skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Bayesian seismic inversion based on rock-physics 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 post-stack 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 rock-physics 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 » 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.« less

Authors:
 [1];  [2]; ;  [1];  [3];  [4]
  1. Physics Department, Federal University of Santa Catarina, Florianópolis (Brazil)
  2. Department of Geology and Geophysics, University of Wyoming, Laramie (United States)
  3. Informatic and Statistics Department, Federal University of Santa Catarina, Florianópolis (Brazil)
  4. 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, E-mail: leandrop.fgr@gmail.com, Grana, Dario, Santos, Marcio, Figueiredo, Wagner, Roisenberg, Mauro, and Schwedersky Neto, Guenther. Bayesian seismic inversion based on rock-physics 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, E-mail: leandrop.fgr@gmail.com, Grana, Dario, Santos, Marcio, Figueiredo, Wagner, Roisenberg, Mauro, & Schwedersky Neto, Guenther. Bayesian seismic inversion based on rock-physics 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, E-mail: leandrop.fgr@gmail.com, Grana, Dario, Santos, Marcio, Figueiredo, Wagner, Roisenberg, Mauro, and Schwedersky Neto, Guenther. Mon . "Bayesian seismic inversion based on rock-physics 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 rock-physics prior modeling for the joint estimation of acoustic impedance, porosity and lithofacies},
author = {Passos de Figueiredo, Leandro, E-mail: 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 post-stack 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 rock-physics 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 rock-physics relationships to link seismic attributes to electrical properties. Ideally, we can connect them through reservoir parameters (e.g., porosity and water saturation) by developing physical-based 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 rock-physics models and for subsequently obtaining reliable estimates of the associated parameters. The use of improper rock-physics models and the inaccuracy of the estimates of model parameters may cause misleading inversion results. Conversely, it is easy tomore » derive statistical relationships among seismic and electrical attributes and reservoir parameters from distant borehole logs. In this study, we develop a Bayesian model to jointly invert seismic AVA and CSEM data for reservoir parameter estimation using statistical rock-physics models; the spatial dependence of geophysical and reservoir parameters are carried out by lithotypes through Markov random fields. We apply the developed model to a synthetic case, which simulates a CO{sub 2} monitoring application. We derive statistical rock-physics relations from borehole logs at one location and estimate seismic P- and S-wave velocity ratio, acoustic impedance, density, electrical resistivity, lithotypes, porosity, and water saturation at three different locations by conditioning to seismic AVA and CSEM data. Comparison of the inversion results with their corresponding true values shows that the correlation-based statistical rock-physics models provide significant information for improving the joint inversion results.« less
  • 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 » strictly a problem of information reduction. This can be achieved in several stages. Standard deviations can be computed as overall uncertainty measures of the parameters, when the shape of the a posteriori pdf is not too far from Gaussian. Covariance and related matrices give more detailed information. An eigenvalue or principle component analysis allows the inspection of essential linear combinations of the parameters. The relative contributions of a priori information and data to the solution can be elegantly studied. Results in this paper are especially worked out for the non-linear Gaussian case. Comparisons with other approaches are given. The procedures are illustrated with a simple two-parameter inverse problem.« less
  • Astochastic joint-inversion approach for estimatingreservoir-fluid saturations and porosity is proposed. The approachcouples seismic amplitude variation with angle (AVA) and marinecontrolled-source 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 single-value estimates provided by deterministic methods, theapproach gives a probability distribution for any unknown parameter ofinterest, such as reservoir-fluid 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 » demonstrated using synthetic andfield data sets. Results show that joint inversion using seismic and EMdata gives better estimates of reservoir parameters than estimates fromeither geophysical data set used in isolation. Moreover, a moreinformative prior leads to much narrower predictive intervals of thetarget parameters, with mean values of the posterior distributions closerto logged values.« less
  • 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 layer-based, representing pseudo acoustic impedance logs at each common mid-point (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 » and removal of wavelet sidelobes. The result is a reduction in tuning effects, and improved resolution. This, in turn, results in an enhanced interpretation and more accurate quantification of reservoir properties.« less
  • Physics-based 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 1-D steady state conduit model for effusive eruptions including equilibrium crystallization and gas transport through the conduit and compare with the quasi-steady dome growth phase of Mount St. Helens in 2005. Viscosity increase resulting from pressure-dependent 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 » magma permeabilities due to their impact on magma density. Including both lateral and vertical gas transport reveals competing effects that produce nonmonotonic behavior in the mass flux when increasing magma permeability. Using this physics-based model in a Bayesian inversion, we link data sets from Mount St. Helens such as extrusion flux and earthquake depths with petrological data to estimate unknown model parameters, including magma chamber pressure and water content, magma permeability constants, conduit radius, and friction along the conduit walls. Even with this relatively simple model and limited data, we obtain improved constraints on important model parameters. We find that the magma chamber had low (<5 wt %) total volatiles and that the magma permeability scale is well constrained at ~10 –11.4m 2 to reproduce observed dome rock porosities. Here, compared with previous results, higher magma overpressure and lower wall friction are required to compensate for increased viscous resistance while keeping extrusion rate at the observed value.« less