测绘学报 ›› 2015, Vol. 44 ›› Issue (9): 965-972.doi: 10.11947/j.AGCS.2015.20150110

• 大地测量学与导航 • 上一篇    下一篇

利用GOCE模拟观测反演重力场的Torus法

刘焕玲1,2, 文汉江2, 徐新禹1, 朱广彬3   

  1. 1. 武汉大学测绘学院, 湖北 武汉 430079;
    2. 中国测绘科学研究院地理空间信息工程国家测绘地理信息局重点实验室, 北京 100830;
    3. 国家测绘地理信息局卫星测绘应用中心, 北京 101300
  • 收稿日期:2015-03-03 修回日期:2015-06-09 出版日期:2015-09-24 发布日期:2015-09-24
  • 作者简介:刘焕玲(1986—),女,博士生,主要研究方向为卫星大地测量。E-mail:liuhuanling1986@126.com
  • 基金资助:
    国家973计划(2013CB733302);国家863计划(2013AA122502);国家自然科学基金(41274031;41404014;41204007;41574019);地球空间环境与大地测量教育部重点实验室开放基金(14-02-07);地理空间信息工程国家测绘地理信息局重点实验室经费(777142108)

Torus Approach in Gravity Field Determination from Simulated GOCE Observations

LIU Huanling1,2, WEN Hanjiang2, XU Xinyu1, ZHU Guangbin3   

  1. 1. School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China;
    2. Key Laboratory of Geo-informatics of SBSM, Chinese Academy of Surveying and Mapping, Beijing 100830, China;
    3. Satellite Surveying and Mapping Application Center, National Administration of Surveying, Mapping and Geoinformation, Beijing 101300, China
  • Received:2015-03-03 Revised:2015-06-09 Online:2015-09-24 Published:2015-09-24
  • Supported by:
    The National Basic Research Program of China (973 Program) (No.2013CB733302) The National High-techn Research and Development Program of China (863 Program) (No.2013AA122502) The National Natural Science Foundation of China (Nos.41274031 41404014 41204007 41574019) The Open Fund of Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, China (No.14-02-07) Key Laboratory of Geo-informatics of NASG (No.777142108)

摘要: 在介绍Torus方法反演地球重力场模型的基本原理和方法的基础上,基于圆环面上均匀分布的卫星引力梯度模拟观测值解算了200阶次的地球重力场模型,在无误差情况下,Torus方法解算模型的阶误差RMS小于10-16,验证了该方法的严密性。利用61 d GOCE卫星轨道上无误差的模拟引力梯度观测值解算了200阶次的地球重力场模型,分析了格网化误差、极空白对解算精度的影响,迭代3次后,在不考虑低次系数情况下,模型的大地水准面阶误差和累积误差均较小,最大值仅为0.022 mm和0.099 mm。在沿轨卫星引力梯度模拟数据中加入5 mE/Hz1/2的白噪声,基于Torus方法和空域最小二乘法解算了200阶次的地球重力场模型,Torus方法的精度略低于空域最小二乘法的精度,在不考虑低次项的情况下,两种方法解算模型的大地水准面阶误差最大值分别为1.58 cm和1.45 cm,累积误差最大值分别为6.37 cm和5.55 cm。但由于采用了二维快速傅里叶技术和块对角最小二乘法,极大地提高了计算效率。本文数值结果说明Torus方法是一种独立有效的方法,可用于GOCE任务海量卫星引力梯度观测值反演重力场的快速解算。

关键词: GOCE, 地球重力场, Torus方法, 极空白

Abstract: The principle of Torus method in Earth's gravity field determination from GOCE satellite gradiometry data is discussed. The Earth's gravity field model complete to degree and order 200 is recovered using simulated satellite gradiometry observations on a Torus grid, and the degree error RMS is smaller than 10-16, which shows the effectiveness of Torus approach. The gravity field model is also resolved using the simulated satellite gradiometry observations given on GOCE orbits of 61 days. The influences of interpolation and polar gaps are analyzed. Without considering the low-order coefficients the geoid degree errors and cumulative errors are very small after three iterations. The maximums of them are only 0.022 mm and 0.099 mm. The white noise with PSD 5 mE/Hz1/2 is added to the simulated observations and the gravity field model complete to degree and order 200 is also computed. The model is compared with that model which is derived using space-wise LS method and the same observations. It shows that the precision of Torus is slightly lower. Without considering the low-order coefficients the maximum geoid degree errors of Torus and space-wise LS method are 1.58 cm and 1.45 cm, and the maximum cumulative geoid errors are 6.37 cm and 5.55 cm, respectively. But the computational efficiency of Torus is greatly improved by using the two-dimensional FFT and the block-diagonal least-squares adjustment. The numerical results show that Torus method is independent and valid. Meanwhile fast resolution of gravity field based on massive amount of GOCE satellite gradiometry observations is feasible.

Key words: GOCE, earth's gravity field model, Torus method, polar gaps

中图分类号: