发明名称 一种适用于大规模风电接入的电网输电裕度控制方法
摘要 本发明提供了一种适用于大规模风电接入的电网输电裕度控制方法,所述方法包括步骤1:采用混合整数规划法进行安全约束机组组合优化计算,获取满足目标函数和约束条件的常规发电机组的开停状态和出力计划;步骤2:采用概率潮流方法量化未来多时段风电随机性对输电裕度的影响,得到电网支路或输电断面的潮流概率分布范围;步骤3:判断潮流概率分布范围是否满足置信区间的要求,若不满足则返回步骤1,重新进行机组组合优化计算,调整所述出力计划。与现有技术相比,本发明提供的一种适用于大规模风电接入的电网输电裕度控制方法,充分考虑了风电的随机性,完善了短期安全约束调度计划在风电接入地区的使用效果,增强了短期计划的可执行性。
申请公布号 CN104079000B 申请公布日期 2016.04.20
申请号 CN201410332243.5 申请日期 2014.07.14
申请人 国家电网公司;中国电力科学研究院 发明人 张传成;丁强;董炜;戴赛;崔辉;许丹;黄国栋;蔡帜;朱翊;闫会;门德月
分类号 H02J3/38(2006.01)I;H02J3/46(2006.01)I 主分类号 H02J3/38(2006.01)I
代理机构 北京安博达知识产权代理有限公司 11271 代理人 徐国文
主权项 一种适用于大规模风电接入的电网输电裕度控制方法,其中包括调整常规发电机组的开停状态和出力计划来控制电网输电裕度,其特征在于,所述方法包括:步骤1:用混合整数规划法进行安全约束机组组合优化计算,获取满足目标函数和约束条件的常规发电机组的开停状态和出力计划;步骤2:用基于半不变量和Gram‑Charlier级数展开的概率潮流方法,计算风电接入时未来多时段电网支路或输电断面的潮流概率分布范围;以及步骤3:判断所述潮流概率分布范围是否满足置信区间的要求,若不满足则返回步骤1,重新进行所述安全约束机组组合优化计算,调整所述出力计划;所述步骤1中所述目标函数为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>F</mi><mo>=</mo><mi>min</mi><mrow><mo>&lsqb;</mo><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>T</mi></msub></munderover><mrow><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>t</mi></msub></munderover><mrow><msub><mi>B</mi><mi>i</mi></msub><mrow><mo>(</mo><mrow><msub><mi>P</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mi>t</mi></mrow><mo>)</mo></mrow></mrow></mrow><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>T</mi></msub></munderover><mrow><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>t</mi></msub></munderover><mrow><msub><mi>Cu</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mrow></mrow><mo>&rsqb;</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000910118320000011.GIF" wi="861" he="155" /></maths>其中,N<sub>T</sub>为发电机组数,i为发电机组序号,N<sub>t</sub>为时段数,t为时段序号,每15分钟为一个时段,则时段数N<sub>t</sub>=96;P<sub>i</sub>(t)为t时段发电机组i的有功出力值,Cu<sub>i</sub>(t)为t时段发电机组i的启机费用,B<sub>i</sub>(P<sub>i</sub>(t),t)为t时段发电机组i的发电费用;所述t时段发电机组i的有功出力值P<sub>i</sub>(t)的计算公式为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>P</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>u</mi><mrow><mi>i</mi><mo>,</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msub><munder><mi>P</mi><mo>&OverBar;</mo></munder><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>2</mn></mrow><msub><mi>N</mi><mi>L</mi></msub></munderover><mrow><mo>(</mo><mrow><msub><mi>u</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>s</mi><mi>e</mi><mi>g</mi></mrow></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000910118320000012.GIF" wi="1070" he="147" /></maths>其中,u<sub>i,j</sub>(t)为t时段发电机组i的分段费用微增率曲线j段的标志位,<img file="FDA0000910118320000013.GIF" wi="146" he="78" />为t时段发电机组i的分段费用微增率曲线j段的初始出力值,P<sub>i,j</sub>(t)为t时段发电机组i在分段费用微增率曲线j段的出力值与初始出力值<img file="FDA0000910118320000014.GIF" wi="139" he="70" />的差值,N<sub>L</sub>为分段费用微增率曲线的段数,<u>P<sub>i</sub></u>(t)为t时段发电机组i的出力下限值;所述差值P<sub>i,j</sub>(t)的限制条件为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&GreaterEqual;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mi>u</mi><mrow><mi>i</mi><mo>,</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mrow><mo>(</mo><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mn>1</mn></mrow><mrow><mi>s</mi><mi>e</mi><mi>g</mi></mrow></msubsup><mo>(</mo><mi>t</mi><mo>)</mo><mo>-</mo><munder><msub><mi>P</mi><mi>i</mi></msub><mo>&OverBar;</mo></munder><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mi>u</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mrow><mo>(</mo><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>s</mi><mi>e</mi><mi>g</mi></mrow></msubsup><mo>(</mo><mi>t</mi><mo>)</mo><mo>-</mo><msubsup><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn></mrow><mrow><mi>s</mi><mi>e</mi><mi>g</mi></mrow></msubsup><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000910118320000015.GIF" wi="693" he="262" /></maths>所述标志位<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>L</mi></msub></munderover><msub><mi>u</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&le;</mo><mn>1</mn><mo>;</mo></mrow>]]></math><img file="FDA0000910118320000016.GIF" wi="271" he="141" /></maths>若标志位<img file="FDA0000910118320000017.GIF" wi="228" he="142" />则发电机组为开机状态,若标志位<img file="FDA0000910118320000018.GIF" wi="236" he="150" />为停机状态。
地址 100031 北京市西城区西长安街86号