Simulation of Fluid Flow on Fractures and Implications for Reactive Transport Simulations

The diverse numerical simulation techniques employed to predict fluid flow properties of fractures yield differing results which limits their applicability for reactive transport simulations. Basically the fluid flow simulation techniques can be divided in two groups: (i) techniques that yield average fluid flow characteristics and (ii) techniques that produce space-resolved properties. These differences may have substantial impact on the reactive transport simulations but may also depend on the fracture characteristics. For this reason, a sensitivity analysis of the geometrical properties of fractures on the fluid flow properties is conducted and evaluated with respect to their impact on reactive transport modeling. Although employing space-resolved simulation techniques, the results of the tests show average values for permeability and fluid velocity that are comparable to previous studies that used other simulation techniques. Observed fluid flow channeling appears to be related to fracture surfaces matching and anisotropy. However, average flow velocities at potential sites for reactive transport differ up to a factor of five from the average ones for the entire fracture. Furthermore, extreme values at reactive transport sites may differ even more and the flow may be directed against the applied pressure gradient. For studies concerned with simulation of reactive transport, these deviations are crucial and should be explicitly considered in the calculations. Hence space-resolved fluid flow simulations should be employed for the simulation of reactive transport.