A mountain-scale thermal-hydrologic model for simulating fluid flow and heat transfer in unsaturated fractured rock

A multidimensional, mountain-scale, thermal-hydrologic (TH) numerical model is presented for investigating unsaturated flow behavior in response to decay heat from the proposed radioactive waste repository in the Yucca Mountain unsaturated zone (UZ), The model, consisting of both two-dimensional (2-D) and three-dimensional (3-D) representations of the UZ repository system, is based on the current repository design, drift layout, thermal loading scenario, and estimated current and future climate conditions. This mountain-scale TH model evaluates the coupled TH processes related to mountain-scale UZ flow. It also simulates the impact of radioactive waste heat release on the natural hydrogeological system, including heat-driven processes occurring near and far away from the emplacement tunnels or drifts. The model simulates predict thermally perturbed liquid saturation, gas- and liquid-phase fluxes, and water and rock temperature elevations, as well as the changes in water flux driven by evaporation/condensation processes and drainage between drifts. These simulations provide insights into mountain-scale thermally perturbed flow fields under thermal loading conditions.

An active region model for capturing fractal flow patterns in unsaturated soils: model development

Preferential flow commonly observed in unsaturated soils allows rapid movement of solute from the soil surface or vadose zone to the groundwater, bypassing a significant volume of unsaturated soil and increasing the risk of groundwater contamination. A variety of evidence indicates that complex preferential patterns observed from fields are fractals. In this study, we developed a relatively simple active region model to incorporate the fractal flow pattern into the continuum approach. In the model, the flow domain is divided into active and inactive regions. Flow occurs preferentially in the active region (characterized by fractals), and inactive region is simply bypassed. A new constitutive relationship (the portion of the active region as a function of saturation) was derived. The validity of the proposed model is demonstrated by the consistency between field observations and the new constitutive relationship.

Field investigation into unsaturated flow and transport in a fault: model analyses

Results of a fault test performed in the unsaturated zone of Yucca Mountain, Nevada, were analyzed using a three-dimensional numerical model. The fault was explicitly represented as a discrete feature and the surrounding rock was treated as a dual-continuum (fracture-matrix) system. Model calibration against seepage and water-travel-velocity data suggests that lithophysal cavities connected to fractures can considerably enhance the effective fracture porosity and therefore retard water flow in fractures. Comparisons between simulation results and tracer concentration data also indicate that matrix diffusion is an important mechanism for solute transport in unsaturated fractured rock. We found that an increased fault-matrix and fracture-matrix interface areas were needed to match the observed tracer data, which is consistent with previous studies. The study results suggest that the current site-scale model for the unsaturated zone of Yucca Mountain may underestimate radionuclide transport time within the unsaturated zone, because an increased fracture-matrix interface area and the increased effective fracture porosity arising from lithophysal cavities are not considered in the current site-scale model.

A triple-continuum approach for modeling flow and transport processes in fractured rock

This paper presents a triple-continuum conceptual model for simulating flow and transport processes in fractured rock. Field data collected from the unsaturated zone of Yucca Mountain, a repository site of high-level nuclear waste, show a large number of small-scale fractures. The effect of these small fractures has not been considered in previous modeling investigations within the context of a continuum approach. A new triple-continuum model (consisting of matrix, small-fracture, and large-fracture continua) has been developed to investigate the effect of these small fractures. This paper derives the model formulation and discusses the basic triple-continuum behavior of flow and transport processes under different conditions, using both analytical solutions and numerical approaches. The simulation results from the site-scale model of the unsaturated zone of Yucca Mountain indicate that these small fractures may have an important effect on radionuclide transport within the mountain.

Conceptual evaluation of the potential role of fractures in unsaturated processes at Yucca Mountain

A wide array of field observations, in situ testing, and rock and water sampling (and subsequent analyses) within the unsaturated zone (UZ) of Yucca Mountain demonstrate the importance of fractures to flow and transport in the welded tuffs. The abundance of fractures and the spatial variability in their hydraulic properties, along with the heterogeneity within lithologic formations, make evaluation of unsaturated processes occurring within the potential repository horizon complex. Fracture mapping and field testing show that fractures are well connected, yet considerable variation is seen within and between units comprising the potential repository horizon with regard to fracture trace length, spacing, permeability, and capillarity. These variations have important implications for the distribution and movement of water and solutes through the unsaturated zone. Numerical models designed to assess such phenomena as unsaturated flow, transport, and coupled thermal-hydrological processes each require their own conceptual model for fracture networks, in order to identify the subset of all fractures that is relevant to the particular study. We evaluate several process-dependent conceptual models for fractures and identify the relevant fracture subsets related to these processes.

Modeling thermal-hydrological response of the unsaturated zone at Yucca Mountain, Nevada, to thermal load at a potential repository

This paper presents a numerical study on the response of the unsaturated zone (UZ) system of Yucca Mountain to heat generated from decaying radioactive wastes emplaced at the proposed repository. The modeling study is based on the current thermal-hydrological (TH) mountain-scale model, which uses a locally refined 2D north-south cross-section and dual-permeability numerical approach. The model provides a prediction of the mountain-scale TH response under the thermal-load scenario of 1.45 kW/m, while accounting for future climatic changes and the effects of drift ventilation. The TH simulation results show that ventilation of the repository drifts has a large impact on thermal-hydrologic regimes and moisture-flow conditions at the repository. In both cases, with and without ventilation, the TH model predicts dry or reduced liquid saturation near the drifts for over 1000 years, during which liquid flux through the drifts is reduced to either zero or less than the ambient flux. Without ventilation, the model predicts higher temperatures at the repository, but no major moisture redistribution in the UZ except in the areas very near the heated drifts.

Evolution of the unsaturated zone testing at Yucca Mountain

The evaluation of the Yucca Mountain site has evolved from intensive surface-based investigations in the early 1980s to current focus on testing in underground drifts. Different periods of site characterization activities and prominent issues concerning the unsaturated zone (UZ) are summarized. Data collection activities have evolved from mapping of faults and fractures to estimation of percolation through tuff layers, and to quantification of seepage into drifts. Evaluation of discrete flow paths in drifts has led to fracture-matrix interaction and matrix diffusion tests over different scales. The effects of tuff interfaces and local faults are evaluated in fractured-welded and porous-nonwelded units. Mobilization of matrix water and redistribution of moisture are measured in thermal tests. Lessons learned from underground tests are used to focus on processes needed for additional quantification. Migration through the drift shadow zone and liquid flow through faults are two important issues that have evolved from current knowledge.

Parallel computing simulation of fluid flow in the unsaturated zone of Yucca Mountain, Nevada

This paper presents the application of parallel computing techniques to large-scale modeling of fluid flow in the unsaturated zone (UZ) at Yucca Mountain, Nevada. In this study, parallel computing techniques, as implemented into the TOUGH2 code, are applied in large-scale numerical simulations on a distributed-memory parallel computer. The modeling study has been conducted using an over-1-million-cell three-dimensional numerical model, which incorporates a wide variety of field data for the highly heterogeneous fractured formation at Yucca Mountain. The objective of this study is to analyze the impact of various surface infiltration scenarios (under current and possible future climates) on flow through the UZ system, using various hydrogeological conceptual models with refined grids. The results indicate that the 1-million-cell models produce better resolution results and reveal some flow patterns that cannot be obtained using coarse-grid modeling models.

Estimation of percolation flux from borehole temperature data at Yucca Mountain, Nevada

Temperature data from the unsaturated zone (UZ) at Yucca Mountain are analyzed to estimate percolation-flux rates and overall heat flux. A multilayer, one-dimensional analytical solution is presented for determining percolation flux from temperature data. Case studies have shown that the analytical solution agrees very well with results from the numerical code, TOUGH2. The results of the analysis yield percolation fluxes in the range from 0 to 20 mm/year for most of the deep boreholes. This range is in good agreement with the results of infiltration studies at Yucca Mountain. Percolation flux for the shallower boreholes, however, cannot be accurately determined from temperature data alone because large gas flow in the shallow system alters the temperature profiles. Percolation-flux estimates for boreholes located near or intersecting major faults are significantly higher than those for other boreholes. These estimates may be affected by gas flow in the faults.

Analysis of flow behavior in fractured lithophysal reservoirs

This study develops a mathematical model for the analysis of pressure behavior in fractured lithophysal reservoirs. The lithophysal rock is described as a tri-continuum medium, consisting of fractures, rock matrices, and cavities. In the conceptual model, fractures have homogeneous properties throughout and interact with rock matrices and cavities that have different permeabilities and porosities. Global flow occurs through the fracture network only, while rock matrices and cavities contain the majority of fluid storage and provide fluid drainage to the fractures. Interporosity flows between the triple media are described using a pseudosteady-state concept and the system is characterized by interporosity transmissivity ratios and storativity ratio of each continuum. Pressure behavior is analyzed by examining the pressure drawdown curves, the derivative plots, and the effects of the characteristic parameters. Typical pressure responses from fractures, matrices, and cavities are represented by three semilog straight lines; the transitions by two troughs below the stabilization lines in the derivative plots. The analytical solution to the proposed model is further verified using a numerical simulation. The analytical model has also been applied to a published field-buildup well test and is able to match the pressure buildup data.

Development of discrete flow paths in unsaturated fractures at Yucca Mountain

We have carried out numerical modeling studies to investigate the development of discrete-fracture flow paths and flow-focusing phenomena in the unsaturated rock of the potential repository horizon at Yucca Mountain, Nevada. These studies are based on two- and three-dimensional (2-D and 3-D) numerical models using site-specific parameters. The 2-D and 3-D models use high-resolution spatial discretization to explicitly include effects of discrete fractures with stochastically developed fracture permeabilities and a continuum approach. The permeability field is generated based on air permeability measurements at various scales. For most of the cases considered, uniform infiltration with different average rates (1-500 mm/year) is prescribed at the top of the model, while variability in outflow at the bottom of the model is used to evaluate the degree of flow focusing. In addition, scenarios involving nonuniform infiltration at the top boundary, different permeability correlation lengths and different flow-allocation schemes were analyzed. The modeling results obtained from all of the cases showed a remarkably similar flow-focusing pattern at the repository horizon. Furthermore, tracer transport simulation results also revealed additional features of focused flow and transport through the fracture network.

On the physics of unstable infiltration, seepage, and gravity drainage in partially saturated tuffs

To improve understanding of the physics of dynamic instabilities in unsaturated flow processes within the Paintbrush nonwelded unit (PTn) and the middle nonlithophysal portion of the Topopah Spring welded tuff unit (TSw) of Yucca Mountain, we analyzed data from a series of infiltration tests carried out at two sites (Alcove 4 and Alcove 6) in the Exploratory Studies Facility (ESF), using analytical and empirical functions. The analysis of infiltration rates measured at both sites showed three temporal scales of infiltration rate: (1) a macro-scale trend of overall decreasing flow, (2) a meso-scale trend of fast and slow motion exhibiting three-stage variations of the flow rate (decreasing, increasing, and [again] decreasing flow rate, as observed in soils in the presence of entrapped air), and (3) micro-scale (high frequency) fluctuations. Infiltration tests in the nonwelded unit at Alcove 4 indicate that this unit may effectively dampen episodic fast infiltration events; however, well-known Kostyakov, Horton, and Philip equations do not satisfactorily describe the observed trends of the infiltration rate. Instead, a Weibull distribution model can most accurately describe experimentally determined time trends of the infiltration rate. Infiltration tests in highly permeable, fractured, welded tuff at Alcove 6 indicate that the infiltration rate exhibits pulsation, which may have been caused by multiple threshold effects and water-air redistribution between fractures and matrix. The empirical relationships between the extrinsic seepage from fractures, matrix imbibition, and gravity drainage versus the infiltration rate, as well as scaling and self-similarity for the leading edge of the water front are the hallmark of the nonlinear dynamic processes in water flow under episodic infiltration through fractured tuff. Based on the analysis of experimental data, we propose a conceptual model of a dynamic fracture flow and fracture-matrix interaction in fractured tuff, incorporating the time-dependent processes of water redistribution in the fracture-matrix system.

Flow and transport in the drift shadow in a dual-continuum model

The current concept for high-level radioactive waste disposal at Yucca Mountain is for the waste to be placed in underground tunnels (or drifts) in the middle of a thick unsaturated zone. Flow modeling and field testing have shown that not all flow encountering a drift will seep into the drift. The underlying reason for the diversion of unsaturated flow around a drift is that capillary forces in the fractures and matrix prevent water entry into the drift unless the capillary pressure in the rock decreases sufficiently to allow for gravity forces to overcome the capillary barrier. As a result of the capillary barrier effect, flow tends to be diverted around the drift, affecting the flow pattern beneath the drift. For some distance beneath the drift, water saturation and flux are reduced. This drift shadow zone is much more pronounced in the fractures than in the matrix due to dominance of gravity over capillary forces in the fractures. Moving downward, away from the drift, the shadow zone asymptotically re-equilibrates to the undisturbed flow conditions due to capillary forces. The behavior of radionuclide transport in this zone of reduced flow is investigated here because this will affect the amount of time required for radionuclides to penetrate the unsaturated zone. The delay of radionuclide movement in the geosphere is one aspect of the potential repository system that could limit public exposure to radioactive waste. The behavior of flow and transport is calculated using a two-dimensional, drift-scale dual-permeability model extending to nine drift diameters below the potential waste emplacement drift. The flow model is first compared with an analytical model for a single continuum. Then, the dual-continuum flow model is investigated with respect to drift-scale and mountain-scale property sets. Transport calculations are performed for a wide range of flow conditions and for different aqueous radionuclides and colloids. Findings indicate that transport times for dissolved or colloidal material released from a drift without seepage are several orders of magnitude longer than if the releases occurred in the undisturbed flow field. Furthermore, the calculations indicate that the transport rate for radionuclides released in the drift shadow is relatively insensitive to flow rates in the fractures, but is sensitive to the flow rate in the matrix.

Effect of heterogeneity in fracture permeability on the potential for liquid seepage into a heated emplacement drift of the potential repository

A numerical model was used to investigate the effect of spatial variability in fracture permeability on liquid seepage and moisture distribution in the vicinity of a waste emplacement drift in the unsaturated zone (UZ) of Yucca Mountain. The model is based on a two-dimensional, cross-sectional, dual-permeability model of the unsaturated zone at Yucca Mountain and uses a stochastic approach to investigate the effect of small-scale heterogeneous features. The studies were conducted using one uniform fracture permeability case, three realizations of stochastically generated fracture permeability, one discrete permeability feature case, and one increased ambient liquid flux case. In all cases, the models predict that completely dry drift conditions will develop above and below the drift in 10-100 years and remain dry for 1000-2000 years. During this period, the models predict no seepage into drifts, although liquid flux above the drifts and within the drift pillars may increase by up to two orders of magnitude above ambient flux. This is because the heat released by the emplaced waste is sufficient to vaporize liquid flux of one to two orders of magnitude higher than present-day ambient flux for over 1000 years. The results also show that unsaturated zone thermal-hydrological (TH) models with uniform layer permeability can adequately predict the evolution of seepage and moisture distribution in the rock mass surrounding the repository drifts. The models further show that although variability in fracture permeability may focus and enhance liquid flow in regions of enhanced liquid saturation (due to condensation above the drifts), vaporization and vapor diffusion can maintain a dry environment within the drifts for thousands of years.

Preliminary 3-D site-scale studies of radioactive colloid transport in the unsaturated zone at Yucca Mountain, Nevada

The U.S. Department of Energy (DOE) is actively investigating the technical feasibility of permanent disposal of high-level nuclear waste in a repository to be situated in the unsaturated zone (UZ) at Yucca Mountain (YM), Nevada. In this study we investigate, by means of numerical simulation, the transport of radioactive colloids under ambient conditions from the potential repository horizon to the water table. The site hydrology and the effects of the spatial distribution of hydraulic and transport properties in the Yucca Mountain subsurface are considered. The study of migration and retardation of colloids accounts for the complex processes in the unsaturated zone of Yucca Mountain, and includes advection, diffusion, hydrodynamic dispersion, kinetic colloid filtration, colloid straining, and radioactive decay. The results of the study indicate that the most important factors affecting colloid transport are the subsurface geology and site hydrology, i.e., the presence of faults (they dominate and control transport), fractures (the main migration pathways), and the relative distribution of zeolitic and vitric tuffs. The transport of colloids is strongly influenced by their size (as it affects diffusion into the matrix, straining at hydrogeologic unit interfaces, and transport velocity) and by the parameters of the kinetic-filtration model used for the simulations. Arrival times at the water table decrease with an increasing colloid size because of smaller diffusion, increased straining, and higher transport velocities. The importance of diffusion as a retardation mechanism increases with a decreasing colloid size, but appears to be minimal in large colloids.

Characterization of flow and transport processes within the unsaturated zone of Yucca Mountain, Nevada, under current and future climates

This paper presents a large-scale modeling study characterizing fluid flow and tracer transport in the unsaturated zone of Yucca Mountain, Nevada, a potential repository site for storing high-level radioactive waste. The study has been conducted using a three-dimensional numerical model, which incorporates a wide variety of field data and takes into account the coupled processes of flow and transport in the highly heterogeneous, unsaturated fractured porous rock. The modeling approach is based on a dual-continuum formulation of coupled multiphase fluid and tracer transport through fractured porous rock. Various scenarios of current and future climate conditions and their effects on the unsaturated zone are evaluated to aid in the assessment of the proposed repository's system performance using different conceptual models. These models are calibrated against field-measured data. Model-predicted flow and transport processes under current and future climates are discussed.