TOUGH2 User's Guide Version 2

TOUGH2 is a numerical simulator for nonisothermal flows of multicomponent, multiphase fluids in one, two, and three-dimensional porous and fractured media. The chief applications for which TOUGH2 is designed are in geothermal reservoir engineering, nuclear waste disposal, environmental assessment and remediation, and unsaturated and saturated zone hydrology. TOUGH2 was first released to the public in 1991; the 1991 code was updated in 1994 when a set of preconditioned conjugate gradient solvers was added to allow a more efficient solution of large problems. The current Version 2.0 features several new fluid property modules and offers enhanced process modeling capabilities, such as coupled reservoir-wellbore flow, precipitation and dissolution effects, and multiphase diffusion. Numerous improvements in previously released modules have been made and new user features have been added, such as enhanced linear equation solvers, and writing of graphics files. The T2VOC module for three-phase flows of water, air and a volatile organic chemical (VOC), and the T2DM module for hydrodynamic dispersion in 2-D flow systems have been integrated into the overall structure of the code and are included in the Version 2.0 package. Data inputs are upwardly compatible with the previous version. Coding changes were generally kept to a minimum, and were only made as needed to achieve the additional functionalities desired. TOUGH2 is written in standard FORTRAN77 and can be run on any platform, such as workstations, PCs, Macintosh, mainframe and supercomputers, for which appropriate FORTRAN compilers are available. This report is a self-contained guide to application of TOUGH2 to subsurface flow problems. It gives a technical description of the TOUGH2 code, including a discussion of the physical processes modeled, and the mathematical and numerical methods used. Illustrative sample problems are presented along with detailed instructions for preparing input data.

[1]  A. Herbert,et al.  Coupled groundwater flow and solute transport with fluid density strongly dependent upon concentration , 1988 .

[2]  D. W. Peaceman Interpretation of well-block pressures in numerical reservoir simulation with nonsquare grid blocks and anisotropic permeability , 1983 .

[3]  J. L. Yanosik,et al.  A Nine-Point, Finite-Difference Reservoir Simulator for Realistic Prediction of Adverse Mobility Ratio Displacements , 1979 .

[4]  George J. Moridis,et al.  TOUGH2 software qualification , 1996 .

[5]  S. W. Webb,et al.  Pore-scale modeling of enhanced vapor diffusion in porous media , 1997 .

[6]  K. Aziz,et al.  The Grid Orientation Effect in Reservoir Simulation , 1991 .

[7]  S. W. Webb Gas-phase diffusion in porous media: Evaluation of an advective- dispersive formulation and the dusty-gas model including comparison to data for binary mixtures , 1996 .

[8]  S. W. Webb,et al.  Review of enhanced vapor diffusion in porous media , 1998 .

[9]  Curtis M. Oldenburg,et al.  Higher-order differencing for phase-front propagation in geothermal systems , 1998 .

[10]  Robert W. Gillham,et al.  Finite-element analysis of the transport of water and solutes in title-drained soils , 1979 .

[11]  K. H. Coats,et al.  Geothermal Reservoir Modelling , 1977 .

[12]  G. W. Thomas,et al.  Principles of Hydrocarbon Reservoir Simulation , 1981 .

[13]  George J. Moridis,et al.  T2SOLV: An enhanced package of solvers for the TOUGH2 family of reservoir simulation codes , 1997 .

[14]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[15]  N. E. Edlefsen,et al.  Thermodynamics of soil moisture , 1943 .

[16]  G. J. Moridis,et al.  Flow and transport simulations using T2CG1, a package of conjugate gradient solvers for the TOUGH2 family of codes , 1995 .

[17]  P. Witherspoon,et al.  Study of two-phase concurrent flow of steam and water in an unconsolidated porous medium , 1985 .

[18]  M. O'Sullivan A Similarity Method for Geothermal Well Test Analysis , 1981 .

[19]  R. Corsi,et al.  Behavior of the Bagnore Steam/CO2 Geothermal Reservoir, Italy , 1980 .

[20]  P. Vaughan,et al.  Analysis of Permeability Reduction During Flow of Heated, Aqueous Fluid Through Westerly Granite , 1987 .

[21]  L. Klinkenberg The Permeability Of Porous Media To Liquids And Gases , 2012 .

[22]  E. A. Mason,et al.  Gas Transport in Porous Media: The Dusty-Gas Model , 1983 .

[23]  Karsten Pruess,et al.  Density-driven flow of gas in the unsaturated zone due to the evaporation of volatile organic compounds , 1989 .

[24]  M. C. Leverett,et al.  Capillary Behavior in Porous Solids , 1941 .

[25]  C. F. Curtiss,et al.  Molecular Theory Of Gases And Liquids , 1954 .

[26]  M. Reeves,et al.  Theory and implementation for SWIFT II. The Sandia waste-isolation flow and transport model for fractured media. Release 4. 84 , 1986 .

[27]  Robert W. Potter,et al.  The volumetric properties of aqueous sodium chloride solutions from 0 degrees to 500 degrees C at pressures up to 2000 bars based on a regression of the available literature data , 1975 .

[28]  T. Narasimhan,et al.  AN INTEGRATED FINITE DIFFERENCE METHOD FOR ANALYZING FLUID FLOW IN POROUS MEDIA , 1976 .

[29]  Karsten Pruess,et al.  A similarity solution for two‐phase water, air, and heat flow near a linear heat source in a porous medium , 1992 .

[30]  K. Aziz,et al.  Petroleum Reservoir Simulation , 1979 .

[31]  Y. Mualem A New Model for Predicting the Hydraulic Conductivity , 1976 .

[32]  J. Quirk,et al.  Permeability of porous solids , 1961 .

[33]  Y. Tsang,et al.  On thermohydrologic conditions near high‐level nuclear wastes emplaced in partially saturated fractured tuff: 1. Simulation studies with explicit consideration of fracture effects , 1990 .

[34]  Gudmundur S. Bodvarsson,et al.  FLUID AND HEAT FLOW IN GAS-RICH GEOTHERMAL RESERVOIRS , 1985 .

[35]  K. H. Coats,et al.  Effects of grid type and difference scheme on pattern steamflood simulation results , 1982 .

[36]  D. W. Peaceman Fundamentals of numerical reservoir simulation , 1977 .

[37]  P. Milly Moisture and heat transport in hysteretic, inhomogeneous porous media: A matric head‐based formulation and a numerical model , 1982 .

[38]  D. R. Fokkema,et al.  BiCGstab(ell) for Linear Equations involving Unsymmetric Matrices with Complex Spectrum , 1993 .

[39]  Iain S. Duff,et al.  MA28 --- A set of Fortran subroutines for sparse unsymmetric linear equations , 1980 .

[40]  Curtis M. Oldenburg,et al.  EOS7R: Radionuclide transport for TOUGH2 , 1995 .

[41]  K. Pruess,et al.  Modeling studies of gas movement and moisture migration at Yucca Mountain, Nevada , 1990 .

[42]  J. E. Warren,et al.  The Behavior of Naturally Fractured Reservoirs , 1963 .

[43]  G. Marsily Quantitative Hydrogeology: Groundwater Hydrology for Engineers , 1986 .

[44]  P. Vinsome,et al.  A Simple Method For Predicting Cap And Base Rock Heat Losses In' Thermal Reservoir Simulators , 1980 .

[45]  T. N. Narasimhan,et al.  A PRACTICAL METHOD FOR MODELING FLUID AND HEAT FLOW IN FRACTURED POROUS MEDIA , 1985 .

[46]  N. B. Vargaftik Tables on the thermophysical properties of liquids and gases: In normal and dissociated states , 1975 .

[47]  T. L. Jones,et al.  Enhancement of thermal water vapor diffusion in soil , 1984 .

[48]  K. Pruess,et al.  The simulator TOUGH2/EWASG for modelling geothermal reservoirs with brines and non-condensible gas , 1997 .

[49]  T. Narasimhan,et al.  On fluid reserves and the production of superheated steam from fractured, vapor‐dominated geothermal reservoirs , 1982 .

[50]  William R. Walker,et al.  Studies of heat transfer and water migration in soils. Final report , 1981 .

[51]  L. A. Richards Capillary conduction of liquids through porous mediums , 1931 .

[52]  E. Michaelides Thermodynamic properties of geothermal fluids , 1981 .