A laboratory study of hydraulic fracturing at the brittle-ductile transition.

Developing high-enthalpy geothermal systems requires a sufficiently permeable formation to extract energy through fluid circulation. Injection experiments above water's critical point have shown that fluid flow can generate a network of highly conductive tensile cracks. However, what remains unclear is the role played by fluid and solid rheology on the formation of a dense crack network. The decrease of fluid viscosity with temperature and the thermally activated visco-plasticity in rock are expected to change the deformation mechanisms and could prevent the formation of fractures. To isolate the solid rheological effects from the fluid ones and the associated poromechanics, we devise a hydro-fracture experimental program in a non-porous material, polymethyl methacrylate (PMMA). In the brittle regime, we observe rotating cracks and complex fracture patterns if a non-uniform stress distribution is introduced in the samples. We observe an increase of ductility with temperature, hampering the propagation of hydraulic fractures close to the glass transition temperature of PMMA, which acts as a limit for brittle fracture propagation. Above the glass transition temperature, acoustic emission energy drops of several orders of magnitude. Our findings provide a helpful guidance for future studies of hydro-fracturing of supercritical geothermal systems.

How Equivalent Are Equivalent Porous Media?

Geoenergy and geoengineering applications usually involve fluid injection into and production from fractured media. Accounting for fractures is important because of the strong poromechanical coupling that ties pore pressure changes and deformation. A possible approach to the problem uses equivalent porous media to reduce the computational cost and model complexity instead of explicitly including fractures in the models. We investigate the validity of this simplification by comparing these two approaches. Simulation results show that pore pressure distribution significantly differs between the two approaches even when both are calibrated to predict identical values at the injection and production wells. Additionally, changes in fracture stability are not well captured with the equivalent porous medium. We conclude that explicitly accounting for fractures in numerical models may be necessary under some circumstances to perform reliable coupled thermohydromechanical simulations, which could be used in conjunction with other tools for induced seismicity forecasting.

Analytical Solution to Assess the Induced Seismicity Potential of Faults in Pressurized and Depleted Reservoirs.

Displaced faults crossing the reservoir could significantly increase the induced earthquake frequency in geo-energy projects. Understanding and predicting the stress variation in such cases is essential to minimize the risk of induced seismicity. Here, we adopt the inclusion theory to develop an analytical solution for the stress response to pore pressure variations within the reservoir for both permeable and impermeable faults with offset ranging from zero to the reservoir thickness. By analyzing fault stability changes due to reservoir pressurization/depletion under different scenarios, we find that (1) the induced seismicity potential of impermeable faults is always larger than that of permeable faults under any initial and injection conditions-the maximum size of the fault undergoing failure is 3-5 times larger for impermeable than for permeable faults; (2) stress concentration at the corners results in the occurrence of reversed slip in normal faults with a normal faulting stress regime; (3) while fault offset has no impact on the slip potential for impermeable faults, the slip potential increases with the offset for permeable faults, which indicates that non-displaced permeable faults constitute a safer choice for site selection; (4) an impermeable fault would rupture at a lower deviatoric stress, and at a smaller pressure buildup than a permeable one; and (5) the induced seismicity potential is overestimated and the injectivity underestimated if the stress arching (i.e., the poromechanical coupling) is neglected. This analytical solution is a useful tool for site selection and for supporting decision making during the lifetime of geo-energy projects.

Sinking CO2 in Supercritical Reservoirs.

Geologic carbon storage is required for achieving negative CO2 emissions to deal with the climate crisis. The classical concept of CO2 storage consists in injecting CO2 in geological formations at depths greater than 800 m, where CO2 becomes a dense fluid, minimizing storage volume. Yet CO2 has a density lower than the resident brine and tends to float, challenging the widespread deployment of geologic carbon storage. Here, we propose for the first time to store CO2 in supercritical reservoirs to reduce the buoyancy-driven leakage risk. Supercritical reservoirs are found at drilling-reachable depth in volcanic areas, where high pressure (p > 21.8 MPa) and temperature (T > 374°C) imply CO2 is denser than water. We estimate that a CO2 storage capacity in the range of 50-500 Mt yr-1 could be achieved for every 100 injection wells. Carbon storage in supercritical reservoirs is an appealing alternative to the traditional approach.

The risks of long-term re-injection in supercritical geothermal systems.

Supercritical geothermal systems are appealing sources of sustainable and carbon-free energy located in volcanic areas. Recent successes in drilling and exploration have opened new possibilities and spiked interest in this technology. Experimental and numerical studies have also confirmed the feasibility of creating fluid conducting fractures in sedimentary and crystalline rocks at high temperature, paving the road towards Enhanced Supercritical Geothermal Systems. Despite their attractiveness, several important questions regarding safe exploitation remain open. We dedicate this manuscript to the first thermo-hydro-mechanical numerical study of a doublet geothermal system in supercritical conditions. Here we show that thermally-induced stress and strain effects dominate the geomechanical response of supercritical systems compared to pore pressure-related instabilities, and greatly enhance seismicity during cold water re-injection. This finding has important consequences in the design of Supercritical Geothermal Systems.

The brittle-ductile transition in active volcanoes.

Contrasting deformation mechanisms precede volcanic eruptions and control precursory signals. Density increase and high uplifts consistent with magma intrusion and pressurization are in contrast with dilatant responses and reduced surface uplifts observed before eruptions. We investigate the impact that the rheology of rocks constituting the volcanic edifice has on the deformation mechanisms preceding eruptions. We propose a model for the pressure and temperature dependent brittle-ductile transition through which we build a strength profile of the shallow crust in two idealized volcanic settings (igneous and sedimentary basement). We have performed finite element analyses in coupled thermo-hydro-mechanical conditions to investigate the influence of static diking on the local brittle-ductile transition. Our results show that in active volcanoes: (i) dilatancy is an appropriate indicator for the brittle-ductile transition; (ii) the predicted depth of the brittle-ductile transition agrees with the observed attenuated seismicity; (iii) seismicity associated with diking is likely to be affected by ductile deformation mode caused by the local temperature increase; (iv) if failure occurs within the edifice, it is likely to be brittle-dilatant with strength and stiffness reduction that blocks stress transfers within the volcanic edifice, ultimately damping surface uplifts.

Tensor Field Visualization using Fiber Surfaces of Invariant Space.

Scientific visualization developed successful methods for scalar and vector fields. For tensor fields, however, effective, interactive visualizations are still missing despite progress over the last decades. We present a general approach for the generation of separating surfaces in symmetric, second-order, three-dimensional tensor fields. These surfaces are defined as fiber surfaces of the invariant space, i.e. as pre-images of surfaces in the range of a complete set of invariants. This approach leads to a generalization of the fiber surface algorithm by Klacansky et al. [16] to three dimensions in the range. This is due to the fact that the invariant space is three-dimensional for symmetric second-order tensors over a spatial domain. We present an algorithm for surface construction for simplicial grids in the domain and simplicial surfaces in the invariant space. We demonstrate our approach by applying it to stress fields from component design in mechanical engineering.