一種會聚區(qū)聲場基于參數(shù)化表征的方法研究?

呂琳 馮達(dá)

（91388部隊91分隊 湛江 524022）

1 引言

會聚區(qū)現(xiàn)象是一種特殊形式的深海聲道波導(dǎo)傳播引起的，根據(jù)Hale的描述［1］，從海面附近聲源發(fā)出的聲波在深海中折射并發(fā)生反轉(zhuǎn)，約在30～35海里以外的范圍折回海面，形成幾公里寬的環(huán)帶狀高聲強區(qū)域，即為會聚區(qū)。會聚區(qū)中的聲波具有較高的強度，這對于在大洋中的遠(yuǎn)程警戒、打擊具有十分重要的意義［2］。潛艇在深?；顒咏?jīng)常在具有會聚區(qū)形成條件的海區(qū)發(fā)現(xiàn)目標(biāo)，由于目標(biāo)方位變化緩慢，這給獲得目標(biāo)距離等其他信息帶來了很大的困難，無法確定目標(biāo)是否為會聚區(qū)目標(biāo)，無法進行有效的攻防決策。關(guān)于會聚區(qū)方面的研究，主要集中在會聚區(qū)環(huán)境特點、會聚區(qū)聲傳播理論等，針對會聚區(qū)目標(biāo)聲場特征進行會聚區(qū)目標(biāo)判定方面的研究較少，近年來有學(xué)者［4］利用聲場的干涉現(xiàn)象進行深海會聚區(qū)和影區(qū)的特性分析，對會聚區(qū)目標(biāo)的判定具有一定的指導(dǎo)意義。本文試圖從會聚區(qū)聲傳播特性的另一個角度入手，分析會聚區(qū)聲場特征，希望通過對聲場的參數(shù)化表征，獲得會聚區(qū)目標(biāo)判定的某些理論依據(jù)。

2 深海波導(dǎo)簡正波聲場特征

深海典型的聲速剖面為“三層結(jié)構(gòu)”型，存在一個聲速最小聲道軸，由于聲線被限制在海面和海底之間，在水中反轉(zhuǎn)時沒有能量損失，故能在深海中進行遠(yuǎn)距離傳播，形成會聚區(qū)現(xiàn)象。簡正波理論是分析和計算聲場的重要理論方法，在描述聲場特征時具有物理概念清晰，數(shù)學(xué)表達(dá)明確等優(yōu)點。

假設(shè)海水為水平無變化介質(zhì)，波動方程在海面和海底邊界條件下，遠(yuǎn)場解析解［5］為

其中p(r,z)表示接收處聲壓，z、zs分別為接收處和聲源深度，r為水平距離，Ψm為特征函數(shù)，krm為特征值。

實際深海波導(dǎo)中，水體中僅存在有限號簡正波，高號簡正波被”截止”［6］，式（1）重新表達(dá)為

其中，N表示簡正波最大號數(shù)，由W.K.B.模態(tài)條件確定。對于典型的深海聲道，A.O.WILLIAMS曾經(jīng)給出了簡單估算N的公式［7］：

假設(shè)聲源頻率為1KHz，水體中的簡正波個數(shù)N≈2400，因此深海中參與傳播的簡正波個數(shù)往往很多。這與淺海有很大不同。由于深海聲場簡正波個數(shù)較多，且不同簇的簡正波具有各不相同的干涉特征，在此首先將典型深海波導(dǎo)（不考慮表面聲道）的簡正波進行分族［8］。對于冬季波導(dǎo)，其常見聲速剖面如圖1所示?？筛鶕?jù)水平波數(shù)對簡正波分簇，當(dāng)聲源信號頻率為100Hz，根據(jù)圖1的參數(shù)，計算得到相鄰號簡正波的干涉跨度如圖2所示。

圖1 冬季典型深海聲速剖面

圖2 相鄰號簡正波的干涉跨度

根據(jù)圖2可知，中間號段“會聚模態(tài)”的干涉跨度為45km～55km，恰好與會聚區(qū)的距離間隔吻合。因此利用簡正波理論對深海會聚區(qū)聲場進行分析，可能為會聚區(qū)目標(biāo)的判定提供一些可靠的判據(jù)。

3 會聚區(qū)聲場的參數(shù)化設(shè)計

在分析會聚區(qū)聲場時，我們在此忽略不參與會聚的簡正波貢獻(xiàn)，令參與會聚的一簇簡正波初始號為m0族內(nèi)共有n0號簡正波，探討會聚區(qū)聲場在水平方向的干涉起伏可將式（2）做如下代換：

由式（3）求和號里面的表達(dá)式可以看出，聲場在距離r處的聲強起伏特征主要由水平方向的波數(shù)差以及深度方向上的本征函數(shù)決定。當(dāng)僅考慮水平方向上的干涉起伏時，聲壓振幅的起伏量主要由波數(shù)差影響，由于任意兩號簡正波存在波數(shù)差，總聲場表現(xiàn)為隨水平距離干涉起伏。由于會聚模態(tài)的相鄰兩號簡正波波數(shù)差隨號數(shù)為非均勻分布，導(dǎo)致在不同的會聚區(qū)內(nèi)各號簡正波干涉貢獻(xiàn)不同，將會表現(xiàn)出不同的特征差異［9～10］。我們擬用一曲線函數(shù)來描述水平波數(shù)差起伏，假定用簡正波號數(shù)為自變量擬合出的相鄰兩號簡正波的水平波數(shù)差為

則

整理后得到

其中

下面我們分別討論曲線函數(shù)中的系數(shù)a，b，c的貢獻(xiàn)，以此描述會聚區(qū)特征相似性和差異性。

假設(shè)式（4）中b=0，c=0，會聚區(qū)內(nèi)模態(tài)km+1-km=a，式（3）可寫為如下形式：

為討論問題方便，暫不考慮式（3）求和號前面的衰減系數(shù)，根據(jù)離散傅里葉變換的周期特性可知，得到水平距離上的聲強［11］分布如圖3所示。

圖3 會聚模態(tài)的波數(shù)差用a擬合后求的聲強分布

從圖3可以看出，當(dāng)只存在系數(shù)a時，會聚區(qū)表現(xiàn)為嚴(yán)格的周期特性，各個會聚區(qū)未出現(xiàn)展寬和差異。

假 設(shè) 式（4）中 c=0，會 聚 區(qū) 內(nèi) 模 態(tài)km+1-km=a+bm為m的一次函數(shù)。此時，式（3）可寫為如下形式：

相當(dāng)于在r域卷積了一項類高斯函數(shù)，這一項使會聚區(qū)展寬，且隨水平距離r增大，會聚區(qū)越來越寬。如圖4所示，由于m2項的存在，會聚區(qū)不再具備嚴(yán)格的周期性，表現(xiàn)為會聚區(qū)隨距離展寬并在會聚區(qū)內(nèi)存在起伏。

假設(shè)式（4）中會聚區(qū)內(nèi)模態(tài)km+1-km=a+bm+cm2為m的二次函數(shù)。此時，式（3）可寫為如下形式：

圖4 會聚模態(tài)的波數(shù)差用a和b擬合后求的聲強分布

m的二次項使會聚區(qū)展寬，m的三次項導(dǎo)致波形不對稱，出現(xiàn)峰側(cè)傾現(xiàn)象［12］。如圖5所示。

圖5 會聚模態(tài)的波數(shù)差用a、b、c擬合后求的聲強分布

以上分析表明，當(dāng)用a，b，c三參數(shù)來描述水平波數(shù)差的起伏，進而刻畫聲場的水平特征時，a，b，c三參數(shù)有明確的物理意義：a參數(shù)主要決定水平干涉跨度，b參數(shù)使會聚區(qū)展寬，c參數(shù)使每一會聚區(qū)在水平方向上不對稱，出現(xiàn)峰側(cè)傾現(xiàn)象。

