Recent development in numerical simulation of enhanced geothermal reservoirs

Abstract

This paper briefly introduces the current state in computer modelling of geothermal reservoir system and then focuses on our research efforts in high performance simulation of enhanced geothermal reservoir system. A novel supercomputer simulation tool has been developing towards simulating the highly non-linear coupled geomechanical-fluid flow-thermal systems involving heterogeneously fractured geomaterials at different spatial and temporal scales. It is applied here to simulate and visualise the enhanced geothermal system (EGS), such as (1) visualisation of the microseismic events to monitor and determine where/how the underground rupture proceeds during a hydraulic stimulation, to generate the mesh using the recorded data for determining the domain of the ruptured zone and to evaluate the material parameters (i.e., the permeability) for the further numerical analysis and evaluation of the enhanced geothermal reservoir; (2) converting the available fractured rock image/fracture data as well as the reservoir geological geometry to suitable meshes/grids and further simulating the fluid flow in the complicated fractures involving the detailed description of fracture dimension and geometry by the lattice Boltzmann method and/or finite element method; (3) interacting fault system simulation to determine the relevant complicated rupture process for evaluating the geological setting and the in-situ reservoir properties; (4) coupled thermo-fluid flow analysis of a geothermal reservoir system for an optimised geothermal reservoir design and management. A few of application examples are presented to show its usefulness in simulating the enhanced geothermal reservoir system.

Keywords

numerical simulation / geothermal / EGS / microseismicity / finite element method / lattice Boltzmann method

References

[1] Bjornsson G, Bodvarsson G. A Survey of Geothermal Reservoir Properties. Geothermics, 1990, 19: 17-27.

[2] Bringemeier D, Wang X, Xing H L, . Modelling of Multiphase Fluid Flow for an Open Pit Development within a Geothermal Active Caldera. Proceedings of the 11th IAEG Congress (International Association for Engineering Geology and the Environment). Auckland, New Zealand, 2010

[3] Brown D, DuTeaux R, Kruger P, . Fluid Circulation and Heat Extraction from Geothermal Reservoirs. Geothermics, 1999, 28: 553-572.

[4] Cox S F, Knackstedt M A, Braun J. Principles of Structural Control on Permeability and Fluid Flow in Hydrothermal Systems. Reviews in Economic Geology, 2001, 14: 1-24.

[5] De Simone S, Vilarrasa V, Carrera J, . Thermal Coupling may Control Mechanical Stability of Geothermal Reservoirs during Cold Water Injection. Physics and Chemistry of the Earth, Parts A/B/C, 2013, 64: 117-126.

[6] Driesner T, Heinrich C A. The System H2O-NaCl. I. Correlation Formulae for Phase Relations in Temperature-Pressure-Composition Space from 0 to 1 000 °C, 0 to 5 000 bar, and 0 to 1 XNaCl. Geochimica et Cosmochimica Acta, 2007, 71: 4880-4901.

[7] Gao J, Xing H L. LBM Simulation of Fluid Flow in Fractured Porous Media with Permeable Matrix. Theor. Appl. Mech. Lett., 2012, 2 032001^(4).

[8] Gao J, Xing H L. High Performance Simulation of Complicated Fluid Flow in 3D Fractured Porous Media Using LBM. 10th International Meeting on High-Performance Computing for Computational Science. Kobe, Japan, 2012

[9] Gao J, Xing H L, Tian Z. Lattice Boltzmann Modeling and Evaluation of Fluid Flow in Heterogeneous Porous Media Involving Multiple Matrix Constituents. Computers & Geosciences, 2013, 62: 198-207.

[10] Held S, Genter A, Kohl T, . Economic Evaluation of Geothermal Reservoir Performance through Modeling the Complexity of the Operating EGS in Soultz-Sous-Forêts. Geothermics, 2014, 51: 270-280.

[11] Ingebritsen S E, Geiger S, Hurwitz S, . Numerical Simulation of Magmatic Hydrothermal Systems. Rev. Geophys., 2010, 48 RG1002.

[12] Jeanne P, Rutqvist J, Vasco D, . Development of a 3D Hydrogeological and Geomechanical Model of an Enhanced Geothermal System Using Microseismic and Ground Deformation Data from a 1-year Injection Program. Proceedings, Thirty-Ninth Workshop on Geothermal Reservoir Engineering Stanford University. Stanford, California, 2014

[13] Li Q, Wei Y N, Liu G, . Combination of CO2 Geological Storage with Deep Saline Water Recovery in Western China: Insights from Numerical Analyses. Applied Energy, 2014, 116: 101-110.

[14] Liu X, Xing H L, Zhang D. Fluid Focusing and Its Link to Vertical Morphological Zonation at the Dajishan Vein-Type Tungsten Deposit, South China. Ore Geology Reviews, 2014, 62: 245-258.

[15] Liu Y H, Xing H L, Guan Z Q. An Indirect Approach for Automatic Generation of Quadrilateral Meshes with Arbitrary Line Constraints. International Journal for Numerical Methods in Engineering, 2011, 87(9): 906-922.

