Experimental and Modelled Reactions of CO 2 and SO 2 with Core from a Low Salinity Aquifer Overlying a Target CO 2 Storage Complex

: CO 2 ‐ induced reactions in low salinity aquifers overlying CO 2 storage sites are of interest to understand potential reactions or impacts in the possible case of a leak. Previous investigations of overlying aquifers in the context of CO 2 storage have focused on pure CO 2 streams, however captured industrial CO 2 streams may contain ancillary gases, including SO 2 , O 2 , NOx, H 2 S, N 2 , etc., some of which may be more reactive than CO 2 when dissolved in formation water. Eight drill cores from two wells in a low salinity sandstone aquifer that overlies a target CO 2 storage complex are characterised for porosity (helium, mercury injection, or micro CT), permeability, and mineral content. The eight Hutton Sandstone cores are variable with porosities of 5.2–19.6%, including carbonaceous mudstones, calcite cemented sandstones, and quartz rich sandstones, common lithologies that may be found generally in overlying aquifers of CO 2 storage sites. A chlorite rich sandstone was experimentally reacted with CO 2 and low concentrations of SO 2 to investigate the potential reactions and possible mineral trapping in the unlikely event of a leak. Micro CT characterisation before and after the reaction indicated no significant change in porosity, although some fines movement was observed that could affect permeability. Dissolved concentrations of Fe, Ca, Mn, Cr, Mg, Rb, Li, Zn, etc., increased during the reaction, including from dissolution of chlorite and trace amounts of ankerite. After ~40 days dissolved concentrations including Fe, Zn, Al, Ba, As and Cr decreased. Chlorite was corroded, and Fe ‐ rich precipitates mainly Fe ‐ Cr oxides were observed to be precipitated on rock surfaces after experimental reaction. Concentrations of Rb and Li increased steadily and deserve further investigation as potential monitoring indicators for a leak. The reaction of chlorite rich sandstone with CO 2 and SO 2 was geochemically modelled over 10 years, with mainly chlorite alteration to siderite mineral trapping 1.55 kg/m 3 of CO 2 and removing dissolved Fe from solution. Kaolinite and chalcedony precipitation was also predicted, with minor pyrite precipitation trapping SO 2, however no changes to porosity were predicted.


Introduction
Owing to international interest in existing and potential CO2 storage sites in deep saline aquifers, the majority of experimental CO2-water-rock reactions have been performed in brines and with pure CO2. There have however been international studies of CO2-water-rock interactions of core from low salinity or fresh water aquifers overlying CO2 storage sites to understand any potential changes in the unlikely event of a leak. Examples of field studies relevant to low salinity aquifer conditions are stock use. In addition the Hutton Sandstone underlies the Walloon coal measures, an important coal seam gas production interval [17,18]. Differences in the groundwater hydrochemistry have been reported, as being fresh in the north and western outcrops with higher salinity to the south and east, especially east of the Burunga Leichhardt fault zone [11,19]. Hydrochemical analyses of groundwater have indicated that a dual permeability is likely in the Hutton Sandstone of the northern Surat Basin near the Mimosa syncline, with groundwater flow limited to ~50 m of its total thickness [20,21]. There is also reported evidence for possible groundwater movement from the Precipice Sandstone to the Hutton Sandstone up the Hutton Wallumbilla fault, and Burunga Leichhardt fault zone [20,21]. Sequence stratigraphy and also non-quantitative well core scale mineral identification have been performed by Hylogger on several wells including West Wandoan 1, GSQ Chinchilla 4, and Woleebee Creek GW4 [22].
Here we quantitatively characterize in detail sandstone, mudstone and calcite cemented cores from two wells in a low salinity aquifer, the Hutton Sandstone, overlying a potential storage reservoir and caprock. One core is characterized before and after experimental reaction with supercritical CO2, SO2 and low salinity water by automated mineral quantification (QEMSCAN), scanning electron microscopy (SEM-EDS) and micro Computed Tomography (CT) to determine potential changes to minerals and porosity. Dissolved metals released in the reaction waters were determined, and longer term geochemical modelling performed to predict potential for mineral trapping of CO2.

