Oxide solution mechanisms in silicate melts: Systematic variations in the activity coefficient of SiO2

Phase equilibria and spectroscopic data are used to develop a simple model for the interaction of various oxide components and molten SiO₂. Network modifying oxides, M_xO_y produce nonbridging oxygens thereby depolymerizing the SiO₂ network. The energetics of nonbridging oxygen formation are least favorable when the field strength of the metal cation is high. This produces relatively strong M-O and Si-O-Si bridging bonds at the expense of weaker Si-O-M bonds (De Jonget al., 1980). This relationship is manifested by an increase in positive deviations from ideality with increasing cation field strength in M_xO_y-SiO₂ systems; the activity coefficient of SiO₂ is inversely correlated with Si-O-M bond strength. Network forming oxides (aluminates, phosphates, titanates. zirconates, etc.) may copolymerize with the SiO₂ network. Mixing on the same quasi-lattice produces solutions which approach ideality. Deviations from ideality in such solutions can be linked to distortions in the SiO₂ network. Discrete anion formers (phosphates, titanates, chromates, zirconates) complex with metal oxides other than SiO₂ to form discrete structural units which do not copolymerize with SiO₂. The SiO₂ network is essentially shielded from the high charge density cations in such systems and unmixing is common. As a result, the relative deviations from ideality in such melts are high. It is important to recognize that oxides such as P₂O₅, TiO₂ and ZrO₂ may act as either network-formers or discrete anion formers depending upon melt composition, and are probably distributed between these two “sites” in most geologically important liquids. The latter structural role is favored in more basic compositions.     

CO2 solubility in dacitic melts equilibrated with H2O-CO2 fluids: Implications for modeling the solubility of CO2 in silicic melts

The solubility of CO₂ in dacitic melts equilibrated with H₂O-CO₂ fluids was experimentally investigated at 1250°C and 100 to 500 MPa. CO₂ is dissolved in dacitic glasses as molecular CO₂ and carbonate. The quantification of total CO₂ in the glasses by mid-infrared (MIR) spectroscopy is difficult because the weak carbonate bands at 1430 and 1530 cm⁻¹ can not be reliably separated from background features in the spectra. Furthermore, the ratio of CO_2,mol/carbonate in the quenched glasses strongly decreases with increasing water content. Due to the difficulties in quantifying CO₂ species concentrations from the MIR spectra we have measured total CO₂ contents of dacitic glasses by secondary ion mass spectrometry (SIMS).At all pressures, the dependence of CO₂ solubility in dacitic melts on x^fluid_CO2,total shows a strong positive deviation from linearity with almost constant CO₂ solubility at x_CO2fluid > 0.8 (maximum CO₂ solubility of 795 ± 41, 1376 ± 73 and 2949 ± 166 ppm at 100, 200 and 500 MPa, respectively), indicating that dissolved water strongly enhances the solubility of CO₂. A similar nonlinear variation of CO₂ solubility with x_CO2fluid has been observed for rhyolitic melts in which carbon dioxide is incorporated exclusively as molecular CO₂ (Tamic et al., 2001). We infer that water species in the melt do not only stabilize carbonate groups as has been suggested earlier but also CO₂ molecules.A thermodynamic model describing the dependence of the CO₂ solubility in hydrous rhyolitic and dacitic melts on T, P, f_CO2 and the mol fraction of water in the melt (x_water) has been developed. An exponential variation of the equilibrium constant K₁ with x_water is proposed to account for the nonlinear dependence of x_CO2,totalmelt on x_CO2fluid. The model reproduces the CO₂ solubility data for dacitic melts within ±14% relative and the data for rhyolitic melts within 10% relative in the pressure range 100-500 MPa (except for six outliers at low x_CO2fluid). Data obtained for rhyolitic melts at 75 MPa and 850°C show a stronger deviation from the model, suggesting a change in the solubility behavior of CO₂ at low pressures (a Henrian behavior of the CO₂ solubility is observed at low pressure and low H₂O concentrations in the melt). We recommend to use our model only in the pressure range 100-500 MPa and in the x_CO2fluid range 0.1-0.95. The thermodynamic modeling indicates that the partial molar volume of total CO₂ is much lower in rhyolitic melts (31.7 cm³/mol) than in dacitic melts (46.6 cm³/mol). The dissolution enthalpy for CO₂ in hydrous rhyolitic melts was found to be negligible. This result suggests that temperature is of minor importance for CO₂ solubility in silicic melts.     

Argon and CO2 on the race track in silicate melts: A tool for the development of a CO2 speciation and diffusion model

