Quantitative CO2 monitoring workflow

Author(s):  
Bastien Dupuy ◽  
Anouar Romdhane ◽  
Peder Eliasson

<p>CO<sub>2</sub> storage operators are required to monitor storage safety during injection with a long-term perspective (Ringrose and Meckel, 2019), implying that efficient measurement, monitoring and verification (MMV) plans are of critical importance for the viability of such projects. MMV plans usually include containment, conformance and contingency monitoring. Conformance monitoring is carried out to verify that observations from monitoring data are consistent with predictions from prior reservoir modelling within a given uncertainty range. Quantitative estimates of relevant reservoir parameters (e.g. pore pressure and fluid saturations) are usually derived from geophysical monitoring data (e.g. seismic, electromagnetic and/or gravity data) and potential prior knowledge of the storage reservoir.</p><p>In this work, we describe and apply a two-step strategy combining geophysical and rock physics inversions for quantitative CO<sub>2</sub> monitoring. Bayesian formulations are used to propagate and account for uncertainties in both steps (Dupuy et al., 2017). We apply our workflow to data from the Sleipner CO<sub>2</sub> storage project, located offshore Norway. At Sleipner, the CO<sub>2</sub> has been injected at approx. 1000 m deep, in the high porosity, high permeability Utsira aquifer sandstone since 1996 with an approximate rate of 1 million tonnes per year. We combine seismic full waveform inversion and rock physics inversion to show that 2D spatial distribution of CO<sub>2</sub> saturation can be obtained. Appropriate and calibrated rock physics models need to take into account the way fluid phases are mixed together (uniform to patchy mixing) and the trade-off effects between pore pressure and fluid saturation. For the Sleipner case, we show that the pore pressure build-up can be neglected and that the derived CO<sub>2</sub> saturation distributions mainly depend on P-wave velocities and on the rock physics model. The CO<sub>2</sub> saturation is larger at the top of the reservoir and the mixing tends to be more uniform. These mixing properties are, however, one of the main uncertainties in the inversion. We discuss the added value of a joint rock physics inversion approach, where multi-physics (electromagnetic, seismic, gravimetry), and multi-parameter inversion can be used to reduce the under-determination of the inverse problem and to better discriminate pressure, saturation, and fluid mixing effects.</p><p>Acknowledgements:</p><p>This publication has been produced with support from the NCCS Centre, performed under the Norwegian research program Centres for Environment-friendly Energy Research (FME). The authors acknowledge the following partners for their contributions: Aker Solutions, Ansaldo Energia, CoorsTek Membrane Sciences, Emgs, Equinor, Gassco, Krohne, Larvik Shipping, Lundin, Norcem, Norwegian Oil and Gas, Quad Geometrics, Total, Vår Energi, and the Research Council of Norway (257579/E20).</p><p>References:</p><p>Dupuy, B., Romdhane, A., Eliasson, P., Querendez, E., Yan, H., Torres, V. A., and Ghaderi, A. (2017). Quantitative seismic characterization of CO<sub>2</sub> at the Sleipner storage site, North Sea. Interpretation, 5(4):SS23–SS42.</p><p>Ringrose, P. S. and Meckel, T. A. (2019). Maturing global CO<sub>2</sub> storage resources on offshore continental margins to achieve 2DS emissions reductions. Scientific Reports, 9(1):1–10.</p>

Geophysics ◽  
2003 ◽  
Vol 68 (5) ◽  
pp. 1569-1579 ◽  
Author(s):  
José M. Carcione ◽  
Hans B. Helle ◽  
Nam H. Pham ◽  
Tommy Toverud

A method is used to obtain pore pressure in shaly sandstones based upon an acoustic model for seismic velocity versus clay content and effective pressure. Calibration of the model requires log data—porosity, clay content, and sonic velocities—to obtain the dry‐rock moduli and the effective stress coefficients as a function of depth and pore pressure. The seismic P‐wave velocity, derived from reflection tomography, is fitted to the theoretical velocity by using pore pressure as the fitting parameter. This approach, based on a rock‐physics model, is an improvement over existing pore‐pressure prediction methods, which mainly rely on empirical relations between velocity and pressure. The method is applied to the Tune field in the Viking Graben sedimentary basin of the North Sea. We have obtained a high‐resolution velocity map that reveals the sensitivity to pore pressure and fluid saturation in the Tarbert reservoir. The velocity map of the Tarbert reservoir and the inverted pressure distribution agree with the structural features of the Tarbert Formation and its known pressure compartments.


Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. MR271-MR283
Author(s):  
Ismael Himar Falcon-Suarez ◽  
Laurence North ◽  
Ben Callow ◽  
Gaye Bayrakci ◽  
Jon Bull ◽  
...  

Seismic and electromagnetic properties generally are anisotropic, depending on the microscale rock fabric and the macroscale stress field. We have assessed the stress-dependent anisotropy of poorly consolidated (porosity of approximately 0.35) sandstones (broadly representative of shallow reservoirs) experimentally, combining ultrasonic (0.6 MHz P-wave velocity, [Formula: see text], and attenuation [Formula: see text]) and electrical resistivity measurements. We used three cores from an outcrop sandstone sample extracted at 0°, 45°, and 90° angles with respect to the visible geologic bedding plane and subjected them to unloading/loading cycles with variations of the confining (20–35 MPa) and pore (2–17 MPa) pressures. Our results indicate that stress field orientation, loading history, rock fabric, and the measurement scale all affect the elastic and electrical anisotropies. Strong linear correlations ([Formula: see text]) between [Formula: see text], [Formula: see text], and resistivity in the three considered directions suggest that the stress orientation similarly affects the elastic and electrical properties of poorly consolidated, high-porosity (shallow) sandstone reservoirs. However, resistivity is more sensitive to pore-pressure changes (effective stress coefficients [Formula: see text]), whereas P-wave properties provide simultaneous information about the confining (from [Formula: see text], with n slightly less than 1) and pore pressure (from [Formula: see text], with n slightly greater than 1) variations. We found n is also anisotropic for the three measured properties because a more intense and rapid grain rearrangement occurs when the stress field changes result from oblique stress orientations with respect to rock layering. Altogether, our results highlighted the potential of joint elastic-electrical stress-dependent anisotropy assessments to enhance the geomechanical interpretation of reservoirs during production or injection activities.


Geophysics ◽  
1998 ◽  
Vol 63 (5) ◽  
pp. 1604-1617 ◽  
Author(s):  
Zhijing Wang ◽  
Michael E. Cates ◽  
Robert T. Langan

A carbon dioxide (CO2) injection pilot project is underway in Section 205 of the McElroy field, West Texas. High‐resolution crosswell seismic imaging surveys were conducted before and after CO2 flooding to monitor the CO2 flood process and map the flooded zones. The velocity changes observed by these time‐lapse surveys are typically on the order of −6%, with maximum values on the order of −10% in the vicinity of the injection well. These values generally agree with laboratory measurements if the effects of changing pore pressure are included. The observed dramatic compressional ([Formula: see text]) and shear ([Formula: see text]) velocity changes are considerably greater than we had initially predicted using the Gassmann (1951) fluid substitution analysis (Nolen‐Hoeksema et al., 1995) because we had assumed reservoir pressure would not change from survey to survey. However, the post‐CO2 reservoir pore fluid pressure was substantially higher than the original pore pressure. In addition, our original petrophysical data for dry and brine‐saturated reservoir rocks did not cover the range of pressures actually seen in the field. Therefore, we undertook a rock physics study of CO2 flooding in the laboratory, under the simulated McElroy pressures and temperature. Our results show that the combined effects of pore pressure buildup and fluid substitution caused by CO2 flooding make it petrophysically feasible to monitor the CO2 flood process and to map the flooded zones seismically. The measured data show that [Formula: see text] decreases from a minimum 3.0% to as high as 10.9%, while [Formula: see text] decreases from 3.3% to 9.5% as the reservoir rocks are flooded with CO2 under in‐situ conditions. Such [Formula: see text] and [Formula: see text] decreases, even if averaged over all the samples measured, are probably detectable by either crosswell or high‐resolution surface seismic imaging technologies. Our results show [Formula: see text] is sensitive to both the CO2 saturation and the pore pressure increase, but [Formula: see text] is particularly sensitive to the pore pressure increase. As a result, the combined [Formula: see text] and [Formula: see text] changes caused by the CO2 injection may be used, at least semiquantitatively, to separate CO2‐flooded zones with pore pressure buildup from those regions without pore pressure buildup or to separate CO2 zones from pressured‐up, non‐CO2 zones. Our laboratory results show that the largest [Formula: see text] and [Formula: see text] changes caused by CO2 injection are associated with high‐porosity, high‐permeability rocks. In other words, CO2 flooding and pore pressure buildup decrease [Formula: see text] and [Formula: see text] more in high‐porosity, high‐permeability samples. Therefore, it may be possible to delineate such high‐porosity, high‐permeability streaks seismically in situ. If the streaks are thick enough compared to seismic resolution, they can be identified by the larger [Formula: see text] or [Formula: see text] changes.


2021 ◽  
pp. 1-47
Author(s):  
Chao Li ◽  
Peng Hu ◽  
Jing Ba ◽  
José M. Carcione ◽  
Tianwen Hu ◽  
...  

Tight-gas sandstone reservoirs of the Ordos Basin of China are characterized by high rock-fragment content, dissimilar pore types and a random distribution of fluids, leading to strong local heterogeneity. We model the seismic properties of these sandstones with the double-double porosity (DDP) theory, which considers water saturation, porosity and the frame characteristics. A generalized seismic wavelet is used to fit the real wavelet and the peak frequency-shift method is combined with the generalized S-transform to estimate attenuation. Then, we establish rock-physics templates (RPTs) based on P-wave attenuation and impedance. We use the log data and related seismic traces to calibrate the RPTs and generate a 3D volume of rock-physics attributes for the quantitative prediction of saturation and porosity. The predicted values are in good agreement with the actual gas production reports, indicating that the method can be effectively applied to heterogeneous tight-gas sandstone reservoirs.


Geophysics ◽  
2021 ◽  
pp. 1-54
Author(s):  
Yijun Wei ◽  
Jing Ba ◽  
José M. Carcione ◽  
Li-Yun Fu ◽  
Mengqiang Pang ◽  
...  

Ultra-deep carbonate reservoirs have high temperatures and pressures, complex pressure/tectonic stress settings and pore structures. These conditions make their seismic detection and characterization difficult, particularly if the signal-to-noise ratio is low, as it is the case in most situations. Moreover, the high risk of deep-drilling exploration makes it impractical to carry out normal logging operations. We propose a temperature-differential pressure-porosity (TPP) rock-physics model based on the Biot-Rayleigh poroelasticity theory to describe the wave response of the reservoir. A preliminary analysis shows that temperature, pressure and porosity are well correlated with wave velocity and attenuation. On the basis of this theory, we built 3D rock-physics templates that account for the effects of TPP on the P-wave impedance, VP/ VS ratio and attenuation. The templates are calibrated with laboratory, well-log and seismic data of the S area (Shuntuoguole uplift, Tarim Basin, Xinjiang, China). Then, the template is used to obtain the properties of the reservoir at seismic frequencies. The predicted results are consistent with the field reports, high temperature, low differential pressure and high porosity, indicating high production rates. The methodology will be useful for the hydrocarbon exploration in ultra-deep carbonate reservoirs.


Geophysics ◽  
2020 ◽  
pp. 1-55
Author(s):  
Wolfgang Weinzierl ◽  
Bernd Wiese

Determining saturation and pore pressure is relevant for hydrocarbon production as well as natural gas and CO2 storage. In this context seismic methods provide spatially distributed data used to determine gas and fluid migration. A method is developed that allows to determine saturation and reservoir pressure from seismic data, more precisely from rock physical attributes that are velocity, attenuation and density. Two rock physical models based on Hertz-Mindlin-Gassmann and Biot-Gassmann are developed. Both generate poroelastic attributes from pore pressure, gas saturation and other rock-physical parameters. The rock physical models are inverted with deep neural networks to derive e.g. saturation, pore pressure and porosity from rock physical attributes. The method is demonstrated with a 65 m deep unconsolidated high porosity reservoir at the Svelvik ridge, Norway. Tests for the most suitable structure of the neural network are carried out. Saturation and pressure can be meaningfully determined under condition of a gas-free baseline with known pressure and data from an accurate seismic campaign, preferably cross-well seismic. Including seismic attenuation increases the accuracy. The training requires hours, predictions just a few seconds, allowing for rapid interpretation of seismic results.


2017 ◽  
Vol 5 (4) ◽  
pp. SS23-SS42 ◽  
Author(s):  
Bastien Dupuy ◽  
Anouar Romdhane ◽  
Peder Eliasson ◽  
Etor Querendez ◽  
Hong Yan ◽  
...  

Reliable quantification of carbon dioxide ([Formula: see text]) properties and saturation is crucial in the monitoring of [Formula: see text] underground storage projects. We have focused on quantitative seismic characterization of [Formula: see text] at the Sleipner storage pilot site. We evaluate a methodology combining high-resolution seismic waveform tomography, with uncertainty quantification and rock physics inversion. We use full-waveform inversion (FWI) to provide high-resolution estimates of P-wave velocity [Formula: see text] and perform an evaluation of the reliability of the derived model based on posterior covariance matrix analysis. To get realistic estimates of [Formula: see text] saturation, we implement advanced rock physics models taking into account effective fluid phase theory and patchy saturation. We determine through sensitivity tests that the estimation of [Formula: see text] saturation is possible even when using only the P-wave velocity as input. After a characterization of rock frame properties based on log data prior to the [Formula: see text] injection at Sleipner, we apply our two-step methodology. The FWI result provides clear indications of the injected [Formula: see text] plume being observed as low-velocity zones corresponding to thin [Formula: see text] filled layers. Several tests, varying the rock physics model and [Formula: see text] properties, are then performed to estimate [Formula: see text] saturation. The results suggest saturations reaching 30%–35% in the thin sand layers and up to 75% when patchy mixing is considered. We have carried out a joint estimation of saturation with distribution type and, even if the inversion is not well-constrained due to limited input data, we conclude that the [Formula: see text] has an intermediate pattern between uniform and patchy mixing, which leads to saturation levels of approximately [Formula: see text]. It is worth noting that the 2D section used in this work is located 533 m east of the injection point. We also conclude that the joint estimation of [Formula: see text] properties with saturation is not crucial and consequently that knowing the pressure and temperature state of the reservoir does not prevent reliable estimation of [Formula: see text] saturation.


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.


Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. B229-B240 ◽  
Author(s):  
Rajive Kumar ◽  
Prashant Bansal ◽  
Bader S. Al-Mal ◽  
Sagnik Dasgupta ◽  
Colin Sayers ◽  
...  

Optimization of production from unconventional reservoirs requires estimates of reservoir properties such as porosity, total organic carbon (TOC) content, clay content, fluid saturation, and fracture intensity. The porosity and TOC content help to determine reservoir quality, and the natural fracture intensity provides information important for the completion strategy. Because shale reservoirs display intrinsic anisotropy due to layering and the partial alignment of clay minerals and kerogen with the bedding plane, the minimum acceptable representation of the anisotropy of naturally fractured shale-gas reservoirs is orthotropy, in which a set of vertical compliant fractures is embedded in a vertical transverse isotropic (VTI) background medium. Full-azimuth seismic data are required to characterize such reservoirs and to invert for the anisotropic elastic properties. Orthotropic inversion uses azimuthally sectored seismic data stacked according to the incident angle. Even for high-fold acquisition, this azimuth/angle grouping can result in low-fold angle stacks. Orthotropic amplitude-variation-with-offset-and-azimuth (AVOAz) inversion requires seismic preconditioning techniques that ensure proper primary amplitude preservation, noise attenuation, and data alignment, and a workflow implemented for the construction of an orthotropic rock-physics model. This model integrates well and core data to estimate reservoir properties using the results of the AVOAz inversion. The seismic inversion results include the P- and S-impedance and parameters quantifying the azimuthal anisotropy. The rock model assumes a VTI kerogen-rich layer, containing aligned vertical fractures, and it uses prestack orthotropic AVOAz inversion results to predict porosity, TOC, and fracture intensity.


2012 ◽  
Vol 2012 ◽  
pp. 1-17 ◽  
Author(s):  
Aaron V. Wandler ◽  
Thomas L. Davis ◽  
Paritosh K. Singh

In mature oil fields undergoing enhanced oil recovery methods, such as CO2injection, monitoring the reservoir changes becomes important. To understand how reservoir changes influence compressional wave (P) and shear wave (S) velocities, we conducted laboratory core experiments on five core samples taken from the Morrow A sandstone at Postle Field, Oklahoma. The laboratory experiments measured P- and S-wave velocities as a function of confining pressure, pore pressure, and fluid type (which included CO2in the gas and supercritical phase). P-wave velocity shows a response that is sensitive to both pore pressure and fluid saturation. However, S-wave velocity is primarily sensitive to changes in pore pressure. We use the fluid and pore pressure response measured from the core samples to modify velocity well logs through a log facies model correlation. The modified well logs simulate the brine- and CO2-saturated cases at minimum and maximum reservoir pressure and are inputs for full waveform seismic modeling. Modeling shows how P- and S-waves have a different time-lapse amplitude response with offset. The results from the laboratory experiments and modeling show the advantages of combining P- and S-wave attributes in recognizing the mechanism responsible for time-lapse changes due to CO2injection.


Sign in / Sign up

Export Citation Format

Share Document