CN115343714A - 一种用于浅海主动声纳目标深度估计的匹配相位处理方法 - Google Patents
一种用于浅海主动声纳目标深度估计的匹配相位处理方法 Download PDFInfo
- Publication number
- CN115343714A CN115343714A CN202210799159.9A CN202210799159A CN115343714A CN 115343714 A CN115343714 A CN 115343714A CN 202210799159 A CN202210799159 A CN 202210799159A CN 115343714 A CN115343714 A CN 115343714A
- Authority
- CN
- China
- Prior art keywords
- field
- target
- interference structure
- propagation vector
- depth estimation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/06—Systems determining the position data of a target
- G01S15/08—Systems for measuring distance only
Abstract
本发明公开了一种用于浅海主动声纳目标深度估计的匹配相位处理方法,包括以下步骤,步骤1,测量场构建,将接收器1和接收器2收到的宽带回波信号进行处理得到测量场的单程传播向量比值宽带干涉结构ηdata(zt,ω);步骤2,拷贝场构建,基于海洋波导环境,利用KRAKEN模型分别计算出接收器1的声场Pr1(zs,zt,zr1,r,ω)和接收器2的声场Pr2(zs,zt,zr1,r,ω),构建拷贝场单程传播向量比值宽带干涉结构ηrplc(zt,ω);步骤3,模糊度面构建,对测量场的单程传播向量比值宽带干涉结构ηdata(zt,ω)和拷贝场单程传播向量比值宽带干涉结构ηrplc(zt,ω)进行匹配相位处理,构建目标深度估计的模糊度面。本发明可实现浅海环境下的主动声纳目标深度估计。
Description
技术领域:
本发明涉及水声领域,具体涉及一种用于浅海主动声纳目标深度估计的匹配相位处理方法。
背景技术:
海洋波导中的目标深度估计一直以来都是水声领域关注的难题之一。从被动目标深度估计研究现状看,已有的深度估计方法主要分为两类,一类是基于匹配场处理(MFP,Matched Field Processing)的深度估计方法;一类为基于匹配模处理(MMP,Matched ModeProcessing)或匹配波束处理(MBP,Matched Beam Processing)的深度估计方法。Bucker(见文献[1]H.P.Bucker.Use of calculated sound fields and matched fielddetection to locate sound sources in shallow water[J].Journal of theAcoustical Society of America,1976,59:368–373.)被认为将匹配场处理公式化的第一人,他通过声场模型生成拷贝场数据与实际接收场数据进行匹配来实现被动目标的距离和深度估计,并引入了模糊度表面的概念。Klemm(见文献[2]Klemm R.Range and depthestimation by line arrays in shallow water[J].Signal Processing,1981,3(4):333-344.)引入了高分辨的广义最大熵(ME,Maximum Entropy)波束形成器来进行深度估计,得到了超过常规MFP方法的性能。李鹏(见文献[3]Li P,Zhang X H,Li L R,etal.Source depth discrimination using wavenumber domain feature with ahorizontal array[J].Applied Acoustics,2020,164(2):107287.)等提出了一种基于水平阵列的f-k模态域分布提取方法,实现基于模态域分布特性的深度辨识。于梦枭见(见文献[4]于梦枭,周士弘,张岩,等.浅海宽带简正模相干/非相干能量比值特征匹配的源深估计[J].声学学报,2020,45(3):309-324.)等提出了一种利用单水听器接收信号中多阶简正模相干项能量与非相干项能量比值特征匹配的目标深度估计方法,消除了声源频谱对被动目标深度估计的影响,提升了深度估计的稳健性。郭良浩[11]等根据近水面声源难以激发低阶模态的物理现象,研究了利用声源波数谱结构和波数位置的不同来分辨水面声源和水下声源。Zheng(见文献[5]Zheng G Y,Zhu F W.Difference factor of vertical beampattern for shallow-water source depth discrimination[J].Acoustics Australia/Australian Acoustical Society,2021,49:105-123.)等人考虑垂直阵波束域对应的水面目标与水下目标激发声场差异,定义了一种隐含目标深度信息的特征-垂直波束谱空间差异因子,用于被动目标深度辨识,并基于Swellex-96数据验证了算法的有效性和宽容性。
与被动目标深度估计研究相比,主动声纳目标深度估计研究相对较少。Premus(见文献[6]Zheng G Y,Zhu F W.Difference factor of vertical beam pattern forshallow-water source depth discrimination[J].Acoustics Australia/AustralianAcoustical Society,2021,49:105-123.)等人将其基于模态能量分布的被动目标深度辨识方法扩展到主动声纳目标深度辨识,给出了典型的仿真分析结果。Han(见文献[7]Han N,Yao S.Discrimination of the active submerged/bottom target based on the totalscintillation index[J].Applied Acoustics,2021,172:107646.)提出了传统模态闪烁指数的变体,为主动声纳目标深度辨识引入了调制闪烁指数和总闪烁指数的数学定义,用于区分浅水波导中水面/水下目标。上述方法的前提是利用足够孔径的垂直或水平接收阵列实现模态滤波,然而对于传统的有源水平拖曳声呐平台,在考虑目标方位的条件下,其等效的水平物理孔径是不足以用于模态滤波的。因此在主动声纳目标深度估计研究方面,匹配场定位技术仍得到了一定的关注。Yang(见文献[8]Yang T C.Method and system forsensing with an active acoustic array[R].DEPARTMENT OF THE NAVY WASHINGTONDC,1996.)提出了利用多途相对幅度的宽带匹配场定位技术,采用基于接收数据特征值分解的模式分解方法实现了对目标深度和距离的估计。Ryan(见文献[9]Goldhahn RA.Waveguide invariant active sonar target detection and depth classificationin shallow water[D].Duke University,2010.)等人提出了波导不变量自适应滤波器技术,通过均匀亚采样目标回波的短时傅里叶变换(STFT,Short Time Fourier Transform)构建频域多快拍目标回波进而估计互谱密度矩阵(SDM,Spectral Density Matrix),用于主动声纳目标深度估计。然而,基于上述方法的有源匹配场处理有以下限制:(1)目标散射会改变目标与声源和接收水听器之间声传播的幅度和相位;(2)目标的散射函数通常是未知的。因此探索一种不依赖与目标散射特性的主动声纳目标深度估计方法是有意义的。
发明内容:
本发明所要解决的技术问题是,提供一种用于浅海主动声纳目标深度估计的匹配相位处理方法,以实现浅海环境下的主动声纳目标深度估计。
本发明的技术解决方案是,提供一种用于浅海主动声纳目标深度估计的匹配相位处理方法,包括以下步骤,
步骤1,测量场构建,将接收器1和接收器2收到的宽带回波信号进行处理得到测量场的单程传播向量比值宽带干涉结构ηdata(zt,ω);
步骤2,拷贝场构建,基于海洋波导环境,利用KRAKEN模型分别计算出接收器1的声场Pr1(zs,zt,zr1,r,ω)和接收器2的声场Pr2(zs,zt,zr1,r,ω),构建拷贝场单程传播向量比值宽带干涉结构ηrplc(zt,ω);
步骤3,模糊度面构建,对测量场的单程传播向量比值宽带干涉结构ηdata(zt,ω)和拷贝场单程传播向量比值宽带干涉结构ηrplc(zt,ω)进行匹配相位处理,构建目标深度估计的模糊度面。
作为优选,步骤3的匹配相位处理,通过计算数据与拷贝场的单程传播向量比值宽带干涉结构的相位特征,来进行相关运算,构建模糊度函数AP(z)如下,
作为优选,通过单个声源发射低频宽带信号,垂直双接收水听器接收目标回波。
采用以上方案后与现有技术相比,本发明具有以下优点:
本发明利用浅海波导中单个声源发射,垂直双水听器接收的低频宽带目标回波,在信号带宽内计算接收的宽带回波的比值作为测量场向量,再与利用声场模型计算的拷贝场向量进行匹配相位特征处理,实现对目标深度的估计,无需考虑目标散射特性对深度估计的影响,且运算量小,有利于实时处理的实现需求。
附图说明:
图1为本发明实施例提供的一种用于浅海主动声纳目标深度估计的匹配相位处理方法流程图;
图2为单声源发射,垂直双水听器接收探测示意图;
图3为Pekeris波导环境参数及声源和接收器布防示意图;
图4为单程传播向量宽带干涉结构的幅度特征;
图5为单程传播向量宽带干涉结构的相位特征;
图6为利用不同简正波阶数贡献的单程传播向量宽带干涉结构进行匹配,深度估计误差与信噪比的关系;
图7为不同目标深度时,匹配相位处理测深算法的模糊度表面和10次模糊度曲线。
具体实施方式:
下面结合附图就具体实施方式对本发明作进一步说明:
如图1-7所示,本实施例提供了一种用于浅海主动声纳目标深度估计的匹配相位处理方法,该方法包括以下步骤:
(1)测量场构建,将接收器1和接收器2收到的宽带回波信号进行处理得到测量场的单程传播向量比值宽带干涉结构ηdata(zt,ω);
(2)拷贝场构建,基于海洋波导环境,利用KRAKEN模型分别计算出接收器1的声场Pr1(zs,zt,zr1,r,ω)和接收器2的声场Pr2(zs,zt,zr1,r,ω),构建拷贝场单程传播向量比值宽带干涉结构ηrplc(zt,ω);
(3)模糊度面构建,对测量场的单程传播向量比值宽带干涉结构ηdata(zt,ω)和拷贝场单程传播向量比值宽带干涉结构ηrplc(zt,ω)进行匹配相位处理,构建目标深度估计的模糊度面。
主动声纳目标深度估计一直是水声领域的一个难题。针对浅海波导中主动声纳目标深度估计问题,本发明提出了通过单个声源发射低频宽带信号,垂直双接收水听器接收目标回波,利用宽带目标回波比值的相位特征进行匹配相位处理的方法,不仅消除了目标散射特性对有源声呐目标深度估计的影响,而且仅需计算信道传递函数,运算量小,有利于实时处理的实现需求。
本发明是这样实现的:
图1给出了单个声源、两个垂直接收水听器以及目标之间分布的示意图(侧视图),假设目标到声源与两个垂直接收水听器的距离均为r。
由简想正波理论可知,对于水平不变的分层波导,深度为zs的单位强度声源在距离r、深度zt处激发的声场可表示为若干简正波的叠加:
距离r、深度zt处的单位强度目标激发的声场在接收深度zri处为:
式中,krm,krn表示本征值,ψm为对应本征值krm的本征函数,ψn为对应本征值krn的本征函数,M为简正波阶数,r为目标到声源和接收水听器的距离,zs为声源深度,zri为第i个接收水听器的深度,zt为目标的深度,ρ为海水介质密度,ω为角频率。
在有源声呐目标深度估计中,目标复杂的散射特性往往是影响深度估计的重要因素。事实上,目标的散射函数T(ω,θi,θo)为频率、入射角度θi、散射角度θo的函数,考虑到声波经浅海远距离的传播,声源激发声信号传播至目标,入射的掠射角θi很小,近似于水平入射。并且经由目标散射后能远距离传播至不同深度接收水听器,散射的掠射角θo也很小。因此,在本文讨论的参数设置下做出如下假设,忽略目标的散射函数T(ω,θi,θo)中入射的掠射角θi、散射的掠射角θo的影响,将目标散射函数简化为T(ω)。
假设单个声源、两个垂直接收水听器与目标之间的水平距离无明显差异(近似相等),则由声源发射声波S(ω)传播至目标p(r,zs,zt,ω),经由目标散射T(ω)后,再传播至接收水听器p(r,zri,zt,ω),接收水听器接收声场可表示为:
其中,S(ω)表示声源发射声波的频谱,T(ω)表示目标的散射函数。
考虑到仅由同一声源发射,深度不同的两个垂直水听器接收的宽带目标回波Pri(zs,zt,zri,r,ω)中项B(zs,zt,r,ω)是一致的,由此可以利用接收到的宽带目标回波消除目标散射的影响,
将η(zt,ω)定义为目标回波的单程传播向量比值宽带干涉结构,可见上式只与接收水听器深度、目标深度、目标距离、频率等因素有关,与目标散射特性无关。这意味着,如果我们采取两个不同深度的水听器接收宽带目标回波,便可利用回波的比值消除目标散射的影响。
目标的深度信息隐含在η(zt,ω)随频率的振荡起伏分布中,由于难以总结出目标深度与η(zt,ω)宽带分布的关系,一种可行的方法是将η(zt,ω)作为隐含目标深度的特征,通过匹配特征结构的方法来实现目标深度的估计。
方法一:匹配幅度处理(MAP,Matched Amplitude Procssing)通过计算数据与拷贝场的单程传播向量比值宽带干涉结构的幅度特征,来进行相关运算,构建模糊度函数AM(z)如下,
方法二:匹配相位处理(MPP,Matched Phase Procssing)通过计算数据与拷贝场的单程传播向量比值宽带干涉结构的相位特征,来进行相关运算,构建模糊度函数AP(z)如下,
方法三:匹配幅相处理(MAPP,Matched Amplitude-Phase Procssing)通过直接计算数据与拷贝场的单程传播向量比值宽带干涉结构(包含幅度和相位信息),来进行相关运算,构建模糊度函数A(z)如下,
其中,ηdata(zt,ω),ηrplc(z,ω)分别表示数据和拷贝场的单程传播向量宽带干涉结构,模糊度函数A(z)的峰值对应的深度即为估计的目标深度。
本发明选取Pekeris波导进行仿真分析,利用KRAKEN简正波声场计算程序进行模型仿真,环境参数及声源和接收水听器布放示意图如图2所示。
仿真条件如下:水层深度为100m,水层中声速剖面为等声速,声速为1500m/s,密度为1g/cm3;沉积层深度50m,声速为1700m/s,密度为1.9g/cm3。吸收系数为0.1dB/λ。单个声源深度zs=50m,垂直接收水听器1深度zr1=40m,垂直接收水听器2深度zr2=60m,且假设目标深度zt为待估计参量,目标距离r=30km。发射信号为双曲调频(HFM)信号,信号时间长度4s,频率覆盖400Hz-800Hz,利用KRAKEN模型计算出声源至目标的格林函数p(r,zs,zt,ω)和目标至不同接收水听器的格林函数p(r,zt,zri,ω),则不同垂直接收水听器接收声场可表示为Pri(zs,zt,zri,r,ω)=S(ω)p(r,zs,zt,ω)T(ω)p(r,zt,zri,ω),设置目标的散射函数T(ω)是一个与频率有关的随机复数,S(ω)是发射信号的频谱。对接收水听器接收声场Pri(zs,zt,zri,r,ω)进行傅里叶逆变换得到接收数据的时域波形Pri(zs,zt,zri,r,t)。
信噪比(SNR,Signal-to-Noise Ratio)是评价方法性能的重要参数之一。信噪比的设定方式如下,考虑空间分布的、服从高斯分布的白噪声叠加于模拟的时域信号中,以一定信噪比对接收水听器接收的时域信号加白噪声,有
其中Pri(zs,zt,zri,r,t)与Pri(zs,zt,zri,r,ω)为傅里叶变换对,w(t)为叠加的白噪声。值得注意的是,应设置噪声与信号的带宽和时间长度一致。在后续仿真分析时,所加高斯白噪声的带宽和时间长度均与信号一致,不再重复说明。
对接收数据时域波形Pri(zs,zt,zri,r,t)叠加高斯白噪声,信噪比设置为10dB。利用垂直双接收水听器得到的宽带目标回波数据的比值消去目标散射的影响,得到单程传播向量比值宽带干涉结构。图3和图4分别给出了基于不同阶模态的单程传播向量比值宽带干涉结构的幅度特征和相位特征。可以看出,无论是单程传播向量比值宽带干涉结构的幅度特征还是相位特征,实际的声场由于全部阶模态的作用,振荡变化相对剧烈,而有限的低阶模态构建的声场则是刻画了振荡起伏的包络,即低阶模态刻画包络结构,高阶模态刻画精细结构。当选取的简正波阶数增多时,单程传播向量比值宽带干涉结构在不同深度和频率上的变化更加复杂,对目标深度的敏感性明显增强,说明选取较多的简正波阶数贡献的单程传播向量比值宽带干涉结构进行匹配处理可以提高目标深度估计的性能。
将全部阶简正波贡献的单程传播向量宽带干涉结构作为测量场,再利用KRAKEN计算出不同阶简正波贡献的单程传播向量比值宽带干涉结构作为拷贝场,利用式(7)-式(9)计算不同匹配测深算法的模糊度函数,估计目标深度。仿真条件如下:海深100m,等声速水文,环境参数,声源和接收水听器布放深度以及发射信号参数同上。目标距离r=30km,目标深度zt=50m。设置信噪比从-20dB连续变化至10dB,间隔1dB,每个信噪比下进行100次仿真试验,统计深度估计的平均误差。图5(a)-(c)给出了目标深度zt=50m时,利用5阶、10阶以及全部阶简正波贡献的单程传播向量比值宽带干涉结构进行匹配及不同匹配测深算法下的仿真估计结果。可见,
(1)对于匹配幅度处理测深算法和匹配幅相处理测深算法来说,无论利用有限阶还是全部阶的简正波,深度估计的平均误差随着信噪比增加整体上趋于不变,平均误差整体保持在20m左右,说明基于单程传播向量比值宽带干涉结构的匹配幅度处理测深算法和匹配幅相处理测深算法不适用于进行目标深度估计;
(2)对于匹配相位处理测深算法来说,深度估计平均误差随着信噪比增加整体上趋于变小,且当利用简正波阶数越多时,估计平均误差下降幅度越大。当利用全部阶简正波且信噪比大于-10dB时,估计平均误差整体上在5m之内。说明基于单程传播向量比值宽带干涉结构的匹配相位处理测深算法可用于进行目标深度估计,信噪比越高,估计的平均误差整体上越小,在无参数失配的条件下,利用简正波阶数越多,深度估计性能越好。在后续所有仿真中,选择全部阶模态贡献的匹配相位处理测深算法进行目标深度估计。
利用KRAKEN模型拷贝场划分为201×101的网格,由于在主动声呐场景中,目标距离大致获知,假设目标距离30km,匹配定位搜索范围在距离方向上为29km~31km,距离间隔10m,深度方向上为0m~100m,深度间隔1m,信噪比0dB,图6分别给出了目标距离30km处不同目标深度分别为10m,50m和90m的深度估计的模糊度面和10次模糊度曲线。可见,在不同目标深度的情况下,目标深度都得到了正确的估计,且模糊度曲线中没有其他的峰值。水体中(目标深度50m),近水面(目标深度10m)以及近海底(目标深度90m)目标的主旁瓣比都约为6dB。
以上仅就本发明较佳的实施例作了说明,但不能理解为是对权利要求的限制。凡是利用本发明说明书所做的等效结构或等效流程变换,均包括在本发明的专利保护范围之内。
Claims (3)
1.一种用于浅海主动声纳目标深度估计的匹配相位处理方法,其特征在于:包括以下步骤,
步骤1,测量场构建,将接收器1和接收器2收到的宽带回波信号进行处理得到测量场的单程传播向量比值宽带干涉结构ηdata(zt,ω);
步骤2,拷贝场构建,基于海洋波导环境,利用KRAKEN模型分别计算出接收器1的声场Pr1(zs,zt,zr1,r,ω)和接收器2的声场Pr2(zs,zt,zr1,r,ω),构建拷贝场单程传播向量比值宽带干涉结构ηrplc(zt,ω);
步骤3,模糊度面构建,对测量场的单程传播向量比值宽带干涉结构ηdata(zt,ω)和拷贝场单程传播向量比值宽带干涉结构ηrplc(zt,ω)进行匹配相位处理,构建目标深度估计的模糊度面。
3.根据权利要求1所述的用于浅海主动声纳目标深度估计的匹配相位处理方法,其特征在于:步骤1中,通过单个声源发射低频宽带信号,垂直双接收水听器接收目标回波。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210799159.9A CN115343714A (zh) | 2022-07-06 | 2022-07-06 | 一种用于浅海主动声纳目标深度估计的匹配相位处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210799159.9A CN115343714A (zh) | 2022-07-06 | 2022-07-06 | 一种用于浅海主动声纳目标深度估计的匹配相位处理方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115343714A true CN115343714A (zh) | 2022-11-15 |
Family
ID=83948015
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210799159.9A Pending CN115343714A (zh) | 2022-07-06 | 2022-07-06 | 一种用于浅海主动声纳目标深度估计的匹配相位处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115343714A (zh) |
-
2022
- 2022-07-06 CN CN202210799159.9A patent/CN115343714A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hermand et al. | Inversion of broad-band multitone acoustic data from the YELLOW SHARK summer experiments | |
JP5801527B2 (ja) | 海の魚の個体群を特徴付ける方法および装置 | |
Maguer et al. | Mechanisms for subcritical penetration into a sandy bottom: Experimental and modeling results | |
JP2008545991A5 (zh) | ||
Battle et al. | Geoacoustic inversion of tow-ship noise via near-field-matched-field processing | |
CN110058245B (zh) | 基于云模型的低频主动拖线阵声呐浅海探测效能评估方法 | |
CN107247291A (zh) | 海底浅地层声能衰减模型构建及其两个重要声能衰减特征参数提取方法 | |
Jiang et al. | Estimation of geoacoustic properties of marine sediment using a hybrid differential evolution inversion method | |
CN108398690B (zh) | 一种海底反向散射强度测量方法 | |
CN110763327A (zh) | 一种平面障板条件下单矢量水听器信号检测方法 | |
Kargl et al. | Synthetic aperture sonar imaging of simple finite targets | |
Jagannathan et al. | Scattering from extended targets in range-dependent fluctuating ocean-waveguides with clutter from theory and experiments | |
Wendelboe | Backscattering from a sandy seabed measured by a calibrated multibeam echosounder in the 190–400 kHz frequency range | |
CN117233737A (zh) | 主动目标深度估计方法、装置、计算机设备和存储介质 | |
Li et al. | Geoacoustic inversion by matched-field processing combined with vertical reflection coefficients and vertical correlation | |
CN115343714A (zh) | 一种用于浅海主动声纳目标深度估计的匹配相位处理方法 | |
Schock et al. | Spatial and temporal pulse design considerations for a marine sediment classification sonar | |
Caiti et al. | Parametric sonars for seafloor characterization | |
CN113126029A (zh) | 适用于深海可靠声路径环境的多传感器脉冲声源定位方法 | |
Arvelo Jr | Robustness and constraints of ambient noise inversion | |
RU2510045C2 (ru) | Фазовый гидролокатор бокового обзора | |
Weydert | Measurements of the acoustic backscatter of selected areas of the deep seafloor and some implications for the assessment of manganese nodule resources | |
Mikhalevsky et al. | Continuous wave and M‐sequence transmissions across the Arctic | |
Salin et al. | Measuring the characteristics of backscattering of sound on a rough surface in the near-field zone of a phased array | |
Rajan et al. | Modal inverse techniques for inferring geoacoustic properties in shallow water |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |