Three-dimensional control-volume distributed multi-point flux approximation coupled with a lower-dimensional surface fracture model

A novel cell-centred control-volume distributed multi-point flux approximation (CVD-MPFA) finite-volume formulation is presented for discrete fracture-matrix simulations on unstructured grids in three-dimensions (3D). The grid is aligned with fractures and barriers which are then modelled as lower-dimensional surface interfaces located between the matrix cells in the physical domain. The three-dimensional pressure equation is solved in the matrix domain coupled with a two-dimensional (2D) surface pressure equation solved over fracture networks via a novel surface CVD-MPFA formulation. The CVD-MPFA formulation naturally handles fractures with anisotropic permeabilities on unstructured grids. Matrix-fracture fluxes are expressed in terms of matrix and fracture pressures and define the transfer function, which is added to the lower-dimensional flow equation and couples the three-dimensional and surface systems. An additional transmission condition is used between matrix cells adjacent to low permeable fractures to couple the velocity and pressure jump across the fractures. Convergence and accuracy of the lower-dimensional fracture model is assessed for highly anisotropic fractures having a range of apertures and permeability tensors. A transport equation for tracer flow is coupled via the Darcy flux for single and intersecting fractures. The lower-dimensional approximation for intersecting fractures avoids the more restrictive CFL condition corresponding to the equi-dimensional approximation with explicit time discretisation. Lower-dimensional fracture model results are compared with equi-dimensional model results. Fractures and barriers are efficiently modelled by lower-dimensional interfaces which yield comparable results to those of the equi-dimensional model. Pressure continuity is built into the model across highly conductive fractures, leading to reduced local degrees of freedom in the CVD-MPFA approximation. The formulation is applied to geologically complex fracture networks in three-dimensions. The effects of the fracture permeability, aperture and grid resolution are also assessed with respect to convergence and computational cost.