CN104809339A - 基于有功平衡的scada系统电网有功坏数据识别方法 - Google Patents

Abstract

本发明公开了基于有功平衡的SCADA系统电网有功坏数据识别方法。其具体步骤为：（1）设定采集电网有功数据的时刻，获取该时刻下电网变电站的数量m；（2）获取第f个变电站的总主变容量，设定测量误差和各个有功损耗误差系数，获取所有主变各侧及出线的有功数据；（3）设置变电站的有功数据的有功坏数据度的初始值为0；（4）对有功数据进行有功平衡校验；（5）列出变电站所有有功数据的有功坏数据度，有功坏数据度为3对应的有功数据为有功坏数据；（6）判断f是否等于m，否则令f=f+1，返回步骤（2）；是则结束对电网有功坏数据的识别。本发明的方法，能有效识别SCADA系统中的有功坏数据，防止有功坏数据对电网运行分析的不良影响。

Description

基于有功平衡的SCADA系统电网有功坏数据识别方法

技术领域

本发明涉及电力系统技术领域，特别涉及SCADA系统平台的电网有功坏数据识别技术。

背景技术

SCADA(Supervisory Control And Data Acquisition)系统，即数据采集与监视控制系统发挥着至关重要的作用。SCADA系统是以计算机为基础的DCS与电力自动化监控系统；它应用领域很广，可以应用于电力、冶金、石油、化工、燃气、铁路等领域的数据采集与监视控制以及过程控制等诸多领域。

在电力系统中，SCADA系统应用最为广泛，技术发展也最为成熟。它在远动系统中占重要地位，可以对现场的运行设备进行监视和控制，以实现数据采集、设备控制、测量、参数调节以及各类信号报警等各项功能。

对电网的历史运行情况进行分析，或对未来的负荷进行预测，都需要从SCADA系统中导出历史数据。然而，在数据采集的过程中，测量装置异常或数据传输错误都可能造成坏数据。由于电网设备多、数据量庞大，历史数据中往往存在大量的坏数据，这些坏数据将影响电网的运行分析和规划，使电网的运行分析和规划偏离实际情况，严重时会导致错误的结论。

SCADA系统中采集到的变电站有功数据，往往是调度运行、规划研究中最关心的。因此本文着重研究如何识别有功坏数据，从有功平衡入手，结合有功平衡校验提出了一套切实可行的方法。

发明内容

本发明的目的在于克服现有技术的缺点与不足，提供一种电网SCADA系统有功坏数据识别方法。

本发明的目的通过以下的技术方案实现：

基于有功平衡校验的SCADA系统电网有功坏数据识别方法，其特征在于，包含以下顺序的步骤：

S1.设定采集电网有功数据的历史时刻T，获取T时刻下电网变电站的数量m，并对变电站从1至m进行编号；定义f＝1～m，且f的初始值为1；

S2.获取第f个变电站的总主变容量S_N，设定测量误差系数ε₀、主变有功损耗误差系数ε₁、高压侧出线有功损耗误差系数ε_g、中压侧出线有功损耗误差系数ε_z和低压侧出线有功损耗误差系数ε_d；获取第f个变电站的主变台数n，令l＝1～n，l表示主变的序号；获取第f个变电站所有主变高、中和低压侧的有功数据P_gl、P_zl和P_dl；获取第f个变电站高、中和低压侧的出线数a、b和c，令i＝1～a，i表示高压侧出线的序号，令j＝1～b，j表示中压侧出线的序号，令k＝1～c，k表示低压侧出线的序号；获取第f个变电站高、中和低压侧出线的有功数据P_gxi、P_zxj和P_dxk，以及连接这些出线另外的变电站同一出线的有功数据P′_gxi、P′_zxj和P′_dxk；以上获取的有功数据均为T时刻下电网的运行数据；

S3.设置第f个变电站的有功数据的有功坏数据度为0；所述有功坏数据度包括：主变高、中和低压侧的有功坏数据度H_gl、H_zl和H_dl，主变高、中和低压侧出线的有功坏数据度H_gxi、H_zxj和H_dxk；

