The conventional Markov chain Monte Carlo (MCMC) method is limited to the selected shape and size of proposal distribution and is not easy to start when the initial proposal distribution is far away from the target distribution. To overcome these drawbacks of the conventional MCMC method, two useful improvements in MCMC method, adaptive Metropolis (AM) algorithm and delayed rejection (DR) algorithm, are attempted to be combined. The AM algorithm aims at adapting the proposal distribution by using the generated estimators, and the DR algorithm aims at enhancing the efficiency of the improved MCMC method. Based on the improved MCMC method, a Bayesian amplitude versus offset (AVO) inversion method on the basis of the exact Zoeppritz equation has been developed, with which the P- and S-wave velocities and the density can be obtained directly, and the uncertainty of AVO inversion results has been estimated as well. The study based on the logging data and the seismic data demonstrates the feasibility and robustness of the method and shows that all three parameters are well retrieved. So the exact Zoeppritz-based nonlinear inversion method by using the improved MCMC is not only suitable for reservoirs with strong-contrast interfaces and long-offset ranges but also it is more stable, accurate and anti-noise.
Keywords: Adaptive Metropolis (AM) algorithm; Bayesian AVO inversion; Delayed rejection (DR) algorithm; Exact Zoeppritz; Nonlinear inversion.