一种高分辨率全球重力场模型构建的高性能科学计算方法

梁伟; 徐新禹; 李建成; 李朋远; 黄剑; 朱天林; 王东亚; 苑刚

doi:10.13203/j.whugis20240321

您当前的位置：

首页 >

文章列表页 >

一种高分辨率全球重力场模型构建的高性能科学计算方法

大地测量学与测量工程 | 更新时间：2026-11-05

- 一种高分辨率全球重力场模型构建的高性能科学计算方法
- A High Performance Scientific Computation Method for High‑Resolution Earth’s Gravity Field Modeling
- 武汉大学学报（信息科学版） 2026年51卷第5期页码：935-943
- 作者机构：
  
  1.北京跟踪与通信技术研究所，北京，100094
  2.中南大学地球科学与信息物理学院，湖南长沙，410083
  3.武汉大学测绘学院，湖北武汉，430079
- 作者简介：
  
  梁伟，博士，助理研究员，主要从事空间态势感知、高分辨率重力场模型构建研究。wliang@whu.edu.cn
  苑刚，博士，研究员。yuangang1978@126.com
- 基金信息：
  
  国家自然科学基金(42388102;42192533;42074019)
- DOI：10.13203/j.whugis20240321
  中图分类号： P223
- 收稿：2025-07-15，
  
  纸质出版：2026-05-05
- 稿件说明：
移动端阅览
梁伟,徐新禹,李建成，等.一种高分辨率全球重力场模型构建的高性能科学计算方法[J].武汉大学学报（信息科学版）,2026,51(5):935-943.

LIANG Wei,XU Xinyu,LI Jiancheng,et al.A High Performance Scientific Computation Method for High‑Resolution Earth’s Gravity Field Modeling[J].Geomatics and Information Science of Wuhan University,2026,51(5):935-943.
梁伟,徐新禹,李建成，等.一种高分辨率全球重力场模型构建的高性能科学计算方法[J].武汉大学学报（信息科学版）,2026,51(5):935-943. DOI：10.13203/j.whugis20240321

LIANG Wei,XU Xinyu,LI Jiancheng,et al.A High Performance Scientific Computation Method for High‑Resolution Earth’s Gravity Field Modeling[J].Geomatics and Information Science of Wuhan University,2026,51(5):935-943. DOI：10.13203/j.whugis20240321

摘要

利用多源重力观测数据构建高分辨率全球重力场模型是重力场研究的基本任务之一。将半解析方法引入高分辨率重力场模型构建，同时采用开放多处理（open multi-processing，OpenMP）+信息传递接口（message passing interface，MPI）并行计算思路，设计了一种高分辨率全球重力场模型构建的高性能科学计算方法，该方法可以在大型集群计算系统下分布式运行。数值实验验证了计算方法的计算精度和效率，结果表明：半解析方法可有效减少模型构建的计算量和内存消耗，解算30′分辨率重力场模型的计算效率是传统块对角最小二乘方法的7倍；使用OpenMP+MPI混合并行算法的计算效率是串行程序的18倍（20个计算核心）；利用4个节点的80个中央处理器可在1 d内完成2′高分辨率重力场模型的解算。

Abstract

Objective

The computation task in high resolution Earth's gravity field modeling is very huge

as the dataset and the parameters are dramatically. We focus on high performance and scientific computation method in high‑resolution Earth's gravity field modeling to make model determination feasible.

Methods

The basic methods we used for high performance and computing are the semi-analysis method

open multi-processing (OpenMP) and message passing interface (MPI) parallel computing method. A new high performance and scientific computing method is designed. The semi-analysis method is also known as the lumped-coefficients method. The OpenMP and MPI technique are to make use of the huge-computing machine. And afterward the method is validated by numerical experiments on computing error and computing efficiency.

Results

The semi-analysis can effectively reduce the computation and storage consumption in high‑resolutoin gravity field modeling. The time consumption of OpenMP+MPI hybrid parallel computing is only 1/18.4 that of sequential computing with 20 cores. The proposed computing method with semi-analysis and OpenMP+MPI computing can be used on super computers for Earth's gravity field modeling

and the modeling of coefficients up to degree and order 5 399 can be completed within one day using 80 central processing unit on 4 nodes.

Conclusions

High performance and scientific computing is useful for high-resolution gravity field modeling

the method need to be improved tested for finer resolution field gravity field modeling.

关键词

Keywords

references

PAVLIS N K ， HOLMES S A ， KENYON S C ， et al . The Development and Evaluation of the Earth Gravitational Model 2008 （EGM2008）［J］. Journal of Geophysical Research： Solid Earth ， 2012 ， 117 （ B4 ）： 2011 JB008916.

FÖRSTE C . EIGEN-6C4： The Latest Combined Global Gravity Field Model Including GOCE Data up to Degree and Order 1949 of GFZ Potsdam and GRGS Toulouse ［C］// EGU General Assembly ， Vienna，Austria ， 2014 .

FECHER T ， PAIL R ， GRUBER T ， et al . GOCO05c： A New Combined Gravity Field Model Based on Full Normal Equations and Regionally Varying Weighting ［J］. Surveys in Geophysics ， 2017 ， 38 （ 3 ）： 571 - 590 .

ZINGERLE P ， PAIL R ， GRUBER T ， et al . The Combined Global Gravity Field Model XGM2019e ［J］. Journal of Geodesy ， 2020 ， 94 （ 7 ）： 66 .

