CN110048402B - 一种两阶段电网参数估计方法 - Google Patents

一种两阶段电网参数估计方法 Download PDF

Info

Publication number
CN110048402B
CN110048402B CN201910148284.1A CN201910148284A CN110048402B CN 110048402 B CN110048402 B CN 110048402B CN 201910148284 A CN201910148284 A CN 201910148284A CN 110048402 B CN110048402 B CN 110048402B
Authority
CN
China
Prior art keywords
power grid
parameter
measurement
error
grid parameter
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.)
Active
Application number
CN201910148284.1A
Other languages
English (en)
Other versions
CN110048402A (zh
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.)
State Grid Corp of China SGCC
State Grid Liaoning Electric Power Co Ltd
Electric Power Research Institute of State Grid Liaoning Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Liaoning Electric Power Co Ltd
Electric Power Research Institute of State Grid Liaoning Electric Power Co Ltd
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 State Grid Corp of China SGCC, State Grid Liaoning Electric Power Co Ltd, Electric Power Research Institute of State Grid Liaoning Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Publication of CN110048402A publication Critical patent/CN110048402A/zh
Application granted granted Critical
Publication of CN110048402B publication Critical patent/CN110048402B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

一种两阶段电网参数估计方法,包括以下步骤:建立基于最优量测断面的可疑错误参数检测模型,检测出电网中的可疑错误参数;将可疑错误参数作为增广状态状态变量,建立分段增广状态估计模型,实现可疑错误参数的修正本发明计算精度高:通过分段目标函数模型,降低错误参数之间的相互影响和减小不良数据对参数估计精度的影响,具有很高的参数估计精度;本发明可以在不良数据和多个参数错误同时存在时实现参数错误的检测和修正;可以自动地消除错误参数之间、错误参数与错误量测之间的相互影响,简单便捷。因此本发明具有很好的应用前景。

Description

