Numerical Simulation on Non-Darcy Flow in a Single Rock Fracture Domain Inverted by Digital Images

The influence of rock seepage must be considered in geotechnical engineering, and understanding the fluid flow in rock fractures is of great concern in the seepage effect investigation. This study is aimed at developing a model for inversion of rock fracture domains based on digital images and further study of non-Darcy flow. The visualization model of single rock fracture domain is realized by digital images, which is further used in flow numerical simulation. We further discuss the influence of fracture domain geometry on non-Darcy flow. The results show that it is feasible to study non-Darcy flow in rock fracture domains by inversion based on digital images. In addition, as the joint roughness coefficient (JRC) increases or the fracture aperture decreases, distortion of the fluid flow path increases, and the pressure gradient loss caused by the inertial force increases. Both coefficients of the Forchheimer equation decrease with increasing fracture aperture and increase with increasing JRC. Meanwhile, the critical Reynolds number tends to decrease when JRC increases or the fracture aperture decreases, indicating that the fluid tends to non-Darcy flow. This work provides a reference for the study of non-Darcy flow through rock fractures.

Permeability of matrix-fracture systems under mechanical loading – constraints from laboratory experiments and 3-D numerical modelling

The permeability of single fractures is commonly approximated by the cubic law assumption, which is however only valid under the condition of a single phase laminar flow between parallel plates. Departure from cubic law are related to many features like aperture fluctuations due to fracture surface roughness, relative shear displacement, the amount of flow exchange between the matrix and the fracture itself, etc. In order to quantify constitutive relationships among the aforementioned aspects, we have conducted a flow-through experiment with a porous rock sample (Flechtinger sandstone) containing a single macroscopic fracture. Based on this experiment, we obtained range of variations of intrinsic rock parameters, permeability and stress-strain relationships of the combined matrix-fracture system under hydrostatic loading. From the measured deformation of the matrix-fracture system, we derived the evolution in the mechanical aperture of the fracture. In order to quantify the processes behind the laboratory observations, we carried out coupled hydro-mechanical simulations of the matrix-fracture system. Navier–Stokes flow was solved in the 3-dimensional open rough fracture domain, and back-coupled to the Darcy flow and the poroelastic behaviour of the rock matrix. The results demonstrate that the elastic behaviour and the related permeability alteration of the fracture domain could be captured by the numerical simulation. Furthermore, the stress-strain values obtained in the vicinity of the fracture asperities suggest that inelastic deformation develops at low mechanical load. An attempt was made to quantify the inelastic deformation by using the failure envelope obtained by laboratory experiments (whether tensile, shear, compaction, or a combination of those). However, change in permeability observed in the experiments are significantly larger than that in the simulation showing the importance of plastic deformation during opening and closure of the fracture and its impact on the cubic law approximation.