CN106154037B - 一种基于校验的同步相量自适应计算方法 - Google Patents

一种基于校验的同步相量自适应计算方法 Download PDF

Info

Publication number
CN106154037B
CN106154037B CN201610657464.9A CN201610657464A CN106154037B CN 106154037 B CN106154037 B CN 106154037B CN 201610657464 A CN201610657464 A CN 201610657464A CN 106154037 B CN106154037 B CN 106154037B
Authority
CN
China
Prior art keywords
algorithm
verification
dynamic
result
synchronized phasor
Prior art date
Application number
CN201610657464.9A
Other languages
English (en)
Other versions
CN106154037A (zh
Inventor
苏寅生
唐卓尧
李鹏
徐全
陆超
韩英铎
吴京涛
Original Assignee
中国南方电网有限责任公司
清华大学
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 中国南方电网有限责任公司, 清华大学 filed Critical 中国南方电网有限责任公司
Priority to CN201610657464.9A priority Critical patent/CN106154037B/zh
Publication of CN106154037A publication Critical patent/CN106154037A/zh
Application granted granted Critical
Publication of CN106154037B publication Critical patent/CN106154037B/zh

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R25/00Arrangements for measuring phase angle between a voltage and a current, or between voltages or currents

Abstract

本发明是一种基于校验的同步相量自适应计算方法,包括以下步骤:1)初始化,确定每周波采样点数、信号模型阶数及稳态算法中两数据窗的时间间隔;2)建立稳态和动态算法中电力信号的模型;3)通过两个数据窗的DFT结果,求解电力信号的同步相量;4)通过比较反推值和实测值,校验计算结果是否正确,校验通过结束,未通过进入5);5)采用动态算法计算同步相量;6)再次校验,若校验通过则使用动态算法的结果作为最终结果,若未通过则仍使用稳态算法结果作为最终结果。进而自适应切换计算模式,本发明同时兼顾稳态情况的计算精度及暂态情况的动态性能问题,兼顾了准确性和快速性要求。本发明可以广泛应用于电力系统同步相量计算中。

Description

