电力系统状态估计new.PPT

《电力系统状态估计new.PPT》由会员分享，可在线阅读，更多相关《电力系统状态估计new.PPT（152页珍藏版）》请在金锄头文库上搜索。

1、第二章第二章 电力系统状态估计电力系统状态估计一一概述概述 二二电力系统量测系统的数学描述与可观察性电力系统量测系统的数学描述与可观察性 三三电力系统状态估计的理论与计算方法电力系统状态估计的理论与计算方法 n最小二乘估计及改进最小二乘估计及改进n支路潮流状态估计法支路潮流状态估计法四四不良数据的检测与辨识不良数据的检测与辨识 五五. 网络结线分析及网络结构辨识的基本概念网络结线分析及网络结构辨识的基本概念 1 一一 常规状态估计的概念常规状态估计的概念 根据可获取的根据可获取的量测数据量测数据估算系统内部状估算系统内部状态的方法。态的方法。 由于随机干扰及测量误差的介入，无论是由于随机干扰及

2、测量误差的介入，无论是理想的运动方程或测量方程均不能求出精确的理想的运动方程或测量方程均不能求出精确的状态向量。只有通过统计学的方法加以处理以状态向量。只有通过统计学的方法加以处理以求出对状态向量的估计值。这种方法，称为状求出对状态向量的估计值。这种方法，称为状态估计。态估计。一概述一概述2 状态估计分为动态估计和静态估计两种。状态估计分为动态估计和静态估计两种。动态估计动态估计: :根据运动方程以某一时刻的测量数根据运动方程以某一时刻的测量数据作为初值进行下一个时刻状态量的估计，叫据作为初值进行下一个时刻状态量的估计，叫做动态估计；做动态估计；静态估计静态估计: :仅仅根据某时刻测量数据，确

3、定该仅仅根据某时刻测量数据，确定该时刻的状态量的估计，叫做静态估计。时刻的状态量的估计，叫做静态估计。本章介绍电力系统的静态估计。本章介绍电力系统的静态估计。一概述一概述3二二 电力系统状态估计的必要性电力系统状态估计的必要性 SCADA(SCADA(Supervisory Control And Data Acquisition)装置采集电网中的信息，并通过信息网装置采集电网中的信息，并通过信息网络将采集数据传送至能量控制中心的计算机监控系统。络将采集数据传送至能量控制中心的计算机监控系统。所获得的数据用于一系列应用程序，包括保证系统的所获得的数据用于一系列应用程序，包括保证系统的经济运行及

4、对系统发生设备或线路故障时进行安全性经济运行及对系统发生设备或线路故障时进行安全性评估分析，并最终构成了我们所称的能量管理系统评估分析，并最终构成了我们所称的能量管理系统（EMSEMS）。）。 电力系统状态估计（电力系统状态估计（POWER SYSTEM STATE POWER SYSTEM STATE ESTIMATIONESTIMATION）是）是EMSEMS中保证电力系统实时数据质量的中保证电力系统实时数据质量的重要一环，它为其它应用程序的实现奠定了基础。重要一环，它为其它应用程序的实现奠定了基础。 一概述一概述4n采集数据存在的问题采集数据存在的问题采集的数据是有误差的，不可靠（错误数

5、据）或者采集的数据是有误差的，不可靠（错误数据）或者局部信息不完整。局部信息不完整。模拟量模拟量母线电压、线路功率、负载功率。母线电压、线路功率、负载功率。一般要经过互感器、功率变换器、一般要经过互感器、功率变换器、A/D转换器量化成转换器量化成数字量，并通过通信传送到控制中心。数字量，并通过通信传送到控制中心。开关量开关量断路器、隔离开关等位置信息。断路器、隔离开关等位置信息。由于通信状态定义不一致造成开关位置错误。由于通信状态定义不一致造成开关位置错误。测量装置不全或种类限制。测量装置不全或种类限制。一概述一概述5n电力系统状态估计电力系统状态估计：对给定的系统结构及量测对给定的系统结构及

6、量测配置，在量测量有误差的情况下，估计出系统配置，在量测量有误差的情况下，估计出系统的真实状态的真实状态-各母线上的电压相角与模值及各各母线上的电压相角与模值及各元件上的潮流。元件上的潮流。作用：作用：提高数据精度，去除不良数据提高数据精度，去除不良数据计算出难以测量的电气量，相当于补充了量测量。计算出难以测量的电气量，相当于补充了量测量。 状态估计为建立一个高质量的数据库提供数据信状态估计为建立一个高质量的数据库提供数据信息，以便于进一步实现在线潮流、安全分析及经济调息，以便于进一步实现在线潮流、安全分析及经济调度等功能。度等功能。 一概述一概述6n实现状态估计需要的条件：实现状态估计需要的

7、条件： 1.1.量测冗余度：量测冗余度：量测冗余度是指量测量个数量测冗余度是指量测量个数m与与待估计的状态量个数待估计的状态量个数n之间的比值之间的比值m/n。 系统冗余度越高，对状态估计采用一定的估计方法系统冗余度越高，对状态估计采用一定的估计方法排除不良数据以及消除误差影响就越好。冗余量测的存排除不良数据以及消除误差影响就越好。冗余量测的存在是状态估计可以实现提高数据精度的基础。在是状态估计可以实现提高数据精度的基础。 2. 分析系统分析系统可观性可观性 当收集到的量测量通过量测方程能够覆盖所有母当收集到的量测量通过量测方程能够覆盖所有母线的电压幅值和相角时，则通过状态估计可以得到这些线的

8、电压幅值和相角时，则通过状态估计可以得到这些值，称该系统是可观测的值，称该系统是可观测的一概述一概述7 由于电力系统远动装置的工作情况经常变化，由于电力系统远动装置的工作情况经常变化，当远动信息量严重不足时，状态估计无法工作。当远动信息量严重不足时，状态估计无法工作。因此，在状态估计之前应先进行因此，在状态估计之前应先进行可观察性检验可观察性检验。如果系统中某些部分被判定是不可观察的，无如果系统中某些部分被判定是不可观察的，无法通过状态估计建立实时数据库，则应把它从法通过状态估计建立实时数据库，则应把它从状态估计的计算中退出来，或用增加人工设置状态估计的计算中退出来，或用增加人工设置的虚拟测量

9、或称伪测量数据来使它变成可观察的虚拟测量或称伪测量数据来使它变成可观察的。的。一概述一概述8n协同状态估计工作的是协同状态估计工作的是不良数据的检测与辨识不良数据的检测与辨识，如果有误差很大的，一般没有随机性的数据如果有误差很大的，一般没有随机性的数据（也称不良数据），就应该将它剔除，并重新（也称不良数据），就应该将它剔除，并重新进行状态估计，最终建立起完整的电力系统模进行状态估计，最终建立起完整的电力系统模型。型。 由于状态估计必须在几分钟内完成，因此它由于状态估计必须在几分钟内完成，因此它通常可以跟踪节点负荷的变化规律，在必要时通常可以跟踪节点负荷的变化规律，在必要时可用来提供补充的测量量

10、。因此，状态估计的可用来提供补充的测量量。因此，状态估计的计算结果也可以用于负荷预测。计算结果也可以用于负荷预测。一概述一概述9图图2-12-1电力系统状态估计的功能流程框图电力系统状态估计的功能流程框图 一概述一概述10 三三 状态估计与常规潮流计算的比较状态估计与常规潮流计算的比较图图2-2 2-2 状态估计与潮流计算的比较框图状态估计与潮流计算的比较框图(a)(a)潮流计算；潮流计算；(b)(b)状态估计状态估计一概述一概述11 潮流计算与状态估计的区别潮流计算与状态估计的区别 潮流计算方程式的数目潮流计算方程式的数目等于等于未知数的数目。未知数的数目。而状态估计的测量向量的维数一般而状

11、态估计的测量向量的维数一般大于大于未知状未知状态向量的维数，即方程数的个数多于未知数的态向量的维数，即方程数的个数多于未知数的个数。其中，测量向量可以是节点电压、节点个数。其中，测量向量可以是节点电压、节点注入功率、线路潮流等测量量的任意组合。注入功率、线路潮流等测量量的任意组合。两者求解的数学方法也不同。潮流计算一般两者求解的数学方法也不同。潮流计算一般用牛顿用牛顿- -拉夫逊法求解拉夫逊法求解 个非线性方程组。而个非线性方程组。而状态估计则是根据一定的估计准则，按估计理状态估计则是根据一定的估计准则，按估计理论的方法求解方程组。论的方法求解方程组。一概述一概述12 状态估计中的状态估计中的

12、“估计估计”不意味着不准不意味着不准确，相反，对于实际运行的系统来说，确，相反，对于实际运行的系统来说，不能认为潮流计算是绝对准确的，而状不能认为潮流计算是绝对准确的，而状态估计的值显然更准确。态估计的值显然更准确。状态估计可认为是一种广义潮流，而常状态估计可认为是一种广义潮流，而常规潮流计算是一种狭义潮流，及状态估规潮流计算是一种狭义潮流，及状态估计中计中m=nm=n的特例。的特例。一概述一概述13141516一一 电力系统测量系统的数学描述电力系统测量系统的数学描述n 电力系统的量测量：电力系统的量测量：二电力系统状态估计的数学描述与可观察性二电力系统状态估计的数学描述与可观察性式中，z为

13、量测向量，假设维数为m；Pij为支路ij有功潮流量测量；Qij为支路ij无功潮流量测量；Pi为母线i有功注入功率量测量；Qi为母线i无功注入功率量测量；Vi为母线i的电压幅值量测量。 17n待求状态变量，待求状态变量，如果系统结构与参数都已知，根据状态变量就如果系统结构与参数都已知，根据状态变量就不难求出各支路的有功潮流、无功潮流及所有不难求出各支路的有功潮流、无功潮流及所有节点的注入功率。节点的注入功率。二电力系统状态估计的数学描述与可观察性二电力系统状态估计的数学描述与可观察性 x = 式中，式中，x为状态向量，为状态向量， i为母线为母线i的电压相角的电压相角; Vi为为母线母线i的电压

14、幅值。的电压幅值。18n量测方程式量测方程式 即联系状态向量与测量量向量之间的函即联系状态向量与测量量向量之间的函数关系。在考虑有测量噪声时，它们之数关系。在考虑有测量噪声时，它们之间的关系为间的关系为 (2-1)(2-1) 式中：式中： 为为 维的测量量向量；维的测量量向量； 为测为测量函数向量量函数向量 (2-2)(2-2) 为测量噪声向量，其表达式为为测量噪声向量，其表达式为 (2-3)(2-3)二电力系统状态估计的数学描述与可观察性二电力系统状态估计的数学描述与可观察性19量测方程式：量测方程式： 二电力系统状态估计的数学描述与可观察性二电力系统状态估计的数学描述与可观察性h(x)=

15、节点电压测量方程式：状态变量与支路潮流节点电压测量方程式：状态变量与支路潮流的非线性函数表达式。的非线性函数表达式。注入功率测量方程式：节点注入功率与支路注入功率测量方程式：节点注入功率与支路潮流的非线性函数表达式。潮流的非线性函数表达式。20 表表2-12-1列出五种基本测量方式。第一种列出五种基本测量方式。第一种测量其维数为测量其维数为 ，显然没有冗余度，这在，显然没有冗余度，这在状态估计是不实际的。第五种测量方式具状态估计是不实际的。第五种测量方式具有最高的维数和冗余度，但所需投资太高，有最高的维数和冗余度，但所需投资太高，也是不现实的。因此，实际测量方式是第也是不现实的。因此，实际测量

16、方式是第一到第四的组合。一到第四的组合。二电力系统状态估计的数学描述与可观察性二电力系统状态估计的数学描述与可观察性21表表2-1 五种基本测量方式五种基本测量方式 测量方式的分量方程式 的维数（1）平衡节点除平衡节点外所有节点的注入功率 、 式(2-4) 、(2-5) 、(2-9)（2）（1）加上所有节点的电压模值 式(2-4) 、(2-5) 、(2-9)（3）支路两侧的有功、无功潮流 式(2-6) 、(2-7)（4）（3）加上所有节点的电压模值式(2-6) 、(2-7) 、(2-9)（5）完全的测量系统式(2-4) (2-7)、(2-9)二电力系统状态估计的数学描述与可观察性二电力系统状态

17、估计的数学描述与可观察性22 相应的方程式为相应的方程式为 (2-4)(2-4) (2-5) (2-5) (2-6) (2-6) (2-7) (2-7) (2-8) (2-8) (2-9) (2-9) 二电力系统状态估计的数学描述与可观察性二电力系统状态估计的数学描述与可观察性23图图2-3 2-3 形线路元件模型形线路元件模型二电力系统状态估计的数学描述与可观察性二电力系统状态估计的数学描述与可观察性24 用测量量来估计系统的状态存在若干用测量量来估计系统的状态存在若干不准确的因素，概括起来有以下几点。不准确的因素，概括起来有以下几点。 （1 1）数学模型不完善数学模型不完善。测量数学模型。

18、测量数学模型通常有工程性的近似处理。此外，还存在通常有工程性的近似处理。此外，还存在模型采用参数不精确的问题，另外，网络模型采用参数不精确的问题，另外，网络结构变化时，结构模型不能及时更新。上结构变化时，结构模型不能及时更新。上述问题属于参数不精确的，通常用参数估述问题属于参数不精确的，通常用参数估计方法解决；属于网络结构错误的，则采计方法解决；属于网络结构错误的，则采用网络接线错误的检测与辨识来解决。用网络接线错误的检测与辨识来解决。二电力系统状态估计的数学描述与可观察性二电力系统状态估计的数学描述与可观察性25 （2 2）测量系统的系统误差测量系统的系统误差。这是由于仪。这是由于仪表不精确

19、，通道不完善所引起的。它的特表不精确，通道不完善所引起的。它的特点是误差恒为正或负而没有随机性。一般点是误差恒为正或负而没有随机性。一般这类数据属于不良数据。清除这类误差的这类数据属于不良数据。清除这类误差的方法，主要是依靠提高测量系统的精确性方法，主要是依靠提高测量系统的精确性与可靠性，也可以用软件方法来检测与辨与可靠性，也可以用软件方法来检测与辨识出不良数据，并通过增加测量系统的冗识出不良数据，并通过增加测量系统的冗余度来补救，但这仅是一种辅助手段。余度来补救，但这仅是一种辅助手段。二电力系统状态估计的数学描述与可观察性二电力系统状态估计的数学描述与可观察性26 （3 3）随机误差随机误差

20、。这是测量系统中不可避。这是测量系统中不可避免会出现的。其特点是小误差比大误差出免会出现的。其特点是小误差比大误差出现的概率大，正负误差出现的概率相等，现的概率大，正负误差出现的概率相等，即概率密度曲线对称于零值或误差的数学即概率密度曲线对称于零值或误差的数学期望为零。状态估计式期望为零。状态估计式(2-1)(2-1)和式和式(2-3)(2-3)中中的误差向量的误差向量 就是这种误差。就是这种误差。二电力系统状态估计的数学描述与可观察性二电力系统状态估计的数学描述与可观察性27 测量的随机误差或噪声向量测量的随机误差或噪声向量 是均值是均值为零的高斯白噪声，其概率密度为为零的高斯白噪声，其概率

21、密度为 式中：式中： 是误差是误差 的标准差；方差的标准差；方差 越大越大表示误差大的概率增大。表示误差大的概率增大。二电力系统状态估计的数学描述与可观察性二电力系统状态估计的数学描述与可观察性28 由于误差的概率密度或协方差很难由测由于误差的概率密度或协方差很难由测量或计算确定，因此在实际应用中常用测量或计算确定，因此在实际应用中常用测量设备的误差来代替。测量误差的方差为量设备的误差来代替。测量误差的方差为 (2-(2-11)11) 式中：为仪表测量误差，一般取式中：为仪表测量误差，一般取0.01 0.02；为远动和模数转换的误差，一般取；为远动和模数转换的误差，一般取0.0025 0.00

22、5；为满刻度时的仪表误差；为满刻度时的仪表误差；为规格化因子。为规格化因子。二电力系统状态估计的数学描述与可观察性二电力系统状态估计的数学描述与可观察性29 每个测量量的方差为每个测量量的方差为 。测量误差。测量误差的方差阵，可以写成每个测量误差方差的的方差阵，可以写成每个测量误差方差的对角阵为对角阵为 (2-12)(2-12)二电力系统状态估计的数学描述与可观察性二电力系统状态估计的数学描述与可观察性30二二 电力系统的可观察性电力系统的可观察性 电力系统状态能够被表征的必要条件电力系统状态能够被表征的必要条件是它的可观察性。如果对系统进行有限次是它的可观察性。如果对系统进行有限次独立的观察

23、（测量），由这些观察向量所独立的观察（测量），由这些观察向量所确定的状态是唯一的，就称该系统是可观确定的状态是唯一的，就称该系统是可观察的。卡尔曼最初提出可观察的概念只是察的。卡尔曼最初提出可观察的概念只是在线性系统范围内，在电力系统的问题中在线性系统范围内，在电力系统的问题中可以由式可以由式(2-1)(2-1)的雅可比矩阵的雅可比矩阵 来确定来确定 (2-13)(2-13)二电力系统状态估计的数学描述与可观察性二电力系统状态估计的数学描述与可观察性31 只要只要 阶测量矩阵阶测量矩阵 的秩为的秩为 ，则系，则系统是可观察的，这表示通过测量量可以唯统是可观察的，这表示通过测量量可以唯一地确定系

24、统的状态量，或者说，测量点一地确定系统的状态量，或者说，测量点的数量及其分布可以保证系统是可观察的。的数量及其分布可以保证系统是可观察的。在非线性系统中，可观察性问题虽复杂得在非线性系统中，可观察性问题虽复杂得多，但可观察的一个必要但非充分条件仍多，但可观察的一个必要但非充分条件仍是雅可比矩阵是雅可比矩阵 的秩等于的秩等于 ，每一时刻的，每一时刻的测量量维数至少应与状态量的维数相等。测量量维数至少应与状态量的维数相等。二电力系统状态估计的数学描述与可观察性二电力系统状态估计的数学描述与可观察性32 电力系统测量需要有较大的冗余度。有电力系统测量需要有较大的冗余度。有冗余度的目的是提高测量系统的

25、可靠性和冗余度的目的是提高测量系统的可靠性和提高状态估计的精确度。保证可观察性是提高状态估计的精确度。保证可观察性是测量点布置的最低要求。测量点布置的最低要求。 前面说过，电力系统出现异常大误差的前面说过，电力系统出现异常大误差的数据，称为不良数据。查找出不良数据，数据，称为不良数据。查找出不良数据，并将其剔除是建立实时数据库的基本要求。并将其剔除是建立实时数据库的基本要求。测量具有冗余度则是实现这一工作的基本测量具有冗余度则是实现这一工作的基本条件。条件。二电力系统状态估计的数学描述与可观察性二电力系统状态估计的数学描述与可观察性33状态估计算法分类：状态估计算法分类： 最小二乘法（最小二乘

26、法（静态估计静态估计） 加权最小二乘法加权最小二乘法 快速解耦状态估计快速解耦状态估计 正交变换法正交变换法 支路潮流状态估计法支路潮流状态估计法 递推状态估计（动态估计）递推状态估计（动态估计）三三 状态估计的算法状态估计的算法34一一 基本原理基本原理 静态估计是用一定的统计学准则通过静态估计是用一定的统计学准则通过测量向量测量向量 求出状态向量求出状态向量 ，并使之尽量，并使之尽量接近其真值接近其真值 。 是一个估计值，估计值是一个估计值，估计值与真值之间的误差称为估计误差与真值之间的误差称为估计误差 (2-14)(2-14) 估计误差值估计误差值 是是 维向量。判断估计维向量。判断估计

27、方法的优劣不是根据方法的优劣不是根据 中个别分量的估计中个别分量的估计误差值，而是根据误差值，而是根据 的整个统计特性来决的整个统计特性来决定。定。三最小二乘估计三最小二乘估计35 最小二乘估计是在电力系统状态估计最小二乘估计是在电力系统状态估计中应用最为广泛的方法之一。它以测量中应用最为广泛的方法之一。它以测量值值 和测量估计值和测量估计值 之差的平方和最小之差的平方和最小为目标准则，即为目标准则，即 三最小二乘估计三最小二乘估计36 应用在电力系统，状态估计是按测量值应用在电力系统，状态估计是按测量值 与系统数学模型确定的值与系统数学模型确定的值 的误差平方的误差平方和最小来确定的系统状态

28、和最小来确定的系统状态 ，即目标函，即目标函数为数为 (2-15)(2-15)三最小二乘估计三最小二乘估计37 二二 加权的意义加权的意义 上述方法对于任一个测量分量的误差上述方法对于任一个测量分量的误差 都以相同的机会加进目标函数，都以相同的机会加进目标函数，即它们在目标函数中所占的份额一样。但即它们在目标函数中所占的份额一样。但由于各个测量量的量测由于各个测量量的量测精度不一致精度不一致，因此，因此它们以同样的权重组成目标函数是不合理它们以同样的权重组成目标函数是不合理的。为提高整个估计值的精度，应该使各的。为提高整个估计值的精度，应该使各个量测量各取一个权值，精度高的测量量个量测量各取一

29、个权值，精度高的测量量权大一些，而精度低的测量量权小一些。权大一些，而精度低的测量量权小一些。根据这一原理提出了加权最小二乘准则。根据这一原理提出了加权最小二乘准则。三最小二乘估计三最小二乘估计38 加权最小二乘准则的目标函数为加权最小二乘准则的目标函数为 (2-(2-16)16) 式中：式中： 为一适当选择的正定阵，当为一适当选择的正定阵，当 为为单位阵时单位阵时(2-16)(2-16)就是最小二乘准则。就是最小二乘准则。 假设假设 ， 为式为式(2-12)(2-12)的测量误差方的测量误差方差阵。其中各元素为差阵。其中各元素为 于是目标函数于是目标函数可写成可写成 (2-(2-17)17)

30、三最小二乘估计三最小二乘估计39 三三 最小二乘算法最小二乘算法 为线性函数为线性函数 先假定先假定 是线性向量函数。是线性向量函数。 (2-18)(2-18) 或或 式中：式中： 为为 矩阵，其元素为矩阵，其元素为 。 状态量的值状态量的值 与测量值与测量值 的关系为的关系为三最小二乘估计三最小二乘估计40 按最小二乘准则建立目标函数按最小二乘准则建立目标函数 或或 (2-19) (2-19) 对目标函数求导数并取为零，即对目标函数求导数并取为零，即 (2-20)(2-20) 。三最小二乘估计三最小二乘估计41 亦即亦即 这是一组有这是一组有 个未知数的个未知数的 维方程组，维方程组，联立求

31、解即可求得联立求解即可求得 的最佳估计值的最佳估计值 。三最小二乘估计三最小二乘估计42 写成矩阵方程式的形式，写成矩阵方程式的形式， (2-21)(2-21) ：称为：称为信息矩阵信息矩阵三最小二乘估计三最小二乘估计：称为最佳估计值：称为最佳估计值43 状态估计值的误差：状态估计值的误差： (2-21)(2-21) 测量量的测量量的测量值测量值与与估计值估计值的差称为的差称为残差残差残差残差方程方程称为称为残差灵敏度矩阵残差灵敏度矩阵三最小二乘估计三最小二乘估计44 为非线性函数为非线性函数 以上是在以上是在 为线性函数的前提下讨沦为线性函数的前提下讨沦的。但电力系统的测量函数向量的。但电力

32、系统的测量函数向量 是非是非线性的向量函数，这时无法直接由目标函线性的向量函数，这时无法直接由目标函数数 的极值条件求解的极值条件求解 ，需要用迭代的，需要用迭代的方法求解。方法求解。 1 1 设状态变量的初值为设状态变量的初值为 将将 在在 处线性化，并用泰勒级数在处线性化，并用泰勒级数在附近展开，即附近展开，即 (2-27)(2-27)三最小二乘估计三最小二乘估计45 是函数向量是函数向量 的雅可比矩阵，其的雅可比矩阵，其元素为元素为 (2-28)(2-28) 2 2 目标函数目标函数 略去略去 的高阶项，取目标函数为的高阶项，取目标函数为 (2-29)(2-29) 取取 ，有，有 (2-

33、30)(2-30)三最小二乘估计三最小二乘估计46 3 3 极值条件极值条件即即则则式中式中 由此可得由此可得 (2-31)(2-31)三最小二乘估计三最小二乘估计47 4 4 迭代格式迭代格式 当当 充分接近充分接近 时泰勒级数略去高阶项时泰勒级数略去高阶项后才是足够近似的。用式后才是足够近似的。用式(2-31)(2-31)作逐次迭作逐次迭代，可以得到代，可以得到 。若以。若以 表示迭代序号，表示迭代序号，式式(2-31)(2-31)可以写成可以写成 (2-32)(2-32) (2-33) (2-33)三最小二乘估计三最小二乘估计48 5 5 收敛判据收敛判据 按式按式(2-32)(2-32

34、)和式和式(2-33)(2-33)进行迭代修正，进行迭代修正，直到目标函数接近于最小为止。所采用的直到目标函数接近于最小为止。所采用的收敛判据可以是以下三项中的任一项收敛判据可以是以下三项中的任一项 （1 1） (2-34)(2-34) （2 2） (2-35)(2-35) （3 3） (2-36)(2-36)三最小二乘估计三最小二乘估计49 上三式是三种收敛标准。其中式上三式是三种收敛标准。其中式(2-34)(2-34)表示状态修正量绝对值最大者小于规定的表示状态修正量绝对值最大者小于规定的收敛标准，这是最常用的判据。收敛标准，这是最常用的判据。 可取基可取基准电压模值的准电压模值的 。 满

35、足收敛标准时的满足收敛标准时的 即为最优状态估即为最优状态估计值计值 。此时测量量的估计值是。此时测量量的估计值是 。三最小二乘估计三最小二乘估计506 6 状态估计的计算步骤及程序框图状态估计的计算步骤及程序框图 当当 是是 的非线性函数时，进行状态的非线性函数时，进行状态估计的步骤如下：估计的步骤如下： 1)1)从状态量的初值计算测量函数向量从状态量的初值计算测量函数向量 和雅可比矩阵和雅可比矩阵 。 2)2)由遥测量由遥测量 和和 计算残差计算残差 和目和目标函数标函数 ，并由雅可比矩阵，并由雅可比矩阵 计算信计算信息矩阵息矩阵 和向量和向量 。三最小二乘估计三最小二乘估计51 3) 3

36、)解方程式解方程式(2-32)(2-32)求得状态修正量求得状态修正量 ，并取其中绝对值最大者并取其中绝对值最大者 。 4)4)检查是否达到收敛标准。检查是否达到收敛标准。 5)5)若未达到收敛标准，修改状态量，若未达到收敛标准，修改状态量， 继续迭代计算，直到收敛为止。继续迭代计算，直到收敛为止。 6)6)将计算结果送入不良数据检测与辨识将计算结果送入不良数据检测与辨识入口。入口。三最小二乘估计三最小二乘估计52 图图2-42-4是加权最小二乘估计程序框图，是加权最小二乘估计程序框图，其中框其中框1 1包括输入各测量量的权值。框包括输入各测量量的权值。框1 1的的初值在实际应用中一般取前一次

37、状态估计初值在实际应用中一般取前一次状态估计的电压值，以加快迭代的收敛速度。框的电压值，以加快迭代的收敛速度。框3 3中用现有的状态量中用现有的状态量 （如电压模值与电压（如电压模值与电压相角）计算相角）计算 及其偏导数及其偏导数 。框。框4 4求求解电压模值与相角的修正量，选出解电压模值与相角的修正量，选出 及及 ，供框，供框5 5作收敛检查。框作收敛检查。框6 6转入下一转入下一次迭代并对状态变量作修正。次迭代并对状态变量作修正。三最小二乘估计三最小二乘估计53图图2-4 2-4 加权最小二乘估计框图加权最小二乘估计框图 三最小二乘估计三最小二乘估计入口输入测量信息给定初值计算计算计算解线

38、性方程式(2-30) 求及l = l + 1到不良数据检测与辨识入口l154 四四 信息矩阵（信息矩阵（ 阵）的特点阵）的特点 稀疏性和对称性稀疏性和对称性 因为因为 一般为稀疏矩阵，所以可以一般为稀疏矩阵，所以可以用稀疏矩阵技巧进行求解。以下先讨论这用稀疏矩阵技巧进行求解。以下先讨论这个矩阵的结构，由式个矩阵的结构，由式(2-32)(2-32)可得可得 (2-38)(2-38) 或写成或写成 (2-39)(2-39) 为了求解式为了求解式(2-39)(2-39)，先研究，先研究 阵的特点。阵的特点。三最小二乘估计三最小二乘估计55 阵的元素：阵的元素： 因为因为 是对角阵，所以是对角阵，所以

39、 阵的结构与阵的结构与 的结构一致。由于的结构一致。由于 是稀疏的，而且是稀疏的，而且 和和 换位并不影响换位并不影响 的值，因此的值，因此 阵是阵是 的的对称稀疏矩阵对称稀疏矩阵。 阵的结构与导纳矩阵不一样，阵的结构与导纳矩阵不一样，取决于取决于网络结构与测点的布置。网络结构与测点的布置。三最小二乘估计三最小二乘估计56 的结构与网络结构和测点配置的关系的结构与网络结构和测点配置的关系 三最小二乘估计三最小二乘估计n支路功率测量支路功率测量： 对连接两个节点的支路路，不论在线路哪一侧，也不对连接两个节点的支路路，不论在线路哪一侧，也不论是有功或无功，只要有一个测量就能出现论是有功或无功，只要

