Integrated VTI model building with seismic data, geologic information, and rock-physics modeling — Part 1: Theory and synthetic test

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. C177-C191 ◽  
Author(s):  
Yunyue Li ◽  
Biondo Biondi ◽  
Robert Clapp ◽  
Dave Nichols

Seismic anisotropy plays an important role in structural imaging and lithologic interpretation. However, anisotropic model building is a challenging underdetermined inverse problem. It is well-understood that single component pressure wave seismic data recorded on the upper surface are insufficient to resolve a unique solution for velocity and anisotropy parameters. To overcome the limitations of seismic data, we have developed an integrated model building scheme based on Bayesian inference to consider seismic data, geologic information, and rock-physics knowledge simultaneously. We have performed the prestack seismic inversion using wave-equation migration velocity analysis (WEMVA) for vertical transverse isotropic (VTI) models. This image-space method enabled automatic geologic interpretation. We have integrated the geologic information as spatial model correlations, applied on each parameter individually. We integrate the rock-physics information as lithologic model correlations, bringing additional information, so that the parameters weakly constrained by seismic are updated as well as the strongly constrained parameters. The constraints provided by the additional information help the inversion converge faster, mitigate the ambiguities among the parameters, and yield VTI models that were consistent with the underlying geologic and lithologic assumptions. We have developed the theoretical framework for the proposed integrated WEMVA for VTI models and determined the added information contained in the regularization terms, especially the rock-physics constraints.

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. C205-C218 ◽  
Author(s):  
Yunyue Li ◽  
Biondo Biondi ◽  
Robert Clapp ◽  
Dave Nichols

Velocity model building is the first step of seismic inversion and the foundation of the subsequent processing and interpretation workflow. Velocity model building from surface seismic data only becomes severely underdetermined and nonunique when more than one parameter is needed to characterize the velocity anisotropy. The traditional seismic processing workflow sequentially performs seismic velocity model building, structural imaging/interpretation, and lithologic inversion, modifying the subsurface model in each step without verifications against the previously used data. We have developed an integrated model building scheme that uses all available information: seismic data, geologic structural information, well logs, and rock-physics knowledge. We have evaluated the accuracy of the anisotropic model in the image space, in which structural information is estimated. The lithologic inversion results from well logs and the dynamic seismic information (amplitude versus angle) are also fed back to the kinematic seismic inversion via a cross-parameter covariance matrix, which is a multivariate Gaussian approximation to the numerical distribution modeled from stochastic rock-physics modeling. The procedure of building the rock-physics prior information and the improvements using these extra constraints were tested on a Gulf of Mexico data set. The inverted vertical transverse isotropic model not only better focused the seismic image, but it also satisfied the geologic and rock-physics principles.


2015 ◽  
Author(s):  
Shazia Asim* ◽  
Peimin Zhu ◽  
Tayyab Naseer ◽  
Shabeer Ahmed ◽  
Farrukh Hussain ◽  
...  

2020 ◽  
Vol 70 (1) ◽  
pp. 209-220
Author(s):  
Qazi Sohail Imran ◽  
◽  
Numair Ahmad Siddiqui ◽  
Abdul Halim Abdul Latif ◽  
Yasir Bashir ◽  
...  

Offshore petroleum systems are often very complex and subtle because of a variety of depositional environments. Characterizing a reservoir based on conventional seismic and well-log stratigraphic analysis in intricate settings often leads to uncertainties. Drilling risks, as well as associated subsurface uncertainties can be minimized by accurate reservoir delineation. Moreover, a forecast can also be made about production and performance of a reservoir. This study is aimed to design a workflow in reservoir characterization by integrating seismic inversion, petrophysics and rock physics tools. Firstly, to define litho facies, rock physics modeling was carried out through well log analysis separately for each facies. Next, the available subsurface information is incorporated in a Bayesian engine which outputs several simulations of elastic reservoir properties, as well as their probabilities that were used for post-inversion analysis. Vast areal coverage of seismic and sparse vertical well log data was integrated by geostatistical inversion to produce acoustic impedance realizations of high-resolution. Porosity models were built later using the 3D impedance model. Lastly, reservoir bodies were identified and cross plot analysis discriminated the lithology and fluid within the bodies successfully.


