CN1996047A - 单接收阵合成孔径声纳干涉测量方法 - Google Patents

单接收阵合成孔径声纳干涉测量方法 Download PDF

Info

Publication number
CN1996047A
CN1996047A CNA2006101253421A CN200610125342A CN1996047A CN 1996047 A CN1996047 A CN 1996047A CN A2006101253421 A CNA2006101253421 A CN A2006101253421A CN 200610125342 A CN200610125342 A CN 200610125342A CN 1996047 A CN1996047 A CN 1996047A
Authority
CN
China
Prior art keywords
signal
emission
synthetic
echo
battle 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.)
Granted
Application number
CNA2006101253421A
Other languages
English (en)
Other versions
CN100529794C (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 CNB2006101253421A priority Critical patent/CN100529794C/zh
Publication of CN1996047A publication Critical patent/CN1996047A/zh
Application granted granted Critical
Publication of CN100529794C publication Critical patent/CN100529794C/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/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8902Side-looking sonar
    • G01S15/8904Side-looking sonar using synthetic aperture techniques

Landscapes

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

Abstract

本发明涉及一种单接收阵合成孔径声纳干涉测量方法,该方法采用两个或两个以上发射阵发射同频的正交信号,采用一个接收阵接收信号回波获取合成孔径成像干涉处理所需信息,从而有效降低系统的复杂程度,控制设备的成本。

Description

单接收阵合成孔径声纳干涉测量方法
技术领域
本发明属于合成孔径声纳干涉领域,提出了一种全新概念的单接收阵合成孔径声纳干涉测量方法。
背景技术
在合成孔径声纳成像的应用中,为了实现测深的目的,需要用垂直于平台运动方向的两条或者多条收发射阵进行干涉处理。以往的方法是采用一条发射基阵,两条或多条接收基阵,在同一时刻获取两幅或多幅合成孔径声纳成像数据,然后对图像进行干涉处理,得到目标的高程数据。
这种方法原理如图1所示,信号处理时将来自两个不同接收阵的信号分别进行合成孔径成像,然后再对它们进行干涉处理,就可以得到一幅目标的高程数据了。从目前的应用情况来看,实现这一技术需要克服以下几个方面的要求。
首先,在合成孔径声纳应用中,为了克服测距模糊和空间采样率不足之间的矛盾,往往采用多子阵接收的方法(如图2所示)。因此实用的合成孔径声纳,接收阵长度一般比发射阵长很多,成本也要高很多。将这一技术应用到干涉上,多出的一个或多个接收阵会进一步增加成本。
其次,由于干涉仪合成孔径声纳接收阵元的数目大幅度增加,声纳接收机的复杂程度也会增加,信号采集设备的通道数目也将成倍增长,同时也给数据存储和数据传输提出了比普通合成孔径声纳更高的要求。
此外,干涉合成孔径声纳对参加干涉的两个或多个接收阵各通道的幅度和相位一致性提出了很高的要求。
发明内容
本发明的目的是为了解决上述背景技术存在的不足,提出一种利用一个接收阵接收回波信号来实现单基线或多基线的单接收阵合成孔径声纳干涉测量方法,从而有效降低系统的复杂程度,控制设备的成本。
本发明采用两个或两个以上发射阵发射同频的正交信号,采用一个接收阵接收信号回波来获取合成孔径成像干涉处理所需信息。
上述每个发射阵之间可以具有不同的相位中心,发射阵最好沿垂直于平台运动的方向排列。
若信号Si(t)和Sj(t)满足如下的关系式:
∫ T S i ( t ) S j ( t - τ ) dt = δ ( τ ) , i = j ∫ T S i ( t ) S j ( t - τ ) dt = 0 , i ≠ j 则称Si(t)i=1,2,…n为一组正交信号。其中,δ(τ)是冲击函数。
理想的同频正交信号是指频率相同、互相关函数等于零的信号。但实际应用中难以找到这样理想的正交信号。在本发明中只要其互相关函数的最大值远小于自相关函数最大值的信号都可以视为正交信号。
本发明步骤为:
1.用上述发射阵同时发射相互正交的信号。发射阵发射的信号为Ti,经目标反射后回波信号分别为Sii)。式中τi是由于收发路径不同而产生的不同延迟。
2、用一个接收阵接收这些回波,并对信号回波进行脉冲压缩,同时将来自于不同发射阵的信号回波分离;
3、分别对每个分离后的回波信号进行合成孔径成像处理,得到Fi(T,R),式中,T是合成孔径成像声纳缓变时间坐标,R是声纳作用方向的坐标。
4、将合成孔径成像的结果进行干涉处理,得到目标区域的高度信息。
上述脉冲压缩采用匹配滤波的方法,公式为:Scomii)=∫S(t)×Ti *i-t)dt,其中,S(t)是接收阵接收到的所有发射信号的回波之和,Ti(t)是第i个发射阵发射的信号,τi是由于收发路径不同而产生的不同延迟,Scomii)是脉冲压缩后被分离的信号回波。
如果该方法采用两个发射阵发射同频的正交信号,所述干涉处理是采用上述两个发射阵的回波合成孔径信号F1(T,R),F2(T,R)进行干涉、相位解卷和相位高度变换处理,得到基于单基线的目标区域高度信息。
如果该方法采用三个或三个以上发射阵发射同频的正交信号,所述干涉处理是采用上述三个或三个以上发射阵的回波合成孔径信号Fi1(T,R),Fi2(T,R),…,Fik(T,R)进行干涉、相位解卷和相位高度变换处理,得到基于多基线的目标区域高度信息。
综上所述,本发明利用多个发射阵和一个接收阵的收发阵结构,通过对正交回波信号进行分离和分别成像处理得到多幅来自不同收发路径的合成孔径图像,将这些合成孔径图像进行干涉处理,实现了干涉仪合成孔径声纳。这种方法相对于以往的采用多接收阵的干涉仪合成孔径技术主要有以下优点
一、相对于以往的采用多接收阵的干涉仪合成孔径技术,本发明大幅度地削减了声纳接收阵的成本。一般而言,同样长度的声纳接收阵和声纳发射阵需要花费的成本大致相等,而在多子阵合成孔径声纳技术投入应用之后,合成孔径声纳的接收阵的长度已经远远大于发射阵的长度,在成本上也远高于发射阵。按照以往的做法,如果要实现干涉或者多基线条件下的干涉测量,就有必要装备两条或数条成本较高的接收阵,而如果采用本发明中的办法,我们只需要装配多条成本比较低的发射阵,和一条成本较高的接收阵,这样整体系统的成本也大幅度降低了。
二、相对于以往的采用多接收阵的干涉仪合成孔径技术,本发明有效减少了接收机的通道数,降低接收系统的复杂程度。由于本发明只采用一条接收阵,多个正交信号同时接收,显然接收机的通道数会成倍降低。
三、本发明降低了干涉仪合成孔径声纳对原始数据采集、存储和传输的要求,很大程度上降低了数据采集和存储系统的成本。干涉仪合成孔径声纳对于数据采集、存储和传输的需求非常高,往往几分钟的声纳扫描,就将产生几千兆甚至更多的原始数据,因此数据采集设备和存储设备的成本都非常高。采用单接收阵的干涉仪合成孔径声纳,数据量也将成倍减少,对数据采集和存储设备的技术指标要求也会相应降低。
附图说明
图1干涉合成孔径成像示意图;
图2多子阵接收方法示意图;
图3两发射阵发射两种正交波形T1和T2;
图4基于双发射阵的干涉原理图;
图5N发射阵分别发射N种正交波形T1,T2,…,TN
图6基于N发射阵的干涉原理图。
具体实施方式
实施例1:参见图3,本实施例的声纳阵列由2个发射阵和1个接收阵组成。发射阵分别发送正交信号T1和T2,经目标反射后,接收阵元同时收到回波信号S11)和S22),其中τ1,τ2分别表示经过不同收发路径到达接收阵的信号延迟时间。利用发射信号的正交特性,在进行脉冲压缩时将来自于不同发射阵的信号回波分离,分别对分离后的回波信号进行合成孔径成像处理得到F1(T,R),F2(T,R)。将两个合成孔径成像的结果进行干涉处理,得到两个信号的相位差主值(T,R)=arg(F1(T,R)×conj(F2(T,R))),再通过二维相位解卷,得到两个信号的真实相位差φ(T,R)=(T,R)+2kπ。附图4中绘出了基于双发射阵的干涉测量原理:显然,二维相位解卷后的相位差是由信号延迟时间的差异决定的。
φ = 4 π λ [ ρ 1 - ρ 2 ] - - - ( 1 )
其中
ρ 2 = B 2 + ρ 1 2 - 2 Bρ 1 cos β = B 2 + ρ 1 2 - 2 Bρ 1 sin ( θ - α ) - - - ( 2 )
基于平面波(远场)假设,(1)式可简写成:
φ = 4 π λ B sin ( θ - α ) - - - ( 3 )
定义垂直基线:
B=Bcos(θ-α)    (4)
由几何关系可得:
cos θ = H - h ρ 1 - - - ( 5 )
利用(3)、(4)两式即可得到场景的高度:
h=H-ρ1 cosθ    (6)
实施例2:参见图5,本实施例的声纳阵列由N个发射阵(3个或3个以上)和1个接收阵组成。参与干涉的任意两对收发射阵相位中心之间的距离又称为基线,包含N个发射阵和1个接收阵的干涉合成孔经声纳系统可以有
Figure A20061012534200081
条基线。附图6中绘出了上述情况下基线的排列情况,基线分别为B(1,2),B(2,3),B(1,3),...,B(1,N)。实施例中的基线长度满足B(1,2)<B(2,3)<B(1,3)<,...,<B(1,N)=B0
发射阵分别发送正交信号T1,T2,…Tn,经目标反射后,接收阵元同时收到回波信号S11),S22),…,Snn),其中τ1,τ2,…τn分别表示经过不同收发路径到达接收阵的信号延迟时间。利用发射信号的正交特性,在进行脉冲压缩时将来自于不同发射阵的信号回波分离,分别对分离后的回波信号进行合成孔径成像处理得到F1(T,R),F2(T,R),…,Fn(T,R)。第i和第j个合成孔径成像的结果的干涉相位为i,j(T,R),第1和第N个合成孔径成像的结果的干涉相位为0。根据干涉相位和基线的近似关系有 多基线干涉可以采用最大似然法估计干涉相位。利用多个合成孔径成像结果构造合适的似然函数fMI(0),求解偏微分方程 可以得到0的估计值。与实施例1中的单基线干涉处理相比,①在低信噪比条件下,多基线干涉得到的干涉图具有更高的干涉相位估计精度。②多基线干涉的干涉相位图具有更大的模糊基,相位解卷过程可以获得简化。
多基线干涉处理在相位——高度变换的过程中可以采用迭代法进行处理,先从最小的基线的干涉相位图出发,估计地形高度,然后以此为基准解卷较大基线的干涉相位图,依此类推直至获得最大基线的地形高度图。这种方法既利用了小基线相位解卷较好的优点,又获得了大基线在高度计算时的高精度。

Claims (6)

1、单接收阵合成孔径声纳干涉测量方法,其特征是该方法采用两个或两个以上发射阵发射同频的正交信号,采用一个接收阵接收信号回波获取合成孔径成像干涉处理所需信息。
2、根据权利要求1所述的单接收阵合成孔径声纳干涉测量方法,其特征是所述发射阵沿垂直于平台运动的方向排列。
3、根据权利要求1或2所述的单接收阵合成孔径声纳干涉测量方法,其特征是所述接收阵在接收信号回波后,首先对信号回波进行脉冲压缩,同时将来自于不同发射阵的信号回波分离;然后分别对每个分离后的回波信号进行合成孔径成像处理;最后将合成孔径成像的结果进行干涉处理,得到目标区域的高度信息。
4、根据权利要求3所述的单接收阵合成孔径声纳干涉测量方法,其特征是所述脉冲压缩采用匹配滤波的方法,公式为: S comi ( τ i ) = ∫ S ( t ) × T i * ( τ i - t ) dt , 其中,S(t)是接收阵接收到的所有发射信号的回波之和,Ti(t)是第i个发射阵发射的信号,τ,是由于收发路径不同而产生的不同延迟,Scomii)是脉冲压缩后被分离的信号回波。
5、根据权利要求3所述的单接收阵合成孔径声纳干涉测量方法,其特征是该方法采用两个发射阵发射同频的正交信号,所述干涉处理是采用上述两个发射阵的回波合成孔径信号F1(T,R),F2(T,R)进行干涉、相位解卷和相位高度变换处理,得到基于单基线的目标区域高度信息。
6、根据权利要求3所述的单接收阵合成孔径声纳干涉测量方法,其特征是该方法采用三个或三个以上发射阵发射同频的正交信号,所述干涉处理是采用上述三个或三个以上发射阵的回波合成孔径信号Fi1(T,R),Fi2(T,R),…,Fik(T,R)进行干涉、相位解卷和相位高度变换处理,得到基于多基线的目标区域高度信息。
CNB2006101253421A 2006-12-07 2006-12-07 单接收阵合成孔径声纳干涉测量方法 Active CN100529794C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2006101253421A CN100529794C (zh) 2006-12-07 2006-12-07 单接收阵合成孔径声纳干涉测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2006101253421A CN100529794C (zh) 2006-12-07 2006-12-07 单接收阵合成孔径声纳干涉测量方法

