A critical evaluation of recent developments in hydrogeochemical transport models of reactive multichemical components

Water Resources Research - Tập 25 Số 1 - Trang 93-108 - 1989
Gour‐Tsyh Yeh1, Vishal Tripathi
1Graduate Institute of Applied Geology

Tóm tắt

Many hydrogeochemical models have appeared in recent years for simulating subsurface solute transport. The hydrological transport of solutes can be described by a set of linear partial differential equations, and the chemical equilibria are described by a set of nonlinear algebraic equations. Three approaches are currently used to formulate the problem: (1) the mixed differential and algebraic equation (DAE) approach, (2) the direct substitution approach (DSA), and (3) the sequential iteration approach (SIA). An extremely important consideration in any approach is the choice of primary dependent variables (PDVs). Six types of PDVs have been employed in the existing models: (1) concentrations of all species, (2) concentrations of all component species and precipitated species, (3) total analytical concentrations of aqueous components, (4) total dissolved concentrations of aqueous components, (5) concentrations of aqueous component species, and (6) hybrid concentrations. Because of many possible combinations of PDVs and approaches, many hydrogeochemical transport models for multicomponent systems have been developed. This paper critically evaluates and discusses these models. The discussion and evaluation are conducted in terms of (1) how severe can the constraints be that a model imposes on computer resources, (2) which equilibrium geochemical processes can a model include, and (3) how easily can a model be modified to deal with mixed kinetic and equilibrium reactions. The use of SIA models leads to the fewest constraints on computer resources in terms of central processing unit (CPU) memory and CPU time; both DAE and DSA models require excessive CPU memory and CPU time for realistic two‐ and three‐dimensional problems. Only those models that use the first three types of PDVs can treat the full complement of equilibrium reactions simultaneously. DAE and SIA models can be modified with reasonable ease to handle mixed chemical kinetics and equilibria. DSA models require strenuous efforts to modify for treating mixed chemical kinetics and equilibria. Therefore SIA models using the third type of PDVs are recommended for their practicality and flexibility. DSA and DAE models should remain research tools for one‐dimensional investigations.

Từ khóa


Tài liệu tham khảo

Bjerrum N., 1926, Ionic association, I, Influence of ionic association on the activity of ions at moderate degrees of association, Kgl. Danske Videnskab. Selskb. Math‐fys, Medd., 7, 1

10.1002/aic.690320505

10.1029/WR021i008p01095

Duff I. S. MA28—A set of FORTRAN subroutines for sparse unsymmetric linear equationsRep. R.8730Atomic Energy Res. Estab. Her Majesteys Stn. London 1977.

Freundlich H., 1926, Colloid and Capillary Chemistry

10.1021/cr60056a002

10.1111/j.1745-6584.1979.tb03319.x

Hageman L. A., 1981, Applied Iterative Methods

10.1016/0016-7037(68)90100-2

Helgeson H. C., 1979, Geochemistry of Hydrothermal Ore Deposits, 568

10.1145/1218052.1218054

10.1029/WR021i003p00346

10.1029/WR018i004p01089

10.1016/0309-1708(84)90040-X

10.1016/0022-1694(85)90092-7

Kirkner D. J. A. A.Jennings T. L.T L Theis Multicomponent mass face transport modeling for energy solid wastesFinal Rep. COO‐10253‐5Dep. of Civ. Eng. Univ. of Notre Dame Notre Dame Ind. 1987.

Lapidus L., 1982, Numerical Solution of Partial Differential Equations in Science and Engineering

10.1016/0022-1694(87)90174-0

10.1016/0016-7037(85)90172-3

Miller C. W. CHEMTRN USER's MANUAL LBL‐16152 Lawrence Berkeley Laboratory University of California Berkeley 1983.

10.1029/WR019i002p00381

10.1021/es60060a006

10.1029/WR022i013p01820

Parkhurst D. L., 1980, PHREEQE—A Computer Program for Geochemical Calculations, U. S. Geol. Surv. Water Resour. Invest., 80

Pauling L., 1956, General Chemistry

Petzold L. R. Differential/algebraic equations are not ODE'sRep. SAND81‐8668Sandia Nat. Lab. Albuquerque N.M. 1981.

10.1016/0016-7037(82)90155-7

10.1021/j100679a023

10.1029/WR019i005p01231

10.1029/WR009i005p01332

10.1029/WR019i002p00493

Spiegel M. R., 1963, Theory and Problems of Advanced Calculus

Stumm W., 1981, Aquatic Chemistry An Introduction Emphasizing Chemical Equilibria in Natural Waters

Theis T. L. D. J.Kirkner A. A.Jennings Multi‐solute subsurface transport modeling for energy solid wastes Tech. Progr. RepC00‐10253‐3Dep. of Civ. Eng. Univ. of Notre Dame Notre Dame Ind. 1982.

10.1016/0039-9140(86)80243-0

Truesdell A. H., 1974, WATEQ, A Computer Program for Calculating Chemical Equilibria of Natural Waters, J. Res. U.S. Geol. Surv., 2, 233

10.1029/WR017i005p01517

10.1002/aic.690300222

Wang H. F., 1982, Introduction to Groundwater Modeling

Westall J. C. J. L.Zachary F. M. M.Morel MINEQL: A computer Program for the Calculation of Chemical Equilibrium Composition of Aqueous SystemTech. Note 18 91Dep. of Civ. Eng. MIT Cambridge Mass. 1976.

10.1029/WR021i003p00272

10.1029/WR022i006p00952

Yeh G. T. FEMWATER: A Finite Element Model of WATER Flow through Saturated‐Unsaturated Porus Media First revisionRep. ORNL‐5567/R1Oak Ridge Nat. Lab. Oak Ridge Tenn. 37831 1987.

Yeh G. T., 1986, Mathematical strategies for hydrogeochemical transport models of reactive chemicals, Eos Trans. AGU, 67, 965

Yeh G. T. V. S.Tripathi A Lagrangian‐Eulerian approach to modeling hydrogeochemical transport of multi‐component systemsProceedings International Conference on Groundwater Contamination: Use of Models in Decision‐Making in the European Year of the EnvironmentInt. Groundwater Modeling Amsterdamthe NetherlandsOct. 26–29 1987.

Yeh G. T. V. S.Tripathi HYDROGEOCHEM: A Coupled Model of HYDROlogical Transport and GEOCHEMical Equilibrium of Multi‐component SystemsRep. ORNL‐6371Oak Ridge Nat. Lab. Oak Ridge Tenn. 1988.

Thông tin
Thông tin xuất bản

Water Resources Research

Tập 25 Số 1

93-108

Thông tin tác giả