回转式空预器直接漏风流量的实时估计方法与流程

本发明涉及火力发电技术领域，具体地，涉及回转式空预器直接漏风流量的实时估计方法。

背景技术：

火力发电站为了利用烟气余热，在烟道尾部安装回转式空预器，以预加热进入炉膛的助燃空气。回转式空预器主要由外壳和转子组成，转子包括数以万计的蓄热元件，烟气和空气之间用密封隔板分隔为独立仓室。当转子缓慢旋转时，烟气和空气沿各自仓室逆向流经蓄热元件。蓄热元件在烟气侧吸热，空气侧放热，从而实现烟气向空气的传热。

目前，绝大多数锅炉均采用三分仓回转式空预器(烟气仓、一次风仓、二次风仓)。回转式空预器的漏风分为携带漏风与直接漏风，携带漏风是指在回转式空预器的转子转动过程中，部分驻留在转子内部空间中的空气被携带到下一个分仓；直接漏风是指通过密封间隙漏入烟气侧的空气。其中，回转式空预器的直接漏风会造成以下后果：

1)送风机、引风机电耗增大；

2)进风量减少，燃烧不稳定，锅炉出力下降；

3)排烟热损失增大，锅炉热效率降低。

因此，空预器的直接漏风质量流量是机组运行的一项重要监控指标。

经过对现有技术的检索，中国专利申请号201410080958.6，公开号103886188a，公开了一种基于压差的空预器漏风率的实时估计方法。该方法根据烟气侧与空气侧能量衡算计算瞬时总漏风质量流量，进而得到瞬时漏风面积。然后在日平均漏风面积的基础上，利用加权平均更新平均漏风面积。基于空气侧和烟气侧的压差和平均漏风面积，计算瞬时总漏风质量流量。

但是，上述方法在计算烟气与空气的压差时，将一次风和二次风作为一个整体考虑，而直接漏风本质上是一次风直接漏风和二次风直接漏风之和，因此，其对直接漏风质量流量的估计精度不高。

技术实现要素：

针对现有技术中的缺陷，本发明的目的是提供一种回转式空预器直接漏风流量的实时估计方法。

根据本发明提供的回转式空预器直接漏风流量的实时估计方法，可以包括：

沿气体流动方向和圆周方向，将回转式空预器离散化为多个传热微元；

对各个传热微元建立能量守恒模型，通过所述能量守恒模型得到考虑携带漏风影响的回转式空预器温度分布；其中，所述回转式空预器温度分布包括：回转式空预器的金属壁温度分布、流经回转式空预器的气体温度分布；

根据所述回转式空预器温度分布，实时估计所述回转式空预器的直接漏风流量。

可选地，在沿气体流动方向和圆周方向，将回转式空预器离散化为多个传热微元之前，还包括：

从锅炉集控运行规程中获取回转式空预器的结构参数，所述结构参数包括：转子直径，转子在高温段、中温段、低温段的高度，在高温段、中温段、低温段内传热微元所用材料，一次风仓、二次风仓、烟气仓所占的角度；

从分布式控制系统dcs的实时数据库中，获取给定时刻下所述回转式空预器的当前检测参数，所述当前检测参数包括：烟气入口温度、压力，烟气出口温度、压力，引风机烟气体积流量，一次风入口温度、压力，一次风出口温度、压力，进入磨煤机的一次风体积流量，二次风入口温度、压力，二次风出口温度、压力，送风机的二次风体积流量。

可选地，对各个传热微元建立能量守恒模型，通过所述能量守恒模型得到考虑携带漏风影响的回转式空预器温度分布，包括：

计算气体流经所述回转式空预器所需时间，计算公式如下：

其中：hap为回转式空预器高度，单位m；ugs为气体流速，单位m/s；tgs为气体流经所述回转式空预器所需时间，单位s；

计算所述回转式空预器转子转过的角度，计算公式如下：

其中：ω为回转式空预器转子的转动角速度，单位rad/s；βet为所述回转式空预器转子转过的角度，单位rad；

计算携带漏风体积流量，计算公式如下：

其中：fpa为一次风出口的体积流量，单位m³/s；npa为一次风仓沿圆周方向的离散化份数；βu,pa为一次风仓内一列传热微元对应角度，单位rad；fpa,et为携带漏风体积流量，单位m³/s；

流经烟气仓第一列传热微元的烟气体积流量的计算公式如下：

其中：ffg为烟气出口的体积流量，单位m³/s；nfg为烟气仓沿圆周方向的离散化份数；βu,fg为烟气仓内一列传热微元对应角度，单位rad；ffg,1为流经烟气仓第一列传热微元的烟气体积流量，单位m³/s；

对于烟气仓内任意传热微元(i,j)，其中i、j分别是该传热微元所在行数和列数，若j≠1，则根据能量守恒关系，烟气在传热微元的进、出口的能量之差等于烟气与金属的对流换热量，即满足如下公式：

其中：ρfg为烟气的密度，单位kg/m³；cfg为烟气的比热容，单位kj/(kg·℃)；t'fg(i,j)和t”fg(i,j)分别是传热微元(i,j)的进、出口烟气温度，单位℃；αfg为烟气对流传热系数，kw/(m²·℃)；aconv(i,j)为传热微元(i,j)的对流换热面积，单位m²；tm(i,j)和tm(i,j+1)分别为传热微元(i,j)左侧金属壁温度、右侧金属壁温度，单位℃；

若j＝1，则所述传热微元位于烟气仓内的第一列，根据能量守恒关系，满足如下公式：

其中：t’fg(i,1)为传热微元(i,1)的进口烟气温度，单位℃；t”fg(i,1)为传热微元(i,1)的出口烟气温度，单位℃；ρpa为一次风的密度，单位kg/m³；cpa为一次风的比热容，单位kj/(kg·℃)；t'pa(i,1)和t”pa(i,1)分别为传热微元(i,1)进、出口一次风温度，单位℃；aconv(i,1)为传热微元(i,1)的对流换热面积，单位m²；tm(i,1)和tm(i,2)为传热微元(i,1)左侧金属壁温度、右侧金属壁温度，单位℃；αpa为一次风对流传热系数，单位kw/(m²·℃)。

可选地，根据所述回转式空预器温度分布，实时估计所述回转式空预器的直接漏风流量，包括：

根据一次风的质量守恒和能量守恒，实时估计一次风直接漏风质量流量，一次风质量和能量稳态平衡关系为：

d”pa＝d'pa-dlp

q”pa＝q'pa+qpm-qlp

其中，d”pa为回转式空预器出口一次风的质量流量，单位kg/s；d'pa为回转式空预器入口一次风的质量流量，单位kg/s；dlp为回转式空预器一次风直接漏风的质量流量，单位kg/s；q”pa为回转式空预器出口一次风热量，单位kw；q'pa为回转式空预器入口一次风热量，单位kw；qpm为一次风与金属壁换热量，单位kw；qlp为一次风直接漏风的热量；其中：

q”pa＝d”pah”pa＝(d'pa-dlp)h”pa

q'pa＝d'pah'pa

其中：h”pa为回转式空预器出口一次风的比焓，单位kj/kg；h'pa为回转式空预器入口一次风的比焓，单位kj/kg；αpa为一次风对流传热系数，单位kw/(m²·℃)；t'pa(i,j)和t”pa(i,j)分别是传热微元(i,j)进、出口处的一次风温度，单位℃；m是回转式空预器沿烟气流动方向离散化的份数；

由于转子的蘑菇状变形，直接漏风主要发生在空预器的热端和冷端：

dlp＝dlp,h+dlp,c

其中，dlp,h,dlp,c分别为回转式空预器热端和冷端的一次风直接漏风质量流量，kg/s；则qlp的计算公式如下：

qlp＝dlp,hh”pa+dlp,ch'pa

其中，dlp,h和dlp,c通过求解如下目标函数的最优化问题得到，令

其中：j1为辨识dlp,h和dlp,c的目标函数；

根据二次风的质量守恒和能量守恒，实时估计二次风直接漏风质量流量，二次风质量和能量稳态平衡关系为：

d”sa＝d’sa-dls

q”sa＝q’sa+qsm-qls

其中，d”sa为回转式空预器出口二次风的质量流量，单位kg/s；d'sa为回转式空预器入口二次风的质量流量，单位kg/s；dls为回转式空预器二次风直接漏风的质量流量，单位kg/s；q”sa为回转式空预器出口二次风热量，单位kw；q'sa为回转式空预器入口二次风热量，单位kw；qsm为二次风与金属壁换热量，单位kw；qls为二次风直接漏风的热量，单位kw；其中：

q”sa＝d”sah”sa＝(d’sa-dls)h”sa

q’sa＝d’sah’sa

其中：h”sa为回转式空预器出口二次风的比焓，单位kj/kg；h'sa为回转式空预器入口二次风的比焓，单位kj/kg；αsa为二次风对流传热系数，kw/(m²·℃)；t'sa(i,j)和t”sa(i,j)分别是传热微元(i,j)进、出口处的二次风温度，单位℃；nsa是二次风仓沿圆周方向离散化的份数；

由于转子的蘑菇状变形，直接漏风主要发生在空预器的热端和冷端：

dls＝dls,h+dls,c

其中，dls,h,dls,c分别为回转式空预器热端和冷端的二次风直接漏风质量流量，单位kg/s；则qls的计算公式如下：

qls＝dls,hh”sa+dls,ch’sa

其中，dls,h和dls,c通过求解如下目标函数的最优化问题得到，令

其中：j2为辨识dls,h和dls,c的目标函数。

与现有技术相比，本发明具有如下的有益效果：

本发明提供的回转式空预器直接漏风流量的实时估计方法，通过计算回转式空预器的温度分布；然后根据一次风和二次风短时段稳态质量和能量平衡，估计一次风和二次风的直接漏风质量流量，从而可以实现在线监测回转式空预器漏风状况；通过对回转式空预器的一次风直接漏风和二次风直接漏风进行分离计算，提高直接漏风质量流量的估计精度。

附图说明

通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：

图1为本发明一实施例空预器烟气、空气流向示意图；

图2为本发明一实施例空预器金属壁温度分布示意图；

图3为本发明一实施例流经空预器气体温度分布示意图；

图4为本发明一实施例空预器直接漏风计算结果示意图。

图中：

1-烟气仓；

2-密封区；

3-一次风仓；

4-二次风仓。

具体实施方式

下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。

本发明提供一种回转式空预器直接漏风流量的实时估计方法，包括：

s101、沿气体流动方向和圆周方向，将回转式空预器离散化为多个传热微元。

本实施例中，在沿气体流动方向和圆周方向，将回转式空预器离散化为多个传热微元之前，还可以预先从锅炉集控运行规程中获取回转式空预器的结构参数，所述结构参数包括：转子直径，转子在高温段、中温段、低温段的高度，在高温段、中温段、低温段内传热微元所用的材料，一次风仓、二次风仓、烟气仓所占的角度。从分布式控制系统dcs的实时数据库中，获取给定时刻下所述回转式空预器的当前检测参数，所述当前检测参数包括：烟气入口温度、压力，烟气出口温度、压力，引风机烟气体积流量，一次风入口温度、压力，一次风出口温度、压力，进入磨煤机的一次风体积流量，二次风入口温度、压力，二次风出口温度、压力，送风机的二次风体积流量。

图1为本发明一实施例空预器烟气、空气流向示意图，如图1所示，三分仓回转式空预器包括：烟气仓1、一次风仓3、二次风仓4，烟气和空气之间用密封隔板分隔为密封区2。回转式空预器的漏风分为携带漏风与直接漏风，携带漏风是指在回转式空预器的转子转动过程中，部分驻留在转子内部空间中的空气被携带到下一个分仓。

s102、对各个传热微元建立能量守恒模型，通过所述能量守恒模型得到考虑携带漏风影响的回转式空预器温度分布；其中，所述回转式空预器温度分布包括：回转式空预器的金属壁温度分布、流经回转式空预器的气体温度分布。

本实施例中，计算气体流经所述回转式空预器所需时间，计算公式如下：

其中：hap为回转式空预器高度，单位m；ugs为气体流速，m/s；tgs为气体流经所述回转式空预器所需时间，单位s；

计算所述回转式空预器转子转过的角度，计算公式如下：

其中：ω为回转式空预器转子的转动角速度，单位rad/s；βet为所述回转式空预器转子转过的角度，单位rad；

计算携带漏风体积流量，计算公式如下：

其中：fpa为一次风出口的体积流量，单位m³/s；npa为一次风仓沿圆周方向的离散化份数；βu,pa为一次风仓内一列传热微元对应角度，单位rad；fpa,et为携带漏风体积流量，单位m³/s；

流经烟气仓一列传热微元的烟气体积流量的计算公式如下：

其中：ffg为烟气出口的体积流量，单位m³/s；nfg为烟气仓沿圆周方向的离散化份数；βu,fg为烟气仓内一列传热微元对应角度，单位rad；ffg,1为流经烟气仓第一列传热微元的烟气体积流量，单位m³/s；

对于烟气仓内任意传热微元(i,j)，其中i、j分别是该传热微元所在行数和列数，若j≠1，则根据能量守恒关系，烟气在传热微元的进、出口的能量之差等于烟气与金属的对流换热量，即满足如下公式：

其中：ρfg为烟气的密度，单位kg/m³；cfg为烟气的比热容，单位kj/(kg·℃)；t'fg(i,j)和t”fg(i,j)分别是传热微元(i,j)的进、出口烟气温度，单位℃；αfg为烟气对流传热系数，kw/(m²·℃)；aconv(i,j)为传热微元(i,j)的对流换热面积；tm(i,j)和tm(i,j+1)分别为传热微元(i,j)左侧金属壁温度、右侧金属壁温度，单位℃；

若j＝1，则所述传热微元位于烟气仓内的第一列，根据能量守恒关系，满足如下公式：

其中：t'fg(i,1)为传热微元(i,1)的进口烟气温度，单位℃；t”fg(i,1)为传热微元(i,1)的出口烟气温度，单位℃；ρpa为一次风的密度，单位kg/m³；cpa为一次风的比热容，kj/(kg·℃)；t'pa(i,1)和t”pa(i,1)分别为传热微元(i,1)进、出口一次风温度，单位℃；aconv(i,1)为传热微元(i,1)的对流换热面积，单位m²；tm(i,1)和tm(i,2)为传热微元(i,1)左侧金属壁温度、右侧金属壁温度，单位℃；αpa为一次风对流传热系数，kw/(m²·℃)。

s103、根据所述回转式空预器温度分布，实时估计所述回转式空预器的直接漏风流量。

本实施例中，可以根据一次风的质量守恒和能量守恒，实时估计一次风直接漏风质量流量，一次风质量和能量稳态平衡关系为：

d”pa＝d'pa-dlp

q”pa＝q'pa+qpm-qlp

其中，d”pa为回转式空预器出口一次风的质量流量，kg/s；d'pa为回转式空预器入口一次风的质量流量，单位kg/s；dlp为回转式空预器一次风直接漏风的质量流量，单位kg/s；q”pa为回转式空预器出口一次风热量，单位kw；q'pa为回转式空预器入口一次风热量，kw；qpm为一次风与金属壁换热量，单位kw；qlp为一次风直接漏风的热量，单位kw；其中：

q”pa＝d”pah”pa＝(d'pa-dlp)h”pa

q'pa＝d'pah'pa

其中：h”pa为回转式空预器出口一次风的比焓，单位kj/kg；h'pa为回转式空预器入口一次风的比焓，单位kj/kg；αpa为一次风对流传热系数，单位kw/(m²·℃)；t'pa(i,j)和t”pa(i,j)分别是传热微元(i,j)进、出口处的一次风温度，单位℃；m是回转式空预器沿烟气流动方向离散化的份数；

由于转子的蘑菇状变形，直接漏风主要发生在空预器的热端和冷端：

dlp＝dlp,h+dlp,c

其中，dlp,h,dlp,c分别为回转式空预器一次风热端和冷端的直接漏风质量流量，则qlp的计算公式如下：

qlp＝dlp,hh”pa+dlp,ch'pa

其中，dlp,h和dlp,c通过求解如下目标函数的最优化问题得到，令

其中：j1为辨识dlp,h和dlp,c的目标函数；

进一步地，可以根据二次风的质量守恒和能量守恒，实时估计二次风直接漏风质量流量，二次风质量和能量稳态平衡关系为：

d”sa＝d’sa-dls

q”sa＝q’sa+qsm-qls

其中，d”sa为回转式空预器出口二次风的质量流量，单位kg/s；d'sa为回转式空预器入口二次风的质量流量，单位kg/s；dls为回转式空预器二次风直接漏风的质量流量，单位kg/s；q”sa为回转式空预器出口二次风热量，单位kw；q'sa为回转式空预器入口二次风热量，单位kw；qsm为回转式空预器二次风与金属壁换热量，单位kw；qls为回转式空预器二次风直接漏风的热量，单位kw；其中：

q”sa＝d”sah”sa＝(d’sa-dls)h”sa

q'sa＝d'sah'sa

其中：h”sa为回转式空预器出口二次风的比焓，单位kj/kg；h'sa为回转式空预器入口二次风的比焓，单位kj/kg；αsa为二次风对流传热系数，单位kw/(m²·℃)；t'sa(i,j)和t”sa(i,j)分别是传热微元(i,j)进、出口处的二次风温度，单位℃；nsa是二次风仓沿圆周方向离散化的份数；

由于转子的蘑菇状变形，直接漏风主要发生在空预器的热端和冷端：

dls＝dls,h+dls,c

其中，dls,h,dls,c分别为回转式空预器二次风热端和冷端的直接漏风质量流量，单位kg/s；则qls的计算公式如下：

qls＝dls,hh”sa+dls,ch'sa

其中，dls,h和dls,c通过求解如下目标函数的最优化问题得到，令

其中：j2为辨识dls,h和dls,c的目标函数。

本实施例，通过计算回转式空预器的温度分布；然后根据一次风和二次风短时段稳态质量和能量平衡，估计一次风和二次风的直接漏风质量流量。可以实现在线监测回转式空预器漏风状况，通过对回转式空预器的一次风直接漏风和二次风直接漏风进行分离计算，从而提高直接漏风质量流量的估计精度。

具体地，以某300mw亚临界机组回转式空预器为例，进行详细说明。已知锅炉型号hg-1025/17.3-wm18，为自然循环、一次中间再热汽包炉。该锅炉配备两台空预器，均为三分仓回转式空预器，型号28.5-fi(t)-1780-smr，转速1转/min。空预器分a、b两侧，直接漏风质量流量估计方法相同，下面以a侧空预器为例描述，具体步骤如下：

步骤1：从锅炉的集控运行规程中，获取空预器如下结构参数：空预器转子直径9.97m，转子总高度1.99m，其中高温段、中温段、低温段的高度分别为1.08m、0.56m、0.36m。高温段、中温段传热元件所用材料为碳钢，低温段传热元件所用材料为搪瓷钢。烟气仓、二次风仓和一次风仓的角度为2.88rad、2.01rad、0.61rad。

步骤2：从分布式控制系统(distributedcontrolsystem，dcs)的实时数据库中读取给定时刻下的空预器烟气入口温度、压力，空预器烟气出口温度、压力，引风机烟气体积流量；空预器一次风入口温度、压力，空预器出口一次风温度、压力，进入a、b、c、d4台磨煤机的一次风体积流量；空预器二次风入口温度、压力，空预器出口二次风温度、压力，送风机二次风体积流量。

步骤3：将空预器沿气体流动方向和圆周方向离散化得到大量传热微元，其中烟气仓、二次风仓和一次风仓沿圆周方向离散化为30份、20份、7份，空预器高温段、中温段和低温段均沿气体流动方向离散化为10份。

由于空预器转子转动，每个分仓第一列传热微元内含有流经该分仓的气体以及从上一分仓携带的气体。由于各分仓内气体流速几乎相同，故气体流经空预器所需时间为：

其中：hap为回转式空预器高度，单位m；ugs为气体流速，单位m/s；tgs为气体流经所述回转式空预器所需时间，单位s；

计算所述回转式空预器转子转过的角度，计算公式如下：

其中：ω为回转式空预器转子的转动角速度，单位rad/s；βet为所述回转式空预器转子转过的角度，单位rad；

计算携带漏风体积流量，计算公式如下：

其中：fpa为一次风出口的体积流量，单位m³/s；npa为一次风仓沿圆周方向的离散化份数；βu,pa为一次风仓内一列传热微元对应角度，单位rad；fpa,et为携带漏风体积流量，单位m³/s；

流经烟气仓一列传热微元的烟气体积流量的计算公式如下：

其中：ffg为烟气出口的体积流量，单位m³/s；nfg为烟气仓沿圆周方向的离散化份数；βu,fg为烟气仓内一列传热微元对应角度，单位rad；ffg,1为流经烟气仓第一列传热微元的烟气体积流量，单位m³/s；

对于烟气仓内任意传热微元(i,j)，其中i、j分别是该传热微元所在行数和列数，若j≠1，则根据能量守恒关系，烟气在传热微元的进、出口的能量之差等于烟气与金属的对流换热量，即满足如下公式：

其中：ρfg为烟气的密度，单位kg/m³；cfg为烟气的比热容，单位kj/(kg·℃)；t'fg(i,j)和t”fg(i,j)分别是传热微元(i,j)的进、出口烟气温度，单位℃；αfg为烟气对流传热系数，kw/(m²·℃)；aconv(i,j)为传热微元(i,j)的对流换热面积，单位m²；tm(i,j)和tm(i,j+1)分别为传热微元(i,j)左侧金属壁温度、右侧金属壁温度，单位℃；

若j＝1，则所述传热微元位于烟气仓内的第一列，根据能量守恒关系，满足如下公式：

其中：t'fg(i,1)为传热微元(i,1)的进口烟气温度，单位℃；t”fg(i,1))为传热微元(i,1)的出口烟气温度，单位℃；ρpa为一次风的密度，单位kg/m³；cpa为一次风的比热容，单位kj/(kg·℃)；t'pa(i,1)和t”pa(i,1)分别为传热微元(i,1)进、出口一次风温度，单位℃；

aconv(i,1)为传热微元(i,1)的对流换热面积，单位m²；tm(i,1)和tm(i,2)为传热微元(i,1)左侧金属壁温度、右侧金属壁温度，单位℃；αpa为一次风对流传热系数，kw/(m²·℃)；

考虑到空预器转子转动引起不稳定换热，烟气对流传热系数的计算方法为(参见杨世铭，陶文铨.传热学(第四版).高等教育出版社,2010,259-262)：

其中，k1是烟气仓内烟气对流传热系数的修正因子；λfg是烟气导热系数，单位w/(m·℃)；de是烟气仓的等效直径，单位m；re是雷诺数，pr是普朗特数，二者计算方法如下：

其中，ufg是烟气流速，单位m/s；vfg是烟气黏度，单位pa·s。相应地，分别引入烟气、二次风和一次风对流传热系数的修正因子k1，k2和k3；

传热单元(i,j)在圆周方向相对上一个传热单元的热量增量δq为：

式中，kper是圆周方向温度修正因子；cm是金属比热容，单位kj/(kg·℃)；hu是一个传热微元的高度，单位m；asec,fg是烟气仓的横截面积，单位m²；e是导热波纹板的孔隙率；βfg是烟气仓的角度，单位rad；ρm为金属壁密度，单位kg/m³；

记k＝[kper,k1,k2,k3]，k的辨识通过解算如下目标函数获得：

其中：j为辨识k的目标函数；nsa是二次风仓沿圆周方向离散化的份数；n为回转式空预器沿圆周方向离散化的份数；m是回转式空预器沿烟气流动方向离散化的份数；t”fg(m+1,i)为传热微元(m+1,i)出口烟气温度，单位℃；t”fg,mea为回转式空预器出口烟气温度实测值，单位℃；t”sa(1,i)为传热微元(1,i)出口二次风温度，单位℃；t”sa,mea为回转式空预器出口二次风温度实测值，单位℃；t”pa(1,i)为传热微元(1,i)出口二次风温度，单位℃；t”pa,mea为回转式空预器出口一次风温度实测值，单位℃；

辨识结果k＝[4.881.552.772.39]。将k的值代入，从烟气仓向一次风仓逐列计算气体温度和金属壁温度。当金属壁温度满足如下收敛条件时停止计算气体温度和金属壁温度：

其中：tdev,m为回转式空预器最左侧金属壁与最右侧金属壁温度偏差，单位℃；tm(i,m+1)和tm(i,1)分别为回转式空预器最左侧和最右侧金属壁温度计算值，单位℃。否则，将一次风仓最后一列金属壁温度赋给烟气仓第一列金属壁温度，进行下一轮迭代计算，直到金属壁温度满足收敛条件为止。

步骤4：分别估计一次风、二次风的直接漏风质量流量。

本发明实现了空预器直接漏风质量流量的实时估计，可用于在线监测空预器漏风状况，对回转式空预器安全、经济运行具有潜在应用价值。

以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。