K2GWB: Utility for generating thermodynamic data files for The Geochemist's Workbench® at 0-1000 °C and 1-5000 bar from UT2K and the UNITHERM database

The Geochemists Workbench^(R) (GWB) is a versatile geochemical modelling software package capable of reaction path modelling, equilibrium fluid speciation and activity diagramming, although it is restricted by the availability of thermodynamic (logK) data files with a temperature range of 0-300^oC and pressures on the critical curve for water. The K2GWB (logK to GWB) utility is used to produce a usable Geochemists Workbench data file from a generic logK and system file generated by the UT2K utility complementing the UNITHERM thermodynamic database system. The GWB data files can be generated for pressure-temperature conditions within the bounds of the modified Helgeson-Kirkham-Flowers model (0-1000^oC, 1-5000bar). Data can be generated with any combination of pressure-temperature logK grids required, allowing the Geochemists Workbench to be used to model a wide range of geological scenarios under polythermal and polybaric conditions, as well as the generation of activity diagrams over a wide range of pressures and temperatures. A comparative HCh and Geochemists Workbench model using the same thermodynamic dataset was run to benchmark test the validity of the new data files. Modelling the cooling of 2molal NaCl-H"2O fluid in equilibrium with a ganodiorite-type rock from 600 to 300^oC gave close to identical results between the two codes, validating the thermodynamic data conversion process described here. An investigation of the speciation of Au across a temperature-pressure gradient from 550^oC and 1500bar to 150^oC and saturated vapor curve pressures has been used to illustrate the usefulness of coupling thermodynamic data file generation at wide temperatures and pressures from a central database and the functionality of The Geochemists Workbench.

[1]  A. Stefánsson,et al.  Gold(I) complexing in aqueous sulphide solutions to 500°C at 500 bar , 2004 .

[2]  H. Helgeson,et al.  Thermodynamics of hydrothermal systems at elevated temperatures and pressures , 1969 .

[3]  L. Pankratz,et al.  Thermodynamic properties of elements and oxides , 1982 .

[4]  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 .

[5]  M. H. Lietzke,et al.  THE CALCULATION OF ACTIVITY COEFFICIENTS FROM OSMOTIC COEFFICIENT DATA1 , 1962 .

[6]  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.

[7]  R. Kretz.  Symbols for rock-forming minerals , 1983 .

[8]  M. Reed,et al.  Calculation of multicomponent chemical equilibria and reaction processes in systems involving minerals, gases and an aqueous phase , 1982 .

[9]  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 .

[10]  Y. Shvarov.  Algorithmization of the numeric equilibrium modeling of dynamic geochemical processes , 1999 .

[11]  E. Shock,et al.  Inorganic species in geologic fluids: correlations among standard molal thermodynamic properties of aqueous ions and hydroxide complexes. , 1997, Geochimica et cosmochimica acta.

[12]  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 .

[13]  E. Bastrakov.  Gold metallogenesis at the Lake Cowal Prospect, NSW , 2000 .

[14]  T. J. Wolery,et al.  EQ3NR, a computer program for geochemical aqueous speciation-solubility calculations: Theoretical manual, user`s guide, and related documentation (Version 7.0); Part 3 , 1992 .