CN112269163B - 一种基于坐底单三维矢量水听器的水中声源方位深度协同跟踪方法 - Google Patents
一种基于坐底单三维矢量水听器的水中声源方位深度协同跟踪方法 Download PDFInfo
- Publication number
- CN112269163B CN112269163B CN202011060194.6A CN202011060194A CN112269163B CN 112269163 B CN112269163 B CN 112269163B CN 202011060194 A CN202011060194 A CN 202011060194A CN 112269163 B CN112269163 B CN 112269163B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- sound
- depth
- horizontal
- sound source
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000001228 spectrum Methods 0.000 claims abstract description 85
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 4
- 230000000737 periodic effect Effects 0.000 claims description 9
- 239000013535 sea water Substances 0.000 claims description 8
- 239000002245 particle Substances 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000000926 separation method Methods 0.000 claims description 4
- 230000005236 sound signal Effects 0.000 claims description 3
- 238000010183 spectrum analysis Methods 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 abstract description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于坐底单三维矢量水听器的水中声源方位深度协同跟踪方法。利用放置于海底深度的单支三维矢量水听器,接收声源发出的声信号,生成声压谱和振速谱数据;由水平和垂直的声强互谱获得水平方位角估计值;利用估计值合成水平振速谱;由垂直声强互谱与合成水平声强互谱获得垂直到达角估计值;对声强谱进行高分辨谱估计获得调制频率,并利用调制频率获得声源深度估计结果;对声源在任一时间段的数据进行分段处理,联合绘制不同时段的水平方位角、垂直到达角和深度估计结果的跟踪曲线。本发明解决以往使用大型阵列通过空间扫描和导向空间匹配方法进行水中声源参数估计与跟踪的***复杂问题。
Description
技术领域
本发明属于水中目标跟踪领域;具体涉及一种基于坐底单三维矢量水听器的水中声源方位深度协同跟踪方法。
背景技术
对水下目标的方位、距离和深度等参数进行估计是水下目标监测、探测和跟踪的核心任务。一般是使用大型阵列,但由于***复杂度大,信息处理压力大,***海上操作和布放困难,设计更为简化的***和优秀的跟踪方法需求迫切。
相对于一般的声压水听器只能获得声场的声压信息,三维矢量水听器可共点同步测量得到声场中的声压和质点振速信息,为声场信息的全面感知和获取创造了更有利的条件,也为简化***规模提供了可能的解决方案。利用单只矢量水听器可完成对自由场条件下目标方位的高精度估计,但无法实现对水下目标更为丰富信息的估计,且直接使用与海洋信道环境下的误差较大,这是目前存在的一大问题。
结合海洋信道条件,重新构造矢量水听器接收信号的声压振速声强互谱形式,将声源的水平方位角、垂直到达角和深度信息隐含于声强互谱之中,从而更为有效的利用海洋信道条件和矢量信息,通过解算,可同时获得水平方位角、垂直到达角和深度的估计结果,这一思想既可以解决目前信道条件下矢量水听器方位估计误差大的问题,也可利用信道信息追加获得垂直到达角和深度的估计结果,得到的声源信息更为丰富。与传统的垂直阵处理模式相比,仅使用三维矢量水听器就可以获得与阵列类似的垂直到达角结果,且无需利用信道模态滤波就可获得深度估计结果,这些因素都使得该思路具有极大的新颖性。此外,由于采用了坐底的布放使用模式,***布放后具有极高的隐蔽性,矢量水听器接收***所处海洋环境的背景噪声更低,且固定坐底安装后水听器姿态稳定,接收信号质量更好。
发明内容
本发明提供了一种基于坐底单三维矢量水听器的水中声源方位深度协同跟踪方法,该方法改善了以往使用大型阵列通过空间扫描和导向空间匹配方法进行水中声源参数估计与跟踪的***复杂问题,可获得简化***、高估计精度和高跟踪效率的有益效果。且***布放与海底,稳定性和隐蔽性更优,可为水中目标的跟踪和探测提供丰富的水平方位角、垂直到达角和深度联合信息,具有广阔的应用前景。
本发明通过以下技术方案实现:
一种基于坐底单三维矢量水听器的水中声源方位深度协同跟踪方法,所述水中声源方位深度协同跟踪方法包括以下步骤:
步骤1:利用放置于海底深度的单支三维矢量水听器,接收靠近海面上部声源发出的声信号,生成声压谱和振速谱数据;
进一步的,所述步骤1具体为,
设置矢量水听器所在位置为参考坐标系原点,根据虚源理论,接收点处接收到的由(rs,zs)处声源发出的声波,其声压谱表示为,
其中,P1(rs,zs,ω)和P2(rs,zs,ω)分别表示直达波和海面反射波的声压;声速为c时,波数为k=ω/c;ω∈[ωl ωh]为角频率,ωl和ωh分别是分析频带的下限和上限;S(ω)是声源复谱振幅;η是海面反射系数,有η≈-1;Rs-和Rs+是直达波和海面反射波路径的倾斜距离,即当zb>>zs,仅考虑直达波和海面反射波的相位差,则(1)式表达为:
矢量水听器接收到的水平x方向和y方向的质点振速谱表示为:
其中,ρ为海水密度;φs为声源入射方位角;θs1和θs2分别是直达波和表面反射波的垂直到达角;由于声源和接收器之间的深度间隔大,θs1和θs2都接近θs,公式(4)近似为
矢量水听器接收到的垂直z方向的质点振速谱表示为:
进一步的,所述步骤2具体为,
水平x方向的声强互谱Ix(rs,zs,ω)和y方向声强互谱Iy(rs,zs,ω)分别为:
其中,(公式中有笔误,应该均为P(rs,zs,ω),公式(8)和(9)均已修改。)*表示复共轭算子,
其中,arctan表示反正切运算。
进一步的,所述步骤3具体为,
进一步的,所述步骤4具体为,
合成水平方向声强互谱Ir(rs,zs,ω)和垂直方向声强互谱Iz(rs,zs,ω)分别为:
进一步的,所述步骤5具体为,
利用(2)式中的声压谱近似形式,得到声强谱近似为,
上式表明,与zs、波数k和直达波和海面反射波的加权平均入射角θs有关的周期性调制项通过对声强谱的高分辨谱分析来估计,获得调制周期频率fperiod的估计结果;
本发明的有益效果是:
1.本发明利用水面及水中声源在信道中的多途声场结构,通过单三维矢量水听器接收声压和振速信号,获得声压振速声强的内在参数联系,利用声强互谱估计,获得水平方位角、垂直到达角和深度估计结果,实现对水中声源的联合跟踪。
2.本发明声强互谱估计的抗噪效果良好,适用于一般情况下的水下声场环境,由于水听器布放于海底附近,干扰小,接收***布放稳定性高,可带来更高的接收信号质量。
3.本发明相比于一般使用的大型阵列,该***组成简单、估计参数完整丰富,无需空间扫描跟踪效率高是本专利方法的关键特点,在实际应用中具有较大的优势。
4.本发明仅利用单三维矢量水听器,即可同时获得水平方位角、垂直到达角和深度的估计结果,在***规模、信息丰富程度、估计精度和速度等方面均有明显提升。适用于水下目标监测、探测和跟踪等领域。
附图说明
附图1本发明方法流程图。
附图2本发明布置情况与声线示意图。
附图3本发明的联合跟踪曲线图。
附图4本发明的水平方位角跟踪曲线图。
附图5本发明的垂直到达角跟踪曲线图。
附图6本发明的深度跟踪曲线图。
具体实施方式
下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
一种基于坐底单三维矢量水听器的水中声源方位深度协同跟踪方法,所述水中声源方位深度协同跟踪方法包括以下步骤:
步骤1:利用放置于海底深度的单支三维矢量水听器,接收靠近海面上部声源发出的声信号,生成声压谱和振速谱数据;
步骤5:对步骤2的声强谱进行高分辨谱估计获得调制频率fperiod,并利用fperiod具有与矢量水听器的布放深度zs、波数k和加权平均入射角θs有关的周期性调制关系,获得声源深度估计结果步骤1公式(2)计算声压,步骤5的公式(16)在公式(2)的基础上计算了声强,用这个声强求得深度;
进一步的,所述步骤1具体为,
海水深度为H,矢量水听器的布放深度为zb≈H,忽略海水声速分层带来的声线弯曲,并由于矢量水听器布放于海底附近而忽略海底反射波的影响,在深度为zs的水下声源产生的水下声场是直达波和海表面反射波的组合,如图1所示。
对于靠近海洋上部和底部区域的声源和水听器(声源和水听器之间的深度间隔较大,zb>>zs),仅考虑直达波和海面反射波。设置矢量水听器所在位置为参考坐标系原点,根据虚源理论,接收点处接收到的由(rs,zs)处声源发出的声波,其声压谱表示为,
其中,P1(rs,zs,ω)和P2(rs,zs,ω)分别表示直达波和海面反射波的声压;声速为c时,波数为k=ω/c;ω∈[ωlωh]为角频率,ωl和ωh分别是分析频带的下限和上限;S(ω)是声源复谱振幅;η是海面反射系数,有η≈-1;Rs-和Rs+是直达波和海面反射波路径的倾斜距离,即当zb>>zs,仅考虑直达波和海面反射波的相位差,则(1)式表达为:
矢量水听器接收到的水平x方向和y方向的质点振速谱表示为:
其中,ρ为海水密度;φs为声源入射方位角;θs1和θs2分别是直达波和表面反射波的垂直到达角;由于声源和接收器之间的深度间隔大,θs1和θs2都接近θs,公式(4)近似为
矢量水听器接收到的垂直z方向的质点振速谱表示为:
进一步的,所述步骤2具体为,
水平x方向的声强互谱Ix(rs,zs,ω)和y方向声强互谱Iy(rs,zs,ω)分别为:
其中,*表示复共轭算子,
其中,arctan表示反正切运算。
进一步的,所述步骤3具体为,
进一步的,所述步骤4具体为,
合成水平方向声强互谱Ir(rs,zs,ω)和垂直方向声强互谱Iz(rs,zs,ω)分别为:
进一步的,所述步骤5具体为,
利用(2)式中的声压谱近似形式,得到声强谱近似为,
上式表明,与zs、波数k和直达波和海面反射波的加权平均入射角θs有关的周期性调制项通过对声强谱的高分辨谱分析来估计,获得调制周期频率fperiod的估计结果;
实施例2
海水深度为500m,单三维矢量水听器布放于距离海底1m处。声源频率为200Hz,采样率为3.2kHz,水中声速1480m/s。声源位于水下深度20m处,与水听器的水平方位角为45度,声源由水平距离矢量水听器10km处向15km处移动,对运动过程中的数据进行分段处理,每次快拍为1024,以保证合适的频谱分辨率。在谱级信噪比为10dB条件下,绘制运动过程中水平方位角、垂直到达角和深度联合跟踪曲线(如图2所示),并对估计误差进行分析(如图3、4、5所示)。
对算例效果分析可知:本专利提出的方法,仅利用单支矢量水听器,即可对水中声源的水平方位角、垂直到达角和深度进行联合跟踪,***的构成简单,使用方便。跟踪过程的估计误差与真值进行对比证明,该方法有较高的精度,利用声强互谱计算的运算量小,且无需进行空间扫描,跟踪速度和效率较高。
相比于一般使用的大型阵列,该***组成简单、估计参数完整丰富,无需空间扫描跟踪效率高是本专利方法的关键特点,在实际应用中具有较大的优势。
Claims (5)
1.一种基于坐底单三维矢量水听器的水中声源方位深度协同跟踪方法,其特征在于,所述水中声源方位深度协同跟踪方法包括以下步骤:
步骤1:利用放置于海底深度的单支三维矢量水听器,接收靠近海面上部声源发出的声信号,生成声压谱和振速谱数据;
所述步骤1具体为,
设置矢量水听器所在位置为参考坐标系原点,根据虚源理论,接收点处接收到的由(rs,zs)处声源发出的声波,其声压谱表示为,
其中,P1(rs,zs,ω)和P2(rs,zs,ω)分别表示直达波和海面反射波的声压;声速为c时,波数为k=ω/c;ω∈[ωlωh]为角频率,ωl和ωh分别是分析频带的下限和上限;S(ω)是声源复谱振幅;η是海面反射系数,有η≈-1;Rs-和Rs+是直达波和海面反射波路径的倾斜距离,即当zb>>zs,仅考虑直达波和海面反射波的相位差,则(1)式表达为:
矢量水听器接收到的水平x方向和y方向的质点振速谱表示为:
其中,ρ为海水密度;φs为声源入射方位角;θs1和θs2分别是直达波和表面反射波的垂直到达角;由于声源和接收器之间的深度间隔大,θs1和θs2都接近θs,公式(4)近似为
矢量水听器接收到的垂直z方向的质点振速谱表示为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011060194.6A CN112269163B (zh) | 2020-09-30 | 2020-09-30 | 一种基于坐底单三维矢量水听器的水中声源方位深度协同跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011060194.6A CN112269163B (zh) | 2020-09-30 | 2020-09-30 | 一种基于坐底单三维矢量水听器的水中声源方位深度协同跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112269163A CN112269163A (zh) | 2021-01-26 |
CN112269163B true CN112269163B (zh) | 2023-04-25 |
Family
ID=74338523
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011060194.6A Active CN112269163B (zh) | 2020-09-30 | 2020-09-30 | 一种基于坐底单三维矢量水听器的水中声源方位深度协同跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112269163B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113109817B (zh) * | 2021-03-03 | 2021-11-09 | 中国科学院声学研究所 | 一种矢量水听器布放深度估计方法 |
CN113671442B (zh) * | 2021-07-30 | 2024-07-09 | 青岛海纳水下信息技术有限公司 | 基于矢量水听器技术的水下无人集群导航定位方法 |
CN116593965B (zh) * | 2023-07-19 | 2023-09-29 | 中国海洋大学 | 一种基于远场机会声源的水平阵阵型估计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105158734A (zh) * | 2015-07-09 | 2015-12-16 | 哈尔滨工程大学 | 一种基于阵不变量的单矢量水听器被动定位方法 |
EP3094988A1 (de) * | 2014-01-15 | 2016-11-23 | ThyssenKrupp Marine Systems GmbH | Erkennen von objekten auf see mittels mobiler sonaremitter |
CN108089155A (zh) * | 2017-12-28 | 2018-05-29 | 西北工业大学 | 一种深海环境下单水听器声源被动定位方法 |
CN111142071A (zh) * | 2020-01-08 | 2020-05-12 | 中国科学院声学研究所 | 一种结合半经验公式的单阵元匹配场***声源定位方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2442749B (en) * | 2006-10-12 | 2010-05-19 | Electromagnetic Geoservices As | Positioning system |
JP4968827B2 (ja) * | 2006-11-24 | 2012-07-04 | 国立大学法人東京海洋大学 | 水中測位システムおよび水中測位方法 |
FR2930649B1 (fr) * | 2008-04-24 | 2016-01-22 | Ixsea | Systeme de positionnement acoustique sous-marin |
CN202383285U (zh) * | 2010-12-30 | 2012-08-15 | 中国科学院声学研究所 | 一种网络化水声定位节点*** |
US9383428B2 (en) * | 2013-03-12 | 2016-07-05 | Teledyne Instruments, Inc. | Passive acoustic bearing estimation via ultra short baseline wideband methods |
CN109696657B (zh) * | 2018-06-06 | 2022-10-14 | 南京信息工程大学 | 一种基于矢量水听器的相干声源定位方法 |
CN110824429B (zh) * | 2019-10-28 | 2022-09-13 | 西北工业大学 | 深海环境下利用非同步垂直阵的宽带声源被动定位方法 |
CN111580048B (zh) * | 2020-05-09 | 2020-12-29 | 中国科学院声学研究所 | 一种利用单矢量水听器的宽带声源深度估计方法 |
CN111679248B (zh) * | 2020-05-15 | 2023-04-21 | 黑龙江工程学院 | 一种基于海底水平l型阵列的目标方位和距离联合稀疏重构定位方法 |
-
2020
- 2020-09-30 CN CN202011060194.6A patent/CN112269163B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3094988A1 (de) * | 2014-01-15 | 2016-11-23 | ThyssenKrupp Marine Systems GmbH | Erkennen von objekten auf see mittels mobiler sonaremitter |
CN105158734A (zh) * | 2015-07-09 | 2015-12-16 | 哈尔滨工程大学 | 一种基于阵不变量的单矢量水听器被动定位方法 |
CN108089155A (zh) * | 2017-12-28 | 2018-05-29 | 西北工业大学 | 一种深海环境下单水听器声源被动定位方法 |
CN111142071A (zh) * | 2020-01-08 | 2020-05-12 | 中国科学院声学研究所 | 一种结合半经验公式的单阵元匹配场***声源定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112269163A (zh) | 2021-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112269163B (zh) | 一种基于坐底单三维矢量水听器的水中声源方位深度协同跟踪方法 | |
JP2021181996A (ja) | 合成開口ソナーのためのシステムおよび方法 | |
CN112083404B (zh) | 一种基于多途特征匹配的单矢量水听器声源深度估计方法 | |
CN111580048B (zh) | 一种利用单矢量水听器的宽带声源深度估计方法 | |
CN111458008B (zh) | 基于单固定水听器和单移动声源的海洋地声参数反演方法 | |
CN112505666B (zh) | 一种二维多波束合成孔径声呐目标回波仿真方法 | |
CN112987004B (zh) | 一种浅海环境下基于水平阵列的水面水下目标分类方法 | |
CN109884730B (zh) | 一种基于深海水下滑翔机平台的虚拟阵被动探测方法 | |
CN112098938B (zh) | 一种基于六元锥矢量阵的水声目标降维匹配声场定位方法 | |
WO2012051267A2 (en) | Generating an angle domain common image gather | |
Gerstoft et al. | Adaptive beamforming of a towed array during a turn | |
CN114280541B (zh) | 一种基于深海分布式垂直线列阵的目标被动定位方法 | |
CN104714235A (zh) | 一种低频双矢量水听器阵列的测距方法及*** | |
CN111580079B (zh) | 一种基于单矢量水听器的水下运动目标最近接近距离估计方法 | |
CN116068493A (zh) | 一种深海大深度垂直分布式水听器被动声源定位方法 | |
Wu et al. | Side-scan sonar and sub-bottom profiler surveying | |
CN110907937B (zh) | 一种基于“t”型阵的掩埋物合成孔径三维成像方法 | |
CN113534161B (zh) | 一种用于远距离定位水下声源的波束镜像聚焦方法 | |
CN115561764A (zh) | 一种基于单矢量水听器的运动目标深度估计方法 | |
CN115902849A (zh) | 一种基于波束输出强度重采样的深海声源深度估计方法 | |
CN113126029B (zh) | 适用于深海可靠声路径环境的多传感器脉冲声源定位方法 | |
Li et al. | Shallow water high resolution multi-beam echo sounder | |
Feng et al. | Detection of gas leaks from sea bed using a small circular array | |
Lingsch et al. | Acoustic imagery using a multibeam bathymetric system | |
Guo et al. | Review of research on sonar imaging technology in China |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |