Numerical studies of fluid and heat flow near high-level nuclear waste packages emplaced in partially saturated fractured tuff

We have performed modeling studies on the simultaneous transport of heat, liquid water, vapor, and air in partially saturated fractured porous rock. Formation parameters were chosen as representative of the potential repository horizon in the Topopah Spring Unit of the Yucca Mountain tuffs. The presence of fractures makes the transport problem very complex, both in terms of flow geometry and physics. The numerical simulator ''TOUGH'' used for our flow calculations takes into account most of the physical effects which are important in multi-phase fluid and heat flow. It has provisions for handling the extreme non-linearities which arise in phase transitions, component disappearances, and capillary discontinuities at fracture faces. We model a region around an infinite linear string of nuclear waste canisters, taking into account both the discrete fractures and the porous matrix. From an analysis of the results obtained with explicit fractures, we develop equivalent continuum models which can reproduce the temperature, saturation, and pressure variation, and gas and liquid flow rates of the discrete fracture-porous matrix calculations. The equivalent continuum approach makes use of a generalized relative permeability concept to take into account the fracture effects. This results in a substantial simplification of the flow problem which makes larger scale modeling of complicated unsaturated fractured porous systems feasible. Potential applications for regional scale simulations and limitations of the continuum approach are discussed. 35 refs., 14 figs., 4 tabs