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

第二章 多波束勘测系统工作原理及构造多波束系统是70年代兴起、80年代中、末期又得到飞速发展的一项全新的海底地形精密勘测技术。它是目前爱好的焦点，由于它既有条带测深数据，又同步可获取反映底质属性的回波强度数据(Laurent Hellequin et al.,)。该技术采用广角度定向发射和多通道信息接受，获得水下高密度具有上百个波束的条幅式海底地形数据，彻底变化了老式测深技术概念，使测深原理、勘测措施、外围设备和数据解决技术诸方面都发生了巨大变化，大大提高了海底地形勘测的精度、辨别率和工作效率，实现了测深技术史上的一次革命性突破（李家彪等，）。多波束系统的工作原理与老式的单波束回声测深仪工作原理类似，都是根据声波在水下来回传播的时间与声速的乘积得到距离，从而得到水深。不同的是单波束测深仪一般采用较宽的发射波束（8左右）向船底垂直发射，声传播途径不会发生弯曲，来回的途径最短，能量衰减很小，通过对回声信号的幅度检测拟定信号来回传播的时间，再根据声波在水介质中的平均传播速度计算测量水深。在多波束系统中，换能器配备有一种或者多种换能器单元的阵列，通过控制不同单元的相位，形成多种具有不同指向角的波束，一般只发射一种波束而在接受时形成多种波束。除换能器天底波束外，外缘波束随着入射角的增长，波束在倾斜穿过水层时会发生折射，同步由于多波束沿航迹方向采用较窄的波束角而在垂直航迹方向采用较宽的覆盖角，要获得整个测幅上精确的水深和位置，必须要精确地懂得测量区域水柱的声速剖面和波束在发射和接受时船的姿态和船艏向。因此，多波束测深在系统构成和测量时比单波束测深仪要复杂得多（周兴华等，1999）。2.1 多波束勘测系统的工作原理2.1.1 单波束的形成2.1.1.1 发射阵和波束的形成一种单波束在水中发射后，是球形等幅度传播，因此方向上的声能相等。这种均匀传播称为各向同性传播（isotropic expansion），发射阵也叫各向同性源（isotropic source）。例如，一种小石头扔进池塘时就是这种状况，如图2.7所示。图2.7 波的各向同性传播显然，测深时是不能采用如此的声波的。采用发射基阵就可以产生各向异性的声波。下面简要论述它形成的原理。如果两个相邻的发射器发射相似的各向同性的声信号，声波图将互相重叠和干涉，如图2.8所示。两个波峰或者两个波谷之间的叠加会增强波的能量，波峰与波谷的叠加正好互相抵消，能量为零。一般地，相长干涉发生在距离每个发射器相等的点或者整波长处，而相消干涉发生在相距发射器半波长或者整波长加半波长处。显然，水听器需要放置在相长干涉处。图2.8 相长干涉和相消干涉（Constructive and Destructive Interference）一种典型的声纳，基阵的间距d（图2.8中1、2点的距离）是/2（半波长）。在这种状况下，相长和相消干涉发生时的点位处在最有利的角度（点位与基阵中心的连线与水平线的夹角），相长干涉：= 0, 180，相消干涉：= 90, 270，如图2.9所示。图2.9 两个发射器相距/2时的相长和相消干涉图2.10是两个发射器间距/2时的波束能量图（Beam Pattern），左边为平面图，右边为三维图，从图上可以清晰地看到能量的分布，不同的角度有不同的能量，这就是能量的指向性（directivity）。如果一种发射阵的能量分布在狭窄的角度中，就称该系统指向性高。真正的发射阵由多种发射器构成，有直线阵和圆形阵等。这里只讨论离散直线阵，其他阵列类似可以推导出。如图2.11所示，根据两个发射器的基阵可以推导出多种发射器构成的直线阵的波束图。图2.10 两个发射器间距/2时的波束能量图（Beam Pattern）图2.11 多基元线性基阵的波束图（Beam Pattern）图2.11中，能量最大的波束叫主瓣，侧边的某些小瓣是旁瓣,也是相长干涉的地方，引起了能量的泄漏。旁瓣还也许引起回波，对主瓣的回波产生干扰。旁瓣是不可避免的，可以通过加权的措施减少旁瓣的水平，但是加权后旁瓣水平值减少了，波束却展宽了。主瓣的中心轴叫最大响应轴（maximum response axisMRA），主瓣半功率处（相对于主瓣能量的-3db）角度的两倍就是波束角。发射器越多，基阵越长，则波束角越小，指向性就越高。设基阵的长度为D，则波束角= 50.6/D（2.36）可以看出，减小波长或者增大基阵的长度都可以提高波束的指向性。但是，基阵的长度不也许无限增大，而波长越小，在水中衰减得越快，因此指向性不也许无限提高。2.1.1.2 波束的指向（Beam Steering）换能器如何在指定的方向上发射或者接受声波，称为波束的指向。以水听器接受回波为例。如图2.12，当回波以方向达到接受基阵时，一方面在点3达到，另一方面为点2和点1，则在图2.12 夹角为的回波点2的回波比点3多旅行了距离Ad sin q，点1比点3的回波多旅行了距离B2d sin q，相应的增长的时间为T2=A/c =(d sin)/c （2.37）T1=B/c =(2d sin)/c （2.38）计算出偏移时间后，在基阵中作相应的调节，引入延时，使回波在基阵上正好构成相长干涉，这样就可以使主瓣在指定的方向上，如图2.13所示。 图2.13 引入延时后主瓣方向的偏移 图2.14 多波束的几何构成2.1.2 多波束的形成当接受波束发射出扇形波束后，接受波束按一定的间距（等距离或者等角度）与之相交，就形成了一种个在纵横向的窄波束脚印，如图2.14所示。设水听器共有N个基元，每个基元i记录的回波Si(t)的振幅为A(t),且S(t) = A(t)cos(2pft) （2.39）写成相位的形式为S(t) = A(t)cos(f(t) （2.40）或 （2.41）其中，f(t) = 2pft。多波束系统需测量回波S(t)和相位(t)，然后将模拟接受信号转换为数字信号，采用率一般在13ms之间。所有基元在采样点上的回波和相位值称为时间片（time slice）。在上节中，讨论了基元i相对于第一种基元的距离差，转换为相位差为 （2.42）由（2.41）和（2.42）得 （2.43）其中为第i个基元在角方向接受时的回波，则基阵接受的回波为 （2.44）其中为加权系数。如果规定在一种时间片(time slice)里，由N个基元形成M个指定方向的波束，用矩阵表达为 （2.45）其中，为接受角的波束时的第i个基元的相位差，为。为了在如此短的时间（ms级）完毕计算，必须采用某些迅速算法。这里，引入迅速傅立叶变换(FFT)。式（2.44）类似于傅立叶变换，设 （2.46）得 （2.47）由于k必须为整数，因此的取值受到一定的限制，如d、N在一定值时，同k的关系表2.1所示。表2.1 k和的相应值K012345()02.44.87.29.612.02.1.3 多波束脚印的归位波束脚印的归位是多波束数据解决的核心问题之一。多波束测量的最后成果是得到地理坐标系（或地方系）下的海底地形或者地物，由于多波束采用广角度定向发射、多阵列信号接受和多种波束形成解决等技术，为了更好的拟定波束的空间关系和波束脚印的空间位置，必须一方面定义多波束船体参照坐标系VFS，并根据船体坐标系同地理坐标系LLS之间的关系，将波束脚印的船体坐标转化到地理坐标系（或本地坐标系）和某一高程基准面下的平面坐标和高程。该过程即为波束脚印的归位。船体坐标系原点位于换能器中心，x轴指向航向，z轴垂直向下，y轴指向侧向，与x、z轴构成右手正交坐标系。地理坐标系原点为换能器中心，x轴指向地北子午线，y同x垂直指向东，z与x、y轴构成正交坐标系。归位需要的参数涉及船位、船姿、声速断面、波束达到角和来回程时间。归位过程涉及如下四个环节：（1） 姿态改正。（2） 船体坐标系下波束脚印位置的计算。（3） 波束脚印地理坐标的计算。（4） 海底点高程的计算。为以便波束脚印在船体坐标系下坐标的计算（声线跟踪），现作如下假设：（1） 换能器处在一种平均深度，静、动吃水觉得仅对深度有影响，而对平面坐标没有影响。（2） 觉得波束的往、返程途径重叠。（3） 对于高频发射系统，换能器的航向变化影响可以忽视。波束脚印船体坐标的计算需要用到三个参量，即垂直参照面下的波束达到角、传播时间和声速剖面。由于海水的作用，声束在海水中不是沿直线传播，而是在不同介质层的界面处发生折射，因此声束在海水中的传播途径为一折线。为了得到波束脚印的真实位置，就必须沿着波束的实际传播途径跟踪波束，该过程即为声线跟踪，通过声线跟踪得到波束脚印船体坐标的计算过程被称为声线弯曲改正。为了计算以便，对声速断面作如下假设：(1) 声速断面是精确的，无代表性误差。(2) 声速在波束形成的垂面内变化，不存在侧向变化。(3) 声速在海水中的传播特性遵循Snell法则。(4) 换能器的动吃水引起的声速剖面的变化对深度的计算可以忽视不计。根据上述讨论和假设，波束脚印的计算模型可体现为：Snell法则可描述为： （2.48）将波束的实际传播途径进行微分，则波束脚印在船体坐标系下的点位(x，y，z)可体现为： （2.49）zz: 深度 R:距离 q:波束角c: 声速 t:脉冲长度ln: 中心波束脚印长度ls: 边沿波束脚印长度 q0RRlnlgct/2图2.15 单个波束脚印坐标的计算x换能器(x0,z0)其一级近似式为： （2.50）更精确的公式见2.1.3。波束脚印的船体坐标系拟定后，下一步就可以转化为地理坐标。转换关系为： （2.51）式中，下脚g、gG别代表波束脚印的地理坐标、运用GPS拟定的船体地理坐标，R(h,r,p)为船体坐标系与地理坐标系的旋转关系，航向、横摇和纵摇是三个欧拉角。式(2.49)拟定的深度z仅为换能器面达到海底的垂直距离，测点的实际深度还应当考虑换能器的静吃水hss、动吃水hds、船体姿态对深度的影响ha，若潮位的变化htide是相对于某一深度基准面或者高程基准面拟定的，则波束脚印的高程为： （2.52）换能器的静吃水在测量前或换能器安装后被量定，作为一种常量输入到多波束的数据解决单元中；动吃水是由于船体的运动而产生的，它可通过姿态传感器中的Heaven参数拟定。船体姿态对波束脚印的地理坐标也有一定的影响，可通过姿态传感器的横摇r和纵摇p参数拟定。上述参数的测定及其对波束脚印平面位置和深度的补偿属于纯几何问题，武汉大学的赵建虎博士对此有具体的研究，本文不再赘述。2.2 多波束勘测系统的构成构造2.2.1 多波束的构成多波束系统重要由三个部分构成。第一部分是多波束的主系统，重要涉及换能器阵列，收发器和数据解决、显示和记录单元等；第二部分是辅助系统，涉及定位系统、船姿（横摇、纵摇、起伏和船艏向）测量传感器和测量水柱声速剖面的声速仪；第三部分是数据存贮和后解决系统，涉及数据解决计算机、数据存贮设备和绘图仪等（周兴华等，1999）。数据网GPS声速断面罗经姿态传感器换能器Transceiver操作和检测单元监控器导航监控器后解决实时数据解决工作站数据存储绘图仪打印机声纳影像记录数据存储图 2.1: Simrad EM950/1000 多波束声纳系统构成单元2.2.2 换能器的物理构成换能器是用来作为电声能量转换的重要器件。一般把电能转换成声能的器件称为发射换能器，把水下声能转换成电能的器件称为接受换能器（或水听器），许多积极声纳中采用同一的换能器兼作发射和接受（秦臻，1984）。水声换能器技术由3类构成，一是水声换能器材料；二是水声换能器设计；三是水声换能器制作。水声换能器材料方面波及材料的配备、生成及成型等3方面的技术。水声换能器设计方面波及性能设计及构造设计等两个方面的技术。水声换能器制作波及加工、安装及检测等3方面的技术。这3类水声换能器技术，是获取高性能、高质量水声换能器的保障。这3类水声换能器技术的发展，直接增进了水声换能器的发展，影响着声纳的发展。水深换能器就其所采用的换能器材料大体可分为两大类，一类是磁致伸缩换能器，一类是电致伸缩换能器。磁致伸缩材料涉及着具有水久磁矩的原子，它们以这种方式被耦合到晶格上，使之在加磁场时可以变化晶格的间距，从而相应地变化宏观样品的尺寸。由于磁致伸缩力与晶格变化有关，因此这个力可以和固体材料受机械变形所产生的弹性力相比较。磁致伸缩换能器把大的恒定磁场和较小的交变磁场叠加起来使之线性化，这种线性化的工作方式常常称为压磁，它最适合在声阻抗比较高的介质(如水)中工作（阎福旺等）。但是，此类材料在高频工作时有着高的涡流损耗和磁致损耗，以致电声转换效率减少，并且工作时，还需外加极化偏置。因此，目前多数声纳换能器采用了具有良好机电性能的电致伸缩材料（秦臻，1984）。电致伸缩材料是磁致伸缩材料的电学类比，它有永久性的电矩。电矩与晶格之间这样耦合，使之在外加电场时，尺寸发生变化。事实上这些极化了的材料，由于已经线性化，因而可以用描述压电材料的措施描述它。石英或许是最早实际用于换能器的材料，此后，被广泛用于实验工作、声能学和超声学方面。大概在1950年前后，可以实际应用的电致伸缩材料开始浮现，这些材料可以做成陶瓷，在极化后有极好的压电性质。钛酸钡是此类材料中一方面被广泛使用的材料。锆钛酸铅目前诸多应用中取代了钛酸钡（阎福旺等）。换能器是一种将电能转换成声能及逆过程的装置。不管使用的多波束是哪种类型，换能器是单个的或是一种阵列，都必须要在测量的海底形成一种声照射区（ensonified area），换能器通过发射声脉冲并接受产生的回波获得水深测量值。换能器的大小是根据波束的宽度需要来设计的。波束的宽度被定义为从天底到扩展波面半功率点的角距离，半功率的量值是-3dB。如图2-5所示Wells ,1996。主 瓣波束宽度旁瓣旁瓣换能器图2-5 波束宽度 每个多波束系统均有一种工作频率，换能器是由多种阵元构成的阵列。整个换能器阵列的大小是根据这样一种总的原则拟定的，即波束宽度与横跨孔径的波长成反比。反之，如果已选择了一种需要的波长，孔径的大小就可以用这个原则来拟定de Moustier, 1996。 如果： 需要的波束宽度 =2（=0.035rad ） 选择的频率为 F=100KHz， 声速C=1500m/s 波长 =0.015m 由总的原则： 孔径= 得到： 孔径=0.43m很明显波束宽度越窄，孔径必然也越大。构成这个孔径的阵元之间必须以某种方式隔开，以尽量减少光栅瓣的产生，这可以通过阵元间相隔波长一半的距离来实现，如果阵元间隔不不小于这个距离就也许看到光栅瓣。由阵元构成的换能器，抱负的是阵列看起来象是一种持续的阵元。对于高频，规定阵元之间的间距要比阵元自身小，因此阵元在横行上要交错排列。旁瓣克制通过单个阵元对整个阵列面作用的权重来实现。EM950的换能器（见图2-6a）是一种直径45cm，扇形角为160的扇形圆柱体，由一种换能器完毕发射和接受，它由128个压电陶瓷条构成，每个陶瓷条由5个阵元构成，这5个阵元在前后方向上有固定的权重，在左右方向上间隔为1.25，前后方向上的波束开角为3.3，波束中心与换能器面垂直。图2-6a EM950换能器EM3000的换能器（图2-6b）由不同的阵列来完毕接受和发射。发射阵列有56个环形阵元。接受阵阵列有80个栅条，每个栅条包具有3个圆柱形阵元。接受栅条以半随机方式交错排列，以使沿航迹方向的旁瓣级低于-20dB。图2-6b EM3000换能器根据上面的讨论，换能器阵列是设计发射一种沿航迹方向非常窄而垂直航迹方向宽的脉冲，垂直航迹方向的宽度取决于换能器和换能器阵元的波束图的宽度。一种典型的多波束系统发射的波束宽度沿航迹方向不不小于3，而垂直航迹方向在天底两侧不小于75。对没有自动纵摇稳定的系统发射的波束，接受波束的长度必须至少长于纵摇最大也许的变化才不致丢失数据（图2-7），如EM3000沿航迹方向的接受波束宽度是25，Seabat8111换能器接受波束沿航迹方向的波束宽度是15。在米氏交叉（Mills Cross）这样构造的换能器，规定接受阵列垂直于发射阵列（图2-8）。 最大向上纵摇时声照射区 最大向下纵摇时声照射区规定的最小接受波束长度图2-7 接受波束的最小宽度发射阵列接受阵列沿航迹方向图2-8 米氏交叉2.2.3 多波束的底部检测单元一般，多波束测深系统的回波检测方式有两种，幅度检测和相位检测。当入射角小时，回波幅度高，持续时间短；当入射角变得十分大时，回波幅度低且持续时间长，但波束间的相位差变大，故振幅检测对于中间波束传播时间的检测具有较高的精度，而对边沿波束的检测精度较差，而相位检测正好相反。精密多波束测深系统运用相位检测用于边沿波束检测，振幅检测用于中间波束检测，由系统取舍。这样，在可保证每个波束检测精度的同步，又可保证整个波束的检测精度一致，从而达到波束旅行时高精度测定的目的1 2。BDI（Bearing Direction Indicator方位批示）和WMT（Weighted Mean Time加权平均时间）是两种不同的计算波束达到角和旅行时的措施。BDI试图先定位每个波束回波的方向即达到角，然后再精确计算旅行时；而WMT先固定每个波束中心为达到角，然后再精确计算出每个回波的旅行时3 4 5。为以便后续问题的展开，下面简要简介它们的原理。1、BDI解决措施在一种发射接受周期（ping）内，波束形成采用的FFT解决措施中的数据可表达为矩阵的形式，设一种周期内涉及M个时间片（time slice），每个时间片分别表达为t1，t2，tM，相对的时间周期起点为t0。设有N个波束，每个时间片则可观测N个幅度值，1，2，N，如图1所示。时间片波束角回波幅度一种发射接受周期图1 单ping回波强度矩阵为了避免旁瓣或海底混响等其他因素引起的回波噪声，对于每个时间片，使用一种动态阈值，也叫检测门限，高于阈值的幅度接受，低于阈值的幅度视为噪声。阈值由旁瓣水平和噪声水平决定，每个时间片均不同3。将高于阈值的幅度值进行拟合，求得幅度的极大值，同步记下相应的角度值，就是该时间片相应的DOA，对于如图2所示。这样，每个幅度极大值相应此外两个数据：角度和时间。这3个数反映了波束发射到接受的过程，以击中（hit）表达，即每个击中用这3个数据表达。图2 单时间片中精确角度的估计如图3（取自文献3）所示，在一种发射接受周期内，将所有的击中以角度和时间表达出来，在每个波束范畴内，计算平均时间tmean和方差，保存与tmean相差1倍内的击中，最后根据幅度加权计算出每个波束的DOA和TOA3，图中表达为和。为了问题的简化，本文没有考虑姿态补偿。 （1） （2） 图3 幅度加权平均的DOA和TOA 图4 参与计算的幅度值2、WMT解决措施WMT一方面固定一系列的波束角，这些波束角已经进行了姿态补偿。然后计算出每个波束角的精确的TOA。固定的波束角往往并不在波束轴中心。在每个时间片，都可得到每个固定角的幅度值，在整个发射接受周期内，可得到每个固定角的时序观测的幅度值，使用开始、结束门限和动态阈值，就可得到最后计算TOA需要的幅度值，然后按幅度加权就可得到精确的TOA，如图4（取自文献3）所示，图中粗线即为计算TOA保存的幅度值，处在开始门、结束门外和动态阈值下的均被忽视，最后的TOA计算同式（2）。具体是选择BDI还是WMT解决措施，要根据具体状况看待。对于镜面反射，反射强度大，回波脉冲尖，持续时间短，此时BDI就不太合用，由于很难精确地计算回波的角度，因而采用WMT措施；相反，在非镜面区，反射强度小，回波脉冲平缓，持续时间长，此时WMT就不太合用，由于很难精确地计算回波时间，因而采用BDI措施3 4 5。