The hydraulic efficiency of single fractures: correcting the cubic law parameterization for self-affine surface roughness and fracture closure

Abstract

Abstract. Quantifying the hydraulic properties of single fractures is a fundamental requirement to understand fluid flow in fractured reservoirs. For an ideal planar fracture, the effective flow is proportional to the cube of the fracture aperture. In contrast, real fractures are rarely planar, and correcting the cubic law in terms of fracture roughness has therefore been a subject of numerous studies in the past. Several empirical relationships between hydraulic and mechanical aperture have been proposed based on statistical variations of the aperture field. However, often, they exhibit non-unique solutions, attributed to the geometrical variety of naturally occurring fractures. In this study, a non-dimensional fracture roughness quantification scheme is acquired, opposing effective surface area against relative fracture closure. This is used to capture deviations from the cubic law as a function of quantified fracture roughness, here termed hydraulic efficiencies. For that, we combine existing methods to generate synthetic 3-D fracture voxel models. Each fracture consists of two random, 25 cm2 wide self-affine surfaces with prescribed roughness amplitude, scaling exponent, and correlation length, which are separated by varying distances to form fracture configurations that are broadly spread in the newly formed two-parameter space (mean apertures in submillimeter range). First, we performed a percolation analysis on 600 000 synthetic fractures to narrow down the parameter space on which to conduct fluid flow simulations. This revealed that the fractional amount of contact and the percolation probability solely depend on the relative fracture closure. Next, Stokes flow calculations are performed, using a 3-D finite differences code on 6400 fracture models to compute directional permeabilities. The deviations from the cubic law prediction and their statistical variability for equal roughness configurations were quantified. The resulting 2-D solution fields reveal decreasing cubic law accordance down to 1 % for extreme roughness configurations. We show that the non-uniqueness of the results significantly reduces if the correlation length of the aperture field is much smaller than the spatial extent of the fracture. An equation was provided that predicts the average behavior of hydraulic efficiencies and respective fracture permeabilities as a function of their statistical properties. A model to capture fluctuations around that average behavior with respect to their correlation lengths has been proposed. Numerical inaccuracies were quantified with a resolution test, revealing an error of 7 %. By this, we propose a revised parameterization for the permeability of rough single fractures, which takes numerical inaccuracies of the flow calculations into account. We show that this approach is more accurate compared to existing formulations. It can be employed to estimate the permeability of fractures if a measure of fracture roughness is available, and it can readily be incorporated in discrete fracture network modeling approaches.