李文浩,余 波,胡錦蘅,劉 彬
(西華大學 能源與動力工程學院,成都 610039)
基于小波奇異性變換理論的離心泵空化分析
李文浩,余 波*,胡錦蘅,劉 彬
(西華大學 能源與動力工程學院,成都 610039)
以泵閥實驗臺為實驗平臺對離心泵進行空化實驗,并對不同流量下空化噪聲信號進行采集,獲得了空化性能曲線。運用小波奇異性變換理論對噪聲突變信號進行分析。結果表明,通過突變信號的模極大值線條數統(tǒng)計可對空化程度進行定性、定量分析,與實驗結果相吻合,為離心泵空化狀況檢測提供了可行的方法。
離心泵;空化噪聲;小波變換;奇異性檢測
離心泵運轉時,當葉輪葉片入口附近壓力小于等于液體輸送溫度下的飽和蒸汽壓力時,液體發(fā)生汽化產生氣泡。當氣泡隨液體流到液內壓力較高處時,外面的液體壓力高于氣泡內的汽化壓力,則氣泡會凝結潰滅形成空穴。這種由于由于液體局部壓力降低時,空泡在液體內部的形成、發(fā)展和潰滅的過程叫做空化[1]。
空泡潰滅時會輻射出尖脈沖,并形成沖擊波[2],是一種典型的突變信號,而突變信號中經常攜帶有比較重要的信息[3]。小波變換具有很好的時頻局部化特性,非常適合處理突變信號,這一特點已在很多領域得到了廣泛的應用[4-7]。離心泵空化的小波奇異性檢測方法就是通過檢測空化信號的奇異性來對空化進行檢測。
1.1 小波變換
數學上稱無限次可導函數是光滑的或沒有奇異性,若函數在某處有間斷或某階導數不連續(xù),則稱函數在此處有奇異性,該點就是奇異點。信號的突變點往往包含著重要的信息,一般用Lipschitz指數來描述信號的奇異性。Lipschitzα越大,函數越光滑。描述信號奇異性的一般定義如下:
定義:設有非負整數n(n≤α≤n+1),如果存在著兩個常數C>0和x0>0與一個n階多項式Pn(x),對于x∈(x0-δ,x0+δ)使得
(1)
則稱函數f在點x0是Lipschitzα的。
如果存在一個常數C與一個n階多項式Pn使得對于任意x0 函數f(x)的小波變換用Wf(s,x;Ψ)表示,Ψ為小波基,s尺度。設f(x)∈L2(R),[a,b]為R上的閉區(qū)間,存在0<α<1,對?ε>0,f(x)在(α+ε,b-ε)上是一致Lipschitz指數α的,當且僅當存在一個常數A,使對所有x∈(α+ε,b-ε)和任意的尺度s>0,有 (2) 1.2 小波變換的模極大值 設一個光滑的函數Φx,滿足下面的條件: (3) 以式(3 )的一階導數Ψ1(x)=dΦ(x)/dx作為小波,則有如下關系: (4) 2.1 實驗設計 本次離心泵空化實驗分別測試離心泵在100、92、80 m3/h工況下空化現象不同時期的噪聲信號,然后運用小波奇異性變換對噪聲信號進行處理,探究噪聲信號與離心泵空化發(fā)生程度之間的關系。 2.2 實驗介紹 本次實驗是在泵閥實驗臺上進行,研究對象采用的是清水式離心泵,型號IS100-80-125。該次實驗設備主要包括離心泵、水環(huán)式真空泵、電機、電磁流量計、扭矩測量儀以及水聽器。實驗裝置示意圖見圖1。由圖1可見,水流從水池經A進口流入,流經離心泵,并從B出口流入水池。 具體實驗步驟如下: 1)按照離心泵能量特性實驗方法啟動實驗裝置,并將進出口閥門至全開狀態(tài)。 2)保持進口閥門為全開狀態(tài),關閉出口閥門,使得流量計的讀數為設定流量工況下(做3個流量工況,100 m3/h,92 m3/h,80 m3/h),保持約5 min,測量此時的大氣壓及水溫,并記錄。再用照相機對準有機玻璃拍照,記錄下此時的水流狀態(tài)。 3)同時,利用NoiseA2.10噪聲測試系統(tǒng)(采樣頻率為1 MHz,采樣長度為32 k/幀,保存幀數10幀,參數可修改)連續(xù)采集10次,保存此時的信號。并同時測量此時的流量,離心泵的進出口壓力、轉速、轉矩各10次。采用這10組數據的平均值作為每個測量點的實際實驗參數。 4)逐漸關閉泵的進口閥門,誘發(fā)離心泵發(fā)生空化,調節(jié)出口閥門,保持流量計的讀數為設定流量工況下。重復2)、3)步驟,測量并記錄各個工況下所對應的參數,采集噪聲信號。改變實驗工況,重復上述實驗步驟,測量并記錄實驗數據。 2.3 水聽器的安裝 空化噪聲信號具有微弱、成分復雜、受環(huán)境干擾嚴重等特點,需采用高靈敏度的水聽器對信號進行采集,因此本次實驗采用的水聽器型號為RHS-10,其工作頻率為10 Hz~200 kHz,靈敏度為-210 dB。 圖1 實驗裝置示意圖Fig.1 Test device schematic diagram 水聽器探頭的布置方式對能否有效采集到空化噪聲信號影響很大,并且水聽器的安裝方式對水流的流態(tài)也有很大的影響。為了能夠有效地采集空化噪聲信號,考慮到現有實驗條件,水聽器探頭采用如下安裝方式: 1)在離心泵的進口安裝一段有機玻璃管。 2)將水聽器探頭與銅管相連。 3)在有機玻璃管上開孔,使得孔徑中心距離離心泵進口距離為85 mm。 4)將帶有探頭的銅管安裝在有機玻璃管徑的中心位置,見圖2。 圖2 水聽器的安裝方式Fig.2 Hydrophone installation method 根據上述實驗采集得到的實驗數據,根據式(7)計算離心泵的汽蝕余量,并繪制離心泵的空化性能曲線,見圖3。當汽蝕余量在某一臨界值以上時,空化對揚程的影響很小,當汽蝕余量小于這一臨界值時,離心泵揚程下降很快。 (7) 圖3 3種工況下的空化性能曲線Fig.3 Three kinds of working conditions of cavitation performance curve 在離心泵空化檢測時,采用運算量較小的db2小波作為小波基。對離心泵某一工況下采集的一組數據(圖3)進行5層的db2小波分析,結果見圖4。由圖4可見,各尺度的突變信號都很明顯,可認為這些突變點是對空泡潰滅時輻射出聲波脈沖的響應。 設置閾值,對第5尺度的突變點的模極大值條數進行統(tǒng)計。λ是調整系數,影響檢測的精度和效果,可在3~10選取,這里λ取值為3。圖4分別為3種工況下汽蝕余量所對應的模極大值的數目。離心泵流量為80 m3/h時不同汽蝕余量的空化情況照片見圖5。 以流量為80 m3/h為例(圖6),汽蝕余量為6.84 m時,極大值線數目為0,對照圖7可見有機玻璃管內沒有發(fā)生空化現象。當汽蝕余量減小到5.95 m時,極大值線數目為1,有機玻璃管中開始出現零散的單個氣泡。汽蝕余量降到為3.04 m時,極大值線數目為28,有機玻璃管中出現了許多連續(xù)的大氣泡,離心泵空化現象趨于嚴重。汽蝕余量降至2.4 m以下時,極大值線數目30以上,空化嚴重影響離心泵性能,揚程快速下降,見圖3。極大值線數目隨著離心泵汽蝕余量的減小而增加,它們之間有較好的曲線關系。 圖4 原噪聲信號Fig.4 Original noise signal 圖5 離心泵空化5層小波分析Fig.5 Centrifugal pump cavitation 5 layers of wavelet analysis 提出了一種基于小波奇異性變換理論的離心泵噪聲檢測方法,通過模極大值線的條數統(tǒng)計量與汽蝕余量的關系曲線,發(fā)現模極大值線的條數統(tǒng)計量隨著空化程度的加劇而增加,與有機玻璃管內氣泡變化吻合,證實了該方法能夠很好地反映空化發(fā)展的過程,可定性地判斷空化發(fā)展程度。 圖6 3個流量工況下(80、92、100 m3/h)時模極大值線Fig.6 Modulus maxima line of three traffic conditions(80、92、100 m3/h) 圖7 流量為80 m3/h時不同汽蝕余量的空化情況 Fig.7 Different cavitation NPSH of the flow rate of 80 m3/h [1]潘中永,袁壽其.泵空化基礎[M].江蘇:江蘇大學出版社,2013. Pan Zhongyong, Yuan Shouqi. The pump cavitation foundation [M]. Jiangsu, Jiangsu University Press, 2013. [2]Ross D. 水下噪聲原理 [M]. 北京:海洋出版社, 1983. Ross D. Underwater noise principle [M]. Beijing: Ocean Press, 1983. [3]Mallat S,Hwang W L.Singularity detection and processing with wavelets [J].IEEE Trans on Information Theory, 1992,38(2) : 617-643. [4]梁武科,羅興锜,張彥寧,等. 水力發(fā)電機組振動故障診斷系統(tǒng)中的信號預處理[J].水力發(fā)電學報,2003,22(3): 114-120. Liang Wuke, Luo Xingqi, Zhang Yanning, et al. Signal pretreatment of hydropower unit vibration fault diagnosis system [J]. Journal of Hydroelectric Power, 2003, 22 (3): 114-120. [5]蒲中奇,張 偉,施克仁,等. 水輪機空化的小波奇異性檢測研究[J]. 中國電機工程學報, 2005, 25(8): 105-109. Pu Zhongqi,Zhang Wei, Shi ke-ven, et al. Research on the wavelet singularity detection of hydraulic turbine cavitation [J]. Proceedings of the CSEE, 2005, 25 (8): 105-109. [6]賈清泉, 劉連光, 楊以涵, 等. 應用小波檢測故障突變特性實現配電網小電流故障選線保護[J]. 中國電機工程學報, 2001, 21(10):78-82. Jia Qingquan ,Liu Lianguang, Yang Yinhan, et al. Application of the wavelet to detect mutation characteristics to achieve distribution network fault protection of small current fault line selection [J]. Proceedings of the CSEE, 2001, 21 (10): 78-82. [7]卿 彪,余 波,陳凌平,等.CFD在調壓室設計計算中的可行性研究[J].黑龍江大學工程學報,2015,6(2):7-12. Qing Biao,Yu Bo,Chen Lingping,et al.Feasibilivy study of CFD on the design calculation of surge chamber[J].Journal of Engineering of Heilongjiang University,2015,6(2):7-12. [8]孔玲軍. Matlab小波分析超級學習手冊[M]. 北京:人民郵電出版社, 2014. Kong Lingjun. Matlab wavelet analysis super learning handbook [M]. Beijing: People’s Posts and Telecommunications Publishing House, 2014. Analysis of centrifugal pump cavitation based on wavelet singularity transformation theory LI Wen-Hao, YU Bo, HU Jin-Heng, LIU Bin (Energy and Power Engineering College of Xihua University, Chengdu 610039,China) Using pump valve test as test platform, the centrifugal pump cavitation experiment was carried out, and to collect the cavitation noise signal under different flow rates, cavitation performance curves were obtained. The singularity of noise signals was analyzed by the use of wavelet transform theory. The results showed that the degree of cavitation could be analyzed qualitatively and quantitatively through the modulus maxima of mutations in the signal lines for statistical number, and the results were in good agreement with the experimental results, a feasible method was provided for the detection of cavitation of centrifugal pump. centrifugal pump; cavitation noise; wavelet transform; singularity detection 10.13524/j.2095-008x.2015.04.069 2015-06-15; 2015-07-07 國家自然科學基金資助項目(51379179);四川省制造與自動化重點實驗室開放基金項目(szjj2012-04) 李文浩(1990-),男,江蘇徐州人,碩士研究生,研究方向:動力工程系統(tǒng)優(yōu)化與節(jié)能技術,E-mail:18482117165@163.com;*通訊作者:余 波(1965-),男,四川西昌人,教授,碩士研究生導師,研究方向:水利水電工程及自動化、計算機監(jiān)控與仿真技術等,E-mail:yubo@163.com。 TH311 A 2095-008X(2015)04-0078-052 離心泵空化噪聲實驗
3 實驗數據分析
4 結 論