CN108008187A - 基于变分模态分解的电网谐波检测方法 - Google Patents
基于变分模态分解的电网谐波检测方法 Download PDFInfo
- Publication number
- CN108008187A CN108008187A CN201711297045.XA CN201711297045A CN108008187A CN 108008187 A CN108008187 A CN 108008187A CN 201711297045 A CN201711297045 A CN 201711297045A CN 108008187 A CN108008187 A CN 108008187A
- Authority
- CN
- China
- Prior art keywords
- frequency
- mode
- harmonic wave
- instantaneous frequency
- mode decomposition
- 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
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 36
- 238000001514 detection method Methods 0.000 title claims abstract description 14
- 238000001228 spectrum Methods 0.000 claims abstract description 12
- 238000012544 monitoring process Methods 0.000 abstract description 2
- 238000000034 method Methods 0.000 description 12
- 238000012880 independent component analysis Methods 0.000 description 6
- 239000000284 extract Substances 0.000 description 3
- 238000011217 control strategy Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000004087 circulation Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开一种基于变分模态分解的电网谐波检测方法,是以快速傅里叶变换选择变分模态分解关键参数K,并以频谱中分布的各谐波尖峰处的频率值为模态分解的中心频率初始值,可以有效限制VMD局部最优方向,使最终结果局部最优收敛于目标频率处。同时通过监测模态瞬时频率均值与目标模态瞬时频率之差以及模态瞬时频率方差,逐步增大惩罚因子,既可以得到更加精确的模态分量,又可以有效避免模态混叠、冗余和丢失现象。
Description
技术领域
本发明涉及一种电网谐波检测方法,尤其是一种可提取谐波信号的相位、幅值及频率信息的基于变分模态分解的电网谐波检测方法。
背景技术
随着新能源发电***的大规模应用,电力电子设备已越来越多地被引入到电力***中,从而产生过多电网谐波。电网谐波会导致电力电子设备过流,严重时将造成新能源发电***大规模脱网等重大事故。理论上讲,电力电子设备可以在明确各类电网谐波信号特征全面信息(幅值、相位及频率)的前提下调整控制策略抑制过流,避免发生因保护动作而导致设备脱网的现象。
目前,电网谐波检测方法的研究主要集中在傅里叶变换(Fourier Transform,FT)、小波变换(Wavelet Transform, WT)、经验模态分解(Empirical ModesDecomposition, EMD)以及独立成分分析(Independent Component Analysis, ICA)4种方法。FT方法可将谐波信号特征的频率信息有效提取出来,但丢失了相位、幅值信息,无法得知扰动的发生时间、持续时长以及严重程度。WT方法可以同时获得谐波信号的相位与频率信息,从而得到比FT方法更好的辨识结果。但却具有对噪声敏感、鲁棒性差、频谱泄露等缺点。EMD方法具有无需待测信号的先验信息,可对非平稳信号进行分析等优点,但其缺乏数学理论支撑且所得结果易发生模态混叠现象,使不同的检测结果含有相同的频率信息,不能够对谐波信号进行有效分离。ICA方法能够在源信号未知的情况下,提取电网谐波信号特征信息,是一种功能强大的信号处理方法。但该方法要求输入为多路信号,而电网信号多数为单通道电压信号,因此,采用ICA方法进行电网谐波检测需要先将单通道信号转换为多通道信号,增加了算法复杂度;同时由于ICA的算法原理忽略信号的幅值信息,从而无法得知谐波扰动的严重程度。综上所述,现有方法并不能够将谐波信号的相位、幅值及频率信息全部提取出来。
变分模态分解(Variational Mode Decomposition,VMD)方法由K.Dragomiretskiy和D. Zosso于2014年提出,是一种非迭代的新的模态估计方法,能够自适应地实现信号的频域划分及各分量的有效分离,并且所得结果同时含有幅值、相位及频率信息,已被应用到多个研究领域。
变分模态分解有两个重要参数,即模态函数个数K和惩罚因子。模态函数个数K取值过小时,所得到的模态分量出现了模态混叠现象,不能够得到正确的模态分量;而模态函数个数K取值过大,则会有多余模态分量的产生,导致模态冗余,同时也增大了算法运算时间,降低算法效率。同样,当取值过小时,造成了模态混叠现象,而增大参数,虽可有效增加模态及其中心频率分离精度,但取值过大将出现模态混叠和模态丢失现象。另外,现有变分模态分解中谐波分量的中心频率初始值为零,有时最终结果局部会偏离目标频率。因此,迄今为止未见既可增加模态及其中心频率分离精度,又可避免模态混叠、冗余及丢失的采用VMD进行电网谐波信号特征检测的相关报道。
发明内容
本发明是为了解决现有技术所存在的上述技术问题,提供一种可提取谐波信号的相位、幅值及频率信息的基于变分模态分解的电网谐波检测方法。
本发明的技术解决方案是:一种基于变分模态分解的电网谐波检测方法,其特征在于依次按照如下步骤进行:
a. 取电网电压信号;
b. 对电网电压信号进行FFT变换;
c. 以频谱中分布的谐波个数为模态分解个数K,以频谱中分布的各谐波尖峰处的频率值为模态分解的中心频率初始值;
d. 确定惩罚因子初始值;
e. 进行VMD分解;
f. 计算模态瞬时频率fins;
g. 判断是否满足;
式中:为模态瞬时频率值,为目标模态瞬时频率,为模态瞬时频率均值,所述小于或等于0.01;
否,按步长2000增大惩罚因子,返回e步骤;
是,则进行下一步骤;
h. 计算模态瞬时频率方差;
i. 判断是否满足 ,所述小于或等于2;
否,则按步长2000、1000、500或100继续增大惩罚因子,返回e步骤;
是,则获得模态分量。
本发明是以快速傅里叶变换选择变分模态分解关键参数K,并以频谱中分布的各谐波尖峰处的频率值为模态分解的中心频率初始值,可以有效限制VMD局部最优方向,使最终结果局部最优收敛于目标频率处。同时通过监测模态瞬时频率均值与目标模态瞬时频率之差以及模态瞬时频率方差,逐步增大惩罚因子,既可以得到更加精确的模态分量,又可以有效避免模态混叠、冗余和丢失现象。通过仿真和真实电压信号进行实验,能够将电网中谐波、间谐波以及基波分量准确提取出来,从而为控制器及时调整控制策略抑制过流,避免电力设备因发生保护动作而导致脱网的现象发生。
附图说明
图1是本发明实施例的流程图。
图2是本发明实施例1的仿真信号时域波形图。
图3是本发明实施例1的仿真信号频域波形图。
图4是本发明实施例1第一次分离所得各个分量的时域波形图。
图5是本发明实施例1第一次分离所得各个分量的频域波形图。
图6是本发明实施例1最后一次分离所得各个分量的时域波形图。
图7是本发明实施例1最后一次分离所得各个分量的频域波形图。
图8是本发明实施例2的仿真信号时域波形图。
图9是本发明实施例2的仿真信号频域波形图。
图10是本发明实施例2第一次分离所得各个分量的时域波形图。
图11是本发明实施例2第一次分离所得各个分量的频域波形图。
图12是本发明实施例4的电网电压时域波形图。
图13是本发明实施例4的电网电压频域波形图。
图14是本发明实施例3第一次分离所得各个分量的时域波形图。
图15是本发明实施例3第一次分离所得各个分量的频域波形图。
具体实施方式
实施例1:
本发明的一种基于变分模态分解的电网谐波检测方法如图1所示,依次按照如下步骤进行:
a. 取电网电压信号;
基于MATLAB平台建立电网电压仿真信号,由3个不同频率分量以及高斯白噪声混合而成,其表达式为:
其中,3个分量的频率分别为12Hz,50Hz,150Hz, 是高斯白噪声。N是噪声幅值系数,N=0.2,SNR=7.09dB(只考虑150Hz分量与噪声的比较),
仿真信号时域波形如图2所示;
b. 对电网电压信号进行FFT变换,频域波形如图3所示,图3中Y轴为幅值,X轴为角频率。从图3可见,该仿真信号主要由12Hz,50Hz,150Hz的3个谐波分量组成;
c. 以频谱中分布的谐波个数为模态分解个数K,即K=3;以频谱中分布的各谐波尖峰处的频率值为模态分解的中心频率初始值,即中心频率初始值分别为12Hz,50Hz,150Hz;
d. 确定惩罚因子初始值,令;
e. 进行VMD分解,结果如图4、图5所示。
图4是本发明实施例1第一次分离所得各个分量的时域波形图。
图5是本发明实施例1第一次分离所得各个分量的频域波形图。
从图4、图5可以看出,没有模态冗余和丢失问题,但是150Hz模态分量的瞬时频率仍有较大波动;
f. 计算模态瞬时频率fins;
g. 判断是否满足;
式中:为模态瞬时频率值,为目标模态瞬时频率,为模态瞬时频率均值,所述等于0.01;
否,按步长2000增大惩罚因子,返回e步骤;
是,说明通过VMD方法所得频率为目标频率,则进行下一步骤;
h. 计算模态瞬时频率方差, 3个模态分量的瞬时频率方差分别为0.013,1.17和12.36;
i. 判断是否满足 ,所述等于0.5;
结果并不满足小于0.5的要求,因此,则按步长2000继续增大惩罚因子,此时当波动量逐渐接近阈值时,可减小步长,如取步长为1000,500或100,返回e步骤;
本发明实施例1经过多个循环,直至179800,返回e步骤,即进行VMD分解,……计算模态瞬时频率方差, 3个模态分量的瞬时频率方差分别为0.0022,0.11和0.49,结果满足了小于0.5的要求,说明瞬时频率足够平稳,则获得模态分量。
获得的模态分量如图6、图7所示。
图6是本发明实施例1最后一次分离所得各个分量的时域波形图。
图7是本发明实施例1最后一次分离所得各个分量的频域波形图。
从图6、图7可以看出,模态分量的瞬时频率波动大幅度减小,即获得精确的模态分量,又有效避免模态混叠、冗余和丢失现象。
实施例2:
本发明的一种基于变分模态分解和去除噪声的电网谐波检测方法如图1所示,依次按照如下步骤进行:
a. 取电网电压信号;
基于MATLAB平台建立电网电压仿真信号,由3个不同频率分量以及高斯白噪声混合而成,其表达式为:
其中,3个分量的频率分别为12Hz,50Hz,150Hz, 是高斯白噪声。N是噪声幅值系数,N=0.8,SNR= -18.96dB(只考虑150Hz分量与噪声的比较),
仿真信号时域波形如图8所示;
b. 对电网电压信号进行FFT变换,频域波形如图9所示,图9中Y轴为幅值,X轴为角频率。从图9可见,该仿真信号主要由12Hz,50Hz,150Hz的3个谐波分量组成;
c. 以频谱中分布的谐波个数为模态分解个数K,即K=3;以频谱中分布的各谐波尖峰处的频率值为模态分解的中心频率初始值,即中心频率初始值分别为12Hz,50Hz,150Hz;
d. 选择惩罚因子,2000000;
e. 进行VMD分解;
f. 计算模态瞬时频率fins;
g. 判断是否满足;
式中:为模态瞬时频率值,为目标模态瞬时频率,为模态瞬时频率均值,所述等于0.01;
否,按步长2000增大惩罚因子,返回e步骤;
是,则进行下一步骤;
h. 计算模态瞬时频率方差, 3个模态分量的瞬时频率方差分别为其方差分别为0.0068,1.45以及0.12;
i. 判断是否满足 ,所述等于1.5;
结果满足了小于1.5的要求,则获得模态分量。
获得的模态分量如图10、图11所示。
图10是本发明实施例2第一次分离所得各个分量的时域波形图。
图11是本发明实施例2第一次分离所得各个分量的频域波形图。
从图10、图11可以看出,所得模态分量的瞬时频率基本稳定,即获得精确的模态分量,又有效避免模态混叠、冗余和丢失现象。
实施例3:
本发明的一种基于变分模态分解电网谐波检测方法如图1所示,其特征在于依次按照如下步骤进行:
a. 在风电场采集电网电压信号,时域波形如图12所示;
b. 对电网电压信号进行FFT变换,所得频域波形如图13所示,从图13中可以看出,该电压信号主要含有50Hz基波与10Hz间谐波分量;
c. 以频谱中分布的谐波个数为模态分解个数K,即K=2;以频谱中分布的各谐波尖峰处的频率值为模态分解的中心频率初始值,即中心频率初始值分别为050Hz,10Hz;
d. 选择惩罚因子,令;
e. 进行VMD分解;
f. 计算模态瞬时频率fins;
g. 判断是否满足;
式中:为模态瞬时频率值,为目标模态瞬时频率,为模态瞬时频率均值,所述等于0.01;
否,按步长2000增大惩罚因子,返回e步骤;
是,则进行下一步骤;
h. 计算模态瞬时频率方差, 2个模态分量的瞬时频率方差分别为0.000174,0.081;
i. 判断是否满足 ,所述等于0.1;
结果满足了小于0.1的要求,则获得模态分量。
获得的模态分量如图14、图15所示。
图14是本发明实施例3第一次分离所得各个分量的时域波形图。
图15是本发明实施例3第一次分离所得各个分量的频域波形图。
从图14、图15可以看出,所得到的有效模态分量时域波形图均比较平稳,其中10Hz间谐波分量的幅值为900V,50Hz基波分量幅值为11090V。可将电网电压信号中有效模态分量提出,模态分量的中心频率与原始信号含有谐波分量频率相同。即获得精确的模态分量,又有效避免模态混叠、冗余和丢失现象。
结论:通过仿真和真实电压信号进行实验,有效的验证了本发明的正确性与有效性。
Claims (1)
1.一种基于变分模态分解的电网谐波检测方法,其特征在于依次按照如下步骤进行:
a. 取电网电压信号;
b. 对电网电压信号进行FFT变换;
c. 以频谱中分布的谐波个数为模态分解个数K,以频谱中分布的各谐波尖峰处的频率值为模态分解的中心频率初始值;
d. 确定惩罚因子初始值;
e. 进行VMD分解;
f. 计算模态瞬时频率fins;
g. 判断是否满足;
式中:为模态瞬时频率值,为目标模态瞬时频率,为模态瞬时频率均值,所述小于或等于0.01;
否,按步长2000增大惩罚因子,返回e步骤;
是,则进行下一步骤;
h. 计算模态瞬时频率方差;
i. 判断是否满足 ,所述小于或等于2;
否,则按步长2000、1000、500或100继续增大惩罚因子,返回e步骤;
是,则获得模态分量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711297045.XA CN108008187B (zh) | 2017-12-08 | 2017-12-08 | 基于变分模态分解的电网谐波检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711297045.XA CN108008187B (zh) | 2017-12-08 | 2017-12-08 | 基于变分模态分解的电网谐波检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108008187A true CN108008187A (zh) | 2018-05-08 |
CN108008187B CN108008187B (zh) | 2020-01-10 |
Family
ID=62057805
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711297045.XA Expired - Fee Related CN108008187B (zh) | 2017-12-08 | 2017-12-08 | 基于变分模态分解的电网谐波检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108008187B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109542089A (zh) * | 2018-12-21 | 2019-03-29 | 浙江大学 | 一种基于改进变分模态分解的工业过程非线性振荡检测方法 |
CN110726875A (zh) * | 2019-12-02 | 2020-01-24 | 山东大学 | 一种新能源柔性直流并网暂态谐波检测方法及*** |
CN111081274A (zh) * | 2020-01-03 | 2020-04-28 | 西安电子科技大学 | 基于瞬时频率的带静默间断信号的端点识别方法及应用 |
CN112903296A (zh) * | 2021-01-25 | 2021-06-04 | 燕山大学 | 一种滚动轴承故障检测方法及*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106017926A (zh) * | 2016-05-13 | 2016-10-12 | 山东理工大学 | 基于变模态分解的滚动轴承故障诊断方法 |
CN106019271A (zh) * | 2016-04-27 | 2016-10-12 | 南京理工大学 | 一种基于变分模态分解的多人穿墙时变呼吸信号检测方法 |
JP2017173345A (ja) * | 2012-09-27 | 2017-09-28 | 株式会社ダイヘン | 周波数解析装置、当該周波数解析装置を用いた信号処理装置、および、当該信号処理装置を用いた高周波測定装置 |
-
2017
- 2017-12-08 CN CN201711297045.XA patent/CN108008187B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017173345A (ja) * | 2012-09-27 | 2017-09-28 | 株式会社ダイヘン | 周波数解析装置、当該周波数解析装置を用いた信号処理装置、および、当該信号処理装置を用いた高周波測定装置 |
CN106019271A (zh) * | 2016-04-27 | 2016-10-12 | 南京理工大学 | 一种基于变分模态分解的多人穿墙时变呼吸信号检测方法 |
CN106017926A (zh) * | 2016-05-13 | 2016-10-12 | 山东理工大学 | 基于变模态分解的滚动轴承故障诊断方法 |
Non-Patent Citations (2)
Title |
---|
K.P. SOMAN ET AL.: "Recursive Variational Mode Decomposition Algorithm for Real Time Power Signal Decomposition", 《PROCEDIA TECHNOLOGY》 * |
马增强等: "基于变分模态分解和Teager 能量算子的滚动轴承故障特征提取", 《振动与冲击》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109542089A (zh) * | 2018-12-21 | 2019-03-29 | 浙江大学 | 一种基于改进变分模态分解的工业过程非线性振荡检测方法 |
CN110726875A (zh) * | 2019-12-02 | 2020-01-24 | 山东大学 | 一种新能源柔性直流并网暂态谐波检测方法及*** |
CN111081274A (zh) * | 2020-01-03 | 2020-04-28 | 西安电子科技大学 | 基于瞬时频率的带静默间断信号的端点识别方法及应用 |
CN111081274B (zh) * | 2020-01-03 | 2022-03-25 | 西安电子科技大学 | 基于瞬时频率的带静默间断信号的端点识别方法及应用 |
CN112903296A (zh) * | 2021-01-25 | 2021-06-04 | 燕山大学 | 一种滚动轴承故障检测方法及*** |
CN112903296B (zh) * | 2021-01-25 | 2021-12-14 | 燕山大学 | 一种滚动轴承故障检测方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN108008187B (zh) | 2020-01-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108008187A (zh) | 基于变分模态分解的电网谐波检测方法 | |
Achlerkar et al. | Variational mode decomposition and decision tree based detection and classification of power quality disturbances in grid-connected distributed generation system | |
Song et al. | A fault-location method for VSC-HVDC transmission lines based on natural frequency of current | |
CN108390395B (zh) | 一种在线自适应抑制柔性直流引起的高频振荡方法 | |
WO2019140584A1 (zh) | 一种在线自适应抑制柔性直流引起的高频振荡方法 | |
CN104459398B (zh) | 一种采用二维形态学降噪的电能质量复合扰动识别方法 | |
CN108535613B (zh) | 一种基于组合窗函数的电压闪变参数检测方法 | |
CN107991531A (zh) | 基于变分模态分解和去除噪声的电网谐波检测方法 | |
Prince et al. | Targeted approach to apply masking signal-based empirical mode decomposition for mode identification from dynamic power system wide area measurement signal data | |
Samantaray et al. | Power transformer protection using S-transform with complex window and pattern recognition approach | |
CN110221116A (zh) | 基于加窗插值和解析模式分解的电压闪变包络检测方法 | |
CN104320157A (zh) | 一种电力线双向工频通信上行信号检测方法 | |
Agrawal et al. | Harmonics and inter harmonics estimation of DFIG based standalone wind power system by parametric techniques | |
CN103647550A (zh) | 一种用于动态电压无功补偿的锁相环方法 | |
CN105044458B (zh) | 一种基于tt变换的电力***谐波检测方法 | |
CN102375085B (zh) | 电压骤升或骤降监测方法及使用该方法的监测装置 | |
CN102759659B (zh) | 一种电力***中电气信号谐波瞬时值的提取方法 | |
Becejac et al. | Impact of PMU data errors on modal extraction using matrix pencil method | |
CN108152679A (zh) | 一种基于hht与闪变功率的闪变源定位方法 | |
Li et al. | Power quality disturbance detection based on IEWT | |
CN110763914B (zh) | 一种pwm谐波有效值计算方法 | |
Reza et al. | Frequency adaptive instantaneous power quality analysis using frequency locked loop based Kalman filter technique | |
Katić et al. | Wavelet transform for voltage dips detection in a microgrid with distributed generation | |
CN105044457B (zh) | 一种抗噪声的电力***谐波检测方法 | |
Bukka et al. | Performance Analysis of three phase shunt hybrid active power filter |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200110 |
|
CF01 | Termination of patent right due to non-payment of annual fee |