Rock and Soil Mechanics ›› 2023, Vol. 44 ›› Issue (7): 2143-2150.doi: 10.16285/j.rsm.2022.6230

Previous Articles     Next Articles

Simulation of seepage and heat transfer in 3D fractured rock mass based on fracture continuum method

LIU Dong-dong1, 2, WEI Li-xin2, XU Guo-yuan1, XIANG Yan-yong3   

  1. 1. School of Civil Engineering and Transportation, South China University of Technology, Guangzhou, Guangdong 510641, China 2. Guangzhou Municipal Engineering Design and Research Institute Co. Ltd., Guangzhou, Guangdong 510600, China 3. School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China
  • Online:2023-07-25 Published:2023-10-08
  • About author:LIU Dong-dong, male, born in 1992, PhD candidate, focusing on multiphysics couplings in porous media.

Abstract: The main challenge of simulating the processes of seepage and heat transfer in fractured rock mass is the heterogeneity of rock mass at various scales. In order to balance the efficiency and accuracy of numerical simulation, the two-dimensional fracture continuum method is extended to a three-dimensional one. Based on the depth-first search algorithm, the effective fractures contributing to the permeability of rock blocks are collected. The block permeability tensor is obtained by considering the contribution of effective fractures and rock matrix. The finite element software COMSOL Multiphysics is integrated with Matlab to generate a three-dimensional fracture continuum model that contains multiple blocks with various permeabilities. The simulation results show that due to the combination between stochastic continuum model and discrete fracture model, the fracture continuum model can avoid the complexity of addressing fracture networks and consider the spatial variability of rock mass permeability, thus both the simulation efficiency and accuracy are balanced. As the order of magnitude for the ratio of rock matrix permeability to fracture permeability ranges from 10−4 to 10−6, the discharge calculation error of the effective fracture network model exceeds 5%.

Key words: 3D fractured rock mass, seepage, heat transfer, depth-first search, fracture continuum method