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.

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.