一种两阶段电网参数估计方法
技术领域
本发明涉及电力系统参数辨识技术领域,尤其是涉及一种两阶段参数估计方法。
背景技术
随着电网的不断扩大,电力系统模型变得更加复杂,各类系统参数的数量也爆炸性增加,系统的正常工作对参数的准确性提出了更高的要求。同时,由于周期性电网检修及改造导致参数的人为设置错误和各种设备的偶发性参数漂移,导致基于系统参数的高级应用产生较大计算误差,甚至不收敛。因此,如何对网络参数错误准确识别并修正是电力系统精确建模中的重要问题。
现有参数辨识方法大致可以分为公式法、实测法、试探法、单线路模型法、灵敏度分析法和增广变量法六类。
公式法。公式法的基本思想是按照电力计算手册中的经验公式,通过线路长度、材质等多种因素综合计算。这种计算通常在线路建成初期进行测算,由于线路参数受到多种地形气候和设备老化的影响,因此按照经验公式得到的理论计算值一般只适用于电网规划,并不能满足系统实时运行对线路参数准确性的要求。
实测法。该方法主要利用测量仪器对线路参数进行实地测量,其测量方式分两种:一种是直接测正序和零序参数,另一种则是测量各相参数后,再经相序将相参数转换为序参数。该类方法通常需要进行停电实验或施加异频电源,无法避免复杂的人工记录和工作量,并且无法及时发现发生漂移的参数。
试探法。试探法是在已知可疑参数的基础上,通过特定步长改变参数后观察状态估计指标的变化情况,从而确定参数是欠修正还是过修正,通过步长的调整得到最终使得状态估计指标最好时相应的参数即为最终估计的参数值。该类方法虽然操作简单,但其最优修正步长难以准确获得,无法避免冗余工作。
单线路模型法。该方法是以单条线路的阻抗为未知参数对单条线路进行建模并利用该线路的相关量测数据实现其参数的辨识。此类方法的缺点在于由于中国输电网尚未实现PMU的完全覆盖,限制了其工程应用。
灵敏度分析法。基本思想是利用状态估计结果建立参数误差与量测误差之间的灵敏度矩阵,从而根据量测误差的估计即可得到参数误差的大小。当错误参数过多时该类方法的参数估计精度很容易受到不良数据和错误参数之间的相互影响,严重时会导致参数误修正。
增广变量法。其基本思想是将可疑参数作为增广状态变量进行状态估计。该类方法由于增加了状态量的维数,使得系统的冗余度下降,数值稳定性降低,严重时甚至导致不收敛。
发明内容
发明目的:
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种“检测+辨识”的两阶段参数辨识估计方法。
技术方案:
本发明的目的通过以下技术方案来实现:
一种两阶段参数估计方法,包括以下步骤:
1)建立基于最优量测断面的可疑错误参数检测模型,检测出电网中的可疑错误参数;
2)将可疑错误参数作为增广状态状态变量,建立分段增广状态估计模型,实现可疑错误参数的修正。
1)步骤中,所述的基于最优量测断面的可疑错误参数检测模型为:
Figure BDA0001980729980000021
式中,t为最优量测断面数;Pe为N维网络的错误参数向量;n为所有量测断面总数;λ∑i为错误参数的累积拉格朗日乘子,其计算公式为:
λ∑i=λi,1i,2+…+λi,N
Figure BDA0001980729980000022
式中,xi为第i个量测断面的状态量;Ji(xi,Pe)为第i个量测断面的目标函数;λi,1i,2+…+λi,N为向量λi中的各个元素;
Figure BDA0001980729980000031
为量测方程对量测误差Pe的雅格比矩阵的转置。
2)步骤中所述的分段增广状态估计模型为:
Figure BDA0001980729980000032
式中,hi(x,Pe)为含状态变量和增广状态变量的量测计算方程;k为迭代步骤;zi为第i个量测,ωi为其权重;x为电压幅值和相角构成的状态变量向量;σi为量测标准差;
Figure BDA0001980729980000033
为量测窗宽。
Figure BDA0001980729980000034
的计算公式:
Figure BDA0001980729980000035
式中,m为量测数。
ωi的计算公式:
Figure BDA0001980729980000036
式中,m为量测数。
实施上述的一种两阶段电网参数估计方法所用的两阶段电网参数估计系统,其特征在于:该系统包括电网参数采集单元和模型建立单元;电网参数采集单元与模型建立单元数据连接。
电网参数采集单元采集电网参数信息,模型建立单元建立基于最优量测断面的可疑错误参数检测模型,并检测出电网参数采集单元采集的电网参数信息中的可疑错误参数;
模型建立单元将可疑错误参数作为增广状态状态变量,建立分段增广状态估计模型,实现可疑错误参数的修正。
优点效果:
与现有技术相比,本发明具有以下优点:
一、计算精度高:本发明通过分段目标函数模型,降低错误参数之间的相互影响和减小不良数据对参数估计精度的影响,具有很高的参数估计精度;
二、具有很强的抗不良数据干扰的能力:相较于其他的参数辨识法在不良数据存在时,参数辨识的精度下降甚至不收敛,本发明可以在不良数据和多个参数错误同时存在时实现参数错误的检测和修正;
三、应用前景广阔:通过本发明的可疑错误参数检测方法,可以结合人工经验将临近生命周期的电网设备参数也作为可疑错误参数,从而很好地避免了错误参数的漏检;通过状态估计模型的分段处理,可以自动地消除错误参数之间、错误参数与错误量测之间的相互影响,简单便捷。因此本发明具有很好的应用前景。
附图说明
图1为本发明的方法流程图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
如图1所示,本发明提供一种两阶段参数估计方法,包括按顺序执行的下列步骤:
步骤1)建立基于最优量测断面的可疑错误参数检测模型,检测出电网中的可疑错误参数,最优量测断面检测模型为:
Figure BDA0001980729980000041
式中,t为最优量测断面数;Pe为N维网络的错误参数向量;n为所有量测断面总数;λ∑i为错误参数的累积拉格朗日乘子,其计算公式为:
λ∑i=λi,1i,2+…+λi,N
Figure BDA0001980729980000042
式中,xi为第i个量测断面的状态量;Ji(xi,Pe)为第i个量测断面的目标函数;λi,1i,2+…+λi,N为向量λi中的各个元素;
Figure BDA0001980729980000043
为量测方程对量测误差Pe的雅格比矩阵的转置。
步骤2)将可疑错误参数作为增广状态状态变量,建立分段增广状态估计模型,实现可疑错误参数的修正,分段增广状态估计模型为:
Figure BDA0001980729980000051
式中,hi(x,Pe)为含状态变量和增广状态变量的量测计算方程;k为迭代步骤;zi为第i个量测,ωi为其权重;x为电压幅值和相角构成的状态变量向量;σi为量测标准差;
Figure BDA0001980729980000052
为量测窗宽。
Figure BDA0001980729980000053
的计算公式:
Figure BDA0001980729980000054
式中,m为量测数。
ωi的计算公式:
Figure BDA0001980729980000055
式中,m为量测数。
实施上述的一种两阶段电网参数估计方法所用的两阶段电网参数估计系统,该系统包括电网参数采集单元和模型建立单元;电网参数采集单元与模型建立单元数据连接。
电网参数采集单元采集电网参数信息,模型建立单元建立基于最优量测断面的可疑错误参数检测模型,并检测出电网参数采集单元采集的电网参数信息中的可疑错误参数;
模型建立单元将可疑错误参数作为增广状态状态变量,建立分段增广状态估计模型,实现可疑错误参数的修正。
两阶段参数估计方法实施例
本申请算例在Microsoft Visual C++2010环境下编制了相应程序,采用IEEE9节点系统,以修正和辨识支路阻抗为例验证所提出的方法的可行性和有效性,算例所有量均采用标幺值表示。
9节点算例量测采用完全配置,量测数据使用潮流计算结果。设置的不良数据在潮流计算结果基础上加上该量测的3倍标准差。表1-表4的结果表明,本申请的方法对电抗估计更加敏感,电阻估计和电抗估计的收敛精度设置为0.001和0.0001。
表1(a)中设置了变压器有功
Figure BDA0001980729980000061
和无功
Figure BDA0001980729980000062
以及线路有功P45三个量测错误。由表1(b)可知,当系统中存在不良数据时,特别当不良数据
Figure BDA0001980729980000063
和错误参数(线路2-7)强相关时,参数辨识的结果会受到不良数据的影响导致精度下降,甚至不收敛(由表4(a)和表4(b)),但是,该方法的收敛速度不受到影响而且不必对不良数据做出额外处理就能完成收敛。结果表明本发明可以自适应的消除非强相关量测对错误参数的影响,直接得到参数辨识的结果,而不需要先区分并剔除不良数据和参数错误导致的错误量测。但是,当不良数据与参数错误强相关时,可能导致精度下降甚至不收敛(由表1、表4),这一点需要进一步研究。
表1(a)强相关不良数据设置
Figure BDA0001980729980000065
表1(b)强相关不良数据参数估计
线路 真值 设定值 检测 估计值
L2-7 0+0.0625j 0+0.0645j L2-7 -0.000003+0.062660j
表2(a)非强相关不良数据设置
Figure BDA0001980729980000064
表2(b)非强相关不良数据参数估计
线路 真值 设定值 检测 估计值
<![CDATA[L<sub>2-7</sub>]]> 0+0.0625j 0+0.0645j <![CDATA[L<sub>2-7</sub>]]> -0.000002+0.062503
表3(a)非强相关不良数据设置(r不能忽略)
Figure BDA0001980729980000071
表3(b)非强相关不良数据电阻估计
线路 真值 设定值 检测 估计值
<![CDATA[r<sub>45</sub>]]> 0.01 0.011 <![CDATA[r<sub>45</sub>]]> 0.010559
表4(a)强相关不良数据设置(r不能忽略)
Figure BDA0001980729980000072
表4(b)强相关不良数据电阻辨识结果
线路 真值 设定值 检测 估计值
<![CDATA[r<sub>45</sub>]]> 0.01 0.011 <![CDATA[r<sub>45</sub>]]> 不收敛
综上,与现有技术相比,本发明具有抗不良数据干扰、计算精度高、应用前景广阔等优点。

