Efficient chemical equilibrium calculations for geochemical speciation and reactive transport modelling

Abstract Chemical equilibrium calculations are essential for many environmental problems. It is also a fundamental tool for chemical kinetics and reactive transport modelling, since these applications may require hundreds to billions equilibrium calculations in a single simulation. Therefore, an equilibrium method for such critical applications must be very efficient, robust and accurate. In this work we demonstrate the potential effectiveness of a novel Gibbs energy minimisation algorithm for reactive transport simulations. The algorithm includes strategies to converge from poor initial guesses; capabilities to specify non-linear equilibrium constraints such as pH of an aqueous solution and activity or fugacity of a species; a rigorous phase stability test to determine the unstable phases; and a strategy to boost the convergence speed of the calculations to quadratic rates, requiring only few iterations to converge. We use this equilibrium method to solve geochemical problems relevant to carbon storage in saline aquifers, where aqueous, gaseous and minerals phases are present. The problems are formulated to mimic the ones found in kinetics and transport simulations, where a sequence of equilibrium calculations are performed, each one using the previous solution as the initial guess. The efficiency and convergence rates of the calculations are presented, which require an average of 1–2 iterations. These results indicate that critical applications such as chemical kinetics and reactive transport modelling can potentially benefit by using this multiphase equilibrium algorithm.

[1]  T. Tsuchiya,et al.  On the formulation and theory of the Newton interior-point method for nonlinear programming , 1996 .

[2]  H. Helgeson,et al.  Theoretical prediction of the thermodynamic behavior of aqueous electrolytes at high pressures and temperatures; I, Summary of the thermodynamic/electrostatic properties of the solvent , 1974 .

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

[4]  M. Fernanda P. Costa,et al.  On Minimizing Objective and KKT Error in a Filter Line Search Strategy for an Interior Point Method , 2011, ICCSA.

[5]  Jesús Carrera,et al.  CHEPROO: A Fortran 90 object-oriented module to solve chemical processes in Earth Science models , 2009, Comput. Geosci..

[6]  Dmitrii A. Kulik,et al.  The convex programming minimization of five thermodynamic potentials other than Gibbs energy in geochemical modeling , 2002 .

[7]  F. Morel,et al.  MINEQL-EIR. Chemical Equilibrium Composition of Aqueous Systems , 1986 .

[8]  H. Helgeson,et al.  Theoretical prediction of the thermodynamic behavior of aqueous electrolytes at high pressures and temperatures , 1974 .

[9]  Anthony V. Fiacco,et al.  Nonlinear programming;: Sequential unconstrained minimization techniques , 1968 .

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

[11]  R. Fletcher.  Practical Methods of Optimization , 1988 .

[12]  C. Kelley.  Iterative Methods for Linear and Nonlinear Equations , 1987 .

[13]  Thomas Wagner,et al.  GEM-Selektor geochemical modeling package: revised algorithm and GEMS3K numerical kernel for coupled simulation codes , 2012, Computational Geosciences.

[14]  W. Wagner,et al.  The IAPWS Formulation 1995 for the Thermodynamic Properties of Ordinary Water Substance for General and Scientific Use , 2002 .

[15]  Jesús Carrera,et al.  A mathematical formulation for reactive transport that eliminates mineral concentrations , 1998 .

[16]  J. Allison,et al.  MINTEQA2/PRODEFA2, a geochemical assessment model for environmental systems: Version 3. 0 user's manual , 1991 .

[17]  M. Fernanda P. Costa,et al.  Practical implementation of an interior point nonmonotone line search filter method , 2008, Int. J. Comput. Math..

[18]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[19]  B. Jones,et al.  WATEQ: A COMPUTER PROGRAM FOR CALCULATING CHEMICAL EQUILIBRIA OF NATURAL WATERS , 1973 .

[20]  Geoffrey C. Maitland,et al.  Phase equilibria of (CO2 + H2O + NaCl) and (CO2 + H2O + KCl): Measurements and modeling , 2013 .

[21]  Lorenz T. Biegler,et al.  Line Search Filter Methods for Nonlinear Programming: Motivation and Global Convergence , 2005, SIAM J. Optim..

[22]  C. Steefel,et al.  A coupled model for transport of multiple chemical species and kinetic precipitation/dissolution rea , 1994 .

[23]  J. Weare,et al.  A chemical equilibrium algorithm for highly non-ideal multiphase systems: Free energy minimization , 1987 .

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

[25]  Karsten Pruess,et al.  TOUGHREACT - A simulation program for non-isothermal multiphase reactive geochemical transport in variably saturated geologic media: Applications to geothermal injectivity and CO2 geological sequestration , 2006, Comput. Geosci..

[26]  Dmitrii A. Kulik,et al.  Dual-thermodynamic estimation of stoichiometry and stability of solid solution end members in aqueous-solid solution systems , 2006 .

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

[28]  C. Steefel,et al.  Reactive transport modeling: An essential tool and a new research approach for the Earth sciences , 2005 .

[29]  Allan M.M. Leal,et al.  A robust and efficient numerical method for multiphase equilibrium calculations: Application to CO2–brine–rock systems at high temperatures, pressures and salinities , 2013 .

[30]  Michael A. Saunders,et al.  SNOPT: An SQP Algorithm for Large-Scale Constrained Optimization , 2002, SIAM J. Optim..

