深海海面目標單水聽器被動測距方法與驗證

翁晉寶,楊燕明

(1.國家海洋局第三海洋研究所 海洋聲學與遙感實驗室，福建 廈門 361005)

?

深海海面目標單水聽器被動測距方法與驗證

翁晉寶1,楊燕明1*

(1.國家海洋局第三海洋研究所 海洋聲學與遙感實驗室，福建 廈門 361005)

基于射線理論分析了在深海情況下海面聲源產(chǎn)生聲場的頻率-距離干涉結(jié)構(gòu)，給出了影區(qū)內(nèi)聲場頻率-距離干涉結(jié)構(gòu)的近似理論表達式，分析得到影區(qū)內(nèi)聲場頻域干涉周期隨收發(fā)距離的增加而增大、隨著接收水聽器深度的增加而減小。因此由單水聽器記錄的聲場干涉結(jié)構(gòu)即可實現(xiàn)被動聲源距離估計。在南海深海實驗中觀測到海面寬帶噪聲源在聲場影區(qū)形成的聲場干涉結(jié)構(gòu)，對實驗獲得聲場干涉結(jié)構(gòu)的處理結(jié)果驗證了深海聲場影區(qū)干涉結(jié)構(gòu)用于被動聲源距離估計的有效性。與傳統(tǒng)的匹配場被動定位方法相比，該方法不需要已知海底聲學參數(shù)和大規(guī)模的拷貝場計算。

海面目標；深海聲場影區(qū)；聲場干涉結(jié)構(gòu)；被動測距

1 引言

傳統(tǒng)的水下被動目標定位系統(tǒng)一般需要多個基元對目標實現(xiàn)定位，例如三子陣法[1]和球面內(nèi)插法[2]，這類方法利用聲信號到達時間或相位進行測距，對基陣布放精度要求高，工作頻率高，測距距離近，精度低。近些年來，定位技術(shù)的發(fā)展方向主要有兩大類，匹配場定位技術(shù)[3]和目標運動分析[4]。匹配場定位技術(shù)利用已知的海洋環(huán)境參數(shù)，采用適當?shù)穆晥瞿Ｐ陀嬎隳M的接收信號，然后與實際測量得到的聲信號進行相關(guān)“匹配”，實現(xiàn)對目標聲源定位。這種方法需要解決的主要是環(huán)境、聲場模型和基陣系統(tǒng)失配等問題，在復雜海洋環(huán)境中往往受限。目標運動分析通過一段時間數(shù)據(jù)的信息綜合，利用如方位、頻率和速度等信息來確定目標位置，然而該方法要求觀測平臺必須機動，而且收斂時間長，測距精度也較低。近年來隨著水聲技術(shù)的發(fā)展，利用波導不變量進行目標定位受到了國內(nèi)外研究人員[5—8]的關(guān)注，該類算法的優(yōu)點是對環(huán)境具有更好的寬容性。

在淺海波導中聲波以簡正波疊加的形式傳播，不同號簡正波發(fā)生干涉，使聲場具有穩(wěn)定的頻率-空間干涉結(jié)構(gòu)，表現(xiàn)為聲強的頻率-距離圖上出現(xiàn)明暗相間的干涉條紋。Chuprov[9]用一個簡單的公式定義了波導不變量：

(1)

用于描述干涉條紋斜率與距離和頻率的關(guān)系。波導不變量從一個全新的角度揭示了聲波在海洋中傳播的結(jié)構(gòu)特征，在水聲定位中得到應用。文獻[5—8]描述了利用單陣元進行淺海目標的波導不變量測距方法。Zhao等[10]基于波導不變量，推導出了聲源距離與聲場干涉譜頻域周期的關(guān)系，利用引導源可實現(xiàn)單個水聽器寬帶運動聲源距離估計。戚聿波等[11]基于淺海簡正波水平波數(shù)差與波導不變量之間的關(guān)系，提出了一種適用于水平不變淺海聲波導中接收信號自相關(guān)函數(shù)的頻域卷繞變換算子，可用于被動聲源距離估計。任云等[12]基于聲強的距離-頻率干涉結(jié)構(gòu)，利用對水平陣不同陣元信號的聲強譜進行頻域補償?shù)姆椒ü烙嫴▽Р蛔兞亢途嚯x，在波導不變量和距離其中一項已知的情況下，可以通過該方法估計另外一項。Thode等[13]將波導不變量原理與匹配場處理結(jié)合，一定程度上提高了匹配場定位的性能。

對于深海低頻聲場，根據(jù)射線理論由于經(jīng)多次海底反射的聲線幅度較小可忽略，實際上對于聲場影響較大的聲線數(shù)目不多，便于分析深海聲場干涉結(jié)構(gòu)。2004年吳國清[14]分析了淺海近距離(數(shù)倍海深)范圍內(nèi)，直達聲、海面反射聲和海底反射聲互相干涉引起的譜干涉現(xiàn)象，利用聲線直線近似和聲速均勻近似幾何推導了干涉譜與聲源深度、距離的關(guān)系，從而根據(jù)實測干涉譜反推聲源深度和距離，并利用實驗數(shù)據(jù)實現(xiàn)淺海水平距離500 m范圍內(nèi)水下目標被動測距。翁晉寶等[15]基于射線理論分析當聲源與接收水聽器位于深度較淺的表面附近時聲源激發(fā)的典型深海聲場頻率-距離干涉結(jié)構(gòu)，在第一影區(qū)內(nèi)由固定的4條海底一次反射聲線干涉疊加，在這種情況下聲場干涉結(jié)構(gòu)在固定距離下聲場的頻率干涉周期固定，因此聲場的干涉條紋在頻率域上等間隔分布，在深海影區(qū)目標定位中具有應用前景。此外，翁晉寶等[16]基于深海近距離聲場頻率-距離干涉結(jié)構(gòu)，提出了一種當聲源與接收水聽器分別位于海面與海底附近時匹配聲場干涉周期反演深海沉積層聲速與厚度的方法。

針對深海海洋環(huán)境，本文首先基于射線理論分析了海面聲源產(chǎn)生聲場的頻率-距離干涉結(jié)構(gòu)，得到影區(qū)內(nèi)干涉結(jié)構(gòu)的近似理論表達式，給出聲強頻域干涉周期與收發(fā)距離的變化關(guān)系，研究深海海洋環(huán)境中海面聲源距離被動估計問題，并且利用海上實驗數(shù)據(jù)進行了驗證。

2 深海海面聲源聲場干涉結(jié)構(gòu)分析

在典型深海海洋環(huán)境下，當聲源與接收水聽器位于海水表層時，根據(jù)接收水聽器與聲源的水平距離可以將典型深海聲場在空間上劃分為直達聲作用區(qū)、影區(qū)及會聚區(qū)，直達聲作用區(qū)指聲線未經(jīng)海底反射或未經(jīng)海底附近反轉(zhuǎn)的作用區(qū)域，影區(qū)指直達聲線與反轉(zhuǎn)聲線無法到達的區(qū)域，會聚區(qū)指反轉(zhuǎn)聲線形成的高強度焦散的區(qū)域。典型5 000 m深海聲場聲線如圖1所示，其中聲源深度為100 m、聲速剖面如圖2所示。

圖1 聲場區(qū)域的劃分Fig.1 The division of sound field region

圖2 聲速剖面Fig.2 The sound velocity profile

圖3 對影區(qū)聲場起主要貢獻的聲線Fig.3 The main contributing sound rays to sound field of the shadow zone

根據(jù)文獻[15]，在影區(qū)內(nèi)，假設(shè)多次海底反射聲線由于能量衰減大而對聲場的貢獻可以忽略不計，第一影區(qū)的聲場僅考慮一次海底反射聲線的貢獻，在所示的典型深海條件下，當聲源與接收水聽器位于海水表層時，如圖3所示，一次海底反射聲線包含4條聲線路徑：聲源—海底—接收器、聲源—海面—海底—接收器、聲源—海底—海面—接收器和聲源—海面—海底—海面—接收器，則接收水聽器接收到的聲壓可以近似為：

(2)

式中，S(f)表示聲源的譜級，W表示聲源單位立體角內(nèi)的輻射聲功率，F(xiàn)、R、Vb和t分別表示聲線的聚焦因子、斜距、海底聲壓反射系數(shù)和傳播時間，變量的下標依次代表海底反射次數(shù)和聲線的序號。

圖4 聲源掠射角Fig.4 The source angle

當聲源位于海水表層時，聲線1與聲線2的軌跡基本一致，如圖4所示，聲源掠射角基本一致(不考慮正負)即α11≈α12，聚焦因子F、斜距R和海底聲壓反射系數(shù)Vb也基本一致，令Δt1=t12-t11，則式(2)中聲線1與聲線2的相干疊加部分可以近似表示為：

(3)

α11與α12差別較小，Δt1可以近似表示為[15]：

(4)

由于聲源位于海水表層則聲線3與聲線4的軌跡也基本一致，聲線的聲源掠射角基本一致即α13≈α14，聚焦因子F、斜距R和海底聲壓反射系數(shù)Vb也基本一致，則式(2)中聲線3與聲線4的相干疊加部分可以近似表示為：

(5)

α13與α14差別較小，t14-t13可以近似表示為[15]：

(6)

(7)

令Δt2=t13-t11，則式(2)中4條聲線相干疊加部分可以近似為：

(8)

由于接收水聽器位于海水表層時，聲線1與聲線3的軌跡相近，α11與α13差別較小，則Δt2可以近似表示為[15]：

(9)

綜上所述接收水聽器接收到的聲壓可以近似為：

×(1-ei2πfΔt1)(1-ei2πfΔt2)，

(10)

接收到的聲強可以近似為：

(11)

其中，

Δt1=t12-t11，

(12)

Δt2=t13-t11，

(13)

即聲強隨頻率具有兩個干涉周期，當頻率f滿足

f=m/Δt1(m=1,2,3,…)，

(14)

或者

f=n/Δt2(n=1,2,3,…)，

(15)

時聲強為極小值，在聲強頻率-距離圖上表現(xiàn)為干涉相消的條紋。

當聲源為船舶等海面目標時，由于船舶吃水深度一般僅為幾米，即聲源深度僅有幾米，則聲源—海底—接收器和聲源—海面—海底—接收器的聲線傳播時間差Δt1極小，即頻域干涉周期1/Δt1極大，約數(shù)百赫茲。但是，考慮到船舶噪聲中高頻(大于1 kHz)部分的譜級小于低頻(小于1 kHz)部分的譜級，和海洋波導中高頻聲傳播的衰減比低頻聲傳播的衰減大的特點，以及海洋中本身存在的環(huán)境噪聲，在實際深海影區(qū)接收信號中僅在數(shù)百赫茲的低頻部分可以觀測到船舶噪聲形成的干涉條紋。因此，對于船舶等海面噪聲源，利用單水聽器難以在影區(qū)內(nèi)觀測到式(14)代表的第一種干涉結(jié)構(gòu)。

在實際應用中，接收水聽器的深度往往為幾十米或者數(shù)百米，則聲源—海底—接收器和聲源—海底—海面—接收器的聲線傳播時間差Δt2較大，即頻域干涉周期1/Δt2較小，約幾十赫茲或幾赫茲，便于利用單水聽器在影區(qū)內(nèi)低頻段觀測到式(15)代表的第二種干涉結(jié)構(gòu)。根據(jù)文獻[15]，在固定距離的情況下，Δt2隨接收深度z增加而增大，頻域干涉周期1/Δt2隨接收深度z增加而減小；在固定接收深度的情況下，Δt2隨距離增加而減小，頻域干涉周期1/Δt2隨距離增加而增大。

