Geoelectrical and hydro-chemical monitoring of karst formation at the laboratory scale
Ensuring sustainable strategies to manage water resources in karst reservoirs requires a better understanding of the mechanisms responsible for conduits formation in the rock mass and the development of detection methods for these hydrological and geochemical processes. In this study, we monitored the electrical conductivity of two limestone core samples during controlled dissolution experiments. We interpret the results with a physics-based model describing the porous medium as effective structural parameters that are tortuosity and constrictivity. We obtain that constrictivity is more affected by calcite dissolution compared to tortuosity. Based on our experimental results and data sets from the literature, we show that the characteristic Johnson length is a valuable structural witness of calcite dissolution impact linking electrical and hydrological properties.
Carbonate rocks occupy a vast area on the Earth’s surface (Chen et al., 2017) and constitute reservoirs for key resources such as groundwater (e.g., Kačaroǧlu, 1999; Bakalowicz, 2005), geothermal energy (e.g., Montanari et al., 2017), or fossil energies (e.g., Burchette, 2012). The study of carbonate rocks is also very active because of their potential to serve for carbon dioxide (CO2) geologic sequestration (e.g., Luquot and Gouze, 2009; Cherubini et al., 2019). Carbonate rocks are also well known to be linked to coastal or agricultural issues with contamination, erosion, and landsliding, to civil engineering with risks of cavity presence and collapsing, but also as a commonly exploited noble building material (e.g., Drew et al., 2017; Buckerfield et al., 2020). Investigating carbonate reservoirs is a crucial challenge due to the multi-scale heterogeneity of rock properties and their strong chemical reactivity (e.g., Choquette and Pray, 1970; Lønøy, 2006). These features are responsible for specific processes occurring in carbonate rocks over a wide size range (nm to km), such as groundwater flow and ionic transport in a reactive porous medium (e.g., Ford and Williams, 2013).
Dissolution of carbonate samples caused by CO2 or acid solution injection has already been well studied to understand the formations of wormholes and their complicated relationship with transport properties such as permeability and porosity (e.g., Golfier et al., 2002; Noiriel et al., 2004, 2005; Rötting et al., 2015). However, these experiments were generally based on images analysis, which is an accurate technique for laboratory works but cannot be used on the field. Additionally, in a subsurface context, chemical analysis of the pore water can be quite intrusive, providing only restricted and spatially limited information (e.g., Goldscheider et al., 2008). Thus, studying large-scale heterogeneities such as karst environments can benefit from the use of non-invasive tools such as the ones proposed in hydrogeophysics (e.g., Hubbard and Linde, 2011; Binley et al., 2015). In particular, geoelectrical methods are good candidates to detect the emergence of a sinkhole, to identify infiltration area, or to map ghost-rock features (e.g., Liñán Baena et al., 2009; Chalikakis et al., 2011; Meyerhoff et al., 2014; Kaufmann and Romanov, 2016) since they present a high sensitivity to physical and chemical properties of both porous matrix and interstitial fluids (e.g., Glover, 2015). An increasing amount of work has shown the interest and the effectiveness of geoelectrical methods for laboratory or in situ monitoring of pore space description, hydrological processes, and contaminant transport (e.g., Revil et al., 2012; Garing et al., 2014; Jougnot et al., 2018; Ben Moshe et al., 2021; Sun et al., 2021; Koohbor et al., 2022). Studies on the geoelectrical monitoring of calcite precipitation and dissolution processes are emerging, highlighting their interest for non-intrusive characterizations (e.g., Wu et al., 2010; Zhang et al., 2012; Leroy et al., 2017; Cherubini et al., 2019; Niu and Zhang, 2019; Saneiyan et al., 2019, 2021; Izumoto et al., 2020, 2022; Rembert et al., 2022a). Nevertheless, geophysical methods are indirect, and thus require appropriate models to give a quantitative interpretation (e.g., Binley and Kemna, 2005; Kemna et al., 2012).
The present work attempts to give some answers on how the electrical signal is impacted by conduits formation in limestone due to calcite dissolution, and how the electrical properties can be related to evolving structural parameters. To this end, acid injections are conducted on two core samples of the same limestone at atmospheric conditions and under two different flow rates. The electrical, chemical, and hydrodynamic properties are recorded during the experiments, and the samples are characterized with laboratory and images methods before and after experiments. Evolution of all these properties coupled with the initial conditions is analyzed in order to observe in which proportion they are responsible for the changes induced by the acid flowing through the rock. Then, we focus on the electrical conductivity interpretation and apply a physics-based model that conceptualizes the porous medium as sub-parallel tortuous capillaries which follow a pseudo-fractal size distribution and present sinusoidal variations of their aperture. Finally, we compare the electrical conductivity results of this study with experimental and numerical data sets from the literature.
2.1 Samples properties and petrophysical characterization
The two core samples of this study are part of a wider published data set (Leger and Luquot, 2021; Leger et al., 2022a). They were cored into an Oxfordian crinoid limestone from the Euville quarry in Nancy, France. The two samples are named E04 and E05. The cores have a length of 32 and 31 mm, respectively, and a diameter of 18 mm. They are surrounded with epoxy resin and PVC pipe for a total diameter of 25 mm (i.e., 1 in.).
X-ray diffraction (XRD) on a Bruker D8 Discover and X-ray fluorescence (XRF) measurements show that these rock samples are fully composed of calcium carbonate (CaCO3). To characterize the cores samples' structure, petrophysical measurements are done by non-destructive laboratory methods: porosity ϕ (–), permeability k (m2), and electrical conductivity σ (S m−1) on saturated samples at different salinities to determine the formation factor F (–) and the surface conductivity σs (S m−1), assuming
where σw (S m−1) is the saturating brine conductivity. Note that the surface conductivity can usually be neglected in carbonate material (e.g., Cherubini et al., 2019). The formation factor is an important parameter that can be related to microstructural properties through appropriate petrophysical relationships (e.g., Rembert et al., 2020). Among others, it can be related to the porosity of the medium through Archie (1942) as
where mc (–) is the cementation exponent. Many studies use this relationship to find a unique value of mc to characterize their rock type. From this, for typical carbonate rocks, the cementation exponent is comprised between 2 and 3 (e.g., Lucia, 1983). However, for highly heterogeneous rocks and samples submitted to dissolution, it has been found that Archie's equation is not suitable (e.g., Garing et al., 2014; Niu and Zhang, 2019; Revil et al., 2014). Nevertheless, it has been shown that the formation factor is an interesting proxy to quantify the pore space geometrical properties, such as the tortuosity. Clennell (1997) proposes an extensive review of seminal works on tortuosity, among which the model of Winsauer et al. (1952) that describes the electrical tortuosity τe (–) as
Recent developments in petrophysics tend to develop physics-based models that rely on more accurate conceptualization of the porous medium geometries (e.g., Rembert et al., 2020; Soldi et al., 2017; Vinogradov et al., 2021).
In addition, X-ray microtomographic (XRMT) images of the two core samples are collected with a pixel size of 12 µm and are analyzed with a homemade software to obtain complementary information about the samples structure. In particular, a pore size distribution is obtained from a probabilistic method displaying the chord lengths.
For further details about rock type, calculation methods and protocols, complete results, and discussion on the samples properties, the reader will find the information in Leger and Luquot (2021), Leger et al. (2022a), and Leger et al. (2022b). After dissolution, in order to estimate the impact of the experiments on the core properties, the same laboratory and imaging characterization as before the experiments are conducted on the samples.
2.2 Experimental setup and monitoring
The experimental protocol followed here is the same as the one used in Leger et al. (2022a) and Leger et al. (2022b). Under atmospheric pressure and temperature conditions, acidic solution is injected through the samples. The injected solution is composed of water previously balanced with samples of the same rock type mixed with acetic acid (CH3CO2H) and sodium acetate (CH3CO2Na). The resulting acid concentration is about 10−2 mol L−1 buffered with a pH of 4. Figure 1 displays the homemade experimental device used for the percolation experiments.
The inlet acid solution is stored in a beaker where electrical fluid conductivity and pH are monitored continuously, and is injected through the core sample by a peristaltic pump.
Sensors of absolute and differential pressure continuously record the pressure difference between the sample inlet and outlet. It allows to calculate the evolution of the samples permeability k (m2) during the experiment using Darcy's law,
where Q (m3 s−1) is the flow rate, L (m) is the sample length, μw=0.001 Pa s is the dynamic viscosity of water, A (m2) is the area of fluid injection, and ΔP (Pa) is the differential pressure between the inlet and outlet.
In order to evaluate the impact of the hydraulic conditions, E04 and E05 were flushed with two different flow rates, inducing two different transport conditions. These two experiment conditions are associated to two Péclet numbers Pe (–) defined by
where u (m s−1) is the flow velocity injected in each sample, m2 s−1 is the mean value of diffusion coefficients for ions presented in the fluid at 25 ∘C, and l (m) is the characteristic length, where the value is based on the pore size distribution calculated in Leger and Luquot (2021). The corresponding Péclet numbers are PeE04=5.26 and PeE05=1.97. They correspond to flow rates m3 s−1 and m3 s−1.
The setup is also equipped with a four-wire half bridge array to monitor the water-saturated sample electrical resistance R (Ω) with a CR1000 data logger from Campbell Scientific.
The measured resistance is converted into the conductivity term , where kg (m) is the geometrical coefficient determined from the electrodes shape, size, and distance, using 3D numerical computing on the Electrical Impedance and Diffuse Optical tomography Reconstruction Software (EIDORS). We obtain kg E04=36 mm and kg E05=37 mm.
At the outlet, fluid samplings of 5 mL are continuously collected during day time and punctually during night time. pH, brine conductivity, and cation concentrations are measured for each fluid sample. From calcium (Ca2+) concentration monitoring, we can determine the mass of dissolved calcite through time, and thus obtain the time variations of the rock sample porosity ϕ (–). At time ti (s), we have , where ϕdt (–) is the porosity difference over the time laps (s). By definition, , where Vech=AL (m3) is the core sample volume and Vdt (m3) is the pore volume difference over time laps dt. The latter is calculated as follows:
where (g) is the mass of calcium carbonate over time laps dt, (g m−3) is the density of calcium carbonate, (mol m−3) is the calcium concentration at the sampling time ti, (mol m−3) is the initial calcium concentration, and (g mol−1) is the molecular masses of calcium carbonate.
2.3 Analytical model
2.3.1 Capillary size distribution
To relate the electrical monitoring with the porous medium microstructure, we use the bundle of capillary tubes (BOCT) approach (e.g., Vinogradov et al., 2021). It uses a bundle of parallel tortuous and constrictive capillaries as a conceptualization of the pore structure (e.g., Guarracino et al., 2014; Rembert et al., 2020). The conceptual geometry is represented in Fig. 2.
The BOCT model considers N capillaries that do not intersect and run with the same orientation. The capillaries are confined within a cylindrical model of length L (m) and radius R (m). They are permitted to have different radius r (m), but a single length l. Thus, the tortuosity (–) is an effective property. Each tortuous and constrictive capillary presents a varying radius r(x) (m) defined with the following sinusoidal expression:
where (m) is the average capillary radius, r′ (m) the amplitude of the radius size fluctuation, and λ (m) is the wavelength. The parameter a (–) is the pore radius fluctuation ratio, also called the constriction factor. It is defined by , thus ranging from 0 to 0.5, corresponding to cylindrical pores () and periodically closed pores, respectively.
with the number of capillaries of radius between and . As many geologic materials present a skewed pore size distribution (e.g., Bennion and Griffiths, 1966; Dullien, 1992), we use a pseudo-fractal law to describe the distribution of capillaries throughout the model (Jackson, 2008; Vinogradov et al., 2021),
where is a dimensionless exponent. and (m) are the extreme average pore sizes. They are estimated from the chord length computation using 3D tomography of the samples before and after the dissolution experiments.
Thus, N is calculated for average pore sizes comprised between these two extremes. The normalization coefficient, D=1, since the pore size distribution from the chord lengths estimation is normalized and expressed in percent. The minus sign of −dN implies that the number of pores decreases when the average radius increases (Soldi et al., 2017; Thanh et al., 2019; Yu et al., 2003).
2.3.2 Electrical properties
We follow the approach of Rembert et al. (2020), which consists in the expression of the electrical properties at the pore scale before upscaling them, considering the bundle of capillaries as an equivalent circuit of parallel conductances. This approach is valid for a negligible surface conductivity and leads to
where f (–) is the constrictivity of the porous medium varying from 0 to 1, corresponding to periodically closed capillaries and cylindrical capillaries, respectively. This approach quantifies the effect of pore throats and pore bodies on the medium electrical conductivity.
Rembert et al. (2020) demonstrated that the constriction factor and the tortuosity vary during dissolution process, since it affects the structure of the porous medium. Thus, a(ϕ) and τ(ϕ) are defined following logarithmic laws:
where Pa and Pτ are empirical dimensionless parameters.
2.3.3 Johnson length
The Johnson length Λ (m) represents a characteristic length that corresponds to the dynamically connected pores (Banavar and Schwartz, 1987; Ghanbarian, 2020). Following Johnson et al. (1986), it can be defined from the permeability and the formation factor by
3.1 Hydro-chemical and electrical results from percolation experiments
Figure 3 displays the evolution over time of measured and calculated parameters for the two percolation experiments, with red curves for E04 and blue curves for E05. In detail, there are porosity, permeability, outlet pH values, calcium concentration and core sample conductivity, water conductivity, and normalized formation factor. The initial and final values of the measured parameters obtained from petrophysical measurements are presented in Table 1. They are also represented for the initial and final porosity with error bars, permeability, and normalized formation factor values as diamonds and circles for E04 and E05 in Fig. 3, respectively. One can note that error bars are higher for final porosity values compared to the initial values. This comes from the formation of conduits through the sample due to dissolution, which allows water to leak out and is responsible for the measurement accuracy deterioration.
One can observe that E04 and E05 curves present similar amplitudes and shapes for all the measured and calculated parameters displayed in Fig. 3. However, due to lower flow regime (Leger et al., 2022a), the curves of E05 present retarded or lower variations compared to E04 curves.
The generation of dissolved calcium (Fig. 3d) and the increase of porosity (Fig. 3a), permeability (Fig. 3b), and water conductivity difference (Fig. 3e) indicate that calcite dissolution is the main chemical reaction for both experiments. However, calcium difference shows slow decrease until the end of the experiment, which can also be seen in the pH curves and the water conductivity difference curves (Fig. 3c and e), and indicates that the dissolution regime was initially strong and got weaker over time.
The permeability curves (Fig. 3b) show a small initial increase followed by a stabilization for a few hours, then they increase sharply and reach final values several orders of magnitude higher than the initial values. The initiation of this sharp increase in permeability is interpreted as being due to the breakthrough of the sample by the acid solution, and the corresponding time, defined as the percolation time, is highlighted as a vertical dashed line in Fig. 3 (in red for E04 and blue for E05). The percolation times are 4.3 and 11.9 h for E04 and E05, respectively. For E05, the percolation time is also observable through a slight break in slope on the pH, dissolved calcium, and water conductivity curves.
Before the percolation time, one can observe that E05 permeability curve presents variations attributed to complex effects of grain displacement and transitory pore clogging. After percolation time, both E04 and E05 permeability curves present high-frequency noise corresponding to the limit of detection of the sensors of differential pressure. In order to avoid these instabilities, the curves are filtered. After a comparison with the moving average and Butterworth filters, we choose the Savitsky–Golay filter, since it offers a better adjustment. We observe that the filtered curves are well superimposed on the curves of the measurements, except between 10 and 12 h for the dissolution of E05, because of the transient complex effects. The final permeability of E05 is higher than the last value of the temporal series. This discrepancy is also due to the limit of detection of the sensors of differential pressure. However, we did not take this higher value into account when filtering, to be in accordance with the temporal series.
The sample conductivity increases during both experiments (Fig. 3d). At first glance, this increase could be attributed to the increase of the pore water conductivity (Fig. 3e) instead of calcite dissolution. However, one can observe that the normalized formation factor (Fig. 3f) decreases during the experiment. Thus, the sample conductivity and formation factor variations are linked to calcite dissolution. It can be observed that for sample conductivity and formation factor, the variations are monotonous, with a rapid initial change followed by a gentler slope until the end of the experiment. The initial slopes are correlated to higher calcium concentrations at the beginning of the experiment (Fig. 3d) and are thus due to initial strong dissolution regimes. The sample conductivity and formation factor smooth variations show that they are more impacted by porosity increase than by permeability variations (Fig. 3a and b). Furthermore, the comparison of the results between E04 and E05 shows that with equal water conductivity and close permeability values before the breakthrough, higher porosity induces higher sample conductivity. Compared to the petrophysical measurements, the normalized formation factor curves show good agreement with the measured amplitude difference between the initial and final values, and are represented as the final normalized formation factor (Fig. 3f).
3.2 Applying the pseudo-fractal law to the chord length distribution
Leger et al. (2022a) have published the chord length distribution of E04 and E05 samples before and after percolation in X, Y, and Z directions. Since there are some slight variations depending on the direction, we apply the pseudo-fractal law describing the distribution of capillaries using Eq. (9) on the mean chord length distribution. The adjusted parameter mj is determined using a Monte Carlo optimization scheme. From the results displayed in Fig. 4, the pseudo-fractal law shows a good agreement with the chord length distributions for both samples.
In order to attest of the accuracy of the pseudo-fractal law to describe the chord length distribution, we computed the corresponding error ϵ<0.5 % for each value of the exponent mj. The error ϵ is defined as follows:
where Nd, , and refer to the number of data, the number of capillaries from the data, and the number of capillaries from the model, respectively. The best adjustment is found for values of mj of the same order of magnitude (see legends of Fig. 4), reflecting that E04 and E05 are samples of the same rock type. However, since mj can vary over large ranges, the slight differences obtained for the four different data sets (although accurately determined, see Appendix A), cannot be interpreted as impact of the dissolution on this parameter describing the pore size distribution, while percolation has a huge impact on the rock properties as presented above (see Fig. 3 and Table 1). Thus, the evolution of mj is not a sufficient indicator of the effect of dissolution on the sample.
Moreover, Leger et al. (2022a) have also determined the pore size distribution by centrifugation. Since the size intervals for the determination of the pore number are not regular, we could not use the pseudo-fractal law on these data. Nevertheless, the characterization of pore size distribution by centrifugation has shown that between 46 and 70 % of the pores have a pore size comprised between 0.5 and 12 µm, while the minimal resolution of the chord length is of 12 µm. In addition, a very porous area of the sample only composed of small pores can be interpreted by the tomography software as a single pore over the entire pixel. This means that the chord length distribution does not have the best precision for the smallest pores distribution. However, Leger and Luquot (2021) noticed the pore size distribution from the two methods are similar and consistent, leading to the conclusion that large pores are more likely to control rock parameters on core scale. In addition, this importance of macropores is even more of interest in reactive transport, as it can be seen in fractured media (e.g., Noiriel, 2015; Garcia-Rios et al., 2017).
3.3 Influence of rock structure on electrical signature
We have shown in the previous section that the pore size distributions of both samples were successfully describable using pseudo-fractal laws. Thus, we can apply the relationship relating the formation factor and the porosity developed in Eq. (11). We can see in Fig. 5 the results of the Monte Carlo inversions to adjust the constrictivity f, defined from Eq. (10), and the tortuosity τ. To reflect the structural changes in the porous medium, we consider the parameters as distributions a(ϕ) and τ(ϕ), expressed in Eqs. (12) and (13). We observe that the evolution of f and τ with the porosity is physically plausible as the increase of porosity due to calcite dissolution reflects less complex medium geometries (i.e., less constrictive and less tortuous), described with higher values of f and lower values of τ. The accuracy of the model is evaluated using the error computation of Eq. (15). For both samples it is below 10 %.
The constrictivity varies over different ranges of values between E04 and E05 (Fig. 5c and d). For E05, f is lower by half (i.e., more constrictive) compared to E04, which is in agreement with E04 higher porosity values. In contrast, the evolution with porosity follows roughly the same law for these two samples. This indicates that we are working on two samples of the same rock type with similar pore shapes.
The tortuosity varies little during the experiment for both samples, although it is more pronounced for E05 (Fig. 5e and f). A surprising result is that the tortuosity values are higher for E04 than for E05, while E05 is less porous. This result is in contradiction with the values of initial and final electrical tortuosity τe displayed in Table 1, showing higher values for both samples and higher values for E05 compared to E04. This discrepancy has already been observed in Niu and Zhang (2019), and can be explained with the definition of τe, given in Eq. (3), which combines the impact of the tortuosity and the constrictivity, and is thus artificially higher. Moreover, the main conduits generated by dissolution can be visualized from tomographic images. Figure 6 shows the phase connecting the inlet to the outlet for the two samples, according to two different image treatments. Two segmentations based on a region growing method were applied to isolate two different void phases. The red phase represents the phase connecting the largest pores in the samples from the inlet to the outlet, as already explained and used in Leger et al. (2022a). The blue phase only shows the darker pixels of the images, corresponding to the center of the conduit used by the fluid to cross the samples. In comparison, E05 has a straighter path than E04. Thus, the model parameters adjustment reproduces well the tomography observations and the experimental data.
However, for the lowest porosity data in Fig. 5, the formation factor shows high values, highlighted in purple, that do not fit with the gentle slope of the model. As mentioned earlier, these points correspond to an initial strong dissolution regime that is out of the validity domain of the model (Guarracino et al., 2014; Rembert et al., 2020). Moreover, the points in green in Fig. 5 highlight the data acquired after the percolation, where the samples cannot be considered anymore as representative elementary volume (REV) and where the model is not valid either. Nevertheless, for E04, the acid percolation does not show radical changes in the formation factor and porosity compared to the smooth tendency that follows the initial strong dissolution regime, whereas for E05, the formation factor falls out of range after percolation. These different behaviors are due to the imposed hydrodynamic regimes and their consequences on the dissolution profile. Indeed, the red phase in Fig. 6 shows that the new paths created in E04 are concentrated in a reduced volume of the sample, while for E05, the red phase is extensively spread in the upper part of the sample. This portion of the sample, close to the inlet, was probably dissolved at an early stage of the experiment, therefore creating no apparent effect on the permeability or formation factor before the sample percolation.
3.4 Influence of hydrodynamic properties on rock structure and electrical signature
3.4.1 Impact of the experimental conditions and permeability changes
In this section, we are interested in highlighting the link between the evolution of hydrodynamic and electrical properties. The permeability is plotted as a function of the formation factor and porosity in Fig. 7. We can first note that the permeability varies over the same ranges of values for both samples, while E04 has a larger porosity (and lower formation factor) compared to E05. As discussed in the literature (e.g., Guarracino et al., 2014; Soldi et al., 2017; Rembert et al., 2020), the pore throats play a primary role on the hydrodynamic behavior of porous media. According to the model used in this study, the studied samples follow the same distribution for constrictivity (Fig. 5), which is consistent with the similarity in permeability.
It can be observed that before the breakthrough, the permeability increases little, while the formation factor and the porosity vary strongly. Then, after percolation, the trends are reversed, the permeability increases strongly compared to the lower variation of the porosity and the formation factor. This behavior change is very clear for both samples and indicates that monitoring electrical properties allows us to be sensitive to the impact of dissolution on the porous medium long before the sample is percolated.
In order to locate our experiments in the scientific landscape, we plotted our data with the well known power law relating permeability and porosity k∼ϕn (Kozeny, 1927; Carman, 1937), where n (–) is an empirical exponent. This relationship has been widely used in literature, sometimes by trying to determine the n exponent (Martys et al., 1994), sometimes by taking into account the critical pore size (Gueguen and Dienes, 1989; Thompson et al., 1987), the specific surface with the Kozeny factor (Paterson, 1983; Walsh and Brace, 1984; Fabricius et al., 2007; Weger et al., 2009), or the grain diameter in a non-fractal (Kozeny, 1927; Carman, 1937, 1948, 1956; Bourbie et al., 1987) or in a fractal dimension (Pape et al., 1999; Costa, 2006). The classic power law has been tested here, and the n values obtained are a bit high, but are still consistent with those proposed in the literature on carbonate rocks, known for their natural heterogeneity and also the one created by dissolution (see Leger et al., 2022b, for a synthetic review).
3.4.2 The characteristic pore size defined through the Johnson length
The comparison of the permeability of E04 and E05 reveals that they are samples with similar initial hydrodynamic properties related to their identical rock type (Table 1). However, their difference of porosity leads to different ranges for the formation factor, and the different imposed hydrodynamic regimes lead to higher final permeability for E05. Moreover, the use of a normalized formation factor (Fig. 3) enables to compare trends but does not allow to quantitatively interpret the electrical monitoring in term of rock structure. Thus, as recommended by Niu and Zhang (2019), we use the Johnson length as an effective pore size (Fig. 8), that enables to compare the behavior of both samples E04 and E05 regarding other data sets from the literature (Niu and Zhang, 2019; Vialle et al., 2014).
The series from Niu and Zhang (2019) comes from calcite dissolution simulation of a digital image of a carbonate mudstone (Wellington Formation, Kansas, USA), that has an initial porosity of 0.13. The fluid transport is assumed as advection-dominated (Pe>1), and under this transport condition, referring to Pereira Nunes et al. (2016), they impose a transport‐limited case, related to wormholing, in which the reaction at the solid‐-liquid interface is limited by the diffusion of reactants to and from solid surfaces (PeDa>1, where is the so-called Damköhler number and kr is the local reaction rate).
The data points from Vialle et al. (2014) correspond to bioclastic limestone (Estaillades limestone, France), composed of more than 99 % calcite and of initial porosity of 0.29. The values for Pe ranges from 18.5 (for the smaller grains) to 37 (for the larger ones) and Da=3486. Meaning that for this data set, the dissolution is also transport-controlled and leads to dominant wormhole formation.
One can observe that all the series present the same behavior, but with different trends and on different ranges. These differences do not seem related to initial porosity value or hydrodynamic regime. Indeed, for different Péclet numbers, there are close values and trends for E04 and E05 compared to the other series. In addition, while the carbonate mudstone from Niu and Zhang (2019) presents similar initial porosity value with E05, they do present contrasted Johnson lengths. Thus, the Johnson length is related to the rock type and highlights the difference in the pore structure among the samples.
This study investigates the impact of calcite dissolution on conductivity monitoring through percolation experiments of two crinoidal limestone core samples in different hydrodynamic conditions. We thus demonstrate that the lower the Peclet condition applied to the samples is, the lower the dissolution rate is.
Then, we express the electrical conductivity and the formation factor of the porous medium in terms of effective petrophysical parameters: tortuosity and constrictivity. In this representation, the pore space is conceptualized as a fractal cumulative distribution of tortuous capillaries with a sinusoidal variation of their radius (i.e., periodical pore throats). We verify the assumption of a fractal pore size distribution using the chord length distribution computed from XRMT images. The confrontation of the model with the formation factor monitoring for both samples shows that the model can reproduce structural changes linked to calcite dissolution by setting logarithmic dependence of constriction factor and tortuosity to porosity. These empirical relations allow the model to accurately fit the data within the range of validity of the model, defined as a steady reaction rate obtained after first acid invasion through the sample and before the acid breakthrough. We observe that dissolution leads to less constrictive rather than less tortuous sample, thus designating constrictivity as a property of the pore structure highly impacted by chemical processes.
By means of an upscaling procedure, we also link the formation factor to permeability through the Johnson length. We show on the experimental data of this study and other data sets (experimental and numerical) from the literature that the Johnson length is a relevant parameter that allows comparison of chemical reaction in comparable Peclet–Damköhler conditions on different rock types (crinoidal limestone, carbonate mudstone, and bioclastic limestone).
We believe that the present study contributes to a better understanding of the links between the electrical conductivity measurement and the evolution of the hydrodynamic and structural properties of the porous medium subjected to reactive transport.
The accuracy of the pseudo-fractal computation for the determination of the pore size distribution is controlled through the error ϵ defined in Eq. (15). Since dissolution does not generate notable variations of the pseudo-fractal exponent mj, we confront the computed error to a range of value mj in Fig. A1. We obtain V-shaped curves for each data set, showing that mj is well determined for each sample, before and after percolation experiments.
Samples of the same rock type are available upon request.
ML and LL designed the experiments and carried them out. FR developed the model code and performed the simulations. FR and ML prepared the manuscript with contributions from all the co-authors. DJ and LL supervised the work and provided the resources.
The contact author has declared that none of the authors has any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
We acknowledge the MRI platform member of the national infrastructure France-BioImaging, the labex CEMEB, the labex VOLTAIRE, and the CIFRE PhD fellowship provided by Voxaya and ANRT. We also acknowledge the petrophysics platform members of Geosciences Montpellier Laboratory.
This research has been supported by the Centre National de la Recherche Scientifique (grant no. EC2CO StarTrek), the Agence Nationale de la Recherche (grant nos. ANR-10-LABX-100-01, ANR-10-INBS-04, ANR-10-LABX-0004, and ANR-10-LABX-0020), and the WaterWorks2017 project of the Water Joint Programming Initiative (grant nos. UrbanWat and Maradentro).
This paper was edited by Alberto Guadagnini and reviewed by two anonymous referees.
Ben Moshe, S., Kessouri, P., Erlich, D., and Furman, A.: Geophysically based analysis of breakthrough curves and ion exchange processes in soil, Hydrol. Earth Syst. Sci., 25, 3041–3052, https://doi.org/10.5194/hess-25-3041-2021, 2021. a
Binley, A. and Kemna, A.: DC resistivity and induced polarization methods, Springer, edited by: Rubin, Y. and Hubbard, S. S., Hydrogeophysics, 11, 19–41, https://doi.org/10.1016/B978-0-444-53802-4.00188-3, 2005. a
Binley, A., Hubbard, A. S., Huisman, J. A., Revil, A., Robinson, D. A., Singha, K., and Slater, L. D.: The emergence of hydrogeophysics for improved understanding of subsurface processes over multiple scales, Water Resour. Res., 51, 3837–3866, https://doi.org/10.1002/2015WR017016, 2015. a
Bourbie, T., Coussy, O., Zinszner, B., and Junger, M. C.: Acoustics of Porous Media, Houston, TX, gulf publishing company edn., J. Acoust. Soc. Am., 91, 3080, https://doi.org/10.1121/1.402899, 1987. a
Buckerfield, S. J., Quilliam, R. S., Bussiere, L., Waldron, S., Naylor, L. A., Li, S., and Oliver, D. M.: Chronic urban hotspots and agricultural drainage drive microbial pollution of karst water resources in rural developing regions, Sci. Total Environ., 744, 140898, https://doi.org/10.1016/j.scitotenv.2020.140898, 2020. a
Burchette, T. P.: Carbonate rocks and petroleum reservoirs: a geological perspective from the industry, Geological Society, London, Special Publications, 370, 17–37, https://doi.org/10.1144/SP370.14, 2012. a
Carman, P. C.: Flow of gases through porous media, Butterworths Scientific Publications, London, 182 pp., 1956. a
Chalikakis, K., Guérin, R., Valois, R., and Bosch, F.: Contribution of geophysical methods to karst-system exploration: an overview, Hydogeol. J., 19, 1169–1180, https://doi.org/10.1007/s10040-011-0746-x, 2011. a
Chen, Z., Goldscheider, N., Auler, A. S., and Bakalowicz, M.: The World karst aquifer mapping project: concept, mapping procedure and map of Europe, Hydrogeol. J., 25, 771–785, https://doi.org/10.1007/s10040-016-1519-3, 2017. a
Cherubini, A., Garcia, B., Cerepi, A., and Revil, A.: Influence of CO2 on the electrical conductivity and streaming potential of carbonate rocks, J. Geophys. Res.-Sol. Ea., 124, 10056–10073, https://doi.org/10.1029/2018JB017057, 2019. a, b, c
Choquette, P. W. and Pray, L. C.: Geologic nomenclature and classification of porosity in sedimentary carbonates, AAPG Bulletin, 54, 207–250, https://doi.org/10.1306/5D25C98B-16C1-11D7-8645000102C1865D, 1970. a
Costa, A.: Permeability-porosity relationship: A reexamination of the Kozeny-Carman equation based on a fractal pore-space geometry assumption, Geophys. Res. Lett., 33, L02318, https://doi.org/10.1029/2005GL025134, 2006. a
Drew, D., Lamoreaux, P. E., Coxon, C., Wess, J. W., Slattery, L. D., Bosch, A. P., and Hötzl, H.: Karst hydrogeology and human activities: impacts, consequences and implications: IAH International Contributions To Hydrogeology 20, 5th edn., 338, https://doi.org/10.1201/9780203749692, 2017. a
Dullien, F. A. L.: 1 – Pore structure, in: Porous media (2nd edition), Academic Press, San Diego, 1992, 5–115, https://doi.org/10.1016/B978-0-12-223651-8.50007-9, 1992. a
Garcia-Rios, M., Luquot, L., Soler, J. M., and Cama, J.: The role of mineral heterogeneity on the hydrogeochemical response of two fractured reservoir rocks in contact with dissolved CO2, Appl. Geochem., 84, 202–217, https://doi.org/10.1016/j.apgeochem.2017.06.008, 2017. a
Ghanbarian, B.: Applications of critical path analysis to uniform grain packings with narrow conductance distributions: I. Single-phase permeability, Adv. Water Res., 137, 103529, https://doi.org/10.1016/j.advwatres.2020.103529, 2020. a
Glover, P. W. J.: Geophysical properties of the near surface Earth: electrical properties, Treatise on Geophysics, 11, 89–137, https://doi.org/10.1016/B978-0-444-53802-4.00189-5, 2015. a
Golfier, F., Zarcone, C., Bazin, B., Lenormand, R., Lasseux, D., and Quintard, M.: On the ability of a Darcy-scale model to capture wormhole formation during the dissolution of a porous medium, J. Fluid Mech., 457, 213–254, https://doi.org/10.1017/S0022112002007735, 2002. a
Guarracino, L., Rötting, T., and Carrera, J.: A fractal model to describe the evolution of multiphase flow properties during mineral dissolution, Adv. Water Res., 67, 78–86, https://doi.org/10.1029/WR024i004p00566, 2014. a, b, c
Hubbard, S. S. and Linde, N.: 2.15 – Hydrogeophysics, in: Treatise on Water Science, edited by: Wilderer, P., Elsevier, Oxford, 2, 401–434, https://doi.org/10.1016/B978-0-444-53199-5.00043-9, 2011. a
Izumoto, S., Huisman, J. A., Wu, Y., and Vereecken, H.: Effect of solute concentration on the spectral induced polarization response of calcite precipitation, Geophys. J. Int., 220, 1187–1196, https://doi.org/10.1093/gji/ggz515, 2020. a
Izumoto, S., Huisman, J. A., Zimmermann, E., Heyman, J., Gomez, F., Tabuteau, H., Laniel, R., Vereecken, H., Méheust, Y., and Le Borgne, T.: Pore-scale mechanisms for spectral induced polarization of calcite precipitation inferred from geo-electrical millifluidics, Environ. Sci. Technol., 56, 4998–5008, https://doi.org/10.1021/acs.est.1c07742, 2022. a
Jackson, M. D.: Characterization of multiphase electrokinetic coupling using a bundle of capillary tubes model, J. Geophys. Res.-Sol. Ea., 113, B04201, https://doi.org/10.1029/2007JB005490, 2008. a, b
Johnson, D., Koplik, J., and Schwartz, L.: New pore-size parameter characterizing transport in porous media, Phys. Rev. Lett., 57, 2564–2567, 1986. a
Jougnot, D., Jiménez-Martínez, J., Legendre, R., Le Borgne, T., Méheust, Y., and Linde, N.: Impact of small-scale saline tracer heterogeneity on electrical resistivity monitoring in fully and partially saturated porous media: Insights from geoelectrical milli-fluidic experiments, Adv. Water Res., 113, 295–309, https://doi.org/10.1016/j.advwatres.2018.01.014, 2018. a
Kaufmann, G. and Romanov, D.: Structure and evolution of collapse sinkholes: Combined interpretation from physico-chemical modelling and geophysical field work, J. Hydrol., 540, 688–698, https://doi.org/10.1016/j.jhydrol.2016.06.050, 2016. a
Kemna, A., Binley, A., Cassiani, G., Niederleithinger, E., Revil, A., Slater, L., Williams, K. H., Orozco, A. F., Haegel, F.-H., Hördt, A., Kruschwitz, S., Leroux, V., Titov, K., and Zimmermann, E.: An overview of the spectral induced polarization method for near-surface applications, Near Surf. Geophys., 10, 453–468, https://doi.org/10.3997/1873-0604.2012027, 2012. a
Koohbor, B., Deparis, J., Leroy, P., Ataie-Ashtiani, B., Davarzani, H., and Colombano, S.: DNAPL flow and complex electrical resistivity evolution in saturated porous media: A coupled numerical simulation, J. Contam. Hydrol., 248, 104003, https://doi.org/10.1016/j.jconhyd.2022.104003, 2022. a
Kozeny, J.: Über kapillare Leitung des Wassers im Boden-Aufstieg, Versickerung und Anwendung auf die Bewasserung, Sitzungsberichte der Akademie der Wissenschaften Wien, Mathematisch Naturwissenschaftliche Abteilung, 136, 271–306, http://www.zobodat.at/pdf/SBAWW_136_2a_0271-0306.pdf (last access: 16 January 2023), 1927. a, b, c
Leger, M. and Luquot, L.: Importance of microstructure in carbonate rocks: laboratory and 3D-imaging petrophysical characterization, Appl. Sci., 11, 3784, https://doi.org/10.3390/app11093784, 2021. a, b, c, d
Leger, M., Luquot, L., and Roubinet, D.: Role of mineralogical, structural and hydrodynamic rock properties in conduits formation in three distinct carbonate rock types, Chem. Geol., 607, 121008, https://doi.org/10.1016/j.chemgeo.2022.121008, 2022a. a, b, c, d, e, f, g
Leger, M., Roubinet, D., Jamet, M., and Luquot, L.: Impact of hydro-chemical conditions on structural and hydro-mechanical properties of chalk samples during dissolution experiments, Chem. Geol., 594, 120763, https://doi.org/10.1016/j.chemgeo.2022.120763, 2022b. a, b, c
Leroy, P., Li, S., Revil, A., and Wu, Y.: Modelling the evolution of complex conductivity during calcite precipitation on glass beads, Geophys. J. Int., 209, 123–140, https://doi.org/10.1093/gji/ggx001, 2017. a
Liñán Baena, C., Andreo, B., Mudry, J., and Carrasco Cantos, F.: Groundwater temperature and electrical conductivity as tools to characterize flow patterns in carbonate aquifers: The Sierra de las Nieves karst aquifer, southern Spain, Hydrogeol. J., 17, 843–853, https://doi.org/10.1007/s10040-008-0395-x, 2009. a
Lucia, F.: Petrophysical parameters estimated from visual descriptions of carbonate rocks: a field classification of carbonate pore space, J. Petrol. Technol., 35, 629–637, https://doi.org/10.2118/10073-PA, 1983. a
Luquot, L. and Gouze, P.: Experimental determination of porosity and permeability changes induced by injection of CO2 into carbonate rocks, Chem. Geol., 265, 148–159, https://doi.org/10.1016/j.chemgeo.2009.03.028, 2009. a
Meyerhoff, S. B., Maxwell, R. M., Revil, A., Martin, J. B., Karaoulis, M., and Graham, W. D.: Characterization of groundwater and surface water mixing in a semi confined karst aquifer using time-lapse electrical resistivity tomography, Water Resour. Res., 50, 2566–2585, https://doi.org/10.1002/2013WR013991, 2014. a
Montanari, D., Minissale, A., Doveri, M., Gola, G., Trumpy, E., Santilano, A., and Manzella, A.: Geothermal resources within carbonate reservoirs in western Sicily (Italy): A review, Earth-Sci. Rev., 169, 180–201, https://doi.org/10.1016/j.earscirev.2017.04.016, 2017. a
Niu, Q. and Zhang, C.: Permeability prediction in rocks experiencing mineral precipitation and dissolution: a numerical study, Water Resour. Res., 55, 3107–3121, https://doi.org/10.1029/2018WR024174, 2019. a, b, c, d, e, f, g, h, i
Noiriel, C.: Resolving time-dependent evolution of pore-scale structure, permeability and reactivity using X-ray microtomography, Rev. Mineral. Geochem., 80, 247–285, https://doi.org/10.2138/rmg.2015.80.08, 2015. a
Noiriel, C., Gouze, P., and Bernard, D.: Investigation of porosity and permeability effects from microstructure changes during limestone dissolution, Geophys. Res. Lett., 31, L24603, https://doi.org/10.1029/2004GL021572, 2004. a
Noiriel, C., Bernard, D., Gouze, P., and Thibault, X.: Hydraulic properties and micro-geometry evolution accompanying limestone dissolution by acidic water, Oil Gas Sci. Technol., 60, 177–192, https://doi.org/10.2516/ogst:2005011, 2005. a
Pereira Nunes, J., Blunt, M., and Bijeljic, B.: Pore-scale simulation of carbonate dissolution in micro-CT images, J. Geophys. Res.-Sol. Ea., 121, 558–576, 2016. a
Rembert, F., Jougnot, D., and Guarracino, L.: A fractal model for the electrical conductivity of water-saturated porous media during mineral precipitation-dissolution processes, Adv. Water Res., 145, 103742, https://doi.org/10.1016/j.advwatres.2020.103742, 2020. a, b, c, d, e, f, g
Rembert, F., Jougnot, D., Luquot, L., and Guérin, R.: Interpreting Self-Potential Signal during Reactive Transport: Application to Calcite Dissolution and Precipitation, Water, 14, 1632, https://doi.org/10.3390/w14101632, 2022a. a
Rembert, F., Léger, M., Jougnot, D., and Luquot, L.: Geoelectrical and hydro-chemical monitoring of karst formation at the laboratory scale, Zenodo [data set], https://doi.org/10.5281/zenodo.6908522, 2022b. a
Revil, A., Karaoulis, M., Johnson, T., and Kemna, A.: Some low-frequency electrical methods for subsurface characterization and monitoring in hydrogeology, Hydrogeol. J., 20, 617–658, https://doi.org/10.1007/s10040-011-0819-x, 2012. a
Revil, A., Kessouri, P., and Torres-Verdin, C.: Electrical conductivity, induced polarization, and permeability of the Fontainebleau sandstone, Geophysics, 79, D301–D318, https://doi.org/10.1190/geo2014-0036.1, 2014. a
Rötting, T., Luquot, L., Carrera, J., and Casalinuovo, D. J.: Changes in porosity, permeability, water retention curve and reactive surface area during carbonate rock dissolution, Chem. Geol., 403, 86–98, https://doi.org/10.1016/j.chemgeo.2015.03.008, 2015. a
Saneiyan, S., Ntarlagiannis, D., Ohan, J., Lee, J., Colwell, F., and Burns, S.: Induced polarization as a monitoring tool for in-situ microbial induced carbonate precipitation (MICP) processes, Ecol. Engin., 127, 36–47, https://doi.org/10.1016/j.ecoleng.2018.11.010, 2019. a
Saneiyan, S., Ntarlagiannis, D., and Colwell, F.: Complex conductivity signatures of microbial induced calcite precipitation, field and laboratory scales, Geophys. J. Int., 224, 1811–1824, https://doi.org/10.1093/gji/ggaa510, 2021. a
Soldi, M., Guarracino, L., and Jougnot, D.: A simple hysteretic constitutive model for unsaturated flow, Transport in Porous Media, 120, 271–285, https://doi.org/10.1007/s11242-017-0920-2, 2017. a, b, c
Sun, Z., Mehmani, A., and Torres‐Verdín, C.: Pore‐scale investigation of the electrical resistivity of saturated porous media: flow patterns and porosity efficiency, J. Geophys. Res.-Sol. Ea., 126, e2021JB022608, https://doi.org/10.1029/2021JB022608, 2021. a
Thanh, L., Jougnot, D., Phan, V., and Nguyen, V.: A physically based model for the electrical conductivity of water-saturated porous media, Geophys. J. Int., 219, 866–876, https://doi.org/10.1093/gji/ggz328, 2019. a, b
Vialle, S., Contraires, S., Zinzsner, B., Clavaud, J. B., Mahiouz, K., Zuddas, P., and Zamora, M.: Percolation of CO2‐rich fluids in a limestone sample: Evolution of hydraulic, electrical, chemical, and structural properties, J. Geophys. Res.-Sol. Ea., 119, 2828–2847, https://doi.org/10.1002/2013JB010656, 2014. a, b, c
Vinogradov, J., Hill, R., and Jougnot, D.: Influence of pore size distribution on the electrokinetic coupling coefficient in two-phase flow conditions, Water, 13, 2316, https://doi.org/10.3390/w13172316, 2021. a, b, c
Weger, R. J., Eberli, G. P., Baechle, G. T., Massaferro, J. L., and Sun, Y.-F.: Quantification of pore structure and its effect on sonic velocity and permeability in carbonates, AAPG Bulletin, 93, 1297–1317, https://doi.org/10.1306/05270909001, 2009. a
Winsauer, W. O., Shearin, H. M., J., Masson, P. H., and Williams, M.: Resistivity of brine-saturated sands in relation to pore geometry, AAPG Bulletin, 36, 253–277, https://doi.org/10.1306/3D9343F4-16B1-11D7-8645000102C1865D, 1952. a
Wu, Y., Hubbard, S., Williams, K. H., and Ajo-Franklin, J.: On the complex conductivity signatures of calcite precipitation, J. Geophys. Res., 115, G00G04, https://doi.org/10.1029/2009JG001129, 2010. a
Zhang, C., Slater, L., Redden, G., Fujita, Y., Johnson, T., and Fox, D.: Spectral induced polarization signatures of hydroxide adsorption and mineral precipitation in porous media, Environ. Sci. Technol., 46, 4357–4364, https://doi.org/10.1021/es204404e, 2012. a