首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Based on our previous development of the molecular interaction potential for pure H2O and CO2 [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 CO2-H2O molecules constructed in this study, we carried out more than one thousand molecular dynamics simulations of the PVTx properties of the CO2-H2O 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 CO2, H2O, and CO2-H2O 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.  相似文献   

2.
This work reports the application of thermodynamic models, including equations of state, to binary (salt-free) CH4-H2O fluid inclusions. A general method is presented to calculate the compositions of CH4-H2O inclusions using the phase volume fractions and dissolution temperatures of CH4 hydrate. To calculate the homogenization pressures and isolines of the CH4-H2O 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) CH4-H2O fluid mixtures (Mao S., Duan Z., Hu J. and Zhang D. (2010) A model for single-phase PVTx properties of CO2-CH4-C2H6-N2-H2O-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/.  相似文献   

3.
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 H2O-CO2 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. CO2 is modeled as a linear chain molecule with a constant quadrupole moment, and H2O 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 H2O-CO2 system, is calculated using the theory of Gubbins and Twu (1978). Six parameters for pure H2O and four parameters for pure CO2 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 CO2-H2O 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 CO2 better than other SAFT EOS without a quadrupole contribution. For the CO2-H2O system, our model agrees well with the vapor-liquid equilibrium data from 323-623 K. The average relative deviation for CO2 solubility (expressed in mole fraction) in water is within 6%. Our model can also predict the PVTx properties of CO2-H2O 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.  相似文献   

4.
Recent isopiestic studies of the Fe2(SO4)3-H2SO4-H2O system at 298.15 K are represented with an extended version of Pitzer’s ion interaction model. The model represents osmotic coefficients for aqueous {(1 − y)Fe2(SO4)3 + yH2SO4} mixtures from 0.45 to 3.0 m at 298.15 K and 0.0435 ? y ? 0.9370. In addition, a slightly less accurate representation of a more extended molality range to 5.47 m extends over the same y values, translating to a maximum ionic strength of 45 m. Recent isopiestic data for the system at 323.15 K are represented with the extended Pitzer model over a limited range in molality and solute fraction. These datasets are also represented with the usual “3-parameter” version of Pitzer’s model so that it may be incorporated in geochemical modeling software, but is a slightly less accurate representation of thermodynamic properties for this system. Comparisons made between our ion interaction model and available solubility data display partial agreement for rhomboclase and significant discrepancy for ferricopiapite. The comparisons highlight uncertainty remaining for solubility predictions in this system as well as the need for additional solubility measurements for Fe3+-bearing sulfate minerals. The resulting Pitzer ion interaction models provide an important step toward an accurate and comprehensive representation of thermodynamic properties in this geochemically important system.  相似文献   

5.
The position of the Raman methane (CH4) symmetric stretching band (ν1) over the range 1-650 bar and 0.3-22 °C has been determined using a high-pressure optical cell mounted on a Raman microprobe. Two neon emission lines that closely bracket the CH4 band were collected simultaneously with each CH4 spectrum. The peak position was determined after least squares fitting using a summed Gaussian-Lorentzian method, resulting in a precision of ≈±0.02 cm−1 in peak position determination. The CH4ν1 band position shifts to lower wave number with increasing pressure. At a given pressure, the band shifts to lower wave number with decreasing temperature, and the magnitude of the temperature shift increases with increasing pressure. The relationship between the Raman CH4ν1 band position and temperature and pressure determined here may be used to estimate the internal pressure in natural or synthetic CH4-bearing fluid inclusions. This information, in turn, may be used to determine the density of pure CH4 fluid inclusions and the salinity of CH4-bearing aqueous inclusions.  相似文献   

6.
A unified equation has been derived by using all available data for calculating methane vapor pressures with measured Raman shifts of C-H symmetric stretching band (υ1) in the vapor phase of sample fluids near room temperature. This equation eliminates discrepancies among the existing data sets and can be applied at any Raman laboratory. Raman shifts of C-H symmetric stretching band of methane in the vapor phase of CH4-H2O mixtures prepared in a high-pressure optical cell were also measured at temperatures between room temperature and 200 °C, and pressures up to 37 MPa. The results show that the CH4υ1 band position shifts to higher wavenumber as temperature increases. We also demonstrated that this Raman band shift is a simple function of methane vapor density, and, therefore, when combined with equation of state of methane, methane vapor pressures in the sample fluids at elevated temperatures can be calculated from measured Raman peak positions. This method can be applied to determine the pressure of CH4-bearing systems, such as methane-rich fluid inclusions from sedimentary basins or experimental fluids in hydrothermal diamond-anvil cell or other types of optical cell.  相似文献   

7.
The composition and evolution of a metallic planetary core is determined by the behavior with pressure of the eutectic and the liquidus on the Fe-rich side of the Fe-FeS eutectic. New experiments at 6 GPa presented here, along with existing experimental data, inform a thermodynamic model for this liquidus from 1 bar to at least 10 GPa. Fe-FeS has a eutectic that becomes more Fe-rich but remains constant in T up to 6 GPa. The 1 bar, 3 GPa, and 6 GPa liquidi all cross at a pivot point at 1640 ± 5 K and FeS37 ± 0.5. This liquid/crystalline metal equilibrium is T-x-fixed and pressure independent through 6 GPa. Models of the 1 bar through 10 GPa experimental liquidi show that with increasing P there is an increase in the T separation between the liquidus and the crest of the metastable two-liquid solvus. The solvus crest decreases in T with increasing P. The model accurately reproduces all the experimental liquidi from 1 bar to 10 GPa, as well as reproducing the 0-6 GPa pivot point. The 14 GPa experimental liquidus ( [Chen et al., 2008a] and Chen et al., 2008b) deviates sharply from the lower pressure trends indicating that the 0-10 GPa model no longer applies to this 14 GPa data.  相似文献   

8.
A thermodynamic model is presented to calculate methane solubility, liquid phase density and gas phase composition of the H2O-CH4 and H2O-CH4-NaCl systems from 273 to 523 K (possibly up to 573 K), from 1 to 2000 bar and from 0 to 6 mol kg−1 of NaCl with experimental accuracy. By a more strict theoretical approach and using updated experimental data, this model made substantial improvements over previous models: (1) the accuracy of methane solubility in pure water in the temperature range between 273 and 283 K is increased from about 10% to about 5%, but confirms the accuracy of the Duan model [Duan Z., Moller N., Weare J.H., 1992a. Prediction of methane solubilities in natural waters to high ionic strength from 0 to 250 °C and from 0 to 1600 bar. Geochim. Cosmochim. Acta56, 1451-1460] above 283 K up to 2000 bar; (2) the accuracy of methane solubility in the NaCl aqueous solutions is increased from >12% to about 6% on average from 273 K and 1 bar to 523 K and 2000 bar; (3) this model is able to calculate water content in the gas phase and liquid phase density, which cannot be calculated by previous models; and (4) it covers a wider range of temperature and pressure space. With a simple approach, this model is extended to predict CH4 solubility in other aqueous salt solutions containing Na+, K+, Mg2+, Ca2+, Cl− and , such as seawater and geothermal brines, with excellent accuracy. This model is also able to calculate homogenization pressure of fluid inclusions (CH4-H2O-NaCl) and CH4 solubility in water at gas-liquid-hydrate phase equilibrium. A computer code is developed for this model and can be downloaded from the website: www.geochem-model.org/programs.htm.  相似文献   

9.
The isopiestic method has been used to determine the osmotic coefficients of the binary solutions NaBr-H2O (from 0.745 to 5.953 mol kg−1) and KBr-H2O (from 0.741 to 5.683 mol kg−1) at the temperature t = 50 °C. Sodium chloride solutions have been used as isopiestic reference standards. The isopiestic results obtained have been combined with all other experimental thermodynamic quantities available in literature (osmotic coefficients, water activities, bromide mineral’s solubilities) to construct a chemical model that calculates solute and solvent activities and solid-liquid equilibria in the NaBr-H2O, KBr-H2O and Na-K-Br-H2O systems from dilute to high solution concentration within the 0-300 °C temperature range. The Harvie and Weare [Harvie C., and Weare J. (1980) The prediction of mineral solubilities in naturalwaters: the Na-K-Mg-Ca-Cl-SO4-H2O system from zero to high concentration at 25 °C. Geochim. Cosmochim. Acta44, 981-997] solubility modeling approach, incorporating their implementation of the concentration-dependent specific interaction equations of Pitzer [Pitzer K. (1973) Thermodynamics of electrolytes. I. Theoretical basis and general equations. J. Phys. Chem.77, 268-277] is employed. The model for binary systems is validated by comparing activity coefficient predictions with those given in literature, and not used in the parameterization process. Limitations of the mixed solutions model due to data insufficiencies are discussed. This model expands the variable temperature sodium-potassium model of Greenberg and Moller [Greenberg J., and Moller N. (1989) The prediction of mineral solubilities in natural waters: a chemical equilibrium model for the Na-K-Ca-Cl-SO4-H2O system to high concentration from 0 to 250 °C. Geochim. Cosmochim. Acta53, 2503-2518] by evaluating Br− pure electrolyte and mixing solution parameters and the chemical potentials of three bromide solid phases: NaBr-2H2O (cr), NaBr (cr) and KBr (cr).  相似文献   

10.
Aqueous fluids in sedimentary basins often contain dissolved methane, particularly in petroleum environments. PVTX (Pressure-Volume-Temperature-Composition) reconstructions performed using fluid inclusion data are largely based on the assumption that inclusions do not change from the time of trapping until the present. Many authors, however, consider that fluid inclusions can re-equilibrate, particularly in fragile minerals like calcite. In order to understand this re-equilibration phenomenon in the metamorphic domain, previous experiments have been performed under high PT conditions, but few have been performed at low to medium PT conditions such as those associated with sedimentary burial diagenesis, and no previous studies have examined CH4-bearing aqueous inclusions in calcite.An experimental study of the preservation/modification of CH4-rich synthetic fluid inclusions in calcite during isothermal decompression was conducted. An autoclave was used for accurate PTX control allowing equilibrium between liquid and vapour in the CH4-H2O system. PTX conditions were maintained at four stages of decreasing pressure, with each stage held for 7 days to simulate an isothermal pressure drop. In order of decreasing pressure, the pressure-temperature conditions monitored were 276 ± 10 bar at 180 ± 7 °C, 176 ± 10 bar at 180 ± 7 °C, 76 ± 10 bar at 180 ± 7 °C and 10 ± 3 bar at 180 ± 15 °C. At the end of the experiment, the calcite was recovered and analyzed by microthermometry and Raman microspectroscopy for PTX reconstruction. A careful procedure was adopted to limit re-equilibration of inclusions during analytical procedures. Four types of inclusion shapes and four types of strain patterns were differentiated. Classification of the petrographic strain patterns was carried out. These strain patterns were associated with inclusion stretching and/or leakage regarding CH4, Th and Ph compared to experimental conditions. Factors controlling the preservation or acquisition of strain patterns included the initial shape and size of the inclusion, and the pressure differential (ΔP) between the confining pressure (Pcf) and the internal pressure (Pi) within the inclusion. Most fluid inclusions seemed to be trapped during the first 7 days of the experiment, although few (4%) of these preserved the initial PT conditions of 276 ± 10 bar, whereas 8% preserved the second and third run of PT conditions. Overall, the majority of inclusions (88%) did not reflect accurately the PTX trapping conditions. A petrographic guide to the inclusions is presented here that allows strain identification for PVT reconstructions. Re-equilibration patterns and evidence for preferential methane leakage from aqueous inclusions in calcite are important findings revealed by this study, and may be useful for the reconstruction of post-trapping events in investigations of natural samples, and in other experiments using synthetic inclusions in calcite.  相似文献   

11.
The solubility of molybdenum (Mo) was determined at temperatures from 500 °C to 800 °C and 150 to 300 MPa in KCl-H2O and pure H2O solutions in cold-seal experiments. The solutions were trapped as synthetic fluid inclusions in quartz at experimental conditions, and analyzed by laser ablation inductively coupled plasma mass spectrometry (LA ICPMS).Mo solubilities of 1.6 wt% in the case of KCl-bearing aqueous solutions and up to 0.8 wt% in pure H2O were found. Mo solubility is temperature dependent, but not pressure dependent over the investigated range, and correlates positively with salinity (KCl concentration). Molar ratios of ∼1 for Mo/Cl and Mo/K are derived based on our data. In combination with results of synchrotron X-ray absorption spectroscopy of individual fluid inclusions, it is suggested that Mo-oxo-chloride complexes are present at high salinity (>20 wt% KCl) and ion pairs at moderate to low salinity (<11 wt% KCl) in KCl-H2O aqueous solutions. Similarly, in the pure H2O experiments molybdic acid is the dominant species in aqueous solution. The results of these hydrothermal Mo experiments fit with earlier studies conducted at lower temperatures and indicate that high Mo concentrations can be transported in aqueous solutions. Therefore, the Mo concentration in aqueous fluids seems not to be the limiting factor for ore formation, whereas precipitation processes and the availability of sulfur appear to be the main controlling factors in the formation of molybdenite (MoS2).  相似文献   

12.
The structure of H2O-saturated silicate melts and of silicate-saturated aqueous solutions, as well as that of supercritical silicate-rich aqueous liquids, has been characterized in-situ while the sample was at high temperature (to 800 °C) and pressure (up to 796 MPa). Structural information was obtained with confocal microRaman and with FTIR spectroscopy. Two Al-bearing glasses compositionally along the join Na2O•4SiO2-Na2O•4(NaAl)O2-H2O (5 and 10 mol% Al2O3, denoted NA5 and NA10) were used as starting materials. Fluids and melts were examined along pressure-temperature trajectories of isochores of H2O at nominal densities (from PVT properties of pure H2O) of 0.85 g/cm3 (NA10 experiments) and 0.86 g/cm3 (NA5 experiments) with the aluminosilicate + H2O sample contained in an externally-heated, Ir-gasketed hydrothermal diamond anvil cell.Molecular H2O (H2O°) and OH groups that form bonds with cations exist in all three phases. The OH/H2O° ratio is positively correlated with temperature and pressure (and, therefore, fugacity of H2O, fH2O) with (OH/H2O°)melt > (OH/H2O°)fluid at all pressures and temperatures. Structural units of Q3, Q2, Q1, and Q0 type occur together in fluids, in melts, and, when outside the two-phase melt + fluid boundary, in single-phase liquids. The abundance of Q0 and Q1 increases and Q2 and Q3 decrease with fH2O. Therefore, the NBO/T (nonbridging oxygen per tetrahedrally coordination cations), of melt is a positive function of fH2O. The NBO/T of silicate in coexisting aqueous fluid, although greater than in melt, is less sensitive to fH2O.The melt structural data are used to describe relationships between activity of H2O and melting phase relations of silicate systems at high pressure and temperature. The data were also combined with available partial molar configurational heat capacity of Qn-species in melts to illustrate how these quantities can be employed to estimate relationships between heat capacity of melts and their H2O content.  相似文献   

13.
The solubility and stability of synthetic grossular were determined at 800 °C and 10 kbar in NaCl-H2O solutions over a large range of salinity. The measurements were made by evaluating the weight losses of grossular, corundum, and wollastonite crystals equilibrated with fluid for up to one week in Pt capsules and a piston-cylinder apparatus. Grossular dissolves congruently over the entire salinity range and displays a large solubility increase of 0.0053 to 0.132 molal Ca3Al2Si3O12 with increasing NaCl mole fraction (XNaCl) from 0 to 0.4. There is thus a solubility enhancement 25 times the pure H2O value over the investigated range, indicating strong solute interaction with NaCl. The Ca3Al2Si3O12 mole fraction versus NaCl mole fraction curve has a broad plateau between XNaCl = 0.2 and 0.4, indicating that the solute products are hydrous; the enhancement effect of NaCl interaction is eventually overtaken by the destabilizing effect of lowering H2O activity. In this respect, the solubility behavior of grossular in NaCl solutions is similar to that of corundum and wollastonite. There is a substantial field of stability of grossular at 800 °C and 10 kbar in the system CaSiO3-Al2O3-H2O-NaCl. At high Al2O3/CaSiO3 bulk compositions the grossular + fluid field is limited by the appearance of corundum. Zoisite appears metastably with corundum in initially pure H2O, but disappears once grossular is nucleated. At XNaCl = 0.3, however, zoisite is stable with corundum and fluid; this is the only departure from the quaternary system encountered in this study. Corundum solubility is very high in solutions containing both NaCl and CaSiO3: Al2O3 molality increases from 0.0013 in initially pure H2O to near 0.15 at XNaCl = 0.4 in CaSiO3-saturated solutions, a >100-fold enhancement. In contrast, addition of Al2O3 to wollastonite-saturated NaCl solutions increases CaSiO3 molality by only 12%. This suggests that at high pH (quench pH is 11-12), the stability of solute Ca chloride and Na-Al ± Si complexes account for high Al2O3 solubility, and that Ca-Al ± Si complexes are minor. The high solubility and basic dissolution reaction of grossular suggest that Al may be a very mobile component in calcareous rocks in the deep crust and upper mantle when migrating saline solutions are present.  相似文献   

14.
The structure of silicate melts in the system Na2O·4SiO2 saturated with reduced C-O-H volatile components and of coexisting silicate-saturated C-O-H solutions has been determined in a hydrothermal diamond anvil cell (HDAC) by using confocal microRaman and FTIR spectroscopy as structural probes. The experiments were conducted in-situ with the melt and fluid at high temperature (up to 800 °C) and pressure (up to 1435 MPa). Redox conditions in the HDAC were controlled with the reaction, Mo + H2O = MoO2 + H2, which is slightly more reducing than the Fe + H2O = FeO + H2 buffer at 800 °C and less.The dominant species in the fluid are CH4 + H2O together with minor amounts of molecular H2 and an undersaturated hydrocarbon species. In coexisting melt, CH3 - groups linked to the silicate melt structure via Si-O-CH3 bonding may dominate and possibly coexists with molecular CH4. The abundance ratio of CH3 - groups in melts relative to CH4 in fluids increases from 0.01 to 0.07 between 500 and 800 °C. Carbon-bearing species in melts were not detected at temperatures and pressures below 400 °C and 730 MPa, respectively. A schematic solution mechanism is, Si-O-Si + CH4?Si-O-CH3+H-O-Si. This mechanism causes depolymerization of silicate melts. Solution of reduced (C-O-H) components will, therefore, affect melt properties in a manner resembling dissolved H2O.  相似文献   

15.
High-pressure liquids in the MgO-SiO2-H2O (MSH) system have been investigated at 11 and 13.5 GPa and between 1000 and 1350 °C. A bulk composition more magnesian than the tie-line forsterite-H2O was employed for the study. Rocking multi-anvil experiments were combined with a diamond trap set-up. After termination of the experiments, the liquid trapped in the diamond layer was analysed by laser ablation ICP-MS using the ‘freezing’ technique. At 11 GPa, liquids coexist with one or two of phase A, clinohumite, chondrodite, and forsterite. A marked discontinuity in the evolution of liquid compositions near 1100 °C is observed at 11 GPa. A step of ∼13 wt% H2O and 13 wt% MgO is interpreted to result from overstepping the fluid-saturated solidus reaction mass balanced to 1.00(18) phase A + 1.07(4) fluid = 0.63(15) chondrodite + 1.44(2) melt. At 13.5 GPa liquids coexist with one or two of hydrous wadsleyite, clinohumite, superhydrous B, phase B, and forsterite. The discontinuity in liquid composition is no longer present, indicating that the second critical endpoint of the solidus has been overstepped. Thus, hydrous melts in the Mg-rich part of the MSH system (molar bulk Mg/Si > 2) are chemically distinct from aqueous fluids at pressure up to 11 GPa. Convergence of fluid and melt compositions along the solidus resulting in a supercritical liquid occurs between 11 and 13.5 GPa, at which pressure the entire MSH system becomes supercritical.  相似文献   

16.
The electrical conductivities of aqueous solutions of Li2SO4 and K2SO4 have been measured at 523-673 K at 20-29 MPa in dilute solutions for molalities up to 2 × 10−2 mol kg−1. These conductivities have been fitted to the conductance equation of Turq, Blum, Bernard, and Kunz with a consensus mixing rule and mean spherical approximation activity coefficients. In the temperature interval 523-653 K, where the dielectric constant, ε, is greater than 14, the electrical conductance data can be fitted by a solution model which includes ion association to form , , and , where M is Li or K. The adjustable parameters of this model are the first and second dissociation constants of the M2SO4. For the 673 K and 300 kg m−3 state point where the Coulomb interactions are the strongest (dielectric constant, ε = 5), models with more extensive association give good fits to the data. In the case of the Li2SO4 model, including the multi-ion associate, , gave an extremely good fit to the conductance data.  相似文献   

17.
  相似文献   

18.
We perform first principles molecular dynamics simulations of Mg2SiO4 liquid and crystalline forsterite. On compression by a factor of two, we find that the Grüneisen parameter of the liquid increases linearly from 0.6 to 1.2. Comparison of liquid and forsterite equations of state reveals a temperature-dependent density crossover at pressures of ∼12-17 GPa. Along the melting curve, which we calculate by integration of the Clapeyron equation, the density crossover occurs within the forsterite stability field at P = 13 GPa and T = 2550 K. The melting curve obtained from the root mean-square atomic displacement in forsterite using the Lindemann law fails to match experimental or calculated melting curves. We attribute this failure to the liquid structure that differs significantly from that of forsterite, and which changes markedly upon compression, with increases in the degree of polymerization and coordination. The mean Si coordination increases from 4 in the uncompressed system to 6 upon twofold compression. The self-diffusion coefficients increase with temperature and decrease monotonically with pressure, and are well described by the Arrhenian relation. We compare our equation of state to the available highpressure shock wave data for forsterite and wadsleyite. Our theoretical liquid Hugoniot is consistent with partial melting along the forsterite Hugoniot at pressures 150-170 GPa, and complete melting at 170 GPa. The wadsleyite Hugoniot is likely sub-liquidus at the highest experimental pressure to date (200 GPa).  相似文献   

19.
Solubilities of corundum (Al2O3) and wollastonite (CaSiO3) were measured in H2O-NaCl solutions at 800 °C and 10 kbar and NaCl concentrations up to halite saturation by weight-loss methods. Additional data on quartz solubility at a single NaCl concentration were obtained as a supplement to previous work. Single crystals of synthetic corundum, natural wollastonite or natural quartz were equilibrated with H2O and NaCl at pressure (P) and temperature (T) in a piston-cylinder apparatus with NaCl pressure medium and graphite heater sleeves. The three minerals show fundamentally different dissolution behavior. Corundum solubility undergoes large enhancement with NaCl concentration, rising rapidly from Al2O3 molality (mAl2O3) of 0.0013(1) (1σ error) in pure H2O and then leveling off to a maximum of ∼0.015 at halite saturation (XNaCl ≈ 0.58, where X is mole fraction). Solubility enhancement relative to that in pure H2O, , passes through a maximum at XNaCl ≈ 0.15 and then declines towards halite saturation. Quenched fluids have neutral pH at 25 °C. Wollastonite has low solubility in pure H2O at this P and T(mCaSiO3=0.0167(6)). It undergoes great enhancement, with a maximum solubility relative to that in H2O at XNaCl ≈ 0.33, and solubility >0.5 molal at halite saturation. Solute silica is 2.5 times higher than at quartz saturation in the system H2O-NaCl-SiO2, and quenched fluids are very basic (pH 11). Quartz shows monotonically decreasing solubility from mSiO2=1.248 in pure H2O to 0.202 at halite saturation. Quenched fluids are pH neutral. A simple ideal-mixing model for quartz-saturated solutions that requires as input only the solubility and speciation of silica in pure H2O reproduces the data and indicates that hydrogen bonding of molecular H2O to dissolved silica species is thermodynamically negligible. The maxima in for corundum and wollastonite indicate that the solute products include hydrates and Na+ and/or Cl− species produced by molar ratios of reactant H2O to NaCl of 6:1 and 2:1, respectively. Our results imply that quite simple mechanisms may exist in the dissolution of common rock-forming minerals in saline fluids at high P and T and allow assessment of the interaction of simple, congruently soluble rock-forming minerals with brines associated with deep-crustal metamorphism.  相似文献   

20.
In order to quantify Al transfer in response to fluid-mineral equilibration under evolving metamorphic conditions, isobaric (0.7 GPa) experiments were conducted in the 350-550 °C range. Disequilibrium was induced (1) by holding initially pure water and natural minerals (kyanite + quartz ± muscovite enclosed in a perforated inner capsule) under isothermal conditions and (2) by stepwise temperature variations. In all experiments, secondary Al-bearing phases crystallized in the external tube of a “tube-in-tube” setup (SEM characterization); they are interpreted as witnesses of the evolution of the fluid composition (fluid reaction path). These reaction paths and the subsequent amount of secondary crystallizations were modeled using thermodynamic data from SUPCRT92 and estimates of both starting-mineral dissolution rates and elemental diffusion coefficients from the literature. A major result is that the amount of aluminum transferred to secondary phases is a thousand times larger than the calculated Al concentration in the fluid. Although the crystallization of Al-bearing phases was expected as a response to a temperature decrease, the stepwise temperature increase (20 °C/day) also led to aluminum transfer towards secondary phases. In the course of re-equilibration, the fluid first becomes saturated with respect to aluminosilicates and then reaches silica saturation, due to the low solubility of Al-minerals. Consequently, aluminosilicates partly recrystallize in response to a temperature increase. Crystallization of secondary Al-phases in the external tube implies that aqueous aluminum was efficiently transported from the inner capsule, even in the pure Al2O3-SiO2-H2O system. Therefore, mass balance calculations considering a constant Al reference frame, i.e., postulating Al immobility, should be regarded with caution.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号