论文推荐| 章传银:高程系统定义分析与高精度GNSS代替水准算法

测绘学报2017.9.14

测绘学报

构建与学术的桥梁 拉近与权威的距离

测绘地理信息与导航高端论坛通知(第二号)

高程系统定义分析与高精度GNSS代替水准算法

章传银, 蒋涛, 柯宝贵, 王伟

中国测绘科学研究院, 北京 100830

收稿日期:2017-02-10;修回日期:2017-07-21

基金项目:国家自然科学基金(41374081;41674024)

第一作者简介:章传银(1968-), 男, 博士, 研究员, 研究方向为大地测量学. E-mail: zhangchy@casm.ac.cn

通信作者:蒋涛, E-mail: jiangtao@casm.ac.cn

摘要:从高程系统定义出发,探讨高程基准面的重力等位性质,测试分析不同类型高程系统地面点高程之间的差异,考察GNSS代替水准与实际水准测量成果的一致性,进而提出新的GNSS代替水准算法。主要结论包括:① 当精度要求达到厘米级水平时,正常高的基准面也应是大地水准面。中国国家1985高程基准采用正常高系统,其高程基准面是过青岛零点的大地水准面。② 近地空间中等解析正高面与大地水准面平行,GNSS代替水准能直接测定地面点的解析正高,但正常高系统更有利于描述地势和地形起伏。③ 本文给出的GNSS代替水准测定近地点正常高算法,大地高误差对正常高结果的影响比大地水准面误差大,前者影响约为后者的1.5倍。

关键词:GNSS代替水准 高程系统 高程基准面

The Analysis of Height System Definition and the High Precision GNSS Replacing Leveling Method

ZHANG Chuanyin, JIANG Tao, KE Baogui, WANG Wei

Abstract: Based on the definition of height system, the gravitational equipotential property of height datum surface is discussed in this paper, differences of the heights at ground points that defined in different height systems are tested and analyzed as well. A new method for replacing leveling using GNSS is proposed to ensure the consistency between GNSS replacing leveling and spirit leveling at mm accuracy level. The main conclusions include:① For determining normal height at centimeter accuracy level, the datum surface of normal height should be the geoid. The 1985 national height datum of China adopts normal height system, its datum surface is the geoid passing the Qingdao zero point.② The surface of equi-orthometric height in the near earth space is parallel to the geoid. The combination of GNSS precise positioning and geoid model can be directly used for orthometric height determination. However, the normal height system is more advantageous for describing the terrain and relief.③ Based on the proposed method of GNSS replacing leveling, the errors in geodetic height affect more on normal height result than the errors of geoid model, the former is about 1.5 times of the latter.

Key words: GNSS replacing leveling height system height datum surface

1 若干定义1.1 正(常)高的力学定义

地面点的正高h定义为

(1)式中,gg分别为地面点到大地水准面间的流动点重力和平均重力;WW0 (2)为地面点的实测重力;G为万有引力常数;ρ为地壳密度;∂γ/∂h为正常重力梯度。利用式(2),由水准测量高差和水准路线上的重力值计算得到的正高,称为赫尔默特(Helmert)正高H*[1, 2, 7] (3)为地面点到大地水准面间的平均正常重力。

1.2 解析正高与Molodensky正常高

1.2.1 解析正高

若直接由地球外部重力场数据,按解析延拓法计算地面点到大地水准面间的流动点重力g,则平均重力可按式(4) 确定

(4)

本文将按解析平均重力g′确定的正高暂称为解析正高h′。由式(4) 得解析正高定义如下

(5)

与Helmert正高相比,解析正高可适用于近地空间的非地面点,利用解析正高公式,易于建立不同类型高程之间的重力学关系。

1.2.2 Molodensky正常高

图 1所示,从地面点A向下量取其高程异常ζ后得到似地形面上的Q点(在一定海拔高度的大陆地区,ζ<0,似地形面在地面上方)。令Q到对应地球椭球面上G点之间的平均正常重力为γG,用γ替换式(1) 或式(3) 中A点到大地水准面上O点之间的平均(正常)重力,得到Molodensky理论中重力学定义[3]的正常高hγ图 1 Molodensky正常高的几何意义Fig. 1 Geometric interpretation for Molodensky normal height图选项设U000-UQ (6)式中,U0式(6) 中,通过引入Molodensky基本条件,只是将计算重力位数c的积分变量由难以测定的实际重力g变换为可解析计算正常重力γ,但仍然有c=W0-W

由式(5) 和式(6) 可得

(7)

式中, HA点的大地高;ζA点的高程异常;NA点处的大地水准面高。

由于γγH (8)

本文测试不同类型高程系统地面点高程之间的差异,分析GNSS代替水准与实际水准测量成果的一致性,研究高精度GNSS代替水准算法。

2 高程基准面的重力等位性质

传统的高程基准由水准原点(高程基准零点)和高程控制网实现,水准点正(常)高的起算面称为高程基准面或零高程面。为保证GNSS代替水准测定的正(常)高与水准测量成果一致,要求GNSS代替水准中采用的高程基准面与水准点正(常)高的起算面完全一致。

2.1 由高程系统定义分析高程基准面的等位性质

h00h=c在区域高程基准中,水准点的重力位数c等于高程基准零点重力位WR与该点重力位W之差,水准点高程的起算重力位等于高程基准零点重力位W对于任一地面点,不论是用正高、正常高还是力高来表示其海拔高程,该点的重力位数cc都是唯一的,高程基准面的重力位处处都等于高程基准零点重力位,高程基准面的这种重力等位性质与采用的高程系统类型无关。

总之,不论是正高系统、正常高系统、力高系统,还是重力位数系统,高程基准面是唯一不变的,都是过高程基准零点的大地水准面。

2.2 非地面点的高程异常及其正常高起算面

由式(3) 或式(6) 可知,正常高定义适合非地面点,而由重力方法确定的似大地水准面是相对计算者选定、用于表示地面的特定数字高程模型而言。设在似大地水准面模型计算时用2′×2′数字高程模型代表地面,若某地面点A的正常高为h,由2′×2′数字高程模型内插得到该点在数字高程模型面上的正常高为h0,由对应的2′×2′似大地水准面模型内插得到该点在数字高程模型面上的高程异常为ζ0 (9)

式中,Tδgγ分别为地面点到对应数字高程模型面间流动点的扰动位、扰动重力和正常重力。

中国大陆2′×2′数字高程模型的高程最大值一般只有6 km,而珠峰地区的实际地面高程超过8 km。若令δg=-100 mGal,γ≈10 m/s2,Δh=2000 m,则Δζ=100×10-5/10×2000 m=20 cm。即在珠峰地区,若将2′×2′重力似大地水准面直接作为正常高的基准面,则由GNSS代替水准测定的正常高要比实际正常高大20 cm。

A点不在地面上,如在地面下方或低空高度上,Δζ的绝对值还会进一步增大。因此,将似大地水准面作为正常高的起算面是不严密的。

综上所述,在严格意义上所有类型高程系统的基准面都是重力等位面。中国国家1985高程基准采用正常高系统,高程基准面是过青岛零点的大地水准面。

