Particle Based Discrete Element Modeling of Hydraulic Stimulation of Geothermal Reservoirs, Induced Seismicity and Fault Zone Deformation
0
조회
0
다운로드
2
인용
0
참고문헌
키워드
초록
This numerical study investigates seismicity and fault slip induced by fluid injection in deep geothermal reservoir with pre-existing fractures and fault. Particle Flow Code 2D is used with additionally implemented hydro-mechanical coupled fluid flow algorithm and acoustic emission moment tensor inversion algorithm. The output of the model includes spatio-temporal evolution of induced seismicity (hypocenter locations and magnitudes) and fault deformation (failure and slip) in relation to fluid pressure distribution. The model is applied to a case of fluid injection with constant rates changing in three steps using different fluid characters, i.e. the viscosity, and different injection locations. In fractured reservoir, spatio-temporal distribution of the induced seismicity differs significantly depending on the viscosity of the fracturing fluid. In a fractured reservoir, injection of low viscosity fluid results in larger volume of induced seismicity cloud as the fluid can migrate easily to the reservoir and cause large number and magnitude of induced seismicity in the post-shut-in period. In a faulted reservoir, fault deformation (co-seismic failure and aseismic slip) can occur by a small perturbation of fracturing fluid (<0.1 MPa) can be induced when the injection location is set close to the fault. The presented numerical model technique can practically be used in geothermal industry to predict the induced seismicity pattern and magnitude distribution resulting from hydraulic stimulation of geothermal reservoirs prior to actual injection operation.
Abstract
This numerical study investigates seismicity and fault slip induced by fluid injection in deep geothermal reservoir with pre-existing fractures and fault. Particle Flow Code 2D is used with additionally implemented hydro-mechanical coupled fluid flow algorithm and acoustic emission moment tensor inversion algorithm. The output of the model includes spatio-temporal evolution of induced seismicity (hypocenter locations and magnitudes) and fault deformation (failure and slip) in relation to fluid pressure distribution. The model is applied to a case of fluid injection with constant rates changing in three steps using different fluid characters, i.e. the viscosity, and different injection locations. In fractured reservoir, spatio-temporal distribution of the induced seismicity differs significantly depending on the viscosity of the fracturing fluid. In a fractured reservoir, injection of low viscosity fluid results in larger volume of induced seismicity cloud as the fluid can migrate easily to the reservoir and cause large number and magnitude of induced seismicity in the post-shut-in period. In a faulted reservoir, fault deformation (co-seismic failure and aseismic slip) can occur by a small perturbation of fracturing fluid (<0.1 MPa) can be induced when the injection location is set close to the fault. The presented numerical model technique can practically be used in geothermal industry to predict the induced seismicity pattern and magnitude distribution resulting from hydraulic stimulation of geothermal reservoirs prior to actual injection operation.
참고문헌 (28)
- Altmann J., "Poroelastic contribution to the reservoir stress path", International Journal of Rock Mechanics and Mining Sciences, 47(7), pp.1104~1113 (2010)
- Al-Busaidi A., "Distinct element modeling of hydraulically fractured Lac du Bonnet granite", Journal of Geophysical Research, 110, pp.B06302~ (2005)[DOI]
- Baisch S., "A numerical model for fluid injection induced seismicity at Soultz-sous-Forêts", International Journal of Rock Mechanics Mining Sciences, 47, pp.405~413 (2010)
- Bruel D, "Using the migration of induced seismicity as a constraint for fractured hot dry rock reservoir modelling", International Journal of Rock Mechanics and Mining Sciences, 44, pp.1106~1117 (2007)
- Cornet F. H., "How close to failure is a granite rock mass at a 5 km depth?", International Journal of Rock Mechanics Mining Sciences, 44, pp.47~66 (2007)
- Gutenberg B., "Earthquake magnitude, intensity, energy and acceleration (second paper)", Bulletin of the Seismological Society of America, 46, pp.105~145 (1956)
- Hakimhashemi A., "Forward induced seismic hazard assessment (FISHA) application to synthetic seismicity catalog generated by hydromechanical modeling of fluid injection", Proceedings of Thirty-Eighth Workshop on Geothermal Reservoir Engineering (2013), Stanford University
- Hanks T. C., "A Moment Magnitude Scale", Journal of Geophysical Research, 84, pp.2348~2350 (1979)
- Hazzard J. F., "Moment tensors and micromechanical models", Tectonophysics, 356, pp.181~197 (2002)
- Hazzard J. F., "Dynamic modeling of induced seismicity", International Journal of Rock Mechanics Mining Sciences, 41, pp.1365~1376 (2004)
- Hazzard J. F., "Numerical modeling of seismicity induced by fluid injection in a fractured reservoir. In: Mining and Tunnel Innovation and Opportunity", Proceedings of the 5th North American Rock Mechanics Symposium, pp.1023~1030 (2002)
- Hillis R, "Pore pressure", stress coupling and its implications for seismicity, 2000, pp.31, 448~454
- THM-issues in repository rock-Thermal, mechanical, thermo-mechanical and hydro-mechanical evolution of the rock at the Forsmark and Laxemar sites, "Hökmark H.", SKB-Swedish Nuclear Fuel and Waste Management Co, ~, pp.Technical Report TR-10-23 (2010)
- Itasca Consulting Group Inc., "PFC2D-Particle Flow Code in 2 Dimensions, Version 4.0" (2008)
- Kohl T., "Predictive modeling of reservoir response to hydraulic stimulations at the European EGS site Soultz-sous-Forets", International Journal of Rock Mechanics and Mining Sciences, 44(8), pp.1118~1131 (2007)
- Kraft T., "Enhanced Geothermal Systems: Mitigating Risk in Urban Areas", EOS, 90(32) (2009)
- Kramer S. L, "Geotechnical Earthquake Engineering", pp.653~ (1996), Prentice-Hall
- Majer E. L., "Induced seismicity associated with Enhanced Geothermal Systems", Geothermics, 36, pp.185~222 (2007)
- McClure M., "Investigation of injectioninduced seismicity using a coupled fluid flow and rate and state friction model", Geophysics, 76(6), pp.WC183~WC200 (2011)
- Mukuhira Y., "Characteristics of large-magnitude microseismic events recorded during and after stimulation of a geothermal reservoir at Basel, Switzerland", Geothermics, 45, pp.1~17 (2013)
- Potyondy D. O., "A bonded-particle model for rock", International Journal of Rock Mechanics Mining Sciences, 41, pp.1329~1364 (2004)
- Rutqvist J., "Estimating maximum sustainable injection pressure during geological sequestration of CO2 using coupled fluid flow and geomechanical fault-slip analysis", Energy Conversion and Management, 48, pp.1798~1807 (2007)
- Schoenball M., "Fluid-induced microseismicity in pre-stressed rock masses", Geophysical Journal International, 180, pp.113~119 (2010)
- Yoon J. S., "Simulating fracture and friction of Aue granite under confined asymmetric compressive test using clumped particle model", International Journal of Rock Mechanics Mining Sciences, 49, pp.68~83 (2012)
- Yoon J. S., "Simulation of hydraulic stimulation of fractured reservoir and induced seismicity using discrete element-fracture network model", Proceedings of Thirty-Eighth Workshop on Geothermal Reservoir Engineering (2013), Stanford University
- Zang A., "Stress Field of the Earth’s Crust" (2010), Springer Science + Business Media B.V.
- Zang A., "Fatigue hydraulic fracturing by cyclic reservoir treatment enhances permeability and reduced induced seismicity", Geophysical Journal International, 195, pp.1282~1287 (2013)
- Zhao X., "Numerical modelling of seismicity induced by fluid injection in naturally fractured reservoirs", Geophysics, 76(6), pp.WC167~WC180 (2011)