2002 Denver Annual Meeting (October 27-30, 2002)

Paper No. 14
Presentation Time: 11:30 AM


HOUSEWORTH, J.E.1, FINSTERLE, S.A.1 and BODVARSSON, G.S.2, (1)Earth Sciences Division, Lawrence Berkeley National Lab, 1 Cyclotron Rd, Berkeley, CA 94720, (2)Earth Sciences Division, Lawrence Berkeley National Lab, 1 Cyclotron Road, Berkeley, CA 94720, jehouseworth@lbl.gov

The current concept for high-level radioactive waste disposal at Yucca Mountain is for 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 because of the 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 because of 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 are 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.