浅海内孤立波动态传播过程中声波模态强度起伏规律*

李沁然 孙超 谢磊†

1) (西北工业大学航海学院,西安 710072)

2) (陕西省水下信息技术重点实验室,西安 710072)

内孤立波是一种常见于浅海海域的非线性内波,具有振幅大、周期短和流速强等特点,它通过扰动水体中的温盐分布使声速剖面产生明显的距离依赖性,进而影响水下声传播特性.内波自生成后通常以1 m/s量级的速度传播,运动的内波使声传播路径上的声波模态能量在空间和时间上剧烈起伏.本文定义模态强度为模态系数模值 (模态幅度) 的平方,并用其衡量各阶模态所含声能量的大小.文中基于耦合简正波理论,推导了内波运动时声波模态强度起伏的表达式,将模态强度表征为振荡项和趋势项的线性叠加.以往的工作大多局限于单独从时域或频域研究内波运动时声波模态强度的时变规律,本文则结合短时傅里叶变换在时频平面上揭示了模态强度的起伏机理.理论推导和数值仿真均表明内孤立波使各阶声波模态之间发生能量交换,即模态耦合.内波的动态传播进一步引起模态干涉,这种干涉效应表现为模态强度中的振荡项并使模态强度随时间快速起伏.受模态剥离效应 (不同阶模态之间衰减系数的差异) 的影响,趋势项的幅度随时间不断变化,进而对模态干涉引起的振荡叠加了时变的偏置.模态强度的整体走势和振荡项中各频率分量振幅的时变特征均与模态衰减密切相关.同时,本文使用深度积分声强作为总接收声场强度的度量,研究了内波运动时模态强度起伏对接收阵位置处声能量的影响.结果表明,能量较大且起伏剧烈的模态强度将主导总接收声强的变化.

1 引言

内波是海水内部不同密度水体之间的波动,内孤立波是一种常存在于浅海海域的非线性内波.正压潮流经起伏的海底地形可生成内潮 (具有潮汐频率的内波),内潮在传播过程中受到某些非线性因素影响时其波形逐渐变陡,进一步演变为内孤立波[1,2].浅海内孤立波通常以其较大振幅剧烈扰动水体中的温盐分布,使声速剖面产生明显的距离依赖性,导致能量在声波模态之间发生转移,即模态耦合,进而引起声传播损失起伏[3,4]、声场空时相关性下降[5,6]和信号到达结构畸变[7,8].

Zhou等[9]最早研究了浅海内孤立波引起声场起伏的机制,他们通过分析黄海声学实验数据发现在某些特定频带内传播损失异常增大.后期数值仿真表明,内孤立波使各阶声波模态之间发生能量交换,位于温跃层以下的声源激发的低阶模态能量转移到衰减较大的高阶模态中导致传播损失增大.在特定频带内各阶简正波之间可发生尤为显著的能量转移,这种具有“频率共鸣”特性的模态耦合效应可引起传播损失的明显变化.Zhou等[9]仅研究了静止内孤立波引起声场起伏的机理,而实际上内孤立波自生成后一般以 1 m/s 量级的速度传播.Duda和Preisig[10]分析SWARM95实验[11]数据后发现,在声传播路径上运动的内孤立波可使内波后方声能量在数分钟内起伏达 7 dB,进一步研究表明到达内波前沿的声波模态之间的相对相位随内波位置的改变而快速变化,由此引起的模态干涉导致内波后方模态幅度和声能量随时间剧烈起伏.Rouseff等[12]根据耦合简正波理论推导了声波模态系数随内波位置变化的关系式,发现模态幅度起伏频率与模态干涉长度和内波传播速度有关,数值仿真表明模态幅度振荡曲线呈现出准周期特征.刘进忠等[13]和胡涛等[14]基于内波运动时简正波系数起伏的表达式,从模态幅度时间序列及其频谱中反演得到耦合矩阵和内波传播速度.Finette等[15]利用SWARM95实验中实测的水文数据计算出内波模型所需的环境参数,并基于Dnoidal模型[16]重构了浅海波导中空时演化的内波场及声速剖面,他们利用接收阵位置处模态强度的时变特征解释了传播损失的起伏现象,发现不同内波成分影响下的模态强度起伏程度有明显差异.Yang[17]使用 Dyson级数展开法求解随机耦合模态方程[18]获得模态系数耦合矩阵,利用该耦合矩阵推导出内波运动时接收阵位置处模态系数和深度积分声强的表达式.研究表明,模态干涉导致深度积分声强随内波运动发生起伏,从深度积分声强振荡曲线频谱的谱峰频率可识别出发生干涉的模态阶数.秦继兴等[19]根据随机耦合简正波理论推导了内孤立波移动时声强起伏的表达式,研究了内波引起的声能量变化的规律.结果表明,接收器位置处声强随时间的变化表现出明显的准周期性,不同深度阵元上声强起伏的频谱具有明显的垂直结构,该结构与主导声场起伏的简正波的模态函数形状有关.宋文华等[20]研究发现内波引起的模态耦合使主导声场的简正波成分发生变化,从而影响波导不变量的分布.Song等[21]之后进一步根据耦合简正波理论给出了存在内波时声场强度分布的表达式,发现内波引起的各阶简正波之间的能量交换导致声强中出现耦合成分.波导不变量分布中出现多个峰值,每个峰值的贡献来自于具有相近波导不变量的耦合成分.内孤立波的振幅、特征宽度和位置都会影响波导不变量分布中峰值的能量大小,并且内波位置的变化对分布中峰值的位置有显著影响.

可以看出,上述工作中多数学者仅单独在时域或频域研究了内波运动时声波模态强度 (或幅度)的变化特性.然而模态强度中通常包含多个频率分量,且受到浅海环境中尤为明显的模态衰减的影响,不同阶模态干涉形成的频率成分的时变特征也存在差异,这就导致单独的时域或频域分析均不能全面揭示内波影响下的声波模态强度起伏规律.针对以上问题,本文将模态衰减系数引入到理论推导中,给出了内孤立波运动时接收阵位置处声波模态强度的完整表达式,将模态强度表征为趋势项和振荡项的线性叠加,揭示了声衰减对模态强度起伏的调制作用.文中使用抛物方程 (parabolic equation,PE) 模型计算声场分布,之后通过模态滤波提取各阶模态强度时间序列,并使用短时傅里叶变换对模态强度进行时频分析,同时在时域和频域研究模态强度起伏特性.结果表明,内波运动时模态强度包含由模态干涉引起的多个频率分量,不同频率成分振幅的变化特征和模态强度的整体走势与波导中的模态衰减有关.同时,能量较大且起伏剧烈的模态强度将主导总接收声场强度的时变特性.

2 理论推导与讨论

2.1 模态强度表达式推导

由于内波传播时声波模态能量会随时间剧烈起伏,本文定义模态强度为模态系数模值的平方,并用其描述内波运动时声波模态中所含能量的变化情况.下面将从一个简单的波导模型出发推导模态强度起伏的表达式.本文使用的浅海波导模型如图1 所示,由单频声源发射的声信号穿越内波后到达距离声源rre处的垂直阵.假设内波出现的区域范围为r1

图1 声源与接收器的布放位置及内波运动情况Fig.1.Locations of the source and receivers and the details of internal wave motion.

式中,是第j段内第m阶模态函数,是界面rrj−1右侧第m阶模态系数,即第j段内初始位置处的模态系数,描述了初始模态系数的相位在该距离段内的变化情况.使用汉克尔函数的渐近形式,相位因子可表示为

如前述,垂直阵位置rre处的波导环境不受内波影响,则接收阵位置处简正波的水平波数和模态函数可通过在背景声速剖面下求解模态方程得到.为表达简洁,本文统一使用km,φm(z) 分别表示无内波环境中第m阶简正波的水平波数与模态函数.

在内波扰动区域内相邻距离段之间的分界面上使用声压连续条件,得到界面rrJ右侧的模态系数为

其中,矩阵Cj的元素为

式中,ρj+1(z)是第j+1段内的密度剖面;矩Hj(r)是对角阵,其元素由(2)式给出.现将(4)式代入(3)式,得到垂直阵处的声压为

其中,Rmn是耦合矩阵R的元素.波导内的模态函数满足以下正交性条件[24]:

式中,和分别是某给定波导环境条件下的模态函数和密度剖面,δmn是 Kroneckerδ函数.利用(8)式的正交性条件对垂直阵处的声压p(re,z)进行模态滤波[25],提取出任一阶(如第l阶)模态系数

其中,ρ(z) 为本文波导模型中的密度剖面,并且假设其不随时间和距离变化.

