CN115598627A - 一种收发分置单基地主动声纳界面混响信号仿真方法 - Google Patents

一种收发分置单基地主动声纳界面混响信号仿真方法 Download PDF

Info

Publication number
CN115598627A
CN115598627A CN202211220364.1A CN202211220364A CN115598627A CN 115598627 A CN115598627 A CN 115598627A CN 202211220364 A CN202211220364 A CN 202211220364A CN 115598627 A CN115598627 A CN 115598627A
Authority
CN
China
Prior art keywords
scattering
reverberation
point
receiving
unit
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
Application number
CN202211220364.1A
Other languages
English (en)
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.)
Dalian Huahai Zhikong Electronic Information Technology Co ltd
Original Assignee
Dalian Huahai Zhikong Electronic Information Technology Co ltd
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 Dalian Huahai Zhikong Electronic Information Technology Co ltd filed Critical Dalian Huahai Zhikong Electronic Information Technology Co ltd
Priority to CN202211220364.1A priority Critical patent/CN115598627A/zh
Publication of CN115598627A publication Critical patent/CN115598627A/zh
Pending legal-status Critical Current

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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/523—Details of pulse systems
    • G01S7/526—Receivers
    • G01S7/527—Extracting wanted echo signals
    • 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/003—Bistatic sonar systems; Multistatic sonar systems
    • 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/88—Sonar systems specially adapted for specific applications
    • 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/537—Counter-measures or counter-counter-measures, e.g. jamming, anti-jamming

Landscapes

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

Abstract

本发明涉及水声工程技术领域，提供一种收发分置单基地主动声纳界面混响信号仿真方法，包括：建立收发分置单基地声纳海底、海面散射的最小散射椭球确定方法和后续散射椭球的确定方法；提出基于界面散射椭圆环的散射单元划分方法；建立散射单元边界点、散射单元中心点坐标的计算模型；建立散射单元发射入射角、接收散射角、到收‑发点的距离、混响时延等计算模型；建立海底、海面混响信号的计算模型；更新界面混响计算结束条件。本发明能够计算海底海面散射产生的综合混响信号，极大提高了收发分置单基地工作模式下混响到达时间、混响强度，以及多普勒计算等的精度，增强了混响模型的可靠性。

Description

一种收发分置单基地主动声纳界面混响信号仿真方法
技术领域
本发明涉及水声工程技术领域，尤其涉及一种收发分置单基地主动声纳界面混响信号仿真方法。
背景技术
海洋混响是主动声纳发射信号后由海底、海面，以及海洋水体不均匀性等共同产生的散射信号的总和，是一种不可避免的水声现象。无论在浅海还是深海，由海底和海面散射引起的界面混响是主动声纳的主要干扰信号，对主动目标回波信号检测有着重要影响。混响信号建模和仿真是水声领域的重要课题，它主要有两个方面的用途：一是用于主动声纳抗混响信号处理技术验证和完善；二是用于主动声纳系统推演和评估，评估声纳系统在混响环境下的探测效果。
收发分置的单基地声纳系统是新型主动声纳常用的形式，这种声纳系统的发射阵和接收阵由同一平台搭载，发射单元与接收阵列独立工作，两者相距一定距离(一般100-600米)。因此，它既不同于传统的单基地(发射单阵和接收阵在同一平台上、且共用，位置相同)，也不同于多基地声纳(发射与接收分别位于不同平台，两者距离的配置一般多比较远，几公里以上)。目前，主动声纳混响建模和仿真主要是针对单基地收发合置声纳或双基地(收-发分置距离远)两种情况，为了简化计算，通常将其按收发合置的单基地处理的，严重影响仿真的可靠性。现有混响信号仿真方法主要存在以下问题：
1、收发分置距离可达100-600米的单基地声纳，发射信号到达散射体(海底、海面和水体)的入射角与散射体到达接收点的散射角大不相同，利用常规单基地概念进行混响建模，仅考虑单程双倍几何传播损失，使散射信号强度、散射信号到达时间产生误差，特别在浅海环境下，影响更加严重；
2、浅海主动声纳混响主要由海底散射引起的，海底混响信号的多普勒特性仿真不可忽略，由于浅海环境下收发分置角较大，用收发合置进行混响信号的多普勒建模，会造成混响信号的多普勒仿真误差，会影响推演评估效果可靠性，甚至无法用此对具有抗混响信号处理能力的声纳进行推演评估。
3、双基地声纳混响信号建模方法，因为双基地收-发单元距离是几公里、十几公里，甚至更远，这种建模方法不能用于收发距离较近的分置式单基地。而且，现有的双基地混响信号建模方法，大多没有考虑发射垂直波束角的控制，在垂直方向上是全向发射，这样对混响强度、混响信号到达时间均有较大影响，影响混响信号仿真的可靠性和精度；
4、现有的混响模型都是针对传统的短脉冲声纳(脉冲长度τ一般是几十毫秒至几百毫秒)，进行建模的。对于短脉冲混响计算，散射环之间的距离间隔可以按照cτ/2计算，这样得到的混响信号的瞬时值符合高斯分布、包络符合瑞利分布的特性。然而，新型主动声纳的脉冲长度越来越长，脉冲长度一般在1-10多秒，甚至更长。这种情况下，如果按照现有的方法推算散射环的间隔，得到的混响信号会呈现明显的锯齿状，不能满足瞬时值符合高斯分布、包络符合瑞利分布的混响信号基本特性要求。
鉴于以上现实需求和存在问题，特提出本发明专利申请。
发明内容
本发明主要解决现有技术的声纳混响信号仿真方法存在的上述四项技术问题，提出一种收发分置单基地主动声纳界面混响信号仿真方法，以达到提高混响到达时间、混响强度，增强混响模型的可靠性，保证混响强度的准确性，以及更好地满足混响模型计算精度和实时性等工程化实现的需要的目的。
本发明提供了一种收发分置单基地主动声纳界面混响信号仿真方法，包括：
建立收发分置单基地声纳界面散射椭球模型；
基于所述散射椭球模型确定海底/海面散射的最小散射椭球，并根据所述最小散射椭球得到后续的散射椭球；
将每个所述散射椭球按方位进行划分，形成距离由散射椭圆环确定、方向由方位区域确定的若干散射单元，并计算各所述散射单元的边界点和中心点坐标；
根据所述散射单元的边界点和中心点坐标，计算散射单元的发射入射角、散射角、接收散射角、混响时延，以及各散射单元到收-发点的距离和所在区域面积；
建立多普勒频移模型，计算所述散射单元处散射回波的归一化双基地多普勒频移、时域多普勒因子；
根据所述散射单元的发射入射角、散射角、接收散射角、混响时延、到收-发点的距离及所在区域面积，计算所述散射单元的散射损失；根据所述散射损失和所述多普勒频移、时域多普勒因子，计算所述散射单元的混响信号；
当计算得到的所述混响信号满足混响级低于海洋环境噪声级与基阵指向性DI之差值时，所述散射单元的混响仿真程序结束。
进一步地，所述建立收发分置单基地声纳界面散射椭球模型，包括：
设X轴沿拖曳平台的运动方向为正，发射阵在前、接收阵在后，z轴朝水面方向为正、海底方向为负；T、R分别为发射阵和接收阵的位置，发射阵与接收阵深度分别为hT和hR，收-发站间的斜距离为2L，其连线与水平面的倾斜角为δ。
进一步地，所述基于所述散射椭球模型确定海底/海面散射的最小散射椭球，包括：
基于发射垂直波束为±Ψ/2的限定条件下，对应收发分置单基地声纳，存在一海底/海面无声散射圆，所述无声散射圆是一个以T′为圆心，以Ly为半径的圆，海底无声散射圆的半径表示为：
Ly＝(H-hT)tan(90-Ψ/2)＝(H-hT)ctan(Ψ/2) (1-1)
海面无声散射圆的半径表示为：
Ly＝hTtan(90-Ψ/2)＝hTctan(Ψ/2) (1-2)
以Ly为半径的无声散射圆与T'-R'连线内的交点作为散射点；
发射声到达海底散射点的最小斜距RTmin可表示为：
Figure BDA0003877683210000031
发射声到达海面散射点的最小斜距RTmin可表示为：
Figure BDA0003877683210000032
海底散射点到达接收点的斜距RRmin可表示为：
Figure BDA0003877683210000033
海面散射点到达接收点的斜距RRmin可表示为：
Figure BDA0003877683210000034
其中，2L'为发射点T与接收点R之间的水平距离，即T′、R′之间的距离；
最先到达接收点的的海底/海面散射椭球为最小散射椭球，其半长轴表示为：
Figure BDA0003877683210000041
对应的两个半短轴为：
Figure BDA0003877683210000042
式(4)和式(5)中，a0、b0、c0定义了所述最小散射椭球。
进一步地，所述根据所述最小散射椭球得到后续的散射椭球，包括：
所述后续的散射椭球中第1个散射椭球外环的半长轴和半短轴分别为：
Figure BDA0003877683210000043
第i＝1,2,....个散射椭球外环的半长轴和半短轴分别为：
Figure BDA0003877683210000044
进一步地，所述将每个所述散射椭球按方位进行划分，形成距离由散射椭圆环确定、方向由方位区域确定的若干散射单元，并计算各所述散射单元的边界点和中心点坐标，包括：
所述散射椭球按方位划分的方法是：若某一方位范围内的多普勒偏移量的变化在允许范围内，则该方位区域就确定为一独立散射单元格；
计算各所述散射单元的边界点的方法如下：
散射椭球在海底/海面的截面为xoy平面上的椭圆，椭圆方程为：
Figure BDA0003877683210000045
式中，对于海底混响计算时，通过第i＝1,2,....个散射椭球在海底z＝-(H-H0)的截面来计算；对于海面混响计算时，通过第i＝1,2,....个散射椭球在海底z＝H0的截面来计算；
将每个散射椭圆划分成K个区域，则第i个散射椭圆外环上、第k个散射单元边界点Bik的极角为ρbk，第k个散射单元的坐标为bxik,byik，通过式(9)-(10)来计算：
bxik＝byikcot(ρbk)或byik＝bxiktan(ρbk) (9)
Figure BDA0003877683210000051
计算各所述散射单元的中心点坐标的方法如下：
设Sik为第i个散射椭圆环上、第k个散射单元的散射中心位置，其坐标用(sxik,syki,z)表示，第i＝1,2,....个散射椭圆环中心环的半长轴和半短轴分别为：
Figure BDA0003877683210000052
设第k个散射单元中心的极角为ρk，则Sik的坐标为：
Figure BDA0003877683210000053
或
Figure BDA0003877683210000054
上式(12)的两个公式之一与椭圆方程联合，即可求解所述散射单元的中心点坐标(sxik,syik)，表示为：
Figure BDA0003877683210000055
进一步地，所述根据所述散射单元的边界点和中心点坐标，计算散射单元的发射入射角、散射角、接收散射角、混响时延，以及各散射单元到收-发点的距离和所在区域面积，并计算所述散射单元处散射回波的归一化双基地多普勒频移，包括：
(1)建立散射单元距离及其传播时延计算模型：
坐标为(-Lcosδ，0，Lsinδ)的发射点T到散射单元Sik中心点的距离为：
Figure BDA0003877683210000056
Sik到坐标为(Lcosδ,0,-Lsinδ)的接收点R的距离为：
Figure BDA0003877683210000057
由发射点T到海底散射单元Sik，再到接收点R的传播时间为：
Figure BDA0003877683210000058
由于垂直发射波束的限制，在发射点正下方海底存在声散射影区，该区域声波不能到达，这个区域不发生散射，在后面散射椭球扩展计算散射强度的过程中，不应计入，对应平坦海底，基本判据是：
RTik>RTmin (17)
(2)散射单元的入射角、散射角计算模型：
T点入射到散射单元Sik的垂直入射角为：
海底散射：
Figure BDA0003877683210000061
海面散射：
Figure BDA0003877683210000062
散射单元Sik到接收点R的垂直散射角为：
海底散射：
Figure BDA0003877683210000063
海面散射：
Figure BDA0003877683210000064
(3)各个散射单元所在区域面积的计算模型：
根据计算得到的第i个椭圆环、第k个散射单元的边界点坐标为(bxik,byik)，进一步计算第i个椭圆环、第k个扇区的面积：
Figure BDA0003877683210000065
则散射单元Sik的面积为：
Abik＝Aik-A(i-1)k
Figure BDA0003877683210000066
计算扇区面积时，在同一椭圆按照下标号递增进行，考虑区域的对称性，只需要计算1/4个扇区、面积相同区域直接采用即可。
进一步地，所述建立多普勒频移模型，计算所述散射单元处散射回波的归一化双基地多普勒频移、时域多普勒因子，包括：
根据同一拖曳平台运动情况下，发射站速度、接收站运动速度与平台运动速度v相同，则散射点Sik处散射回波的归一化双基地多普勒频率表示为：
Figure BDA0003877683210000067
时域多普勒因子表示为：
Figure BDA0003877683210000068
式中，αik为T到达散射单元Sik的掠射角、βik为散射单元Sik到R的掠射角；θTik为T'到达散射单元Sik的舷角；θRik为R'到达散射单元Sik的舷角。
进一步地，所述根据所述散射单元的发射入射角、散射角、接收散射角、混响时延、到收-发点的距离及所在区域面积，计算所述散射单元的散射损失，包括：
在深海，经海底/海面散射单元Sik的总散射损失为：
TLrik(t)＝TSbik(t)-20lgRTik(t)-20lg RRik(t)-10-3α(f)(RTik(t)+RRik(t))(24-1)
在浅海，扩展损失服从1.5次方规律，即介于平面波和球面波之间，经海底/海面散射单元Sik的总散射损失为：
TLrik(t)＝TSbik(t)-15lgRTik(t)-15lg RRik(t)-10-3α(f)(RTik(t)+RRik(t))(24-2)
其中，TSbik为散射单元Sik的海底散射强度，表示为：
TSbik＝Sbik+10lg Abik (25-1)
对于海面散射，式(24)中散射单元Sik的TSbik由海面散射强度TSsik代替：
TSsik＝Ssik+10lg Asik (25-2)
式(27)中的Abik和Asik为第i个散射环第k个海底、海面散射单元的面积，Sbik、Ssik分别为海底散射强度、海面散射系数，由(26)式表示：
Sbik＝10lg(μsinαiksinβik) (26-1)
Figure BDA0003877683210000071
式中，β＝158(vf1/3)-0.58，αik是该散射单元入射波的掠射角，βik是该散射单元散射波的掠射角，v表示海面风速，f表示中心频率；
α(f)为海水的声吸收系数，对于中心频率利用Thorp经验公式表示为：
Figure BDA0003877683210000072
散射单元Sik散射信号的散射损失为：
Figure BDA0003877683210000073
进一步地，所述根据所述散射损失和所述多普勒频移、时域多普勒因子，计算所述散射单元的混响信号，包括：
假设发射信号为声源级为SL的信号s(t)，不考虑发射水平指向性，发射垂直指向性按照前面的限定，并考虑散射单元的多普勒影响，T点发射、通过海底/海面散射单元Sik到达接收点R的散射信号可表示为：
Figure BDA0003877683210000074
或：
Figure BDA0003877683210000081
式中，fdik为多普勒频移，Δik为对应的多普勒因子(时域表示)，φik为服从均匀分布的随机相位，在0-2π间均匀分布，τik为信号发射到第Sik散射单元、再返回接收点总声程的往返时延，按式(16)计算；
第i个散射椭圆环的散射信号计算模型，表示为：
Figure BDA0003877683210000082
ABi为第i个散射环产生的散射信号的随机幅度，服从正态分布；
海底/海面散射的混响信号计算模型，表示为：
海底混响信号：
Figure BDA0003877683210000083
海面混响信号：
Figure BDA0003877683210000084
Nb(t)、Ns(t)分别为t时刻海底散射信号、海面散射信号到达接收端的散射椭圆环的总数；
将式(31-1)和(31-2)同时刻的信号相加，即得到海底和海面共同产生的总的混响信号：
y(t)＝yb(t)+ys(t) (32)。
进一步地，所述当计算得到的所述混响信号满足混响级低于海洋环境噪声级的条件时，所述散射单元的混响仿真程序结束，包括：
考虑到主动声纳接收机的阵处理能力，增加阵列处理的指向性DI增益，以使得界面混响与海洋混响衔接更加平滑，当满足式(33)的条件时，所述散射单元的混响仿真程序结束；
Figure BDA0003877683210000085
本发明提供的一种收发分置单基地主动声纳界面混响信号仿真方法，通可以综合得到海底、海面界面散射产生的混响信号，改变了当前将其作为收发合置单基地或多基地进行混响建模的现状，极大改进了这种收发分置单基地工作模式下，提高了混响到达时间、混响强度，以及界面混响信号多普勒仿真计算等的精度，增强了混响模型的可靠性；
本发明基于散射椭球模型，提出了形成收发分置单基地声纳主动混响的最小椭球及其后续散射椭圆环的确定方法，提出了海底(海面)混响信号最早到达时间(混响时延)的确定方法。混响到达时间体现了在接收点到达的混响信号的位置，时延精度直接影响混响信号的时间和能量关系，对任何收发分置系统的混响信号建模都十分重要；
本发明针对垂直发射波束受限制情况，在动态计算计算过程中，提出了一种散射单元是否位于垂直波束范围内的判据。如果在垂直波束范围内，则参与混响强度计算，否则不参与计算，从而更好地保证混响强度准确性；
本发明提出了一种散射椭圆单元划分方法，从保证混响信号瞬时幅度满足高斯分布、包络符合瑞利分布的要求，给出了声传播引起的散射椭球扩展的长轴增量(对于散射椭圆环范围)确定原则，解决了长脉冲引起的界面散射环范围如何确定、无据可依的问题；
本发明提出了每个散射环按方位间隔划分散射单元格的方法，可以更好地符合水平全向发射阵(即：发射阵没有水平波束)的特点，方位间隔可根据界面混响多普勒特性对方位分辨精度和混响信号计算实时性等方面的要求加以设定，更好地满足混响模型计算精度和实时性等工程化实现的需要。
附图说明
图1是本发明的方法流程图。
图1-1是实施例中，海深200米，考虑“无声散射圆”的界面混响图(从左至右依次为模型计算的界面混响强度图、界面混响信号图、混响信号计算的混响强度图)；
图1-2是实施例中，海深200米，没有考虑“无声散射圆”的界面混响图(从左至右依次为界面混响信号图、混响信号计算的混响强度图)；
图1-3是实施例中，海深5000米，考虑“无声散射圆”的界面混响图(从左至右依次为模型计算的界面混响强度图、界面混响信号图、混响信号计算的混响强度图)；
图1-4是实施例中，海深5000米，没有考虑“无声散射圆”的界面混响图(从左至右依次为界面混响信号图、混响信号计算的混响强度图)；
图1-5是实施例中，200米海深海底散射单元边界点(实心点)、中心点(空心圈)分布图；
图1-6是实施例中，200米海深海面散射单元边界点(实心点)、中心点(空心圈)分布图；
图1-7是实施例中，归一化界面混响信号瞬时幅值直方图和高斯分布对比图；
图1-8是实施例中，归一化界面混响信号包络直方图和瑞利分布对比图；
图2-1是实施例中，收发分置海底散射垂直剖面几何示意图；
图2-2是实施例中，收-发分置散射单元三维示意图；
图2-3是实施例中，散射椭圆环与散射单元划分示意图；
图2-4是实施例中，收发分置混响信号计算流程图。
具体实施方式
为使本发明解决的技术问题、采用的技术方案和达到的技术效果更加清楚，下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部内容。
如图1所示，本发明实施例提供的收发分置单基地主动声纳界面混响信号仿真方法包括：
101、建立收发分置单基地声纳界面散射椭球模型；
具体而言，声纳系统的几何配置垂直剖面如图2-1所示。设X轴沿拖曳平台的运动方向为正，发射阵在前、接收阵在后，z轴朝水面方向为正、海底方向为负；T、R分别为发射阵和接收阵的位置，发射阵与接收阵深度分别为hT和hR，收-发站间的斜距离为2L，其连线与水平面的倾斜角为δ。
借鉴单元散射进行混响建模思想，这样的散射点满足“发射-散射单元-接收”的声传播时间相同条件。根据几何学原理，在三维空间，收发分置声纳系统时间同步到达的散射点，是以收-发分置连线中点O为中心、以T和R为焦点的椭球，而该椭球在海底或海面上相交的椭圆，即为海底或海面散射椭圆。假如海洋声速不随声传播方向、距离变化，则从发射点T到散射椭圆上任意点再到接收点R的声传播时间均相同，简称为“时同椭球”。这是收发分置系统建立散射模型的基础。T′、R′表示发站、收站在海底投影平面上的垂足，O'为等时散射椭球在海底投影形成的椭圆的中心。Ψ/2为由发射垂直波束(±Ψ/2)决定的最大掠射角，由T沿该立体角发出的声线，在海底形成了一个以T′为圆心、以|AB|为直径的圆，该圆称为“无声散射圆”，也就是说声能不达到该圆弧内，或者说该圆内所有面积上没有散射发生。收发分置单基地主动声纳三维散射如图2-2所示。
102、基于散射椭球模型确定海底/海面散射的最小散射椭球，并根据最小散射椭球得到后续的散射椭球；
具体而言，(1)确定海底、海面散射的最小散射椭球
确定最小散射椭球是混响建模的重要基础，因为这涉及最先到达散射波的时间，也因此影响后续散射波叠加的信号能量。
基于发射垂直波束为±Ψ/2的限定条件下，对应收发分置单基地声纳，存在一海底/海面无声散射圆，所述无声散射圆是一个以T′为圆心，以Ly为半径的圆，海底无声散射圆的半径表示为：
Ly＝(H-hT)tan(90-Ψ/2)＝(H-hT)ctan(Ψ/2) (1-1)
海面无声散射圆的半径表示为：
Ly＝hTtan(90-Ψ/2)＝hTctan(Ψ/2) (1-2)
以Ly为半径的无声散射圆与T'-R'连线内的交点作为散射点；
发射声到达海底散射点的最小斜距RTmin可表示为：
Figure BDA0003877683210000111
发射声到达海面散射点的最小斜距RTmin可表示为：
Figure BDA0003877683210000112
海底散射点到达接收点的斜距RRmin可表示为：
Figure BDA0003877683210000113
海面散射点到达接收点的斜距RRmin可表示为：
Figure BDA0003877683210000114
其中，2L'为发射点T与接收点R之间的水平距离，即T′、R′之间的距离；
无论收-发单元怎样配置，必定存在一个与海底(或海面)刚刚接触的椭球，最早接触海底(或海面)的只是一个点，也即存在一个最先到达接收点的散射点，该点确定的海底(或海面)散射椭球，实际上就是最小散射椭球。以此接触点为起点的椭球，继续扩展椭球长轴形成的椭球，与海底(或海面)截面就形成了海底(或海面)散射椭圆，以一定半长轴间距Δa扩展的两个椭球在海底(或海面)的散射椭圆就形成了一个散射椭圆环。
最小散射椭球及散射点的确定。最先到达接收点的的海底/海面散射椭球为最小散射椭球，表示为：
Figure BDA0003877683210000121
对应的两个半短轴为：
Figure BDA0003877683210000122
式(4)和式(5)中，a0、b0、c0定义了所述最小散射椭球。
(2)确定散射椭球间隔和后续散射椭球
常规收发合置单基地声纳通常发射的是小于百毫秒级的短脉冲信号，混响建模按照半径间隔cτ/2的散射圆计算。但对于几秒～10多秒长脉冲信号，按此计算的散射圆半径间隔是几公里～十多公里，每个散射环后续相继到达的散射信号时延是几秒～10多秒，由此产生的混响信号会呈明显的锯齿状，混响信号幅度和包络的统计特性不满足混响信号的特征。
本发明通过对浅海、深海混响的仿真研究，对收发分置后续散射椭球环间隔(散射椭球半长轴a间隔)，按照Δa＝0.1c～0.2c米扩展(c为声速)，由此计算得到的混响信号满足声纳混响信号的基本特性要求。
根据最小散射椭球得到后续的散射椭球，包括：后续的散射椭球中第1个散射椭球外环的半长轴和半短轴分别为：
Figure BDA0003877683210000123
第i＝1,2,....个散射椭球外环的半长轴和半短轴分别为：
Figure BDA0003877683210000124
103、将每个散射椭球按方位进行划分，形成距离由散射椭圆环确定、方向由方位区域确定的若干散射单元，并计算各散射单元的边界点和中心点坐标；
具体而言，(1)散射椭球按方位划分的方法是：若某一方位范围内的多普勒偏移量的变化在允许范围内，则该方位区域就确定为一独立散射单元格。因为混响的多普勒右散射体的相对位置(方位、距离)决定的。理论上，方位区域范围越小，计算越精确，单计算量越大，影响计算的实时性。假如，发射阵具有水平波束扫描能力，则方位区域大小由水平扫描波束宽度确定；否则，方位区域大小由多普勒特性对方位变化的敏感性确定，即：若某一方位范围内的多普勒偏移量的变化在允许范围内，则该方位区域就确定为一独立散射单元格。
(2)计算各散射单元的边界点的方法如下：
散射椭球在海底/海面的截面为xoy平面上的椭圆，椭圆方程为：
Figure BDA0003877683210000131
式中，对于海底混响计算时，通过第i＝1,2,....个散射椭球在海底z＝-(H-H0)的截面来计算；对于海面混响计算时，通过第i＝1,2,....个散射椭球在海底z＝H0的截面来计算；
将每个散射椭圆划分成K个区域，则第i个散射椭圆外环上、第k个散射单元边界点Bik的极角为ρbk，第k个散射单元的坐标为bxik,byik，通过式(9)-(10)来计算：
bxik＝byikcot(ρbk)或byik＝bxiktan(ρbk) (9)
Figure BDA0003877683210000132
(3)计算各散射单元的中心点坐标的方法如下：
设Sik为第i个散射椭圆环上、第k个散射单元的散射中心位置，其坐标用(sxik,syki,z)表示，第i＝1,2,....个散射椭圆环中心环的半长轴和半短轴分别为：
Figure BDA0003877683210000133
设第k个散射单元中心的极角为ρk，则Sik的坐标为：
Figure BDA0003877683210000134
或
Figure BDA0003877683210000135
上式(12)的两个公式之一与椭圆方程联合，即可求解所述散射单元的中心点坐标(sxik,syik)，表示为：
Figure BDA0003877683210000136
104、根据散射单元的边界点和中心点坐标，计算散射单元的发射入射角、散射角、接收散射角、混响时延，以及各散射单元到收-发点的距离和所在区域面积；
具体而言，(1)建立散射单元距离及其传播时延计算模型：
坐标为(-Lcosδ，0，Lsinδ)的发射点T到散射单元Sik中心点的距离为：
Figure BDA0003877683210000141
Sik到坐标为(Lcosδ,0,-Lsinδ)的接收点R的距离为：
Figure BDA0003877683210000142
由发射点T到海底散射单元Sik，再到接收点R的传播时间为：
Figure BDA0003877683210000143
由于垂直发射波束的限制，在发射点正下方海底存在声散射影区，该区域声波不能到达，这个区域不发生散射，在后面散射椭球扩展计算散射强度的过程中，不应计入，对应平坦海底，基本判据是：
RTik>RTmin (17)
(2)散射单元的入射角、散射角计算模型：
T点入射到散射单元Sik的垂直入射角为：
海底散射：
Figure BDA0003877683210000144
海面散射：
Figure BDA0003877683210000145
散射单元Sik到接收点R的垂直散射角为：
海底散射：
Figure BDA0003877683210000146
海面散射：
Figure BDA0003877683210000147
(3)各个散射单元所在区域面积的计算模型：
如图2-3所示，根据计算得到的第i个椭圆环、第k个散射单元的边界点坐标为(bxik,byik)，进一步计算第i个椭圆环、第k个扇区的面积：
Figure BDA0003877683210000148
则散射单元Sik的面积为：
Abik＝Aik-A(i-1)k
Figure BDA0003877683210000149
计算扇区面积时，在同一椭圆按照下标号递增进行，考虑区域的对称性，只需要计算1/4个扇区、面积相同区域直接采用即可。
105、建立多普勒频移模型，计算散射单元处散射回波的归一化双基地多普勒频移、时域多普勒因子；
具体而言，根据同一拖曳平台运动情况下，发射站速度、接收站运动速度与平台运动速度v相同，则散射点Sik处散射回波的归一化双基地多普勒频率表示为：
Figure BDA0003877683210000151
时域多普勒因子表示为：
Figure BDA0003877683210000152
式中，αik为T到达散射单元Sik的掠射角、βik为散射单元Sik到R的掠射角；θTik为T'到达散射单元Sik的舷角；θRik为R'到达散射单元Sik的舷角。显然，不同的方位具有不同的多普勒频移。
106、根据散射单元的发射入射角、散射角、接收散射角、混响时延、到收-发点的距离及所在区域面积，计算散射单元的散射损失；根据散射损失和多普勒频移、时域多普勒因子，计算散射单元的混响信号；
具体而言，在深海，经海底/海面散射单元Sik的总散射损失为：
TLrik(t)＝TSbik(t)-20lgRTik(t)-20lg RRik(t)-10-3α(f)(RTik(t)+RRik(t))(24-1)
在浅海，扩展损失服从1.5次方规律，即介于平面波和球面波之间，经海底/海面散射单元Sik的总散射损失为：
TLrik(t)＝TSbik(t)-15lgRTik(t)-15lg RRik(t)-10-3α(f)(RTik(t)+RRik(t))(24-2)
其中，TSbik为散射单元Sik的海底散射强度，表示为：
TSbik＝Sbik+10lg Abik (25-1)
对于海面散射，式(24)中散射单元Sik的TSbik由海面散射强度TSsik代替：
TSsik＝Ssik+10lg Asik (25-2)
式(27)中的Abik和Asik为第i个散射环第k个海底、海面散射单元的面积，Sbik、Ssik分别为海底散射强度、海面散射系数，由(26)式表示：
Sbik＝10lg(μsinaiksinβik) (26-1)
Figure BDA0003877683210000153
式中，β＝158(vf1/3)-0.58，αik是该散射单元入射波的掠射角，βik是该散射单元散射波的掠射角v表示海面风速(单位是节)，f表示中心频率(单位Hz)。
α(f)为海水的声吸收系数，对于中心频率利用Thorp经验公式表示为：
Figure BDA0003877683210000161
散射单元Sik散射信号的散射损失为：
Figure BDA0003877683210000162
假设发射信号为声源级为SL的信号s(t)，不考虑发射水平指向性，发射垂直指向性按照前面的限定，并考虑散射单元的多普勒影响，T点发射、通过海底/海面散射单元Sik到达接收点R的散射信号可表示为：
Figure BDA0003877683210000163
或：
Figure BDA0003877683210000164
式中，fdik为多普勒频移，Δik为对应的多普勒因子(时域表示)，φik为服从均匀分布的随机相位，在0-2π间均匀分布，τik为信号发射到第Sik散射单元、再返回接收点总声程的往返时延，按式(16)计算；
第i个散射椭圆环的散射信号计算模型，表示为：
Figure BDA0003877683210000165
ABi为第i个散射环产生的散射信号的随机幅度，服从正态分布；
海底/海面散射的混响信号计算模型，表示为：
海底混响信号：
Figure BDA0003877683210000166
海面混响信号：
Figure BDA0003877683210000167
Nb(t)、Ns(t)分别为t时刻海底散射信号、海面散射信号到达接收端的散射椭圆环的总数；
将式(31-1)和(31-2)同时刻的信号相加，即得到海底和海面共同产生的总的混响信号：
y(t)＝yb(t)+ys(t) (32)。
107、当计算得到的所述混响信号满足混响级低于海洋环境噪声级与基阵指向性DI之差值时，所述散射单元的混响仿真程序结束。
具体而言，N(t)的确定：随着掠射角越来越小、散射元越来越远，散射强度越来越弱，仿真中并不需要计算垂直开角内所有掠射角度，因此需要判断何时终止计算。
基本准则：是混响级低于海洋环境噪声级(同主动声纳工作频段)，即当由混响掩蔽级为主变为噪声掩蔽级为主时，混响仿真程序结束。考虑到主动声纳接收机的阵处理能力，增加阵列处理的指向性DI增益，以使得界面混响与海洋混响衔接更加平滑，当满足式(33)的条件时，所述散射单元的混响仿真程序结束；
Figure BDA0003877683210000171
总体流程如图2-4。
按照图2-4进行混响信号仿真计算。具有实现过程如下：
1、设置混响计算的背景参数
(1)主动发射信号为频率范围为900Hz-1500Hz的HFM信号、脉冲长度3秒。
(2)发射声源级SL＝210dB，根据声源级计算发射信号声压幅度；
(3)根据采样频率fs＝25000H在，产生信号波形s(t)。
(4)发射垂直开角(±Ψ/2)：全向，±Ψ/2＝±90°；受限：±Ψ/2＝±30°。
(5)海区深度H：(1)浅海200米水深；(2)深海5000米水深。
(6)二级海况，海面风速v＝13.3kn。
(7)发射阵深度hT＝70米，接收阵深度hR＝50米，H0＝(70+50)/2＝60米。
(8)发射-接收点之间的斜距2L＝400米。
(9)扩展椭球半长轴的增量Δa＝150米。
(10)声速c＝1500米/秒。
(11)每个散射单元数量为K＝12。
(12)海底散射反向散射强度10logμ＝-27。
2、确定海底z＝-(H-H0)散射最小散射椭球及其后续散射椭球
(1)按照公式(1)-(5)确定最小散射椭球参数；
(2)按照公式(6)-(8)确定后续第i＝1,2,....个海底散射的椭球参数。
3、对所有海底散射椭球环，进行以下循环计算：
3.1计算海底散射椭圆单元边界点坐标
将海底散射椭圆环，分成K＝12个方位区域，每个区域的方位角大小Δρ＝360/K＝30°，则第i个散射椭圆外环上、第k个散射单元边界点Bik的极角为ρbk＝0°，3°0，60°，…，330°，联立(9)～(10)式，计算得到散射单元格的边界坐标(bxik,byik)。
3.2计算海底第i个散射椭圆环各散射单元的中心点坐标
按照公式(11)，计算第i＝1,2,....个散射椭球的中心环的半长轴和半短轴；
联立(12)-(13)式，计算第i个散射椭圆环上、第k个散射单元(方位区域)的散射中心位置Sik的坐标(sxik,syki)，计算式中的ρk＝ρbk+15°。
3.3计算散射单元距离及其传播时延等参数
(1)散射单元距离及其传播时延计算
根据公式(14)，计算坐标为(-Lcosδ，0，Lsinδ)的发射点T到散射单元Sik中心点的距离RTik。
根据公式(15)，计算Sik到坐标为(Lcosδ,0,-Lsinδ)的接收点R的距离RRik。
根据公式(16)，计算由发射点T到海底散射单元Sik，再到接收点R的传播时间τik。
3.4计算散射单元的入射角、散射角
按照公式(18)、(19)计算，T点入射到散射单元Sik的垂直入射角αik和散射单元Sik到接收点R的垂直散射角βik。
3.5计算散射单元的散射强度
(1)计算散射单元Sik所在区域的面积
根据步骤3.1计算得到的散射单元格边界坐标(bxik,byik)，用公式(20)-(21)计算各个散射单元的面积。
(2)计算第i个散射环第k个散射单元Sik的海底散射强度
根据步骤3.4计算αik和βik，按照公式(26)计算散射单元Sik的散射系数Sbik。
(3)根据散射单元面积Abik和界面散射强度Sbik，按公式(25)计算散射单元的散射强度TSbik。
3.6计算散射单元Sik的多普勒频移
根据公式(22)计算散射点Sik处散射回波的归一化双基地多普勒频率fdik。
3.7散射单元的散射损失计算
(1)根据中心频率f，按照公式(27)计算，海水的声吸收系数α(f)。
(2)根据步骤3.4得到的散射单元距离RTik和RRik，以及步骤3.5得到的TSbik，按照公式(24)计算散射单元Sik以分贝为单位的散射损失TLrik。
假如RTik满足(17)式，则该散射单元、纳入混响计算，否则不计算。
(3)再由公式(28)换算为散射损失的声压幅度系数SBik。
3.8单个散射单元的混响信号计算
(1)按照(29)式计算单个散射单元的混响信号
(2)按照(30)式计算第i个散射环的混响信号
(3)按照(31)式计算t时刻的海底混响信号
3.9判断混响计算是否满足结束条件
根据(32)判断是否满足结束条件。
如果满足结束条件，混响计算结束；否则重复3.1～3.9，继续计算混响信号。
4设置z＝H0，按照第2～3步骤(3.1～3.9)计算海面混响信号。
计算海面混响信号时相应的公式，用海面散射相关公式替换。
5按照相同时间采样点，叠加海面和海底混响信号，形成总的海洋界面混响信号。
6计算结果图：
图1-1和图1-2分别是在海深200米，考虑“无声散射圆”和没有考虑“无声散射圆”的界面混响结果。因为是浅海，海底混响为主导，两种情况下的混响到达时间、混响强度差别不大，比较接近。
图1-3和图1-4分别是在海深5000米，考虑“无声散射圆”和没有考虑“无声散射圆”的界面混响结果。这是深海，先到达的是海面混响、后到达的是海底混响。这两种情况下，海面混响到达时间、混响强度基本相同；海底混响的到达时间和混响强度，差别明显。考虑“无声散射圆”的情况，更加符合仿真计算的实际。
图1-5和图1-6分别是海深200米时的海底、海面散射椭圆环、散射中心位置点、散射区域边界点的坐标位置图。
图1-7和图1-8分别是归一化界面混响信号瞬时幅值直方图和高斯分布对比图、归一化界面混响信号包络直方图和瑞利分布对比图。从两幅图可以看出，由本专利计算得到的界面混响信号瞬时幅值符合高斯分布、混响信号包络符合瑞利分布，满足混响信号的统计特性。
有益效果：
1、本发明建立了收发分置单基地声纳海底、海面散射的最小散射椭球确定方法和后续散射椭球的确定方法。最小散射椭球确定方法解决了收发分置单基地声纳最先达到、后续到达界面散射信号的确定方法问题；后续散射椭球的确定方法可以有限解决长脉冲主动信号混响信号的计算问题，使得产生的混响信号特征更加符合混响的统计特性。两者综合的优点是提高混响信号到达时间精度、多种混响信号叠加的强度精度，以及混响信号的真实性。该发明同样适用于多基地声纳界面混响建模仿真。
2、本发明提出了一种基于界面散射椭圆环的散射单元划分方法，该方法针对水平全向发射主动声纳，将整个水平面分成K个方位区域，这样在每个界面散射椭圆环上就形成了K个散射单元。区域数量K可根据接收阵的波束宽度、界面混响信号的多普勒频移分辨率，以及混响计算的时间代价等综合因素确定。
3、本发明建立了散射单元边界点、散射单元中心点坐标的计算模型，为计算散射单元面积、散射单元多普勒频移等奠定基础。
4、本发明建立散射单元发射入射角、接收散射角、到收-发点的距离、混响时延等计算模型，为各个散射单元的散射损失计算和混响信号计算奠定基础。
5、本发明建立海底、海面混响信号的计算模型。(1)提出了一种基于散射单元中心的收发分置的多普勒频移计算模型，给出了频域和时域表达式，可分别用于窄带信号和宽度信号；(2)给出了收发分置海面散射强度计算公式，可为各种收发分置的单基地或多基地声纳海面混响建模提供支撑。
6、本发明更新了界面混响计算结束条件。本发明考虑到主动声纳接收机的阵处理能力，增加了阵列处理的指向性DI因素，使得界面混响与海洋混响衔接更加平滑。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (10)

