Next Article in Journal
Controls on Mg/Ca Ratios in Recent Stromatolites: Insights from Fluvial Systems in the Iberian Range (Spain)
Next Article in Special Issue
CO2 Dipole Moment: A Simple Model and Its Implications for CO2-Rock Interactions
Previous Article in Journal
Environmental Impact Assessment of the Subsurface in a Former W-Sn Mine: Integration of Geophysical Methodologies
Previous Article in Special Issue
Reaction Characteristics of Two Types of Shale with Supercritical CO2 and Its Potential Impact on Flow-Back Strategies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

CO2 Reaction-Diffusion Experiments in Shales and Carbonates

by
Giordano Montegrossi
1,2,*,
Barbara Cantucci
3,*,
Monica Piochi
4,
Lorenzo Fusi
5,
M. Shahir Misnan
6,
M. Rashad Amir Rashidi
6,
Zainol Affendi Abu Bakar
6,
Zuhar Zahir Tuan Harith
7,
Nabila Hannah Samsol Bahri
7 and
Noorbaizura Hashim
7
1
Centro Nazionale delle Ricerche-Istituto di Geoscienze e Georisorse, Via G. La Pira 4, 50121 Firenze, Italy
2
Consorzio Interuniversitario Nazionale per la Scienza e Tecnologia dei Materiali–INSTM, Via G. Giusti 9, 50121 Firenze, Italy
3
Istituto Nazionale di Geofisica e Vulcanologia, 00143 Rome, Italy
4
Istituto Nazionale di Geofisica e Vulcanologia, Osservatorio Vesuviano, 80125 Naples, Italy
5
Dipartimento di Matematica e Informatica “U. Dini”, Viale Morgagni 67/a, 50134 Firenze, Italy
6
PETRONAS Research Sdn. Bhd, Bandar Baru Bangi 43650, Malaysia
7
AEM Energy Solutions Sdn. Bhd, Kuala Lumpur 50450, Malaysia
*
Authors to whom correspondence should be addressed.
Minerals 2023, 13(1), 56; https://doi.org/10.3390/min13010056
Submission received: 6 November 2022 / Revised: 20 December 2022 / Accepted: 24 December 2022 / Published: 29 December 2022

Abstract

:
The evaluation of caprock integrity and reservoir efficiency is critical for safe CO2 geological storage management. It is therefore important to investigate geochemical reactions between CO2-rich fluids and host rocks and their contribution in retaining CO2 at depth. This study deals with diffusive reaction experiments on shales and carbonate samples cored from an offshore structure in the Malaysian basin, a potential target for CO2-enhanced gas recovery. The aim is to evaluate the CO2 reaction front velocity in a typical shaly caprock and the mineral response of the reservoir. Rock samples were characterized in terms of texture, chemistry, and mineralogy by X-ray diffraction, electron microscopy (SEM), microanalysis (EDS), infrared spectroscopy (FT-IR), rock geochemistry (XRF), and mercury injection capillary permeability (MICP). Performed analyses show mineralogical alteration induced by CO2 as it penetrated into the samples. Carbonate dissolution and weathering of pyrite to form secondary carbonates belonging to siderite-ankerite series were observed along two reaction fronts. Estimated diffusion coefficients of CO2 are two orders of magnitude lower than CO2(aq) molecular diffusion in pure water and from half to an order of magnitude lower than diffusivity computed on unaltered sample, highlighting the important effect of gas–water–rock reactions on the CO2(aq) diffusivities in shales and carbonates. Results obtained in this study provide an insight regarding the effect of geochemical reactions on CO2 transport and represent a further discussion point on the diffusion coefficients.

1. Introduction

CO2 geological storage (CCS) has been recognized as able to efficiently mitigate the climate change induced by greenhouse gasses [1,2,3,4]. In order to ensure a safe storage of CO2 at depth, potential leakage from the reservoir throughout the overburden have to be carefully evaluated [5,6,7,8,9].
In sedimentary basins, shales are generally considered optimal caprock for CCS reservoirs since they are characterized by low permeability and high capillary entry pressure which prevent migration of injected CO2 toward the surface. In addition, they are mainly constituted by siliciclastic and clay minerals, which can react with CO2-rich brines, permanently trapping CO2 as newly formed solid phases [10]. These solid phases are generally present also in sandstones and limestones reservoirs, although in minor amounts. The dissolution of primary minerals and the precipitation of secondary ones control the porosity and the permeability of the hosting rock, ruling the velocity of the CO2 reaction front at the caprock–reservoir interface. Depending on mineral assemblage, fluid chemical composition, and boundary conditions, geochemical reactions may contribute to retain CO2 at depth by self-sealing process or to enhance fluid migration towards the surface by mineral dissolution; and they should be evaluated throughout site-specific analyses.
The characterization of rocks mineralogical and petrophysical alteration due to CO2 injection can be carried out by laboratory analysis and experiments at in-situ condition by using available, site-specific rock samples [11,12,13]. Due to the slow kinetics of some minerals (e.g., siliciclastic and clayey), laboratory experiments allow to evaluate the chemical weathering effect on reservoir and caprock stability only in the short-term, but they can be a useful tool to obtain reliable geochemical input parameters for long-term numerical models, covering uncertainties linked to kinetic rates, mineral reactive surface areas, petrophysical properties of rocks, and diffusion coefficients of supercritical CO2.
In this framework, the objective of this study is to qualitatively evaluate the geochemical reactions induced by CO2 in a carbonate reservoir and on the shaly caprock and quantify the velocity of reaction front by means of reaction-diffusion experiments. A depleted hydrocarbon reservoir, in offshore Malaysia, a potential target for CO2 geological storage [14], was used as a case study for our experiments, focusing on shale overburden, which is often neglected during hydrocarbon exploitation.
The texture, geochemistry, and mineral composition of three selected rock samples were determined before and after alteration experiments, combining several petrophysical techniques: X-ray diffraction, rock geochemistry (XRF), electron microscopy plus energy dispersive system (SEM-EDS), infrared spectroscopy (FT-IR), and mercury injection capillary permeability (MICP).
The lab analysis and diffusion experiments allowed to: (i) characterize the mineralogy of cores; (ii) describe their changes related to geochemical weathering induced by CO2; (iii) evaluate the velocity of CO2 penetration into the samples; and (iv) estimate the CO2 diffusion coefficients. Following our previous methodological approach (e.g., [15]), these results provide information to set and calibrate reliable reactive transport models useful to investigate short-to-long-term consequences of CO2 storage and to evaluate the effect of trapping mechanisms on fluid flow [16].

2. Materials and Methods

2.1. Geological Structure, Rock Formations and Studied Samples

The selected study structure (Figure 1a) is located in the north of Central Luconia Province in the Sarawak Basin, East Malaysia. Sedimentary stratigraphy of the area can be divided in 8 regressive cycles (from the Eocene to present) that are separated by major transgressions [17,18]. In the selected structure, rock formations belong to V-VII sedimentary cycles [19], which are comprised of middle Miocene carbonate reservoir (V cycle) discordantly overlain by Pliocene sediments (VI and VII cycles). The overburden formations are thick, approximately 2000 m, and consist mainly of interbedded silt, silty clay, and shale intervals. They are divided in three main layers based on the contrast in sonic-log velocity: soft sediments right below the seabed, the Seal 6 and below the Seal 4 (Figure 1a). Based on the quick correlation, the presence of silty layers is prominent in Seal 6 while shales dominate Seal 4 with some intercalation of silt. The underlying reservoir is constituted by fractured limestones. Geothermal and geobaric gradients (Figure 1b), measured in well logs, show a high-temperature, high-pressure zone from around 1400 m of depth, roughly corresponding to Seal4 top to the reservoir (133 °C and 315 bar at about 1900 m, top reservoir).
Rock samples employed in diffusion experiments were taken from a deep well, located in the case study structure and provided by operator. Generally speaking, in hydrocarbon exploration very limited attention is reserved for coverture characterization, mainly focusing on the reservoir. Only three rocks samples (C1, C2, and C3), representing Seal 6, Seal 4, and limestone reservoir matrix (Figure 1a), were available for the experiments.

2.2. Sample Preparation and Characterization

Petrophysical analysis of cores was performed before and after the alteration experiments to study textural and mineralogical changes.
In particular, X-ray diffraction (XRPD), X-ray fluorescence (XRF), Rietveld, and electron with energy dispersive system microscopy (SEM-EDS) (Carl Zeiss, Oberkochen, Germany) analyses were carried out before the diffusion experiments to define chemical composition of the samples and semi-quantitative and morphological evaluation of minerals. SEM-EDS analyses were repeated after experiments to evaluate the alteration degree and the velocity of the reaction front. Infrared spectroscopy (FT-IR) (is50, Thermo Fisher Scientific Inc., Waltham, MA, USA) was used to further evaluate the compositional changes, allowing an easy comparison among spectra of pristine and altered materials. Finally, mercury injection capillary pressure (MICP) analyses were also carried out to measure porosity, pore distribution, and tortuosity.
All original samples were massive and well-preserved. After experiments, shale specimens were still cohesive, but after the SEM analyses, preparation, and storage, they progressively lost cohesion (Figure 1d,e), requiring different analytical procedures as described below.
Before starting diffusion experiments, each rock was cut to obtain a representative sample and heated in an oven for five days at 40 °C to avoid a crumbly braking of clays and remove the water present in the samples. This procedure is necessary to evaluate the weight of the sample at dry and saturated conditions before and after the experiments.
A rock portion of initial cores was powdered for XRPD and XRF analysis, whereas another one was cut to create a thin section for SEM-EDS analysis. Finally, a core fraction was devoted to destructive analysis of MICP.
The XRPD analyses were performed on powdered samples by using an X-ray Philips PW 1050/37 diffractometer (PW 1050/37, Philips, Amsterdam, Netherlands), with the X’Pert PRO Philips acquisition system, operating at 40 kV-20 mA, with a Cu anode and graphite monochromator at an interval 2θ of 5–70° and 2–32° for the clay analysis. The powders were mixed with corundum at 2.5% w/w.
XRF was performed on powdered samples with a Rigaku II wavelength dispersive spectrometer (ZSX, Rikaku, Tokyo, Japan) with a Rh anode. The main composition was acquired with a ZSX Rigaku Software Package version 3.50 and quantified with the SQX semi-quantitative elemental analysis routine by the fundamental parameter (FP) method without the use of reference standard samples. In the quantifying process, a sensitivity library (Rigaku analytical system) containing the FP reference values for every element was used. Pure metals and reagents were used for the XRF calibration. The analytic error is 3% w/w for major and 10% w/w for trace elements.
CO2 was determined as loss from carbonates by calcimetric analysis which was used to correct XRF data. The concentration of FeO was calculated by volumetric titration with Potassium dichromate (N: 0.400) (K2Cr2O7) on a quantity of powdered sample (up to 0.5 g), transferred to a 10 mL Pt capsule and solubilized by 10 mL of H2SO4 and 5 mL of HF [20] under a Bunsen burner. Natrium-diphenylamine-sulfonate (solubilized in H2SO4) was used as an indicator. To compute the amount of FeO (in % by weight), the mL of titration solution was multiplied for 0.4 (calibration coefficient for titrating solution).
Rietveld analyses were performed on XRPD data by Maud v2.2 code [21] to obtain the modal mineralogy. Maud v2.2 is a general diffraction/reflectivity analysis program mainly based on the Rietveld method which permits to calculate the reciprocal ratios among present phases and normalize them. Carbonate content was determined by calcimetry with a Dietrich–Fruhling apparatus, including correction for carbonates other than calcite, based on the XRPD Rietveld refinement data. XRF results for total chemical composition were then corrected including carbonate abundances.
Scanning electron microscopy equipped with an energy dispersive system (SEM/EDS) on thin sections (preliminary carbon coated) allowed rock texture and semi-quantitative chemical analysis. Morphological and morphometrical images at high resolution in a back scattered mode and micro-chemical analyses were obtained using a ZEISS EVO MA15 scanning electron microscope with microanalysis (EDS) OXFORD INCA 250 at the M.E.M.A lab. Analytical conditions include a working distance of 9mm, voltage 15–20 kV and spot size for EDS of 2 µm. Error is estimated at 5% for major elements.
After diffusion experiments SEM-EDS analyses were also performed on material positioned on stub to highlight mineralogical changes due to CO2–water interaction. In this case, the most cohesive portion of C1 and C2 specimens was taken carefully considering the orientation with respect to the fluid interaction front during the experiment.
FT-IR analyses were carried out by a Nicolet is50 instrument equipped with a Polaris™ long-life IR source (is50, Thermo Fisher Scientific Inc., Waltham, MA, USA). Spectra were acquired in the 7400–400 cm−1 range on grain fraction taken from the original samples and from the altered zone after the experiments. Furthermore, the original C3 and C1 samples and C3 altered sample allowed analysis as blocks in NIR mode. For the NIR range, sample portions and blocks were positioned on a glass lamina and analyzed in the Integrating Sphere module using a InGaAs detector and a CaF2 beamsplitter; settings included 8 resolution, 1 sample gain, 70 aperture. The rest of spectrum range was acquired on grain portions pressed under a built-in mid- and far-IR capable diamond ATR using the DTGS detector and the KBr beamsplitter; settings included 4 resolution, 1 sample gain, 100 aperture. Hence, 50 scans were acquired in both cases. An air dryer fluxing air up to 34 l/m at the dew point of –75 °C and low (down to 1 ppm) CO2 was used to purge the system. The OMNIC Data Collector 9.11© provided instrument control, data acquisition, and spectra elaboration.
Finally, porosity, pore distribution, and tortuosity were determined by MICP on 1 cm3 specimens from original samples, after drying to 105 °C. Thermo Scientific™ Pascal Mercury Porosimeter (PASCAL 120 and 240) instruments (PASCAL 120 and 240, Thermo Fisher Scientific Inc., Waltham, MA, USA) were used to determine macropores (with diameter ranging from 116 to 3.8 µm) and micropores (from 15 to 0.0074 µm) working up the maximum pressure of 200 MPa and temperature of 22 °C. The samples were first submitted to a vacuum to degas the material and to permit the mercury filling of the dilatometer. The PASCAL instruments were managed by the S.O.L.I.D. (SOLver of Intrusion Data) software (by Thermo Fisher Scientific) to elaborate raw data.

2.3. Experimental Settings

The diffusion experiments allow to define the velocity of reactive reaction front in a rock by placing a solution with high concentration of CO2 in contact with an intact rock for a certain period of time; in this study, for 2–5 days.
Experiment duration was calibrated based for both the small dimension of the specimens and estimated velocity of the diffusion front, to allow a noticeable CO2 diffusion (some mm) while leaving an unaltered core zone, to be able to detect the altered/unaltered boundary. In this short period of time, only fast kinetic reactions (i.e., carbonate, oxides and pyrite) and clay mineral structure changes were expected to occur after weathering, since clay minerals take within weeks to months to fully react [22].
Original carbonate reservoir sample (C3) was cut in two cubes of 1 × 1 × 1 cm each and positioned on the base of the micro-reactor (Figure 2). In this way, the reaction with CO2-rich brine proceeds symmetrically from the external to internal of the cube by means of diffusion only. Shale samples of Seal 6 and Seal 4 (C1 and C2, respectively), due to the clay swelling and sample dissolution in water, were cut in cylinders of about 1 × 3 cm and inserted in a custom-made reaction vessel in high density polytetrafluoroethylene (PTFE) open from two sides to allow the reaction.
The experimental apparatus is a micro-reactor 4590 HP/HT (manufactured by Parr Instrument Company, Moline, IL, USA), equipped with perfluoroelastomer (FFKM) O-ring for temperatures up to 300 °C, 200 bar. In each experiment, the micro-reactor was charged with one specimen and 5.5 mL of MilliQ water (enough to completely submerge the sample). The micro-reactor was initially flushed with pure CO2(g) to get rid of atmospheric air before setting the CO2 partial pressure (PCO2). The specimen within the micro-reactor was completely immersed in water and covered by a thin layer of water, separating the gas phase, above, from the rock below (Figure 2). This thin water layer is crossed by a diffusive CO2 flux controlling reaction front.
Experiments were carried out at temperatures of 75 °C, 105 °C and 145 °C, corresponding to Seal 6 (C1), Seal 4 (C2), and reservoir (C3) core samples at in-situ conditions (Figure 1b) and PCO2 of 60 bars. Experimental pressure, employed for all the experiments, was subject to CO2 cylinder availability, a maximum pressure of 60 bars.
After experiments, the samples were taken out from the apparatus, dried in a 40 °C oven to recover mechanical consistency of the material, extracted from the reaction vessel, and weighted again. The specimens were then cut in half, from the CO2-rich brine exposed side toward the inner part.

