混合量测电力系统多区域分布式状态估计

混合量测电力系统多区域分布式状态估计

0分布式状态估计的基本原理

中国各地区电网的互联将带来新的全球形势评估和计算带来新的挑战。同时，国内电网调

度采用分层控制、分布处理的模式,导致各区域电网数据无法共享，因而传统的整体式状态

估计难以符合现代电网的发展需求

为此，文献［2］提出了一种双层次分布式方法，在第1层次各子区域独立进行状态估计，将边

界信息传递至中央协调侧后，第2层算法协调各区域参考母线的相角以及边界母线的状态

量。文献［3］以辅助问题原理将全网状态估计分解为独立的子问题。文献［4］提出了根据注

入功率对状态量的灵敏度矩阵修正子区域状态量的分解协调算法。文献［5］利用搭接

式分区、外网浓缩等值法，提出了基于黑板模型的分解协调计算机制。文献［6］以节点撕裂

法切分互联电网，根据KKT(Karush-KuhnTucker)条件方程分解协调全网状态估计。文

献［7］以支路切割法将网络分块后分解协调计算。文献［8］将异步迭代

此外,相量测量单元(PMU)的配置促进了分布式状态估计的发展。文献［1475］在每个子区

域的参考母线配置PMU,从而实现各区域参考母线到全回参考母线的转化。文献［16-17］的

分布式方法计及了数据采集与监控(SCADA)系统、PMU泡合量测。文献［18］利用PMU量测将

状态估计转化为线性最小二乘问题，以交替方向乘子法

鉴于当前(或未来)电网的量测配置,本文的分布式状态估计恻重于SCADA系统、PMU混合量

测。同时,考虑到非线性SCADA系统量测影响了各子区域内部状态估计以及相邻区域间协

调计算的收敛性,本文研究了SCADA系统、PMU混合量测下的双线性

1基于admm的线性分布估计

1.1边界植物增加

扩展子区域法分区原理如羽1所示，各扩展的子区域除原有的区域外，还包含r相邻了区域

的边界母线。

由图1可以看此相邻子区域k与子区域1边界部分为母线a,b,c,d,则扩展子区域k包含

了子区域1中的边界母线b和d,同样扩展子区域1包含了子区域k中的边界母线a和Co

边界母线的注入功率量测一般难以处理

定义x为全网的状态变量,x

式中:N

1.2加权最小二乘估计

假设各子区域量测量与状态量呈如下线性关系：

式中:Z

基于广泛应用于工程实践的加权最小二乘（WLS）估计，计及边界状态约束（式（1））的子

区域k的状态估计目标函数为：

式中：f（x

由于子区域k的状态估计目标函数（式（3））的等式约束包含了子区域1的边界状态量x

引入中间辅助变量x

1.3admm求解

作为一种分布式凸优化算法,ADMM在电力系统分布式计算中已有应用

式中:w为状态变量;h（・）和g（•）为凸优化函数;A和B为恒定矩阵;c为恒定向量。

式⑸的伸缩增广拉格朗口函数为：

式中：P为惩罚因子;u为伸缩拉格朗日乘子。

式⑹中含状态变量x和，及伸缩拉格朗日乘子u,采用ADMM进行求解,x,w,u可通过解耦

分步计算出:

式中：t为迭代次数。

同理,基于ADMM求解式（41,扩展子区域k的状态量｛x

式中:u

在式（10）至式（12）中,ADW/I以类似高斯一赛德尔的方法，将第t次迭代值用于求解第t+1次

状态最，其中式（10）求解扩展子区域k的状态变量x

需要说明的是，引入中间辅助变量x

式中:u

由附录A式（A5）可知，在ADMM迭代过程中，式（13）恒成立。因而通过引入中间辅助变量

推广到全网各子区域,基于上述方法求解，各相邻子区域间在迭代过程中仅需交换边界母线

的状态量信息，无需中央协调侧,最大限度地保留了各子区域的独立性,简化了通信机制。

然而，需要注意的是,ADYM必须在严格凸优化的条件才能可靠收敛

2scada系统与pmu的混合体积模型

2.1基于scada系统测量的二元模型

基于SCADA系统量测,通过引入中间变量，可将传统的非线性WLS估计转化为三阶段问题

式中:z为量测量;e

2.1.1油路两端投资

1）对于网络中的每条支路:定义如下变量：

式中：i和j为支路两端母线;V

2）对于网络中的每条母线,定义如下变量：

SCADA系统量测量z包括支路功率量测P

式中

2.1.2第2阶段的非线性变换

