发明名称 基于配电网故障失负荷恢复值的配电网运行安全评估方法
摘要 本发明涉及一种基于配电网故障失负荷恢复值的运行安全评估方法,属于电力系统调度自动化与电网仿真技术领域。该方法对配电网进行N-1故障扫描,每一故障下对失电负荷进行恢复优化计算得到能够恢复的失电负荷值;扫描完成后,累加每次恢复的失电负荷得到总的恢复失电负荷值。根据总失电负荷值、总的可能恢复的失负荷值、总的恢复负荷值,建立配电网运行安全指标。这其中,故障失负荷恢复优化是考虑了有功、无功、电压及辐射状运行约束的配电网详细恢复优化模型,该混合整数二次约束规划问题采用分支定界法求解;本发明基于上述优化模型,制定的配电网运行安全评估的指标,能够反映出当前配电网运行条件和网络结构的安全性。
申请公布号 CN103366315A 申请公布日期 2013.10.23
申请号 CN201310314232.X 申请日期 2013.07.24
申请人 国家电网公司;清华大学;辽宁省电力有限公司沈阳供电公司 发明人 吴文传;张伯明;郭昆亚;孙宏斌;陈珂宁;王英男;黄哲洙
分类号 G06Q50/06(2012.01)I 主分类号 G06Q50/06(2012.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 廖元秋
主权项 1.一种基于配电网故障失负荷恢复值的配电网运行安全评估方法,其特征在于,该方法包括以下步骤:1)构建带混合整数二次约束的优化模型:1.1)构建初始恢复优化模型如下所示:<maths num="0001"><![CDATA[<math><mrow><mi>Obj</mi><mo>.</mo><mi>Min</mi><munder><mi>&Sigma;</mi><msub><mrow><mi>i</mi><mo>&Element;</mo><mi>&Phi;</mi></mrow><mi>out</mi></msub></munder><mrow><mo>(</mo><msubsup><mi>L</mi><mrow><mi>p</mi><mo>,</mo><mi>i</mi></mrow><mn>0</mn></msubsup><mo>-</mo><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>x</mi><mi>ij</mi></msub><msub><mi>P</mi><mi>ij</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>xij∈{0,1},i,j∈Φall(2)<maths num="0002"><![CDATA[<math><mrow><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>x</mi><mi>ij</mi></msub><msub><mi>P</mi><mi>ij</mi></msub><mo>=</mo><msubsup><mi>L</mi><mrow><mi>p</mi><mo>,</mo><mi>i</mi></mrow><mn>0</mn></msubsup><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mi>all</mi></msub><mo>-</mo><msub><mi>&Phi;</mi><mi>loss</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>x</mi><mi>ij</mi></msub><msub><mrow><mo>.</mo><mi>Q</mi></mrow><mi>ij</mi></msub><mo>=</mo><msubsup><mi>L</mi><mrow><mi>p</mi><mo>,</mo><mi>i</mi></mrow><mn>0</mn></msubsup><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mi>all</mi></msub><mo>-</mo><msub><mi>&Phi;</mi><mi>loss</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><munder><mi>&Sigma;</mi><mrow><mi>i</mi><mo>&NotEqual;</mo><mi>j</mi></mrow></munder><msub><mi>x</mi><mi>ij</mi></msub><mo>=</mo><msub><mi>N</mi><mi>node</mi></msub><mo>-</mo><msub><mi>N</mi><mi>root</mi></msub><mo>,</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mi>all</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>x</mi><mi>ij</mi></msub><msub><mi>P</mi><mi>ij</mi></msub><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mi>loss</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msubsup><mi>L</mi><mrow><mi>p</mi><mo>,</mo><mi>i</mi></mrow><mn>0</mn></msubsup><mo>-</mo><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>x</mi><mi>ij</mi></msub><msub><mi>P</mi><mi>ij</mi></msub><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mi>loss</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><mfrac><mrow><msubsup><mi>P</mi><mi>ij</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>Q</mi><mi>ij</mi><mn>2</mn></msubsup></mrow><msubsup><mi>V</mi><mi>i</mi><mn>2</mn></msubsup></mfrac><mo>&le;</mo><msubsup><mover><mi>I</mi><mo>&OverBar;</mo></mover><mi>ij</mi><mn>2</mn></msubsup><mo>,</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mi>all</mi></msub><mo>,</mo><mi>i</mi><mo>&NotEqual;</mo><mi>j</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msubsup><mi>V</mi><mi>i</mi><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>V</mi><mi>j</mi><mn>2</mn></msubsup><mo>=</mo><mn>2</mn><mrow><mo>(</mo><msub><mi>R</mi><mi>ij</mi></msub><msub><mi>P</mi><mi>ij</mi></msub><mo>+</mo><msub><mi>Q</mi><mi>ij</mi></msub><msub><mi>X</mi><mi>ij</mi></msub><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><msubsup><mi>R</mi><mi>ij</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>X</mi><mi>ij</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mfrac><mrow><msubsup><mi>P</mi><mi>ij</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>Q</mi><mi>ij</mi><mn>2</mn></msubsup></mrow><msubsup><mi>V</mi><mi>i</mi><mn>2</mn></msubsup></mfrac><mo>,</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mrow><mi>all</mi><mo>,</mo><mi>i</mi><mo>&NotEqual;</mo></mrow></msub><mi>j</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0009"><![CDATA[<math><mrow><munder><msub><mi>V</mi><mi>i</mi></msub><mo>&OverBar;</mo></munder><mo>&le;</mo><msub><mi>V</mi><mi>i</mi></msub><mo>&le;</mo><msub><mover><mi>V</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mi>all</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,N<sub>node</sub>为节点总个数;N<sub>branch</sub>为支路总个数;N<sub>root</sub>为馈线根节点个数;Φ<sub>loss</sub>为故障隔离后失电的负荷节点集合;Φ<sub>all</sub>为网络中所有节点集合;j∈i为表明与节点i有连接关系的节点;<img file="FDA000035607668000113.GIF" wi="72" he="66" />和<img file="FDA000035607668000114.GIF" wi="70" he="69" />为故障前节点i处的负荷有功和无功数值;x<sub>ij</sub>为模型中支路的连接状态,0和1分别表示断开及闭合;P<sub>ij</sub>和Q<sub>ij</sub>为支路ij上的有功与无功潮流数值;<img file="FDA000035607668000110.GIF" wi="54" he="93" />为支路ij的安全电流限值;R<sub>ij</sub>和X<sub>ij</sub>为支路ij的电阻与电抗值;V<sub>i</sub>,<img file="FDA000035607668000111.GIF" wi="80" he="74" /><img file="FDA000035607668000112.GIF" wi="51" he="81" />分别为节点i的电压幅值、电压上限、电压下限;该初始优化模型中,式(1)为目标函数,要求失负荷尽可能地少;式(2)-(10)为约束;式(2)为反映支路开闭情况的01变量;(3)-(4)为功率平衡约束;(3-5)共同组成了辐射状约束;(6)-(7)为失电节点恢复功率约束;(8)为线路热稳定约束;(9)-(10)为支路两端电压关系及节点电压约束;1.2)对初始优化模型中的部分约束形式进行处理:a)将约束(8)改进为:<maths num="0010"><![CDATA[<math><mrow><msubsup><mi>P</mi><mi>ij</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>Q</mi><mi>ij</mi><mn>2</mn></msubsup><mo>&le;</mo><msub><mi>x</mi><mi>ij</mi></msub><mo>.</mo><msubsup><mover><mi>I</mi><mo>&OverBar;</mo></mover><mi>ij</mi><mn>2</mn></msubsup><mo>.</mo><msubsup><mi>V</mi><mn>0</mn><mn>2</mn></msubsup><mo>,</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mi>all</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>'</mo><mo>)</mo></mrow></mrow></math>]]></maths>其中的<img file="FDA00003560766800022.GIF" wi="146" he="87" />为常数,<img file="FDA00003560766800023.GIF" wi="62" he="87" />为支路电流上限,V<sub>0</sub>为节点预设电压,都是常数;式(8)改进后,式(3)-(4)和(6)-(7)中不同变量相乘的形式:x<sub>ij</sub>·P<sub>ij</sub>与x<sub>ij</sub>·Q<sub>ij</sub>也被分别表示为P<sub>ij</sub>和Q<sub>ij</sub>:<maths num="0011"><![CDATA[<math><mrow><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>P</mi><mi>ij</mi></msub><mo>=</mo><msubsup><mi>L</mi><mrow><mi>p</mi><mo>,</mo><mi>i</mi></mrow><mn>0</mn></msubsup><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mi>all</mi></msub><mo>-</mo><msub><mi>&Phi;</mi><mi>loss</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>'</mo><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0012"><![CDATA[<math><mrow><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>Q</mi><mi>ij</mi></msub><mo>=</mo><msubsup><mi>L</mi><mrow><mi>q</mi><mo>,</mo><mi>i</mi></mrow><mn>0</mn></msubsup><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mi>all</mi></msub><mo>-</mo><msub><mi>&Phi;</mi><mi>loss</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>'</mo><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0013"><![CDATA[<math><mrow><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>P</mi><mi>ij</mi></msub><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mi>loss</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>'</mo><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0014"><![CDATA[<math><mrow><msubsup><mi>L</mi><mrow><mi>p</mi><mo>,</mo><mi>i</mi></mrow><mn>0</mn></msubsup><mo>-</mo><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>P</mi><mi>ij</mi></msub><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mi>loss</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>'</mo><mo>)</mo></mrow></mrow></math>]]></maths>使所述优化模型变得可求解;b)约束(9)改为:U<sub>i</sub>-U<sub>j</sub>=2(P<sub>ij</sub>R<sub>ij</sub>+Q<sub>ij</sub>X<sub>ij</sub>)(9’)i,j∈Φ<sub>all</sub>,U<sub>i</sub>≥0,U<sub>j</sub>≥0其中<img file="FDA00003560766800028.GIF" wi="180" he="78" />U<sub>i</sub>≥0、<img file="FDA00003560766800029.GIF" wi="182" he="85" />U<sub>j</sub>≥0;c)对支路两端电压关系的约束(9’)进一步处理:引入变量Mij:M<sub>ij</sub>=(1-x<sub>ij</sub>)·M<sub>0</sub>(11)其中的M<sub>0</sub>是给定的一个足够大的正数,进而(9’)改进为:U<sub>i</sub>-U<sub>j</sub>≤M<sub>ij</sub>+2(P<sub>ij</sub>R<sub>ij</sub>+Q<sub>ij</sub>X<sub>ij</sub>)U<sub>i</sub>-U<sub>j</sub>≥-M<sub>ij</sub>+2(P<sub>ij</sub>R<sub>ij</sub>+Q<sub>ij</sub>X<sub>ij</sub>);(9”)i,j∈Φ<sub>all</sub>,U<sub>i</sub>≥0,U<sub>j</sub>≥0d.)对约束(3)‐(5)进行处理,以解决零注入节点产生的问题:给所有零注入节点增设一个足够小的负荷ε,即对于没有失电的节点,如果其初始负荷为0,则设<img file="FDA00003560766800031.GIF" wi="214" he="86" />而对于失电节点,约束(6)改为:<maths num="0015"><![CDATA[<math><mrow><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>P</mi><mi>ij</mi></msub><mo>&GreaterEqual;</mo><mi>&epsiv;</mi><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mi>loss</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>'</mo><mo>)</mo></mrow></mrow></math>]]></maths>通过步骤1.2)改进后原初始优化模型成为带线性目标、带混合整数二次约束的优化模型;该优化模型以(1)为优化目标,约束为(2)(3’)(4’)(5)(6’)(7’)(8’)(9’’)(11);2)基于配电网故障失负荷恢复值对配电网的运行安全进行评估:具体步骤如下:2.1)网络的初始化:将配电网简化为可开闭的线路,和带有负荷的节点;2.2)令i=1;2.3)设置第i条线路为故障线路;2.4)故障隔离:将故障线路i隔离,故障隔离后,配电网网络分为两个部分,一个部分是带电部分,另一个部分是包括一个或多个失电区域的失电部分;失电区域进一步分为可能恢复的失电区域和不可能恢复的失电区域;由带电部分和失电部分的可能恢复的失电区域的线路和节点获得用于计算的网络结构、总的失电负荷数值和可能恢复的失电负荷数值;2.5)求解失负荷恢复带混合整数二次约束的优化模型将步骤2.1)-2.4)中得到的线路和节点网络结构、节点负荷值、线路电流限额和节点电压限额作为带混合整数二次约束的优化模型的输入,该网络结构仅包括带电部分和可能恢复的失电部分;用分支定界法求解带混合整数二次约束的优化模型,得到电网故障失负荷恢复值;即得到当前第i条线路故障情况下的实际恢复失电负荷值;2.6)“N-1”故障扫描判断:若i=N,N为线路总个数,则故障扫描已完成,进行步骤2.7);若i&lt;N,故障扫描未完成,即还有其它线路需要计算,令i=i+1,回到步骤2.3);2.7)计算配电网安全评估结果:定义如下配电网运行安全指标作为配电网安全评估结果:<maths num="0016"><![CDATA[<math><mrow><mi>Index</mi><mn>2</mn><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>L</mi><mrow><mi>restorable</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>L</mi><mrow><mi>outage</mi><mo>,</mo><mi>i</mi></mrow></msub></mrow></math>]]></maths><maths num="0017"><![CDATA[<math><mrow><mi>Index</mi><mn>2</mn><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>L</mi><mrow><mi>restorable</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>L</mi><mrow><mi>outage</mi><mo>,</mo><mi>i</mi></mrow></msub></mrow></math>]]></maths>其中,<img file="FDA00003560766800034.GIF" wi="206" he="142" />为故障扫描中各次总失电负荷之和;<img file="FDA00003560766800035.GIF" wi="245" he="141" />为故障扫描中各次可能恢复的失电负荷之和;<img file="FDA00003560766800036.GIF" wi="225" he="142" />为由恢复算法计算得到的各次实际能够恢复的失电负荷之和;由Index1和Index2共同表征配电网的运行安全性,是最终的安全评估结果。
地址 100031 北京市西城区西长安街86号