多波束勘测系统工作基础学习知识原理及其结构

1、第二章多波束勘测系统工作原理及结构 多波束系统是70年代兴起、80年代中、末期又得到飞速发展的一项全新的海底地形 精密勘测技术。它是当前兴趣的焦点，因为它既有条带测深数据，又同时可获取反映底 质属性的回波强度数据(Laurent Hellequin et al.,2003)。该技术采取广角度定向发射和多通 道信息接收，获得水下高密度具有上百个波束的条幅式海底地形数据，彻底改变了传统 测深技术概念，使测深原理、勘测方法、外围设备和数据处理技术诸方面都发生了巨大 变化，大大提高了海底地形勘测的精度、分辨率和工作效率，实现了测深技术史上的一 次革命性突破(李家彪等，2000)。多波束系统的工作原理与

2、传统的单波束回声测深仪工 作原理类似，都是根据声波在水下往返传播的时间与声速的乘积得到距离，从而得到水 深。不同的是单波束测深仪一般采用较宽的发射波束(8 左右)向船底垂直发射，声传 播路径不会发生弯曲，来回的路径最短，能量衰减很小，通过对回声信号的幅度检测确 定信号往返传播的时间，再根据声波在水介质中的平均传播速度计算测量水深。在多波 束系统中，换能器配置有一个或者多个换能器单元的阵列，通过控制不同单元的相位， 形成多个具有不同指向角的波束，通常只发射一个波束而在接收时形成多个波束。除换 能器天底波束外，外缘波束随着入射角的增加，波束在倾斜穿过水层时会发生折射，同 时由于多波束沿航迹方向采用

3、较窄的波束角而在垂直航迹方向采用较宽的覆盖角，要获 得整个测幅上精确的水深和位置，必须要精确地知道测量区域水柱的声速剖面和波束在 发射和接收时船的姿态和船艏向。因此，多波束测深在系统组成和测量时比单波束测深 仪要复杂得多(周兴华等，1999 )。 2.1多波束勘测系统的工作原理 2.1.1单波束的形成 2.1.1.1发射阵和波束的形成 一个单波束在水中发射后，是球形等幅度传播，所以方向上的声能相等。这种均匀 传播称为各向同性传播(isotropic expansion)，发射阵也叫各向同性源(isotropic source)。 例如，一个小石头扔进池塘时就是这种情况，如图2.7所示。 XTT

4、HtirMaiif 图2.7波的各向同性传播 显然，测深时是不能采用如此的声波的。采用发射基阵就可以产生各向异性的声波。 下面简要叙述它形成的原理。 如果两个相邻的发射器发射相同的各向同性的声信号，声波图将互相重叠和干涉， 如图2.8所示。两个波峰或者两个波谷之间的叠加会增强波的能量，波峰与波谷的叠加正 好互相抵消，能量为零。一般地，相长干涉发生在距离每个发射器相等的点或者整波长 处，而相消干涉发生在相距发射器半波长或者整波长加半波长处。显然，水听器需要放 置在相长干涉处。 3C *柘上干円 图2.8 相长干涉和相消干涉(Constructive and Destructive Interfe

5、rence ) 一个典型的声纳，基阵的间距 d （图2.8中1、2点的距离）是入/2 （半波长）。在这种情况 下，相长和相消干涉发生时的点位处于最有利的角度（点位与基阵中心的连线与水平线 的夹角），相长干涉：B = 0, 180，相消干涉：B = 90, 270，如图2.9所示。 图2.9两个发射器相距入/2时的相长和相消干涉 图2.10是两个发射器间距 入12时的波束能量图(Beam Pattern),左边为平面图，右边 为三维图，从图上可以清楚地看到能量的分布，不同的角度有不同的能量，这就是能量 的指向性(directivity )。如果一个发射阵的能量分布在狭窄的角度中，就称该系统指向性

6、 高。真正的发射阵由多个发射器组成，有直线阵和圆形阵等。这里只讨论离散直线阵， 其它阵列类似可以推导出。如图2.11所示，根据两个发射器的基阵可以推导出多个发射器 组成的直线阵的波束图。 图2.10两个发射器间距 入/2时的波束能量图(Beam Pattern) 图2.11多基元线性基阵的波束图( Beam Pattern) 图2.11中，能量最大的波束叫主瓣，侧边的一些小瓣是旁瓣，也是相长干涉的地方， 引起了能量的泄漏。旁瓣还可能引起回波，对主瓣的回波产生干扰。旁瓣是不可避免的， 可以通过加权的方法降低旁瓣的水平，但是加权后旁瓣水平值降低了，波束却展宽了。 主瓣的中心轴叫最大响应轴(maxi

7、mum response axis MRA),主瓣半功率处(相对于主 瓣能量的-3db )角度的两倍就是波束角。发射器越多，基阵越长，则波束角越小，指向 性就越高。设基阵的长度为 D,则波束角 e = 50.6 xx /D(2.36) 可以看出，减小波长或者增大基阵的长度都可以提高波束的指向性。但是，基阵的长度 不可能无限增大，而波长越小，在水中衰减得越快，所以指向性不可能无限提高。 2.1.1.2 波束的指向(Beam Steering) 换能器怎样在指定的方向上发射或者接收声波，称为波束的指向。以水听器接收回 波为例。如图2.12，当回波以B方向到达接收基阵时，首先在点3到达，其次为点2和

8、点1, 相应的增加的时间为 (2.37) T2=A/c =(d sin 0 )/c T1 = B/c =(2d sin 0 )/c(2.38) 计算出偏移时间后，在基阵中作相应的调整，弓I入延时，使回波在基阵上正好构成 相长干涉，这样就可以使主瓣在指定的方向上，如图2.13所示。 图2.13引入延时后主瓣方向的偏移 图2.14多波束的几何构成 2.1.2多波束的形成 当接收波束发射出扇形波束后，接收波束按一定的间距(等距离或者等角度)与之 相交，就形成了一个个在纵横向的窄波束脚印，如图2.14所示。 设水听器共有N个基元，每个基元i记录的回波S(t)的振幅为A(t),且 写成相位的形式为 S(

9、t) = A(t)cos(2pft) S(t) = A(t)cos(f(t) (2.39) (2.40) S(t) A(t)ej (t) (2.41) 其中，f(t) = 2pft。 多波束系统需测量回波S(t)和相位0 (t)，然后将模拟接收信号转换为数字信号，采用 率一般在13ms之间。所有基元在采样点上的回波和相位值称为时间片(time slice )。 在上节中，讨论了基元 i相对于第一个基元的距离差，转换为相位差为 id sin (2.42) 由(2.41)和(2.42) 得 j( i id sin ) Bi( ) Ae 2 j id sin Aej i e j( id sin (2

10、.43) 其中Bi()为第i个基元在角B方向接收时的回波，则基阵接收的回波为 B() N 1j(id sin ) 0 (2.44) 其中S为加权系数。如果要求在一个时间片 (time slice)里，由N个基元形成 M个指定方向的 波束，用矩阵表示为 B( B( 1) 2 ) D11 D21 D12 D22 D1N D2N 51 52 (2.45) B( DM2 D MN SN 其中，Dij为接收角j的波束时的第i个基元的相位差，为 e 时间(ms级)完成计算，必须采用一些快速算法。这里，引入快速傅立叶变换 式(2.44)类似于傅立叶变换，设 j ( id sin j ) 。为了在如此短的 (

11、FFT)。 B( N hiJ晋 i 0 (2.46) sin (2.47) k同k的关 由于k必须为整数，所以k的取值受到一定的限制，如d、N在一定值时, 系表2.1所示。 表2.1 k和&的对应值 K 0 1 2 3 4 5 0 ( ) 0 2.4 4.8 7.2 9.6 12.0 2.1.3多波束脚印的归位 波束脚印的归位是多波束数据处理的关键问题之一。多波束测量的最终成果是得到 地理坐标系（或地方系）下的海底地形或者地物，由于多波束采用广角度定向发射、多 阵列信号接收和多个波束形成处理等技术，为了更好的确定波束的空间关系和波束脚印 的空间位置，必须首先定义多波束船体参考坐标系VFS,并根

12、据船体坐标系同地理坐标 系LLS之间的关系，将波束脚印的船体坐标转化到地理坐标系（或当地坐标系）和某一 高程基准面下的平面坐标和高程。该过程即为波束脚印的归位。船体坐标系原点位于换 能器中心，x轴指向航向，z轴垂直向下，y轴指向侧向，与x、z轴构成右手正交坐标系。 地理坐标系原点为换能器中心，x轴指向地北子午线，y同x垂直指向东，z与x、y轴构 成正交坐标系。 归位需要的参数包括船位、船姿、声速断面、波束到达角和往返程时间。归位过程 包括如下四个步骤： （1）姿态改正。 （2）船体坐标系下波束脚印位置的计算。 （3）波束脚印地理坐标的计算。 （4）海底点高程的计算。 为方便波束脚印在船体坐标系

13、下坐标的计算（声线跟踪），现作如下假设： （1）换能器处于一个平均深度，静、动吃水认为仅对深度有影响，而对平面坐标没有影 响。 （2）认为波束的往、返程路径重合。 （3）对于高频发射系统，换能器的航向变化影响可以忽略。 波束脚印船体坐标的计算需要用到三个参量，即垂直参考面下的波束到达角、传播 时间和声速剖面。由于海水的作用，声束在海水中不是沿直线传播，而是在不同介质层 的界面处发生折射，因此声束在海水中的传播路径为一折线。为了得到波束脚印的真实 位置，就必须沿着波束的实际传播路径跟踪波束，该过程即为声线跟踪，通过声线跟踪 得到波束脚印船体坐标的计算过程被称为声线弯曲改正。为了计算方便，对声速断

14、面作 如下假设： （1）声速断面是精确的，无代表性误差。 （2）声速在波束形成的垂面内变化，不存在侧向变化。 （3）声速在海水中的传播特性遵循 Snell法则。 换能器的动吃水引起的声速剖面的变化对深度的计算可以忽略不计。 根据上述讨论和假设，波束脚印的计算模型可表达为： Snell法则可描述为: Co sin Ci sin 1 Cn sin n (2.48) (x, y, z)可表 将波束的实际传播路径进行微分，则波束脚印在船体坐标系下的点位 达为： (2.49) z z0C(z)cos( (z)dz x xoC(z)sin( (z)dz 其一级近似式为: CoT zZ0 xX。 y o c

15、os CoTp 2 sin o (2.5。) 更精确的公式见 2.1.3。 波束脚印的船体坐标系确定后，下一步就可以转化为地理坐标。转换关系为: x xox R(h,r,p)( 2.51) y g yo gGy 式中，下脚g、gG别代表波束脚印的地理坐标、利用GPS确定的船体地理坐标，R(h,r,p) 为船体坐标系与地理坐标系的旋转关系，航向、横摇和纵摇是三个欧拉角。 式(2.49)确定的深度z仅为换能器面到达海底的垂直距离，测点的实际深度还应该考 虑换能器的静吃水 hss、动吃水hds、船体姿态对深度的影响ha，若潮位的变化htide是相对 于某一深度基准面或者高程基准面确定的，则波束脚印的

16、高程为： zg z Aide hss hds ha（ 2.52） 换能器的静吃水在测量前或换能器安装后被量定，作为一个常量输入到多波束的数据处 理单元中；动吃水是由于船体的运动而产生的，它可通过姿态传感器中的Heaven参数确 定。船体姿态对波束脚印的地理坐标也有一定的影响，可通过姿态传感器的横摇r和纵 摇p参数确定。上述参数的测定及其对波束脚印平面位置和深度的补偿属于纯几何问题， 武汉大学的赵建虎博士对此有详细的研究，本文不再赘述。 2.2多波束勘测系统的组成结构 2.2.1多波束的组成 多波束系统主要由三个部分组成。第一部分是多波束的主系统，主要包括换能器阵 列，收发器和数据处理、显示和记

17、录单元等；第二部分是辅助系统，包括定位系统、船 姿（横摇、纵摇、起伏和船艏向）测量传感器和测量水柱声速剖面的声速仪；第三部分 是数据存贮和后处理系统，包括数据处理计算机、数据存贮设备和绘图仪等（周兴华等， 1999）。 数据网 GPS 声速断 罗经 姿态传感器 监控器 操作和检 导航监控器 后处理 操作和检 I二测单元 测单元 ransceiver 数据存储 实时数据处理工作站 声纳影像记录 换能器 数据存储 绘图仪 打印机 图2.1: Simrad EM950/1000 多波束声纳系统组成单元 2.2.2换能器的物理构成 换能器是用来作为电声能量转换的重要器件。通常把电能转换成声能的器件称为

18、发 射换能器，把水下声能转换成电能的器件称为接收换能器（或水听器），许多主动声纳中 采用同一的换能器兼作发射和接收（秦臻，1984）。 水声换能器技术由 3类组成，一是水声换能器材料；二是水声换能器设计；三是水 声换能器制作。水声换能器材料方面涉及材料的配置、生成及成型等3方面的技术。水 声换能器设计方面涉及性能设计及结构设计等两个方面的技术。水声换能器制作涉及加 工、安装及检测等 3方面的技术。这3类水声换能器技术，是获取高性能、高质量水声 换能器的保障。这 3类水声换能器技术的发展，直接促进了水声换能器的发展，影响着 声纳的发展。 水深换能器就其所采用的换能器材料大致可分为两大类，一类是磁

19、致伸缩换能器， 一类是电致伸缩换能器。磁致伸缩材料包含着具有水久磁矩的原子，它们以这种方式被 耦合到晶格上，使之在加磁场时能够改变晶格的间距，从而相应地改变宏观样品的尺寸。 由于磁致伸缩力与晶格变化有关，所以这个力可以和固体材料受机械变形所产生的弹性 力相比较。磁致伸缩换能器把大的恒定磁场和较小的交变磁场叠加起来使之线性化，这 种线性化的工作方式常常称为压磁，它最适合在声阻抗比较高的介质 （如水）中工作（阎福 旺等）。不过，这类材料在高频工作时有着高的涡流损耗和磁致损耗，以致电声转换效率 降低，而且工作时，还需外加极化偏置。因此，目前多数声纳换能器采用了具有良好机 电性能的电致伸缩材料(秦臻，

20、1984)。电致伸缩材料是磁致伸缩材料的电学类比，它有 永久性的电矩。电矩与晶格之间这样耦合，使之在外加电场时，尺寸发生变化。实际上 这些极化了的材料，由于已经线性化，因而能够用描述压电材料的方法描述它。石英或 许是最早实际用于换能器的材料，此后，被广泛用于实验工作、声能学和超声学方面。 大约在1950年前后，能够实际应用的电致伸缩材料开始出现，这些材料能够做成陶瓷， 在极化后有极好的压电性质。钛酸钡是这类材料中首先被广泛使用的材料。锆钛酸铅现 在很多应用中取代了钛酸钡(阎福旺等) 。 换能器是一个将电能转换成声能及逆过程的装置。不管使用的多波束是哪种类型， 换能器是单个的或是一个阵列，都必须

21、要在测量的海底形成一个声照射区(en so nified area)，换能器通过发射声脉冲并接收产生的回波获得水深测量值。换能器的大小是根据 波束的宽度需要来设计的。波束的宽度被定义为从天底到扩展波面半功率点的角距离， 半功率的量值是-3dB。如图2-5所示Wells ,1996。 波束宽度 图2-5波束宽度 每个多波束系统都有一个工作频率，换能器是由多个阵元组成的阵列。整个换能器 阵列的大小是根据这样一个总的原则确定的，即波束宽度与横跨孔径的波长成反比。反 之，如果已选择了一个需要的波长，孔径的大小就可以用这个原则来确定de Moustier, 1996。 如果： 需要的波束宽度 3 =2(

22、 =0.035rad ) 选择的频率为F=100KHz , 声速C=1500m/s 波长 入=匸=0.015m 由总的原则：孔径= 得到： 孔径=0.43m 很明显波束宽度越窄，孔径必定也越大。构成这个孔径的阵元之间必须以某种方式 隔开，以尽可能减少光栅瓣的产生，这可以通过阵元间相隔波长一半的距离来实现，如 果阵元间隔小于这个距离就可能看到光栅瓣。由阵元组成的换能器，理想的是阵列看起 来象是一个连续的阵元。对于高频，要求阵元之间的间距要比阵元本身小，因此阵元在 横行上要交错排列。旁瓣抑制通过单个阵元对整个阵列面作用的权重来实现。EM950的换 能器（见图2-6a ）是一个直径45cm,扇形角为

23、160的扇形圆柱体，由一个换能器完成 发射和接收，它由128个压电陶瓷条组成，每个陶瓷条由 5个阵元组成，这 5个阵元在 前后方向上有固定的权重， 在左右方向上间隔为 1.25 ，前后方向上的波束开角为 3.3 ， 波束中心与换能器面垂直。 图2-6a EM950换能器 EM3000的换能器（图2-6b ）由不同的阵列来完成接收和发射。发射阵列有56个环 形阵元。接收阵阵列有 80个栅条，每个栅条包含有 3个圆柱形阵元。接收栅条以半随机 方式交错排列，以使沿航迹方向的旁瓣级低于-20dB。 图2-6b EM3000换能器 根据上面的讨论，换能器阵列是设计发射一个沿航迹方向非常窄而垂直航迹方向宽

24、 的脉冲，垂直航迹方向的宽度取决于换能器和换能器阵元的波束图的宽度。一个典型的 多波束系统发射的波束宽度沿航迹方向小于3，而垂直航迹方向在天底两侧大于75 对没有自动纵摇稳定的系统发射的波束，接收波束的长度必须至少长于纵摇最大可能的 变化才不致丢失数据（图2-7 ），如EM3000沿航迹方向的接收波束宽度是25,Seabat8111 换能器接收波束沿航迹方向的波束宽度是15。在米氏交叉（Mill s Cross）这样结构 的换能器，要求接收阵列垂直于发射阵列（图 2-8 ）。 图2-8米氏交叉 2.2.3多波束的底部检测单元 一般，多波束测深系统的回波检测方式有两种，幅度检测和相位检测。当入射

25、角小 时，回波幅度高，持续时间短；当入射角变得十分大时，回波幅度低且持续时间长，但 波束间的相位差变大，故振幅检测对于中间波束传播时间的检测具有较高的精度，而对 边沿波束的检测精度较差，而相位检测正好相反。精密多波束测深系统利用相位检测用 于边沿波束检测，振幅检测用于中间波束检测，由系统取舍。这样，在可保证每个波束 检测精度的同时，又可保证整个波束的检测精度一致，从而达到波束旅行时高精度测定的目的1【2。 BDI （ Bearing Direction Indicator 方位指示）和WMT （ Weighted Mean Time 加权平 均时间）是两种不同的计算波束到达角和旅行时的方法。B

26、DI试图先定位每个波束回波的 方向即到达角，然后再精确计算旅行时；而 WMT先固定每个波束中心为到达角，然后再 精确计算出每个回波的旅行时3 4 5。为方便后续问题的展开，下面简要介绍它们的原理。 1、BDI处理方法 在一个发射接收周期（ping）内，波束形成采用的FFT处理方法中的数据可表示为矩 阵的形式，设一个周期内包含 M个时间片（time slice），每个时间片分别表示为 h, t2, tM,相对的时间周期起点为t。设有N个波束，每个时间片则可观测 N个幅度值，0 1, 0 2, 0 N,如图1所示。 A2N A3N AMN 波束角 回波幅度 1 AI1 A21A31 AM 1 1 A11 A21A31 AM 1 1 A11 A21A31 AM 1 时间片 totit2t3t m 一个发射接收周期 * 图1单ping