為了驗證上述理論分析，首先給出由BELLHOP[17]程序計算得到的海面目標產(chǎn)生的典型深海聲場干涉結(jié)構(gòu)。海水深度為5 000 m，聲源深度為5 m，接收深度為250 m，海水聲速如圖2所示，海底聲速為1 600 m/s，海底密度為1.8 g/cm3，海底吸聲系數(shù)為0.8 dB/λ，聲源帶寬為50～250 Hz。圖5給出了收發(fā)距離在4～55 km內(nèi)的聲場頻率-距離干涉結(jié)構(gòu)，圖6給出了收發(fā)距離在4～55 km內(nèi)的1/Δt2，1/Δt2在第一影區(qū)內(nèi)隨著距離的增加而單調(diào)增大，其中聲線傳播時間t11和t13由BELLHOP程序計算得到，圖7給出了根據(jù)式(15)預測的聲場干涉相消條紋，圖8給出了將圖7預測的干涉相消條紋與圖5仿真得到的聲場頻率-距離干涉結(jié)構(gòu)進行比對，從圖8可以看出兩者完全一致，即式(15)可以準確的預測影區(qū)內(nèi)干涉相消條紋。圖9給出了在同樣的條件下當接收深度為125 m時的聲場干涉結(jié)構(gòu)，將圖9和圖5進行對比，可以看出聲場頻域干涉周期隨接收深度的增加而減小。

圖5 深海聲場頻率-距離干涉結(jié)構(gòu)(接收深度250 m)Fig.5 The sound field frequency-range interference patterns in deep water when the receiver located at a depth of 250 m

圖6 聲場頻域干涉周期1/Δt2Fig.6 The sound field interference cycle in the frequency domain

圖7 根據(jù)式(15)預測的聲場干涉相消條紋Fig.7 The sound field stripes of interference cancellation cycle predicted by equation (15)

圖9 深海聲場頻率-距離干涉結(jié)構(gòu)(接收深度125 m)Fig.9 The sound field frequency-range interference pat-terns in deep water when the receiver located at a depth of 125 m

3 單水聽器被動聲源距離估計方法

通過上面的理論分析與仿真結(jié)果可知，影區(qū)內(nèi)海面目標形成聲場干涉結(jié)構(gòu)的頻域干涉周期隨距離的增加而單調(diào)增大。因此在已知接收水聽器深度、海底地形和海水聲速剖面的前提下，通過仿真計算得到如圖6所示的頻域干涉周期隨水平距離的變化情況，之后只要能獲得影區(qū)內(nèi)海面目標形成聲場干涉結(jié)構(gòu)的實測頻域干涉周期，通過兩者之間的匹配便可以估計海面目標與接收水聽器的水平距離。在實際情況中，影區(qū)內(nèi)的傳播損失隨著影區(qū)號數(shù)的增加而不斷增大，一般海面目標形成的聲場干涉結(jié)構(gòu)由于目標聲源級的限制往往僅在第一影區(qū)內(nèi)能被觀測到。具體實施步驟如下：

(1)海面目標形成聲場干涉結(jié)構(gòu)的判斷

將接收水聽器布放于海水表層(例如幾十米或數(shù)百米)持續(xù)接收記錄寬帶噪聲信號，對寬帶噪聲信號進行連續(xù)時頻分析，觀察時頻譜圖中數(shù)百赫茲內(nèi)的低頻部分，如果存在明顯的單一的明暗相間的條紋結(jié)構(gòu)，且相鄰暗條紋的頻域間隔不隨頻率變化，則判斷在第一影區(qū)內(nèi)存在海面目標。

如果在同一頻帶內(nèi)存在兩種明暗相間的條紋結(jié)構(gòu)疊加在一起，且兩種條紋結(jié)構(gòu)的頻域間隔都不隨頻率變化，則可能是第一影區(qū)內(nèi)存在水下目標，或者第一影區(qū)內(nèi)存在兩個海面目標，難以進行判斷，暫不討論。

(2)頻域干涉周期的提取

當判斷在第一影區(qū)內(nèi)存在海面目標后，時頻譜圖上存在著明顯的單一的明暗相間的條紋結(jié)構(gòu)，且相鄰暗條紋的頻域間隔不隨頻率變化，因此可以在每個時刻提取包含N個亮條紋的頻域間隔ΔfN(t)，之后通過Δf=ΔfN/N得到每個時刻的頻域干涉周期Δf(t)，即通過平均計算減少誤差。

