CN104049144B - 一种滤除衰减直流分量的同步相量测量的实现方法 - Google Patents
一种滤除衰减直流分量的同步相量测量的实现方法 Download PDFInfo
- Publication number
- CN104049144B CN104049144B CN201410298815.2A CN201410298815A CN104049144B CN 104049144 B CN104049144 B CN 104049144B CN 201410298815 A CN201410298815 A CN 201410298815A CN 104049144 B CN104049144 B CN 104049144B
- Authority
- CN
- China
- Prior art keywords
- delta
- rsqb
- lsqb
- centerdot
- omega
- 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
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明涉及电力系统同步相量测量技术领域,特别是一种滤除衰减直流分量的同步相量测量方法。基于传统离散傅里叶(DFT)算法,用复合梯形公式和线性插值法分别计算三个连续数据窗基波分量的实部和虚部,再计算实部和虚部的误差估计值,利用所得误差估计值对原实部和虚部进行修正,滤除电力系统中衰减直流分量给同步相量测量方法带来的影响。该方法在电力系统静态条件和动态条件下都拥有较高的同步相量测量精度和抗噪性,并且滤除了电力系统中衰减直流分量带来的影响,提高了同步相量测量方法在电力系统发生故障时的测量精度。
Description
技术领域
本发明涉及电力系统同步相量测量技术领域,特别是一种滤除衰减直流分量的同步相量测量方法。
背景技术
随着全球电力市场和区域电网工程的发展,电网的运行环境变得日渐复杂,新的安全隐患也随之产生,在电力系统动态条件下对电网进行实时地监控变得尤为重要。目前,广域测量系统(WAMS)作为一种新型、高效电网动态监控系统,为电网动态安全监控提供了新的技术保障。WAMS采集的数据具有同步性,同时,采集数据的处理结果要求相量数据带有精确的时标,所以同步相量测量技术是WAMS得以实现的基础,且同步相量测量方法是同步相量测量技术的核心部分,测量的精度将直接对电力系统故障分析、继电保护和稳定控制等应用的准确性造成影响。因此,同步相量测量方法逐渐成为近几年来研究的热点。
目前,电力系统相量的测量方法有过零点检测法、瞬时值法、Prony法、小波变换法、卡尔曼滤波法和DFT法等。过零点检测法容易受随机噪声和信号过零点处谐波影响而且实时性不好,易受系统动态条件的影响,测量误差不易控制。瞬时值法要求信号的波形为标准工频正弦波,对输入波形要求较高,通用性不强,且计算量较大。Prony法无法反映出系统动态条件下的非平稳性,且噪声对其拟合的结果影响很大,当噪声信噪比小于40dB,会得到不正确的结果。小波变换法不能在以额定频率为中心频率的窄带信号的同步相量测量中得到很好的体现,且由于小波分析运算复杂会影响同步相量测量的速度以及增加了数字信号处理器(DSP)的负担。卡尔曼滤波法在同步相量的相角测量中,存在着较大的误差,无法满足同步相量相角测量的要求。DFT法在电力系统动态条件下,如信号频率发生偏移时,会由于频率泄露的原因而产生较大的误差。但是由于在抑制谐波上拥有的优势,使现在绝大多数实际的同步相量测量方法都是以DFT法为基础的。
发明内容
本发明的目的在于提供一种滤除衰减直流分量的同步相量测量方法。该方法不仅能在电力系统静态﹑动态条件下都拥有较高的同步相量测量精度和抗噪性,而且对电力系统中衰减直流分量带来的误差能进行有效的滤除,提高同步相量测量方法在电力系统发生故障时的测量精度。
为实现上述目的,本发明的技术方案是:用积分公式和欧拉公式将电力信号转换为复数形式,用连续傅里叶积分表示出电力信号复数形式中的基波分量,采用复合梯形公式和线性插值法分别计算三个连续数据窗基波分量的实部和虚部,再计算由于衰减直流分量而产生的相量实部和虚部的误差估计值,利用所得误差估计值对原实部和虚部进行修正。其具体步骤如下:
步骤1:给定电力信号x(t):
式中,X0为直流分量;τ为衰减时间常数;X(n),分别为n次谐波的幅值和初相角,其中fs为采样频率,f0为额定频率;
步骤2:将电力信号x(t)用傅里叶级数表示:
式中,傅里叶系数
步骤3:设数据窗个数为M,且每个数据窗长度为N+2;设m为数据窗序号,取m=1;用积分公式将电力信号x(t)傅里叶级数表达式中的傅里叶系数转换为如下离散形式,从离散形式中可获得各采样点的值:
步骤4:用欧拉公式将电力信号的傅里叶级数形式转换为复数形式:
式中,x(t)包含多次谐波,而系数n=1,2,3...;
步骤5:用连续傅里叶积分表示电力信号复数形式中的基波分量:
式中,ΔN为分数,为当实际工作频率f偏离额定频率f0时,数据窗长度发生的偏移,f通过频率跟踪方法获得;
步骤6:用复合梯形公式和线性插值法分别计算第m、m+1、m+2三个连续数据窗中基波分量连续傅里叶积分表达式中相量的实部和虚部;
式中,Δt为采样时间间隔,ω0为额定角频率;
步骤7:设通过计算得出由于衰减直流分量而产生的相量实部误差估计值Δδa和虚部误差估计值Δδb:
式中,δ1=cos(ωΔt),δ2=sin(ωΔt),ω=2πf,f为实际工作频率,f通过频率跟踪方法获得;
步骤8:将与Δδa,Δδb相减,即对原相量实部和虚部进行修正,得到滤除了衰减直流分量所产生的误差的相量实部a1和虚部a2:
步骤9:令m=m+1,数据窗序号m如果大于数据窗个数M,则结束同步相量测量,否则转到步骤6继续同步相量的测量。
相较于现有技术,本发明有以下有益效果:
1、在电力系统静态条件和动态条件下都拥有较高的同步相量测量精度并具有良好的同步相量测量抗噪性。
2、滤除了电力系统中衰减直流分量带来的误差,提高了同步相量测量方法在电力系统发生故障时对同步相量测量的精度。
附图说明
图1是本发明实施例的工作流程图。
图2是τ=0.01,t≥0.06s和频率从49Hz跳变45Hz时,三种相量测量方法的幅值误差比较。
图3是τ=0.01,t≥0.06s和频率从49Hz跳变45Hz时,三种相量测量方法的相角误差比较。
图4是τ=0.01,t≥0.06s和频率从49Hz跳变45Hz时,三种相量测量方法的TVE值比较。
具体实施方式
该滤除衰减直流分量的同步相量测量方法结合图1进行说明,并对电力系统电力信号用积分公式和欧拉公式将其转换为复数形式,用连续傅里叶积分表示出复数形式中的基波分量,采用复合梯形公式和线性插值法分别计算三个连续数据窗基波分量实部和虚部的误差估计值,利用所得误差估计值对原实部和虚部进行修正。得到的效果图示于附图,具体步骤如下:
步骤1:给定电力信号x(t):
式中,X0为直流分量,X0=1;τ为衰减时间常数,τ=0.01;X(n),分别为n次谐波的幅值和初相角;其中fs为采样频率,fs=2000Hz,f0为额定频率,f0=50Hz,
步骤2:将电力信号x(t)用傅里叶级数表示:
式中,傅里叶系数
步骤3:设数据窗个数为M,M=47,每个数据窗长度为N+2=42。设m为数据窗序号,取m=1。用积分公式将电力信号x(t)傅里叶级数表达式中的傅里叶系数转换为如下离散形式,从离散形式中可获得各采样点的值:
步骤4:用欧拉公式将电力信号的傅里叶级数形式转换为复数形式:
式中,x(t)包含多次谐波,而系数n=1,2,3...;
步骤5:用连续傅里叶积分表示电力信号复数形式中的基波分量:
式中,ΔN为分数,为当实际工作频率f偏离额定频率f0时,数据窗长度发生的偏移,f通过频率跟踪方法获得;
步骤6:用复合梯形公式计算基波分量中的c1′和c1″部分:
通过线性插值法得到分数点采样点X(N+ΔN)估计值,代入c1″中
得到电力信号复数形式中的基波分量
分别计算第m、m+1、m+2三个连续数据窗中基波分量连续傅里叶积分表达式中相量的实部和虚部。
式中,Δt为采样时间间隔,ω0为额定角频率;
步骤7:设通过计算得出由于衰减直流分量而产生的相量实部误差估计值Δδa和虚部误差估计值Δδb:
式中,δ1=cos(ωΔt),δ2=sin(ωΔt),ω=2πf,f为实际工作频率,f通过频率跟踪方法获得;
步骤8:将与Δδa,Δδb相减,即对原相量实部和虚部进行修正,得到滤除了衰减直流分量所产生的误差的相量实部a1和虚部a2:
步骤9:令m=m+1,数据窗序号m如果大于数据窗个数M,则结束同步相量测量,否则转到步骤6继续同步相量的测量。
以上是本发明的较佳实施例,凡依本发明技术方案所作的改变,所产生的功能作用未超出本发明技术方案的范围时,均属于本发明的保护范围。
Claims (1)
1.一种滤除衰减直流分量的同步相量测量方法,其特征在于:
步骤1:给定电力信号x(t):
式中,X0为直流分量;τ为衰减时间常数;X(n),分别为n次谐波的幅值和初相角,其中fs为采样频率,f0为额定频率;ω=2πf,f为实际工作频率,f通过频率跟踪方法获得;
步骤2:将电力信号x(t)用傅里叶级数表示:
式中,傅里叶系数 ω0=2πf0;
步骤3:设数据窗个数为M,且每个数据窗长度为N+2;设m为数据窗序号,取m=1;用求积分公式将电力信号x(t)傅里叶级数表达式中的傅里叶系数转换为如下离散形式,从而可获得a0、an、bn的值:
式中,a0、an、bn为电力信号x(t)傅里叶级数表达式中的傅里叶系数;
步骤4:用欧拉公式将电力信号的傅里叶级数形式转换为复数形式:
式中,x(t)包含多次谐波,而系数n=1,2,3…,与步骤1所表示的n含义相同;
步骤5:用连续傅里叶积分表示电力信号复数形式中的基波分量:
式中,ΔN为分数,为当实际工作频率f偏离额定频率f0时,数据窗长度发生的偏移,f通过频率跟踪方法获得,Δt为采样时间间隔;
步骤6:用复合梯形公式和线性插值法分别计算第m、m+1、m+2三个连续数据窗中基波分量连续傅里叶积分表达式中相量的实部和虚部;
式中,Δt为采样时间间隔,ω0为额定角频率;
步骤7:设通过计算得出由于衰减直流分量而产生的相量实部误差估计值Δδa和虚部误差估计值Δδb:
式中,δ1=cos(ωΔt),δ2=sin(ωΔt),ω=2πf,f为实际工作频率,f通过频率跟踪方法获得;
步骤8:将与Δδa,Δδb相减,即对原相量实部和虚部进行修正,得到滤除了衰减直流分量所产生的误差的相量实部a1和虚部a2:
步骤9:令m=m+1,数据窗序号m如果大于数据窗个数M,则结束同步相量测量,否则转到步骤6继续同步相量的测量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410298815.2A CN104049144B (zh) | 2014-06-27 | 2014-06-27 | 一种滤除衰减直流分量的同步相量测量的实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410298815.2A CN104049144B (zh) | 2014-06-27 | 2014-06-27 | 一种滤除衰减直流分量的同步相量测量的实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104049144A CN104049144A (zh) | 2014-09-17 |
CN104049144B true CN104049144B (zh) | 2017-07-04 |
Family
ID=51502237
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410298815.2A Active CN104049144B (zh) | 2014-06-27 | 2014-06-27 | 一种滤除衰减直流分量的同步相量测量的实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104049144B (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104502707A (zh) * | 2015-01-06 | 2015-04-08 | 福州大学 | 一种基于三次样条插值的电力系统同步相量测量方法 |
CN104569689B (zh) * | 2015-01-22 | 2017-06-06 | 福州大学 | 一种基于极坐标系的电力系统同步相量测量方法 |
CN104678152B (zh) * | 2015-02-06 | 2018-01-16 | 中国铁路总公司 | 利用频域信息滤除衰减直流的电力信号同步相量测量方法 |
CN104655965B (zh) * | 2015-02-13 | 2017-07-07 | 中国农业大学 | 一种电力系统中相量测量方法 |
CN106156489B (zh) * | 2016-06-23 | 2019-01-04 | 清华大学 | 一种电力系统同步相量快速计算方法 |
CN106018960B (zh) * | 2016-07-13 | 2019-05-03 | 东北电力大学 | 一种基于压缩传感的同步相量测量方法 |
CN106154037B (zh) * | 2016-08-11 | 2019-04-02 | 中国南方电网有限责任公司 | 一种基于校验的同步相量自适应计算方法 |
CN107064630B (zh) * | 2017-03-31 | 2019-11-12 | 许继集团有限公司 | 一种电力系统频率测量方法及装置 |
CN107144734B (zh) * | 2017-05-15 | 2019-09-27 | 北京理工大学 | 一种适用于pmu的配电网高精度相量测量方法 |
CN109061345B (zh) * | 2018-07-13 | 2020-09-18 | 许继集团有限公司 | 适用于电力系统的有效值测量方法与系统 |
CN109374970B (zh) * | 2018-10-24 | 2021-08-13 | 南方电网科学研究院有限责任公司 | 实时校验的同步相量测量方法、装置、设备及存储介质 |
CN109444525B (zh) * | 2018-10-29 | 2020-10-09 | 南京国电南自维美德自动化有限公司 | 变电站测控装置交流采样数据有效值及功率因数计算方法 |
CN109521275B (zh) * | 2018-11-23 | 2020-09-25 | 南方电网科学研究院有限责任公司 | 一种同步相量确定方法、系统、装置及可读存储介质 |
CN109521273B (zh) * | 2018-11-23 | 2020-07-03 | 南方电网科学研究院有限责任公司 | 一种同步相量测量方法、系统及装置 |
CN109521274B (zh) * | 2018-11-23 | 2020-09-25 | 南方电网科学研究院有限责任公司 | 一种同步相量测量方法、系统、装置及可读存储介质 |
CN110018346A (zh) * | 2019-04-19 | 2019-07-16 | 南京工程学院 | 一种谐波状态估计的量测配置方法 |
CN110568254A (zh) * | 2019-09-27 | 2019-12-13 | 宁夏凯晨电气集团有限公司 | 一种准确检测故障电流中衰减直流分量参数的方法 |
CN111398731A (zh) * | 2020-03-09 | 2020-07-10 | 西南交通大学 | 基于多频率-泰勒模型滤除衰减直流的动态相量测量方法 |
CN111624441B (zh) * | 2020-06-10 | 2022-02-18 | 国网四川省电力公司电力科学研究院 | 一种低频振荡影响下pmu量测误差分析方法 |
CN111948454B (zh) * | 2020-08-11 | 2021-10-08 | 山东大学 | 一种抗衰减直流分量干扰的同步相量测量方法及系统 |
CN113109621B (zh) * | 2021-03-09 | 2022-05-24 | 华南理工大学 | 滤除故障信号中衰减直流分量的方法、系统、装置及介质 |
CN116577594B (zh) * | 2023-07-10 | 2023-11-14 | 广东电网有限责任公司佛山供电局 | 一种继电保护装置的故障监测方法和系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100476443C (zh) * | 2006-10-10 | 2009-04-08 | 浙江大学 | 一种同步时标可自由设定的广域电网相量同步测量装置 |
CN202153239U (zh) * | 2011-07-18 | 2012-02-29 | 武汉钢铁(集团)公司 | 差动相量检测装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003254900A (ja) * | 2002-02-28 | 2003-09-10 | Japan Science & Technology Corp | 画像同期測定方法 |
-
2014
- 2014-06-27 CN CN201410298815.2A patent/CN104049144B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100476443C (zh) * | 2006-10-10 | 2009-04-08 | 浙江大学 | 一种同步时标可自由设定的广域电网相量同步测量装置 |
CN202153239U (zh) * | 2011-07-18 | 2012-02-29 | 武汉钢铁(集团)公司 | 差动相量检测装置 |
Non-Patent Citations (2)
Title |
---|
基于同步相量测量单元的直流附加控制器研究;黄莹等;《中国电机工程学报》;20040930;第24卷(第9期);全文 * |
滤除衰减直流分量的动态同步相量测量算法;麦瑞坤等;《中国电机工程学报》;20101105;第30卷(第31期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104049144A (zh) | 2014-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104049144B (zh) | 一种滤除衰减直流分量的同步相量测量的实现方法 | |
CN103869162B (zh) | 一种基于时域准同步的动态信号相量测量方法 | |
CN102435844B (zh) | 一种频率无关的正弦信号相量计算方法 | |
CN103245832B (zh) | 基于快速s变换的谐波时频特性参数估计方法及分析仪 | |
CN101587146B (zh) | 一种谐波和间谐波参数的检测方法 | |
CN104597320B (zh) | 一种适用于多个频率交流信号计算的方法 | |
CN105548739B (zh) | 一种避雷器运行状态信号处理方法 | |
CN102323477A (zh) | 一种电网间谐波测量方法及测量仪 | |
Jain et al. | An adaptive time-efficient technique for harmonic estimation of nonstationary signals | |
Chen et al. | Improved VMD-FRFT based on initial center frequency for early fault diagnosis of rolling element bearing | |
CN104020352A (zh) | 一种适用于m类pmu单元的同步相量测量方法 | |
CN106597229A (zh) | 35kV以上变电设备绝缘在线监测系统的工作方法 | |
CN104678152A (zh) | 利用频域信息滤除衰减直流的电力信号同步相量测量方法 | |
CN106324340B (zh) | 一种同步相量和频率测量动态性能的方法 | |
Warichet et al. | Considerations about synchrophasors measurement in dynamic system conditions | |
CN106501602A (zh) | 一种基于滑窗频谱分离的基波参数测量方法 | |
CN104777356A (zh) | 一种基于神经网络的实时高精度谐波检测方法 | |
CN104502707A (zh) | 一种基于三次样条插值的电力系统同步相量测量方法 | |
CN108333426A (zh) | 基于傅氏算法的电力系统频率测量方法 | |
CN109061345B (zh) | 适用于电力系统的有效值测量方法与系统 | |
CN104502703A (zh) | 基于频域动态模型的电力信号同步相量测量方法 | |
CN103543331B (zh) | 一种计算电信号谐波和间谐波的方法 | |
CN104360156A (zh) | 一种电力系统频率偏移时信号相角的测量方法 | |
CN202102051U (zh) | 一种电网间谐波测量仪 | |
Park et al. | Advanced frequency estimation technique using gain compensation |
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 |