HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Full Text Full Text (PDF) Services Similar articles in this journal Similar articles in Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Web of Science (43) Citing Articles via Google Scholar Google Scholar Articles by Xu, T. Articles by Pruess, K. Search for Related Content Articles by Xu, T. Articles by Pruess, K. GeoRef GeoRef Citation

Modeling Multiphase Non-isothermal Fluid Flow and Reactive Geochemical Transport in Variably Saturated Fractured Rocks: 1. Methodology

Earth Sciences Division, Lawrence Berkeley National Laboratory, University of California, Berkeley, California 94720

Reactive fluid flow and geochemical transport in unsaturated fractured rocks have received increasing attention for studies of contaminant transport, ground- water quality, waste disposal, acid mine drainage remediation, mineral deposits, sedimentary diagenesis, and fluid-rock interactions in hydrothermal systems. This paper presents methods for modeling geochemical systems that emphasize: (1) involvement of the gas phase in addition to liquid and solid phases in fluid flow, mass transport, and chemical reactions; (2) treatment of physically and chemically heterogeneous and fractured rocks, (3) the effect of heat on fluid flow and reaction properties and processes, and (4) the kinetics of fluid-rock interaction. The physical and chemical process model is embodied in a system of partial differential equations for flow and transport, coupled to algebraic equations and ordinary differential equations for chemical interactions. For numerical solution, the continuum equations are discretized in space and time. Space discretization is based on a flexible integral finite difference approach that can use irregular gridding to model geologic structure; time is discretized fully implicitly as a first-order finite difference. Heterogeneous and fractured media are treated with a general multiple interacting continua method that includes double-porosity, dual-permeability, and multi-region models as special cases. A sequential iteration approach is used to treat the coupling between fluid flow and mass transport on the one hand, chemical reactions on the other. Applications of the methods developed here to variably saturated geochemical systems are presented in a companion paper (part 2, this issue).

This article has been cited by other articles:

				J. M. Mogollon, I. L'Heureux, A. W. Dale, and P. Regnier Methane gas-phase dynamics in marine sediments: A model study Am J Sci, March 1, 2009; 309(3): 189 - 220. [Abstract] [Full Text] [PDF]

				P. Audigane, I. Gaus, I. Czernichowski-Lauriol, K. Pruess, and T. Xu Two-dimensional reactive transport modeling of CO2 injection in a saline aquifer at the Sleipner site, North Sea Am J Sci, September 1, 2007; 307(7): 974 - 1008. [Abstract] [Full Text] [PDF]

				G. D. Jones and Y. Xiao Geothermal convection in the Tengiz carbonate platform, Kazakhstan: Reactive transport models of diagenesis and reservoir quality AAPG Bulletin, August 1, 2006; 90(8): 1251 - 1272. [Abstract] [Full Text] [PDF]

				J. Bartels, C. Clauser, M. Kuhn, H. Pape, and W. Schneider Reactive flow and permeability prediction -- numerical simulation of complex hydrogeothermal problems Geological Society, London, Special Publications, January 1, 2005; 240(1): 133 - 151. [Abstract] [PDF]

				T. Xu, E. Sonnenthal, N. Spycher, K. Pruess, G. Brimhall, and J. Apps Modeling Multiphase Non-Isothermal Fluid Flow and Reactive Geochemical Transport in Variably Saturated Fractured Rocks: 2. Applications to Supergene Copper Enrichment and Hydrothermal Flows Am J Sci, January 1, 2001; 301(1): 34 - 59. [Abstract] [Full Text] [PDF]