Thermodynamic modeling of binary CH4-H2O fluid inclusions

This work reports the application of thermodynamic models, including equations of state, to binary (salt-free) CH₄-H₂O fluid inclusions. A general method is presented to calculate the compositions of CH₄-H₂O inclusions using the phase volume fractions and dissolution temperatures of CH₄ hydrate. To calculate the homogenization pressures and isolines of the CH₄-H₂O inclusions, an improved activity-fugacity model is developed to predict the vapor-liquid phase equilibrium. The phase equilibrium model can predict methane solubility in the liquid phase and water content in the vapor phase from 273 to 623 K and from 1 to 1000 bar (up to 2000 bar for the liquid phase), within or close to experimental uncertainties. Compared to reliable experimental phase equilibrium data, the average deviation of the water content in the vapor phase and methane solubility in the liquid phase is 4.29% and 3.63%, respectively. In the near-critical region, the predicted composition deviations increase to over 10%. The vapor-liquid phase equilibrium model together with the updated volumetric model of homogenous (single-phase) CH₄-H₂O fluid mixtures (Mao S., Duan Z., Hu J. and Zhang D. (2010) A model for single-phase PVTx properties of CO₂-CH₄-C₂H₆-N₂-H₂O-NaCl fluid mixtures from 273 to 1273 K and from 1 to 5000 bar. Chem. Geol.275, 148-160), is applied to calculate the isolines, homogenization pressures, homogenization volumes, and isochores at specified homogenization temperatures and compositions. Online calculation is on the website: http://www.geochem-model.org/.     

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.     

Equation of state of the H2O, CO2, and H2O-CO2 systems up to 10 GPa and 2573.15 K: Molecular dynamics simulations with ab initio potential surface

Based on our previous development of the molecular interaction potential for pure H₂O and CO₂ [Zhang, Z.G., Duan, Z.H. 2005a. Isothermal-isobaric molecular dynamics simulations of the PVT properties of water over wide range of temperatures and pressures. Phys. Earth Planet Interiors149, 335-354; Zhang, Z.G., Duan, Z.H. 2005b. An optimized molecular potential for carbon dioxide. J. Chem. Phys.122, 214507] and the ab initio potential surface across CO₂-H₂O molecules constructed in this study, we carried out more than one thousand molecular dynamics simulations of the PVTx properties of the CO₂-H₂O mixtures in the temperature-pressure range from 673.15 to 2573.15 K up to 10.0 GPa. Comparison with extensive experimental PVTx data indicates that the simulated results generally agree with experimental data within 2% in density, equivalent to experimental uncertainty. Even the data under the highest experimental temperature-pressure conditions (up to 1673 K and 1.94 GPa) are well predicted with the agreement within 1.0% in density, indicating that the high accuracy of the simulation is well retained as the temperature and pressure increase. The consistent and stable predictability of the simulation from low to high temperature-pressure and the fact that the molecular dynamics simulation resort to no experimental data but to ab initio molecular potential makes us convinced that the simulation results should be reliable up to at least 2573 K and 10 GPa with errors less than 2% in density. In order to integrate all the simulation results of this study and previous studies [Zhang and Duan, 2005a, 2005b] and the experimental data for the calculation of volumetric properties (volume, density, and excess volume), heat properties, and chemical properties (fugacity, activity, and possibly supercritical phase separation), an equation of state (EOS) is laboriously developed for the CO₂, H₂O, and CO₂-H₂O systems. This EOS reproduces all the experimental and simulated data covering a wide temperature and pressure range from 673.15 to 2573.15 K and from 0 to 10.0 GPa within experimental or simulation uncertainty.     

Equations of state for the NaCl-H2O-CH4 system and the NaCl-H2O-CO2-CH4 system: Phase equilibria and volumetric properties above 573 k

An equation of state (EOS) based on thermodynamic perturbation theory is presented for the NaCl-H₂O-CH₄ system. This equation consistently reproduces PvTX properties and phase equilibria with an accuracy close to that of data in the temperature, pressure and concentration ranges from 648 K to 873 K, 0 to 2500 bar and up to 2.37 mol % NaCl. Good agreement with recent ternary immiscibility data from 673 K to 873 K suggests that the EOS may provide accurate predictions for NaCl concentrations as high as 40 mol %. We could not find any experimental data above 873 K that can be used to validate the predictions of the EOS inside the ternary. However, parameters for the mixed ternary system were established from parameters evaluated for pure and binary systems and accurate combination rules. Therefore, predictions in the ternary should be reliable to the high temperatures and pressures where the EOS for the lower order systems are valid (about 1300 K and 5000 bar). Using the same combining approach, an EOS for the quaternary NaCl-H₂O-CO₂-CH₄ is constructed on the basis of parameters from our earlier model for the NaCl-H₂O-CO₂ system and the present NaCl-H₂O-CH₄ model. This suggests that predictions of the quaternary EOS are reliable also to about 1300 K and 5000 bar.     

Molecular dynamics simulation of the CH4 and CH4-H2O systems up to 10 GPa and 2573 K

Based on our previous study of the intermolecular potential for pure H₂O and the strict evaluation of the competitive potential models for pure CH₄ and the ab initio fitting potential surface across CH₄-H₂O molecules in this study, we carried out more than two thousand molecular dynamics simulations for the PVTx properties of pure CH₄ and the CH₄-H₂O mixtures up to 2573 K and 10 GPa. Comparison of 1941 simulations with experimental PVT data for pure CH₄ shows an average deviation of 0.96% and a maximum deviation of 2.82%. The comparison of the results of 519 simulations of the mixtures with the experimental measurements reveals that the PVTx properties of the CH₄-H₂O mixtures generally agree with the extensive experimental data with an average deviation of 0.83% and 4% in maximum, which is equivalent to the experimental uncertainty. Moreover, the maximum deviation between the experimental data and the simulation results decreases to about 2% as temperature and pressure increase, indicating that the high accuracy of the simulation is well retained in the high temperature and pressure region.After the validation of the simulation method and the intermolecular potential models, we systematically simulated the PVTx properties of this binary system from 673 K and 0.05 GPa to 2573 K and 10 GPa. In order to integrate all the simulation results and the experimental data for the calculation of thermodynamic properties, an equation of state (EOS) is developed for the CH₄-H₂O system covering 673-2573 K and 0.01-10 GPa. Isochores for compositions <4 mol% CH₄ up to 773 K and 600 MPa are also determined in this paper. The program for the EOS can be downloaded from www.geochem-model.org/programs.htm.     

Methane and CO2 sorption and desorption measurements on dry Argonne premium coals: pure components and mixtures

Sorption and desorption behaviour of methane, carbon dioxide, and mixtures of the two gases has been studied on a set of well-characterised coals from the Argonne Premium Coal Programme. The coal samples cover a maturity range from 0.25% to 1.68% vitrinite reflectance. The maceral compositions were dominated by vitrinite (85% to 91%). Inertinite contents ranged from 8% to 11% and liptinite contents around 1% with one exception (Illinois coal, 5%). All sorption experiments were performed on powdered (−100 mesh), dry coal samples.Single component sorption/desorption measurements were carried out at 22 °C up to final pressures around 51 bar (5.1 MPa) for CO₂ (subcritical state) and 110 bar (11 MPa) for methane.The ratios of the final sorption capacities for pure CO₂ and methane (in molar units) on the five coal samples vary between 1.15 and 3.16. The lowest ratio (1.15) was found for the North Dakota Beulah-Zap lignite (VR_r=0.25%) and the highest ratios (2.7 and 3.16) were encountered for the low-rank coals (VR_r 0.32% and 0.48%) while the ratio decreases to 1.6–1.7 for the highest rank coals in this series.Desorption isotherms for CH₄ and CO₂ were measured immediately after the corresponding sorption isotherms. They generally lie above the sorption isotherms. The degree of hysteresis, i.e. deviation of sorption and desorption isotherms, varies and shows no dependence on coal rank.Adsorption tests with CH₄/CO₂ mixtures were conducted to study the degree of preferential sorption of these two gases on coals of different rank. These experiments were performed on dry coals at 45 °C and pressures up to 180 bar (18 MPa). For the highest rank samples of this sequence preferential sorption behaviour was “as expected”, i.e. preferential adsorption of CO₂ and preferential desorption of CH₄ were observed. For the low rank samples, however, preferential adsorption of CH₄ was found in the low pressure range and preferential desorption of CO₂ over the entire pressure range.Follow-up tests for single gas CO₂ sorption measurements consistently showed a significant increase in sorption capacity for re-runs on the same sample. This phenomenon could be due to extraction of volatile coal components by CO₂ in the first experiment. Reproducibility tests with methane and CO₂ using fresh sample material in each experiment did not show this effect.     

CO2-H2O mixtures in the geological sequestration of CO2. I. Assessment and calculation of mutual solubilities from 12 to 100°C and up to 600 bar

