A reaction-transport model for calcite precipitation and evaluation of infiltration fluxes in unsaturated fractured rock

The percolation flux in the unsaturated zone (UZ) is an important parameter addressed in site characterization and flow and transport modeling of the potential nuclear-waste repository at Yucca Mountain, NV, USA. The US Geological Survey (USGS) has documented hydrogenic calcite abundances in fractures and lithophysal cavities at Yucca Mountain to provide constraints on percolation fluxes in the UZ. The purpose of this study was to investigate the relationship between percolation flux and measured calcite abundances using reactive transport modeling. Our model considers the following essential factors affecting calcite precipitation: (1) infiltration, (2) the ambient geothermal gradient, (3) gaseous CO(2) diffusive transport and partitioning in liquid and gas phases, (4) fracture-matrix interaction for water flow and chemical constituents, and (5) water-rock interaction. Over a bounding range of 2-20 mm/year infiltration rate, the simulated calcite distributions capture the trend in calcite abundances measured in a deep borehole (WT-24) by the USGS. The calcite is found predominantly in fractures in the welded tuffs, which is also captured by the model simulations. Simulations showed that from about 2 to 6 mm/year, the amount of calcite precipitated in the welded Topopah Spring tuff is sensitive to the infiltration rate. This dependence decreases at higher infiltration rates owing to a modification of the geothermal gradient from the increased percolation flux. The model also confirms the conceptual model for higher percolation fluxes in the fractures compared to the matrix in the welded units, and the significant contribution of Ca from water-rock interaction. This study indicates that reactive transport modeling of calcite deposition can yield important constraints on the unsaturated zone infiltration-percolation flux and provide useful insight into processes such as fracture-matrix interaction as well as conditions and parameters controlling calcite deposition.     

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.     

Tritium and 36Cl as constraints on fast fracture flow and percolation flux in the unsaturated zone at Yucca Mountain

An analysis of tritium and 36Cl data collected at Yucca Mountain, Nevada suggests that fracture flow may occur at high velocities through the thick unsaturated zone. The mechanisms and extent of this "fast flow" in fractures at Yucca Mountain are investigated with data analysis, mixing models and several one-dimensional modeling scenarios. The model results and data analysis provide evidence substantiating the weeps model [Gauthier, J.H., Wilson, M.L., Lauffer, F.C., 1992. Proceedings of the Third Annual International High-level Radioactive Waste Management Conference, vol. 1, Las Vegas, NV. American Nuclear Society, La Grange Park, IL, pp. 891-989] and suggest that fast flow in fractures with minimal fracture-matrix interaction may comprise a substantial proportion of the total infiltration through Yucca Mountain. Mixing calculations suggest that bomb-pulse tritium measurements, in general, represent the tail end of travel times for thermonuclear-test-era (bomb-pulse) infiltration. The data analysis shows that bomb-pulse tritium and 36Cl measurements are correlated with discrete features such as horizontal fractures and areas where lateral flow may occur. The results presented here imply that fast flow in fractures may be ubiquitous at Yucca Mountain, occurring when transient infiltration (storms) generates flow in the connected fracture network.     

Flow and transport in unsaturated fractured rock: effects of multiscale heterogeneity of hydrogeologic properties

The heterogeneity of hydrogeologic properties at different scales may have different effects on flow and transport processes in a subsurface system. A model for the unsaturated zone of Yucca Mountain, Nevada, is developed to represent complex heterogeneity at two different scales: (1) layer scale corresponding to geologic layering and (2) local scale. The layer-scale hydrogeologic properties are obtained using inverse modeling, based on the available measurements collected from the Yucca Mountain site. Calibration results show a significant lateral and vertical variability in matrix and fracture properties. Hydrogeologic property distributions in a two-dimensional, vertical cross-section of the site are generated by combining the average layer-scale matrix and fracture properties with local-scale perturbations generated using a stochastic simulation method. The unsaturated water flow and conservative (nonsorbing) tracer transport through the cross-section are simulated for different sets of matrix and fracture property fields. Comparison of simulation results indicates that the local-scale heterogeneity of matrix and fracture properties has a considerable effect on unsaturated flow processes, leading to fast flow paths in fractures and the matrix. These paths shorten the travel time of a conservative tracer from the source (repository) horizon in the unsaturated zone to the water table for small fractions of total released tracer mass. As a result, the local-scale heterogeneity also has a noticeable effect on global tracer transport processes, characterized by an average breakthrough curve at the water table, especially at the early arrival time of tracer mass. However, the effect is not significant at the later time after 20% tracer mass reaches the water table. The simulation results also verify that matrix diffusion plays an important role in overall solute transport processes in the unsaturated zone at Yucca Mountain.     

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.     

Field tracer-transport tests in unsaturated fractured tuff.

This paper presents the results of a field investigation in the unsaturated, fractured welded tuff within the Exploratory Studies Facility (ESF) at Yucca Mountain, NV. This investigation included a series of tests during which tracer-laced water was released into a high-permeability zone within a horizontal injection borehole. The tracer concentration was monitored in the seepage collected in an excavated slot about 1.6 m below the borehole. Results showed significant variability in the hydrologic response of fractures and the matrix. Analyses of the breakthrough curves suggest that flow and transport pathways are dynamic, rather than fixed, and related to liquid-release rates. Under high release rates, fractures acted as the predominant flow pathways, with limited fracture-matrix interaction. Under low release rates, fracture flow was comparatively less dominant, with a noticeable contribution from matrix flow. Observations of tracer concentrations rebounding in seepage water, following an interruption of flow, provided evidence of mass exchange between the fast-flowing fractures and slow- or non-flowing regions. The tests also showed the applicability of fluorinated benzoate tracers in situations where multiple tracers of similar physical properties are warranted.     

Flow and transport in fractured tuff at Yucca Mountain: numerical experiments on fast preferential flow mechanisms

Recent discovery of bomb-related ³⁶Cl at depth in fractured tuff in the unsaturated zone at the Yucca Mountain candidate high-level waste (HLW) repository site has called into question the usual modeling assumptions based on the equivalent continuum model (ECM). A dual continuum model (DCM) for simulating transient flow and transport at Yucca Mountain is developed. In order to ensure properly converged flow solutions, which are used in the transport simulation, a new flow solution convergence criteria is derived. An extensive series of simulation studies is presented which indicates that rapid movement of solute through the fractures will not occur unless there are intense episodic infiltration events. Movement of solute in the environs of the repository is enhanced if the properties of the tuff layer at the repository horizon are modified from current best-estimate values. Due to a large advective–dispersive coupling between the matrix and fractures, the matrix acts as a major buffer which inhibits rapid transport along the fractures. Consequently, fast movement of solutes through the fractures to the repository depth can only be explained if the matrix–fracture coupling term is significantly reduced from a value that would be calculated on the basis of data currently available.     

Tritium and Cl as constraints on fast fracture flow and percolation flux in the unsaturated zone at Yucca Mountain

An analysis of tritium and ³⁶Cl data collected at Yucca Mountain, Nevada suggests that fracture flow may occur at high velocities through the thick unsaturated zone. The mechanisms and extent of this “fast flow” in fractures at Yucca Mountain are investigated with data analysis, mixing models and several one-dimensional modeling scenarios. The model results and data analysis provide evidence substantiating the weeps model [Gauthier, J.H., Wilson, M.L., Lauffer, F.C., 1992. Proceedings of the Third Annual International High-level Radioactive Waste Management Conference, vol. 1, Las Vegas, NV. American Nuclear Society, La Grange Park, IL, pp. 891–989] and suggest that fast flow in fractures with minimal fracture–matrix interaction may comprise a substantial proportion of the total infiltration through Yucca Mountain. Mixing calculations suggest that bomb-pulse tritium measurements, in general, represent the tail end of travel times for thermonuclear-test-era (bomb-pulse) infiltration. The data analysis shows that bomb-pulse tritium and ³⁶Cl measurements are correlated with discrete features such as horizontal fractures and areas where lateral flow may occur. The results presented here imply that fast flow in fractures may be ubiquitous at Yucca Mountain, occurring when transient infiltration (storms) generates flow in the connected fracture network.     

Modeling tracer diffusion in fractured and unfractured, unsaturated, porous media

A proposed tracer diffusion test for the Exploratory Shaft Facility at Yucca Mountain, NV, is modeled. For the proposed test, a solution containing conservative tracers will be introduced into a borehole in the geologic medium of interest. The tracers will diffuse and advect from the saturated source region into the unsaturated matrix in the surrounding tuff. After some time, the borehole is to be overcored, and tracer concentrations in the fluid will be measured in the core as a function of distance from emplacement. The data will be used to evaluate diffusive behavior and to derive effective diffusion coefficients for the tracers in the specific tuff. Numerical simulations are used to study the effects of effective diffusion coefficient, porosity, saturation, and fracturing on tracer transport. Results are reported for numerical simulations of tests in the Topopah Spring Member and the Tuff of Calico Hills, which have significantly different porosities and saturations. The simulations make the following predictions: The spread of tracer during the test will be sensitive to the effective diffusion coefficient of the tracer. Tracer will diffuse farther in the Topopah Spring Member than in the Tuff of Calico Hills because of the former's lower porosity and saturation. Tracer transport by advection into the Topopah Spring Member will be greater than that into the Tuff of Calico Hills because of capillary effects. While advection will be a significant mechanism for tracer penetration into the Topopah Spring tuff, it will be less significant for tracer penetration into the Calico Hills tuff. The proximity of a single vertical fracture to the source region determines its effects on tracer transport, especially if the fracture diverts fluid flowing from the source region into the matrix.     

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.     

Incorporating layer- and local-scale heterogeneities in numerical simulation of unsaturated flow and tracer transport

This study characterizes layer- and local-scale heterogeneities in hydraulic parameters (i.e., matrix permeability and porosity) and investigates the relative effect of layer- and local-scale heterogeneities on the uncertainty assessment of unsaturated flow and tracer transport in the unsaturated zone of Yucca Mountain, USA. The layer-scale heterogeneity is specific to hydrogeologic layers with layerwise properties, while the local-scale heterogeneity refers to the spatial variation of hydraulic properties within a layer. A Monte Carlo method is used to estimate mean, variance, and 5th, and 95th percentiles for the quantities of interest (e.g., matrix saturation and normalized cumulative mass arrival). Model simulations of unsaturated flow are evaluated by comparing the simulated and observed matrix saturations. Local-scale heterogeneity is examined by comparing the results of this study with those of the previous study that only considers layer-scale heterogeneity. We find that local-scale heterogeneity significantly increases predictive uncertainty in the percolation fluxes and tracer plumes, whereas the mean predictions are only slightly affected by the local-scale heterogeneity. The mean travel time of the conservative and reactive tracers to the water table in the early stage increases significantly due to the local-scale heterogeneity, while the influence of local-scale heterogeneity on travel time gradually decreases over time. Layer-scale heterogeneity is more important than local-scale heterogeneity for simulating overall tracer travel time, suggesting that it would be more cost-effective to reduce the layer-scale parameter uncertainty in order to reduce predictive uncertainty in tracer transport.     

Solute transport in crystalline rocks at Aspö--I: geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted Asp? Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175]. The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments. While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours-days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption K(d)s are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.     

Solute transport in crystalline rocks at Äspö — I: Geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted Äspö Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175].The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments.While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours–days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption K_ds are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.     

Modelling the migration of contaminants through variably saturated dual-porosity, dual-permeability chalk

In the Hesbaye region in Belgium, tracer tests performed in variably saturated fissured chalk rocks presented very contrasting results in terms of transit times, according to artificially controlled water recharge conditions prevailing during the experiments. Under intense recharge conditions, tracers migrated across the partially or fully saturated fissure network, at high velocity in accordance with the high hydraulic conductivity and low effective porosity (fracture porosity). At the same time, a portion of the tracer was temporarily retarded in the almost immobile water located in the matrix. Under natural infiltration conditions, the fissure network remained inactive. Tracers migrated downward through the matrix, at low velocity in relation with the low hydraulic conductivity and the large porosity of the matrix. Based on these observations, Brouyère et al. (2004a) [Brouyère, S., Dassargues, A., Hallet, V., 2004a. Migration of contaminants through the unsaturated zone overlying the Hesbaye chalky aquifer in Belgium: a field investigation, J. Contam. Hydrol., 72 (1-4), 135-164, doi: 10.1016/j.conhyd.2003.10.009] proposed a conceptual model in order to explain the migration of solutes in variably saturated, dual-porosity, dual-permeability chalk. Here, mathematical and numerical modelling of tracer and contaminant migration in variably saturated fissured chalk is presented, considering the aforementioned conceptual model. A new mathematical formulation is proposed to represent the unsaturated properties of the fissured chalk in a more dynamic and appropriate way. At the same time, the rock water content is partitioned between mobile and immobile water phases, as a function of the water saturation of the chalk rock. The groundwater flow and contaminant transport in the variably saturated chalk is solved using the control volume finite element method. Modelling the field tracer experiments performed in the variably saturated chalk shows the adequacy and usefulness of the new conceptual, mathematical and numerical model.     

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.     

Migration of a water pulse through fractured porous media

Contaminant transport from waste-disposal sites is strongly affected by the presence of fractures and the degree of fracture matrix interaction. Characterization of potential contaminant plumes at such sites is difficult, both experimentally and numerically. Simulations of water flow through fractured rock were performed to examine the penetration depth of a large pulse of water entering such a system. Construction water traced with lithium bromide was released during the excavation of a tunnel at Yucca Mountain, Nevada, which is located in an unsaturated fractured tuff formation. Modeling of construction-water migration is qualitatively compared with bromide-to-chloride ratio (Br/Cl) data for pore-water salts extracted from drillcores. The influences of local heterogeneities in the fracture network and variations in hydrogeologic parameters were examined by sensitivity analyses and Monte Carlo simulations. The simulation results are qualitatively consistent with the observed Br/Cl signals, although these data may only indicate a minimum penetration depth, and water may have migrated farther through the fracture network.     

Commentary: assessment of past infiltration fluxes through Yucca Mountain on the basis of the secondary mineral record-is it a viable methodology?

Two papers recently published in the Journal of Contaminant Hydrology by Marshall et al. [Marshall, B.D., Neymark, L.A., Peterman, Z.E., 2003. Estimation of past seepage volumes from calcite distribution in the Topopah Spring Tuff, Yucca Mountain, Nevada. J. Contam. Hydrol. 62-63, 237-247] and Xu et al. [Xu, T., Sonnenthal, E., Bodvarsson, G., 2003. A reaction-transport model for calcite precipitation and evaluation of infiltration fluxes in unsaturated fractured rock. J. Contam. Hydrol. 64, 113-127] attempt to assess past volumes of seepage and infiltration fluxes through the vadose zone of Yucca Mountain, Nevada, on the basis of the modeling of the spatial distribution of secondary calcite. In this commentary, we argue that the employed methodology is not viable. In addition, the thermal boundary conditions used in simulations do not correspond to the temperatures of the mineral forming fluids established on the basis of the fluid inclusion studies.     

Kinematic analysis of solute mass flows in rock fractures with spatially random parameters

Field data of physical properties in heterogeneous crystalline bedrock, like porosity and fracture aperture, is associated with uncertainty that can have a significant impact on the analysis of solute transport in rock fractures. Solutions to the central temporal moments of the residence time probability density function (PDF) are derived in a closed form for a solute Dirac pulse. The solutions are based on a model that takes into account advection along the fracture plane, diffusion into the rock matrix and sorption kinetics in the rock matrix. The most relevant rock properties including fracture aperture and several matrix properties as well as flow velocity are assumed to be spatially random along transport pathways. The mass transport is first solved in a general form along one-dimensional pathways, but the results can be extended to multi-dimensional flows simply by substituting the expected travel time for inert water parcels. Based on data obtained with rock samples taken at Asp? Hard Rock Laboratory in Sweden, the solutions indicate that the heterogeneity of the rock properties contributes to increasing significantly both the variance and the skewness of the residence time probability density function for a pulse travelling in a fracture. The Asp? data suggests that the bias introduced in the variance of the residence time PDF by neglecting the effect of heterogeneity of the rock properties on the radionuclide migration is very large for fractures thinner than a few tenths of a millimetre.     

Modeling of non-reactive solute transport in fractured clayey till during variable flow rate and time

Fractures and biopores can act as preferential flow paths in clay aquitards and may rapidly transmit contaminants into underlying aquifers. Reliable numerical models for assessment of groundwater contamination from such aquitards are needed for planning, regulatory and remediation purposes. In this investigation, high resolution preferential water-saturated flow and bromide transport data were used to evaluate the suitability of equivalent porous medium (EPM), dual porosity (DP) and discrete fracture/matrix diffusion (DFMD) numerical modeling approaches for assessment of flow and non-reactive solute transport in clayey till. The experimental data were obtained from four large undisturbed soil columns (taken from 1.5 to 3.5 m depth) in which biopores and channels along fractures controlled 96-99% of water-saturated flow. Simulating the transport data with the EPM effective porosity model (FRACTRAN in EPM mode) was not successful because calibrated effective porosity for the same column had to be varied up to 1 order of magnitude in order to simulate solute breakthrough for the applied flow rates between 11 and 49 mm/day. Attempts to simulate the same data with the DP models CXTFIT and MODFLOW/MT3D were also unsuccessful because fitted values for dispersion, mobile zone porosity, and mass transfer coefficient between mobile and immobile zones varied several orders of magnitude for the different flow rates, and because dispersion values were furthermore not physically realistic. Only the DFMD modeling approach (FRACTRAN in DFMD mode) was capable to simulate the observed changes in solute transport behavior during alternating flow rate without changing values of calibrated fracture spacing and fracture aperture to represent the macropores.     

Transport of radon gas into a tunnel at Yucca Mountain--estimating large-scale fractured tuff hydraulic properties and implications for the operation of the ventilation system

Radon gas concentrations have been monitored as part of the operation of a tunnel (the Exploratory Studies Facility-ESF) at Yucca Mountain to ensure worker safety. The objective of this study was to examine the potential use of the radon data to estimate large-scale formation properties of fractured tuffs. This objective was examined by developing a numerical model, based upon the characteristics of the ESF and the Topopah Spring welded (TSw) tuff unit, capable of predicting radon concentrations for prescribed ventilation conditions. The model was used to address two specific issues. First, it was used to estimate the permeability and porosity of the fractures in the TSw at the length scale of the ESF and extending tens of meters into the TSw, which surrounds the ESF. Second, the model was used to understand the mechanism leading to radon concentrations exceeding a specified level within the ESF. The mechanism controlling radon concentrations in the ESF is a function of atmospheric barometric fluctuations being propagated down the ESF along with ventilated air flow and the slight suction induced by the ventilation exhaust fans at the South Portal of the ESF. These pressure fluctuations are dampened in the TSw fracture continuum according to its permeability and porosity. Consequently, as the barometric pressure in the ESF drops rapidly, formation gases from the TSw are pulled into the ESF, resulting in an increase in radon concentrations. Model calibration to both radon concentrations measured in the ESF and gas-phase pressure fluctuations in the TSw yielded concurrent estimates of TSw fracture permeability and porosity of 1 x 10(-11) m2 and 0.00034, respectively. The calibrated model was then used as a design tool to predict the effect of adjusting the current ventilation-system operation strategy for reducing the probability of radon gas concentrations exceeding a specified level.