Numerical investigation of the flow characteristics and permeability of 2D irregular columnar jointed rock masses

PurposeThe objective of this paper is to provide a better understanding of the effect of irregular columnar jointed structure on the permeability and flow characteristics of rock masses.Design/methodology/approachAn efficient numerical procedure is proposed to investigate the permeability and fluid flow in columnar jointed rock masses (CJRMs), of which the columnar jointed networks are generated by a modified constrained centroid Voronoi algorithm according to the field statistical results. The fractures are represented explicitly by using the lower-dimensional zero thickness elements. And the modeling scheme is validated by a benchmark test for flow in fractured porous media. The effective permeability and representative elementary volume (REV) size of CJRMs are estimated using finite element method (FEM). The influences of joint density and variation coefficient of columnar joint structure on the permeability of the rock mass are discussed.FindingsThe simulation results indicate that the permeability is scale-dependent and tends to be stable with increase of model size. The hydraulic REV size is determined as 3.5 m for CJRMs in the present study. Moreover, the joint density is a dominant factor affecting the permeability of CJRMs. The average permeability of columnar jointed structures increases linearly with the joint density under the same REV size, while the influence from the coefficient of variation can be neglected.Originality/valueThe present paper investigates the REV size of the CJRMs and the effect of joint parameters on the permeability. The proposed method and the results obtained are useful on understanding the hydraulic characteristic of the irregular CJRMs in engineering projects.