Publications (2)

Publication Number Publication Date
CN1996047A true CN1996047A (zh) 2007-07-11
CN100529794C CN100529794C (zh) 2009-08-19

Family

ID=38251189

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2006101253421A Active CN100529794C (zh) 2006-12-07 2006-12-07 单接收阵合成孔径声纳干涉测量方法

Country Status (1)

Country Link
CN (1) CN100529794C (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101825709A (zh) * 2009-12-08 2010-09-08 中国科学院声学研究所 一种水下高分辨率侧视声学成像系统与方法
CN103149274A (zh) * 2013-01-28 2013-06-12 中国科学院声学研究所 一种混凝土缺陷检测方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101825709A (zh) * 2009-12-08 2010-09-08 中国科学院声学研究所 一种水下高分辨率侧视声学成像系统与方法
CN101825709B (zh) * 2009-12-08 2012-07-25 中国科学院声学研究所 一种水下高分辨率侧视声学成像方法
CN103149274A (zh) * 2013-01-28 2013-06-12 中国科学院声学研究所 一种混凝土缺陷检测方法
CN103149274B (zh) * 2013-01-28 2015-10-07 中国科学院声学研究所 一种混凝土缺陷检测方法

Also Published As

Publication number Publication date
CN100529794C (zh) 2009-08-19

Similar Documents

Publication Publication Date Title
CN107957574B (zh) 基于ifft和混合匹配追踪的时分地基mimo滑坡雷达成像方法
CN105699945B (zh) 频控阵mimo雷达系统中的波形优化设计方法
US9261592B2 (en) Method and system for through-the-wall imaging using compressive sensing and MIMO antenna arrays
CN102866401B (zh) 一种基于mimo技术的三维成像方法
CN110297233B (zh) Lfmcw阵列雷达信号并行流水化处理方法
CN106680815B (zh) 基于张量稀疏表示的mimo雷达成像方法
CN106291540A (zh) 一种基于doa估计的多输入多输出探地雷达逆向投影目标成像方法
CN109765562A (zh) 一种三维前视声像声纳系统和方法
CN102176008A (zh) 一种三维地层成像的相控方位滤波方法
CN105911544A (zh) 一种基于压缩感知技术的扩频穿墙雷达成像方法
CN103454630A (zh) 一种基于多阵元发射技术的超宽带三维成像方法
CN113253223B (zh) 基于步进频信号的非平稳杂波抑制的目标检测方法
CN103245943A (zh) 一种使用线性调频信号的mimo阵列扇扫成像方法
CN100419455C (zh) 合成孔径声纳运动误差和水声信道相位误差的测量方法
CN106125078A (zh) 一种水下多维声学成像系统和方法
CN104062663B (zh) 一种多波束海底浅地层剖面探测设备
CN100529794C (zh) 单接收阵合成孔径声纳干涉测量方法
CN103217671B (zh) 色噪声环境下的多输入多输出雷达收发角度快速估计方法
CN114814817A (zh) 运动目标速度解模糊方法、装置、电子设备及存储介质
CN103605112A (zh) 多发多收干涉合成孔径雷达时频二维信号波形设计方法
CN103926586B (zh) 一种使用发射子阵的mimo阵列测深方法
CN106093926A (zh) 双边带频率分集阵列雷达系统及其目标定位方法
CN110596707B (zh) 基于多快拍图像联合的mimo雷达三维成像方法
CN102590861B (zh) 基于接收阵列的地震波束形成方法
JP2020530562A (ja) センサアレイ撮像装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant