Integrated geophysical and hydromechanical assessment for CO2 storage: shallow low permeable reservoir sandstones

Geological reservoirs can be structurally complex and can respond to CO2 injection both geochemically and geomechanically. Hence, predicting reservoir formation behaviour in response to CO2 injection and assessing the resulting hazards are important prerequisites for safe geological CO2 storage. This requires a detailed study of thermal-hydro-mechanical-chemical coupled phenomena that can be triggered in the reservoir formation, most readily achieved through laboratory simulations of CO2 injection into typical reservoir formations. Here, we present the first results from a new experimental apparatus of a steady-state drainage flooding test conducted through a synthetic sandstone sample, simulating real CO2 storage reservoir conditions in a shallow (?1 km), low permeability ?1mD, 26% porosity sandstone formation. The injected pore fluid comprised brine with CO2 saturation increasing in steps of 20% brine/CO2 partial flow rates up to 100% CO2 flow. At each pore fluid stage, an unload/loading cycle of effective pressure was imposed to study the response of the rock under different geomechanical scenarios. The monitoring included axial strains and relative permeability in a continuous mode (hydromechanical assessment), and related geophysical signatures (ultrasonic P-wave and S-wave velocities Vp and Vs, and attenuations Qp?1 and Qs?1, respectively, and electrical resistivity). On average, the results showed Vp and Vs dropped ?7% and ?4% respectively during the test, whereas Qp?1 increased ?55% and Qs?1 decreased by ?25%. From the electrical resistivity data, we estimated a maximum CO2 saturation of ?0.5, whereas relative permeability curves were adjusted for both fluids. Comparing the experimental results to theoretical predictions, we found that Gassmann's equations explain Vp at high and very low CO2 saturations, whereas bulk modulus yields results consistent with White and Dutta–Odé model predictions. This is interpreted as a heterogeneous distribution of the two pore fluid phases, corroborated by electrical resistivity tomography images. The integration of laboratory geophysical and hydromechanical observations on representative shallow low-permeable sandstone reservoir allowed us to distinguish between pure geomechanical responses and those associated with the pore fluid distribution. This is a key aspect in understanding CO2 injection effects in deep geological reservoirs associated with carbon capture and storage practices

Water saturation effects on P-wave anisotropy in synthetic sandstone with aligned fractures

The seismic properties of rocks are known to be sensitive to partial liquid or gas saturation, and to aligned fractures. P-wave anisotropy is widely used for fracture characterization and is known to be sensitive to the saturating fluid. However, studies combining the effect of multiphase saturation and aligned fractures are limited even though such conditions are common in the subsurface. An understanding of the effects of partial liquid or gas saturation on P-wave anisotropy could help improve seismic characterization of fractured, gas bearing reservoirs. Using octagonal-shaped synthetic sandstone samples, one containing aligned penny-shaped fractures and the other without fractures, we examined the influence of water saturation on P-wave anisotropy in fractured rocks. In the fractured rock, the saturation related stiffening effect at higher water saturation values is larger in the direction across the fractures than along the fractures. Consequently, the anisotropy parameter ‘?’ decreases as a result of this fluid stiffening effect. These effects are frequency dependent as a result of wave-induced fluid flow mechanisms. Our observations can be explained by combining a frequency-dependent fractured rock model and a frequency-dependent partial saturation model

Modelling ultrasonic laboratory measurements of the saturation dependence of elastic modulus: new insights and implications for wave propagation mechanisms

