Pore-scale dissolution mechanisms in calcite-CO2-brine systems: The impact of non-linear reaction kinetics and coupled ion transport

We simulate two sets of dissolution experiments in which CO2-saturated solutions are injected into calcite formations. We explore the impact of non-linear reaction kinetics and charge-coupled ion transport in systems representing different levels of flow and mineralogical complexity. First, we flow CO2-saturated water and brine through cylindrical channels drilled through solid calcite cores and compare directly with experimental dissolution rates. We find that simulations using a linear saturation model match experimental results much better than the batch-reactor-derived non-linear saturation model. The use of a coupled diffusion model causes only a very small increase in the overall dissolution rate compared to a single diffusion coefficient, due to the increase in transport rates of reaction products, particularly the highly charged Ca2+ ion. We also determine the relative importance of the two calcite dissolution pathways, with H+ and H2CO3, and conclude that the H2CO3 – calcite reaction is by far the more dominant, in contrast with common assumptions in the literature. Then, we compare to the experiments of Menke et al. (2015) in which CO2-saturated brine was injected into a microporous Ketton carbonate, and compare dissolution rates over time. We find that including non-linear saturation behaviour markedly changes the simulated dissolution rate, by up to a factor of 0.7 in the case of the experimentally derived saturation model of Anabaraonye (2017), however neither case matches the experimental result which is several times slower than the simulation. Including the effects of coupled ion transport lead to virtually no change in overall dissolution rate due to the convection dominated behaviour. The model also shows differences in the trend of the dissolution rate over time observed in Menke et al, with an approximately linear relationship with time compared to the experimental square-root dependence on time. We conclude that the geochemical model may need to include other effects such as dissolution inside microporous regions.