中文 | English
Hefei University of Technology Logo J. Hefei Univ. Tech. (Nat. Sci.)

Navigation

Computation of marine gravity gradient models based on Stokes' planar approximation function and GPU parallel computing

Download Full Text

基于斯托克斯平面近似函数与 GPU 并行的海洋重力梯度模型计算

Journal Information

J. Hefei Univ. Tech. (Nat. Sci.), Feb.2026, Vol.49 No.2: 253-259

DOI: 10.3969/j.issn.1003-5060.2026.02.018

Authors

BU Jingyu $ ^{1} $, YE Zhourun $ ^{1} $, LIANG Xinghui $ ^{2} $, LIU Jinzhao $ ^{3} $, LIU Lintao $ ^{2} $, WANG Jiachen $ ^{1} $

(1. School of Civil and Hydraulic Engineering, Hefei University of Technology, Hefei 230009, China; 2. State Key Laboratory of Geodesy and Earth's Dynamics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430077, China; 3. The First Monitoring and Application Center, China Earthquake Administration, Tianjin 300180, China)

Abstract and Keywords

Abstract: Compared to other gravitational field elements, the disturbing gravity gradients can reflect more high-frequency information produced by changing irregular Earth. When calculating the disturbing gravity gradient, due to the complexity of the Stokes integral, the integrand is complex, making it difficult to directly apply Newton-Leibniz formula, and the excessive data makes the computation too time-consuming. In order to effectively solve this problem, this paper uses the Gaussian numerical integration method to calculate the disturbing gravity gradient, and uses the compute unified device architecture (CUDA) to realize parallel computing on the graphics processing unit (GPU) side during the calculation process. According to the Laplace's equation, the accuracy of the calculation results can be verified. This paper selects the gravity anomaly data of sea level within $ 3^{\circ} \times 2^{\circ} $ of a certain sea area for calculation, and the results show that the Gaussian numerical integration and CUDA parallel computing method can provide accurate calculation results and improve the computational efficiency.

Keywords: disturbing gravity gradient; gravity anomaly; compute unified device architecture(CUDA) parallel computing; graphics processing unit(GPU); Gaussian numerical integration

作者信息

卜靖宇 $ ^{1} $，叶周润 $ ^{1} $，梁星辉 $ ^{2} $，刘金钊 $ ^{3} $，柳林涛 $ ^{2} $，王嘉琛 $ ^{1} $

(1. 合肥工业大学 土木与水利工程学院, 安徽 合肥 230009; 2. 中国科学院精密测量科学与技术创新研究院 大地测量与地球动力学国家重点实验室, 湖北 武汉 430077; 3. 中国地震局第一监测中心, 天津 300180)

摘要和关键词

摘要: 相对于其他重力场元素,扰动重力梯度能更多地反映变化的不规则地球产生的高频信息。在计算扰动重力梯度时,由于斯托克斯积分较为复杂导致被积函数复杂难以直接用牛顿-莱布尼茨公式计算、且计算的数据量过于庞大导致计算耗时过长。为有效解决该问题,文章使用高斯数值积分解决被积函数复杂的问题,同时利用统一计算设备架构(compute unified device architecture, CUDA)在计算过程中实现了在图形处理器(graphics processing unit, GPU)端的并行计算,根据拉普拉斯方程可以检验计算结果的准确性,并且选取了某海域 $ 3^{\circ}\times2^{\circ} $范围海平面的重力异常数据进行计算。结果表明,使用高斯数值积分以及CUDA并行计算的方法,提供准确计算结果的同时也提高了计算效率。

关键词: 扰动重力梯度；重力异常；CUDA 并行计算；图形处理器（GPU）；高斯数值积分 中图分类号：P223 文献标志码：A 文章编号：1003-5060（2026）02-0253-07

Funding

国家重点研发计划资助项目（2024YFB3908104）；国家自然科学基金资助项目（42430101；42574005）

Previous Next
Back to Contents

Account

Submission Login Reviewer Login Editor Login