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

C02-CH4流体与金成矿作用：以阿尔泰山南缘和穆龙套金矿为例

造山型金矿的成矿作用与H20-CO2流体有着密切的联系。然而对阿尔泰山南缘和穆龙套金矿的流体包裹体研究表明，无水的CO2-CH4流体在中亚成矿域中一些金矿床中具有同样重要意义。阿尔泰山南缘萨热阔布金矿包裹体的Xch4达0.20～0.23，穆龙套金矿的XCH4为0.07～0.23。CH4扩大了流体不混溶的范围，有利于对Au的富集沉淀。CO2流体在Au成矿中的重要作用至少包括了三方面的意义，即：缓冲流体PH值范围、提高流体中的Au含量并使其维持与还原硫的络合作用进行迁移；扩大超临界流体的温度范围；增加流体不混溶的区域。CH4的加入扩大了流体不混溶的范围，有利于对Au的富集沉淀。     

Review of the systematics of CO2–H2O fluid inclusions

Aqueous solutions that contain volatile (gas) components are one of the most important types of fluid in the Earth's crust. The record that such fluids have left in the form of fluid inclusions in minerals provides a wealth of insight into the geochemical and petrologic processes in which the fluids participated. This article reviews the systematics of CO₂–H₂O fluid inclusions as a starting point for interpreting the chemically more complex systems. The phase relations of the binary are described with respect to a qualitative P–T–X model, and isoplethic–isochoric paths through this model are used to explain the equilibrium and non-equilibrium behaviour of fluid inclusions during microthermometric heating and cooling. The P–T–X framework is then used to discuss the various modes of fluid inclusion entrapment, and how the resulting assemblage textures can be used to interpret the P–T conditions, phase states, and evolution paths of the parent solutions. Finally, quantitative methods are reviewed by which bulk molar volume and composition of CO₂–H₂O fluid inclusions can be determined from microthermometric observations of phase transitions.     

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.     

Methane: An equation of state with application to the ternary system H2O-CO2-CH4

The following hardsphere modified Redlich-Kwong (HSMRK) equation of state was obtained by least squares fitting to available P-V-T data for methane (P in bars; T in Kelvins; v in cm³ mol^?1; b = 60.00 cm³ mol^?1; R = 83.14 cm³barmol^?1K^?1): $\text{P}\text{RT(1 + y + y}{}^2\text{?y}{}^3\text{v(1?y)}{}^3\text{)}\text{-}\text{c(T) + d(T)}\text{v}\text{+}\text{e(T)}\text{v}{}^2\text{/v(v + b)T}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$$\text{y =}\text{b}\text{4v}$c(T) = 13.403 × 10⁶ + (9.28 × 10⁴)T + 2.7 T²d(T) = 5.216 × 10⁹ ? (6.8 × 10⁶)T + (3.28 × 10³)T²e(T) = (?2.3322 × 10¹¹) + (6.738 × 10⁸)T + (3.179 × 10⁵)T² For the P-T range of experimental data used in the fit (50 to 8600 bars and from 320 to 670 K), calculated volumes and fugacity coefficients for CH₄ relative to experimentally determined volumes and fugacity coefficients have average percent deviations of 0.279 and 1.373, respectively. The HSMRK equation, which predicts linear isochores over a wide P-T range, should yield reasonable estimates of fugacity coefficients for CH₄ to pressures and temperatures well outside the P-T range of available P-V-T data. Calculations for the system H₂O-CO₂-CH₄, using the HSMRK equations for H₂O and CO₂ of Kerrick and Jacobs (1981) and the HSMRK equation for CH₄ of this study, indicate that compared to the binary H₂O-CO₂ system, small amounts of CH₄ in the ternary system H₂O-CO₂-CH₄ slightly increases the activity of H₂O, and significantly decreases the activity of CO₂.     

A general algorithm for the O abundance correction to C/C determinations from CO2 isotopologue measurements, including CO2 characterised by ‘mass-independent’ oxygen isotope distributions

It is widely recognised that a significant limitation to the ultimate precision of carbon stable isotope ratio measurements, as obtained from dual-inlet mass spectrometric measurements of CO₂ isotopologue ion abundances at m/z 44, 45, and 46, is the correction for interference from ¹⁷O-bearing molecular ions. Two long-established, alternative procedures for determining the magnitude of this correction are in widespread use (although only one has IAEA approval); their differences lead to small but potentially significant discrepancies in the magnitude of the resulting correction. Furthermore, neither approach was designed to accommodate oxygen three-isotope distributions which do not conform to terrestrial mass-dependent behaviour. Stratospheric CO₂, for example, contains a strongly ‘mass-independent’ oxygen isotope composition. A new strategy for determining the ¹⁷O-bearing ion correction is presented, for application where the oxygen three-isotope characteristics of the analyte CO₂ are accurately known (or assigned) in terms of the slope λ of the three-isotope fractionation line and the ordinate axis intercept 10³ ln(1 + k) on a 10³ ln(1 + δ¹⁷O) versus 10³ ln(1 + δ¹⁸O) plot. At the heart of the approach is the relationship between ¹⁷R, which is the ¹⁷O/¹⁶O ratio of the sample CO₂, and other assigned or empirically determined parameters needed for the δ¹³C evaluation:

     

