CN101924364A - 变电站-调度中心两级分布式电网的非线性状态估计方法 - Google Patents

Abstract

本发明涉及一种变电站-调度中心两级分布式电网的非线性状态估计方法，属于电力***运行和控制技术领域。首先在变电站中利用来自PMU、RTU的多源量测，采用变电站三相无阻抗非线性多源状态估计方法进行变电站状态估计，同步辨识并剔除变电站内开关开合状态坏数据和模拟量坏数据，获得各变电站的高可靠性的三相拓扑结构、高精度的三相状态估计结果及三相不平衡度并实时上传到调度中心，同时没有参与变电站状态估计的变电站选择部分原始量测量上传到调度中心，调度中心利用上传的正序或单相量对电网进行拓扑分析及非线性单相状态估计，有效监视全网三相不平衡程度和非全相运行的情况，并实现了传统状态估计到两级分布式状态估计的过渡。

Description

变电站-调度中心两级分布式电网的非线性状态估计方法

技术领域

本发明涉及一种变电站-调度中心两级分布式电网的非线性状态估计方法，属于电力***运行和控制技术领域。

背景技术

能量管理***(Energy management system，以下简称EMS)是基于计算机的现代电力***的调度自动化***，其任务是对电力***进行实时采集、监视、分析、优化和控制决策。电力***状态估计是EMS的基础和核心环节，状态估计是利用从电力***中采集的实时量测信息，排除错误信息，计算出完整、一致和可信的电力***实时变量，保证EMS控制决策的正确性。

传统状态估计在电力调度中心实施，利用远方终端单元(Remote Terminal Unit，以下简称RTU)采集遥测、遥信生数据并通过数据采集与监控***(Supervisory Control AndData Acquisition，以下简称SCADA)上传到调度中心，集中完成全网拓扑分析和状态估计。由于传送到调度中心的信息量的先天不足，通过传统状态估计模型和算法的改进，已无法从根本上解决调度中心自动化基础数据的准确性问题，拓扑错误、非线性迭代发散、大误差等导致的集中式状态估计的不可用，已成为世界范围内调度中心高级应用实用化中的一个瓶颈问题。造成上述问题的一个根本原因在于信息分布和处理的不合理。一方面，调度中心的信息过于集中。全网电网模型(含设备参数、静态拓扑、一次单线图)的建立都要在调度中心完成，工作量随着电网规模的日益增大而显著增加。同时，调度中心的维护人员不可能对电网的每个细节都非常熟悉，潜在的出错可能性非常大，这种参数错误或拓扑错误将湮没在庞大的电网模型信息中，难以排查。另一方面，调度中心的信息在局部又不够冗余。为了避免海量信息传输和存储，厂站侧通常只将调度中心最关心的部分信息上送，因此调度中心无法获得变电站内部的大量冗余的量测采集信息。由于量测冗余度不足，致使拓扑检错和坏数据辨识成为困扰调度中心多年的难点问题。此外，这种集中式的调度中心一旦受到灾难打击，很容易导致全部功能的瘫痪，难以自愈。

基于全球定位***(global positioning system，GPS)的相量测量单元(phasormeasurement unit，以下简称PMU)能够给电力***提供高精度高采样频率的相量量测值，包括各测点的复电压、复电流、功率及频率量测，从而使电力***线性状态估计成为可能。论文《Tao Yang，Hongbin Sun，Anjan Bose.Two-level PMU-based Linear State Estimator.IEEEPES Power Systems Conference & Exposition，Seattle，USA，March 16-18，2009》提出了一种基于PMU的变电站——调度中心两级分布式线性状态估计方法，即首先在各个变电站中利用相量量测单元进行变电站本地的线性状态估计，同时辨识出模拟量坏数据和开关开合状态坏数据，获得熟数据。然后，通过通信网络将熟数据传送给调度中心，实现全网的线性状态估计。这种方法有效提高了全网状态估计的可靠性，特别是提高了拓扑检错的能力。但是，这种方法只采用PMU的量测进行线性状态估计，因此仅仅适用于装有PMU的变电站，对绝大部分无PMU的变电站无法适用。即便对于装有PMU的变电站，由于没有利用RTU量测，降低了量测冗余度和量测错误辨识能力。与此同时，该方法对网络结构进行了简化，假设***运行在三相平衡条件下，在变电站中进行单相状态估计并将单相估计结果上传给调度中心，调度中心无法监控电网三相不平衡程度及非全相运行的情况。

出于成本的考虑，在未来相当长的一段时间内，电力***不可能装设足够多的PMU，以满足全***可观测性的要求，因此一种比较可行的方案是把PMU和原有的RTU构成混合测量***，利于多源的量测数据对***进行混合的非线性状态估计，既提高量测冗余度，又具有普遍适用性。同时，实际电力***运行时，通常三相不完全平衡，结构不完全对称。例如：三相参数不对称、负荷不对称、非全相运行和检修及电气化铁路都会造成电力***的三相不平衡。电力***三相不平衡时产生的负序电流和谐波对电力设备的危害性很大，如果不采取措施，会对电力***运行造成负面影响，还可能引起继电保护装置的误动作。因此，调度中心人员需了解电网的三相不平衡度，发现非全相运行的状态，做出快速正确的调度决策，并及时维修故障部分。同时，三相不平衡也是传统状态估计精度不高的一个重要原因。因此，变电站中状态估计应采用三相模型，并向调度中心上传三相不平衡度和非全相运行信息，便于监控电网三相不平衡程度及非全相运行的情况。

本申请人曾经申请过专利申请号为201010140811.3，发明名称为《一种变电站三相无阻抗非线性多源状态估计方法》，该方法发明了利用变电站内高冗余的三相多源同步实时数据，建立三相无阻抗节点-支路模型，通过三相解耦的非线性状态估计，剔除模拟量坏数据和开关状态坏数据，得到变电站内正确的开关状态和熟数据。

发明内容

本发明的目的是提出一种变电站-调度中心两级分布式电网的非线性状态估计方法，克服已有技术中存在的缺点，以实现在变电站中同时利用来自PMU、RTU等多源量测量进行变电站三相无阻抗非线性多源状态估计，同步辨识并剔除变电站内开关开合状态坏数据和模拟量坏数据，并在调度中心采集各变电站状态估计后的熟数据及没有参与变电站状态估计的老变电站的原始量测数据，对全电网进行拓扑分析及非线性单相状态估计，有效监视全网三相不平衡程度和非全相运行的情况，并实现了传统状态估计到两级分布式状态估计的过渡。

本发明提出的变电站-调度中心两级分布式电网的非线性状态估计方法，包括以下步骤：

(1)设T1为电网的变电站的状态估计周期，T2为电网的调度中心的状态估计周期，T1＜T2；采集各变电站中各电压等级下来自电网的RTU和PMU的多源量测数据，来自RTU的量测数据为：变电站内的各三相电压幅值、三相电流幅值、三相有功功率、三相无功功率以及变电站内各开关各相的开关开合状态，来自PMU的量测数据为：变电站内的各三相复电压、三相复电流、三相有功功率和三相无功功率；

(2)设电网变电站的各电压等级下各开关三相为闭合，形成由连通的无阻抗开关支路组成的开关岛，在各开关岛内分别进行无阻抗电压状态估计和无阻抗功率状态估计，并根据无阻抗功率状态估计结果分别进行模拟量坏数据的辨识和开关开合状态坏数据辨识，剔除坏数据，根据新的开关开合状态，形成变电站内各电压等级下的新开关岛，设新开关岛有M个，则

其中，K为该变电站内电压等级个数，N_k为第k个电压等级下开关岛个数，对各新开关岛进行状态估计，得到变电站内各电压等级下各开关的三相开合状态

各节点的三相复电压

各开关上三相有功功率

和无功功率

以及各节点三相注入有功功率和无功功率

其中

分别为与变电站内的变压器、线路及容抗器相连的出线节点的三相注入有功功率和无功功率，

分别为除出线节点以外的物理节点的三相注入有功功率和无功功率，上标

表示三相；

(3)以步骤(2)的变电站内各电压等级下各新开关岛分别为一条母线，则共有M条母线，母线的复电压为通过求解

得到变电站内第m个新开关岛的母线复电压其中，

为上述各节点三相复电压

的第i个值，n为第m个新开关岛内的节点数，则母线的a、b线电压相量为

为母线a相电压相量，

为母线b相电压相量；

(4)根据步骤(2)的变电站三相状态估计结果，通过求解

得到第m个新开关岛内各开关上三相复电流

通过求解

得到第m个新开关岛内各节点的三相注入复电流

其中

为与变电站内的变压器、线路及容抗器相连的出线节点的三相注入复电流，

为除出线节点以外的物理节点的三相注入复电流，出线节点a相的电流相量为

(5)根据步骤(3)中变电站内各母线三相复电压

得到变电站内各母线三相电压的正序分量

和负序分量

