Experimental and numerical analysis of flow through a natural rough fracture subject to normal loading.

Fractured crystalline rocks have been chosen or are under consideration by several countries as host rock formations for deep geological repositories for spent nuclear fuel. In such geological formations, flow and solute transport are mostly controlled by a network of connected natural fractures, each of them being characterised by internal heterogeneity, also denoted as roughness. Fractures are, in turn, subject to variable load caused by various factors, such as the presence of thick ice sheets formed during glaciation periods. Understanding how coupled hydro-mechanical (HM) processes affect flow and transport at the scale of a single natural fracture is crucial for a robust parameterisation of large-scale discrete fracture network models, which are not only used for nuclear waste disposal applications but are also of interest to problems related to geothermics, oil and gas production or groundwater remediation. In this work, we analyse and model an HM experiment carried out in a single natural fracture and use the results of both, the experimental and the modelling work, to get insights into fundamental questions such as the applicability of local cubic law or the effect of normal load on channeling. The initial fracture aperture was obtained from laser scanning of the two fracture surfaces and an equivalent initial aperture was then defined by moving the two fracture surfaces together and comparing the results obtained using a Navier-Stokes based computational fluid dynamics (CFD) model with the experimental flowrate obtained for unloaded conditions. The mechanical effect of the different loading stages was simulated using a high-resolution contact model. The different computed fracture apertures were then used to run groundwater flow simulations using a modified Reynolds equation. The results show that, without correction, local cubic law largely overestimates flowrates. Instead, we show that by explicitly acknowledging the difference between the mechanical aperture and the hydraulic aperture and setting the latter equal to 1/5 of the former, cubic law provides a very reasonable approximation of the experimental flowrates over the entire loading cycle. A positive correlation between fluid flow channeling and normal load is also found.

Transport of oxygen into granitic rocks: Role of physical and mineralogical heterogeneity.

The rock matrix of granites is expected to be an important buffer against the dispersion of contaminants, e.g. radionuclides, and against the ingress of oxygenated glacial meltwater. The influence of matrix heterogeneity on O2 diffusive transport is assessed here by means of numerical experiments based on a micro-Discrete Fracture Network (micro-DFN) representation of the diffusion-available pore space along with random realisations of idealized biotite grains, to simulate the heterogeneous nature of granitic rocks. A homogeneous-based analytical solution is also presented and used to assess possible deviations of the numerical experiments from the assumption of homogeneity. The analytical solution is also used to test upscaled values of mineral surface area. The numerical experiments show that the matrix behaves as a composite system, with the coexistence of fast and slow diffusive pathways. This behavior is more evident at low Damköhler numbers. Our interpretation of the numerical experiments points out the importance to properly characterise the heterogeneity of the rock matrix.

Microtomography-based Inter-Granular Network for the simulation of radionuclide diffusion and sorption in a granitic rock.

Field investigation studies, conducted in the context of safety analyses of deep geological repositories for nuclear waste, have pointed out that in fractured crystalline rocks sorbing radionuclides can diffuse surprisingly long distances deep into the intact rock matrix; i.e. much longer distances than those predicted by reactive transport models based on a homogeneous description of the properties of the rock matrix. Here, we focus on cesium diffusion and use detailed micro characterisation data, based on micro computed tomography, along with a grain-scale Inter-Granular Network model, to offer a plausible explanation for the anomalously long cesium penetration profiles observed in these in-situ experiments. The sparse distribution of chemically reactive grains (i.e. grains belonging to sorbing mineral phases) is shown to have a strong control on the diffusive patterns of sorbing radionuclides. The computed penetration profiles of cesium agree well with an analytical model based on two parallel diffusive pathways. This agreement, along with visual inspection of the spatial distribution of cesium concentration, indicates that for sorbing radionuclides the medium indeed behaves as a composite system, with most of the mass being retained close to the injection boundary and a non-negligible part diffusing faster along preferential diffusive pathways.

Continuum-based DFN-consistent numerical framework for the simulation of oxygen infiltration into fractured crystalline rocks.

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.

A new method for the interpretation of pumping tests in leaky aquifers.

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.