Transdimensional Bayesian inversion of time-domain airborne EM data

To reduce the dependence of EM inversion on the choice of initial model and to obtain the global minimum, we apply transdimensional Bayesian inversion to time-domain airborne electromagnetic data. The transdimensional Bayesian inversion uses the Monte Carlo method to search the model space and yields models that simultaneously satisfy the acceptance probability and data fitting requirements. Finally, we obtain the probability distribution and uncertainty of the model parameters as well as the maximum probability. Because it is difficult to know the height of the transmitting source during flight, we consider a fixed and a variable flight height. Furthermore, we introduce weights into the prior probability density function of the resistivity and adjust the constraint strength in the inversion model by changing the weighing coefficients. This effectively solves the problem of unsatisfactory inversion results in the middle high-resistivity layer. We validate the proposed method by inverting synthetic data with 3% Gaussian noise and field survey data.

[1]  Changchun Yin,et al.  3D inversion for multipulse airborne transient electromagnetic data , 2016 .

[2]  E. Auken,et al.  A comparison of helicopter-borne electromagnetics in frequency- and time-domain at the Cuxhaven valley in Northern Germany , 2009 .

[3]  Greg Hodges,et al.  Simulated annealing for airborne EM inversion , 2007 .

[4]  Stan E. Dosso,et al.  Non-linearity in Bayesian 1-D magnetotelluric inversion , 2011 .

[5]  James Macnae,et al.  Evaluating EM waveforms by singular‐value decomposition of exponential basis functions , 1998 .

[6]  P. Green Reversible jump Markov chain Monte Carlo computation and Bayesian model determination , 1995 .

[7]  Hiroshi Amano,et al.  2.5‐D inversion of frequency‐domain electromagnetic data generated by a grounded‐wire source , 2002 .

[8]  Anandaroop Ray,et al.  Bayesian inversion of marine CSEM data from the Scarborough gas field using a transdimensional 2-D parametrization , 2014 .

[9]  Alberto Malinverno,et al.  A Bayesian criterion for simplicity in inverse problem parametrization , 2000 .

[10]  Yin Chang,et al.  The full-time electromagnetic modeling for time-domain airborne electromagnetic systems , 2013 .

[11]  Misac N. Nabighian,et al.  Electromagnetic Methods in Applied Geophysics: Voume 1, Theory , 1988 .

[12]  B. Minsley A trans-dimensional Bayesian Markov chain Monte Carlo algorithm for model assessment using frequency-domain electromagnetic data , 2011 .

[13]  Mrinal K. Sen,et al.  2-D resistivity inversion using spline parameterization and simulated annealing , 1996 .

[14]  A. Malinverno Parsimonious Bayesian Markov chain Monte Carlo inversion in a nonlinear geophysical problem , 2002 .

[15]  Greg Hodges,et al.  MULTIPULSE – high resolution and high power in one TDEM system , 2015 .

[16]  Anandaroop Ray,et al.  Bayesian inversion of marine CSEM data with a trans‐dimensional self parametrizing algorithm , 2012 .

[17]  M. Sambridge,et al.  Seismic tomography with the reversible jump algorithm , 2009 .

[18]  Rhys Hawkins,et al.  Trans-dimensional Bayesian inversion of airborne electromagnetic data for 2D conductivity profiles , 2017 .

[19]  Jun Lin,et al.  PC-based artificial neural network inversion for airborne time-domain electromagnetic data , 2012, Applied Geophysics.

[20]  Yin Changchun,et al.  Weighted Laterally-Constrained Inversion of Time-Domain Airborne Electromagnetic Data , 2016 .

[21]  A. Malinverno,et al.  Receiver function inversion by trans‐dimensional Monte Carlo sampling , 2010 .

[22]  Whitney Trainor-Guitton,et al.  Stochastic inversion for electromagnetic geophysics: Practical challenges and improving convergence efficiency , 2011 .

[23]  X Shi MULTISCALE GENETIC ALGORITHM AND ITS APPLICATION IN MAGNETOTELLURIC SOUNDING DATA INVERSION , 2000 .