Abstract. Predicting effective permeabilities of fractured rock masses is a key component of reservoir modelling. This is often realized with the discrete fracture network (DFN) method, where single-phase incompressible fluid flow is modelled in discrete representations of individual fractures in a network. Depending on the overall number of fractures, this can result in significant computational costs. Equivalent continuum models (ECM) provide an alternative approach by subdividing the fracture network into a grid of continuous medium cells, over which hydraulic properties are averaged for fluid flow simulations. While this has the advantage of lower computational costs and the possibility to include matrix properties, choosing the right cell size for discretizing the fracture network into an ECM is crucial to provide accurate flow results and conserve anisotropic flow properties. Whereas several techniques exist to map a fracture network onto a grid of continuum cells, the complexity related to flow in fracture intersections is often ignored. Here, numerical simulations of Stokes-flow in simple fracture intersections are utilized to analyze their effect on permeability. It is demonstrated that intersection lineaments oriented parallel to the principal direction of flow increase permeability in a process termed intersection flow localization (IFL). We propose a new method to generate ECM's that includes this effect with a directional pipe flow parametrization: the fracture-and-pipe model. Our approach is tested by conducting resolution tests with a massively parallelized Darcy-flow solver, capable of representing the full permeability anisotropy for individual grid cells. The results suggest that as long as the cell size is smaller than the minimal fracture length and larger than the maximal hydraulic aperture of the considered fracture network, the resulting effective permeabilities and anisotropies are resolution-independent. Within that range, ECM's are applicable to upscale flow in fracture networks, which reduces computational expenses for numerical permeability predictions of fractured rock masses. Furthermore, incorporating the off-diagonal terms of the individual permeability tensors into numerical simulations results in an improved representation of anisotropy in ECM's that was previously reserved for the DFN method.
Correlations between Geometric Properties and Permeability of 2D Fracture Networks
