Effects of inherent spatial variability of rock properties on the thermo-hydro-mechanical responses of a high-level radioactive waste repository

Abstract

In France, a deep geological disposal for high level radioactive waste (HLW) and intermediate-level long-lived radioactive waste (ILW-LL) called Cigéo is planned to be constructed in a deep Callovo-Oxfordian claystone (COx) formation. The heat released from the HLW packages leads to a temperature increase in the host rock, which causes a pore pressure increase essentially due to the difference between the thermal expansion of the pore water and the solid skeleton. The low permeability of the COx and its relative rigidity inhibit the discharge of the induced pressure build-up. Moreover, thermal loading also induces thermo-mechanical stresses within the formation. The French National Radioactive Waste Management Agency (Andra) participated in the DECOVALEX-2019 project (Task E) to model the Thermo-Hydro-Mechanical (THM) responses of a benchmark exercise inspired of a HLW repository based on the French concept. This paper aims at studying the THM behavior of the HLW repository by taking into account the inherent spatial variability of the properties of the host rock through thermo-poro-elastic numerical simulations. Only the most influential parameters were selected to perform the spatial variability analysis by using Monte Carlo simulations coupled with the Random Finite Element Method. The selection of the most influential parameters was based on two parametric sensitivity analyses consisting in computing the Sobol indices, which determine the contribution of each input parameter to the overall model output variance. The temporal evolution of the Terzaghi effective stress in the far field is defined as indicator of the influence of each parameter on the THM response.