1、浅析长江口深水航道治理二期工程 GPS (B 级)控制网布设技术上海达华测绘公司 万 军摘要:本文对长江口深水航道二期工程 GPS 控制局域网的结构特点、外业观测成果、平差方法、高程异常网等技术作一番简析。关键词:浅析 GPS 控制网、布设、技术一、前言长江口深水航道治理工程是我国有史以来国内最大的航道治理工程。由于工程范围大、远离陆域,定位和定位精度控制也成为航道治理工程中一大难题,尤其是南、北导堤施工测设定位需要高精度定位。于 1996 年向长江口航道建设有限公司(业主)提出了一系列应用 GPS 定位方案,并提出建立覆盖施工区域的 GPS 工程局域网(B 级)的方案。经测绘专家多次论证和评
2、审,获得了很好评价。在 1997 年 12 月至 2001 年 6 月 分别完成一、二期 GPS 工程局域网,尤其是二期工程 GPS 局域控制网覆盖整个长江口施工区域,其覆盖范围达 800 平方公里,而且该网还包括了高程异常网(HACN) 。目前,长江口施工区已有50 多台 THALES 公司的 DNSP 6502 双频接收机分布在各工段进行定位,如疏浚测量定位、铺排定位、400 多吨的半圆体抛设定位、导堤变形监测等等。从最近的测量资料获悉,南导堤处两段导堤(由两个不同时工单位施工)合拢后,其两条导堤中心轴线仅仅偏差 2cm,这表明 GPS 网对定位精度有了最佳控制,各次各幅的测图的拼接也没有
3、出现断裂,拼接得相当好,这一些都能足以表明长江口深水航道治理工程 GPS 局域网(以下简称:GPS 网)的布设是成功的。为了进一步掌握 GPS 定位技术在水运工程中的应用,本文将对 GPS 网的布设技术进行简析,作一番探讨。二、GPS 控制网建立的三大特点1、 用长基线构建 GPS 骨架网根据长江口沿岸的地形走势、观测环境与条件;交通和通讯的状况、河口水文状况,建立 GPS 网要比陆域或城市建立 GPS 控制网更具特殊性、复杂性,如长江口河口处南、北两岸的宽度达 100 多公里,航道进口端至横沙岛长达 60 多公里,在这样大的范围建立GPS 控制网没有先例可鉴。面对这种情况,提出了采用长基线构
4、成 GPS 骨架网,根据有关理论与公式计算,基线愈长,其相对误差愈小,反之愈大,因此采用长基线作为骨架网,也可满足业主提出“基线相对精度 10-710-8”的要求。为了加强网的整体强度,在长边骨架网之内,布设了加密网,从而加强了网的整体强度,长江口深水航道治理二期工程 GPS控制网的结构如图(1)所示。图(1)长江口深水航道治理工程二期 GPS 局域控制网2、 选用佘山 IGS 站佘山 IGS 站是国家 GPS 网 A 级点,是 IGS 长期跟踪站,它在 WGS-84 系统内具有很高绝对位置精度,充分利用这个重要资源,对建立长江口 GPS 网有着重大的作用。如图(2)所示,从佘山 IGS 站于
5、长江口 GPS 网各控制点进行联测,只要保证全网基线相对精度 10-710-8 的条件下,则可使得网中各点在 WGS84 坐标系中的点位中误差与佘山 IGS 站的点位中误差相匹配,由于佘山点位绝对精度优于0.2m,这样网中各点的 WGS84 坐标绝对精度可达0.2m,大地高差基本上可以控制在 0.3ppm 之内,这是一个很重要的技术特点,使得在基线解算、坐标转换中获得了具有精度较高的 WGS84 坐标的起算点。同时,在计算中采用了 IGS 提供的精密星历,以保证测得的大地高值的稳定性,从而实在大地高系统转换为正常高的过程中,能给与较好的精度保证,为提高整网的精度打下了很好的基础。图(2)从佘山
6、 IGS 站引测网中各点示意图3、重合点的选择在建立 GPS 网中,如何选择重合点也是一个重要的技术,其关系到转换源的坐标精度是否达到需要的精度,所以,在这次长江口建立 GPS 网中,共选用了 17 个点,其中 4个点为国家一等点作为转换源,控制点点号和点名如表(1)所示。控制点点号和点名表 表(1)点号 点名 备注1 花鸟山控制点 国家一等点2 绿华山控制点3 大戢山控制点 国家一等点4 牛皮礁控制点5 北槽中控制点6 江夏村控制点7 倪家鸿控制点 城建二等点8 芦潮港水位站9 横沙岛小白楼10 横沙岛水手搂11 长兴岛水位站12 陈家镇控制点 国家一等点13 陈家宅桥控制点 国家一等点14
7、 吴淞水位站水准点16 大戢山 GPS A 级点17 佘山 IGS 跟踪站 国家 GPS A 级点三 、外业成果本次 GPS 网外业观测使用了 8 台测量型双频 GPS 接收机,其中有 THALES 6502、ASHTECH Z12 SURVEYER GPS 接收机,其标称精度为 ,外业观)15(Dpm测工作按预先编制的技术设计书的要求与规定进行。由于 GPS 网是采用长基线骨架网组成,基线平均长度在 50km 以上,最长基线达140km,而且要求的精度达到 10-7或更高,这样的精度要求,是一般商用软件所不能达到的,因此,选用 GAMIT 软件。GAMIT 软件是麻省理工学院研制的集 GPS
8、 定位定轨于一体的高精度 GPS 数据分析与处理软件,是目前国际上著名的定位定轨软件之一。我国武汉测绘科技大学、同济大学在消化、吸收 GAMIT 软件的基础上作了移植到微机上工作,为高精度GPS 网的数据处理提供了相应的软件支持。表(3)列出了 GAMIT 软件计算结果,其空间最大中误差,Mx=0.0021m My=0.0021m Mz=0.0021m. GAMIT 软件计算结果表 表(3)序号 基线 Dx Dy Dz 闭合差 环长 Km 基线长 PPm1 0001-0006 0.0173 -0.0093 -0.0076 0.0211 253.3 126.65 0.1672 0001-0007
9、 -0.003 0.0046 -0.0031 0.0063 173.1 86.55 0.0733 0001-0008 0.0091 -0.0025 -0.0082 0.0125 160.5 80.25 0.1564 0001-0010 0.0046 -0.0061 -0.008 0.0111 188.7 94.35 0.1185 0001-0015 .-0.0021 0.0124 -0.0417 0.0436 194.5 97.25 0.4486 0001-SHAO 0.0042 0.0032 -0.0048 0.0071 288.5 144.25 0.0497 0006-0007 -0.02
10、03 0.0139 0.0045 0.025 141.9 70.95 0.3528 0006-0008 -0.0081 0.0068 -0.0007 0.0106 197.9 98.95 0.1079 0006-0010 -0.0127 0.0032 -0.0003 0.0131 98.3 49.15 0.26710 0006-0015 -0.0194 0.0217 -0.0341 0.0448 115.6 57.8 0.77511 0006-SHAO -0.013 0.0125 0.0029 0.0183 186.9 93.45 0.19612 0007-0008 0.0121 -0.007
11、1 -0.0052 0.015 56.2 28.1 0.53413 0007-0010 0.0076 -0.0108 -0.0049 0.0141 43.6 21.8 0.64714 0007-0015 0.0009 0.0078 -0.0386 0.0394 29.9 14.95 2.63515 0007-SHAO 0.0073 -0.0014 -0.0017 0.0076 118.8 59.4 0.12816 0008-0010 -0.0046 -0.0036 0.0004 0.0059 99.6 49.8 0.11817 0008-0015 -0.0112 0.0149 -0.0335
12、0.0383 84.9 42.45 0.90218 0008-SHAO -0.0049 0.0058 0.0035 0.0084 135.2 67.6 0.12419 0010-0015 -0.0067 0.0185 0.0337 0.039 22.2 11.1 3.51420 0010-SHAO -0.0003 0.0094 0.0032 0.0099 129.1 64.55 0.15321 0015-SHAO 0.0063 -0.0091 0.037 0.0386 109.9 54.95 0.702为了验证用 GAMIT 软件计算结果的可靠性,给中国测绘研究院用 BERNESE 软件进行单
13、独计算,计算结果如表(4)所示。GAMIT 软件与 BERNESE 软件计算成果比较:见表(4)BERNESE 软件计算结果表- 6 -BERNESE 软件计算结果表 表(4)序号 基线 Dx Dy Dz 闭合差 环长(Km) 基线长(Km) PPm1 0001-0006 0.0169 -0.0117 -0.0126 00241 253.3 126.65 0.1902 0001-0007 -0.0084 -0.0155 -0.0091 0.0198 173.1 86.55 0.2293 0001-0008 0.0009 -0.0189 -0.0115 0.0221 160.5 80.25 0.
14、2754 0001-0010 0.0067 -0.0234 -0.0204 0.0318 188.7 94.35 0.3375 0001-0015 -0.0052 0.0034 -00456 0.046 194.5 97.25 0.4736 0001-SHAO 0.0137 -0.0342 -0.0276 0.046 288.5 144.25 0.3197 0006-0007 -0.0253 -0.0038 0.0035 0.0258 141.9 70.95 0.3648 0006-0008 -0.016 -0.0072 0.0011 0.0176 197.9 98.95 0.1789 000
15、6-0010 -0.0102 -0.0117 -0.0078 0.0174 98.3 49.15 0.35410 0006-0015 -0.0221 0.0151 -0.033 0.0425 115.6 57.8 0.73511 0006-SHAO -0.0032 -0.0225 -0.015 0.0272 186.9 93.45 0.29112 0007-0008 0.0093 -0.0034 -0.0024 0.0102 56.2 28.1 0.36313 0007-0010 0.0151 -0.0079 -0.0113 0.0204 43.6 21.8 0.93614 0007-0015
16、 0.0032 0.0189 -0.0365 0.0412 29.9 14.95 2.75615 0007-SHAO 0.0221 -0.0187 -0.0185 0.0344 118.8 59.4 0.57916 0008-0010 0.058 -0.0045 -0.0089 0.0115 99.6 49.8 0.23117 0008-0015 -0.0061 0.0223 -0.0341 0.0412 84.9 42.45 0.97118 0008-SHAO 0.0128 -0.0153 -0.0161 0.0256 135.2 67.6 0.37919 0010-0015 -0.0119
17、 0.0268 -0.0252 0.0387 22.2 11.1 3.48620 0010-SHAO 0.007 -0.0108 -0.0072 0.0147 129.1 64.55 0.22821 0015-SHAO 0.0189 -0.0376 0.018 0.0458 109.9 54.95 0.833将两者计算结果进行比较,如表(5)所示:GAMIT 软件和 BERNESE 软件计算结果比较表 表(5)Site DX DY DZ Dn De Dh1 0.0101 0.0101 0.0026 0.0006 -0.014 0.0033 -0.0128 0.0262 0.0203 0.002
18、5 -0.0031 0.0354 -0.0082 0.0235 0.0164 0.0015 -0.0057 0.0295 0.0042 0.0016 -0.0009 -0.0003 -0.0044 -0.0016 -0.007 0.0229 0.0097 -0.0039 -0.0061 0.0257 0.0004 0.0008 -0.0006 -0.0007 -0.0007 08 -0.0002 0.0011 0.0015 0.0007 -0.0005 0.00210 -0.0102 0.0205 0.0115 -0.002 -0.0021 -0.02615 -0.0143 0.0317 0.
19、0179 -0.0026 -0.0046 0.03816 -0.0041 0.0685 0.0414 0.0046 -0.033 0.073- 7 -SHAO 0 0 0 0 0 02 -0.0015 0.0152 0.0054 -0.0023 -0.007 0.0159 -0.0162 0.026 0.0118 -0.0058 0 0.03211 -0.0067 0.0141 0.0054 -0.0034 -0.0017 0.01612 -0.0015 0.0067 0.0021 -0.0016 -0.0023 0.00613 -0.0059 0.0148 0.0049 -0.0039 -0
20、.0028 0.01614 -0.0196 0.0261 0.0205 0.0006 0.0031 0.038由表可见:GAMIT 和 BERNESE 坐标直接比较,两者最大差异:Dx -0.0196mDy0.0685m Dz0.0414m此外,还进行复测基线较差计算,异步环闭合差计算,它们是 GPS 网精度的技术指标之一,其计算结果见表(6)、表(7)所示。复测基线较差精度统计表 表(6)基线名 测站名 118 计算结果 120 计算结果 差值(m) 平均值(m) m 基m基2相对精度(ppm)0001-0006 花鸟山江夏村 126674.2080 0.0038126674.1797 0.
21、00300.0283 126674.1938 0.133 0.3724 2.910-82.310 -80001-0007 花鸟山倪家鸿 86569.6309 0.0036 86569.6327 0.0030 -0.0018 86569.6318 0.091 0.25662 4.1510-83.410 -80001-0008 花鸟山芦潮港 80271.6794 0.0040 80271.6635 0.0025 0.0159 80271.6714 0.085 0.2397 4.910-83.110 -80001-0010 花鸟山横沙基站 94337.5828 0.0034 94337.5758 0
22、.0029 0.007 94337.5793 0.099 0.2792 3.610-83.110 -80001-0015 花鸟山三甲港 97234.0667 0.0036 97234.0715 0.0030 -0.0048 97234.0691 0.102 0.2856 3.710-83.010 -80001-SHAO 花鸟山佘山 144242.2834 0.0049144242.2594 0.00330.024 144242.2714 0.149 0.4172 3.310-82.310 -80006-0007 江夏村倪家鸿 70949.9702 0.0027 70949.9629 0.002
23、6 0.0073 70949.9665 0.714 0.0138 3.810-83.110 -80006-0008 江夏村芦潮港 98930.9089 0.0029 98930.9069 0.0027 0.002 98930.9079 0.103 0.2884 2.910-80006-0010 江夏村横沙基站 49153.7159 0.0026 49153.7150 0.0025 0.0009 49153.7154 0.054 0.1512 5.010-80006-0015 江夏村三甲港 57817.0308 0.0028 57816.9968 0.0026 0.034 57817.0138
24、0.062 0.1736 4.410-80006-SHAO 江夏村佘山 93463.7680 0.004393463.7719 0.0029-0.0039 93463.7700 0.098 0.2744 4.610-8310-80007-0008 倪家鸿芦潮港 28077.0495 0.0025 28077.0569 0.0025 -0.0074 27077.0532 0.033 0.0929 8.910-70007-0010 倪家鸿横沙基站 21824.4829 0.0024 21824.4763 0.0024 0.0066 21824.4796 0.026 0.0728 1.010-700
25、07-0015 倪家鸿三甲港 14948.7874 0.0023 14948.8105 0.0023 -0.0231 14948.7990 0.019 0.0532 1.510-70007-SHAO 倪家鸿佘山 59377.6982 0.0040 59377.6724 0.0027 0.0258 59377.6853 0.064 0.1792 6.710-84.510 -60008-0010 芦潮港横沙基站 49781.9209 0.002649781.9202 0.00250.0007 49781.9205 0.054 0.1512 5.210-80008-0015 芦潮港三甲港 42438
26、.3990 0.0026 42438.4359 0.0025 -0.0369 42438.4174 0.047 0.1316 5.810-80008-SHAO 芦潮港佘山 67594.9734 0.0040 67594.9673 0.0032 0.0061 67594.9703 0.073 0.2044 5.910-84.710 -80010-0015 横沙基站三甲港 11098.8909 0.0023 11098.8609 0.0023 0.0300 11098.8759 0.016 0.0448 2.010-70010-SHAO 横沙基站佘山 64538.0916 0.0040 64538
27、.0719 0.0027 0.0197 64538.0817 0.069 0.1932 6.010-84.210 -80015-SHAO 三甲港佘山 5494.3232 0.004054941.3167 0.00270.0065 54941.3200 0.059 0.1652 7.210-84.910 -8- 8 -异步环闭合差计算表 表(7)异 步 环 闭合差(m) 环路长度 (m)限差(m) 基n3相对精度(ppm)花鸟山陈家镇江夏村倪家鸿花鸟山 0.022 295016.5526 0.248 0.07花鸟山大戢山北槽中江夏村 花鸟山 0.022 282886.5152 0.238 0.0
28、7大戢山江夏村 倪家鸿北槽中大戢山 0.027 251155.7861 0.214 0.107北槽中江夏村倪家鸿花鸟山北槽中 0.021 292570.5934 0.246 0.07江夏村倪家鸿花鸟山大戢山江夏村 0.022 217057.8826 0.187 0.101倪家鸿花鸟山大戢山北槽中倪家鸿 0.0064 208202.4625 0.181 0.03复测基线的较差能反映出 GPS 测量的总体精度,计算结果表明,基线相对中误差均在 2.310-88.910 -8 之间,短基线也在 1.010-72.0 10-7 之间,复测基线的较差全部小于 (m 基 (5mm+1ppmD) ) ,从而
29、表明了 GPS 测量的总体精度是较高的。异基2步环闭合差均在限差之内。四、高程异常网的建立针对水运工程测量的要求,开发与利用 GPS 三维定位功能,将测量定位、测深、验潮(即包括测深归算与改正)即 RTK GPS 接收集于一体。同时也是为了适应于国际疏浚新技术的要求,于是在长江口施工区域布设高程异常网,将测得的大地高按 H 大 =h 正 + 计算出高程异常,然后进行测深归算与改正,获得图载水深,是一个现代 GPS 高程在水运工程中的应用。长江口高程异常网是利用长江口各水位站,如吴淞水位站、牛皮礁水位站、长兴水位站、横沙水位站、北槽中水位站、大戢山水位站、绿华山水位站、芦潮港水位站等构成。如图(
30、)所示。在这些水位站上有较完整的几何水准点、水尺零点,这对求算各站的高程异常值提供良好的基础条件,此外为了求得陆域较高精确的高程异常值,进行了几何水准联测,如芦潮港水文站等与基岩点进行二等水准联测,获得精确的水准高。由高程系统理论可知,大地高 H 与正常高 h 关系为: 式中 为高程异H常,即似大地水准面至参考椭球面之间的距离,所以,要使 GPS 高程 H 在工程测量中应用,必须将其转换为正常高 h,而其中的关键是求定点的高程异常 ,也就是说 GPS 高程拟合,实质上就是求定的 值。长江口各点高程异常变化趋势如图(4)所示。根据观测数据处理的结果表明:- 9 -图(4)长江口各点高程异常值变化
31、趋势示意图佘山至花鸟山之间的高程异常值之间差( )为 4.824,佘山至吴淞之间的高程异常值之间的差( )为 1.179,横沙至北槽中之间的高程异常值之间的差( )为 0.54,北槽中至牛皮礁之间的高程异常值( )为 1.342 等。从图(4)中可看出长江口区高程异常的趋势变化,所以在 GPS 高程拟合计算中,既考虑了高程异常的趋势性,即长江口的高程异常不是一个“抛物单曲面” ,而是呈鞍形和波浪形。若用一个二阶多项式函数计算,势必其结果出现较大误差,而在凹凸面变换处残差更大,内插也不会令人满意。因此,采用加权平均值的方法。加权平均值的思路:由于似大地水准面起伏和复杂,因此,高程异常值的分布也是
32、复杂的。实际情况表明,相邻点只见到高程异常值的变化基本上是一个连续的渐变过程,所以当用一点周围各点的高程异常值估算该电高程异常时,周围点距距该点越近,其高程异常值对该点高程异常的影响越大,距该点的距离越远,影响越少,给予这种思路便提出了权的问题,即认为在用某点周围若干点的高程异常推算该点高程异常时,周围各点的高程异常值推算时所占的权重与各点到该点的距离成反比。假设有一点 A,其周围有几个已知点的高程异常值( ) ,这些点到 A 点水平距离为 Di,则各点对 A 点高程异常影响的权重iPi 可用 来描述。PiA/为此,利用大网中 6 到 7 个已知水准点的正常高和高程异常值推算位于已知点所围成的
33、多边形之内的待定点的高程异常值,利用加权平均值求得的高程异常值求得的高程异常值准确度要比其他方法要略高一些。这就是在长江口控制网中用加权平均值法将 GPS 大地高程转换正常高的重要应用条件之一。几年来的使用发现,GPS 拟合高程精度在很大程度上取决于测取高程异常变化的密度、点数及点的分布情况,基于 800km2,东西跨度达 60km 的测区,可以进一步考虑分割成若干小区,建立数字高程异常模拟,使拟合残差达到较高精度,这样高程异常网的实用性更强,应用范围更广。五、平差计算模型与分析- 10 -本次 GPS 网平差计算采用了由同济大学王解先老师提出的一种简便实用的 GPS 网平差模型,该模型是以
34、GPS 基线向量作为观测值,而平差参数直接选为我们感兴趣的高斯平面坐标和大地高,这样将给其它附加计算带来明显的方便,如地面归心测量和其他地面上的常规观测数据,另外,如在法方程中消去大地高参数,就可变成平差模型,如在法方程中消去高斯坐标参数,就可变成高程平差模型。(1)平差模型以下标 I 表示点号, 表示 GPS 观测点的空间直角坐标, 表示TZYXR R由 GAMIT 软件求出的基线向量(空间坐标差) ,则第 j 条基线的误差方程可表示为:(1)12iijjv式中的 V 为改正数,i1 和 i2 为该基线的端点点号, (1)式中的权阵为 ,它是jPGAMIT 解算软件求得的基线协因数阵的逆阵。
35、以 表示 GPS 点的空间直角坐标近似值,将(1)式中 R 化为:0R(2)R0将(2)式中的空间直角坐标改正数 化为高斯坐标和大地高的改正数 的形式:Tggzyxrr(3))(测得的控制网成果必须符合到已有的国家 54 坐标系中,一般而言,直接由 GPS 观测结果根据选定的中央子午线投影出来的高斯坐标与当地已有坐标之间存在平移、旋转、缩放,各点与旋转参考点坐标差的关系可由下式表示:(4)00)(1yxkyxgg上式中的 和 为缩放系数和旋转角, 和 为当地坐标系中的坐标,k为旋转参考点的高斯坐标,它在旋转前后坐标不变,它可以是 GPS 网中的任意0yx一个已知点。(5)cossini将(2)
36、 、 (3) 、 (4) 、 (5)带入(1)式,得: 1122 0cossini0cssiio igiigi HyxrRHyxrRVj kyxyxr iggi 2002 0sinco)(- 11 -(6)iiggil lkyxyxrR 1001 0sincoi)(6)式简写为:(7)ljkErDVj 12(7)式的权阵仍为 ,常数项为:jP(8)0012kRiijj式中的 和 为缩放和旋转角的近似值。0k由(7)式组成法方程,在所有的基线处理完毕后,得到最后的法方程:(9)CxN其维数为 (n 为 GPS 网的总点数) ,N 和 C 为最后法方程的系数和常数项,未知23数 包括 n 个未知点
37、的平面坐标和高程改正数、1 个尺度因子和 1 个旋转角参数。将已知点高斯坐标对应的行列从(9)式中消去,并选择一个大地高的点(一般可选为(4)式中的旋转参考点) ,将其对应的行列也从(9)式中消去。解(9)式便得最后的结果。(2)这样的模型将用到空间坐标和大地高的偏导数,而这个偏导数 的数学表达方式十rR分复杂,为此引入数值方法来解决该偏导数的计算。因为空间坐标与高斯坐标之间存在着一一对应关系,所以其在理论上是正确的 ,便将(3)式改写为:(10)igii zyxazyx32311式中的下标 i 表示 GPS 控制网中的 i 点。由 i 点的标准高斯坐标和大地高近似值 ,算出相应的空间直角Ti
38、giii Hyxr00坐标值 ,然后将 变化一个微小量 d,算出变化后的空间直角Tiiii ZYXR00gi坐标值,并求出其变化量 ,从而求得(10)中的:(11)dZada31211,同理,通过对 和 的微小变化,可求出偏导数的其它值。0giyiH按导数的定义,当 d 趋近于零时,以上方法求得的导数是严密的,但这在数值计算上是不可能的,也是没有必要的,因为平差计算总是要迭代的,偏导数准确到前两位就已足够,通过一些实际计算,d 取 0.1 米是比较适当的,它既不使计算机计算时出现数值问题,也使得求出的偏导数具有足够的精度。- 12 -六、结束语通过上述的技术简析,反映出 GPS 定位技术在水运
39、工程中应用的新特点,也正因为由于这些新特点的作用,使长江口区二期工程 GPS 网获得成功。也为我国水运工程测量提供宝贵而丰富的新技术。同时,我们也应当看到,在长江口的现代 GPS 高程测量仅是跨出了可喜的一步,面对长江口如此大的区域,如何获得高程异常、高精度,还有许多事仍须进一步研究,如在求出各点的高程异常后,将其以一定的数据结构存储,建立高程异常数学模型,并利用虚拟技术将该地区的高程异常变化直观的反映出来。不仅如此,而且在通过 GPS 验潮,将测深归算与改正、测船的姿态改正应用中仍有很大的潜力。这一些是仍为继续为长江口深水航道工程测量技术发展的基本要求。随着长江口深水治理工程与上海洋山深水港建设日益加快,无论是目前的建设,还是今后的运营,都面临着构建数字化港口、数字航道信息更新、获取的新课题。特别是,近几年来,飞速发展的 RTK-DGPS 技术,如 RTK-VRS,将成为长江口洋山港航道空间管理获得航道地理空间信息的重要手段之一。因此,建立 RTK 虚拟基准站网将可能替代现有的 GPS 控制网。