CN102032956B - 锅炉水冷壁吸热量的实时测量方法 - Google Patents

Abstract

一种锅炉系统技术领域的锅炉水冷壁吸热量的实时测量方法。包括步骤如下：步骤1：根据锅炉运行规程，获得锅炉参数；步骤2：由IF97公式求出τ时刻下绘制出水冷壁热量测量值Q_dτ随时间变化的趋势曲线；步骤3：由汽包水位模型和水冷壁汽液循环比模型求得汽包和水冷壁中的汽水比例，由下降管流量型求出下降管出口工质的流量和焓值；步骤4：由汽包模型求出汽包惯性常数Ip，并求出水冷壁吸热量实时测量Q。本发明能用于火力电站的水冷壁吸热量在线测量，还能用于燃煤热值的在线辨识以及锅炉效率的实时测量。

Description

锅炉水冷壁吸热量的实时测量方法

技术领域

本发明涉及的是一种锅炉系统技术领域的测量方法，具体地说，是一种锅炉水冷壁吸热量的实时测量方法。

背景技术

目前测量火力电站锅炉水冷壁吸热量的主要方法是基于辐射传热模型进行推算。例如张腾飞等人在《中国电机工程学报》(2003年10月第23卷第10期第215-219页)发表了“炉膛辐射传热数学模型及其仿真”一文，提出一种基于辐射传热数学模型的水冷壁热量的测量方法。它以辐射传热计算的假想面模型为基础，建立了炉内传热集总参数模型，推导出计算真实壁面和假想壁面有效辐射力的矩阵表达式，然后通过求解介质的能量方程获得了炉内的一维温度分布以及各区段水冷壁的吸热量。

然而，辐射传热模型不仅有大量的经验参数，并且需要炉膛传热面积、内部结构等信息，尤其需要知道炉膛温度分布。近年出现了一些测量炉膛温度的方法，如余兵等人在《热力发电》(2007年第36卷第6期第97-100页)发表了“300MW机组锅炉三维温度场可视化监控系统的研究”一文，它利用计算机图像处理技术设计了一套电站锅炉双炉膛三维温度场监测系统，实现了对炉膛燃烧状况实时监测和诊断，是一种锅炉三维温度场的可视化监测方法。但是，目前国内绝大多数火力电站并没有投运这样的设备。也就是说，对于大量火力电站而言，炉膛温度事实上是没有有效测量值的，从而无法应用辐射传热模型。

发明内容

本发明针对现有技术的不足和缺陷，提出了一种锅炉水冷壁吸热量的实时测量方法。与传统方法相比，本发明无需炉膛温度或温度分布的测量值，而是根据汽水系统的实时运行状况，利用汽水系统质量和能量动态衡算的方法，测量出锅炉水冷壁的实时吸热量。

本发明通过以下技术方案实现：

本发明包括步骤如下：

步骤1：根据锅炉运行规程，获得以下锅炉参数：汽包内径D₁、汽包外径D₂、汽包长度L、下降管和水冷壁高度h、水冷壁直径d、水冷壁根数n、下降管根数n_x、下降管内径d_xn、下降管外径d_xw、金属密度ρ_metal、金属比热c_m。在τ时刻，从厂级监控信息系统的实时数据库读取该时刻下的汽包压力P_t、给水流量Fl_w、主蒸汽流量Fl_s、汽包水位hqb、排污量pw、给水温度T_sm、机组负荷指令P_fh。

步骤2：由IF97公式求出τ时刻下、汽包压力为P_t时的水密度ρ_w、汽密度ρ_s、水焓值i_w、汽焓值i_s，汽包压力变化d_pd时的水焓值变化量

汽焓值变化量

水密度变化量

汽密度变化量

工质饱和温度变化量d_Ts。参见：W.Wagner在《Journal of Engineeringfor Gas Turbines and Power》(JANUARY 2000，Vol.122/151)发表的“The IAPWSIndustrial Formulation 1997 for the Thermodynamic Properties of Water and Steam”一文中给出的公式(以下简称IF97公式)。

步骤3：由汽包水位模型和水冷壁汽液循环比模型求得汽包和水冷壁中的汽水比例，由下降管流量模型求出下降管出口工质的流量和焓值。

步骤4：由汽包模型求出汽包惯性常数Ip，并求出水冷壁吸热量实时测量值Q。

所述的τ时刻，是指在任意运行时刻τ，以时问先后为序，分别获得τ＝t、t+Δt、t+2Δt、……t+nΔt时刻下相应的水冷壁热量值Q_dt、Q_d(t+Δt)、Q_d(t+2Δt)、……、Q_d(t+nΔt)，并绘制出水冷壁热量测量值Q_dτ随时问变化的趋势曲线；

所述汽包水位模型为：

$R_1=\frac{D_1}{2}$

$V_1=\frac{4\mathrm{\π}{R}_1^3}{3}+\mathrm{\π}{R}_1^{2}L$

${\mathrm{hqb}}_{\mathrm{cc}}=R_1+0.051+\frac{\mathrm{hqb}}{1000}$

$S_1=R_1^2\·\arccos (1-\frac{{\mathrm{hqb}}_{\mathrm{cc}}}{R_1})-\sqrt{2R_1\·{\mathrm{hqb}}_{\mathrm{cc}}-{\mathrm{hqb}}_{\mathrm{cc}}^2(R_1-{\mathrm{hqb}}_{\mathrm{cc}})}$

${\mathrm{vw}}_{\mathrm{qb}}=\mathrm{\π}(R_1\·{\mathrm{hqb}}_{\mathrm{cc}}^2-\frac{{\mathrm{hqb}}_{\mathrm{cc}}^3}{3})+S_1\·L$

vs_qb＝V₁-vw_qb

式中：R₁为汽包内半径，V₁为汽包体积，hqb_cc为修正后的汽包水位高度，S₁为汽包圆柱部分截面积，vw_qb为汽包中水的体积，vs_qb为汽包中汽的体积；

所述水冷壁汽液循环比模型为：

${\mathrm{hvw}}_{\mathrm{slb}}=0.3-\frac{P_{\mathrm{fh}}-{\mathrm{mid}}_{\mathrm{fh}}}{7\·{\mathrm{mid}}_{\mathrm{fh}}}$

${\mathrm{kvs}}_{\mathrm{slb}}=0.7+\frac{P_{\mathrm{fh}}-{\mathrm{mid}}_{\mathrm{fh}}}{7\·{\mathrm{mid}}_{\mathrm{fh}}}$

其中：mid_fh为火力电站机组负荷中位值，kvw_slb为水冷壁中水体积比例，kvs_slb为水冷壁中汽体积比例；

所述的下降管流量模型为：

$V_{\mathrm{xj}}=\mathrm{\π}{\left(\frac{d_{\mathrm{xn}}}{2}\right)}^{2}h\·n_x$

$V_2=\mathrm{\π}{\left(\frac{d}{2}\right)}^{2}h\·n$

$\mathrm{Dxje}=750\mathrm{\π}{\left(\frac{d_{\mathrm{xn}}}{2}\right)}^2{n}_x$

Vsh_cc＝V₂·kvs_slb

${\mathrm{Dxj}}_{\mathrm{cc}}=\frac{\mathrm{Dxje}}{\frac{0.15V_2}{4{\mathrm{Vsh}}_{\mathrm{cc}}}+0.85}$

$\mathrm{Hxj}=\frac{0.5{\mathrm{Fl}}_w\·i_s+({\mathrm{Dxj}}_{\mathrm{cc}}-0.5{\mathrm{Fl}}_w)i_w}{{\mathrm{Dxj}}_{\mathrm{cc}}}$

${\mathrm{Mm}}_{\mathrm{xjg}}=(\mathrm{\π}{\left(\frac{d_{\mathrm{xw}}}{2}\right)}^2-\mathrm{\π}{\left(\frac{d_{\mathrm{xn}}}{2}\right)}^2)h\·n_x\·{\mathrm{\ρ}}_{\mathrm{metal}}$

$d_{\mathrm{Hxjh}}={\mathrm{Dxj}}_{\mathrm{cc}}\·(\mathrm{Hxj}-{\mathrm{Hxjh}}_{\mathrm{cc}})-V_{\mathrm{xj}}\·{\mathrm{Hxjh}}_{\mathrm{cc}}\·d_{{\mathrm{\ρ}}_w}-\frac{c_m\·{\mathrm{Mm}}_{\mathrm{xjg}}\·d_{\mathrm{Ts}}}{V_{\mathrm{xj}}\·{\mathrm{\ρ}}_w}$

