一种电力系统状态的双曲余弦型抗差状态估计方法
【技术领域】
[0001]本发明属于电力系统调度自动化领域，特别涉及一种电力系统状态的双曲余弦型 抗差状态估计方法。
【背景技术】
[0002] 电力系统状态估计是能量管理系统的基础和核心。现在几乎每一个大型调度中心 都安装了状态估计器，状态估计已成为电网安全运行的基石。自1970国外学者首次提出状 态估计以来，人们对状态估计的研究和应用已经有40多年的历史了，这期间涌现出了各种 各样的状态估计方法。
[0003] 目前，在国内外应用最为广泛的状态估计是加权最小二乘法(Weighted least squares JLShWLS模型简洁，求解容易，但是其抗差性很差。为了增强抗差性，一般有两种 方法。第一种是在WLS估计之后加入不良数据辨识环节，例如最大正则化残差检验法(LNR) 或估计辨识方法等;另一种是采用抗差状态估计方法。目前，国内外学者已经提出的抗差状 态估计方法（Robust state estimation)包括加权最小绝对值估计(Weighted least absolute value，WLAV)、非二次准则法（QL、QC等）、以合格率最大为目标的状态估计 (Maximum normal measurement rate, MM1R)以及指数型目标函数状态估计（Maximum exponential square，MES)等。但是这些抗差状态估计方法的估计性能仍有待提高。
【发明内容】
[0004] 本发明的目的是提出一种电力系统状态的双曲余弦型抗差状态估计方法，其特征 在于，该方法是基于抗差性好、计算效率高的双曲余弦型抗差状态估计方法;包括如下步 骤：
[0005] 步骤A.提取电力系统节点注入有功和无功、支路有功和无功功率，以及节点电压 幅值参数;以此建立双曲余弦型抗差状态估计模型；
[0006] 步骤B.利用原-对偶内点算法，对所述双曲余弦型抗差状态估计模型求解。
[0007] 所述步骤A中建立双曲余弦型抗差状态估计模型为：
s.t. g(x) = 0，r = z-h(x)，其中：zeRm为量测矢量，包括节点注入有功和无功、支路有功和 无功以及节点电压幅值量测;xeRn为状态矢量，包括节点电压幅值和平衡节点除外的其他 各个节点相角;h:R n-Rm为由状态矢量到量测矢量的非线性映射;Γι为残差矢量r的第i个元 素;为零注入功率等式约束; Wi为第i个量测量的权重，σ为窗宽参数。
[0008] 所述步骤Β为利用原-对偶内点算法，对双曲余弦型抗差状态估计模型求解，包括：
[0009] 步骤Β1:令X为平启动状态变量;选择= = 设置收敛判据ε = 1〇Λ置 迭代计数器k=0;
[0010] 步骤B2:求解修正方程，得到[dxT drT (1λτ d3TT];
[0011] 步骤B3,修正变i
[0012] 步骤Μ:判断是否收敛，若max(dx) < ε，则转步骤B6，否则进入步骤B5;
[0013] 步骤Β5:令迭代计数器k = k+l，进入步骤Β2;
[0014] 步骤B6:输出最优解，结束。
[0015] 所述步骤B2包括：
[0 016 ] 步骤B 21:形成量测方程以及零注入功率约束对应的雅克比矩阵Η = ?⑷/ &及 (? = %(_/&，形成量测方程以及零注入功率约束对应的海森矩阵V2h(x)及V2g( x)，其中h (X)为状态矢量到量测矢量的映射，即为量测估计值;z = h(X)+e为量测方程，z和X如步骤A 所示，ε eRm为量测误差向量:h(x)的具体表达式介绍如下:对于节点i的电压幅值量测，Vi = Vl;对于注入有功量测，
}意思为与节点i相联的所有 节点;Gij，Bij分别为节点导纳矩阵第i行第j列的实部和虚步.V、和占 i和书占 电压幅值;θυ为节点i和节点j的相角差;对于注入无功量测
对于支路i j的有功量测，今=_vf(I + %) - ^名一08- νΛ' \ sin ，其中gsi为支路i j的首端对 地电导，其中gij和bij分别为支路i j的串联电导和串联电纳；对于支路i j的无功量测， % = + \) + v, VA. cos & - v, v,.gy. sin &，其中bsi为支路i j的首端对地电纳;g ( X ) = 0为零注 入功率约束；
[0017]步骤B22:引入拉格朗日函数 [0018]
[0019 ] 式中：λ e Rc及JT e Rm为拉格朗日乘子矢量；
[0020]计算
，
[0021 ] W为RmXm的对角阵，其对角元素为
[0022] 步骤B23:求解方i
得到[dxT drT (1λτ dJiT] 0
[0023] 本发明的有益效果是双曲余弦型抗差状态估计方法在估计过程中可有效抑制包 括一致性不良数据在内的多个不良数据，显示了良好的抗差性，并具有很高的计算效率，非 常适宜于实际工程应用。
【具体实施方式】
[0024]本发明提出一种电力系统状态的双曲余弦型抗差状态估计方法，下面结合实施例 详细描述本发明。
[0025] 所述电力系统状态的双曲余弦型抗差状态估计方法(Hyperbolic cosine state estimation, COSH)包括下列步骤：
[0026] 步骤A:提供双曲余弦型抗差状态估计（Hyperbolic cosine state estimation, COSH)模型。
[0027] 具体地，本发明提出的COSH的模型如下所示
[0028]
(1)
[0029] s.t. g(x)=0 (2)
[0030] r = z~h(x) (3)
[0031] 式中：ZeRm为量测矢量，常包括节点注入有功和无功、支路有功和无功以及节点 电压幅值量测等;为包括节点电压幅值和相角的状态矢量(平衡节点相角除外）；h:R n -Rm为由状态矢量到量测矢量的非线性映射;η是残差矢量r的第i个元素;gUhRn-F为 零注入功率等式约束;wi为第i个量测量的权重， σ为窗宽参数。
[0032] 步骤Β:利用原-对偶内点算法，对所述双曲余弦型抗差状态估计模型求解。
[0033] (1)C0SH模型的求解方法
[0034]注意到C0SH模型（1)~(3)是一个含有等式约束和不等式约束的最优化问题，适宜 用原-对偶内点算法进行求解。为使本领域技术人员更好地理解本发明，首先给出详细的推 导过程如下：
[0035]引入拉格朗日函数
[0036]
(4)
[0037] 式中：为拉格朗日乘子矢量。
[0038]为取得最优值，根据KKT条件，可得
[0039]
[0040]
[0041]
[0042]
[0043]
[0044] 以上方程由牛顿法求解可得
[0045] [ V2g (X) λ-V2h (X) π ] dx+GTdA-HTd3i = -Lx (9)
[0046] Gdx = -Lx (10)
[0047] -Hdx-dr = -Li (11)
[0048]
(12)
[0049] 式（12)可表示为矩阵形式，为
[0050] Wdr+d3i = -Lr (13)
[0051] 式中：W为RmXm的对角阵，其对角元素为
[0052]根据式(9)、（10)、（11)及（13)，可得修正方程为
[0053]
(14)
[0054] 求解式（14)可得[dxT drT (1λτ dJiT]，则迭代即可持续进行。
[0055] (2)C0SH模型的求解步骤
[0056]在介绍C0SH模型的求解推导过程之