CN106501848B - 一种隧道掘进过程中隐性断层超前物探方法 - Google Patents
一种隧道掘进过程中隐性断层超前物探方法 Download PDFInfo
- Publication number
- CN106501848B CN106501848B CN201611025673.8A CN201611025673A CN106501848B CN 106501848 B CN106501848 B CN 106501848B CN 201611025673 A CN201611025673 A CN 201611025673A CN 106501848 B CN106501848 B CN 106501848B
- Authority
- CN
- China
- Prior art keywords
- microseismic
- fault
- caused
- tunneling
- events
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 62
- 230000005641 tunneling Effects 0.000 title claims abstract description 44
- 238000012544 monitoring process Methods 0.000 claims abstract description 32
- 230000004913 activation Effects 0.000 claims abstract description 30
- 238000001914 filtration Methods 0.000 claims description 6
- 238000013459 approach Methods 0.000 claims description 4
- 230000002238 attenuated effect Effects 0.000 claims description 2
- 238000012937 correction Methods 0.000 claims description 2
- 238000002474 experimental method Methods 0.000 claims description 2
- 238000010276 construction Methods 0.000 abstract description 4
- 238000001994 activation Methods 0.000 description 23
- 238000005516 engineering process Methods 0.000 description 7
- 239000011435 rock Substances 0.000 description 4
- 238000005553 drilling Methods 0.000 description 3
- 230000035939 shock Effects 0.000 description 3
- 208000013201 Stress fracture Diseases 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 206010063385 Intellectualisation Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
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
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21D—SHAFTS; TUNNELS; GALLERIES; LARGE UNDERGROUND CHAMBERS
- E21D9/00—Tunnels or galleries, with or without linings; Methods or apparatus for making thereof; Layout of tunnels or galleries
- E21D9/003—Arrangement of measuring or indicating devices for use during driving of tunnels, e.g. for guiding machines
-
- 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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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
-
- 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/65—Source localisation, e.g. faults, hypocenters or reservoirs
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Environmental & Geological Engineering (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geochemistry & Mineralogy (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种隧道掘进过程中隐性断层超前物探方法,基于微震监测系统,通过分析微震事件的空间分布与隧道掘进工作面之间的关系,判断出断层活化引起的微震事件;选取断层活化引起的微震事件进行分析,画出断层活化引起的微震事件的能量密度云图,从而判断隐性断层的倾角和走向。该方法施工简单、准确率高、预报距离较长,具有极强的经济价值和社会意义。
Description
技术领域
本发明属于隐性断层探测技术领域,具体涉及一种隧道掘进过程中隐性断层超前物探方法。
背景技术
在工程尺度下,隐性断层活化对突水、隧道(井巷)支护等都有很大的影响。因此,开展隧道(巷道)掘进过程中隐性断层超前物探方法工作对于工程断层活化引起的相关工程灾害的预警具有十分重要的意义。由于隐性断层往往规模较小,地质上不易被揭露,因此隧道(巷道)掘进过程中往往需要对隐性断层超前物探。现价段,对隐性断层超前物探方法主要有三维地震、超前钻孔和声波法。
大量实践表明,三维地震勘探技术在大型断层超前物探中的准确性较好,在小规模断层超前物探中的准确性较差。超前钻孔施工比较复杂,同时容易诱发瓦斯突出、突水等灾害,超前钻孔的长度一般在20m到50m之间,其超前预报的距离较短。声波法最大的缺点在于,近场准确性较高,超过一定距离(50m)其准确性大大降低。
近年来,微震监测技术作为一种岩体微破裂三维空间监测技术得到了迅速发展。微震主要监测岩体脆性材料在外力作用下,发生微破裂的时间、空间、强度以及相关震源参数,相对传统的位移或应力监测技术而言,微震监测技术具有如下几个重要特点:(1)监控范围广;(2)实现了监测的自动化、信息化和智能化;(3)监测仪器设备高集成性、小体积、多通道、高灵敏度;(4)支持自动监测和信息远程传输;(5)更有利于保证监测系统的长期运行。大量实践表明,微震监测技术是监测岩体微破坏及与微岩体破坏相关灾害问题的一种强有力手段。而断层在扰动应力场作用下,在活化过程中必然产生大量微震事件,因此如何将微震监测技术应用到隧道掘进过程中隐性断层超前物探中,是目前需要解决的问题。
发明内容
本发明的目的是解决上述问题,提供一种施工简单、准确性高的隧道掘进过程中隐性断层超前物探方法。
为解决上述技术问题,本发明的技术方案是:一种隧道掘进过程中隐性断层超前物探方法,包括以下步骤:
步骤一:在被监测隧道建立微震监测系统;
步骤二:根据微震监测系统监测到的实时波形数据依次进行滤波、微震事件判定、微震事件定位以及震源参数的确定;
步骤三:分析微震事件的空间分布与隧道掘进工作面之间的关系,判断哪些是隧道掘进引起的微震事件,哪些是断层活化引起的微震事件;
步骤四:选取断层活化引起的微震事件进行分析,采用微震事件定位信息和震源参数画出断层活化引起的微震事件的能量密度云图,并判断隐性断层的倾角和走向。
优选地,所述步骤一中微震监测系统包括设于已掘进的隧道一侧,距离工作面300m范围内的检波器,检波器布置在已掘进的隧道的两帮,在隧道走向方向和深度方向呈交错布置,检波器与数据采集系统相连;工作面每推进100m,隧道内布置的检波器整体向前移动100m。
优选地,所述数据采集系统依次与信号接收器、数字转化器以及数据采集与分析系统相连。
优选地,所述步骤三中,隧道掘进引起的微震事件分布于工作面前方50m以内;而当前方有断层,当工作面距离断层200m时,断层区域附近区域便开始产生微震事件,随着掘进工作面与断层距离的靠近,微震事件越来越多,直至掘进工作面通过断层200m距离,以此作为根据,判断哪些是隧道掘进引起的微震事件,哪些是断层活化引起的微震事件。
优选地,所述步骤四中,通过二维插值手段画出断层活化引起的微震事件的能量密度云图。
优选地,所述步骤四中,隐性断层引起的微震事件的能量密度云图呈带状分布,并以此判断隐性断层的倾角和走向。
本发明的有益效果是:本发明所提供的一种隧道掘进过程中隐性断层超前物探方法,基于微震监测系统,通过分析微震事件的空间分布与隧道掘进工作面之间的关系,判断出断层活化引起的微震事件;选取断层活化引起的微震事件进行分析,画出断层活化引起的微震事件的能量密度云图,从而判断隐性断层的倾角和走向。该方法施工简单、准确率高、预报距离较长,具有极强的经济价值和社会意义。
附图说明
图1是本发明隧道掘进过程中隐性断层超前物探方法的流程图;
图2是本发明检波器安装位置示意图;
图3是本发明微震监测系统的拓扑图;
图4是实施例一中掘进与断层活化引起微震事件数的关系图;
图5是实施例一中2015年2月1日断层活化引起微震事件的能量密度云图;
图6是实施例一中2015年3月1日断层活化引起微震事件的能量密度云图;
图7是实施例一中2015年4月1日断层活化引起微震事件的能量密度云图;
图8是实施例一中2015年5月1日断层活化引起微震事件的能量密度云图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步的说明:
如图1至图3所示,本发明的一种隧道掘进过程中隐性断层超前物探方法,包括以下步骤:
步骤一:在被监测隧道建立微震监测系统;微震监测系统包括设于已掘进的隧道一侧,距离工作面300m范围内的检波器,检波器布置在已掘进的隧道的两帮,在隧道走向方向和深度方向呈交错布置,检波器与数据采集系统相连;数据采集系统依次与信号接收器、数字转化器以及数据采集与分析系统相连。工作面每推进100m,隧道内布置的检波器整体向前移动100m。
步骤二:微震监测系统实时将监测到的波形数据传到地面的数据采集与分析系统,根据微震监测系统监测到的实时波形数据依次进行滤波、微震事件判定、微震事件定位以及震源参数的确定;
滤波主要是指将一些有规律性的干扰波形信号过滤掉,提高微震事件判定的效率和准确度。微震监测系统可能监测到的波形信号包括敲击实验信号、干扰信号和微震信号,这三种波形型号存在典型区别,微震信号可以看到明显的P波起震,然后衰减,紧接着S波起震,然后衰减的完整过程,微震事件持续的事件一般不超过300ms,微震主频率一般不超过200,并且至少被四个传感器接收到信号,以此作为微震事件判定的依据。
所述微震事件定位,首先确定每个传感器的微震事件波形的P波起震位置和S波起震位置,然后采用相对定位法、单纯形法及其混合定位法、震源位置和台站联合校正法、Geiger及各种改进法、最小二乘法、台偶时差法、EHB法以及双重残差法中的一种或几种方法的组合对微震事件进行定位,从而确定微震事件发生的平面坐标、深度以及发生时间;所述震源参数的确定,在微震事件发生的平面坐标、深度以及发生时间确定的条件下,通过分析波形,采用地震学的基本模型,从而确定微震事件的能量和矩震级。
步骤三:分析微震事件的空间分布与隧道掘进工作面之间的关系,隧道掘进引起的微震事件分布于工作面前方50m以内;而当前方有断层,当工作面距离断层200m时,断层区域附近区域便开始产生微震事件,随着掘进工作面与断层距离的靠近,微震事件越来越多,直至掘进工作面通过断层200m距离,以此作为根据,判断哪些是隧道掘进引起的微震事件,哪些是断层活化引起的微震事件。
步骤四:选取断层活化引起的微震事件进行分析,采用微震事件定位信息和震源参数,通过二维插值手段画出断层活化引起的微震事件的能量密度云图(包括俯视图和侧视图),隐性断层引起的微震事件的能量密度云图(包括俯视图和侧视图)呈带状分布,并以此判断隐性断层的倾角和走向。
以下通过隧道掘进过程中的具体的实施例对上述隐性断层超前物探方法作进一步说明:
实施例一
一种基于该微震监测系统的隧道掘进过程中隐性断层超前物探方法,包括以下步骤:
步骤一:如图3所示,在被监测隧道建立微震监测系统;图中S1、S2、S3、S4、S5和S6为交错布置的检波器。
步骤二:微震监测系统24小时不间断采集工作面的微震监测数据,采集过程中,应对工作面进度以及工作面地质情况做好详细记录。根据微震监测系统监测到的实时波形数据依次进行滤波、微震事件判定、微震事件定位以及震源参数的确定,监测到的微震事件数随着时间的变化如图4所示,隧道掘进进度如表1所示:
表1掘进隧道掘进进度表
从图4可以看出,微震事件从1月中旬开始产生,表1中工作面从1月中旬开始掘进,两者之间具有对应关系。图4中从4月上旬至8月下旬,构造影响微震区的微震事件数呈现出稳定状态,8月下旬以后微震事件数继续增加(由其他工作面开采引起)。
步骤三:分析微震事件的空间分布与隧道掘进工作面之间的关系,隧道掘进引起的微震事件分布于工作面前方50m以内;而当前方有断层,当工作面距离断层200m时,断层区域附近区域便开始产生微震事件,随着掘进工作面与断层距离的靠近,微震事件越来越多,直至掘进工作面通过断层200m距离,以此作为根据,判断哪些是隧道掘进引起的微震事件,哪些是断层活化引起的微震事件。
步骤四:选取断层活化引起的微震事件进行分析,采用微震事件定位信息和震源参数画出断层活化引起的微震事件的能量密度云图,如图5至图8所示,依次为2015年2月1日、2015年3月1日、2015年4月1日和2015年5月1日的断层活化引起微震事件的能量密度云图,并判断隐性断层的倾角和走向。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (2)
1.一种隧道掘进过程中隐性断层超前物探方法,其特征在于,包括以下步骤:
步骤一:在被监测隧道建立微震监测系统;
步骤二:根据微震监测系统监测到的实时波形数据依次进行滤波、微震事件判定、微震事件定位以及震源参数的确定;微震监测系统监测到的波形信号包括敲击实验信号、干扰信号和微震信号,微震信号可以看到P波起震,然后衰减,紧接着S波起震,然后衰减的完整过程,微震事件持续的事件不超过300ms,微震主频率不超过200,并且至少被四个传感器接收到信号,以此作为微震事件判定的依据;所述微震事件定位,首先确定每个传感器的微震事件波形的P波起震位置和S波起震位置,然后采用相对定位法、单纯形法及其混合定位法、震源位置和台站联合校正法、Geiger及各种改进法、最小二乘法、台偶时差法、EHB法以及双重残差法中的一种或几种方法的组合对微震事件进行定位,从而确定微震事件发生的平面坐标、深度以及发生时间;所述震源参数的确定,在微震事件发生的平面坐标、深度以及发生时间确定的条件下,通过分析波形,采用地震学的基本模型,从而确定微震事件的能量和矩震级;
步骤三:分析微震事件的空间分布与隧道掘进工作面之间的关系,判断哪些是隧道掘进引起的微震事件,哪些是断层活化引起的微震事件;
步骤四:选取断层活化引起的微震事件进行分析,采用微震事件定位信息和震源参数画出断层活化引起的微震事件的能量密度云图,并判断隐性断层的倾角和走向;
所述步骤一中微震监测系统包括设于已掘进的隧道一侧,距离工作面300m范围内的检波器,检波器布置在已掘进的隧道的两帮,在隧道走向方向和深度方向呈交错布置,检波器与数据采集系统相连;工作面每推进100m,隧道内布置的检波器整体向前移动100m;
所述步骤三中,隧道掘进引起的微震事件分布于工作面前方50m以内;而当前方有断层,当工作面距离断层200m时,断层区域附近区域便开始产生微震事件,随着掘进工作面与断层距离的靠近,微震事件越来越多,直至掘进工作面通过断层200m距离,以此作为根据,判断哪些是隧道掘进引起的微震事件,哪些是断层活化引起的微震事件;
所述步骤四中,通过二维插值手段画出断层活化引起的微震事件的能量密度云图;
所述步骤四中,隐性断层引起的微震事件的能量密度云图呈带状分布,并以此判断隐性断层的倾角和走向。
2.根据权利要求1所述的隧道掘进过程中隐性断层超前物探方法,其特征在于:所述数据采集系统依次与信号接收器、数字转化器以及数据采集与分析系统相连。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611025673.8A CN106501848B (zh) | 2016-11-15 | 2016-11-15 | 一种隧道掘进过程中隐性断层超前物探方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611025673.8A CN106501848B (zh) | 2016-11-15 | 2016-11-15 | 一种隧道掘进过程中隐性断层超前物探方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106501848A CN106501848A (zh) | 2017-03-15 |
CN106501848B true CN106501848B (zh) | 2020-11-10 |
Family
ID=58327569
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611025673.8A Expired - Fee Related CN106501848B (zh) | 2016-11-15 | 2016-11-15 | 一种隧道掘进过程中隐性断层超前物探方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106501848B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108415066B (zh) * | 2017-12-13 | 2020-03-24 | 北方工业大学 | 一种隧道施工地质灾害预报方法 |
CN108194142A (zh) * | 2017-12-25 | 2018-06-22 | 河北煤炭科学研究院 | 隐性裂隙带探测方法 |
CN108490485B (zh) * | 2018-06-20 | 2024-04-05 | 成都理工大学 | 双线隧道微震监测系统及其定位精度评估方法 |
CN108931822B (zh) * | 2018-09-17 | 2020-08-28 | 贵州省水利水电勘测设计研究院有限公司 | 一种地下全填充溶洞勘探方法 |
CN109441455B (zh) * | 2019-01-12 | 2020-06-30 | 黄捷军 | 一种隧道工程安全施工开挖方法 |
CN109779635B (zh) * | 2019-02-02 | 2021-01-05 | 中建铁投轨道交通建设有限公司 | 一种隧道工程安全施工开挖方法 |
CN111397593B (zh) * | 2020-02-28 | 2023-09-22 | 中国海洋大学 | 一种轨道交通装备导航数据处理系统及方法 |
CN111636859B (zh) * | 2020-07-09 | 2022-08-16 | 中煤科工集团重庆研究院有限公司 | 基于微破裂波检测的煤岩随钻自识别方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102506993A (zh) * | 2011-11-21 | 2012-06-20 | 大同煤矿集团有限责任公司 | 煤矿井下围岩微震检测方法 |
CN104153814A (zh) * | 2014-06-23 | 2014-11-19 | 中国矿业大学 | 一种冲击地压的微震多参量预警方法 |
CN105652311A (zh) * | 2016-01-18 | 2016-06-08 | 陕西煤业化工技术研究院有限责任公司 | 一种监测底板突水的微震监测方法 |
-
2016
- 2016-11-15 CN CN201611025673.8A patent/CN106501848B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN106501848A (zh) | 2017-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106501848B (zh) | 一种隧道掘进过程中隐性断层超前物探方法 | |
Parker et al. | Active‐source seismic tomography at the Brady geothermal field, Nevada, with dense nodal and fiber‐optic seismic arrays | |
AU2007272702B2 (en) | Time lapse marine seismic surveying employing interpolated multicomponent streamer pressure data | |
Picozzi et al. | Interferometric analysis of strong ground motion for structural health monitoring: the example of the L’Aquila, Italy, seismic sequence of 2009 | |
Zollo et al. | Bayesian estimation of 2-DP-velocity models from active seismic arrival time data: imaging of the shallow structure of Mt Vesuvius (Southern Italy) | |
CN103336297B (zh) | 微破裂向量扫描方法 | |
Bindi et al. | Site effects observed in alluvial basins: the case of Norcia (Central Italy) | |
Papadimitriou et al. | The Santorini Volcanic Complex: A detailed multi-parameter seismological approach with emphasis on the 2011–2012 unrest period | |
CN105765408A (zh) | 用于分析位于地下矿场巷道上方的层的地质结构以及相对应力变化的方法及系统 | |
Stabile et al. | A comprehensive approach for evaluating network performance in surface and borehole seismic monitoring | |
Naoi et al. | Twenty thousand aftershocks of a very small (M 2) earthquake and their relation to the mainshock rupture and geological structures | |
Satoh et al. | Three-dimensional finite-difference waveform modeling of strong motions observed in the Sendai basin, Japan | |
Juhlin et al. | Reflection seismic studies over the end-glacial Burträsk fault, Skellefteå, Sweden | |
Moriya et al. | Delineation of large localized damage structures forming ahead of an active mining front by using advanced acoustic emission mapping techniques | |
Ugalde et al. | Passive seismic monitoring of an experimental CO2 geological storage site in Hontomín (Northern Spain) | |
Langet et al. | Joint focal mechanism inversion using downhole and surface monitoring at the Decatur, Illinois, CO2 injection site | |
Odum et al. | Multisource, high-resolution seismic-reflection imaging of Meeman-Shelby fault and a possible tectonic model for a Joiner Ridge–Manila High stepover structure in the upper Mississippi Embayment region | |
CN104749630B (zh) | 构建微地震监测速度模型的方法 | |
Zeng et al. | Turning a telecom fiber‐optic cable into an ultradense seismic array for rapid postearthquake response in an urban area | |
D’Alessandro et al. | Montana Seismic Network Performance: an evaluation through the SNES method | |
Xu et al. | Optimal design of microseismic monitoring networking and error analysis of seismic source location for rock slope | |
Ahmadi et al. | Revealing the deeper structure of the end-glacial Pärvie fault system in northern Sweden by seismic reflection profiling | |
Jousset et al. | Seismic tomography in Reykjanes, SW Iceland | |
Quiros et al. | Aftershock imaging with dense arrays (AIDA) following the M w 4.0 Waterboro earthquake of 16 October 2012 Maine, USA | |
Loginov et al. | Joint processing of surface and underground microseismic monitoring data in hard mineral mining |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |
Granted publication date: 20201110 Termination date: 20201115 |
|
CF01 | Termination of patent right due to non-payment of annual fee |