多波束水深测量数据处理及思考.docx

多波束水深测量数据处理及思考.docx

《多波束水深测量数据处理及思考.docx》由会员分享，可在线阅读，更多相关《多波束水深测量数据处理及思考.docx（21页珍藏版）》请在冰点文库上搜索。

多波束水深测量数据处理及思考

多波束水深测量数据处理及思考

顾顺隆（上海海事局海测大队）

摘要：

根据多波束测深及其辅助传感器测量原理，对多波束数据处理方法进行了探讨，并对各种误差影响进行了思索。

关键词：

多波束坐标系传感器系统偏差归算误差

1前言

随着科技的发展，水深测量方法已经从测杆、测绳测深发展到超声波测深；从电火花测深记录纸模拟输出到数字化输出、断面图像文件（无纸化）输出；从单波束线测深发展到多波束面测深；从六分仪定位发展到GPS定位；从人工测量方式发展到自动化数据采集。

伴随着测深技术的进步，我们的测量能力和测绘成果的质量得到了空前提高。

多波束测深仪是单波束测深仪的发展，其显著特点是能一次发射和接收一列波束，在测量船速等控制得当的情况下，多波束测深仪可以对一个区域进行全覆盖的面测深，这是对单波束测深技术的巨大改进。

多波束测深技术的应用，不仅提高了外业测量的作业效率，同时也提高了野外探测能力；不仅为水深测量提供了方便，也为沉船或沉物的寻找、暗礁或浅点的探查提供了方便。

随着多波束测深技术的全面推广应用，大家赞叹着技术进步带来的便利，但同时也出现了一些令人疑惑甚至令人担忧的问题，需要从实践和理论两方面进行分析研究，如，单波束测深仪或旁侧声呐测到的某些目标，多波束却未能测出。

问题的存在，引人思考：

面对技术进步，不仅要引进新技术，还要能够及时消化和吸收；不仅提高设备的技术水平，还要同时提高作业人员的技术素养。

以下是对多波束数据处理等的初略探讨。

2多波束系统的坐标系

多波束测深系统是多传感器组成的综合系统，除了多波束测深仪和定位仪器外，还包括测定船舶航向的电罗经、测定纵摇横摇的姿态传感器及测定上下起伏的涌浪滤波器等辅助测量设备。

只有当它们相互匹配时，才能正常开展测量工作。

为了方便叙述及公式推导，本文定义和使用如下的坐标系统：

2·1测量船坐标系（XbYbZb）

右手系，往船艏方向看，向左为X轴正向，与船艏艉线平行为Y轴正向，垂直向下为Z轴正向，坐标原点一般选在多波束探头发射中心。

多波束探头、姿态传感器、电罗经、GPS等的安装位置及其坐标系均依测量船坐标系确定。

（见图1）

2·2多波束探头坐标系（XcYcZc）

与多波束测深扇面垂直，并指向船艏为Y轴正向，向左为X轴正向，中心波束方向为Z轴正向，坐标原点在多波束探头发射中心，三轴与测量船坐标系对应轴平行。

2·3姿态传感器坐标系（XdYdZd）

三轴与测量船坐标系对应轴同向平行。

2·4电罗经坐标系（XeYeZe）

三轴与测量船坐标系对应轴同向平行。

2·5水平坐标系（XaYaZa）

与平静时的测量船坐标系重合，XaYa平面与当地水平面平行。

水深点纵摇改正、横摇改正和航向倾斜改正等，均参照该坐标系进行。

2·6测量坐标系（XY）

一般使用高斯坐标系，多波束水深经航向归算后，最后纳入该坐标系。

3横摇与纵摇改正

多波束测深系统进行水深测量时，因受风浪的影响，测量船不免产生摇摆，导致测量船坐标系的倾斜，使测得的水深出现系统偏移，从而需要同时进行船舶姿态测量，并对水深进行倾斜改正。

姿态传感器可以测量船舶的横摇（roll）与纵摇（pitch）（如图2），依此可进行横摇改正与纵摇改正。

改正时，测量船坐标系（XbYbZb）绕直线l旋转ρ即可与水平坐标系（XaYaZa）重合。

所以有

