[发明专利]考虑零注入约束的修正快速解耦电力系统状态估计方法

权利要求书

1.一种考虑零注入约束的修正快速解耦电力系统状态估计方法，其特征在于该方法包括以下步骤：

(1)建立一个含等式约束的电力状态估计模型：

min12Σi=1m(zi-hi(x))2]]>

s.t        c(x)＝0

等式约束为：使零注入节点的节点注入功率为0，以c(x)＝0表示，其中x是电力系统的状态变量，包括零注入节点的电压幅值向量V_z、零注入节点的电压相角向量θ_z、非零注入节点的电压幅值向量V_n和非零注入节点的电压相角向量θ_n，z_i是第i号量测的实时测量值，h_i(x)是第i号量测的实时测量方程，m为量测个数；

(2)根据电力系统当前的拓扑结构和网络参数，形成电力系统当前的节点导纳矩阵，并计算以下系数矩阵F

F=-BzzGzzGzz-Bzz-1BznGznGzn-Bzn]]>

其中矩阵G_zz和B_zz分别是节点导纳矩阵中零注入节点对应的对角子阵的实部和虚部，矩阵G_zn和B_zn分别是节点导纳矩阵中零注入节点与非零注入节点的交叉非对角子阵对应的实部和虚部；

(3)根据电力系统当前的拓扑结构和网络参数，形成快速解耦估计有功迭代的系数矩阵A和无功迭代的系数矩阵R，分别为：

A=HaTWaHa]]>

R=HrTWrHr]]>

其中H_a，H_r分别是电力系统状态估计的有功迭代和无功迭代中使用的定常量测雅可比矩阵，W_a，W_r分别为电力系统状态估计的有功量测和无功量测的权重矩阵，上标T表示转置；

(4)设置电力系统状态估计的节点电压幅值的初值，记为V⁽⁰⁾，节点电压相角的初值为θ⁽⁰⁾，并设置迭代次数，记为k＝0；

(5)第k次迭代得到节点电压幅值向量V^(k)和电压相角向量θ^(k)，将θ^(k)中零注入节点的电压相角向量和非零注入节点的电压相角向量分别记为和将V^(k)中零注入节点的电压幅值向量和非零注入节点的电压幅值向量分别记为和

(6)保持非零注入节点的相角子向量不变，重新计算零注入节点对应的电压幅值向量相角子向量计算方法如下：

通过下式计算零注入节点的状态变量向量

xz(k)=ΦFΦ-1xn(k)]]>

其中矩阵F为步骤(2)中的计算结果，Φ为复数理论中用直角坐标表示的复数到用极坐标表示的复数的转换映射；Φ^-1表示Φ的逆映射。是非零注入节点的状态变量，由非零注入节点的电压幅值向量和相角向量组成，计算得到的是零注入节点的状态变量向量，由零注入节点的电压幅值向量和相角向量组成，由向量就可以得到和

Φ的表达式是：

V=e2+f2]]>

θ=arctan(ef)]]>

Φ^-1的表达式是：

e＝Vcosθ

f＝Vsinθ

其中e，f是用直角坐标表示的节点电压的实部和虚部，V，θ是用极坐标表示的节点电压的幅值和相角；。

(7)根据当前的零注入节点电压幅值向量和相角向量以及非零注入节点的电压幅值向量和相角向量计算第k次迭代的有功迭代的右手向量α^(k)，计算公式为：

α(k)=-HaTWar(V~z(k),θ~z(k),Vn(k),θn(k))]]>

其中r为状态估计第k次迭代的残差向量，可由当前的零注入节点电压幅值向量和相角向量以及非零注入节点的电压幅值向量和相角向量计算得到；

(8)求解如下的快速解耦估计有功迭代修正方程，得到第k次迭代节点相角向量的修正量Δθ^(k)

AΔθ^(k)＝α^(k)

(9)根据设定的收敛精度ε，对上述修正量Δθ^(k)进行判断，若max|Δθ^(k)|＜ε，且无功迭代已收敛，则状态估计计算完成；若max|Δθ^(k)|＜ε，且无功迭代未收敛，则标记当前有功迭代已收敛，并进行步骤(12)；若max|Δθ^(k)|≥ε，则进行步骤(10)；通常取ε＝0.0001；

(10)根据第k次迭代的节点电压相角修正量Δθ^(k)，将Δθ^(k)中零注入节点的电压相角修正量和非零注入节点的电压相角修正量分别记为和得到第k+1次迭代的电压相角变量为：θn(k+1)=θn(k)+Δθn(k);]]>

(11)根据上述第k+1次迭代的电压相角变量利用步骤(6)的方法，重新计算零注入节点对应的电压幅值向量相角子向量，记为和

(12)根据当前的零注入节点电压幅值向量和相角向量以及非零注入节点的电压幅值向量和相角向量计算第k次迭代的无功迭代的右手向量β^(k)，计算公式为：

β(k)=-HrTWrr(V~z(k)′,θ~z(k)′,Vn(k),θn(k+1))]]>

(13)求解以下的快速解耦估计无功迭代修正方程，得到第k次迭代节点电压幅值向量的修正量ΔV^(k)

RΔV^(k)＝β^(k)；

(14)根据设定的收敛精度ε，对上述修正量ΔV^(k)进行判断，若max|ΔV^(k)|＜ε，且有功迭代已收敛，则状态估计计算完成；若max|ΔV^(k)|＜ε，且有功迭代未收敛，则标记当前无功迭代已收敛，并进行步骤(15)，若max|ΔV^(k)|≥ε，则进行步骤(5)；

(15)第k次迭代的节点电压幅值修正量为ΔV^(k)，将ΔV^(k)中零注入节点的电压幅值修正量和非零注入节点的电压幅值修正量分别记为和k＝k+l，进行步骤(5)。