CN106154037B - 一种基于校验的同步相量自适应计算方法 - Google Patents
一种基于校验的同步相量自适应计算方法 Download PDFInfo
- 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
- 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
- 238000012795 verification Methods 0.000 title claims abstract description 29
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 27
- 238000004364 calculation method Methods 0.000 title claims abstract description 15
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 53
- 238000005070 sampling Methods 0.000 claims abstract description 5
- 238000000034 method Methods 0.000 claims description 34
- 238000012937 correction Methods 0.000 claims description 7
- 230000035772 mutation Effects 0.000 claims description 6
- 230000003542 behavioural effect Effects 0.000 claims description 3
- 230000015556 catabolic process Effects 0.000 claims description 3
- 235000013399 edible fruits Nutrition 0.000 claims 1
- 230000005611 electricity Effects 0.000 claims 1
- 230000001052 transient effect Effects 0.000 abstract description 4
- 238000000605 extraction Methods 0.000 description 10
- 238000005259 measurement Methods 0.000 description 10
- 238000012360 testing method Methods 0.000 description 7
- 230000003044 adaptive effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000007689 inspection Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R25/00—Arrangements for measuring phase angle between a voltage and a current or between voltages or currents
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Phase Differences (AREA)
- Supply And Distribution Of Alternating Current (AREA)
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)=θ0+θ1t;动态算法中为更好反映信号的动态特征,假设相量模型中幅值和相角均为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)=θ0+θ1t。动态算法中为更好反映信号的动态特征,假设相量模型中幅值和相角均为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)=θ0+θ1t;动态算法中为更好反映信号的动态特征,假设相量模型中幅值和相角均为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)对应采样点实测值。
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 Active CN106154037B (zh) | 2016-08-11 | 2016-08-11 | 一种基于校验的同步相量自适应计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106154037B (zh) |
Families Citing this family (6)
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 | 国家电网公司 | 实现全时序电网量测数据同步采样及处理的测控方法 |
CN109374970B (zh) * | 2018-10-24 | 2021-08-13 | 南方电网科学研究院有限责任公司 | 实时校验的同步相量测量方法、装置、设备及存储介质 |
CN109521274B (zh) * | 2018-11-23 | 2020-09-25 | 南方电网科学研究院有限责任公司 | 一种同步相量测量方法、***、装置及可读存储介质 |
CN109995067B (zh) * | 2019-01-15 | 2022-02-08 | 东北电力大学 | 一种直流配电网换流站多功能备用电源控制方法 |
CN109991473B (zh) * | 2019-04-18 | 2021-09-24 | 南方电网科学研究院有限责任公司 | 导线电流相量的测量方法、测量装置及同步相量测量装置 |
Citations (7)
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的同步相量测量方法 |
-
2016
- 2016-08-11 CN CN201610657464.9A patent/CN106154037B/zh active Active
Patent Citations (7)
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)
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 |
---|---|---|
CN106154037B (zh) | 一种基于校验的同步相量自适应计算方法 | |
CN103245832B (zh) | 基于快速s变换的谐波时频特性参数估计方法及分析仪 | |
CN107271768B (zh) | 一种最小二乘拟合动态频率测量方法 | |
US10288667B2 (en) | Method and system for fault detection and faulted line identification in power systems using synchrophasors-based real-time state estimation | |
Abdelsalam et al. | Characterization of power quality disturbances using hybrid technique of linear Kalman filter and fuzzy-expert system | |
CN106093593B (zh) | 电容式电力设备介质损耗在线检测方法和*** | |
CN1892240B (zh) | 一种对避雷器泄漏电流中阻性电流值的在线监测方法 | |
CN106646140B (zh) | 基于测量波阻抗的高压直流输电线路区内外故障识别方法 | |
CN103401238B (zh) | 一种基于总体测辨法的电力负荷建模方法 | |
WO2016138764A1 (zh) | 一种提升pmu同步相量测量精度的方法 | |
CN103995178A (zh) | 一种基于时频聚集特性准则s变换的电压暂降检测方法 | |
CN102288804B (zh) | 一种基于正交变换的避雷器阻性电流计算方法 | |
WO2016065959A1 (zh) | 中性点不接地的10kV***中铁磁谐振的诊断方法 | |
CN106324340A (zh) | 一种同步相量和频率测量动态性能的方法 | |
CN104808035A (zh) | 一种基于广义双曲s变换的电压暂降检测方法 | |
CN109284933A (zh) | 一种基于数理统计的电子式互感器状态评估***及方法 | |
CN106199349B (zh) | 基于校准方波信号传递比的变压器局部放电区域判断方法 | |
CN102495264B (zh) | 基于信息熵的电压跌落持续时间检测方法 | |
Ma et al. | Harmonic and interharmonic analysis of mixed dense frequency signals | |
US20150094966A1 (en) | Quadrature-based voltage events detection method | |
Kolosok et al. | Detection of systematic errors in PMU measurements by the power system state estimation methods | |
CN106772193B (zh) | 一种利用电流互感器频率特性测量装置的测量方法 | |
Ghazizadeh-Ahsaee | Accurate NHIF locator utilizing two-end unsynchronized measurements | |
CN108982954B (zh) | 适用于馈线终端的计算相电压幅值与相位的方法及*** | |
CN106156489B (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 |