[发明专利]一种基于最大指数绝对值目标函数的抗差状态估计方法

说明书

技术领域

本发明属于电力系统调度自动化领域，具体涉及一种基于最大指数绝对值目标函数的抗差状态估计方法。

背景技术

电力系统状态估计是能量管理系统的基础和核心。现在几乎每一个大型调度中心都安装了状态估计器，状态估计已成为电网安全运行的基石。自1970国外学者首次提出状态估计以来，人们对状态估计的研究和应用已经有40多年的历史了，这期间涌现出了各种各样的状态估计方法。

目前，在国内外应用最为广泛的状态估计是加权最小二乘法(Weighted least squares,WLS)。WLS模型简洁，求解容易，但是其抗差性很差。为了增强抗差性，一般有两种方法。第一种是在WLS估计之后加入不良数据辨识环节，例如最大正则化残差检验法(LNR)或估计辨识方法等；另一种是采用抗差状态估计方法(Robust state estimation)。目前，国内外学者已经提出的抗差状态估计方法包括加权最小绝对值估计(Weighted least absolute value,WLAV)、非二次准则法(QL、QC等)、以合格率最大为目标的状态估计(Maximum normal measurement rate,MNMR)以及指数型目标函数状态估计(Maximum exponential square,MES)等。但是这些抗差状态估计方法普遍存在计算效率不够高的特点，从而在一定程度上影响了它们在实际系统中的应用。

发明内容

本发明旨在至少在一定程度上解决上述技术问题之一或至少提供一种有用的商业选择。为此，本发明的一个目的在于提出一种抗差性好、计算效率高的基于最大指数绝对值目标函数的抗差状态估计方法(Maximum exponential absolute value state estimation,MEAV)。

根据本发明实施例的基于最大指数绝对值目标函数的抗差状态估计方法，包括步骤：步骤A.提供基于最大指数绝对值目标函数的抗差状态估计基本模型；步骤B.对所述基于最大指数绝对值目标函数的抗差状态估计基本模型引进辅助变量，变换得到基于最大指数绝对值目标函数的抗差状态估计等价模型；以及步骤C.利用原-对偶内点算法，对所述基于最大指数绝对值目标函数的抗差状态估计等价模型求解。

在本发明的一个实施例中，所述基于最大指数绝对值目标函数的抗差状态估计基本模型为：s.t.g(x)=0，r=z-h(x)，其中：z∈R^m为量测矢量，包括节点注入有功和无功、支路有功和无功以及节点电压幅值量测；x∈Rⁿ为状态矢量，包括节点电压幅值和平衡节点除外的其他各个节点相角；h:Rⁿ→R^m为由状态矢量到量测矢量的非线性映射；r_i为残差矢量r的第i个元素；g(x):Rⁿ→R^c为零注入功率等式约束；w_i为第i个量测量的权重，同是为窗宽参数。

在本发明的一个实施例中，所述步骤B包括：引进非负松弛变量u,v∈R^m，变换得到的所述基于最大指数绝对值目标函数的抗差状态估计等价模型为：s.t.g(x)=0，z-h(x)-u+v＝0，u,v≥0。