Time-lapse seismic tomography of an underground mining zone

Abstract In this study, we determined P-wave tomographic images of the Yongshaba deposit, an underground mining zone in Guizhou Province, China, by inverting arrival-time data of micro-seismic events and blasts recorded by a passive seismic array consisting of 28 sensors during January to April 2014 using an event location technique and a travel-time tomography method that can handle complex seismic discontinuities. The damping parameter used in the damped least-squares method is determined with the help of trade-off curve. To reveal internal P-wave velocity changes of the study area under severe mining activities, a time-lapse data partition scheme is used. To assess the human influence on the underground structure, we introduce a parameter to measure the relative velocity changes in different periods. The resolution of the tomographic images and the robustness of the obtained features are examined by conducting a series of checkerboard resolution tests and a restoring resolution test. Our tomographic results obtained from the entire data set reveal a prominent low-velocity (low-V) zone and many obvious high-velocity (high-V) zones. These features match well with the geological setting and the excavation plan carried out in the mine, according to in-site surveys. The low-V zone may reflect empty volumes, stress releases and rock breaking and cracking caused by mining activities, whereas the high-V zones are probably the consequences of local stress concentrations caused by regional stress redistribution. The time-lapse tomographic images may reveal the process of stress concentrations caused by rock breaking due to the continuing excavations in the entire observation period, indicating that the mining activities influenced the rock property and caused complex changes of the underground structure.

[1]  M. Jackson,et al.  Geochemical zoning of volcanic chains associated with Pacific hotspots , 2011 .

[2]  Tao Xu,et al.  Microseismicity Induced by Fault Activation During the Fracture Process of a Crown Pillar , 2015, Rock Mechanics and Rock Engineering.

[3]  Keiiti Aki,et al.  Determination of the three‐dimensional seismic structure of the lithosphere , 1977 .

[4]  H. Shiobara,et al.  South Pacific mantle plumes imaged by seismic observation on islands and seafloor , 2009 .

[5]  Akira Hasegawa,et al.  Tomographic imaging of P and S wave velocity structure beneath northeastern Japan , 1992 .

[6]  Clifford H. Thurber,et al.  A fast algorithm for two-point seismic ray tracing , 1987 .

[7]  S. Vinciguerra,et al.  Relating seismic velocities, thermal cracking and permeability in Mt. Etna and Iceland basalts , 2005 .

[8]  D. F. Scott,et al.  3-D tomographic imaging of anomalous conditions in a deep silver mine , 1995 .

[10]  J. W. Cole,et al.  Physical property relationships of the Rotokawa Andesite, a significant geothermal reservoir rock in the Taupo Volcanic Zone, New Zealand , 2014, Geothermal Energy.

[11]  G. Cho,et al.  Long-Wavelength Elastic Wave Propagation Across Naturally Fractured Rock Masses , 2014, Rock Mechanics and Rock Engineering.

[12]  Dapeng Zhao,et al.  Tomographic imaging of the Cascadia subduction zone: Constraints on the Juan de Fuca slab , 2015 .

[13]  D. Zhao.  Multiscale Seismic Tomography , 2015 .

[14]  Dazhao Song,et al.  Numerical simulation of rock-burst relief and prevention by water-jet cutting , 2014 .

[15]  Robert G. Clapp,et al.  Time‐lapse seismic noise correlation tomography at Valhall , 2014 .

[16]  Zhouchuan Huang,et al.  Mechanism of the 2011 Tohoku-oki earthquake (Mw 9.0) and tsunami: Insight from seismic tomography , 2013 .

[17]  D. F. Scott,et al.  Temporal imaging of mine-induced stress change using seismic tomography , 1997 .

[18]  Xin Liu,et al.  P and S wave tomography of Japan subduction zone from joint inversions of local and teleseismic travel times and surface-wave data , 2016 .

[19]  D. Zhao,et al.  Seismic evidence for a metastable olivine wedge in the subducting Pacific slab under Japan Sea , 2008 .

[20]  Roland Martin,et al.  Multiple Scattering of Elastic Waves by Subsurface Fractures and Cavities , 2006 .

[21]  M. Oda,et al.  Material property changes in granitic rock during long-term immersion in hot water , 1995 .

[22]  O. P. Mishra,et al.  Influence of fluids and magma on earthquakes: seismological evidence , 2002 .

[23]  T. Mukerji,et al.  In situ identification of high vertical stress areas in an underground coal mine panel using seismic refraction tomography , 2015 .

[24]  Frederik Tilmann,et al.  Coalescence microseismic mapping , 2013 .

[25]  Sijing Wang,et al.  Engineering geology, ground surface movement and fissures induced by underground mining in the Jinchuan Nickel Mine , 2004 .

[26]  Nong Zhang,et al.  Microseismic multi-parameter characteristics of rockburst hazard induced by hard roof fall and high stress concentration , 2015 .

[27]  M. Beiki,et al.  Application of genetic programming to predict the uniaxial compressive strength and elastic modulus of carbonate rocks , 2013 .

[28]  H. Kanamori,et al.  Tomography of the Source Area of the 1995 Kobe Earthquake: Evidence for Fluids at the Hypocenter? , 1996, Science.

[29]  Xin Liu,et al.  Depth-varying azimuthal anisotropy in the Tohoku subduction channel , 2017 .

[30]  Koichi Sassa,et al.  Seismic attenuation tomography and its application to rock mass evaluation , 1996 .

[31]  Sanford Ballard,et al.  Efficient and Accurate Calculation of Ray Theory Seismic Travel Time Through Variable Resolution 3D Earth Models , 2009 .

[32]  G. Pavlis,et al.  The mixed discrete‐continuous inverse problem: Application to the simultaneous determination of earthquake hypocenters and velocity structure , 1980 .

[33]  Xibing Li,et al.  Rockburst characteristics and numerical simulation based on a strain energy density index: A case study of a roadway in Linglong gold mine, China , 2017 .

[34]  Xibing Li,et al.  Fracture Evolution Around a Cavity in Brittle Rock Under Uniaxial Compression and Coupled Static–Dynamic Loads , 2018, Rock Mechanics and Rock Engineering.

[35]  Yixian Xu,et al.  P wave anisotropic tomography of the Alps , 2017 .

[36]  A. Malehmir,et al.  Seismic characterization of the Grängesberg iron deposit and its mining-induced structures, central Sweden , 2015 .

[37]  Clifford H. Thurber,et al.  Earthquake locations and three‐dimensional crustal structure in the Coyote Lake Area, central California , 1983 .

[38]  Remko Scharroo,et al.  Generic Mapping Tools: Improved Version Released , 2013 .

[39]  S. Ivansson,et al.  Seismic borehole tomography—Measurement system and field studies , 1986, Proceedings of the IEEE.

[40]  Karl Berteussen,et al.  A microseismic experiment in Abu Dhabi, United Arab Emirates: implications for carbonate reservoir monitoring , 2013, Arabian Journal of Geosciences.

[41]  Xibing Li,et al.  Seismic attenuation tomography of the source zone of the 2016 Kumamoto earthquake (M 7.3) , 2017 .

[42]  Wu Cai,et al.  Application of seismic velocity tomography in underground coal mines: A case study of Yima mining area, Henan, China , 2014 .

[43]  John E. Peterson,et al.  Applications of algebraic reconstruction techniques to crosshole seismic data , 1985 .

[44]  V. Oye,et al.  Automated microearthquake location using envelope stacking and robust global optimization , 2010 .

[45]  B. Biondi,et al.  Target-oriented joint least-squares migration/inversion of time-lapse seismic data sets , 2010 .

[46]  N. Umino,et al.  Tomographic Imaging outside a Seismic Network: Application to the Northeast Japan Arc , 2007 .

[47]  B. Kennett,et al.  Joint seismic tomography for bulk sound and shear wave speed in the Earth's mantle , 1998 .

[48]  S. Ivansson,et al.  Seismic borehole tomography—Theory and computational methods , 1986, Proceedings of the IEEE.

[49]  M. Monjezi,et al.  Prediction of the strength and elasticity modulus of granite through an expert artificial neural network , 2015, Arabian Journal of Geosciences.

[50]  S. Reneau,et al.  Fault interaction and along-strike variation in throw in the Pajarito fault system, Rio Grande rift, New Mexico , 2009 .

[51]  F. Waldhauser,et al.  Streaks, multiplets, and holes: High‐resolution spatio‐temporal behavior of Parkfield seismicity , 2004 .

[52]  M. Downey,et al.  Geotomography for the delineation of coal seam structure , 1987 .

[53]  Peng Yan,et al.  Dynamic response of rock mass induced by the transient release of in-situ stress , 2012 .

[54]  M. Ge.  Comparison of Least Squares and Absolute Value Methods in Ae/MS Source Location: A Case Study , 1997 .

[55]  Hiroshi Takenaka,et al.  A TWO-POINT, THREE-DIMENSIONAL SEISMIC RAY TRACING USING GENETIC ALGORITHMS , 1999 .

[56]  Influence of Dip and Velocity Heterogeneity on Reverse- and Normal-Faulting Rupture Dynamics and Near-Fault Ground Motions , 2007 .

[57]  Keiiti Aki,et al.  Determination of three‐dimensional velocity anomalies under a seismic array using first P arrival times from local earthquakes: 1. A homogeneous initial model , 1976 .

[58]  Xibing Li,et al.  Impact of blasting parameters on vibration signal spectrum: Determination and statistical evidence , 2015 .

[59]  Long-Jun Dong,et al.  Locating single-point sources from arrival times containing large picking errors (LPEs): the virtual field optimization method (VFOM) , 2016, Scientific Reports.

[60]  Hu He,et al.  Rockburst hazard determination by using computed tomography technology in deep workface , 2012 .

[61]  Linming Dou,et al.  Quantitative analysis of seismic velocity tomography in rock burst hazard assessment , 2015, Natural Hazards.

[62]  Xibing Li,et al.  Determination of the minimum thickness of crown pillar for safe exploitation of a subsea gold mine based on numerical modelling , 2013 .

[63]  D. Eberhart‐Phillips,et al.  Three-dimensional velocity structure in northern California Coast Ranges from inversion of local earthquake arrival times , 1986 .

[64]  Hiroshi Ogasawara,et al.  Nucleation process of an M2 earthquake in a deep gold mine in South Africa inferred from on‐fault foreshock activity , 2015 .

[65]  M. Toksöz,et al.  An improved method for hydrofracture induced microseismic event detection and phase picking , 2010 .

[66]  Kourosh Shahriar,et al.  Studying the stress redistribution around the longwall mining panel using passive seismic velocity tomography and geostatistical estimation , 2013, Arabian Journal of Geosciences.

[67]  K. Luxbacher,et al.  Three-dimensional time-lapse velocity tomography of an underground longwall panel , 2008 .

[68]  David N. Dewhurst,et al.  Stress anisotropy and velocity anisotropy in low porosity shale , 2011 .

[69]  J. Mori,et al.  High resolution seismic velocity structure around the Yamasaki fault zone of southwest Japan as revealed from travel-time tomography , 2013, Earth, Planets and Space.

[70]  Gene Simmons,et al.  Velocity of compressional waves in various minerals at pressures to 10 kilobars , 1964 .

[71]  D. Stromeyer,et al.  Mining-Induced Stress Transfer and Its Relation to a $$\text{M}_w$$Mw 1.9 Seismic Event in an Ultra-deep South African Gold Mine , 2015 .

[72]  A. Lurka.  Location of high seismic activity zones and seismic hazard assessment in Zabrze Bielszowice coal mine using passive tomography , 2008 .

[73]  P. Bois,et al.  Well-to-well seismic measurements , 1972 .

[74]  Stefan Buske,et al.  Seismic travel-time and attenuation tomography to characterize the excavation damaged zone and the surrounding rock mass of a newly excavated ramp and chamber , 2014 .

[75]  Hiroshi Ogasawara,et al.  Steady activity of microfractures on geological faults loaded by mining stress , 2015 .

[76]  M. Baan,et al.  Dynamic triggering of microseismicity in a mine setting , 2015 .

[77]  Dapeng Zhao,et al.  P and S wave attenuation tomography of the Japan subduction zone , 2017 .