High fluid pressure and triggered earthquakes in the enhanced geothermal system in Basel, Switzerland
|
|
- Wilfred Blair
- 6 years ago
- Views:
Transcription
1 JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 117,, doi: /2011jb008980, 2012 High fluid pressure and triggered earthquakes in the enhanced geothermal system in Basel, Switzerland Toshiko Terakawa, 1 Stephen A. Miller, 2 and Nicholas Deichmann 3 Received 28 October 2011; revised 15 May 2012; accepted 23 May 2012; published 20 July [1] We analyzed 118 well-constrained focal mechanisms to estimate the pore fluid pressure field of the stimulated region during the fluid injection experiment in Basel, Switzerland. This technique, termed focal mechanism tomography (FMT), uses the orientations of slip planes within the prevailing regional stress field as an indicator of the fluid pressure along the plane at the time of slip. The maximum value and temporal change of excess pore fluid pressures are consistent with the known history of the wellhead pressure applied at the borehole. Elevated pore fluid pressures were concentrated within 500 m of the open hole section, which are consistent with the spatiotemporal evolution of the induced microseismicity. Our results demonstrate that FMT is a robust approach, being validated at the meso-scale of the Basel stimulation experiment. We found average earthquake triggering excess pore fluid pressures of about 10 MPa above hydrostatic. Overpressured fluids induced many small events (M < 3) along faults unfavorably oriented relative to the tectonic stress pattern, while the larger events tended to occur along optimally oriented faults. This suggests that small-scale hydraulic networks, developed from the high pressure stimulation, interact to load (hydraulically isolated) high strength bridges that produce the larger events. The triggering pore fluid pressures are substantially higher than that predicted from a linear pressure diffusion process from the source boundary, and shows that the system is highly permeable along flow paths that allow fast pressure diffusion to the boundaries of the stimulated region. Citation: Terakawa, T., S. A. Miller, and N. Deichmann (2012), High fluid pressure and triggered earthquakes in the enhanced geothermal system in Basel, Switzerland, J. Geophys. Res., 117,, doi: /2011jb Introduction [2] Hydraulic stimulation is one of the fundamental methods to increase the permeability at depth in the enhanced geothermal system (EGS). From 2 to 8 December, 2006, a private/public consortium in Basel, Switzerland (Figure 1), attempted to stimulate a geothermal reservoir by injecting approximately 11,500 m 3 of water at high pressure into a 5 km deep well[häring et al., 2008]. Thousands of induced seismic events were recorded by a six-sensor borehole array installed at depths between 317 m and 2740 m around the well. At a peak wellhead pressure of almost 30 MPa, a seismic event of M L 2.6 occurred, that surpassed the safety threshold and which resulted in a termination of the injection. It was followed a few hours later by an M L 3.4 event. Seismic activity in the stimulated rock volume declined rapidly over the following 1 Earthquake and Volcano Research Center, Graduate School of Environmental Studies, Nagoya University, Nagoya, Japan. 2 Geodynamics/Geophysics, Steinmann-Institute, University of Bonn, Bonn, Germany. 3 Swiss Seismological Service, ETH-Zürich, Zürich, Switzerland. Corresponding author: T. Terakawa, Earthquake and Volcano Research Center, Nagoya University, Nagoya , Japan. (terakawa@seis.nagoya-u.ac.jp) American Geophysical Union. All Rights Reserved /12/2011JB three weeks, but increased again in January and February 2007, including three additional events with M L >3. [3] Determining the pore fluid pressure field in the stimulated region of EGS is essential for the long-term management of geothermal resources. Focal mechanisms of seismic events reflect the pore fluid pressure and the friction coefficient on the fault as well as the tectonic stress pattern, because fault strength is controlled by Coulomb failure criterion (Coulomb friction law). Several studies [e.g., Zoback, 1992b; Cornet and Yin, 1995; Plenefisch and Bonjer, 1997; Rivera and Kanamori, 2002] estimated a uniform stress pattern from a variety of focal mechanisms of seismic events, and investigated levels of the pore fluid pressure and the friction coefficient using the following two end-member models. One endmember model assumes hydrostatic pore fluid pressure, and attributes focal mechanism variations to variations in friction coefficients of pre-existing faults. The other end-member model assumes a constant standard friction coefficient (0.6), as suggested by Byerlee s law [Byerlee, 1978]. This model attributes the observed variation of focal mechanisms to fault strength heterogeneity controlled by variations in pore fluid pressures acting on faults. Using Byerlee s law as well as Bayesian statistical inference and Akaike s Bayesian Information Criterion (ABIC) [Akaike, 1980], Terakawa et al. [2010] developed an analysis technique, termed focal mechanisms tomography (FMT), to estimate the 3D pore fluid 1of15
2 Figure 1. Topographic map of Switzerland and surroundings with regional seismic stations that contributed data to the focal mechanisms. For the locations of the seismic stations in the region of Basel see Deichmann and Giardini [2009] or Deichmann and Ernst [2009]. pressure field from focal mechanisms. Application of this technique to the 2009 L Aquila earthquake in Italy showed that overpressured fluid reservoirs at hypocentral depths likely contributed to that sequence. The results of Terakawa et al. [2010] were consistent with other seismic studies that showed the involvement of high pressure fluids [Lucente et al., 2010], but there was no independent way to validate the pore fluid pressures determined using FMT. [4] The Basel fluid-injection experiment provided an excellent opportunity to validate the FMT approach because of the well-constrained regional stress field, the known fluid pressure history applied at the open section of the borehole, and numerous well-constrained focal mechanisms with relocated hypocenters [e.g., Häring et al., 2008; Valley and Evans, 2006, 2009; Deichmann and Giardini, 2009]. In this paper, we apply FMT to the Basel experiment to estimate the pore fluid pressure in the stimulated region and quantitatively demonstrate the validity of the FMT method. We also discuss how the results provide additional insight into the role of high pressure fluids and their interaction with tectonic stress to earthquake generation. 2. Methodology of Focal Mechanism Tomography [5] Focal mechanism tomography is an analysis technique to estimate pore fluid pressures at depth from focal mechanisms of seismic events [Terakawa et al., 2010]. The primary assumptions in this approach are; 1) fault strength is controlled by the Coulomb failure criterion with a constant standard friction coefficient [Byerlee, 1978], and 2) seismic slip occurs in the direction of the resolved shear traction acting on pre-existing faults [Wallace, 1951; Bott, 1959]. We further assume that 3) the vertical stress is the weight of the overburden, and 4) seismic slip on optimally oriented faults relative to the prevailing regional stress pattern occurs under hydrostatic pressure. The following sections describe basic ideas and the procedure of the method Basic Ideas [6] We assume that earthquakes are governed by the Coulomb failure criterion (assumption 1). The fault strength t s is described by t s ¼ mðs n P f Þ where s n and P f are the normal stress (compression is positive) and the pore fluid pressure, respectively, and m is the friction coefficient. Within a uniform stress pattern and hydrostatic pore fluid pressure, the variety of earthquake focal mechanisms can be attributed to variations in friction coefficients of pre-existing faults. Conversely, in a uniform stress pattern with constant friction coefficients, variations in focal mechanisms can be attributed to fault strength heterogeneity controlled by variations in pore fluid pressures ð1þ 2of15
3 Figure 2. Theory of FMT in 3D Mohr diagram. The horizontal and vertical axes show normal and shear stresses; s 1, s 2, and s 3 are the maximum, intermediate, and minimum compressive principal stresses; P f, P h, r w, g, z, and m are the pore fluid pressure, the hydrostatic pore fluid pressure, density of water, gravity, depth of hypocenter, and the friction coefficient. A stress state acting on a plane (s 0 : normal stress, t 0 : shear stress) is shown with a point within the shadow zone surrounded by the three Mohr circles. The solid line labeled H 0 shows the fault strength controlled by Coulomb failure criterion under hydrostatic pressure. Assuming that rock tensile strength is negligible, the fluid pressure can increase to the minimum principal stress, and the fault strength is shown by the solid line labeled H max. Elevated pore fluid pressure reduces fault strength and triggers seismic slip on faults unfavorably oriented relative to the regional stress pattern. acting on pre-existing faults. Based on laboratory experiments of rock friction [e.g., Byerlee, 1978] and in situ stress measurements in deep boreholes [e.g., Zoback and Healy, 1992; Zoback and Townend, 2001], heterogeneity of friction coefficients within the rock mass is normally be substantially smaller than the pore fluid pressure heterogeneity in the study region. In FMT we estimate pore fluid pressure with a constant friction coefficient. [7] The theory of FMT can be explained using the 3D Mohr diagram (Figure 2). This diagram represents the stress state on any plane within a uniform stress field, designated by a point in the region surrounded by the three Mohr circles. The horizontal and vertical coordinates of the point show the normal and shear stresses acting on the plane, while the relative position of the point shows the orientation of the fault relative to the stress pattern [e.g., Zoback, 2007]. Since seismic slip occurs when shear stress reaches the fault strength, the intersection of the Mohr-Coulomb failure line passing through this point and the horizontal axis shows the pore fluid pressure P f activating the event. Thus elevated fluid pressure initiates seismic slip on unfavorably oriented faults relative to the regional stress pattern even in the absence of shear stress build-up, and the fault orientation within that stress field is a measure of pore fluid pressure. It should be noted that pore fluid pressure activating seismic slip must be theoretically less than the least principal stress, otherwise hydraulic fracture occurs [e.g., Zoback, 1992b; Rice, 1992; Hardebeck and Michael, 2004] The Procedure of FMT [8] In the first step we estimate the tectonic stress pattern, which is usually represented by the 3 directions of the principal stress axes and the R value, R = (s 1 s 2 )/ (s 1 s 3 ) (i.e., the relative dimensions of the Mohr circles), through stress inversion of seismic events [e.g., Gephart and Forsyth, 1984; Michael, 1984, 1987; Horiuchi et al., 1995; Terakawa and Matsu ura, 2008]. The signs of s 1, s 2, and s 3 show maximum, intermediate, and minimum (compressive) principal stresses, respectively. [9] In the second step we identify which of the two nodal planes is a geometrically possible slip plane. We calculate expected moment tensors from the orientation (strike and dip angles) of both nodal planes [Terakawa and Matsu ura, 2008], based on the assumption 2 (the Wallace-Bott hypothesis). We measure the closeness of the observed and expected moment tensors (CT) with inner tensor product as defined by Michael [1987]: where m o ij X 3 i¼1 X 3 m o ij me ij j¼1 CT ¼ vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffivffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ; ð2þ u X 3 X t 3 u X 3 X t 3 and m e ij i¼1 ðm o ij Þ2 j¼1 i¼1 ðm e ij Þ2 j¼1 are components of the observed and expected moment tensors. The CT value ranges from 1 CT 1, where +1 indicates that the observed moment tensor exactly coincides with the expected moment tensor, and 1 indicating that the observed moment tensor is opposite to that expected. We choose the nodal plane with the higher CT value as the true fault plane. It should be noted that this criterion for selection of true fault planes depends only on the stress pattern but not on the absolute stress tensor. [10] In order to determine absolute pore fluid pressures, we need to know all the six components of the stress tensor. In the third step we determine the deviatoric and isotropic stress levels (i.e., the actual diameters and the center of the largest Mohr circle, respectively) with assumptions 3 and 4. We can determine the isotropic level by assumption 3 that the vertical stress is the weight of the overburden. Assumption 4 requires that the Mohr-Coulomb failure line (under hydrostatic pore fluid pressure) must connect with the largest Mohr circle, which enables a determination of the diameters of the Mohr circles. [11] Finally, from Cauchy s relation and given an absolute stress tensor and fault orientations (determined from focal mechanisms in step 2), we can estimate normal and shear stresses acting on the fault plane. Pore fluid pressures are then determined by solving equation (1). Determining pore fluid pressure from this method results in a pore fluid pressure measurement at the discrete points of the hypocenters. By applying an inversion technique developed by Yabuki and Matsu ura [1992] to the excess (above hydrostatic) pore fluid pressures p n (n =1,, N) of the data set, we can estimate the 3D excess pore fluid pressure field as a continuous function with estimation errors (see Appendix). This technique is based on Bayesian statistic inference and Akaike s Bayesian Information Criterion (ABIC) [Akaike, 1980]. 3. FMT Analysis in the Basel Geothermal Reservoir [12] We applied the FMT method to focal mechanisms based on first-motion polarities (Figure 3) determined for 3of15
4 Figure 3. Hypocenter maps of 118 induced events and focal mechanisms of the 42 events with ML 1.7. (a) Map of epicenters (gray dots). Focal mechanisms are represented as lower-hemisphere stereographic projections, scaled with magnitude. The color of the focal spheres corresponds to focal depth. (b) Vertical distribution of hypocenters (E-W section). (c) Vertical distribution of hypocenters (N-S section). The color of the dots corresponds to focal depth. The focal depth is measured from sea level. 118 events with local magnitudes (M L ) between 1.0 and 3.4, which occurred from 3 December 2006 to 30 November The focal mechanisms were computed as described in Deichmann and Ernst [2009]. According to the classification of Zoback [1992a], of the 118 focal mechanisms, 101 are strike-slip, 5 are oblique strike-slip with a normal component, and 12 are normal faulting mechanisms (Figure 4b). Hypocentral locations of the events analyzed in this study were determined using a master-event relocation technique. As documented in Deichmann and Giardini [2009] mean and standard deviations of the calculated relative location errors for the x, y, and z coordinates are 49 4m,54 6 m, and 72 8 m. In Section 3.1, following the procedure described in Section 2.2, we first determined the tectonic stress field from the data set. Next, in Section 3.2, we estimated the pore fluid pressure distribution at depth by examining the fault orientation relative to the tectonic stress pattern Tectonic Stress Field in the Basel Geothermal Reservoir [13] Using the computer code of Terakawa et al. [2010], from the data set we searched the best fitting stress pattern which minimizes the difference between the expected and observed moment tensors. We chose the nodal planes with the higher CT values between the observed and expected moment tensors as the true fault planes (Table 1). [14] Figure 4a shows the tectonic stress pattern and the estimation errors of the three principal stress axes estimated in the present study. We evaluated the estimation errors of the stress pattern using a bootstrap method. The number of re-sampling was 500. The stress pattern is characterized by strike-slip faulting, where the maximum and minimum principal stress axes are almost horizontal in the directions of N38 WandN48 E and the value of R is The 67% confidence regions of maximum and minimum horizontal principal stress axes are N30 W-N45 WandN41 E-N59 E. R-values in the confidence region vary from 0.3 to 0.6. This stress pattern is entirely consistent with that deduced by Valley and Evans [2006, 2009] from borehole measurements (maximum horizontal principal stress: N36 W 14,minimum horizontal principal stress: N54 E 14 ). For the Basel data set, the average CT is 0.98, which corresponds to an average misfit angle between observed and expected slip vectors of about 9.7 degrees. Seismic events with high CT values can be explained as reactivation of pre-existing faults in a uniform stress field, while those with low CT may be indicative of local stress heterogeneities and/or very weak fault strengths. The horizontal orientation of the principle stress axes coincides to within 4 degrees with the average azimuth of the P- and T-axes (Figure 4b). [15] With the known stress pattern, and assumed friction coefficient of 0.85 [Byerlee, 1978] and rock density of 2540 kg/m 3 [Häring et al., 2008], we determined the absolute stress tensor (Figure 5). The absolute values of s 1, s 2, and s 3 at the depth of 5 km are 144 MPa, 117 MPa and Figure 4. The tectonic stress pattern and focal mechanisms of the triggered seismicity. (a) The tectonic stress pattern. (b) Stereographic projection (equal area, lower hemisphere) of the orientation of the P-axes (white circles) and T-axes (black circles) of the 118 focal mechanisms of the triggered Basel events. The diagonal lines correspond to the average azimuths of the P- and T-axes. (c) Symmetric polar histogram of the strike of all nodal planes, together with the overall strike of the macroscopic structure defined by the hypocenter locations. The arrows in Figures 4b and 4c point in the direction of the maximum and minimum horizontal compressive stress estimated in the present study. The red line in Figure 4c shows the orientation of the macroscopic plane of microseismic cloud. 4of15
5 Table 1. Focal Mechanism Parameters a Event Number Lat Lon Z M Plane 1 f 1 CT 1 Pf1 Plane 2 f 2 CT 2 Pf /63/ /86/153/ /62/ /84/152/ /62/ /84/151/ /84/ /73/5/ /63/ /86/152/ /90/ /89/0/ /84/ /77/5/ /82/ /63/7/ /84/ /78/6/ /61/ /85/151/ /42/ /52/ 112/ /78/ /89/ 12/ /46/ /49/ 116/ /88/ /66/ 2/ /47/ /68/ 132/ /44/ /47/ 105/ /65/ /84/154/ /45/ /51/ 115/ /41/ /53/ 112/ /44/ /50/ 114/ /66/ /84/155/ /52/ /46/ 123/ /83/ /57/ 7/ /44/ /50/ 114/ /70/ /74/ 20/ /49/ /45/ 114/ /68/ /62/ 24/ /44/ /50/ 114/ /87/ /75/3/ /79/ /66/ 168/ /71/ /65/ 20/ /71/ /66/ 19/ /79/ /58/12/ /87/ /80/ 2/ /45/ /47/ 108/ /49/ /45/ 114/ /70/ /76/159/ /86/ /77/176/ /72/ /63/ 19/ /83/ /88/172/ /60/ /79/149/ /85/ /88/174/ /88/ /80/178/ /78/ /72/11/ /80/ /88/ 10/ /72/ /63/ 19/ /84/ /80/5/ /87/ /81/177/ /75/ /66/15/ /87/ /81/177/ /80/ /73/170/ /79/ /77/ 10/ /68/ /86/157/ /74/ /75/ 16/ /76/ /50/ 162/ /87/ /72/ 2/ /83/ /72/173/ /82/ /90/7/ /80/ /87/10/ /80/ /78/170/ /75/ /51/ 160/ /74/ /51/ 160/ /82/ /57/9/ /62/ /77/152/ /80/ /87/ 10/ /87/ /77/ 2/ /77/ /75/ 13/ /84/ /80/5/ /87/ /89/177/ /88/ /80/178/ /87/ /58/ 176/ of15
6 Table 1. (continued) Event Number Lat Lon Z M Plane 1 f 1 CT 1 Pf1 Plane 2 f 2 CT 2 Pf /74/ /64/ 16/ /82/ /80/8/ /82/ /84/171/ /76/ /73/166/ /74/ /78/ 16/ /87/ /85/ 3/ /89/ /77/0/ /86/ /86/4/ /87/ /73/2/ /84/ /74/174/ /76/ /44/ 20/ /87/ /87/ 3/ /80/ /85/169/ /87/ /73/2/ /85/ /82/ 4/ /81/ /73/171/ /70/ /76/159/ /76/ /46/ 19/ /76/ /73/166/ /75/ /46/ 19/ /80/ /49/ 12/ /75/ /46/ 19/ /59/ /81/149/ /74/ /60/18/ /84/ /79/174/ /76/ /59/15/ /74/ /58/17/ /83/ /72/173/ /88/ /82/178/ /85/ /85/4/ /65/ /84/154/ /79/ /86/ 11/ /90/ /79/180/ /84/ /72/6/ /83/ /86/ 7/ /87/ /84/ 176/ /80/ /87/ 169/ /80/ /87/ 169/ /81/ /86/ 170/ /71/ /72/ 20/ /85/ /75/ 174/ /81/ /80/171/ /72/ /78/ 161/ /86/ /89/175/ /78/ /85/ 167/ /73/ /65/17/ /73/ /65/17/ a Z: focal depth (km); M: local magnitude; Plane1 and Plane 2: strike/dip/rake; f 1 and f 2 : expected rake angles from the best fit stress pattern; CT1 and CT2: CT values between theoretical and observed moment tensors, Pf1 and Pf2: excess pore fluid pressures (MPa). 69 MPa, respectively. These values are consistent with those determined by Häring et al. [2008], whose values of s 1, s 2, and s 3 at the same depth are 160 MPa, 122 MPa and 84 MPa. One of the reasons of the slight difference between the two studies is that, contrary to our approach, Häring et al. [2008] assumed s 2 to be perfectly vertical. The other reason is that, in the absence of any evidence for the occurrence of hydrofracking during the injection, Häring et al. [2008] constrained the value of s 3 to be slightly higher than the maximum injection pressure at the bottom of the well The Distribution of the Pore Fluid Pressure in the Basel Fluid-Injection Experiment [16] Given the absolute stress tensor and true fault planes, we estimated pore fluid pressures necessary to trigger the observed seismic events. Assuming errors of focal mechanisms to be within degrees, we chose a tight constraint of CT for events to be included in the analysis to estimate the pore fluid pressure field. This corresponds to misfit angles between the observed and expected slip vectors less than about 15 degrees. Figure 6a shows the distribution of 96 excess pore fluid pressures plotted in the 3D Mohr diagram. The maximum value of the excess pore fluid pressures is about 31 MPa (Event 46, Table 1), which is consistent with the maximum wellhead pressure. Figures 6b and 6c show the temporal and spatiotemporal evolutions of excess pore fluid pressures, respectively. Figure 6b shows that temporal change of excess pore fluid pressures before the shut-in is consistent with the history of wellhead pressures. Figure 6c shows that overpressured fluids are concentrated within 500 m from the injection point. 6of15
7 estimation errors (the standard deviation) of 2 5 MPa. Figure 8 shows that the hypocenters mainly migrated upward in the direction of S18 E. Figure 5. The depth distribution of absolute stresses. We assume that the density of rock is 2540 kg/m 3 [Häring et al., 2008]. Seismic events used in the present study occurred in the depth range from 3.77 km to 4.93 km. [17] The hypocenters of the events are localized along a near-vertical plane that strikes approximately N18 W (Figures 3a and 4c). We took this macroscopic plane of the microseismic cloud with lengths of 2800 m and 3600 m along strike and in depth, respectively, as the model region. We divided it into subsections, and distributed 165 (11 15) bi-cubic B splines (basis functions) with 200 m equally spaced local support (grid interval) to represent the distribution of the pore fluid pressure on this macroscopic plane. In this case, since the variation of pore fluid pressures within each grid interval is represented by the superposition of four local cubic functions, the theoretical limit of spatial resolution is 50 m. These values are smaller than the average resolution length expected from the distribution of seismic events (Figure 3), and so the present choice of grid intervals will not affect the results of inversion analysis. Then, determining the best estimates of the expansion coefficients (model parameters) and the covariance matrix with the Yabuki-Matsu ura inversion formula [Yabuki and Matsu ura, 1992], we estimated the pore fluid pressure field as a continuous function with estimation errors. Figure 7a shows the distribution of excess pore fluid pressures on the macroscopic plane, and Figure 7b shows the estimation errors. Figures 7a and 7b show that elevated pore fluid pressures with a peak of 17 MPa are concentrated within about 500 m of the open hole section with the 4. Discussion [18] In FMT we estimate the pore fluid pressure field from a given set of focal mechanisms of seismic events by examining the fault orientation relative to the tectonic stress field, based on an end-member model that the friction coefficient is constant for each event. This assumption is suggested by laboratory experiments of rock friction [e.g., Byerlee, 1978] and in situ stress measurements in deep boreholes [e.g., Zoback and Healy, 1992; Zoback and Townend, 2001]. The other end-member model that pore fluid pressure is hydrostatic everywhere, and the variation in focal mechanisms can attributed to variations in friction coefficients does not apply in our case because large volume of fluids were injected at high pressures, so we can assume that the effect of possible variations of friction coefficients is small in comparison. [19] In the present study we make the simplifying assumption that the tectonic stress pattern is spatiotemporally uniform in the geothermal reservoir, although pore fluid pressures should be coupled with the stress fields. The most influential stress perturbation seems to be stress changes due to seismic events triggered by increase of pore fluid pressures, although an excess pore fluid pressure will deform porous rock and cause stress perturbations. From seismic observations, the stress drop in the source region is scale independent, and it is at most 5 MPa on average [e.g., Kanamori and Anderson, 1975; Goertz-Allmann et al., 2011]. As shown in studies on Coulomb stress changes [e.g., King et al., 1994; Stein, 1999], stress changes due to a seismic event rapidly vanish with distance from the source. The stress perturbation due to seismic events is at most 10% of the absolute stress level at the depth of 4 5 km, and therefore the uniform stress pattern assumption is reasonable. In the Basel case, more than 80% of the slip vectors were consistent with those expected under the estimated uniform stress pattern, but in order to apply FMT to seismic data in a broader region, we have to consider spatial variation of stress fields. The formulation described in Section 2 is also applicable to the heterogeneous stress field, and it would be useful to combine FMT with stress inversion methods that can infer spatial variation of stress fields [e.g., Terakawa and Matsu ura, 2008, 2010]. [20] The excess pore fluid pressure for each event in the data set (Figures 6a 6c), and the distribution of pore fluid pressures (Figure 7a) are consistent with the conditions of the wellhead pressure of the injection experiment [Häring et al., 2008] and the spatiotemporal evolution of the induced microseismicity. In addition to the analysis in Section 3.2, we performed nine additional numerical experiments to examine the robustness of FMT. The conditions of these analyses are listed in Table 2. In Experiments 1 3 we estimated pore fluid pressure fields, varying the friction coefficients (Figures 9a 9c). These results show that the shape of the inferred fluid reservoir does not depend on the friction coefficient, although the absolute value of pore fluid pressures (in the upper part of the reservoir) tends to increase slightly with increasing friction coefficients. These 7of15
8 Figure 6 8of15
9 Figure 7. The distribution of pore fluid pressures formed by the fluid injection in Basel. (a) The excess pore fluid pressure field projected onto the macroscopic plane defined by the overall hypocenter distribution. The circles show hypocenters of events projected on the plane, whose excess pore fluid pressures are included in the analysis to estimate the pore fluid pressure field. The color and the size of circles are the same as in Figure 6a. (b) The distribution of estimation errors of pore fluid pressures. The size of the circles is the same as in Figure 7a. results show that FMT can reproduce reliable pore fluid pressure distributions at depth. As discussed in the next paragraph, reliability of results in the deeper parts of the reservoir is unreliable because of limited data (one data point) in this region. [21] Some of the excess pore pressures in Figures 6a 6c exceed the maximum wellhead pressure (29.6 MPa). This may be caused by observation errors of focal mechanisms and/or incorrect selection of true fault planes. In order to check the error propagation, we applied FMT to three data sets of focal mechanisms, each of which is rotated from the original focal mechanism around an arbitrary vector. In Experiments 4 and 5 we determined rotation angles using uniform random numbers. Rotation angles vary from 5 degrees to 5 degrees, and from 10 degrees to 10 degrees in Experiments 4 and 5, respectively. In Experiment 6 we determined rotation angles using a normal distribution (mean: 0 degree, standard deviation: 5 degrees). Figures 10a 10c show the distributions of excess pore fluid pressures in the three cases. In all the cases we can see two large fluid pockets with peak pore fluid pressures of about 20 MPa, indicating that the pore fluid pressure field inferred from FMT is robust. It is expected that the distribution of excess pore fluid pressures (Figure 7a) would not be affected by random errors of the data set because we statistically obtained this result using an inversion technique [Yabuki and Matsu ura, 1992]. However, the deepest fluid pocket in Figure 7a does not appear in Experiments 5 and 6 (Figures 10b 10c) because there are few data in the vicinity of the deepest parts of the stimulated rock volume. The estimation error of pore fluid pressures in the deeper part (> 4.5 km) is larger than in the region of the two large shallower pockets of high fluid pressure (Figure 7b). The existence of the deepest concentration of high fluid pressure may be less reliable. [22] We determined true fault planes by choosing the higher CT of the observed and expected moment tensors. As long as the B-axis does not coincide with that of the intermediate Figure 6. The pore fluid pressures formed by the fluid injection in Basel. (a) The excess pore fluid pressure for each event (with CT 0.965) in the 3D Mohr diagram. The horizontal and vertical axes show normal and shear stresses which are normalized by the maximum shear stresses, defined by (s 1 s 3 )/2. s 1, s 2, s 3 and P h show positions of the 3 principal stresses and the hydrostatic pressure in the horizontal axis. The color indicates the excess pore fluid pressure, and the size scales with magnitude. The line labeled H 0 indicates fault strength under hydrostatic pressure, and the other two lines (H 1,andH 2 ) under supra-hydrostatic pore fluid pressures. The inlet figure is the histogram of excess pore fluid pressures. (b) The temporal evolution of pore fluid pressures. The diamonds show excess pore pressures, and the solid line shows the history of the wellhead pressure of the injection experiment. (c) The spatiotemporal evolution of pore fluid pressures. The color circles are the same as in Figure 6a. 9of15
10 Figure 8. The spatiotemporal evolution of seismicity. Excess pore fluid pressure for each event is plotted on the excess pore fluid pressure field in Figure 7a. The color and size of the circles are the same as in Figure 6a. (a) 3 6 December 2006, (b) 7 December 2006, (c) 8 December 2006, (d) 9 24 December, and (e) Jan. to Nov principal stress, this criterion can identify the true fault plane [Gephart, 1985]. However, some incorrect selections may be caused by focal mechanism errors and/or by local stress-field heterogeneities. To determine the effect of choosing CT values that are close to each other, we used a criterion that for CT values within 0.05 (33 events), the true fault plane was chosen at random between the two possibilities (Experiment 7). The criterion of 0.05 was used because that corresponds to the uncertainty in the focal mechanisms. The distribution of pore fluid pressures (Figure 11) has a few fluid pockets. This feature is similar to that in Figure 7a, although we randomly chose true fault planes for a third of the data set. This indicates that the selection of true fault planes using CT does not result in systematic errors, and so the pore fluid pressure field inferred through FMT is unbiased. Table 2. Parameters of the Numerical Experiments a Data Set T. F. P. m Basis Function Figure The standard analysis original CT Figure 7a Experiment 1 original CT Figure 9a Experiment 2 original CT Figure 9b Experiment 3 original CT Figure 9c Experiment 4 original + random errors (5 ) CT Figure 10a Experiment 5 original + random errors (10 ) CT Figure 10b Experiment 6 original + normal random errors CT Figure 10c Experiment 7 original CT + random Figure 11 Experiment 8 original CT Figure 12a Experiment 9 original CT Figure 12b a Data set: The data set used in the present study; T. F. P: The criterion for selecting the true fault plane; m: The friction coefficient; Basis function: The number of basis functions; Figure: Figure number. 10 of 15
11 Figure 9. Dependence of the pore fluid pressure distribution on the friction coefficient: (a) m = 0.4 (Experiment 1), (b) m = 0.6 (Experiment 2), and (c) m = 1.0 (Experiment 3). The color and size of the circles are the same as in Figure 6a. Figure 10. Influences from the error propagation on the distribution of pore fluid pressures. We show the excess pore fluid pressure fields inferred from three data sets of (a) Experiment 4, (b) Experiment 5, and (c) Experiment 6. The data sets consist of focal mechanisms each of which was rotated around an arbitrary vector. Rotation angles were determined by uniform random numbers in Experiment 4 ( 5 q 5 ) and Experiment 5 ( 10 q 10 ), and by normal random numbers in Experiment 6 (mean: 0 degree, standard deviation: 5 degrees). The color and size of the circles are the same as in Figure 6a. 11 of 15
12 Figure 11. Influences from incorrect selection of true fault planes with CT on the distribution of pore fluid pressures. We show the excess pore fluid pressure field of Experiment 7 (in text). [23] The distribution of pore fluid pressures peaks (17 MPa) in a large fluid pocket south of the well, and the peak value is 56% of the applied maximum wellhead pressure. There should be a pressure gradient from borehole to the triggering location with the peak value, and so the estimate seems to be rational. But on the other hand this may be related to resolution of the present analysis. In order to check dependence of the number of the basis function (model parameters), we estimated the distributions of the pore fluid pressure using 513 (19 27) and 99 (9 11) bi-cubic B splines with 100 m and 300 m equally spaced local support, respectively (Experiments 8 and 9). The estimated pore fluid pressure fields are independent of the number of the basis function (Figures 12a 12b), although the peak value of pore fluid pressures (18 MPa in Experiment 8, 16 MPa in Experiment 9) was slightly higher with smaller size of basis function. This can be improved with additional data that would result in higher resolution of the pore pressure field. [24] Figures 6a, 6b, and 7a indicate that many events occurred at intermediate pore fluid pressures with an average of about 10 MPa 6 MPa. Such earthquake triggering pressures are significantly larger (by several orders of magnitude) than earthquake triggering pressures inferred from linear diffusion models of porous rock media [e.g., Dinske and Shapiro, 2010; Goertz-Allmann et al., 2011]. This is significant because it shows that the system is highly permeable along flow paths, allowing fast pressure transfer and substantially elevated fluid pressures within the stimulated region. The large change in local permeability at the onset of slip [e.g., Miller and Nur, 2000] should be included in Figure 12. Dependence of the number of the basis functions on the distribution of pore fluid pressures. (a) Experiment 8, and (b) Experiment 9. We used 513 and 99 basis functions in Experiments 8 and 9, respectively. 12 of 15
13 Figure 13. Roles of pore fluid pressure and tectonic stress. (a) Relationship between magnitudes of seismic events and pore fluid pressures (P f ) activating the events. The vertical axis shows the over-pressure coefficient, C, which is defined by C =(P f P h )/(s 3 P h ), where P h,ands 3 are the hydrostatic pressure and the minimum principal stress. (b) Relationship between magnitudes of seismic events and shear stress acting on fault planes. The vertical axis shows shear stress normalized by the maximum shear stress. numerical modeling efforts of EGS systems. In addition, the observation that the injected fluids mainly migrated upward in the direction of S18 E (Figure 8) indicates that future modeling efforts include the effective stress-dependence of permeability as used in previous models of fluid-rock interactions [Rice, 1992; Miller et al., 2004]. [25] The strike of the macroscopic plane of the microseismic cloud formed by the injection experiment is roughly N18 W, although only a small fraction of the available focal mechanisms have nodal planes whose strike is more or less parallel to it (Figure 4c). The orientation is very consistent with that of one of the optimally oriented planes for shear faulting with the friction coefficient of 0.85 (N13 W). This indicates that the macroscopic plane formed in response to shear faulting due to overpressured fluids. Although formation of the macroscopic plane appears to be controlled by the tectonic stress field, the evolution of this system includes many individual seismic events occurring along unfavorably oriented faults due to elevated pore fluid pressures (Figure 6a). [26] Figure 13a shows the relationship between event magnitudes and the degree of overpressure, which is defined by the overpressure coefficient, C =(P f P h )/(s 3 P h ), and Figure 13b shows the relationship between event magnitudes and shear stress normalized by the maximum shear stress at hypocentral depth. These figures show that many small events (M < 3) occurred on unfavorably oriented faults with relatively low values of shear stress, which we attribute to triggering by elevated fluid pressure. However, the 3 largest events occurred along faults with relatively high shear stress and with overpressure coefficients indicative of near hydrostatic pore pressure at the time of failure. Because of insufficient data set, an unambiguous conclusion is difficult, but we note that a similar tendency was found in the case of the 2009 L Aquila earthquake (auxiliary material), although the critical magnitude of events triggered by fluids in L Aquila is larger than that in Basel [Terakawa et al., 2010]. 1 It has also been reported for natural earthquakes that focal mechanisms of large earthquakes directly represent the regional tectonic stress field, although those of small earthquakes are diverse [e.g., England and Jackson, 1989; Amelung and King, 1997; Terakawa and Matsu ura, 2010]. [27] The implication of this observation is that small-scale triggering by reduction of fault strength due to overpressured fluids provides a mechanism for an evolving interconnected network, but whether events grow into larger ones or not appears to be controlled by tectonic stress field. The overpressured fluid reservoir may play a role as the nonasperity, slip on which would cause tectonic loading at the asperity around them [Kanamori, 1981]. How this evolving interconnected network influences larger-scale seismicity requires additional study. 5. Conclusions [28] We estimated the distribution of pore fluid pressures developed within the reservoir of the enhanced geothermal system in Basel Switzerland. The injection of high-pressure fluids induced thousands of earthquakes, and we applied focal mechanism tomography to 118 well-constrained focal mechanisms to infer the pore fluid pressure field at depth. The inferred pore fluid pressure distribution is consistent with the known fluid pressure history applied at the borehole. Our analysis demonstrates that focal mechanism tomography is a viable method for inferring pore fluid pressures, and their distribution, at depth. Our results show that triggering pressures of about 10 MPa are orders of magnitude larger than would be inferred using a simple linear diffusion model. The results demonstrate that high fluid pressures trigger earthquakes on non-optimally oriented faults, but that larger events are primarily controlled by shear failure in response to the tectonic stresses. The complex interaction between high pore pressure and shear 1 Auxiliary materials are available in the HTML. doi: / 2011JB of 15
14 failure [e.g., Miller et al., 1996] requires that future modeling efforts include evolving crack networks coupled to pressure diffusion and fluid flow. A further coupling between the evolving fluid pressure field and its poro-elastic consequences will be the focus of future studies. Appendix A [29] Determining pore fluid pressure from FMT results in a fluid pressure measurement at the discrete points of the hypocenters. In order to estimate the 3D excess fluid pressure field we apply an inversion technique developed by Yabuki and Matsu ura [1992] to the excess (above hydrostatic) fluid pressures p n (n =1,, N) of the data set. First, we represent the excess fluid pressure field p(x) by the superposition of a finite number of tri-cubic B-splines: pðxþ ¼ XM m¼1 a m F m ðxþ; ða1þ where F m are the tri-cubic B-splines, and a m (m =1,, M) are the expansion coefficients (model parameters). Then, we obtain system linear equations to be solved for the expansion coefficients a m (m =1,, M), 0 1 p 1 ðx 1 Þ p n ðx n Þ C A ¼ 0 10 F 1 ðx 1 Þ F M ðx 1 Þ B A@ F 1 ðx N Þ F M ðx N Þ a 1 a M 1 0 C B A e 1 e N 1 C A ða2þ where N is the number of data, x n is the position vector of the hypocenter of the n-th event, and e n is the observation error for p n. These observation equations can be rewritten in vector form as p ¼ Fa þ e; ða3þ where p is a N 1 data vector, a is a M 1 model parameter vector, F is a N M coefficient matrix, and e is a N 1 error vector. Assuming the errors to be Gaussian with zero mean and a covariance s 2, we obtain a stochastic model that relates the data p with the model parameter a as h i pðp a; s 2 Þ¼ð2ps 2 Þ N=2 exp ðp FaÞ T ðp FaÞ=2s 2 : ða4þ [30] In addition to observed data we can use the physical consideration that the pore fluid pressure field must be smooth to some degree because fluids are connected within a porous network. We thus impose prior constraint on the roughness of the pore fluid pressure field, which is expressed in the form of probability density function with an unknown scale factor r 2 as pða; r 2 Þ¼ð2pr 2 Þ M=2 k G k 1=2 expð a T Ga=2r 2 Þ ða5þ where k G k denotes the absolute value of the product of (nonzero) eigenvalues of G. Here G is a M M roughness matrix defined in Terakawa and Matsu ura [2008]. Combining the prior constraint pða; r 2 Þ in equation (A5) with the stochastic model pðp a; s 2 Þ in equation (A4) by Bayes rule, and introduce a new hyper-parameter a 2 (=s 2 /r 2 ) instead of r 2, we can obtain a hierarchic, highly flexible model controlled by the hyper-parameters s 2 and a 2 : pða; s 2 ; a 2 pþ ¼cð2ps 2 Þ ðmþnþ=2 ða 2 Þ M=2 k G k 1=2 exp½ sðaþ=2s 2 Š ða6þ with sðaþ ¼ðp FaÞ T ðp FaÞþa 2 a T Ga: ða7þ [31] Here, it should be noted that the hyper-parameter a 2 controls the relative weight of observed data to the prior information. Now, our problem is to find the values of a, s 2 and a 2 that maximize the posterior pdf in equation (A6) for a given data set of excess fluid pressures. To find the optimum values of hyper-parameters, we can use Akaike s Bayesian Information Criterion (ABIC) [Akaike, 1980]. The explicit expression of ABIC is given in Yabuki and Matsu ura [1992] as ABICða 2 Þ¼Nlog sða Þ M log a 2 þ log F T F þ a 2 G þ C ; ða8þ where a* denotes a maximum likelihood solution for a given a 2. Once the optimum value ^a 2, which minimizes ABIC in equation (A8), has been found, the best estimates ^a of the model parameters a and the covariance matrix C(^a) of their estimation errors are given by ^a ¼ðF T F þ ^a 2 GÞ 1 F T p Cð^aÞ ¼^s 2 ðf T F þ ^a 2 GÞ 1 : ða9þ ða10þ where ^s 2 is the optimum value of the hyper-parameter s 2, and it s obtained by: ^s 2 ¼ sð^aþ=n: ða11þ [32] The best estimate of the excess fluid pressure at an arbitrary point x o can be calculated by substituting ^a =[^a m ] in equation (A9) into equation (A1) as ^pðx o Þ¼ XM m¼1 ^a m F m ðx o Þ; and so the variance of its estimation errors is given by var½^pðx o ÞŠ ¼ XM l¼1 X M k¼1 F l ðx o ÞC lk ð^aþf k ðx o Þ; ða12þ ða13þ where C lk ð^aþ denotes the lk-element of the covariance matrix Cð^aÞ in equation (A10). We can use the square root of variances (the standard error) to measure the uncertainty of estimated excess fluid pressure fields. [33] Acknowledgments. We thank GeoPower Basel for permission to publish these results. We also thank the Editor Robert Nowack, one of the reviewers Thomas Plenefisch, and the anonymous Associate Editor and reviewer. Thoughtful comments from them were helpful to improve the manuscript. This work was supported by grant-in aid for scientific research C ( ). 14 of 15
Yusuke Mukuhira. Integration of Induced Seismicity and Geomechanics For Better Understanding of Reservoir Physics
Integration of Induced Seismicity and Geomechanics For Better Understanding of Reservoir Physics Yusuke Mukuhira Postdoctoral Fellow (JSPS research fellow) Department of Earth, Atmospheric, and Planetary
More informationFrom observations of microseismic source parameters to reservoir geomechanics
From observations of microseismic source parameters to reservoir geomechanics B. Goertz-Allmann 1, V. Gischig, C. Bachmann, and S. Wiemer AGIS workshop, 27 th November 2012 1 now at Presentation outline
More informationConstraining of Focal Mechanisms of Induced Seismicity Using Borehole Logging Information
PROCEEDINGS, 43rd Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 12-14, 2018 SGP-TR-213 Constraining of Focal Mechanisms of Induced Seismicity Using Borehole
More informationThe role of Coulomb stress changes for injection induced seismicity: The Basel enhanced geothermal system
GEOPHYSICAL RESEARCH LETTERS, VOL. 4, 72 77, doi:1.129/212gl54147, 213 The role of Coulomb stress changes for injection induced seismicity: The Basel enhanced geothermal system Flaminia Catalli, 1 Men
More informationTensor character of pore pressure/stress coupling in reservoir depletion and injection
Tensor character of pore pressure/stress coupling in reservoir depletion and injection Müller, B., Altmann, J.B., Müller, T.M., Weißhardt, A., Shapiro, S., Schilling, F.R., Heidbach, O. Geophysical Institute
More informationInterpretation of Microseismic Events of Large Magnitudes Collected at Cooper Basin, Australia and at Basel, Switzerland
Proceedings World Geothermal Congress 2010 Bali, Indonesia, 25-29 April 2010 Interpretation of Microseismic Events of Large Magnitudes Collected at Cooper Basin, Australia and at Basel, Switzerland Yusuke
More informationMigration of Shut-in Pressure and its Effect to Occurrence of the Large Events at Basel Hydraulic Stimulation
PROCEEDINGS, 41st Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 22-24, 2016 SGP-TR-209 Migration of Shut-in Pressure and its Effect to Occurrence of the
More informationBrittle Deformation. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm
Lecture 6 Brittle Deformation Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm WW Norton, unless noted otherwise Brittle deformation EarthStructure (2 nd
More informationMicroseismic monitoring of borehole fluid injections: Data modeling and inversion for hydraulic properties of rocks
GEOPHYSICS, VOL. 68, NO. 2 (MARCH-APRIL 2003); P. 685 689, 5 FIGS. 10.1190/1.1567239 Short Note Microseismic monitoring of borehole fluid injections: Data modeling and inversion for hydraulic properties
More informationMaterial is perfectly elastic until it undergoes brittle fracture when applied stress reaches σ f
Material is perfectly elastic until it undergoes brittle fracture when applied stress reaches σ f Material undergoes plastic deformation when stress exceeds yield stress σ 0 Permanent strain results from
More informationChapter 6. Conclusions. 6.1 Conclusions and perspectives
Chapter 6 Conclusions 6.1 Conclusions and perspectives In this thesis an approach is presented for the in-situ characterization of rocks in terms of the distribution of hydraulic parameters (called SBRC
More informationAssessing the solution quality of the earthquake location problem
Semesterarbeit HS 2008 Department of Earth Sciences ETH Zurich Assessing the solution quality of the earthquake location problem Sabrina Schönholzer Supervisor: Technical Supervisor: Prof. Dr. E. Kissling
More informationEstimating energy balance for hydraulic fracture stimulations: Lessons Learned from Basel
Estimating energy balance for hydraulic fracture stimulations: Lessons Learned from Basel David W. Eaton*, Department of Geoscience, University of Calgary, Calgary, Canada eatond@ucalgary.ca and Neda Boroumand,
More informationSUPPLEMENTAL INFORMATION
GSA DATA REPOSITORY 2013310 A.M. Thomas et al. MOMENT TENSOR SOLUTIONS SUPPLEMENTAL INFORMATION Earthquake records were acquired from the Northern California Earthquake Data Center. Waveforms are corrected
More informationMEASUREMENT OF HYDRAULICALLY ACTIVATED SUBSURFACE FRACTURE SYSTEM IN GEOTHERMAL RESERVOIR BY USING ACOUSTIC EMISSION MULTIPLET-CLUSTERING ANALYSIS
MEASUREMENT OF HYDRAULICALLY ACTIVATED SUBSURFACE FRACTURE SYSTEM IN GEOTHERMAL RESERVOIR BY USING ACOUSTIC EMISSION MULTIPLET-CLUSTERING ANALYSIS HIROKAZU MORIYA 1, HIROAKI NIITSUMA 1 and ROY BARIA 2
More informationAutomatic Moment Tensor Analyses, In-Situ Stress Estimation and Temporal Stress Changes at The Geysers EGS Demonstration Project
PROCEEDINGS, 42nd Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 13-15, 2017 SGP-TR-212 Automatic Moment Tensor Analyses, In-Situ Stress Estimation and
More informationMonte-Carlo Simulations of EGS Stimulation Phase with a 3-D Hybrid Model Dimitrios Karvounis and Stefan Wiemer
Monte-Carlo Simulations of EGS Stimulation Phase with a 3-D Hybrid Model Dimitrios Karvounis and Stefan Wiemer 11.03.2015 Enhanced Geothermal System in Basel In Basel: A vertical well was drilled until
More informationThe Frictional Regime
The Frictional Regime Processes in Structural Geology & Tectonics Ben van der Pluijm WW Norton+Authors, unless noted otherwise 1/25/2016 10:08 AM We Discuss The Frictional Regime Processes of Brittle Deformation
More informationON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL
1 Best Practices in Physics-based Fault Rupture Models for Seismic Hazard Assessment of Nuclear ON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL Hideo AOCHI
More informationAnalysis of Microseismic Events from a Stimulation at Basel, Switzerland
GRC Transactions, Vol. 31, 2007 Analysis of Microseismic Events from a Stimulation at Basel, Switzerland Hiroshi Asanuma 1, Yusuke Kumano 1, Akito Hotta 1, Ulrich Schanz 2, Hiroaki Niitsuma 1, and Markus
More informationMaterials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation.
Materials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation. In the process zone, stress amplitudes are poorly determined and much
More informationTectonics. Lecture 12 Earthquake Faulting GNH7/GG09/GEOL4002 EARTHQUAKE SEISMOLOGY AND EARTHQUAKE HAZARD
Tectonics Lecture 12 Earthquake Faulting Plane strain 3 Strain occurs only in a plane. In the third direction strain is zero. 1 ε 2 = 0 3 2 Assumption of plane strain for faulting e.g., reverse fault:
More informationAnalytical Stress Modeling for Mine related Microseismicity
Analytical Stress Modeling for Mine related Microseismicity Himanshu Barthwal and Mirko van der Baan University of Alberta Summary Microseismicity is recorded in an underground mine by a network of 7 boreholes
More informationFocal Mechanism Analysis of a Multi-lateral Completion in the Horn River Basin
Focal Mechanism Analysis of a Multi-lateral Completion in the Horn River Basin Paige Snelling*, Cameron Wilson, MicroSeismic Inc., Calgary, AB, Canada psnelling@microseismic.com Neil Taylor, Michael de
More informationModeling pressure response into a fractured zone of Precambrian basement to understand deep induced-earthquake hypocenters from shallow injection
Modeling pressure response into a fractured zone of Precambrian basement to understand deep induced-earthquake hypocenters from shallow injection S. Raziperchikolaee 1 and J. F. Miller 1 Abstract Analysis
More informationSHale gas. in relations to technological activity at The Geysers geothermal field, California
SHale gas Exploration and Exploitation induced Risks Temporal static stress drop variations in relations to technological activity at The Geysers geothermal field, California Monika Staszek, Beata Orlecka-Sikora,
More informationof other regional earthquakes (e.g. Zoback and Zoback, 1980). I also want to find out
4. Focal Mechanism Solutions A way to investigate source properties of the 2001 sequence is to attempt finding well-constrained focal mechanism solutions to determine if they are consistent with those
More informationCharacterization of Induced Seismicity in a Petroleum Reservoir: A Case Study
Characterization of Induced Seismicity in a Petroleum Reservoir: A Case Study Edmond Sze, M. Nafi Toksöz, and Daniel R. Burns Earth Resources Laboratory Dept. of Earth, Atmospheric and Planetary Sciences
More informationReservoir Geomechanics and Faults
Reservoir Geomechanics and Faults Dr David McNamara National University of Ireland, Galway david.d.mcnamara@nuigalway.ie @mcnamadd What is a Geological Structure? Geological structures include fractures
More informationA new methodology for estimating field horizontal stress from microseismic focal mechanisms
A new methodology for estimating field horizontal stress from microseismic focal mechanisms Alireza Agharazi 1*, Peter Duncan 1 and Michael Thornton 1 introduce a new methodology for the estimation of
More informationThis paper was prepared for presentation at the Unconventional Resources Technology Conference held in San Antonio, Texas, USA, 1-3 August 2016.
URTeC: 2461621 Determining Maximum Horizontal Stress With Microseismic Focal Mechanisms Case Studies in the Marcellus, Eagle Ford, Wolfcamp Alireza Agharazi*, MicroSeismic Inc. Copyright 2016, Unconventional
More information! EN! EU! NE! EE.! ij! NN! NU! UE! UN! UU
A-1 Appendix A. Equations for Translating Between Stress Matrices, Fault Parameters, and P-T Axes Coordinate Systems and Rotations We use the same right-handed coordinate system as Andy Michael s program,
More informationMicroseismic monitoring is a valuable
SPECIAL SECTION: M i c r o s e i s m i c moment tensors: A path to understanding frac growth ADAM BAIG and TED URBANCIC, Engineering Seismology Group Canada monitoring is a valuable tool in understanding
More information3D Finite Element Modeling of fault-slip triggering caused by porepressure
3D Finite Element Modeling of fault-slip triggering caused by porepressure changes Arsalan Sattari and David W. Eaton Department of Geoscience, University of Calgary Suary We present a 3D model using a
More informationAbstracts ESG Solutions
Abstracts ESG Solutions 2015-2016 For more information, please contact Melissa Hoy, Technical Marketing Coordinator at melissa.hoy@esgsolutions.com Visit us online at www.esgsolutions.com Abstract #1 Fracture
More informationCritical Borehole Orientations Rock Mechanics Aspects
Critical Borehole Orientations Rock Mechanics Aspects By R. BRAUN* Abstract This article discusses rock mechanics aspects of the relationship between borehole stability and borehole orientation. Two kinds
More informationInsights gained from the injection-induced seismicity in the southwestern Sichuan Basin, China
6 th Int. Symp. on In-Situ Rock Stress RS03 RS03-078 0- August 03, Sendai, Japan Insights gained from the injection-induced seismicity in the southwestern Sichuan Basin, China Xinglin Lei a * and Shengli
More informationF021 Detetection of Mechanical Failure During Hyraulic Fracturing Through Passive Seismic Microseismic Monitoring
F021 Detetection of Mechanical Failure During Hyraulic Fracturing Through Passive Seismic Microseismic Monitoring A. De La Pena* (Microseismic Inc.), L. Eisner (Microseismic Inc.), M.P. Thornton (Microseismic
More informationHigh-resolution analysis of microseismicity related to hydraulic stimulations in the Berlín Geothermal Field, El Salvador
High-resolution analysis of microseismicity related to hydraulic stimulations in the Berlín Geothermal Field, El Salvador Grzegorz Kwiatek 1, Fatih Bulut 1 Marco Bohnhoff 1, Georg Dresen 1 1) GFZ Potsdam,
More informationAftershocks are well aligned with the background stress field, contradicting the hypothesis of highly heterogeneous crustal stress
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2010jb007586, 2010 Aftershocks are well aligned with the background stress field, contradicting the hypothesis of highly heterogeneous crustal stress
More informationOverview of moment-tensor inversion of microseismic events
Overview of moment-tensor inversion of microseismic events Thomas S. Eyre 1 and Mirko van der Baan 1 Downloaded 08/11/15 to 142.244.191.52. Redistribution subject to SEG license or copyright; see Terms
More informationMagnitude, scaling, and spectral signature of tensile microseisms
Magnitude, scaling, and spectral signature of tensile microseisms David W. Eaton Department of Geoscience, University of Calgary Summary The spatial dimensions and rupture characteristics of microseismic
More informationDETAILED IMAGE OF FRACTURES ACTIVATED BY A FLUID INJECTION IN A PRODUCING INDONESIAN GEOTHERMAL FIELD
PROCEEDINGS, Thirty-Fourth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 9-11, 2009 SGP-TR-187 DETAILED IMAGE OF FRACTURES ACTIVATED BY A FLUID INJECTION
More informationAn Investigation on the Effects of Different Stress Regimes on the Magnitude Distribution of Induced Seismic Events
An Investigation on the Effects of Different Stress Regimes on the Magnitude Distribution of Induced Seismic Events Afshin Amini, Erik Eberhardt Geological Engineering, University of British Columbia,
More informationAftershock Sequences Modeled with 3-D Stress Heterogeneity and Rate-State Seismicity Equations: Implications for Crustal Stress Estimation
Pure Appl. Geophys. 167 (2010), 1067 1085 Ó 2010 The Author(s) This article is published with open access at Springerlink.com DOI 10.1007/s00024-010-0093-1 Pure and Applied Geophysics Aftershock Sequences
More informationInduced Seismicity AP Note
Note Some data used in this study belongs to Geopower Basel AG. Data and concepts published in this study may not be further used without written approval of Geopower Basel AG. SERIANEX AP 3000 AP 3000
More informationThe geomechanical significance of clay in geothermal reservoirs
The geomechanical significance of clay in geothermal reservoirs Carola Meller, Emmanuel Gaucher, Thomas Kohl KIT Institute of Applied Geosciences, Div. Geothermal Research, Adenauerring 20b, 76131 Karlsruhe,
More informationRole of lithological layering on spatial variation of natural and induced fractures in hydraulic fracture stimulation
Role of lithological layering on spatial variation of natural and induced fractures in hydraulic fracture stimulation Vincent Roche *, Department of Physics, University of Alberta, Edmonton roche@ualberta.ca
More informationSEEKING ANDERSON S FAULTING IN SEISMICITY: A CENTENNIAL CELEBRATION
SEEKING ANDERSON S FAULTING IN SEISMICITY: A CENTENNIAL CELEBRATION Bernard Célérier 1 Received 14 August 2007; accepted 3 April 2008; published 1 October 2008. [1] The reactivation of faults with near-optimal
More informationSOURCE MODELING OF RECENT LARGE INLAND CRUSTAL EARTHQUAKES IN JAPAN AND SOURCE CHARACTERIZATION FOR STRONG MOTION PREDICTION
SOURCE MODELING OF RECENT LARGE INLAND CRUSTAL EARTHQUAKES IN JAPAN AND SOURCE CHARACTERIZATION FOR STRONG MOTION PREDICTION Kimiyuki Asano 1 and Tomotaka Iwata 2 1 Assistant Professor, Disaster Prevention
More informationGeomechanical Analysis of Hydraulic Fracturing Induced Seismicity at Duvernay Field in Western Canadian Sedimentary Basin
Geomechanical Analysis of Hydraulic Fracturing Induced Seismicity at Duvernay Field in Western Canadian Sedimentary Basin Suvrat P Lele 1, Timothy Tyrrell 2, Ganeswara R Dasari 1, William A Symington 1
More informationFULL MOMENT TENSOR ANALYSIS USING FIRST MOTION DATA AT THE GEYSERS GEOTHERMAL FIELD
PROCEEDINGS, Thirty-Eighth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2013 SGP-TR-198 FULL MOMENT TENSOR ANALYSIS USING FIRST MOTION DATA AT
More informationSHALE GAS AND HYDRAULIC FRACTURING
SHALE GAS AND HYDRAULIC FRACTURING JAMES VERDON FRONTIERS IN EARTH SCIENCES (EASC M0016) WEDNESDAY 29.11.2011 Tight gas reservoirs have permeability of 0.1mD or less difficult to produce at economic rates.
More informationA mechanical model of the San Andreas fault and SAFOD Pilot Hole stress measurements
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L15S13, doi:10.1029/2004gl019521, 2004 A mechanical model of the San Andreas fault and SAFOD Pilot Hole stress measurements Jean Chéry Laboratoire Dynamique de la
More informationRotation of the Principal Stress Directions Due to Earthquake Faulting and Its Seismological Implications
Bulletin of the Seismological Society of America, Vol. 85, No. 5, pp. 1513-1517, October 1995 Rotation of the Principal Stress Directions Due to Earthquake Faulting and Its Seismological Implications by
More informationMicroseismic Geomechanical Modelling of Asymmetric Upper Montney Hydraulic Fractures
Microseismic Geomechanical Modelling of Asymmetric Upper Montney Hydraulic Fractures Drew Chorney, Byungtark Lee, Shawn Maxwell (IMaGE) Summary Geomechanical modelling is a powerful tool to quantitatively
More informationMonitoring induced microseismic events usually
SPECIAL M i c r SECTION: o s e i s m M i ci c r o s e i s m i c Beyond the dots in the box: microseismicity-constrained fracture models for reservoir simulation Leo Eisner, Sherilyn Williams-Stroud, Andrew
More informationHijiori HDR Reservoir Evaluation by Micro-Earthquake Observation
GRC Transactions, Vol. 38, 2014 Hijiori HDR Reservoir Evaluation by Micro-Earthquake Observation Hideshi Kaieda Central Research Institute of Electric Power Industry, Abiko, Chiba, Japan Keywords HDR,
More informationVerification of the asperity model using seismogenic fault materials Abstract
Verification of the asperity model using seismogenic fault materials Takehiro Hirose*, Wataru Tanikawa and Weiren Lin Kochi Institute for Core Sample Research/JAMSTEC, JAPAN * Corresponding author: hiroset@jamstec.go.jp
More informationModel Inversion for Induced Seismicity
Model Inversion for Induced Seismicity David Castiñeira Research Associate Department of Civil and Environmental Engineering In collaboration with Ruben Juanes (MIT) and Birendra Jha (USC) May 30th, 2017
More informationKinematic inversion of pre-existing faults by wastewater injection-related induced seismicity: the Val d Agri oil field case study (Italy)
Kinematic inversion of pre-existing faults by wastewater injection-related induced seismicity: the Val d Agri oil field case study (Italy) Buttinelli M., Improta L., Bagh S., Chiarabba C. 1/10 The Val
More informationCHANGES OF COULOMB FAILURE STRESS DUE TO DISLOCATIONS DURING STIMULATION OF GPK2 AT SOULTZ-SOUS-FORÊTS
PROCEEDINGS, Thirty-Seventh Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 30 - February 1, 2012 SGP-TR-194 CHANGES OF COULOMB FAILURE STRESS DUE TO DISLOCATIONS
More informationProduction-induced stress change in and above a reservoir pierced by two salt domes: A geomechanical model and its applications
Production-induced stress change in and above a reservoir pierced by two salt domes: A geomechanical model and its applications Peter Schutjens, Jeroen Snippe, Hassan Mahani, Jane Turner, Joel Ita and
More informationPossibility of reservoir induced seismicity around three gorges dam on Yangtze river
Int. J. Rock Mech. & Min. Sci. Vol. 34, No. 3-4, 1997 To cite this paper: Int. J. RockMech. &Min. Sci. 34:34, Paper No. 076 Possibility of reservoir induced seismicity around three gorges dam on Yangtze
More informationHydraulic Characteristics of the Basel 1 Enhanced Geothermal System
GRC Transactions, Vol. 33, 9 Hydraulic Characteristics of the Basel 1 Enhanced Geothermal System Florentin Ladner and Markus O. Häring Geothermal Explorers International Ltd, Switzerland Keywords Geothermal,
More informationA review of friction laws and their application for simulation of microseismicity prior to hydraulic fracturing
A review of friction laws and their application for simulation of microseismicity prior to hydraulic fracturing Jiyang Ye, Mirko Van Der Baan (Email: jiyang1@ualberta.ca, Mirko.VanderBaan@ualberta.ca)
More informationJ.V. Herwanger* (Ikon Science), A. Bottrill (Ikon Science) & P. Popov (Ikon Science)
29829. One 4D geomechanical model and its many applications J.V. Herwanger* (Ikon Science), A. Bottrill (Ikon Science) & P. Popov (Ikon Science) Main objectives (i) Field case study demonstrating application
More informationObservation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story
Observation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story Petr Kolinsky 1, Leo Eisner 1, Vladimir Grechka 2, Dana Jurick 3, Peter Duncan 1 Summary Shear
More informationStress Orientations Obtained from Earthquake Focal Mechanisms: What Are Appropriate Uncertainty Estimates?
Bulletin of the Seismological Society of America, 91, 2, pp. 250 262, April 2001 Stress Orientations Obtained from Earthquake Focal Mechanisms: What Are Appropriate Uncertainty Estimates? by Jeanne L.
More informationPassive seismic monitoring in unconventional oil and gas
Passive seismic monitoring in unconventional oil and gas Michael Kendall, James Verdon, Alan Baird, Anna Stork and Philip Usher Bristol University Microseismicity Projects (BUMPS) Microseismicity and
More informationHDR PROJECT SOULTZ: HYDRAULIC AND SEISMIC OBSERVATIONS DURING STIMULATION OF THE 3 DEEP WELLS BY MASSIVE WATER INJECTIONS
PROCEEDINGS, Thirty-Second Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 22-24, 27 SGP-TR-183 HDR PROJECT SOULTZ: HYDRAULIC AND SEISMIC OBSERVATIONS DURING
More informationThis paper was prepared for presentation at the Unconventional Resources Technology Conference held in Denver, Colorado, USA, August 2014.
URTeC 1922263 Utilizing Ant-tracking to Identify Slowly Slipping Faults in the Barnett Shale Noha Sameh Farghal* and Mark D. Zoback, Stanford University, Stanford, CA, USA Copyright 2014, Unconventional
More informationExercise: concepts from chapter 6
Reading: Fundamentals of Structural Geology, Chapter 6 1) The definition of the traction vector (6.7) relies upon the approximation of rock as a continuum, so the ratio of resultant force to surface area
More informationTectonic Seismogenic Index of Geothermal Reservoirs
Tectonic Seismogenic Index of Geothermal Reservoirs C. Dinske 1, F. Wenzel 2 and S.A. Shapiro 1 1 Freie Universität Berlin 2 KIT Karlsruhe November 27, 2012 Introduction M max reservoir location -2.0 Barnett
More informationEarth Tides Can Trigger Shallow Thrust Fault Earthquakes
Earth Tides Can Trigger Shallow Thrust Fault Earthquakes Elizabeth S. Cochran, 1 * John E. Vidale, 1 Sachiko Tanaka 2 1 Department of Earth and Space Sciences and Institute of Geophysics and Planetary
More informationLocation uncertainty for a microearhquake cluster
Analysis of location uncertainty for a microearhquake cluster: A case study Gabriela Melo, Alison Malcolm, Oleg Poliannikov, and Michael Fehler Earth Resources Laboratory - Earth, Atmospheric, and Planetary
More informationRELOCATION OF THE MACHAZE AND LACERDA EARTHQUAKES IN MOZAMBIQUE AND THE RUPTURE PROCESS OF THE 2006 Mw7.0 MACHAZE EARTHQUAKE
RELOCATION OF THE MACHAZE AND LACERDA EARTHQUAKES IN MOZAMBIQUE AND THE RUPTURE PROCESS OF THE 2006 Mw7.0 MACHAZE EARTHQUAKE Paulino C. FEITIO* Supervisors: Nobuo HURUKAWA** MEE07165 Toshiaki YOKOI** ABSTRACT
More informationFinite element modelling of fault stress triggering due to hydraulic fracturing
Finite element modelling of fault stress triggering due to hydraulic fracturing Arsalan, Sattari and David, Eaton University of Calgary, Geoscience Department Summary In this study we aim to model fault
More informationCoupling between deformation and fluid flow: impacts on ore genesis in fracture-controlled hydrothermal systems
Coupling between deformation and fluid flow: impacts on ore genesis in fracture-controlled hydrothermal systems Stephen F Cox Research School of Earth Sciences The Australian National University INTRODUCTION
More informationPREDICTIVE MODELING OF INDUCED SEISMICITY: NUMERICAL APPROACHES, APPLICATIONS, AND CHALLENGES
PREDICTIVE MODELING OF INDUCED SEISMICITY: NUMERICAL APPROACHES, APPLICATIONS, AND CHALLENGES Mark McClure Assistant Professor Petroleum and Geosystems Engineering The University of Texas at Austin Overview
More informationPreliminary slip model of M9 Tohoku earthquake from strongmotion stations in Japan - an extreme application of ISOLA code.
Preliminary slip model of M9 Tohoku earthquake from strongmotion stations in Japan - an extreme application of ISOLA code. J. Zahradnik 1), F. Gallovic 1), E. Sokos 2) G-A. Tselentis 2) 1) Charles University
More informationJihoon Kim, George J. Moridis, John Edmiston, Evan S. Um, Ernest Majer. Earth Sciences Division, Lawrence Berkeley National Laboratory 24 Mar.
TOUGH+ROCMECH for the Analysis of coupled Flow, Thermal, Geomechanical and Geophysical Processes Code Description and Applications to Tight/Shale Gas Problems Jihoon Kim, George J. Moridis, John Edmiston,
More informationMicroearthquake Focal Mechanisms
Microearthquake Focal Mechanisms A Tool for Monitoring Geothermal Systems By Bruce R. Julian (U. S. Geological Survey - Menlo Park, CA) and Gillian R. Foulger (University of Durham - Durham, United Kingdom)
More information2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
STRUCTURE OF THE KOREAN PENINSULA FROM WAVEFORM TRAVEL-TIME ANALYSIS Roland Gritto 1, Jacob E. Siegel 1, and Winston W. Chan 2 Array Information Technology 1 and Harris Corporation 2 Sponsored by Air Force
More informationEarthquake Mechanisms and Stress Field
Earthquake Mechanisms and Stress Field Václav Vavryčuk* Institute of Geophysics, Czech Academy of Sciences, Praha 4, Czech Republic Synonyms Dynamics and mechanics of faulting; Earthquake source observations;
More informationPEAT SEISMOLOGY Lecture 12: Earthquake source mechanisms and radiation patterns II
PEAT8002 - SEISMOLOGY Lecture 12: Earthquake source mechanisms and radiation patterns II Nick Rawlinson Research School of Earth Sciences Australian National University Waveform modelling P-wave first-motions
More informationDetermination of Maximum Horizontal Field Stress from Microseismic Focal Mechanisms A Deterministic Approach
ARMA 16-691 Determination of Maximum Horizontal Field Stress from Microseismic Focal Mechanisms A Deterministic Approach Agharazi, A MicroSeismic Inc., Houston, Texas, USA Copyright 2016 ARMA, American
More information** New Energy and Industrial Technology Development Organization (NEDO) Keisoku Co., Ltd. (CK)
et Micro-Earthquake Data Processing and Analysis System (MEPAS), a Software for Geothermal Applications Shin-ichi Shigeo * Takashi Japan Metals and Chemicals Co., Ltd. ** New Energy and Industrial Technology
More informationAddressing the risks of induced seismicity in sub-surface energy operations
Addressing the risks of induced seismicity in sub-surface energy operations Richard Porter a, Alberto Striolo a, Haroun Mahgerefteh a, Joanna Faure Walker b a Department of Chemical Engineering, University
More informationProperties of large ruptures and the dynamical influence of fluids on earthquakes and faulting
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 107, NO. B9, 2182, doi:10.1029/2000jb000032, 2002 Properties of large ruptures and the dynamical influence of fluids on earthquakes and faulting S. A. Miller 1 Geology
More informationDevelopments on Microseismic Monitoring and Risk Assessment of Large-scale CO 2 Storage
Developments on Microseismic Monitoring and Risk Assessment of Large-scale CO 2 Storage Bettina Goertz-Allmann NORSAR CCS Technical Workshop, RITE, Tokyo, January 16, 2019 Outline Introduction: Induced
More informationFault Reactivation Predictions: Why Getting the In-situ Stresses Right Matters
Fault Reactivation Predictions: Why Getting the In-situ Stresses Right Matters Pat McLellan, M.Sc., P.Eng. Principal Consultant McLellan Energy Advisors Inc. Calgary, Alberta May 8, 2015 Presented at the
More informationSummary so far. Geological structures Earthquakes and their mechanisms Continuous versus block-like behavior Link with dynamics?
Summary so far Geodetic measurements velocities velocity gradient tensor (spatial derivatives of velocity) Velocity gradient tensor = strain rate (sym.) + rotation rate (antisym.) Strain rate tensor can
More informationDynamic source inversion for physical parameters controlling the 2017 Lesvos earthquake
Dynamic source inversion for physical parameters controlling the 2017 Lesvos earthquake F. Kostka 1, E. Sokos 2, J. Zahradník 1, F. Gallovič 1 1 Charles University in Prague, Faculty of Mathematics and
More informationGain information on the development of Engineered Geothermal Systems (EGS) and enhance our understanding of long-term reservoir behaviour.
1 GEOTHERMAL ICL was involved in the first commercial Enhanced Geothermal System (EGS) project in the EU and for over 15 years has undertaken international research programs for the monitoring of rock
More informationEarthquake patterns in the Flinders Ranges - Temporary network , preliminary results
Earthquake patterns in the Flinders Ranges - Temporary network 2003-2006, preliminary results Objectives David Love 1, Phil Cummins 2, Natalie Balfour 3 1 Primary Industries and Resources South Australia
More informationINFLUENCE OF LOCAL PERTURBATION ON REGIONAL STRESS AND ITS IMPACT ON THE DESIGN OF MAJOR UNDERGROUND STRUCTURE IN HYDROELECTRIC PROJECT
INFLUENCE OF LOCAL PERTURBATION ON REGIONAL STRESS AND ITS IMPACT ON THE DESIGN OF MAJOR UNDERGROUND STRUCTURE IN HYDROELECTRIC PROJECT *D.S. Subrahmanyam National Institute of Rock Mechanics, Bangalore
More informationSeismological monitoring of the GRT1 hydraulic stimulation (Rittershoffen, Alsace, France)
Seismological monitoring of the GRT1 hydraulic stimulation (Rittershoffen, Alsace, France) Mohamed Boubacar 1, Olivier Lengliné 2 1. University of Quebec, INRS 2. University of Strasbourg, EOST Summary
More informationTheory - SBRC fundamentals
Chapter 2 Theory - SBRC fundamentals The approach proposed to estimate hydraulic properties of rocks using microseismicity is called Seismicity Based Reservoir Characterization (SBRC). It uses a spatio-temporal
More informationGrimsel In-situ Stimulation and Circulation experiment: First results SCCER Annual meeting , Birmensdorf, Switzerland
Grimsel In-situ Stimulation and Circulation experiment: First results SCCER Annual meeting 14 15.09.2017, Birmensdorf, Switzerland J. Doetsch, F. Amann, V. Gischig, M. Jalali, H. Krietsch, B. Valley, C.
More information