0 }; t, O; ]) H0 R7 e7 J
! V t- t$ E$ d# i1 W0 v+ C0 u0 b( e# J ⒈ 试验方法 $ {' Z8 s) a3 q" ^% b
结合部队年度生产作业，按照GB12327-1998《海道测量规范》的要求进行测线布设和海上作业，将采集记录的测量数据分别按GPS测高模式和传统人工验潮模式两种方式处理，对两种模式所获得的数值成果、图件成果进行定性和定量分析，同时将GPS测高模式所获得的数据反推出潮位数据，与验潮站人工观测数据进行比对，从而检核GPS测高模式水深测量技术的正确性，以及所获取测量成果的可靠性、精度及其误差分布特性。 , ?3 x- H7 ] b' G
本文海上试验所完成的水深测线轨迹如图3所示。 7 c1 d5 X4 `$ |4 H- W, n6 n3 ~% m
% d' }& V9 G, \$ {2 d2 g$ a ⒉ GPS测高反推潮位与人工潮位比较分析
7 T2 ~2 n* Y+ d E. Q 根据每一测点处GPS天线的大地高可由下式反推出测点潮位，该推算潮位中包含了测深仪的动态吃水变化。 : |8 I) J8 b4 ]* f+ V
hi=H-Hop-(ha-测深仪静态吃水) : t, f- I% y G' {
⑹式中,各符号的含义同式⑶～⑸。
4 b1 L$ i+ |2 O ?, [- m: t7 g 图4为海上试验期间8月31日GPS测高推算出的潮位与沙埕、海尾两验潮站人工验潮潮位曲线图。 # t7 G/ @% |4 S4 c5 J6 Y
( U0 P3 ?4 }4 x6 C
同时，我们对8月30日、31日、9月2日这三天由两种方法计算的潮位数据进行了数值比较分析，其结果列于下表。
; `' m) j/ |- B9 Q+ ^- Y! g 表1 GPS测高反推的潮位与人工验潮值的比较结果
" i0 p$ `$ J/ P4 b/ \/ j
' v& e* h3 T% [2 s 6 K. U7 {7 Q9 h4 g- {. h7 `
日期/时段
5 B2 j! o3 [$ {7 m+ X( W/ U+ k, ^" { 最大互差(m) : h! r( T+ P9 b8 X
最小互差(m)
; `4 M% _; U- ?' e, _, J 平均差(m)
# q+ n8 v0 G0 O3 K% s 均方差(m)
' y/ V4 s$ N* X' g+ U 标准差(m)
! p+ M8 U" e5 X/ L 8月30日
; Q3 Q6 Q) R# S/ @8 S 0.506
( A, c3 {! ~" N( ]' ~2 B D -0.017
$ q, G. o2 P" X: i3 |4 |3 J 0.182
4 |" u7 A- B l 0.194
: m+ `" ?% F% M# L' r) G 0.068
- x1 M3 T* C* l 8月31日
% r( r* R5 ?5 U 0.468
1 V% _; F" I: u2 g/ u! v& T9 g' X -0.094
+ t. c ]7 y5 \* Q8 a7 d% K 0.101
/ _9 ]6 D5 I4 A9 y! b: Z 0.172 4 }. C, X% C1 h, U1 p2 o+ Z) {) D" }$ F9 }0 _
0.139
/ E b/ u' ~$ A' \* g& p* { 9月02日
; s, t: ^" G; F1 a 0.434
$ d4 c3 \: v: X. ]: M# I. w+ a -0.052
7 S$ S7 q$ a4 W# } 0.186
3 D& q. s4 u" t Q$ i 0.197
4 a# U+ B; B: m% L 0.065
/ c) O! q, {" @* v- _7 S+ l ) P% g$ v$ m2 a( h+ |7 z- _
) ~! v; q7 E% {8 A" B; s2 _
从图4及表1可知，由GPS测高数据反推出的潮位曲线与两验潮站的人工验潮曲线总体上一致，但由于测区位于两验潮站中间，在理论上测高曲线应该位于两验潮站潮位曲线中间，而实际上与理论曲线相比存在10～18cm的系统差，根据有关文献的研究结果，这正是船只的动态吃水及涌浪等因素可能造成的影响量值，这也从另一个侧面说明，基于高精度GPS测高的水深测量作业模式能够有效地消除船只动态吃水等因素对水深测量成果的影响。 + K! W5 M# U# g: Q' T: k
⒊ GPS测高作业模式与传统作业模式内符合精度及精度分布规律比较 7 w. m% z& {9 w1 Z" z4 A
本次海上试验，上半测区航道水深变化比较剧烈，下半测区水深变化比较平缓，因此我们将试验数据分成两个区块进行数值分析比较，数值分析的主要方法是：通过计算主检测线交叉点水深不符值的统计规律来进行。 % f/ Z/ z$ A4 n
按GPS测高模式处理得到的上半部测区主检测线交叉点不符值分布规律见表2。 % p7 ^0 q" q @) |5 `/ C' y
表2 GPS测高模式处理的上半部测区主检测线不符值分布规律统计表
# o, i! _5 G- a! A. m+ P ! c% f! R( v8 s# u
/ u6 ^" U9 J/ R- J" U9 ~ 不符值范围
/ {1 V1 h+ s6 q$ T) O 点数
0 O" d% x% D5 f& z8 e* b 所占比例
S, t; Y* b, E! x" Y) V7 ? 0.0≤|dh|≤0.1 % p5 P6 H. Q* F2 Q1 z
119 ) {9 z6 c9 H; M0 }- j9 U
35.20%
) ^, U+ b" ^8 Q* M' M1 l# Q 0.1<|dh|≤0.2
" X" M- h9 d, m# L9 Z$ F 63
) k4 ?* ` z+ Z6 H3 G: P# I 18.64% v3 K, X2 O, t, P
0.2<|dh|≤0.3
6 y: E5 t) v' @" b! t 51 # q) C% v' Y! f4 m* i- J9 B2 j2 `
15.09% % N. j# U+ Z/ f& U: Q; \# o
0.3<|dh|≤0.4
9 [( R L) x8 e' |4 N: L- Q 38 ' o+ P$ N; }, j! }( Q" ]
11.24%
. k1 W9 ?: e' Q, l5 h* b; S- } 0.4<|dh|≤0.5
o! D! z( ^$ X 22
* t, d2 |8 D; z# k1 U 6.51% 6 n7 }4 w' E. f: z: D& q/ Q# [
0.5<|dh|≤1.0
; ]. Y! B0 M! W: }+ v) f4 a 29 - \& h, t- A: J. o! i+ u$ D
8.58%
( O# X7 a6 E2 e! L' ]9 @. _, m% n/ R 1.0<|dh|≤2.0
4 T- s3 J6 X. F( M 15
5 f1 D8 S! h' {& R5 |% ?( j- U 4.44% % f" k4 g2 l. A3 C) e+ n! n7 C
2.0<|dh|
6 C6 ]! O1 N3 J4 I$ a 1 + I; ^* y, U/ K7 `; u2 V' N
0.30%
! n1 Q. E4 U# o( j+ |, s$ V 合计
$ I# S9 J: X( ^* z6 X 338
8 b' Y0 T9 _; p& i 100%
3 F+ r% }+ A5 l& t l/ O7 V& P 7 W. b5 j- B- T, U/ X0 O8 F
+ e. l/ B% d5 o: E GB12327-1998《海道测量规范》规定，对主检测线交叉点不符值进行系统误差及粗差检验，其主、检不符值限差为：水深0～20m时为0.5m；水深20～30m时为0.6m；水深30～50m时为0.7m；水深50～100m时为1.5m；水深大于100m时为水深的3%。同时还规定超限的点数不得超过参加比对总点数的15%。从上面的主检测深线深度比较表可知，按GPS测高模式获得的上半部测区水深成果完全满足GB123217-1998《海道测量规范》的质量要求。 ! |. F4 C( q/ I' [ v
采用沙埕、海尾两验潮站潮位数据，按传统作业模式处理得到的上半部测区主检测线交叉点水深不符值分布规律见表3。
$ e5 N6 X2 V, Q/ I8 k 表3 传统作业模式处理的上半部测区主检测线不符值分布规律统计表 7 k% w% Q- b O0 o
3 R% L/ U( P6 ^+ f! p5 |! w6 o
) `4 k6 M5 y2 A 不符值范围
) L5 |4 p3 u" K" e( W 点数 % |0 m5 m, ]' R! l
所占比例 : A% w# ~1 A! r3 z1 T/ H
0.0≤|dh|≤0.1
! e. Z' e8 L' k$ U 115 ) i: e/ O8 }5 e, k0 b3 R7 M% Z
34.02% " D% U" K" e/ X* D5 \
0.1<|dh|≤0.2 . `; V8 _- V" @; a3 N, G
64 5 b' t6 h6 `3 d& r7 H5 m3 Z
18.93%
% M- Y& j; j1 l# P, T 0.2<|dh|≤0.3
7 G" Y; Q" D, H l# d 59
7 T" P! ~; H4 g1 H j* R K 17.46% ( V' u7 a# C f2 G( w
0.3<|dh|≤0.4
( P, r1 L: _3 T! O V 35
" e; }& H3 X# M# w+ c: U 10.36%
$ c8 E: ~* S/ z6 r( G3 N# G8 W+ C 0.4<|dh|≤0.5
) f+ J& o1 u- E# r0 ? 17
: N& C1 b K' S4 w8 V; v: U; z4 H 5.03% + l5 k: N. P* i1 g- s. ^
0.5<|dh|≤1.0 , _, ^) J( m% i: R6 X
34 * r3 [/ P V% T
10.06% * u( [3 w2 r- b% h
1.0<|dh|≤2.0
+ B" I& v7 N" |( u1 O 13
$ t: r0 { B/ @ 3.84% - H P% a* r- y8 Z+ K/ ]7 R% h
2.0<|dh| - X. Q) ~$ C) c4 a/ j9 H
1 3 l5 {% E5 z0 }$ s+ n
0.30% 8 _; j5 _/ v% u# z# i w$ x" u
合计
: O. h6 Q. v! n$ J {$ I" j$ w 338 - ] X, S7 J* a' b+ B. v# m
100.00% 6 v1 |& z* ~/ Q3 c! _" w# ^, W
6 C) n. p7 a1 ]; H- Q0 v+ i2 w
8 {3 t' \' j* ~/ ~ 按GPS测高模式处理得到的下半部测区主检测线交叉点不符值分布规律见表4。
2 d" i O7 m/ ~; ?& u2 A 表4 按GPS测高模式处理的下半部测区主检测线不符值分布规律统计表
! o k. x n* i* q4 r9 H1 {4 u( S' W 7 u! Y0 e) ~6 z, v
1 a& a# @* H) {/ q) G7 X
不符值范围 $ H, t. s. b2 s" x7 E
点数 ( a: w0 _8 T) J& a6 n; H8 y
所占比例 * s- X# T, m( [
0.0≤|dh|≤0.1 9 v# }! O2 x8 u$ g* j6 P
59
$ e* V% C; U# f4 \* T7 \0 F 59%
1 ]& Y! d8 K/ E 0.1<|dh|≤0.2
/ e8 n8 Q- P' K u: n% A 27 % Z; G" W2 A+ E. c" C$ C% q. t
27%
. [3 g" z% | P9 V 0.2<|dh|≤0.3
* x0 X; i: z6 o- E5 T 9
6 N! e2 B, t# I' w7 e% @# j9 c 9% : w* ^% O' U/ W* P3 X, G
0.3<|dh|≤0.4
% [/ J( z2 P+ N3 y+ k6 y5 d! c 4
5 M$ {6 f" h2 e4 { 4%
" o2 [8 B2 p9 m! p 0.4<|dh|≤0.5 2 ?" l$ s. A1 J: i
0
A' y5 z: U9 t 0% 5 r$ U! W- y" ~/ X% f% c
0.5<|dh|≤1.0
; U) ~4 K9 @$ A% i 1 1 M0 ?$ ]6 Y1 X# d- {4 U/ p4 T# L V* B
1%
5 a( h5 B1 _" C+ W+ I% d+ ]& r 1.0<|dh|≤2.0 6 H. C) R0 `2 ~ Z2 \% y, M
0 - L$ k+ T0 F7 ]6 p7 w
0% 1 e" j6 U- T x4 R% h8 c9 r
2.0<|dh| # K0 x" G/ c& h: e+ S& {. H$ X
0 : h0 q: n/ t1 v5 u# {* K8 ^
0% , U/ F; c* K5 ^) @/ }0 f6 Y8 a
合计 + @% r* f3 b S0 P; K4 p
100
" M6 B) W: l6 K+ \0 \4 ? 100%
* m) T; Y3 e: i+ C, x . J! n5 O" F9 p9 s! `6 W0 i
5 Z! s& d9 H! Y
采用沙埕、海尾两验潮站潮位数据，按传统作业模式处理得到的下半部测区主检测线交叉点水深不符值分布规律见表5。 % c& x1 U" K* F3 \" U
表5 按传统作业模式处理的下半部测区主检测线不符值分布规律统计表
+ |- W3 x( Y% D, q5 p; G+ Y: P * Z" H& }, v- X; x
! u( e) X7 |6 W, c( j. I( }: l( y! k" B2 z 不符值范围 8 ^) i% K. L/ p* h4 T
点数 ! w9 u0 D6 X" W5 ~* W
所占比例 / b' r1 g2 J; X8 j8 d
0.0≤|dh|≤0.1 6 v# i5 w# ]5 A' h, J
52
7 Z+ u9 U# b. N% @. W3 }7 b 52%
. ~7 B$ o r0 }! [! v! w 0.1<|dh|≤0.2 5 |' z, ]4 h: \5 @3 T
30 5 |. B b6 X: L
30%
% @) h8 U2 d S! x3 y* c 0.2<|dh|≤0.3 [+ B8 ^9 b4 T; G
11 & `7 r6 H2 R4 c1 k# Q; Q! C" P
11%
8 H+ u" m }3 n9 {$ d 0.3<|dh|≤0.4 ; |& [% g* S4 P5 p2 U; a$ f
6 2 x& S9 Y6 D- h+ k: J4 x
6%
* p/ _$ P* u8 ~, t 0.4<|dh|≤0.5 0 c8 p. k* i, o2 g
0 # }' I0 ]6 _, F! k# Y+ n$ h, u
0%
1 y6 E4 K# e8 u6 U8 X) I! M4 ~5 \ 0.5<|dh|≤1.0
* y+ w3 D5 D- o9 w 1 7 C" d: E; h3 J+ s
1%
% b1 c4 y7 w6 G* G$ r 1.0<|dh|≤2.0
! g6 b6 G# J \+ l 0
0 O [, {$ k ^' b* W0 l. d- Z 0%
6 y7 z4 R8 }+ E9 X, w 2.0<|dh|
' O( _& ?( c4 d0 n 0 : a1 P& N) W4 ]& K+ I' V: \
0% 8 M; H" j# z: K a9 Z9 B+ M; I
合计
# O- L; B& s* t/ x1 e1 ~- P 100
/ C/ v% V" {% O- Z" r, [. Y3 Z1 B 100%
; K1 F" a9 a3 f7 U1 r
! W. Q6 K8 w2 v, t: z3 X: f ' G1 }4 F3 F+ O, `% m
根据表2～表5的统计分析，我们可以得出如下结论：GPS测高作业模式获得的水深成果，其内符合精度与传统作业模式获得的水深测量成果内符合精度相当，而且两种作业模式得到的测量成果交叉点不符值的统计规律一致，这充分证明了基于GPS测高的水深测量作业模式原理的可靠性。 h1 [8 K/ P' s/ U* V9 @
⒋ GPS测高作业模式与传统作业模式成果图外部检核比较 5 l; `+ J1 O- b5 s. Z5 Q
为了便于比较，我们在进行图形显示时，有意将测高处理数据的水深成果注记进行了移位，实质上两套成果的测点坐标是重合的。见图5～6。 ; C- H$ R3 ~, C, ?
图5中，红色为传统作业模式成果，黑色为GPS测高作业模式成果，等值线间隔为0.5m，很显然，传统作业模式处理的水深成果普遍浅于GPS测高作业模式处理获得的水深成果。
. s0 \5 ~, q6 @* S9 t' \" G 
' Y4 y9 P8 h0 v5 H% k' k" y @
图5 GPS测高作业模式成果与传统作业模式成果局部比较图 - x. c% _* w! R5 o0 Z" V7 C- H" h% m
+ e- }# h+ Y# J9 {! Q6 r7 P
图6 两种处理方式得到的水深成果的叠加等值线图
+ u- \) p! ]+ z4 l% ~0 R, [ 我们对8月30日、8月31日、9月2日三天的测量数据，将传统作业模式处理得到的成果与GPS测高作业模式处理得到的成果进行了全面比较，其比较结果列于表6。
: a, k! [/ V4 I$ O" h+ H 表6 GPS测高作业模式与传统作业模式处理的水深成果比较
( |# o% r V' L& B+ @' R% ^ " o7 g. X) u5 q8 X3 |" I. x+ m
) A1 u, r' j; ^3 e- g" F' q
日期 ; @& {) W5 C6 P. |& n" O
最小差(m) T3 z, R6 r8 w( n
(传统-GPS) * b, I2 U* ^3 c$ P% g
最大差(m) ' d3 L$ f' t+ s5 E4 M/ F. S' \
(传统-GPS)
3 I/ _. Q7 U" _ X Q- O: H0 m+ z 平均差(m) % K! f* p% x& `2 l; [+ a
(传统-GPS) * i' Z' k: V) K
均方差(m) $ A4 q- W; e- ]) w; e3 s. Y( D
标准差(m)
; X* Q/ e9 C: B& ? 比较测点数 , X% c- P# k4 j
8月30日 1 i% P% w8 l: O Y' a1 j( l
-0.430 * S6 q1 K W: {4 \1 M5 R
0.346 $ n+ f8 O3 |+ r9 {! C& Q1 D$ [1 ?, | n, H
-0.1640 4 ~* [2 {6 Q: |6 D
0.1801
/ R. E5 V/ M4 b/ R 0.0527
* H: c3 X: A" Z* ] 7443 ! l: n) ?3 x7 a, X0 \8 B
8月31日 7 g* K; C% E. P% Q8 [( D. \7 |
-0.744
N& U: E& K2 }1 G9 o4 P 0.268 $ K" ^5 P( ~/ I
-0.0702 : `9 u( t2 {5 |; R
0.0970 $ \6 y. v! ~1 R
0.0473 ; \, h6 ]2 m- j
5091
; N( t- U3 j3 m/ Y 9月2日 5 V% f$ A9 y! z8 U
-0.323
9 {) `* j) a y; {) K, `6 c0 R0 Z) f5 t 0.397
( k2 X2 O9 [. i- [% K7 f9 Y; ] -0.0977
2 t6 Z, e0 I. Z- l$ h9 o/ E 0.1217 " l% K+ O6 }! w; F
0.0513
9 J5 R: i6 ?9 `6 |1 ~ 5750 4 D% r2 G* w- _3 s8 ]
& t$ J5 \ f; h! f* e
9 E) \( t8 K' b9 d0 b" { 从表6可知，按GPS测高作业模式获得的成果水深与传统作业模式获得的成果水深存在7～16cm的系统偏差，如前所述，该系统差主要是由传统作业模式不能准确完成测量船动态吃水改正引起的。 $ @* A. Q* @. T$ U4 V
6 Q+ |3 g7 |9 f/ I: e1 q
- J9 d# ?- L: q4 [' C6 v2 [; T+ D4 O( a7 d7 D
% e; ?6 Y' G ^. G
% \ }$ i$ h) r0 Z& `/ N |