$\left\{\begin{array}{c}{\stackrel{\·}{U}}_1=\frac{1}{3}({\stackrel{\·}{U}}^a+a{\stackrel{\·}{U}}^b+a^2{\stackrel{\·}{U}}^c)\\ {\stackrel{\·}{U}}_2=\frac{1}{3}({\stackrel{\·}{U}}^a+a^2{\stackrel{\·}{U}}^b+a{\stackrel{\·}{U}}^c)\end{array}\right.$

其中，

计算得到变电站内各母线的三相复电压正序分量幅值|U₁|和负序分量幅值|U₂|：

$\left\{\begin{array}{c}\left|U_1\right|=\sqrt{U_{1\mathrm{real}}^2+U_{1\mathrm{imag}}^2}\\ \left|U_2\right|=\sqrt{U_{2\mathrm{real}}^2+U_{2\mathrm{imag}}^2}\end{array}\right.$

其中，U_1real和U_1imag为变电站内各母线三相复电压正序分量

的实部和虚部，U_2real和U_2imag为负序分量

的实部和虚部，通过求解得到变电站内各母线的三相电压不平衡度；

(6)由步骤(4)变电站内各出线节点的三相注入复电流

计算得到各出线节点三相电流的正序分量

和负序分量

$\left\{\begin{array}{c}{\stackrel{\·}{I}}_1=\frac{1}{3}({\stackrel{\·}{I}}_{\mathrm{LN}}^a+a{\stackrel{\·}{I}}_{\mathrm{LN}}^b+a^2{\stackrel{\·}{I}}_{\mathrm{LN}}^c)\\ {\stackrel{\·}{I}}_2=\frac{1}{3}({\stackrel{\·}{I}}_{\mathrm{LN}}^a+a^2{\stackrel{\·}{I}}_{\mathrm{LN}}^b+a{\stackrel{\·}{I}}_{\mathrm{LN}}^c)\end{array}\right.$

其中，

计算得到各出线节点三相电流的正序分量幅值

和负序分量幅值

$\left\{\begin{array}{c}\left|{\stackrel{\·}{I}}_1\right|=\sqrt{I_{1\mathrm{real}}^2+I_{1\mathrm{imag}}^2}\\ \left|{\stackrel{\·}{I}}_2\right|=\sqrt{I_{2\mathrm{real}}^2+I_{2\mathrm{imag}}^2}\end{array}\right.$

其中，I_1real和I_1imag分别为正序分量

的实部和虚部，I_2real和I_2imag分别为负序分量

的实部和虚部，通过求解

得到变电站内各出线节点三相电流不平衡度；

(7)通过如下公式：

$\left\{\begin{array}{c}{P}_{\mathrm{CB}}=P_{\mathrm{CB}}^a+P_{\mathrm{CB}}^b+P_{\mathrm{CB}}^c\\ Q_{\mathrm{CB}}=Q_{\mathrm{CB}}^a+Q_{\mathrm{CB}}^b+Q_{\mathrm{CB}}^c\\ P_{\mathrm{ND}}=P_{\mathrm{ND}}^a+P_{\mathrm{ND}}^b+P_{\mathrm{ND}}^c\\ Q_{\mathrm{ND}}=Q_{\mathrm{ND}}^a+Q_{\mathrm{ND}}^b+Q_{\mathrm{ND}}^c\end{array}\right.,$

计算变电站内各开关三相总有功功率P_CB和三相总无功功率Q_CB，以及变电站内各节点的三相注入总有功功率

和三相总无功功率

其中P_LN、Q_LN分别为变电站内各出线节点的三相总有功功率和无功功率，P_nd、Q_nd分别为除出线节点以外的物理节点的三相总有功功功率和无功功率；

(8)计算变电站内各开关的三相总开合状态z_CB，其中第i个开关的三相总开合状态计算公式如下：

$z_{\mathrm{CB},i}=\left\{\begin{array}{c}0,\mathrm{if}{z}_{\mathrm{CB},i}^a=z_{\mathrm{CB},i}^b=z_{\mathrm{CB},i}^c=0;\\ 1,\mathrm{if}{z}_{\mathrm{CB},i}^a=z_{\mathrm{CB},i}^b=z_{\mathrm{CB},i}^c=1;\\ 2,\mathrm{alarm}\end{array}\right.$

其中，0表示开关的三相都闭合，1表示开关的三相都断开，2表示开关的三相开合状态不一致，处于报警状态，并向调度中心发送非全相运行状态的报警信号；

(9)当到达电网的变电站的状态估计周期T1时，重复步骤(1)-(9)，进行周期性变电站状态估计，当到达电网的调度中心的状态估计周期T2时，进行以下步骤；

(10)电网的调度中心从各变电站采集量测量熟数据z_se，其中，当调度中心需要正序熟数据对电网进行状态估计时，

其中，

表示正序熟数据，当调度中心需要单相熟数据对电网进行状态估计时，

表示单相熟数据，同时调度中心从各变电站采集以下变电站状态估计结果：变电站的各母线三相电压不平衡度ε_U、各出线节点三相电流不平衡度ε_I以及各开关的三相总开合状态z_CB；

(11)调度中心采集电网中不参与变电站状态估计的各变电站原始量测量z_meas，其中，当调度中心需要正序量测量对电网进行状态估计时，

其中，

表示正序量测量，|U′₁|为RTU采集到的变电站内各电气母线电压幅值，|I′₁|为RTU采集到的各出线节点的电流幅值，P′、Q′分别为RTU采集到的各出线节点的三相总的有功功率和无功功率，为PMU采集到的各电气母线正序复电压，

分别为PMU采集到的各出线节点的正序复电流，P′_PMU、Q′_PMU分别为PMU采集到的各出线节点的三相总有功功率和无功功率，当调度中心需要单相熟数据对电网进行状态估计时，

其中，

表示单相量测量，|U′_ab|为RTU采集到的变电站内各电气母线a、b线电压幅值，|I′_a|为RTU采集到的变电站内各出线节点的a相电流幅值，为PMU采集到的变电站内各电气母线a、b线电压相量，

为PMU采集到的变电站内各出线节点的a相电流相量；

(12)当调度中心的第i个开关的三相总开合状态z_CBi等于2时，即为非全相运行状态时，处于报警状态；设置一个三相电流不平衡度阈值，当变电站内各出线节点的三相复电流不平衡度ε_I大于此阈值时，处于报警状态；设置一个三相电压不平衡度阈值，当变电站内各母线的三相电压不平衡度ε_U大于此阈值时，处于报警状态；利用各变电站状态估计结果，对电网进行传统的电力***拓扑分析和状态估计，传统状态估计目标函数如下：

J(x)＝[z-h(x)]^TR^-1[z-h(x)]

其中，状态量x＝(v₁，...v_m，θ₁，...，θ_m)^T表示电网各母线的电压幅值和相角，量测矢量

包括变电站状态估计得到的量测量熟数据z_se及没有进行状态估计的变电站的原始量测量z_meas；

(13)重复步骤(1)-(13)，进行变电站-调度中心两级分布式电网的非线性状态估计。

本发明提出的变电站-调度中心两级分布式电网的非线性状态估计方法，利用已有的在变电站内进行三相无阻抗非线性多源状态估计，进而将变电站内三相状态估计结果进行转化与处理，得到变电站状态估计后的熟数据、三相不平衡度和开关状态信息并上传到调度中心，对于一些不能进行变电站状态估计的老变电站也有针对性地选择一部分原始量测数据上传到调度中心，调度中心利用上传的熟数据及原始量测对全网进行状态估计，同时监视电网三相不平衡度及非全相运行的情况，并实现了传统状态估计到两级分布式状态估计的过渡。

本发明提出的变电站-调度中心两级分布式电网的非线性状态估计方法，与传统的调度中心集中式的状态估计相比，其优点是：由于站内量测冗余度高(包括PMU和RTU不同测量装置的量测量)，采集速度快，因此，变电站状态估计可靠性高、精度高，尤其是变电站内开关上三相量测信息丰富，拓扑检错能力强，可显著减少拓扑错误。另外，由于站内网络规模小，计算速度快，实时性很高。经过变电站状态估计，上传预过滤后的熟数据到调度中心，显著减少了全网的拓扑错误和模拟量坏数据，提高了调度中心状态估计的可靠性和精度，同时变电站三相状态估计也为调度中心提供了各变电站三相不平衡度和非全相运行信息，提高了调度中心可靠地监视全网三相不平衡度和非全相运行的能力。对于不能进行站内状态估计的老变电站有针对性地选择部分原始量测上传给调度中心，实现了传统状态估计到两级分布式状态估计的过渡。

附图说明

图1是本发明方法的流程框图。

图2是本发明方法的一个实施例中采用的IEEE14节点模型。

图3是本发明方法的一个实施例中采用的变电站开关支路-节点模型。

具体实施方式

本发明提出的变电站-调度中心两级分布式电网的非线性状态估计方法，其流程框图如图1所示，包括以下步骤：

(1)设T1为电网的变电站的状态估计周期，T2为电网的调度中心的状态估计周期，T1＜T2；采集各变电站中各电压等级下来自电网的RTU和PMU的多源量测数据，来自RTU的量测数据为：变电站内的各三相电压幅值、三相电流幅值、三相有功功率、三相无功功率以及变电站内各开关各相的开关开合状态，来自PMU的量测数据为：变电站内的各三相复电压、三相复电流、三相有功功率和三相无功功率；

(2)设电网变电站的各电压等级下各开关三相为闭合，形成由连通的无阻抗开关支路组成的开关岛，在各开关岛内分别进行无阻抗电压状态估计和无阻抗功率状态估计，并根据无阻抗功率状态估计结果分别进行模拟量坏数据的辨识和开关开合状态坏数据辨识，剔除坏数据，根据新的开关开合状态，形成变电站内各电压等级下的新开关岛，设新开关岛有M个，则

其中，K为该变电站内电压等级个数，N_k为第k个电压等级下开关岛个数，对各新开关岛进行状态估计，得到变电站内各电压等级下各开关的三相开合状态

各节点的三相复电压

各开关上三相有功功率

和无功功率

以及各节点三相注入有功功率

和无功功率

其中

分别为与变电站内的变压器、线路及容抗器相连的出线节点的三相注入有功功率和无功功率，

分别为除出线节点以外的物理节点的三相注入有功功率和无功功率，上标

表示三相；

(3)以步骤(2)的变电站内各电压等级下各新开关岛分别为一条母线，则共有M条母线，母线的复电压为

通过求解

得到变电站内第m个新开关岛的母线复电压其中，

为上述各节点三相复电压

的第i个值，n为第m个新开关岛内的节点数，则母线的a、b线电压相量为

为母线a相电压相量，

为母线b相电压相量；

(4)根据步骤(2)的变电站三相状态估计结果，通过求解得到第m个新开关岛内各开关上三相复电流

通过求解得到第m个新开关岛内各节点的三相注入复电流

其中

为与变电站内的变压器、线路及容抗器相连的出线节点的三相注入复电流，

为除出线节点以外的物理节点的三相注入复电流，出线节点a相的电流相量为

(5)根据步骤(3)中变电站内各母线三相复电压

得到变电站内各母线三相电压的正序分量

和负序分量

$\left\{\begin{array}{c}{\stackrel{\·}{U}}_1=\frac{1}{3}({\stackrel{\·}{U}}^a+a{\stackrel{\·}{U}}^b+a^2{\stackrel{\·}{U}}^c)\\ {\stackrel{\·}{U}}_2=\frac{1}{3}({\stackrel{\·}{U}}^a+a^2{\stackrel{\·}{U}}^b+a{\stackrel{\·}{U}}^c)\end{array}\right.$

其中，

计算得到变电站内各母线的三相复电压正序分量幅值|U₁|和负序分量幅值|U₂|：

$\left\{\begin{array}{c}\left|U_1\right|=\sqrt{U_{1\mathrm{real}}^2+U_{1\mathrm{imag}}^2}\\ \left|U_2\right|=\sqrt{U_{2\mathrm{real}}^2+U_{2\mathrm{imag}}^2}\end{array}\right.$

其中，U_1rel和U_1imag为变电站内各母线三相复电压正序分量

的实部和虚部，U_2real和U_2imag为负序分量

的实部和虚部，通过求解得到变电站内各母线的三相电压不平衡度；

(6)由步骤(4)变电站内各出线节点的三相注入复电流

计算得到各出线节点三相电流的正序分量

和负序分量

$\left\{\begin{array}{c}{\stackrel{\·}{I}}_1=\frac{1}{3}({\stackrel{\·}{I}}_{\mathrm{LN}}^a+a{\stackrel{\·}{I}}_{\mathrm{LN}}^b+a^2{\stackrel{\·}{I}}_{\mathrm{LN}}^c)\\ {\stackrel{\·}{I}}_2=\frac{1}{3}({\stackrel{\·}{I}}_{\mathrm{LN}}^a+a^2{\stackrel{\·}{I}}_{\mathrm{LN}}^b+a{\stackrel{\·}{I}}_{\mathrm{LN}}^c)\end{array}\right.$

其中，

计算得到各出线节点三相电流的正序分量幅值

和负序分量幅值

$\left\{\begin{array}{c}\left|{\stackrel{\·}{I}}_1\right|=\sqrt{I_{1\mathrm{real}}^2+I_{1\mathrm{imag}}^2}\\ \left|{\stackrel{\·}{I}}_2\right|=\sqrt{I_{2\mathrm{real}}^2+I_{2\mathrm{imag}}^2}\end{array}\right.$

其中，I_1real和I_1imag分别为正序分量

的实部和虚部，I_2real和I_2imag分别为负序分量

的实部和虚部，通过求解

得到变电站内各出线节点三相电流不平衡度；

(7)通过如下公式：

$\left\{\begin{array}{c}{P}_{\mathrm{CB}}=P_{\mathrm{CB}}^a+P_{\mathrm{CB}}^b+P_{\mathrm{CB}}^c\\ Q_{\mathrm{CB}}=Q_{\mathrm{CB}}^a+Q_{\mathrm{CB}}^b+Q_{\mathrm{CB}}^c\\ P_{\mathrm{ND}}=P_{\mathrm{ND}}^a+P_{\mathrm{ND}}^b+P_{\mathrm{ND}}^c\\ Q_{\mathrm{ND}}=Q_{\mathrm{ND}}^a+Q_{\mathrm{ND}}^b+Q_{\mathrm{ND}}^c\end{array}\right.,$

计算变电站内各开关三相总有功功率P_CB和三相总无功功率Q_CB，以及变电站内各节点的三相注入总有功功率和三相总无功功率

其中P_LN、Q_LN分别为变电站内各出线节点的三相总有功功率和无功功率，P_nd、Q_nd分别为除出线节点以外的物理节点的三相总有功功功率和无功功率；

(8)计算变电站内各开关的三相总开合状态z_CB，其中第i个开关的三相总开合状态计算公式如下：

$z_{\mathrm{CB},i}=\left\{\begin{array}{c}0,\mathrm{if}{z}_{\mathrm{CB},i}^a=z_{\mathrm{CB},i}^b=z_{\mathrm{CB},i}^c=0;\\ 1,\mathrm{if}{z}_{\mathrm{CB},i}^a=z_{\mathrm{CB},i}^b=z_{\mathrm{CB},i}^c=1;\\ 2,\mathrm{alarm}\end{array}\right.$

其中，0表示开关的三相都闭合，1表示开关的三相都断开，2表示开关的三相开合状态不一致，处于报警状态，并向调度中心发送非全相运行状态的报警信号；

(9)当到达电网的变电站的状态估计周期T1时，重复步骤(1)-(9)，进行周期性变电站状态估计，当到达电网的调度中心的状态估计周期T2时，进行以下步骤；

(10)电网的调度中心从各变电站采集量测量熟数据z_se其中，当调度中心需要正序熟数据对电网进行状态估计时，

其中，

表示正序熟数据，当调度中心需要单相熟数据对电网进行状态估计时，

表示单相熟数据，同时调度中心从各变电站采集以下变电站状态估计结果：变电站的各母线三相电压不平衡度ε_U、各出线节点三相电流不平衡度ε_I以及各开关的三相总开合状态z_CB；

(11)调度中心采集电网中不参与变电站状态估计的各变电站原始量测量z_meas，其中，当调度中心需要正序量测量对电网进行状态估计时，其中，

表示正序量测量，|U′₁|为RTU采集到的变电站内各电气母线电压幅值，|I′₁|为RTU采集到的各出线节点的电流幅值，P′、Q′分别为RTU采集到的各出线节点的三相总的有功功率和无功功率，

为PMU采集到的各电气母线正序复电压，

分别为PMU采集到的各出线节点的正序复电流，P′_PMU、Q′_PMU分别为PMU采集到的各出线节点的三相总有功功率和无功功率，当调度中心需要单相熟数据对电网进行状态估计时，

其中，

表示单相量测量，|U′_ab|为RTU采集到的变电站内各电气母线a、b线电压幅值，|I′_a|为RTU采集到的变电站内各出线节点的a相电流幅值，

为PMU采集到的变电站内各电气母线a、b线电压相量，为PMU采集到的变电站内各出线节点的a相电流相量；

(12)当调度中心的第i个开关的三相总开合状态z_CBi等于2时，即为非全相运行状态时，处于报警状态；设置一个三相电流不平衡度阈值，当变电站内各出线节点的三相复电流不平衡度ε_I大于此阈值时，处于报警状态；设置一个三相电压不平衡度阈值，当变电站内各母线的三相电压不平衡度ε_U大于此阈值时，处于报警状态；利用各变电站状态估计结果，对电网进行传统的电力***拓扑分析和状态估计，传统状态估计目标函数如下：

J(x)＝[z-h(x)]^TR^-1[z-h(x)]

其中，状态量x＝(v₁，...v_m，θ₁，...，θ_m)^T表示电网各母线的电压幅值和相角，量测矢量

包括变电站状态估计得到的量测量熟数据z_sc。及没有进行状态估计的变电站的原始量测量z_meas；

(13)重复步骤(1)-(13)，进行变电站-调度中心两级分布式电网的非线性状态估计。

以下结合附图，介绍本发明的一个实施例。

使用本发明方法时，建立的变电站状态估计的模型如下：

应用图2所示的IEEE14节点模型。在这个模型中，由变压器连接的所有母线都在同一个变电站里，图中母线5和母线6组成变电站5，母线4和母线9组成变电站4，其余各母线分别为一个变电站，对各变电站建立开关三相模型。母线5和母线6组成的变电站5开关支路-节点模型如图3所示。

在本实施例中，以潮流计算结果作为真值来对***进行状态估计，同时，为了更好地模拟真实的实时***量测，在量测真值上加入了高斯噪声，以使三相量测在数值上有所差异。这里假设变电站中各开关节点上只有一组来自数据采集与监控***的功率和电流幅值量测。实际变电站中，可能有来自不同量测装置，如PMU和RTU的几组精度不同的多源量测数据，量测冗余度将更大。

1.)实验一：只有模拟量坏数据

对于传统的状态估计算法，变电站中没有进行变电站状态估计，只是将各出线节点的三相总的有功无功量测数据直接上传给调度中心。调度中心只能建立全***的单相模型进行单相状态估计，不能对A、B、C各相坏数据进行辨识。采用两级分布式状态估计算法时，通过变电站三相状态估计，能有效剔除各相坏数据，对量测数据进行预过滤。

在图3所示的变电站5中，假设所有变电站内所有开关开合状态均闭合，且开关开合状态没有坏数据。设定线路T1、T5、T7(即图2中的支路1-5，2-5，4-5)三相总的有功功率和无功功率量测为坏数据。其中，T1上A相的有功无功量测为坏数据，B，C相没有坏数据；T5上B相有功无功量测为坏数据，A，C相没有坏数据；T7上C相有功无功量测为坏数据，A，B相没有坏数据。其余各量测均不是坏数据。经过变电站状态估计，T1、T5、T7上的有功无功量测坏数据在变电站级就已经被剔除，不会影响调度中心全网状态估计的可靠性。

各母线电压幅值和相角估计结果如表1所示。从估计值与真值之间的量测估计误差可以看出，采用传统的状态估计方法，电压幅值和相角估计值与真值偏差较大。这是因为去路T1、T5、T7上三个强相关的功率量测都被当作坏数据给剔除了，降低了调度中心进行全电网状态估计时的量测量冗余度，从而得到误差较大的估计结果。采用两级状态估计方法时，首先利用变电站中高度冗余的量测进行变电站三相状态估计，分别剔除T1的A相、T5的B相、T7的C相上的功率坏数据，再将正确三相总功率估计值送往调度中心进行状态估计，从而得到精度较高的估计结果。

表1只有模拟量坏数据时母线电压估计结果比较

表2给出了在图2中母线1、2、3、4、5、6组成的局部网络中，采用两种方法得到的各母线的注入功率和各支路功率估计结果的比较。可以看出，采用传统状态估计方法将得到量测估计误差较大的功率估计结果。

表2只有模拟量坏数据时母线注入功率和支路功率估计结果比较(p.u)

2.)实验二：变电站中有一个开关开合状态坏数据时

在母线5和母线6组成的变电站5中(见图3)，假设开关1和6正确的开合状态为断开，其它开关状态为闭合，并且开关3上功率量测为坏数据，开关开合状态量测也为坏数据，即开关3开合状态量测为断开。如果这些坏数据直接传送到调度中心，那么错误的开关状态将导致错误的拓扑分析，即开关2为一个开关岛，开关4、5为另一个开关岛，即图2中母线5将被分解为两条母线(母线1和母线15)，从而得到与真值偏差较大的状态估计结果。采用本发明方法首先在变电站进行三相状态估计后，能在变电站内同时有效地检测与辨识出模拟量坏数据和开关开合状态坏数据并剔除，得到开关3正确的开合状态应为闭合，将正确的开关开合状态和开关岛信息送往调度中心，通过调度中心对全网进行拓扑分析得到全网母线个数应该仍为14个，从而得到可靠的全网状态估计结果。

表3给出了这种情况下两种估计方法的电压估计结果比较，可以看出，当开关上同时存在模拟量坏数据和开关开合状态坏数据时，传统的状态估计方法不能同时辨识出这两种坏数据，经拓扑分析得到错误的母线个数。采用本发明方法，能在变电站级状态估计中同时辨识出这两种坏数据，并将正确的拓扑连接关系送往调度中心，通过调度中心级状态估计得到全网正确的拓扑结构和电压大小。表4给出了两种方法下的局部网络(母线1、2、3、4、5、6及其相关支路)支路潮流和节点注入功率状态估计结果的比较。可以看出，采用传统的状态估计方法，错误的拓扑结构将会导致量测估计误差较大的估计结果。

表3同时存在模拟量坏数据和开关状态坏数据时的调度中心电压状态估计结果

表4同时存在模拟量坏数据和开关状态坏数据时母线注入功率和支路功率估计结果比较(P.U.)

3.)实验三：非全相运行情况

假设图2中支路1-5，即支路T1的A相发生断线故障，处于非全相稳态运行状态，从而导致母线5上的负荷L12的三相不平衡。但T1上的A相开关状态量测却为闭合，即坏数据。在传统的估计方法中，调度中心只能采集到三相总的遥测遥信量测，不能监视电网三相不平衡情况。由于此时全网在非全相运行状态下达到了一种新的稳态，因此，与全相运行时相比，此时的三相总的功率分布发生了变化，但开关开合状态坏数据的存在使调度中心并不能发现电网处于非全相运行状态。而在两级状态估计方法中，在变电站内能够采集到三相分相量测数据，从而得到分相估计结果，当出线节点三相电流不平衡度大于电流不平衡度阈值时，向调度中心发出报警信息，方便调度人员对全网状态进行监控并对故障部分及时进行维修。表5给出了采用传统状态估计方法时局部功率估计结果和运行状态以及采用两级状态估计方法时变电站内各出线节点的功率估计结果。从表5可以看出，传统状态估计的功率估计结果并不能反映电网正处于非全相运行状态。采用两级状态估计方法时，经变电站级状态估计可知，图3中出线节点2(支路1-5)及出线节点5(负荷L12)的三相各相间的功率状态估计结果不平衡较大。变电站级出线节点的三相复电流估计结果及三相电流不平衡度如表6所示。从表6中可以看出，上述节点2及节点5三相电流不平衡度很大，即支路T1和负荷L12三相电流不平衡度很大，需向调度中心发送报警信息。

表5传统状态估计局部功率估计结果和运行状态以及两级状态估计变电站级功率估计结果(p.u.)

表6变电站级节点注入电流估计结果及三相不平衡度(p.u.)

Claims (1)

1.一种变电站-调度中心两级分布式电网的非线性状态估计方法，其特征在于该方法包括以下步骤：

(1)设T1为电网的变电站的状态估计周期，T2为电网的调度中心的状态估计周期，T1＜T2；采集各变电站中各电压等级下来自电网的远方终端单元和相量量测单元的多源量测数据，来自远方终端单元的量测数据为：变电站内的各三相电压幅值、三相电流幅值、三相有功功率、三相无功功率以及变电站内各开关各相的开关开合状态，来自相量量测单元的量测数据为：变电站内的各三相复电压、三相复电流、三相有功功率和三相无功功率；

(2)设电网变电站的各电压等级下各开关三相为闭合，形成由连通的无阻抗开关支路组成的开关岛，在各开关岛内分别进行无阻抗电压状态估计和无阻抗功率状态估计，并根据无阻抗功率状态估计结果分别进行模拟量坏数据的辨识和开关开合状态坏数据辨识，剔除坏数据，根据新的开关开合状态，形成变电站内各电压等级下的新开关岛，设新开关岛有M个，则

其中，K为该变电站内电压等级个数，N_k为第k个电压等级下开关岛个数，对各新开关岛进行状态估计，得到变电站内各电压等级下各开关的三相开合状态

各节点的三相复电压

各开关上三相有功功率

和无功功率以及各节点三相注入有功功率

和无功功率

其中

分别为与变电站内的变压器、线路及容抗器相连的出线节点的三相注入有功功率和无功功率，

分别为除出线节点以外的物理节点的三相注入有功功率和无功功率，上标

表示三相；

(3)以步骤(2)的变电站内各电压等级下各新开关岛分别为一条母线，则共有M条母线，母线的复电压为

通过求解得到变电站内第m个新开关岛的母线复电压

其中，为上述各节点三相复电压

的第i个值，n为第m个新开关岛内的节点数，则母线的a、b线电压相量为

为母线a相电压相量，为母线b相电压相量；

(4)根据步骤(2)的变电站三相状态估计结果，通过求解

得到第m个新开关岛内各开关上三相复电流

通过求解

得到第m个新开关岛内各节点的三相注入复电流

其中

为与变电站内的变压器、线路及容抗器相连的出线节点的三相注入复电流，

为除出线节点以外的物理节点的三相注入复电流，出线节点a相的电流相量为

(5)根据步骤(3)中变电站内各母线三相复电压得到变电站内各母线三相电压的正序分量

和负序分量

$\left\{\begin{array}{c}{\stackrel{\·}{U}}_1=\frac{1}{3}({\stackrel{\·}{U}}^a+a{\stackrel{\·}{U}}^b+a^2{\stackrel{\·}{U}}^c)\\ {\stackrel{\·}{U}}_2=\frac{1}{3}({\stackrel{\·}{U}}^a+a^2{\stackrel{\·}{U}}^b+a{\stackrel{\·}{U}}^c)\end{array}\right.$

其中，

计算得到变电站内各母线的三相复电压正序分量幅值|U₁|和负序分量幅值|U₂|：

$\left\{\begin{array}{c}\left|U_1\right|=\sqrt{U_{1\mathrm{real}}^2+U_{1\mathrm{imag}}^2}\\ \left|U_2\right|=\sqrt{U_{2\mathrm{real}}^2+U_{2\mathrm{imag}}^2}\end{array}\right.$

其中，U_1real和U_1imag为变电站内各母线三相复电压正序分量

的实部和虚部，U_2real和U_2imag为负序分量

的实部和虚部，通过求解

得到变电站内各母线的三相电压不平衡度；

(6)由步骤(4)变电站内各出线节点的三相注入复电流

计算得到各出线节点三相电流的正序分量和负序分量

$\left\{\begin{array}{c}{\stackrel{\·}{I}}_1=\frac{1}{3}({\stackrel{\·}{I}}_{\mathrm{LN}}^a+a{\stackrel{\·}{I}}_{\mathrm{LN}}^b+a^2{\stackrel{\·}{I}}_{\mathrm{LN}}^c)\\ {\stackrel{\·}{I}}_2=\frac{1}{3}({\stackrel{\·}{I}}_{\mathrm{LN}}^a+a^2{\stackrel{\·}{I}}_{\mathrm{LN}}^b+a{\stackrel{\·}{I}}_{\mathrm{LN}}^c)\end{array}\right.$

其中，

计算得到各出线节点三相电流的正序分量幅值

和负序分量幅值

$\left\{\begin{array}{c}\left|{\stackrel{\·}{I}}_1\right|=\sqrt{I_{1\mathrm{real}}^2+I_{1\mathrm{imag}}^2}\\ \left|{\stackrel{\·}{I}}_2\right|=\sqrt{I_{2\mathrm{real}}^2+I_{2\mathrm{imag}}^2}\end{array}\right.$

其中，I_1real和I_1imag分别为正序分量

的实部和虚部，I_2real和I_2imag分别为负序分量的实部和虚部，通过求解

得到变电站内各出线节点三相电流不平衡度；

(7)通过如下公式：

$\left\{\begin{array}{c}{P}_{\mathrm{CB}}=P_{\mathrm{CB}}^a+P_{\mathrm{CB}}^b+P_{\mathrm{CB}}^c\\ Q_{\mathrm{CB}}=Q_{\mathrm{CB}}^a+Q_{\mathrm{CB}}^b+Q_{\mathrm{CB}}^c\\ P_{\mathrm{ND}}=P_{\mathrm{ND}}^a+P_{\mathrm{ND}}^b+P_{\mathrm{ND}}^c\\ Q_{\mathrm{ND}}=Q_{\mathrm{ND}}^a+Q_{\mathrm{ND}}^b+Q_{\mathrm{ND}}^c\end{array}\right.,$

计算变电站内各开关三相总有功功率P_CB和三相总无功功率Q_CB，以及变电站内各节点的三相注入总有功功率

和三相总无功功率

其中P_LN、Q_LN分别为变电站内各出线节点的三相总有功功率和无功功率，P_nd、Q_nd分别为除出线节点以外的物理节点的三相总有功功功率和无功功率；

(8)计算变电站内各开关的三相总开合状态z_CB，其中第i个开关的三相总开合状态计算公式如下：

$z_{\mathrm{CB},i}=\left\{\begin{array}{c}0,\mathrm{if}{z}_{\mathrm{CB},i}^a=z_{\mathrm{CB},i}^b=z_{\mathrm{CB},i}^c=0;\\ 1,\mathrm{if}{z}_{\mathrm{CB},i}^a=z_{\mathrm{CB},i}^b=z_{\mathrm{CB},i}^c=1;\\ 2,\mathrm{alarm}\end{array}\right.$

其中，0表示开关的三相都闭合，1表示开关的三相都断开，2表示开关的三相开合状态不一致，处于报警状态，并向调度中心发送非全相运行状态的报警信号；

(9)当到达电网的变电站的状态估计周期T1时，重复步骤(1)-(9)，进行周期性变电站状态估计，当到达电网的调度中心的状态估计周期T2时，进行以下步骤；

(10)电网的调度中心从各变电站采集量测量熟数据z_se。，其中，当调度中心需要正序熟数据对电网进行状态估计时，

其中，

表示正序熟数据，当调度中心需要单相熟数据对电网进行状态估计时，

表示单相熟数据，同时调度中心从各变电站采集以下变电站状态估计结果：变电站的各母线三相电压不平衡度ε_U、各出线节点三相电流不平衡度ε_I以及各开关的三相总开合状态z_CB；

(11)调度中心采集电网中不参与变电站状态估计的各变电站原始量测量z_meas，其中，当调度中心需要正序量测量对电网进行状态估计时，

其中，

表示正序量测量，|U′₁|为远方终端单元采集到的变电站内各电气母线电压幅值，|I′₁|为远方终端单元采集到的各出线节点的电流幅值，P′、Q′分别为远方终端单元采集到的各出线节点的三相总的有功功率和无功功率，

为相量量测单元采集到的各电气母线正序复电压，分别为相量量测单元采集到的各出线节点的正序复电流，P′_PMU、Q′_PMU分别为相量量测单元采集到的各出线节点的三相总有功功率和无功功率，当调度中心需要单相熟数据对电网进行状态估计时，

其中，

表示单相量测量，|U′_ab|为远方终端单元采集到的变电站内各电气母线a、b线电压幅值，|I′_a|为远方终端单元采集到的变电站内各出线节点的a相电流幅值，

为相量量测单元采集到的变电站内各电气母线a、b线电压相量，

为相量量测单元采集到的变电站内各出线节点的a相电流相量；

(12)当调度中心的第i个开关的三相总开合状态z_CB1等于2时，即为非全相运行状态时，处于报警状态；设置一个三相电流不平衡度阈值，当变电站内各出线节点的三相复电流不平衡度ε_I大于此阈值时，处于报警状态；设置一个三相电压不平衡度阈值，当变电站内各母线的三相电压不平衡度ε_U大于此阈值时，处于报警状态；利用各变电站状态估计结果，对电网进行传统的电力***拓扑分析和状态估计，传统状态估计目标函数如下：

J(x)＝[z-h(x)]^TR^-1]T[z-h(x)]

其中，状态量x＝(v₁，...v_m，θ₁，...，θ_m)^T表示电网各母线的电压幅值和相角，量测矢量包括变电站状态估计得到的量测量熟数据z_se及没有进行状态估计的变电站的原始量测量z_meas；

(13)重复步骤(1)-(13)，进行变电站-调度中心两级分布式电网的非线性状态估计。

Images