CN106918809B - 快速干涉合成孔径声纳原始回波时域仿真方法 - Google Patents

快速干涉合成孔径声纳原始回波时域仿真方法 Download PDF

Info

Publication number
CN106918809B
CN106918809B CN201710265617.XA CN201710265617A CN106918809B CN 106918809 B CN106918809 B CN 106918809B CN 201710265617 A CN201710265617 A CN 201710265617A CN 106918809 B CN106918809 B CN 106918809B
Authority
CN
China
Prior art keywords
distance
echo
synthetic aperture
time domain
array
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
CN201710265617.XA
Other languages
English (en)
Other versions
CN106918809A (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.)
Naval University of Engineering PLA
Original Assignee
Naval University of Engineering PLA
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 Naval University of Engineering PLA filed Critical Naval University of Engineering PLA
Priority to CN201710265617.XA priority Critical patent/CN106918809B/zh
Publication of CN106918809A publication Critical patent/CN106918809A/zh
Application granted granted Critical
Publication of CN106918809B publication Critical patent/CN106918809B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/52004—Means for monitoring or calibrating

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (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

快速干涉合成孔径声纳原始回波时域仿真方法
技术领域
本发明涉及干涉合成孔径声纳原始回波的仿真，具体地指一种快速干涉合成孔径声纳原始回波时域仿真方法。
背景技术
干涉合成孔径声纳(InSAS)是在合成孔径声纳基础上增加一幅(或多幅)接收基阵，通过比相测高的方法得到场景的高度信息，从而得到场景的三维图像。InSAS兼备了合成孔径声纳分辨率与成像距离和工作频率无关的优点和干涉测深精度高的优点。可以用于海底测量、水下考古、海底石油勘探和搜寻水下失落物体等，也能提高水雷的识别能力。
干涉合成孔径声纳信号处理包括合成孔径成像处理和干涉信号处理，信号处理过程非常复杂。为了验证各部分信号处理算法的正确性，在试验之前充分测试程序，必须开展原始回波的仿真方法研究。原始信号的计算过程十分复杂，特别是对于大场景的原始回波仿真，效率显得至关重要，研究快速的回波仿真方法意义重大。
目前干涉合成孔径声纳原始回波仿真方法主要有时域仿真和频域仿真两类。时域仿真完全按照干涉合成孔径声纳的工作过程进行仿真，是最精确的仿真方法，但是计算量也非常大，对计算机的性能要求很高。频域仿真方法效率比时域仿真要高，但是精度低。如何在保持仿真精度的前提下提高速度是一个技术难题。
发明内容
本发明的目的是为了解决现有技术中干涉合成孔径声纳原始回波时域仿真效率低的问题，从而提供一种快速的干涉合成孔径声纳原始回波时域仿真方法。为了实现以上目的，本发明采用如下方案：
根据仿真的场景数据、声基阵位置和姿态数据、系统运行的几何参数，计算每一反射点的反射系数和延迟距离，最后直接在时域内生成每一接收阵元的原始回波信号。
具体地，包括：
1)生成仿真所用的三维场景数据。将场景划分为小三角形面元，面元尺寸要大于声波波长，小于成像分辨率。计算出场景中每个三角形面元的反射系数和遮挡效应。对于场景的距离向、方位向和高度信息都加入随机偏移量。
2)根据声纳的工作参数载体速度、脉冲重复频率、基阵距底高度、基线倾角确定声基阵的空间位置。
3)参考信号二维数组的生成。发射信号的时宽乘以距离向采样频率作为二维数组的列数，将一个距离采样时间间隔T细分成NS等份，每份记为Ts，NS作为二维数组的行数，二维数组的每行起始时间相差Ts。
4)场景中每个面元的回波计算。按照基阵在方位向采样点的位置，计算发射阵到面元的距离和接收子阵到面元的距离，两者之和的一半作为延迟距离。延迟距离减去最近采样距离再除以距离向采样间隔的整数部分作为整数平移量，即回波信号中距离向采样起点，小数部分乘以NS后的整数部分作为参考信号在二维数组中的索引。
5)将场景中每个面元的回波进行叠加作为最终的回波结果。
6)采用同样方法分别对主、辅接收阵进行计算，获得原始干涉合成孔径声纳原始回波数据。
本发明方法能够有效避免冗余计算，大大提高原始回波计算效率，有利于合成孔径成像算法和干涉信号处理算法验证。
附图说明
图1为本发明快速干涉合成孔径声纳原始回波时域仿真方法的流程示意图。
图2为参考信号示意图。
图3为平移量计算示意图。
图4为平移量小数部分计算示意图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步的详细说明。
本发明的流程如附图1所示，下面结合附图和具体实施方式对本发明做进一步的详细说明。
(1)计算待仿真的三维场景数据。首先确定需要仿真的场景的大小，根据成像分辨率和波长的大小确定出面元尺寸，面元的大小要符合两个要求：第一，面元要大于声波的波长，第二，面元要小于成像分辨率。然后计算出整个场景在方位向和距离向上划分的面元数Ny，Nx。由于真实的场景并不是平坦的，有高低起伏，为了更真实的反应场景信息，在方位向和距离向以及场景高度上都加入随机量，随机量小于面元的尺寸。本实施例中，随机偏移量分为距离向随机偏移量、方位向随机偏移量和高度随机偏移量，其中距离向随机偏移量和方位向随机偏移量分别控制在半个距离向和方位向分辨率范围内，高度随机偏移量控制在半个波长范围内。
最后，将三角形的中点作为面元的中心，本实施例中所用三角形的中点为三角形三个顶点的坐标和的均值即为中点。
(2)计算场景中各个面元的反射系数。首先计算入射角，即发射阵与面元中心连线向量与面元的法向量的夹角，然后根据兰伯特定律计算出反向散射系数。最后，考虑遮挡效应，将位于阴影中的点处的反射系数设置为0。
(3)仿真基阵位置和姿态信息。加入姿态信息，可以使仿真的基阵位置更准确，更接近真实情况，因此本实施例需要考虑基阵的横滚、俯仰、横摇、纵荡、侧摆和升沉六个运动参数，以及x,y,z的空间位置信息。具体地，声基阵的位置根据载体速度、脉冲重复频率、基阵距底高度、基线倾角确定，声基阵姿态数据根据六自由度运动仿真参数确定。
(4)预先生成的参考信号是一个二维数组，如附图2所示。其行数NTp＝Tp·rFs，其中Tp是发射信号的时宽，rFs是距离向采样频率。参考信号的行数NS等于一个距离向采样间隔的划分数，即将一个距离向采样时间间隔T细分成1000等份，每份记为Ts，参考信号的每行相差一个Ts。参考信号的表达式为
Figure BDA0001275953130000041
其中K是线性调频斜率，
Figure BDA0001275953130000042
Figure BDA0001275953130000043
是参考信号列索引，ia＝0,1,…999是参考信号行索引。
(5)计算场景中每个面元的回波。根据基阵的位置计算发射阵到面元的距离r1和接收子阵到面元的距离r2，延迟距离
Figure BDA0001275953130000044
(6)计算当前延迟距离r对应的参考信号的平移量，平移量计算如附图3所示，计算平移量的整数部分
Figure BDA0001275953130000045
是最近采样距离，Δr是距离向一个采样间隔。NRstrat即对应该点回波信号的距离向采样点的起始点。
(7)计算平移量的小数部分，如附图4所示，将一个距离采样间隔Δr细分成NS等份，每份为
Figure BDA0001275953130000046
每份对应的时间记为Ts，计算从NRstrat到r处共有多少个Ts，计算式为
Figure BDA0001275953130000047
it即对应着二维参考信号中的行数。例如将距离向上一个时间采样间隔T分为1000等分，每份记为Ts，生成1000个发射参考信号，将其存放在一个二维数组中，参考信号每行的起始时间分别为0,1*Ts,…,999*Ts。时间采样间隔的选取与计算机资源有关，采样间隔分割数选取越大，仿真回波效果越好。
对于延迟距离为r的回波反射点，其回波的起始下标为(r-Rmin)/Δr的整数部分，其中Rmin为最近采样距离，回波在参考信号二维数组中的行数为(r-Rmin)/Δr的小数部分乘以NS(1000)后的整数部分。
(8)将参考信号的it行数据的开头平移到回波信号的NRstrat处，其结尾对应着回波信号的NRstrat+NTp处。
(9)对方位向每个脉冲和每个采样点都进行循环，计算回波的叠加，这样就完成了一个面元的仿真。
(10)对场景中的每个面元都进行上述过程，然后进行回波叠加，得到最终的仿真回波数据。
(11)对主、辅接收阵回波采用同样方式计算，获得干涉合成孔径声纳回波的原始数据数据。

Claims (8)

1.一种快速干涉合成孔径声纳原始回波时域仿真方法，其特征在于：根据仿真的场景数据、声基阵位置和姿态数据、系统运行的几何参数，计算每一反射点的反射系数和延迟距离，最后直接在时域内生成每一接收阵元的原始回波信号；计算任意一点原始回波时，将延时距离转换为距离向上的平移量；距离向上的平移量分为整部平移部分与小数平移部分；延迟距离减去最近采样距离然后再除以距离采样间隔得到的整数部分作为整数平移量，小数部分作为小数偏移量部分；对于延迟距离为r的回波反射点，其回波的起始下标为(r-Rmin)/Δr的整数部分，其中Rmin为最近采样距离，回波在参考信号二维数组中的行数为(r-Rmin)/Δr的小数部分乘以NS后的整数部分。
2.根据权利要求1所述快速干涉合成孔径声纳原始回波时域仿真方法，其特征在于：仿真场景的小面元的距离向和方位向间隔要小于合成孔径成像结果的距离向和方位向分辨率，大于波长，并且仿真场景中每一点的位置加入随机偏移量。
3.根据权利要求2所述快速干涉合成孔径声纳原始回波时域仿真方法，其特征在于：随机偏移量包括距离向随机偏移量、方位向随机偏移量和高度随机偏移量，其中距离向随机偏移量和方位向随机偏移量分别控制在半个距离向和方位向分辨率范围内，高度随机偏移量控制在半个波长范围内。
4.根据权利要求1所述快速干涉合成孔径声纳原始回波时域仿真方法，其特征在于：声基阵的位置根据载体速度、脉冲重复频率、基阵距底高度、基线倾角确定，声基阵姿态数据根据六自由度运动仿真参数确定。
5.根据权利要求1所述快速干涉合成孔径声纳原始回波时域仿真方法，其特征在于：系统运行几何参数需要满足仿真场景照射条件，将仿真场景置于声纳照射范围内。
6.根据权利要求1所述快速干涉合成孔径声纳原始回波时域仿真方法，其特征在于：首先计算入射角，即发射阵与面元中心连线向量与面元的法向量的夹角，然后根据兰伯特定律计算出反射系数；延迟距离的计算以仿真三角形面元的中心作为参考点，分别计算发射阵和接收阵到面元中心的距离，然后将两个距离之和作为总延迟距离。
7.根据权利要求6所述快速干涉合成孔径声纳原始回波时域仿真方法，其特征在于：小数部分取3位有效数据。
8.根据权利要求1所述快速干涉合成孔径声纳原始回波时域仿真方法，其特征在于：将一个距离采样间隔Δr细分成NS等份，每份为
Figure FDA0003916307690000021
每份对应的时间记为Ts，计算从NRstart到r处共有多少个Ts，计算式为
Figure FDA0003916307690000022
it即对应着二维参考信号中的行数。
CN201710265617.XA 2017-04-21 2017-04-21 快速干涉合成孔径声纳原始回波时域仿真方法 Active CN106918809B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710265617.XA CN106918809B (zh) 2017-04-21 2017-04-21 快速干涉合成孔径声纳原始回波时域仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710265617.XA CN106918809B (zh) 2017-04-21 2017-04-21 快速干涉合成孔径声纳原始回波时域仿真方法

