风成噪声是海洋环境中最广泛存在的环境噪声, 是浅海环境中被动声纳系统关心频段范围内的主要噪声成分和限制声纳性能发挥的主要因素^[1]. 过去几十年里相继开展了一系列有关风成噪声理论建模和特性研究的工作^[2-5], 以期建立合适的数学物理模型对风成噪声的特性进行预报和利用, 以达到降低噪声对声纳工作性能的影响并提升声纳处理性能的目的.

风成噪声建模是风成噪声特性预报的基础, 包括风成噪声源建模和声传播建模两方面. 代表性的风成噪声模型有Cron和Sherman^[6]提出的C-S模型^[6]以及由Kuperman和Ingenito^[7]提出的K-I模型. 前者将风成噪声源建模为分布于无穷大海表上的指向性声源, 后者将其建模为位于水面下方一无穷大平面上的单极子声源, 尽管两种建模方式略有不同, 但可以证明它们在特定条件下是等价的^[8]. C-S模型和K-I模型最主要的区别在于它们对水下声传播的处理方式不同. C-S模型假设声波沿直线传播, 且不与海底发生任何作用, 这种假设使其应用场景较为局限, 通常只适用于均匀深海环境. K-I模型则直接采用波动理论来处理声传播问题, 无论在浅海或深海环境中均适用, 因而也成为后续研究中最广泛使用的风成噪声模型. 基于K-I模型, 对风成噪声的空间相关性、噪声级垂直分布特性、空间指向性等进行了大量分析, 并发现了“噪声凹槽”现象(即在风成噪声垂直方向响应图的水平方向上总是存在一个凹槽), 此外, 还深入分析了声速剖面、底质特性参数等环境参数以及海面起伏、风速、表面声道等对有关噪声特性的影响^[2,7-12].

风成噪声的空间相关特性是阵列处理算法设计的关键. 以往研究通过将基于K-I模型推导得到的风成噪声空间协方差矩阵引入接收信号建模中, 并结合统计阵列信号处理理论推导得到了一系列基于垂直阵接收数据的目标检测、定位及辐射声功率估计算法, 很大程度上提高了风成噪声背景下的被动声纳处理性能^[13-18]. 然而, 这些研究大多关注的只是算法设计本身. 与目标辐射声波相同, 风成噪声也是通过源辐射、波导传播到达接收阵, 因此, 它也具有与前者相类似的多模态结构. 即根据简正波理论^[5,19], 声源激励的声场可以分解为若干阶简正波之和的形式, 各阶简正波对声场的贡献与声源深度上该阶简正波对应的模深函数幅度成正比. 由于风成噪声源位于海面附近, 对应的低阶模态幅度非常小, 因此, 风成噪声很少对应低阶模态、其大部分能量都来源于中高阶模态(这也是“噪声凹槽”出现的本质原因), 尤其是在负声速梯度环境中^[20]. 而另一方面, 信号场对应的模态结构与声源深度息息相关; 随声源深度的变化, 风成噪声场与信号场两者之间的模态结构差异势必发生变化, 进而导致算法性能也出现一定的变化. 这一变化为何? 是否存在一些特定的变化规律? 以上问题的答案对于声纳设备的使用以及目标机动隐蔽具有重要的指导意义, 但至今未见报道.

考虑到无论是对于目标检测或是参数估计, 计入了垂直阵采样声强、背景噪声功率、阵处理增益的阵列信噪比(signal-to-noise ratio, SNR)是决定被动声纳处理性能最关键的物理量^[13-15,21], 本文研究聚焦于风成噪声背景下垂直阵采样SNR随声源深度的变化.本文在K-I模型的基础上, 对典型夏季浅海环境中风成噪声的模态结构进行研究, 探究风成噪声空间协方差矩阵特征向量与垂直阵模态采样之间的对应关系, 在此基础之上, 结合辐射声场的模态结构, 进一步对SNR展开理论推导, 揭示SNR随声源深度的变化规律.

一简谐声源在浅海环境中可以激发出若干阶简正波. 每一阶简正波对应了一阶特征函数和特征值, 两者分别被称为模深函数(后文简称模态)和水平波数. 模态的形状和阶数, 以及水平波数均由海洋环境和声源频率决定^[5].

