Changes in hydrodynamic, structural and geochemical properties in carbonate rock samples due to reactive transport

Reactive transport plays an important role in the development of a wide range of both anthropic and natural processes affecting geological media. To predict the consequences of reactive transport processes on structural and hydrodynamic properties of a porous media at large time and spatial scales, numerical modeling is a powerful tool. Nevertheless, such models, to be realistic, need geochemical, structural and hydrodynamic data inputs representative of the studied reservoir or material. Here, we present an experimental study coupling traditional laboratory measurements and percolation experiments in order to obtain the parameters that define rock heterogeneity, which can be altered during the percolation of a reactive fluid. In order to validate the experimental methodology and identify the role of the initial heterogeneities on the localization of the reactive transport processes, we used three different limestones with different petrophysical characteristics. We tracked the changes of geochemical, structural and hydrodynamic parameters in these samples induced by the percolation of an acid fluid by measuring, before and after the percolation experiment, petrophysical and hydrodynamic properties of the rocks.Peer ReviewedPostprint (published version

Flow simulation of artificially induced microfractures using digital rock and lattice boltzmann methods

Microfractures have great significance in the study of reservoir development because they are an effective reserving space and main contributor to permeability in a large amount of reservoirs. Usually, microfractures are divided into natural microfractures and induced microfractures. Artificially induced rough microfractures are our research objects, the existence of which will affect the fluid-flow system (expand the production radius of production wells), and act as a flow path for the leakage of fluids injected to the wells, and even facilitate depletion in tight reservoirs. Therefore, the characteristic of the flow in artificially induced fractures is of great significance. The Lattice Boltzmann Method (LBM) was used to calculate the equivalent permeability of artificially induced three-dimensional (3D) fractures. The 3D box fractal dimensions and porosity of artificially induced fractures in Berea sandstone were calculated based on the fractal theory and image-segmentation method, respectively. The geometrical parameters (surface roughness, minimum fracture aperture, and mean fracture aperture), were also calculated on the base of digital cores of fractures. According to the results, the permeability lies between 0.071–3.759 (dimensionless LB units) in artificially induced fractures. The wide range of permeability indicates that artificially induced fractures have complex structures and connectivity. It was also found that 3D fractal dimensions of artificially induced fractures in Berea sandstone are between 2.247 and 2.367, which shows that the artificially induced fractures have the characteristics of self-similarity. Finally, the following relations were studied: (a) exponentially increasing permeability with increasing 3D box fractal dimension, (b) linearly increasing permeability with increasing square of mean fracture aperture, (c) indistinct relationship between permeability and surface roughness, and (d) linearly increasing 3D box fractal dimension with increasing porosity

Time-lapse cross-hole electrical resistivity tomography (CHERT) for monitoring seawater intrusion dynamics in a Mediterranean aquifer

Surface electrical resistivity tomography (ERT) is a widely used tool to study seawater intrusion (SWI). It is noninvasive and offers a high spatial coverage at a low cost, but its imaging capabilities are strongly affected by decreasing resolution with depth. We conjecture that the use of CHERT (cross-hole ERT) can partly overcome these resolution limitations since the electrodes are placed at depth, which implies that the model resolution does not decrease at the depths of interest. The objective of this study is to test the CHERT for imaging the SWI and monitoring its dynamics at the Argentona site, a well-instrumented field site of a coastal alluvial aquifer located 40¿km NE of Barcelona. To do so, we installed permanent electrodes around boreholes attached to the PVC pipes to perform time-lapse monitoring of the SWI on a transect perpendicular to the coastline. After 2 years of monitoring, we observe variability of SWI at different timescales: (1) natural seasonal variations and aquifer salinization that we attribute to long-term drought and (2) short-term fluctuations due to sea storms or flooding in the nearby stream during heavy rain events. The spatial imaging of bulk electrical conductivity allows us to explain non-monotonic salinity profiles in open boreholes (step-wise profiles really reflect the presence of freshwater at depth). By comparing CHERT results with traditional in situ measurements such as electrical conductivity of water samples and bulk electrical conductivity from induction logs, we conclude that CHERT is a reliable and cost-effective imaging tool for monitoring SWI dynamics.This work was funded by the project CGL2016-77122-C2-1-R/2-R of the Spanish Government. This project also received funding from the European Commission, Horizon 2020 research and innovation programme (Marie Sklodowska-Curie (grant no. 722028)). The author Albert Folch is a Serra Húnter Fellow.Peer ReviewedPostprint (published version

Time-lapse cross-hole electrical resistivity tomography (CHERT) for monitoring seawater intrusion dynamics in a Mediterranean aquifer

Surface electrical resistivity tomography (ERT) is a widely used tool to study seawater intrusion (SWI). It is noninvasive and offers a high spatial coverage at a low cost, but its imaging capabilities are strongly affected by decreasing resolution with depth. We conjecture that the use of CHERT (cross-hole ERT) can partly overcome these resolution limitations since the electrodes are placed at depth, which implies that the model resolution does not decrease at the depths of interest. The objective of this study is to test the CHERT for imaging the SWI and monitoring its dynamics at the Argentona site, a well-instrumented field site of a coastal alluvial aquifer located 40 km NE of Barcelona. To do so, we installed permanent electrodes around boreholes attached to the PVC pipes to perform time-lapse monitoring of the SWI on a transect perpendicular to the coastline. After 2 years of monitoring, we observe variability of SWI at different timescales: (1) natural seasonal variations and aquifer salinization that we attribute to long-term drought and (2) short-term fluctuations due to sea storms or flooding in the nearby stream during heavy rain events. The spatial imaging of bulk electrical conductivity allows us to explain non-monotonic salinity profiles in open boreholes (step-wise profiles really reflect the presence of freshwater at depth). By comparing CHERT results with traditional in situ measurements such as electrical conductivity of water samples and bulk electrical conductivity from induction logs, we conclude that CHERT is a reliable and cost-effective imaging tool for monitoring SWI dynamics

Heletz experimental site overview, characterization and data analysis for CO2 injection and geological storage

International audienceThis paper provides an overview of the site characterization work at the Heletz site, in preparation to scientifically motivated CO2 injection experiments. The outcomes are geological and hydrogeological models with associated medium properties and baseline conditions. The work has consisted on first re-analyzing the existing data base from ∼40 wells from the previous oil exploration studies, based on which a 3-dimensional structural model was constructed along with first estimates of the properties. The CO2 injection site is located on the saline edges of the Heletz depleted oil field. Two new deep (>1600 m) wells were drilled within the injection site and from these wells a detailed characterization program was carried out, including coring, core analyses, fluid sampling, geophysical logging, seismic survey, in situ hydraulic testing and measurement of the baseline pressure and temperature. The results are presented and discussed in terms of characteristics of the reservoir and cap-rock, the mineralogy, water composition and other baseline conditions, porosity, permeability, capillary pressure and relative permeability. Special emphasis is given to petrophysical properties of the reservoir and the seal, such as comparing the estimates determined by different methods, looking at their geostatistical distributions as well as changes in them when exposed to CO2

Localization and dissolution rates of two different carbonate rocks: role of heterogeneity and mineralogy

International audienc

Quantification expérimentale à l'échelle mésoscopique des processus réactionnels dans le cadre de l'injection de CO2 dans des roches carbonatées et silicatées

In order to minimize CO2 atmospheric concentration, a solution consists in sequestrating CO2 in geological reservoirs. To estimate long term risks, it is necessary to quantify the couplings between reaction processes as well as structural and hydrodynamical modifications. We realised two experimental benches enabling injecting CO2-enriched-brine in conditions corresponding to in situ storage (T an experimental protocol using X-Ray microtomography and fluid and rock analyses in order to measure the variations of physical and chemical parameters.The study of carbonated reservoirs near the injection well, allows quantifying different k-phi relationships depending on the dissolution processes and triggered by the local fluid chemical composition and initials conditions. Away from the injection well, we observe carbonate precipitation decreasing the permeability. The study of fracturated caprock samples shows that alternative percolation of CO2-enriched-brine and CO2 gas increases the fracture permeability. The study of silicated rocks indicates carbonate precipitation in zeolite sandstone and sintered dunite grains. Nevertheless, in zeolite sandstone we also observe the precipitation of clay particles located in the uid pathways which decrease strongly the permeability.Afin de minimiser la concentration en CO2 dans l'atmosphère une solution consiste à le séquestrer dans les réservoirs géologiques. Pour évaluer les risques à long terme, il est nécessaire de quantifier les couplages entre les processus réactionnels et les modifications structurales et hydrodynamiques. Dans cette optique, nous avons construit deux dispositifs expérimentaux permettant d'injecter des saumures chargées en CO2 dans les conditions de stockage (T < 200 °C et P < 200 bar) et mis en oeuvre un protocole expérimental utilisant la microtomographie RX et l'analyse des roches et des fuides pour quantifier les variations des paramètres physiques et chimiques. Dans le cas des réservoirs carbonatés, on observe, près du puits d'injection, une forte variabilité des relations k-phi en fonction du régime de dissolution contrôlé par la chimie locale du fluide, ainsi que par les conditions initiales. A plus grande distance, on observe des processus de précipitation diminuant fortement la perméabilité. L'étude des roches de couverture fracturées (argilite) a montré qu'une percolation alternée d'une saumure chargée en CO2 et de CO2 gaz augmentait la perméabilité de fracture. Dans le cas des roches silicatées, la précipitation de carbonates est mise en évidence aussi bien dans les grès à zéolites que dans les frittés de dunite de San Carlos. Cependant, dans les grès à zéolites, la précipitation d'une phase argileuse est observée dans les chemins d'écoulement et engendre une forte diminution de la perméabilité

Calculating structural and geometrical parameters by laboratory measurements and X-ray microtomography: A comparative study applied to a limestone sample before and after a dissolution experiment

The aim of this study is to compare the structural, geometrical and transport parameters of a limestone rock sample determined by X-ray microtomography (XMT) images and laboratory experiments. Total and effective porosity, pore-size distribution, tortuosity, and effective diffusion coefficient have been estimated. Sensitivity analyses of the segmentation parameters have been performed. The limestone rock sample studied here has been characterized using both approaches before and after a reactive percolation experiment. Strong dissolution process occurred during the percolation, promoting a wormhole formation. This strong heterogeneity formed after the percolation step allows us to apply our methodology to two different samples and enhance the use of experimental techniques or XMT images depending on the rock heterogeneity. We established that for most of the parameters calculated here, the values obtained by computing XMT images are in agreement with the classical laboratory measurements. We demonstrated that the computational porosity is more informative than the laboratory measurement. We observed that pore-size distributions obtained by XMT images and laboratory experiments are slightly different but complementary. Regarding the effective diffusion coefficient, we concluded that both approaches are valuable and give similar results. Nevertheless, we concluded that computing XMT images to determine transport, geometrical, and petrophysical parameters provide similar results to those measured at the laboratory but with much shorter durations