Search

Natural fractures are characterized by high internal heterogeneity. This internal variability is the cause of flow channeling, which in turn leads to contaminant transport taking place primarily along the high‐velocity channels. Mass exchange between the high‐velocity channels and the low‐velocity zones has the potential to enhance contaminant retention, due to solute diffusion into the low‐velocity zones and subsequent exposure to additional surface area for diffusion into the bordering rock matrix. Here, we derive a random walk particle tracking method for heterogeneous fractures, which includes an additional term to account for the aperture gradient. The method takes into account advection, diffusion in the fracture and matrix diffusion. The developed numerical framework is applied to assess the effect of low‐velocity zones in rough self‐affine fractures. The results show that diffusion into low‐velocity zones has a visible but modest impact on contaminant retention. The magnitude of this impact does not change considerably, regardless of whether diffusion into the rock matrix is considered in the model, and increases for a decreasing average Péclet number of the fracture. Key Points: Natural fractures are highly heterogeneous, comprising flowing channels and lower‐velocity zonesWe study the effects that diffusion into low‐velocity zones has in contaminant transport through rough fracturesAccounting for diffusion in low‐velocity zones has a relatively modest impact on contaminant retention [ABSTRACT FROM AUTHOR]

The REPRO-TDE test was performed at a depth of about 400 m in the ONKALO underground research facility in Finland. Synthetic groundwater containing radionuclide tracers [tritiated water tracer (HTO), 36Cl, 22Na, 133Ba, and 134Cs] was circulated for about 4 years in a packed-off interval of the injection borehole. Tracer activities were additionally monitored in two observation boreholes. The test was the subject of a modeling exercise by the SKB GroundWater Flow and Transport of Solutes Task Force. Eleven teams participated in the exercise, using different model concepts and approaches. Predictive model calculations were based on laboratory-based information concerning porosities, diffusion coefficients, and sorption partition coefficients. After the experimental results were made available, the teams were able to revise their models to reproduce the observations. General conclusions from these back-analysis calculations include the need for reduced effective diffusion coefficients for 36Cl compared to those applicable to HTO (anion exclusion), the need to implement weaker sorption for 22Na compared to results from laboratory batch sorption experiments, and the observation of large differences between the theoretical initial concentrations for the strongly sorbing 133Ba and 134Cs, and the first measured values a few hours after tracer injection. Different teams applied different concepts, concerning mainly the implementation of isotropic versus anisotropic diffusion, or the possible existence of borehole disturbed zones around the different boreholes. The role of microstructure was also addressed in two of the models. [ABSTRACT FROM AUTHOR]

We use a realistic fracture along with pore‐scale reactive transport simulations to investigate the impact of fracture filling minerals and water‐mineral reactions on in‐plane groundwater flow patterns and on contaminant uptake by matrix diffusion. The results of the simulations show that sparsely distributed fracture filling calcite patches provide a significant buffering capacity against the infiltration of acidic water. The advance of the geochemical perturbation is characterized by channeling, which is partly smeared out when low Peclet numbers are used. Mineral patches have a significant impact on groundwater flow patterns as shown by the computed distributions of groundwater travel time, which are characterized by long tails. The ability of the fracture to retain contaminants by uptake due to matrix diffusion is investigated by means of two metrics: the Eulerian quantity denoted as F‐factor and the Lagrangian variable known as transport resistance. F‐factor is shown to provide an average single estimate of the bulk retention capacity of the fracture whereas transport resistance is characterized by a distribution, whose very long tail is related to streamlines forced to travel through very narrow throats. A reduction of up to ∼70% in the mean travel time and ∼85% in the mean transport resistance is observed as a consequence of geochemical dissolution reactions, which also lead to an overall reduction of their variability. Key Points: Fracture internal variability and infilling minerals lead to channelingSparsely distributed calcite mineral patches provide a significant buffering capacity against the infiltration of acidic waterMineral dissolution leads to a general decrease of contaminant retention by matrix diffusion [ABSTRACT FROM AUTHOR]

