Numerical modeling of porous flow in fractured rock and its applications in geothermal energy extraction

Abstract

Understanding the characteristics of hydraulic fracture, porous flow and heat transfer in fractured rock is critical for geothermal power generation applications, and numerical simulation can provide a powerful approach for systematically and thoroughly investigating these problems. In this paper, we present a fully coupled solid-fluid code using discrete element method (DEM) and lattice Boltzmann method (LBM). The DEM with bonded particles is used to model the deformation and fracture in solid, while the LBM is used to model the fluid flow. The two methods are two-way coupled, i.e., the solid part provides a moving boundary condition and transfers momentum to fluid, while the fluid exerts a dragging force to the solid. Two widely used open source codes, the ESyS_Particle and the OpenLB, are integrated into one code and paralleled with Message Passing Interface (MPI) library. Some preliminary 2D simulations, including particles moving in a fluid and hydraulic fracturing induced by injection of fluid into a borehole, are carried out to validate the integrated code. The preliminary results indicate that the new code is capable of reproducing the basic features of hydraulic fracture and thus offers a promising tool for multiscale simulation of porous flow and heat transfer in fractured rock.

Keywords

discrete element method / lattice Boltzmann method / hydraulic fracturing / geothermal energy extraction / multiscale modelling

References

[1] Abe S, Place D, Mora P. A Parallel Implementation of the Lattice Solid Model for the Simulation of Rock Mechanics and Earthquake Dynamics. Pure Appl. Geophys., 2004, 161: 2265-2277.

[2] Bataille A, Genthon P, Rabinowicz M, . Modeling the Coupling between Free and Forced Convection in a Vertical Permeable Slot: Implications for the Heat Production of an Enhanced Geothermal System. Geothermics, 2006, 35: 654-682.

[3] Blocher M G, Zimmermann G, Moeck I, . 3D Numerical Modeling of Hydrothermal Processes during the Lifetime of a Deep Geothermal Reservoir. Geofluids, 2010, 10: 406-421.

[4] Chen S, Doolen G. Lattice Boltzmann Method for Fluid Flows. Ann. Rev. Fluid Mech., 1998, 30: 329-364.

[5] Cundall P A, Strack O D L. Discrete Numerical Model for Granular Assemblies. Geotechnique, 1979, 29: 47-65.

[6] Gong B, Liang H, Xin S, . Effect of Water Injection on Reservoir Temperature during Power Generation in Oil Fields. Proceedings of 36th Workshop on Geothermal Reservoir Engineering, 2011 Stanford: Stanford University, 2011.

[7] Hunsweck M J, Shen Y, Lew A J. A Finite Element Approach to the Simulation of Hydraulic Fractures with Lag. Int. J. Numer. Anal. Meth. Geomech., 2013, 37: 993-1015.

[8] Mora P, Place D. A Lattice Solid Model for the Nonlinear Dynamics of Earthquakes. Int. J. Mod. Phys., 1993, 4: 1059-1074.

[9] Mora P, Place D. Simulation of the Frictional Stick-Slip Instability. Pure Appl. Geophys., 1994, 143: 61-87.

[10] Pang Z H, Hu S B, Wang J Y. A Roadmap to Geothermal Energy Development in China. Science & Technology Review, 2012, 30(32): 18-24.

[11] Secchi S, Schrefler B A. A Method for 3-D Hydraulic Fracturing Simulation. Int. J. Fract., 2012, 178: 245-258.

[12] Wang Y C. A New Algorithm to Model the Dynamics of 3-D Bonded Rigid Bodies with Rotations. Acta Geotechnica, 2009, 4: 117-127.

[13] Wang Y C, Abe S, Latham S, . Implementation of Particle-Scale Rotation in the 3D Lattice Solid Model. Pure Appl. Geophys., 2006, 163: 1769-1785.

[14] Wang Y C, Alonso-Marroquin F. A New Discrete Element Model: Particle Rotation and Parameter Calibration. Granular Matter, 2009, 11: 331-343.

[15] Wang Y C, Mora P. Elastic Properties of Regular Lattices. J. Mech. Phys. Solids, 2008, 56: 3459-3474.

[16] Wang Y C, Mora P. Xing H L. ESyS_Particle: A New 3-D Discrete Element Model with Single Particle Rotation. Advances in Geocomputing, 2009, 183-228.

[17] Wang Y C, Xue S, Xie J. Li Y G. Discrete Element Method and Its Applications in Earthquake and Rock Fracture Modeling. Imaging, Modeling and Assimilation in Seismology, 2012 Beijing: China High Education Press

[18] Zhang X, Jeffrey R G, Thiercelin M. Mechanics of Fluid-Driven Fracture Growth in Naturally Fractured Reservoirs with Simple Network Geometries. J. Geophys. Res., 2009, 114 B12406

[19] Zhou L, Hou M Z. A New Numerical 3D-Model for Simulation of Hydraulic Fracturing in Consideration of Hydro-Mechanical Coupling Effects. International Journal of Rock Mechanics & Mining Sciences, 2013, 60: 370-380.

[20] Zou Q, He X. On Pressure and Velocity Boundary Conditions for the Lattice Boltzmann BGK Model. Phys. Fluids, 1997, 9: 1591-1598.