CN112881799A - 基于fft变换地电力***谐波检测方法 - Google Patents

基于fft变换地电力***谐波检测方法 Download PDF

Info

Publication number
CN112881799A
CN112881799A CN202110061016.3A CN202110061016A CN112881799A CN 112881799 A CN112881799 A CN 112881799A CN 202110061016 A CN202110061016 A CN 202110061016A CN 112881799 A CN112881799 A CN 112881799A
Authority
CN
China
Prior art keywords
frequency
harmonic
phase
power system
formula
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.)
Pending
Application number
CN202110061016.3A
Other languages
English (en)
Inventor
鄢烈奇
唐登平
明东岳
李俊
姚鹏
王尚鹏
夏天
李君�
庞博
王信
王琪
汪应春
余鹤
雷鸣
马奔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Metering Center of State Grid Hubei Electric Power Co Ltd
Original Assignee
Wuhan Glory Road Intelligent Technology Co ltd
Metering Center of State Grid Hubei Electric Power Co Ltd
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 Wuhan Glory Road Intelligent Technology Co ltd, Metering Center of State Grid Hubei Electric Power Co Ltd filed Critical Wuhan Glory Road Intelligent Technology Co ltd
Priority to CN202110061016.3A priority Critical patent/CN112881799A/zh
Publication of CN112881799A publication Critical patent/CN112881799A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum analysis; Fourier analysis

Landscapes

  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明提出了基于FFT变换地电力***谐波检测方法,将频率变化率设置为前次和当次测得的基波频率进行差商运算的结果,用本实施例的频率变化率去修正当次测量的归一化频率校正量,并通过频率校正公式、幅值校正公式和相位校正公式计算修正后的频率、幅值、相位,最终得到的频率、幅值、相位更符合实际结果,进而可以适应频率动态变化的电网环境;将前次测得的频率用于本次频率变化率的计算,从而使计算时间缩短了一半。

Description

基于FFT变换地电力***谐波检测方法
技术领域
本发明涉及谐波检测技术领域,尤其涉及基于FFT变换地电力***谐波检测方法。
背景技术
谐波不但使电网效率降低,还会形成谐振环流,严重时候甚至会烧坏电机。常用傅里叶分析法抑制谐波。FFT分析法将检测到地一个周期地谐波信号通过快速FFT进行分解,得到各次谐波地幅值和相位,将拟消除地谐波分量通过带通滤波器或傅里叶变换器得出所需地误差信号,再将该误差信号进行傅里叶反变换,即可得补偿信号。这种方法不能分开有功电流和无功电流,并且要求严格地同步采样和等间隔采样,否则将会出现栅栏效应和频谱泄露现象,从而影响检测地精确度。
非同步采样下的测量方法从采样方式上可以分为定速率采样方法和自适应采样方法,而现有的定速率采样方法不能根据电网频率的变化实时调整采样频率以实现同步采样,并且当电网频率发生动态变化时,其测量精度急剧降低,因此,为了解决上述问题,本发明提供了基于FFT变换地电力***谐波检测方法,通过相位差校正法对两段加窗地时域序列做FFT并利用对应峰值谱线地相位差来进行频率、幅值和相位的校正,以得到同步数据,并且可以适应频率动态变化地电网环境。
发明内容
有鉴于此,本发明提出了基于FFT变换地电力***谐波检测方法,通过相位差校正法对两段加窗地时域序列做FFT并利用对应峰值谱线地相位差来进行频率、幅值和相位的校正,以得到同步数据,并且可以适应频率动态变化的电网环境。
本发明的技术方案是这样实现的:本发明提供了基于FFT变换地电力***谐波检测方法,包括以下步骤:
S1、对采集的非同步数据进行FFT变换,得到频谱序列,由频谱序列推导归一化频率校正量;
S2、基于频率变化率修正归一化频率校正量,将修正后的归一化频率校正量代入频率校正公式、幅值校正公式和相位校正公式中得出K次谐波的频率、幅值和相位。
在以上技术方案的基础上,优选的,S1具体包括以下步骤:
S101、确定谐波信号的相位,并将该谐波信号在时域上向左平移时间长度为t0的长度,确定平移后的相位;
S102、将平移前与平移后的相位相减得到两段信号的相位差,根据相位差推得归一化频率校正量。
在以上技术方案的基础上,优选的,S102中平移前与平移后的相位相减得到两段信号的相位差为:
Δφk=2πfkt0
式中,Δφk为两段信号的相位差;fk为k次谐波的频率;t0为第二段信号平移的时间长度。
进一步优选的,S2中频率变化率为当次测得的基波频率减去前次测得的基波频率再除以两次测量的时间间隔的结果。
进一步优选的,S2中修正后的归一化频率校正量为:
Figure BDA0002902353180000031
式中,Δmk为修正后的归一化频率校正量;N为各段信号的采样点数;L为第二段信号平移的点数;fk为k次谐波的频率;fs为采样率;mk为归一化峰值谱线频率。
进一步优选的,频率校正公式为:
Figure BDA0002902353180000032
进一步优选的,幅值校正公式为:
Figure BDA0002902353180000033
式中,Ak为k次谐波的幅值;Amk为谐波信号频谱中峰值谱线的幅值;Wnorm(Δmk)为Δmk的傅里叶变换函数。
进一步优选的,相位校正公式为:
Figure BDA0002902353180000034
式中,
Figure BDA0002902353180000035
为k次谐波的相位;Ik和Rk分别为归一化频率为mk时Wnorm(Δmk)函数的虚部和实部。
本发明的基于FFT变换地电力***谐波检测方法相对于现有技术具有以下有益效果:
(1)将频率变化率设置为前次和当次测得的基波频率进行差商运算的结果,用本实施例的频率变化率去修正当次测量的归一化频率校正量,并通过频率校正公式、幅值校正公式和相位校正公式计算修正后的频率、幅值、相位,最终得到的频率、幅值、相位更符合实际结果,进而可以适应频率动态变化的电网环境;
(2)将前次测得的频率用于本次频率变化率的计算,从而使计算时间缩短了一半。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明基于FFT变换地电力***谐波检测方法的流程图;
图2为本发明基于FFT变换地电力***谐波检测方法中步骤S1的具体流程图。
具体实施方式
下面将结合本发明实施方式,对本发明实施方式中的技术方案进行清楚、完整地描述,显然,所描述的实施方式仅仅是本发明一部分实施方式,而不是全部的实施方式。基于本发明中的实施方式,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式,都属于本发明保护的范围。
对于定速率采样,频域上的方法需要推导出相应的频率、幅值、相位校正公式,各校正公式均与归一化频率校正量有关,但现有方法在推导归一化频率校正量计算公式过程中均未考虑到电网频率变化率不为零时的情况,因而在频率动态过程中测量精度严重降低,为了解决上述问题,如图1所示,本实施例提供了基于FFT变换地电力***谐波检测方法,包括以下步骤:
S1、对采集的非同步数据进行FFT变换,得到频谱序列,由频谱序列推导归一化频率校正量;如图2所示,具体包括以下步骤:
S101、确定谐波信号的相位,并将该谐波信号在时域上向左平移时间长度为t0的长度,确定平移后的相位;
S102、将平移前与平移后的相位相减得到两段信号的相位差,根据相位差推得归一化频率校正量。其中,平移前与平移后的相位相减得到两段信号的相位差为:Δφk=2πfkt0;式中,Δφk为两段信号的相位差;fk为k次谐波的频率;t0为第二段信号平移的时间长度。根据相位差推得归一化频率校正量可以采用现有技术实现,在此不再累述。本实施例是对计算出的归一化频率校正量进一步校正,使其适用于频率动态过程。
S2、基于频率变化率修正归一化频率校正量,将修正后的归一化频率校正量代入频率校正公式、幅值校正公式和相位校正公式中得出K次谐波的频率、幅值和相位。
其中,频率变化率为当次测得的基波频率减去前次测得的基波频率再除以两次测量的时间间隔的结果。由于在电力谐波测量技术领域中,频率变化率是一个很重要的参数,因此,本实施例中,频率变化率为前次和当次测得的基波频率进行差商运算的结果。用本实施例的频率变化率去修正当次测量的归一化频率校正量,并通过频率校正公式、幅值校正公式和相位校正公式计算修正后的频率、幅值、相位。本实施例中,将前次测得的频率用于本次频率变化率的计算,从而使计算时间缩短了一半。
进一步优选的,修正后的归一化频率校正量为:
Figure BDA0002902353180000061
式中,Δmk为修正后的归一化频率校正量;N为各段信号的采样点数;L为第二段信号平移的点数;fk为k次谐波的频率;fs为采样率;mk为归一化峰值谱线频率。
进一步优选的,频率校正公式为:
Figure BDA0002902353180000062
进一步优选的,幅值校正公式为:
Figure BDA0002902353180000063
式中,Ak为k次谐波的幅值;Amk为谐波信号频谱中峰值谱线的幅值;Wnorm(Δmk)为Δmk的傅里叶变换函数。
进一步优选的,相位校正公式为:
Figure BDA0002902353180000064
式中,
Figure BDA0002902353180000065
为k次谐波的相位;Ik和Rk分别为归一化频率为mk时Wnorm(Δmk)函数的虚部和实部。
本实施例的有益效果为:将频率变化率设置为前次和当次测得的基波频率进行差商运算的结果,用本实施例的频率变化率去修正当次测量的归一化频率校正量,并通过频率校正公式、幅值校正公式和相位校正公式计算修正后的频率、幅值、相位,最终得到的频率、幅值、相位更符合实际结果,进而可以适应频率动态变化的电网环境;
本实施例中,将前次测得的频率用于本次频率变化率的计算,从而使计算时间缩短了一半。
以上所述仅为本发明的较佳实施方式而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.基于FFT变换地电力***谐波检测方法,其特征在于:包括以下步骤:
S1、对采集的非同步数据进行FFT变换,得到频谱序列,由频谱序列推导归一化频率校正量;
S2、基于频率变化率修正归一化频率校正量,将修正后的归一化频率校正量代入频率校正公式、幅值校正公式和相位校正公式中得出K次谐波的频率、幅值和相位。
2.如权利要求1所述的基于FFT变换地电力***谐波检测方法,其特征在于:所述S1具体包括以下步骤:
S101、确定谐波信号的相位,并将该谐波信号在时域上向左平移时间长度为t0的长度,确定平移后的相位;
S102、将平移前与平移后的相位相减得到两段信号的相位差,根据相位差推得归一化频率校正量。
3.如权利要求2所述的基于FFT变换地电力***谐波检测方法,其特征在于:所述S102中平移前与平移后的相位相减得到两段信号的相位差为:
Δφk=2πfkt0
式中,Δφk为两段信号的相位差;fk为k次谐波的频率;t0为第二段信号平移的时间长度。
4.如权利要求1所述的基于FFT变换地电力***谐波检测方法,其特征在于:所述S2中频率变化率为当次测得的基波频率减去前次测得的基波频率再除以两次测量的时间间隔的结果。
5.如权利要求4所述的基于FFT变换地电力***谐波检测方法,其特征在于:所述S2中修正后的归一化频率校正量为:
Figure FDA0002902353170000021
式中,Δmk为修正后的归一化频率校正量;N为各段信号的采样点数;L为第二段信号平移的点数;fk为k次谐波的频率;fs为采样率;mk为归一化峰值谱线频率。
6.如权利要求5所述的基于FFT变换地电力***谐波检测方法,其特征在于:所述频率校正公式为:
Figure FDA0002902353170000022
7.如权利要求5所述的基于FFT变换地电力***谐波检测方法,其特征在于:所述幅值校正公式为:
Figure FDA0002902353170000023
式中,Ak为k次谐波的幅值;Amk为谐波信号频谱中峰值谱线的幅值;Wnorm(Δmk)为Δmk的傅里叶变换函数。
8.如权利要求5所述的基于FFT变换地电力***谐波检测方法,其特征在于:所述相位校正公式为:
Figure FDA0002902353170000024
式中,
Figure FDA0002902353170000025
为k次谐波的相位;Ik和Rk分别为归一化频率为mk时Wnorm(Δmk)函数的虚部和实部。
CN202110061016.3A 2021-01-18 2021-01-18 基于fft变换地电力***谐波检测方法 Pending CN112881799A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110061016.3A CN112881799A (zh) 2021-01-18 2021-01-18 基于fft变换地电力***谐波检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110061016.3A CN112881799A (zh) 2021-01-18 2021-01-18 基于fft变换地电力***谐波检测方法

Publications (1)

Publication Number Publication Date
CN112881799A true CN112881799A (zh) 2021-06-01

Family

ID=76048800

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110061016.3A Pending CN112881799A (zh) 2021-01-18 2021-01-18 基于fft变换地电力***谐波检测方法

Country Status (1)

Country Link
CN (1) CN112881799A (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106918741A (zh) * 2017-03-02 2017-07-04 浙江大学 应用于频率宽范围波动电网的自适应采样的相位差校正法
CN106970264A (zh) * 2017-03-02 2017-07-21 浙江大学 一种考虑电网频率变化率的改进相位差校正法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106918741A (zh) * 2017-03-02 2017-07-04 浙江大学 应用于频率宽范围波动电网的自适应采样的相位差校正法
CN106970264A (zh) * 2017-03-02 2017-07-21 浙江大学 一种考虑电网频率变化率的改进相位差校正法

Similar Documents

Publication Publication Date Title
CN105738696B (zh) 全相位时移相位差频率估计方法及装置
CN110837001B (zh) 一种电力***中谐波和间谐波的分析方法与装置
CN101261292A (zh) 基于5项Rife-Vincent(I)窗双谱线插值FFT的基波与谐波检测方法
Wen et al. Hanning self-convolution window and its application to harmonic analysis
CN110954746A (zh) 一种基于四项Nuttall余弦窗的六插值FFT算法
CN109521275B (zh) 一种同步相量确定方法、***、装置及可读存储介质
CN105137180B (zh) 基于六项余弦窗四谱线插值的高精度谐波分析方法
CN102818930B (zh) 一种高精度快速计算电力谐波参数的方法
CN108896944B (zh) 一种同步测量装置实验室校准仪及其同步相量测量方法
CN105353215A (zh) 基于Nuttall窗四谱线插值FFT的谐波检测方法
CN108535613A (zh) 一种基于组合窗函数的电压闪变参数检测方法
CN106405229A (zh) 一种基波和谐波电能计量方法
CN102495285B (zh) 对称窗函数功率重心估计电力谐波参数的方法
CN114460527B (zh) 关联度延拓Hilbert移相电子式互感器校验仪溯源方法和***
CN109444539B (zh) 一种基于克拉克变换的同步相量测量方法
CN104020350B (zh) 一种克服频率摄动的电压基波分量检测方法
CN112881799A (zh) 基于fft变换地电力***谐波检测方法
CN106053936B (zh) 一种获取电学信号瞬时频率的方法及***
CN117234176A (zh) 一种基于扫频试验的控制回路频域辨识方法
CN106324342A (zh) 一种基于查表的谐波检测方法
CN114487589A (zh) 电网宽频信号自适应测量方法、装置及***
CN111579868B (zh) 一种高次谐波的测量方法及装置
CN115372696A (zh) 基于递归傅里叶级数的电网电压频率计算方法和装置
CN113129912A (zh) 一种单音信号的检测方法
CN108982966A (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
TA01 Transfer of patent application right

Effective date of registration: 20220921

Address after: No. 488, QingWang Road, Donghu eco tourism scenic spot, Wuhan City, Hubei Province, 430000

Applicant after: Marketing Service Center (metering center) of State Grid Hubei Electric Power Co.,Ltd.

Address before: No. 488, QingWang Road, Donghu eco tourism scenic spot, Wuhan City, Hubei Province, 430000

Applicant before: Marketing Service Center (metering center) of State Grid Hubei Electric Power Co.,Ltd.

Applicant before: WUHAN GLORY ROAD INTELLIGENT TECHNOLOGY Co.,Ltd.

TA01 Transfer of patent application right
RJ01 Rejection of invention patent application after publication

Application publication date: 20210601

RJ01 Rejection of invention patent application after publication