考虑到内孤立波在海洋中的传播过程,现假设内波以速度v从垂直阵附近向声源匀速运动,此时内波前沿位置rr1和后沿位置rrJ均随时间T变化,且分别表示为

式中,表示模态衰减系数.为突出接收模态系数与声源位置处模态系数之间的关系,将rr0处的模态系数展开为

式中,bn与第n阶模态激励[27]成正比:

其中,zzs是声源深度.上文已假设内孤立波仅存在于区域r1

由(14)式得到垂直阵处第l阶模态强度为

利用欧拉公式和以下复数运算法则:

(15) 式可进一步整理为

式中,D1Re[(Rlnbn)(Rlsbs)∗],D2Im[(Rlnbn)(Rlsbs)∗] .对 (17) 式继续使用三角恒等式,则第l阶接收模态强度最终化简为

式 中,θlnsarctan(D2/D1),−π/2≤θlns≤π/2 ;H[x]是Heaviside阶跃函数:

至此,得到了内波运动时声波模态强度的表达式,下面将围绕(18)式做进一步分析.

2.2分析与讨论

(18)式给出的第l阶接收模态强度可看作直流成分 DCl(t)和交流成分 ACl(t) 的线性叠加(此处的“直流成分”和“交流成分”借用了电气工程学科中的专业术语):

两种成分的表达式分别为

其中,DCln(T)和 AClns(T) 别是直流分量和交流分量.如 (21)式和(22)式所示,它们各自组成了直流成分 DCl(T)和交流成分 ACl(T) .下面将根据理论推导简要讨论内波动态传播时声波模态强度的起伏机理,以及模态强度中直流成分和交流成分的时变规律.

内波出现在声传播路径上时,其扰动区域内相邻距离段之间的模态函数不再满足(8)式给出的正交性条件,(6)式表明此时耦合矩阵R的非对角线元素非零.从(22)式余弦函数前的模值项可以看出,第n,s阶入射模态的部分能量转移进第l阶接收模态中,发生模态耦合.内波的动态传播引起入射模态之间相对相位φns(krn−krs)()的快速变化,从而使(22)式中的余弦项振荡,发生模态干涉.模态干涉分量构成了模态强度中的交流成分,也是导致内波运动过程中声波模态强度起伏的原因.

单个交流分量即模态干涉分量 AClns(T) 的振荡频率fns|(krn−krs)v/(2π)|取决于模态干涉长度L2π/(krn−krs) 和内波传播速度v.通过观察(22)式中调制余弦振荡的指数函数可以发现,接收模态衰减系数αl与入射模态衰减系数αn,αs之和的相对大小决定了对应的模态干涉分量AClns(T)振幅的变化特征.当第l阶接收模态衰减系数的2倍大于第n,s阶入射模态衰减系数之和时(2αl>αn+αs),交流分量 AClns(T) 的振幅随时间减小,反之则随时间增大.

模态强度中的直流成分 DCl(T) 在交流成分随时间的快速振荡上叠加了时变的偏置,影响了模态强度的整体趋势.与交流分量类似,直流分量DCln(T)的变化趋势也取决于接收模态衰减系数αl与入射模态衰减系数αn的相对大小.从(21)式中指数函数的自变量可以看出,当第l阶接收模态衰减系数大于第n阶入射模态衰减系数时(αl>αn),直流分量 DCln(T) 随时间减小,反之则随时间增大.

从以上讨论可以看出,不同入射模态干涉形成的交流分量振荡幅度的时变特征和振荡频率都有差异,并且所有直流分量的变化趋势将共同决定模态强度的整体走向,以上复杂性质决定了单从时域或频域均无法全面揭示模态强度的起伏规律.本文将对仿真得到的模态强度时间序列进行时频分析,上文讨论到的模态强度起伏特性均会体现在时频平面内.

3 仿真实验与分析

3.1 仿真环境参数

图2(a)给出了无内波时的背景声速剖面(与距离r、时间T无关),其表达式为

图2 波导环境声速剖面 (a) 背景声速剖面;(b) 内波扰动后的声速剖面,白色虚线标注了内波波形Fig.2.Sound speed profile in the waveguide:(a) Background sound speed profile;(b) sound speed profile perturbed by internal waves,with the waveform of internal waves indicated by white dashed lines.

式 中,z00m为海面,z115m至z230m 为跃变层深度范围,H80m 为波导深度.上下两层的声速分别为c11522m/s,c21481m/s,Δ(c1−c2)/(z1−z2)为中间跃变层声速梯度.海底为底质特性均匀的半空间,海底声速为 1600 m/s,密度为1.5 g/cm3,衰减系数为0.5 dB/λ.仿真中假设一组由 3 个双曲正割形孤子组成的内孤立波包对区域r1

式中,孤子振幅An15,12,10 m,n1,2,3 ;特征宽度Λn100,130,160 m,n1,2,3 .相邻孤子间距分别为d2−d1600 m,d3−d2500 m .受内波扰动后声速剖面的表达式为

本文以波包中首个孤子波峰的位置表示内波位置rp,图2(b)给出了位于rp26 km 的内波扰动后的声速剖面,白色虚线标注了内孤立波包波形.声场计算频率为 200Hz,声源位于跃层以下的zs37m 深度处以充分激发衰减较小的低阶模态.分别使用IFDPE[28]和KRAKEN[29]软件计算声场和模态函数.

下面首先探讨内孤立波引起的声波模态耦合效应,揭示模态强度变化特征对内波位置的高度敏感性.接下来再对内波运动时声波模态强度时间序列进行时频分析,同时在时域和频域分析模态强度的起伏现象.最后使用深度积分声强度量总接收声场强度,分析模态强度起伏对接收阵位置处声能量变化的影响.

3.2 模态耦合效应及其对内波位置的敏感性分析

从(4)式可以看出,(5)式定义的模态系数耦合矩阵表征了各阶声波模态在区域r1

式中,M8 为背景波导环境中传播模态的个数,其余计算参数分别为r00.2,r10.5,rJ2.7,re2.72km,重复以上步骤 1)— 3)直到耦合矩阵按列填充完毕.

图3(a)和图3(b)分别给出了区域r1

图3 模态系数耦合矩阵的模值 |R| (a) 无内波;(b) 有内波Fig.3.Modulus of the coupling matrix |R| :(a) Internal waves are not present in the waveguide;(b) internal waves are present in the waveguide.

2.2 节中讨论到入射模态之间的相对相位φns(krn−krs)()在内波运动过程中快速变化,从而使交流成分中的余弦项振荡,这种模态干涉效应导致内波后方的模态强度出现起伏.所以模态耦合的最终结果不仅取决于耦合矩阵,还与内波位置密切相关.在整个声传播过程中,耦合矩阵可以看作是内孤立波在当前波导环境中所固有的声学效应,该矩阵表征了各阶模态之间的能量转移关系.耦合矩阵中非对角元素的出现表明声能量可以从入射模态转移进接收模态,模态干涉则进一步影响接收模态所分配到的能量大小.从(22)式可以看出,当第n,s阶入射模态发生相长干涉时,第l阶接收模态强度出现极大值;反之当入射模态发生相消干涉时,接收模态强度出现极小值,所以模态干涉也将影响模态耦合过程的最终结果,即内波后方声能量在各阶模态中的分布情况.

图4、图5、图6分别给出了无内波、内波位于rp3km和内波位于rp2.8 km 时的内波波形(分别见图4—图6的子图 (a))以及前 5 阶模态强度随距离的变化 (分别见图4—图6的子图 (b),几何扩展效应已从结果中移除).由于图4是无内波的情况,图4(a)中内波幅度恒为η0 m.图5(b)和图6(b)中阴影覆盖的区域是内波出现的范围.图4(b)表明无内波时各阶模态绝热传播,大部分声能量集中在第1,2阶模态中,这是因为位于跃层下方的声源充分激发了低阶模态,且高阶模态传播时受到较大衰减.对比图4(b)、图5(b)和图6(b)可以看出,在内波扰动的区域内发生了模态耦合,模态强度在内波波形斜率较大(声速剖面的距离依赖性较强)的区域变化剧烈,在内波波峰和波谷(声速剖面的距离依赖性较弱)附近变化平缓,经过内波后声能量在各阶模态之间重新分配.对比图5(b)和图6(b)还可以发现,内波分别位于rp3km和rp2.8 km时的模态强度变化情况有很大差异.如内波位于rp3 km时,与图4(b)所示的绝热传播情况相比,第 3—5 阶模态经过内波后能量明显上升,并且第1,3阶和第4阶模态强度相近;内波位于rp2.8 km时,第4阶模态经过内波后其能量与无内波情况基本持平,但此时远低于第1,3阶模态强度.

