CN110048402B - 一种两阶段电网参数估计方法 - Google Patents
一种两阶段电网参数估计方法 Download PDFInfo
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
-
- 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
- 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
- Y02E60/00—Enabling 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)步骤中,所述的基于最优量测断面的可疑错误参数检测模型为:
式中,t为最优量测断面数;Pe为N维网络的错误参数向量;n为所有量测断面总数;λ∑i为错误参数的累积拉格朗日乘子,其计算公式为:
λ∑i=λi,1+λi,2+…+λi,N
2)步骤中所述的分段增广状态估计模型为:
ωi的计算公式:
实施上述的一种两阶段电网参数估计方法所用的两阶段电网参数估计系统,其特征在于:该系统包括电网参数采集单元和模型建立单元;电网参数采集单元与模型建立单元数据连接。
电网参数采集单元采集电网参数信息,模型建立单元建立基于最优量测断面的可疑错误参数检测模型,并检测出电网参数采集单元采集的电网参数信息中的可疑错误参数;
模型建立单元将可疑错误参数作为增广状态状态变量,建立分段增广状态估计模型,实现可疑错误参数的修正。
优点效果:
与现有技术相比,本发明具有以下优点:
一、计算精度高:本发明通过分段目标函数模型,降低错误参数之间的相互影响和减小不良数据对参数估计精度的影响,具有很高的参数估计精度;
二、具有很强的抗不良数据干扰的能力:相较于其他的参数辨识法在不良数据存在时,参数辨识的精度下降甚至不收敛,本发明可以在不良数据和多个参数错误同时存在时实现参数错误的检测和修正;
三、应用前景广阔:通过本发明的可疑错误参数检测方法,可以结合人工经验将临近生命周期的电网设备参数也作为可疑错误参数,从而很好地避免了错误参数的漏检;通过状态估计模型的分段处理,可以自动地消除错误参数之间、错误参数与错误量测之间的相互影响,简单便捷。因此本发明具有很好的应用前景。
附图说明
图1为本发明的方法流程图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
如图1所示,本发明提供一种两阶段参数估计方法,包括按顺序执行的下列步骤:
步骤1)建立基于最优量测断面的可疑错误参数检测模型,检测出电网中的可疑错误参数,最优量测断面检测模型为:
式中,t为最优量测断面数;Pe为N维网络的错误参数向量;n为所有量测断面总数;λ∑i为错误参数的累积拉格朗日乘子,其计算公式为:
λ∑i=λi,1+λi,2+…+λi,N
步骤2)将可疑错误参数作为增广状态状态变量,建立分段增广状态估计模型,实现可疑错误参数的修正,分段增广状态估计模型为:
ωi的计算公式:
实施上述的一种两阶段电网参数估计方法所用的两阶段电网参数估计系统,该系统包括电网参数采集单元和模型建立单元;电网参数采集单元与模型建立单元数据连接。
电网参数采集单元采集电网参数信息,模型建立单元建立基于最优量测断面的可疑错误参数检测模型,并检测出电网参数采集单元采集的电网参数信息中的可疑错误参数;
模型建立单元将可疑错误参数作为增广状态状态变量,建立分段增广状态估计模型,实现可疑错误参数的修正。
两阶段参数估计方法实施例
本申请算例在Microsoft Visual C++2010环境下编制了相应程序,采用IEEE9节点系统,以修正和辨识支路阻抗为例验证所提出的方法的可行性和有效性,算例所有量均采用标幺值表示。
9节点算例量测采用完全配置,量测数据使用潮流计算结果。设置的不良数据在潮流计算结果基础上加上该量测的3倍标准差。表1-表4的结果表明,本申请的方法对电抗估计更加敏感,电阻估计和电抗估计的收敛精度设置为0.001和0.0001。
表1(a)中设置了变压器有功和无功以及线路有功P45三个量测错误。由表1(b)可知,当系统中存在不良数据时,特别当不良数据和错误参数(线路2-7)强相关时,参数辨识的结果会受到不良数据的影响导致精度下降,甚至不收敛(由表4(a)和表4(b)),但是,该方法的收敛速度不受到影响而且不必对不良数据做出额外处理就能完成收敛。结果表明本发明可以自适应的消除非强相关量测对错误参数的影响,直接得到参数辨识的结果,而不需要先区分并剔除不良数据和参数错误导致的错误量测。但是,当不良数据与参数错误强相关时,可能导致精度下降甚至不收敛(由表1、表4),这一点需要进一步研究。
表1(a)强相关不良数据设置
表1(b)强相关不良数据参数估计
线路 | 真值 | 设定值 | 检测 | 估计值 |
L2-7 | 0+0.0625j | 0+0.0645j | L2-7 | -0.000003+0.062660j |
表2(a)非强相关不良数据设置
表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不能忽略)
表3(b)非强相关不良数据电阻估计
线路 | 真值 | 设定值 | 检测 | 估计值 |
<![CDATA[r<sub>45</sub>]]> | 0.01 | 0.011 | <![CDATA[r<sub>45</sub>]]> | 0.010559 |
表4(a)强相关不良数据设置(r不能忽略)
表4(b)强相关不良数据电阻辨识结果
线路 | 真值 | 设定值 | 检测 | 估计值 |
<![CDATA[r<sub>45</sub>]]> | 0.01 | 0.011 | <![CDATA[r<sub>45</sub>]]> | 不收敛 |
综上,与现有技术相比,本发明具有抗不良数据干扰、计算精度高、应用前景广阔等优点。
Claims (7)
6.实施权利要求1所述的一种两阶段电网参数估计方法所用的两阶段电网参数估计系统,其特征在于:该系统包括电网参数采集单元和模型建立单元;电网参数采集单元与模型建立单元数据连接。
7.根据权利要求6所述的两阶段电网参数估计系统,其特征在于:电网参数采集单元采集电网参数信息,模型建立单元建立基于最优量测断面的可疑错误参数检测模型,并检测出电网参数采集单元采集的电网参数信息中的可疑错误参数;
模型建立单元将可疑错误参数作为增广状态状态变量,建立分段增广状态估计模型,实现可疑错误参数的修正。
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110429587A (zh) * | 2019-07-19 | 2019-11-08 | 国网辽宁省电力有限公司大连供电公司 | 一种两阶段电网参数估计方法 |
Citations (6)
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 | 国电南瑞科技股份有限公司 | 一种基于拓扑分区的多断面联合参数估计方法 |
-
2019
- 2019-02-28 CN CN201910148284.1A patent/CN110048402B/zh active Active
Patent Citations (6)
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)
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 |