壓氣機喘振聲音信號的快速傅里葉變換和小波變換分析

曹昳劼,臧樹升,葛 冰

(上海交通大學葉輪機械研究所,上海200240)

壓氣機發生喘振時會伴隨特有的噪聲,但是這種噪聲很可能被壓氣機本身的噪聲所淹沒。有經驗的操作人員在生產現場雖然可以根據聽到的噪聲中判斷壓氣機的運行是否正常,當壓氣機喘振時雖然可以分辨出喘振噪聲甚至聲源的位置,但是這種建立個人的經驗基礎上的判斷帶有主觀性,例如不同的人判斷的結果往往不一致,有時很難及時做出正確的決策。隨著計算機模式識別技術的提高,目前已經可以采用聲學測量和信號分析等方法,根據噪聲判斷壓氣機的喘振故障。例如居太亮等人利用麥克風聲音陣列的特性,建立了較為準確的麥克風陣列近場寬帶三維聲源定位算法[1];陳露明等人在對麥克風的采樣噪聲信號進行處理,提出了一種在DSP平臺上的硬件和軟件解決方案[2]。M.Morini等人對一個小型的六級軸流式壓氣機在穩定工況下以及變工況運行時的振動和聲音信號進行了對比研究,指出聲音信號在壓氣機監測上具有一定的可行性,并對聲音信號在壓氣機具體工況的測量及信號處理的方法進行了一定的闡述[3]。所有的這些監測方法關鍵是噪聲信號的識別和特征信號的提取,但是沒有對喘振一壓氣機運行時的危險工況進行特別的研究和討論。本文將在此基礎上,對壓氣機在進入喘振時的聲音進行研究和處理,通過消除噪聲干擾信號來得到喘振時聲音的特征信號,為實際生產中使用聲音信號監測壓氣機狀態以及故障診斷提供了良好的理論基礎和依據。

1 試驗裝置與數據采集系統

本文的試驗裝置由壓氣機,變頻電機和減速齒輪箱組成,俯視圖如圖1所示。

圖1 試驗裝置俯視圖

試驗裝置采用三個PCB公司的130D20麥克風傳感器收集壓氣機運行發出的聲音信息,傳感器的頻率響應為10 Hz到15 kHz,布置在壓氣機的周圍;采用庫力特公司的XTL-190M系列的壓力傳感器測量壓氣機出口的動態壓力;采用LDS Nicolet公司的Genesis系列GEN7t高速采集儀器進行聲音和動態壓力信號的采集,采樣頻率為1 k/s。試驗通過改變電機頻率調整壓氣機的轉速,通過閥門開度調整壓氣機的工況,壓氣機運行過程的壓力、溫度、流量和轉速等信號則通過PLC數據采集系統進行記錄。

2 試驗數據及分析

試驗主要是通過改變電機頻率以調節壓氣機轉速以及改變閥門開度以調節流量的方法來控制壓氣機的實際運行工況。

試驗以壓氣機的特性曲線為依據,從壓氣機的低工況開始做起在保證閥門開度不變的情況下通過逐步提高壓氣機的轉速來使其進入喘振區域。接著增大閥門的開度使流量增大,這時壓氣機的運行工況將逐漸遠離喘振線,然后保持閥門開度不變,再次提高轉速使得壓氣機的運行工況再次進入喘振區域,具體的試驗過程如表1所示。

在整個試驗測量的過程中,壓氣機的運行工況點有兩次進入了喘振區域。圖2所示是壓氣機運行工況兩次進入喘振區域的總體圖像,圖中的第一個通道顯示的信號為動態壓力傳感器的壓力信號,第二至第四的三個通道所顯示的信號為PCB的聲音傳感器采得的聲音信號。圖中所顯示的橫坐標為時間,縱坐標為壓力信號的幅值。從圖中可以清晰地看到,聲音信號在整個過程中有兩個明顯的變化的部分,這表征著壓氣機進入了喘振區域。

表1 壓氣機運行工況表

圖2 聲音信號整體圖像

圖3所示是壓氣機第一次進入喘振時聲音信號的局部圖,從圖中可以看到在進入喘振區域后,雖然存在著很多的噪聲干擾,但是進入喘振時的聲音信號存在一定的波動,與喘振之前的信號比較,進入喘振時的聲音信號出現了較大的振蕩波形的。

試驗給出了兩次喘振,每次喘振和喘振前的幅值計算結果如表2所示。

表2 兩次喘振的信號幅值 V

圖3 聲音信號局部圖像

從表中可以看到,每次喘振時,聲音傳感器所測得的幅值都發生了明顯的變化,聲音信號的幅值出現大幅的增加,增加的比例最大的超過180%。

3 信號的進一步處理

為了得到壓氣機運行工況進入喘振區域時聲音的頻譜特征信號,本文將采用快速傅里葉變換和小波變換對收集的噪聲信號進行處理分析,并將喘振特征信號從噪聲中提取出來。

3.1 快速傅里葉變換

傅里葉分析可以分析出采集所得的時域信號的頻譜和能量譜特征,本文把原始數據中每次壓氣機進入喘振和進入喘振前5 s的聲音信號截取出來進行對比處理,然后對截取的信號進行fft變換,第一次接近喘振邊界前5 s和喘振時,1#聲音傳感器所測得信號的fft變換見圖4和圖5。

圖4和圖5表征的是聲音信號在不同頻率段的振幅情況,可以看到壓氣機在喘振時,信號的fft變換在50 Hz之前頻率段出現了明顯了波峰,可以認為在第一次喘振的的頻率在50 Hz以內,喘振前和喘振時大于50 Hz信號幅值雖然有些差別但是波形相似,可以認定為背景噪聲信號。

圖4 第一次喘振前數據fft變換圖

圖5 第一次喘振時fft變換圖

圖6和圖7是第二次接近喘振前和喘振時,2#聲音傳感器所采集到得信號進行fft變換圖,與第一次喘振一樣,也可以認為喘振的頻率在50 Hz以內。

圖6 第三次喘振前數據fft變換圖

圖7 第三次喘振時數據fft變換圖

對照圖5和圖7以及表2可以看到,壓氣機發生喘振時不同傳感器所采得的信號幅值有一定的差別,但是在fft變換所得到的頻譜中并沒有明顯的差別,可以確定第一次喘振的頻率在50 Hz以內。

根據1#聲音傳感器的數據對第一次喘振的聲音信號進行功率譜的分析,結果如圖8所示。

從圖中可以看到,在50 Hz以下的信號頻率尤其是在20 Hz以下的低頻區域出現了明顯的能量聚集。

3.2 小波變換

小波分析的目的是對采得信號進行去噪處理,由于喘振發生在低頻區域,為了防止低頻成分的影響和滲入所以選擇中心頻率較高的小波基,即采用Daubechies(dbN)小波系中db8系的小波對信號進行處理,并使用7層信號分解,保證噪處理有著較好的效果。圖9是對2號聲音傳感器所采得在第一次進入喘振時的聲音信號進小波變換后得到的圖形,其中橫坐標顯示的為時間,縱坐標顯示的為振幅。

圖8 第一次喘振信號能量譜

圖中d1到d7分別為第一層到底七層小波分析的,a7為小波分析的近似信號。從a7信號中可以看出,開始一段信號較為平穩,表明壓氣機的運行點還沒接近喘振區,信號突然出現了劇烈的波動,表明已經進入了喘振區,可以確認,對低頻區域的聲音信號使用db8系列小波7層分解可以較好的去除噪聲,提取出聲音的特征信號。

圖9 第一次喘振信號小波分解

4 結論

通過此次壓氣機特性的試驗和分析結果,可以得到以下結論。

(1)在壓氣機運行工況進入喘振時,聲音信號有著明顯的變化,所以可以通過聲音信號識別壓氣機是否進入喘振狀態。

(2)壓氣機進入喘振狀態下的聲音信號的特征頻率范圍在50 Hz以下,根據相應的功率譜可以更準確地確定發生喘振的頻率,但振幅比功率譜的變化更為顯著。

(3)對進入喘振區的低頻區聲音信號進行小波分解時使用db8系列的小波進行7層分解可以有效地對聲音的特征信號進行提取,為生產中利用聲音信號監測壓氣機狀態,以及為故障的診斷提供了有價值的依據。

[1]居太亮,彭啟宗.基于任意麥克風陣列的近場聲源三維定位算法研究[J].信號處理,2007,23(2):231-234.

[2]陳露明.嘈雜背景下的聲源定位及語音分離實現技術的研究[D].成都:電子科技大學,2005.

[3]MORINI M,PINELLI M,VENTURINI M.Acoustic and Vibrational Analyses on a Multi-Stage Compressor for Unstable Behavior Precursor Identification[C]//ASME Turbo Expo 2007:Power for Land,Sea,and Air(GT2007),2007.