一种基于校验的同步相量自适应计算方法
技术领域
本发明涉及电力系统自动测量技术领域,是一种基于校验的同步相量自适应计算方法,特别涉及一种高精度和快速响应的同步相量测量方法。
背景技术
近年来,以同步相量测量装置(phasor measurement unit,PMU)为基础的广域测量系统在电力系统动态过程监视、在线辨识、安全稳定分析以及广域控制等领域中得到广泛的应用。随着广域测量系统应用研究的不断深入,PMU装置对同步相量测量的要求越来越高,其相量算法的准确性和快速性将直接影响到相关应用功能的可靠性。
已有基于DFT(Discrete Fourier Transform,离散傅里叶变换)的修正算法,相对于传统DFT算法,较大的提高了计算精度,但对突变等动态过程响应速度有限,且在幅值时刻变化时难以满足精度要求。而已有基于频域动态模型的算法,利用同一数据窗不同频点滤波器的响应来修正DFT的估计结果,提高了对突变等动态过程的响应速度,但此算法存在噪声和谐波抑制能力弱,且在故障过程中精度低的问题。
发明内容
本发明的目的在于提供一种基于校验的同步相量自适应计算方法,该方法对稳态和动态情况,分别设计了稳态算法和动态算法两种不同计算模式,通过反推校验实现二者间的自适应切换,以同时兼顾稳态情况的计算精度及暂态情况的动态性能要求。本发明方法可在不同情况下,通过校验实现计算模式自适应切换,可以兼顾快速响应和高精度的要求。
本发明采用的技术方案是:本发明的基于校验的同步相量自适应计算方法,包括以下步骤:
(1)初始化,确定每周波采样点数N,稳态和动态算法的电力信号模型中幅值和相角的阶数K,稳态算法中两数据窗的间隔l点;
(2)信号建模,该方法采用复信号P(t)表示电力信号的动态相量为P(t)=a(t)ej θ(t),电力信号x(t)表示为:式中:a(t)和θ(t)分别表示电力信号幅值和相角的多项式;f0为额定频率;稳态算法模型假设计算周期内幅值恒定,频率偏差恒定,即a(t)=a,θ(t)=θ01t;动态算法中为更好反映信号的动态特征,假设相量模型中幅值和相角均为K阶模型,即
(3)通过两个数据窗数据的DFT(Discrete Fourier Transform,离散傅里叶变换)结果,求解电力信号的同步相量,
模型参数:
校正后相角:
校正后幅值:
式中:θc和θcl为校正前两数据窗分别计算角度,ac为校正前DFT计算的幅值,
(4)通过比较反推值和实测值,校验计算结果是否正确,如果校验通过则结束,若未通过则说明当前可能处于突变等动态过程,进入步骤(5);
(5)假设电力信号为K阶模型,采用已有动态算法求解同步相量,如通过基于同一数据窗各谐波含量计算同步相量的方法;
(6)同步骤(4),校验动态算法计算结果是否正确,若校验通过则使用动态算法的结果作为最终结果,如果校验未通过则可能故障发生在本数据窗内等动态算法同样无法处理的情况,为避免波动较大仍使用稳态算法结果。
上述步骤(4)的具体地校验过程为:通过计算得到的相量,反推各采样点的理论计算值,并与实际测得的点进行比较,若总偏差较小,则校验通过。
本发明针对稳态算法,由校正后θ0、θ1和a,反推各抽取点理论计算值,抽取点的密度根据实际情况确定。
本发明抽取点的理论计算值为:
式中:ni(i=1,2,···,L)为所抽取点,且1≤ni≤N;xc(ni)为所抽取点的理论计算值;
理论计算值与实测值之间的总偏差为:
式中:xm(ni)对应采样点实测值。
本发明由于采取以上技术方案,与现有技术相比,具有以下优点:本发明在同时考虑PMU算法快速性和准确性要求下,提出一种基于校验的同步相量自适应计算方法,在稳态和动态情况下自适应采用不同的计算模式。本发明测量方法对稳态和动态情况,分别设计了稳态算法和动态算法两种不同计算模式,通过反推校验实现二者间的自适应切换,从而同时兼顾稳态情况的计算精度及暂态情况的动态性能要求。本发明针对已有算法无法同时兼顾稳态情况的计算精度及暂态情况的动态性能问题,本发明方法可在不同情况下,通过校验实现计算模式自适应切换,可以兼顾快速响应和高精度的要求。
附图说明
图1为本发明方法的流程示意图。
图2为本发明方法实施例中的算法流程示意图。
图3为本发明方法的仿真测试中阶跃响应示意图。
具体实施方式
以下结合附图对本发明进行详细的描述。本发明提出来的基于校验的同步相量自适应计算方法,可以采用多种硬件方案来实现,在此不再赘述。本发明所提出的测量算法流程如图1所示,其中PMU算法的流程如图2所示。测量方法包括以下步骤:
(1)初始化,确定每周波采样点数N,稳态和动态算法的电力信号模型中幅值和相角的阶数K,稳态算法中两数据窗的间隔l点。
(2)信号建模,该方法采用复信号P(t)表示电力信号的动态相量为P(t)=a(t)ej θ(t),电力信号x(t)可以表示为:式中:a(t)和θ(t)分别表示电力信号幅值和相角的多项式;f0为额定频率。稳态算法模型假设计算周期内幅值恒定,频率偏差恒定,即a(t)=a,θ(t)=θ01t。动态算法中为更好反映信号的动态特征,假设相量模型中幅值和相角均为K阶模型,即
(3)通过两个数据窗数据的DFT结果,求解电力信号的同步相量。
模型参数:
校正后相角:
校正后幅值:
式中:θc和θcl为校正前两数据窗分别计算角度,ac为校正前DFT计算的幅值,
(4)通过比较反推值和实测值,校验计算结果是否正确,如果校验通过则结束,若未通过则说明当前可能处于突变等动态过程,进入(5);
具体地校验过程为:通过计算得到的相量,反推各采样点的理论计算值,并与实际测得的点进行比较,若总偏差较小,则校验通过。针对稳态算法,由校正后θ0、θ1和a,反推各抽取点理论计算值,抽取点的密度根据实际情况确定。抽取点的理论计算值为:
式中:ni(i=1,2,···,L)为所抽取点,且1≤ni≤N;xc(ni)为所抽取点的理论计算值。
理论计算值与实测值之间的总偏差为:
式中:xm(ni)对应采样点实测值。
(5)假设电力信号为K阶模型,采用已有动态算法求解同步相量,如通过基于同一数据窗各谐波含量计算同步相量的方法。
(6)同步骤(4),校验动态算法计算结果是否正确,若校验通过则使用动态算法的结果作为最终结果。如果校验未通过则可能故障发生在本数据窗内等动态算法同样无法处理的情况,为避免波动较大仍使用稳态算法结果。
为进一步说明本发明方法,下面对本发明方法进行仿真测试。仿真过程中算法的采样率为6400Hz等间隔采样。
1、频率偏差测试
为了测试相量测量算法在频率偏离额定频率时的性能,国标《电力系统同步相量测量装置检测规范》规定频率测量范围是45Hz至55Hz,在基波频率偏离额定值5Hz时,电压、电流幅值测量误差改变量应小于额定频率时测量误差极限值的100%,相角测量误差改变量应不大于1°。表1为本发明方法在频率偏离额定频率5Hz时的测试结果。可以看出,本发明方法的量测精度高于标准要求。
表1频率偏差为5Hz时相量测试结果
最大误差 平均绝对误差 均方根误差
角度误差 -0.02° 0.0174° 0.0175°
幅值误差 0.014% 0.0092% 0.0097%
2、为测试算法对突变等动态过程的响应性能,参考国标《电力系统同步相量测量装置检测规范》,施加90°相角阶跃信号:
本发明方法测试结果如图3所示,可以看出,本发明方法的阶跃响应时间为20ms明显快于国标规定的30ms。

Claims (2)

1.一种基于校验的同步相量自适应计算方法,其特征在于包括以下步骤:
(1)初始化,确定每周波采样点数N,稳态和动态算法的电力信号模型中幅值和相角的阶数K,稳态算法中两数据窗的间隔l点;
(2)信号建模,该方法采用复信号P(t)表示电力信号的动态相量为P(t)=a(t)ejθ(t),电力信号x(t)表示为:式中:a(t)和θ(t)分别表示电力信号幅值和相角的多项式;f0为额定频率;稳态算法模型假设计算周期内幅值恒定,频率偏差恒定,即a(t)=a,θ(t)=θ01t;动态算法中为更好反映信号的动态特征,假设相量模型中幅值和相角均为K阶模型,即
(3)通过两个数据窗数据的DFT(Discrete Fourier Transform,离散傅里叶变换)结果,求解电力信号的同步相量,
模型参数:
校正后相角:
校正后幅值:
式中:θc和θcl为校正前两数据窗分别计算角度,ac为校正前DFT计算的幅值,
(4)通过比较反推值和实测值,校验计算结果是否正确,如果校验通过则结束,若未通过则说明当前可能处于突变等动态过程,则进入步骤(5);
(5)假设电力信号为K阶模型,采用已有动态算法求解同步相量,所述动态算法是通过基于同一数据窗各谐波含量计算同步相量的方法;
(6)同步骤(4),校验动态算法计算结果是否正确,若校验通过则使用动态算法的结果作为最终结果,如果校验未通过则可能故障发生在本数据窗内等动态算法同样无法处理的情况,为避免波动较大仍使用稳态算法结果;
上述步骤(4)的具体地校验过程为:通过计算得到的相量,反推各采样点的反推值,并与实测值进行比较,若总偏差较小,则校验通过;
针对稳态算法,由校正后θ0、θ1和a,反推各采样点的反推值,采样点的密度根据实际情况确定。
2.根据权利要求1所述的基于校验的同步相量自适应计算方法,其特征在于采样点的反推值为:
式中:ni(i=1,2,···,L)为所采样点,且1≤ni≤N;xc(ni)为所采样点的反推值;
反推值与实测值之间的总偏差为:
式中:xm(ni)对应采样点实测值。
CN201610657464.9A 2016-08-11 2016-08-11 一种基于校验的同步相量自适应计算方法 CN106154037B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610657464.9A CN106154037B (zh) 2016-08-11 2016-08-11 一种基于校验的同步相量自适应计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610657464.9A CN106154037B (zh) 2016-08-11 2016-08-11 一种基于校验的同步相量自适应计算方法

Publications (2)

Publication Number Publication Date
CN106154037A CN106154037A (zh) 2016-11-23
CN106154037B true CN106154037B (zh) 2019-04-02

Family

ID=57329984

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610657464.9A CN106154037B (zh) 2016-08-11 2016-08-11 一种基于校验的同步相量自适应计算方法

Country Status (1)

