一种高效的基于逆Stokes公式的海域重力异常计算方法  

Efficient method of calculating the marine gravity anomalies based on the inverse Stokes formula

作　　者：王纯 孟小红[1] 王俊[1] 方圆[1] 干永康 WANG Chun;MENG Xiao-hong;WANG Jun;FANG Yuan;GAN Yong-kang(School of Geophysics and Information Technology,China University of Geosciences(Beijing),Beijing 100083,China)

机构地区：[1]中国地质大学(北京)地球物理与信息技术学院,北京100083

出　　处：《地球物理学进展》2020年第6期2384-2389,共6页Progress in Geophysics

基　　金：中国地质调查局项目(DD20191006);国家重点研发项目课题(2016YFC0303001)联合资助。

摘　　要：利用卫星测高数据解算重力场是获取大范围海域重力异常的主要途径.然而,由于涉及到全球范围的多重积分,传统串行计算方法效率低下.针对上述问题,本文研究提出一种高效的海域重力异常计算方法.该快速方法利用逆Stokes公式进行重力场解算,同时针对其中不同点位多重循环计算的效率低下问题提出了一种基于多核CPU技术的并行化改进方法,该方法通过将程序合理地划分为CPU不同线程,在运算过程中合理调配计算机资源以完成海域重力异常的高效并行计算.文章利用EGM2008模型提供的大地水准面数据进行了快速算法测试,并将解算出的重力异常与船测重力异常进行了比对,表明该快速算法能够准确高效地求取出研究区重力异常,计算效率相比较于传统串行算法提升30%以上.本文提出的并行化方法对解算大范围海域重力异常和精化地球重力场模型具有积极意义.Marine gravity anomalies over a large area are usually calculated from the satellite altimetry data. However, numerous conventional sequential methods are time-consuming due to the multiple integration at a global scale. To handle this issue, this paper proposes a fast method for calculating the gravity anomalies from the satellite altimetry data, which is based on the inverse Stokes formula. The proposed fast method employs the multi-core Central Processing Unit technology for parallel computing, with which the most inefficient part of the conventional method can be accelerated. Computer resources are efficiently allocated and seized upon according to divide the program reasonably into multi-threads of CPU in the process of operation, and the parallel operation of marine gravity anomaly can be completed. The geoid data provided by the EGM2008 model are utilized to inverse the research area gravity anomalies to test the fast method. And the corresponding, gravity anomalies provided by the shipborne measurement are used for verifying correctness of the fast method. Numerous testing results show that the fast method can calculate the gravity anomalies in the study area accurately with much less time, which is of positive significance for the satellite altimetry data inversion and refining the Earth gravity field model, and the computation efficiency can be improved by more than 30% compared with the conventional sequential method.

关 键 词：逆Stokes公式 大地水准面高 重力异常 多核CPU加速 

分 类 号：P631[天文地球—地质矿产勘探] P738[天文地球—地质学]