发明名称 含风电系统机组的备用配置优化方法
摘要 本发明公开了一种含风电系统机组的备用配置优化方法。本发明方法提出了一种利用风电和负荷联合概率密度函数的积分式计算一定正负旋转备用容量下的停电损失和弃风损失两种风险成本的计算方法;给出系统总正负旋转备用容量分配到不同煤耗机组之间的最优分配方法;给出综合考虑风险成本和发电成本,电力系统最经济运行时系统总正负旋转备用容量、各机组发电计划及备用分配容量的决策方法。本发明适用于已接入风电场的电力系统,可为系统决策出兼顾经济性和安全性的总正负旋转备用容量、各机组出力计划及最优正负旋转备用分配容量,同时备用容量在各机组的分配考虑了机组爬坡约束和网络安全约束,保证了各机组提供的备用完全可用。
申请公布号 CN103956773B 申请公布日期 2016.08.31
申请号 CN201410196885.7 申请日期 2014.05.09
申请人 河海大学 发明人 赵晋泉;唐洁
分类号 H02J3/38(2006.01)I 主分类号 H02J3/38(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 柏尚春
主权项 一种含风电系统机组的备用配置优化方法,其包括以下步骤:步骤A、利用风电和负荷两种随机变量的联合概率密度函数的积分式计算当前正负旋转备用容量下的停电损失和弃风损失两种风险成本,具体按照以下方法:步骤A1、用正态分布模型描述风电和负荷的概率分布,其表达式如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>f</mi><msub><mi>P</mi><mi>W</mi></msub></msub><mrow><mo>(</mo><msub><mi>P</mi><mi>W</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msqrt><mrow><mn>2</mn><mi>&pi;</mi></mrow></msqrt><msub><mi>&sigma;</mi><mi>W</mi></msub></mrow></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msup><mrow><mo>(</mo><msub><mi>P</mi><mi>W</mi></msub><mo>-</mo><msub><mover><mi>P</mi><mo>&OverBar;</mo></mover><mi>W</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mrow><mn>2</mn><msubsup><mi>&sigma;</mi><mi>W</mi><mn>2</mn></msubsup></mrow></mfrac></mrow></msup></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>f</mi><mrow><msub><mi>&Delta;P</mi><mi>L</mi></msub></mrow></msub><mrow><mo>(</mo><msub><mi>&Delta;P</mi><mi>L</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msqrt><mrow><mn>2</mn><mi>&pi;</mi></mrow></msqrt><msub><mi>&sigma;</mi><mi>L</mi></msub></mrow></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mfrac><mrow><msup><msub><mi>&Delta;p</mi><mi>L</mi></msub><mn>2</mn></msup></mrow><mrow><mn>2</mn><msubsup><mi>&sigma;</mi><mi>L</mi><mn>2</mn></msubsup></mrow></mfrac></mrow></msup></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001012396040000011.GIF" wi="582" he="351" /></maths>式中:P<sub>W</sub>为实际风电功率;<img file="FDA0001012396040000012.GIF" wi="63" he="63" />为风功率预测值;ΔP<sub>L</sub>为负荷预测误差;σ<sub>W</sub>、σ<sub>L</sub>分别为风电功率和负荷预测误差的标准差;步骤A2、通过卷积公式求出二者联合分布Z=P<sub>W</sub>‑ΔP<sub>L</sub>的概率密度函数,其表达式如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>f</mi><mi>Z</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><msub><mi>P</mi><mi>N</mi></msub></msubsup><msub><mi>f</mi><msub><mi>P</mi><mi>W</mi></msub></msub><mrow><mo>(</mo><msub><mi>P</mi><mi>W</mi></msub><mo>)</mo></mrow><msub><mi>f</mi><mrow><msub><mi>&Delta;P</mi><mi>L</mi></msub></mrow></msub><mrow><mo>(</mo><mi>z</mi><mo>-</mo><msub><mi>P</mi><mi>W</mi></msub><mo>)</mo></mrow><msub><mi>dP</mi><mi>W</mi></msub></mrow>]]></math><img file="FDA0001012396040000013.GIF" wi="702" he="95" /></maths>式中:P<sub>N</sub>为风电额定功率;步骤A3、利用风电功率和负荷预测误差的联合概率密度函数,获得在当前正负旋转备用值下的停电期望和弃风期望;电力系统使用全部的正旋转备用,仍不能弥补系统功率的不足,就会产生停电,其中,停电期望的表达式如下:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>L</mi><mrow><mi>l</mi><mi>o</mi><mi>s</mi><mi>s</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>=</mo><munder><mo>&Integral;</mo><mrow><msubsup><mi>R</mi><mi>t</mi><mi>U</mi></msubsup><mo>+</mo><mi>Z</mi><mo>-</mo><msub><mover><mi>P</mi><mo>&OverBar;</mo></mover><mi>W</mi></msub><mo>&lt;</mo><mn>0</mn></mrow></munder><msub><mi>&Delta;P</mi><mrow><mi>l</mi><mi>o</mi><mi>s</mi><mi>s</mi></mrow></msub><msub><mi>f</mi><mi>Z</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mi>d</mi><mi>z</mi></mrow>]]></math><img file="FDA0001012396040000014.GIF" wi="582" he="126" /></maths>式中:ΔP<sub>loss</sub>为功率缺额,有<img file="FDA0001012396040000015.GIF" wi="467" he="63" /><img file="FDA0001012396040000016.GIF" wi="66" he="66" />为系统在时段t的正旋转备用,L<sub>loss,t</sub>为在正备用值为<img file="FDA0001012396040000017.GIF" wi="64" he="71" />时的停电期望;电力系统使用全部的负旋转备用,仍不能平衡系统功率的盈余,则会产生弃风,其中,弃风期望的表达式如下:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>W</mi><mrow><mi>l</mi><mi>o</mi><mi>s</mi><mi>s</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>=</mo><munder><mo>&Integral;</mo><mrow><msubsup><mi>R</mi><mi>t</mi><mi>D</mi></msubsup><mo>-</mo><mrow><mo>(</mo><mi>Z</mi><mo>-</mo><msub><mover><mi>P</mi><mo>&OverBar;</mo></mover><mi>W</mi></msub><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn></mrow></munder><msub><mi>&Delta;W</mi><mrow><mi>l</mi><mi>o</mi><mi>s</mi><mi>s</mi></mrow></msub><msub><mi>f</mi><mi>Z</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mi>d</mi><mi>z</mi></mrow>]]></math><img file="FDA0001012396040000018.GIF" wi="626" he="127" /></maths>式中:ΔW<sub>loss</sub>为弃风功率,有<img file="FDA0001012396040000021.GIF" wi="418" he="63" /><img file="FDA0001012396040000022.GIF" wi="69" he="62" />为系统在时段t的负备用;W<sub>loss,t</sub>为在负备用值为<img file="FDA0001012396040000023.GIF" wi="67" he="69" />时的弃风期望;步骤A4、根据停电期望和停电成本求出停电损失,同时根据弃风期望和弃风成本求出弃风损失,表达式如下:<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><msub><mi>F</mi><mrow><mi>U</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>=</mo><msub><mi>c</mi><mi>L</mi></msub><msub><mi>L</mi><mrow><mi>l</mi><mi>o</mi><mi>s</mi><mi>s</mi><mo>,</mo><mi>t</mi></mrow></msub></mtd></mtr><mtr><mtd><mrow><msub><mi>F</mi><mrow><mi>D</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>=</mo><msub><mi>c</mi><mi>W</mi></msub><msub><mi>W</mi><mrow><mi>l</mi><mi>o</mi><mi>s</mi><mi>s</mi><mo>,</mo><mi>t</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001012396040000024.GIF" wi="286" he="150" /></maths>式中:c<sub>L</sub>为停电的单位成本;F<sub>U,t</sub>为停电损失;c<sub>W</sub>为弃风的单位成本;F<sub>D,t</sub>为弃风损失;步骤B、考虑备用的可用性和节煤效益,给出系统总正负旋转备用容量确定时,备用在不同煤耗机组的分配方法;步骤C、将风险成本加入发电成本目标函数中,通过对系统总正负旋转备用容量进行寻优,找到总成本最小时的系统总正负备用需求、各机组的发电计划和正负旋转备用分配容量。
地址 211100 江苏省南京市江宁区佛城西路8号