Introduction to this special section: The role of geophysics in a net-zero-carbon world

Human activities are changing the earth's climate, causing increasingly disruptive social and ecological impacts. These impacts can be reduced if global carbon dioxide (CO2) emissions reach net zero in the near future. A net-zero-carbon world can be achieved by using energy more efficiently and responsibly; transitioning toward energy sources, products, and services that minimize greenhouse gas release; and implementing existing and novel technologies to remove and store CO2 from the atmosphere

Experimental study of geophysical and transport properties of salt rocks in the context of underground energy storage

Artificial caverns in salt rock formations play an important role in the net-zero energy transition challenge, both for covering short-term fluctuations in energy demand and serving as safe locations for long-term underground gas storage both for hydrogen and natural gas. Geophysical tools can serve for monitoring geomechanical changes in the salt cavern during selection and development, and during gas storage/extraction activities, but the use of common geophysical monitoring techniques has been very limited in this area. Here, we present experimental work on physical and transport properties of halite rocks within the energy storage context and assess the potential of seismic and electromagnetic data to monitor gas storage activities in salt formations. First, we analysed the stress-dependency of the elastic and transport properties of five halite rocks to improve our understanding on changes in the geological system during gas storage operations. Second, we conducted two dissolution tests, using cracked and intact halite samples, monitored with seismic (ultrasonic P- and S-waves velocities and their attenuation factors) and electromagnetic (electrical resistivity) sources to evaluate (i) the use of these common geophysical sensing methods to remotely interpret caverning development and (ii) the effect of structural discontinuities on rock salt dissolution. Elastic properties and permeability showed an increasing trend towards rock sealing and mechanical enhancement with increasing pressure for permeabilities above 10−21 m2, with strong linear correlations up to 20 MPa. In the dissolution tests, the ultrasonic waves and electrical resistivity showed that the presence of small structural discontinuities largely impacts the dissolution patterns. Our results indicate that seismic and electromagnetic methods might help in the selection and monitoring of the caverning process and gas storage operations, contributing to the expected increase in demand of large-scale underground hydrogen storage

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

Alteration of ultrasonic signatures by stress-induced changes in hydro-mechanical properties of fractured rocks

In this study, we evaluate the changes in ultrasonic signatures (i.e., frequency content, velocity, amplitude, and attenuation) due to stress-induced alteration of fracture aperture/permeability. Flow-through experiments were performed on artificially-fractured phyllite specimens along with the concurrent measurements of ultrasonic signatures under different stress conditions. Increasing pore pressure led to fracture opening, as indicated by increases in both mechanical and hydraulic apertures. In addition, we observed that increase in confining pressure (and decrease in pore pressure) led to increases in ultrasonic velocities, ultrasonic amplitudes, and fracture specific stiffness, and decrease in ultrasonic attenuations. It was found that time-frequency partitioning depends on hydraulic aperture. The higher frequency band, for both P- and S-waves, was insensitive to the changes in stress conditions; the lower band was sensitive to the changes in stress conditions, as long as the hydraulic aperture was changing. Three-Element rheological and Power-Law models successfully predicted the time-dependent fracture displacement, with the former being more accurate at higher levels of pore pressures

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

Introduction to special section on the rock physics contribution to the energy transition challenge

Rock physics connects with geophysics, petrophysics and geomechanics to adequately characterize geological reservoirs, optimize monitoring operations in the field, interpret in situ and laboratory test data, and develop accurate predictive models for extraction/injection activities. The application of rock physics is crucial to achieving net-zero carbon emissions worldwide, as we need to combine large-scale mitigation technologies like carbon capture usage and storage, together with an increasing use of renewables such as geothermal and underground hydrogen storage (UHS)

CO2‐Brine substitution effects on ultrasonic wave propagation through sandstone with oblique fractures

Seismic monitoring of injected CO2 plumes in fractured storage reservoirs relies on accurate knowledge of the physical mechanisms governing elastic wave propagation, as described by appropriate, validated rock physics models. We measured laboratory ultrasonic velocity and attenuation of P and S waves, and electrical resistivity, of a synthetic fractured sandstone with obliquely aligned (penny‐shaped) fractures, undergoing a brine‐CO2 flow‐through test at simulated reservoir pressure and temperature. Our results show systematic differences in the dependence of velocity and attenuation on fluid saturation between imbibition and drainage episodes, which we attribute to uniform and patchy fluid distributions, respectively, and the relative permeability of CO2 and brine in the rock. This behavior is consistent with predictions from a multifluid rock physics model, facilitating the identification of the dispersive mechanisms associated with wave‐induced fluid flow in fractured systems at seismic scales

Shale distribution effects on the joint elastic–electrical properties in reservoir sandstone

We investigated the effect of shale distribution on the joint elastic wave and electrical properties of shaly reservoir sandstones using a dataset of laboratory measurements on 75 brine-saturated (35 g/L salinity) rock samples (63 samples from the literature, 12 newly measured samples). All the data were collected using the ultrasonic (700 kHz) pulse-echo measurement technique for P- and S-wave velocities (Vp, Vs), attenuations (Qp−1, Qs−1), and a four-electrode method for resistivity under elevated hydrostatic confining pressures between 10 and 50 MPa (pore fluid pressure 5 MPa). The distribution of volumetric shale content was classified by comparing the calculated dry P-wave modulus to the modified Upper Hashin–Shtrikman bound for quartz and air mixtures, assuming pore-filling shale. This scheme in particular allowed us to distinguish between pore-filling and load-bearing shale distributions according to idealized definitions, which provides new insight into the joint ultrasonic properties and resistivity behaviour for shaly sandstones. In resistivity–velocity space, the resistivity of load-bearing shale increases with increasing velocity which form a more distinct trend with steeper gradient compared to those for partial pore-filling shale and clean sandstones. Moreover, the pore-filling shale trend straddles the clean sandstone trend and meets the load-bearing shale trend between 100 and 150 apparent formation factors. In resistivity–attenuation space, the highest attenuations exist when the volumetric shale content is close to the frame porosity (for Qp−1 in particular), at the transition between pore-filling and load-bearing shales. The results will inform the development of improved rock physics models to aid reservoir characterization from geophysical remote sensing, particularly for joint seismic and controlled source electromagnetic surveys

Transport properties of saline CO2 storage reservoirs with unconnected fractures from brine-CO2 flow-through tests

CO2 storage in fractured reservoirs may lead to fast CO2 flow through interconnected fracture networks; but the role of isolated fractures on brine-CO2 multiphase flow systems remains unclear. We present the results of a brine-CO2 flow-through experiment in which we assess the change in transport properties of a synthetic sandstone plug (a surrogate of a saline siliciclastic CO2 reservoir) containing non-connected fractures aligned 45° from its axis. The test was performed at 40 MPa of constant hydrostatic confining pressure and ~11 MPa of pore pressure, at room temperature (~19.5 °C), using pure liquid-CO2 and 35 g L−1 NaCl salt solution. The injected CO2-brine volume fraction was increased from 0 to 1 in 0.2 units-steps (drainage). Upon achievement of the maximum CO2 saturation (SCO2 ~0.6), the plug was flushed-back with the original brine (imbibition). During the test, we monitored simultaneously pore pressure, temperature, axial and radial strains, and bulk electrical resistivity. The fractured sample showed lower values of cross- and end-points in the relative permeability curves to CO2 compared to non-fractured samples, from comparable experiments performed at similar pressure and brine salinity conditions, but different temperature. Our results suggest that a non-connected fracture network affects the mobility of the individual phases, favouring the trapping of CO2 in the porous medium and improving the storage efficiency of the reservoir. These evidences show the need of a better understanding of fracture connectivity prior to discard fractured reservoirs as unsuitable geological formations for CO2 storage

Uplift and exposure of serpentinized massifs: Modeling differential serpentinite diapirism and exhumation of the Troodos Mantle Sequence, Cyprus

Serpentinized mantle peridotites form prominent mountains, including the highest elevations of the Troodos ophiolite in Cyprus (Mount Olympus, 1,952 m), but to date, only qualitative mechanisms have been proposed to explain the uplift of mantle rocks to high altitudes. Serpentinization reactions between mantle rocks and water result in profound changes to the rheology and physical properties of peridotites including significant density reduction (∼900 kg/m3). Field observations, density measurements, and isostatic uplift and erosional modeling provide new constraints on the contribution of serpentinization to the uplift of the Troodos Mantle Sequence. Different serpentinization styles have resulted in two distinct serpentinite domains with contrasting densities. Our modeling shows that the Troodos Mountains can form within the geologically constrained uplift time frame (∼5.5 Myr) exclusively through partial serpentinization reactions. We interpret the serpentinite domains as two nested diapirs that formed due to different extents of serpentinization and density reduction. Differential uplift and exhumation have decoupled the two serpentinite diapirs from the originally overlying ocean crustal rocks. Once at high altitudes the incursion of meteoric water reinforced coupled deformation-alteration-recrystallization processes in the shallow subsurface producing a localized low density completely serpentinized diapir. A second decoupling between the contrasting serpentinite diapirs results in localized differential uplift and exhumation, extruding deep materials to the east of Mount Olympus. Application of our modeling to other serpentinite massifs (e.g., St. Peter and St. Paul Rocks, New Idria, California) highlights the contribution of isostasy to the uplift of serpentinized massifs