小区域似大地水准面精化方法的分析-analysis of refining method of quasi - geoid in small area.docxVIP

1绪论1.1研究背景GPS技术的出现给地球科学以及测绘科学注入了新的生命力,实现了三维(几何)坐标定位技术的突破,带来了一场地球科学技术体系的革命,极大地推动了地球科学的发展,随着GPS定位技术的广泛应用，人们已经能够在10-7至10-9的精度量级上，简捷而经济地获得所测点位的平面位置，但是却一直未能以相应的精度解求点的高程(海拔高)。原因是尽管GPS能给出高精度的大地高，却由于没有一个具有相应精度和高分辨率的(似)大地水准面模型，至使在GPS大地高至GPS海拔高的转换中精度严重丢失。为此，目前包括我国在内的国际大地测量学界都在致力于研究区域性高分辨率高精度似大地水准面或大地水准面的建立[1,2,3]。由GPS定位获取的大地高直接解求厘米级精度的正高或正常高已成为一种趋势，其实现可以取代传统的高程测量方法(如三角高程测量和低等级水准测量)，使GPS高程得到广泛应用，进而减轻外业工作量，大大提高工作效率，从而带来直接的经济效益。它将为基础测绘、数字中国地理空间基础框架、区域沉降监测、海洋科学、气象预报、地学研究、交通、水利等多学科研究与应用提供必要的测绘服务，同时精确的似大地水准面模型可以对“数字中国”、“数字城市”等工程的建设起到“加速器”的作用。因此研究精化似大地水准面的理论和方法对我国来说具有特别重要的科学意义以及巨大的社会效益和经济效益。1.2国外似大地水准面的现状和发展20世纪90年代以来，世界各国和各地区先后研制和推出了各自的新一代（似）大地水准面模型，但各个国家和地区的局部大地水准面模型在分辨率和精度方面差异较大。如美国的GEOID99、欧洲的EGG97和加拿大的GSD95都是当地高精度的（似）大地水准面模型。下面就简要介绍写国外似大地水准面的现状和发展情况。整个欧洲地区大地水准面的计算始于20世纪80年初，第一代欧洲重力大地水准面EGG1和EAGG1精度为几分米，分辨率约20km，其中EAGG1比EGG1多了天文重力水准资料，精度略优于后者[4]。从1994年开始，先后推出EGG94、95、96和97序列欧洲重力似大地水准面解(欧洲采用正常高程系统)。这些重力似大地水准面均按移去-恢复方法计算，位模型在EGG94-96中采用GEM-T2、OSU91A以及0SU91A与JGM3的联合，在EGG97中采用EGM96，短波分量采用残差地形模型(RTM)归算法，地面重力数据分辨率优于10km，残差重力异常用快速最小二乘推估法形成1.0′×1.5′的格网数据。相应残差高程异常采用Wenzel的最小二乘频谱组合法计算，由位系数和重力值的误差阶方差定义谱权，其中重力值的零阶阶方差取4mGal2。所得重力似大地水准面与GPS水准确定的似大地水准面的拟合，应用了含趋势项的最小二乘配置法，其中趋势项简单模拟为坐标原点平移量?x、?y、?z对大地水准面高的影响，这一趋势项表现为一个倾斜平面。最新的EGG97以1.0′×1.5′格网表示【5】，分别与德国的GPS网NDS92(300km)，法国的GPS网RBF(1000km)以及横贯欧洲的GPS导线(3000km)进行比较，表明中、长波系统误差为±8.0cm，短波误差信号为±1.3cm。美国在20世纪90年代先后推出了GEOID90、GEOID93和G9501【1】区域大地水准面模型，三个模型计算方法基本相同。由于Rapp(拉普)给出了由OSU91A位系数直接计算高程异常的公式，在G9501的计算中，首先是按莫洛金斯基级数计算高程异常，再转换为大地水准面高。地面重力值的归算采用顾及局部地形改正的布格改正，应用一种“张力连续样条”方法内插形成3′×3′的Bouguer异常格网值，再对这一格网值利用格网平均方法恢复移去的Bouguer片，形成Faye(法耶)异常格网值，其值近似于地形面重力异常加G1项改正，由法耶异常值中移去OSU91A的模型重力异常，对残差重力异常按斯托克斯公式求解，得到残差高程异常，加上OSU91A模型高程异常得最后解，再利用去掉局部地形改正项的布格异常和高程数据，按一阶近似将高程异常转换为大地水准面高。G9501与美国的GPS联邦基准网(FBN)和州建合作基准网(CBN)的GPS水准大地水准面作了比较和拟合，包括1889点。消去系统偏差后的G9501与GPS水准大地水准面之差的经验协方差函数具有L=500km的相关长度和C0=(0.0185)2m2的方差，用此残差由30′×30′格网按最小二乘配置法推估残余误差信号，其中经迭代确定了观测值的噪声方差为6.52cm2，对G9501进行改正得到一个校正大地水准面模型G9501C。对该模型的残差再作一个经验协方差函数，得L=40km，C0=(0.026)2m2。以上表明了两个不同波段的误差源，一个约500km的长波误差源，一个