Fluid Flow Behavior of Sheared Rough Fractures Subjected to Different Stress State

The hydraulic properties of fractures are greatly affected by the stress. Knowing the fluid flow behavior of fractures is of great importance to underground engineering construction and environmental safety. The main purpose of this paper is to study the fluid flow characteristics of rough fractures under different stress states. First, rough fracture surfaces were generated by using the corrected successive random addition (SRA) algorithm. Then, the sheared fracture models subjected to different stress condition were obtained under the boundary condition of constant normal stiffness (CNS). Finally, the hydraulic characteristics of the three-dimensional rough rock fractures were analyzed by numerically solving the full Navier–Stokes equation. It has been found that (1) the aperture of fractures all obeys the Gaussian distribution. The dilatancy effect is gradually obvious and aperture becomes larger with the increase of shear displacement. (2) When the initial normal stress increases, the contact area of fracture becomes larger and the reverse flow can be observed around the contact area. (3) The relationship between hydraulic gradient and flowrate exhibits nonlinearity which can be described by the Forchheimer’s law. The linear coefficient a and the nonlinear coefficient b gradually decrease with the increase of shear displacement and finally stabilize. The values of a and b are reduced by 1–2 and 1–3 orders of magnitude respectively during the shear. The critical Reynolds number increases with the increase of shear displacement and decrease as the initial normal stress increases.