发明名称 一种集群风电场有功功率控制方法
摘要 本发明公开了一种集群风电场有功功率控制方法,是一种计算省(自治区)级调度机构下发的有功功率分配到各风电场的计划值的方法,采用本发明提出的有功功率分配计算方法可以确定下一控制周期各大集群风电场的有功出力计划值,通过实时计算电网总接纳风电能力,从而保证了在电力系统安全稳定运行的前提下,各集群风电场出力最大化。本发明考虑电网总接纳能力的变化、各风电场当前出力、风电预测情况以及负荷预测情况等多方面因素,每固定周期计算一次各风场计划,并下发至各风场,该计划为上限值,低于该计划风场可自由发电,高于该计划需要申请。本发明最大限度的提高风电场的发电效率,并保证电力系统最大化接纳风电。
申请公布号 CN103219751B 申请公布日期 2015.04.29
申请号 CN201310164592.6 申请日期 2013.05.08
申请人 东南大学 发明人 高赐威;陈曦寒
分类号 H02J3/46(2006.01)I 主分类号 H02J3/46(2006.01)I
代理机构 南京瑞弘专利商标事务所(普通合伙) 32249 代理人 杨晓玲
主权项 一种集群风电场有功功率控制方法,其特征在于,该方法包括以下步骤:1)集群风电有功功率控制系统的通信终端从风电场调度中心采集当前第i个风电场的出力<img file="FDA0000622520430000011.GIF" wi="134" he="115" />第i个风电场的最大出力为<img file="FDA0000622520430000012.GIF" wi="153" he="107" />第i个风电场功率预测装置预测的下一控制周期风电场的有功出力预测值<img file="FDA0000622520430000013.GIF" wi="143" he="133" />从省(自治区)级调度机构采集电力系统当前负荷数据<img file="FDA0000622520430000014.GIF" wi="100" he="95" />当前下发给各风电场的有功功率计划值<img file="FDA0000622520430000015.GIF" wi="137" he="99" />以及负荷预测系统预测的下一控制周期系统负荷预测值<img file="FDA0000622520430000016.GIF" wi="138" he="126" />其中t表示当前控制周期,t+1表示下一控制周期,i为风电场中的序号,W表示风电场数据,L表示负荷数据;2)根据以下目标函数和约束条件,通过混合整数规划法计算得到1到n个风电场下一控制周期的并网系数<img file="FDA0000622520430000017.GIF" wi="559" he="100" />所述目标函数为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>max</mi></mtd><mtd><mi>F</mi><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>K</mi><mi>i</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>|</mo><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>plan</mi></mrow><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>W</mi></mrow><mi>t</mi></msubsup><mo>|</mo></mrow></mfrac></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000622520430000018.GIF" wi="804" he="339" /></maths>式中:<img file="FDA0000622520430000019.GIF" wi="122" he="98" />为第i个风电场下一控制周期并网系数,n为风电场的个数,<img file="FDA00006225204300000110.GIF" wi="150" he="108" />为第i个风电场下一控制周期的有功出力计划值,即<img file="FDA00006225204300000111.GIF" wi="536" he="135" /><img file="FDA00006225204300000112.GIF" wi="102" he="96" />为当前第i个风电场的出力;所述约束条件包括系统有功功率平衡约束、备用约束、火电机组出力上下限约束、火电机组开停机时间约束、火电机组爬坡约束、线路潮流约束、风电场出力约束、风电场功率变化约束;所述系统有功功率平衡约束如下式所示,表示各个时刻开机火电机组的出力之和应等于负荷需求与风电场出力之差:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><munderover><mi>&Sigma;</mi><mi>j</mi><mi>m</mi></munderover><msubsup><mi>P</mi><mi>j</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&CenterDot;</mo><msubsup><mi>U</mi><mi>j</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mover><msubsup><mi>P</mi><mi>L</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>^</mo></mover><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msubsup><mi>K</mi><mi>i</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&times;</mo><mover><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>W</mi></mrow><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>^</mo></mover><mo>)</mo></mrow><mo>,</mo><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA00006225204300000113.GIF" wi="1154" he="173" /></maths>式中:m为火电机组数,n为风电场数,<img file="FDA00006225204300000114.GIF" wi="121" he="116" />为机组j在下一控制周期的有功出力,<img file="FDA0000622520430000021.GIF" wi="120" he="105" />为火电机组j在下一控制周期的启停状态,<img file="FDA0000622520430000022.GIF" wi="207" he="102" />即关机状态,<img file="FDA0000622520430000023.GIF" wi="199" he="104" />为开机状态,<img file="FDA0000622520430000024.GIF" wi="115" he="126" />为下一控制周期时段的负荷预测值,<img file="FDA0000622520430000025.GIF" wi="117" he="134" />为第i个风电场下一控制周期的风电场出力预测值;所述备用约束如下式所示,开机火电机组所提供的备用应能够满足负荷与风电预测双重误差所带来的影响:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><munderover><mi>&Sigma;</mi><mi>j</mi><mi>m</mi></munderover><msub><mi>P</mi><mrow><mi>j</mi><mo>,</mo><mi>max</mi></mrow></msub><mo>&CenterDot;</mo><msubsup><mi>U</mi><mi>j</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>P</mi><mi>D</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&GreaterEqual;</mo><msubsup><mi>R</mi><mi>L</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msubsup><mi>K</mi><mi>i</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&times;</mo><msubsup><mi>R</mi><mrow><mi>i</mi><mo>,</mo><mi>W</mi></mrow><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>,</mo><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000622520430000026.GIF" wi="1291" he="163" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>P</mi><mi>D</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>-</mo><munderover><mi>&Sigma;</mi><mi>j</mi><mi>m</mi></munderover><msub><mi>P</mi><mrow><mi>j</mi><mo>,</mo><mi>max</mi></mrow></msub><mo>&CenterDot;</mo><msubsup><mi>U</mi><mi>j</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&le;</mo><msubsup><mi>R</mi><mi>L</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msubsup><mi>K</mi><mi>i</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&times;</mo><msubsup><mi>R</mi><mrow><mi>i</mi><mo>,</mo><mi>W</mi></mrow><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>,</mo><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000622520430000027.GIF" wi="1281" he="159" /></maths>式中:P<sub>j,min</sub>为第j个机组最大出力,P<sub>j,max</sub>为第j个机组最小出力,<img file="FDA0000622520430000028.GIF" wi="68" he="77" />为t时段的负荷减去t时段的风电出力的净负荷;<img file="FDA0000622520430000029.GIF" wi="68" he="80" />为t时段负荷预测误差所要求的备用,<img file="FDA00006225204300000210.GIF" wi="76" he="80" />为t时段风电预测误差所要求的备用;所述火电机组出力上下限约束如下式所示:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>j</mi><mo>,</mo><mi>max</mi></mrow></msub><mo>&CenterDot;</mo><msubsup><mi>U</mi><mi>j</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&le;</mo><msubsup><mi>P</mi><mi>j</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&le;</mo><msub><mi>P</mi><mrow><mi>j</mi><mo>,</mo><mi>max</mi></mrow></msub><mo>&CenterDot;</mo><msubsup><mi>U</mi><mi>j</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>,</mo><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA00006225204300000211.GIF" wi="991" he="110" /></maths>火电机组开停机时间约束如下式所示,因机组不能频繁启停,机组在一定时间内开停次数应限制在设定范围以内;<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mrow><mo>(</mo><msubsup><mi>X</mi><mi>j</mi><mi>on</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>T</mi><mi>j</mi><mi>on</mi></msubsup><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><msubsup><mi>U</mi><mi>j</mi><mi>t</mi></msubsup><mo>-</mo><msubsup><mi>U</mi><mi>j</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA00006225204300000212.GIF" wi="1023" he="101" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mrow><mo>(</mo><msubsup><mi>X</mi><mi>j</mi><mi>off</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>T</mi><mi>j</mi><mi>off</mi></msubsup><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><msubsup><mi>U</mi><mi>j</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>U</mi><mi>j</mi><mi>t</mi></msubsup><mo>)</mo></mrow><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA00006225204300000213.GIF" wi="1061" he="108" /></maths>式中,<img file="FDA00006225204300000214.GIF" wi="122" he="103" />和<img file="FDA00006225204300000215.GIF" wi="121" he="111" />分别表示火电机组j到当前累计开机时间和累计关机时间,<img file="FDA00006225204300000216.GIF" wi="140" he="105" /><img file="FDA00006225204300000217.GIF" wi="113" he="117" />分别为机组j的最小运行时间和最小停运时间;所述火电机组爬坡约束如下式所示,机组出力在前后两个时段应满足机组爬坡约束:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msubsup><mi>P</mi><mi>j</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>P</mi><mi>j</mi><mi>t</mi></msubsup><mo>&le;</mo><msub><mi>R</mi><mi>u</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA00006225204300000218.GIF" wi="462" he="109" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msubsup><mi>P</mi><mi>j</mi><mi>t</mi></msubsup><mo>-</mo><msubsup><mi>P</mi><mi>j</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&le;</mo><msub><mi>R</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA00006225204300000219.GIF" wi="461" he="103" /></maths>式中,R<sub>u</sub>(j)和R<sub>d</sub>(j)分别代表火电机组j的上坡速率限值和下坡速率限值;所述线路潮流约束如下式所示,采用直流潮流来表示网络约束对机组组合的作用:<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msubsup><mi>P</mi><mi>k</mi><mi>min</mi></msubsup><mo>&le;</mo><msubsup><mi>P</mi><mi>k</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&le;</mo><msubsup><mi>P</mi><mi>k</mi><mi>max</mi></msubsup></mrow>]]></math><img file="FDA0000622520430000031.GIF" wi="518" he="108" /></maths>式中,<img file="FDA0000622520430000032.GIF" wi="142" he="106" />和<img file="FDA0000622520430000033.GIF" wi="132" he="101" />分别为线路k的最小传输容量和最大传输容量;所述风电场出力约束如下式所示:<maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>W</mi></mrow><mi>min</mi></msubsup><mo>&le;</mo><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>plan</mi></mrow><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&le;</mo><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>W</mi></mrow><mi>max</mi></msubsup></mrow>]]></math><img file="FDA0000622520430000034.GIF" wi="499" he="106" /></maths>式中:<img file="FDA0000622520430000035.GIF" wi="124" he="92" />为第i个风电场的最大出力;<img file="FDA0000622520430000036.GIF" wi="122" he="104" />为第i个风电场的最小出力;所述风电场功率变化约束如下式所示,是根据《风电场接入电网技术规定》,确保风电场功率变化率不超过电网调度部门的给定值:<maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><mo>|</mo><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>plan</mi></mrow><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>W</mi></mrow><mi>t</mi></msubsup><mo>|</mo><mo>&le;</mo><msub><mi>&Delta;P</mi><mrow><mi>i</mi><mo>,</mo><mi>rule</mi></mrow></msub></mrow>]]></math><img file="FDA0000622520430000037.GIF" wi="548" he="108" /></maths>式中:ΔP<sub>i,rule</sub>为电网调度部门的规定的第i个风电场输出功率变化率给定值;3)根据下式计算下一控制周期第i个风电场有功出力计划值:<maths num="0013" id="cmaths0013"><math><![CDATA[<mrow><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>W</mi></mrow><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>W</mi></mrow><mi>max</mi></msubsup><mo>,</mo></mtd><mtd><msubsup><mi>K</mi><mi>i</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&times;</mo><mover><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>W</mi></mrow><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>^</mo></mover><mo>&GreaterEqual;</mo><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>W</mi></mrow><mi>max</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>K</mi><mi>i</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&times;</mo><mover><mrow><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>W</mi></mrow><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>,</mo></mrow><mo>^</mo></mover></mtd><mtd><msubsup><mi>K</mi><mi>i</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&times;</mo><mrow><mover><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>W</mi></mrow><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>^</mo></mover><mo>&lt;</mo><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>W</mi></mrow><mi>max</mi></msubsup></mrow></mtd></mtr></mtable></mfenced><mo>.</mo></mrow>]]></math><img file="FDA0000622520430000038.GIF" wi="881" he="254" /></maths>
地址 211189 江苏省南京市江宁区东南大学路2号