混流式水輪機(jī)模型空化狀態(tài)與聲發(fā)射信號特征關(guān)系試驗(yàn)

劉 忠， 鄒淑云， 陳 瑩， 李志鵬

(1．長沙理工大學(xué) 能源與動力工程學(xué)院，長沙 410114； 2．中國水利水電科學(xué)研究院，北京 100038)

?

混流式水輪機(jī)模型空化狀態(tài)與聲發(fā)射信號特征關(guān)系試驗(yàn)

劉 忠1， 鄒淑云1， 陳 瑩2， 李志鵬1

(1．長沙理工大學(xué) 能源與動力工程學(xué)院，長沙 410114； 2．中國水利水電科學(xué)研究院，北京 100038)

為了獲取水輪機(jī)空化狀態(tài)下聲發(fā)射信號的特征參數(shù)及其變化規(guī)律，在混流式水輪機(jī)模型空化試驗(yàn)基礎(chǔ)上，采集了不同空化狀態(tài)時典型測點(diǎn)的聲發(fā)射信號，基于傳統(tǒng)分析方法提取了聲發(fā)射信號的事件計(jì)數(shù)、振鈴計(jì)數(shù)、頻率中心和電壓有效值等特征參數(shù)，分析了這些參數(shù)隨水輪機(jī)空化系數(shù)的變化關(guān)系.結(jié)果表明：隨著水輪機(jī)空化系數(shù)的逐漸減小，從無空化、空化初生到完全空化狀態(tài)，聲發(fā)射信號的各特征參數(shù)發(fā)生了較為明顯的規(guī)律性變化，整體上，事件計(jì)數(shù)不斷減小，振鈴計(jì)數(shù)、頻率中心和電壓有效值等不斷增大，變化幅度隨空化狀態(tài)所處階段有所不同.

混流式水輪機(jī)； 聲發(fā)射信號； 特征提?。?空化試驗(yàn)

水輪機(jī)內(nèi)發(fā)生空化時，容易產(chǎn)生高頻脈沖和輻射噪聲，引起水體壓力波動，并誘發(fā)水輪機(jī)能量特性下降、穩(wěn)定性惡化，嚴(yán)重時會造成流道表面破壞、部件壽命縮短[1].微型射流或沖擊波作用于葉片和流道內(nèi)壁，會使表面原子晶格產(chǎn)生錯位變形，發(fā)出0.02～1 MHz頻率范圍的聲發(fā)射(Acoustic Emission，AE)信號[2-3].對伴隨空化而生的聲發(fā)射信號進(jìn)行采集和分析，已被證實(shí)是一種可行的水輪機(jī)空化性能監(jiān)測方式[4].

在原型水輪機(jī)空化監(jiān)測研究方面，國內(nèi)外學(xué)者研究了不同空化狀態(tài)下聲發(fā)射信號特征與水輪機(jī)工作參數(shù)之間的關(guān)系.Escaler 等[5]對比了單機(jī)容量為11 MW的混流式水輪機(jī)在不同導(dǎo)葉開度下聲發(fā)射信號的頻譜特征，以及單機(jī)容量為65 MW的混流式水輪機(jī)在不同負(fù)荷下聲發(fā)射信號均方根的變化曲線.王輝斌等[6]提出了“脈沖重復(fù)率”和“聲強(qiáng)烈度” 2種聲發(fā)射信號的特征指標(biāo)，研究了東江水電站2號混流式水輪機(jī)的空化聲強(qiáng)烈度隨負(fù)荷變化的趨勢.Cencc 等[7]研究了抽水蓄能機(jī)組工作在水泵狀態(tài)時不同流量系數(shù)下聲發(fā)射信號均方根和空化系數(shù)之間的關(guān)系.Queiroz 等[8]分析了單機(jī)容量為175 MW的水輪機(jī)在空蝕修復(fù)前后的聲發(fā)射信號特征，建立了初步的判斷標(biāo)準(zhǔn).然而，原型水輪機(jī)的空化性能受水力設(shè)計(jì)參數(shù)、安裝條件和工作參數(shù)等的綜合影響，加上原型水輪機(jī)的真實(shí)空化狀態(tài)在現(xiàn)場無法準(zhǔn)確判斷，上述研究結(jié)果在水輪機(jī)空化監(jiān)測方面的適用性還有待商榷.

基于模型試驗(yàn)研究水輪機(jī)空化的聲發(fā)射監(jiān)測技術(shù)越來越受到關(guān)注.Rus等[9]研究了某軸流式水輪機(jī)模型的空化系數(shù)與聲發(fā)射信號均方根的變化關(guān)系.張俊華等[10]研究了某軸流式水輪機(jī)模型的空化系數(shù)與聲發(fā)射信號頻譜特征的變化關(guān)系.田浩等[11]研究了空化噪聲強(qiáng)度和噪聲脈沖計(jì)數(shù)與空化系數(shù)的關(guān)系.Schmidt 等[12]基于聲發(fā)射技術(shù)監(jiān)測模型水泵水輪機(jī)的前緣空化，分析了聲發(fā)射信號特征與圖像觀察結(jié)果之間的關(guān)系.薛延剛等[13]提出了基于希爾伯特-黃變換(HHT)的空化信號分析方法，并通過水輪機(jī)模型空化試驗(yàn)進(jìn)行了檢驗(yàn).但上述研究或多或少地存在選定的試驗(yàn)工況過于單一、空化類型與聲發(fā)射信號特征對應(yīng)關(guān)系過于籠統(tǒng)等問題.

因此，筆者在混流式水輪機(jī)模型空化試驗(yàn)基礎(chǔ)上，針對特定空化類型選擇試驗(yàn)工況，采集并提取聲發(fā)射信號的典型特征，分析混流式水輪機(jī)模型空化狀態(tài)變化過程中聲發(fā)射信號特征的變化規(guī)律.

1 試驗(yàn)過程

混流式水輪機(jī)模型空化試驗(yàn)在國內(nèi)某具有國際先進(jìn)水平、綜合精度<±0.2%的閉式水輪機(jī)模型試驗(yàn)臺[14]上進(jìn)行.試驗(yàn)過程中，采用與試驗(yàn)臺配套的水輪機(jī)性能測試系統(tǒng)采集、保存并計(jì)算得到反映水輪機(jī)模型能量特性及空化性能的參數(shù)，如流量、壓差、轉(zhuǎn)速和力矩等.采用閃頻儀觀測水輪機(jī)流態(tài).采用自行開發(fā)的聲發(fā)射信號采集與分析系統(tǒng)[15]采集來自不同測點(diǎn)的聲發(fā)射信號.其中，2個聲發(fā)射傳感器分別布置在導(dǎo)葉拐臂上和轉(zhuǎn)輪下環(huán)底部附近(如圖1所示)，以保證傳感器距離空化發(fā)生部位盡可能近，且與聲源之間的金屬厚度盡可能小.聲發(fā)射信號采集裝置的A/D轉(zhuǎn)換器分辨率為16 位，采樣頻率設(shè)置為2.0 MHz，帶通濾波頻率范圍為20～500 kHz，閾值電壓取0.05 V，所得信號分別用AE1和AE2表示.

圖1 聲發(fā)射傳感器布置

為了避免不同空化類型之間的相互干擾，在該混流式水輪機(jī)模型綜合試驗(yàn)基礎(chǔ)上，選取偏離設(shè)計(jì)工況且對應(yīng)于不同導(dǎo)葉開度a0和單位轉(zhuǎn)速n11的6組試驗(yàn)工況點(diǎn)(如圖2所示)，以便集中研究間隙空化過程中聲發(fā)射信號的變化規(guī)律.

每組工況點(diǎn)確定后，在試驗(yàn)中保持水頭和單位轉(zhuǎn)速不變，通過逐步抽真空改變尾水箱的壓力來實(shí)現(xiàn)水輪機(jī)空化狀態(tài)的調(diào)整.不同的單位流量、單位轉(zhuǎn)速和空化系數(shù)構(gòu)成一個測試工況，每個測試工況重復(fù)采集10組空化聲發(fā)射信號.將水輪機(jī)模型效率降低1%時的空化系數(shù)定義為臨界空化系數(shù)σ1，三個轉(zhuǎn)輪葉片上出現(xiàn)空泡判斷為空化初生，相應(yīng)的空化系數(shù)為初生空化系數(shù)σi.將空化試驗(yàn)結(jié)果與聲發(fā)射信號特征提取結(jié)果進(jìn)行對比分析.

圖2 試驗(yàn)工況示意圖

2 試驗(yàn)數(shù)據(jù)分析

2.1 各試驗(yàn)工況下的空化特性曲線

根據(jù)試驗(yàn)數(shù)據(jù)計(jì)算及閃頻儀觀測，得到各試驗(yàn)工況點(diǎn)的空化特性曲線和特征空化系數(shù)，分別如圖3和表1所示.

圖3 各試驗(yàn)工況點(diǎn)的空化特性曲線

工況點(diǎn)(σi，σ1)工況點(diǎn)(σi，σ1)n111?a01(0．1082，0．0554)n111?a02(0．1482，0．0753)n112?a01(0．1251，0．0582)n111?a02(0．1590，0．0691)n113?a01(0．1170，0．0685)n111?a02(0．1248，0．0848)

2.2 聲發(fā)射信號參數(shù)提取

針對空化試驗(yàn)中產(chǎn)生的連續(xù)型聲發(fā)射信號，采用傳統(tǒng)分析方法提取了其基本特征參數(shù)[16].

(1)事件計(jì)數(shù). 當(dāng)脈沖衰減波的包絡(luò)線超過閾值電壓時，產(chǎn)生一個矩形脈沖，此脈沖即稱為一次事件.它反映了聲發(fā)射事件的總量頻度.

(2)振鈴計(jì)數(shù). 聲發(fā)射信號超過閾值電壓而產(chǎn)生振蕩脈沖的個數(shù)稱為振鈴計(jì)數(shù).它反映聲發(fā)射信號的強(qiáng)度和頻度.一個事件可以包含若干個振鈴計(jì)數(shù).

(3)頻率中心. 在基于快速傅里葉變換得到的信號幅頻特性曲線上取上限頻率和下限頻率的幾何平均值，稱為頻率中心.它反映聲發(fā)射信號的主要頻率范圍.

(4)電壓有效值. 該值反映聲發(fā)射信號的能量，與頻率中心一樣，不受閾值電壓的影響，其表達(dá)式為：

(1)

式中：xrms為電壓有效值；n為信號電壓時間序列x的采樣個數(shù).

2.3 聲發(fā)射信號特征參數(shù)隨空化系數(shù)的變化關(guān)系

提取各試驗(yàn)工況點(diǎn)采集到的聲發(fā)射信號的事件計(jì)數(shù)、振鈴計(jì)數(shù)、頻率中心和電壓有效值等特征參數(shù).按相同n11、不同a0的工況點(diǎn)分類繪制特征參數(shù)隨水輪機(jī)空化系數(shù)變化的曲線，如圖4～圖6所示.

由圖4～圖6可以看出，隨著空化系數(shù)的逐漸減小， 2路聲發(fā)射信號AE1和AE2的各特征參數(shù)均發(fā)生了較明顯的變化，且具有一定的規(guī)律性.整體上，事件計(jì)數(shù)不斷減小，振鈴計(jì)數(shù)、頻率中心和電壓有效值均不斷增大，變化幅度隨空化狀態(tài)演變有所不同.這期間水輪機(jī)經(jīng)歷了從無空化到空化初生、再到臨界空化(已完全形成)的過程.

以圖4中單位轉(zhuǎn)速n111、導(dǎo)葉開度a01的工況點(diǎn)為例，隨著空化系數(shù)的減小，聲發(fā)射信號AE2的事件計(jì)數(shù)從1 281逐漸減小到47，振鈴計(jì)數(shù)從13 404逐漸增大到18 046，頻率中心從144 kHz增大到179 kHz，電壓有效值從0.125 V增大到2.492 V.

聲發(fā)射信號主要是由水體流動、設(shè)備旋轉(zhuǎn)以及水流與流道內(nèi)壁和水輪機(jī)轉(zhuǎn)輪葉片的相互作用產(chǎn)生的.在不斷減小尾水真空壓力但還未發(fā)生空化時，水流流態(tài)相對穩(wěn)定，聲發(fā)射信號的能量和頻度不斷增大，表現(xiàn)為事件計(jì)數(shù)減小，振鈴計(jì)數(shù)、頻率中心和電壓有效值增大.但就電壓有效值而言，增幅很小.

在空化初生(σ=0.108 2)時，由圖像觀測到水體中出現(xiàn)少許空泡.空泡的出現(xiàn)破壞了水流的平穩(wěn)性，局部出現(xiàn)氣液兩相流，水流內(nèi)部以及水流與流道邊壁和轉(zhuǎn)輪葉片之間的作用變得復(fù)雜.聲發(fā)射信號的能量和頻度快速增大，表現(xiàn)為事件計(jì)數(shù)繼續(xù)減小，振鈴計(jì)數(shù)、頻率中心和電壓有效值繼續(xù)增大.電壓有效值的增幅較之前明顯(從0.557 V增大到1.138 V).

隨著工況條件的繼續(xù)變化，水中空泡的數(shù)量不斷增加，伴隨著空泡的產(chǎn)生、成長、相互擠壓以及潰滅等，氣液兩相流程度加劇.聲發(fā)射信號的能量和頻度繼續(xù)增大，表現(xiàn)為事件計(jì)數(shù)繼續(xù)減小，頻率中心繼續(xù)向高頻方向移動，振鈴計(jì)數(shù)和電壓有效值繼續(xù)增大.就電壓有效值而言，增幅尤為明顯.到臨界空化(σ=0.055 4)時，空化已完全形成，電壓有效值達(dá)到最大(2.492 V).

(a) 事件計(jì)數(shù)-σ

(b) 振鈴計(jì)數(shù)-σ

(c) 頻率中心-σ

(d) 電壓有效值-σ

(a) 事件計(jì)數(shù)-σ

(b) 振鈴計(jì)數(shù)-σ

(c) 頻率中心-σ

(d) 電壓有效值-σ

(a) 事件計(jì)數(shù)-σ

(b) 振鈴計(jì)數(shù)-σ

(c) 頻率中心-σ

(d) 電壓有效值-σ

此后，隨著空化進(jìn)一步加劇，水中空泡的存在及其相互作用對聲發(fā)射信號的傳播產(chǎn)生了一定程度的衰減作用，表現(xiàn)為振鈴計(jì)數(shù)、頻率中心和電壓有效值均有小幅下降的趨勢.

盡管水輪機(jī)不同部位的聲發(fā)射信號特征參數(shù)的變化規(guī)律基本一致，但在相同工況下其數(shù)值還是存在一定差異的.其原因在于聲發(fā)射源發(fā)出的信號會經(jīng)過水體、流道固體邊壁、水輪機(jī)金屬部件等到達(dá)聲發(fā)射傳感器安裝位置，其傳播路徑不同，傳播過程中聲發(fā)射波的折射、散射和衰減程度也不同，從而造成不同測點(diǎn)處聲發(fā)射信號的差異.因此，聲發(fā)射傳感器位置和數(shù)量對空化狀態(tài)監(jiān)測結(jié)果的影響不容忽視.

3 結(jié) 論

