Implementing the Variability of Crystal Surface Reactivity in Reactive Transport Modeling

Current reactive transport model (RTM) uses transport control as the sole arbiter of differences in reactivity. For the simulation of crystal dissolution, a constant reaction rate is assumed for the entire crystal surface as a function of chemical parameters. However, multiple dissolution experiments conirmed the existence of an intrinsic variability of reaction rates, spanning two to three orders of magnitude. Modeling this variance in the dissolution process is vital for predicting the dissolution of minerals in multiple systems. Novel approaches to solve this problem are currently under discussion. Critical applications
include reactions in reservoir rocks, corrosion of materials, or contaminated soils. The goal of this study is to provide an algorithm for multi-rate dissolution of single crystals, to discuss its software implementation, and to present case studies illustrating the diference between the single rate and multi-rate dissolution models. This improved model approach is applied to a set of test cases in order to illustrate the diference between the new model
and the standard approach. First, a Kossel crystal is utilized to illustrate the existence of critical rate modes of crystal faces, edges, and corners. A second system exempliies the efect of multiple rate modes in a reservoir rock system during calcite cement dissolution in a sandstone. The results suggest that reported variations in average dissolution rates can be explained by the multi-rate model, depending on the geometric conigurations of the crystal surfaces.