Probabilistic inversion of airborne electromagnetic data for basement conductors

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. E389-E400 ◽  
Author(s):  
Juerg Hauser ◽  
James Gunning ◽  
David Annetts

Probabilistic inversion of airborne electromagnetic data is often approximated by a layered earth using a computationally efficient 1D kernel. If the underlying framework accounts for prior beliefs on spatial correlation, the inversion will be able to recover spatially coherent interfaces and associated uncertainties. Greenfield exploration using airborne electromagnetic data, however, often seeks to identify discrete economical targets. In mature exploration provinces, such bodies are frequently obscured by thick, conductive regolith, and the response of such economic basement conductors presents a challenge to any layered earth inversion. A well-known computationally efficient way to approximate the response of a basement conductor is to use a thin plate. Here we have extended a Bayesian parametric bootstrap approach, so that the basement of a spatially varying layered earth can contain a thin plate. The resulting Bayesian framework allowed for the inversion of basement conductors and associated uncertainties, but more importantly, the use of model selection concepts to determine if the data supports a basement conductor model or not. Recovered maps of basement conductor probabilities show the expected patterns in uncertainty; for example, a decrease in target probability with increasing depth. Such maps of target probabilities generated using the thin plate approximation are a potentially valuable source of information for the planning of exploration activity, such as the targeting of drillholes to confirm the existence of a discrete conductor in a greenfield exploration scenario. We have used a field data set from northwest Queensland, Australia, to illustrate how the approach allowed inversion for a basement conductor and related uncertainties in a spatially variable layered earth, using the information from multiple survey lines and prior beliefs of geology.

Geophysics ◽  
2015 ◽  
Vol 80 (2) ◽  
pp. E135-E146 ◽  
Author(s):  
Juerg Hauser ◽  
James Gunning ◽  
David Annetts

Probabilistic 1D inversions of airborne electromagnetic (AEM) surveys allow an exhaustive search of model space for each station, but they often assume that there is no spatial correlation between neighboring stations. This can result in abrupt transverse model discontinuities when attempting to construct a 3D model. In contrast to this, fully spatially regularized deterministic inversions can take spatial correlation between 1D models into account, but they do not explore the model space sufficiently to be able to evaluate model robustness. The Bayesian parametric bootstrap (BPB) approach that we developed is a practical compromise between computationally expensive exhaustive search techniques and computationally efficient deterministic inversions. Using a 1D kernel, we inverted for the interfaces, layer properties, and related uncertainties, taking lateral spatial correlations and additional prior information into account. Numerical examples revealed that a BPB technique was likely to explore the model space sufficiently for nonpathological situations. Using a subset of a large AEM survey collected in northwest Australia for aquifer mapping, we show how the BPB approach can be used to produce a spatially coherent map of the base of the Broome sandstone aquifer. The recovered uncertainties, which are likely to be one of the main sources of uncertainty in any groundwater model, exhibited the well-known increase in uncertainty of a depth to interface with increasing depth to the interface.


2020 ◽  
Vol 224 (1) ◽  
pp. 543-557
Author(s):  
Thomas M Hansen

SUMMARY Probabilistic inversion methods, typically based on Markov chain Monte Carlo, exist that allow exploring the full uncertainty of geophysical inverse problems. The use of such methods is though limited by significant computational demands, and non-trivial analysis of the obtained set of dependent models. Here, a novel approach, for sampling the posterior distribution is suggested based on using pre-calculated lookup tables with the extended rejection sampler. The method is (1) fast, (2) generates independent realizations of the posterior, and (3) does not get stuck in local minima. It can be applied to any inverse problem (and sample an approximate posterior distribution) but is most promising applied to problems with informed prior information and/or localized inverse problems. The method is tested on the inversion of airborne electromagnetic data and shows an increase in the computational efficiency of many orders of magnitude as compared to using the extended Metropolis algorithm.


Geophysics ◽  
2003 ◽  
Vol 68 (4) ◽  
pp. 1211-1223 ◽  
Author(s):  
Haoping Huang ◽  
Douglas C. Fraser

Inversion of airborne electromagnetic (EM) data for a layered earth has been commonly performed under the assumption that the magnetic permeability of the layers is the same as that of free space. The resistivity inverted from helicopter EM data in this way is not reliable in highly magnetic areas because magnetic polarization currents occur in addition to conduction currents, causing the inverted resistivity to be erroneously high. A new algorithm for inverting for the resistivity, magnetic permeability, and thickness of a layered model has been developed for a magnetic conductive layered earth. It is based on traditional inversion methodologies for solving nonlinear inverse problems and minimizes an objective function subject to fitting the data in a least‐squares sense. Studies using synthetic helicopter EM data indicate that the inversion technique is reasonably dependable and provides fast convergence. When six synthetic in‐phase and quadrature data from three frequencies are used, the model parameters for two‐ and three‐layer models are estimated to within a few percent of their true values after several iterations. The analysis of partial derivatives with respect to the model parameters contributes to a better understanding of the relative importance of the model parameters and the reliability of their determination. The inversion algorithm is tested on field data obtained with a Dighem helicopter EM system at Mt. Milligan, British Columbia, Canada. The output magnetic susceptibility‐depth section compares favorably with that of Zhang and Oldenburg who inverted for the susceptibility on the assumption that the resistivity distribution was known.


Geophysics ◽  
2015 ◽  
Vol 80 (6) ◽  
pp. K25-K36 ◽  
Author(s):  
Michael S. McMillan ◽  
Christoph Schwarzbach ◽  
Eldad Haber ◽  
Douglas W. Oldenburg

1990 ◽  
Author(s):  
Clyde Bergeron ◽  
Terrence L. Morris ◽  
Juliette W. Ioup

Geophysics ◽  
2002 ◽  
Vol 67 (2) ◽  
pp. 492-500 ◽  
Author(s):  
James E. Reid ◽  
James C. Macnae

When a confined conductive target embedded in a conductive host is energized by an electromagnetic (EM) source, current flow in the target comes from both direct induction of vortex currents and current channeling. At the resistive limit, a modified magnetometric resistivity integral equation method can be used to rapidly model the current channeling component of the response of a thin-plate target energized by an airborne EM transmitter. For towed-bird transmitter–receiver geometries, the airborne EM anomalies of near-surface, weakly conductive features of large strike extent may be almost entirely attributable to current channeling. However, many targets in contact with a conductive host respond both inductively and galvanically to an airborne EM system. In such cases, the total resistive-limit response of the target is complicated and is not the superposition of the purely inductive and purely galvanic resistive-limit profiles. Numerical model experiments demonstrate that while current channeling increases the width of the resistive-limit airborne EM anomaly of a wide horizontal plate target, it does not necessarily increase the peak anomaly amplitude.


Sign in / Sign up

Export Citation Format

Share Document