复数域加权最小二乘法在基于相量量测的状态估计中的应用Application of complex field weighted least square algorithm in state estimation based on phasor measurement
ABSTRACT：State estimation based on phasor measurement can improve the accuracy of state. The widely placed phase measurement uints(PMU)in high-voltage network provide a substantial foundation for state estimation only using PMU.Studying state estimation only using PMU is great significance. A measurement equation in complex domain is established, and the complex field weighted least square algorithm(CWLS), which needs no iteration to achieve estimation results, is introduced to solve the model. Research of phasor error characteristics in complex domain and the robustness of CWLS algorithm are conducted,so as an effective method of choosing weight in CWLS as well as bad data identification. The simulation results show the correctness of the model and the CWLS algorithm proposed is effective and strongly adaptive.
KEY WORDS：phasor measurement unit；phasor error characteristics; linear state estimation；weight；improved normalized residual
摘要：基于相量量测的状态估计能提高状态估计的精度，而相量测量单元在高电压等级网络中的广泛配置为电力系统全相量量测状态估计提供了坚实的硬件基础，研究全相量量测的状态估计意义重大。本文建立了基于全相量量测的状态估计量测方程，在此基础上建立了复数加权最小二乘法(Complex Field Weighted Least Square，CWLS)状态估计模型。本文还对复数域下量测量的误差特性和CWLS算法的抗差性进行了研究，同时详细分析了CWLS算法中权重的选取和不良数据的辨识。仿真结果表明了所建模型的合理性以及CWLS 算法的有效性和适应性。
关键词：相量测量单元；相量误差特性；线性状态估计；权重；改进标准化残差
0引言
自从上世纪70年代 F.C.Schweppe等人提出电力系统状态估计以来[1-3]，经过四十多年的发展和应用，状态估计在可观测性分析、估计算法、不良数据及拓扑错误的检测和辨识等方面取得了很大的进展[4-8]，状态估计的实用化进程也不断得以推进。
传统状态估计基于数据采集与监控（Supervisory Control and Data Acquisition，SCADA）系统提供的数据进行估计。近年来，PMU及建立在PMU基础上的广域测量系统（Wide Area Measurement System，WAMS）不断发展，从而为状态估计提供了新的数据源。PMU有较高的采样速率和精度，能够直接测量电压和电流的相角，且采样数据都打上了高精度的时标，因此数据的精度和速度相比SCADA系统有着巨大的提升。
自PMU装备电力系统以来，如何将PMU 量测数据应用到状态估计中是一个比较热点的研究方向，国内外一大批科研工作者对此进行了大量的研究，取得了丰富的成果[9-13]。这些研究大多是将PMU数据和SCADA数据进行混合，从而建立基于混合量测的估计模型，但PMU和SCADA量测的采样周期不同，且SCADA采样没有时标，难以保证两个系统采样数据的同步性，从而大大影响了估计结果的精度。
目前基于全PMU数据进行状态估计的研究还不太多。文献[14]首次建立了全部节点电压相量和部分支路电流相量可量测情况下的状态估计模型，并利用加权最小二乘法（Weighted Least Square，WLS）求解该模型，但没有讨论WLS 法应用于复数域的特殊性，也没有考虑权重的选取问题。文献[15]提出只要网络中PMU的配置满足可观性，即可建立线性状态估计模型。文献[16]将电压相量和电流相量分别用实部和虚部来表示，从而建立了基于直角坐标系的实数形式的线性状态估计模型，但这使得模型中系数矩阵的维数相比相量模型扩大了一倍，增加了计算规模，同时对实部和虚部间不相关的假设也不成立。
目前，PMU在国内电网中的配置已经超过1100台，500kV及以上电压等级已基本完成配置，正在向220kV电压等级推进，海南电网已率先实现了220kV厂站的全PMU覆盖。PMU 的大量配置为进行全PMU状态估计奠定了的基
础。在这种情况下，研究基于全PMU 数据的状态估计建模方法以及相应的不良数据检测和辨识方法，具有十分重要的意义。
本文建立了PMU 配置满足可观性条件下的电力系统状态估计模型，包括直角坐标形式、极坐标形式和相量形式。从数值稳定性和效率方面考虑，推荐使用相量形式的量测方程，并建立复数域加权最小二乘法（Complex Field Weighted Least Square ，CWLS ）对该模型进行求解，得到了状态量的无需迭代的、线性的估计公式。文章还对PMU 相量的误差特性、CWLS 法中权重的选取以及CWLS 法的抗差性进行了研究，并用改进标准化残差法（Improved Normlized Residual ，INR ）对不良数据进行辨识。仿真验证了本文方法的正确性和有效性。
1 基于全PMU 数据的线性量测方程
对于一个由n 个节点、b 条支路构成的电力网络，当PMU 量测使网络可观测时，即可建立全PMU 估计模型，并可利用PMU 量测的冗余度，提高状态变量的估计精度。若PMU 的量测精度很高，则可直接将节点电压相量量测作为该节点的状态，从而实现在线状态测量。然而PMU 量测的误差不可避免，且可能存在不良数据，因而将PMU 电压相量量测直接作为系统的状态并不合适。充分利用PMU 节点电压以及支路电流相量量测数据，建立相应的状态估计模型，并采用合适的算法求解模型，是构建全PMU 状态估计模型的基础。
设网络中共配置了k 个PMU ()k n ≤，这些PMU 可给出k 个节点电压相量量测以及l 条支路的电流相量量测(2)l b ≤。
PMU 的量测量分别为电压幅值、相角以及电流幅值、相角。将电压幅值U 和电压相角θ组合构成电压复相量j Ue θ，将电流幅值I 和电流相角δ组合构成电流复相量j Ie δ。则可得到如下形式的量测方程：
U m T mea
m T I ??
????==+????????????
εU U Z I I ε (1)
式中，m U 、m I 是节点电压相量量测和支路
电流相量量测构成的矢量，其中
12
12k T
j j j m k U e U e U e θθθ??=??U ，
1212l
T
j j j m l I e I e I e δδδ??=??I ；T U 、T I 是电
压相量真值和电流相量真值构成的矢量；
U ε、I ε是相应的误差矢量，上述矢量中的
元素均为复数。
将支路电流相量T I 用节点电压相量U 来表
示，其中T cT ??
=????
U U U ，cT U 表示未配置PMU 节点的电压相量真值。考虑T I 的一个分量pq I ，对应于从节点p 流向节点q 的电流相量。对于支路pq ，不管其是否含有变压器、移相器，均可用如图1所示的π型等值电路表示。不妨假设其包含变压器和移相器，变压器的变比为1:K ，移相器的移相角度为θ，线路的电阻、电抗和充电电纳分别为r 、x 、b ，由基本电路定律可得：
y p
y q
I pq
I qp
y pq
图1 支路的π型等值电路
()()00pq p pq p pq q qp qp p q qp q I y y V y V I y V y y V ???
????
=+-?
?
?=-++?
(2)
其中，02p jb
y =
，()1pq j y r jx Ke θ
-=+，022q jb
y K =
，()()
*1qp j y r jx Ke θ-=+。 由(2)式可知，可以通过支路导纳矩阵将支
路电流相量T I 用节点电压相量U 表出。设支路导纳矩阵为b Y ，则b Y 的维数为l n ?，且是一个
稀疏矩阵，其中的非零元素定义如下：
()()()0,,-b i ij b ij Y i i y y Y i j y ?=+?
?=??
当ij I 有电流量测 (3)
因此有T =b I Y U
(4)
将等式(4)代入到等式(1)得量测方程：
12U mea I ????=+=+????
????
b b εI
0Z U A U εY Y ε (5) 式中，I 表示k k ?阶单位矩阵，0表示
()k n k ?-阶零矩阵，[]12=b b b Y Y Y ，1b Y 表示l k
?阶矩阵，2b Y 表示()l n k ?-阶矩阵。
式(5)建立了量测量mea Z 和状态量U 间的关系，其中mea Z 和U 均是复矢量。在网络结构和量测配置不变的情况下，系数矩阵A 为常稀疏复矩阵，维数为()l k n +?。因此mea Z 和U 间是线性关系，(5)即为基于PMU 数据的量测方程。
(5)的直角坐标形式和极坐标形式分别如式(6)和(7)所示：
,,11
mea r r i r r mea i i
r i i -????????
=+????????????????=+Z ΑA U εZ A A U εB U ε (6)
222
m
mea
a δδ????
??????????==
+????????????????
??=+IU Ιm
m a m m m a m θU θI 0U 0I U U Z εI U I C U H ΗH H ε (7) 其中,mea r Z 、,mea i Z 分别表示mea Z 的实部和虚部，r A 、i A 分别表示A 的实部和虚部，r U 、i U 分别表示U 的实部和虚部，r ε、i ε分别表示ε的实部和虚部。m
m
U 、a m
U 表示m U 的幅值和相角，
m m I 、a m I 表示m I 的幅值和相角，m U 、a U 表示U
的幅值和相角，IU H 、ΙθΗ、δU H 、δθH 表示电
流幅值、相角对电压幅值、相角的雅克比矩阵，
是非常数矩阵。
2 复相量的误差特性
2.1 误差定义
一般来说PMU 量测数据的幅值和相角是相
互独立的，且一般假设误差均值为0。若PMU 电压幅值量测的标准差为u v ，相角量测的标准差为u Vθ，电流幅值量测的标准差为u I ，相角量测的标准差为u Iθ。则误差复相量如图2所示
U k U kT
θkT
ε2 ε1
图2 电压量测复相量和真值复相量间的关系
对量测方程(5)中mea Z 的任一电压分量k U ?
有
()()
21kT k j j k k kT kT Uk
U U e U e U θεθεε?
+?
==+=+ (8)
其中，kT U 表示电压k U ?
幅值的真值，kT θ表示电压k U ?
相角的真值。化简(8)得：
()221e e e kT
j j j Uk kT
kT U
U θ
εεεε=+- (9)
(9)即是量测方程中任一电压量测量k U ?
的量测误差。在复数领域，一般用距离来表征两个相量间的差别，对图2所示的量测量k U ?
和真值
kT U ?
，两者间的距离为：
2122121
(,)()e 2()(1cos )j k kT Uk kT kT kT
kT d U U U U U U εεεεεε
??
==+-=+-+ (10)
因此可以用式(10)来用来表征复相量间的绝对误差。
2.2 正常量测情况下的期望和方差
假设幅值和相角服从相互独立的零均值高斯分布，则对某个量测绝对误差Uk ε，可求其期望和方差。
()()22
2
1e e e e e 1kT kT j j j Uk kT kT j j kT E E U U U E θεεθεεε??=+-????=-??
(11)
令()222e (cos sin )j E E j a jb εεε=+=+，则上式可简化为()e (1)kT j Uk kT E U a jb θε=-+，从而可以得到Uk ε的期望和方差：
()22(1)Uk kT E U a b ε=-+ (12)
()()()2
2
2
2
222
2222()2(1)(1)(1)
Uk Uk Uk kT v kT v kT D E E U a u U a b u U a b εεε??=-??
=-+--+=+-- (13)
由(12)和(13)可知，量测误差的期望和相角
的误差有关，与幅值误差无关；而方差同幅值和相角的量测误差均有关。由于a 和b 的值无法用解析方法求得，在实际应用中可用统计分析方法来估计()Uk E ε和()Uk D ε的值。
假设PMU 电压相角量测的标准差u Vθ分别为0.02o 、0.1o 、0.5o 、1o 、3o ，对每种标准差u Vθ分别对2ε进行10000000次蒙特卡洛采样，对得到的2ε值分别计算()Uk E ε和()Uk D ε，然后求其平均值，得到的结果如表1所示。
PMU 相角量测的误差标准差一般可以达到0.02o ，最差的情况下误差也不会超过1o 。从表1中对()Uk E ε和()Uk D ε进行的仿真统计分析可以得出在PMU 量测的误差范围内，误差随机变量的平均值约为0，方差基本为222V V kT
u u U θ+。因此在实际应用中，可以认为()0Uk E ε≈，
()222
Uk V V kT D u u U θε≈+。
表1 ()Uk E ε和()Uk D ε的数值计算值
测量误
差 ()Uk E ε ()Uk D ε
2
V u θ
kT U
22+V kT u U
0.02o 1.4021*10-7
1.2188*10-7 1.2185*10-7 0.1o 1.7326*10-6 3.0458*10-6 3.0462*10-6 0.5o 3.8675*10-5 7.6128*10-5 7.6154*10-5 1o 1.5230*10-4 3.0444*10-4 3.0462*10-4 3o
1.3703*10-3
2.7385*10-3
2.7416*10-3
对电流相量量测误差可做类似的分析，可得
()0Ik E ε≈，()222
Ik I I kT D u u I θε≈+。
2.3 有不良数据时的特性
研究表明，PMU 量测数据中，幅值的精度
一般较高，而相角的精度相对差一点。在准稳态情况下，电压幅值的波动远低于1%，一般不超过0.4%，而相角的波动相对较大，出现错误数据的概率也相对大些。不妨设某次采样的数据中幅值是准确的，而此时的相角数据是错误的，错误幅度为1o ，即21o ε=。此时可得绝对误差Uk ε为：
(
)
1e 1 1.75%o
j Uk kT kT U U ε=-=? (14)
可见即使只有1o 的相角误差，就使得绝对误差高达幅值的1.75 %，大大超过幅值的正常量测误差，因此绝对误差对相角误差的灵敏度非常高。
3 复数域加权最小二乘法
WLS 是电力系统状态估计的实用方法之一，但WLS 没有抗差性，在实际应用时，在WLS 估计之后需要加上一个不良数据辨识环节。 WLS 是应用在量测方程中的各量为实数的基础上的。而建立的方程(5)是基于向量的复数线性方程；基于直角坐标的方程(6)虽是实数线性方程，但其关系矩阵B 的维数为
()222l k n +?，比矩阵
A 大一倍，由于计算中涉及关系矩阵的求逆，将增加计算量，另外将PMU 量测转换为实部和虚部时，还将引入转换误差，
同时实部和虚部之间是相关的，不再相互独立；基于极坐标的方程(7)是实数非线性方程，需要进行迭代计算，存在收敛性问题。因此本文选择方程(5)作为全PMU 状态估计的量测方程，并将WLS 引入复数域建立复数域加权最小二乘法CWLS 求解之。
当量测方程中的各量是复数时，由于量测误差矢量ε是一个复矢量，基于极大似然估计原则，CWLS 的目标函数是使得ε的模的加权平方和最小，即：
()[][]
*
*
min T
T
T
mea mea J ===--U εW εεW ε
Z A U W Z A U (15)
(15)式中W 是权重矩阵，为()l k +阶实对角方阵。推导后可以得到状态量U 的估计结果如下：
*
*
1
T
T est mea
mea
-??=????=U A W A A W Z D Z (16)
模型(15)是一个二次规划问题，可保证得到全局最优解，而(16)式不需要迭代，通过一次运算就可以得到状态量的估计值。因此不存在迭代计算中常见的收敛性和快速性等问题。实际上，在网络结构和量测配置不变的情况下，式中的D 为常量，因此状态量的估计结果可以通过直接计算得到。实际计算中只需先形成支路导纳矩阵b Y ，然后形成矩阵A ，之后得到权重矩阵W ，即可完成估计计算。
由于CWLS 法是以绝对误差的加权平方和最小为目标的，因此CWLS 法对相角错误数据有较大的抗差性，仿真算例的结果也证明了这一点。
4 权重矩阵W 的选取
权重矩阵一般是对角矩阵，其对角线元素同量测误差的方差有关，在用CWLS 进行估计的过程中，可以通过2种方法获取权重矩阵W ：定权重法和自适应调整权重法。
定权重一般是通过对量测误差的先验知识以及一定的误差分布假设（通常假设为高斯分布），给相应的量测一个固定的权重，一般取量测误差方差的倒数作为相应量测的权重。根据第
2节中关于误差特性的介绍，可以得到定权重情况下，权重矩阵的各元素如下：
()()()()122211
222,,V v kT
I I kT u u U W i i R i i u u I θθ---?+?==??????+?
(17) 实际应用时，可以用测量值代替式(17)中的
真值。
自适应调整权重主要是依据量测残差和量测误差间的关系来完成的。利用CWLS 完成估计计算后，可得量测残差r 为：
**
1mea est
T T -=-????=-=????????
r Z Z I A A W A A W εS ε
(18) 其中S 为残差灵敏度矩阵，是实数矩阵。在网络结构和量测配置不变的情况下，S 矩阵为定常矩阵。
残差r 的协方差为：
*
*
*
()()()T
T
T T
T
T
E E E =====r R rr S εεS S εεS SRS SR
(19)
R 为误差协方差矩阵，
是实对角矩阵。式(19)通过灵敏度矩阵S 建立起了残差协方差r R 和误差协方差R 之间的关系，因此通过r R 和S ，经式(20)可以求得R 中的元素：
(,)(,)(,)
r R i i R i i S i i =
(20)
对R 中的元素求倒数，即可得到相应的权重，从而可以对权重矩阵中的元素进行自适应的调整，使权重同实际的量测误差相匹配。
5 不良数据的检测和辨识
最大标准残差（Largest Normlized Residual ，
LNR ）法是应用WLS 方法进行状态估计后，一种实用的不良数据检测、辨识方法。针对CWLS 方法的特性，对LNR 法进行改进，可以将其用于CWLS 法估计后的不良数据检测和辨识中。
据式(18)，可以将复误差矢量ε写作：
()()=++-=-εr I S εS εI S ε
(21)
两边取模值，可得绝对误差
()+=-=+εS εI S εr t
(22)
矢量ε所在的向量空间为l k R +，定义l k R +上的内积为,T <>=αβαW β。据此定义，CWLS 的目标函数即为矢量ε在l k R +上的内积。由式(22)可知绝对误差矢量ε可以分解为两个分量，即残差的模矢量=r S ε和矢量()=-t I S ε，可以证明这两个矢量的内积为零，这两个矢量在l k R +空间上是相互正交的。因此残差的模矢量r 只是绝对误差矢量ε的一部分。
假设第i 个量测数据是不良数据，即i c =ε，其余分量为零。此时有i c =r S ，()i c =-t I S ，其中i S 、()i -I S 表示矩阵的第i 列构成的列矢量，对第i 个量测，定义指标：
()(),,i i i i i c c K c c <-->
=
<>
I S I S S S (23)
定义改进残差（Improved Residual ，IR ）ir 来更好的表征误差矢量ε：
21i i i K =+ir r (24)
其中i r 表示残差模矢量的第i 个分量，i
ir 表示改进残差的第i 个分量，并用ir 进行不良数据的检测和辨识。
6 仿真研究
这一节将在IEEE 节点系统上，对前面提出的全PMU 状态估计的CWLS 方法进行仿真研究和验证。
6.1 估计结果评价
在电力系统状态估计的仿真研究中，可用状
态估计值与真值之差的绝对值均值来表征状态估计结果的优劣，其表征式如下所示：
,11
11M
N
average
ij estimate j i j x x M
N
σ==??
=- ???
∑∑
(25)
其中，σaverage 表示M 次蒙特卡洛分析中，估计误差绝对值的平均值，N 表示系统状态变量的个数，x ij,estimate 表示第i 次蒙特卡洛分析中第j 个状态量的估计值，x j 表示第j 个状态量的真值。对每一次蒙特卡洛分析，量测值是通过在真值基础上加随机白噪声得到的，白噪声的分布由量测量的误差分布来确定。
6.2 测量误差
PMU 测量的误差来源主要有：PT 、CT 测量误差、光缆通道传输误差、A/D 转换误差和算法误差。PT 、CT 测量误差在工频下一般是固定的，通道传输误差与通道的长度有关，这两种误差一般可以通过外部校准设备进行补偿。本文主要考虑由A/D 转换误差和算法误差造成的量测不确定度。本文的量测误差是根据国标对PMU 量测的最低标准来设定的，如表2所示。
表2 最大量测误差
PMU 量测
电压幅值 电压相角 电流幅值 电流相角
0.4% 0.5o 0.4% 1o
6.3 无不良数据时估计结果
对IEEE300节点系统，选节点257为相角参考节点，在全部300个节点上均装设PMU ，可以测量300个电压相量和822个电流相量。对该标准系统，进行1000次蒙特卡洛试验，对每一次蒙特卡洛试验得到的量测数据，分别进行估计计算，得到状态量的估计值。其中电压幅值单位是标幺值，相角单位是度。表3列出了IEEE300节点系统在给定误差下的σaverage 及完成状态估计的计算时间，仿真的硬件条件是Intel Core i5 2.27GHz CPU ，内存为2GB 。表4列出了用该方法对IEEE14节点系统、IEEE118节点系统、IEEE300节点系统和波兰3012节点系统进行1000次仿真所用的平均时间。分析表3和表4的仿真结果可以得出以下几点结论：
（1）本文提出的模型和方法能快速有效的提高测量值的准确度，较好的完成状态估计计算，尤其是大幅度的提高了相角的精度。
（2）完成状态估计所需时间同系统规模大致为线性关系；这是因为全PMU 状态估计模型为线性模型，估计算法不需要经过迭代即可完成。
表3 给定误差下估计的σaverage 和时间
σaverage 幅值误差/10-3 角度误差/度 时间/ms
量测 1.77 0.30 0.18
估计
1.58
0.06
表4 不同规模系统估计所需时间
系统 时间/ms IEEE14节点 0.22 IEEE118节点 1.09 IEEE300节点 2.38 波兰3012节点
29.10
表5 假设及实际的误差水平
PMU 量测
电压幅值 电压相角 电流幅值 电流相角 假设误差 0.2% 0.4o 0.3% 0.6o 低于假设误差 0.2% 0.2o 0.3% 0.2o 高于假设误差
0.4%
0.5o
0.4%
1o 为了比较当测量系统实际的误差水平同假设水平不一致时，本文提出的状态估计方法的适应性。以上述300节点系统为对象，以表5中列出的假设及实际误差水平为条件，进行仿真计算，得到表6所示的结果。
从表5和6中可以看出，当量测量的实际误差水平同假设的误差水平不同时，不管是实际误差水平高于或者低于假设，本文提出的方法都能够有效的提高测量值的准确度，表现出较强的适应性。
表6 算法的适应性结果
状态估计结果
误差水平 低于假设误差时 高于假设误差时
幅值/10-3 角度/度
幅值/10-3 角度/度 估计前 0.88 0.12 1.76 0.30 估计后
0.53 0.02
1.47 0.05
6.4 有不良数据时估计结果
对IEEE9节点系统，选节点1为相角参考
节点，在全部节点上均装设PMU ，可以测量9个电压相量和18个电流相量。在电压量测2U 、
5U 、6U 和8U 处分别设置1%幅值坏数据、1o 相
角坏数据、5o 相角坏数据，15o 相角坏数据，利用CWLS 算法进行估计计算后，可以得到表7所示的结果，其中坏数据及其估计结果在图中用粗体显示。
表7 多个坏数据时状态估计结果
真值
量测 估计 序号 幅值
角度 幅值 角度 幅值 角度 1 1.000
1.001
0 1.022
2 1.000 9.669 2.012 9.65
3 1.038 9.625 3 1.000 4.771 0.999
4.987
1.027 4.819
4 0.987 -2.407 0.988 -2.589 1.007 -2.303
5 0.975 -4.017 0.975 -6.967 0.994 -3.854
6 1.003 1.926 1.005 8.942 1.029 2.059
7 0.986 0.622 0.984
0.532
1.010 0.856
8 0.996 3.799 0.996 22.586 1.024 3.970 9
0.958 -4.350 0.958 -4.593 0.976 -4.115 从表7中可以看出，即使不经过任何坏数据
的处理，CWLS 对相角的估计结果仍然是可接受的，8U 处15o 的相角坏数据都能够被很好的抑制，体现了CWLS 法对相角误差极强的抗差性。 各电压相量量测的改进标准化残差如表8所示。
表8 各相量量测的改进标准化残差
序号
1
2
3
4
5
6
7
8
9
INR 6.9 166 10.6 7.5 20 41 8.2 109 6.3
从中可以看出，含有坏数据的量测的改进标准残差INR 远大于正常量测的INR ，从而可以有效的对坏数据进行辨识。
7 结论
本文建立了基于复数形式的全相量量测状态估计量测方程，在此基础上建立了复数域加权最小二乘法对该模型进行求解，通过直接计算得到状态估计值，而无需进行迭代。
本文将估计误差绝对值的平均值以及估计计算的时间作为评价标准，通过蒙特卡洛仿真方法，对IEEE300节点系统及其他节点系统进行了仿真分析，结果表明了建立的全相量量测状态估计模型及算法的正确性和有效性。该方法可以有效提高量测量的准确性，尤其是角度的准确性，同时极大的缩短了估计运算的时间，对300节点的系统只需2.4ms 即可完成估计计算。同时对实际量测误差水平同假设不一致的情况进行了仿真，结果表明了该方法在这种情况下依然有效，体现了较强的适应性和抗差性。下一步将重
点研究CWLS算法在网络拓扑和参数存在错误情况下的特性。
参考文献
[1]. Schweppe F.C.，D.B. Rom，Power System Static-State
Estimation, Part II: Approximate Model. Power Apparatus and Systems, IEEE Transactions on, 1970. PAS-89(1): p. 125-130.
[2]. Schweppe F.C. and J. Wildes, Power System Static-State
Estimation, Part I: Exact Model. Power Apparatus and Systems, IEEE Transactions on, 1970. PAS-89(1): p. 120-125.
[3]. Schweppe F.C., Power System Static-State Estimation, Part
III:Implementation. Power Apparatus and Systems, IEEE Transactions on, 1970. PAS-89(1): p. 130-135.
[4]. 刘莉,翟登辉与姜新丽,电力系统不良数据检测与辨识方法的
现状与发展.电力系统保护与控制,2010(5):第143-147+152页.
LIU Li, ZHAI Deng-hui, JIANG Xin-li. Current situation and
development of the methods on bad-data detection and
identification of power system[J]. Power System Protection
and Control, 2010, 38 (5):143-147+152.
[5]. 亓俊健等,电力系统抗差状态估计研究综述.电工电能新技术,
2011(3): 第59-64页.
QI Jun-jian,et al.A Review of Power System Robust State Estimation[J].Advanced Technology of Electric Engineering and Energy,2011(3):59-64.
[6]. 高奥,兰华与李晋,电力系统状态估计算法综述与研究趋势.
东北电力大学学报, 2011(1): 第52-56页.
GAO Ao,LAN Hua,LI Jin.Review and Development on State Estimation of Electrical Power Systems[J].Journal of Northeast Dianli University,2011(1):52-56.
[9]. 丁军策,蔡泽祥与王克英,基于广域测量系统的状态估计研究
综述.电力系统自动化, 2006(7): 第98-103页.
DING Jun-ce , CAI Ze-xiang , WANGKe-ying.An overview
of state estimation based on wide area measurement
system.Automation of Electric Power Systems,2006,30
(7):98-103.
[10]. 秦晓辉,毕天姝与杨奇逊,计及PMU的混合非线性状态估计
新方法. 电力系统自动化, 2007(4): 第28-32页.
Qin Xiao-hui,Bi Tian-shu,Yang Qi-xun.A New Method for Hybrid Non-linear State Estimation with PMU[J]. Automation of Electric Power Systems,2007(4):28-32
[11].李卫国等, 基于广域测量系统的电力系统状态估计. 电网技
术, 2008(S2): 第124-127页.
Li Wei-guo,et al.Power System State Estimation Based on Wide-Area Measurement System.Power SystemTechnology,2008(S2):124-127.
[12]. P hadke, A.G., et al., Recent developments in state estimation
with phasor measurements, in 2009 IEEE/PES Power Systems Conference and Exposition, PSCE 2009. 2009: Seattle, WA, United states.
[13]. 李大路等, 利用PMU数据提高电力系统状态估计精度的方
法. 电网技术, 2009(3): 第74-78页.
Li Da-lu,et al.An Algorithm Improving Precision of Power System State Estimation by PMU Data[J]. Power System Technology,2009(3):74-78.
[14]. P hadke, A.G., J.S. Thorp and K.J. Karimi, State Estimlatjon
with Phasor Measurements. Power Systems, IEEE Transactions on, 1986. 1(1): p. 233-238.
[15]. H ongga, Z. A New State Estimation Model of Utilizing PMU
Measurements. in Power System Technology, 2006. PowerCon 2006. International Conference on. 2006.
[16]. 李强等, 基于相量量测的电力系统线性状态估计. 电力系统
自动化, 2005(18): 第24-28页.
Li Qiang,et al.Power System Linear State Estimation Based on Phasor Measurement[J]. Automation of Electric Power Systems,2005(18):24-28.