Numerical simulation of gas flow through unsaturated fractured rock at Yucca Mountain, Nevada

Abstract

Numerical analysis is used to identify the physical phenomena associated with barometrically driven gas (air and water vapor) flow through unsaturated fractured rock at Yucca Mountain, Nevada. Results from simple finite difference simulations indicate that for a fractured rock scenario, the maximum velocity of air out of an uncased 10 cm borehole is 0.002 m s{sub {minus}1}. An equivalent porous medium (EPM) model was incorporated into a multiphase, multicomponent simulator to test more complex conceptual models. Results indicate that for a typical June day, a diurnal pressure wave propagates about 160 m into the surrounding Tiva Canyon hydrogeologic unit. Dry air that enters the formation evaporates water around the borehole which reduces capillary pressure. Multiphase countercurrent flow develops in the vicinity of the hole; the gas phase flows into the formation while the liquid phase flows toward the borehole. The effect occurs within 0.5 m of the borehole. The amount of water vapor leaving the formation during 1 day is 900 cm{sup 3}. This is less than 0.1% of the total recharge into the formation, suggesting that the barometric effect may be insignificant in drying the unsaturated zone. However, gas phase velocities out of the borehole (3 m s{sup {minus}1}), indicating that observed flow rates from wells along the east flank of Yucca Mountain were able to be simulated with a barometric model.