Nancy Møller

The prediction of mineral solubilities in natural waters: The Na-K-Mg-Ca-H-Cl-SO4-OH-HCO3-CO3-CO2-H2O system to high ionic strengths at 25°C

The mineral solubility model of Harvie and Weare (1980) is extended to the eight component system, Na-K-Mg-Ca-H-Cl-SO4-OH-HCO3-CO3-CO2-H2O at 25°C to high concentrations. The model is based on the semi-empirical equations of Pitzer (1973) and co-workers for the thermodynamics of aqueous electrolyte solutions. The model is parameterized using many of the available isopiestic, electromotive force, and solubility data available for many of the subsystems. The predictive abilities of the model are demonstrated by comparison to experimental data in systems more complex than those used in parameterization. The essential features of a chemical model for aqueous electrolyte solutions and the relationship between pH and the equilibrium properties of a solution are discussed.

An equation of state for the CH4-CO2-H2O system: I. Pure systems from 0 to 1000°C and 0 to 8000 bar

An equation of state (EOS) for the CH4-CO2-H2O system covering a wide T-P range has been developed. In this article the new EOS is presented and applied to the pure endmembers. The equation is similar to that of Lee and Kesler (1975) and contains fifteen parameters. It is used with a mixing rule in the following article to provide a thermodynamic model for the mixed system. Though the parameters are evaluated from the PVT data in the temperature range from 0 to 450°C for CH4, from 0 to 1000°C for CO2 and H2O, and for pressures from 0 to 3500 bar, comparison of this EOS with a large amount of experimental data in the pure systems indicates that predictions for temperatures and pressures from 0 to 1000°C and 0 to 8000 bar (or slightly above) are very nearly within experimental uncertainty. The EOS can describe both the gas and the liquid phases of the endmember systems with similar accuracy. Fugacity coefficients are derived and compiled. In this paper mixing is considered using ideal mixing based on the endmember fugacities (Amagats rule). It is shown that such an approach leads to quite accurate predictions for high temperatures and low pressures.

The prediction of methane solubility in natural waters to high ionic strength from 0 to 250°C and from 0 to 1600 bar

Abstract A model for the solubility of methane in brines (0–6 m) for temperatures from 0 to 250°C and for pressures from 0 to 1600 bar (or slightly above) is presented. The model is based on Pitzer phenomenology for the liquid phase and a highly accurate equation of state recently developed for the vapor phase. Comparison of model predictions with experimental data indicates that they are within experimental uncertainty. Most experimental data sets are consistent within errors of about 7%. Although the parameters were evaluated from binary and ternary data, the model accurately predicts methane solubility in much more complicated systems like seawater and Salton geothermal brines. Application to fluid inclusion analysis is discussed. Minimum trapping pressures are calculated given the composition and homogenization temperature.

The prediction of mineral solubilities in natural waters: A chemical equilibrium model for the Na-Ca-Cl-SO4-H2O system, to high temperature and concentration

This paper describes a chemical equilibrium model for the Na-Ca-Cl-SO4-H2O system which calculates solubilities from 25°C to 250°C and from zero to high concentration (I ~ 18. m) within experimental uncertainty. The concentration and temperature dependence of the model were established by fitting available activity (solubility, osmotic and emf) data. A single ion complex, CaSO04, which increases in strength with temperature, is included explicitly in the model. n nThe validation of model accuracy by comparison to laboratory and field solubility data is included. Applications of the model are also given. Phase diagrams constructed for the Na-Ca-Cl-SO4-H2O system and predicted solubilities of anhydrite and hemihydrate in concentrated seawater at high temperature are in very good agreement with the data. Calculations of the temperature of gypsum-anhydrite coexistence as a function of water activity are compared to reported values, and are used to estimate the composition-temperature relation for gypsum-anhydrite transition in a natural brine evaporation. A preliminary model for barite solubility in sodium chloride solutions at high temperature (100°C to 250°C), based on this parameterization of the CaSO4-NaCl-H2O system, gives good agreement with the data.

Equation of state for the NaClH2OCO2 system: prediction of phase equilibria and volumetric properties

An equation of state (EOS) has been developed for the NaClue5f8H2Oue5f8CO2 system which consistently predicts various properties including PVTX, immiscibility or phase equilibria, solubilities, and activities with an accuracy close to that of experimental data from 300 to about 1000°C and 0–6000 bar with NaCl concentrations to about 30 wt% of NaCl (relative to NaCl + H2O) or to about 50 wt% with less accuracy. The EOS predicts that excess volumes can be over 30% of the total volume under some T-P conditions. Adding NaCl to the H2Oue5f8CO2 system dramatically increases the T-P range of immiscibility. The immiscibility field is minimal around 400–500°C. Above or below this temperature, it expands for a constant pressure. A moderately saline brine can evolve into a very saline brine by phase separation at high temperatures. The presence of NaCl can substantially decrease the activity of H2O and increase that of CO2, thus affecting decarbonation and dehydration reactions. Compared to the EOS of Bowers and Helgeson (1983a), the EOS of this study is more reliable in the calculation of volumetric properties particularly in the low pressure range. In addition, BH EOS cannot predict phase equilibria.

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☆

A chemical equilibrium model is described which is used to calculate solubilities within experimental uncertainties in the Na-K-Ca-Cl-SO4-H2O system from zero to high ionic strength and from 0 to 250°C. This model is an expansion of the variable temperature Na-Ca-Cl-SO4-H2O model of Moller (1988). It is parameterized by fitting available osmotic and solubility data in all common ion systems involving the potassium ion: Na-K-Cl-H2O, Na-K-SO4-H2O, K-Cl-SO4-H2O, K-Ca-Cl-H2O, and K-Ca-SO4-H2O. Limitations of the model due to data insufficiencies are discussed and comparisons of model calculations with the available data are given. Model predictions for solubility in the complex reciprocal systems, Na-K-Cl-SO4-H2O and K-Ca-Cl-SO4-H2O, are compared with experiment. Data for the two systems were available only in the temperature ranges 0–100°C and 25–55°C, respectively. The phase diagram predicted for the halite-saturated Na-K-Ca-Cl-SO4-H2O quinary system at 100°C is also presented. This model will be used to extend the Harvie et al. (1984) model for the Na-K-Ca-Mg-Cl-SO4-CO2-H2O seawater system to high temperature.

An equation of state for the CH4-CO2-H2O system: II. Mixtures from 50 to 1000°C and 0 to 1000 bar

An equation of state (EOS) for mixtures in the CH4-CO2-H2O system has been developed. The model is based on the highly accurate endmember EOS presented in the previous article and on an empirical mixing rule. The mixing rule is based on an analogy with high order contributions to the virial expansion for mixtures. Comparison with experimental data indicates that the mixed system EOS can predict both phase equilibria and volumetric properties for the binaries with accuracy close to that of the experimental data for a temperature range from 50 to 1000°C and a pressure range from 0 to 1000 bar (or to 3000 bar with less accuracy). n nFor temperatures below the critical point of water, there is very little PVTX (density) data. However, even for temperatures for which sufficient data exists we found that parameterization from PVTX data alone did not lead to a free energy that would accurately predict liquid-vapor equilibria. On the other hand, using this data alone we obtain a free energy that predicts both liquid-vapor equilibria and the PVTX properties of the binaries with roughly experimental accuracy. n nThere are very few data inside the ternary. However, the mixing rule contains third order parameters, which require evaluation from ternary mixtures. For the single temperature for which we have data, the adjustment of one parameter gives good prediction of phase equilibrium in the ternary. The resulting EOS predicts that the presence of a small amount of CO2 can significantly affect the solubility of CH4. The application of this EOS to the study of fluid inclusions is discussed. The presence of CH4 in CO2rich fluid inclusions can significantly affect the predicted trapping pressure.

The prediction of mineral solubilities in natural waters: A chemical equilibrium model for the NaKCaMgClSO4H2O system at temperatures below 25°C☆

Abstract A low temperature thermochemical model for the system Naue5f8Kue5f8Caue5f8Mgue5f8Clue5f8SO4ue5f8H2O is presented. Aqueous species and standard chemical potentials of solid-solution reactions are modeled from published data for binary and ternary solutions. The temperature range below 25°C (to near −60°C) is emphasized, although the model parameters are fitted to merge smoothly with those of higher temperature models at temperatures between 25 and 100°C. Binary and ternary specific ion interaction terms vary independently with temperature and are modeled using freezing point depression and mineral solubility measurements. The standard chemical potential of the ice-water reaction is fitted independent of the model (from vapor pressure and free energy data). Remaining standard chemical potentials of solidsolution reactions are fitted along with the specific ion interaction terms. Model predictions are tested against published data for minerals formed and brine compositions obtained by chilling seawater to the eutectic (about −54°C). The model predicts the sequence of solid phases observed to precipitate from chilled seawater (mice-mirabilite-hydrohalite-sylvite-MgCl2 · 12H2Oue5f8antarcticite). For all but mirabilite model temperatures are within the uncertainty of the measured temperature. The compositions of brines predicted by the model also closely follow the observed compositions. The model allows accurate predictions of the freezing points of simple and complex solutions in the system. Low temperature phase equilibria and mineral solubilities may also be predicted. The model may be used to determine the composition of brines in fluid inclusions in the multicomponent system based on low temperature phase equilibria.

A general equation of state for supercritical fluid mixtures and molecular dynamics simulation of mixture PVTX properties

Abstract A general Equation of State (EOS), which we previously developed for pure nonpolar systems, is extended to polar systems such as water and to mixtures in this study. This EOS contains only two parameters for each pure component and two additional parameters for each binary mixture (no higher order parameters are needed for more complicated mixture systems). The two mixing parameters can be eliminated for nonaqueous mixtures with a slight loss of accuracy in both total mole volume and in excess volume (or nonideal mixing). Comparison with a large amount of experimental PVTX data in pure systems (including H2O) and in the mixtures, H2O-CO2, CO2-N2, CH4-CO2, and N2-CO2-CH4 results in an average error of 1.6% in density. Comparison with commonly used EOS for supercritical fluids shows that the EOS of this study covers far more T-P-X space with higher accuracy. We believe that it is accurate from supercritical temperature to 2000 K and from 0 to 25,000 bar or higher with an average error in density of less than 2% for both pure members and mixtures in the system H2O-CO2-CH4-N2-CO-H2-O2-H2S-Ar and possibly with additional gases. Comparison with the published simulated data suggests that this EOS is approximately correct up to 300,000 bar and 2800 K. We also simulated the PVTX properties of a number of supercritical fluid mixtures using molecular dynamics (MD) simulation. These results and those of other authors are well predicted by the EOS of this study.

Molecular dynamics simulation of PVT properties of geological fluids and a general equation of state of nonpolar and weakly polar gases up to 2000 K and 20,000 bar

The PVT properties of CH[sub 4] from 30-360 cm[sup 3] mol and roughly from 273-2,000 K and from 100-20,000 bar have been simulated by molecular dynamics using Lennard-Jones potentials. The simulated results compare with data within 1.5% in volume. Using these simulated values and experimental PVT data, an equation of state (EOS) was developed. Because of the choice of potential, a simple scaling generalizes the EOS to predict the supercritical PVT properties of CO[sub 2], N[sub 2], CO, H[sub 2], O[sub 2], and Cl[sub 2] within an average error of about 1.5%.