多波束勘测系统工作原理与结构

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

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

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

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

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

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

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

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

9、t )cos(2p ft )(2.40 )S(t) = A(t)cos(f( t)(2.41 )S(t) = A(t)ej：其中，f( t) = 2p ft 。多波束系统需测量回波 S（t）和相位0 （t），然后将模拟接收信号转换为数字信号，采 用率一般在13ms之间。所有基元在采样点上的回波和相位值称为时间片（time slice ）。(2.42 )在上节中，讨论了基元i相对于第一个基元的距离差，转换为相位差为 id sin 二由（2.41 ）和（2.42 ）得Bi（小扁心肚 2%玉s= S&吩(2.43 )其中BiL）为第i个基元在角B方向接收时的回波，则基阵接收的回波为N 二j (竺扁

10、sin B(2.44 )B(R 八sQei =9（time slice ）里，由N个基元形成M个指定方向的波束，用矩阵表示为DuD12D1NS1 )Bd)D21D22D2NXS2lD m 1DM2Dmn JSn丿其中S为加权系数。如果要求在一个时间片(2.45 )j （绍d sin 日）一其中，Dij为接收角可的波束时的第i个基元的相位差，为 e 。为了在如此短的时间（ms 完成计算，必须采用一些快速算法。这里，引入快速傅立叶变换（FFT）。式（2.44 ）类似于傅立叶变换，设(2.46 )(2.47 )BQ)，hieFi=0由于k必须为整数，所以 入的取值受到一定的限制，如 d、N在一定值时

11、， 耳同k的关 系表2.1所示。表2.1 k和&的对应值K0123450 ( ) 0 2.4 4.8 7.2 9.6 12.0 2.1.3多波束脚印的归位波束脚印的归位是多波束数据处理的关键问题之一。多波束测量的最终成果是得到 地理坐标系（或地方系）下的海底地形或者地物，由于多波束采用广角度定向发射、多 阵列信号接收和多个波束形成处理等技术，为了更好的确定波束的空间关系和波束脚印 的空间位置，必须首先定义多波束船体参考坐标系VFS并根据船体坐标系同地理坐标系LLS之间的关系，将波束脚印的船体坐标转化到地理坐标系（或当地坐标系）和某一高程基准面下的平面坐标和高程。该过程即为波束脚印的归位。船体坐

12、标系原点位于换能器 中心，x轴指向航向，z轴垂直向下，y轴指向侧向，与x、z轴构成右手正交坐标系。地 理坐标系原点为换能器中心，x轴指向地北子午线，y同x垂直指向东，z与x、y轴构成正交坐标系。归位需要的参数包括船位、船姿、声速断面、波束到达角和往返程时间。归位过程 包括如下四个步骤：（1）姿态改正。（2）船体坐标系下波束脚印位置的计算。（3）波束脚印地理坐标的计算。（4）海底点高程的计算。为方便波束脚印在船体坐标系下坐标的计算（声线跟踪），现作如下假设：（1）换能器处于一个平均深度，静、动吃水认为仅对深度有影响，而对平面坐标没有影响。（2）认为波束的往、返程路径重合。（3）对于高频发射系统，

13、换能器的航向变化影响可以忽略。波束脚印船体坐标的计算需要用到三个参量，即垂直参考面下的波束到达角、传播 时间和声速剖面。由于海水的作用，声束在海水中不是沿直线传播，而是在不同介质层 的界面处发生折射，因此声束在海水中的传播路径为一折线。为了得到波束脚印的真实 位置，就必须沿着波束的实际传播路径跟踪波束，该过程即为声线跟踪，通过声线跟踪 得到波束脚印船体坐标的计算过程被称为声线弯曲改正。为了计算方便，对声速断面作 如下假设：（1）声速断面是精确的，无代表性误差。（2）声速在波束形成的垂面内变化，不存在侧向变化。（3）声速在海水中的传播特性遵循Snell法则。 换能器的动吃水引起的声速剖面的变化对

14、深度的计算可以忽略不计。根据上述讨论和假设，波束脚印的计算模型可表达为：Snell法则可描述为:C0C1(2.48 )(x, y，z)可表将波束的实际传播路径进行微分，则波束脚印在船体坐标系下的点位 达为：z = z0 亠 iC(z) cos(v(z)dzx=XoC(z)s in (v(z)dz(2.49)z =ZoCTp cos%2旦siE2(2.50)更精确的公式见 2.1.3 。波束脚印的船体坐标系确定后，下一步就可以转化为地理坐标。转换关系为:X = X0R(h,r,p)X_y g_y0 gG一y(2.51)式中，下脚g、G别代表波束脚印的地理坐标、利用GPS确定的船体地理坐标，R(h

15、,r,p)为船体坐标系与地理坐标系的旋转关系，航向、横摇和纵摇是三个欧拉角。式(2.49)确定的深度z仅为换能器面到达海底的垂直距离，测点的实际深度还应该 考虑换能器的静吃水hss、动吃水hds、船体姿态对深度的影响ha，若潮位的变化htide是相对于某一深度基准面或者高程基准面确定的，则波束脚印的高程为：Zg 二 Z hide hss hds - ha（ 2.52）换能器的静吃水在测量前或换能器安装后被量定，作为一个常量输入到多波束的数据处 理单元中；动吃水是由于船体的运动而产生的，它可通过姿态传感器中的Heaven参数确定。船体姿态对波束脚印的地理坐标也有一定的影响，可通过姿态传感器的横摇

16、r和纵摇p参数确定。上述参数的测定及其对波束脚印平面位置和深度的补偿属于纯几何问题， 武汉大学的赵建虎博士对此有详细的研究，本文不再赘述。 2.2多波束勘测系统的组成结构2.2.1多波束的组成多波束系统主要由三个部分组成。第一部分是多波束的主系统，主要包括换能器阵 列，收发器和数据处理、显示和记录单元等；第二部分是辅助系统，包括定位系统、船 姿（横摇、纵摇、起伏和船艏向）测量传感器和测量水柱声速剖面的声速仪；第三部分 是数据存贮和后处理系统，包括数据处理计算机、数据存贮设备和绘图仪等（周兴华等， 1999）。监控器GPS操作和检测单元声速断罗经ransceiver姿态传感器换能器导航监控器实时

17、数据处理工作站数据存储声纳影像记录后处理数据存储绘图仪打印机图2.1: Simrad EM950/1000 多波束声纳系统组成单元2.2.2换能器的物理构成换能器是用来作为电声能量转换的重要器件。通常把电能转换成声能的器件称为发 射换能器，把水下声能转换成电能的器件称为接收换能器（或水听器），许多主动声纳中采用同一的换能器兼作发射和接收（秦臻，1984）。水声换能器技术由 3类组成，一是水声换能器材料；二是水声换能器设计；三是水 声换能器制作。水声换能器材料方面涉及材料的配置、生成及成型等3方面的技术。水声换能器设计方面涉及性能设计及结构设计等两个方面的技术。水声换能器制作涉及加 工、安装及检

18、测等 3方面的技术。这3类水声换能器技术，是获取高性能、高质量水声 换能器的保障。这 3类水声换能器技术的发展，直接促进了水声换能器的发展，影响着 声纳的发展。水深换能器就其所采用的换能器材料大致可分为两大类，一类是磁致伸缩换能器， 一类是电致伸缩换能器。磁致伸缩材料包含着具有水久磁矩的原子，它们以这种方式被 耦合到晶格上，使之在加磁场时能够改变晶格的间距，从而相应地改变宏观样品的尺寸。 由于磁致伸缩力与晶格变化有关，所以这个力可以和固体材料受机械变形所产生的弹性 力相比较。磁致伸缩换能器把大的恒定磁场和较小的交变磁场叠加起来使之线性化，这 种线性化的工作方式常常称为压磁，它最适合在声阻抗比较

19、高的介质（如水）中工作（阎福旺等）。不过，这类材料在高频工作时有着高的涡流损耗和磁致损耗，以致电声转换效率降低，而且工作时，还需外加极化偏置。因此，目前多数声纳换能器采用了具有良好 机电性能的电致伸缩材料(秦臻，1984)。电致伸缩材料是磁致伸缩材料的电学类比，它有永久性的电矩。电矩与晶格之间这样耦合，使之在外加电场时，尺寸发生变化。实际 上这些极化了的材料，由于已经线性化，因而能够用描述压电材料的方法描述它。石英 或许是最早实际用于换能器的材料，此后，被广泛用于实验工作、声能学和超声学方面。 大约在1950年前后，能够实际应用的电致伸缩材料开始出现，这些材料能够做成陶瓷， 在极化后有极好的压

20、电性质。钛酸钡是这类材料中首先被广泛使用的材料。锆钛酸铅现 在很多应用中取代了钛酸钡(阎福旺等)。换能器是一个将电能转换成声能及逆过程的装置。不管使用的多波束是哪种类型，换能器是单个的或是一个阵列，都必须要在测量的海底形成一个声照射区(en so nifiedarea)，换能器通过发射声脉冲并接收产生的回波获得水深测量值。换能器的大小是根据 波束的宽度需要来设计的。波束的宽度被定义为从天底到扩展波面半功率点的角距离，半功率的量值是-3dB。如图2-5所示Wells ,1996。换能器波束宽度图2-5波束宽度每个多波束系统都有一个工作频率，换能器是由多个阵元组成的阵列。整个换能器 阵列的大小是根

21、据这样一个总的原则确定的，即波束宽度与横跨孔径的波长成反比。反 之，如果已选择了一个需要的波长，孔径的大小就可以用这个原则来确定de Moustier,1996。如果： 需要的波束宽度 3 =2( =0.035rad )选择的频率为F=100KHz , 声速C=1500m/s波长入=匸=0.015m由总的原则：孔径=-p入得到： 孔径=0.43m很明显波束宽度越窄，孔径必定也越大。构成这个孔径的阵元之间必须以某种方式隔开，以尽可能减少光栅瓣的产生，这可以通过阵元间相隔波长一半的距离来实现，如 果阵元间隔小于这个距离就可能看到光栅瓣。由阵元组成的换能器，理想的是阵列看起 来象是一个连续的阵元。对

22、于高频，要求阵元之间的间距要比阵元本身小，因此阵元在 横行上要交错排列。旁瓣抑制通过单个阵元对整个阵列面作用的权重来实现。EM950的换能器（见图2-6a ）是一个直径45cm,扇形角为160的扇形圆柱体，由一个换能器完成 发射和接收，它由128个压电陶瓷条组成，每个陶瓷条由 5个阵元组成，这 5个阵元在 前后方向上有固定的权重， 在左右方向上间隔为 1.25 ，前后方向上的波束开角为 3.3 ， 波束中心与换能器面垂直。图2-6a EM950换能器EM3000的换能器（图2-6b ）由不同的阵列来完成接收和发射。发射阵列有56个环形阵元。接收阵阵列有 80个栅条，每个栅条包含有 3个圆柱形阵

23、元。接收栅条以半随机 方式交错排列，以使沿航迹方向的旁瓣级低于-20dB。图2-6b EM3000换能器根据上面的讨论，换能器阵列是设计发射一个沿航迹方向非常窄而垂直航迹方向宽 的脉冲，垂直航迹方向的宽度取决于换能器和换能器阵元的波束图的宽度。一个典型的 多波束系统发射的波束宽度沿航迹方向小于3，而垂直航迹方向在天底两侧大于75 对没有自动纵摇稳定的系统发射的波束，接收波束的长度必须至少长于纵摇最大可能的 变化才不致丢失数据（图2-7 ），如EM3000沿航迹方向的接收波束宽度是25,Seabat8111换能器接收波束沿航迹方向的波束宽度是15。在米氏交叉（Mill s Cross）这样结构的

24、换能器，要求接收阵列垂直于发射阵列（图 2-8 ）。最大向上纵摇时声照射区要求的最小接收 波束长度图2-8米氏交叉2.2.3多波束的底部检测单元一般，多波束测深系统的回波检测方式有两种，幅度检测和相位检测。当入射角小 时，回波幅度高，持续时间短；当入射角变得十分大时，回波幅度低且持续时间长，但 波束间的相位差变大，故振幅检测对于中间波束传播时间的检测具有较高的精度，而对 边沿波束的检测精度较差，而相位检测正好相反。精密多波束测深系统利用相位检测用 于边沿波束检测，振幅检测用于中间波束检测，由系统取舍。这样，在可保证每个波束 检测精度的同时，又可保证整个波束的检测精度一致，从而达到波束旅行时高精

25、度测定的目的1 2BDI （Bearing Direction Indicator方位指示）和 WM（ Weighted Mean Time 加权平均时间）是两种不同的计算波束到达角和旅行时的方法。BDI试图先定位每个波束回波的方向即到达角，然后再精确计算旅行时；而WM先固定每个波束中心为到达角，然后再精确计算出每个回波的旅行时3 4 5。为方便后续问题的展开，下面简要介绍它们的原理。1、BDI处理方法在一个发射接收周期（ping ）内，波束形成采用的FFT处理方法中的数据可表示为矩 阵的形式，设一个周期内包含M个时间片（time slice ），每个时间片分别表示为t1,t2，, tM,相对的时间周期起点为t。设有N个波束，每个时间片则可观测 N个幅度值，0 1, 0 2, 0 N,如图1所示。波束角回波幅度AI1A21A31AMAI1A21A31AMA11A21A31AMA1NA2NA3NAMNItoTt1Tt2Tt3Tt