Publications (2)

Publication Number Publication Date
CN106918809A CN106918809A (zh) 2017-07-04
CN106918809B true CN106918809B (zh) 2022-12-13

Family

ID=59567325

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710265617.XA Active CN106918809B (zh) 2017-04-21 2017-04-21 快速干涉合成孔径声纳原始回波时域仿真方法

Country Status (1)

Country Link
CN (1) CN106918809B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111142083B (zh) * 2020-02-27 2022-05-03 西北核技术研究院 一种短脉冲非相参雷达中频回波构建方法
CN112505667B (zh) * 2020-11-19 2022-07-15 哈尔滨工程大学 一种二维声呐基阵运动姿态自校准方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101598788A (zh) * 2008-06-03 2009-12-09 中国科学院声学研究所 合成孔径声纳信号的快速仿真方法
CN101907707A (zh) * 2010-07-02 2010-12-08 哈尔滨工程大学 用于多波束合成孔径声呐的组合声基阵
EP2660622A2 (en) * 2012-09-03 2013-11-06 Institute of Electronics, Chinese Academy of Sciences Imaging method and device in FMCW SAR
CN103593877A (zh) * 2013-11-07 2014-02-19 清华大学 合成孔径声纳图像的仿真方法及系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101598788A (zh) * 2008-06-03 2009-12-09 中国科学院声学研究所 合成孔径声纳信号的快速仿真方法
CN101907707A (zh) * 2010-07-02 2010-12-08 哈尔滨工程大学 用于多波束合成孔径声呐的组合声基阵
EP2660622A2 (en) * 2012-09-03 2013-11-06 Institute of Electronics, Chinese Academy of Sciences Imaging method and device in FMCW SAR
CN103593877A (zh) * 2013-11-07 2014-02-19 清华大学 合成孔径声纳图像的仿真方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
多子阵合成孔径声呐波数域算法不均匀采样问题研究;刘维等;《声学学报(中文版)》;20090515(第03期);全文 *

