非線性內(nèi)波引起的淺海混響強度異常震蕩

高博，宋文華，王寧

(中國海洋大學 信息科學與工程學院，山東 青島 266100)

主動聲吶工作在淺海中時，海底散射強度大，散射機理復雜，與目標回波極為類似的海底混響中的雜波干擾更是導致虛警概率升高，限制主動聲吶探測性能的重要干擾。明確淺海海底雜波的形成機理是區(qū)分真假目標，抑制雜波干擾的重要理論基礎。海底雜波已經(jīng)成為淺?；祉懷芯康臒狳c問題，許多學者探討了不同的淺海環(huán)境中雜波干擾的形成機制。淺海雜波散射源有以下幾類：海底界面以及海底淺層底質(zhì)中的離散的強散射體[1](如海底山、海底沉船)、非連續(xù)的海底底質(zhì)結(jié)構(gòu)[2]、水體中魚群[3]、和特殊的海洋動力學現(xiàn)象[4](如非線性孤子內(nèi)波)等。文獻[5]通過仿真計算首次給出了孤子內(nèi)波會導致類似目標回波的強雜波干擾。并從射線理論角度，分析了聲波入射到孤子內(nèi)波后產(chǎn)生強折射，從而以更大的入射角與海底發(fā)生散射，散射強度增大，因此孤子內(nèi)波可能會產(chǎn)生10 dB左右的雜波干擾。文獻[4]發(fā)現(xiàn)雜波干擾過后的平均混響強度會比產(chǎn)生雜波前的平均混響強度略高，而其機理并未給出。另外，淺海信道中可能產(chǎn)生雜波的機制繁雜，如何確定某個雜波是水體起伏引起的也是亟需解決的問題。淺海信道具有典型的時變空變特性，移動的非線性孤子內(nèi)波引起的水體起伏對混響強度的影響是本文的重點關注內(nèi)容。本文為了進一步揭示內(nèi)波引起的水體起伏對淺?；祉憦姸鹊挠绊懀苑蔷€性孤子內(nèi)波為例，將孤子內(nèi)波導致的聲模態(tài)間的能量耦合嵌入到簡正波混響理論，提出了一種考慮水體耦合效應的耦合簡正波混響模型。

1 非線性內(nèi)波引起的混響與雜波理論

本文選取淺海Pekeris波導，平均海深為70 m，上層海水的平均聲速為1 500 m/s，液態(tài)半無限海底聲速為1 700 m/s。水體和液態(tài)海底密度比為1.6，海底聲吸收系數(shù)為0.4 dB/λ，海水中的聲速梯度是黃海夏季常見的典型孤子內(nèi)波環(huán)境，在15～25 m處有強的溫躍層，具體的海水中聲速分布如圖1所示。聲源接收器的配置形式為收發(fā)合置，布放深度為35 m，仿真所用信號為中心頻率3 kHz的窄帶信號。

圖1 內(nèi)波結(jié)構(gòu)和聲速梯度Fig.1 The soliton wave structure and the sound speed profile

為了仿真需要，假設孤子內(nèi)波的波包沿著聲傳播路徑移動，且移動中孤子波的波形保持不變。孤子內(nèi)波的描述采用經(jīng)典的KdV方程解的形式：

(1)

式中：rx為內(nèi)波的中心位置；Δ為內(nèi)波的寬度。內(nèi)波引起的強雜波主要是受聲波遇到內(nèi)波的第1個波包影響[4]。雖KdV方程解并非最優(yōu)，但本文的主要目的是闡述在聲波的傳播路徑上遇到孤子波包后，引起的聲場各模態(tài)簡的能量耦合現(xiàn)象。采用KdV方程解的形式刻畫內(nèi)波結(jié)構(gòu)能夠滿足本文的需求。

1.1 引入內(nèi)波耦合矩陣的簡正波混響修正模型

對于淺海收發(fā)合置的混響，包含傳播和海底散射2個過程，為了刻畫聲波在傳播過程中遇到強孤子內(nèi)波導致的聲場各模態(tài)間的能量耦合過程，本文在經(jīng)典簡正波混響模型基礎上，借助耦合簡正波理論刻畫水體不均勻?qū)е碌哪B(tài)耦合，表征孤子內(nèi)波引起的簡正波間的能量交換，揭示內(nèi)波引起的強雜波的物理機制。

淺?；祉懙睦碚撃Ｐ捅姸郲6-11]，利用經(jīng)典的簡正波混響理論，海底混響的聲壓ps(r,zr)為：

(2)

(3)

(3)

當聲波的前向傳播路徑中遭遇孤子內(nèi)波引起的水文躍變，在入射到海底發(fā)生海底散射前會產(chǎn)生強烈的模態(tài)間耦合效應。假設內(nèi)波波包存在的水平距離范圍是r1

A(r)=E(r,r2)Λ(r2,r1)E(r1,0)A(0)

(4)

式中：Λ(r2,r1)表征在距離r1

(5)

依據(jù)式(5)得到混響強度仿真結(jié)果如圖2所示。為了進一步驗證基于水體耦合矩陣的簡正波混響理論的準確性，探討有無內(nèi)波以及內(nèi)波移動后引起的混響雜波和混響強度的畸變規(guī)律。本文分別仿真了無內(nèi)波情況、有內(nèi)波且內(nèi)波移動到距聲源不同水平位置的情況。結(jié)果如圖3所示。通過對比得出，當沒有孤子內(nèi)波引起的水文躍變時，海底混響強度正常衰減，不會產(chǎn)生海底雜波；而當孤子內(nèi)波出現(xiàn)在不同位置時，與內(nèi)波位置相對應的時間上會產(chǎn)生雜波。需注意：1)在仿真時假設內(nèi)波在移動過程中保持波形不變，即內(nèi)波導致的模態(tài)間的耦合能量交換的強度不變，因此產(chǎn)生海底雜波的聲強也是一致的，出現(xiàn)的時間隨著內(nèi)波位置發(fā)生了相對應的變化；2)海底散射后的混響聲場與前向傳播聲場相比為小量，故只需考慮前向傳播過程中水體起伏引起的能量耦合過程，因此式(5)中只在前向傳播的n階模態(tài)水平因子中考慮了耦合矩陣影響，以此提高算法的運算效率。

圖2 距聲源3 km處的孤子內(nèi)波引起的海底雜波Fig.2 Simulated RL when the distance between a source and the NIW is 3 km

圖3 無孤子內(nèi)波和有孤子內(nèi)波其位于不同位置時的非相干混響強度Fig.3 Simulated incoherent RL when there is no soliton wave and when the soliton wave is existing at different position

1.2 非線性內(nèi)波對前向傳播聲場的耦合過程

在聲波傳播路徑上的非線性孤子內(nèi)波會導致各模態(tài)間的能量耦合，耦合的強度由內(nèi)波的垂向幅度η和水平寬度決定。為了表征該耦合過程，將孤子內(nèi)波結(jié)構(gòu)分為3段，如圖4中的r1、rmid、r2所示。其中r1為聲波剛剛進入內(nèi)波，rmid為聲波進入到內(nèi)波的中心位置，r2為聲波基本經(jīng)過了一個完整的內(nèi)波波包。圖5則分別給出了與r1、rmid、r2位置處相對應的耦合矩陣，各仿真結(jié)果中只展示了前100階模態(tài)間的能量轉(zhuǎn)換?？梢园l(fā)現(xiàn)：1)聲波剛進入內(nèi)波的r1位置時，水文環(huán)境開始發(fā)生變化，相應的圖5(a)中可以看出能量開始向高號簡正波轉(zhuǎn)移，高號簡正波能量出現(xiàn)增大的現(xiàn)象；2)當前向傳播的聲場進入到rmid位置時，即孤子內(nèi)波波包的中心位置，此處的水文躍變最為劇烈，大量的能量轉(zhuǎn)移到高號簡正波(如圖6所示)；3)等傳播到r2位置，即將完全經(jīng)過內(nèi)波波包時，除了部分高號簡正波的能量消逝掉，和入射內(nèi)波前相比，仍出現(xiàn)了大量的高階模態(tài)(如圖6所示)。在前向傳播路徑上有孤子內(nèi)波出現(xiàn)時，水體變化引起強烈的耦合效應，大量能量轉(zhuǎn)移到高號模態(tài)，海底混響的入射角增大，海底散射強度增加，從而產(chǎn)生類目標的強雜波。仿真結(jié)果表明，當內(nèi)波中心位置在距離聲源3 km處時，考慮了水體耦合效應后的相干混響強度和非相干混響強度的衰減規(guī)律均在混響時間為4 s處(對應于空間上的3 km處的內(nèi)波位置)產(chǎn)生了類似目標回波的強雜波。

