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

doi:10.11947/j.AGCS.2015.20150110

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

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

刘焕玲^1,2, 文汉江², 徐新禹¹, 朱广彬³

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 Huanling^1,2, WEN Hanjiang², XU Xinyu¹, ZHU Guangbin³

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)

摘要/Abstract

摘要： 在介绍Torus方法反演地球重力场模型的基本原理和方法的基础上,基于圆环面上均匀分布的卫星引力梯度模拟观测值解算了200阶次的地球重力场模型,在无误差情况下,Torus方法解算模型的阶误差RMS小于10^-16,验证了该方法的严密性。利用61 d GOCE卫星轨道上无误差的模拟引力梯度观测值解算了200阶次的地球重力场模型,分析了格网化误差、极空白对解算精度的影响,迭代3次后,在不考虑低次系数情况下,模型的大地水准面阶误差和累积误差均较小,最大值仅为0.022 mm和0.099 mm。在沿轨卫星引力梯度模拟数据中加入5 mE/Hz^1/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/Hz^1/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

中图分类号:

P223

刘焕玲, 文汉江, 徐新禹, 朱广彬. 利用GOCE模拟观测反演重力场的Torus法[J]. 测绘学报, 2015, 44(9): 965-972.

LIU Huanling, WEN Hanjiang, XU Xinyu, ZHU Guangbin. Torus Approach in Gravity Field Determination from Simulated GOCE Observations[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(9): 965-972.

参考文献

[1] ESA. Gravity Field and Steady-state Ocean Circulation Mission. Reports for Mission Selection of the Four Candidate Earth Explorer Core Missions[R]. ESA Publications Division, ES SP-1233(1), 1999.
[2] DRINKWATER M R, HAAGMANS R, MUZI D, et al. The GOCE Gravity Mission: ESA's First Core Earth Explorer[R]. Proceedings of the 3rd International GOCE User Workshop, Frascati, Italy, ESA Special Publication, SP-627, ISBN 92-9092-938-3, 2006: 1-8.
[3] RUMMEL R, VAN GELDEREN M, KOOP R, et al. Spherical Harmonic Analysis of Satellite Gradiometry[M]. Netherlands Geodetic Commission: Publications on Geodesy, New Series 39, 1993.
[4] KOOP R. Global Gravity Field Modeling Using Satellite Gravity Gradiometry[M]. Netherlands Geodetic Commission: Publications on Geodesy, New Series 38, 1993.
[5] PAIL R, BRUINSMA S, MIGLIACCIO F, et al. First GOCE Gravity Field Models Derived by Three Different Approaches[J]. Journal of Geodesy, 2011, 85(11): 819-843.
[6] SNEEUW N J. A Semi-analytical Approach to Gravity Field Analysis from Satellite Observations[D]. Munich, Germany: Institut für Astronomische und Physikalische Geodsie, Technische Universitt München, 2000.
[7] SNEEUW N J. Space-wise, Time-wise, Torus and Rosborough Representations in Gravity Field Modeling[J]. Space Science Reviews, 2003, 108(1-2): 37-46.
[8] KLEES R, DITMAR P. The Performance of the Time-wise Semi-analytical Inversion of Satellite Gravity Gradients[M]//DM J, SCHWARZ K P. Vistas for Geodesy in the New Millennium, International Association of Geodesy Symposia. Berlin Heidelberg: Springer, 2001, 125: 253-258.
[9] SCHUH W D, PAIL R, PLANK G. Assessment of Different Numerical Solution Strategies for Gravity Field Recovery[C]// Proceedings of the 1st International GOCE User Workshop, ESA WPP-188, 87-95, ESA/ESTEC, 2001.
[10] PAIL R, WERMUTH M. GOCE SGG and SST Quick-look Gravity Field Analysis[J]. Advances in Geosciences, 2003, 1: 5-9.
[11] PAIL R, PLANK G. GOCE Gravity Field Processing Strategy[J]. Studia Geophysica et Geodaetica, 2004, 48(2): 289-309.
[12] XU Chen. The Torus-based Semi-analytical Approach in Spaceborne Gravimetry[D]. Calgary: Department of Geomatics Engineering, University of Calgary, 2008.
[13] KAULA W M. Theory of Satellite Geodesy: Applications of Satellites to Geodesy[M]. Waltham Massachusetts: Blaisdell Publishing Company, 1966.
[14] SCHRAMA E J O. The Role of Orbit Errors in Processing of Satellite Altimeter Data[M]. Netherlands Geodetic Commission: Publications on Geodesy, New Series 33, 1989.
[15] PAIL R, PLANK G. Assessment of Three Numerical Solution Strategies for Gravity Field Recovery from GOCE Satellite Gravity Gradiometry Implemented on a Parallel Platform[J]. Journal of Geodesy, 2002, 76(8): 462-474.
[16] ESA. GOCE HPF: GOCE Level 2 Product Data Handbook[R]. Technical Note, GO-MA-HPF-GS-0110, 2010.
[17] CAPITAINE N, WALLACE P T, MCCARTHY D D. Expressions to Implement the IAU 2000 Definition of UT1[J]. Astronomy & Astrophysics, 2003, 406(3): 1135-1149.
[18] PAVLIS N K, HOLMES S A, KENYON S C, et al. An Earth Gravitational Model to Degree 2160: EGM2008[J]. EGU General Assembly, 2008, 10: 13-18.
[19] ESA. GOCE L1b Products User Handbook[R]. Technical Note, GOCE-GSEG-EOPGTN-06-0137, 2006.
[20] SNEEUW N J, VAN GELDEREN M. The Polar Gap[M]//SANSÓ F, RUMMEL R. Geodetic Boundary Value Problems in View of the One Centimeter Geoid. Lecture Notes in Earth Sciences. Berlin Heidelberg: Springer, 1997, 65: 559-568.
[21] XU Xinyu, LI Jiancheng, JIANG Weiping, et al. Simulation Study for Recovering GOCE Satellite Gravity Model Based on Space-wise LS Method[J]. Acta Geodaetica et Cartographica Sinica, 2011, 40(6): 697-702. (徐新禹, 李建成, 姜卫平, 等. 基于空域最小二乘法求解 GOCE 卫星重力场的模拟研究[J]. 测绘学报, 2011, 40(6): 697-702.)
[22] RUDOLPH S, KUSCHE J, IIK K H. Investigations on the Polar Gap Problem in ESA's Gravity Field and Steady-state Ocean Circulation Explorer Mission (GOCE)[J]. Journal of Geodynamics, 2002, 33(1-2): 65-74.
[23] XU Xinyu, LI Jiancheng, WANG Zhengtao, et al. The Simulation Research on the Tikhonov Regularization Applied in Gravity Field Determination of GOCE Satellite Mission[J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(5): 465-470. (徐新禹, 李建成, 王正涛, 等. Tikhonov正则化方法在GOCE重力场求解中的模拟研究[J]. 测绘学报, 2010, 39(5): 465-470.)
[24] ZHU Guangbin, LI Jiancheng, WEN Hanjiang, et al. Slepian Localized Spectral Analysis of the Determination of the Earth's Gravity Field Using Satellite Gravity Gradiometry Data[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(1): 1-7. (朱广彬, 李建成, 文汉江, 等. 卫星重力梯度数据确定地球重力场的Slepian局部谱分析方法[J]. 测绘学报, 2012, 41(1): 1-7.)
[25] BAUR O, SNEEUW N, GRAFAREND E W. Methodology and Use of Tensor Invariants for Satellite Gravity Gradiometry[J]. Journal of Geodesy, 2008, 82(4-5): 279-293.

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

Torus Approach in Gravity Field Determination from Simulated GOCE Observations

PDF

可视化

被引次数

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[1]	潘娟霞, 邹贤才. GOCE引力梯度内部校准方法[J]. 测绘学报, 2022, 51(2): 192-200.
[2]	陈鑑华, 张兴福, 沈云中, 陈秋杰, 李伟超. GOCE卫星重力梯度观测值的时变重力场变化改正及影响分析[J]. 测绘学报, 2021, 50(3): 324-332.
[3]	刘焕玲, 文汉江, 徐新禹, 赵永奇, 蔡剑青. GOCE实测数据反演高阶重力场模型的Torus方法[J]. 测绘学报, 2020, 49(8): 965-973.
[4]	黄炎, 王庆宾, 冯进凯, 邢志斌, 范雕, 谭勖立, 吕明昊. 局部地形改正快速计算的GPU并行的棱柱法[J]. 测绘学报, 2020, 49(11): 1430-1437.
[5]	瞿庆亮, 常晓涛, 于胜文, 朱广彬. 利用先验重力场模型求定GOCE卫星重力梯度仪校准参数[J]. 测绘学报, 2019, 48(2): 176-184.
[6]	邹贤才. GOCE卫星重力梯度校准与无阻尼控制效果分析[J]. 测绘学报, 2018, 47(3): 291-297.
[7]	田家磊, 李新星, 吴晓平, 邢志斌. 超高阶重力场模型最小二乘快速实现[J]. 测绘学报, 2018, 47(11): 1437-1445.
[8]	祝竺, 赵艳彬, 廖鹤, 涂海波, 张国万, 魏小刚. 星载原子干涉技术用于地球重力场测量及其精度评估[J]. 测绘学报, 2017, 46(9): 1088-1097.
[9]	游为. 利用局部高空间分辨率大气数据计算GRACE大气去混频模型[J]. 测绘学报, 2017, 46(3): 316-324.
[10]	孟祥超, 万晓云, 于锦海, 朱永超, 冯炜. 利用引力梯度数据计算径向梯度的优化方法[J]. 测绘学报, 2016, 45(7): 775-781.
[11]	周浩, 罗志才, 钟波, 陆飚. 利用最小二乘直接法反演卫星重力场模型的MPI并行算法[J]. 测绘学报, 2015, 44(8): 833-839.
[12]	叶周润, 柳林涛, 梁星辉. Moho面扰动重力梯度信息的提取[J]. 测绘学报, 2015, 44(6): 609-615.
[13]	苏勇, 范东明, 游为. 利用GOCE卫星轨道数据恢复地球重力场模型方法的分析[J]. 测绘学报, 2015, 44(2): 142-149.
[14]	徐新禹, 赵永奇, 魏辉, 吴汤婷. 利用先验重力场模型对GOCE卫星重力梯度观测值进行校准分析[J]. 测绘学报, 2015, 44(11): 1196-1201.
[15]	李新星吴晓平李姗姗刘晓刚崔志伟. 块对角最小二乘方法在确定全球重力场模型中的应用[J]. 测绘学报, 2014, 43(8): 778-785.