一种深海环境下单水听器声源被动定位方法技术

本发明专利技术涉及一种深海环境下单水听器声源被动定位方法，利用单个水听器接收的宽带信号，不需要多个水听器组成大孔径的阵列，不需要水听器间的信号同步，同时在被动接收方式下不会暴露自身位置即可实现在深海环境下的声源定位。本发明专利技术将无法区分的六个相对时间延迟按照从小到大的顺序排列，构成相对时延向量，用于声源定位。本发明专利技术将测量时间延迟向量与预测时间延迟向量相减得到的误差向量求1范数并求倒数对声源定位。本方法用于夏季南海2680米海深环境，声源为带宽2.9kHz‑3.1kHz的线性调频信号，声源距离13.8公里，深度49米，完成对声源进行定位，距离相对误差2.9％，深度相对误差22.5％。

【技术实现步骤摘要】
一种深海环境下单水听器声源被动定位方法
本专利技术属于水声探测、声呐技术等领域，涉及水中的声源定位方法，尤其涉及一种深海环境下单水听器声源被动定位方法。
技术介绍
水下声源定位方法通常是通过提取单个水听器(传感器)或多个按照一定方式排列的水听器阵列接收到信号所包含的信息来重构声源位置信息。比较常见的主要是匹配场处理方法，但该方法需要较大的孔径和阵元数来采集足够的信息，才能对声源进行定位，并且这类方法计算量大、受环境失配影响较大，不适用于深海环境下的声源定位。在深海环境下，声波的传播路径主要有两种，分别是汇聚区模式和海底弹射模式。利用海底弹射模式的信号时延特征，重构声源位置信息，可以减小计算量，并减小对水听器数量的需求。现有的深海声源定位方法多数需要水听器组成阵列采集足够的信息，或需要工作在主动模式下。公布号为CN107202988A的专利涉及一种利用航船不停移动接收到声学释放器应答的距离信息进行深海潜标定位的方法，该方法简便易行，但是应用范围小，仅限于对短距离潜标系统定位；公布号为CN104793212A的专利涉及一种通过水听器阵列主动探测目标位置信息的方法，主动探测方式容易暴露自身位置；公布号为CN106707240A的专利涉及一种提取水听器阵列接收到信号时延信息对声源深度进行估计的方法，该方法计算量小，但是需要同步的较多数量的水听器组成阵列，并且只能得到声源的深度信息。
技术实现思路
要解决的技术问题为了避免现有技术的不足之处，本专利技术提出一种深海环境下单水听器声源被动定位方法，解决现有深海环境下声源定位方法需要水听器数量多、需要主动发声等问题。技术方案一种深海环境下单水听器声源被动定位方法，其特征在于步骤如下：步骤1：将单个水听器采用潜标或者浮标的方式布放在海中，水听器接收到海底弹射模式下的来自四条路径的宽带信号为：海底反射信号B、海面-海底反射信号SB、海底-海面反射信号BS、海面-海底-海面反射信号SBS；不同路径的信号的传播时间记为TB、TSB、TBS、TSBS；通过水听器的深度传感器获取水听器的深度z；所述深度z指水听器与海平面的距离；步骤2：对接收的宽带信号进行自相关处理，得到四条途径到达信号之间的六个相对时间延迟；自相关处理：其中：s(t)为接收信号，T是信号长度，(·)*表示信号的复共轭，τ是积分变量；对信号自相关输出取包络，并归一化C(τ)＝|R(τ)+j*Hilbert(R(τ))|其中：j等于Hilbert(·)表示信号的希尔伯特变换，|·|表示绝对值；得到四条路径在自相关输出的正半轴形成六个峰值，表示六个相对时间延迟；当TBS＞TSB时，这六个峰值对应的时延分别记为τSB-B，τBS-B，τSBS-B，τSBS-SB，τSBS-BS，τBS-SB，当TSB＞TBS时，τBS-SB是负数，用τSB-BS替代；所述的相对时延表示为：步骤3：将六个时间延迟从小到大排序，构成一个时间延迟测量向量：τm(r,d)＝[sort(τSB-B,τBS-B,τSBS-B,τSBS-BS,τSBS-SB,|τBS-SB|)]其中：r是声源与接收水听器的水平距离，d是声源深度，sort(·)表示将时延从小到大排列，|·|表示绝对值；由于未知TBS和TSB的大小，采用绝对值表示；步骤4：利用代价函数求解声源位置其中：||·||1表示向量的1范数，是某个搜索网格点的深度，是某个搜索网格点距离接收水听器的水平距离；所述的搜索网格点为：对感兴趣的区域的深度，以深度维度网格划分步长进行深度网格划分，以距离维度网格划分步长进行网格划分；所述τp是根据环境信息，假设声源在某个搜索网格点时，用声场计算模型BELLHOP进行计算，得到的时间延迟预测向量：步骤5：找到代价函数最大值点以及最大值点对应的网格点的深度和距离接收水听器的水平距离，该深度和水平距离接近声源的深度和距离接收水听器的水平距离。所述感兴趣的区域为目标可能存在的区域。所述目标可能存在的感兴趣的距离范围为0-20公里，深度范围为0-500米。所述深度维度网格划分步长1米。所述距离维度网格划分步长为0.01公里。有益效果本专利技术提出的一种深海环境下单水听器声源被动定位方法，利用单个水听器接收的宽带信号，不需要多个水听器组成大孔径的阵列，不需要水听器间的信号同步，同时在被动接收方式下不会暴露自身位置即可实现在深海环境下的声源定位。本专利技术将无法区分的六个相对时间延迟按照从小到大的顺序排列，构成相对时延向量，用于声源定位。本专利技术将测量时间延迟向量与预测时间延迟向量相减得到的误差向量求1范数并求倒数对声源定位。本方法用于夏季南海2680米海深环境，声源为带宽2.9kHz-3.1kHz的线性调频信号，声源距离13.8公里，深度49米，完成对声源进行定位，距离相对误差2.9％，深度相对误差22.5％。附图说明图1为深海环境下声速不随深度变化时四条理想声线传播路径；图2为实验测量声速剖面，海底地形以及实验配置图；图3为实验数据、数据自相关与自相关部分放大图；图4为感兴趣范围内声源定位结果。具体实施方式现结合实施例、附图对本专利技术作进一步描述：步骤1：将单个水听器采用潜标或者浮标的方式布放在海中，通过深度传感器获取水听器的深度信息。所述的深度z指水听器与海平面的距离，海洋深度为H，在深海环境下，声源深度为d，声源与水听器的水平距离为r，当d＜H/2时，水听器能够更好的接收到海底弹射模式下的来自四条路径的信号。来自四条路径的信号分别是：海底反射信号(B)、海面-海底反射信号(SB)、海底-海面反射信号(BS)、海面-海底-海面反射信号(SBS)。不同路径的信号具有不同的传播时间，分别记为TB、TSB、TBS、TSBS。步骤2：对于接收到的宽带信号，通过对信号进行自相关处理可以在未知信号先验信息的情况下得到四条途径到达信号之间的六个相对时间延迟。所述的接收信号设为s(t)，信号的自相关可以表示为其中T是信号长度，(·)*表示信号的复共轭，τ是积分变量。为了更精确的确定信号时延，对信号自相关输出取包络，并归一化C(τ)＝|R(τ)+j*Hilbert(R(τ))|其中j等于Hilbert(·)表示信号的希尔伯特变换，|·|表示绝对值。四条路径在自相关输出的正半轴形成六个峰值，当TBS＞TSB时，这六个峰值对应的时延分别记为τSB-B，τBS-B，τSBS-B，τSBS-SB，τSBS-BS，τBS-SB，当TSB＞TBS时，τBS-SB是负数，用τSB-BS替代。所述的相对时延表示为：步骤3：时间延迟与声源深度和距离有关，当声源深度为d，水平距离为r时，将六个时间延迟从大到小排序，构成一个时间延迟测量向量。所述的六个相对时延在未知声源位置时无法确定相对大小，即不能将六个峰值一一对应到六个相对时延，为了利用相对时延信息重构声源位置，将由自相关输出峰值得到的六个相对时延数值按照从小到大排序，构成时延测量向量。τm(r,d)＝[sort(τSB-B,τBS-B,τSBS-B,τSBS-BS,τSBS-SB,|τBS-SB|)]其中sort(·)表示将时延从小到大排列，|·|表示绝对值。由于未知TBS和TSB的大小，采用绝对值表示。步骤4：根据环境信息用声场计算模型BELLHOP进本文档来自技高网...

【技术保护点】
一种深海环境下单水听器声源被动定位方法，其特征在于步骤如下：步骤1：将单个水听器采用潜标或者浮标的方式布放在海中，水听器接收到海底弹射模式下的来自四条路径的宽带信号为：海底反射信号B、海面‑海底反射信号SB、海底‑海面反射信号BS、海面‑海底‑海面反射信号SBS；不同路径的信号的传播时间记为TB、TSB、TBS、TSBS；通过水听器的深度传感器获取水听器的深度z；所述深度z指水听器与海平面的距离；步骤2：对接收的宽带信号进行自相关处理，得到四条途径到达信号之间的六个相对时间延迟；自相关处理：

【技术特征摘要】
1.一种深海环境下单水听器声源被动定位方法，其特征在于步骤如下：步骤1：将单个水听器采用潜标或者浮标的方式布放在海中，水听器接收到海底弹射模式下的来自四条路径的宽带信号为：海底反射信号B、海面-海底反射信号SB、海底-海面反射信号BS、海面-海底-海面反射信号SBS；不同路径的信号的传播时间记为TB、TSB、TBS、TSBS；通过水听器的深度传感器获取水听器的深度z；所述深度z指水听器与海平面的距离；步骤2：对接收的宽带信号进行自相关处理，得到四条途径到达信号之间的六个相对时间延迟；自相关处理：其中：s(t)为接收信号，T是信号长度，(·)*表示信号的复共轭，τ是积分变量；对信号自相关输出取包络，并归一化C(τ)＝|R(τ)+j*Hilbert(R(τ))|其中：j等于Hilbert(·)表示信号的希尔伯特变换，|·|表示绝对值；得到四条路径在自相关输出的正半轴形成六个峰值，表示六个相对时间延迟；当TBS＞TSB时，这六个峰值对应的时延分别记为τSB-B，τBS-B，τSBS-B，τSBS-SB，τSBS-BS，τBS-SB，当TSB＞TBS时，τBS-SB是负数，用τSB-BS替代；所述的相对时延表示为：步骤3：将六个时间延迟从小到大排序，构成一个时间延迟测量向量：τm(r,d)＝[sort(τSB-B,τBS-B,τSBS-...

【专利技术属性】
技术研发人员：杨益新，刘文旭，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西,61