3D Inversion of Magnetic Amplitude Data with Sparseness Constraint

Author(s):  
Mohammad Rezaie
Geophysics ◽  
2010 ◽  
Vol 75 (1) ◽  
pp. L1-L11 ◽  
Author(s):  
Yaoguo Li ◽  
Sarah E. Shearer ◽  
Matthew M. Haney ◽  
Neal Dannemiller

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.


Geophysics ◽  
2018 ◽  
Vol 83 (5) ◽  
pp. J75-J84 ◽  
Author(s):  
Camriel Coleman ◽  
Yaoguo Li

Three-dimensional inversion plays an important role in the quantitative interpretation of magnetic data in exploration problems, and magnetic amplitude data can be an effective tool in cases in which remanently magnetized materials are present. Because amplitude data are typically calculated from total-field anomaly data, the error levels must be characterized for inversions. Lack of knowledge of the error in amplitude data hinders the ability to properly estimate the data misfit associated with an inverse model and, therefore, the selection of the appropriate regularization parameter for a final model. To overcome these challenges, we have investigated the propagation of errors from total-field anomaly to amplitude data. Using parametric bootstrapping, we find that the standard deviation of the noise in amplitude data is approximately equal to that of the noise in total-field anomaly data when the amplitude data are derived from the conversion of total-field data to three orthogonal components. We then illustrate how the equivalent source method can be used to estimate the error in total-field anomaly data when needed. The obtained noise estimate can be applied to amplitude inversion to recover an optimal inverse model by applying the discrepancy principle. We test this method on synthetic and field data and determine its effectiveness.


2019 ◽  
Author(s):  
Zelin Li* ◽  
Changli Yao

Geophysics ◽  
2015 ◽  
Vol 80 (2) ◽  
pp. B13-B22 ◽  
Author(s):  
Marcelo Leão-Santos ◽  
Yaoguo Li ◽  
Roberto Moraes

Strong hydrothermal alteration modifies rock physical properties in iron oxide-copper-gold deposits (IOCGs) and may result in characteristic signatures detectable in geophysical surveys. Magnetic data are commonly used in characterizing orebodies, and 3D inversions are often used to assist in interpretations. In areas with strong remanence and self-demagnetization, the total magnetization can have directions different from the inducing field direction. This deviation precludes the use of traditional inversion methods. Magnetic amplitude inversion offers one solution to this challenge because the amplitude data are weakly dependent on the magnetization direction. In addition, the low magnetic latitude also imposes difficulty in amplitude data calculation due to the instability in the component conversion in the wavenumber domain. To formulate a practical approach, we present a case study on applying the magnetic amplitude inversion to the Furnas southeast IOCG deposit at the low magnetic latitude in Carajás Mineral Province, Brazil, and demonstrate that the approach can reliably recover an interpretable distribution of effective magnetic susceptibility and identify massive magnetite from hydrothermal alterations associated with the high-grade ore.


2020 ◽  
Vol 68 (2) ◽  
pp. 365-375
Author(s):  
Zelin Li ◽  
Changli Yao

Geophysics ◽  
2021 ◽  
Vol 86 (1) ◽  
pp. J1-J11
Author(s):  
Marlon C. Hidalgo-Gato ◽  
Valéria C. F. Barbosa ◽  
Vanderlei C. Oliveira

We have developed an inversion method to recover the depth and the total magnetization intensity of the basement under a sedimentary basin using the amplitude of the magnetic anomaly vector (amplitude data). Because the amplitude data are weakly dependent on the magnetization direction, our method is suitable for interpreting areas with remanent magnetization. Our method assumes constant magnetized basement rocks overlain by nonmagnetic sediments. To overcome the inherent ambiguity of potential field data, we assume knowledge of the average depth of the basement and use it as a constraint to regularize the inversion. A sensitivity analysis with synthetic data shows the weak dependency of the magnetic amplitude inversion on the magnetization direction. Different combinations of magnetization directions recover the interface separating sediments from basement rocks. Test on field data over the Foz do Amazonas Basin, Brazil, recovers the shape of the basement relief without any knowledge about the magnetization intensity and direction. The estimated basement relief reveals a smooth basement framework with basement highs in the central part of the area. In a regional-scale perspective, the deeper and constant estimated basement relief at the northernmost limit of the area may suggest changing in crustal domains from a hyperextended continental crust to homogeneous oceanic crust.