1.一种收发分置单基地主动声纳界面混响信号仿真方法，其特征在于，所述方法包括：
建立收发分置单基地声纳界面散射椭球模型；
基于所述散射椭球模型确定海底/海面散射的最小散射椭球，并根据所述最小散射椭球得到后续的散射椭球；
将每个所述散射椭球按方位进行划分，形成距离由散射椭圆环确定、方向由方位区域确定的若干散射单元，并计算各所述散射单元的边界点和中心点坐标；
根据所述散射单元的边界点和中心点坐标，计算散射单元的发射入射角、散射角、接收散射角、混响时延，以及各散射单元到收-发点的距离和所在区域面积；
建立多普勒频移模型，计算所述散射单元处散射回波的归一化双基地多普勒频移、时域多普勒因子；
根据所述散射单元的发射入射角、散射角、接收散射角、混响时延、到收-发点的距离及所在区域面积，计算所述散射单元的散射损失；根据所述散射损失和所述多普勒频移、时域多普勒因子，计算所述散射单元的混响信号；
当计算得到的所述混响信号满足混响级低于海洋环境噪声级与基阵指向性DI之差值时，所述散射单元的混响仿真程序结束。
2.根据权利要求1所述的收发分置单基地主动声纳界面混响信号仿真方法，其特征在于，所述建立收发分置单基地声纳界面散射椭球模型，包括：
设X轴沿拖曳平台的运动方向为正，发射阵在前、接收阵在后，z轴朝水面方向为正、海底方向为负；T、R分别为发射阵和接收阵的位置，发射阵与接收阵深度分别为hT和hR，收-发站间的斜距离为2L，其连线与水平面的倾斜角为δ。
3.根据权利要求2所述的收发分置单基地主动声纳界面混响信号仿真方法，其特征在于，所述基于所述散射椭球模型确定海底/海面散射的最小散射椭球，包括：
基于发射垂直波束为±Ψ/2的限定条件下，对应收发分置单基地声纳，存在一海底/海面无声散射圆，所述无声散射圆是一个以T′为圆心，以Ly为半径的圆，海底无声散射圆的半径表示为：
Ly＝(H-hT)tan(90-Ψ/2)＝(H-hT)ctan(Ψ/2) (1-1)
海面无声散射圆的半径表示为：
Ly＝hTtan(90-Ψ/2)＝hTctan(Ψ/2) (1-2)
以Ly为半径的无声散射圆与T'-R'连线内的交点作为散射点；
发射声到达海底散射点的最小斜距RTmin可表示为：
Figure FDA0003877683200000021
发射声到达海面散射点的最小斜距RTmin可表示为：
Figure FDA0003877683200000022
海底散射点到达接收点的斜距RRmin可表示为：
Figure FDA0003877683200000023
海面散射点到达接收点的斜距RRmin可表示为：
Figure FDA0003877683200000024
其中，2L'为发射点T与接收点R之间的水平距离，即T′、R′之间的距离；
最先到达接收点的的海底/海面散射椭球为最小散射椭球，其半长轴表示为：
Figure FDA0003877683200000025
对应的两个半短轴为：
Figure FDA0003877683200000026
式(4)和式(5)中，a0、b0、c0定义了所述最小散射椭球。
4.根据权利要求3所述的收发分置单基地主动声纳界面混响信号仿真方法，其特征在于，所述根据所述最小散射椭球得到后续的散射椭球，包括：
所述后续的散射椭球中第1个散射椭球外环的半长轴和半短轴分别为：
Figure FDA0003877683200000027
第i＝1,2,....个散射椭球外环的半长轴和半短轴分别为：
Figure FDA0003877683200000028
5.根据权利要求4所述的收发分置单基地主动声纳界面混响信号仿真方法，其特征在于，所述将每个所述散射椭球按方位进行划分，形成距离由散射椭圆环确定、方向由方位区域确定的若干散射单元，并计算各所述散射单元的边界点和中心点坐标，包括：
所述散射椭球按方位划分的方法是：若某一方位范围内的多普勒偏移量的变化在允许范围内，则该方位区域就确定为一独立散射单元格。
设第i个散射椭圆的单元格边界点坐标为(bxi,byi)，计算各所述散射单元的边界点坐标的方法如下：
散射椭球在海底/海面的截面为xoy平面上的椭圆，椭圆方程为：
Figure FDA0003877683200000031
式中，对于海底混响计算时，通过第i＝1,2,....个散射椭球在海底z＝-(H-H0)的截面来计算；对于海面混响计算时，通过第i＝1,2,....个散射椭球在海底z＝H0的截面来计算；
将每个散射椭圆划分成K个区域，则第i个散射椭圆外环上、第k个散射单元边界点Bik的极角为ρbk，第k个散射单元的坐标为bxik,byik，通过式(9)-(10)来计算：
bxik＝byikcot(ρbk)或byik＝bxiktan(ρbk) (9)
Figure FDA0003877683200000032
计算各所述散射单元的中心点坐标的方法如下：
设Sik为第i个散射椭圆环上、第k个散射单元的散射中心位置，其坐标用(sxik,syki,z)表示，第i＝1,2,....个散射椭圆环中心环的半长轴和半短轴分别为：
Figure FDA0003877683200000033
设第k个散射单元中心的极角为ρk，则Sik的坐标为：
Figure FDA0003877683200000034
或
Figure FDA0003877683200000035
上式(12)的两个公式之一与椭圆方程联合，即可求解所述散射单元的中心点坐标(sxik,syik)，表示为：
Figure FDA0003877683200000036
6.根据权利要求5所述的收发分置单基地主动声纳界面混响信号仿真方法，其特征在于，所述根据所述散射单元的边界点和中心点坐标，计算散射单元的发射入射角、散射角、接收散射角、混响时延，以及各散射单元到收-发点的距离和所在区域面积，并计算所述散射单元处散射回波的归一化双基地多普勒频移，包括：
(1)建立散射单元距离及其传播时延计算模型：
坐标为(-Lcosδ，0，Lsinδ)的发射点T到散射单元Sik中心点的距离为：
Figure FDA0003877683200000041
Sik到坐标为(Lcosδ,0,-Lsinδ)的接收点R的距离为：
Figure FDA0003877683200000042
由发射点T到海底散射单元Sik，再到接收点R的传播时间为：
Figure FDA0003877683200000043
由于垂直发射波束的限制，在发射点正下方海底存在声散射影区，该区域声波不能到达，这个区域不发生散射，在后面散射椭球扩展计算散射强度的过程中，不应计入，对应平坦海底，基本判据是：
RTik>RTmin (17)
(2)散射单元的入射角、散射角计算模型：
T点入射到散射单元Sik的垂直入射角为：
海底散射：
Figure FDA0003877683200000044
海面散射：
Figure FDA0003877683200000045
散射单元Sik到接收点R的垂直散射角为：
海底散射：
Figure FDA0003877683200000046
海面散射：
Figure FDA0003877683200000047
(3)各个散射单元所在区域面积的计算模型：
根据计算得到的第i个椭圆环、第k个散射单元的边界点坐标为(bxik,byik)，进一步计算第i个椭圆环、第k个扇区的面积：
Figure FDA0003877683200000048
则散射单元Sik的面积为：
Figure FDA0003877683200000051
计算扇区面积时，在同一椭圆按照下标号递增进行，考虑区域的对称性，只需要计算1/4个扇区、面积相同区域直接采用即可。
7.根据权利要求6所述的收发分置单基地主动声纳界面混响信号仿真方法，其特征在于，所述建立多普勒频移模型，计算所述散射单元处散射回波的归一化双基地多普勒频移、时域多普勒因子，包括：
根据同一拖曳平台运动情况下，发射站速度、接收站运动速度与平台运动速度v相同，则散射点Sik处散射回波的归一化双基地多普勒频率表示为：
Figure FDA0003877683200000052
时域多普勒因子表示为：
Figure FDA0003877683200000053
式中，αik为T到达散射单元Sik的掠射角、βik为散射单元Sik到R的掠射角；θTik为T'到达散射单元Sik的舷角；θRik为R'到达散射单元Sik的舷角。
8.根据权利要求7所述的收发分置单基地主动声纳界面混响信号仿真方法，其特征在于，所述根据所述散射单元的发射入射角、散射角、接收散射角、混响时延、到收-发点的距离及所在区域面积，计算所述散射单元的散射损失，包括：
在深海，经海底/海面散射单元Sik的总散射损失为：
TLrik(t)＝TSbik(t)-20lgRTik(t)-20lgRRik(t)-10-3α(f)(RTik(t)+RRik(t)) (24-1)
在浅海，扩展损失服从1.5次方规律，即介于平面波和球面波之间，经海底/海面散射单元Sik的总散射损失为：
TLrik(t)＝TSbik(t)-15lgRTik(t)-15lgRRik(t)-10-3α(f)(RTik(t)+RRik(t)) (24-2)
其中，TSbik为散射单元Sik的海底散射强度，表示为：
TSbik＝Sbik+10lgAbik (25-1)
对于海面散射，式(24)中散射单元Sik的TSbik由海面散射强度TSsik代替：
TSsik＝Ssik+10lgAsik (25-2)
式(27)中的Abik和Asik为第i个散射环第k个海底、海面散射单元的面积，Sbik、Ssik分别为海底散射强度、海面散射系数，由(26)式表示：
Sbik＝10lg(μsinαiksinβik) (26-1)
Figure FDA0003877683200000061
式中，β＝158(vf1/3)-0.58，αik是该散射单元入射波的掠射角，βik是该散射单元散射波的掠射角，v表示海面风速，f表示中心频率；
α(f)为海水的声吸收系数，对于中心频率利用Thorp经验公式表示为：
Figure FDA0003877683200000062
散射单元Sik散射信号的散射损失为：
Figure FDA0003877683200000063
9.根据权利要求8所述的收发分置单基地主动声纳界面混响信号仿真方法，其特征在于，所述根据所述散射损失和所述多普勒频移、时域多普勒因子，计算所述散射单元的混响信号，包括：
假设发射信号为声源级为SL的信号s(t)，不考虑发射水平指向性，发射垂直指向性按照前面的限定，并考虑散射单元的多普勒影响，T点发射、通过海底/海面散射单元Sik到达接收点R的散射信号可表示为：
Figure FDA0003877683200000064
或：
Figure FDA0003877683200000065
式中，fdik为多普勒频移，Δik为对应的多普勒因子(时域表示)，φik为服从均匀分布的随机相位，在0-2π间均匀分布，τik为信号发射到第Sik散射单元、再返回接收点总声程的往返时延，按式(16)计算；
第i个散射椭圆环的散射信号计算模型，表示为：
Figure FDA0003877683200000066
ABi为第i个散射环产生的散射信号的随机幅度，服从正态分布；
海底/海面散射的混响信号计算模型，表示为：
海底混响信号：
Figure FDA0003877683200000071
海面混响信号：
Figure FDA0003877683200000072
Nb(t)、Ns(t)分别为t时刻海底散射信号、海面散射信号到达接收端的散射椭圆环的总数；
将式(31-1)和(31-2)同时刻的信号相加，即得到海底和海面共同产生的总的混响信号：
y(t)＝yb(t)+ys(t) (32)。
10.根据权利要求9所述的收发分置单基地主动声纳界面混响信号仿真方法，其特征在于，所述当计算得到的所述混响信号满足混响级低于海洋环境噪声级的条件时，所述散射单元的混响仿真程序结束，包括：
考虑到主动声纳接收机的阵处理能力，增加阵列处理的指向性DI增益，以使得界面混响与海洋混响衔接更加平滑，当满足式(33)的条件时，所述散射单元的混响仿真程序结束；
Figure FDA0003877683200000073
CN202211220364.1A 2022-10-08 2022-10-08 一种收发分置单基地主动声纳界面混响信号仿真方法 Pending CN115598627A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211220364.1A CN115598627A (zh) 2022-10-08 2022-10-08 一种收发分置单基地主动声纳界面混响信号仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211220364.1A CN115598627A (zh) 2022-10-08 2022-10-08 一种收发分置单基地主动声纳界面混响信号仿真方法

Publications (1)

Publication Number Publication Date
CN115598627A true CN115598627A (zh) 2023-01-13

Family

ID=84845071

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211220364.1A Pending CN115598627A (zh) 2022-10-08 2022-10-08 一种收发分置单基地主动声纳界面混响信号仿真方法

Country Status (1)

Country Link
CN (1) CN115598627A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116299381A (zh) * 2023-05-19 2023-06-23 中国海洋大学 叠加海底山杂波耦合效应的深海混响快速预报方法及系统

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116299381A (zh) * 2023-05-19 2023-06-23 中国海洋大学 叠加海底山杂波耦合效应的深海混响快速预报方法及系统
CN116299381B (zh) * 2023-05-19 2023-08-15 中国海洋大学 叠加海底山杂波耦合效应的深海混响快速预报方法及系统

Similar Documents

Publication Publication Date Title
Denbigh Swath bathymetry: Principles of operation and an analysis of errors
US8107320B2 (en) Autonomous sonar system and method
EP3449279A1 (en) Acoustic doppler system and method
CN108107436A (zh) 一种基于可靠声路径的水下目标主动分类与定位方法
CN109061654B (zh) 一种深海环境下单圆环阵主动三维定位方法
CN110456361A (zh) 大范围海底地貌遥测的双基地声学成像方法
US7738319B2 (en) Determining angles of arrival using multipaths
US8400875B2 (en) Active sonar system and active sonar method using a pulse sorting transform
CN115598627A (zh) 一种收发分置单基地主动声纳界面混响信号仿真方法
Farmer et al. Intermediate range fish detection with a 12-kHz sidescan sonar
RU2700278C1 (ru) Способ определения местоположения подводного объекта
CN207689518U (zh) 三维空间来风的风速测量系统
CA2774758A1 (en) Method and device for measuring a profile of the ground
CN110441779B (zh) 一种多声纳浮标分布式协同定位方法
Qin et al. The 3D imaging for underwater objects using SAS processing based on sparse planar array
Ainslie et al. Echo, reverberation, and echo-to-reverberation ratio for a short pulse in a range-dependent Pekeris waveguide
CN110865359A (zh) 一种基于接收信号强度的水声测距方法
Abraham et al. Sonar Systems and the Sonar Equation
CN111650159B (zh) 一种海面反向散射强度测量方法
Keary et al. Simulation of the correlation velocity log using a computer based acoustic model
Xue et al. Integrated design of velocity measurement sonar system based on ACL and phased array ADL
Alexandrou Computer simulations of bathymetric sidescan sonar
Xiang-jian et al. An Imaging Algorithm for High-speed Side-scan Sonar Based on Multi-beam Forming Technology
CN118011372A (zh) 一种深水mbes海底回波信号仿真方法及系统
Dong et al. Modeling of Space-Time Reverberation under High-speed Vertical Gradient Motion

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