N2-CH4（CO2）混合气体在线标样制备及其拉曼定量因子测定

利用混合气体的标准样品对激光拉曼探针进行标定,可以快速准确地对包裹体中的无机及有机气相组分进行定量分析。而常用的商用钢瓶装混合气体标样,存在费用高、气体组成单一固定等缺点。本文设计了一套在线标样制备装置,提出一种在线配置不同浓度和压力条件下混合气体标样的方法。利用高纯度(纯度99.999%)的N2、CH4以及CO2钢瓶气,经过在线混合增压,在5 MPa和10 MPa条件下制备了N2摩尔分数为30%、50%和70%的N2-CH4以及N2-CO2混合气体在线标样。该方法制备的标样与70%N2+30%CO2的商用钢瓶气标样对比表明,CO2与N2的拉曼相对峰高以及相对峰面积值的误差在4%以内,具有较高的准确度和重现性。通过不同压力和浓度条件下CH4以及CO2的拉曼相对定量因子测定表明,气体的相对定量因子在5~10 MPa压力条件下与压力及组成无关。地质样品应用结果表明,本方法可以方便、灵活、准确地按任意比例将两瓶及两瓶以上纯气体钢瓶样品进行混合及增压,为激光拉曼标定、气体组成原位测量等提供了一种新的技术思路。     

A high temperature equation of state for the H2O-CaCl2 and H2O-MgCl2 systems

An equation of state (EOS) is developed for salt-water systems in the high temperature range. As an example of the applications, this EOS is parameterized for the calculation of density, immiscibility, and the compositions of coexisting phases in the CaCl₂-H₂O and MgCl₂-H₂O systems from 523 to 973 K and from saturation pressure to 1500 bar. All available volumetric and phase equilibrium measurements of these binaries are well represented by this equation. This EOS is based on a Helmholtz free energy representation constructed from a reference system containing hard-sphere and polar contributions plus an empirical correction. For the temperature and pressure range in this study, the electrolyte solutes are assumed to be associated. The water molecules are modeled as hard spheres with point dipoles and the solute molecules, MgCl₂ and CaCl₂, as hard spheres with point quadrupoles. The free energy of the reference system is calculated from an analytical representation of the Helmholtz free energy of the hard-sphere contributions and perturbative estimates of the electrostatic contributions. The empirical correction used to account for deviations of the reference system predictions from measured data is based on a virial expansion. The formalism allows generalization to aqueous systems containing insoluble gases (CO₂, CH₄), alkali chlorides (NaCl, KCl), and alkaline earth chlorides (CaCl₂, MgCl₂). The program of this model is available as an electronic annex (see EA1 and EA2) and can also be downloaded at: http://www.geochem-model.org/programs.htm.     

The stability of UO2OH+ and UO2[HPO4]2−2 complexes at 25°C

The cumulative association constant (β₂) for the geochemically important aqueous complex UO₂[HPO₄]^2?₂ has been determined by potentiometric titration in Na₂HPO₄-UO₂(NO₃)₂ solutions in the pH range 3.9–4.7, at ionic strengths below 0.024 molal with the Newton-Raphson method used to compute β₂ from the chemical analytical data. Based on 25 measurements we obtain logβ₂ = 18.3 ± 0.2 at 25°C. From the same experiments we compute that the association constant of UO₂OH⁺ is 8.9 ± 0.1, in disagreement with the value of 8.3 ± 0.3 for this constant given by Baes and Mesmer (1976).     

低硅埃达克岩和高硅埃达克岩问题

低硅埃达克岩(low-SiO2adakites,LSA,SiO2<60%)和高硅埃达克岩(high-SiO2adakites,HSA,SiO2>60%)的术语是Martin等(2005)提出来的。本文不赞同上述分类,因为,他们的低硅埃达克岩是赞岐岩而不是埃达克岩。赞岐岩是幔源的,埃达克岩是壳源的,虽然埃达克岩可以与地幔混合形成高镁的埃达克岩,但仍然是壳源的。不应当把壳源的岩浆和幔源的岩浆混淆起来。     

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.     

The bearing of the system CaMgSi2O6CaAl2SiO6CaFeAlSiO6 on fassaitic pyroxene

The system CaMgSi₂O₆CaAl₂SiO₆CaFeAlSiO₆ has been studied in air at 1 atm. The phase assemblage at subsolidus temperatures in the CaMgSi₂O₆-rich portion is Cpx + An + Mel and that in the CaMgSi₂O₆-poor portion Cpx + An + Mel + Sp. At subsolidus temperatures the sigle-phase field of clinopyroxene increases with an increase in the CaFeAlSiO₆ component of the system. The Al₂O₃ content of clinopyroxene, however, continues to increase beyond the single-phase field and attains at least 16.04 wt.% Al₂O₃ with 3.9 wt.% Fe₂O₃. The stability field of fassaite in the system over a range of pressures and oxygen fugacities has been estimated from data in the literature as well as the present data. The CaFeAlSiO₆ content of fassaite is dependent on oxygen fugacity, but is not influenced by pressure. The stability field is strongly influenced by oxygen fugacity at low and high pressure, and decreases with decreasing oxygen fugacity. Clinopyroxenes in both volcanic and metamorphic rocks from various localities, when plotted on the CaMgSi₂O₆CaAl₂SiO₆CaFeAlSiO₆ triangle, show that there is no compositional gap between diopside and fassaitic pyroxene in metamorphic rocks, and that the fassaitic pyroxene in alkalic rocks becomes richer in both CaAl₂SiO₆ and CaFeAlSiO₅ components as crystallization proceeds. These results agree with those obtained in the experimental study.     

Actinide crystal-liquid partitioning for clinopyroxene and Ca3(PO4)2

Using fission and alpha track radiography techniques, we have measured partition coefficients (D) for the actinide elements Th, U and Pu between diopsidic clinopyroxene, whitlockite [β-Ca₃ (PO₄)₂] and silicate liquid at 20kbar. Equilibrium partitioning at the crystal-liquid interface is assumed, and corrections for actinide zoning have been applied to the measured D values. Reproducibility for both actinide and minor element D values is carefully examined as a criterion for crystal-liquid interface equilibrium. The data are mostly compatible with interface equilibrium except for experiments at high cooling rates ( ? 30 deg/hr). Partition coefficients for Th/U/Pu of about 0.002/0.002/0.06 are measured for clinopyroxene and 1.2/0.5/3.4 for whitlockite. At an oxygen fugacity of 10^?8.5, Pu is much more readily incorporated into the crystalline phases than is U or Th because of the importance of trivalent Pu. The D_Pu(cpx) is similar to D(cpx) of the light rare earths supporting the concept of Pu/(rare earth) dating.     

K-H3OFe3(SO4)2(OH)6铁矾的XRD研究

根据X射线衍射(XRD)分析发现: A Fe₃(SO₄)₂(OH)₆(A=K⁺、H₃O⁺)系列铁钒的XRD数据十分相近,难以用XRD区别,需通过能谱(EDS)辅助分析,才能区分此类铁矾。另外,此类铁矾的003和107面网间距d随K⁺含量增大而增大,且呈一元三次方程的关系;而033和220面网间距d随K⁺含量增大而减小,呈一元二次方程的关系。对该现象从铁矾晶体结构方面进行解释:K⁺、H₃O⁺离子位于较大空隙中,且沿着Z轴方向排列,当K⁺、H₃O⁺离子之间相互替换时,会导致该铁矾晶体结构在Z轴方向有较明显的变化。     

Laboratory study of self-potential signals during releasing of CO2 and N2 plumes

Carbon Capture Sequestration (CCS) projects require, for safety reasons, monitoring programmes focused on surveying gas leakage on the surface. Generally, these programmes include detection of chemical tracers that, once on the surface, could be associated with CO₂ degassing. We take a different approach by analysing feasibility of applying electrical surface techniques, specifically Self-Potential. A laboratory-scale model, using water-sand, was built for simulating a leakage scenario being monitored with non-polarisable electrodes. Electrical potentials were measured before, during and after gas injection (CO₂ and N₂) to determine if gas leakage is detectable. Variations of settings were done for assessing how the electrical potentials changed according to size of electrodes, distance from electrodes to the gas source, and type of gas. Results indicated that a degassing event is indeed detectable on electrodes located above injection source. Although the amount of gas could not be quantified from signals, injection timespan and increasing of injection rate were identified. Even though conditions of experiments were highly controlled contrasting to those usually found at field scale, we project that Self-Potential is a promising tool for detecting CO₂ leakage if electrodes are properly placed.     

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.     

Thermochemistry of glasses and liquids in the systems CaMgSi2O6-CaAl2Si2O8-NaAlSi3O8, SiO2-CaAl2Si2O8-NaAlSi3O8 and SiO2-Al2O3-CaO-Na2O

Enthalpies of solution in 2PbO· B₂O₃ at 712°C have been measured for glasses in the systems albite anorthite diopside, NaAlO₂-SiO₂, Ca_0.5AlO₂-SiO₂ and albite-anorthite-quartz. The systems albite-anorthite and diopside-anorthite show substantial negative enthalpies of mixing, albite-diopside shows significant positive heats of mixing. For compositions up to NaAlO₂ = 0.42 (which includes the subsystem albite-silica) the system NaAlO₂-SiO₂ shows essentially zero heats of mixing. A negative ternary excess heat of mixing is found in the plagioclase-rich portion of the albite-anorthite-diopside system. The join Si₄O₈-CaAl₂Si₂O₈ shows small but significant heats of mixing. In albite-anorthite-quartz. ternary glasses, the ternary excess enthalpy of mixing is positive.Based on available heat capacity data and appropriate consideration of the glass transition, the enthalpy of the crystal-glass transition (vitrification) is a serious underestimate of the enthalpy of the crystal-liquid transition (fusion) especially when the melting point, T_f, is many hundreds of degrees higher than the glass transition temperature, T_g. On the other hand, the same heat capacity data suggest that the enthalpies of mixing in albite-anorthite-diopside liquids are calculated to be quite similar to those in the glasses. The enthalpies of mixing observed in general support the structural models proposed by Taylor and Brown (1979a, b) and others for the structure of aluminosilicate glasses.     

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.     

The rise of trees and their effects on Paleozoic atmospheric CO2 and O2

The rise of large vascular plants during the mid-Paleozoic brought about a major increase in the rates of weathering of silicate minerals that induced a drop in the level of atmospheric CO₂ and contributed, via the atmospheric greenhouse effect, to global cooling and the initiation of the most long lived and a really extensive glaciation of the past 550 million years. Sedimentary burial of the microbiologically resistant remains of the plants resulted during the Permo-Carboniferous in both further lowering of CO₂ and in elevation of atmospheric O₂. Evidence of changes in CO₂ and O₂ are provided by mathematical models, studies of paleosols, fossil plants, fossil insects, and the effects of modern plants on silicate weathering, and by laboratory studies of the effects of changes in O₂ on plants and insects. To cite this article: R.A. Berner, C. R. Geoscience 335 (2003).     

Groundwater N speciation and redox control of organic N mineralization by O2 reduction to H2O2

Excess N from agriculture induces eutrophication in major river systems and hypoxia in coastal waters throughout the world. Much of this N is from headwaters far up the watersheds. In turn, much of the N in these headwaters is from ground-water discharge. Consequently, the concentrations and forms of N in groundwater are important factors affecting major aquatic ecosystems; despite this, few data exist for several species of N in groundwater and controls on speciation are ill-defined. Herein, we report N speciation for a spring and well that were selected to reflect agricultural impacts, and a spring and well that show little to no agricultural-N impact. Samples were characterized for NO₃⁻, NO₂⁻, N₂O, NH₄⁺, urea, particulate organic N(), and dissolved organic N(). These analytes were monitored in the agricultural spring for up to two years along with other analytes that we reported upon previously. For all samples, when oxidized N was present, the dominant species was NO₃⁻ (88-98% of total fixed N pool) followed by (<4-12%) and only trace fractions of the other N analytes. In the non-agriculturally impacted well sample, which had no quantifiable NO₃⁻ or dissolved O₂, comprised the dominant fraction (68%) followed by NH₄⁺ (32%), with only a trace balance comprised of other N analytes. Water drawn from the well, spring and a wetland situated in the agricultural watershed also were analyzed for dissolved N₂ and found to have a fugacity in excess of that of the atmosphere. H₂O₂ was analyzed in the agricultural spring to evaluate the O₂/H₂O₂ redox potential and compare it to other calculated potentials. The potential of the O₂/H₂O₂ couple was close in value to the NO₃⁻/NO₂⁻ couple suggesting the important role of H₂O₂ as an O₂-reduction intermediate product and that O₂ and NO₃⁻ are reduced concomitantly. The O₂/H₂O₂ and NO₃⁻/NO₂⁻ couples also were close in value to a cluster of other inorganic N and Fe couples indicating near partial equilibrium among these species. Urea mineralization to NO₂⁻ was found to approach equilibrium with the reduction of O₂ to H₂O₂. By modeling as amide functional groups, as justified by recent analytical work, similar thermodynamic calculations support that mineralization to NO₂⁻ proceeds nearly to equilibrium with the reduction of O₂ to H₂O₂ as well. This near equilibration of redox couples for urea- and -oxidation with O₂-reduction places these two couples within the oxidized redox cluster that is shared among several other couples we have reported previously. In the monitored agricultural spring, [NO₃⁻] was lower in the summer than at other times, whereas [N₂O] was higher in the summer than at other times, perhaps reflecting a seasonal variation in the degree of denitrification reaction progress. No other N analytes were observed to vary seasonally in our study. In the well having no agricultural-N impact, C_org/N_org = 5.5, close to the typical value for natural aqueous systems of about 6.6. In the agricultural watershed C_org/N_org varied widely, from ∼1.2 to ?9.     

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.     

A hydrothermal investigation of the system FeO, Fe2O3, TiO2

Phase relations in the system FeO---Fe₂O₃---TiO₂, at temperatures ranging between 300°C and 700°C, have been investigated experimentally with special refference to the reaction Fe₃O₄ + TiO₂ = Fe₂O₃ + FeTiO₃. Pressure was varied between 500 and 2000 bars but its effect was negligible. Magnetite and rutile are the stable assemblage at temperatures above 550 dgC, and hematite and ilmenite are stable for lower temperatures. The equilibrium oxygen fugacity is estimated to be 10^−17.5 bars at equilibrium temperature. It is suggested that intermediate hematite-ilmenite solid solutions are inhomogeneous, consisting of ‘domains’ of hematite and ilmenite. The ‘domains’ are too small to be resolved by X-ray diffraction techniques. The top of the solvus curve in the hematite-ilmenite solution corresponds to a temperature of 660°C. Regular solution theory is not applicable to the solid solution.