一种基于声阵的海洋声传播损失实时测量方法技术

本发明专利技术提出一种基于声阵的海洋声传播损失测量方法，包括利用水声接收阵系统实时实现传播损失测量，并且在阵元域的信噪比小于6dB情况下，仍可通过声接收阵的波束域输出结果，实时实现不同频段水下声传播损失的计算。与传统的基于浮标的静态传播损失测量技术相比，基于声阵的海洋声传播损失实时测量技术可以实时进行阵元域和波束域的多频段的水下声传播损失测量。并且在阵元域信噪比低于6dB而无法进行传播损失测量的情况下，利用声基阵的阵增益，通过波束域可以实现多频段的水下声传播损失测量。

【技术实现步骤摘要】
一种基于声阵的海洋声传播损失实时测量方法
本专利技术属于水声信号处理
，具体涉及一种海洋声传播损失的实时测量方法。
技术介绍
随着水下目标隐身技术的发展，作为探测设备的主/被动声纳，其低频、大孔径、大功率特征愈加明显。主/被动拖曳线列阵声纳已成为最主要的声纳探测设备，它具有远离本艇干扰、深度可变、低频大孔径、主动声源级高等优点，已成为探测安静型潜艇的最重要设备。为充分发挥声纳设备，尤其是主被动拖曳线列阵声纳在各种海区的使用效能，获取工作海区的海洋声学环境十分必要。现代海洋斗争的特点，已由冷战时期的深海对抗转为既兼顾深海同时又兼顾近岸复杂浅水域，对声纳的探测及目标参数估计提出了更高的性能要求。因此，对复杂多变海洋声学环境的充分了解，直接关系到声纳的性能预报、最佳效能发挥、设计参数优化。水声信号的传播衰减直接关系声纳作用距离，传播衰减与工作频段、声速剖面、介质/界面吸收/反射等因素密切相关，测量并获取传播衰减曲线，可以为声纳性能的预报和设计提供依据，将实测曲线与模型估计曲线同步比对，可提高声纳性能预报的准确度。
技术实现思路
本专利技术所要解决的技术问题是：利用水声接收阵系统实时实现传播损失测量，在阵元域的信噪比小于6dB情况下，利用声阵的阵增益，通过声接收阵的波束域输出结果，实时实现不同频段水下声传播损失的计算。为解决以上技术问题，本专利技术是通过以下技术方案实现的：一种基于声阵的海洋声传播损失实时测量方法，其特征在于：采用水平拖曳线列阵实现，所述水平拖曳线列阵系统由16段声学段组成，各声学段由其内置的水听器阵列、前置与传输部分一起，完成声信号的声电转换、前置预处理、信号采集和传输功能；阵元数为184元，为实现波束域不同频段传播损失计算，采取布阵方式是：前64个声通道，阵元间距1.5m，中间24个声通道，阵元间距24m，最后96个声通道，阵元间距0.4m；声基阵总长758.4m；用于单阵元传播损失实时测量，具体包括以下步骤：步骤一：对实时接收的数据进行数据累计，累积8批数据得到x1(i,n)和x1_1(i,n)，x1(i,n)和x1_1(i,n)分别为累积8批得到的当前测量数据和之前数据；步骤二：对接收的数据进行傅立叶变化处理：其中，16384为傅立叶变换的长度，处理数据有87.5％的重叠，即每次更新2048点，频率分辨力为0.7324Hz；步骤三：分频段进行能量计算：步骤四：对数据进行能量判断，以确定有无信号存在；在计算传播损失时，如果是脉冲信号，需进行信噪比判断：如果Est(fb)-Ent(fb)>Th，Th为门限，则输出传播损失；否则不输出；如果是连续波，则不需要进行能量判断；步骤五：进行传播损失计算：TL(fb)＝SL(fb)－Est(fb)－HS－AG(5)其中SL为声源级，事先标定；HS表示水听器灵敏度，单位为dB；AG表示放大倍数，单位为dB。一种基于声阵的海洋声传播损失实时测量方法，其特征在于：用于波束域传播损失实时测量，利用64元细线阵或96元粗线阵输出信号测量阵波束传播损失；接收数据表示为x1_64(i,n)，i＝0,……63；x1_96(i,n)，i＝0,……95，其中，i表示水听器通道号，x1_64(i,n)表示前64个水听器接收的阵数据，x1_96(i,n)表示后96个水听器接收的阵数据；n＝0，……，2048为时间采样序列；具体包括以下步骤：步骤一：对实时接收的数据进行滤波：其中h1(n)为带通滤波器，滤波器通带为50-500Hz，过渡带为[1050Hz]和[500600Hz]，通带波纹小于等于0.5dB，阻带衰减大于等于40dB；h2(n)为带通滤波器，滤波器通带为500-2000Hz，过渡带为[400500Hz]和[20002100Hz]，通带波纹小于等于0.5dB，阻带衰减大于等于40dB；步骤二：对接收的数据进行傅立叶变换处理：步骤三：分别对两个阵的接收数据进行波束形成处理，在余弦坐标系中覆盖0°～180°水平空间：其中0≤l≤128，k＝17,...,85；其中0≤l≤128，k＝85,...,341；数据重叠50％，基元间距d1＝1.5m,d2＝0.4m,c为声速；步骤四：波束选取；根据显控输入的人工选取方位或自动选取的方位lmax选择信号波束；自动选取方位过程如下：公式(12)表示对100～500Hz频段信号计算波束图；公式(13)表示对500～2000Hz频段信号计算波束图；利用公式(14)判断选取的方位；步骤五：对选定的lmax波束数据进行逆傅立叶变换处理：其中，2048为逆傅立叶变换的长度；步骤六：对人工选取方位或自动选取的方位数据累积8批，得到的Bftime1_64(n)和Bftime1_96(n)，分别为指定方位的累积8批得到的当前测量数据，Bftime1_1_64(n)和Bftime1_1_96(n)分别为指定方位的累积8批得到的之前噪声数据；步骤七：对目标方位数据累积进行傅立叶变换处理：其中，16384为傅立叶变换的长度，处理数据有87.5％的重叠，即每次更新2048点；频率分辨力为0.7324Hz；步骤八：对波束形成并积累后的数据进行分频段能量计算：步骤九：对接收信号进行能量判断；如果是脉冲信号需进行信噪比判断：如果Est_64(fb)-Ent_64(fb)>Th，Est_96(fb)-Ent_96(fb)>Th，Th为门限，则输出传播损失；否则，不输出；如果是连续波信号不需要进行能量判断；步骤十：传播损失计算：TL_64(fb)＝SL(fb)－Est_64(fb)－HS－AG(25)TL_96(fb)＝SL(fb)－Est_96(fb)－HS－AG(26)其中SL为配合舰船发射声源或爆炸声源的声源级，事先标定；HS表示水听器灵敏度，单位为dB；AG表示放大倍数，单位为dB。本专利技术能带来以下有益效果：1、基于声阵的海洋声传播损失实时测量技术可以实时进行阵元域和波束域的多频段的水下声传播损失测量。2、基于声阵的海洋声传播损失实时测量技术在阵元域信噪比低于6dB而无法进行传播损失测量的情况下，利用声基阵的阵增益，通过波束域可以实现多频段的水下声传播损失测量。附图说明图1本专利技术的系统原理框图；图2海上试验630-800Hz传播损失实时测量结果；图3海上试验1000-1250Hz传播损失实时测量结果；图4海上试验1600-2000Hz传播损失实时测量结果。具体实施方式下面结合具体实施例子和附图对本专利技术做进一步的描述：本专利技术采用水平拖曳线列阵实现，由16段声学段组成，各声学段由其内置的水听器阵列、前置与传输部分一起，完成声信号的声电转换、前置预处理、信号采集和传输功能；阵元数为184元，为实现波束域不同频段传播损失计算，采取布阵方式是：前64个声通道，阵元间距1.5m，中间24个声通道，阵元间距24m，最后96个声通道，阵元间距0.4m；声基阵总长758.4m。图1是本专利技术的系统原理框图，从图中可以看出本专利技术的实施过程：对水平拖曳线列阵的接收阵数据，首先进行放大、滤波、进行A/D转换，将接收的模拟信号进行数字采样。然后分别进行阵元域和波束域的传播损失计算，最后根据现场测量情况，操作者可以选择进行阵元域或波束域的传播损失测量结果输出。1本文档来自技高网...

【技术保护点】
一种基于声阵的海洋声传播损失实时测量方法，其特征在于：用于单阵元传播损失实时测量，具体包括以下步骤：步骤一：对实时接收的数据进行数据累计，累积8批数据得到x1(i,n)和x1_1(i,n)，x1(i,n)和x1_1(i,n)分别为累积8批得到的当前测量数据和之前数据；步骤二：对接收的数据进行傅立叶变化处理：Y1(i,k)=FFT16384{x1(i,n)}i=0,1,2---(1)]]>Y1_1(i,k)=FFT16384{x1_1(i,n)}i=0,1,2---(2)]]>Nfft＝16384为FFT长度，处理数据有87.5％的重叠，即每次更新2048点，频率分辨力为0.7324Hz；步骤三：分频段进行能量计算：Est(fb)=13Σi=02(1NkΣk=kLkH|Y1(i,k)|2)fb=1,2,...,18---(3)]]>Ent(fb)=13Σi=02(1NkΣk=kLkH|Y1_1(i,k)|2)fb=1,2,...,18---(4)]]>步骤四：对数据进行能量判断，以确定有无信号存在；在计算传播损失时，如果是脉冲信号，需进行信噪比判断：如果Est(i,kn)‑Ent(i,kn)>Th，Th为门限，则输出传播损失；否则不输出；如果是连续波，则不需要进行能量判断；步骤五：进行传播损失计算：TL(fb)＝SL(fb)－Est(fb)－HS－AG    (5)其中SL为声源级，事先标定；HS表示水听器灵敏度，单位为dB；AG表示放大倍数，单位为dB。...

【技术特征摘要】
1.一种基于声阵的海洋声传播损失实时测量方法，其特征在于：用于单阵元传播损失实时测量，具体包括以下步骤：步骤一：对实时接收的数据进行数据累计，累积8批数据得到x1(i,n)和x1_1(i,n)，x1(i,n)和x1_1(i,n)分别为累积8批得到的当前测量数据和之前数据；步骤二：对接收的数据进行傅立叶变化处理：其中，16384为傅立叶变换的长度，处理数据有87.5％的重叠，即每次更新2048点，频率分辨力为0.7324Hz；步骤三：分频段进行能量计算：步骤四：对数据进行能量判断，以确定有无信号存在；在计算传播损失时，如果是脉冲信号，需进行信噪比判断：如果Est(fb)-Ent(fb)＞Th，Th为门限，则输出传播损失；否则不输出；如果是连续波，则不需要进行能量判断；步骤五：进行传播损失计算：TL(fb)＝SL(fb)－Est(fb)－HS－AG(5)其中SL为声源级，事先标定；HS表示水听器灵敏度，单位为dB；AG表示放大倍数，单位为dB。2.一种基于声阵的海洋声传播损失实时测量方法，其特征在于：用于波束域传播损失实时测量，利用64元细线阵或96元粗线阵输出信号测量阵波束传播损失；接收数据表示为x1_64(i,n)，i＝0,……63；x1_96(i,n)，i＝0,……95，其中，i表示水听器通道号，x1_64(i,n)表示前64个水听器接收的阵数据，x1_96(i,n)表示后96个水听器接收的阵数据；n＝0，……，2048为时间采样序列；具体包括以下步骤：步骤一：对实时接收的数据进行滤波：其中h1(n)为带通滤波器，滤波器通带为50-500Hz，过渡带为[1050Hz]和[500600Hz]，通带波纹小于等于0.5dB，阻带衰减大于等于40dB；h2(n)为带通滤波器，滤波器通带为500-2000Hz，过渡带为[400500Hz]和[20002100Hz]，通带波纹小于等于0.5dB，阻带衰减大于等于40dB；步骤二：对接收的数据进行傅立叶变换处理：1步骤三：分别对两个阵的接收数据进行波束形成处理，在余弦坐标系中覆盖0°～180°水平空间：其中0≤l≤128，k＝17,...,85；其中0≤l≤128，k＝85...

【专利技术属性】
技术研发人员：邹丽娜，王忠康，谷新禹，
申请(专利权)人：中国船舶重工集团公司第七一五研究所，
类型：发明
国别省市：浙江;33