CN112987004A - 一种浅海环境下基于水平阵列的水面水下目标分类方法 - Google Patents

Abstract

本发明属于水声工程、海洋工程技术等领域，具体涉及一种浅海环境下基于水平阵列的水面水下目标分类方法，该方法先利用岸基水平阵列实际测量获取目标辐射噪声时序信号，然后分别对各阵元接收信号进行常规波束形成和频率自适应的最优权重阵列波束形成处理，得到两种情况下的目标辐射噪声波束形成宽带时频域干涉谱图，最后基于两干涉谱图中的干涉条纹差异实现水上水下目标分类。本发明的有益结果是：该方法利用固定于海底的水平阵列，对接收的目标辐射声信号展开不同的波束形成处理完成目标的水面水下分辨，将浅海声场时频干涉结构特征融入目标识别过程，处理流程简单，实现运动声源快速分类，目标判别过程简明、结果直观。

Description

一种浅海环境下基于水平阵列的水面水下目标分类方法

技术领域

本发明属于水声工程、海洋工程技术等领域，具体涉及一种浅海环境下基于水平阵列的水面水下目标分类方法，该方法基于不同模态的声场干涉条纹差异性来分类水面和水下目标，适用于岸基水平阵列、浅海负跃层环境，目标深度在跃层上的水面目标和跃层下的水下目标。

背景技术

在水声领域，水面水下目标的分辨一直是一个亟待解决的问题。传统方法有利用匹配场处理和采用垂直阵列基于信道特性的信息处理方法等，该类方法由于利用海洋信道特性，对声环境参数要求高，声信息采集装备布放复杂等缺点。近年来发展的基于图像处理和机器学习的目标分类方法、基于目标辐射噪声特征分析和深度学习的目标分类方法等，将机器学习等先进的计算机技术引入水下目标信息处理，在目标样本充足时能较好的训练出分类目标的典型特征。但实际中，由于水下目标的数据样本少，该类方法估计的目标分类结果精确度一般不高。

在浅海环境中，目标的声传播特性可利用简正波模型来表征，因而不同深度声源辐射声场会由不同的模态组成。事实上，由于浅海环境通常为典型的负跃层环境(如我国东海、黄海的夏季环境)，导致海底岸基阵在浅海环境中接收到的水面目标声场主要由SRBR模态(海面海底反射模态)组成，而水下目标声场主要由nSRBR模态(非海面海底反射模态)中一种特殊的模态，即RBR模态(折射-海底反射模态)组成。采用波导不变量来描述SRBR模态和nSRBR模态时可以发现，SRBR模态波导不变量约为1，而nSRBR模态的波导不变量为负数或大于1。另一方面，不同取值的波导不变量对应在距离-频率谱图中展现的干涉条纹斜率方向存在显著差异。

此外，相对于传统方法采用的垂直阵列，水平阵列在工程上更加容易布放，且布放后阵形更加固定，这使得阵列波束形成的结果更加准确，有助于更方便地获得目标辐射噪声的时间-频率谱图。

发明内容

基于上述分析，如果基于水平阵列有效利用水面水下目标在模态上的差异性，就能够相对于传统方法更加方便、直观、准确地实现水面水下目标分类。为了拓展岸基水平阵列水面水下目标分辨能力，本发明提出一种浅海环境下基于水平阵列的水面水下目标分类方法，先利用岸基水平阵列实际测量获取目标辐射噪声时序信号，然后分别对各阵元接收信号进行常规波束形成和频率自适应的最优权重阵列波束形成处理，得到两种情况下的目标辐射噪声波束形成宽带时频域干涉谱图，最后基于两干涉谱图中的干涉条纹差异实现水上水下目标分类。该方法主要利用浅海环境中水面和水下目标声源辐射信号时，水平阵列接收到的目标声场模态在的距离-频率谱图中干涉条纹的角度差，实现水面水下目标的快速分类判别。

本发明采用的技术方案为，所述水平阵列布设于浅于100m深的海底、阵元个数N为大于30的奇数、水平阵列长度大于200m，海区水体声速剖面为负跃层，步骤如下：

步骤1：实际测量获取目标辐射声信号时间段内海区水体声速剖面最大值和海底声速值。

步骤2：实际测量获取N元水平阵列各阵元接收的目标辐射噪声时序信号声压向量。

步骤3：利用常规波束形成方法(详见文献：杜选民等著，声纳阵列信号处理技术，电子工业出版社，2018:52-53。)，得到被跟踪运动目标辐射噪声时间-频率谱图。

步骤4：利用频率自适应的最优权重阵列处理方法(详见文献：徐国军等，频率自适应最优权重阵列干涉条纹处理技术，声学学报，2017，42(3)：257-266。)，得到被跟踪运动目标辐射噪声时间-频率谱图。

步骤5：将步骤3和步骤4得到的被跟踪运动目标辐射噪声时间-频率谱图分别进行Radon变换，得到两幅Radon变换后的距离-角度图。

步骤6：将步骤5中得到的两幅Radon变换后的距离-角度图中各像素点分别沿距离维相加，得到两组角度估计向量，读出向量最大值对应的角度作为两个干涉条纹估计角度值。

步骤7：比较步骤6中两个干涉条纹估计角度值的角度差，当角度差小于5°时，判别目标为水面目标；当角度差大于5°时，判别为水下目标。

本发明的有益结果是：

该方法利用固定于海底的水平阵列，对接收的目标辐射声信号展开不同的波束形成处理完成目标的水面水下分辨，将浅海声场时频干涉结构特征融入目标识别过程，处理流程简单，实现运动声源快速分类，目标判别过程简明、结果直观。

附图说明

图1是本发明的实现流程图；

图2是Radon变换示意图；

图3是海区声速剖面及单水听器接收声场距离-频率谱。其中图3(a)为声速剖面，横坐标表示声速，单位为“m/s”，纵坐标为深度，单位为“m”；图3(b)和图3(c)分别为深度10m、60m声源的单水听器接收声场距离-频率谱，横坐标均为距离，单位为“m”，纵坐标为频率，单位为“Hz”。

图4为深度10m声源的常规波束形成处理结果。其中图4(a)为常规波束形成距离-频率谱，横坐标表示距离，单位为“m”，纵坐标为频率，单位为“Hz”；图4(b)为Radon变换结果，横坐标表示距离，单位为“m”，纵坐标为角度，单位为“°”；图4(c)为根据干涉条纹估计的角度结果，横坐标为角度，单位为“°”，纵坐标表示幅度，单位为“1”。

图5为深度10m声源的频率自适应最优权重波束形成处理结果。其中图5(a)为频率自适应波束形成距离-频率谱，横坐标表示距离，单位为“m”，纵坐标为频率，单位为“Hz”；图5(b)为Radon变换结果，横坐标表示距离，单位为“m”，纵坐标为角度，单位为“°”；图5(c)为根据干涉条纹估计的角度结果，横坐标为角度，单位为“°”，纵坐标表示幅度，单位为“1”。

图6为深度60m声源的常规波束形成处理结果。其中图6(a)为常规波束形成距离-频率谱，横坐标表示距离，单位为“m”，纵坐标为频率，单位为“Hz”；图6(b)为Radon变换结果，横坐标表示距离，单位为“m”，纵坐标为角度，单位为“°”；图6(c)为根据干涉条纹估计的角度结果，横坐标为角度，单位为“°”，纵坐标表示幅度，单位为“1”。

图7为深度60m声源的频率自适应最优权重波束形成处理结果。其中图7(a)为频率自适应波束形成距离-频率谱，横坐标表示距离，单位为“m”，纵坐标为频率，单位为“Hz”；图7(b)为Radon变换结果，横坐标表示距离，单位为“m”，纵坐标为角度，单位为“°”；图7(c)为根据干涉条纹估计的角度结果，横坐标为角度，单位为“°”，纵坐标表示幅度，单位为“1”。

具体实施方式

图1是本发明的实现流程图。具体实施方式包括以下步骤：

步骤1：实际测量获取目标辐射声信号时间段内海区水体声速剖面最大值和海底声速值。

所述海区水体声速剖面最大值为海区水体各深度上声速最大值，所述海底声速值为海底沉积层声速值。

步骤2：实际测量获取N元水平阵列各阵元接收的目标辐射噪声时序信号声压向量。

步骤3：利用常规波束形成方法(详见文献：杜选民等著，声纳阵列信号处理技术，电子工业出版社，2018:52-53。)，得到被跟踪运动目标辐射噪声时间-频率谱图。

利用常规波束形成方法得到波束形成后的目标辐射噪声时间-频率谱图的过程为：

水平阵列波束响应函数可表示为：

其中

为波束响应，ω为角频率，θ_T,θ_S,r_t分别为导向方位角、目标方位角和目标与水平阵参考阵元(不失一般性地，在本具体实施方式中设为第一阵元)之间的距离，当目标运动时该距离随时间t变化，v(k_t)为各阵元导向矢量，W＝(w₁,w₂,…w_N)^T为各阵元加权权向量，T表示向量转置，P^H为各阵元接收的目标辐射噪声时序信号声压的复共轭转置向量。当导向方位与目标方位一致时，

即为在t时刻目标的输出波束在频率f时的结果，其中f＝ω/(2π)。

对阵列分析进行常规波束形成处理时，令各阵元权矢量为w₁＝w₂＝…＝w_N＝1。波束处理过程中，首先获取某时刻t的目标辐射噪声，计算目标方位与导向方位一致即θ_T＝θ_S时

的值，f为频率，最后遍历所有的时刻t和频率f，可得到波束形成后的目标辐射噪声时间-频率谱图。

步骤4：利用频率自适应的最优权重阵列处理方法(详见文献：徐国军等，频率自适应最优权重阵列干涉条纹处理技术，声学学报，2017，42(3)：257-266。)，得到被跟踪运动目标辐射噪声时间-频率谱图。

频率自适应的最优权重阵列处理方法中权矢量计算过程为：

频率自适应的最优权重阵列处理方法计算过程与步骤3基本一致，仅各阵元加权权值不同，即权向量W中各元素不均等于1，而是依赖于频率值和阵列信息。考虑到权重的对称性和阵元位置的对称性，仅需设计前一半阵元权矢量(后一半阵元权矢量与前一半对称)，令其为W⁽¹⁾＝(w₁,w₂,…w_(N+1)/2)^T，则有部分阵元权矢量计算公式如下：

其中c、A、a各向量元素个数均为(N+1)/2，各变量具体计算公式如下：

a＝(2 2 … 2 1)^T，

有l＝2π/d，d为阵元间距；

c_max、c_seafloor分别为海水水体最大声速值、海底声速值；d_i＝(i-(N+1)/2)d,i＝1,2,…,(N-1)/2；

所述步骤3和步骤4获得的时间-频率谱图中，频带范围取300-500Hz，每隔1Hz取一个频点；时间维总长取5分钟以上，每1秒为一个时间点。

通过该步骤，可以实现声场模态的滤波，输出声场中的SRBR模态。

步骤5：将步骤3和步骤4得到的被跟踪运动目标辐射噪声时间-频率谱图分别进行Radon变换，得到两幅Radon变换后的距离-角度图。

对时间-频率谱图中干涉条纹角度的获取利用Radon变换来实现。将待处理图形上的一条线(图2中的线P)变换成距离-角度图上的一个点p，线的积分值转变为点的幅值，变换后图2中横轴为角度、纵轴为距离。如图2，角度

(取值为0°～180°)为直线P相对于y轴的偏角，ρ为线P与变换后的距离-角度图中的距离一般取图形中点(即参考点)相对于直线的距离，显然图中多条斜率相同的条纹，经Radon变换后会在变换图中同一个角度上不同距离位置上呈现很多点，获得Radon变换后的距离-角度图。

步骤6：将步骤5中得到的两幅Radon变换后的距离-角度图中各像素点分别沿距离维相加，得到两组角度估计向量，读出向量最大值对应的角度作为两个干涉条纹估计角度值。

步骤7：比较步骤6中两个干涉条纹估计角度值的角度差，当角度差小于5°时，判别目标为水面目标；当角度差大于5°时，判别为水下目标。

应用实例：

设定海区环境为典型黄海夏季声速剖面为负跃层环境，0-20m为跃层上声学层、20-75m为负跃层、75-100m为跃层下弱负梯度，100-150m为沉积层。声速剖面如图3(a)所示。海区声速最大值和海底声速值分别为1539.6m/s和1572m/s。声源分别位于10m和60m，位布设于海底水平阵列的端射方向，开始时刻距离阵列最近阵元5km，远离阵列方向航行至8km，在200-300Hz频率范围内，单基元接收的声场距离-频度谱图如图3(b)、图3(c)所示。

图4为声源10米深时，海底水平阵列利用常规波束形成方法的处理结果，采用Radon变换对波束处理的声场结构干涉条纹处理得到条纹角度为80度。图5给出了声源10米深时，海底水平阵列利用频率自适应最优权重波束形成方法的处理结果，采用Radon变换对波束处理的声场结构干涉条纹处理得到条纹角度为79.5度。显然，两种波束形成方法处理后得到的条纹角度差小于5度，判别该目标为水面目标。

图6、图7分别为声源60米深时，海底水平阵列利用常规波束形成方法和频率自适应最优权重波束形成方位的处理结果。常规波束形成输出距离-频率谱图中估计的条纹角度为100度；而频率自适应最优权重波束形成输出距离-频率谱图中干涉条纹角度估计为79度。显然，两种波束形成方法处理后得到的条纹角度差大于5度，判别该目标为水下目标。

Claims (3)

1.一种浅海环境下基于水平阵列的水面水下目标分类方法，该方法针对水平阵列布设于浅于100m深的海底、阵元个数N为大于30的奇数、水平阵列长度大于200m，其特征在于，该方法包括以下步骤：

步骤1：实际测量获取目标辐射声信号时间段内海区水体声速剖面最大值和海底声速值；

步骤2：实际测量获取N元水平阵列各阵元接收的目标辐射噪声时序信号声压向量；

步骤3：利用常规波束形成方法，得到被跟踪运动目标辐射噪声时间-频率谱图；

步骤4：利用频率自适应的最优权重阵列处理方法，得到被跟踪运动目标辐射噪声时间-频率谱图；

步骤5：将步骤3和步骤4得到的被跟踪运动目标辐射噪声时间-频率谱图分别进行Radon变换，得到两幅Radon变换后的距离-角度图；

步骤6：将步骤5中得到的两幅Radon变换后的距离-角度图中各像素点分别沿距离维相加，得到两组角度估计向量，读出向量最大值对应的角度作为两个干涉条纹估计角度值；

步骤7：比较步骤6中两个干涉条纹估计角度值的角度差，当角度差小于5°时，判别目标为水面目标；当角度差大于5°时，判别为水下目标。

2.一种根据权利要求1所述浅海环境下基于水平阵列的水面水下目标分类方法，其特征在于：所述海区水体声速剖面最大值为海区水体各深度上声速最大值，所述海底声速值为海底沉积层声速值。

3.一种根据权利要求1所述浅海环境下基于水平阵列的水面水下目标分类方法，其特征在于：利用常规波束形成方法得到波束形成后的目标辐射噪声时间-频率谱图的过程为：

水平阵列波束响应函数可表示为：

其中

为波束响应，ω为角频率，θ_T,θ_S,r_t分别为导向方位角、目标方位角和目标与水平阵参考阵元之间的距离，当目标运动时该距离随时间t变化，v(k_t)为各阵元导向矢量，W＝(w₁,w₂,…w_N)^T为各阵元加权权向量，T表示向量转置，P^H为各阵元接收的目标辐射噪声时序信号声压的复共轭转置向量；当导向方位与目标方位一致时，

即为在t时刻目标的输出波束在频率f时的结果，其中f＝ω/(2π)；

对阵列分析进行常规波束形成处理时，令各阵元权矢量为w₁＝w₂＝…＝w_N＝1；波束处理过程中，首先获取某时刻t的目标辐射噪声，计算目标方位与导向方位一致即θ_T＝θ_S时

的值，f为频率，最后遍历所有的时刻t和频率f，可得到波束形成后的目标辐射噪声时间-频率谱图。

Images