ABSTRACT
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.
ABSTRACT
Water bodies provide essential ecosystem services linked to morphometric features that might differ between natural lakes and reservoirs. We use the HydroLAKES global dataset to quantitatively compare large (area > 1 km2) reservoirs and natural lakes in terms of scaling exponents between morphometric measures (volume, area, shore length). These exponents are further compared to those expected from geometrical assumptions and constraints. Lakes cover a larger range of volumes for the same range of surface areas than reservoirs, and have a larger volume-area scaling exponent. The volume-area scaling exponent for reservoirs (but not natural lakes) and the area-shore length exponent for all water bodies follow the predictions for self-affine surfaces. Land cover and terrain influence the scaling relations more for lakes than for reservoirs. These morphometric differences may be used to model the impact of reservoirs and lakes on hydrological processes and associated ecosystem services at regional to global scales.
