Andrea Lisjak

Numerical Modelling of the Anisotropic Mechanical Behaviour of Opalinus Clay at the Laboratory-Scale Using FEM/DEM

The Opalinus Clay (OPA) is an argillaceous rock formation selected to host a deep geologic repository for high-level nuclear waste in Switzerland. It has been shown that the excavation damaged zone (EDZ) in this formation is heavily affected by the anisotropic mechanical response of the material related to the presence of bedding planes. In this context, the purpose of this study is twofold: (i) to illustrate the new developments that have been introduced into the combined finite-discrete element method (FEM/DEM) to model layered materials and (ii) to demonstrate the effectiveness of this new modelling approach in simulating the short-term mechanical response of OPA at the laboratory-scale. A transversely isotropic elastic constitutive law is implemented to account for the anisotropic elastic modulus, while a procedure to incorporate a distribution of preferentially oriented defects is devised to capture the anisotropic strength. Laboratory results of indirect tensile tests and uniaxial compression tests are used to calibrate the numerical model. Emergent strength and deformation properties, together with the simulated damage mechanisms, are shown to be in strong agreement with experimental observations. Subsequently, the calibrated model is validated by investigating the effect of confinement and the influence of the loading angle with respect to the specimen anisotropy. Simulated fracture patterns are discussed in the context of the theory of brittle rock failure and analyzed with reference to the EDZ formation mechanisms observed at the Mont Terri Underground Research Laboratory.

Influence of pre-existing discontinuities and bedding planes on hydraulic fracturing initiation

Pressure-driven fracturing, also known as hydraulic fracturing, is a process widely used for developing geothermal resources, extracting hydrocarbons from unconventional reservoirs such as tight sandstone and shale formations, as well as for preconditioning the rock-mass during deep mining operations. While the overall process of pressure-driven fracturing is well understood, a quantitative description of the process is difficult due to both geologic and mechanistic uncertainties. Among them, the simulation of fractures growing in a complex heterogeneous medium is associated with computational difficulties. Experimental evidence based on micro-seismic monitoring clearly demonstrates the important influence of rock mass fabric on hydraulic fracture development, and the interaction between fluid-driven fractures and pre-existing discontinuities. However, these components are not well accounted for by standard numerical approaches. Thus, the design of hydraulic fracturing operations continues to be based on simplified models whereby the rock mass is treated as a homogeneous continuum. The purpose of this paper is to present the preliminary results obtained using the combined finite-discrete element technology to study the interaction between fluid driven fractures and natural rock mass discontinuities.

Hybrid Finite-Discrete Element Simulation of the EDZ Formation and Mechanical Sealing Process Around a Microtunnel in Opalinus Clay

The analysis and prediction of the rock mass disturbance around underground excavations are critical components of the performance and safety assessment of deep geological repositories for nuclear waste. In the short term, an excavation damaged zone (EDZ) tends to develop due to the redistribution of stresses around the underground openings. The EDZ is associated with an increase in hydraulic conductivity of several orders of magnitude. In argillaceous rocks, sealing mechanisms ultimately lead to a partial reduction in the effective hydraulic conductivity of the EDZ with time. The goal of this study is to strengthen the understanding of the phenomena involved in the EDZ formation and sealing in Opalinus Clay, an indurated claystone currently being assessed as a host rock for a geological repository in Switzerland. To achieve this goal, hybrid finite-discrete element method (FDEM) simulations are performed. With its explicit consideration of fracturing processes, FDEM modeling is applied to the HG-A experiment, an in situ test carried out at the Mont Terri underground rock laboratory to investigate the hydro-mechanical response of a backfilled and sealed microtunnel. A quantitative simulation of the EDZ formation process around the microtunnel is first carried out, and the numerical results are compared with field observations. Then, the re-compression of the EDZ under the effect of a purely mechanical loading, capturing the increase of swelling pressure from the backfill onto the rock, is considered. The simulation results highlight distinctive rock failure kinematics due to the bedded structure of the rock mass. Also, fracture termination is simulated at the intersection with a pre-existing discontinuity, representing a fault plane oblique to the bedding orientation. Simulation of the EDZ re-compression indicates an overall reduction of the total fracture area as a function of the applied pressure, with locations of ineffective sealing associated with self-propping of fractures. These results are consistent with hydraulic testing data revealing a negative correlation between pressure values and an increase in the EDZ transmissivity.

Hydro-mechanical evolution of the EDZ as transport path for radionuclides and gas: insights from the Mont Terri rock laboratory (Switzerland)

The excavation damaged zone (EDZ) around the backfilled underground structures of a geological repository represents a release path for radionuclides, which needs to be addressed in the assessment of long-term safety. Additionally, the EDZ may form a highly efficient escape route for corrosion and degradation gases, thus limiting the gas overpressures in the backfilled repository structures. The efficiency of this release path depends not only on the shape and extent of the EDZ, but also on the self-sealing capacity of the host rock formation and the prevailing state conditions, such as in situ stresses and pore pressure. The hydro-mechanical and chemico-osmotic phenomena associated with the formation and temporal evolution of the EDZ are complex, thus precluding a detailed representation of the EDZ in conventional modelling tools for safety assessment. Therefore, simplified EDZ models, able to mimic the safety-relevant functional features of the EDZ in a traceable manner are required. In the framework of the Mont Terri Project, a versatile modelling approach has been developed for the simulation of flow and transport processes along the EDZ with the goal of capturing the evolution of hydraulic significance of the EDZ after closure of the backfilled underground structures. The approach draws on both empirical evidence and experimental data, collected in the niches and tunnels of the Mont Terri rock laboratory. The model was benchmarked with a data set from an in situ self-sealing experiment at the Mont Terri rock laboratory. This paper summarises the outcomes of the benchmark exercise that comprises relevant empirical evidence, experimental data bases and the conceptual framework for modelling the evolution of the hydraulic significance of the EDZ around a backfilled tunnel section during the entire re-saturation phase.

A pragmatic approach to abstract the excavation damaged zone around tunnels of a geological radioactive waste repository: application to the HG-A experiment in Mont Terri

Abstract The excavation damaged zone (EDZ) around the backfilled tunnels of a geological repository represents a possible release path for radionuclides, corrosion and degradation gases that needs to be adequately addressed by safety assessment (SA) modelling tools. The hydromechanical phenomena associated with the creation and temporal evolution of the EDZ are of high complexity, precluding detailed representations of the EDZ in conventional SA. Thus, simplified EDZ models mimicking the safety-relevant features of the EDZ are required. In this context, a heuristic modelling approach has been developed to represent the creation and evolution of the EDZ in an abstracted and simplified manner. The key features addressed are the stochastic character of the excavation-induced fracture network and the self-sealing processes associated with the re-saturation after backfilling of the tunnels. The approach has been applied to a range of generic repository settings to investigate the impact of repository depth and in situ conditions on the hydraulic significance of the EDZ after repository closure. The model has been benchmarked with a dataset from a self-sealing experiment at the Mont Terri underground rock laboratory (URL), demonstrating the ability of the approach to mimic the evolution of the hydraulic significance of the EDZ during the re-saturation phase.