王　瀟，金江明，盧奐采

（1.浙江工業(yè)大學(xué) 機(jī)械工程學(xué)院 特種裝備制造與先進(jìn)加工技術(shù)教育部/浙江省重點(diǎn)實(shí)驗(yàn)室，杭州 310014；　2.浙江省信號(hào)處理重點(diǎn)實(shí)驗(yàn)室，杭州 310014）

?

浮標(biāo)水聲信號(hào)中海浪聲的分離方法研究

王瀟1，金江明2，盧奐采1

（1.浙江工業(yè)大學(xué) 機(jī)械工程學(xué)院 特種裝備制造與先進(jìn)加工技術(shù)教育部/浙江省重點(diǎn)實(shí)驗(yàn)室，杭州 310014；2.浙江省信號(hào)處理重點(diǎn)實(shí)驗(yàn)室，杭州 310014）

針對(duì)浮標(biāo)系統(tǒng)測量的水聲信號(hào)開展研究，使用信號(hào)處理的方法對(duì)水聽器信號(hào)中的海浪聲成份進(jìn)行分離。首先對(duì)于16路水聽器的浮標(biāo)水聲信號(hào)，依據(jù)采樣頻率、采樣時(shí)間等參數(shù)將水聽器采集到的原始數(shù)字信號(hào)轉(zhuǎn)換成音頻信號(hào)，利用人耳區(qū)別水聽器信號(hào)中的各種噪聲成份；然后在時(shí)域內(nèi)，從原始水聽器信號(hào)中截取各種噪聲時(shí)段進(jìn)行頻譜分析，得到它們的特征頻率；最后設(shè)計(jì)合理的數(shù)字濾波器進(jìn)行去噪處理。試驗(yàn)結(jié)果表明水聽器信號(hào)中海浪聲成份得到有效抑制。

聲學(xué)；海浪聲；水聽器；特征頻率；數(shù)字濾波器

海浪聲是所有水聲信號(hào)在采集過程中不可避免的一種噪聲干擾，它會(huì)污染水聽器采集的原始信號(hào)，影響后期水聲信號(hào)特征的提取，降低水聲設(shè)備在定位或是環(huán)境監(jiān)測方面的精度，所以對(duì)水聽器原始信號(hào)中的海浪聲的分離研究具有重要的意義。為了有效地將海浪聲成份從原始水聽器信號(hào)中分離出去，就必須先分析原始水聽器信號(hào)中各種噪聲成份，并提取它們的特征。國內(nèi)的楊向峰、王盼盼等［1-2］使用功率譜估計(jì)的方法，對(duì)艦船［3］輻射的水聲信號(hào)特征提取進(jìn)行了研究；國外許多研究者使用短時(shí)傅里葉變換（STSF）對(duì)時(shí)變的水聲信號(hào)進(jìn)行分析［4-5］，實(shí)現(xiàn)分類識(shí)別的目的；李新欣［6］和裴善報(bào)等［7］使用小波變換對(duì)鯨類聲信號(hào)和水下連續(xù)爆炸［8］聲信號(hào)的特征進(jìn)行了提取；在2009年，李秀坤等［9］使用希爾伯特-黃變換（HHT），對(duì)水下目標(biāo)特征進(jìn)行了提取?？傊?，很多國內(nèi)外專家在對(duì)水下目標(biāo)特性的提取問題上，有較深入和完整的研究；但是，如何依據(jù)提取的特性來分離水聲信號(hào)中的各種噪聲成份，尤其是海浪聲，有待進(jìn)一步仔細(xì)的研究。文中將利用人耳來辨識(shí)水聽器信號(hào)中的噪聲成份，然后在時(shí)域截取各種噪聲時(shí)段，分析它們的頻率特性，最后針對(duì)海浪聲的特征頻率設(shè)計(jì)合理的數(shù)字濾波器，干凈地分離掉水聽器原始信號(hào)中的海浪聲成份。

1　浮標(biāo)水聲信號(hào)中噪聲成份分析

海洋環(huán)境很復(fù)雜，有來自海面上的噪聲，也有來自水下的噪聲，故使用浮標(biāo)水聲測量系統(tǒng)進(jìn)行海洋環(huán)境監(jiān)測時(shí)，會(huì)將干擾噪聲采集到原始信號(hào)中，這對(duì)浮標(biāo)水聲測量系統(tǒng)的監(jiān)測精度有很大的影響。針對(duì)水聽器采集到的原始數(shù)字信號(hào)進(jìn)行噪聲成份分析，通過采樣頻率、采樣時(shí)間等參數(shù)將水聽器采集到的原始數(shù)字信號(hào)轉(zhuǎn)換成音頻信號(hào)，利用人耳辨識(shí)水聽器信號(hào)中的噪聲成份，尋找出原始水聽器信號(hào)中的海浪聲成份。

通過分析發(fā)現(xiàn)一段原始的水聽器信號(hào)，由三種噪聲成份和有用的信號(hào)成份構(gòu)成，這三種噪聲成份分別是海浪聲、雷聲和機(jī)械撞擊聲，有用的信號(hào)成份指浮標(biāo)研究的目標(biāo)對(duì)象輻射的聲信號(hào)。假定原始水聽器信號(hào)為L，海浪聲成份為L1，雷聲成份為L2，機(jī)械撞擊聲成份為L3，以及有用的信號(hào)成份為Q。則原始水聽器信號(hào)L可表示為

現(xiàn)選取浮標(biāo)水聲測量系統(tǒng)的16路水聽器中第1路水聽器信號(hào)為分析對(duì)象。圖1是第1路水聽器的原始時(shí)域信號(hào)，信號(hào)的時(shí)長為60 s，采樣頻率為20 kHz，橫坐標(biāo)是時(shí)間（s），縱坐標(biāo)是聲壓（Pa）。

圖1　第1路水聽器的原始時(shí)域信號(hào)

在圖1中，海浪聲貫穿整個(gè)60 s的原始信號(hào)。圖中標(biāo)示出了三種典型噪聲成份發(fā)生時(shí)段，其中包含有6段海浪聲、6段雷聲、6段機(jī)械撞擊聲（圖1中橫向大括號(hào)標(biāo)出），這里的6段海浪聲是指海浪聲占總水聲信號(hào)中的主要成份（下文的海浪聲均表示這個(gè)意思，不另作解釋）。將這三種噪聲成份的具體數(shù)據(jù)統(tǒng)計(jì)在表1中。

通過圖1和表1發(fā)現(xiàn)，海浪聲成份的一些特征已經(jīng)被辨識(shí)出來。并且海浪聲和雷聲交替發(fā)生，海浪聲的幅值為3 Pa～4 Pa，雷聲的幅值為6 Pa～7 Pa，然而機(jī)械撞擊聲的幅值變化較大，從最小的7.6 Pa到最大的45 Pa不等。在下節(jié)中將對(duì)這三種噪聲成份進(jìn)行頻率特性分析。

表1　第1路水聽器的原始時(shí)域信號(hào)中各噪聲成份分析統(tǒng)計(jì)表

2　三種噪聲成份的頻率特性分析

在上一節(jié)中，分析了浮標(biāo)水聽器信號(hào)的噪聲成份。第1路水聽器的原始時(shí)域信號(hào)中主要包含三種噪聲成份，分別是海浪聲、雷聲和機(jī)械撞擊聲。為了從水聽器信號(hào)中抑制海浪聲成份引起的干擾，對(duì)海浪聲、雷聲和機(jī)械撞擊聲這三種噪聲成份的頻率特性進(jìn)行了分析。

首先，選取海浪聲進(jìn)行頻率特性分析，從60 s的原始信號(hào)中截取19.39 s～23.99 s處的水聲信號(hào)，幅值約為3 Pa。圖2為該段水聲信號(hào)的時(shí)域圖和頻譜分析圖。

