1.一种电力系统状态估计方法，其特征在于，包括以下步骤：S1、在电网中配置PMU，同时对电网进行划分，具体如下：S11、在发电机和一级、二级负荷节点处安装PMU，其中，一级、二级负荷的定义参见《民用建筑电气设计规范》；

S12、根据S11所述PMU的配置情况将电网划分为PMU可观测区域和PMU不可直接观测但SCADA可观测区域，其中，PMU可观测区域记作区域O，PMU不可直接观测但SCADA可观测区域记作区域U，NO为区域O的节点数，NU为区域U的节点数，电网中总节点数为N＝NO+NU；

S2、读取电网的数据信息，所述数据信息包括：区域O的网络参数、拓扑结构和线路阻抗，区域U的网络参数、拓扑结构和线路阻抗；

S3、根据S2所述电网中的数据信息并结合形成电力系统中相应的节点导纳矩阵和支路-节点关联矩阵；

S4、读取电网的SCADA量测配置信息,所述SCADA量测配置信息包括：节点电压幅值量测、节点电流幅值量测、节点功率注入量测和线路潮流量测；

0 0

S5、计算插值矩阵H，所述插值矩阵H 由S3所得的节点导纳矩阵经插值处理得到，具体插值处理方法如下：S51、电力系统通用节点电流方程为[Ibus]＝[Ybus]·[Ebus]，其中，Ibus是节点注入电流向量，Ybus是节点导纳矩阵，Ebus是节点电压向量；

S52、根据S12所述电网划分将节点导纳矩阵分为区域O节点自导纳矩阵YOO，区域O节点互导纳矩阵YOU，区域U的自导纳矩阵YUU，连接区域O和区域U的互导纳矩阵YUO；

S53、得到节点电流方程 其中，IO是区域O中节点注入电流，EO是区域O中的节点电压，IU是区域U中的节点注入电流，EU是区域U中的节点电压；

S54、将IU等效转化为一个负荷导纳向量YU，即，当区域U中的节点数目为NU，那么区域U中各个节点处的电流注入可以表示为 其中，Si为节点i处的复功率，Pi为节点i处的有功功率，Qi为节点i处的无功功率，[YU]＝2

[YL]＝[I/E]＝[(Pi-jQi)/|Ei|],i＝1,2,...,NU，YL为一个NU×NU的对角矩阵；

S55、将S54所得电流注入表达式[IU]代入S53所述节点电流方程，得到然后计算得到YUO·EO+YT·EU＝0，其中，YT＝YUU+YL；

0 -1 -1

S56、得到系统的插值矩阵H ＝-YT ·YUO，其中，矩阵YT ·YUO是稀疏的，所述矩阵-1YT ·YUO维数是NU×NO；

S6、根据区域O中基于PMU量测插值得到的先验信息(xk)并利用改进加权最小二乘法估计出区域U中各节点的节点状态向量 并由此计算出各节点处的注入复功率，所述区域U的节点处的注入复功率计算方法如下：S61、电力系统在第k次采样时的状态xk可以用zk＝h(xk)+vk表示，其中，x是状态向量，zk为量测向量，h(·)是m维非线性函数向量，vk是服从正态分布的随机白噪声，即vk～N(0,Rk)，Rk是量测误差的方差；

S62、考虑到PMU采样速率通常比SCADA快很多，在SCADA量测更新的间隔期间通过PMU插值得到的区域U状态可以等效为区域U的先验信息；

S63、对于区域U，以节点电压为状态变量为例，根据S62所述区域U的先验信息有如下的优化目标函数：其中, 是k时刻的先验状态，在两次SCADA量测采样间隔期间， 可由式得到，Pk是NU×NU阶的误差协方差矩阵，对上面的目标函数进行优化求解可得： 其中， Wk是对角

权重矩阵， Q是h(x)关于x的雅克比矩阵，所述权重Wk是由新的自适应量测权重分配函数来计算， 其中，λ是一个大于0的自然数；

S64、得到电力系统状态向量的迭代等式xk+1＝xk+Δx；

