一种隧道围岩微震源定位方法技术

本发明专利技术公开了一种隧道围岩微震源定位方法，该方法基于启发式类算法中的引力搜索法GSA，在隧道掌子面后方布置至少4个微震监测传感器，采集传感器接收范围内的岩体破裂的波形信号；建立隧道空间坐标系并准确测定传感器空间坐标，拾取各个传感器上波形信号的观测到时；并设定观测到时与计算到时的累积绝对差值为微震源定位的目标函数，最终通过计算目标函数满足终止准则来搜索微震源位置，本发明专利技术能够在较快时间内获取较准确的微震源位置，能够满足现场工程对定位精度需要，提高了微震震源定位精度，方法经济实用，操作简单，适用于各种交通、水利水电等隧道工程微震监测。

A method for locating microseismic source of surrounding rock in tunnel

The invention discloses a method for locating micro-seismic source in surrounding rock of tunnel. Based on the gravitational search method GSA of heuristic algorithm, at least four micro-seismic monitoring sensors are arranged behind the face of tunnel face to collect the waveform signals of rock mass rupture within the range of sensor reception, establish the space coordinate system of tunnel and accurately determine the sensor. The space coordinate picks up the observation time of the waveform signal on each sensor and sets the accumulative absolute difference between the observation time and the calculation time as the target function for locating the microseismic source. Finally, the position of the microseismic source can be searched by calculating the target function to satisfy the termination criterion. The invention can obtain a more accurate microseismic source in a relatively fast time. The method is economical, practical and easy to operate. It is suitable for micro-earthquake monitoring in various traffic, water conservancy and Hydropower tunnels.

【技术实现步骤摘要】
一种隧道围岩微震源定位方法
本专利技术涉及微震监测领域，尤其是一种基于启发式算法GSA的隧道围岩微震源定位方法，适用于各种交通、水利水电等隧道工程微震监测。
技术介绍
岩体在外界扰动的影响下，内部会产生微裂隙并以弹性波的形式释放应变能，微裂隙不断发育伴随着弹性波在岩体内迅速传播与释放，这种弹性波被称之为微震。微震监测技术就是基于弹性波的解译来分析岩体内部微裂纹扩展以及岩体稳定性的监测方法。微震源定位是微震监测技术的核心，它是利用微震传感器记录的微震波形信号、到时数据和微震波波速反演微震事件的空间坐标和发震时刻。目前，在微震源定位研究中，震源定位原理主要分为两大类：一类是基于到时不同理论的震源定位方法，另一类则是基于三轴传感器的震源定位方法。基于到时不同理论发展起来的震源定位方法种类繁多，是应用最广的一类震源定位方法，例如经典的Geiger法、Thurber法、单纯形定位算法、双重残差法等震源定位方法。但是，经典的Geiger法对初始条件的依赖性较大，在迭代过程中存在失稳发散问题；Thurber法虽然引入二阶偏导数提高了算法的稳定性，但同时也大大增加了计算量；单纯形定位算法与Powell等直接算法无须求解方程组，而是直接进行多维搜索，具有收敛速度快，简单且易实现等优点，但此类算法过于依赖初值的选取，与最小二乘法一样易陷入局部极小值。目前的现有技术还存在获取难度高、工作量较大、误差较大和效率不高的缺陷。因此，研究一种在较快时间内获取较准确的微震源位置的方法十分有意义，能够很大程度上保证微震源定位精度。
技术实现思路
本专利技术为解决上述技术问题是提供一种能够在较快时间内获取较准确的隧道围岩微震源定位方法。一种隧道围岩微震源定位方法一种隧道围岩微震源定位方法，其特征在于，包括以下步骤：a.在隧道施工掌子面后方布置n个微震检测传感器，n≥4；b.收集微震检测传感器接收范围内的岩体破裂的波形信号，定义所有传感器的观测到时与计算到时的累积绝对差值为震源波速反演的目标函数；c.假定微震定位空间有多个粒子，并假定各个粒子的质量与多维位置向量，；d.以各粒子的多维位置向量计算目标函数，记录得到的最小目标函数值及对应最小目标函数值的位置；e.判断最小函数值是否小于规定量值ε，若最小目标函数值小于规定量值ε，则对应最小目标函数值的位置为震源，若最小目标函数值不小于规定量值ε，则更新各粒子的位置向量，重新计算目标函数值，至最小目标函数值小于规定量值ε为止。进一步，更新假定粒子的位置向量时采用引力搜索算法GSA(GravitationalSearchAlgorithm)进行更新。引力搜索算法GSA(GravitationalSearchAlgorithm)是基于万有引力定律和牛顿第二定律的种群优化算法，万有引力是自然界4种基础力之一。在自然界中，万有引力的作用无处不在，使得任意一个粒子都会与其它的粒子相互吸引而不断的靠近，即较大质量的粒子能吸引较小质量的粒子，较大质量粒子可代表较准确的微震源位置，从而最终获取到最准确的微震源位置。因此，利用基于自然界物理法则的引力搜索算法能实现对微震源的准确定位与获取。进一步，所述的规定量值为ε，其范围为1e-4到1e-5。如上所述的一种隧道围岩微震源定位方法，包括如下步骤：a在隧道施工掌子面后方布置至少4个微震监测传感器。b建立隧道空间坐标系，采集传感器接收范围内的岩体破裂的波形信号，定义所有传感器的观测到时与计算到时的累积绝对差值为震源波速反演的目标函数，所述目标函数的计算公式如下：式中，fit为到时的累积绝对差值，n为监测传感器数量；ti为第i个传感器的观测到时，上标p,s为P波或S波，t0为震源的初始发震时刻，Ri/V为计算走时，Ri为微震源位置(x1，x2，x3)与第i个传感器位置(xi，yi，zi)的距离，V代表微震波在传播路径上的速度。c假定微震定位空间有N个粒子(代表计算震源)。在初始时刻，每个粒子有质量Mi和多维位置向量Xi：M＝(M1,...,Mi,...,MN)，(i＝1,2,...,N)(2)式2代表有N个具有质量的粒子，式3代表每个粒子具有n维的数值(若仅为三维空间则退化为三维数值)，其中，表示第i个粒子在第d维的数值，且有上下限值，即d将各粒子的多维位置向量X(包含震源的位置和微震波速)带入式1，获得各粒子的目标函数值。记录历史循环目标函数的最小值fitbest及对应粒子的多维位置Xbest。e判断当前目标函数最小值是否满足终止条件(残差值是否小于规定量值ε)，即是否fitbest<ε。如果fitbest<ε则输出最优微震源位置，如果fitbest＞ε执行下一步骤，所述ε的范围为1e-4到1e-5，ε取值越小，获取的震源定位精度就越高，但需要运算时间也会越久，因此ε的具体取值与工程实际情况相关。f计算粒子间相互作用的引力。在第k次迭代，定义为在d维度上粒子i受到粒子j作用的引力：其中，Maj(k)和Mpi(k)分别为主动粒子j的惯性质量和被动粒子i的惯性质量，ε为小量值常量。G(k)为引力系数函数，满足如下：其中，G0与α为已知值；k为当前迭代次数，K为迭代总次数。Rij(k)为粒子i和粒子j的欧式距离：Rij(k)＝||Xi(k),Xj(k)||2(6)第d维上第i个粒子受到其他所有粒子引力作用的总和为：式中，randj为[0，1]之间的随机数，为引力的总和添加随机。在每一次迭代中，每个粒子都会更新惯性质量。惯性质量根据目标函数值计算，粒子惯性质量越大，表明越接近最优值，也表明对其他粒子有更大的吸引力。根据以下公式更新粒子的惯性质量Mi：Mai＝Mpi＝Mii＝Mi,(i＝1,2,...,N)(8)其中，fiti(k)为粒子i在第k次迭代的目标函数值大小。对于求解最小目标函数值问题，best为N个粒子目标函数的最小值，worst为N个粒子目标函数的最大值，定义如下：g计算每个粒子的加速度和速度。根据牛顿第二定理，第d维上粒子i的加速度和速度为：其中，Mi(k)为当前时刻粒子i的惯性质量。h在每一次迭代中，每个粒子都会根据以下公式更新位置：i根据各粒子更新后的位置向量重新计算目标函数值，并判断是否fitbest<ε，若满足则退出循环，输出全体粒子目标函数的最小值fitbest和对应的最优震源位置向量Xbest；否则继续执行迭代循环。作为优选，G0为100，α为20。采用本专利技术所述的方法对隧道围岩微震源进行定位，具备以下有益效果：(1)提供了一种新的方法用于微震源定位，且该方法与传统的微震源定位方法相比，经济实用、操作简单、参数设置较少、收敛速度快，且能够很好的和各种优化问题相结合。(2)本专利技术所采用的方法的通用性比较强，对问题中不确定的信息具有一定的适应能力。(3)本专利技术所采用的方法能够在较快时间内获取较准确的微震源位置，能够满足现场工程对定位精度需要，提高了微震震源定位精度。附图说明图1是万有引力作用示意图；图2是本专利技术提供的隧道围岩微震源定位方法的流程图。具体实施方式下面结合附图及具体实施内容对本专利技术进一步说明。以某高地应力高速公路双洞隧道安装微震监测系统为例，依照图2流程获取微震源位置。在该高地应力高速公路双洞隧道安本文档来自技高网...

【技术保护点】
1.一种隧道围岩微震源定位方法，其特征在于，包括以下步骤：a.在隧道施工掌子面后方布置n个微震检测传感器，n≥4；b.收集微震检测传感器接收范围内的岩体破裂的波形信号，定义所有传感器的观测到时与计算到时的累积绝对差值为震源波速反演的目标函数；c.假定微震定位空间有多个粒子，并假定各个粒子的质量与多维位置向量；d.以各粒子的多维位置向量计算目标函数，记录得到的最小目标函数值及对应最小目标函数值的位置；e.判断最小函数值是否小于规定量值ε，若最小目标函数值小于规定量值ε，则对应最小目标函数值的位置为震源，若最小目标函数值不小于规定量值ε，则更新各粒子的位置向量，重新计算目标函数值，至最小目标函数值小于规定量值ε为止。

【技术特征摘要】
1.一种隧道围岩微震源定位方法，其特征在于，包括以下步骤：a.在隧道施工掌子面后方布置n个微震检测传感器，n≥4；b.收集微震检测传感器接收范围内的岩体破裂的波形信号，定义所有传感器的观测到时与计算到时的累积绝对差值为震源波速反演的目标函数；c.假定微震定位空间有多个粒子，并假定各个粒子的质量与多维位置向量；d.以各粒子的多维位置向量计算目标函数，记录得到的最小目标函数值及对应最小目标函数值的位置；e.判断最小函数值是否小于规定量值ε，若最小目标函数值小于规定量值ε，则对应最小目标函数值的位置为震源，若最小目标函数值不小于规定量值ε，则更新各粒子的位置向量，重新计算目标函数值，至最小目标函数值小于规定量值ε为止。2.根据权利要求1所述的一种隧道围岩微震源定位方法，其特征在于：所述步骤e中更新各粒子的位置向量采用引力搜索算法进行更新。3.根据权利要求1所述的一种隧道围岩微震源定位方法，其特征在于：所述的规定量值为ε，其范围为1e-4到1e-5。4.根据权利要求1或2或3所述的一种隧道围岩微震源定位方法，其特征在于，包括以下步骤：a.在隧道施工掌子面后方布置n个微震监测传感器，n≥4；b.建立隧道空间坐标系，采集传感器接收范围内的岩体破裂的波形信号，定义所有传感器的观测到时与计算到时的累积绝对差值为震源波速反演的目标函数，所述目标函数的计算公式如下：式1中，fit为到时的累积绝对差值，n为监测传感器数量；ti为第i个传感器的观测到时，上标p,s为P波或S波，t0为震源的初始发震时刻，Ri/V为计算走时，Ri为微震源位置与第i个传感器位置的距离，V代表微震波在传播路径上的速度；c.假定微震定位空间有N个粒子，在初始时刻，每个粒子有质量Mi和多维位置向量Xi：M＝(M1,...,Mi,...,MN)，(i＝1,2,...,N)(2)式2代表有N个具有质量的粒子，式3代表每个粒子具有n维的数值(若仅为三维空间则退化为三维数值)，其中，表示...

【专利技术属性】
技术研发人员：马春驰，李天斌，张航，韩瑀萱，周雄华，王剑锋，
申请(专利权)人：成都理工大学，
类型：发明
国别省市：四川,51