Transport and Reactive Flow Modelling Using A Particle Tracking Method Based on Continuous Time Random Walks

Mixing and reactive processes have a large impact on the viability of enhanced oil and gas recovery projects that involve acid stimulation and CO2 injection. To achieve a successful design of the injection schemes an accurate understanding of the interplay between pore structure, flow and reactive transport is necessary. Dependent on transport and reactive conditions, this complex coupling can also be dependent on initial rock heterogeneity across a variety of scales. To address these issues, we devise a new method to study transport and reactive flow in porous media at multiple scales. The transport model is based on an efficient Particle Tracking Method based on Continuous Time Random Walks (CTRW-PTM) on a lattice. Transport is modelled using an algorithm described in Rhodes and Blunt (2006) and Srinivasan et al. (2010); this model is expanded to enable for reactive flow predictions in subsurface rock undergoing a first-order fluid/solid chemical reaction. The reaction-induced alteration in fluid/solid interface is accommodated in the model through changes in porosity and flow field, leading to time dependent transport characteristics in the form of transit time distributions which account for rock heterogeneity change. This also enables the study of concentration profiles at the scale of interest. Firstly, we validate transport model by comparing the probability of molecular displacement (propagators) measured by Nuclear Magnetic Resonance (NMR) with our modelled predictions for concentration profiles. The experimental propagators for three different porous media of increasing complexity, a beadpack, a Bentheimer sandstone and a Portland carbonate, show a good agreement with the model. Next, we capture the time evolution of the propagators distribution in a reactive flow experiment, where hydrochloric acid is injected into a limestone rock. We analyse the time-evolving non-Fickian signatures for the transport during reactive flow and observe an increase in transport heterogeneity at latter times, representing the increase in rock heterogeneity. Evolution of transit time distribution is associated with the evolution of concentration profiles, thus highlighting the impact of initial rock structure on the reactive transport for a range of Pe and Da numbers.