一个N元垂直阵可以对各阶模态进行采样, 得到的模态信息通常表示为如下$N \times M$维的模态矩阵:

其中, M表示感兴趣频率声源在给定浅海环境中激发的模态数, ${\varPsi _m}(z)$表示第m阶模态($m = 1, $$ 2, \cdots, M$), ${z_n}$为第n个阵元所在深度($n = 1, 2, \cdots, $$ N$), ${\varPsi _m}({z_n})$为第n个阵元对第m阶模态的采样, ${\boldsymbol{\varPsi }}$的第m列$ {{\boldsymbol{\psi }}_m} = {\left( {{\varPsi _m}({z_1}), {\varPsi _m}(z_2), \cdots, {\varPsi _m}({z_N})} \right)^{\rm{T}}} $称为第m阶模态采样.

根据传播距离的远近, 模态可分为传播模态与泄露模态, 前者传播距离远, 对应水平波数虚部小(虚部对应能量衰减); 后者传播距离近, 对应水平波数虚部大. 记传播模态阶数为${M_{\rm{p}}}$, 相应地, ${\boldsymbol{\varPsi }}$中前${M_{\rm{p}}}$列表示传播模态采样, 后$M - {M_{\rm{p}}}$列表示泄露模态采样.

在理论分析中, 假定垂直阵孔径足够大、可以采集到各阶模态的完整信息. 此时, 根据简正波理论, 各阶模态采样近似单位正交^[5,22]:

式中, $ {\delta _{mn}} $表示狄拉克函数, 且仅当$m = n$时, ${\delta _{mn}}$等于1, 否则为0, 上标H表示共轭转置, 水体密度$\rho = 1$ g/cm³及深度采样间隔1 m被省略.

风成噪声通常假设由均匀分布在水面下一无穷大平面上的单极子声源产生^[7]. 对于该平面上的某个圆环微元, 其对应的第i号和第j号阵元上风成噪声的空间互谱密度为$ v{q^2}p(r, {z_i}){p^*}(r, {z_j})2{\text{π}}r{\rm{d}}r $, 其中, v表示噪声源密度, q表示单位噪声源强度, $ p(r, {z_i}) $, $ p(r, {z_j}) $分别表示与垂直阵相距r的微元在深度${z_i}$, ${z_j}$上的声场格林函数, 上标“*”表示共轭. 假定噪声源平面内各微元统计独立, 则第i, j号阵元上风成噪声的互谱密度$ C({z_i}, {z_j}) $为

其中, R为噪声源平面半径, 理论上为无穷大. 将(3)式中的声场格林函数表示成依深度变化的波数积分的形式, 进一步可以推导得到$ C({z_i}, {z_j}) $的近似理论表达式^[7]:

其中, ${z_0}$为噪声源平面所在深度, 一般不大于波长的四分之一, $ {\varPsi _m}({z_0}) $为噪声平面深度上第m阶模态的幅值, $ {\kappa _m} $和${\alpha _m}$分别表示第m阶模态对应水平波数的实部与虚部, 即衰减系数.

按照(4)式将风成噪声协方差矩阵${{\boldsymbol{C}}_{{\rm{sn}}}}$中的各元素进行展开(下标“sn”为“surface noise”的缩写), 可以得到${{\boldsymbol{C}}_{{\rm{sn}}}}$的模态域形式:

其中, ${{\boldsymbol{\varLambda }}_{\boldsymbol{\varPsi }}}$表示M维的对角阵

注意到各阶模态采样在采样完整情况下是正交归一化的, 因此, (5)式恰好对应了风成噪声协方差矩阵${{\boldsymbol{C}}_{{\rm{sn}}}}$的特征分解形式:

