内燃机IFOA-VMD振动信号降噪预处理方法研究

李遥平,岳应娟,蔡艳平,王 旭,凤 林

(火箭军工程大学, 西安 710025)

摘要:为克服内燃机振动信号高度耦合、信噪比低,针对果蝇算法本身存在固定步长在搜索后期精度不高和气味判别函数不合理,提出一种改进的果蝇算法,以变分模态算法分解后的单分量模态信号最小功率谱熵作为评判函数,通过IFOA迭代搜索得出VMD分解层数k和惩戒函数α的最优组合,通过参数优化变分模态算法完成内燃机振动信号的降噪预处理。仿真结果表明,改进果蝇算法兼顾全局搜索能力,有效提升局部搜寻精度,参数优化的结果与仿真情况相符。将IFOA-VMD用于内燃机缸盖振动信号降噪预处理,通过平滑伪魏格纳-维利分布进行振动信号时频表征,可明显提高不同工况信号的信噪比,比经验模态分解方法降噪效果更好。

关键词:内燃机;振动信号;降噪预处理;改进果蝇算法;变分模态分解

内燃机作为一种机械装备主要的动力部件,其一旦发生故障,必然导致整个动力装置出现安全性和可靠性问题。所以,及时发现内燃机存在的故障,对于提高装备维修保养能力和效率具有重要意义。然而,内燃机是非常复杂的机械系统,目前主流的方法是通过分析其振动信号进行故障诊断,但振动信号为典型的非线性、非平稳信号,夹杂着大量的噪声信号[1]。因此对振动信号进行降噪预处理,减少其耦合度及提高信噪比,为故障诊断重要的一环。

目前,常用的振动信号降噪预处理方法有EMD分解[2]、LMD分解[3]、小波变换[4]等,但这些方法都存在一定的局限和不足。例如,EMD分解存在模态混叠及分解迭代条件缺乏标准等问题;小波变换由于符合测不准原理,因此在高频段处频率分辨率较低。2014年,Dragomiretskiy[5]等提出了一种完全非递归的信号分解方法——变分模态分解(VMD),该方法能够自适应地将复杂信号分解为多个本征模态函数。VMD一经提出,便广泛应用于一些领域,如文献[6]通过变分模态分解与形态差值滤波器结合,有效提取了轴承早期故障;文献[7]将变分模态分解和奇异谱分析结合,对微震信号进行二次降噪,更加有效地去除微震信号中的噪声干扰。

然而,由于VMD针对不同设备的振动信号设置参数存在区别,且参数设置不当对信号降噪预处理效果影响较大。因此,许多学者在此方面做了大量研究,文献[8]通过粒子群优化算法搜索变分模态分解算法的最佳参数组合,并将优化后的算法运用于提取轴承振动信号微弱特征频率信息,实现了滚动轴承早期故障的准确判别;文献[9]以合成峭度倒数和信息熵乘积最小值为评判准则选取VMD参数的最优组合,并将优化后的VMD与PSO-SVM结合,较好地实现轴承的故障诊断。尽管前人提出优化VMD的方法已经较为成熟,但依然需要不断的开发新的方法以达到更快的搜索速度和更好的搜索精度。基于此,笔者提出一种改进的果蝇优化算法以经VMD分解后的IMF分量局部最小功率谱熵为准则函数,迭代搜寻VMD的2个参数kα的最优组合。

果蝇优化算法(fruit fly optimization algorithm,FOA)是2011年6月由学者潘文超提出的一种具有全局性优化的智能算法[10]。由于该算法具有易理解、参数调节少、计算量小、搜索能力较强且搜索精度较高等优点,已成功应用于许多领域,例如文献[11]利用改进果蝇算法(IFOA)优化支持向量机(SVM)的方法来建立凝汽器真空预测模型,实现对凝汽器真空的准确预测;文献[12]利用多群果蝇优化算法,解决了多功能异构无人机协同任务规划问题;文献[13]通过改进果蝇算法优化BP神经网络,并应用于接地网腐蚀速率的预测,提高了预测的精度。但原始的果蝇优化算法仍存在步长在迭代搜索末期精度不高和气味判别函数设置不合理问题,因此,提出改进的果蝇优化算法并将其应用于VMD最优参数组合的选择上具有重要的研究意义。为了提高算法的优化性能,本文从迭代搜索步长设置及味道浓度判定函数两个方面对FOA进行改进,提出一种改进果蝇优化算法(improved fruit fly optimization algorithm,IFOA),使算法全局和局部寻优性能都得到提升并增加负值空间搜索能力,后将其成功应用于VMD的2个参数kα优化,通过仿真和实验验证都取得较好的效果。

1 改进果蝇优化(IFOA)算法

1.1 果蝇优化算法

果蝇优化算法(fruit fly optimization algorithm,FOA)是根据果蝇觅食行为而创新性提出的一种具有全局性寻找最优解的群智能优化算法。主要思想:首先,果蝇通过超强的嗅觉感知空气中的气味分子,后迅速地向食物源的方向飞去;然后,其他果蝇个体通过敏锐的视觉第一时间飞向同伴聚集的位置,共同寻找食物源具体位置;同时,防止判断错误情况发生,少部分果蝇种群随机向其他方向飞去。如此,反复的进行迭代更新直至到达最大迭代次数,具体步骤如下[10]

步骤1初始参数设置

设置FOA算法初始参数:最大迭代次数(Maxgen)、种群数量(Sizepop)、果蝇群体位置范围(R)、果蝇搜索步长(L)和每个果蝇种群包含的果蝇个体数量(dim)。

步骤2随机生成果蝇个体

Xaxis=R×rand()

Yaxis=R×rand()

(1)

步骤3果蝇个体被赋予随机的方向和搜索步长

Xi=Xaxis+L×rand()

Yi=Yaxis+L×rand()

(2)

步骤4计算果蝇个体与坐标原点的距离

(3)

步骤5计算味道浓度判定值

Si=1/Di

(4)

步骤6计算味道浓度值

Smilli=function(kα)

(5)

步骤7选择味道浓度最小的果蝇个体

[bestsmell,bestindex]=min(Smelli)

(6)

步骤8保留最佳的味道浓度和最优的果蝇个体位置

(7)

步骤9进行迭代寻优,判断是否到达最大迭代次数,否则进行步骤3~步骤8,到达后输出最优果蝇个体,即为寻优结果,与之对应的为VMD的参数k和α。

1.2 果蝇优化算法的改进

1) 搜索步长的优化

将搜索步长L代入到式(2)中可以看出,如果L设置过大果蝇个体后期搜索范围较大则寻优结果精度过低,反之L设置过小则寻优时间较长,甚至陷入局部最优。针对此问题,笔者提出一种改进的搜索步长方法——定变递进步长,具体如下:

(8)

式中: g为实时迭代次数,G为最大迭代次数,L为果蝇固定搜索步长,Lg为此次搜索步长,Lg-1为上次搜索步长。搜索步长L的具体变化曲线如图1所示。将式(8)代入到式(2)中可推导出果蝇个体分布搜索式(9)。

(9)

通过图1可清晰的观察到改进后的搜索步长在迭代前2/G次数为固定大步长,后2/G次数采取迭代递减步长方式。由式(9)可知,前半段(gG/2)时,L=1,则果蝇个体XY方向的搜索步长为此时步长取值较大且固定,因此果蝇个体搜索范围相对较大,可较大提高搜索效率,同时有效避免陷入局部最优;后半段(g>G/2)时,Lg=Lg-1-2/G×Lg-1,则果蝇个体XY方向的搜索步长为此时步长取值为迭代递减,因此果蝇个体搜索范围逐渐减小,结果值精度也逐渐提高。采取迭代步长的优点在于保持步长变化的平滑性,防止出现步长突变而导致陷入局部最优解问题的出现,同时也避免了最后一次迭代步长为零而导致算法出现bug情况。若未采用迭代递减方式而是采用定量递减的步长变化方式,则后半段每次步长减少量为固定步长L的2/G倍,最后一次迭代步长正好为零;正因为采用迭代方式,所以每次步长减少量为上一次步长的2/G倍,最后一次迭代步长为((G-2)/G)G/2×L,而不可能为零,因此该改进算法收敛且有效提高搜索精度。

图1 搜索步长与迭代次数关系曲线
Fig.1 Relation curve between search step size
and iteration times

2) 优化气味浓度判定函数

从式(3)和(4)可看出,算法中气味浓度判定函数采用倒欧式距离计算式,在实际应用在存在两个问题,一是计算出的气味浓度判定值都为正数,未考虑负数值的情况,无法解决负面空间问题,缺乏应对高纬度和复杂多变情况的求解能力;二是若采用倒欧式距离计算式,在搜索的前期判定值变化较大,存在陷入局部最优解的风险,同时步长(L)变化与判定值成反函数关系,不利于控制判断值的改变。综合以上两点问题,对气味浓度判定值计算式进行改进具体如式(10)所示:

Si=Di×sgn(Xi×Yi)

(10)

式(10)中采用欧氏距离与sgn函数的结合,在保证搜索步长与判定值同向的同时,也增加了对负值空间的搜索能力,扩展了算法的应用场景,提升了算法有效应对高纬度和复杂多变问题的解决能力。

2 变分模态分解(VMD)算法

为变分模态分解过程主要为2个过程:变分问题的构造和求解。

1) 变分问题的构造

变分问题实质就是迭代搜寻k个模态分量(IMF),且使得所有模态分量的带宽为最小。假设需处理的信号分解成k个模态分量,则可构造约束变分模型如下:

(11)

(12)

式中,δ(t)为Dirac函数;*为卷积;uk(k=1,2,…,K)为分解得到的模态分量;ωk为各个模态分量的中心频率;j为虚数;t为时间;h为输入信号。

2) 变分问题的求解

VMD算法将惩罚因子α和拉格朗日乘子λ(t)引入其中,从而使约束问题转化为非约束问题,以求得模型的最优解。增广后的拉格朗日函数Lp表达式为:

(13)

VMD算法是采取乘法算子交替的方法,通过不断迭代更新ωkn+1λn+1ukn+1,当满足迭代条件时停止循环,输出扩展拉格朗日函数的‘鞍点’,即为约束变分模型的最优解,需处理的信号便分解为k个模态分量。具体过程可查看文献[5]和文献[14]。

3 IFOA-VMD算法实现

IFOA优化VMD分解层数k和惩戒函数α时,引入局部最小功率谱熵[15]作为评判kα组合最优的标准。功率谱熵的大小反应分解后的IMF分量中包含噪音的多少,功率谱熵大则IMF分量中含有的故障相关的冲击信号不明显,混合较多的噪音;反之功率谱熵小IMF分量中包含的故障冲击特征信号明显,很清楚的反应在波形中,因此选用局部最小功率谱熵作为适应度函数是行之有效的,其计算表达式为:

(14)

在运用KFOA优化VMD算法具体问题中,相关参数的设置需根据实际情况而定,其中优化的参数为2个,因此dim设为2,分别用nS(i,1)对应α值,mS(i,2)对应k值。由于α值范围是1 000~10 000,所有设n=1 000;k值范围是1~10,所有设m=1,其算法流程如图2所示。

图2 IFOA-VMD算法流程框图
Fig.2 Flow chart of IFOA-VMD algorithm

4 仿真及实测验证

4.1 仿真验证

采用多分量合成信号X(t)作为仿真信号,其表达式如下:

X(t)=cos(40×π×t)+0.25cos(120×π×t)+

0.125cos(200×π×t)+0.

randn(size(cos(40×π×t)))

(15)

式中t=0.001∶1(间隔为0.001),randn为产生随机噪音的函数。分别运用定步长、变步长和定变递进步长这3种不同搜索步长的果蝇算法,进行VMD算法的参数优化。图3采用L=1定步长搜索,可发现在迭代15次后搜索精度没有变化,功率谱熵一直在2.6附近徘徊,是由于步长过大而导致精度较低;图4采用L=0.1定步长搜索,从图中可清晰看出在循环7次后,迭代结果基本未发生变化,且与仿真信号K=4实际不相符,发生此情况的原因为步长过小提前陷入局部最优解。

图3 搜索步长L=1时VMD参数优化及局部
最小功率谱熵迭代曲线
Fig.3 VMD parameter optimization and local minimum power
spectral entropy iteration curve when search step L=1

图4 搜索步长L=0.1时VMD参数优化及
局部最小功率谱熵迭代曲线
Fig.4 VMD parameter optimization and local minimum power
spectral entropy iteration curve when search step L=0.1

图5采用L=G-g/G均匀递减步长搜寻,迭代28次后,输出结果较好且与实际相符,仍存在循环次数多、寻优时间长的问题;图6步长采用式(8),从参数寻优结果可发现,在循环23次后,输出结果基本与实际相符且精度较好,可得出结论:搜索步长改进后的果蝇算法性能上优于定步长和均匀变步长果蝇算法,在保证全局寻优基础上,进一步提升局部寻优能力。

图5 搜索步长L=G-g/G时VMD参数优化及
局部最小功率谱熵迭代曲线
Fig.5 VMD parameter optimization and local minimum power
spectral entropy iteration curve when search step L=G-G/g

图6 定变递进步长时VMD参数优化及
局部最小功率谱熵迭代曲线
Fig.6 Constant variation progressive long-term VMD
parameter optimization and local minimum power
spectral entropy iteration curve

4.2 实验验证

1) 内燃机实验工况

本实验以6135型柴油机为研究对象,平台主要由控制台、传动轴、电机和柴油机4部分组成,如图7所示。预处理数据为内燃机第2缸盖的表面振动信号,采样频率为25 kHz,转速为1 500 r/min,采集过程中,内燃机保持空载运行。共设置气门间隙状况4种,具体如表1所示。其中气门间隙0.06 mm,0.3 mm和0.5 mm分别与排气阀气门间隙较小、正常和较大相对应,开口表示为气阀上开长(4 mm)×宽(1 mm)的缺口来模拟严重漏气。

图7 实验平台示意图
Fig.7 Experimental platform

表1 4种实验工况设置(mm)
Table 1 Setting of four experimental conditions

工况进气口排气口正常(1)0.300.30气门间隙较小(2)0.300.06气门间隙较大(3)0.300.50气门开口(4)0.30开口4×1

2) 内燃机缸盖振动信号预处理分析

首先运用改进的果蝇优化算法分别对4种工况下的振动信号进行迭代搜寻最优VMD参数kα组合,其原始振动信号时、频域曲线如图8~图11。由于果蝇优化算法中最大迭代次数(G)和果蝇种群数量(sizepop)这2个参数的设置一般对优化结果影响不大,因此可根据经验分别设为G=30,sizepop=30。具体参数优化结果如表2。

图8 4种工况的时、频曲线
Fig.8 Time frequency curves of four working conditions

表2 4种实验工况下果蝇算法优化结果
Table 2 Optimization results of Drosophila algorithm
under four experimental conditions

工况kα局部极小包络熵Epmin正常(1)42 6435.76气门间隙较小(2)42 7836.32气门间隙较大(3)42 8066.55气门开口(4)42 8636.76

因篇幅所限,只展示内燃机正常工况下的参数优化结果,如图9所示。

图9 正常工况果蝇算法参数优化曲线
Fig.9 Optimization curve of Drosophila algorithm parameters
under normal working conditions

利用参数优化后的VMD算法与SPWVD算法结合,表示出各工况下的振动信号时频图像:

图10 4种工况时频图
Fig.10 Time frequency diagram of four working conditions

图11 4种工况时频图
Fig.11 Time frequency diagram of four working conditions

同样,利用EMD算法与SPWVD算法相结合,其时频图像如下:

进行对比分析,可以明显发现各工况的振动信号经VMD处理后时域图像更加清晰,无信号重叠,尤其从图17和图20可看出振动信号图像发生明显的干扰重叠,这是由于VMD算法分解信号更加充分,而EMD算法不能将振动信号与噪音干扰信号完全分离开来,因此参数优化VMD算法相比EMD算法降噪效果更好。

5 结论

1) 内燃机缸盖产生的振动信号内夹杂着大量的噪声干扰信号,若直接对原始振动信号进行相关分析,难于得出理想分析结果。运用IFOA-VMD分解方法对原始信号进行降噪预处理,可取得较好的降噪作用,为故障分析打下基础。

2) 通过仿真对比分析,改进后的果蝇优化算法相比于定步长和均匀递减步长果蝇优化算法,在迭代效率和搜索精度都具有明显优势,既保证了全局搜索性又具备局部收敛性。气味浓度判定值的优化进一步扩大算法搜索空间,同时保证步长变化方向与气味浓度判定值一致。

