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.     

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.     

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.     

Field investigation into unsaturated flow and transport in a fault: model analyses

Results of a fault test performed in the unsaturated zone of Yucca Mountain, Nevada, were analyzed using a three-dimensional numerical model. The fault was explicitly represented as a discrete feature and the surrounding rock was treated as a dual-continuum (fracture-matrix) system. Model calibration against seepage and water-travel-velocity data suggests that lithophysal cavities connected to fractures can considerably enhance the effective fracture porosity and therefore retard water flow in fractures. Comparisons between simulation results and tracer concentration data also indicate that matrix diffusion is an important mechanism for solute transport in unsaturated fractured rock. We found that an increased fault-matrix and fracture-matrix interface areas were needed to match the observed tracer data, which is consistent with previous studies. The study results suggest that the current site-scale model for the unsaturated zone of Yucca Mountain may underestimate radionuclide transport time within the unsaturated zone, because an increased fracture-matrix interface area and the increased effective fracture porosity arising from lithophysal cavities are not considered in the current site-scale model.     

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.     

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.     

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.     

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.     

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.     

Analysis of a mesoscale infiltration and water seepage test in unsaturated fractured rock: Spatial variabilities and discrete fracture patterns

A mesoscale (21 m in flow distance) infiltration and seepage test was recently conducted in a deep, unsaturated fractured rock system at the crossover point of two underground tunnels. Water was released from a 3 mx4 m infiltration plot on the floor of an alcove in the upper tunnel, and seepage was collected from the ceiling of a niche in the lower tunnel. Significant temporal and (particularly) spatial variabilities were observed in both measured infiltration and seepage rates. To analyze the test results, a three-dimensional unsaturated flow model was used. A column-based scheme was developed to capture heterogeneous hydraulic properties reflected by these spatial variabilities observed. Fracture permeability and van Genuchten alpha parameter [van Genuchten, M.T., 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 44, 892-898] were calibrated for each rock column in the upper and lower hydrogeologic units in the test bed. The calibrated fracture properties for the infiltration and seepage zone enabled a good match between simulated and measured (spatially varying) seepage rates. The numerical model was also able to capture the general trend of the highly transient seepage processes through a discrete fracture network. The calibrated properties and measured infiltration/seepage rates were further compared with mapped discrete fracture patterns at the top and bottom boundaries. The measured infiltration rates and calibrated fracture permeability of the upper unit were found to be partially controlled by the fracture patterns on the infiltration plot (as indicated by their positive correlations with fracture density). However, no correlation could be established between measured seepage rates and density of fractures mapped on the niche ceiling. This lack of correlation indicates the complexity of (preferential) unsaturated flow within the discrete fracture network. This also indicates that continuum-based modeling of unsaturated flow in fractured rock at mesoscale or a larger scale is not necessarily conditional explicitly on discrete fracture patterns.     

Comparison of two micro-analytical methods for detecting the spatial distribution of sorbed Pu on geologic materials

Subsurface transport of groundwater contaminants is greatly influenced by chemical speciation, precipitation and sorption processes at the mineral-water interface. The retardation of contaminants is often greatest at boundaries between minerals and in fractures and pore spaces. The investigation of the spatial distribution of sorbed contaminants along these boundaries requires micro-analytical techniques. The sorption of dissolved Pu(V) on a natural zeolitic tuff from Yucca Mountain (NV, USA) was examined using microautoradiography (MAR), X-ray diffraction (XRD), electron microprobe (EM) techniques, and synchrotron-based micro-X-ray fluorescence (micro-SXRF). The tuff contained a heterogeneous distribution of zeolites and trace quantities of smectites, Fe oxides (hematite), and Mn oxides (rancieite), which are present as fracture fill and pore space materials. Micro-SXRF studies showed that Pu is mostly associated with bodies of smectite plus Mn oxides, which were typically elevated in Ce, Ga, Nb, Pb, Y, Ca, Ti, and Zn. Sorbed Pu was not associated with Fe-rich bodies, which were enriched in Cl and Rb. Results of the MAR studies were complementary to that of the micro-SXRF studies in that Pu was associated with similar elements in the tuff. Indirect detection of Pu by EM or micro-SXRF (by analyzing Ag developed on the MAR photoemulsion) was a more sensitive method for detecting lower levels of sorbed Pu than the direct detection of sorbed Pu via micro-SXRF in the absence of the photoemulsion.     

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.     

Characterization of fracture aperture field heterogeneity by electrical resistance measurement

We use electrical resistance measurements to characterize the aperture field in a rough fracture. This is done by performing displacement experiments using two miscible fluids of different electrical resistivity and monitoring the time variation of the overall fracture resistance. Two fractures have been used: their complementary rough walls are identical but have different relative shear displacements which create "channel" or "barrier" structures in the aperture field, respectively parallel or perpendicular to the mean flow velocity U(→). In the "channel" geometry, the resistance displays an initial linear variation followed by a tail part which reflects the velocity contrast between slow and fast flow channels. In the "barrier" geometry, a change in the slope between two linear zones suggests the existence of domains of different characteristic aperture along the fracture. These variations are well reproduced analytically and numerically using simple flow models. For each geometry, we present then a data inversion procedure that allows one to extract the key features of the heterogeneity from the resistance measurement.     

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.     

Simulation and analysis of solute transport in 2D fracture/pipe networks: the SOLFRAC program

The Time Domain Random Walk (TDRW) method has been recently developed by Delay and Bodin [Delay, F. and Bodin, J., 2001. Time domain random walk method to simulate transport by advection-dispersion and matrix diffusion in fracture networks. Geophys. Res. Lett., 28(21): 4051-4054.] and Bodin et al. [Bodin, J., Porel, G. and Delay, F., 2003c. Simulation of solute transport in discrete fracture networks using the time domain random walk method. Earth Planet. Sci. Lett., 6566: 1-8.] for simulating solute transport in discrete fracture networks. It is assumed that the fracture network can reasonably be represented by a network of interconnected one-dimensional pipes (i.e. flow channels). Processes accounted for are: (1) advection and hydrodynamic dispersion in the channels, (2) matrix diffusion, (3) diffusion into stagnant zones within the fracture planes, (4) sorption reactions onto the fracture walls and in the matrix, (5) linear decay, and (6) mass sharing at fracture intersections. The TDRW method is handy and very efficient in terms of computation costs since it allows for the one-step calculation of the particle residence time in each bond of the network. This method has been programmed in C++, and efforts have been made to develop an efficient and user-friendly software, called SOLFRAC. This program is freely downloadable at the URL (labo.univ-poitiers.fr/hydrasa/intranet/telechargement.htm). It calculates solute transport into 2D pipe networks, while considering different types of injections and different concepts of local dispersion within each flow channel. Post-simulation analyses are also available, such as the mean velocity or the macroscopic dispersion at the scale of the entire network. The program may be used to evaluate how a given transport mechanism influences the macroscopic transport behaviour of fracture networks. It may also be used, as is the case, e.g., with analytical solutions, to interpret laboratory or field tracer test experiments performed in single fractures.     

Mass removal of chlorinated ethenes from rough-walled fractures using permanganate

