Abstract
We have extended reactive flow simulation in pore-network models to include geometric changes in the medium from dissolution effects. These effects include changes in pore volume and reactive surface area, as well as topological changes that open new connections. The computed changes were based upon a mineral map from an X-ray computed tomography image of a sandstone core. We studied the effect of these changes on upscaled (pore-scale to core-scale) reaction rates and compared against the predictions of a continuum model. Specifically, we modeled anorthite and kaolinite reactions under acidic flow conditions during which the anorthite reactions remain far from equilibrium (dissolution only), while the kaolinite reactions can be near-equilibrium. Under dissolution changes, core-scale reaction rates continuously and nonlinearly evolved in time. At higher injection rates, agreement with predictions of the continuum model degraded significantly. For the far-from-equilibrium reaction, our results indicate that the ability to correctly capture the heterogeneity in dissolution changes in the reactive mineral surface area is critical to accurately predict upscaled reaction rates. For the near-equilibrium reaction, the ability to correctly capture the heterogeneity in the saturation state remains critical. Inclusion of a Nernst-Planck term to ensure neutral ionic currents under differential diffusion resulted in at most a 9% correction in upscaled rates.