楊 波,黃 倩,付 強(qiáng)*,朱榮生
(1.江蘇大學(xué) 流體機(jī)械技術(shù)研究中心,江蘇 鎮(zhèn)江 212001;2.中國(guó)核電工程有限公司,北京 100840;3.核電泵及裝置智能診斷運(yùn)維聯(lián)合實(shí)驗(yàn)室,江蘇 鎮(zhèn)江 212013)
在我國(guó)化工、能源等領(lǐng)域中,曾接連出現(xiàn)過(guò)多起離心泵機(jī)械故障[1],其故障類型主要包括:機(jī)械密封失效、葉輪磨損導(dǎo)致的不平衡、軸承溫度異常偏高、齒輪箱異物、軸端密封導(dǎo)致泄漏等。目前,對(duì)離心泵各類故障進(jìn)行處理,以及對(duì)機(jī)組進(jìn)行日常維護(hù)已經(jīng)成為了臥式離心泵健康管理的重要組成部分[2]。
而對(duì)于臥式離心泵傳統(tǒng)的機(jī)械故障進(jìn)行必要的診斷也極為重要。目前,國(guó)內(nèi)外對(duì)于離心泵振動(dòng)信號(hào)的降噪研究已有很多。
HUANG N E等人[3]提出了一種基于經(jīng)驗(yàn)?zāi)B(tài)分解(empirical mode decomposition,EMD)算法的降噪方法,即將振動(dòng)故障信號(hào)分解成各階的本征模函數(shù)(IMF)分量,再?gòu)闹刑崛≌駝?dòng)信號(hào)中的故障特征。但在含有大量噪聲的背景下,使用EMD算法提取故障特征的效果會(huì)受到嚴(yán)重影響,存在模態(tài)混疊和末端效應(yīng)等問(wèn)題,導(dǎo)致提取的故障特征存在特征不明顯、誤差大或失真等問(wèn)題。
為了解決EMD分解的問(wèn)題,YEH J R等人[4]提出了互補(bǔ)集合經(jīng)驗(yàn)?zāi)B(tài)分解(CEEMD)降噪方法,通過(guò)添加噪聲的方法解決模態(tài)混疊等問(wèn)題。雖然CEEMD能解決EMD分解存在的問(wèn)題,但想要在強(qiáng)噪聲的背景下,提取出退化初期時(shí)振動(dòng)信號(hào)的微弱特征還有很大困難。因此,急需要對(duì)CEEMD分解后的分量做進(jìn)一步研究。
小波閾值降噪的方法是由DONOHO D L[5,6]提出的,該方法因?yàn)橛?jì)算量小而得到了廣泛應(yīng)用;但該方法本身存在缺陷,即小波硬閾值函數(shù)不連續(xù),降噪后可能會(huì)產(chǎn)生振蕩。小波軟閾值雖然具有較好的連續(xù)性,但處理后的小波系數(shù)和真實(shí)小波系數(shù)存在偏差,重構(gòu)信號(hào)時(shí)誤差增大,會(huì)導(dǎo)致精度下降。
因此,選取合適的小波閾值尤為重要。原磊明等人[7]提出了將小波閾值降噪用于單相接地故障時(shí)的電流信號(hào)中,發(fā)現(xiàn)其去噪效果明顯提高,得到了有效包含故障特征的向量集合。
對(duì)于臥式離心泵進(jìn)行故障診斷時(shí),模型分類算法的準(zhǔn)確率尤為重要。采用最鄰近(KNN)算法,對(duì)小樣本和少數(shù)據(jù)樣本進(jìn)行分類時(shí),其無(wú)需預(yù)估參數(shù),算法步驟較為簡(jiǎn)單。丁正生等人[8]提出通過(guò)聚類及數(shù)據(jù)重構(gòu),改進(jìn)KNN算法分類速度慢的缺點(diǎn),用改進(jìn)后的KNN算法對(duì)文本進(jìn)行向量化處理,再用特征向量對(duì)數(shù)據(jù)進(jìn)行分類,分類準(zhǔn)確率提高明顯。
綜上所述,筆者提出一種基于CEEMD改進(jìn)小波閾值去噪與優(yōu)化KNN算法的臥式離心泵機(jī)械故障診斷方法;即先通過(guò)CEEMD分解,計(jì)算出相關(guān)系數(shù),得到不相關(guān)向量,再用改進(jìn)小波閾值去噪,輸入到優(yōu)化后的KNN算法進(jìn)行故障分類,以實(shí)現(xiàn)對(duì)臥式離心泵的故障診斷。
該處筆者采用CEEMD結(jié)合改進(jìn)小波閾值的算法,對(duì)采集到的4種故障信號(hào)進(jìn)行分解,再用改進(jìn)小波閾值進(jìn)行去噪,最后,通過(guò)優(yōu)化后的KNN算法進(jìn)行故障分類。
故障診斷流程圖如圖1所示。

圖1 故障診斷流程圖
圖1中,首先對(duì)振動(dòng)信號(hào)進(jìn)行CEEMD分解,得到分解后的IMFs,計(jì)算其和原信號(hào)的相關(guān)系數(shù),得出不相關(guān)分量;再用改進(jìn)后的小波閾值去噪分解不相關(guān)分量,對(duì)高頻系數(shù)進(jìn)行去噪,將去噪后的不相關(guān)分量重構(gòu)高頻信號(hào);最后,使用粗糙懲罰法對(duì)相關(guān)分量進(jìn)行平滑處理,合成相關(guān)分量與不相關(guān)分量,得到了降噪信號(hào)。
對(duì)降噪后的信號(hào)進(jìn)行時(shí)頻特征向量提取,選取相應(yīng)的特征值構(gòu)建特征集,用于構(gòu)建模型測(cè)試訓(xùn)練,通過(guò)不斷地微調(diào)權(quán)重,對(duì)KNN算法進(jìn)行優(yōu)化,最終達(dá)到故障診斷的目的。
2.1.1 CEEMD原理
CEEMD算法的原理是,信號(hào)在EMD[9]分解之前添加成對(duì)的正、負(fù)白噪聲,以減輕EMD分解過(guò)程中出現(xiàn)的模態(tài)混疊、末端效應(yīng)及殘余噪聲的問(wèn)題[10]。
針對(duì)臥式離心泵不同工況下振動(dòng)信號(hào)分解后IMFs的選取重構(gòu)十分關(guān)鍵,直接影響故障信號(hào)的信噪比。合理選擇出關(guān)鍵的IMF分量,不僅可以提高故障信號(hào)的信噪比,還有助于信號(hào)的降噪。
2.1.2 相關(guān)系數(shù)法計(jì)算
在真實(shí)信號(hào)w(n)的基礎(chǔ)上,添加干擾e(n),取得加噪后的信號(hào)[11]:
y(n)=w(n)+e(n)
(1)
CEEMD去噪法其原理為將信號(hào)的相關(guān)分量進(jìn)行重構(gòu):
(2)
式中:imft(n)—y(n)分解出的第i個(gè)量;M—分解的IMF分量個(gè)數(shù);res(n)—解殘差。
重構(gòu)信號(hào)的另一種表示方法為:
(3)

(4)
式中:L—IMF分量的長(zhǎng)度;m—β(m)首次小于常數(shù)C所對(duì)應(yīng)的值。
β(m)呈逐漸下降的趨勢(shì),直到達(dá)到首個(gè)最小值。通過(guò)選取常數(shù)C計(jì)算得到kth的值,即相關(guān)分量首次出現(xiàn)的位置:
kth=arglast1≤m≤M{β(m)≥C}+1
(5)
式中:last1≤m≤M—矩陣中最后一個(gè)滿足條件的值。
在該處C取0.75,可以確定kth的值[12],則可推出前項(xiàng)為不相關(guān)分量,其余項(xiàng)相反。
2.1.3 對(duì)不相關(guān)分量、相關(guān)分量的處理
使用CEEMD對(duì)振動(dòng)信號(hào)進(jìn)行分解,通常認(rèn)為噪聲存在于不相關(guān)分量中,但不相關(guān)分量中存在有效信號(hào),使用CEEMD去噪法對(duì)信號(hào)去噪會(huì)丟失有效信號(hào),去噪效果不佳,使重構(gòu)信號(hào)失真。故筆者使用改進(jìn)小波閾值去噪法對(duì)不相關(guān)分量進(jìn)行去噪,使用粗糙懲罰法作為平滑濾波算法,對(duì)信號(hào)相關(guān)分量進(jìn)行平滑處理。
粗糙懲罰法的基本原理是在最小二乘法的原理上增加一個(gè)粗糙平滑項(xiàng)作為平滑函數(shù)[13],即:
(6)

粗糙平滑算法可以使擬合函數(shù)更加平滑,同時(shí)保證結(jié)果在平滑和失真間保持平衡。
2.2.1 理論分析
小波閾值降噪方法的核心在于閾值函數(shù)的構(gòu)建。閾值函數(shù)的不同表明對(duì)系數(shù)的估計(jì)方法的不同。閾值的大小也關(guān)系到降噪效果的好壞,只有選取合理的閾值才能使得降噪效果明顯,且有效信號(hào)不會(huì)丟失。
傳統(tǒng)的小波閾值為小波硬閾值函數(shù)和軟閾值函數(shù)[14],其表達(dá)式為:
(7)
硬閾值和軟閾值雖然在實(shí)踐中得到了廣泛的應(yīng)用,但該方法本身還存在一些缺陷,如容易出現(xiàn)不連續(xù)的間斷點(diǎn)和信號(hào)失真等問(wèn)題。
為了克服小波軟、硬閾值方法的缺點(diǎn),筆者在選取小波閾值函數(shù)的問(wèn)題上做出改進(jìn)。
新改進(jìn)的小波閾值表達(dá)式為:
(8)