Claims (7)

1.一种两阶段电网参数估计方法,其特征在于:该方法包括以下步骤:
1)建立基于最优量测断面的可疑错误参数检测模型,检测出电网中的可疑错误参数;
2)将可疑错误参数作为增广状态状态变量,建立分段增广状态估计模型,实现可疑错误参数的修正;
2)步骤中所述的分段增广状态估计模型为:
Figure FDA0004059811450000011
式中,hi(x,Pe)为含状态变量和增广状态变量的量测计算方程;k为迭代步骤;zi为第i个量测,ωi为其权重;x为电压幅值和相角构成的状态变量向量;σi为量测标准差;
Figure FDA0004059811450000012
为量测窗宽。
2.根据权利要求1所述的一种两阶段电网参数估计方法,其特征在于:1)步骤中,所述的基于最优量测断面的可疑错误参数检测模型为:
Figure FDA0004059811450000013
式中,t为最优量测断面数;Pe为N维网络的错误参数向量;n为所有量测断面总数;λ∑i为错误参数的累积拉格朗日乘子。
3.根据权利要求2所述的一种两阶段电网参数估计方法,其特征在于:λ∑i计算公式为:
λ∑i=λi,1i,2+…+λi,N
Figure FDA0004059811450000014
式中,xi为第i个量测断面的状态量;Ji(xi,Pe)为第i个量测断面的目标函数;λi,1i,2+…+λi,N为向量λi中的各个元素;
Figure FDA0004059811450000015
为量测方程对量测误差Pe的雅格比矩阵的转置。
4.根据权利要求1所述的一种两阶段电网参数估计方法,其特征在于:
Figure FDA0004059811450000016
的计算公式:
Figure FDA0004059811450000021
式中,m为量测数。
5.根据权利要求1所述的一种两阶段电网参数估计方法,其特征在于:ωi的计算公式:
Figure FDA0004059811450000022
式中,m为量测数。
6.实施权利要求1所述的一种两阶段电网参数估计方法所用的两阶段电网参数估计系统,其特征在于:该系统包括电网参数采集单元和模型建立单元;电网参数采集单元与模型建立单元数据连接。
7.根据权利要求6所述的两阶段电网参数估计系统,其特征在于:电网参数采集单元采集电网参数信息,模型建立单元建立基于最优量测断面的可疑错误参数检测模型,并检测出电网参数采集单元采集的电网参数信息中的可疑错误参数;
模型建立单元将可疑错误参数作为增广状态状态变量,建立分段增广状态估计模型,实现可疑错误参数的修正。
CN201910148284.1A 2018-12-31 2019-02-28 一种两阶段电网参数估计方法 Active CN110048402B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN2018116517393 2018-12-31
CN201811651739 2018-12-31

