孟春霞，馬忠成

(1. 水下測控技術(shù)重點(diǎn)實(shí)驗(yàn)室，遼寧 大連 116013；2. 大連測控技術(shù)研究所，遼寧 大連 116013)

為了描述水下聲源并研究其輻射的聲場，經(jīng)常使用全向點(diǎn)輻射器模型，稱其為單極子模型。但實(shí)際上聲源在一定程度上都具有指向性，即聲源在介質(zhì)中產(chǎn)生的聲壓幅度和相位與觀測點(diǎn)的位置有關(guān)。在同一水平面內(nèi)，沿船一周測得的輻射噪聲分布曲線，明顯是不均勻的，船艏和船艉方向上要小一些，前者是因?yàn)榇w對螺旋槳噪聲的屏蔽，后者則因?yàn)槲擦鞯钠帘嗡斐蒣1]。因此，采用單極子模型難以對聲場特性準(zhǔn)確描述。近年來，研究者提出了一些輻射噪聲模型和聲場計(jì)算方法。Бреховских等[2]利用簡正波展開研究了不同波導(dǎo)和不同邊界條件下單極子輻射器的聲場。Хауга等[3]研究了平行分層波導(dǎo)中有限尺寸的指向性聲源的簡正波理論。Кулакова等[4-6]的工作中研究了含有指向性和無指向性聲源的波導(dǎo)中各種局部情況下利用簡正波展開的特點(diǎn)。被積函數(shù)在積分區(qū)域內(nèi)不僅有極點(diǎn)還有分支點(diǎn)，在近場除了考慮簡正波離散譜外，還必須考慮旁側(cè)波形式的連續(xù)譜[7-8]。Кузнецов等[9-10]對水下指向性聲源聲場特性進(jìn)行了持續(xù)的研究，對指向性聲源聲場特性的數(shù)學(xué)物理問題進(jìn)行了細(xì)致闡述。孟春霞[11]研究了在船舶的艏部、舯后部和艉部，用不同權(quán)系數(shù)的多極子來描述船舶的水下輻射噪聲場，并利用漁政船水下輻射噪聲數(shù)據(jù)對多極子模型的有效性進(jìn)行了驗(yàn)證。羅建[12]把艦船作為具有主機(jī)、輔機(jī)和螺旋槳3個(gè)輻射噪聲亮點(diǎn)的體積目標(biāo)，研究了其輻射噪聲的空間指向性的重構(gòu)方法。

本文推導(dǎo)出多極子聲源聲場勢函數(shù)，并針對能夠遠(yuǎn)距離傳輸?shù)牡皖l段，深入探討波導(dǎo)中不同階次多極子聲特性的衰減規(guī)律。

1 多極子聲場勢函數(shù)的數(shù)學(xué)表示

令有限尺寸并且具有任意形狀的單頻聲源位于均勻無限空間中。在聲源的內(nèi)部選擇任意一點(diǎn)O，這點(diǎn)是半徑為r0的球面S0的中心，使得聲源完全處于球面內(nèi)，如圖1所示。

圖1 多極子聲源模型概念的基礎(chǔ)Fig.1 Basis of the model concept for multi-pole sources

(1)

Δψ(M)+k2ψ(M)=0,M∈Ω

(2)

在邊界上具有確定的值：

(3)

式中：Δ是拉普拉斯算子；k=ω/c是波數(shù)；c是介質(zhì)的聲速；M是研究的區(qū)域Ω或它的邊界S0上的點(diǎn)。在無窮遠(yuǎn)處所求的勢函數(shù)ψ滿足索莫菲爾德條件。

在利用指向性點(diǎn)狀多極子輻射器模型時(shí)，省略時(shí)間因子e-iωt，無界均勻空間中聲場勢函數(shù)為：

(4)

式(4)是與有限尺寸的實(shí)際聲源等效的指向性點(diǎn)狀輻射器的聲場勢函數(shù)。在包含原始聲源的任意球面以外，這種點(diǎn)狀輻射器的聲場與有限尺寸的原始聲源的聲場處處吻合。模型輻射器的方向特性完全由系數(shù)Cnm確定，Cnm為復(fù)數(shù)，-n≤m≤n，該參數(shù)的物理意義是構(gòu)成展開式的初級球面多極子的力矩。對于式(4)來說，利用特殊函數(shù)的積分表達(dá)式，可以得到多極子聲源聲場的積分表達(dá)式為平面波之和的展開形式：

(5)

當(dāng)L=0時(shí)，這個(gè)表達(dá)式實(shí)際上是無指向性單極子輻射器的勢函數(shù)；當(dāng)L=2時(shí)，相當(dāng)于研究到四極子，其中單極子有1種形式，偶極子有3種形式，四極子有5種形式，這時(shí)，指向性點(diǎn)狀輻射器是包含18個(gè)未知參數(shù)的模型。在聲學(xué)正問題中，在研究各個(gè)多極子非相干分量的特性時(shí)，為簡便起見假設(shè)它們的多極子矩等于1，即Cnm=1。求解聲學(xué)逆問題，實(shí)際上就是根據(jù)實(shí)測的輻射聲壓值，對等效源模型中的各個(gè)多極子系數(shù)Cnm進(jìn)行求解。

1.1 多極子勢函數(shù)的簡正波近似

考慮海底和海面的聲學(xué)邊界條件，經(jīng)過推導(dǎo)，可得到Pekeris波導(dǎo)中多極子聲源聲場勢函數(shù)的簡正波近似形式的表達(dá)式：

(6)

(7)

聲場勢函數(shù)可以表示成非相干分量和相干分量之和的形式，非相干分量確定聲場的平均聲級，相干分量確定聲場在平均聲級附近的振蕩。多極子聲場勢函數(shù)幅度平方非相干分量的表達(dá)式為：

(8)

(xl/kh))2exp(-δrl2)

(9)

式(9)是單個(gè)選取的多極子聲場勢函數(shù)的幅度平方的非相干分量。

近場δr<<1時(shí)，所有號數(shù)為l的簡正波，exp(-δrl2)≈1，在中等距離exp(-δrl2)不為1；在遠(yuǎn)距離δr>1時(shí)，聲場中的指數(shù)項(xiàng)迅速衰減。利用不同的近似條件可以得到不同區(qū)域聲場的近似表達(dá)式，針對實(shí)際需求，能夠方便地對不同區(qū)域的聲傳播規(guī)律進(jìn)行估計(jì)。

1.2 多極子勢函數(shù)的圍線積分表示

圍線積分方法直接進(jìn)行數(shù)值積分，能給出聲場的精確解。利用圍線積分的方法將式(1)表示為：

(10)

式中：ψH(R,θ,φ)為無界空間中的聲場勢函數(shù)，由式(1)確定；ψF(R,θ,φ)為波導(dǎo)中來自絕對軟上邊界的聲波對總聲場的貢獻(xiàn)。

(11)

2 聲場近似解適用條件仿真

淺海波導(dǎo)通?？梢钥醋魇撬骄鶆虿▽?dǎo)，能夠用Pekeris波導(dǎo)來描述。仿真計(jì)算條件如下：海水深度60 m，水中聲速1 500 m/s，水體下方為液態(tài)半無限空間，海底聲速1 800 m/s，海底密度2 700 kg/m3,海底無量綱的吸收系數(shù)為0.1。圖2為對聲源位于海面和接收器位于海底的情況，以及聲源和接收器都位于海水中部的情況進(jìn)行仿真。

圖2 單極子勢函數(shù)非相干分量及其近場近似解Fig.2 Regular components of monopole potential function and its approximate solutions at near field

波導(dǎo)的聲學(xué)參數(shù)不變，設(shè)聲源深度5 m、接收器深度55 m，圖3給出單極子聲源頻率60 Hz和100 Hz時(shí)勢函數(shù)非相干分量與中等距離近似解及遠(yuǎn)場近似解的對比。

圖3 非相干分量精確解與近似解對比Fig.3 Exact solutions of regular components of monopole potential function and its approximate solutions

從圖3可以看出：對于水深60 m的情況來說，勢函數(shù)遠(yuǎn)場近似解適用距離大于5 km，而且頻率越高，適用距離越大?？梢娭械染嚯x近似解的適用條件很寬泛。

3 多極子衰減規(guī)律仿真

3.1 中遠(yuǎn)距離聲場勢函數(shù)衰減規(guī)律仿真

在中遠(yuǎn)距離上，與簡正波相比旁側(cè)波對聲場的貢獻(xiàn)可以忽略，聲場表現(xiàn)為簡正波疊加的形式。在式(1)中，n=1,m=-1為y偶極子；n=2，m=-2為有2個(gè)水平軸的yy四極子；n=1，m=1為有1個(gè)水平軸1個(gè)垂直軸的yz四極子，n=1，m=0為z偶極子；n=2，m=0為有2個(gè)垂直軸的zz四極子。其他類型的多極子在聲信號觀測角度上勢函數(shù)為0。

在研究聲學(xué)逆問題時(shí)，利用聲場勢函數(shù)與聲壓之間的關(guān)系：

(12)

并進(jìn)一步與實(shí)測的聲壓值進(jìn)行比對。