Materials and Methods
Drill cores were sampled from the Hutton Sandstone of the GSQ Chinchilla 4 well (latitude -26.72722, and longitude of 150.2014 decimal degrees, approximately 10 km SSE of Miles), and the West Wandoan 1 (WW1) well (latitude -26.181622, longitude 149.812422, approximately 19 km south west of the town of Wandoan), in Queensland, Australia. The stratigraphic column and a map showing the well locations in given in supplementary material Figure S1. The West Wandoan 1 well was drilled for a CO2 storage feasibility study of a demo scale injection into the Precipice Sandstone. Core air permeability, mercury intrusion porosimetry (MICP), and Helium porosity of the Chinchilla 4 well cores described here were performed at Weatherford. MICP, He pychnometry, permeability, and X-ray diffraction of West Wandoan 1 core was performed at the University of Queensland (UQ) by methods reported previously [10,23,24]. Scanning electron microscopy (SEM) and Energy Dispersive Spectroscopy (EDS) in back scatter detection mode (BSE) with both a low-vacuum JEOL 6460LA environmental SEM (Peabody, MA, USA), and Hitachi TM3030 (Tokyo, Japan) with a Bruker EDS (Billerica, MA, USA) was performed on polished thin sections and blocks. Core total acid digestions or lithium metaborate fusions and loss on ignition were performed in the UQ Environmental Geochemistry Laboratory.
Micro CT was performed on sub-plugs to visualise and calculate pore space. QEMSCAN is an automated mineral quantification based on SEM-EDS that also provides visualisation of mineral associations and grain sizes. These were performed on polished sub-plug sections by FEI Australia (Canberra, Australia) (and more recently at the Australian National University) as described in detail previously [11,25,26]. Briefly, for the WW1 724.1 m sample that was reacted, a 3 mm diameter subplug was digitally characterised in 3D by X-ray micro computed tomography (μCT) with a voxel size of 2.2 µm. The 3D μCT images were registered into perfect geometric alignment with higherresolution 2D SEM images and automated quantified SEM-EDS (QEMSCAN®) mineral maps were produced of a polished sub-plug slice from a trimmed end. The sub-plug was reacted, before being imaged again after reaction, and the two sets of before and after reaction images registered to one another to characterise the changes.
Long term batch experiments were performed in Parr reactors (Moline, IL, USA) at 120 bar and 60 °C for approximately 10 weeks with the reactors described in detail previously [12]. A low salinity water (100 ml of 1500 ppm NaCl) was added with the WW1 724.1 m rock core (15 mm cube, sub-plug and offcut) at a water:rock ratio of 7 by mass. Data on the in situ formation water chemistry was not available, however measurements in other regions of the aquifer range from ~ 282-1863 mg/kg NaCl with variable bicarbonate alkalinity ~136-733 mg/kg [19]. The simplified synthetic formation water chemistry was chosen to both be in the range of possible formation water salinities, and to be consistent with previous studies for comparison. Reactors were purged with N2 and pressurised to 120 bar at 60 °C for 5 days to provide a baseline water-rock soak chemistry. This time was similar to previous work given time constraints, however it is possible that full equilibration of water and rock may not have occurred. After fluid sampling, N2 was depressurised and reactors filled with 160 ppm SO2 in a balance of CO2. After another 18 days sodium bicarbonate was added (to 105 ppm) to approximate reservoir buffering. Note that this likely introduced some O2 from air contact, however O2 or air may be present anyway in injected CO2 streams especially from oxyfuel firing. Fluid was sampled during the experiments, with pH and electrical conductivity immediately measured with a TPS WP81 meter and probes. Aliquots were diluted 20 times and acidified to 2% ultra-pure nitric acid for analysis of ions in the UQ Environmental Geochemistry Laboratory by ICP-OES (Optima 3300 DV ICP-OES, Perkin Elmer, Waltham, MA, USA) with an error of ~5% for major elements. Trace element concentrations were measured by ICP MS (7900 ICP-MS, Agilent, Santa Clara, CA, USA with collision cell) with errors of less than 10%. Total inorganic and organic carbon, alkalinity, and sulphate were determined on selected unacidified samples (total organic carbon analyser, and ion chromatography performed at ALS environmental), and also NOx and phosphate at the end of the reaction (UQ, Lachat QuikChem8500 Flow Injection Analyser). A blank experiment without core was also performed to determine if cations were leached from the reactor.
Kinetic geochemical models were constructed for the reaction of the WW1 724.1 m Hutton Sandstone for up to 10 years from the characterization data, with the input minerals given in Supplementary material. The upscaled reactive surface areas of minerals used are also given in supplementary material, with the water-rock ratio based on the porosity and an equilibrated water chemistry based on formation water measurements [19]. The general methods and mineral kinetic and thermodynamic parameters have been published in detail previously for other rock reactions [11,12,[26][27][28][29]. Briefly, geochemical models were run in the react module of Geochemist work Bench (GWB) version 9, using the EQ3/6 database, with minerals input via mineral script files to describe acid neutral and basic mechanisms [29][30][31]. CO2 fugacity was calculated at 12 MPa and 60 °C from Duan and Sun (2003), with SO2 gas added by mass, models were also run with the CO2 fugacity at half of the full fugacity to test the effect on pH [32]. Saturated minerals were allowed to precipitate based on observations of experiments and natural analogue systems, e.g., the carbonates siderite and ankerite/dolomite have been observed to precipitate in natural systems along with kaolinite, smectites, silica, and pyrite [33][34][35][36].

Core Characterisation
The Hutton Sandstone cores from the two wells show a variety of pore throat distributions ( Figure 1A). The Chin 4828.76 m core with large pore throats corresponds to the coarsest grained sandstone with visible open pores (Figure 2A), and Chin 4835.48 m with a relatively high clay content has the majority in the range 0.01-0.1 μm ( Figure 1A, Figure 2B, Figure 3). The Chin 4828.76 m sandstone corresponds to a region with high measured air permeabilities up to 1228 mD ( Figure 1B). Table 1 compares the corresponding MICP porosities with those by μCT or helium (He) porosities (along with Table 2). The He porosities tend to give the highest estimates likely owing to the small size of the He molecule, with MICP and μCT values generally in reasonable agreement with each other. Selected µCT tomograms and SEM images of the Chinchilla 4 well cores are shown in Figure 2 and Figure 3. The 799.5 m core for example is calcite cemented, the 835.48 m core contains 10% chlorite which along with kaolinite fills and rims pores. QEMSCAN images and quantified mineral components by QEMSCAN or XRD are given in Supplementary material Figure S2, Table S1 and  Table S2, along with the volume % of X-ray distinct components Table S3. The WW1 cores are equally variable in lithology, including chlorite rich and calcite cemented sandstones, two also contain coal ( Figure 4, and Supplementary Material Figure S3).

Hutton Sandstone WW1 724.1 m Before and After Reaction
The chlorite rich sandstone WW1 724.1 m was reacted with CO2 and SO2 since chlorite is known to alter to siderite and ankerite mineral trapping CO2 in natural systems. The µCT image of the subplug before and after reaction is shown in Figure 5 along with the mineral segmentation image. The porosity did not change measurably after reaction (Table 2). SEM and QEMSCAN of a sub-plug slice before reaction show that chlorite tended to rim porosity, the sub-plug also contained ~7% plagioclase and K-feldspar ( Figure 6, Figure 7). After reaction the mineral content by QEMSCAN did not appear to change significantly (Figure 7). Some movement of fine material was however observed in pores ( Figure 7, and supplementary material Figure S4, S5, S6). SEM-EDS of a core block before and after reaction is shown in Figure 8. Along with the minerals identified in QEMSCAN of the sub-plug, other minerals present included Fe-Mg or Fe-Ti altered micas, phosphates containing rare earth elements, and coal. Zircon crystals were present in a band through the core, and chlorite was Fe-rich also containing Mg and Mn, S signatures were present in illite. After reaction only alteration to chlorite surfaces and loss of Fe and Ca signatures from chlorite and trace calcite/ankerite on illite were observed, with precipitation of Fe-Cr-Ni-oxides on chlorite. Rock mass decreased only slightly from 13.68 to 13.62 g after reaction.
Whole rock digest data for total metal content in the core is given in supplementary material Table S5

Water Chemistry during Reaction of WW1 724.1 m
The measured ex situ pH decreased from 7.5 to 5.2 after the rock -water soak period, and then to 4.8 after CO2-SO2 addition to reactors. The pH then increased slightly to 5.41 by the end of the reaction (Figure 9). The measured electrical conductivity decreased slightly from 2.79 to 2.48 ms/cm and increased to 3.57 ms/cm by the end of the experiment (Supplementary Material Figure S7). The pH varied from 3.54 to 3.93 during the blank experiment without rock, and conductivity from 2.09 to 2.59 ms/cm indicating some pH buffering by minerals in the experiment containing the rock core. After CO2 -SO2 injection, several dissolved elements increased including Ca, Mn, Mg, Fe, Cr, As, Pb, Rb, Ti, Tl, Al, Zn, K, Si, S (Figure 9, Figure 10, Supplementary Material Figure S7, S8). Fe-Mg-(Mn)chlorite corrosion was observed directly in SEM images, and contributed to dissolved Fe, Mg, Mn, Si, and Al. Ni, Zn, Ti and Li can also substitute into the chlorite structure and may have been partly sourced from chlorite dissolution. Dissolution of trace amounts of calcite or ankerite/siderite would contribute to the Ca, Mg, Mn, and Fe. Tl has been shown elsewhere to substitute into and be sourced from carbonate mineral dissolution in CO2-water-rock reactions of calcite or dolomite [3]. The gradual increases in K and Rb indicate minor corrosion of feldspars or illite continuing over the reaction timescale. Several elements including Fe, As, Zn, Al, Ca, had a decreasing trend after ~1000 h, with dissolved concentrations of Cr, Pb, Al, Ba, decreasing significantly before ~500 h. These were likely either adsorbed or precipitated onto rock surfaces e.g., in the precipitated oxide minerals observed or as surface coatings. Barite precipitation was not directly observed however its low solubility, the decreasing dissolved Ba, and the presence of dissolved S indicates its likely precipitation. Barite precipitation has been directly observed in other CO2-SO2-water rock reactions of calcite cemented sandstone where higher concentrations of SO2 (0.16%) were used [25,37]. The rock surfaces had a brown coloration after reaction supporting precipitation of Fe-containing minerals as surface coatings that may not have been visible in SEM images (Supplementary Material Figure S9). Dissolved Ni was variable and appeared to increase at the end of reaction with the blank experiment indicating some potential contribution from the reactor. Dissolved Cr concentration at the end of the blank experiment was higher than the experiment with the rock present also indicating some contribution from the reactor. The reason for the lower Cr concentration at the end of the experiment with rock present may be owing to the higher pH (therefore less reactor corrosion), or likely the precipitation of Cr containing minerals seen on the rock surface. The dissolved concentrations of Pb, U and As were relatively low at less than 8, 1 and 1 µg/kg, respectively, during the reaction.  Total S was 4.6 mg/kg after the N2-low salinity water-rock soak indicating some trace pyrite or sphalerite may have reacted or adsorbed S released, this may have also contributed to the decrease in pH. Sphalerite is acid reactive, and reaction would also have contributed to the increase in dissolved Zn and metals such as Pb and As. While it was not directly observed in 724.1 m core here, trace amounts of sphalerite have been observed sporadically in other Hutton Sandstone sections e.g., 800.7 m. Total S increased to 14.4 after CO2-SO2 injection and reached 18.5 mg/kg at the end of reaction. Alkalinity was in the bicarbonate form at 358 mg/kg, and Cl was 845 mg/kg on experiment termination. This is lower that might be expected from the initial salinity and indicates either minor salt precipitation, or adsorption onto clays and surfaces. Dissolved total carbon, total organic carbon, inorganic carbon and sulphate measured during reaction are shown in Supplementary Material Figure S7.

Geochemical Model
The kinetic geochemical model prediction of reaction of Hutton Sandstone with a mineral composition of the 724.1 m core is shown over 10 years in Figure 11. The alteration of chlorite to siderite, kaolinite and chalcedony is mainly predicted. Dissolution of albite, K-feldspar, siderite, ankerite, calcite and precipitation of pyrite and ankerite is also predicted. Pyrite was predicted to precipitate in Hutton Sandstone WW1 724.1 m, sourcing S from the co-injected SO2. In the current reaction of Hutton Sandstone WW1 724.1 m the predicted altering mineral volumes are relatively small, with only 25 cm 3 of (original 1279 cm 3 Supplementary Material Table S4) chlorite dissolved, and 13 cm 3 of siderite precipitated with 12 cm 3 of kaolinite, and with no net change in volume (or porosity). A pH of 5.2 was predicted after 10 years, with the concentration of dissolved Fe only 4.18 mg/kg as Fe was mainly sequestered with CO2 as siderite. Net mineral trapping was predicted as 1.55 kg/m 3 CO2 as siderite after 10 y (based on the method of Watson and Gibson-Poole ) [38]. Additional modelling scenarios for Hutton Sandstone WW1 724.1 m allowing precipitation of smectite are shown in Supplementary material ( Figure S10), where precipitation of smectite occurs mainly replacing kaolinite precipitation.

Discussion
Overall the drill cores characterized from the Hutton Sandstone varied in mineral content, metal content, porosity, pore throat sizes, and permeability. Air permeability measurements of the Chinchilla 4 Hutton Sandstone core showed a high permeability section of ~50 m which is in agreement with the findings of [20,21], that suggested groundwater flow may be restricted to a ~50 m section of the Hutton Sandstone near the Mimosa syncline. Here a chlorite rich Hutton Sandstone, WW1 724.1 m was reacted with CO2 and SO2 in a low salinity water. Porosity measured by micro CT did not change measurably on reaction, although fines movement was observed. No increase to porosity would be favorable to avoid further leakage or migration of CO2. In contrast, an increase to porosity of 1.1% on pure CO2 reaction of calcite cemented Chinchilla 4 Hutton Sandstone 799 m at 60 °C was reported elsewhere [8]. In that case, however, dissolution of 2 vol% calcite cementing framework grains had caused the increased porosity. Two previous studies have reported the reaction of a calcite cemented Hutton Sandstone core from WW1 800.83 m with 1500 mg/kg NaCl and CO2 or CO2 containing 0.16% SO2 and 2% O2 at 60 °C [11,12]. They also reported calcite cement dissolution, and a porosity increase from 5.2 to 11.3% after CO2 reaction; or from 3.8 to 7.5% after CO2-SO2-O2 reaction. These studies indicate that the lithology of the section of the overlying aquifer receiving a potential leak may influence its extent of propagation, with carbonate cemented lithologies a higher potential risk than chlorite rich sandstone.
Precipitation of gypsum or anhydrite had been observed or predicted in CO2 reactivity studies with SO2 ± O2 reaction of rock containing calcite cements as a source of calcium [27,[39][40][41][42]. The sandstone reacted here however had a small amount of calcium containing minerals, so a lower dissolved Ca concentration, therefore gypsum/anhydrite did not reach saturation and was not predicted to precipitate. In addition the experiment performed here was at 60 °C, below the stability for anhydrite. A feldspar rich sandstone from the Chinchilla 4 well 868 m was also recently reacted for 16 days with water and CO2 containing 0.16% SO2 and 2% O2, with reported corrosion of siderite and chlorite, precipitation of Fe-oxides with Cr and Ni signatures, and an Al-sulphate mineral jarosite [29]. While Fe-oxides were observed to precipitate in the current reaction reported here, no sulphate mineral precipitation was observed, in contrast FeS/pyrite was observed and predicted to form owing to the lack of co-injected O2. Precipitation of pyrite was observed elsewhere in a natural analogue system of CO2 and S co-sequestration in the Madison Limestone of the Moxa Arch, Wyoming, at current reservoir temperatures of 90-110 °C [34,43]. There a high natural CO2-rich reservoir with dissolved aqueous sulphur complexes contained pyrite, native sulphur, and anhydrite cementing pore space in the drill core. These were shown by the authors to be an example endpoint for CO2 and S co-sequestration. Anhydrite precipitation could be expected to dominate over gypsum precipitation at those higher temperatures (above 60 °C) in the presence of dissolved Ca from the dolomite, and dissolved S from H2S. The prediction in the work presented here of pyrite precipitation with SO2 coinjection is in agreement with that work, although the temperature of the current study is much lower.
During the experimental reaction of the chlorite rich Hutton Sandstone WW1 724.1 m reported here for the reaction with CO2 and SO2, dissolved Fe was somewhat correlated with Mg (R 2 = 0.71), Fe was also correlated with Mn (R 2 = 0.82) as the majority of these were sourced from the chlorite dissolution. The correlations also reflect subsequent precipitation of metals including Fe and Mn, likely into Fe-oxyhydroxides. Fe was also somewhat correlated with Ba (R 2 = 0.73), and dissolved Ca with Fe (R 2 = 0.69). This indicates that the dissolution of trace amounts of ankerite also contributed to dissolved cations and was likely the source of Ba. The correlations of Fe with Ca and Ba also reflect later decreases in concentrations through precipitation indicating their incorporation into similar precipitating minerals. Dissolved Rb was somewhat correlated with Li (R 2 = 0.69) indicating a similar source mineral, likely plagioclase. The increase and subsequent decrease of dissolved Zn was also correlated with As (R 2 = 0.72) likely from trace amounts of sulphides dissolving and re-precipitating. Published experimental reactions of core overlying CO2 storage sites have mainly included calcite cemented sandstone, dolomite or limestones [2][3][4]8]. In those cases the majority of dissolved cations or metals were released from calcite (or dolomite) dissolution and reflected metals substituted into the calcite structure including Ba, Sr, Mn, Mg. The current study has shown that Fe, Mg and Mn are also released from reactive clays such as chlorite in overlying sandstone aquifers. The Fe, Mg and Mn are then available for mineral trapping over longer time periods. After 10 years, mineral trapping of CO2 was predicted to be 1.55 kg/m 3 in the form of siderite in models for our current study at 60 °C. Predictions of CO2 mineral trapping have been reported over longer time scales elsewhere. After 1000 y of CO2 reaction, the quartz rich Precipice Sandstone had a predicted mineral trapping from the reaction of trace amounts of chlorite to form siderite trapping only 1.24-1.30 kg/m 3 CO2 [44]. No mineral trapping was predicted before 30 y reaction in that case. Up to 2.57 kg/m 3 CO2 was predicted to be mineral trapped as siderite or ankerite after 30 y CO2 reaction of chlorite and plagioclase rich mudstones and sandstones of the Evergreen Formation caprock [44]. Watson and Gibson-Poole determined that 34.3 and 231.7 kg/m 3 of net CO2 was trapped as siderite, ankerite and calcite in the quartz-rich Waare Sandstone reservoir and the chlorite/berthierine-rich Flaxman Formation respectively [38]. These much higher mineral trapping amounts were estimated from observations of natural mineral trapping in a CO2-rich well vs a low CO2 well (natural analogue system) reacted over unconstrained geological time scales (thought to be since ~5 Ka-1 Ma), with the CO2 sourced from magmatic activity, and at higher present day temperature ranges ~ 75-116 °C. The main process mineral trapping CO2 was alteration of Fe-rich chlorite/bertherine clay to form siderite and ankerite. Clay minerals were additionally observed to have precipitated in that natural system, mainly kaolinite from alteration of plagioclase, in agreement with the current study, and additionally smectite/illite from alteration of K-feldspar. The predictions of mineral trapping of CO2 as siderite from reaction of Fe-rich chlorite clay in the Hutton Sandstone 724.1 m presented here are reasonable given the above and other studies of natural systems [33,36,[45][46][47]. A reactive transport modelling study of an arkose (20 vol% plagioclase), saline reservoir predicted Ca-Na-plagioclase (oligoclase) and chlorite alteration to ankerite, dawsonite, and siderite on reaction with CO2 and 1% SO2 at 75 °C. Predicted mineral trapping was 40-50 kg/m 3 over 1000-10,000 y, with dawsonite predicted to be formed from the Na supplied by plagioclase dissolution. SO2 was trapped as alunite, anhydrite and pyrite. In general, although the studies above were at different temperatures, mineral trapping as siderite and ankerite has been observed in reservoirs with Fe-bearing reactive minerals present e.g., Fe-rich chlorite; with dolomite or calcite precipitated from Mg and Ca-rich source minerals such as Mg-chlorite or Ca-plagioclase. Siderite and ankerite are able to precipitate at lower pH than dolomite and calcite, therefore if Fe is present in reducing conditions these ferroan carbonates tend to precipitate first [33]. In sandstone reservoirs, the presence of chlorite and Ca-rich plagioclase are generally associated with higher mineral trapping capacities. Precipitation of dawsonite has been predicted or observed in a few systems with, for example, high Na-plagioclase content and persistent high CO2 partial pressures.

Conclusions
This study has shown that low salinity aquifers overlying CO2 storage sites may be very variable in porosity, permeability and mineral content. The response to a possible CO2 leakage is therefore dependent on the lithologies receiving the leak. A chlorite rich sandstone showed no measurable increase in micro CT porosity when reacted with CO2 and low concentrations of SO2 over experimental timescales. This is favorable to avoid an increase in porosity and potential further leakage or migration. Fines movement however was observed in experiments and has the potential to plug permeability. The likelihood of this occurring in a reservoir could be tested in future with flow-through or core flood experiments at a range of flow rates to simulate different potential leakage scenarios. Reaction of mainly Fe-rich chlorite and minor amounts of carbonates, plagioclase and sulphides were observed via changes in the experimental water chemistry. Dissolved elements increased in concentration after CO2-SO2 injection, however several including Fe, Cr, Al, Zn, Ba, As, Pb subsequently decreased with Fe-oxide precipitation in the experiment. Concentrations of Rb and Li increased steadily in the experiment and deserve further investigation as potential indicators for monitoring a leak. Simulations indicated that siderite may mineral trap CO2 in Fe-chlorite rich rocks after 5 to 10 years, with pyrite trapping dissolved S.
Future work is suggested to react different lithologies of core overlying potential CO2 sites (e.g., mudstones, sandstones, carbonate cemented core), such as the Hutton sandstone, and to directly compare reactions with pure CO2 or CO2 containing SO2, NOx or O2 to understand more broadly the potential impacts to porosity, permeability or water chemistry in the unlikely event that a leak were to occur. The use of both batch reactors and flow-through experiments is also suggested to determine if minor gases in CO2 streams such as SO2, NOx or O2 are transported or depleted by reaction with formation water and minerals when moving through the core subsurface.
Supplementary Materials: The following are available online at www.mdpi.com/link, Figure S1: Stratigraphic column and map, Figure S2: QEMSCAN images of a section of sub-plugs, Table S1: QEMSCAN minerals (%), Table S2: XRD % semi quantitative mineral components, Figure S3: Photo of coal in Hutton Sandstone core, Figure S4: QEMSCAN selected areas of WW1 724.1 m, Figure S5: Pre-reaction tomogram and area, Figure S6: Post-reaction tomogram and movement of fines, Figure S7, and S8: Water chemistry during reaction of WW1 724.1 m, Figure S9: Photo of the rock surface before and after reaction with brown coloration. Table S3: Volume percentages of X-ray distinct components, Table S4: Geochemical model input. Figure S10: Geochemical models over 5 years, Table S5: Rock core acid digest total metal content or lithium metaborate fusion major elements and loss on ignition data.
Author Contributions: J.P. and S.G. conceived and designed the experiments; J.P. performed the experiments; J.P. and G.D. analyzed the data; S.S. analyzed micro CT and QEMSCAN before and after reaction; J.P. wrote the paper.