波导中多极子声源低频声场特性衰减规律

2020-12-15孟春霞马忠成

哈尔滨工程大学学报 2020年10期

孟春霞，马忠成

(1. 水下测控技术重点实验室，辽宁大连 116013；2. 大连测控技术研究所，辽宁大连 116013)

为了描述水下声源并研究其辐射的声场，经常使用全向点辐射器模型，称其为单极子模型。但实际上声源在一定程度上都具有指向性，即声源在介质中产生的声压幅度和相位与观测点的位置有关。在同一水平面内，沿船一周测得的辐射噪声分布曲线，明显是不均匀的，船艏和船艉方向上要小一些，前者是因为船体对螺旋桨噪声的屏蔽，后者则因为尾流的屏蔽所造成[1]。因此，采用单极子模型难以对声场特性准确描述。近年来，研究者提出了一些辐射噪声模型和声场计算方法。Бреховских等[2]利用简正波展开研究了不同波导和不同边界条件下单极子辐射器的声场。Хауга等[3]研究了平行分层波导中有限尺寸的指向性声源的简正波理论。Кулакова等[4-6]的工作中研究了含有指向性和无指向性声源的波导中各种局部情况下利用简正波展开的特点。被积函数在积分区域内不仅有极点还有分支点，在近场除了考虑简正波离散谱外，还必须考虑旁侧波形式的连续谱[7-8]。Кузнецов等[9-10]对水下指向性声源声场特性进行了持续的研究，对指向性声源声场特性的数学物理问题进行了细致阐述。孟春霞[11]研究了在船舶的艏部、舯后部和艉部，用不同权系数的多极子来描述船舶的水下辐射噪声场，并利用渔政船水下辐射噪声数据对多极子模型的有效性进行了验证。罗建[12]把舰船作为具有主机、辅机和螺旋桨3个辐射噪声亮点的体积目标，研究了其辐射噪声的空间指向性的重构方法。

本文推导出多极子声源声场势函数，并针对能够远距离传输的低频段，深入探讨波导中不同阶次多极子声特性的衰减规律。

1 多极子声场势函数的数学表示

令有限尺寸并且具有任意形状的单频声源位于均匀无限空间中。在声源的内部选择任意一点O，这点是半径为r0的球面S0的中心，使得声源完全处于球面内，如图1所示。

图1 多极子声源模型概念的基础Fig.1 Basis of the model concept for multi-pole sources

(1)

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

(2)

在边界上具有确定的值：

(3)

式中：Δ是拉普拉斯算子；k=ω/c是波数；c是介质的声速；M是研究的区域Ω或它的边界S0上的点。在无穷远处所求的势函数ψ满足索莫菲尔德条件。

在利用指向性点状多极子辐射器模型时，省略时间因子e-iωt，无界均匀空间中声场势函数为：

(4)

式(4)是与有限尺寸的实际声源等效的指向性点状辐射器的声场势函数。在包含原始声源的任意球面以外，这种点状辐射器的声场与有限尺寸的原始声源的声场处处吻合。模型辐射器的方向特性完全由系数Cnm确定，Cnm为复数，-n≤m≤n，该参数的物理意义是构成展开式的初级球面多极子的力矩。对于式(4)来说，利用特殊函数的积分表达式，可以得到多极子声源声场的积分表达式为平面波之和的展开形式：

(5)

当L=0时，这个表达式实际上是无指向性单极子辐射器的势函数；当L=2时，相当于研究到四极子，其中单极子有1种形式，偶极子有3种形式，四极子有5种形式，这时，指向性点状辐射器是包含18个未知参数的模型。在声学正问题中，在研究各个多极子非相干分量的特性时，为简便起见假设它们的多极子矩等于1，即Cnm=1。求解声学逆问题，实际上就是根据实测的辐射声压值，对等效源模型中的各个多极子系数Cnm进行求解。

1.1 多极子势函数的简正波近似

考虑海底和海面的声学边界条件，经过推导，可得到Pekeris波导中多极子声源声场势函数的简正波近似形式的表达式：

(6)

(7)

声场势函数可以表示成非相干分量和相干分量之和的形式，非相干分量确定声场的平均声级，相干分量确定声场在平均声级附近的振荡。多极子声场势函数幅度平方非相干分量的表达式为：

(8)

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

(9)

式(9)是单个选取的多极子声场势函数的幅度平方的非相干分量。

近场δr<<1时，所有号数为l的简正波，exp(-δrl2)≈1，在中等距离exp(-δrl2)不为1；在远距离δr>1时，声场中的指数项迅速衰减。利用不同的近似条件可以得到不同区域声场的近似表达式，针对实际需求，能够方便地对不同区域的声传播规律进行估计。

1.2 多极子势函数的围线积分表示

围线积分方法直接进行数值积分，能给出声场的精确解。利用围线积分的方法将式(1)表示为：

(10)

式中：ψH(R,θ,φ)为无界空间中的声场势函数，由式(1)确定；ψF(R,θ,φ)为波导中来自绝对软上边界的声波对总声场的贡献。

(11)

2 声场近似解适用条件仿真

浅海波导通常可以看作是水平均匀波导，能够用Pekeris波导来描述。仿真计算条件如下：海水深度60 m，水中声速1 500 m/s，水体下方为液态半无限空间，海底声速1 800 m/s，海底密度2 700 kg/m3,海底无量纲的吸收系数为0.1。图2为对声源位于海面和接收器位于海底的情况，以及声源和接收器都位于海水中部的情况进行仿真。

图2 单极子势函数非相干分量及其近场近似解Fig.2 Regular components of monopole potential function and its approximate solutions at near field

波导的声学参数不变，设声源深度5 m、接收器深度55 m，图3给出单极子声源频率60 Hz和100 Hz时势函数非相干分量与中等距离近似解及远场近似解的对比。

图3 非相干分量精确解与近似解对比Fig.3 Exact solutions of regular components of monopole potential function and its approximate solutions

从图3可以看出：对于水深60 m的情况来说，势函数远场近似解适用距离大于5 km，而且频率越高，适用距离越大。可见中等距离近似解的适用条件很宽泛。

3 多极子衰减规律仿真

3.1 中远距离声场势函数衰减规律仿真

在中远距离上，与简正波相比旁侧波对声场的贡献可以忽略，声场表现为简正波叠加的形式。在式(1)中，n=1,m=-1为y偶极子；n=2，m=-2为有2个水平轴的yy四极子；n=1，m=1为有1个水平轴1个垂直轴的yz四极子，n=1，m=0为z偶极子；n=2，m=0为有2个垂直轴的zz四极子。其他类型的多极子在声信号观测角度上势函数为0。

在研究声学逆问题时，利用声场势函数与声压之间的关系：

(12)

并进一步与实测的声压值进行比对。

波导的声学参数不变，图4(a)为声源深度5 m，接收器深度6 m频率100 Hz不同类型多极子声场势函数非相干分量衰减曲线。图4(b)为声源深度5 m，接收器深度28 m，频率100 Hz时不同类型多极子势函数非相干分量衰减曲线。为了便于比较，将图中曲线沿着垂直轴进行了平移。

图4 频率100 Hz时多极子势函数衰减曲线Fig.4 Attenuation curves of multi-poles source potential function f=100 Hz

从图4可以看出：在中远距离上，根据声场势函数衰减规律，多极子可以分为2组，第1组包括：单极子、水平偶极子、具有2个水平轴的四极子、具有2个垂直轴的四极子;第2组包括：垂直偶极子与具有一个水平轴一个垂直轴的四极子。这样，在远场分析实际声源的指向性时，可以忽略垂直面上的指向性，并且远场条件下，实际水下声场的多极子模型可以简化。

3.2 近场衰减规律仿真

在声源特性精细研究和观测的近场，接收器与声源之间的距离在300～400 m到1～2 km。在这个距离上旁侧波对声场的作用不能忽视，此时简正波形式的声场近似解无法给出声场在空间上的细致结构，因此采用围线积分方法研究多极子的声场分布。声压接收方向为沿着x轴方向，此时y轴水平偶极子的声场为0，因此只研究x轴水平偶极子的声场衰减规律。

声源深度36 m，接收器深度98 m，波导厚度为100 m。海底声学参数与前述相同。图5分别给出了声源频率100 Hz和50 Hz时单极子与水平偶极子和垂直偶极子声场特性的对比结果。

在图5中，声信号观测方向为水平偶极子极大值方向，因此声场看不出方位的变化。从图5可以看出：相同波导条件下，在近距离上，水平偶极子与单极子的衰减规律相同，垂直偶极子比单极子和水平偶极子都衰减的快，传播损失较大。对于同一种多极子来说，频率越高，声信号衰减速度越快。

图5 单极子与及垂直偶极子声场特性对比Fig.5 Comparison of sound field characteristics between monopole and horizontal dipole and vertical dipole

水平偶极子声源是最常见的多极子声源，设波导的声学参数不变，水平偶极子声源沿着直线运动，研究声源与接收器之间正横距离分别为50、100和200 m时，不同频率水平偶极子声源声场衰减规律，仿真结果如图6所示。

从图6可以看出：水平偶极子的声场具有指向性，但是在海洋信道滤波器的作用下水平偶极子声场的8字型指向性发生了畸变。在水平距离较小的情况下，接收器与声源之间的正横距离对声场的指向性非常明显；随着水平距离的增大，接收器与声源之间正横距离对声场的影响减弱。相同频率不同接收深度的声场也存在差异，这是由于在波导中到达不同接收位置的声线能量不同。