Omid K. Mahabadi

Y-GUI : A graphical user interface and pre-processor for the combined finite-discrete element code, Y2D, incorporating material heterogeneity

Numerical modelling of a discontinuous medium has gained much popularity in recent decades. The combined finite-discrete element method (FEM/DEM) is a state-of-the-art numerical modelling technique pioneered in the mid-1990s. Y2D is a robust two-dimensional FEM/DEM research code developed by Munjiza in 2004. The major limitations of this code are (1) the lack of a graphical user interface (GUI) meaning that all pre-processing has to be made directly on an ASCII input file and (2) the inability of dealing with heterogeneous media. This contribution presents the first GUI and pre-processor, known as Y-GUI, developed for Y2D and the implementation of a new algorithm that allows for the use of heterogeneous materials. In the text all major FEM/DEM concepts are described, together with the main features available in the Y-GUI. The use of Y-GUI is presented in detail and some of its functionalities, including the heterogeneity module to be used to randomly assign materials to a mesh, are introduced. At the end of the manuscript, four case studies, including Brazilian tests of a homogeneous and a layered rock sample and a rock avalanche, are presented.

Influence of microscale heterogeneity and microstructure on the tensile behavior of crystalline rocks

This study investigates the influence of microscale heterogeneity and microcracks on the failure behavior and mechanical response of a crystalline rock. The thin section analysis for obtaining the microcrack density is presented. Using micro X-ray computed tomography (μCT) scanning of failed laboratory specimens, the influence of heterogeneity and, in particular, biotite grains on the brittle fracture of the specimens is discussed and various failure patterns are characterized. Three groups of numerical simulations are presented, which demonstrate the role of microcracks and the influence of μCT-based and stochastically generated phase distributions. The mechanical response, stress distribution, and fracturing process obtained by the numerical simulations are also discussed. The simulation results illustrate that heterogeneity and microcracks should be considered to accurately predict the tensile strength and failure behavior of the sample.

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.

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.