发明名称 一种考虑动态频率约束的含高渗透率光伏电源的孤立电网机组组合优化方法
摘要 本发明涉及一种考虑动态频率约束的含高渗透率光伏电源的孤立电网机组组合优化方法,其在机组组优化合模型中考虑动态频率特性约束,并且让光伏电源减载参与调频来增强系统的调频能力。本发明首先基于含光伏电源调频的孤立电网在故障后动态频率最低点的表达式,并公开了光伏电源参与调频时所需最小调频容量关于动态频率最低点的表达式,用以确定机组组合模型在优化过程留有充足但不过量的光伏电源调频容量。并建立了含动态频率最低点以及光伏电源调频容量约束的机组组合优化模型。本发明提出基于分解思想产生优化割与Benders分解一起对所提优化模型进行内外双层优化的方法,有效降低了问题求解复杂度。
申请公布号 CN105303267A 申请公布日期 2016.02.03
申请号 CN201510819729.6 申请日期 2015.11.20
申请人 武汉大学;湖北工业大学 发明人 林涛;叶婧;陈汝斯;徐遐龄;洪雯;毕如玉
分类号 G06Q10/04(2012.01)I;G06Q50/06(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 武汉科皓知识产权代理事务所(特殊普通合伙) 42222 代理人 鲁力
主权项 一种考虑动态频率约束的含高渗透率光伏电源的孤立电网机组组合优化方法,其特征在于,基于一个光伏电源并网下考虑动态频率约束机组组合的优化规划模型,该模型基于以下目标函数以及约束条件:所述目标函数定义为系统发电成本最低,基于:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>F</mi><mo>=</mo><mi>min</mi><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><mo>{</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo>&lsqb;</mo><msub><mi>U</mi><mrow><mi>i</mi><mi>h</mi></mrow></msub><mi>f</mi><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>i</mi><mi>h</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>U</mi><mrow><mi>i</mi><mi>h</mi></mrow></msub><mrow><mo>(</mo><mo>-</mo><msub><mi>U</mi><mrow><mi>i</mi><mrow><mo>(</mo><mi>h</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></msub><mo>)</mo></mrow><msub><mi>S</mi><mrow><mi>i</mi><mi>h</mi></mrow></msub><mo>&rsqb;</mo><mo>}</mo></mrow>]]></math><img file="FDA0000853511620000011.GIF" wi="716" he="159" /></maths>   式一式中N为常规机组的台数,T为研究周期内小时数,P<sub>ih</sub>为常规机组i在时段h的输出有功功率,f(P<sub>ih</sub>)为常规机组i的运行成本,U<sub>ih</sub>为常规机组i在时段h的启停机状态,U<sub>ih</sub>=1表示运行,U<sub>ih</sub>=0表示停机,用式表示为<img file="FDA0000853511620000012.GIF" wi="502" he="95" />a<sub>i</sub>、b<sub>i</sub>、c<sub>i</sub>为成本函数的系数;S<sub>ih</sub>为机组i在时段h的启动成本;所述约束条件为:常规机组发电、光伏发电与负荷功率平衡约束:<img file="FDA0000853511620000013.GIF" wi="476" he="151" />式二;系统旋转备用容量<img file="FDA0000853511620000014.GIF" wi="494" he="140" />式三;<img file="FDA0000853511620000015.GIF" wi="518" he="143" />式四;机组有功出力上下限约束P<sub>i</sub><sup>min</sup>≤P<sub>ih</sub>≤P<sub>i</sub><sup>max</sup>式五;机组有功功率爬坡速率约束D<sub>i</sub>≤P<sub>ih</sub>‑P<sub>i(h‑1)</sub>≤L<sub>i</sub>式六;最小启停时间约束<img file="FDA0000853511620000016.GIF" wi="205" he="78" />式七、<img file="FDA0000853511620000017.GIF" wi="208" he="82" />式八;动态频率最低点约束f<sub>min‑h</sub>≥49.0Hz式九;光伏调频备用容量约束<img file="FDA0000853511620000018.GIF" wi="383" he="79" />式十;其中<img file="FDA0000853511620000019.GIF" wi="247" he="82" />分别为第j个光伏电源、第h时段的预测值和调度值,;P<sub>Lh</sub>为第h时段负荷预测值大小;P<sub>i</sub><sup>max</sup>、P<sub>i</sub><sup>min</sup>分别为机组i的出力上、下限;γ<sub>j</sub>为第j个光伏电源预测误差百分比;L<sub>i</sub>、D<sub>i</sub>为第i台机组的爬坡上、下限;<img file="FDA0000853511620000021.GIF" wi="221" he="87" />分别为第i台机组在h时段连续开、停时间;UT<sub>i</sub>、DT<sub>i</sub>分别为第i台机组最小开、停时间;f<sub>min‑h</sub>为第h时段,系统发电功率损失后动态频率的最小值;<img file="FDA0000853511620000022.GIF" wi="121" he="70" />为第j个光伏电源、在h时段的最小调频容量;优化方法具体包括以下步骤:步骤1:获取电力系统的每个发电机组的机组特性数据P<sub>i</sub><sup>max</sup>、P<sub>i</sub><sup>min</sup>、L<sub>i</sub>、D<sub>i</sub>、UT<sub>i</sub>、DT<sub>i</sub>;负荷预测数据P<sub>ih</sub>;光伏电源出力预测数据<img file="FDA0000853511620000023.GIF" wi="137" he="94" />γ<sub>j</sub>;步骤2:将优化规划模型分解成3个相互制约的子问题,分别为无动态频率约束以及光伏调频容量约束的机组组合问题SP1、动态频率校验问题SP2、光伏调频备用容量校验问题SP3;步骤3:计算无动态频率约束以及光伏调频容量约束的机组组合问题SP1,问题SP1的目标函数为式一,约束条件为式二至式八,求解子问题SP1获得当前机组启停方案;步骤4:子问题SP2为动态频率最低点即约束式九的校验;将步骤3所得机组启停方案、以及出力代入检测问题SP2;如果SP2问题无越限情况,则进入下一步;如果检测到某个时段有动态频率最低点越限,且该时段系统发电机单位调节功率有增加的可能,则产生优化割,并在SP1中增加一个新的约束条件,然后计算SP1,重复步骤4和步骤5使SP1和SP2反复相互迭代,直至SP2问题满足;步骤5:将步骤3、步骤4迭代收敛最终所得机组启停方案、机组出力作为光伏调频容量检测问题SP3运算的初始值,计算SP3子问题的目标函数,若SP3子问题目标函数为0,则代表光伏调频容量校验子问题SP3无越限,所得机组启停方案、机组出力为系统最优方案;若越限,则进入下一步;步骤6:将步骤5中得到的非0惩罚变量值进行标记,并在SP1中增加一个新的约束条件,然后计算SP1,重复步骤4、5、6;使SP1、SP2、SP3反复相互迭代,直至得到系统最优方案。
地址 430072 湖北省武汉市武昌区珞珈山武汉大学