We have analysed the kinetics of Argon and CO₂ diffusion in simplified iron free rhyolitic to hawaiitic melts using the diffusion couple technique. The concentration distance profiles of Ar and CO₂ were measured with electron microprobe analysis and Fourier Transform Infrared Spectroscopy, respectively. Error functions were fitted to the symmetrical concentration distance profiles to extract the diffusion coefficients.In the temperature range 1373 to 1773 K the activation energies for Ar diffusion range from 169 ± 20 to 257 ± 62 kJ mol⁻¹. Ar diffusivity increases exponentially with the degree of depolymerisation. In contrast, the mobility of total CO₂, that is identical to Ar mobility in rhyolitic melt, keeps constant with changing bulk composition from rhyolite to hawaiite. CO₂ speciation at 1623 K and 500 MPa was modeled for the range of compositions studied using the diffusion data of Ar and total CO₂ in combination with network former diffusion calculated from viscosity data. Within error this model is in excellent agreement with CO₂ speciation data extrapolated from temperatures near the glass transition temperature for dacitic melt composition. This model shows that even in highly depolymerised hawaiitic and tholeiitic melts molecular CO₂ is a stable species and contributes 70 to 80% to the total CO₂ diffusion, respectively.     

Carbon dioxide and argon diffusion in silicate melts: Insights into the CO2 speciation in magmas

To investigate the influence of temperature and composition on the diffusivities of dissolved carbon dioxide and argon in silicate melts, diffusion experiments were performed at magmatic pressure and temperature conditions in (a) albite melts with excess Na₂O (0-8.6 wt%) and a constant Si/Al ratio of 3, and (b) albite₇₀quartz₃₀ to jadeite melts with decreasing SiO₂ content and a constant Na/Al ratio of 1. We obtained diffusion coefficients at 500 MPa and 1323-1673 K. In the fully polymerized system Ab₇₀Qz₃₀ - Jd, the change in composition only has a weak effect on bulk CO₂ diffusivity, but Ar diffusivity increases clearly with decreasing SiO₂ content. In the system Ab + Na₂O, bulk CO₂ and Ar diffusivity increase significantly with gradual depolymerisation. The relatively small change in composition on molar basis in the depolymerized system leads to a significantly larger change in diffusivities compared to the fully polymerized Ab₇₀Qz₃₀-Jd join. Within error, activation energies for bulk CO₂ and Ar diffusion in both systems are identical with decreasing silica content (Ab + Na₂O: 159 ± 25 kJ mol⁻¹ for bulk CO₂ and 130 ± 8 kJ mol⁻¹ for Ar; Ab₇₀Qz₃₀-Jd: 163 ± 16 kJ mol⁻¹ for bulk CO₂ and 148 ± 15 kJ mol⁻¹ for Ar) even though this results in depolymerisation in one system and not the other.Although there is a variation in CO₂ speciation with changing composition as observed in quenched glasses, it has previously established that this is not a true representation of the species present in the melt, with the ratio of molecular CO₂ to carbonate decreasing during quenching. Thus, diffusion coefficients for the individual CO₂ species cannot be directly derived by measuring molecular CO₂ and CO₃^2- concentration-distance profiles in the glasses. To obtain diffusivities of individual CO₂ species, we have made two assumptions that (1) inert Ar can be used as a proxy for molecular CO₂ diffusion characteristics as shown by our previous work and (2) the diffusivity of CO₃²⁻ can be calculated assuming it is identical to network forming components (Si⁴⁺ and Al³⁺). This is derived from viscosity data (Eyring eqn.) and suggests that CO₃²⁻ diffusion would be several orders of magnitude slower than molecular CO₂ diffusion.The systematics of measured bulk CO₂ diffusivity rates and comparison with the Ar proxy all suggest that the faster molecular CO₂ species is much more dominant in melts than measurements on resulting quenched glasses would suggest. This study has confirmed an observation of surprisingly consistent bulk CO₂ diffusivity across a range of natural compositions were Ar diffusivity significantly increases. This is consistent with an actual increase in molecular CO₂ mobility (similar to Ar) that is combined with an increase in the proportion of the slower carbonate in the melt.These results demonstrate that the CO₂ diffusion and speciation model provides an insight into the transport processes in the melt and is promising and an alternative tool to in situ speciation measurements at magmatic conditions, which at the moment are technically extremely difficult. We present the first high pressure high temperature in situ MIR spectra of a CO₂ bearing albitic glass/melt suggesting that molecular CO₂ is a stable species at high temperature, which is qualitatively consistent with the modelled CO₂ speciation data.     

The role of P2O5 in silicate melts

Phase equilibria data in the systems SiO₂-P₂O₅, P₂O₅-M_xO_y, and P₂O₅-M_xO_y-SiO₂ are employed in conjunction with Chromatographic and spectral data to investigate the role of P₂O₅ in silicate melts. Such data indicate that the behavior of P₂O₅ is complex. P₂O₅ depolymerizes pure SiO₂ melts by entering the network as a four-fold coordinated cation, but polymerizes melts in which an additional metal cation other than silicon is present. The effect of this polymerization is apparent in the widening of the granite-ferrobasalt two-liquid solvus. In this complex system P₂O₅ acts to increase phase separation by further enrichment of the high charge density cations Ti, Fe, Mg, Mn, Ca, in the ferrobasaltic liquid. P₂O₅ also produces an increase in the ferrobasalt-granite REE liquid distribution coefficients. These distribution coefficients are close to 4 in P₂O₅-free melts, but close to 15 in P₂O₅-bearing melts.The dual behavior of P₂O₅ is explained in a model which requires complexing of phosphate anions (analogous to silicate anions) and metal cations in the melt. This interaction destroys Si-O-M-O-Si bonds polymerizing the melt. The higher concentration of Si-O-M-O-Si bond complexes in immiscible ferrobasaltic liquids relative to their conjugate immiscible granite liquids explains the partitioning of P₂O₅ into the ferrobasaltic liquid.     

Diffusion in silicate melts: II. Multicomponent diffusion in CaOAl2O3SiO2 at 1500°C and 1 GPa

Chemical diffusion profiles in molten CaO---Al₂O₃---SiO₂ have been measured over a large range of compositions at 1500°C and l GPa. The diffusion profiles have been inverted for effective binary diffusion coefficients (EBDCs) and for the chemical diffusion matrix. The EBDCs are shown to depend strongly on both composition and direction of diffusion in composition space. The dependence of EBDCs on direction in composition space, which for the system studied here can be as large as a factor of seven, severely limits the applicability of EBDCs to interdiffusion in any direction other than the one used to derive the EBDCs.

The chemical diffusion matrix for molten CaO---Al₂O₃---SiO₂ was determined using diffusion profiles from two or three mutually orthogonal diffusion couples in the ternary composition space. All features of the diffusion profiles shown in this work can be reproduced by representing the chemical fluxes in the three-component system as a linear combination of concentration gradients via a 2 × 2 diffusion matrix. Chemical diffusion in molten CaO---Al₂O₃---SiO₂ shows clear evidence of strong diffusive coupling among the components. This can be seen in the uphill diffusion profiles of components that were initially uniform, in the fact that the apparent rate of diffusion of some components is a strong function of direction in composition space, and most quantitatively in the magnitude of off-diagonal elements of the diffusion matrix relative to the magnitude of the diagonal elements. SiO₂ for example, is found to be strongly coupled with CaO in relatively silicic melts, whereas Al₂O₃ is strongly coupled with CaO in less silicic melts. Furthermore, the coupling of CaO with either Al₂O₃ or Si02 reverses sign between more and less polymerized compositions. Interdiffusion profiles in natural melts have numerous features that suggest similar coupling between Al₂O₃ and CaO and between SiO₂ and CaO.     

Oxygen isotope exchange between H2O and CO2 at elevated CO2 pressures: Implications for monitoring of geological CO2 storage

Traditionally, the application of stable isotopes in Carbon Capture and Storage (CCS) projects has focused on δ¹³C values of CO₂ to trace the migration of injected CO₂ in the subsurface. More recently the use of δ¹⁸O values of both CO₂ and reservoir fluids has been proposed as a method for quantifying in situ CO₂ reservoir saturations due to O isotope exchange between CO₂ and H₂O and subsequent changes in δ¹⁸O_H2O values in the presence of high concentrations of CO₂. To verify that O isotope exchange between CO₂ and H₂O reaches equilibrium within days, and that δ¹⁸O_H2O values indeed change predictably due to the presence of CO₂, a laboratory study was conducted during which the isotope composition of H₂O, CO₂, and dissolved inorganic C (DIC) was determined at representative reservoir conditions (50 °C and up to 19 MPa) and varying CO₂ pressures. Conditions typical for the Pembina Cardium CO₂ Monitoring Pilot in Alberta (Canada) were chosen for the experiments. Results obtained showed that δ¹⁸O values of CO₂ were on average 36.4 ± 2.2‰ (1σ, n = 15) higher than those of water at all pressures up to and including reservoir pressure (19 MPa), in excellent agreement with the theoretically predicted isotope enrichment factor of 35.5‰ for the experimental temperatures of 50 °C. By using ¹⁸O enriched water for the experiments it was demonstrated that changes in the δ¹⁸O values of water were predictably related to the fraction of O in the system sourced from CO₂ in excellent agreement with theoretical predictions. Since the fraction of O sourced from CO₂ is related to the total volumetric saturation of CO₂ and water as a fraction of the total volume of the system, it is concluded that changes in δ¹⁸O values of reservoir fluids can be used to calculate reservoir saturations of CO₂ in CCS settings given that the δ¹⁸O values of CO₂ and water are sufficiently distinct.     

CO2 solubility and speciation in intermediate (andesitic) melts: the role of H2O and composition

We determined total CO₂ solubilities in andesite melts with a range of compositions. Melts were equilibrated with excess C-O(-H) fluid at 1 GPa and 1300°C then quenched to glasses. Samples were analyzed using an electron microprobe for major elements, ion microprobe for C-O-H volatiles, and Fourier transform infrared spectroscopy for molecular H₂O, OH⁻, molecular CO₂, and CO₃²⁻. CO₂ solubility was determined in hydrous andesite glasses and we found that H₂O content has a strong influence on C-O speciation and total CO₂ solubility. In anhydrous andesite melts with ∼60 wt.% SiO₂, total CO₂ solubility is ∼0.3 wt.% at 1300°C and 1 GPa and total CO₂ solubility increases by about 0.06 wt.% per wt.% of total H₂O. As total H₂O increases from ∼0 to ∼3.4 wt.%, molecular CO₂ decreases (from 0.07 ± 0.01 wt.% to ∼0.01 wt.%) and CO₃²⁻ increases (from 0.24 ± 0.04 wt.% to 0.57 ± 0.09 wt.%). Molecular CO₂ increases as the calculated mole fraction of CO₂ in the fluid increases, showing Henrian behavior. In contrast, CO₃²⁻ decreases as the calculated mole fraction of CO₂ in the fluid increases, indicating that CO₃²⁻ solubility is strongly dependent on the availability of reactive oxygens in the melt. These findings have implications for CO₂ degassing. If substantial H₂O is present, total CO₂ solubility is higher and CO₂ will degas at relatively shallow levels compared to a drier melt. Total CO₂ solubility was also examined in andesitic glasses with additional Ca, K, or Mg and low H₂O contents (<1 wt.%). We found that total CO₂ solubility is negatively correlated with (Si + Al) cation mole fraction and positively correlated with cations with large Gibbs free energy of decarbonation or high charge-to-radius ratios (e.g., Ca). Combining our andesite data with data from the literature, we find that molecular CO₂ is more abundant in highly polymerized melts with high ionic porosities (>∼48.3%), and low nonbridging oxygen/tetrahedral oxygen (<∼0.3). Carbonate dominates most silicate melts and is most abundant in depolymerized melts with low ionic porosities, high nonbridging oxygen/tetrahedral oxygen (>∼0.3), and abundant cations with large Gibbs free energy of decarbonation or high charge-to-radius ratio. In natural silicate melt, the oxygens in the carbonate are likely associated with tetrahedral and network-modifying cations (including Ca, H, or H-bonds) or a combinations of those cations.     

Interaction between sulphide and H2O in silicate melts

Reaction between dissolved water and sulphide was experimentally investigated in soda-lime-silicate (NCS) and sodium trisilicate (NS3) melts at temperatures from 1000 to 1200 °C and pressures of 100 or 200 MPa in internally heated gas pressure vessels. Diffusion couple experiments were conducted at water-undersaturated conditions with one half of the couple being doped with sulphide (added as FeS or Na₂S; 1500-2000 ppm S by weight) and the other with H₂O (∼3.0 wt.%). Additionally, two experiments were performed using a dry NCS glass cylinder and a free H₂O fluid. Here, the melt was water-saturated at least at the melt/fluid interface. Profiling by electron microprobe (sulphur) and infrared microscopy (H₂O) demonstrate that H₂O diffusion in the melts is faster by 1.5-2.3 orders of magnitude than sulphur diffusion and, hence, H₂O can be considered as a rapidly diffusing oxidant while sulphur is quasi immobile in these experiments.In Raman spectra a band at 2576 cm⁻¹ appears in the sulphide - H₂O transition zone which is attributed to fundamental S-H stretching vibrations. Formation of new IR absorption bands at 5025 cm⁻¹ (on expense of the combination band of molecular H₂O at 5225 cm⁻¹) and at 3400 cm⁻¹ was observed at the front of the in-diffusing water in the sulphide bearing melt. The appearance and intensity of these two IR bands is correlated with systematic changes in S K-edge XANES spectra. A pre-edge excitation at 2466.5 eV grows with increasing H₂O concentration while the sulphide peak at 2474.0 eV decreases in intensity relative to the peak at 2477.0 eV and the feature at 2472.3 eV becomes more pronounced (all energies are relative to the sulphate excitation, calibrated to 2482.5 eV). The observations by Raman, IR and XANES spectroscopy indicate a well coordinated S²⁻ - H₂O complex which was probably formed in the glasses during cooling at the glass transition. No oxidation of sulphide was observed in any of the diffusion couple experiments. On the contrary, XANES spectra from experiments conducted with a free H₂O fluid show complete transformation of sulphide to sulphate near the melt surface and coexistence of sulphate and sulphide in the center of the melt. This can be explained by a lower H₂O activity in the diffusion couple experiments or by the need of a sink for hydrogen (e.g., a fluid which can dissolve high concentration of hydrogen) to promote oxidation of sulphide by H₂O via the reaction S²⁻ + 4H₂O = SO₄²⁻ + 4H₂. Sulphite could not be detected in any of the XANES spectra implying that this species, if it exists in the melt, it is a subordinate or transient species only.     

Solubility and solution mechanism of H2O in alkali silicate melts and glasses at high pressure and temperature

The solubility behavior of H₂O in melts in the system Na₂O-SiO₂-H₂O was determined by locating the univariant phase boundary, melt = melt + vapor in the 0.8-2 GPa and 1000°-1300°C pressure and temperature range, respectively. The NBO/Si-range of the melts (0.25-1) was chosen to cover that of most natural magmatic liquids. The H₂O solubility in melts in the system Na₂O-SiO₂-H₂O (X_H2O) ranges between 18 and 45 mol% (O = 1) with (∂X_H2O/∂P)_T∼14-18 mol% H₂O/GPa. The (∂X_H2O/∂P)_T is negatively correlated with NBO/Si (= Na/Si) of the melt. The (∂X_H2O/∂T)_P is in the −0.03 to +0.05 mol% H₂O/°C range, and is negatively correlated with NBO/Si. The [∂X_H2O/∂(NBO/Si)]_P,T is in the −3 to −8 mol% H₂O/(NBO/Si) range. Melts with NBO/Si similar to basaltic liquids (∼0.6-∼1.0) show (∂X_H2O/∂T)_P<0, whereas more polymerized melts exhibit (∂X_H2O/∂T)_P>0. Complete miscibility between hydrous melt and aqueous fluid occurs in the 0.8-2 GPa pressure range for melts with NBO/Si ≤0.5 at T >1100°C. Miscibility occurs at lower pressure the more polymerized the melt.     

Na2O solubility in CaO-MgO-SiO2 melts

The sodium solubility in silicate melts in the CaO-MgO-SiO₂ (CMS) system at 1400 °C has been measured by using a closed thermochemical reactor designed to control alkali metal activity. In this reactor, Na_(g) evaporation from a Na₂O-xSiO₂ melt imposes an alkali metal vapor pressure in equilibrium with the molten silicate samples. Because of equilibrium conditions in the reactor, the activity of sodium-metal oxide in the molten samples is the same as that of the source, i.e., aNa₂O_(sample) = aNa₂O_(source). This design also allows to determine the sodium oxide activity coefficient in the samples. Thirty-three different CMS compositions were studied. The results show that the amount of sodium entering from the gas phase (i.e., Na₂O solubility) is strongly sensitive to silica content of the melt and, to a lesser extent, the relative amounts of CaO and MgO. Despite the large range of tested melt compositions (0 < CaO and MgO < 40; 40 < SiO₂ < 100; in wt%), we found that Na₂O solubility is conveniently modeled as a linear function of the optical basicity (Λ) calculated on a Na-free basis melt composition. In our experiments, γNa₂O_(sample) ranges from 7 × 10⁻⁷ to 5 × 10⁻⁶, indicating a strongly non-ideal behavior of Na₂O solubility in the studied CMS melts (γNa₂O_(sample) ? 1). In addition to showing the effect of sodium on phase relationships in the CMS system, this Na₂O solubility study brings valuable new constraints on how melt structure controls the solubility of Na in the CMS silicate melts. Our results suggest that Na₂O addition causes depolymerization of the melt by preferential breaking of Si-O-Si bonds of the most polymerized tetrahedral sites, mainly Q⁴.     

硅酸盐和金属氧化物矿物氧同位素组成的CO2激光氟化分析

我室采用MIR-10型CO2激光器,在一种富BrF5的氛围中使激光对硅酸盐和氧化物矿物样品加热形成O2,经多次纯化后用5A的分子筛吸收,再直接送至气体质谱仪进行氧同位素比值测定.这个实验流程与传统方法相比的改进不仅在使用激光加热技术及样品的放置上,而且在直接采用O2而不是CO2进行质谱测定.采用O2进行直接分析的优点不仅避免了向CO2转化过程中的潜在同位素分馏,而且能够得到样品的δ17O值,因此为宇宙样品分析提供了可能.CO2激光氟化技术的优点是所需样品量小(可低达1～2 mg),因此能够分析微小岩石区域或单矿物晶体内的氧同位素分布.同时,激光可以达到非常高的温度(＞4000K),因此能够对某些难熔矿物(如锆石、蓝晶石、橄榄石等)进行氧同位素分析.     

PbO-SiO2 melts: structure and thermodynamics of mixing

Thermodynamic properties of PbO-SiO₂ melts, obtained from published data and calculated from freezing point depressions, reflect the gradual polymerization of silicate anions in the melt as the $\text{SiO}{}_2\text{PbO}$ ratio is increased. The free energy of mixing curve at 1000°C has a minimum at 40 mole % SiO₂ and is convex-upward between 72 and 98 mole % SiO₂. The latter is an indication of metastable liquid immiscibility. The free energy minimum is correlated with the maximum in the distribution of nonbridging oxygens in the melt. In SiO₂-poor melts, the activities of PbO and SiO₂ (pure liquid standard states) show sharp negative deviations from ideality. The PbO activity reflects the paucity of free oxygen species in the melt whereas the SiO₂ activity reflects the depolymerized state of the silicate anions. In more SiO₂-rich melts, the activity of SiO₂ shows a positive deviation from ideality which is qualitatively correlated to a polymerization parameter. The heat of mixing term has a minimum of ?2000 cal at 35 mole % SiO₂ and a maximum of +200 cal at 90 mole % SiO₂. The minimum is associated with the exothermic heat effect obtained during the reaction (O⁰) + (O^2?) = 2(O^?), whereas the maximum corresponds to the endothermic heat effect obtained when coordination polyhedra of oxygens form around the Pb cation. The entropy of mixing curve has the same form but is systematically smaller than a theoretical curve calculated on the assumption of random mixing of oxygen species. The discrepancy is due to the entropy loss obtained by the clustering of oxygen species to form complex silicate species.     

Caprock interaction with CO2: A laboratory study of reactivity of shale with supercritical CO2 and brine

Crushed rock from two caprock samples, a carbonate-rich shale and a clay-rich shale, were reacted with a mixture of brine and supercritical CO₂ (CO₂–brine) in a laboratory batch reactor, at different temperature and pressure conditions. The samples were cored from a proposed underground CO₂ storage site near the town of Longyearbyen in Svalbard. The reacting fluid was a mixture of 1 M NaCl solution and CO₂ (110 bar) and the water/rock ratio was 20:1. Carbon dioxide was injected into the reactors after the solution had been bubbled with N₂, in order to mimic O₂-depleted natural storage conditions. A control reaction was also run on the clay-rich shale sample, where the crushed rock was reacted with brine (CO₂-free brine) at the same experimental conditions. A total of 8 batch reaction experiments were run at temperatures ranging from 80 to 250 °C and total pressures of 110 bar (∼40 bar for the control experiment). The experiments lasted 1–5 weeks.Fluid analysis showed that the aqueous concentration of major elements (i.e. Ca, Mg, Fe, K, Al) and SiO₂ increased in all experiments. Release rates of Fe and SiO₂ were more pronounced in solutions reacted with CO₂–brine as compared to those reacted with CO₂-free brine. For samples reacted with the CO₂–brine, lower temperature reactions (80 °C) released much more Fe and SiO₂ than higher temperature reactions (150–250 °C). Analysis by SEM and XRD of reacted solids also revealed changes in mineralogical compositions. The carbonate-rich shale was more reactive at 250 °C, as revealed by the dissolution of plagioclase and clay minerals (illite and chlorite), dissolution and re-precipitation of carbonates, and the formation of smectite. Carbon dioxide was also permanently sequestered as calcite in the same sample. The clay-rich shale reacted with CO₂–brine did not show major mineralogical alteration. However, a significant amount of analcime was formed in the clay-rich shale reacted with CO₂-free brine; while no trace of analcime was observed in either of the samples reacted with CO₂–brine.     

GEOCARBSULF: A combined model for Phanerozoic atmospheric O2 and CO2

A model for the combined long-term cycles of carbon and sulfur has been constructed which combines all the factors modifying weathering and degassing of the GEOCARB III model [Berner R.A., Kothavala Z., 2001. GEOCARB III: a revised model of atmospheric CO₂ over Phanerozoic time. Am. J. Sci. 301, 182-204] for CO₂ with rapid recycling and oxygen dependent carbon and sulfur isotope fractionation of an isotope mass balance model for O₂ [Berner R.A., 2001. Modeling atmospheric O₂ over Phanerozoic time. Geochim. Cosmochim. Acta65, 685-694]. New isotopic data for both carbon and sulfur are used and new feedbacks are created by combining the models. Sensitivity analysis is done by determining (1) the effect on weathering rates of using rapid recycling (rapid recycling treats carbon and sulfur weathering in terms of young rapidly weathering rocks and older more slowly weathering rocks); (2) the effect on O₂ of using different initial starting conditions; (3) the effect on O₂ of using different data for carbon isotope fractionation during photosynthesis and alternative values of oceanic δ¹³C for the past 200 million years; (4) the effect on sulfur isotope fractionation and on O₂ of varying the size of O₂ feedback during sedimentary pyrite formation; (5) the effect on O₂ of varying the dependence of organic matter and pyrite weathering on tectonic uplift plus erosion, and the degree of exposure of coastal lands by sea level change; (6) the effect on CO₂ of adding the variability of volcanic rock weathering over time [Berner, R.A., 2006. Inclusion of the weathering of volcanic rocks in the GEOCARBSULF model. Am. J. Sci.306 (in press)]. Results show a similar trend of atmospheric CO₂ over the Phanerozoic to the results of GEOCARB III, but with some differences during the early Paleozoic and, for variable volcanic rock weathering, lower CO₂ values during the Mesozoic. Atmospheric oxygen shows a major broad late Paleozoic peak with a maximum value of about 30% O₂ in the Permian, a secondary less-broad peak centered near the Silurian/Devonian boundary, variation between 15% and 20% O₂ during the Cambrian and Ordovician, a very sharp drop from 30% to 15% O₂ at the Permo-Triassic boundary, and a more-or less continuous rise in O₂ from the late Triassic to the present.     

Melting relationships in the systems CaOCO2 and MgOCO2 to 33 kilobars

The melting temperatures of calcite and magnesite in the presence of excess CO₂ have been measured using Ag₂C₂O₄ in sealed capsules m a piston-cylinder apparatus. At 27 kbar, 11.5 wt % CO₂ dissolves in molten CaCO₂, depressing the freezing temperature from 1610 to 1505°C; and 6.5 wt % CO₂ dissolves in molten MgCO₃, depressing the freezing temperature from 1590 to 1510°C. The eutectic between calcite and lime was located at 1385°C at 27 kbar. These and other new results, combined with previously published data, permit completion of PT diagrams for the systems CaO-CO₂ and MgO-CO₂ from 1 bar to 35 kbar. The dissociation curve for each carbonate terminates at an invariant point where melting begins, at 40 bars and 1230°C for CaO-CO₂ and 23 kbar and 1550°C for MgO-CO₂ The differences between the two systems are explained by the different solubilities of CO₂ in the invariant liquids consequent upon the large pressure difference between the locations of these two invariant points. The results show that the temperatures for the beginning of melting of carbonates in the asthenosphere are lowered by about 100°C in the presence of CO₂.     

Selection of coals of different maturities for CO2 Storage by modelling of CH4 and CO2 adsorption isotherms

CO₂ injection in unmineable coal seams could be one interesting option for both storage and methane recovery processes. The objective of this study is to compare and model pure gas sorption isotherms (CO₂ and CH₄) for well-characterised coals of different maturities to determine the most suitable coal for CO₂ storage. Carbon dioxide and methane adsorption on several coals have been investigated using a gravimetric adsorption method. The experiments were carried out using both CO₂ and CH₄ pure gases at 25 °C from 0.1 to 5 MPa (1 to 50 bar). The experimental results were fitted using Temkin's approach but also with the corrected Langmuir's and the corrected Tóth's equations. The two last approaches are more accurate from a thermodynamical point of view, and have the advantage of taking into account the fact that experimental data (isotherms) correspond to excess adsorption capacities. These approaches allow better quantification of the adsorbed gas. Determined CO₂ adsorption capacities are from 0.5 to 2 mmol/g of dry coal. Modelling provides also the affinity parameters of the two gases for the different coals. We have shown these parameters determined with adsorption models could be used for classification and first selection of coals for CO₂ storage. The affinity ratio ranges from a value close to 1 for immature coals to 41 for high rank coals like anthracites. This ratio allows selecting coals having high CO₂ adsorption capacities. In our case, the modelling study of a significant number of coals from various ranks shows that anthracites seem to have the highest CO₂ storage capacities. Our study provides high quality affinity parameters and values of CO₂ and CH₄ adsorption capacities on various coals for the future modelling of CO₂ injection in coal seams.     

深部咸含水层CO2注入下压力抬升与流体迁移的大尺度模拟研究:以松辽盆地为例

地下深部封存CO2已经被公认是人类削减温室气体排放的一条有效而又科学的途径。深部咸含水层CO2地质封存因封存潜力巨大,技术可行,且已有实际的工程运行,因而备受关注。松辽盆地是中国潜在的CO2储存场地之一,选择松辽盆地为大尺度模拟研究对象,选取姚家组砂岩层为储层,选取嫩江组泥岩为盖层,运用TOUGH-MP并行计算代码建立了覆盖整个松辽盆地的三维地质模型,在中央凹陷区开展大尺度CO2注入模拟研究,包括CO2运移、储存、地层压力提升以及储存安全性等问题。模拟结果表明:持续注入100a后形成的CO2羽远小于产生的压力积聚区影响范围。注入产生的压力抬升将在注入停止后迅速消散,不会对区域地层压力和浅层地下水系统产生显著影响。在千年之内注入的CO2将随着时间持续,逐渐溶解于水中,而不会因盖层微弱的渗透性而逃逸。     

火山岩吸附CO2气的成藏潜力及实例分析

火山岩的脱气实验和对昌德东CO2气藏气源的分析结果表明加热火山岩到250℃时,脱出挥发分总量为0.0299～0.0790mL/g,其中CO2脱出量为0.0218～0.0706mL/g(0.429～1.387wt%);挥发组分以CO2为主,还含有H2、CO、CH4等还原性气体,以及少量低碳烷烃,CO2含量和总烃呈现反比关系;基性岩的CO2脱出量、脱出率高于中、酸性岩;CO2脱出量与岩石碱质含量正相关.松辽盆地北部昌德东CO2气藏成藏模式为"自生自储",成藏CO2气主要来自深部被火山岩吸附的气.随岩浆上升,在岩浆冷凝成火山岩的过程中被吸附于火山岩的节理、劈理和晶体位错之中的CO2气,连同火山岩包体中的残留气,成为高纯CO2气藏的主要补给源,并非地幔气体沿大断裂上来直接充注成藏.     

Mitigation of asphaltics deposition during CO2 flood by enhancing CO2 solvency with chemical modifiers

CO₂ injected in the reservoir of McElroy Field, TX, for a CO₂ flood was in the supercritical state. Supercritical CO₂ fluid is capable of extracting light and intermediate hydrocarbons from rocks but is unable to extract heavy hydrocarbons and asphaltics. Therefore, plugging of asphaltics in reservoir rocks and a consequent reduction in injectivity and recovery may result when CO₂ only is used in enhanced oil recovery. By adding common solvents as chemical modifiers, the flooding fluid shows marked improvement in solvency for heavy components of crudes due to its increased density and polarity. Numerous supercritical CO₂ fluid extractions of dolomite rock from the Grayburg Formation containing known amount of spiked McElroy crude oil have been carried out to evaluate extraction efficiencies of CO₂ and CO₂ with chemical modifiers at various temperatures and pressures. All experiments show that extraction efficiency increases with increasing CO₂ pressure but decreases with increasing temperature. Addition of chemical modifiers to CO₂ also shows improved extraction efficiency and reduced asphaltic deposits. Under the pressure and temperature similar to McElroy reservoir conditions; chemically modified CO₂ yielded almost 3 times as much oil extracts as those in extractions with CO₂ only. It also reduced the asphaltics content in extracted rocks to nearly one half; indicating its potential for mitigating asphaltics plugging of formation rocks