姚美娟　馬 力　鹿力成　郭圣明

（1中國科學(xué)院聲學(xué)研究所　北京　100190）

（2中國科學(xué)院水聲環(huán)境特性重點實驗室　北京　100190）

（3中國科學(xué)院大學(xué)　北京　100049）

利用寬帶引導(dǎo)聲源重構(gòu)拷貝聲場的目標(biāo)聲源定位?

姚美娟1，2，3?馬 力1，2鹿力成1，2郭圣明1，2

（1中國科學(xué)院聲學(xué)研究所北京100190）

（2中國科學(xué)院水聲環(huán)境特性重點實驗室北京100190）

（3中國科學(xué)院大學(xué)北京100049）

為了避開傳統(tǒng)的匹配場目標(biāo)定位技術(shù)對環(huán)境先驗知識的依賴性，提出了在均勻淺海環(huán)境中只知道少量環(huán)境參數(shù)的情況下，利用垂直接收陣和不同距離上的兩枚寬帶引導(dǎo)聲源重構(gòu)聲場對目標(biāo)聲源進行定位的一種方法。這種方法主要基于簡正波估計和聲場重構(gòu)兩種關(guān)鍵技術(shù)，同時省去了匹配場定位技術(shù)中大量的拷貝聲場計算。數(shù)值仿真主要采用線性Bartlett匹配處理器分析了目標(biāo)定位效果，在信噪比高于10 dB的情況下，定位效果良好。

引導(dǎo)聲源，簡正波估計，聲場重構(gòu)，目標(biāo)定位

1　引言

由于復(fù)雜的聲環(huán)境使得聲納技術(shù)受到明顯制約，目標(biāo)定位技術(shù)在水聲研究領(lǐng)域中受到越來越多的關(guān)注。傳統(tǒng)的匹配場處理技術(shù)（Matched field processing，MFP）出現(xiàn)于上世紀(jì)80年代，在目標(biāo)探測和三維定位等方面的應(yīng)用取得了很大進展，大大推動了被動聲納技術(shù)的發(fā)展。但匹配場技術(shù)對環(huán)境參數(shù)和聲場模型的依賴性極高，容易受到環(huán)境參數(shù)和聲場模型失配的影響。為了克服匹配場定位技術(shù)的不足，減少定位效果對海洋環(huán)境參數(shù)失配和聲場模型失配的影響，人們力求尋找一種新的不需要精確的海洋環(huán)境先驗知識的定位技術(shù)?！疤摂M接收”的概念于上世紀(jì)90年代末出現(xiàn)于水聲學(xué)中，它通過對垂直陣接收到的未知目標(biāo)聲源和已知引導(dǎo)聲源信號的相關(guān)處理，來得到在引導(dǎo)聲源位置處對未知聲源的虛擬接收聲場，從而消除原先垂直接收陣和未知聲源之間的海洋環(huán)境變化影響。Martin Siderius等人利用虛擬接收的概念，研究了淺海環(huán)境中的多途補償問題，適當(dāng)調(diào)節(jié)引導(dǎo)聲源的位置，可以在很大程度上消除垂直接收陣和虛擬接收陣之間的環(huán)境影響［1］。引導(dǎo)聲源的引入達到了對垂直接收陣和引導(dǎo)聲源之間的海洋環(huán)境未知補償效果。當(dāng)存在引導(dǎo)聲源時（在實際使用中，引導(dǎo)聲源一般選取為人為布放聲源或者機會聲源），可以有多種途徑來實現(xiàn)聲源定位：（1）根據(jù)垂直陣接收到的引導(dǎo)聲源聲場反演得到所需要的海洋環(huán)境參數(shù)，再利用聲場算法模型計算得到拷貝聲場，進行匹配聲場定位，同時給出目標(biāo)聲源的距離和深度估計；（2）利用文獻［1］等所給的虛擬接收技術(shù)，結(jié)合虛擬聲場在頻率-頻率偏移平面上的干涉條紋或者直接分析實際聲場在頻率-距離平面上的干涉條紋來獲得目標(biāo)聲源距離估計，這種方法原則上不能直接得到目標(biāo)聲源深度；本文中將提出一種新的基于簡正波分離技術(shù)和聲場重構(gòu)技術(shù)的目標(biāo)定位技術(shù)。利用聲場重構(gòu)技術(shù)與前面所述的與第一種方法相比，可以節(jié)約甚至避免拷貝聲場計算過程，與第二種方法（虛擬接收技術(shù)）相比，可以同時給出聲源距離和深度估計。

聲場重構(gòu)技術(shù)是一種新穎的水聲學(xué)技術(shù)，啟發(fā)于S.C.Walker［2］等人提出時間反轉(zhuǎn)聚焦在深度上的偏移，主要從聲場位置信息和聲場相位信息兩個方向重構(gòu)聲場。基于聲場重構(gòu)技術(shù)的目標(biāo)定位技術(shù)，首先從引導(dǎo)聲源接收聲場中估計簡正波模態(tài)函數(shù)。近年來，隨著人們對水聲信道認(rèn)識的加深和信息處理技術(shù)的發(fā)展，水聲簡正波分離技術(shù)得到了發(fā)展，Westwood等人提出了在沒有任何先驗信息的條件下，利用垂直線列陣（Vertical line array，VLA）測量到的聲壓數(shù)據(jù)，可以直接提取簡正波的信息［3］。Kuperman等人利用相同的方法，通過處理環(huán)境噪聲的VLA數(shù)據(jù)提取了簡正波的信息，并且提出了在已知聲速剖面的條件下，利用打靶法得到水平波數(shù)［4］。在提取得到簡正波模態(tài)函數(shù)的基礎(chǔ)上重構(gòu)得到拷貝聲場，再結(jié)合傳統(tǒng)的匹配場處理方法進行目標(biāo)聲源定位。

2　理論方法

2.1匹配場聲源定位

傳統(tǒng)的匹配場定位方法，是利用已知的（或反演獲得的）海洋環(huán)境參數(shù)信息，基于聲信號多途傳播特點采用傳播模型來計算拷貝聲場（向量），并與實際的陣接收聲場（向量）進行“匹配”，來實現(xiàn)對目標(biāo)聲源的定位。目前主要采用的匹配場處理算法，按照其權(quán)向量是否依賴測量數(shù)據(jù)，可以劃分為線性和自適應(yīng)匹配場處理兩大類［5］。垂直接收陣匹配場處理的線性Bartlett處理器的輸出為

其中R（rs,zs）是估計得到的目標(biāo)聲源陣接收聲場po的互譜密度矩陣，權(quán)系數(shù)v（r,z）取拷貝聲場向量p（r,z），上角標(biāo)符號H表示向量的共軛轉(zhuǎn)置，（r,z）是目標(biāo)聲源定位空間搜索點位置。在真實聲源位置（rs,zs）處D（r,z）取最大值，從而決定出目標(biāo)聲源位置?？梢钥吹?，對于傳統(tǒng)的匹配場聲源定位來說，關(guān)鍵性的工作就是如何從接收聲場中估計得到互譜矩陣以及如何計算拷貝聲場。

拷貝聲場計算，可以在海洋環(huán)境參數(shù)已知的情況下，利用合適的聲場算法模型來計算得到的，也可以利用接收到的已知聲源信號，來估計海洋環(huán)境參數(shù)或者是拷貝聲場計算直接可用的簡正波參數(shù)（包括簡正波本征值、本征函數(shù)和衰減系數(shù)等）。相比較而言，前者容易受到環(huán)境參數(shù)和聲場算法模型不準(zhǔn)確性的影響，引起失配問題，并且通過模型計算匹配聲場也有相當(dāng)大的計算量；后者直接從現(xiàn)場數(shù)據(jù)中提取所需要的環(huán)境參數(shù)信息或者是聲場計算所需要的簡正波參數(shù)，盡管也會存在參數(shù)提取不準(zhǔn)確性的失配可能，但可以很大程度上降低這種風(fēng)險，并且可以顯著降低匹配聲場計算量甚至于不需要直接計算。

2.2利用聲場重構(gòu)技術(shù)的拷貝聲場計算

聲場重構(gòu)技術(shù)是一種新穎的水聲學(xué)技術(shù)，利用聲場重構(gòu)技術(shù)的拷貝聲場計算以及進一步的目標(biāo)聲源定位，得益于S.C.Walker［2］等人提出的時間反轉(zhuǎn)聚焦在深度上調(diào)焦技術(shù)的啟發(fā)，類似于用測量的拷貝聲場來進行匹配場定位，具有較好的寬容性。

注意到在水平分層介質(zhì)中，垂直陣接收時拷貝聲場（向量）計算可以寫成

其中本征函數(shù)矩陣Φ和聲源位置信息的系數(shù)向量b分別為

