Permeability Prediction in Rocks Experiencing Mineral Precipitation and Dissolution: A Numerical Study

In this study, we focus on the electrical tortuosity-based permeability model k = r(eff)(2)/8F (r(eff) is an effective pore size, and F is the formation factor) and analyze its applicability to rocks experiencing mineral precipitation and dissolution. Two limiting cases of advection-dominated water-rock reactions are simulated, that is, the reaction-limited and transport-limited cases. At the pore scale, the two precipitation/dissolution patterns are simulated with a geometrical model and a phenomenological model. The fluid and electric flows in the rocks are simulated by directly solving the linear Stokes equation and Laplace equation on the representative elementary volume of the samples. The numerical results show that evolutions of k and F differ significantly in the two limiting cases. In general, the reaction-limited precipitation/dissolution would result in a smooth variation of k and F, which can be roughly modeled with a power function of porosity phi with a constant exponent. In contrast, the transport-limited precipitation/dissolution mostly occurs near the pore throats where the fluid velocity is high. This induces a sharp change in k and F despite a minor variation in phi. The commonly used power laws with constant exponents are not able to describe such variations. The results also reveal that the electrical tortuosity-based permeability prediction generally works well for rocks experiencing precipitation/dissolution if r(eff) can be appropriately estimated, for example, with the electrical field normalized pore size Lambda. The associated prediction errors are mainly due to the use of electrical tortuosity, which might be considerably larger than the true hydraulic tortuosity.