在圖2（a）中，海浪聲時(shí)域信號(hào)的包絡(luò)面，近似一個(gè)簡諧波，其周期約為T=1.8 s，近似頻率為f=0.56 Hz。此時(shí)，觀察圖2（b）發(fā)現(xiàn)，海浪聲的特征頻率發(fā)生在0.6 Hz為中心的寬帶內(nèi)，與之前估算的頻率相似，故確定水聽器信號(hào)中含有的海浪聲的特征頻率為0.6 Hz為中心的寬帶，帶寬為1 Hz左右。

同理，分析水聽器信號(hào)中雷聲成份的特征頻率。從60 s的原始信號(hào)中截取17.4 s～19.46 s處的雷聲信號(hào)，其幅值約為6 Pa。圖3為該段雷聲信號(hào)的時(shí)域圖和頻譜分析圖。

從圖3（a）中可以發(fā)現(xiàn)，雷聲發(fā)生時(shí)段的時(shí)域數(shù)據(jù)無明顯的周期性，故不會(huì)出現(xiàn)固定的特征頻率，雷聲的特征頻率可能會(huì)在某一寬頻段出現(xiàn)。由于這一時(shí)段雷聲和海浪聲同時(shí)存在，根據(jù)上面對(duì)海浪聲特征頻率的分析，圖3（b）中0.6 Hz附近的寬帶屬于海浪聲的特征頻率，故水聽器信號(hào)中雷聲的頻率范圍約發(fā)生在5 Hz～600 Hz之間。

最后，對(duì)水聽器信號(hào)中機(jī)械撞擊聲成份的特征頻率進(jìn)行分析。從60 s的原始信號(hào)中截取41.2 s附近含有的機(jī)械撞擊聲信號(hào)，其幅值約為44 Pa。圖4為該時(shí)段含有機(jī)械撞擊聲的水聲信號(hào)的時(shí)域圖和頻譜分析圖。

通過觀察圖4（a）發(fā)現(xiàn)，機(jī)械撞擊聲類似一個(gè)脈沖信號(hào)，故它是寬頻帶信號(hào)。返觀圖4（b），并結(jié)合前面對(duì)海浪聲和雷聲的分析可知，頻率在0.6 Hz附近的寬帶的水聲信號(hào)成份，是海浪聲的貢獻(xiàn)。機(jī)械撞擊聲的特征頻率范圍為8 Hz～1 000 Hz，與雷聲的特征頻率范圍有所重疊。

3　海浪聲成份的分離與消除方法

3.1海浪聲分離方法

回顧前面對(duì)水聽器信號(hào)中包含的各類噪聲成份的分析，可知水聽器信號(hào)中的噪聲成份主要由海浪聲、雷聲和機(jī)械撞擊聲組成。其中，海浪聲的特征頻率比較固定，在0.6 Hz左右，雷聲和機(jī)械撞擊聲的頻率范圍，主要集中在5 Hz～1 000 Hz之間。

為了分離出海浪聲，消除其對(duì)水聽器信號(hào)的污染，依據(jù)海浪聲的頻率特征，以及它與雷聲和機(jī)械撞擊聲的頻率范圍基本沒有重疊的情況，可對(duì)含有海浪聲的水聽器信號(hào)加入Butterworth數(shù)字帶通濾波器，帶通濾波器的中心頻率設(shè)計(jì)在0.6 Hz，帶通寬度為1.18 Hz，這樣可提取出純頻的海浪聲信號(hào)，如圖5（a）所示。然后，將原始的水聽器記錄信號(hào)減掉純頻的海浪聲信號(hào)，則可以得到分離出掉海浪聲后的水聽器信號(hào)，如圖5（b）所示。

3.2海浪聲分離處理

首先，選取水聲信號(hào)中發(fā)生在19.6 s～23.53 s處的海浪聲信號(hào)。對(duì)這一時(shí)段的信號(hào)，進(jìn)行中心頻率為0.6 Hz、帶通寬度為1.18 Hz的帶通濾波，得到如圖5（a）所示的海浪聲信號(hào)，其中細(xì)線代表水聽器原始記錄信號(hào)，粗線代表經(jīng)過分離處理的信號(hào)。最后，從原始的水聽器記錄信號(hào)中，分離掉帶通濾波器獲得的海浪聲信號(hào)，得到如圖5（b）所示的消除了海浪聲后的信號(hào)。在此信號(hào)中，海浪聲基本被完全分離掉，水聲信號(hào)不再隨海浪聲波動(dòng)。

圖2　發(fā)生在19.39 s～23.99 s處的海浪聲信號(hào)時(shí)域圖和頻譜分析圖

圖3　發(fā)生在17.4 s～19.46 s處的雷聲信號(hào)時(shí)域圖和頻譜分析圖

圖4　發(fā)生在41.2 s附近的機(jī)械撞擊聲信號(hào)時(shí)域圖和頻譜分析圖

圖5　海浪聲分離效果圖

根據(jù)以上的海浪聲分離方法，對(duì)第1路水聽器信號(hào)進(jìn)行海浪聲信號(hào)的分離處理，將會(huì)得到兩幅包含有兩條對(duì)比曲線的時(shí)域圖，如圖6和圖7所示。其中，圖6包含水聽器原始記錄信號(hào)和分離出的海浪聲信號(hào)，圖7包含水聽器原始記錄信號(hào)和分離出海浪聲后的水聽器信號(hào)。這兩幅時(shí)域圖中，圖6反映出分離出的海浪聲，很好地吻合了水聽器原始記錄信號(hào)上下包絡(luò)線的中心軌跡。圖7反映出海浪聲基本從水聽器原始記錄信號(hào)中被分離出去。

圖6　第1路水聽器的原始記錄信號(hào)與分離出的海浪聲信號(hào)

圖7　第1路水聽器的原始記錄信號(hào)與分離掉海浪聲后的信號(hào)

3.3海浪聲分離結(jié)果總結(jié)

通過對(duì)第1路水聽器信號(hào)進(jìn)行海浪聲的分離處理和分析，發(fā)現(xiàn)該分離方法可以較精確地將海浪聲從水聽器原始記錄信號(hào)中分離出來，并保留水聽器原始記錄信號(hào)中的其他聲音成份。這樣，達(dá)到了研究所要求的從水聽器原始信號(hào)中，分離海浪聲成份的目的?，F(xiàn)給出以聲壓級(jí)形式表示的海浪聲成份分離前和分離后的各噪聲成份聲壓幅值的計(jì)算公式，見式（2）。其中，P是以聲壓級(jí)形式表示的海浪聲成份分離前和分離后的各噪聲成份的聲壓幅值，單位是dB；pe是海浪聲成份分離前和分離后的各噪聲成份時(shí)域信號(hào)的聲壓幅值，單位是Pa；pref是參考聲壓，它在水中的取值為1 μPa。

經(jīng)過計(jì)算，將以聲壓級(jí)形式表示的海浪聲成份分離前和分離后的聲壓幅值統(tǒng)計(jì)到表2中。

表2　第1路水聽器信號(hào)進(jìn)行數(shù)字濾波前后各類特征聲音幅值對(duì)比統(tǒng)計(jì)/dB

4　結(jié)語

針對(duì)擁有16路水聽器的浮標(biāo)水聲測量系統(tǒng)，選取第1路水聽器信號(hào)為研究對(duì)象，分析水聽器信號(hào)中包含的噪聲成份，并依據(jù)各個(gè)噪聲成份的頻率特性，提出分離水聽器信號(hào)中海浪聲成份的方法。研究結(jié)論如下：

（1）利用人耳辨識(shí)出水聽器原始信號(hào)中包含有三種噪聲成份，即海浪聲、雷聲和機(jī)械撞擊聲。

（2）通過研究三種噪聲成份的頻率特性，發(fā)現(xiàn)海浪聲的特征頻率為0.6 Hz左右；而雷聲和機(jī)械撞擊聲的特征頻率為一寬頻帶，并有所重疊，主要集中在5 Hz～1 000 Hz之間。

（3）水聽器原始記錄信號(hào)經(jīng)過設(shè)計(jì)的數(shù)字濾波器，可以很好地分離出不同時(shí)段的海浪聲，從而可實(shí)現(xiàn)抑制和消除海浪聲的目的。經(jīng)統(tǒng)計(jì)，在水聽器原始記錄信號(hào)中，海浪聲可以從原來的130 dB下降到了0 dB，達(dá)到了對(duì)水聽器原始記錄信號(hào)中海浪聲成份抑制的目的。

［1］楊向鋒，張效民，孫繼紅.艦船輻射噪聲功率譜特征提取方法研究［J］.魚雷技術(shù)，2006，14（1）：35-38.

［2］王盼盼，張國軍，關(guān)凌綱，等.艦船輻射噪聲特征線譜提取方法研究［J］.海洋技術(shù)，2010，29（3）：47-50.

［3］申志偉，李榮，李志華.特征相似的艦船噪聲模型可信性分析［J］.噪聲與振動(dòng)控制，2012，32（2）：159-162.

［4］ALI PEZESHKI，MAHMOOD R.AZIMI-SADJADI and LOUIS L.SCHARF.Undersea target classification using canonical correlation analysis［J］.IEEE Journal of Oceanic Engineering，2007，32（4）:948-955.

［5］CORNEL IOANA，ANDRé QUINQUIS and YANN STEPHAN.Feature extraction from underwater signals using time-frequency warping operators［J］.IEEE Journal of Oceanic Engineering，2006，31（3）:628-645.

［6］李新欣.船舶及鯨類聲信號(hào)特征提取和分類識(shí)別研究［D］.哈爾濱：哈爾濱工程大學(xué)，2012.

［7］裴善報(bào)，劉榮忠，郭銳.基于小波變換的水下連續(xù)爆炸聲信號(hào)特征分析［J］.爆炸與沖擊，2015，35（4）：520-526.

［8］翟紅波，李芝絨，張玉磊，等.典型艦船艙室水下爆炸響應(yīng)特性［J］.噪聲與振動(dòng)控制，2016，36（1）：60-64.

［9］李秀坤，謝磊，秦宇.應(yīng)用希爾伯特黃變換的水下目標(biāo)特征提?。跩］.哈爾濱工程大學(xué)學(xué)報(bào)，2009，30（5）：542-546.

Separation of Ocean Wave Sound from the Underwater Sound Signals Collected by the Buoy

WANGXiao1，JIN Jiang-ming2，LU Huan-cai1

（1.Key Laboratory of E and M，Ministry of Education and Zhejiang Province，Zhejiang University of Technology，Hangzhou 310014，China；2.Key Laboratory of Signal Processing of Zhejiang Province，Hangzhou 310014，China）

The underwater sound signals collected by the buoys are studied.A signal processing method is applied to separate the ocean wave sound from the original measured signals with hydrophones.Each measured hydrophone signal of the 16 channels is transformed into audio signal based on the parameters such as sampling frequency，sampling time and so on.Then，the audio signal is played and classified into different segments with different noise characteristics by human ears. These segments with different noise characteristics are extracted from the original time-domain signals for frequency spectral analysis，and their characteristic frequencies are obtained.Finally，the digital filter is designed to de-noise the ocean wave sound.Experimental results show that the ocean wave sound components are effectively suppressed from the measured original signals by hydrophones.

acoustics；ocean wave sound；hydrophone；characteristic frequency；digital filter

TB565

ADOI編碼：10.3969/j.issn.1006-1335.2016.04.035

1006-1355（2016）04-0165-04+192

2016-03-04

王瀟（1987-），男，江蘇省泰興市人，碩士生，主要研究方向?yàn)樗曅盘?hào)處理。E-mail:wangxiao_svlab@163.com

盧奐采（1962-），女，博士生導(dǎo)師。E-mail:huancailu@zjut.edu.cn