Fault normal stress (σ_n) changes dynamically during earthquakes. However, the impact of these changes on fault strength is poorly understood. We explore the effects of rapidly varying σ_n by conducting rotary-shear experiments on simulated fault gouges at 1 μm/s, under well-drained, hydrothermal conditions. Our results show both elastic and anelastic (time-dependent but recoverable) changes in gouge layer thickness in response to step changes and sinusoidal oscillations in σ_n. In particular, we observe dilation associated with marked weakening during ongoing σ_n-oscillations at frequencies >0.1 Hz. Moreover, recovery of shear stress after such oscillations is accompanied by transient (anelastic) compaction. We propose a microphysically based friction model that explains most of the observations made, including the effects of temperature and step versus sinusoidal perturbation modes. Our results highlight that σ_n-oscillations above a specific frequency threshold, controlled by the loading regime and frictional properties of the fault, may enhance seismic hazards.

@en

Laboratory studies suggest that seismogenic rupture on faults in carbonate terrains can be explained by a transition from high friction, at low sliding velocities (V), to low friction due to rapid dynamic weakening as seismic slip velocities are approached. However, consensus on the controlling physical processes is lacking. We previously proposed a microphysically based model (the “Chen–Niemeijer–Spiers” [CNS] model) that accounts for the (rate-and-state) frictional behavior of carbonate fault gouges seen at low velocities characteristic of rupture nucleation. In the present study, we extend the CNS model to high velocities (1 mm/s ≤ V ≤ 10 m/s) by introducing multiple grain-scale deformation mechanisms activated by frictional heating. As velocity and hence temperature increase, the model predicts a continuous transition in dominant deformation mechanisms, from frictional granular flow with partial accommodation by plasticity at low velocities and temperatures, to grain boundary sliding with increasing accommodation by solid-state diffusion at high velocities and temperatures. Assuming that slip occurs in a localized shear band, within which grain size decreases with increasing velocity, the model results capture the main mechanical trends seen in high-velocity friction experiments on room-dry calcite-rich rocks, including steady-state and transient aspects, with reasonable quantitative agreement and without the need to invoke thermal decomposition or fluid pressurization effects. The extended CNS model covers the full spectrum of slip velocities from earthquake nucleation to seismic slip rates. Since it is based on realistic fault structure, measurable microstructural state variables, and established deformation mechanisms, it may offer an improved basis for extrapolating lab-derived friction data to natural fault conditions.

@en

Rock materials show dramatic dynamic weakening in large-displacement (m), high-velocity (∼1 m/s) friction experiments, providing a mechanism for the generation of large, natural earthquakes. However, whether such weakening occurs during induced M3-4 earthquakes (dm displacements) is unknown. We performed rotary-shear experiments on simulated fault gouges prepared from the source-, reservoir- and caprock formations present in the seismogenic Groningen gas field (Netherlands). Water-saturated gouges were subjected to a slip pulse reaching a peak circumferential velocity of 1.2–1.7 m/s and total displacements of 13–20 cm, at 2.5–20 MPa normal stress. The results show 22%–81% dynamic weakening within 5–12 cm of slip, depending on normal stress and gouge composition. At 20 MPa normal stress, dynamic weakening from peak friction coefficients of 0.4–0.9 to 0.19–0.27 was observed, probably through thermal pressurization. We infer that similar effects play a key role during induced seismic slip on faults in the Groningen and other reservoir systems.

@en

A promising option for storing large-scale quantities of green gases (e.g., hydrogen) is in subsurface rock salt caverns. The mechanical performance of salt caverns utilized for long-term subsurface energy storage plays a significant role in long-term stability and serviceability. However, rock salt undergoes non-linear creep deformation due to long-term loading caused by subsurface storage. Salt caverns have complex geometries and the geological domain surrounding salt caverns has a vast amount of material heterogeneity. To safely store gases in caverns, a thorough analysis of the geological domain becomes crucial. To date, few studies have attempted to analyze the influence of geometrical and material heterogeneity on the state of stress in salt caverns subjected to long-term loading. In this work, we present a rigorous and systematic modeling study to quantify the impact of heterogeneity on the deformation of salt caverns and quantify the state of stress around the caverns. A 2D finite element simulator was developed to consistently account for the non-linear creep deformation and also to model tertiary creep. The computational scheme was benchmarked with the already existing experimental study. The impact of cyclic loading on the cavern was studied considering maximum and minimum pressure that depends on lithostatic pressure. The influence of geometric heterogeneity such as irregularly-shaped caverns and material heterogeneity, which involves different elastic and creep properties of the different materials in the geological domain, is rigorously studied and quantified. Moreover, multi-cavern simulations are conducted to investigate the influence of a cavern on the adjacent caverns. An elaborate sensitivity analysis of parameters involved with creep and damage constitutive laws is performed to understand the influence of creep and damage on deformation and stress evolution around the salt cavern configurations. The simulator developed in this work is publicly available at https://gitlab.tudelft.nl/ADMIRE_Public/Salt_Cavern.

