Numerical Simulation of Non-Matching Rough Fracture Seepage

Abstract

Natural rock fractures often exhibit non-matching characteristics at certain scales, leading to uneven aperture distributions that significantly affect fluid flow. This study investigates the impact of the mismatch between the upper and lower surfaces on the flow through three-dimensional rough fractures. By applying fractal theory, a rough upper surface of the fracture is generated, and different degrees of mismatch are introduced by adding random noise to this surface. This approach enables the construction of a variety of three-dimensional rough fracture flow models. Numerical simulations, which involve directly solving the Navier-Stokes equations, are used to simulate flow through a rough single fracture, assessing the effects of various degrees of mismatch between the surfaces. The study also examines how the inclusion of the matrix alters flow characteristics. The results demonstrate that the Forchheimer equation accurately describes the nonlinear flow behavior in fractures with different degrees of mismatch. The increased mismatch intensifies the uneven distribution of fracture apertures, causing the flow velocity to shift from uniform to discrete and the streamlines to become increasingly curved. The overall tortuosity of the flow path increases and the formation of ‘concave’ and ‘convex’ areas leads to vortex zones, promoting nonlinear seepage. The correlation between both viscous and inertial permeability with the degree of mismatch is negative, whereas the impact of matrix permeability on the flow capacity of the fracture shows a positive correlation with a mismatch.