Induced earthquakes are still highly unpredictable, and often caused by variations in pore fluid pressure. Monitoring and understanding the mechanisms of fluid-induced fault slip is essential for seismic risk mitigation and seismicity forecasting. Fluid-induced slip experiments were performed on critically stressed faulted sandstone samples, and the evolution of the actively sent ultrasonic waves throughout the experiment was measured. Two different fault types were used: smooth saw-cut fault samples at a 35º angle, and a rough fault created by in-situ faulting of the samples. Variations in the seismic slip velocity and friction along the fault plane were identified by the coda of the ultrasonic waves. Additionally, ultrasonic amplitudes show precursory signals to laboratory fault reactivation. Our results show that small and local variations in stress before fault failure can be inferred using coda wave interferometry for time-lapse monitoring, as coda waves are more sensitive to small perturbations in a medium than direct waves. Hence, these signals can be used as precursors to laboratory fault slip and to give insight into reactivation mechanisms. Our results show that time-lapse monitoring of coda waves can be used to monitor local stress changes associated with fault reactivation in this laboratory setting of fluid-induced fault reactivation. This is a critical first step towards a method for continuous monitoring of natural fault zones, contributing to seismic risk mitigation of induced and natural earthquakes.

Fracture networks are abundant in subsurface applications (e.g., geothermal energy production, CO2 sequestration). Fractured reservoirs often have a very complex structure, making modeling flow and transport in such networks slow and unstable. Consequently, this limits our ability to perform uncertainty quantification and increases development costs and environmental risks. This study provides an advanced methodology for simulation based on Discrete Fracture Model (DFM) approach. The preprocessing framework results in a fully conformal, uniformly distributed grid for realistic 2D fracture networks at a required level of precision. The simplified geometry and topology of the resulting network are compared with input (i.e., unchanged) data to evaluate the preprocessing influence. The resulting mesh-related parameters, such as volume distributions and orthogonality of control volume connections, are analyzed. Furthermore, changes in fluid-flow response related to preprocessing are evaluated using a high-enthalpy two-phase flow geothermal simulator. The simplified topology directly improves meshing results and, consequently, the accuracy and efficiency of numerical simulation. The main novelty of this work is the introduction of an automatic preprocessing framework allowing us to simplify the fracture network down to required level of complexity and addition of a fracture aperture correction capable of handling heterogeneous aperture distributions, low connectivity fracture networks, and sealing fractures. The graph-based framework is fully open-source and explicitly resolves small-angle intersections within the fracture network. A rigorous analysis of changes in the static and dynamic impact of the preprocessing algorithm demonstrates that explicit fracture representation can be computationally efficient, enabling their use in large-scale uncertainty quantification studies.

Using active ultrasonic source survey data, Coda-wave Decorrelation (CWD) time-lapse imaging during the triaxial compression of Whitby Mudstone cores provides a 3-D description of the evolution and redistribution of inelastic strain concentrations. Acoustic Emissions (AEs) monitoring is also performed between any two consecutive surveys. From these data, we investigate the impact of initial water saturation $S_w$ on the onset, growth, and reactivation of inelastic deformation, compared to the post-deformation fracture network extracted from X-ray tomography scans. Our results indicate for the applied strain-rate and degree of initial water saturation, and within the frequency range of our ultrasonic transducers (0.1 to 1 MHz), that inelastic strain localisation and propagation in the Whitby Mudstone does not radiate AEs of sufficient magnitude to be detected above the average noise level. This is true for both the initial onset of inelasticity (strain localisation), and during macroscopic failure. In contrast, the CWD results indicate the onset of what is interpreted as localised regions of inelastic strain at less than fifty percent of the peak differential stress the Whitby Mudstone can sustain. The seemingly aseismic nature of these clay-rich rocks suggests the gradual development of inelastic strain, from the microscopic diffuse damage, up until the macroscopic shear failure.

Pore-pressure depletion in sandstone reservoirs is well known to cause both elastic and inelastic compaction, often resulting in notable surface subsidence and induced seismicity. Recent studies indicate that in such cases inelastic strain, which is often neglected in geomechanical models, represents a significant proportion of the total strain throughout reservoir production. While there has been considerable effort to quantify the proportion of continuous inelastic deformation from the mechanical response of laboratory samples, there has been little focus to date on the associated acoustic response throughout compaction. With this in mind, we employ three coda-wave based processing methods for the active source monitoring of ultrasonic velocity, scattering power, and intrinsic/scattering attenuation during the pore-pressure depletion of core samples from the Slochteren sandstone reservoir in the Groningen gas field (the Netherlands). Our results corroborate previous studies suggesting that initially, inelastic deformation occurs primarily along intergranular boundaries, with intergranular cracking developing towards the end of depletion and particularly for the highest porosity samples. Furthermore, analysis of Biot type intrinsic attenuation indicates that this compaction occurs in several stages of predominately intergranular closure, transitioning into predominantly intergranular slip/cracking, and eventually porosity-dependent intragranular cracking. We demonstrate how this segmentation of pore-pressure driven compaction can be used to characterise differences in sample properties, and monitor the evolution of microstructural inelastic deformation throughout depletion. We further discuss the feasibility of in/cross-borehole monitoring of reservoir compaction, for both improved geo-mechanical modelling and early warning detection of induced seismicity.

Predicting stress changes in the subsurface leading to failure or seismicity remains challenging. Developing a robust monitoring method can help the prediction and thus mitigation of natural hazards. Ultrasonic transmission experiments were performed on Red Pfaelzer sandstones to investigate the forecasting potential to failure at different confining pressures. The forecasting potential for failure of the energy of the direct and coda wave, the transmissivity, Q-factor, coda wave decorrelation coefficient, and velocity change by coda wave interferometry are investigated and compared. Our results show the failure of the tested samples can be forecasted from 40 to 70% of the failure point. Small differences are visible in the precursors between the tested confining pressures, but as the trends are very similar, a robust prediction of failure can be made by combining the various analyses techniques. In this paper, we propose a traffic light forecasting system using active acoustic monitoring which is applicable for forecasting failure at various depths and or stress conditions, for a better prediction of small stress-induced changes in the subsurface and thus mitigation of failure (and seismicity) in the subsurface.

Our understanding of the past behavior of the geomagnetic field arises from magnetic signals stored in geological materials, e.g. (volcanic) rocks. Bulk rock samples, however, often contain magnetic grains that differ in chemistry, size and shape; some of them record the Earth’s magnetic field well, others are unreliable. The presence of a small amount of adverse behaved magnetic grains in a sample may already obscure important information on the past state of the geomagnetic field. Recently it was shown that it is possible to determine magnetizations of individual grains in a sample by combining X-ray computed tomography and magnetic surface scanning measurements. Here we establish this new Micromagnetic Tomography (MMT) technique and make it suitable for use with different magnetic scanning techniques, and for both synthetic and natural samples. We acquired reliable magnetic directions by selecting subsets of grains in a synthetic sample, and we obtained rock-magnetic information of individual grains in a volcanic sample. This illustrates that MMT opens up entirely new venues of paleomagnetic and rock-magnetic research. MMT’s unique ability to determine the magnetization of individual grains in a nondestructive way allows for a systematic analysis of how geological materials record and retain information on the past state of the Earth’s magnetic field. Moreover, by interpreting only the contributions of known magnetically well-behaved grains in a sample MMT has the potential to unlock paleomagnetic information from even the most complex, crucial, or valuable recorders that current methods are unable to recover.