Bayesian Seismic Inversion for the Joint Estimation of Facies and Elastic Properties of the Lula Field

Author(s):  
L. De Figueiredo ◽  
B. Barbosa Rodrigues ◽  
D. Grana ◽  
M. Roisenberg
2017 ◽  
Vol 336 ◽  
pp. 128-142 ◽  
Author(s):  
Leandro Passos de Figueiredo ◽  
Dario Grana ◽  
Marcio Santos ◽  
Wagner Figueiredo ◽  
Mauro Roisenberg ◽  
...  

Geophysics ◽  
2017 ◽  
Vol 82 (4) ◽  
pp. M55-M65 ◽  
Author(s):  
Xiaozheng Lang ◽  
Dario Grana

We have developed a seismic inversion method for the joint estimation of facies and elastic properties from prestack seismic data based on a geostatistical approach. The objectives of our inversion methodology are to sample from the posterior distribution of seismic properties and to simultaneously classify the lithology conditioned by seismic data. The inversion algorithm is a sequential Gaussian mixture inversion based on Bayesian linearized amplitude variation with offset inverse theory and sequential geostatistical simulations. The stochastic approach to the inversion allows generating multiple elastic models that match the seismic data. To mathematically represent the multimodal behavior of elastic properties due to their variations within different lithologies, we adopt a Gaussian mixture distribution for the prior model of the elastic properties and we use the prior probability of the facies as weights of the Gaussian components of the mixture. The solution of the inverse problem is achieved by deriving the explicit analytical expression of the posterior distribution of the elastic properties and facies and by sampling from this distribution according to a spatial correlation model. The inversion methodology has been validated using well logs and synthetic seismic data with different noise levels, and it is then applied to a real 3D seismic data set in North Sea.


2021 ◽  
Author(s):  
Zahid U. Khan ◽  
◽  
Mona Lisa ◽  
Muyyassar Hussain ◽  
Syed A. Ahmed ◽  
...  

The Pab Formation of Zamzama block, lying in the Lower Indus Basin of Pakistan, is a prominent gas-producing sand reservoir. The optimized production is limited by water encroachment in producing wells, thus it is required to distinguish the gas-sand facies from the remainder of the wet sands and shales for additional drilling zones. An approach is adopted based on a relation between petrophysical and elastic properties to characterize the prospect locations. Petro-elastic models for the identified facies are generated to discriminate lithologies in their elastic ranges. Several elastic properties, including p-impedance (11,600-12,100 m/s*g/cc), s-impedance (7,000-7,330 m/s*g/cc), and Vp/Vs ratio (1.57-1.62), are calculated from the simultaneous prestack seismic inversion, allowing the identification of gas sands in the field. Furthermore, inverted elastic attributes and well-based lithologies are incorporated into the Bayesian framework to evaluate the probability of gas sands. To better determine reservoir quality, bulk volumes of PHIE and clay are estimated using elastic volumes trained on well logs employing Probabilistic Neural Networking (PNN), which effectively handles heterogeneity effects. The results showed that the channelized gas-sands passing through existing well locations exhibited reduced clay content and maximum effective porosities of 9%, confirming the reservoir's good quality. Such approaches can be widely implemented in producing fields to completely assess litho-facies and achieve maximum production with minimal risk.


2020 ◽  
Vol 39 (2) ◽  
pp. 92-101
Author(s):  
Arturo Contreras ◽  
Andre Gerhardt ◽  
Paul Spaans ◽  
Matthew Docherty

Multiple state-of-the-art inversion methods have been implemented to integrate 3D seismic amplitude data, well logs, geologic information, and spatial variability to produce models of the subsurface. Amplitude variation with angle (AVA) deterministic, stochastic, and wave-equation-based amplitude variation with offset (WEB-AVO) inversion algorithms are used to describe Intra-Triassic Mungaroo gas reservoirs located in the Carnarvon Basin, Western Australia. The interpretation of inverted elastic properties in terms of lithology- and fluid-sensitive attributes from AVA deterministic inversion provides quantitative information about the geomorphology of fluvio-deltaic sediments as well as the delineation of gas reservoirs. AVA stochastic inversion delivers higher resolution realizations than those obtained from standard deterministic methods and allows for uncertainty analysis. Additionally, the cosimulation of petrophysical parameters from elastic properties provides precise 3D models of reservoir properties, such as volume of shale and water saturation, which can be used as part of the static model building process. Internal multiple scattering, transmission effects, and mode conversion (considered as noise in conventional linear inversion) become useful signals in WEB-AVO inversion. WEB-AVO compressibility shows increased sensitivity to residual/live gas discrimination compared to fluid-sensitive attributes obtained with conventional inversions.


Geophysics ◽  
2010 ◽  
Vol 75 (5) ◽  
pp. 75A165-75A176 ◽  
Author(s):  
Miguel Bosch ◽  
Tapan Mukerji ◽  
Ezequiel F. Gonzalez

