Effects of complex formation in flowing fluids on the hydrothermal solubilities of minerals as a function of fluid pressure and temperature in the critical and supercritical regions of the system H2O

Abstract Consideration of Darcy's law and its analogs for open channel flow, together with the thermodynamics of hydrolysis reactions in hydrothermal systems indicates that either fluid pressure and/or isothermal mineral solubilities or both may decrease or increase in the direction of fluid flow, depending on the volume of reaction and the permeability or aperture, cross-sectional area, and angle of flow as a function of distance along the flow path. Recent progress in theoretical geochemistry has led to improved equations of state which can be used to calculate the standard partial molal thermodynamic properties of both charged and neutral inorganic and organic aqueous species at pressures and temperatures to 1000°C and 5 kb ( Tanger and Helgeson , 1988; Shock and Helgeson , 1988, 1990; Shock , et al., 1989, 1992; Sverjensky et al., 1992). Thermodynamic properties generated from these revised equations of state for the hydrolysis of minerals to form aqueous complexes at high pressures and temperatures indicate that the signs of the standard partial molal volume, enthalpy, and heat capacity of reaction depend primarily on the number of ligands in the complexes, as well as their charge. If polyligand complexes and/or certain neutral aqueous species appear on the right side of the reaction, the isobaric and isothermal partial derivatives of the logarithm of the equilibrium constant (logK) at PSAT may tend toward infinity and negative infinity, respectively, as fluid pressure and temperature increase in the liquid phase region and approach the critical point of H2O. This behavior results in positive values of (∂logK/∂T)p and negative values of (∂logK/∂P)T at supercritical pressures and temperatures. For example, thermodynamic calculations indicate that values of log K for reactions representing hydrothermal sulfide solubilities in the acid pH range where the predominant sulfide species is H2S(aq) decrease with increasing fluid pressure to an increasing degree with increasing temperature, which is consistent with experimental data reported by Hemley et al. (1986, 1992). In contrast, (∂logK/∂T)p and (∂logK/∂P)T in the supercritical region may be negative and positive, respectively, for reactions representing sulfide solubilities in hydrothermal solutions with higher pHs where HS− predominates over H2S, but only if the chloride concentration is low. The opposite may be the case in concentrated alkali chloride solutions, regardless of the pH. Similar calculations indicate that log K for the incongruent reaction of K-feldspar or other aluminosilicates with supercritical hydrothermal solutions to form quartz and A1(OH)4−increases monotonically with increasing fluid pressure at constant temperature. However, the log K values maximize with increasing temperature at all pressures to at least ~3 kb, which is not true of the solubility of quartz. In contrast, values of log K for analogous reactions written in terms of Al3+ or Al (OH)2+ minimize with increasing temperature at constant pressure.

[1]  E. Oelkers,et al.  Calculation of the thermodynamic properties of aqueous species at high pressures and temperatures. Effective electrostatic radii, dissociation constants and standard partial molal properties to 1000 °C and 5 kbar , 1992 .

[2]  G. Robinson,et al.  Hydrothermal ore-forming processes in the light of studies in rock-buffered systems; I, Iron-copper-zinc-lead sulfide solubility relations , 1992 .

[3]  D. Norton,et al.  Transport phenomena in hydrothermal systems; cooling plutons , 1977 .

[4]  K. K. Kelley,et al.  AN EQUATION FOR THE REPRESENTATION OF HIGH-TEMPERATURE HEAT CONTENT DATA1 , 1932 .

[5]  W. B. Kamb The thermodynamic theory of nonhydrostatically stressed solids , 1961 .

[6]  James W. Johnson,et al.  Theoretical prediction of hydrothermal conditions and chemical equilibria during skarn formation in porphyry copper systems , 1985 .

[7]  Everett L. Shock,et al.  Calculation of the thermodynamic and transport properties of aqueous species at high pressures and temperatures: Standard partial molal properties of organic species , 1990 .

[8]  H. Helgeson,et al.  Theoretical prediction of the thermodynamic behavior of aqueous electrolytes by high pressures and temperatures; IV, Calculation of activity coefficients, osmotic coefficients, and apparent molal and standard and relative partial molal properties to 600 degrees C and 5kb , 1981 .

[9]  D. Norton Metasomatism and permeability , 1988 .

[10]  H. Helgeson,et al.  UNIFIED DESCRIPTION OF INCONGRUENT REACTIONS AND MINERAL SOLUBILITIES AS A FUNCTION OF BULK COMPOSITION AND SOLUTION pH IN HYDROTHERMAL SYSTEMS , 1991 .

[11]  J. K. Costain,et al.  Hydroseismicity—A hypothesis for the role of water in the generation of intraplate seismicity , 1987 .

[12]  Stephen R. Brown,et al.  Fluid flow through rock joints: The effect of surface roughness , 1987 .

[13]  R. Berman,et al.  Heat capacity of minerals in the system Na2O-K2O-CaO-MgO-FeO-Fe2O3-Al2O3-Sio2-TiO2-H2O-CO2: representation, estimation, and high temperature extrapolation , 1985 .

[14]  E. Oelkers,et al.  Calculation of activity coefficients and degrees of formation of neutral ion pairs in supercritical electrolyte solutions , 1991 .

[15]  Roger Powell,et al.  An enlarged and updated internally consistent thermodynamic dataset with uncertainties and correlations: the system K2O–Na2O–CaO–MgO–MnO–FeO–Fe2O3–Al2O3–TiO2–SiO2–C–H2–O2 , 1990 .

[16]  J. Bear Dynamics of Fluids in Porous Media , 1975 .

[17]  H. Helgeson Prediction of the thermodynamic properties of electrolytes at high pressures and temperatures , 1981 .

