(1. State Key Laboratory of Hydraulics and Mountain River Engineering，Sichuan University，Chengdu，Sichuan 610065，China；
2. China Three Gorges Corporation，Beijing 100038，China)
摘要快速精准的微震震源定位方法是微震监测技术更好地发挥岩体稳定性预测预警作用的基础。精确的初至波走时计算是提高微震定位精度的关键。针对实际工程中带有空洞和速度分区的复杂波速岩体，引入速度快、精度高和对复杂模型适应性强的多模板快速行进法multi-stencils fast marching methods(MSFM)用于初至波走时计算。该方法利用坐标旋转的方式生成新的计算模板，使网格对角邻点参与计算，提高了the Fast Marching Method(FMM)在对角方向的计算精度。首先，对比分析了一、二阶FMM和MSFM计算精度和效率，表明二阶MSFM具有更高的计算精度；其次，在分层速度岩体和带空洞岩体中分别采用单一速度模型和二阶MSFM计算初至波走时，通过与解析解对比发现，二阶MSFM较单一速度模型的计算绝对误差平均减小了97.65%和95.18%；然后，建立速度分层且带有隧洞的岩体模型，验证了二阶MSFM算法对复杂速度模型适应性极强的特点；最后，提出了基于MSFM的复杂速度岩体微震定位算法，并将其应用到白鹤滩水电站左岸边坡微震定位。通过对现场采集的4个爆破事件进行定位试算，得到使定位误差平均值最小(9.6 m)的相对最优层速度组合并用于考虑空洞的速度模型构建，爆破事件走时正演验证了本速度模型的可靠性。对2015年5月的128个微震事件分别采用单一速度模型和提出的定位算法进行定位，通过对比定位效果，发现后者定位的事件较前者在岩体主要损伤区(层内错动带LS331和LS337)附近有更明显呈条带状分布的聚集现象。研究表明MSFM算法作为走时正演方法具有很好的应用价值，提出的基于MSFM的微震定位方法能够有效地提高震源定位精度。
Abstract：A rapid and accurate microseismic(MS) source location method is the foundation of MS monitoring to play a better forecasting and early-warning role for rock mass stability. The accurate travel time calculation of first arrival wave is the key to improve the location accuracy of MS source. Aiming at rock mass with complex velocities in actual engineering，like voids filled with air and partitioned velocities，the multi-stencils fast marching methods(MSFM) was introduced to calculate the travel time of first arrival wave，which is fast and with high precision and strong applicability for complex model. In this method，a new calculation stencil is generated by means of coordinate rotation to make the diagonal adjacent points of the mesh calculated，and the precision of the fast marching method(FMM) is improved in diagonal direction. Firstly，the accuracy and efficiency of the first or second order FMM and MSFM were compared and analyzed，showing that the second order MSFM has higher accuracy. Secondly，for rock mass with partitioned velocities model and one with a void，both single velocity model and the second order MSFM were used to calculate the travel time. Through comparing them to analytical solution，the average absolute errors of the second order MSFM were reduced by 97.65% and 95.18% compared with the single velocity model. Then，the rock mass model with partitioned velocities model and tunnels was built，and it was proved that the second order MSFM has strong applicability for complex model. Finally，the MSFM-based MS source location algorithm of rock mass with complex velocities was proposed，and applied to Baihetan hydropower station left bank slope of MS location. By trial location of four blasting events collected in fields，the relative optimal layer velocities combination was obtained，which had the minimum mean value(9.6 m) of location errors and was used to establish the void-considered velocity model. The reliability of this velocity model was validated by the travel time forward calculation of a blasting event. The 128 MS events in May 2015 were located by both the single velocity model and the proposed location algorithm，by comparing the location results，the located events by the latter were distributed as belts more obviously near the main damage zone(LS331 and LS337) than those by the former. The study concludes that the MSFM has very valuable application as a travel time calculation method，the proposed MSFM-based MS source location method can improve the source location precision effectively.
[1] XU N W，TANG C A，LI L C，et al. Microseismic monitoring and stability analysis of the left bank slope in Jinping first stage hydropower station in southwestern China[J]. International Journal of Rock Mechanics and Mining Sciences，2011，48(6)：950–963.
[2] 姜福兴，杨淑华，成云海，等. 煤矿冲击地压的微地震监测研究[J]. 地球物理学报，2006，49(5)：1 511–1 516.(JIANG Fuxing，YANG Shuhua，CHENG Yunhai，et al. A study on microseismic monitoring of rock burst in coal mine[J]. Chinese Journal of Geophysics，2006，49(5)：1 511–1 516.(in Chinese))
[3] 陆菜平，窦林名，王耀峰，等. 坚硬顶板诱发煤体冲击破坏的微震效应[J]. 地球物理学报，2010，53(2)：450–456.(LU Caiping，DOU Linming，WANG Yaofeng，et al. Microseismic effect of coal material rock burst failure induced by hard roof[J]. Chinese Journal of Geophysics，2010，53(2)：450–456.(in Chinese))
[4] 潘一山，赵扬锋，官福海，等. 矿震监测定位系统的研究及应用[J]. 岩石力学与工程学报，2007，26(5)：1 002–1 011.(PAN Yishan，ZHAO Yangfeng，GUAN Fuhai，et al. Study on rockburst monitoring and orientation system and its application[J]. Chinese Journal of Rock Mechanics and Engineering，2007，26(5)：1 002–1 011.(in Chinese))
[5] 唐礼忠，杨承祥，潘长良. 大规模深井开采微震监测系统站网布置优化[J]. 岩石力学与工程学报，2006，25(10)：2 036–2 042.(TANG Lizhong，YANG Chengxiang，PAN Changliang. Optimization of microseismic monitoring network for large-scale deep well mining[J]. Chinese Journal of Rock Mechanics and Engineering，2006，25(10)：2 036–2 042.(in Chinese))
[6] 李庶林，尹贤刚，郑文达，等. 凡口铅锌矿多通道微震监测系统及其应用研究[J]. 岩石力学与工程学报，2005，24(12)：2 048– 2 053.(LI Shulin，YIN Xiangang，ZHENG Wenda，et al. Research of multi-channel microseismic monitoring system and its application to Fankou lead-zinc mine[J]. Chinese Journal of Rock Mechanics and Engineering，2005，24(12)：2 048–2 053.(in Chinese))
[7] 陈际经，胡静云，李庶林. 柿竹园矿地下特大爆破的微震监测技术研究[J]. 矿业研究与开发，2011，31(2)：59–62.(CHEN Jijing，HU Jingyun，LI Shulin，et al. Microseismic monitoring for the underground extra large-scale blasting in Shizhuyuan Mine[J]. Mining Research and Development，2011，31(2)：59–62.(in Chinese))
[8] 宋维琪，徐奔奔，喻志超，等. 基于各向异性分析的微地震震源矢量场重建和裂缝解释[J]. 地球物理学报，2015，58(2)：656–663. (SONG Weiqi，XU Benben，YU Zhichao，et al. Reconstruction of the micro-seismic source vector field and fissure interpretation based on the anisotropy analysis[J]. Chinese Journal of Geophysics，2015，58(2)：656–663.(in Chinese))
[9] XU N，TANG C，LI H，et al. Excavation-induced microseismicity：microseismic monitoring and numerical simulation[J]. Journal of Zhejiang University Science A，2012，13(6)：445–460.
[10] 徐奴文，唐春安，周 钟，等. 岩石边坡潜在失稳区域微震识别方法[J]. 岩石力学与工程学报，2011，30(5)：893–900.(XU Nuwen，TANG Chunan，ZHOU Zhong，et al．Identification method of potential failure regions of rock slope using microseismic monitoring technique[J]. Chinese Journal of Rock Mechanics and Engineering，2011，30(5)：893–900.(in Chinese))
[11] 李 彪，戴 峰，徐奴文，等. 深埋地下厂房微震监测系统及其工程应用[J]. 岩石力学工程学报，2014，33(增1)：3 375–3 383. (LI Biao，DAI Feng，XU Nuwen，et al. Microseismic monitoring system and its engineering applications of deep-buried underground powerhouse[J]. Chinese Journal of Rock Mechanics and Engineering，2014，33(Supp.1)：3 375–3 383.(in Chinese))
[12] 张伯虎，邓建辉，周志辉，等. 大岗山水电站厂房断层控制区域微震监测分析[J]. 岩土力学，2012，33(增2)：213–218.(ZHANG Bohu，DENG Jianhui，ZHOU Zhihui，et al. Analysis of monitoring microseism in areas controlled by faults near powerhouse in Dagangshan hydropower station[J]. Rock and Soil Mechanics，2012，33(Supp.2)：213–218.(in Chinese))
[13] TRIFU C I，SHUMILA V. Geometrical and inhomogeneous raypath effects on the characterization of open-pit seismicity[C]// 44th US Rock Mechanics Symposium and 5th US-Canada Rock Mechanics Symposium. [S.l.]：American Rock Mechanics Association，2010.
[14] 巩思园，窦林名，马小平，等. 煤矿矿震定位中异向波速模型的构建与求解[J]．地球物理学报，2012，55(5)：1 757–1 763.(GONG Siyuan，DOU Linming，MA Xiaoping，et al. Study on the construction and solution technique of anisotropic velocity model in the location of coal mine tremor[J]. Chinese Journal of Geophysics，2012，55(5)：1 757–1 763.(in Chinese))
[15] FENG G L，FENG X T，CHEN B R，et al. Sectional velocity model for microseismic source location in tunnels[J]. Tunnelling and Underground Space Technology，2015，45：73–83.
[16] UM J，THURBER C. A fast algorithm for two-point seismic ray tracing[J]. Bulletin of the Seismological Society of America，1987，77(3)：972–986.
[17] LANGAN R T，LERCHE I，CUTLER R T. Tracing of rays through heterogeneous media：An accurate and efficient procedure[J]. Geophysics，1985，50(9)：1 456–1 465.
[18] 马争鸣，李衍达. 二步法射线追踪[J]. 地球物理学报，1991，34(4)：501–508.(MA Zhengming，LI Yanda. Two-step ray tracing method[J]. Chinese Journal of Geophysics，1991，34(4)：501–508.(in Chinese))
[19] 高尔根，徐果明. 二维速度随机分布逐步迭代射线追踪方法[J]. 地球物理学报，1996，39(增1)：302–308.(GAO Ergen，XU Guoming. A new kind of step by step iterative ray-tracing method[J]. Chinese Journal of Geophysics，1996，39(Supp.1)：302–308.(in Chinese))
[20] 李 强. 基于程函方程快速行进法地震走时成像方法研究[硕士学位论文][D]. 西安：长安大学，2012.(LI Qiang. Seismic traveltime tomography by using a Fast Marching Method based on solving Eikonal equation[M. S. Thesis][D]. Xi?an：Chang?an University，2012.(in Chinese))
[21] VIDALE J. Finite-difference calculation of travel times[J]. Bulletin of the Seismological Society of America，1988，78(6)：2 062–2 076.
[22] VIDALE J E. Finite-difference calculation of traveltimes in three dimensions[J]. Geophysics，1990，55(5)：521–526.
[23] NELSON G D，VIDALE J E. Earthquake locations by 3D finite-difference travel times[J]. Bulletin of the Seismological Society of America，1990，80(2)：395–410.
[24] SETHIAN J A. A fast marching level set method for monotonically advancing fronts[J]. Proceedings of the National Academy of Sciences，1996，93(4)：1 591–1 595.
[25] VAN TRIER J，SYMES W W. Upwind finite-difference calculation of traveltimes[J]. Geophysics，1991，56(6)：812–821.
[26] PODVIN P，LECOMTE I. Finite difference computation of traveltimes in very contrasted velocity models：a massively parallel approach and its associated tools[J]. Geophysical Journal International，1991，105(1)：271–284.
[27] KIM S. ENO-DNO-PS：a stable，second-order accuracy eikonal solver[C]//1999 SEG Annual Meeting. [S.l.]：Society of Exploration Geophysicists，1999：1 747–1 750.
[28] CHOPP D L. Some improvements of the fast marching method[J]. SIAM Journal on Scientific Computing，2001，23(1)：230–244.
[29] ZHAO H. A fast sweeping method for Eikonal equations[J]. Mathematics of Computation，2005，74：603–627.
[30] 李庆春，李永博，叶 佩. 分区多步快速行进法射线路径计算方法[C]// 中国地球物理学会第二十八届年会. [S.l.]：[s.n.]，2012：570.(LI Qingchun，LI Yongbo，YE Pei. The Raypath calculating method of multistage fast marching method[C]// The 28th annual meeting of Chinese Geophysical Society. [S.l.]：[s.n.]，2012：570.(in Chinese))
[31] HASSOUNA S M，FARAG A A. Multistencils fast marching methods：A highly accurate solution to the Eikonal equation on cartesian domains[J]. Pattern Analysis and Machine Intelligence，IEEE Transactions on，2007，29(9)：1 563–1 574.
[32] 刘 斌. 地震学原理与应用[M]. 合肥：中国科学技术大学出版社，2009：100–104.(LIU Bin. Principles and applications of seismology[M]. Hefei：Press of University of Science and Technology of China，2009：100–104.(in Chinese))
[33] 卢回忆，刘伊克，常 旭. 基于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))
[34] DIJKSTRA E W. A note on two problems in connexion with graphs[J]. Numerische mathematik，1959，1(1)：269–271.
[35] 董陇军，李夕兵，唐礼忠，等. 无需预先测速的微震震源定位的数学形式及震源参数确定[J]. 岩石力学与工程学报，2011，30(10)：2 057–2 067.(DONG Longjun，LI Xibing，TANG Lizhong，et al．Mathematical functions and parameters for microseismic source location without pre-measuring speed[J]. Chinese Journal of Rock Mechanics and Engineering，2011，30(10)：2 057–2 067.(in Chinese))
[36] 华东勘测设计研究院有限公司. 金沙江白鹤滩水电站左岸坝肩开挖(834～600 m)工程施工招标文件参考资料——工程地质[R]. [S.l.]：[s.n.]，2013.(HydroChina Huadong Engineering Corporation Limited. Bidding documents reference of the left bank abutment excavation engineer construction(834–600 m) in Baihetan hydropower station[R]. [S.l.]：[s.n.]，2013.(in Chinese))
[37] 华东勘测设计研究院有限公司. 金沙江白鹤滩水电站左岸拱轴线工程地质剖面图[R]. [S.l.]：[s.n.]，2015.(HydroChina Huadong Engineering Corporation Limited. The engineering geological profile of arch axis in Baihetan hydropower station[R]. [S.l.]：[s.n.]，2015.(in Chinese))
2017, Vol. 36 
