动力学法的卫星重力反演算法特点与改进设想

doi:10.11947/j.AGCS.2017.20170380

测绘学报 ›› 2017, Vol. 46 ›› Issue (10): 1308-1315.doi: 10.11947/j.AGCS.2017.20170380

动力学法的卫星重力反演算法特点与改进设想

沈云中

同济大学测绘与地理信息学院, 上海 200092

收稿日期:2017-07-04 修回日期:2017-08-11 出版日期:2017-10-20 发布日期:2017-10-26
作者简介:沈云中(1962-),男,教授,研究方向为大地测量数据处理及其卫星重力和卫星定位应用。E-mail:yzshen@tongji.edu.cn
基金资助:
国家自然科学基金（41474017）；中国科学院战略性先导科技专项（B类）（XDB23030100）；高分遥感测绘应用示范系统资助项目

Algorithm Characteristics of Dynamic Approach-based Satellite Gravimetry and Its Improvement Proposals

SHEN Yunzhong

College of Surveying and Geo-informatics, Tongji University, Shanghai 200092, China

Received:2017-07-04 Revised:2017-08-11 Online:2017-10-20 Published:2017-10-26
Supported by:
The National Natural Science Foundation of China(No. 41474017);The Strategic Priority Research Program of the Chinese Academy of Sciences(No. XD1323030100);The Program of Application and Demonstration System of High Resolution Remote Sensing in Surveying and Mapping.

摘要/Abstract

摘要： 根据卫星轨道计算的积分公式，导出了以参考轨道为初值的线性化解算地球重力场的观测方程，给出了其系数矩阵的积分计算公式，阐明了动力学法本质上是观测值相对于参考轨道的线性摄动方法，因此其变分方程力模型参数的偏导数初值必定为0。在此公式的基础上，分析了动力学法观测方程的主要特点，即线性化误差随轨道弧段增长而快速增大，其观测方程的性质也随弧段增长而变差，且积分计算误差将是下一代重力卫星数据处理的重要瓶颈问题。提出了进一步提高动力学法重力反演精度的方法，主要归结为：改进以几何轨道为初值的线性化方法以减小线性化误差，改变参数化方式以改善观测方程的性质，综合应用解析公式与数值积分公式以提高轨道计算精度。

关键词: 卫星重力反演, 动力学法, 轨道积分, 解析轨道, 线性化方法

Abstract: By using the integration equation for computing satellite orbit, this paper derives the observational equation for computing gravity field model linearized with respect to the reference orbit, provides the integration equation of calculating the design matrices of the observational equation, and clarifies that the dynamic approach is in principle the perturbation method relative to the reference orbit, therefore the initial partial derivatives with respect to the force model parameters must be zero. Based on the derived formulae this paper analyzes the main characteristics of dynamic approach-based observational equation, i.e. the linearization error will rapidly increase and the property of observational equation becomes worse as the integration arc extends longer, and the numerical integration error will be the bottle-neck problem for the data processing of next generation of satellite gravity exploration. Then this paper proposes the methods for improving the accuracy of gravity recovery, which can be summarized as that refining the linearization method relative kinematic orbit to reduce the linearization error, modifying parameterization method to improve the property of observational equation, and combined using analytic formula and numerical integration formula to increase the accuracy of orbit computation.

Key words: satellite gravimetry, dynamic approach, orbit integration, analytic orbit, linearization method

中图分类号:

P223

沈云中. 动力学法的卫星重力反演算法特点与改进设想[J]. 测绘学报, 2017, 46(10): 1308-1315.

SHEN Yunzhong. Algorithm Characteristics of Dynamic Approach-based Satellite Gravimetry and Its Improvement Proposals[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(10): 1308-1315.

参考文献

[1] TAPLEY B D, BETTADPUR S, WATKINS M, et al. The Gravity Recovery and Climate Experiment:Mission Overview and Early Results[J]. Geophysical Research Letters, 2004, 31(9):L09607.
[2] 许厚泽, 陆洋, 钟敏, 等. 卫星重力测量及其在地球物理环境变化监测中的应用[J]. 中国科学(地球科学):2012, 42(6):843-853. XU Houze, LU Yang, ZHONG Min, et al. Satellite Gravity and Its Application to Monitoring Geophysical Environmental Change[J]. Scientia Sinica (Terrae), 2012, 42(6):843-853.
[3] CHEN J L, WILSON C R, TAPLEY B D. Satellite Gravity Measurements Confirm Accelerated Melting of Greenland Ice Sheet[J]. Science, 2006, 313(5795):1958-1960. DOI:10.1126/science.1129007.
[4] FENG Wei, ZHONG Min, LEMOINE J M, et al. Evaluation of Groundwater Depletion in North China Using the Gravity Recovery and Climate Experiment (GRACE) Data and Ground-based Measurements[J]. Water Resources Research, 2013, 49(4):2110-2118.
[5] DAHLE C, FLECHTNER F, GRUBER C, et al. GFZ GRACE Level-2 Processing Standards Document for Level-2 Product Release 0005[R]. Scientific Technical Report STR 12/02. Potsdam:Deutsches GeoForschungsZentrum GFZ, 2012. DOI:10.2312/GFZ.b103-1202-25,GFZ.
[6] BETTADPUR S. Gravity Recovery and Climate Experiment UTCSR Level-2 Processing Standards Document for Level-2 Product Release 0005[M]. Austin:Center for Space Research, University of Texas, 2012.
[7] WATKINS M, YUAN D N. GRACE JPL Level-2 Processing Standards Document for Level-2 Product Release 05[M]. GRACE 327-744(v 5.0). Pasadena, USA:Jet Propulsion Laboratory, 2012.
[8] MAYER-GVRR T, BEHZADPOUR S, ELLMER M, et al. ITSG-Grace2016-Monthly and Daily Gravity Field Solutions from GRACE[M/OL]. GFZ Data Services, 2016, http://doi.org/10.5880/icgem.2016.007.
[9] SAVE H, BETTADPUR S, TAPLEY B D. High-resolution CSR GRACE RL05 Mascons[J]. Journal of Geophysical Research:Solid Earth, 2016, 121(10):7547-7569.
[10] FÖRSTE C, BRUINSMA S, ABRIKOSOV O, et al. EIGEN-6S4 a Time-variable Satellite-only Gravity Field Model to D/O 300 Based on LAGEOS, GRACE and GOCE Data from the Collaboration of GFZ Potsdam and GRGS Toulouse[R].[S.l.]:GFZ, 2016,
[11] ZHOU Hao, LUO Zhicai, ZHOU Zebing, et al. HUST-GRACE2016s:A New GRACE Static Gravity Field Model Derived from a Modified Dynamic Approach over a 13-year Observation Period[J]. Advances in Space Research, 2017, 60(3):597-611.
[12] MAYER-GVRR T, EICKER A, ILK K H. ITG-Grace02s:A GRACE Gravity Field Derived from Short Arcs of the Satellites Orbit[C]//Proceedings of the 1st International Symposium of the International Gravity Field Service "Gravity Field of the Earth". Istanbul:Harita Dergisi, 2007. 193-198.
[13] CHEN Qiujie, SHEN Yunzhong, ZHANG Xingfu, et al. Monthly Gravity Field Models Derived from GRACE Level 1B Data Using a Modified Short-arc Approach[J]. Journal of Geophysical Research, 2015, 120(3):1804-1819.
[14] CHEN Qiujie, SHEN Yunzhong, CHEN Wu, et al. An Improved GRACE Monthly Gravity Field Solution By Modeling the Non-conservative Acceleration and Attitude Observation Errors[J]. Journal of Geodesy, 2016, 90(6):503-523.
[15] DITMAR P, KUZNETSOV V, VAN ECK VAN DER SLUIJS A A, et al. 'DEOS_CHAMP-01C_70':A Model of the Earth's Gravity Field Computed from Accelerations of the CHAMP Satellite[J]. Journal of Geodesy, 2006, 79(10-11):586-601.
[16] LIU X, DITMAR P, SIEMES C, et al. DEOS Mass Transport Model (DMT-1) Based on GRACE Satellite Data:Methodology and Validation[J]. Geophysical Journal International, 2010, 181(2):769-788.
[17] CHEN Qiujie, SHEN Yunzhong, CHEN Wu, et al. A Modified Acceleration-based Monthly Gravity Field Solution from GRACE Data[J]. Geophysical Journal International, 2015, 202(2):1190-1206.
[18] 徐天河, 杨元喜. 利用CHAMP卫星几何法轨道恢复地球重力场模型[J]. 地球物理学报, 2005, 48(2):288-293. XU Tianhe, YANG Yuanxi. CHAMP Gravity Field Recovery Using Kinematic Orbits[J]. Chinese Journal of Geophysics, 2005, 48(2):288-293.
[19] 王长青, 许厚泽, 钟敏, 等. 利用动力学方法解算GRACE时变重力场研究[J]. 地球物理学报, 2015, 58(3):756-766. WANG Changqing, XU Houze, ZHONG Min, et al. An Investigation on GRACE Temporal Gravity Field Recovery Using the Dynamic Approach[J]. Chinese Journal of Geophysics, 2015, 58(3):756-766.
[20] 王正涛, 李建成, 姜卫平, 等. 基于GRACE卫星重力数据确定地球重力场模型WHU-GM-05[J]. 地球物理学报, 2008, 51(5):1364-1371. WANG Zhengtao, LI Jiancheng, JIANG Weiping, et al. Determination of Earth Gravity Field Model WHU-GM-05 Using GRACE Gravity Data[J]. Chinese Journal of Geophysics, 2008, 51(5):1364-1371.
[21] 肖云, 夏哲仁, 王兴涛. 用GRACE星间速度恢复地球重力场[J]. 测绘学报, 2007, 36(1):19-25. DOI:10.3321/j.issn:1001-1595.2007.01.004. XIAO Yun, XIA Zheren, WANG Xingtao. Recovering the Earth Gravity Field from Inter-satellite Range-rate of GRACE[J]. Acta Geodaetica et Cartographica Sinica, 2007, 36(1):19-25. DOI:10.3321/j.issn:1001-1595.2007.01.004.
[22] 宁津生, 钟波, 罗志才, 等. 基于卫星加速度恢复地球重力场的去相关滤波法[J]. 测绘学报, 2010, 39(4):331-337, 343. NING Jinsheng, ZHONG Bo, LUO Zhicai, et al. Decorrelation Filtering Methods for Gravity Field Recovery Based on Satellite Acceleration Data[J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(4):331-337, 343.
[23] XU Peiliang. Position and Velocity Perturbations for the Determination of Geopotential from Space Geodetic Measurements[J]. Celestial Mechanics and Dynamical Astronomy, 2008, 100(3):231-249.
[24] SHEN Y, CHEN Q, HSU H, et al. A Modified Short Arc Approach for Recovering Gravity Field Model[C]//Oral Presentation at the GRACE Science Team Meeting. Austin:Center of Space Research, University of Texas, 2013.
[25] FLECHTNER F, WATKINS M, MORTON P, et al. Status of the GRACE Follow-on Mission[C]//Joint GSTM/SPP Final Colloquium. Potsdam, Germany:[s.n.], 2012.
[26] GRUBER T, MURBÖCK M, PAIL R, et al. Earth System Mass Transport Mission 2 e. Motion2-proposal for an Earth Explorer 9 Mission[C]//Asia Oceania Geoscience Society 13th Annual Meeting, Session SE01. Beijing:AOGS, 2016.
[27] LOOMIS B D, NEREM R S, LUTHCKE S B. Simulation Study of A Follow-on Gravity Mission to GRACE[J]. Journal of Geodesy, 2012, 86(5):319-335. DOI:10.1007/s00190-011-0521-8.
[28] 刘晓刚, 肖云, 李婧, 等. 低低卫——卫跟踪模式中卫星轨道高度和星间距离的指标设计论证[J]. 地球物理学进展, 2013, 28(5):2247-2255. LIU Xiaogang, XIAO Yun, LI Jing, et al. Demonstration on the Indexes Design of Satellite Orbit Height and Inter-satellite Range in the Low-low Satellite-to-satellite Tracking Mode[J]. Progress in Geophysics, 2013, 28(5):2247-2255.
[29] 冉将军, 许厚泽, 沈云中, 等. 新一代GRACE重力卫星反演地球重力场的预期精度[J]. 地球物理学报, 2012, 55(9):2898-2908. RAN Jiangjun, XU Houze, SHEN Yunzhong, et al. Expected Accuracy of the Global Gravity Field for Next GRACE Satellite Gravity Mission[J]. Chinese Journal of Geophysics, 2012, 55(9):2898-2908.
[30] XU Peiliang. Zero Initial Partial Derivatives of Satellite Orbits with Respect to Force Parameters Violate the Physics of Motion of Celestial Bodies[J]. Science in China Series D:Earth Sciences, 2009, 52(4):562-566.
[31] 陈秋杰, 沈云中, 张兴福, 等. 基于GRACE卫星数据的高精度全球静态重力场模型[J]. 测绘学报, 2016, 45(4):396-403. DOI:10.11947/j.AGCS.2016.20150422. CHEN Qiujie, SHEN Yunzhong, ZHANG Xingfu, et al. GRACE Data-based High Accuracy Global Static Earth's Gravity Field Model[J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(4):396-403. DOI:10.11947/j.AGCS.2016.20150422.
[32] XU Guochang. Orbits[M]. Berlin:Springer, 2008.
[33] ELLMER M, MAYER-GVRR T. High Precision Dynamic Orbit Integration for Spaceborne Gravimetry in View of GRACE Follow-on[J]. Advances in Space Research, 2017, 60(1):1-13.

动力学法的卫星重力反演算法特点与改进设想

Algorithm Characteristics of Dynamic Approach-based Satellite Gravimetry and Its Improvement Proposals

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 3

编辑推荐

Metrics

本文评价

[1]	袁俊军, 李凯, 唐成盼, 周善石, 胡小工, 曹建峰. 面向精密位置服务的低轨卫星轨道预报精度分析[J]. 测绘学报, 2022, 51(5): 640-647.
[2]	范磊, 李敏, 宋伟伟, 施闯, 王成. 多系统GNSS卫星轨道快速积分方法[J]. 测绘学报, 2016, 45(S2): 93-100.
[3]	陈秋杰, 沈云中, 张兴福, 陈武, 许厚泽. 基于GRACE卫星数据的高精度全球静态重力场模型[J]. 测绘学报, 2016, 45(4): 396-403.