陳　發(fā)，李永勝，趙　罡，呂林夏，陳亞林

（1.中國船舶重工集團公司 第705研究所，陜西 西安，710075；2.水下信息與控制重點實驗室，陜西 西安，710075）

基于實際水文條件的點目標回波仿真

陳發(fā)1，2，李永勝1，趙罡1，呂林夏1，陳亞林1

（1.中國船舶重工集團公司 第705研究所，陜西 西安，710075；2.水下信息與控制重點實驗室，陜西 西安，710075）

為了分析水文條件對主動水聲探測系統(tǒng)信號的影響，首先利用射線聲學理論和BELLHOP模型，求解實測海洋水文數(shù)據(jù)下的本征聲線，然后選取主要聲線所攜帶的幅度、時延、相位等信息構建點目標回波信號模型。仿真分析給出了發(fā)射單頻正弦脈沖下的點目標回波仿真波形，并與實測波形進行了對比。計算結果表明，基于水文條件的點目標回波模型可以很好的模擬水文條件對水聲信號的影響。

水文條件；目標回波；本征聲線

0　引言

主動水聲探測系統(tǒng)通過發(fā)射信號，接收目標反射回來的信號探測水下目標［1］。由于目標信息包含在回波中，所以通過處理接收信號可以判斷目標的存在，利用回波信息還可以測定目標距離、速度、方位、航向等運動要素，但是水聲系統(tǒng)的性能會受到海洋環(huán)境的影響。水聲信道的低信息傳輸率、色散效應、界面影響、聲速剖面、以及多途等特性都會對接收信號產(chǎn)生影響。我國周邊海域多為淺海，海洋環(huán)境參數(shù)復雜多變。為了提高系統(tǒng)的環(huán)境適應性，有必要研究海洋信道對主動水聲探測信號的影響情況。

水聲信道對主動水聲探測信號的影響主要體現(xiàn)在2個方面：一是由于信道頻散和多途導致的波形衰減；二是參數(shù)的時變特性引起波形起伏。這些都會導致接收信號相對發(fā)射信號發(fā)生波形畸變和模糊，具體表現(xiàn)為幅度起伏、時頻擴展、去相關性。信道參數(shù)種類多樣，其中影響較大的是水文條件，也就是海洋中的聲速分布。

近年來，國內(nèi)外學者開展了海洋環(huán)境參數(shù)對水聲探測系統(tǒng)作用距離的影響研究［2］，方法多集中于求解經(jīng)典聲吶方程，對于海洋環(huán)境參數(shù)多使用平面波混響級RL和傳播損失TL代替。由于這2個參數(shù)是平均值、經(jīng)驗值，并不能反映信道的空變特性。文中利用海洋實測聲速分布數(shù)據(jù)，結合BELLHOP射線聲學模型，計算本征聲線［3］，然后選取合適的聲線信息（幅度、相位和延時），構建點目標回波，采用混合模型仿真混響，并疊加高斯白噪聲，從而得到接收信號。通過與實測數(shù)據(jù)對比分析，驗證了模型的正確性和可行性。

1　淺海聲傳播模型

1.1射線聲學

海水介質(zhì)中的聲傳播在數(shù)學上可以用波動方程的解來描述［4］，這個解是在特定問題的相應邊界條件和介質(zhì)條件下得到的。波動方程是聲壓p對坐標x，y z，和時間t的偏微分方程式中，c為聲速，且為坐標的函數(shù)，通常考慮為僅與深度有關。

有2種方法求解上述方程，一種是簡正波理論，利用名為簡正波的特征函數(shù)描述聲傳播，每一個特征函數(shù)都是方程的一個解，將這些簡正波疊加起來，以滿足邊界條件和聲源條件；另一種是射線聲學，利用聲線圖描述傳播圖像，聲場分布更加形象化。2種方法應用場合不同，前者通常用于低頻近場，后者只適用于高頻。由于目前的主動水聲探測系統(tǒng)使用的頻率通常為幾十kHz，所以通常使用聲線模型來求解探測信號聲場。

射線聲學的2個基本方程為

式（2）為強度方程，描述聲能量隨聲線傳播的情況；式（3）為程函方程，描述聲線的走向，該方程不僅可以導出聲線的方向，還可以求得聲線的軌跡和傳播時間。

1.2BELLHOP模型

射線聲學模型有很多便于計算機實現(xiàn)的算法，BELLHOP是其中較為常用的，它由Porter和Bucker從地震勘探常用的高斯波束追蹤法演變而來［5-6］，該算法解決了聲影區(qū)的問題，特別適合高頻、距離相關問題的求解。假設條件為點源在聲速僅與深度有關的柱形波導中。

圓柱坐標系下的標準射線方程為

其中：r=r（s），為弧長的函數(shù)；c（r，z）為聲速。

聲線的曲率和寬度由以下輔助方程推導而來

聲線定義為

加入初始條件后，可以求得聲場

其中，α為掠射角，該式利于計算機實現(xiàn)。輸入水文條件（聲速剖面、海底特性等）和求解參數(shù)，可以使用BELLHOP算法計算得到聲線圖。本征聲線為經(jīng)過發(fā)射點和接收點之間的聲線，該算法還可以計算每根本征聲線所攜帶的信息：幅度衰減系數(shù)、時延和相位變化量。

2　點目標回波模型

從大多數(shù)應用場合來看，水聲信道可以看成是緩慢時變的相干多途信道［7］。淺海試驗數(shù)據(jù)已經(jīng)表明，信道至少在20 s內(nèi)是穩(wěn)定的。由于海洋介質(zhì)中的聲速不是恒定的，加上界面的反射效應，聲波在海水中的傳播是多途的。不同路徑到達的信號相互干涉疊加形成接收信號。

聲波是微振幅波，滿足疊加定理，因此水聲信道對信號的影響可以使用線性網(wǎng)絡的系統(tǒng)函數(shù)來描述。假設聲傳播過程中不存在色散現(xiàn)象，系統(tǒng)函數(shù)可以描述為

其中多途路徑數(shù)為N，假設發(fā)射信號為s（t），則接收到的回波信號為

利用BELLHOP模型可計算出每根聲線的幅度衰減系數(shù)、時延和相位偏轉(zhuǎn)，故，式（10）簡化為

式中：AT為點目標強度；Ai，τi，θi為第i條路徑的幅度衰減系數(shù)、時延和相位偏移量，其中時延考慮了目標運動的影響。理論上能夠到達接收點的聲線是無數(shù)多的，為了計算的方便和實際情況的需要，可以只計算那些衰減量比較小的聲線。

同時到達接收陣列的還有混響和噪聲?；祉懯呛Ｑ笾猩⑸潴w的散射成分總和，具體可以劃分為體積混響和界面混響。噪聲包括海洋環(huán)境噪聲和系統(tǒng)自噪聲。那么接收信號可以表示為

其中：srev（t）為體積混響、海面混響、海底混響之和；snoise（t）為海洋環(huán)境噪聲。

3　仿真分析

3.1信道參數(shù)

這里采用海洋實測聲速剖面計算信道參數(shù)，海深H=150 m ，其相對海面聲速分布如圖1所示，可見該海域聲速呈負梯度分布。

設計4組仿真試驗，其基本參數(shù)為聲源深度Hs=70 m ，接收深度Hr=65 m ，發(fā)射單頻脈沖信號中心頻率為f0=20 kHz，采樣頻率fs=128 kHz ，計算波束數(shù)Num=101，垂直發(fā)射角為φ=±10o。假設海底是平坦的，底質(zhì)為淤泥粘土，其特性參考理論值［8］：海底聲速cs=1535 ms ，海底介質(zhì)密度ρb=143gcm3。每組試驗的探測距離和目標徑向速度，如表1所示。

下面給出第1組試驗的信道參數(shù)計算結果。由于入射聲場本征聲線較少，下面給出入射時的信道沖擊響應（實際上對應的就是每根本征聲線所攜帶的幅度衰減和時延信息）。

系統(tǒng)使用收發(fā)合置換能器，因此還需要計算從目標反射回去的本征聲線信息，才能構建目標回波。首先利用圖2中的信道沖擊響應，結合式（11）構建點目標反射波，以此合成信號作為發(fā)射信號，交換Hr和Hs，重置垂直發(fā)射角為φ=±20o，其他參數(shù)不變，求解本征聲線，見圖3。

圖1　淺海相對聲速剖面Fig.1　Profile of relative sound speed in shallow sea

表1　仿真參數(shù)Table 1　Parameters of simulation

圖2　聲源發(fā)射信號時的信道沖擊響應Fig.2　Impulse response of channel to sound source transmitting signal

圖3　目標反射到接收器的主要聲線Fig.3　Main sound ray reflected from target to receiver

從圖3可以看出，到達目標的聲線除了直達聲之外，還有經(jīng)海面和海底反射的聲線，接收信號為不同途徑到達信號的疊加。

反射信道沖擊響應如圖4所示?？梢娦诺罌_擊響應不僅與信道有關，還和換能器、目標深度、發(fā)射角有關。注意：以上信道參數(shù)，是在假設信道不存在頻散條件下求得的；程序中相位偏轉(zhuǎn)信息已經(jīng)通過換算到幅度衰減上。根據(jù)圖2和圖4中的信道沖擊響應參數(shù)，就可以利用模型完成點回波重建。

圖4　目標反射信號時的信道沖擊響應Fig.4　Impulse response of channel when target reflects signal

3.2回波仿真對比分析

下面通過比較發(fā)射信號、回波信號的特征，分析信道對信號的影響情況。因為實際海洋環(huán)境噪聲并非高斯白噪聲，淺海聲傳播會受到海底介質(zhì)特性的影響，但實際海底介質(zhì)特性與模型假設相差較大，所以仿真波形與實測波形在細節(jié)上存在差異。這里主要考慮信號的脈寬和頻域特性，并分析信號的去相關性。海上試驗示意見圖5。

圖5　海上試驗環(huán)境示意圖Fig.5　Schematic of sea trial environment

發(fā)射信號為單頻矩形脈沖，脈寬T=60 ms（信號持續(xù)時間均小于信道時間相關半徑，滿足模型假設條件）。發(fā)射信號如圖6所示。

因為水聽器實際接收到的信號還包括混響和噪聲，且混響為主要干擾。所以為了更好的與實測數(shù)據(jù)進行對比，還需要進行混響建模仿真。

混響伴隨著發(fā)射信號而產(chǎn)生，這里采用點散射和單元散射混合模型計算混響，模型假設條件為收發(fā)合置聲吶位置固定，聲波在介質(zhì)中沿直線傳播，忽略二次散射，信號為窄帶信號，計算得到的體積、界面混響及總混響如圖7所示?？梢娀祉懗掷m(xù)時間較長，且隨時間衰減很快。

圖6　發(fā)射矩形脈沖信號時頻波形圖Fig.6　Time-frequency waveform of transmitting rectangular pulse signal

圖7　仿真得到的體積、海面、海底混響及總混響波形Fig.7　Waveforms of volume，sea surface，sea bottom reverberation and total reverberation obtained by simulation

利用信道沖擊響應參數(shù)和式（11）構建點目標回波，這里考慮了目標的徑向運動帶來的多普勒頻移。將發(fā)射回波、混響、噪聲疊加起來，得到接收信號。由于發(fā)射信號載頻較高，且為窄帶特性，為了減少系統(tǒng)硬件和軟件的復雜性，可對信號的復包絡進行采樣。這里對信號先進行正交解調(diào)［9］，然后對2個正交分量進行采樣，采樣率為fs′=4 kHz，從而減少了數(shù)據(jù)量，得到的目標回波見圖8（這里只給出了第1組試驗的結果圖，其余組均只給出了分析結果）。為方便后面與實測回波進行對比，圖8中的時域波形只給出了目標反射回波部分；信號頻譜主要由回波和混響組成，且回波相對混響有頻移，頻移量約等于多普勒。

海上試驗接收到的信號時頻分布參見圖9。為了觀察的方便和對比的需要，圖中的時域圖只給出了目標回波出現(xiàn)的區(qū)間。從實測接收波形的頻譜可以看出，混響與回波的中心頻率之差基本上等于目標的多普勒，偏差來源于散射體、換能器的運動。實測與仿真混響頻譜相差較大，這是因為海底介質(zhì)特性采用了理論值，所以海底混響與仿真存在差異。對比實測波形和仿真回波，兩者在細節(jié)上存在較大差別，下面從回波脈寬、帶寬、中心頻率、互相關系數(shù)等參數(shù)進行研究。

圖8　仿真正交采樣回波時頻波形Fig.8　Simulated time-frequency waveform of orthogonal sample echo

圖9　實測正交采樣回波時頻波形Fig.9　Measured time-frequence waveform of orthogonal sample echo

實測和仿真回波參數(shù)相對誤差計算結果如表2所示。表中的頻率為正交采樣后的頻譜中心頻率，帶寬定義為峰值左右第1個零點之間的差值；求解發(fā)射信號與仿真回波之間互相關系數(shù)的積分時間為回波持續(xù)時間。從表2可以看出，信道對信號進行了去相關，這會影響信號的相關檢測性能。

表2　仿真與實測回波誤差分析Table 2　Error analysis between simulated and measured echoes

從表2可以看出，仿真得到的回波信號頻率誤差較小，而帶寬和脈寬誤差在15％以內(nèi)，誤差來自海底散射體的影響。利用水文條件構建的點目標回波模型，能夠很好地反映聲信道對發(fā)射信號的變換：幅度相位起伏、時域頻域擴展、去相關等。脈寬和帶寬計算結果表明，利用該模型來反映信道特性是可行的，該模型可以比較準確地反映信道特性。

4　結束語

為了探索水文條件對水聲探測系統(tǒng)信號的影響，文中采用實測淺海聲速數(shù)據(jù)，利用BELLHOP模型計算本征聲線，然后利用主要聲線信息重建點目標回波。將其與混響、噪聲疊加起來就得到水聽器的接收信號。

仿真結果表明，雖然海底、噪聲、混響等特性與仿真模型有較大差異，但通過分析信道對信號主要的影響因子，結果表明，該模型可以很好地反映信道特性。同時與實測波形的對比、計算結果驗證了該模型的可行性?；谒臈l件的點目標回波模型理論簡單，容易實現(xiàn)，可為今后進一步研究水文條件對水聲探測系統(tǒng)檢測性能的影響打下基礎。

［1］田坦，劉國枝，孫大軍.聲納技術［M］.哈爾濱：哈爾濱工程大學出版社，2006.

［2］何心怡，錢東，王光宇，等.水文條件對魚雷聲自導作用距離的影響［J］.魚雷技術，2007，15（5）：33-36.

［3］李永平，毛衛(wèi)寧.基于分層海洋射線傳播模型的接收信號仿真［J］.聲學與電子工程，2003（1）：12-15.

［4］Urick R J.水聲原理［M］.洪申，譯.哈爾濱：哈爾濱船舶工程學院出版社，1990.

［5］Alec J Duncan，Amos L Maggi.A Consistent，User Friendly Interface for Running a Variety of Underwater Acoustic Propagation Codes［J］.Proceeding of Acoustic，2006，20（22）：471-487.

［6］Poter M B，Bucker H P.Gaussian Beam Tracing for Computing Ocean Acoustic Fields［J］.The Journal of the Acoustical Society of America，1987，82（4）：1349-1359.

［7］惠俊英，生雪莉.水下聲信道［M］.北京：國防工業(yè)出版社，2007.

［8］劉伯勝，雷家煜.水聲學原理［M］.哈爾濱：哈爾濱工程大學，2010.

［9］潘仲明.隨機信號分析與最優(yōu)估計理論［M］.長沙：國防科技大學出版社，2007.

（責任編輯：許妍）

Simulation of Point Target Echo Based on the Actual
Hydrological Conditions

CHEN Fa1，2，LI Yong-sheng1，ZHAO Gang1，Lü Lin-xia1，CHEN Ya-lin1
（1.The 705th Research Institute，China Shipbuilding Industry Corporation，Xi′an 710075，China；2.Science and Technology on Underwater Information and Control Laboratory，Xi′an 710075，China）

To analyze the influences of hydrological conditions on the signal of an active hydroacoustic detection system，the ray acoustic theory and BELLHOP models are employed to solve the eigenray from measured oceanic hydrological data.A point target echo model is built with the information of amplitude，time-delay，and phase，which are carried by the main eigenray.The point target echo waveform is obtained from numerical simulation in the case of transmitting single-frequency sine impulse，and is compared with the measured waveform.Computational results show that the point target echo model based on hydrological conditions is applicable to simulate the influences of hydrological conditions on underwater sound signal.

hydrological conditions；target echo；eigenray

TJ630.34；TN911

A

1673-1948（2015）04-0311-05

2015-05-25；

2015-06-03.

陳發(fā)（1988-），男，在讀碩士，研究方向為水聲信號處理.