其中, ${\boldsymbol{\varLambda }} = {\rm{diag}}({\lambda _1}, {\lambda _2}, \cdots, {\lambda _M})$为由M个非零特征值${\lambda _m}$按照由大到小顺序排列而成的对角阵(${\lambda _1} > $$ {\lambda _2} > \cdots > {\lambda _M}$), ${\boldsymbol{U}} = ({{\boldsymbol{u}}_1}, {{\boldsymbol{u}}_2}, \cdots, {{\boldsymbol{u}}_M})$为$N \times M$维的特征矩阵, 第m列${{\boldsymbol{u}}_m}$表示与特征值${\lambda _m}$相对应的$N \times 1$维特征向量. 可见, ${{\boldsymbol{C}}_{{\rm{sn}}}}$的每一阶特征向量对应了一阶模态采样, 对应的特征值则与噪声源平面深度上该模态幅度的强度${\left| {{\varPsi _m}({z_0})} \right|^2}$成正比. 在负声速梯度环境中, 由于低阶模态在水面附近的幅值趋于0, 而中高阶模态更为显著, 因此, ${{\boldsymbol{C}}_{{\rm{sn}}}}$的小特征值对应的特征向量主要表征低阶模态采样, 而大特征值对应的特征向量则主要表征中高阶模态采样.

例如, 图1给出了一典型的负声速梯度浅海三层波导, 包括100 m深的水体、2.5 m厚的沉积层和半空间. 海面声速为1500 m/s, 海底声速为1480 m/s, 声速随深度线性变化; 沉积层和半空间为等声速层, 声速分别为1650 m/s和1700 m/s, 两层的介质密度和压缩波衰减系数均相同, 分别为1.80 g/cm³和0.35 dB/λ. 采用KRAKENC声场计算软件^[23], 可以计算得到200 Hz声源在该波导中可以激发出18阶模态, 其中前13阶为传播模态、后5阶为泄露模态. 图2给出了这18阶模态的形状. 可以看到, 低阶模态在水面附近的幅值趋近于0, 而中高阶模态幅值则较为显著, 如在2 m深度上, 前3阶模态的幅值依次为0.0003, 0.0036, 0.01, 而第8, 9, 10阶模态幅值分别为0.06, 0.07, 0.08.

图3所示为垂直阵覆盖全水深时(孔径100 m、阵元间隔1 m)${{\boldsymbol{C}}_{{\rm{sn}}}}$的归一化特征值及其特征向量与各阶模态采样的相关系数$\left| {{\boldsymbol{u}}_m^{\rm{H}}{{\boldsymbol{\psi }}_n}} \right|$. 噪声协方差矩阵中各元素由(3)式计算得到, 其中, $R = 100 \;{\rm{km}}$, ${z_0} = 1.875\;{\rm{m}}$(距海面四分之一个波长). 从图3(a)和(b)可以看到, 该矩阵包含有趋近于0的小特征值, 与其对应的特征向量表征低阶模态采样, 并且, 由于阶数越低的模态在水面附近幅值越小, 因此, 特征值越小的特征向量表征的也是阶数越低的模态采样, 如最小、次小和第13阶特征值对应的特征向量分别表征1, 2, 3阶模态采样, 而大特征值对应的特征向量表征的则是第8—11阶模态采样.

由于在损耗型海底中模态一般都存在拖尾, 而低阶模态的拖尾可以忽略, 高阶模态则较为显著, 因此, 覆盖全水深的垂直阵只能采集到前者的完整信息, 相应地, 如图3(b)所示, 低阶模态采样与特征向量之间几乎一一对应, 如前5阶, 而中高阶模态采样则可能对应多个特征向量. 此外, 尽管泄露模态在水面附近的幅值也较大, 但由于这部分模态在传播过程中衰减快、对风成噪声的贡献较小(即(6)式中衰减系数${\alpha _m}$显著大于其余传播模态), 因此, 这部分模态采样表征的也主要是${{\boldsymbol{C}}_{{\rm{sn}}}}$小特征值对应的特征向量.

根据水声信道传播理论, 声源、海洋声信道、阵元在一定时间范围内可看作是一个分置式的线性时不变系统: 输入为声源辐射声信号, 输出为阵元采样的辐射声信号, 系统传递函数为阵元所在空间位置上由辐射声源激发的声场格林函数, 其是声源频率f、声源位置S(包括声源深度${z_{\rm{s}}}$和声源距离${r_{\rm{s}}}$)和环境参数集${\boldsymbol{\varPhi }}$(包括水体深度、声速剖面、地声参数等)的函数. 记第n号阵元所在空间位置上的声场格林函数为${g_n}(f, {\boldsymbol{S}}, {\boldsymbol{\varPhi }})$, 相应地, $N \times 1$维的垂直阵采样声场频域快拍$ {\boldsymbol{r}}(f) $可以表示为

其中, $s(f)$表示辐射声信号复幅度; ${\boldsymbol{G}}(f, {\boldsymbol{S}}, {\boldsymbol{\varPhi }}) \triangleq $$ \left\{ {{g_n}(f, {\boldsymbol{S}}, {\boldsymbol{\varPhi }})} \right\}_{n = 1}^N$为$N \times 1$维的格林函数向量, 也称为信号波前; ${\boldsymbol{n}}(f)$表示$N \times 1$维的背景噪声向量, 假设服从均值为0、协方差阵为${{\boldsymbol{C}}_{{\rm{nn}}}}(f)$的复高斯分布, 即${\boldsymbol{n}}(f)\sim\mathcal{C}\mathcal{N}({\rm{0}}, {{\boldsymbol{C}}_{{\rm{nn}}}}(f))$ ^[15-18].

本文考虑简谐声源以及垂直阵位于声源远场的情况. 根据简正波理论, 并省略有关参数对频率的依赖, 水平不变浅海波导中各阵元上的声场格林函数可以表示为一系列简正波之和的形式^[5]:

由于泄露模态经过远程传播后对声场的贡献几乎为零, 因此(9)式中仅保留传播模态. 将信号波前中的所有元素均表示为(9)式的形式, 可以得到:

其中, ${\boldsymbol{d}} = {({d_1}, {d_2}, \cdots, {d_{{M_{\rm{p}}}}})^{\rm{T}}}$表示${M_{\rm{p}}} \times 1$维的模态系数向量:

(10)式和(11)式表明, 垂直阵采样到的目标辐射声信号为各阶模态采样的线性组合, 各阶模态采样对信号波前的贡献与声源深度上该模态的幅值成正比.

背景噪声通常由风成噪声与水听器噪声两部分构成, 由于彼此不相关, 因此, 背景噪声协方差矩阵可以表示为各自协方差矩阵之和的形式:

其中, ${{\boldsymbol{C}}_{{\rm{hn}}}} \triangleq \sigma _{\rm{w}}^2{{\boldsymbol{I}}_N}$为水听器噪声空间协方差矩阵, 即各水听器自噪声互不相关且功率$\sigma _{\rm{w}}^2$相等, 下标“hn”为水听器噪声“hydrophone noise”的缩写; 风成噪声协方差阵由(3)式给出. 为衡量风成噪声与水听器噪声的相对比重, 定义:

其中, ${\rm{tr}}( \cdot )$表示矩阵的迹, α表示的是风成噪声功率在各阵元上的平均与水听器噪声功率的比值. 随着声纳技术的不断发展, 当前水听器自噪声已大幅下降、声纳系统背景噪声大多以环境噪声为主. 例如, 某型水听器的自噪声谱级在200 Hz附近约为38 dB (参考值: 1 ${\rm{μ P}}{{\rm{a}}^2}/{\rm{Hz}}$)^[24], 比0.5级海况下的表面噪声谱级约低10 dB、比8级海况条件下的表面噪声谱级约低40—50 dB^[25] (α大致在10⁴—10⁵量级). 当海况进一步升高, α的取值还将进一步增大.

对于(8)式给出的垂直阵采样声场, SNR定义如下^[13,21]

其中, ${\boldsymbol{G}} \triangleq {\boldsymbol{G}}({\boldsymbol{S}}, {\boldsymbol{\varPhi }})$, $ {\left| s \right|^2}{{\boldsymbol{G}}^{\rm{H}}}{\boldsymbol{G}} $为垂直阵采样的感兴趣目标辐射声总声强. (14)式给出的SNR不仅计入了垂直阵采样声强、背景噪声功率, 同时还计入了阵处理增益, 其本质上是波束形成器所能达到的最大输出信噪比^[21]、是决定阵处理性能最关键的物理量; 对于匹配滤波器、模态空间检测器、能量检测器以及声源位置已知时的最大似然辐射声功率估计器等, 它甚至是唯一的性能决定因素^[13-15]. 一般情况下, SNR越大, 阵处理性能越好, 如对目标的检测概率更高、参数估计的均方误差越小等. 本研究的目的是结合风成噪声和信号波前的模态结构, 研究负声速梯度环境中SNR随声源深度的变化规律. 注意到${\left| s \right|^2}$为常数, 对此研究没有影响.

将(12)式代入(14)式中, 并记${{\bar{\boldsymbol{C}}}_{{\rm{sn}}}} = {{\boldsymbol{C}}_{{\rm{sn}}}}/{\rm{tr}}({{\boldsymbol{C}}_{{\rm{sn}}}})$, 可得:

将${\bar{\boldsymbol{C}}_{{\rm{sn}}}}$替换为其特征分解的形式, 即${\boldsymbol{U\bar \Lambda }}{{\boldsymbol{U}}^{\rm{H}}}$, 其中${\boldsymbol{\bar \varLambda }} \;=\; {\boldsymbol{\varLambda }}/{\rm{tr}}({\boldsymbol{\varLambda }}) \;=\; {\rm{diag}}({\bar \lambda _1}, \;{\bar \lambda _2}, \cdots,\; {\bar \lambda _N})$, ${\bar \lambda _n} = $$ {\lambda _n}/{\rm{tr}}({\boldsymbol{\varLambda }})$, 并根据矩阵求逆引理, 可以推导得到:

其中省略了常数因子$1/\sigma _{\rm{w}}^2$. 由2.2节分析可知, 风成噪声协方差阵的各特征向量对应各阶模态采样, 因此, 信号波前在这些特征向量张成空间上的投影应等于其本身(${\boldsymbol{U}}{{\boldsymbol{U}}^{\rm{H}}}{\boldsymbol{G}} \approx {\boldsymbol{G}}$), 即${{\boldsymbol{G}}^{\rm{H}}}{\boldsymbol{U}}{{\boldsymbol{U}}^{\rm{H}}}{\boldsymbol{G}} \approx $$ \left\| {\boldsymbol{G}} \right\|_2^2$. 将其代入(16)式替换被减数, 并同时利用(10)式对信号波前进行展开:

由于模态采样与特征向量之间是一一对应的, 且各阶模态采样之间是正交归一化的, 因此, (17)式可以化简为

其中, 下标${n_m}$表示第m阶模态采样对应的风成噪声协方差矩阵特征值的序号.

(18)式表明, 风成噪声背景下的SNR可以表示为各阶模态系数强度$ |{d_m}{|^2} $的加权求和, 而各模态系数强度对应的权值:

与该模态采样对应的风成噪声协方差矩阵${{\boldsymbol{C}}_{{\rm{sn}}}}$的特征值成反比. 如2.2的分析, 在负声速梯度环境中由于低阶模态采样对应的特征值趋近于0, 中高阶模态对应的特征值更为显著, 因此, 低阶模态对应的权值更高、而中高阶模态对应的权值则相对很小. 例如, 如后文所示, 对于200 Hz频率、图1所示的典型夏季浅海波导, 1阶和2阶模态对应的权值远大于其余模态. 进一步将(11)式代入(18)式中可得:

从(20)式可总结得到SNR随声源深度的变化具有如下规律:

1) SNR随声源深度的变化可以近似为低阶模态幅度强度$ |{\varPsi _m}(z){|^2} $随深度变化的线性叠加, 且模态阶数越低、贡献越大. 需要指出的是, 尽管当声源位于海面附近时中高阶模态对SNR事实上也具有一定的贡献, 但由于这部分SNR相比其他声源深度上的SNR(尤其是水体下半部分)非常小, 因此, 这部分模态对SNR随声源深度的总体变化规律的影响很小.

2)由于低阶模态在水面附近的幅值非常小、在水体下半部分较大, 特别是1阶模态, 因此, 在同一距离条件下, 相比放置于海面附近的情况, 声源放置于水体下半部分时SNR总是更高(放置于不同深度时声源级相同); 并且, 由于更高阶的模态在1阶模态峰值点下方均具有局部最大值点, 因此, SNR最大值出现的声源深度通常对应了1阶模态峰值点下方某一深度.

3)当声速负梯度较强时, 可以验证, 1阶模态采样对应的风成噪声协方差矩阵特征值一般比2阶模态采样对应的特征值至少小2—3个数量级. 此时, 若风成噪声也非常强, 如特别地, $\alpha \to + \infty $, 可以得到:

将其代入(20)式, 有:

即, 此时SNR随声源深度的变化近似可由1阶模态幅度强度随深度的变化独立表征. 相应地, 最大SNR出现的声源深度即为1阶模态峰值点所在深度.

以上分析假定了所有模态采样与风成噪声协方差矩阵特征向量一一对应, 而由2.2节分析可知, 在损耗型海底的浅海环境中这对中高阶模态采样较为苛刻. 但由于这部分模态对SNR随声源深度的总体变化规律影响很小, 因此, 以上结论在该环境中仍然适用.

本节在典型负声速梯度浅海环境中对SNR随声源深度的变化规律进行仿真研究. 仿真中使用的浅海环境、声源频率与2.2节保持一致. 首先在垂直阵覆盖全水深、对各阶模态采样较为完整的情况下(阵列孔径为100 m、阵元间距为1 m)开展研究, 而后分析风成噪声强度、水体声速梯度、阵列孔径带来的影响. 不同声源位置处的SNR均由其定义(14)式计算得到, 信号波前由KRAKEN声场软件计算得到^[23]. 如无其他说明, α设定为10, 对应于风成噪声平均功率比水听器自噪声功率高10 dB的情况. 为体现SNR随声源深度的变化规律, 不同声源距离的SNR均关于该距离上的最大SNR进行归一化处理.

图4(a)和(b)分别给出了不同声源距离上归一化SNR随声源深度的变化关系以及${r_{\rm{s}}} = 10\;{\rm{km}}$上的结果. 由图4可知, SNR随声源深度的变化在不同距离上近似相同(这并不说明SNR本身与声源距离无关, 由(20)式可知, SNR与r_s总体成反比), 当声源深度由1 m增大到85 m时, 归一化SNR由0.0007增大到1, 而后随声源深度的增大迅速减小. 在不同距离上, 最大SNR均出现在声源位于85 m深度附近, 该深度略大于1阶模态峰值点对应深度82 m; 并且当声源位于水体下半部分时SNR总是高于其位于海面附近的情况, 例如, 当${z_{\rm{s}}} \leqslant 20\;{\rm{m}}$时的SNR (归一化值不超过0.063)至少要比声源位于85 m深度附近时低12 dB ($10{\log _{10}}(1/0.063)$).

为了更好理解以上SNR随声源深度的变化结果, 图5(a)给出了风成噪声协方差矩阵各特征向量关于SNR的归一化权值(关于最小特征向量/1阶模态采样对应权值进行了归一化). 结合图3可知, 前3阶模态采样对应的归一化权值分别为1, 0.4, 0.05, 而其余特征向量/模态采样对应的权值几乎可以忽略, 相应地, SNR随声源深度的变化可由前2阶模态幅度强度随深度的变化来表征(最小两阶特征向量). 图5(b)给出的各特征向量对SNR的贡献$ {w_m}\left\| {{{\boldsymbol{G}}^{\rm{H}}}{{\boldsymbol{u}}_m}} \right\|_2^2 $随声源深度的变化关系(${r_{\rm{s}}} = $$ 10\;{\rm{km}}$)对此进行进一步的验证. 对比图2和图5(b)也可以看到, 最小两阶特征向量对SNR贡献随声源深度的变化分别与1阶和2阶模态函数的形状相吻合.

根据上述分析, 图4所示结果不难理解: 1)可以验证, 前两阶模态对应的指数衰减${{\rm{e}} ^{ - 2{\alpha _m}{r_{\rm{s}}}}}$在不同声源距离上近似保持相对不变, 因此, 由(20)式可知, SNR随声源深度的变化在不同声源距离上近似相同; 2)尽管1阶模态幅度在0—20 m深度范围内随深度变化不大、趋近于0, 但由于第2阶模态幅度在此范围内随深度存在一定程度的增大, 因此, SNR也将随声源深度在此范围内增大而增大; 随后, 在两阶模态的共同作用下, SNR随声源深度的增大而持续增大、并在深度增大到1阶模态峰值点下方时, 即85 m深度上, 达到最大; 在靠近海底的深度上, 由于1, 2阶模态幅度均随深度增大迅速减小, 因此, SNR随声源深度的增大也迅速减小.

图6(a)进一步给出了α取20, 50, 100, 1000时归一化SNR随声源深度的变化结果. 根据(19)式可以计算得到, 风成噪声越强将导致关于SNR的权值越来越集中在1阶模态上, 因此, 在这几种α取值情况下SNR随声源深度的变化规律也与声源距离近似无关, 这里只给出${r_{\rm{s}}} = 10\;{\rm{km}}$的结果. 由于更高阶模态对SNR的贡献是SNR随声源深度在海面附近增大而增大以及使最大SNR对应声源深度偏离1阶模态峰值点的主要原因, 因此, 从图6(a)可以看到, α由20增大到1000, 声源位于海面附近时的归一化SNR逐渐变小, 最大SNR对应的声源深度也由84 m减小至82 m. 特别当$\alpha = 1000$时, 最大SNR对应的声源深度恰好是1阶模态峰值点所在深度(82 m). 这是因为此时SNR随声源深度的变化可由1阶模态幅度强度随深度的变化近似独立表征(前两阶模态采样对SNR的权值分别为1, 0.01): 图6(b)对两者进行了对比, 可以看到它们是近乎重合的. 进一步增大α将不会对以上结果产生明显影响.

需要指出的是, 当风成噪声较弱时, 部分中高阶模态对SNR将产生较大的贡献(权值变大). 相比低阶模态, 由于这部分模态随声源距离的增大衰减更快, 因此, 此时SNR随声源深度的变化规律在不同声源距离上将略有不同. 图7给出了$\alpha = 1$时不同声源距离上归一化SNR随声源深度的变化结果(实际中α通常大于此值, 这里主要关注α变化对SNR随声源深度变化规律的距离无关性的影响趋势). 可以看到, 由于中高阶模态具有幅值变化频率高、在水面附近幅值增长快的特点, 因此, 在距离声源较近的情况下垂直阵采样SNR随声源深度在海面附近增大而增大的速度要比距离较远时更快, 同时还伴随着一定的振荡.

水体声速梯度定义为海底声速与海面声速的差值与水深之比, 记为g, 其表征声速随深度变化的快慢, 很大程度上决定了模态的形状. 本节研究水体声速梯度对SNR随声源深度变化的影响. 固定海面声速为1500 m/s, 海底声速分别取1495, 1493, 1490, 1470 m/s, 以此生成4种不同的负声速梯度环境: $g = - 0.05$, $ - 0.07$, $ - 0.10$和–0.30 s^–1.

图8给出了这4种环境中不同声源距离上SNR随声源深度的变化结果. 从中可以观察到与图4(a)相类似的结果, 但由于在弱负梯度环境中低阶模态与中高阶模态在水面附近的幅值差异并不大、部分中高阶模态对SNR也存在一定的贡献, 因此, 如图8(a)所示, 当g = –0.05 s^–1时, SNR随声源深度的变化在不同声源距离上有一定差异, 但对比图8(b)—(d)可知, 随着负梯度的增强, 由于低阶模态与中高阶模态在水面附近的幅值差异越来越大, SNR随声源深度变化将主要受低阶模态的影响, 而这部分模态对应的指数衰减项在不同声源距离上通常差别不大、甚至近似相同, 因此, SNR随声源深度的变化将呈现出越来越好的距离无关性.

图9对这4种环境中${r_{\rm{s}}} = 10\;{\rm{km}}$上SNR随声源深度的变化结果进行对比. g由–0.05 s^–1减小到$ - 0.30\;{{\rm{s}}^{ - 1}}$, 声源位于海面附近时的归一化SNR越来越小、最大SNR相对其优势越来越明显, 同时, 最大SNR出现的声源深度也越来越大, 由71 m增大到89 m. 这是因为声速负梯度越强, 低阶模态在海面附近的幅值将越趋近于0, 同时, 1阶模态峰值点所在深度也将越大, 在g = –0.05 s^–1, g = –0.07 s^–1, g = –0.10 s^–1和g = –0.30 s^–1环境中该深度依次为69, 73, 76, 85 m.

