The convex programming minimization of five thermodynamic potentials other than Gibbs energy in geochemical modeling

An approach is suggested for computation of complete or metastable equilibrium states in the multicomponent, multiphase, and multiaggregate systems with parameters of state other than fixed temperature T and pressure P. It consists in computation of a series of isobaric-isothermal equilibrium states by Gibbs energy minimization. Thus, the general problem of minimization of thermodynamic potentials other than Gibbs energy G(T,P) is reduced to a sequence of simpler problems of two- or three-criteria parametric minimization of either of five thermodynamic potentials [isochoric-isothermal or Helmholtz energy A(T,V); isochoric-isentropic or internal energy at isochoric conditions U(S,V); isobaric-isentropic or enthalpy H(P,S); negative entropy at isobaric conditions and fixed enthalpy -S(P,H); and negative entropy at isochoric conditions and fixed internal energy -S(V,U)]. The computational scheme is simple, efficient, and reliable, because in any case it uses the IPM (Interior Points Method) algorithm for searching minimum of the non-linear total Gibbs energy function whereby one- or two-sided restrictions may be imposed on some or all dependent components in the system (Karpov and others, 1997). The new approach provides a rigorous basis for the formulation and numerical solution of the problems of chemical mass transfer under factors of state other than isobaric-isothermal. The developed algorithms and their implementation in Selektor-C program package can be used in a wide variety of scientific and practical applications. The numerical tests and examples illustrate the new technique on a model of irreversible oil-water interaction in hydrocarbon reservoirs (system Ca-Mg-C-O-H).

[1]  E. Shock,et al.  Organic acids in hydrothermal solutions: standard molal thermodynamic properties of carboxylic acids and estimates of dissociation constants at high temperatures and pressures. , 1995, American journal of science.

[2]  J. Kiefer.  Optimum Experimental Designs , 1959 .

[3]  J. Fisher.  Distribution and occurrence of aliphatic acid anions in deep subsurface waters , 1987 .

[4]  T. Barth.  Organic acids and inorganic ions in waters from petroleum reservoirs, Norwegian continental shelf: a multivariate statistical analysis and comparison with American reservoir formation waters , 1991 .

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

[6]  N. Hubbard,et al.  Short-chain aliphatic acid anions in deep subsurface brines: A review of their origin, occurrence, properties, and importance and new data on their distribution and geochemical implications in the Palo Duro Basin, Texas , 1987 .

[7]  J. Prausnitz,et al.  Thermodynamic properties of supercritical fluids and their mixtures at very high pressures , 1973 .

[8]  Craig M. Bethke,et al.  Geochemical Reaction Modeling: Concepts and Applications , 1996 .

[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]  Everett L. Shock,et al.  Prediction of the thermodynamic properties of aqueous metal complexes to 1000°C and 5 kb , 1997 .

[11]  E. Shock,et al.  Petroleum, oil field waters, and authigenic mineral assemblages Are they in metastable equilibrium in hydrocarbon reservoirs☆ , 1993 .

[12]  R. Reid,et al.  The Properties of Gases and Liquids , 1977 .

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

[14]  M. S. Bazaraa,et al.  Nonlinear Programming , 1979 .

[15]  Everett L. Shock,et al.  Metal-organic complexes in geochemical processes: Calculation of standard partial molal thermodynamic properties of aqueous acetate complexes at high pressures and temperatures , 1993 .

[16]  Y. Kharaka,et al.  Thermal decarboxylation of acetic acid: Implications for origin of natural gas , 1983 .

[17]  J. B. Rapp,et al.  Short chain aliphatic acid anions in oil field waters and their contribution to the measured alkalinity , 1975 .

[18]  D. S. Korzhinskiĭ.  Physicochemical basis of the analysis of the paragenesis of minerals , 1959 .

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

[20]  Y. Kharaka,et al.  Aliphatic Acid Anions in Oil-Field Waters--Implications for Origin of Natural Gas , 1978 .

[21]  Mark S. Ghiorso,et al.  Chemical mass transfer in magmatic processes IV. A revised and internally consistent thermodynamic model for the interpolation and extrapolation of liquid-solid equilibria in magmatic systems at elevated temperatures and pressures , 1995 .

[22]  S. E. Drummond,et al.  Thermal decarboxylation of acetate. Part I. The kinetics and mechanism of reaction in aqueous solution , 1986 .

[23]  Sanford Gordon,et al.  CALCULATION OF COMPLEX CHEMICAL EQUILIBRIA , 1968 .

[24]  Warren D. Seider,et al.  Computation of Phase and Chemical Equilibrium: A Review , 1980 .

[25]  Byung-Ik Lee,et al.  A generalized thermodynamic correlation based on three‐parameter corresponding states , 1975 .

[26]  Konstantin V. Chudnenko,et al.  Modeling chemical mass transfer in geochemical processes; thermodynamic relations, conditions of equilibria and numerical algorithms , 1997 .

[27]  D. Crerar,et al.  Organic complexes in hydrothermal systems , 1988 .

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

[29]  Mokhtar S. Bazaraa,et al.  Nonlinear Programming: Theory and Algorithms , 1993 .

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

[31]  M. Ghiorso.  Modeling magmatic systems; thermodynamic relations , 1987 .

[32]  J. Hanor,et al.  Distribution of dissolved volatile fatty acids in some Louisiana oil field brines , 1986 .

[33]  D. R. Waldbaum.  Temperature Changes associated with Adiabatic Decompression in Geological Processes , 1971, Nature.