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

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 4x up to 100x 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.