四川冕宁木落稀土矿床稀土矿物中富CO2和H2流体包裹体研究

木落稀土矿床位于四川省冕宁县境内，其成矿与喜山期岩浆碳酸岩有关。通过对矿床中主要稀土矿物氟碳铈矿中流体包裹体的岩相学、包裹体显微测温分析和包裹体成分的LRM分析等，对成矿流体的特征、演化及稀土矿化过程进行了讨论，结果表明与稀土成矿有关的流体为富CO₂和SO₂ 的中高温、高压超临界流体。温度、压力降低和流体不混溶是造成稀土矿物沉淀的主要机制。     

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.     

Numerical model to determine the composition of H2O-NaCl-CaCl2 fluid inclusions based on microthermometric and microanalytical data

Natural fluids approximated by the H₂O-NaCl-CaCl₂ system are common in a wide range of geologic environments, including sedimentary basins associated with hydrocarbon occurrences and MVT deposits, submarine hydrothermal systems, and other metamorphic, magmatic and hydrothermal environments. We present a comprehensive numerical model and Microsoft® Excel©-based computer program to determine the compositions of fluid inclusions in the H₂O-NaCl-CaCl₂ system based on microthermometric and microanalytical data. The model consists of six polynomial correlation equations that describe liquid salinity as a function of NaCl/CaCl₂ ratio and melting temperature on each of the ice, hydrohalite, halite, antarcticite, CaCl₂·4H₂O and CaCl₂·2H₂O vapor-saturated liquidus surfaces. The cotectic and peritectic boundaries are determined from the intersections of the liquidus surfaces. The model is implicitly internally consistent and topologically correct.The model expands upon the compositional range of applicability and the data types that can be used for compositional determination. It reproduces experimental data for all compositions that lie within the H₂O-NaCl-CaCl₂·4H₂O compositional triangle in the H₂O-NaCl-CaCl₂ system and yields accurate reproductions of the H₂O-NaCl and H₂O-CaCl₂ binaries. Furthermore, in comparison to previously published models, the one presented here eliminates systematic errors, wavy isotherms and cotectic and peritectic curves with local “bumps.”     

宁洱火山区壳内岩浆囊现今温度:来自CO_2-CH_4碳同位素地质温度计的估计

温度是岩浆囊的重要物理参数,获取温度参数并监测其变化对更好地理解岩浆囊的物理化学性质和行为、评价火山的活动性和喷发危险性具有重要的理论和现实意义。火山岩浆-水热系统丰富的含碳气体,这些含碳气体间碳同位素温度计为我们解决休眠火山区地下岩浆囊的温度问题提供了可能。我们对宁洱火山区地表地热异常强度最高区的4个温泉点CO_2、CH_4碳同位素组成进行了2年的连续观测和2次平行观测,利用观测数据对CO_2、CH_4气体的岩浆来源进行了甄别,对CO_2-CH_4间碳同位素分馏是否发到平衡进行了判断,然后利用Horita CO_2-CH_4碳同位素平衡分馏方程计算了宁洱-通关火山区逸出气体的源区温度。结果表明:宁洱火山区4个样点中,1个样点的CO_2和CH_4均为非岩浆成因,其他3个样点的主要属于岩浆来源。剔除了非岩浆来源的数据后的δ-Δ图解法判断CO_2-CH_4间碳同位素分馏达到了平衡。3个样点的CO_2、CH_4碳气体源区温度分别为425~475℃、941~995℃和1179~1578℃。结合最上部地壳温度场结果,我们认为,宁洱火山区现今存在2个壳内岩浆房,分别位于宁洱火山以南和通关火山以北,宁洱岩浆房的温度为1179~1578℃,通关岩浆房的温度为941~995℃。2个岩浆房的现今温度已达到流纹岩浆(600~900℃)、安山岩(800~1100℃)和玄武岩浆(1000~1250℃)的形成温度。思普地震带空间上密集的6级地震丛集活动可以用宁洱岩浆房的高活动性来解释。δ-Δ图解法判断CO_2-CH_4间碳同位素分馏平衡准则应修正为:在保持两拟合直线的斜率符号相反的条件下,δ~(13)CCO_2-ΔCO_2-CH_4拟合直线和δ~(13)CCH_4-ΔCO_2-CH_4拟合直线应相交于δ轴附近截距差Δb≤0.16处。     

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.     

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.     

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.     

