李黎明，張微之，季子良，景 嚴(yán)，龐國(guó)富，王 力，宋三明

（1.中國(guó)科學(xué)院沈陽(yáng)自動(dòng)化研究所 機(jī)器人學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，遼寧 沈陽(yáng) 110016；2.中國(guó)科學(xué)院機(jī)器人與智能制造創(chuàng)新研究院，遼寧 沈陽(yáng) 110169；3.中國(guó)科學(xué)院大學(xué)，北京 100049；4.沈陽(yáng)工業(yè)大學(xué)人工智能學(xué)院，遼寧 沈陽(yáng) 110870；5.沈陽(yáng)化工大學(xué)信息工程學(xué)院，遼寧 沈陽(yáng) 110142）

0 引言

水下機(jī)動(dòng)目標(biāo)在航行時(shí)會(huì)向水中輻射噪聲。因此，可以通過水聽器接收到的噪聲信號(hào)對(duì)航行器的槳葉數(shù)、螺旋槳轉(zhuǎn)速、噸位數(shù)以及推進(jìn)器類型等運(yùn)動(dòng)和物理特性等進(jìn)行監(jiān)測(cè)[1]。通常情況下，采用艦船搭載舷側(cè)陣或拖曳陣采集目標(biāo)噪聲[2-3]。隨著水下無(wú)人系統(tǒng)技術(shù)的日漸成熟，采用水下機(jī)器人如AUV搭載水聽器（陣列）對(duì)目標(biāo)進(jìn)行探測(cè)開始成為可能（圖1）。

但是，基于AUV的移動(dòng)探測(cè)任務(wù)至少面臨如下2方面的挑戰(zhàn)：1）由于本體尺寸的原因，AUV搭載的舷側(cè)陣尺寸和孔徑要小得多；2）由于能源和動(dòng)力的限制，拖曳陣的陣元數(shù)規(guī)模和孔徑不僅相對(duì)較小，而且更加靠近本體，導(dǎo)致目標(biāo)噪聲淹沒在本體噪聲中。因此，AUV移動(dòng)探測(cè)的一個(gè)核心任務(wù)是如何消除或者降低本體噪聲的影響，提高目標(biāo)探測(cè)和識(shí)別的效率。

針對(duì)移動(dòng)平臺(tái)的自噪聲去除問題，現(xiàn)有的研究主要基于自適應(yīng)濾波和帶通濾波的思想展開。文獻(xiàn)[3]、文獻(xiàn)[4]分別提出了基于自適應(yīng)噪聲消除和自適應(yīng)波束形成的算法，用于移動(dòng)AUV拖曳平臺(tái)和拖曳船的本體噪聲消除。文獻(xiàn)[5]根據(jù)平臺(tái)噪聲和目標(biāo)噪聲之間的頻譜差異進(jìn)行帶通濾波，去除本體噪聲，提高目標(biāo)探測(cè)能力。文獻(xiàn)[6]針對(duì)矢量舷側(cè)陣，根據(jù)本體噪聲的方向選擇性設(shè)置通、阻帶得到空域矩陣，通過空間濾波抑制本體噪聲的干擾。這些研究均將本體噪聲視為噪聲成分，通過抵消或抑制策略消除或減弱其影響。但是，從聲源的角度看，本體噪聲也是一個(gè)聲源，因此自噪聲抑制問題也可以視為聲源分離（一般稱為盲源分離）任務(wù)。盲源分離算法在水聲信號(hào)處理中主要用于多聲源場(chǎng)景下的分離任務(wù)。Smaragdis等提出了一種基于短時(shí)傅里葉譜的Infomax盲源分離算法，在2個(gè)聲源卷積混合條件下實(shí)現(xiàn)了聲源分離[7]。Gaeta 等對(duì)噪聲信號(hào)進(jìn)行局部白化處理，提出了基于神經(jīng)網(wǎng)絡(luò)的盲源分離算法，實(shí)現(xiàn)了淺水信道盲源分離[8]。在考慮多水聲信道多途效應(yīng)的基礎(chǔ)上，文獻(xiàn)[9]開展了多輸入多輸出通道模型的盲源分離研究，在頻域上采用最小化二階統(tǒng)計(jì)量對(duì)源信號(hào)進(jìn)行估計(jì)。截止目前，水聲信號(hào)盲源分離算法多基于模型驅(qū)動(dòng)，根據(jù)混合噪聲模型設(shè)定合適的分離準(zhǔn)則，通過優(yōu)化求解分離矩陣獲得單一聲源數(shù)據(jù)。且現(xiàn)有研究多基于陣列信號(hào)，基于單水聽器信號(hào)的盲源分離較為少見。本文將在單水聽器和聲源分離的框架下開展本體噪聲去除研究。

濾波法需要目標(biāo)噪聲與干擾聲源和背景噪聲之間具有統(tǒng)計(jì)獨(dú)立性，當(dāng)本體噪聲或目標(biāo)噪聲具有非常相似的倍頻特性且存在部分重疊時(shí)，濾波法的性能將急劇下降，為此，本文提出了一種基于U-Net[10]深度網(wǎng)絡(luò)的本體噪聲去除算法。首先，對(duì)水聽器采集到的噪聲信號(hào)或仿真信號(hào)作短時(shí)傅里葉變換（Short Time Fourier Transform，STFT）提取混合噪聲的時(shí)頻特征；其次，利用U-Net網(wǎng)絡(luò)對(duì)時(shí)頻譜的幅度譜（以下簡(jiǎn)稱頻譜）進(jìn)行分割，即從混合譜中提取出目標(biāo)噪聲的頻譜；最后，結(jié)合相位信息，經(jīng)逆變換（Inverse STFT，ISTFT）獲得分離后純凈的目標(biāo)噪聲。訓(xùn)練深度網(wǎng)絡(luò)需要大量標(biāo)定好的數(shù)據(jù)，但真實(shí)噪聲數(shù)據(jù)往往比較稀缺，因此，為了測(cè)試算法的可行性，本文首先采用仿真數(shù)據(jù)集開展水下噪聲分離研究。后續(xù)，將進(jìn)一步通過湖試和海試數(shù)據(jù)對(duì)網(wǎng)絡(luò)進(jìn)行優(yōu)化，提高系統(tǒng)的實(shí)用性和準(zhǔn)確性。

本文的主要貢獻(xiàn)如下：

1）面向移動(dòng)探測(cè)場(chǎng)景，在AUV本體噪聲與水中目標(biāo)噪聲頻譜特性和倍頻關(guān)系極其相似的情形下，提出了一種可行的本體噪聲抑制方案；

圖2 AUV本體噪聲去除算法流程圖Fig.2 Flow chart of AUV self-radiated noise removal algorithm

2）基于單通道水聽器信號(hào)，開展了水下目標(biāo)噪聲仿真和基于 U-Net深度學(xué)習(xí)網(wǎng)絡(luò)的水下目標(biāo)噪聲分離研究，實(shí)現(xiàn)了目標(biāo)噪聲與AUV本體噪聲的分離，大幅度提高了目標(biāo)噪聲的信噪比；

3）通過一系列仿真實(shí)驗(yàn)，測(cè)試了算法在不同轉(zhuǎn)速、槳葉數(shù)以及信噪比條件下的性能。

1 算法處理流程

1.1 噪聲信號(hào)仿真

為了測(cè)試算法的可行性，并獲得深度網(wǎng)絡(luò)的初始權(quán)值，我們利用現(xiàn)有的艦船噪聲模型生成仿真數(shù)據(jù)，合成訓(xùn)練網(wǎng)絡(luò)所需要的混合噪聲和目標(biāo)噪聲。文獻(xiàn)[11]指出，從聲源角度來說，艦船輻射噪聲主要由水動(dòng)力噪聲、機(jī)械噪聲和螺旋槳噪聲3大部分組成，這幾部分決定了組成艦船輻射噪聲頻譜的成分主要有線譜、連續(xù)譜和調(diào)制譜。這里選用與文獻(xiàn)[12]中相一致的方法生成輻射噪聲仿真數(shù)據(jù)，其時(shí)域信號(hào)表述為

式中，gX(t)，a(t)，gl(t)分別為連續(xù)譜分量、調(diào)制波形和線譜分量。圖3（a）和圖3（b）分別給出了一段仿真信號(hào)的時(shí)域波形和頻譜。

圖3 艦船輻射噪聲仿真Fig.3 An example for the ship-radiated noise simulation

基于上述模型，通過設(shè)定船只的槳葉數(shù)和轉(zhuǎn)速，我們可以獲取到不同轉(zhuǎn)速、不同槳葉數(shù)條件下艦船噪聲的時(shí)域仿真信號(hào)。為模擬混合噪聲，引入線性瞬時(shí)混合模型。記目標(biāo)船只輻射噪聲信號(hào)為st(t)；母船輻射噪聲信號(hào)為ss(t)；則混合噪聲的時(shí)域表達(dá)式記為sm(t)=ss(t)+ηst(t)，η表示目標(biāo)噪聲的相對(duì)權(quán)值。訓(xùn)練數(shù)據(jù)由混合噪聲和與之對(duì)應(yīng)的目標(biāo)噪聲信號(hào)組成，采樣頻率Fs=8 192 Hz，訓(xùn)練數(shù)據(jù)包括3 276段混合噪聲和目標(biāo)噪聲，其中驗(yàn)證集占20%，單段信號(hào)的時(shí)長(zhǎng)2 s左右。

1.2 時(shí)頻特征提取

在水下環(huán)境中，當(dāng)目標(biāo)距離水聽器較遠(yuǎn)時(shí)，其傳播過程中的能量衰減會(huì)很明顯，這就導(dǎo)致目標(biāo)噪聲的頻譜在混合噪聲的頻譜中不夠突出，淹沒在了背景噪聲之中。通過時(shí)頻變換，可以通過統(tǒng)計(jì)的方式對(duì)特征頻譜進(jìn)行增強(qiáng)。短時(shí)傅里葉譜是最原始的頻譜特征。出于簡(jiǎn)化特征提取和測(cè)試網(wǎng)絡(luò)可行性的考慮，本文直接通過STFT變換提取短時(shí)傅里葉譜。式（2）給出了STFT的定義：

式中：s(t)為源信號(hào)；w(t)為窗函數(shù)。變換窗長(zhǎng)取8 192，步長(zhǎng)為64，單段信號(hào)頻譜幀數(shù)取 128，最終頻譜大小取為4 096× 128。

1.3 U-Net網(wǎng)絡(luò)

本文將U-Net網(wǎng)絡(luò)用于本體噪聲去除任務(wù)，圖4給出了網(wǎng)絡(luò)的詳細(xì)結(jié)構(gòu)。降采樣部分采用卷積層和池化層來減小輸入圖像的尺寸并增加特征通道數(shù)，從而方便捕獲全局信息；升采樣部分采用轉(zhuǎn)置卷積來增加圖像尺寸并減少特征通道數(shù)，用于精確定位。

圖4 U-Net網(wǎng)絡(luò)結(jié)構(gòu)圖Fig.4 Network structure of U-Net

U-Net網(wǎng)絡(luò)配置如下：降采樣部分的二維卷積每層卷積核大小為3× 3、步長(zhǎng)為2，激活函數(shù)采用ReLU，引入批歸一化策略加速網(wǎng)絡(luò)的收斂；升采樣的每層使用與降采樣層對(duì)稱的卷積核大小和卷積步長(zhǎng)從而保證輸入輸出維度一致，其中卷積層采用的是反卷積，或稱轉(zhuǎn)置卷積。我們舍棄了文獻(xiàn)[10]中的 Dropout層，因?yàn)榕鷼w一化的引入可以在一定程度上抑制過擬合的出現(xiàn)，不需要額外引入 Dropout。網(wǎng)絡(luò)經(jīng)最后的 Sigmoid層輸出時(shí)頻掩模。

