CN104036435A - 一种电力系统状态估计的量测量权重设置方法 - Google Patents

Abstract

本发明公开了一种电力系统状态估计的量测量权重设置方法，包括以下步骤：根据支路数据统计各节点的出线数α_j；判断当前量测量是否为电压幅值量测量；判断当前量测量是否为支路功率量测量；判断当前量测量是否为节点注入功率量测量；设置当前量测量的类型因子为k_wi＝c/α_j；设置当前量测量权重为由于本发明设置权重时同时考虑量测量误差和量测量类型，节点注入功率量测量的类型因子k_w为支路功率量测量的类型因子除以该节点的出线数。实施算例表明本发明可以明显提高电力系统状态估计的精度。对IEEE30系统，各量测量均存在1％误差时，采用本发明计算，状态估计的精度比传统方法提高7.4％。

Description

一种电力系统状态估计的量测量权重设置方法

技术领域

本发明涉及一种电力系统的状态估计方法，特别是一种加权最小二乘法状态估计的量测量权重设置方法。

背景技术

状态估计是电力系统调度部门所使用的能量管理系统(EMS)的重要组成部分。随着电网规模的逐年加大，电力系统的自动化水平也越来越高。为了保证电力系统安全经济地运行，要求现代化的电力调度系统能准确、迅速、全面地掌握电力系统实际的运行状态，用来分析当前系统的运行状态，预测系统的运行趋势，并能及时地对运行中发生的各种问题，提出处理方法，制定下一步的操作方案。

为了监视电力系统的运行状态，在电力系统的各个厂站(发电厂和变电站)都装设了测量单元来获取各种电气量的实时数据，这些被称为量测数据的实时数据通过各种通信手段传送到电网调度中心。由于量测数据存在一些问题，因此无法直接依据量测数据来判断电力系统状态。量测数据有如下两个缺点：

(1)量测数据不全。量测数据并不包含电网的全部数据，仅仅采集了电网的部分数据。

(2)量测数据不准。电网调度中心得到的实时数据是通过传感器、变换器等远动装置传送进来的，数据采集、变换和通信每个环节都有误差，并且存在不同程度的干扰，因而这些数据不是完全可靠的，甚至存在很大的误差。

由于存在以上缺点，直接使用这些量测数据分析系统的状态显然是不可靠的，不能满足电网调度的要求。量测数据不能直接使用，必须加工处理，补齐不足数据，剔除误差很大的坏数据，减少误差影响，才能使用。这个加工处理的过程就是状态估计。

加权最小二乘法是状态估计最基本、最常用的算法，最小二乘法就是根据最小二乘准则使目标函数取得极小值。

加权最小二乘法状态估计的原理如下：

令电力系统量测方程组为

z＝h(x) (1)

式中，z为m维量测量向量；x为n维状态变量向量(由节点电压幅值和电压相角组成)；h(x)为m维量测量的计算值函数向量，由状态变量求得，是非线性函数。

量测量包括支路有功功率量测量、支路无功功率量测量、节点注入有功功率量测量、节点注入无功功率量测量和节点电压幅值量测量。

由于量测方程的个数m大于状态变量个数n，即m>n，方程一般是无解的，即不存在满足式(1)中所有方程的解。需要寻找一组状态变量，使方程的加权残差(量测量与量测量的计算值的差)最小。

令目标函数为

$J\left(x\right)={(z-h\left(x\right))}^T{R}^{-1}(z-h\left(x\right))=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{R}_i^{-1}{(z_i-h_i\left(x\right))}^2=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{R}_i^{-1}{r_i}^2---\left(2\right)$

式中，R^-1表示量测量权重矩阵；r_i＝z_i-h_i(x)称为第i个量测量的残差。

传统的量测量权重设置方法一般是选择量测量的方差倒数为权重值，即

$R_i^{-1}=1/{\mathrm{\σ}}_i^2---\left(3\right)$

式中，为第i个量测量的方差。

最小二乘法就是要使得目标函数J(x)最小。J(x)的偏导为零是求取J(x)极值的必要条件，即

$\frac{\mathrm{dJ}\left(x\right)}{\mathrm{dx}}=-{2H}^T\left(x\right)R^{-1}(z-h\left(x\right))=0---\left(4\right)$

式中，H(x)是量测量雅可比矩阵，由量测量的计算值函数h(x)对x的偏导数构成。

式(4)是一个非线性方程组，可以用牛顿法求解。

令f(x)＝H^T(x)R^-1(z-h(x))＝0，将x在x₀点附近进行泰勒展开，并取至线性项，得

$f\left(x_0\right)+\frac{\∂f\left(x\right)}{{\∂x}^T}{|}_{x_0}\mathrm{\Δx}=0---\left(5\right)$

式中，

$\frac{\∂f\left(x\right)}{{\∂x}^T}=\frac{\∂}{{\∂x}^T}\left[H^T\left(x\right)R^{-1}(z-h\left(x\right))\right]\≈-H^T\left(x\right)R^{-1}H\left(x\right)---\left(6\right)$

式(6)代入式(5)，得状态估计的法方程

H^TR^-1HΔx＝H^TR^-1(z-h(x₀)) (7)

式中，方程系数部分H^TR^-1H称为信息矩阵。

解法方程，得

Δx＝[H^TR^-1H]^-1H^TR^-1(z-h(x₀)) (8)

用Δx对初值x₀进行修正，得x为

x＝x₀+Δx (9)

由于对f(x)做一阶泰勒展开，当x₀十分接近x的真值时，一阶泰勒展开才能足够精确，但是实际上很难取到足够接近x真值的初值x₀，所以应该反复迭代来逐步修正x使它逼近x的真值。迭代公式为

Δx^(k)＝[H^T(x^(k))R^-1H(x^(k))]^-1H^T(x^(k))R^-1(z-h(x^(k))) (10)

x^(k+1)＝x^(k)+Δx^(k) (11)

按式(10)和式(11)进行迭代修正，直到为止(ε为收敛精度)，这时得到的就是最佳估计值，使目标函数J(x^(k))取得最小值。

如图1所示，加权最小二乘法状态估计方法，包括以下步骤：

步骤1：读取网络数据和量测数据z；

步骤2：按平启动设置电压初值，即所有节点电压幅值为1.0，电压相角为0.0；

步骤3：令迭代计数k＝1；

步骤4：计算量测量的计算值函数对x的偏导数的雅可比矩阵H(x^(k))、量测量的计算值函数h(x^(k))、信息矩阵H^TR^-1H、方程右端向量H^TR^-1(z-h(x^(k)))；

步骤5：解法方程(H^TR^-1H)Δx＝H^TR^-1(z-h(x^(k)))求Δx^(k)及max|Δx_i|；

步骤6：判断max|Δx_i|是否小于收敛精度ε，如果小于ε，结束；否则执行步骤7；

步骤7：令x^(k+1)＝x^(k)+Δx^(k)、k＝k+1，返回到步骤4进行下一次迭代。

上述计算方法中，量测量的权重设置对状态估计计算精度的影响很大，传统的量测量权重设置方法一般是选择量测量方差的倒数为权重值，实践中发现误差相同的不同类型的量测量对状态估计的影响程度有很大不同。设置权重时不考虑量测量类型对状态估计精度的影响必然影响状态估计的精度。中国专利ZL201310314260.1中披露了一种电力系统状态估计的量测量权重设置方法，状态估计的精度有了较大提高，但仍然有进一步提高的余地。

发明内容

为解决现有技术存在的上述问题，本发明要设计一种电力系统状态估计的量测量权重设置方法，进一步提高电力系统状态估计的精度。

为了实现上述目的，本发明的技术方案如下：一种电力系统状态估计的量测量权重设置方法，将节点注入功率量测量的类型因子设置为支路功率量测量的类型因子除以该节点的出线数。具体包括以下步骤：

A、根据支路数据统计各节点的出线数α_j；

B、读取量测数据z；

C、设置当前序号i＝1；

D、判断当前量测量是否为电压幅值量测量，如果不是电压幅值量测量转至步骤F；

E、设置当前量测量的类型因子为k_wi＝a，取a＝0.8～1.2；

F、判断当前量测量是否为支路功率量测量，如果不是支路功率量测量转至步骤H；

G、设置当前量测量的类型因子为k_wi＝c，取c＝8.5～11.5；

H、判断当前量测量是否为节点注入功率量测量，如果不是节点注入功率量测量转至步骤J；

I、取出节点注入功率量测量对应节点的出线数α_j，并设置当前量测量的类型因子为k_wi＝c/α_j；

J、设置当前量测量权重为

K、令i＝i+1；

L、判断i是否大于量测量个数m，如果i不大于m，则转至步骤D；否则结束。

与现有技术相比，本发明具有以下有益效果：

本发明设置权重时同时考虑量测量误差和量测量类型，节点注入功率量测量的类型因子k_w为支路功率量测量的类型因子除以该节点的出线数。实施算例表明本发明提出的量测量权重设置方法可以明显提高电力系统状态估计的精度。对IEEE30系统，各量测量均存在1％误差时，采用本发明方法设置量测量权重，状态估计的精度比传统方法提高7.4％，比专利ZL201310314260.1方法提高2.7％；各量测量均存在5％误差时，采用本发明方法设置量测量权重，状态估计的精度比传统方法提高7.8％，比专利ZL201310314260.1方法提高1.7％。

附图说明

本发明共有附图4张，其中：

图1是现有技术的加权最小二乘法状态估计流程图。

图2是节点注入功率量测量权重设置方法的证明用图。

图3是本发明的量测量权重设置流程图。

图4是本发明的量测量权重设置算例IEEE30系统接线图。

具体实施方式

下面结合附图对本发明作进一步地说明。

结合图2节点注入功率量测量权重设置方法的证明用图，对本发明的节点注入功率量测量权重设置方法的依据论证如下：

图2中节点i通过支路l_ij、l_ik、l_il分别与节点j、k、l三个节点相连，即节点i的出线数为3。为了论述方便，假设量测量P_i、P_ij、P_ik、P_il的编号排在前面，为前4个量测量；状态量的修正量ΔV_i、ΔV_j、ΔV_k、ΔV_l、Δθ_i、Δθ_j、Δθ_k、Δθ_l的编号也排在前面，为前8个状态量的修正量；并设各量测量的方差都相同，所以权重也相同，这里为分析方便设R^-1＝1。

节点注入有功功率与支路有功功率存在如下关系：

P_i＝P_ij+P_ik+P_il (12)

状态估计的法方程为

H^TR^-1HΔx＝H^TR^-1r (13)

由于设R^-1＝1，式(13)简化为

H^THΔx＝H^Tr (14)

下面分别分析节点i的注入有功功率量测量的误差对状态估计的影响，支路l_ij、l_ik、l_il的支路有功功率量测量的误差对状态估计的影响。假设原来的量测量没有误差，状态量已经计算出来，定义这时的量测量和状态量数据为基态，此时式(14)中残差向量r的元素都为0，信息矩阵H^TH是个定值。分析量测量误差影响时，以此基态为基础。如果某些量测量出现了误差，式(14)中残差向量r的对应元素不为0，信息矩阵H^TH不变，还是原来的值，右端向量的数值则会发生变化。下面考察一下式(14)的右端向量的数值变化情况。

(1)节点i的注入有功功率量测量的误差对状态估计的影响

假设各量测量原来均无误差，当只有一个节点i的注入有功功率量测量存在误差r₀时，分析法方程的右端向量的值。式(14)中除r₁≠0外，其他量测量的残差都为0，法方程的右端向量的值为

$H^{T}r=\left[\begin{array}{ccccc}\∂P_i/\∂V_i& \∂P_{\mathrm{ij}}/\∂V_i& \∂P_{\mathrm{ik}}/\∂V_i& \∂P_{\mathrm{il}}/\∂V_i& \·\·\·\\ \∂P_i/\∂V_j& \∂P_{\mathrm{ij}}/\∂V_j& 0& 0& \·\·\·\\ \∂P_i/\∂V_k& 0& \∂P_{\mathrm{ik}}/\∂V_k& 0& \·\·\·\\ \∂P_i/\∂V_l& 0& 0& \∂P_{\mathrm{il}}/\∂V_l& \·\·\·\\ \∂P_i/\∂{\mathrm{\θ}}_i& \∂P_{\mathrm{ij}}/\∂{\mathrm{\θ}}_i& \∂P_{\mathrm{ik}}/\∂{\mathrm{\θ}}_i& \∂P_{\mathrm{il}}/\∂{\mathrm{\θ}}_i& \·\·\·\\ \∂P_i/\∂{\mathrm{\θ}}_j& \∂P_{\mathrm{ij}}/\∂{\mathrm{\θ}}_j& 0& 0& \·\·\·\\ \∂P_i/\∂{\mathrm{\θ}}_k& 0& \∂P_{\mathrm{ik}}/\∂{\mathrm{\θ}}_k& 0& \·\·\·\\ \∂P_i/\∂{\mathrm{\θ}}_l& 0& 0& \∂P_{\mathrm{il}}/\∂{\mathrm{\θ}}_l& \·\·\·\\ \·& & & & \\ \·& 0& 0& 0& \·\·\·\\ \·& & & & \end{array}\right]\left[\begin{array}{c}{r}_0\\ 0\\ 0\\ 0\\ 0\\ 0\\ 0\\ 0\\ 0\end{array}\right]=r_0\left[\begin{array}{c}\∂P_i/\∂V_i\\ \∂P_i/\∂V_j\\ \∂P_i/\∂V_k\\ \∂P_i/\∂V_l\\ \∂P_i/\∂{\mathrm{\θ}}_i\\ \∂P_i/\∂{\mathrm{\θ}}_j\\ \∂P_i/\∂{\mathrm{\θ}}_k\\ \∂P_i/\∂{\mathrm{\θ}}_l\\ 0\end{array}\right]---\left(15\right)$

由式(15)可见，当节点i的注入有功功率量测量存在误差时，法方程的右端向量中有8个元素不为0，分别与本节点的有功功率函数对本节点及相连节点的电压幅值和电压相角的偏导有关。

(2)支路l_ij的支路有功功率量测量的误差对状态估计的影响

假设各量测量原来均无误差，当从节点i出发的1条支路l_ij的支路有功功率量测量存在误差r₀时，分析法方程的右端向量的值。式(14)中除r₂≠0外，其他量测量的残差都为0，法方程的右端向量的值为

$H^{T}r=\left[\begin{array}{ccccc}\∂P_i/\∂V_i& \∂P_{\mathrm{ij}}/\∂V_i& \∂P_{\mathrm{ik}}/\∂V_i& \∂P_{\mathrm{il}}/\∂V_i& \·\·\·\\ \∂P_i/\∂V_j& \∂P_{\mathrm{ij}}/\∂V_j& 0& 0& \·\·\·\\ \∂P_i/\∂V_k& 0& \∂P_{\mathrm{ik}}/\∂V_k& 0& \·\·\·\\ \∂P_i/\∂V_l& 0& 0& \∂P_{\mathrm{il}}/\∂V_l& \·\·\·\\ \∂P_i/\∂{\mathrm{\θ}}_i& \∂P_{\mathrm{ij}}/\∂{\mathrm{\θ}}_i& \∂P_{\mathrm{ik}}/\∂{\mathrm{\θ}}_i& \∂P_{\mathrm{il}}/\∂{\mathrm{\θ}}_i& \·\·\·\\ \∂P_i/\∂{\mathrm{\θ}}_j& \∂P_{\mathrm{ij}}/\∂{\mathrm{\θ}}_j& 0& 0& \·\·\·\\ \∂P_i/\∂{\mathrm{\θ}}_k& 0& \∂P_{\mathrm{ik}}/\∂{\mathrm{\θ}}_k& 0& \·\·\·\\ \∂P_i/\∂{\mathrm{\θ}}_l& 0& 0& \∂P_{\mathrm{il}}/\∂{\mathrm{\θ}}_l& \·\·\·\\ \·& & & & \\ \·& 0& 0& 0& \·\·\·\\ \·& & & & \end{array}\right]\left[\begin{array}{c}0\\ r_0\\ 0\\ 0\\ 0\\ 0\\ 0\\ 0\\ 0\end{array}\right]=r_0\left[\begin{array}{c}\∂P_{\mathrm{ij}}/\∂V_i\\ \∂P_{\mathrm{ij}}/\∂V_j\\ 0\\ 0\\ \∂P_{\mathrm{ij}}/\∂{\mathrm{\θ}}_i\\ \∂P_{\mathrm{ij}}/\∂{\mathrm{\θ}}_j\\ 0\\ 0\\ 0\end{array}\right]---\left(16\right)$

由式(16)可见，当支路l_ij的支路有功功率量测量存在误差时，法方程的右端向量中有4个元素不为0，分别与支路有功功率函数对支路两端节点的电压幅值和电压相角的偏导有关。

对支路l_ik、l_il的分析也能得出类似的结论。

(3)支路l_ij、l_ik、l_il的支路有功功率量测量的误差对状态估计的影响

假设各量测量原来均无误差，当从节点i出发的3条支路的支路有功功率量测量都存在相同误差r0时，分析法方程的右端向量的值。式(14)中除r₂≠0、r₃≠0、r₄≠0外，其他量测量的残差都为0，法方程的右端向量的值为

$H^{T}r=\left[\begin{array}{ccccc}\∂P_i/\∂V_i& \∂P_{\mathrm{ij}}/\∂V_i& \∂P_{\mathrm{ik}}/\∂V_i& \∂P_{\mathrm{il}}/\∂V_i& \·\·\·\\ \∂P_i/\∂V_j& \∂P_{\mathrm{ij}}/\∂V_j& 0& 0& \·\·\·\\ \∂P_i/\∂V_k& 0& \∂P_{\mathrm{ik}}/\∂V_k& 0& \·\·\·\\ \∂P_i/\∂V_l& 0& 0& \∂P_{\mathrm{il}}/\∂V_l& \·\·\·\\ \∂P_i/\∂{\mathrm{\θ}}_i& \∂P_{\mathrm{ij}}/\∂{\mathrm{\θ}}_i& \∂P_{\mathrm{ik}}/\∂{\mathrm{\θ}}_i& \∂P_{\mathrm{il}}/\∂{\mathrm{\θ}}_i& \·\·\·\\ \∂P_i/\∂{\mathrm{\θ}}_j& \∂P_{\mathrm{ij}}/\∂{\mathrm{\θ}}_j& 0& 0& \·\·\·\\ \∂P_i/\∂{\mathrm{\θ}}_k& 0& \∂P_{\mathrm{ik}}/\∂{\mathrm{\θ}}_k& 0& \·\·\·\\ \∂P_i/\∂{\mathrm{\θ}}_l& 0& 0& \∂P_{\mathrm{il}}/\∂{\mathrm{\θ}}_l& \·\·\·\\ \·& & & & \\ \·& 0& 0& 0& \·\·\·\\ \·& & & \end{array}\right]\left[\begin{array}{c}0\\ r_0\\ r_0\\ r_0\\ 0\\ 0\\ 0\\ 0\\ 0\end{array}\right]=r_0\left[\begin{array}{c}\∂P_{\mathrm{ij}}/\∂V_i+\∂P_{\mathrm{ik}}/\∂V_i+\∂P_{\mathrm{il}}/\∂V_i\\ \∂P_{\mathrm{ij}}/\∂V_j\\ \∂P_{\mathrm{ik}}/\∂V_k\\ \∂P_{\mathrm{il}}/\∂V_l\\ \∂P_{\mathrm{ij}}/\∂{\mathrm{\θ}}_i+\∂P_{\mathrm{ik}}/\∂{\mathrm{\θ}}_i+\∂P_{\mathrm{il}}/\∂{\mathrm{\θ}}_i\\ \∂P_{\mathrm{ij}}/\∂{\mathrm{\θ}}_j\\ \∂P_{\mathrm{ik}}/\∂{\mathrm{\θ}}_k\\ \∂P_{\mathrm{il}}/\∂{\mathrm{\θ}}_l\\ 0\end{array}\right]---\left(17\right)$

由式(17)可见，当从节点i出发的3条支路l_ij、l_ik、l_il的支路有功功率量测量都存在相同误差时，法方程的右端向量中有8个元素不为0，分别与这些支路有功功率函数对支路两端节点的电压幅值和电压相角的偏导有关。

(4)误差对状态估计影响的分析结论

由式(12)可知

$\left\{\begin{array}{c}\∂P_i/\∂V_i=\∂P_{\mathrm{ij}}/\∂V_i+\∂P_{\mathrm{ik}}/\∂V_i+\∂P_{\mathrm{il}}/\∂V_i\\ \∂P_i/\∂V_j=\∂P_{\mathrm{ij}}/\∂V_j\\ \∂P_i/\∂V_k=\∂P_{\mathrm{ik}}/\∂V_k\\ \∂P_i/\∂V_l=\∂P_{\mathrm{il}}/\∂V_l\\ \∂P_i/\∂{\mathrm{\θ}}_i=\∂P_{\mathrm{ij}}/\∂{\mathrm{\θ}}_i+\∂P_{\mathrm{ik}}/\∂{\mathrm{\θ}}_i+\∂P_{\mathrm{il}}/\∂{\mathrm{\θ}}_i\\ \∂P_i/\∂{\mathrm{\θ}}_j=\∂P_{\mathrm{ij}}/\∂{\mathrm{\θ}}_j\\ \∂P_i/\∂{\mathrm{\θ}}_k=\∂P_{\mathrm{ik}}/\∂{\mathrm{\θ}}_k\\ \∂P_i/\∂{\mathrm{\θ}}_l=\∂P_{\mathrm{il}}/\∂{\mathrm{\θ}}_l\end{array}\right.---\left(18\right)$

对比式(15)和式(16)，考虑式(18)，可知式(16)的非零元素比式(15)多，且式(16)的第1列和第4列的元素也比式(15)的对应元素大。说明在相同误差的情况下，节点注入有功功率量测量的误差对状态估计精度的影响比支路有功功率量测量的误差对状态估计精度的影响大。

对比式(15)和式(17)，考虑式(18)，可知式(15)和式(17)是相等的，说明节点i的注入有功功率量测量的误差对状态估计的影响与从该节点出发的所有支路的有功功率量测量存在相同误差时对状态估计的影响是相同的。

节点注入有功功率量测量的误差对状态估计的影响与节点的出线数有关，出线越多，影响越大。因此可以根据节点的出线数设置节点注入有功功率量测量的权重。

同理节点注入无功功率量测量也可得出相同的结论。

下面采用图1和图3所示的算法对图4所示的IEEE30系统进行状态估计。计算时以潮流计算结果为量测量的真值和状态量的真值，采用3种方法分别针对5种工况进行状态估计。专利ZL201310314260.1的量测量类型因子均取典型值，即a＝1，b＝3，c＝10；本发明的量测量类型因子也取典型值，即a＝1，c＝10，节点j的注入功率量测量的类型因子为k_wi＝c/α_j，分析时数据单位取标幺值。

5种工况分别如下：

工况1：所有量测量增加各自量测量的1％，即所有量测量都存在1％误差；

工况2：所有量测量增加各自量测量的2％，即所有量测量都存在2％误差；

工况3：所有量测量增加各自量测量的3％，即所有量测量都存在3％误差；

工况4：所有量测量增加各自量测量的4％，即所有量测量都存在4％误差；

工况5：所有量测量增加各自量测量的5％，即所有量测量都存在5％误差。

3种方法分别如下：

方法1：传统的以量测量方差倒数作为量测量权重的方法；

方法2：专利ZL201310314260.1的设置量测量权重方法；

方法3：本发明提出的设置量测量权重方法。

计算结果如表1，表中状态估计总误差为状态估计误差平方和，即其中z_0i为第i个量测量的真值；h_i(x)为第i个量测量的计算值。本发明方法精度提高率的计算规则：与方法1比较时以方法1的状态估计总误差为基准，与方法2比较时则以方法2的状态估计总误差为基准，计算结果见表1。

表1 状态估计结果表

表1可见，与传统的量测量权重设置方法相比，本发明提出的量测量权重设置方法的状态估计的精度明显提高；与专利ZL201310314260.1的量测量权重设置方法相比，本发明提出的量测量权重设置方法的状态估计的精度也有一定的提高。

本发明可以采用任何一种编程语言和编程环境实现，如C语言、C++、FORTRAN、Delphi等。开发环境可以采用VisualC++、BorlandC++Builder、VisualFORTRAN等。

Claims (1)

1.一种电力系统状态估计的量测量权重设置方法，其特征在于：包括以下步骤：

A、根据支路数据统计各节点的出线数α_j；

B、读取量测数据z；

C、设置当前序号i＝1；

D、判断当前量测量是否为电压幅值量测量，如果不是电压幅值量测量转至步骤F；

E、设置当前量测量的类型因子为k_wi＝a，取a＝0.8～1.2；

F、判断当前量测量是否为支路功率量测量，如果不是支路功率量测量转至步骤H；

G、设置当前量测量的类型因子为k_wi＝c，取c＝8.5～11.5；

H、判断当前量测量是否为节点注入功率量测量，如果不是节点注入功率量测量转至步骤J；

I、取出节点注入功率量测量对应节点的出线数α_j，并设置当前量测量的类型因子为k_wi＝c/α_j；

J、设置当前量测量权重为

K、令i＝i+1；

L、判断i是否大于量测量个数m，如果i不大于m，则转至步骤D；否则结束。