CN110231778B - 一种通用的uuv水下目标探测仿真方法及*** - Google Patents

一种通用的uuv水下目标探测仿真方法及*** Download PDF

Info

Publication number
CN110231778B
CN110231778B CN201910374580.3A CN201910374580A CN110231778B CN 110231778 B CN110231778 B CN 110231778B CN 201910374580 A CN201910374580 A CN 201910374580A CN 110231778 B CN110231778 B CN 110231778B
Authority
CN
China
Prior art keywords
detection
target
active
passive
propagation loss
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
Application number
CN201910374580.3A
Other languages
English (en)
Other versions
CN110231778A (zh
Inventor
孙飞虎
高善国
张扬帆
李宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Acoustics CAS
Original Assignee
Institute of Acoustics CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Acoustics CAS filed Critical Institute of Acoustics CAS
Priority to CN201910374580.3A priority Critical patent/CN110231778B/zh
Publication of CN110231778A publication Critical patent/CN110231778A/zh
Application granted granted Critical
Publication of CN110231778B publication Critical patent/CN110231778B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/66Sonar tracking systems
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B17/00Systems involving the use of models or simulators of said systems
    • G05B17/02Systems involving the use of models or simulators of said systems electric

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Automation & Control Theory (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种通用的UUV水下目标探测仿真方法和***,所述方法包括:接收选择的探测模式;所述探测模式包括:主动探测和被动探测;如果选择的是主动探测,则开启主动探测,获取目标探测结果;否则,开启被动探测,获取目标探测结果。本发明的方法和***能够提供UUV能力建设的仿真研究服务,通过对典型的水声探测装备和探测过程的建模,提供标准化、通用性水下目标探测模型。

Description

一种通用的UUV水下目标探测仿真方法及***
技术领域
本发明涉及UUV目标探测建模领域,尤其涉及一种通用的UUV水下目标探测仿真方法及***。
背景技术
无人水下航行器(UUV)是一种海上力量倍增器,拥有广泛而重要的军事、科研用途,在未来的海洋探测中发挥不可替代的作用。UUV搭载多种传感器,可以在复杂海况、恶劣环境下完成水下警戒、跟踪、勘探、中继通信等水下航行任务。随着技术的进步,UUV拥有更多的潜在应用领域,尤其是对人员要求严格的深海。
随着UUV技术的日臻成熟,为了满足水下无人探测、自主能力建设的要求, UUV开始装备多种声纳探测传感器,主要以安装在载壳上的载壳共形声传感器为主,国海军UUVBluefin 21艉部安装的主动声传感器,而美国的Progeny***公司准备在无人平台载体上加装矢量传感器用于水下探测。另外,还有一些实验性的探测声传感器装备,如SACLANTCEN和MIT合作在Odyssey UUV安装的前鼻阵声传感器,全阵共8个阵元间隔10cm工作频段在8-16kHz。2008年美国NURC将SLITA拖曳阵声传感器(31mm外径)安装在Ocean Explorer AUV上,试验UUV拖曳声传感器性能,以期完成其整个水下探测网络的设想。
在主被动多传感器探测感知方面,美国水下持续跟踪无人艇(ACTUV)的两翼安装雷声公司的MS3中频远程主被动壳体声传感器,是第五代艏部和壳体主被动声传感器***(改进SQS-56***),主动有效探测距离可以达到18km,同时具备主动和被动目标搜索能力,并通过使用光纤水听器,降低了***功耗和重量。2015年,在密西西比州海域,成功对1km外的水下目标进行了持续跟踪。但是,当前应用与 UUV平台的主被动声探测技术的成熟度不高,距离真正的实用还有差距。目前有必要抓紧开展基于仿真的通用探测方法研究,以加快装备实用化的进程。
发明内容
本发明的目的在于克服上述技术缺陷,采用的技术方案为:从仿真运行的角度,梳理UUV水声探测过程,主要理清主动声纳、被动声纳的工作逻辑,按照平台的控制要求合理地进行主动、被动探测。UUV水声探测通用模型的工作流程包括参数装订、主被动开关机控制、自主探测上报等基本流程。
为了实现上述目的,本发明提出了一种通用的UUV水下目标探测仿真方法,所述方法包括:
接收选择的探测模式;所述探测模式包括:主动探测和被动探测;
如果选择的是主动探测,则开启主动探测,获取目标探测结果;
否则,开启被动探测,获取目标探测结果。
作为上述方法的一种改进,所述主动探测具体包括:
发射主动探测信号,发射范围为左舷-15°~-115°,右舷15°~115°;
基于声场分析计算主动探测传播损失值;
基于主动声纳方程计算主动探测优质因素;
根据主动探测传播损失值和主动探测优质因素判断目标是否在探测范围内,如果在探测范围内,计算主动探测目标信息。
作为上述方法的一种改进,基于声场分析计算主动探测传播损失,具体包括:
在目标到探测源的方向上,从距离目标的距离r=1m位置处 开始,按照距离目标的距离按照距离步长Δr递进,依次计算m个位置上的传播损失值:l1,l2,...,lm;其中m为传播损失值的个数,与所选取的距离有关;
Figure BDA0002051207150000021
其中,1≤i≤m;p(1)表示距离目标距离1m处的声压值;p(r+iΔr,z)为r+iΔr 处的声压:
Figure BDA0002051207150000022
其中,z表示目标距离海平面的深度;
Figure BDA0002051207150000023
c0表示典型的相位速度,ω表示频率;αj,nj,n表示耦合系数,n表示有理函数的项数,用于逼近指数函数;
Figure BDA0002051207150000024
其中k表示波束数目,ρ表示海水密度, p(r+(i-1)Δr,z)为r+(i-1)Δr处的声压;
根据m个位置上的传播损失值,插值生成距离-传播损失值曲线。
作为上述方法的一种改进,所述主动探测优质因素为:
Figure BDA0002051207150000031
其中,FOM为主动探测优质因素;SL指的发射声源级;TS表示目标强度,由目标正横方向的目标强度和方位角计算得出;NL表示噪声级,包括环境噪声和平台自噪声;GS和GT分别表示主动声纳***的空间增益和时间增益;DT为主动声纳识别系数。
作为上述方法的一种改进,所述根据主动探测传播损失和主动探测优质因素判断目标是否在探测范围内,如果在探测范围内,计算主动探测目标信息,具体包括:
依据下述探测目标判定准则判定目标是否在探测范围内:
Figure BDA0002051207150000032
其中,TL为根据距离-传播损失值曲线得到的目标位置的传播损失值;Δa表示主动探测的工程裕量,单位为db;
如果目标在探测范围内,则计算主动探测目标信息:
目标距离为:
Figure BDA0002051207150000033
目标方位为:
Figure BDA0002051207150000034
其中,为Er声纳主动探测距离误差,Ed为目标方位估计误差;(xt,yt)表示仿真中的目标的经纬度,(xz,yz)表示UUV中探测载荷的经纬度,
Figure BDA0002051207150000035
表示由经纬度到距离的转换因子,randn表示[0,1]范围内的随机数。
作为上述方法的一种改进,所述被动探测具体包括:
基于声场分析计算被动探测传播损失值;
基于被动声纳方程计算被动探测优质因素;
根据被动探测传播损失值和被动探测优质因素判断目标是否在探测范围内,如果在探测范围内,计算被动探测目标信息。
作为上述方法的一种改进,所述基于声场分析计算被动探测传播损失值:
在目标到探测源的方向上,从距离目标的距离r=1m位置处开始,按照距离目标的距离按照距离步长Δr递进,依次计算m个位置上的传播损失值:l1,l2,...,lm;其中m为传播损失值的个数,与所选取的距离有关;
Figure BDA0002051207150000041
其中,1≤i≤m;p(1)表示距离目标距离1m处的声压值;p(r+iΔr,z)为r+iΔr 处的声压:
Figure BDA0002051207150000042
其中,z表示目标距离海平面的深度;
Figure BDA0002051207150000043
c0表示典型的相位速度,ω表示频率;αj,nj,n表示耦合系数,n表示有理函数的项数,用于逼近指数函数;
Figure BDA0002051207150000044
其中k表示波束数目,ρ表示海水密度, p(r+(i-1)Δr,z)为r+(i-1)Δr处的声压;
根据m个位置上的传播损失值,插值生成距离-传播损失值曲线。
作为上述方法的一种改进,所述被动探测优质因素为:
FOMp=SLp-NL+GSp+GTp-DTp
其中,FOMp为被动探测优质因素;SLp指的目标辐射噪声源级;GSp和GTp分别表示被动声纳***的空间增益和时间增益;DTp表示被动声纳识别系数。
作为上述方法的一种改进,所述根据被动探测传播损失和被动探测优质因素判断目标是否在探测范围内,如果在探测范围内,计算被动探测目标信息,具体包括:
依据下述探测目标判定准则判定目标是否在探测范围内:
Figure BDA0002051207150000045
其中,TLp为根据距离-传播损失值曲线得到的目标位置的传播损失值;Δp表示被动探测的工程裕量,单位为db;被动探测范围为左舷-15°~-165°,右舷15°~165°;
如果目标在探测范围内,则计算被动探测目标信息:
目标方位
Figure BDA0002051207150000051
为:
Figure BDA0002051207150000052
其中,(xt,yt)表示仿真中的目标的经纬度,(xz,yz)表示UUV中探测载荷的经纬度,
Figure BDA0002051207150000053
为声纳被动探测目标方位估计误差,randn表示[0,1]范围内的随机数。
本发明还提供了一种通用的UUV水下目标探测仿真***,所述***包括:探测模式选择模块、主动探测模块和被动探测模块;
所述探测模式选择模块,用于接收选择的探测模式;所述探测模式包括:主动探测和被动探测;如果选择的是主动探测,则开启主动探测模块,获取目标探测结果;否则,开启被动探测模块,获取目标探测结果;
所述主动探测模块,用于发射主动探测信号,基于声场分析计算主动探测传播损失值;基于主动声纳方程计算主动探测优质因素;根据主动探测传播损失值和主动探测优质因素判断目标是否在探测范围内,如果在探测范围内,计算主动探测目标信息;
所述被动探测模块,用于基于声场分析计算被动探测传播损失值;基于被动声纳方程计算被动探测优质因素;根据被动探测传播损失值和被动探测优质因素判断目标是否在探测范围内,如果在探测范围内,计算被动探测目标信息。
本发明的优势在于:
本发明的方法及***能够提供UUV能力建设的仿真研究服务,通过对典型的水声探测装备和探测过程的建模,提供标准化、通用性水下目标探测模型。
附图说明
图1为本发明的通用的UUV水下目标探测仿真方法的流程图;
图2为本发明的主被动协同探测时UUV探测结果列表;
图3为本发明的UUV探测目标方位瀑布图。
具体实施方式
下面结合附图和实例对本发明提出的方法及***进行详细描述。
实施例1
如图1所示,通用的UUV水下目标探测过程主要包括:参数装订、主被动开关机控制、自主探测上报等基本流程。从仿真运行的角度,梳理UUV水声探测过程,主要理清主动声纳、被动声纳的工作逻辑,按照平台的控制要求合理地进行主动、被动探测。
(1)参数装订
开始阶段,UUV水声探测通用模型进行初始化参数装订,通过外部模块装订目标初始信息、UUV平台初始信息、声纳探测初始信息等,同时装载环境信息。
(2)开机控制
UUV平台给出主动工作指令,水下目标探测开启主动探测;若没有开启主动探测,探测模型则开启纯被动探测。
(3)水声探测
按照主动探测方式或被动探测方式,水下目标探测开展探测工作过程。通过模拟声传播的过程,分析探测信号给出探测目标信息。
(4)探测结果上报
当连续多个环节均探测到稳定目标信号时,分析目标信息,上报UUV平台探测结果。
(5)关机控制
UUV平台给出关机指令,探测模型停止工作。
步骤1,初始化参数装订;
步骤2,判断主动探测开机命令是否开启,如开启则转到步骤3,否则转向步骤 5;
步骤3,开启主动探测,上报探测结果;
步骤4,判断是否接收到主动关机命令,如接收到,则关机,否则转到步骤2;
步骤5,判断被动探测开机命令是否开启,如开启则转到步骤6,否则转到步骤 2;
步骤6,开启被动探测,上报探测结果;
步骤7,判断是否接收到被动关机命令,如接收到,则关机,否则转到步骤2。
首先,本发明设计了适用于UUV的主动探测过程。
第一步,基于声场分析计算主动探测传播损失。
针对随距离变化情况下的声场分析,采用抛物线方程方法进行分析。假定r表示距目标点源的水平距离,z表示距离海平面的深度,θ表示方位角,p表示声压。
Figure BDA0002051207150000071
其中,ρ表示海水密度,
Figure BDA0002051207150000072
表示波束数目,ω表示频率,c表示声速,β是衰减因子。进一步,可以得到:
Figure BDA0002051207150000073
Figure BDA0002051207150000074
其中,
Figure BDA0002051207150000075
且c0表示典型的相位速度。由此,可以获得:
Figure BDA0002051207150000076
计算上式,可以获得迭代解如下:
Figure BDA0002051207150000077
其中,Δr是距离步长,n表示有理函数的项数,用于逼近指数函数,αj,nj,n表示耦合系数,其大小通过约束有理函数的准确性和稳定性来定义。
在目标到探测源的方向上,从距离目标的距离r=1m位置处 开始,按照距离目标的距离按照距离步长Δr递进,依次计算m个位置上的传播损失值:l1,l2,...,lm;其中m为传播损失值的个数,与所选取的距离有关;
Figure BDA0002051207150000078
其中,1≤i≤m;p(1)表示距离目标距离1m处的声压值;p(r+iΔr,z)为r+iΔr 处的声压:
Figure BDA0002051207150000079
其中,z表示目标距离海平面的深度;
Figure BDA0002051207150000081
c0表示典型的相位速度,ω表示频率;αj,nj,n表示耦合系数,n表示有理函数的项数,用于逼近指数函数;
Figure BDA0002051207150000082
其中k表示波束数目,ρ表示海水密度, p(r+(i-1)Δr,z)为r+(i-1)Δr处的声压;
根据m个位置上的传播损失值,插值生成距离-传播损失值曲线。
第二步,基于主动声纳方程计算优质因素。
其特征在于根据主动声纳方程,计算优质因素FOM如下:
Figure BDA0002051207150000083
其中,SL指的发射声源级;TS表示目标强度,由目标正横方向的目标强度(目标模型给出)和方位角计算得出;NL表示噪声级,包括环境噪声和平台自噪声组成; GS和GT分别表示主动声纳***的空间增益和时间增益;DT表示主动声纳识别系数,又称为检测阈。
(1)目标强度
根据目标和UUV声纳平台的经纬度,计算出相对于目标艇艏艇艉方向的方位角γ,则目标强度TS为:
TS=TS0(16.17-2.98cos2γ-3.083cos6γ)/22.233
其中,TS0为目标正横的目标强度,γ为声波入射舷角。对于收发合置的声纳,声脉冲信号的入射方向和目标艇体的夹角会影响回波信号的强度,不同入射舷角时目标的目标强度不同。
(2)海洋环境噪声级
大风大浪的恶劣条件使海洋环境噪声明显增大,从而显著降低声纳的检测性能。作为衡量海洋环境好坏的参数,海况等级按照从好至坏通常分为0~9级。
以海况和频率为参数的海洋环境噪声级,对于浅海,计算噪声谱级的近似公式为:
NL=10lgf-1.7+6S+55
式中f为频率(KHz),S为海况等级,S=0,1,2,...,9。
(3)检测阈
假设主动声纳采用匹配滤波法检测目标回波信号,则:
Figure BDA0002051207150000091
其中,T为主动声纳发射脉冲宽度,d为检测指数,根据选定的检测概率和虚警概率从接收机工作特性(ROC)曲线上查得。
(4)空间增益
GS=10lgN
其中,N表示阵元数。
(5)时间增益
主动探测通过匹配滤波处理获得时间处理增益,脉冲时间T=0.1s、B=2000Hz,可以获得的时间增益为:
GT=10lg(T·B)
其中,T表示脉冲时间,B表示脉冲频率。
第三步,输出主动探测目标信息。
首先假定声纳主动探测距离误差Er,目标方位估计误差Ed。探测目标判定准则为,目标的传播损失小于声纳的优质因素减去裕量,且目标不在探测盲区内。
Figure BDA0002051207150000092
其中,TL为根据距离-传播损失值曲线得到的目标位置的传播损失值;Δa表示主动探测的工程裕量,单位为db。
假定主动信号的发射范围为左舷-15°~-115°,右舷15°~115°。连续4个探测周期内,探测到目标后转入跟踪状态。则输出目标信息如下:
目标距离
Figure BDA0002051207150000093
目标方位
Figure BDA0002051207150000094
其中,(xt,yt)表示目标的经纬度,(xz,yz)表示UUV中探测载荷的经纬度,
Figure BDA0002051207150000095
表示由经纬度到距离(单位:米)的转换因子,randn表示[0,1]范围内的随机数。
其次,本发明设计了适用于UUV的被动目标探测过程。
第一步,基于声场分析计算被动探测传播损失TLp
这一过程与主动目标探测过程一致。
第二步,基于被动声纳方程计算优质因素。
根据被动声纳方程计算优质因素如下:
FOMp=SLp-NL+GSp+GTp-DTp
其中,SLp指的目标辐射噪声源级;GSp和GTp分别表示被动声纳***的空间增益和时间增益;DTp表示被动声纳识别系数,又称为检测阈。
(1)声源级
1)水面目标声源级
SLp=60lgν+91lg Tp-20lg f+20lg D+33.5
其中v为航速,单位节,Tp为排水量,单位吨,f为频率,单位kHz,D为距离,单位码。
2)水下目标声源级
SLp=60lgνm+91lg Tp-20lg f+35.8
其中,νm为水下目标速度,Tp为排水量,单位吨,f为噪声频率(kHz)。
(2)环境噪声级
与主动探测计算方式一致。
(3)检测阈
Figure BDA0002051207150000101
其中d为检测指数,WSL为谱线宽度或FFT分析的频率间隔,Ttotal为LOFAR 图总的观测时间。
(4)空间增益
GSp=10lgN
其中,N表示有效阵元数。频段中间,根据孔径计算有效阵元数换算。
(5)时间增益
被动探测通过能量积累获得时间处理增益,可以获得的时间增益为:
GTp=5lg(T/Δτ)
其中,T表示脉冲时间,Δτ表示脉冲积累时间。
第三步,输出被动探测目标信息。
首先假定声纳被动探测目标方位估计误差
Figure BDA0002051207150000113
探测目标判定准则为,目标的传播损失小于声纳的优质因素减去裕量,且目标不在探测盲区内。
Figure BDA0002051207150000111
其中,TLp为根据距离-传播损失值曲线得到的目标位置的传播损失值,与TL的计算方法相同;Δp表示被动探测的工程裕量,单位为db。
假定被动探测范围为左舷-15°~-165°,右舷15°~165°。连续4个探测周期内,探测到目标后转入跟踪状态。则输出目标信息如下:
目标方位
Figure BDA0002051207150000112
本发明适用于利用仿真技术实现高效的UUV探测过程,有助于UUV自主能力的仿真研究。图2给出了利用仿真数据驱动本发明方法给出的UUV探测目标的结果,图3给出了UUV探测过程中的目标方位瀑布图。从仿真结果可以看出,本发明的方法具备良好的通用性,能够模拟UUV典型水声装备的工作过程,具备简单高效的特点。
实施例2
本发明的实施例2提供了一种通用的UUV水下目标探测仿真***,所述***包括:探测模式选择模块、主动探测模块和被动探测模块;
所述探测模式选择模块,用于接收选择的探测模式;所述探测模式包括:主动探测和被动探测;如果选择的是主动探测,则开启主动探测模块,获取目标探测结果;否则,开启被动探测模块,获取目标探测结果;
所述主动探测模块,用于发射主动探测信号,基于声场分析计算主动探测传播损失值;基于主动声纳方程计算主动探测优质因素;根据主动探测传播损失值和主动探测优质因素判断目标是否在探测范围内,如果在探测范围内,计算主动探测目标信息;
所述被动探测模块,用于基于声场分析计算被动探测传播损失值;基于被动声纳方程计算被动探测优质因素;根据被动探测传播损失值和被动探测优质因素判断目标是否在探测范围内,如果在探测范围内,计算被动探测目标信息。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (8)

1.一种通用的UUV水下目标探测仿真方法,所述方法包括:
接收选择的探测模式;所述探测模式包括:主动探测和被动探测;
如果选择的是主动探测,则开启主动探测,获取目标探测结果;
否则,开启被动探测,获取目标探测结果;
所述主动探测具体包括:
发射主动探测信号,发射范围为左舷-15°~-115°,右舷15°~115°;
基于声场分析计算主动探测传播损失值;
基于主动声纳方程计算主动探测优质因素;
根据主动探测传播损失值和主动探测优质因素判断目标是否在探测范围内,如果在探测范围内,计算主动探测目标信息;
所述基于声场分析计算主动探测传播损失值,具体包括:
在目标到探测源的方向上,从距离目标的距离r=1m位置处开始,按照距离目标的距离按照距离步长Δr递进,依次计算m个位置上的传播损失值:l1,l2,...,lm;其中m为传播损失值的个数,与所选取的距离有关;
Figure FDA0002386270730000011
其中,1≤i≤m;p(1)表示距离目标距离1m处的声压值;p(r+iΔr,z)为r+iΔr处的声压:
Figure FDA0002386270730000012
其中,z表示目标距离海平面的深度;
Figure FDA0002386270730000013
c0表示典型的相位速度,ω表示频率;αj,nj,n表示耦合系数,n表示有理函数的项数,用于逼近指数函数;
Figure FDA0002386270730000014
其中k表示波束数目,ρ表示海水密度,p(r+(i-1)Δr,z)为r+(i-1)Δr处的声压;
根据m个位置上的传播损失值,插值生成距离-传播损失值曲线。
2.根据权利要求1所述通用的UUV水下目标探测仿真方法,其特征在于,所述主动探测优质因素为:
Figure FDA0002386270730000021
其中,FOM为主动探测优质因素;SL指的发射声源级;TS表示目标强度,由目标正横方向的目标强度和方位角计算得出;NL表示噪声级,包括环境噪声和平台自噪声;GS和GT分别表示主动声纳***的空间增益和时间增益;DT为主动声纳识别系数。
3.根据权利要求2所述的通用的UUV水下目标探测仿真方法,其特征在于,所述根据主动探测传播损失和主动探测优质因素判断目标是否在探测范围内,如果在探测范围内,计算主动探测目标信息,具体包括:
依据下述探测目标判定准则判定目标是否在探测范围内:
Figure FDA0002386270730000022
其中,TL为根据距离-传播损失值曲线得到的目标位置的传播损失值;Δa表示主动探测的工程裕量,单位为db;
如果目标在探测范围内,则计算主动探测目标信息:
目标距离为:
Figure FDA0002386270730000023
目标方位为:
Figure FDA0002386270730000024
其中,为Er声纳主动探测距离误差,Ed为目标方位估计误差;(xt,yt)表示仿真中的目标的经纬度,(xz,yz)表示UUV中探测载荷的经纬度,
Figure FDA0002386270730000025
表示由经纬度到距离的转换因子,randn表示[0,1]范围内的随机数。
4.根据权利要求1所述通用的UUV水下目标探测仿真方法,其特征在于,所述被动探测具体包括:
基于声场分析计算被动探测传播损失值;
基于被动声纳方程计算被动探测优质因素;
根据被动探测传播损失值和被动探测优质因素判断目标是否在探测范围内,如果在探测范围内,计算被动探测目标信息。
5.根据权利要求4所述通用的UUV水下目标探测仿真方法,其特征在于,所述基于声场分析计算被动探测传播损失值:
在目标到探测源的方向上,从距离目标的距离r=1m位置处开始,按照距离目标的距离按照距离步长Δr递进,依次计算m个位置上的传播损失值:l1,l2,...,lm;其中m为传播损失值的个数,与所选取的距离有关;
Figure FDA0002386270730000031
其中,1≤i≤m;p(1)表示距离目标距离1m处的声压值;p(r+iΔr,z)为r+iΔr处的声压:
Figure FDA0002386270730000032
其中,z表示目标距离海平面的深度;
Figure FDA0002386270730000033
c0表示典型的相位速度,ω表示频率;αj,nj,n表示耦合系数,n表示有理函数的项数,用于逼近指数函数;
Figure FDA0002386270730000034
其中k表示波束数目,ρ表示海水密度,p(r+(i-1)Δr,z)为r+(i-1)Δr处的声压;
根据m个位置上的传播损失值,插值生成距离-传播损失值曲线。
6.根据权利要求5所述的通用的UUV水下目标探测仿真方法,其特征在于,所述被动探测优质因素为:
FOMp=SLp-NL+GSp+GTp-DTp
其中,FOMp为被动探测优质因素;SLp指的目标辐射噪声源级;NL表示噪声级;GSp和GTp分别表示被动声纳***的空间增益和时间增益;DTp表示被动声纳识别系数。
7.根据权利要求6所述的通用的UUV水下目标探测仿真方法,其特征在于,所述根据被动探测传播损失和被动探测优质因素判断目标是否在探测范围内,如果在探测范围内,计算被动探测目标信息,具体包括:
依据下述探测目标判定准则判定目标是否在探测范围内:
Figure FDA0002386270730000041
其中,TLp为根据距离-传播损失值曲线得到的目标位置的传播损失值;Δp表示被动探测的工程裕量,单位为db;被动探测范围为左舷-15°~-165°,右舷15°~165°;
如果目标在探测范围内,则计算被动探测目标信息:
目标方位
Figure FDA0002386270730000042
为:
Figure FDA0002386270730000043
其中,(xt,yt)表示仿真中的目标的经纬度,(xz,yz)表示UUV中探测载荷的经纬度,
Figure FDA0002386270730000044
为声纳被动探测目标方位估计误差,randn表示[0,1]范围内的随机数。
8.一种通用的UUV水下目标探测仿真***,其特征在于,所述***包括:探测模式选择模块、主动探测模块和被动探测模块;
所述探测模式选择模块,用于接收选择的探测模式;所述探测模式包括:主动探测和被动探测;如果选择的是主动探测,则开启主动探测模块,获取目标探测结果;否则,开启被动探测模块,获取目标探测结果;
所述主动探测模块,用于发射主动探测信号,基于声场分析计算主动探测传播损失值;基于主动声纳方程计算主动探测优质因素;根据主动探测传播损失值和主动探测优质因素判断目标是否在探测范围内,如果在探测范围内,计算主动探测目标信息;
所述被动探测模块,用于基于声场分析计算被动探测传播损失值;基于被动声纳方程计算被动探测优质因素;根据被动探测传播损失值和被动探测优质因素判断目标是否在探测范围内,如果在探测范围内,计算被动探测目标信息;
所述基于声场分析计算主动探测传播损失值,具体包括:
在目标到探测源的方向上,从距离目标的距离r=1m位置处开始,按照距离目标的距离按照距离步长Δr递进,依次计算m个位置上的传播损失值:l1,l2,...,lm;其中m为传播损失值的个数,与所选取的距离有关;
Figure FDA0002386270730000051
其中,1≤i≤m;p(1)表示距离目标距离1m处的声压值;p(r+iΔr,z)为r+iΔr处的声压:
Figure FDA0002386270730000052
其中,z表示目标距离海平面的深度;
Figure FDA0002386270730000053
c0表示典型的相位速度,ω表示频率;αj,nj,n表示耦合系数,n表示有理函数的项数,用于逼近指数函数;
Figure FDA0002386270730000054
其中k表示波束数目,ρ表示海水密度,p(r+(i-1)Δr,z)为r+(i-1)Δr处的声压;
根据m个位置上的传播损失值,插值生成距离-传播损失值曲线。
CN201910374580.3A 2019-05-07 2019-05-07 一种通用的uuv水下目标探测仿真方法及*** Active CN110231778B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910374580.3A CN110231778B (zh) 2019-05-07 2019-05-07 一种通用的uuv水下目标探测仿真方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910374580.3A CN110231778B (zh) 2019-05-07 2019-05-07 一种通用的uuv水下目标探测仿真方法及***

Publications (2)

Publication Number Publication Date
CN110231778A CN110231778A (zh) 2019-09-13
CN110231778B true CN110231778B (zh) 2020-06-16

Family

ID=67860415

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910374580.3A Active CN110231778B (zh) 2019-05-07 2019-05-07 一种通用的uuv水下目标探测仿真方法及***

Country Status (1)

Country Link
CN (1) CN110231778B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110908393A (zh) * 2019-10-31 2020-03-24 中国矿业大学 基于探通一体化的水下无人航行器编队协同方法
CN112362153A (zh) * 2020-09-15 2021-02-12 中国人民解放军军事科学院国防科技创新研究院 基于uuv平台的低频主动水声探测***及方法
CN112630724B (zh) * 2020-10-30 2022-11-08 哈尔滨工程大学 一种适用于uuv平台的高分辨目标方位估计方法
CN112596412A (zh) * 2020-12-11 2021-04-02 中国科学院沈阳自动化研究所 一种多auv仿真平台

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3492347B2 (ja) * 2001-11-28 2004-02-03 三菱重工業株式会社 水中音響伝播特性模擬のための装置、方法及びプログラム
CN103257347A (zh) * 2012-02-17 2013-08-21 中国人民解放军海军装备研究院舰艇作战***论证研究所 一种实际使用环境下的声纳作用距离指标的获取方法
CN104679008A (zh) * 2013-11-27 2015-06-03 中国科学院沈阳自动化研究所 一种auv自主寻底控制方法
CN104730528A (zh) * 2013-12-19 2015-06-24 中国科学院声学研究所 一种水声多目标自主检测与方位跟踪方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3492347B2 (ja) * 2001-11-28 2004-02-03 三菱重工業株式会社 水中音響伝播特性模擬のための装置、方法及びプログラム
CN103257347A (zh) * 2012-02-17 2013-08-21 中国人民解放军海军装备研究院舰艇作战***论证研究所 一种实际使用环境下的声纳作用距离指标的获取方法
CN104679008A (zh) * 2013-11-27 2015-06-03 中国科学院沈阳自动化研究所 一种auv自主寻底控制方法
CN104730528A (zh) * 2013-12-19 2015-06-24 中国科学院声学研究所 一种水声多目标自主检测与方位跟踪方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
典型水面舰艇对潜声探测能力研究;马旭勃;《中国优秀硕士学位论文全文数据库·工程科技Ⅱ辑》;20130315(第3期);6-8、15-18、36 *
浅海声速剖面对吊放声纳探测距离的影响;曲晓慧等;《测试技术学报》;20110630;第25卷(第6期);471-476 *

