The potential response of the hydrate reservoir in the South Shetland Margin, Antarctic Peninsula, to ocean warming over the 21st century

In the South Shetland Margin (SSM), Antarctic Peninsula, a bottom-simulating reflector indicates the presence of hydrate between ca. 500 and 3000 m water depth (mwd). The cold seabed temperatures allow hydrate stability at shallower water depths. During the past five decades, the Antarctic Peninsula has been warming up faster than any other part of the Southern Hemisphere, and long-term ocean warming could affect the stability of the SSM hydrate reservoir at shallow waters. Here, we model the transient response of the SSM hydrate reservoir between 375 and 450 mwd to ocean warming for the period 1958–2100. For the period 1958–2010, seabed temperatures are given by oceanographic measurements in the area, and for 2010–2100 by two temperature scenarios represented by the observed trends for the periods 1960–2010 (0.0034°C y−1) and 1980–2010 (0.023°C y−1). Our results show no hydrate-sourced methane emissions for an ocean warming rate at the seabed of 0.0034 °C y−1. For a rate of 0.023°C y−1, emissions start in 2028 at 375 mwd and extend to 442 mwd at an average rate of about 0.91 mwd y−1, releasing ca. 1.13×103 mol y−1 of methane per metre along the margin by 2100. These emissions originate from dissociation at the top of the hydrate layer, a physical process that steady-state modelling cannot represent. Our results are speculative on account of the lack of direct evidence of a shallow water hydrate reservoir, but they illustrate that the SSM is a key area to observe the effects of ocean warming-induced hydrate dissociation in the coming decades

A quick-look method for initial evaluation of gas hydrate stability below subaqueous permafrost

Many studies demonstrated the coexistence of subaqueous permafrost and gas hydrate. Subaqueous permafrost could be a factor affecting the formation/dissociation of gas hydrate. Here, we propose a simple empirical approach that allows estimating the steady-state conditions for gas hydrate stability in the presence of subaqueous permafrost. This approach was derived for pressure, temperature, and salinity conditions typical of subaqueous permafrost in marine (brine) and lacustrine (freshwater) environments

Spatial and temporal evolution of rifting and continental breakup in the Eastern Black Sea Basin revealed by long‐offset seismic reflection data

The age and distribution of the synrift and early postrift infill records the spatial and temporal distribution of extension and breakup processes in a rift basin. The Eastern Black Sea Basin (EBSB) is thought to have formed by back‐arc extension during Cretaceous to Early Cenozoic time. However, a lack of direct constraints on its deep stratigraphy leaves uncertainties over the time, duration, and location for rifting and breakup processes in the basin. Here we use the enhanced imaging provided by 2‐D long‐offset seismic reflection profiles to analyze the deep structural and stratigraphic elements of the EBSB. Based on these elements, we infer the presence of two distinct Late Cretaceous synrift units, recording initial extension (rift stage 1) over the continental highs (Shatsky Ridge and the Mid Black Sea High), followed by strain localization along the major basin‐bounding faults and rift migration toward the basin axis (rift stage 2). Overlying these units, Palaeocene(?)‐Eocene and Oligocene units show a synkinematic character in the NW, with evidence for ongoing extension until Oligocene time. Toward the SE, these sequences are instead postkinematic, directly overlaying a basement emplaced during breakup. We interpret the Palaeocene(?)‐Oligocene units to record the time spanning from the initiation of breakup (Late Cretaceous‐Palaeocene, in the SE) to the end of extension (Oligocene, in the NW). The first ubiquitously postrift infill is the Lower Miocene Maykop Formation. Our results highlight the along‐strike temporal variability of extension and breakup processes in the EBSB

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

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)

Laboratory observations of frequency-dependent ultrasonic P-wave velocity and attenuation during methane hydrate formation in Berea sandstone

Knowledge of the effect of methane hydrate saturation and morphology on elastic wave attenuation could help reduce ambiguity in seafloor hydrate content estimates. These are needed for seafloor resource and geohazard assessment, as well as to improve predictions of greenhouse gas fluxes into the water column. At low hydrate saturations, measuring attenuation can be particularly useful as the seismic velocity of hydrate-bearing sediments is relatively insensitive to hydrate content. Here, we present laboratory ultrasonic (448–782 kHz) measurements of P-wave velocity and attenuation for successive cycles of methane hydrate formation (maximum hydrate saturation of 26 per cent) in Berea sandstone. We observed systematic and repeatable changes in the velocity and attenuation frequency spectra with hydrate saturation. Attenuation generally increases with hydrate saturation, and with measurement frequency at hydrate saturations below 6 per cent. For hydrate saturations greater than 6 per cent, attenuation decreases with frequency. The results support earlier experimental observations of frequency-dependent attenuation peaks at specific hydrate saturations. We used an effective medium rock-physics model which considers attenuation from gas bubble resonance, inertial fluid flow and squirt flow from both fluid inclusions in hydrate and different aspect ratio pores created during hydrate formation. Using this model, we linked the measured attenuation spectral changes to a decrease in coexisting methane gas bubble radius, and creation of different aspect ratio pores during hydrate formation

Reactive transport modelling insights into CO2 migration through sub-vertical fluid flow structures

Sub-vertical geological structures that cut through the overburden, usually called chimneys or pipes, are common in sedimentary basins. Chimneys behave as conduits that hydraulically connect deep strata with the overburden and seabed. Hence, if stored CO2 migrates to a sufficiently high permeability chimney the risk of CO2 leakage at the seabed increases. Despite the possible negative effects these structures may have on the integrity of CO2 storage sites, little is known about (i) their effective permeability distribution, controlled by the combined role of fractures and matrix, and (ii) feedback mechanisms between porosity-permeability, CO2 reactivity and mineralogy within them. Reactive transport modelling is used to perform 2D axisymmetric radial simulations of geological systems containing chimneys. CO2 saturations of 10%, 30% and 50% are imposed on a cell located next to the symmetry axis at the base of the model. Under hydrostatic conditions, CO2 reaches the seabed, at 500 m above the injection point, in less than 100 yr using injected CO2 saturations at or above 30% and with overburden isotropic permeabilities and chimney vertical permeabilities above 10−14  m2. Vertical fractures with apertures larger than 0.05 mm for volume fractions below 1% are sufficient to sustain such high vertical permeabilities in the chimney with a relatively high cap rock matrix permeability of 10−16 m2. Over 100 yr of CO2 injection, changes in porosity and permeability due to mineral precipitation/dissolution are negligible. For this time scale, in systems containing chimneys sufficiently far away from the injection well, the risk of CO2 leakage at the seabed is primarily controlled by the pre-existing hydrogeological state of the system

Core-scale geophysical and hydromechanical analysis of seabed sediments affected by CO2 venting

Safe offshore Carbon Capture Utilization and Storage (CCUS) includes monitoring of the subseafloor, to identify and assess potential CO2 leaks from the geological reservoir through seal bypass structures. We simulated CO2-leaking through shallow marine sediments of the North Sea, using two gravity core samples from ∼1 and ∼2.1 m below seafloor. Both samples were subjected to brine−CO2 flow-through, with continuous monitoring of their transport, elastic and mechanical properties, using electrical resistivity, permeability, P-wave velocity and attenuation, and axial strains. We used the collected geophysical data to calibrate a resistivity-saturation model based on Archie’s law extended for clay content, and a rock physics for the elastic properties. The P-wave attributes detected the presence of CO2 in the sediment, but failed in providing accurate estimates of the CO2 saturation. Our results estimate porosities of 0.44 and 0.54, a background permeability of ∼10−15 and ∼10-17 m2, and maximum CO2 saturation of 18 % and 10 % (±5 %), for the sandier (shallower) and muddier (deeper) sample, respectively. The finer-grained sample likely suffered some degree of gas-induced fracturing, exhibiting an effective CO2 permeability increase sharper than the coarser-grained sample. Our core-scale multidisciplinary experiment contributes to improve the general interpretation of shallow sub-seafloor gas distribution and migration patterns

Experimental assessment of pore fluid distribution and geomechanical changes in saline sandstone reservoirs during and after CO2 injection

Responsible CO2 geosequestration requires a comprehensive assessment of the geomechanical integrity of saline reservoir formations during and after CO2 injection. We assessed the geomechanical effects of CO2 injection and post-injection aquifer recharge on weakly cemented, synthetic-sandstone (38% porosity) sample in the laboratory under dry and brine-saturated conditions, before and after subjecting the sample to variable pore pressure brine-CO2 flow-through tests (∼170 h). We measured ultrasonic P- and S-wave velocities (Vp, Vs) and attenuations, electrical resistivity and volumetric strain (εv). Vs was found to be an excellent indicator of mechanical deformation during CO2 injection; Vp gives mechanical and pore fluid distribution information, allowing quantification of the individual contribution of both phenomena when combined with resistivity. Abrupt strain recovery during imbibition suggests that aquifer recharge after ceasing CO2 injection might affect the geomechanical stability of the reservoir. Static and dynamic parameters indicate the sample experienced minor geomechanical changes during CO2 exposure, with an increase of Δεv <3% and a drop in ΔVs ∼1%. In contrast, due to brine-induced hydro-mechanical alteration, Δεv increased by ∼10% and ΔVs by ∼6%. This study provides a multiparameter, thermo-hydro-mechanical-chemical database needed to validate monitoring tools and simulators, for prediction of the geomechanical behaviour of CO2 storage reservoirs