In situ chemical oxidation (ISCO) employing permanganate is an emerging technology that has been successful at enhancing mass removal from DNAPL source zones in unconsolidated media at the pilot-scale. The focus of this study was to evaluate the applicability of flushing a permanganate solution across two single vertical fractures in a laboratory environment to remove free phase DNAPL. The fracture experiments were designed to represent a portion of a larger fractured aquifer system impacted by a near-surface DNAPL spill over a shallow fractured rock aquifer. Each fracture was characterized by hydraulic and tracer tests, and the aperture field for one of the fractures was mapped using a co-ordinate measurement machine. Following DNAPL emplacement, a series of water and permanganate flushes were performed. To support observations from the fracture experiments, a set of batch experiments was conducted. The data from both fracture experiments showed that the post-oxidation effluent concentration was not impacted by the oxidant flush; however, changes in the aperture distribution, flow field, and flow rate were observed. These changes resulted in a significant decrease to the mass loading from the fractures, and were attributed to the build-up of oxidation by-products (manganese oxides and carbon dioxide) within the fracture which was corroborated by the batch experiment data and visual examination of the walls of one fracture. These results provide insight into the potential impact that a permanganate solution and oxidation by-products can have on the aperture distribution within a fracture and on DNAPL mass transfer rates. A permanganate flush or injection completed within a fractured rock aquifer may lead to the development of an insoluble product adjacent to the DNAPL which results in the reduction or complete elimination of advective regions near the DNAPL and reduces mass transfer rates. This outcome would have significant implications on the plume generating potential of the remaining DNAPL.     

Hydraulic anisotropy characterization of pneumatic-fractured sediments using azimuthal self potential gradient

The pneumatic fracturing technique is used to enhance the permeability and porosity of tight unconsolidated soils (e.g. clays), thereby improving the effectiveness of remediation treatments. Azimuthal self potential gradient (ASPG) surveys were performed on a compacted, unconsolidated clay block in order to evaluate their potential to delineate contaminant migration pathways in a mechanically-induced fracture network. Azimuthal resistivity (ARS) measurements were also made for comparative purposes. Following similar procedures to those used in the field, compressed kaolinite sediments were pneumatically fractured and the resulting fracture geometry characterized from strike analysis of visible fractures combined with strike data from optical borehole televiewer (BHTV) imaging. We subsequently injected a simulated treatment (electrolyte/dye) into the fractures. Both ASPG and ARS data exhibit anisotropic geoelectric signatures resulting from the fracturing. Self potentials observed during injection of electrolyte are consistent with electrokinetic theory and previous laboratory results on a fracture block model. Visual (polar plot) analysis and linear regression of cross plots show ASPG lobes are correlated with azimuths of high fracture strike density, evidence that the ASPG anisotropy is a proxy measure of hydraulic anisotropy created by the pneumatic fracturing. However, ARS data are uncorrelated with fracture strike maxima and resistivity anisotropy is probably dominated by enhanced surface conduction along azimuths of weak ‘starter paths’ formed from pulverization of the clay and increases in interfacial surface area. We find the magnitude of electrokinetic SP scales with the applied N₂ gas pressure gradient (ΔPN₂) for any particular hydraulically-active fracture set and that the positive lobe of the ASPG anomaly indicates the flow direction within the fracture network. These findings demonstrate the use of ASPG in characterizing the effectiveness of (1) pneumatic fracturing and (2) defining likely flow directions of remedial treatments in unconsolidated sediments and rock.     

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.     

Determination of the flow-wetted surface in fractured media

Diffusion and sorption in the rock matrix are important retardation mechanisms for radionuclide transport in fractured media. For the conditions existing in a deep repository in crystalline rock, interaction with the rock matrix is controlled by the water flowrate in the fractures and the surface area in contact with the flowing water (the so-called "flow-wetted surface" (FWS)). The flow-wetted surface may be determined from the frequency of open fractures intersecting a borehole. The choice of packer distance used in these hydraulic measurements is crucial, however, since several open fractures may be found in one packer interval. The use of a packer distance that is too large may result in a considerable underestimation of the flow-wetted surface. This is especially important in zones with a high frequency of open fractures (fracture zones) where a small packer distance is a fundamental requirement. A large volume of hydraulic data has been compiled in Sweden from measurements using quite small packer distances. Over the last decade, the most common packer distance used for the hydraulic tests has been 3 m, although some new measurements using a shorter packer distance have also been performed. In several cases, the resolution of these measurements has been less than 0.5 m. All these data have been analysed in detail. From these data, the flow-wetted surface has been calculated and compared with the flow-wetted surface estimated in earlier studies. The results show the importance of using a small packer distance for carrying out borehole transmissivity measurements.