PUBLICATIONS. Geophysical Research Letters. Linking microearthquakes to fracture permeability change: The role of surface roughness
|
|
- Susan Owens
- 5 years ago
- Views:
Transcription
1 PUBLICATIONS RESEARCH LETTER Key Points: A plausible scale dependence for channeling flow through joints and faults is formulated A quantitative relation between the moment magnitude of MEQs and the fracture permeability change is obtained Potential of in situ MEQ data that allow inverse mapping of permeability change is suggested Correspondence to: T. Ishibashi, takuya.ishibashi@aist.go.jp Citation: Ishibashi, T., N. Watanabe, H. Asanuma, and N. Tsuchiya (2016), Linking microearthquakes to fracture permeability change: The role of surface roughness, Geophys. Res. Lett., 43, , doi: / 2016GL Received 20 APR 2016 Accepted 6 JUL 2016 Accepted article online 11 JUL 2016 Published online 25 JUL American Geophysical Union. All Rights Reserved. Linking microearthquakes to fracture permeability change: The role of surface roughness Takuya Ishibashi 1, Noriaki Watanabe 2, Hiroshi Asanuma 1, and Noriyoshi Tsuchiya 2 1 Fukushima Renewable Energy Institute, National Institute of Advanced Industrial Science and Technology, Koriyama, Japan, 2 Graduate School of Environmental Studies, Tohoku University, Sendai, Japan Abstract Despite its importance, the relation between microearthquakes (MEQs) and changes in hydraulic properties during hydraulic stimulation of a fractured reservoir has rarely been explored, and it is still not well understood. To investigate this relation, we first formulate a plausible scale dependence, where fracture length and shear displacement are variables, for channeling flow through heterogeneous aperture distributions for joints and faults. By combining this formulation with the concept of the seismic moment, we derive quantitative relations between the moment magnitude (M w ) of MEQs and the fracture permeability change in the directions orthogonal to (k fault, /k joint ) and parallel to (k fault, /k joint ) the shear displacement, in the form k fault; =k joint ¼ 116:410 0:46Mw and k fault;===k joint ¼ 13:110 0:46Mw. Despite the simplicity of the derivation, these relations have the potential to explain the results of field experiments on hydraulic stimulation, such as the enhanced geothermal systems at Soultz-sous-Fôret and Basel. 1. Introduction Traditionally, fractured reservoirs such as geothermal or hydrocarbon reservoirs are stimulated hydraulically for the purpose of improving or maintaining their permeability [Evans et al., 2005; Häring et al., 2008]. During this stimulation, a massive amount of pressurized fluid is injected into the reservoir, and as a result, preexistent fractures undergo slip which induces microearthquakes (MEQs) [Ellsworth, 2013; Guglielmi et al., 2015]. In situ data on such MEQs are generally recorded, since they can be good indicators of active processes within the reservoir [Majer et al., 2007]. The mechanisms involved in (hydro-)shear slip of rock fractures have been widely investigated and such slip is generally explained in terms of the Mohr- Coulomb (MC) failure criterion. Specifically, changes in mechanical properties, e.g., based on the rate and state friction (RSF) model, have been actively studied [Marone, 1998], and by coupling the RSF concept with the MC failure criterion, McClure and Horne [2012] have recently been successful in reproducing the diversity of timings and magnitudes of MEQs induced by hydraulic stimulation of a fractured system. This has led to a better understanding of the link between MEQs and the mechanical properties of rock fractures. In contrast, the relation between MEQs and changes in hydraulic properties has rarely been explored and is still poorly understood. Recently, however, the link between MEQs and the change in permeability in fractures and faults has become one of the hottest topics in geo-engineering [McClure and Horne, 2012] and seismology [Shelly et al., 2013], since understanding this relation can potentially open the door to novel approaches for coupling hydraulic and seismic effects. Since the permeability of a fractured system has been shown to change during or after hydraulic stimulation [Evans et al., 2005; Häring et al., 2008], it is natural to explore a link with MEQs in such systems. In the present study, we focus specifically on the role of the surface topography (roughness) of rock fractures for the hydraulic properties of the fractured system. For a fracture consisting of two opposing rock faces, such roughness can produce asperities that give rise to self-propping of the fracture during a MEQ, and this is highly likely to change the hydraulic properties of the fracture or reservoir [Preisig et al., 2015]. Here the roughness of natural fractures and faults can be characterized by a self-affine fractal nature regardless of their scale (length) [Renard et al., 2013]. Thus, the distribution of apertures formed by these rough surfaces remains heterogeneous regardless of shear displacement [Ameli et al., 2013; Watanabe et al., 2009]. Considering such characteristics of rock fractures, we can determine the significance of geometrical properties of aperture distribution and resulting fluid flow through rock fractures with natural (field) scales. Ishibashi et al. [2015] recently proposed a prediction method for aperture distributions of beyond-laboratory-scale rock ISHIBASHI ET AL. LINKING MEQS TO PERMEABILITY CHANGE 7486
2 fractures at in situ reservoir conditions. However, the relation between MEQs and changes in fracture permeability has still not been clarified. In the present study, our main focus is achieving a quantitative relation between the magnitude of a MEQ and the change in permeability during hydraulic stimulation of a fractured system. Based on the results of numerical modeling, we derive the scale dependences of mechanical and hydraulic properties of subsurface rock fractures (joint/fault), where the shear displacement of the fracture is treated as a variable parameter rather than having a fixed value. Consequently, the dependence of mechanical and hydraulic properties of rock fractures on scale (length) and shear displacement can be successively determined. By introducing the concept of the seismic moment, which relates the magnitude of a MEQ to the fracture scale, we then obtain the relation between the MEQ magnitude and the permeability change. 2. Methods Based on the method developed by Ishibashi et al. [2015], heterogeneous aperture distributions for effective confining stresses of up to ~100 MPa (corresponding to depths of ~6 7 km in a continental crust) are simulated for rock fractures with various combinations of length l (m) and shear displacement δ (m). In their method, a pair of synthetic fractal surfaces were placed together, so that the fractures had a scaleindependent contact area (see Ishibashi et al. [2015] for details). The contact area of a fracture having shear displacement was 42%, whereas that of fracture having no shear displacement was 59%. In our study, shear displacements (δ) of more than ~0.4 mm were considered. This limitation was because of the change from joint to fault fracture characteristics for δ between 0 mm and ~0.4 mm [Chen et al., 2000]. Following Ishibashi et al. [2015], for the rough fracture surfaces (hanging wall and footwall), a fractal dimension, a standard deviation for a referential length of 200 mm, and a mismatch length were set as 2.3 mm, 1.3 mm, and 0.7 mm. The proportion of the power spectral density (PSD) of the initial aperture with a single contact point to the PSD of the footwall surface geometry, R(f), was expressed as Rf ðþ¼ef 6: ðln fþ 3 2: ðln fþ 2 þ2: ðln fþþ5: g, where f is the spatial frequency. These values and relationships are essential for numerical modeling of fracture surface topographies and were computed based on fractures in Inada granite (Ibaraki, Japan). Using the aforementioned parameter values, we numerically paired surface topographies of square fractures with fracture lengths of 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, and 0.6 m, using a 250 μm square-shaped grid. Heterogeneous aperture structures of fractures with shear displacement (i.e., faults) and without shear displacement (i.e., joint) were then modeled over various scales. To imitate natural faults, δ/l was kept constant, which corresponds to a linear increase in shear displacement with fracture length [Kanamori and Anderson, 1975]. However, since δ/l inherently depends on both rock properties and tectonic environment, values of , , , and were considered. Thus, a total of 35 aperture distributions were considered in this study. For each distribution, the geometrical average and standard deviation for nonzero aperture sizes were first calculated. Unidirectional fluid flow was then simulated by solving the Reynolds equation [Brown, 1987; Brush and Thomson, 2003]. In the simulations, flow both orthogonal and parallel to the imposed shear displacement direction was considered, since it has been reported that there is a large degree of hydraulic anisotropy associated with hydraulic fracturing in fractured reservoirs [Nemoto et al., 2009]. Based on the fluid flow simulations, changes in the fracture permeability and flow area which corresponds to the region of the channeling flow paths within a fracture were evaluated for different fracture lengths and shear displacements. Finally, the scale dependence of fracture flow was combined with equations defining the scale of earthquakes [Aki and Richards, 2002; Hanks and Kanamori, 1979], and the relation between MEQs and the change in the fracture permeability was investigated. 3. Results and Discussion 3.1. Scale Dependence of Channeling Flow Through Rock Fracture Regardless of the fracture length or the shear displacement, rock fractures contain a large number of contacting points, and the distribution of nonzero apertures is skewed and has a long tail, i.e., a lognormal appearing ISHIBASHI ET AL. LINKING MEQS TO PERMEABILITY CHANGE 7487
3 Figure 1. Representative results for channeling flow within heterogeneous aperture distribution for model fractures (a) with no shear displacement and (b) with constant shear displacement of δ/l = distribution. Hence, preferential flow paths are always formed through the heterogeneous aperture distribution (i.e., channeling flow). Figure 1 shows typical aperture maps (gray scale) and flow rate maps (color scale) for (a) a joint (δ/l=0) and (b) a fault (δ/l= ); for both fractures, channeling flow can clearly be seen. These results are completely different than those for conventional parallel-plate fractures. Note that the flow rates are normalized, normalized flow rates of less than 0.01 are colorless, and the shear displacement is in the direction orthogonal to the fluid flow. As shown in Figure 2a, the average aperture size of joints/faults increases with fracture length. In addition, the average aperture size for faults with a constant length increases with shear displacement (i.e., δ/l), as a result of the shear dilation effect [Chen et al., 2000; Watanabe et al., 2009]. In the case of a joint, the average aperture size can be approximated by e m;joint ¼ 1:310 1 l 0:10 ; (1) where e m, joint and l are the average aperture size for the joint (mm) and the fracture length (m), respectively. For a fault with a constant δ/l, the average aperture size can be approximated by e m;fault ¼ α 1 l β 1 ; (2) Figure 2. Predicted fracture-length dependence of (a) average aperture size and (b) standard deviation, for fractures with no shear displacement (joints) and fractures with constant shear displacement to fracture length ratios (faults). Curves from equations (1) and (2) are also shown. ISHIBASHI ET AL. LINKING MEQS TO PERMEABILITY CHANGE 7488
4 Figure 3. Predicted fracture-length dependence of (a) permeability in direction orthogonal to shear displacement, k ; (b) permeability in direction parallel to shear displacement, k // ; (c) permeability anisotropy ratio k /k ; and (d) flow area for joints and faults. Curves from equations (4) and (6) are also shown. where e m, fault is the average aperture size for the fault (mm) and α 1 and β 1 are a preexponential factor and an exponent, respectively. Curves from equation (2) are also shown in Figure 2a, and they yield α 1 = (δ/l) 0.59 and β 1 = Equation (2) can then be rearranged to give e m;fault ¼ 1:310 δ 0:59 l 0:71 : (3) l Equations (1) and (3) indicate that the scale dependence of e m, joint is weak, whereas that of e m, fault is strong. This discrepancy is due to differences in the degree of matedness, which is defined by how well matched the opposing fracture surfaces are [Olsson and Barton, 2001]. As shown in Figure 2b, the standard deviation is approximately 3, and its dependence on l or δ/l is ambiguous. The predicted fracture-length dependence of hydraulic properties is shown in Figure 3. Figure 3a shows that fracture permeability (k ) in the direction orthogonal to the shear displacement increases with fracture length for both joints and faults. For joints, the permeability can be approximated by k joint ¼ 9: l 0:16 ; (4) where k joint is the joint permeability (m 2 ). For faults with a constant δ/l, the permeability can be approximated by k fault; ¼ α 2 l β 2 ; (5) where k fault, is the fault permeability (m 2 ) and α 2 and β 2 are a preexponential factor and an exponent, respectively. Curves from equation (5) yield α 2 = (δ/l) 1.18 and of β 2 = Equation (5) can then be rearranged to give k fault; ¼ 2:310 6 δ 1:18 l 1:08 : (6) l Figure 3b shows fracture permeability (k // ) in the direction parallel to the shear displacement. The line and the shaded region represent the k values shown in Figure 3a. Figure 3c shows the permeability anisotropy ratio, ISHIBASHI ET AL. LINKING MEQS TO PERMEABILITY CHANGE 7489
5 Figure 4. (a) Model used to link microearthquakes to permeability change during hydraulic stimulation, and estimated relations between moment magnitude and (b) fracture length and shear displacement, (c) average aperture size, and (d) fracture permeability. which is defined as k /k //. For the joint, this ratio is about 1 regardless of the fracture length, indicating that the permeability anisotropy is low. For faults, the ratio varies from 2 to 40, indicating a high degree of permeability anisotropy. Since the permeability anisotropy ratio is poorly correlated with fracture length or shear displacement (Figure 3c), we evaluate the arithmetic mean value of permeability ratio as 8.9. The fault permeability (k fault, // ) in the direction parallel to the shear displacement is then arranged as k fault;== ¼ k fault; 8:9 ¼ δ 1:18 2:610 7 l 1:08 : (7) l The permeability anisotropy may be caused by contacting asperities grouped along the direction orthogonal to the shear displacement that partially block the fluid flow [Nemoto et al., 2009]. As shown in Figure 3d, the flow area is limited to about 5 20% of the fracture plane for fluid flow in the direction orthogonal to the introduced shear displacement. This is typical of fluid flow within rock fractures, matching the results of field investigations [Tsang and Neretnieks, 1998] Relation Between Microearthquakes and Change in Fracture Permeability Having established the dependence of the average aperture size and fracture permeability on the fracture length and shear displacement, we introduce the concept of seismic moment. This concept relates the scale of an earthquake to the fracture length, allowing a quantitative relationship between the change in permeability and the magnitude of MEQs. As shown in Figure 4a, we consider a simple configuration of hydroshear failure for the MEQs generated during stimulation of fractured reservoirs. A shear displacement of δ (m) is uniformly applied to a square fracture plane with a fracture length of l (m), although in reality the shear displacement would have a nonuniform distribution [Preisig et al., 2015]. For such a MEQ, the seismic moment is defined as [Aki and Richards, 2002] M 0 ¼ μ 0 Aδ; (8) where M 0 is the seismic moment (N m), μ 0 is the shear modulus for the host rock, and A is the surface area of the fracture (m 2 ). This seismic moment can be converted into the equivalent moment magnitude as log M 0 ¼ 1:5M w þ 9:1; (9) ISHIBASHI ET AL. LINKING MEQS TO PERMEABILITY CHANGE 7490
6 where M w is the moment magnitude [Hanks and Kanamori, 1979]. As long as the stress drop has a unique value, A and δ can be estimated as functions of M w using empirical scaling relations [Zoback and Gorelick, 2012] and assuming a static stress drop. According to Leonard [2010], for strike-slip earthquakes of relatively small magnitude, the seismic moment is related to the fracture surface area as log M 0 ¼ 1:5log A þ 6:09: (10) Because the fracture surface area is equivalent to the square of the fracture length (i.e., A = l 2 ), equation (10) can be rearranged to yield The shear displacement is given by l ¼ 10 0:5Mwþ1:00 : (11) log δ ¼ 0:5log A 4:43 ¼ log l 4:43; (12) which indicates that δ/l is constant (δ/l= ). The small value of δ/l is due to the fact that a fault experiences a single slip when using the scaling for fracture surface area (equations (10)) and shear displacement (equation (12)). Since these scaling of strike-slip earthquakes (i.e., equations (10) and (12)) are virtually identical to those of dip-slip earthquakes, the following discussions are generally acceptable for both strikeslip and dip-slip earthquakes. Thus, M w, l, and δ can be quantitatively correlated. For a fracture that has not experienced shear failure (i.e., a joint prior to hydraulic stimulation), we assume zero shear displacement. This allows the initial values for average aperture size (e m, joint ) and fracture permeability (k joint ) for a fracture with a given length to be estimated using equations (1) and (4). In contrast, for a fracture that has experienced a MEQ (i.e., a fault after hydraulic stimulation), the average aperture size (e m, fault ), the permeability (k fault, ) in the direction orthogonal to the shear displacement, and the permeability (k fault,// ) in the direction parallel to the shear displacement can be estimated using equations (3), (6), and (7), respectively. Here the shear displacement is limited to values larger than ~0.4 mm (M w > 0 from equations (11) and (12)). In this manner, the quantitative change in average aperture size and fracture permeability can be expressed in terms of M w for a MEQ. By combining equations (1), (4), (11), and (12), the change in the average aperture size is e m;fault =e m;joint ¼ 1:010 0:35Mw : (13) Similarly, by combining equations (3), (6), (11), and (12), the corresponding change in the fracture permeability in the orthogonal direction to the shear displacement can be formulated as k fault; =k joint ¼ 116:410 0:46Mw : (14) Finally, by combining equations (3), (7), (11), and (12), the change in the fracture permeability in the parallel direction to the shear displacement is k fault;===k joint ¼ 13:110 0:46Mw : (15) The generality of these expressions is clearly limited, because they can be influenced by rock type, the roughness of fracture surface, fault gouge formation, and the absolute value of the static stress drop. Next, l and δ are estimated for a MEQ with a M w of 0 2, and the average aperture size and the fracture permeability are calculated before and after the earthquake. Figure 4b shows the relation between l and M w and between δ and M w. Figures 4c and 4d, respectively, show the change in the average aperture size and the fracture permeability due to MEQs with a specific M w. MEQs with M w 0 cause an increase in the average aperture size due to shear dilation. In contrast, the impact of shear dilation is expected to be small or absent for MEQs of M w < 0, and this is consistent with previously reported experimental results [Chen et al., 2000; Javadi et al., 2014]. As a result, the increase in permeability due to a MEQ with M w < 0 is expected to be very small or zero, so that during stimulation of a fractured reservoir, such increases would only be expected for zones experiencing seismic events with M w 0. Although the mechanisms of fracture shear are much more complicated than is assumed here (e.g., a heterogeneous slip distribution along the fracture), we feel that the proposed approach is valuable as a first step toward linking MEQs to fracture permeability change. ISHIBASHI ET AL. LINKING MEQS TO PERMEABILITY CHANGE 7491
7 3.3. Comparison With Results of Field Experiments We next compare the permeability increases predicted using the proposed method with the results of actual field experiments. According to Evans et al. [2005], the hydraulic stimulation for creating an enhanced geothermal system (EGS) at Soultz-sous-Fôret, France, caused several MEQ events with local magnitudes of , and as a result, the reservoir transmissivity increased from ~ to ~ m 2 (approximately 200 times). Similar transmissivity gains (400 times) were reported at the EGS project in Basel, Switzerland, where the local magnitude of observed MEQs was 2.5 [Häring et al., 2008]. In comparison, for a MEQ with M w = 2.5, equations (14) and (15), respectively, predict permeability gains of 1640 and 185 times. Since these equations represent extreme cases with regard to the slip direction, the actual permeability increase would be expected to be between these values. This is consistent with the results of the field experiments, and thus supports the validity of the approach proposed here. In future work, we will attempt to refine this method, for example, by considering uncertainties in the stress drop [Zoback and Gorelick, 2012] or the difference between the characteristics of a single rock fracture and those of a discrete rock fracture network. The results of the present study suggest that the relationship between earthquake dynamics and crustal fluids is likely influenced by the in situ characteristics of the Earth s crust, including the stress field and the fault geometry (size and morphologic characteristics), as well as by the characteristics of the crustal fluid. This relationship can be simply formulated and constrained by considering both scale-independent parameters (ratio of asperity contact and flow area) and scale-dependent parameters (average aperture, fracture permeability, and shear displacement). This raises the possibility of new interpretations of subsurface hydraulic properties based on seismic data. 4. Conclusions We explore the linkage between MEQs and fracture permeability change during hydraulic stimulation of a fracture reservoir. Plausible scale dependencies were systematically formulated for channeling flow within heterogeneous aperture structures for the case of joints (fractures having no shear displacement) and faults (fractures having shear displacement), for different fracture lengths and shear displacements. These formulations were then combined with the concept of the seismic moment, which connects the scale of the earthquake with the fracture length. Consequently, expressions were obtained for the relation between the increase in fracture permeability and the moment magnitude of MEQs during hydraulic stimulation, in the form k fault; =k joint ¼ 116:410 0:46Mw and k fault;===k joint ¼ 13:110 0:46Mw. This may allow rough inverse mapping of changes in fracture permeability using in situ MEQ data and possibly lead to novel views into fluid migration phenomenon within the Earth s crust. This is also intimately connected to the development of fractured reservoirs of hydrocarbon/geothermal or could lead to breakthrough for earthquake mechanisms. Acknowledgments The authors thank Steven Ingebritsen for his thoughtful comments that significantly improved the manuscript. This work also benefitted from discussions with Derek Elsworth, Yi Fang, Atsushi Okamoto, and Nobuo Hirano. The present study was supported in part by JSPS through Postdoctoral Fellowships for Research Abroad, (to T.I.), a Grant-in-Aid for Specially Promoted Research, (to N.T.), and by METI, Japan through International Research Program for Innovative Energy Technology (to H.A.). The data to reproduce the work are available from the corresponding author: Takuya Ishibashi (takuya.ishibashi@aist.go.jp). References Aki, K., and P. G. Richards (2002), Quantitative Seismology, 2nd ed., 700 pp., Univ. Sci. Books, Sausalito, Calif. Ameli, P., J. E. Elkhoury, and R. L. Detwiler (2013), High-resolution fracture apertuire mapping using optical profilometry, Water. Resour. Res., 49, 1 7, doi: /wrcr Brown, S. R. (1987), Fluid flow thorough rock joints: The effect of surface roughness, J. Geophys. Res., 92(B2), , doi: / JB029iB02p Brush, D. J., and N. R. Thomson (2003), Fluid flow in synthetic rough-walled fractures: Navier-Stokes, Stokes, and local cubic law assumptions, Water Resour. Res., 39(4), 1085, doi: /2002wr Chen, Z., S. P. Narayan, Z. Yang, and S. S. Rahman (2000), An experimental investigation of hydraulic behavior of fractures and joints in granitic rock, Int. J. Rock Mech. Min. Sci., 37, Ellsworth, W. L. (2013), Injection-induced earthquakes, Science, 341(6142), doi: /science Evans, F. E., A. Genter, and J. Sausse (2005), Permeability creation and damage due to massive fluid injections into granite at 3.5 km at Soultz; 1. Borehole observations, J. Geophys. Res., 110, B04203, doi: /2004jb Guglielmi, Y., F. Cappa, J. F. Avouac, P. Henry, and D. Elsworth (2015), Seismicity triggered by fluid injection-induced aseismic slip, Science, 348(6240), Hanks, T. C., and H. Kanamori (1979), A moment magnitude scale, J. Geophys. Res., 84(B5), , doi: /jb084ib05p Häring, M. O., U. Schanz, F. Ladner, and B. C. Dyer (2008), Characterisation of the Basel 1 enhanced geothermal system, Geothermics, 37, Ishibashi, T., N. Watanabe, N. Hirano, A. Okamoto, and N. Tsuchiya (2015), Beyond-laboratory-scale prediction for channeling flows through subsurface rock fractures with heterogeneous aperture distributions revealed by laboratory evaluation, J. Geophys. Res. Solid Earth, 120, , doi: /2014jb Javadi, M., M. Sharifzadeh, K. Shahriar, and Y. Mitani (2014), Critical Reynolds number for nonlinear flow through rough-walled fractures: The role of shear process, Water Resour. Res., 50, , doi: /2013wr ISHIBASHI ET AL. LINKING MEQS TO PERMEABILITY CHANGE 7492
8 Kanamori, H., and D. L. Anderson (1975), Theoretical basis of some empirical relations in seismology, Bull. Seismol. Soc. Am., 65, Leonard, M. (2010), Earthquake fault scaling: Self-consistent relating of rupture length, width, average displacement, and moment release, Bull. Seismol. Soc. Am., 100(5A), Majer, E. L., R. Baria, M. Stark, S. Oates, J. Bommer, B. Smith, and H. Asanuma (2007), Induced seismicity associated with enhanced geothermal systems, Geothermics, 36, Marone, C. (1998), Laboratory-derived friction laws and their application to seismic faulting, Annu. Rev. Earth Planet. Sci., 26, McClure, M. W., and R. N. Horne (2012), Investigation of injection-induced seismicity using a coupled fluid flow and rate/state friction model, Geophysics, 76(6), Nemoto, K., N. Watanabe, N. Hirano, and N. Tsuchiya (2009), Direct measurement of contact area and stress dependence of anisotropic flow through rock fracture with heterogeneous aperture distribution, Earth Planet. Sci. Lett., 281, Olsson, R., and N. Barton (2001), An improved model for hydromechanical coupling during shearing of rock joints, Int. J. Rock Mech. Min. Sci., 38, Preisig, G., E. Eberhardt, V. Gischig, V. Roche, M. V. Baan, B. Valley, P. K. Kaiser, D. Duff, and R. Lowther (2015), Development of connected permeability in massive crystalline rocks through hydraulic fracture propagation and shearing accompanying fluid injection, Geofluids, 15, Renard, F., T. Candela, and E. Bouchaud (2013), Constant dimensionality of fault roughness from the scale of micro-fracture to the scale of continents, Geophys. Res. Lett., 40, 83 87, doi: /2012gl Shelly, D. R., D. P. Hill, F. Massin, J. Farrell, R. B. Smith, and T. Taira (2013), A fluid-driven earthquake swarm on the margin of the Yellowstone caldera, J. Geophys. Res. Solid Earth, 118, , doi: /jgrb Tsang, C. F., and I. Neretnieks (1998), Flow channeling in heterogeneous fractured rocks, Rev. Geophys., 36(2), , doi: / 97RG Watanabe, N., N. Hirano, and N. Tsuchiya (2009), Diversity of channeling flow in heterogeneous aperture distribution inferred from integrated experimental-numerical analysis on flow through shear fracture in granite, J. Geophys. Res., 114, B04208, doi: / 2008JB Zoback, M. D., and S. M. Gorelick (2012), Earthquake triggering and large-scale geologic storage of carbon dioxide, Proc. Natl. Acad. Sci. U.S.A., 109(26), 10,164 10,168. ISHIBASHI ET AL. LINKING MEQS TO PERMEABILITY CHANGE 7493
Mapping the Preferential Flow Paths within a Fractured Reservoir
Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Mapping the Preferential Flow Paths within a Fractured Reservoir Takuya Ishibashi 1, Noriaki Watanabe 2, Tetsuya Tamagawa
More informationPredicting the Channeling Flows through Fractures at Various Scales
Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Predicting the Channeling Flows through Fractures at Various Scales Noriaki Watanabe 1, Takuya Ishibashi 2 and Noriyoshi
More informationRelative Permeability Measurement and Numerical Modeling of Two-Phase Flow Through Variable Aperture Fracture in Granite Under Confining Pressure
GRC Transactions, Vol. 36, 2012 Relative Permeability Measurement and Numerical Modeling of Two-Phase Flow Through Variable Aperture Fracture in Granite Under Confining Pressure Noriaki Watanabe, Keisuke
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 informationGeoFlow: A novel model simulator for prediction of the 3-D channeling flow in a rock fracture network
WATER RESOURCES RESEARCH, VOL. 48, W07601, doi:10.1029/2011wr011226, 2012 GeoFlow: A novel model simulator for prediction of the 3-D channeling flow in a rock fracture network Takuya Ishibashi, 1 Noriaki
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 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 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 informationdoi: /2006WR005411
doi: 10.1029/2006WR005411 WATER RESOURCES RESEARCH, VOL. 44, W06412, doi:10.1029/2006wr005411, 2008 Determination of aperture structure and fluid flow in a rock fracture by high-resolution numerical modeling
More informationYusuke 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 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 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 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 informationThermo-Hydro-Mechanical modeling of EGS using COMSOL Multiphysics
PROCEEDINGS, Fourtieth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 26-28, 2015 SGP-TR-204 Thermo-Hydro-Mechanical modeling of EGS using COMSOL Multiphysics
More informationStress-Permeability Relationships in Low Permeability Systems: Application to Shear Fractures
PROCEEDINGS, Thirty-Ninth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 24-26, 2014 SGP-TR-202 Stress-Permeability Relationships in Low Permeability Systems:
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 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 informationExperimental Investigation of Injection-driven Shear Slip and Permeability Evolution in Granite for EGS Stimulation
PROCEEDINGS, 42nd Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 13-15, 217 SGP-TR-212 Experimental Investigation of Injection-driven Shear Slip and Permeability
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 informationPUBLICATIONS. Geophysical Research Letters. Postinjection Normal Closure of Fractures as a Mechanism for Induced Seismicity
PUBLICATIONS Geophysical Research Letters RESEARCH LETTER Key Points: Normal closure of stimulated fractures after the termination of injection enhances postinjection seismicity Processes are strongly
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 informationThermal Performance Implications of Flow Orientation Relative to Fracture Shear Offset in Enhanced Geothermal Systems
PROCEEDINGS, 44th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2019 SGP-TR-214 Thermal Performance Implications of Flow Orientation Relative to
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 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 informationNew Concept Discrete Fracture Network Model Simulator, GeoFlow, and Three Dimensional Channeling Flow in Fracture Network
Proceedings World Geotheral Congress 2010 Bali, Indonesia, 25-29 April 2010 New Concept Discrete Fracture Network Model Siulator, GeoFlow, and Three Diensional Channeling Flow in Fracture Network Nobuo
More informationswisstopo, Swiss Geological Survey, Wabern, Switzerland, Lawrence Berkeley National Laboratory, Earth & Environmental Sciences Area, Berkeley, USA
Aseismic fault slip and leakage preceding an earthquake induced during an in-situ fault reactivation experiment in the Opalinus Clay, Mont Terri rock laboratory, Switzerland Christophe Nussbaum 1, Yves
More informationDepth variation of coseismic stress drop explains bimodal earthquake magnitude-frequency distribution
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L24301, doi:10.1029/2008gl036249, 2008 Depth variation of coseismic stress drop explains bimodal earthquake magnitude-frequency distribution
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 informationInvestigation of Injection-Triggered Slip on Basement Faults: Role of Fluid Leakoff on Post Shut-In Seismicity
Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Investigation of Injection-Triggered Slip on Basement Faults: Role of Fluid Leakoff on Post Shut-In Seismicity Jack Norbeck
More informationDi#erences in Earthquake Source and Ground Motion Characteristics between Surface and Buried Crustal Earthquakes
Bull. Earthq. Res. Inst. Univ. Tokyo Vol. 2+,**0 pp.,/3,00 Di#erences in Earthquake Source and Ground Motion Characteristics between Surface and Buried Crustal Earthquakes Paul Somerville* and Arben Pitarka
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 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 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 informationSPATIAL AND TEMPORAL DISTRBUTION OF LARGER SEISMIC EVENTS AT EUROPEAN AND AUSTRALIAN HDR SITES
PROCEEDINGS, Thirty-First Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 30-February 1, 2006 SGP-TR-179 SPATIAL AND TEMPORAL DISTRBUTION OF LARGER SEISMIC
More informationSeismicity-permeability coupling in the behavior of gas shales, CO 2 storage and deep geothermal energy
DOI 10.1007/s40948-017-0051-9 ORIGINAL ARTICLE Seismicity-permeability coupling in the behavior of gas shales, CO 2 storage and deep geothermal energy Y. Fang. C. Wang. D. Elsworth. T. Ishibashi Received:
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 3,500 108,000 1.7 M Open access books available International authors and editors Downloads Our
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 informationThe Coso Geothermal Area: A Laboratory for Advanced MEQ Studies for Geothermal Monitoring
The Coso Geothermal Area: A Laboratory for Advanced MEQ Studies for Geothermal Monitoring Bruce R. Julian U. S. Geological Survey, Menlo Park, CA 94025 USA julian@usgs.gov Gillian R. Foulger Dept. Earth
More informationScaling relations of seismic moment, rupture area, average slip, and asperity size for M~9 subduction-zone earthquakes
GEOPHYSICAL RESEARCH LETTERS, VOL. 4, 7 74, doi:1.12/grl.976, 213 Scaling relations of seismic moment, rupture area, average slip, and asperity size for M~9 subduction-zone earthquakes Satoko Murotani,
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 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 informationDevelopment of connected permeability in massive crystalline rocks through hydraulic fracture propagation and shearing accompanying fluid injection
Geofluids (2014) doi: 10.1111/gfl.12097 Development of connected permeability in massive crystalline rocks through hydraulic fracture propagation and shearing accompanying fluid injection G. PREISIG 1,
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 informationModeling Microseismic Activity in the Newberry Enhanced Geothermal System
PROCEEDINGS, 42nd Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 13-15, 2017 SGP-TR-212 Modeling Microseismic Activity in the Newberry Enhanced Geothermal
More informationExtending the magnitude range of seismic reservoir monitoring by Utilizing Hybrid Surface Downhole Seismic Networks
Extending the magnitude range of seismic reservoir monitoring by Utilizing Hybrid Surface Downhole Seismic Networks Gisela Viegas*, ESG, Kingston, Ontario, Canada Gisela.Fernandes@esgsolutions.com and
More informationMULTIPLE-WELL, MULTIPLE-RESERVOIR, LONG TERM THERMAL MODELING AT SOULTZ EGS SITE
PROCEEDINGS, Thirty-Eighth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2013 SGP-TR-198 MULTIPLE-WELL, MULTIPLE-RESERVOIR, LONG TERM THERMAL MODELING
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 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 informationThe Pennsylvania State University. The Graduate School. College of Earth and Mineral Sciences ANALYSIS OF INDUCED SEISMICITY AND HEAT TRANSFER
The Pennsylvania State University The Graduate School College of Earth and Mineral Sciences ANALYSIS OF INDUCED SEISMICITY AND HEAT TRANSFER IN GEOTHERMAL RESERVOIRS BY COUPLED SIMULATION A Dissertation
More informationFault Specific, Dynamic Rupture Scenarios for Strong Ground Motion Prediction
Fault Specific, Dynamic Rupture Scenarios for Strong Ground Motion Prediction H. Sekiguchi Disaster Prevention Research Institute, Kyoto University, Japan Blank Line 9 pt Y. Kase Active Fault and Earthquake
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 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 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 informationMathematical Modelling of a Fault Slip Induced by Water Injection
Mathematical Modelling of a Fault Slip Induced by Water Injection T. S. Nguyen, 1 J.Rutqvist 2 and Y. Gugliemi 2 1 Canadian Nuclear Safety Commission 2 Lawrence Berkeley National Laboratory ComGeo IV Symposium
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 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 informationEstimation of fracture flow parameters through numerical analysis of hydromechanical pressure pulses
WATER RESOURCES RESEARCH, VOL. 44,, doi:10.1029/2008wr007015, 2008 Estimation of fracture flow parameters through numerical analysis of hydromechanical pressure pulses Frédéric Cappa, 1 Yves Guglielmi,
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 informationSize effect on aperture and permeability of a fracture as estimated in large synthetic fractures
International Journal of Rock Mechanics & Mining Sciences 43 (2006) 726 755 www.elsevier.com/locate/ijrmms Size effect on aperture and permeability of a fracture as estimated in large synthetic fractures
More informationIntegrating Lab and Numerical Experiments to Investigate Fractured Rock
Integrating Lab and Numerical Experiments to Investigate Fractured Rock Bradford H. Hager Director, Earth Resources Laboratory and Cecil and Ida Green Professor Department of Earth, Atmospheric and Planetary
More informationFracture induced shear wave splitting in a source area of triggered seismicity by the Tohoku-oki earthquake in northeastern Japan.
Fracture induced shear wave splitting in a source area of triggered seismicity by the Tohoku-oki earthquake in northeastern Japan Masahiro Kosuga 1 1. Corresponding Author. Professor, Graduate School of
More informationInduced seismicity in Enhanced Geothermal Systems (EGS) in Alsace, France. Jean Schmittbuhl 1
Induced seismicity in Enhanced Geothermal Systems (EGS) in Alsace, France Jean Schmittbuhl 1 Olivier Lengliné 1, Mohamed Boubacard 1, Lucie Lamourette 1, Lilas Vivin 1, Nicolas Cuenot 2 1 EOST-IPGS, Université
More informationSynthetic Seismicity Models of Multiple Interacting Faults
Synthetic Seismicity Models of Multiple Interacting Faults Russell Robinson and Rafael Benites Institute of Geological & Nuclear Sciences, Box 30368, Lower Hutt, New Zealand (email: r.robinson@gns.cri.nz).
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 informationModeling Injection-Induced Seismicity with the Physics-Based Earthquake Simulator RSQSim
Modeling Injection-Induced Seismicity with the Physics-Based Earthquake Simulator RSQSim by James H. Dieterich, Keith B. Richards-Dinger, and Kayla A. Kroll INTRODUCTION Although the phenomenon of earthquakes
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 informationSimulation of earthquake rupture process and strong ground motion
Simulation of earthquake rupture process and strong ground motion Takashi Miyatake (1) and Tomohiro Inoue (2) (1) Earthquake Research Institute, University of Tokyo, Yayoi, Bunkyo, Tokyo, 113-0032, Japan
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lecture 20, 30 Nov. 2017 www.geosc.psu.edu/courses/geosc508 Seismic Spectra & Earthquake Scaling laws. Seismic Spectra & Earthquake Scaling laws. Aki, Scaling law
More informationAmélie Neuville (1,2), Renaud Toussaint (2), Eirik Flekkøy (1), Jean Schmittbuhl (2)
Amélie Neuville (1,2), Renaud Toussaint (2), Eirik Flekkøy (1), Jean Schmittbuhl (2) Thermal exchanges between a hot fractured rock and a cold fluid Deep geothermal systems Enhanced Geothermal Systems
More informationInjection-Triggered Seismicity: An Investigation of Porothermoelastic Effects Using a Rate-and- State Earthquake Model
PROCEEDINGS, Fourtieth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 26-28, 2015 SGP-TR-204 Injection-Triggered Seismicity: An Investigation of Porothermoelastic
More informationEffect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics
Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics Naoyuki Kato (1) and Tomowo Hirasawa (2) (1) Geological
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 informationNumerical study on multi-scaling earthquake rupture
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L266, doi:1.129/23gl1878, 24 Numerical study on multi-scaling earthquake rupture Hideo Aochi Institut de Radioprotection et de Sûreté Nucléaire, France Satoshi Ide
More informationShort Note Fault Slip Velocities Inferred from the Spectra of Ground Motions
Bulletin of the Seismological Society of America, Vol. 99, No. 2A, pp. 876 883, April 2009, doi: 10.1785/0120080008 Short Note Fault Slip Velocities Inferred from the Spectra of Ground Motions by N. Ani
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 informationGeothermal Induced Seismicity : Summary of International Experience IEA-GIA Environmental Mitigation Workshop Taupo, June 2012
Geothermal Induced Seismicity : Summary of International Experience IEA-GIA Environmental Mitigation Workshop Taupo, 15-16 June 2012 1 Chris Bromley, Wairakei Research Centre, Taupo, New Zealand 1.0 Background
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 information3D HM-DEM model for Hydro-Fracturing
3D HM-DEM model for Hydro-Fracturing E. Papachristos, F.V. Donzé & B. Chareyre Laboratoire Sols, Solides, Structures, Grenoble, France efthymios.papachristos@3sr-grenoble.fr,, frederic.donze@3srgrenoble.fr,
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 informationHigh fluid pressure and triggered earthquakes in the enhanced geothermal system in Basel, Switzerland
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 117,, doi:10.1029/2011jb008980, 2012 High fluid pressure and triggered earthquakes in the enhanced geothermal system in Basel, Switzerland Toshiko Terakawa, 1 Stephen
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 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 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 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 informationSHEAR-SLIP ANALYSIS IN MULTIPHASE FLUID-FLOW RESERVOIR ENGINEERING APPLICATIONS USING TOUGH-FLAC
PROCEEDINGS, TOUGH Symposium 6 Lawrence Berkeley National Laboratory, Berkeley, California, May 5 7, 6 SHEAR-SLIP ANALYSIS IN MULTIPHASE FLUID-FLOW RESERVOIR ENGINEERING APPLICATIONS USING TOUGH-FLAC Jonny
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 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 informationDISCRETE FRACTURE NETWORK MODELLING OF HYDRAULIC FRACTURING IN A STRUCTURALLY CONTROLLED AREA OF THE MONTNEY FORMATION, BC
DISCRETE FRACTURE NETWORK MODELLING OF HYDRAULIC FRACTURING IN A STRUCTURALLY CONTROLLED AREA OF THE MONTNEY FORMATION, BC Steve Rogers Golder Associates Ltd Pat McLellan McLellan Energy Advisors Inc Gordon
More informationThe Influence of Shear and Deviatoric Stress on the Evolution of Permeability in Fractured Novaculite and Diorite
Igor Faoro, Andre Niemeijer, Chris Marone, Derek Elsworth The Influence of Shear and Deviatoric Stress on the Evolution of Permeability in Fractured Novaculite and Diorite Department of Energy and Geo-Environmental
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 information3D modeling of the cycle of a great Tohoku oki earthquake, considering frictional behavior at low to high slip velocities
GEOPHYSICAL RESEARCH LETTERS, VOL. 38,, doi:10.1029/2011gl049308, 2011 3D modeling of the cycle of a great Tohoku oki earthquake, considering frictional behavior at low to high slip velocities B. Shibazaki,
More informationCOUPLED HYDRO-MECHANICAL MODELLING OF THE GPK3 RESERVOIR STIMULATION AT THE EUROPEAN EGS SITE SOULTZ-SOUS-FORÊTS
PROCEEDINGS, Thirtieth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 31-February 2, 25 SGP-TR-176 COUPLED HYDRO-MECHANICAL MODELLING OF THE GPK3 RESERVOIR
More informationA THREE-DIMENSIONAL STOCHASTIC FRACTURE NETWORK MODEL FOR GEOTHERMAL RESERVOIR STIMULATION
PROCEEDINGS, Thirty-Sixth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 31 - February 2, 2011 SGP-TR-191 A THREE-DIMENSIONAL STOCHASTIC FRACTURE NETWORK
More informationSupplementary Materials for
advances.sciencemag.org/cgi/content/full/4/3/eaao4915/dc1 Supplementary Materials for Global variations of large megathrust earthquake rupture characteristics This PDF file includes: Lingling Ye, Hiroo
More informationWhat allows seismic events to grow big?: Insights from fault roughness and b-value analysis in stick-slip experiments
What allows seismic events to grow big?: Insights from fault roughness and b-value analysis in stick-slip experiments T. H. W. Goebel 1, G. Kwiatek 2, T. W. Becker 3, G. Dresen 2, E.E. Brodsky 1 1 UC Santa
More informationNumerical simulation of seismic cycles at a subduction zone with a laboratory-derived friction law
Numerical simulation of seismic cycles at a subduction zone with a laboratory-derived friction law Naoyuki Kato (1), Kazuro Hirahara (2) and Mikio Iizuka (3) (1) Earthquake Research Institute, University
More informationSheared Fracture Conductivity
PROCEEDINGS, Thirty-Ninth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 24-26, 2014 SGP-TR-202 Sheared Fracture Conductivity Ravindra Bhide 1, Tyler Gohring
More informationKeywords: induced seismicity, fault heterogeneity, rate-and-state friction, poroelasticity, earthquake modeling
POCEEDINGS, 41st Workshop on Geothermal eservoir Engineering Stanford University, Stanford, California, February 22-24, 2016 SGP-T-209 evisiting Stimulation Mechanism at Fenton Hill and an Investigation
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 informationIntegration of seismic and fluid-flow data: a two-way road linked by rock physics
Integration of seismic and fluid-flow data: a two-way road linked by rock physics Abstract Yunyue (Elita) Li, Yi Shen, and Peter K. Kang Geologic model building of the subsurface is a complicated and lengthy
More information