发明名称 基于同步相角测量装置的电力系统动态状态估计方法
摘要 本发明涉及电力系统运行与控制技术领域。本发明公开了一种基于同步相角测量装置的电力系统动态状态估计方法。本发明的技术方案包括步骤:A、读取电力系统当前网络参数和网络拓扑结构,并由此形成节点导纳矩阵和支路-节点关联矩阵;B、根据电力系统网络拓扑结构建立等效电路,配置电力系统量测函数和PMU,系统的量测包括节点电压幅值量测、节点电流幅值量测、节点功率注入量测和节点潮流量测;C、在扩展卡尔曼滤波基础上对系统进行动态状态估计;D、收敛条件判断。本发明中PMU的引入能够为系统提供实时准确的电压、相角等量测信息,保证系统获得更高的量测冗余度,进而提高状态估计的精度。具有鲁棒性好,状态估计精度高,收敛性好的特点。
申请公布号 CN103326358A 申请公布日期 2013.09.25
申请号 CN201310238251.9 申请日期 2013.06.17
申请人 西南交通大学 发明人 张葛祥;赵俊博
分类号 H02J3/00(2006.01)I;G06Q50/06(2012.01)I 主分类号 H02J3/00(2006.01)I
代理机构 成都宏顺专利代理事务所(普通合伙) 51227 代理人 李顺德
主权项 1.基于同步相角测量装置的电力系统动态状态估计方法,其特征在于,包括以下步骤:A、读取电力系统当前网络参数和网络拓扑结构,并由此形成节点导纳矩阵和支路-节点关联矩阵;B、根据电力系统网络拓扑结构建立等效电路,配置电力系统量测函数和PMU,系统的量测包括节点电压幅值量测、节点电流幅值量测、节点功率注入量测和节点潮流量测;C、在扩展卡尔曼滤波基础上对系统进行动态状态估计一个电力系统的状态可以由一个包含一系列复杂电压幅值和相角的n维状态向量x来表示,系统通常使用状态估计每隔几分钟或者一定的采样间隔更新一次状态;在已知量测z<sub>k</sub>后,系统在第k次采样时的状态x<sub>k</sub>由下式表示:z<sub>k</sub>=h(x<sub>k</sub>)+v<sub>k</sub>其中,h(·)表示m维非线性函数向量;v<sub>k</sub>是服从正态分布的随机白噪声,即v<sub>k</sub>~N(0,R<sub>k</sub>),R<sub>k</sub>是量测误差的方差;对系统进行动态状态估计,主要包括参数辨识、状态预测和状态滤波:参数辨识:电力系统运行状态的变化可由以下线性化的准稳态模型来描述:x<sub>k+1</sub>=F<sub>k</sub>x<sub>k</sub>+G<sub>k</sub>+w<sub>k</sub>其中,F<sub>k</sub>,G<sub>k</sub>是n维非零对角动态模型参数矩阵;F<sub>k</sub>是状态转移矩阵;G<sub>k</sub>是控制向量;w<sub>k</sub>是服从正态分布的随机白噪声,即w<sub>k</sub>~N(0,Q<sub>k</sub>),Q<sub>k</sub>是系统模型误差的方差;参数F<sub>k</sub>,G<sub>k</sub>由Holt’s两参数线性指数平滑法求得;状态预测:一个含n节点系统的状态<img file="FDA00003353620000011.GIF" wi="90" he="86" />和协方差矩阵M<sub>k+1</sub>预测结果为:<maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>x</mi><mo>~</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>F</mi><mi>k</mi></msub><msub><mover><mi>x</mi><mo>^</mo></mover><mi>k</mi></msub><mo>+</mo><msub><mi>G</mi><mi>k</mi></msub></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>M</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>F</mi><mi>k</mi></msub><msub><mi>P</mi><mi>k</mi></msub><msubsup><mi>F</mi><mi>k</mi><mi>T</mi></msubsup><mo>+</mo><msub><mi>Q</mi><mi>k</mi></msub></mrow></math>]]></maths>其中,<img file="FDA00003353620000014.GIF" wi="66" he="86" />是一个(2n-1)×1维状态向量的估计值,P<sub>k</sub>是误差协方差关联矩阵;状态滤波:假设已经获得一组电力系统的实时量测值z<sub>k</sub>,则通过对预测的状态向量<img file="FDA00003353620000015.GIF" wi="90" he="80" />进行滤波可得到新的状态估计向量<img file="FDA00003353620000016.GIF" wi="110" he="92" />因此在k+1时刻,状态估计的目标函数为:<maths num="0003"><![CDATA[<math><mrow><msub><mi>J</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>[</mo><mi>z</mi><mo>-</mo><mi>h</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>]</mo></mrow><mi>T</mi></msup><msup><mi>R</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>[</mo><mi>z</mi><mo>-</mo><mi>h</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>]</mo><mo>+</mo><mo>[</mo><mi>x</mi><mo>-</mo><mover><mi>x</mi><mo>~</mo></mover><mo>]</mo><msup><mi>M</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>[</mo><mi>x</mi><mo>-</mo><mover><mi>x</mi><mo>~</mo></mover><mo>]</mo></mrow></math>]]></maths>其中,R为量测误差矩阵,W=R<sup>-1</sup>是一个对角权重矩阵,它的对角元素由每个量测的标准差决定;D、收敛条件判断若<img file="FDA00003353620000021.GIF" wi="344" he="145" />算法收敛,输出各个节点剔除不良数据之后系统更准确的电压幅值和相角,否则转步骤C。
地址 610031 四川省成都市二环路北一段111号