目標被動測距的引導源方法研究*

李東東

(杭州應用聲學研究所 杭州 310023)

?

目標被動測距的引導源方法研究*

李東東

(杭州應用聲學研究所 杭州 310023)

基于波導不變量的原理，討論了一種可用于中程目標定位的引導源被動測距法。通過單水聽器提取的聲場信息，結合波導不變量理論和Hough變換的圖形處理方法，實現對目標的被動測距。

波導不變量; 引導源; 被動測距

Class Number TP391

1 引言

海洋環境是一個復雜的波導環境。海洋波導中不同階之間的簡正波，存在加強和減弱的干涉，即聲場存在干涉現象。淺海低頻聲場干涉結構被證明是精細及較穩定的[1]。

淺海環境中的聲源定位由于環境本身的復雜性變得十分困難。波導不變量原理最初由俄羅斯學者S.D.Chuprov[2]提出，用于描述相干性中的確定性成分，可以用于表征和解釋波導聲場中的干涉圖像和色散特性等波導環境信息，并可以反演聲場中聲源的位置信息。

目前計算淺海中波導不變量的方法有：頻移補償[3]、利用 LOFAR 譜圖的二維傅里葉變換脊[4]、利用 Hough 變換[5]，及深海環境下基于WKBZ[6]等。文獻[7]中將波導不變量和虛擬接收器結合，構建水平虛擬接收陣對目標聲源進行距離定位；文獻[8]中首次提出采用水平陣和引導聲源實現目標測距；文獻[9]針對實際海洋環境噪聲較復雜的情況,提出利用Gabor濾波進行降噪處理來改善目標聲源的定位效果；文獻[10]為了實現引導源目標定位算法在低信噪比聲源信號條件下的準確定位, 提出了利用領域平均法對算法進行降噪處理的方法；文獻[11]對于實現引導源目標定位算法在低信噪比聲源信號條件下的準確定位, 提出了利用領域平均法對算法進行降噪處理的方法。

本文主要是研究淺海聲場中波導不變量的提取方法，并討論了一種可用于中程被動定位的引導聲源被動測距法。該方法主要是基于單個水聽器，結合波導不變量以及Hough變換圖像處理的方法，實現被動目標的距離估計和跟蹤定位。

2 基本原理

2.1 波導不變量提取原理

低頻聲場具有穩定的干涉結構，在輻射連續譜的寬帶聲源的接收聲功率譜圖上，可以觀察到穩定的干涉條紋。在高頻情況，干涉圖可以通過射線理論進行預報；中、低頻時，有效的理論是簡正波理論，由于干涉條紋復雜，采用波導不變量進行分析較方便。

波導不變量，反映了距離、頻率和干涉條紋斜率的關系，描述了聲場的頻散特性和相消相長的干涉結構。在Pekeris波導中波導不變量β≈1。

波導不變量提取方法就是在低頻聲場上提取的LOFAR圖利用Hough變換，結合與波導不變量有關的干涉條紋方程提取波導不變量[12]。

如圖1所示。假定目標聲源作勻速直線運動，其線速度為v，接收水聽器位于坐標原點，目標聲源相對于接收點的最近通過距離為r0，目標聲源相對于接收點的最近通過時刻為t0。其中α為航向角(規定與x軸正方向的夾角)，θ為目標方位角。

圖1 目標運動幾何關系圖

首先推導第一種干涉條紋公式：

由目標運動幾何關系可知：

(1)

(2)

對式(2)求導:

(3)

由波導不變量的定義可知：

(4)

為關于頻率f的方程：

(5)

(6)

由上面可以得到：

(7)

整理積分后，得到：

(8)

若是已知干涉條紋頂點頻率f0，則f=f(τ)，f0=f(0)，對上式取反對數可得到：

(9)

推廣到一般情況，則有：

(10)

式(10)推導出的是第一種干涉條紋方程，表示的干涉條紋是一族廣義雙曲線，該方程只適用于干涉條紋頂點已知的情況。若是干涉條紋未出現頂點，則上述方程式不再適用。下面推導不受頂點限制的第二種干涉條紋方程。

第二種干涉條紋公式：

(11)

由圖1目標運動的幾何關系可知：

(12)

(13)

即第二種干涉條紋方程。若條紋有干涉頂點，則在過頂點的情況下，第二種干涉條紋方程為

(14)

在過頂點條件下，第二種干涉條紋方程與第一種干涉條紋方程一樣，都表示一簇廣義雙曲線；對于遠場目標的LOFAR圖，第二種干涉條紋方程則近似表示一簇直線。

2.2 波導不變量提取仿真

仿真條件：采用Pekeris模型，海深55m，聲速1500m/s，海水密度為1024kg/m3，海底介質聲速1749m/s，其密度1941kg/m3，海底無吸收。接收水聽器深度：30m。目標等深航行，聲源深度：4m，航行速度：12m/s，航向角為30°,最近通過距離：1200m，選取時間段0s～300s，最近通過時間為觀測時刻150s。處理頻段：300Hz～1000Hz。采用KRAKENC進行簡正波模型的聲場建模[13]。仿真流程圖如圖2所示。

圖2 參數估計仿真流程圖

首先是分析最近通過距離的情況，根據建模后的聲場得到目標信號的LOFAR圖如圖3所示。按照第一種干涉條紋方程，進行Hough變換，得到圖4。提取出的信道波導不變量為0.97，目標速度與接收水聽器距艦船水平最近距離比v/r0為0.0100。

圖3 目標信號的LOFAR圖

圖4 Hough變換(第一種干涉條紋方程)

按照第二種干涉條紋方程進行Hough變換，結果如圖5。提取信道波導不變量β為0.96，航向角為30°。

圖5 Hough變換(第二種干涉條紋方程)

按照已知頂點條件下第二種干涉條紋方程進行Hough變換，得到圖6。提取信道波導不變量β為0.96，航向角為30°。與Pekeris模型中理論指出的β≈1相比，β偏差為4%。

圖6 Hough變換(已知頂點條件下第二種干涉條紋方程)

若是在沒有頂點通過的情況下，首先是在沒有通過最近距離的情況下得到的聲場的LOFAR圖以及Hough變換結果，如圖7所示。提取波導不變量約為0.98，與在Pekeris模型中理論值相比，偏差為2%，航向角為30°。

圖7 未知頂點情況下按照第一種干涉條紋方程仿真圖

3 引導源被動定位測距原理

3.1 基本原理

美國物理學者Aaron[14]將提出了一種新的水下定位方法，該算法不需要知道確切的海洋環境的具體參數，而是根據接收到的引導聲源信號，與接收的目標聲源信號進行處理，從而求得目標聲源的距離，有效地克服了環境的時變性的影響。

本節通過水聽器接收到來自引導聲源和目標聲源的聲場信號，將得到的LOFAR圖進行Hough變換，提取出聲場干涉條紋斜率，帶入到推導出的測距公式中，求得目標聲源的距離。

圖8中給出了引導源被動測距聲源和接收器的位置關系圖。假設海水深度為H，矢量水聽器的深度為h，目標聲源深度為z0，且目標聲源距離接收水聽器的水平距離為r1，引導聲源的深度以及距接收水聽器的水平距離分別為zg和rs。

圖8 引導聲源被動測距的空間布放圖

在引導聲源的LOFAR圖上的干涉條紋滿足方程：

(15)

其中，fs表示引導聲源的頻率，rs表示引導聲源相對接收點的距離，fgs為引導聲源的參考頻率，rgs為引導聲源相對接收點的參考距離，βs為波導不變量。

同樣，目標聲源的干涉條紋方程為：

(16)

此外，根據波導不變量理論可知：

(17)

假設引導源和目標所處的海洋環境相同，則可以認為：

βs≈β1

(18)

則根據式(17)和式(18)可得：

(19)

其中，dfs/drs和df1/dr1分別為引導源和目標干涉條紋的斜率。

在同一海區，假設海底是平整的前提下，利用引導聲源的聲場LOFAR圖估計出干涉條紋的斜率和波導不變量，結合目標聲源的LOFAR圖，按照式(19)，即可以測得運動目標的距離。

3.2 仿真研究

首先利用單水聽器接收一定時間內的引導聲源的寬帶噪聲，獲得引導聲源的LOFAR圖。根據干涉條紋方程，采用Hough變換對引導聲源的LOFAR圖進行處理，分別提取出干涉條紋斜率和波導不變量。然后再進一步結合目標聲源的LOFAR圖，估計目標聲源的距離。

仿真條件：設海平面為0m，海深為55m，引導聲源的深度為4m，引導聲源以10m/s的速度航行，接收水聽器距離引導聲源的水平最近距離P=10m，假設引導聲源的位置為[-1000m:10m:-400m]，負數表示目標向接收水聽器運動。目標聲源深度為4m，目標源以15m/s的速度航行，接收水聽器距離目標聲源的水平最近距離P=100m，設目標聲源的位置為[-3600m:15m:-2700m]。接收水聽器的深度為30m，頻率下限為300Hz，頻率上限為2000Hz，采樣率為10kHz。

圖9 引導聲源的LOFAR圖

圖10 目標聲源LOFAR圖

圖11 目標聲源距離隨時間變化圖

已知目標聲源做勻速運動，以20s處的距離為例，此時目標的真實距離是3300m，目標的估計距離是3386m，誤差為2.16%。

若是將仿真條件中的引導源深度改為5m，其他條件不變，仿真引導聲源被動測距如圖12所示。同樣以20秒處為例，目標真實距離是3300m，估計距離是3475m，誤差為5.23%。

圖12 目標聲源位置隨時間的變化圖

4 結語

文中主要研究了基于波導不變量的引導源被動測距。首先是由兩種干涉條紋方程，結合Hough變換提取出聲場中的波導不變量。然后提出一種基于單水聽器通過接收引導聲源和目標聲源的聲場信息，結合波導不變量的知識，就可以實現在海洋環境未知的情況下估計出目標的距離，并進行了真實值和估計值對比驗證。

[1] 凌青,生雪莉,袁延藝,等.運動目標低頻矢量寬帶譜干涉結構及其應用研究[J].聲學技術,2012,31(4): 366-369.

[2] S.D.Chuprov. Interference structure of a Sound Field in a Layered Ocean[M]. Acoustics of the Ocean: Current State, Edited by L.M. Brekhovskih and I.B. Andreevoi, Moscow,1982:71-91.

[3] 蘇曉星,張仁和,李風華.用頻移補償方法計算聲場的波導不變量[J].聲學技術,2007,26(6):1073-1076.

[4] 安良,王志強,陸佶人.利用LOFAR譜圖的二維傅里葉變換脊計算波導不變量[J].電子與信息學報,2008,30(12):2930-2933.

[5] 田愛玲,劉福臣,周士弘.利用Hough變換提取波導不變量[J].聲學與電子工程,2009(4):22-24.

[6] 崔寶龍,笪良龍,徐國軍,等.基于WKBZ方法計算深海環境下的波導不變量[J].聲學技術,2013,32(4):426-430.

[7] 鹿力成.海底變化環境下目標距離定位[J].水聲及物理聲學,2007,26(4):65-68.

[8] 趙振東,王大治,王好忠,等.波導不變量原理在目標測距中的應用[J].聲學技術,2009,28(2):45-46.

[9] 徐國軍,笪良龍,張林,等.基于引導源目標定位算法的降噪處理研究[J].聲學技術,2009,28(5):97-98.

[10] 徐國軍,笪良龍,李玉陽,等.低信噪比條件下引導源目標定位算法[J].聲學技術,2010,29(3):336-339.

[10] 何青海,笪良龍,韓梅,等.聲場干涉條紋圖像的降噪研究與性能仿真[J].火力與指揮控制,2012,37(6):130-133.

[11] 余赟.淺海低頻聲場干涉結構及其應用研究[D].哈爾濱:哈爾濱工程大學,2010:124-130.

[12] 溫韶娟.中近程目標定位方法研究[D].哈爾濱:哈爾濱工程大學,2011:32-35.

[13] MICHAEL B P.The KRAKEN normal mode program (DRAFT)[M]. LaSpezia: SACLANT Undersea Research Centre,1991.

[14] Aaron M Thode. Source ranging with minimal environmental information using a virtual receiver and waveguide invariant theory[J]. J.A.S.A,2000,108:1582-1594.

Guide Source of Target Passive Ranging

LI Dongdong

(Hangzhou Research Institute of Applied Acoustics, Hangzhou 310023)

Based on the principle of the waveguide invariant, a method for passive ranging which can be used to the source medium range target positioning is discussed. Combined the acoustic field information extracted by a single hydrophone with the theory of the waveguide invariant and the method of Hough transform, the target passive location is realized.

waveguide invariant, guide source, passive ranging

2016年5月10日,

2016年6月16日

李東東,男,碩士研究生,助理工程師,研究方向：水聲信號處理。

TP391

10.3969/j.issn.1672-9730.2016.11.033