Also Published As

Publication number Publication date
CN110231778A (zh) 2019-09-13

Similar Documents

Publication Publication Date Title
CN110231778B (zh) 一种通用的uuv水下目标探测仿真方法及***
US10725149B1 (en) System and method for autonomous joint detection-classification and tracking of acoustic signals of interest
EP2304467B1 (en) Method and apparatus for detection and classification of a swimming object
KR101740157B1 (ko) 소노부이의 최적운용 위치 및 수심 결정 방법
CN110703203A (zh) 基于多声学波浪滑翔机的水下脉冲声定位***
WO2016017358A9 (ja) 船舶特性推定装置、及び自動操舵装置
CN113311388B (zh) 一种水下机器人的超短基线定位***
CN110703202B (zh) 基于多声学波浪滑翔机和水面无人艇的水下脉冲声定位***
CN110058245B (zh) 基于云模型的低频主动拖线阵声呐浅海探测效能评估方法
CN101762823B (zh) 基于单向水声测距的海洋地震拖缆位置测量***
US5790474A (en) Active sonar for under-ice conditions
CN110333369B (zh) 基于水面gps校正的uuv的dvl测速***及自适应去噪方法
Terracciano et al. Passive bearing estimation using a 2-D acoustic vector sensor mounted on a hybrid autonomous underwater vehicle
CN110806760B (zh) 一种无人水下航行器的目标跟踪控制方法
CN109632258B (zh) 一种基于矢量传感器的收发分离的海洋内波声学检测方法
CN110412583B (zh) 基于多rov协同的水下不明运动目标侦测***
CN117250961A (zh) 无人移动靶船协同编队控制方法及***
CN114455042B (zh) 一种基于水下滑翔机的智能水声探测***
Avcioglu et al. Three dimensional volume coverage in multistatic sonar sensor networks
US20230213599A1 (en) Inverse estimation-based radius calculation method and system for ferromagnetic target detection
CN114852274B (zh) 一种舰船多物理场特征预报评估***
CN115100841B (zh) 一种基于电磁波的海底预置无人***远程激活方法
KR102345739B1 (ko) 표적 탐지를 할 수 있는 능동 소노부이 및 수동 소노부이쌍을 전시하는 방법, 그리고 이를 실행시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 매체
RU2735630C1 (ru) Система шумопеленгования гидроакустического комплекса подводной лодки
Klauson et al. Ship source strength estimation 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
GR01 Patent grant
GR01 Patent grant