基于人工神經(jīng)網(wǎng)絡(luò)的離心泵空化故障診斷

摘 要：為了研究離心泵空化故障診斷問題，探究人工神經(jīng)網(wǎng)絡(luò)在該問題上的預(yù)測效果，通過數(shù)值仿真的方法對離心泵流場進行模擬，采集不同狀態(tài)下流場內(nèi)各點壓力值及工作點作為輸入特征，以旋轉(zhuǎn)區(qū)域氣體體積分?jǐn)?shù)為標(biāo)簽特征，對離心泵空化狀態(tài)進行神經(jīng)網(wǎng)絡(luò)建模。使用LSTM和一維卷積網(wǎng)絡(luò)處理時序數(shù)據(jù)，并在特征提取階段添加正則化損失函數(shù)以保證網(wǎng)絡(luò)稀疏性。最終模型在測試集上的分類任務(wù)準(zhǔn)確率達到95%以上，能夠有效地對離心泵空化程度進行診斷。

關(guān)鍵詞：離心泵；空化；人工神經(jīng)網(wǎng)絡(luò)；自編碼器

中圖分類號：TP277.3 文獻標(biāo)志碼：B 文章編號：1671-5276（2024）04-0151-03

Centrifugal Pump Cavitation Fault Diagnosis Based on Artificial Neural Network

LIU Chaoxi¹， YE Zhifeng¹， WANG Bin¹， YAN Shebin²

（1. Nanjing University of Aeronautics and Astronautics， Nanjing 2100016， China;2. AECC Guizhou Honglin Aro-engine Control Technology Co.， Ltd.， Guiyang 550009， China）

Abstract：In order to study the cavitation fault diagnosis of centrifugal pump and explore the prediction effect of artificial neural network concerning the fault diagnosis， the flow field of centrifugal pump was simulated by numerical simulation， the pressure values and working points of various points in the downstream field in different states were collected as the input characteristics and the volume fraction of the gas in the rotating region was taken as the label characteristic to conduct the neural network modeling for the cavitation state of the centrifugal pump.The LSTM and one-dimensional convolutional network were used to process the time series data， and the regularization loss function was added in feature extraction stage to ensure network sparsity.The accuracy rate of classification task of the trained model on the test set exceeded 95%， which can effectively diagnose the cavitation degree of centrifugal pump.

Keywords：centrifugal pump; cavitation; artificial neural network; autoencoder

0 引言

在航空航天領(lǐng)域中，離心泵由于其結(jié)構(gòu)簡單、質(zhì)量較輕以及能承受較高轉(zhuǎn)速等優(yōu)點，在燃油、液壓系統(tǒng)中獲得廣泛應(yīng)用。而在其日常運行中，空化故障是其最常發(fā)生的故障之一。在離心泵內(nèi)部流場中，當(dāng)局部的壓力低于液體的汽化壓力時，液體會氣化而形成氣泡，并在流到高壓區(qū)后發(fā)生破裂，對周圍產(chǎn)生沖擊，這種現(xiàn)象就稱為空化。泵空化時由于內(nèi)部的氣泡會占據(jù)流道，干擾正常的流動，影響液泵的加壓性能，并且氣泡不斷地產(chǎn)生與破滅會對壁面持續(xù)產(chǎn)生沖擊，使得葉片發(fā)生損傷或腐蝕^[1]。由于空化現(xiàn)象機理復(fù)雜，涉及到湍流、兩相流等多個方面，且對離心泵的正常工作具有很大的影響，因此其故障診斷方法一直受到廣大學(xué)者的重視。

目前工業(yè)領(lǐng)域故障診斷主要有3種方法，基于失效物理模型、基于統(tǒng)計可靠性和基于數(shù)據(jù)驅(qū)動的故障診斷方法^[2]。對于實際的復(fù)雜系統(tǒng)而言，直接使用其物理模型進行理論上的分析推導(dǎo)來得到一個可靠的模型十分困難。此時基于數(shù)據(jù)驅(qū)動的故障診斷方法就可以發(fā)揮重要作用。文獻[3]則將壓力信號和診斷信號兩種數(shù)據(jù)進行結(jié)合建模，實現(xiàn)了基于SVM的航空燃油泵故障識別。通過仿真、實驗等多種方法獲取正常或異常工作狀態(tài)下的數(shù)據(jù)，繞過系統(tǒng)的機理，直接對數(shù)據(jù)進行建模，最終形成一個只依靠數(shù)據(jù)的映射模型，來對系統(tǒng)的故障進行診斷。文獻[4]以齒輪泵為研究對象，對其4種常見的故障進行了小波變換+SVM建模，最終建立了一個由4個SVM構(gòu)成的診斷模型，能夠?qū)Πㄕ９ぷ鞯?種狀態(tài)進行有效區(qū)分。文獻[5]對振動信號進行經(jīng)驗?zāi)Ｊ椒纸?，得到了多個固有模態(tài)函數(shù)（IMF），并將IMF中的前6個送入廣義回歸神經(jīng)網(wǎng)絡(luò)中進行故障診斷，最終分類正確率達到了97.5%。

考慮到實際空化原理十分復(fù)雜，對其進行純理論分析顯得十分困難，基于CFD軟件的數(shù)值模擬方法和實驗成為了研究空化現(xiàn)象的主要方法^[6]。本文以旋轉(zhuǎn)區(qū)域的氣體體積分?jǐn)?shù)做為離心泵空化程度的度量并將之作為神經(jīng)網(wǎng)絡(luò)輸出，使用流場中監(jiān)測點的特征以及工作點參數(shù)作為輸入，建立離心泵空化模型，研究神經(jīng)網(wǎng)絡(luò)對離心泵空化程度的診斷效果。

相較于一般的使用神經(jīng)網(wǎng)絡(luò)建模，本文充分考慮到了不同數(shù)據(jù)類型的處理方式，分別對時序型數(shù)據(jù)、數(shù)值型數(shù)據(jù)設(shè)計了網(wǎng)絡(luò)層并對其進行特征融合，以使網(wǎng)絡(luò)建模思路可以在其他更廣泛的問題上得到應(yīng)用。此外，本文對使用的稀疏自編碼器進行了改進，通過修改稀疏性限制函數(shù)來實現(xiàn)更好的正則化，并將自編碼器直接置于大的預(yù)測網(wǎng)絡(luò)之中，使得自編碼器在訓(xùn)練過程中可得到預(yù)測結(jié)果方向的梯度，優(yōu)化了訓(xùn)練過程。

1 離心泵流場數(shù)值仿真

為了獲得離心泵在不同狀態(tài)下的內(nèi)部流場特征，獲取空化故障診斷模型的訓(xùn)練數(shù)據(jù)，本文使用CFD軟件對流場進行了數(shù)值仿真。與實驗相比，數(shù)值仿真不僅簡單易行，且可以獲取離心泵內(nèi)部氣體分布情況，直接對離心泵空化的程度進行分析。

1.1 數(shù)值仿真方法

對流場中的湍流計算選用了標(biāo)準(zhǔn)k-ε模型。該模型由Launder和Spalding提出，計算簡單，魯棒性好，且在大范圍湍流的情況下具有較好的預(yù)測效果。此外，對于空化模型選擇了Signhal全空化模型^[7]。

1.2 物理模型與計算條件

選擇某型單吸離心泵為研究對象，仿真軟件為Pumplinx。如圖1所示，根據(jù)區(qū)域劃分笛卡兒網(wǎng)格，最終網(wǎng)格數(shù)量為35萬。

離心泵額定工作狀態(tài)為：轉(zhuǎn)速9 000r/min，進口壓力0.2MPa，出口流量50 000L/h。為了得到離心泵在不同空化狀態(tài)下的流場特征，通過在不同轉(zhuǎn)速與流量的工況下減小進口壓力，以逐步誘發(fā)離心泵空化。

2 離心泵空化診斷建模

在得到離心泵在不同狀態(tài)下的流場數(shù)據(jù)后，使用人工神經(jīng)網(wǎng)絡(luò)來構(gòu)建泵的空化程度與所選特征的映射模型，并對神經(jīng)網(wǎng)絡(luò)模型在離心泵空化診斷上的有效性進行驗證。

2.1 數(shù)據(jù)集構(gòu)建

基于本文上一節(jié)的仿真方案，選定神經(jīng)網(wǎng)絡(luò)模型輸入數(shù)據(jù)為：流場內(nèi)各部位如進出口、葉片流道、轉(zhuǎn)動軸附近預(yù)先設(shè)置的檢測點壓力、離心泵工作轉(zhuǎn)速、進口壓力和出口流量。

對于離心泵空化程度的表達，選用了轉(zhuǎn)動區(qū)域的氣體體積分?jǐn)?shù)為判定標(biāo)準(zhǔn)。由于體積分?jǐn)?shù)的大小隨著選定區(qū)域的總體積大小變化而變化，而且在一定小范圍內(nèi)出現(xiàn)的波動幾乎不會對流場產(chǎn)生影響，尤其在空化初生階段，往往在流場內(nèi)要積累有一定量的氣體才會對外表現(xiàn)出空化的特征。因此，本文對仿真結(jié)果進行等級劃分，將氣體體積分?jǐn)?shù)劃定為4個等級來表達空化的程度，將問題轉(zhuǎn)化為典型的分類問題。

2.2 神經(jīng)網(wǎng)絡(luò)模型

模型的輸入數(shù)據(jù)有兩種類型：一是監(jiān)測點壓力值，按照軟件設(shè)置的非穩(wěn)態(tài)仿真時間間隔進行采集，得到的是典型的多維時序數(shù)據(jù)，格式為二維張量；二是工作點狀態(tài)，屬于數(shù)值型數(shù)據(jù)，格式為一維張量。

網(wǎng)絡(luò)整體結(jié)構(gòu)如圖2所示，其中，Predict為預(yù)測網(wǎng)絡(luò)，負(fù)責(zé)輸出最終結(jié)果。編碼器Encoder和解碼器Decoder網(wǎng)絡(luò)組成了一個自編碼器網(wǎng)絡(luò)。本文采取了稀疏自編碼器^[8]的思想，不直接限制中間層編碼h的維數(shù)，而是對其實際表達信息量進行限制，并使用一種新的稀疏性限制損失函數(shù)：對中間層編碼h的各維度分量進行了限制。在神經(jīng)網(wǎng)絡(luò)中，一般認(rèn)為某個分量等于0，便認(rèn)為其基本不攜帶信息。對于中間層各分量，如果其低于某個閾值，便將其置為0，而高于閾值的置為1，這就使得中間層的有效輸出維度受到了抑制，實現(xiàn)了稀疏性限制。

2.3 網(wǎng)絡(luò)訓(xùn)練

在工作點中隨機取15個作為測試集，15個作為驗證集，剩余的56個為訓(xùn)練集。由于每個工作點都截取了3段數(shù)據(jù)，因此訓(xùn)練集最終包含有168條數(shù)據(jù)，驗證集和測試集各45條。將訓(xùn)練集數(shù)據(jù)送入未經(jīng)訓(xùn)練的網(wǎng)絡(luò)，使用Adam算法對網(wǎng)絡(luò)進行迭代，并在每輪結(jié)束后再驗證集上進行準(zhǔn)確率測試。

2.4 網(wǎng)絡(luò)訓(xùn)練結(jié)果分析

圖3和圖4為神經(jīng)網(wǎng)絡(luò)在訓(xùn)練集和驗證集上分類損失函數(shù)的變化過程，橫軸為迭代次數(shù)，縱軸為網(wǎng)絡(luò)輸出與實際標(biāo)簽的交叉熵?fù)p失值?？梢钥闯?，損失函數(shù)隨著迭代次數(shù)的增加呈現(xiàn)出了顯著的下降趨勢，而在訓(xùn)練末期，訓(xùn)練集上損失函數(shù)雖仍有下降，但驗證集上卻開始有上升，出現(xiàn)了過擬合趨勢，網(wǎng)絡(luò)過度偏向于訓(xùn)練集的數(shù)據(jù)本身，忽略了輸入數(shù)據(jù)的深層特征，使得網(wǎng)絡(luò)僅僅能適合于訓(xùn)練集數(shù)據(jù)而不具有足夠的泛用性。因此最終選取驗證集上損失最小時刻的模型為最終結(jié)果。

在神經(jīng)網(wǎng)絡(luò)建模流程中，由于測試集在一開始便獨立劃分出去，并未參與模型的訓(xùn)練與調(diào)整，因此可以將其視為全新的數(shù)據(jù)，基于此數(shù)據(jù)的預(yù)測性能也可看做是模型的最終效果。將此前劃分的45條測試集作為網(wǎng)絡(luò)輸入進行預(yù)測，最終獲得的結(jié)果如表1所示，其準(zhǔn)確率為95.6%，達到了預(yù)期效果。

3 結(jié)語

1）人工神經(jīng)網(wǎng)絡(luò)可以有效地建立離心泵流場特征與空化程度之間的對應(yīng)模型，最終準(zhǔn)確率達到了95%以上。因此，使用神經(jīng)網(wǎng)絡(luò)來診斷離心泵空化情況的方法是有效的，在未來的仿真模擬或者實機運行中可以發(fā)揮作用。

2）在網(wǎng)絡(luò)中嵌入了稀疏自編碼器，并對稀疏性損失函數(shù)進行優(yōu)化，能夠有效地對帶有噪聲的數(shù)據(jù)進行降噪，并且在實現(xiàn)信息壓縮的情況下防止過擬合。

3）針對實際數(shù)據(jù)中最常見的兩種類型：時間序列和數(shù)值型特征，設(shè)計了不同的網(wǎng)絡(luò)類型對其進行處理，并且使用concatenate層將二者的信息進行了融合。對于沒有出現(xiàn)的分類型特征，也可以參考數(shù)值型特征使用全連接或其他網(wǎng)絡(luò)層進行處理，基本可包含所有可能出現(xiàn)的數(shù)據(jù)類型。通過多種類輸入、多網(wǎng)絡(luò)融合的神經(jīng)網(wǎng)絡(luò)設(shè)計思路，使該網(wǎng)絡(luò)以及其建模思想獲得更廣的應(yīng)用場景。

參考文獻：

[1] LU J X，YUAN S Q，ZHOU B L，et al. Research on the characteristics under the condition of asymmetric cavitation in a centrifugal pump[C]//2014 ISFMFE - 6th International Symposium on Fluid Machinery and Fluid Engineering. Wuhan. London：IET，2015：1-8.

[2] 李晗，蕭德云. 基于數(shù)據(jù)驅(qū)動的故障診斷方法綜述[J]. 控制與決策，2011，26（1）：1-9，16.

[3] 陳俊柏，劉勇智，陳勇，等. 基于EMD和SVM的機載燃油泵故障診斷研究[J]. 北京航空航天大學(xué)學(xué)報，2021，47（8）：1687-1696.

[4] 毋文峰，蘇勛家，陳小虎，等. 基于小波包和支持向量機的液壓泵故障診斷[J]. 機床與液壓，2011，39（9）：146-147，154.

[5] AZIZI R，ATTARAN B，HAJNAYEB A，et al. Improving accuracy of cavitation severity detection in centrifugal pumps using a hybrid feature selection technique[J]. Measurement，2017，108：9-17.

[6] 程文潔，張飛宇. 低溫泵內(nèi)空化流動演化規(guī)律的數(shù)值仿真[J]. 低溫與超導(dǎo)，2021，49（8）：78-81，92.

[7] SINGHAL A K，ATHAVALE M M，LI H Y，et al. Mathematical basis and validation of the full cavitation model[J]. Journal of Fluids Engineering，2002，124（3）：617-624.

[8] NG A. Sparse autoencoder[J]. CS294A Lecture Notes，2011， 72（1）：1-19.

收稿日期：2023-01-30