zn是垂直接收陣陣元深度，km和δm分別為簡正波本征值和衰減系數(shù)?？梢灾溃灰軌驈膮⒖悸曉唇邮章晥鲋刑崛〉玫奖菊骱瘮?shù)矩陣Φ，再重構(gòu)出不同空間位置（r,z）處系數(shù)bm（r,z）或am（r），就可以方便的重構(gòu)得到所需要的拷貝聲場向量p（r,z）。

2.2.1簡正波本征函數(shù)獲取

在有引導(dǎo)聲源的情況下，從垂直陣接收到的寬帶引導(dǎo)聲源聲場中就可以提取得到對應(yīng)的簡正波本征函數(shù)。

對垂直接收陣接收得到的引導(dǎo)聲源時域聲壓信號序列作離散傅里葉變換，得到頻域的復(fù)聲壓信號，然后取某一窄帶（ω1：ωNω）范圍內(nèi)不同深度的復(fù)聲壓組成矩陣

Sg（ωj）和zg分別是引導(dǎo)聲源頻譜和深度，zi（i=1,2,···,Nz）對應(yīng)Nz個不同的垂直陣接收水聽器深度，km、δm和φm（z）分別是第m號簡正波的水平波數(shù)、衰減系數(shù)和本征函數(shù)。

簡正波本征函數(shù)φm（zi|ω）原則上也是頻率的函數(shù)。當(dāng)頻帶范圍（ω1：ωNω）比較窄時，簡正波本征函數(shù)變化不大，可以認(rèn)為是與頻率無關(guān)的常量，φm（zi）。這樣注意到簡正波本征函數(shù)的正交歸一化特性，它相當(dāng)于是矩陣Q中列向量的一組正交基函數(shù)。為了獲得這組正交基函數(shù)，對矩陣Q作SVD分解有

在垂直陣滿陣接收的情況下，V是與簡正波系數(shù)相關(guān)的矩陣，容易知道正交矩陣U的列向量與各號簡正波本征函數(shù)是成比例［6］，由此這樣我們就可以構(gòu)造出所需要的簡正波本征函數(shù)矩陣Φ。當(dāng)垂直陣不是滿陣時，由于簡正波本征函數(shù)在接收陣覆蓋深度范圍不是完全正交的，不能如（8）式這樣簡單的換算獲得，需要采用更為精致的方法。

2.2.2聲場重構(gòu)技術(shù)計算拷貝聲場

獲得簡正波本征函數(shù)矩陣Φ之后，再利用簡正波分解技術(shù)，可以從引導(dǎo)聲源接收聲場中提取得到對應(yīng)的簡正波系數(shù)

其中ω∈［ω1：ωNω］，

寬帶參考聲源一般可以選取所在海區(qū)的水面船只，距離和深度都是能夠獲得作為已知參數(shù)。這樣首先從提取得到的本征函數(shù)矩陣Φ，插值得到參考聲源深度和拷貝聲場計算深度上的本征函數(shù)值φm（zg）和φm（z），再結(jié)合（9）式對接收聲場的簡正波分解結(jié)果換算得到

注意到

利用兩個不同距離rg1，rg2上接收到的參考聲源聲場（暫時認(rèn)為參考聲源深度是一樣沒變化的），可以構(gòu)造出

以及

結(jié)合前面提取到的本征函數(shù)矩陣之后，不難寫出所需的拷貝聲場矢量

得出拷貝聲場向量后，代入（1）式就可得到定位結(jié)果。注意聲場重構(gòu)過程中的源信號譜級差異等方面引入的常系數(shù)比例系數(shù)差異，并不會影響到最后的匹配場定位結(jié)果（因為有歸一化處理）。

鑒于搜索點（r,z）空間網(wǎng)格點不宜太粗，這樣要求兩個引導(dǎo)聲源距離點差值Δr不要太大，要滿足網(wǎng)格點間距要求。

3　數(shù)值計算

3.1環(huán)境模型介紹

至少需要兩枚寬帶引導(dǎo)聲源對目標(biāo)聲源聲場進行重構(gòu)。如圖1是利用寬帶引導(dǎo)聲源重構(gòu)聲場對目標(biāo)聲源進行定位的水聲環(huán)境模型示意圖。垂直陣接收到來自于引導(dǎo)聲源1，2以及目標(biāo)聲源的聲信號，利用垂直陣接收到的引導(dǎo)聲源1，2的聲信號對目標(biāo)聲源處的聲場進行重構(gòu)。聲速剖面如圖1中所示，各水聲參數(shù)如表1中所示。

圖1　水聲環(huán)境模擬示意圖Fig.1　Sketch of underwater acoustic environment for simulation

