基于变步长加速搜索的微震源定位方法

微震定位方法是微震监测技术的重要组成部分,其关键是定位震源位置。利用空间网格划分并计算网格交点目标函数值,对微震定位目标函数二维及三维空间分布进行了分析,并据此获取了目标函数连续且极小值唯一、单轴收敛范围逐步减小、各轴收敛范围不一的规律。利用以上规律及模式搜索法、网格搜索法的优缺点,探索出了基于连续比较模块、变步长模块、加速模块的变步长加速搜索法。通过模拟算例与工程数据下收敛稳定性、结果精确度、计算速度以及参数初始值影响程度4个指标的效果对比,结果表明：模拟算例下,对比模拟退火算法、遗传算法,变步长加速搜索法的目标函数值标准差、定位误差标准差、波速误差标准差均为0;该算法的定位误差平均值分别为其余二者的0.7%、1.9%;该算法的计算时间平均值分别为其余二者的6.9%、33.2%。该算法单独更改各参数对定位误差的影响在0.005～0.025 m之间;减小搜索步长下限可有效提高结果精确度,并增加相应的计算时间。在规定初至到时与目标函数模型及检波器位置坐标下,搜索算法对定位精度无实质影响。     

基于快速行进迎风线性插值的微震定位算法研究

准确的微震定位方法是微震监测技术进行岩体稳定性分析的基础,而初至波走时计算方法是决定微震定位精度重要因素之一。引入一种由快速行进法(fast marching methods, FMM)与线性插值射线追踪法(LTI ray tracing)改进而来的线性插值算法,实现了在复杂岩体中相对准确的微震初至走时计算,并在均一速度模型和Marmousi模型中得到了验证。在此基础上,采用常微分方程数值算法中的单步高精度Runge-Kutta法,得到了震源点和传感器之间的射线路径,并在分层和带空洞的岩体模型中验证了其合理性。选取最小二乘法为目标函数,通过寻找目标函数值最小的方式来获得最优节点实现微震定位。在所建立的加入地下硐室的Marmousi岩体模型中,对于预设震源点的定位结果误差分别为1.41、1.00、2.83、1.41、1.41m,验证了所提方法在复杂岩体速度模型中能够实现较准确的微震定位,为进一步提高工程岩体的微震定位精度提供了新思路。     

微震震源定位新方法——混合差定位法

当微地震事件个数较多时,传统双差定位算法的数据存储量和计算量巨大,难以适应微震震源的定位。提出一种新的微震震源定位算法——混合差定位法。这种算法综合了双差定位算法和Geiger定位法的优点,不仅数据的存储量和计算量大大减小,且与双差定位相比,其定位结果不受初始震源丛的质心位置的影响。模型数据和实际微地震数据测试结果验证了该方法的有效性和实用性。     

基于微震矩张量的矿山围岩破坏机制分析

微震是矿山岩体损伤的重要信号,并被用于矿山围岩稳定性分析。为应用矩张量理论,研究矿山开采围岩破坏机制,对理论运用中的关键问题提出了解决办法:(1)通过引入适于矿山岩体应力状态的矩张量三分量分解模型,解决了由于矿山应力状态不确定带来的矩张量分解方程组求解的不确定问题;(2)针对Ohstu提出的矩张量岩体破坏模式判据的不足之处进行了讨论和修正,使其适于矿山岩体工程岩体破坏模式分析。利用改进后的矩张量分析方法对冬瓜山铜矿的部分巷道围岩破坏机制进行研究,印证了冬瓜山巷道围岩破坏是一种以剪切破坏为主,多因素综合作用的破坏形式的结论。当巷道穿过不同岩层接触带时,微震事件在接触带附近积聚,表明巷道围岩出现剪切破坏带并形成冒落区;位于单一岩层中的巷道围岩中的微震活动较为分散,破坏原因多为围岩应力集中,破坏形式表现为局部性围岩冒落。     

井中微震监测数据处理方法研究与应用

井中微震监测的目的是确定压裂过程产生的裂缝的几何形态,采集到具有微震事件的原始记录信噪比低,纵横波初至拾取困难,直接影响到微震的定位及后续的裂缝形态描述。通过分析但由于微震信息较弱,因此微震事件的记录特点,结合实际应用过程中遇到的问题对处理方法进行了研究,包括去噪处理、微震事件识别、纵横波初至拾取、极化分析、模型建立及校正、微震定位等。研究表明,带通滤波可有效法除噪音;长、短时窗能量比法可精确的记录微震事件及初至时间自动拾取;利用纵波尖端曲线图可以确定微震源方向;速度模型的建立应以校正炮时距曲线拟合理论曲线为原则;利用Geiger和网格搜索法可精确得到震源定位结果。     

基于差分进化算法的微震定位

针对传统基于走时微震定位方法存在的定位精度依赖于走时精度问题,提出了改进定位目标函数的差分进化微震定位方法。基于传统微震P波走时定位方法,在定位目标函数中加入了误差加权系数,通过降低初至拾取精度低的微震信号道对定位误差的影响,减弱了初至拾取精度对微震定位的影响;利用差分进化算法,通过种群初始化、变异、交叉和选择等操作,求解改进的微震定位目标函数,实现微震定位;对模型数据和实测数据进行试处理,结果表明：改进的目标函数能明显降低误差较大分量对总误差的贡献;提出的改进定位目标函数的差分进化微震定位方法,在部分微震信号道走时拾取误差较大甚至错误的情况下,能很大程度减弱拾取误差对定位精度的影响,可以获得较好的定位结果。     

深埋隧道岩爆位置偏转及其微震活动特征

川藏铁路拉林段控制性工程巴玉隧道,施工过程中岩爆频发且呈现邻近区域岩爆位置不一现象。为此,通过微震监测系统原位观测了岩爆位置偏转过程中的岩体破裂演化规律,进一步探讨了岩爆位置发生偏转的原因。研究结果表明:岩爆位置偏转段两侧微震信息呈不均匀分布特征,微震活动(事件数、能量和视体积等参数)随岩爆偏转向岩爆发生位置集核偏转,两者具有良好的空间相关性。岩爆位置偏转段岩体损伤区及视应力的演化特征说明岩爆位置偏转与局部应力集中区变化相关,应力集中程度大的区域更易发生岩爆。相关研究成果对岩爆位置偏转预警及防治具有一定借鉴和参考意义。     

基于EMS微震参数的岩爆预警方法及探讨

一系列深埋长大隧道工程在我国西部涌现,深埋隧道围岩的岩爆成因机制及监测预警成为亟待解决的课题。微震监测作为岩体微破裂空间监测技术在隧道工程岩爆预警上发挥了重要作用。由于微震数据存在的波动性（时间序列）、离散性（空间分布）和误差,利用震源参数随施工循环的演化特征评估岩爆发育过程并非足够稳定、有效。基于3种微震监测易于获取且紧密联系的震源参数（地震能量、地震矩和视应力）,扩展建立了用于岩爆评估与预警的EMS方法。该方法定义了地震能量、地震矩、视应力3种评价指标,并获得了3种评价指标间的关系;将3种评价指标用于岩爆预测,并进行了岩爆等级评估;同时还给出了微震路径与震源参数空间的概念。最后,将EMS方法应用于巴陕高速米仓山隧道的岩爆灾害实例,对岩爆里程段微震事件簇（微震事件聚集区域）内岩体累进破裂和岩爆发育过程进行了评估,实现了简易有效的微震数据日常处理与岩爆灾害评估预警。     

缓倾地层微震定位算法探讨及其数值验证

倾斜层状地层中速度模型对于震源精确定位具有重要意义。基于水平层状地层两点间的快速射线追踪法,分析了应力波在倾斜地层中的传播路径及规律。结合网格搜索法,提出了针对倾斜地层的微震定位算法,并与传统的基于单一速度模型的定位算法进行对比分析;探索性地采用了颗粒流理论构建地层数值模型,通过颗粒间的相互作用模拟应力波传播进行定位算法有效性验证。研究结果表明：（1）对于倾斜地层,基于简化弹性波等速传播规律的定位算法精度较差,单一速度模型不能满足复杂介质的定位要求;（2）采用倾斜分层速度模型并改进应力波传播路径的计算方法能有效提高倾斜地层的震源定位精度,同时分析了地层倾角以及地层数对算法精度的影响,表明定位误差与倾角、地层数成正比;（3）通过减小网格搜索法的网格尺寸可在一定程度上提高定位精度,合理的尺寸大小和搜索策略有助于网格搜索法在实际工程中的应用。建立的倾斜地层震源定位及其验证方法可为进一步研究复杂地层定位监测技术提供重要的理论和技术支撑。     

基于稳健模拟退火-单纯形混合算法的微震定位研究

为提高矿井微震的定位精度,通过模拟试验对比分析常规定位及新方法在定位求解中的适用范围及优缺点,研究反演未知变量数目、检波器的密度、波速等对各自定位方法的影响。基于单纯形法的优点,并结合模拟退火法的全局收敛性,提出采用稳健的模拟退火-单纯形法进行微震定位。研究发现,该方法不必直接搜索发震时刻,只需通过稳健方法预先估算,利用较少的检波器就可进行微震定位,且基本不受初始波速影响,从而大大减少搜索步数,有效地提高定位精度及收敛速度。以柿竹园矿人工爆破定位测试试验为例,重点分析比较模拟退火-单纯形法、稳健模拟退火-单纯形法及其他定位方法的空间绝对距离及目标函数误差,认为检波器数量较多且可获得理想的P波信息时,应用非线性最小二乘法能求出可信的结果;其余定位方法的三维定位误差及目标函数误差相差不大,但稳健模拟退火-单纯形法目标函数误差最小,稳定性最高,可在矿井微震定位中推广使用。     

基于STA/LTA岩石破裂微震信号实时识别算法及工程应用

微震监测获取的数据中通常混有大量的非岩石破裂信号,该类信号目前主要通过人工经验进行识别与滤除,这消耗了大量的宝贵时间,严重影响灾害的防治和救援效率。对大量微震信号进行分析,发现STA/LTA算法在信号实时触发后能大致表征波形振幅和频率的变化,岩石破裂信号和非岩石破裂信号在延迟位置处R值具有差异性。基于此,提出了岩石破裂微震信号实时识别算法。新算法应用到白鹤滩水电站地下厂房、红透山和阿舍勒铜矿深部采场3个工程,岩石破裂事件识别的准确率分别是85.98%、92.45%和91.06%,非岩石破裂事件滤除的准确率分别是72.06%、83.11%和49.87%。该算法使基于岩石破裂微震信息的岩石工程灾害自动分析与预警成为可能,具有重要意义。     

太宁隧道围岩类别的再认识

通过论证比较，按照隧道围岩主要工程地质特征，围岩岩体结构及结构面特征，水文地质特征，围岩弹性纵波波速，变形系数，泊松比，内摩擦角，软化系数等指标，重新划分和确定了围岩类别，分析了产生围岩类别偏差的原因，以期对同类工程勘察起到借鉴和参考作用。     

深埋隧洞微震活动区与岩爆的相关性研究

基于锦屏二级水电站深埋引水隧洞和排水洞大量微震监测数据及上百个不同等级岩爆实例,研究了深埋隧洞微震活动区与岩爆之间的关系。研究结果表明：（1）微震活动分布范围主要介于掌子面后方3倍洞径至前方1.5倍洞径之间,而岩爆高发区位于掌子面后方3倍洞径以内,表明岩爆高发区与微震事件主要分布范围相吻合;（2）隧洞工程岩爆潜在风险重点关注区域是掌子面后方3倍洞径已开挖范围,以及掌子面前方1.5倍洞径施工范围;（3）微震事件及岩爆分布呈区域性集结特点,其中一部分岩爆发生于微震事件集结区内部,另一部分岩爆发生于微震事件集结区边缘,这是岩体破坏过程中所固有的现象,与微震事件集结区边缘局部应力集中密切相关。     

浅谈坑道围岩稳定性监测

结合延吉市八道金矿坑探工程,利用FLAC3D软件建立坑道三维实体模型,对坑道围岩的应力分布规律、变形位移大小及塑性区分布等方面进行围岩稳定性的综合分析。然后根据分析结果对施工坑道围岩形变大的部位进行监测,为坑道施工安全及支护方案提供了一定的参考依据。     

基于? -SVR算法的隧道围岩位移演化规律预测

针对目前广泛应用的灰色理论、遗传算法(GA)和人工神经网络(ANN)等方法预测隧道稳定性的缺陷，提出应用稳健性能较好的? -SVR（support vector machine）算法对非对称连拱隧道围岩位移演化规律进行预测研究。应用加速混合遗传算法搜索? -SVR最优参数，以提高? -SVR的预测能力。将预测结果与灰色理论、BP神经网络预测结果进行比较，显示? -SVR算法学习和预测精度高。     

基于ANP-KL-TOPSIS法的铁路隧道围岩安全性评价

为客观评价隧道围岩安全状态,结合围岩岩体结构及地质特征等影响因素,建立隧道围岩安全评价指标体系。用ANP法应用Super Decision软件建立指标因素关联图并计算出各指标权重。在传统TOPSIS法的基础上结合相对熵原理,用Kull-back-Leibler距离计算贴近度,构建隧道围岩安全评价模型。运用该模型对张吉怀铁路兰新乡隧道进口段七处隧道段的围岩进行安全性评价,得出各隧道段围岩安全等级,评价结果切合实际,对隧道施工有一定的指导意义。该评价模型简单可操作,可运用于隧道围岩安全性评价研究中。     

强震区隧道软弱围岩洞口段桩-筏抗震措施的作用效果分析

为进一步提高地震时软弱围岩铁路隧道安全性及稳定性,以玉磨铁路新平隧道工程作为依托,利用FLAC3D数值模拟软件,对强震区隧道软弱围岩洞口段桩-筏抗震措施效果进行研究,对比分析了无加固措施、措施一（桩-拱结构）及措施二(桩-筏-拱结构)二衬的结构位移及内力,研究结果表明：措施一和措施二的竖向位移分别减小了1.35%、1.09%;措施一的边墙收敛几无变化（增大了0.44%）;由于桩-筏-拱结构的整体稳定性较好,措施二的边墙收敛减小了30.49%;距离洞口最近的两个危险截面1和2处,措施一的最小安全系数提高了1.43%、6.71%,措施二的最小安全系数提高了145.91%、143.72%;综合位移及内力分析,措施二的抗震效果优于措施一,建议新平隧道洞口段采用桩-筏-拱结构进行抗震加固。