（1）

其中，

（2）

α为横摇，β为纵摇。

声线处理时，该公式不很方便，为此，下面给出另一种形式的转换公式：

首先，水平坐标系（XaYaZa）绕Za轴旋转偏角ε，得坐标系（X'aY'aZ'a），并要求Y'a轴与Yb轴在XaYa上的投影重合（见图3），有

（3）

分别绕Y'（Yb）和X'a轴旋转，得

（4）

其中，

（5）

α为横摇，β为纵摇。

偏转角ε可按下面的公式直接求取：

（6）

其中，ρ和κ按

（2）式求取（参见图4）。

4航向归算

多波束系统测向一般采用光纤罗经或陀螺罗经，它们都依据地球自转进行航向测定，但就测量原理来讲，有着本质的不同。

在航向归算前，先要弄清它们测得的是什么量，需要进行怎样的变换处理。

4·1光纤罗经

光纤陀螺原理基于光学的萨格钠克（Sagnac）效应：

当光束在一个环形的通道中前进时，如果环形通道本身具有一个转动速度，则在不同的前进方向上，光学环路的光程相对于环路在静止时的光程都会产生变化。

利用这种光程的变化，光纤陀螺可以测量光纤环路的转动速度。

光纤罗经安装三个相互垂直的光纤陀螺，可分别测出地球自转角速度在这三个垂直方向的分量值，如果三个光纤陀螺仪的旋转轴分别与测量船坐标系的三根轴相对应平行，则由所测得的角速度分量便可确定船舶航向。

公式如下：

（7）

而

（8）

其中：

ωe为地球自转的角速度，φ为当地纬度，ω1为ωe的水平方向分量，ω1为ωe的垂直方向分量。

这是船舶平稳水平时的航向计算公式，当船舶倾斜时，三个光纤陀螺测出的是倾斜坐标系三轴的角速度分量值，一般可借助于光纤罗经仪自带的姿态传感器所测出的船舶的横摇、纵摇，利用公式

（1）进行转换，得到水平坐标系的Xa、Ya轴方向的两水平分量ωx和ωy，再由（7）式求得水平坐标系的Ya轴的真方位A。

4·2陀螺罗经

陀螺罗经是基于物理学角动量原理（或称动量矩原理）的一种导航仪器。

根据角动量原理，当高速旋转的陀螺转子，没有外力矩作用时，转动轴将保持方向不变，这就是陀螺仪的定轴性；而有外力矩作用时，转动轴将朝外力矩的方向转动，这就是陀螺仪的进动性。

陀螺罗经利用了陀螺仪的定轴性和进动性，由地球自转时产生相应的控制力矩（如重力力矩），使陀螺主轴自动找北；外加使用阻尼力矩，使陀螺主轴稳定指北。

从结构上看，陀螺罗经一般由灵敏部分，随动部分和固定部分三部分组成。

灵敏部分为陀螺球，内有陀螺转子、陀螺电机等，其作用是指北。

随动部分包括随动支架、方位齿轮等，随动部分的作用之一是跟踪陀螺球，使主罗经刻度盘零度线与陀螺主轴方向保持一致；对于数码刻度盘，陀螺球的航向检测不用传统的机械电器式的传感器、齿轮、同步电机等，而是采用数码检测信号。

固定部分包括罗经桌或罗经箱体等，其作用除了支承陀螺球、随动部分等外，是要利用其上设置的罗经基线获取航向。

从对不同类型的陀螺罗经的结构分析，不难看出，船舶航向是陀螺主轴和罗经基线在主罗经刻度盘平面上投影线间的夹角。

陀螺罗经工作时，主罗经刻度盘的可以绕垂直轴旋转，而垂直轴安装于固定部分位置相对不变，因此，当船舶倾斜时，主罗经刻度盘随之倾斜。

于是可得结论：

陀螺罗经测定的是倾斜面上的航向。

如果，陀螺罗经未进行倾斜修正，必要时可按以下面的公式进行处理。

在水平坐标系（X'aY'aZ'a）中，陀螺主轴向量为{x'a，y'a，z'a}，在倾斜坐标系中，陀螺主轴向量为{xb，yb，zb}，令z'a=0，xb=sinA'，yb=cosA'，根据公式（4）可得，