表1　仿真環(huán)境水聲參數(shù)Tabel 1 Underwater acoustics parameters of the simulation

3.2數(shù)值計算結(jié)果分析

利用寬帶引導(dǎo)聲源進行SVD分解估計簡正波函數(shù)時，選用寬帶引導(dǎo)聲源的頻率段為210～250 Hz，對此寬帶聲場進行SVD分解，每隔1 Hz計算一次以估計簡正波模態(tài)函數(shù)；同時在中心頻率230 Hz下直接計算得到簡正波模態(tài)函數(shù)，將這兩者進行比較以驗證所估計的簡正波模態(tài)函數(shù)的準(zhǔn)確性。當(dāng)引導(dǎo)聲源深度處于某號簡正波的節(jié)點時，該號簡正波不能被激發(fā)出來，所以將引導(dǎo)聲源近海面布放。引導(dǎo)聲源1和引導(dǎo)聲源2布放的位置分別為（11 km，3.5 m）和（10 km，3.5 m），目標(biāo)聲源位置（18 km，18 m）。這里看到引導(dǎo)聲源1，2以及目標(biāo)聲源距離之間的關(guān)系為r-rg1=7?（rg1-rg2）；利用SVD分解方法估計簡正波本征函數(shù)時，聲場要滿足遠場條件，這里在較遠距離處（15 km，3.5 m）放置一枚寬帶引導(dǎo)聲源以估計簡正波模態(tài)函數(shù)。

如圖2，虛線的簡正波模態(tài)函數(shù)是直接計算的結(jié)果，實線的簡正波模態(tài)函數(shù)是利用寬帶引導(dǎo)聲源進行SVD分解的估計結(jié)果，實線與虛線幾乎完全吻合，完美的驗證了利用寬帶引導(dǎo)聲源進行SVD分解方法提取簡正波的理論。簡正波估計準(zhǔn)確，是后續(xù)聲場重構(gòu)以進行目標(biāo)聲源定位的前提。

圖2　未加噪聲時簡正波模態(tài)函數(shù)Fig.2 Normal mode function with no noise

利用提取的簡正波函數(shù)以及垂直陣接收到的引導(dǎo)聲源1，2的聲場信息的衰減和相位變化重構(gòu)聲場，再將重構(gòu)聲場與目標(biāo)聲源聲場匹配從而對目標(biāo)聲源進行定位。如圖3是線性Bartlett處理器輸出模糊度平面圖和surf立體圖。標(biāo)注“o”符號的位置為目標(biāo)聲源的真實位置（18 km，18 m），模糊度平面中標(biāo)注“+”符號的位置為線性Bartlett處理器的輸出峰值為（17.8 km，19 m），處理器輸出結(jié)果為水平距離定位誤差為1.11%，深度上定位誤差3.33%，定位效果良好。

圖3　未加噪聲時線性Bartlett處理器目標(biāo)定位結(jié)果Fig.3 Linear Bartlett processor localization results with no noise

在實際應(yīng)用中，垂直接收陣接收到的聲場都含有一定的噪聲場，下面將引導(dǎo)聲源與目標(biāo)聲源的聲場中均加入信噪比為15 dB的噪聲，分析目標(biāo)定位結(jié)果。如圖4，噪聲的干擾使SVD分解方法提取得到的簡正波模態(tài)函數(shù)與真實的模態(tài)函數(shù)值之間有一定的偏差。如圖5，目標(biāo)聲源真實位置為（18 km，18 m），自適應(yīng)匹配處理器的輸出峰值位置為（18.3 km，19 m），處理器輸出結(jié)果為水平距離定位誤差為1.67%，深度上定位誤差為3.33%，定位效果良好。

圖4　信噪比為15 dB簡正波模態(tài)函數(shù)Fig.4 Normal mode function with 15 dB signalto-noise-ratio

圖5　信噪比為15 dB時線性Bartlett處理器目標(biāo)定位結(jié)果Fig.5 Linear Bartlett processor localization results with 15 dB signal-to-noise-ratio

聲場中噪聲的加入對簡正波模態(tài)獲取和目標(biāo)定位都產(chǎn)生了一定的影響，將引導(dǎo)聲源與目標(biāo)聲源的聲場中的信噪比降低至10 dB，下面分析目標(biāo)定位結(jié)果。

如圖6，信噪比的降低使SVD分解方法提取得到的簡正波模態(tài)函數(shù)與真實的模態(tài)函數(shù)值之間的偏差進一步增大（相對于圖4）。因此由提取的簡正波模態(tài)函數(shù)重構(gòu)得到的聲場與真實聲場間的偏差也將進一步增大，這必然增加定位的誤差，如圖7，線性Bartlett處理器的輸出峰值位置為（18.4 km，20 m），處理器輸出結(jié)果為水平距離定位誤差為2.22%，深度上定位誤差為6.67%，信噪比的降低使得定位誤差變大。

圖6　信噪比為10 dB時簡正波模態(tài)函數(shù)Fig 6 Normal mode function with 10 dB signalto-noise-ratio

圖7　信噪比為10 dB時線性Bartlett處理器目標(biāo)定位結(jié)果Fig.7 Linear Bartlett processor localization results with 10 dB signal-to-noise-ratio

4　結(jié)論

利用寬帶引導(dǎo)聲源重構(gòu)拷貝聲場進行目標(biāo)定位的方法，主要基于簡正波獲取技術(shù)和聲場重構(gòu)技術(shù)兩種關(guān)鍵技術(shù)。它與傳統(tǒng)的匹配場定位技術(shù)相比具有可以省去匹配聲場建模和大量聲場計算的優(yōu)勢，與利用虛擬接收方法估計聲場干涉條紋的斜率對目標(biāo)測距的方法相比可以同時給出目標(biāo)聲源的距離和深度。本文通過仿真實驗，分析了在淺海海域中利用兩枚寬帶引導(dǎo)聲源進行目標(biāo)定位的效果。在信噪比為10 dB以上的情況下，水平距離定位誤差和深度定位誤差在5%到10%左右，定位效果良好。隨著信噪比的降低，定位誤差會增大，下一步工作將研究如何在較低信噪比下獲得較好的定位效果。另外，聲場重構(gòu)技術(shù)中相位信息的獲取在工程應(yīng)用中受到限制，這也是下一步工作的方向。

［1］SIDERIUS M，JACKSON D R，ROUSEFF D，et al.Multipath compensation in shallow water environments using a virtual receiver［J］.J.Acoust.Soc.Am.，1997，102（6）：3439-3449.

［2］WALKER S C，ROUX P，KUPERMAN W A.Focal depth shifting of a time reversal mirror in a rangeindependent waveguide［J］.J.Acoust.Soc.Am.，2005，118（3）：1341-1347.

［3］NEILSEN T B，WESTWOOD E K.Extraction of acoustic normal mode depth functions using vertical line array data［J］.J.Acoust.Soc.Am.，2002，111（2）：748-756.

［4］HURSKY P，HODGKISS W S，KUPERMAN W A. Extracting modal structure from vertical array ambient noise data in shallow water［J］.J.Acoust.Soc.Am.，1995，98（5）：2971.

［5］楊坤德.水聲陣列信號的匹配場處理［M］.西安：西北工業(yè)大學(xué)出版社，2008：165-166.

［6］ZHAO Zhendong，WANG Ning，GAO Dazhi，et al. Broadband source ranging in shallow water using the ?-interference spectrum［J］.Chin.Phys.Lett.，2010，27（6）：064301-1-064301-4.

Object source localization using the broadband guide source to reconstruct the replica field

YAO Meijuan1,2,3MA Li1,2LU Licheng1,2GUO Shengming1,2
（1 Institute of Acoustics，Chinese Academy of Sciences，Beijing 100190，China）
（2 Key Laboratory of Underwater Acoustics Environment，Chinese Academy of Sciences，Beijing 100190，China）
（3 University of Chinese Academy of Sciences，Beijing 100049，China）

A new method of object source localization，on the condition that a few environment parameters are known in range-independent shallow water environment，by using the vertical receiver array and the broadband guide source to reconstruct acoustic field，is presented to avoid the strong dependence on environment parameters of traditional matched field processing localization technique.This method is mainly based on normal modes estimation technique and acoustic field reconstruction technique and also can avoid the complicated computation of replica field for matched field processing localization.Linear Bartlett processor is used to analyze the simulation results.Desired results are obtained when the signal-noise ratio is higher than 10 dB. Key wordsGuide source，Normal mode estimation，Acoustic field reconstruction，Source localization

O427.1

A

1000-310X（2015）02-0135-07

10.11684/j.issn.1000-310X.2015.02.007

2014-05-08收稿；2014-09-15定稿

?國家自然科學(xué)基金項目（10774156）

姚美娟（1988-），女，山東臨沂人，博士研究生，研究方向：聲學(xué)。

E-mail：yyyshui@126.com