S4.对S2中的有功数据进行有功平衡校验；

S5.列出第f个变电站高、中和低压侧及其出线的所有有功坏数据度H_gl、H_zl、H_dl、H_gxi、H_zxj和H_dxk，有功坏数据度为3对应的有功数据，均为有功坏数据；

S6.判断f是否等于m，若f不等于m，令f＝f+1，返回步骤S2；若f等于m，结束对电网有功坏数据的识别。

所述的有功坏数据，是指从SCADA系统中采集到的变电站主变高、中和低压侧及其各条出线的有功功率数据中，由于测量装置损坏、通信设备中断或存储出错等原因，与真实值的误差超出了电网运行分析和规划所允许的误差范围的异常运行数据；有功坏数据会导致电网运行分析和规划偏离实际情况，甚至导致错误的结论，因此需要进行识别和剔除。

所述的有功坏数据度，是用来判断变电站主变和出线的有功数据是否为有功坏数据的条件，当某个有功坏数据度为3时，则判断其对应的有功数据是有功坏数据。

所述的有功平衡校验，是用来计算有功坏数据度的方法，有功平衡校验具体如下：

A.对l从1到n，判断式(1)是否成立：

|P_gl+P_zl+P_dl-ε₁S_N|<ε₀S_N    (1)

若不成立，则将H_gl、H_zl和H_dl分别加1；若成立，则H_gl、H_zl和H_dl保持原值；

B.对i从1到a，判断式(2)是否成立：

|P_gxi+P′_gxi-ε_gS_N|<ε₀S_N    (2)

若不成立，则将H_gxi加1；若成立，则H_gxi保持原值；

C.对j从1到b，判断式(3)是否成立：

|P_zxj+P′_zxj-ε_zS_N|<ε₀S_N    (3)

若不成立，则将H_zxj加1；若成立，则H_zxj保持原值；

D.对k从1到c，判断式(4)是否成立：

|P_dxk+P′_dxk-ε_dS_N|<ε₀S_N    (4)

若不成立，则将H_dxk加1；若成立，则H_dxk保持原值；

E.判断式(5)是否成立：

       $|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gl}}+\underset{i=1}{\overset{a}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gxi}}|<{\mathrm{\&epsiv;}}_0{S}_N---\left(5\right)$

若不成立，则对l从1到n，将H_gl加1，对i从1到a，将H_gxi加1；若成立，则H_gl和H_gxi保持原值；

F.判断式(6)是否成立：

       $|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zl}}+\underset{j=1}{\overset{b}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zxj}}|<{\mathrm{\&epsiv;}}_0{S}_N---\left(6\right)$

若不成立，则对l从1到n，将H_zl加1，对j从1到b，将H_zxj加1；若成立，则H_zl和H_zxj保持原值；

G.判断式(7)是否成立：

       $|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dl}}+\underset{k=1}{\overset{c}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dxk}}|<{\mathrm{\&epsiv;}}_0{S}_N---\left(7\right)$

若不成立，则对l从1到n，将H_dl加1，对k从1到c，将H_dxk加1；若成立，则H_dl和H_dxk保持原值；

H.判断式(8)是否成立：

       $|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gl}}-\underset{i=1}{\overset{a}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gxi}}^{\&prime;}+a{\mathrm{\&epsiv;}}_g{S}_N|<{\mathrm{\&epsiv;}}_0{S}_N---\left(8\right)$

若不成立，则对l从1到n，将H_gl加1；若成立，则H_gl保持原值；

I.判断式(9)是否成立：

       $|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zl}}-\underset{j=1}{\overset{b}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zxj}}^{\&prime;}+b{\mathrm{\&epsiv;}}_z{S}_N|<{\mathrm{\&epsiv;}}_0{S}_N---\left(9\right)$

若不成立，则对l从1到n，将H_zl加1；若成立，则H_zl保持原值；

