On Water Flow in Hot Fractured Rock -- A Sensitivity Study on theImpact of Fracture-Matrix Heat Transfer

Dual-continuum models have been widely used in modeling flowand transport in fractured porous rocks. Among many other applications,dual-continuum approaches were utilized in predictive models of thethermal-hydrological conditions near emplacement tunnels (drifts) atYucca Mountain, Nevada, the proposed site for a radioactive wasterepository in the U.S. In unsaturated formations such as those at YuccaMountain, the magnitude of mass and heat exchange between the twocontinua fracture network and matrix is largely dependent on the flowcharacteristics in the fractures, because channelized finger-type flowstrongly reduces the interface area between the matrix surfaces and theflowing liquid. This effect may have important implications, for example,during the time period that the fractured rock near the repository driftswould be heated above the boiling point of water. Depending on themagnitude of heat transfer from the matrix, water percolating down thefractures will either boil off in the hot rock region above drifts or maypenetrate all the way to the drift walls and possibly seep into the opencavities. In this paper, we describe a sensitivity analysis using avariety of approaches to treat fracture-matrix interaction in athree-dimensional dual-continuum setting. Our simulation example is alaboratory heater experiment described in the literature that providesevidence of rapid water flow in fractures, leading to drift seepagedespite above-boiling conditions in the adjacent fractured rock. Theexperimental finding can only be reproduced when the interface area forheat transfer between the matrix and fracture continua is reduced toaccount for flow channeling