九、电力系统状态估计(二).ppt

《九、电力系统状态估计(二).ppt》由会员分享，可在线阅读，更多相关《九、电力系统状态估计(二).ppt（51页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、东南大学电气工程系东南大学电气工程系九电力系统状态估计（二）版权所有东南大学电气工程系东南大学电气工程系内容提要内容提要n基本概念n不良数据的检测n 检测n加权残差检测n标准化残差检测n算例n不良数据辨识n残差搜索辨识n零残差辨识n估计辨识版权所有东南大学电气工程系东南大学电气工程系 一、基本概念一、基本概念一、基本概念一、基本概念n不良数据是指误差大于某一标准（如从统计学角度，大于3-10倍标准方差）的量测数据。n只有排除不良数据才能得到正确的状态估计结果，这一过程称为不良数据检测与辨识过程。n对SCADA原始量测数据的状态估计结果进行检查，判断是否存在不良数据并指出具体可疑量测数据的过程称

2、之为不良数据检测。n对检测出的可疑数据验证真正不良数据的过程称之为不良数据的辨识。版权所有东南大学电气工程系东南大学电气工程系n状态估计辨识不良数据的能力来自于量测系统的冗余度n能够估计出全部状态量的量测系统具有可观测性n而去掉不良数据仍保持可观测性的量测系统具有可辨识性。版权所有东南大学电气工程系东南大学电气工程系残差方程残差方程n假设按最小二乘估计算法获得状态估计 后，则量测估计 为：式中，是m维量测估计向量。n将上式在状态真值x附近线形化：式中：h(x)是m维量测真值向量 是n维状态估计误差向量；H是mn阶雅可比矩阵。(1)版权所有东南大学电气工程系东南大学电气工程系n状态估计误差：n量

3、测估计误差方程阵：(4)(5)n量测估计误差：(3)n对线性情况，状态估计值为(2)版权所有东南大学电气工程系东南大学电气工程系n定义：残差r为量测向量与量测估计向量之差 r n将量测方程式和（4）式代入（6）式，可得：以及 W为mm阶残差灵敏度矩阵I为单位矩阵(6)(7)(8)版权所有东南大学电气工程系东南大学电气工程系W为mm阶残差灵敏度矩阵(7)式就是残差方程，它建立了残差与量测误差之间的关系 n残差r的协方差矩阵为(7)(9)n量测残差：r为量测向量与量测估计向量之差 n量测误差：v为量测向量与量测真实值向量之差 版权所有东南大学电气工程系东南大学电气工程系WW特点特点nW阵是幂等阵，

4、W2 WnW阵秩为K（K=mn），由于其秩小于其维数m，无法求逆。nW矩阵不对称，但是WR是对称的。n残差灵敏度矩阵W与电网结构参数以及量测分布有直接关系n一般而言残差灵敏度矩阵具有对角占优特性，有利于可疑数据检测，因为对应于不良数据的残差项最大，即根据残差大小就可以检测不良数据。版权所有东南大学电气工程系东南大学电气工程系小结小结观测值观测值真实值真实值估计值估计值估计误差估计误差残差残差误差误差状态-x-量测zh(x)v=z-h(x)r 版权所有东南大学电气工程系东南大学电气工程系加权残差加权残差加权残差加权残差 n将残差方程用加权残差表示为：n其中：rw为m维加权残差向量，Ww为mm阶加

5、权残差灵敏度矩阵，vw为m维加权量测误差向量。且分别有：n加权残差灵敏度矩阵具有n幂等性 Ww2=Wwn对称性 WwT=Ww (10)(11)(12)(13)n加权残差rw的协方差矩阵为：版权所有东南大学电气工程系东南大学电气工程系标准化残差标准化残差标准化残差标准化残差 n记对角阵D为：DdiagWR=diagr n标准化残差的定义如下：n或者有：其中，是矩阵r的第i个对角元素。n定义标准化残差灵敏度矩阵WN为 n则标准化残差方程：(14)(15)(16)(17)(18)版权所有东南大学电气工程系东南大学电气工程系二、不良数据的检测二、不良数据的检测二、不良数据的检测二、不良数据的检测n不良

