Pore-scale modeling of acid etching in a carbonate fracture

Acid fracturing technique has been widely used in the oil and gas industry for improving the carbonate reservoir permeability. In recent years this chemical stimulation technique is borrowed from the oil and gas industry, employed in the enhanced geothermal systems at Groß Schönebeck, Germany (Zimmermann et al., 2010), and at Soultz-sous-Forêts, France (Portier et. al., 2009). In principle, acid fracturing technique utilizes strong acids that react with acid-soluble rock matrix to non-uniformly etch the fracture surfaces. The permeability-enhancing effect depends upon the degree of surface irregularity after pore-scale acidizing which is affected by the compositional heterogeneity of the reacting rock matirx, fracture aperture heterogeneity, and flow and transport heterogeneity. In order to have an insight into these impacts on the acid etching process with the final goal of determining optimum operating conditions (e.g., acid type and acid injection rate), pore-scale acid-fracturing model is needed. The core components of the pore-scale acid-fracturing model consist in tracking the motion of the fluid-matrix boundary surface induced by acid etching. To date, a number of front tracking approaches (e.g., local remeshing technique, embedded boundary method, immersed boundary method, and level-set method) have been proposed by many researchers in order for moving boundary problems. Each approach has its pros and cons. In this work, we propose employing the phase-field approach as an alternative to the existing front tracking approaches to capture the physically sharp concentration discontinuities across the liquid-solid interface. The developed pore-scale acid-fracturing model includes the Stokes-Brinkmann equations for fluid flow in the fracture-matrix system, the multi-component reactive transport equation for transport of solute species in the rough-walled fracture, and the phase-field equation for the reaction-driven motion of the fluid-matrix boundary surface. The simulation results show that the developed pore-scale acid-fracturing model enables to track recession of carbonate fracture surface by acid etching and to capture the solute concentration jump (w.r.t., Ca2+, H+, and HCO3−) across the solid-liquid interface. Reference Zimmermann, G., Moeck, I. and Blöcher, G., 2010. Cyclic waterfrac stimulation to develop an enhanced geothermal system (EGS) — conceptual design and experimental results. Geothermics, 39(1), pp.59-69. Portier, S., Vuataz, F.D., Nami, P., Sanjuan, B. and Gérard, A., 2009. Chemical stimulation techniques for geothermal wells: experiments on the three-well EGS system at Soultz-sous-Forêts, France. Geothermics, 38(4), pp.349-359.