Using the migration of the induced seismicity as a constraint for fractured Hot Dry Rock reservoir modelling

The challenge for Hot Dry Rock technology is to develop a heat exchanger in deep hot rock masses and to circulate a fluid in order to extract its energy to be used at the ground surface. The present day strategy is to take advantage of natural fractures that pre-exist at these depths and to improve their hydraulic properties. The extension of the area with enhanced properties must then be evaluated so that the best locations for further boreholes can be proposed. To date, this development procedure is based on fluid injection at high rates, forcing hydro-mechanical interactions to take place along pressurised fractures. These pore-pressure-driven mechanisms are accompanied by seismic activity. Assuming the validity of poro-elastic theory in the fractured host rock, some authors have derived the virgin hydraulic diffusivity of the fractured reservoir from the analysis of the spatio-temporal growth of the induced seismicity. The present work is aimed at verifying this approach using a numerical code to solve directly for hydro-mechanical interactions in random fracture networks. Our approach assumes that the seismic activity is controlled by a Coulomb shear criterion and we show how the interpretation of spreading rate of the modelled shear activity in a given network coincides with the upscaled virgin hydraulic diffusivity of the same fracture network, calculated from an independent numerical procedure at the reservoir scale. Therefore, it is shown that the direct analysis of the seismicity migration is appropriate to give reliable estimates of virgin hydraulic and mechanic parameters. These parameters can then be used for performing any further quantitative analysis of a reservoir open to the far field. This is of importance, as fluid mass-balance in multi-well exploitation systems will include the contribution from areas surrounding the stimulated zone.

[1]  D. Bruel Impact of Induced Thermal Stresses During Circulation Tests in an Engineered Fractured Geothermal Reservoir: Example of the Soultz-Sous-Forêts European Hot Fractured Rock Geothermal Project, Rhine Graben, France , 2002 .

[2]  F. Cornet,et al.  Integrated stress determination by joint inversion of hydraulic tests and focal mechanisms , 1994 .

[3]  F. H. Cornet,et al.  Analysis of induced seismicity for fault zone identification , 1993 .

[4]  R. J. Pine,et al.  Downward migration of shearing in jointed rock during hydraulic injections , 1984 .

[5]  Jean-Raynald de Dreuzy,et al.  Hydraulic properties of two‐dimensional random fracture networks following a power law length distribution: 1. Effective connectivity , 2001 .

[6]  F. Cornet Comment on 'Large-scale in situ permeability tensor of rocks from induced microseismicity' by S. A. Shapiro, P. Audigane and J.-J. Royer , 2000 .

[7]  Paul Segall,et al.  Joint formation in granitic rock of the Sierra Nevada , 1983 .

[8]  J. Mallet,et al.  3D Fracture Modeling Integrating Geomechanics and Geologic Data. , 2004 .

[9]  Pascal Audigane,et al.  Large-scale in situ permeability tensor of rocks from induced microseismicity , 1999 .

[10]  Olivier Bour,et al.  On the connectivity of three‐dimensional fault networks , 1998 .

[11]  Tomáš Fischer,et al.  Slip‐generated patterns of swarm microearthquakes from West Bohemia/Vogtland (central Europe): Evidence of their triggering mechanism? , 2005 .

[12]  T. Cladouhos,et al.  Are fault growth and linkage models consistent with power-law distributions of fault lengths? , 1996 .

[13]  A.C.G. Nagelhout,et al.  Investigating fault slip in a model of an underground gas storage facility , 1997 .

[14]  Michael Fehler,et al.  Stress control of seismicity patterns observed during hydraulic fracturing experiments at the Fenton Hill hot dry rock geothermal energy site, New Mexico , 1987 .

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

[16]  Serge A. Shapiro,et al.  Microseismic monitoring of borehole fluid injections: Data modeling and inversion for hydraulic properties of rocks , 2003 .

[17]  J. Thovert,et al.  Macroscopic permeability of three-dimensional fracture networks with power-law size distribution. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[18]  M. Nepveu,et al.  Correlation between hydrocarbon reservoir properties and induced seismicity in the Netherlands , 2006 .

[19]  H. Kanamori,et al.  A moment magnitude scale , 1979 .

[20]  S. L. Crouch Solution of plane elasticity problems by the displacement discontinuity method. I. Infinite body solution , 1976 .

[21]  H. Takahashi,et al.  Progress toward a stochastic rock mechanics model of engineered geothermal systems , 1996 .

[22]  D. Bruel Modelling Hydraulic Jacking-tests On a Pre-existing Fracture System , 1999 .

[23]  H. H. Einstein 5th International Congress on Rock Mechanics , 1983 .

[24]  N. Odling,et al.  Scaling of fracture systems in geological media , 2001 .

[25]  Hiroaki Niitsuma,et al.  Detailed Fracture System of the Soultz-sous-Forêts HDR Field Evaluated Using Microseismic Multiplet Analysis , 2002 .

[26]  P. R. La Pointe,et al.  Derivation of parent fracture population statistics from trace length measurements of fractal fracture populations , 2002 .

[27]  T. Yamashita Simulation of seismicity due to fluid migration in a fault zone , 1998 .

[28]  Hiroshi Asanuma,et al.  Creation of an HDR reservoir at 5000 m depth at the European HDR project , 2006 .

[29]  William L. Ellsworth,et al.  Monitoring velocity variations in the crust using earthquake doublets: An application to the Calaveras Fault, California , 1984 .

[30]  M. Fehler,et al.  More Than Cloud New techniques for characterizing reservoir structure using induced seismicity , 2001 .

[31]  Philippe Davy,et al.  On the Frequency-Length Distribution of the San Andreas Fault System , 1993 .

[32]  Y. Okada Internal deformation due to shear and tensile faults in a half-space , 1992, Bulletin of the Seismological Society of America.

[33]  Serge A. Shapiro,et al.  Characterization of fluid transport properties of reservoirs using induced microseismicity , 2002 .

[34]  A. Rejeb,et al.  Hydromechanical effects of shaft sinking at the Sellafield site , 2001 .

[35]  N. Barton,et al.  Strength, deformation and conductivity coupling of rock joints , 1985 .

[36]  C. Fidelibus,et al.  Derivation of equivalent pipe network analogues for three‐dimensional discrete fracture networks by the boundary element method , 1999 .

[37]  Nicolas Cuenot,et al.  Characterization of fluid transport properties of the Hot Dry Rock reservoir Soultz-2000 using induced microseismicity , 2004 .

[38]  T. Tamagawa,et al.  Numerical Simulation of Hydraulic Shearing in Fractured Reservoir , 2005 .

[39]  Gudmundur S. Bodvarsson,et al.  Hydraulic conductivity of rock fractures , 1996 .

[40]  M. Fehler,et al.  Spatial Correlation of Seismic Slip at the HDR-Soultz Geothermal Site: Qualitative Approach , 2000 .

[41]  S. Talebi Seismicity Associated with Mines, Reservoirs and Fluid Injections , 1998 .