发明名称 一种综合节能优化调度日计划生成方法
摘要 本发明涉及一种综合节能优化调度日计划生成方法,包括通过煤耗在线监测系统获得机组供电煤耗信息,通过烟气排放在线监测系统获得脱硫效率信息,通过热电机组在线监测系统获得供热负荷信息,通过资源综合利用在线监测系统获得资源综合利用信息;将供电煤耗信息、脱硫效率信息、供热负荷信息和资源综合利用信息,结合输电网损信息构建综合节能减排数学模型;采用多目标线性规划算法优化求解,获得次日计划自动生成。本发明对供电煤耗、脱硫效率、供热负荷、资源综合利用、电网网损等因素采用了不同的处理办法,在目标函数中考虑机组的供电煤耗、脱硫效率以及电网网损,在约束条件中考虑供热负荷、资源综合利用,从而实现节能调度的日计划自动生成。
申请公布号 CN102930351B 申请公布日期 2016.01.20
申请号 CN201210419520.7 申请日期 2012.10.26
申请人 安徽省电力公司;安徽南瑞继远软件有限公司 发明人 杜贵和;王正风;李端超;陈实;宋祥春;朱刚刚;马金辉;胡健;汪春燕;唐超
分类号 G06Q10/04(2012.01)I;G06Q50/06(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 合肥天明专利事务所 34115 代理人 吴娜
主权项 一种综合节能优化调度日计划生成方法,该方法包括下列顺序的步骤:(1)通过煤耗在线监测系统获得机组供电煤耗信息,通过烟气排放在线监测系统获得脱硫效率信息,通过热电机组在线监测系统获得供热负荷信息,通过资源综合利用在线监测系统获得资源综合利用信息;(2)将上述供电煤耗信息、脱硫效率信息、供热负荷信息和资源综合利用信息,结合输电网损信息构建综合节能减排数学模型;(3)采用多目标线性规划算法优化求解,获得次日计划自动生成;所述的基于供电煤耗信息、脱硫效率信息、供热负荷信息、资源综合利用信息和输电网损信息构建的综合节能减排数学模型,其目标函数如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mi>min</mi></mtd><mtd><mrow><mi>F</mi><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo>&lsqb;</mo><msub><mi>&lambda;</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>C</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>p</mi><mo>(</mo><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow><mo>)</mo><mo>)</mo></mrow><mo>&rsqb;</mo><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mi>&Delta;P</mi><mi>t</mi></msub><mo>&CenterDot;</mo><mi>&beta;</mi></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000771575470000011.GIF" wi="1700" he="171" /></maths>其中,T为调度期间的时段数;N为机组数;p(i,t)为机组i在时段t的发电出力;Ci(p(i,t))为机组i在t时的能耗;△P<sub>t</sub>表示t时段的网损,β表示电量与煤耗的折算系数;λ<sub>i</sub>表示机组脱硫状态系数;所述的目标函数满足系统负荷平衡约束如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>p</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>p</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000771575470000012.GIF" wi="1501" he="164" /></maths>式中,p(i,t)为机组i在时段t的发电出力,p<sub>d</sub>(t)为t时的系统发电口径净负荷;所述的目标函数满足机组旋转备用约束如下:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mover><mi>p</mi><mo>&OverBar;</mo></mover><mi>r</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&GreaterEqual;</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>r</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>&GreaterEqual;</mo><munder><msub><mi>p</mi><mi>r</mi></msub><mo>&OverBar;</mo></munder><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000771575470000013.GIF" wi="1578" he="174" /></maths>式中,r(i,t)为机组i在t时提供的旋转备用;<img file="FDA0000771575470000014.GIF" wi="168" he="91" /><u>p<sub>r</sub></u>(t)分别表示系统t时的上调旋转备用需求与下调旋转备用需求;所述的目标函数满足电网安全约束,其中满足的支路潮流约束如下:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><munder><msub><mi>p</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&OverBar;</mo></munder><mo>&le;</mo><msub><mi>p</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&le;</mo><mover><msub><mi>p</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000771575470000015.GIF" wi="1355" he="130" /></maths>式中,<img file="FDA0000771575470000016.GIF" wi="114" he="100" /><u>p<sub>ij</sub></u>分别表示支路ij的潮流功率上、下限值;满足的联络线断面潮流约束如下:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><munder><msub><mi>P</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&OverBar;</mo></munder><mo>&le;</mo><msub><mi>p</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&le;</mo><mover><msub><mi>P</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000771575470000021.GIF" wi="1392" he="124" /></maths>式中,<img file="FDA0000771575470000022.GIF" wi="114" he="115" /><u>p<sub>ij</sub></u>分别表示联络线断面ij的潮流功率上、下限值;所述的目标函数满足机组运行约束如下:p<sub>i,min</sub>u(i,t)≤p(i,t)≤p<sub>i,max</sub>u(i,t)   (6)式中,p(i,t)为机组i在时段t的发电出力,p<sub>i,min</sub>u(i,t),p<sub>i,max</sub>u(i,t)分别表示发电机组i在时段t的运行状态输出功率的上、下限;所述的目标函数满足固定出力约束,固定出力约束是依据电网运行的特点以及机组的特性确定所需要满足的约束条件;对于热定电机组,依据热电机组的供热负荷需求,确定发电计划;对于综合利用机组,若为煤矸石型综合利用机组,当前日低位发热量≤3500kCal/kg时,按固定出力发电,否则参与常规火电优化调度确定出力;若为燃用生活垃圾型综合利用机组,当前日燃料比≥80%时,按固定出力发电,否则参与常规火电优化调度确定出力;对于风电、太阳能新能源机组,按预测出力确定发电计划;所述的目标函数满足最小运行时间和最小停运时间约束如下:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mo>(</mo><msubsup><mi>V</mi><mrow><mi>t</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi>o</mi><mi>n</mi></mrow></msubsup><mo>-</mo><msubsup><mi>T</mi><mi>i</mi><mrow><mi>min</mi><mo>_</mo><mi>o</mi><mi>n</mi></mrow></msubsup><mo>)</mo><mo>&CenterDot;</mo><mo>(</mo><mi>u</mi><mo>(</mo><mi>i</mi><mo>,</mo><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo><mo>-</mo><mi>u</mi><mo>(</mo><mi>i</mi><mo>,</mo><mi>t</mi><mo>)</mo><mo>)</mo><mo>&GreaterEqual;</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000771575470000023.GIF" wi="1536" he="106" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mo>(</mo><msubsup><mi>V</mi><mrow><mi>t</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi>o</mi><mi>f</mi><mi>f</mi></mrow></msubsup><mo>-</mo><msubsup><mi>T</mi><mi>i</mi><mrow><mi>min</mi><mo>_</mo><mi>o</mi><mi>f</mi><mi>f</mi></mrow></msubsup><mo>)</mo><mo>&CenterDot;</mo><mo>(</mo><mi>u</mi><mo>(</mo><mi>i</mi><mo>,</mo><mi>t</mi><mo>)</mo><mo>-</mo><mi>u</mi><mo>(</mo><mi>i</mi><mo>,</mo><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo><mo>)</mo><mo>&GreaterEqual;</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000771575470000024.GIF" wi="1537" he="99" /></maths>式中,T<sub>imin_on</sub>和T<sub>imin_off</sub>分别为机组i的最小开机时间和最小停机时间;V<sub>t,ion</sub>和V<sub>t,ioff</sub>分别为机组i在t时段之前的连续开机和停机时间;u(i,t)表示机组i在时段t的启停方式,运行或停止,1表示运行,0表示停止;所述的目标函数满足机组调节备用约束:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mover><msubsup><mi>p</mi><mi>r</mi><mo>&prime;</mo></msubsup><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&GreaterEqual;</mo><munder><mo>&Sigma;</mo><mrow><mi>i</mi><mo>&Element;</mo><msub><mi>I</mi><mi>g</mi></msub></mrow></munder><msup><mi>r</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>&GreaterEqual;</mo><munder><msubsup><mi>p</mi><mi>r</mi><mo>&prime;</mo></msubsup><mo>&OverBar;</mo></munder><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000771575470000025.GIF" wi="1517" he="162" /></maths>式中,p'<sub>r</sub>(t)、<u>p</u><sub>r</sub><u>′</u>(t)分别表示系统t时的AGC上调备用需求和下调备用需求,r'(i,t)表示机组i在t时提供的AGC备用;所述的目标函数满足机组最大开停次数约束:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><mi>y</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>N</mi><mi>s</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000771575470000031.GIF" wi="1527" he="185" /></maths>式中,Ns表示调度期内最大开停次数;y(i,t)为机组i在时段t是否有停机到开机状态变化的标志。
地址 230061 安徽省合肥市包河区黄山路9号