A numerical manifold method model for analyzing fully coupled hydro-mechanical processes in porous rock masses with discrete fractures

被引:70
作者
Hu, Mengsu [1 ,2 ]
Rutqvist, Jonny [2 ]
Wang, Yuan [1 ,2 ]
机构
[1] Hohai Univ, Coll Civil & Transportat Engn, Nanjing 210098, Jiangsu, Peoples R China
[2] Lawrence Berkeley Natl Lab, Energy Geosci Div, Berkeley, CA 94720 USA
关键词
Discrete fractures; Full coupling; Direct and indirect couplings; Numerical manifold method; Non-conforming mesh; FLUID-FLOW; GEOLOGICAL MEDIA; DEFORMATION; SIMULATION; BEHAVIOR; STRESS;
D O I
10.1016/j.advwatres.2017.02.007
中图分类号
TV21 [水资源调查与水利规划];
学科分类号
081501 ;
摘要
In this study, a numerical manifold method (NMM) model was developed for fully coupled analysis of hydro-mechanical (HM) processes in porous rock masses with discrete fractures. Using an NMM two cover-mesh system of mathematical and physical covers, fractures are conveniently discretized by dividing the mathematical cover along fracture traces to physical cover, resulting in a discontinuous model on a non-conforming mesh. In this model, discrete fracture deformation (e.g. open and slip) and fracture fluid flow within a permeable and deformable porous rock matrix are rigorously considered. For porous rock, direct pore-volume coupling was modeled based on an energy-work scheme. For mechanical analysis of fractures, a fracture constitutive model for mechanically open states was introduced. For fluid flow in fractures, both along-fracture and normal-to-fracture fluid flow are modeled without introducing additional degrees of freedom. When the mechanical aperture of a fracture is changing, its hydraulic aperture and hydraulic conductivity is updated. At the same time, under the effect of coupled deformation and fluid flow, the contact state may dynamically change, and the corresponding contact constraint is updated each time step. Therefore, indirect coupling is realized under stringent considerations of coupled HM effects and fracture constitutive behavior transfer dynamically. To verify the new model, examples involving deformable porous media containing a single and two sets of fractures were designed, showing good accuracy. Last, the model was applied to analyze coupled HM behavior of fractured porous rock domains with complex fracture networks under effects of loading and injection. Published by Elsevier Ltd.
引用
收藏
页码:111 / 126
页数:16
相关论文
共 34 条
[1]  
An X., 2012, INT J COMPUT METHODS, V9
[2]   General theory of three-dimensional consolidation [J].
Biot, MA .
JOURNAL OF APPLIED PHYSICS, 1941, 12 (02) :155-164
[3]   A numerical model for thermo-hydro-mechanical coupling in fractured rock [J].
Bower, KM ;
Zyvoloski, G .
INTERNATIONAL JOURNAL OF ROCK MECHANICS AND MINING SCIENCES, 1997, 34 (08) :1201-1211
[4]  
Chen GQ, 1998, INT J NUMER METH ENG, V43, P685, DOI 10.1002/(SICI)1097-0207(19981030)43:4<685::AID-NME442>3.0.CO
[5]  
2-7
[6]   Development of three-dimensional numerical manifold method for jointed rock slope stability analysis [J].
He, L. ;
An, X. M. ;
Ma, G. W. ;
Zhao, Z. Y. .
INTERNATIONAL JOURNAL OF ROCK MECHANICS AND MINING SCIENCES, 2013, 64 :22-35
[7]  
Hu M., 2016, ACTA GEOTECHNICA
[8]   A practical model for fluid flow in discrete-fracture porous media by using the numerical manifold method [J].
Hu, Mengsu ;
Rutqvist, Jonny ;
Wang, Yuan .
ADVANCES IN WATER RESOURCES, 2016, 97 :38-51
[9]   On continuous and discontinuous approaches for modeling groundwater flow in heterogeneous media using the Numerical Manifold Method: Model development and comparison [J].
Hu, Mengsu ;
Wang, Yuan ;
Rutqvist, Jonny .
ADVANCES IN WATER RESOURCES, 2015, 80 :17-29
[10]  
Itasca Consulting Group, 2013, 3DEC ADV 3 DIM DIST