 点击上方“溪流之海洋人生”即可订阅哦多波束测深波束归位计算过程涉及参数较多，同时受海洋复杂环境的影响，Ping周期换能器姿态存在动态变化，进一步加剧了波束归位计算的复杂性，而姿态改正作为波束归位中的重要环节之一，其产生的影响对于后续多波束精细化处理至关重要。 5 u8 v% F( i- S7 v6 {
对于多波束瞬时姿态影响，国内外学者进行了大量的研究。2001年，赵建虎等分析了测船瞬时姿态各参数的受动因素，并给出各角度所带来误差的计算公式；2009年，阳凡林等考虑到水下声线弯曲对波束归位的影响，提出瞬时波束向量是安装偏差和姿态影响的综合产物，并据此推导了顾及安装偏差角和姿态角影响的瞬时波束向量计算的数学模型；2015年，于家成等基于载体姿态效应分别对位于水平海底和倾斜海底的波束脚印进行改正，最终得到精度更高的水深图；2019年，金绍华等考虑姿态和声线弯曲对测量影响的耦合效应，推导了海底测深点位置和水深归算模型；陈若婷等、肖元弼等针对横摇变化导致测深数据误差问题，分别于2013年和2021年提出了基于横摇稳定和利用回波时间反演横摇角的方法，提高了波束脚印位置的计算精度；HARER等剖析了多波束测深过程中的各项误差源，推导了测深点水深及平面位置的误差估计模型，并提出了用于多波束测深不确定度计算的HGM（Hare-Godin-Mayer）模型。然而，上述研究成果均将多波束发射和接收换能器近似看作一体（本文称之为收发同体模型），仅利用波束发射或接收时刻的单一姿态进行波束归位，未顾及Ping周期内的换能器姿态变化，此外，为保障海底全覆盖，纵摇和艏摇实时补偿技术也应用于多扇区多波束系统，使得发射导向角成为波束入射向量的组成部分，在此影响下，收发换能器Ping周期内姿态变化的影响进一步加剧。2004年，加拿大新布伦瑞克大学教授HUGHES CLARKE JE利用波束传播过程中涉及的所有参数提出一种波束归位方法———虚拟同心（VCCA）模型，模型将多波束收发换能器分开考虑，可更为精确地计算水深结果。 
) @+ _$ E3 c* ~& K: X 本图来自网络
5 i4 a3 B8 M; q$ O, r 综合现有研究，正确分析Ping周期内换能器姿态所产生的影响可为多波束数据精细化处理提供一定的参考，因此，本文依据多波束测深原理，并基于VCCA模型，利用模拟数据和实测数据分析波束收发过程内姿态变化，以及在发射导向角共同作用下对波束点水深的影响。 + M) p+ g! Z- a
一、瞬时姿态影响分析  ! Q4 l- M1 d( i; `! k$ Y( n' Y
波束发射时刻的姿态主要影响发射波束能量与海底平面的相交过程，导致其偏离理想的位置；波束接收时刻的姿态主要影响接收换能器对波束脚印回波信息的接收，导致信号接收时换能器理想状态被打破，进而使得测得的波束入射角存在误差。 ; e, o5 L5 A& M) [: u. |8 T
对比收发同体模型和VCCA模型可知，前者在计算波束入射向量时利用了收发换能器整体的安装偏差、波束入射角，以及接收时刻横摇和发射时刻纵摇、航向，而后者使用了收发换能器各自安装偏差角、收发时刻两换能器瞬时姿态及波束发射与波束入射角等参数。结合波束形成和波束导向原理可知，发射时刻横摇和接收时刻纵摇对波束能量在海底投影影响微小，因此本节对发射时刻纵摇、航向和接收时刻横摇、航向，以及在发射导向角共同作用下其对波束入射向量计算产生的影响进行分析。
: e/ i. ~9 a$ i) ?* n ⒈接收时刻横摇 , n8 ` t z2 l5 C
若将发射和接收换能器看作整体，其换能器声学中心可假定位于发射换能器与接收换能器米尔斯交叉点附近。如图1所示，波束接收（BR）时刻横摇会对接收波束能量产生明显偏转，当波束入射角增大时，波束在海底投影形成的曲线曲率增加，开口变窄，在无其他角度影响时，其角度偏差大小为横摇值，波束脚印实际位置由a点变化为b点；在发射导向角的共同影响下，波束脚印实际位置变化为c点，且该误差随发射导向角和波束入射角的增大而逐渐扩大。  ' ~, { Q/ ?$ R! }! J- o
图1 接收时刻横摇影响示意图 ! e2 _& N/ v7 D. B
结合两模型，此时，在横摇和发射导向角的共同作用下，单位波束入射向量计算偏差值如下。 
% N2 F" p7 R2 p" ~ 式中，r为接收时刻横摇角；Tsteer为发射导向角；Rsteer为波束入射角。 ! x- V7 v+ e: C: D; t3 F5 M
⒉发射时刻纵摇  , F/ n2 Z5 A$ ~& g0 ]% U
图2 发射时刻纵摇影响示意图
: ^/ y0 ?0 t& q$ ` K: N 波束发射（BT）时刻纵摇使得发射换能器绕换能器坐标系Y轴在XOZ面内发生旋转，造成发射波束能量沿航迹线方向发生前后位移，进而产生大小为纵摇值的角度偏移。如图2所示，受波束发射时刻纵摇影响，波束脚印实际位置由a点变化至b点，在发射导向角的共同作用下变化为c点。
+ p q" u) \' t% Y* d7 ~& T1 [ 结合两模型，此时在纵摇和发射导向角的共同作用下，单位波束入射向量计算偏差值如下。  8 a) ^3 y5 c- P2 J2 H) P3 Z. |
式中，p为接收时刻纵摇角。 K5 `" z1 X: F0 \& ~; Y1 c5 F5 t
⒊发射时刻航向、接收时刻航向 ) ` y9 h7 B6 z4 X( e8 h
受测量环境复杂因素影响，测船实际航向方向并不会与测线方向一致，且在Ping周期内波束发射与接收存在一定差别，进而使得波束扇面绕换能器坐标系Z轴在水平面上产生扭动。如图3所示，受波束发射时刻航向影响，发射波束能量产生大小约为航向值的角度偏转，波束脚印实际位置由a点变化为b点；受接收时刻航向影响，波束脚印位置由a点变化为c点；若同时顾及发射与接收时刻航向影响，波束脚印实际位置由a点变化为d点。同时可知，当存在发射导向角的共同作用时，上述两时刻航向产生的误差会随发射导向角增大而逐渐扩大。 
0 Q) i- S+ `& A k! M4 u) ~ 图3 发射与接收时刻航向影响示意图 5 {# S; ?3 y8 l9 I; a) r. D& D9 R
结合两模型，在收发时刻航向和发射导向角的共同影响下，单位波束入射向量计算偏差值如下（为简化书写，部分表达式以A、B代替）。 
0 E) S, G1 m" C$ g3 [ 式中，h1、h2分别为发射时刻和接收时刻航向；δ为瞬时发射向量和瞬时接收向量的夹角；A、B及δ计算公式如下。 
: Z3 o: Q0 U' s 综上，波束发射与接收时刻姿态对测深结果产生不可忽视的影响，尤其是两时刻的航向变化会对测深点X、Y、Z坐标均产生影响，在发射导向的作用下姿态对测深结果的影响进一步扩大，随水深增加测深结果（尤其是边缘波束）极大可能会超出规范要求；同时也可知，对于收发同体模型，其与实际波束脚印位置间的计算误差来源于接收时刻航向、发射导向及各角度间可能存在的耦合误差。 " d! t4 l/ o' m8 Q
二、Ping周期内姿态变化影响分析
( O1 Y1 F% [: C1 L, I D/ I 结合上一节可知，波束发射与接收时刻的某一瞬时姿态会对入射向量产生影响且在存在发射导向角情况下影响进一步增大。因此，针对Ping周期内姿态变化所产生的影响，本节使用收发同体和VCCA两种模型，首先利用MATLAB进行单一姿态角变化影响的模拟实验；然后选取实测数据，通过对比数据结果分析综合姿态变化，以及在发射导向作用下对水深计算结果的影响。 ' X' }, q- \0 P8 K$ E3 l
⒈模拟实验分析
) C$ k, F& z! t% I5 |" I 本节分别分析波束发射与接收过程中横摇、纵摇和航向变化，以及在发射导向角作用下对水深计算结果的影响。设置波束入射角为10°和60°分别代表中央波束与边缘波束，姿态角度在0°～3°内变化，以发射导向角0°和3°设置两组数据，同时为简化运算过程，除实验测试参数外，安装偏差角和其余姿态角均假设为0°，数据参数如表1所示。 - b( r$ F. j; |* l
表1 模拟实验参数设置 
* V% _9 ~! w6 `& W: r 利用VCCA模型对表1数据进行处理，计算水深相对变化，结果如图4所示。其中，图4（a）和图4（b）分别表示发射导向角为0°和3°时各姿态角度变化所产生的水深相对变化。 
3 i7 m r/ { ^- b1 D; V 图4 姿态变化对水深的影响 s% L4 j ]4 a0 d, O* l
从图4可知，①相较于纵摇和航向变化，横摇变化对水深结果影响更大，且对边缘波束影响远大于中央波束，即便是微小的角度变化也会使得边缘波束易超出水深测量精度要求；②纵摇变化对中央波束和边缘波束影响近似，但随发射导向角数值而逐渐对边缘波束产生更大影响；③对于中央波束，除横摇变化外，其受纵摇变化影响较航向变化影响大，边缘波束则反之；④当存在发射导向角时，会加剧姿态变化对水深结果的影响，进一步影响测深精度；⑤由第1节可知，收发同体模型误差来源于接收时刻航向及发射导向角，在两者共同影响下，中央波束和边缘波束的水深相对变化约为0.2%和1.9%，若同时顾及横摇影响，水深相对误差扩大为1.1%和9.2%。 8 L H' m4 o/ `7 E+ N
⒉实测数据分析
% G+ y/ {! J. u 选取无发射导向作用的浅水数据、带发射导向作用的浅水数据和带发射导向作用的深水数据，三组数据信息如表2所示。三组数据地形图如图5所示。其中，图5（a）、图5（b）、图5（c）分别对应表2中1、2、3三组数据中用于计算结果比对的测线地形图，红框区域为用于精度评定的与检查测线重叠部分。 ' h8 F* T/ M/ H' N* ^; Y# D# V
表2 实测数据参数   6 d- l5 g3 Y6 V, v$ b
图5 选取测线数据光照地形图 % P8 ]9 ?9 V1 o; x, n2 y
利用收发同体模型和VCCA模型对三组数据进行处理，统计两算法水深偏差，并利用GlobalMapper软件提取部分中央波束与边缘波束的水深偏差，同时对交叉重叠部分数据进行精度评定，进一步说明姿态变化和发射导向角的影响。其中，对比1、2组数据可以分析在存在发射导向角的情况下姿态变化所产生的影响；对比2、3组数据可以分析姿态变化随水深增加所产生的影响；对比1、3组数据可以分析在发射导向角的作用下，姿态变化随水深增加而产生的影响。 9 w6 X: ` [ P, t
⑴无导向浅水数据
* A. u! `5 V! z& z( E% E/ D6 | 利用两模型对第1组数据进行处理，计算两模型水深差值并提取部分中央波束和边缘波束结果，并对结果进行统计，如表3和图6所示（统计结果为差值的绝对值）。其中，图6（a）为两算法水深计算差值，图6（b）为提取的部分中央波束与边缘波束水深差值。 ) j; I# m" x$ _% Z# B
表3 无导向作用浅水区两种模型计算对比分析结果   ! L/ j0 S% W" T1 J! e# H) P
图6 收发同体模型与VCCA模型水深计算差值
& J3 ^! C4 r: F6 G 结合表3和图6可知，对于无导向作用的浅水数据，两算法水深误差在0～0.54m，均值约为0.01m，约有98%误差分布在0.2m内；中央波束误差与边缘波束误差均值分别约为0.005m和0.007m，误差相对较小。对测线重叠部分数据进行精度评定，并参照水运工程测量规范中相关精度要求进行统计，结果如表4所示，此时两模型精度结果近似，姿态变化产生的影响并不明显。 [. E$ K+ U+ X J# g. ]) B4 F
表4 无导向作用浅水区两种模型计算结果精度评定  8 e% J P( X6 w& `5 f
⑵带导向作用浅水数据 & j, t- d; A) {$ F
利用两模型对第2组数据进行处理，计算两模型水深差值并提取部分中央波束和边缘波束结果，并对结果进行统计，如表5和图7所示（统计结果为差值的绝对值）。参照《水运工程测量规范》（JTS131—2012）中相关精度要求进行精度评定，结果如表6所示。 V2 G4 g& ~+ q3 y8 i- X- u. C: Z
表5 带导向作用浅水区两种模型计算对比分析结果  fill=%23FFFFFF%3E%3Crect x=249 y=126 width=1 height=1%3E%3C/rect%3E%3C/g%3E%3C/g%3E%3C/svg%3E) - h; S6 h- p0 w$ J2 A
图7 收发同体模型与VCCA模型水深计算差值
$ l" ~! Y6 j* E' Q/ D 表6 有导向作用浅水区两种模型计算结果精度评定 / C! ?" F' _" W$ b r! n
对水深差值进行统计并结合图7可知，对于带导向浅水数据，两算法水深误差在0～18.43m，均值为0.47m，约有89.3%数据误差分布在1.0m内，且存在部分较大误差的数据（约1.4%数据误差大于2.0m）；选取的中央波束与边缘波束误差均值分别约为0.38m和0.57m，水深误差在1.0m内的数据占比分别约为94.8%和83.6%。同时由表6可知，收发同体模型深度不符值中误差及均值较VCCA模型大，且存在部分超出2%水深精度的点。对比上节数据可知，第2组数据姿态变化较第1组较小，但在发射导向角的综合作用下，姿态变化所产生的影响凸显。 7 ^9 `' W& |0 d, N3 P% Q3 q
⑶带导向作用深水数据
/ ?9 w6 S- y# I, x5 N 同样利用两模型对第3组数据进行处理，计算两模型水深差值并提取部分中央波束和边缘波束结果，统计结果如表7和图8所示（统计结果为差值的绝对值）。对重叠区域进行精度评定，结果如表8所示。 * r0 Z" Y- R/ l+ L# {9 p( x8 z
表7 带导向作用深水区两种模型计算对比分析结果 8 o$ ?6 C6 R+ [4 Y- z0 [0 v
图8 收发同体模型与VCCA模型水深计算差值
+ T7 p# B3 @7 w( L; z( z6 o 表8 带导向作用深水区两种模型计算结果精度评定 2 B6 e- Q$ o# @& I6 l! O
对水深差值进行统计并结合图8可知，对于带导向的深水数据，两算法水深误差在0～212.09m，均值为9.23m，约有77.4%误差分布在10.0m内，约13.0%数据误差大于20.0m，且存在有极大误差的数据；中央波束误差与边缘波束误差均值分别约为1.39m和6.47m，水深误差在10.0m内的数据占比分别约为99.2%和83.3%。由表8可知，收发同体模型深度不符值中误差及均值远大于VCCA模型，大于2%水深的超限点也进一步增多。对比本文上两节数据可知，随水深增加，姿态变化影响进一步明显并扩大，且在发射导向角的共同作用下，最终影响数据处理结果。 ! [$ r4 P* L+ c
本文实测数据选用R2Sonic2024和EM122多波束系统采集数据，其标称测量结果符合IHOS-44特级精度要求，对应三组数据分别约为0.33m、1.5m、15m。对比前小节可知，除无导向浅水区数据符合要求外，在发射导向和姿态的综合影响下，当水深大于200m，存在部分数据超出标称精度要求的情况。 ; _2 V7 H g$ n4 ]$ F
三、结论
$ r0 b/ N- N; g 本文从理论、模拟实验和实测数据三方面出发，对多波束Ping周期内的姿态变化及其在发射导向角的共同作用下对测深点产生的影响进行了分析，结论如下。
# P6 K# Q( U G& O6 _2 O ⑴姿态角度通过改变发射与接收时刻的波束能量来产生影响，其中，发射时刻纵摇、接收时刻横摇及两时刻的航向会对最终水深结果产生较大影响。
( M& [! t- t# w3 O) {2 G ⑵在单独考虑某一姿态角度变化下，横摇变化对水深结果的影响最大，且边缘波束受影响远大于中央波束，在横摇和发射导向的共同作用下，其水深相对误差可达到9.2%。
4 ~5 \# U5 R) e7 F ⑶浅水环境下，对于集成安装的常规多波束测深仪（如R2Sonic2024、SeaBat7125等）采集的数据，由于其Ping率较高，姿态变化产生的影响可近似忽略，但当使用船体固定安装的深水或全海深型多波束进行浅水区作业、设备具有发射导向作用或水深相对较深（大于200m）时应考虑其影响以得到更精细化的数据。
; I- j Z. V5 V; M- M1 V ⑷中深水及深水环境下，测量结果受姿态变化影响明显且迅速增大，此外，随着水深的增加和海况条件的影响，波束往返时间可达数秒甚至十数秒（特别对于深水或全海深型多波束系统，在不理想海况条件下千米水深Ping周期也可能达十几秒），姿态变化明显的同时若存在发射导向角，则影响会更加明显，极易导致水深结果超出规范要求。 % r! `( Z, Y* _" @- W0 N
⑸在数据采集过程中，应根据实际情况调整船速以及开角，当姿态变化明显时可适当降低船速，减小开角，以减少测点间距，同时保证航向和船艏向差别控制在小范围之内。
4 [2 @9 n! Q' N) I$ Y6 h2 s' k 1
: I8 ?. {, b+ ^1 J* W% L8 L END # O) v' m0 J+ U( p) F. }
1 7 Z' |7 @2 A ]' v( {; u
【作者简介】文/吴冬强 王晓明 卜宪海 于宗泽，分别来自自然资源部海底矿产资源重点实验室中国地质调查局广州海洋地质调查局、中交星宇科技有限公司、山东科技大学测绘与空间信息学院。第一作者吴冬强，男，1994出生，硕士，助理工程师，主要从事多波束数据采集与处理；通讯作者卜宪海，男，1990出生，博士，讲师，主要从事海底地形地貌测量研究。本文为基金项目，国家自然科学基金资助项目（42072324）、国家自然科学基金青年科学基金项目（42204049）、山东省自然科学基金青年项目（ZR2022QD008）。文章来自《海洋技术学报》（2023年第5期），参考文章略，版权归出版单位与作者所有，用于学习与交流，转载也请备注由“溪流之海洋人生”微信公众平台编辑与整理。 , w9 `5 L& I6 v2 I1 {# {4 H# H0 q
相关阅读推荐 海洋论坛▏浅水多波束测深声呐关键技术剖析海洋论坛▏几种探测方法在海洋地形测绘中的应用海洋论坛▏利用多波束获取水下桩体点云的方法研究经验交流▏多波束测深系统在海上风电地形测量中的应用与探讨海洋论坛▏海洋测绘中的测深技术分析海洋论坛▏多波束水柱图像联合降噪方法研究 3 v3 `4 Z4 z9 l! ^! ` ~9 R7 `
公众号
9 B8 F, j W# ^ 溪流之海洋人生
& R7 {! d8 ]) u: j3 r 微信号▏xiliu92899 $ ]9 X% j6 C9 H6 Y7 s2 L T
用专业精神创造价值
3 I4 q' s7 d( V7 [/ O1 r 用人文关怀引发共鸣 您的关注就是我们前行的动力 $ S: s$ `2 e1 {& L* h
投稿邮箱▏191291624@qq.com
' H) j0 W! U2 ]6 [. q7 F 【温馨提示】：由于微信修改了推送规则，没有经常点“在看”的，会慢慢收不到推送。如果你还想看到“溪流之海洋人生”的相关信息，请将“溪流之海洋人生”加为星标。 3 J9 S% v% N, J) W9 y
$ v- u( l3 l- _0 b; o M1 J6 T9 S
4 E; a4 U1 O) D1 Q& I ]; g
$ ^( g2 j3 S. o* v
- s* G7 A: c3 V% { | 
