 点击上方“溪流之海洋人生”即可订阅哦水道地形测量作为河道资料的主要获取方式，其模式已基本定型于利用GNSS测定点的平面位置，利用测深仪测定点的水深，附之以瞬时潮位或水位资料获得高程，以此确定点的三维坐标故水位资料是生产过程控制的关键环节。有关河道观测的水位的获取及推算论述很多，如黄佑生等提出了由上游往下游推求的水面线试算迭代模型；周儒夫等提出了水道水位节点布控和水下测点高程推算方法等；也有采用GNSS测量值作为水位的方法等。 
8 u% Z$ B" D$ Q3 B l 但针对具体河段如何选取合理的水位获取方式和推算模型方面的文章较少，为了验证这些不同水位推算方法对水位变化复杂条件下的工程河段水道地形测量精度的影响，笔者选择了宜昌河段的某一区域进行了试验研究，采用不同的软硬件、不同的水位获取方法及推算模型进行对比试验，为水位变化复杂条件下工程河段水道地形观测提供了理论及实践依据，具有较强的针对性、实用性。  . w5 \4 ^, K% A0 C9 }! J6 F
一、水位推算模型与原理  , M; K' b- W/ x" ^
⒈线性内插推算模型
" H( `; t/ s: ~& _6 F/ n8 ] 对于恒定均匀流量的天然河道，其水面线计算原理基于一维能量方程，然后逐个断面采用直接步进法推求，计算公式如下：
( e/ F$ R3 v7 G Z2＝Z1＋hf＋hj＋(a1u12/2g)－(a2u22/2g) ⑴
3 @. n3 D3 {0 a8 T- f4 A n 式中：Z1、Z2分别为河道下游、上游断面水位；a1、a2分别为河道下游、上游断面流量系数；u1、u2分别为河道下游、上游断面流速；g为重力加速度；hf、hj分别为上下游河道断面之间的沿程水头损失和局部水头损失。 6 ]; M7 B! c$ ?" @+ T) q
实际项目生产中，特别是水位不变或均匀线性变化的情况，往往按照上下游（测量区域起止点）断面水位落差和断面布设情况，采用平均落差逐个进行推算。 Z6 n9 F9 W' e' e, n3 U
⒉时差法水位推算
7 c- ^# J0 N% }" _1 b) m% r 大多数测量作业区域的河道水位变化为非线性、非均匀变化，故水位推算也不能简单采用线性内插的方法，这时可以采用时差法水位推算模型。主要原理如下：首先在水道观测一岸的开、收工位置布设水位节点（存在横比降测区可以在河道另一岸合理位置加设水位观测水尺）；从水下开始前到结束后按照规定的频次或水位变化的情况进行水位观测。然后根据测量作业布置的中心线法或三点法（也叫区域法）进行水位推算，称之为时差法水位推算。根据不同的条件，可采用二步内插法、平面内插法、距离加权法等进行水位推算，水道观测规范（SL257-2017）给出了具体的推算公式与要求，这里不再赘述。 9 n8 a( b2 L; u* k
⒊GNSS潮位法推算
6 t# r/ E" b# y. x2 G 以上两种方法中，潮位采用了人工或自计水位进行推算。在水道测量时，有时水位通过GNSS测量获得，称为GNSS三维水道观测技术，即利用GNSS观测值作为潮位的模式进行水道地形测量。 + \+ M( i% T0 x# E8 B
若设H为WGS84参考椭球下大地高，L为GNSS天线高，d为测深系统探头静吃水，T为基于当地基准的当地潮位改正，S为测深系统的测量水深值，ζ为高程异常值，h为基于当地的河底高程。通过相互关系可知： 6 q6 {) l# M( d) D3 Y: @& Z
h＝S＋D－T ⑵
2 `8 x+ Q t( U6 _" Z5 M T＝H－ζ－L ⑶
, \" a) g2 W/ q3 H" l! ?! h h＝s＋d＋L﹣(H﹣ζ) ⑷
, u# h' ^/ j, K5 Y' G 若当地基准面为 1985国家高程基准面时，ζ为高程异常，此时H﹣ζ＝H85高程。则式⑷可转换得到：
- B# Z1 r+ D9 s+ F, J2 R1 x: t h＝s＋d＋L﹣H85 ⑸
}" v" l: r5 e+ }4 W( L 根据项目外业生产可知：S由测深仪实时测得水深，d＋L为固定值（钢卷尺丈量），若 GNSS接收机实时采集到1985国家高程基准面下的正常高，便可实时测得水下1985国家高程基准的河底高程。
: x5 g7 s1 w! d' Q 二、试验布置及数据处理
5 m8 }* V$ ?5 C' v" A ⒈试验实施 0 n" [7 W' \; f- w# B
本试验选取河段为三峡工程坝下游宜昌河段，一是本河段受水利工程影响水位变化大，且变化过程经常呈非线性变化；二是河道基本比较顺直、河床稳定，底质为山体礁岩构成，测深仪施测过程中其模拟回声效果良好。项目实施区域示意图如1所示（示意图进行了缩放与旋转调整）。  : N2 J% x3 H X/ G+ P% F5 s2 z
图1 项目实施河段、控制分布及断面布设意图 7 l( V4 _% ^2 L* q# a
⑴控制网布设及转换参数求取 6 E g. q4 g5 O1 @3 { S/ u
在长约20km的带状测区内，上中下合理选取3个国家C级控制点，通过测量及解算，求取国家坐标转换七参数，平面及高程转换残差满足CH/T2009-2010标准要求，并对转换参数进行实地检查与校核工作。 $ O; z K9 }* H! H7 G4 n
⑵水位控制基点测量
8 I5 r6 N& M+ m' n- u# ?5 g 在转换参数控制区域内，选择长约5km左右范围的测量作业区域起止埋设2个水位基点，采用四等水准观测精度联测国家二等水准成果，用作水位施测起算点（河段基本顺直，根据长期资料积累可知，本段不存在横比降情况）。
! g9 N; `, h% d) Y0 \2 M% O( S' ? ⑶区域水道地形测量 # ]8 B f0 N" M
①GNSS基准站及流动站设置：在水下测量区域中间，架设GNSS基准站，设备系统配置中输入当地坐标系统及求取的七参数，以准确得到当地实时平面坐标与当地基准高程值，流动站接收到基站信息，进行实时改正，并选择已知点进行精度检校。 % h+ [' m5 A/ I0 r0 Y
②水位控制测量：在水下测量区域的首尾进行水位控制测量，精度达到五等水准精度，同时设立人工水尺和潮位自计仪；在水下测量期间，采用人工和自计两种方式记录水位变化值。
% `8 Y: x! H- J8 {. n' O9 b ③水道地形测量：水下作业前先将测深仪各部件连接、设置、改正及检校。将GNSS天线与换能器采用一体杆连接，确保其在同一垂直平面上，并量取GNSS天线相位中心（或天线底部）到换能器底部之间的距离，同时精确量取回声仪探头吃水深度。 . ]! M- {5 M6 X+ v1 Z: l$ |/ b3 z
本次试验采用3m×10m铁船作为测量船（此船长期在本区域实施测量工作，稳定性较好），在头部约三分之一处左边安装一套中海达测深系统，单独数据采集与处理；右边安装海鹰Hy1601测深系统，海鹰Hy1601作为水深采集单元，将水下和平面数据分别输入Hypack软件和亿点通海洋测量软件，实现全野外一体化、实时同步水道地形测量与数据采集。 , R# C( j; B6 P* u5 J8 W. w
水道地形测量采用横断面法，断面方向应大致与水流方向垂直，断面间距为180m，测点间距取5m。作业过程中应用Hypack导航测量软件现场严格控制船的走向，随时关注行船的平稳性、GNSS跟踪的卫星数量、卫星图形精度、基站信号和模式及观测精度等情况。三种软件自动记录流动站测得的水面点的平面位置、高程，同步瞬时记录测深仪测得水深值。 : Y0 r) {2 A8 E) h0 C7 ?
④为了验证单波束GNSS验潮与船的姿态关系，在试验过程中，将一体化罗经与运动传感器Octans的改正参数输入了亿点通海洋测量软件。
5 r B* f7 I3 l. s6 k% u2 C# ]0 P ⒉数据处理 O3 P( X- h" ?+ J; a
本次试验数据处理总共得到以下几种数据形式： % d. S! @' r9 Y4 h7 X
⑴采用在水道地形测量开始及结束的位置观测的水位值，按照断面个数及水位落差进行逐个断面水位推算，即线性内插法。
! K+ s S% j1 w7 K( i: { ⑵采用Hypack软件按照时间和距离的中心线法（时差法）内插推算模式推算每个断面的水位。 - }: l0 K9 P: n( u: B9 O9 U
⑶分别采用中海达测深系统自带软件和亿点通海洋测量软件进行处理，得出验潮模式下（时差法）每个断面和每个定标点的水位，以及GNSS验潮模式下的每个断面和每个定标点水位。
9 S1 s2 [& L% V; W 三、水位获取与改正分析
& o( ]- r( f+ R F ⒈试验辅助观测分析
" Q6 Z8 }: u, t ⑴人工水位观测
4 {! `+ q, i, W/ X6 @1 O; O. { 测量作业期间通过人工和自计的方式对现场水位进行了观测，其变化情况如图2所示。其测量作业期间水位是先上涨、后回落、再上涨的一个变化过程且开始阶段涨幅较大，整个作业期间水位变化达0.5m左右，是典型的工程调节模式下的非线性变化过程。 
4 T; u5 G4 U; c4 _; M! t& D# \+ x 图2 测量作业期间上下游水位变化情况图 1 F, J* t! {0 ~! e0 |
⑵声速剖面观测 7 a9 G& n; A+ y0 V# {0 |
在测量期间，试验组于上午、中午采用声速剖面仪对测区内进行了两次垂直声速测量，为了便于水深声速改正，测量结果见图3所示。通过观测得知，无论在时间和空间上，其声速变化较小，其影响可以忽略。  : ^* X1 P0 v: g8 |/ j- {
图3 声速剖面测量示意图 g. A6 P' A) T) ~) e0 I* Z
⒉试验数据对比分析 * b2 d* |" b1 d( f' q7 g0 @
⑴不同测深系统测深对比 6 w& b# t& Z3 c8 [: J
为了验证不同测深系统的水深观测值的是否存在差异性，试验组挑选了水下变化情况复杂的SY16号断面进行对比分析，如图4所示。  8 @& i9 p- u% F" f5 a# m1 S
图4 不同设备同步水深观测对比图
1 w' W) I. w2 f0 q% Y0 Y" ` 通过查看现场的模拟回声及后处理回放数据对比分析可知，国内两套测深设备系统观测数据基本无差异，满足水道地形测量生产的需要。 % y* O n5 Y$ {5 ^
⑵不同水位获取与改正方式对比分析
- g9 d6 c& U4 N" a 本次试验采用基于Hypack时差法（时间和距离的中心线法）推算的水位作为基准值，将传统断面个数线性内插法、中海达测深系统和亿点通海洋测量软件处理人工（自计）实测验潮数据(时差法)、GNSS验潮数据分别作对比分析，其结果见图5及表1所示。  * @6 L( g7 o, H/ p$ ] S3 n
图5 各种水位推算方法与对比分析图
2 W; C( Y( v7 |& K" S7 Z 表1 不同水位值与基准值互差分布
$ Z1 K% p* V u/ \- @7 y 设备/软硬件
" x: ^6 o: z% ~% ]% W2 X a 中海达测深系统 9 l% |# \4 t9 w, @
亿点通海洋测量
1 {3 N. G' V {" _: i; F5 ?& E 软件 " ~9 W8 ?' E9 _; S
统计 a0 O' r3 l( U! _- a( n8 c4 q& Q
时差法 * F" O& N0 I0 c$ I+ x! s( U
GNSS潮位 h/ A, l, I" ~$ Y
时差法
: V" v, L V# n/ p GNSS潮位
# [" q% n, s, g* p# `( y/ S ＜0.05 % D: U7 m1 @" A, l; [
26个
2 k" W* l; c. V# z 4个
n7 M& n5 R+ G2 M3 b: A$ n! c 26个 ' \6 [5 Y& V: {8 ], }
21个
+ W; e" t% z5 }# y3 ^2 Y 0.05～0.1
' ?" N5 ^* N; {) u$ }3 G 0个
J" A) A5 `5 c. o0 I! | 19个 5 y) B6 p$ f- l! N% N; I& {
0个 - R2 v* J; C/ [* j
5个
3 e4 M8 _; M/ e5 T& z- e ＞0.01 , M& r. V: g$ t2 Z, o8 S1 `
0个
; M' b6 R: ]7 c2 [ 3个
: [. G4 {0 [5 n$ N0 U$ }' i 0个
3 Y) T5 U; P* c 0个
- Y. }0 w, A" |+ v 最大值
% J2 V5 N# b) a' g ]; `, R4 o8 {( c 0.049m
$ d9 |+ M! } Y& n 0.130m & g: c5 ?2 r+ e5 H4 {6 E! u' Z; Y
0.02m
" u, ]2 I: v7 ]5 K( m 0.068m
. \6 {1 C" n s2 {& I. \2 g 为了验证姿态改正对GNSS潮位值的影响情况，实验组对接入姿态仪的亿点通海洋测量软件，在数据后处理时分别进行姿态改正和不改正两种情况数据分析，取其中一个断面进行对比，其结果见图6所示。 
& H: A6 M) o& \# A/ {' `: x+ J 图6 亿点通海洋测量软件进行姿态改正与不改正之间的比较
- P" o/ a' M+ } 从图5、图6和表1可知： ) \6 Y# |4 T3 Q% ~0 F
①工程调节河段，其水位呈现非线性变化的情况下，不宜采用在水道地形测量开始及结束的位置观测的水位，按照断面个数及水位开始和结束时的落差进行逐个断面水位推算的这种线性内插方法。
0 T6 j7 k% Z" S ②不同软件处理系统对数据进行处理，在时差法上对比可知，其互差均小于5cm，故采用人工（自计）验潮值通过时差法进行水位改正，精度最高且符合实际情况。
; x* B! W' [# ? ③中海达测深系统（未接入姿态改正）GNSS验潮值经过粗差剔除及多次滤波处理后，按断面取平均与基准值对比，互差最大值为0.13m；亿点通海洋测量软件（接入姿态仪）GNSS验潮值进行了姿态改正、粗差剔除及多次滤波处理后，按断面取平均值与基准值对比，互差最大值为0.068m，且小于5cm的达80.7%；说明采用GNSS潮位接入姿态改正是很有必要的。
* j1 P) t6 M1 [9 L7 o ④根据《水运工程测量规范》，在水深≤20m时，要求水深测量精度σh≤0.2m，按照水位改正精度为σh的1/3，即≤0.07m的要求，表明接入姿态仪后，采用合理的处理方法，GNSS三维水深测量可以达到此精度要求。
& N- @0 H( t6 z2 P: A+ ] 四、结论与建议 3 L4 [, H0 P) r9 z6 N3 Q0 d/ E
通过试验及分析表明，不同的水位获取与推算方式，其值会出现较大差别；如果采用的方法不合理，甚至会出现水位值错误，从而导致水道地形测量精度不高的情况发生。所以在水利工程调节影响、水位变化复杂的河段进行水道测量时，水位获取与推算方面建议如下：
' _& U' v1 L. C1 V; d# l' F- S! o ⑴结合水道地形测量的项目用途、成图比例尺、测区情况等因素合理选取水位的获取与推算方式，特别是基于非线性变化的河段，其人工（自计）验潮模式下的水位推算方式（时差法，中心线法或区域法）最为可靠，但要注意水位零点的误差传播。 7 i. K) Z3 u. B: }6 I6 Y9 `
⑵通过对比分析，采用断面平均的GNSS验潮方式（有姿态改正），满足《水运工程测量规范》相关要求，可以根据项目用途及要求合理选取GNSS验潮方式。
$ L; [& O! x; \8 Q) p$ L. ?/ F ⑶试验表明，姿态改正能有效提高GNSS验潮精度。但是在数据处理过程中，要应用软件剔除粗差，选择适宜的GNSS 滤波方式及次数。如果是长江这种带状河段，可以优先选取断面平均法作为最终的断面水位。 6 M4 f [! |, h/ V) v! P
总之，受水利工程调节影响河段水位获取与改正方法对水道测量精度的影响较大，在类似河段进行水下地形测量时，要充分分析河段特性，选取合理的水位获取及推算方法。 5 q, r9 j0 G/ J# @
1 " T6 M7 V: M/ |# f0 y. g
END
& k- w) Y* b, N5 |7 s 1
) _0 U+ O7 E' Q7 ]* a 【作者简介】文/解祥成 郭文周，均来自长江水利委员会水文局荆江水文水资源勘测局；第一作者解祥成，1981年出生，男，任职于水文荆江局河道勘测中心，大学本科、工程硕士，长期从事河道勘测工作；本文为科技项目，国家重点研发计划专题《堤防工程风险识别和监测预警技术研究（2017YFC1502604）》的分项研究内容；文章来自《水资源研究》（2019年第6期），参考文献略，用于学习与交流，版权归作者及出版社共同拥有，转载也请备注由“溪流之海洋人生”微信公众号编辑整理。    8 [' |8 L. d. v8 P% q4 W; ]
相关阅读推荐 论文专区▏海道测量规范中水位控制部分修订的要点海测技术▏顾及地形变化的多测线粗差剔除方法在单波束水深测量数据采集和处理一体化中的应用论文专区▏海底地形测量成果的质量检核评估（四）：自容式压力验潮仪验潮零点漂移处理实践海测技术▏GNSS支持下的无验潮测深模式优化海测技术▏内陆水域水下地形测量中声速改正方法探讨技术交流▏海洋测绘中的余水位应用     * z; [/ ^' }, p+ w7 V
公众号 0 \# ^" I/ Q6 W; ~8 a! @0 r
溪流之海洋人生
3 L$ T' h W3 {' T$ c8 x 微信号▏xiliu92899 
5 i' _: A/ W3 l9 l. W( j 用专业精神创造价值 * @1 Q" [% t( Q! y; {
用人文关怀引发共鸣 您的关注就是我们前行的动力
+ f% A$ y( k0 q0 W. j- K 投稿邮箱▏12163440@qq.com  
- q( f6 p7 _1 {* l' \/ {7 K1 X3 y5 k& G( L) E: [* \2 S2 J
' L8 W; f. }8 K$ w, v/ H7 Z
5 L1 C: R# F0 N
5 W1 i& f7 U3 |$ E | 