6、数据的检测一般是通过检查目标函数是否大大偏离正常值或残差是否超过正常值来反映的，常用有三种方法：n 检测n加权残差检测n标准化残差检测 版权所有东南大学电气工程系东南大学电气工程系n当正常情况下测量为正态分布时，由概率论可知，就是c2分布。其数学期望和方差可以分别由(19)式的展开式求出。n由此可见，为vwz的二次型。检测检测n(1)先假定电力系统没有不良数据存在，此时加权残差写为rwz，加权测量误差写为vwz，于是目标函数为由式（10）得：(19)版权所有东南大学电气工程系东南大学电气工程系n求 数学期望n方差 检测检测推得：n根据量测误差矢量vz为0均值的正态分布及其分量间两两独立的假设。

7、并考虑到（13）式。版权所有东南大学电气工程系东南大学电气工程系n随着自由度的增加，c2(K)越来越接近于正态分布，当K30时，可以用相应的正态分布来代替c2分布，的标准化随机变量形式为：检测检测n式中：K为测量的冗余度，亦即是c2分布的自由度；为K阶自由度的c2分布的随机变量，可写为 c2（K）N(0，1)，当K30(20)(21)版权所有东南大学电气工程系东南大学电气工程系将上式代入n此时加权测量误差向量可写为 n(2)假定在电力系统的测量量中，第i个量是不良数据，且假定其值为ai，于是测量的误差向量应改写为式中：ei是m维列向量，其中仅i元素为1，其余元素均为0，即(22)(23)检测检

8、测版权所有东南大学电气工程系东南大学电气工程系n上式右侧的第一项即为 ，是c2分布；第二项是0均值分布；第三项为常数，所以 的数学期望与方差分别为：检测检测n于是含一个不良数据的目标函数为(24)n当K30时，式24右侧第一项趋于正态分布，此时整个 也趋于正态分布。用标准化随机变量形式表示为：(25)(26)版权所有东南大学电气工程系东南大学电气工程系 检测的检测的 H H0 0和和H H1 1两种假设两种假设1.H0假设：如果 (为检验阈值)，则无不良数据，H0属真2.H1假设：如果 (为检验阈值)，则有不良数据，H1属真n由式（24）可以看出，当存在不良数据后，目标函数就急剧增大，利用这一

9、特性可以检验不良数据，其具体方法是用H0和H1两种假设的假设性检验方法，内容如下：版权所有东南大学电气工程系东南大学电气工程系H0：，则无不良数据，H0属真H1：，则有不良数据，H1属真 检测的检测的 H H0 0和和H H1 1两种假设两种假设n当确定了阈值后，如果H0属真而拒绝H0，接受H1，则是误报警。其出现的概率是pe，又称伪警概率。n如果H0不真而接受了H0，拒绝H1，则是漏报警。其出现的概率是pd，又称漏检概率。n例如，若pe0.05，即：如果K30，则可以由给定的N（0，1）正态分布表查出相应的值为1.645。版权所有东南大学电气工程系东南大学电气工程系 加权残差检测加权残差检测

10、n由残差定义可知，当测量值z是符合正态分布的随机变量时，其估计值可以认为等于均值，所以残差也是一个按正态分布的随机变量。n又由于加权残差的权值是相应测量的标准差的倒数，因而加权残差也符合正态分布n所以利用加权残差同样可以用假设性检验的方法来检测不良数据版权所有东南大学电气工程系东南大学电气工程系 加权残差检测加权残差检测n由于是正态分布，故Erwz,i0，那么Ww的对角元素就是加权残差的方差n亦即加权残差rwz,i为下列形式的正态分布随机变量：(27)于是加权残差阀值可定为:(28)n在通常情况下，如果规定Pe0.005，则正常的加权残差取值范围：版权所有东南大学电气工程系东南大学电气工程系

