章修龍,劉貴杰,寧東紅
(中國海洋大學(xué) 機(jī)電工程系,山東 青島 266100)
水下采油樹結(jié)構(gòu)復(fù)雜,在高溫、高壓的惡劣海洋環(huán)境中作業(yè)容易發(fā)生故障。深海設(shè)備檢修維護(hù)非常困難,維修成本極高[1],一旦在生產(chǎn)中出現(xiàn)故障,會給日常生產(chǎn)造成很大的損失。
閥門作為水下采油樹中的關(guān)鍵裝置,工作時處于高溫、高壓的惡劣環(huán)境中[2],受到介質(zhì)的長期沖刷;且高溫閥門關(guān)閉后會迅速冷卻,使密封面出現(xiàn)漏孔或裂紋,易導(dǎo)致泄漏,影響水下采油樹的正常工作[3]。
閥門受損后產(chǎn)生的泄漏故障,會影響水下油氣開采設(shè)備的正常工作,造成重大的經(jīng)濟(jì)損失。因此,為保障水下采油樹健康、高效地正常運行,對水下采油樹閥門泄漏故障進(jìn)行診斷具有重要的意義。
不少學(xué)者已經(jīng)對閥門進(jìn)行了故障診斷的研究,并提出了不少方法。寧方立等人[4]為了檢測輸氣管道閥門泄漏,提出了基于改進(jìn)卷積神經(jīng)網(wǎng)絡(luò)(convolutional neural network, CNN)的閥門泄漏超聲信號識別方法,提高了閥門故障的識別性能。楊海榮等人[5]提出了一種基于主元分析和支持向量多分類器的故障診斷方法,實現(xiàn)了對閥門故障快速、準(zhǔn)確檢測與診斷目的。
在水下環(huán)境中,為克服接觸式故障診斷方法中存在的安裝和維護(hù)困難問題,可采用聲學(xué)傳感器采集的水聲信號對故障進(jìn)行診斷。在該研究領(lǐng)域中,XIE Ying-chun等人[6]利用數(shù)字水聽器收集的聲學(xué)泄漏信號檢測了海底生產(chǎn)系統(tǒng)管道泄漏,并利用經(jīng)驗?zāi)B(tài)分解和希爾伯特-黃變換獲得的信號,準(zhǔn)確識別了管道泄漏。
然而,聲音在水下的傳輸過程中,由于吸收了水下環(huán)境背景噪聲,聲傳感器吸收了大量的噪聲成分,會對最終的診斷結(jié)果造成重大影響。因此,需要尋找一種合適的降噪方法,從收集的原始聲信號中去除噪聲成分。
傳統(tǒng)的信號處理方法通常是使用基于經(jīng)驗的方法對信號進(jìn)行預(yù)處理。例如,ZHOU Peng等人[7]利用集合經(jīng)驗?zāi)B(tài)分解(ensembleempirical mode decomposition, EEMD)對信號進(jìn)行降噪預(yù)處理,獲得了清晰的故障特征信號,提高了故障診斷的準(zhǔn)確性。JIANG Wei-xiong等人[8]通過使用自適應(yīng)噪聲的集成經(jīng)驗?zāi)B(tài)分解(complementary EEMD with adaptive noise, CEEMDAN)方法,實現(xiàn)了從原始聲壓信號中去除部分噪聲的目的。
然而,當(dāng)噪聲信號與水下聲音信號處于同一頻帶范圍時,這些方法無法取得有效的降噪效果。當(dāng)噪聲信號與水聲信號存在頻率混疊(即二者有部分頻率處于相同的范圍內(nèi))時,雖然采用基于時頻域的降噪方法[9],可以實現(xiàn)去除噪聲的目的(比如,小波變換[10]和S-變換[11,12]等),但是基于這些方法去除含有混疊噪聲信號時,對閾值函數(shù)的選擇要求極高,使其在實際應(yīng)用中受到限制。
近年來,基于深度學(xué)習(xí)的降噪方法受到了業(yè)界的廣泛關(guān)注,同時該類方法也確實可用于學(xué)習(xí)水聲數(shù)據(jù)的特征[13]。但由于復(fù)雜的海洋環(huán)境和信號傳輸過程中的水聲噪聲,需要找到一種自適應(yīng)智能降噪方法,以最大限度地減少數(shù)字水聽器收集的水聲數(shù)據(jù)中的噪聲成分,提高故障診斷的準(zhǔn)確性。
針對上述問題,筆者通過研究萊特希爾聲比擬理論、降噪的基本理論及方法,結(jié)合數(shù)值分析與實驗驗證,提出一種基于深度神經(jīng)網(wǎng)絡(luò)(deep neural network, DNN)的閥門故障診斷方法,以降低水聲信號數(shù)據(jù)中噪聲成分,提高閥門故障診斷的準(zhǔn)確率。
閥門泄露主要由密封結(jié)構(gòu)被破壞引起,而閥門密封面未關(guān)嚴(yán)導(dǎo)致縫隙泄漏,以及閥門出現(xiàn)漏孔導(dǎo)致泄漏,都是較為常見的泄漏方式。以縫隙泄漏和漏孔泄漏為例,流體在經(jīng)過縫隙或漏孔高速噴出后,進(jìn)入到一個足夠大的空間內(nèi),由于不再受邊界的限制而進(jìn)一步擴(kuò)散流動,該現(xiàn)象稱為射流。
液體從泄漏口高速噴射出來,所產(chǎn)生的噴流噪聲可用軸對稱自由射流表示,如圖1所示。
圖1中,閥門漏孔泄漏口的直徑為D,當(dāng)液體從泄漏口高速噴射出去時(射流的分布可分為混合區(qū)、過渡區(qū)和充分發(fā)展區(qū)),混合區(qū)的液體流速最快,只有少量湍流的產(chǎn)生,其內(nèi)部的速度可以看作與泄漏口的流速近似相同;過渡區(qū)中有大量的湍流產(chǎn)生,該部分區(qū)域的液體流速隨著距離的增加而減小,混合區(qū)與過渡區(qū)之間的湍流也是射流噪聲生成的主要原因;過渡區(qū)之后是寬廣的充分發(fā)展區(qū),該區(qū)域的流速相對較小,且隨著距離的增加,流速呈遞減趨勢。
閥門泄漏發(fā)聲機(jī)理實質(zhì)上是屬于射流產(chǎn)生的噪聲,而射流噪聲通常是基于萊特希爾聲比擬方法來處理的。萊特希爾(Lighthill)[14,15]將聲場分為聲源場(近場)和遠(yuǎn)場兩部分。聲源場處,高速噴射的流體產(chǎn)生湍流,湍流中的部分能量轉(zhuǎn)換成聲能,聲能經(jīng)海水介質(zhì)傳播到遠(yuǎn)場。
為方便計算,筆者將聲源點y作為坐標(biāo)原點建立笛卡爾坐標(biāo)系,如圖2所示。
圖2中,r表示近場聲源與遠(yuǎn)場之間的距離,遠(yuǎn)場與坐標(biāo)軸i之間的夾角記為α,坐標(biāo)軸j與坐標(biāo)軸k之間的夾角記為β。
液體噴流產(chǎn)生的聲場遵循有源聲場的波動方程,萊特希爾聲比擬方法的基本理論如下:
根據(jù)納維-斯托克斯方程導(dǎo)出的基本公式,考慮靜止介質(zhì)及有源聲場的情況,假設(shè)介質(zhì)密度為ρ,聲速為c,則其連續(xù)方程可表示為:
(1)
式中:vi—質(zhì)點速度矢量。
vi的動量方程可以表示為:
(2)
式中:pij—微元體表面上的應(yīng)力張量。
結(jié)合式(1,2)可得:
(3)
式(1,3)消去ρvi,可得萊特希爾方程:
(4)
令p=c2ρ,則式(4)可表示為:
(5)
從式(4,5)可以看出:流體右邊存在聲源,表明其存在傳播現(xiàn)象,可以計算出遠(yuǎn)場的聲壓。
由于液體噴流產(chǎn)生的聲場遵循有源聲場的波動方程,大量的研究結(jié)果表明,射流產(chǎn)生的聲源屬于四極子聲源,該類型聲源在結(jié)構(gòu)上完全對稱。四極子聲源中兩個偶極子之間距離為d,且存在應(yīng)力pij,即存在應(yīng)力張量Tij。
因此,四級子聲源聲波方程的解為:
(6)
由于r=|x-y|,式(6)經(jīng)過微分計算后可得:
(7)
(8)
式(8)中,括號內(nèi)第一項代表遠(yuǎn)場的聲壓,即:
(9)
由式(9)可以看出,遠(yuǎn)場聲壓的大小與近場至遠(yuǎn)場之間的距離成反比,與頻率成正相關(guān)關(guān)系。
1.2.1 降噪基本理論
數(shù)字水聽器采集的信號類型為遠(yuǎn)場的聲壓信號,從水下環(huán)境中收集的原始聲壓信號與環(huán)境背景噪聲信號混合在一起。因此,為了獲得準(zhǔn)確的故障診斷結(jié)果,需要對采集到的原始聲壓信號進(jìn)行去噪。
將時域中帶有噪聲的閥門泄漏信號轉(zhuǎn)換為時頻域后,在時頻域中,筆者將帶有噪聲的閥門泄漏信號標(biāo)記為:
Y(t,f)=X(t,f)+N(t,f)
(10)
式中:X(t,f)—閥門泄漏產(chǎn)生的不含噪聲成分的時頻域信號,也是該實驗降噪處理后希望得到的目標(biāo)信號;N(t,f)—復(fù)雜的環(huán)境背景噪聲的時頻域信號。
(11)
1.2.2 DNN模型
受到卷積降噪自動編碼器的啟發(fā),基于降噪的基本理論,筆者建立用于水聲信號降噪處理的DNN模型結(jié)構(gòu),如圖3所示。
該神經(jīng)網(wǎng)絡(luò)設(shè)計為由一系列的二維卷積層、批量歸一化層和整流線性單元(rectified linear unit, ReLU)層組成。輸入的幅值頻譜,通過這些二維卷積層進(jìn)行處理和轉(zhuǎn)換。通過使用2×2的步幅,減少了網(wǎng)絡(luò)反卷積層的特征空間,但卷積濾波器尺寸保持在3×3。卷積處理和下采樣操作可以實現(xiàn)網(wǎng)絡(luò)的特征提取功能,并壓縮特征矩陣,使下一次卷積操作的參數(shù)數(shù)量最小化。
在該神經(jīng)網(wǎng)絡(luò)的反卷積層中,這些層被用來重建特征,并對降噪信號的輸出幅值譜生成高維非線性映射。神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)在最后一層中使用一個回歸函數(shù),用來降低其輸出結(jié)果與輸入目標(biāo)之間的均方誤差。因為,隨著設(shè)計的神經(jīng)網(wǎng)絡(luò)層數(shù)的增多、深度的增加,神經(jīng)網(wǎng)絡(luò)的準(zhǔn)確度會出現(xiàn)飽和現(xiàn)象,可能還會出現(xiàn)梯度消失的問題。因此,該網(wǎng)絡(luò)使用跨層連接將卷積層與反卷積層連接起來,以實現(xiàn)加快訓(xùn)練的收斂性,達(dá)到高質(zhì)量局部最優(yōu)目的。
模型的詳細(xì)參數(shù)如表1所示。
1.2.3 原始水聲信號數(shù)據(jù)降噪處理流程
筆者提出的基于DNN的降噪方法的基本流程為:
(1)對其進(jìn)行數(shù)據(jù)預(yù)處理。通過使用短時傅里葉變換(short-time Fourier transform, STFT),將時域信號轉(zhuǎn)換為時頻域信號,使其成為神經(jīng)網(wǎng)絡(luò)可以訓(xùn)練學(xué)習(xí)的數(shù)據(jù)類型;
(2)將不含噪聲的閥門泄漏水聲信號和含噪聲的閥門泄漏水聲信號的幅值頻譜,輸入神經(jīng)網(wǎng)絡(luò)中,通過神經(jīng)網(wǎng)絡(luò)對含噪聲與不含噪聲的閥門泄漏水聲信號進(jìn)行學(xué)習(xí)訓(xùn)練,對含噪聲的閥門泄漏水聲信號進(jìn)行降噪處理;
表1 改進(jìn)的DNN模型的詳細(xì)參數(shù)
(3)將降噪后的信號與噪聲信號的幅值時頻譜輸出;
(4)通過使用短時傅里葉逆變換,將時頻域信號轉(zhuǎn)換為時域信號,得到筆者所需要的降噪后的閥門泄漏水聲信號數(shù)據(jù)。
在其訓(xùn)練和測試過程中,模型輸入的是含噪聲的閥門泄漏水聲信號幅值頻譜,和不含噪聲的閥門泄漏水聲信號幅值頻譜,輸出的是降噪后的閥門泄漏水聲信號的幅值頻譜,與降噪后得到的噪聲信號的幅值頻譜。
基于DNN降噪方法流程圖如圖4所示。
水下采油樹閥門泄漏的故障診斷方法主要有3個步驟,分別為原始數(shù)據(jù)獲取、原始數(shù)據(jù)預(yù)處理及降噪、故障診斷。
基于深度卷積神經(jīng)網(wǎng)絡(luò)進(jìn)行水聲信號降噪處理的,水下采油樹閥門泄漏故障診斷模型的流程,如圖5所示。
筆者使用COMSOL軟件對聲場進(jìn)行仿真處理,以確定數(shù)字水聽器的位置,并通過數(shù)字水聽器獲取原始水聲數(shù)據(jù),然后將采集好的原始水聲數(shù)據(jù)按照故障類型標(biāo)記進(jìn)行分類,對其使用STFT將時域信號轉(zhuǎn)變?yōu)闀r頻域信號,再利用設(shè)計訓(xùn)練好的水聲信號降噪改進(jìn)的DNN處理該時頻域信號,信號經(jīng)過降噪處理后,最后將其輸入至用于故障分類的卷積神經(jīng)網(wǎng)絡(luò)中,診斷其故障類型。
筆者采用數(shù)字水聲器采集水下采油樹閥門泄漏產(chǎn)生的水聲信號,還應(yīng)用Lighthill理論、水聲傳播理論和COMSOL軟件進(jìn)行數(shù)值模擬分析,獲得了數(shù)字水聽器的最優(yōu)位置。
在測試水下采油樹閥門泄漏的遠(yuǎn)場聲壓信號時,以圖2所建立的坐標(biāo)系為基礎(chǔ),將坐標(biāo)系中的x用水聽器代替,以確定距離和方位(圖中,α和β分別為俯仰角補角和方位角)。
數(shù)字水聽器采集到的原始水聲信號數(shù)據(jù)中含有大量的環(huán)境背景噪聲成分,信號的信噪比較高,對故障診斷的結(jié)果產(chǎn)生很大的影響。因此,需要在進(jìn)行故障診斷之前對原始水聲信號數(shù)據(jù)進(jìn)行處理。
由于數(shù)字水聽器采集到的原始水聲信號為時域信號,無法對其進(jìn)行頻域分析,因此,筆者采用STFT將時域的原始閥門泄漏水聲信號轉(zhuǎn)換為時頻域的信號,便于對其進(jìn)行頻域分析,再將轉(zhuǎn)換后的時頻域信號作為DNN的輸入數(shù)據(jù),進(jìn)行降噪處理。
采用基于神經(jīng)網(wǎng)絡(luò)的故障診斷方法進(jìn)行診斷時,需要對處理后的數(shù)據(jù)進(jìn)行特征提取,再分類標(biāo)記,然后對其進(jìn)行故障診斷。而卷積神經(jīng)網(wǎng)絡(luò)在卷積處理過程中,卷積層可以自主進(jìn)行特征提取,因此,對神經(jīng)網(wǎng)絡(luò)進(jìn)行一定的改進(jìn)有望實現(xiàn)該目的。
受一維殘差卷積神經(jīng)網(wǎng)絡(luò)的啟發(fā)[16],筆者設(shè)計了一種用于故障診斷的卷積神經(jīng)網(wǎng)絡(luò)(其不需要大量的標(biāo)記圖像,可以無監(jiān)督、自適應(yīng)地提取降噪后的閥門泄漏水聲信號特征,實現(xiàn)故障診斷目標(biāo))。
卷積神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)如圖6所示。
該實驗在中國海洋大學(xué)工程訓(xùn)練實驗中心的水池中進(jìn)行。實驗水池由鋼架、高強(qiáng)度防水布組成,其長、寬和深度分別為2.00 m、1.00 m和0.98 m。該實驗用的閥門類型為#DN10,以液體侵蝕和機(jī)械故障引起的閥門泄漏作為典型的故障類型。
實驗設(shè)計了4個不同漏孔直徑的閥門作為實驗組,閥門的漏孔直徑分別為d1=0 mm,d2=3.00 mm,d3=4.00 mm,d4=5.00 mm。
實驗中,閥門的出水壓力由一個恒壓水箱提供,它可以調(diào)節(jié)出水口的壓力。由于水箱的壓力限制,通過設(shè)計使用更小內(nèi)徑的管道來提高閥門出水口的壓力,將閥門前部的水壓最高提高到0.40 MPa,因此,該實驗將閥門泄漏口之前的水壓設(shè)置為0.10 MPa,0.20 MPa,0.30 MPa。
實驗中設(shè)計的9種不同的工況類型如表2所示。
表2 不同壓力及孔徑下的故障類型
為了模擬真正的實驗環(huán)境,筆者設(shè)計了流動面積大小為長度2.00 m、寬度1.00 m、深度0.98 m的矩形域,與實驗的實際尺寸相匹配,分別表示向量坐標(biāo)系中x軸、y軸、z軸的坐標(biāo)。
閥門模型安裝在流場底部,內(nèi)外流體密度分別為1 000 kg/m3。z方向的上邊界由聲音軟邊界條件控制,用于模擬池中空氣介質(zhì)與水介質(zhì)之間的界面。
筆者將其他邊界作為聲音硬邊界條件,以最小化池的混響對聲壓的影響。
多體動力學(xué)和聲學(xué)響應(yīng)分析結(jié)果如圖7所示。
圖7(a)為0.30 MPa故障條件下的聲壓場分布。由于水和空氣接觸表面的散射效應(yīng)和池體邊界的反射效應(yīng),筆者發(fā)現(xiàn)聲場分布圖中的聲壓在泄漏孔處有很強(qiáng)的波動;聲壓沿y軸方向迅速下降,沿x軸方向逐漸減弱,聲壓峰值出現(xiàn)在混合區(qū)和湍流區(qū)之間;
圖7(b)為相同故障條件下流場的速度分布。液體泄漏產(chǎn)生的射流沿x軸波動最大,流速在泄漏孔附近達(dá)到峰值后迅速下降。該現(xiàn)象表明,在混合區(qū)和湍流區(qū)之間的邊界附近是聲場的最強(qiáng)值區(qū)域。因此,該處是安裝數(shù)字水聽器,以提高原始聲壓信號信噪比的最佳位置。
筆者將閥門固定在實驗池底部,閥門泄漏處朝向數(shù)字水聽器,數(shù)字水聽器的安裝方向與水平面垂直,角度為90°,安裝位置處于混合區(qū)和湍流區(qū)之間的邊界附近。
實驗現(xiàn)場如圖8所示。
實驗現(xiàn)場可分為實驗平臺數(shù)據(jù)采集區(qū)及實驗區(qū)。當(dāng)閥門泄漏發(fā)生時,泄漏產(chǎn)生的水聲信號在水中傳播之后,由遠(yuǎn)場的數(shù)字水聽器進(jìn)行收集,并通過線纜傳輸至電腦。傳輸過程中,使用聲學(xué)信號放大器對水聲信號進(jìn)行放大處理,增強(qiáng)信號的信噪比。
水聲信號數(shù)據(jù)由實驗采集得到。實驗中,使用構(gòu)建好的改進(jìn)DNN來降低原始聲壓數(shù)據(jù)中的噪聲成分,將采集到的時域聲壓信號轉(zhuǎn)換為時頻圖,并利用轉(zhuǎn)換好的時頻圖來訓(xùn)練神經(jīng)網(wǎng)絡(luò),并測試其性能。
筆者將采集到的信號進(jìn)行分割,得到2 400個不含噪聲成分的閥門泄漏水聲頻譜,作為泄漏信號樣本(該部分信號具有高信噪比,需要10個不同類型的噪聲數(shù)據(jù)作為噪聲信號樣本);在2 400個泄漏信號樣本中,筆者隨機(jī)選取80%作為訓(xùn)練樣本,10%作為驗證樣本,10%作為測試樣本。
在訓(xùn)練過程中,為了確保神經(jīng)網(wǎng)絡(luò)能夠識別泄漏信號和噪聲信號的特征,筆者對訓(xùn)練泄漏信號樣本進(jìn)行反復(fù)迭代,每次迭代中選擇一個噪聲樣本,并將其疊加到泄漏信號樣本上,以降低泄漏信號的信噪比,通過連續(xù)訓(xùn)練,積累訓(xùn)練數(shù)據(jù)集。
在預(yù)測過程中,計算和記錄歸一化均方誤差。筆者將泄漏信號的振幅和含有噪聲的泄漏信號的振幅輸入神經(jīng)網(wǎng)絡(luò)中,進(jìn)行獨立學(xué)習(xí),該神經(jīng)網(wǎng)絡(luò)的預(yù)測目標(biāo)為泄漏信號的振幅(即期望的輸出結(jié)果)。在驗證過程中,利用驗證樣本調(diào)整網(wǎng)絡(luò)參數(shù),防止神經(jīng)網(wǎng)絡(luò)出現(xiàn)過擬合現(xiàn)象,以實現(xiàn)更好的降噪性能。
利用測試樣本去分析網(wǎng)絡(luò)的降噪性能,在實驗中,筆者用采集到的由水下設(shè)備機(jī)械振動引起的噪聲(簡稱振動噪聲)波形來測試該網(wǎng)絡(luò),其實驗結(jié)果如圖9所示。
圖9中(1)圖表示的是無噪聲的原始閥門泄漏聲壓信號的頻譜圖及波形圖。對泄漏信號樣本進(jìn)行迭代的過程中,將該振動噪聲樣本疊加到泄漏信號樣本上,得到如圖9中(2)圖所示的頻譜圖及波形圖。
含有振動噪聲的復(fù)合信號經(jīng)過神經(jīng)網(wǎng)絡(luò)降噪處理后,降噪后的水聲信號其結(jié)果如圖9中(3)圖所示,降噪后的噪聲信號其結(jié)果如圖9中(4)圖所示。
通過對比降噪后的閥門泄漏水聲信號與降噪后的振動噪聲可以發(fā)現(xiàn):降噪后的閥門泄漏水聲信號,與添加噪聲之前的閥門泄漏水聲數(shù)據(jù)的波形圖類似,由于振動噪聲波形的幅值、周期性、頻率等特征明顯,有利于提取噪聲特征。
該結(jié)果表明:該神經(jīng)網(wǎng)絡(luò)可以成功地將含有振動噪聲的復(fù)合信號分解為振動噪聲和降噪信號;該降噪方法能有效地降低復(fù)合信號中的噪聲成分,較好地恢復(fù)原始信號。因此,該方法能有效地從復(fù)合信號中分離出振動類的噪聲。
為了進(jìn)一步驗證該降噪方法的有效性,筆者將基于深度卷積神經(jīng)網(wǎng)絡(luò)的降噪方法與其他降噪方法進(jìn)行對比分析[17]。
筆者選擇全卷積神經(jīng)網(wǎng)絡(luò)(fully convolutional neutral network, FCNN)[18]與常規(guī)的卷積神經(jīng)網(wǎng)絡(luò)[19]作為對比方法組進(jìn)行比較測試。其中,全卷積神經(jīng)網(wǎng)絡(luò)為所有層都是卷積層的神經(jīng)網(wǎng)絡(luò),比卷積神經(jīng)網(wǎng)絡(luò)計算速度更快、效率更高,有助于網(wǎng)絡(luò)性能的提升。
筆者從收集到的具有環(huán)境背景噪聲的聲壓信號中,分別選擇100~600個數(shù)據(jù)樣本進(jìn)行實驗,每次遞增100個數(shù)據(jù)樣本。
為了避免樣本選擇對降噪結(jié)果的影響,筆者采用控制變量法,對相同條件下的不同方法都使用相同的樣本;此外,增加樣本數(shù)量,同時保證神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)和參數(shù)保持不變。
為了更好地衡量該方法降噪效果的優(yōu)劣,筆者還設(shè)置了一組對照組,將其與上述3種方法進(jìn)行對比,即通過設(shè)置含有環(huán)境背景噪聲的原始閥門泄漏聲壓信號(簡稱含噪聲的信號)作為對照組,將降噪后得到的信號其歸一化均方誤差值變化進(jìn)行對比,結(jié)果如圖10所示。
由圖10可知:隨著樣本數(shù)量的增加,未經(jīng)處理的含噪聲的信號歸一化均方誤差值穩(wěn)定在0.50附近;
經(jīng)過降噪處理后,基于FCNN的降噪方法所得到的信號歸一化均方誤差值由0.499 2最低降至0.011 2,降幅達(dá)97.76%;基于CNN的降噪方法所得到的信號歸一化均方誤差值由0.499 2最低降至0.015 6,降幅達(dá)96.88%;筆者所提出的降噪方法所得到的信號歸一化均方誤差值由0.499 2最低降至0.011 0,降幅達(dá)97.80%。
歸一化均方誤差值在該實驗中用于衡量降噪后信號與目標(biāo)信號之間的相似度,數(shù)值越小代表二者之間的相似度越高。
因此,在歸一化均方誤差值這一評價指標(biāo)中,筆者提出的降噪方法的表現(xiàn)優(yōu)于其他兩種降噪方法,這也證明經(jīng)過降噪后的信號波形更接近目標(biāo)信號,具有更好的精度和更強(qiáng)的魯棒性。
通過對9種不同工況下的含噪聲的原始水聲信號進(jìn)行降噪處理,筆者得到了9種降噪后不同類型的水聲故障數(shù)據(jù),它們分別對應(yīng)于故障1至故障9,并將其存儲在一個數(shù)據(jù)集中(該數(shù)據(jù)集包含4 500個數(shù)據(jù)組。每種故障類型包含500個數(shù)據(jù)組(每個數(shù)據(jù)組包含1 024個數(shù)據(jù)點),隨機(jī)選取其中的80%為訓(xùn)練樣本,10%為測試樣本,10%為驗證樣品);
然后,筆者將3 600個訓(xùn)練樣本納入CNN模型進(jìn)行訓(xùn)練,用450個測試樣本評估模型的泛化能力,用450個驗證樣本對模型的參數(shù)進(jìn)行調(diào)整;
經(jīng)過參數(shù)調(diào)整處理后,得到了網(wǎng)絡(luò)模型,模型中E—卷積層+批量歸一化層+整流線性單元層;F—卷積層+整流線性單元層。
模型的詳細(xì)參數(shù)如表3所示。
表3 卷積神經(jīng)網(wǎng)絡(luò)模型的詳細(xì)參數(shù)
對含環(huán)境背景噪聲的原始閥門泄漏水聲信號進(jìn)行降噪處理后,得到故障診斷所需要的數(shù)據(jù)集,數(shù)據(jù)集被分成了訓(xùn)練集、測試集和驗證集[20]。
分組后,筆者將這些數(shù)據(jù)集輸入至搭建好的卷積神經(jīng)網(wǎng)絡(luò)模型中,按照最小批量尺寸為128,學(xué)習(xí)率為0.001,下降周期為5,下降因子為0.1,迭代次數(shù)為300進(jìn)行訓(xùn)練。
經(jīng)過300次迭代之后,筆者從3 600個訓(xùn)練樣本中隨機(jī)選取250個樣本,從450個測試樣本中隨機(jī)選取50個樣本,得到最終的診斷結(jié)果,如圖11所示。
從圖11可以看出:訓(xùn)練和測試的卷積神經(jīng)網(wǎng)絡(luò)故障診斷結(jié)果均顯示出最優(yōu)診斷率為100%,訓(xùn)練樣本的診斷準(zhǔn)確率均高于測試樣本。
最后,筆者采用基于改進(jìn)的DNN用于降噪處理的故障診斷方法,其診斷準(zhǔn)確率為98.89%。結(jié)合故障診斷結(jié)果及各故障類型的信號特點,發(fā)現(xiàn)泄漏壓力小、泄漏孔徑小的故障工況易被誤診。
筆者將用于降噪處理改進(jìn)的DNN故障診斷結(jié)果,與其他兩種對照組的降噪處理方法的故障診斷結(jié)果進(jìn)行對比,結(jié)果如表4所示。
表4 診斷結(jié)果對比
從表4中可以看出:經(jīng)過改進(jìn)的DNN處理后的數(shù)據(jù)輸入至同一卷積神經(jīng)網(wǎng)絡(luò)進(jìn)行故障診斷,其診斷率比未經(jīng)過處理的信號提升了46.42%;相比于其他兩種方法,其分別提升了5.17%和2.55%。
筆者通過研究閥門泄漏基本理論,建立了閥門泄漏模型,求解出了遠(yuǎn)場聲壓;基于降噪的基本理論和方法,建立了DNN模型用于水聲降噪處理,對降噪結(jié)果進(jìn)行了有效性驗證;還將處理后的信號輸入至卷積神經(jīng)網(wǎng)絡(luò)中進(jìn)行故障診斷。
筆者結(jié)合數(shù)值模擬、實驗驗證得出了如下結(jié)論:
(1)綜合使用閥門泄漏理論,包括軸對稱自由射流模型、萊特希爾聲比擬理論和四極子聲源遠(yuǎn)場聲壓求解算法,建立了閥門泄漏射流聲場模型,并利用COMSOL數(shù)值仿真軟件,完成了對數(shù)字水聽器位置的確定,為聲學(xué)實驗開展奠定了基礎(chǔ);
(2)受卷積降噪自動編碼器的啟發(fā),基于降噪的基本理論,筆者設(shè)計并使用DNN對輸入的信號進(jìn)行了降噪處理,實驗證明,該網(wǎng)絡(luò)可以有效地去除噪聲成分;除此之外,利用CNN作為故障分類器,實現(xiàn)了對水下采油樹閥門泄漏故障的準(zhǔn)確診斷,準(zhǔn)確率為98.89%,相比于其他兩種方法,該方法的故障診斷準(zhǔn)確率更高。
在之后的研究中,筆者擬進(jìn)一步探索多噪聲環(huán)境背景下該降噪模型的降噪效果;并開展海試,獲取真實海況下的原始聲學(xué)數(shù)據(jù);并基于該數(shù)據(jù)遠(yuǎn)程診斷水下閥門的故障類型,優(yōu)化神經(jīng)網(wǎng)絡(luò)模型結(jié)構(gòu)參數(shù),進(jìn)一步驗證該降噪模型在實際工況環(huán)境下的有效性和正確性。