（9）

其中，A'是陀螺罗经未经倾斜修正的航向，其余参数同（4）式，求得x'a，y'a后，再由下式

（10）

求得Y'a轴的真方位A1，Ya轴的真方位

（11）

偏转角ε按（6）式计算。

航向倾斜改正估算，假定只存在横摇，即β=0，由（9）式不难得到：

（12）

此时，航向倾斜改正最大值出现在45°、135°、225°或315°航向上，α=10º，倾斜改正约有°；α=15º，倾斜改正约有1°。

5声线计算

多波束声线弯曲计算采用二维模式，即声速只在垂直深度上发生变化，而在水平方向上均匀。

5·1公式

对于连续变化的声速，有

（13）

其中，t为计算点单程旅行时的时间，θ是计算点声线切线与垂线之间的夹角，v为计算点声速，D为换能器到海底的水深，X为侧向中心距。

对于离散的声速，有

（14）

其中，Δt为各层单程旅行时的时间，θ是各层入射角，v为各层声速，Δz是各层厚，D为换能器到海底的水深，X为侧向中心距。

θ值按斯涅耳定律计算

（15）

θ1是波束入射角，v1为换能器表面声速。

5·2水深及侧向中心距表

多波束水深测量每Ping波束少则几十，多则上千，为减少重复计算，可根据垂线声速预先计算不同入射角和不同旅行时间所对应的水深及中心偏距，并制成表。

再利用该表，并根据各波束实际的入射角和旅行时间插值计算出水深及偏距。

插值计算一般采用线性内插，如图5是对应（ti，ti+1）及（θi，θi+1）的四个水深和偏距，令

，

（16）

有

（17）

（18）

5·3线性内插误差估算

1.假定垂线声速不变，并取λ=，φ=，得

（19）

其中，l为波束斜距，δθ=（θi+1-θi）/2。

对于边缘波束，取l=200米，δθ=1°，θ=75°，有Vx=0.03米，VD=0.01米；对于近中央波束，取l=200米，δθ=1°，θ=1°，有Vx=0米，VD=0.03米。

2.

一般垂线声速是变化的，假定在声波旅行的某个时段（∆t）范围内，声速变化为∆v，则该范围内最大声速与最小声速所对应的声线长度之差的公式是Δl<Δv×Δt（即声线长度的变化）。

取∆t=2毫秒，∆v=15米/秒，则在步长中点处声线长度的变化Δl<Δv×Δt=厘米，即对垂直波束水深的影响小于厘米。

由上可见，当入射角步长取2°，旅行时间取2毫秒时，线性内插对偏距和水深的影响不大于5厘米。

5·4空间波束入射角计算

波束入射角是波束中心线与水平面法线（垂直轴Za）之间的夹角；波束单元入射角（换能器各波束单元的入射角）是波束中心线与换能器垂直轴之间的夹角，Y轴方向逆时针为正。

由于测量时存在横摇和纵摇，一般波束单元入射角不能直接用作波束入射角，而需要经过必要的换算。

如图6不难得到，波束入射角可按下式计算，

（20）

而

（21）

其中，γ为波束单元入射角，α'、β'、b'按（5）式计算。

根据波束入射角θ及旅行时间t，利用水深和中心偏距表即可求得水深D和中心偏距X，再按以下公式，

（22）

及（3）式便得到水平坐标系坐标。

波束入射角误差对水深点的影响，可按下式估算，

（23）

其中，l为声线斜距，θ为波束入射角，δθ为波束入射角误差。

6系统偏差改正

多波束测深系统各传感器坐标系的对应坐标轴应当相互平行，但由于安装偏差的存在，这点很难保证。

电罗经和姿态传感器的安装偏差尚可通过静态测定进行测定或校正，多波束探头的安装偏差却很难进行这样的静态测定。

如果不对系统安装偏差进行有效处理，将极大影响多波束测深系统的测量质量。

为此，根据多波束测深系统的特点，设计了横摇综合偏差、纵摇综合偏差和艏向综合偏差等偏差测定方法。

6·1横摇综合偏差

横摇综合偏差测定宜选择在平坦水域进行，因为在平坦水域，航向偏差只影响到水深点位；纵摇偏差尽管还影响到水深，但是对各水深点是等幅度的。

在航向偏差和纵摇偏差的影响下，水深还是平坦的。

因此，在平坦水域可突现横摇偏差的影响，测定后应先行改正。

测定采用同测线往返测量的方式。

6·2纵摇综合偏差

纵摇综合偏差测定宜选择在水深有明显变化的水域进行，如，有一定坡度的斜坡（航向垂直于斜坡）。

测定也采用同测线往返测量的方式，通过比较同水深点航向上的点位变化，确定纵摇偏差。

采用同测线往返测量的方式时，航向偏差使多波束探头扇面偏转一角度，如无纵摇和横摇的影响，在同一船位处，两Ping平行，将测得相同的水深。

因此，采用同测线往返测量可突出纵摇偏差的影响，宜在航向偏差前进行测定。

6·3艏向综合偏差

艏向综合偏差测定也宜选择在水

深有明显变化的水域进行，有明显的目标点更佳。

测定采用两平行测线同向测量的方式，通过比较重叠带内同水深点（目标点）航向上的点位变化，确定艏向偏差。

横摇综合偏差和纵摇综合偏差包括多波束探头和运动姿态传感器安装偏差的共同影响；艏向综合偏差则包括多波束探头和电罗经安装偏差的共同影响。

除了进行上述几项偏差测定外，对于采用GPS等设备定位的，还应先进行定位时延测定。

综合偏差改正适用于姿态传感器与光纤罗经一体化的设备，一般姿态传感器与光纤罗经两者的坐标系具有很好的一致性，即对应轴平行。

综合偏差的测定，实际上，相当于确定了当α=0、β=0时，多波束探头坐标系（XcYcZc）相对于水平坐标系（XaYaZa）的偏转，如图7所示，也建立了（XcYcZc）与（XdYdZd）及（XeYeZe）的空间关系。

当α、β≠0时，姿态传感器坐标系（XdYdZd）相对于水平坐标系（XaYaZa）偏转，多波束探头坐标系（XcYcZc）随之进一步偏转。

采用下面的方法，能十分有效地对横摇综合偏差等进行改正，把多波束探头坐标归算到水平坐标。

在此转换中，（XdYdZd）及（XeYeZe）担当了（XaYaZa）与（XcYcZc）的中间过渡。

求出α、β≠0时Yc轴在XaYa上的投影

取

，

，

。

其中，

（24）

利用

（1）式，得到（xa，ya，za），并由下式得坐标转偏角ε'。

（25）

求出α、β≠0时Zc轴相对于（XaYaZa）的倾斜

取xb=δa，yb=δb，zb=1。

由

（1）式得到（xa，ya，za），并由下式计算

（26）

（27）

由（20）、（21）、（22）各式进行声线计算，得（x'a，y'a，z'a），坐标偏转角ε'转换，可在以后的水深点航向转换时一并处理。

对于姿态传感器与电罗经分离的设备来说，应先对姿态传感器和电罗经坐标系偏差进行静态测定，再对多波束探头进行横摇偏差、纵摇偏差和航向偏差动态测定。

这些测定应以测量船坐标系（XbYbZb）为参照，并采用响应的方法进行数据处理。

在此，值得注意的是如何确定（XdYdZd）与（XeYeZe）的空间关系，即姿态传感器坐标系Yd轴的方向如何？

罗经坐标系是否倾斜？

7水深点测量坐标转换

对于一体化的设备，通过静态测定，对α、β及A的观测值取算术平均，即得α、β及A的系统偏差，分别为δαd、δβd及δAe（船航向与A均值之差），如图8所示，与综合偏差相类似，但处理方式有所不同。

7·1定位仪点位归算

3.根据设备安装文件，对定位仪点位进行α、β及A的系统偏差改正，取xb=Xps，yb=Yps，zb=Zps，α=δαd，β=δβd，ε=-δAe，其中，（Xps，Yps，Zps）为定位仪在测量船坐标系中坐标（相对于多波束探头坐标系（XcYcZc）原点），由（3）、（4）和（5）式得（xa，ya，za）；