(1)在混流式水輪機(jī)模型空化試驗(yàn)的基礎(chǔ)上，采集并提取了試驗(yàn)工況下的聲發(fā)射信號特征參數(shù)，分析了聲發(fā)射信號特征參數(shù)隨空化狀態(tài)演變的變化關(guān)系，進(jìn)一步促進(jìn)了水輪機(jī)空化監(jiān)測技術(shù)的研究.

(2)隨著水輪機(jī)空化系數(shù)的逐漸減小，聲發(fā)射信號的各特征參數(shù)均發(fā)生了較為明顯的變化，且具有一定的規(guī)律性.整體上，事件計(jì)數(shù)不斷減小，振鈴計(jì)數(shù)、頻率中心和電壓有效值均不斷增大，變化幅度隨空化狀態(tài)所處階段有所不同.

[1] KUMAR P, SAINI R P. Study of cavitation in hydro turbines-a review[J]. Renewable and Sustainable Energy Reviews, 2010, 14(1): 374-383.

[2] 楊晶, 李錄平, 饒洪德, 等. 基于聲發(fā)射檢測的閥門泄漏故障模式診斷技術(shù)研究[J]. 動力工程學(xué)報(bào), 2013, 33(6): 455-460, 483.

YANG Jing, LI Luping, RAO Hongde,etal. Diagnosis of valve leakage fault patterns based on acoustic emission detection[J]. Journal of Chinese Society of Power Engineering, 2013, 33(6): 455-460, 483.

[3] 劉忠, 鄒淑云, 晉風(fēng)華, 等. 聲發(fā)射技術(shù)在水電機(jī)組狀態(tài)監(jiān)測與故障診斷中的應(yīng)用研究綜述[J]. 水利水電技術(shù), 2011, 42(2): 49-51.

LIU Zhong, ZOU Shuyun, JIN Fenghua,etal. Overview on application of acoustic emission technique to operation monitoring and fault diagnosis of hydropower generating unit[J]. Water Resources and Hydropower Engineering, 2011, 42(2): 49-51.

[4] HUSIN S, ADDALI A, MBA D. Feasibility study on the use of the acoustic emission technology for monitoring flow patterns in two phase flow[J]. Flow Measurement and Instrumentation, 2013, 33: 251-256.

[5] ESCALER X, EGUSQUIZA E, FARHAT M,etal. Detection of cavitation in hydraulic turbines[J]. Mechanical Systems and Signal Processing, 2006, 20(4): 983-1007.

[6] 王輝斌, 吳長利, 鄒桂麗, 等. 水輪機(jī)故障聲學(xué)診斷技術(shù)研究及其應(yīng)用探討[J]. 大電機(jī)技術(shù), 2010(3): 45-50, 54.

WANG Huibin, WU Changli, ZOU Guili,etal. Acoustic fault diagnosis technology research and application in hydraulic turbine[J]. Large Electric Machine and Hydraulic Turbine, 2010(3): 45-50, 54.

[8] QUEIROZ F R, MEDEIRO E B, DONADON L V,etal. Acoustic emission evaluation of cavitation induced hydraulic turbine wear[C]//18th International Congress on Sound and Vibration 2011. Brazil, USA: International Institute of Acoustics and Vibrations, 2480-2484.

[9] RUS T, DULAR M,IROK B,etal. An investigation of the relationship between acoustic emission, vibration, noise, and cavitation structures on a Kaplan turbine[J]. Journal of Fluids Engineering, 2007, 129(9): 1112-1122.

[10] 張俊華, 張偉, 蒲中奇, 等. 軸流轉(zhuǎn)槳式水輪機(jī)空化程度聲信號辨識研究[J]. 中國電機(jī)工程學(xué)報(bào), 2006, 26(8): 72-76.

ZHANG Junhua, ZHANG Wei, PU Zhongqi,etal. Research on the cavitation identification of Kaplan turbine using acoustic signals[J]. Proceedings of the CSEE, 2006, 26(8): 72-76.

