Stress waves, known as acoustic emissions (AEs), are released by localized inelastic deformation events during the progressive failure of brittle rocks. Although several numerical models have been developed to simulate the deformation and damage processes of rocks, such as non-linear stress-strain behaviour and localization of failure, only a limited number have been capable of providing quantitative information regarding the associated seismicity. Moreover, the majority of these studies have adopted a pseudo-static approach based on elastic strain energy dissipation that completely disregards elastodynamic effects. This paper describes a new AE modelling technique based on the combined finite-discrete element method (FEM/DEM), a numerical tool that simulates material failure by explicitly considering fracture nucleation and propagation in the modelling domain. Given the explicit time integration scheme of the solver, stress wave propagation and the effect of radiated seismic energy can be directly captured. Quasi-dynamic seismic information is extracted from a FEM/DEM model with a newly developed algorithm based on the monitoring of internal variables (e.g. relative displacements and kinetic energy) in proximity to propagating cracks. The AE of a wing crack propagation model based on this algorithm are cross-analysed by traveltime inversion and energy estimation from seismic recordings. Results indicate a good correlation of AE initiation times and locations, and scaling of energies, independently calculated with the two methods. Finally, the modelling technique is validated by simulating a laboratory compression test on a granite sample. The micromechanical parameters of the heterogeneous model are first calibrated to reproduce the macroscopic stress-strain response measured during standard laboratory tests. Subsequently, AE frequency-magnitude statistics, spatial clustering of source locations and the evolution of AE rate are investigated. The distribution of event magnitude tends to decay as power law while the spatial distribution of sources exhibits a fractal character, in agreement with experimental observations. Moreover, the model can capture the decrease of seismic b value associated with the macrorupture of the rock sample and the transition of AE spatial distribution from diffuse, in the pre-peak stage, to strongly localized at the peak and post-peak stages, as reported in a number of published laboratory studies. In future studies, the validated FEM/DEM-AE modelling technique will be used to obtain further insights into the micromechanics of rock failure with potential applications ranging from laboratory-scale microcracking to engineering-scale processes (e.g. excavations within mines, tunnels and caverns, petroleum and geothermal reservoirs) to tectonic earthquakes triggering.
a b s t r a c tThe Opalinus Clay formation is currently being investigated as a potential host rock for the deep geological disposal of radioactive waste in Switzerland. Recently, a test tunnel was excavated at the Mont Terri underground rock laboratory (URL) as part of a long-term research project (''Full-scale Emplacement (FE) experiment'') aimed at studying the thermo-hydro-mechanical (THM) effects induced by the presence of an underground repository. The objective of this paper is twofold. Firstly, the results of the rock mass monitoring programme carried out during the construction of the 3 m diameter, 50 m long FE tunnel are presented, with particular focus on the short-term deformation response. The deformation measurements, including geodetic monitoring of tunnel wall displacements, radial extensometers and longitudinal inclinometers, indicate a strong directionality in the excavation response. Secondly, the deformational behaviour observed in the field is analyzed using a hybrid finite-discrete element (FDEM) analysis to obtain further insights into the formation of the excavation damaged zone (EDZ). The FDEM simulation using the Y-Geo code is calibrated based on the average short-term response observed in the field. Deformation and strength anisotropy are captured using a transversely isotropic, linear elastic constitutive law and cohesive elements with orientation-dependent strength parameters. Overall, a good agreement is obtained between convergences measured in the field and numerical results. The simulated EDZ formation process highlights the importance of bedding planes in controlling the failure mechanisms around the underground opening. Specifically, failure initiates due to shearing of bedding planes critically oriented with respect to the compressive circumferential stress induced around the tunnel. Slippageinduced rock mass deconfinement then promotes extensional fracturing in the direction perpendicular to the bedding orientation. The simulated fracture pattern is consistent with previous experimental evidence from the Mont Terri URL.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.