图4 无内波时模态强度随距离的变化 (a) 内波波形;(b) 模态强度变化Fig.4.Modal intensity variation with range when internal waves are not present in the waveguide:(a) Waveform of internal waves;(b) modal intensity variation.

图5 内波位于 rp=3km 时模态强度随距离的变化 (a) 内波波形;(b) 模态强度变化Fig.5.Modal intensity variation with range when internal waves are at rp=3km :(a) Waveform of internal waves;(b) modal intensity variation.

图6 内波 位于 rp=2.8km 时模态强 度随距离的 变化(a) 内波波形;(b) 模态强度变化Fig.6.Modal intensity variation with range when internal waves are at rp=2.8km :(a) Waveform of internal waves;(b) modal intensity variation.

从以上对模态耦合效应的分析可知,模态强度的距离变化特征对内波位置的改变高度敏感,因此有必要深入研究内波运动时模态强度的起伏规律.

3.3 内波传播过程中模态强度起伏特性的时频分析

3.3.1 交流成分起伏分析

以下仿真中假设内孤立波包从起点rp26 km以v1m/s的速度匀速运动到终点rp1 km,对波导环境观测的快拍间隔为5 s,共产生5001个环境快拍.然后使用每组快拍下的声速分布进行声场计算,接下来对位于re30 km处垂直阵上的声压进行模态滤波得到模态强度时间序列,最后对该时间序列进行短时傅里叶变换得到其时频图.图7和图8分别给出了第1,3阶接收模态强度的时频分析结果,子图(a),(b)和(c)分别是归一化的模态强度时间序列、模态强度时频图和模态强度频谱,由第n,s阶入射模态干涉形成的频率分量对应的谱峰在图中标注为“Moden&s”.

从图7(a)和图8(a)可以看出,内波在声传播路径上运动时,接收阵处的声波模态能量剧烈起伏,模态强度振荡曲线中包含多组频率分量.(22) 式表明交流分量AClns(t) 的振荡频率为

结合(27)式和图7(c)、图8(c)中各谱峰的位置不难发现,由阶数间隔较大的第n,s阶入射模态干涉形成的交流分量AClns(t) 具有较高频率.从(22)式还可以看出,影响频谱中各谱峰幅度的因素十分复杂,包括受声源深度影响的模态激励(体现在|bn|,|bs|中)、由有损海底导致的模态衰减(体现在exp(−αn),exp(−αs) 中)、由内孤立波引起的模态耦合(体现在|Rln|,|Rls|中)和指数函数对各频率成分在时间上的调制.对比图7(a)和图8(a)可以发现,不同阶模态强度的起伏程度存在显著差异,这具体表现为第1阶模态强度起伏程度随时间增大而第3阶模态强度的振荡幅度随时间减小.3.2节中讨论到模态衰减系数随着模态阶数升高而单调增大,则对于第l1 阶接收模态,条件2αl|l=1<αn+αs对所有入射模态(第n,s阶模态)均成立,由(22)式可知此时第1阶接收模态强度中的所有交流分量的振幅均随时间增大.从图7(b)的时频平面中也可以看出,各频率分量的能量随时间增强,这最终导致第1阶模态强度的起伏程度逐渐增大.对于第l3 阶接收模态,条件2αl|l=3>αn+αs对图8(c)中标注的入射模态(即第n1,s2阶、第n2,s3阶、第n1,s3阶、第n2,s4阶和第n1,s4阶入射模态)均成立.从图8(b)的时频图中可以看出,以上入射模态干涉形成的频率分量的能量随时间减小,在时域中则体现为这些交流分量的振幅逐渐减小.注意到图8(c)中标注的这些频率分量的能量远大于其他频率成分,从而将主导模态强度的变化特征,导致第3阶模态强度的起伏程度变弱.

图7 第 1 阶接收模态强度时频分析结果 (a) 时间序列;(b) 时频图;(c) 频谱Fig.7.Time-frequency analysis results of the first received mode:(a) Time series;(b) time-frequency plot;(c) frequency spectrum.

图8 第 3 阶接收模态强度时频分析结果 (a) 时间序列;(b) 时频图;(c) 频谱Fig.8.Time-frequency analysis results of the third received mode:(a) Time series;(b) time-frequency plot;(c) frequency spectrum.

