CN104122291B - 超超临界火电机组水冷壁向工质传热速率的实时辨识方法 - Google Patents

Abstract

本发明提供了一种超超临界火电机组水冷壁向工质传热速率的实时辨识方法，该方法具体步骤包括：1：根据锅炉运行设计规程，建立水冷壁的结构参数数据库，并从厂级监控信息系统的实时数据库中读取相关实时数据；2：根据工质物性参数库，结合现场厂级监控信息系统服务器信息，实时计算给定时刻水冷壁内工质的物性参数；3：由水冷壁内工质质量、能量和动量动态衡算模型求得水冷壁各个分段向工质的传热速率，加和后获得水冷壁向工质的总传热速率。水冷壁传热速率是燃烧状况和炉膛内结焦程度的重要表征，可用于燃烧状态监控、低氮燃烧和智能吹灰控制，还可用于入炉煤热值的在线辨识和锅炉效率等的在线估计。

Description

超超临界火电机组水冷壁向工质传热速率的实时辨识方法

技术领域

本发明涉及一种用于超超临界火电机组水冷壁向工质传热速率的实时辨识方法，具体地说，是一种基于水冷壁内工质质量、能量和动量动态衡算的工质传热速率的实时辨识方法，属于火力发电控制领域。

背景技术

燃煤电站水冷壁内工质吸热量的主要监控方法是通过机组负荷和煤-电转换效率反推炉膛总发热量，再将总发热量按一定比例(通常为45-50％)折算为水冷壁内工质的吸热量。

然而，实际电厂生产过程中，最终的煤-电转换效率受煤质、设备状态以及机组负荷等诸多因素影响，时变性很强，所以通过煤-电转换效率反推出的炉膛总发热量通常与实际总发热量有较大偏差；另一方面，水冷壁吸热量占炉膛总发热量的比例也并非固定，而是与燃烧工况、燃煤热值、过量空气系数等操作参数相关，故目前大部分生产现场对水冷壁吸热量只能粗略估计，无法精确地跟踪，这对燃烧状态监控十分不利。

对现有技术的检索发现，中国专利申请号201010553886.4，公开日2011-04-27，涉及了一种锅炉水冷壁对其内部工质传热速率的测量方法，该方法将整个亚临界机组蒸发系统(包括汽包、下降管和水冷壁)内的工质分为饱和水和饱和汽，根据汽包压力和机组负荷的变化实时计算二者的质量比例变化，结合实时计算得工质比焓进一步获得整个蒸发系统内工质蓄积的总能量变化率，最后通过进出口工质焓差和工质蓄积的总能量变化率在线计算获得水冷壁对其内部工质传热速率。但是该方法只针对亚临界机组的蒸发系统进行建模和实时计算，而超超临界机组蒸发系统的结构、工质成分和汽水循环流程与亚临界机组差别较大，导致该方法无法应用于超超临界机组。

发明内容

本发明针对现有技术存在的上述不足，提出了一种基于水冷壁内工质质量、能量和动量平衡的水冷壁工质吸热速率的实时辨识方法。

本发明是通过以下技术方案实现的：本发明基于工质物性参数数据库和蒸发系统设备结构参数数据库，根据工质相态和对象设备具体结构结合水冷壁的实时运行工况信息，获得动态变化的工质物性参数；再利用水冷壁内工质质量、能量和动量动态衡算的方法，对整个水冷壁内工质能量蓄积量变化速率进行实时估计，进而结合水冷壁进出口工质的实时测量值，计算水冷壁进出口焓差，最后通过工质能量蓄积量变化速率和水冷壁进出口焓差在线估计水冷壁向工质的传热速率。

本发明所述的超超临界火电机组水冷壁向工质传热速率的实时辨识方法，具体包括如下步骤：

步骤1：根据锅炉运行设计规程，建立水冷壁的结构参数数据库，所述参数数据库包含：水冷壁换热管道沿工质流动方向的总长度L；换热管数量n；换热管内径d；并从厂级监控信息系统的实时数据库读取当前时刻(定义为k)的水冷壁给水压力p₀、给水温度T₀、给水质量流量D₀，水冷壁出口集箱处过热蒸汽压力p₃、过热蒸汽温度T₃、过热蒸汽质量流量D₃；

步骤2：根据可覆盖全操作工况的工质物性参数数据库，结合现场厂级监控信息系统服务器上传的工质状态参数即温度压力，实时计算该当前时刻(定义为k)水冷壁给水密度ρ₀和比焓h₀、水冷壁出口集箱处过热蒸汽密度ρ₃和比焓h₃；

步骤3：由水冷壁内水相段工质质量、能量、动量动态衡算模型计算得到水相段向工质的传热速率Q_W,由水冷壁内蒸发段工质质量、能量、动量动态衡算模型计算得到蒸发段向工质的传热速率Q_E,由水冷壁内过热段工质质量、能量、动量动态衡算模型计算得到过热段向工质的传热速率Q_S，然后加和得到整个水冷壁向工质的总传热速率Q_T，实现超超临界火电机组水冷壁向工质传热速率的实时辨识。

优选地，以所述给定时刻t₀为计时起点，Δt为计算步长，使用水冷壁向工质传热速率实时辨识方法，按时间先后次序，分别确定t＝t₀、t＝t₀+Δt、t＝t₀+2Δt、……、t＝t₀+nΔt时刻下相应的工质传热速率值Q_T(k)、Q_T(k+Δt)、Q_T(k+2Δt)、……、Q_T(k+nΔt)，并获得水冷壁向工质传热速率Q_T随时间变化的曲线。

优选地，步骤2中：

所述密度ρ₀和比焓h₀均通过当前时刻(定义为k)给水压力p₀和温度T₀计算获得；

所述密度ρ₃和比焓h₃通过当前时刻(定义为k)水冷壁出口集箱处过热蒸汽温度T₃和压力p₃计算获得。

优选地，步骤3中，所述水相段工质质量动态衡算模型具体为：

$A(L_W\frac{{\mathrm{d\ρ}}_W}{dt}+{\mathrm{\ρ}}_W\frac{{\mathrm{dL}}_W}{dt})=D_0-(D_1-{\mathrm{A\ρ}}_1\frac{{\mathrm{dL}}_W}{dt})---\left(1\right)$

其中ρ_W为水相段工质平均密度，ρ_W由水相段两端即水冷壁入口界面和水相段与蒸发段界面处的工质密度的算术平均值替代；A为水冷壁管流通截面积，由水冷壁换热管数量n和管内径d计算获得；L_W为水相段长度；D₀为通过水冷壁入口界面的工质质量流量，即给水质量流量；D₁为通过水相段与蒸发段界面的工质质量流量；ρ₁为水相段与蒸发段界面处工质密度，水相段与蒸发段界面处工质为饱和水，故ρ₁就是对应压力下的饱和水密度；随着机组负荷变化，水相段与蒸发段界面的位置是时变的，故水相段的长度和工质容积也是时变的，式(1)中Aρ₁(dL_W/dt)即为水相段长度变化引起的水相段工质质量的附加变化率；

所述水相段工质能量动态衡算模型具体形式是：

$A(L_W{h}_W\frac{{\mathrm{d\ρ}}_W}{dt}+{\mathrm{\ρ}}_W{h}_W\frac{{\mathrm{dL}}_W}{dt}+{\mathrm{\ρ}}_W{L}_W\frac{{\mathrm{dh}}_W}{dt})=D_0{h}_0-(D_1-{\mathrm{A\ρ}}_1\frac{{\mathrm{dL}}_W}{dt})h+Q_W---\left(2\right)$

其中h_W为水相段工质平均比焓，h_W由水相段两端即水冷壁入口界面和水相段与蒸发段界面处的工质比焓的算术平均值替代；h₀为水相段入口界面处工质比焓；h₁为水相段与蒸发段界面处工质比焓，水相段与蒸发段界面处工质为饱和水，故h₁就是对应压力下的饱和水比焓；Q_W为水冷壁对水相段工质的传热速率；其余变量含义同式(1)；

所述水相段工质动量动态衡算模型具体形式是：

$p_0-p_1=\frac{{\mathrm{\ϵ}}_W{D_W}^2{L}_W}{{\mathrm{\ρ}}_W}---\left(3\right)$