Seismic time-lapse techniques are a valuable tool used to estimate the mobilization and distribution of stored CO2 in depleted reservoirs. The success of these techniques depends on knowing the seismic properties of partially saturated rocks with accuracy. It is commonplace to use controlled laboratory-scale experiments to determine how the fluid content impacts on their properties. In this work, we measure the ultrasonic P- and S-wave velocities of a set of synthetic sandstones of about 30% porosity. Using an accurate method, we span the entire saturation range of an air-water system. We show that the rocks’ elastic behaviour is consistent with patchy saturation and squirt flow models but observe a discontinuity at around 90% gas saturation which can be interpreted in two very different ways. In one interpretation, the responsible mechanism is frequency-dependent squirt-flow that occurs in narrow pores that are preferentially saturated. An equally plausible mechanism is the change of the mobile fluid in the pores once they are wetted. Extrapolated to seismic frequencies, our results imply that the seismic properties of rocks may be affected by the wetting effect with an impact on the interpretation of field data but would potentially be unaffected by the squirt flow effect. This provides strong motivation to conduct laboratory-scale experiments with partially saturated samples at lower frequency or, ideally, a range of frequencies in the seismo-acoustic range

Effects of aligned fractures on the response of velocity and attenuation ratios to water saturation variation: A laboratory study using synthetic sandstones

P-wave-to-S-wave ratios are important seismic characterization attributes. Velocity ratios are sensitive to the petrophysical properties of rocks and to the presence of gas. Attenuation ratios have also been shown to be sensitive to the presence of partial liquid/gas saturation. The relationship between liquid/gas saturation and P-wave and S-wave ratios has been used to distinguish gas-saturated rocks from liquid-saturated rocks. Aligned fractures are common in the Earth's crust and cause seismic anisotropy and shear wave splitting. However, most existing relationships between partial gas/liquid saturation and P-wave and S-wave ratios are for non-fractured rocks. We present experimental results comparing the effects of changing water saturation on Qs/Qp versus Vp/Vs ratios between a non-fractured rock and one containing fractures aligned parallel to wave propagation direction. We also study the effects of aligned fractures on the response of Vp/Vs to changing water saturation using synthetic fractured sandstones with fractures aligned at 45o and parallel to the wave propagation direction. The results suggest that aligned fractures could have significant effects on the observed trends, some of which may not be obvious. Fractures aligned parallel to wave propagation could change the response of Qs/Qp versus Vp/Vs ratios to water saturation from previously reported trends. Shear wave splitting due to the presence of aligned fractures results in two velocity ratios (Vp/Vs1 and Vp/Vs2). The fluid independence of shear wave splitting for fractures aligned parallel to wave propagation direction means the difference between Vp/Vs1 and Vp/Vs2 is independent of water saturation. For fractures aligned at oblique angles, shear wave splitting can be sensitive to water saturation and consequently be frequency dependent, which can lead to fluid and frequency-dependent differences between Vp/Vs1 and Vp/Vs2. The effect of aligned fractures on Vp/Vs ratios not only depends on the fracture effects on both P-wave and S-wave velocities but also on the effects of water saturation distribution on the rock and fracture stiffness, and hence on the P-wave and S-wave velocities. As such, these effects can be frequency dependent due to wave-induced fluid flow. A simple modelling study combining a frequency-dependent fractured rock model, and a frequency-dependent partial saturation model was used to gain valuable interpretations of our experimental observations and possible implications, which would be useful for field seismic data interpretation

Theoretical Derivation of a Brie-like Fluid Mixing Law

Prediction of the velocity of acoustic waves in partially saturated rocks is very important in geophysical applications. The need to accurately predict acoustic velocities has resulted in a widespread popularity of Brie's effective fluid mixing law. This empirical model together with Gassmann's formula are used routinely in fluid substitution problems in petroleum geophysics and seismic monitoring of carbon capture and storage. Most attempts to justify Brie's model have been focused on interpretation in terms of patchy saturation models and attaching meaning to the Brie parameter in terms of the patch size. In this paper, using a microstructural description of the rock and a parameter relating to capillary pressure, we calculate an effective fluid modulus that is very similar to Brie's law. The fluid mixing law we propose is independent of frequency and has a solid theoretical foundation. This proposed law produces analytically harmonic and arithmetic averaging at the endpoints. Our results indicate that Brie-like behaviour may not necessarily be related to frequency- and patch-size- dependent phenomena

Comparison of stress-dependent geophysical, hydraulic and mechanical properties of synthetic and natural sandstones for reservoir characterization and monitoring studies

