发明名称 一种电站锅炉烟气软测量方法
摘要 本发明提供一种基于最小二乘支持向量机及在线更新的电站锅炉烟气软测量方法,属于热工技术和人工智能交叉技术领域。该方法选择电站锅炉有关运行和状态参数作为模型的输入,要预测的烟气成分含量作为模型的输出,选取历史运行数据作为初始训练样本,利用最小二乘支持向量机方法建立烟气排放的初始模型。另外,基于对烟气排放时变特性的分析,提出了基于样本替换和样本追加的更新策略,并采用删减样本和增加样本两种模式以增量的形式来实现参数的求解和模型的更新。本发明提出的最小二乘支持向量机及在线更新软测量方法随着过程特性的变化自适应地改进模型性能,能够实现对烟气排放的精确预测,对电站锅炉的安全和优化运行有重要的意义。
申请公布号 CN103728879B 申请公布日期 2016.08.17
申请号 CN201410024169.0 申请日期 2014.01.20
申请人 华北电力大学 发明人 吕游;杨婷婷;刘吉臻
分类号 G05B13/00(2006.01)I 主分类号 G05B13/00(2006.01)I
代理机构 北京金智普华知识产权代理有限公司 11401 代理人 巴晓艳
主权项 一种基于LSSVM及在线更新的电站锅炉烟气软测量方法,其特征在于,所述方法采用LSSVM构建初始烟气排放模型,然后利用样本追加和样本替换来实现模型的增量更新,能够对锅炉烟气各成分进行精确地预测,所述方法包括以下步骤:步骤1:收集初始训练样本来构建LSSVM模型;步骤2:利用此模型对烟气成分含量进行预测;步骤3:当实际的烟气成分含量的传感器测量值y<sub>q</sub>采集到后,计算样本(x<sub>q</sub>,y<sub>q</sub>)的预测误差Er;步骤4:判断预测误差:若Er&gt;Δ,Δ为误差阈值,则执行步骤5,否则执行步骤9;步骤5:从历史运行数据中选取距新采样样本(x<sub>q</sub>,y<sub>q</sub>)最近的样本点(x<sub>k</sub>,y<sub>k</sub>),其中<img file="FDA0000992235660000011.GIF" wi="527" he="107" />步骤6:根据新采样样本(x<sub>q</sub>,y<sub>q</sub>)与其最近的样本(x<sub>k</sub>,y<sub>k</sub>)之间的距离进行判断并确定更新类型;步骤7:根据确定的更新类型,对步骤1获得的初始LSSVM模型进行增量更新,根据步骤6进行样本增加或样本替换来更新;步骤8:所述步骤1中LSSVM模型可描述为以下优化问题,<maths num="0001"><math><![CDATA[<mrow><munder><mi>min</mi><mrow><mi>w</mi><mo>,</mo><mi>b</mi><mo>,</mo><mi>&xi;</mi></mrow></munder><mi>J</mi><mrow><mo>(</mo><mrow><mi>w</mi><mo>,</mo><mi>&xi;</mi></mrow><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>w</mi><mi>T</mi></msup><mi>w</mi><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>&gamma;</mi><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><msub><mi>&xi;</mi><mi>i</mi></msub><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000992235660000012.GIF" wi="1825" he="162" /></maths><img file="FDA0000992235660000013.GIF" wi="1201" he="113" />其中,J为目标函数,<img file="FDA0000992235660000015.GIF" wi="92" he="55" />是核空间映射函数,w为权重向量,γ为惩罚系数,ξ<sub>i</sub>为误差变量,b为偏差,上标T表示矩阵的转置;为解此优化问题,定义Lagrange函数如下:<img file="FDA0000992235660000014.GIF" wi="1847" he="154" />其中,α=[α<sub>1</sub>,…,α<sub>n</sub>]<sup>T</sup>为Lagrange乘子;利用Lagrange函数对各变量求偏导,并令导数值为零可得到:<img file="FDA0000992235660000021.GIF" wi="948" he="617" />消去中间变量w和ξ<sub>i</sub>,将其转化为求解线性方程组:<maths num="0002"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>0</mn></mtd><mtd><msup><mover><mi>I</mi><mo>&RightArrow;</mo></mover><mi>T</mi></msup></mtd></mtr><mtr><mtd><mover><mi>I</mi><mo>&RightArrow;</mo></mover></mtd><mtd><mrow><mi>&Omega;</mi><mo>+</mo><mn>1</mn><mo>/</mo><mi>&gamma;</mi><mi>I</mi></mrow></mtd></mtr></mtable></mfenced><mfenced open = "[" close = "]"><mtable><mtr><mtd><mi>b</mi></mtd></mtr><mtr><mtd><mi>a</mi></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mi>y</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000992235660000022.GIF" wi="1327" he="166" /></maths>其中y=[y<sub>1</sub>,…,y<sub>n</sub>]<sup>T</sup>,<img file="FDA0000992235660000025.GIF" wi="283" he="73" />I为n×n阶单位矩阵,Ω={Ω<sub>kl</sub>|k,l=1,…,n},且<img file="FDA0000992235660000023.GIF" wi="774" he="98" />定义为核函数;通过求解方程组得到α和b的值为:<maths num="0003"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>b</mi><mo>=</mo><mfrac><mrow><msup><mover><mi>I</mi><mo>&RightArrow;</mo></mover><mi>T</mi></msup><msup><mi>H</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>y</mi></mrow><mrow><mover><mi>I</mi><mo>&RightArrow;</mo></mover><msup><mi>H</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mover><mi>I</mi><mo>&RightArrow;</mo></mover></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mi>a</mi><mo>=</mo><msup><mi>H</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mrow><mi>y</mi><mo>-</mo><mover><mi>I</mi><mo>&RightArrow;</mo></mover><mi>b</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000992235660000024.GIF" wi="1927" he="322" /></maths>用求得的新的模型的特征矩阵H<sup>‑1</sup>计算得到相应的模型参数α和b,实现对烟气软测量模型h(x)的更新;步骤9:判断测试样本是否结束,若结束则停止程序;否则执行步骤2,对下一样本进行预测。
地址 102206 北京市昌平区回龙观镇华北电力大学