其中ε_W为定值，与水冷壁管道几何结构相关，由国标GB/T2624-93获得；D_W为水相段工质平均质量流量，D_W由水相段两端即水冷壁入口界面和水相段与蒸发段界面处的工质质量流量D₀和D₁的算术平均值替代；p₀为水冷壁入口界面处工质压力；p₁为水相段与蒸发段界面处工质压力；其余变量含义同式(1)；

与水相段类似，所述蒸发段工质的质量、能量、动量动态衡算模型，即公式(4)～(6)：

蒸发段相应变量的下标是E：

$A(L_E\frac{{\mathrm{d\ρ}}_E}{dt}+{\mathrm{\ρ}}_E\frac{{\mathrm{dL}}_E}{dt})=(D_1-{\mathrm{A\ρ}}_1\frac{{\mathrm{dL}}_W}{dt})-\[D_2-{\mathrm{A\ρ}}_2\frac{d(L_E+L_W)}{dt}\]---\left(4\right)$

$A(L_E{h}_E\frac{{\mathrm{d\ρ}}_E}{dt}+{\mathrm{\ρ}}_E{h}_E\frac{{\mathrm{dL}}_E}{dt}+{\mathrm{\ρ}}_E{L}_E\frac{{\mathrm{dh}}_E}{dt})=(D_1-{\mathrm{A\ρ}}_1\frac{{\mathrm{dL}}_W}{dt})h_1-\[D_2-{\mathrm{A\ρ}}_2\frac{d(L_E+L_W)}{dt}\]h_2+Q_E---\left(5\right)$

$p_1-p_2=\frac{{\mathrm{\ϵ}}_E{D_E}^2{L}_E}{{\mathrm{\ρ}}_E}---\left(6\right)$

所述过热段工质的质量、能量、动量动态衡算模型，即公式(7)～(9)：

过热段相应变量的下标是S：

$A(L_S\frac{{\mathrm{d\ρ}}_S}{dt}+{\mathrm{\ρ}}_S\frac{{\mathrm{dL}}_S}{dt})=(D_2-{\mathrm{A\ρ}}_2\frac{{\mathrm{dL}}_S}{dt})-D_3---\left(7\right)$

$A(L_E{h}_E\frac{{\mathrm{d\ρ}}_E}{dt}+{\mathrm{\ρ}}_E{h}_E\frac{{\mathrm{dL}}_E}{dt}+{\mathrm{\ρ}}_E{L}_E\frac{{\mathrm{dh}}_E}{dt})=(D_2-{\mathrm{A\ρ}}_2\frac{{\mathrm{dL}}_S}{dt})h_2-D_3{h}_3+Q_S---\left(8\right)$

$p_2-p_3=\frac{{\mathrm{\ϵ}}_S{D_S}^2{L}_S}{{\mathrm{\ρ}}_S}---\left(9\right)$

式(1)至(9)是线性互不相关的9维方程组，其中水冷壁总长固定为L，故L_W、L_E、L_S中只需求解其中两个即可；水冷壁入口界面(即水冷壁给水)处的工质状态变量(压力、温度和质量流量)由物理测点实时获得；水冷壁出口界面(即水冷壁出口集箱)处工质压力和温度由物理测点实时获得，工质质量流量视为与主蒸汽质量流量相等；由于水相段与蒸发段界面和蒸发段与过热段界面处的工质分别是饱和水和饱和蒸汽，故这两处的工质温度也是压力的函数，通过压力值计算获得；所有分界面处工质物性参数均根据国际水和水蒸气热力性质工业公式1997(IF97)，通过对应分界面处的工质温度、压力加以计算；所有分段内工质物性参数均以两端分界面处物性参数的算术平均值替代。

至此，式(1)至(9)构成的方程组中，待求解变量分别为：L_W、L_E、p₁、p₂、D₁、D₂、Q_W、Q_E、Q_S。

待求解变量数量与方程组维数相同，故解算方程组即获得各段水冷壁向工质的传热速率：Q_W、Q_E、Q_S。进一步，水冷壁向工质的总传热速率Q_T＝Q_W+Q_E+Q_S。

当水冷壁内工质压力超过22Mpa时，将不存在蒸发段，方程组中，式(4)至(6)被剔除；待求解变量中L_E、Q_E被剔除；p₁和p₂合并为p_W/S，D₁和D₂合并为D_W/S；不过，由于工质处于超临界状态，水相段与过热段(W段和S段)分界面处工质的温度不再是压力的函数，而是一个自变量，故待求解变量中需增加分界面处工质温度T_W/S，因此，待求解变量总数为6个：L_W、p_W/S、D_W/S、T_W/S、Q_W、Q_S。

待求解变量数量仍与方程组维数相同，故解算方程组即获得各段水冷壁向工质的传热速率：Q_W、Q_S。则水冷壁向工质的总传热速率Q_T＝Q_W+Q_S。

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

本发明给出了在满足水冷壁内工质质量、能量和动量守恒关系下的水冷壁向工质传热速率的实时辨识值，计算水冷壁向工质传热速率所需的测点数据均从厂级监控信息系统的实时数据库直接读取，现场不需要额外增加测点，仅需在已有的控制系统中增加相应的软件计算模块，成本低。水冷壁向工质传热速率的实时辨识是迭代进行的，迭代中仅需要保存上一步(即k-1时刻)的结果，内存占用少，实时性强。同时，借助可覆盖全工况的工质物性参数数据库，该方法在不同工况和负荷条件下都可以准确跟踪水冷壁向工质传热速率的变化，而设备结构参数数据库的建立则使该方法可以非常方便地推广到同类型但设备结构不同的机组。可用于指导锅炉的燃烧优化，扩展应用范围广。

附图说明

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

图1为本发明水冷壁分段结构示意图，其中0表示水冷壁工质入口界面，1表示水相段和蒸发段分界面，2表示蒸发段与过热段分界面，3表示水冷壁工质出口界面；

图2工质焓值与温度、压力的关系示意图；

图3工质密度与温度、压力的关系示意图。

具体实施方式

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

如图1所示，以某1000MW燃煤机组为例，该直流锅炉烟道中的换热器主要包括过热器、再热器、省煤器三大部分。锅炉型号为DG-3000/26.15-II1型，锅炉为超超临界、直流、一次中间再热、前后墙对冲火焰燃烧方式、双拱单炉膛、平衡通风、尾部双烟道、烟气挡板调温、固态排渣、露天布置、全钢架悬吊式直流炉。

基于水冷壁内工质质量、能量和动量平衡的水冷壁工质吸热速率的实时辨识方法详细步骤如下：

步骤1：根据锅炉运行规程，建立主要设备结构参数数据库，该数据库包含以下固定参数信息：水冷壁换热管道沿工质流动方向的总长度L为38.1m，换热管总数778根，换热管单管内径23mm；

并从厂级监控信息系统的实时数据库中读取相关实时数据，如在当前时刻k的运行工况下实时数据如下：

水冷壁给水压力29.02MPa、给水温度297℃、给水质量流量3033t/h、水冷壁出口集箱处过热蒸汽压力27.4MPa、过热蒸汽温度411℃、过热蒸汽质量流量3033t/h。

本实施例中，所述水冷壁分段结构和各分段之间分界面的编号如图1所示。

步骤2：根据可覆盖全操作工况的工质物性参数数据库，结合现场厂级监控信息系统服务器上传的工质状态参数即温度压力，实时计算该给定时刻下水冷壁给水密度ρ₀和比焓h₀、水冷壁出口集箱处过热蒸汽密度ρ₃和比焓h₃。上述工质物性参数库输出结果如图2、图3所示，其中：图2为工质焓值与温度、压力的关系示意图，图3为工质密度与温度、压力的关系示意图。

本实施例中，所述工质物性参数数据库(即图2、图3)即IF97公式，是指根据水和水蒸气热力性质工业公式(IAPWS-IF97)开发的具有可并行调用的、区域自动判别、批处理运算等特点的用于在线计算的工质物性参数数据库，可参考文献：王旭辉，于彤，惠兆宇，袁景淇，用于火电全范围仿真的工质物性参数数据库，控制工程，2011；18:131-133。

步骤3：由水冷壁内水相段工质质量、能量、动量动态衡算模型计算得到水相段向工质的传热速率Q_W,由水冷壁内蒸发段工质质量、能量、动量动态衡算模型计算得到蒸发段向工质的传热速率Q_E,由水冷壁内过热段工质质量、能量、动量动态衡算模型计算得到过热段向工质的传热速率Q_S，然后加和得到整个水冷壁向工质的总传热速率Q_T。

