A Discrete Fracture Network Model for Coupled Variable-Density Flow and Dissolution with Dynamic Fracture Aperture Evolution

Abstract

Fluid flow and mass transfer processes in some fractured aquifers are negligible in the low-permeability rock matrix and occur mainly in the fracture network. In this work, we consider coupled variable-density flow (VDF) and mass transport with dissolution in discrete fracture networks (DFNs). These three processes are ruled by nonlinear and strongly coupled partial differential equations (PDEs) due to the (i) density variation induced by concentration and (ii) fracture aperture evolution induced by dissolution. In this study, we develop an efficient model to solve the resulting system of nonlinear PDEs. The new model leverages the method of lines (MOL) to combine the robust finite volume (FV) method for spatial discretization with a high-order method for temporal discretization. A suitable upwind scheme is used on the fracture network to eliminate spurious oscillations in the advection-dominated case. The time step size and the order of the time integration are adapted during simulations to reduce the computational burden while preserving accuracy. The developed VDF-DFN model is validated by simulating saltwater intrusion and dissolution in a coastal fractured aquifer. The results of the VDF-DFN model, in the case of a dense fracture network, show excellent agreement with the Henry semi-analytical solution for saltwater intrusion and dissolution in a coastal aquifer. The VDF-DFN model is then employed to investigate coupled flow, mass transfer and dissolution for an injection/extraction well pair problem. This test problem enables an exploration of how dissolution influences the evolution of the fracture aperture, considering both constant and variable dissolution rates.