三维Tesseroid网格模型重力异常正演方法及并行算法  

The approach to gravity forward calculation of 3D Tesseroid mesh model and its parallel algorithm

作　　者：王博[1] 郭良辉[1,2] 崔亚彤 王祥 WANG Bo;GUO Liang-Hui;CUI Ya-Tong;Wang Xiang(School of Geophysics and Information Technology,China University of Geosciences(Beijing),Beijing 100083,China;State Key Laboratory of Geological Processes and Mineral Resources,China Universityof Geosciences(Beijing),Beijing 100083,China)

机构地区：[1]中国地质大学(北京)地球物理与信息技术学院,北京100083 [2]地质过程与矿产资源国家重点实验室中国地质大学(北京),北京100083

出　　处：《物探与化探》2021年第6期1597-1605,共9页Geophysical and Geochemical Exploration

基　　金：国家自然科学基金面上项目(41774098,41974101);中央高校基本科研业务费专项资金;地质过程与矿产资源国家重点实验室科技部专项经费。

摘　　要：三维网格模型的正演计算是重力资料反演的基础。高精度、高效率的正演有利于提高反演解释质量。针对大尺度、地表观测面研究区的高精度、高效率重力正演问题,本文给出球坐标系三维Tesseroid网格模型重力异常正演方法及并行算法。其中,正演算法采用改进的高斯—勒让德积分法实现大尺度、地表观测面的重力异常高精度计算,并行算法采用基于OpenMP的MATLAB任务并行算法实现高效率计算。理论模型和华东岩石圈三维模型数据试验验证了本文方法的有效性。本文方法为高效的大尺度重力场模拟和三维反演提供技术支撑。The forward modeling of a 3 D mesh model is the basis of gravity data inversion.High precision and high efficiency forward modeling is helpful to the improvement of the quality of inversion interpretation.In order to solve the problem of high precision and high efficiency gravity forward modeling based on a large-scale surface observation area,this paper presents the gravity anomaly forward modeling method and parallel algorithm of a 3 D Tesseroid mesh model in the spherical coordinate system.The forward modeling uses the improved Gauss-Legendre Quadrature integration method to realize the high-precision gravity anomaly calculation based on a large-scale surface observation area,and also uses the MATLAB task parallel algorithm based on OpenMP to realize the high-efficiency forward modeling.The test on the 3 D theoretical model and the Eastern China lithospheric model has verified the validity of the proposed method.This method can provide technical support for efficient large-scale gravity field simulation and 3 D inversion.

关 键 词：球坐标系 Tesseroid 重力正演 并行计算 地表观测面 

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