0 引　言

水下目标辐射噪声的测量对水下移动平台的发展至关重要，是发展低辐射噪声水下移动平台的重要参考条件。对水下目标辐射噪声进行测量，需保证测量阵列在一定带宽范围内有恒定的主瓣宽度，且在浅海条件下进行测量时，还需考虑到浅海波导界面反射能量对测量结果的影响。随着目标辐射噪声的逐年降低，测量系统经历了单阵元、垂直阵到体积阵的不断演变^[1-3]。美国等国家率先将双锥阵等体积阵结构应用到噪声测量中，双锥阵具有较宽的水平向主瓣宽度和较窄的俯仰向主瓣宽度，有利于确保目标在波束主瓣内同时抑制界面反射的干扰，能较好满足噪声测量的需求^[4]。

恒定束宽波束形成器是宽带噪声测量的基础^5]。现有方法利用嵌套阵、聚焦变换等方法实现恒定束宽，但大多应用于均匀线阵。应用于双锥阵的嵌套阵方法对阵元数量要求高，布放难度大。本文通过建立双锥阵噪声测量的信号模型，利用双锥阵作为测量阵列，同时引入一种低旁瓣稳健的恒定束宽波束形成方法，有效解决噪声测量中的多径干扰抑制和二维稳健恒定束宽波束形成问题。

1 双锥阵稳健宽带波束形成 1.1 双锥阵噪声测量信号模型

在浅海波导环境中，双锥阵列第 $l$ 个阵元接收到的频率为 $f$ 的复声压信号可以表示为：

$ {p_l}\left( f \right) = S\left( f \right)G\left( {{r_l},{d_l},{z_s},f} \right) + n\left( f \right)。$

(1)

式中： $S\left( f \right)$ 为复声源幅度； $ G\left( {{r_l},{d_l},{z_s},f} \right) $ 为频域格林函数； $ n\left( f \right) $ 为噪声； ${r_l}$ 和 ${d_l}$ 分别为声源与第 $l$ 个阵元之间的距离和阵元深度； $ {z_s} $ 为声源深度。根据射线理论，浅海波导的格林函数可表示为：

$ g\left( {{r_l},{d_l},{z_s},t} \right) = \sum\limits_{i = 1}^N {{A_i}\left( {{r_l},{d_l},{z_s}} \right)\delta } \left[ {t - {\tau _i}\left( {{r_l},{d_l},{z_s}} \right)} \right] 。$

(2)

式中； $ {A_i} $ 和 ${\tau _i}$ 分别为第 $i$ 个特征声线的幅度和传播时间； $\delta \left[ \cdot \right]$ 为冲激函数。不同多途声线的到达时间和到达角度不同，在噪声测量中，需保留直达波而抑制其他多途信号。

考虑如图1所示的双锥阵结构，假设接收阵列最上方阵元深度大于声源深度且声源位于 $y$ 轴方向，则每个阵元与辐射声源之间的相对距离及每个阵元的接收深度可表示为：

$ {d_l} = {d_0} - {z_l} ，$

(3)

$ {r_l} = \sqrt {x_l^2 + y_l^2 + r_0^2 - 2\sqrt {\left( {x_l^2 + y_l^2} \right){r_0}} \cos {\theta _l}}。$

(4)

式中： $ {\theta _l} = \arctan \left( {{{ - {y_l}} \mathord{\left/ {\vphantom {{ - {y_l}} {{x_l}}}} \right. } {{x_l}}}} \right) $ ； ${r_0}$ 和 ${d_0}$ 分别为参考点的距离和深度； $\left\{ {{x_l},{y_l},{z_l}} \right\}$ 分别为各阵元相对参考点的坐标。

1.2 二维稳健恒定束宽波束形成

为了在俯仰方向有较低的旁瓣，在水平方向同时具有较宽的3 dB主瓣宽度和较低的旁瓣，本文将低旁瓣稳健波束形成方法扩展到二维场景中，通过分别设计俯仰角和水平角的旁瓣约束达到旁瓣控制的目的。

