起伏海底掩埋目標聲散射特性數(shù)值仿真

于福建, 王 斌, 張培珍

?

起伏海底掩埋目標聲散射特性數(shù)值仿真

于福建1, 王 斌2, 張培珍3

(1.中國船舶工業(yè)系統(tǒng)工程研究院 水聲對抗技術重點實驗室, 北京, 100094; 2.上海交通大學 海洋工程國家重點實驗室, 上海, 200240; 3.廣東海洋大學 電子與信息工程學院, 廣東 湛江, 524088)

基于COMSOL多物理場平臺下的有限元數(shù)值計算方法, 完成收發(fā)分置條件下掩埋彈性目標聲散射數(shù)值計算。運用高斯譜法建立起伏海底模型, 給出粗糙海底起伏高度、相關長度、掠射角等參數(shù)對掩埋目標聲場特性的影響。結果表明, 在單一方向入射多點線陣上接收, 計算得到收發(fā)分置的信混比呈現(xiàn)起伏特性: 平坦或小起伏海底前向波形振蕩劇烈且呈現(xiàn)準周期性, 隨著海底粗糙程度增加, 起伏的規(guī)律性逐漸消失。當聲波垂直照射海底時, 海底起伏的相關長度變化對散射能量瑞利方差影響較大, 而幾乎不改變?nèi)鹄怠Ｐ比肷鋾r, 相關長度和起伏高度均是影響瑞利方差和均值的主要因素。研究成果為海底掩埋目標聲探測提供了參考依據(jù)。

起伏海底; 掩埋目標; 聲散射; 收發(fā)分置

0 引言

沉底水雷、海底管路、海底沉物等的探測和識別是水聲技術中的一大難題, 其原因是靜止目標往往和海底界面渾然一體, 目標回波與海底界面的強混響難以區(qū)分。非均勻海底的散射和海底透射是影響目標聲波散射特性的主要因素, 特別是當目標和界面粗糙度的均方高度以及波長可比擬時, 粗糙海底混響以及沉積層的衰減特性嚴重影響了目標的探測和識別。針對目標體散射和隨機粗糙面散射這2類目標, 已有的研究成果將2類散射問題分別單獨考慮, 或者將海底視為平坦界面條件下, 分析界面附近目標特性[1-3]。然而, 真實的海底不是絕對光滑的, 可以看成符合特定概率密度分布的隨機粗糙表面, 采用統(tǒng)計參數(shù)來描述, 如均方根高度、相關長度和分維數(shù)等。用于解決粗糙海底散射問題的經(jīng)典方法有基于Rayleigh假設的微擾法[4]、Kirchhoff近似法[5]、小斜率近似理論[6-7], 在流-固界面條件下將微擾理論和Kirhchoff近似結合的雙尺度方法[8]、消光定理[9], 以及雙基地海底散射模型仿真[10]等方法。至今為止, 針對平坦界面下目標散射問題取得了一定的進展, 用于解決非平滑海底界面附近目標高效目標散射的求解尚有待研究。文中基于COMSOL多物理場平臺, 構建海底和目標聯(lián)合的二維聲場模型, 解釋海底起伏對掩埋目標聲散射特性的影響。得到隨入射角、起伏參數(shù)變化的收發(fā)分置目標散射能量的統(tǒng)計變化規(guī)律。

1 掩埋目標聲散射數(shù)值模型

假定無限長柱體目標掩埋于海底, 柱體半徑為1 m, 入射聲波頻率為1 kHz, 計算聲場區(qū)域面積為60 m×30 m, 掩埋深度距離海底平面2 m。接收信號換能器位于距離海底4 m的線陣。二維掩埋彈性目標聲場模型如圖1所示。其中圖1(a)為幾何視圖, 圖1(b)為COMSOL平臺構建的聲固耦合模型。物理場中背景場、沙層及目標區(qū)域采用三角網(wǎng)格剖分, 最大網(wǎng)格為波長的1/6, 周圍完美匹配層有上下2種材質構成, 海底邊界上方為水, 下方為沙。采用6層的矩形掃掠網(wǎng)格。

考慮到指向性聲源可以大幅度提高信混比, 入射聲波采用錐形波[11], 其表達式為

其中

假定海底為非彈性且密度均勻介質, 目標及周圍介質聲學參數(shù)見表1。

表1 目標及周圍介質聲學參數(shù)

為了模擬更為真實的海底場景, 利用高斯譜法產(chǎn)生一維高斯型粗糙海底, 設置界面的起伏高度()和相關長度()。建立起伏海底下掩埋目標二維聲場模型。圖2給出了根據(jù)式(1)繪制的入射波背景場, 其中=11.25 m, 海底起伏參數(shù):=0.158 m;=0.5 m。

2 收發(fā)分置散射特性仿真分析

圖3~圖5依次給出垂直入射和入射掠射角為60°時的平坦海底, 弱起伏海底(起伏參數(shù):=0.158 m;=0.5 m)、強起伏海底(起伏參數(shù):=1.58 m;=5 m)條件下海底散射聲壓分布。

針對不同海底狀況, 依次改變?nèi)肷渎由浣菫?0°, 50°, 60°, 70°, 80°, 90°, 計算與海底垂直距離4 m的線陣(圖1)不同位置上的信混比, 如圖6所示?？梢钥闯? 當海底平坦或起伏較小時, 線陣不同位置上的信混比在前向方位的幅度起伏范圍相比后向更大, 且隨距離變化呈現(xiàn)準周期性波動, 掠射角越小, 信混比起伏周期越小振蕩越快。這一特性隨著海底粗糙程度的增加逐漸消失。

3 聲散射特性統(tǒng)計特性分析

固定掠射角, 隨機生成粗糙界面100次, 求解圖1所示線陣上目標散射強度。為了得到受起伏界面影響的目標聲散射特性的統(tǒng)計特性, 將分析線陣上各點能量之和的概率分布。即相關長度、起伏高度、掠射角固定, 可得到根據(jù)100次計算結果得到1條關于散射能量的概率分布密度函數(shù), 如圖7所示?？梢钥闯? 斜入射和垂直發(fā)射聲波, 得到的散射能量服從瑞利分布。當垂直照射海底時, 海底起伏的相關長度變化對瑞利方差影響較大, 而幾乎不改變?nèi)鹄?。斜入射時, 相關長度和起伏高度均是影響瑞利方差和均值的主要因素。

4 結束語

文中以二維聲場中彈性目標為例, 針對起伏海底下方目標聲散射特性的數(shù)值計算, 分析不同方向掠射角照射條件下的聲場分布特性, 并得出收發(fā)分置掩埋目標回波受粗糙界面影響產(chǎn)生的起伏規(guī)律。相應的成果可以拓展應用于復雜背景下掩埋目標的探測與識別。由于聲場計算區(qū)域較大, 特別是海底起伏加大了邊界處的網(wǎng)格剖分自由度的數(shù)量, 導致基于有限元的網(wǎng)格剖分技術計算量巨大, 難以在高頻、寬帶、三維目標散射特性預報的實際工程中實現(xiàn)。特別是隨著起伏的增強, 僅通過單一的散射數(shù)據(jù)不足以描述聲場特性, 此時需要整合目標散射、海底散射及面體耦合散射模型, 進一步對海底掩埋目標散射特征的統(tǒng)計分析。

[1] Zampolli M, Tesei A, Jensen F B, et al.A Computationally Efficient Finite Element Model with Perfectly Matched Layers Applied to Scattering from Axially Symmetric Objects[J].Journal of the Acoustical Society of America, 2007, 122(3): 1472-1485.

[2] Zampolli M, Tesei A, Canepa G, et al.Computing the Far Field Scattered or Radiated by Objects Inside Layered Fluid Media Using Approximate Green’s Functions[J].Journal of the Acoustical Society of America, 2008, 123 (6): 4051-4058.

[3] Zampolli M, Espana A L, Williams K L.Azimuthal Fourier Mode Decomposition of Generic Incident Fields in the Frequency Domain: Computational Model and Results[J].Journal of the Acoustical Society of America, 2011, 130(4): 2330-2330.

[4] Sot-CresPo J M, Nieto-Vesperinas M, Friberg A T.Scattering from Slighttly Rough Random Surfaces: A Detailed Study on the Validity of the Small Perturbation Method[J] Journal of the Optical Society of America.1990, 7(7): 1185-1201.

[5] Collaro A, Francesehetti G, Migliaccio M, et al.Gaussian Rough Surfaces and Kirchhoff Approximation[J].IEEE Transactions on Antennas and Propagation, 1999, 47(2): 392-398.

[6] 陳小泉, 馬忠成.小斜率近似方法分析粗糙界面聲散射問題[J].聲學技術, 2009, 28(6): 725-728.Chen Xiao-quan, Ma Zhong-cheng.The Small-Slope Approximation Specular Reflection Coefficient Rough Surface Acoustic Scattering[J].Technical Acoustics, 2009, 28(6): 725-728.

[7] Thorsos E, Broschat S.An Investigation of the Small Slope Approximation for Scattering from Rough Surfaces.Part I.Theory[J].Journal of the Acoustical Society of America, 1995 , 97(4) : 2082-2093

[8] 周士弘, 田玲愛.粗糙海洋界面散射對低頻聲傳播的影響分析[J].聲學技術, 2008, 27(5): 80-81.Zhou Shi-hong, Tian Ling-ai.Effect of Scattering from Rough Sea Surface and Bottom Boundaries on Low-Fre-que-ncy Sound Propagation[J].Technical Acoustics, 2008, 27(2): 80-81

[9] 郁高坤, 彭臨慧.聲波在粗糙液-液分界面上的散射[J].中國科學G輯, 2008, 38(2): 156-164.Yu Gao-kun, Peng Ling-hui.Scattering of Sound Wave from Rough Liquid-Liquid Interface[J].Science in China Series G-Physics, Mechanics & Astronomy, 2008, 38(2): 156-164.

[10] 李楠, 張明敏.雙基地海底散射模型仿真與試驗研究[J].海洋技術, 2009, 28(4): 48-51.Li nan, Zhang ming-min.Study on Simulation and Experiment of the Bistatic Bottom Scattering Model[J].Ocean Technology, 2009, 28(4): 48-51

[11] Isakson M J, Chotiros N P.Finite Element Modeling of Acoustic Scattering from Fluid and Elastic Rough Interfaces[J].IEEE Journal of Oceanic Engineering, 2014, 40(2): 1-10.

Numerical Simulation on Acoustic Scattering Characteristics of Targets Buried in Fluctuating Seabed

YU Fu-jian1, WANG Bin2, ZHANG Pei-zhen3

(1.Science and Technology on Underwater Acoustic Antagonizing Laboratory, Systems Engineering Research Institute of China State Shipbuilding Corporation, Beijing 100094, China; 2.State Key Laboratory of Ocean Engineering, Shanghai Jiao Tong University, Shanghai 200240, China; 3.College of Electronic and Information Engineering, Guangdong Ocean University, Zhanjiang 524088, China)

The acoustic scattering of buried elastic target under the bistatic condition is numerically calculated with the finite element method(FEM) on the COMSOL multi-physical field platform.A fluctuating seabed model is established by Gaussian spectrum method, and then the effects of fluctuation height, correlation length, grazing angle and other parameters on the acoustic characteristic of buried targets are obtained.The results show that the echo-to-reverberation ratio of the bistatic is fluctuant when sound enters in single direction while is received on a multi-point linear array.That is, for flat or gently fluctuating seabed, the forward waveform oscillates violently and presents quasi-periodicity, but the regularity of the echo-to-reverberation ratio fluctuation gradually disappears with the increase of seabed roughness.When a sound wave irradiates the seabed vertically, the variation of the correlation length of seabed fluctuation has a significant influence on the Rayleigh variance of the scattering energy, but has little effect on Rayleigh mean.If a sound wave enters obliquely, the correlation length and fluctuation height are the main factors affecting the Rayleigh variance and mean value.These results may provide a reference for acoustic detection of buried targets in seabed.

fluctuating seabed; ried target; acoustic scattering; bistatic

TJ617; O427.2

A

2096-3920(2018)06-0533-04

10.11993/j.issn.2096-3920.2018.06.004

2018-07-19;

2018-08-01.

國家重點研發(fā)計劃項目資助(2016YFC1400206); 水聲對抗技術重點實驗室項目資助(SSDKKFJJ-20160206).

于福建(1990-), 男, 碩士, 工程師, 主要研究方向為水聲工程.

于福建, 王斌, 張培珍, 等.起伏海底掩埋目標聲散射特性數(shù)值仿真[J].水下無人系統(tǒng)學報, 2018, 26(6): 533-536.

(責任編輯: 陳 曦)