Investigation of magnetic inversion methods in highly magnetic environments under strong self-demagnetization effect

Geophysics ◽  
2017 ◽  
Vol 82 (6) ◽  
pp. J83-J97 ◽  
Author(s):  
Richard A. Krahenbuhl ◽  
Yaoguo Li
Keyword(s):  
Numerical Experiments ◽  
Large Field ◽  
Magnetic Data ◽  
Complex Structures ◽  
Magnetization Direction ◽  
Data Set ◽  
Inversion Methods ◽  
Amplitude Data ◽  
Magnetic Inversion ◽  
Complex Geology

We investigate self-demagnetization effects on magnetic data and develop a comparison of two existing inversion methods as they apply to quantitative interpretation of such data in highly magnetic environments. We begin by evaluating the effect on magnetization direction when susceptibility is a scalar and increases from low values into the realm of self-demagnetization. We show through numerical experiments that susceptibility values of greater than 0.1 SI lead to significant self-demagnetization effects. Second, we show that conventional inversion can perform well for interpretation of self-demagnetization problems with simple source geometries. However, as the geometry becomes more complex in realistically complex problems, this approach can produce poor results and a more robust technique is required. Our numerical experiments indicate that directly inverting amplitude data, which can be derived from total-field magnetic anomaly data and are weakly dependent on magnetization direction, produces superior results when interpreting data from areas with complex geology and high magnetic susceptibilities. We conclude by evaluating the application of our preferred approach on a large field data set exhibiting strong self-demagnetization, multiple source bodies, and complex structures.

Adaptive sampling of potential-field data: A direct approach to compressive inversion

Geophysics ◽  
2014 ◽  
Vol 79 (1) ◽  
pp. IM1-IM9 ◽  
Author(s):  
Nathan Leon Foks ◽  
Richard Krahenbuhl ◽  
Yaoguo Li
Keyword(s):  
Potential Field ◽  
Adaptive Sampling ◽  
Field Data ◽  
Computational Cost ◽  
Large Field ◽  
Magnetic Data ◽  
Data Sampling ◽  
Data Types ◽  
Data Set ◽  
Potential Field Data

Compressive inversion uses computational algorithms that decrease the time and storage needs of a traditional inverse problem. Most compression approaches focus on the model domain, and very few, other than traditional downsampling focus on the data domain for potential-field applications. To further the compression in the data domain, a direct and practical approach to the adaptive downsampling of potential-field data for large inversion problems has been developed. The approach is formulated to significantly reduce the quantity of data in relatively smooth or quiet regions of the data set, while preserving the signal anomalies that contain the relevant target information. Two major benefits arise from this form of compressive inversion. First, because the approach compresses the problem in the data domain, it can be applied immediately without the addition of, or modification to, existing inversion software. Second, as most industry software use some form of model or sensitivity compression, the addition of this adaptive data sampling creates a complete compressive inversion methodology whereby the reduction of computational cost is achieved simultaneously in the model and data domains. We applied the method to a synthetic magnetic data set and two large field magnetic data sets; however, the method is also applicable to other data types. Our results showed that the relevant model information is maintained after inversion despite using 1%–5% of the data.

Comprehensive approaches to 3D inversion of magnetic data affected by remanent magnetization

Geophysics ◽  
2010 ◽  
Vol 75 (1) ◽  
pp. L1-L11 ◽  
Author(s):  
Yaoguo Li ◽  
Sarah E. Shearer ◽  
Matthew M. Haney ◽  
Neal Dannemiller
Keyword(s):  
Mineral Exploration ◽  
Three Dimensional ◽  
Canadian Arctic ◽  
Magnetization Vector ◽  
Magnetic Data ◽  
Inversion Algorithm ◽  
3D Inversion ◽  
Magnetization Direction ◽  
Direct Inversion ◽  
Amplitude Data

Three-dimensional (3D) inversion of magnetic data to recover a distribution of magnetic susceptibility has been successfully used for mineral exploration during the last decade. However, the unknown direction of magnetization has limited the use of this technique when significant remanence is present. We have developed a comprehensive methodology for solving this problem by examining two classes of approaches and have formulated a suite of methods of practical utility. The first class focuses on estimating total magnetization direction and then incorporating the resultant direction into an inversion algorithm that assumes a known direction. The second class focuses on direct inversion of the amplitude of the magnetic anomaly vector. Amplitude data depend weakly upon magnetization direction and are amenable to direct inversion for the magnitude of magnetization vector in 3D subsurface. Two sets of high-resolution aeromagnetic data acquired for diamond exploration in the Canadian Arctic are used to illustrate the methods’ usefulness.

Stepped inversion of magnetic data

Geophysics ◽  
2007 ◽  
Vol 72 (3) ◽  
pp. L21-L30 ◽  
Author(s):  
Soraya Lozada Tuma ◽  
Carlos Alberto Mendonça
Keyword(s):  
Magnetic Anomaly ◽  
Shape Function ◽  
Source Parameters ◽  
Real Data ◽  
Magnetic Data ◽  
Total Field ◽  
Magnetization Direction ◽  
Magnetization Intensity ◽  
Magnetic Inversion ◽  
Total Gradient

We present a three-step magnetic inversion procedure in which invariant quantities with respect to source parameters are inverted sequentially to give (1) shape cross section, (2) magnetization intensity, and (3) magnetization direction for a 2D (elongated) magnetic source. The quantity first inverted (called here the shape function) is obtained from the ratio of the gradient intensity of the total-field anomaly to the intensity of the anomalous vector field. For homogenous sources, the shape function is invariant with source magnetization and allows reconstruction of the source geometry by attributing an arbitrary magnetization to trial solutions. Once determined, the source shape is fixed and magnetization intensity is estimated by fitting the total gradient of the total-field anomaly (equivalent to the amplitude of the analytic signal of magnetic anomaly). Finally, the source shape and magnetization intensity are fixed and the magnetization direction is determined by fitting the magnetic anomaly. As suggested by numerical modeling and real data application, stepped inversion allows checking whether causative sources are homogeneous. This is possible because the shape function from inhomogeneous sources can be fitted by homogeneous models, but a model obtained in this way fits neither the total gradient of the magnetic anomaly nor the magnetic anomaly itself. Such a criterion seems effective in recognizing strongly inhomogeneous sources. Stepped inversion is tested with numerical experiments, and is used to model a magnetic anomaly from intrusive basic rocks from the Paraná Basin, Brazil.

Separation of regional and residual magnetic field data

Geophysics ◽  
1998 ◽  
Vol 63 (2) ◽  
pp. 431-439 ◽  
Author(s):  
Yaoguo Li ◽  
Douglas W. Oldenburg
Keyword(s):  
Magnetic Field ◽  
Field Data ◽  
Power Spectra ◽  
Magnetic Data ◽  
Inversion Algorithm ◽  
Large Area ◽  
Data Set ◽  
Regional Sources ◽  
Magnetic Inversion ◽  
Susceptibility Model

We present a method for separating regional and residual magnetic fields using a 3-D magnetic inversion algorithm. The separation is achieved by inverting the observed magnetic data from a large area to construct a regional susceptibility distribution. The magnetic field produced by the regional susceptibility model is then used as the regional field, and the residual data are obtained by simple subtraction. The advantages of this method of separation are that it introduces little distortion to the shape of the extracted anomaly and that it is not affected significantly by factors such as topography and the overlap of power spectra of regional and residual fields. The proposed method is tested using a synthetic example having varying relative positions between the local and regional sources and then using a field data set from Australia. Results show that the residual field extracted using this method enables good recovery of target susceptibility distribution from inversions.

Inversion of magnetic anomaly on rugged observation surface in the presence of strong remanent magnetization

Geophysics ◽  
2014 ◽  
Vol 79 (2) ◽  
pp. J11-J19 ◽  
Author(s):  
Shu-Ling Li ◽  
Yaoguo Li
Keyword(s):  
Remanent Magnetization ◽  
Country Rock ◽  
Magnetization Vector ◽  
Magnetic Data ◽  
Magnetization Direction ◽  
Equivalent Source ◽  
Total Magnetization ◽  
Amplitude Inversion ◽  
Amplitude Data ◽  
And Inversion

We study the inversion of magnetic data acquired over a rugged observation surface and where the buried source bodies have strong remanent magnetization that leads to unknown total magnetization directions. These factors pose significant challenges for processing and inversion of such data. To tackle the challenges from both a rugged observation surface and an unknown magnetization direction, we propose a strategy through the joint use of the equivalent source technique and 3D amplitude inversion to obtain 3D magnetization strength. We use equivalent source processing to calculate the amplitude data in the space domain because the use of the wavenumber-domain method is invalid due to large variations in the data elevation. We then carried out an amplitude inversion to generate a 3D subsurface distribution of the magnitude of the total magnetization vector. The results from a synthetic example and aeromagnetic data in Daye Mine in China showed that this approach is effective and images the magnetic units whose contact zones with the limestone country rock host the mineralization. The method is general and can be applied to a variety of cases with similar challenges.

Magnetic inversion to recover the subsurface block structures based on L1 norm and total variation regularization

2021 ◽  
Author(s):  
Mitsuru Utsugi
Keyword(s):  
Total Variation ◽  
Magnetic Data ◽  
Preconditioned Conjugate Gradient ◽  
Inversion Method ◽  
Total Variation Regularization ◽  
Subsurface Structure ◽  
L1 Norm ◽  
Magnetic Inversion ◽  
Original Definition ◽  
Primal Dual

Summary This paper presents a new sparse inversion method based on L1 norm regularization for 3D magnetic data. In isolation, L1 norm regularization yields model elements which are unconstrained by the input data to be exactly zero, leading to a sparse model with compact and focused structure. Here, we complement the L1 norm with a penalty minimizing total variation, the L1 norm of the model gradients; it is expected that the sharp boundaries of the subsurface structure are not compromised by incorporating this penalty. Although this penalty is widely used in the geophysical inversion studies, it is often replaced by an alternative quadratic penalty to ease solution of the penalized inversion problem; in this study, the original definition of the total variation, i.e., form of the L1 norm of the model gradients, is used. To solve the problem with this combined penalty of L1 norm and total variation, this study introduces alternative direction method of multipliers (ADMM), which is a primal-dual optimization algorithm that solves convex penalized problems based on the optimization of an augmented Lagrange function. To improve the computational efficiency of the algorithm to make this method applicable to large-scale magnetic inverse problems, this study applies matrix compression using the wavelet transform and the preconditioned conjugate gradient method. The inversion method is applied to both synthetic tests and real data, the synthetic tests demonstrate that, when subsurface structure is blocky, it can be reproduced almost perfectly.

Constrained 3D inversion of magnetic data with structural orientation and borehole lithology: A case study in the Macheng iron deposit, Hebei, China

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. B121-B133 ◽  
Author(s):  
Shida Sun ◽  
Chao Chen ◽  
Yiming Liu
Keyword(s):  
Remanent Magnetization ◽  
Magnetic Data ◽  
Iron Deposit ◽  
Reference Field ◽  
Total Field ◽  
Equivalent Source ◽  
Structural Orientation ◽  
Amplitude Data ◽  
Ore Bodies

We have developed a case study on the use of constrained inversion of magnetic data for recovering ore bodies quantitatively in the Macheng iron deposit, China. The inversion is constrained by the structural orientation and the borehole lithology in the presence of high magnetic susceptibility and strong remanent magnetization. Either the self-demagnetization effect caused by high susceptibility or strong remanent magnetization would lead to an unknown total magnetization direction. Here, we chose inversion of amplitude data that indicate low sensitivity to the direction of magnetization of the sources when constructing the underground model of effective susceptibility. To reduce the errors that arise when treating the total-field anomaly as the projection of an anomalous field vector in the direction of the geomagnetic reference field, we develop an equivalent source technique to calculate the amplitude data from the total-field anomaly. This equivalent source technique is based on the acquisition of the total-field anomaly, which uses the total-field intensity minus the magnitude of the reference field. We first design a synthetic model from a simplified real case to test the new approach, involving the amplitude data calculation and the constrained amplitude inversion. Then, we apply this approach to the real data. The results indicate that the structural orientation and borehole susceptibility bounds are compatible with each other and are able to improve the quality of the recovered model to obtain the distribution of ore bodies quantitatively and effectively.

scholarly journals A fast methodology for large-scale focusing inversion of gravity and magnetic data using the structured model matrix and the 2-D fast Fourier transform

2020 ◽  
Vol 223 (2) ◽  
pp. 1378-1397
Author(s):  
Rosemary A Renaut ◽  
Jarom D Hogue ◽  
Saeed Vatankhah ◽  
Shuang Liu
Keyword(s):  
Fourier Transform ◽  
Fast Fourier Transform ◽  
Linear Systems ◽  
Large Scale ◽  
Surface Measurement ◽  
Magnetic Data ◽  
Uniform Grid ◽  
Data Sets ◽  
Inversion Algorithm ◽  
Data Set

