CN109870695B - 基于深海海底反射声的非合作目标多节点水声定位方法 - Google Patents
基于深海海底反射声的非合作目标多节点水声定位方法 Download PDFInfo
- Publication number
- CN109870695B CN109870695B CN201910130546.1A CN201910130546A CN109870695B CN 109870695 B CN109870695 B CN 109870695B CN 201910130546 A CN201910130546 A CN 201910130546A CN 109870695 B CN109870695 B CN 109870695B
- Authority
- CN
- China
- Prior art keywords
- sound source
- sound
- array
- delay difference
- time delay
- 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
Images
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
基于深海海底反射声的非合作目标多节点水声定位方法,涉及水声定位领域。本发明是为了解决现有的水声定位方法中定位信号经常受海面海底反射信号的干扰,如果定位节点只有一个水听器,无法区分信号是直达波还是反射波,将会引起时延差估计错误,影响最终的水声定位结果的问题。将多个基阵作为定位节点布放在海底,监测多个定位节点在水中的坐标;采用波束形成方法测量多种声波信号到达基阵的俯仰角度,利用垂直直线阵波束形成将各基阵接收到的声源直达波和反射波分离开,然后求互相关得到时延,再进行位置解算,从而得到声源位置。它用于准确获得声源位置。
Description
技术领域
本发明涉及一种基于深海海底反射声的非合作目标多节点水声定位方法,属于水声定位领域。
背景技术
非合作目标的水声定位技术是水声研究的一个重要问题。通常的做法是通过测量目标信号到达各定位节点之间的时延差,进而求得目标与各定位节点之间的距离差。在深海的条件下,已有的方法均是针对直达声进行分析,并结合声速剖面进行声线弯曲修正。但是在实际的水声定位时,信号经常收到海面海底反射信号的干扰,如果定位节点只有一个水听器,无法区分信号是直达波还是反射波,将会引起时延差估计的错误,影响最终的水声定位结果。
发明内容
本发明是为了解决现有的水声定位方法中定位信号经常受海面海底反射信号的干扰,如果定位节点只有一个水听器,无法区分信号是直达波还是反射波,将会引起时延差估计错误,影响最终的水声定位结果的问题。现提供基于深海海底反射声的非合作目标多节点水声定位方法。
基于深海海底反射声的非合作目标多节点水声定位方法,所述方法包括以下步骤:
步骤一、将多个基阵作为定位节点布放在海底,监测多个定位节点在水中的坐标;
当声源目标发出辐射声时,每个基阵均接收到一种或者多种声波信号,采用波束形成方法测量多种声波信号到达基阵的俯仰角度,根据得到的俯仰角度区分直达波、一次海面反射波和一次海底反射波;
步骤二、对每两个基站的直达波做互相关、一次海面反射波做互相关或一次海底反射波做互相关,得到声源信号到达每两个基站之间的时延差;
步骤三、根据声源信号到达每两个基站之间的延时差、多个定位节点在水中的坐标和双曲定位方程,获得声源水平位置的初值,根据声源水平位置的初值和声源目标的深度,获得声源位置的初值;
步骤四、根据声速剖面信号和snell定律,得到声源从声源位置的初值到各基阵之间的时延差,将该时延差与步骤二中得到的时延差做差,得到时延差误差,根据时延差误差和偏微分方程,得到声源位置的偏差,用该偏差不断修正声源位置的初值,直至声源位置的偏差小于设定值,从而得到声源位置。
本发明的有益效果为:
该定位方法是在深海多途条件下,声源位置固定,利用垂直直线阵波束形成将各基阵接收到的声源直达波和反射波分离开,然后求互相关得到时延,再进行位置解算,从而得到声源位置。本申请采用多个基阵,利用波束形成方法,将直达波和反射波分离,再将分离后的直达波做互相关或对反射波做互相关,得到声源信号到达每两个基站之间的时延差;采用此方式获得的时延差准确,因此,使最终获得的声源位置准确。
附图说明
图1为具体实施方式一所述的基于深海海底反射声的非合作目标多节点水声定位方法的流程图;
图2为多节点定位几何配置示意图;
图3为单定位节点的组成;
图4为定位算法实现的流程图;
图5为深海声线传播图。
具体实施方式
具体实施方式一:参照图1至图5具体说明本实施方式,本实施方式所述的基于深海海底反射声的非合作目标多节点水声定位方法,所述方法包括以下步骤:
步骤一、将多个基阵作为定位节点布放在海底,监测多个定位节点在水中的坐标;
当声源目标发出辐射声时,每个基阵均接收到一种或者多种声波信号,采用波束形成方法测量多种声波信号到达基阵的俯仰角度,根据得到的俯仰角度区分直达波、一次海面反射波和一次海底反射波;
步骤二、对每两个基站的直达波做互相关、一次海面反射波做互相关或一次海底反射波做互相关,得到声源信号到达每两个基站之间的时延差;
步骤三、根据声源信号到达每两个基站之间的延时差、多个定位节点在水中的坐标和双曲定位方程,获得声源水平位置的初值,根据声源水平位置的初值和声源目标的深度,获得声源位置的初值;
步骤四、根据声速剖面信号和snell定律,得到声源从声源位置的初值到各基阵之间的时延差,将该时延差与步骤二中得到的时延差做差,得到时延差误差,根据时延差误差和偏微分方程,得到声源位置的偏差,用该偏差不断修正声源位置的初值,直至声源位置的偏差小于设定值,从而得到声源位置。
本实施方式中,对于每个基阵,直达波与反射波到达基阵的角度均可以测出来。以直达波为例,再次进行波束形成,只对其到达基阵的角度方向做补偿,可以得到该方向上的信号,即直达波信号;同样地,可以得到反射波信号,这样就将直达波与反射波分离开来。
具体实施方式二:本实施方式是对具体实施方式一所述的基于深海海底反射声的非合作目标多节点水声定位方法作进一步说明,本实施方式中,基阵数量为4个,4个基阵均为垂直直线阵。
本实施方式中,为实现本发明的目的,声源距水面较近,各定位节点需要布放在距离海面较远的位置。
为实现本发明的目的,定位节点由对定位用的信号有一定方位分辨能力的垂直线阵组成。
具体实施方式三:本实施方式是对具体实施方式一或二所述的基于深海海底反射声的非合作目标多节点水声定位方法作进一步说明,本实施方式中,步骤一中,采用波束形成方法测量多种声波信号到达基阵的俯仰角度的过程为:
每个垂直直线阵利用波束形成方法测量直达波、一次海面反射波和一次海底反射波的俯仰角度的方式均相同,以一个垂直直线阵为例,利用该垂直阵作波束形成可得到直达波、一次海面反射波或一次海底反射波到达接收基阵的角度为:
利用垂直直线阵对接收到的来自K个不同方向的信号进行常规波束形成,假设垂直直线阵的阵元数为N,阵元间距为d,第n号阵元接收到的信号为Xn(t),n=1,2,…,N:
式中,sk(t)表示第k个信源产生信号的时间序列,τkn1为第k个信源的信号在第n通道和第1通道的时延差,表示为:
式中,c为水中声速,θk为第k个信源信号的入射方向与水平面的夹角,向上为正,向下为负,其取值范围是-90°~90°;
假定信号入射角为θ,则信号在第n通道和第1通道的时延差为:
对各通道进行时延补偿,平均得到θ方向波束输出:
当θ=θk时,波束对准θk方向,B(θ,t)中第k个信源信号输出达到最大,定义B(θ,t)的功率谱为:
则当θ=θk时,Bp(θ)为局部最大,因此通过计算不同θ方向的波束功率输出,取局部最大值方向对应的角度为声波信号的俯仰角度。
本实施方式中,每个信源信号输出都对应一个功率谱,从得到的多个功率谱中选取多个不同的峰值(最大值),每个峰值对应一个声波信号的俯仰角度,从而得到多种声波信号的俯仰角度。
具体实施方式四:本实施方式是对具体实施方式一或二所述的基于深海海底反射声的非合作目标多节点水声定位方法作进一步说明,本实施方式中,步骤三中,获得声源位置的初值的具体过程为:
4个垂直直线阵的坐标分别为(x1,y1,z1)、(x2,y2,z2)、(x3,y3,z3)、(x4,y4,z4),声源坐标为(x,y,z),
根据双曲定位方程:
获得声源水平位置(x,y)的初值,
式中,t21为声源信号到达一号基阵与二号基阵的时延差,t31为声源信号到达一号基阵与三号基阵的时延差,t41为声源信号到达一号基阵与四号基阵的时延差。
具体实施方式五:本实施方式是对具体实施方式四所述的基于深海海底反射声的非合作目标多节点水声定位方法作进一步说明,本实施方式中,步骤三中,根据声源水平位置的初值和声源目标的深度,获得声源位置的初值的具体过程为:
根据snell定律,获得不同深度下声源目标到各定位节点的直达波、一次海面反射波或一次海底反射波的接收角度,该角度与步骤一中测量出的俯仰角度做差,差值最小时的深度值作为声源目标的深度。
具体实施方式六:本实施方式是对具体实施方式五所述的基于深海海底反射声的非合作目标多节点水声定位方法作进一步说明,本实施方式中,步骤四中,得到声源位置的具体过程为:
根据声速剖面信号和snell定律,获得声源从声源位置的初值直达基阵与声源位置的初值经一次海底反射达到基站的时延差t21′,t31′,t41′,将该时延差与步骤二中得到的时延差做差,得到时延差误差Δt21,Δt31,Δt41,将该时延差误差Δt21,Δt31,Δt41带入偏微分方程:
得到声源水平位置(x,y)的偏差(Δx,Δy),
将声源位置偏差加在声源位置初值上,利用snell定律求时延差,进而得到新的声源位置偏差,重复此过程,直到y的偏差小于设定值,从而得到声源位置。
实施例1:
基于深海海底反射声的非合作目标多节点水声定位方法,所述包括以下步骤:
步骤一、将由垂直直线阵构成的的定位节点由无人船布放到一定的海深处,并用超短基线定位***和深度计实时监测节点在水中的位置和阵型;
各定位节点利用垂直直线阵对非合作目标的信号进行检测,检测到信号后对信号进行方位估计,根据公式1至公式5提取海底反射信号的俯仰方位;
步骤二、在反射信号的俯仰方位进行波束形成,提取海底一次反射信号波形,对各节点的海底一次反射信号波形进行相关,获得时延差;
步骤三、根据定位方程和偏微分方程,获得公式6中的位置偏差与时延差误差的关系方程,假设海域声速均匀,声速值为1500米/秒,将从该关系方程中解算出的声源水平位置作为初值,初值为(1515,2965.5);
步骤四、根据声速剖面信号和snell定律,求不同深度下目标信号到定位节点的直达声与一次反射声的接收角度,与步骤一中实测的俯仰方位值比较,相差最小的深度值为目标深度,求得的声源深度为194米,根据声源水平位置的初值和声源目标的深度,获得声源位置的初值,初值位置为(1515,2965.5,194);
利用声速剖面,根据snell定律,求得该深度的声源位置到各定位节点经一次海底反射的时延值,得到信号从目标到各节点之间的时延差,并与步骤二中实测的时延差比较得到时延差误差;
将时延差误差代入公式7中的偏微分方程求得位置偏差,利用该位置偏差对初值进行修正和更新,重复步骤四操作,直到y的偏差小于1米,该位置就是目标位置估计结果。
Claims (4)
1.基于深海海底反射声的非合作目标多节点水声定位方法,其特征在于,所述方法包括以下步骤:
步骤一、将多个基阵作为定位节点布放在海底,监测多个定位节点在水中的坐标;
当声源目标发出辐射声时,每个基阵均接收到一种或者多种声波信号,采用波束形成方法测量多种声波信号到达基阵的俯仰角度,根据得到的俯仰角度区分直达波、一次海面反射波和一次海底反射波;
步骤二、对每两个基站的直达波做互相关、一次海面反射波做互相关或一次海底反射波做互相关,得到声源信号到达每两个基站之间的时延差;
步骤三、根据声源信号到达每两个基站之间的延时差、多个定位节点在水中的坐标和双曲定位方程,获得声源水平位置的初值,根据声源水平位置的初值和声源目标的深度,获得声源位置的初值;
步骤四、根据声速剖面信号和snell定律,得到声源从声源位置的初值到各基阵之间的时延差,将该时延差与步骤二中得到的时延差做差,得到时延差误差,根据时延差误差和偏微分方程,得到声源位置的偏差,用该偏差不断修正声源位置的初值,直至声源位置的偏差小于设定值,从而得到声源位置;
步骤一中,采用波束形成方法测量多种声波信号到达基阵的俯仰角度的过程为:
每个垂直直线阵利用波束形成方法测量直达波、一次海面反射波和一次海底反射波的俯仰角度的方式均相同,以一个垂直直线阵为例,利用该垂直阵作波束形成可得到直达波、一次海面反射波或一次海底反射波到达接收基阵的角度为:
利用垂直直线阵对接收到的来自K个不同方向的信号进行常规波束形成,假设垂直直线阵的阵元数为N,阵元间距为d,第n号阵元接收到的信号为Xn(t),n=1,2,…,N:
式中,sk(t)表示第k个信源产生信号的时间序列,τkn1为第k个信源的信号在第n通道和第1通道的时延差,表示为:
式中,c为水中声速,θk为第k个信源信号的入射方向与水平面的夹角,向上为正,向下为负,其取值范围是-90°~90°;
假定信号入射角为θ,则信号在第n通道和第1通道的时延差为:
对各通道进行时延补偿,平均得到θ方向波束输出:
当θ=θk时,波束对准θk方向,B(θ,t)中第k个信源信号输出达到最大,定义B(θ,t)的功率谱为:
则当θ=θk时,Bp(θ)为局部最大,因此通过计算不同θ方向的波束功率输出,取局部最大值方向对应的角度为声波信号的俯仰角度;
步骤四中,得到声源位置的具体过程为:
根据声速剖面信号和snell定律,获得声源从声源位置的初值直达基阵与声源位置的初值经一次海底反射达到基站的时延差t21′,t31′,t41′,将该时延差与步骤二中得到的时延差做差,得到时延差误差Δt21,Δt31,Δt41,将该时延差误差Δt21,Δt31,Δt41带入偏微分方程:
得到声源水平位置(x,y)的偏差(Δx,Δy),
将声源水平位置偏差加在声源位置初值上,利用snell定律求时延差,进而得到新的声源位置偏差,重复此过程,直到y的偏差小于设定值,从而得到声源位置。
2.根据权利要求1所述的基于深海海底反射声的非合作目标多节点水声定位方法,其特征在于,基阵数量为4个,4个基阵均为垂直直线阵。
4.根据权利要求3所述的基于深海海底反射声的非合作目标多节点水声定位方法,其特征在于,步骤三中,根据声源水平位置的初值和声源目标的深度,获得声源位置的初值的具体过程为:
根据snell定律,获得不同深度下声源目标到各定位节点的直达波、一次海面反射波或一次海底反射波的接收角度,该角度与步骤一中测量出的俯仰角度做差,差值最小时的深度值作为声源目标的深度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910130546.1A CN109870695B (zh) | 2019-02-21 | 2019-02-21 | 基于深海海底反射声的非合作目标多节点水声定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910130546.1A CN109870695B (zh) | 2019-02-21 | 2019-02-21 | 基于深海海底反射声的非合作目标多节点水声定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109870695A CN109870695A (zh) | 2019-06-11 |
CN109870695B true CN109870695B (zh) | 2023-02-03 |
Family
ID=66919093
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910130546.1A Active CN109870695B (zh) | 2019-02-21 | 2019-02-21 | 基于深海海底反射声的非合作目标多节点水声定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109870695B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110411480B (zh) * | 2019-08-27 | 2022-10-21 | 哈尔滨工程大学 | 一种复杂海洋环境下水下机动平台声学导航误差预测方法 |
CN111142071B (zh) * | 2020-01-08 | 2021-11-09 | 中国科学院声学研究所 | 一种结合半经验公式的单阵元匹配场***声源定位方法 |
CN112051577B (zh) * | 2020-07-14 | 2024-06-11 | 中国船舶重工集团公司第七一五研究所 | 一种利用分段式垂直阵实现深海中距目标测距及定深方法 |
CN112782704B (zh) * | 2020-12-30 | 2023-04-25 | 中国人民解放军91550部队 | 一种海面声源被动定位声波覆盖特性及定位精度分析方法 |
CN113189543B (zh) * | 2021-04-27 | 2023-07-14 | 哈尔滨工程大学 | 一种基于运动补偿鲁棒主成分分析的干扰抑制方法 |
CN113534161B (zh) * | 2021-05-31 | 2023-07-04 | 中国船舶重工集团公司第七一五研究所 | 一种用于远距离定位水下声源的波束镜像聚焦方法 |
CN113820717B (zh) * | 2021-08-23 | 2022-10-21 | 哈尔滨工程大学 | 一种基于负梯度波导下warping变换的单矢量传感器无源测距方法 |
CN116953616B (zh) * | 2023-07-17 | 2024-05-28 | 哈尔滨工程大学 | 一种对空中目标的跨冰层三维声学定位方法和装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101335574A (zh) * | 2008-07-30 | 2008-12-31 | 哈尔滨工程大学 | 多途信道中的声屏蔽方法 |
CN102253364A (zh) * | 2011-04-20 | 2011-11-23 | 电子科技大学 | 一种基于分布式照射源的非合作目标被动定位方法 |
CN102809744A (zh) * | 2012-05-15 | 2012-12-05 | 哈尔滨工程大学 | 水下远场目标被动定位过程中的近场强干扰源抑制方法 |
CN108181651A (zh) * | 2017-12-25 | 2018-06-19 | 西北工业大学 | 一种基于可靠声路径的深海地声参数反演方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4266669B2 (ja) * | 2003-02-28 | 2009-05-20 | 日本電気株式会社 | バイスタティック方位検出システム及び検出方法 |
-
2019
- 2019-02-21 CN CN201910130546.1A patent/CN109870695B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101335574A (zh) * | 2008-07-30 | 2008-12-31 | 哈尔滨工程大学 | 多途信道中的声屏蔽方法 |
CN102253364A (zh) * | 2011-04-20 | 2011-11-23 | 电子科技大学 | 一种基于分布式照射源的非合作目标被动定位方法 |
CN102809744A (zh) * | 2012-05-15 | 2012-12-05 | 哈尔滨工程大学 | 水下远场目标被动定位过程中的近场强干扰源抑制方法 |
CN108181651A (zh) * | 2017-12-25 | 2018-06-19 | 西北工业大学 | 一种基于可靠声路径的深海地声参数反演方法 |
Non-Patent Citations (2)
Title |
---|
基于三元阵声纳的强干扰抑制技术研究;吕曜辉等;《中国电子科学研究院学报》;20170831;第12卷(第4期);第406-409页 * |
基于元启发算法的纯方位被动定位方法;赵伟康等;《水下无人***学报》;20181231;第26卷(第6期);第623-627段 * |
Also Published As
Publication number | Publication date |
---|---|
CN109870695A (zh) | 2019-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109870695B (zh) | 基于深海海底反射声的非合作目标多节点水声定位方法 | |
CN108828522B (zh) | 一种利用垂直阵lcmv波束形成的水下目标辐射噪声测量方法 | |
CN108226933B (zh) | 一种基于条纹干涉结构的深海宽带目标深度估计方法 | |
CN103713285B (zh) | 基于信息融合的分布式米波阵列雷达测角方法 | |
CN107664758B (zh) | 基于长基线或超短基线组网的深海导航定位***及方法 | |
CN112782704B (zh) | 一种海面声源被动定位声波覆盖特性及定位精度分析方法 | |
CN110703203A (zh) | 基于多声学波浪滑翔机的水下脉冲声定位*** | |
CN109239686B (zh) | 一种针对分布式mimo雷达目标定位的发射机和接收机布局方法 | |
CN107976660A (zh) | 弹载多通道雷达超低空目标分析与多径回波建模方法 | |
CN111142071B (zh) | 一种结合半经验公式的单阵元匹配场***声源定位方法 | |
CN112540348A (zh) | 一种基于空间尺度上的声线修正算法在长基线水声定位***上的应用 | |
AU2008250691B2 (en) | Method for locating a sound-projecting target | |
CN113126030B (zh) | 基于宽带声场干涉结构的深海直达声区目标深度估计方法 | |
CN110471032B (zh) | 一种水下目标被动定位的方法 | |
AU2010297455B2 (en) | Method and device for measuring a profile of the ground | |
CN114563778B (zh) | 基于时延冗余测量的阵形校准方法 | |
CN109541557A (zh) | 一种雷达天线保护通道生成方法 | |
US8300501B2 (en) | Supercavitating projectile tracking system and method | |
Xu et al. | A Simulation method for USBL Localization Performance Analysis using a Stereo Array | |
RU2649899C1 (ru) | Способ измерения угловых координат групповых низколетящих целей | |
CN112180328B (zh) | 一种基于交叉定位及特征搜索的被动定位方法 | |
Yayu et al. | Research on location of underwater sound source target in deep sea sound field based on bellhop model | |
RU2820807C1 (ru) | Способ определения координат шумящих объектов с использованием вертикально развитых бортовых антенн гидроакустических комплексов | |
CN112881981B (zh) | 多种失配情况下无线传感器阵列空间处理增益损失的方法 | |
RU2761955C1 (ru) | Способ определения высоты полета низколетающей цели моноимпульсной РЛС сопровождения |
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 |