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

广州市重力异常计算及其在似大地水准面精化中应用与探讨 　　摘 要：地区似大地水准面精化工作中一个重要环节即是测定并计算该地区的重力异常数据，本文以广州市为例，探讨该地区重力异常的计算方法，及其在采用移去-恢复法进行似大地水准面精化工作中的应用，并对其中若干问题展开探讨。 　　关键词：重力异常 似大地水准面 移去-恢复法 　　中图分类号：P22 文献标识码：A 文章编号：1672-3791（2013）04（c）-0042-02 　　一般有以下几何方法可以确定四大地水准面：天文水准、卫星测高、GPS水准，还有重力法，和将几何与重力联合法。现在一般使用联合法精化局部地区四大地水准面。这种方法采用的是现将通过GPS水准的得到的高精度但分辨率较低的几何大地水准面作为控制，再拟合通过重力学法确定的高分辨但精度低的重力大地水准面，这样来精化局部似大地水准面[1]。 　　地区似大地水准面精化工作中一个重要环节即是测定并计算该地区的重力异常数据，该数据的精确性和可靠性直接影响到似大地水准面的精度。2005年，“广州市高精度三维控制网的建立及似大地水准面确定”项目利用全球高阶地球重力场模型、实测重力数据和GPS水准数据及高分辨率数字地形模型（DTM），确定了广州市1 km和2 km分辨率的空间重力异常及厘米级似大地水准面。计算范围为：在WGS-84格网坐标下，纬度从22度33分9.550秒至23度59分49.549747秒，经度从112度54分18.893701秒至114度6分38.894秒，南北方向160 km，东西方向124 km。为此，本文依托以该项目的实际数据，以广州市为例，探讨该地区重力异常的计算方法，对其数据在广州市似大地水准面精化工作中的应用展开分析，并对其中若干问题进行探讨。 　　1 重力场模型的选取 　　选择本地区的高阶全球重力场模型作为参考重力场，在实际计算局部大地水准面的时候是非常必要的。到现在为止同阶次模型（360阶）EGM96被公认为是最好的，原武汉测绘科技大学自行研制的WDM94是360阶全球重力场模型，所以，对这两个重力场模型进行实测高精度重力和GPS水准数据比较和评价，用来作为广州计算高精度大地水准面选取最适合本地区的参考模型。 　　利用GPS测定的椭球高和精??水准测量测定的正常高可得到实测似大地水准面高（或高程异常）为： 　　这里GM为地心引力常数；为计算点的正常重力；a为参考椭球的长半径；、和r分别为计算点的地心纬度、经度和向径；和为完全规格化位系数；为完全规格化缔合Legendre函数；nmax为计算模型的最大阶数，在本项目的计算中，EGM96和WDM94均取为360。 　　以WGS-84为参考椭球，则实测重力观测值可归算为大地水准面上的空间重力异常，实际上为我国1956黄海高程基准面上的空间重力异常，即： 　　2 实测值与模型计算值之间的比较分析 　　2.1 源数据选取 　　采用的重力数据包括：收集到的广州地区及周边3838个陆地重力数据，其密度为每6平方公里一个点。 　　GPS水准数据：在广州建立了由123个点组成的GPS二等水准控制网，其平均间距约14 km，并按二等水准测量的要求与原有一、二等水准点进行了联测，正常高为1985国家高程基准。 　　数字地面模型（DTM）数据：高分辨率的数字地形模型是计算高分辨率高精度大地水准面的重要信息，为此，利用广东省国土厅信息中心提供的100 m分辨率的DTM，分别建立了广州及其周边地区分辨率为100 m、500 m、1 km和2 km的DTM，覆盖范围为：在WGS-84格网坐标下，纬度从22度25分51.6547秒至24度2分45.9047秒，经度从112度49分17.1189秒至114度11分28.6192秒，南北方向179 km，东西方向141 km。 　　2.2 实测大地水准面与模型计算值之间的比较 　　表1列出了利用123个GPS水准数据计算的值统计结果。从表1可以看出，利用EGM96、WDM94计算大地水准面高的精度（标准差）分别为±0.1271 m和±0.2119 m，而大地水准面高差的绝对精度则分别为±0.154 m和0.258 m。 　　一般说来，尽可能地消除系统偏差后的模型大地水准面与GPS水准的符合精度反映了重力场模型的实际精度。因此，首先采用五参数模型消除用模型计算的大地水准面中的系统偏差，即： 　　其中为未知参数；和分别为大地纬度和经度；为随机噪声。 　　采用最小二乘法解算模型（8）中的未知参数，并且根据模型显著性检验的结果来决定采用三参数、五参数或五参数模型。消除系统偏差后，由地球重力场模型计算的大地水准面高和大地水准面高差的统计结果分别列于表1。由表1可知，消除系统偏差后模型大地水准面的精度得