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].

Parametric Sensitivity to Capillary Entry Pressure in Two-Phase Water–Gas Flow Models: Deep Geologic Disposal of Radioactive Waste

Gas-entry pressure impact on the evaluation of hydrogen migration at different scales of a deep geological disposal of radioactive waste

Two phase partially miscible flow and transport modeling in porous media: application to gas migration in a nuclear waste repository

Modelling and Numerical Simulation of Gas Migration in a Nuclear Waste Repository

Numerical modeling of gas injection induced cracking with unsaturated hydromechanical process in the context of radioactive waste disposal

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

Diffusion and Gas Flow Dynamics in Partially Saturated Smectites

Multiphase flow and reactive transport benchmark for radioactive waste disposal

Coupled multiphase flow and viscoelastic mechanics modeling of gas injection in a compacted bentonite buffer

A Gas Flow Model for Layered Landfills with Vertical Extraction Wells.

Role of capillary pressure on carbon dioxide sequestration: mathematical analysis

The application of high-resolution 3D seismic data to model the distribution of mechanical and hydrogeological properties of a potential host rock for the deep storage of radioactive waste in France

Evaluation of Gas Migration and Rock Damage Characteristics for Underground Nuclear Waste Storage Based on a Coupled Model

Soft Interface Instability and Gas Flow Channeling in Low-Permeability Deformable Media

Isothermal Dual-Phase Flow Modeling to Assess the Impact of Gas Collection on Geotechnical and Hydraulic Performance of Landfills

Reactive transport modeling for assessing coupled hydrogeochemical processes at interfaces in deep geological repositories: from the laboratory to the real world

Pore‐Scale Modeling of Carbon Dioxide and Hydrogen Transport During Geologic Gas Storage

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

A dual porosity model of high-pressure gas flow for geoenergy applications

Seepage Characteristics of Mixed-Wettability Porous Media on the Phase-Field Model

Numerical Assessments of Flow and Advective Transport Uncertainty for Performance Measures of Radioactive Waste Geological Disposal in Fractured Rocks