Forward Induced Seismic Hazard Assessment (FISHA) Based on a Hydromechanical Coupled Fracture Mechanics Earthquake Model

The occurrence of an Mw 3.2 induced seismic event in 2006 during the injection of fluid at the Basel geothermal site in Switzerland initiated an on-going discussion on the potential hazard and risk of hydraulic stimulation in Europe. So far, development of mitigation strategies of induced Seismic Events of Economic Concern (SEECo) has become an important issue. Statistical methods have been developed to understand the induced seismicity from which appropriate strategy for mitigating effect of SEECo can be established. However, all these methods require catalogues of monitored seismicity. Therefore, they cannot deliver a priori strategies for stimulation treatment. Here we present the Forward Induced Seismic Hazard Assessment (FISHA). This procedure converts the output of a hydromechanical coupled and fracture mechanics based model into probabilistic seismic hazard of induced seismicity in terms of time-dependent occurrence rate of induced seismic events. FISHA is applied to the simulation results of various stimulation scenarios, where injection flow rate, duration and style of injection, as well as reservoir parameters are varied. The geothermal reservoir is modelled using a discrete element fracture network that is coupled to a routine of viscous fluid flow in porous media generating flow-driven failure of rock matrix and pre-existing fractures in Mode I (tensile) and Mode II (shear). The output is a synthetic seismicity catalogue. This catalogue includes hypocentre, occurrence time and magnitude of induced seismic events. The synthetic catalogue is then used to determine the magnitude completeness Mc and to estimate the time-dependent parameters a and b of the Gutenberg-Richter frequency-magnitude relation using statistical methods. These parameters are then applied to estimate the time-dependent occurrence rate of SEECo for each stimulation scenario. In this way, we can evaluate the stimulation scenarios in terms of hazard according to the corresponding rate of SEECo. The advantage of FISHA is then the assessment of the occurrence rate of SEECo resulting from various stimulation scenarios, which gives a priori information to operator how to mitigate the risk of occurrence of SEECo.

[1]  R. P. Young,et al.  Moment tensors and micromechanical models , 2002 .

[2]  R. P. Young,et al.  Dynamic modelling of induced seismicity , 2004 .

[3]  Jeoung Seok Yoon,et al.  Numerical investigation on optimized stimulation of intact and naturally fractured deep geothermal reservoirs using hydro-mechanical coupled discrete particles joints model , 2014 .

[4]  Gottfried Grünthal,et al.  Induced seismicity related to geothermal projects versus natural tectonic earthquakes and other types of induced seismic events in Central Europe , 2014 .

[5]  Friedemann Wenzel,et al.  Probability of earthquake occurrence and magnitude estimation in the post shut-in phase of geothermal projects , 2011, Journal of Seismology.

[6]  O. Stephansson,et al.  Hydro-Mechanical Coupled Discrete Element Modeling of Geothermal Reservoir Stimulation and Induced Seismicity , 2013 .

[7]  J. Dieterich A constitutive law for rate of earthquake production and its application to earthquake clustering , 1994 .

[8]  Fabio Moia,et al.  A survey of the induced seismic responses to fluid injection in geothermal and CO2 reservoirs in Europe , 2012 .

[9]  G. Grünthal,et al.  Forward induced seismic hazard assessment: application to a synthetic seismicity catalogue from hydraulic stimulation modelling , 2014, Journal of Seismology.

[10]  C. Tsang,et al.  Estimating maximum sustainable injection pressure during geological sequestration of CO2 using coupled fluid flow and geomechanical fault-slip analysis , 2006 .

[11]  Arno Zang,et al.  Forward modelling of seismicity rate changes in georeservoirs with a hybrid geomechanical–statistical prototype model , 2014 .

[12]  Thomas Kohl,et al.  Predictive modeling of reservoir response to hydraulic stimulations at the European EGS site Soultz-sous-Forêts , 2007 .

[13]  Jeoung Seok Yoon,et al.  Fatigue hydraulic fracturing by cyclic reservoir treatment enhances permeability and reduces induced seismicity , 2013 .

[14]  Jeoung Seok Yoon,et al.  SIMULATION OF HYDRAULIC STIMULATION OF FRACTURED RESERVOIR AND INDUCED SEISMICITY USING DISCRETE ELEMENT-FRACTURE NETWORK MODEL , 2013 .

[15]  Julian J. Bommer,et al.  Induced seismicity associated with Enhanced Geothermal Systems , 2007 .

[16]  Reinhard Jung,et al.  A numerical model for fluid injection induced seismicity at Soultz-sous-Forêts , 2010 .

[17]  Markus Häring,et al.  Characterisation of the Basel 1 enhanced geothermal system , 2008 .

[18]  R. Horne,et al.  Investigation of injection-induced seismicity using a coupled fluid flow and rate/state friction model , 2011 .

[19]  M. King Hubbert,et al.  Mechanics of Hydraulic Fracturing , 1972 .

[20]  Jacques R. Ernst,et al.  Earthquake focal mechanisms of the induced seismicity in 2006 and 2007 below Basel (Switzerland) , 2009 .

[21]  O. Heidbach,et al.  Fluid-induced microseismicity in pre-stressed rock masses , 2010 .

[22]  Domenico Giardini,et al.  Earthquakes Induced By the Stimulation of an Enhanced Geothermal System Below Basel (Switzerland) , 2009 .

[23]  Serge A. Shapiro,et al.  Seismotectonic state of reservoirs inferred from magnitude distributions of fluid-induced seismicity , 2012, Journal of Seismology.

[24]  G. Grünthal,et al.  FISHA-FORWARD INDUCED SEISMIC HAZARD ASSESSMENT APPLICATION TO SYNTHETIC SEISMICITY CATALOG GENERATED BY HYDRULIC STIMULATION MODELING , 2013 .

[25]  J. Bommer,et al.  Control of hazard due to seismicity induced by a hot fractured rock geothermal project , 2006 .

[26]  J. Altmann,et al.  Poroelastic contribution to the reservoir stress path , 2010 .

[27]  T. Bérard,et al.  How close to failure is a granite rock mass at a 5 km depth , 2007 .

[28]  Ivan G. Wong,et al.  Protocol for Addressing Induced Seismicity Associated with Enhanced Geothermal Systems , 2012 .

[29]  Friedemann Wenzel,et al.  Seismogenic index and magnitude probability of earthquakes induced during reservoir fluid stimulations , 2010 .

[30]  J. Labuz,et al.  Mohr–Coulomb Failure Criterion , 2012, Rock Mechanics and Rock Engineering.

[31]  Dominique Bruel,et al.  Using the migration of the induced seismicity as a constraint for fractured Hot Dry Rock reservoir modelling , 2007 .