J.判断式(10)是否成立：

       $|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dl}}-\underset{k=1}{\overset{c}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dxk}}^{\&prime;}+c{\mathrm{\&epsiv;}}_d{S}_N|<{\mathrm{\&epsiv;}}_0{S}_N---\left(10\right)$

若不成立，则对l从1到n，将H_dl加1；若成立，则H_dl保持原值；

K.判断式(11)是否成立：

       $|\underset{i=1}{\overset{a}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gxi}}+\underset{j=1}{\overset{b}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zxj}}+\underset{k=1}{\overset{c}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dxk}}-{\mathrm{\&epsiv;}}_1{S}_N|<{\mathrm{\&epsiv;}}_0{S}_N---\left(11\right).$

若不成立，则对i从1到a，将H_gxi加1，对j从1到b，将H_zxj加1，对k从1到c，将H_dxk加1；若成立，则H_gxi、H_zxj和H_dxk保持原值。

在进行上述的有功平衡校验时，变电站高、中和低压侧及其出线的所有有功数据均分别进行了3次校验，当某个有功数据都不满足这3次校验时，才可判断其为有功坏数据。因此，有功坏数据度为3对应的有功数据，才能被判断为有功坏数据。

本发明与现有技术相比，具有如下优点和有益效果：

(1)本发明给出了具体的识别SCADA系统电网有功坏数据的方法，为剔除历史有功坏数据，提供了方便、有效的识别方法。

(2)本发明所提供的识别方法，基于有功平衡和误差理论，并提出了评价指标有功坏数据度，通过对主变高、中和低压侧以及出线的有功坏数据度进行了三次计算，最终识别有功坏数据，结果更可靠。

附图说明

图1为本发明所述的基于有功平衡校验的SCADA系统电网有功坏数据识别方法的流程图；

图2为图1所述方法的计算实施例变电站示意图。

具体实施方式

下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。

某220kV变电站由2台变压器组成，高压侧有2条出线，中和低压侧分别有3条出线，其示意图如图2所示。

结合图1规划流程，基于有功平衡校验的电网SCADA系统有功坏数据识别方法，包括以下步骤：

(1)设定采集电网有功数据的历史时刻T，获取T时刻下电网变电站的数量m，并对变电站从1至m进行编号；定义f＝1～m，且f的初始值为1；

本实施例中，采集电网有功数据的历史时刻为2013年2月11日5:00，变电站的个数为1。

(2)获取第f个变电站的总主变容量S_N，设定测量误差系数ε₀、主变有功损耗误差系数ε₁、高压侧出线有功损耗误差系数ε_g、中压侧出线有功损耗误差系数ε_z和低压侧出线有功损耗误差系数ε_d；获取第f个变电站的主变台数n，令l＝1～n，l表示主变的序号；获取第f个变电站所有主变高、中和低压侧的有功数据P_gl、P_zl和P_dl；获取第f个变电站高、中和低压侧的出线数a、b和c，令i＝1～a，i表示高压侧出线的序号，令j＝1～b，j表示中压侧出线的序号，令k＝1～c，k表示低压侧出线的序号；获取第f个变电站高、中和低压侧出线的有功数据P_gxi、P_zxj和P_dxk，以及连接这些出线另外的变电站同一出线的有功数据P′_gxi、P′_zxj和P′_dxk；

本实施例中，S_N＝480MVA，设定ε₀＝0.01，ε₁＝0.002，ε_g＝0.002，ε_z＝0.001，ε_d＝0.001；n＝2，所有主变高、中和低压侧的有功数据如表1所示；a＝2，b＝3，c＝3，主变高、中和低压侧的出线以及连接这些出线另外的变电站同一出线的有功数据如表2所示。

表1主变的有功数据(单位：MW)

      

表2主变出线及对端出线的有功数据(单位：MW)

      

(3)设置第f个变电站的有功坏数据度为0；所述有功坏数据度包括：主变高、中和低压侧的有功坏数据度H_gl、H_zl和H_dl，主变高、中和低压侧出线的有功坏数据度H_gxi、H_zxj和H_dxk。

(4)对(2)中的有功数据进行有功平衡校验；