Publications (2)

Publication Number Publication Date
CN110048402A CN110048402A (zh) 2019-07-23
CN110048402B true CN110048402B (zh) 2023-04-07

Family

ID=67274247

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910148284.1A Active CN110048402B (zh) 2018-12-31 2019-02-28 一种两阶段电网参数估计方法

Country Status (1)

Country Link
CN (1) CN110048402B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110429587A (zh) * 2019-07-19 2019-11-08 国网辽宁省电力有限公司大连供电公司 一种两阶段电网参数估计方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102280877A (zh) * 2011-07-25 2011-12-14 清华大学 一种多量测断面的电力系统不良支路参数辨识方法
WO2011156799A2 (en) * 2010-06-11 2011-12-15 Abb Research Ltd. Detecting state estimation network model data errors
CN103944165A (zh) * 2014-04-28 2014-07-23 国电南瑞科技股份有限公司 一种大电网参数辨识估计方法
CN104270119A (zh) * 2014-09-22 2015-01-07 衢州学院 基于非线性未知随机偏差的两阶段容积卡尔曼滤波方法
CN104836223A (zh) * 2014-11-14 2015-08-12 浙江大学 电网参数错误与不良数据协同辨识与估计方法
CN104992010A (zh) * 2015-06-25 2015-10-21 国电南瑞科技股份有限公司 一种基于拓扑分区的多断面联合参数估计方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011156799A2 (en) * 2010-06-11 2011-12-15 Abb Research Ltd. Detecting state estimation network model data errors
CN102280877A (zh) * 2011-07-25 2011-12-14 清华大学 一种多量测断面的电力系统不良支路参数辨识方法
CN103944165A (zh) * 2014-04-28 2014-07-23 国电南瑞科技股份有限公司 一种大电网参数辨识估计方法
CN104270119A (zh) * 2014-09-22 2015-01-07 衢州学院 基于非线性未知随机偏差的两阶段容积卡尔曼滤波方法
CN104836223A (zh) * 2014-11-14 2015-08-12 浙江大学 电网参数错误与不良数据协同辨识与估计方法
CN104992010A (zh) * 2015-06-25 2015-10-21 国电南瑞科技股份有限公司 一种基于拓扑分区的多断面联合参数估计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
拉格朗日乘子法电力系统网络参数错误辨识的应用;郭烨;《中国电机工程学报》;20131231;第43-48页 *
非线性输出误差模型的两阶段递推辨识算法;贾杰;《南昌航空大学学报(自然科学版)》;20141231;第24-30页 *

Also Published As

Publication number Publication date
CN110048402A (zh) 2019-07-23

Similar Documents

Publication Publication Date Title
CN107016236B (zh) 基于非线性量测方程的电网假数据注入攻击检测方法
CN107643507B (zh) 一种基于电网线路运行误差远程校准的精益线损分析与管控方法
CN104836223B (zh) 电网参数错误与不良数据协同辨识与估计方法
CN107453484B (zh) 一种基于wams信息的scada数据校准方法
CN102778631B (zh) 基于分压补偿的传感电缆泄漏检测精确定位方法
CN104992010A (zh) 一种基于拓扑分区的多断面联合参数估计方法
CN111046327A (zh) 适用于低频振荡与次同步振荡辨识的Prony分析方法
CN103324858A (zh) 配电网三相潮流状态估计方法
CN110048402B (zh) 一种两阶段电网参数估计方法
CN104537581B (zh) 采用模糊相似度匹配的电压暂降源在线定位方法
CN100478650C (zh) 改善发电设备故障检测的方法和装置
CN102636706B (zh) 一种电网中参数误差支路的辨识方法
CN111999691B (zh) 一种计量传感器装置的误差校准方法和误差校准装置
CN116520080A (zh) 基于时差矩阵数据拟合的配电网多分支故障行波定位方法和装置
CN112989587B (zh) 一种电容式电压互感器劣化成因在线分析方法及系统
CN105407496B (zh) 一种识别无线传感器网络中错误测量值的方法
CN109901022A (zh) 基于同步量测数据的配电网区域定位方法
CN113189513A (zh) 一种基于纹波的冗余电源均流状态识别方法
CN109193639B (zh) 一种电力系统抗差估计方法
CN114942402B (zh) 一种异常电能表定位方法及系统
CN116718161A (zh) 基于强化学习算法的架空线路弧垂状态估计系统及方法
CN115795222A (zh) 基于同步相量测量的电网不良参数辨识和校正方法
CN106355029B (zh) 一种配电系统关键量测解耦快速辨识的方法
CN105139286A (zh) 基于无功不平衡量的scada数据有效性评估方法
CN111797545B (zh) 一种基于实测数据的风电机组偏航折减系数计算方法

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
GR01 Patent grant
GR01 Patent grant