发明名称 一种基于无功网损的电力系统支路电抗参数辨识方法
摘要 本发明涉及一种基于无功网损的电力系统支路电抗参数辨识方法,属于电力系统调度自动化与电网仿真技术领域。首先通过电力系统量测数据,计算电力系统每条支路输送功率的容量与无功网损的量测值;并根据电力系统状态估计计算结果,得到每条支路输送功率的容量与无功网损的估计值;将量测值与估计值相减,得到条支路输送功率的容量与无功网损的残差。比较每条支路输送功率容量残差与无功网损残差的符号,如果符号相反,则说明该条支路的电抗参数需要进行校核。本发明方法的计算量少,实现起来简单方便;本方法对不良电抗的定位精确,由于无功网损和支路容量之间的关系是本支路电抗数值最直接的体现,能够实现对错误参数的准确定位。
申请公布号 CN102156770B 申请公布日期 2012.11.14
申请号 CN201110007359.8 申请日期 2011.01.14
申请人 江西省电力公司;清华大学 发明人 张伯明;吴文传;郭玉金;姚诸香;孙宏斌;郭烨;郭庆来
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 罗文群
主权项 1.一种基于无功网损的电力系统支路电抗参数辨识方法,其特征在于该方法包括以下步骤:(1)根据电力系统中量测系统的量测值,计算得到电网中每条支路的无功网损量测值和容量量测值,计算公式如下:<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>Ql</mi><mi>k</mi><mi>meas</mi></msubsup><mo>=</mo><mo>|</mo><msubsup><mi>Q</mi><mi>kij</mi><mi>meas</mi></msubsup><mo>+</mo><msubsup><mi>Q</mi><mi>kji</mi><mi>meas</mi></msubsup><mo>|</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msubsup><mi>S</mi><mi>k</mi><mi>meas</mi></msubsup><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msqrt><msup><mrow><mo>(</mo><msubsup><mi>P</mi><mi>kij</mi><mi>meas</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msubsup><mi>Q</mi><mi>kij</mi><mi>meas</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msqrt><msup><mrow><mo>(</mo><msubsup><mi>P</mi><mi>kji</mi><mi>meas</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msubsup><mi>Q</mi><mi>kji</mi><mi>meas</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow></math>]]></maths>其中上标meas表示量测值,下标k表示第k号支路,i、j分别为支路的两端节点,<img file="FDA00001666986900013.GIF" wi="140" he="60" />表示第k号支路的无功网损量测值,<img file="FDA00001666986900014.GIF" wi="113" he="59" />表示第k号支路的容量量测值,<img file="FDA00001666986900015.GIF" wi="313" he="73" />分别表示第k号支路由端节点i流向端节点j的有功功率量测值和无功功率量测值,<img file="FDA00001666986900016.GIF" wi="158" he="73" /><img file="FDA00001666986900017.GIF" wi="121" he="67" />表示第k号支路由端节点j流向端节点i的有功功率量测值和无功功率量测值;(2)根据电力系统中量测系统的量测值及电力系统网络模型和参数,进行状态估计,即求解如下的状态估计优化模型:<maths num="0003"><![CDATA[<math><mrow><mi>min</mi><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mfrac><mn>1</mn><msub><mi>R</mi><mi>i</mi></msub></mfrac><msup><mrow><mo>(</mo><msub><mi>z</mi><mi>i</mi></msub><mo>-</mo><msub><mi>h</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow></math>]]></maths>s.t    c(x)=0其中R<sub>i</sub>、z<sub>i</sub>、h<sub>i</sub>分别为第i号量测的量测标准差、量测值和量测函数;m为系统量测数,x为电力系统状态变量,由节点电压的幅值和相角组成,c(x)=0表示零注入等式约束,对以上的电力系统状态估计优化模型进行求解,得到电力系统状态估计结果;(3)根据上述电力系统状态估计结果,计算得到电网中每条支路的无功网损估计值和容量估计值,计算公式如下:<maths num="0004"><![CDATA[<math><mrow><msubsup><mi>Ql</mi><mi>k</mi><mi>se</mi></msubsup><mo>=</mo><mo>|</mo><msubsup><mi>Q</mi><mi>kij</mi><mi>se</mi></msubsup><mo>+</mo><msubsup><mi>Q</mi><mi>kji</mi><mi>se</mi></msubsup><mo>|</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msubsup><mi>S</mi><mi>k</mi><mi>se</mi></msubsup><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msqrt><msup><mrow><mo>(</mo><msubsup><mi>P</mi><mi>kij</mi><mi>se</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msubsup><mi>Q</mi><mi>kij</mi><mi>se</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msqrt><msup><mrow><mo>(</mo><msubsup><mi>P</mi><mi>kji</mi><mi>se</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msubsup><mi>Q</mi><mi>kji</mi><mi>se</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow></math>]]></maths>其中上标se表示状态估计值,下标k表示第k号支路,i,j分别为支路两端节点,<img file="FDA000016669869000111.GIF" wi="94" he="60" />表示第k号支路的无功网损状态估计值,<img file="FDA00001666986900021.GIF" wi="67" he="59" />表示第k号支路的容量状态估计值,<img file="FDA00001666986900022.GIF" wi="195" he="67" />分别表示第k号支路由端节点i流向端节点j的有功功率状态估计值和无功功率状态估计值,<img file="FDA00001666986900023.GIF" wi="212" he="67" />表示第k号支路由端节点j流向端节点i的有功功率状态估计值和无功功率状态估计值;(4)根据步骤(1)和步骤(3)的计算结果,计算电网中每条支路容量和无功网损的残差,计算公式如下:<maths num="0006"><![CDATA[<math><mrow><msubsup><mi>Ql</mi><mi>k</mi><mi>rad</mi></msubsup><mo>=</mo><msubsup><mi>Ql</mi><mi>k</mi><mi>meas</mi></msubsup><mo>-</mo><msubsup><mi>Ql</mi><mi>k</mi><mi>se</mi></msubsup></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><msubsup><mi>S</mi><mi>k</mi><mi>rad</mi></msubsup><mo>=</mo><msubsup><mi>S</mi><mi>k</mi><mi>meas</mi></msubsup><mo>-</mo><msubsup><mi>S</mi><mi>k</mi><mi>se</mi></msubsup></mrow></math>]]></maths>其中上标rad表示残差,<img file="FDA00001666986900026.GIF" wi="116" he="68" />和<img file="FDA00001666986900027.GIF" wi="91" he="68" />分别表示k号支路的无功网损残差和容量残差;(5)对电网中的每条支路进行判断,若<img file="FDA00001666986900028.GIF" wi="358" he="68" />则该条支路的电抗参数存在错误,需要进行校核,若<img file="FDA00001666986900029.GIF" wi="357" he="68" />则该条支路的电抗参数不存在错误。
地址 330077 江西省南昌市湖滨东路666号