Country Link
CN (1) CN106154037B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107144734B (zh) * 2017-05-15 2019-09-27 北京理工大学 一种适用于pmu的配电网高精度相量测量方法
CN107656151B (zh) * 2017-08-10 2020-12-29 国家电网公司 实现全时序电网量测数据同步采样及处理的测控方法
CN109374970A (zh) * 2018-10-24 2019-02-22 南方电网科学研究院有限责任公司 实时校验的同步相量测量方法、装置、设备及存储介质
CN109521274B (zh) * 2018-11-23 2020-09-25 南方电网科学研究院有限责任公司 一种同步相量测量方法、系统、装置及可读存储介质
CN109991473A (zh) * 2019-04-18 2019-07-09 南方电网科学研究院有限责任公司 导线电流相量的测量方法、测量装置及同步相量测量装置

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1477401A (zh) * 2003-07-18 2004-02-25 清华大学 高精度的同步相量测量方法
CN101261293A (zh) * 2007-03-08 2008-09-10 国网南京自动化研究院 基于自适应滤波器的电力稳态信号跟踪测量法
CN103543335A (zh) * 2013-10-30 2014-01-29 国家电网公司 一种同步相量测量方法
CN104049144A (zh) * 2014-06-27 2014-09-17 福州大学 一种滤除衰减直流分量的同步相量测量的实现方法
CN104502703A (zh) * 2014-12-13 2015-04-08 西南交通大学 基于频域动态模型的电力信号同步相量测量方法
CN104793053A (zh) * 2015-04-22 2015-07-22 福州大学 一种基于dft的同步相量相角测量方法
WO2015157989A1 (zh) * 2014-04-18 2015-10-22 华北电力大学 一种适用于p类相量测量单元pmu的同步相量测量方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1477401A (zh) * 2003-07-18 2004-02-25 清华大学 高精度的同步相量测量方法
CN101261293A (zh) * 2007-03-08 2008-09-10 国网南京自动化研究院 基于自适应滤波器的电力稳态信号跟踪测量法
CN103543335A (zh) * 2013-10-30 2014-01-29 国家电网公司 一种同步相量测量方法
WO2015157989A1 (zh) * 2014-04-18 2015-10-22 华北电力大学 一种适用于p类相量测量单元pmu的同步相量测量方法
CN104049144A (zh) * 2014-06-27 2014-09-17 福州大学 一种滤除衰减直流分量的同步相量测量的实现方法
CN104502703A (zh) * 2014-12-13 2015-04-08 西南交通大学 基于频域动态模型的电力信号同步相量测量方法
CN104793053A (zh) * 2015-04-22 2015-07-22 福州大学 一种基于dft的同步相量相角测量方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
一种估计同步发电机功角的新方法;谢小荣;《中国电机工程学报》;20031130;第23卷(第11期);106-110
基于频域动态模型的同步相量测量算法;符玲;《中国电机工程学报》;20150320;第35卷(第6期);1371-1378

Also Published As

Publication number Publication date
CN106154037A (zh) 2016-11-23

Similar Documents

Publication Publication Date Title
Zhan et al. A Clarke transformation-based DFT phasor and frequency algorithm for wide frequency range
CN102435844B (zh) 一种频率无关的正弦信号相量计算方法
CN102255291B (zh) 一种基于贝杰龙模型的交流输电线路纵联保护的测后模拟方法
CN103869162B (zh) 一种基于时域准同步的动态信号相量测量方法
CN102998545B (zh) 一种变压器绕组工作状态的在线监测方法
CN101806832B (zh) 一种低频率信号的频率测量方法
Abdelsalam et al. Characterization of power quality disturbances using hybrid technique of linear Kalman filter and fuzzy-expert system
Shaw et al. A Kalman-filter spectral envelope preprocessor
CN104049144A (zh) 一种滤除衰减直流分量的同步相量测量的实现方法
CN102135569B (zh) 基于波动量法的用户侧谐波发射水平实用化的估计方法
Dasgupta et al. Line parameter estimation using phasor measurements by the total least squares approach
Karimzadeh et al. Method for determining utility and consumer harmonic contributions based on complex independent component analysis
CN104569481B (zh) 瓦斯继电器油流流速采集系统和重瓦斯整定值校验方法
CN102749521A (zh) 一种电力系统谐波阻抗计算方法
Mai et al. Dynamic phasor and frequency estimator for phasor measurement units
CN103439567B (zh) 一种moa阻性电流在线监测系统的工作方法
Ji et al. Disturbance detection, location and classification in phase space
CN102830319B (zh) 一种氧化锌避雷器绝缘状态带电检测装置及方法
CN1892240B (zh) 一种对避雷器泄漏电流中阻性电流值的在线监测方法
CN102163844B (zh) 基于相量测量装置的电力系统状态检测方法
CN103353550A (zh) 一种测量电力系统信号频率及谐波参数的方法
WO2013079933A1 (en) Fault location in power distribution systems
CN105841762B (zh) 超声波水表的流量计量方法和系统
Ren et al. An adaptive phasor estimator for power system waveforms containing transients
CN103941072B (zh) 一种基于实数强跟踪滤波的电力信号突变参数测量方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant