The role of fractures on coupled dissolution and precipitation patterns in carbonate rocks
Introduction
The importance of fractures in the transport of contaminants in groundwater systems is well-known. While the hydraulic conductivity of fractures is usually orders of magnitude larger than that of the host rock matrix, dissolution and precipitation processes can modify significantly the physical and chemical properties of fractured porous media. Qualitative and quantitative description of these coupled processes, and their interaction with both anthropogenic and natural pollutants, is fundamental to management and protection of water resources, as well as to the analysis of geological formations and to processes of interest in petroleum and chemical engineering.
A specific instance in which coupled precipitation and dissolution phenomena occur arises in oil reservoir exploitation. Acid is often injected in an attempt to increase hydraulic conductivity in enhanced oil recovery operations, as well as to reduce the resistance to fluid flow during oil well maintenance procedures. However, some dissolution products can precipitate and ultimate clog the pore space. For example, injection of HF acid into sandstone led to partial pore clogging by precipitated silica [3]. Similarly, Rege and Fogler [13] showed that injection of HCl into carbonate rock results in significant precipitation of ferric hydroxide. Another situation for which coupled precipitation and dissolution phenomena can occur is associated with waste rock deposits from mining operations.
Several theoretical and experimental studies have provided partial understanding of the evolution of hydraulic conductivity caused by precipitation or dissolution processes in single fractures [5], [7], [8], [10], [12] and fractured media [17], [18]. However, only a few studies have analyzed effects of coupled precipitation and dissolution on the temporal evolution of fluid flow and solute transport; moreover, these studies have focused on porous rocks without fractures. Rege and Fogler [13] examined both theoretically and experimentally a system of concurrent dissolution of calcium carbonate and precipitation of ferric hydroxide. They demonstrated that periodic fluctuations in the overall permeability, and the frequency of these fluctuations, are related to the interplay between geochemical reactions and flow rates. Singurindy and Berkowitz [15] investigated the evolution of hydraulic conductivity in uniform carbonate rock due to competition among flow, dissolution of calcium carbonate and precipitation of gypsum (CaSO4 · 2H2O; hereafter, we use the term gypsum). Their experimental study presented a range of injected H+/ ratios and flow rates in which the effective hydraulic conductivity displays oscillatory patterns. This study provided a conceptual picture of carbonate rock evolution, describing a variety of possible system behaviors (e.g., system clogging by precipitation dominance or wormhole formation by dissolution dominance), and a detailed description of the factors controlling the oscillations. A parallel study [16] investigated the process of dedolomitization, which couples dolomite dissolution with calcite precipitation and magnesium liberation into solution. Here too, the interplay among flow, precipitation and dissolution processes led to oscillations in the temporal and spatial evolution of effective hydraulic conductivity and porosity.
In spite of the interest of these problems, few analyses exist and the available literature provides only partial understanding of hydraulic conductivity evolution by coupled flow, precipitation and dissolution in fractured rocks. A model simulation describing dissolution and precipitation of minerals during water-rock interaction in fractured granite is given by Sausse et al. [14]. The permeability calculation of this model is based on description only of the fracture network, with the microcrack (fracture) and matrix permeabilities being neglected. Another study of coupled precipitation and dissolution processes in fractured rocks is based on a multicomponent reactive transport model in fractured dolostone [1]. This 1D model accounts for dedolomitization during diffusion from a fracture towards the rock matrix and during advective flow along a fracture. Together, these processes were seen to cause changes in fracture/matrix volume and porosity.
The main aim of this current experimental study is to investigate the influence of fractures on the evolution of hydraulic conductivity structure, caused by the interplay of flow, dissolution of calcium carbonate, and precipitation of gypsum, in carbonate rocks. Specifically, we use laboratory flow cells to determine effects of fracture orientation, fracture wall roughness, fluid flow rate and chemistry, and coupled dissolution/precipitation reaction mechanisms on overall patterns of hydraulic conductivity and porosity evolution. To separate the relative effects of these parameters, we work with flow experiments in quasi-two-dimensional (2D) rock fractures, in three-dimensional (3D) intact rock cores, and in 3D rock cores containing different fracture configurations. While the flow and geochemical conditions considered here do not necessarily correspond to specific natural situations, the dynamic behaviors we investigate are relevant to, e.g., acid injection for enhanced oil recovery [13] and neutralization of sulfuric acid wastes [4].
Section snippets
Rock samples and fluid input solutions
Calcareous sandstone consisting of 97.5% CaCO3 and 2.5% SiO2, collected from the Rosh-Hanikra coastal area of northern Israel, was used as the host rock. Input solutions of two acids, hydrochloric and sulfuric, were injected into the rock samples as reacting fluid. Coupled calcium carbonate dissolution and gypsum precipitation occurred according to the following reactions [9]:
Calcium carbonate dissolution is a mass-transfer limited process
Experimental results
As outlined above, we first examine flow dynamics in the quasi-2D fracture systems. We then consider intact and fractured 3D core samples.
Conclusions
A range of evolutionary patterns of overall hydraulic conductivity occurs in systems for which fluid flow is accompanied by coupled precipitation and dissolution processes. Fractures present in such systems play a particularly important role. We have carried out a series of laboratory experiments which allow us to delineate, qualitatively and in part quantitatively, the effect of fractures on flow and dissolution/precipitation patterns. In the systems we considered, oscillations in the
Acknowledgement
We thank Dr. Ishai Dror, Simon Emmanuel, Dr. Eugenia Klein and Dr. Eytan Sass for assistance and useful discussions. The US-Israel Binational Science Foundation (Contract No. 1999142) and British Petroleum International are thanked for financial support.
References (17)
- et al.
The genesis of dedolomites: a discussion based on reactive transport modeling
J Hydrol
(1998) - et al.
Gypsum overgrowths passivate calcite to acid attack
J Coll Interf Sci
(1997) - et al.
Experimental study of syntaxial vein growth during lateral fluid flow in transmitted lights: first results
J Struct Geol
(2002) - et al.
Evolution of crack permeability during fluid–rock interaction. Example of the Brezouard granite (Vosges, France)
Tectonophysics
(2001) - et al.
Multicomponent reactive transport in discrete fractures: 1. Controls on reaction front geometry
J Hydrol
(1998) Precipitation of hydrated silica from spent hydrofluoric acid: how much of a problem is it
J Petr Technol
(1986)- et al.
Geochemical engineering in chalk: neutralization of waste sulfuric acid
Ann Soc Geol du Nord
(1997) - et al.
Precipitation and dissolution of reactive solutes in fractures
Water Resour Res
(1998)
Cited by (72)
Phase diagram and permeability evolution for dissolving vertical fractures in a gravity field
2024, Advances in Water ResourcesPermeability evolution in open fractures during precipitation and dissolution - A phase-field study
2023, Advances in Water ResourcesCoupled dissolution-precipitation reactions as a potential method for mitigating contaminant transport in carbonate aquifers
2021, Advances in Water ResourcesScale-independent relationship between permeability and resistivity in mated fractures with natural rough surfaces
2021, GeothermicsCitation Excerpt :Channeling flow has also been revealed through field observations (Ishibashi et al., 2015; Tsang and Neretnieks, 1998). In geothermal systems, the channeling flow controls the heat-exchange performance of the hot rock mass and the fluid–rock chemical interactions of dissolution and precipitation (Hawkins et al., 2018; Neuville et al., 2010; Okoroafor and Horne, 2019; Sausse, 2002; Singurindy and Berkowitz, 2005). Thus, in addition to fracture permeability, changes in the channeling flow with changes in fracture aperture need to be clarified for a better understanding of hydraulic activities in geothermal systems.
Experimental and numerical study of bimolecular reactive transport in a single rough-wall fracture
2021, Journal of Hydrology