姚美娟　馬 力　鹿力成　郭圣明

（1中國科學院聲學研究所　北京　100190）

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

（3中國科學院大學　北京　100049）

利用寬帶引導聲源重構拷貝聲場的目標聲源定位?

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

（1中國科學院聲學研究所北京100190）

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

（3中國科學院大學北京100049）

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

引導聲源，簡正波估計，聲場重構，目標定位

1　引言

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

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

2　理論方法

2.1匹配場聲源定位

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

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

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

2.2利用聲場重構技術的拷貝聲場計算

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

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

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

zn是垂直接收陣陣元深度，km和δm分別為簡正波本征值和衰減系數?？梢灾?，只要能夠從參考聲源接收聲場中提取得到本征函數矩陣Φ，再重構出不同空間位置（r,z）處系數bm（r,z）或am（r），就可以方便的重構得到所需要的拷貝聲場向量p（r,z）。

2.2.1簡正波本征函數獲取

在有引導聲源的情況下，從垂直陣接收到的寬帶引導聲源聲場中就可以提取得到對應的簡正波本征函數。

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

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

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

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

2.2.2聲場重構技術計算拷貝聲場

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

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

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

注意到

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

以及

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

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

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

3　數值計算

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

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

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

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

3.2數值計算結果分析

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

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

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

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

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

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

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

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

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

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

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

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

4　結論

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

［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è)大學出版社，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定稿

?國家自然科學基金項目（10774156）

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

E-mail：yyyshui@126.com