Abstract. In this contribution, we present a review of scientific research results that address seismo-hydromechanically coupled processes relevant for the development of a sustainable heat exchanger in low-permeability crystalline rock and introduce the design of the In situ Stimulation and Circulation (ISC) experiment at the Grimsel Test Site dedicated to studying such processes under controlled conditions. The review shows that research on reservoir stimulation for deep geothermal energy exploitation has been largely based on laboratory observations, large-scale projects and numerical models. Observations of full-scale reservoir stimulations have yielded important results. However, the limited access to the reservoir and limitations in the control on the experimental conditions during deep reservoir stimulations is insufficient to resolve the details of the hydromechanical processes that would enhance process understanding in a way that aids future stimulation design. Small-scale laboratory experiments provide fundamental insights into various processes relevant for enhanced geothermal energy, but suffer from (1) difficulties and uncertainties in upscaling the results to the field scale and (2) relatively homogeneous material and stress conditions that lead to an oversimplistic fracture flow and/or hydraulic fracture propagation behavior that is not representative of a heterogeneous reservoir. Thus, there is a need for intermediate-scale hydraulic stimulation experiments with high experimental control that bridge the various scales and for which access to the target rock mass with a comprehensive monitoring system is possible. The ISC experiment is designed to address open research questions in a naturally fractured and faulted crystalline rock mass at the Grimsel Test Site (Switzerland). Two hydraulic injection phases were executed to enhance the permeability of the rock mass. During the injection phases the rock mass deformation across fractures and within intact rock, the pore pressure distribution and propagation, and the microseismic response were monitored at a high spatial and temporal resolution.

Fracture network geometry is crucial for transport in hard rock aquifers, but it can only be approximated in models. While fracture orientation, spacing, and intensity can be obtained from borehole logs, core images, and outcrops, the characterization of in situ fracture network geometry requires the interpretation of spatially distributed hydraulic and transport experiments. In this study, we present a novel concept using a transdimensional inversion method (reversible jump Markov Chain Monte Carlo, rjMCMC) to invert a two-dimensional cross-well discrete fracture network (DFN) geometry from tracer tomography experiments. The conservative tracer transport is modeled via a fast finite difference model neglecting matrix diffusion. The proposed DFN inversion method iteratively evolves DFN variants by geometry updates to fit the observed tomographic data evaluated by the Metropolis-Hastings-Green acceptance criteria. A main feature is the varying dimensions of the inverse problem, which allows for the calibration of fracture geometries and numbers. This delivers an ensemble of thousands of DFN realizations that can be utilized for probabilistic identification of fractures in the aquifer. In the presented hypothetical and outcrop-based case studies, cross sections between boreholes are investigated. The procedure successfully identifies major transport pathways in the investigated domain and explores equally probable DFN realizations, which are analyzed in fracture probability maps and by multidimensional scaling.

Abstract. To characterize the stress field at the Grimsel Test Site (GTS) underground rock laboratory, a series of hydrofracturing and overcoring tests were performed. Hydrofracturing was accompanied by seismic monitoring using a network of highly sensitive piezosensors and accelerometers that were able to record small seismic events associated with metre-sized fractures. Due to potential discrepancies between the hydrofracture orientation and stress field estimates from overcoring, it was essential to obtain highprecision hypocentre locations that reliably illuminate fracture growth. Absolute locations were improved using a transverse isotropic P-wave velocity model and by applying joint hypocentre determination that allowed for the computation of station corrections. We further exploited the high degree of waveform similarity of events by applying cluster analysis and relative relocation. Resulting clouds of absolute and relative located seismicity showed a consistent east-west strike and 70 • dip for all hydrofractures. The fracture growth direction from microseismicity is consistent with the principal stress orientations from the overcoring stress tests, provided that an anisotropic elastic model for the rock mass is used in the data inversions. The σ 1 stress is significantly larger than the other two principal stresses and has a reasonably welldefined orientation that is subparallel to the fracture plane; σ 2 and σ 3 are almost equal in magnitude and thus lie on a circle defined by the standard errors of the solutions. The poles of the microseismicity planes also lie on this circle towards the north. Analysis of P-wave polarizations suggested double-couple focal mechanisms with both thrust and normal faulting mechanisms present, whereas strike-slip and thrust mechanisms would be expected from the overcoring-derived stress solution. The reasons for these discrepancies can be explained by pressure leak-off, but possibly may also involve stress field rotation around the propagating hydrofracture. Our study demonstrates that microseismicity monitoring along with high-resolution event locations provides valuable information for interpreting stress characterization measurements.

No abstract

Abstract. As part of the In-situ Stimulation and Circulation (ISC) experiment, hydraulic fracturing (HF) tests were conducted in a moderately fractured crystalline rock mass at the Grimsel Test Site (GTS), Switzerland. The aim of these injection tests was to improve our understanding of processes associated with high-pressure fluid injection. A total of six HF experiments were performed in two inclined boreholes, where the surrounding rock mass was accessed with twelve observation boreholes, which allow high-resolution monitoring of fracture fluid pressure, strain and micro-seismicity in an exceptionally well-characterized rock mass. A similar injection protocol was used for all six experiments to investigate the complexity of the fracture propagation processes. At the borehole scale, these processes involved newly created tensile fractures intersecting the injection interval while at the cross-hole scale, the natural network of fractures dominated the propagation process. The six HF experiments can be divided into two groups based on their injection location (i.e., south or north to a brittle ductile shear zone), their similarity of injection pressures and their response to deformation and pressure propagation. The injection tests performed in the south connect upon propagation to the brittle ductile shear zone. Thus, the shear zone acts as a dominant drain and a constant pressure boundary. The experiments executed north of the shear zone, show smaller injection pressures and larger backflow during bleed-off phases. From a seismic perspective, the injection tests show high variability in seismic response independent of the location of injection. For two injection experiments, we observe re-orientation of the seismic cloud as the fracture propagated away from the wellbore. In both cases, the main propagation direction is normal to the minimum principal stress direction. The re-orientation during propagation is interpreted to be related to a strong stress heterogeneity and the intersection of natural fractures striking different than the propagating hydraulic fracture. The seismic activity was limited to about 10 m radial distance from the injection point. In contrast, strain and pressure signals reach further into the rock mass indicating that the process zone around the injection point is larger than the zone illuminated by seismic signals. Furthermore, strain signals indicate not just single fracture openings but also the propagation of multiple fractures. Transmissivities of injection intervals increase about 2–4 orders of magnitudes.