注意到低阶模态在海面附近的幅值越趋近于0, 同时也意味着越来越多的低阶模态将对SNR产生更大的影响(特征值越趋近于0、权值越大), 因此, 随着负梯度的增强, 最大SNR对应的声源深度将越来越偏离1阶模态峰值点对应深度, 例如, 在g = –0.05 s^–1时, 最大SNR对应声源深度偏离1阶模态峰值点对应深度$71 - 69 = 2\;{\rm{m}}$, 而在g = –0.20 s^–1和g = －0.30 s^–1两种情况下, 这一深度偏离分别为$85 - 82 = 3\;{\rm{m}}$和$89 - 85 = 4\;{\rm{m}}$.

前述分析均假设垂直阵孔径足够大、覆盖全水深, 本节分析孔径不足、阵列仅覆盖部分水体带来的影响. 不失一般性, 考虑g = –0.20 s^–1.

图10给出了阵列覆盖10—100 m时不同声源距离上归一化SNR随声源深度的变化结果. 对比图10和图4(a)可以发现, 该结果与垂直阵覆盖全水深时近似相同. 事实上, 可以验证, 此时1阶和2阶模态采样分别对应了风成噪声协方差矩阵最小特征值和第16阶特征值对应的特征向量, 这两阶模态采样之所以在此情况下仍能与特征向量一一对应, 是因为它们在水面附近的幅值非常小, 即便垂直阵没有覆盖这部分水体也仍能采集到这些模态的完整信息, 而这两阶模态对SNR的贡献与阵列覆盖全水深时近似相同—— 归一化权值分别为1和0.40——其余模态采样/特征向量对SNR的贡献均非常小. 相应地, 在此阵列布放条件下SNR随声源深度的变化与阵列覆盖全水深时理应相同.

图11进一步给出了当阵列仅覆盖水体下半部分时(50—100 m)的结果. 从中仍可以观察到与前述类似的结果, 但不同的是, 此时SNR-${z_{\rm{s}}}$关于声源距离呈现近似周期性的变化趋势、最大SNR对应的声源深度在1阶模态峰值点(82 m)所在深度附近振荡.

图12给出了此时风成噪声协方差矩阵特征向量与各阶模态采样的对应关系以及各特征向量关于SNR的归一化权值(关于第11阶特征向量对应的权值做归一化). 由于阵列孔径小到一定程度时各阶模态采样之间将具有近似的复共线性, 即某几阶模态采样可由其余模态采样近似线性表示^[26,27], 因此, 此时所有模态采样可由前11阶特征向量近似表征. 其中, 第11阶、10阶特征向量对SNR具有最大的权值、且同时对应了前两阶模态, 相应地, 这两阶特征值对SNR的贡献均为两阶模态系数的相干求和, 如对于第11阶特征向量, 有:

包含于模态系数中复指数项${{\rm{e}}^{ - {\rm{j}}{k_{{\rm{r}}, m}}{r_{\rm{s}}}}}$的叠加即是最大SNR对应声源深度关于声源距离呈现出近似周期性变化的原因.

以上分析仅讨论了垂直阵覆盖下半部分水体的情况. 当阵列覆盖上半部分水体时, 由于低阶模态主要信息丢失、风成噪声协方差矩阵的特征向量与模态采样之间不再具有上述对应关系, 因此, 以上SNR随声源深度的变化规律也将消失.

本文研究了受风成噪声影响, 典型夏季浅海环境中垂直阵SNR随声源深度的变化关系: 1)在一定条件下, SNR随声源深度的变化与声源距离近似无关; 2)当风成噪声越强时(一定范围内), 同一距离上最大SNR相对于声源位于海面附近时的SNR优势越明显、最大SNR对应的声源深度将越靠近1阶模态峰值点, 并且SNR随声源深度的变化规律将呈现出越来越好的距离无关性; 3)当声速负梯度越强时, 上述SNR优势也越明显、SNR随声源深度的变化规律也将呈现出更好的距离无关性, 但最大SNR对应的声源深度将越来越偏离1阶模态峰值点; 4)该变化规律尽管是在垂直阵覆盖全水深条件下研究得到的, 但事实上只要垂直阵覆盖水体下半部分、对低阶模态采样较为完整时即可观察到. 只不过当阵列孔径小到一定程度时, 模态采样之间的近似复共线性将导致该变化规律随声源距离呈现出近似周期性的变化, 而最大SNR出现的声源深度也将在1阶模态峰值点对应深度附近作周期振荡.