In situ experiments and field‐scale characterization studies have pointed out that, in fractured crystalline media, groundwater flow is highly channelized. This implies that, at the scale of a single fracture, only part of the fracture surface area is in contact with flowing water, while the rest of in‐plane water is essentially stagnant and can be accessed by solutes via molecular diffusion. Despite their importance for contaminant retention, to date, there are no numerical or analytical approaches that could be used to assess the implication of stagnant water zones on solute transport in realistic large‐scale Discrete Fracture Network‐based models. Here, we present an efficient and flexible algorithm for the simulation of transport in fractured rock with diffusion into stagnant water and rock matrix. The algorithm is a generalization of a previously developed numerical framework for time domain particle tracking in sparsely fractured rock. The key of the generalization is that total time in fracture (τf) is first evaluated using a Monte Carlo sampling and then a second sampling is performed conditioned on τf. The algorithm has been successfully validated against existing independent solutions and the implication of diffusion into stagnant water and secondary diffusion into the matrix has been assessed for a realistic modeling scenario. The results show that, due to diffusion into stagnant water, contaminants are more strongly retarded. This increased retention is more significant for sorbing species, as a larger number of sorption sites is accessible. A high sensitivity to the flowing channel/stagnant water zone geometry has also been observed. Key Points: In fractured crystalline rocks flow is highly channelizedDiffusion into stagnant water and secondary diffusion into the matrix provide additional retentionWe present an efficient algorithm for the simulation of transport with diffusion into stagnant water and rock matrix [ABSTRACT FROM AUTHOR]

Antonellini et al. (2017, https://doi.org/10.1002/2016WR019864) presented a detailed analysis of alteration halos observed in outcrops of a fractured rock system. We point out here the need for a critical discussion on the low value of dissolved oxygen concentration used in the boundary water of their analytical models. Moreover, we argue that the model used by Antonellini et al. (2017) to evaluate the time for fracture filling is not consistent with the hypotheses used by the same authors to analyze the alteration halos. We propose here an alternative explanation along with an alternative model for fracture filling by precipitation of hydroxide minerals. Unlike the model used by Antonellini et al. (2017), in the model presented here, the time for fracture filling, defined as the time when the relative fracture aperture is equal to a given minimum threshold, does not depend on the initial fracture aperture. Key Points: Interpretation of alteration halos is significantly sensitive to oxygen concentration in the infiltrating waterIn the absence of organic‐rich sediments, oxygen content in the infiltrating water should be in equilibrium with atmospheric conditionsAn alternative model for the evaluation of the time for fracture filling is proposed [ABSTRACT FROM AUTHOR]

We present an enhanced continuum-based approach for the modelling of groundwater flow coupled with reactive transport in crystalline fractured rocks. In the proposed formulation, flow, transport and geochemical parameters are represented onto a numerical grid using Discrete Fracture Network (DFN) derived parameters. The geochemical reactions are further constrained by field observations of mineral distribution. To illustrate how the approach can be used to include physical and geochemical complexities into reactive transport calculations, we have analysed the potential ingress of oxygenated glacial-meltwater in a heterogeneous fractured rock using the Forsmark site (Sweden) as an example. The results of high-performance reactive transport calculations show that, after a quick oxygen penetration, steady state conditions are attained where abiotic reactions (i.e. the dissolution of chlorite and the homogeneous oxidation of aqueous iron(II) ions) counterbalance advective oxygen fluxes. The results show that most of the chlorite becomes depleted in the highly conductive deformation zones where higher mineral surface areas are available for reactions. [ABSTRACT FROM AUTHOR]

In the framework of safety assessment studies for geological disposal, large scale reactive transport models are powerful inter-disciplinary tools aiming at supporting regulatory decision making as well as providing input to repository engineering activities. Important aspects of these kinds of models are their often very large temporal and spatial modelling scales and the need to integrate different non-linear processes (e.g., mineral dissolution and precipitation, adsorption and desorption, microbial reactions and redox transformations). It turns out that these types of models may be computationally highly demanding. In this work, we present a Lagrangian-based framework, denoted as FASTREACT, that aims at solving multi-component-reactive transport problems with a computationally efficient approach allowing complex modelling problems to be solved in large spatial and temporal scales. The tool has been applied to simulate radionuclide migration in a synthetic heterogeneous transmissivity field and the results have been successfully compared with those obtained using a standard Eulerian approach. Finally, the same geochemical model has been coupled to an ensemble of realistic three-dimensional transport pathways to simulate the migration of a set of radionuclides from a hypothetical repository for spent nuclear fuel to the surface. The results of this modelling exercise, which includes key processes such as the exchange of mass between the conductive fractures and the matrix, show that FASTREACT can efficiently solve large-scale reactive transport models. [ABSTRACT FROM AUTHOR]

Waste rock piles are an outcome of open pit and underground mining operations. Unprocessed low-grade rock is disposed of in piles from several meters high to 100 m+ high. Waste rock piles may still contain sufficient concentrations of metals to be a potential source of pollution. The evaluation of the potential risk involves properly characterizing flow through these piles under unsaturated conditions. The main flow characteristic of the piles is the presence of a large range of grain and pore sizes. Based on data from an instrumented rock pile located in Saskatchewan Canada, unsaturated flow through the pile is modeled as a linear system after separating a fast and a slow component. Water reaching the base of the pile is monitored by 16 contiguous zero-tension lysimeters. The fast component, flowing through macropores, is assumed to be released instantaneously, while the slow component is simulated using a linear-reservoir model that assumes the presence of an interconnected porous matrix. An empirical transfer function (TF) is computed as the ratio of the spectra of signals between the output (basal outflow) and the input (rainfall time series). Determination of a parametric transfer function model provides information on the characteristic time of water storage in the matrix and on the fraction of the water within each subsection of the experimental pile that is channeled through the macropores. An analysis of the output signal at different support scales is performed, indicating the nonlinearity of the macropore fraction scaling processes. [ABSTRACT FROM AUTHOR]

Hydrologists routinely analyze pumping test data using conventional interpretation methods that are based on the assumption of homogeneity and that, consequently, yield single estimates of representative flow parameters. However, natural subsurface formations are intrinsically heterogeneous, and hence, the flow parameters influencing the drawdown vary as the cone of depression expands in time. In this paper a novel procedure for the analysis of pumping tests in heterogeneous confined aquifers is developed. We assume that a given heterogeneous aquifer can be represented by a homogeneous system whose flow parameters evolve in time as the pumping test progresses. At any point in time, the interpreted flow parameters are estimated using the ratio of the drawdown and its derivative observed at that particular time. The procedure is repeated for all times, yielding time-dependent estimates of transmissivity Ti( t) and storativity, Si( t). Based on the analysis of the sensitivity of drawdown to inhomogeneities in the T field, the time-dependent interpreted transmissivity values are found to be a good estimate of Tg( r), the geometric mean of the transmissivity values encompassed within a progressively increasing radius r from the well. The procedure is illustrated for Gaussian heterogeneous fields with ln( T) variances up to a value of 2. The impact of the separation distance between the pumping well and observation point on data interpretation is discussed. The results show that information about the spatial variability of the transmissivity field can be inferred from time-drawdown data collected at a single observation point. [ABSTRACT FROM AUTHOR]

Pumping tests are routinely interpreted from the analysis of drawdown data and their derivatives. These interpretations result in a small number of apparent parameter values which lump the underlying heterogeneous structure of the aquifer. Key questions in such interpretations are (1) what is the physical meaning of those lumped parameters and (2) whether it is possible to infer some information about the spatial variability of the hydraulic parameters. The system analyzed in this paper consists of an aquifer separated from a second recharging aquifer by means of an aquitard. The natural log transforms of the transmissivity, ln T, and the vertical conductance of the aquitard, ln C, are modeled as two independent second-order stationary spatial random functions (SRFs). The Monte Carlo approach is used to simulate the time-dependent drawdown at a suite of observation points for different values of the statistical parameters defining the SRFs. Drawdown data at each observation point are independently used to estimate hydraulic parameters using three existing methods: (1) the inflection-point method, (2) curve-fitting, and (3) the double inflection-point method. The resulting estimated parameters are shown to be space dependent and vary with the interpretation method since each method gives different emphasis to different parts of the time-drawdown data. Moreover, the heterogeneity in the pumped aquifer or the aquitard influences the estimates in distinct manners. Finally, we show that, by combining the parameter estimates obtained from the different analysis procedures, information about the heterogeneity of the leaky aquifer system may be inferred. [ABSTRACT FROM AUTHOR]

A novel methodology for the interpretation of pumping tests in leaky aquifer systems, referred to as the double inflection point (DIP) method, is presented. The method is based on the analysis of the first and second derivatives of the drawdown with respect to log time for the estimation of the flow parameters. Like commonly used analysis procedures, such as the type-curve approach developed by Walton (1962) and the inflection point method developed by Hantush (1956) , the mathematical development of the DIP method is based on the assumption of homogeneity of the leaky aquifer layers. However, contrary to the two methods developed by Hantush and Walton, the new method does not need any fitting process. In homogeneous media, the two classic methods and the one proposed here provide exact results for transmissivity, storativity, and leakage factor when aquifer storage is neglected and the recharging aquifer is unperturbed. The real advantage of the DIP method comes when applying all methods independently to a test in a heterogeneous aquifer, where each method yields parameter values that are weighted differently, and thus each method provides different information about the heterogeneity distribution. Therefore, the methods are complementary and not competitive. In particular, the combination of the DIP method and Hantush method is shown to lead to the identification of contrasts between the local transmissivity in the vicinity of the well and the equivalent transmissivity of the perturbed aquifer volume. [ABSTRACT FROM AUTHOR]

The article offers brief information on the Steel Connections eToolKit digital design instrument for steelwork connections from the Southern African Institute of Steel Construction (SAISC).