CN106154037A - 一种基于校验的同步相量自适应计算方法 - Google Patents
一种基于校验的同步相量自适应计算方法 Download PDFInfo
- Publication number
- CN106154037A CN106154037A CN201610657464.9A CN201610657464A CN106154037A CN 106154037 A CN106154037 A CN 106154037A CN 201610657464 A CN201610657464 A CN 201610657464A CN 106154037 A CN106154037 A CN 106154037A
- Authority
- CN
- China
- Prior art keywords
- verification
- algorithm
- synchronized phasor
- dynamic
- result
- 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.)
- Granted
Links
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
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 (4)
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),校验动态算法计算结果是否正确,若校验通过则使用动态算法的结果作为最终结果,如果校验未通过则可能故障发生在本数据窗内等动态算法同样无法处理的情况,为避免波动较大仍使用稳态算法结果。
2.根据权利要求1所述的基于校验的同步相量自适应计算方法,其特征在于上述步骤(4)的具体地校验过程为:通过计算得到的相量,反推各采样点的理论计算值,并与实际测得的点进行比较,若总偏差较小,则校验通过。
3.根据权利要求1所述的基于校验的同步相量自适应计算方法,其特征在于针对稳态算法,由校正后θ0、θ1和a,反推各抽取点理论计算值,抽取点的密度根据实际情况确定。
4.根据权利要求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 true CN106154037A (zh) | 2016-11-23 |
CN106154037B 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) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107144734A (zh) * | 2017-05-15 | 2017-09-08 | 北京理工大学 | 一种适用于pmu的配电网高精度相量测量方法 |
CN107656151A (zh) * | 2017-08-10 | 2018-02-02 | 国家电网公司 | 实现全时序电网量测数据同步采样及处理的测控方法 |
CN109374970A (zh) * | 2018-10-24 | 2019-02-22 | 南方电网科学研究院有限责任公司 | 实时校验的同步相量测量方法、装置、设备及存储介质 |
CN109521274A (zh) * | 2018-11-23 | 2019-03-26 | 南方电网科学研究院有限责任公司 | 一种同步相量测量方法、系统、装置及可读存储介质 |
CN109991473A (zh) * | 2019-04-18 | 2019-07-09 | 南方电网科学研究院有限责任公司 | 导线电流相量的测量方法、测量装置及同步相量测量装置 |
CN109995067A (zh) * | 2019-01-15 | 2019-07-09 | 东北电力大学 | 一种直流配电网换流站多功能备用电源控制方法 |
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 |
---|
符玲: "基于频域动态模型的同步相量测量算法", 《中国电机工程学报》 * |
谢小荣: "一种估计同步发电机功角的新方法", 《中国电机工程学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107144734A (zh) * | 2017-05-15 | 2017-09-08 | 北京理工大学 | 一种适用于pmu的配电网高精度相量测量方法 |
CN107144734B (zh) * | 2017-05-15 | 2019-09-27 | 北京理工大学 | 一种适用于pmu的配电网高精度相量测量方法 |
CN107656151A (zh) * | 2017-08-10 | 2018-02-02 | 国家电网公司 | 实现全时序电网量测数据同步采样及处理的测控方法 |
CN109374970A (zh) * | 2018-10-24 | 2019-02-22 | 南方电网科学研究院有限责任公司 | 实时校验的同步相量测量方法、装置、设备及存储介质 |
CN109521274A (zh) * | 2018-11-23 | 2019-03-26 | 南方电网科学研究院有限责任公司 | 一种同步相量测量方法、系统、装置及可读存储介质 |
CN109521274B (zh) * | 2018-11-23 | 2020-09-25 | 南方电网科学研究院有限责任公司 | 一种同步相量测量方法、系统、装置及可读存储介质 |
CN109995067A (zh) * | 2019-01-15 | 2019-07-09 | 东北电力大学 | 一种直流配电网换流站多功能备用电源控制方法 |
CN109995067B (zh) * | 2019-01-15 | 2022-02-08 | 东北电力大学 | 一种直流配电网换流站多功能备用电源控制方法 |
CN109991473A (zh) * | 2019-04-18 | 2019-07-09 | 南方电网科学研究院有限责任公司 | 导线电流相量的测量方法、测量装置及同步相量测量装置 |
CN109991473B (zh) * | 2019-04-18 | 2021-09-24 | 南方电网科学研究院有限责任公司 | 导线电流相量的测量方法、测量装置及同步相量测量装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106154037B (zh) | 2019-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106154037B (zh) | 一种基于校验的同步相量自适应计算方法 | |
CN103245832B (zh) | 基于快速s变换的谐波时频特性参数估计方法及分析仪 | |
CN106646140B (zh) | 基于测量波阻抗的高压直流输电线路区内外故障识别方法 | |
CN103401238B (zh) | 一种基于总体测辨法的电力负荷建模方法 | |
CN102435815B (zh) | 一种moa阻性电流在线监测系统的工作方法 | |
CN106324340B (zh) | 一种同步相量和频率测量动态性能的方法 | |
CN110007141B (zh) | 一种基于电压电流谐波相似度的谐振点检测方法 | |
CN102508026B (zh) | 一种电能质量谐波分析仪的谐波分析方法 | |
CN102288873A (zh) | 一种基于平波电感元件性能方程测后模拟识别直流输电线路区内外故障的方法 | |
CN103439566B (zh) | 一种moa阻性电流测试仪的工作方法 | |
CN104808035A (zh) | 一种基于广义双曲s变换的电压暂降检测方法 | |
CN102621388A (zh) | 基于同步时域信号的输电线路集中参数在线确定方法 | |
CN102955068A (zh) | 一种基于压缩采样正交匹配追踪的谐波检测方法 | |
CN105138843A (zh) | 电力系统采样飞点检测及其修复方法 | |
CN106054023A (zh) | 一种输电线路单端测距中估计两侧系统阻抗的方法 | |
CN102495264B (zh) | 基于信息熵的电压跌落持续时间检测方法 | |
CN106772193B (zh) | 一种利用电流互感器频率特性测量装置的测量方法 | |
CN102590690B (zh) | 一种测后模拟识别直流输电线路区内外故障的方法 | |
CN108982954B (zh) | 适用于馈线终端的计算相电压幅值与相位的方法及系统 | |
CN108051682B (zh) | 一种单相整流器系统阻抗模型的验证方法 | |
CN103630744B (zh) | Pmu相角检测方法及系统 | |
CN106156489B (zh) | 一种电力系统同步相量快速计算方法 | |
CN102435816B (zh) | 一种moa阻性电流测试仪的工作方法 | |
CN104730384A (zh) | 基于不完全s变换的电能扰动识别与定位方法 | |
CN105004920A (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 |