1.一种SOA架构的电网实时抗差状态估计方法，采用一种SOA架构的电网实时抗差状态估计系统，系统包括部署在云计算平台上对状态估计系统进行设置和结果显示的用户端，用户端与电网系统相连接，接收来自电力系统的量测值和网架拓扑信息，部署在云计算平台上用于实现状态估计逻辑计算的服务端，服务端包括根据电网量测值和网架拓扑信息进行状态估计计算的处理单元，以及供用户端接收数据的服务接口单元，用户端与服务接口通过网络相连接，其特征是包括以下步骤：S1.用户端获取电网量测值和网架拓扑信息；

S2.根据用户对量测值的选择调整作为状态估计计算的原始数据；

S3.根据原始数据分析电网在电气上是否为单一区域运行情况，若是，则直接将原始数据处理生成一个区域的标准格式数据文件，若不是，则按照多区域处理办法生成合并计算的标准格式数据文件；计算标准格式数据文件的具体过程包括：S31.根据网架拓扑连接关系识别电网在电气关系上是一个电气区域，还是相互独立的多个电气区域，即是否为多区域独立运行情况；

S32.若是一个区域，则按常规方法将原始数据处理为标准格式数据文件，进入步骤S4；

若为多个区域，则将每个区域的原始数据分别处理为标准格式数据文件；

S33.以节点最多的电气区域为基础，保持该区域所有节点的编号不变，其他区域的节点号依次按顺序递增，并重新按照标准格式数据文件进行处理，直到将所有区域的标准格式数据合并为一个多平衡节点的标准格式数据文件；

S34.合并后，按常规方法将原始数据处理为标准格式数据文件；

S4.将标准格式的数据文件发送给服务端；

S5.服务端读取数据文件中的数据进行状态估计计算，计算结果供用户端进行调用。

2.根据权利要求1所述的一种SOA架构的电网实时抗差状态估计方法，其特征是步骤S1中获取的电网量测值包括节点注入功率量测值、节点电压量测值、支路功率量测值、零注入功率量测值；网架拓扑信息包括电网节点数、平衡节点号、节点连接关系、线路和变压器支路阻抗参数、接地支路参数。

3.根据权利要求1所述的一种SOA架构的电网实时抗差状态估计方法，其特征是步骤S2的具体过程包括：

S21.将电网中所有可以选择的量测数据点在用户端上显示并供用户选择；

S22.根据用户选择的量测数据点并结合默认选择的节点量测数据点，与上一次进行状态估计计算的量测数据点进行比较；

S23.将上一次进行状态估计计算的量测数据点作为原始数据，若比较为增加了量测数据点，将增加的数据点添加到原始数据中，若比较为减少了量测数据点，将减少的数据点从原始数据中删除，若比较为没有改变，则保持原始数据量测点不变。

4.根据权利要求1或2所述的一种SOA架构的电网实时抗差状态估计方法，其特征是步骤5中状态估计计算的具体过程包括：

S51.读取数据文件，根据状态估计的非线性量测方程z＝g(x)+ε

采用最小绝对值法进行计算，最小绝对值模型为其中z为m维量测向量，即三类量测值中用于计算的m各量测值组成的向量，z＝[Pi,Qi；

Vi；Pij,Qij]，Pi、Qi分别为节点注入有功率和无功功率量测值，Vi为节点电压幅值量测值，Pij、Qij为支路功率量测值；x为n维状态向量，即x＝[V,δ]，V、δ分别为电网节点电压幅值和相角，共n个；g(x)为m维量测函数向量，为量测值z中m个量测值对应的m个函数表达式向量，即为估计值， 为节点注入功率估计值， 为节点电压估计值， 为支路功率估计值；ε为m维量测残差，h(x)为节点零节点注入功率量测值，w为m维权重向量；T为矩阵运算符号，表示对向量的转置；

S52.初始化参数，设定迭代次数k＝0，kmax＝20中心参数σ＝0.1，l＝u＝λ＝[1,…,1]T∈Rm，y＝0，α＝β＝w，收敛精度e＝10-5，状态变量x采用平启动，即V＝1，δ＝0，其中l,u为松弛变量，(l,u)∈Rm≥0，y、λ、α、β为拉格朗日乘子向量，y∈Rc，λ、α、β∈Rm，Rc、Rm分别为c维和m维实数向量；

S53.判断是否k

S54.计算互补间隙Gap＝αTl+βTu，若Gap

S55.计算扰动因子

S56.根据最小绝对值模型获得修正方程，求解修正方程得到修正量Δx，Δλ，Δy，Δα，Δβ，Δl，Δu；

S57.求出原始对偶步长

i＝1,…,m；

S58.修正原始、对偶变量

将k＝k+1，返回步骤S53 。

5.根据权利要求1所述的一种SOA架构的电网实时抗差状态估计方法，其特征是步骤S56 修正量获取的过程包括：

S561.等价转化最小绝对值模型，获得最小值等价模型引入松弛变量(l,u)∈Rm≥0将不等式转化为等式，并进行等价转换，得到S562.构建拉格朗日函数：

L＝wT(l+u)-yTh(x)-λT(z-g(x)+l-u)-αTl-βTu  (3)求取拉格朗日函数极限，导出下一阶最优性条件：Ly≡h(x)＝0     (5)

Lλ≡z-g(x)+l-u＝0     (6)Ll≡w-λ-α＝0    (7)

Lu≡w+λ-β＝0   (8)

Lα≡diag(α)diag(l)e＝0   (9)Lβ≡diag(β)diag(u)e＝0   (10)其中，e＝[1,…,1]T∈Rm，引入扰动参数μ＞0，对互补项(9)、(10)实施扰动松弛得到：S563.对以上非线性方程采用牛顿法求解，线性化得到：-Δλ-Δα＝-Ll    (16)

Δλ-Δβ＝-Lu    (17)

S564.变量α，β，λ选取合适初始值，使得Ll＝Lu＝0，由公式(16)、(17)、(18)、(19)求得以下关系式：求得Δα，Δβ，Δl，Δu；

S565.将Δl、Δu代入公式(15)，得到：公式中S＝diag(α)-1diag(l)+diag(β)-1diag(u)；

由公式(13)、(14)、(21)得到修正方程：公式中

求解修正方程得到Δx，Δλ，Δy。