ALBERT

Description: We performed a series of 12 hydraulic stimulation experiments in a 20m×20m×20m foliated, crystalline rock volume intersected by two distinct fault sets at the Grimsel Test Site, Switzerland. The goal of these experiments was to improve our understanding of stimulation processes associated with high-pressure fluid injection used for reservoir creation in enhanced or engineered geothermal systems. In the first six experiments, pre-existing fractures were stimulated to induce shear dilation and enhance permeability. Two types of shear zones were targeted for these hydroshearing experiments: (i) ductile ones with intense foliation and (ii) brittle–ductile ones associated with a fractured zone. The second series of six stimulations were performed in borehole intervals without natural fractures to initiate and propagate hydraulic fractures that connect the wellbore to the existing fracture network. The same injection protocol was used for all experiments within each stimulation series so that the differences observed will give insights into the effect of geology on the seismo-hydromechanical response rather than differences due to the injection protocols. Deformations and fluid pressure were monitored using a dense sensor network in boreholes surrounding the injection locations. Seismicity was recorded with sensitive in situ acoustic emission sensors both in boreholes and at the tunnel walls. We observed high variability in the seismic response in terms of seismogenic indices, b values, and spatial and temporal evolution during both hydroshearing and hydrofracturing experiments, which we attribute to local geological heterogeneities. Seismicity was most pronounced for injections into the highly conductive brittle–ductile shear zones, while the injectivity increase on these structures was only marginal. No significant differences between the seismic response of hydroshearing and hydrofracturing was identified, possibly because the hydrofractures interact with the same pre-existing fracture network that is reactivated during the hydroshearing experiments. Fault slip during the hydroshearing experiments was predominantly aseismic. The results of our hydraulic stimulations indicate that stimulation of short borehole intervals with limited fluid volumes (i.e., the concept of zonal insulation) may be an effective approach to limit induced seismic hazard if highly seismogenic structures can be avoided.

Description: In this contribution we present a review of scientific research results that address seismo-hydro-mechanical 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 study 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 hydro-mechanical processes that would enhance process understanding in a way that aids future stimulation design. Small scale laboratory experiments provide a 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 over-simplistic fracture flow and/or hydraulic fracture propagation behaviour that is not representative for 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. Only few intermediate-scale hydro-shearing and hydro-fracturing experiments have recently been performed in a densely instrumented rock mass. No such measurements have been performed on faults in crystalline basement rocks. The In-situ Stimulation and Circulation (ISC) experiment currently performed in a naturally fractured and faulted crystalline rock mass at the Grimsel Test Site (Switzerland) is designed to address open research questions, which could not be investigated in the required detail so far. Two hydraulic injection phases were executed to enhance the permeability of the rock mass: a hydro-shearing phase and then a hydraulic fracturing phase. During the injection phases the rock mass deformation across fractures and within intact rock, the pore pressure distribution and propagation and the micro-seismic response were monitored at a high spatial and temporal resolution.

Description: 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.

Description: 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 high-precision 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 well-defined 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.

Description: To characterize the stress field at the Grimsel Test Site (GTS) underground rock laboratory a series of hydrofracturing test and overcoring test were performed. Hydrofracturing was accompanied by seismic monitoring using a network of highly sensitive piezo sensors and accelerometers that were able to record small seismic events associated with decimeter-sized fractures. Due to potential discrepancies between the hydro-fracture orientation and stress field estimates from overcoring, it was essential to obtain high-precision hypocenter locations that reliably illuminate fracture growth. Absolute locations were improved using a transverse isotropic P-wave velocity model and by applying joint hypocenter determination that allowed 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 hydro-fractures. The fracture growth direction from microseismicity is consistent with the principal stress orientations from the overcoring stress tests provided an anisotropic elastic model for the rock mass is used in the data inversions. σ1 is significantly larger than the other two principal stresses, and has a reasonably well-defined 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. The trace of the hydraulic fracture imaged at the borehole wall show that they initiated within the foliation plane, which differs in orientation from the microseismicity planes. Thus, fracture initiation was most likely influenced by a foliation-related strength anisotropy. 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 are not well understood, but may 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.

Description: AlpArray is a large collaborative seismological project in Europe that includes more than 50 research institutes and seismological observatories. At the heart of the project is the collection of top-quality seismological data from a dense network of broadband temporary seismic stations, in compliment to the existing permanent networks, that ensures a homogeneous station coverage of the greater Alpine region. This Alp Array Seismic Network (AASN) began operation in January 2016 and will have a duration of at least 2 years. In this work we report the Swiss contribution to the AASN, we concentrate on the site selection process, our methods for stations installation, data quality and data management. We deployed 27 temporary broadband stations equipped with STS-2 and Trillium Compact 120 s sensors. The deployment and maintenance of the temporary stations across 5 countries is managed by ETH Zurich and it is the result of a fruitful collaboration between five institutes in Europe.

Description: The success of geological carbon storage depends on the assurance of permanent containment for injected carbon dioxide (CO2) in the storage formation at depth. One of the critical elements of the safekeeping of CO2 is the sealing capacity of the caprock overlying the storage formation despite faults and/or fractures, which may occur in it. In this work, we present an ongoing injection experiment performed in a fault hosted in clay at the Mont Terri underground rock laboratory (NW Switzerland). The experiment aims to improve our understanding of the main physical and chemical mechanisms controlling (i) the migration of CO2 through a fault damage zone, (ii) the interaction of the CO2 with the neighboring intact rock, and (iii) the impact of the injection on the transmissivity in the fault. To this end, we inject CO2-saturated saline water in the top of a 3 m thick fault in the Opalinus Clay, a clay formation that is a good analog of common caprock for CO2 storage at depth. The mobility of the CO2 within the fault is studied at the decameter scale by using a comprehensive monitoring system. Our experiment aims to close the knowledge gap between laboratory and reservoir scales. Therefore, an important aspect of the experiment is the decameter scale and the prolonged duration of observations over many months. We collect observations and data from a wide range of monitoring systems, such as a seismic network, pressure temperature and electrical conductivity sensors, fiber optics, extensometers, and an in situ mass spectrometer for dissolved gas monitoring. The observations are complemented by laboratory data on collected fluids and rock samples. Here we show the details of the experimental concept and installed instrumentation, as well as the first results of the preliminary characterization. An analysis of borehole logging allows for identifying potential hydraulic transmissive structures within the fault zone. A preliminary analysis of the injection tests helped estimate the transmissivity of such structures within the fault zone and the pressure required to mechanically open such features. The preliminary tests did not record any induced microseismic events. Active seismic tomography enabled sharp imaging the fault zone.

Description: In this study, we analyzed a large seismological dataset from temporary and permanent networks in the southern and eastern Alps to establish high-precision hypocenters and 1-D VP and VP/VS models. The waveform data of a subset of local earthquakes with magnitudes in the range of 1–4.2 ML were recorded by the dense, temporary SWATH-D network and selected stations of the AlpArray network between September 2017 and the end of 2018. The first arrival times of P and S waves of earthquakes are determined by a semi-automatic procedure. We applied a Markov chain Monte Carlo inversion method to simultaneously calculate robust hypocenters, a 1-D velocity model, and station corrections without prior assumptions, such as initial velocity models or earthquake locations. A further advantage of this method is the derivation of the model parameter uncertainties and noise levels of the data. The precision estimates of the localization procedure is checked by inverting a synthetic travel time dataset from a complex 3-D velocity model and by using the real stations and earthquakes geometry. The location accuracy is further investigated by a quarry blast test. The average uncertainties of the locations of the earthquakes are below 500 m in their epicenter and ∼ 1.7 km in depth. The earthquake distribution reveals seismicity in the upper crust (0–20 km), which is characterized by pronounced clusters along the Alpine frontal thrust, e.g., the Friuli-Venetia (FV) region, the Giudicarie–Lessini (GL) and Schio-Vicenza domains, the Austroalpine nappes, and the Inntal area. Some seismicity also occurs along the Periadriatic Fault. The general pattern of seismicity reflects head-on convergence of the Adriatic indenter with the Alpine orogenic crust. The seismicity in the FV and GL regions is deeper than the modeled frontal thrusts, which we interpret as indication for southward propagation of the southern Alpine deformation front (blind thrusts).

Description: The recent development of ambient noise tomography, in combination with the increasing number of permanent seismic stations and dense networks of temporary stations operated during passive seismic experiments, provides a unique opportunity to build the first high-resolution 3-D shear wave velocity (vS) model of the entire crust of the Bohemian Massif (BM). This paper provides a regional-scale model of velocity distribution in the BM crust. The velocity model with a cell size of 22 km is built using a conventional two-step inversion approach from Rayleigh wave group velocity dispersion curves measured at more than 400 stations. The shear velocities within the upper crust of the BM are ∼0.2 km s−1 higher than those in its surroundings. The highest crustal velocities appear in its southern part, the Moldanubian unit. The Cadomian part of the region has a thinner crust, whereas the crust assembled, or tectonically transformed in the Variscan period, is thicker. The sharp Moho discontinuity preserves traces of its dynamic development expressed in remnants of Variscan subductions imprinted in bands of crustal thickening. A significant feature of the presented model is the velocity-drop interface (VDI) modelled in the lower part of the crust. We explain this feature by the anisotropic fabric of the lower crust, which is characterised as vertical transverse isotropy with the low velocity being the symmetry axis. The VDI is often interrupted around the boundaries of the crustal units, usually above locally increased velocities in the lowermost crust. Due to the north-west–south-east shortening of the crust and the late-Variscan strike-slip movements along the north-east–south-west oriented sutures preserved in the BM lithosphere, the anisotropic fabric of the lower crust was partly or fully erased along the boundaries of original microplates. These weakened zones accompanied by a velocity increase above the Moho (which indicate an emplacement of mantle rocks into the lower crust) can represent channels through which portions of subducted and later molten rocks have percolated upwards providing magma to subsequently form granitoid plutons.

Description: The region around the town Albstadt, SW Germany, was struck by four damaging earthquakes with magnitudes greater than 5 during the last century. These earthquakes occurred along the Albstadt Shear Zone (ASZ), which is characterized by more or less continuous microseismicity. As there are no visible surface ruptures that may be connected to the fault zone, we study its characteristics by its seismicity distribution and faulting pattern. We use the earthquake data of the state earthquake service of Baden-Württemberg from 2011 to 2018 and complement it with additional phase picks beginning in 2016 at the AlpArray and StressTransfer seismic networks in the vicinity of the ASZ. This extended data set is used to determine new minimum 1-D seismic vp and vs velocity models and corresponding station delay times for earthquake relocation. Fault plane solutions are determined for selected events, and the principal stress directions are derived. The minimum 1-D seismic velocity models have a simple and stable layering with increasing velocity with depth in the upper crust. The corresponding station delay times can be explained well by the lateral depth variation of the crystalline basement. The relocated events align about north–south with most of the seismic activity between the towns of Tübingen and Albstadt, east of the 9∘ E meridian. The events can be separated into several subclusters that indicate a segmentation of the ASZ. The majority of the 25 determined fault plane solutions feature an NNE–SSW strike but NNW–SSE-striking fault planes are also observed. The main fault plane associated with the ASZ dips steeply, and the rake indicates mainly sinistral strike-slip, but we also find minor components of normal and reverse faulting. The determined direction of the maximum horizontal stress of 140–149∘ is in good agreement with prior studies. Down to ca. 7–8 km depth SHmax is bigger than SV; below this depth, SV is the main stress component. The direction of SHmax indicates that the stress field in the area of the ASZ is mainly generated by the regional plate driving forces and the Alpine topography.