A.对l从1到n，判断式(1)是否成立：

|P_gl+P_zl+P_dl-ε₁S_N|<ε₀S_N    (1)若不成立，则将H_gl、H_zl和H_dl分别加1；若成立，则H_gl、H_zl和H_dl保持原值；本实施例中，

ε₀S_N＝0.01×480＝4.8MW

|P_g1+P_z1+P_d1-ε₁S_N|＝|60.4-54.1-8.0-0.002×480|＝2.66MW<ε₀S_N

|P_g2+P_z2+P_d2-ε₁S_N|＝|32.3-58.6-4.2-0.002×480|＝31.46MW>ε₀S_N

因此，H_g1＝H_zl＝H_dl＝0，H_g2＝H_z2＝H_d2＝1；

B.对i从1到a，判断式(2)是否成立：

|P_gxi+P′_gxi-ε_gS_N|<ε₀S_N    (2)若不成立，则将H_gxi加1；若成立，则H_gxi保持原值；

本实施例中，

|P_gx1+P′_gx1-ε_gS_N|＝|173.5-176.6-0.002×480|＝4.06MW<ε₀S_N

|P_gx2+P′_gx2-ε_gS_N|＝|-304.6+305.7-0.002×480|＝0.14MW<ε₀S_N

因此，H_gx1＝H_gx2＝0；

C.对j从1到b，判断式(3)是否成立：

|P_zxj+P′_zxj-ε_zS_N|<ε₀S_N    (3)若不成立，则将H_zxj加1；若成立，则H_zxj保持原值；

本实施例中，

|P_zx1+P′_zx1-ε_zS_N|＝|32.5-31.8-0.001×480|＝0.22MW<ε₀S_N

|P_zx2+P′_zx2-ε_zS_N|＝|60.8-58.4-0.001×480|＝1.92MW<ε₀S_N

|P_zx3+P′_zx3-ε_zS_N|＝|-20.2-18.6-0.001×480|＝39.28MW>ε₀S_N

因此，H_zx1＝H_zx2＝0，H_zx3＝1；

D.对k从1到c，判断式(4)是否成立：

|P_dxk+P′_dxk-ε_dS_N|<ε₀S_N    (4)若不成立，则将H_dxk加1；若成立，则H_dxk保持原值；

本实施例中，

|P_dx1+P′_dx1-ε_dS_N|＝|4.6-4.5-0.001×480|＝0.38MW<ε₀S_N

|P_dx2+P′_dx2-ε_dS_N|＝|3.8-3.9-0.001×480|＝0.58MW<ε₀S_N

|P_dx3+P′_dx3-ε_dS_N|＝|3.6-3.2-0.001×480|＝0.08MW<ε₀S_N

因此，H_dx1＝H_dx2＝H_dx3＝0；

E.判断式(5)是否成立：

       $|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gl}}+\underset{i=1}{\overset{a}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gxi}}|<{\mathrm{\&epsiv;}}_0{S}_N---\left(5\right)$

若不成立，则对l从1到n，将H_gl加1，对i从1到a，将H_gxi加1；若成立，则H_gl和H_gxi保持原值；

本实施例中，

       $|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gl}}+\underset{i=1}{\overset{a}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gxi}}|=|60.4+32.3+173.5-304.6|=38.6\mathrm{MW}>{\mathrm{\&epsiv;}}_0{S}_N$

因此，H_g1＝H_gx1＝H_gx2＝1，H_g2＝2；

F.判断式(6)是否成立：

       $|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zl}}+\underset{j=1}{\overset{b}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zxj}}|<{\mathrm{\&epsiv;}}_0{S}_N---\left(6\right)$

若不成立，则对l从1到n，将H_zl加1，对j从1到b，将H_zxj加1；若成立，则H_zl和H_zxj保持原值；

本实施例中，

       $|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zl}}+\underset{j=1}{\overset{b}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zxj}}|=|-54.1-58.6+32.5+60.8-20.2|=39.6\mathrm{MW}>{\mathrm{\&epsiv;}}_0{S}_N$

因此，H_z1＝H_zx1＝H_zx2＝1，H_z2＝H_zx3＝2；

G.判断式(7)是否成立：

       $|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dl}}+\underset{k=1}{\overset{c}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dxk}}|<{\mathrm{\&epsiv;}}_0{S}_N---\left(7\right)$

若不成立，则对l从1到n，将H_dl加1，对k从1到c，将H_dxk加1；若成立，则H_dl和H_dxk保持原值；

本实施例中，

       $|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dl}}+\underset{k=1}{\overset{c}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dxk}}|=|-8.0-4.2+4.6+3.8+3.6|=0.2\mathrm{MW}>{\mathrm{\&epsiv;}}_0{S}_N$

因此，H_d1＝H_dx1＝H_dx2＝H_dx3＝0，H_d2＝1；

H.判断式(8)是否成立：

       $|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gl}}-\underset{i=1}{\overset{a}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gxi}}^{\&prime;}+a{\mathrm{\&epsiv;}}_g{S}_N|<{\mathrm{\&epsiv;}}_0{S}_N---\left(8\right)$

若不成立，则对l从1到n，将H_gl加1；若成立，则H_gl保持原值；

本实施例中，

       $\begin{array}{c}|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gl}}-\underset{i=1}{\overset{a}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gxi}}^{\&prime;}+a{\mathrm{\&epsiv;}}_g{S}_N|=|60.4+32.3-(-176.6+305.7)+2\&times;0.002\&times;480|\\ =34.48\mathrm{MW}>{\mathrm{\&epsiv;}}_0{S}_N\end{array}$

因此，H_g1＝2，H_g2＝3；

I.判断式(9)是否成立：

       $|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zl}}-\underset{j=1}{\overset{b}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zxj}}^{\&prime;}+b{\mathrm{\&epsiv;}}_z{S}_N|<{\mathrm{\&epsiv;}}_0{S}_N---\left(9\right)$

若不成立，则对l从1到n，将H_zl加1；若成立，则H_zl保持原值；

本实施例中，

       $\begin{array}{c}|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zl}}-\underset{j=1}{\overset{b}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zxj}}^{\&prime;}+b{\mathrm{\&epsiv;}}_z{S}_N|=|-54.1-58.6-(-31.8-58.4-18.6)+3\&times;0.001\&times;480|\\ =2.46\mathrm{MW}<{\mathrm{\&epsiv;}}_0{S}_N\end{array}$

因此，H_z1＝1，H_z2＝2；

J.判断式(10)是否成立：

       $|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dl}}-\underset{k=1}{\overset{c}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dxk}}^{\&prime;}+c{\mathrm{\&epsiv;}}_d{S}_N|<{\mathrm{\&epsiv;}}_0{S}_N---\left(10\right)$

若不成立，则对l从1到n，将H_dl加1；若成立，则H_dl保持原值；

本实施例中，

       $\begin{array}{c}|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dl}}-\underset{k=1}{\overset{c}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dxk}}^{\&prime;}+c{\mathrm{\&epsiv;}}_d{S}_N|=|-8.0-4.2-(-4.5-3.9-3.2)+3\&times;0.001\&times;480|\\ =0.84\mathrm{MW}<{\mathrm{\&epsiv;}}_0{S}_N\end{array}$

因此，H_d1＝0，H_d2＝1；

K.判断式(11)是否成立：

       $|\underset{i=1}{\overset{a}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gxi}}+\underset{j=1}{\overset{b}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zxj}}+\underset{k=1}{\overset{c}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dxk}}-{\mathrm{\&epsiv;}}_1{S}_N|<{\mathrm{\&epsiv;}}_0{S}_N---\left(11\right).$

若不成立，则对i从1到a，将H_gxi加1，对j从1到b，将H_zxj加1，对k从1到c，将H_dxk加1；若成立，则H_gxi、H_zxj和H_dxk保持原值。

本实施例中，

       $\begin{array}{c}|\underset{i=1}{\overset{a}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gxi}}+\underset{j=1}{\overset{b}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zxj}}+\underset{k=1}{\overset{c}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dxk}}-{\mathrm{\&epsiv;}}_1{S}_N|=\left|\begin{array}{c}173.5-304.6+32.5+60.8-20.2+4.6+3.8+3.6\\ -0.002\&times;480\end{array}\right|\\ =46.96\mathrm{MW}>{\mathrm{\&epsiv;}}_0{S}_N\end{array}$

因此，H_gx1＝H_gx2＝H_zx1＝H_zx2＝2，H_zx3＝3，H_dx1＝H_dx2＝H_dx3＝1。

(5)列出第f个变电站高、中和低压侧及其出线的所有有功坏数据度H_gl、H_zl、H_dl、H_gxi、H_zxj和H_dxk，有功坏数据度为3对应的有功数据，均为有功坏数据；

本实施例中，

H_d1＝0；H_z1＝H_d2＝H_dx1＝H_dx2＝H_dx3＝1；H_g1＝H_z2＝H_gx1＝H_gx2＝H_zx1＝H_zx2＝2；H_g2＝H_zx3＝3；其中H_g2和H_zx3的值为3，由此可见，P_g2＝32.3MW和P_zx3＝﹣20.2MW为有功坏数据。

可见，采用本发明所提出的基于有功平衡校验的SCADA系统电网有功坏数据识别方法，能有效地识别SCADA系统采集的有功坏数据，对更好地分析电网历史运行情况提供较为准确的数据支撑。

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他任何未背离本发明的精神实质和原理下所作的修改、修饰、替代、组合、简化，均应为等效的置换方式，都应包含在本发明的保护范围之内。

Claims (3)

1.基于有功平衡的SCADA系统电网有功坏数据识别方法，其特征在于，包含以下步骤：

S1.设定采集电网有功数据的历史时刻T，获取T时刻下电网变电站的数量m，并对变电站从1至m进行编号；定义f＝1～m，且f的初始值为1；

S2.获取第f个变电站的总主变容量S_N，设定测量误差系数ε₀、主变有功损耗误差系数ε₁、高压侧出线有功损耗误差系数ε_g、中压侧出线有功损耗误差系数ε_z和低压侧出线有功损耗误差系数ε_d；获取第f个变电站的主变台数n，令l＝1～n，l表示主变的序号；获取第f个变电站所有主变高、中和低压侧的有功数据P_gl、P_zl和P_dl；获取第f个变电站高、中和低压侧的出线数a、b和c，令i＝1～a，i表示高压侧出线的序号，令j＝1～b，j表示中压侧出线的序号，令k＝1～c，k表示低压侧出线的序号；获取第f个变电站高、中和低压侧出线的有功数据P_gxi、P_zxj和P_dxk，以及连接这些出线另外的变电站同一出线的有功数据P′_gxi、P′_zxj和P′_dxk；以上获取的有功数据均为T时刻下电网的运行数据；

S3.设置第f个变电站的有功数据的有功坏数据度为0；所述有功坏数据度包括：主变高、中和低压侧的有功坏数据度H_gl、H_zl和H_dl，主变高、中和低压侧出线的有功坏数据度H_gxi、H_zxj和H_dxk；

S4.对步骤S2中的有功数据进行有功平衡校验；

S5.列出第f个变电站高、中和低压侧及其出线的所有有功坏数据度H_gl、H_zl、H_dl、H_gxi、H_zxj和H_dxk，有功坏数据度为3对应的有功数据，均为有功坏数据；

S6.判断f是否等于m，若f不等于m，令f＝f+1，返回步骤S2；若f等于m，结束对电网有功坏数据的识别。

2.根据权利要求1所述的基于有功平衡校验的电网SCADA系统有功坏数据识别方法，其特征在于，所述的有功坏数据度是用来判断变电站主变和出线的有功数据是否为有功坏数据的条件，当某个有功坏数据度为3时，则判断其对应的有功数据是有功坏数据。

3.根据权利要求1所述的基于有功平衡校验的SCADA系统电网有功坏数据识别方法，其特征在于，所述的有功平衡校验，是用来计算有功坏数据度的方法，有功平衡校验具体如下：

A.对l从1到n，判断式(1)是否成立：

|P_gl+P_zl+P_dl-ε₁S_N|<ε₀S_N    (1)

若不成立，则将H_gl、H_zl和H_dl分别加1；若成立，则H_gl、H_zl和H_dl保持原值；

B.对i从1到a，判断式(2)是否成立：

|P_gxi+P′_gxi-ε_gS_N|<ε₀S_N    (2)

若不成立，则将H_gxi加1；若成立，则H_gxi保持原值；

C.对j从1到b，判断式(3)是否成立：

|P_zxj+P′_zxj-ε_zS_N|<ε₀S_N    (3)

若不成立，则将H_zxj加1；若成立，则H_zxj保持原值；

D.对k从1到c，判断式(4)是否成立：

|P_dxk+P′_dxk-ε_dS_N|<ε₀S_N    (4)

若不成立，则将H_dxk加1；若成立，则H_dxk保持原值；

E.判断式(5)是否成立：

$|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gl}}+\underset{i=1}{\overset{a}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gxi}}|<{\mathrm{\&epsiv;}}_0{S}_N---\left(5\right)$

若不成立，则对l从1到n，将H_gl加1，对i从1到a，将H_gxi加1；若成立，则H_gl和H_gxi保持原值；

F.判断式(6)是否成立：

$|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zl}}+\underset{j=1}{\overset{b}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zxj}}|<{\mathrm{\&epsiv;}}_0{S}_N---\left(6\right)$

若不成立，则对l从1到n，将H_zl加1，对j从1到b，将H_zxj加1；若成立，则H_zl和H_zxj保持原值；

G.判断式(7)是否成立：

$|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dl}}+\underset{k=1}{\overset{c}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dxk}}|<{\mathrm{\&epsiv;}}_0{S}_N---\left(7\right)$

若不成立，则对l从1到n，将H_dl加1，对k从1到c，将H_dxk加1；若成立，则H_dl和H_dxk保持原值；

H.判断式(8)是否成立：

$|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gl}}-\underset{i=1}{\overset{a}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gxi}}^{\&prime;}+{\mathrm{a\&epsiv;}}_g{S}_N|<{\mathrm{\&epsiv;}}_0{S}_N---\left(8\right)$

若不成立，则对l从1到n，将H_gl加1；若成立，则H_gl保持原值；

I.判断式(9)是否成立：

$|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zl}}-\underset{j=1}{\overset{b}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zxj}}^{\&prime;}+{\mathrm{b\&epsiv;}}_z{S}_N|<{\mathrm{\&epsiv;}}_0{S}_N---\left(9\right)$

若不成立，则对l从1到n，将H_zl加1；若成立，则H_zl保持原值；

J.判断式(10)是否成立：

$|\underset{l=1}{\overset{n}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dl}}-\underset{k=1}{\overset{c}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dxk}}^{\&prime;}+{\mathrm{c\&epsiv;}}_d{S}_N|<{\mathrm{\&epsiv;}}_0{S}_N---\left(10\right)$

若不成立，则对l从1到n，将H_dl加1；若成立，则H_dl保持原值；

K.判断式(11)是否成立：

$|\underset{i=1}{\overset{a}{\mathrm{\&Sigma;}}}{P}_{\mathrm{gxi}}+\underset{j=1}{\overset{b}{\mathrm{\&Sigma;}}}{P}_{\mathrm{zxj}}+\underset{k=1}{\overset{c}{\mathrm{\&Sigma;}}}{P}_{\mathrm{dxk}}-{\mathrm{\&epsiv;}}_1{S}_N|{\mathrm{\&epsiv;}}_0{S}_N---\left(11\right)$

若不成立，则对i从1到a，将H_gxi加1，对j从1到b，将H_zxj加1，对k从1到c，将H_dxk加1；若成立，则H_gxi、H_zxj和H_dxk保持原值。