Also Published As

Publication number Publication date
CN106918809A (zh) 2017-07-04

Similar Documents

Publication Publication Date Title
CN106772365B (zh) 一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法
CN107976660B (zh) 弹载多通道雷达超低空目标分析与多径回波建模方法
AU2005201859B2 (en) Method of seismic source monitoring using modeled source signatures with calibration functions
CN112505666B (zh) 一种二维多波束合成孔径声呐目标回波仿真方法
EP3078991B1 (en) Method for swell effect and mis-tie correction in high-resolution seismic data using multi-beam echo sounder data
CN104101864B (zh) 基于eof分解的导航x波段雷达海浪参数反演算法
JP5950835B2 (ja) シーン内の反射物を再構成するシステム
CN105445718B (zh) 一种基于阵列重构的分布式多载舰超视距雷达的doa估计方法
CN106199601B (zh) 基于粗数字高程图的InSAR绝对相位模糊估计方法
CN108845325A (zh) 拖曳线列阵声纳子阵误差失配估计方法
Young et al. Resolving the internal and basal geometry of ice masses using imaging phase-sensitive radar
CN103760546A (zh) 一种雷达用低空目标波达方向估计方法
US11448755B2 (en) System and techniques for split-aperture beamforming
CN110907938B (zh) 一种近场的快速下视合成孔径三维成像方法
CN106918809B (zh) 快速干涉合成孔径声纳原始回波时域仿真方法
CN206546434U (zh) 一种水下多维声学成像系统
WO2009110096A1 (ja) 合成開口処理システム及び合成開口処理方法
CN102129068A (zh) 一种对合成孔径声纳系统相位误差估计的测试系统及方法
CN117848285A (zh) 海底地形测绘方法、设备、测绘系统及存储介质
CN110412585B (zh) 一种基于mvdr的下视合成孔径三维成像方法及系统
CN110850396B (zh) 一种应用于深海黑匣子搜探定位系统电模拟器及其轨迹生成方法
CN109884621B (zh) 雷达高度表回波相参积累方法
CN103984007A (zh) 定向地震波延时参数优化设计方法
JP5443891B2 (ja) 合成開口ソーナー
CN112902931B (zh) 无人船测深数据与定位数据之间延迟的测定和消除方法

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