波導(dǎo)的聲學(xué)參數(shù)不變，圖4(a)為聲源深度5 m，接收器深度6 m頻率100 Hz不同類型多極子聲場勢函數(shù)非相干分量衰減曲線。圖4(b)為聲源深度5 m，接收器深度28 m，頻率100 Hz時(shí)不同類型多極子勢函數(shù)非相干分量衰減曲線。為了便于比較，將圖中曲線沿著垂直軸進(jìn)行了平移。

圖4 頻率100 Hz時(shí)多極子勢函數(shù)衰減曲線Fig.4 Attenuation curves of multi-poles source potential function f=100 Hz

從圖4可以看出：在中遠(yuǎn)距離上，根據(jù)聲場勢函數(shù)衰減規(guī)律，多極子可以分為2組，第1組包括：單極子、水平偶極子、具有2個(gè)水平軸的四極子、具有2個(gè)垂直軸的四極子;第2組包括：垂直偶極子與具有一個(gè)水平軸一個(gè)垂直軸的四極子。這樣，在遠(yuǎn)場分析實(shí)際聲源的指向性時(shí)，可以忽略垂直面上的指向性，并且遠(yuǎn)場條件下，實(shí)際水下聲場的多極子模型可以簡化。

3.2 近場衰減規(guī)律仿真

在聲源特性精細(xì)研究和觀測的近場，接收器與聲源之間的距離在300～400 m到1～2 km。在這個(gè)距離上旁側(cè)波對聲場的作用不能忽視，此時(shí)簡正波形式的聲場近似解無法給出聲場在空間上的細(xì)致結(jié)構(gòu)，因此采用圍線積分方法研究多極子的聲場分布。聲壓接收方向?yàn)檠刂鴛軸方向，此時(shí)y軸水平偶極子的聲場為0，因此只研究x軸水平偶極子的聲場衰減規(guī)律。

聲源深度36 m，接收器深度98 m，波導(dǎo)厚度為100 m。海底聲學(xué)參數(shù)與前述相同。圖5分別給出了聲源頻率100 Hz和50 Hz時(shí)單極子與水平偶極子和垂直偶極子聲場特性的對比結(jié)果。

在圖5中，聲信號觀測方向?yàn)樗脚紭O子極大值方向，因此聲場看不出方位的變化。從圖5可以看出：相同波導(dǎo)條件下，在近距離上，水平偶極子與單極子的衰減規(guī)律相同，垂直偶極子比單極子和水平偶極子都衰減的快，傳播損失較大。對于同一種多極子來說，頻率越高，聲信號衰減速度越快。

圖5 單極子與及垂直偶極子聲場特性對比Fig.5 Comparison of sound field characteristics between monopole and horizontal dipole and vertical dipole

水平偶極子聲源是最常見的多極子聲源，設(shè)波導(dǎo)的聲學(xué)參數(shù)不變，水平偶極子聲源沿著直線運(yùn)動(dòng)，研究聲源與接收器之間正橫距離分別為50、100和200 m時(shí)，不同頻率水平偶極子聲源聲場衰減規(guī)律，仿真結(jié)果如圖6所示。

從圖6可以看出：水平偶極子的聲場具有指向性，但是在海洋信道濾波器的作用下水平偶極子聲場的8字型指向性發(fā)生了畸變。在水平距離較小的情況下，接收器與聲源之間的正橫距離對聲場的指向性非常明顯；隨著水平距離的增大，接收器與聲源之間正橫距離對聲場的影響減弱。相同頻率不同接收深度的聲場也存在差異，這是由于在波導(dǎo)中到達(dá)不同接收位置的聲線能量不同。

圖6 水平偶極子不同正橫距離時(shí)聲場對比Fig.6 Comparison of sound field characteristics between cross section for horizontal dipole

4 結(jié)論

1)在不同的傳播距離上，可以采用不同的近似公式來描述多極子勢函數(shù)非相干分量的衰減規(guī)律。

2)在中遠(yuǎn)距離上，按照勢函數(shù)衰減規(guī)律分類，多極子可以分為2類：單極子、水平偶極子、具有2個(gè)水平軸的四極子和具有2個(gè)垂直軸的四極子；垂直偶極子及具有1個(gè)水平軸和1個(gè)垂直軸的四極子。

3)在近場水平方向上，單極子與水平偶極子傳播規(guī)律相同，垂直偶極子比單極子衰減速度快。海洋波導(dǎo)中隨著傳播距離的增大，多極子的指向性減弱。

4)本文的研究工作為描述水下聲源的聲場提供了精確的理論模型，利用該模型，能夠有效地預(yù)報(bào)海洋波導(dǎo)中低頻水下聲源輻射聲信號在不同區(qū)域的傳輸規(guī)律，具有重要的實(shí)用價(jià)值。