Abstract:In the field of radioactive waste confinement, the question of gas transfers in clay formations is a crucial issue [1]. A certain amount of gas, such as Hydrogen may be generated during the exploitation phase in the nearfield by the deterioration of the metal components of the system. Since the host medium is characterised by a very low permeability, the mechanisms of gas transport by advection and diffusion within the groundwater remain insufficient to evacuate the gas generated and a free gas phase is formed. If the gas pressure keeps increasing and reaches the minimum principal stress acting on the rock, micro-fractures, known as preferential gas pathways could develop through the rock mass [2], and affect the clay barrier integrity. There is a growing body of experimental evidences [3, 4] that separation planes such as bedding planes or pre-existing fractures and heterogeneities in clay-rich materials represent preferred weaknesses for the process of opening discrete gas-filled pathways. Capturing the related transport mechanisms therefore requires to go from macroscopic to microscopic scale. In this context, the adopted numerical approach would ideally take into account the effect of each constituent of the material microstructure on the macroscopic gas flow. Yet, direct modelling of the entire microstructure using small-scale models is usually not possible due to the huge computational costs it would require at the scale of a repository. Conversely, indirect modelling of the behaviour of all the micro-constituents by collective closed-form macroscopic constitutive equations using large-scale models [5] has also limits in terms of assumptions formulation and parameters identification. Hence, the use of a multi-scale approach that models the micro-scale effects explicitely on their specific length scale and couples their homogenized effects to the macro-scale is proposed in the present work. Based on a periodicity assumption of the microstructure, the physical and geometrical properties of the microstructure are embedded on a Representative Element Volume (REV) which contains a detailed model of the microstructure constituents, i.e. the pore network, the bedding planes and the bridging planes, while eluding a complex description of the whole microstructure over the domain. From this local description of the material, the macroscopic response to loading can be derived using numerical homogenisation techniques [6]. The multi-scale model proposed here has been implemented in the finite element code Lagamine, and validated against a well established and documented macro-scale THM coupled model [7]. The size and structure of the REV comes from experimental data acquired from scanning images in clay material, like for instance the Boom Clay formation which is envisaged as a potential host rock for a deep geological disposal in Belgium. Such scans reveal the opening of large aperture with a repeated distance in-between, which make it possible to extrapolate a physical idealisation of the microstructure (Figure 1a), built with one of these fractures corresponding to the bedding planes (red) and the pore network substituted by an assembly of tubes (green). Narrow-aperture fissures, known as the bridging planes (blue) were also identified to contribute to the flow normal to the predominant direction of fissure [4], by connecting the natural discontinuities which were initially closed. These micro-scale constituents are defined in such a way as to satisfy the conditions of pore size distribution, macroporosity and intrinsic permeability. For the mechanical behaviour, a simple linear elastic constitutive law was introduced for the tube network while a non-linear stress-strain relationship was used for the interfaces defining the bedding and bridging planes. Concerning the hydraulic behaviour, gas diffusive fluxes were considered and modelled by a Fick’s law while a channel flow based on the Navier-Stokes equations was used for the multiphase flow in the different constituents. The model has been subsequently applied to simulate a gas injection test parallel and perpendicular to the bedding of initially saturated samples of Boom Clay [4]. This analysis provides a rather good agreement with the experimental results in terms of injection and outflow pressure response, outflow volume and average axial strain along the sample height. In addition, it allows to simulated the creation of a preferential flow pathway along the sample axis (Figure 1b, top), which serves as basis to numerically reproduce the development of random pathway through the sample in plane strain state (Figure 1b, bottom), and aims to improve the mechanistic understanding of the gas transport processes at play in clayey barriers [8].

Diffusion Experiments in Callovo-Oxfordian Clay from the Meuse/Haute-Marne URL, France. Experimental Setup and Data Analyses

Numerical Model of A Long-Term in Situ Diffusion and Retention (dr) Experiment in Opalinus Clay

DIGIT: An In Situ Experiment for Studying the Diffusion of Water and Solutes under Thermal Gradient in the Toarcian Clay Rock at the Tournemire Underground Research Laboratory: Part 1—Goals, Scoping Calculations, Installation and First Results under Unheated Conditions

A laboratory investigation of thermally induced pore pressures in the Callovo-Oxfordian Claystone

In-situ X-ray fluorescence to investigate iodide diffusion in opalinus clay: Demonstration of a novel experimental approach

Anion Diffusion in Compacted Clays by Pore‐Scale Simulation and Experiments

Solute transport properties of compacted Ca-bentonite used in FEBEX project

Characterization of heterogeneities in the sandy facies of the Opalinus Clay (Mont Terri underground rock laboratory, Switzerland)

RETRACTED ARTICLE: Experimental study and numerical prediction of HTO and 36 Cl − diffusion in radioactive waste at Téguline Clay

Treatment of boundary conditions in through-diffusion: A case study of (85)Sr(2+) diffusion in compacted illite

@Modelling of failure and fracture development of the Callovo-Oxfordian claystone during an in-situ heating experiment associated with geological disposal of high-level radioactive waste

Influence of salinity gradients on the diffusion of water and ionic species in dual porosity clay samples

Overview of radionuclide migration experiments in the HADES Underground Research Facility at Mol (Belgium)

Predictive Modeling of a Simple Field Matrix Diffusion Experiment Addressing Radionuclide Transport in Fractured Rock. is It So Straightforward?

A single-site reactive transport model of Cs + for the in situ diffusion and retention (DR) experiment

multi-scale model to study gas transport processes in clay materials

Studying radionuclide migration on different scales: the complementary roles of laboratory and in situ experiments

Modeling of a radionuclide transport experiment in granitic rock matrix at the grimsel test site (Switzerland). The role of advection

Investigation of Se(IV) Diffusion in Compacted Tamusu Clay by Capillary Method

The Influence of Sedimentary Heterogeneity on the Diffusion of Radionuclides in the Sandy Facies of Opalinus Clay at the Field Scale

Light reflection visualization to determine solute diffusion into clays