发明名称 一种基于多参数预测的火电机组协调控制方法
摘要 本发明公开了一种基于多参数预测的火电机组协调控制方法,该方法针对锅炉被控对象大滞后、响应过程缓慢等特性,通过对给水流量、给煤量、汽机阀门开度的合理估计,及主汽压力和分离器出口温度的偏差预测,分别构成主汽压力控制回路和分离器出口温度控制回路,并于负荷控制回路相协调,构成火电机组协调控制系统。本发明有效提升了火电机组协调控制系统的控制品质,在快速响应负荷变动的同时,有效抑制负荷变动时的主汽压力波动,减小主汽压力和设定值之间的偏差,抑制分离器出口温度波动,减少分离器出口温度超温情况的发生。同时,提高锅炉控制系统的稳定性,和机组的运行效率。
申请公布号 CN102608911A 申请公布日期 2012.07.25
申请号 CN201210066442.7 申请日期 2012.03.14
申请人 东南大学 发明人 陈雨亭;吕剑虹;秦文炜;崔晓波
分类号 G05B11/42(2006.01)I 主分类号 G05B11/42(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 柏尚春
主权项 1.一种基于多参数预测的火电机组协调控制方法,其特征在于该控制方法包括如下两个控制回路:主汽压力控制回路、分离器出口温度控制回路;其中主汽压力控制回路包括如下结构:保持预测器、线性增量预测器、主汽压力模型、主汽压力控制器;分离器出口温度控制回路包括如下结构:保持预测器、线性增量预测器、分离器出口温度模型、分离器出口温度控制器;线性增量预测器为按增量线性递减的方式求取输出,即当该预测器当前输入为u<sub>z</sub>(k)时,线性增量预测器输出为Y<sub>z</sub>(k)=[y<sub>z</sub>(k+1)y<sub>z</sub>(k+2)...y<sub>z</sub>(k+N)],其中<maths num="0001"><![CDATA[<math><mrow><msub><mi>y</mi><mi>z</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>1</mn><mo>-</mo><msup><mi>&beta;</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msup></mrow><mrow><mn>1</mn><mo>-</mo><mi>&beta;</mi></mrow></mfrac><mo>[</mo><msub><mi>u</mi><mi>z</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>u</mi><mi>z</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>]</mo><mo>+</mo><msub><mi>u</mi><mi>z</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></math>]]></maths>j=1,2,3...N,β为递减系数,N为预测参数时间长度,k为当前时刻,z为增量预测器的代号;保持预测器以当前该预测器输入作为其输出,即当该预测器当前输入为u<sub>b</sub>(k)时,保持预测器输出为Y<sub>b</sub>(k)=[u<sub>b</sub>(k)u<sub>b</sub>(k)...u<sub>b</sub>(k)]<sub>1×N</sub>,N为预测参数时间长度,b为保持预测器的代号;在以上结构基础上,基于多参数预测的火电机组协调控制方法包括如下步骤:1).分别以给水流量Fw、给煤量Fu、汽机阀门开度Tm为阶跃量,获取主汽压力Msp、分离器出口温度Tsp、负荷Ne的阶跃响应值,并记录在历史数据库中;2).通过最小二乘系统辨识法对历史数据库中的数据进行拟合,获取如下传递函数:以Fw为输入、Msp为输出的给水-主汽压力传递函数<img file="FDA0000143312960000012.GIF" wi="107" he="56" />以Fu为输入、Msp为输出的给煤-主汽压力传递函数<img file="FDA0000143312960000013.GIF" wi="107" he="56" />以Tm为输入、Msp为输出的汽机阀门-主汽压力传递函数<img file="FDA0000143312960000014.GIF" wi="107" he="50" />以Fw为输入、Ne为输出的给水-负荷传递函数<img file="FDA0000143312960000015.GIF" wi="90" he="56" />以Fu为输入、Ne为输出的给煤-负荷传递函数<img file="FDA0000143312960000016.GIF" wi="90" he="56" />以Tm为输入、Ne为输出的汽机阀门-负荷传递函数<img file="FDA0000143312960000017.GIF" wi="90" he="50" />以Fw为输入、Tsp为输出的给水-分离器出口温度传递函数<img file="FDA0000143312960000018.GIF" wi="95" he="60" />以Fu为输入、Tsp为输出的给煤-分离器出口温度传递函数<img file="FDA0000143312960000021.GIF" wi="95" he="60" />以Tm为输入、Tsp为输出的汽机阀门-分离器出口温度传递函数<img file="FDA0000143312960000022.GIF" wi="95" he="54" />3).对<img file="FDA0000143312960000023.GIF" wi="85" he="56" />进行双线性离散化得离散化传递函数<img file="FDA0000143312960000024.GIF" wi="217" he="139" />其中<img file="FDA0000143312960000025.GIF" wi="185" he="62" />为给水-主汽压力离散分子多项式,<img file="FDA0000143312960000026.GIF" wi="184" he="62" />为给水-主汽压力离散分母多项式,q为超前因子;4).构造丢番图方程:<maths num="0002"><![CDATA[<math><mrow><mn>1</mn><mo>=</mo><msubsup><mi>E</mi><mi>N</mi><mi>msp</mi></msubsup><mrow><mo>(</mo><msup><mi>q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><msubsup><mi>A</mi><mi>fw</mi><mi>msp</mi></msubsup><mrow><mo>(</mo><msup><mi>q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mi>&Delta;</mi><mo>+</mo><msup><mi>q</mi><mrow><mo>-</mo><mi>N</mi></mrow></msup><msubsup><mi>F</mi><mi>N</mi><mi>msp</mi></msubsup><mrow><mo>(</mo><msup><mi>q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>得给水-主汽压力分离第一多项式<img file="FDA0000143312960000028.GIF" wi="186" he="56" />与给水-主汽压力分离第二多项式<img file="FDA0000143312960000029.GIF" wi="212" he="56" />5).令<maths num="0003"><![CDATA[<math><mrow><msup><mi>G</mi><mi>msp</mi></msup><mrow><mo>(</mo><msup><mi>q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mo>=</mo><msubsup><mi>E</mi><mi>N</mi><mi>msp</mi></msubsup><mrow><mo>(</mo><msup><mi>q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><msubsup><mi>B</mi><mi>fw</mi><mi>msp</mi></msubsup><mrow><mo>(</mo><msup><mi>q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>+</mo><msub><mi>N</mi><mi>bfw</mi></msub></mrow></munderover><msubsup><mi>g</mi><mi>i</mi><mi>msp</mi></msubsup><msup><mi>q</mi><mrow><mo>-</mo><mi>i</mi></mrow></msup><mo>,</mo></mrow></math>]]></maths>其中G<sup>msp</sup>(q<sup>-1</sup>)为主汽压力前向多项式,<img file="FDA00001433129600000211.GIF" wi="80" he="50" />为主汽压力前向多项式系数,N<sub>bfw</sub>为给水-主汽压力离散分子多项式<img file="FDA00001433129600000212.GIF" wi="185" he="62" />的阶次;6).构造主汽压力控制系数<maths num="0004"><![CDATA[<math><mrow><msup><mi>&alpha;</mi><mi>msp</mi></msup><mo>=</mo><msup><mrow><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>i</mi><mo>=</mo><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mrow><mo>(</mo><msubsup><mi>g</mi><mi>i</mi><mi>msp</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mi>&gamma;</mi><mi>msp</mi></msup><mo>]</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>,</mo></mrow></math>]]></maths>其中γ<sup>msp</sup>为主汽压力控制量影响分量系数;7).令当前为k时刻,当前给煤量Fu(k)、汽机阀门开度Tm(k)经过线性增量预测器、给水流量Fw(k)经过保持预测器计算后,未来给水流量、给煤量、汽机阀门开度预测值为:Fw(k+j)=Fw(k-1)j=1,2,3...N<maths num="0005"><![CDATA[<math><mrow><mi>Fu</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>1</mn><mo>-</mo><msup><msub><mi>&beta;</mi><mi>fu</mi></msub><mrow><mi>j</mi><mo>+</mo><mn>2</mn></mrow></msup></mrow><mrow><mn>1</mn><mo>-</mo><msub><mi>&beta;</mi><mi>fu</mi></msub></mrow></mfrac><mi>&Delta;Fu</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mi>Fu</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>j=1,2,3...N<img file="FDA00001433129600000215.GIF" wi="752" he="128" />j=1,2,3...N,其中β<sub>fu</sub>和β<sub>tm</sub>分别为给煤控制量衰减系数与汽机阀门开度衰减系数且均为小于1大于0的常数;8).依据未来给水流量、给煤量、汽机阀门开度预测值和对应传递函数<img file="FDA00001433129600000216.GIF" wi="112" he="56" /><img file="FDA00001433129600000217.GIF" wi="85" he="56" />和<img file="FDA00001433129600000218.GIF" wi="108" he="50" />经过主汽压力模型计算得未来主汽压力预测值Msp<sub>t</sub>(k+j)    j=1,2,3...N;9).取主汽压力设定值为Msp<sub>r</sub>经过保持预测器计算后与未来主汽压力预测值作差,并由主汽压力控制器计算得主汽压力控制回路给水控制量为<maths num="0006"><![CDATA[<math><mrow><msub><mi>Fw</mi><mi>c</mi></msub><mo>=</mo><msup><mi>&alpha;</mi><mi>msp</mi></msup><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mo>[</mo><msub><mi>Msp</mi><mi>r</mi></msub><mo>-</mo><msub><mi>Msp</mi><mi>t</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>i</mi><mo>)</mo></mrow><mo>]</mo><msubsup><mi>g</mi><mi>i</mi><mi>msp</mi></msubsup><mo>;</mo></mrow></math>]]></maths>10).将Fw<sub>c</sub>作为给水流量的目标值,传送到给水泵执行机构中,通过改变给水流量控制主汽压力;a)给水流量增大,给煤量不变,锅炉过热段减少,主汽压力减小;b)给水流量减少,给煤量不变,锅炉过热段增长,主汽压力增大;从而控制锅炉主汽压力;11).对<img file="FDA0000143312960000032.GIF" wi="72" he="60" />进行双线性离散化得离散化传递函数<img file="FDA0000143312960000033.GIF" wi="205" he="139" />其中<img file="FDA0000143312960000034.GIF" wi="171" he="62" />为给煤-分离器出口温度离散分子多项式,<img file="FDA0000143312960000035.GIF" wi="169" he="62" />为给煤-分离器温度离散分母多项式;12).构造丢番图方程,<maths num="0007"><![CDATA[<math><mrow><mn>1</mn><mo>=</mo><msubsup><mi>E</mi><mi>N</mi><mi>tsp</mi></msubsup><mrow><mo>(</mo><msup><mi>q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><msubsup><mi>A</mi><mi>fu</mi><mi>tsp</mi></msubsup><mrow><mo>(</mo><msup><mi>q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mi>&Delta;</mi><mo>+</mo><msup><mi>q</mi><mrow><mo>-</mo><mi>N</mi></mrow></msup><msubsup><mi>F</mi><mi>N</mi><mi>tsp</mi></msubsup><mrow><mo>(</mo><msup><mi>q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>得给煤-分离器出口温度分离第一多项式<img file="FDA0000143312960000037.GIF" wi="173" he="56" />与给煤-分离器出口温度分离第二多项式<maths num="0008"><![CDATA[<math><mrow><msubsup><mi>F</mi><mi>N</mi><mi>tsp</mi></msubsup><mrow><mo>(</mo><msup><mi>q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>13).令<maths num="0009"><![CDATA[<math><mrow><msup><mi>G</mi><mi>tsp</mi></msup><mrow><mo>(</mo><msup><mi>q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mo>=</mo><msubsup><mi>E</mi><mi>N</mi><mi>tsp</mi></msubsup><mrow><mo>(</mo><msup><mi>q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><msubsup><mi>B</mi><mi>fu</mi><mi>tsp</mi></msubsup><mrow><mo>(</mo><msup><mi>q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>+</mo><msub><mi>N</mi><mi>bfu</mi></msub></mrow></munderover><msubsup><mi>g</mi><mi>i</mi><mi>tsp</mi></msubsup><msup><mi>q</mi><mrow><mo>-</mo><mi>i</mi></mrow></msup><mo>,</mo></mrow></math>]]></maths>其中G<sup>tsp</sup>(q<sup>-1</sup>)为分离器出口温度前向多项式,<img file="FDA00001433129600000310.GIF" wi="66" he="54" />为分离器出口温度前向多项式系数,N<sub>bfu</sub>为给煤-分离器出口温度离散分子多项式<img file="FDA00001433129600000311.GIF" wi="171" he="62" />的阶次;14).构造分离器出口温度控制系数<img file="FDA00001433129600000312.GIF" wi="554" he="137" />其中γ<sup>tsp</sup>为分离器出口温度控制量影响分量系数;15).令当前为k时刻,当前给水流量为Fw(k)、汽机阀门开度为Tm(k)经过线性增量预测器、给煤量Fu(k)经过保持预测器计算后,得未来给水流量、给煤量、汽机阀门开度预测值为Fu(k+j)=Fu(k-1)    j=1,2,3...N<maths num="0010"><![CDATA[<math><mrow><mi>Fw</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>1</mn><mo>-</mo><msup><msub><mi>&beta;</mi><mi>fw</mi></msub><mrow><mi>j</mi><mo>+</mo><mn>2</mn></mrow></msup></mrow><mrow><mn>1</mn><mo>-</mo><msub><mi>&beta;</mi><mi>fw</mi></msub></mrow></mfrac><mi>&Delta;Fw</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mi>Fw</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>j=1,2,3...N<img file="FDA0000143312960000042.GIF" wi="752" he="128" />j=1,2,3...N,其中β<sub>fw</sub>为给水控制量衰减系数,且为小于1大于0的常数。16).依据未来给水流量、给煤量、汽机阀门开度预测值和对应传递函数<img file="FDA0000143312960000043.GIF" wi="100" he="60" /><img file="FDA0000143312960000044.GIF" wi="72" he="60" />和<img file="FDA0000143312960000045.GIF" wi="72" he="54" />经过分离器出口温度模型计算,得未来分离器出口温度预测值Tsp<sub>t</sub>(k+j)   j=1,2,3...N;17).取分离器出口温度设定值为Tsp<sub>r</sub>,经过保持预测器计算后与未来分离器出口温度预测值作差,并由分离器出口温度控制器计算得给煤控制量为<maths num="0011"><![CDATA[<math><mrow><msub><mi>Fu</mi><mi>c</mi></msub><mo>=</mo><msup><mi>&alpha;</mi><mi>tsp</mi></msup><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mo>[</mo><msub><mi>Tsp</mi><mi>r</mi></msub><mo>-</mo><msub><mi>Tsp</mi><mi>t</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>i</mi><mo>)</mo></mrow><mo>]</mo><msubsup><mi>g</mi><mi>i</mi><mi>tsp</mi></msubsup><mo>;</mo></mrow></math>]]></maths>18).将Fu<sub>c</sub>作为给煤量的目标值,传送到给煤机执行机构中,通过改给煤量控制分离器出口温度;a)给煤机转速上升,给煤量增加,锅炉过热段增长,过热蒸汽出口温度上升,分离器出口温度上升;b)给煤机转速下降,给煤量减少,锅炉过热段减少,过热蒸汽出口温度下降,分离器出口温度下降;从而控制锅炉分离器出口温度。
地址 210096 江苏省南京市四牌楼2号