Estimating saturation and density changes caused by CO2 injection at Sleipner — Using time-lapse seismic amplitude-variation-with-offset and time-lapse gravity

2017 ◽  
Vol 5 (2) ◽  
pp. T243-T257 ◽  
Author(s):  
Martin Landrø ◽  
Mark Zumberge

We have developed a calibrated, simple time-lapse seismic method for estimating saturation changes from the [Formula: see text]-storage project at Sleipner offshore Norway. This seismic method works well to map changes when [Formula: see text] is migrating laterally away from the injection point. However, it is challenging to detect changes occurring below [Formula: see text] layers that have already been charged by some [Formula: see text]. Not only is this partly caused by the seismic shadow effects, but also by the fact that the velocity sensitivity for [Formula: see text] change in saturation from 0.3 to 1.0 is significantly less than saturation changes from zero to 0.3. To circumvent the seismic shadow zone problem, we combine the time-lapse seismic method with time-lapse gravity measurements. This is done by a simple forward modeling of gravity changes based on the seismically derived saturation changes, letting these saturation changes be scaled by an arbitrary constant and then by minimizing the least-squares error to obtain the best fit between the scaled saturation changes and the measured time-lapse gravity data. In this way, we are able to exploit the complementary properties of time-lapse seismic and gravity data.

2020 ◽  
Author(s):  
F. Greco ◽  
D. Carbone ◽  
F. Cannavò ◽  
A. A. Messina ◽  
G. Siligato

AbstractGravity changes associated with volcanic processes occur over a wide range of time scales, from minutes to years and with magnitudes between a few and a few hundred microGal. High-precision instruments are needed to detect such small signals and both time-lapse surveys along networks of stations, and continuous measurements at single points, are accomplished. Continuous volcano gravimetry is mostly carried out through relative gravimeters, either superconducting instruments, providing higher quality data, or the more widely used spring meters. On the other hand, time-lapse surveys can be carried out with relative (spring) gravimeters, that measure gravity differences between pairs of stations, or by absolute gravimeters, capable of measuring the absolute value of the gravitational acceleration at the observation point. Here we present the state-of-the-art of terrestrial gravity measurements to monitor and study active volcanoes and the possibilities of new gravimeters that are under development. In particular, we present data from a mini array of three iGrav superconducting gravimeters (SGs) at Mount Etna (the first network of SGs ever installed on an active volcano). A comparison between continuous gravity measurements recorded through the iGrav#016 superconducting gravimeter at Serra La Nave station (1730 m a.s.l.) and absolute gravity data collected with the Microg LaCoste FG5#238 gravimeter in the framework of repeated campaigns is also presented. Furthermore, we introduce the Horizon 2020 NEWTON-g project (New Tools for Terrain Gravimetry), funded under the FET-OPEN Research and Innovation Actions call, Work Programme 2016–2017 (Grant Agreement No 801221). In the framework of this project, we aim to develop a field-compatible gravity imager, including an array of low-costs Micro-Electro-Mechanical Systems (MEMS)-based relative gravimeters, anchored on an absolute quantum gravimeter. After the design and production phases, the gravity imager will be field-tested at Mt. Etna (Italy) during the last 2 years of the project.


Geophysics ◽  
2007 ◽  
Vol 72 (1) ◽  
pp. C9-C17 ◽  
Author(s):  
Aaron Wandler ◽  
Brian Evans ◽  
Curtis Link

Information on time-lapse changes in seismic amplitude variation with offset (AVO) from a reservoir can be used to optimize production. We designed a scaled physical model experiment to study the AVO response of mixtures of brine, oil, and carbon dioxide at pressures of 0, 1.03, and [Formula: see text]. The small changes in density and velocity for each fluid because of increasing pressure were not detectable and were assumed to lie within the error of the experiment. However, AVO analysis was able to detect changes in the elastic properties between fluids that contained oil and those that did not. When the AVO response was plotted in the AVO intercept-gradient domain, fluids containing oil were clearly separated from fluids not containing oil. This was observed in the AVO response from both the top and base of the fluids in the physical model. We then compared the measured AVO response with the theoretical AVO response given by the Zoeppritz equations. The measured and theoretical AVO intercept responses for the top fluid reflection agree well, although the AVO gradients disagree slightly. For the fluid base reflection, the measured and theoretical responses are in close agreement.


Geophysics ◽  
2008 ◽  
Vol 73 (6) ◽  
pp. WA155-WA161 ◽  
Author(s):  
Håvard Alnes ◽  
Ola Eiken ◽  
Torkjell Stenvold

Thirty seafloor gravity stations have been placed above the carbon dioxide [Formula: see text] injection site and producing gas reservoir at the Sleipner Øst Ty field. Gravity and depth measurements from 2002 and 2005 reveal vertical changes of the permanently deployed benchmarks, probably caused by seafloor erosion and biologic activity (fish). The original gravity data have been reprocessed, resulting in slightly different gravity-change values compared with earlier published results. Observed gravity changes are caused by height variances, gas production and water influx in the Ty Formation, and [Formula: see text] injection in the Utsira Formation. Simultaneous matches to models for these effects have been made. The latest simulation model of the Ty Formation was fitted by permitting a scale factor, and the gravity contribution from the [Formula: see text] plume was determined by using the plume geometry as observed in 4D seismic data and varying the average density. The best-fit vertical gravity gradient is [Formula: see text], and the response from the Ty Formation suggests more water influx than expected in the presurvey simulation model. The best-fit average density of [Formula: see text] is [Formula: see text]. Estimates of the reservoir temperature combined with the equation of state for [Formula: see text] indicate an upper bound on [Formula: see text] density of [Formula: see text]. The gravity data suggest a lower bound of [Formula: see text] at 95% confidence.


2021 ◽  
Author(s):  
Pankaj Kumar Tiwari ◽  
Debasis Priyadarshan Das ◽  
Parimal Arjun Patil ◽  
Prasanna Chidambaram ◽  
Prasanna Kumar Chandran ◽  
...  

Abstract CO2 sequestration in depleted carbonate reservoir stipulate incorporation of comprehensive and trailblazing monitoring technologies. 4D time-lapse seismic is sine qua non for Monitoring, Measurement and Verification (MMV) planning to demonstrate the migration of CO2 plume within geological storage. An ingenious, adaptive and site specific MMV plan for monitoring CO2 plume is paramount to minimize possible subsurface and project integrity risks. Integration of dynamic simulation with seismic forward modeling aggrandize the capabilities of 4D seismic in CO2 sequestration projects. Depleted carbonate reservoir has been thoroughly studied and its geomechanical and geochemical modeling results were coupled into dynamic simulation. Reservoir porosity and fluid properties along with CO2 saturation and injection pressure distribution within each reservoir level were generated. The dynamic simulation results were integrated with seismic forward modeling to demonstrate the CO2 plume migration and its impact on seismic amplitude. Fluid acoustic properties were computed for carbonate reservoir using FLAG method. Selection of wells was based on availability of superior quality acoustic logs as well as those representing the reservoir best. Gassmann fluid substitution exercise was carried using dry rock modeling. Several scenarios were generated, and results were analyzed to demonstrate the effect of CO2 saturation and pressure build-ups within reservoir on the seismic amplitude due to continuous CO2 injection. Synthetic seismic AVO gathers were generated for angles ranging from 5 to 50 degree. Near, Mid and Far seismic amplitude response at the top of carbonate reservoir were analyzed with respect to in-situ condition for each scenario. Results reveal that CO2 saturation as low as 25 - 30% in depleted carbonate reservoir can be distinguished from 4D time-lapse seismic. With continuous CO2 injection, the reservoir pressure increases and this in turn controls the properties of both in-situ and injected fluids. The gradual changes in fluid properties and their impact on bulk acoustic properties of reservoir were modeled to assess the feasibility of using 4D seismic as a predictive tool for detection of localized and provincial pressure build-ups. Modeling results show that although observed changes in amplitude on synthetic gathers were subtle, it is expected that 4D seismic with high signal-to-noise ratio possibly be able to image such localized pressure build-ups. To monitor CO2 plume migration as well as localized pressure build-ups, we recommend acquiring multi-azimuth (MAZ) surface seismic in combination with 3D DAS-VSP for superior subsurface imaging. The integrated modeling approach ensures that 4D Seismic in subsurface CO2 plume monitoring is robust. Monitoring pressure build-ups from MAZ surface seismic and 3D DAS-VSP will reduce the associated risks.


Geophysics ◽  
2004 ◽  
Vol 69 (1) ◽  
pp. 25-36 ◽  
Author(s):  
Shelley J. Ellison ◽  
Matthias G. Imhof ◽  
Cahit Çoruh ◽  
Alan D. Fuqua ◽  
Stephen C. Henry

The objective of this case study is to predict whether 5 years of water‐flood production from a thinly layered Gulf of Mexico reservoir will change its seismic amplitude‐variation‐with‐offset (AVO) response in a detectable manner. Density and velocity profiles were computed from in situ wireline logs for 100% oil, gas, and brine saturations and for a 5‐year prediction that was based on a fluid‐flow and production simulation. Analytical AVO curves for simple half‐space models did not match AVO curves extracted from synthetic seismograms computed with a full‐waveform layer‐stack algorithm. Several different amplitude corrections were tried to reduce the AVO curves from the synthetic data to the analytical ones, but, ultimately, none was deemed satisfactory. Instead, AVO change attributes based on relative changes, polarity changes, or ratios were used. Attributes based on the change of AVO gradient were perceived to be most diagnostic of the water flood, but they were also overly sensitive to interference noise and amplitude correction errors. For field data from the study area, a large decrease in intercept magnitude may be the best indicator of the waterfront.


2021 ◽  
Vol 95 (2) ◽  
Author(s):  
Mirjam Bilker-Koivula ◽  
Jaakko Mäkinen ◽  
Hannu Ruotsalainen ◽  
Jyri Näränen ◽  
Timo Saari

AbstractPostglacial rebound in Fennoscandia causes striking trends in gravity measurements of the area. We present time series of absolute gravity data collected between 1976 and 2019 on 12 stations in Finland with different types of instruments. First, we determine the trends at each station and analyse the effect of the instrument types. We estimate, for example, an offset of 6.8 μgal for the JILAg-5 instrument with respect to the FG5-type instruments. Applying the offsets in the trend analysis strengthens the trends being in good agreement with the NKG2016LU_gdot model of gravity change. Trends of seven stations were found robust and were used to analyse the stabilization of the trends in time and to determine the relationship between gravity change rates and land uplift rates as measured with global navigation satellite systems (GNSS) as well as from the NKG2016LU_abs land uplift model. Trends calculated from combined and offset-corrected measurements of JILAg-5- and FG5-type instruments stabilized in 15 to 20 years and at some stations even faster. The trends of FG5-type instrument data alone stabilized generally within 10 years. The ratio between gravity change rates and vertical rates from different data sets yields values between − 0.206 ± 0.017 and − 0.227 ± 0.024 µGal/mm and axis intercept values between 0.248 ± 0.089 and 0.335 ± 0.136 µGal/yr. These values are larger than previous estimates for Fennoscandia.


Geophysics ◽  
2008 ◽  
Vol 73 (1) ◽  
pp. E1-E5 ◽  
Author(s):  
Lev Vernik

Seismic reservoir characterization and pore-pressure prediction projects rely heavily on the accuracy and consistency of sonic logs. Sonic data acquisition in wells with large relative dip is known to suffer from anisotropic effects related to microanisotropy of shales and thin-bed laminations of sand, silt, and shale. Nonetheless, if anisotropy parameters can be related to shale content [Formula: see text] in siliciclastic rocks, then I show that it is straightforward to compute the anisotropy correction to both compressional and shear logs using [Formula: see text] and the formation relative dip angle. The resulting rotated P-wave sonic logs can be used to enhance time-depth ties, velocity to effective stress transforms, and low-frequency models necessary for prestack seismic amplitude variation with offset (AVO) inversion.


2017 ◽  
Vol 7 (1) ◽  
Author(s):  
K. Breili ◽  
R. Hougen ◽  
D. I. Lysaker ◽  
O. C. D. Omang ◽  
B. Tangen

AbstractThe Norwegian Mapping Authority (NMA) has recently established a new gravity laboratory in Ny-Ålesund at Svalbard, Norway. The laboratory consists of three independent pillars and is part of the geodetic core station that is presently under construction at Brandal, approximately 1.5 km north of NMA’s old station. In anticipation of future use of the new gravity laboratory, we present benchmark gravity values, gravity gradients, and final coordinates of all new pillars. Test measurements indicate a higher noise level at Brandal compared to the old station. The increased noise level is attributed to higher sensitivity to wind.We have also investigated possible consequences of moving to Brandal when it comes to the gravitational signal of present-day ice mass changes and ocean tide loading. Plausible models representing ice mass changes at the Svalbard archipelago indicate that the gravitational signal at Brandal may differ from that at the old site with a size detectable with modern gravimeters. Users of gravity data from Ny-Ålesund should, therefore, be cautious if future observations from the new observatory are used to extend the existing gravity record. Due to its lower elevation, Brandal is significantly less sensitive to gravitational ocean tide loading. In the future, Brandal will be the prime site for gravimetry in Ny-Ålesund. This ensures gravity measurements collocated with space geodetic techniques like VLBI, SLR, and GNSS.


2020 ◽  
Vol 55 (3) ◽  
pp. 100-117
Author(s):  
Viktor Szabó ◽  
Dorota Marjańska

AbstractGlobal satellite gravity measurements provide unique information regarding gravity field distribution and its variability on the Earth. The main cause of gravity changes is the mass transportation within the Earth, appearing as, e.g. dynamic fluctuations in hydrology, glaciology, oceanology, meteorology and the lithosphere. This phenomenon has become more comprehensible thanks to the dedicated gravimetric missions such as Gravity Recovery and Climate Experiment (GRACE), Challenging Minisatellite Payload (CHAMP) and Gravity Field and Steady-State Ocean Circulation Explorer (GOCE). From among these missions, GRACE seems to be the most dominating source of gravity data, sharing a unique set of observations from over 15 years. The results of this experiment are often of interest to geodesists and geophysicists due to its high compatibility with the other methods of gravity measurements, especially absolute gravimetry. Direct validation of gravity field solutions is crucial as it can provide conclusions concerning forecasts of subsurface water changes. The aim of this work is to present the issue of selection of filtration parameters for monthly gravity field solutions in RL06 and RL05 releases and then to compare them to a time series of absolute gravimetric data conducted in quasi-monthly measurements in Astro-Geodetic Observatory in Józefosław (Poland). The other purpose of this study is to estimate the accuracy of GRACE temporal solutions in comparison with absolute terrestrial gravimetry data and making an attempt to indicate the significance of differences between solutions using various types of filtration (DDK, Gaussian) from selected research centres.


Author(s):  
K. Jourde ◽  
D. Gibert ◽  
J. Marteau

Abstract. This paper examines how the resolution of small-scale geological density models is improved through the fusion of information provided by gravity measurements and density muon radiographies. Muon radiography aims at determining the density of geological bodies by measuring their screening effect on the natural flux of cosmic muons. Muon radiography essentially works like medical X-ray scan and integrates density information along elongated narrow conical volumes. Gravity measurements are linked to density by a 3-D integration encompassing the whole studied domain. We establish the mathematical expressions of these integration formulas – called acquisition kernels – and derive the resolving kernels that are spatial filters relating the true unknown density structure to the density distribution actually recovered from the available data. The resolving kernels approach allows to quantitatively describe the improvement of the resolution of the density models achieved by merging gravity data and muon radiographies. The method developed in this paper may be used to optimally design the geometry of the field measurements to perform in order to obtain a given spatial resolution pattern of the density model to construct. The resolving kernels derived in the joined muon/gravimetry case indicate that gravity data are almost useless to constrain the density structure in regions sampled by more than two muon tomography acquisitions. Interestingly the resolution in deeper regions not sampled by muon tomography is significantly improved by joining the two techniques. The method is illustrated with examples for La Soufrière of Guadeloupe volcano.


Sign in / Sign up

Export Citation Format

Share Document