Rocks can undergo fatigue failure when subjected to cyclic mechanical, hydraulic, or thermal loadings, or a combination of these. Therefore, accounting for possible fatigue damage is important for subsurface engineering projects, such as the cyclic stimulation of geothermal reservoirs. However, existing models do not simultaneously account for degradation of both tensile strength and stiffness under varying-amplitude loading and coupled thermo-hydro-mechanical (THM) conditions. To address this, a new cohesive zone model is developed to account for the effect of fatigue on tensile strength and stiffness. The model is then used within the framework of zero-thickness interface elements to simulate the response of pre-existing or new fractures. Hydraulic and thermal processes are included in both the cohesive interface elements and the continuum elements, allowing the consideration of coupled thermo-hydro-mechanical processes. The fatigue damage variable is set to evolve with the number and magnitude of cycles according to Palmgren-Miner's rule. The proposed method is validated against three laboratory tests from the literature, including cyclic Brazilian test, cyclic hydraulic fracturing test and cyclic thermal stimulation test. All three validation results show that the fatigue damage or reduced breakdown pressure can be well reproduced. Mesh sensitivity based on the simulation of the Brazilian test, in which interface elements are inserted in-between all the continuum elements, highlights the influence of the mesh orientation and mesh density on the simulation results. In addition, stabilisation of the method is demonstrated by increasing the mechanical viscosity, which must be used with care to avoid predicting a longer fatigue life. The ability of the method to handle varying-amplitude cyclic loading is demonstrated by the simulation of a synthetic cyclic loading scheme based on the Brazilian test. The proposed method can be used to support the design of cyclic thermal stimulation campaigns for geothermal (or other) reservoirs, by being able to simulate the reduction in strength due to fatigue, and thus reducing stimulation pressures needed.

Soft stimulation technologies have been proposed as a means to reduce the breakdown pressure and mitigate the risk of induced seismicity during geothermal reservoir stimulation. Yet, the underlying mechanisms remain poorly understood due to the complexity of the coupled thermo-hydro-mechanical (THM) processes. In this work, a fully coupled THM model is developed to evaluate and compare the performance of different stimulation scenarios (monotonic, stepwise injection rate, cyclic injection rate or temperature, and stepwise combined with cyclic injection rate stimulation) on a synthetic, highly permeable reservoir with near-borehole clogging. Simulation results show that stepwise injection rate stimulation yields the most favourable outcomes, followed by the stepwise injection rate combined with cyclic injection rate stimulation. On the other hand, fatigue effects are seen to play a negligible role in the improved performance since the tensile stress at the fracture tip is relaxed with the continuous fracture growth. In addition, cyclic injection temperature stimulation is generally neither better nor worse than monotonic stimulation, but has slightly different characteristics, creating more local damage controlled by the period of the injection cycle. Cyclic injection rate stimulation can slightly reduce the peak pressure, compared with monotonic stimulation, but only when the injection rate is low. The reduction in peak pressure occurs due to the combination of thermally-induced stresses associated with cooling and incremental damage rather than any influence of fatigue. Stepwise or low-frequency cyclic injection rate stimulation are suggested rather than a high-frequency cyclic injection rate stimulation, while injection with cyclic temperatures is suggested when more local damage is wanted.

This study presents a thermo-hydro-mechanical framework to model hydrothermal systems within a simplified faulted synthetic reservoir, replicating current production scenarios in The Netherlands and Germany. The reservoir, composed of porous and permeable sandstone, and the confining layer, made of porous but less permeable shale, undergoes a process where cold water is injected and hot water is extracted. A fault, situated 750 meters from the injection well, is investigated to examine the conditions when fault slip could occur. Various fault and formation stiffnesses are modeled to assess their impact on fault stability. Our analysis reveals that stress changes induced by hydrothermal operations can lead to fault reactivation, with the stiffness contrast between the reservoir and confining layers playing a significant role in when and where fault reactivation can occur. Stiffer confining layers lead to reactivation occurring more closely associated with the passage of the cooling front. In contrast, a stiffer reservoir results in greater and more gradual stress changes, making reactivation more closely related to the total volume of cooled rock.

This study presents a fully coupled hydro-mechanical framework for modeling hydraulic shearing in a mesoscale reservoir located at the Grimsel Test Site, Switzerland. The experiment was conducted on a ductile–brittle fault embedded in low-permeable granite. We observe that normal fracture opening increases flow channel recoverably, while fracture sliding locks asperities leading to a non-recoverable increase in flow. To couple these processes, we use a poro-elasto-plastic constitutive framework and employ a permeability function that depends on several parameters, such as dilation angle, in-situ stresses, residual aperture and maximum aperture. Our results capture the recorded pressure responses well, and indicate that the permeability changes by one order of magnitude during the experiment.