3 k0 v( Q6 P- n$ _# D% f
摘要& f7 v7 t% g& e% o+ J7 R
GNSS PPK ( post processing kinematic)技术属动态后处理技术，不受电台无线电传播距离的限制，有效作用距离可达 80km。为解决沿海长航道水深测量需在海上设置定点验潮站的难题，提出基于 Hypack2014 的 GNSS PPK三维水深测量方法。该方法无需进行水位提取与水位拟合，提高了水深测量精度; 且采用现有商业软件 TBC 和 Hypack，避免了繁琐的二次开发。介绍 GNSS PPK定位原理及 GNSS PPK 三维水深测量工作原理，阐述基于 Hypack 的 GNSS PPK三维水深测量的实施方法，在沿海长航道水深测量方面具有推广价值。
2 O- w R3 U3 ]" l( I+ i' g/ u 正文
" D5 R: M% Y! ?6 T' G7 H3 l 对于沿海长航道水深测量，当测区超出岸边水位站有效控制范围时，可在海上设立临时水位 站 ( 修建验潮站、抛压力式验潮仪)、潮位推算， 或直接进行 GPS ＲTK三维水深测量。但海上修建验潮站的难度大、成本高，因此较少采用;抛压力式验潮仪则会面临验潮仪丢失或由于仪器在海底滑动导致采集的水位数据不可用的情况，且不适用于作业频率较高的疏浚工程水深测量; 潮位推算精度受限; ＲTK受电台无线电传播距离的限制，只能解决基线长度为 20 km范围内的水位测量问题 。鉴于上述方法均存在难以克服的缺点，本文提出基于 Hypack 2014 的 GNSS PPK三维水深测量方法。
. b4 Y( K' F$ F8 x+ v; U 1.基本原理# G4 g. h$ p$ M2 K7 {
GNSS PPK 定位原理GNSS PPK技术属动态后处理技术，不受电台无线电传播距离的限制，有效作用距离可以达到80km 。定位原理: 利用进行同步观测的 1 台基准站接收机和至少 1台流动接收机对卫星的载波相位观测量; 事后在计算机中利用 GPS处理软件进行线性组合，形成虚拟的载波相位观测量值，确定接收机之间厘米级的相对位置;然后进行坐标转换得到流动站在地方坐标系中的三维坐标 。
$ {5 @2 N: }3 c; r6 v9 ]: X( x GNSS PPK 三维水深测量原理三维水深测量原理见图 1。图中 T 为水位; A 为 GNSS 天线的大地高; B为水面到海底的距离; C S 为海图基准面以下的水深; D 为动吃水; K为似大地水准面到海图基准面的距离; N 为参考椭球面到似大地水准面的距离。 7 D6 `4 R; @& _+ f. n
& @- f+ C5 ^1 M" p0 d8 C4 r5 Y9 _ 水位T = －T = N－K －A －H －D ( 1) 8 ~, N* O# n+ c4 N9 M" O
图载水深C S = B+D +T = B+N－K －A －H ( 2)
* X0 ^0 l/ ?, h* k5 b' k- K 可见，通过 GNSS PPK 三维定位技术获得流动站的平面坐标和大地高 A，在定位系统中央处理器的驱动下，同步实施采集水深数据B，并根据流动站天线至水面的高度 H 及深度基准面与大地高的关系 N-K，进行综合数据处理，即可获得测点的平面位置和图载水深 C S，实现 GNSS PPK 三维水深测量。
6 |& O) c$ K+ o' f+ x( R* W 在 GNSS PPK 三维水深测量外业数据采集过程中，采用 GNSS PPK进行导航定位，采集的平面定位坐标和高程坐标是绝对定位 ( 单点定位) ，其精度相对较低。需在内业处理过程中，用经后处理的精确三维坐标逐一替换单点定位数据，达到精确定位和水位改正，从而不需要建立水位站 进行验潮。 {. u8 o; {2 g0 q E% N' e. }3 [
2.GNSS PPK 三维水深测量
3 A: q) T8 S! e3 e# J0 L h5 Z 2.1仪器安置7 O3 l' k1 A) E, i V# A5 M) I1 S
2.1.1架设基准站
- `2 B( O S* ~: W, r* W% k 定位设备均采用美国 Trimble Ｒ8 GNSS接收机。在岸上已知控制点上架设基准站，测量类型为FastStatic，记录间隔宜为 1 s。 : V" }( `1 W5 N
2.1.2流动站设置
% U$ ~# Q. ?; }2 Y' J 流动站安置在测量船上较高处，测量类型为后处理动态测量，记录间隔应与基准站保持一致， 宜为 1 s。 ( }4 m% R3 y6 l# p
2.2测深系统安置! t0 s R9 L9 }" d# }. B, U+ p
按常规作业的要求， 安置测深系统的其他设备。
0 j3 I2 c7 v9 R 2.2.1Hypack 中的设置添加设备在 Hypack 软件中接入 GNSS PPK 时，设备驱动文件是 GPS. dll。若GNSS信号失锁，经过后处理的三维坐标精度较低，甚至无法进行后处理，不能满足三维水深测量精度要求，严重则导致水深测量的补测或 重测。为便于实时跟踪 GNSS状态， 及时发现GNSS 失锁情况，须在 Hypack软件的硬件中进行相应设置。在 Survey Devices 选项中，选择 Setup，打开GPS Setup 对话框。根据 GNSS 输出的 NEMA 数据标准， 选择对应格式 (NAME 2. 1、NEMA 3. 0 或自定义) ，并勾选所有可选项 ( GNSS 不固定，则报警并暂停记录)，如图 2 ～ 3 所示。
. [* |6 w, F1 u7 | 4 C' v! w, Z1 u$ [+ l
图 2 GPS 一般设置 " H/ J/ _+ _' t
" y. o4 r `+ @& e
Figure :
; e5 J' D3 B) _4 _3 V8 R1 C 图 3 GPS 状态码设置 ( F) m/ p" j6 y, ~
在 Offsets 选项中，Starboard 和 Forward 是相对于测量船原点 (一般以测深仪换能器位置为原点) 的位置， 分别以向右、向前为正。Vertical 是GNSS 接收机天线到水线的距离，向下为正，即图 1中的 H 值。 高程转换设置GNSS 测量的是基于 WGS84 椭球面的大地高， 而我国法定高程系统为正常高。要实现WGS84 椭球面与深度基准面 ( 沿海港口和内河感潮河段采用理论最低潮面)之间的高程转换，需要通过 2 步来实现: 第 1步是高程基准面从椭球面到似大地水准面之间的转换，即大地高向正常高的转换， 图1 中的 N 值; 第 2步是高程基准面从似大地水准面向深度基准面的转换，即正常高向海图高的 转换，图1 中的 K 值。根据不同的高程转换方法，进行相应的设置( 表 1) 。在完成测深系统各仪器设备安置、GNSS PPK初始化后，测试测深系统通讯情况。 $ z! r( ]) S7 P3 R
2.3外业数据采集
' R; ]. x, E* Q/ V7 @5 l+ z GNSS PPK 三维水深测量数据采集包括: 采用 Hypack 软件进行实时采集测深数据、定位数据 ( 单点定位)和姿态数据等。具体操作与ＲTK 三维水深测量完全相同;流动站主机和手簿GNSS PPK 数据采集。手簿开始记录，流动站主机开始自动记录;基准站主机数据采集。Figure :3.内业数据处理
8 {3 @. f w e) {% G. ~ 3.1GNSS PPK 数据处理
$ j" x8 q1 M& |; s% V; {0 {5 |# C GNSS PPK 测量作业结束后，要采用 Trimble Data Transfer软件及时将基准站和流动站的数据一一对应地进行下载、保存。GNSS PPK数据需采用美国 Trimble GPS 接收机配备的商用软件 TBC ( Trimble BusinessCenter) 进行数据处理。处理过程主要包括: ) M5 o% q1 V) e: J
新建项目，选择所需坐标系统。若没有所需坐标系统可供选择， 应在Coordinate System Manger 中新建完成。导入观测数据。Trimble 系列接收机观测数据可直接导入基准站及流动站的 DAT文件，若采用其他系列接收机观测数据，必须事先将文件格式转换为Ｒinex 格式才能使用。输入相应的天线类型、测量方法和仪器高度。基线处理。进行基线解算，并删除无法进行基线处理或精度超标的基线，保存得到固定解的基线。平差。更改坐标系统，约束基准站，进行约束平差。数据导出。在导出时选择 “轨迹 ( CSV)文件导出器” 直接导出成相应的CSV 文件，用于替换原始数据中的单点定位数据。3.2GNSS PPK 数据与水深数据联合处理数据采集软件 Hypack 存储的原始数据中定位数据是 GNSS PPK 单点定位的GGA数据，其平面精度和高程精度均不能满足水深测量要求，所以必须将经后处理的GNSS PPK 高精度平面定位数据和高程数据逐一替换外业采集的原始 GGA数据。处理过程如下:打开 Hypack 软件中单波束编辑器。打开水深测量原始数据索引文件( * . log) ，选择需要进行数据处理的原始数据文件。设置读取参数: 取消 “Ｒecalculate ＲTK Tides Using Project Geodesy”，选择 “ＲTK Tide Method”，平均水位数据周期默认 30 s，也可修改。重要的是，必须选择“Use Post-Processed GPS Correction File” 来实现经后处理的 GNSS PPK 高精度平面定位数据和高程数据逐一替换外业采集的原始 GGA数据。根据TBC 导出的GNSS PPK 数据文件( * . csv)中每个定位点数据 ( 每行数据) 的格式，正确设置各 种 参 数， 匹 配导 入数据的格式和时区( 图 4) 。
* i6 q! u! d. P' m/ x: O+ u8 e Figure :
. h4 F( c& _1 A 图 4 导入 GNSS PPK 数据设置界面 . f: A4 D- H: w, {* E# t6 p
4.精度分析& ]$ k3 l- V# S! d7 S* Z" `( |
4.1内部符合精度% W3 U: r! l3 M+ o' I7 T
在Hypack 平台下，使用GNSS PPK 和GDCOＲS (广东省连续运行卫星定位服务系统) 同时进行水深测量。GDCOＲS 采集基于 WGS84的大地高，测深数据进行人工验潮水位改正。经过统计，GNSS- PPK 和广东省COＲS 的水深测量图进行水深重合的比对，误差小于等于 0. 4 m 占 100.0%，最大0. 39 m。GNSS PPK 三维水深测量成果内符合精度满足现行《水运工程测量规范》要求，比对结果见表 2。
. P; N: c' V, {* |; e 表 2 GNSS PPK 和 GDCOＲS 水深测量图比较 8 y5 ]2 C5 M: [( P# D. ^
2 U: E; ~3 Y$ H5 E
0 ^ A t w! u$ N
. e( A% w2 y; w: _& C 不符值分布区间/m
# Q$ a, ]8 D9 w8 @5 h8 T: o 测点/个
! _8 d* X1 Q% \ 占比例/% 6 p4 W' C; x- W2 f
( －0. 4，－0. 3］
1 l) K3 [) t3 y+ ], ]8 N u, M& C) S 3
5 Y2 p* P' T3 @9 r5 J 0. 1
4 d' M; C% i% P$ y: G ( －0. 3，－0. 2］ , f6 c8 ]' g/ w3 t
2 ; a$ h/ H. ~. o) Q! u
0. 1
% A4 Q; N: [1 H ( －0. 2，－0. 1］
- M2 r/ `' w/ Q% G) U 21 9 t) N( j9 J9 s% G
0. 7 . v' e+ _/ U8 `) Y3 v
( －0. 1，0］
2 P: J5 i. v" R/ Z) A& Z 715 9 A7 S; @% `. a6 G' ]
24. 1 & N+ B4 Z, `) G
( 0，0. 1］ 6 |5 e/ j p; h- {. c7 Q
2 184 1 W6 ?7 h5 J7 p8 x5 \
73. 5
X1 E: Y% [' D% n; P ( 0. 1，0. 2 ］ ' G8 L& G- p- f8 S H# O, Q- T
39
- k5 c4 a$ B# a M 1. 3
% ]+ p4 t$ {: B ( 0. 2，0. 3］
/ r. |9 O4 ?. `+ Z* w; C/ a 3
, J1 w. {! N3 A3 m7 u6 B' a 0. 1
# N+ c- [4 o ^& G1 Y ( 0. 3，0. 4］ ' {* a2 J* X% T7 _% F4 k
4 # i9 }- L* @; m3 E. `
0. 1
, N0 g# h+ e. r- k. Y 其他
+ F% f5 O: W, |) ?2 n3 J O 0
$ H7 k+ _" Y- |* R( h" R 0. 0 + a& A/ E" e) _4 Y; v8 E
+ U8 s* T6 R7 w* o
* @, i- l" n$ g! j ) Z9 c0 r/ Q" o! ]" K
注: 不符值= GNSS PPK 水深－GDCOＲS 水深。
% U$ a% q' p' ~3 y# r( } 4.2外部符合精度
$ _: I w! N2 j; B/ e% c 在广州港南沙港池进行 GNSS PPK三维水深测量，与广东省海测大队的水深测量图进行水深重合比 对，误差小于等于 0.4 m 占 97. 8%，最大 0. 59 m。GNSS PPK 三维水深测量成果外符合精度满足现行《水运工程测量规范》要求，比对结果见表 3。
+ e( \" U& B, ~6 y5 r/ D 表 3 GNSS PPK 和广东省海测大队水深测量图比较 A. a8 |: e3 }' c" }
+ S% @! t3 F r/ S0 @
& |9 o( K! d9 n q; z/ C 8 v4 e$ y' G5 l
不符值分布区间/m
( r# V7 }8 a: v 测点/个
, F+ M! F7 M4 T1 D 占比例/%
1 B* ~/ L# L6 K6 G+ q2 O j9 ? ( －0. 4，－0. 3］ : w6 K6 Z" Z: P; T
1
# l( M, L% V6 X- ~ 0. 7
& j! y0 V) T8 a+ T ^# a) y+ |( {) F2 y ( －0. 3，－0. 2］
9 p G# \( m+ E3 q 2
. b# j" f) x5 E9 y5 u* ~ 1. 4
% t! {+ l/ l e' C* Q, ]2 Z. L ( －0. 2，－0. 1］ ! F0 ]1 }/ N" p
12
7 P) W1 q/ c4 n; K! E8 L# ] 8. 7 * ~5 U+ c, K* R8 Q- x7 r4 P! ?
( －0. 1，0］ & ^$ h; U; I& r, ~2 h
39
9 n( w! v4 N& i- G 28. 3 3 {+ f2 s2 x' L* r+ S
( 0，0. 1］ + R. `6 _, l# I+ [9 A+ ^
55 0 z, Y3 K v/ h9 d6 M! n" y
39. 9
* H, E' s/ z/ Z5 z' \ ( 0. 1，0. 2 ］
* _: D3 J6 Q% F1 `! s' [ 22 1 g! C/ l4 j: E
15. 9
7 H R# N# s* B: i4 m" _: H, R ( 0. 2，0. 3］
1 k5 ]4 [4 m2 h' [% S: } 4 3 R' \0 u% I C6 ^3 U7 W4 w
2. 9
. n6 s2 ^; c% s! A: Q. T/ m ( 0. 3，0. 4］ , L+ ~, m" y O& V/ U7 ^
1 ' t6 \, {( B' t, }, S8 X
0. 7 9 n+ \: g5 ^9 B# c% [
其他
4 C0 z% j1 G* r2 e' L 2 2 W7 B$ `9 I+ h% V& M ?$ k
1. 4
! q9 b+ s' D4 n* \" Q0 }7 R 合计
) D4 P- j9 h* {9 g! F 138 ) ], \3 @" }7 e5 @
100. 0
+ s) r7 C6 f+ B4 Q: L# u& D/ \: P v8 X
D) C; H5 Q/ c# R+ p7 w: c* Z 8 @ t& @/ W/ N4 @4 ]: W( u; M
6 {' n- U; l: }- C& i
4.3与现有水深测量方法相比的优势# l: a- Q" e7 c3 Q! b4 G/ W- j; c G
与 ＲTK 三维水深测量相比，无需 GPS ＲTK 设备所需的无线数据链，突破了 GPSＲTK 设备 10 ～ 15 km 有效作用距离的限制，能够在 80 km的范围内实施高精度三维水深测量作业。与传统的人工验潮水深测量相比，由公式( 2)中看出，无需进行水位观测工作、涌浪补偿和测 深仪换能器动吃水改正数测定，基于高精度的PPK 测高技术模式优于传统的水深测量作业模式。与采用 GPS PPK水位测量技术的水深测量作业相比，无需进行水位提取与水位曲线拟合，而是进行实时水位改正，避免了水位曲线拟合所引起的误差，进一步提高了测深精度;平面定位采用高精度的 PPK 定位数据，避免采用精度相对较低的 DGPS，可适用于大比例尺测图。 % m' H* w7 r; ~
5.结语
6 @$ D5 O- d$ F 提出的方法实现了基于 Hypack 软件的 GNSS PPK 导航定位、数据采集;内业数据处理过程中， 将 GNSS PPK处理后的高精度平面定位数据和高程数据逐一替换外业采集的原始 GGA数据，完成定位数据 ( 含高程信息) 与姿态数据、测深数据的融合，并最终实现 GNSS PPK 三维水深测量， 且避免了繁琐的二次开发。
) c- `# x( g/ F7 }0 R! F 基于Hypack 软件的GNSS PPK 三维水深测量技术，适用于沿海长航道、超出 ＲTK有效作业范围的测区。该技术已成功应用于湄洲湾航道三期工程、厦门港主航道扩建工程四期工程和港珠澳大桥工程沉管隧道工程，解决了沿海长航道水深测量的难题。
( w; t, A. r, O9 U. V! _
0 k% m3 {7 n5 F8 ?. ^# ]' ]7 W- l Y) R& }( H$ \+ d
0 f J1 _0 R, e. \! t: o
$ W1 O: z/ ^& U6 R" W6 \) Q* p1 \ | 
