Numerical analysis of thermo-hydro-mechanical (THM) processes in the clay based material

Clay formations are investigated worldwide as potential host rock for the deep geological disposal of high-level radioactive waste (HLW). Usually bentonite is preferred as the buffer and backfill material in the disposal system. In the disposal of HLW, heat emission is one of the most important issues as it can generate a series of complex thermo-hydro-mechanical (THM) processes in the surrounding materials and thus change the material properties. In the context of safety assessment, it is important to understand the thermally induced THM interactions and the associated change in material properties. In this work, the thermally induced coupled THM behaviours in the clay host rock and in the bentonite buffer as well as the corresponding coupling effects among the relevant material properties are numerically analysed. A coupled non-isothermal Richards flow mechanical model and a non-isothermal multiphase flow model were developed based on the scientific computer codes OpenGeoSys (OGS). Heat transfer in the porous media is governed by thermal conduction and advective flow of the pore fluids. Within the hydraulic processes, evaporation, vapour diffusion, and the unsaturated flow field are considered. Darcy’s law is used to describe the advective flux of gas and liquid phases. The relative permeability of each phase is considered. The elastic deformation process is modelled by the generalized Hooke’s law complemented with additional strain caused by swelling/shrinkage behaviour and by temperature change. In this study, special attention has been paid to the analysis of the thermally induced changes in material properties. The strong mechanical and hydraulic anisotropic properties of clay rock are described by a transversely isotropic mechanical model and by a transversely isotropic permeability tensor, respectively. The thermal anisotropy is described by adoption of the bedding-orientation-dependent thermal conductivity. The dependency of the thermal conductivity on the degree of water saturation, the dependency of the thermal effects on the water retention behaviour, and the dependency of the effects of the pore pressure variation on the permeability and the anisotropic swelling/shrinkage behaviour have been intensively analysed and the corresponding numerical models to consider those coupling effects have been developed. The developed numerical model has been applied to simulate the laboratory and in situ heating experiments on the bentonite and clay rock at different scales. Firstly the laboratory heating experiment on Callovo-Oxfordian Clay (COX) and the laboratory long-term heating and hydration experiment on MX80 pellets were simulated. Based on the knowledge from the numerical analysis of the laboratory experiments, a 1:2 scale in situ heating experiment of an integrated system of the bentonite engineered barrier system (EBS) in the Opalinus Clay host rock was simulated. All the relevant operation phases were considered in the modelling. Besides, the modelling was extended to 50 years after the heat shut-down with the aim of predicting the long-term behaviours. Additionally, variation calculations were carried out to investigate the effects of the storage capacity of the Opalinus Clay on the thermally induced hydraulic response. In the long-term modelling, the effects of different saturated water permeabilities of buffer material on the resaturation process were analysed. Based on the current researches and model developments, the observed THM behaviours of the bentonite buffer and the clay rock, that is, the measured evolution of temperature, pore pressure, humidity, swelling pressure, and so on in the laboratory and in situ experiments can be reproduced and interpreted well. It is proved that by using both a non-isothermal multiphase flow model and a non-isothermal Richards flow model combined with the corresponding thermal and mechanical models, the major THM behaviours can be captured. It is validated that the developed model is able to simulate the relevant coupled THM behaviours of clayey material under the well-defined laboratory conditions as well as under the complex natural disposal conditions