-6

S65、当Δx＜10 ，则算法收敛，输出当前电力系统各个节点的电压和相角值，进而估计出区域U中各节点处的注入复功率；

S7、判断电力系统是否出现较大扰动，如果 则进入S10，如果则进入S8；

0

S8、更新插值矩阵，插值矩阵的更新式子为H+ΔH，其中，ΔH为电力系统在运行点发生0

变化造成的插值矩阵H 的偏移量，所述插值矩阵更新方法如下：

0 0 -1

S81、在给定的参考运行点‘0’处有：H ＝-(YT) ·YUO，其中，所述给定的参考运行点‘0’为稳态时电力系统的初始运行状态；

* 2 -1 -1

S82、电力系统的运行点可以由H＝-(YUU+diag[(Si)/|Ei|]) ·YUO＝-YT ·YUO确定，即有YTH＝-YUO，对此等式两边求导可得ΔYT·H+YT·ΔH＝0，从而可以得-1 * 2到ΔH ＝ -YT ·ΔYT·H由 于YT＝ YUU+diag[(Si)/|Ei|]，那 么 对 YT求 导 可 得：YUU是常数矩阵,求导后为0；

S83、如果初始运行点相对应的复功率为 定义复功率偏移量为：，其中， 为当前时刻的注入复功率；那么插值矩阵的偏移量可以表示为从而插值矩阵的更新式子为H0+ΔH；

S84、将S6中已经估计出的U区域各个节点的注入复功率值和电压值带入插值矩阵的更新式子中就得到了更新后的插值矩阵；

S9、通过插值矩阵和式子

对区域U各

节点状态进行计算，得到当前时刻区域U各个节点的电压，其中，EU为区域U中的节点电压，EO为区域O中的节点电压， 和 分别为初始运行点和当前时刻的注入复功率， 和Ei分别为初始运行点和当前时刻的电压；

S10、判断收敛性，当目前电力系统量测采样次数k＜TH2时，算法收敛，则输出区域U中各个节点的电压，当k≥TH2时，则进入S5，其中，TH2为预先设定的系统量测总采样次数。

2.根据权利要求1所述的一种电力系统状态估计方法，其特征在于：将输电线等效为典型π型电路，S4所述SCADA量测采用的量测函数包括：典型π型等效电路不含非变压器支路时节点i的有功注入量测函数典型π型等效电路不含非变压器支路时节点i的无功注入量测函数典型π型等效电路不含非变压器支路时节点i到节点j的注入有功函数典型π型等效电路不含非变压器支路时节点i到节点j的注入无功潮流函数典型π型等效电路不含非变压器支路时节点i到节点j的线路电流幅值量测函数典 型π型 等 效电 路 含 变 压器 支 路 时 节 点i的 有 功注 入 量 测 函数典 型π型 等 效电 路 含 变 压器 支 路 时 节 点i的 无 功注 入 量 测 函数典型π型等效电路含变压器支路时节点i到节点j的注入有功函数典型π型等效电路含变压器支路时节点i到节点j的注入无功潮流函数典型π型等效电路含变压器支路时节点j到节点i的注入有功函数典型π型等效电路含变压器支路时节点j到节点i的注入无功潮流函数典型π型等效电路含非变压器支路时节点i到节点j的线路电流幅值量测函数其中，Vi为节点i的电压幅值，Vj为节点j的电压幅值，θi为节点i的相角，θj为节点j的相角，θij＝θi-θj为节点i和节点j的相角差，Gij+jBij为导纳矩阵的第i行第j列元素，gij+jbij为节点i到节点j间的序导纳，gsi+jbsi为节点i到节点j间的并联导纳，K为变压器非标准变比，bT为变压器标准侧的电纳。

3.根据权利要求1所述的一种电力系统状态估计方法，其特征在于：S63所述λ＝1。

4.根据权利要求1所述的一种电力系统状态估计方法，其特征在于：S8所述运行点发生变化是由拓扑变化即线路断开或者是由节点负荷突变造成。