@en

The strength properties of fault rocks at shearing rates spanning the transition from crystal-plastic flow to frictional slip play a central role in determining the distribution of crustal stress, strain, and seismicity in tectonically active regions. We review experimental and microphysical modelling work, which is aimed at elucidating the processes that control the transition from pervasive ductile flow of fault rock to rate-and-state-dependent frictional (RSF) slip and to runaway rupture, carried out at Utrecht University in the past 2 decades or so. We address shear experiments on simulated gouges composed of calcite, halite-phyllosilicate mixtures, and phyllosilicate-quartz mixtures performed under laboratory conditions spanning the brittle-ductile transition. With increasing shear rate (or decreasing temperature), the results consistently show transitions from (1) stable velocity-strengthening (v-strengthening) behaviour, to potentially unstable v-weakening behaviour, and (2) back to v strengthening. Sample microstructures show that the first transition seen at low shear rates and/or high temperatures represents a switch from pervasive, fully ductile deformation to frictional sliding involving dilatant granular flow in localized shear bands where intergranular slip is incompletely accommodated by creep of individual mineral grains. A recent microphysical model, which treats fault rock deformation as controlled by competition between ratesensitive (diffusional or crystal-plastic) deformation of individual grains and rate-insensitive sliding interactions between grains (granular flow), predicts both transitions well. Unlike classical RSF approaches, this model quantitatively reproduces a wide range of (transient) frictional behaviours using input parameters with direct physical meaning, with the latest progress focusing on incorporation of dynamic weakening processes characterizing co-seismic fault rupture. When implemented in numerical codes for crustal fault slip, the model offers a single unified framework for understanding slip patch nucleation and growth to critical (seismogenic) dimensions, as well as for simulating the entire seismic cycle.

@en

As part of a study to investigate methods to enhance pore/crack connectivity between the shale matrix and the induced fractures, we have investigated the matrix microstructure of an exposed analogue of the Jurassic Posidonia shales in the Dutch sub-surface. A combination of Precision-Ion- Polishing and Scanning Electron Microscopy has been used to image the in-situ porosity and mineralogy in shale samples from Whitby (UK), which are an analogue for the Dutch Posidonia shale. First results show a fine-grained mudstone with cm-sub-mm scale stratification. The section at Whitby can be divided into a clay matrix dominated upper half and a coarser grained, calcite-rich, lower half of the section. Commonly occurring minerals are pyrite, calcite (fossils, grains and cement), quartz, mica and dolomite. Organic matter content varies from 0 -2 % in the calcite dominated layers to 5 -18% in the clay matrix dominated layers. The most porous phases are the clay matrix and calcite fossils. This microstructural study shows which intervals within Posidonia shale contain the largest porosity and organic matter contents for sweet spot analyses and forms a basis for future research on enhancing connectivity between the pores and induced fractures.@en

A modular pore-scale model is developed to assess the response of wellbore cement to geological storage of CO₂. Numerical formulations for modeling of solute transport are presented and a methodology for coupling with geochemical processes is discussed, which includes: (1) advective and diffusive fluxes of solutes within the pore space, (2) aqueous phase speciation, (3) mineral dissolution-precipitation kinetics, and (4) the subsequent changes in pore space geometry. A Complex Pore Network Model (CPNM) is used to discretize the continuum porous structure as a network of pore bodies and pore throats, both with finite volumes. CPNM allows for a distribution of pore coordination numbers ranging between 1 and 26. This topological property, together with a geometrical distribution of pore sizes, enables the microstructure of porous media to be mimicked. For each pore element, transport of solute is calculated by solving the governing mass balance equations. Chemical reaction of the fluid phase with the main reactive solid components (portlandite and calcite) is incorporated through coupling with a geochemical reactive simulator. Average values and properties are obtained by integration over a large number of pores. Using this approach, we investigate how chemical reaction between water-bearing wellbore cement and supercritical CO₂ can create a distribution of porosity in a direction parallel to the CO₂ concentration gradient and transport path, at 50°C. The dynamics of this process involve interaction between diffusion dominated mass transport and the kinetics of dissolution and precipitation of portlandite and/or calcite. Simulation of unconfined chemical degradation, in a fluid of constant composition, shows development of different regions: (1) a zone adjacent to the inlet face, which is characterized by an increase in porosity due to extensive dissolution, (2) a carbonation zone with decreased porosity, (3) the carbonation front which made a thin layer with the lowest porosity due to calcium carbonate precipitation, and (4) dissolution zone. These results are in agreement with laboratory observations under similar conditions. This provides confidence that this pore-scale approach can ultimately be applied to model the progress of coupled CO₂ transport and cement degradation at critical points along the length of cemented wellbore sections at CO₂ storage sites.

@en