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.