We present a pore scale model capable of simulating fluid/fluid reactive transport on images of porous media from first principles. We use a streamline-based particle tracking method for simulating flow and transport, while for reaction to occur, both reactants must be within a diffusive distance of each other during a time-step. We assign a probability of reaction (Pr), as a function of the reaction rate constant (kr) and the diffusion length. Firstly, we validate our model for reaction against analytical solutions for the bimolecular reaction (A+B→C) in a free fluid. Then, we simulate transport and reaction in a beadpack to validate the model through predicting the fluid/fluid reaction experimental results provided by Gramling et al. (2002). Our model accurately predicts the experimental data, as it takes into account the degree of incomplete mixing present at the sub-pore (image voxel) level, in contrast to advection-dispersion-reaction equation (ADRE) model that over-predicts pore scale mixing. Finally, we show how our model can predict dynamic changes in the reaction rate accurately accounting for the local geometry, topology and flow field at the pore scale. We demonstrate the substantial difference between the predicted early-time reaction rate in comparison to the ADRE model.

We perform direct numerical simulation using a pore-scale fluid-fluid reactive transport model (Alhashmi et al. in J Contam Hydrol 179:171-181, 2015. doi:10.1016/j. jconhyd.2015 to investigate the impact of pore structure heterogeneity on the effective reaction rate in different porous media. We simulate flow, transport, and reaction in three pore-scale images: a beadpack, Bentheimer sandstone, and Doddington sandstone for a range of transport and reaction conditions. We compute the reaction rate, velocity distributions, and dispersion coefficient comparing them with the results for non-reactive transport. The rate of reaction is a subtle combination of the amount of mixing and spreading that cannot be predicted from the non-reactive dispersion coefficient. We demonstrate how the flow field heterogeneity affects the effective reaction rate. Dependent on the intrinsic flow field heterogeneity characteristic and the flow rate, reaction may: (a) occur throughout the zones where both resident and injected particles exist (for low Péclet number and a homogeneous flow field), (b) preferentially occur at the trailing edge of the plume (for high Péclet number and a homogeneous flow field), or (c) be disfavored in slow-moving zones (for high Péclet number and a heterogeneous flow field).