CN109521221B - 一种钻爆法施工硬岩隧道微震波波速实时获取方法 - Google Patents
一种钻爆法施工硬岩隧道微震波波速实时获取方法 Download PDFInfo
- Publication number
- CN109521221B CN109521221B CN201811501964.9A CN201811501964A CN109521221B CN 109521221 B CN109521221 B CN 109521221B CN 201811501964 A CN201811501964 A CN 201811501964A CN 109521221 B CN109521221 B CN 109521221B
- Authority
- CN
- China
- Prior art keywords
- wave
- microseismic
- wave velocity
- sensor
- velocity
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P3/00—Measuring linear or angular speed; Measuring differences of linear or angular speeds
- G01P3/64—Devices characterised by the determination of the time taken to traverse a fixed distance
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提出一种钻爆法施工硬岩隧道微震波波速实时获取方法,具体流程为:已知坐标的点进行单孔单药包***;提取每个传感器监测到的该信号的P波监测到时;计算第i+1号传感器和第i号传感器的P波监测到时之差与计算到时之差的残差;求解微震P波波速的适应值函数Q;求得微震P波波速;求得S波波速的估计值;进行初步定位,并提取P波和S波监测到时,将P波和S波监测到时绘制在直角坐标系中,并进行直线拟合,求得实际P波和S波的波速比,进而求得微震S波波速;本发明将微震P波波速和S波波速分开求解,提高S波波速反演精度的同时,不会增加人力物力的消耗,方法经济简洁,操作可行,实时快速获取岩体中微震波波速,提高微震监测效果。
Description
技术领域
本发明属于隧道微震监测技术领域,具体涉及一种钻爆法施工硬岩隧道微震波波速实时获取方法。
背景技术
微震监测是指通过监测岩体破裂产生的震动或其它物体震动,对监测对象的破坏状况、安全状况等作出评价,从而为预警和防控灾害提供依据的成套设备和技术。微震震源定位是微震监测的基础,它是利用微震***记录的微震波形信息、岩体中微震波传播速度,包括P波波速和S波波速,求解微震事件的空间坐标、发震时刻并计算相关微震参数,如微震能量、震级、视体积、震源机制等。因此,微震波波速直接影响微震源的定位精度和微震参数的计算准确度,进而影响微震信息的解译和监测效果。
目前,获取微震波波速的方法主要有室内试验法和现场测试法。室内试验法主要是指在室内通过对岩块进行声波试验而获取岩块波速,然后乘以一定的系数得到岩体的波速。由于现场地质条件复杂,室内试验得到的岩块波速与现场岩体波速之间并不是简单的线性关系。因此,通过室内试验法获得的岩体波速并不能准确反映现场岩体波速。现场测试法主要是指现场通过在已知空间坐标的点进行单孔单药包***,然后利用微震监测***中各传感器监测到该***产生的微震波的P波和S波到时经过反演获得岩体的波速。然而,通常情况下,相对于正常的微震事件,***产生的微震波S波较弱,很难准确识别S波监测到时。因此,采用现场定点***反演获得的微震S波波速误差较大,如果将其作为岩体波速输入,进行微震源定位和相关微震参数计算,必然带来较大的误差。
中国专利公开号CN102096093A,发明名称“一种利用微震点作为震源计算矿区地震波传播速度的方法”,该发明针对矿区地震波传播速度的计算,利用多个传感器测得的微震点产生的地震波时间确定微震点的坐标估计,然后利用微震点的坐标估计值来计算矿区地震波传播速度。由于微震点坐标估计值精度不高,因此该法所获取的地震波传播速度存在较大误差。中国专利公开号CN103697999A,发明名称“一种高应力硬岩TBM施工隧道微震波速实时获取方法”,该发明针对高应力硬岩TBM施工隧道,利用掌子面附近的塌方、岩爆或岩体开裂作为坐标已知的微震源,然后进行波速反演以获得微震波波速。该法只有在施工过程中发生了宏观岩体破坏的洞段可以应用,而没有发生宏观岩体破坏的洞段难以实施。中国专利公开号CN104502964A,发明名称“一种基于空间几何关系的获得微震波速的方法”,该发明提供了一种基于传感器和未知震源的空间几何关系,利用迭代算法构建等效微震波速的方法,但该法只适用于微震事件分布在传感器阵列内的情况,而对于隧道工程,微震事件主要分布在传感器阵列外,该法难以有效应用。中国专利公开号CN104406681A,发明名称“一种实时确定微震波速的测试方法”,该发明利用已知震源的人工***通过反演实时确定微震波波速,但该方法只对P波波速给出了反演方法,没有给出S波波速的反演方法。
发明内容
针对上述问题,本发明的目的在于提供一种钻爆法施工硬岩隧道微震波波速实时获取方法,用于解决硬岩钻爆法施工隧道微震监测中波速难以实时准确获取的问题,提高微震源定位精度和微震参数计算准确度,从而保证微震信息被正确解译并提高微震监测效果。
一种钻爆法施工硬岩隧道微震波波速实时获取方法,包括如下步骤:
步骤1:在硬岩钻爆法施工隧道掌子面附近已知坐标的点进行单孔单药包***,掌子面附近定义为:掌子面前方10m至掌子面后方10m,将该已知坐标点作为微震源,利用掌子面后方布置的微震监测***捕捉***产生的微震信号,用于微震P波波速反演;
式中,(x,y,z)为微震源的空间坐标,(xi,yi,zi)为第i号传感器的空间坐标,VP为微震P波波速,其中当i=n时有:
步骤4:求解微震P波波速的适应值函数Q描述为:
步骤5:当Q等于或趋于零时,求得微震P波波速VP;
步骤6:求得微震P波波速VP后,根据岩体中P波和S波的波速比估计值为求得S波波速的估计值VS′。然后利用求得的微震P波波速VP和微震S波波速的估计值VS′对隧道开挖过程中产生的微震事件进行初步定位,选取初步定位坐标在掌子面前方10m至掌子面后方10m范围内,且S波监测到时明显的微震事件进行分析,提取每个传感器监测到的该微震事件的P波和S波监测到时,根据传感器编号分别记录为:
步骤7:第i号传感器的P波和S波监测到时分别满足如下关系:
式中,t0为微震事件的发震时间,Ri为微震源与第i个传感器之间的距离,VS为微震S波波速;
步骤8:公式(4)消除Ri后得到第i号传感器的S波和P波监测到时之差为:
步骤9:以Ti P为横坐标,(Ti S-Ti P)为纵坐标,将所有传感器记录的P波和S波监测到时以散点的形式绘制在直角坐标系中,然后对这些点进行直线拟合,由公式(5)可知为该直线的斜率,即实际的微震P波和微震S波的波速比,微震P波波速步骤4已经求出,则由直线斜率求得微震S波波速。
在隧道施工过程中,随着掌子面的推进,围岩地质情况不断发生变化,微震波波速也将不断发生变化,因此,根据现场施工情况,重复步骤1~步骤9实时动态地获取随掌子面的推进而不断变化的微震波波速。
有益技术效果:
(1)将微震P波波速和S波波速分开求解,避免了采用***事件同时求解P波波速和S波波速时,S波监测到时难以准确拾取而导致的S波波速的求解误差较大的问题。
(2)充分利用钻爆法施工隧道掌子面附近围岩中频繁发生的微破裂事件,选用S波监测到时明显的事件进行S波波速反演,在提高S波波速反演精度的同时,不会增加人力物力的消耗。
(3)直接利用掌子面后方已有的微震监测***进行波速反演,方法经济简洁,操作可行,随着掌子面的推进,能够实时快速获取岩体中微震波波速,解决了隧道开挖过程中岩体微震波波速频繁变化的问题,提高了微震监测效果。
附图说明
图1为本发明实施例的一种钻爆法施工硬岩隧道微震波波速实时获取方法流程图;
图2为本发明实施例的S波监测到时明显的波形图;
图3为本发明实施例的(Ti S-Ti P)随Ti P变化的散点分布及直线拟合。
具体实施方式
下面结合附图和具体实施实例对发明做进一步说明:一种钻爆法施工硬岩隧道微震波波速实时获取方法,如图1所示,包括如下流程:
步骤1:在硬岩钻爆法施工隧道微震监测过程中,微震传感器阵列一般布置于掌子面后方一定距离的洞壁围岩中,且微震传感器阵列一般跟随掌子面向前移动。根据隧道中的测量基准点,可以用全站仪、测斜仪和钢卷尺测得各传感器的空间坐标。掌子面附近任意点的空间坐标也能快速准确地获得。在掌子面附近选择一个点作为微震源,在该点进行单孔单药包***,利用微震监测***捕捉***产生的微震信号,该信号可用于微震P波波速反演;
2017年6月27日,该隧道掌子面开挖至桩号SK194+835处。掌子面后方布置有8个微震传感器,传感器阵列的中心与掌子面间的距离约为100m,各微震传感器的相对坐标按照传感器编号依次记录为(1243.8,6.5,-11.9),(1242.5,5.3,-6.5),(1243.8,-6.1,-7.5),(1242.9,-6.9,-12.0),(1213.5,5.9,-11.2),(1212.6,5.3,-6.9),(1212.6,-6.2,-7.4)和(1212.2,-6.8,-11.2)。
为进行微震波波速反演,施工人员在掌子面附近选取相对空间坐标为(1354.5,0.0,-10.7)的点于2017年6月27日18时00分进行了单孔单药包***;
步骤2:***产生的微震信号震级高,一般能触发监测***中所有的传感器,且***信号P波监测到时非常明显,可以准确地提取每个传感器的P波监测到时,根据传感器编号将各传感器记录的P波监测到时分别记录为:n为传感器个数,本实施例n=8;
微震监测***于2017年6月27日18时00分02秒监测到该***信号。经过分析提取每个传感器捕捉到的该***信号的P波监测到时,按照传感器编号依次记录为(18:00:02.595583),(18:00:02.596749),(18:00:02.595583),(18:00:02.595583),(18:00:02.600833),(18:00:02.600833),(18:00:02.600666),和(18:00:02.600833)。
式中,(x,y,z)为微震源的空间坐标,(xi,yi,zi)为第i号传感器的空间坐标,VP为微震P波波速,其中当i=n时有:
步骤4:求解微震P波波速的适应值函数Q可以描述为:
步骤5:当Q等于或趋于零时,求得的VP即为微震P波波速;
将震源点坐标、传感器坐标和各传感器记录的***信号的P波监测到时代入公式(1)~(3),反演获得该洞段微震P波波速为VP=5995m/s。
步骤6:根据地震方面的研究可知,理论上岩体中P波和S波的波速比估计值为因此在获得微震P波波速后可求得S波波速的估计值为在硬岩钻爆法施工过程中掌子面附近会产生较多的微震事件,利用求得的微震P波波速和S波波速的估计值对这些微震事件进行初步定位。然后选取初步定位坐标在掌子面附近且S波监测到时明显的微震事件进行分析,提取每个传感器监测到的该微震信号的P波和S波监测到时,根据传感器编号将各传感器记录的P波监测到时分别记录为:将各传感器记录的S波监测到时分别记录为:
步骤7:第i号传感器的P波和S波监测到时分别满足如下关系:
式中,t0为微震事件的发震时间,Ri为微震源与第i个传感器之间的距离,VS为微震S波波速;
步骤8:公式(4)前后两式相减可消除Ri,从而得到第i号传感器的S波和P波监测到时之差为:
根据反演获得的P波波速和S波波速的估计值对6月27日微震***捕捉到的微震事件进行初步定位,其中2017年6月27日18时04分04秒监测到的微震事件初步定位在掌子面后方4m,且其S波监测到时明显,如图2所示,经过分析提取每个传感器捕捉到的该微震信号的P波和S波监测到时,按照传感器编号依次记录为(18:04:04.399458,18:04:04.411958),(18:04:04.399958,18:04:04.413291),(18:04:04.399458,18:04:04.412125),(18:04:04.399458,18:04:04.412291),(18:04:04.404645,18:04:04.420979),(18:04:04.404812,18:04:04.422145),(18:04:04.404478,18:04:04.420978)和(18:04:04.404646,18:04:04.421146)。
步骤9:以Ti P为横坐标,(Ti S-Ti P)为纵坐标,将所有传感器记录的P波和S波监测到时以点的形式绘制在直角坐标系中,然后对这些点进行直线拟合,该直线的斜率记为k,由公式(5)可知由此可求得实际的P波和S波的波速比为再结合步骤步骤1~步骤4求得的微震P波波速VP可求得微震S波波速
以各传感器捕捉到的该信号的P波监测到时为横坐标,以对应的S波监测到时与P波监测到时差为纵坐标,在直角坐标系中绘制(Ti S-Ti P)随Ti P变化的散点分布图,然后进行直线拟合,结果如图3所示。该拟合直线的斜率为0.7612,由式(5)可知,所以微震S波波速为VS=5995÷1.7612=3404m/s。
在隧道施工过程中,随着掌子面的推进,新揭露的围岩地质情况不断发生变化,而且微震传感器阵列也在跟随掌子面向前移动,因此微震波波速也将不断发生变化。所以实际监测过程中需要根据现场揭露的围岩地质情况及微震传感器阵列的布置情况,重复步骤1~步骤9实时动态地获取随掌子面的推进而不断变化的微震波波速。
Claims (1)
1.一种钻爆法施工硬岩隧道微震波波速实时获取方法,其特征在于,包括如下步骤:
步骤1:在硬岩钻爆法施工隧道掌子面附近已知坐标的点进行单孔单药包***,掌子面附近定义为:掌子面前方10m至掌子面后方10m,将该已知坐标点作为微震源,利用掌子面后方布置的微震监测***捕捉***产生的微震信号,用于微震P波波速反演;
式中,(x,y,z)为微震源的空间坐标,(xi,yi,zi)为第i号传感器的空间坐标,VP为微震P波波速,其中当i=n时有:
步骤4:求解微震P波波速的适应值函数Q描述为:
步骤5:当Q等于或趋于零时,求得微震P波波速VP;
步骤6:求得微震P波波速VP后,根据岩体中P波和S波的波速比估计值为求得S波波速的估计值VS′,然后利用求得的微震P波波速VP和S波波速的估计值VS′对隧道开挖过程中产生的微震事件进行初步定位,选取初步定位坐标在掌子面前方10m至掌子面后方10m范围内,且S波监测到时明显的岩石破裂微震信号进行分析,提取每个传感器监测到的岩石破裂微震信号的P波监测到时和S波监测到时,根据传感器编号将各传感器记录的P波监测到时分别记为:将各传感器记录的S波监测到时分别记为:
步骤7:第i号传感器的P波和S波监测到时分别满足如下关系:
式中,t0为微震事件的发震时间,Ri为微震源与第i个传感器之间的距离,VS为微震S波波速;
步骤8:公式(4)消除Ri后得到第i号传感器的S波和P波监测到时之差为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811501964.9A CN109521221B (zh) | 2018-12-10 | 2018-12-10 | 一种钻爆法施工硬岩隧道微震波波速实时获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811501964.9A CN109521221B (zh) | 2018-12-10 | 2018-12-10 | 一种钻爆法施工硬岩隧道微震波波速实时获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109521221A CN109521221A (zh) | 2019-03-26 |
CN109521221B true CN109521221B (zh) | 2020-09-29 |
Family
ID=65795026
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811501964.9A Active CN109521221B (zh) | 2018-12-10 | 2018-12-10 | 一种钻爆法施工硬岩隧道微震波波速实时获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109521221B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110107308B (zh) * | 2019-04-30 | 2020-08-25 | 中国科学院武汉岩土力学研究所 | 一种发生在隧洞tbm刀盘内岩爆的等级确定方法 |
CN112130204B (zh) * | 2020-09-02 | 2022-03-04 | 武汉大学 | 一种适用于深埋隧洞围岩波速的现场测试方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5109362A (en) * | 1990-10-22 | 1992-04-28 | Shell Oil Company | Remote seismic sensing |
JPH06324160A (ja) * | 1993-01-12 | 1994-11-25 | Kajima Corp | 震度予測システム |
CN101025084A (zh) * | 2006-02-20 | 2007-08-29 | 中国石油大学(北京) | 一种随钻预测钻头底下地层孔隙压力的方法 |
CN101251603A (zh) * | 2008-04-12 | 2008-08-27 | 中国石油集团西北地质研究所 | X、z两分量波场矢量合成方法 |
US7457195B2 (en) * | 2003-02-08 | 2008-11-25 | Schlumberger Technology Corporation | Estimating the time of arrival of a seismic wave |
CN101329405A (zh) * | 2007-06-20 | 2008-12-24 | 中国石油天然气集团公司 | 一种简单的多参数地震反演方法 |
CN102435980A (zh) * | 2011-09-15 | 2012-05-02 | 中南大学 | 一种基于解析求解的声发射源或微震源定位方法 |
CN105510959A (zh) * | 2015-11-30 | 2016-04-20 | 山东康威通信技术股份有限公司 | 一种隧道的振源类型识别方法及振源位置定位方法 |
CN105785436A (zh) * | 2016-03-17 | 2016-07-20 | 北京矿冶研究总院 | 一种矿用微震监测方法 |
CN106680871A (zh) * | 2016-12-05 | 2017-05-17 | 中南大学 | 一种基于p波到时与波形振幅的震源定位方法 |
-
2018
- 2018-12-10 CN CN201811501964.9A patent/CN109521221B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5109362A (en) * | 1990-10-22 | 1992-04-28 | Shell Oil Company | Remote seismic sensing |
JPH06324160A (ja) * | 1993-01-12 | 1994-11-25 | Kajima Corp | 震度予測システム |
US7457195B2 (en) * | 2003-02-08 | 2008-11-25 | Schlumberger Technology Corporation | Estimating the time of arrival of a seismic wave |
CN101025084A (zh) * | 2006-02-20 | 2007-08-29 | 中国石油大学(北京) | 一种随钻预测钻头底下地层孔隙压力的方法 |
CN101329405A (zh) * | 2007-06-20 | 2008-12-24 | 中国石油天然气集团公司 | 一种简单的多参数地震反演方法 |
CN101251603A (zh) * | 2008-04-12 | 2008-08-27 | 中国石油集团西北地质研究所 | X、z两分量波场矢量合成方法 |
CN102435980A (zh) * | 2011-09-15 | 2012-05-02 | 中南大学 | 一种基于解析求解的声发射源或微震源定位方法 |
CN105510959A (zh) * | 2015-11-30 | 2016-04-20 | 山东康威通信技术股份有限公司 | 一种隧道的振源类型识别方法及振源位置定位方法 |
CN105785436A (zh) * | 2016-03-17 | 2016-07-20 | 北京矿冶研究总院 | 一种矿用微震监测方法 |
CN106680871A (zh) * | 2016-12-05 | 2017-05-17 | 中南大学 | 一种基于p波到时与波形振幅的震源定位方法 |
Non-Patent Citations (1)
Title |
---|
汶川地震科学钻探3号井孔附近微震震源参数确定及其意义;叶庆东 等;《地球物理学报》;20170731;第60卷(第7期);第2716-2732页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109521221A (zh) | 2019-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xiao et al. | ISRM suggested method for in situ microseismic monitoring of the fracturing process in rock masses | |
Feng et al. | ISRM suggested method for in situ acoustic emission monitoring of the fracturing process in rock masses | |
CN109613616B (zh) | 一种浅层岩溶隧道地质综合立体探测方法 | |
CN108957521B (zh) | 一种用于隧道长距离三维超前地质预报方法 | |
CN106501848B (zh) | 一种隧道掘进过程中隐性断层超前物探方法 | |
CN202837558U (zh) | 一种地下溶洞地震跨孔ct探测及层析成像装置 | |
CN104267442A (zh) | 一种用于煤矿井下的瞬变电磁拟地震探测方法 | |
CN113189644B (zh) | 一种微震震源定位方法及*** | |
CN113153430B (zh) | 巷道围岩破坏声发射定位与波速成像监测及灾变预警方法 | |
CN109521221B (zh) | 一种钻爆法施工硬岩隧道微震波波速实时获取方法 | |
NO341717B1 (no) | Stakking av seismiske støydata for å analysere mikroseismiske hendelser | |
CN104374828A (zh) | 一种隐患探测的超声波层析成像方法 | |
CN112485823A (zh) | 高效综合超前地质预报方法 | |
CN108121010B (zh) | 基于孔巷联合的井下独头巷道槽波超前探测方法及*** | |
CN114994791B (zh) | 一种用于评估井地一体微震监测***监测能力的方法 | |
Huang et al. | Source identification of microseismic events in underground mines with interferometric imaging and cross wavelet transform | |
CN113820743A (zh) | 一种微震定位方法 | |
CN111352151B (zh) | 一种井中虚源法 | |
CN208334667U (zh) | 一种用于软岩隧道长距离三维超前地质预报装置 | |
CN103984006A (zh) | 一种全断面观测***的隧道超前地质探测方法 | |
CN112130204B (zh) | 一种适用于深埋隧洞围岩波速的现场测试方法 | |
JPH10311880A (ja) | 切羽前方探査システムおよびその方法およびトンネル掘削方法 | |
CN113050085A (zh) | 一种超前地质预报的方法 | |
CN110579799A (zh) | 一种等旅行时间间隔的地震采集观测方法及*** | |
CN111290022A (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 |