LIANG W ， LI J C ， XU X Y ， et al . A High-Resolution Earth's Gravity Field Model SGG-UGM-2 from GOCE， GRACE， Satellite Altimetry， and EGM2008 ［J］. Engineering ， 2020 ， 6 （ 8 ）： 860 - 878 .

COOLEY J W ， TUKEY J W . An Algorithm for the Machine Calculation of Complex Fourier Series ［J］. Mathematics of Computation ， 2019 ， 90 ： 297 - 301 .

COLOMBO O L . Numerical Methods for Harmonic Analysis On the Sphere ［R］. Columbus ： The Ohio State University ， 1981 .

SNEEUW N . Global Spherical Harmonic Analysis by Least-Squares and Numerical Quadrature Methods in Historical Perspective ［J］. Geophysical Journal International ， 1994 ， 118 （ 3 ）： 707 - 716 .

PAVLIS N K . Modeling and Estimation of a Low Degree Geopotential Model from Terestrial Gravity Data ［R］. Columbus ： The Ohio State University Columbus ， 1988 .

RUMMEL R ， VAN GELDEREN M ， KOOP R ， et al . Spherical Harmonic Analysis of Satellite Gradio‑ metry ［M］. New York ： Nederlandse Commissie voor Geodesie ， 1993 .

SNEEUW N J . A Semi-analytical Approach to Gravity Field Analysis from Satellite ［D］. Munich， Germany ： Technische Universität München ， 2000 .

PREIMESBERGER T ， PAIL R . GOCE Quick-Look Gravity Solution： Application of the Semianalytic Approach in the Case of Data Gaps and Non-Repeat Orbits ［J］. Studia Geophysica et Geodaetica ， 2003 ， 47 （ 3 ）： 435 - 453 .

徐新禹 . 卫星重力梯度及卫星跟踪卫星数据确定地球重力场的研宄［D］. 武汉：武汉大学 . 2008 .

XU Xinyu . Study of Determining the Earth’s Gravity Field from Satellite Gravity Gradient and Satellite-to-Satellite Tracking Data ［D］. Wuhan ： Wuhan University ， 2008 .

CHAPMAN B ， JOST G ， PAS R . Using OpenMP： Portable Shared Memory Parallel Programming ［M］. Massachusetts ： MIT Press ， 2007 .

GROPP W ， LUSK E ， SKJELLUM A . Using MPI- Portable Parallel Programming with the Message-Passing Interface ［M］. Massachusetts ： MIT Press ， 1999 .

GROPP W ， LUSK E ， THAKUR R . Using MPI-2 Advanced features of the Message-Passing Interface ［M］. Massachusetts ： MIT Press ， 1999 .

周浩，钟波，罗志才，等 . OpenMP并行算法在卫星重力场模型反演中的应用［J］. 大地测量与地球动力学， 2011 ， 31 （ 5 ）： 123 - 127 .

ZHOU Hao ， ZHONG Bo ， LUO Zhicai ， et al . Application of Parallel Algorithms Based on Openmp to Satellite Gravity Field Recovery ［J］. Journal of Geodesy and Geodynamics ， 2011 ， 31 （ 5 ）： 123 - 127 .

MAYER-GÜRR T ， BEHZADPOUR S ， EICKER A ， et al . GROOPS： A Software Toolkit for Gravity Field Recovery and GNSS Processing ［J］. Computers & Geosciences ， 2021 ， 155 ： 104864 .

王鸿琰，关雪峰，吴华意 . 一种面向CPU/GPU异构环境的协同并行空间插值算法［J］. 武汉大学学报（信息科学版）， 2017 ， 42 （ 12 ）： 1688 - 1695 .

WANG Hongyan ， GUAN Xuefeng ， WU Huayi . A Collaborative Parallel Spatial Interpolation Algorithmon Oriented Towards the Heterogeneous CPU/GPU System ［J］. Geomatics and Information Scien‑ce of Wuhan University ， 2017 ， 42 （ 12 ）： 1688 - 1695 .

陈正生，吕志平，崔阳，等 . 大规模GNSS数据的分布式处理与实现［J］. 武汉大学学报（信息科学版）， 2015 ， 40 （ 3 ）： 384 - 389 .

CHEN Zhengsheng ， LÜ Zhiping ， CUI Yang ， et al . Implementation of Distributed Computing with Large-Scale GNSS Data ［J］. Geomatics and Information Science of Wuhan University ， 2015 ， 40 （ 3 ）： 384 - 389 .

BROCKMANN J M ， ROESE-KOERNER L ， SCHUH W D . A Concept for the Estimation of High-Degree Gravity Field Models in a High Performance Computing Environment ［J］. Studia Geophysica et Geodaetica ， 2014 ， 58 （ 4 ）： 571 - 594 .

BLACKFORD L S ， CHOI J ， CLEARY A ， et al . ScaLAPACK Users' Guide ［M］. Philadelphia ： Society for Industrial and Applied Mathematics ， 1997 .

HOFMANN-WELLENHOF B ， MORITZ H . Physical Geodesy ［M］. Vienna ： Springer ， 2006 .

浏览量

下载量

CSCD

文章被引用时，请邮件提醒。

提交

工具集

关联资源

呼包鄂榆城市群草地脆弱性时空演变及因子可解释性分析

土壤水分机器学习遥感反演：进展与挑战

关于智慧海图的思考

时空智能学：理论与应用

几何与频域双约束的单步扩散海底地形超分辨率重建方法