Abstract:An efficient algorithm for the Lp-norm joint inversion of gravity and magnetic data using the cross-gradient constraint is presented. The presented framework incorporates stabilizers that use Lp-norms (0 <= p <= 2) of the model parameters, and/or the gradient of the model parameters. The formulation is developed from standard approaches for independent inversion of single data sets, and, thus, also facilitates the inclusion of necessary model and data weighting matrices, for example, depth weighting and hard constraint matrices. Using the block Toeplitz Toeplitz block structure of the underlying sensitivity matrices for gravity and magnetic models, when data are obtained on a uniform grid, the blocks for each layer of the depth are embedded in block circulant circulant block matrices. Then, all operations with these matrices are implemented efficiently using 2-D fast Fourier transforms, with a significant reduction in storage requirements. The nonlinear global objective function is minimized iteratively by imposing stationarity on the linear equation that results from applying linearization of the objective function about a starting model. To numerically solve the resulting linear system, at each iteration, the conjugate gradient algorithm is used. This is improved for large scale problems by the introduction of an algorithm in which updates for the magnetic and gravity parameter models are alternated at each iteration, further reducing total computational cost and storage requirements. Numerical results using a complicated 3-D synthetic model and real data sets obtained over the Galinge iron-ore deposit in the Qinghai province, north-west (NW) of China, demonstrate the efficiency of the presented algorithm.

Joint Gravity and Gravity Gradient Inversion for Subsurface Object Detection

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

Joint Inversion of Gravity and Gravity Gradient Data Using Smoothed L0 Norm Regularization Algorithm with Sensitivity Matrix Compression

Stable solution of normalized total gravity gradient

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

Joint Inversion of Gravity Gradient Tensor Data Based on L1 and L2 Norms

A Joint Inversion Algorithm for the Establishment of High-Accuracy 3-D Marine Gravity Field

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

Inversion of Gravity Anomaly Based on AlexNet Deep Neural Network

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

Joint inversion of full-tensor gravity gradiometry data based on source growing

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

A Density-Weighted and Cross-Gradient Constrained Joint Inversion Method of Gravity and Vertical Gravity Gradient Data in Spherical Coordinates and Its Application to Lunar Data

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

Cross-Gradient Joint Inversion of DC Resistivity and Gravity Gradient Data: A Multi-Disciplinary Approach for Geoscience, Heritage, and the Built Environment

Gravity Inversion in Spherical Coordinates With Dynamic Re-Weighting Matrix

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

3D cooperative inversion of airborne magnetic and gravity gradient data using deep learning techniques

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

A Deep Learning Gravity Inversion Method Based on a Self-Constrained Network and Its Application

Joint Gravity and Magnetic Inversion Using CNNs' Deep Learning