3.3.2 直流成分变化分析

对比图7(a)和图8(a)可以看出,第1阶和第3阶模态强度的整体走势也存在显著差异,第1阶模态强度曲线随时间略微向上偏移,而第3阶模态强度曲线明显向零值偏移.图9中的蓝色、红色曲线分别给出了从第1,3阶模态强度时频图的 0 Hz处提取出的值,图9反映出的直流成分时变特征与上述现象相符.如前文所述,模态强度中振荡的交流成分会被叠加上时变的直流偏置,(21)式表明模态剥离(不同阶模态之间衰减系数的差异)是这种偏置具有时变特征的原因.对于第l1 阶接收模态,条αl|l=1<αn对除第1阶以外的所有入射模态(第n2—8 阶模态)均成立,由(21)式可知这些入射模态对应的直流分量 DCln(T)|l=1,n=2—8均随时间增大.第1阶入射模态对应的直流分量DCln(T)|l=1,n=1为常数,从而对直流成分的时变特征没有贡献.以上原因导致直流成分 DCl(T)|l=1呈上升趋势,第1阶接收模态强度向上偏移.

图9 第1阶 (蓝色曲线) 和第3阶 (红色曲线) 接收模态强度时频图中0 Hz成分随时间的变化Fig.9.Variation of the 0 Hz-component in the time-frequency plane of the first (blue line) and the third (red line)modal intensity.

第l3 阶接收模态的情况则比较复杂,因为条件αl|l=3>αn仅对第n1,2 阶入射模态成立,对第n4—8阶入射模态则有αl|l=3<αn,这意味着直流分量 DCln(T)|l=3,n=1,2随时间减小,而直流分量 DCln(T)|l=3,n=4—8随时间增大,但图8(a)和图9表明第3阶模态强度中的直流成分DCl(T)|l=3整体仍表现出下降趋势.由于直流成分中的各直流分量无法像交流分量一样被分解到时频平面上,下面将通过考察第3阶模态强度中直流分量的时不变幅度来衡量各分量在直流成分中的贡献.因为(21)式表明第3阶入射模态(n3)对应的直流分量 DCln(T)|l=3,n=3为常数,其对直流成分的时变特征没有贡献,故下文将不再讨论该分量的影响.从(21)式得到单个直流分量 DCln(T)|l=3的时不变幅度为

从(28)式可以看出,直流分量的时不变幅度与模态激励(体现在|bn|中)、模态耦合(体现在|Rln|l=3中)和模态衰减(体现在exp(−2αn)中)这3种效应有关.耦合矩阵模值已在图3(b)中给出,图10(a)和 图10(b) 分别给出了各阶模态归一化的激励程度和衰减系数.从图10(a)可以看出,声源处各阶传播模态均受到充分激发,并且第1,2阶模态的激励程度大于其余模态.图3(b)表明矩阵元素|Rln|l=3,n=1,|Rln|l=3,n=2,|Rln|l=3,n=4的值大于同行中其他元素,则经过内波后第1,2阶和第4阶模态的能量更容易转移进第3阶模态.从图10(b)可以看出,模态衰减系数随模态阶数升高而单调增加,高阶模态(第 4—8 阶模态)在传播过程中比低阶模态(第1,2阶模态)受到更大衰减.以上现象说明,虽然声波在穿越过内波后会有部分能量从第1,2阶和第4阶模态耦合进第3阶模态,但受模态衰减影响,(28)式表明第n4—8 阶入射模态对应的直流分量 DCln(T)|l=3,n=4—8在直流成分中的贡献将小于第n1,2 阶入射模态对应的直流分量 DCln(T)|l=3,n=1,2.

图10 各阶传播模态归一化的激励程度和衰减系数:(a) 模态激励;(b) 模态衰减Fig.10.Normalized excitation and attenuation of each propagating mode:(a) Modal excitation;(b) modal attenuation.

图11(a)给出了第l3 阶接收模态强度中直流分量的归一化时不变幅度Bln|l=3.可以明显看出,第n1,2 阶入射模态对应的直流分量时不变幅度远大于其余分量的幅度,这与上文得到的结论一致.直流分量 DCln(T)|l=3,n=1,2将主导直流成分DCl(T)|l=3的变化趋势,即该直流成分的幅度将随时间减小.图11(b)给出了利用(21)式重构出的第3阶模态强度的直流成分 DCl(T)|l=3,可以看出基于理论公式得到的直流成分确实呈下降趋势,符合上文的推测.

图11(a) 第 l=3 阶接收模态强度中各直流分量的归一化时不变幅度 Bln|l=3;(b) 第 l=3 阶接收模态强度中归一化的直流成分Fig.11.(a) Normalized time-invariant amplitude of each DC component in the third modal intensity;(b) normalized DC part in the third modal intensity.

3.4 声波模态强度起伏对总接收声强的影响

上文同时在时域和频域分析了内波传播过程中声波模态强度的起伏特性,而与模态强度相比,声场强度在实际中受到更广泛的关注,下面将研究动态内波场中声波模态强度起伏对接收阵位置处声强的影响.

式中,p(re,z;T) 为T时刻下垂直阵接收到的声压.利用声场的简正波表达式

(28) 式可展开为

海底中的声衰减使得声波模态函数φ(z) 为复函数.本文仅考虑在远场有主要贡献的传播模态,这些简正波衰减系数较小,其模态函数的虚部相比于实部可被忽略[26],故在(31)式中使用了近似条件φ(z)φ∗(z).根据(8)式给出的正交性条件,(31)式最终化简为

(32)式表明接收阵处的深度积分声强可表示为到达该基阵的各阶声波模态强度之和.

为考察各阶模态强度对深度积分声强的贡献,图12(a)和图12(b)分别给出了接收阵位置处所有传播模态的强度和深度积分声强随时间的变化情况(几何扩展已从结果中移除),表1总结了图12中模态强度和深度积分声强起伏的均值和标准差.表1中的均值反映了各声学量的平均水平,结合图12(a)可以看出,第1,2阶模态强度远高于其余(第 3—8 阶)模态,从而在接收阵位置处的总声强中占比最大.表1中的标准差反映了各项声学量的起伏程度,第1,2阶模态强度标准差分别高达0.0144,0.0079,且明显高于其余模态强度.图12(a)也表明第1阶模态强度起伏最为剧烈,而第2阶模态强度次之.由于第1阶模态的能量在总接收声强中占比最大且起伏最为显著,因此其将主导深度积分声强的变化特征.对比图12(a)和图12(b)可以看出,深度积分声强变化与第1阶模态强度起伏呈现高度正相关,并且经计算,两者的相关系数高达0.87.从以上讨论可以看出,深度积分声强的变化由各阶模态强度控制,能量最大且起伏最为剧烈的模态将主导总接收声强的变化.

图12(a) 各阶传播模态强度随时间的变化;(b) 深度积分声强随时间的变化Fig.12.(a) Temporal variations of modal intensities of propagating modes;(b) temporal variation of the depth-integrated intensity.

表1 深度积分声强和各阶模态强度起伏的均值和标准差Table 1.Mean values and the standard deviations of the depth-integrated intensity and the modal intensities.

4 结论

本文基于耦合简正波理论,推导了内孤立波运动时声波模态强度的表达式,并结合短时傅里叶变换同时在时域和频域研究了声波模态强度的时变特性.研究表明:1)内孤立波可引起显著的模态耦合效应,声波在穿越内波时各阶模态之间发生能量交换.内孤立波的动态传播导致模态干涉,从而使内波后方模态强度随时间起伏;2)模态强度由交流成分和直流成分线性叠加组成,交流成分由模态干涉引起,表现为模态强度中随时间快速振荡的部分.直流成分对这种振荡叠加了时变的偏置,这种时变特征由模态剥离(不同阶模态之间衰减系数的差异)引起;3)模态强度中各交流分量的频率取决于模态干涉长度和内波传播速度,各频率分量的振幅随时间的变化趋势取决于对应入射模态(到达内波前沿的声波模态)和接收模态(到达接收阵的声波模态)之间衰减系数的差异.具体来说,当接收模态衰减系数的2倍大于入射模态衰减系数之和时,由入射模态干涉形成的频率分量的振幅随时间减小,反之则随时间增大;4)与模态干涉形成的频率分量类似,当接收模态衰减系数大于入射模态衰减系数时,对应的直流分量幅度随时间减小,反之则随时间增大;5)声波模态强度起伏对总接收声场强度的变化有明显影响.深度积分声强(总接收声强的度量)可表示为接收阵位置处各阶模态强度的线性叠加,能量较大且起伏剧烈的模态强度将主导深度积分声强的时变特征.