CN108802814A - 一种隧道围岩微震波速的获取方法 - Google Patents
一种隧道围岩微震波速的获取方法 Download PDFInfo
- Publication number
- CN108802814A CN108802814A CN201810634022.1A CN201810634022A CN108802814A CN 108802814 A CN108802814 A CN 108802814A CN 201810634022 A CN201810634022 A CN 201810634022A CN 108802814 A CN108802814 A CN 108802814A
- Authority
- CN
- China
- Prior art keywords
- wave
- velocity
- particle
- microseism
- tunnel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000012544 monitoring process Methods 0.000 claims abstract description 22
- 239000002360 explosive Substances 0.000 claims abstract description 13
- 238000009825 accumulation Methods 0.000 claims abstract description 6
- 230000000694 effects Effects 0.000 claims abstract description 5
- 239000002245 particle Substances 0.000 claims description 72
- 239000011435 rock Substances 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims description 7
- 238000010845 search algorithm Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 4
- 230000003993 interaction Effects 0.000 claims description 2
- 210000002415 kinetochore Anatomy 0.000 claims description 2
- 238000004880 explosion Methods 0.000 abstract description 6
- 238000004422 calculation algorithm Methods 0.000 abstract description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract description 2
- 238000005422 blasting Methods 0.000 abstract 1
- 238000012360 testing method Methods 0.000 description 7
- 238000013459 approach Methods 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 238000009533 lab test Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000001960 triggered effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000005293 physical law Methods 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/288—Event detection in seismic signals, e.g. microseismics
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种隧道围岩微震波速的获取方法,该方法基于启发式类算法中的引力搜索法GSA,在隧道掌子面后方布置至少4个微震监测传感器,建立隧道空间坐标系并准确测定传感器空间坐标,利用人工***事件触发传感器震动响应,拾取各个传感器上波形信号的观测到时;并设定观测到时与计算到时的累积绝对差值为微震波速反演的目标函数,最终通过计算目标函数满足终止准则来搜索微震波波速,本发明能够利用隧道施工的***事件实时准确获取隧道围岩的微震波速,保证了微震定位精度;该方法不需要专门的***活动确定微震波速,经济实用,操作简单,适用于各种交通、水利水电等隧道工程微震监测。
Description
技术领域
本发明涉及微震监测领域,尤其是一种基于启发式算法GSA的隧道围岩微震波速获取方法,适用于各种交通、水利水电等隧道工程微震监测。
背景技术
微震是指岩体应力超过自身强度时,岩体内部产生裂纹错动、开裂和破坏,累积的能量以弹性波形式释放并产生震动。微震震源定位是微震监测技术的核心,它是利用微震传感器记录的微震波形信号、到时数据和微震波波速反演微震事件的空间坐标和发震时刻。准确的微震波波速对于提高震源定位精度和稳定性至关重要,同时也是较难精确确定的震源定位输入参数。
目前,微震震源定位过程中,主要还是通过室内试验和现场勘探来获取微震波波速。室内试验是指在室内通过对岩块试样进行声波试验而获取岩块的波速;现场勘探主要指在微震监测区域内利用现场试验反演微震波波速。由于微震监测区域地质构造和介质特性的复杂多变性,现场勘探只能反映微震监测区域中某些区域的微震波波速,并不能完全代表监测区域内所有路径下的微震波波速;而室内声波试验中所用的岩块并不能代表现场岩体的性质,将所测试的波速作为微震波速用于微震监测必然造成较大的定位误差。
现有技术中,中国专利CN103697999A公开了一种高应力硬岩TBM施工隧道微震波速实时获取方法,在高应力硬岩TBM施工隧道掌子面后方岩体内布置至少4个微震传感器,采集掌子面附近岩体破裂产生的震动信号,将掌子面附近发生的岩体破坏事件作为微震源,测量微震源及微震传感器坐标,提取微震传感器监测到的微震源震动信号到时,反演获取微震波速。该方法反演的微震波速基于天然微震事件,由于岩体破裂位置未知,可能造成反演的微震波速存在较大误差。
现有技术中,中国专利CN104406681B公开的一种实时确定微震波波速的测试方法,利用微震监测***自带的微震传感器和已知震源的人工***实时确定微震波波速。在微震监测***已有的空间坐标系内准确测定各微震传感器和人工***事件的空间坐标;拾取微震监测***中各微震传感器监测到的微震波到时,将第一个被触发的微震传感器到时作为比较对象,分别采用其它微震传感器和第一个被触发微震传感器的到时差以及它们相对于人工***震源的距离差实时反演微震波波速。随着采掘工程的进行。但是,该方法仅能获取单一的微震波速,不能获取微震波速的向量或不同传播路径上的波速,可能出现应用该微震波速进行定位时的较大误差。
现有技术中,中国专利CN104502964A公开了一种基于空间几何关系的获得微震波速的方法,根据震源点和基站的空间几何关系,构建微震等效波速V的不等式,构建形如开口向上的二次函数的变化规律的时差适应函数,提出一个单调单向搜索算法,求出最佳波速。但是,该方法仅能获取单一的微震波速,不能获取微震波速的向量或不同传播路径上的波速,可能出现应用该微震波速进行定位时的较大误差。
目前的现有技术还存在获取难度高、工作量较大、误差较大和效率不高的缺陷。因此,研究一种能够综合利用现场已知震源的人工***事件和已有的微震监测***实时测试微震波波速的方法十分有意义,能够很大程度上保证微震震源定位精度。
发明内容
本发明为解决上述技术问题是提供一种能够综合利用现场已知震源的人工***事件和已有的微震监测***实时测试微震波波速的方法。
本发明所采用的技术方案是,一种隧道围岩微震波速的获取方法,其特征在于,包括如下步骤:
a.在隧道施工掌子面后方布置n个微震检测传感器,n≥4;
b.采用人工***事件为震源位置,收集微震检测传感器接收范围内的岩体破裂的波形信号,定义所有传感器的观测到时与计算到时的累积绝对差值为震源波速反演的目标函数;
c.假定微震定位空间有多个粒子,并假定各个粒子的质量与多维波速向量;
d.以各粒子的多维波速向量计算目标函数,记录得到的最小目标函数值及对应最小目标函数值的多维波速向量;
e.判断最小函数值是否小于规定量值ε,若最小目标函数值小于规定量值ε,则对应最小目标函数值的多维波速向量为波速,若最小目标函数值不小于规定量值ε,则更新各粒子的多维波速向量,重新计算目标函数值,至最小目标函数值小于规定量值ε为止。
进一步,所述步骤e中各粒子的多维波速向量的更新采用引力搜索算法GSA(Gravitational Search Algorithm)进行。
引力搜索算法GSA(Gravitational Search Algorithm)是基于万有引力定律和牛顿第二定律的种群优化算法,万有引力是自然界4种基础力之一。在自然界中,万有引力的作用无处不在,使得任意一个粒子都会与其它的粒子相互吸引而不断的靠近,即较大质量的粒子能吸引较小质量的粒子,较大质量粒子可代表较准确的微震波速模型,从而最终获取到最准确的微震波速。因此,利用基于自然界物理法则的引力搜索算法能实现对微震源的准确定位与获取。
进一步,所述的规定量值ε,其范围为1e-4到1e-5。
作为优选,如上所述的一种隧道围岩微震波速的获取方法,包括如下步骤:
a在隧道施工掌子面后方布置至少4个微震监测传感器;
b采用人工***事件,以人工***事件位置为震源,建立隧道空间坐标系,定义所有传感器的观测到时与计算到时的累积绝对差值为震源波速反演的目标函数,所述目标函数的计算公式如下:
式中,fit为到时的累积绝对差值,n为传感器数量;ti为第i个传感器的观测到时,上标p,s为P波或S波,t0为震源的初始发震时刻,Ri/V为计算走时,Ri代表人工***事件位置(x0,y0,z0)与第i个传感器位置(xi,yi,zi)的距离;vi代表震源与第i个传感器传播路径上的微震波速,对于所有的传感器则有微震波速向量V;
c假定隧道空间有N个粒子,每个粒子代表一种可能的微震波速模型,在初始时刻,每个粒子有质量Mi和n维波速向量Vi;
M=(M1,...,Mi,...,MN),(i=1,2,...,N) (2)
式2代表有N个具有质量的粒子,式3代表每个粒子具有n维的波速。其中,表示第i个粒子在第d维的波速,并有上下限值,即
d将式(2)带入式(1)中计算各粒子的目标函数值fit,并记录历史循环的最小目标函数值Fbest及其对应粒子波速向量Vbest;
e判断当前最小目标函数值Fbest是否满足终止准则,所述终止准则为最小目标函数值Fbest是否小于规定量值ε,即是否Fbest<ε。如果Fbest<ε则输出最优微震波波速Vbest;如果Fbest>ε则执行下一步骤,所述ε的范围为1e-4到1e-5,ε取值越小,获取的微震波波速精度就越高,但需要运算时间也会越久,因此ε的具体取值与工程实际情况相关;
f计算粒子间相互作用的引力,在第k次迭代,定义为在d维度上粒子i受到粒子j作用的引力,计算公式如下:
式中,Maj(k)为主动粒子j的惯性质量,Mpi(k)为被动粒子i的惯性质量,为粒子j在第d维度上的波速,为粒子i在第d维度上的波速,ε为小量值常量,Rij(k)为粒子i和粒子j的欧式距离,G(k)为引力系数函数,G(k)和Rij(k)满足如下公式:
式中,G0与α为确定值,k为当前迭代次数,K为迭代总次数;随着当前迭代次数逐渐增加到总迭代次数,引力系数逐渐衰减;
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个粒子目标函数的最大值,worst(t)和best(t)的定义如下:
g计算每个粒子的加速度。根据牛顿第二定理,第d维上粒子i的的加速度为:
式中,Fi d为当前迭代次数粒子i受到的引力总和,Mi(k)为粒子i的惯性质量;
h在每一次迭代中,每个粒子都会根据以下公式更新粒子速度和微震波速
i将步骤h中各粒子更新的微震波速带入式(1)计算目标函数值fit,并判断是否Fbest<ε。若满足则退出循环,输出全体粒子目标函数的最小值Fbest和对应的最优微震波波速Vbest;若不满足,则继续执行迭代循环。
进一步,G0可取100,α可取20。
本发明的有益效果是:
(1)采用施工活动***事件作为微震震源反演微震波波速,人工***事件的空间坐标能够采用全站仪精确测量,从而避免了因震源空间坐标误差导致的微震波波速误差,提高了微震波波速的精度。
(2)本方法不需要进行专门***获取微震波速,经济实用,操作简单,参数设置较少,收敛速度快,且能够很好的和各种优化问题相结合,且通用性比较强,对问题中不确定的信息具有一定的适应能力。
(3)能够实时获取微震波波速,对微震震源定位中采用的微震波波速进行实时调整,从而解决了由于隧道施工导致的监测区域中微震波波速可能发生变化的问题,保证了微震波波速的实时准确性,提高了微震震源定位精度。
附图说明
图1是万有引力作用示意图。
图2是本发明隧道围岩微震波速获取方法流程图。
具体实施方式
下面结合附图与实施例对本发明做进一步说明。
以某高地应力高速公路双洞隧道安装微震监测***为例,依照图2流程获取微震波速。在该高地应力高速公路双洞隧道安装微震监测***中,该监测***含有9个微震传感器,分别布置在左右双洞掌子面后方,其中滞后洞布置3个,超前洞布置6个,分2个断面布置。由此开展针对隧道围岩内部破裂发生的实时监测,捕捉破裂触发的微震波形信号,拾取微震波形信号在每个传感器上的观测到时。在所建立的隧道空间坐标系内采用全站仪精确测得各传感器空间坐标(见表1)。2017年3月9日15时37分30秒在微震监测范围内进行了人工***试验,采用全站仪准确测量该人工***事件空间坐标为M(1008,770,1010)。同时,提取对各微震传感器所接受的微震波形到时(见表1)。
表1各微震传感器坐标
定义所有传感器的观测到时与计算到时的累积绝对差值为震源波速反演的目标函数,其计算公式如下:
在上述基础上,基于启发式算法GSA获取隧道围岩微震波速,其参数选取如下:K=1000,G0=100,α=20。最终根据上述发明,计算的最小目标函数Fbest=5.44E-4,对应的粒子波速向量Vbest=(5601,5566,4178),即该实例反演获取的最优微震波波速v1=5601(对应传感器S1、S2和S3),v2=5566(对应传感器S4、S5和S6),v3=4178(对应传感器S7、S8和S9)。
以上实例仅用以说明本发明的技术方案而非限制,本领域的普通技术人员应当理解,本发明的技术方案进行修改或同等替换,而不脱离本发明方案的精神和范围,均应覆盖在本发明中。
Claims (5)
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个传感器位置的距离;vi代表震源与第i个传感器传播路径上的微震波速,对于所有的传感器则有微震波速向量V;
c假定隧道空间有N个粒子,每个粒子代表一种可能的微震波速模型,在初始时刻,每个粒子有质量Mi和n维波速向量Vi;
M=(M1,...,Mi,...,MN),(i=1,2,...,N) (2)
式2代表有N个具有质量的粒子,式3代表每个粒子具有n维的波速,其中,表示第i个粒子在第d维的波速,并有上下限值,即
d将式(3)带入式(1)中计算各粒子的目标函数值fit,并记录历史循环的最小目标函数值Fbest及其对应粒子波速向量Vbest;
e判断当前最小目标函数值Fbest是否满足终止准则,所述终止准则为最小目标函数值Fbest是否小于规定量值ε,如果Fbest<ε则输出最优微震波波速Vbest;如果Fbest>ε则执行下一步骤;
f计算粒子间相互作用的引力,在第k次迭代,定义为在d维度上粒子i受到粒子j作用的引力,计算公式如下:
式中,Maj(k)为主动粒子j的惯性质量,Mpi(k)为被动粒子i的惯性质量,为粒子j在第d维度上的波速,为粒子i在第d维度上的波速,ε为小量值常量,,Rij(k)为粒子i和粒子j的欧式距离,G(k)为引力系数函数,G(k)和Rij(k)满足如下公式:
式中,G0与α为确定值,k为当前迭代次数,K为迭代总次数;;
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个粒子目标函数的最大值,worst(t)和best(t)的定义如下:
g计算每个粒子的加速度,根据牛顿第二定理,第d维上粒子i的加速度为:
式中,Fi d为当前迭代次数粒子i受到的引力总和,Mi(k)为粒子i的惯性质量;
h在每一次迭代中,每个粒子都会根据以下公式更新粒子速度和微震波速
i将步骤h中各粒子更新的微震波速带入式(1)计算目标函数值fit,并判断是否Fbest<ε,若满足则退出循环,输出全体粒子目标函数的最小值Fbest和对应的最优微震波波速Vbest;若不满足,则继续执行迭代循环。
5.根据权利要求4所述的一种隧道围岩微震波波速的获取方法,其特征在于:G0为100,α为20。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810634022.1A CN108802814B (zh) | 2018-06-20 | 2018-06-20 | 一种隧道围岩微震波速的获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810634022.1A CN108802814B (zh) | 2018-06-20 | 2018-06-20 | 一种隧道围岩微震波速的获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108802814A true CN108802814A (zh) | 2018-11-13 |
CN108802814B CN108802814B (zh) | 2019-10-25 |
Family
ID=64083839
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810634022.1A Expired - Fee Related CN108802814B (zh) | 2018-06-20 | 2018-06-20 | 一种隧道围岩微震波速的获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108802814B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109782356A (zh) * | 2019-02-25 | 2019-05-21 | 西南大学 | 基于能量网格搜索的井下微震监测传感器最优布设方法 |
CN110333530A (zh) * | 2019-06-25 | 2019-10-15 | 广东石油化工学院 | 一种新的微震事件检测方法和*** |
CN113885074A (zh) * | 2021-09-09 | 2022-01-04 | 中煤科工开采研究院有限公司 | 微震监测能量校核方法 |
CN117572503A (zh) * | 2024-01-15 | 2024-02-20 | 湖北神龙工程测试技术有限公司 | 一种岩石波速测试方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120120762A1 (en) * | 2009-05-25 | 2012-05-17 | Masahiro Kamata | Methods and systems for seismic signal detection |
CN105022031A (zh) * | 2015-07-03 | 2015-11-04 | 四川大学 | 一种区域岩体微震震源的分层速度定位方法 |
CN105652315A (zh) * | 2016-01-19 | 2016-06-08 | 山东大学 | 地下工程施工过程中多通道声发射监测***及定位方法 |
CN106154334A (zh) * | 2015-04-13 | 2016-11-23 | 中石化石油工程地球物理有限公司胜利分公司 | 基于网格搜索的井下微地震事件实时反演定位方法 |
CN107562537A (zh) * | 2017-08-21 | 2018-01-09 | 中南大学 | 一种基于万有引力搜索的云计算任务调度方法 |
CN107884822A (zh) * | 2017-11-13 | 2018-04-06 | 北京矿冶研究总院 | 一种提高矿用微震震源定位精度的方法 |
-
2018
- 2018-06-20 CN CN201810634022.1A patent/CN108802814B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120120762A1 (en) * | 2009-05-25 | 2012-05-17 | Masahiro Kamata | Methods and systems for seismic signal detection |
CN106154334A (zh) * | 2015-04-13 | 2016-11-23 | 中石化石油工程地球物理有限公司胜利分公司 | 基于网格搜索的井下微地震事件实时反演定位方法 |
CN105022031A (zh) * | 2015-07-03 | 2015-11-04 | 四川大学 | 一种区域岩体微震震源的分层速度定位方法 |
CN105652315A (zh) * | 2016-01-19 | 2016-06-08 | 山东大学 | 地下工程施工过程中多通道声发射监测***及定位方法 |
CN107562537A (zh) * | 2017-08-21 | 2018-01-09 | 中南大学 | 一种基于万有引力搜索的云计算任务调度方法 |
CN107884822A (zh) * | 2017-11-13 | 2018-04-06 | 北京矿冶研究总院 | 一种提高矿用微震震源定位精度的方法 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109782356A (zh) * | 2019-02-25 | 2019-05-21 | 西南大学 | 基于能量网格搜索的井下微震监测传感器最优布设方法 |
CN110333530A (zh) * | 2019-06-25 | 2019-10-15 | 广东石油化工学院 | 一种新的微震事件检测方法和*** |
CN113885074A (zh) * | 2021-09-09 | 2022-01-04 | 中煤科工开采研究院有限公司 | 微震监测能量校核方法 |
CN113885074B (zh) * | 2021-09-09 | 2024-04-02 | 中煤科工开采研究院有限公司 | 微震监测能量校核方法 |
CN117572503A (zh) * | 2024-01-15 | 2024-02-20 | 湖北神龙工程测试技术有限公司 | 一种岩石波速测试方法 |
CN117572503B (zh) * | 2024-01-15 | 2024-03-26 | 湖北神龙工程测试技术有限公司 | 一种岩石波速测试方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108802814B (zh) | 2019-10-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108802814B (zh) | 一种隧道围岩微震波速的获取方法 | |
CN108717201B (zh) | 一种隧道围岩微震源定位方法 | |
Golubev et al. | Numerical simulation of seismic activity by the grid-characteristic method | |
CN105022031A (zh) | 一种区域岩体微震震源的分层速度定位方法 | |
CN106840381A (zh) | 一种***振动衰减效应的测试方法 | |
CN106680869B (zh) | 微地震事件的检测和定位方法与装置 | |
CN109738940B (zh) | 一种存在空区条件下的声发射/微震事件定位方法 | |
CN109597125B (zh) | 一种基于p波到时与最大振幅波形的微震源定位方法 | |
CN109991658B (zh) | 一种基于“震源-台站”速度模型的微地震事件定位方法 | |
CN105242328B (zh) | 古热岩石圈厚度的确定方法及装置 | |
CN113189644B (zh) | 一种微震震源定位方法及*** | |
CN112305591B (zh) | 隧道超前地质预报方法、计算机可读存储介质 | |
CN107290787A (zh) | 一种地震次声同址台站的监测信号关联方法 | |
CN104749630B (zh) | 构建微地震监测速度模型的方法 | |
D'Alessandro et al. | Retrieval of ocean bottom and downhole seismic sensors orientation using integrated MEMS gyroscope and direct rotation measurements | |
CN110398775A (zh) | 隧道突涌水灾害微震事件信号波动初至拾取方法及*** | |
CN103364823A (zh) | 震动源实时定位与分析*** | |
CN109061723B (zh) | 一种隧洞岩爆孕育过程的微震源高精度定位方法及*** | |
CN105022091B (zh) | 一种无预测速的远场震源快速定位方法 | |
CN109709602A (zh) | 一种远探测声波偏移成像方法、装置及*** | |
CN103941281A (zh) | 一种分布式矿震检测方法及装置 | |
CN105301638A (zh) | 一种提取风化层底界面的方法和装置 | |
CN109521221B (zh) | 一种钻爆法施工硬岩隧道微震波波速实时获取方法 | |
US11079468B2 (en) | Detection of acoustic events | |
CN112525407B (zh) | 一种冲击波单体多传感器阵列测试方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191025 Termination date: 20210620 |