发明名称 一种配电网最小化开关操作次数的线性恢复控制方法
摘要 本发明涉及一种配电网最小化开关操作次数的线性恢复控制方法,属于电力系统调度自动化领域。该方法建立最大化恢复更多失电负荷和最小化开关操作次数的多目标优化模型,通过等效变换建立了线性多目标函数,通过选择开关状态、节点电压平方、支路有功、无功作为变量、引入大M数的方式建立了恢复控制的线性约束条件,从而将主流的高度非线性的混合整数二次模型转换为混合整数线性规划模型,本方法可以快速求解,解决了传统模型求解成功率低、可解性差、耗时长的问题,可以满足在线恢复控制的要求。
申请公布号 CN105160438A 申请公布日期 2015.12.16
申请号 CN201510638864.0 申请日期 2015.09.29
申请人 贵州电网公司信息通信分公司;北京清大高科系统控制有限公司 发明人 吴忠;王玮;初祥祥;纪元;王鹏;张克贤;汤磊;尹佳;杨耀
分类号 G06Q10/04(2012.01)I;G06Q50/06(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 廖元秋
主权项 一种配电网最小化开关操作次数的线性恢复控制方法,包括以下步骤:1)将配电网中包含的节点总个数、支路总个数、馈线根节点个数分别记为N<sub>node</sub>、N<sub>branch</sub>、N<sub>root</sub>,该支路包括线路及线路串联的开关;故障隔离后失电的负荷节点集合记为Φ<sub>out</sub>;网络中所有可用的节点集合记为Φ<sub>all</sub>;用j∈i代表节点j与节点i有连接关系;故障前节点i处的负荷有功、无功分别记为<img file="FDA0000816312700000011.GIF" wi="208" he="79" />支路i‑j的优化后开关状态记为x<sub>ij</sub>,当前即优化前的开关状态为x′<sub>ij</sub>,x<sub>ij</sub>和x′<sub>ij</sub>取值为0或者1,安全电流限值记为<img file="FDA0000816312700000012.GIF" wi="146" he="78" />是支路的铭牌参数,有功与无功分别记为P<sub>ij</sub>、Q<sub>ij</sub>,电阻与电抗值分别记为R<sub>ij</sub>、X<sub>ij</sub>;节点i的电压幅值及其上、下限分别记为u<sub>i</sub>、<img file="FDA0000816312700000013.GIF" wi="239" he="69" />和<img file="FDA0000816312700000014.GIF" wi="46" he="63" />是事先指定的常数,节点i电压的平方记为U<sub>i</sub>,M<sub>0</sub>是一个较大的常数,取10000以上;2)建立如式(1)所示目标函数,该目标函数包括两个子目标,第一个子目标为尽可能多的恢复失电负荷,第二个子目标为开关操作次数尽可能少:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>minw</mi><mn>1</mn></msub><munder><mo>&Sigma;</mo><mrow><mi>i</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow></msub></mrow></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><mo>&Sigma;</mo><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>P</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>w</mi><mn>2</mn></msub><munder><mo>&Sigma;</mo><mrow><mi>i</mi><mi>j</mi><mo>&Element;</mo><msub><mi>N</mi><mrow><mi>b</mi><mi>r</mi><mi>a</mi><mi>n</mi><mi>c</mi><mi>h</mi></mrow></msub></mrow></munder><mo>|</mo><mrow><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>-</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow><mo>&prime;</mo></msubsup></mrow><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000816312700000015.GIF" wi="1391" he="117" /></maths>式中w<sub>1</sub>为第一个子目标的权重因子,w<sub>2</sub>为第二个子目标的权重因子,其中w<sub>1</sub>+w<sub>2</sub>=1,w<sub>1</sub>、w<sub>2</sub>均为非负实数;3)开关状态x′<sub>ij</sub>按照0为断开,1为闭合定义,将式(1)的目标函数等效转化成式(2)的优化线性目标函数:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>minw</mi><mn>1</mn></msub><munder><mo>&Sigma;</mo><mrow><mi>i</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow></msub></mrow></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><mo>&Sigma;</mo><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>P</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>w</mi><mn>2</mn></msub><mrow><mo>(</mo><munder><mo>&Sigma;</mo><mrow><mi>i</mi><mi>j</mi><mo>&Element;</mo><msub><mi>N</mi><mrow><mi>b</mi><mi>r</mi><mi>a</mi><mi>n</mi><mi>c</mi><mi>h</mi></mrow></msub><mo>,</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow><mo>&prime;</mo></msubsup><mo>=</mo><mn>1</mn></mrow></munder><mn>1</mn><mo>-</mo><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>+</mo><munder><mo>&Sigma;</mo><mrow><mi>i</mi><mi>j</mi><mo>&Element;</mo><msub><mi>N</mi><mrow><mi>b</mi><mi>r</mi><mi>a</mi><mi>n</mi><mi>c</mi><mi>h</mi></mrow></msub><mo>,</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow><mo>&prime;</mo></msubsup><mo>=</mo><mn>0</mn></mrow></munder><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000816312700000016.GIF" wi="1542" he="118" /></maths>4)建立的如式(3)~(16)所示线性约束条件:x<sub>ij</sub>∈{0,1},i,j∈Φ<sub>all</sub>        (3)<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><munder><mo>&Sigma;</mo><mrow><mi>i</mi><mo>&NotEqual;</mo><mi>j</mi></mrow></munder><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><msub><mi>N</mi><mrow><mi>n</mi><mi>o</mi><mi>d</mi><mi>e</mi></mrow></msub><mo>-</mo><msub><mi>N</mi><mrow><mi>r</mi><mi>o</mi><mi>o</mi><mi>t</mi></mrow></msub><mo>,</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mrow><mi>a</mi><mi>l</mi><mi>l</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000816312700000017.GIF" wi="1254" he="110" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><munder><mo>&Sigma;</mo><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>P</mi><mrow><mi>i</mi><mi>j</mi></mrow></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><mrow><mi>a</mi><mi>l</mi><mi>l</mi></mrow></msub><mo>-</mo><msub><mi>&Phi;</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000816312700000018.GIF" wi="1231" he="116" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><munder><mo>&Sigma;</mo><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>P</mi><mrow><mi>i</mi><mi>j</mi></mrow></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><mrow><mi>a</mi><mi>l</mi><mi>l</mi></mrow></msub><mo>-</mo><msub><mi>&Phi;</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000816312700000019.GIF" wi="1226" he="112" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><munder><mo>&Sigma;</mo><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>Q</mi><mrow><mi>i</mi><mi>j</mi></mrow></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><mrow><mi>a</mi><mi>l</mi><mi>l</mi></mrow></msub><mo>-</mo><msub><mi>&Phi;</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00008163127000000110.GIF" wi="1228" he="106" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><munder><mo>&Sigma;</mo><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>P</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00008163127000000111.GIF" wi="1158" he="107" /></maths><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msubsup><mi>L</mi><mrow><mi>p</mi><mo>,</mo><mi>i</mi></mrow><mn>0</mn></msubsup><mo>-</mo><munder><mo>&Sigma;</mo><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>P</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000816312700000021.GIF" wi="1211" he="123" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msubsup><mi>L</mi><mrow><mi>p</mi><mo>,</mo><mi>i</mi></mrow><mn>0</mn></msubsup><mo>-</mo><munder><mo>&Sigma;</mo><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>P</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>&Phi;</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000816312700000022.GIF" wi="1212" he="131" /></maths><maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mo>-</mo><mn>0.7</mn><mo>*</mo><msub><mover><mi>I</mi><mo>&OverBar;</mo></mover><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&le;</mo><msub><mi>P</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&le;</mo><mn>0.7</mn><mo>*</mo><msub><mover><mi>I</mi><mo>&OverBar;</mo></mover><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000816312700000023.GIF" wi="1186" he="99" /></maths><maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><mo>-</mo><mn>0.7</mn><mo>*</mo><msub><mover><mi>I</mi><mo>&OverBar;</mo></mover><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&le;</mo><msub><mi>Q</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&le;</mo><mn>0.7</mn><mo>*</mo><msub><mover><mi>I</mi><mo>&OverBar;</mo></mover><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000816312700000024.GIF" wi="1194" he="101" /></maths>M<sub>ij</sub>=(1‑x<sub>ij</sub>)M<sub>0</sub>             (13)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>Xi<sub>j)</sub>,i,j∈Φ<sub>all</sub>       (14)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>),i,j∈Φ<sub>all</sub>       (15)<maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><msubsup><mover><mi>u</mi><mo>&OverBar;</mo></mover><mi>i</mi><mn>2</mn></msubsup><mo>&GreaterEqual;</mo><msubsup><mi>U</mi><mi>i</mi><mn>2</mn></msubsup><mo>&GreaterEqual;</mo><msup><msub><munder><mi>u</mi><mo>&OverBar;</mo></munder><mi>i</mi></msub><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000816312700000025.GIF" wi="1098" he="92" /></maths>5)求解由式(2)代表的优化线性目标函数及式(3)~(16)线性约束条件形成的模型,该模型为混合整数线性规划,求解得到各个支路开关的优化状态,判断支路i‑j的优化开关状态x<sub>ij</sub>,若不等于优化前的开关x′<sub>ij</sub>,则进行开关变位操作;遍历所有支路找到优化前后开关状态不同的支路,其开关变位操作后的状态构成恢复负荷的开关序列。
地址 550003 贵州省贵阳市瑞金南路38号