scholarly journals Monitoring geological storage of CO2 using a new rock physics model

2022 ◽  
Vol 12 (1) ◽  
Author(s):  
Manzar Fawad ◽  
Nazmul Haque Mondol

AbstractTo mitigate the global warming crisis, one of the effective ways is to capture CO2 at an emitting source and inject it underground in saline aquifers, depleted oil and gas reservoirs, or in coal beds. This process is known as carbon capture and storage (CCS). With CCS, CO2 is considered a waste product that has to be disposed of properly, like sewage and other pollutants. While and after CO2 injection, monitoring of the CO2 storage site is necessary to observe CO2 plume movement and detect potential leakage. For CO2 monitoring, various physical property changes are employed to delineate the plume area and migration pathways with their pros and cons. We introduce a new rock physics model to facilitate the time-lapse estimation of CO2 saturation and possible pressure changes within a CO2 storage reservoir based on physical properties obtained from the prestack seismic inversion. We demonstrate that the CO2 plume delineation, saturation, and pressure changes estimations are possible using a combination of Acoustic Impedance (AI) and P- to S-wave velocity ratio (Vp/Vs) inverted from time-lapse or four-dimensional (4D) seismic. We assumed a scenario over a period of 40 years comprising an initial 25 year injection period. Our results show that monitoring the CO2 plume in terms of extent and saturation can be carried out using our rock physics-derived method. The suggested method, without going into the elastic moduli level, handles the elastic property cubes, which are commonly obtained from the prestack seismic inversion. Pressure changes quantification is also possible within un-cemented sands; however, the stress/cementation coefficient in our proposed model needs further study to relate that with effective stress in various types of sandstones. The three-dimensional (3D) seismic usually covers the area from the reservoir's base to the surface making it possible to detect the CO2 plume's lateral and vertical migration. However, the comparatively low resolution of seismic, the inversion uncertainties, lateral mineral, and shale property variations are some limitations, which warrant consideration. This method can also be applied for the exploration and monitoring of hydrocarbon production.

2006 ◽  
Author(s):  
Kyle Spikes ◽  
Jack Dvorkin ◽  
Gary Mavko

2018 ◽  
Vol 6 (4) ◽  
pp. SM1-SM8 ◽  
Author(s):  
Tingting Zhang ◽  
Yuefeng Sun

Fractured zones in deeply buried carbonate hills are important because they often have better permeability resulting in prolific production than similar low-porosity rocks. Nevertheless, their detection poses great challenge to conventional seismic inversion methods because they are mostly low in acoustic impedance and bulk modulus, hardly distinguishable from high-porosity zones or mudstones. A proxy parameter of pore structure defined in a rock-physics model, the so-called Sun model, has been used for delineating fractured zones in which the pore structure parameter is relatively high, whereas the porosity is low in general. Simultaneous seismic inversion of the pore structure parameter and porosity proves to be difficult and nontrivial in practice. Although the pore structure parameter is well-defined at locations where density, P-, and S-velocity are known from logs, estimation of P- and S-velocity information, especially density information from prestack seismic data is rather challenging. A three-step iterative inversion method, which uses acoustic, gradient, and elastic impedance from angle-stacked seismic data as input to the rock-physics model for calculating porosity and bulk and shear pore structure parameters simultaneously, is proposed and implemented to solve this problem. The methodology is successfully tested with well logs and seismic data from a deeply buried carbonate hill in the Bohai Bay Basin, China.


2019 ◽  
Vol 38 (5) ◽  
pp. 358-365 ◽  
Author(s):  
Colin M. Sayers ◽  
Sagnik Dasgupta

This paper presents a predictive rock-physics model for unconventional shale reservoirs based on an extended Maxwell scheme. This model accounts for intrinsic anisotropy of rock matrix and heterogeneities and shape-induced anisotropy arising because the dimensions of kerogen inclusions and pores are larger parallel to the bedding plane than perpendicular to this plane. The model relates the results of seismic amplitude variation with offset inversion, such as P- and S-impedance, to the composition of the rock and enables identification of rock classes such as calcareous, argillaceous, siliceous, and mixed shales. This allows the choice of locations with the best potential for economic production of hydrocarbons. While this can be done using well data, prestack inversion of seismic P-wave data allows identification of the best locations before the wells are drilled. The results clearly show the ambiguity in rock classification obtained using poststack inversion of P-wave seismic data and demonstrate the need for prestack seismic inversion. The model provides estimates of formation anisotropy, as required for accurate determination of P- and S-impedance, and shows that anisotropy is a function not only of clay content but also other components of the rock as well as the aspect ratio of kerogen and pores. Estimates of minimum horizontal stress based on the model demonstrate the need to identify rock class and estimate anisotropy to determine the location of any stress barriers that may inhibit hydraulic fracture growth.


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.


2020 ◽  
Vol 223 (3) ◽  
pp. 1610-1629
Author(s):  
Gil Hetz ◽  
Akhil Datta-Gupta ◽  
Justyna K Przybysz-Jarnut ◽  
Jorge L Lopez ◽  
D W Vasco

SUMMARY Our limited knowledge of the relationship between changes in the state of an aquifer or reservoir and the corresponding changes in the elastic moduli, that is the rock physics model, hampers the effective use of time-lapse seismic observations for estimating flow properties within the Earth. A central problem is the complicated dependence of the magnitude of time-lapse changes on the saturation, pressure, and temperature changes within an aquifer or reservoir. We describe an inversion methodology for reservoir characterization that uses onset times, the calendar time of the change in seismic attributes, rather than the magnitude of the changes. We find that onset times are much less sensitive than magnitudes to the rock physics model used to relate time-lapse observations to changes in saturation, temperature and fluid pressure. We apply the inversion scheme to observations from daily monitoring of enhanced oil recovery at the Peace River field in Canada. An array of 1492 buried hydrophones record seismic signals from 49 buried sources. Time-shifts for elastic waves traversing the reservoir are extracted from the daily time-lapse cubes. In our analysis 175 images of time-shifts are transformed into a single map of onset times, leading to a substantial reduction in the volume of data. These observations are used in conjunction with bottom hole pressure data to infer the initial conditions prior to the injection, and to update the reservoir permeability model. The combination of a global and local inversion scheme produces a collection of reservoir models that are best described by three clusters. The updated model leads to a nearly 70 percent reduction in seismic data misfit. The final set of solutions successfully predict the observed normalized pressure history during the soak and flow-back into the wells between 82 and 175 days into the cyclic steaming operation.


Sign in / Sign up

Export Citation Format

Share Document