2019 ◽  
Vol 7 (3) ◽  
pp. SG11-SG22 ◽  
Author(s):  
Heather Bedle

Gas hydrates in the oceanic subsurface are often difficult to image with reflection seismic data, particularly when the strata run parallel to the seafloor and in regions that lack the presence of a bottom-simulating reflector (BSR). To address and understand these imaging complications, rock-physics modeling and seismic attribute analysis are performed on modern 2D lines in the Pegasus Basin in New Zealand, where the BSR is not continuously imaged. Based on rock-physics and seismic analyses, several seismic attribute methods identify weak BSR reflections, with the far-angle stack data being particularly effective. Rock modeling results demonstrate that far-offset seismic data are critical in improving the imaging and interpretation of the base of the gas hydrate stability zone. The rock-physics modeling results are applied to the Pegasus 2009 2D data set that reveals a very weak seismic reflection at the base of the hydrates in the far-angle stack. This often-discontinuous reflection is significantly weaker in amplitude than typical BSRs associated with hydrates. These weak far-angle stack BSRs often do not appear clearly in full stack data, the most commonly interpreted seismic data type. Additional amplitude variation with angle (AVA) attribute analyses provide insight into identifying the presence of gas hydrates in regions lacking a strong BSR. Although dozens of seismic attributes were investigated for their ability to reveal weak reflections at the base of the gas hydrate stability zone, those that enhance class 2 AVA anomalies were most effective, particularly the seismic fluid factor attribute.


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.


2020 ◽  
Vol 8 (2) ◽  
pp. T275-T291 ◽  
Author(s):  
Kenneth Bredesen ◽  
Esben Dalgaard ◽  
Anders Mathiesen ◽  
Rasmus Rasmussen ◽  
Niels Balling

We have seismically characterized a Triassic-Jurassic deep geothermal sandstone reservoir north of Copenhagen, onshore Denmark. A suite of regional geophysical measurements, including prestack seismic data and well logs, was integrated with geologic information to obtain facies and reservoir property predictions in a Bayesian framework. The applied workflow combined a facies-dependent calibrated rock-physics model with a simultaneous amplitude-variation-with-offset seismic inversion. The results suggest that certain sandstone distributions are potential aquifers within the target interval, which appear reasonable based on the geologic properties. However, prediction accuracy suffers from a restricted data foundation and should, therefore, only be considered as an indicator of potential aquifers. Despite these issues, the results demonstrate new possibilities for future seismic reservoir characterization and rock-physics modeling for exploration purposes, derisking, and the exploitation of geothermal energy as a green and sustainable energy resource.


Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. O53-O63 ◽  
Author(s):  
Ran Bachrach

Sediment porosity and saturation affect bulk modulus, shear modulus, and density. Consequently, estimating hydrocarbon saturation and reservoir porosity from seismic data is a joint estimation problem: Uncertainty in porosity will lead to errors in saturation prediction, and vice versa. Porosity and saturation can be jointly estimated using stochastic rock-physics modeling and formal Bayesian estimation methodology. Knowledge of shear impedance reduces the uncertainty in porosity and thus also reduces uncertainty in saturation estimation. This study investigates joint estimation of porosity and saturation by using rock-physics, stochastic modeling, and Bayesian estimation theory to derive saturation and porosity maps of expected pay sands. In the field example, the uncertainty in porosity, quantified by the standard deviation (STD) associated with the posterior probability density function (pdf), derived from inversion of seismic data is much less than the uncertainty in the derived saturation. For a typical case, the STD associated with saturation is [Formula: see text] while porosity STD is about 1.34 porosity units given seismic-derived inversion attributes with reasonable accuracy. Comparison of these numbers with prior estimates showed that inversion of seismic data decreased the uncertainty in porosity to 15% of the prior uncertainty while saturation uncertainty was only reduced to 92% of the prior uncertainty. Although these results may vary from one location to another, the methodology is general and can be applied to other locations.


Sign in / Sign up

Export Citation Format

Share Document