Synthetic rock samples can offer advantages over natural rock samples when used for laboratory rock physical properties studies, provided their success as natural analogues is well understood. The ability of synthetic rocks to mimic the natural stress dependency of elastic wave, electrical and fluid transport properties is of primary interest. Hence, we compare a consistent set of laboratory multi‐physics measurements obtained on four quartz sandstone samples (porosity range 20–25%) comprising two synthetic and two natural (Berea and Corvio) samples, the latter used extensively as standards in rock physics research. We measured simultaneously ultrasonic (P‐ and S‐wave) velocity and attenuation, electrical resistivity, permeability and axial and radial strains over a wide range of differential pressure (confining stress 15–50 MPa; pore pressure 5–10 MPa) on the four brine saturated samples. Despite some obvious physical discrepancies caused by the synthetic manufacturing process, such as silica cementation and anisotropy, the results show only small differences in stress dependency between the synthetic and natural sandstones for all measured parameters. Stress dependency analysis of the dry samples using an isotropic effective medium model of spheroidal pores and penny‐shaped cracks, together with a granular cohesion model, provide evidence of crack closure mechanisms in the natural sandstones, seen to a much lesser extent in the synthetic sandstones. The smaller grain size, greater cement content, and cementation under oedometric conditions particularly affect the fluid transport properties of the synthetic sandstones, resulting in lower permeability and higher electrical resistivity for a similar porosity. The effective stress coefficients, determined for each parameter, are in agreement with data reported in the literature. Our results for the particular synthetic materials that were tested suggest that synthetic sandstones can serve as good proxies for natural sandstones for studies of elastic and mechanical properties, but should be used with care for transport properties studies

Estimating gas saturation in a thin layer by using frequency-dependent amplitude versus offset modelling

Various models have been proposed to link partial gas saturation to seismic attenuation and dispersion, suggesting that the reflection coefficient should be frequency-dependent in many cases of practical importance. Previous approaches to studying this phenomenon have typically been limited to single interface models. Here we propose a modelling technique which allows us to incorporate frequency-dependent reflectivity into convolutional modelling. With this modelling framework, seismic data can be synthesized from well logs of velocity, density, porosity and water saturation. This forward modelling could act as a basis for inversion schemes aimed at recovering gas saturation variations with depth. We present a Bayesian inversion scheme for a simple thin layer case and a particular rock physics model, and show that although the method is very sensitive to prior information and constrains, gas saturation and layer thickness can both theoretically be estimated in the case of interfering reflections

CO2-brine flow-through on an Utsira Sand core sample: Experimental and modelling. Implications for the Sleipner storage field

Sleipner (North Sea) is the world’s first commercial-scale carbon capture and storage (CCS) project, active since 1996, with ∼17 million tonnes of CO2 stored. The main reservoir, Utsira Sand, constitutes an ideal host formation of exceptionally high porosity-permeability and large lateral extent. However, the extensive seismic time-lapse, gravity and electromagnetic monitoring surveys deployed at Sleipner have not been well-supported by laboratory measurements. Here, we investigate the geophysical and geomechanical response of an Utsira core sample for the first time, using controlled inflation/depletion cycles at variable CO2-to-brine fractional flow rates. Ultrasonic P-wave velocities and attenuations are measured together with electrical resistivity (converted into CO2-saturation), along with continuous axial and radial strain monitoring. Ultrasonic velocity and attenuation data were simultaneously inverted and results extrapolated to field-scale seismic-frequencies using a new rock physics theory, which combines patchy fluid distribution and squirt flow effects. It provides a velocity-saturation relationship of practical importance to CO2 plume monitoring. Furthermore, by combining ultrasonic and deformation data, we report empirical relations between pore pressure changes and geomechanical effects in the reservoir, for different saturation ranges. Our dataset complements and constrains existing geophysical monitoring surveys at Sleipner and, more generally, improves the understanding of shallow weakly-cemented sand reservoirs