基于簡正波模態(tài)的淺海聲傳播的最佳深度規(guī)律研究

竇雨芮，周其斗，譚 路

(海軍工程大學 艦船與海洋學院，湖北 武漢 430033)

0 引 言

在水下探測和通信中，聲波有著不可替代的優(yōu)勢[1]，而聲吶通過對聲信號進行轉換和處理實現對水下目標的偵察和定位，聲吶能否有效發(fā)揮其作用不僅與聲吶自身性能有著密不可分的關系[2]，與聲吶所處的海洋環(huán)境同樣息息相關。聲吶所在區(qū)域的海洋環(huán)境作為聲信號的載體對聲波在其中傳播有較大的影響。Jensen[3]關于海洋環(huán)境對聲信號的影響多集中在水平方向上聲能量的衰減，對深度方向上聲能量分布有較少研究。由于海洋聲信道上下分別受到海面和海底的制約，其邊界條件的不相同，導致聲能量在深度方向上的分布并不均勻，當傳播距離一定時，不同接收點深度的聲傳播損失差異可達20 dB，因此研究不同條件下深度方向上的聲傳播損失最小值的位置即最佳深度對聲吶探測和水聲通信有一定幫助。

國內外學者對最佳深度做了一定研究：Weston[4-5]對聲場的平均聲強進行了研究，通過計算等聲速梯度條件下深度方向上的距離平均聲強，得出聲強在聲源深度和聲源對稱深度存在峰值，且峰值比平均聲強大50%；Gershfeld[6]通過對聲傳播損失和信噪比的計算得出淺海聲傳播的最佳深度與聲源點深度相同的結論，并研究了聲源位置、海底類型、聲速梯度和聲源頻率對最佳深度的影響；Ferla[7]基于前人結論對聲源深度未知時的最佳深度位置進行了研究；王曉宇等[8]分別對理想Pekeris波導和實際淺海分層波導情況下的水平線列陣最佳布放深度進行研究，得出最佳布放深度應與聲源深度相同。

現有的研究[9-10]只針對最佳深度的現象進行了歸納和總結，并沒有探索影響最佳深度的根本原因，且得出的結論存在一定的局限性。本文對影響最佳深度的因素進行了研究，通過對理想液體波導中簡正波階數對最佳深度影響的研究，得出簡正波模態(tài)疊加是產生最佳深度的根本原因，并較為全面地總結了不同環(huán)境下最佳深度的規(guī)律。

1 理論公式

淺海聲傳播問題可以近似采用單位強度的諧和點聲源在水平均勻的分層波導中的響應來模擬，任意一點的聲壓p滿足如下波動方程[11]：

式中：p(r,z,t)為與水平距離r、深度z和時間t相關的聲壓；c(z),ρ(z)分別為與深度相關的聲速和密度，s(t)代表點聲源強度。

在柱坐標系下上述方程轉化為：

計算區(qū)域的邊界條件、海水深度和c(z)均和r沒有關系，因此可以進行如下設置：p(r,z)=Φ(R)Ψ(Z)，并將此表達式代入上述方程整理后得：

上述方程加號左右兩側分別是只關于r和z的函數，因此使等式滿足的唯一辦法就是方括號內的值同時為常數，選取這一常數值為進行分離：

由于海面即z=0被設置為壓力釋放邊界，海底z=D在足夠深處被設置為絕對硬邊界，則函數Ψ 同時滿足：

Ψm(z)表示當常數值為krm時的函數。

利用Ψm(z)的正交性，在對Ψm(z)進行歸一化后聲壓的任意解可以表示為解空間里模式函數之和，即聲壓可以用以下形式表示：

將上述表達式代入聲壓波動方程中并進行簡化得到：

運用模式函數的正交性質得到：

上述方程為標準貝塞爾方程，其解采用漢克爾函數形式：

由于聲壓解滿足無窮遠輻射條件，所以解選用第一類漢克爾函數，省去時間項 e xp(-iwt)，并用漢克爾函數展開式代入得到：

在海洋聲場計算中，聲壓并不總是有意義的，當聲壓為復數時，其表示壓力幅值和相位，為更清楚地描述聲場，采用聲傳播損失描述：

其中：p0為在無窮大的自由空間中單位強度的點聲源產生的聲壓，。

將其代入聲傳播損失表達式中：

上式的傳播損失是基于不同模態(tài)簡正波之間相位疊加的結果，如果忽略簡正波之間的相位差異，只考慮簡正波攜帶的能量多少，通過這種方式疊加則產生非相干損失。Jensen通過對比不同模態(tài)的非相干疊加與1/3倍頻層平均后的聲傳播損失測量結果，發(fā)現非相干疊加的聲傳播損失計算結果和1/3倍頻層平均的測量結果一致。驗證說明了非相干疊加的計算結果在描述傳播規(guī)律時可以代替相干疊加。下面給出非相干疊加的計算公式：

2 最佳深度的影響因素研究

2.1 模型參數

選取如圖1所示的海洋環(huán)境進行計算，單位強度的諧和點源位于海水層z0處，海水和海底層密度和聲速均保持不變，海底為半無限液體海底，聲信號在傳播過程中分別受到海面和海底邊界條件的影響，海底參數的選取會對聲傳播有一定影響，本文采用Hamilton[12-13]總結的大陸架附近的海底參數進行計算，部分參數選取如表1所示。

圖1 Pekeris波導示意圖Fig.1 A Pekeris waveguide with fluid bottom

表1 不同類型的海底參數Tab.1 Parameters for different basement types

聲信號在實際海水中傳播時會產生能量損失，這種能量損失不僅包括擴展損失，還包括海底衰減損失和海水衰減損失。試驗證明，聲能量在海底衰減的多少與頻率相關，且頻率越高，能量衰減越大，不同類型海底與頻率相關的衰減系數見表1。海水中的衰減系數與頻率同樣存在一定關系，本文采用Thorp[14]總結的海水衰減系數與頻率的關系進行計算。

2.2 模型計算

淺海聲場比深海聲場有著更為復雜的傳播狀況，聲波在其中傳播受到多種因素的影響，本文主要考慮聲源頻率、聲源深度、海水深度和海底參數對最佳深度的影響。

2.1.1 聲源參數對最佳深度的影響

不同頻率的點聲源在海水中激發(fā)的聲場不同，同一頻率的聲源在不同位置也會激發(fā)不同的聲場。為研究聲源參數對最佳深度的影響，本文選取點聲源頻率Freq 分別為 100 Hz和 800 Hz，海水深度為 100 m，海底參數為泥，分別計算放置于10，30和50 m處的點聲源激發(fā)聲場的非相干聲傳播損失場，結果如圖2所示。

從上述現象可以得出：聲源頻率和聲源深度對最佳深度有很大影響，在低頻時，無論聲源深度為多少，在近場和遠場時最佳深度始終在海水層中部，在高頻時，最佳深度位于聲源深度和聲源的對稱深度，且頻率越大，距離聲源點水平距離越近，該最佳深度現象就越明顯。

2.2.2 海水參數對最佳深度的影響

海水參數中聲速剖面的存在對聲場有很大影響，張旭芝[15]已經對此做了一定研究，本文主要研究等聲速情況下的最佳深度。為研究海水深度對最佳深度的影響，選取海水深度為100，150和200 m，聲源頻率為100 Hz，聲源點深度分別設置為10，15和20 m，其余參數保持不變，得到聲傳播損失隨距離和深度變化如圖3所示。

從圖3可以看出：海水深度不同導致靠近海面處放置的點聲源激發(fā)的聲場也不同，水深100 m，聲源深度10 m條件下，聲傳播損失隨距離和深度變化的偽彩圖在近場和遠場均呈現為弧形，該現象說明海水層中間處為最佳深度；水深為150 m，聲源深度15 m條件下，近場的聲傳播損失出現弓形現象，說明在近場處最佳深度位于聲源深度和聲源的對稱深度時，隨著傳播距離的增大，最佳深度向中間移動，傳播至遠場變?yōu)榛⌒?，最佳深度位于海水中間處；水深200 m，聲源深度20 m時，在近場和遠場均出現弓形現象，隨著傳播距離的增大，弓形曲線的峰值向海水中部移動。

圖2 聲源深度分別為 10、30 和 50 m 的聲傳播損失場Fig.2 Transmission loss field for fixed sources at 10,30,and 50 m depth

從上述現象可以得出海水深度對最佳深度有很大影響，海水深度越大，近場處的最佳深度位于聲源深度和聲源的對稱深度的現象就越明顯，當聲傳播較遠距離后，最佳深度由聲源深度和其對稱深度向海水中部移動。

2.2.3 海底參數對最佳深度的影響

淺海邊界條件對聲傳播的影響較深海更為顯著，不同海底參數下的聲傳播也不同，本文選取表1中泥、極細砂和細砂三類海底進行計算，研究海底參數對最佳深度的影響。聲源頻率為100 Hz，聲源深度為10 m，其余參數均與圖1一致。得到聲傳播損失偽彩圖如圖4所示。

由圖4可知：隨著海底參數的變化，最佳深度的現象變化相對較小，海底類型為細砂時才出現相對明顯的最佳深度位于聲源深度和聲源的對稱深度的現象?？梢缘贸龊５讌祵ψ罴焉疃扔绊懴鄬^小。

本節(jié)研究聲源參數、海水參數和海底參數對最佳深度的影響，從上述結論初步得出聲源頻率、聲源深度和海水深度對最佳深度的影響最大，與Gershfeld和Ferla得出結論不同的是，最佳深度并不總是出現在聲源深度和聲源的對稱深度，在頻率較低、海水深度較小情況下，最佳深度出現在海水層中間處。

圖3 海水深度分別為 100，150 和 200 m 時的聲傳播損失場Fig.3 Transmission loss field for water depth of 100,150, and 200 m

3 簡正波模態(tài)主導的最佳深度研究

由波動理論可知，海水深度和聲源頻率是影響海水中簡正波階數的直接因素，因此為探究影響最佳深度的根本原因，本文研究了由海水深度和聲源頻率主導的簡正波模態(tài)階數對最佳深度的影響，通過對理想液體波導（見圖5）中簡正波階數對最佳深度影響的研究來驗證簡正波疊加是產生最佳深度的根本原因，并總結最佳深度的相關規(guī)律。

3.1 理想液體波導中的最佳深度研究

理想液體波導是簡化了的實際海洋環(huán)境，在波導中海面和海底均設置為全反射邊界，聲波在界面處只有反射沒有透射，只存在水平擴展損失。這種邊界條件的設置對遠距離傳播的情況是合理的，這是因為對遠場起主要貢獻的是掠射角較小的聲波，而小掠射角在實際海底處的反射類似全反射邊界處的反射。

圖4 海底類型分別為泥、極細砂和細砂時的聲傳播損失場Fig.4 Transmission loss field for basement type as silt, very fine sand and fine sand

圖5 理想液體波導示意圖Fig.5 Idealized ocean waveguide model

3.1.1 理想液體波導中的簡正波

理想液體波導中的聲場是由各階簡正波疊加得到的，每一階簡正波的幅值與海水深度密切相關，為得到各階簡正波與深度相關的位移幅值，下面對理想液體波導中的位移勢函數的表達式進行推導，其計算模型和參數如圖6所示。

圖6 理想波導中的前5階模態(tài)隨深度變化圖Fig.6 Depth dependence of the first 5 normal modes in ideal waveguide

由波動理論可知：聲場的位移勢函數 ψ 同樣滿足波動方程，在圖6所示的柱坐標系下的理想液體波導中，其滿足以下的波動方程：

其中Sw為點聲源強度，運用積分變化法得出總聲場的積分表達式為：

則kr可以表示為。

理想波導的位移勢函數的無限積分形式可以轉換為無限個留數之和：

kzm代表某階簡正波的垂直波數。只有當krm為實數時，簡正波才會對聲場真正做出貢獻，因此簡正波階數m必須滿足因此式（18）又可以表示為：

M滿足從式（19）中可以看出當點聲源位置確定時，式中 s in(kzmz0)為定值，在深度方向上，位移勢函數只與z有關，且與z成正弦函數關系，由于kzm同樣滿足式（17），下圖給出水深100 m，聲源深度10 m，頻率100 Hz情況下的前5階簡正波模態(tài)與深度相關的位移幅值：

3.1.2 理想液體波導中簡正波階數對最佳深度的影響

為研究理想液體波導中簡正波階數對最佳深度的影響，本文計算了在圖6條件下頻率為100 Hz，海水中的簡正波階數分別為1，2，5，10時的聲傳播損失偽彩圖和不同階簡正波模態(tài)圖，結果如圖7所示。

可以看出：當波導中只有1階簡正波時，聲場的最佳深度和簡正波最大幅值均出現在50 m處；當波導中存在2階簡正波時，聲場的最佳深度和前2階簡正波疊加后幅值最大處吻合；隨著簡正波階數的增多，最佳深度從50 m處分別向海面和海底移動，波導中對稱地出現最佳深度，當海水中簡正波階數足夠多時，最佳深度穩(wěn)定在聲源深度和其對稱深度處。

下面解釋了在理想波導中為什么在聲源的對稱深度同樣會出現最佳深度的現象：

任意接收點深度和其對稱深度滿足深度之和為海水深度D，設接收點深度為zr，其對稱深度為D-zr，接收點深度為zr時的位移勢函數：

其對稱接收點深度的位移勢函數為：

由上述推導可知在理想液體波導條件下聲傳播損失的大小在深度方向上關于對稱，因此最佳深度也總是對稱存在的。

Buckingham[16]經理論推導解釋了在海水中簡正波足夠多時，最佳深度位于聲源深度和其對稱深度的原因，本文從簡正波疊加的角度來解釋上述現象：在聲源產生的所有簡正波中，各階簡正波的量級并不相同，其量級大小與該階簡正波在聲源深度處的幅度成正比，如圖7中的第5階簡正波在與聲源深度和聲源的對稱深度相同時具有最大幅值，所以第5階簡正波的量級最大，疊加后的聲場在聲源深度和聲源的對稱

深度能量最高，所以其為最佳深度。

圖7 理想波導中簡正波階數為1，2，5，10時的聲傳播損失場和簡正波模態(tài)圖Fig.7 Transmission loss field and normal modes with the number of 1,2,5,and 10 in ideal waveguide

總結本節(jié)結果得出以下結論：在聲源位置一定情況下，簡正波階數時影響最佳深度的唯一因素，當海水中只有1階簡正波時，最佳深度和簡正波幅值最大處吻合，都位于海水層中間處，隨著簡正波個數的不斷增多，最佳深度也不斷地向海面和海底移動，當簡正波個數足夠多時，最佳深度均位于聲源深度和聲源的對稱深度。

4 結 語

本文運用簡正波方法對淺海聲傳播的最佳深度進行研究，通過研究最佳頻率的影響因素和不同環(huán)境下簡正波階數對最佳深度的影響，得出以下結論：

1）聲源深度、聲源頻率和海水深度對最佳深度的影響最大，聲源頻率和海水深度通過影響海水中簡正波階數來影響最佳深度，聲源深度通過影響簡正波的量級來影響最佳深度。

2）在聲源位置一定時，最佳深度只與簡正波階數有關，當海水只有1階簡正波時，最佳深度與第1階簡正波幅值最大處位于相同位置，隨著簡正波階數增大，最佳深度從中間處分別向海面和海底移動，當簡正波個數足夠多時，最佳深度穩(wěn)定在聲源深度和聲源的對稱深度。