Markov chain Monte Carlo for seismic facies classification

Dario Grana,Leandro de Figueiredo,Klaus Mosegaard
DOI: https://doi.org/10.1190/geo2022-0442.1
IF: 3.264
2023-04-05
Geophysics
Abstract:Seismic facies classification aims to predict a facies model, or a set of facies models, from measured seismic data. We focus on stochastic classification methods to estimate the probability distribution of facies conditioned on seismic data. Bayesian classification methods based on analytical solutions are generally applied to seismically inverted elastic attributes or petrophysical properties rather than measured seismic data and they do not include spatial correlation models in the prior distribution. Instead, iterative stochastic methods such as Monte Carlo acceptance-rejection sampling and Approximate Bayesian Computation can be directly applied to seismic data and can be implemented with complex prior models; however, their convergence is generally slow as it requires the simulation of a large number of facies models from the prior distribution. We propose an efficient implementation of a Markov chain Monte Carlo (MCMC) approach that adopt complex prior models, such as multiple point statistics (MPS) simulations based on a training image (TI), to generate geologically realistic facies realizations. The novelty of the approach is that the proposal distribution of the proposed MCMC method is based on a gradual deformation algorithm, where the proposal distribution is expressed as a linear combination of the indicator variables of the previously accepted model and a random prior simulation, according to a uniform random weight. The synthetic examples show accurate facies predictions with a success rate of approximately 80% for well log facies classification based on elastic attributes and approximately 60% for seismic facies classification based on seismic data. The proposed inversion is compared to a MCMC method with prior models sampled from a first order Markov chain (MC) and Bayesian facies classification.
geochemistry & geophysics
What problem does this paper attempt to address?