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.

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.

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.

Acid treatments aimed at reducing formation breakdown pressure are becoming increasingly popular in stimulating deep and ultra-deep geothermal and natural gas reservoirs, where high pumping pressures are typically required for fracturing. This technique effectively reduces surface pumping pressures, ensuring safe operations. Through multiscale experiments and mechanistic analysis, this study reveals the fundamental mechanisms underlying breakdown pressure reduction via acid preconditioning. Key findings include: (1) High breakdown pressures arise from high in-situ stress, low porosity and permeability, and engineering contamination. (2) Acid-induced mineral dissolution triggers dual effects-pore structure evolution enhances reservoir permeability (facilitating subsequent fracturing fluid imbibition and pressure transmission), while mechanical property degradation substantially weakens rock resistance to fracturing. (3) Comparative analysis of HCl, organic acids, and chelating agents demonstrates that high-temperature reservoirs benefit from low-corrosivity chelating agents (e.g., GLDA) or organic acid systems combined with low-concentration HCl, achieving optimal dissolution efficiency while ensuring wellbore integrity. (4) A multiscale laboratory evaluation framework was established to integrate experimental data for optimizing acid formulations and post-acid fracturing strategies. This paper provides mechanistic insights, acid system selection criteria, and experimental methodologies for breakdown pressure reduction in deep carbonate reservoirs, offering significant engineering value for achieving safe and efficient reservoir stimulation.

Distributed acoustic sensing (DAS) that uses optical fibres as sensing units is attracting increasing interest for micro-seismic monitoring in geothermal projects. Standard optical fibres provide one-component measurements along the fibre and this pose challenges in determining certain characteristics of the source, such as its azimuth and its full moment tensor. Full source characteristics can be obtained via offset downhole measurements and/or measurements from horizontal well sections but these come with substantial extra costs. This paper proposes a single-well dual-cable DAS configuration to reduce the need for drilling additional wells or sections, where two DAS cables are assumed to be positioned within a single vertical well at opposite sides of the well. Synthetic DAS signals are generated by an open-source code that assumes plane-layered media and are used to study the feasibility of the dual-cable DAS for localising a seismic source and resolving its moment tensor. A localisation procedure is presented, and a sensitivity analysis of localisation accuracy is conducted with respect to source parameters and noise levels. In addition, an analysis is performed to assess the resolvability of the moment tensor components from the dual-cable DAS configuration. Results suggest the source location can be fully determined, yet low signal-to-noise ratio and azimuth close to 0∘ (North, aligned with the two cables) lead to a decrease in accuracy. The full moment tensor can be resolved only if the epicentral distance is 5 m or less, while non-double-couple components can be reliably resolved with an epicentral distance up to 20 m, showing improvement compared to installations with a single cable. Consequently, near-borehole failures, regardless of the source mechanisms, can be characterised within an epicentral distance of 5 m. With epicentral distance increasing, resolvability of the mix-mode failures is reduced first, followed by the resolvability of the pure shear or tensile failures, which depends on the azimuth. Overall, the results demonstrate that a single-well dual-cable configuration has the potential for monitoring and understanding near-borehole micro-seismic events induced during geothermal reinjection and stimulation operations.

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.

In geothermal projects, reinjection of produced water has been widely applied for disposing wastewater, supplying heat exchange media and maintaining reservoir pressure. Accordingly, it is a key process for environmental and well performance assessment, which partly controls the success of projects. However, the injectivity, a measure of how easily fluids can be reinjected into reservoirs, is influenced by various processes throughout installation and operation. Both injectivity decline and enhancement have been reported during reinjection operations, while most current studies tend to only focus on one aspect. This review aims to provide a comprehensive discussion on how the injectivity can be influenced during reinjection, both positively and negatively. This includes a detailed overview of the different clogging mechanisms, in which decreasing reservoir temperature plays a major role, leading to injectivity decline. Strategies to avoid and recover from injectivity reduction are also introduced. Followed is an overview of mechanisms underlying injectivity enhancement during reinjection, wherein re-opening/shearing of pre-existing fractures and thermal cracking have been identified as the main contributors. In practice, nevertheless, mixed-mechanism processes play a key role during reinjection. Finally, this review provides an outlook on future research directions that can enhance the understanding of injectivity-related issues.