发明名称 一种二次规划形式的电力系统抗差状态估计方法
摘要 本发明提出一种二次规划形式的电力系统抗差状态估计方法,包括以下步骤:S1,提出二次规划形式的电力系统抗差状态估计基本模型;S2,对二次规划形式的电力系统抗差状态估计基本模型进行转化,得到其对偶模型,对偶模型同样具有二次规划形式;以及S3,对对偶模型求解。本发明实施例的二次规划形式的电力系统抗差状态估计方法,具有良好的抗差性,并且所提出的模型属于二次规划问题,从数学上可保证得到全局最优解,计算效率更高,具有很好的工程应用前景。
申请公布号 CN103701115B 申请公布日期 2015.10.28
申请号 CN201310596140.5 申请日期 2013.11.22
申请人 清华大学;华北电力大学 发明人 刘锋;陈艳波;梅生伟;黄少伟;马进;崔静思
分类号 H02J3/00(2006.01)I;G06F17/50(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 张大威
主权项 一种二次规划形式的电力系统抗差状态估计方法,其特征在于,包括以下步骤:S1,提出二次规划形式的电力系统抗差状态估计基本模型,其中,所述二次规划形式的电力系统抗差状态估计基本模型为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><munder><mrow><mi>M</mi><mi>i</mi><mi>n</mi></mrow><mrow><mi>X</mi><mo>,</mo><mi>&sigma;</mi><mo>,</mo><mi>&xi;</mi><mo>,</mo><msup><mi>&xi;</mi><mo>*</mo></msup></mrow></munder><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>X</mi><mi>T</mi></msup><mi>X</mi><mo>+</mo><mfrac><mi>B</mi><mn>2</mn></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msubsup><mi>&sigma;</mi><mi>i</mi><mn>2</mn></msubsup><mo>+</mo><mi>C</mi><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>i</mi></msub><mo>+</mo><msubsup><mi>&xi;</mi><mi>i</mi><mo>*</mo></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000778673000000011.GIF" wi="841" he="140" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><msup><mi>X</mi><mi>T</mi></msup><msubsup><mi>J</mi><mi>i</mi><mi>T</mi></msubsup><mo>+</mo><msub><mi>&sigma;</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub><mo>&le;</mo><mi>&epsiv;</mi><mo>+</mo><msub><mi>&xi;</mi><mi>i</mi></msub><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>m</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msup><mi>X</mi><mi>T</mi></msup><msubsup><mi>J</mi><mi>i</mi><mi>T</mi></msubsup><mo>-</mo><msub><mi>&sigma;</mi><mi>i</mi></msub><mo>&le;</mo><mi>&epsiv;</mi><mo>+</mo><msubsup><mi>&xi;</mi><mi>i</mi><mo>*</mo></msubsup><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>m</mi></mrow></mtd></mtr><mtr><mtd><mrow><munder><mi>X</mi><mo>&OverBar;</mo></munder><mo>&le;</mo><mi>X</mi><mo>&le;</mo><mover><mi>X</mi><mo>&OverBar;</mo></mover></mrow></mtd></mtr><mtr><mtd><mrow><mi>&xi;</mi><mo>,</mo><msup><mi>&xi;</mi><mo>*</mo></msup><mo>&GreaterEqual;</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000778673000000012.GIF" wi="867" he="331" /></maths>其中:X∈R<sup>N+2b</sup>为中间状态变量,N为网络中所有节点的总数目,b为网络中所有支路的数目;m为量测量的总个数;B&gt;0;σ<sub>i</sub>为量测误差矢量σ的第i个分量;C&gt;0为惩罚因子;J<sub>i</sub>∈R<sup>m×(N+2b)</sup>为常数雅可比矩阵;<img file="FDA0000778673000000013.GIF" wi="376" he="81" /><img file="FDA0000778673000000014.GIF" wi="430" he="80" />y<sub>i</sub>∈R<sup>m</sup>为量测矢量,包括节点电压幅值的平方、支路有功、支路无功、注入有功、注入无功和支路电流幅值的平方;ε&gt;0为某个预先设定的参数;<img file="FDA00007786730000000111.GIF" wi="69" he="64" />和<img file="FDA00007786730000000112.GIF" wi="62" he="59" />分别代表X的下界和上界;S2,对所述二次规划形式的电力系统抗差状态估计基本模型进行转化,得到所述二次规划形式的电力系统抗差状态估计基本模型的对偶模型;以及S3,对所述二次规划形式的电力系统抗差状态估计基本模型的对偶模型求解,得到状态变量估计值。
地址 100084 北京市海淀区100084-82信箱