Two-scale continuum model for simulation of wormholes in carbonate acidization

• Published in: AIChE journal Vol. 51; no. 12; pp. 3231 - 3248
• Main Authors: Panga, Mohan K. R., Ziauddin, Murtaza, Balakotaiah, Vemuri
• Format: Journal Article
• Language: English
• Published: Hoboken Wiley Subscription Services, Inc., A Wiley Company 01.12.2005; Wiley Subscription Services; American Institute of Chemical Engineers
• Subjects: Applied sciences; carbonate acidizing; Chemical engineering; Dissolution; Exact sciences and technology; Heat and mass transfer. Packings, plates; Particle physics; porous media; reactive dissolution; Reactors; wormholes; Mixing; Porous medium; Forecast model; Modeling; Transport process; Heterogeneity; Pore size; Mass transfer coefficient; Reaction mechanism; Rate constant; Diffusion; carbonate acidizing; wormholes; Schmidt number; Permeability; Dissolution; Dispersion; Mass transfer; Chemical reaction; Numerical simulation; Property structure relationship; reactive dissolution; Kinetics; Porosity; porous media
• ISSN: 0001-1541; 1547-5905
• DOI: 10.1002/aic.10574

A two‐scale continuum model is developed to describe transport and reaction mechanisms in reactive dissolution of a porous medium, and used to study wormhole formation during acid stimulation of carbonate cores. The model accounts for pore level physics by coupling local pore‐scale phenomena to macroscopic variables (Darcy velocity, pressure and reactant cup‐mixing concentration) through structure‐property relationships (permeability‐porosity, average pore size‐porosity, and so on), and the dependence of mass transfer and dispersion coefficients on evolving pore scale variables (average pore size and local Reynolds and Schmidt numbers). The gradients in concentration at the pore level caused by flow, species diffusion and chemical reaction are described using two concentration variables and a local mass‐transfer coefficient. Numerical simulations of the model on a two‐dimensional (2‐D) domain show that the model captures the different types of dissolution patterns observed in the experiments. A qualitative criterion for wormhole formation is developed and it is given by Λ ∼ O(1), where Λ = . Here, keff is the effective volumetric dissolution rate constant, DeT is the transverse dispersion coefficient, and uo is the injection velocity. The model is used to examine the influence of the level of dispersion, the heterogeneities present in the core, reaction kinetics and mass transfer on wormhole formation. The model predictions are favorably compared to laboratory data. © 2005 American Institute of Chemical Engineers AIChE J, 2005