4.再进行倾斜转换，取xb=xa，yb=ya，zb=za，α=αd，β=βd，其中，（xa，ya，za）为上段计算所得，αd、βd为横摇和为纵摇观测值。

由

（1）式得（xa，ya，za）；

5.再利用下式，对以上转换结果进行航向变换，

（28）

其中，

（29）

A为真航向观测值，γ为高斯坐标收敛角，（x'，y'）为经航向旋转后的定位仪点坐标。

7·2水深点坐标转换

经过声线归算后，得到各水深点在水平坐标系的点位（x'a，y'a，z'a），利用下式对其转换，

（30）

其中，

（31）

A为真航向观测值，γ为高斯坐标收敛角，ε'为按（25）计算的坐标偏转角，（x'，y'）为经航向旋转后的定位仪点坐标，（x，y）为水深点高斯坐标，（X，Y）为测量船定位坐标值。

7·3起伏和潮位改正

经过声线归算后，得到是各水深点在探头以下的水深za，另外，测量时除了潮位变化外，还有测量船的上下起伏。

因此，需要进行以下处理，

（32）

其中，d为深度基准面下水深，za为声线归算后水深，h1为船舶吃水（包括静态吃水和动态吃水），h为潮位改正，hv为涌浪改正。

对于涌浪改正，一个明显的感觉：

在船的不同位置，上下起伏是不同的，船中后部较小，船头较大。

这种不同，不仅引起了船舶的倾斜，同时也使得在船的不同位置有不同的涌浪值，因此，需要对涌浪观测值进行倾斜改正。

如图9所示未考虑α、β及A的系统偏差，有则应对涌浪观测值进行相应改正（参见a））。

a）根据设备安装文件，对涌浪观测值进行α、β及A的系统偏差改正，取xb=Xhv，yb=Yhv，zb=Zhv，α=δαd，β=δβd，ε=-δAe，其中，（Xhv，Yhv）为涌浪滤波器在测量船坐标系中坐标（相对于多波束探头坐标系（XcYcZc）原点），h'v为涌浪观测值（向上为正）。

由（3）、（4）和（5）式得（xa，ya，za）；

b）再进行倾斜转换，取xb=xa，yb=ya，zb=za，α=αd，β=βd，其中

（xa，ya，za）为上段计算所得，αd、βd为横摇和纵摇观测值。

由

（1）式得（xa，ya，za），其中，hv=za，就是测量船坐标系原点处的涌浪值，也即多波束探头发射中心的涌浪改正值。

7·4α、β及A系统偏差的影响估算

由前述可知，对于一体化的设备，经过综合偏差及α、β及A的系统偏差改正，若各项偏差测定精确度足够，则设备安装偏差对水深点的影响将被全面消除。

综合偏差改正全面消除了各设备安装偏差对水深点归位的影响；α、β及A的系统偏差改正消除了这些偏差对定位仪点位归算和涌浪改正值的影响。

下面对α、β

及A系统偏差的影响进行估算：

6.α、β及A系统偏差对（x'，y'）影响的估算

设αd=10º，βd=5º，δαd=1º，δβd=1º，δAe=1º，Xps=0m，Yps=-10m，Zps=-10m，得x'=-1.97m，y'=-10.97m；若取δαd=0º，δβd=0º，δAe=0º，其余不变，得x'=-1.66m，y'=-10.82m，点位偏移0.35m。

7.α、β及A系统偏差对hv的影响估算

设αd=10º，βd=5º，A=1º，δαd=1º，δβd=1º，δAe=1º，Xhv=0m，Yhv=-10m，h'v=-0.50m，得hv=0.57m；取δαd=0º，δβd=0º，δAe=0º，其余不变，得hv=0.57m，Δhv=0.20m。

由上可见，尽管α、β及A系统偏差对（x'，y'）影响不大，但当涌浪滤波器与多波束探头距离较大时，姿态传感器系统偏差对涌浪改正值的影响还是相当可观的，必要时，应进行静态测定。

8多波束系统误差来源

