基于快速譜峭度的沖擊能量法在水輪機空化評價中的應(yīng)用研究

林家洋，張民威，蘇疆東，劉 藝，王振鑫

(1.福建水口發(fā)電集團有限公司，福建福州350004; 2.北京中元瑞訊科技有限公司，北京100085)

0 引 言

水輪機轉(zhuǎn)輪流道內(nèi)及其過流部件的局部表面，經(jīng)常會發(fā)生空化而后引起輕微的空蝕，并在表面產(chǎn)生少量蝕點，進而會形成蜂窩狀的空蝕，甚至有的會使葉片穿孔或掉邊，會嚴重威脅到機組的安全運行[1]?？瘴g空化是導致水輪機運行效率下降、結(jié)構(gòu)破壞以及使機組壽命縮短的主要因素之一，深入研究空化問題，尋找能準確評價轉(zhuǎn)輪空化狀態(tài)的方法有著積極的意義[2]。

1 空化空蝕產(chǎn)生機理及空化信號特點

1.1 空化空蝕產(chǎn)生機理

空化發(fā)生在水輪機轉(zhuǎn)輪流道中局部壓力下降到臨界壓力(一般接近汽化壓力)時，水中氣核發(fā)展成為氣泡，從而使液相流體的連續(xù)性遭到破壞，變?yōu)楹瑲獾亩嗔?。氣泡中主耍充滿著液體的蒸汽以及從溶液中析出的氣體。當這些氣泡進入壓力較低的區(qū)域時，發(fā)育成長為較大的氣泡，當氣泡隨水流運動到壓力較高區(qū)域，氣泡將迅速凝縮并潰滅。因此，空化包括了氣泡的積聚、流動、分裂到潰滅的整個過程[3]。

水輪機空蝕主要發(fā)生在流速較高的轉(zhuǎn)輪通道區(qū)域和尾水管肘管。水輪機中的空蝕大體可以分為：翼型空蝕、局部空蝕、間隙空蝕和空腔空蝕。對于混流式水輪機，主要有翼型空蝕和空腔空蝕。混流式水輪機的翼型空蝕主要是由于葉片的局部低壓區(qū)形成的?？涨豢瘴g一般發(fā)生在低水頭、低負荷時，尾水管空蝕渦帶內(nèi)部含有大量的氣穴，當氣穴崩潰時，將使轉(zhuǎn)輪的泄水錐及尾水管的錐管段和肘管段遭受空蝕破壞[4]。對于軸流式水輪機，一般空蝕發(fā)生在轉(zhuǎn)輪葉片和水輪機導葉尾部、底環(huán)以及尾水管進口處[5]。

對于空化的破壞機理多年來國內(nèi)外學者進行了大量理論探索和實驗研究，可以歸納為：機械作用、化學腐蝕作用、電化學作用和熱力學作用等。

(1)機械作用理論。機械作用理論研究者認為過流壁面產(chǎn)生空化破壞是由于汽泡潰滅時產(chǎn)生微射流和沖擊波的強大沖擊作用所致。Hammitt通過計算和實測得出，游移型空泡潰滅時，近壁處微射流速度可達70～180 m/s；在過流壁面產(chǎn)生的沖擊壓力可高達705 MPa；微射流直徑約為2～3μm，沖坑直徑為2～20 μm，表面受到微射流沖擊次數(shù)約為100～1 000次/(s·cm2)：沖擊脈沖作用時間每次只有幾微秒，這樣高的沖擊作用將直接破壞物體表面而形成蝕坑，較小沖擊力的反復作用則引起物體表面疲勞破壞[6]。

(2)化學腐蝕理論。一般說來，化學腐蝕作用常常與機械空化作用互相促進，空化加速腐蝕，腐蝕也加速空化，二者聯(lián)合作用造成更嚴重的壁面破壞。

(3)電化學理論。在空泡潰滅時的高溫高壓作用下，金屬晶粒中形成熱電偶，冷熱端之間存在電位差，對金屬表面產(chǎn)生電解作用，造成電化學腐蝕。

(4)熱力學作用理論。當氣泡高速受壓后，汽相高速凝結(jié)，從而放出大量的熱，足以使金屬融化造成損壞。

目前，機械作用是較為公認的造成空化破壞的主要原因。

1.2 空化信號的特點

研究表明，水輪機中的空化發(fā)展到比較嚴重的程度時，才會導致尾水管壓力脈動信號特征的變化、機組振動和噪聲增強以及水輪機效率下降，而在空化發(fā)生的初始階段，液體沖產(chǎn)生的絕大多數(shù)空化的體積較小，空泡潰火時較短，輻射出髙頻率的沖擊脈沖聲信號，主要是超聲波信號。因此，監(jiān)測空化超聲波信有利于及早發(fā)現(xiàn)水輪機中的空化[1]。有關(guān)研究表明，空泡在固體表面潰滅時，像在彈性物體上潰滅一樣，產(chǎn)生振動。但其聲脈沖持續(xù)時間很短，是典型的沖擊信號，具有很寬頻譜，而同樣的空泡在液流中潰滅時，其聲沖擊脈沖持續(xù)時間增大，而頻帶變窄，此時，空化輻射出的超聲信號主要分布在30～120 kHz的頻率范圍內(nèi)。這個頻段很難由機械振動等低頻聲源產(chǎn)生，可以說是主要由空化產(chǎn)生的，避免了背景噪聲的干擾[7]。另外，隨著空化進一步發(fā)生時，液體中多數(shù)的空泡體積也隨著增大，空泡潰滅時間增加，輻射出的聲信號頻率逐漸向低頻區(qū)域過渡。因此，不同空化狀態(tài)下的引起的脈沖輻射的頻率并不固定。

綜上，可以看出，空化信號具有如下的特點:①空泡潰滅時產(chǎn)生的聲音是一些列的沖擊脈沖信號，沖擊脈沖信號的頻次以及脈沖間隔并不固定，是一個典型的隨機脈沖信號；②單次脈沖引起的聲音幅射頻率并不固定。

因此，通過檢測識別空泡潰滅時發(fā)出的高頻沖擊脈沖信號的頻次密度以及沖擊脈沖的幅值能量強度是識別轉(zhuǎn)輪空化程度的重要手段。

然而，傳統(tǒng)的方式是通過在水輪機的適當位置安放振動聲響傳感器，以測量在適當頻帶內(nèi)的噪聲強度，計算反映水輪機空化脈沖性能的噪聲脈沖頻率[8,9]。這種方法的最大問題是難以識別出混雜在各種振動、噪聲信號中的空化沖擊脈沖信號，而且由于不同空化狀態(tài)下脈沖輻射頻率是變化的，這就更增加了有效識別空化沖擊脈沖信號的難度[10,11]。

針對上述特點，本文提出采用基于快速譜峭度，結(jié)合數(shù)字包絡(luò)進行空化信號沖擊脈沖檢測，進而利用識別出的沖擊脈沖信號進行空化狀態(tài)評價的方法。

2 基于快速譜峭度的沖擊能量空化信號評價方法

2.1 譜峭度

譜峭度法(Kurtogram)是一個共振解調(diào)方法中帶通濾波器中心頻率和帶寬參數(shù)確定問題的有效方法，廣泛用于旋轉(zhuǎn)機械故障診斷[12]；Kurtogmm以濾波后時域信號的峭度值作為濾波效果度量指標，進而確定最優(yōu)的濾波器頻帶范圍。

峭度是描述波形尖峰度的一個指標，對沖擊信號非常敏感，可以對信號非平穩(wěn)強弱進行評判。但早期故障中信噪比很低，峭度作為一個全局指標不能正確地反應(yīng)出故障類型，不適用于強噪聲干擾環(huán)境下的故障檢測。而譜峭度方法實際上是計算每一條譜線上信號的非平穩(wěn)特性來確定非平穩(wěn)特征所在的頻率。研究表明該方法能夠很好地提取出被噪聲淹沒的故障特征信號，同時還能表示故障相應(yīng)頻率的峭度值。

