Reactive Transport Modeling in Porous Fractured Media: Application to Weathering Processes

Abstract

Complex permeability distribution is frequently encountered in weathered rocks with inherited heterogeneities, such as networks of discontinuities that impact fluid flow and water‐rock interactions. Reactive transport modeling is a powerful tool for analyzing fractured media's chemical weathering dynamics. In this study, two different ways to model the fracture‐matrix interaction were explored. A 1D dual‐porosity simulation of reactive transport is designed to study the global impact of fractures on the weathering progression front, compared with a Single Porosity (S.P.) (unfractured) reference model. Then a 2D discrete model using a Discrete Fracture Matrix (DFM) representation is built to highlight the importance of medium geometry and fracture connectivity on weathering. The relevance of such models was evaluated in the case of weathering processes affecting ultramafic units in New Caledonia. Comparing the average behavior of the dual and S.P. models reveals a similar trend in both global mineralogical evolutions. Yet differences occur when confronting the weathering rates. By extending the weathering over a greater depth, the fractures favor a smoother transition between the saprolite horizon and the oxide ore. Furthermore, a partition of the newly formed phyllosilicates occurs in fractures and matrix: Ni‐phyllosilicates preferentially precipitate in fractures, while Mg‐phyllosilicates remain trapped within the matrix due to differences in residence time and liquid‐to‐solid ratio. The DFM modeling highlights the effect of the connectivity of a complex fracture network on the dissolution‐precipitation processes. Depending on the permeability and geometry, fluid flow is modified, favoring the formation of weathering heterogeneities typically observed on the field.