Robust 3D Joint Inversion of Gravity and Magnetic Data: A High-Performance Computing Approach

Abstract

One of the fundamental challenges in geophysics is the calculation of distribution models for physical properties in the subsurface that accurately reproduce the measurements obtained in the survey and are geologically plausible in the context of the study area. This is known as inverse modeling. Performing a 3D joint inversion of multimodal geophysical data is a computationally intensive task. Additionally, since it involves a modeling process, finding a solution that matches the desired characteristics requires iterative calculations, which can take days or even weeks to obtain final results. In this paper, we propose a robust numerical solution for 3D joint inversion of gravimetric and magnetic data with Gramian-based structural similarity and structural direction constraints using parallelization as a high-performance computing technique, which allows us to significantly reduce the total processing time based on the available Random-Access Memory (RAM) and Video Random-Access Memory (VRAM)and improve the efficiency of interpretation. The solution is implemented in the high-level programming languages Fortran and Compute Unified Device Architecture (CUDA) Fortran, capable of optimal resource management while being straightforward to implement. Through the analysis of performance and computational costs of serial, parallel, and hybrid implementations, we conclude that as the inversion domain expands, the processing speed could increase from 4× up to 100× times faster, rendering it particularly advantageous for applications in larger domains. We tested our algorithm with two synthetic data sets and field data, showing better results than standard separate inversion. The proposed method will be useful for joint geological and geophysical interpretation of gravimetric and magnetic data used in exploration geophysics for example minerals, ore, and petroleum search and prospecting. Its application will significantly increase the reliability of physical-geological models and accelerate the process of data processing.

Fast Magnetization Vector Inversion Method with Undulating Observation Surface in Spherical Coordinate for Revealing Lunar Weak Magnetic Anomaly Feature

“…Thus, a three-dimensional magnetic vector inversion system in a spherical coordinate system is needed. When the inversion data amount is large, in order to improve the practicality of the inversion algorithm, computer hardware (GPU and CPU) can be utilized for parallel computation to increase the inversion speed [10][11][12][13]. Due to the large amount of large region satellite data and the slow inversion speed, Zhao et al [14,15] proposed a highly efficient gravity inversion method based on the block-Toeplitz Toeplitz-block (BTTB) structure and the fast Fourier transform (FFT) method for spherical coordinate systems that reduces the inversion time and improves the practicality of satellite data inversion.…”