Evaluating the feasibility of CO₂ geologic sequestration requires the use of pressure-temperature-composition (P-T-X) data for mixtures of CO₂ and H₂O at moderate pressures and temperatures (typically below 500 bar and below 100°C). For this purpose, published experimental P-T-X data in this temperature and pressure range are reviewed. These data cover the two-phase region where a CO₂-rich phase (generally gas) and an H₂O-rich liquid coexist and are reported as the mutual solubilities of H₂O and CO₂ in the two coexisting phases. For the most part, mutual solubilities reported from various sources are in good agreement. In this paper, a noniterative procedure is presented to calculate the composition of the compressed CO₂ and liquid H₂O phases at equilibrium, based on equating chemical potentials and using the Redlich-Kwong equation of state to express departure from ideal behavior. The procedure is an extension of that used by King et al. (1992), covering a broader range of temperatures and experimental data than those authors, and is readily expandable to a nonideal liquid phase. The calculation method and formulation are kept as simple as possible to avoid degrading the performance of numerical models of water-CO₂ flows for which they are intended. The method is implemented in a computer routine, and inverse modeling is used to determine, simultaneously, (1) new Redlich-Kwong parameters for the CO₂-H₂O mixture, and (2) aqueous solubility constants for gaseous and liquid CO₂ as a function of temperature. In doing so, mutual solubilities of H₂O from 15 to 100°C and CO₂ from 12 to 110°C and up to 600 bar are generally reproduced within a few percent of experimental values. Fugacity coefficients of pure CO₂ are reproduced mostly within one percent of published reference data.     

利用量子力学势能数据预测温度、压力、盐度和沉积物孔径对甲烷水合物形成和分解的影响

介绍一个预测不同温度、压力、盐度和沉积物毛细管孔径条件下甲烷水合物 溶液 气体多相平衡模型。该模型以Van der Waals和 Platteeuw热力学模型、量子力学从头算粒子相互作用势能、DMW 92状态方程和Pitzer电解质理论为基础,能在很宽广温压范围内预测温度、压力、盐度和毛细管力对甲烷水合物形成和分解的影响。通过对比本模型的预测结果与实验数据,可知本模型能够准确地预测海水和多孔介质中甲烷水合物的相平衡条件。对于一定盐度下多孔介质中甲烷水合物的形成温压条件的在线计算可浏览: www.geochem model.org/models.htm。     

Re-equilibration of CO2 fluid inclusions at controlled hydrogen fugacities

Abstract Natural, pure CO₂ inclusions in quartz and olivine (c. Fo₉₀) were exposed to controlled f_H2 conditions at T= 718–728°C and P_total= 2 kbar; their compositions were monitored (before and after exposures) by microsampling Raman spectroscopy (MRS) and microthermometry. In both minerals exposed at the graphite–methane buffer (f_H2= 73 bar), fluid speciations record the diffusion of hydrogen into the inclusions. In quartz, room-temperature products in euhedral isolated (EI type) inclusions are carbonic phases with molar compositions of c. CO₂(60) + CH₄(40) plus graphite (Gr) and H₂O, whereas anhedral inclusions along secondary fractures (AS type) are Gr-free and contain H₂O plus carbonic phases with compositions in the range c. CO₂(60) + CH₄(40) to CO₂(10) + CH₄(90). EI type inclusions in olivine evolved to c. CO₂(90–95) + CH₄(5–10) without Gr, whereas AS type inclusions have a range of compositions from CO₂(90) + CH₄(10) ± Gr to CH₄(50) + H₂(50) ± Gr; neither H₂O nor any hydrous species was detected by optical microscopy or MRS in the olivine-hosted products. Differences in composition between and among the texturally distinct populations of inclusions in both minerals probably arise from variations in initial fluid densities, as all inclusions apparently equilibrated with the ambient f_H2. These relations suggest that compositional variability among inclusions in a given natural sample does not require the entrapment of multiple generations of fluids. In addition, the absence of H₂O in the olivine-hosted inclusions would require the extraction of oxygen from the fluids, in which case re-equilibration mechanisms may be dependent on the composition and structure of the host mineral. Many of the same samples were re-exposed to identical P–T conditions using Ar as the pressure medium, yielding ambient f_H2= 0.06 bar. In most inclusions, the carbonic fluids returned to pure CO₂ and graphite persisted in the products. Reversal of the mechanisms from the prior exposure at f_H2= 73 bar did not occur in any inclusions but the AS types in olivine, in which minor CO₂ was produced at the expense of CH₄ and/or graphite. The observed non-reversibility of previous mechanisms may be attributed to: (1) slower fluid–solid reactions compared to reactions in the homogeneous fluid phase; (2) depressed activities of graphite due to poor ordering; and/or (3) low ambient f_O2 at the conditions of the second run.     

Chemical and isotopic equilibrium between CO2 and CH4 in fumarolic gas discharges: Generation of CH4 in arc magmatic-hydrothermal systems

The chemical and isotopic composition of fumarolic gases emitted from Nisyros Volcano, Greece, and of a single gas sample from Vesuvio, Italy, was investigated in order to determine the origin of methane (CH₄) within two subduction-related magmatic-hydrothermal environments.Apparent temperatures derived from carbon isotope partitioning between CH₄ and CO₂ of around 340°C for Nisyros and 470°C for Vesuvio correlate well with aquifer temperatures as measured directly and/or inferred from compositional data using the H₂O-H₂-CO₂-CO-CH₄ geothermometer. Thermodynamic modeling reveals chemical equilibrium between CH₄, CO₂ and H₂O implying that carbon isotope partitioning between CO₂ and CH₄ in both systems is controlled by aquifer temperature.N₂/³He and CH₄/³He ratios of Nisyros fumarolic gases are unusually low for subduction zone gases and correspond to those of midoceanic ridge environments. Accordingly, CH₄ may have been primarily generated through the reduction of CO₂ by H₂ in the absence of any organic matter following a Fischer-Tropsch-type reaction. However, primary occurrence of minor amounts of thermogenic CH₄ and subsequent re-equilibration with co-existing CO₂ cannot be ruled out entirely. CO₂/³He ratios and δ¹³C_CO2 values imply that the evolved CO₂ either derives from a metasomatized mantle or is a mixture between two components, one outgassing from an unaltered mantle and the other released by thermal breakdown of marine carbonates. The latter may contain traces of organic matter possibly decomposing to CH₄ during thermometamorphism.     

Stable carbon isotope composition and concentrations of CO2 and CH4 in the deep catotelm of a peat bog

Vertical profiles of concentration and C-isotopic composition of dissolved methane and carbon dioxide were observed over 26 months in the catotelm of a deep (6.5 m) peat bog in Switzerland. The dissolved concentrations of these gases increase with depth while CO₂ predominates over CH₄ (CO₂ ca. 5 times CH₄). This pattern can be reproduced by a reaction-advection-ebullition model, where CO₂ and CH₄ are formed in a ratio of 1:1. The less soluble methane is preferentially lost via outgassing (bubbles). The isotopic fractionation between CO₂ and CH₄ also increases with depth, with α_C values ranging from 1.045 to 1.075. The isotopic composition of the gases traces the passage of respiration-derived CO₂ (from the near surface) through a shallow zone with methanogenesis of low isotopic fractionation (splitting of fermentation-derived acetate). This solution then moves through the catotelm, where methanogenesis occurs by CO₂ reduction (large isotopic fractionation). In the upper part of the catotelm the C-13-depleted respiration-derived CO₂ pool buffers the isotopic composition of CO₂; the δ¹³C of CO₂ increases only slowly. At the same time strongly depleted CH₄ is formed as CO₂ reduction consumes the depleted CO₂. In the lower part of the catotelm, the respiration-derived CO₂ and shallow CH₄ become less important and CO₂ reduction is the dominant source of CO₂ and CH₄. Now, the δ¹³C values of both gases increase until equilibrium is reached with respect to the isotopic composition of the substrate. Thus, the δ¹³C values of methane reach a minimum at intermediate depth, and the deep methane has δ¹³C values comparable to shallow methane. A simple mixing model for the isotopic evolution is suggested. Only minor changes of the observed patterns of methanogenesis (in terms of concentration and isotopic composition) occur over the seasons. The most pronounced of these is a slightly higher rate of acetate splitting in spring.     

Characterization of CO2CH4H2O fluid inclusions by microthermometry and laser Raman microprobe spectroscopy: Inferences for clathrate and fluid equilibria

Microthermometry (MT) and laser Raman microprobe (LRM) spectroscopy (at room temperature and at about 0°C) were done on 33 synthetically produced CO₂CH₄H₂O fluid inclusions in quartz (from R. Bodnar and M. Sterner). At room temperature, the inclusions consist of an aqueous liquid, a CO₂CH₄ supercritical carbonic fluid and (in most cases) graphite. In all these inclusions, the melting temperature for solid CO₂ is less than that for the homogenization of the vapor bubble in the carbonic fluid.A method is described whereby MT data for CO₂CH₄H₂O inclusions can be projected within the CO₂CH₄ binary phase diagram to infer CO₂:CH₄ ratios in the carbonic fluid (<2 to > 35 mole% CH₄ in the inclusions under study). This method takes into account the formation of CO₂CH₄ clathrate hydrate during MT analysis. Unless clathrate formation is properly considered, errors arise in the determination of the bulk CO₂CH₄ ratio. For the inclusions in our study, these errors are on the order of 5 to 8 mole% CH₄. Our interpretation of the MT data indicates that CH₄ is preferentially partitioned into the clathrate over the coexisting carbonic fluid, in contradiction to the prediction from Parrish and Prausnitz's (1972) model for clathrate equilibria. Comparison of LRM analyses on the bulk carbonic fluid (in the absence of clathrate) and the residual carbonic fluid (in the presence of clathrate) confirm the preferential partitioning of CH₄ into the clathrate. LRM analyses of the clathrate itself indicate that CH₄ occupies both types of cage sites in the clathrate structure, whereas CO₂ may only occupy one site. Two by-products of the combined LRM and MT analyses of the same inclusions are derivation of empirical ratios of Raman quantification factors for high-density CO₂CH₄ fluids and the ability to determine CO₂:CH₄ ratios of inclusions whose MT data lie near the critical region for CO₂CH₄. Thus, the joint use of LRM and MT techniques provides information that could not be obtained by either technique alone.     

Prediction of vapor-liquid equilibrium and PVTx properties of geological fluid system with SAFT-LJ EOS including multi-polar contribution. Part I: Application to H2O-CO2 system

Molecular based equations of state (EOS) are attractive because they can take into account the energetic contribution of the main types of molecular interactions. This study models vapor-liquid equilibrium (VLE) and PVTx properties of the H₂O-CO₂ binary system using a Lennard-Jones (LJ) referenced SAFT (Statistical Associating Fluid Theory) EOS. The improved SAFT-LJ EOS is defined in terms of the residual molar Helmholtz energy, which is a sum of four terms representing the contributions from LJ segment-segment interactions, chain-forming among the LJ segments, short-range associations and long-range multi-polar interactions. CO₂ is modeled as a linear chain molecule with a constant quadrupole moment, and H₂O is modeled as a spherical molecule with four association sites and a dipole moment. The multi-polar contribution to Helmholtz energy, including the dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole contribution for H₂O-CO₂ system, is calculated using the theory of Gubbins and Twu (1978). Six parameters for pure H₂O and four parameters for pure CO₂ are needed in our model. The Van der Waals one-fluid mixing rule is used to calculate the Lennard-Jones energy parameter and volume parameter for the mixture. Two or three binary parameters are needed for CO₂-H₂O mixtures, which are evaluated from phase equilibrium data of the binary system. Comparison with the experimental data shows that our model represents the PVT properties of CO₂ better than other SAFT EOS without a quadrupole contribution. For the CO₂-H₂O system, our model agrees well with the vapor-liquid equilibrium data from 323-623 K. The average relative deviation for CO₂ solubility (expressed in mole fraction) in water is within 6%. Our model can also predict the PVTx properties of CO₂-H₂O mixtures up to 1073 K and 3000 bar. The good performance of this model indicates that: (1) taking account of the multi-polar contribution explicitly improves the agreement of calculated properties with experimental data at high temperatures and high pressures, (2) the molecular-based EOS with just a few parameters fit to data in the sub-critical region can predict the thermodynamic properties of fluids over a wide range of P-T conditions.     

An algorithm for finding composition,molar volume and isochors of CO2-CH4 fluid inclusions from Th and Tfm (for Th < Tfm)

A modified Redlich-Kwong equation of state is used to calculate the solubility of CO₂ in methane at various temperatures and pressures. From the solubility of CO₂ in CH₄ at the triple point and at final melting (T_h < T_fm), and the molar volume of solid CO₂, the volume of solid at the triple point, and the molar volume of the inclusion can be calculated using a mass balance. The pressure at the melting point is calculated from the equation of state.The algorithm predicts composition, molar volume, pressure at final melting and the isochor pressure (for a given temperature of trapping) for CO₂-CH₄ fluid inclusions for the case T_h < T_fm, given T_h, T_fm and experimental data on P_h and d_co2 (solid) at T_h.     

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.     

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压力条件下与压力及组成无关。地质样品应用结果表明,本方法可以方便、灵活、准确地按任意比例将两瓶及两瓶以上纯气体钢瓶样品进行混合及增压,为激光拉曼标定、气体组成原位测量等提供了一种新的技术思路。     

汶川M_s 8.0地震破裂带CO_2、CH_4、Rn和Hg脱气强度

地震活动断裂带能够向大气释放大量的温室气体、放射性气体和有毒气体(CO_2、CH_4、Rn和Hg),并对大气环境的影响产生复杂的影响。利用静态暗箱法,对汶川M_s8.0地震破裂带CO_2、Rn和Hg脱气强度进行实地测量,并计算了CO_2和Hg脱气对大气的年贡献量。结果表明:(1)破裂带土壤气中CO_2、CH_4、Rn和Hg异常浓度最大值分别可以达到7.98%、2.38%、524.30k Bq/m~3和161.00ng/m~3;破裂带CO_2、Rn和Hg脱气平均通量是34.95g·m~(-2)d~(-1)、36.11m Bq·m~(-2)s~(-1)和26.56ng·m~(-2)h~(-1),最大值分别达到259.23g·m~(-2)d~(-1)、580.35m Bq·m~(-2)s~(-1)和387.67ng·m~(-2)h~(-1);(2)汶川Ms8.0地震破裂带向大气脱气的CO_2年贡献量是0.95Mt,Hg的年贡献量是15.94kg。汶川Ms8.0地震破裂带破裂CO_2、CH_4、Rn和Hg等的脱气强度,不仅与破裂带渗透率有关,还与断裂带浅部存在的气藏、煤层以及磷矿层等气体源有重要的联系。     

Flue gas and pure CO2 sorption properties of coal: A comparative study

Presently many research projects focus on the reduction of anthropogenic CO₂ emissions. It is intended to apply underground storage techniques such as flue gas injection in unminable coal seams. In this context, an experimental study has been performed on the adsorption of pure CO₂ and preferential sorption behavior of flue gas. A coal sample from the Silesian Basin in Poland (0.68% V Rr), measured in the dry and wet state at 353 K has been chosen for this approach. The flue gas used was a custom class industrial flue gas with 10.9% of CO₂, 0.01% of CO, 9% of H₂, 3.01% of CH₄, 3.0% of O₂, 0.106% of SO₂ and nitrogen as balance.Adsorption isotherms of CO₂ and flue gas were measured upto a maximum of 11 MPa using a volumetric method. Total excess sorption capacities for CO₂ on dry and wet Silesia coal ranged between 1.9 and 1.3 mmol/g, respectively. Flue gas sorption capacities on dry and wet Silesia coal were much lower and ranged between 0.45 and 0.2 mmol/g, respectively, at pressures of 8 MPa. The low sorption capacity of wet coal has resulted from water occupying some of the more active adsorption sites and hence reducing the heterogeneity of adsorption sites relative to that of dry coal. Desorption tests with flue gas were conducted to study the degree of preferential sorption of the individual components. These experiments indicate that CO₂ is by far the prefered sorbing component under both wet and dry conditions. This is followed by CH₄. N₂ adsorbs very little on the coal in the presence of CO₂ and CH₄. It is also observed that the adsorption of CO₂ onto coal is not significantly hindered by the addition of other gases, other than dilution effect of the pressure.In addition to the sorption experiments, the density of the flue gas mixture has been determined up to 18 MPa at 318 K. A very good precision of these measurements were documented by volumetric methods.     

Carbon-13 exchange between CO2 and CH4 under geothermal conditions

On the basis of recently reported data on the kinetics of carbon-13 exchange between CO₂ and CH₄ at temperatures above 500°C, first order rate constants $\text{log k = 11.16?10,190/T}$ were derived allowing variations in Δ, the difference in the isotopic composition of coexisting CO₂ and CH₄, to be evaluated as a function of initial composition and cooling rate of the rising geothermal fluid. Observed Δ-values in geothermal discharges are likely to represent frozen in compositions attained after minimum residence times of 20 ka at 400°C or 10 Ma at 300°C. The carbon-13 contents of any biogenic gases are unlikely to have been affected by thermal re-equilibration at temperatures below 200°C. The chemical equilibrium involving CO₂ and CH₄ can be expected to proceed about a hundred times faster than isotopic equilibration.     

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.