1887

Abstract

Summary

We implement a transdimensional Bayesian approach to solve the 1D elastic full-waveform inversion (FWI) in which the reflectivity algorithm constitutes the forward modelling. In this approach the number of model parameters (i.e. the number of layers) is treated as an unknown, and a reversible jump Markov Chain Monte Carlo algorithm is used to sample the variable-dimension model space. We also treat the noise standard deviation as an unknown parameter to be solved for, thus letting the algorithm infer the appropriate level of data-fitting. A Parallel tempering strategy and a delayed rejection updating scheme are used to improve the efficiency of the probabilistic sampling. We focus the attention to synthetic data inversions, with the aim to draw general conclusions about the suitability of our approach for pre-stack inversion of reflection seismic data. Our tests prove that the implemented inversion algorithm provides a parsimonious solution and successfully estimates model uncertainty, noise level, model dimensionality and elastic parameters. Our experiments also demonstrate that there is a trade-off between property uncertainty and location uncertainty: A strong elastic contrast determines high uncertainty in the model property values, but low uncertainty in the location of the elastic discontinuity.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.201901339
2019-06-03
2020-04-08
Loading full text...

Full text loading...

References

  1. Bodin, T., and Sambridge, M.
    (2009). Seismic tomography with the reversible jump algorithm. Geophysical Journal International, 178(3), 1411–1436.
    [Google Scholar]
  2. Bodin, T., Sambridge, M., Tkalčić, H., Arroucau, P., Gallagher, K., and Rawlinson, N.
    (2012). Transdimensional inversion of receiver functions and surface wave dispersion. Journal of Geophysical Research: Solid Earth, 117, B02301.
    [Google Scholar]
  3. Mandolesi, E., Ogaya, X., Campanyà, J., and Agostinetti, N. P.
    (2018). A reversible-jump Markov chain Monte Carlo algorithm for 1D inversion of magnetotelluric data. Computers & Geosciences, 113, 94–105.
    [Google Scholar]
  4. Sambridge, M.
    (2014). A parallel tempering algorithm for probabilistic sampling and multimodal optimization. Geophysical Journal International, 196(1), 357–374.
    [Google Scholar]
  5. Sen, M. K., and Biswas, R.
    (2017). Transdimensional seismic inversion using the reversible jump Hamiltonian Monte Carlo algorithm. Geophysics, 82(3), R119–R134.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.201901339
Loading
/content/papers/10.3997/2214-4609.201901339
Loading

Data & Media loading...

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