40、有一个测量就能出现aij元素。元素。n节点注入功率测量：节点注入功率测量：对节点对节点i的有功或无功注入的测量值，不仅与节点的有功或无功注入的测量值，不仅与节点i的的状态量有关，而且还与同节点状态量有关，而且还与同节点i有直接连接的相邻节有直接连接的相邻节点的状态量有关。点的状态量有关。n3.节点电压测量节点电压测量 节点节点i的电压测量值仅在的电压测量值仅在H阵阵i列有非零元素，在列有非零元素，在A阵阵 中也只影响相应的中也只影响相应的i行对角元行对角元57 2 2节点注入功率测量节点注入功率测量 节点节点 的有功或无功注入的测量值，不的有功或无功注入的测量值，不仅与节点仅与节点 的状态量有

41、关，而且还与同节的状态量有关，而且还与同节点点 直接连接的相邻节点的状态量有关。直接连接的相邻节点的状态量有关。 对于图对于图2-52-5例子，在例子，在 阵中，相应于节阵中，相应于节点点 注入测量的行（设为注入测量的行（设为 行）的行）的 列以列以及与及与 相关的各节点（如相关的各节点（如 、）的列均、）的列均为非零元素，即为非零元素，即 、 、 、 为非零元素，为非零元素，即相应的阵为即相应的阵为三最小二乘估计三最小二乘估计58图图2-5 2-5 节点注入对节点注入对H H阵影响示意图阵影响示意图iekj三最小二乘估计三最小二乘估计59 根据式根据式(2-38)(2-38)可看出，这一测量

42、值，在可看出，这一测量值，在 阵中将使阵中将使 、 、 六个非对角元六个非对角元发生变化（由于发生变化（由于 是对称阵，这里仅列出是对称阵，这里仅列出下三角部分）并成为非零元素（这时下三角部分）并成为非零元素（这时 、 、均非零）。、均非零）。它的作用相当于在它的作用相当于在 、 、 、 、 、 六条支路上装有测量，而实际上图六条支路上装有测量，而实际上图2-52-5中以中以虚线表示的线路是不存在的。虚线表示的线路是不存在的。 三最小二乘估计三最小二乘估计60 根据上述，对于图根据上述，对于图2-6(a)2-6(a)的网络与测的网络与测点布置情况，其点布置情况，其 阵的结构如图阵的结构如图2-

43、6(b)2-6(b)所所示，其中列号为节点号，亦即该节点的状示，其中列号为节点号，亦即该节点的状态量电压模值与电压相角的序号。态量电压模值与电压相角的序号。 图中有图中有9 9个测量量，个测量量，7 7个状态量。个状态量。三最小二乘估计三最小二乘估计61图图2-6 2-6 信息矩阵的结构示意图信息矩阵的结构示意图(a)(a)系统网络示意图系统网络示意图;(b);(b)H H矩阵矩阵;(c);(c)A A矩阵矩阵;(d);(d)A A矩阵网络示意图矩阵网络示意图电压测量；电压测量；支路功率测量；支路功率测量；注入功率测量注入功率测量三最小二乘估计三最小二乘估计432（a）1（d）4321（c）

44、（b）P12 Q12P23 Q23P34 Q34V1P3 Q3 62 由式由式(2-38)(2-38)， ， 阵结构如图阵结构如图2-2-6(c)6(c)所示。用图所示。用图2-6(c)2-6(c)的关联关系可以绘的关联关系可以绘出出 阵的线图阵的线图2-6(d)2-6(d)。 比较图比较图2-6(a)2-6(a)与图与图2-6(d)2-6(d)可见，可见，n 凡没有配置支路功率测量，且其两侧又无凡没有配置支路功率测量，且其两侧又无注入功率的，其注入功率的，其 阵的阵的 。n 如果在节点如果在节点 有注入功率测量，则与有注入功率测量，则与 有有关联的各节点间就形成一闭合的回路。关联的各节点间就

45、形成一闭合的回路。 三最小二乘估计三最小二乘估计63四静态最小二乘估计的改进四静态最小二乘估计的改进思想：思想：n有功与无功的分解。有功与电压模值，有功与无功的分解。有功与电压模值，无功与电压相角间联系很弱。无功与电压相角间联系很弱。减少内存，提高每次迭代速度。减少内存，提高每次迭代速度。但增加迭代次数但增加迭代次数n信息矩阵常数化信息矩阵常数化进行一次因子分解进行一次因子分解。对。对角化角化提高计算效率提高计算效率一、快速解耦（分解）状态估计一、快速解耦（分解）状态估计64 数学模型数学模型 对极坐标形式的电力系统加权最小二乘对极坐标形式的电力系统加权最小二乘状态估计基本算法进行简化。状态估

46、计基本算法进行简化。 将状态变量按节点电压相角和模值分别将状态变量按节点电压相角和模值分别排列，即排列，即 将测量量按有功和无功分别排列，即将测量量按有功和无功分别排列，即 四静态最小二乘估计的改进四静态最小二乘估计的改进65 式中：表示支路有功潮流、节点有功注入式中：表示支路有功潮流、节点有功注入测量量向量；表示支路无功潮流、节点无测量量向量；表示支路无功潮流、节点无功注入、节点电压模值的测量向量。功注入、节点电压模值的测量向量。 雅可比矩阵可表示为雅可比矩阵可表示为 (2-40)(2-40) 同时，对角权矩阵也相应地按有功和无功同时，对角权矩阵也相应地按有功和无功分别排列，即分别排列，即

47、(2-41)(2-41)四静态最小二乘估计的改进四静态最小二乘估计的改进66 信息矩阵可以写成信息矩阵可以写成 (2-42)(2-42) 在高压电网中，有功主要取决于节点电在高压电网中，有功主要取决于节点电压相角，无功主要取决于节点电压模值。压相角，无功主要取决于节点电压模值。即即四静态最小二乘估计的改进四静态最小二乘估计的改进67 因此，可引入第一项简化假设：因此，可引入第一项简化假设： 这样，这样， 矩阵变为准对角阵：矩阵变为准对角阵： 式式(2-42)(2-42)可以转化为对角矩阵可以转化为对角矩阵四静态最小二乘估计的改进四静态最小二乘估计的改进68 如再假定各支路电阻远远小于电抗，支如

48、再假定各支路电阻远远小于电抗，支路两端的相角差很小，各节点电压模值接路两端的相角差很小，各节点电压模值接近于参考节点电压，即近于参考节点电压，即 ， ， ， 这样有这样有 取支路电抗倒数（不计变压器非标准取支路电抗倒数（不计变压器非标准变比及线路对地电容的影响）；变比及线路对地电容的影响）； 取支路导纳的虚部（电压测量的取支路导纳的虚部（电压测量的 元元素取素取 ）。）。四静态最小二乘估计的改进四静态最小二乘估计的改进69 于是信息矩阵就成为常数矩阵，不必在于是信息矩阵就成为常数矩阵，不必在迭代过程中修改。迭代过程中修改。 (2-43)(2-43) 对于修正方程右边项的处理不同于常规对于修正方

49、程右边项的处理不同于常规潮流计算。计算经验表明：潮流计算。计算经验表明： 矩阵元素采矩阵元素采用上述两项假设比准确计算更有利于收敛用上述两项假设比准确计算更有利于收敛性的改善，提高迭代计算的速度。性的改善，提高迭代计算的速度。四静态最小二乘估计的改进四静态最小二乘估计的改进70 迭代的修正方程式可以写成迭代的修正方程式可以写成 (2-44)(2-44) (2-45) (2-45) 展开为展开为 (2-46)(2-46) (2-47) (2-47)四静态最小二乘估计的改进四静态最小二乘估计的改进71 其中：其中： (2-48)(2-48) (2-49) (2-49) 式中：式中： 为节点电压相角

50、的向量，为节点电压相角的向量， 为节点为节点电压模值的向量。电压模值的向量。 方程式方程式(2-44)(2-44)和式和式(2-45)(2-45)的方法，称为的方法，称为快速解耦状态估计算法。快速解耦状态估计算法。四静态最小二乘估计的改进四静态最小二乘估计的改进72 当有功测量的维数为当有功测量的维数为 ，无功测量的维，无功测量的维数为数为 时，状态量时，状态量 、 的维数是网络节的维数是网络节点数中减去平衡节点的状态量数，分别为点数中减去平衡节点的状态量数，分别为 、 ，于是，于是 是是 阶的，阶的， 是是 阶的，阶的， 是是 阶常数对称矩阵，阶常数对称矩阵， 是是 阶常数阶常数对称矩阵，对

51、称矩阵， 是是 维向量，维向量， 是是 维向量。维向量。四静态最小二乘估计的改进四静态最小二乘估计的改进73 为进一步加快速度，可对式为进一步加快速度，可对式(2-44)(2-44)和式和式(2-45)(2-45)右边也做类似简化。这种方法，也右边也做类似简化。这种方法，也称为模分解估计算法。其简化式为称为模分解估计算法。其简化式为 (2-50)(2-50) (2-51) (2-51)四静态最小二乘估计的改进四静态最小二乘估计的改进74 支路潮流状态估计法是早期一种较成支路潮流状态估计法是早期一种较成功的算法。用这种算法进行状态估计功的算法。用这种算法进行状态估计仅需仅需支路潮流测量量，支路潮

52、流测量量，在状态估计计算时在状态估计计算时将将支支路功率路功率转换成转换成支路两端电压差支路两端电压差的量，最后的量，最后得到与基本加权最小二乘估计相类似的迭得到与基本加权最小二乘估计相类似的迭代修正公式。代修正公式。 由于这种方法只用支路测量量，所以由于这种方法只用支路测量量，所以也称也称“唯支路法唯支路法”。又由于这种方法需用。又由于这种方法需用支路测量量转变成支路两端电压差的量，支路测量量转变成支路两端电压差的量，所以又称所以又称“量测量变换法量测量变换法”。五支路潮流状态估计法五支路潮流状态估计法75 一一 数学模型数学模型 支路潮流测量量支路潮流测量量 ，表示连接节点，表示连接节点

53、、的支路的支路 上测量到的复功率。若应用加权上测量到的复功率。若应用加权最小二乘的算式，其目标函数为最小二乘的算式，其目标函数为 测量量向量测量量向量 是支路复潮流，以是支路复潮流，以 表示；表示； 测量函数向量测量函数向量 用用 表示，它是状态表示，它是状态向量向量- -节点复电压节点复电压 的函数。的函数。五支路潮流状态估计法五支路潮流状态估计法76 这样，状态估计目标函数为这样，状态估计目标函数为 (2-65)(2-65) 式中，式中， 是是 维实数对角矩阵。维实数对角矩阵。 若该支路若该支路 两端的节点电压为两端的节点电压为 ，则，则该支路两端的电压差为该支路两端的电压差为 (2-66

54、)(2-66) 则支路则支路 、的测量方程可写为：、的测量方程可写为： (2-67)(2-67)五支路潮流状态估计法五支路潮流状态估计法77 式中：式中： 为支路为支路 的的 侧复功率；、为图侧复功率；、为图2-112-11中等值电路的参数。中等值电路的参数。n 对于线路：对于线路：n 对于变压器：对于变压器：五支路潮流状态估计法五支路潮流状态估计法78图图2-11 2-11 输电线与变压器等值电路输电线与变压器等值电路 (a)(a)输电线输电线;(b);(b)变压器变压器五支路潮流状态估计法五支路潮流状态估计法79n潮流测量量潮流测量量 与经转换得到的与之相应的与经转换得到的与之相应的支路电

55、压差的关系为支路电压差的关系为 (2-68)(2-68) 由此可求出以支路潮流测量值表示的由此可求出以支路潮流测量值表示的该支路电压差该支路电压差“测量值测量值”： 五支路潮流状态估计法五支路潮流状态估计法80线路线路 侧：侧： 线路线路 侧：侧： 变压器变压器 侧：侧： 变压器变压器 侧：侧：五支路潮流状态估计法五支路潮流状态估计法81 若定义这个线路电压差的向量为功率若定义这个线路电压差的向量为功率测量变换来的电压差向量测量变换来的电压差向量 ，则写成矩，则写成矩阵形式后为阵形式后为 (2-69)(2-69) 即即 (2-(2-70)70) 式中：支路功率测量向量式中：支路功率测量向量 与

56、对应的支路与对应的支路电压差值向量电压差值向量 是是 维的；为维的；为 阶的阶的对角矩阵对角矩阵. .五支路潮流状态估计法五支路潮流状态估计法82 的元素为的元素为 (2-71)(2-71)五支路潮流状态估计法五支路潮流状态估计法83 为复数向量，其元素为为复数向量，其元素为 (2-72)(2-72) 上两式中：上两式中： 表示测点号。表示测点号。五支路潮流状态估计法五支路潮流状态估计法84n 当支路功率的估计值为当支路功率的估计值为 ，与之相应的，与之相应的支路电压值为支路电压值为 ，它们之间的关系也可，它们之间的关系也可以用式以用式(2-70)(2-70)表示表示 (2-73)(2-73)

57、 将式将式(2-70)(2-70)、式、式(2-73)(2-73)代入式代入式(2-65)(2-65)得得 五支路潮流状态估计法五支路潮流状态估计法85 至此，通过测量量的变换，原来以支至此，通过测量量的变换，原来以支路潮流表示的目标函数已经化成了以支路潮流表示的目标函数已经化成了以支路电压差来表示的目标函数。路电压差来表示的目标函数。五支路潮流状态估计法五支路潮流状态估计法上式中，电压差值的估计值向量上式中，电压差值的估计值向量 可以分别用平衡节点电压可以分别用平衡节点电压 与其余节与其余节点的电压向量点的电压向量 表示，即表示，即86 (2-75)(2-75) 式中：为测量点所在支路与节点

58、的关联矩式中：为测量点所在支路与节点的关联矩阵。由于平衡节点的状态是给定的。所阵。由于平衡节点的状态是给定的。所以可以把以可以把 阵写成阵写成 (2-(2-76)76)五支路潮流状态估计法五支路潮流状态估计法87 A A矩阵各行矩阵各行+1+1和和-1-1表示对应于每一个表示对应于每一个支路测量点，各有两个非零元素。若测支路测量点，各有两个非零元素。若测量点在支路量点在支路ijij的的i i侧，则侧，则i i列为列为+1 , j+1 , j列列为为-1-1。若在。若在j j侧，则反之；其余元素均为侧，则反之；其余元素均为零。零。 当线路两端均有测量点时，此线路将当线路两端均有测量点时，此线路将

59、在在A A中出现两次。当节点数为中出现两次。当节点数为 时，为时，为 阶矩阵。阶矩阵。五支路潮流状态估计法五支路潮流状态估计法88于是式于是式(2-75)(2-75)可写成可写成 (2-77)(2-77) 将将(2-77)(2-77)、式、式(2-76)(2-76)代入式代入式(2-74)(2-74)得得 式中，式中， 为为 阶实对角矩阵阶实对角矩阵. .五支路潮流状态估计法五支路潮流状态估计法89 若以若以 表示测量点号，其元素为表示测量点号，其元素为 在求目标函数最小化时，可以在求目标函数最小化时，可以假定假定在电在电力系统运行中电压变化不大，因此力系统运行中电压变化不大，因此 可以可以取

60、为取为常数矩阵常数矩阵。 五支路潮流状态估计法五支路潮流状态估计法90 于是于是, ,极值条件：极值条件： 因为因为 所以所以 (2-79)(2-79) 由上式求出的解是次优解，因为电由上式求出的解是次优解，因为电压差是由线路功率计算出的，而功率与压差是由线路功率计算出的，而功率与节点电压是相关的，需要迭代。节点电压是相关的，需要迭代。 五支路潮流状态估计法五支路潮流状态估计法假定：假定： 不依赖于节点电压不依赖于节点电压。91迭代格式为：迭代格式为：(2-80)(2-80) 若己知上一次迭代的若己知上一次迭代的 值，就可以求值，就可以求解下一次迭代的解下一次迭代的 值。式中的值。式中的 由上

61、一由上一次迭代求出：次迭代求出： 式中：式中： 和和 与与 有关。有关。五支路潮流状态估计法五支路潮流状态估计法92 二二 算法特点算法特点 支路潮流估计法与最小二乘估计法的支路潮流估计法与最小二乘估计法的差别在于：差别在于： 1.1.信息矩阵信息矩阵 中，是测点与节点的关中，是测点与节点的关联矩阵，是常数对角阵，信息矩阵的结联矩阵，是常数对角阵，信息矩阵的结构与节点导纳矩阵结构完全相同，因而构与节点导纳矩阵结构完全相同，因而稀疏程度高，程序设计方便。稀疏程度高，程序设计方便。 2.2.信息矩阵的常数化，可在网络结构和信息矩阵的常数化，可在网络结构和测点配置不变时只进行一次因子分解，测点配置不

62、变时只进行一次因子分解，从而节约每次迭代的计算量，同快速分从而节约每次迭代的计算量，同快速分解法一样。解法一样。五支路潮流状态估计法五支路潮流状态估计法93 3. 3.实部和虚部的迭代用同一个实数信实部和虚部的迭代用同一个实数信息矩阵息矩阵 ，既节约内存，又节约矩阵分，既节约内存，又节约矩阵分解计算时间解计算时间。 4.4.只能处理支路潮流测量量，而不能处只能处理支路潮流测量量，而不能处理其他形式的测量量，如节点注入功率。理其他形式的测量量，如节点注入功率。在实际系统中，因不能充分利用全部测在实际系统中，因不能充分利用全部测量量而降低了状态估计结果的可信度。量量而降低了状态估计结果的可信度。五

63、支路潮流状态估计法五支路潮流状态估计法5 解出的解出的 是待求量，而不是修正量。是待求量，而不是修正量。94 支路潮流状态估计程序框图，如图支路潮流状态估计程序框图，如图2-12 2-12 所示。其计算步骤如下：所示。其计算步骤如下： (1)(1)给定节点电压向量的初值给定节点电压向量的初值 ，可取所，可取所有节点电压与平衡节点的电压相同。有节点电压与平衡节点的电压相同。 (2)(2)利用测量量计算支路电压差值利用测量量计算支路电压差值 (3)(3)利用迭代方程式利用迭代方程式(2-80)(2-80)求解求解 。 (4)(4)重复步骤重复步骤(2)(2)、(3) (3) ，直至符合收敛，直至符

64、合收敛条件条件 (2-(2-81)81)五支路潮流状态估计法五支路潮流状态估计法95图图2-12支路潮流状态估计程序框图支路潮流状态估计程序框图 五支路潮流状态估计法五支路潮流状态估计法96三种方法综合比较：三种方法综合比较：n加权最小二乘法：估计质量好，收敛性好，加权最小二乘法：估计质量好，收敛性好，但内存需要大，计算时间长，适用于小型电但内存需要大，计算时间长，适用于小型电力系统。力系统。n快速分解法：估计质量好，收敛性好，计算快速分解法：估计质量好，收敛性好，计算速度快，程序复杂，是一种实用方法。速度快，程序复杂，是一种实用方法。n量测量变换法：计算速度快，内存需要少，量测量变换法：计算

65、速度快，内存需要少，但仅适用于支路型量测系统中。但仅适用于支路型量测系统中。五支路潮流状态估计法五支路潮流状态估计法97七不良数据的检测与辨识七不良数据的检测与辨识 一一 不良数据的检测不良数据的检测 电力系统的测量信息如果误差不大，测电力系统的测量信息如果误差不大，测量系统的配置恰当，则用一般的状态估计量系统的配置恰当，则用一般的状态估计方法可以得到满意的实时数据库。如果调方法可以得到满意的实时数据库。如果调度中心收到的远动测量数据具有度中心收到的远动测量数据具有异常大异常大的误差的误差，则常规状态估计算法无法估计出，则常规状态估计算法无法估计出正确的数值，影响电力系统的实时调度管正确的数值

66、，影响电力系统的实时调度管理。理。98七不良数据的检测与辨识七不良数据的检测与辨识 电力系统中测量系统的标准误差电力系统中测量系统的标准误差 大约大约为正常测量范围的为正常测量范围的 ，误差大于，误差大于 的测量值就为不良数据，但在实用中由于的测量值就为不良数据，但在实用中由于达不到这个标准，所以通常把误差达到达不到这个标准，所以通常把误差达到 以上的数据作为不良数据。当以上的数据作为不良数据。当电力系统出现不良数据时，需要通过检测电力系统出现不良数据时，需要通过检测与辨识的方法处理，以满足状态估计计算与辨识的方法处理，以满足状态估计计算对测量数据的要求。对测量数据的要求。99正正态态分分布布

67、图图100七不良数据的检测与辨识七不良数据的检测与辨识 检测是判定是否存在不良数据，而辨检测是判定是否存在不良数据，而辨识则是为了寻找出哪一个数据是不良数识则是为了寻找出哪一个数据是不良数据，以便进行剔除或补充据，以便进行剔除或补充。 不良数据的出现，会在目标函数不良数据的出现，会在目标函数 中中得到反映，使它大大偏离正常值。因此，得到反映，使它大大偏离正常值。因此，可以根据对可以根据对 的检测来确定不良数据的的检测来确定不良数据的是否存在。是否存在。101七不良数据的检测与辨识七不良数据的检测与辨识 目标函数如式目标函数如式(2-16)(2-16)所示，其中所示，其中 项项可用残差可用残差

68、表示表示 (2-107)(2-107) 测量误差为测量误差为 ，则残差可写成，则残差可写成 ( :( :拗拗) ) (2-108)(2-108)式中：式中： 是是残差灵敏度矩阵残差灵敏度矩阵。式。式(2-108)(2-108)也也就是前述的式就是前述的式(2-25)(2-25)，称为，称为残差方程残差方程，它，它表示了表示了残差与测量误差残差与测量误差间的关系。间的关系。102七不良数据的检测与辨识七不良数据的检测与辨识 下面定义加权残差下面定义加权残差 (2-109)(2-109) 再定义加权测量误差再定义加权测量误差 (2-110)(2-110) 引入上述定义后，残差方程可以写成引入上述定

69、义后，残差方程可以写成 (2-111)(2-111) 式中：为加权残差灵敏度。其表示式为式中：为加权残差灵敏度。其表示式为 (2-112)(2-112)103七不良数据的检测与辨识七不良数据的检测与辨识 采用加权残差灵敏度，从数学运算方面采用加权残差灵敏度，从数学运算方面可以带来一些方便。例如：可以带来一些方便。例如： 是不对称的，是不对称的，而而 是对称的，所以有是对称的，所以有 (2-113)(2-113) (2-114) (2-114)以及加权残差的协方差阵为以及加权残差的协方差阵为 (2-115)(2-115) 不良数据的检测一般不良数据的检测一般通过检查目标函数通过检查目标函数是否远

70、离正常值或残差是否超过正常值来是否远离正常值或残差是否超过正常值来反映反映。104七不良数据的检测与辨识七不良数据的检测与辨识不良数据的检测常用的方法有三种：不良数据的检测常用的方法有三种：n 检测法检测法n加权残差检测法加权残差检测法n标准化残差检测法标准化残差检测法以下分别介绍。以下分别介绍。105 检测法检测法 (1)(1)假定电力系统没有不良数据，加权假定电力系统没有不良数据，加权残差为残差为 ，加权测量误差为，加权测量误差为 。于是目标。于是目标函数为函数为 将式将式(2-111)(2-111)和式和式(2-114)(2-114)代入上式得代入上式得 (2-116)(2-116) 可

71、见可见 为为 的二次型。正常情况下测的二次型。正常情况下测量为正态分布时量为正态分布时 是是 - -分布。其数学期分布。其数学期望和方差可以分别由式望和方差可以分别由式(2-116)(2-116)的展开式的展开式求出求出。七不良数据的检测与辨识七不良数据的检测与辨识106 的数学期望和方差为的数学期望和方差为 式中：式中： 为测量冗余度，即为测量冗余度，即 - -分布的自分布的自由度；由度； 为为 阶自由度的阶自由度的 - -分布的随机分布的随机变量，可写为变量，可写为 (2-117)(2-117) 随着自由度的增大，随着自由度的增大， 越来越接近正态越来越接近正态分布，当分布，当 时，可以用

72、相应的正态分时，可以用相应的正态分布代替布代替 - -分布。分布。七不良数据的检测与辨识七不良数据的检测与辨识107的标准化随机变量形式为的标准化随机变量形式为: :(2-118)(2-118)( )( )七不良数据的检测与辨识七不良数据的检测与辨识108卡方分布图卡方分布图109正正态态分分布布图图110卡方分布图卡方分布图111 (2) (2)假定在电力系统的测量量中，第假定在电力系统的测量量中，第 个个量是值为量是值为 的不良数据，于是测量的误差的不良数据，于是测量的误差向量为：向量为： (2-119)(2-119) 式中：式中： 是是 维向量，其中仅维向量，其中仅 元素为元素为1 1，

73、其余元素均为其余元素均为0 0，即，即 (2-120)(2-120) 此时加权测量误差向量为此时加权测量误差向量为 将上式代入式将上式代入式(2-111)(2-111)得得 于是含一个不良数据时的目标函数将不于是含一个不良数据时的目标函数将不同于式同于式(2-116)(2-116)，而是，而是 (2-(2-121)121)七不良数据的检测与辨识七不良数据的检测与辨识112 式式(2-121)(2-121)右侧第一项即为右侧第一项即为 ，是，是 - -分布；第二项是分布；第二项是0 0均值的正态分布。第三均值的正态分布。第三项为常数。所以项为常数。所以 的数学期望与方差分的数学期望与方差分别为别

74、为 (2-122)(2-122) (2-123) (2-123) 时，式时，式(2-121)(2-121)右侧第一项趋于正右侧第一项趋于正态分布，此时整个态分布，此时整个 也趋于正态分布。也趋于正态分布。的标准化随机变量形式为的标准化随机变量形式为 (2-(2-124)124)七不良数据的检测与辨识七不良数据的检测与辨识113 比较式比较式(2-121)(2-121)可以看出，存在不良数可以看出，存在不良数据后，目标函数据后，目标函数 急剧增大。利用这一急剧增大。利用这一特性可以检测不良数据，具体方法是用特性可以检测不良数据，具体方法是用 和和 两种假设性检验方法。内容如下：两种假设性检验方法

75、。内容如下： 1)1) 假设：如假设：如 ( ( 为检验阈值为检验阈值) )，则没，则没有不良数据，有不良数据， 属真。属真。 2)2) 假设：如假设：如 ( ( 为检验阈值为检验阈值) )，则有，则有不良数据，不良数据， 属真。属真。七不良数据的检测与辨识七不良数据的检测与辨识114 当确定了阈值当确定了阈值 后，如某次采样后，如某次采样 就就认为认为 属真。这时可能犯属真。这时可能犯第一类错误第一类错误，即，即 属真而拒绝了属真而拒绝了 ，接受了，接受了 。这类错误。这类错误称误报警，其出现的概率为称误报警，其出现的概率为 称为称为伪警概伪警概率率。上述检验结果也可能犯。上述检验结果也可能

76、犯第二类错误第二类错误，即即 不真而接受了不真而接受了 ，拒绝了，拒绝了 。这类错。这类错误称漏报，出现的概率为误称漏报，出现的概率为 称为称为漏检概率漏检概率。这两类错误的概率由阈值这两类错误的概率由阈值 确定，一般漏确定，一般漏检概率越小，伪警概率就越大，反之亦然。检概率越小，伪警概率就越大，反之亦然。七不良数据的检测与辨识七不良数据的检测与辨识115 为了减少这两类错误，通常将概率范为了减少这两类错误，通常将概率范围取为围取为 。若。若 且且 ，则可由给定的则可由给定的 正态分布表查到相应正态分布表查到相应的的 值为值为1.645 1.645 。七不良数据的检测与辨识七不良数据的检测与辨

77、识116七不良数据的检测与辨识七不良数据的检测与辨识117 加权残差检测法加权残差检测法 由残差定义可知由残差定义可知, ,残差也是一个按正态残差也是一个按正态分布的随机变量。分布的随机变量。又由于加权残差的权值又由于加权残差的权值是相应测量量标准差的倒数，因而是相应测量量标准差的倒数，因而加权残加权残差也符合正态分布。差也符合正态分布。所以利用加权残差同所以利用加权残差同样也可以用假设性检验的方法来检测不良样也可以用假设性检验的方法来检测不良数据。数据。七不良数据的检测与辨识七不良数据的检测与辨识118 由于是正态分布，故由于是正态分布，故 ， 的对角元的对角元素就是加权残差的方差素就是加权

78、残差的方差 即即 为正态分布随机变量为正态分布随机变量 通常测量情况下，若规定伪警概率为通常测量情况下，若规定伪警概率为 ，则正常的加权残差取值范围为，则正常的加权残差取值范围为 (2-125)(2-125) 于是加权残差阈值可定为于是加权残差阈值可定为 (2-126)(2-126)七不良数据的检测与辨识七不良数据的检测与辨识119 加权残差加权残差 检测是将逐维残差按假设性检测是将逐维残差按假设性检验的方法来进行。检验的方法来进行。 (2-127)(2-127) 式中：式中：七不良数据的检测与辨识七不良数据的检测与辨识120 标准化残差检测法标准化残差检测法 除了除了 检测法与加权残差检测法

