Home > Research > Publications & Outputs > Efficient multiscale imaging of subsurface resi...

Associated organisational unit

Electronic data

  • Tso et al 2021 EKI (author accepted copy)

    Rights statement: This is a pre-copy-editing, author-produced PDF of an article accepted for publication in Geophysical Journal International following peer review. The definitive publisher-authenticated version Chak-Hau Michael Tso, Marco Iglesias, Paul Wilkinson, Oliver Kuras, Jonathan Chambers, Andrew Binley, Efficient multiscale imaging of subsurface resistivity with uncertainty quantification using ensemble Kalman inversion, Geophysical Journal International, Volume 225, Issue 2, May 2021, Pages 887–905, https://doi.org/10.1093/gji/ggab013 is available online at: https://academic.oup.com/gji/article-abstract/225/2/887/6081097

    Accepted author manuscript, 2.37 MB, PDF document

    Available under license: CC BY-NC: Creative Commons Attribution-NonCommercial 4.0 International License

Links

Text available via DOI:

View graph of relations

Efficient multiscale imaging of subsurface resistivity with uncertainty quantification using ensemble Kalman inversion

Research output: Contribution to Journal/MagazineJournal articlepeer-review

Published
  • C.-H.M. Tso
  • M. Iglesias
  • P. Wilkinson
  • O. Kuras
  • J. Chambers
  • A. Binley
Close
<mark>Journal publication date</mark>31/05/2021
<mark>Journal</mark>Geophysical Journal International
Issue number2
Volume225
Number of pages19
Pages (from-to)887-905
Publication StatusPublished
Early online date11/01/21
<mark>Original language</mark>English

Abstract

Electrical resistivity tomography (ERT) is widely used to image the Earth's subsurface and has proven to be an extremely useful tool in application to hydrological problems. Conventional smoothness-constrained inversion of ERT data is efficient and robust, and consequently very popular. However, it does not resolve well sharp interfaces of a resistivity field and tends to reduce and smooth resistivity variations. These issues can be problematic in a range of hydrological or near-surface studies, for example mapping regolith-bedrock interfaces. While fully Bayesian approaches, such as those using Markov chain Monte Carlo sampling, can address the above issues, their very high computation cost makes them impractical for many applications. Ensemble Kalman inversion (EKI) offers a computationally efficient alternative by approximating the Bayesian posterior distribution in a derivative-free manner, which means only a relatively small number of 'black-box' model runs are required. Although common limitations for ensemble Kalman filter-type methods apply to EKI, it is both efficient and generally captures uncertainty patterns correctly. We propose the use of a new EKI-based framework for ERT which estimates a resistivity model and its uncertainty at a modest computational cost. Our EKI framework uses a level-set parametrization of the unknown resistivity to allow efficient estimation of discontinuous resistivity fields. Instead of estimating level-set parameters directly, we introduce a second step to characterize the spatial variability of the resistivity field and infer length scale hyperparameters directly. We demonstrate these features by applying the method to a series of synthetic and field examples. We also benchmark our results by comparing them to those obtained from standard smoothness-constrained inversion. Resultant resistivity images from EKI successfully capture arbitrarily shaped interfaces between resistivity zones and the inverted resistivities are close to the true values in synthetic cases. We highlight its readiness and applicability to similar problems in geophysics.

Bibliographic note

This is a pre-copy-editing, author-produced PDF of an article accepted for publication in Geophysical Journal International following peer review. The definitive publisher-authenticated version Chak-Hau Michael Tso, Marco Iglesias, Paul Wilkinson, Oliver Kuras, Jonathan Chambers, Andrew Binley, Efficient multiscale imaging of subsurface resistivity with uncertainty quantification using ensemble Kalman inversion, Geophysical Journal International, Volume 225, Issue 2, May 2021, Pages 887–905, https://doi.org/10.1093/gji/ggab013 is available online at: https://academic.oup.com/gji/article-abstract/225/2/887/6081097