3.1 地面点正常高h*任意地面点的正常高h* (10)式中,∂γ/∂h=0.308 6×10-5s-2为正常重力梯度。在一定海拔高度的大陆地区,ζ<0,因此似地形面在地面上方,由式(10) 可得:hγζ=-50 m,h=6000 m,则h-hγ=-0.308 6×10-6×50×6000 m≈-0.09 m。Molodensky正常高h利用Eigen6C4地球重力场模型分别按Molodensky正常高定义hγ=(W0G*0*计算2′×2′数字高程模型面上的正常高值,得到hγ*结果见图 2图 2 地面点Molodensky正常高h与正常高h*之差Fig. 2 Differences between Molodensky normal height hγ (11)式中,γγi+1分别为水准点ii+1在地球椭球面上的正常重力(单位为mGal);hi为水准点i的概略高程;hm=(hi+hi+1)/2;γ式(11) 中,γ=γi-0.308 6hi/2,是水准点i处似地形面到地球椭球面间正常重力平均值的线性近似,可见在中国国家水准网成果中,水准点的正常高属于Molodensky正常高系统。

3.2 解析正高h′与H-N的差别测试分析

在中国大陆东部、中部北方和南方、西部北方和南方各取1点,分别从大地水准面起算,以大地高每增加100 m为考察点,用Eigen6C4重力场模型计算考察点的重力位W、大地水准面的重力位W00-W/g′。以h=H-Nx轴,以h与解析正高h′之差h-h′为y轴绘图,见图 3图 3 GNSS直接代替水准测定的正高h与解析正高h′的差别Fig. 3 Differences between orthometric height h determined by GNSS replacing leveling and actual orthometric height h图选项由图 3可知,在整个中国大陆地区,GNSS直接代替水准测定的正高h=H-N与解析正高h′相差不到2 mm。忽略由地球重力场模型计算其对应大地水准面上重力位W0的计算误差,则有:h′=H-N。可见,等解析正高面与大地水准面平行,GNSS代替水准能直接测定地面点的解析正高。 (12)

式中,Δg为地面点到大地水准面间的平均空间重力异常。

令Δg=-100 mGal,h=6000 m,则h-hγ=100×10-6×6000 m=0.6 m。在中国大陆东部、中部北方和南方、西部北方和南方各取1点,由Eigen6C4地球重力场模型按定义计算Molodensky正常高hγ0G和有关模型参数。令h=H-N,以hx轴,以h与Molodensky正常高hγ之差h-hγy轴绘图,见图 4图 4 Molodensky正常高h与(H-N)的差别Fig. 4 Differences between Molodensky normal height h图 4可知,h-h随高度增大而增大,在大地水准面附近接近零,在中国西部北方超过0.5 m;在中国东部地区,h-hγ一般小于0.15 m。

实际上,由式(7) 可得

(13)

4 不同类型高程的水准面不平行性4.1 GNSS代替水准中水准面不平行改正的通用形式

设GNSS定位测得地面点大地高为H,地面点到大地水准面的大地高差为h=H-N,由地球重力场知识可知,该点的重力位W=T+U=ζγ+U,其中,TζU分别为地面点的扰动位、高程异常和正常重力位,顾及通用高程定义

(15)

类比水准测量中正常高的水准面不平行改正,这里仍将∈作为GNSS代替水准中通用高程的水准面不平行改正。

GG=ghd(力高中G为常量,力高系统可看成是与重力位数系统等价的一种特殊类型正常高系统),都需要增加水准面不平行改正∈。由于大地水准面重力位W0已知,地面点高程异常ζ可由大地水准面模型按Possion积分计算,正常重力位U、正常重力γ和平均正常重力γ*G)可严密计算,因此在GNSS代替水准中,正常高(hG)水准面不平行改正∈能由GNSS定位和大地水准面模型按式(15) 计算。的水准面不平行改正也难以准确计算。

若用GNSS代替水准传递高差,也需对高差增加水准面不平行改正Δ∈。即

(16)

4.2 地面点正高与正常高力学性质的差异

图 2可知,等解析正高面相互之间平行,且都与大地水准面平行;但由图 3可以看出,等正常高面相互之间不平行,等正常高面与大地水准面也不平行。利用等高程面上的重力位差大小和空间分布情况可以评价高程系统的力学性质。图 5是用Eigen6C4重力场模型计算的3000 m等正高面的重力位(图 5(a))和3000 m等正常高面的重力位(图 5(b))。

图 5可知,等正高面的重力位(已扣除了平均值)随空间变化复杂,等正常高面沿卯酉圈方向具有等位性质。考察卯酉圈上正常高等于3000 m的两点AB,其重力位差ΔWAB≈0.308 6×10-5ΔNAB≈9.2×10-3ΔNAB m22NAB为两点处大地水准面高之差),是小量。图 5(b)中各卯酉圈方向的重力位标准差均小于0.03 m22。可见,正常高系统更有利于描述地势和地形起伏。图 5 重力位对比图Fig. 5 Gravitational potential contrast图选项


5 GNSS代替水准测定正常高的新算法5.1 GNSS代替水准测定近地点正常高的统一算法

当忽略大陆水负荷、非潮汐海平面和非潮汐大气负荷等变化引起的大地水准面形变(GRACE监测结果显示,非潮汐大地水准面形变达到厘米量级)时,大地水准面高不随时间变化,是稳态的。本文只讨论稳态大地水准面情况下GNSS代替水准测定正常高算法。

设GNSS测定近地点A(可以是地面点,也可以是非地面点)的大地高为H,忽略厘米量级的地心运动影响,则GNSS代替水准测定正常高h* (17)式中,h=hhγγα*γG为区域高程基准位差。 (18)为大地水准面重力位;γU分别为A点正常重力与正常重力位;ζA点高程异常,能由大地水准面模型按Possion公式精确计算。GNSS大地高HζγU自变量的形式代入。的高精度计算。

5.2 高精度GNSS代替水准测定正常高的误差分析

对GNSS代替水准测定正常高的通用算法式(18) 进行变分,省略γ上标,顾及高程基准位差ΔW00r (19)

式中,δγδU分别为A点正常重力变分和正常重力位变分;δγ为平均正常重力变分;δζA点高程异常变分。

对于误差分析目的,不失一般性,可采用如下线性近似

(20)对式(20) 进一步进行变分,顾及,得 (21)将式(21) 代入式(19),并顾及γ,得 (22)≈0.308 6×10-5s-2ζ=-70 m,γ-≈9.8 m2,因此式(22) 右边第2项可以忽略,即 (23)式(23) 中,GNSS代替水准测定正常高(hh)的误差δhα由两项构成:大地高误差δH和大地水准面误差δN,误差影响因子分别为τH≈1.5(<1.5) 和τ≈1(<1)。这表明,1倍大地高误差会导致近1.5倍的正常高误差,1倍大地水准面误差会导致近1倍的正常高误差。在GNSS代替水准测定力高d (24)

6 小结

本文从高程系统定义出发,探讨高程基准面的重力等位性质,测试分析不同类型高程系统地面点高程之间的差异,考察GNSS代替水准与实际水准测量成果的一致性,进而提出新的GNSS代替水准算法。主要结论包括:

(1) 当精度要求达到厘米级水平时,正常高的基准面也应是大地水准面。中国国家1985高程基准采用正常高系统,其高程基准面是过青岛零点的大地水准面。

(2) 近地空间中等解析正高面与大地水准面平行,GNSS代替水准能直接测定地面点的解析正高,但正常高系统更有利于描述地势和地形起伏。

(3) 本文给出的GNSS代替水准测定近地点正常高算法,大地高误差对正常高结果的影响比大地水准面误差大,前者影响约为后者的1.5倍。

致谢: 本文是在许厚泽院士、宁津生院士、晁定波教授、魏子卿院士、杨元喜院士、李建成院士、吴晓平教授、黄谟涛高工、李斐教授、翟国君高工、张传定教授、申文斌教授、罗志才教授与几位审稿专家的许多建设性意见基础上逐步完善的,在此特别致谢。

【引文格式】章传银,蒋涛,柯宝贵,等。高程系统定义分析与高精度GNSS代替水准算法[J]. 测绘学报,2017,46(8):945-951. DOI: 10.11947/j.AGCS.2017.20170058

更多精彩内容:

李志林:博士研究生学习--从技能到智慧的全面提升

2017年测绘科技进步奖评选结果公告(附名单)

泰晤士发布2018世界大学排行榜,看看有多少地学院校上榜?

重磅| 地学领域的44个国家重点实验室,都在这里了!

《测绘学报》2017年第8期网刊发布

关于坐标系统的知识要点

书讯| 《海量精细点云数据组织与管理》王晏民 等 著

2017年国家自然科学基金立项数和立项金额(完整名单)出炉!

2017测绘地理信息前沿技术论坛(无人机测绘技术及应用专场)在西宁召开

武汉大学“最奢侈测绘基础课”!众院士坚持了20年!

李德仁: 测绘遥感领域十个前沿技术都在这里了!

论文推荐| 王乐洋:Partial EIV模型的非负最小二乘方差分量估计

科技部发布“重大自然灾害监测预警与防范”重点专项2017年度项目申报指南

《战狼2》撤侨到底发生在哪?破译电影里一幅近代地图的秘密

2017无人机测绘技术及应用论坛在线报名已开通

惊艳的骗局:原来你每天坐的地铁线路图这样欺骗你的双眼!

权威 | 专业 | 学术 | 前沿

微信投稿邮箱 | song_qi_fan@163.com

欢迎加入《测绘学报》作者QQ群: 297834524

最新评论

    查看更多评论
    加载失败,请点击重试
    已加载全部评论
    测绘学报

    关注《测绘学报》官方一点号,及时掌握测绘地理信息领域前沿科研信息,随时随地浏览优秀论文摘要,为读者、作者、专家提供一个自由互动的学术交流平台。

    +关注