专利名称：一种锅炉在线燃烧优化的模型更新方法
技术领域：
本发明属于信息控制技术领域，涉及到增量学习技术，特别是涉及一 种锅炉在线燃烧优化的模型更新方法。
背景技术：
电站锅炉的燃烧优化是节能减排的重要技术手段，其目标是通过调整 锅炉配风、给煤等运行参数而获得高效率或低污染排放的运行状态。锅炉 的很多运行参数如配风情况、给煤情况及氧量等的搭配对锅炉燃烧状态都 有直接的影响，不同的配风、给煤等操作参数的配置会直接导致不同的锅 炉效率及污染气体的排放的情况。对于给定的电站锅炉，在一定的负荷条 件下，对于不同的燃烧状态的指标存在一种最优的操作参数配置方案，能 够使燃烧状态最优化。然而，锅炉的操作参数间有着非常复杂的耦合关系， 要找到最优的操作参数的配置是非常困难。随着科学技术的不断进步，锅 炉运行自动化程度不断提高，但是锅炉燃烧优化问题一直没有很好的得到 解决。
目前锅炉燃烧优化的研究热点是通过数据挖掘，在大量不同的实际运 行参数组合中，应用机器学习的方法，挖掘出运行参数与锅炉燃烧指标间 的关系模型，再利用优化算法结合挖掘出的模型进行锅炉的在线燃烧优化。 这种方法节省人力物力，而且可以找到比人工实验更优的参数配置，但是 由于锅炉设备的特性随着使用时间的增长会有所改变，而且燃用煤种有时 也会变化，因此如何保证模型能够快速、高效的更新以适应新的情况成为 了这种方法的瓶颈问题。该问题与建模方法、样本数据选取及更新策略等 都有很大关系。
发明内容
本发明的目的是针对锅炉燃烧优化中的瓶颈问题，提出一种兼顾模型 历史学习结果与新的变化情况的模型更新方法。
本发明具体是利用超出己存在模型预测误差限度的数据建立新的模 型，然后再利用锅炉燃烧数据和优化算法，寻找最优的新模型和已有模型 的权重系数，利用最优的权重系数将新模型和已有模型结合在一起，共同 对新的锅炉工况进行预测，实现模型更新。该方法克服了传统更新方法中 将已有模型完全放弃，不能利用己有模型的学习结果的缺点，充分利用了 己有模型的学习结果，大大縮短了模型更新的计算工作量和时间。
本发明的技术方案是通过利用已有模型预测超出误差限度的数据作样 本，建立新的模型，并将新模型与己有模型相结合等手段，确立了一种锅 炉在线燃烧优化的模型更新方法，利用该方法可以快速、高效的实现模型本发明方法的步骤包括
步骤(1)建立已存在模型的预测错误数据库。根据具体锅炉燃烧情况和 对模型预测精度的要求，设定模型的允许预测误差限s，在采集数据时， 判断模型预测值与实际运行值之间的误差与允许预测误差限S的大小，如
果预测误差大于s，即i^-kI〉^其中^为模型预测值，r为实际运行数
据，则将超限数据存入预测错误数据库中，以备模型更新之用。
步骤(2)建立新模型。原模型需要更新时，利用预测错误数据库中的数
据作为训练样本，样本可以表示为^,J^:，其中、表示第z组作为输入数
据的锅炉运行参数向量，X表示第Z组作为输出参数的表征锅炉燃烧状态的 参数，采用支持向量机算法建立新的模型，核函数选为径向基函数
〖"，、)=00,) 00,) = exp
义,-_X;
-(x)为映射函数，设所求的目标函数为/(x,) = w—(义,)+ 6， /(义,)为模
型输出的锅炉燃烧指标预测值，w为权重系数向量，6为截距。引入松弛因 子^》0和DO和允许拟合误差e，模型可以通过在约束 _y,. — w. — 6 S ￡■ +《 〔x,) + Zj-乂 " +《
《20 "l,…，iV，条件下，最小化
《2 0
min 《,《*)二 * w. w+cZ《+《*
获得，其中常数O0为惩罚系数。该最小化问题为一个凸二次规划问题，
引入拉格朗日函数
丄(一， , a V,卢)=垂w , +《(纟+ 。 —1>, [乂.—(《+ e + ))]
_z [《+s+) 1 ] - z ("+")
其中a,, ;》o， z,，<》o，为拉格朗日乘数。
在鞍点处，函数L是关于w，b， ^， H勺极小点，也是",，《，y,乂极大点.
最小化问题转化为求其对偶问题的最大化问题。
拉格朗日函数L在鞍点处是关于W， b， L，〖，，炎小点，得
5丄=0 —》《,,<) = 0
会丄=0 —C_ X =0
可得拉格朗日函数的对偶函数:
07
此时，
/(义)=Z o, - )《"a )+6
按照库恩-塔克(KKT)条件定理，在鞍点有下式成: [",["《1 + /(义,)]=0
,*[s +《+ x-/(J,)] = 0
由上式可见，"/《",和"/都不会同时为非零，可得:
Kj, =0
从上式可求出b，获得模型。清空预测错误数据库，以备下一次更新 采集数据之用。
步骤(3)确定新模型和已有模型的权重。采集最新锅炉不同运行工况下
的数据作为检验样本，应用原模型预测与新模型预测加权平均的方法，对 检验数据进行预测，即}',= 化+/ 化，其中y,.为第，组工况的优化目标预测 值，化为新模型预测值，D。为原有模型预测值，"为新模型预测值权重系
数，y9为原有模型的预测权重系数，且《 + / = ]， 《与/ 的确定采用粒子群 算法迭代寻优确定，定义粒子群算法初始位置Z向量的各维分量，分别为新模型权重"和原有模型权重/ ，目标函数为mini/ = Z|X-i;|，其中《为 第/组工况的实际运行值。当J达到了最小、达到设定值或完成迭代次数时， 训练完成，获得新模型和己有模型的权重系数。
步骤(4)将原有模型与新模型按最优的权重比例相结合，构成新的模 型，即/) = 1) +^)。，其中D为更新后的模型，实现锅炉燃烧优化模型的更新。
本发明提出的模型更新方法充分利用了已有的模型所包含的有用信 息，大大减少了模型更新的工作量，提高了模型更新的效率，满足了锅炉 燃烧在线优化的实际要求，保证了锅炉燃烧优化的实时性和准确性。
具体实施例方式
一种锅炉在线燃烧优化的模型更新方法，具体步骤是
(1) 建立已存在模型的预测错误数据库。根据具体锅炉燃烧情况和对模 型预测精度的要求，设定模型的允许预测误差限S，在采集数据时，判断 模型预测值与实际运行值之间的误差与允许预测误差限S的大小，如果预 测误差大于S，即|^-">^其中K为模型预测值，K为实际运行数据，
则将超限数据存入预测错误数据库中，以备模型更新之用。
(2) 建立新模型。原模型需要更新时，利用预测错误数据库中的数据作
为训练样本，样本可以表示为^,J^:，其中:c,表示第/组作为输入数据的 锅炉运行参数向量，y,.表示第/组作S输出参数的表征锅炉燃烧状态的参 数，采用支持向量机算法建立新的模型，核函数选为径向基函数
〖(x,，x,)=《x,)—(> ) = exp (11 '/11 )
0(x)为映射函数，设所求的目标函数为/(x,) = >v—".) + Z>， /(x,)为模 型输出的锅炉燃烧指标预测值，w为权重系数向量，6为截距。引入松弛因 子1、》0和DO和允许拟合误差e，模型可以通过在约束
min 《,《*) = * w. >v+cZ《+《*
获得，其中常数DO为惩罚系数。该最小化问题为一个凸二次规划问题，
引入拉格朗曰函数
WXI,) + 6-乂 Sf+《■
条件下，最小化:
7丄ov》，《，，,《， 《*,=* w ,+(《+，)—|>,[乂's+/" ))]
-Z [《+" /(a ) i ] - Z ("'.+ )
其中《,, ;>o， r,，r;^0，为拉格朗日乘数。
在鞍点处，函数L是关于w， b，〖,，〖，*的极小点，也是a,，^,r,乂极大点:
最小化问题转化为求其对偶问题的最大化问题。
拉格朗日函数L在鞍点处是关于W，b， " H及小点，得:￡丄=0 — w = g" _ O(x,.)
可得拉格朗日函数的对偶函数
-z ( ,+)乂.
此时，
/(力=z ("' - a)+6
按照库恩-塔克(KKT)条件定理，在鞍点有下式成立
+《-y,+/(x,)] = 01 ,["《+乂'_/(1,)] = 0
由上式可见，"，和"/都不会同时为非零，可得f《7, = 0 1" =0
从上式可求出b，获得模型。清空预测错误数据库，以备下一次更新 采集数据之用。
(3) 确定新模型和已有模型的权重。采集最新锅炉不同运行工况下的数
据作为检验样本，应用原模型预测与新模型预测加权平均的方法，对检验
数据进行预测，即乃= 马+"化，其中少,为第/组工况的优化目标预测值，A, 为新模型预测值，D。为原有模型预测值，a为新模型预测值权重系数，p为 原有模型的预测权重系数，且《 + ^ = 1，"与/ 的确定采用粒子群算法迭代 寻优确定，定义粒子群算法初始位置z向量的各维分量，分别为新模型权
重a和原有模型权重^，目标函数为mini/ = Z|X-"，其中"为第/组工况
的实际运行值。当J达到了最小、达到设定值或完成迭代次数时，训练完 成，获得新模型和己有模型的权重系数。
(4) 将原有模型与新模型按最优的权重比例相结合，构成新的模型，即 Z^"A+^D。，其中D为更新后的模型，实现锅炉燃烧优化模型的更新。
权利要求
1、一种锅炉在线燃烧优化的模型更新方法，其特征在于该方法的步骤包括步骤(1). 建立已存在模型的预测错误数据库根据具体锅炉燃烧情况和对模型预测精度的要求，设定模型的允许预测误差限δ，在采集数据时，判断模型预测值与实际运行值之间的误差与允许预测误差限δ的大小，如果预测误差大于δ，即|Vc-Vs|>δ，则将超限数据存入预测错误数据库中，其中Vc为模型预测值，Vs为实际运行数据；步骤(2). 建立新模型原模型需要更新时，利用预测错误数据库中的数据作为训练样本，样本表示为其中xi表示第i组作为输入数据的锅炉运行参数向量，yi表示第i组作为输出参数的表征锅炉燃烧状态的参数，采用支持向量机算法建立新的模型，核函数选为径向基函数φ(x)为映射函数，设所求的目标函数为f(xi)＝w·φ(xi)+b，f(xi)为模型输出的锅炉燃烧指标预测值，w为权重系数向量，b为截距；w、b由在条件求解而得；式中引入松弛因子ξ*i和ξi，ξ*i≥0和ξi≥0，ε为允许拟合误差；步骤(3). 确定新模型和已有模型的权重；采集锅炉不同运行工况下的数据作为检验样本，应用原模型预测与新模型预测加权平均的方法，对检验数据进行预测，即yi＝αDn+βDo，其中yi为第i组工况的优化目标预测值，Dn为新模型预测值，Do为原有模型预测值，α为新模型预测值权重系数，β为原有模型的预测权重系数，且α+β＝1；α与β的确定采用粒子群算法迭代寻优确定，定义粒子群算法初始位置Z向量的各维分量，分别为新模型权重α和原有模型权重β，目标函数为minJ＝∑|yi-Yi|，其中Yi为第i组工况的实际运行值，当J达到了最小、达到设定值或完成迭代次数时，训练完成，获得新模型和已有模型的权重系数；步骤(4). 将原有模型与新模型按最优的权重比例相结合，构成新的模型，即D＝αDn+βDo，其中D为更新后的模型，实现锅炉燃烧优化模型的更新。
全文摘要
本发明涉及一种锅炉在线燃烧优化的模型更新方法。本发明方法利用超出已存在模型预测误差限度的数据建立新的模型，然后再利用锅炉燃烧数据和优化算法，寻找最优的新模型和已有模型的权重系数，利用最优的权重系数将新模型和已有模型结合在一起，共同对新的锅炉工况进行预测，实现模型更新。本发明方法克服了传统更新方法中将已有模型完全放弃，不能利用已有模型的学习结果的缺点，充分利用了已有模型的学习结果，大大缩短了模型更新的计算工作量和时间。
文档编号G05B19/418GK101498458SQ200910096410
公开日2009年8月5日 申请日期2009年3月2日 优先权日2009年3月2日
发明者张日东, 王建中, 王春林, 铭 葛, 薛安克 申请人:杭州电子科技大学