scholarly journals 3D magnetic amplitude inversion in the presence of self-demagnetization and remanent magnetization

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. J69-J82 ◽  
Author(s):  
Boxin Zuo ◽  
Xiangyun Hu ◽  
Yi Cai ◽  
Shuang Liu
Keyword(s):  
Remanent Magnetization ◽  
Large Scale ◽  
Magnetic Data ◽  
Inversion Problem ◽  
Aeromagnetic Data ◽  
Inversion Algorithm ◽  
Rapid Convergence ◽  
Magnetization Direction ◽  
Amplitude Inversion ◽  
Nonlinear Conjugate Gradient

We have developed a general 3D amplitude inversion algorithm for magnetic data in the presence of self-demagnetization and remanent magnetization. The algorithm uses a nonlinear conjugate gradient (NLCG) scheme to invert the amplitude of the magnetic anomaly vector within a partial differential equation framework. Three quantities— the amplitude of the anomalous magnetic field, the analytic signal, and the normalized source strength, defined as the amplitudes of magnetic data that are weakly dependent on the magnetization direction — are inverted to recover the 3D distribution of the subsurface magnetic susceptibility. Numerical experiments indicate that our NLCG amplitude inversion algorithm has a rapid convergence rate that provides a reasonable inversion solution in the absence of knowing the total magnetization direction. High-resolution aeromagnetic data collected from the Pea Ridge iron oxide-apatite-rare earth element deposit, southeast Missouri, USA, are used to illustrate the efficacy of our amplitude inversion algorithm. This algorithm is generally applicable for tackling the large-scale inversion problem in the presence of self-demagnetization and remanent magnetization.

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.

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.

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.

Practical considerations in the use of edge detectors for geologic mapping using magnetic data

Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. J1-J8 ◽  
Author(s):  
Mark Pilkington ◽  
Victoria Tschirhart
Keyword(s):  
Magnetic Field ◽  
Detection Methods ◽  
Magnetic Data ◽  
Horizontal Gradient ◽  
Recent Effort ◽  
Aeromagnetic Data ◽  
Geologic Mapping ◽  
Magnetization Direction ◽  
Model Studies ◽  
Edge Detectors

Locating the edges of magnetized sources provides a fundamental tool in the geologic interpretation of magnetic field data. Much recent effort has been expended on developing improvements to existing edge-detection methods, resulting in purported increases in accuracy and continuity along edges, reduction of noise effects, and limiting the influences of variable depth to source, magnetization direction, and source dip. These endeavors are valuable and provide interpreters with a wider range of tools to carry out geologic interpretations of aeromagnetic data. Nevertheless, survey parameters such as flight height and line spacing impose limits on the quality of edge locations that can be achieved. Using model studies, we quantify the effects that source size, depth, and interference between sources have on calculated edge locations. Based on the known behavior of established edge detectors, we found that many of the newer approaches offer limited advantages over older methods. Consequently, we studied an example of field mapping of geologic contacts in the Canadian Shield, supported by aeromagnetic data, using calculation of a standard edge detector: the horizontal gradient magnitude of the total magnetic field or TF-hgm. Calculated edge locations estimated from this method appear sufficiently accurate and continuous to provide a solid basis on which the mapping campaign was based and executed successfully.

A New Method of Cross-Correlation by Magnetic Dipole for Estimating Magnetization Direction under the Influence of Remanent Magnetization

2014 ◽  
Vol 644-650 ◽  
pp. 3459-3462 ◽  
Author(s):  
Lei Shi ◽  
Liang Hui Guo ◽  
Feng Yi Guo
Keyword(s):  
Correlation Coefficient ◽  
Magnetic Dipole ◽  
Remanent Magnetization ◽  
Cross Correlation ◽  
New Method ◽  
Magnetic Data ◽  
Dipole Source ◽  
Total Field ◽  
Magnetization Direction ◽  
Total Magnetization

Processing and interpretation of magnetic data usually require information of total magnetization direction. However, under the effects of remanent magnetization, total magnetization direction is different from induced magnetization direction, which makes data processing and interpretation complexity. In this paper, we present a new method by cross-correlation of magnetic dipole source for determination of magnetization direction from relatively isolated and approximate equiaxial-shape magnetic total field anomaly. This method calculates cross-correlation coefficient between observed magnetic total field anomaly and theoretical magnetic total field anomaly caused by a magnetic dipole source, by using a set of varying parameters of positions and total magnetization direction of dipole source for trial and error. The corresponding magnetization direction of maximum correlation coefficient is regarded as estimated total magnetization direction. Test on synthetic data indicates that this method reliably and effectively estimates the magnetization direction from relatively isolated and approximate equiaxial-shape magnetic total field anomaly.

Estimating the Magnetization Direction of Sources through Correlation between Reduced-to-Pole Anomaly and Normalized Source Strength

2014 ◽  
Vol 644-650 ◽  
pp. 3793-3796
Author(s):  
Liang Hui Guo ◽  
Rui Gao ◽  
Guo Li Zhang
Keyword(s):  
Geomagnetic Field ◽  
Remanent Magnetization ◽  
Cross Correlation ◽  
Synthetic Data ◽  
Field Direction ◽  
Magnetic Data ◽  
Source Strength ◽  
Magnetization Direction ◽  
New Approach ◽  
Total Magnetization

Under the effects of remanent magnetization, total magnetization direction is different from geomagnetic field direction, which makes magnetic data processing and interpretation complexity. In this paper, we present a new approach for estimating the total magnetization direction of sources via cross-correlation between the reduced-to-pole anomaly and the normalized source strength (who is less sensitive to remanent magnetization). The geomagnetic field direction is used to calculated the normalized source strength, while various assumed total magnetization directions are used to calculated the RTP anomalies. The maximum correlation between the RTP anomalies and the normalized corresponds to the estimated total magnetization direction. Test on synthetic data showed that the new approach is simple and effective.

Mitigating remanent magnetization effects in magnetic data using the normalized source strength

Geophysics ◽  
2013 ◽  
Vol 78 (3) ◽  
pp. J25-J32 ◽  
Author(s):  
Mark Pilkington ◽  
Majid Beiki
Keyword(s):  
Magnetic Field ◽  
Remanent Magnetization ◽  
Field Data ◽  
Magnetic Data ◽  
Magnetization Distribution ◽  
Source Strength ◽  
Inversion Algorithm ◽  
Data Set ◽  
Magnetic Field Data ◽  
The Magnetic Field

We have developed an approach for the interpretation of magnetic field data that can be used when measured anomalies are affected by significant remanent magnetization components. The method deals with remanent effects by using the normalized source strength (NSS), a quantity calculated from the eigenvectors of the magnetic gradient tensor. The NSS is minimally affected by the direction of remanent magnetization present and compares well with other transformations of the magnetic field that are used for the same purpose. It therefore offers a way of inverting magnetic data containing the effects of remanent magnetizations, particularly when these are unknown and are possibly varying within a given data set. We use a standard 3D inversion algorithm to invert NSS data from an area where varying remanence directions are apparent, resulting in a more reliable image of the subsurface magnetization distribution than possible using the observed magnetic field data directly.

scholarly journals Seismic reflectivity inversion using an L1-norm basis-pursuit method and GPU parallelisation

2020 ◽  
Author(s):  
Ruo Wang ◽  
Yanghua Wang ◽  
Ying Rao
Keyword(s):  
Large Scale ◽  
Seismic Inversion ◽  
Computational Time ◽  
Basis Pursuit ◽  
Inversion Problem ◽  
Inversion Algorithm ◽  
L1 Norm ◽  
The Matrix ◽  
Seismic Reflectivity ◽  
Reflectivity Inversion

Abstract Seismic reflectivity inversion problem can be formulated using a basis-pursuit method, aiming to generate a sparse reflectivity series of the subsurface media. In the basis-pursuit method, the reflectivity series is composed by large amounts of even and odd dipoles, thus the size of the seismic response matrix is huge and the matrix operations involved in seismic inversion are very time-consuming. In order to accelerate the matrix computation, a basis-pursuit method-based seismic inversion algorithm is implemented on Graphics Processing Unit (GPU). In the basis-persuit inversion algorithm, the problem is imposed with a L1-norm model constraint for sparsity, and this L1-norm basis-pursuit inversion problem is reformulated using a linear programming method. The core problems in the inversion are large-scale linear systems, which are resolved by a parallelised conjugate gradient method. The performance of this fully parallelised implementation is evaluated and compared to the conventional serial coding. Specifically, the investigation using several field seismic data sets with different sizes indicates that GPU-based parallelisation can significantly reduce the computational time with an overall factor up to 145. This efficiency improvement demonstrates a great potential of the basis-pursuit inversion method in practical application to large-scale seismic reflectivity inversion problems.

Application of magnetic amplitude inversion in exploration for volcanic units in a basin environment

Geophysics ◽  
2012 ◽  
Vol 77 (5) ◽  
pp. B219-B225 ◽  
Author(s):  
Yaoguo Li ◽  
Zhanxiang He ◽  
Yunxiang Liu
Keyword(s):  
Remanent Magnetization ◽  
Volcanic Rocks ◽  
Wiener Filter ◽  
Inversion Method ◽  
Magnetization Direction ◽  
Amplitude Inversion ◽  
Final Interpretation ◽  
Local Geology ◽  
Magnetic Amplitude

We have performed a case study on the use of magnetic amplitude inversion in imaging volcanic rocks that are buried in a sedimentary basin and have strong remanent magnetization. The application arises in exploration for natural gas hosted in volcanic rocks in basin environments. The weak anomaly associated with the volcanic units and the presence of remanent magnetization therein pose major challenges in the magnetic interpretation. We first use a Wiener filter based on an ensemble analysis to examine the depth characteristics of the anomalies, and then we use an inversion-based signal separation to extract the anomaly for final interpretation. We apply an amplitude inversion method to recover the distribution of effective susceptibility in the absence of the knowledge about the total magnetization direction. The result effectively identifies the volcanic rock units at large depths and the imaged distribution of these units is consistent with information from drillholes and local geology.

Identifying remanent magnetization effects in magnetic data

Geophysics ◽  
1993 ◽  
Vol 58 (5) ◽  
pp. 653-659 ◽  
Author(s):  
Walter R. Roest ◽  
Mark Pilkington
Keyword(s):  
Magnetic Anomaly ◽  
Remanent Magnetization ◽  
Close Agreement ◽  
Magnetic Anomalies ◽  
Analytic Signal ◽  
Magnetic Data ◽  
Horizontal Gradient ◽  
Impact Structure ◽  
Magnetization Direction ◽  
Additional Constraints

Remanent magnetization can have a significant influence on the shape of magnetic anomalies in areas that are generally characterized by induced magnetization. Since modeling of magnetic anomalies is nonunique, additional constraints on the direction of magnetization are useful. A method is proposed here to study the possible contribution of remanent magnetization to a particular anomaly, by comparing two functions that are calculated directly from the observations: (1) the amplitude of the analytic signal, and (2) the horizontal gradient of pseudogravity. From the amplitude and relative position of maxima in these derived quantities, we infer the deviation of the magnetization direction from that of the ambient field. The approach is applied to the magnetic anomaly in the center of the Manicouagan impact structure (Canada). Our results, based only on the magnetic anomaly observations, are in close agreement with constraints on the direction of remanent magnetization from rock samples.

Sign in / Sign up

Export Citation Format

Share Document