[11] 田浩, 于石生. 水輪機(jī)空化的超聲監(jiān)測技術(shù)研究[J]. 無損檢測, 2003, 25(5): 250-251, 278.

TIAN Hao, YU Shisheng. Research on ultrasonic detection of hydraulic turbine cavitation[J]. Nondestructive Testing, 2003, 25(5): 250-251, 278.

[12] SCHMIDT H, KIRSCHNER O, RIEDELBAUCH S. Cavitation measurements on a pump-turbine model[J]. Journal of Physics: Conference Series, 2015, 656(1): 012071.

[13] 薛延剛, 王瀚. 基于HHT的水輪機(jī)空化信號研究[J]. 水力發(fā)電學(xué)報(bào), 2015, 34(5): 147-151.

XUE Yan'gang, WANG Han. Investigation on turbine cavitation signals analysis based on Hilbert-Huang transform[J]. Journal of Hydroelectric Engineering, 2015, 34(5): 147-151.

[14] 馬素萍, 朱雷, 張海平, 等. 抽水蓄能機(jī)組水泵水輪機(jī)模型測試[J]. 水利水電技術(shù), 2012, 43(4): 94-98.

MA Suping, ZHU Lei, ZHANG Haiping,etal. Model test on pump-turbine of pumped-storage hydropower unit[J]. Water Resources and Hydropower Engineering, 2012, 43(4): 94-98.

[15] 鄒淑云, 劉忠, 李志鵬, 等. 離心泵空化狀態(tài)與聲發(fā)射信號特征關(guān)系試驗(yàn)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2014, 45(3): 45-48.

ZOU Shuyun, LIU Zhong, LI Zhipeng,etal. Experiment on relationships between centrifugal pump cavitation conditions and acoustic emission signal characteristics[J]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(3): 45-48.

[16] GROSSE C U, OHTSU M. Acoustic emission testing[M]. Berlin, Germany: Springer Science & Business Media, 2008.

Experiments on the Relationship Between Cavitation Status and Acoustic Emission Signal Features for a Francis Turbine Model

LIU Zhong1, ZOU Shuyun1, CHEN Ying2, LI Zhipeng1

(1.School of Energy and Power Engineering, Changsha University of Science and Technology,Changsha410114,China;2.ChinaInstituteofWaterResourcesandHydropowerResearch,Beijing100038,China)

To obtain the features of acoustic emission signals and their variation laws under the condition of hydraulic turbine cavitation, a series of cavitation experiments were conducted for a Francis turbine model, during which the acoustic emission signals were collected at typical measuring points under different cavitation conditions, while following features were extracted in a conventional analysis way, such as the event count, hit count, frequency center and the root mean square value of voltages, etc., so as to analyze the effects of cavitation coefficient on above features. Results show that with the decrease of turbine cavitation coefficient, the cavitation state changes from none to incipient, and further to fully-developed status; correspondingly, the features of acoustic emission signals vary in an obvious regularity; as a whole, the event count decreases gradually, while other three features increase, and their variation amplitudes depend on the stage in a cavitation process.

Francis turbine; acoustic emission signal; feature extraction; cavitation experiment

2016-01-25

2016-03-21

國家自然科學(xué)基金資助項(xiàng)目(51309034)；湖南省教育廳科學(xué)研究優(yōu)秀青年資助項(xiàng)目(14B004)；可再生能源電力技術(shù)湖南省重點(diǎn)實(shí)驗(yàn)室(長沙理工大學(xué))開放基金資助項(xiàng)目(2015ZNDL008)

劉 忠(1978-)，男，湖北漢川人，副教授，博士，研究方向?yàn)樗C(jī)械狀態(tài)監(jiān)測與故障診斷. 電話(Tel．)：0731-85258410；E-mail：drliu20003@126.com．

1674-7607(2016)12-1017-06

TK73

A 學(xué)科分類號：470.30