Hxj_glb＝Hxj_cc+d_Hxjh

式中：V_xj为下降管体积，V₂为水冷壁内体积，Dxje为额定负荷下下降管入口工质流量，Vsh_cc为水冷壁中水体积，Dxj_cc为下降管入口工质流量，Hxj为下降管入口工质焓值，Mm_xjg为下降管金属质量，Hxjh_cc为上一步下降管出口焓值，d_Hxjh为该步下降管焓值变化量，Hxjh_glb为当前下降管出口焓值。

所述的汽包模型为：

${\mathrm{Mm}}_{\tan k}=(\frac{4\mathrm{\π}{R}_2^3}{3}-\frac{4\mathrm{\π}{R}_1^3}{3}+(\mathrm{\π}{R}_2^2-\mathrm{\π}{R}_1^2)L)\·{\mathrm{\ρ}}_{\mathrm{metal}}$

del_i＝1083+3.99P_t-4.67T_sm

$\mathrm{Ip}=\mathrm{Vxj}\·{\mathrm{\ρ}}_w\·d_{i_w}+V_{\mathrm{xj}}\·{\mathrm{Hxjh}}_{\mathrm{cc}}\·d_{{\mathrm{\ρ}}_w}+c_m({\mathrm{Mm}}_{\tan k}+{\mathrm{Mm}}_{\mathrm{xjh}})d_{\mathrm{Ts}}$

$+({\mathrm{\ρ}}_w\·d_{i_w}+\frac{r_q\·{\mathrm{\ρ}}_s\·d_{{\mathrm{\ρ}}_w}}{{\mathrm{\ρ}}_w-{\mathrm{\ρ}}_s}(V_2\·{\mathrm{kvw}}_{\mathrm{slb}}+{\mathrm{vw}}_{\mathrm{qb}})+({\mathrm{\ρ}}_s\·d_{i_s}+\frac{r_q\·{\mathrm{\ρ}}_w\·d_{{\mathrm{\ρ}}_s}}{{\mathrm{\ρ}}_w-{\mathrm{\ρ}}_s}(V_2\·{\mathrm{kvs}}_{\mathrm{slb}}+{\mathrm{vs}}_{\mathrm{qb}}))$

$Q=d_{\mathrm{pd}}\·\mathrm{Ip}-(\frac{r_q\·{\mathrm{\ρ}}_s}{{\mathrm{\ρ}}_s-{\mathrm{\ρ}}_s}-{\mathrm{del}}_i)\·{\mathrm{Fl}}_w+\frac{r_q\·{\mathrm{\ρ}}_w}{{\mathrm{\ρ}}_w-{\mathrm{\ρ}}_s}{\mathrm{Fl}}_s+\frac{r_q\·{\mathrm{\ρ}}_s}{{\mathrm{\ρ}}_w-{\mathrm{\ρ}}_s}\mathrm{pw}$

式中：Mm_tank为汽包金属质量，Mm_xjh为下降管的金属质量，R₁、R₂分别为汽包的内径和外径，del_i为给水欠焓，r_q为汽化潜热，Ip为汽包惯性常数，Q为水冷壁吸收热量。

本发明与以往的基于辐射传热模型的测量方法相比，该方法给出了在满足锅炉汽水系统质量和能量守恒关系下的水冷壁吸热量的实时测量，并不需要炉膛温度、传热面积和其它结构等信息。在每一时刻的水冷壁吸热量测量中，都将锅炉汽水系统的测点信息统一起来，得到相应时刻下的水冷壁吸热量。在下一时刻，基于同样的质量和能量守恒关系可测量出水冷壁吸热量，它与前一时刻水冷壁吸热量的差值可真实、准确的反映出水冷壁吸热量随时问的变化量。作为一种实时测量方法，水冷壁吸热量所需的测点数据均从厂级监控信息系统的实时数据库直接读取，现场不需要额外增加测点、探头或昂贵的仪表，只需在已有的控制系统中增加相应的软件模块即可，成本低。水冷壁吸热量的软测量是迭代进行的，迭代中仅需要保存上一步的结果，内存占用少，实时性强。可将本方法加载到火力电站控制系统中，用于指导锅炉的燃烧调节和优化控制，进一步完善整个机组性能，扩展应用范围广。

本发明能用于火力电站的水冷壁吸热量在线测量。若辅以合适的模型进一步测量后续换热器设备的吸热量、烟气损失热量、燃煤机械热损失和化学不完全燃烧热损失等，则它还能用于燃煤热值的在线辨识以及锅炉效率的实时测量。

附图说明

图1锅炉汽水系统结构示意图

图2水冷壁热量与炉膛燃烧总热量随时问变化的曲线示意图

其中：水冷壁吸热量(实时测量值)与炉膛燃烧总热量(实测值)随时问变化的曲线。

具体实施方式

下面结合附图对本发明的实施例作详细说明：本实施例在本发明技术方案前提下进行实施，给出了详细的实施方式和具体操作过程，但本发明的保护范围不限于下述的实施例。

如图1所示，以某300MW火力电站为例，该自然循环锅炉汽水系统有下降管、水冷壁和汽包三大部分，锅炉在任意运行时刻τ，以时问先后为序，分别获得τ＝t、t+Δt、t+2Δt、……t+nΔt时刻下相应的水冷壁热量值Q_dt、Q_d(t+Δt)、Q_d(t+2Δt)、……、Q_d(t+nΔt)，并绘制出水冷壁热量测量值Q_dτ随时问变化的趋势曲线。实现基于汽水系统质量和能量守恒的锅炉水冷壁吸热量的软测量方法。锅炉型号为HG-1025/17.3-WM18型，锅炉为亚临界、自然循环、一次中问再热、“W”火焰燃烧方式、双拱单炉膛、平衡通风、尾部双烟道、烟气挡板调温、固态排渣、露天布置、全钢架悬吊式汽包炉；机组配置容量均为15％BMCR的高、低压两级串级旁路。

本实施例步骤如下：

步骤1：根据锅炉运行规程，获得以下固定参数：

汽包内径1.83m、汽包外径2.19m、汽包长度15.8m、下降管和水冷壁高度48.0m、水冷壁直径0.0527m、水冷壁根数624、下降管根数4、下降管内径0.325m、下降管外径0.406m、金属密度7.8×10³kg/m³、金属比热0.502kJ/kg℃。

从厂级监控信息系统的实时数据库中读取相关实时数据，如在τ＝t时刻的运行工况下，主要的实时数据如下：

汽包压力17.6MPa、给水流量236.6kg/s、主蒸汽流量233.6kg/s、汽包水位-133.5mm、排污量5.0kg/s、给水温度292.3℃、机组负荷指令278.5MW。

步骤2：由IF97公式，求出τ＝t时刻、汽包压力为P_t时的水密度ρ_w、汽密度ρ_s、水焓值i_w、汽焓值i_s，汽包压力变化d_pd时的水焓值变化量

汽焓值变化量

水密度变化量

汽密度变化量

工质饱和温度变化量d_Ts。

步骤3：由汽包水位模型和水冷壁汽液循环比模型求得汽包和水冷壁中的汽水比例，由下降管流量模型求出下降管出口工质的流量和焓值。

步骤4：由汽包模型求出汽包惯性常数Ip，并求出水冷壁吸热量实时测量值Q。

步骤5：选取时问问隔Δt＝5s，记录τ＝t时刻起一天的水冷壁吸热量变化趋势，在每一时刻τ＝t、t+Δt、t+2Δt、……、t+17280Δt下，重复步骤1-4，获得相应时刻的水冷壁热量软测量值Q_dt、Q_d(t+Δt)、Q_d(t+2Δt)、……、Q_{d(t+17280Δt)}，并绘制出水冷壁热量测量值Qdτ随时问变化的趋势曲线。示例运行结果如图2所示。作为参考，其中的炉膛燃烧总热量是用给煤量乘以燃煤热值标准值得到的。

本实施例能用于火力电站的水冷壁吸热量在线测量，还能用于燃煤热值的在线辨识以及锅炉效率的实时测量。

Claims (6)

1.一种锅炉水冷壁吸热量的实时测量方法，其特征在于，包括步骤如下：

步骤1：根据锅炉运行规程，获得汽包内径D₁、汽包外径D₂、汽包长度L、下降管和水冷壁高度h、水冷壁直径d、水冷壁限数n、下降管根数n_x、下降管内径d_xn、下降管外径d_xw、金属密度ρ_metal和金属比热c_m，在τ时刻，从厂级监控信息系统的实时数据库读取该时刻下的汽包压力P_t、给水流量Fl_w、主蒸汽流量Fl_s、汽包水位hqb、排污量pw、给水温度T_sm、机组负荷指令P_fh；

步骤2：由IF97公式求出τ时刻下，汽包压力为P_t时的水密度ρ_w、汽密度ρ_s、水焓值i_w、汽焓值i_s，汽包压力变化d_pd时的水焓值变化量

汽焓值变化量

水密度变化量

汽密度变化量

工质饱和温度变化量d_Ts；

步骤3：由汽包水位模型和水冷壁汽液循环比模型求得汽包和水冷壁中的汽水比例，由下降管流量模型求出下降管出口工质的流量和焓值；

步骤4：由汽包模型求出汽包惯性常数Ip，并求出水冷壁吸热量实时测量Q。

2.根据权利要求1所述的锅炉水冷壁吸热量的实时测量方法，其特征是，所述的τ时刻，是指任意运行时刻τ，以时间先后为序，分别获得τ＝t、t+Δt、t+2Δt、……t+nΔt时刻下相应的水冷壁热量值Q_dt、Q_d(t+Δt)、Q_d(t+2Δt)、……、Q_d(t+nΔt)，并绘制出水冷壁热量测量值Qdτ随时间变化的趋势曲线。

3.根据权利要求1所述的锅炉水冷壁吸热量实时测量方法，其特征是，所述的汽包水位模型为：

$R_1=\frac{D_1}{2}$

$V_1=\frac{4\mathrm{\π}{R}_1^3}{3}+\mathrm{\π}{R}_1^{2}L$

${\mathrm{hqb}}_{\mathrm{cc}}=R_1+0.051+\frac{\mathrm{hqb}}{1000}$

$S_1=R_1^2\·\arccos (1-\frac{{\mathrm{hqb}}_{\mathrm{cc}}}{R_1})-\sqrt{2R_1\·{\mathrm{hqb}}_{\mathrm{cc}}-{\mathrm{hqb}}_{\mathrm{cc}}^2(R_1-{\mathrm{hqb}}_{\mathrm{cc}})}$

${\mathrm{vw}}_{\mathrm{qb}}=\mathrm{\π}(R_1\·{\mathrm{hqb}}_{\mathrm{cc}}^2-\frac{{\mathrm{hqb}}_{\mathrm{cc}}^3}{3})+S_1\·L$

vs_qb＝V₁-vw_qb

式中：R₁为汽包内半径，V₁为汽包体积，hqb_cc为修正后的汽包水位高度，S₁为汽包圆柱部分截面积，vw_qb为汽包中水的体积，vs_qb为汽包中汽的体积。

4.根据权利要求1所述的锅炉水冷壁吸热量实时测量方法，其特征是，所述的水冷壁汽液循环比模型为：

${\mathrm{hvw}}_{\mathrm{slb}}=0.3-\frac{P_{\mathrm{fh}}-{\mathrm{mid}}_{\mathrm{fh}}}{7\·{\mathrm{mid}}_{\mathrm{fh}}}$

${\mathrm{kvs}}_{\mathrm{slb}}=0.7+\frac{P_{\mathrm{fh}}-{\mathrm{mid}}_{\mathrm{fh}}}{7\·{\mathrm{mid}}_{\mathrm{fh}}}$

式中：mid_fh为火力电站机组负荷中位值，kvw_slb为水冷壁中水体积比例，kvs_slb为水冷壁中汽体积比例。

5.根据权利要求1所述的锅炉水冷壁吸热量实时测量方法，其特征是，所述的下降管流量模型为：

$V_{\mathrm{xj}}=\mathrm{\π}{\left(\frac{d_{\mathrm{xn}}}{2}\right)}^{2}h\·n_x$

$V_2=\mathrm{\π}{\left(\frac{d}{2}\right)}^{2}h\·n$

$\mathrm{Dxje}=750\mathrm{\π}{\left(\frac{d_{\mathrm{xn}}}{2}\right)}^2{n}_x$

Vsh_cc＝V₂·kvs_slb

${\mathrm{Dxj}}_{\mathrm{cc}}=\frac{\mathrm{Dxje}}{\frac{0.15V_2}{4{\mathrm{Vsh}}_{\mathrm{cc}}}+0.85}$

$\mathrm{Hxj}=\frac{0.5{\mathrm{Fl}}_w\·i_s+({\mathrm{Dxj}}_{\mathrm{cc}}-0.5{\mathrm{Fl}}_w)i_w}{{\mathrm{Dxj}}_{\mathrm{cc}}}$

${\mathrm{Mm}}_{\mathrm{xjg}}=(\mathrm{\π}{\left(\frac{d_{\mathrm{xw}}}{2}\right)}^2-\mathrm{\π}{\left(\frac{d_{\mathrm{xn}}}{2}\right)}^2)h\·n_x\·{\mathrm{\ρ}}_{\mathrm{metal}}$

$d_{\mathrm{Hxjh}}={\mathrm{Dxj}}_{\mathrm{cc}}\·(\mathrm{Hxj}-{\mathrm{Hxjh}}_{\mathrm{cc}})-V_{\mathrm{xj}}\·{\mathrm{Hxjh}}_{\mathrm{cc}}\·d_{{\mathrm{\ρ}}_w}-\frac{c_m\·{\mathrm{Mm}}_{\mathrm{xjg}}\·d_{\mathrm{Ts}}}{V_{\mathrm{xj}}\·{\mathrm{\ρ}}_w}$

Hxj_glb＝Hxjh_cc+d_Hxjh

式中：V_xj为下降管体积，V₂为水冷壁内体积，Dxje为额定负荷下下降管入口工质流量，Vsh_cc为水冷壁中水体积，Dxj_cc为下降管入口工质流量，Hxj为下降管入口工质焓值，Mm_xjg为下降管金属质量，Hxjh_cc为上一步下降管出口焓值，d_Hxjh为该步下降管焓值变化量，Hxjh_glb为当前下降管出口焓值。

6.根据权利要求1所述的锅炉水冷壁吸热量实时测量方法，其特征是，所述的汽包模型为：

${\mathrm{Mm}}_{\tan k}=(\frac{4\mathrm{\π}{R}_2^3}{3}-\frac{4\mathrm{\π}{R}_1^3}{3}+(\mathrm{\π}{R}_2^2-\mathrm{\π}{R}_1^2)L)\·{\mathrm{\ρ}}_{\mathrm{metal}}$

del_i＝1083+3.99P_t-4.67T_sm

$\mathrm{Ip}=\mathrm{Vxj}\·{\mathrm{\ρ}}_w\·d_{i_w}+V_{\mathrm{xj}}\·{\mathrm{Hxjh}}_{\mathrm{cc}}\·d_{{\mathrm{\ρ}}_w}+c_m({\mathrm{Mm}}_{\tan k}+{\mathrm{Mm}}_{\mathrm{xjh}})d_{\mathrm{Ts}}$

$+({\mathrm{\ρ}}_w\·d_{i_w}+\frac{r_q\·{\mathrm{\ρ}}_s\·d_{{\mathrm{\ρ}}_w}}{{\mathrm{\ρ}}_w-{\mathrm{\ρ}}_s}(V_2\·{\mathrm{kvw}}_{\mathrm{slb}}+{\mathrm{vw}}_{\mathrm{qb}})+({\mathrm{\ρ}}_s\·d_{i_s}+\frac{r_q\·{\mathrm{\ρ}}_w\·d_{{\mathrm{\ρ}}_s}}{{\mathrm{\ρ}}_w-{\mathrm{\ρ}}_s}(V_2\·{\mathrm{kvs}}_{\mathrm{slb}}+{\mathrm{vs}}_{\mathrm{qb}}))$

$Q=d_{\mathrm{pd}}\·\mathrm{Ip}-(\frac{r_q\·{\mathrm{\ρ}}_s}{{\mathrm{\ρ}}_s-{\mathrm{\ρ}}_s}-{\mathrm{del}}_i)\·{\mathrm{Fl}}_w+\frac{r_q\·{\mathrm{\ρ}}_w}{{\mathrm{\ρ}}_w-{\mathrm{\ρ}}_s}{\mathrm{Fl}}_s+\frac{r_q\·{\mathrm{\ρ}}_s}{{\mathrm{\ρ}}_w-{\mathrm{\ρ}}_s}\mathrm{pw}$

式中：Mm_tank为汽包金属质量，Mm_xjh为下降管的金属质量，R₂分别为汽包的内径和外径，del_i为给水欠焓，r_q为汽化潜热，Ip为汽包惯性常数，Q为水冷壁吸收热量。

Images