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

Pressure-varying CO 2 distribution affects the ultrasonic velocities of synthetic sandstones

We performed a novel experiment in which three synthetic sandstones – manufactured using a common method but having different porosities – were saturated with brine and progressively flooded with CO2 under constant confining pressure. The fluid pressure was varied around the critical pressure of CO2 and repeated measurements were made of resistivity, in order to assess the saturation, and elastic wave velocity during the flood. The measured saturated bulk moduli were higher than those predicted by the Gassmann–Wood theory, but were consistent with behaviour described by a recently derived poroelastic model which combines “patch” and “squirt” effects. Measurements on two of the samples followed a patch-based model while those on the highest porosity sample showed evidence of squirt-flow behaviour. Our analysis suggests that the appropriate fluid mixing law is pressure dependent, which is consistent with the notion that the effective patch size decreases as fluid pressure is increased. We derive simple empirical models for the patch dependence from fluid pressure which may be used in seismic modelling and interpretation exercises relevant to monitoring of CO2 injection

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

Controlled-source electromagnetic and seismic delineation of sub-seafloor fluid flow structures in a gas hydrate province, offshore Norway

Deep sea pockmarks underlain by chimney-like or pipe structures that contain methane hydrate are abundant along the Norwegian continental margin. In such hydrate provinces the interaction between hydrate formation and fluid flow has significance for benthic ecosystems and possibly climate change. The Nyegga region, situated on the western Norwegian continental slope, is characterized by an extensive pockmark field known to accommodate substantial methane gas hydrate deposits. The aim of this study is to detect and delineate both the gas hydrate and free gas reservoirs at one of Nyegga's pockmarks. In 2012, a marine controlled-source electromagnetic (CSEM) survey was performed at a pockmark in this region, where high-resolution three-dimensional seismic data were previously collected in 2006. Two-dimensional CSEM inversions were computed using the data acquired by ocean bottom electrical field receivers. Our results, derived from unconstrained and seismically constrained CSEM inversions, suggest the presence of two distinctive resistivity anomalies beneath the pockmark: a shallow vertical anomaly at the underlying pipe structure, likely due to gas hydrate accumulation, and a laterally extensive anomaly attributed to a free gas zone below the base of the gas hydrate stability zone. This work contributes to a robust characterization of gas hydrate deposits within sub-seafloor fluid flow pipe structures

Elastic and electrical properties and permeability of serpentinites from Atlantis Massif, Mid-Atlantic Ridge

Serpentinized peridotites co-exist with mafic rocks in a variety of marine environments including subduction zones, continental rifts and mid-ocean ridges. Remote geophysical methods are crucial to distinguish between them and improve the understanding of the tectonic, magmatic and metamorphic history of the oceanic crust. But, serpentinite peridotites exhibit a wide range of physical properties that complicate such a distinction. We analyzed the ultrasonic P- and S-wave velocities (Vp, Vs) and their respective attenuation (Qp−1, Qs−1), electrical resistivity and permeability of four serpentinized peridotite samples from the southern wall of the Atlantis Massif, Mid-Atlantic Ridge, collected during International Ocean Discovery Program (IODP) Expedition 357. The measurements were taken over a range of loading-unloading stress paths (5 - 45 MPa), using ∼1.7 cm length, 5 cm diameter samples horizontally extracted from the original cores drilled on the seafloor. The measured parameters showed variable degrees of stress dependence, but followed similar trends. Vp, Vs, resistivity and permeability show good inter-correlations, while relationships that included Qp−1 and Qs−1 are less clear. Resistivity showed high contrast between highly serpentinized ultramafic matrix (> 50 Ω m) and mechanically/geochemically altered (magmatic/hydrothermal-driven alteration) domains (< 20 Ω m). This information together with the elastic constants (Vp/Vs ratio and bulk moduli) of the samples allowed us to infer useful information about the degree of serpentinization and the alteration state of the rock, contrasted by petrographic analysis. This study shows the potential of combining seismic techniques and controlled source electromagnetic surveys for understanding tectono-magmatic processes and fluid pathways in hydrothermal systems

The influence of gas hydrate morphology on reservoir permeability and geophysical shear wave remote sensing

We show that direct estimates of the permeability of hydrate-bearing geological formations are possible from remote measurements of shear wave velocity (Vs) and attenuation (Qs−1). We measured Vs, Qs−1 and electrical resistivity at time intervals during methane hydrate formation in Berea sandstone using a laboratory ultrasonic pulse-echo system. We observed that Vs and Qs−1 both increase with hydrate saturation Sh, with two peaks in Qs−1 at hydrate saturations of around 6% and 20% that correspond to changes in gradient of Vs. We implemented changes in permeability with hydrate saturation into well-known Biot-type poro-elastic models for two- and three-phases for low (Sh 12%) hydrate saturations respectively. By accounting for changes in permeability linked to hydrate morphology, the models were able to describe the Vs and Qs−1 observations. We found that the first Qs−1 peak is caused by a reduction of permeability during hydrate formation associated with a transition from pore-floating to pore-bridging hydrate morphology; similarly, the second Qs−1 peak is caused by a permeability reduction associated with a transition from pore-bridging hydrate morphology to an interlocking network of hydrate in the pores. We inverted for permeability using our poro-elastic models from Vs and Qs−1. This inverted permeability agrees with permeability obtained independently from electrical resistivity. We demonstrate a good match of our models to shear wave data at 200 Hz and 2 kHz frequencies from the literature, indicating the general applicability of the models

Experimental assessment of the stress-sensitivity of combined elastic and electrical anisotropy in shallow reservoir sandstones

Seismic and electromagnetic properties are generally 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, VP, and attenuation 1/QP) 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 (R2 > 0.9) between VP, 1/QP, 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 n > 1), whereas P-wave properties provide simultaneous information about the confining (from VP, with n slightly less than 1) and pore pressure (from 1/QP, 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

Geophysical early warning of salt precipitation during geological carbon sequestration

Sequestration of industrial carbon dioxide (CO2) in deep geological saline aquifers is needed to mitigate global greenhouse gas emissions; monitoring the mechanical integrity of reservoir formations is essential for effective and safe operations. Clogging of fluid transport pathways in rocks from CO2-induced salt precipitation reduces injectivity and potentially compromises the reservoir storage integrity through pore fluid pressure build-up. Here, we show that early warning of salt precipitation can be achieved through geophysical remote sensing. From elastic P- and S-wave velocity and electrical resistivity monitoring during controlled laboratory CO2 injection experiments into brine-saturated quartz-sandstone of high porosity (29%) and permeability (1660 mD), and X-ray CT imaging of pore-scale salt precipitation, we were able to observe, for the first time, how CO2-induced salt precipitation leads to detectable geophysical signatures. We inferred salt-induced rock changes from (i) strain changes, (ii) a permanent ~ 1.5% decrease in wave velocities, linking the geophysical signatures to salt volume fraction through geophysical models, and (iii) increases of porosity (by ~ 6%) and permeability (~ 7%). Despite over 10% salt saturation, no clogging effects were observed, which suggests salt precipitation could extend to large sub-surface regions without loss of CO2 injectivity into high porosity and permeability saline sandstone aquifers