3) 实验结果与基于EMD预处理方法相比,基于IFOA-VMD预处理方法生成的振动信号时频图像更加清晰,信号更加纯净,无明显信号干扰重叠情况,反映后者在去除干扰信号方面具有更优越的性能。

参考文献:

[1] 蔡艳平,徐光华,张恒,等基于内燃机振动信号的可视化识别诊断[J].振动与冲击,2019,38(24):150-157.

Cai Y P,Xu G H,Zhang H,et al.Visual recognition and diagnosis based on internal combustion engine vibration signal[J].Vibration and shock,2019,38(24):150-157.

[2] 王剑,梁斌,薛飞.一种基于白噪声EMD分解统计特性的水声信号去噪方法[J].四川兵工学报,2015,36(06):131-134.

Wang J,Liang B,Xue F.A denoising method of underwater acoustic signal based on EMD decomposition statistical characteristics of white noise[J].Sichuan Journal of ordnance industry,2015,36(06):131-134.

[3] 郭枫.局域能量逐层提取的LMD模态混叠抑制方法[J].兵器装备工程学报,2016,37(12):155-158.

Guo F.LMD mode aliasing suppression method for layer by layer extraction of local energy[J].Journal of ordnance and equipment engineering,2016,37(12):155-158.

[4] 柳振民,彭宗尧,郭力闻.小波变换在火箭遥测数据误码剔除中的应用[J].兵器装备工程学报,2017,38(05):94-98.

Liu Z M,Peng Z Y,Guo L W.Application of wavelet transform in error elimination of rocket telemetry data[J].Journal of Ordnance Engineering,2017,38(05):94-98.

[5] DRAGOMIRETSKIY K,ZOSSO D.Variational mode decomposition[J].IEEE Transactions on Signal Processing,2014,62(3):531-544.

[6] 艾澍海,张寿明.基于VMD和形态差值滤波器的特征提取算法[J].传感器与微系统,2018,37(09):151-154.

AI Shu H,Zhang S M.Feature extraction algorithm based on VMD and morphological difference filter[J].Sensors and Microsystems,2018,37(09):151-154.

[7] 苏梦哲,李发本,卢才武,等.基于VMD-SSA的微震信号降噪[J/OL].地球物理学进展:1-13[2020-11-09].http://kns.cnki.net/kcms/detail/11.2982.P.20201109.1008.024.html.

Su M Z,Li F B,Lu C W,et al.Noise reduction of microseismic signals based on vmd-ssa[J/OL].Progress in Geophysics:1-13[2020-11-09].http://kns.cnki.net/kcms/detail/11.2982.P.20201109.1008.024.html.

[8] 唐贵基,王晓龙.参数优化变分模态分解方法在滚动轴承早期故障诊断中的应用[J].西安交通大学学报,2015,49(05):73-81.

Tang G J,Wang X L.Application of parameter optimization variational modal decomposition method in early fault diagnosis of rolling bearing[J].Journal of Xi’an Jiaotong University,2015,49(05):73-81.

[9] 刘臻,彭珍瑞.信息熵和合成峭度优化的VMD和PSO-SVM的轴承故障诊断[J/OL].机械科学与技术:1-8[2020-11-09].https://doi.org/10.13433/j.cnki.1003-8728.20200239.

Liu Z,Peng Z R.Bearing fault diagnosis based on VMD and PSO-SVM optimized by information entropy and synthetic kurtosis[J/OL].Mechanical science and technology:1-8[2020-11-09].https://doi.org/10.13433/j.cnki.1003-8728.20200239.

[10] Pan W T.A new fruit fly optimization algorithm:taking the financial distress model as an example[J].KnowledgeBased Systems,2012,26:69-74.

[11] 葛晓霞,赵舒莹,肖洪闯,等.基于改进果蝇算法优化SVM的凝汽器真空预测[J].热能动力工程,2020,35(11):39-45.

Ge X X,Zhao S Y,Xiao H H,et al.Condenser vacuum prediction based on improved Drosophila algorithm optimization SVM[J].Thermal power engineering,2020,35(11):39-45.

[12] Aviation-Unmanned Aerial Vehicle; New Findings from Harbin Institute of Technology in the Area of Unmanned Aerial Vehicle Reported(Solving the Multi-Functional Heterogeneous UAV Cooperative Mission Planning Problem Using Multi-Swarm Fruit Fly Optimization Algorithm).2020:172.

[13] 程宏伟,高莲,于虹,等.基于改进果蝇算法优化BP神经网络的接地网腐蚀速率预测模型[J/OL].电测与仪表:1-10[2020-12-08].http://kns.cnki.net /kcms/detail/23.1202.TH.20200925.1733.006.html.

Cheng H W,Gao L,Yu H,et al.Corrosion rate prediction model of grounding grid based on improved Drosophila algorithm optimized BP neural network[J/OL].Electrical measurement and instrumentation:1-10[2020-12-08].http://kns.cnki.net/kcms/detail/23.1202.TH.20200925.1733.006.html.

[14] LI Changxi.Fault diagnosis of high-speed train bogie based on variational mode decomposition[D].Chengdu:Southwest Jiaotong University,2018:25-38.

[15] 闫晓玲,董世运,徐滨士.基于最优小波包Shannon熵的再制造电机转子缺陷诊断技术[J].机械工程学报,2016,52(04):7-12.

Yan X L,Dong S Y,Xu B S.Rotor defect diagnosis technology of remanufactured motor based on optimal wavelet packet Shannon entropy[J].Journal of mechanical engineering,2016,52(04):7-12.

Research on Noise Reduction Pretreatment Method of IFOA-VMD Vibration Signal of Internal Combustion Engine

LI Yaoping, YUE Yingjuan, CAI Yanping, FENG Lin

(Rocket Force University of Engineering, Xi’an 710025, China)

Abstract: In order to solve the problems of high coupling and low signal-to-noise ratio of internal combustion engine vibration signals, a noise reduction preprocessing method of internal combustion engine vibration signals based on the combination of improved Drosophila algorithm (IFOA) and variational modal algorithm (VMD) was proposed. Aiming at the problems of fixed step size, low precision in the later stage of search and unreasonable odor discrimination function of Drosophila algorithm, an improved Drosophila algorithm was proposed. Then, the minimum power spectral entropy of single component modal signal decomposed by variational modal algorithm was used as the evaluation function, and the VMD decomposition layers k and punishment function were obtained through IFOA iterative search α The noise reduction preprocessing of internal combustion engine vibration signal was completed by parameter optimization variational modal algorithm. The simulation results show that the improved Drosophila algorithm not only takes into account the global search ability, but also effectively improves the local search accuracy. The results of parameter optimization are consistent with the simulation. IFOA-VMD was applied to the noise reduction preprocessing of internal combustion engine cylinder head vibration signal, and then the vibration signal was characterized by smoothing pseudo Wigner Ville distribution. It can be clearly found that the signal-to-noise ratio of signals under different working conditions is significantly improved, and the noise reduction effect is better than that of empirical mode decomposition method.

Key words internal combustion engine; vibration signal; noise reduction pretreatment; improved drosophila algorithm; variational mode decomposition

收稿日期:2020-01-06;修回日期: 2020-12-19

基金项目:陕西省自然科学基金项目(2019JQ-712)

作者简介:李遥平(1994—),男,硕士研究生,E-mail:838881651@qq.com。

通信作者:岳应娟(1973—),女,博士,教授,E-mail:yingjuanyue@163.com。

doi: 10.11809/bqzbgcxb2021.09.046

本文引用格式:李遥平,岳应娟,蔡艳平,等.内燃机IFOA-VMD振动信号降噪预处理方法研究[J].兵器装备工程学报,2021,42(09):286-292.

Citation format:LI Yaoping, YUE Yingjuan, CAI Yanping, et al.Research on Noise Reduction Pretreatment Method of IFOA-VMD Vibration Signal of Internal Combustion Engine[J].Journal of Ordnance Equipment Engineering,2021,42(09):286-292.

中图分类号:TK418;TN911.6

文献标识码:A

文章编号:2096-2304(2021)09-0286-07

科学编辑 胡绍林 博士(广东石油化工学院教授、博导)责任编辑 唐定国