Numerical modelling of spatially and temporally distributed on-fault induced seismicity: implication for seismic hazards

International Journal of Coal Science & Technology

Research

Hani Mitri¹  Atsushi Sainoki²  Adam Karl Schwartzkopff³  Lishuai Jiang⁴ 

[1]Department of Mining and Materials Engineering, McGill University, Montreal, Canada

[2]Faculty of Advanced Science and Technology, Kumamoto University, Kumamoto, Japan

[3]School of Engineering, The University of Western Australia, Perth, Australia

[4]State Key Laboratory of Mining Disaster Prevention and Control, Shandong University of Science and Technology, Shandong, China

关键词: Induced seismicity;    Fracture network;    Fault damage zone;    Crack tensor theory;   

DOI  :  10.1007/s40789-022-00560-7

 received in 2022-03-29, accepted in 2022-12-30,  发布年份 2022

来源: Springer

【 摘 要 】

Induced seismicity is strongly related to various engineering projects that cause anthropogenic in-situ stress change at a great depth. Hence, there is a need to estimate and mitigate the associated risks. In the past, various simulation methods have been developed and applied to induced seismicity analysis, but there is still a fundamental difference between simulation results and field observations in terms of the spatial distribution of seismic events and its frequency. The present study aims to develop a method to simulate spatially distributed on-fault seismicity whilst reproducing a complex stress state in the fault zone. Hence, an equivalent continuum model is constructed, based on a discrete fracture network within a fault damage zone, by employing the crack tensor theory. A fault core is simulated at the center of the model as a discontinuous plane. Using the model, a heterogeneous stress state with stress anomalies in the fault zone is first simulated by applying tractions on the model outer boundaries. Subsequently, the effective normal stress on the fault plane is decreased in a stepwise manner to induce slip. The simulation result is validated in terms of the b-value and other seismic source parameters, hence demonstrating that the model can reproduce spatially and temporally distributed on-fault seismicity. Further analysis on the parameters shows the variation of frequency-magnitude distribution before the occurrence of large seismic events. This variation is found to be consistent with field observations, thus suggesting the potential use of this simulation method in evaluating the risk for seismic hazards in various engineering projects.

【 授权许可】

CC BY   
© The Author(s) 2023

【 图 表 】

41116_2022_35_Article_IEq580.gif

41116_2022_35_Article_IEq578.gif

41116_2022_35_Article_IEq576.gif

41116_2022_35_Article_IEq573.gif

Fig. 42

41116_2022_35_Article_IEq563.gif

41116_2022_35_Article_IEq561.gif

41116_2022_35_Article_IEq560.gif

41116_2022_35_Article_IEq558.gif

41116_2022_35_Article_IEq550.gif

41116_2022_35_Article_IEq537.gif

41116_2022_35_Article_IEq530.gif

Fig. 2

41116_2022_35_Article_IEq518.gif

41116_2022_35_Article_IEq503.gif

41116_2022_35_Article_IEq496.gif

Fig. 21

41116_2022_35_Article_IEq473.gif