3D Modelling of Excavation-Induced Anisotropic Responses of Deep Drifts at the Meuse/Haute-Marne URL

Abstract

Understanding the excavation induced fractured zone (EFZ) around drifts is paramount in the context of the deep geological disposal for nuclear waste since fractures can introduce pathways for the migration of radionuclides.Drifts in the Meuse/Haute-Marne Underground Research Laboratory (URL) have been essentially excavated following the two main directions of major and minor horizontal stresses. Field observations on the two drifts GCS (parallel to major horizontal stress direction) and GED (parallel to minor horizontal stress direction) in the URL show anisotropic shapes of EFZ around drifts through both orientations and anisotropic convergences. These anisotropic responses resulted from the inherent and/or induced anisotropies of the host rock; as well as the anisotropic stress field.This study focuses on 3D numerical modelling of excavation induced anisotropic responses including shape and extent of EFZ, and short term convergences of drifts. The main assumption is that the failure of claystone material is due to fracturing along weakness planes (ubiquitous joints) and the failure of the rock matrix. The ubiquitous joint failure is represented by perfectly plastic models for both tensile and shear yield functions. Their orientation is determined from the stress state based on the fracture mechanics, which includes tensile, longitudinal splitting and shear (conjugate planes) cracks. The rock matrix is assumed to be elastoplastic with hardening, softening and residual behaviours. Confining pressure dependency for the post-peak behaviour with a brittle-ductile transition is taken into account for the rock matrix. The proposed model is implemented into a commercial numerical software FLAC 3D . The main features of the implemented model are shown by the simulation of laboratory triaxial compression tests, as well as field observation within the URL. In particular, comparisons between 3D simulations of GCS and GED drifts with in situ observations shows promising results, which demonstrates advances of present model advances with respect to existing models.

Numerical simulation of underground excavations in an indurated clay using non-local regularisation. Part 1: formulation and base case

“…Different approaches have been adopted for the simulation of excavations, including localised deformations, in the context of nuclear waste disposal. Perhaps, the most extended technique is the local second gradient model (van den Eijnden et al, 2017;Pardoen et al, 2015;Pardoen and Collin, 2017;Salehnia et al, 2015) although other approaches such as coupling between elastoplastic rock matrix and weakness planes (ubiquitous joints) (Souley et al, 2020); discrete (or semi-discrete) approaches such as the extended rigid block spring method (Yao et al, 2017), or the combined finite-discrete element method (Lisjak et al, 2014) have also been employed. This paper presents a detailed simulation, using a different approach, of the GCS drift, a tunnel excavated parallel to the σ H direction at the MHM URL main level.…”

Scalable block preconditioners for saturated thermo-hydro-mechanics problems

“…Each linear solve requires less than 9 outer iterations and takes roughly 30 s. Given the total number of 130 Newton's iterations, the simulation is achieved in 3800 s. If the Schur preconditioner is used, the total solution time reaches 4700 s, increasing by 24%. We emphasize that the use of the proposed preconditioner allows for a simulation whose size is incomparable to other simulations in the specific field of radioactive waste disposal [4,53,54].…”

Constitutive Model for Thermal Compaction of Clayey Geomaterials and Application to COx Claystone

“…where 𝜀 𝑠 𝑝 (the internal flow variable) is the plastic distortion (or shear plastic strain), varies between a minimum value  0 and a maximum value  m . 𝛽(𝜀 𝑠 𝑝 ) is the dilatancy rate (or coefficient) that its evolution is based on Souley et al (2022) and given by:…”