Simulation of Reactive Transport in Fractured Porous Media

Abstract

Numerical simulations of reactive transport in fractured porous media require the solution of coupled physical and chemical processes that depend on the fractures. Such coupled processes are described by a system of nonlinear partial differential-algebraic equations, while strong heterogeneities characterise fractures. This paper presents an approach to simulate single-phase flow and non-isothermal reactive transport with mineral dissolution and precipitation in fractured porous media. Our numerical solution strategy is based on two ingredients. First, the model equations consist of coupled partial differential equations for the fluid flow, heat transfer and solute transport and nonlinear algebraic equations representing the chemical reactions. Second, fractures are explicitly represented and treated as lower-dimensional objects. The partial differential equations are discretised using finite-volume methods, and at each time step, we solve a nonlinear system of equations using Newton’s method. With numerical simulations, we illustrate our model’s ability to accurately describe the two-way interaction between coupled multi-physical processes and two- and three-dimensional porous media with intersecting fractures.