Investigating the Influence of Mechanical Anisotropy on the Fracturing Behaviour of Brittle Clay Shales with Application to Deep Geological Repositories

Clay shales are currently being assessed as possible host rock formations for the deep geological disposal of radioactive waste. However, one main concern is that the favourable long-term isolation properties of the intact rock mass could be negatively affected by the formation of an excavation damaged zone (EDZ) around the underground openings. The goal of this thesis was to develop and validate a mechanically-based model of the deformation and failure process of a clay shale, namely Opalinus Clay (OPA), with particular focus on the influence of anisotropy on the short-term response of circular tunnels. To achieve this goal, a hybrid continuum-discontinuum numerical approach was developed in combination with new field measurements from the Mont Terri underground research laboratory. The response of OPA during the excavation of a full-scale emplacement (FE) test tunnel was characterized by geodetic monitoring of wall displacements, radial extensometers and longitudinal inclinometers. The deformation measurements indicated strong directionality induced by the combined effect of in-situ stress field and bedding planes striking parallel to the tunnel axis, with the most severe deformation occurring in the direction approximately perpendicular to the material layering. Computer simulations were conducted using a newly-extended combined finite-discrete element method (FEM/DEM), a numerical technique which allows the explicit simulation of brittle fracturing and associated seismicity. The numerical experimentation firstly focused on the laboratory-scale analysis of failure processes (e.g., acoustic activity) in brittle rocks, and on the role of strength and modulus anisotropy in the failure behaviour of OPA in tension and compression. The fracturing behaviour of unsupported circular excavations in laminated rock masses was then analyzed under different in-situ stress conditions. Lastly, the proposed modelling methodology was applied to the aforementioned FE tunnel to obtain unique insights into the EDZ formation process around emplacement tunnels for nuclear waste. The calibrated numerical model suggested delamination along bedding planes and subsequent extensional fracturing as key mechanisms of the damage process potentially leading to buckling and spalling phenomena. Overall, the research findings are already having a direct impact on the preliminary design process of an underground repository for nuclear waste in Switzerland.