[18]  W. D'angelo,et al.  Effect of pressure on ore mineral solubilities under hydrothermal conditions , 1986 .

[19]  Peter C. Lichtner,et al.  Continuum model for simultaneous chemical reactions and mass transport in hydrothermal systems , 1985 .

[20]  J. Gibbs On the equilibrium of heterogeneous substances , 1878, American Journal of Science and Arts.

[21]  R. Wood,et al.  A relation between the critical properties of aqueous salt solutions and the heat capacity of the solutions near the critical point using a single-fluid corresponding-states theory☆ , 1982 .

[22]  J. Gates,et al.  Experimental evidence for the remarkable behavior of the partial molar heat capacity at infinite dilution of aqueous electrolytes at the critical point , 1982 .

[23]  H. Green On the thermodynamics of non-hydrostatically stressed solids , 1980 .

[24]  E. Oelkers,et al.  Calculation of dissociation constants and the relative stabilities of polynuclear clusters of 1:1 electrolytes in hydrothermal solutions at supercritical pressures and temperatures , 1993 .

[25]  H. Helgeson,et al.  Theoretical prediction of thermodynamic properties of aqueous electrolytes at high pressures and temperatures. III. Equation of state for aqueous species at infinite dilution , 1976 .

[26]  Lester Haar Nbs/Nrc Steam Tables , 1984 .

[27]  K. Pitzer Dielectric constant of water at very high temperature and pressure. , 1983, Proceedings of the National Academy of Sciences of the United States of America.

[28]  Y. A. Kozlovsky KOLA SUPER-DEEP: INTERIM RESULTS AND PROSPECTS , 1982 .

[29]  P. Lichtner The quasi-stationary state approximation to coupled mass transport and fluid-rock interaction in a porous medium , 1988 .

[30]  J. R. Fisher,et al.  Simultaneous evaluation and correlation of thermodynamic data , 1976 .

[31]  E. C. Beutner Slaty cleavage and related strain in Martinsburg Slate, Delaware Water Gap, New Jersey , 1978 .

[32]  E. Oelkers,et al.  Calculation of the thermodynamic and transport properties of aqueous species at high pressures and temperatures: dissociation constants for supercritical alkali metal halides at temperatures from 400 to 800.degree.C and pressures from 500 to 4000 bar , 1988 .

[33]  H. Helgeson Organic/inorganic reactions in metamorphic processes , 1991 .

[34]  Everett L. Shock,et al.  Calculation of the thermodynamic and transport properties of aqueous species at high pressures and temperatures: Correlation algorithms for ionic species and equation of state predictions to 5 kb and 1000°C , 1988 .

[35]  P. A. Witherspoon,et al.  Hydromechanical behavior of a deformable rock fracture subject to normal stress , 1981 .

[36]  H. Taylor,et al.  Stable isotopic evidence for large-scale seawater infiltration in a regional metamorphic terrane; the Trois Seigneurs Massif, Pyrenees, France , 1985 .

[37]  S. Saxena,et al.  An equation for the heat capacity of solids , 1987 .

[38]  D. Norton,et al.  Transport phenomena in hydrothermal systems: the nature of porosity , 1977 .

[39]  J. Hunt,et al.  Hydrothermal ore-forming processes in the light of studies in rock-buffered systems; II, Some general geologic applications , 1992 .

[40]  Chin-Fu Tsang,et al.  Channel model of flow through fractured media , 1987 .

[41]  E. Oelkers,et al.  SUPCRT92: a software package for calculating the standard molal thermodynamic properties of minerals, gases, aqueous species, and reactions from 1 to 5000 bar and 0 to 1000 ° C , 1992 .

[42]  E. Oelkers,et al.  Triple-ion anions and polynuclear complexing in supercritical electrolyte solutions , 1990 .

[43]  H. Taylor,et al.  Quantitative Simulation of the Hydrothermal Systems of Crystallizing Magmas on the Basis of Transport Theory and Oxygen Isotope Data: An analysis of the Skaergaard Intrusion , 1979 .

[44]  T. L. Tour,et al.  Fluid participation in deep fault zones: Evidence from geological, geochemical, and 18O/16O relations , 1984 .

[45]  Y. Tsang,et al.  The Effect of Tortuosity on Fluid Flow Through a Single Fracture , 1984 .

[46]  P. Lichtner,et al.  Fluid flow and mineral reactions at high temperatures and pressures , 1987, Journal of the Geological Society.

[47]  James W. Johnson,et al.  Critical phenomena in hydrothermal systems; state, thermodynamic, electrostatic, and transport properties of H 2 O in the critical region , 1991 .

[48]  H. Helgeson Errata II; Thermodynamics of minerals, reactions, and aqueous solutions at high pressures and temperatures , 1985 .

[49]  J. Sengers,et al.  An Improved Representative Equation for the Dynamic Viscosity of Water Substance , 1980 .

[50]  H. Helgeson Thermodynamics of minerals, reactions, and aqueous solutions at high pressures and temperatures; errata , 1982 .

[51]  H. Helgeson,et al.  Calculation of the thermodynamic properties of aqueous silica and the solubility of quartz and its polymorphs at high pressures and temperatures , 1977 .

[52]  W. Brace Permeability of crystalline rocks: New in situ measurements , 1984 .

[53]  E. Shock,et al.  Erratum to Geochim. Cosmochim.: E. L. Shock and H. C. Helgeson: Cosmochimica Acta 52, pp. 2009–2036 , 1989 .

[54]  H. Helgeson,et al.  Calculation of the Thermodynamic and Transport Properties of Aqueous Species at High Pressures and Temperatures; Revised Equations of State for the Standard Partial Molal Properties of Ions and Electrolytes , 1988, American Journal of Science.