(3)根據(jù)海洋環(huán)境參數(shù)仿真計算頻域干涉周期

根據(jù)布放海域采集的海水聲速剖面、海深和接收深度，假定聲源深度為5 m，利用BELLHOP程序計算得到第一影區(qū)內(nèi)不同距離處聲源-海底-接收器聲線到達時間t11(x)和聲源—海底—海面—接收器聲線到達時間t13(x)，之后根據(jù)Δf2=1/(t13-t11)計算第二種頻域干涉周期Δf2(x)隨水平距離的變化情況。

(4)匹配頻域干涉周期估計聲源距離

將提取得到的每個時刻的頻域干涉周期Δf(t)，與根據(jù)海洋環(huán)境參數(shù)仿真計算得到的每個距離的第二種頻域干涉周期Δf2(x)進行匹配，得到海面目標與接收水聽器的水平距離x(t)隨時間的變化情況。

4 實驗數(shù)據(jù)分析

實驗數(shù)據(jù)來自于在南海海域進行的聲學實驗，由錨定于海底的一條等間距32陣元垂直線列陣潛標記錄噪聲信號，垂直陣列陣元間隔30 m、覆蓋深度為188～1 118 m，由固定于附近實驗船上的船舶自動識別系統(tǒng)(Automatic Identification System, 簡稱AIS)記錄經(jīng)過接收點附近的海面船只信息。圖10給出了實驗海域聲速剖面，水深3 495 m。

圖10 實驗海域聲速剖面Fig.10 The sound velocity profile

圖11、圖12和圖13分別給出了陣元1、5和10接收的當天20:00-20:10記錄的噪聲信號的時頻分析結(jié)果，3個陣元所處深度依次為188 m、308 m和458 m，從這3張圖中可以看出該時間段內(nèi)聲場具有兩種干涉條紋，在70～160 Hz頻段內(nèi)干涉條紋的頻域干涉周期隨著時間增大逐漸變小，說明70～160 Hz頻段內(nèi)干涉條紋對應的是接近接收點的目標，在160～230 Hz頻段內(nèi)干涉條紋的頻域干涉周期隨著時間增大逐漸變大，說明160～230 Hz頻段內(nèi)干涉條紋對應的是遠離接收點的目標。其中，在70～160 Hz頻段內(nèi)，第1陣元和第5陣元的干涉條紋較為清晰；在160～230 Hz頻段內(nèi)，第5陣元和第10陣元的干涉條紋較為清晰。另外，圖11～13中干涉條紋的頻域干涉周期明顯隨深度的增加而減小，與以上理論分析是一致的。

圖11 第1陣元(188 m)接收噪聲信號的時頻分析Fig.11 The time-frequency analysis of the received noise signal in channel 1(188 m)

圖12 第5陣元(308 m)接收噪聲信號的時頻分析Fig.12 The time-frequency analysis of the received noise signal in channel 5 (308 m)

圖13 第10陣元(458 m)接收噪聲信號的時頻分析Fig.13 The time-frequency analysis of the received noise signal in channel 10 (458 m)

圖12第5陣元接收噪聲信號的時頻分析中，在70～160 Hz頻段和160～230 Hz頻段都存在明顯的單一的明暗相間的條紋結(jié)構(gòu)，且相鄰暗條紋的頻域間隔不隨頻率變化，考慮到接收水聽器深度為308 m，可以判斷該條紋結(jié)構(gòu)由聲源—海底—接收器、聲源—海底—海面—接收器干涉形成。圖14給出了由圖12中70～160 Hz頻段和160～230 Hz頻段分別直接提取的4條暗條紋，由于該條紋結(jié)構(gòu)為明暗相間的干涉條紋，單條亮條紋的寬度即為兩條相鄰暗條紋的頻域間隔，其中70～160 Hz頻段的兩條暗條紋中間間隔9條亮條紋、160～230 Hz頻段的兩條暗條紋中間間隔5條亮條紋，因此可以通過平均計算得到相鄰暗條紋的頻域間隔，即頻域干涉周期，圖15給出了根據(jù)圖14計算的兩個頻段內(nèi)干涉條紋的頻域干涉周期，其中70～160 Hz頻段干涉周期在9～592 s內(nèi)從6.0 Hz減小到4.7 Hz、160～230 Hz頻段干涉周期在9～592 s內(nèi)從8.6 Hz增大到9.4 Hz。