譜峭度的計算過程是首先確定一個一定時間長度的窗口，并且在窗口內(nèi)對含噪信號進行Fourier變換得到其頻譜， 沿時間軸移動窗口，可以得到不同時段的頻譜。然后，對不同頻帶的譜統(tǒng)計其峭度， 得到的結(jié)果即為譜峭度。若頻帶選擇適當并且避開了干擾噪聲， 則可得到有用信號分量的峭度，該峭度與在時域得到的峭度一樣反映相應(yīng)信號分量的分布情況， 對信號的畸變敏感。譜峭度方法可以避開噪聲干擾在某些頻段上通過峭度檢測其不確定性，從而得到故障特征信息。

假定信號為

(1)

式中，H(t，ω)為被分析信號x(t)的時頻復包絡(luò)，可采用短時傅里葉變換計算得到。

根據(jù)譜的階矩定義，譜峭度可表示為

(2)

式中，C4y(ω)為信號y(t)的四階譜累積量；S(ω)為譜瞬時矩。

假設(shè)實測振動或聲音信號v(t)為

v(t)=x(t)+N(t)

(3)

式中，x(t)為有效的故障特征信號；N(t)為噪聲信號。實際運行過程，空泡潰滅過程產(chǎn)生沖擊脈沖，引起整個系統(tǒng)的共振(包括水體擊金屬結(jié)構(gòu))，所得故障信號x(t)的通用模型為

(4)

式中，h(t)為單個脈沖的脈沖響應(yīng)；Ak和τk分別為各個脈沖的幅值和發(fā)生時間。

從信號處理的角度看，譜峭度可以解釋為：理想濾波器組的輸出在頻率ω上計算得到的峭度值，于是有

(5)

式中，Kv(ω)為實測信號的譜峭度；Kx(ω)為故障特征信號的譜峭度；ρ(ω)為信噪比倒數(shù)。因此，某一頻率處信噪比高，則表示該處實測信號的譜峭度近似等于故障信號的譜峭度，從而可以找出最優(yōu)濾波頻帶。

2.2 譜峭度空化脈沖沖擊信號檢測

正如前文所述，譜峭度可以很好地分析非穩(wěn)定過程，如瞬時信號。而高度非穩(wěn)態(tài)的瞬時信號的峭度數(shù)值取決于頻率分辨率(Δf)。同時由于每一種瞬變時變現(xiàn)象對應(yīng)著一種最優(yōu)的頻率帶B(f，Δf)。因此，在實際的分析過程中，應(yīng)該找到最優(yōu)的頻率與頻率分辨率的信息，從而在這個區(qū)間內(nèi)，峭度值達到最大值，即可以找到相關(guān)的瞬態(tài)信息。

采用譜峭度方法對空化沖擊信號的檢測，帶通濾波器必須選擇包含3倍以上的空化噪聲頻率。正如本文所說，在超聲段頻率，空化信號主要集中在30～120 kHz，因此需選擇帶通濾波器在20～360 kHz。

本文提出的采用快速譜峭度方法對空化沖擊信號的沖擊脈沖檢測方法流程描述如下：

(1)根據(jù)空化噪聲信號超聲段的主要頻率范圍，結(jié)合采樣頻率確定最大分層數(shù)Nlevel。

(2)設(shè)定初始帶通濾波器范圍在20～360 kHz，采用1/3-二進濾波方式，按照確定的分解層數(shù)對信號進行濾波，根據(jù)式(2)計算各濾波信號的包絡(luò)譜帶通濾波峭度值K值，并生成直觀的改進譜峭度圖。最后選擇最大譜峭度值對應(yīng)頻帶的最優(yōu)帶通濾波器Bb(f，Δf)。

(3)選擇Bb(f，Δf)進行帶通濾波，通過數(shù)字包絡(luò)解調(diào)，求解空化信號的包絡(luò)信號。

(4)通過空化信號的包絡(luò)信號識別統(tǒng)計沖擊脈沖信號的頻次和強度。

詳細算法流程如圖1所示。

圖1 譜峭度空化脈沖沖擊信號檢測法流程

以下幾個圖是利用本方法對某軸流轉(zhuǎn)槳式機組的空化信號進行檢測的一個實例。其中圖2是該機組在低負荷區(qū)采集到的原始的空化聲音信號，在該負荷區(qū)存在著較為嚴重的空腔空化和葉形空化，然而由于原始信號中混疊有大量的振動、碰撞、噪聲等信號，直接地從原始信號中識別出空泡潰滅導致的沖擊脈沖信號幾乎不可能，或者難以準確識別和統(tǒng)計。