[16] Liu Y, Xing H L. Gurgenci H, Weber R D. Automatic Meshing and Construction of a 3D Reservoir System: From Visualization towards Simulation. Proceedings of the 2010 Australian Geothermal Energy Conference. Nov. 17-19, Adelaide, 2010

[17] Liu Y, Xing H L. A Boundary Focused Quadrilateral Mesh Generation Algorithm for Multi-Material Structures. Journal of Computational Physics, 2013, 232: 516-528.

[18] O’Sullivan M J, Pruess K, Lippmann M J. State of the Art of Geothermal Reservoir Simulation. Geothermics, 2001, 30: 395-429.

[19] Rutqvist J. Pre-Stimulation Coupled THM Modeling Related to the Northwest Geysers EGS Demonstration Project. 38th Annual Workshop on Geothermal Reservoir Engineering, Stanford, 2014

[20] Rybach L. The Future of Geothermal Energy and Its Challenges. Proceedings World Geothermal Congress 2010. 25–29 April, 2010, Bali, 2010

[21] Sanyal S K, Butler S J, Swenson D, . Review of the State-of-the-Art of Numerical Simulation of Enhanced Geothermal Systems. Proceedings World Geothermal Congress 2000, 2000, 3853-3858.

[22] Sukop M C, Thorne D T. Lattice Boltzmann Modeling: An Introduction for Geoscientists and Engineers, 2007

[23] Tester J W, Anderson B J, Batchelor A S, . The Future of Geothermal Energy-Impact of Enhanced Geothermal Systems (EGS) on the United States in the 21st Century, 2006, 358.

[24] Tian Z, Xing H L, Tan Y, . A Coupled Lattice Boltzmann Model for Simulating Reactive Transport in CO₂ Injection. Physica A: Statistical Mechanics and Its Applications, 2014, 403: 155-164.

[25] Wyborn D. The Innamincka EGS Project—10 Years of Operation, 2013

[26] Xing H L. Progress Report: Supercomputer Simulation of Hot Fractured Rock Geothermal Reservoir Systems. ESSCC/ACcESS Technical Report, 2008, 1-48.

[27] Xing H L. Finite Element Simulation of Transient Geothermal Flow in Extremely Heterogeneous Fractured Porous Media. Journal of Geochemical Exploration, 2014, 144(PartA): 168-178.

[28] Xing H L, Liu Y. Automated Quadrilateral Mesh Generation for Digital Image Structure. Theoretical and Applied Mechanics Letters, 2011, 16: 061001.1-061001.3.

[29] Xing H L, Liu Y. Mesh Generation for 3D Geological Reservoirs with Arbitrary Stratigraphic Surface Constraints. Procedia Computer Science, 2014, 29: 897-909.

[30] Xing H L, Makinouchi A. Finite Element Analysis of Sandwich Friction Experimental Model of Rocks. Pure and Applied Geophysics, 2002, 159: 1985-2009.

[31] Xing H L, Makinouchi A. Finite-Element Modelling of Multibody Contact and Its Application to Active Faults. Concurrency and Computation: Practice and Experience, 2002, 14: 431-450.

[32] Xing H L, Makinouchi A. Three Dimensional Finite Element Modelling of Thermomechanical Frictional Contact between Finite Deformation Bodies Using R-Minimum Strategy. Computer Methods in Applied Mechanics and Engineering, 2002, 191: 4193-4214.

[33] Xing H L, Makinouchi A, Mora P. Finite Element Modelling of Interacting Fault System. Physics of the Earth and Planetary Interiors, 2007, 163: 106-121.

[34] Xing H L, Mora P. Construction of an Intraplate Fault System Model of South Australia, and Simulation Tool for the iSERVO Institute Seed Project. Pure and Applied Geophysics, 2006, 163: 2297-2316.

[35] Xing H L, Mora P, Makinouchi A. A Unified Frictional Description and It’s Application to the Simulation of Frictional Instability Using the Finite Element Method. Philosophical Magazine, 2006, 86(21/22): 3453-3475.

[36] Xing H L, Yu W, Zhang J. Xing H L. 3D Mesh Generation in Geocomputing. Advances in Geocomputing, 2009 Berlin & Heidelberg: Springer-Verlag, 27-64.

[37] Xing H L, Zhang J. Finite Element Modelling of Nonlinear Deformation of Rate-Dependent Materials Using a R-Minimum Strategy. Acta Geotechnica., 2009, 4(2): 139-148.

[38] Xiong Y, Fakcharoenphol P, Winterfeld P, . Coupled Geomechanical and Reactive Geochemical Model for Fluid and Heat Flow: Application for Enhanced Geothermal Reservoir. SPE Reservoir Characterization and Simulation Conference and Exhibition, Society of Petroleum Engineers, 2013

[39] Zeng Y C, Wu N Y, Su Z, . Numerical Simulation of Electricity Generation Potential from Fractured Granite Reservoir through a Single Horizontal Well at Yangbajing Geothermal Field. Energy, 2014, 65: 472-487.