测绘学报 ›› 2015, Vol. 44 ›› Issue (8): 833-839.doi: 10.11947/j.AGCS.2015.20140396

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

利用最小二乘直接法反演卫星重力场模型的MPI并行算法

周浩1, 罗志才1,2, 钟波1,2, 陆飚1   

  1. 1. 武汉大学测绘学院, 湖北 武汉 430079;
    2. 地球空间环境与大地测量教育部重点实验室, 湖北 武汉 430079
  • 收稿日期:2014-07-23 修回日期:2015-04-20 出版日期:2015-09-20 发布日期:2015-09-02
  • 作者简介:周浩(1987-),男,博士生,主要从事物理大地测量和卫星重力学研究。E-mail:zhouhao2013@whu.edu.cn
  • 基金资助:
    国家973计划(2013CB733302);国家自然科学基金(41104014;41131067);地理信息工程国家重点实验室开放基金(SKLGIE2013-M-1-3);大地测量与地球动力学国家重点实验室开放基金(SKLGED2015-1-3-E);地球空间环境与大地测量教育部重点实验室开放基金(13-02-05)

MPI Parallel Algorithm in Satellite Gravity Field Model Inversion on the Basis of Least Square Method

ZHOU Hao1, LUO Zhicai1,2, ZHONG Bo1,2, LU Biao1   

  1. 1. School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China;
    2. Key Lab of Geospace Environment and Geodesy, Ministry of Education, Wuhan University, Wuhan 430079, China
  • Received:2014-07-23 Revised:2015-04-20 Online:2015-09-20 Published:2015-09-02
  • Supported by:
    The National Basic Research Program of China(973 Program)(No. 2013CB733302);The National Natural Science Foundation of China(Nos. 41104014;41131067);The Open Research Fund Program of the State Key Laboratory of Geo-information Engineering(No. SKLGIE2013-M-1-3) The Open Research Fund Program of the State Key Laboratory of Geodesy and Earth's Dynamics(No. SKLGED2015-1-3-E) The Open Research Fund Program of the Key Laboratory of Geospace Environment and Geodesy, Ministry of Education(No. 13-02-05)

摘要: 针对海量卫星重力数据反演高阶次地球重力场模型的密集型计算任务与高内存耗用问题,基于MPI实现了最小二乘直接法恢复高阶次位系数的并行算法。引入并行读写、分块存储与分块计算等方式完成了设计矩阵的构建、法方程的形成与求解等密集型计算任务的并行算法,数值计算结果表明三者的并行相对效率峰值可分别达到95%、68%、63%。利用GOCE轨道跟踪和径向扰动重力梯度数据(共518 400个历元)分别反演了120、240阶次地球重力场模型,计算时间仅为40 min、7 h,内存耗用峰值仅为290 MB、1.57 GB;采用与GOCE同等噪声水平的观测数据恢复的重力场模型精度与GOCE已发布模型的解算精度相一致,联合GRACE和GOCE的解算模型能够实现二者独立信息的频谱互补,表明本文方法可高效稳定地恢复高阶次地球重力场模型。

关键词: MPI, 并行算法, 卫星重力场模型, Gauss-Jordan算法, GOCE

Abstract: In order to solve the intensive computing tasks and high memory demand problem in satellite gravity field model inversion on the basis of huge amounts of satellite gravity observations, the parallel algorithm for high truncated order and degree satellite gravity field model inversion with least square method on the basis of MPI was introduced. After analyzing the time and space complexity of each step in the solving flow, the parallel I/O, block-organized storage and block-organized computation algorithm on the basis of MPI are introduced to design the parallel algorithm for building design matrix, establishing and solving normal equation, and the simulation results indicate that the parallel efficiency of building design matrix, establishing and solving normal equation can reach to 95%, 68%and 63% respectively. In addition, on the basis of GOCE simulated orbits and radial disturbance gravity gradient data(518 400 epochs in total), two earth gravity models truncated to degree and order 120, 240 are inversed, and the relative computation time and memory demand are only about 40 minutes and 7 hours, 290 MB and 1.57 GB respectively. Eventually, a simulation numerical calculation for earth gravity field model inversion with the simulation data, which has the equivalent noise level with GRACE and GOCE mission, is conducted. The accuracy of inversion model has a good consistent with current released model, and the combined mode can complement the spectral information of each individual mission, which indicates that the parallel algorithm in this paper can be applied to inverse the high truncated degree and order earth gravity model efficiently and stably.

Key words: MPI, parallel algorithm, satellite gravity field model, Gauss-Jordan method, GOCE

中图分类号: