Fluid extraction from sandstone reservoirs leads to reservoir compaction, potentially inducing surface subsidence and seismicity, as observed in the Groningen Gas Field, Netherlands. Such compaction is partly elastic, but can additionally be caused by instantaneous plastic and rate/time-dependent processes, such as subcritical crack growth, meaning that compaction may continue even if production is stopped. Despite the need to evaluate the impact of post-abandonment reservoir behavior (>10–100 years), few mechanism-based, rate/time-dependent compaction laws exist. Compaction due to grain breakage, either via critical or subcritical crack growth, is driven by tensile stresses acting on surface and volume flaws. We performed high-resolution 3D linear elastic finite element method simulations on simplified grain assemblies to investigate the effect of stress–strain boundary conditions, porosity and mineralogical variations on grain-scale stress fields. Our simulations showed tensile stress concentrations at grain contact edges and on pore walls, which increased in magnitude with increasing aggregate porosity and local porosity variation. The fraction of surface area with tensile stresses sufficient to extend flaws with a size up to 30μm showed a clear correlation with compactive yield envelopes for the Groningen reservoir sandstone. This suggests that compactive failure is related to the probability of pre-existing surface flaws, falling in a pore surface region where the Griffith criterion is satisfied. A preliminary, time-independent failure probability model, using the observed tensile stress distribution, qualitatively predicts a non-linear increase in grain cracking during deviatoric loading, and suggests a new route to predict sandstone compaction through brittle grain failure.

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.

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.

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.

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.

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.