圖14 第5陣元接收噪聲信號的干涉條紋Fig.14 The interference stripes of the received noise signal in channel 5

圖15 第5陣元接收噪聲信號的頻域干涉周期Fig.15 The interference cycle in the frequency domain of the received noise signal in channel 5

圖16給出了根據(jù)實測海洋環(huán)境計算的頻域干涉周期隨收發(fā)距離的變化情況，頻域干涉周期在5～25 km內(nèi)從3.0 Hz增大到11.0 Hz。表1給出了根據(jù)圖15實測頻域干涉周期與圖16仿真頻域干涉周期的對比估計出的兩個海面目標的水平距離，70～160 Hz頻段內(nèi)目標水平距離在9～592 s內(nèi)從14.27 km減小到11.00 km，160～230 Hz頻段內(nèi)目標水平距離在9～592 s內(nèi)從20.35 km增大到22.06 km。

圖16 根據(jù)實測海洋環(huán)境計算的頻域干涉周期Fig.16 The calculated interference cycle in the frequency domain based on practical marine environment

Tab.1 The estimated range of targets based on the frequency-domain interference cycle

表2給出了9～592 s內(nèi)AIS系統(tǒng)提供的附近海面所有船舶的水平距離信息，該時間段內(nèi)實驗海域內(nèi)存在4個目標，目標1的水平距離從16.61 km減小到12.10 km，目標1的船速為12.3 kn、航向為209°(正北為0°，順時針旋轉(zhuǎn))，目標1逐步靠近接收點，目標2的水平距離從21.09 km增大到21.73 km，目標2的船速為14.1 kn、航向為214°，目標2逐步遠離接收點，目標3和4的水平距離都保持在30 km以上．表3給出了該時間段內(nèi)根據(jù)AIS系統(tǒng)提供的水上移動通信業(yè)務(wù)標識碼(Maritime Mobile Service Identify，簡稱MMSI)查詢得到的附近海面船舶的船長、船寬和總噸，目標1和2的規(guī)模遠遠大于目標3和4。

根據(jù)9～592 s內(nèi)估計海面目標水平距離(表1)和實測海面目標水平距離(表2)的對比，以及目標1、2的船速、航向信息，可以判斷產(chǎn)生70～160 Hz頻段內(nèi)干涉條紋的是目標1、產(chǎn)生160～230 Hz頻段內(nèi)干涉條紋的是目標2，如圖17所示。圖18給出了目標1和2的估計距離與實測距離的對比，目標1和2的距離估計平均誤差為11.3%和1.4%，對于目標1的距離估計誤差大于對于目標2的距離估計誤差。圖12中70～160 Hz頻段內(nèi)干涉條紋與160～230 Hz頻段內(nèi)干涉條紋相比，干涉條紋間隔較小較模糊不利于干涉條紋的提取。

表2 實驗海域內(nèi)海面船舶與接收點的水平距離

表3 實驗海域內(nèi)海面船舶的船長、船寬和總噸

5 結(jié)論

當海面目標位于深海聲場第一影區(qū)內(nèi)，由4條固定的海底一次反射聲線干涉疊加，本文給出了干涉結(jié)構(gòu)的近似理論表達式，分析得到在影區(qū)內(nèi)聲場頻域干涉周期隨收發(fā)距離的增加而增大、隨著接收水聽器深度的增加而減小?；陬l域干涉周期與收發(fā)距離的關(guān)系，提出了一種單水聽器的被動聲源測距方法，利用南海深海海域獲得的噪聲測量實驗數(shù)據(jù)和AIS數(shù)據(jù)對方法的有效性進行了驗證。相對于吳國清[14]分析淺海近距離(數(shù)倍海深)范圍內(nèi)直達聲、海面反射聲和海底反射聲互相干涉引起的譜干涉現(xiàn)象，利用聲線直線近似和聲速均勻近似幾何推導了干涉譜與聲源深度、距離的關(guān)系，從而實現(xiàn)淺海近距離目標測距測深，本文對于深海聲影區(qū)4條海底反射聲線形成的聲場干涉結(jié)構(gòu)進行近似推導，得到聲強在頻率域僅有兩種干涉周期，對于影區(qū)海面目標由于頻帶限制僅可觀測到第二種干涉結(jié)構(gòu)，頻域干涉周期隨距離單調(diào)增大，通過匹配聲場頻域干涉周期實現(xiàn)21 km遠距離海面目標的被動測距，利用船舶自動識別系統(tǒng)驗證了測距結(jié)果。另外，當影區(qū)內(nèi)存在多個聲源目標時，不同聲源形成的聲場干涉結(jié)構(gòu)相互疊加，限制本文所提方法的使用，有待進一步改善和發(fā)展這種方法。

圖17 實驗海域船舶分布與海面目標距離估計的對比(第592 s)Fig.17 The comparison between the practical distribution of vessels and the estimated results of target range

圖18 目標1和2的估計距離與實測距離的對比Fig.18 The comparison between the estimated range and the practical range of target 1 and 2

致謝：感謝參與海上實驗人員，是他們的辛勤勞動為本文提供可靠的實驗數(shù)據(jù)。

[1] 李啟虎. 聲吶信號處理引論[M]. 北京: 海洋出版社, 1985.

Li Qihu. The Introduction of Sonar Signal Processing[M]. Beijing: China Ocean Press, 1985.

[2] Smith J O, Abel J S. Closed-form least-square source location estimation from range-difference measurements[J]. IEEE Transactions on Acoustics, Speech, and Signal Processing, 1987, 35(12): 1661-1669.

[3] Baggeroer A B, Kuperman W A, Mikhalevsky P N. An overview of matched field methods in ocean acoustics[J]. IEEE Journal of Oceanic Engineering, 1993, 18(4): 401-424.

[4] Nardone S C, Lindgren A G, Gong K F. Fundamental properties and performance of conventional bearings-only target motion analysis[J]. IEEE Transactions on Automatic Control, 1984, 29(9): 775-787.

[5] Tao H, Hickman G, Krolik J L, et al. Single hydrophone passive localization of transiting acoustics sources[C]//Proceedings of OCEANS 2007，IEEE. Scotland, 2007：1-3.

[6] Kevin L C, Henrik S. Robust passive range estimation using the waveguide invariant[J]. Journal of the Acoustical Society of America, 2010, 127(5): 2780-2789.

[7] Rakotonarivo S T, Kuperman W A. Model-independent range localization of a moving source in shallow water[J]. Journal of the Acoustical Society of America, 2012, 132(4): 2218-2223.

[8] Lee S, Makris N C. The array invariant[J]. Journal of the Acoustical Society of America, 2006, 119(1): 336-351.

[9] Chuprov S D. Interference structure of a sound field in a layered ocean[M]//Ocean Acoustics: Current State. Moscow: Nauka, 1982: 71-79.

[10] Zhao Z D, Wang N, Gao D Z. Broadband source ranging in shallow water using the Ω-interference spectrum[J]. Chinese Physics Letters, 2010, 27(6): 110-113.

[11] 戚聿波, 周士弘, 張仁和, 等. 一種基于β-warping變換算子的被動聲源距離估計方法[J]. 物理學報, 2015, 64(7): 074301-1-074301-6.

Qi Yubo, Zhou Shihong, Zhang Renhe, et al. A passive source ranging method using the waveguide-invariant-warping operator[J]. Acta Physica Sinica, 2015, 64(7): 074301-1-074301-6.

[12] 任云, 戚聿波, 周士弘, 等. 聲強譜頻移補償?shù)牟▽Р蛔兞亢途嚯x估計方法[J]. 聲學學報, 2014, 39(3): 319-323.