11、加权残差检测加权残差检测n加权残差rw检测是逐维残差按假设性检验的方法来进行的。此时H0属真，接受H0 此时H0不真，接受H1(29)版权所有东南大学电气工程系东南大学电气工程系n将式 代入上式。可得上式右端矩阵的对角元素均为1，故有：标准化残差检测标准化残差检测n在正常测量条件下，具有下列关系：(30)(31)n当Pe0.005时，得到第i个标准化残差的检测阀值：(32)版权所有东南大学电气工程系东南大学电气工程系 标准化残差检测标准化残差检测n逐维残差的标准化检测方法是：此时H0属真，接受H0 此时H0不真，接受H1(33)版权所有东南大学电气工程系东南大学电气工程系几种检测方法比较几种检

12、测方法比较n以上三种检测方法的共同特点是利用采样的残差信息来检测出不良数据，其检测的效果与阀值的选择有关。n 检测法是一种总体型的检测，它能知道不良数据是否存在，但不能知道哪个是不良数据。在系统规模较大及冗余度大的情况下，个别不良数据对 的影响相对减小，亦即式（24）右侧的第三项相对减小，从而使灵敏度降低。版权所有东南大学电气工程系东南大学电气工程系几种检测方法比较几种检测方法比较nrW与rN检测法与系统大小无关，它取决于WW或WN的对角元素。当测量系统完善，冗余度K越大，则对角元素越占优势，检测不良数据越灵敏。n在冗余度为m/n23时，rN法要比rW法在灵敏度方面更优越，但是rN法需付出计算

13、D的代价n在冗余度更高时这两种方法的效果相近。版权所有东南大学电气工程系东南大学电气工程系几种检测方法比较几种检测方法比较nrW与rN法有时除了不良数据点的残差呈现超过检测阀值外，还有一些正常测点的残差也超过阀值，这种现象称为残差污染。n在多个不良数据情况下，由于相互作用会导致部分或全部不良数据测点上的残差近于正常残差的现象，这称为残差淹没。版权所有东南大学电气工程系东南大学电气工程系几种检测方法比较几种检测方法比较n在应用rW与rN检验时，增加测量可使W矩阵的对角元素增大，同时使非对角元素减小。n采用伪测量可以加强残差灵敏度矩阵的对角元素优势，可以有效地削弱单个不良数据情况下的残差污染和多个

14、不良数据情况下的残差淹没现象。n为了减少由于增加伪测量所导致的计算时间增长，应只在真正薄弱的某些测点增加相应的局部测量量。版权所有东南大学电气工程系东南大学电气工程系 例题例题n如图所示的三母线电力系统，支路电抗和节点注入有功功率如图所示。n选取P1、P2、P12、P13、P23作为用于状态估计的量测量估计计算公式版权所有东南大学电气工程系东南大学电气工程系例题例题n则其加权残差平方和 服从 分布，在本例中其自由度k=3，取 从 分布表中查得 ，加权平方和 的估计值计算如下：n假定测量得到的量测量向量z=-1.98,0.502,-0.596-1.404,-0.097T则计算状态量估计值及量测估

15、计值。从而残差向量为：估计结果可以接受 版权所有东南大学电气工程系东南大学电气工程系例题例题n现假设另一组量测量数据z=-1.98,0.502,-0.1-1.404,-0.097T，重新进行上述计算：状态量估计值：量测量估计值：残差向量：n加权残差平方和 21.16 n此值大于11.35，因此，在99%的置信度下原始量测数据中可能存在着坏数据。估计结果不能接受版权所有东南大学电气工程系东南大学电气工程系例题例题n现在我们要计算标准化残差，并将标准化残差最大的那个量测量判为坏数据，将其滤除，再重新进行估计。标准化残差：其中D=diag(I-HG-1HTR-1)R 本例标准化残差：n第三个量测量的

16、标准化残差的绝对值较大，有可能为坏数据。现将其滤除后重新进行估计 版权所有东南大学电气工程系东南大学电气工程系例题例题n去除第三个量测后，量测量数据z=-1.98,0.502,-1.404,-0.097T，重新进行上述计算：状态量估计值：量测量估计值：残差向量：n新自由度k=2，仍取 从 分布表中查得 ，新的加权平方和 的估计值计算结果为0.0245，远小于9.21，可见将原先的第三个量测量确为坏数据，将其滤除后重新进行估计得到了较为可信的结果。版权所有东南大学电气工程系东南大学电气工程系三、不良数据的辨识三、不良数据的辨识n不良数据的辩识方法较多，主要有：n残差搜索辩识法；n零残差辩识法；n

17、估计辩识法；版权所有东南大学电气工程系东南大学电气工程系 残差搜索辨识法残差搜索辨识法n过程：进行估计计算，量测按残差大小排队，去掉残差最大的量测。不断重复上述过程直至所有残差大的量测。n通常分rW 和rN法n残差搜索法一般只适用于单个不良数据的辨识，或弱相关的多个不良数据的辨识。对于强相关的可疑数据的情况，则由于需搜索次数过多而难以有效。n残差搜索法在确定一个残差大 的可疑数据并将它暂时排除后，需重作状态估计以确定排除的是否真为不良数据。n程序简单，但要进行多次状态估计，需时多版权所有东南大学电气工程系东南大学电气工程系 残残差差搜搜索索法法流流程程框框图图版权所有东南大学电气工程系东南大学

18、电气工程系 零残差辩识法零残差辩识法n不将残差绝对值大的测量量从状态估计中排除，而是在迭代过程中减小它的权值，亦即减小它在状态估计中的影响，使得最终能获得精确的状态估计量。n零残差辩识法：n将残差大的量测权重降低为0，参与到估计计算公式中。(34)其中的S为mm阶对角矩阵，对应于不良数据的对角元素为0，其余为1版权所有东南大学电气工程系东南大学电气工程系 零残差辩识法特点零残差辩识法特点n不需要重新设定量测模型进行状态估计，计算速度快n与估计辨识法比，不需要存放残差灵敏度矩阵，节省内存n程序简单，能配合各种状态估计算法，适应性强n缺点：量测系统弱，不良数据过多或过大，收敛速度下降。版权所有东南

19、大学电气工程系东南大学电气工程系 估计辩识法基础估计辩识法基础n估计辩识法的基础是残差方程(7)：n重要特点及论断：nW矩阵无逆，秩等于K=m-n，即残差方程式中只有K个独立线性方程，相应地，只可以估计出K个未知数。n假定在误差矢量v中有某K个分量，把它做为状态量，就可能利用量测多余信息把这K个误差分量估计出来。n如果不良数据超过K个，就不可能被辨识。因为这意味着用K个独立方程去解超过K个的独立变量。n即使不良数据个数P小于K个，也不一定可以辨识。要看P的测点组合。版权所有东南大学电气工程系东南大学电气工程系 估计辩识法（估计辩识法（1 1）n根据可观察性的概念，如果有p个不良数据从m个测量量

20、中移去，余下的量测量不能保证系统的可观察性，则这p个量也就不可能辨识。n如果在一次测量中可疑数据有s个，可靠数据有t个，测量总数为mts，则p个不良数据也包括在可疑数据s中，于是有：Ws、Wt分别为W阵中对应于数据s和t部分的ms及mt子阵；(35)vs、vt分别为可疑数据误差向量与可靠数据误差向量。版权所有东南大学电气工程系东南大学电气工程系 估计辩识法（估计辩识法（2 2）n对于t维正常误差向量vt有：(36)式中：Rt为tt阶对角阵n在可疑数据误差向量vs中可能含有正常测量误差，但肯定含有p个不良数据，可以建立一个以可疑数据误差向量为变量的目标函数。版权所有东南大学电气工程系东南大学电气

21、工程系 估计辩识法（估计辩识法（3 3）n建立以可疑数据误差向量为变量的目标函数。：式中：G-1为mm阶的正定加权阵，可以取G为 的对角阵，其中的元素为：(37)(38)n式中：wt,ij为矩阵Wt中第i行第j列元素；t,j2为v t,j的方差nvs的加权最小二乘目标可通过对（37）式的导数等于零来求得版权所有东南大学电气工程系东南大学电气工程系 估计辩识法（估计辩识法（4 4）n利用对（37）式的导数等于0，得：n式中：是可疑数据误差的估计值，它可以用来判断那些分量为不良数据，哪些数据为正常测量误差。从而达到辨识不良数据的目的。上式也可称为可疑数据的估计方程。得(39)n在求出可疑数据误差的

22、估计值后，可以求出状态估计修正量。版权所有东南大学电气工程系东南大学电气工程系 估计辩识法（估计辩识法（5 5）n状态估计的修正量计算方法如下：n 正常情况下状态估计误差表达式如下：n当存在s个不良数据时，其测量误差向量为vb，这时状态估计误差表达式为 式中：为含有不良数据时求出的状态估计值。n将上述两式相减，可得(40)(41)(42)n式中：vsm为一个m维向量，其中对应于s个不良数据测点的相应元素等于不良数据的真值vsi，而其余元素均为零。版权所有东南大学电气工程系东南大学电气工程系 估计辩识法（估计辩识法（6 6）n总结上述过程：n在检测出存在不良数据后，应用完善的可疑数据检测系统检测

23、出s个可疑数据n按照式（38）求出G矩阵。n按照式（39）求出vs的加权最小二乘估计n根据 ，将求出的vs的估计值 代入vsm向量的相应元素，利用公式（42）对原来状态估计进行修正。修正后的估计值为：n修正：(43)版权所有东南大学电气工程系东南大学电气工程系 估计辩识法（估计辩识法（7 7）n上述修正方法是建立在系统线性化基础上的，并假定不良数据的存在对H阵没有影响。n估计辩识法可以快速处理多个不良数据，无需重新状态估计计算，但在准确辩识可疑数据方面含有缺陷。n逐次型估计辩识法：结合残差搜索法中准确辩识可疑数据的优势和总体型估计辩识法中快速处理多个不良数据的能力。版权所有东南大学电气工程系东

24、南大学电气工程系四、其它四、其它n状态估计接口SCADARTNETDPFLOWCASE管理AGCOPFRTCABUSLFSSLFnSCADA：状态估计从SCADA取实时遥信遥测数据，并把计算结果和遥测质量信息反送SCADAnAGC：状态估计从AGC取厂站发电计划和机组计划，当量测不足时，可用计划数据作为伪量测版权所有东南大学电气工程系东南大学电气工程系nSSLF：状态估计向超短期负荷预计送系统发电量和系统总负荷；nOPF：状态估计向最优潮流送实时方式；nRTCA：状态估计向实时预想事故分析送实时方式；nCASE管理：状态估计向CASE管理送整点方式。状态估计实时运行时，每到整点保存实时方式，但覆盖昨天此整点保存的方式。在调度员潮流CASE管理中可见到由状态估计保存的最近24个整点方式及保存时间。nDPFLOW：状态估计向调度员潮流送实时方式；nBUSLF：状态估计向母线负荷预测送负荷模型系数，量测不足时取母线负荷预测值；版权所有东南大学电气工程系东南大学电气工程系n非二次准则方法nRobust estimation版权所有东南大学电气工程系东南大学电气工程系结 束n版权所有