1. 引言
由于浅海波导环境复杂多变,声信号在海洋中传播受其影响很大,海面和海底的反射、散射产生了多途干涉效应,使声波呈多模态传播。因此在估计目标的方位或者深度时,如果可以实现模态的有效分离或得到水平波数,对浅海波导中的目标探测与定位有极大的助益。Vincent E. Premus [1] 等人基于垂直阵对模态进行分离,然而垂直阵的姿态在水中很难保持,应用有很大的局限性。Florent Le Courtois [2] 等人利用水平阵实现了对模态的分离,但是需要阵列的端射方向朝向目标,而线阵的方位估计性能在端射方向最差。因此水平阵的模态分离与方位估计存在着矛盾。
在现代声纳系统中,圆弧阵列因其特有的阵列结构具有360˚均匀的指向性,并且不存在左右舷模糊的缺陷 [3] 。因此在圆阵上应用模态分离,不会与方位估计产生冲突。
2. 均匀圆阵
如图1所示,设平面离散均匀间隔的圆阵半径为r,为方便表示,将阵元按顺时针方向表示为G1···GN,设圆心o通过阵元G1的方向选作0˚方向。阵元的个数为N,则相邻两个阵元之间的夹角为
。设信号从
方向入射,那么到达圆心o的信号假定为Acos2πft。第i个阵元Gi所接收到的信号为
,其中
为Gi相对于圆心o的延时 [4] ,即
(1)
这里c代表水中的声速。c/f用波长λ来表示,其中
(2)
为了在理想方向形成波束,本文将Gi信号加以延时
,得到信号为
![](//html.hanspub.org/file/1-2670123x16_hanspub.png)
Figure 1. The element phase difference of m modal normal mode wave
图1. 第m阶模态简正波的阵元间相位差
(3)
其中,
,
文献 [3] [5] [6] 给出了均匀圆阵的阵列流行导向矢量(即每个阵元的相移矢量):
; (4)
,又因为时间
,因此
,令
,则
(5)
根据波束形成原理 [7] ,当信号传播到各基元时,由于声程差的缘故,每个基元输出的信号是有差异的。我们将这差异进行人为的补偿到每个信号都相同,则叠加后的结果很容易得到目标的方位。类似地,我们在搜索时不搜索方位,而是改成搜索(6)式中k对应的相位,则可以实现圆阵的简正波模态分解。
3. 简正波原理与模态分离
3.1. 简正波原理与模态分离
文献 [8] [9] [10] 给出了简正波的绝热近似理论,可以得出接收声压谱为:
(6)
其中zs和z分别为声源深度和接收阵布放深度,r为声源和接收阵之间的距离,ψn是依赖于深度变化的第n阶简正波的模式函数,
为第n阶简正波的水平波数,
为声源的幅频响应,ρ为介质密度,N为当前频率下的波导中存在的简正波阶数 [11] 。由(6)式可以得出,水平阵各个阵元接收到的声压可以看作是声源产生的各阶简正波在接收点处叠加的结果,各阶简正波的相位随距离的变化是独立的,依赖于每一阶简正波的水平波数变化。因此利用简正波的相位信息,就可以对接收到声压信号进行模态域上的分离。
3.2. 模态域波束形成原理
在远场条件下,第k个阵元接收到的数据与原点之间的时间延时为
。R为圆阵的
半径,c为声源在该水域的声速,φk是第k个阵元与x轴的夹角,
是信号到达角,相位差为
。则第m个阵元所接收到的声压可表示为:
(7)
其中
为频点f处到达第一个阵元的声压。考虑到每阶简正波的传播方式,上式可以改写为:
(8)
根据简正波理论可知,声波在浅海波导中传播时是以简正波的形式进行传播的。一个声信号可以分解为多阶简正波,而简正波之间可以看作是独立传播的。本文把每一阶简正波当作一个独立的声信号,如图1所示,因为每一阶简正波的水平波数不同,所以声波到达水平圆阵的每个阵元的相位延迟τk也是不同的。根据圆阵波束形成原理 [12] ,按照x阶简正波所对应的水平波数进行相位补偿时,其它阶简正波分量并不能同相叠加。首先根据圆阵对目标信号进行方位估计,得到信号方位信息后,本文做波数域上的波束形成,当搜索的波数等于某阶简正波的水平波数,即会出现该阶简正波对应的波束输出的峰值。各阶简正波波束输出之间的幅度比与模态强度比是一致的 [11] 。波束输出
表达式如下:
(9)
其中M为其它模态的残留分量,当阵元数较多时,叠加输出
主要由第x阶的模态分量决定,不考虑其它模态残留量M产生的影响。我们就得到了各阶简正波的能量分布与对应的水平波数,并实现了在波数域上的模态分离。为了便于观察模态分离的效果,本文定义
为模态强度,作为模态域波束形成后的幅度的参考值。式中
为深度z对应的模式函数,经过kracken仿真得到。模态强度B公式如下:
(10)
4. 仿真与实验结果
4.1. 波导模型
本文应用Pekeris波导,基于kracken软件进行仿真。本文获取了简正波模态,以经过声源与水听器连线的所在直线做海洋垂直剖面,该剖面以海面为r轴,过声源做垂直于海平面的z轴,可以得到声源坐标为(0, zs),水听器坐标为(r, z)。浅海波导条件如图2所示,波导深度D = 100 m,海水中的介质密度为ρ = 1000 kg/m3,声速为1500 m/s,海底的介质密度为1500 kg/m3,声速为1600 m/s,海底衰减为0.2 dB/λ。
4.2. 模态分离仿真实验结果及结论
图3所示是差异化条件下所做仿真实验结果,差异化条件在对应图下均有标注。
图3(a)所示是在声源目标频率为100 Hz条件下,圆阵阵元数为256,半径为500 m,入射角度为50˚,
![](//html.hanspub.org/file/1-2670123x40_hanspub.png)
Figure 2. The spreading model of normal mode wave in shallow sea and its simulation conditions
图2. 浅海波导模型及仿真条件
(a) DR = 500 m, Ds = 5 m, DZ = 95 m
(b) R = 1000 m, Ds = 5 m, DZ = 95 m
Figure 3. The simulation conclusion and reference chart of normal mode wave modal separation based on circle array
图3. 圆阵的简正波模态分离仿真结果与参考图
接收阵与声源距离为10 km,声源深度为5 m,接收阵深度为95 m时的经过模态域波束形成后的模态分离结果。图3(b)所示是在声源目标频率为100 Hz条件下,圆阵阵元数为256,半径为1000 m,入射角度为50˚,接收阵与声源距离为10 km,声源深度为5 m,接收阵深度为95 m时的经过模态域波束形成后的模态分离结果。图3(c)是同等条件下经过kracken软件直接仿真得到的结果,作为模态幅度的参考值。仿真结果显示:经过模态分离后的各阶简正波幅度与kracken软件得到的模态强度参考值是对应的,但是第二阶模态从图上观察未与第一阶模态分开。这是因为模态分离要求阵列孔径较大,而第一阶第二阶模态间的波数间隔较小,所以在当前仿真条件下,只能区分1,3,4,5阶简正波。并根据峰值的位置得到1,3,4,5阶水平波数为:0.4167、0.4109、0.4035、0.3958,与参考水平波数0.4179、0.4100、0.4028、0.3938基本上是一致的,均匀圆阵的5阶模态分离效果比较良好。
图3(d)所示是在声源目标频率为100 Hz条件下,圆阵阵元数为256,半径为500 m,入射角度为50˚,接收阵与声源距离为10 km,声源深度为10 m,接收阵深度为50 m时的经过模态域波束形成后的模态分离结果。图3(e)所示是在声源目标频率为100 Hz条件下,圆阵阵元数为256,半径为1000米,入射角度为50˚,接收阵与声源距离为10 km,声源深度为10 m,接收阵深度为50 m时的经过模态域波束形成后的模态分离结果。图3(f)是同等条件下经过kracken软件直接仿真得到的结果,作为模态幅度的参考值。仿真结果显示:经过模态分离后的各阶简正波幅度与kracken软件得到的模态强度参考值是对应的,但是第5阶模态没有被明显分离开来。这是因为模态分离第5阶能量较小受旁瓣影响较大,所以在当前仿真条件下,只能区分1,2,3,4阶简正波。并根据峰值的位置得到1,2,3,4阶水平波数为:0.4178、0.4142、0.4103、0.4023,与参考水平波数0.4179、0.4150、0.4100、0.4028基本上是一致的,均匀圆阵的5阶模态分离效果比较理想。
基于此,我们得出如下结论:本文在已知海洋环境条件、目标的距离范围和目标声源的深度条件下,基于均匀圆阵进行模态域分离,取得了理想的分离效果。在同等条件下,圆阵半径越大,则模态分离效果越好,声源和接收阵的深度差越小,则分离效果越好。基于本文结论,我们就可以通过分离模态进行匹配,从而可以较为准确地估计目标的深度。
致谢
感谢我的导师章新华教授,他严谨细致、一丝不苟的作风一直是我工作、学习中的榜样;他循循善诱的教导和不拘一格的思路给予我无尽的启迪。
感谢我的师兄李鹏,这篇学术论文的每个实验细节和每个数据,都离不开你的细心指导。而你开朗的个性和宽容的态度,帮助我能够很快的融入我们这个实验室。在此,我还要感谢在一起愉快的度过研究生生活的师兄师弟们,正是由于你们的帮助和支持,我才能克服一个一个的困难和疑惑,直至本文的顺利完成。