Ren Yun, Qi Yubo, Zhou Shihong, et al. Waveguide invariant and range estimation based on frequency-shift-compensation of underwater acoustic intensity spectrograms[J]. Acta Acustica, 2014, 39(3): 319-323.

[13] Thode A M, Kuperman W A, D’Spain G L, et al. Localization using Bartlett matched-field processor sidelobes[J]. Journal of the Acoustical Society of America, 2000, 107(1): 278-286.

[14] 吳國清. 線譜非平穩(wěn)性分析和利用干涉譜測距[J]. 聲學學報, 2004, 29(4): 363-368.

Wu Guoqing. Nonstationarity analysis of line spectrum and range estimation by spectrum interference[J]. Acta Acustica, 2004, 29(4): 363-368.

[15] 翁晉寶, 李風華, 郭永剛. 典型深海聲場頻率-距離干涉結(jié)構(gòu)分析及實驗研究[J]. 聲學學報, 2016, 41(3): 330-342.

Weng Jinbao, Li Fenghua, Guo Yonggang. The sound field frequency-range interference patterns in deep water: theory and experiment[J]. Acta Acustica, 2016, 41(3): 330-342.

[16] 翁晉寶, 李風華, 郭永剛. 深海近距離聲場頻率-距離干涉結(jié)構(gòu)反演海底聲學參數(shù)[J]. 聲學學報, 2015, 40(2): 207-215.

Weng Jinbao, Li Fenghua, Guo Yonggang. Geoacoustic inversion based on near-field frequency-range interference pattern in deep water[J]. Acta Acustica, 2015, 40(2): 207-215.

[17] Porter M B. The BELLHOP Manual and User’s Guide: PRELIMINARY DRAFT[OL]. http://oalib.hlsresearch.com/Rays/HLS-2010-1.pdf, 2011-01-31/2016-01-07.

A passive sea-surface source ranging method by single hydrophone in deep sea

Weng Jinbao1，Yang Yanming1

(1.OceanLaboratoryofAcousticsandRemoteSensing,ThirdInstituteofOceanography,StateOceanicAdministration,Xiamen361005,China)

The sound field frequency-range interference patterns of sea-surface source in deep water are discussed based on ray theory. Approximate formulae of the sound field frequency-range interference patterns in the shadow zone are provided. In the shadow zone, the interference cycle decreases with the increase of the depth of the receiver and increases with the increase of the propagation range. As a result, the source range can be extracted from the interference cycle of the sound intensity using a single hydrophone. The sound field interference patterns caused by sea-surface broadband noise sources in the shadow zone have been observed in a deep sea experiment in the South China Sea. The availability of the shadow zone sound field interference structure in the passive source range estimation is verified by experimental data. Compared with conventional matched field passive source localization method, this method doesn’t need seafloor acoustic parameters and large scale computing of replica fields.

sea-surface source; shadow zone of deep water; sound field interference structure; passive ranging

10.3969/j.issn.0253-4193.2017.01.004

2016-01-07；

2016-07-13。

國家海洋局第三海洋研究所基本科研業(yè)務(wù)費專項資金(海三科2016016)；國家重點研發(fā)計劃(2016YFC1400103)；國家自然科學基金(61601132，41606116)；福建省自然科學基金(2016J01019)。

翁晉寶(1989—)，男，福建省晉江市人，從事海洋聲學研究。E-mail:wengjinbao@tio.org.cn

*通信作者：楊燕明，男，研究員，主要從事海洋聲學研究。E-mail:yangyanming@tio.org.cn

P733.23

A

0253-4193(2017)01-0036-10

翁晉寶，楊燕明. 深海海面目標單水聽器被動測距方法與驗證[J]. 海洋學報, 2017, 39(1): 36-45,

Weng Jinbao, Yang Yanming. A passive sea-surface source ranging method by single hydrophone in deep sea[J]. Haiyang Xuebao, 2017, 39(1): 36-45, doi: 10.3969/j.issn.0253-4193.2017.01.004