基于变分模态分解的电网谐波检测方法与流程

本发明涉及一种电网谐波检测方法，尤其是一种可提取谐波信号的相位、幅值及频率信息的基于变分模态分解的电网谐波检测方法。

背景技术：

随着新能源发电系统的大规模应用，电力电子设备已越来越多地被引入到电力系统中，从而产生过多电网谐波。电网谐波会导致电力电子设备过流，严重时将造成新能源发电系统大规模脱网等重大事故。理论上讲，电力电子设备可以在明确各类电网谐波信号特征全面信息（幅值、相位及频率）的前提下调整控制策略抑制过流，避免发生因保护动作而导致设备脱网的现象。

目前，电网谐波检测方法的研究主要集中在傅里叶变换（fouriertransform,ft）、小波变换（wavelettransform,wt）、经验模态分解（empiricalmodesdecomposition,emd）以及独立成分分析（independentcomponentanalysis,ica）4种方法。ft方法可将谐波信号特征的频率信息有效提取出来，但丢失了相位、幅值信息，无法得知扰动的发生时间、持续时长以及严重程度。wt方法可以同时获得谐波信号的相位与频率信息，从而得到比ft方法更好的辨识结果。但却具有对噪声敏感、鲁棒性差、频谱泄露等缺点。emd方法具有无需待测信号的先验信息，可对非平稳信号进行分析等优点，但其缺乏数学理论支撑且所得结果易发生模态混叠现象，使不同的检测结果含有相同的频率信息，不能够对谐波信号进行有效分离。ica方法能够在源信号未知的情况下，提取电网谐波信号特征信息，是一种功能强大的信号处理方法。但该方法要求输入为多路信号，而电网信号多数为单通道电压信号，因此，采用ica方法进行电网谐波检测需要先将单通道信号转换为多通道信号，增加了算法复杂度；同时由于ica的算法原理忽略信号的幅值信息，从而无法得知谐波扰动的严重程度。综上所述，现有方法并不能够将谐波信号的相位、幅值及频率信息全部提取出来。

变分模态分解（variationalmodedecomposition,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。可将电网电压信号中有效模态分量提出，模态分量的中心频率与原始信号含有谐波分量频率相同。即获得精确的模态分量，又有效避免模态混叠、冗余和丢失现象。

结论：通过仿真和真实电压信号进行实验，有效的验证了本发明的正确性与有效性。