Geothermal energy plays a vital role in decarbonizing electricity and heat supply. Effective utilization of geothermal resources hinges on identifying or generating permeable reservoir zones and understanding how effective pressure variations affect fluid circulation and reservoir properties by poroelastic deformation. Hydrothermal alteration can modify the petrophysical properties of geothermal reservoir rocks, which may increase or decrease its productivity. Understanding these alteration effects is essential to predict and optimize long‐term sustainable geothermal operations. Here, we investigate the impact of hydrothermal alteration on poroelastic and hydraulic properties of diverse lithologies in a series of deformation tests performed at several confining (0–80 MPa) and pore pressure (10–30 MPa) levels. Experimental results of hydrothermally altered dikes and phyllites obtained from the Blue Mountain geothermal field (Nevada, USA) are compared to thermally cracked La Peyratte granite (France) and correlated with petrophysical properties, mineral composition, and microstructures. Argillic alteration of dikes increases porosity and storage capacity but lowers thermal conductivity and increases pore compressibility. Conversely, silicate precipitation in phyllites increases stiffness and thermal conductivity but also reduces porosity and permeability. Experimentally determined effective pressure coefficients range from 0.1 to 0.9, differ for permeability and volumetric strain and decrease with increasing effective pressure. The presence of compliant microcracks and crack‐like pores significantly increases the stress sensitivity of La Peyratte granite and silicified phyllites. This study demonstrates how thermal and chemical alteration impacts poromechanical and petrophysical characteristics of geothermal targets, which ultimately govern reservoir stability and subsidence, induced seismicity as well as fluid and heat extraction efficiency during geothermal operations.

Dynamic failure in the laboratory is commonly preceded by many foreshocks which accompany premonitory aseismic slip. Aseismic slip is also thought to govern earthquake nucleation in nature, yet, foreshocks are rare. Here, we examine how heterogeneity due to different roughness, damage and pore pressures affects premonitory slip and acoustic emission characteristics. High fluid pressures increase stiffness and reduce heterogeneity which promotes more rapid slip acceleration and shorter precursory periods, similar to the effect of low geometric heterogeneity on smooth faults. The associated acoustic emission activity in low-heterogeneity samples becomes increasingly dominated by earthquake-like double-couple focal mechanisms. The similarity of fluid pressure increase and roughness reduction suggests that increased stress and geometric homogeneity may substantially shorten the duration of foreshock activity. Gradual fault activation and extended foreshock activity is more likely observable on immature faults at shallow depth.

Unconsolidated, undrained triaxial deformation tests were performed on sandy facies Opalinus Clay at 50 MPa confining pressure to characterize the effect of water and microfabric orientation on the deformation behavior, mechanical properties, and P-wave velocity evolution. Dry and wet (≈ 8 and > 95% initial water saturation, respectively) samples with 12.6 ± 0.4 vol% porosity were deformed parallel and perpendicular to the bedding direction at a constant strain rate of 5 × 10 –6 s −1 . Dry samples revealed semi-brittle behavior and exhibited strain localization at failure, while deformation was more ductile at saturated conditions, promoting stable, slow faulting. Peak strength, Young’s modulus, and number of cumulative acoustic emissions decreased significantly for wet samples compared to dry samples; the opposite was observed for Poisson’s ratio. P-wave velocity anisotropy was significantly altered by differential stress, primarily due to the interplay between pore and fracture closure and stress-induced microcrack formation. For samples that were deformed perpendicular to bedding, we observed a reduction and reversal of P-wave velocity anisotropy with increasing differential stress, whereas anisotropy of parallel samples increased. The results suggest that water saturation reduces the pressure at the brittle-ductile transition and that the elastic properties and anisotropy of sandy facies Opalinus Clay can be significantly altered in an anisotropic stress field, e.g., adjacent to fault zones or tunnel excavations. Changes in elastic anisotropy are primarily controlled by the orientation between the pre-existing microfabric and the maximum principal stress direction, stress magnitude, and the degree of water saturation.

<p>Many reservoir rocks of productive geothermal energy resources display low porosity and matrix permeability. Therefore, to enhance fluid flow, fault zones and natural fracture networks are increasingly targeted for geothermal energy exploitation that are hydraulically connected to geothermal wells by stimulating the reservoir units. To this end, fluid is injected into the reservoir, which is generally believed to reduce effective stress and induce minor slip along stressed faults. Fluid injection can also lead to induced microseismicity and remotely-triggered earthquakes at great distances from the target reservoirs. The Blue Mountain geothermal field produces the largest seismic activity during maintenance shutdowns of injection and production requiring additional mechanisms such as poroelastic stress effects. In order to improve seismic hazard assessment and the understanding of induced seismicity around injection wells, we explore the coupling between matrix permeability, fault zone hydrology and mechanical behavior.</p><p>One main goal of this work is to provide insight into the scale-dependence of permeability by comparing laboratory results with previous permeability measurements using tidal responses in three different idle wells in Blue Mountain. We present a series of laboratory experiments performed on rock samples collected from the DB2 well at the Blue Mountain geothermal site in Humboldt County, Nevada, USA. The geothermal field benefits from the intersection of two W- and NW-directed normal faults resulting in high permeability of the geothermal reservoir production zone controlled by a brittle damage zone. Samples were obtained from two different lithologies, both of Triassic age, that constitute the reservoir. The first set of samples are low-porosity, (0.4 vol%) quartz-dominated (~50 – 60 wt%) phyllites, which exhibit zones with pronounced fracturing and elevated porosity (3.8 vol%). The second set of samples are felsic intrusive rocks with moderate to high porosity (7 – 15 vol%) due to strong hydrothermal alteration and clay mineral formation. Samples from both lithologies were selected from different sections of the damage zone showing varying degrees of faulting, from intact to highly brecciated, containing mineralized veins. We determine flow and poroelastic properties of cylindrical samples with a length of 2 cm and a diameter of 5 cm, subjected to stepwise cyclic variation of pore (<40 MPa) and confining pressure (<45 MPa). At each pressure step, we measure volumetric strain changes to derive the bulk modulus and effective stress coefficient, and use steady-state or pore pressure oscillation to determine permeability.</p><p>In addition to the tidal response and laboratory results, we developed a high-resolution seismicity catalog based on more than three years of continuous waveforms records from 2016 to 2019. We performed template-matching and differential travel-time inversions and use the resulting seismic events with magnitudes between 0.7 – 2.7 to search for seismicity migration patterns associated with discrete injection events. Integration of field and laboratory results can improve the characterization of the permeability structure of the fault zone at Blue Mountain and help to understand the mechanisms that trigger seismic events during production shutdown as well as the role of poroelastic stress release.</p>

The Opalinus Clay (OPA) formation is considered a suitable host rock candidate for nuclear waste storage. However, the sealing integrity and long-term safety of OPA are potentially compromised by pre-existing natural or artificially induced faults. Therefore, characterizing the mechanical behavior and microscale deformation mechanisms of faults and the surrounding rock is relevant for predicting repository damage evolution. In this study, we performed triaxial tests using saw-cut samples of the shaly and sandy facies of OPA to investigate the influence of pressure and mineral composition on the deformation behavior during fault reactivation. Dried samples were hydrostatically pre-compacted at 50 MPa and then deformed at constant strain rate, drained conditions and confining pressures ( p c ) of 5–35 MPa. Mechanical data from triaxial tests was complemented by local strain measurements to determine the relative contribution of bulk deformation and fault slip, as well as by acoustic emission (AE) monitoring, and elastic P-wave velocity measurements using ultrasonic transmissions. With increasing p c , we observe a transition from brittle deformation behavior with highly localized fault slip to semi-brittle behavior characterized by non-linear strain hardening with increasing delocalization of deformation. We find that brittle localization behavior is limited by p c at which fault strength exceeds matrix yield strength. AEs were only detected in tests performed on sandy facies samples, and activity decreased with increasing p c . Microstructural analysis of deformed samples revealed a positive correlation between increasing p c and gouge layer thickness. This goes along with a change from brittle fragmentation and frictional sliding to the development of shear zones with a higher contribution of cataclastic and granular flow. Friction coefficient at fault reactivation is only slightly higher for the sandy ( µ ~ 0.48) compared to the shaly facies ( µ ~ 0.4). Slide-hold-slide tests performed after ~ 6 mm axial shortening suggest stable creeping and long-term weakness of faults at the applied conditions. Our results demonstrate that the mode of fault reactivation highly depends on the present stress field and burial history.

We experimentally determined the hydraulic properties of fractures within various rock types, focusing on a variety of Variscan rocks. Flow-through experiments were performed on slate, graywacke, quartzite, granite, natural fault gouge, and claystone samples containing an artificial fracture with a given roughness. For slate samples, the hydraulic transmissivity of the fractures was measured at confining pressures, pc, at up to 50 MPa, temperatures, T, between 25 and 100 °C, and differential stress, σ, acting perpendicular to the fracture surface of up to 45 MPa. Fracture transmissivity decreases non-linearly and irreversibly by about an order of magnitude with increasing confining pressure and differential stress, with a slightly stronger influence of pc than of σ. Increasing temperature reduces fracture transmissivity only at high confining pressures when the fracture aperture is already low. An increase in the fracture surface roughness by about three times yields an initial fracture transmissivity of almost one order of magnitude higher. Fractures with similar surface roughness display the highest initial transmissivity within slate, graywacke, quartzite and granite samples, whereas the transmissivity in claystone and granitic gouge material is up to several orders of magnitude lower. The reduction in transmissivity with increasing stress at room temperature varies with composition and uniaxial strength, where the deduction is lowest for rocks with a high fraction of strong minerals and associated high brittleness and strength. Microstructural investigations suggest that the reduction is induced by the compaction of the matrix and crushing of strong asperities. Our results suggest that for a given surface roughness, the fracture transmissivity of slate as an example of a target reservoir for unconventional EGS, is comparable to that of other hard rocks, e.g., granite, whereas highly altered and/or clay-bearing rocks display poor potential for extracting geothermal energy from discrete fractures.

Studying the mechanical properties of argillaceous rocks is of major interest in geoscience. For example, these rocks are important in engineering applications such as being suitable cap-rocks for the geological storage of carbon dioxide and potential host rocks for the storage of nuclear waste. Furthermore, argillaceous rocks are encountered in different natural settings such as accretionary wedges or fault zones. As a result of their sedimentary and diagenetic history clay rich rocks are often characterised by multiscale textural anisotropy and compositional heterogeneity resulting in anisotropic mechanical and hydraulic properties. Here, we studied the anisotropic deformation behaviour of Opalinus Clay, collected from the Mont Terri underground laboratory, which is the envisaged host rock formation for nuclear waste disposal in Switzerland. We used the sandy facies of Opalinus Clay, characterized by an irregular wavy lamination of quartz-rich and carbonate-cemented lenses with clay-rich interlayers. Unconsolidated-cylindrical samples cored at 0°, 45° and 90° to the macroscopically visible bedding were deformed in undrained constant strain rate experiments using a Paterson-type deformation apparatus. For each orientation, tests were performed at dry conditions varying either confining pressure (in the range of 50 - 100 MPa), temperature (25 - 200 °C) or strain rate (1*10-3 - 5*10-6 s-1) to study the influence of testing condition and sample orientation on the deformation behaviour. In addition, we deformed a set of back saturated samples at fixed conditions of 50 MPa, 100 °C and 5*10-4 s-1 to investigate the effect of water content on the material strength. The results show semi-brittle deformation with low yield strength and strain weakening behaviour, in which strain is localized in sub-millimetre to millimetre-wide shear zones at all conditions. Increasing water content reduces, whereas increasing confining pressure increases the peak strength. Samples that were deformed parallel to bedding orientation exhibit the highest strength compared to samples with an orientation of 90° and 45° to bedding. Only for the latter orientations a weak correlation was found between temperature and failure behaviour. The variation of strain rate shows no clear influence for all orientations. Within this test series, there appears to be a potentially greater influence of the porosity on the peak strength for 45° and 90° oriented samples. Clay rich layers seem to have a strong influence on localization and formation of shear zones, in particular for samples oriented at 45° and 90° to bedding. This observation was confirmed by electron microscopy performed on broad ion beam polished surfaces of deformed sample material. Our experiments reveal that water content, sample orientation with respect to bedding and confining pressure are the most important factors influencing the peak strength of the sandy facies of Opalinus Clay, whereas compositional heterogeneity is responsible for the localization behaviour.

Extensive groundwater extraction leads to the drawdown of the groundwater table, which can increase the stress of soils and lead to land subsidence. Analyses of land subsidence generally require numerical models based on the poroelasticity theory. FEFLOW is a well-known commercial finite element software that is one of the most widely used for groundwater, heat transfer, and transport modelling. However, it is not capable of modelling land subsidence. Here, we introduce a plug-in called SUB+ and a pre-processing software PSUB that can be used for fully coupled land subsidence analyses using the poroelasticity theory. The plug-in and the software are developed using C++ with Qt-Framework. The Eigen and Intel MKL libraries are used for matrix operations and the sparse matrix solver. Although FEFLOW is closed-source, SUB+ and PSUB are open-source and available on GitHub.