1887
1st Australasian Exploration Geoscience Conference – Exploration Innovation Integration
  • ISSN: 2202-0586
  • E-ISSN:

Abstract

We have developed an algorithm and released open-source code for the 1D inversion of magnetotelluric data. The algorithm uses trans-dimensional Markov chain Monte Carlo techniques to solve for a probabilistic conductivity-depth model.

The inversion of each station employs multiple Markov Chains in parallel to generate an ensemble of millions of conductivity models that adequately fit the data given the assigned noise levels. The trans-dimensional aspect of the inversion means that the number of layers in the conductivity model is solved for rather than being predetermined and kept fixed. Each Markov chain increases and decrease the number of layers in the model and the depths of the interfaces as it samples.

Once the ensemble of models is generated, its statistics are analysed to assess the posterior probability distribution of the conductivity at any particular depth, as well as the number of layers and the depths of the interfaces. This stochastic approach gives a thorough exploration of model space and a more robust estimation of uncertainty than deterministic methods allow.

The method’s application to cover thickness estimation is discussed with synthetic and real examples. Inversion of complex impedance tensor and also derived apparent resistivity/phase data are both demonstrated. It is found that the more pronounced layer boundaries allow more straightforward interpretation of cover thickness than that from deterministic smooth model inversions. It is concluded that thickness estimates compare favourably with borehole stratigraphic logs in most cases, and that the method is a useful addition to a range of cover thickness estimation tools.

Loading

Article metrics loading...

/content/journals/10.1071/ASEG2018abT5_1F
2018-12-01
2026-01-20
Loading full text...

Full text loading...

References

  1. Bodin, T. and Sambridge, M., 2009. Seismic tomography with the reversible jump algorithm. Geophysical Journal International178(3), 1411-1436.
  2. Bodin, T., Sambridge, M., Tkalcic, H., Arroucau, P., Gallagher, K. and Rawlinson, N., 2012. Transdimensional inversion of receiver functions and surface wave dispersion. Journal of Geophysical Research117(B02301). doi: 10.1029/2011JB008560.
  3. Brodie, R. C. and Sambridge, M., 2012. Transdimensional Monte Carlo inversion of AEM data. ASEG Extended Abstracts2012(1), 4.
  4. Constable, S. C., Parker, R. L. and Constable, C. G., 1987. Occam’s inversion; a practical algorithm for generating smooth models from electromagnetic sounding data. Geophysics52(3), 289-300.
  5. Hawkins, R., 2013. iEarth web page for rj-McMC. Online: http://www.iearth.org.au/codes/rj-MCMC.
  6. Malinverno, A., 2002. Parsimonious Bayesian Markov chain Monte Carlo inversion in a nonlinear geophysical problem. Geophysical Journal International151(3), 675-688.
  7. Minsley, B. J., 2011. A trans-dimensional Bayesian Markov chain Monte Carlo algorithm for model assessment using frequency-domain electromagnetic data. Geophysical Journal International187(1), 252-272.
  8. Ray, A. and Key, K., 2012. Bayesian inversion of marine CSEM data with a trans-dimensional self parametrizing algorithm. Geophysical Journal International191(3), 1135-1151. doi: 10.1111/j.1365-246X.2012.05677.x.
  9. Roach, I. C., McPherson, A. A., Clarke, J. D. A., Brodie, R. C., Doublier, M. P., Armistead, S. E., Skirrow, R. G. and Main, P. T. (ed), 2015. Southern Thomson Orogen VTEMplus AEM Survey: Using airborne electromagnetics as an UNCOVER aplication. Geoscience Australia. Record 2015/29. Online: https://d28rz98at9flks.cloudfront.net/83844/Rec2015_029.pdf.
/content/journals/10.1071/ASEG2018abT5_1F
Loading
  • Article Type: Research Article
Keyword(s): inversion; Magnetotelluric; Monte Carlo; trans-dimensional
This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error