CN105257339B - 掘进工作面多参量综合监测预警方法 - Google Patents
掘进工作面多参量综合监测预警方法 Download PDFInfo
- Publication number
- CN105257339B CN105257339B CN201510786271.9A CN201510786271A CN105257339B CN 105257339 B CN105257339 B CN 105257339B CN 201510786271 A CN201510786271 A CN 201510786271A CN 105257339 B CN105257339 B CN 105257339B
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- monitoring
- stress
- mtr
- 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
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
- E21F17/18—Special adaptations of signalling or alarm devices
Abstract
本发明公开了一种掘进工作面多参量综合监测预警方法,属于煤矿动力灾害监测预警领域,所述方法包括:步骤1:将掘进工作面分区监测;步骤2:在掘进巷道内安设煤体应力、地音和微震、锚杆锚索应力监测系统中的至少两种监测系统;步骤3:获取各监测预警指标;步骤4:将各监测预警指标分别赋予一定的权重并累加,得到综合预警指标。本发明根据掘进工作面冲击地压发生机理,提出了掘进工作面分区监测的思想,将各监测指标在不同分区内的变化情况耦合,以多参量综合监测指标为预警参数,实现掘进工作面冲击危险性的综合预警。本发明解决了现有技术中仅依靠掘进工作面冲击地压单指标监测的方法无法达到预测预警掘进工作面冲击地压危险的问题。
Description
技术领域
本发明涉及煤矿动力灾害监测预警领域,特别是指一种掘进工作面多参量综合监测预警方法。
背景技术
近年来,掘进工作面冲击地压事故在矿山开采中屡见不鲜,冲击范围有些达到100m以上,致使巷道设备损毁,人员伤亡,对煤矿造成了巨大的损失。由于煤岩体的流变特性,工作面掘进后需要经历一定的时间段,才能形成相对稳定的区域。实际过程中,由于矿山的快速掘进,采掘空间形成后,煤岩体内地质构造的不均匀性和煤体应力分布的不均衡性,会使工作面前方煤体内集聚的应力来不及充分释放,新揭露的煤体还未形成卸压稳定的状态,致使煤体失稳,造成冲击地压事故。
对于冲击地压灾害的监测预警,目前在回采工作面较为成熟,掘进工作面仍相对薄弱。国内外已经提出了钻屑监测法、煤体应力监测法、电磁辐射监测法、地音和微震监测法等多种方法。
1)钻屑监测法:通过监测煤层钻孔排粉量变化规律以及相关的动力效应,了解煤体应力状态,以此预测冲击地压危险性,是目前最常使用的监测方法。该方法具有因人而异的操作误差和不能连续监测的缺点;
2)煤体应力监测法:通过对煤体内采动应力进行连续监测,从应力场角度对煤岩体冲击危险性进行预测和评价,该方法实现了煤体采动应力变化量的连续监测。该方法监测范围较小,且监测结果对于自发型冲击地压的监测效果较好,对于诱发型冲击地压监测效果较差;
3)电磁辐射监测法:对煤岩体破裂过程中向外辐射的电磁强度和脉冲数进行监测,以此判断煤岩体受载程度和破裂强度,得到冲击危险程度。该方法受井下各种电信号的影响,所得结果具有不确定性;
4)地音和微震监测法:对煤岩体破裂过程中释放的震动信号进行监测,可分析得到煤岩体不同破裂阶段,了解煤岩体整体破坏情况和能量释放情况,实现冲击的预测和预警。该方法只监测破裂产生的震动信号,且对于诱发型冲击地压的监测效果较好,对于自发型冲击地压监测效果较差。
当前综合机械化掘进速度普遍较快,现场围岩破裂信息量不足,无法掌握掘进巷道围岩活动区与围岩稳定区的分布情况,无法了解围岩破裂程度及活动趋势;掘进工作面开挖空间小,应力扰动区域相对较小,应力转移范围有限,可以监测到的应力场变化幅度不大;巷道掘进一段距离后,无法实时了解围岩支护情况,缺少对支护结构工作状态、支承压力及围岩安全性的实时评价;各种监测方法的预警指标及判据相对独立,无法统一表征巷道的冲击危险性,也很难从矿压理论上揭示各指标之间存在的内在联系。
对于冲击地压灾害的预警,现有监测方法都有各自本身的缺点;另外,现有单个监测方法只对冲击地压发生过程中某一种现象进行监测,仅用单一参量监测指标,很难达到准确预警的目的。鉴于此,有些学者提出了灰色关联度分析法、综合指数法、模糊综合评价法、神经网络分析法等综合评价冲击危险性,但这些方法的基础是模糊数学,所得到的结论均属于宏观评价,对于局部区域危险程度的判定误差较大。
发明内容
本发明要解决的技术问题是提供一种准确性高,能够实现掘进工作面冲击危险性综合预警的掘进工作面多参量综合监测预警方法。
为解决上述技术问题,本发明提供技术方案如下:
一种掘进工作面多参量综合监测预警方法,包括:
步骤1:掘进工作面分区监测
将掘进工作面划分为两个危险区域:动态演化区和蠕变隐患区,其中,动态演化区是自掘进工作面至后方60m范围,蠕变隐患区是自掘进工作面后方60m至后方260m范围;
步骤2:布置测点
在掘进巷道内安设煤体应力监测系统、地音和微震监测系统、锚杆锚索应力监测系统中的至少两种监测系统,其中,煤体应力测点自掘进工作面后方25m处开始布置,布设间距为25m,共布置8-12组;地音和微震监测系统传感器自掘进工作面后方30m处开始布置,布设间距为30m,共布置6-10组;锚杆锚索应力测点自掘进工作面后方50m处开始布置,布设间距为50m,共布置4-6组;随着掘进工作面向前推进,各监测系统的测点相应也需要进行前移;
步骤3:获取煤体应力监测预警指标、地音和微震监测预警指标、锚杆锚索应力监测预警指标;
步骤4:将所获取的煤体应力监测预警指标、地音和微震监测预警指标、锚杆锚索应力监测预警指标分别赋予一定的权重并累加,得到所述动态演化区的综合预警指标和蠕变隐患区的综合预警指标。
本发明具有以下有益效果:
本发明根据掘进工作面冲击地压发生机理,以长期现场监测为基础,提出了掘进工作面分区监测的思想,将各监测指标在不同分区内的变化情况耦合,以多参量综合监测指标为预警参数,实时在线监测掘进工作面冲击地压危险性,实现掘进工作面冲击危险性的综合预警。本发明解决了现有技术中仅依靠掘进工作面冲击地压单指标监测的方法无法达到预测预警掘进工作面冲击地压危险的问题。
附图说明
图1为本发明的掘进工作面多参量综合监测系统架构示意图;
图2为本发明中的测点布置和移组方案示意图;
图3为本发明中的煤体应力增幅预警指标隶属图;
图4为本发明中的掘进工作面多参量综合监测预警系统运算流程示意图;
图5为本发明中的掘进巷道冲击危险性实时综合监测预警判别流程示意图。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚,下面将结合附图及具体实施例进行详细描述。
本发明提供一种掘进工作面多参量综合监测预警方法,监测子参量由主要子监测参量(煤体应力监测参量、地音和微震监测参量)和辅助子监测参量(锚杆锚索应力监测参量)组成,系统架构可以参考图1来设置。
在图1中,11为监测数据服务器,12为实时监控,13为数据分析服务器,14为光纤环网,15为顶板,16为锚杆锚索应力测点,17为地音传感器,18为迎头(即掘进工作面),19为钻孔应力计,20为无线应力测点,21为无线数据采集子站,22为底板,23为监测分站,24为光缆,25为127V供电电缆,26为电缆,27为CAN总线。
需要说明的是,图1中的系统架构仅作参考,本领域技术人员在本发明的原理的基础上可以采用各种可行的架构方案。
本发明具体实现步骤包括:
步骤1:掘进工作面分区监测
根据掘进工作面冲击地压发生机理、围岩应力与支护结构耦合作用,将掘进工作面划分为两个危险区域:动态演化区和蠕变隐患区。
第一个为动态演化区,自掘进工作面至后方60m范围,此区域处于巷道开挖初期应力调整和重新分布阶段,煤岩体破裂及应力重新分布过程中产生的地音信号和煤体应力信号较其它监测信号变化更为显著;
第二个为蠕变隐患区,自掘进工作面后方60m至后方260m范围,此区域围岩应力与支护结构逐渐趋于相对稳定状态,巷道在高应力环境作用下仍会发生蠕变,煤体应力信号和锚杆锚索监测信号较其它监测信号变化更为显著。
步骤2:布置测点
在掘进巷道内安设煤体应力监测系统、地音和微震监测系统、锚杆锚索应力监测系统中的至少两种监测系统。
其中,煤体应力测点自掘进工作面后方25m处开始布置,布设间距为25m,共布置8-12组,优选布置10组;
地音和微震监测系统传感器自掘进工作面后方30m处开始布置,布设间距为30m,共布置6-10组,优选布置8组;
锚杆锚索应力测点自掘进工作面后方50m处开始布置,布设间距为50m,共布置4-6组,优选布置5组;
各监测系统优选布置在监测掘进工作面及后方260m范围,包括动态演化区及蠕变隐患区,随着掘进工作面向前推进,各监测系统的测点相应也需要进行移组,测点布置和移组方案可参考图2。
在图2中,1为锚杆锚索应力测点,2为地音传感器,3为迎头(即掘进工作面),4为深孔应力测点,5为浅孔应力测点,6为无线数据采集子站,7为光缆,8为电缆,9为CAN总线,10为移组的应力测点,11为移组的传感器,12为监测分站。
步骤3:获取煤体应力监测预警指标、地音和微震监测预警指标、锚杆锚索应力监测预警指标;
本步骤中,煤体应力监测预警指标、地音和微震监测预警指标、锚杆锚索应力监测预警指标可按照现有技术中的常规方式计算得到,也可按照本发明后续给出的方式计算得到。
步骤4:将所获取的煤体应力监测预警指标、地音和微震监测预警指标、锚杆锚索应力监测预警指标分别赋予一定的权重并累加,得到所述动态演化区的综合预警指标和蠕变隐患区的综合预警指标。
本发明中,掘进工作面分区域监测实现了根据不同区域内冲击机理的不同,合理选择了各监测系统的布设位置,并将各监测子参量按一定权重累加,充分发挥了各监测系统的优势,提高预警结果的科学性和准确性。
本发明根据掘进工作面冲击地压发生机理,以长期现场监测为基础,提出了掘进工作面分区监测的思想,将各监测指标在不同分区内的变化情况耦合,以多参量综合监测指标为预警参数,实时在线监测掘进工作面冲击地压危险性,实现掘进工作面冲击危险性的综合预警。本发明解决了现有技术中仅依靠掘进工作面冲击地压单指标监测的方法无法达到预测预警掘进工作面冲击地压危险的问题。
本发明能够对掘进巷道冲击地压事故进行有效的预警,提高现场管理和施工人员的工作效率,指导危险区卸压。
本发明中,各监测预警指标优选采用发明人优化、升级后的如下方法计算得到。具体的:
1、对于步骤3中的煤体应力监测预警指标,优选通过以下步骤获得:
步骤31:计算煤体应力值大小指标IY;
步骤32:计算煤体应力增幅预警指标IDY;
步骤33:计算煤体应力监测预警指标I1=a1·IY+b1·IDY,其中(a1,b1)∈[0.4,0.6]且a1+b1=1。
上述步骤31-32可以采用现有技术中的常规计算方法计算得到煤体应力值大小指标IY、煤体应力增幅预警指标IDY。上述步骤33实现了IY和IDY两者较佳的耦合。
优选的,步骤31具体为:
其中,Y为监测区域内最大应力测点数值;
煤体应力值大小指标IY采用区域内最大应力测点数值,作为总应力预警指标。并与设定的临界值对比,确定危险性。
例如:动态演化区内有2个应力测点,一个班8小时内的最大应力测点数值为8.5MPa,而设定的测点预警值如下表:
则本班动态演化区总应力预警等级为中等冲击危险。
本例中:IY=0.5Y-3.5=0.75。
优选的,步骤32具体为:
其中,α为系数,P为临界值,两者一一对应,P=1/α。α的取值与应力测点稳定初装力值有关,取值关系如下表所示。
DY为区域内煤体应力最大应力测点的增幅值,由当前监测区域内最大应力测点的应力与测点稳定初装力之差除以测点稳定初装力得到,具体的,DY=100(F’-F)/F,F’为当前监测区域内最大应力测点的应力,F为测点稳定初装力。
煤体应力增幅预警指标是在总应力达到某个临界值后才起作用,总应力小时不起作用。因此,根据不同的“稳定初装应力”,设定不同的增幅预警隶属函数。
例如:动态演化区测点稳定初装力为4.5MPA,如图3所示,可以选择上述4MPA的曲线,8小时内如果应力增加到6MPA,则DY=100(6-4.5)/4.5=33.3(实际上是33.3%),代入上式,得到:
IDY=0.01×33.3=0.333。
本例中,对于煤体应力监测预警指标I1,假设a1=b1=0.5,则I1=0.5(IY+IDY)=0.5(0.75+0.333)=0.54。
2、对于步骤3中的地音和微震监测预警指标,优选通过以下步骤获得:
步骤31’:计算地音和微震事件数异常率指标INb;
步骤32’:计算地音和微震能量异常率指标IEb;
步骤33’:计算地音和微震监测预警指标I2=a2·IEb+b2·INb,其中(a2,b2)∈[0.4,0.6]且a2+b2=1。
上述步骤31’-32’可以采用现有技术中的常规计算方法计算得到地音和微震事件数异常率指标INb、地音和微震能量异常率指标IEb。上述步骤33’实现了INb和IEb两者较佳的耦合。
优选的,步骤31’具体为:
其中,Nb=100(N’-N)/N,Nb为地音和微震事件数异常率比值,N为正常阶段一段时间内的有效事件数,N’为当前相同长度时间段内的有效事件数;
优选的,步骤32’具体为:
其中,Eb=100(E’-E)/E,Eb为地音和微震能量异常率比值,E为正常阶段一段时间内的能量之和,E’为当前相同长度时间段内的能量之和。
例如:采用某区域内的地音或微震探头数据,计算每8小时(一个班)的有效事件数和能量之和,并与正常阶段(一般为前7天、相同班次)的平均数相比,得到当班的异常率比值。通过一段时间的监测,得到比较准确的“正常”有效事件数和能量和的值,作为掘进工作面本监测区的标准值。
设正常阶段(一般为前7天、相同班次)的有效事件平均数为N,能量之和为E,当班的有效事件数为N’,能量之和为E’,则地音和微震事件数异常率比值为:
Nb=100(N’-N)/N;
地音和微震能量异常率比值为:
Eb=100(E’-E)/E;
在地音和微震监测表达冲击危险性中,能量的权重略大于事件次数。
进一步例如:前7天早班动态演化区事件数N平均为20个,能量和E平均为50000,本早班动态演化区事件数N’为30个,能量和E’为60000。则事件数异常率比值为:Nb=100(N’-N)/N=50;能量异常率比值为:Eb=100(E’-E)/E=20;Nb和Eb的计算结果分别为:INb=0.01Nb=0.5,IEb=0.01Eb=0.2。
取a2=0.6,取b2=0.4,则本班动态演化区的地音和微震监测预警指标I2为:
I2=0.6IEb+0.4INb=0.6*0.5+0.4*0.2=0.38。
3、对于步骤3中的锚杆锚索应力监测预警指标,优选通过以下步骤获得(与上述煤体应力监测预警指标的算法相同):
步骤31”:计算锚杆锚索应力值大小指标I”Y;
步骤32”:计算锚杆锚索应力增幅预警指标I”DY;
步骤33”:计算锚杆锚索应力监测预警指标I3=a3·I”Y+b3·I”DY,其中(a3,b3)∈[0.4,0.6]且a3+b3=1。
上述步骤31”-32”可以采用现有技术中的常规计算方法计算得到锚杆锚索应力值大小指标I”Y、锚杆锚索应力增幅预警指标I”DY。上述步骤33”实现了I”Y和I”DY两者较佳的耦合。
优选的,步骤31”具体为:
其中,Y为监测区域内最大应力测点数值;
优选的,步骤32”具体为:
其中,α为系数,P为临界值,两者一一对应,P=1/α。α的取值与应力测点稳定初装力值有关,取值关系如下表所示。
DY为区域内锚杆锚索应力最大值测点的增幅值,由当前监测区域内最大应力测点的应力与测点稳定初装力之差除以测点稳定初装力得到,具体的,DY=100(F’-F)/F,F’为当前监测区域内最大应力测点的应力,F为测点稳定初装力。
此监测参量数据处理方法按照煤体应力监测指标的数据处理方法,处理锚杆锚索应力值大小指标和应力增幅预警指标,最后得到掘进工作面某监测区锚杆锚索应力监测预警指标I3。
例如,假设锚杆锚索应力值大小指标为0.4,应力增幅预警指标为0.5,取a3=b3=0.5,则I3=a3·0.4+b3·0.5=0.45。
4、对于步骤4中的综合预警指标,优选通过以下步骤获得:
I0=K1I1+K2I2+K3I3,
其中,I0为动态演化区或蠕变隐患区的综合预警指标,I1为该区对应的煤体应力监测预警指标,K1为该区对应的煤体应力监测预警指标权重系数,I2为该区对应的地音和微震监测预警指标,K2为该区对应的地音和微震监测预警指标权重系数,I3为该区对应的锚杆锚索应力监测预警指标,K3为该区对应的锚杆锚索应力监测预警指标权重系数;
根据掘进工作面不同分区内冲击地压发生机理的差异及顶板性质差异,各权重系数取值如下表:
另外,由于监测现场环境复杂,可能存在某子监测系统无监测数据的情况。
煤体应力监测系统监测数据缺失时,各权重系数取值如下表:
地音和微震监测系统监测数据缺失时,各权重系数取值如下表:
锚杆锚索应力监测系统监测数据缺失时,各权重系数取值如下表:
I0计算完成后,对比危险性等级对照表(下表),可得多参量综合监测危险等级,整体运算流程图可参考图4。
根据以上实例,假设该工作面动态演化区顶板性质属于坚硬顶板,则:I0=0.35·I1+0.4·I2+0.25·I3=0.35×0.54+0.4×0.38+0.25×0.45=0.4535。
经对照,运用掘进工作面多参量综合监测预警方法,得到掘进面动态演化区多参量综合危险等级I0=0.4535,属于弱冲击危险等级。
本发明中,为进一步增强监测预警结果的可靠性,步骤4之后还可以包括:
步骤5:根据所得到的综合预警指标,实时判断巷道是否处于危险阶段及危险程度,同时确定出现场需要钻屑法检验的危险区域,依据现场实测钻屑指标对上述巷道的综合预警指标进行验证,最终确定现场危险特征。
由此可以给出现场监测系统的具体输出参量及综合判别过程如图5所示。本步骤5中,钻屑法为本领域中公知的监测方法,此处不再赘述。
综上,通过本发明方法,可以实现煤矿掘进作业中冲击危险性实时判定、蠕变型冲击危险性判定、支护质量与围岩活动趋势评定、钻屑法实施具体参数(时间间隔、钻孔间距和区域范围)的确定、巷道冲击危险区划分。
(1)实时预警掘进工作面冲击危险性:通过监测到的围岩破裂变化程度和趋势、围岩应力变化量和趋势、锚杆(索)支护力变化趋势给出各自在监测位置的危险程度,耦合各参量后给出掘进中是否具有冲击危险性。
(2)判定蠕变型冲击危险性:分析地音信号有效监测时间内频次和平均能量、钻孔应力有效监测时间内应力增量、锚杆索应力变化量是否趋于稳定,通过综合分析其各参量的危险程度及权重判断巷道是否会发生蠕变型冲击。
(3)确定钻屑法实施的具体参量(时间间隔、钻孔间距和区域范围):通过综合判断巷道围岩的稳定性、煤岩体应力变化趋势和支护结构的稳定性,初步划分巷道冲击危险区域和程度,根据此参数确定现场是否需要进行煤粉钻检验以及检测参量。
(4)划分巷道冲击危险区:对比分析综合判定的掘进巷道冲击危险性和煤粉钻检测值,最终划定巷道冲击危险区范围,为后期回采防冲提供科学依据。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (8)
1.一种掘进工作面多参量综合监测预警方法,其特征在于,包括:
步骤1:掘进工作面分区监测
将掘进工作面划分为两个危险区域:动态演化区和蠕变隐患区,其中,动态演化区是自掘进工作面至后方60m范围,蠕变隐患区是自掘进工作面后方60m至后方260m范围;
步骤2:布置测点
在掘进巷道内安设煤体应力监测系统、地音和微震监测系统、锚杆锚索应力监测系统中的至少两种监测系统,其中,煤体应力测点自掘进工作面后方25m处开始布置,布设间距为25m,共布置8-12组;地音和微震监测系统传感器自掘进工作面后方30m处开始布置,布设间距为30m,共布置6-10组;锚杆锚索应力测点自掘进工作面后方50m处开始布置,布设间距为50m,共布置4-6组;随着掘进工作面向前推进,各监测系统的测点相应也需要进行前移;
步骤3:获取煤体应力监测预警指标、地音和微震监测预警指标、锚杆锚索应力监测预警指标;
步骤4:将所获取的煤体应力监测预警指标、地音和微震监测预警指标、锚杆锚索应力监测预警指标分别赋予一定的权重并累加,得到所述动态演化区的综合预警指标和蠕变隐患区的综合预警指标;
所述步骤3中的煤体应力监测预警指标通过以下步骤获得:
步骤31:计算煤体应力值大小指标IY;
步骤32:计算煤体应力增幅预警指标IDY;
步骤33:计算煤体应力监测预警指标I1=a1·IY+b1·IDY,其中(a1,b1)∈[0.4,0.6]且a1+b1=1;
所述步骤31具体为:
<mrow>
<msub>
<mi>I</mi>
<mi>Y</mi>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>Y</mi>
<mo><</mo>
<mn>7</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.5</mn>
<mi>Y</mi>
<mn>3.5</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>7</mn>
<mo>&le;</mo>
<mi>Y</mi>
<mo>&le;</mo>
<mn>9</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>Y</mi>
<mo>></mo>
<mn>9</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,Y为监测区域内最大应力测点数值;
所述步骤32具体为:
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>D</mi>
<mi>Y</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>D</mi>
<mi>Y</mi>
<mo>&le;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&alpha;</mi>
<mo>&CenterDot;</mo>
<mi>D</mi>
<mi>Y</mi>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo><</mo>
<mi>D</mi>
<mi>Y</mi>
<mo><</mo>
<mi>P</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>D</mi>
<mi>Y</mi>
<mo>&GreaterEqual;</mo>
<mi>P</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,α为系数,P为临界值,两者一一对应,P=1/α,α的取值与应力测点稳定初装力值有关,取值关系如下表所示:
DY为区域内煤体应力最大应力测点的增幅值,由当前监测区域内最大应力测点的应力与测点稳定初装力之差除以测点稳定初装力得到,具体的,DY=100(F’-F)/F,F’为当前监测区域内最大应力测点的应力,F为测点稳定初装力。
2.根据权利要求1所述的方法,其特征在于,所述步骤3中的地音和微震监测预警指标通过以下步骤获得:
步骤31’:计算地音和微震事件数异常率指标INb;
步骤32’:计算地音和微震能量异常率指标IEb;
步骤33’:计算地音和微震监测预警指标I2=a2·IEb+b2·INb,其中(a2,b2)∈[0.4,0.6]且a2+b2=1。
3.根据权利要求2所述的方法,其特征在于,
所述步骤31’具体为:
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>N</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>N</mi>
<mi>b</mi>
<mo>&le;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.01</mn>
<mi>N</mi>
<mi>b</mi>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo><</mo>
<mi>N</mi>
<mi>b</mi>
<mo><</mo>
<mn>100</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>N</mi>
<mi>b</mi>
<mo>&GreaterEqual;</mo>
<mn>100</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,Nb=100(N’-N)/N,Nb为地音和微震事件数异常率比值,N为正常阶段一段时间内的有效事件数,N’为当前相同长度时间段内的有效事件数;
所述步骤32’具体为:
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>E</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>E</mi>
<mi>b</mi>
<mo>&le;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.01</mn>
<mi>E</mi>
<mi>b</mi>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo><</mo>
<mi>E</mi>
<mi>b</mi>
<mo><</mo>
<mn>100</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>E</mi>
<mi>b</mi>
<mo>&GreaterEqual;</mo>
<mn>100</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,Eb=100(E’-E)/E,Eb为地音和微震能量异常率比值,E为正常阶段一段时间内的能量之和,E’为当前相同长度时间段内的能量之和。
4.根据权利要求3所述的方法,其特征在于,所述步骤3中的锚杆锚索应力监测预警指标通过以下步骤获得:
步骤31”:计算锚杆锚索应力值大小指标I”Y;
步骤32”:计算锚杆锚索应力增幅预警指标I”DY;
步骤33”:计算锚杆锚索应力监测预警指标I3=a3·I”Y+b3·I”DY,其中(a3,b3)∈[0.4,0.6]且a3+b3=1。
5.根据权利要求4所述的方法,其特征在于,
所述步骤31”具体为:
<mrow>
<msub>
<msup>
<mi>I</mi>
<mrow>
<mo>,</mo>
<mo>,</mo>
</mrow>
</msup>
<mi>Y</mi>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>Y</mi>
<mo><</mo>
<mn>50</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.1</mn>
<mi>Y</mi>
<mo>-</mo>
<mn>5</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>50</mn>
<mo>&le;</mo>
<mi>Y</mi>
<mo>&le;</mo>
<mn>60</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>Y</mi>
<mo>></mo>
<mn>60</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,Y为监测区域内最大应力测点数值;
所述步骤32”具体为:
<mrow>
<msub>
<msup>
<mi>I</mi>
<mrow>
<mo>,</mo>
<mo>,</mo>
</mrow>
</msup>
<mrow>
<mi>D</mi>
<mi>Y</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>D</mi>
<mi>Y</mi>
<mo>&le;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&alpha;</mi>
<mo>&CenterDot;</mo>
<mi>D</mi>
<mi>Y</mi>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo><</mo>
<mi>D</mi>
<mi>Y</mi>
<mo><</mo>
<mi>P</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>D</mi>
<mi>Y</mi>
<mo>&GreaterEqual;</mo>
<mi>P</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
2
其中,α为系数,P为临界值,两者一一对应,P=1/α,α的取值与应力测点稳定初装力值有关,取值关系如下表所示:
DY为区域内锚杆锚索应力最大值测点的增幅值,由当前监测区域内最大应力测点的应力与测点稳定初装力之差除以测点稳定初装力得到,具体的,DY=100(F’-F)/F,F’为当前监测区域内最大应力测点的应力,F为测点稳定初装力。
6.根据权利要求5所述的方法,其特征在于,
所述步骤4进一步为:
I0=K1I1+K2I2+K3I3,
其中,I0为动态演化区或蠕变隐患区的综合预警指标,I1为煤体应力监测预警指标,K1为煤体应力监测预警指标权重系数,I2为地音和微震监测预警指标,K2为地音和微震监测预警指标权重系数,I3为锚杆锚索应力监测预警指标,K3为锚杆锚索应力监测预警指标权重系数;
根据掘进工作面不同分区内冲击地压发生机理的差异及顶板性质差异,各权重系数取值如下表:
煤体应力监测系统监测数据缺失时,各权重系数取值如下表:
地音和微震监测系统监测数据缺失时,各权重系数取值如下表:
锚杆锚索应力监测系统监测数据缺失时,各权重系数取值如下表:
7.根据权利要求6所述的方法,其特征在于,
I0取值所对应的危险性等级如下表:
8.根据权利要求1-7中任一所述的方法,其特征在于,所述步骤4之后还包括:
步骤5:根据所得到的综合预警指标,实时判断巷道是否处于危险阶段及危险程度,同时确定出现场需要钻屑法检验的危险区域,依据现场实测钻屑指标对上述巷道的综合预警指标进行验证,最终确定现场危险特征。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510786271.9A CN105257339B (zh) | 2015-11-16 | 2015-11-16 | 掘进工作面多参量综合监测预警方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510786271.9A CN105257339B (zh) | 2015-11-16 | 2015-11-16 | 掘进工作面多参量综合监测预警方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105257339A CN105257339A (zh) | 2016-01-20 |
CN105257339B true CN105257339B (zh) | 2017-10-13 |
Family
ID=55097228
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510786271.9A Active CN105257339B (zh) | 2015-11-16 | 2015-11-16 | 掘进工作面多参量综合监测预警方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105257339B (zh) |
Families Citing this family (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105759009B (zh) * | 2016-01-28 | 2017-08-01 | 山东大学 | 一种基于微震监测的地下工程岩体锚注效果的检测方法 |
CN105626150B (zh) * | 2016-02-04 | 2017-12-08 | 山东大学 | 一种基于微震监测的沿空掘进巷道动态监测与稳定评价方法 |
CN105759010B (zh) * | 2016-02-04 | 2017-11-07 | 山东大学 | 一种采动影响巷道动态监测与稳定评价方法 |
US10884154B2 (en) | 2018-02-26 | 2021-01-05 | University Of Science And Technology Beijing | Monitoring and forewarning method for coal-rock dynamic disasters based on electromagnetic radiation and earth sound |
CN109281673B (zh) * | 2018-10-11 | 2021-01-26 | 中国铁建重工集团股份有限公司 | 一种掘进机远程监控系统的预警方法及系统 |
CN109556775A (zh) * | 2018-11-08 | 2019-04-02 | 魏全德 | 一种监测冲击地压危险性的方法 |
CN110043242A (zh) * | 2019-03-11 | 2019-07-23 | 北京安科兴业科技股份有限公司 | 一种基于应力高频采集的冲击地压监测预警方法 |
CN109991313B (zh) * | 2019-03-11 | 2024-04-19 | 北京安科兴业科技股份有限公司 | 声发射监测装置以及利用该装置的声发射监测方法 |
CN110043316A (zh) * | 2019-04-15 | 2019-07-23 | 北京安科兴业科技股份有限公司 | 一种矿山复合动力灾害综合监测预警方法 |
CN110107358A (zh) * | 2019-05-01 | 2019-08-09 | 山东科技大学 | 一种巷道掘进期间的大范围岩爆危险性预测方法 |
CN110596750B (zh) * | 2019-09-19 | 2021-01-26 | 东北大学 | 一种煤矿井下锚杆预紧安装微震传感器的装配器 |
CN110779574B (zh) * | 2019-10-30 | 2020-11-13 | 北京科技大学 | 一种煤岩动力灾害多系统多参量集成综合预警方法及系统 |
CN111220305B (zh) * | 2020-01-16 | 2021-09-07 | 山东爱拓软件开发有限公司 | 一种煤岩体应力增幅监测方法及系统 |
CN111476420B (zh) * | 2020-04-08 | 2023-07-11 | 中煤能源研究院有限责任公司 | 一种微震监测冲击地压预警指标优化方法 |
CN111648826B (zh) * | 2020-05-26 | 2022-05-13 | 陕西彬长孟村矿业有限公司 | 一种煤矿冲击地压预测预警系统及方法 |
CN112231936B (zh) * | 2020-11-11 | 2023-11-14 | 中国矿业大学 | 一种多参量影响下沿空留巷稳定性关键参数的确定方法 |
CN112443357A (zh) * | 2020-11-25 | 2021-03-05 | 攀枝花学院 | 矿山井巷围岩稳定性监测智能预警方法 |
CN112963202B (zh) * | 2021-02-05 | 2023-05-12 | 中煤科工开采研究院有限公司 | 冲击地压监测预警方法及系统 |
CN113642865A (zh) * | 2021-07-30 | 2021-11-12 | 深圳市中金岭南有色金属股份有限公司凡口铅锌矿 | 一种矿山风险等级的确定方法、装置及服务器 |
CN113759097B (zh) * | 2021-09-07 | 2023-06-16 | 重庆大学 | 基于煤矿巷道围岩应力在线监测系统的应力状态分析方法 |
CN113515802B (zh) * | 2021-09-14 | 2021-12-07 | 四川交达预应力工程检测科技有限公司 | 基于机器学习的锚固临界值检测方法、系统及存储介质 |
CN114233386B (zh) * | 2021-12-09 | 2024-01-16 | 北京安科兴业矿山安全技术研究院有限公司 | 一种基于多参量风险判识数据库的煤矿灾害预警方法 |
CN114687737B (zh) * | 2022-04-15 | 2023-08-11 | 中国中煤能源股份有限公司 | 冲击地压巷道的支护方法 |
CN115452216B (zh) * | 2022-09-23 | 2023-04-07 | 安徽理工大学 | 一种基于互联网的煤矿巷道光纤智能在线监测系统 |
CN116085053B (zh) * | 2023-04-12 | 2023-08-01 | 平顶山天安煤业股份有限公司 | 一种基于传感技术的工作面矿压在线监测系统 |
CN116591777B (zh) * | 2023-05-16 | 2023-10-24 | 中国矿业大学 | 多场多源信息融合的冲击地压智能化监测预警装置及方法 |
CN117128044B (zh) * | 2023-08-28 | 2024-04-02 | 浙江华东测绘与工程安全技术有限公司 | 一种地下洞室施工围岩稳定安全在线预警方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1463915A1 (ru) * | 1987-07-08 | 1989-03-07 | Научно-производственное объединение "Сибцветметавтоматика" | Устройство дл контрол степени удароопасности горных пород |
SU1712601A1 (ru) * | 1989-05-12 | 1992-02-15 | Научно-производственное объединение "Сибцветметавтоматика" | Устройство дл определени степени удароопасности горных пород |
CN102628373B (zh) * | 2012-04-27 | 2014-04-23 | 天地科技股份有限公司 | 煤矿冲击地压分源综合预警方法 |
CN102644482B (zh) * | 2012-05-18 | 2014-04-02 | 河南大有能源股份有限公司 | 冲击地压预测预警方法 |
CN103306722B (zh) * | 2013-06-21 | 2015-08-05 | 中国矿业大学 | 一种冲击危险区的微震多维信息综合区域探测评价方法 |
CN103362551A (zh) * | 2013-07-23 | 2013-10-23 | 中国矿业大学 | 一种冲击矿压的综合指数评估方法 |
CN104153814A (zh) * | 2014-06-23 | 2014-11-19 | 中国矿业大学 | 一种冲击地压的微震多参量预警方法 |
-
2015
- 2015-11-16 CN CN201510786271.9A patent/CN105257339B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105257339A (zh) | 2016-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105257339B (zh) | 掘进工作面多参量综合监测预警方法 | |
CN103233777B (zh) | 一种顶板安全矿压监测方法及装置 | |
CN111058855B (zh) | 一种盾构下穿的结构物的形变控制方法及评估系统 | |
CN105785471A (zh) | 一种矿井预开采煤层的冲击危险性评价方法 | |
CN106680894B (zh) | 一种基于c/s架构的隧洞地质超前预报方法 | |
Zhang et al. | Criterion and classification for strain mode rockbursts based on five-factor comprehensive method | |
CN107387166A (zh) | 回采工作面煤层底板破坏深度实时监测预警系统及方法 | |
CN103308946B (zh) | 一种基于炮孔钻进信息的隧道超前地质预报方法 | |
CN102913285B (zh) | 矿井顶板危险报警的方法 | |
CN105759010B (zh) | 一种采动影响巷道动态监测与稳定评价方法 | |
CN109268072A (zh) | 煤矿底板突水灾害智能、实时预测预警的大数据云平台 | |
CN110043316A (zh) | 一种矿山复合动力灾害综合监测预警方法 | |
CN107764323B (zh) | 一种多参数封闭区危险性实时监测预警方法 | |
CN102979071A (zh) | 深基坑受力稳定性远程智能监测及三维预警方法与系统 | |
Coşar | Application of rock mass classification systems for future support design of the Dim Tunnel near Alanya | |
CN106401652A (zh) | 一种隧道远程在线遥测系统串联式组网方法 | |
CN113375753A (zh) | 一种煤矿综采工作面采动对地下水影响的监测与分析方法 | |
CN201857953U (zh) | 采煤工作面超前支承压力监测装置 | |
Kim et al. | Evaluation of monitoring items for adverse ground conditions in subsea tunneling | |
Zhang et al. | Stability of a roadway below a coal seam under dynamic pressure: a case study of the 11123 floor gas drainage roadway of a mine in huainan, China | |
Xiao et al. | Effect of open-pit blasting vibrations on a Hanging-Wall slope: A Case Study of the Beizhan Iron Mine in China | |
CN105698853B (zh) | 一种工作面回采期间围岩锚注支护效果的检测方法 | |
Cheng et al. | Vibration behavior during underground drilling based on an innovative measurement method and the application | |
CN113958367A (zh) | 一种地下矿山软岩巷道地压安全预警监测方法 | |
Quick et al. | Aicha‐Mauls on the Brenner Base Tunnel–status of the works and results |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |