Seismic reservoir characterization uses pre-stack reflection seismic data to describe the spatial variability of subsurface properties around the target zone. In this work, I implement a non linear, target-oriented, inversion algorithm that exploits the amplitude versus angle (AVA) variations of the seismic reflections to estimate in a single-step both the litho-fluid facies and the petrophysical properties. The inversion is casted in a Bayesian framework and given the non linearity of the forward modeling, it is formulated in terms of a Markov Chain Monte Carlo (MCMC) algorithm (Sambridge and Moosegard, 2002) in order to produce accurate and unbiased uncertainty estimations. I apply this algorithm in a clastic reservoir located in offshore Nile Delta where the reservoir zone is gas saturated and hosted in sand channels surrounded by shale sequences. A linear empirical rock physics model (RPM; see Aleardi et al. 2016) is used to link the petrophysical properties to the elastic parameters, whereas the non-linear Zoeppritz equations relate such elastic properties to the observed AVA response. The exact Zoeppritz equations allow me to take advantage of the long offset seismic acquisition and to consider a wide range of incidence angles (0 and 60 degrees) in the inversion. The Gaussian mixture (GM) distribution used to describe the a-priori information about the petrophysical properties takes into consideration the multimodality and the correlation that characterize the distribution of these properties in the reservoir zone. In the field data application the close match between the outcomes of the MCMC algorithm and the well log information demonstrates the applicability of the method and the reliability of the final results. In the following discussion, F indicates the litho-fluid facies that are shale, brine sand and gas sand, E represents the elastic properties that are P- wave, S-wave velocities (Vp and Vs, respectively) and density, R indicates the petrophysical properties that are water saturation (Sw), porosity (φ) and shaliness (Sh), whereas d is the observed data that is the AVA response pertaining to the top of the interpreted reservoir extracted for each considered CMP gather.
A Markov Chain Monte Carlo algorithm for litho-fluid facies prediction and petrophysical property estimation: An application for reservoir characterization in offshore Nile Delta
ALEARDI, MATTIA
2016-01-01
Abstract
Seismic reservoir characterization uses pre-stack reflection seismic data to describe the spatial variability of subsurface properties around the target zone. In this work, I implement a non linear, target-oriented, inversion algorithm that exploits the amplitude versus angle (AVA) variations of the seismic reflections to estimate in a single-step both the litho-fluid facies and the petrophysical properties. The inversion is casted in a Bayesian framework and given the non linearity of the forward modeling, it is formulated in terms of a Markov Chain Monte Carlo (MCMC) algorithm (Sambridge and Moosegard, 2002) in order to produce accurate and unbiased uncertainty estimations. I apply this algorithm in a clastic reservoir located in offshore Nile Delta where the reservoir zone is gas saturated and hosted in sand channels surrounded by shale sequences. A linear empirical rock physics model (RPM; see Aleardi et al. 2016) is used to link the petrophysical properties to the elastic parameters, whereas the non-linear Zoeppritz equations relate such elastic properties to the observed AVA response. The exact Zoeppritz equations allow me to take advantage of the long offset seismic acquisition and to consider a wide range of incidence angles (0 and 60 degrees) in the inversion. The Gaussian mixture (GM) distribution used to describe the a-priori information about the petrophysical properties takes into consideration the multimodality and the correlation that characterize the distribution of these properties in the reservoir zone. In the field data application the close match between the outcomes of the MCMC algorithm and the well log information demonstrates the applicability of the method and the reliability of the final results. In the following discussion, F indicates the litho-fluid facies that are shale, brine sand and gas sand, E represents the elastic properties that are P- wave, S-wave velocities (Vp and Vs, respectively) and density, R indicates the petrophysical properties that are water saturation (Sw), porosity (φ) and shaliness (Sh), whereas d is the observed data that is the AVA response pertaining to the top of the interpreted reservoir extracted for each considered CMP gather.File | Dimensione | Formato | |
---|---|---|---|
A Markov chain Monte Carlo algorithm.pdf
accesso aperto
Tipologia:
Versione finale editoriale
Licenza:
Tutti i diritti riservati (All rights reserved)
Dimensione
3.11 MB
Formato
Adobe PDF
|
3.11 MB | Adobe PDF | Visualizza/Apri |
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.