在實際應(yīng)用中，如果可求解波數(shù)差曲線函數(shù)表達(dá)式中的a，b，c，能夠描述聲場水平特性，找出會聚區(qū)與非會聚區(qū)時聲場參數(shù)特征因子的差異，從而可為會聚區(qū)的判定提供參考標(biāo)準(zhǔn)。

4 聲場的參數(shù)化普適性分析

為進一步考核聲場的參數(shù)化普適性，現(xiàn)研究在不同水文條件下該參數(shù)化系數(shù)的穩(wěn)定性。這里我們僅僅給出在前述典型聲速剖面下的聲場計算。利用波數(shù)差曲線函數(shù)中的系數(shù)a，b，c，求出會聚模態(tài)下的簡正波水平波數(shù)，以此計算聲場，并與kraken計算得到的聲場分布比較，以此驗證聲場參數(shù)化的可行性，如圖6所示。

圖6 聲場參數(shù)化與kraken計算聲強分布對比

從圖6可以看出，用聲場參數(shù)化方法計算出的聲場傳播損失和用kraken計算出的結(jié)果吻合良好，這說明聲場參數(shù)化方法可很好地描述會聚區(qū)聲場特征。

5 結(jié)語

本文從另外一個角度研究了會聚區(qū)聲傳播特性，提出利用參數(shù)化表征會聚區(qū)聲場特征結(jié)構(gòu)，通過上述的理論分析和計算，采用聲場參數(shù)化特征因子a，b，c，可以很好地描述會聚區(qū)聲場特征，且具有明確的物理意義。在實際應(yīng)用中如能夠根據(jù)會聚區(qū)聲場特征，獲得特征因子a，b，c，表征會聚區(qū)與非會聚區(qū)聲場特征差異，可為會聚區(qū)目標(biāo)判定提供判定依據(jù)，這也是后續(xù)將要深入研究的問題。

參考文獻(xiàn)

［1］ F.E.Hale.Long Range Sound Propagation in the deep Ocean［J］.J.Acoust.Soc.Am.1961，33（4）：456-464.

［2］ R.J.Urick，G.R.Lund.Coherence of Convergence Zone Sound［J］.J.Acoust.Soc.Am.1968，43（4）：723-729.

［3］張仁和，何怡.大洋遠(yuǎn)程聲傳播理論及其應(yīng)用［J］.自然科學(xué)進展，1994，4（6）：670-676.

［4］李倩倩，李整林，張仁和.利用波導(dǎo)不變性解釋深海中的干涉現(xiàn)象［J］. 聲學(xué)技術(shù)，2011，30（3）：13-15.

［5］Finn B.Jensen，William A.Kuperman，Michael B.Por?ter，Henrik Schmidt.Computational Ocean Acoustics（Sec?ond Edition）［M］.New York：Springer，2011.

［6］張仁和.水下聲道中的反轉(zhuǎn)點匯聚區(qū)簡正波理論［J］.聲學(xué)學(xué)報，1980，1：28-41.

［7］A.O.Williams.Normal-mode methods in propagation of underwater sound［M］.New York： R.W.B.Stephens，1970.

［8］薛福茹.深海聲道會聚區(qū)的研究［J］.艦船科學(xué)技術(shù)，1982（1）：48-60.

［9］張晶晶，羅博.深海會聚區(qū)特征參數(shù)計算與分析［J］.聲學(xué)與電子工程，2017（3）：8-11.

［10］范培勤，笪良龍，李玉陽.深海會聚區(qū)特征參數(shù)計算方法研究［J］. 海洋技術(shù)學(xué)報，2012，31（4）：23-25.

［11］楊剛，韓樹平，黃銳.深海會聚區(qū)目標(biāo)強度分析［J］.艦船科學(xué)技術(shù)，2016，38（23）：131-134.

［12］張旭，張永剛.聲速垂直結(jié)構(gòu)變化引起的會聚區(qū)偏移［J］. 海洋科學(xué)進展，2010，28（3）：311-317.