发明名称 微电网运行方法
摘要 本发明公开了一种微电网运行方法,该微电网中的电源包括风电、火电、小水电群,该方法包括以下步骤:对微电网建立经济运行策略问题目标函数模型;利用模拟退火算法和粒子群算法分别求解目标函数模型的值,比较所求出的两个值并协同优化,计算出目标函数模型的最优值,确立微电网的运行模式;按照上一步所求出的最优值确立微电网的运行模式,运行电网。本发明应用对象为含有风电、火电、小水电群作为电源供给的微电网,设计了一种电力系统经济调度优化运行算法,着重分析了小水电群应采用的运行策略,应用于求解微网独立运行状态下的经济运行模式,实现该类型微电网的平稳可靠、高效益运行。
申请公布号 CN103580061B 申请公布日期 2015.05.20
申请号 CN201310516199.9 申请日期 2013.10.28
申请人 贵州电网公司电网规划研究中心;清华大学 发明人 郑竞宏;沈欣炜;朱守真;李金霞
分类号 H02J3/38(2006.01)I;G06Q10/04(2012.01)I;G06Q50/06(2012.01)I 主分类号 H02J3/38(2006.01)I
代理机构 北京卓言知识产权代理事务所(普通合伙) 11365 代理人 王茀智;龚清媛
主权项 一种微电网运行方法,其特征在于,所述微电网中的电源包括风电、火电、小水电群,该方法包括以下步骤:A.对微电网建立经济运行策略问题目标函数模型;B.利用模拟退火算法和粒子群算法分别求解目标函数模型的值,比较所求出的两个值并协同优化,计算出目标函数模型的最优值,确立微电网的运行模式;和C.按照步骤B所求出的最优值确立微电网的运行模式,运行电网;将所述微电网中的电源划分为:出力可调、作为电压频率参考源的电源S<sub>VF</sub>;出力可调但不作为S<sub>VF</sub>电源的电源S<sub>AG</sub>;和出力不可调的电源S<sub>UG</sub>;所述经济运行策略问题目标函数模型中:目标函数为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>obj</mi><mi>min</mi><mo>&CenterDot;</mo><munder><mi>&Sigma;</mi><mi>t</mi></munder><mo>[</mo><munder><mi>&Sigma;</mi><mrow><mi>i</mi><mo>&Element;</mo><msub><mi>S</mi><mi>VF</mi></msub></mrow></munder><mrow><mo>(</mo><msub><mi>a</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><msubsup><mi>P</mi><mrow><mi>VFi</mi><mo>,</mo><mi>t</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msub><mi>a</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><msub><mi>P</mi><mrow><mi>VFi</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>+</mo><msub><mi>a</mi><mrow><mn>0</mn><mi>i</mi></mrow></msub><mo>)</mo></mrow><mo>]</mo><mo>,</mo></mrow>]]></math><img file="FDA0000677593190000011.GIF" wi="795" he="120" /></maths>其中P<sub>VFi,t</sub>为控制变量,表示S<sub>VF</sub>电源在第t时段的出力,a<sub>2i</sub>、a<sub>1i</sub>、a<sub>0i</sub>为已知参数,表示S<sub>VF</sub>电源耗量特性曲线参数,约束条件包括:节点有功潮流方程和节点无功潮流方程<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>VFi</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>Di</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>-</mo><msub><mi>P</mi><mrow><mi>AGi</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>-</mo><msub><mi>U</mi><mi>UGi</mi></msub><msub><mi>P</mi><mrow><mi>UGi</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>+</mo><msub><mi>V</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>V</mi><mrow><mi>j</mi><mo>,</mo><mi>t</mi></mrow></msub><mrow><mo>(</mo><msub><mi>G</mi><mrow><mi>ij</mi><mo>,</mo><mi>t</mi></mrow></msub><mi>cos</mi><msub><mi>&theta;</mi><mrow><mi>ij</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>+</mo><msub><mi>B</mi><mrow><mi>ij</mi><mo>,</mo><mi>t</mi></mrow></msub><mi>sin</mi><msub><mi>&theta;</mi><mrow><mi>ij</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000677593190000012.GIF" wi="1432" he="145" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>Q</mi><mrow><mi>Ri</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>=</mo><msub><mi>Q</mi><mrow><mi>Di</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>+</mo><msub><mi>V</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>V</mi><mrow><mi>j</mi><mo>,</mo><mi>t</mi></mrow></msub><mrow><mo>(</mo><msub><mi>G</mi><mrow><mi>ij</mi><mo>,</mo><mi>t</mi></mrow></msub><mi>sin</mi><msub><mi>&theta;</mi><mrow><mi>ij</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>-</mo><msub><mi>B</mi><mrow><mi>ij</mi><mo>,</mo><mi>t</mi></mrow></msub><mi>cos</mi><msub><mi>&theta;</mi><mrow><mi>ij</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000677593190000013.GIF" wi="1036" he="148" /></maths>有功出力约束和无功出力约束<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><munder><mi>P</mi><mo>&OverBar;</mo></munder><mi>VFi</mi></msub><mo>&le;</mo><msub><mi>P</mi><mrow><mi>VFi</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&le;</mo><msub><mover><mi>P</mi><mo>&OverBar;</mo></mover><mi>VFi</mi></msub><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>S</mi><mi>VF</mi></msub></mrow>]]></math><img file="FDA0000677593190000014.GIF" wi="561" he="95" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><munder><mi>Q</mi><mo>&OverBar;</mo></munder><mi>Ri</mi></msub><mo>&le;</mo><msub><mi>Q</mi><mrow><mi>Ri</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&le;</mo><msub><mover><mi>Q</mi><mo>&OverBar;</mo></mover><mi>Ri</mi></msub><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>S</mi><mi>R</mi></msub></mrow>]]></math><img file="FDA0000677593190000015.GIF" wi="514" he="95" /></maths>机组爬坡约束‑r<sub>Gi,dn</sub>T≤P<sub>VFi,t</sub>‑P<sub>VFi,t‑1</sub>≤r<sub>Gi,up</sub>T,i∈S<sub>VF</sub>节点电压约束<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><munder><msub><mi>V</mi><mi>i</mi></msub><mo>&OverBar;</mo></munder><mo>&le;</mo><msub><mi>V</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&le;</mo><msub><mover><mi>V</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi></mrow>]]></math><img file="FDA0000677593190000016.GIF" wi="480" he="103" /></maths>线路传输约束<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><munder><mi>P</mi><mo>&OverBar;</mo></munder><mi>ij</mi></msub><mo>&le;</mo><msub><mi>P</mi><mrow><mi>ij</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&le;</mo><msub><mover><mi>P</mi><mo>&OverBar;</mo></mover><mi>ij</mi></msub><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&Element;</mo><msub><mi>S</mi><mi>Cb</mi></msub></mrow>]]></math><img file="FDA0000677593190000021.GIF" wi="605" he="100" /></maths>第t时段线路(i,j)传输功率P<sub>ij,t</sub>为<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>ij</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>=</mo><msub><mi>V</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>V</mi><mrow><mi>j</mi><mo>,</mo><mi>t</mi></mrow></msub><mrow><mo>(</mo><msub><mi>G</mi><mrow><mi>ij</mi><mo>,</mo><mi>t</mi></mrow></msub><mi>cos</mi><msub><mi>&theta;</mi><mrow><mi>ij</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>+</mo><msub><mi>B</mi><mrow><mi>ij</mi><mo>,</mo><mi>t</mi></mrow></msub><mi>sin</mi><msub><mi>&theta;</mi><mrow><mi>ij</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000677593190000022.GIF" wi="826" he="144" /></maths>其中,Q<sub>Ri,t</sub>,U<sub>UGi</sub>为控制变量,具体地:Q<sub>Ri,t</sub>表示第t时段节点i可调无功电源出力,U<sub>UGi</sub>为并网标识符,取值1时表示并网发电,取值为0时表示不并网发电,其中,P<sub>UGi,t</sub>、P<sub>Di,t</sub>、Q<sub>Di,t</sub>、G<sub>ij,t</sub>、B<sub>ij,t</sub>、<u>P</u><sub>VFi</sub>、<img file="FDA0000677593190000023.GIF" wi="113" he="76" />Q<sub>Ri</sub>、<img file="FDA0000677593190000024.GIF" wi="109" he="82" />r<sub>Gi,dn</sub>、r<sub>Gi,up</sub>、<u>V<sub>i</sub></u>、<img file="FDA0000677593190000025.GIF" wi="80" he="83" /><u>P</u><sub>ij</sub>、<img file="FDA0000677593190000026.GIF" wi="89" he="73" />参为已知参数,具体地:P<sub>UGi,t</sub>表示第t时段节点i上属于UG电源的有功出力,P<sub>Di,t</sub>和Q<sub>Di,t</sub>分别表示第t时段节点i的有功无功负荷,G<sub>ij,t</sub>、B<sub>ij,t</sub>分别表示第t时段线路(i,j)的电导、电纳,<u>P</u><sub>VFi</sub>和<img file="FDA0000677593190000027.GIF" wi="86" he="68" />表示节点i上S<sub>VF</sub>电源有功出力的约束上下限,<u>Q</u><sub>Ri</sub>和<img file="FDA0000677593190000028.GIF" wi="76" he="77" />分别表示节点i可调无功电源出力的约束上下限,r<sub>Gi,dn</sub>和r<sub>Gi,up</sub>分别表示机组的上、下爬坡速率,<u>V<sub>i</sub></u>和<img file="FDA0000677593190000029.GIF" wi="54" he="82" />分别表示第t时段节点i的电压的约束上下限,<u>P</u><sub>ij</sub>和<img file="FDA00006775931900000210.GIF" wi="64" he="73" />分别表示第t时段线路(i,j)传输功率的约束上下限,其中,P<sub>AGi,t</sub>、V<sub>i,t</sub>、θ<sub>ij,t</sub>、P<sub>ij,t</sub>为依从变量,具体地,V<sub>i,t</sub>表示第t时段节点i的电压,θ<sub>ij,t</sub>表示第t时段节点i、j之间相角差,P<sub>ij,t</sub>表示第t时段线路(i,j)传输的有功功率,其中,P<sub>AGi,t</sub>表示第t时间段节点i上的库容式小水电有功出力,其中,S<sub>VF</sub>为可调有功出力电源的集合、S<sub>R</sub>为无功电源的集合,n为仿真系统节点个数,S<sub>Cb</sub>是约束线路的集合。
地址 550003 贵州省贵阳市瑞金南路38号