根据图1所示的双锥阵结构，假设俯仰角 $\theta $ 为入射方向与 $z$ 轴的夹角，方位角 $\phi $ 为入射方向在 $xy$ 平面投影后与 $y$ 轴的夹角，则双锥阵各阵元相对参考阵元的时延差可表示为：

$ \Delta {\tau _l} = \frac{{u \cdot {e_l}}}{c}。$

(5)

式中： $u = \left( {\sin \theta \cos \phi ,\sin \theta \sin \phi ,\cos \theta } \right)$ 为来波方向的单位方向向量； $ {e_l} = \left( {{x_l},{y_l},{z_l}} \right) $ 为第 $l$ 个阵元的位置坐标； $c$ 为参考声速。双锥阵对每个方位的向量为：

$ a\left( {\theta ,\phi } \right) = {\left[ {{e^{ - j2{\rm{\pi}} f\Delta {\tau _1}}},{e^{ - j2{\rm{\pi}} f\Delta {\tau _2}}}, \cdots {e^{ - j2{\rm{\pi}} f\Delta {\tau _L}}}} \right]^{\text{T}}}。$

(6)

对于二维低旁瓣稳健波束形成方法，通过观察方向无失真响应、最大旁瓣值最小以及权重向量的范数约束，实现低旁瓣稳健波束形成，保证其他方位的多途干扰被有效抑制，降低其他方位噪声源对目标噪声测量的影响，从而提高测量准确度。由于水平角和俯仰角主瓣宽度和旁瓣抑制的需求存在差异，需要对2 个方向的主瓣和旁瓣分别设计。因此，求解波束形成加权向量 $w$ 的约束优化问题可表示为：

$ \begin{gathered} \mathop {\min }\limits_w \mu ，\\ {\rm{ s.t.}}{\text{ }}a\left( {\theta _0},{\phi _0} \right)^Hw = 1 ，\\ {\text{ }}\left\| a{{\left( {\theta _i},\phi \right)}^H}w \right\| \leqslant \mu ,\left( \left| {{\theta _i} - {\theta _0}} \right| > {\Delta _\theta } \right) ，\\ {\text{ }}\left\| {a{{\left( {\theta ,{\phi _j}} \right)}^H}w} \right\| \leqslant \frac{\mu }{\lambda },\left( {\left| {{\phi _j} - {\phi _0}} \right| > {\Delta _\phi }} \right) ，\\ {\text{ }}\left\| w \right\| \leqslant \delta 。\\ \end{gathered} $

(7)

通过引入额外的权重参数 $\lambda $ ，实现俯仰角和水平角的旁瓣分别设计，以满足噪声测量的需求。

该二维低旁瓣稳健波束形成方法在设计频率处有很好的旁瓣抑制效果，且水平角和俯仰角的主瓣宽度均能满足要求，但在其他频率处的效果不好。通过将设计频率处的波束图作为参考波束，利用聚焦变换实现恒定束宽。

聚焦变换利用一个聚焦矩阵 $ T $ 将不同频率的信号数据聚焦到参考频率 $ {f_0} $ 上。本文采用旋转信号子空间方法作为聚焦处理器，使得聚焦后的方向响应向量与参考频率处的误差最小：

$ \begin{gathered} \mathop {\min }\limits_{{T_{{f_k}}}\left( {\theta ,\phi } \right)} {\left\| {{a_{{f_0}}}\left( {\theta ,\phi } \right) - {T_{{f_k}}}\left( {\theta ,\phi } \right){a_{{f_k}}}\left( {\theta ,\phi } \right)} \right\|_F}，\\ {\rm{s.t.}}{\text{ }}{T_{{f_k}}}^H\left( {\theta ,\phi } \right){T_{{f_k}}}\left( {\theta ,\phi } \right) = I 。\\ \end{gathered} $

(8)

因此，不同频率处的波束图分别为：

$ {B_{{f_0}}}\left( {\theta ,\phi } \right) = \left| {a_{{f_0}}^H\left( {{\theta _0},{\phi _0}} \right){w_{{f_0}}}\left( {\theta ,\phi } \right)} \right| ，$

(9)

$ {B_{{f_k}}}\left( {\theta ,\phi } \right) = \left| {a_{{f_k}}^H\left( {{\theta _0},{\phi _0}} \right){T_{{f_k}}}\left( {\theta ,\phi } \right){w_{{f_k}}}\left( {\theta ,\phi } \right)} \right|。$

(10)

结合两式即可获得恒定束宽的波束输出。

2 仿真分析 2.1 测试条件

为实现真实浅海波导环境下对200～1000 Hz范围内的宽带信号的测量，基于BELLHOP模型，对双锥阵恒定束宽方法进行测试。

环境参数如下：海水深度150 m，测量距离200 m，接收深度100 m，目标深度50 m。假设水下目标直径12 m，长度100 m，测量距离200 m，则垂直方向 $ - 3$ dB主瓣宽度最小要求 $ \arcsin \left( {6/200} \right) \times 2 = 3^\circ $ ，水平方向 $ - 3$ dB主瓣宽度最小要求 $ \arctan \left( {50/200} \right) \times 2 = 28^\circ $ 。垂直主瓣要求较窄，容易满足，而水平主瓣要求较宽，需要进行波束设计。此外，浅海波导中声线传播情况如图2所示，存在较强的界面反射信号，直达波、海面反射波和一次海底反射波到达角度分别为 $76^\circ $ 、 $52^\circ $ 、 $127^\circ $ ，需要对除直达波外的多途信号进行抑制。

2.2 稳健恒定束宽波束形成方法

考虑如图3所示的双锥阵结构，共9层8列，层高为0.5d，最大半径为3d，d为最小频率对应的半波长。假定波束形成器的主瓣指向直达波方向 $76^\circ $ ，处理的频率为200～1000 Hz。

图4为本文方法在俯仰方向和水平方向的宽带波束图，图5为各频率处波束图的叠加，不同频率处的波束图与参考频率处基本一致，旁瓣低于 $ - 40 $ dB，能有效抑制多径信号，且水平方向 $ - 3$ dB主瓣宽度为 $ 29^\circ $ ，满足结构噪声源测量需求。

2.3 水下目标辐射噪声测量性能分析

仿真测量信号为声源信号经BELLHOP仿真信道传输后的阵列信号，信噪比为5 dB。声源信号波形和频谱及浅海波导信道的时频特征如图6所示，可以看到，在浅海波导中，受界面反射的影响，信道特征在时域上呈现明显的多途特性，而在频域则体现为周期震荡，造成单水听器接收信号失真。

图7为浅海波导下，多频信号的双锥阵恒定束宽波束输出结果。可以看到，在浅海波导中，受界面反射的影响，单阵元信号的时频波形都存在明显的失真。常规波束形成方法的主瓣宽度随处理频率增加而减小，因此在波束输出中体现出明显的低通滤波器效应，且随频率变化有明显起伏，表明其不能有效滤除界面反射信号的干扰。本文提出的测量方法能够保证信号波形不失真，因此更准确地测量噪声能量。单水听器、常规波束形成以及本文方法的测量结果误差分别为7.1 dB、4.1 dB、0.7 dB，表明该方法在噪声测量中的可靠性。

3 结　语

本文分析了浅海波导对噪声测量的影响，建立了浅海波导环境中的双锥阵信号模型。通过引入低旁瓣稳健的恒定束宽波束形成方法，在水平和垂直2 个方位均获得低旁瓣的恒定束宽波束图，有利于浅海波导噪声测量中缓解界面反射信号对测量结果的干扰。仿真结果表明，本文提出的二维恒定束宽波束形成方法能有效应用到双锥阵噪声测量系统中，相比单阵元和利用常规波束形成的测量方法，测量准确性更高。