基于双水听器lofar谱图分析的目标深度分辨方法

基于双水听器lofar谱图分析的目标深度分辨方法
【专利摘要】本发明提供一种基于双水听器LOFAR谱图分析的目标深度分辨方法，包括以下步骤：1)利用深度分别处于跃层上下的两个水听器接收到的声场数据，计算得到两个深度处接收声强的LOFAR谱图；2)在得到的两个LOFAR图中分别选择一个矩形区域，两个矩形区域的坐标范围相同，对该矩形区域进行二维傅里叶变换，然后在变换域内按角度遍历的方式计算声强的积分值，得到干涉条纹角度的概率密度函数；3)计算两个概率密度函数的期望值，进而计算两个期望值之差；设置深度分辨判决门限角度值；4)根据角度期望值之差与判决门限的比较结果判断目标的深度。本发明不需要海洋环境先验信息，计算量小，能分辨具有跃层的浅海环境中的目标深度。
【专利说明】
基于双水听器LOFAR谱图分析的目标深度分辨方法
技术领域
[0001 ]本发明涉及一种目标深度分辨方法。
【背景技术】
[0002] 船只等水面声源会给水下目标监测设备带来虚警，实现对水面声源和水下目标的 有效区分是声纳设备的关键能力之一。
[0003] 现有的目标深度分辨方法主要基于匹配场处理方法以及更稳健的一些扩展方法， 例如模态闪烁和模态滤波等[Baggeroer( 1988)，Premus( 1999,2004)]。匹配场处理方法基 于确定的传播模型，将实际测量的声场数据与由传播模型计算出的拷贝场作互相关，求得 一个性能表面，估计出目标的距离和深度。基于匹配场处理的定位方法多数采用垂直阵或 水平阵，以获得足够多的声场空间信息。匹配场方法需要已知准确的海洋环境传播参数等 先验信息，以及海底底质的声学特性(例如沉积层声速、密度和衰减系数)等。模态闪烁和模 态滤波类方法则要求在海洋环境中布放多个水听器垂直阵列或水平阵列。在实际应用中， 这两类方法主要存在两个方面的缺点，一是采用的水听器阵列系统比较复杂，成本较高且 海上布放较为困难，例如由于海流的影响垂直阵会发生倾斜，并且该倾斜量一般很难精确 得到。一是匹配场处理方法对模型误差敏感，由于实际水声环境信息的难以获取及其不确 定性使得基于匹配场处理方法的目标距离和深度估计容易出现失配问题。

【发明内容】

