2005年珠峰测高中的平面测量数据的获取与处理

2005年珠峰测高中的平面测量数据的获取与处理2005年珠峰测高中的平面测量数据主要通过珠峰地区3个网进行测定。这3个网分别是GPS监测网、GPS控制网和GPS联测网。GPS监测网以青藏周边地区的17个GPS连续运行站为框架点，GPS控制网以GPS监测网点为基准，ePS联测网由包括珠峰峰顶1个GPS点和地面交会珠峰的7个GPS点构成(其中定日点系珠峰地区临时性GPS连续运行站)。珠峰地区3个GPS网的数据进行处理时，GPS联测网则以GPS控制网为基准。所有这些GPS网的数据进行处理时，采用IGS精密星历、WGS一84坐标系统和GLOBKl0．2版软件进行整体平差；坐标框架采用ITRF2000，历元采用平均瞬时观测历元。2005年珠峰测高的高程控制数据的获取与处理分为2个部分，一是从相距珠峰约160km的国家一等水准点“I萨拉40基”开始的高程传递，用二等、三等水准将我国黄海高程基准通过该点传递至测定珠峰高程的地面测站。各种等级的水准路线长度累计约400km。此外，在水准无法工作的地区，还施测了测距高程导线以传递高程，总长20．5km。2005年珠峰测高中峰顶的高程必须由珠峰地区的大地水准面起算。因此，2005年珠峰测高的高程控制数据的获取与处理的第二个部分工作是精化该地区的局域大地水准面。这项工作的基础是推算具有较高分辨率的栅格重力值，由于珠峰地区地形崎岖，重力资料稀缺，因此，在归算和推估珠峰地区重力时，必须充分利用地形资料和均衡模型。经过试算比较，决定在2005珠峰重力推估中，采用爱黎一海斯卡宁均衡模型和34km的均衡低偿深度。在利用栅格数字地形模型(DTM)方面，国内部分主要采用国家测绘局珠峰地区的1”×l”DTM，国外部分以SRTM3的DSM为基础，结合GToP030(3∥×3矿)的DTM。最后计算了珠峰地区两种栅格(3矿×30”和2．5’×2．5’)形式的重力空间异常。利用这一栅格重力资料，采用移去恢复法，精化与EGM96相应的珠峰地区的重力大地水准面。为了将这一精化后的珠峰地区的重力大地水准面定位于我国黄海高程基准，由此选用珠峰地区的44个GPS／水准点与该重力大地水准面进行拟合，作为在2005年珠峰测高计算中最终采用的似大地水准面。经过外部检测其精度为±0．09m。此次测量主要采用GPs技术和经典大地测量技术综合进行珠峰高程测量。因此在珠峰峰顶上树立了测量觇标，以供三角高程测量精确照准之用之外，在该觇标上还架设了GPs天线和激光测距的反射棱镜，供GPS测量和激光测距仪照准之用。经典大地测量技术测定珠峰峰顶(即峰顶测量觇标处，以下同)雪面大地高时，采用单向三角高程结合激光测距这二种技术进行，与GPS技术测得的珠峰峰顶雪面大地高的值取权平均值，最终求得珠峰顶雪面大地高。顾及该处的大地水准面，由此求得珠峰顶雪面的正高(海拔高)为8847．93m，该处相应的冰雪层厚度为3．50m。根据上面的计算，2005年珠峰测高中珠峰峰顶岩面正高(海拔高)的国家采用值为8844．43m，精度为±0．21m。参考文献[1]常吉庆.历次珠峰高程测量情况简介[J].测绘通报,2005,(10):2-6.[2]张燕平,张鹏,岳建利等.2005珠穆朗玛峰高程测量[J].地理信息世界,2005,3(4):48-52.[3]郭春喜,王斌,程传录等.珠穆朗玛峰高程测量[J].地球科学与环境学报,2009,31(1):106-110.2005年的珠峰高程复测是一项GPS、精密水准测量、大地水准面精化、重力测量、雷达探测、气象探测以及地学理论等大地测量各种相关技术的国家基础测绘项目，经数据处理取得了重要的可靠成果。以上是我国军民大地测量工作者和登山队员艰辛劳动的结晶35年来所获得的历次珠峰高程测量的成果，在此谨向他们致以崇高的敬意和谢意。珠峰的测量分两个阶段第一阶段：海拔5600米之前———水准测量法“从拉孜(位于西藏自治区西南部,为此次测量起点)到5600米的珠峰半山坡都要使用水准测量法测量高度，但使用这种传统而精确方法测量每天只能行进4公里。”总体来说，测量珠峰行动分为两大部分，一部分是从拉孜出发行进500公里到达珠峰5600米处；另一部分从5600米起使用觇标、GPS、雪深雷达综合测量峰顶高度。”据介绍，所谓水准测量法是指从最初的基准点———青岛海拔基本面开始，每35米设立一个标杆，用水准测量仪计算两边标杆尺度的变化，在水平视线内测出两根标尺之间的高差，然后通过这种在两个相距70米标杆间设立水准仪的测量方法一站一站地将高差累加起来，直到珠峰5600米处。“这种方法能够把误差控制在每公里0.5毫米以内，是目前世界上最精确的方法。不过好在拉孜县的海拔早在1997年就已经测定，此次测量珠峰我们只需从拉孜出发行程500公里就能到达珠峰5600米处的下一个测绘站了。”第二阶段：海拔5600米以后———6点联测确保精度在用传统水准测量法到达5600米高程之后，测量工作将转入第二阶段，测量人员将直接进行珠峰山体测量。这一阶段的测量将是由测量人员在观测点通过观测登山队员立到珠峰顶上的觇标，通过计算最终得出珠峰山体高度。为了提高测量精度，本次珠峰测量一共在珠峰脚下部下了6个观测点，届时观测队员将进行6点联测。据悉这种多角度测量是测量精度的可靠保证。1975年我国进行第一次大规模珠峰测量时，觇标只具有角度测量功能，而这次设计的觇标添加了专门的棱镜，增加了测距功能。这也进一步提高了测量的精度。针对珠峰顶部冰雪层厚度年际变化大的问题，此次珠峰测量还首次动用了冰雪深雷达。据悉，这将帮助测量人员搞清珠峰顶部冰雪层与岩石层之间的关系。这将帮助测量人员掌握珠峰顶岩石层的高度。这个高度将不会随着冰雪层的变化而变化。据透露，此次测量，登山队员还将在珠峰8300米处的一块坚固岩石上竖立一根永久性的觇标。测量人员将每年定期观测这根觇标，以掌握珠峰山体的年际变化。观测数据将成为“珠峰究竟仍在长高还是正在坍塌”问题的最有力证据。而由于珠峰地区地质构造极为复杂，珠峰山体的重力线并不是一根直线。此次测量，科研人员还大量进行了不同地点的重力测量，这些工作将描述出珠峰山体重力线弯曲的具体情形，并为最终修正珠峰高度数据提供科学依据。而觇标是放置于峰顶的一种测量标志。此次珠峰测量，科研人员将在珠峰下的观测点通过观测设在珠峰顶部的觇标最终计算出珠峰高度。而这根将矗立在世界最高峰顶上的觇标是地道的北京造。1)GPS监测网(如图1所示)自青海格尔木至珠峰地区，跨越冈底斯、喜马拉雅构造带布设GPs监测网，并联测其万方数据他地壳运动GPS监测网点，主要目的是对青藏板块运动进行监测和研究，同时将珠峰高程测量的控制网纳入到国家整体坐标基准框架和地壳运动整体监测网中。2)珠峰GPs控制网(如图2所示)在珠峰及邻近地区布设，在选点和观测模式等方面区别于GPS监测网。作为本次珠峰地区的控制网，控制点沿水准路线布设，联系GPs监测网，并形成GPS／水准点，辅以加密重力测量等手段精确求定珠峰高程测量的起算点并服务于大地水准面精化。3)水准测量(如图3所示)由最新一期国家一等水准网中的I萨拉40基点通过二等和三等水准测量、高程导线测量将国家高程基准准确地传递至珠峰高程测量起算点。。4)；：；重力观测(如霆碍所零》i—i：；|¨』些四图22005年珠穆朗玛峰GPs控制网展点图(共32点)用于局部重力场优化、垂线偏差计算、大地水准面精化以及水准测量改正计算等。5)峰顶测量(如图5所示)主要包括GPS测量以及雷达探测确定峰顶高度及冰雪层厚度、采用测距和三角高程交会测量方法确定峰顶的位置和高度等。峰局部重力场精化3．1基础数据整理与处理珠峰地区局部重力场精化[1—3所采用的基础资料，主要包括以下几个方面：(1)以珠峰地区SRTM3(航天雷达测绘任务获得的3”地形数据)为基础，结合1：50000DEM和全球GTOP030数据，综合确定珠峰地区分辨率为3”X3”的数字地形模型。(2)收集与采用了2501个重力点成果，包括中(国)尼(泊尔)重力剖面成果和中国在珠峰地区施测的重力成果。(3)采用重力点成果、3”×3”DEM数据，完成了30”×30”和2．5’×2．51平均空间异常的计算。(4)采用44个GPS水准点作为拟合珠峰地区重力似大地水准面的重控制数据。3．2珠峰地区均衡模型及均衡深度确定由于珠峰地区重力资料较少，因此在重力归算和推估中应采用合适的均衡模型和抵偿均衡深度。1975年计算珠峰地区均衡改正时，因受当时计算条件等限制，采用的是海福特模板，均衡抵偿深度采用113．7km。本次计算珠峰地区均衡改正时，拟采用目前运用较为广泛的爱黎模型，采用快速算法进行积分。为此采用普拉特模型和爱黎模型分别进行了试算，并对两种模型的计算结果进行了比较，同时通过试算确定了爱黎模型在珠峰地区应采用的均衡深度。在两种均衡模型的试算中，主要采用登山路线上5个重力点的地形均衡异常的计算结果进行比较分析。采用普拉特均衡模型计算时，均衡低偿深度采用113．7km，积分半径采用166．7km，对于爱黎模型，选择了34、35、37、40km4个均衡深度分别计算了登山路线上5个重力点的地形均衡异常。从两种均衡模型的几种试算结果可以看出，均衡模型选用爱黎模型，均衡深度采用34km时，登山路线上5个重力点的均衡异常的数值小且变化均匀，因此在计算珠峰地区均衡改正时，均衡模型选择爱黎模型，均衡深度采用34km。3．3珠峰地区似大地水准面的确定利用珠峰地区加密重力点成果、高分辨率数字地形模型、高阶次的地球重力场模型及分布较均匀、现势性较好的GPS水准成果，采用重力法(Molodell—sky原理)及移去一恢复技术完成了珠峰地区分辨率为2．5’x2．5’的高精度似大地水准面精化工作。计算时选择EGM96、WDM94、IGC如5B、DQM2000D和CG03C5种参考重力场模型进行了试算比较，确定了最终采用的似大地水准面模型。(1)选择E∞他6、WDM94、IGG05B、DQM2000D和CG03C作为重力场模型，分别完成了珠峰地区模型似大地水准面和模型平均空间异常计算。(2)由实测的平均空间异常、模型平均空间异常以及前面计算的地形改正，计算剩余法耶异常，用剩余法耶异常计算区域重力似大地水准面。(3)在采用重力法计算区域重力似大地水准面时，考虑到地球重力场模型的阶次和积分半径大小对似大地水准面的计算精度有着重要影响，因此，采用Molodensky公式，每个地球重力场模型选用360阶次，按5km间隔，选择从20km到100km的积分半径进行区域重力似大地水准面的试算。经比较分析，采用55km积分半径，选择EGM96，wDM94、IGG05B、DQM2000D和CG03C作为参考重力场模型，分别完成了珠峰地区区域重力似大地水准面的计算。(4)区域重力似大地水准面采用的是GRS80椭球，而GPS水准网则使用的是WGS84椭球，加上重力基准等因素的影响，使得GPS水准与重力似大地水准面存在一定的差异，对此利用44个GPS水准点，采用二元二次多项式拟合法，分别将EGM96、WDM94、IGG05B、DQM2000D和CG03C作为参考重力场模型计算的区域重力似大地水准面拟合适配于该地区GPS水准网，由此获得了相应的拟合似大地水准面成果及有关的精度信息。(5)利用拟合似大地水准面的格网成果，通过内插方法求得各GPS水准点的内插高程异常值，然后将内插高程异常值同GPS水准点的已知高程异常值进行比较，从而计算拟合似大地水准面的高程异常残差中误差。从计算结果发现，分别选择EGM96、WDM94、IGG05B、DQM2000D和IG03C作为参考重力场模型确定的珠峰地区拟合似大地水准面的精度比较接近，考虑到EGM96地球重力场模型在全球的广泛应用，并结合该模型在中国的应用实践，选择EGM96作为参考重力场模型确定的珠峰地区拟合似大地水准面作为最终似大地水准面。利用最终确定的似大地水准面，根据珠峰峰顶坐标，采用Shepard插值法完成了珠峰顶高程异常等：珠穆朗玛峰高程测量109值的推算，求得