多波束系统需由多传感器协同进行水深测量，观测值多，误差来源也多，从参数测定的准确性，到工作环境的优劣，都会影响水深成果的质量。

内部来源有设备的观测精度、仪器的结构设计等；外部来源有风浪的影响、各种干扰信号等。

8·1观测误差

横摇观测误差、纵摇观测误差、起伏观测误差、航向观测误差、波束旅行时间误差、声速误差等，一般以设备标称精度来衡量，或通过设备率定确定。

这些误差属偶然误差，按误差传播律，最终反映在水深误差和点位误差中。

多波束水深测量是一种动态的测量，各设备的观测精度会受到现场环境条件的很大影响，尤其是有活动部件的观测设备，环境条件差时，观测精度会明显下降，例如，风浪过大能引起涌浪滤波器测量值的严重失真。

因此，系统应在各设备的限差范围内工作，才能保证质量。

8·2偏差测定误差

多波束系统的偏差测定主要包括姿态传感器、罗经偏差的静态测定和多波束探头综合偏差的动态测定。

另外，还有定位时延测定（动态测定）。

显然，静态测定的精度较高，能以设备的标称精度来衡量，但如若不进行此项改正，在一定条件下，也会造成水深误差；综合偏差是用波束进行目标测定，而每个波束都具有相当大的波束角，就好比经纬仪测角时，不是用十字丝照准目标，而是用手电筒光柱去照准目标，每次观测的精度很低，因此，必须通过大量的观测才可能达到所需精度。

8·3起伏变化和摇摆变化的影响

多波束水深测量每Ping需接收很多波束，从接收第一个回波到最后一个回波，有一个时间差，在这一接收时段内，起伏和摇摆都会发生变化，可能影响水深及其点位归算的准确性。

以30米水深，单边3倍水深宽度测量为例，边缘波束和中央波束接收的时间差约为80毫秒。

假定起伏变化为振幅0.50米，周期5秒的正弦曲线，则在80毫秒的接收时段内，起伏最大变化近5厘米；假定横摇变化为振幅10°，周期5秒的正弦曲线，则在80毫秒的接收时段内，横摇最大变化近1°，对边缘波束水深的影响高达1.66米。

摇摆变化不仅对水深产生影响，也对起伏的倾斜改正产生影响。

起伏变化和摇摆变化应采用适当的模型，进行恰当的处理。

8·4粗差处理

应能准确识别和处理观测值的粗差，否则将严重影响水下地形的真实性。

水深观测值的粗差主要通过相邻波束和相邻Ping的水深变化情况加以判别；姿态观测值的粗差则要通过观测值过程线变化来判断。

8·5计算误差

主要是近似方法使用不得当，如，声线计算时，入射角和旅行时间步长过大，可能带入较大的计算误差，可以通过改进计算方法加以解决。

8·6波束单元入射角观测误差

单波束水深测量是以波束脚印中，最快回波作为该波束的水深观测值。

如果多波束也采用同样方式测量水深，则可能使水深出现很大误差，对波束单元入射角大的波束尤甚。

因此，多波束水深测量应当能测定波束脚印内的最浅点，确定其声线和的波束单元入射角。

在对水深敏感的浅水航行水域进行测量时，要严格控制作业环境的影响，要分析水深与设备参数设置之间的关系，确保多波束系统的探测能力。

除了上述外，观测值的时间同步问题、测量船的振动等都会影响水深的质量。

9结束语

在大家尽情利用多波束带来的高效和便利的同时，要充分关注多波束存在的问题。

应当知道，多波束测深系统并没有想象的那么完美，本文前言中提的问题只是例证之一，此问题有两种可能，一是系统未测到目标，这是系统探测能力的缺陷，需要从硬件方面改进；二是系统测到了目标，但数据处理不恰当，如，水深值假信号的过滤不恰当，水深归算系统偏差改正不到位等，这需要对数采软件或后处理软件进行改进，并且，数据处理软件应具备质量检查及水深恢复功能。

通过对多波束数据处理过程等的了解，可以更深入地理解多波束，希望本文能对大家有用。

参考文献

[1]李家标等著《多波束勘测原理技术与方法》1999年海洋出版社

[2]《船用陀螺罗经》1993年