Study on the quantitative definition of blasting damage zone scope based on wave velocity field inversion
FU Shuaiyang1，2，LI Haibo1，2，WU Di1，2，WANG Ben1，2，LIU Liwang3，LI Xiaofeng1，2
(1. State Key Laboratory of Geomechanics and Geotechnical Engineering，Institute of Rock and Soil Mechanics，Chinese Academy of Sciences，Wuhan，Hubei 430071，China；2. University of Chinese Academy of Sciences，Beijing 100049，China；3. School of Mechanical Engineering，Nanjing University of Science and Technology，Nanjing，Jiangsu 210094，China)
Abstract：Determining the scope of the rock mass blasting damage zone is important for evaluating the stability of underground engineering surrounding rock under blasting and controlling bedrock damage under excavation. A method based on wave velocity field inversion was proposed to determine the rock blasting damage zone scope quantitatively. This method was based on multistencils fast marching methods(MSFM) and simultaneous iterative reconstructive technique(SIRT). The effectiveness of the proposed method was verified by comparing the inversed wave velocity field with the known one. Blasting and wave velocity measurement simulations were carried out on the self-developed finite-discrete element solver OpenFDEM. The waveform data were processed using a self-programmed post-processing program，and thus the post-blasting wave velocity field was reconstructed. Based on the relationship between the changing rate of wave velocity and the degree of rock mass damage in the current code，the areas with different wave velocity changing rates were compared with the areas of cracks. The scopes of the blasting damage zone with different damage degrees were quantified using the ratio of the equivalent radius to the blast hole radius. 10% and 15% of the wave velocity changing rate were used as the thresholds for defining slightly and completely damage of the rock mass，respectively. The ratio of the equivalent radius of the slightly damage zone and the damage zone to the radius of the blast hole is about 11.86 and 7.2，respectively. As the blast load increases，the shape of the damage zones with different damage degree becomes irregular and the equivalent radius increases. The study results are expected to provide a reference for determining the scope of the blasting damage zone in the field.
傅帅旸1，2，李海波1，2，吴 迪1，2，王 犇1，2，刘黎旺3，李晓锋1，2. 基于波速场反演的爆破损伤区范围界定研究[J]. 岩石力学与工程学报, 2024, 43(S1): 3257-3266.
FU Shuaiyang1，2，LI Haibo1，2，WU Di1，2，WANG Ben1，2，LIU Liwang3，LI Xiaofeng1，2. Study on the quantitative definition of blasting damage zone scope based on wave velocity field inversion. , 2024, 43(S1): 3257-3266.
[1] 夏 祥，李海波，李俊如，等. 岭澳核电站二期工程基岩爆破安全阈值分析[J]. 岩土力学，2008，29(11)：2 946–2 951.(XIA Xiang，LI Haibo，LI Junru，et al. Research on vibration safety threshold for rock under blasting excavation[J]. Rock and Soil Mechanics，2008，29(11)：2 946–2 951.(in Chinese))
[2] YANG J H，YAO C，JIANG Q H，et al. 2D numerical analysis of rock damage induced by dynamic in-situ stress redistribution and blast loading in underground blasting excavation[J]. Tunnelling and Underground Space Technology，2017，70：221–232.
[3] LI H B，XIA X，LI J C，et al. Rock damage control in bedrock blasting excavation for a nuclear power plant[J]. International Journal of Rock Mechanics and Mining Sciences，2011，48(2)：210–218.
[4] 中华人民共和国行业标准编写组. SL47—2020 水工建筑物岩石地基开挖施工技术规范[S]. 北京：中国水利水电出版社，2021.(The Professional Standards Compilation Group of the People?s Republic of China. SL47—2020 Technical specification for excavation construction of rock-foundation of hydraulic structures[S]. Beijing：China Water Power Press，2021.(in Chinese))
[5] ZENG Y Q，LI H B，XIA X，et al. Blast-induced rock damage control in Fangchenggang nuclear power station，China[J]. Journal of Rock Mechanics and Geotechnical Engineering，2018，10(5)：914–923.
[6] 夏 祥，李俊如，李海波，等. 广东岭澳核电站爆破开挖岩体损伤特征研究[J]. 岩石力学与工程学报，2007，26(12)：2 510–2 516. (XIA Xiang，LI Junru，LI Haibo，et al. Study on damage characteristics of rock mass under blasting load in Ling?Ao nuclear power station，Guangdong province[J]. Chinese Journal of Rock Mechanics and Engineering，2007，26(12)：2 510–2 516.(in Chinese))
[7] 梁生贤，王 桥，焦彦杰，等. LSQR法在位场反演中的分析与评价[J]. 物探与化探，2019，43(2)：359–366.(LIANG Shengxian，WANG Qiao，JIAO Yanjie，et al. Analysis and evaluation of the potential field inversion using LSQR method[J]. Geophysical and Geochemical Exploration，2019，43(2)：359–366.(in Chinese))
[8] GÖKTÜRKLER G，BALKAYA Ç. Traveltime tomography of crosshole radar data without ray tracing[J]. Journal of Applied Geophysics，2010，72(4)：213–224.
[9] 石林珂，孙懿斐. 声波层析成像技术[J]. 岩石力学与工程学报，2003，22(1)：122–125.(SHI Linke，SUN Yifei. Technique of acoustic tomography[J]. Chinese Journal of Rock Mechanics and Engineering，2003，22(1)：122–125.(in Chinese))
[10] HAACH V G，RAMIREZ F C. Qualitative assessment of concrete by ultrasound tomography[J]. Construction and Building Materials，2016，119：61–70.
[11] BECK A，TEBOULLE M. A fast iterative shrinkage-thresholding algorithm for linear inverse problems[J]. SIAM Journal on Imaging Sciences，2009，2(1)：183–202.
[12] 姚旭龙，张艳博，孙 林，等. 基于区域相关性的岩石损伤声发射探测与成像方法研究[J]. 岩石力学与工程学报，2017，36(9)：2 113– 2 123.(YAO Xulong，ZHANG Yanbo，SUN Lin，et al. Research on acoustic emission inspecting and imaging method for rock damage based on regional correlation[J]. Chinese Journal of Rock Mechanics and Engineering，2017，36(9)：2 113–2 123.(in Chinese))
[13] 张艳博，王科学，姚旭龙，等. 基于波速场成像技术的岩石损伤评价研究[J]. 岩石力学与工程学报，2019，38(12)：2 404–2 417. (ZHANG Yanbo，WANG Kexue，YAO Xulong，et al. Rock damage evaluation based on wave velocity field imaging technology[J]. Chinese Journal of Rock Mechanics and Engineering，2019，38(12)：2 404–2 417.(in Chinese))
[14] HE T M，ZHAO Q，HA J，et al. Understanding progressive rock failure and associated seismicity using ultrasonic tomography and numerical simulation[J]. Tunnelling and Underground Space Technology，2018，81：26–34.
[15] 巩思园，田鑫元，郑有雷，等. 弹性波CT反演识别煤岩体空区实验研究[J]. 采矿与安全工程学报，2020，37(4)：759–766.(GONG Siyuan，TIAN Xinyuan，ZHENG Youlei，et al. Empty area recognition technology of coal and rock mass by elastic wave CT inversion[J]. Journal of Mining and Safety Engineering，2020，37(4)：759–766.(in Chinese))
[16] HOSSEINI N，ORAEE K，SHAHRIAR K，et al. Passive seismic velocity tomography on longwall mining panel based on simultaneous iterative reconstructive technique(SIRT)[J]. Journal of Central South University，2012，19：2 297–2 306.
[17] DUAN C L，YAN C H，XU B T，et al. Crosshole seismic CT data field experiments and interpretation for karst caves in deep foundations[J]. Engineering Geology，2017，228：180–196.
[18] SAGONG M，PARK C S，LEE B，et al. Cross-hole seismic technique for assessing in situ rock mass conditions around a tunnel[J]. International Journal of Rock Mechanics and Mining Sciences，2012，53：86–93.
[19] 蒋元建，王李管，彭平安. 复杂不规则地质体的速度模型反演[J]. 煤炭学报，2018，43(增2)：420–425.(JIANG Yuanjian，WANG Liguan，PENG Ping?an. Inversion of velocity model for complex irregular geology[J]. Journal of China Coal Society，2018，43(Supp.2)：420–425.(in Chinese))
[20] 卢回忆，刘伊克，常 旭. 基于MSFM的复杂近地表模型走时计算[J]. 地球物理学报，2013，56(9)：3 100–3 108.(LU Huiyi，LIU Yike，CHANG Xu. MSFM-based travel-times calculation in complex near-surface model[J]. Chinese Journal of Geophysics，2013，56(9)：3 100–3 108.(in Chinese))
[21] DONG L J，TONG X J，MA J. Quantitative investigation of tomographic effects in abnormal regions of complex structures[J]. Engineering，2021，7(7)：1 011–1 022.
[22] 熊 炜，周曾辉，余开彪，等. 基于弯曲路径的混凝土超声层析成像及其改进[J]. 岩土力学，2011，32(2)：629–634.(XIONG Wei，ZHOU Zenghui，YU Kaibiao，et al. Concrete ultrasonic tomography imaging and improvement based on curved path[J]. Rock and Soil Mechanics，2011，32(2)：629–634.(in Chinese))
[23] TRIVINO L F，MOHANTY B. Assessment of crack initiation and propagation in rock from explosion-induced stress waves and gas expansion by cross-hole seismometry and FEM-DEM method[J]. International Journal of Rock Mechanics and Mining Sciences，2015，77：287–299.
[24] DONG L J，PEI Z W，XIE X，et al. Early identification of abnormal regions in rock-mass using traveltime tomography[J]. Engineering，2023，22：191–200.
[25] ZIMA B，K?DRA R. Detection and size estimation of crack in plate based on guided wave propagation[J]. Mechanical Systems and Signal Processing，2020，142：106788.
[26] 邴 琦，孙章庆，韩复兴，等. 地震波射线追踪方法综述——方法，分类，发展现状与趋势[J]. 地球物理学进展，2020，35(2)：536– 547.(BING Qi，SUN Zhangqing，HAN Fuxing，et al. Review on the seismic wave ray tracing：methods，classification，developmental present and trend[J]. Progress in Geophysics，2020，35(2)：536– 547.(in Chinese))
[27] HASSOUNA M S，FARAG A A. Multistencils fast marching methods：A highly accurate solution to the eikonal equation on cartesian domains[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence，2007，29(9)：1 563–1 574.
[28] QU X X，LIU S X，WANG F. A new ray tracing technique for crosshole radar traveltime tomography based on multistencils fast marching method and the steepest descend method[C]// Proceedings of the 15th International Conference on Ground Penetrating Radar. Brussels：[s. n.]，2014：503–508.
[29] BANADAKI M D，MOHANTY B. Numerical simulation of stress wave induced fractures in rock[J]. International Journal of Impact Engineering，2012，40：16–25.
[30] WU D，LI H B，FUKUDA D，et al. Development of a finite-discrete element method with finite-strain elasto-plasticity and cohesive zone models for simulating the dynamic fracture of rocks[J]. Computers and Geotechnics，2023，156：105271.
[31] DUVALL W I. Strain-wave shapes in rock near explosions[J]. Geophysics，1953，18(2)：310–323.
[32] YANG R S，DING C X，YANG L Y，et al. Visualizing the blast-induced stress wave and blasting gas action effects using digital image correlation[J]. International Journal of Rock Mechanics and Mining Sciences，2018，112：47–54.
[33] BANADAKI M M. Stress-wave induced fracture in rock due to explosive action[M]. Toronto：University of Toronto，2010：219.
[34] 魏 东，陈 明，叶志伟，等. 基于应变率相关动力特性的岩体爆破破坏区范围研究[J]. 工程科学与技术，2021，53(1)：67–74. (WEI Dong，CHEN Ming，YE Zhiwei，et al. Study on blasting failure zone based on rate-dependent dynamic characteristics of rock mass[J]. Advanced Engineering Sciences，2021，53(1)：67–74.(in Chinese))
2024, Vol. 43 
