Discrimination between Pressure and Fluid Saturation Changes from Time-lapse Seismic Data at a CO2 Storage Site

Author(s):  
A. Ivanova ◽  
C. Yang ◽  
J. Kummerow ◽  
S. Lueth ◽  
C. Juhlin
Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. C81-C92 ◽  
Author(s):  
Helene Hafslund Veire ◽  
Hilde Grude Borgos ◽  
Martin Landrø

Effects of pressure and fluid saturation can have the same degree of impact on seismic amplitudes and differential traveltimes in the reservoir interval; thus, they are often inseparable by analysis of a single stacked seismic data set. In such cases, time-lapse AVO analysis offers an opportunity to discriminate between the two effects. We quantify the uncertainty in estimations to utilize information about pressure- and saturation-related changes in reservoir modeling and simulation. One way of analyzing uncertainties is to formulate the problem in a Bayesian framework. Here, the solution of the problem will be represented by a probability density function (PDF), providing estimations of uncertainties as well as direct estimations of the properties. A stochastic model for estimation of pressure and saturation changes from time-lapse seismic AVO data is investigated within a Bayesian framework. Well-known rock physical relationships are used to set up a prior stochastic model. PP reflection coefficient differences are used to establish a likelihood model for linking reservoir variables and time-lapse seismic data. The methodology incorporates correlation between different variables of the model as well as spatial dependencies for each of the variables. In addition, information about possible bottlenecks causing large uncertainties in the estimations can be identified through sensitivity analysis of the system. The method has been tested on 1D synthetic data and on field time-lapse seismic AVO data from the Gullfaks Field in the North Sea.


Geophysics ◽  
2015 ◽  
Vol 80 (2) ◽  
pp. WA1-WA13 ◽  
Author(s):  
Lisa A. N. Roach ◽  
Donald J. White ◽  
Brian Roberts

Two 3D time-lapse seismic surveys were acquired in 2012 and 2013 at the Aquistore [Formula: see text] storage site prior to the start of [Formula: see text] injection. Using these surveys, we determined the background time-lapse noise at the site and assessed the feasibility of using a sparse areal permanent receiver array as a monitoring tool. Applying a standard processing sequence to these data, we adequately imaged the reservoir at 3150–3350 m depth. Evaluation of the impact of each processing step on the repeatability revealed a general monotonic increase in similarity between the data sets as a function of processing. The prestack processing sequence reduced the normalized root mean squared difference (nrms) from 1.13 between the raw stacks to 0.13 after poststack time migration. The postmigration cross-equalization sequence further reduced the global nrms to 0.07. A simulation of the changes in seismic response due to a range of [Formula: see text] injection scenarios suggested that [Formula: see text] was detectable within the reservoir at the Aquistore site provided that zones of greater thickness than 6–13 m have reached [Formula: see text] saturations of greater than 5%.


Author(s):  
J. Kasahara ◽  
A. Kato ◽  
M. Takanashi ◽  
Y. Hasada ◽  
S. Lüth ◽  
...  
Keyword(s):  

2012 ◽  
Vol 84 ◽  
pp. 14-28 ◽  
Author(s):  
Monika Ivandic ◽  
Can Yang ◽  
Stefan Lüth ◽  
Calin Cosma ◽  
Christopher Juhlin

Geophysics ◽  
2019 ◽  
Vol 85 (1) ◽  
pp. M1-M13 ◽  
Author(s):  
Yichuan Wang ◽  
Igor B. Morozov

For seismic monitoring injected fluids during enhanced oil recovery or geologic [Formula: see text] sequestration, it is useful to measure time-lapse (TL) variations of acoustic impedance (AI). AI gives direct connections to the mechanical and fluid-related properties of the reservoir or [Formula: see text] storage site; however, evaluation of its subtle TL variations is complicated by the low-frequency and scaling uncertainties of this attribute. We have developed three enhancements of TL AI analysis to resolve these issues. First, following waveform calibration (cross-equalization) of the monitor seismic data sets to the baseline one, the reflectivity difference was evaluated from the attributes measured during the calibration. Second, a robust approach to AI inversion was applied to the baseline data set, based on calibration of the records by using the well-log data and spatially variant stacking and interval velocities derived during seismic data processing. This inversion method is straightforward and does not require subjective selections of parameterization and regularization schemes. Unlike joint or statistical inverse approaches, this method does not require prior models and produces accurate fitting of the observed reflectivity. Third, the TL AI difference is obtained directly from the baseline AI and reflectivity difference but without the uncertainty-prone subtraction of AI volumes from different seismic vintages. The above approaches are applied to TL data sets from the Weyburn [Formula: see text] sequestration project in southern Saskatchewan, Canada. High-quality baseline and TL AI-difference volumes are obtained. TL variations within the reservoir zone are observed in the calibration time-shift, reflectivity-difference, and AI-difference images, which are interpreted as being related to the [Formula: see text] injection.


Geophysics ◽  
2014 ◽  
Vol 79 (6) ◽  
pp. B243-B252 ◽  
Author(s):  
Peter Bergmann ◽  
Artem Kashubin ◽  
Monika Ivandic ◽  
Stefan Lüth ◽  
Christopher Juhlin

A method for static correction of time-lapse differences in reflection arrival times of time-lapse prestack seismic data is presented. These arrival-time differences are typically caused by changes in the near-surface velocities between the acquisitions and had a detrimental impact on time-lapse seismic imaging. Trace-to-trace time shifts of the data sets from different vintages are determined by crosscorrelations. The time shifts are decomposed in a surface-consistent manner, which yields static corrections that tie the repeat data to the baseline data. Hence, this approach implies that new refraction static corrections for the repeat data sets are unnecessary. The approach is demonstrated on a 4D seismic data set from the Ketzin [Formula: see text] pilot storage site, Germany, and is compared with the result of an initial processing that was based on separate refraction static corrections. It is shown that the time-lapse difference static correction approach reduces 4D noise more effectively than separate refraction static corrections and is significantly less labor intensive.


2016 ◽  
Author(s):  
Masashi Nakatsukasa ◽  
Isao Kurosawa ◽  
Ayato Kato ◽  
Mamoru Takanashi ◽  
Don J White ◽  
...  

Geophysics ◽  
2005 ◽  
Vol 70 (3) ◽  
pp. O1-O11 ◽  
Author(s):  
Alexey Stovas ◽  
Martin Landrø

We investigate how seismic anisotropy influences our ability to distinguish between various production-related effects from time-lapse seismic data. Based on rock physics models and ultrasonic core measurements, we estimate variations in PP and PS reflectivity at the top reservoir interface for fluid saturation and pore pressure changes. The tested scenarios include isotropic shale, weak anisotropic shale, and highly anisotropic shale layers overlaying either an isotropic reservoir sand layer or a weak anisotropic sand layer. We find that, for transverse isotropic media with a vertical symmetry axis (TIV), the effect of weak anisotropy in the cap rock does not lead to significant errors in, for instance, the simultaneous determination of pore-pressure and fluid-saturation changes. On the other hand, changes in seismic anisotropy within the reservoir rock (caused by, for instance, increased fracturing) might be detectable from time-lapse seismic data. A new method using exact expressions for PP and PS reflectivity, including TIV anisotropy, is used to determine pressure and saturation changes over production time. This method is assumed to be more accurate than previous methods.


Sign in / Sign up

Export Citation Format

Share Document