79、外，检测法与加权残差检测法外，有时可以采用标准化残差检测方法以取得有时可以采用标准化残差检测方法以取得更理想的效果。更理想的效果。 标准化残差的定义为标准化残差的定义为 (2-128)(2-128) 其中其中 （p112p112）(2-129)(2-129) 于是于是 (2-130)(2-130) 式中：式中： 是矩阵是矩阵 的第的第 个对角元素。个对角元素。七不良数据的检测与辨识七不良数据的检测与辨识121 由残差方程式由残差方程式(2-110)(2-110)可以写出可以写出标准化标准化残差方程式残差方程式为为 (2-131)(2-131) 式中：式中： 为为标准化残差灵敏度矩阵标准化残差灵

80、敏度矩阵。 在正常测量条件下，在正常测量条件下， 具有下列关系具有下列关系 (2-132)(2-132) 将式将式(2-129)(2-129)代入上式，可得上式右端代入上式，可得上式右端矩阵的对角元素均为矩阵的对角元素均为1 1 ，故有，故有 (2-133)(2-133)七不良数据的检测与辨识七不良数据的检测与辨识122 当当 时，得到第时，得到第 个标准化残差个标准化残差的检测阈值为的检测阈值为 (2-134)(2-134) 逐维残差的标准化残差检测方法为逐维残差的标准化残差检测方法为 (2-135)(2-135) 式中：式中： 为第为第 个标准化残差分量。个标准化残差分量。七不良数据的检测

81、与辨识七不良数据的检测与辨识123 以上三种检测方法的共同特点是利用以上三种检测方法的共同特点是利用采样的残差信息来检测不良数据，其检测采样的残差信息来检测不良数据，其检测效果与阈值的选择有关，当阈值较低时，效果与阈值的选择有关，当阈值较低时，检测不良数据的能力较强，但是过低的阈检测不良数据的能力较强，但是过低的阈值又会使误检率增大。值又会使误检率增大。 检测法是一种检测法是一种总体型的检测总体型的检测，它能，它能测知不良数据是否存在，但不能知道哪一测知不良数据是否存在，但不能知道哪一个是不良数据。在系统个是不良数据。在系统规模较大及冗余度规模较大及冗余度大大的情况下，个别不良数据对的情况下，

82、个别不良数据对 的影响的影响相对减小，亦即式相对减小，亦即式(2-121)(2-121)右侧的第三顶右侧的第三顶相对减小，从而使检测的灵敏度较低。相对减小，从而使检测的灵敏度较低。七不良数据的检测与辨识七不良数据的检测与辨识124 与与 检测法与系统大小无关，它取决检测法与系统大小无关，它取决于于 或或 的对角元素。当测量系统完善，的对角元素。当测量系统完善，冗余度冗余度 越大，则对角元素越占优势，检越大，则对角元素越占优势，检测不良数据越灵敏。在冗余度为测不良数据越灵敏。在冗余度为 时，时， 法比法比 法在灵敏度方面更优越，但法在灵敏度方面更优越，但是是 法需付出计算法需付出计算 的代价，在

83、冗余度更的代价，在冗余度更高时这两种方法的效果相近。高时这两种方法的效果相近。 与与 法在法在单个不良数据单个不良数据时一般可取得时一般可取得理想的效果，但理想的效果，但有时除了不良数据点的残有时除了不良数据点的残差超过检测阈值外，一些正常测点的残差差超过检测阈值外，一些正常测点的残差也超过阈值，这种现象称为也超过阈值，这种现象称为残差污染残差污染。七不良数据的检测与辨识七不良数据的检测与辨识125 有多个不良数据时有多个不良数据时，由于相互作用可能由于相互作用可能导致部分或全部不良数据测点上的残差近导致部分或全部不良数据测点上的残差近于正常残差现象，这称为于正常残差现象，这称为残差淹没残差淹

84、没。残差。残差污染和残差淹没使不良数据点模糊，导致污染和残差淹没使不良数据点模糊，导致辨识不良数据的困难。辨识不良数据的困难。 在应用在应用 或或 检验时，增加测量可使检验时，增加测量可使 矩阵的对角元素增大，同时使其非对角元矩阵的对角元素增大，同时使其非对角元素减小。素减小。七不良数据的检测与辨识七不良数据的检测与辨识126不良数据的辨识法不良数据的辨识法七不良数据的检测与辨识七不良数据的检测与辨识残差搜索辨识法残差搜索辨识法估计辨识法估计辨识法 加权残差搜索法加权残差搜索法标准化残差搜索法标准化残差搜索法 128 二二 不良数据的残差搜索辨识法不良数据的残差搜索辨识法 对不良数据辨识的基本

85、思路是：检测出对不良数据辨识的基本思路是：检测出不良数据后，设法找出这个不良数据并不良数据后，设法找出这个不良数据并在测量向量中将其排除，然后重新进行在测量向量中将其排除，然后重新进行状态估计。状态估计。 假设在检测中发现有不良数据。最简单假设在检测中发现有不良数据。最简单的辨识方法是在的辨识方法是在 个测量量中去掉第一个个测量量中去掉第一个测量量，再用不良数据检测法检查余下测量量，再用不良数据检测法检查余下的的 个中是否存在不良数据。如果个中是否存在不良数据。如果 个测量的个测量的 值与值与 个时的个时的 值差不多，值差不多，则表示第一个测量量是正常量，应予以则表示第一个测量量是正常量，应予

86、以恢复。恢复。七不良数据的检测与辨识七不良数据的检测与辨识129 然后试第二个测量量，直到找出不良然后试第二个测量量，直到找出不良数据为止。如果存在两个不良数据，则数据为止。如果存在两个不良数据，则应试探每次去掉两个测量量的各种组合。应试探每次去掉两个测量量的各种组合。这种方法试探的次数非常多这种方法试探的次数非常多, ,而且每次试而且每次试探都要进行状态估计，因此问题的关键探都要进行状态估计，因此问题的关键在于如何减少试探的次数。在于如何减少试探的次数。 残差搜索辨识法，即用残差绝对值由大残差搜索辨识法，即用残差绝对值由大到小排队来逐维作试探，通常分为到小排队来逐维作试探，通常分为 与与 。

87、七不良数据的检测与辨识七不良数据的检测与辨识130 加权残差搜索法加权残差搜索法 加权残差搜索法是按加权残差搜索法是按 大小排队，逐大小排队，逐维试探。加权残差可以写成维试探。加权残差可以写成 (2-136)(2-136) 式中：式中： 为有不良数据为有不良数据 时时 维加权残差维加权残差向量；向量； 为没有不良数据时的为没有不良数据时的 维加权残维加权残差向量；差向量； 为为 矩阵的第矩阵的第 个列向量；为个列向量；为出现在第出现在第 点上的不良数据值。点上的不良数据值。七不良数据的检测与辨识七不良数据的检测与辨识131 略去正常残差，式略去正常残差，式(2-136)(2-136)可以写成可

88、以写成 (2-137)(2-137) 由于需研究的是由于需研究的是 的大小与其排队次序的大小与其排队次序问题。由式问题。由式(2-137)(2-137)可见，与可见，与 的值是影的值是影响因素之一，此外响因素之一，此外 因此，排队顺序可用上列系数值分别除因此，排队顺序可用上列系数值分别除以以 、 来确定。来确定。七不良数据的检测与辨识七不良数据的检测与辨识132 当当 比比 小，则小，则 法排队次序的提前较法排队次序的提前较 法更明显。对于注入功率较小，而穿越法更明显。对于注入功率较小，而穿越功率较大的节点，由于功率较大的节点，由于 与与 分别与注分别与注入功率与穿越功率成线性关系，所以用入功

89、率与穿越功率成线性关系，所以用法排队效果较好。法排队效果较好。七不良数据的检测与辨识七不良数据的检测与辨识133 标准化残差搜索法标准化残差搜索法 标准化残差搜索法是按标准化残差搜索法是按 排队，逐维排队，逐维试探。标准化残差可以写成试探。标准化残差可以写成 (2-139)(2-139) 式中：为有一个不良数据式中：为有一个不良数据 时的时的 维标准维标准化残差向量，化残差向量， 为正常测量时的为正常测量时的 维标准维标准化残差向量；化残差向量； 的含义见式的含义见式(2-129)(2-129)；为；为 阵的第阵的第 个列向量。个列向量。 略去正常残差略去正常残差 ，式，式(2-139)(2-

90、139)可写成可写成 (2-140)(2-140)七不良数据的检测与辨识七不良数据的检测与辨识134 由式由式(2-132)(2-132)和式和式(2-133)(2-133)得得 (2-141)(2-141) 这表示随机变量这表示随机变量 的自相关系数为的自相关系数为1 1。根据概率论，在式根据概率论，在式(2-132)(2-132)中，、中，、 之之间的互相关系数的绝对值恒小于或等于间的互相关系数的绝对值恒小于或等于1 1。 (2-142)(2-142) 与式与式(2-133)(2-133)比较，可得比较，可得 (2-143)(2-143)七不良数据的检测与辨识七不良数据的检测与辨识135

91、可见，在单个不良数据时，按可见，在单个不良数据时，按 大小排大小排队，不良数据点的标准化残差绝对值总是队，不良数据点的标准化残差绝对值总是排在前面。计及正常残差排在前面。计及正常残差 的影响后也排的影响后也排在前面，亦即按在前面，亦即按 法只需搜索法只需搜索1 12 2次即可次即可辨识成功。辨识成功。 残差搜索法只适用于单个不良数据的辨残差搜索法只适用于单个不良数据的辨识，或弱相关的多个不良数据的辨识。对识，或弱相关的多个不良数据的辨识。对于强相关的多个不良数据，则由于需搜索于强相关的多个不良数据，则由于需搜索次数过多而难以奏效。次数过多而难以奏效。七不良数据的检测与辨识七不良数据的检测与辨识

