CN109283511B - 一种宽覆盖多波束接收基阵校准方法 - Google Patents

一种宽覆盖多波束接收基阵校准方法 Download PDF

Info

Publication number
CN109283511B
CN109283511B CN201811017290.5A CN201811017290A CN109283511B CN 109283511 B CN109283511 B CN 109283511B CN 201811017290 A CN201811017290 A CN 201811017290A CN 109283511 B CN109283511 B CN 109283511B
Authority
CN
China
Prior art keywords
array
data
transducer
angle
rotation
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
CN201811017290.5A
Other languages
English (en)
Other versions
CN109283511A (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.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201811017290.5A priority Critical patent/CN109283511B/zh
Publication of CN109283511A publication Critical patent/CN109283511A/zh
Application granted granted Critical
Publication of CN109283511B publication Critical patent/CN109283511B/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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52004Means 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

一种宽覆盖多波束接收基阵校准方法
技术领域
本发明属于多波束接收基阵校准领域,具体涉及一种宽覆盖多波束接收基阵校准方法。
背景技术
围绕海洋的科学研究、资源开发、工程建设以及军事等活动中,通常都需要准确地获取所关注区域内的海底地形地貌信息作为基础资料与支撑依据。多波束测深声纳已成为国内外海洋科学研究、海底资源开发、海洋工程建设等海洋活动中最主要的海洋调查勘测仪器之一。
多波束测深声纳的发展逐渐朝着宽覆盖、高精度、一体化测量方向发展,所以对于接收换能器基阵的精度要求也逐步提高。对于宽覆盖多波束测深***,接收基阵低掠射波束角度下的测深精度逐渐降低,并且呈非线性变化,波束角度偏移严重。所以,为了达到宽覆盖情况下的高精度测量,对于接收基阵的测量校准就是必不可少的。
传统的换能器基阵测量方法只是单纯的将各接收阵元并联,信号无指向性叠加,测量组合基阵的自然指向性。这种测量方法并不能与多波束测深机理相吻合,多波束测深算法需要首先进行波束形成,因此更关心基阵在各波束角度上的指向性输出,自然指向性不能表征此特性。另外,传统方法要求基阵处于远场平面波假设下,各阵元信号才能够等效成无指向性直接叠加,当基元工作频率较低时远场距离较远,需要水池尺寸巨大。传统的基阵校准方法只是单纯的校正换能器特性,不考虑与之配套的信号接收电路特性,信号接收电路特性一般单独测量,这种测量方法对于通道数较少的***影响不大,而对于具有上百路接收通道,电路幅度、相位一致性各不相同的多波束测深***明显不适用。在进行指向性测量时,一般需要在大尺度水池下使换能器处于远场条件下,进行多自由度伸缩杆的旋转,每个角度的测量都需要旋转、暂停螺杆,这将会引入空程误差,对于大角度的指向性范围测量影响不大,但对于高精度要求的多波束声纳波束指向性测量会带来较大的角度偏差。
基于以上问题,本发明提出了一种宽覆盖多波束接收基阵校准方法,将多波束测深声纳接收机电路与换能器基阵作为整体测量,经***充分预热后排除电路温度漂移影响,使***处于日常工作状态下,不单独测量各分机特性。测量过程采取自动单向无停顿旋转方法,通过信号同步线控制数据采集,不需要严格找寻中心零点角度,通过后续算法处理校正初始测量位置角度,排除机械旋转空程误差。采用近场聚焦波束形成算法进行基阵校准,能够在较小尺寸的消声水池中使待测基阵处于近场条件下进行测量,计算得到基阵波束指向性曲线和波束角度误差曲线。本发明提出的方法能够精细的校准多波束声纳接收基阵波束指向性,有效地提升多波束测深声纳宽覆盖角度下的测深精度。
发明内容
本发明的目的是在于提供一种宽覆盖多波束接收基阵校准方法。
该方法包含以下步骤:
(1)根据接收换能器基阵工作频率、换能器尺寸选取相应尺寸的消声水池,通过声速仪记录下测量时刻水池的实时声速值,计算基阵布放位置和布放方位向;
(2)根据选取的合适尺寸的消声水池以及计算得到的基阵布放位置和布放方位向,将待测换能器基阵按照布放方位向安装在多自由度伸缩杆上,测量待测换能器基阵与水平方向夹角,基阵入水放置在计算好的布放位置处,将一个全指向性发射换能器放置在计算好的发射位置处,调节布放深度使其对准待测基阵等效声中心;
(3)根据调节好的布放深度,将发射换能器与数据控制存储单元通过同步信号线相连,打开信号源设置工作参数,功率放大器开机预热,打开数据控制存储单元,多波束测深***空采预热,使用示波器实时监测发射换能器发射波形与待测基阵接收波形;
(4)预热完成后将待测基阵旋转至垂直于发射换能器角度,设置多自由度伸缩杆旋转方向与旋转速度,开启多波束测深***采集软件,多自由度伸缩杆顺时针自动旋转的同时采集软件自动将对应位置的回波信号存储至数据控制存储单元,由数据控制存储单元按照秒脉冲产生同步信号,发射换能器和待测基阵采集***按照秒脉冲信号同步工作、记录数据;
(5)单方向测量完成后设置参数,再进行一次逆时针自动旋转测量,采集软件自动将对应位置的回波信号存储至数据控制存储单元;
(6)数据采集完成后,通过接收基阵自校准算法计算出待测基阵的接收波束指向性图、波束角度误差曲线,对双程往返测试数据求平均,降低测量误差。
将多波束测深声纳接收机电路与换能器基阵作为整体测量,不单独测量各分机特性,测量过程采取自动单向无停顿旋转方法,通过信号同步线控制数据采集,排除机械旋转空程误差,采用近场聚焦波束形成算法进行基阵校准,在较小尺寸的消声水池计算基阵波束指向性曲线和波束角度误差曲线。
步骤(1)所述发射换能器与待测基阵之间的间距,按下式计算:
Figure BDA0001786082780000021
其中D为待测基阵孔径,λ为信号波长,需要根据当前工作频率和声速计算得到。
步骤(2)所述将发射换能器调节布放深度使其对准待测基阵等效声中心的方法为:
在待测基阵等效声中心处,一般为线阵几何中心,放置一个水听器,发射换能器不断发射脉冲信号,观察水听器信号输出,同时调节发射换能器深度,选取接收信号最强位置作为发射换能器布放位置。
步骤(4)所述信号源发射信号参数与自动旋转参数控制方法为:
发射信号采用待测基阵中心频率,发射脉宽选择0.1ms长度正弦波,信号发射间隔选取1s,多自由度伸缩杆采取自动匀速旋转方式,旋转速度0.2°/s,旋转范围[-90°90°],由数据控制存储单元按照秒脉冲产生同步信号,发射换能器和待测基阵按照秒脉冲信号同步工作、记录数据,不需要严格与旋转角度对应,旋转误差通过后续算法自校正,总共采集900ping数据。
步骤(6)所述接收基阵自校准算法包含以下步骤:
(6.1)对所有数据进行近场聚焦波束形成处理,根据波束角度和布放距离预估回波时间范围,进行数据压缩、整理,减小计算数据量;
(6.2)对压缩后的波束形成输出数据选取波束形成输出极值并记录,在所有波束输出极值中寻找最大值位置D=argmax{D(n)}作为正对向发射换能器的中心波束位置;
(6.3)对中心波束位置所对应的数据文件再进行分辨率更高的近场聚焦波束形成,找出波束输出最大值位置对应的波束到达角度,观察与0°的偏差,综合测量出的基阵安装误差角度,修正旋转起始位置和终止位置角度范围,给出各旋转位置理论对应波束角;
(6.4)利用步骤(6.1)中得到的各位置回波数据的波束形成结果最大值,找出该数据对应的波束形成输出能量最强点所对应的回波角度,与步骤(6.3)中计算的理论值做差得到波束角度误差曲线;
(6.5)利用步骤(6.2)中计算得到的各位置波束形成输出极值,经角度误差曲线修正后,选取回波极值-3dB范围作为波束覆盖宽度,得到接收波束指向性图。
本发明的有益效果在于:
将多波束测深声纳接收机电路与换能器基阵作为整体测量基阵指向性,经***充分预热后排除电路温度漂移影响,使***处于日常工作状态下,不需单独测量各分机特性。测量过程采取自动单向无停顿旋转方法,通过信号同步线控制数据采集,采样时刻精准对标,不需要严格找寻旋转中心零点角度,排除机械旋转空程误差。采用近场聚焦波束形成算法进行基阵校准,可以在较小尺寸的消声水池计算基阵波束指向性曲线和波束角度误差曲线。此发明能够有效地表征宽覆盖多波束接收基阵对各角度回波的响应能力,将接收换能器基阵与信号调理电路一体化测量,更真实的反映接收***整体的信号响应能力,可以广泛地应用在多波束接收基阵校准领域。
附图说明
图1为***布放示意图。
图2为基阵旋转及波束指向性图。
图3为数据压缩后波束形成输出。
图4为基准位置波束输出。
图5为波束指向性输出图。
图6为波束误差曲线图。
具体实施方式
下面结合附图1对本发明的具体实施措施进行描述:
本发明公开了一种宽覆盖多波束接收基阵校准方法,将多波束测深声纳接收机电路与换能器基阵作为整体测量基阵指向性,经***充分预热后排除电路温度漂移影响,使***处于日常工作状态下,不需单独测量各分机特性。测量过程采取自动单向无停顿旋转方法,通过信号同步线控制数据采集,采样时刻精准对标,不需要严格找寻旋转中心零点角度,排除机械旋转空程误差。采用近场聚焦波束形成算法进行基阵校准,可以在较小尺寸的消声水池计算基阵波束指向性曲线和波束角度误差曲线。此发明能够有效地表征宽覆盖多波束接收基阵对各角度回波的响应能力,将接收换能器基阵与信号调理电路一体化测量,更真实的反映接收***整体的信号响应能力,可以广泛地应用在多波束接收基阵校准领域。
将多波束测深声纳接收机电路与换能器基阵作为整体测量,不单独测量各分机特性。测量过程采取自动单向无停顿旋转方法,通过信号同步线控制数据采集,排除机械旋转空程误差。采用近场聚焦波束形成算法进行基阵校准,在较小尺寸的消声水池计算基阵波束指向性曲线,包含以下几个步骤:
步骤(1):根据接收换能器基阵工作频率、换能器尺寸选取相应尺寸的消声水池,通过声速仪记录下测量时刻水池的实时声速值,计算基阵布放位置和布放方位向。
步骤(2):将待测换能器基阵按照布放方位向安装在多自由度伸缩杆上,测量待测换能器基阵与水平方向夹角,基阵入水放置在计算好的布放位置处。将一个全指向性发射换能器放置在计算好的发射位置处,调节布放深度使其对准待测基阵等效声中心。
步骤(3):将发射换能器与数据控制存储单元通过同步信号线相连,打开信号源设置工作参数,功率放大器开机预热。打开数据控制存储单元,多波束测深***空采预热,使用示波器实时监测发射换能器发射波形与待测基阵接收波形。
步骤(4):预热完成后将待测基阵旋转至垂直于发射换能器角度,设置多自由度伸缩杆旋转方向与旋转速度,开启多波束测深***采集软件,多自由度伸缩杆顺时针自动旋转的同时采集软件自动将对应位置的回波信号存储至数据控制存储单元。由数据控制存储单元按照秒脉冲产生同步信号,发射换能器和待测基阵采集***按照秒脉冲信号同步工作、记录数据。
步骤(5):单方向测量完成后设置参数,再进行一次逆时针自动旋转测量,采集软件自动将对应位置的回波信号存储至数据控制存储单元。
步骤(6):数据采集完成后,通过接收基阵自校准算法计算出待测基阵的接收波束指向性图、波束角度误差曲线,对双程往返测试数据求平均,降低测量误差。
本发明还可以包括:
1、发射换能器与待测基阵之间的间距:
Figure BDA0001786082780000051
其中D为待测基阵孔径,λ为信号波长,需要根据当前工作频率和声速计算得到
2、信号源发射信号参数与自动旋转参数控制方法包含但不限于:
发射信号采用待测基阵中心频率,发射脉宽选择0.1ms长度正弦波,发射间隔1s。多自由度伸缩杆采取自动匀速旋转方式,旋转速度0.2°/s,旋转范围[-90°90°]。发射信号与采集单元通过信号线进行同步,不需要严格与旋转角度对应,旋转误差通过后续算法自校正,总共采集900ping数据。
3、接收基阵自校准算法包含以下几个步骤:
步骤(1):对所有数据进行近场聚焦波束形成处理,根据波束角度和布放距离预估回波时间范围,进行数据压缩、整理,减小计算数据量。
步骤(2):对压缩后的波束形成输出数据选取波束形成输出极值并记录,在所有波束输出极值中寻找最大值位置D=argmax{D(n)}作为正对向发射换能器的中心波束位置。
步骤(3):对中心波束位置所对应的数据文件再进行分辨率更高的近场聚焦波束形成,找出波束输出最大值位置对应的波束到达角度,观察与0°的偏差,综合测量出的基阵安装误差角度,修正旋转起始位置和终止位置角度范围,给出各旋转位置理论对应波束角。
步骤(4):利用步骤(1)中得到的各位置数据的波束形成结果最大值,找出该数据对应的波束形成输出最强点对应的回波角度,与步骤(3)中计算的理论值做差即可得到波束角度误差曲线。
步骤(5):利用步骤(2)中计算得到的各位置波束形成输出极值,经角度误差曲线修正后,选取回波极值-3dB范围作为波束覆盖宽度,得到接收波束指向性图。对双程往返测试数据求平均,降低测量误差。
本发明公开了一种宽覆盖多波束接收基阵校准方法,将多波束测深声纳接收机电路与换能器基阵作为整体测量,经***充分预热后排除电路温度漂移影响,使***处于日常工作状态下,不单独测量各分机特性。测量过程采取自动单向无停顿旋转方法,通过信号同步线控制数据采集,不需要严格找寻中心零点角度,排除机械旋转空程误差。采用近场聚焦波束形成算法进行基阵校准,在较小尺寸的消声水池计算基阵波束指向性曲线和波束角度误差曲线。
结合附图1说明仪器布放结构,利用声速仪记录下测量时刻水池的实时声速值,根据待测换能器基阵尺寸和工作频率、实时声速计算基阵的远场距离
Figure BDA0001786082780000061
由于本发明使用近场聚焦波束形成技术,可以不处于远场平面波假设下,但是为了避免混响造成的测量误差,选取
Figure BDA0001786082780000062
作为布放间距。
将待测换能器基阵安装在多自由度伸缩杆上,尽量保持安装水平,同时正对对向伸缩杆。使用水平仪测量待测换能器与水平方向夹角,将换能器入水至3m处。在对向伸缩杆上放置全指向性发射换能器,连接好同步线、信号发生器与功率放大器。持续发射短脉冲,在待测换能器端使用水听器检测接收到的直达波信号,不断调节发射换能器深度,直至接收回波值最大,此处即为等效声中心位置。
将信号源与功率放大器调节至接收换能器日常工作功率与参数,同时待测换能器基阵开始空采数据,所有电子部分开始预热。这样做的目的是使电路部分处于额定工作状态,排除电路温度漂移干扰。
结合附图2,完成预热后的***旋转至垂直于发射换能器的角度,发射信号采用待测基阵中心频率,发射脉宽选择0.1ms长度正弦波,发射间隔1s。多自由度伸缩杆采取自动匀速旋转方式,旋转速度0.2°/s,旋转范围[-90°90°]。由数据控制存储单元按照秒脉冲产生同步信号,发射换能器和待测基阵按照秒脉冲信号同步工作、记录数据。旋转过程中信号采集位置不需要严格与旋转角度对应,初始位置误差通过后续算法自校正,总共采集900ping数据。
完成顺时针方向测量后,切换数据文件,在顺时针测量终止位置处反向自动旋转换能器,再进行一次测量并存储数据。
数据采集完成后,通过接收基阵自校准算法计算出待测基阵的接收波束指向性图、波束角度误差曲线。
接收基阵自校准算法包含以下几个步骤:
步骤(1):对所有数据进行近场聚焦波束形成处理,根据波束角度和布放距离预估回波时间范围,进行数据压缩、整理,减小计算数据量。
对于接收到的数据文件,按照探测ping号切割成900个小数据文件,近场聚焦波束形成算法对数据进行处理,利用公式:
Figure BDA0001786082780000063
即可计算出近场聚焦波束形成输出,其中波长λ需要按照记录下的测量现场声速实时计算。按照波束角度计算回波斜距
Figure BDA0001786082780000071
Figure BDA0001786082780000072
范围内的数据提取出来作为应用的数据。
步骤(2):对压缩后的波束形成输出数据选取波束形成输出极值并记录,在所有波束输出极值中寻找最大值位置D=argmax{D(n)}作为正对向发射换能器的中心波束位置。按照同相叠加原理,只有在换能器基阵正前方位置的波束的同相叠加最为充分,同时聚焦传输距离最近,波束输出形成最强,记录此位置波束号,作为中心波束位置。
步骤(3):对中心波束位置所对应的数据文件再进行N=1024点的近场聚焦波束形成,在角度轴上更细的划分出波束到达角度,找出波束输出最大值位置对应的波束到达角度,观察与0°的偏差,综合测量出的基阵安装误差角度,修正旋转起始位置和终止位置角度范围,新计算旋转的起始位置与终止位置。由此,可以根据伸缩杆转速与ping号计算出每一组数据的真实空间位置,给出各旋转位置理论对应波束角。
步骤(4):利用步骤(1)中得到的各位置数据的波束形成结果最大值,找出该数据对应的波束形成输出最强点对应的回波角度,与步骤(3)中计算的理论值做差即可得到波束角度误差曲线。将此曲线记录下来,既可作为多波束测深***波束到达角度的修正曲线。这条曲线反映的是换能器与接收电路共同作用下的多波束测深声纳接收基阵性能,是综合所有影响因素下的校正曲线,可以为测深算法性能评估提供参考。
步骤(5):利用步骤(2)中计算得到的各位置波束形成输出极值,经角度误差曲线修正后,选取回波极值-3dB范围作为波束覆盖宽度,得到接收波束指向性图。对双程往返测试数据求平均,降低测量误差。
多波束测深声纳实际工作过程中首先就需要将数据变换到波束域,因此直接的基阵自然指向性图并不能真实的反映接收基阵的性能,经过本发明的算法处理后,能够给出真实的***波束角度响应性能,能够有效地评估基阵对外侧波束角度情况下的信号接收能力。
下面结合本方法的一种较佳实现范例进行更详细的说明:
针对某型多波束测深***,工作频率f0=200kHz,阵元数目N=100,经计算得到基阵远场距离R=37.5m,因此选取布放距离L=10m。按照作用斜距计算回波到达时间,对数据进行近场聚焦波束形成计算,数据压缩后的输出如图3所示,其中的亮点部分即为每个时间片对应的波束输出极值位置。
找寻到各极值中的最大值位置后,对基准位置的数据进行更为细致的近场聚焦波束形成,波束输出如图4所示,选取波束输出最大值位置的角度与中心原点比较,即可得到安装偏差值,之后对测量起始位置和终止位置角度进行修正。
经角度修正后的极值,选取回波-3dB范围作为波束覆盖宽度,得到接收波束指向性图如图5所示,波束角度误差曲线如图6。
经过处理可以发现,对于基阵外侧波束的偏差能够达到0.4°,这将对多波束测深声纳的探测效果产生极大影响,经计算在水深50m情况下,基阵波束角度偏差引入的计算精度偏差达到近40cm,所以通过实验数据处理可以发现,进行多波束测深声纳的接收基阵校准是非常有必要的。
一种宽覆盖多波束接收基阵校准方法,将多波束测深声纳接收机电路与换能器基阵作为整体测量,不单独测量各分机特性。测量过程采取自动单向无停顿旋转方法,通过信号同步线控制数据采集,排除机械旋转空程误差。采用近场聚焦波束形成算法进行基阵校准,在较小尺寸的消声水池计算基阵波束指向性曲线和波束角度误差曲线。
包含以下几个步骤:
步骤(1):根据接收换能器基阵工作频率、换能器尺寸选取相应尺寸的消声水池,通过声速仪记录下测量时刻水池的实时声速值,计算基阵布放位置和布放方位向。
步骤(2):将待测换能器基阵按照布放方位向安装在多自由度伸缩杆上,测量待测换能器基阵与水平方向夹角,基阵入水放置在计算好的布放位置处。将一个全指向性发射换能器放置在计算好的发射位置处,调节布放深度使其对准待测基阵等效声中心。
步骤(3):将发射换能器与数据控制存储单元通过同步信号线相连,打开信号源设置工作参数,功率放大器开机预热。打开数据控制存储单元,多波束测深***空采预热,使用示波器实时监测发射换能器发射波形与待测基阵接收波形。
步骤(4):预热完成后将待测基阵旋转至垂直于发射换能器角度,设置多自由度伸缩杆旋转方向与旋转速度,开启多波束测深***采集软件,多自由度伸缩杆顺时针自动旋转的同时采集软件自动将对应位置的回波信号存储至数据控制存储单元。由数据控制存储单元按照秒脉冲产生同步信号,发射换能器和待测基阵采集***按照秒脉冲信号同步工作、记录数据。
步骤(5):单方向测量完成后设置参数,再进行一次逆时针自动旋转测量,采集软件自动将对应位置的回波信号存储至数据控制存储单元。
步骤(6):数据采集完成后,通过接收基阵自校准算法计算出待测基阵的接收波束指向性图、波束角度误差曲线,对双程往返测试数据求平均,降低测量误差。
发射换能器与待测基阵之间的间距:
Figure BDA0001786082780000081
其中D为待测基阵孔径,λ为信号波长,需要根据当前工作频率和声速计算得到。
将发射换能器调节布放深度使其对准待测基阵等效声中心的方法为:
在待测基阵等效声中心处(一般为线阵几何中心)放置一个水听器,发射换能器不断发射脉冲信号,观察水听器信号输出。同时调节发射换能器深度,选取接收信号最强位置作为发射换能器布放位置。
信号源发射信号参数与自动旋转参数控制方法包含但不限于:
发射信号采用待测基阵中心频率,发射脉宽选择0.1ms长度正弦波。为了排除信号混响干扰,提高测量精度,信号发射间隔选取1s。多自由度伸缩杆采取自动匀速旋转方式,旋转速度0.2°/s,旋转范围[-90°90°]。由数据控制存储单元按照秒脉冲产生同步信号,发射换能器和待测基阵按照秒脉冲信号同步工作、记录数据,不需要严格与旋转角度对应,旋转误差通过后续算法自校正,总共采集900ping数据。
接收基阵自校准算法包含以下几个步骤:
步骤(1):对所有数据进行近场聚焦波束形成处理,根据波束角度和布放距离预估回波时间范围,进行数据压缩、整理,减小计算数据量。
步骤(2):对压缩后的波束形成输出数据选取波束形成输出极值并记录,在所有波束输出极值中寻找最大值位置D=argmax{D(n)}作为正对向发射换能器的中心波束位置。
步骤(3):对中心波束位置所对应的数据文件再进行分辨率更高的近场聚焦波束形成,找出波束输出最大值位置对应的波束到达角度,观察与0°的偏差,综合测量出的基阵安装误差角度,修正旋转起始位置和终止位置角度范围,给出各旋转位置理论对应波束角。
步骤(4):利用步骤(1)中得到的各位置回波数据的波束形成结果最大值,找出该数据对应的波束形成输出能量最强点所对应的回波角度,与步骤(3)中计算的理论值做差即可得到波束角度误差曲线。
步骤(5):利用步骤(2)中计算得到的各位置波束形成输出极值,经角度误差曲线修正后,选取回波极值-3dB范围作为波束覆盖宽度,得到接收波束指向性图。

Claims (6)

1.一种宽覆盖多波束接收基阵校准方法,其特征在于,该方法包含以下步骤:
(1)根据接收换能器基阵工作频率、换能器尺寸选取相应尺寸的消声水池,通过声速仪记录下测量时刻水池的实时声速值,计算基阵布放位置和布放方位向;
(2)根据选取的合适尺寸的消声水池以及计算得到的基阵布放位置和布放方位向,将待测换能器基阵按照布放方位向安装在多自由度伸缩杆上,测量待测换能器基阵与水平方向夹角,基阵入水放置在计算好的布放位置处,将一个全指向性发射换能器放置在计算好的发射位置处,调节布放深度使其对准待测基阵等效声中心;
(3)根据调节好的布放深度,将发射换能器与数据控制存储单元通过同步信号线相连,打开信号源设置工作参数,功率放大器开机预热,打开数据控制存储单元,多波束测深***空采预热,使用示波器实时监测发射换能器发射波形与待测基阵接收波形;
(4)预热完成后将待测基阵旋转至垂直于发射换能器角度,设置多自由度伸缩杆旋转方向与旋转速度,开启多波束测深***采集软件,多自由度伸缩杆顺时针自动旋转的同时采集软件自动将对应位置的回波信号存储至数据控制存储单元,由数据控制存储单元按照秒脉冲产生同步信号,发射换能器和待测基阵采集***按照秒脉冲信号同步工作、记录数据;
(5)单方向测量完成后设置参数,再进行一次逆时针自动旋转测量,采集软件自动将对应位置的回波信号存储至数据控制存储单元;
(6)数据采集完成后,通过接收基阵自校准算法计算出待测基阵的接收波束指向性图、波束角度误差曲线,对双程往返测试数据求平均,降低测量误差。
2.根据权利要求1所述的一种宽覆盖多波束接收基阵校准方法,其特征在于:将多波束测深声纳接收机电路与换能器基阵作为整体测量,不单独测量各分机特性,测量过程采取自动单向无停顿旋转方法,通过信号同步线控制数据采集,排除机械旋转空程误差,采用近场聚焦波束形成算法进行基阵校准,在较小尺寸的消声水池计算基阵波束指向性曲线和波束角度误差曲线。
3.根据权利要求1所述的一种宽覆盖多波束接收基阵校准方法,其特征在于:
步骤(1)所述发射换能器与待测基阵之间的间距,按下式计算:
Figure FDA0001786082770000011
其中D为待测基阵孔径,λ为信号波长,需要根据当前工作频率和声速计算得到。
4.根据权利要求3所述的一种宽覆盖多波束接收基阵校准方法,其特征在于:步骤(2)所述将发射换能器调节布放深度使其对准待测基阵等效声中心的方法为:
在待测基阵等效声中心处,一般为线阵几何中心,放置一个水听器,发射换能器不断发射脉冲信号,观察水听器信号输出,同时调节发射换能器深度,选取接收信号最强位置作为发射换能器布放位置。
5.根据权利要求4所述的一种宽覆盖多波束接收基阵校准方法,其特征在于:步骤(4)所述信号源发射信号参数与自动旋转参数控制方法为:
发射信号采用待测基阵中心频率,发射脉宽选择0.1ms长度正弦波,信号发射间隔选取1s,多自由度伸缩杆采取自动匀速旋转方式,旋转速度0.2°/s,旋转范围[-90°90°],由数据控制存储单元按照秒脉冲产生同步信号,发射换能器和待测基阵按照秒脉冲信号同步工作、记录数据,不需要严格与旋转角度对应,旋转误差通过后续算法自校正,总共采集900ping数据。
6.根据权利要求5所述的一种宽覆盖多波束接收基阵校准方法,其特征在于:步骤(6)所述接收基阵自校准算法包含以下步骤:
(6.1)对所有数据进行近场聚焦波束形成处理,根据波束角度和布放距离预估回波时间范围,进行数据压缩、整理,减小计算数据量;
(6.2)对压缩后的波束形成输出数据选取波束形成输出极值并记录,在所有波束输出极值中寻找最大值位置作为正对向发射换能器的中心波束位置,最大值位置按下式计算:
D=argmax{D(n)}
上式中,D(n)表示第n波束的输出极值;
(6.3)对中心波束位置所对应的数据文件再进行分辨率更高的近场聚焦波束形成,找出波束输出最大值位置对应的波束到达角度,观察与0°的偏差,综合测量出的基阵安装误差角度,修正旋转起始位置和终止位置角度范围,给出各旋转位置理论对应波束角;
(6.4)利用步骤(6.1)中得到的各位置回波数据的波束形成结果最大值,找出该数据对应的波束形成输出能量最强点所对应的回波角度,与步骤(6.3)中计算的理论值做差得到波束角度误差曲线;
(6.5)利用步骤(6.2)中计算得到的各位置波束形成输出极值,经角度误差曲线修正后,选取回波极值-3dB范围作为波束覆盖宽度,得到接收波束指向性图。
CN201811017290.5A 2018-09-01 2018-09-01 一种宽覆盖多波束接收基阵校准方法 Active CN109283511B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811017290.5A CN109283511B (zh) 2018-09-01 2018-09-01 一种宽覆盖多波束接收基阵校准方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811017290.5A CN109283511B (zh) 2018-09-01 2018-09-01 一种宽覆盖多波束接收基阵校准方法

Publications (2)

Publication Number Publication Date
CN109283511A CN109283511A (zh) 2019-01-29
CN109283511B true CN109283511B (zh) 2022-07-29

Family

ID=65184294

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811017290.5A Active CN109283511B (zh) 2018-09-01 2018-09-01 一种宽覆盖多波束接收基阵校准方法

Country Status (1)

Country Link
CN (1) CN109283511B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109839627B (zh) * 2019-03-08 2022-12-13 哈尔滨工程大学 一种水池条件下声呐探测距离量值计量校准***及方法
CN110361720B (zh) * 2019-07-31 2021-04-02 武昌理工学院 一种声呐收声板控制方法
CN110865333B (zh) * 2019-11-19 2021-08-24 浙江大学 洋流影响下水下滑翔机单信标无源声学定位方法
CN111208522A (zh) * 2020-01-15 2020-05-29 上海船舶电子设备研究所(中国船舶重工集团公司第七二六研究所) 岸基高频多波束图像声纳***
CN113495096B (zh) * 2020-06-18 2024-04-09 中控全世科技(杭州)有限公司 一种基于fid-tvoc检测仪的校准方法
CN112327280B (zh) * 2020-10-27 2022-05-17 哈尔滨工程大学 一种水声均匀线阵阵元编号校正方法
CN112505667B (zh) * 2020-11-19 2022-07-15 哈尔滨工程大学 一种二维声呐基阵运动姿态自校准方法
CN113153232B (zh) * 2021-03-11 2024-04-02 浙江大学 一种基于小型阵列的海底天然气水合物气泡泄露定位装置及方法
CN114113314B (zh) * 2021-11-25 2024-03-22 厦门大学 一种用于滨海核电站的海洋致灾生物水声监测***及方法
CN117630894B (zh) * 2024-01-24 2024-04-12 山东省科学院海洋仪器仪表研究所 一种多波束声呐平均阵元间距校准方法
CN118074771B (zh) * 2024-04-18 2024-06-21 智慧尘埃(上海)通信科技有限公司 一种通信感知一体化方法、装置、设备及存储介质

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1489442A (en) * 1974-02-21 1977-10-19 Westinghouse Electric Corp Sonar apparatus
US5530678A (en) * 1994-12-05 1996-06-25 Alliant Techsystems Inc. Real-time calibration acoustic array
WO2001082806A1 (en) * 2000-04-21 2001-11-08 Txsonics Ltd. Systems and methods for testing and calibrating a focused ultrasound transducer array
CN1818710A (zh) * 2006-03-15 2006-08-16 哈尔滨工程大学 8基元超短基线定位***相位测量***及其校准方法
CN101644913A (zh) * 2009-08-20 2010-02-10 中国科学院声学研究所 一种水下授时和同步方法及其***
CN103777177A (zh) * 2013-11-08 2014-05-07 中国船舶重工集团公司七五○试验场 一种基于宽带信号时延检测的超短基线水下目标定位方法
CN104181540A (zh) * 2014-09-05 2014-12-03 哈尔滨工程大学 一种基于编码信号的同时多焦点浅水多波束接收动态聚焦***
CN104267403A (zh) * 2014-09-05 2015-01-07 哈尔滨工程大学 一种浅水多波束测深仪用快速动态聚焦方法
CN104787261A (zh) * 2015-04-14 2015-07-22 重庆交通大学 挂舷式内河离港船舶吃水深度测量装置及其控制方法
CN105891835A (zh) * 2016-06-16 2016-08-24 北京海卓同创科技有限公司 一种实时动态聚焦波束形成方法及***
CN107037418A (zh) * 2016-11-02 2017-08-11 哈尔滨工程大学 一种低频超指向性多极子矢量阵的有界空间校正方法
WO2017190006A1 (en) * 2016-04-29 2017-11-02 R2Sonic, Llc Sonar data compression
CN107703501A (zh) * 2017-09-28 2018-02-16 哈尔滨工程大学 一种校准多波束声呐测深和分辨率的装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8532951B2 (en) * 2009-12-22 2013-09-10 Delphinus Medical Technologies, Inc. Method for calibrating a transducer array

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1489442A (en) * 1974-02-21 1977-10-19 Westinghouse Electric Corp Sonar apparatus
US5530678A (en) * 1994-12-05 1996-06-25 Alliant Techsystems Inc. Real-time calibration acoustic array
WO2001082806A1 (en) * 2000-04-21 2001-11-08 Txsonics Ltd. Systems and methods for testing and calibrating a focused ultrasound transducer array
CN1818710A (zh) * 2006-03-15 2006-08-16 哈尔滨工程大学 8基元超短基线定位***相位测量***及其校准方法
CN101644913A (zh) * 2009-08-20 2010-02-10 中国科学院声学研究所 一种水下授时和同步方法及其***
CN103777177A (zh) * 2013-11-08 2014-05-07 中国船舶重工集团公司七五○试验场 一种基于宽带信号时延检测的超短基线水下目标定位方法
CN104181540A (zh) * 2014-09-05 2014-12-03 哈尔滨工程大学 一种基于编码信号的同时多焦点浅水多波束接收动态聚焦***
CN104267403A (zh) * 2014-09-05 2015-01-07 哈尔滨工程大学 一种浅水多波束测深仪用快速动态聚焦方法
CN104787261A (zh) * 2015-04-14 2015-07-22 重庆交通大学 挂舷式内河离港船舶吃水深度测量装置及其控制方法
WO2017190006A1 (en) * 2016-04-29 2017-11-02 R2Sonic, Llc Sonar data compression
CN105891835A (zh) * 2016-06-16 2016-08-24 北京海卓同创科技有限公司 一种实时动态聚焦波束形成方法及***
CN107037418A (zh) * 2016-11-02 2017-08-11 哈尔滨工程大学 一种低频超指向性多极子矢量阵的有界空间校正方法
CN107703501A (zh) * 2017-09-28 2018-02-16 哈尔滨工程大学 一种校准多波束声呐测深和分辨率的装置

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
Evaluation of target angular position algorithms for multi-beam fishery sonars;Vatnehol, SN and Ona, E;《JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA》;20170223;第141卷(第2期);第1074-1083页 *
Minghui Li ; Michelle McGuire ; Kwok Shun Ho ; Gordon Hayward.Array element failure correction for robust ultrasound beamforming and imaging.《2010 IEEE International Ultrasonics Symposium》.2011,第29-32页. *
Ying Jiang ; Wen Xu ; Xiang Pan.Near-field beamforming for a Multi-Beam Echo Sounder: Approximation and error analysis.《OCEANS》.2010,第1-4页. *
多波束测深声呐多通道接收机的设计与实现;尹家亮;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20180315;第C036-300页 *
多波束测深声纳技术研究新进展;李海森,周天,徐超;《声学技术》;20130415;第32卷(第2期);第73-80页 *
多波束测深***换能器的安装校准分析;张海涛,唐秋华,周兴华,丁继胜,刘建立;《海洋通报》;20090215;第28卷(第1期);第102-107页 *
矢量阵近场聚焦波束形成技术研究;李萌;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20100615(第6期);第C028-38页 *

Also Published As

Publication number Publication date
CN109283511A (zh) 2019-01-29

Similar Documents

Publication Publication Date Title
CN109283511B (zh) 一种宽覆盖多波束接收基阵校准方法
CN100520441C (zh) 8基元超短基线定位***相位测量***及其校准方法
CN102066921A (zh) 利用超声波来无损检测受检件的方法以及执行该方法的装置
CN109991590B (zh) 一种在有限空间压力罐内测试换能器低频发射特性的***与方法
Lanzoni et al. High-resolution calibration of a multibeam echo sounder
CN108226884B (zh) 一种有源相控阵雷达天线近场校准的方法
Lanzoni et al. A method for field calibration of a multibeam echo sounder
CN102087357B (zh) 传感器阵列回波方向估计及多波束回波测深底检测方法
CN110275139B (zh) 一种基于旋转式基元复用的超短基线定位***及方法
CN113534161B (zh) 一种用于远距离定位水下声源的波束镜像聚焦方法
CN109752705B (zh) 高频水声阵列性能参数测量方法及***、设备及存储介质
CN113670433B (zh) 一种能提高水听器阵列阵元在线灵敏度测量精度的方法
CN109375197B (zh) 一种小尺寸矢量阵低频散射校正方法
CN111735436A (zh) 一种基于3条以上均匀分布多波束数据的海底地形数据校验方法
Lauterbach et al. Procedure for the accurate phase calibration of a microphone array
Li et al. Calibration of multibeam echo sounder transducer array based on focused beamforming
Lanzoni et al. Calibration of multibeam echo sounders: a comparison between two methodologies
CN114624660A (zh) 天线发射方向图、接收方向图以及波束方向图测试方法
JP2011247624A (ja) 水中探知装置及び水中探知方法
CN112327280B (zh) 一种水声均匀线阵阵元编号校正方法
CN117630894B (zh) 一种多波束声呐平均阵元间距校准方法
CN112526464B (zh) 一种基于多通道雷达实测数据估计方位波束宽度的方法
RU2815608C1 (ru) Гидроакустическая навигационная система с ультракороткой измерительной базой
CN113687303B (zh) 甲板定位***与主***
CN117872274A (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