发明名称 一种考虑风电预测误差的机组有功实时调度方法
摘要 本发明公开了一种考虑风电预测误差的电力系统有功实时经济调度方法。使用拉普拉斯分布描述风电预测误差随机分布特性,以经济性最优与弃风最小为双目标,引入机会约束条件,基于拉丁超立方采样的概率潮流方法求解系统有功潮流概率分布,建立考虑风电预测误差的缓冲机组有功最优分配模型,并采用改进的遗传算法求解该机会约束规划模型。本发明在实时调度中考虑风电预测误差,能够避免风电出力偏差给系统带来的功率失衡和潮流越限,保障系统的安全运行。此外系统在消纳风电预测误差的同时,提高了风电接纳能力。
申请公布号 CN104810863A 申请公布日期 2015.07.29
申请号 CN201510236708.1 申请日期 2015.05.11
申请人 东南大学 发明人 蒋平;张文婷;霍雨翀
分类号 H02J3/48(2006.01)I 主分类号 H02J3/48(2006.01)I
代理机构 南京瑞弘专利商标事务所(普通合伙) 32249 代理人 严巧巧
主权项 一种考虑风电预测误差的机组有功实时调度方法,其特征在于:包括顺序执行的以下步骤:步骤一、读取待调度电力系统的网络数据和机组参数,以及最新滚动调度后电力系统中的机组出力;步骤二、设定缓冲机组有功最优分配模型目标函数为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>min</mi><mrow><mo>(</mo><munder><mi>&Sigma;</mi><mrow><mi>v</mi><mo>&Element;</mo><msub><mi>N</mi><mi>rt</mi></msub></mrow></munder><msubsup><mi>a</mi><mi>v</mi><mi>rt</mi></msubsup><mo>&CenterDot;</mo><mi>&Delta;</mi><msubsup><mi>P</mi><mi>v</mi><mi>rt</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><munder><mi>&Sigma;</mi><mrow><mi>v</mi><mo>&Element;</mo><msub><mi>N</mi><mi>wind</mi></msub></mrow></munder><msub><mi>&kappa;</mi><mi>v</mi></msub><mo>&CenterDot;</mo><mi>&Delta;</mi><msubsup><mi>P</mi><mi>v</mi><mi>loss</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000714599080000011.GIF" wi="1253" he="117" /></maths>(1)式中:v表示机组包括缓冲机组和风电机组,N<sub>rt</sub>是缓冲机组的集合,N<sub>wind</sub>是风电机组的集合,<img file="FDA0000714599080000012.GIF" wi="67" he="74" />是缓冲机组发电出力单位调整成本,<img file="FDA0000714599080000013.GIF" wi="174" he="77" />是缓冲机组在t时刻初的出力调整量,<img file="FDA0000714599080000014.GIF" wi="198" he="78" />是风电机组在t时刻初的弃风量,κ<sub>v</sub>为弃风权重因子;步骤三、获取下一周期的风电预测出力及总负荷预测信息,并设置阈值P<sub>wind</sub>;利用下一周期的风电预测出力减去当前周期的风电出力获得风电总出力预测增量,判断风电总出力预测增量是否大于P<sub>wind</sub>;若大于等于P<sub>wind</sub>,则下一周期电力系统选择弃风;若小于P<sub>wind</sub>,则下一周期电力系统选择全部接纳风电出力;步骤四、设定<img file="FDA0000714599080000015.GIF" wi="170" he="79" />和<img file="FDA0000714599080000016.GIF" wi="195" he="82" />为决策变量,组成待求染色体;步骤五、随机产生初代染色体,基于拉丁超立方采样获取风电预测误差样品,通过概率潮流方法得到电力系统有功潮流概率分布,验证其是否满足如下4个约束条件;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><munder><mi>&Sigma;</mi><mrow><mi>v</mi><mo>&Element;</mo><msub><mi>N</mi><mi>rt</mi></msub></mrow></munder><mi>&Delta;</mi><msubsup><mi>P</mi><mi>v</mi><mi>rt</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><munder><mi>&Sigma;</mi><mrow><mi>v</mi><mo>&Element;</mo><msub><mi>N</mi><mi>wind</mi></msub></mrow></munder><mrow><mo>(</mo><mi>&Delta;</mi><msubsup><mi>P</mi><mi>v</mi><mi>wind</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>&Delta;</mi><msubsup><mi>P</mi><mi>v</mi><mi>loss</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mi>&Delta;</mi><msub><mi>P</mi><mi>load</mi></msub><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000714599080000017.GIF" wi="1533" he="157" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>Pr</mi><mo>{</mo><munder><msub><mi>P</mi><mi>Tk</mi></msub><mo>&OverBar;</mo></munder><mo>&le;</mo><msub><mi>P</mi><mi>Tk</mi></msub><mrow><mo>(</mo><msubsup><mi>&epsiv;</mi><mi>v</mi><mi>f</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>&le;</mo><munder><mover><msub><mi>P</mi><mi>Tk</mi></msub><mo>&OverBar;</mo></mover><mo>&OverBar;</mo></munder><mo>}</mo><mo>&GreaterEqual;</mo><mn>1</mn><mo>-</mo><msubsup><mi>&theta;</mi><mn>1</mn><mi>rt</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000714599080000018.GIF" wi="1430" he="96" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mo>-</mo><mi>&Delta;</mi><munder><msubsup><mi>P</mi><mi>v</mi><mi>rt</mi></msubsup><mo>&OverBar;</mo></munder><mo>&le;</mo><mi>&Delta;</mi><msubsup><mi>P</mi><mi>v</mi><mi>rt</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&le;</mo><mi>&Delta;</mi><msubsup><mi>P</mi><mi>v</mi><mi>rt</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000714599080000019.GIF" wi="1436" he="90" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>Pr</mi><mo>{</mo><mo>-</mo><munder><mi>&Sigma;</mi><mrow><mi>v</mi><mo>&Element;</mo><msub><mi>N</mi><mi>rt</mi></msub></mrow></munder><mi>&Delta;</mi><munder><msubsup><mi>P</mi><mi>v</mi><mi>rt</mi></msubsup><mo>&OverBar;</mo></munder><mo>&le;</mo><munder><mi>&Sigma;</mi><mrow><mi>v</mi><mo>&Element;</mo><msub><mi>N</mi><mi>rt</mi></msub></mrow></munder><mi>&Delta;</mi><msubsup><mi>P</mi><mi>v</mi><mi>rt</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><munder><mi>&Sigma;</mi><mrow><mi>v</mi><mo>&Element;</mo><msub><mi>N</mi><mi>wind</mi></msub></mrow></munder><mrow><mo>(</mo><msubsup><mi>&epsiv;</mi><mi>v</mi><mi>f</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>&Delta;</mi><msubsup><mi>P</mi><mi>v</mi><mi>loss</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>&le;</mo><munder><mi>&Sigma;</mi><mrow><mi>v</mi><mo>&Element;</mo><msub><mi>N</mi><mi>rt</mi></msub></mrow></munder><mi>&Delta;</mi><mover><msubsup><mi>P</mi><mi>v</mi><mi>rt</mi></msubsup><mo>&OverBar;</mo></mover><mo>}</mo><mo>&le;</mo><mn>1</mn><mo>-</mo><msubsup><mi>&theta;</mi><mn>2</mn><mi>rt</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00007145990800000110.GIF" wi="1662" he="115" /></maths>式中:(2)式中:<img file="FDA00007145990800000111.GIF" wi="214" he="75" />是风电机组预测出力的增量,ΔP<sub>load</sub>为总负荷预测的增量,总负荷预测增量为下一周期总负荷预测信息减去当前周期总负荷;(3)式中:<img file="FDA00007145990800000112.GIF" wi="125" he="75" />为风电机组的风电预测误差量,P<sub>Tk</sub>为t时刻实时调度风电预测误差后的线路断面k的有功功率概率分布,<img file="FDA00007145990800000113.GIF" wi="68" he="77" />为线路断面k的有功功率传输上限,<img file="FDA00007145990800000114.GIF" wi="65" he="69" />为线路断面k的有功功率传输下限,Pr{·}表示概率,<img file="FDA00007145990800000115.GIF" wi="63" he="67" />为给定的潮流约束置信水平,且<img file="FDA00007145990800000116.GIF" wi="268" he="79" />该式需对电网中任一断面k都成立;(4)式中:<img file="FDA00007145990800000117.GIF" wi="126" he="101" />为缓冲机组出力调整的下限,且<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>&Delta;</mi><munder><msubsup><mi>P</mi><mi>v</mi><mi>rt</mi></msubsup><mo>&OverBar;</mo></munder><mo>=</mo><mi>min</mi><mrow><mo>(</mo><msubsup><mi>P</mi><mi>v</mi><mi>rt</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>P</mi><mrow><mi>v</mi><mo>,</mo><mi>min</mi></mrow></msub><mo>,</mo><mi>&Delta;</mi><msubsup><mi>P</mi><mi>v</mi><mi>dn</mi></msubsup><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA00007145990800000118.GIF" wi="755" he="91" /></maths>其中P<sub>v,min</sub>为缓冲机组出力下限,<img file="FDA00007145990800000119.GIF" wi="129" he="87" />为缓冲机组出力的最大下调量,<img file="FDA00007145990800000120.GIF" wi="116" he="92" />为缓冲机组出力调整的上限,且<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>&Delta;</mi><mover><msubsup><mi>P</mi><mi>v</mi><mi>rt</mi></msubsup><mo>&OverBar;</mo></mover><mo>=</mo><mi>min</mi><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>v</mi><mo>,</mo><mi>max</mi></mrow></msub><mo>-</mo><msubsup><mi>P</mi><mi>v</mi><mi>rt</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mi>&Delta;</mi><msubsup><mi>P</mi><mi>v</mi><mi>up</mi></msubsup><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000714599080000021.GIF" wi="751" he="96" /></maths>其中P<sub>v,max</sub>为缓冲机组出力上限,<img file="FDA0000714599080000022.GIF" wi="130" he="86" />为缓冲机组出力的最大上调量;(5)式中:<img file="FDA0000714599080000023.GIF" wi="85" he="82" />为给定的旋转备用约束置信水平;若满足所有约束条件,则将对应的初代染色体放入初代种群;否则重新产生染色体并在新的染色体的基础上重复上述采样及验证方法,将满足约束条件的对应的染色体放入初代种群中,直至达到种群规模N;步骤六:经由改进遗传算法进行保留、选择、交叉、变异、复制操作,采用基于序的评价函数,对新产生种群进行检验以确认其满足步骤五中的4个约束条件,直到改进遗传算法达到指定进化代数d或满足以下收敛性条件:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><munder><mi>&Sigma;</mi><mrow><mi>v</mi><mo>&Element;</mo><mi>N</mi></mrow></munder><mo>|</mo><msubsup><mi>&eta;</mi><mi>v</mi><mi>d</mi></msubsup><mo>-</mo><msubsup><mi>&eta;</mi><mi>v</mi><mrow><mi>d</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>|</mo><mo>&le;</mo><mi>&epsiv;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000714599080000024.GIF" wi="1030" he="144" /></maths>(1)式中:<img file="FDA0000714599080000025.GIF" wi="76" he="87" />为染色体,d为进化代数,ε为收敛判据;步骤七:将最终代染色体中的最优个体作为优化后的缓冲机组出力调整量和弃风量,并按照优化后缓冲机组出力调整量和弃风量进行实时调度。
地址 211189 江苏省南京市江宁区东南大学路2号