[0004] 为了克服现有技术的不足，本发明提供一种基于双水听器L0FAR(Low Frequency Analysis Record)谱图分析的目标深度分辨方法。该方法利用了具有跃层的浅海环境中， 波导不变量的值会随着声源和接收点深度的变化而变化这一物理特性，又因为波导不变量 的值对应了接收声场L0FAR谱图中干涉条纹的斜率，因此对分别处于跃层上下两个深度处 的L0FAR谱图中的干涉条纹进行斜率分析，根据对斜率角度差值大小的判断实现对目标深 度的分辨。
[0005] 本发明解决其技术问题所采用的技术方案包括以下步骤：
[0006] 1)利用深度分别处于跃层上下的两个水听器接收到的声场数据，计算得到两个深 度处接收声强的频率-时间平面，即其L0FAR图；
[0007] 2)在步骤1得到的两个L0FAR图中分别选择一个矩形区域，两个矩形区域的时间范 围和频率范围分别相同，对该矩形区域进行二维傅里叶变换，然后在变换域内按角度遍历 的方式分别计算声强的积分值，积分值的大小反映了概率的高低，得到干涉条纹角度的概 率密度函数；
[0008] 3)计算两个概率密度函数的期望值，进一步计算两个期望值之差;设置深度分辨 判决门限角度值；
[0009] 4)根据角度期望值之差与判决门限的比较结果判断目标的深度:若期望值之差大 于判决门限，则认为目标深度较深，处于跃层之下；若期望值之差小于判决门限，则认为目 标深度较浅，处于跃层之上。
[0010] 所述步骤1)，具体实现如下：
[0011] 分别将两个水听器接收到的原始信号的采样序列分成连续的若干帧，每帧N个采 样点，根据具体情况，帧间可有部分重叠。对每帧信号采样样本L(n)作归一化和中心化处 理，归一化处理的目的是使接收信号的幅度(或方差)在时间上均匀，中心化处理是为了使 样本的均值为零。
[0012] 归一化处理
[0013] 中心化处理
[0014] 对处理后的信号x(n)作短时傅里叶变换，得到L0FAR谱图，记为Ik(?，t)，k=l，2， 其中I表示声强，《表示频率，t表示时间，下标k代表水听器编号。
[0015]所述步骤2)，具体实现如下：
[0016] 在步骤1得到的两个L0FAR图Ik(c〇，〇沽=1，2中分别选择一个矩形区域数据窗， 该数据窗的频率范围为《min< W < ?max，时间范围为tmin<t<tmax，并且满足Wmid= ( ?max+ Wmin)/2> > ( 〇max_ Wmin)和 rmid = (rmax+rmin)/2> > (rmax-rmin);对该矩形区域进行二维离 散傅里叶变换，翻频率域的L0FAR图为6 (M,v) = ￡艺/42"(-^二i，2，其 /=/| 中u = 0，l，…，M-l，v = 0，l，…，N-1;将Fk(u，v)从直角坐标系转换到极坐标系中表示为Kk(9， P)，有Q = arctan(v/u)，p = ;然后在变换域内按角度遍历的方式分别计算声强的 积分值积分值的大小反映了概率的高低，分别得到两个L0FAR图 P 中干涉条纹角度的概率密度函数
[0017] 所述步骤3)，具体实现如下：
[0018] 计算两个概率密度函数的期望值.
k=l，2，K为离散角度值的 个数。进一步计算两个期望值之差Dra= lEi-EsI ;设置深度分辨判决门限角度值Dth。
[0019] 所述步骤4)，具体实现如下
[0020] 根据角度期望值之差Dca与判决门限Dth的比较结果判断目标的深度:若期望值之 差大于判决门限D M>Dth，则认为目标深度较深，处于跃层之下；若期望值之差小于判决门 限Dra<D th，则认为目标深度较浅，处于跃层之上。
[0021]本发明的有益效果是:本发明的基本原理和具体实现得到了简正波模型的验证， 结果表明：
[0022] 1)当海洋环境中存在负声速梯度或者跃层时，波导不变量的值将会随着声源和接 收点的变化而变化，进而导致接收声场LOFAR谱图中干涉条纹的斜率发生变化。
[0023] 2)当声源位于跃层以上的较浅深度时，跃层上下两水听器的L0FAR谱图干涉条纹 斜率比较一致；当声源位于跃层以下的较深深度时，跃层上下两水听器的L0FAR谱图干涉条 纹的斜率相差较大。
[0024] 3)基于二维傅里叶变换方法的L0FAR谱图干涉条纹角度计算方法效率高，且能有 效分辨声源深度。
【附图说明】
[0025] 图1是声源处于不同深度时，声线分布情况示意图。其中，（a)是声源位于跃层以上 的情况，（b)是声源位于跃层以下的情况；
[0026] 图2是声源深度为10米时，两水听器接收声场L0FAR图。其中，（a)是水听器1，（b)是 水听器2;
[0027]图3是声源深度为50米时，两水听器接收声场L0FAR图。其中，（a)是水听器1，（b)是 水听器2;
[0028] 图4是声源位于不同深度时，两水听器L0FAR图中干涉条纹角度的概率密度分布情 况。其中，（a)是声源深度为10米的情况，（b)是声源深度为50米的情况。
【具体实施方式】
[0029] 下面结合附图和实施例对本发明进一步说明，本发明包括但不仅限于下述实施 例。
[0030] 本发明的主要内容有：
[0031] 1、本发明首次提出了采用二维傅里叶变换进行声场干涉条纹角度分析来实现目 标深度分辨的方法。
[0032] 2、给出了声场干涉条纹角度概率密度函数的计算方法。
[0033] 3、利用计算机仿真分析对比了不同深度声场干涉条纹的角度，并进行了目标深度 分辨。
[0034] 本发明的技术方案可以分成如下步骤：
[0035] 1)利用深度分别处于跃层上下的两个水听器接收到的声场数据，计算得到两个深 度处接收声强的频率-时间平面，即其L0FAR图；
[0036] 2)在步骤1得到的两个L0FAR图中分别选择一个矩形区域，两个矩形区域的时间范 围和频率范围分别相同，对该矩形区域进行二维傅里叶变换，然后在变换域内按角度遍历 的方式分别计算声强的积分值，积分值的大小反映了概率的高低，得到干涉条纹角度的概 率密度函数；
[0037] 3)计算两个概率密度函数的期望值，进一步计算两个期望值之差;设置深度分辨 判决门限角度值；
[0038] 4)根据角度期望值之差与判决门限的比较结果判断目标的深度:若期望值之差大 于判决门限，则认为目标深度较深，处于跃层之下；若期望值之差小于判决门限，则认为目 标深度较浅，处于跃层之上。
[0039] 5)性能分析。利用计算机仿真验证了本发明方法的深度分辨性能。
[0040]以下对本发明的每个步骤作进一步的详细说明：
[0041 ]所述步骤1)，具体实现如下：
[0042] 分别将两个水听器接收到的原始信号的采样序列分成连续的若干帧，每帧N个采 样点，根据具体情况，帧间可有部分重叠。对每帧信号采样样本L(n)作归一化和中心化处 理，归一化处理的目的是使接收信号的幅度(或方差)在时间上均匀，中心化处理是为了使 样本的均值为零。
[0043] 归一化处理！ (!)
[0044] 中心化处理： （2)
[0045] 对处理后的信号x(n)作短时傅里叶变换，得到L0FAR谱图，记为Ik(?，t)，k=l，2， 其中I表示声强，《表示频率，t表示时间，下标k代表水听器编号。
[0046]所述步骤2)，具体实现如下：
[0047] 在步骤1得到的两个L0FAR图Ik(c〇，〇沽=1，2中分别选择一个矩形区域数据窗， 该数据窗的频率范围为《min< W < ?max，时间范围为tmin<t<tmax，并且满足
[0048] Wmid-( 〇 max+〇 min)/2〉〉（ W max-〇 min) ( 3 )
[0049] rmid= (rmax+rmin)/2> > (rmax-rmin) (4)
[0050] 对该矩形区域进行二维离散傅里叶变换，得到频率域的LOFAR图为
[0051] pkM = E Z7, {〇),t)e-JMualM+vtlN\k = 1,2 (5) 0)=0)] t=t\
[0052] 其中u = 0,1，…，M-l，v = 0,1，…，N-l;将Fk(u，v)从直角坐标系转换到极坐标系中 表示为Fk(Q，P)，有Q = arctan(vAi) ^ = yju2 +V2 ;然后在变换域内按角度遍历的方式分别 计算声强的积分值= 2积分值的大小反映了概率的高低，分别得到两
P
个L0FAR图中干涉条纹角度的概率密度函数 [0053]所述步骤3)，具体实现如下：
[0054]计算两个概率密度函数的期望僅 ，k=l，2，K为离散角度值的 个数。进一步计算两个期望值之差Dra= lEi-EsI ;设置深度分辨判决门限角度值Dth。
[0055]所述步骤4)，具体实现如下
[0056]根据角度期望值之差Dca与判决门限Dth的比较结果判断目标的深度:若期望值之 差大于判决门限DM>Dth，则认为目标深度较深，处于跃层之下；若期望值之差小于判决门 限D ca<Dth，则认为目标深度较浅，处于跃层之上。实际应用中，判决门限可以根据感兴趣目 标的速度、频带等特性灵活设置，以达到最优的深度分辨目的。
[0057] 所述步骤5)，具体实现如下：
[0058] 以存在负声速梯度的典型浅海环境为例，给出了本发明的实施实例。
[0059]实例参数设置如下:假设海水中声速剖面以及目标和水听器的布放情况如图1中 所示。其中，海水混合层深度为15m，声速为1540m/s;温跃层厚度为10m;声速经过温跃层后， 变为1510m/s，并随着深度的增加不再变化。水听器1和水听器2的布放深度分别为10m和 50m。假设目标辐射150-500HZ的宽带随机信号，速度为5m/s，距离水听器最近点距离为lkm， 水听器开始记录的时间为目标经过最近点后200秒。
[0060] 本实例分为两种情况:情况1，目标深度为10m，如图1(a)所示;情况2,目标深度为 50m，如图1(b)所示。
[0061] 本实例的具体实现过程如下：
[0062] 1)获取两个水听器布放位置处的L0FAR图
[0063] 根据图1中的海洋环境模型，利用KRAKEN程序计算水听器1和水听器2接收的声场 L0FAR谱图。情况1的声场L0FAR谱图如图2所示，情况2的声场L0FAR谱图如图3所示。
[0064] 2)计算L0FAR图中干涉条纹角度的概率密度函数
[0065] 对L0FAR图中时间大于500秒，频率大于200Hz的矩形区域进行二维傅里叶变换(矩 形区域的选择，只影响角度的分散程度，不影响最终判定结果），在变换后的极坐标系中按 角度遍历的方式分别计算声强的积分值，得到L0FAR图中干涉条纹角度的概率密度函数。情 况1的概率密度函数如图4(a)所示，情况2的概率密度函数如图4(b)所示。
[0066] 3)计算两水听器概率密度函数期望值之差，并设置判决门限
[0067]情况1中，计算得到水听器1对应的干涉条纹角度期望值为70.1°，水听器2对应的 干涉条纹角度期望值为67.9°，可得其角度差为Dc；ai = 2.2°。
[0068]情况2中，计算得到水听器1对应的干涉条纹角度期望值为67.9°，水听器2对应的 干涉条纹角度期望值为57.5°，可得其角度差为0。32 = 10.4°。
[0069]将判决门限设置为Dth = 5°。
[0070] 4)分辨目标的深度
[0071 ] 情况1中，Dral<Dth，判断目标深度较浅，处于跃层之上。
[0072] 情况2中，Dra2>Dth，判断目标深度较深，处于跃层之下。
[0073] 5)性能分析。
[0074]在负跃层的声速剖面下，声场简正波可以划分为水面反射-海底反射简正波 (SRBR，即反射简正波)和水体反转-海底反射简正波(RBR，即反转简正波)两大类，不同类型 的简正波对应的波导不变量典型值见表1。当反射简正波占主导地位时，波导不变量的值为 +1(略小于1)，当反转简正波占主导地位时，波导不变量的值为负值或者大于1的正值。
[0075]表1负跃层波导中简正波的分类及其对应的波导不变量典型值
[0077]如图1(a)所示，当声源深度较浅位于跃层以上时，激发的简正波类型为水面反射- 海底反射SRBR型。水听器1和水听器2接收到的均为SRBR型简正波，对应的波导不变量值均 为略小于1的正值，因而形成的干涉条纹角度也较为一致。
[0078]如图1(b)所示，当声源深度较深位于跃层以下时，激发的低阶简正波多为RBR型 (图中虚线示意），高阶简正波多为SRBR型（图中实线示意）。水听器1接收到的为SRBR型简正 波，对应的波导不变量值为略小于1的正数。水听器2接收到的主要为RBR型简正波，对应的 波导不变量值大于1或者为负数。此时，水听器1和2接收到声场干涉条纹角度将会有较大的 不同。
[0079]因此，通过判断两个不同深度处水听器的L0FAR图中干涉条纹的斜率可以实现对 目标深度的分辨。实际应用中，判决门限可以根据感兴趣目标的速度、频带等特性灵活设 置，以达到最优的深度分辨目的。
【主权项】
1. 基于双水听器LOFAR谱图分析的目标深度分辨方法，其特征在于包括下述步骤： 1) 利用深度分别处于跃层上下的两个水听器接收到的声场数据，计算得到两个深度处 接收声强的频率-时间平面，即其L0FAR图； 2) 在步骤1得到的两个L0FAR图中分别选择一个矩形区域，两个矩形区域的时间范围和 频率范围分别相同，对该矩形区域进行二维傅里叶变换，然后在变换域内按角度遍历的方 式分别计算声强的积分值，积分值的大小反映了概率的高低，得到干涉条纹角度的概率密 度函数； 3) 计算两个概率密度函数的期望值，进一步计算两个期望值之差;设置深度分辨判决 门限角度值； 4) 根据角度期望值之差与判决门限的比较结果判断目标的深度:若期望值之差大于判 决门限，则认为目标深度较深，处于跃层之下；若期望值之差小于判决门限，则认为目标深 度较浅，处于跃层之上。2. 根据权利要求1所述的基于双水听器L0FAR谱图分析的目标深度分辨方法，其特征在 于所述的步骤1)具体实现如下： 分别将两个水听器接收到的原始信号的采样序列分成连续的若干帧，每帧N个采样点， 根据具体情况，帧间可有部分重叠;对每帧信号采样样本L(n)作归一化和中心化处理，归一 化处理的目的是使接收信号的幅度(或方差)在时间上均匀，中心化处理是为了使样本的均 值为零；对处理后的信号x(n)作短时傅里叶变换，得到L0FAR谱图，记为Ik(?，〇汰=1，2，其中1 表示声强，ω表示频率，t表示时间，下标k代表水听器编号。3. 根据权利要求1所述的基于双水听器L0FAR谱图分析的目标深度分辨方法，其特征在 于所述的步骤2)具体实现如下： 在得到的两个L0FAR图11{(?，〇沽=1，2中分别选择一个矩形区域数据窗，该数据窗的 频率祀1 围为 CO min< CO < CO max，时间:围为tmin< t < tmax，并且俩足 〇 mid - ( 〇 max+ 〇 min) / 2〉 > (ω max- ω min )和rmid = ( rmax+rmin ) /2 > > ( rmax-rmin );对该矩形区域进行二维离散傅里叶变 换，得到频率域的L0FAR图为4=1，2，其中11 = 0， 1，…，M-1，v = 0,1，···，N-1;将Fk(u，v)从直角坐标系转换到极坐标系中表示为Fk(9，P)，有Θ =arctan (v/u);然后在变换域内按角度遍历的方式分别计算声强的积分值'积分值的大小反映了概率的高低，分别得到两个L0FAR图中干涉 条纹角度的概率密度函数4. 根据权利要求1所述的基于双水听器LOFAR谱图分析的目标深度分辨方法，其特征在 于所述的步骤3)具体实现如下： 计算两个概率密度函数的期望值>k=l，2，K为离散角度值的个数； 进一步计算两个期望值之差Dra= lEi-EsI ;设置深度分辨判决门限角度值Dth。5. 根据权利要求1所述的基于双水听器LOFAR谱图分析的目标深度分辨方法，其特征在 于:所述的步骤4)根据角度期望值之差与判决门限D th的比较结果判断目标的深度:若期 望值之差大于判决门限则认为目标深度较深，处于跃层之下；若期望值之差小于 判决门限D ra<Dth，则认为目标深度较浅，处于跃层之上。
【文档编号】G01S7/539GK105929385SQ201610223586
【公开日】2016年9月7日
【申请日】2016年4月6日
【发明人】杜金燕, 李晔, 崔浩
【申请人】山东省科学院海洋仪器仪表研究所