发明名称 一种基于母管制的多机组汽轮机蒸汽流量实时分配方法
摘要 一种基于母管制的多机组汽轮机蒸汽流量实时分配方法,利用数据挖掘技术处理母管制机组在各台汽轮机单独运行时的历史稳态数据,建立各台汽轮机的功率和进汽热流量函数关系计算模型,在机组多台汽轮机同时运行时,将实时采集到的各台汽轮机的负荷、母管蒸汽压力、蒸汽温度数据输入计算模型,计算模型会输出对应当前控制边界条件下各台汽轮机进汽热流量,以计算出各台汽轮机需要分配的蒸汽流量;同时以定期更新计算模型的方式考虑各台汽轮机随运行时间变化其特性变化的影响。
申请公布号 CN105201565A 申请公布日期 2015.12.30
申请号 CN201510640476.6 申请日期 2015.09.30
申请人 西安西热电站信息技术有限公司 发明人 王智微;李虹锐
分类号 F01D17/00(2006.01)I 主分类号 F01D17/00(2006.01)I
代理机构 西安智大知识产权代理事务所 61215 代理人 刘国智
主权项 一种基于母管制的多机组汽轮机蒸汽流量实时分配方法,其特征在于,包括以下步骤:1)收集母管制机组单台汽轮机运行工况的历史稳态数据设电厂总计n台汽轮机,在电厂的信息系统里查询第j台汽轮机单独运行的工况,其中1≤j≤n,导出第j台汽轮机的历史稳态运行数据,提取历史稳态运行数据中的负荷<img file="FDA0000814863070000011.GIF" wi="76" he="82" />和该负荷下对应的母管蒸汽压力<img file="FDA0000814863070000012.GIF" wi="136" he="88" />蒸汽温度<img file="FDA0000814863070000013.GIF" wi="133" he="85" />每台锅炉出口蒸汽流量<img file="FDA0000814863070000014.GIF" wi="115" he="84" />数据;2)建立各台汽轮机功率和进汽热流量函数关系的计算模型,运用分类、回归数据挖掘技术,建立有监督的机器自学习模型,作为该台汽轮机功率和进汽热流量函数关系计算模型:2‑1)将提取出来的负荷<img file="FDA0000814863070000015.GIF" wi="70" he="80" />和该负荷下对应的母管蒸汽压力<img file="FDA0000814863070000016.GIF" wi="131" he="84" />蒸汽温度<img file="FDA0000814863070000017.GIF" wi="128" he="84" />每台锅炉出口蒸汽流量<img file="FDA0000814863070000018.GIF" wi="114" he="85" />数据随机分开,70%数据用作计算模型的自学习,30%数据作为检验,自学习数据记为<img file="FDA0000814863070000019.GIF" wi="644" he="118" />检验数据记为<img file="FDA00008148630700000110.GIF" wi="672" he="124" />2‑2)对数据进一步处理,已知压力P、温度T和水蒸气热物性表计算出蒸汽焓值<img file="FDA00008148630700000111.GIF" wi="248" he="102" />和第j台汽轮机在负荷<img file="FDA00008148630700000112.GIF" wi="70" he="84" />下的进汽热流量<img file="FDA00008148630700000113.GIF" wi="302" he="94" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>Q</mi><mrow><mo>(</mo><mi>p</mi><mi>e</mi><mo>)</mo></mrow><mi>j</mi></msubsup><mo>=</mo><msubsup><mi>h</mi><mrow><mo>(</mo><mi>P</mi><mi>e</mi><mo>)</mo></mrow><mi>j</mi></msubsup><mo>&CenterDot;</mo><msubsup><mi>&Sigma;D</mi><mrow><mi>b</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>(</mo><mi>P</mi><mi>e</mi><mo>)</mo></mrow></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00008148630700000114.GIF" wi="589" he="107" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>Q</mi><mrow><mo>(</mo><mi>p</mi><mi>e</mi><mo>)</mo></mrow><msup><mi>j</mi><mo>&prime;</mo></msup></msubsup><mo>=</mo><msubsup><mi>h</mi><mrow><mo>(</mo><mi>P</mi><mi>e</mi><mo>)</mo></mrow><msup><mi>j</mi><mo>&prime;</mo></msup></msubsup><mo>&CenterDot;</mo><msubsup><mi>&Sigma;D</mi><mrow><mi>b</mi><mo>,</mo><mi>i</mi></mrow><msup><mrow><mo>(</mo><mi>P</mi><mi>e</mi><mo>)</mo></mrow><mo>&prime;</mo></msup></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00008148630700000115.GIF" wi="591" he="106" /></maths>2‑3)建立结构为1‑5‑1的BP神经网络,即输入层有1个节点,隐含层有5个节点,输入层有1个节点,节点传递函数采用对数S型传递函数logsig,训练函数采用Levenberg_Marquardt的BP算法训练函数trainlm,网络学习函数采用带动量项的BP学习规则learngdm,性能分析函数采用均方差性能分析函数mse,将所有自学习数据<img file="FDA0000814863070000021.GIF" wi="250" he="91" />使用最大最小法进行数据归一化,最大最小法函数形式如下:x'<sub>k</sub>=(x<sub>k</sub>‑x<sub>min</sub>)/(x<sub>max</sub>‑x<sub>min</sub>)  (3)式中:x'<sub>k</sub>—归一化后的值;x<sub>k</sub>—归一化前的原始值;x<sub>min</sub>—原始值中的最小值;x<sub>max</sub>—原始值中的最大值;将归一化后的自学习数据输入BP神经网络模型进行训练,设置训练终止条件为拟合结果误差小于等于0.0001,验证检查结果大于等于10,训练结束得到第j台汽轮机负荷<img file="FDA0000814863070000022.GIF" wi="73" he="81" />和进汽热流量<img file="FDA0000814863070000023.GIF" wi="114" he="91" />的输入输出函数关系,即第j台汽轮机负荷和热耗的计算模型:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>Q</mi><mrow><mo>(</mo><mi>p</mi><mi>e</mi><mo>)</mo></mrow><mi>j</mi></msubsup><mo>=</mo><msup><mi>f</mi><mi>j</mi></msup><mrow><mo>(</mo><msub><mi>p</mi><mi>e</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000814863070000024.GIF" wi="456" he="91" /></maths>2‑4)将检验数据<img file="FDA0000814863070000025.GIF" wi="260" he="93" />使用自学习数据的归一化参数即自学习数据中的最大,最小值进行归一化处理,输入到第j台汽轮机负荷和进汽热流量的计算模型,得到第j台汽轮机在负荷为<img file="FDA0000814863070000026.GIF" wi="89" he="86" />时进汽热流量的模型计算值<img file="FDA0000814863070000027.GIF" wi="180" he="88" /><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>Q</mi><mrow><mo>(</mo><mi>p</mi><mi>e</mi><mo>)</mo></mrow><mi>j</mi></msubsup><mo>,</mo><mi>c</mi><mo>=</mo><msup><mi>f</mi><mi>j</mi></msup><mrow><mo>(</mo><msubsup><mi>p</mi><mi>e</mi><msup><mi>j</mi><mo>&prime;</mo></msup></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000814863070000028.GIF" wi="501" he="94" /></maths>2‑5)记事件K为:<img file="FDA0000814863070000029.GIF" wi="465" he="205" />若P(K)≥95%,则称计算模型收敛正确;如果P(K)<95%,则称计算模型收敛错误,调整BP神经网络的隐含层节点数,或者增加隐含层层数,以提高BP神将网络的非线性函数拟合能力,重新执行步骤2)即重新学习、检验,直至计算模型正确收敛;3)多台汽轮机运行时实时计算应分配给的计算蒸汽流量机组多台汽轮机运行时,将实时采集到的各台汽轮机的负荷<img file="FDA0000814863070000031.GIF" wi="74" he="80" />数据归一化后输入计算模型,计算模型输出对应当前控制边界条件下各台汽轮机的进汽热流量计算值<img file="FDA0000814863070000032.GIF" wi="101" he="81" />已知母管蒸汽压力P、蒸汽温度T,计算出当前控制边界条件下各台汽轮机的计算蒸汽流量D<sub>t,j,c</sub>;4)根据计算蒸汽流量D<sub>t,j,c</sub>分配各台汽轮机蒸汽流量通过系统实时采集到的各台锅炉进入母管的蒸汽流量<img file="FDA0000814863070000033.GIF" wi="118" he="107" />计算母管蒸汽的总流量<img file="FDA0000814863070000034.GIF" wi="215" he="106" />然后利用模型输出的第j台汽轮机的计算蒸汽流量D<sub>t,j,c</sub>占所有运行汽轮机计算蒸汽流量之和∑D<sub>t,j,c</sub>的比例来得到第j台汽轮机实际分配到的蒸汽流量D<sub>t,j</sub>:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>D</mi><mrow><mi>t</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mfrac><msub><mi>D</mi><mrow><mi>t</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>c</mi></mrow></msub><mrow><msub><mi>&Sigma;D</mi><mrow><mi>t</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>c</mi></mrow></msub></mrow></mfrac><mo>&CenterDot;</mo><msubsup><mi>&Sigma;D</mi><mrow><mi>b</mi><mo>,</mo><mi>i</mi></mrow><msub><mi>p</mi><mi>e</mi></msub></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA0000814863070000035.GIF" wi="658" he="161" /></maths>
地址 710032 陕西省西安市兴庆路136号