发明名称 基于网损等值负荷模型的改进直流最优潮流方法
摘要 本发明公开了一种基于网损等值负荷模型的改进直流最优潮流(MDCOPF)。虽然传统直流最优潮流(DCOPF)有着很高的计算效率,但是其模型中的有功平衡方程忽略了支路电阻的影响,计算结果精度偏低。本发明针对DCOPF的不足之处,通过网损迭代的方式得到等值负荷的大小,并建立网损等值负荷模型,提出了一种基于网损等值负荷模型的MDCOPF。经过多个算例的仿真测试,结果表明与传统DCOPF相比,本发明所提方法有着更高的计算精度和计算效率,并且对于病态电网具有较强的处理能力,证明了本发明的正确性和有效性。
申请公布号 CN104794531A 申请公布日期 2015.07.22
申请号 CN201510153568.1 申请日期 2015.04.01
申请人 河海大学 发明人 卫志农;何天雨;孙国强;孙永辉;高沁;臧海祥;朱瑛;向育鹏;李海欣
分类号 G06Q10/04(2012.01)I;G06Q50/06(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 李玉平
主权项 一种基于网损等值负荷模型的改进直流最优潮流方法,其特征在于,所述方法依次按以下步骤实现:(1)获得电力系统的网络参数信息;(2)程序初始化;包括:对算法中的状态变量设置初值、设置迭代计数器k=0、最大迭代次数K<sub>max</sub>、收敛精度ε、形成节点导纳矩阵B、设置网损等值负荷初值P<sub>equ</sub>=0;(3)计算互补间隙Gap,判断其是否满足精度要求,若满足,则输出最优解,结束循环,否则,继续;(4)计算等式约束、不等式约束的雅可比矩阵▽<sub>x</sub>h(x)、<img file="FDA0000692637300000011.GIF" wi="223" he="94" />计算目标函数、等式约束、不等式约束的海森矩阵<img file="FDA0000692637300000012.GIF" wi="713" he="98" />以及各常数项L′<sub>x</sub>、L<sub>y</sub>、L<sub>w</sub>、<img file="FDA0000692637300000013.GIF" wi="99" he="78" />并根据以下方程求解总变量x、拉格朗日乘子y和松弛变量u、w的增量Δx、Δy、Δu、Δw;(5)确定原始变量和对偶变量的迭代步长:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&alpha;</mi><mi>p</mi></msub><mo>=</mo><mn>0.9995</mn><mi>min</mi><mo>{</mo><munder><mi>min</mi><mi>i</mi></munder><mrow><mo>(</mo><mfrac><mrow><mo>-</mo><msub><mi>u</mi><mi>i</mi></msub></mrow><mrow><mi>&Delta;</mi><msub><mi>u</mi><mi>i</mi></msub></mrow></mfrac><msub><mo>|</mo><mrow><mi>&Delta;</mi><msub><mi>u</mi><mi>i</mi></msub><mo>&lt;</mo><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><mn>1</mn><mo>}</mo></mrow>]]></math><img file="FDA0000692637300000014.GIF" wi="830" he="220" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&alpha;</mi><mi>d</mi></msub><mo>=</mo><mn>0.9995</mn><mi>min</mi><mo>{</mo><munder><mi>min</mi><mi>i</mi></munder><mrow><mo>(</mo><mfrac><mrow><mo>-</mo><msub><mi>w</mi><mi>i</mi></msub></mrow><mrow><mi>&Delta;</mi><msub><mi>w</mi><mi>i</mi></msub></mrow></mfrac><msub><mo>|</mo><mrow><mi>&Delta;</mi><msub><mi>w</mi><mi>i</mi></msub><mo>></mo><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><mn>1</mn><mo>}</mo></mrow>]]></math><img file="FDA0000692637300000015.GIF" wi="835" he="219" /></maths>(6)按照下式更新所有变量和拉格朗日乘子:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msup></mtd></mtr><mtr><mtd><msup><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msup></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup></mtd></mtr><mtr><mtd><msup><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup></mtd></mtr></mtable></mfenced><mo>+</mo><msub><mi>&alpha;</mi><mi>p</mi></msub><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Delta;x</mi></mtd></mtr><mtr><mtd><mi>&Delta;u</mi></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000692637300000016.GIF" wi="605" he="198" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msup></mtd></mtr><mtr><mtd><msup><mi>w</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msup></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup></mtd></mtr><mtr><mtd><msup><mi>w</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup></mtd></mtr></mtable></mfenced><mo>+</mo><msub><mi>&alpha;</mi><mi>d</mi></msub><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Delta;y</mi></mtd></mtr><mtr><mtd><mi>&Delta;w</mi></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000692637300000017.GIF" wi="627" he="198" /></maths>(7)提取总变量x中的节点电压相角的信息,储存在θ中,并根据以下公式更新网损等值负荷P<sub>equ</sub>和改进直流最优潮流中的等式约束方程h(x):<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>P</mi><mi>ij</mi><mo>&prime;</mo></msubsup><mo>=</mo><mfrac><mrow><msub><mi>&theta;</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&theta;</mi><mi>j</mi></msub></mrow><msub><mi>x</mi><mi>ij</mi></msub></mfrac></mrow>]]></math><img file="FDA0000692637300000018.GIF" wi="278" he="173" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>loss</mi><mo>,</mo><mi>ij</mi></mrow></msub><mo>=</mo><msubsup><mi>I</mi><mrow><mi>r</mi><mo>,</mo><mi>ij</mi></mrow><mn>2</mn></msubsup><msub><mi>r</mi><mi>ij</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>S</mi><mi>ij</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mn>2</mn></msup><msub><mi>r</mi><mi>ij</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><mi>&alpha;</mi><msubsup><mi>P</mi><mi>ij</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mn>2</mn></msup><msub><mi>r</mi><mi>ij</mi></msub></mrow>]]></math><img file="FDA0000692637300000021.GIF" wi="782" he="111" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>equ</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>=</mo><munder><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><mrow><mi>j</mi><mo>&NotEqual;</mo><mi>i</mi></mrow></munder><mfrac><msub><mi>P</mi><mrow><mi>loss</mi><mo>,</mo><mi>ij</mi></mrow></msub><mn>2</mn></mfrac></mrow>]]></math><img file="FDA0000692637300000022.GIF" wi="365" he="195" /></maths>h(x)=ΔP=P<sub>G</sub>‑P<sub>D</sub>+B·θ‑P<sub>equ</sub>其中:θ<sub>i</sub>、θ<sub>j</sub>为节点i、j的相角;x<sub>ij</sub>、r<sub>ij</sub>为线路ij的电抗、电阻;P′<sub>ij</sub>、S′<sub>ij</sub>、P<sub>loss,ij</sub>为线路ij的有功功率、视在功率、有功功率损耗;I<sub>r,ij</sub>为流过电阻r<sub>ij</sub>的电流大小;α是线路视在功率与有功功率的比例因子;P<sub>equ,i</sub>为节点i的网损等值负荷;P<sub>G</sub>、P<sub>D</sub>为节点注入有功功率、有功负荷向量;h(x)为等式约束。(8)判断迭代次数是否小于最大迭代次数K<sub>max</sub>,若是,则令迭代次数加1,返回(3),否则,输出“计算不收敛”,结束程序。
地址 211100 江苏省南京市江宁区佛城西路8号