第2阶段非线性变换为等维数变换：

定义变量v=[a

2.1.3第3阶段线性测量模型

定义状态变量x=[”，。],其中d和o分别为由&

2.2参考汽车到全网参考母线

相比于传统的SCADA系统量测,PMU可直接测量母线电压相角。在分布式状态估计中，可在

每个区域的参考母线配置PUU量测，实现各区域参考母线到全网参考母线的转化

式中：（）

假设全网参考母线在区域1,则区域1第3阶段v与x的线性关系仍由式（24）描述;对于其

余子区域,v与x的线性关系为：

此时x包含参考母线的相角。

1）估计精度:分布式状态估计结果需与整体式方法相同，或者在允许的工程误差内，估计精

度尽可能接近整体式方法，

2）计算效率:相比于整体式方法，分布式方法的一大优点在于提高了状态估计的计算效率。

3）通信复杂度：区域间交换的信息量应尽可能少，以减少通信延时，降低通信复杂度。

4）收敛的鲁棒性:分布式方法需在不同的分区、拓扑、运行状态、量测配置下可靠收敛。

5）不良数据分析：分布式方法应保留整体式方法的不良数据辨识能力，易于计算正则化残差,

特别是对于边界区域。

下文将从上述5个方面,将本文所提方法与文献［17］的现有方法进行比较分析工

4.1节点电压幅值与管路相角估计误差

在工程误差允许的范围内，分布式方法的估计精度应尽可能接近整体式方法。表1和表2

为相比于整体式方法（即MS估计），分布式方法（本文方法以及文献［17］的现有方法）的节

点电压幅值与支路相角的估计误差。

由表1和表2可知,本文方法的估计精度要高于现有方法,原因在于本文方法无需中央协调

侧,从数学意义上讲,当式（1）成立时,分布式方法的估计结果与整体式方法是等价的;而现

有分布式方法各子区域、协调侧分开求解,其估计结果一般为次优解,难以与整体式方法保

持一致

4.2各节点侧求解时间

本文方法的计算时间等价于各子区域计算时间的最大值，文献［17］的现有方法计算时间除

分区侧计算时间外，还包括通信时间、协调侧求解时间。IEEE118节点系统、实际电网

的整体式与分。式力法的计算时间如表3所示。

由表3可知,相比于整体式方法,分布式方法降低了系统的求解规模，因而两种分布式方法

的计算效率均有所提高。比外，由于本文引入了双线性理论，保持了雅可比矩阵恒定，并且

以ADMM求解,简化了通信机制,无需中央协调侧，因而本文方法计算效率也高于现有方法。

4.3相邻区域边界母件

两种分布式方法的通信复杂度比较如表4所示。本文方法基于扩展子区域法分区，扩展子

区域除包含原来的子区域外,还包含了相邻区域的边界母线，易于处理边界母线的注入量测,

以ADMM迭代求解,相邻区域仅需交换边界母线的状态量信息，无需中央协调侧。而现有方

法首先由各子区域内部执行状态估计，然后将边界母线的估计值与协方差（需要由复杂的矩

阵求逆获得）一起传递至中央协调侧，由中央协调侧计算边界状态量并协调各子区域的参考

相角，因而现有方法通信复杂度更高。

4.4不良数据处理能力

文献［17］中现有的分布式方法包含分区侧、协调侧两步状态估计，收敛性能与整体式方法

相当;且易于以最大正则化残差法检测辨识不良数据，保留了整体式方法的不良数据处理能

力。

本文方法的收敛性由ADMM保证,在严格的凸优化条件下,ADMM能够可靠地收敛至全局最优

解,在本文的测试算例中，第1、第3阶段的分布式状态估计，IEEE118节点系统需迭代4%

次，实际电网需迭代7〜8次;对于不良数据的检测辨识,传统的最大正则化残差法也同样适

用于本文方法。此外，由于抗差估计能够自动抑制不良数据对估计精度的影响，避免了繁琐

的不良数据检测辨识程序(残差灵敏度矩阵计算复杂度较高，难以适应大规模电力系统的发

展需求

5计算方法的适用性

本文研究了SCADA系统和PMU混合量测下的电力系统分布式状态估计,具有如下特点：1双

线性理论与ADMM的结合，保证了分布式状态估计能够快速可靠的收敛，同时也降低了计算

复杂度与通信复杂度;2各子区域的分布式计算中仅需交换边界母线的状态量信息,无需中

央协调侧，最大限度地保留了各子区域的独立性;3相比于整体式力法，在工程误差允许的范

围内，当边界母线状态量相等的