SUMMARY We discuss the focusing inversion of potential field data for the recovery of sparse subsurface structures from surface measurement data on a uniform grid. For the uniform grid, the model sensitivity matrices have a block Toeplitz Toeplitz block structure for each block of columns related to a fixed depth layer of the subsurface. Then, all forward operations with the sensitivity matrix, or its transpose, are performed using the 2-D fast Fourier transform. Simulations are provided to show that the implementation of the focusing inversion algorithm using the fast Fourier transform is efficient, and that the algorithm can be realized on standard desktop computers with sufficient memory for storage of volumes up to size n ≈ 106. The linear systems of equations arising in the focusing inversion algorithm are solved using either Golub–Kahan bidiagonalization or randomized singular value decomposition algorithms. These two algorithms are contrasted for their efficiency when used to solve large-scale problems with respect to the sizes of the projected subspaces adopted for the solutions of the linear systems. The results confirm earlier studies that the randomized algorithms are to be preferred for the inversion of gravity data, and for data sets of size m it is sufficient to use projected spaces of size approximately m/8. For the inversion of magnetic data sets, we show that it is more efficient to use the Golub–Kahan bidiagonalization, and that it is again sufficient to use projected spaces of size approximately m/8. Simulations support the presented conclusions and are verified for the inversion of a magnetic data set obtained over the Wuskwatim Lake region in Manitoba, Canada.

Magnetization vector imaging for borehole magnetic data based on magnitude magnetic anomaly

Geophysics ◽  
2013 ◽  
Vol 78 (6) ◽  
pp. D429-D444 ◽  
Author(s):  
Shuang Liu ◽  
Xiangyun Hu ◽  
Tianyou Liu ◽  
Jie Feng ◽  
Wenli Gao ◽  
...  
Keyword(s):  
Conjugate Gradient Method ◽  
Conjugate Gradient ◽  
Gradient Method ◽  
Magnetization Vector ◽  
Real Data ◽  
Magnetic Anomalies ◽  
Magnetic Data ◽  
Preconditioned Conjugate Gradient ◽  
Magnetization Direction ◽  
Magnetization Intensity

Remanent magnetization and self-demagnetization change the magnitude and direction of the magnetization vector, which complicates the interpretation of magnetic data. To deal with this problem, we evaluated a method for inverting the distributions of 2D magnetization vector or effective susceptibility using 3C borehole magnetic data. The basis for this method is the fact that 2D magnitude magnetic anomalies are not sensitive to the magnetization direction. We calculated magnitude anomalies from the measured borehole magnetic data in a spatial domain. The vector distributions of magnetization were inverted methodically in two steps. The distributions of magnetization magnitude were initially solved based on magnitude magnetic anomalies using the preconditioned conjugate gradient method. The preconditioner determined by the distances between the cells and the borehole observation points greatly improved the quality of the magnetization magnitude imaging. With the calculated magnetization magnitude, the distributions of magnetization direction were computed by fitting the component anomalies secondly using the conjugate gradient method. The two-step approach made full use of the amplitude and phase anomalies of the borehole magnetic data. We studied the influence of remanence and demagnetization based on the recovered magnetization intensity and direction distributions. Finally, we tested our method using synthetic and real data from scenarios that involved high susceptibility and complicated remanence, and all tests returned favorable results.

Inversion of airborne geophysics over the DO-27/DO-18 kimberlites — Part 1: Potential fields

2017 ◽  
Vol 5 (3) ◽  
pp. T299-T311 ◽  
Author(s):  
Sarah G. R. Devriese ◽  
Kristofer Davis ◽  
Douglas W. Oldenburg
Keyword(s):  
Magnetic Data ◽  
Gravity Gradiometry ◽  
3D Inversion ◽  
Potential Field Data ◽  
Electromagnetic Data ◽  
Magnetic Inversion ◽  
Geophysical Surveys ◽  
Airborne Geophysics ◽  
Joint Interpretation ◽  
Geologic Model

The Tli Kwi Cho (TKC) kimberlite complex contains two pipes, called DO-27 and DO-18, which were discovered during the Canadian diamond exploration rush in the 1990s. The complex has been used as a testbed for ground and airborne geophysics, and an abundance of data currently exist over the area. We have evaluated the historical and geologic background of the complex, the physical properties of interest for kimberlite exploration, and the geophysical surveys. We have carried out 3D inversion and joint interpretation of the potential field data. The magnetic data indicate high susceptibility at DO-18, and the magnetic inversion maps the horizontal extent of the pipe. DO-27 is more complicated. The northern part is highly magnetic and is contaminated with remanent magnetization; other parts of DO-27 have a low susceptibility. Low densities, obtained from the gravity and gravity gradiometry data, map the horizontal extents of DO-27 and DO-18. We combine the 3D density contrast and susceptibility models into a single geologic model that identifies three distinct kimberlite rock units that agree with drilling data. In further research, our density and magnetic susceptibility models are combined with information from electromagnetic data to provide a multigeophysical interpretation of the TKC kimberlite complex.

Sign in / Sign up

Export Citation Format

Share Document