There are various approaches for quantitative estimation of reservoir properties from seismic inversion. A general Bayesian formulation for the inverse problem can be implemented in two different work flows. In the sequential approach, first seismic data are inverted, deterministically or stochastically, into elastic properties; then rock-physics models transform those elastic properties to the reservoir property of interest. The joint or simultaneous work flow accounts for the elastic parameters and the reservoir properties, often in a Bayesian formulation, guaranteeing consistency between the elastic and reservoir properties. Rock physics plays the important role of linking elastic parameters such as impedances and velocities to reservoir properties of interest such as lithologies, porosity, and pore fluids. Geostatistical methods help add constraints of spatial correlation, conditioning to different kinds of data and incorporating subseismic scales of heterogeneities.


Geophysics ◽  
2010 ◽  
Vol 75 (3) ◽  
pp. O21-O37 ◽  
Author(s):  
Dario Grana ◽  
Ernesto Della Rossa

A joint estimation of petrophysical properties is proposed that combines statistical rock physics and Bayesian seismic inversion. Because elastic attributes are correlated with petrophysical variables (effective porosity, clay content, and water saturation) and this physical link is associated with uncertainties, the petrophysical-properties estimation from seismic data can be seen as a Bayesian inversion problem. The purpose of this work was to develop a strategy for estimating the probability distributions of petrophysical parameters and litho-fluid classes from seismics. Estimation of reservoir properties and the associated uncertainty was performed in three steps: linearized seismic inversion to estimate the probabilities of elastic parameters, probabilistic upscaling to include the scale-changes effect, and petrophysical inversion to estimate the probabilities of petrophysical variables andlitho-fluid classes. Rock-physics equations provide the linkbetween reservoir properties and velocities, and linearized seismic modeling connects velocities and density to seismic amplitude. A full Bayesian approach was adopted to propagate uncertainty from seismics to petrophysics in an integrated framework that takes into account different sources of uncertainty: heterogeneity of the real data, approximation of physical models, measurement errors, and scale changes. The method has been tested, as a feasibility step, on real well data and synthetic seismic data to show reliable propagation of the uncertainty through the three different steps and to compare two statistical approaches: parametric and nonparametric. Application to a real reservoir study (including data from two wells and partially stacked seismic volumes) has provided as a main result the probability densities of petrophysical properties and litho-fluid classes. It demonstrated the applicability of the proposed inversion method.


2021 ◽  
Author(s):  
Pavlo Kuzmenko ◽  
Rustem Valiakhmetov ◽  
Francesco Gerecitano ◽  
Viktor Maliar ◽  
Grigori Kashuba ◽  
...  

Abstract The seismic data have historically been utilized to perform structural interpretation of the geological subsurface. Modern approaches of Quantitative Interpretation are intended to extract geologically valuable information from the seismic data. This work demonstrates how rock physics enables optimal prediction of reservoir properties from seismic derived attributes. Using a seismic-driven approach with incorporated prior geological knowledge into a probabilistic subsurface model allowed capturing uncertainty and quantifying the risk for targeting new wells in the unexplored areas. Elastic properties estimated from the acquired seismic data are influenced by the depositional environment, fluid content, and local geological trends. By applying the rock physics model, we were able to predict the elastic properties of a potential lithology away from the well control points in the subsurface whether or not it has been penetrated. Seismic amplitude variation with incident angle (AVO) and azimuth (AVAZ) jointly with rock-derived petrophysical interpretations were used for stochastical modeling to capture the reservoir distribution over the deep Visean formation. The seismic inversion was calibrated by available well log data and by traditional structural interpretation. Seismic elastic inversion results in a deep Lower Carboniferous target in the central part of the DDB are described. The fluid has minimal effect on the density and Vp. Well logs with cross-dipole acoustics are used together with wide-azimuth seismic data, processed with amplitude control. It is determined that seismic anisotropy increases in carbonate deposits. The result covers a set of lithoclasses and related probabilities: clay minerals, tight sandstones, porous sandstones, and carbonates. We analyzed the influence of maximum angles determination for elastic inversion that varied from 32.5 to 38.5 degrees. The greatest influence of the far angles selection is on the density. AI does not change significantly. Probably the 38,5 degrees provides a superior response above the carbonates. It does not seem to damage the overall AVA behavior, which result in a good density outcome, as higher angles of incidence are included. It gives a better tie to the wells for the high density layers over the interval of interest. Sand probability cube must always considered in the interpretation of the lithological classification that in many cases may be misleading (i.e. when sand and shale probabilities are very close to each other, because of small changes in elastic parameters). The authors provide an integrated holistic approach for quantitative interpretation, subsurface modeling, uncertainty evaluation, and characterization of reservoir distribution using pre-existing well logs and recently acquired seismic data. This paper underpins the previous efforts and encourages the work yet to be fulfilled on this subject. We will describe how quantitative interpretation was used for describing the reservoir, highlight values and uncertainties, and point a way forward for further improvement of the process for effective subsurface modeling.


Sign in / Sign up

Export Citation Format

Share Document