本实施例中，步骤3中：

所述水相段工质质量动态衡算模型具体为：

$A(L_W\frac{{\mathrm{d\ρ}}_W}{dt}+{\mathrm{\ρ}}_W\frac{{\mathrm{dL}}_W}{dt})=D_0-(D_1-{\mathrm{A\ρ}}_1\frac{{\mathrm{dL}}_W}{dt})---\left(1\right)$

其中ρ_W为水相段工质平均密度，ρ_W由水相段两端即水冷壁入口界面和水相段与蒸发段界面处的工质密度的算术平均值替代；A为水冷壁管流通截面积，由水冷壁换热管数量n和管内径d计算获得；L_W为水相段长度；D₀为通过水冷壁入口界面的工质质量流量，即给水质量流量；D₁为通过水相段与蒸发段界面的工质质量流量；ρ₁为水相段与蒸发段界面处工质密度，水相段与蒸发段界面处工质为饱和水，故ρ₁就是对应压力下的饱和水密度；随着机组负荷变化，水相段与蒸发段界面的位置是时变的，故水相段的长度和工质容积也是时变的，式(1)中Aρ₁(dL_W/dt)即为水相段长度变化引起的水相段工质质量的附加变化率；

所述水相段工质能量动态衡算模型具体形式是：

$A(L_W{h}_W\frac{{\mathrm{d\ρ}}_W}{dt}+{\mathrm{\ρ}}_W{h}_W\frac{{\mathrm{dL}}_W}{dt}+{\mathrm{\ρ}}_W{L}_W\frac{{\mathrm{dh}}_W}{dt})=D_0{h}_0-(D_1-{\mathrm{A\ρ}}_1\frac{{\mathrm{dL}}_W}{dt})h_1+Q_W---\left(2\right)$

其中h_W为水相段工质平均比焓，h_W由水相段两端即水冷壁入口界面和水相段与蒸发段界面处的工质比焓的算术平均值替代；h₀为水相段入口界面处工质比焓；h₁为水相段与蒸发段界面处工质比焓，水相段与蒸发段界面处工质为饱和水，故h₁就是对应压力下的饱和水比焓；Q_W为水冷壁对水相段工质的传热速率；其余变量含义同式(1)；

所述水相段工质动量动态衡算模型具体形式是：

$p_0-p_1=\frac{{\mathrm{\ϵ}}_W{D_W}^2{L}_W}{{\mathrm{\ρ}}_W}---\left(3\right)$

其中ε_W为定值，与水冷壁管道几何结构相关，由国标GB/T2624-93获得；D_W为水相段工质平均质量流量，D_W由水相段两端即水冷壁入口界面和水相段与蒸发段界面处的工质质量流量D₀和D₁的算术平均值替代；p₀为水冷壁入口界面处工质压力；p₁为水相段与蒸发段界面处工质压力；其余变量含义同式(1)；

与水相段类似，所述蒸发段工质的质量、能量、动量动态衡算模型，即公式(4)～(6)：

蒸发段相应变量的下标是E：

$A(L_E\frac{{\mathrm{d\ρ}}_E}{dt}+{\mathrm{\ρ}}_E\frac{{\mathrm{dL}}_E}{dt})=(D_1-{\mathrm{A\ρ}}_1\frac{{\mathrm{dL}}_W}{dt})-\[D_2-{\mathrm{A\ρ}}_2\frac{d(L_E+L_W)}{dt}\]---\left(4\right)$

$A(L_E{h}_E\frac{{\mathrm{d\ρ}}_E}{dt}+{\mathrm{\ρ}}_E{h}_E\frac{{\mathrm{dL}}_E}{dt}+{\mathrm{\ρ}}_E{L}_E\frac{{\mathrm{dh}}_E}{dt})=(D_1-{\mathrm{A\ρ}}_1\frac{{\mathrm{dL}}_W}{dt})h_1-\[D_2-{\mathrm{A\ρ}}_2\frac{d(L_E+L_W)}{dt}\]h_2+Q_E---\left(5\right)$

$p_1-p_2=\frac{{\mathrm{\ϵ}}_E{D_E}^2{L}_E}{{\mathrm{\ρ}}_E}---\left(6\right)$

所述过热段工质的质量、能量、动量动态衡算模型，即公式(7)～(9)：

过热段相应变量的下标是S：

$A(L_S\frac{{\mathrm{d\ρ}}_S}{dt}+{\mathrm{\ρ}}_S\frac{{\mathrm{dL}}_S}{dt})=(D_2-{\mathrm{A\ρ}}_2\frac{{\mathrm{dL}}_S}{dt})-D_3---\left(7\right)$

$A(L_E{h}_E\frac{{\mathrm{d\ρ}}_E}{dt}+{\mathrm{\ρ}}_E{h}_E\frac{{\mathrm{dL}}_E}{dt}+{\mathrm{\ρ}}_E{L}_E\frac{{\mathrm{dh}}_E}{dt})=(D_2-{\mathrm{A\ρ}}_2\frac{{\mathrm{dL}}_S}{dt})h_2-D_3{h}_3+Q_S---\left(8\right)$

$p_2-p_3=\frac{{\mathrm{\ϵ}}_S{D_S}^2{L}_S}{{\mathrm{\ρ}}_S}---\left(9\right)$

(1)至(9)是线性互不相关的9维方程组，其中水冷壁总长固定为L，故L_W、L_E、L_S中只需求解其中两个即可；水冷壁入口界面(即水冷壁给水)处的工质状态变量(压力、温度和质量流量)由物理测点实时获得；水冷壁出口界面(即水冷壁出口集箱)处工质压力和温度由物理测点实时获得，工质质量流量视为与主蒸汽质量流量相等；由于水相段与蒸发段界面和蒸发段与过热段界面处的工质分别是饱和水和饱和蒸汽，故这两处的工质温度也是压力的函数，可通过压力值计算获得；所有分界面处工质物性参数均根据国际水和水蒸气热力性质工业公式1997(IF97)，通过对应分界面处的工质温度、压力加以计算；所有分段内工质物性参数均以两端分界面处物性参数的算术平均值替代。

至此，式(1)至(9)构成的方程组中，待求解变量分别为：L_W、L_E、p₁、p₂、D₁、D₂、Q_W、Q_E、Q_S。

待求解变量数量与方程组维数相同，故解算方程组即可获得各段水冷壁向工质的传热速率：Q_W、Q_E、Q_S。进一步，水冷壁向工质的总传热速率Q_T＝Q_W+Q_E+Q_S。

当水冷壁内工质压力超过22Mpa时，将不存在蒸发段，方程组中，式(4)至(6)被剔除；待求解变量中L_E、Q_E被剔除；p₁和p₂合并为p_W/S，D₁和D₂合并为D_W/S；不过，由于工质处于超临界状态，水相段与过热段(W段和S段)分界面处工质的温度不再是压力的函数，而是一个自变量，故待求解变量中需增加分界面处工质温度T_W/S，因此，待求解变量总数为6个：L_W、p_W/S、D_W/S、T_W/S、Q_W、Q_S。

待求解变量数量仍与方程组维数相同，故解算方程组即可获得各段水冷壁向工质的传热速率：Q_W、Q_S。则水冷壁向工质的总传热速率Q_T＝Q_W+Q_S。

本实施例选取时间间隔Δt＝5s，计算t＝t₀为计时起点的24小时水冷壁向工质传热速率，即在时刻t＝t₀、t＝t₀+5、t＝t₀+10、……、t＝t₀+86400，分别重复步骤1-3，获得相应的水冷壁向工质传热速率实时辨识Q_T(t₀)、Q_T(t₀+5)、Q_T(t₀+10)、……、Q_T(t₀+86400)，继而得到水冷壁向工质传热速率随时间变化曲线。

以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变形或修改，这并不影响本发明的实质内容。

Claims (2)

1.一种超超临界火电机组水冷壁向工质传热速率的实时辨识方法，其特征在于，包括如下步骤：

步骤1：根据锅炉运行设计规程，建立水冷壁的结构参数数据库，所述参数数据库包含：水冷壁换热管道沿工质流动方向的总长度L，换热管数量n，换热管内径d；从厂级监控信息系统的实时数据库读取当前时刻k的水冷壁给水压力p₀、给水温度T₀、给水质量流量D₀，水冷壁出口集箱处过热蒸汽压力p₃、过热蒸汽温度t₃、过热蒸汽质量流量D₃；

步骤2：根据可覆盖全操作工况的工质物性参数数据库，结合现场厂级监控信息系统实时提供的工质状态即温度、压力，实时计算该给定时刻下水冷壁给水密度ρ₀和比焓h₀，以及水冷壁出口集箱处过热蒸汽密度ρ₃和比焓h₃；

步骤3：由水冷壁内水相段工质质量、能量、动量动态衡算模型计算得到水相段向工质的传热速率Q_W,由水冷壁内蒸发段工质质量、能量、动量动态衡算模型计算得到蒸发段向工质的传热速率Q_E,由水冷壁内过热段工质质量、能量、动量动态衡算模型计算得到过热段向工质的传热速率Q_S，然后加和得到整个水冷壁向工质的总传热速率Q_T，实现超超临界火电机组水冷壁向工质传热速率的实时辨识；

所述水相段工质质量动态衡算模型具体为：

$A(L_W\frac{{\mathrm{d\ρ}}_W}{dt}+{\mathrm{\ρ}}_W\frac{{\mathrm{dL}}_W}{dt})=D_0-(D_1-{\mathrm{A\ρ}}_1\frac{{\mathrm{dL}}_W}{dt})---\left(1\right)$

其中ρ_W为水相段工质平均密度，ρ_W由水相段两端即水冷壁入口界面和水相段与蒸发段界面处的工质密度的算术平均值替代；A为水冷壁管流通截面积，由水冷壁换热管数量n和管内径d计算获得；L_W为水相段长度；D₀为通过水冷壁入口界面的工质质量流量，即给水质量流量；D₁为通过水相段与蒸发段界面的工质质量流量；ρ₁为水相段与蒸发段界面处工质密度，水相段与蒸发段界面处工质为饱和水，故ρ₁就是对应压力下的饱和水密度；随着机组负荷变化，水相段与蒸发段界面的位置是时变的，故水相段的长度和工质容积也是时变的，式(1)中Aρ₁(dL_W/dt)即为水相段长度变化引起的水相段工质质量的附加变化率；

所述水相段工质能量动态衡算模型具体形式是：

$A(L_W{h}_W\frac{{\mathrm{d\ρ}}_W}{dt}+{\mathrm{\ρ}}_W{h}_W\frac{{\mathrm{dL}}_W}{dt}+{\mathrm{\ρ}}_W{L}_W\frac{{\mathrm{dh}}_W}{dt})=D_0{h}_0-(D_1-{\mathrm{A\ρ}}_1\frac{{\mathrm{dL}}_W}{dt})h_1+Q_W---\left(2\right)$

其中h_W为水相段工质平均比焓，h_W由水相段两端即水冷壁入口界面和水相段与蒸发段界面处的工质比焓的算术平均值替代；h₀为水相段入口界面处工质比焓；h₁为水相段与蒸发段界面处工质比焓，水相段与蒸发段界面处工质为饱和水，故h₁就是对应压力下的饱和水比焓；Q_W为水冷壁对水相段工质的传热速率；其余变量含义同式(1)；

所述水相段工质动量动态衡算模型具体形式是：

$p_0-p_1=\frac{{\mathrm{\ϵ}}_W{D_W}^2{L}_W}{{\mathrm{\ρ}}_W}---\left(3\right)$

其中ε_W为定值，与水冷壁管道几何结构相关，由国标GB/T2624-93获得；D_W为水相段工质平均质量流量，D_W由水相段两端即水冷壁入口界面和水相段与蒸发段界面处的工质质量流量D₀和D₁的算术平均值替代；p₀为水冷壁入口界面处工质压力；p₁为水相段与蒸发段界面处工质压力；其余变量含义同式(1)；

与水相段类似，所述蒸发段工质的质量、能量、动量动态衡算模型，即公式(4)～(6)：

蒸发段相应变量的下标是E：

$A(L_E\frac{{\mathrm{d\ρ}}_E}{dt}+{\mathrm{\ρ}}_E\frac{{\mathrm{dL}}_E}{dt})=(D_1-{\mathrm{A\ρ}}_1\frac{{\mathrm{dL}}_W}{dt})-\[D_2-{\mathrm{A\ρ}}_2\frac{d(L_E+L_W)}{dt}\]---\left(4\right)$

$A(L_E{h}_E\frac{{\mathrm{d\ρ}}_E}{dt}+{\mathrm{\ρ}}_E{h}_E\frac{{\mathrm{dL}}_E}{dt}+{\mathrm{\ρ}}_E{L}_E\frac{{\mathrm{dh}}_E}{dt})=(D_1-{\mathrm{A\ρ}}_1\frac{{\mathrm{dL}}_W}{dt})h_1-\[D_2-{\mathrm{A\ρ}}_2\frac{d(L_E+L_W)}{dt}\]h_2+Q_E---\left(5\right)$

$p_1-p_2=\frac{{\mathrm{\ϵ}}_E{D_E}^2{L}_E}{{\mathrm{\ρ}}_E}---\left(6\right)$

所述过热段工质的质量、能量、动量动态衡算模型，即公式(7)～(9)：

过热段相应变量的下标是S：

$A(L_S\frac{{\mathrm{d\ρ}}_S}{dt}+{\mathrm{\ρ}}_S\frac{{\mathrm{dL}}_S}{dt})=(D_2-{\mathrm{A\ρ}}_2\frac{{\mathrm{dL}}_S}{dt})-D_3---\left(7\right)$

$A(L_E{h}_E\frac{{\mathrm{d\ρ}}_E}{dt}+{\mathrm{\ρ}}_E{h}_E\frac{{\mathrm{dL}}_E}{dt}+{\mathrm{\ρ}}_E{L}_E\frac{{\mathrm{dh}}_E}{dt})=(D_2-{\mathrm{A\ρ}}_2\frac{{\mathrm{dL}}_S}{dt})h_2-D_3{h}_3+Q_S---\left(8\right)$

$p_2-p_3=\frac{{\mathrm{\ϵ}}_S{D_S}^2{L}_S}{{\mathrm{\ρ}}_S}---\left(9\right)$

式(1)至(9)是线性互不相关的9维方程组，其中水冷壁总长固定为L，故L_W、L_E、L_S中只需求解其中两个即可；水冷壁入口界面即水冷壁给水处的工质状态变量由物理测点实时获得；水冷壁出口界面即水冷壁出口集箱处工质压力和温度由物理测点实时获得，工质质量流量视为与主蒸汽质量流量相等；由于水相段与蒸发段界面和蒸发段与过热段界面处的工质分别是饱和水和饱和蒸汽，故这两处的工质温度也是压力的函数，通过压力值计算获得；所有分界面处工质物性参数均根据国际水和水蒸气热力性质工业公式1997(IF97)，通过对应分界面处的工质温度、压力加以计算；所有分段内工质物性参数均以两端分界面处物性参数的算术平均值替代；

至此，式(1)至(9)构成的方程组中，待求解变量分别为：L_W、L_E、p₁、p₂、D₁、D₂、Q_W、Q_E、Q_S；

待求解变量数量与方程组维数相同，故解算方程组即获得各段水冷壁向工质的传热速率：Q_W、Q_E、Q_S，水冷壁向工质的总传热速率Q_T＝Q_W+Q_E+Q_S；

当水冷壁内工质压力超过22Mpa时，将不存在蒸发段，方程组中，式(4)至(6)被剔除；待求解变量中L_E、Q_E被剔除；p₁和p₂合并为p_W/S，D₁和D₂合并为D_W/S；不过，由于工质处于超临界状态，水相段与过热段分界面即W段和S段处工质的温度不再是压力的函数，而是一个自变量，故待求解变量中需增加分界面处工质温度T_W/S，因此，待求解变量总数为6个：L_W、p_W/S、D_W/S、T_W/S、Q_W、Q_S；

待求解变量数量仍与方程组维数相同，故解算方程组即获得各段水冷壁向工质的传热速率：Q_W、Q_S，则水冷壁向工质的总传热速率Q_T＝Q_W+Q_S。

2.根据权利要求1所述的超超临界火电机组水冷壁向工质传热速率的实时辨识方法，其特征在于，步骤2中：

所述密度ρ₀和比焓h₀均由当前时刻k给水压力p₀和温度T₀计算获得；

所述密度ρ₃和比焓h₃由当当前时刻k水冷壁出口集箱处过热蒸汽温度T₃和压力p₃计算获得。