Bayesian seismic inversion based on rockphysics prior modeling for the joint estimation of acoustic impedance, porosity and lithofacies
Abstract
We propose a Bayesian approach for seismic inversion to estimate acoustic impedance, porosity and lithofacies within the reservoir conditioned to poststack seismic and well data. The link between elastic and petrophysical properties is given by a joint prior distribution for the logarithm of impedance and porosity, based on a rockphysics model. The well conditioning is performed through a background model obtained by well log interpolation. Two different approaches are presented: in the first approach, the prior is defined by a single Gaussian distribution, whereas in the second approach it is defined by a Gaussian mixture to represent the well data multimodal distribution and link the Gaussian components to different geological lithofacies. The forward model is based on a linearized convolutional model. For the single Gaussian case, we obtain an analytical expression for the posterior distribution, resulting in a fast algorithm to compute the solution of the inverse problem, i.e. the posterior distribution of acoustic impedance and porosity as well as the facies probability given the observed data. For the Gaussian mixture prior, it is not possible to obtain the distributions analytically, hence we propose a Gibbs algorithm to perform the posterior sampling and obtain several reservoir model realizations, allowing anmore »
 Authors:
 Physics Department, Federal University of Santa Catarina, Florianópolis (Brazil)
 Department of Geology and Geophysics, University of Wyoming, Laramie (United States)
 Informatic and Statistics Department, Federal University of Santa Catarina, Florianópolis (Brazil)
 Petrobras Research Center, Rio de Janeiro (Brazil)
 Publication Date:
 OSTI Identifier:
 22622286
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 336; Other Information: Copyright (c) 2017 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 97 MATHEMATICAL METHODS AND COMPUTING; ACOUSTICS; ALGORITHMS; COMPARATIVE EVALUATIONS; COMPUTERIZED SIMULATION; DATA ANALYSIS; DISTRIBUTION; GAUSS FUNCTION; IMPEDANCE; INTERPOLATION; MONTE CARLO METHOD; POROSITY; PROBABILITY; ROCKS; STACKS
Citation Formats
Passos de Figueiredo, Leandro, Email: leandrop.fgr@gmail.com, Grana, Dario, Santos, Marcio, Figueiredo, Wagner, Roisenberg, Mauro, and Schwedersky Neto, Guenther. Bayesian seismic inversion based on rockphysics prior modeling for the joint estimation of acoustic impedance, porosity and lithofacies. United States: N. p., 2017.
Web. doi:10.1016/J.JCP.2017.02.013.
Passos de Figueiredo, Leandro, Email: leandrop.fgr@gmail.com, Grana, Dario, Santos, Marcio, Figueiredo, Wagner, Roisenberg, Mauro, & Schwedersky Neto, Guenther. Bayesian seismic inversion based on rockphysics prior modeling for the joint estimation of acoustic impedance, porosity and lithofacies. United States. doi:10.1016/J.JCP.2017.02.013.
Passos de Figueiredo, Leandro, Email: leandrop.fgr@gmail.com, Grana, Dario, Santos, Marcio, Figueiredo, Wagner, Roisenberg, Mauro, and Schwedersky Neto, Guenther. Mon .
"Bayesian seismic inversion based on rockphysics prior modeling for the joint estimation of acoustic impedance, porosity and lithofacies". United States.
doi:10.1016/J.JCP.2017.02.013.
@article{osti_22622286,
title = {Bayesian seismic inversion based on rockphysics prior modeling for the joint estimation of acoustic impedance, porosity and lithofacies},
author = {Passos de Figueiredo, Leandro, Email: leandrop.fgr@gmail.com and Grana, Dario and Santos, Marcio and Figueiredo, Wagner and Roisenberg, Mauro and Schwedersky Neto, Guenther},
abstractNote = {We propose a Bayesian approach for seismic inversion to estimate acoustic impedance, porosity and lithofacies within the reservoir conditioned to poststack seismic and well data. The link between elastic and petrophysical properties is given by a joint prior distribution for the logarithm of impedance and porosity, based on a rockphysics model. The well conditioning is performed through a background model obtained by well log interpolation. Two different approaches are presented: in the first approach, the prior is defined by a single Gaussian distribution, whereas in the second approach it is defined by a Gaussian mixture to represent the well data multimodal distribution and link the Gaussian components to different geological lithofacies. The forward model is based on a linearized convolutional model. For the single Gaussian case, we obtain an analytical expression for the posterior distribution, resulting in a fast algorithm to compute the solution of the inverse problem, i.e. the posterior distribution of acoustic impedance and porosity as well as the facies probability given the observed data. For the Gaussian mixture prior, it is not possible to obtain the distributions analytically, hence we propose a Gibbs algorithm to perform the posterior sampling and obtain several reservoir model realizations, allowing an uncertainty analysis of the estimated properties and lithofacies. Both methodologies are applied to a real seismic dataset with three wells to obtain 3D models of acoustic impedance, porosity and lithofacies. The methodologies are validated through a blind well test and compared to a standard Bayesian inversion approach. Using the probability of the reservoir lithofacies, we also compute a 3D isosurface probability model of the main oil reservoir in the studied field.},
doi = {10.1016/J.JCP.2017.02.013},
journal = {Journal of Computational Physics},
number = ,
volume = 336,
place = {United States},
year = {Mon May 01 00:00:00 EDT 2017},
month = {Mon May 01 00:00:00 EDT 2017}
}

Joint inversion of seismic AVA and CSEM data requires rockphysics relationships to link seismic attributes to electrical properties. Ideally, we can connect them through reservoir parameters (e.g., porosity and water saturation) by developing physicalbased models, such as Gassmann’s equations and Archie’s law, using nearby borehole logs. This could be difficult in the exploration stage because information available is typically insufficient for choosing suitable rockphysics models and for subsequently obtaining reliable estimates of the associated parameters. The use of improper rockphysics models and the inaccuracy of the estimates of model parameters may cause misleading inversion results. Conversely, it is easy tomore »

Bayesian estimation in seismic inversion
A parameter estimation or inversion procedure is incomplete without an analysis of uncertainties in the results. In the fundamental approach of Bayesian parameter estimation, the a posteriori probability density function (pdf) is the solution to the inverse problem. It is the product of the a priori pdf, containing a priori information on the parameters, and the likelihood function, which represents the information from the data. The maximum of the a posteriori pdf is usually taken as a point estimate of the parameters. The shape of this pdf, however, gives the full picture of uncertainty in the parameters. Uncertainty analysis ismore » 
Reservoirparameter identification using minimum relativeentropybased Bayesian inversion of seismic AVA and marine CSEMdata
Astochastic jointinversion approach for estimatingreservoirfluid saturations and porosity is proposed. The approachcouples seismic amplitude variation with angle (AVA) and marinecontrolledsource electromagnetic (CSEM) forward models into a Bayesianframework, which allows for integration of complementary information. Toobtain minimally subjective prior probabilities required for the Bayesianapproach, the principle of minimum relative entropy (MRE) is employed.Instead of singlevalue estimates provided by deterministic methods, theapproach gives a probability distribution for any unknown parameter ofinterest, such as reservoirfluid saturations or porosity at variouslocations. The distribution means, modes, and confidence intervals can becalculated, providing a more complete understanding of the uncertainty inthe parameter estimates. The approach ismore » 
Seismic inversion and modeling for porosity detection
The inversion of seismic reflection data to acoustic impedance rock property data presents opportunities for improved interpretation leading to a significant reduction of risk in exploration and development drilling. While the traditional seismic section represents boundaries between impedance contrasts, seismic inversion output is layerbased, representing pseudo acoustic impedance logs at each common midpoint (CMP). The earth is correctly represented as a complex of geologic layers, and in this form, is immediately recognizable to both geophysicists and geologists. Modern inversion methods also incorporate wave shaping techniques in the inversion process. This facilitates both a correction to zero phase and the estimationmore » 
Constraining the Magmatic System at Mount St. Helens (20042008) Using Bayesian Inversion With PhysicsBased Models Including Gas Escape and Crystallization
Physicsbased models of volcanic eruptions track conduit processes as functions of depth and time. When used in inversions, these models permit integration of diverse geological and geophysical data sets to constrain important parameters of magmatic systems. We develop a 1D steady state conduit model for effusive eruptions including equilibrium crystallization and gas transport through the conduit and compare with the quasisteady dome growth phase of Mount St. Helens in 2005. Viscosity increase resulting from pressuredependent crystallization leads to a natural transition from viscous flow to frictional sliding on the conduit margin. Erupted mass flux depends strongly on wall rock andmore »