Integrated gravity and gravity gradient 3-D inversion can effectively improve the resolution of inversion results. This study aims to find a better way to extract geological information from gravity and gravity gradient data and to obtain the recovered model with a higher resolution. We present three different methods for integrating gravity and gravity gradient data in inversion: directly integrating gravity and gravity gradient data (DIGG), the sequential inversion of gravity and gravity gradient data (SIGG), and integrating gravity and gravity gradient components using a weighting matrix (IGGW). In addition, we propose a new discrete method for smooth inversion. We discretize the stabilizing functional into two parts. One is the square of the 2-norm of the model parameter, and the other one is the square of the Frobenius norm of the gradient of the model parameter. We adopt a conjugate gradient algorithm (CGA) to solve the objective function. We used the DIGG, SIGG, and IGGW inversion methods in a dip-slab model with multiple anomalous bodies. The inversion results show that by integrating gravity and gravity gradient in inversion, more reasonable results are obtained than by inverting a single data type. Among the presented methods, the IGGW method achieved the best performance in using information contained in gravity and gravity gradient data, followed by the SIGG method, and finally, by the DIGG method. We used the methods on actual data from Vinton salt dome, southwestern Louisiana, USA. The obtained results support the abovementioned conclusions.

The subsurface three-dimensional (3D) model of density distribution is obtained by solving an under-determined linear equation that is established by gravity data. Here, we describe a new fast gravity inversion method to recover a 3D density model from gravity data. The subsurface will be divided into a large number of rectangular blocks, each with an unknown constant density. The gravity inversion method introduces a stabiliser model norm with a depth weighting function to produce smooth models. The depth weighting function is combined with the model norm to counteract the skin effect of the gravity potential field. As the numbers of density model parameters is NZ (the number of layers in the vertical subsurface domain) times greater than the observed gravity data parameters, the inverse density parameter is larger than the observed gravity data parameters. Solving the full set of gravity inversion equations is very time-consuming, and applying a new algorithm to estimate gravity inversion can significantly reduce the number of iterations and the computational time. In this paper, a new symmetric successive over-relaxation (SSOR) iterative conjugate gradient (CG) method is shown to be an appropriate algorithm to solve this Tikhonov cost function (gravity inversion equation). The new, faster method is applied on Gaussian noise-contaminated synthetic data to demonstrate its suitability for 3D gravity inversion. To demonstrate the performance of the new algorithm on actual gravity data, we provide a case study that includes ground-based measurement of residual Bouguer gravity anomalies over the Humble salt dome near Houston, Gulf Coast Basin, off the shore of Louisiana. A 3D distribution of salt rock concentration is used to evaluate the inversion results recovered by the new SSOR iterative method. In the test model, the density values in the constructed model coincide with the known location and depth of the salt dome.

Recent advances in drone-based aeromagnetic survey technology face limitations due to existing interference models and calibration methods cannot adequately compensate for the complex additional magnetic interference on in-cabin magnetometers generated by the on-board electronic (OBE) systems, which require the magnetometer to be mounted on a stinger or wing far from the cabin, restricting drone design and size. To solve this, we propose a compensation method for in-cabin magnetometers. Specifically, in the preprocessing stage, we propose a Recursive Empirical Wavelet Transform Reconstruction algorithm (REWT-R) for more precise current denoising. Then, an Iterative Partial Correlation Screening algorithm (IPCS) is introduced to screen current data for model calibration, avoiding multicollinearity issues. Finally, we eliminate the inaccurate bandpass filtering operation in the calibration stage and improve the loss function with a decomposed weighting scheme, enhancing its global optimization and generalization capabilities. Experiments on public datasets validated the effectiveness of the proposed method at each stage. Meanwhile, compared to the latest methods, our proposed in-cabin magnetic interference compensation method reduces the RMSE by 15 times.