[31]  Jorge Nocedal,et al.  An Interior Point Algorithm for Large-Scale Nonlinear Programming , 1999, SIAM J. Optim..

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

[33]  Yu. V. Shvarov,et al.  HCh: New potentialities for the thermodynamic simulation of geochemical systems offered by windows , 2008 .

[34]  Karsten Pruess,et al.  CO2-H2O mixtures in the geological sequestration of CO2. I. Assessment and calculation of mutual solubilities from 12 to 100°C and up to 600 bar , 2003 .

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

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

[37]  G. Eriksson,et al.  ChemSage—A computer program for the calculation of complex chemical equilibria , 1990 .

[38]  James J. Morgan,et al.  Numerical method for computing equilibriums in aqueous chemical systems , 1972 .

[39]  Edite Manuela da G. P. Fernandes,et al.  Assessing the potential of interior point barrier filter line search methods: nonmonotone versus monotone approach , 2011 .

[40]  Ionel Michael Navon,et al.  VARIATM—A FORTRAN program for objective analysis of pseudostress wind fields using large-scale conjugate-gradient minimization , 1991 .

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

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

[43]  Gunnar Eriksson,et al.  FactSage thermochemical software and databases - recent developments , 2009 .

[44]  John H. Weare,et al.  An equation of state for the CH4-CO2-H2O system: I. Pure systems from 0 to 1000°C and 0 to 8000 bar , 1992 .

[45]  Dominique Orban,et al.  Global minimization of the Gibbs energy of multicomponent systems Involving the presence of order/disorder phase transitions , 2013, American Journal of Science.

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

[47]  Stefan Ulbrich,et al.  A globally convergent primal-dual interior-point filter method for nonlinear programming , 2004, Math. Program..

[48]  C. Bethke.  Geochemical and Biogeochemical Reaction Modeling , 2007 .

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

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

[51]  C. Steefel,et al.  Approaches to modeling of reactive transport in porous media , 1996 .

[52]  C. Lemaréchal,et al.  The watchdog technique for forcing convergence in algorithms for constrained optimization , 1982 .

[53]  Karsten Pruess,et al.  TOUGHREACT User's Guide: A Simulation Program for Non-isothermal Multiphase Reactive geochemical Transport in Variable Saturated Geologic Media , 2004 .

[54]  C. Steefel,et al.  A new kinetic approach to modeling water-rock interaction: The role of nucleation, precursors, and Ostwald ripening , 1990 .

[55]  Olaf Kolditz,et al.  Modeling reactive transport in non-ideal aqueous–solid solution system , 2009 .

[56]  Gunnar Eriksson,et al.  FactSage thermochemical software and databases , 2002 .

[57]  R. Vanderbei.  LOQO:an interior point code for quadratic programming , 1999 .

[58]  Thomas H. Brown,et al.  The computation of chemical equilibrium in complex systems containing non-ideal solutions , 1987 .

[59]  Lorenz T. Biegler,et al.  On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming , 2006, Math. Program..

[60]  T. J. Wolery,et al.  EQ3/6, a software package for geochemical modeling of aqueous systems: Package overview and installation guide (Version 7.0) , 1992 .

[61]  Dmitrii A. Kulik,et al.  MODELLING CHEMICAL EQUILIBRIUM PARTITIONING WITH THE GEMS-PSI CODE , 2004 .

[62]  Hande Y. Benson,et al.  INTERIOR-POINT METHODS FOR NONCONVEX NONLINEAR PROGRAMMING: JAMMING AND COMPARATIVE NUMERICAL TESTING , 2000 .

[63]  James W. Ball,et al.  WATEQ4F - a personal computer Fortran translation of the geochemical model WATEQ2 with revised data base , 1987 .

[64]  Nicolas Spycher,et al.  Fugacity coefficients of H2, CO2, CH4, H2O and of H2O- CO2-CH4 mixtures: A virial equation treatment for moderate pressures and temperatures applicable to calculations of hydrothermal boiling , 1988 .

[65]  Sven Leyffer,et al.  Nonlinear programming without a penalty function , 2002, Math. Program..

[66]  Robert J. Vanderbei,et al.  LOQO User’s Manual – Version 4.05 , 2006 .

[67]  Zhenhao Duan,et al.  An improved model calculating CO2 solubility in pure water and aqueous NaCl solutions from 273 to 533 K and from 0 to 2000 bar , 2003 .

[68]  Zhenhao Duan,et al.  An improved model for the calculation of CO2 solubility in aqueous solutions containing Na+, K+, Ca2+, Mg2+, Cl−, and SO42− , 2006 .

[69]  David L. Parkhurst,et al.  Description of input and examples for PHREEQC version 3: a computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations , 2013 .

[70]  Lorenz T. Biegler,et al.  Line Search Filter Methods for Nonlinear Programming: Local Convergence , 2005, SIAM J. Optim..

[71]  S. H. Storey,et al.  The computation of chemical equilibria , 1970 .

[72]  D. L. Parkhurst,et al.  User's guide to PHREEQC (Version 2)-a computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations , 1999 .

[73]  William R. Smith,et al.  Chemical Reaction Equilibrium Analysis: Theory and Algorithms , 1982 .