Chemistry of fluids from a natural analogue for a geological CO2 storage site (Montmiral, France): Lessons for CO2–water–rock interaction assessment and monitoring

Chemical and isotope studies of natural CO₂ accumulations aid in assessing the chemical effects of CO₂ on rock and thus provide a potential for understanding the long-term geochemical processes involved in CO₂ geological storage. Several natural CO₂ accumulations were discovered during gas and oil exploration in France’s carbogaseous peri-Alpine province (south-eastern France) in the 1960s. One of these, the Montmiral accumulation at a depth of more than 2400 m, is currently being exploited. The chemical composition of the water collected at the wellhead has changed in time and the final salinity exceeds 75 g/L. These changes in time can be explained by assuming that the fraction of the reservoir brine in the recovered brine–CO₂–H₂O mixture varies, resulting in variable proportions of H₂O and brine in the sampled water. The proportions can be estimated in selected samples due to the availability of gas and water flowrate data. These data enabled the reconstruction of the chemical and isotope composition of the brine. The proportions of H₂O and brine can also be estimated from isotope (δ²H, δ¹⁸O) composition of collected water and δ¹⁸O of the sulfates or CO₂. The reconstituted brine has a salinity of more than 85 g/L and, according to its Br⁻ content and isotope (δ²H, δ¹⁸O, δ³⁴S) composition, originates from an evaporated Triassic seawater that underwent dilution by meteoric water. The reconstitution of the brine’s chemical composition enabled an evaluation of the CO₂–water–rock interactions based on: (1) mineral saturation indices; and (2) comparison with initial evaporated Triassic seawater. Dissolution of K- and SO₄-containing minerals such as K-feldspar and anhydrite, and precipitation of Ca and Mg containing minerals that are able to trap CO₂ (carbonates) are highlighted. The changes in concentration of these elements in the brine, which are attributed to CO₂ interactions, illustrate the relevance of monitoring the water quality at future industrial CO₂ storage sites.     

Solubility of Fe-ettringite (Ca6[Fe(OH)6]2(SO4)3 · 26H2O)

The solubility of Fe-ettringite (Ca₆[Fe(OH)₆]₂(SO₄)₃ · 26H₂O) was measured in a series of precipitation and dissolution experiments at 20 °C and at pH-values between 11.0 and 14.0 using synthesised material. A time-series study showed that equilibrium was reached within 180 days of ageing. After equilibrating, the solid phases were analysed by XRD and TGA while the aqueous solutions were analysed by ICP-OES (calcium, sulphur) and ICP-MS (iron). Fe-ettringite was found to be stable up to pH 13.0. At higher pH-values Fe-monosulphate (Ca₄[Fe(OH)₆]₂(SO₄) · 6H₂O) and Fe-monocarbonate (Ca₄[Fe(OH)₆]₂(CO₃) · 6H₂O) are formed. The solubilities of these hydrates at 25 °C are:      

阿勒泰恰夏铜矿床的富CO2流体与矿床成因

恰夏铜矿床位于新疆阿尔泰山南缘克兰火山-沉积盆地内,赋矿地层主要为下泥盆统康布铁堡组上亚组变质岩系.脉状铜矿化主要特征为:早期顺层石英脉,呈脉状或透镜状沿变质片理分布,有星点状黄铁矿产出;晚期含铜黄铁矿-石英脉,斜切变质围岩,黄铜矿以浸染状分布于石英脉裂隙中.石英脉中流体包裹体主要为富CO2包裹体,其次为水溶液包裹体,同时含有少量的碳质流体包裹体.显微测温研究表明,早期顺层石英中原生富CO2流体包裹体,CO2三相点温度(tm,CO2)集中在-61.5～-57.5℃,CO2部分均一温度(th,CO2)集中在25～27℃,完全均一温度(th,tot)集中于223～280℃,流体密度为0.82～0.90 g/cm3;含铜黄铁矿-石英脉中原生富CO2包裹体的tm,CO2集中于-61.5～-58.7℃,th,CO2集中在23.5～28.7℃,th,tot集中在230 ～ 310℃,流体密度0.81～0.86 g/cm3.成矿流体为中高温、中低盐度、富CO2的CO2-H2O-NaCl±CH4±N2体系.恰夏铜矿脉状铜矿化的成矿流体特征与造山型金矿床的流体包裹体特征类似,结合矿床产出的地质背景、控矿构造特征,认为脉状铜矿化的成因与造山-变质热液有关,是阿尔泰山南缘晚泥盆世一二叠纪造山-变质作用的产物.SRXRF测试富CO2流体包裹体中金属微量元素,显示其富集Au,可能表明富CO2流体对金的富集起到一定作用.     

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.