圖4 孤子內(nèi)波的結(jié)構(gòu)示意Fig.4 Geometry of the nonlinear internal wave structure

用更精準的非線性方程刻畫出更精細的內(nèi)波結(jié)構(gòu)，有助于得到更為準確的耦合矩陣和海底混響雜波結(jié)構(gòu)，但本文重點關注的是利用模態(tài)耦合矩陣分析混響雜波的物理機制和過程。圖5的結(jié)果表明，本文提出的基于內(nèi)波耦合矩陣的簡正波混響模型可以詳細的表征孤子內(nèi)波引起的海底混響雜波的耦合機理和物理過程。

2 移動內(nèi)波引起的混響強度起伏

在圖3的算例中，由不同位置的孤子內(nèi)波引起的不同時刻的混響雜波到達接收點之前，海底混響強度是一致的。而類目標的雜波出現(xiàn)后的平均混響強度是明顯不同的。與沒有內(nèi)波的情況相比，不同位置的內(nèi)波會引起雜波后的混響強度增加或減小。以第8 s的混響時間為例，即對應于距離聲源6 km位置處的海底散射產(chǎn)生的混響。此時聲波已經(jīng)過內(nèi)波的波包結(jié)構(gòu)。圖3中所示內(nèi)波距離聲源3.2 km時，海底雜波結(jié)構(gòu)后面的混響強度要高于沒有內(nèi)波時的情況。當內(nèi)波距離聲源3 km和3.5 km時，海底雜波結(jié)構(gòu)到達接收點后的混響強度均不同程度的弱于沒有內(nèi)波時的混響強度。

圖5 利用耦合矩陣表示處各模態(tài)間的能量耦合Fig.5 The energy redistributed at position by couple matrix

圖6 移動內(nèi)波引起的混響強度震蕩Fig.6 The fluctuation of RL induced by moving internal wave

為了進一步解釋海底雜波結(jié)構(gòu)后面的混響強度異?，F(xiàn)象，進行了如下環(huán)境的仿真：假設聲源在35 m深度處，每隔10 s發(fā)射一次中心頻率為3 kHz 脈寬為0.1 s的CW窄帶脈沖。設2 000 s內(nèi)發(fā)射200個信號，同時有圖1所示的孤子內(nèi)波在2 000 s內(nèi)以0.5 m/s的速度遠離聲源，即孤子內(nèi)波開始出現(xiàn)的位置由距聲源3 km移動至4 km。運動過程中，設孤子波的波包結(jié)構(gòu)不變。則依據(jù)前文理論，類目標的海底雜波結(jié)構(gòu)出現(xiàn)的初始時刻也從第4 s連續(xù)變化至5.33 s，如圖7所示。

圖7 隨著移動內(nèi)波第8 s時刻的混響強度震蕩Fig.7 The fluctuation of RL at the moment of the 8 s

將圖7中第8 s固定為觀測時間點，有移動內(nèi)波時該觀測時間點的非相干混響強度定義為RLs，無孤子內(nèi)波時，定義該觀測時間點的非相干混響強度為RL0，在內(nèi)波移動的2 000 s內(nèi)該觀測時間點的混響強度變化量可以表述為：ΔRL=RLs-RL0。圖8為ΔRL的變化規(guī)律。

