GOCE实测数据反演高阶重力场模型的Torus方法

doi:10.11947/j.AGCS.2020.20200044

测绘学报 ›› 2020, Vol. 49 ›› Issue (8): 965-973.doi: 10.11947/j.AGCS.2020.20200044

GOCE实测数据反演高阶重力场模型的Torus方法

刘焕玲^1,2,3, 文汉江^1,2, 徐新禹⁴, 赵永奇⁴, 蔡剑青⁵

1. 中国测绘科学研究院, 北京 100036;
2. 自然资源部测绘科学与地球空间信息技术重点实验室, 北京 100036;
3. 城市空间信息工程北京市重点实验室, 北京 100038;
4. 武汉大学测绘学院, 湖北武汉 430079;
5. 德国斯图加特大学大地测量研究所, 斯图加特 70174

收稿日期:2020-02-13 修回日期:2020-05-29 发布日期:2020-08-25
通讯作者: 徐新禹 E-mail:xyxu@sgg.whu.edu.cn
作者简介:刘焕玲(1986-),女,助理研究员,主要研究方向为卫星大地测量、物理大地测量。E-mail:liuhl@casm.ac.cn
基金资助:
高分辨率对地观测系统重大专项（42-Y20A09-9001-17/18）；国家自然科学基金（41574019；41774020）；DAAD战略合作伙伴与专题研究合作网络项目（57421148）；中国测绘科学研究院基本科研业务费（AR1922）；城市空间信息工程北京市重点实验室经费资助项目（2019211）；武汉大学地球空间环境与大地测量教育部重点实验室开放基金（18-01-05）

Torus approach in high degree gravity field model determination from GOCE satellite gradiometry observations

LIU Huanling^1,2,3, WEN Hanjiang^1,2, XU Xinyu⁴, ZHAO Yongqi⁴, CAI Jianqing⁵

1. Chinese Academy of Surveying and Mapping, Beijing 100036, China;
2. Key Laboratory of Surveying and Mapping Science and Geospatial Information Technology, Ministry of Natural Resources, Beijing 100036, China;
3. Beijing Key Laboratory of Urban Spatial Information Engineering, Beijing 100038, China;
4. School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China;
5. Institute of Geodesy Universit't Stuttgart, Stuttgart 70174, Germany

Received:2020-02-13 Revised:2020-05-29 Published:2020-08-25
Supported by:
The Major Project of High-Resolution Earth Observation System(No. 42-Y20A09-9001-17/18);The National Natural Science Foundation of China (Nos. 41574019;41774020);The German Academic Exchange Service (DAAD) Thematic Network Project (No. 57421148);The Basal Research Fund of Chinese Academy of Surveying and Mapping (No. AR1922);Beijing Key Laboratory of Urban Spatial Information Engineering (No. 2019211);The Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, Wuhan University (No. 18-01-05)

摘要/Abstract

摘要： 不同于当前广泛使用的空域法、时域法、直接解法，本文尝试采用Torus方法处理GOCE实测数据，利用71 d的GOCE卫星引力梯度数据反演了200阶次GOCE地球重力场模型，实现了对参考模型的精化。首先，采用Butterworth零相移滤波方法加移去—恢复技术，处理引力梯度观测值中的有色噪声，并利用泰勒级数展开和Kriging方法对GOCE卫星引力梯度数据进行归算和格网化，计算得到了名义轨道上格网点处的引力梯度数据。然后，利用2D-FFT技术和块对角最小二乘方法处理名义轨道上数据，获得了200阶次的GOCE地球重力场模型GOCE_Torus。利用中国和美国的GPS/水准数据进行外部检核结果说明，GOCE_Torus与ESA发布的同期模型的精度相当；GOCE_Torus模型与200阶次的EGM2008模型相比，在美国区域精度相当，但在中国区域精度提高了4.6 cm，这充分体现了GOCE卫星观测数据对地面重力稀疏区的贡献。Torus方法拥有快速高精度反演卫星重力场模型的优势，可以在重力梯度卫星的设计、误差分析及在轨快速评估等方面得到充分应用。

关键词: GOCE, 地球重力场, Torus方法, 精化

Abstract: Currently, space-wise method, time-wise method and direct-wise method are widely used to deal with the GOCE satellite real observations. Different from these conventional methods, this paper proposes one Torus approach. As one refinement of the reference model, the GOCE earth gravity field model complete to 200 d/o is recovered using 71-day GOCE satellite gradiometry observations with Torus approach. Firstly, a Butterworth filter with zero-phase combining with remove-restore is used to deal with the colored noise in GOCE satellite gravitational gradient observations. These observations are reduced from the real orbits to nominal orbits using Taylor series expansion and Kriging is used for gridding on the nominal orbits. Then, 2D-FFT and the block-diagonal least-square adjustment are performed to process the gridded observations, so that the GOCE earth gravity field model complete to 200 d/o named as GOCE_Torus can be realized efficiently. GOCE_Torus is assessed with the GPS/leveling data in China and USA, respectively. GOCE_Torus and the first generation model released by ESA have a similar accuracy. In USA, EGM2008 complete to 200 d/o reveals a similar accuracy with GOCE_Torus. But in China, GOCE_Torus is improved by 4.6 cm due to the contribution of GOCE gravity gradiometry data. With the application of Torus approach the satellite gravity field model can be recovered quickly and accurately, and therefore this approach is useful in mission design, error analysis and in-orbit evaluation of gravity gradient satellites.

Key words: GOCE, earth's gravity field, Torus approach, refinement

中图分类号:

P223

刘焕玲, 文汉江, 徐新禹, 赵永奇, 蔡剑青. GOCE实测数据反演高阶重力场模型的Torus方法[J]. 测绘学报, 2020, 49(8): 965-973.

LIU Huanling, WEN Hanjiang, XU Xinyu, ZHAO Yongqi, CAI Jianqing. Torus approach in high degree gravity field model determination from GOCE satellite gradiometry observations[J]. Acta Geodaetica et Cartographica Sinica, 2020, 49(8): 965-973.

参考文献

[1] RUMMEL R, VAN GELDEREN M, KOOP R, et al. Spherical harmonic analysis of satellite gradiometry[M]. Delft:Netherlands Geodetic Commission, 1993.
[2] KOOP R. Global gravity field modelling using satellite gravity gradiometry[M]. Delft:Netherlands Geodetic Commission, 1993.
[3] 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.
[4] 李建成, 徐新禹, 赵永奇, 等. 由GOCE引力梯度张量不变量确定卫星重力模型的半解析法[J]. 武汉大学学报(信息科学版), 2016, 41(1):21-26. LI Jiancheng, XU Xinyu, ZHAO Yongqi, et al. Approach for determining satellite gravity model from GOCE gravitational gradient tensor invariant observations[J]. Geomatics and Information Science of Wuhan University, 2016, 41(1):21-26.
[5] YU Jinhai, WAN Xiaoyun. Recovery of the gravity field from GOCE data by using the invariants of gradient tensor[J]. Science China Earth Sciences, 2013, 56(7):1193-1199.
[6] 苏勇, 范东明. 利用短弧积分法和GOCE轨道数据确定地球重力场模型[J]. 地球物理学进展, 2014, 29(5):2072-2076. SU Yong, FAN Dongming. Gravity field modeling by using the short-arc integral approach and GOCE orbits[J]. Progress in Geophysics, 2014, 29(5):2072-2076.
[7] 刘晓刚, 吴杉, 王献民, 等. GOCE卫星重力测量中有色噪声滤波器设计[J]. 地球物理学进展, 2012, 27(3):856-860. LIU Xiaogang, WU Shan, WANG Xianmin, et al. Design of color noise filter in GOCE satellite gravimetry mission[J]. Progress in Geophysics, 2012, 27(3):856-860.
[8] 万晓云, 于锦海, 曾艳艳. GOCE引力梯度的频谱分析及滤波[J]. 地球物理学报, 2012, 55(9):2909-2916. WAN Xiaoyun, YU Jinhai, ZENG Yanyan. Frequency analysis and filtering processing of gravity gradients data from GOCE[J]. Chinese Journal of Geophysics, 2012, 55(9):2909-2916.
[9] SNEEUW N. A semi-analytical approach to gravity field analysis from satellite observations[D]. Munich, Germany:Technische Universität München, 2000.
[10] SNEEUW N. Space-wise, time-wise, torus and rosborough representations in gravity field modelling[J]. Space Science Reviews, 2003, 108(1-2):37-46.
[11] XU Chen. The torus-based semi-analytical approach in spaceborne gravimetry[D]. Calgary:University of Calgary, 2008.
[12] 刘焕玲, 文汉江, 徐新禹, 等. 利用GOCE模拟观测反演重力场的Torus法[J]. 测绘学报, 2015, 44(9):965-972. DOI:10.11947/j.AGCS.2015.20150110. LIU Huanling, WEN Hanjiang, XU Xinyu, et al. Torus approach in gravity field determination from simulated GOCE observations[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(9):965-972. DOI:10.11947/j.AGCS.2015.20150110.
[13] HWANG C, KAO Y C. Spherical harmonic analysis and synthesis using FFT:application to temporal gravity variation[J]. Computers & Geosciences, 2006, 32(4):442-451.
[14] GRÜBER C. Kugelfunktionen und analyse heterogener schweredaten im spektralbereich[D]. Berlin, Germany:Technische Universität, 2008.
[15] GRÜBER C, NOVÁK P, SEBERA J. FFT-based high-performance spherical harmonic transformation[J]. Studia Geophysica et Geodaetica, 2011, 55(3):489-500.
[16] BOUMAN J, RISPENS S, GRÜBER T, et al. Preprocessing of gravity gradients at the GOCE high-level processing facility[J]. Journal of Geodesy, 2009, 83:659-678.
[17] 刘焕玲, 赵永奇, 文汉江等. GOCE卫星引力梯度数据滤波方法研究[J]. 测绘科学, 2019(6):66-71. LIU Huanling, ZHAO Yongqi, WEN Hanjiang, et al. Study of filtering method for GOCE satellite gradiometry observations[J]. Science of Surveying and Mapping, 2019(6):66-71.
[18] YI W, RUMMEL R, GRÜBER T. Gravity field contribution analysis of GOCE gravitational gradient components[J]. Studia Geophysica et Geodaetica, 2013, 57(2):174-202.
[19] FUCHS M J, BOUMAN J. Rotation of GOCE gravity gradients to local frames[J]. Geophysical Journal International, 2011, 187(2):743-753.
[20] VAN GELDEREN M, KOOP R. The use of degree variances in satellite gradiometry[J]. Journal of Geodesy, 1997, 71(6):337-343.
[21] SNEEUW N, VAN GELDEREN M. The polar gap[J]. Geodetic Boundary Value Problems in View of the One Centimeter Geoid, 1997.
[22] KAULA W M. Theory of satellite geodesy[M]. Waltham Massachusetts, Blaisdell Publishing Company, 1966.
[23] METZLER B, PAIL R. GOCE data processing:the spherical cap regularization approach[J]. Studia Geophysica et Geodaetica, 2005, 49(4):441-462.
[24] 徐新禹,李建成,王正涛等. Tikhonov正则化方法在GOCE重力场求解中的模拟研究[J]. 测绘学报,2010, 35(5):465-470. XU Xinyu, LI Jianchen, 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, 35(5):465-470.
[25] KAULA W M. Theory of satellite geodesy[M]. Waltham, MA:Blaisdell Publishing Company, 1966.
[26] METZLER B, PAIL R. GOCE data processing:the spherical cap regularization approach[J]. Studia Geophysica et Geodaetica, 2005, 49(4):441-462.
[27] 徐新禹, 李建成, 王正涛, 等. Tikhonov正则化方法在GOCE重力场求解中的模拟研究[J]. 测绘学报, 2010, 39(5):465-470. 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.

GOCE实测数据反演高阶重力场模型的Torus方法

Torus approach in high degree gravity field model determination from GOCE satellite gradiometry observations

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[1]	潘娟霞, 邹贤才. GOCE引力梯度内部校准方法[J]. 测绘学报, 2022, 51(2): 192-200.
[2]	陈鑑华, 张兴福, 沈云中, 陈秋杰, 李伟超. GOCE卫星重力梯度观测值的时变重力场变化改正及影响分析[J]. 测绘学报, 2021, 50(3): 324-332.
[3]	邓非, 陈欣, 颜青松, 曲英杰. 带线约束的摄影测量网格变分精化算法[J]. 测绘学报, 2020, 49(4): 469-479.
[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]	刘宸, 刘长建, 冯绪, 许岭峰, 杜莹. 适用于不同尺度区域的Klobuchar-like电离层模型[J]. 测绘学报, 2016, 45(S2): 54-63.
[11]	孟祥超, 万晓云, 于锦海, 朱永超, 冯炜. 利用引力梯度数据计算径向梯度的优化方法[J]. 测绘学报, 2016, 45(7): 775-781.
[12]	耿弘毅, 龚志辉, 高超, 杨戬峰, 唐可. 天绘一号三线阵相机在轨几何参数精化[J]. 测绘学报, 2016, 45(3): 282-290.
[13]	刘焕玲, 文汉江, 徐新禹, 朱广彬. 利用GOCE模拟观测反演重力场的Torus法[J]. 测绘学报, 2015, 44(9): 965-972.
[14]	周浩, 罗志才, 钟波, 陆飚. 利用最小二乘直接法反演卫星重力场模型的MPI并行算法[J]. 测绘学报, 2015, 44(8): 833-839.
[15]	叶周润, 柳林涛, 梁星辉. Moho面扰动重力梯度信息的提取[J]. 测绘学报, 2015, 44(6): 609-615.