一种小斜视多接收阵合成孔径声呐距离多普勒成像算法

吕金华，唐扶光，赵 煦，吴浩然

(1. 武汉船舶职业技术学院电气与电子工程学院，湖北 武汉 430050；

2. 武汉轻工大学电气与电子工程学院，湖北 武汉 430023；

3. 海军工程大学电子工程学院，湖北 武汉 430033)

0 引 言

合成孔径声呐（synthetic aperture sonar，SAS）利用小尺寸基阵沿运动方向做匀速直线运动来合成大的孔径基阵，获得沿运动方向（横向）的高分辨率[1]。海流和载体平台的非对称等因素，可能导致合成孔径声呐出现斜视。另外，由于多子阵合成孔径声呐的斜视角会导致时延误差以及声呐多普勒效应比雷达显著，即使斜视角很小，也会导致合成孔径声呐图像散焦。

斜视合成孔径声呐与斜视合成孔径雷达信号模型最大的不同有两点：一是由于水中声速低，造成方位频率不模糊和距离不模糊存在矛盾。为了解决这个矛盾。合成孔径声呐的接收阵一般采用多子阵配置[2]。在现有的斜视合成孔径成像算法中，仅有单接收阵斜视合成孔径成像算法，包括单站斜视合成孔径雷达成像算法和双基单接收阵斜视合成孔径成像算法。其中单站斜视合成孔径雷达成像算法有：RDA[3]，CSA[4]，omega-K 算法[5]以及它们的修正算法[6]和斜视子孔径算法[7]，双基单接收阵斜视合成孔径成像算法有：双基RDA[8]、双基CSA[9]和双基omega-K 算法[10]。二是由于声呐平台航速和水下声速可比拟，导致在一般合成孔径雷达中常用的“停走停”假设在合成孔径声呐上不适用[11-12]。

因此，斜视单基站成像算法不能直接用于斜视多接收阵 SAS。本文提出一种小斜视多接收阵 SAS 距离多普勒成像算法，并通过仿真实验证明了算法的有效性和正确性。

1 斜视多接收阵SAS 模型

1.1 精确距离史模型

图 1 斜视多接收阵SAS 几何模型Fig. 1 SAS geometric model of slanting multisensor array

由于在斜视情况下发射阵、接收阵与目标三者不位于一个平面内，需要在三维直角坐标系内建立多子阵斜视合成孔径声呐模型。定义声呐运动方向为轴方向，垂直向下的方向为轴方向，用右手准则确定 y轴方向。设以发射阵位于点的时刻为慢变时间t的起点， 此时波束中心射线通过点目标rsinθr,h) i)tE(vt,0,0)，第个接收子阵位于点。在 时刻，发射阵移动到 点发射信号，那么此时发射阵与点目标的距离为：

经过t*i(t;r) 后， 第i 个接收子阵移动到了(())点，才接收到点目标的回波信号，那么此时第个接收子阵与点目标的距离为：

其中：

由于接收和发射的间隔时间 等于声波通过点目标P 的精确距离史 所需要的时间，因此也叫点目标的精确时延史，还可以表示为：

由式（5）和式（6）解得

其中：

1.2 修正距离史

将式（11）代入式（5），得到近似距离史为：

斜视SAS 的中接收阵排列方向与运动方向存在一个固定的角度偏差，使得斜视多子阵合成孔径声呐与正侧视多子阵合成孔径声呐单基近似过程不一样。如图2 所示，首先将每个接收阵元投影至方位轴上，然后取接收阵的投影位置和发射阵位置的中点作为单基采样点。由于接收子阵在方位轴上的投影长度为，为了满足方位均匀采样，脉冲重复频率（PRF）需要调整为。

图 2 斜视时发射阵位置和接收阵的空间位置以及方位均匀采样时相位中心的位置Fig. 2 Spatial position of transmitting array and receiving array as well as position of phase center during azimuth uniform sampling in the case of squint

由以上分析可知，点目标P 的斜视修正距离史为：

其中：Ri(t;r) 表示修正量，写为

ΔR(r;di)表示修正量，写为

其中：第1 项为平方项，包括非停走停模式引入采样点在方位上的移动距离和?方位向上投影长度；第2 项表示斜视条件下接收子阵引入的偏移量。

通过比较式（13）和精确距离史式（5），得到波长归一的距离史误差为：

在图3 中，距离史误差都是在波束边沿处最大，随着距离增大而增大，随着斜视角增大而增大。为了在不同斜视角下定量地比较距离史误差，分别对图3进行测量，得到距离史误差的最大变化量分别为0.0105λ 0.0248λ 0.0537λ 和0.1128λ从测量结果看，斜视角在0°～4.4°范围内的距离史误差小于 ，满足成像需要。当斜视角增加到6.8°时，距离史误差超过 的1 倍。因此式（13）能够满足小斜视条件下的成像要求。

1.3 信号模型

2 小斜视多接收阵SAS 距离多普勒成像算法

本文提出的小斜视多接收阵SAS 距离多普勒成像算法，实现过程如图4 所示，包括多接收阵信号处理和斜视单基站RDA 两大部分。

图 3 不同斜视角下的距离史误差Fig. 3 Distance history errors at different oblique angles

2.1 斜视多接收阵信号单基处理

从式（13）可知，相比斜视单基合成孔径信号，斜视多接收阵SAS 信号包含带来的相位和时延，以及收发分置项。这3 项的处理对应斜视多接收阵信号的单基过程。

由式（18）和式（13），可以得到ΔR(r;di)对应的相位补偿因子 为：

由于ΔR(r;di)是弱距离依赖的，对应的时延可以用参考距离 （选为测绘带中心）上的时延代替。在距离向频域通过相位相乘的方式完成时延补偿，对应的时延补偿因子为：

图 4 算法流程图Fig. 4 Algorithm flow chart

假设完成时延补偿后的信号变换至二维时域，得到

对于多接收阵SAS 来说，单个接收阵信号是方位欠采样的。为了获得满足奈奎斯特采样定理的方位信号，需要对多接收阵信号进行方位重构。由于所有的接收阵几乎同时接收同一个点目标的回波信号，必须将不同的接收阵信号沿着方位时间轴错开不同的时间间隔，才能等效为一个接收阵沿着方位时间轴在不同的位置接收回波信号，从而得到方位向采样点数增加的倍数与阵元个数相等，采样频率提高倍数与阵元个数相等的方位重构信号。通过对式（21）所示的多接收阵信号按照逐阵元逐脉冲的排列，得到方位重构的信号为：

其中：

此时斜视多接收阵SAS 信号式（22）能够等效为斜视单基合成孔径信号。

2.2 斜视单基RDA

为了推导斜视单基RDA，利用驻定相位原理将式（22）变换至二维频域，得到

将式（25）对 进行泰勒级数展开fr，保留至f2r，得到

式（26）中第1 项表示距离向匹配滤波项，包含发射信号调制项和二次距离压缩项，其中二次距离压缩项是弱距离依赖的，一般可以用参考距离上对应的二次距离压缩调频斜率代替整个场景的调频斜率；第2 项表示距离徙动项，是距离徙动的来源；第3 项表示方位调制项，是方位匹配滤波的来源；第4 项是“非停走停”模式引起的方位向线性走动量。

从式（26）的第1 项可得实现距离向脉压和二次距离压缩的相位函数为：

从式（26）的第2 项可得距离徙动校正量为：

完成距离徙动校正后，从式（26）可以看出方位向脉压和“非停走停”引起的线性走动量可以合并，用一个相位函数实现。该相位函数写为

最后，进行方位向逆傅里叶变换，即可得到成像结果。

3 实验验证

为了验证本文提出的小斜视多接收阵RDA 的有效性，进行仿真实验。多接收阵合成孔径声呐的系统参数如表1 所示，点目标的时延由精确时延公式（7）给出，仿真结果如图5 所示。

表 1 系统仿真参数Tab. 1 System simulation parameters

当斜视角为2°时，比较图5（a）和图5（d）发现与斜视角为0°的成像结果很相似，都对场景中的目标进行了很好的聚焦。不同之处是图5（d）出现了由于斜视造成的图像几何形变。为更详细比较图5 中点目标的成像质量，将点目标从图5（a）和图5（b）中提取出来，并画出距离向和方位向剖面，分别如图5（b）、图5（c）、图5（e）和图5（f）所示。然后分别对目标的方位向和距离向IRW（Impulse Response Width）及PSLR（Peak Sidelobe Ratio）和ISLR（Integration Sidelobe Ratio）进行测量，结果如表2 所示。从图5（b）、图5（c）、图5（e）和图5（f）和表2 的对比结果可知，本文提出的小斜视多接收阵RDA 消除了斜视对成像结果的影响，得到了和正侧视相同的成像结果。

4 结 语

为解决斜视多子阵合成孔径声呐成像问题，本文首先建立斜视多子阵合成孔径声呐精确的几何模型，给出了精确的时延史。为推导成像算法，对该时延史进行2 次近似，分析了近似误差，结果表明在窄波束小斜视角情形下，总的近似误差满足成像要求。然后给出了相应的信号模型。在算法推导部分，借鉴经典RDA 算法。提出小斜视角多子阵合成孔径声呐RDA算法。最后通过计算机仿真实验证明了算法的有效性和正确性。

图 5 不同斜视角时点目标的仿真结果Fig. 5 Simulation results of point targets at different oblique angles

表 2 图像质量比较Tab. 2 Comparison of image quality