3. Results

3.1. Cores Characterization

Original petrographic features of core samples were defined at various scales by XRPD, XRF, FT-IR, and SEM-EDS analyses.
XRPD combined with Rietveld refinement (Figure 3) allowed a semi-quantitative mineralogical composition normalized to 100 (Table 1). Analytical results show that the shale and shaly silt overburden (C1 and C2 samples) is constituted by ca. 36 and 51% by weight of quartz and plagioclase (anorthite by Rietveld analysis), respectively. Clay minerals (illite, kaolinite and chlorite) contribute from 31 to 40% w/w with minor mica (muscovite) content of 11% w/w. Finally, low amounts of carbonates (calcite, dolomite, and siderite) are present (C1 sample = 7% w/w, C2 sample = 14% w/w). If we compare Seal 6 and Seal 4 samples, this latter shows the major contents of carbonate (14%) and clay minerals (40% w/w) and the minor contents of quartz and anorthite (sand = 36% w/w). The reservoir (C3 sample) is constituted by limestone (76% w/w), mainly calcite and dolomite, subordinate quartz (13% w/w), anorthite and mica, and a very low abundance of clay minerals (10% w/w; illite, chlorite and kaolinite).
FT-IR spectra (Figure 4 and Figure 5) corroborate the XRPD results indicating carbonate and clay prevalence in the C3 and in the C1 samples, respectively. Actually, the C3 spectrum is dominated by the signals of calcite with typical spikes of absorption bands at 1396 (𝜈3 CO3−2asymmetrical stretch), 871 (𝜈2 CO3−2asymmetrical bend) and 711 (𝜈4 CO3−2aymmetrical bend) cm−1 [23]. A small contribution from the substitution of other (Mg, Fe, Mn) cations can be suspected by comparing the smoothed peak at 5154 cm−1 of our sample with those of siderite (HS271.3B) in the USGS database (https://www.usgs.gov/labs/spectroscopy-lab (accessed on 1 November 2022)). Considering the USGS database, the signals at 7092 and 7246 cm−1 in the carbonate sample highlight the scarce presence of mica and clays, although the doublet at 7246 and 7.092 cm−1 can be tentatively attributed to the kaolinite occurrence. Conversely, for the investigated C1 shale, the illite determines the peaks sequence in the NIR frequency (7092, 5236, 4504, 4255 cm−1) and the shape in the OH-stretching region. In this later range, the mica contribution can be also revealed at ca. 3645 cm−1, as well as the triplet at ca. 2656, 2923, and 2853 cm−1. According to [24,25], several peaks at ca. 1030, 1000, and 984 cm−1 (Si-O stretching), at ca. 933, 906, and 870 cm−1 (vibration of AlAlOH and AlFeOH groups) and at ca. 543, 520, and 506 cm−1 (Al-O and Si-O deformation) can be attributed to clays (mostly illite/montmorillonite and kaolinite). Notably, for C1, the triplet at 2957 + 2922 + 2850 cm−1 and the smoothed signals at 1421 cm−1 point to muscovite being richest and a significant occurrence of (ankerite, dolomite, siderite) carbonates, respectively. Based on the peak at 3619 and the hamp at 3415 cm−1, it is demonstrated that montmorillonite/illite is the major clay in shales.
XRF results are reported in Table 2, also considering CO2 loss correction. The overburden shales mostly contain SiO2 (38–41% w/w) and Al2O3 (12–16% w/w), with minor amounts of K2O (2.6% w/w) and Fe2O3 (5.2% w/w). Boron and magnesium oxides are also present. Moreover, for XRPD, these oxides suggest the presence of quartz and clay minerals rich in potassium and iron, such as illite, chlorite, and muscovite. In the C3 sample, the most abundant oxide is CaO (40.2% w/w) with minor amounts of SiO2 (13% w/w), Al2O3 (6% w/w), Fe2O3 (1.6% w/w), and SO3 (1.33% w/w). The presence of SO3 is likely due to pyrite contribution, which was recognized in XRPD but not quantified due to the very low amount.
SEM-EDS investigation of Seal 6 shale (C1 sample) at a low magnification scale (Figure 6a) shows the presence of some phenocrysts of muscovite (mean 24 × 33 µm), Fe-rich dolomite-ankerite (mean 24 × 37 µm), and anatase (mean 6 × 12 µm), which are dispersed in a dominant finer matrix. The matrix is cohesive and generally low porous. The total porosity, measured by MICP, is 8.7%, with an asymmetrical distribution of pores and mean frequency of 0.019 μm. Tortuosity is 2.174. The zoom on the matrix (Figure 6b,c) reveals that it is constituted by microcrysts of quartz (mean 1.53 × 2.6 µm), secondary illite (mean 0.7 µm), chlorite (mean 0.93 × 1.33 µm), and kaolinite (mean 0.62 × 3.19 µm). The clay minerals exhibit both spherical and plate like shape. Illite and chlorite are generally rich in iron ion, in variable content ranging from 0.51% to 7.03%, in agreement with marine or hypersaline environments where the sedimentation rate is low. Indeed, the sources of K (for illite) and Fe ions are seawater and detrital minerals, respectively [26]. Cubic siderite-dolomite microcrystals (mean 1.11 × 1.99 µm) are also present. Primary plagioclase (5.5 µm, spectrum 2) and quartz are visible in Figure 6c.
Figure 6. SEM-EDS data of C1 sample in thin section. (a) General aspect at low magnification scale at the analytical Site 1. S1: muscovite with quartz; S2 and S3: solid solution of Fe-rich dolomite and ankerite; S4: anatase. (b) Texture and mineralogy of the analytical Site S2. S1: Fe-rich illite; S2 and S3: quartz; S4: illite. (c) Texture and mineralogy of the analytical Site 3. S1: solid solution of siderite-ankerite; S2: kaolinite with adsorbed Na from underlying primary plagioclase; S3: solid solution of Fe-rich chlorite and illite; S4: kaolinite. (d) EDS spectrum of S1 in panel (a) compatible with mixture of muscovite and quartz. Mineral compositions are related to spectrum for site in Table 3.
Figure 6. SEM-EDS data of C1 sample in thin section. (a) General aspect at low magnification scale at the analytical Site 1. S1: muscovite with quartz; S2 and S3: solid solution of Fe-rich dolomite and ankerite; S4: anatase. (b) Texture and mineralogy of the analytical Site S2. S1: Fe-rich illite; S2 and S3: quartz; S4: illite. (c) Texture and mineralogy of the analytical Site 3. S1: solid solution of siderite-ankerite; S2: kaolinite with adsorbed Na from underlying primary plagioclase; S3: solid solution of Fe-rich chlorite and illite; S4: kaolinite. (d) EDS spectrum of S1 in panel (a) compatible with mixture of muscovite and quartz. Mineral compositions are related to spectrum for site in Table 3.
Minerals 13 00056 g006
Table 3. Semi-quantitative analysis of C1 sample in atoms percentage.
Table 3. Semi-quantitative analysis of C1 sample in atoms percentage.
SpectrumMineralONaMgAlSiKCaTiMnFe
1-Site1muscovite + quartz64.700.64-5.5924.154.92----
2-Site1Fe-dolomite + ankerite74.57-6.64-0.52-13.38-1.213.69
3-Site1Fe-dolomite + ankerite73.96-6.16-0.60-15.03--4.25
4-Site1anatase68.700.67 *-0.25 *0.55 *--29.83--
1-Site2Fe-rich illite 65.391.122.076.3819.250.870.49--4.42
2-Site2quartz65.45--0.31 *34.24-----
3-Site2quartz66.55--0.38 *33.07-----
4-Site2Illite64.911.230.566.0325.930.83---0.51
1-Site3siderite-dolomite68.06-7.621.82 *2.92 *-1.79--17.78
2-Site3kaolinite on primary plagioclase68.354.65-7.7917.560.171.49---
3-Site3Fe-chlorite+illite68.842.853.444.8111.50-1.55--7.03
4-Site3kaolinite66.302.17-5.6024.52-1.23--0.19
*: Contamination from close by minerals. Refer to Figure 6 for details on spectrum and site number.
SEM analysis of Seal 4 (C2 sample) at low magnification (Figure 7a) shows a slightly coarser and more porous structure than C1 sample. Based on MICP, this sample has a total porosity of 9.2%, with an average pore distribution of 0.15 μm. Tortuosity is 2.164. At the above spatial scale, only pyrite (mean 10 µm) and quartz macro crystals (mean 4.8 × 12.7 µm) are clearly distinguishable. An enlarged magnification (Figure 7b) allows further distinguishing cubic crystals of dolomite-siderite (mean 2.8 µm), small crystals of anatase (mean 1.6 µm) (Table 4) and, in the matrix, illite, both spherical (mean 2.15 µm) and plate shape (mean 1.18 × 5.05 µm), kaolinite (mean 0.71 × 4.83 µm), and quartz (mean 12 × 19 µm). The zoom on the matrix (Figure 7c, Table 4) shows that it is mostly constituted by secondary illite-chlorite rich in Fe, both spherical (mean 1.62 µm) and plate shape (mean 1.03 × 2.38 µm), kaolinite (mean 0.28 × 1.35 µm), muscovite (mean 0.84 × 2.3 µm), and micro crystals of dolomite-siderite (mean 2.5 µm).
SEM analysis of carbonate reservoir sample (C3) at a low magnification scale (Figure 8a) shows the presence of primary calcite (Table 5) and clay with a plate-like structure (mean 0.5 × 2.23 µm). The porosity, as measured by MICP, is 9.4%, close to the C2 sample but C3 has a more symmetrical distribution, with average porous dimension ranging from 0.0723 to 0.0792 μm. Tortuosity is 2.125. A fossil re-crystallized by small crystals of pyrite (mean 1 µm) is present in the center of the picture. The zoom on the matrix (Figure 8b,c) reveals that it is constituted by secondary calcite (mean 2.19 µm), kaolinite (mean 0.8 × 6.4 µm), illite (mean 1.2 × 4.75 µm), muscovite (mean 0.52 × 5.7 µm), and quartz (mean 0.92 × 4.2 µm). Rare crystals of pyrite (3–5 µm) and sylvite are also present.

3.2. Diffusion Experiments

The reaction front due to fluid–rock interaction was evaluated by SEM-EDS analysis after two and five days of interaction inside the micro-reactor. As expected, due to the short interaction time, the results indicate limited sample changes that mostly concerned with carbonates and pyrite.
SEM-EDS investigations were conducted in two modes: (i) classical electron microscope image (e.g., Figure 9), and (ii) chemical map (cameo image) by chemical fake-coloring (Figure 10). The color scale was arbitrarily changed to highlight the reaction front, thus not reflecting the real abundance of the elements. The first analysis aimed to recover chemical changes and to spot secondary minerals, whereas the cameo images were used to define the boundary between unaltered sample and the weathering due to CO2 reaction.
A general feature observed on the SEM images after the experiments is the depletion in the pyrite microcrystals and a relative increased abundance of carbonates mostly in the form of micron sized cubes (Figure 9). These carbonates, still pertaining to the dolomite-siderite-ankerite series, generally resulted depleted in Ca in favor of Mg and Fe at the EDS, exhibiting an evident rhombohedral habit (C3 sample, Figure 11) and believed to represent secondary phases. Figure 10a (C1 sample) shows this phenomenon, with a relative enrichment of Fe and Mg (red and blue colors, respectively) with respect to Ca (green color) in the portion of the specimen which reacted with CO2. In Figure 10b (C2 sample) a relative enrichment of Mg (green) is visible within the reaction area. The Mg relative enrichment and false scale color jeopardize the original major abondance of Fe and K (visible on the other side of reaction front) of C2 with respect to C1 samples. The reaction front of C3 sample (Figure 10c) is slightly evident since Ca and Mg are the major elements constituting carbonates. The green color represents the iron content, which is more evident in the unaltered sample portion, thus allowing the front identification by the color intensity. The different coloring across the reaction front is due to the weathering process, which occurs as CO2(aq) reached the minerals, and its effect on the coloring is not easily predictable, since it could be due to both matrix depletion after weathering and/or subsequent deposition of leached elements.
The distribution of minerals and elements at mm-scale allows to define two reaction fronts, the first corresponding to the area of pyrite depletion and the second coinciding with the precipitation of secondary carbonates.
The distance of these fronts from the outer part of the specimen allows defining the CO2 penetration depth (e.g., Figure 9) after two and five days. These distances were measured by SEM ruler, whose instrumental accuracy is within 1 µm, but considering the difficulties in localizing the alteration front and its unevenness, we can reasonably estimate the experimental error to be approximately 10 µm. The measured penetration front and the depletion front length are reported in Table 6.
FT-IR spectra in Figure 4 and Figure 5 indicate no significant variation in shales subjected to diffusion experiments, whereas after two days of hydrothermal action, the C3 carbonate significantly modified its infrared response. In particular, several peaks appeared in the OH- stretching region and between 1500 and 500 cm−1 of the experimental C3 sample, i.e., at ca. 3695, 3648, 3633, 3618, 932–910, 540–470 cm−1(Figure 4), pointing to kaolinite plus illite/muscovite [24,25,27]. The doublet at 3695 and 3618 cm−1 is diagnostic for the kaolin phase and is corroborated by absorption at 933–929 and 906 cm−1 (OH-bending/vibrations). The illite/muscovite is indicated by signals at ca. 1030, 912, and 520 cm−1. The spectra also show a span of both 𝜈3 CO32− asymmetrical stretch and bend bands, the first band assuming a bell-shaped between 1384 and 1434 cm−1 and the second slightly distorting towards higher and lower wavenumbers (Figure 12a,b). Following [23], this indicates Mg, Fe carbonate and, in our case, may suggest calcite dissolution in agreement with the experimental sample aspect under the optical microscope (Figure 12c). Notably, shales are slightly different, maintaining the more accentuated peaks at 3651, 3619, and 932 cm−1 for the kaolinite abundance in sample C2, and for C1, the triplet at 2957 + 2922 + 2850 cm−1 for muscovite being richest, and the smoothed signals at 1421 cm−1 for the occurrence of ankerite, dolomite, and siderite carbonates, in addition to the illite (3622 and 3415 cm−1). However, in these experimental samples, the infrared signals at 800, 777, and 747 cm−1 (Figure 4) show faint but sensible variations that would be interest to OH bending vibrations/interactions with cations, thus pointing to possible leeching processes. The 1000 and 900 cm−1 bands, related to Si-O-Si stretching and Al-Al-OH bending [28], can be tentatively attributed to structural modification of the illite clays or muscovite with possibly consequent cation losses.

4. Discussion

The performed study aimed to explore the penetration depth and time of CO2 during hydrothermal experiments. Despite the limitations presented by the analytical strategy and experimental samples, the obtained results show fast reactive minerals alteration.
A few (two and five) days allows the weathering of only fast reactive minerals, i.e., carbonate, sulfides, and sulfates [29], and a possible leaching from clays. In fact, we observed major alteration in the C3 carbonate sample and minor changes in the C1 and C2 shales, in all cases mainly concerning carbonates and pyrite. Based on the both SEM-EDS and FT-IR results, ankerite is more represented in altered samples as results of dissolution of calcite or of carbonate dissolution and recrystallization. In this latter case, carbonate traps Fe+Mg>Ca with the iron released by disappeared pyrite, although we can guess that some iron is originated from the leaching of Fe-rich minerals. An oxidation process during the experiments, which could be responsible of pyrite depletion, is unlikely since all oxygen was removed by using disareated water and flushing the cell with pure CO2. Moreover, no sulfate phase was observed. The carbonate dissolution determines a relative sample enrichment in clay (kaolinite over montmorillonite and illite). In the shales, the FT-IR spectra suggest that experiments produced modification of the illite/muscovite structure, likely releasing cations.
The mineralogical changes in the matrix after two and five days are attributed to the action of CO2 penetrated in the samples. Considering geochemical reactions, the distance of the penetration front and depletion front from the outer edge of the specimens was defined, with an estimated experimental error of approximately 10 µm (Table 6).
The penetration length L (m) over the time t (s) may be used to roughly compute the apparent diffusion coefficient Draw (m2 s−1) according to Equation (1) (Table 6). This rough evaluation should be conducted only for the penetration length, while the second front observed in the diffusion experiment is related to mineral depletion and should be better described by geochemical modeling.
D r a w = L 2 t
Equation (1) is a simple relationship between the only two parameters, L and t, we can measure by SEM images after the experiments. Equation (1) derives from the analytical solution of diffusion equation in one dimension and has been already used to evaluate the apparent diffusion coefficient in cements after experiments with CO2 by [30,31,32].
The proposed approach is different from the computation of the effective diffusivity of CO2 in a porous media, De (m2 s−1), considering the CO2 diffusion coefficient in water (Dw, m2 s−1), porosity (φ), and tortuosity (τ) (Equation (2)).
Indeed, in our case, the diffusivity values computed by Equation (1) are ruled by CO2 induced reactions, whereas those computed by Equation (2) (Table 6) can refer only to unreacted sample, and a more complex model accounting for geochemical reactions should be performed to obtain an accurate computation instead of Equation (2). The difference in diffusive transport properties between unaltered and altered rocks can be significant, as shown by [33] for shale caprocks.
D e = ϕ τ 2 D w  
where Draw and De computed for C1, C2, and C3 samples are reported in Table 6, where Dw values were taken from [34] at 75 °C, 105 °C, and 145 °C, and φ and τ values were measured by MICP tests on original core samples (see Section 3.1).
Estimated apparent diffusion coefficients for both two- and five-day experiments (with five-day experiments that should prove to be more reliable) range from 2.80 × 10−11 m2 s−1 to 5.18 × 10−11 m2 s−1 for shale samples, and from 3.11 × 10−11 to 3.36 × 10−11 m2 s−1 for carbonate. The experimental error ranges from 1.6 × 10−13 m2 s−1 to 3.5 × 10−13 m2 s−1 and its impact on diffusion coefficients may be considered limited.
Draw values of shales and carbonate are of the same magnitude, with those of shales being slightly higher than carbonate for two-day experiments, whereas after five days, Draw is higher in C3 than C1 and C2. [35] showed that changes in transport properties in limestones could be affected also by a redistribution of calcite within the sample, as well as weathering of siliciclastic and clay minerals. Computed Draw values are in agreement with previously published data for shale samples, ranging from 3.08 × 10−11 up to 7.80 × 10−11 m2 s−1 [35,36,37,38,39] (Table 7). Differences in the values, within the same order of magnitude, can be due to several factors. First of all, depending on mineral assemblage, dissolution and precipitation may change pore network characteristics. [33] demonstrated a De increasing in shales from 5 × 10−12 m2 s−1 to 1.9 × 10−11 m2 s−1 due to redox-sensitive and carbonate reactions. Moreover, CO2 sorption processes on shales may retard the diffusion front depending on organic materials, clay minerals, and micropore structure [36]. Furthermore, the effective diffusivity depends on experimental temperature and PCO2, since CO2 diffusion in water increases with them [34].
Computed effective diffusion coefficients are 9.77 × 10−11 m2 s−1 for C1, 1.53 × 10−11 m2 s−1 for C2, and 2.51 × 10−10 m2 s−1 for C3 sample. These values are from half to an order of magnitude higher than Draw diffusivity ones, highlighting the important contribution of chemical reactions in retarding CO2 diffusivity throughout hosting rock since CO2(aq) is consumed by chemical reaction and adsorption.
Not many experimental laboratories deal with the complex topic of measuring the effect of geochemical reactions on the CO2 diffusion coefficient. This is mostly because the results (Table 7) are two orders of magnitude lower than CO2(aq) molecular diffusion in pure water, mainly due to pore network changes and chemical reactions.

5. Conclusions

Diffusive reaction experiments were carried out on two shales and a carbonate rock sample belonging to an offshore structure in the Malaysian basin with the goals to evaluate the geochemical reactions induced by CO2 and quantify the CO2 reaction front velocity in the caprock and reservoir at the respective in-situ temperature conditions. Detailed laboratory analysis, performed before and after diffusion experiments, allowed to characterize the texture, chemistry, and mineral association of rock samples.
Results obtained using SEM-EDS illustrate a general behavior of samples with the alteration of pyrite in favor of carbonate secondary minerals belonging to siderite-ankerite series, along two identified reaction fronts.
FT-IR spectra clearly show the carbonate dissolution/precipitation processes, and point out illite structure changes, likely leading to the loss of cations, in shale samples.
The diffusion coefficients estimated by the penetration length of CO2 are two orders of magnitude lower than CO2(aq) molecular diffusion in pure water.
Draw coefficients have the same magnitude order of the coefficient reported for similar samples in the recent literature, and from half to an order of magnitude lower than De computed on an unaltered sample. These differences are mainly imputed to gas–water–rock reactions and adsorption process in shales, which retard the diffusion front.
Results obtained in this study provide an insight concerning the mineralogy of Seal 6, Seal 4, and reservoir formations and represent a discussion point on the obtained diffusion coefficients. Carbon dioxide that migrates from the reservoir into the cap rock reacts with minerals, hence reducing the slow, diffusion-driven migration.
Obtained information here reported may be also useful to set up future numerical simulations of reactive transport flow in geological structures potentially suitable for CO2 storage.

Author Contributions

Conceptualization, G.M., B.C., M.P., L.F. and Z.Z.T.H.; methodology, G.M., B.C. and M.P.; laboratory experiments and analyses, G.M. and M.P.; investigation, G.M., B.C. and M.P.; data curation, B.C. and M.P.; writing—original draft preparation, B.C. and M.P.; writing—review and editing, B.C., M.P. and G.M.; project administration, N.H.S.B. and N.H.; funding acquisition, M.S.M., M.R.A.R. and Z.A.A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by PETRONAS Research Sdn. Bhd for the Project “Methodology of Mitigation CO2 Leakages (Overburden Study)”.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank E. Pecchioni for her help during XRPD analyses and PETRONAS Research Sdn Bhd for the cooperation and help in this work. The two anonymous reviewers are thanked for their critical and thoughtful evaluation of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Metz, B. Carbon Dioxide Capture and Storage: Special Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2005. [Google Scholar]
  2. Metz, B.; Davidson, O.R.; Bosh, P.; Dave, R.; Meyer, L.A. Contribution of Working Group III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, 2007; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2007. [Google Scholar]
  3. IEA Greehouse Gas R&D Programme. Natural and Industrial Analugues for Geological Storage of Carbon Dioxide; IEA Greehouse Gas R&D Programme: Cheltenham, UK, 2009. [Google Scholar]
  4. Pachauri, R.K.; Meyer, L.A. IPCC, 2014: Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; IPCC: Geneva, Switzerland, 2014. [Google Scholar]
  5. Pruess, K.; García, J. Multiphase Flow Dynamics during CO2 Disposal into Saline Aquifers. Environ. Geol. 2002, 42, 282–295. [Google Scholar] [CrossRef]
  6. Rutqvist, J.; Tsang, C.-F. A Study of Caprock Hydromechanical Changes Associated with CO2-Injection into a Brine Formation. Environ. Geol. 2002, 42, 296–305. [Google Scholar] [CrossRef]
  7. Damen, K.; Faaij, A.; van Bergen, F.; Gale, J.; Lysen, E. Identification of Early Opportunities for CO2 Sequestration—Worldwide Screening for CO2-EOR and CO2-ECBM Projects. Energy 2005, 30, 1931–1952. [Google Scholar] [CrossRef]
  8. Voltattorni, N.; Sciarra, A.; Caramanna, G.; Cinti, D.; Pizzino, L.; Quattrocchi, F. Gas Geochemistry of Natural Analogues for the Studies of Geological CO2 Sequestration. Appl. Geochem. 2009, 24, 1339–1346. [Google Scholar] [CrossRef]
  9. Balashov, V.N.; Guthrie, G.D.; Lopano, C.L.; Hakala, J.A.; Brantley, S.L. Reaction and Diffusion at the Reservoir/Shale Interface during CO2 Storage: Impact of Geochemical Kinetics. Appl. Geochem. 2015, 61, 119–131. [Google Scholar] [CrossRef]
  10. Gunter, W.D.; Perkins, E.H.; McCann, T.J. Aquifer Disposal of CO2-Rich Gases: Reaction Design for Added Capacity. Energy Convers. Manag. 1993, 34, 941–948. [Google Scholar] [CrossRef]
  11. Shafeen, A.; Carter, T. Geological Sequestration of Greenhouse Gases. In Environmentally Conscious Fossil Energy Production; Kutz, M., Elkamel, A., Eds.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2009; pp. 207–241. ISBN 978-0-470-43274-7. [Google Scholar]
  12. Amann-Hildenbrand, A.; Bertier, P.; Busch, A.; Krooss, B.M. Experimental Investigation of the Sealing Capacity of Generic Clay-Rich Caprocks. Int. J. Greenh. Gas Control. 2013, 19, 620–641. [Google Scholar] [CrossRef]
  13. Shogenov, K.; Shogenova, A.; Vizika-Kavvadias, O.; Nauroy, J.-F. Experimental Modeling of CO2-Fluid-Rock Interaction: The Evolution of the Composition and Properties of Host Rocks in the Baltic Region: EXPERIMENTAL MODELING OF CO2 STORAGE. Earth Space Sci. 2015, 2, 262–284. [Google Scholar] [CrossRef]
  14. Masoudi, R.; Jalil, M.A.A.; Press, D.; Lee, K.H.; Tan, C.P.; Anis, L.; Darman, N.; Othman, M. An Integrated Reservoir Simulation-Geomechanical Study on Feasibility of CO2 Storage in M4 Carbonate Reservoir, Malaysia. In Proceedings of the Society of Petroleum Engineers—International Petroleum Technology Conference 2012, IPTC 2012, Bangkok, Thailand, 7–9 February 2012; Volume 3, pp. 2583–2600. [Google Scholar] [CrossRef]
  15. Piochi, M.; Cantucci, B.; Montegrossi, G.; Currenti, G. Hydrothermal Alteration at the San Vito Area of the Campi Flegrei Geothermal System in Italy: Mineral Review and Geochemical Modeling. Minerals 2021, 11, 810. [Google Scholar] [CrossRef]
  16. Cantucci, B.; Montegrossi, G.; Buttinelli, M.; Vaselli, O.; Scrocca, D.; Quattrocchi, F. Geochemical Barriers in CO2 Capture and Storage Feasibility Studies. Transp. Porous Med. 2015, 106, 107–143. [Google Scholar] [CrossRef]
  17. Doust, H. Geology and Exploration History of Offshore Central Sarawak. In Energy Resources of the Pacific Region; American Association of Petroleum Geologists: Tulsa, OK, USA, 1981; ISBN 978-1-62981-180-2. [Google Scholar]
  18. Vahrenkamp, V.C. Miocene Carbonates of the Luconia Province, Offshore Sarawak: Implications for Regional Geology and Reservoir Properties from Strontium-Isotope Stratigraphy. BGSM 1998, 42, 1–13. [Google Scholar] [CrossRef]
  19. Ho, K.F. Stratigraphic Framework for Oil Exploration in Sarawak. Bull. Geol. Soc. Malays. 1978, 10, 1–13. [Google Scholar] [CrossRef]
  20. Shapiro, L.; Brannock, W.W. Rapid Analysis of Silicates, Carbonates and Phosphate Rocks; US Government Printing Office: Washington, DC, USA, 1962.
  21. Lutterotti, L.; Bortolotti, M.; Ischia, G.; Lonardelli, I.; Wenk, H.R. Rietveld Texture Analysis from Diffraction Images. Z. Fur Krist. Suppl. 2007, 1, 125–130. [Google Scholar] [CrossRef]
  22. Wigand, M.; Carey, J.W.; Schütt, H.; Spangenberg, E.; Erzinger, J. Geochemical Effects of CO2 Sequestration in Sandstones under Simulated in Situ Conditions of Deep Saline Aquifers. Appl. Geochem. 2008, 23, 2735–2745. [Google Scholar] [CrossRef]
  23. Kim, Y.; Caumon, M.-C.; Barres, O.; Sall, A.; Cauzid, J. Identification and Composition of Carbonate Minerals of the Calcite Structure by Raman and Infrared Spectroscopies Using Portable Devices. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2021, 261, 119980. [Google Scholar] [CrossRef]
  24. Madejová, J.; Komadel, P. Baseline studies of the clay minerals society source clays: Infrared methods. Clays Clay Miner. 2001, 49, 410–432. [Google Scholar] [CrossRef]
  25. Madejová, J. FTIR Techniques in Clay Mineral Studies. Vib. Spectrosc. 2003, 31, 1–10. [Google Scholar] [CrossRef]
  26. Meunier, A.; El Albani, A. The Glauconite-Fe-Illite-Fe-Smectite Problem: A Critical Review. Terra Nova 2007, 19, 95–104. [Google Scholar] [CrossRef]
  27. Vaculíková, L.; Plevová, E. Identification of Clay Minerals and Micas in Sedimentary Rocks. Acta Geodyn. Geomater. 2005, 2, 167–175. [Google Scholar]
  28. Di Pietro, S.A.; Emerson, H.P.; Katsenovich, Y.P.; Johnson, T.J.; Francis, R.M.; Mason, H.E.; Marple, M.A.; Sawvel, A.M.; Szecsody, J.E. Solid Phase Characterization and Transformation of Illite Mineral with Gas-Phase Ammonia Treatment. J. Hazard. Mater. 2022, 424, 127657. [Google Scholar] [CrossRef]
  29. Palandri, J.L.; Kharaka, Y.K. A Compilation of Rate Parameters of Water-Mineral Interaction Kinetics for Application to Geochemical Modeling; Geological Survey: Menlo Park, CA, USA, 2004.
  30. Kutchko, B.G.; Strazisar, B.R.; Dzombak, D.A.; Lowry, G.V.; Thaulow, N. Degradation of Well Cement by CO2 under Geologic Sequestration Conditions. Environ. Sci. Technol. 2007, 41, 4787–4792. [Google Scholar] [CrossRef] [PubMed]
  31. Huet, B.M.; Prevost, J.H.; Scherer, G.W. Quantitative Reactive Transport Modeling of Portland Cement in CO2-Saturated Water. Int. J. Greenh. Gas Control. 2010, 4, 561–574. [Google Scholar] [CrossRef]
  32. Hernández-Rodríguez, A.; Orlando, A.; Montegrossi, G.; Huet, B.; Virgili, G.; Vaselli, O. Experimental Analysis on the Carbonation Rate of Portland Cement at Room Temperature and CO2 Partial Pressure from 1 to 51 Bar. Cem. Concr. Compos. 2021, 124, 104271. [Google Scholar] [CrossRef]
  33. Kampman, N.; Busch, A.; Bertier, P.; Snippe, J.; Hangx, S.; Pipich, V.; Di, Z.; Rother, G.; Harrington, J.F.; Evans, J.P.; et al. Observational Evidence Confirms Modelling of the Long-Term Integrity of CO2-Reservoir Caprocks. Nat. Commun. 2016, 7, 12268. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Cadogan, S.P.; Maitland, G.C.; Trusler, J.P.M. Diffusion Coefficients of CO2 and N2 in Water at Temperatures between 298. 15 K and 423.15 K at Pressures up to 45 MPa. J. Chem. Eng. Data 2014, 59, 519–525. [Google Scholar] [CrossRef] [Green Version]
  35. Wollenweber, J.; Alles, S.a.; Kronimus, A.; Busch, A.; Stanjek, H.; Krooss, B.M. Caprock and Overburden Processes in Geological CO2 Storage: An Experimental Study on Sealing Efficiency and Mineral Alterations. Energy Procedia 2009, 1, 3469–3476. [Google Scholar] [CrossRef] [Green Version]
  36. Busch, A.; Alles, S.; Gensterblum, Y.; Prinz, D.; Dewhurst, D.N.; Raven, M.D.; Stanjek, H.; Krooss, B.M. Carbon Dioxide Storage Potential of Shales. Int. J. Greenh. Gas Control. 2008, 2, 297–308. [Google Scholar] [CrossRef]
  37. Busch, A.; Alles, S.; Krooss, B.M.; Stanjek, H.; Dewhurst, D. Effects of Physical Sorption and Chemical Reactions of CO2 in Shaly Caprocks. Energy Procedia 2009, 1, 3229–3235. [Google Scholar] [CrossRef] [Green Version]
  38. Wang, J.G.; Peng, Y. Numerical Modeling for the Combined Effects of Two-Phase Flow, Deformation, Gas Diffusion and CO2 Sorption on Caprock Sealing Efficiency. J. Geochem. Explor. 2014, 144, 154–167. [Google Scholar] [CrossRef]
  39. Wang, Z.; Li, Y.; Liu, H.; Zeng, F.; Guo, P.; Jiang, W. Study on the Adsorption, Diffusion and Permeation Selectivity of Shale Gas in Organics. Energies 2017, 10, 142. [Google Scholar] [CrossRef]
Figure 1. (a) Simplified stratigraphic log of the selected structure with cores location and housing. Depths (m) are referred to sea level. (b) Geothermal and geobaric gradient as measured in drilled wells; (c) carbonate C3 sample as its original form (larger block) and related charge after experiment (smaller cube); (d) shale C1 sample after 2 days of experiment; (e) shale C1 sample after 5 days of experiment.
Figure 1. (a) Simplified stratigraphic log of the selected structure with cores location and housing. Depths (m) are referred to sea level. (b) Geothermal and geobaric gradient as measured in drilled wells; (c) carbonate C3 sample as its original form (larger block) and related charge after experiment (smaller cube); (d) shale C1 sample after 2 days of experiment; (e) shale C1 sample after 5 days of experiment.
Minerals 13 00056 g001
Figure 2. Workflow of diffusion experiments. Carbonate specimens were cut in cubes, whereas shales were cored in cylinders and hosted in a custom-made reaction vessel in high density polytetrafluoroethylene (PTFE).
Figure 2. Workflow of diffusion experiments. Carbonate specimens were cut in cubes, whereas shales were cored in cylinders and hosted in a custom-made reaction vessel in high density polytetrafluoroethylene (PTFE).
Minerals 13 00056 g002
Figure 3. XRPD spectra and Rietveld refinement of rock samples. (a) C1 shale; (b) C2 shale; (c) C3 carbonate. Quantitative data results are in Table 1.
Figure 3. XRPD spectra and Rietveld refinement of rock samples. (a) C1 shale; (b) C2 shale; (c) C3 carbonate. Quantitative data results are in Table 1.
Minerals 13 00056 g003
Figure 4. FT-IR spectra (a) in the 3850–2780 cm−1 and (b) in the 1780–470 cm−1 range of C1, C2 and C3 granular samples pressed under the ATR module with indication of main signals (vertical dotted lines with arrows). Aliquots after experiments are those reporting in sample names “2d” and “5d” for 2 days and 5 days of hydrothermal interaction, respectively. Montmorillonite and muscovite are from the USGS database. Upper boxes indicate the absorption bands for discriminant minerals.
Figure 4. FT-IR spectra (a) in the 3850–2780 cm−1 and (b) in the 1780–470 cm−1 range of C1, C2 and C3 granular samples pressed under the ATR module with indication of main signals (vertical dotted lines with arrows). Aliquots after experiments are those reporting in sample names “2d” and “5d” for 2 days and 5 days of hydrothermal interaction, respectively. Montmorillonite and muscovite are from the USGS database. Upper boxes indicate the absorption bands for discriminant minerals.
Minerals 13 00056 g004
Figure 5. FT-IR spectra of C1, C2 and C3 samples in the 7400–5000 cm−1 range by Integrating Sphere module with indication of main signals (vertical dotted lines with arrows). C1 and C3 are the original shale and carbonate samples, respectively; the remaining sample names indicate aliquots after experiments (2d and 5d for 2 days and 5 days of hydrothermal interaction, respectively). C1, C3 and C3_2d were analyzed as blocks.
Figure 5. FT-IR spectra of C1, C2 and C3 samples in the 7400–5000 cm−1 range by Integrating Sphere module with indication of main signals (vertical dotted lines with arrows). C1 and C3 are the original shale and carbonate samples, respectively; the remaining sample names indicate aliquots after experiments (2d and 5d for 2 days and 5 days of hydrothermal interaction, respectively). C1, C3 and C3_2d were analyzed as blocks.
Minerals 13 00056 g005
Figure 7. SEM-EDS data of C2 sample in thin section. (a) General aspect at low magnification scale showing a dominant fine clayey matrix that includes larger crystals at the analytical Site 1. S1: pyrite; S2–S4: quartz. (b) Texture and mineralogy of the analytical Site 2. S1: dolomite-siderite; S2: illite, S3: anatase; S4: dolomite-siderite. (c) Texture and mineralogy of the analytical Site 3. S1: kaolinite; S2: dolomite-siderite; S3: muscovite; S4: illite-chlorite. (d) EDS spectrum of Site 1 in panel (a), S1: pyrite. Mineral compositions are related to spectrum for site in Table 4.
Figure 7. SEM-EDS data of C2 sample in thin section. (a) General aspect at low magnification scale showing a dominant fine clayey matrix that includes larger crystals at the analytical Site 1. S1: pyrite; S2–S4: quartz. (b) Texture and mineralogy of the analytical Site 2. S1: dolomite-siderite; S2: illite, S3: anatase; S4: dolomite-siderite. (c) Texture and mineralogy of the analytical Site 3. S1: kaolinite; S2: dolomite-siderite; S3: muscovite; S4: illite-chlorite. (d) EDS spectrum of Site 1 in panel (a), S1: pyrite. Mineral compositions are related to spectrum for site in Table 4.
Minerals 13 00056 g007
Figure 8. SEM-EDS data of C3 sample in thin section. (a) General aspect at low magnification scale at the analytical Site 2. S1: pyrite; S2 and S3: primary calcite; S4: kaolinite-illite. (b) Texture and mineralogy of the Site 1. S1: pyrite; S2: calcite; S3: kaolinite-illite; S4: quartz. (c) Texture and mineralogy of the Site 3. S1: kaolinite; S2: illite; S3: pyrite; S4: muscovite. (d) EDS spectrum of Site S4 in panel (c) indicating mixture of kaolinite-illite chemistry. Mineral compositions are related to spectrum for site in Table 5.
Figure 8. SEM-EDS data of C3 sample in thin section. (a) General aspect at low magnification scale at the analytical Site 2. S1: pyrite; S2 and S3: primary calcite; S4: kaolinite-illite. (b) Texture and mineralogy of the Site 1. S1: pyrite; S2: calcite; S3: kaolinite-illite; S4: quartz. (c) Texture and mineralogy of the Site 3. S1: kaolinite; S2: illite; S3: pyrite; S4: muscovite. (d) EDS spectrum of Site S4 in panel (c) indicating mixture of kaolinite-illite chemistry. Mineral compositions are related to spectrum for site in Table 5.
Minerals 13 00056 g008
Figure 9. SEM images of sample C1 after 2 days of experiments, oriented with respect to the CO2 entry. (a) Low magnification view. The two white lines represent, from left to the right, the sample outer edge with CO2 richer fluids, the pyrite and other Fe-rich mineral dissolution and the secondary minerals reaction front and, on the right the zone of deposition. (b) Zoom on the matrix with the formation of Fe-richer carbonate (siderite) as a secondary mineral. Images taken from a piece on a stub.
Figure 9. SEM images of sample C1 after 2 days of experiments, oriented with respect to the CO2 entry. (a) Low magnification view. The two white lines represent, from left to the right, the sample outer edge with CO2 richer fluids, the pyrite and other Fe-rich mineral dissolution and the secondary minerals reaction front and, on the right the zone of deposition. (b) Zoom on the matrix with the formation of Fe-richer carbonate (siderite) as a secondary mineral. Images taken from a piece on a stub.
Minerals 13 00056 g009
Figure 10. SEM-EDS results of C1, C2 and C3 samples. (a) Cameo image of C1 sample after 2 days of experiments. RGB coloring for Fe, Ca and Mg, respectively. (b) Cameo image of C2 sample after 5 days of experiments. RGB coloring for K, Mg and Fe, respectively. (c) Cameo image of C3 sample after 5 days of experiments. Fe abundance in green color. Colors distribution highlight the reaction front (yellow line). Images taken from a piece on a stub.
Figure 10. SEM-EDS results of C1, C2 and C3 samples. (a) Cameo image of C1 sample after 2 days of experiments. RGB coloring for Fe, Ca and Mg, respectively. (b) Cameo image of C2 sample after 5 days of experiments. RGB coloring for K, Mg and Fe, respectively. (c) Cameo image of C3 sample after 5 days of experiments. Fe abundance in green color. Colors distribution highlight the reaction front (yellow line). Images taken from a piece on a stub.
Minerals 13 00056 g010
Figure 11. Sample C3, after 5 days of experiments. Details of the depletion zone with carbonate with higher abundance of Fe and Mg (and/or Mn). Images taken from a piece on a stub.
Figure 11. Sample C3, after 5 days of experiments. Details of the depletion zone with carbonate with higher abundance of Fe and Mg (and/or Mn). Images taken from a piece on a stub.
Minerals 13 00056 g011
Figure 12. FT-IR spectra zoomed in (a) in the 1500–1040 cm−1 and (b) in the 920–800 cm−1 range of C3 samples (c) with textural changes between pristine (in red) and different weathered portions (blue, green and magenta lines) after 2 days of experiments. Bands displacement from [23].
Figure 12. FT-IR spectra zoomed in (a) in the 1500–1040 cm−1 and (b) in the 920–800 cm−1 range of C3 samples (c) with textural changes between pristine (in red) and different weathered portions (blue, green and magenta lines) after 2 days of experiments. Bands displacement from [23].
Minerals 13 00056 g012
Table 1. Semi-quantitative mineralogical composition of selected cores (weight %) after Rietveld refinement, normalized to 100.
Table 1. Semi-quantitative mineralogical composition of selected cores (weight %) after Rietveld refinement, normalized to 100.
MineralsSeal 6Seal 4Reservoir
C1C2C3
Quartz35.9631.2710.04
Siderite4.299.640.64
Dolomite2.053.1611.63
Calcite1.071.0763.94
Ankerite--<1
Plagioclase14.984.322.05
Muscovite10.8310.781.57
Anatase<1--
Pyrite<1<1<1
Chlorite3.041.161.59
Kaolinite5.4314.144.70
Illite22.3524.463.85
-: not present.
Table 2. XRF analysis (oxides weight %) of C1, C2 and C3 rock samples.
Table 2. XRF analysis (oxides weight %) of C1, C2 and C3 rock samples.
OxidesC1C2C3C1-Corrected C2-CorrectedC3-Corrected
B2O32.41- 2.111.61-1.41
Na2O0.810.230.250.540.150.16
MgO1.401.100.810.930.740.54
Al2O319.0922.898.9012.7715.315.95
SiO261.4957.3119.5741.1238.3213.09
P2O50.150.150.110.100.100.07
SO30.500.901.990.340.601.33
Cl0.040.020.240.030.010.16
K2O3.854.251.032.572.840.69
CaO0.821.4160.050.550.9440.15
TiO21.281.130.420.850.760.28
V2O5-- 0.03--0.02
Cr2O30.050.030.020.030.020.01
MnO0.120.090.030.080.060.02
Fe2O37.7610.252.395.196.851.60
NiO0.0130.020.010.010.010.01
CuO0.010.010.010.010.060.01
ZnO0.020.020.010.010.02<0.01
Rb2O0.030.030.01--<0.001
Ga2O3-0.01----
As2O30.010.01----
SrO0.030.040.600.0170.030.40
ZrO20.050.03----
Nb2O5-<0.01----
BaO0.090.071.400.060.050.93
PbO--0.014--0.01
CO2---3.556.5534.716
Corrected: the compositions were corrected for the contribution of CO2 loss from carbonates determined by calcimetric analysis.
Table 4. Semi-quantitative analysis of C2 sample in atoms percentage.
Table 4. Semi-quantitative analysis of C2 sample in atoms percentage.
SpectrumMineralONaMgAlSiKSCaTiFe
2-Site1quartz68.48---31.52-----
3-Site1quartz68.45--0.39 *31.16-----
4-Site1quartz66.96---33.04-----
1-Site2dolomite-siderite74.03-9.901.16 *1.76 *0.20 *-2.06-10.90
2-Site2illite68.570.710.3011.8415.832.39---0.36
3-Site2anatase73.52--4.78 *9.44 *0.65 *--11.340.26 *
4-Site2dolomite-siderite73.38-9.691.19 *1.29 *--2.11-12.34
1-Site3kaolinite68.560.42 *0.54 *11.1117.001.86 *---0.51 *
2-Site3dolomite-siderite74.49-8.402.42 *3.43 *--1.86-9.41
3-Site3muscovite68.340.44 *0.31 *8.8019.871.56--0.210.46
4-Site3illite-chlorite70.14-2.477.9514.481.37-0.510.252.83
*: Contamination from close minerals. Refer to Figure 6 for details on spectrum and site number.
Table 5. Semi-quantitative analysis of C3 sample in atoms percentage.
Table 5. Semi-quantitative analysis of C3 sample in atoms percentage.
SpectrumMineralONaMgAlSiKPSCaClTiFe
2-Site1calcite79.92--1.85 *2.85 *0.32 *-0.28*14.79---
3-Site1kaolinite - illite74.550.49 *0.39 *7.0614.591.10 *--1.45 *0.18 *-0.19 *
4-Site1quartz73.69-0.27 *1.4323.490.20 *--0.92 *---
2-Site2calcite77.71---0.34 *---21.94---
3-Site2calcite74.84-------24.860.30 *--
4-Site2kaolinite-illite69.350.45 *0.59 *11.0114.472.06-0.46 *1.31--0.31
1-Site3kaolinite 75.92--9.5313.070.20 *--1.17 *--0.11
2-Site3illite 71.240.270.519.7514.881.350.28-1.48--0.23
4-Site3muscovite72.98-0.877.8813.072.10--2.38 *-0.260.46
*: Contamination from close minerals. Refer to Figure 6 for details on spectrum and site number.
Table 6. Diffusion coefficients obtained from experiments.
Table 6. Diffusion coefficients obtained from experiments.
Temperature, °CPenetration Length, mm Depletion Length, mmDraw, m2s−1De, m2s−1
C1
2 days752.99 ± 0.010.67 ± 0.015.17 ± 0.035 × 10−119.77 × 10−11
5 days753.48 ± 0.011.48 ± 0.012.80 ± 0.016 × 10−11
C2
2 days1052.60 ± 0.010.66 ± 0.013.91 ± 0.03 × 10−111.53 × 10−10
5 days1053.50 ± 0.011.24 ± 0.012.84 ± 0.016 × 10−11
C3
2 days1452.41 ± 0.011.44 ± 0.013.36 ± 0.028 × 10−112.51 × 10−10
5 days1453.67 ± 0.010.94 ± 0.013.12 ± 0.017 × 10−11
Table 7. Diffusion coefficients for shales and marlstones reported in literature.
Table 7. Diffusion coefficients for shales and marlstones reported in literature.
T, °CDe, m2 s−1Reference
503.08 × 10−11[36]
504.81 × 10−11
28 7.80 × 10−11 − 1.2 × 10−10[35]
1.20 × 10−11[38]
203 × 10−11[33]
4–8 × 10−11[39]
T, °CD_CO2(aq), m2 s−1Reference
252.233 × 10−09[34]
503.718 × 10−09
755.391 × 10−09
1007.521 × 10−09
15012.33 × 10−09
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Montegrossi, G.; Cantucci, B.; Piochi, M.; Fusi, L.; Misnan, M.S.; Rashidi, M.R.A.; Abu Bakar, Z.A.; Tuan Harith, Z.Z.; Bahri, N.H.S.; Hashim, N. CO2 Reaction-Diffusion Experiments in Shales and Carbonates. Minerals 2023, 13, 56. https://doi.org/10.3390/min13010056

AMA Style

Montegrossi G, Cantucci B, Piochi M, Fusi L, Misnan MS, Rashidi MRA, Abu Bakar ZA, Tuan Harith ZZ, Bahri NHS, Hashim N. CO2 Reaction-Diffusion Experiments in Shales and Carbonates. Minerals. 2023; 13(1):56. https://doi.org/10.3390/min13010056

Chicago/Turabian Style

Montegrossi, Giordano, Barbara Cantucci, Monica Piochi, Lorenzo Fusi, M. Shahir Misnan, M. Rashad Amir Rashidi, Zainol Affendi Abu Bakar, Zuhar Zahir Tuan Harith, Nabila Hannah Samsol Bahri, and Noorbaizura Hashim. 2023. "CO2 Reaction-Diffusion Experiments in Shales and Carbonates" Minerals 13, no. 1: 56. https://doi.org/10.3390/min13010056

APA Style

Montegrossi, G., Cantucci, B., Piochi, M., Fusi, L., Misnan, M. S., Rashidi, M. R. A., Abu Bakar, Z. A., Tuan Harith, Z. Z., Bahri, N. H. S., & Hashim, N. (2023). CO2 Reaction-Diffusion Experiments in Shales and Carbonates. Minerals, 13(1), 56. https://doi.org/10.3390/min13010056

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop