Single-step probabilistic inversion of 3D seismic data of a carbonate reservoir in Southwest Iran

Abstract

We use a sampling-based Markov chain Monte Carlo method to invert seismic data directly for porosity and quantify its uncertainty distribution in a hard-rock carbonate reservoir in Southwest Iran. The noise that remains on seismic data after the processing flow is correlated with the bandwidth in the range of the seismic wavelet. Hence, to account for the inherent correlated nature of the band-limited seismic noise in the probabilistic inversion of real seismic data, we assume the estimated seismic wavelet as a suitable proxy for capturing the coupling of noise samples. In contrast to the common approach of inserting a delta function on the main diagonal of the covariance matrix, we insert the seismic wavelet on its main diagonal. We also calibrate an empirical and a semi-empirical inclusion-based rock-physics model to characterize the rock-frame elastic moduli via a lithology constrained fitting of the parameters of these models, i.e. the critical porosity and the pore aspect ratio. These calibrated rock-physics models are embedded in the inversion procedure to link petrophysical and elastic properties. In addition, we obtain the pointwise critical porosity and pore aspect ratio, which can potentially facilitate the interpretation of the reservoir for further studies. The inversion results are evaluated by comparing with porosity logs and an existing geological model (porosity model) constructed through a geostatistical simulation approach. We assess the consistency of the geological model through a geomodel-to-seismic modeling approach. The results confirm the performance of the probabilistic inversion in resolving some thin layers and reconstructing the observed seismic data. We present the applicability of the proposed sampling-based approach to invert 3D seismic data for estimating the porosity distribution and its associated uncertainty for four subzones of the reservoir. The porosity time maps and the facies probabilities obtained via porosity cut-offs indicate the relative quality of the reservoir’s subzones.

Assessment of the impact of noise magnitude and bandwidth variations on a probabilistic inversion of reflection seismic data

“…As showcased by Heidari et al. (2022), the results of this work can directly be applied in industry‐led studies, where the noise model is a challenging aspect in setting up the probabilistic inversion of real seismic data.…”

“…Subsequently, we present a thorough discussion on the impact of different noise models on the results of the probabilistic inversion of 2D realistic synthetic seismic data. As showcased by Heidari et al (2022), the results of this work can directly be applied in industry-led studies, where the noise model is a challenging aspect in setting up the probabilistic inversion of real seismic data.…”