远离大陆海岛的高程传递_柯宝贵

2011年第12期柯宝贵，等：远离大陆海岛的高程传递

中图分类号：P216

文献标识码：B

3

0911（2011）12-0003-02文章编号：0494-

远离大陆海岛的高程传递

柯宝贵，章传银，张利明

（中国测绘科学研究院，北京100830）

Height Transmission Far from the Mainland

KE Baogui ，ZHANG Chuanyin ，ZHANG Liming

摘要：根据舟山庙子湖岛短时验潮站上GPS 点的大地高，结合EGM2008重力场模型计算短时验潮站的重力大地水准面高，得到

GPS 点的海拔高。利用同步水位观测技术得到GPS 点多年平均海面的1985国家高程。将GPS 海拔高与其1985国家高程比较，发现两者相差不超过10cm 。研究结果表明，在没有前期观测数据的情况下，同步水位观测技术和GPS /重力移去-恢复法均可进行远离大陆海岛礁的高程传递。

关键词：重力场模型；1985国家高程；高程传递；同步水位观测

一、引言二、远离大陆海岛高程的确定方法

1．移去-恢复法

利用移去-恢复法确定大地水准面，就是将大地水准面高分解为中长波分量和残余分量两部分。

中长波分量用地球重力场参考模型（EGM2008）计算，通常称为模型大地水准面高，文中以N GM 表示；同时利用参考模型计算模型重力异常，文中以Δg GM 表示。首先从观测的重力异常中扣除Δg GM 后，利用Stokes 公式计算大地水准面的残余分量，此过程为移去的过程；然后将模型重力大地水准面高与残余分量相加即可得到研究区域的大地水准面高，此即恢复过程。计算方法如式（1）

R

N =N GM +（Δg －Δg GM ）S （ψ）d σ

4πγσ

我国是世界上海岛较多的国家之一，海岛礁是海洋经济的重要领域，是优化海洋经济发展布局的重要载体，对其开发利用具有广阔的发展前景。因

我国远离大陆海岛礁的大地基准和高为历史原因，

程基准并未与陆地基准统一起来。而海岛礁的测

如何有效地进行海岛高程传递一绘有其自身特点，

直是海洋测绘中的难点问题

［1］

。目前，国家测绘地

理信息局正在组织开展海岛礁测绘工作，其目的之

一是构建陆海统一高程/深度基准，实现电子地图与电子海图的无缝拼接。

海岛高程基准传递有4种传统方法：静力水准

GPS 水准法及常规大地测量法。李法、动力水准法、建成比较分析了静力水准法、动力水准法及常规大

［2］

地测量法3种方法的优缺点；章传银等对EGM2008模型在中国内地的适用性进行了研究［3］；张利明研究了利用基于重力等位面的方法进行跨

［4］

海高程传递的问题。本文基于在庙子湖岛建立的短时验潮站，利用GPS /重力移去-恢复法实现远离大陆海岛的高程传递，得到庙子湖岛短时验潮站并与基于同步水位观测技术的1985国家的海拔高，

GPS /高程基准传递方法进行对比分析。结果表明，重力移去-恢复法可进行远离大陆海岛的高程传递，

实现陆海高程基准的统一，为庙子湖岛的各种测绘信息采集提供统一的1985国家高程基准。

（1）

N GM 称为模型大地水准面高；Δg GM 称为模型重式中，

力异常。

N GM

GM a =∑r γn =2r

N max

N max

()∑（C

m =0

n

n

nm

cos m λ+

（2）

nm

S nm sin m λ）P nm （cos θ）

Δg GM

GM

=2∑（n －1）r n =2

()∑（C

a r

m =0

n

n

cos m λ+

（3）

S nm sin m λ）P nm （cos θ）

（r ，式中，θ，λ）为计算点的球坐标（地心距离、余纬和经度）；GM 是地心引力常数；a 是地球椭球长半轴；C nm 和S nm 为完全规格化位系数；P nm （cos θ）为完全规格化缔合Legendre 函数；N max 为位模型球谐展

05-16收稿日期：2011-2009AA121402）；国家自然科学基金资助项目（40974016，41004004）；中国测绘基金项目：国家高技术发展计划资助项目（2009AA121405，

科学研究院基本业务费资助项目

作者简介：柯宝贵（1982—），男，湖北黄冈人，博士，主要从事大地测量与地球动力学研究工作。

开的最高阶数；γ为计算点的正常重力；σ为单位球面；ψ为计算点到积分面元之间的角距，式（1）中的S （ψ）为Stokes 函数，其表达式为

1ψ

S （ψ）=－6sin +1－5cos ψ－

2ψ

sin

2

3cos ψln sin

(

ψψ+sin 222

)

（4）

图1

同步改正法

GPS 测量可以得到精确的大地高。根据关系式h =H －N （式中，h 为观测点海拔高程；H 为大地高；

N 为大地水准面起伏），即得到利用海岛GPS 观测点的大地高H 和地球重力扰动位模型得到观测点的海拔高h 。

2．CGCS2000下的大地水准面转换

由于EGM2008地球重力场模型与GPS 大地高（CGCS2000）所用的参考椭球不一致，故需要将EGM2008参考椭球面的物理量转换到CGCS2000的参考椭球面上，利用EGM2008参考椭球所用椭球参数，按照式（6）进行转换

a

d h =－W d a +（1－f ）sin 2φd f

W

三、观测数据分析

1．实际观测数据

构建了3个临为了实现陆海大地基准的统一，

时GPS 基准站，并与IGS 站或浙江CORS 网络联网。第一个建在庙子湖岛南边庙里村，第二个建在东福

GPS 基准站24h 第三个建在朱家尖。试验期间，山，

连续观测不少于60d ，采样间隔为1s 。GPS 临时基

准站的精确坐标选择周边IGS 站以及浙江省的5～6个CORS 站联合解算。文中以庙子湖岛GPS 基准站进行静态大地基准控制，观测到临时验潮站及其他两个验潮站的大地高如表1所示。

表1

点

名

3个验潮站大地高

经度

/（°' ″）…02．69091…54．75103…59．56794

大地高/m19．04919．48218．288

（5）

式中，a 是参考椭球的长半轴；f 是参考椭球的扁－e sin φ；e 为椭球的第率；φ是大地纬度；W =d f =f 0－f ，d a 与d f 分别为椭一偏心率；d a =a 0－a ，

球的长半轴改正和扁率改正；d h 为参考椭球转换引

起的高程变化。

3．基于同步水位观测技术的1985国家高程基准传递

笔者利舟山庙子湖岛缺乏长期的验潮资料，

用两个月时间的短期验潮资料，不足以计算多年平均海面，仅能得到有观测数据期间内的短期平均海面距水尺零点的高度。本文采用同步改正法确定短期站的多年平均海面，如图1所示。同步水位观测的基本假设是长期验潮站与短时验潮站两个站的多年平均海面及短期平均海面高度一致。考虑到长期验潮站的多年平均海面已知，则短期验潮站的多年平均海面从它的水尺零点起算的高度，由式（6）给出

h =h 3+(h 1－h 2)

（6）

h 为多年平均海面与临时验潮站水尺零点之式中，

h 值的大小也可表示为水尺零间的距离，也就是说，点的1985国家高程；h 1为多年平均海面与长期验

潮站水尺零点之间的距离；h 2

为短期平均海面与长期验潮站水尺零点之间的距离；h 3为短期平均海面与临时验潮站水尺零点之间的距离。

纬度/（°' ″）…29．36101…24．87963…35．65791

南码头岱山验潮站沈家门验潮站

2．验潮数据观测

在舟山庙子湖岛南码头获得了2010年8月1日—2010年9月30日两个月的潮位观测数据。分别与岱山、沈家门长期验潮站进行同步观测，采样间隔为1min ，验潮水尺观测不少于1个月（台风等恶劣天气除外），每天3h ，每10min 读取一个数。庙子湖岛1985国家高程水准点与验潮水尺进行三等水准联测，误差不超过±3cm 。

利用同步水位观测的方法，由沈家门长期验潮站推算的庙子湖岛临时验潮站水尺零点的1985国家高程为－3．011m ，而从岱山长期验潮站推算的庙子湖岛临时验潮站水尺零点的1985国家高程为－3．049m 。取两者平均值，临时验潮站水尺零点的1985国家高程为－3．03m 。对庙子湖岛1985国家高程水准点与临时验潮站水尺零点进行三等水准联测，得到1985国家高程水准点与验潮水尺零点距

（下转第32页）

体，水汽总体变化趋势较平稳。水汽含量底层变化平稳，上层变化比较大。结果证明了利用GPS 观测信息层析三维水汽特性的可行性、适用性与先进性，在具有一定密度的GPS 网点信息量的沿海地区，可作为水汽分析、天气预报及气候预测的新手段。

参考文献：

［1］宋淑丽，朱文耀，程宗颐，等．GPS 信号斜路径方向水

．天文学报，2004，45（3）：汽含量的计算方法［J ］338-346．

［2］李国翠，李国平，景华，等．华北三站地基GPS 反演的

．气象科学，2008，28（4）：大气可降水量及其特征［J ］

图10

3测站的综合三维水汽含量

389-394．

［3］毕研盟，毛节泰，刘晓阳，等．应用地基GPS 遥感倾斜

．地球物理学报，2006，路径方向大气水汽总量［J ］49（2）：335-342．

［4］曹玉静，刘晶淼，梁宏，等．基于地基GPS 层析大气水

．自然资源学报，2010，汽资源的方法研究［J ］25（10）：1786-1796．

［5］王勇，何荣，杨彬云，等．GPS 反演的可降水量与降水

J ］．测绘科学，2010，35（5）：80-82．的对比分析研究［

［6］刘敏，郭鹏，叶其欣，等．上海地区地基GPS 水汽三维

．天文学报，2010，51（3）：层析技术和初步应用［J ］299-308．

五、结束语

本文在连云港沿海区域灌云、赣榆、燕尾港3个

JSCORS 站点基础上，以2010年4月1日为例，提出Matlab 工具及层析技术并实施了运用GAMIT 软件、

反演区域三维水汽特征的研究。结果表明，连云港海岸带大气水汽总体变化趋势是底层大气水汽含量较高，随着海拔高度的增加，水汽含量相应降低，而且在同一海拔高度的不同时间，水汽含量变化也不相同。局部地区水汽变化幅度较大，但不会影响总

櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂櫂（上接第4页）得到GPS 点的海拔高以及利用同力大地水准面高，离为5．904m 。从而得到庙子湖岛上1985国家高程水准点（南码头）的高程为2．874m 。

3．数据对比分析由表2可以看出，庙子湖岛南码头的海拔高为2．962m ，而由同步验潮的方法得到南码头水准点的高程为2．874m ，两者相差0．088m ，未超过10cm 。没有前期观测数据的情况下，同步水位观测技术和GPS /重力移去-恢复法均可进行远离大陆海岛礁的高程传递。

表2点

基于移去-恢复法得到3个验潮站的海拔高m 名

大地水准面高16．49814．96315．351

CGCS2000的大地水准面高

16．08714．55214．939

海拔高2．9624．9303．348

步水位观测技术得到GPS 点多年平均海面的1985

国家高程两种技术，在缺乏前期测绘基础信息的情况下，使用两种技术所确定的高程基准，具有一定的参考价值，可为远离大陆的海岛提供公共可行的高程基准值。

参考文献：

［1］李建成．我国现代高程测定关键技术若干问题的研究

．武汉大学学报：信息科学版，2007，及进展［J ］32（11）：980-987．

［2］李建成，姜卫平．长距离跨海高程基准传递方法的研

．武汉大学学报：信息科学版，2001，26（6）：究［J ］514-517．

［3］章传银，郭春喜，陈俊勇，等．EGM2008地球重力场模

．测绘学报，2009，型在中国大陆适用性分析［J ］38（4）：283-289．

［4］张利明．与陆地一致的海岛（礁）高程传递方法：中国，

[1**********]［P ］．2010-06-12．

南码头岱山验潮站沈家门验潮站

四、结束语

利用GPS 大地高与EGM2008重力场模型的重