Geophysics ◽  
2020 ◽  
Vol 85 (3) ◽  
pp. B63-B75
Author(s):  
Kaijun Xu ◽  
Yaoguo Li

We present a case study on imaging volcanic units in gas exploration by constraining magnetic amplitude inversions using magnetotelluric (MT) sounding data at sparse locations. Magnetic data can be effective in mapping volcanic units because they have remanent magnetization and significant susceptibility contrast with surrounding rocks. Although magnetic data can identify the lateral distribution of volcanic units, they often have difficulties in defining the depth extent. For this reason, additional structural constraints from other geophysical methods can often help improve the vertical resolution. Among the independent geophysical methods, MT data can provide the needed structural information at a low cost. We have investigated an approach to combine a set of sparse MT soundings with magnetic amplitude data to image the distribution of volcanics in a basin environment. We first use a blocky 1D MT inversion based on Ekblom norm to obtain the structural constraint, and then we perform a constrained 3D magnetic amplitude inversion to recover the distribution of effective susceptibility by incorporating the structural information from MT soundings. We determine that even a small number of MT stations (e.g., 20) in a [Formula: see text] area is sufficient to drastically improve the magnetic amplitude inversion. Our results indicate that magnetic amplitude inversion with structural constraint from MT soundings form a practical and cost-effective means to map the lateral and vertical distribution of volcanics.


2015 ◽  
Author(s):  
Zelin Li ◽  
Changli Yao ◽  
Yuanman Zheng ◽  
Xiaohong Meng

Entropy ◽  
2021 ◽  
Vol 23 (6) ◽  
pp. 676
Author(s):  
Dimitrios Z. Politis ◽  
Stelios M. Potirakis ◽  
Yiannis F. Contoyiannis ◽  
Sagardweep Biswas ◽  
Sudipta Sasmal ◽  
...  

In this work we present the statistical and criticality analysis of the very low frequency (VLF) sub-ionospheric propagation data recorded by a VLF/LF radio receiver which has recently been established at the University of West Attica in Athens (Greece). We investigate a very recent, strong (M6.9), and shallow earthquake (EQ) that occurred on 30 October 2020, very close to the northern coast of the island of Samos (Greece). We focus on the reception data from two VLF transmitters, located in Turkey and Israel, on the basis that the EQ’s epicenter was located within or very close to the 5th Fresnel zone, respectively, of the corresponding sub-ionospheric propagation path. Firstly, we employed in our study the conventional analyses known as the nighttime fluctuation method (NFM) and the terminator time method (TTM), aiming to reveal any statistical anomalies prior to the EQ’s occurrence. These analyses revealed statistical anomalies in the studied sub-ionospheric propagation paths within ~2 weeks and a few days before the EQ’s occurrence. Secondly, we performed criticality analysis using two well-established complex systems’ time series analysis methods—the natural time (NT) analysis method, and the method of critical fluctuations (MCF). The NT analysis method was applied to the VLF propagation quantities of the NFM, revealing criticality indications over a period of ~2 weeks prior to the Samos EQ, whereas MCF was applied to the raw receiver amplitude data, uncovering the time excerpts of the analyzed time series that present criticality which were closest before the Samos EQ. Interestingly, power-law indications were also found shortly after the EQ’s occurrence. However, it is shown that these do not correspond to criticality related to EQ preparation processes. Finally, it is noted that no other complex space-sourced or geophysical phenomenon that could disturb the lower ionosphere did occur during the studied time period or close after, corroborating the view that our results prior to the Samos EQ are likely related to this mainshock.


Preview ◽  
2011 ◽  
Vol 2011 (152) ◽  
pp. 29-34 ◽  
Author(s):  
Glenn Wilson ◽  
Martin Čuma ◽  
Michael S. Zhdanov

Sign in / Sign up

Export Citation Format

Share Document