CN103793746A - 一种燃煤电站锅炉过热器模型参数辨识方法 - Google Patents

Abstract

本发明公开了一种燃煤电站锅炉过热器模型参数辨识方法，包括以下步骤：建立锅炉过热器的非线性集总参数模型，确定已知参数和需要辨识的参数；对已知参数进行坏点处理和数据平滑处理，得到已知参数实时数据库；建立n个不同负荷段的训练样本数据库，并对参数进行归一化处理；针对每个训练数据库建立一个对应的RBF神经网络模型，并使n个RBF神经网络模型并联组成混合网络；提取当前时刻过热器***的实际测量值，对RBF神经网络模型进行在线参数辨识；每隔时间间隔t，对RBF神经网络模型的进行更新处理。本发明避免了常规过热器模型中采用固定参数的缺点，实现了模型参数的实时辨识，保证了模型参数的辨识精度。

Description

一种燃煤电站锅炉过热器模型参数辨识方法

技术领域

本发明涉及信息处理技术领域，尤其涉及一种燃煤电站锅炉过热器模型参数辨识方法。

背景技术

对锅炉建立数学模型是得到其动态特性的常用方法，模型建立的准确度与锅炉的整体优化控制息息相关。随着锅炉机组容量的增加，所需检测和控制参数的增多，使得锅炉结构也逐渐变得复杂，给模型的建立和参数的确定带来了新的挑战。

过热器作为锅炉的重要组成部分，其输出变量主蒸汽温度和主蒸汽压力是***的重要监控参量，因此对其建立精确的数学模型是对上述关键参量进行优化控制的前提。而模型的精确程度又与模型中参数的选取精度直接相关，故提高模型中的参数辨识精度具有重要意义。现有过热器的数学模型大多数为非线性的集总参数模型，模型中的参数大多是通过对实际***进行不同工况的测试试验获得的，所得参数只能反映有限试验点附近的***运行状况，难以适应***大范围变负荷运行时的参数变化以及实际运行中各种扰动的影响。因此，如何根据实际***的运行情况对模型的参数进行准确辨识并能够定时更新具有重要的实用价值。

经对现有文献检索发现，《中国电机工程学报》上的“锅炉单相受热面链式建模新方法”一文中提出了一种集总参数链式建模新方法，对于锅炉省煤器、过热器、再热器等单相受热管的典型分布参数对象，在多段集总参数模型的基础上，提出了集总参数链式建模新方法，并建立了单相受热管的链式结构的集总参数数学模型。该方法建模简单，计算工作量小，物理意义明确。但是该技术所得到的参数为基于设计参数计算所得，在实际机组运行过程中，难免出现一定的偏差，精确度有所不足，缺乏一定的泛化能力。

又经检索发现，《热能动力工程》上的“锅炉过热器***的动态仿真模型”基于机理模型很难再现锅炉过热器***的复杂动态特性，利用机理模型为主导，动态神经网络在线修正，显著地提高了仿真模型的精度。其仿真结果表明，此建模方法为复杂大***的动态仿真提供了一个理想的建模途径。但是该技术中提到的模型复杂程度不够，对对象的描述不够精确。

因此，本领域的技术人员致力于开发一种针对锅炉过热器模型的参数辨识方法，能够实现在线优化辨别锅炉***的实时变化情况。

发明内容

有鉴于现有技术的上述缺陷，本发明所要解决的技术问题是提供一种基于数据驱动的燃煤电站过热器模型参数辨识方法，通过机组中和模型参数相关的实际运行数据，建立与所需辨识参数之间的神经网络辨识模型，从而能够根据锅炉***的实时变化情况进行模型参数的在线优化辨识，具有数据容易获取、对机组正常运行影响小等优点，且避免了模型中采用固定参数的缺点，使得模型能够较好地跟踪实际***的运行情况，为***的控制优化打下坚实的基础。

为实现上述目的，本发明提供了一种燃煤电站锅炉过热器模型参数辨识方法，包括以下步骤：

a）建立锅炉过热器的非线性集总参数模型，确定已知参数和需要辨识的参数；

b）对步骤a）中所述已知参数依次进行坏点处理和数据平滑处理，对处理后的已知参数建立已知参数实时数据库；

c）对步骤b）中所述已知参数实时数据库进行n段划分处理，建立n个不同负荷段的训练样本数据库，将每个负荷段的训练样本数据库按照负荷区间均匀的分为若干训练样本子集，并对每个训练样本子集中的参数进行归一化处理；

其中，n≥1。

d）针对步骤c）中的每个所述训练样本数据库建立一个对应的径向基函数神经网络模型，并使n个径向基函数（Radial Basis Function，RBF）神经网络模型并联组成混合网络；

e）提取当前时刻过热器***的实际测量值，对径向基函数神经网络模型进行在线参数辨识，得到当前时刻需要辨识的参数的数值；

f）每隔时间间隔t，对径向基函数神经网络模型进行更新处理，由新的混合网络对需要辨识的参数进行在线参数辨识。

进一步地，所述步骤b）中，对所述已知参数进行坏点处理是指：通过多项式滑动拟合方法来判断坏点，并采用七点二阶前推差分算式对坏点进行剔除。

进一步地，所述步骤b）中，对所述已知参数进行数据平滑处理是采用七点加权滤波法实现。

进一步地，所述步骤c）中，所述已知参数实时数据库进行n段划分处理是指：确定分段个数为n，在保证每段数据库个数相等或近似相等的情况下，将已知参数实时数据库按照负荷大小分为n段，得到n个不同负荷段的训练样本数据库。

进一步地，所述步骤c）中，所述归一化处理是指：通过计算获得训练样本子集中的参数的归一化值。所述归一化值的计算如下：

$X^{*}=\frac{X-X_{\min }}{X_{\max }-X_{\min }}$

其中：为训练样本子集中的参数，X_max是X所在训练样本子集中的最大参数，X_min是X所在训练样本子集中的最小参数，X^*是X的归一化值。

进一步地，所述步骤d）中，所述建立一个对应的径向基函数神经网络模型包括以下步骤：

第一步、确定径向基函数神经网络模型的输出量；

第二步、从所述锅炉过热器的非线性集总参数模型中查找与所述输出量有关系的已知参数，选取径向基函数神经网络模型的输入量；

第三步、采用径向基函数神经网络训练方法对所述输入量和所述输出量进行训练，得到径向基函数神经网络模型。

其中，第一步中的输出量为需要辨识的参数。

其中，第三步中，从已知参数中去除已知参数中非独立参量和难以获取的参量，即可获得径向基函数神经网络模型的输入量。

进一步地，所述步骤e）中，所述在线参数辨识是指：对当前时刻的锅炉过热器的实际测量值依次进行坏点处理、数据平滑处理和归一化处理，并根据当前负荷的大小将处理的实际测量值输入相应负荷段的径向基函数神经网络模型中，从而得到当前时刻需要辨识的参数的数值。

进一步地，所述步骤f）中，所述更新处理，包括以下步骤：

第一步、将过去t时间内的过热器***的实际测量数据依次进行坏点处理、数据平滑处理和归一化处理，得到过去t时间内的处理后的实际测量数据；

第二步、根据过去t时间内每一时刻的负荷，得到过去t时间内每一时刻处理后的实际测量数据所属的训练样本数据库，并进一步得到过去t时间内每一时刻处理后的实际测量数据所属的训练样本数据子集；

第三步、将属于同一训练样本数据子集的处理后的实际测量数据按从前向后的顺序，依次替换原有的样本数据，得到更新后的训练样本数据库；

第四步、针对更新后的训练样本数据库建立一个对应的新的径向基函数神经网络模型。

在本发明的一个较佳实施例中，一种燃煤电站锅炉过热器模型参数辨识方法，包括以下步骤：建立锅炉过热器的非线性集总参数模型，确定已知参数和需要辨识的参数；对已知参数进行坏点处理和数据平滑处理，得到已知参数实时数据库；建立n个不同负荷段的训练样本数据库，并对参数进行归一化处理；针对每个训练数据库建立一个对应的RBF神经网络模型，并使n个RBF神经网络模型并联组成混合网络；提取当前时刻过热器***的实际测量值，对RBF神经网络模型进行在线参数辨识；每隔时间间隔t，对RBF神经网络模型的进行更新处理。

由此可见，与现有技术相比，本发明提出的过热器模型的参数辨识方法基于数据驱动，具有数据容易获取，无需对机组进行试验等优点。本发明通过将负荷分段建立由多个RBF神经网络模型组合的混合神经网络进行参数辨识，减少了单个神经网络泛化能力的不足，提高了过热器模型参数在不同负荷下的辨识精度；在各个负荷段内按照步长进一步划分多个训练样本数据子集，保证了数据样本选择的均衡性，进一步提高了模型的泛化能力。此外，本发明建立了神经网络辨识模型的在线定时更新机制，以反映过热器***在长期运行中其模型参数的变化情况，为模型输出对实际***输出的长期在线跟踪提供了良好的基础。本发明避免了常规过热器模型中采用固定参数的缺点，实现了模型参数的实时辨识，保证了模型参数的辨识精度。

具体实施方式

在本发明的一个实施例中，锅炉过热器满足下面的条件：

A、外部热流沿管长及圆周均匀分布；

B、管壁金属只考虑径向传热，不考虑轴向传热；

C、管内工质不可压缩，且流动截面上参数均匀一致；

D、忽略减温器动态特性；

E、以过热器出口参数为集总参数模型的代表参数,阻力集中在入口。

针对本实施例的锅炉过热器，一种燃煤电站锅炉过热器模型参数辨识方法，包括以下步骤：

第一步，建立锅炉过热器的非线性集总参数模型，确定已知参数和需要辨识的参数。

本实施例的非线性集总参数模型，具体是：

$D_1+D_{\mathrm{jw}}-D_2=V_g\frac{d{\mathrm{\ρ}}_2}{\mathrm{dt}}$

$\begin{array}{cc}{D}_1{H}_1+D_{\mathrm{jw}}{H}_{\mathrm{jw}}-D_2{H}_2+Q_2=V_g\frac{d\left({\mathrm{\ρ}}_2{H}_2\right)}{\mathrm{dt}}& Q_s-Q_2=C_m{M}_g\frac{d{T}_g}{\mathrm{dt}}\end{array}$

$Q_2={\mathrm{KD}}_2^{\mathrm{\α}}(T_g-\frac{T_1+T_2}{2})$

$D_1={{\mathrm{kP}}_1}^{\mathrm{\β}}\sqrt{P_1-P_2}$

其中：D₁是过热器入口蒸汽流量；D₂是过热器出口蒸汽流量；D_jw是减温水流量；H₁是过热器入口蒸汽焓；H₂是过热器出口蒸汽焓；H_jw是减温水焓；V_g是过热器容积；ρ₂是过热器出口蒸汽密度；Q₂是过热器金属管壁传给蒸汽工质热量；Q_s是烟气传给过热器金属管壁热量；C_m是过热器金属密度；M_g是过热器金属质量；T_g是过热器管壁温度；K是换热系数；α是指数；T₁是过热器入口蒸汽温度；T₂是过热器出口蒸汽温度；k是系数；β是指数；P₁是过热器入口蒸汽压力；P₂是过热器出口蒸汽压力。

本实施例中需要辨识的参数包括：传热系数K、系数k及指数β；其余已知参数从设计资料获取，或由现场的测量仪器测得。

本实施例中：额定蒸汽流量D₂是909.6t/h，额定减温水流量D_jw是15t/h，过热器容积V_g是120m³，过热器金属密度C_m是7900kg/m³，指数C_m是0.8，金属比热容是460J/(kg.℃)，额定蒸汽压力C_m是14.5Mpa。

第二步，对已知参数进行坏点处理和数据平滑处理，得到处理后的已知参数，从而得到已知参数实时数据库。

所述的坏点处理是通过多项式滑动拟合方法来判断坏点，并采用七点二阶前推差分算式对坏点进行剔除。

所述的七点二阶前推差分算式，具体是：当第k个点满足下面公式时，则将第k个点剔除：

$|y_k-{\hat{y}}_k|>2.2\sqrt{\underset{k=i-6}{\overset{i}{\mathrm{\Σ}}}{(y_k-{\hat{y}}_k)}^2/6},$

其中： $\left\{\begin{array}{c}{\hat{y}}_1=({32y}_1+{15y}_2+{3y}_3-{4y}_4-{6y}_5-{3y}_6+{5y}_7)/42\\ {\hat{y}}_2=({5y}_1+{4y}_2+{3y}_3+{2y}_4+y_5-y_7)/14\\ {\hat{y}}_3=(y_1+{3y}_2+{4y}_3+{4y}_4+{3y}_5+y_6-{2y}_7)/14\\ {\hat{y}}_4=(-{2y}_1+{3y}_2+{6y}_3+{7y}_4+{6y}_5+{3y}_6-{2y}_7)/21\\ {\hat{y}}_5=(-{2y}_1+y_2+{3y}_3+{4y}_4+{4y}_5+{3y}_6-{2y}_7)/14\\ {\hat{y}}_6=(-y_1+y_2+{2y}_4+{3y}_5+{4y}_6+{5y}_7)/14\\ {\hat{y}}_i=({5y}_{i-1}-{3y}_{i-5}-{6y}_{i-4}-{4y}_{i-3}+{3y}_{i-2}+{15y}_{i-1}+{32y}_i)/42\end{array}\right.,$

i=7,8...，y_i为实测数据，为插值数据，

所述的数据平滑处理是采用七点加权滤波法实现的，具体是：

${\hat{y}}_m=({0.025y}_{m-3}+{0.05y}_{m-2}+{0.075y}_{m-1}+0.7y_m+{0.075y}_{m+1}+0.05y_{m+2}+0.025y_{m+3},$

其中：

为滤波后的数值，y_m为时刻m的实际测量值。

为了后面计算的需要，本实施例中由模型公式：

$D_1+D_{\mathrm{jw}}-D_2=V_g\frac{d{\mathrm{\ρ}}_2}{\mathrm{dt}}$

$D_1={{\mathrm{kP}}_1}^{\mathrm{\β}}\sqrt{P_1-P_2}$

得： $D_1+D_{\mathrm{jw}}-D_2=V_g(\frac{{\∂\mathrm{\ρ}}_2}{{\∂p}_2}\frac{{\mathrm{d\ρ}}_2}{\mathrm{dt}}+\frac{{\∂\mathrm{\ρ}}_2}{{\∂T}_2}\frac{{\mathrm{dT}}_2}{\mathrm{dt}})$

将其变换为： $D_1=D_2-D_{\mathrm{jw}}+V_g(\frac{{\∂\mathrm{\ρ}}_2}{\∂p_2}\frac{{\mathrm{d\ρ}}_2}{\mathrm{dt}}+\frac{{\∂\mathrm{\ρ}}_2}{{\∂T}_2}\frac{{\mathrm{dT}}_2}{\mathrm{dt}})$

则由上式及运行数据可以求得过热器入口蒸汽流量，为了便于辨识，将过热器入口蒸汽流量公式两边取对数，得：

$\log \left(D_1\right)=\log \left(k\right)+\mathrm{\β}\log \left(P_1\right)+\log \sqrt{P_1-P_2}$

整理得： $\log \left(D_1\right)-\log \left(\sqrt{P_1-P_2}\right)=\log \left(k\right)+\mathrm{\β}\log \left(P_1\right)$

令

x=log(P₁)，a=log(k)，β=b则有式：

y=ax+b。

由上述讨论，通过实际运行数据可以求得(x,y)数据对，再利用最小二乘法即可辨识出a和b，进而求得k和β，参数K同理可得。

第三步，对已知参数实时数据库进行n段划分处理，建立n个不同负荷段的训练样本数据库，并将每个负荷段的训练样本数据库按照负荷区间S均匀的分为若干训练样本子集，并对每个训练样本子集中的参数进行归一化处理。

所述的n段划分处理，是指：确定分段个数为n，在保证每段数据库个数尽量相等的情况下，将已知参数实时数据库按照负荷大小分为n段，从而得到n个不同负荷段的训练样本数据库。

本实施例中n=4，分为85%以上负荷段的训练样本数据库，70%～85%负荷段的训练样本数据库，55%～70%负荷段的训练样本数据库，55%以下负荷段的训练样本数据库。

所述归一化处理,是：

$X^{*}=\frac{X-X_{\min }}{X_{\max }-X_{\min }}$

其中：为训练样本子集中的参数，X_max是X所在训练样本子集中的最大参数，X_min是X所在训练样本子集中的最小参数，X^*是X的归一化值。

本实施例进一步将每个负荷段的训练样本数据库按照0.5%的总负荷为步长均匀的分为若干训练样本子集。

第四步，针对每个训练数据库建立一个对应的RBF神经网络模型，并使n个RBF神经网络模型并联组成混合网络。

所述的建立一个对应的RBF神经网络模型，包括以下步骤：

1）确定RBF神经网络模型的输出量，即需要辨识的参数；

2）从锅炉过热器的非线性集总参数模型中查找与需要辨识的参数有关系的已知参数，去除已知参数中非独立和难以获得的量，从而得到RBF神经网络模型的输入量；

3）采用现有的RBF神经网络训练方法对输入量和输出量进行训练，得到RBF神经网络模型。

本实施例中RBF神经网络模型的输入量为：P₁、T_g、T₁、D₁和D₂，输出量为：传热系数K、系数k及指数β。

第五步，提取当前时刻过热器***的实际测量值，对RBF神经网络模型进行在线参数辨识，得到此时需要辨识的参数的数值。

所述的在线参数辨识，是：对当前时刻过热器***的实际测量值依次进行坏点处理、数据平滑处理和归一化处理，并根据当前负荷的大小将处理的实际测量值输入相应负荷段的RBF神经网络模型中，从而得到此时需要辨识的参数的数值。

第六步，每隔6小时，对RBF神经网络模型进行更新处理，得到更新后的RBF神经网络模型，并使更新后的RBF神经网络模型并联组成新的混合网络，由新的混合网络对需要辨识的参数进行在线参数辨识。

所述的更新处理，包括以下步骤：

1）将过去6小时内的过热器***的实际测量数据依次进行坏点处理、数据平滑处理和归一化处理，得到过去t时间内的处理后的实际测量数据；

2）根据过去6小时内每一时刻的负荷，得到过去6小时内每一时刻处理后的实际测量数据所属的训练样本数据库，并进一步得到过去6小时内每一时刻处理后的实际测量数据所属的训练样本数据子集；

3）将属于同一训练样本数据子集的处理后的实际测量数据按从前向后的顺序，依次替换原有的样本数据，得到更新后的训练样本数据库；

4）针对更新后的训练样本数据库建立一个对应的新的RBF神经网络模型。

本实施例为了能够准确地辨识出不同负荷工况下模型参数的变化，根据***的负荷进行分段，将差别较大的工况分别建立多个相应的神经网络辨识模型。按分段建立参数辨识模型，一方面可以减少单个神经网络模拟整个***大范围工况运行时泛化能力的不足，从而提高模型参数在不同负荷下的辨识精度；另一方面又减少了每个神经网络的规模和样本数量，提高了网络训练速度，为在线辨识的实时性提供了相应保证。在建立各个神经网络模型前，需要对实际采集的运行数据进行适当的选择和预处理，以进一步提高神经网络的泛化能力和辨识精度。此外，为每个神经网络辨识模型建立了在线更新训练样本的机制，并定时用实测数据替换训练样本库中原有数据进行神经网络的训练，以更好地反映在时间统计意义下模型参数的变化情况，从而使得模型与实际***能够相匹配，为***的控制与优化提供基础。

以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。

Claims (8)

1.一种燃煤电站锅炉过热器模型参数辨识方法，其特征在于，包括以下步骤：

a）建立锅炉过热器的非线性集总参数模型，确定已知参数和需要辨识的参数；

b）对步骤a）中所述已知参数依次进行坏点处理和数据平滑处理，对处理后的已知参数建立已知参数实时数据库；

c）对步骤b）中所述已知参数实时数据库进行n段划分处理，建立n个不同负荷段的训练样本数据库，将每个负荷段的训练样本数据库按照负荷区间均匀的分为若干训练样本子集，并对每个训练样本子集中的参数进行归一化处理；

d）针对步骤c）中的每个所述训练样本数据库建立一个对应的径向基函数神经网络模型，并使n个径向基函数神经网络模型并联组成混合网络；

e）提取当前时刻过热器***的实际测量值，对径向基函数神经网络模型进行在线参数辨识，得到当前时刻需要辨识的参数的数值；

f）每隔时间间隔t，对径向基函数神经网络模型进行更新处理，由新的混合网络对需要辨识的参数进行在线参数辨识。

2.如权利要求1所述的燃煤电站锅炉过热器模型参数辨识方法，其特征在于，所述步骤b）中，对所述已知参数进行坏点处理是指：通过多项式滑动拟合方法来判断坏点，并采用七点二阶前推差分算式对坏点进行剔除。

3.如权利要求1所述的燃煤电站锅炉过热器模型参数辨识方法，其特征在于，所述步骤b）中，对所述已知参数进行数据平滑处理是采用七点加权滤波法实现。

4.如权利要求1所述的燃煤电站锅炉过热器模型参数辨识方法，其特征在于，所述步骤c）中，所述已知参数实时数据库进行n段划分处理是指：确定分段个数为n，在保证每段数据库个数相等或近似相等的情况下，将已知参数实时数据库按照负荷大小分为n段，得到n个不同负荷段的训练样本数据库。

5.如权利要求1所述的燃煤电站锅炉过热器模型参数辨识方法，其特征在于，所述步骤c）中，所述归一化处理是指：通过计算获得训练样本子集中的参数的归一化值。

6.如权利要求1所述的燃煤电站锅炉过热器模型参数辨识方法，其特征在于，所述步骤d）中，所述建立一个对应的径向基函数神经网络模型包括以下步骤：

第一步、确定径向基函数神经网络模型的输出量；

第二步、从所述锅炉过热器的非线性集总参数模型中查找与所述输出量有关系的已知参数，选取径向基函数神经网络模型的输入量；

第三步、采用径向基函数神经网络训练方法对所述输入量和所述输出量进行训练，得到径向基函数神经网络模型。

7.如权利要求1所述的燃煤电站锅炉过热器模型参数辨识方法，其特征在于，所述步骤e）中，所述在线参数辨识是指：对当前时刻的锅炉过热器的实际测量值依次进行坏点处理、数据平滑处理和归一化处理，并根据当前负荷的大小将处理的实际测量值输入相应负荷段的径向基函数神经网络模型中，从而得到当前时刻需要辨识的参数的数值。

8.如权利要求1所述的燃煤电站锅炉过热器模型参数辨识方法，其特征在于，所述步骤f）中的所述更新处理，包括以下步骤：

第一步、将过去t时间内的过热器***的实际测量数据依次进行坏点处理、数据平滑处理和归一化处理，得到过去t时间内的处理后的实际测量数据；

第二步、根据过去t时间内每一时刻的负荷，得到过去t时间内每一时刻处理后的实际测量数据所属的训练样本数据库，以及过去t时间内每一时刻处理后的实际测量数据所属的训练样本数据子集；

第三步、将属于同一训练样本数据子集的处理后的实际测量数据按从前向后的顺序，依次替换原有的样本数据，得到更新后的训练样本数据库；

第四步、针对更新后的训练样本数据库建立一个对应的新的径向基函数神经网络模型。