Quasi 3D transdimensional Markov-chain Monte Carlo for seismic impedance inversion and uncertainty analysis

Accurate estimation of subsurface properties plays an important role in successful hydrocarbon exploration, and a variety of different types of inversion schemes are used to infer earth properties such as velocity or density by analyzing the surface seismic. The Markov-chain Monte Carlo (MCMC) stochastic approach is widely used to estimate subsurface properties. We have used a transdimensional form of MCMC, reversible jump MCMC (RJMCMC), to estimate seismic impedance, which allows the inference of the number of interfaces as well as the interface location and layer impedances. Estimating the uncertainty quantitatively is also very important when performing stochastic inversion. Therefore, the goal of this paper is to apply the transdimensional method to obtain a 3D seismic impedance model and to quantify uncertainty in impedance and interface locations. We also measured the speedup of the proposed algorithm by applying data and task parallelism. To demonstrate the performance and reliability of the proposed RJMCMC impedance inversion, we used seismic data from the E-segment of the Norne field in Norwegian Sea. The results of the quasi 3D transdimensional MCMC approach, which independently inverts data from each common-depth-point location, indicate high velocity contrasts near gas-oil contacts and high uncertainty in impedance near discontinuities. Also, the cross section of the impedance uncertainty volume helps to identify the location of a high-contrast boundary corresponding to the location of the possible gas reservoir. The proposed uncertainty measure can serve as an attribute to identify important reservoir features.