一种电站锅炉炉膛结渣分区段实时软测量方法与流程

技术特征：

1.一种电站锅炉炉膛结渣分区段实时软测量方法，其特征在于，包括以下步骤：

步骤1，将炉膛按燃烧性质由下到上依次分为主燃烧区、燃尽区和换热区三个区段。对炉膛进行相关数据采集。

步骤2，根据步骤1中采集的相关数据的分别计算燃料理论燃烧温度T_th、炉膛整体辐射吸收减弱系数k_a、炉膛整体火焰综合黑度ε_syn、炉膛整体黑度和炉膛整体水冷壁热有效系数ψ_f，为炉膛分区段辐射传热计算作准备。

步骤3，假设主燃烧区出口烟温T_f1″，列出主燃烧区热平衡方程，计算主燃烧区水冷壁热有效系数ψ₁。

步骤4，列出炉膛内主燃烧区热负荷方程，对假设的主燃烧区出口烟温T_f1″进行校核，若满足条件则进行步骤5，若不满足则重复步骤3到4。

步骤5，假设燃尽区的出口烟气温度T_f2″，列出燃尽区热平衡方程，计算燃尽区水冷壁热有效系数ψ₂。

步骤6，列出炉膛内燃尽区热负荷方程，对假设的燃尽区出口烟温T_f2″进行校核，若满足条件则进行步骤7，若不满足则重复步骤5到6。

步骤7，列出换热区热平衡方程，计算换热区水冷壁热有效系数ψ₃，输出计算出的不同区段水冷热有效系数ψ₁、ψ₂和ψ₃作为不同区段的结渣监测指标。

2.根据权利要求1所述的电站锅炉炉膛结渣分区段实时软测量方法，其特征在于：所述步骤1中采集的相关数据包括锅炉炉膛结构及设计参数、入炉煤质数据和锅炉实时运行参数；锅炉炉膛结构及设计参数包括炉膛整体传热面积、不同区段的传热面积、有效容积、计算高度、上下排燃烧器布置高度差、燃烧器平均布置高度、出口烟窗面积、炉膛出口对半辐射受热面的辐射热有效系数；入炉煤质数据通过煤质分析获得，主要包括煤的元素分析、工业分析和热值分析；锅炉实时运行参数通过电厂DCS系统采集，主要测点包括锅炉燃煤量、一次风占总风量比例、二次风占总风量比例、一次风进出口风温、二次风进出口风温、水冷壁工质流量、水冷壁进出口工质压力、水冷壁进出口工质温度和炉膛出口烟气温度。

3.根据权利要求1所述的电站锅炉炉膛结渣分区段实时软测量方法，其特征在于：所述步骤2中根据采集数据的分别计算燃料理论燃烧温度T_th、炉膛整体辐射吸收减弱系数k_a、炉膛整体火焰综合黑度ε_syn、炉膛整体黑度和炉膛整体水冷壁热有效系数ψ_f的方法如下：

步骤2-1，根据理论的热空气焓和理论的冷空气焓得到进入炉内空气热量；

$Q_k=({\mathrm{\α}}_f^{\′\′}-{\mathrm{\Δ\α}}_f-{\mathrm{\Δ\α}}_{pcs})I_{hk}^0+({\mathrm{\Δ\α}}_f+{\mathrm{\Δ\α}}_{pcs})I_{ck}^0$

Q_k为随单位质量燃料带入炉内的空气的热量，且此热量含漏风的热量；为理论的热空气焓，即空预器出口一次风和二次风混合温度；为理论的冷空气焓，即空预器进口一次风和二次风混合温度，α_f"为炉膛出口过量空气系数，Δα_f为炉膛漏风系数，Δα_pcs为制粉系统的漏风系数。；

根据入炉内空气热量、单位质量燃料带入炉内的热量、化学未完全燃烧热损失、机械未完全燃烧热损失以及其他热损失得到燃料带入炉内有效热量；

$Q_f^{ef}=Q_r\frac{100-q_3-q_4-q_6}{100-q_4}+Q_k$

为炉内有效热量，即为理论燃烧温度T_th对应焓值，在获得炉内有效热量后通过烟气焓温表利用插值法查取燃料理论燃烧温度T_th；Q_r为单位质量燃料带入炉内的热量，通常等于燃料收到基低位发热量，q₃为化学未完全燃烧热损失，q₄为机械未完全燃烧热损失，q₆为其他热损失，

步骤2-2：根据收到基灰分以及单位燃料燃烧产生烟气质量得到烟气中灰分颗粒的质量浓度；

${\mathrm{\μ}}_{ash,m}=\frac{A_{ar}{\mathrm{\α}}_{fa}}{100G_y};$

μ_ash,m为烟气中灰分颗粒的质量浓度，A_ar为收到基灰分，α_fa为飞灰系数；G_y为单位燃料燃烧产生烟气质量；

根据收到基碳元素、干燥无灰基挥发份、最上排燃烧器布置高度、最下排燃烧器布置高度以及水冷壁的吸热表面积得到烟气中焦炭颗粒的容积浓度；

${\mathrm{\μ}}_{cok,v}=\frac{5.5C_{ar}(10+q_4)}{(100+V_{daf})V_y}(1+\frac{h_t-h_{un}}{h_f});$

μ_cok,v为烟气中焦炭颗粒的容积浓度，C_ar为收到基碳元素，V_daf为干燥无灰基挥发份，h_t为最上排燃烧器布置高度，h_un为最下排燃烧器布置高度，H_f为水冷壁的吸热表面积；

根据水蒸汽占烟气容积份额、水蒸汽和二氧化物之和占烟气容积份额、炉内辐射层有效厚度以及炉膛整体出口烟温得到三原子气体减弱系数；

$k_{g}r=\[\frac{0.78+1.6r_{H_{2}O}}{\sqrt{rs}}-0.1\](1-0.37\frac{T_f^{\′\′}}{1000})r;$

k_gr为三原子气体减弱系数；为水蒸汽占烟气容积份额，r为水蒸汽和二氧化物之和占烟气容积份额，S为炉内辐射层有效厚度，T_f″为炉膛整体出口烟温；

根据烟气中灰分颗粒的质量浓度以及灰分颗粒的平均粒径得到灰分减弱系数；

$k_{ash}{\mathrm{\μ}}_{ash}=\frac{5330{\mathrm{\μ}}_{ash,m}}{{\left(T_f^{\′\′}{d}_{ash}\right)}^{2/3}}\[1-\frac{0.65}{1+0.0177/{\left({\mathrm{\μ}}_{ash,m}S\right)}^2}\];$

k_ashμ_ash为灰分颗粒减弱系数，μ_ash,m为烟气中灰分颗粒的质量浓度，d_ash为灰分颗粒的平均粒径，

根据焦炭颗粒的平均粒径、烟气中焦炭颗粒的容积浓度得到焦炭颗粒减弱系数；

$k_{cok}{\mathrm{\μ}}_{cok}=\frac{10{\mathrm{\μ}}_{cok,v}}{{\left(T_f^{\′\′}{d}_{cok}\right)}^{2/3}};$

k_cokμ_cok为焦炭颗粒减弱系数，d_cok为焦炭颗粒的平均粒径，

根据三原子气体减弱系数、灰分颗粒减弱系数以及焦炭颗粒减弱系数得到辐射吸收减弱系数；

k_a＝k_gr+k_ashμ_ash+k_cokμ_cok

k_a为炉膛整体辐射吸收减弱系数，

步骤2-3，根据炉膛整体辐射吸收减弱系数、炉内辐射层有效厚度得到炉膛实际火焰黑度；

${\mathrm{\ϵ}}_f=1-e^{-k_{a}S};$

ε_f为炉膛实际火焰黑度；

根据炉膛整体辐射吸收减弱系数、炉内辐射层有效厚度以及炉膛实际火焰黑度得到炉膛火焰综合黑度：

${\mathrm{\ϵ}}_{syn}=\frac{{\mathrm{\ϵ}}_f}{0.32k_a{\mathrm{R\ϵ}}_f+1};$

ε_syn为考虑了火焰辐射强度因介质吸收而减弱的火焰综合黑度；

步骤2-4：根据炉膛火焰综合黑度建立得到炉膛整体黑度关系式：

${\mathrm{\ϵ}}_f^{syn}=\frac{{\mathrm{\ϵ}}_{syn}}{{\mathrm{\ϵ}}_{syn}+(1-{\mathrm{\ϵ}}_{syn}){\mathrm{\ψ}}_f};$

为炉膛整体黑度，ψ_f为炉膛整体水冷壁热有效系数；

根据理论燃烧温度、主燃烧区出口烟气温度、炉膛火焰最高温度位置的相对高度以及炉内烟气在理论燃烧温度至炉膛出口温度区间内的平均热容建立炉膛辐射传热关系式：

T_th为理论燃烧温度，x_m为炉膛火焰最高温度位置的相对高度，σ₀为玻尔兹曼常数，φ为保热系数，B_j为计算燃烧量，为炉内烟气在理论燃烧温度至炉膛出口温度区间内的平均热容；

联立炉膛整体黑度关系式和炉膛辐射传热关系式两个关系式，求出炉膛整体水冷壁热有效系数：

${\mathrm{\ψ}}_f=\frac{{\mathrm{Z\ϵ}}_{syn}}{{\mathrm{\ϵ}}_{syn}-(1-{\mathrm{\ϵ}}_{syn})Z};$

其中，

求出炉膛炉膛整体水冷壁热有效系数后，再计算出炉膛的整体黑度

4.根据权利要求1所述的电站锅炉炉膛结渣分区段实时软测量方法，其特征在于：所述步骤3中主燃烧区热平衡方程如下：

$\begin{array}{c}{B}_j(\frac{100}{100-q_4}{\mathrm{\β}}_{cr}{Q}_r+Q_k-Q_6)\\ =B_j{I}_{f1}^{\′\′}+{\mathrm{\σ}}_0{\mathrm{\ϵ}}_f^{syn}{\mathrm{\ψ}}^{\′\′}{\left(T_{f1}^{\′\′}\right)}^4{F}_{c1}+{\mathrm{\σ}}_0{\mathrm{\ϵ}}_f^{syn}{\mathrm{\ψ}}_1{\left(T_{f1}^{\′\′}\right)}^4{H}_{f1}\end{array}$

根据主燃烧区热平衡式，计算主燃烧区水冷壁热有效系数ψ₁。

其中，T_f1″为主燃烧区出口烟气温度。I_f1″为主燃烧区出口烟气焓值(根据T_f1″查取烟气焓温表获得)。β_cr为燃料的燃尽率，可查阅锅炉手册。ψ"为对上面区段的辐射热有效系数，一般取0.1。F_c1主燃烧区出口炉膛截面积。H_f1为主燃烧区的水冷壁传热面积。ψ₁为主燃烧区水冷壁热有效系数。

5.所述步骤4中炉膛内主燃烧区热负荷方程如下：

$\stackrel{\‾}{q}=B_j(Q_f^{ef}-I_f^{\′\′})/H_f$

$\stackrel{\‾}{q_I}=\frac{{\∫}_0^{h_1}\left({\mathrm{\η}}_h\stackrel{\‾}{q}\right)dh}{h_1}$

主燃烧区出口烟温T_f1″校核式：

$|\stackrel{\‾}{q_I}-{\mathrm{\σ}}_0{\mathrm{\ϵ}}_f^{syn}{\mathrm{\ψ}}_1{\left(T_{f1}^{\′\′}\right)}^4|/\stackrel{\‾}{q_I}\≤2\%$

其中，为炉膛整体平均热负荷。η_h为炉膛内局部热负荷系数，可做试验确定或查锅炉说明书。为主燃烧区的平均热负荷。h₁为主燃烧区出口到冷灰斗中线高度。

6.所述步骤5中燃尽区热平衡方程如下：

$\begin{array}{c}{B}_j{\mathrm{\Δ\β}}_{cr}{Q}_r+B_j{I}_{f1}^{\′\′}+{\mathrm{\σ}}_0{\mathrm{\ϵ}}_f^{syn}{\mathrm{\ψ}}^{\′\′}{\left(T_{f1}^{\′\′}\right)}^4{F}_{c1}\\ =B_j{I}_{f2}^{\′\′}+{\mathrm{\σ}}_0{\mathrm{\ϵ}}_f^{syn}{\mathrm{\ψ}}^{\′\′}{\left(T_{f2}^{\′\′}\right)}^4{F}_{c2}+{\mathrm{\σ}}_0{\mathrm{\ϵ}}_f^{syn}{\mathrm{\ψ}}_2{T}_{g2}^4{H}_{f2}\end{array}$

根据燃尽区平均烟气温度计算燃尽区水冷壁热有效系数ψ₂。

其中，T_g2为燃尽区的烟气平均温度。T_f2″为燃尽区出口烟气温度。I_f2″为燃尽区出口烟气焓值，kJ/kg(根据T_f2″查取烟气焓温表获得)。Δβ_cr为主燃烧区燃料的未燃尽率。F_c2燃尽区出口炉膛截面积。H_f2为燃尽区的水冷壁传热面积。ψ₂为燃尽区水冷壁热有效系数。

7.所述步骤6中炉膛内燃尽区热负荷方程如下：

$\stackrel{\‾}{q_{II}}=\frac{{\∫}_{h_1}^{h_2}\left({\mathrm{\η}}_h\stackrel{\‾}{q}\right)dh}{h_2-h_1}$

燃尽区出口烟温T_f2″校核式：

$|\stackrel{\‾}{q_{II}}-{\mathrm{\σ}}_0{\mathrm{\ϵ}}_f^{syn}{\mathrm{\ψ}}_2{T}_{g2}^4|/\stackrel{\‾}{q_{II}}\≤2\%$

其中，为燃尽区的平均热负荷。h₂为燃尽区出口到冷灰斗中线高度。

8.所述步骤7中换热区热平衡方程如下：

$\begin{array}{c}{B}_j{I}_{f2}^{\′\′}+{\mathrm{\σ}}_0{\mathrm{\ϵ}}_f^{syn}{\mathrm{\ψ}}^{\′\′}{\left(T_{f2}^{\′\′}\right)}^4{F}_{c2}\\ =B_j{I}_f^{\′\′}+{\mathrm{\σ}}_0{\mathrm{\ϵ}}_f^{syn}{\mathrm{\ψ}}_p{\left(T_f^{\′\′}\right)}^4{F}_{c3}+{\mathrm{\σ}}_0{\mathrm{\ϵ}}_f^{syn}{\mathrm{\ψ}}_3{T}_{g3}^4{H}_{f3}\end{array}$

根据换热区平均烟气温度计算换热区水冷壁热有效系数ψ₃。

其中，T_g3为换热区的烟气平均温度。I_f″为燃尽区出口烟气焓值(根据T_f″查取烟气焓温表获得)。ψ_p为炉膛出口对半辐射受热面的辐射热有效系数，等于炉膛整体水冷壁热有效系数ψ_f乘以经验系数，经验系数可按锅炉设计值选取。F_c3炉膛出口截面积。H_f3为换热区的水冷壁传热面积。ψ₃为换热区水冷壁热有效系数。