結(jié)合圖6和圖7可以發(fā)現(xiàn)，移動內(nèi)波引起混響雜波結(jié)構(gòu)到達之后的某固定時間點的混響強度起伏ΔPL呈現(xiàn)出一種不規(guī)則的準周期震蕩結(jié)構(gòu)。為了探究其物理機制，對該結(jié)構(gòu)進行了傅里葉譜分析。圖8給出了其歸一化的頻譜結(jié)構(gòu)，當前向傳播的聲波經(jīng)過移動內(nèi)波后，固定觀測時間點的混響強度震蕩ΔPL的主要頻率分量與某對簡正波間的干涉跨度有關，其關系式為fmn=vkmn/2π。其中，kmn=km-kn為第m階和第n階簡正波的水平波數(shù)差，v表示內(nèi)波的移動速度。

圖8 混響強度隨時間起伏的頻域譜Fig.8 Spectrum of the RL fluctuation shown

本文的仿真算例中，孤子波以0.5 m/s的速度運動，該移動內(nèi)波導致了海底雜波到達后的混響強度產(chǎn)生不規(guī)則的準周期震蕩，其震蕩結(jié)構(gòu)的前4階頻率分量，分別由波數(shù)差為2、3、4、9的簡正波之間的干涉周期決定，如圖8所示。內(nèi)波移動造成前向傳播的聲場與內(nèi)波發(fā)生耦合的位置隨時間變化，由此導致了不同模態(tài)之間相位差發(fā)生變化。孤子內(nèi)波作為一個移動的散射體，引起的模態(tài)間的干涉也存在相應的起伏震蕩，而起伏震蕩的周期則和某對模態(tài)之間的干涉跨度相關，對應的時間周期等于干涉跨度除以內(nèi)波速度。

移動孤子內(nèi)波導致的混響強度在某個確定時刻產(chǎn)生不規(guī)則震蕩現(xiàn)象，即觀測的混響起伏是雜波結(jié)構(gòu)到達后同一個混響時刻的非相干混響強度的起伏，這個物理過程相當于固定的水平位置觀測內(nèi)波移動帶來的相位起伏。類似的現(xiàn)象在孤子內(nèi)波環(huán)境下的聲傳播起伏研究中已有論述，相關理論結(jié)果和物理機制得到了SWARM實驗數(shù)據(jù)的驗證[12]。而內(nèi)波的移動可以導致混響雜波和雜波后混響強度的異常震蕩現(xiàn)象尚未有文獻報道。而文獻[4]中所提的雜波結(jié)構(gòu)后混響強度增大，只是本文所述震蕩現(xiàn)象中某個時刻的特例。雜波結(jié)構(gòu)后緊跟的混響強度增加還是減小，取決于前向傳播的聲波與孤子內(nèi)波發(fā)生耦合的位置。利用本文中所述的耦合簡正波混響理論可以解釋該現(xiàn)象的機理。

3 結(jié)論

1)淺海環(huán)境中非線性內(nèi)波引起的模態(tài)間的能量耦合效應，是淺?；祉戨s波的一種形成機制。當聲波進入內(nèi)波至內(nèi)波波包的中心位置處，能量由低號模態(tài)向高號模態(tài)轉(zhuǎn)移，產(chǎn)生大量的高階模態(tài)，海底散射強度增大，海底雜波強度也達到峰值。而聲波完全穿過內(nèi)波，高號簡正波的能量又向低號轉(zhuǎn)移，只有部分高號簡正波得以存活，因此雜波結(jié)構(gòu)也回落至正常的海底混響強度。

2)移動非線性內(nèi)波會使得海底雜波結(jié)構(gòu)到達之后的混響強度發(fā)生不規(guī)則的準周期震蕩現(xiàn)象。震蕩頻率與內(nèi)波的速度，內(nèi)波引起某對模態(tài)間的干涉跨度有關，具體的關系是為fnm=vkmn/2π。這種雜波后的混響強度震蕩現(xiàn)象有望用來區(qū)分移動目標(孤子內(nèi)波、AUV等)和固定目標(沉船、海底山、海底沉積層結(jié)構(gòu)等)產(chǎn)生的雜波。

3)由非線性內(nèi)波對聲場的耦合效應可知，孤子內(nèi)波引起的海底混響雜波結(jié)構(gòu)、到達時刻分別和內(nèi)波的精細結(jié)構(gòu)、位置有關。因此，移動孤子內(nèi)波導致的混響雜波和混響起伏可以成為獲取內(nèi)波速度、距離的信息載體。