圖2 原始空化聲音信號

圖3是采用譜峭度計算出來的各分層分解的帶通頻帶對應(yīng)的峭度示意圖，其中不同顏色代表了不同的峭度值，顏色越紅峭度值越大。從圖中可以直觀地觀察到，當帶通濾波器選定在[66 666.664 Hz，75 000.0 0 Hz]時，峭度值達到最大值，因此[66 666.664 Hz，75 000.00 Hz]是該空化信號的最優(yōu)帶通濾波器。

圖3 譜峭度計算示意

圖4是將最采用最優(yōu)帶通濾波器對原始空化信號進行濾波以后的空化信號時域波形，從波形中可以清晰地觀察到多個高幅值的沖擊脈沖信號。

圖4 用最優(yōu)帶通濾波器濾波后的空化波形信號

圖5 是采用數(shù)字包絡(luò)解調(diào)技術(shù)對濾波后的空化聲音波形信號生成的包絡(luò)波形[13]，通過包絡(luò)波形信號可以更準確的對沖擊脈沖識別，可以檢測出脈沖出現(xiàn)的頻次、時刻以及沖擊的幅值。

圖5 帶通濾波后空化聲音信號的包絡(luò)波形

2.3 沖擊脈沖能量空化狀態(tài)評價方法

對空化聲音信號通過快速譜峭度計算獲得最優(yōu)帶通濾波器濾波，進而利用數(shù)字包絡(luò)解調(diào)算法獲得空化聲音的包絡(luò)波形，通過對包絡(luò)信號的識別可以獲得沖擊脈沖的數(shù)量、以及每個沖擊脈沖的時間寬度以及脈沖的幅值大小。

設(shè)空化聲音信號的包絡(luò)波形函數(shù)為xe(t)，那么假定第i個沖擊脈沖從t=t0開始，t=t1結(jié)束，如圖6所示。那么定義Epi為該沖擊脈沖的能量，即

圖6 單個空化聲音沖擊脈沖作能量積分示意圖

(6)

也就是Epi對應(yīng)了包絡(luò)波形函數(shù)xe(t)在[t0，t1]之間的面積。根據(jù)定義，Epi的值不僅與脈沖的幅值有關(guān)系，而且與脈沖的總時間寬度有關(guān)系，脈沖的幅值越高，Epi也越大；脈沖持續(xù)的時間越長，Epi也越大。因此通過Epi，能正確反映出不同空化階段，單個空化脈沖能量特性的變化。

設(shè)Lp_T為在給定時間內(nèi)(t=ΔT)通過上述方法檢測到的脈沖的總數(shù)量，Epi為給定時間內(nèi)第i個沖擊脈沖的能量，那么可定義

(7)

(8)

(9)

式中，Ep為單位時間內(nèi)的空化沖擊脈沖的總沖擊能量值；Ep_ave為單位時間內(nèi)的平均空化沖擊脈沖能量值；Lp單位時間內(nèi)的空化沖擊脈沖重復率。

綜上，根據(jù)轉(zhuǎn)輪空化發(fā)生的機理，當空化發(fā)生時，空泡潰滅將會產(chǎn)生大量的隨機沖擊脈沖，因此從測量的空化聲音信號中可以提取沖擊出現(xiàn)的頻率，和累計單位時間內(nèi)的沖擊總能量，通過這兩個評估指標就可以評價轉(zhuǎn)輪的空化狀態(tài)了。

一方面，空泡潰滅的沖擊脈沖頻次Lp越高，說明空泡越多，對轉(zhuǎn)輪葉片、尾水管等部位的破壞也越大。另外一方面，單位時間內(nèi)的沖擊總能量Ep越大，說明空泡越多而且脈沖的幅值越大，或者脈沖持續(xù)的時間越寬，對水輪機過流部件的破壞也越大。因此，借助于單位時間內(nèi)的沖擊總能量Ep和沖擊脈沖頻次Lp可以實現(xiàn)對水輪空化狀態(tài)和發(fā)展水平的分析評價。

3 算例分析

近年來，福建某水電廠軸流轉(zhuǎn)槳式機組轉(zhuǎn)輪空化日益嚴重，在機組轉(zhuǎn)輪輪轂、槳葉出水邊上下部、進水邊上部、轉(zhuǎn)輪室里襯等部位產(chǎn)生了空蝕，造成水輪機轉(zhuǎn)輪的流線破壞、機械構(gòu)件損壞等，影響了機組的安全穩(wěn)定運行。為了驗證本分析評價方法的的有效性，在該電站某機組上尾水管門外壁安裝了超聲寬頻傳感器，測量了該機組在不同負荷下的水聲信號，并利用本評價方法對測量數(shù)據(jù)進行了分析，詳細如下。

圖7是該機組尾水門超聲信號的頻譜陣列圖(濾除可聞聲段的信號，保留20 kHz以上的超聲信號)，右列坐標為機組負荷，橫坐標為頻率(Hz)，縱坐標為空化聲音幅值。從上圖看到，該機組在低負荷區(qū)(30 MW以下)，空化聲音信號主要集中在30～120 kHz 之間，而在高負荷區(qū)(110 MW以上)，這個頻段的聲音信號分布很小。

圖7 空化聲音信號的頻譜陣列示意

經(jīng)與該機組模型試驗對比和早期研究可知，該機組在低負荷區(qū)存在較為強烈的空腔空化和翼型空化；而在高負荷區(qū)，該機組空腔空化則減小、變?nèi)?。對比測量的聲學信號可知，在低負荷區(qū)產(chǎn)生的30～120 kHz信號，正是由于該機組的空腔空化和和翼型空化引起的聲音信號。根據(jù)其頻率分布特點，可以認為，該空腔空化和和翼型空化引起的空化沖擊脈沖信號引起的聲音主要分布在30～120 kHz之間。

圖8是采用譜峭度方法尋求最優(yōu)帶通濾波之后低負荷區(qū)空化聲音信號，從該圖中可以觀察到非常明顯的大幅值的沖擊脈沖信號。而在高負荷區(qū)，沖擊脈沖信號幅值和頻次則小很多，證明了采用該方法可以識別出有效的脈沖沖擊信號。

圖8 采用譜峭度進行最優(yōu)帶通濾波之后的空化聲音信號

表1為不同負荷區(qū)下空化沖擊脈沖能量及頻次統(tǒng)計結(jié)果，圖9、10分別為不同負荷區(qū)下空化沖擊脈沖重復率分布、能量分布。從表1以及圖9、10可以看出，在低負荷區(qū)單位時間內(nèi)的沖擊脈沖重復率最高達到7 000次/s以上，反觀高負荷區(qū)，沖擊脈沖重復率則低于3 000次/s以下；另外，在低負荷區(qū)單位時間內(nèi)的沖擊脈沖能量超過300 Pa·s；而在高負荷區(qū)，沖擊脈沖能量值不高于40 Pa·s。從而可以看出，通過沖擊脈沖重復率和沖擊脈沖能量可以較好的反映出水輪機的空化狀態(tài)和變化水平。

表1 不同負荷區(qū)下空化沖擊脈沖能量及頻次統(tǒng)計表(部分數(shù)據(jù))

圖9 不同負荷區(qū)下空化沖擊脈沖重復率分布

圖10 不同負荷區(qū)下空化沖擊脈沖能量的分布

4 結(jié) 語

本文研究了基于譜峭度的沖擊能量法在空化狀態(tài)評價的應(yīng)用，并通過實際測量數(shù)據(jù)的對比分析證明了該方法的有效性和實用價值。從實際意義來說，如果能找到有效的空化狀態(tài)評價方法，不僅可以實現(xiàn)對機組的早期故障預(yù)警，而且，也可為指導檢修提供重要的參考意義。然而，受限于測量數(shù)據(jù)、測量手段的限制，目前本方法上目前尚不能更精確地區(qū)分各類空化類型，亟需對更多的測試數(shù)據(jù)和和不同類型機組、不同工況下的數(shù)據(jù)進行分析和挖掘處理，以便修正和完善改進分析評價的更精確模型。另外，探索利用機器學習以及深度機器學習進行空化聲音信號的分類與識別以及空蝕的發(fā)展程度也將是一個重要的方向。