賀國(guó)，曹玉良，明廷鋒，蘇永生

(1.海軍工程大學(xué) 管理工程系，湖北 武漢 430033; 2.海軍工程大學(xué) 動(dòng)力工程學(xué)院，湖北 武漢 430033)

基于改進(jìn)倍頻帶特征的離心泵空化狀態(tài)識(shí)別

賀國(guó)1，曹玉良2，明廷鋒2，蘇永生2

(1.海軍工程大學(xué) 管理工程系，湖北 武漢 430033; 2.海軍工程大學(xué) 動(dòng)力工程學(xué)院，湖北 武漢 430033)

空化狀態(tài)識(shí)別是離心泵狀態(tài)監(jiān)測(cè)的難點(diǎn)之一。對(duì)離心泵的空化進(jìn)行了試驗(yàn)研究，采集了三種轉(zhuǎn)速時(shí)泵殼上兩個(gè)位置處的振動(dòng)信號(hào)，根據(jù)離心泵的特點(diǎn)對(duì)標(biāo)準(zhǔn)倍頻帶進(jìn)行改進(jìn)，基于改進(jìn)倍頻帶和標(biāo)準(zhǔn)倍頻帶構(gòu)建了振動(dòng)信號(hào)的特征向量，利用BP神經(jīng)網(wǎng)絡(luò)對(duì)離心泵的四類空化狀態(tài)進(jìn)行識(shí)別。研究表明：改進(jìn)倍頻帶比標(biāo)準(zhǔn)倍頻帶更能有效地提取離心泵空化振動(dòng)信號(hào)的特征；利用任意一處振動(dòng)信號(hào)的頻帶特征都能夠有效地識(shí)別揚(yáng)程降低大于3%的嚴(yán)重空化狀態(tài)，聯(lián)合利用兩處振動(dòng)信號(hào)的頻帶特征能夠有效地提高對(duì)正常狀態(tài)和空化程度較輕狀態(tài)的識(shí)別率。

離心泵； 空化狀態(tài)識(shí)別； 振動(dòng)信號(hào)； 頻帶特征； 神經(jīng)網(wǎng)絡(luò)； 倍頻帶

空化會(huì)導(dǎo)致流體機(jī)械的性能下降，引起振動(dòng)和噪聲，并導(dǎo)致過(guò)流部件的腐蝕破壞[1-2]。為了防止空化造成破壞，應(yīng)及時(shí)識(shí)別空化并進(jìn)行預(yù)警。如何有效地進(jìn)行空化的識(shí)別和監(jiān)測(cè)是研究的難點(diǎn)之一，國(guó)內(nèi)外很多學(xué)者對(duì)此進(jìn)行了研究。Escaler等對(duì)水輪機(jī)模型的振動(dòng)信號(hào)和瞬態(tài)壓力信號(hào)進(jìn)行了頻譜分析[3]。Hashmi等利用離心泵的瞬時(shí)角速度對(duì)空化進(jìn)行監(jiān)測(cè)[4]。Chini等發(fā)現(xiàn)某些頻率的聲壓級(jí)可以作為空化初生的特征，提出了一種基于噪聲頻譜的空化監(jiān)測(cè)方法[5]。Stopa等利用電機(jī)電流的變化計(jì)算離心泵轉(zhuǎn)矩的變化，設(shè)計(jì)了一種監(jiān)測(cè)離心泵空化狀態(tài)的計(jì)算機(jī)輔助工具[6]。然而，上述研究都主要是利用統(tǒng)計(jì)方法或者傅里葉變換，人工分析空化信號(hào)的異常，從而判斷泵內(nèi)空化狀態(tài)。

隨著智能識(shí)別技術(shù)的發(fā)展，人工神經(jīng)網(wǎng)絡(luò)在空化監(jiān)測(cè)和狀態(tài)識(shí)別方面已得到了一定的應(yīng)用。Mi等利用人工神經(jīng)網(wǎng)絡(luò)對(duì)豎直槽道內(nèi)的兩相流進(jìn)行了識(shí)別研究[7-8]。Giorgi等利用一個(gè)三層神經(jīng)網(wǎng)絡(luò)對(duì)孔口的空化狀態(tài)進(jìn)行了預(yù)測(cè)，發(fā)現(xiàn)其能夠識(shí)別溫度對(duì)空化的影響[9]。蘇永生等基于振動(dòng)和噪聲信號(hào)，利用BP神經(jīng)網(wǎng)絡(luò)和徑向基神經(jīng)網(wǎng)絡(luò)對(duì)噴水推進(jìn)泵的空化狀態(tài)進(jìn)行了識(shí)別研究[10]。然而，如何更高效地進(jìn)行空化特征提取和分類識(shí)別，仍是當(dāng)前空化監(jiān)測(cè)研究的難點(diǎn)問(wèn)題之一。求均值、求均方根和傅里葉變換等方法在特征提取中被廣泛應(yīng)用，然而空化特征分布在很寬的頻帶范圍內(nèi)，不同頻帶的特征不同，這些方法無(wú)法反映出信號(hào)在不同頻帶內(nèi)的特征。因此，本文基于離心泵空化試驗(yàn)采集的振動(dòng)信號(hào)，提出了一種結(jié)合頻帶特征和BP神經(jīng)網(wǎng)絡(luò)的離心泵空化狀態(tài)分類識(shí)別方法，將標(biāo)準(zhǔn)倍頻帶和改進(jìn)倍頻帶分別與BP神經(jīng)網(wǎng)絡(luò)相結(jié)合，對(duì)比分析了這兩種方法對(duì)離心泵四類空化狀態(tài)的分類識(shí)別效果。

1 離心泵空化實(shí)驗(yàn)

1.1 試驗(yàn)設(shè)備與過(guò)程

離心泵試驗(yàn)臺(tái)如圖1所示，主要包括兩臺(tái)離心泵、透明塑料管、水箱、球閥、蝶閥、電動(dòng)蝶閥和流量計(jì)等設(shè)備，除此以外，試驗(yàn)過(guò)程中還用到了變頻器、電動(dòng)控制臺(tái)、壓力傳感器、振動(dòng)傳感器、計(jì)算機(jī)和信號(hào)采集器。試驗(yàn)所用離心泵的比轉(zhuǎn)速為130，進(jìn)口直徑為126 mm，出口直徑為100 mm，葉片數(shù)為6片。

圖1 離心泵試驗(yàn)臺(tái)Fig.1 Experimental lab of centrifugal pump

振動(dòng)信號(hào)是狀態(tài)監(jiān)測(cè)中容易測(cè)量的信號(hào)之一，已被廣泛地應(yīng)用在流體機(jī)械空化監(jiān)測(cè)的研究中[11]。為了監(jiān)測(cè)空化所導(dǎo)致的泵殼的振動(dòng)，在泵殼的兩側(cè)V1和V2處安裝了兩個(gè)ICP型加速度傳感器，其安裝的方式和位置如圖2所示；信號(hào)采集使用的是BK 3560C振動(dòng)采集器，最高有效頻率fc為25.6kHz。

在進(jìn)行空化模擬試驗(yàn)時(shí)，離心泵機(jī)組2、球閥、蝶閥1和蝶閥2都保持關(guān)閉狀態(tài)，由離心泵1、流量計(jì)和電動(dòng)蝶閥1、2及其管路構(gòu)成開式試驗(yàn)系統(tǒng)。

在試驗(yàn)開始前，首先清洗水箱，注入干凈的清水，一次試驗(yàn)過(guò)程中保證水溫不發(fā)生變化。啟動(dòng)離心泵后先等待離心泵工作穩(wěn)定，然后對(duì)信號(hào)進(jìn)行頻譜分析，確定無(wú)異常后再進(jìn)行空化模擬試驗(yàn)。在試驗(yàn)過(guò)程中，使用變頻器調(diào)節(jié)離心泵的轉(zhuǎn)速，保證離心泵的轉(zhuǎn)速穩(wěn)定；通過(guò)減小進(jìn)口電動(dòng)蝶閥1的開度和降低水箱內(nèi)水面的高度來(lái)減小離心泵進(jìn)口的有效汽蝕余量，從而使離心泵內(nèi)部發(fā)生空化，在減小進(jìn)口電動(dòng)蝶閥1開度的同時(shí)增大出口電動(dòng)蝶閥2的開度以保持流量穩(wěn)定。

圖2 振動(dòng)傳感器Fig.2 Vibration sensors

1.2 試驗(yàn)結(jié)果

離心泵進(jìn)口有效汽蝕余量NPSHa和揚(yáng)程H的計(jì)算公式為

(1)

(2)

式中:Ps為進(jìn)出口斷面的絕對(duì)壓力，vin、vout為進(jìn)出口的流速，Pa為大氣壓，Pv為試驗(yàn)水溫下水的飽和蒸汽壓，Δh為進(jìn)出口的位置差。利用定轉(zhuǎn)速定流量法進(jìn)行了三種轉(zhuǎn)速的空化模擬試驗(yàn)，轉(zhuǎn)速1 475r/min時(shí)，離心泵的流量為100m3/h；轉(zhuǎn)速1 200r/min時(shí)，流量為80m3/h；轉(zhuǎn)速900r/min時(shí)，流量為54m3/h。試驗(yàn)測(cè)得的離心泵的空化性能曲線如圖3所示，從圖中可以看出，轉(zhuǎn)速1 475r/min時(shí)，揚(yáng)程最高；轉(zhuǎn)速900r/min時(shí)，揚(yáng)程最低；三種轉(zhuǎn)速時(shí)離心泵的揚(yáng)程都隨著有效汽蝕余量的減小而逐漸降低。

2 識(shí)別方法

2.1 提取頻帶特征

空化信號(hào)的特征分布在很寬的頻帶范圍內(nèi)，不同頻帶的特征不同，McKee等將倍頻帶特征與人工神經(jīng)網(wǎng)絡(luò)相結(jié)合進(jìn)行了水泵的空化狀態(tài)識(shí)別研究，并取得了較好的效果[12]。因此，本文也利用倍頻帶分析振動(dòng)信號(hào)的頻帶特征，設(shè)fl、fu分別為頻帶的下限和上限，fc為頻帶中心，則

fu=2mfl

(3)

(4)

當(dāng)m=1時(shí)為標(biāo)準(zhǔn)倍頻帶，根據(jù)國(guó)際標(biāo)準(zhǔn)ISO 532，標(biāo)準(zhǔn)倍頻帶是將人耳聽覺范圍內(nèi)的頻率劃分成10個(gè)固定的頻帶，10個(gè)頻帶的中心頻率分別為31.5、63、125、250、500、1 000、2 000、4 000、8 000 和16 000 Hz[13]。利用切比雪夫1型數(shù)字濾波器對(duì)振動(dòng)信號(hào)進(jìn)行濾波，得到10個(gè)頻帶的分信號(hào)，分別求取這些分信號(hào)的有效值，以這些分信號(hào)的有效值構(gòu)建特征向量。因此，利用標(biāo)準(zhǔn)倍頻帶進(jìn)行分析時(shí)，V1和V2處振動(dòng)信號(hào)的特征向量都為10維。

圖3 離心泵的空化性能曲線Fig.3 Cavitation performance curves of centrifugal pump

然而，離心泵為旋轉(zhuǎn)機(jī)械，不同離心泵的工作轉(zhuǎn)速不同，同一離心泵在不同工作環(huán)境時(shí)轉(zhuǎn)速也不同。為了適用于不同的轉(zhuǎn)速，將第2個(gè)頻帶的中心頻率設(shè)為離心泵的旋轉(zhuǎn)頻率f0，將低于第2頻帶的所有頻率設(shè)為第1頻帶，將最高頻帶的上限設(shè)為采集器的最高有效頻率fc，將其他頻帶的中心都設(shè)為轉(zhuǎn)頻的倍頻，改進(jìn)倍頻帶的頻帶分布如表1所示。

表1 改進(jìn)后的倍頻帶

與標(biāo)準(zhǔn)倍頻帶相比，改進(jìn)后的倍頻帶有兩個(gè)優(yōu)點(diǎn)：1)轉(zhuǎn)頻信號(hào)在旋轉(zhuǎn)機(jī)械中非常重要，而標(biāo)準(zhǔn)倍頻帶的下上限和中心頻率都是給定的，不能夠很好地體現(xiàn)旋轉(zhuǎn)機(jī)械的特征；改進(jìn)后的倍頻帶，第1頻帶只包含了轉(zhuǎn)頻的亞諧波，其他頻帶則圍繞在轉(zhuǎn)頻的倍頻周圍。2)標(biāo)準(zhǔn)倍頻帶的范圍局限于人類的聽覺范圍，只劃分了10個(gè)頻帶，而改進(jìn)后的倍頻帶提高了所分析頻率的下限和上限，使所分析的頻率范圍變?yōu)?～fc，文中fc=25.6kHz。在研究過(guò)程中，頻帶數(shù)k由fc和轉(zhuǎn)頻f0確定：

2k-3/2f0

(5)

根據(jù)離心泵的轉(zhuǎn)頻確定每個(gè)頻帶的上限和下限后，利用切比雪夫1型數(shù)字濾波器對(duì)第1頻帶進(jìn)行低通濾波，對(duì)其他頻帶進(jìn)行帶通濾波，得到頻率范圍不同的多個(gè)分信號(hào)，然后分別求取這些分信號(hào)的有效值，以這些分信號(hào)的有效值構(gòu)建特征向量。

改進(jìn)倍頻帶的頻帶數(shù)和頻帶范圍都隨著離心泵轉(zhuǎn)速的變化而變化，轉(zhuǎn)速1 475r/min和1 200r/min時(shí)，頻帶數(shù)都為12個(gè)，因此V1和V2處振動(dòng)信號(hào)的特征向量都為12維；轉(zhuǎn)速900 r/min時(shí)，頻帶數(shù)為13個(gè)，此時(shí)V1和V2處振動(dòng)信號(hào)的特征向量都為13維。

2.2 神經(jīng)網(wǎng)絡(luò)識(shí)別

人工神經(jīng)網(wǎng)絡(luò)是在生物神經(jīng)網(wǎng)絡(luò)的啟示下建立的數(shù)據(jù)處理模型，具有較強(qiáng)的自適應(yīng)性和學(xué)習(xí)能力，其中BP神經(jīng)網(wǎng)絡(luò)的應(yīng)用最為廣泛，因此本文也利用BP神經(jīng)網(wǎng)絡(luò)對(duì)離心泵的空化狀態(tài)進(jìn)行識(shí)別研究，識(shí)別的流程如圖4所示。

圖4 空化狀態(tài)識(shí)別流程Fig.4 Procedure of cavitation pattern recognition

在建立神經(jīng)網(wǎng)絡(luò)模型時(shí)，由特征向量的維度確定神經(jīng)網(wǎng)絡(luò)輸入層的節(jié)點(diǎn)數(shù)，由空化狀態(tài)的類別確定輸出層的節(jié)點(diǎn)數(shù)，設(shè)定隱含層節(jié)點(diǎn)數(shù)與輸入層節(jié)點(diǎn)數(shù)相等，由此就確定了神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)。

對(duì)于標(biāo)準(zhǔn)倍頻帶，當(dāng)只利用V1或V2處振動(dòng)信號(hào)的頻帶特征進(jìn)行識(shí)別時(shí)，神經(jīng)網(wǎng)絡(luò)輸入層和隱含層的節(jié)點(diǎn)數(shù)為10個(gè)；當(dāng)同時(shí)利用V1和V2的特征向量進(jìn)行聯(lián)合識(shí)別時(shí)，則將V1的特征向量和V2的特征向量合并構(gòu)成20維的特征向量，此時(shí)神經(jīng)網(wǎng)絡(luò)輸入層和隱含層的節(jié)點(diǎn)數(shù)為20個(gè)。

對(duì)于改進(jìn)倍頻帶，由于特征向量的維度隨著頻帶數(shù)的變化而變化，因此神經(jīng)網(wǎng)絡(luò)的輸入層和隱含層的節(jié)點(diǎn)數(shù)也隨著頻帶數(shù)的變化而改變。轉(zhuǎn)速1 475 r/min和1 200 r/min時(shí)，當(dāng)只利用V1或V2處的振動(dòng)信號(hào)進(jìn)行識(shí)別時(shí)，神經(jīng)網(wǎng)絡(luò)輸入層和隱含層的節(jié)點(diǎn)數(shù)為12個(gè)；當(dāng)利用V1和V2的特征向量進(jìn)行聯(lián)合識(shí)別時(shí)，神經(jīng)網(wǎng)絡(luò)輸入層和隱含層的節(jié)點(diǎn)數(shù)為24個(gè)。轉(zhuǎn)速900 r/min時(shí)，當(dāng)只利用V1或V2處的振動(dòng)信號(hào)進(jìn)行分析時(shí)，神經(jīng)網(wǎng)絡(luò)輸入層和隱含層的節(jié)點(diǎn)數(shù)為13個(gè)；當(dāng)同時(shí)利用V1和V2的特征向量進(jìn)行聯(lián)合識(shí)別時(shí)，神經(jīng)網(wǎng)絡(luò)輸入層和隱含層的節(jié)點(diǎn)數(shù)為26個(gè)。

在訓(xùn)練神經(jīng)網(wǎng)絡(luò)時(shí)，從每類空化狀態(tài)所對(duì)應(yīng)的樣本中隨機(jī)抽取60%的樣本作為訓(xùn)練樣本，剩余40%的樣本作為測(cè)試樣本。之后將訓(xùn)練樣本輸入到神經(jīng)網(wǎng)絡(luò)，利用誤差反向傳播算法，不斷地迭代調(diào)整網(wǎng)絡(luò)的權(quán)值和參數(shù)，直至收斂。最后將測(cè)試樣本輸入到訓(xùn)練好的網(wǎng)絡(luò)，比較網(wǎng)絡(luò)輸出的空化狀態(tài)與真實(shí)的空化狀態(tài)之間的差別，從而得到神經(jīng)網(wǎng)絡(luò)識(shí)別空化狀態(tài)的正確率。

3 空化狀態(tài)識(shí)別

3.1 空化狀態(tài)分類

揚(yáng)程和效率的變化是離心泵內(nèi)發(fā)生空化最明顯的標(biāo)志之一，Johann指出必須汽蝕余量NPSH1、

NPSH3可以作為空化狀態(tài)的標(biāo)準(zhǔn)[1]。而NPSH1、NPSH3分別表示離心泵揚(yáng)程降低1%和3%時(shí)的汽蝕余量，因此也可以以揚(yáng)程降低的程度為標(biāo)準(zhǔn)。在研究過(guò)程中，正常情況為狀態(tài)1，揚(yáng)程降低小于1%為狀態(tài)2，揚(yáng)程降低在1%～3%為狀態(tài)3，揚(yáng)程降低大于3%為狀態(tài)4。在利用神經(jīng)網(wǎng)絡(luò)進(jìn)行狀態(tài)識(shí)別時(shí)，四類空化狀態(tài)的標(biāo)簽設(shè)為[1 0 0 0]、[0 1 0 0]、[0 0 1 0]、[0 0 0 1]。

3.2 四類空化狀態(tài)識(shí)別

利用標(biāo)準(zhǔn)倍頻帶與改進(jìn)倍頻帶提取V1和V2處的振動(dòng)信號(hào)的頻帶特征，利用BP神經(jīng)網(wǎng)絡(luò)對(duì)離心泵的四類空化狀態(tài)進(jìn)行識(shí)別研究。研究過(guò)程中，根據(jù)所利用信號(hào)的不同又分為3類，即只利用V1處信號(hào)、只利用V2處信號(hào)和聯(lián)合利用V1、V2兩處信號(hào)三種情況。轉(zhuǎn)速1 475 r/min時(shí)，四類空化狀態(tài)的樣本數(shù)分別為100、200、200和200個(gè)，利用BP神經(jīng)網(wǎng)絡(luò)對(duì)四類空化狀態(tài)的識(shí)別率如表2所示。轉(zhuǎn)速1 200 r/min時(shí)，四類空化狀態(tài)的樣本數(shù)分別為100、300、200和200個(gè)，識(shí)別率如表3所示。轉(zhuǎn)速900 r/min時(shí)，四類空化狀態(tài)的樣本數(shù)分別為100、200、200和300個(gè)，識(shí)別率如表4所示。其中，狀態(tài)i(i=1,2,3,4)的識(shí)別率ri和平均識(shí)別率R的計(jì)算式為

(6)

(7)

式中：Ni表示測(cè)試樣本中屬于空化狀態(tài)i的樣本的數(shù)量，ni表示Ni中被正確識(shí)別出來(lái)的樣本的數(shù)量(ni≤Ni)。

表2 轉(zhuǎn)速為1 475 r/min時(shí)空化狀態(tài)識(shí)別的正確率

表3 轉(zhuǎn)速為1 200 r/min時(shí)空化狀態(tài)識(shí)別的正確率

表4 轉(zhuǎn)速為900 r/min時(shí)空化狀態(tài)識(shí)別的正確率

從表2～4的對(duì)比分析中可以看出，無(wú)論是只利用V1處振動(dòng)信號(hào)的頻帶特征或是只利用V2處振動(dòng)信號(hào)的頻帶特征，還是聯(lián)合利用V1和V2的頻帶特征，改進(jìn)倍頻帶的平均識(shí)別率都比標(biāo)準(zhǔn)倍頻帶的平均識(shí)別率高，說(shuō)明改進(jìn)倍頻帶比標(biāo)準(zhǔn)倍頻帶更能有效地提取離心泵空化振動(dòng)信號(hào)的特征。主要有兩方面原因：一是改進(jìn)倍頻帶的頻率范圍為0～fc，包括了全頻帶的信息，而標(biāo)準(zhǔn)倍頻帶只包含了人耳聽覺范圍內(nèi)的頻率信息，存在信息泄漏。二是，改進(jìn)倍頻帶的頻帶數(shù)、頻帶中心和頻帶上下限都隨著轉(zhuǎn)速的變化而變化，第1頻帶只包含了轉(zhuǎn)頻的亞諧波，其余頻帶都圍繞在轉(zhuǎn)頻及其倍頻周圍，能夠更好地提取離心泵的特征。

無(wú)論是標(biāo)準(zhǔn)倍頻帶還是改進(jìn)倍頻帶對(duì)狀態(tài)4的識(shí)別都接近100%，說(shuō)明結(jié)合振動(dòng)信號(hào)的頻帶特征和BP神經(jīng)網(wǎng)絡(luò)能夠有效地識(shí)別揚(yáng)程降低大于3%的嚴(yán)重空化狀態(tài)。主要是因?yàn)樵谶@種空化狀態(tài)時(shí)，離心泵內(nèi)部空化劇烈，大量的空泡不斷地產(chǎn)生和破滅，對(duì)泵殼產(chǎn)生劇烈沖擊，導(dǎo)致泵殼的振動(dòng)與其他狀態(tài)存在明顯不同，所以無(wú)論是利用標(biāo)準(zhǔn)倍頻帶還是改進(jìn)倍頻帶都能將其正確識(shí)別出來(lái)。

對(duì)于正常狀態(tài)1和空化程度很輕的狀態(tài)2，這兩種狀態(tài)非常相近、空化特征不明顯，而利用改進(jìn)倍頻帶對(duì)這兩種狀態(tài)的識(shí)別率也明顯高于標(biāo)準(zhǔn)倍頻帶，進(jìn)一步表明改進(jìn)倍頻帶比標(biāo)準(zhǔn)倍頻帶在提取離心泵空化特征方面更有效。

與只利用V1或V2處振動(dòng)信號(hào)的頻帶特征相比，聯(lián)合利用V1和V2處振動(dòng)信號(hào)的頻帶特征時(shí)平均識(shí)別率有明顯提高，不僅對(duì)強(qiáng)空化時(shí)的狀態(tài)4有很高的識(shí)別率，而且對(duì)正常狀態(tài)1和輕微空化的狀態(tài)2的識(shí)別率也都超過(guò)了90%，說(shuō)明聯(lián)合利用兩處信號(hào)比單獨(dú)利用一處信號(hào)更能有效地識(shí)別離心泵的空化狀態(tài)，主要是因?yàn)閂1處和V2處的振動(dòng)信號(hào)都只捕捉到了離心泵空化的部分特征，而兩處信號(hào)所包含的特征更加全面。

4 結(jié)論

1)改進(jìn)倍頻帶比標(biāo)準(zhǔn)倍頻帶更能有效地提取離心泵空化振動(dòng)信號(hào)的特征，改進(jìn)倍頻帶與神經(jīng)網(wǎng)絡(luò)相結(jié)合時(shí)空化狀態(tài)的識(shí)別率較高。

2)無(wú)論是標(biāo)準(zhǔn)倍頻帶還是改進(jìn)倍頻帶對(duì)狀態(tài)4的識(shí)別都接近100%，說(shuō)明結(jié)合振動(dòng)信號(hào)的頻帶特征和BP神經(jīng)網(wǎng)絡(luò)能夠有效地識(shí)別揚(yáng)程降低大于3%的嚴(yán)重空化狀態(tài)。另外，對(duì)于正常狀態(tài)1和空化程度較輕的狀態(tài)2，改進(jìn)倍頻帶對(duì)這兩種狀態(tài)的識(shí)別率也明顯高于標(biāo)準(zhǔn)倍頻帶。

3)對(duì)于正常狀態(tài)1和空化程度較輕的狀態(tài)2，單獨(dú)利用V1或V2處的振動(dòng)信號(hào)時(shí)識(shí)別率較低，而聯(lián)合利用V1和V2處振動(dòng)信號(hào)時(shí)識(shí)別率有明顯提高。

[1] JOHANN F G．Centrifugal pumps [M]. Berlin Heidelberg：Springer, 2014: 287-292．

[2] KIM K H, FRANC J P, CHAHINE G, et al. Advanced experimental and numerical techniques for cavitation erosion prediction [M]. Dordrecht: Springer, 2014: 3-6.

[3] ESCALER X, FARHAT M, AUSONI P, et al. Cavitation monitoring of hydroturbines: tests in a francis turbine model [C] // Sixth International Symposium on Cavitation, CAV2006. Wageningen, Netherlands: 2006: 1-5.

[4] AL-HASHMI S, GU F, LI Y, et al. Cavitation detection of a centrifugal pump using instantaneous angular speed [C] // ASME 7th Biennial Conference on Engineering Systems Design and Analysis. American Society of Mechanical Engineers. Manchester, UK, 2004: 185-190.

[5] CHINI S F, RAHIMZADEH H, BAHRAMI M. Cavitation detection of a centrifugal pump using noise spectrum [C]// ASME 2005 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference. American Society of Mechanical Engineers. [S.l.], 2005: 13-19.

[6] STOPA M M, CARDOSO F B J, MARTINEZ C B. Detection of incipient cavitation phenomenon in a centrifugal pump [C]//2012 IEEE Industry Applications Society Annual Meeting. Las Vegas, USA, 2012: 1-6.

[7] MI Y, ISHII M, TSOUKALAS L H. Vertical two-phase flow identification using advanced instrumentation and neural networks [J]. Nuclear engineering and design, 1998, 184(2): 409-420.

[8] MI Y, ISHII M, TSOUKALAS L H. Flow regime identification methodology with neural networks and two-phase flow models [J]. Nuclear engineering and design, 2001, 204(1): 87-100.

[9] DE-GIORGI M G, BELLO D, FICARELLA A. An artificial neural network approach to investigate cavitating flow regime at different temperatures [J]. Measurement, 2014, 47: 971-981.

[10] 蘇永生，王永生，段向陽(yáng). 基于多傳感器數(shù)據(jù)融合的噴水推進(jìn)泵空化分類識(shí)別[J]. 振動(dòng)與沖擊，2012，31(18): 93-96. SU Yongsheng, WANG Yongsheng, DUAN Xiangyang. Classification of cavitation in a waterjet pump based onmulti-sensor data fusion [J]. Journal of vibration and shock, 2012, 31(18): 93-96.

[11] 明廷鋒，曹玉良，賀國(guó)，等. 流體機(jī)械空化檢測(cè)研究進(jìn)展[J]. 武漢理工大學(xué)學(xué)報(bào)：交通科學(xué)與工程版， 2016， 40(2)： 219-227. MING Tingfeng, CAO Yuliang, HE Guo, et al. Research Progress of Cavitation Detection of Fluid Machinery[J]. Journal of Wuhan University of Technology： transportation science & engineering, 2016， 40(2)： 219-227.

[12] MCKEE K K, FORBES G L, MAZHAR I, et al. A vibration cavitation sensitivity parameter based on spectral and statistical methods [J]. Expert systems with applications, 2015, 42: 67-78.

[13] ISO 532. Method for calculating loudness level [S]. International organization for standardization, Switzerland, 1975.

本文引用格式：

賀國(guó)，曹玉良，明廷鋒，等. 基于改進(jìn)倍頻帶特征的離心泵空化狀態(tài)識(shí)別[J]. 哈爾濱工程大學(xué)學(xué)報(bào)， 2017， 38(8): 1263-1267,1302.

HE Guo, CAO Yuliang, MING Tingfeng, et al. Cavitation state recognition of centrifugal pump based on features of modified octave bands[J]. Journal of Harbin Engineering University， 2017， 38(8): 1263-1267,1302.

Cavitation state recognition of centrifugal pump based on features of modified octave bands

HE Guo1, CAO Yuliang2, MING Tingfeng2, SU Yongsheng2

(1.Department of Management Science, Naval University of Engineering, Wuhan 430033, China; 2.College of Power Engineering, Naval University of Engineering, Wuhan 430033, China)

Recognizing cavitation states is one of the challenges involved in monitoring the conditions of centrifugal pumps. First, the cavitation of a centrifugal pump was studied through experiments and the vibration signals at two points on the pump casing were collected at three rotation speeds. Subsequently, the standard octave bands were modified based on the characteristics of the centrifugal pump and the feature vectors of the vibration signals were extracted by using the modified and standard octave bands. Finally, four cavitation states of the centrifugal pump were recognized using a BP neural network. The experimental results demonstrated that the modified octave bands were more effective than the standard octave bands in terms of extracting the features of the cavitation vibration signals of the centrifugal pump. Heavy cavitation state with pump head decreased by 3% can be recognized from the band features of the vibration signals at either point. Using the band features of the vibration signals at the two points, the accuracy of recognizing the non-cavitation state and slight cavitation state was effectively improved.

centrifugal pump; cavitation states recognition; vibration signal; band feature; neural network; octave band

2016-11-06.

日期：2017-04-28.

國(guó)家自然科學(xué)基金資助項(xiàng)目(51306205)；湖北省自然科學(xué)基金項(xiàng)目(2015CFB700)；海軍工程大學(xué)博士生創(chuàng)新基金項(xiàng)目(4142C15K).

賀國(guó)(1965-)，男，教授，博士生導(dǎo)師; 曹玉良(1988-)，男，博士研究生.

曹玉良，E-mail：yuliangc@126.com.

10.11990/jheu.201611027

U664.33

A

1006-7043(2017)08-1263-05

網(wǎng)絡(luò)出版地址：http://www.cnki.net/kcms/detail/23.1390.u.20170428.1438.038.html