92、136 残差搜索法在确定一个残差大的可疑数残差搜索法在确定一个残差大的可疑数据并将它暂时排除后，需重作状态估计以据并将它暂时排除后，需重作状态估计以确定排除的是否真为不良数据。因此，需确定排除的是否真为不良数据。因此，需进行多次状态估计，在大系统中会耗时过进行多次状态估计，在大系统中会耗时过多，但这种方法程序简单，占用内存少，多，但这种方法程序简单，占用内存少，对状态估计程序的适应性好。对状态估计程序的适应性好。 与与 辨识法可以用相同的程序流程框辨识法可以用相同的程序流程框图，如图图，如图2-142-14所示。图中程序说明如下。所示。图中程序说明如下。七不良数据的检测与辨识七不良数据的检测与

93、辨识137 框框1 1 进行一次采样，形成测点集合进行一次采样，形成测点集合 ，在，在 中可疑数据集合为中可疑数据集合为 ，不良数据集合为，不良数据集合为 。 框框2 2 用用 进行状态估计。计算出进行状态估计。计算出 或或 和和 。 框框3 3 进行进行 检测，若检测，若 无不良数据，无不良数据，输出状态估计的各项数据，程序转到出口。输出状态估计的各项数据，程序转到出口。若若 ，进入框，进入框4 4 。 框框4 4 将现有将现有 的保留在的保留在 中。中。 框框5 5 将将 中各测点按中各测点按 或或 的大小排队。的大小排队。 框框6 6 将将 中排在前面的测点号中排在前面的测点号 送到送到

94、 ，并将并将 从从 清除。清除。七不良数据的检测与辨识七不良数据的检测与辨识138 框框7 7 在总的测点集合在总的测点集合 中扣除中扣除 的测量，的测量，形成新的形成新的 集合。集合。 框框8 8 用新的测点集合用新的测点集合 作状态估计，计算作状态估计，计算 或或 及及 。 框框9 9 进行进行 检测，若已无不良数据，则程检测，若已无不良数据，则程序转出口，否则到框序转出口，否则到框10 10 。 框框10 10 比较比较 与与 ，若有，若有 显著减小，显著减小，则表示则表示 点是不良数据，转框点是不良数据，转框4 4，再对，再对 集合用框集合用框8 8的的 或或 重新排队；若重新排队；若

95、 无明无明显减小。则显减小。则 不是不良数据，予以恢复，不是不良数据，予以恢复，转框转框11 11 。七不良数据的检测与辨识七不良数据的检测与辨识139 框框1111将将 从从 中清除，不必重新排队，只需取中清除，不必重新排队，只需取框框6 6中排在第二位的测点进行试探。中排在第二位的测点进行试探。 在上述不良数据的残差搜索辨识中，若不在上述不良数据的残差搜索辨识中，若不将残差绝对值大的测量量从状态估计中排除，将残差绝对值大的测量量从状态估计中排除，而是在迭代过程中减小它的权值，即减小它在而是在迭代过程中减小它的权值，即减小它在状态估计中的影响，使得最终能获得最精确的状态估计中的影响，使得最终

96、能获得最精确的状态估计量。这种修改加权最小二乘目标函数状态估计量。这种修改加权最小二乘目标函数的方法在辨识单个不良数据或多个不产生残差的方法在辨识单个不良数据或多个不产生残差淹没的不良数据时是非常有效的，但在出现残淹没的不良数据时是非常有效的，但在出现残差淹没时则难以奏效，而且它对测量系统的要差淹没时则难以奏效，而且它对测量系统的要求较高。求较高。七不良数据的检测与辨识七不良数据的检测与辨识140 三三 不良数据的估计辨识法不良数据的估计辨识法 不良数据的估计辨识是我国发展起来不良数据的估计辨识是我国发展起来的一种新方法，这种方法具有较好的辨识的一种新方法，这种方法具有较好的辨识多个不良数据的

97、功能，实时性也较好，因多个不良数据的功能，实时性也较好，因此近年来已得到了进一步的发展。此近年来已得到了进一步的发展。 设系统有设系统有 个状态量，个状态量， 个测量量，多个测量量，多余测量信息余测量信息 。根据可观察性的概念，。根据可观察性的概念，如果有如果有 个不良数据从测量量中移去，余个不良数据从测量量中移去，余下的测量量不能保证系统的可观察性，即下的测量量不能保证系统的可观察性，即不能作出不能作出 的估计，则这的估计，则这 个量就不可能个量就不可能辨识。辨识。七不良数据的检测与辨识七不良数据的检测与辨识142 利用各开关的遥信及线路潮流的遥测利用各开关的遥信及线路潮流的遥测检查其是否对

98、应来发现可疑的遥信。然检查其是否对应来发现可疑的遥信。然后用搜索法根据可能存在的网络结构形后用搜索法根据可能存在的网络结构形式，通过状态估计检查式，通过状态估计检查 是否小于门槛是否小于门槛值来判别该类结构形式是否正确。值来判别该类结构形式是否正确。 八八. .电力系统网络拓扑分析及网络结构辨识的基本概念电力系统网络拓扑分析及网络结构辨识的基本概念143 表表 2-2 2-2 结构错误检测结果分类表结构错误检测结果分类表八八. .电力系统网络拓扑分析及网络结构辨识的基本概念电力系统网络拓扑分析及网络结构辨识的基本概念类别类别SWITSWITMEASMEAS检测结果检测结果正常合正常合0 0正常

99、正常误识别开断误识别开断1 1可疑可疑误识别闭合误识别闭合0 00 0可疑可疑正常开断正常开断1 10 0正常正常144 利用上述网络结构错误的检测结果，可利用上述网络结构错误的检测结果，可以确定可能有几种正确的网络结构，并以以确定可能有几种正确的网络结构，并以此网络结构为准进行此网络结构为准进行 检验。然后以检验检验。然后以检验结果来辨识正确的结构方式。结果来辨识正确的结构方式。 若若 (2-156)(2-156) 则网络结构错误已被测出，并得到了正确则网络结构错误已被测出，并得到了正确的结构方式。的结构方式。 若若 (2-157)(2-157) 则表示仍需对另一个可能的网络结构进行则表示仍

100、需对另一个可能的网络结构进行估计辨识。直至找到正确的网络结构或几估计辨识。直至找到正确的网络结构或几种可能的网络结构全部被辨识完为止。种可能的网络结构全部被辨识完为止。八八. .电力系统网络拓扑分析及网络结构辨识的基本概念电力系统网络拓扑分析及网络结构辨识的基本概念145 这种方法的优点是比较简单，对测量这种方法的优点是比较简单，对测量量配置要求也较低。但在多个结构错误量配置要求也较低。但在多个结构错误与不良数据同时出现时可能出现失败。与不良数据同时出现时可能出现失败。因此，在一些简单电力系统中，因此，在一些简单电力系统中，网络结网络结构的辨识通常是仅比较遥测信号是否一构的辨识通常是仅比较遥测

101、信号是否一致单独完成的。致单独完成的。八八. .电力系统网络拓扑分析及网络结构辨识的基本概念电力系统网络拓扑分析及网络结构辨识的基本概念146 电力系统状态估计理论，在电力系统状态估计理论，在2020世纪世纪7070年代初期就已经确立。但实际应用则是在年代初期就已经确立。但实际应用则是在近二十几年中才不断地取得完善，其中以近二十几年中才不断地取得完善，其中以最小二乘法应用最为广泛，并且为了进一最小二乘法应用最为广泛，并且为了进一步改善其功效，陆续又出现了解耦法与正步改善其功效，陆续又出现了解耦法与正交变换法。交变换法。 1 1）目前最小二乘估计存在的主要问题之）目前最小二乘估计存在的主要问题之

102、一是，在某些情况下可能出现病态，导致一是，在某些情况下可能出现病态，导致难以求解。难以求解。状态估计小结状态估计小结147 2 2）本章的解算方法，是假定测量信息是本章的解算方法，是假定测量信息是同时得到的同时得到的，而实际上数据采集系统各个，而实际上数据采集系统各个数据的时刻是不同的，这种情况对大多数数据的时刻是不同的，这种情况对大多数状态估计影响不大，仍可满足要求。但对状态估计影响不大，仍可满足要求。但对不是取自一个数据采集系统的数据不是取自一个数据采集系统的数据, ,时间偏时间偏差问题就会非常严重。目前此类问题尚未差问题就会非常严重。目前此类问题尚未很好解决。很好解决。状态估计小结状态估

103、计小结148 3 3）电力系统网络结构和测量配置，是）电力系统网络结构和测量配置，是实现状态估计的重要条件实现状态估计的重要条件。当测量不足时。当测量不足时系统是不可观察的，无法进行状态估计，系统是不可观察的，无法进行状态估计，所以应在状态估计之前先进行可观察检验。所以应在状态估计之前先进行可观察检验。状态估计小结状态估计小结149 4 4）在不良数据识别与辨识方面）在不良数据识别与辨识方面，对于，对于仅存在单个不良数据的情况，用假设性检仅存在单个不良数据的情况，用假设性检验法检查是否存在不良数据，并以搜索法验法检查是否存在不良数据，并以搜索法进行辨识一般是有效的。进行辨识一般是有效的。状态估

104、计小结状态估计小结150 5) 5)当系统中存在关联的多个不良数据时，当系统中存在关联的多个不良数据时，最大的测量残差不一定对应于有大误差的最大的测量残差不一定对应于有大误差的测量。小测量残差的量也可能是不良数据，测量。小测量残差的量也可能是不良数据，问题十分复杂。不良数据检测和辨识的方问题十分复杂。不良数据检测和辨识的方法还在继续不断地得到开发。法还在继续不断地得到开发。 状态估计小结状态估计小结151电力工程系电力工程系Department of Electrical EngineeringNorth China Electric Power UniversityThanks Http ee152