CN113595058A - 配电网状态估计量测数据的处理方法 - Google Patents

配电网状态估计量测数据的处理方法 Download PDF

Info

Publication number
CN113595058A
CN113595058A CN202110471786.5A CN202110471786A CN113595058A CN 113595058 A CN113595058 A CN 113595058A CN 202110471786 A CN202110471786 A CN 202110471786A CN 113595058 A CN113595058 A CN 113595058A
Authority
CN
China
Prior art keywords
data
measurement
node
pmu
phasor
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202110471786.5A
Other languages
English (en)
Inventor
何西
董恒
陈孟娜
刘柏伸
贺飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hunan Institute of Technology
Original Assignee
Hunan Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hunan Institute of Technology filed Critical Hunan Institute of Technology
Priority to CN202110471786.5A priority Critical patent/CN113595058A/zh
Publication of CN113595058A publication Critical patent/CN113595058A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J13/00—Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network
    • H02J13/00001—Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network characterised by the display of information or by user interaction, e.g. supervisory control and data acquisition systems [SCADA] or graphical user interfaces [GUI]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/003—Load forecast, e.g. methods or systems for forecasting future load demand
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/22—Flexible AC transmission systems [FACTS] or power factor or reactive power compensating or correcting units
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/40—Display of information, e.g. of data or controls

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Human Computer Interaction (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

配电网状态估计量测数据的处理方法，涉及电力***安全技术领域，包括微型同步相量测量装置D‑PMU的量测数据和高级量测体系AMI的量测数据的处理方法，通过对微型同步相量测量装置D‑PMU的量测数据和高级量测体系AMI的量测数据的处理，处理后的微型同步相量测量装置D‑PMU的量测数据和高级量测体系AMI的量测数据与配网SCADA***的数据构成一个混合量测***来实现配电网状态估计。通过处理D‑PMU的量测数据，解决雅可比矩阵病态构造的问题；对AMI量测的低压节点负荷进行超短期预测，得到了同一时间断面注入功率数据。本发明还提出了一种抗虚假数据攻击方法，基于上述的配电网状态估计量测数据的处理方法，通过伪D‑PMU量测数据来提高量测冗余度，达到抗虚假数据攻击的目的。

Description

配电网状态估计量测数据的处理方法
技术领域
本发明涉及电力***安全技术领域，特别涉及一种配电网状态估计量测数据的处理方法。
背景技术
智能配电网通过布置智能量测终端进行数据采集，为状态估计，故障诊断、新能源接入、孤岛控制等高级应用提供支撑，而其中状态估计是基础，只有得到整个***的电压、电流、功率数据才能进行更为复杂的分析、判断与控制。
相对于输电网状态估计，配电网状态估计中存在的最大问题就是量测不足。目前，我国配网馈线自动化率不足30％，大部分节点和线路缺乏量测。随着新型量测技术的发展，比如微型同步相量测量装置(D-PMU)，高级量测体系(AMI)等，但是如何最大化利用这些量测成为实现配网状态估计的关键。
另一方面，配网状态估计还容易受到虚假数据注入攻击(FDIA)的影响。攻击者通过向原量测值z中注入事先设定好的虚假数据a，使得输入数据变成zbad＝z+a，估计状态量将由原来的
Figure BDA0003045715660000011
变成
Figure BDA0003045715660000012
文献证明，只要满足条件a＝Hc，目标函数
Figure BDA0003045715660000013
将等于
Figure BDA0003045715660000014
也就是说，虽然量测值被篡改，但由于注入的虚假数据满足特定的条件，使得坏数据检测的目标函数没变，因此无法将虚假数据检测出来。
现有技术仅在配电网“关键节点”配置D-PMU,且状态估计算法通常是采用极坐标的形式对量测函数和雅可比矩阵进行操作，D-PMU的电压和电流量测的原始值也是极坐标形式。然而，对于线路较短的配电网来说，并联导纳通常可以忽略，从而导致电流相量在求偏微分时分母为零，出现方程病态构造。而传统的配网状态估计大多利用历史负荷数据对节点注入功率进行预测，但这种方法准确率不高。
发明内容
针对上述现有技术所存在的缺陷，本发明要解决的技术问题是提出了一种配电网状态估计量测数据的处理方法，该处理方法可解决雅可比矩阵病态构造的问题以及无法获取同一时间断面节点注入功率的问题。
为解决上述技术问题，本发明采用的技术方案如下：配电网状态估计量测数据的处理方法，包括微型同步相量测量装置D-PMU的量测数据和高级量测体系AMI的量测数据的处理方法，通过对所述微型同步相量测量装置D-PMU的量测数据和高级量测体系AMI的量测数据进行处理，处理后的微型同步相量测量装置D-PMU的量测数据和高级量测体系AMI的量测数据与配网SCADA***的数据构成一个混合量测***来实现配电网状态估计。
处理微型同步相量测量装置D-PMU的量测数据的步骤是：先确定微型同步相量测量装置D-PMU的安装节点，采集安装节点的电压相量数据和线路的电流相量数据；将采集到的线路的电流相量数据通过直角坐标的形式进行处理来获取雅可比矩阵中与电流相量的相关项；采集到的安装节点的电压相量数据用于计算伪D-PMU量测，从而提高量测冗余度。
进一步来说，处理高级量测体系AMI的量测数据的步骤是：先确定操作点的数据与实际数据的最大时延；再采用超短期预测方法对操作点实际数据进行预测，从而得到同一时间断面注入功率数据。
进一步来说，安装节点的电压相量数据为V＝V∠δ，其中，V表示安装节点的电压，δ表示安装节点的电压相量的角度。
进一步来说，线路的电流相量数据为：
I＝I∠θ
或I＝(V安装节点-V安装节点的相邻节点)·Y (1)
上式(1)中，I表示安装节点与安装节点的相邻节点之间线路的电流，θ表示线路电流的角度，Y为导纳矩阵。
进一步来说，通过直角坐标形式处理的电流相量数据表示为：
I12.a＝Y11(V1.a-V2.a)+Y12(V1.b-V2.b)+Y13(V1.c-V2.c) (2)
上式(2)中，I12.a表示安装节点与安装节点的相邻节点之间A相的线路电流，V1.a、V1.b、V1.c表示安装节点的A、B、C相的电压，V2.a、V2.b、V2.c表示安装节点的相邻节点的A、B、C相的电压。
进一步来说，通过安装节点的电压相量数据计算与安装节点相邻的节点的电压相量：
Figure BDA0003045715660000031
与安装节点相邻的节点的电压相量V2即伪D-PMU量测，上式(3)中I12.a、I12.b、I12.c表示安装节点与安装节点的相邻节点之间线路A、B、C相的电流，z表示量测值，V1.a、V1.b、V1.c表示安装节点的A、B、C相的电压；
以A相为例，伪D-PMU量测在直角坐标系中表示为：
Figure BDA0003045715660000041
上式(4)中，上标“r”表示相量的实部，“i”表示虚部；
基于不确定性传播原理分析伪D-PMU量测的标准差，所述伪D-PMU量测的标准差表示为：
Figure BDA0003045715660000042
Figure BDA0003045715660000043
上式(5)、(6)中，x＝[V1.a,δ1.a,I12.a,θ12.a,I12.b,θ12.b,I12.c,θ12.c]，σ(x(k))表示微型同步相量测量装置D-PMU各个量测数据的标准差。
最后通过计算所述伪D-PMU量测的标准差确定计算半径。
进一步来说，超短期预测方法是将单个用户负荷数据看成正弦信号，正弦信号采用傅里叶级数形式的谐波分量模型表示。
进一步来说，超短期预测表示为：
Figure BDA0003045715660000044
上式(7)中，l′(last)表示最近一次采样的高级量测体系AMI的量测数据，参数an和bn通过历史数据计算得到，T表示操作点与最近一次数据更新的时间间隔，数据更新频率为96次/天，N表示谐波分量阶数。
本发明所取得有益效果在于：通过处理微型同步相量量测装置D-PMU的量测数据，解决雅可比矩阵病态构造的问题；基于谐波分量分解方法对高级量测体系AMI量测的低压节点负荷进行超短期预测，得到了同一时间断面注入功率数据。通过对D-PMU临近节点或线路量测的电压、电流和功率进行计算，得到“伪D-PMU量测数据”，基于不确定性传播原理确定计算半径，提高了量测冗余度，从而达到抗虚假数据攻击的目的。
附图说明
图1为本实施例中考虑互感的三相两节点网络；
图2为本实施例中的低压配电节点示意图；
图3为本实施例中的AMI量测时序图；
图4为本实施例中的攻击区域划分的模型；
图5为本实施例中的IEEE 123节点测试***的网络拓扑图；
图6为本实施例中的四种方法对节点55的有功功率和无功功率超短期预测；
图7为本实施例中存在虚假数据攻击的情况下的伪D-PMU量测对于状态估计精度的影响；
图8为本实施例中存在虚假数据攻击的情况下伪D-PMU量测对于状态估计一致性的影响。
具体实施方式
为了便于本领域技术人员的理解，下面结合实施例与附图对本发明作进一步的说明，实施方式提及的内容并非对本发明的限定。
配电网状态估计量测数据的处理方法，包括微型同步相量测量装置D-PMU的量测数据和高级量测体系AMI的量测数据的处理方法，通过对微型同步相量测量装置D-PMU的量测数据和高级量测体系AMI的量测数据的处理，处理后的微型同步相量测量装置D-PMU的量测数据和高级量测体系AMI的量测数据与配网SCADA***的数据构成一个完整的混合量测***，实现配电网状态估计，如图1所示。
实施例1
微型同步相量测量装置D-PMU量测数据的处理；现有技术中出于经济性考虑，仅在配电网“关键节点”配置D-PMU。如果采用***各节点的电压幅值及相角作为状态量，那么对于微型同步相量测量装置D-PMU测量的电压相量数据不需要做任何处理，只需对电流相量数据进行处理。
一般的状态估计算法采用极坐标的形式对量测函数和雅克比矩阵进行操作，微型同步相量测量装置D-PMU的电压和电流量测的原始值也是极坐标形式。然而，对于线路较短的配电网来说，并联导纳通常可以忽略，从而导致电流相量在求偏微分时分母为零，出现方程病态构造。
处理微型同步相量测量装置D-PMU的量测数据的具体步骤如下：先确定微型同步相量测量装置D-PMU的安装节点，采集安装节点的电压相量数据和线路的电流相量数据；将采集到的线路的电流相量数据通过直角坐标的形式进行处理来获取雅可比矩阵中与电流相量的相关项；采集到的安装节点的电压相量数据用于计算伪D-PMU量测，从而提高量测冗余度。
具体以如图1所示的考虑互感的三相两节点电路为例，假设微型同步相量测量装置D-PMU安装于节点1，因此可以得到节点1的电压相量数据和线路1-2的电流相量数据，安装节点的电压相量为V1＝V1∠δ1,线路电流相量为I12＝I12∠θ12或I12＝(V1-V2)·Y，其中Y为导纳矩阵。
接着采用直角坐标的形式对微型同步相量测量装置D-PMU采集的电流相量进行处理，以A相为例进行说明，电流相量I12可表示为
I12a＝Y11(V1a-V2a)+Y12(V1b-V2b)+Y13(V1c-V2c) (1-1)
因此，获取的雅可比矩阵中与电流相量的相关项可以表示为：
Figure BDA0003045715660000071
Figure BDA0003045715660000072
Figure BDA0003045715660000073
Figure BDA0003045715660000074
上式(1-2)至(1-5)中，
Figure BDA0003045715660000075
和
Figure BDA0003045715660000076
分别表示线路电流相量I12的实部和虚部。
综上所述，当电流相量采用直角坐标表示时，在求解雅可比矩阵的过程中，偏微分方程更加简洁，且不会出现矩阵病态构造的问题。
实施例2
高级量测体系AMI量测数据的处理；随着电能计量表等值智能终端在用户侧的普及，可以得到更为准确的有功和无功注入功率数据。用户侧电能计量表的数据采集频率约为15分钟至30分钟一次。由于这些装置自身的采样时延和网络传送时延，从AMI数据采样到控制中心接收到数据，中间存在较大的延时。本文中考虑较为极端的情况，假设存在2分钟的延时。
在本实施例中由于低压节点距离用户较近，假设线路的功率损耗为零，因此在低压节点的注入功率等于各个高级量测体系AMI的量测数据之和。从图2可以看出，每个低压节点(n)由其所辖的多个终端节点构成(n1，n2，n3)。在一个典型的城市配电网中，每个终端节点下面有100-150个用户。低压节点i在时刻t的有功和无功注入功率可以分别表示成：
Figure BDA0003045715660000081
Figure BDA0003045715660000082
不同于微型同步相量测量装置D-PMU的量测数据，高级量测体系AMI数据不带时标信息，且由于数据传输带宽限制，各个计量表计不能接收同一时刻的传输数据。如图3所示，各个高级量测体系AMI的数据采样与数据接收间隔都是15分钟，采用与接收时延为2分钟。假设一个低压节点含有3个高级量测体系AMI，M1，M2和M3。如果想知道在时间2(分钟为计量单位)此节点的功率数据，对于M1来说，时间2数据恰好更新，可以直接使用更新数据，但对于M2来说，最近一次更新位于时间-12，而时间-12对应的更新时间为时间-14。同样对于M3来说，距离时间2最近的一次数据更新是在时间-9，其对应的采样时间是时间-11。可以看到，要想得到此低压节点在时间2的功率数据，需要用到数分钟之前的采样数据。从M2可以看到，最多可能需要用到16分钟之前的数据。
由于操作点的数据与实际数据的最大时延为16分钟，采用超短期预测的方法对操作点实际数据进行预测。超短期预测方法是将单个用户负荷数据看成正弦信号，采用傅里叶级数形式的谐波分量模型对其表示，超短期预测表示为：
Figure BDA0003045715660000091
式(2-3)中，l′(last)表示最近一次采样的AMI数据，T表示操作点与最近一次数据更新的时间间隔，数据更新的频率为96次/天，N表示谐波分量阶数，在本实施例中对不同的高级量测体系AMI的数据采用不同的谐波分量阶数。
同样以图3为例来说明基于谐波分量分解的超短期预测方法。对于M1，M2和M3来说，操作点即所希望知道的时间点数据，用OP表示与最近采样点即与操作点最接近的上一次采样时间点，用LS表示的时间间隔分别是2，16和13。由于谐波分量阶数的最大值为5且为现有技术中已知的，因此通过不同的LS与OP间隔设置不同的谐波分量阶数，如表1所示。
表1不同间隔对应的谐波分量阶数
Figure BDA0003045715660000101
因此对于高级量测体系AMI数据M1，M2和M3，在时间2的超短期预测可以表示为：
M1：l(2)＝l′(2) (2-4)
Figure BDA0003045715660000102
Figure BDA0003045715660000103
参数an和bn可以通过历史数据计算得到。进一步的，由于这些超短期预测数据呈正态分布，其误差跟LS与OP的时间间隔相关，因此我们可以得到预测数据的标准差为：
Figure BDA0003045715660000104
式(2-7)中，σ2(last)和σ2(last-1)表示最近两次更新的方差，INlength表示间隔时间。
为了充分利用微型同步相量测量装置D-PMU数据，对与微型同步相量测量装置D-PMU安装节点相邻的安装节点的电压进行计算，并分析其误差。在如图1所示中，假设节点1为微型同步相量测量装置D-PMU的安装节点，节点2的电压相量称为“伪D-PMU量测”，计算节点2的电压相量V2：
Figure BDA0003045715660000111
式(2-8)中，I12.a、I12.b、I12.c表示安装节点与安装节点的相邻节点之间A、B、C相的线路电流，z表示量测值，V1.a、V1.b、V1.c表示安装节点的A、B、C相的电压。
然后以A相为例，基于经典的不确定性传播原理对其标准差进行分析，V2.a在直角坐标系中可表示为：
Figure BDA0003045715660000112
式中，
Figure BDA0003045715660000113
表示相量的实部，
Figure BDA0003045715660000114
表示虚部；对实部和虚部分别展开：
Figure BDA0003045715660000115
Figure BDA0003045715660000116
式(2-10)、(2-11)中，δ和θ分别表示电压和电流的角度，r和x分别表示电阻和电抗。从V2.a的实部和虚部展开式可以看出，节点2电压相量的实部和虚部由安装于节点1的微型同步相量测量装置D-PMU的各个量测数据决定。
Figure BDA0003045715660000117
与
Figure BDA0003045715660000118
可以表示为：
Figure BDA0003045715660000119
Figure BDA00030457156600001110
根据不确定性传播原理，节点2的电压相量V2的标准差表示为：
Figure BDA0003045715660000121
Figure BDA0003045715660000122
上式(2-14)、(2-15)中，σ(x(k))是微型相量测量装置D-PMU各个量测数据的标准差，可通过数据手册查到。x＝[V1.a,δ1.a,I12.a,θ12.a,I12.b,θ12.b,I12.c,θ12.c]，通过计算伪D-PMU量测的标准差，可以确定计算半径，即当伪量测标准差高于预设值时停止计算，从而确保伪量测的准确性。对于电流和功率也可以采用同样的方法来计算其对应的伪量测，在这里不在赘述。
实施例3
抗虚假数据注入攻击数学证明；从直观上理解，如果在状态估计算法中增加精度和准确性都很高的D-PMU数据，坏数据将得到“稀释”，从而在某种程度上增加算法的鲁棒性。下面对其进行数学证明。
对于一个信息有限的攻击者来说，需要将FDIA控制在某一特定的区域。如图4所示，将整个网络分为攻击区A和非攻击区N，区域A和N之间通过联络线连接。在发生虚假数据注入攻击之后，估计残差变为：
r′＝min||z′-H(x′)x′||2 (3-1)
式(3-1)中，z′＝[z′AzN]T，z′A表示攻击区量测，zN表示非攻击区量测。为了使攻击不被***检测出来，新的残差不能大于原来正常量测状态下的残差。
本实施例中，定义一个指标“攻击复杂度C”，通过攻击复杂度C来定量分析***的抗虚假数据攻击能力。攻击复杂度包括两个方面：(1)攻击者需要篡改的量测数据；(2)所篡改的数据与真实值的偏差。为了不失一般性，状态估计采用高斯-牛顿法进行迭代计算，迭代增量
Figure BDA0003045715660000131
为：
Figure BDA0003045715660000132
其中Hk为雅克比矩阵，
Figure BDA0003045715660000133
对上式进行线性化和代数操作可以得到：
Figure BDA0003045715660000134
Figure BDA0003045715660000135
上式(3-4)中，K叫做“帽矩阵”。帽矩阵的行与量测的冗余度相关：量测冗余度越高，帽矩阵中非零行就越多，算法在噪声和数据缺失的情况下表现较好。
由于帽矩阵的引入，量测残差r可以表示为：
Figure BDA0003045715660000136
虚假注入数据za可表示为za＝z+a。a是雅克比矩阵的线性组合，即a＝Hc。为了绕过坏数据检测***，a必须落在I–K的零空间内。攻击复杂度可以表示为：
min||a||0＝min||Hc||0 (3-6)
min||a||1＝min||Hc||1 (3-7)
st.Σi Hkici＝1
其中||Hc||0表示向量a中的非零行，||a||1＝Σi|ai|表示向量a的幅值。
将雅克比矩阵分解为
Figure BDA0003045715660000141
其中H1表示由前文中的伪D-PMU量测得到的雅克比矩阵，H2表示由其它量测得到的雅克比矩阵。增广矩阵H1将提高攻击复杂度的下限，这意味着完成一次成功的虚假数据注入攻击需要更多的量测被篡改，且注入的虚假量测偏差值要更大。
实施例4
最后结合测试***对本发明进行仿真测试；本实施例中采用IEEE123节点测试***进行仿真测试，网络拓扑如图5所示。此***包含123个节点，85个不平衡负载，4个并联电容器，1个变压器和11个开关，可以代表一个典型的中小型城市配电台区。仿真平台参数为：Intel Core i7 4.2GHz处理器，32GB内存，Win7终极版操作***，仿真软件为MatlabR2019a。
负荷超短期预测仿真测试情况：通过前述分析可知，***在每个时间断面均有准确的负荷节点数据，关键在于对包含n(n通常大于50)个高级量测体系AMI量测的低压节点负荷进行超短期预测。在此对节点55进行仿真测试。假设节点55包含100个高级量测体系AMI数据即用户，该数据是以真实值为期望，方差为0.05的高斯分布。
考虑以下四种情况：
CASE1：将最近一次AMI的量测数据直接相加得到所需要操作点的负荷数据；
CASE2：基于负荷趋势的新型超短期负荷预测法所提的基于负荷趋势的超短期负荷预测方法；
CASE3：基于谐波分量分析方法进行负荷预测，谐波分量阶数全部取3；
CASE4：基于谐波分量分析方法进行负荷预测，谐波分量阶数按照表1选取。
综上述4种方法对节点55的有功功率和无功功率超短期预测如图6所示。从图中可以看出，所提的基于谐波分量模型的预测方法对于低压节点的功率超短期预测效果较好。将CASE3与CASE4进行比较可以发现，对于谐波分量阶数的选取应根据操作点(OP)与最近采样点(LS)的时间间隔适当选取，选取的标准如表1所示。直接将最近一次AMI量测数据相加得到的节点负荷数据与真实值相差最大。
具体的抗虚假数据攻击能力测试情况：在如图5所示的IEEE 123节点测试***中，总共有18个D-PMU量测，8个SCADA量测点，设备安装位置在图上已标记。通过前述的计算，可以得到85个“伪D-PMU量测”，从而将量测冗余度从1.50提高到1.87。这些伪量测包括电压、电流和功率伪量测数据。值得注意的是，在计算伪D-PMU量测时，要对伪量测的标准差进行计算，当标准差高于预设值时停止计算，从而确保伪量测的准确性。
假设在SCADA量测节点18和67，线路57-60，低压终端节点46和71发生了虚假数据注入攻击，考虑两种情况：
CASE1：采用传统的加权最小二乘法(WLS)，包含85个伪D-PMU量测；
CASE2：采用传统的加权最小二乘法(WLS)，仅包含直接量测数据，不包含85个伪D-PMU量测。
如图7所示的节点58的电压和相角变化情况中可以看出，攻击从0.8秒时刻开始，到3秒时刻结束。在存在虚假数据攻击的情况下，传统的加权最小二乘法(WLS)对电压节点幅值和相角预测存在较大的误差(幅值误差最大值为2％，相角误差最大值为4％)。加入85个高精度的伪D-PMU量测之后，***的冗余度增加了25％，由于坏数据得到了“稀释”，状态估计的精度得到大幅度提升，尤其是对于相角的预测，其误差几乎为零。
通过对安装节点电压相量的估计，再加上线路参数以及其它量测数据，可以计算出线路潮流。表2分析了5条线路有功功率和无功功率的情况，包括潮流的实际值，量测值，采用WLS方法的预测值，所提方法的预测值，以及两种方法预测精度的差异(数据为正即表示所提方法预测精确度优于WLS法)。在表中可以看出加入伪D-PMU量测数据之后，线路的潮流计算精度得到提高。
表2五条线路潮流比较(kW/kVar)
Figure BDA0003045715660000161
Figure BDA0003045715660000171
图8表示存在虚假数据注入攻击的情况下，两种算法的一致性比较。状态估计算法的一致性用归一化的状态误差平方表示：
Figure BDA0003045715660000172
上式(4-1)中，
Figure BDA0003045715660000173
表示估计的误差协方差矩阵。把两种方法执行400次，可以看到所提方法的一致性指数集中在100-200之间，而不加伪D-PMU量测的一致性指数在200-400之间。
本发明通过处理微型同步相量量测装置D-PMU的电流相量数据，解决雅可比矩阵病态构造的问题；基于谐波分量分解方法对高级量测体系AMI量测的低压节点负荷进行超短期预测，得到了同一时间断面注入功率数据。通过对D-PMU临近节点或线路量测的电压、电流和功率进行计算，得到“伪D-PMU量测数据”，基于不确定性传播原理确定计算半径，从而提高量测冗余度，达到抗虚假数据攻击的目的。
上述实施例为本发明较佳的实现方案，除此之外，本发明还可以其它方式实现，在不脱离本技术方案构思的前提下任何显而易见的替换均在本发明的保护范围之内。
为了让本领域普通技术人员更方便地理解本发明相对于现有技术的改进之处，本发明的一些附图和描述已经被简化，并且为了清楚起见，本申请文件还省略了一些其它元素，本领域普通技术人员应该意识到这些省略的元素也可构成本发明的内容。

Claims (9)

1.配电网状态估计量测数据的处理方法，其特征在于：包括微型同步相量测量装置D-PMU的量测数据和高级量测体系AMI的量测数据的处理方法，通过对所述微型同步相量测量装置D-PMU的量测数据和高级量测体系AMI的量测数据进行处理，处理后的微型同步相量测量装置D-PMU的量测数据和高级量测体系AMI的量测数据与配网SCADA***的数据构成一个混合量测***来实现配电网状态估计。
2.根据权利要求1所述的配电网状态估计量测数据的处理方法，其特征在于：处理所述微型同步相量测量装置D-PMU的量测数据的步骤是：先所述确定微型同步相量测量装置D-PMU的安装节点，采集所述安装节点的电压相量数据和线路的电流相量数据；将采集到的线路的电流相量数据通过直角坐标的形式进行处理来获取雅可比矩阵中与电流相量的相关项；采集到的安装节点的电压相量数据用于计算伪D-PMU量测，从而提高量测冗余度。
3.根据权利要求1所述的配电网状态估计量测数据的处理方法，其特征在于：处理所述高级量测体系AMI的量测数据的步骤是：先确定操作点的数据与实际数据的最大时延；再采用超短期预测方法对操作点实际数据进行预测，从而得到同一时间断面注入功率数据。
4.根据权利要求2所述的配电网状态估计量测数据的处理方法，其特征在于：所述安装节点的电压相量数据为V＝V∠δ，其中，V表示安装节点的电压，δ表示安装节点的电压相量的角度。
5.根据权利要求2所述的配电网状态估计量测数据的处理方法，其特征在于：所述线路电流为：
I＝I∠θ
或I＝(V安装节点-V安装节点的相邻节点)·Y (1)
式(1)中，I表示安装节点与安装节点的相邻节点之间的线路电流，θ表示线路电流的角度，Y为导纳矩阵。
6.根据权利要求2所述的配电网状态估计量数据的处理方法，其特征在于：通过所述直角坐标形式处理的电流相量数据表示为：
I12.a＝Y11(V1.a-V2.a)+Y12(V1.b-V2.b)+Y13(V1.c-V2.c) (2)
式(2)中，I12.a表示安装节点与安装节点的相邻节点之间A相的线路电流，V1.a、V1.b、V1.c表示安装节点的A、B、C相的电压，V2.a、V2.b、V2.c表示安装节点的相邻节点的A、B、C相的电压。
7.根据权利要求2所述的配电网状态估计量测数据的处理方法，其特征在于：通过所述安装节点的电压相量数据计算与安装节点相邻的节点的电压相量为：
Figure RE-FDA0003126015690000021
与所述安装节点相邻的节点的电压相量V2即伪D-PMU量测，上式(3)中I12.a、I12.b、I12.c表示安装节点与安装节点的相邻节点之间线路A、B、C相的电流，z表示量测值，V1.a、V1.b、V1.c表示安装节点的A、B、C相的电压；
伪D-PMU量测在直角坐标系中表示为：
Figure RE-FDA0003126015690000022
上式(4)中，上标“r”表示相量的实部，“i”表示虚部；
基于不确定性传播原理分析伪D-PMU量测的标准差，所述伪D-PMU量测的标准差表示为：
Figure RE-FDA0003126015690000031
Figure RE-FDA0003126015690000032
上式(5)、(6)中，x＝[V1.a,δ1.a,I12.a,θ12.a,I12.b,θ12.b,I12.c,θ12.c]，σ(x(k))表示微型同步相量测量装置D-PMU各个量测数据的标准差；
最后通过计算所述伪D-PMU量测的标准差确定计算半径。
8.根据权利要求3所述的配电网状态估计量测数据的处理方法，其特征在于：所述超短期预测方法是将单个用户负荷数据看成正弦信号，所述正弦信号采用傅里叶级数形式的谐波分量模型表示。
9.根据权利要求3所述的配电网状态估计量测数据的处理方法，其特征在于：所述超短期预测表示为：
Figure RE-FDA0003126015690000033
上式(7)中，l′(last)表示最近一次采样的高级量测体系AMI的量测数据，参数an和bn通过历史数据计算得到，T表示操作点与最近一次数据更新的时间间隔，数据更新频率为96次/天，N表示谐波分量阶数。
CN202110471786.5A 2021-04-29 2021-04-29 配电网状态估计量测数据的处理方法 Pending CN113595058A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110471786.5A CN113595058A (zh) 2021-04-29 2021-04-29 配电网状态估计量测数据的处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110471786.5A CN113595058A (zh) 2021-04-29 2021-04-29 配电网状态估计量测数据的处理方法

Publications (1)

Publication Number Publication Date
CN113595058A true CN113595058A (zh) 2021-11-02

Family

ID=78243007

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110471786.5A Pending CN113595058A (zh) 2021-04-29 2021-04-29 配电网状态估计量测数据的处理方法

Country Status (1)

Country Link
CN (1) CN113595058A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114113898A (zh) * 2021-11-29 2022-03-01 大连海事大学 一种基于多源量测数据的配电网网损分析方法及***
CN117039893A (zh) * 2023-10-09 2023-11-10 国网天津市电力公司电力科学研究院 配电网状态确定方法、装置及电子设备
CN117060592A (zh) * 2023-08-30 2023-11-14 国网山东省电力公司德州供电公司 基于多域协同的电网数据同步校准方法及***

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114113898A (zh) * 2021-11-29 2022-03-01 大连海事大学 一种基于多源量测数据的配电网网损分析方法及***
CN114113898B (zh) * 2021-11-29 2023-11-14 大连海事大学 一种基于多源量测数据的配电网网损分析方法及***
CN117060592A (zh) * 2023-08-30 2023-11-14 国网山东省电力公司德州供电公司 基于多域协同的电网数据同步校准方法及***
CN117060592B (zh) * 2023-08-30 2024-02-23 国网山东省电力公司德州供电公司 基于多域协同的电网数据同步校准方法及***
CN117039893A (zh) * 2023-10-09 2023-11-10 国网天津市电力公司电力科学研究院 配电网状态确定方法、装置及电子设备
CN117039893B (zh) * 2023-10-09 2024-01-26 国网天津市电力公司电力科学研究院 配电网状态确定方法、装置及电子设备

Similar Documents

Publication Publication Date Title
Dehghanpour et al. A survey on state estimation techniques and challenges in smart distribution systems
Usman et al. Validation of a PMU‐based fault location identification method for smart distribution network with photovoltaics using real‐time data
CN113595058A (zh) 配电网状态估计量测数据的处理方法
Deng et al. False data injection attacks against state estimation in power distribution systems
CN109818349B (zh) 一种基于多维状态矩阵滑动匹配的电网鲁棒状态预测方法
Andreoni et al. Tri-objective optimal PMU placement including accurate state estimation: The case of distribution systems
Zanni et al. Probabilistic assessment of the process-noise covariance matrix of discrete Kalman filter state estimation of active distribution networks
Wang et al. Transmission lines positive sequence parameters estimation and instrument transformers calibration based on PMU measurement error model
Zhang et al. Topology identification method of distribution network based on smart meter measurements
Jamuna et al. Optimal placement of PMU and SCADA measurements for security constrained state estimation
Veerakumar et al. PMU-based real-time distribution system state estimation considering anomaly detection, discrimination and identification
Skok et al. Hybrid state estimation model based on PMU and SCADA measurements
Li et al. Dynamic topology awareness in active distribution networks under DG uncertainties using GMM-PSEs and KL divergence
Wiel et al. Identification of topology changes in power grids using phasor measurements
Rakpenthai et al. Measurement placement for power system state estimation using decomposition technique
Ghaedi et al. Modified WLS three-phase state estimation formulation for fault analysis considering measurement and parameter errors
Goh et al. Dynamic estimation of power system stability in different Kalman filter implementations
Janssen et al. Monitoring, protection and fault location in power distribution networks using system-wide measurements
Kirincic et al. A two‐step hybrid power system state estimator
Karvelis et al. State estimation using SCADA and PMU measurements for networks containing classic HVDC links
Jiang et al. Estimating electromechanical oscillation modes from synchrophasor measurements in bulk power grids using FSSI
Yang et al. Three‐phase feeder parameter estimation using radial basis function neural networks and multi‐run optimisation method with bad data preparation
Farantatos et al. Implementation of a 3-phase state estimation tool suitable for advanced distribution management systems
Ashok et al. Secondary network parameter estimation for distribution transformers
Williams et al. Off-line tracking of series parameters in distribution systems using AMI data

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination