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 3D Joint Inversion of Gravity and Magnetic Data Based on Cross Gradient Constraint

An Efficient Cross-Gradient Joint Inversion Algorithm for Gravity and Magnetic Data Using a Sequential Strategy

3-D Joint Inversion of Gravity and Magnetic Data Using Data-Space and Truncated Gauss–Newton Methods

3-D data-space joint inversion of gravity and magnetic data using a correlation-analysis constraint

An Efficient Alternating Algorithm for the Lp-Norm Cross-Gradient Joint Inversion of Gravity and Magnetic Data Using the 2-D Fast Fourier Transform

Three-Dimensional Joint Inversion Of Gravity And Gravity Gradient Data Based On Data Space And Sparse Constraints

Alternating Joint Inversion of Controlled-Source Electromagnetic and Seismic Data Using the Joint Total Variation Constraint.

Three-dimensional joint inversion of gravity and magnetic data using Gramian constraints and Gauss-Newton method

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

Gravity-magnetic Cross-Gradient Joint Inversion by the Cyclic Gradient Method.

Joint Inversion of 3D Gravity and Magnetic Data under Undulating Terrain Based on Combined Hexahedral Grid

Cross-gradients joint 3D inversion with applications to gravity and magnetic data

Joint Inversion of Multiphysical Parameters Based on a Combination of Cosine Dot-Gradient and Joint Total Variation Constraints.

3-D Joint Inversion of Gravity and Magnetic Data Based on a Deep Learning Network With Automatic Recognition of Structural Similarity

Generalized L$_p$-norm joint inversion of gravity and magnetic data using cross-gradient constraint

Joint Inversion Method of Gravity and Magnetic Analytic Signal Data with Adaptive Unstructured Tetrahedral Subdivision.

Joint Inversion of Gravity and Gravity Gradient Data Based on Cross-Gradient Function

Efficiently Implementing and Balancing the Mixed Lp-Norm Joint Inversion of Gravity and Magnetic Data

Large-scale focusing joint inversion of gravity and magnetic data with Gramian constraint

An efficient joint inversion strategy for 3D seismic travel time and DC resistivity data based on cross-gradient structure constraint

3D joint inversion of magnetotelluric and gravity data based on local correlation constraints