網(wǎng)絡(luò)訓(xùn)練的目標(biāo)是獲取準(zhǔn)確的時(shí)頻掩模。任意給定混合噪聲的頻譜，經(jīng)過“編碼–解碼”之后輸出掩模矩陣的維度與輸入一致，將時(shí)頻掩模與混合噪聲頻譜做點(diǎn)積便可獲得去除本體噪聲后的頻譜。記輸入混合信號(hào)的頻譜為M，期望輸出目標(biāo)噪聲的頻譜為T。訓(xùn)練過程的優(yōu)化器選用Adam，損失函數(shù)為輸出頻譜與期望頻譜差值的L1,1范數(shù)：

式中：f(M,Θ)為輸入頻譜；權(quán)重參數(shù)為Θ時(shí)對(duì)應(yīng)輸出的掩模；⊙表示按位相乘。

2 評(píng)價(jià)指標(biāo)

混合信號(hào)經(jīng)短時(shí)傅里葉變換后獲得的頻譜作為輸入特征，網(wǎng)絡(luò)對(duì)輸入頻譜進(jìn)行編解碼后獲得時(shí)頻掩模矩陣，其與輸入混合頻譜做點(diǎn)積之后即可獲得期望的目標(biāo)噪聲頻譜，再經(jīng)ISTFT即可獲取到時(shí)域目標(biāo)信號(hào)。

為了對(duì)本體噪聲去除效果進(jìn)行評(píng)估，我們引入盲分離中常用定量評(píng)價(jià)指標(biāo)SDR[13]，用以評(píng)估分離結(jié)果的優(yōu)劣[14]。記se為經(jīng)算法處理過后輸出的噪聲信號(hào)，其可以看作目標(biāo)信號(hào)與多種干擾信號(hào)之間的混合，即

式中：starget表示目標(biāo)噪聲；einterf,enoise,eartif分別為干擾、噪聲和偽影誤差。SDR定義為

分離前后目標(biāo)噪聲SDR的差值可以反映本體噪聲的去除效果，該指標(biāo)記為

式中：sr為目標(biāo)噪聲的參考值；sm為混合噪聲。這里SDR的具體計(jì)算基于mir_eval工具箱[16]。進(jìn)一步，為了衡量整體的分離效果，定義 GNSDR（Global NSDR）：

3 實(shí)驗(yàn)

為了測(cè)試上述算法在本體噪聲去除中的可行性和有效性，我們開展了以下3個(gè)實(shí)驗(yàn)：

1）母船與被探測(cè)船只槳葉數(shù)相同、轉(zhuǎn)速不同情形下的本體噪聲去除實(shí)驗(yàn)。該實(shí)驗(yàn)考察網(wǎng)絡(luò)分離2個(gè)不同轉(zhuǎn)速目標(biāo)的能力。在槳葉數(shù)相同的情況下，兩船轉(zhuǎn)速越接近，其頻譜的重合度越高，因此分離難度會(huì)隨轉(zhuǎn)速差的大小發(fā)生變化。

2）母船及被探測(cè)船只轉(zhuǎn)速均固定、槳葉數(shù)不相同情況下的本體噪聲去除實(shí)驗(yàn)。該實(shí)驗(yàn)考察網(wǎng)絡(luò)分離2個(gè)不同槳葉數(shù)目標(biāo)的能力。轉(zhuǎn)速固定時(shí)，隨槳葉數(shù)的變化目標(biāo)的頻譜特征也會(huì)發(fā)生變化，且槳葉數(shù)差異越大，其頻譜特征變換也越大。

3）固定頻率、固定槳葉數(shù)，目標(biāo)噪聲信噪比不同情形下的本體噪聲去除實(shí)驗(yàn)。該實(shí)驗(yàn)考察網(wǎng)絡(luò)分離2個(gè)噪聲強(qiáng)度不同的目標(biāo)的能力。當(dāng)目標(biāo)與母船相距較遠(yuǎn)時(shí)，目標(biāo)噪聲在傳播過程中會(huì)有一定程度的能量損失，合成混合信號(hào)時(shí)可通過調(diào)整η值改變混合噪聲中目標(biāo)噪聲的信噪比。

3.1 不同轉(zhuǎn)速差下分離效果評(píng)估

圖5給出了一個(gè)混合頻譜分割結(jié)果的例子。由于高頻部分包含的頻譜信息較少，這里只展示了512 Hz以下的頻譜。圖5（a）為混合噪聲的頻譜，混合噪聲由槳葉數(shù)均為Nblades=3的2段噪聲信號(hào)疊加而成，其中母船的轉(zhuǎn)速為RS=4r/s，目標(biāo)轉(zhuǎn)速為RT=13r/s，目標(biāo)噪聲信噪比0。圖5（b）為本文算法取得的結(jié)果，圖5（c）為參考目標(biāo)噪聲的頻譜。分析圖5不難看出來，在槳葉數(shù)相同、轉(zhuǎn)速差為dR=9r/s的情況下，本體噪聲的頻譜成分得到了有效去除。

圖5 算法效果示例Fig.5 An example for target noise separation

頻譜基頻對(duì)應(yīng)目標(biāo)螺旋槳的轉(zhuǎn)速。從圖5可以看出，分離混合噪聲頻譜時(shí)，由于目標(biāo)和本體2大主要聲源均為艦船輻射噪聲，均具有倍頻特性且存在部分重疊。因此，本體噪聲去除的質(zhì)量與兩船的轉(zhuǎn)速差直接相關(guān)。圖6給出了不同轉(zhuǎn)速差下本文提出的算法在本體噪聲去除上的性能，此時(shí)兩船槳葉數(shù)為NS=NT=3葉，給定母船轉(zhuǎn)速RS=4r/s，目標(biāo)轉(zhuǎn)速?gòu)腞T={12,13,14,15,16}r/s各值中選取從而設(shè)定不同轉(zhuǎn)速差，信噪比設(shè)為0。給定轉(zhuǎn)速差下各生成20組測(cè)試數(shù)據(jù)，驗(yàn)證算法性能。

圖6 不同轉(zhuǎn)速差情況下本體噪聲去除性能Fig.6 Performance of self-noise removal under different rotational speed difference

從圖6可看出，在轉(zhuǎn)速差從dR=8r/s增加到dR=10r/s過程中，GNSDR指標(biāo)明顯提升，即轉(zhuǎn)速差越大，分離效果越好。在轉(zhuǎn)速差繼續(xù)增大至dR=12r/s過程中，分離指標(biāo)值逐漸下降，這是因?yàn)榇藭r(shí)目標(biāo)噪聲基頻愈發(fā)靠近本體噪聲的倍頻，導(dǎo)致二者頻譜產(chǎn)生重疊，給分離帶來了一定挑戰(zhàn)。

3.2 不同槳葉數(shù)差下分離效果評(píng)估

不同槳葉數(shù)對(duì)應(yīng)不同的葉頻特征，兩船轉(zhuǎn)速固定時(shí)，槳葉數(shù)對(duì)頻譜的影響主要體現(xiàn)在葉頻上。圖7給出了不同槳葉數(shù)差下本文算法去除本體噪聲的性能，實(shí)驗(yàn)在兩船轉(zhuǎn)速固定條件下展開。本體和目標(biāo)轉(zhuǎn)速分別對(duì)應(yīng)為RS=8r/s和RT=11r/s，母船槳葉數(shù)NS=3葉，目標(biāo)槳葉數(shù)在NT={4,5,6,7}各值中選取從而設(shè)定不同槳葉數(shù)差，信噪比設(shè)為0。給定槳葉數(shù)差條件下，分別生成20個(gè)片段的測(cè)試數(shù)據(jù)。圖7的結(jié)果表明，隨著槳葉數(shù)差的增加，本體噪聲的分離效果不斷改善。這是因?yàn)殡S著槳葉數(shù)差的增加，兩噪聲成分的葉頻特征的區(qū)分愈發(fā)明顯。

圖7 不同槳葉數(shù)差條件下本體噪聲去除性能Fig.7 Performance of self-noise removal under different blade number difference

3.3 不同信噪比條件下分離效果評(píng)估

在給定探測(cè)場(chǎng)景下，當(dāng)目標(biāo)和本體船只槳葉數(shù)和轉(zhuǎn)速確定后，影響分離效果的另一主要因素當(dāng)屬目標(biāo)信號(hào)的信噪比。當(dāng)混合譜中二者強(qiáng)度相差不大時(shí)，聲源分離更容易開展，而實(shí)際場(chǎng)景中目標(biāo)信號(hào)強(qiáng)度一般較低。因此，對(duì)二者進(jìn)行混合時(shí)，可以設(shè)定不同的信噪比，驗(yàn)證低信噪比條件下算法去除本體噪聲的性能。圖 8給出了目標(biāo)信號(hào)不同信噪比條件下算法去除本體噪聲的性能。圖9給出了信噪比為–10 dB時(shí)分離效果的一個(gè)實(shí)例，圖9（a）表示的混合譜中目標(biāo)信號(hào)幾乎被完全淹沒，圖9（b）為本文算法輸出的目標(biāo)噪聲的頻譜，其與圖9（c）表示的目標(biāo)噪聲真實(shí)頻譜在特征譜線處基本一致。

圖8 不同信噪比條件下分割效果的對(duì)比圖Fig.8 Performance of self-noise removal under different signal-to-noise ratio

圖9 目標(biāo)信號(hào)信噪比為-10 dB條件下本體噪聲去除效果示例Fig.9 An example of target noise separation when the SNR of target noise is -10dB

實(shí)驗(yàn)結(jié)果表明，在目標(biāo)信號(hào)信噪比不低于–10 dB的條件下，應(yīng)用本文算法可以取得良好的本體噪聲去除效果。需要指出的是，網(wǎng)絡(luò)訓(xùn)練階段所用數(shù)據(jù)中目標(biāo)信號(hào)信噪比均為0（η=1），這也表明本文算法在低信噪比條件下具有一定的魯棒性。

4 結(jié)束語(yǔ)

本文面向水下移動(dòng)探測(cè)場(chǎng)景，開展了AUV本體噪聲去除算法研究。提取噪聲信號(hào)的時(shí)頻譜，輸入U(xiǎn)-Net深度學(xué)習(xí)網(wǎng)絡(luò)，實(shí)現(xiàn)了單水聽器噪聲信號(hào)下的盲源分離，初步達(dá)到了去除本體噪聲的目標(biāo)。不同實(shí)驗(yàn)條件下的仿真結(jié)果表明，在母船與目標(biāo)船只特征存在細(xì)微差異的情形下，本文的算法可以很好地對(duì)本體噪聲進(jìn)行去除，且該算法在低信噪比條件下表現(xiàn)良好。

盡管實(shí)驗(yàn)結(jié)果初步驗(yàn)證了算法的可行性，但是無(wú)論是特征提取、網(wǎng)絡(luò)結(jié)構(gòu)還是數(shù)據(jù)集都還需要開展進(jìn)一步的研究。后續(xù)，將主要在以下3個(gè)方面開展工作：

1）盡管 GNSDR指標(biāo)可以在一定程度上對(duì)盲源分離效果做定量說明，但是實(shí)際應(yīng)用中我們無(wú)法獲取目標(biāo)的真實(shí)噪聲數(shù)據(jù)，無(wú)法使用該指標(biāo)對(duì)結(jié)果進(jìn)行評(píng)估。因此，在后續(xù)工作中我們考慮將采集到的混合噪聲在去除背景噪聲前后分別輸入到識(shí)別模塊中，通過對(duì)比識(shí)別精度進(jìn)一步驗(yàn)證本體噪聲去除效果。

2）基頻信息是表征目標(biāo)特性的重要參數(shù)，其對(duì)聲源分離有著一定指導(dǎo)作用。后續(xù)工作中將引入多任務(wù)學(xué)習(xí)策略，實(shí)現(xiàn)聲源分離與基頻檢測(cè)的同步完成。

3）多目標(biāo)情形是AUV移動(dòng)探測(cè)中不可避免的場(chǎng)景，該條件下的盲源分離算法也是之后研究的重點(diǎn)。后續(xù)考慮使用聚類或及獨(dú)立成分分析的方法處理多目標(biāo)場(chǎng)景下的盲源分離任務(wù)。