广州市重力异常计算及其在似大地水准面精化中的应用与探讨

广州市重力异常计算及其在似大地水准面精化中的应用与探讨

摘要：地区似大地水准面精化工作中一个重要环节即是测定并计算该地区的重力异常数据，本文以广州市为例，探讨该地区重力异常的计算方法，及其在采用移去-恢复法进行似大地水准面精化工作中的应用，并对其中若干问题展开探讨。

关键词：重力异常似大地水准面移去-恢复法

1引言

一般有以下几何方法可以确定四大地水准面：天文水准、卫星测高、GPS水准，还有重力法，和将几何与重力联合法。现在一般使用联合法精化局部地区四大地水准面。这种方法采用的是现将通过GPS 水准的得到的高精度但分辨率较低的几何大地水准面作为控制，再拟合通过重力学法确定的高分辨但精度低的重力大地水准面，这样来精化局部似大地水准面。[1]

地区似大地水准面精化工作中一个重要环节即是测定并计算该地区的重力异常数据，该数据的精确性和可靠性直接影响到似大地水准面的精度。2005年，“广州市高精度三维控制网的建立及似大地水准面确定”项目利用全球高阶地球重力场模型、实测重力数据和GPS 水准数据及高分辨率数字地形模型（DTM），确定了广州市1km和2km分辨率的空间重力异常及厘米级似大地水准面。计算范围为：在WGS-84格网坐标下，纬度从22度33分9.550秒至23度59分

49.549747秒，经度从112度54分18.893701秒至114度6分38.894秒，南北方向160km，东西方向124km。为此，本文依托以该项目的实际数据，以广州市为例，探讨该地区重力异常的计算方法，对其数据在广州市似大地水准面精化工作中的应用展开分析，并对其中若干问题进行探讨。

2重力场模型的选取

选择本地区的高阶全球重力场模型作为参考重力场，在实际计算局部大地水准面的时候是非常必要的。到现在为止同阶次模型（360阶）EGM96被公认为是最好的，原武汉测绘科技大学自行研制的WDM94是360阶全球重力场模型，所以，对这两个重力场模型进行实测高精度重力和GPS水准数据比较和评价，用来作为广州计算高精度大地水准面选取最适合本地区的参考模型。

利用GPS测定的椭球高和精密水准测量测定的正常高可得到实测似大地水准面高(或高程异常)为：

由地球重力场模型可按下式计算任意一点的大地水准面高(Heiskanen and Moritz, 1967)，

这里GM 为地心引力常数；为计算点的正常重力；a是参考椭

球的长半径；、和r分别是计算点的地心纬度、经度和向径；和为完全规格化位系数；是完全规格化缔合Legendre函数；nmax是计算模型的最大阶数，在本项目的计算中，EGM96和WDM94均取为360。

以WGS-84为参考椭球，则实测重力观测值可归算为大地水准面上的空间重力异常，实际上为我国1956黄海高程基准面上的空间重力异常，即

其中为实测绝对重力值，为正常重力，为空间改正，并且

这里和分别为计算点的大地纬度和正常高。

利用地球重力场模型计算空间重力异常的公式为：

式中各量的含义同（2）式。

因此，由（3）和（6）两式可计算它们的差值为：

3实测值与模型计算值之间的比较分析

3.1源数据选取

采用的重力数据包括：收集到的广州地区及周边3838个陆地重力数据，其密度为每6平方公里一个点。

GPS水准数据：在广州建立了由123个点组成的GPS二等水准控制网，其平均间距约14公里，并按二等水准测量的要求与原有一、二等水准点进行了联测，正常高为1985国家高程基准。

数字地面模型（DTM）数据：高分辨率的数字地形模型是计算高分辨率高精度大地水准面的重要信息，为此，利用广东省国土厅信息中心提供的100m分辨率的DTM，分别建立了广州及其周边地区分辨率为100m、500m、1km和2km的DTM，覆盖范围为：在WGS-84格网坐标下，纬度从22度25分51.6547秒至24度2分45.9047秒，经度从112度49分17.1189秒至114度11分28.6192秒，南北方向179km，东西方向141km。

3.2实测大地水准面与模型计算值之间的比较

表1列出了利用123个GPS水准数据计算的值统计结果。从表1可以看出，利用EGM96、WDM94计算大地水准面高的精度（标准差）分别为?0.1271m和?0.2119m，而大地水准面高差的绝对精度则分别为?0.154m和?0.258m。

一般说来，尽可能地消除系统偏差后的模型大地水准面与GPS

水准的符合精度反映了重力场模型的实际精度。因此，首先采用五参数模型消除用模型计算的大地水准面中的系统偏差，即

其中为未知参数，和分别为大地纬度和经度，为随机噪声。

采用最小二乘法解算模型（8）中的未知参数，并且根据模型显著性检验的结果来决定采用三参数、五参数或五参数模型。消除系统偏差后，由地球重力场模型计算的大地水准面高和大地水准面高差的统计结果分别列于表1。由表1可知，消除系统偏差后模型大地水准面的精度得到显著提高，用EGM96和WDM94计算大地水准面高的绝对精度（标准差）分别为? 0.0486m、?0.0438m，而大地水准面高差的绝对精度则分别为?0.069m、?0.063m，这说明在利用EGM96和WDM94来表示广州局部重力大地水准面非常接近。

表1 由地球重力场模型计算的大地水准面高与GPS水准测定的大地水准面高之比较结果

注：before是指没有消除系统偏差的结果，而after则指消除系统偏差后的结果，以下同。

3.3 实测重力异常数据与模型计算值之间的比较

表2给出了的统计结果。从该表可以看出，由EGM96和WDM94计算空间重力异常的精度（标准差）分别为?9.8674mGal、?8.9819mGal。这些结果表明，WDM94稍优于EGM96重力场模型。

综上所述，利用地球重力场模型EGM96和WDM94来表示广州局部重力场相当接近，此次计算我们选取EGM96作为计算广州似大地水准面的参考重力场模型，以下计算都是基于EGM96模型进行的。

表 2 由实测重力异常与地球重力场模型计算的重力异常的比较结果（单位：mGal）

4格网空间重力异常的计算

4.1 计算方法

本文计算格网空间重力异常使用的是移去-恢复方法，也就是先利用高阶地球重力场模型计算出中场波重力异常，在离散重力点的重力异常中将其消除，通过DTM计算的地形改正（短波重力异常），这样得到离散重力点的残差重力异常，再进行残差重力异常的拟合，之后得到网格残差重力异常，最后在格网残差重力异常中恢复重力场模型和DTM的贡献，最终获得格网空间重力异常。

离散重力点的残差重力异常为：