skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Modeling non-isothermal multiphase multi-species reactive chemical transport in geologic media

Abstract

The assessment of mineral deposits, the analysis of hydrothermal convection systems, the performance of radioactive, urban and industrial waste disposal, the study of groundwater pollution, and the understanding of natural groundwater quality patterns all require modeling tools that can consider both the transport of dissolved species as well as their interactions with solid (or other) phases in geologic media and engineered barriers. Here, a general multi-species reactive transport formulation has been developed, which is applicable to homogeneous and/or heterogeneous reactions that can proceed either subject to local equilibrium conditions or kinetic rates under non-isothermal multiphase flow conditions. Two numerical solution methods, the direct substitution approach (DSA) and sequential iteration approach (SIA) for solving the coupled complex subsurface thermo-physical-chemical processes, are described. An efficient sequential iteration approach, which solves transport of solutes and chemical reactions sequentially and iteratively, is proposed for the current reactive chemical transport computer code development. The coupled flow (water, vapor, air and heat) and solute transport equations are also solved sequentially. The existing multiphase flow code TOUGH2 and geochemical code EQ3/6 are used to implement this SIA. The flow chart of the coupled code TOUGH2-EQ3/6, required modifications of the existing codes and additional subroutines needed are presented.

Authors:
; ; ;
Publication Date:
Research Org.:
Lawrence Livermore National Lab., CA (United States)
Sponsoring Org.:
USDOE Office of Energy Research, Washington, DC (United States)
OSTI Identifier:
589221
Report Number(s):
LBNL-40504
ON: DE98052709; TRN: 98:001859
DOE Contract Number:
AC03-76SF00098
Resource Type:
Technical Report
Resource Relation:
Other Information: PBD: 1 Jul 1997
Country of Publication:
United States
Language:
English
Subject:
54 ENVIRONMENTAL SCIENCES; 58 GEOSCIENCES; MULTIPHASE FLOW; POROUS MATERIALS; FLOW MODELS; E CODES; T CODES; COMPUTERIZED SIMULATION; ENVIRONMENTAL TRANSPORT; GROUND WATER; RADIONUCLIDE MIGRATION; LAWRENCE BERKELEY LABORATORY; RESEARCH PROGRAMS

Citation Formats

Tianfu Xu, Gerard, F., Pruess, K., and Brimhall, G.. Modeling non-isothermal multiphase multi-species reactive chemical transport in geologic media. United States: N. p., 1997. Web. doi:10.2172/589221.
Tianfu Xu, Gerard, F., Pruess, K., & Brimhall, G.. Modeling non-isothermal multiphase multi-species reactive chemical transport in geologic media. United States. doi:10.2172/589221.
Tianfu Xu, Gerard, F., Pruess, K., and Brimhall, G.. 1997. "Modeling non-isothermal multiphase multi-species reactive chemical transport in geologic media". United States. doi:10.2172/589221. https://www.osti.gov/servlets/purl/589221.
@article{osti_589221,
title = {Modeling non-isothermal multiphase multi-species reactive chemical transport in geologic media},
author = {Tianfu Xu and Gerard, F. and Pruess, K. and Brimhall, G.},
abstractNote = {The assessment of mineral deposits, the analysis of hydrothermal convection systems, the performance of radioactive, urban and industrial waste disposal, the study of groundwater pollution, and the understanding of natural groundwater quality patterns all require modeling tools that can consider both the transport of dissolved species as well as their interactions with solid (or other) phases in geologic media and engineered barriers. Here, a general multi-species reactive transport formulation has been developed, which is applicable to homogeneous and/or heterogeneous reactions that can proceed either subject to local equilibrium conditions or kinetic rates under non-isothermal multiphase flow conditions. Two numerical solution methods, the direct substitution approach (DSA) and sequential iteration approach (SIA) for solving the coupled complex subsurface thermo-physical-chemical processes, are described. An efficient sequential iteration approach, which solves transport of solutes and chemical reactions sequentially and iteratively, is proposed for the current reactive chemical transport computer code development. The coupled flow (water, vapor, air and heat) and solute transport equations are also solved sequentially. The existing multiphase flow code TOUGH2 and geochemical code EQ3/6 are used to implement this SIA. The flow chart of the coupled code TOUGH2-EQ3/6, required modifications of the existing codes and additional subroutines needed are presented.},
doi = {10.2172/589221},
journal = {},
number = ,
volume = ,
place = {United States},
year = 1997,
month = 7
}

Technical Report:

Save / Share:
  • Coupled modeling of subsurface multiphase fluid and heat flow, solute transport and chemical reactions can be used for the assessment of mineral alteration in hydrothermal systems, waste disposal sites, acid mine drainage remediation, contaminant transport, and groundwater quality. A comprehensive non-isothermal multi-component reactive fluid flow and geochemical transport simulator, TOUGHREACT, has been developed. A wide range of subsurface thermo-physical-chemical processes is considered under various thermohydrological and geochemical conditions of pressure, temperature, water saturation, and ionic strength. The program can be applied to one-, two- or three-dimensional porous and fractured media with physical and chemical heterogeneity. The model can accommodate anymore » number of chemical species present in liquid, gas and solid phases. A variety of equilibrium chemical reactions are considered, such as aqueous complexation, gas dissolution/exsolution, and cation exchange. Mineral dissolution/precipitation can proceed either subject to local equilibrium or kinetic conditions. Changes in porosity and permeability due to mineral dissolution and precipitation can be considered. Linear adsorption and decay can be included. For the purpose of future extensions, surface complexation by double layer model is coded in the program. Xu and Pruess (1998) developed a first version of a non-isothermal reactive geochemical transport model, TOUGHREACT, by introducing reactive geochemistry into the framework of the existing multi-phase fluid and heat flow code TOUGH2 (Pruess, 1991). Xu, Pruess, and their colleagues have applied the program to a variety of problems such as: (1) supergene copper enrichment (Xu et al, 2001), (2) caprock mineral alteration in a hydrothermal system (Xu and Pruess, 2001a), and (3) mineral trapping for CO{sub 2} disposal in deep saline aquifers (Xu et al, 2003b and 2004a). For modeling the coupled thermal, hydrological, and chemical processes during heater tests at proposed nuclear waste disposal site at Yucca Mountain (Nevada), Sonnenthal and Spycher (2000) and Spycher et al. (2003) enhanced TOUGHREACT on (1) high temperature geochemistry, (2) mineral reactive surface area calculations, and (3) porosity and permeability changes due to mineral alteration. On the other hand, Pruess et al. (1999) updated the TOUGH2 simulator to TOUGH2 V2. The present version of TOUGHREACT was developed by introducing the work of Sonnenthal and Spycher (2000) to the original work of Xu and Pruess (1998), and by replacing TOUGH2 (Pruess, 1991) by TOUGH2 V2 (Pruess et al, 1999). The TOUGHREACT program makes use of ''self-documenting'' features. It is distributed with a number of input data files for sample problems. Besides providing benchmarks for proper code installation, these can serve as self-teaching tutorial in the use of TOUGHREACT, and they provide templates to help jump-start new applications. The fluid and heat flow part of TOUGHREACT is derived from TOUGH2 V2, so in addition to the current manual, users must have manual of the TOUGH2 V2 (Pruess et al., 1999). The present version of TOUGHREACT provides the following different TOUGH2 fluid property or ''EOS'' (equation-of-state) modules: (1) EOS1 for water, or two waters with typical applications to hydrothermal problems, (2) EOS2 for multiphase mixtures of water and CO{sub 2} also with typical applications to hydrothermal problems, (3) EOS3 for multiphase mixtures of water and air with typical applications to vadose zone and nuclear waste disposal problems, (4) EOS4 that has the same capabilities as EOS3 but with vapor pressure lowering effects due to capillary pressure, (5) EOS9 for single phase water (Richards. equation) with typical applications to ambient reactive geochemical transport problems, (6) ECO2 for multiphase mixtures of water, CO{sub 2} and NaCl with typical applications to CO{sub 2} disposal in deep brine aquifers.« less
  • Coupled modeling of subsurface multiphase fluid and heat flow, solute transport, and chemical reactions can be applied to many geologic systems and environmental problems, including geothermal systems, diagenetic and weathering processes, subsurface waste disposal, acid mine drainage remediation, contaminant transport, and groundwater quality. TOUGHREACT has been developed as a comprehensive non-isothermal multi-component reactive fluid flow and geochemical transport simulator to investigate these and other problems. A number of subsurface thermo-physical-chemical processes are considered under various thermohydrological and geochemical conditions of pressure, temperature, water saturation, and ionic strength. TOUGHREACT can be applied to one-, two- or three-dimensional porous and fractured mediamore » with physical and chemical heterogeneity. The code can accommodate any number of chemical species present in liquid, gas and solid phases. A variety of equilibrium chemical reactions are considered, such as aqueous complexation, gas dissolution/exsolution, and cation exchange. Mineral dissolution/precipitation can take place subject to either local equilibrium or kinetic controls, with coupling to changes in porosity and permeability and capillary pressure in unsaturated systems. Chemical components can also be treated by linear adsorption and radioactive decay. The first version of the non-isothermal reactive geochemical transport code TOUGHREACT was developed (Xu and Pruess, 1998) by introducing reactive geochemistry into the framework of the existing multi-phase fluid and heat flow code TOUGH2 (Pruess, 1991). TOUGHREACT was further enhanced with the addition of (1) treatment of mineral-water-gas reactive-transport under boiling conditions, (2) an improved HKF activity model for aqueous species, (3) gas species diffusion coefficients calculated as a function of pressure, temperature, and molecular properties, (4) mineral reactive surface area formulations for fractured and porous media, and (5) porosity, permeability, and capillary pressure changes owing to mineral precipitation/dissolution (Sonnenthal et al., 1998, 2000, 2001; Spycher et al., 2003a). Subsequently, TOUGH2 V2 was released with additional EOS modules and features (Pruess et al., 1999). The present version of TOUGHREACT includes all of the previous extensions to the original version, along with the replacement of the original TOUGH2 (Pruess, 1991) by TOUGH2 V2 (Pruess et al., 1999). TOUGHREACT has been applied to a wide variety of problems, some of which are included as examples, such as: (1) Supergene copper enrichment (Xu et al., 2001); (2) Mineral alteration in hydrothermal systems (Xu and Pruess, 2001a; Xu et al., 2004b; Dobson et al., 2004); (3) Mineral trapping for CO{sub 2} disposal in deep saline aquifers (Xu et al., 2003b and 2004a); (4) Coupled thermal, hydrological, and chemical processes in boiling unsaturated tuff for the proposed nuclear waste emplacement site at Yucca Mountain, Nevada (Sonnenthal et al., 1998, 2001; Sonnenthal and Spycher, 2000; Spycher et al., 2003a, b; Xu et al., 2001); (5) Modeling of mineral precipitation/dissolution in plug-flow and fracture-flow experiments under boiling conditions (Dobson et al., 2003); (6) Calcite precipitation in the vadose zone as a function of net infiltration (Xu et al., 2003); and (7) Stable isotope fractionation in unsaturated zone pore water and vapor (Singleton et al., 2004). The TOUGHREACT program makes use of 'self-documenting' features. It is distributed with a number of input data files for sample problems. Besides providing benchmarks for proper code installation, these can serve as a self-teaching tutorial in the use of TOUGHREACT, and they provide templates to help jump-start new applications. The fluid and heat flow part of TOUGHREACT is derived from TOUGH2 V2, so in addition to the current manual, users must have the manual of the TOUGH2 V2 (Pruess et al., 1999). The present version of TOUGHREACT provides the following TOUGH2 fluid property or 'EOS' (equation-of-state) modules: (1) EOS1 for water, or two waters with typical applications to hydrothermal problems, (2) EOS2 for multiphase mixtures of water and CO{sub 2} also with typical applications to hydrothermal problems, (3) EOS3 for multiphase mixtures of water and air with typical applications to vadose zone and nuclear waste disposal problems, (4) EOS4 that has the same capabilities as EOS3 but with vapor pressure lowering effects due to capillary pressure, (5) EOS9 for single phase water (Richards equation) with typical applications to ambient temperature and pressure reactive geochemical transport problems, and (6) ECO2N for multiphase mixtures of water, CO{sub 2} and NaCl with typical applications to CO{sub 2} disposal in deep brine aquifers.« less
  • Coupled modeling of subsurface multiphase fluid and heat flow, solute transport and chemical reactions can be used for the assessment of acid mine drainage remediation, mineral deposition, waste disposal sites, hydrothermal convection, contaminant transport, and groundwater quality. Here they present a numerical simulation model, TOUGHREACT, which considers non-isothermal multi-component chemical transport in both liquid and gas phases. A wide range of subsurface thermo-physical-chemical processes is considered. The model can be applied to one-, two- or three-dimensional porous and fractured media with physical and chemical heterogeneity. The model can accommodate any number of chemical species present in liquid, gas and solidmore » phases. A variety of equilibrium chemical reactions is considered, such as aqueous complexation, gas dissolution/exsolution, cation exchange, and surface complexation. Mineral dissolution/precipitation can proceed either subject to local equilibrium or kinetic conditions. The coupled model employs a sequential iteration approach with reasonable computing efficiency. The development of the governing equations and numerical approach is presented along with the discussion of the model implementation and capabilities. The model is verified for a wide range of subsurface physical and chemical processes. The model is well suited for flow and reactive transport in variably saturated porous and fractured media. In the second of this two-part paper, three applications covering a variety of problems are presented to illustrate the capabilities of the model.« less
  • Coupled modeling of subsurface multiphase fluid and heat flow, solute transport and chemical reactions can be used for the assessment of acid mine drainage remediation, waste disposal sites, hydrothermal convection, contaminant transport, and groundwater quality. We have developed a comprehensive numerical simulator, TOUGHREACT, which considers non-isothermal multi-component chemical transport in both liquid and gas phases. A wide range of subsurface thermo-physical-chemical processes is considered under various thermohydrological and geochemical conditions of pressure, temperature, water saturation, and ionic strength. The code can be applied to one-, two- or three-dimensional porous and fractured media with physical and chemical heterogeneity.