基于单固定水听器和单移动声源的海洋地声参数反演方法技术

本发明专利技术公开了一种基于单固定水听器和单移动声源的海洋地声参数反演方法。所述反演方法包括：利用单固定水听器采集由移动合作声源产生的水声信号；利用声传播的互易性，利用声源移动带来的空间增益，构建虚拟合成阵列接收信号；利用KRAKEN前向模型，在多普勒情况下，计算拷贝声场；对接收信号进行傅立叶变换，在指定频率上，使用频率相干的BARTLETT相关器作为目标函数，计算拷贝声场和测量声场的匹配程度；利用差分进化算法，获取使得匹配程度达到最大值的最优地声参数估计。本发明专利技术提供的单固定水听器和单移动声源的海洋地声参数反演方法，能够准确估计海洋地声参数，并提供了资源需求最小化和设备部署便捷化的解决方案。

【技术实现步骤摘要】
基于单固定水听器和单移动声源的海洋地声参数反演方法
本专利技术涉及传感器(水听器)采集和海洋探测领域，尤其涉及一种基于单固定水听器和单移动声源的海洋地声参数反演方法。技术背景海洋地质方面，海底的沉积层一般被建模为流体，在沉积层底部是硬质的半无限空间，如岩石等，这被称为海洋基底。同样地，一些海洋环境底部则是无沉积层覆盖的基底，介质需要建模为弹性介质，同时支持压缩波和剪切波的传播；海底结构方面，在浅海环境中，声线在水体和海底交界处发生反射，声线传播的多径效应明显。因此，海底的结构和地质、地声特性直接影响声的传播。这些环境参数对于声波传播损失(TransmissionLoss,TL)的预测和声纳性能的评估起着非常重要的影响。地声参数包括沉积层顶部、底部声速、密度、衰减，基底声速、密度、衰减等等，其获取方法主要为直接测量法和间接反演算法。由于直接测量法难度大、开销高，且无法获取全面的地声参数，往往被限制在沉积层的表面几米，间接地通过测量声信号，反演出地声参数，逐渐被学界与业界关注。通过地声参数反演，可以快速地反演出一片海域的地声模型参数，能够最大程度上避免在海底直接测量的困难和危险。目前常用的地声参数反演方法主要为如下几种：(1)到达时间分析方法；(2)模态色散技术；(3)匹配场反演。为获取地声参数反演算法必要的测量声场信息，通常利用垂直阵，水平阵，拖曳阵，分布式水下传感网络等等来获取测量声场。声源通常根据反演区域的规模，和不同频率信号在沉积层的穿透能力，选择低频至中频的信号(反演规模由小及大)。通常情况下，水下设备昂贵，部署难度大，在实际过程中，需要考虑如下几个问题：(1)部署成本，包括使用水听器数量，阵列规模，声源数量等等；(2)部署难度，包括是否锚底等等；(3)信号采集方式，包括同步采集等等。
技术实现思路
针对现有技术存在的缺陷，本专利技术提供了一种基于单固定水听器和单移动声源的海洋地声参数反演方法。该方法利用单水听器和单声源，极大程度上降低了设备成本；同时，通过船只拖曳，或使用水下自助航行器(AUV)搭载声源的方式，进一步降低了阵列拖曳或锚节点布放等部署成本。在部署成本和部署难度都维持在较低水平的同时，获取了阵列接收信号相同的效果。同时反演的优化问题中，采用了频率相干的目标函数，摆脱了时间同步的问题，使得不同时刻移动声源的空间增益得以充分利用。本专利技术的目的是通过以下技术方案来实现的：一种基于单固定水听器和单移动声源的海洋地声参数反演方法，其包括以下步骤：步骤1：利用单固定水听器采集由移动合作声源产生的水声信号；步骤2：利用声传播的互易性，利用声源移动带来的空间增益，构建虚拟合成阵列接收信号；步骤3：利用KRAKEN前向模型，结合移动声源的多普勒效应，计算多普勒效应下的拷贝声场；步骤4：匹配场反演：对接收信号进行傅立叶变换，在指定若干频率上，使用频率相干的BARTLETT相关器作为目标函数，计算拷贝声场和测量声场的匹配程度，进而利用差分进化算法，获取使得匹配程度达到最大的最优地声参数估计。基于上述技术方案，其中的各步骤还可以采用如下具体方式实现。作为优选，所述步骤1具体实现方法为：在待探测水域内，使用船只或水下自助航行器搭载中频声源按规划航迹航行，并以相同时间轴记录航迹信息与发射波形；该发射波形作为反演所用的移动合作声源；单固定水听器用以采集由合作声源产生的水声信号，采集时间覆盖移动声源的整个航行时间。进一步的，所述步骤1中，所述单固定水听器悬挂固定于浮标、船只或水下自助航行器上，或锚定于海底实现固定。作为优选，所述步骤2具体实现方法为：利用所述的航迹信息，选取M个同一直线上不同时刻的声源位置，获取该M个时刻后一定时间内由固定水听器采集的大小为M×N的信号矩阵y，作为接收声场向量，其中N为数据长度；根据声传播的互易性，交换接收端和发射端的位置，声线沿原传播的路径，反向传播；航迹上任意时刻，声源所处位置发射信号，由水听器接收，均视为水听器所处位置发出的声信号由声源所处位置的传感器接收；由此在所述M个时刻位于不同位置的移动声源，构成M个时刻移动声源对应位置上的M个传感器接收到来自于固定水听器位置的虚拟声源所发射的信号；M个位置的接收信号，共同组成虚拟合成阵列接收信号，即一个M元线阵的测量水声信号y。作为优选，所述步骤3具体包括以下子步骤：步骤3-1：首先在静止情况下，利用KRAKEN前向模型，计算接收频率、发射频率下的模深度函数以及波数、群速度，得到：{Ψn(zr；ωr),Ψn(zs；ωs),krn,urn}其中Ψn(zr；ωr)为接收频率下的模深度函数，Ψn(zs；ωs)为发射频率下的模深度函数，krn为波数，urn为群速度；zr为接收节点与水面的距离，zs为发送节点与水面的距离，n代表模的阶数；步骤3-2：利用步骤3-1中得到的krn和urn，结合运动声源的径向速度vs计算多普勒情况下的修正波数kn：步骤3-3：基于修正波数kn计算接收频率ωr的对应多普勒偏移发射频率ωs：其中，为修正波数kn下计算得到的多普勒偏移发射频率，vr为接收节点的速度；步骤3-4：通过发射信号的频谱在指定频率ωs处的幅度和相位信息结合运动声源的位置信息r＝r0+vst，计算声压场ψ(r,z,ωr)，其中r0为该时刻运动声源的瞬时位置，t为运动声源的运动时长，声压场计算公式如下：其中为汉克尔函数的近似，i为虚数，Ψn(z；ωs)为接收节点的模深度函数的近似；最终得到的ψ(r,z,ωr)即为多普勒修正下的拷贝声压场。步骤4：匹配场反演：对接收信号进行傅立叶变换，在指定若干频率上，使用频率相干的BARTLETT相关器作为目标函数，计算拷贝声场和测量声场的匹配程度，进而利用差分进化算法，获取使得匹配程度达到最大的包括沉积层顶部、底部声速、密度、衰减，基底声速、密度、衰减等最优地声参数估计。作为优选，所述步骤4具体实现方法为：对接收声场数据y的每一路信号做傅里叶变换，选取指定频率构成频域上的1×M维接收声场向量d＝[d1,d2,…,dM]，记发射信号指定频点的幅度和相位信息组成的对角矩阵为S，d(m)为待反演参数向量m通过前向模型获得的拷贝场向量，则频率相干的目标函数为：其中dn(m)为拷贝场向量d(m)的第n维；基于所述目标函数，利用差分进化算法，获取使得拷贝声场和测量声场匹配程度达到最大的最优地声参数估计，所述地声参数包括沉积层顶部、底部声速、密度、衰减，以及基底声速、密度、衰减。本专利技术的有益效果是：(1)本专利技术利用单水听器和单声源的地声参数反演方法，极大程度上降低了探测地声参数过程中的设备成本、部署成本(2)本专利技术方法利用声传播的互易性，利用移动声源的空间信息，获取了与使用阵列接收的声场数据进行反演相同的效果。(3本文档来自技高网...

【技术保护点】
1.一种基于单固定水听器和单移动声源的海洋地声参数反演方法，其特征在于，包括以下步骤：/n步骤1：利用单固定水听器采集由移动合作声源产生的水声信号；/n步骤2：利用声传播的互易性，利用声源移动带来的空间增益，构建虚拟合成阵列接收信号；/n步骤3：利用KRAKEN前向模型，结合移动声源的多普勒效应，计算多普勒效应下的拷贝声场；/n步骤4：匹配场反演：对接收信号进行傅立叶变换，在指定若干频率上，使用频率相干的BARTLETT相关器作为目标函数，计算拷贝声场和测量声场的匹配程度，进而利用差分进化算法，获取使得匹配程度达到最大的最优地声参数估计。/n

【技术特征摘要】
1.一种基于单固定水听器和单移动声源的海洋地声参数反演方法，其特征在于，包括以下步骤：
步骤1：利用单固定水听器采集由移动合作声源产生的水声信号；
步骤2：利用声传播的互易性，利用声源移动带来的空间增益，构建虚拟合成阵列接收信号；
步骤3：利用KRAKEN前向模型，结合移动声源的多普勒效应，计算多普勒效应下的拷贝声场；
步骤4：匹配场反演：对接收信号进行傅立叶变换，在指定若干频率上，使用频率相干的BARTLETT相关器作为目标函数，计算拷贝声场和测量声场的匹配程度，进而利用差分进化算法，获取使得匹配程度达到最大的最优地声参数估计。


2.根据权利要求1所述的一种基于单固定水听器和单移动声源的海洋地声参数反演方法，其特征在于，所述步骤1具体实现方法为：
在待探测水域内，使用船只或水下自助航行器搭载中频声源按规划航迹航行，并以相同时间轴记录航迹信息与发射波形；该发射波形作为反演所用的移动合作声源；单固定水听器用以采集由合作声源产生的水声信号，采集时间覆盖移动声源的整个航行时间。


3.根据权利要求2所述的一种基于单固定水听器和单移动声源的海洋地声参数反演方法，其特征在于，所述步骤1中，所述单固定水听器悬挂固定于浮标、船只或水下自助航行器上，或锚定于海底实现固定。


4.根据权利要求2所述的一种基于单固定水听器和单移动声源的海洋地声参数反演方法，其特征在于，所述步骤2具体实现方法为：
利用所述的航迹信息，选取M个同一直线上不同时刻的声源位置，获取该M个时刻后一定时间内由固定水听器采集的大小为M×N的信号矩阵y，作为接收声场向量，其中N为数据长度；根据声传播的互易性，交换接收端和发射端的位置，声线沿原传播的路径，反向传播；航迹上任意时刻，声源所处位置发射信号，由水听器接收，均视为水听器所处位置发出的声信号由声源所处位置的传感器接收；由此在所述M个时刻位于不同位置的移动声源，构成M个时刻移动声源对应位置上的M个传感器接收到来自于固定水听器位置的虚拟声源所发射的信号；M个位置的接收信号，共同组成虚拟合成阵列接收信号，即一个M元线阵的测量水声信号y。


5.根据权利要求4所述的一种基于单固定水听器和单移动声源的海洋地声参数反演方法，其特征在于，所述步骤3具体包括以下子步骤：
步骤3-1：首先...

【专利技术属性】
技术研发人员：潘翔，沈一宁，姜景宁，章杰，楼晨露，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江;33