On Vaporizing Water Flow in Hot Sub-Vertical Rock Fractures

Water injection into unsaturated fractured rock at above-boiling temperatures gives rise to complex fluid flow and heat transfer processes. Examples include water injection into depleted vapor-dominated geothermal reservoirs, and emplacement of heat-generating nuclear wastes in unsaturated fractured rock. We conceptualize fractures as two-dimensional heterogeneous porous media, and use geostatistical techniques to generate synthetic permeability distributions in the fracture plane. Water flow in hot high-angle fractures is simulated numerically, taking into account the combined action of gravity, capillary, and pressure forces, and conductive heat transfer from the wall rocks which gives rise to strong vaporization. In heterogeneous fractures boiling plumes are found to have dendritic shapes, and to be subject to strong lateral flow effects. Fractures with spatially-averaged homogeneous permeabilities tend to give poor approximations for vaporization behavior and liquid migration patterns. Depending on water flow rates, rock temperature, and fracture permeability, liquid water can migrate considerable distances through fractured rock that is at above-boiling temperatures and be only partially vaporized.

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

[2]  G. Matheron The intrinsic random functions and their applications , 1973, Advances in Applied Probability.

[3]  P. Witherspoon,et al.  Porous media equivalents for networks of discontinuous fractures , 1982 .

[4]  Y. Tsang,et al.  On Two-Phase Relative Permeability and Capillary Pressure of Rough-Walled Rock Fractures , 1989 .

[5]  T. N. Narasimhan,et al.  Aperture correlation of a fractal fracture , 1988 .

[6]  D. Tritton,et al.  Physical Fluid Dynamics , 1977 .

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

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

[9]  S. Wheatcraft,et al.  Gravity-driven infiltration instability in initially dry nonhorizontal fractures , 1994 .

[10]  A. Woods,et al.  The instability of a vaporization front in hot porous rock , 1994, Nature.

[11]  Sylvie Gentier,et al.  Morphologie et comportement hydromécanique d'une fracture naturelle dans un granite sous contrainte normale : étude expérimentale et théorique , 1986 .

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

[13]  G. I. Barenblatt,et al.  Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks [strata] , 1960 .

[14]  LawrenceD . Ramspott,et al.  The constructive use of heat in an unsaturated tuff repository , 1991 .

[15]  Shaun D. Fitzgerald,et al.  The vaporization of a liquid front moving through a hot porous rock , 1993, Journal of Fluid Mechanics.

[16]  Eva Hakami Water flow in single rock joints , 1988 .

[17]  I. J. Hall,et al.  Fracture and matrix hydrologic characteristics of tuffaceous materials from Yucca Mountain, Nye County, Nevada , 1984 .

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

[19]  T. N. Narasimhan,et al.  Multidimensional numerical simulation of fluid flow in fractured porous media , 1982 .

[20]  M. O'Sullivan,et al.  Reinjection studies of vapor-dominated systems , 1982 .

[21]  R. Ababou,et al.  Implementation of the three‐dimensional turning bands random field generator , 1989 .

[22]  J. S. Y. Wang,et al.  Validity of cubic law for fluid flow in a deformable rock fracture. Technical information report No. 23 , 1979 .

[23]  Y. Tsang,et al.  On thermohydrologic conditions near high-level nuclear wastes emplaced in partially saturated fractured tuff. 2. Effective continuum approximation , 1990 .

[24]  Injection recovery factors in various areas of the southeast Geysers, California , 1995 .

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

[26]  Scher,et al.  Simulation and theory of two-phase flow in porous media. , 1992, Physical review. A, Atomic, molecular, and optical physics.

[27]  P. A. Witherspoon,et al.  Hydraulic And Mechanical Properties of Natural Fractures In Low Permeability Rock , 1987 .

[28]  R. Celati,et al.  Modeling studies of cold water injection into fluid-depleted, vapor-dominated geothermal reservoirs , 1986 .

[29]  Yu-Shu Wu,et al.  An analytical solution for heat transfer at a boiling front moving through a porous medium , 1987, International Journal of Heat and Mass Transfer.

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

[31]  Alan L. Flint,et al.  Spatial distribution of potential near surface moisture flux at Yucca Mountain , 1994 .

[32]  C. Tsang,et al.  Geohydrologic Data and Models of Rainier Mesa and Their Implications to Yucca Mountain , 1993 .

[33]  R. W. Spengler,et al.  Geologic character of tuffs in the unsaturated zone at Yucca Mountain, southern Nevada , 1982 .

[34]  L. J. Pyrak-Nolte,et al.  The fractal geometry of flow paths in natural fractures in rock and the approach to percolation , 1989 .

[35]  J. J. Nitao,et al.  The analysis of repository-heat-driven hydrothermal flow at Yucca Mountain , 1993 .

[36]  J. W. Mercer,et al.  Role of the unsaturated zone in radioactive and hazardous waste disposal , 1983 .

[37]  Abbas Firoozabadi,et al.  Capillary Pressure in Fractured Porous Media (includes associated papers 21892 and 22212 ) , 1990 .

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

[39]  G. Bodvarsson,et al.  Thermal Effects of Reinjection in Geothermal Reservoirs With Major Vertical Fractures , 1984 .

[40]  Y. Tsang,et al.  MODELING OF STRONGLY HEAT-DRIVEN FLOW IN PARTIALLY SATURATED FRACTURED POROUS MEDIA , 1984 .

[41]  Thomas A. Buscheck,et al.  The implications of episodic nonequilibrium fracture-matrix flow on site suitability and total system performance , 1992 .

[42]  K. Pruess,et al.  TOUGH2-A General-Purpose Numerical Simulator for Multiphase Fluid and Heat Flow , 1991 .

[43]  K. Pruess,et al.  Thermal Modeling for a Potential High-level Nuclear Waste Repository at Yucca Mountain, Nevada , 1994 .

[44]  K. Pruess On water seepage and fast preferential flow in heterogeneous,unsaturated rock fractures , 1998 .

[45]  Shaun D. Fitzgerald,et al.  Viscous fingers in superheated geothermal systems , 1994 .

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

[47]  O. Dahan,et al.  Water Recharge and Solute Transport Through the Vadose Zone of Fractured Chalk Under Desert Conditions , 1995 .

[48]  G. Dagan Time‐dependent macrodispersion for solute transport in anisotropic heterogeneous aquifers , 1988 .

[49]  J. Gale,et al.  Comparison of coupled fracture deformation and fluid flow models with direct measurements of fracture pore structure and stress-flow properties , 1987 .

[50]  A. Mantoglou,et al.  The Turning Bands Method for simulation of random fields using line generation by a spectral method , 1982 .

[51]  Stefan Finsterle,et al.  Determining permeability of tight rock samples using inverse modeling , 1997 .

[52]  D. G. Wilder Alternative Strategies—A Means for Saving Money and Time on the Yucca Mountain Project , 1993 .

[53]  Y. Tsang,et al.  Modeling of strongly heat-driven flow processes at a potential high-level nuclear waste repository at Yucca Mountain, Nevada , 1993 .

[54]  K. Pruess,et al.  Numerical simulation of water injection into vapor-dominated reservoirs , 1995 .

[55]  Chin-Fu Tsang,et al.  Flow and Contaminant Transport in Fractured Rock , 1993 .

[56]  E. A. Klavetter,et al.  A continuum model for water movement in an unsaturated fractured rock mass , 1988 .

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

[58]  O. Phillips Infiltration of a liquid finger down a fracture into superheated rock , 1996 .

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

[60]  E. Kwicklis,et al.  Numerical investigation of steady liquid water flow in a variably saturated fracture network , 1993 .

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

[62]  J. C. Jaeger,et al.  Conduction of Heat in Solids , 1952 .

[63]  Karsten Pruess,et al.  Two‐Phase Flow Visualization and Relative Permeability Measurement in Natural Rough‐Walled Rock Fractures , 1995 .

[64]  Stephen R. Brown,et al.  Broad bandwidth study of the topography of natural rock surfaces , 1985 .

[65]  Roland Lenormand,et al.  Role Of Roughness And Edges During Imbibition In Square Capillaries , 1984 .

[66]  K. Pruess Grid orientation and capillary pressure effects in the simulation of water injection into depleted vapor zones , 1991 .

[67]  T. Narasimhan,et al.  7 – Unsaturated Flow in Fractured Porous Media , 1993 .

[68]  S. Enedy,et al.  Reservoir response to injection in the Southeast Geysers , 1991 .

[69]  Effective continuum approximation for modeling fluid and heat flow in fractured porous tuff: Nevada Nuclear Waste Storage Investigations Project , 1988 .