CN112377257B - 一种基于微震监测的工作面采动超前影响范围确定方法 - Google Patents
一种基于微震监测的工作面采动超前影响范围确定方法 Download PDFInfo
- Publication number
- CN112377257B CN112377257B CN202011154941.2A CN202011154941A CN112377257B CN 112377257 B CN112377257 B CN 112377257B CN 202011154941 A CN202011154941 A CN 202011154941A CN 112377257 B CN112377257 B CN 112377257B
- Authority
- CN
- China
- Prior art keywords
- earthquake
- mining
- advanced
- mine
- total
- 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
- 238000005065 mining Methods 0.000 title claims abstract description 61
- 238000012544 monitoring process Methods 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 23
- 239000003245 coal Substances 0.000 claims abstract description 21
- 238000000605 extraction Methods 0.000 claims abstract description 16
- 238000005192 partition Methods 0.000 claims description 20
- 238000010606 normalization Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000003672 processing method Methods 0.000 claims description 3
- 239000007787 solid Substances 0.000 abstract description 3
- 239000011435 rock Substances 0.000 description 10
- 229910052500 inorganic mineral Inorganic materials 0.000 description 4
- 239000011707 mineral Substances 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000001808 coupling effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK 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
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Geochemistry & Mineralogy (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于微震监测的工作面采动超前影响范围确定方法,利用矿井微震监测系统采集工作面回采过程中的自然矿震信号,确定震源位置及能量;按一定数据步距选择若干回采阶段,统计各阶段矿震分布范围确定矿震统计区间;将矿震统计区间内实体煤区域按一定间距进行等分形成多个超前分区,统计各阶段各个超前分区内矿震总能量和总频次并进行叠加;以超前分区中线距回采阶段中线距离为横坐标、相应超前分区矿震总能量和总频次为纵坐标,绘制矿震总能量和总频次变化曲线;将两曲线进行归一化处理,并在两条归一化曲线上确定工作面矿震能量和矿震频次超前影响范围,在两者范围内选择横坐标最大值为工作面采动超前影响范围。
Description
技术领域
本发明涉及一种基于微震监测的工作面采动超前影响范围确定方法,属于煤矿安全技术领域。
背景技术
我国煤炭资源赋存条件复杂,近年来随开采深度的增加,采场应力环境不断恶化,冲击地压、煤与瓦斯突出等动力灾害频繁发生,严重威胁了矿井安全生产及人员生命安全。现场监测和统计结果表明,工作面回采过程中超前两巷一定范围内的巷道会发生不同程度的顶板下沉、鼓帮及底鼓现象,冲击地压、煤与瓦斯突出等动力事件亦多发生于采煤工作面及其超前两巷一定范围内。工作面超前采动影响范围的确定对工作面尤其是深部具有冲击地压和煤与瓦斯突出危险的工作面采掘布局、卸压解危措施的实施、超前支护和限员等至关重要。目前,采煤工作面超前影响范围确定方法主要有理论分析、数值计算和现场监测(锚杆索受力监测、应力计监测)等。理论分析和数值计算方法往往对现场条件进行了大量简化处理,与实际复杂赋存条件和力学性质差异较大;现场锚杆(索)受力监测、应力计监测等方法受装备自身及安装环境影响较大,往往达不到理想效果。工作面矿震、冲击地压及煤与瓦斯突出等显现是煤岩体静载和工作面采动动载耦合作用的结果,而传统理论分析、数值计算及现场应力监测方法只是对煤岩体静载应力的监测,未能实现对工作面动载的有效监测。
目前,较多矿井尤其具有冲击地压等灾害危险矿井均安装了微震监测系统,能够实现矿井微震事件的有效监测,微震监测技术已是冲击地压等动力灾害研究的一个重要手段。但是目前还没有利用微震监测系统进行采煤工作面采动超前影响范围的确定方法。
发明内容
针对上述现有技术存在的问题,本发明提供一种基于微震监测的工作面采动超前影响范围确定方法,无需额外增加监测设备,能有效监测煤岩体静载和工作面采动动载对工作面耦合影响,最终能精确的确定采煤工作面采动超前影响范围。
为了实现上述目的,本发明采用的技术方案是:一种基于微震监测的工作面采动超前影响范围确定方法,具体步骤为:
(1)利用矿井微震监测系统采集工作面回采过程中的多个自然矿震信号,并将各个自然矿震信号采用已知处理方法确定各自的矿震震源位置及其能量;
(2)先设定数据步距为c(通过设定数据步距能在整个工作面回采长度内均匀选择若干个回采阶段,保证后续统计结果能够覆盖整个工作面回采长度),然后根据该数据步距选择多个工作面回采距离为a的回采阶段,利用步骤(1)确定的数据统计每个回采阶段中各自的矿震超前分布范围,取其中的矿震超前最大值为矿震统计区间L;
(3)将矿震统计区间L中的实体煤区域(即矿震统计区间去除一个回采区域,即L-a)按照间距为b进行等分,获得n个超前分区;分别统计步骤(2)选择的每个回采阶段落入n个超前分区内的矿震总能量和矿震总频次,并将各个回采阶段在每个超前分区内的矿震总能量和矿震总频次进行叠加,具体计算公式为:
式中,En为各个超前分区总能量,J;Qn为各个超前分区总频次,个;ein为各个回采阶段在每个超前分区叠加总能量,J;qin为各个回采阶段在每个超前分区叠加总频次;m为回采阶段的数量;
(4)计算各个超前分区中线分别距回采阶段中线的距离,以各个超前分区中线距回采阶段中线的距离为横坐标,以各个超前分区内矿震总能量和总频次为纵坐标,分别绘制出超前工作面矿震总能量和总频次变化曲线,其中距离的计算公式为:
d0=0即回采阶段中线本身为零。
式中,di为各个超前分区中线距回采阶段中线的距离,m;a为回采距离,m;b为超前分区的间距,m;
(5)将步骤(4)绘制出的两条曲线数据除以各自纵坐标最大值进行归一化处理,分别得出矿震总能量归一化分布曲线和矿震总频次归一化分布曲线,具体公式为:
Qmax=max(Qi),i=0,1,2,…,n
Emax=max(Ei),i=0,1,2,…,n
式中,Qmax为各个超前分区中的总频次最大值,J;Emax为各个超前分区中的总能量最大值,个;R(Qi)为各个超前分区总频次除以总频次最大值的相对值;R(Ei)为各个超前分区总能量除以总能量最大值的相对值;
取两条归一化分布曲线中相对值大于等于10%的数据确定为工作面矿震能量和矿震频次超前影响范围,在工作面矿震能量和矿震频次超前影响范围内,选择两条归一化分布曲线中横坐标最大值确定为工作面采动超前影响范围。
进一步,所述步骤(2)选择的回采阶段数量不少于3个,回采距离a为5~20m,超前分区的间距b为1~10m。
进一步,所述步骤(1)采集的自然矿震信号为通过已知方法将干扰波剔除后的矿震信号。
与现有技术相比,本发明利用井下已有的矿井微震监测系统,先采集工作面回采过程中的自然矿震信号,确定震源位置及能量;按一定数据步距选择若干回采阶段,统计各阶段矿震分布范围,取范围最大值为矿震统计区间;将矿震统计区间内实体煤区域按一定间距进行等分形成多个超前分区,统计各阶段各个超前分区内矿震总能量和总频次并进行叠加;以超前分区中线距回采阶段中线距离为横坐标、相应超前分区矿震总能量和总频次为纵坐标,绘制矿震总能量和总频次变化曲线;将两曲线数据除以各自纵坐标最大值进行归一化处理,两条归一化曲线分别取相对值大于等于10%的部分确定为工作面矿震能量和矿震频次超前影响范围,在两者范围内选择横坐标最大值为工作面采动超前影响范围。本发明无需额外增加监测设备,能有效监测煤岩体静载和工作面采动动载对工作面耦合影响,最终能精确的确定采煤工作面采动超前影响范围。
附图说明
图1是本发明确定工作面采动超前影响范围的计算过程示意图;
图2是实施例1绘制出的超前工作面矿震总能量分布曲线;
图3是实施例1绘制出的超前工作面矿震总频次分布曲线;
图4是实施例1得出的矿震总能量归一化分布曲线
图5是实施例1得出的矿震总频次归一化分布曲线。
具体实施方式
下面将对本发明作进一步说明。
实施例1:如图1所示,本实施例针彬长矿区某矿204工作面具有强冲击危险,工作面已回采1050m,发生过冲击显现,矿井已安装微震监测系统,采用本发明确定其采动超前影响范围,具体步骤为:
(1)利用该矿SOS微震监测系统监测204工作面回采过程中工作面及临近区域8000余次2次方以上矿震信号,并将各个自然矿震信号剔除干扰信号后采用已知处理方法确定各自的矿震震源位置及其能量;
(2)先设定数据步距为50m,在204工作面选择20个回采距离为10m的回采阶段,利用步骤(1)确定的数据统计每个回采阶段中各自的矿震超前分布范围,最大值为超前500m,确定本实施例中矿震统计区间L=500m;
(3)将矿震统计区间L中的实体煤区域(即矿震统计区间去除一个回采区域,即L-a)为490m,按照间距为5m进行等分,获得98个超前分区;分别统计步骤(2)选择的每个回采区域落入98个超前分区内的矿震总能量和矿震总频次,并将各个回采区域在每个超前分区内的矿震总能量和矿震总频次进行叠加,具体计算公式为:
式中,En为各个超前分区总能量,J;Qn为各个超前分区总频次,个;ein为各个回采区域在每个超前分区叠加总能量,J;qin为各个回采区域在每个超前分区叠加总频次;m为回采区域的数量;
(4)计算98个超前分区分别距10m回采阶段中线的距离,分别为7.5m、12.5m、···、497.5m,回采区域距离为0m,以此为横坐标,以(3)中相应分区内矿震总能量和总频次为纵坐标,分别绘制超前区域矿震总能量和总频次变化曲线,如图2和图3所示,其中距离的计算公式为:
d0=0即回采阶段中线本身为零。
式中,di为各个超前分区中线距回采区域中线的距离,m;a为回采距离,m;b为超前分区的间距,m;
(5)确定两曲线纵坐标最大值,分别为2.3×106J和44个,将两曲线各个数据分别除以各自的纵坐标最大值进行归一化处理,分别得出矿震总能量归一化分布曲线和矿震总频次归一化分布曲线,如图4和图5所示,
通过图4和图5,取两条归一化分布曲线中相对值大于等于10%的数据确定为工作面矿震能量和矿震频次超前影响范围(即横坐标最大值),分别为112.5m和132.5m,选择两者中的最大值132.5m确定为工作面采动超前影响范围。
Claims (3)
1.一种基于微震监测的工作面采动超前影响范围确定方法,其特征在于,具体步骤为:
(1)利用矿井微震监测系统采集工作面回采过程中的多个自然矿震信号,并将各个自然矿震信号采用已知处理方法确定各自的矿震震源位置及其能量;
(2)先设定数据步距为c,然后根据该数据步距选择多个工作面回采距离为a的回采阶段,利用步骤(1)确定的数据统计每个回采阶段中各自的矿震超前分布范围,取其中的矿震超前最大值为矿震统计区间L;
(3)将矿震统计区间L中的实体煤区域按照间距为b进行等分,获得n个超前分区;分别统计步骤(2)选择的每个回采阶段落入n个超前分区内的矿震总能量和矿震总频次,并将各个回采阶段在每个超前分区内的矿震总能量和矿震总频次进行叠加,具体计算公式为:
式中,En为各个超前分区总能量,J;Qn为各个超前分区总频次,个;ein为各个回采阶段在每个超前分区叠加总能量,J;qin为各个回采阶段在每个超前分区叠加总频次;m为回采阶段的数量;
(4)计算各个超前分区中线分别距回采阶段中线的距离,以各个超前分区中线距回采阶段中线的距离为横坐标,以各个超前分区内矿震总能量和总频次为纵坐标,分别绘制出超前工作面矿震总能量和总频次变化曲线,其中距离的计算公式为:
d0=0
式中,di为各个超前分区中线距回采阶段中线的距离,m;a为回采距离,m;b为超前分区的间距,m;
(5)将步骤(4)绘制出的两条曲线数据除以各自纵坐标最大值进行归一化处理,分别得出矿震总能量归一化分布曲线和矿震总频次归一化分布曲线,具体公式为:
Qmax=max(Qi),i=0,1,2,…,n
Emax=max(Ei),i=0,1,2,…,n
式中,Qmax为各个超前分区中的总频次最大值,J;Emax为各个超前分区中的总能量最大值,个;R(Qi)为各个超前分区总频次除以总频次最大值的相对值;R(Ei)为各个超前分区总能量除以总能量最大值的相对值;
取两条归一化分布曲线中相对值大于等于10%的数据确定为工作面矿震能量和矿震频次超前影响范围,在工作面矿震能量和矿震频次超前影响范围内,选择两条归一化分布曲线中横坐标最大值确定为工作面采动超前影响范围。
2.根据权利要求1所述的一种基于微震监测的工作面采动超前影响范围确定方法,其特征在于,所述步骤(2)选择的回采阶段数量不少于3个,回采距离a为5~20m,超前分区的间距b为1~10m。
3.根据权利要求1所述的一种基于微震监测的工作面采动超前影响范围确定方法,其特征在于,所述步骤(1)采集的自然矿震信号为通过已知方法将干扰波剔除后的矿震信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011154941.2A CN112377257B (zh) | 2020-10-26 | 2020-10-26 | 一种基于微震监测的工作面采动超前影响范围确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011154941.2A CN112377257B (zh) | 2020-10-26 | 2020-10-26 | 一种基于微震监测的工作面采动超前影响范围确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112377257A CN112377257A (zh) | 2021-02-19 |
CN112377257B true CN112377257B (zh) | 2021-08-20 |
Family
ID=74576212
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011154941.2A Active CN112377257B (zh) | 2020-10-26 | 2020-10-26 | 一种基于微震监测的工作面采动超前影响范围确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112377257B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113279813B (zh) * | 2021-05-31 | 2022-02-18 | 北京科技大学 | 一种基于动载的巷道支护方法 |
CN113434819B (zh) * | 2021-06-24 | 2023-08-01 | 中国矿业大学 | 工作面采动对采空区矿震活动影响时间和距离的确定方法 |
CN113962466B (zh) * | 2021-10-27 | 2024-04-30 | 陕西正通煤业有限责任公司 | 基于微震的冲击地压矿井褶曲区域冲击危险动态预测方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103294926B (zh) * | 2013-06-21 | 2016-03-23 | 中国矿业大学 | 一种冲击矿压时序监测预警方法 |
CN103291364B (zh) * | 2013-06-21 | 2015-08-05 | 中国矿业大学 | 一种冲击矿压的微震多维信息综合时序预警方法 |
CN105673075A (zh) * | 2016-01-13 | 2016-06-15 | 中国矿业大学(北京) | 一种煤岩动力灾害多参量无线监测综合预警技术与方法 |
CN106154307B (zh) * | 2016-06-21 | 2018-11-20 | 中国矿业大学 | 一种煤岩冲击失稳模式的微震识别方法 |
CN110118991B (zh) * | 2019-05-16 | 2020-06-23 | 中国矿业大学 | 一种基于微震损伤重构的采动应力评估方法 |
-
2020
- 2020-10-26 CN CN202011154941.2A patent/CN112377257B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112377257A (zh) | 2021-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112377257B (zh) | 一种基于微震监测的工作面采动超前影响范围确定方法 | |
CN109611143B (zh) | 顶板水压致裂多参量综合监测系统及致裂效果判别方法 | |
CN209761499U (zh) | 顶板水压致裂多参量综合监测系统 | |
Zhang et al. | Microseismicity induced by fault activation during the fracture process of a crown pillar | |
Zhang et al. | Assessment of rockburst risk in deep mining: an improved comprehensive index method | |
CN104653226A (zh) | 一种基于应力梯度的煤矿冲击地压危险区域的划分方法 | |
CN112377258B (zh) | 一种临近采区同采工作面间采动影响程度的评价方法 | |
Liu et al. | A case study of collapses at the Yangshan tunnel of the Coal Transportation Channel from the Western Inner Mongolia to the Central China | |
CN113279813A (zh) | 一种基于动载的巷道支护方法 | |
CN103606019A (zh) | 基于时空关系的矿山采空区覆岩沉降动态预测方法 | |
Yuan-hui et al. | The spatial-temporal evolution law of microseismic activities in the failure process of deep rock masses | |
CN110761841B (zh) | 一种基于微震监测的矿井群工作面开采互扰范围计算方法 | |
CN115169948A (zh) | 采煤工作面覆岩离层突水风险预测方法及安全采矿方法 | |
CN115220092A (zh) | 一种确定工作面超前冲击危险范围的微震统计方法 | |
CN109555502B (zh) | 一种高瓦斯煤巷预裂增透的工业性试验方法 | |
Chen et al. | Reliability assessment of the hydraulic fracturing process in coal mine based on the analysis of micro-seismic source parameters | |
Wang et al. | Microseismicity evolution related to two extremely intense rockbursts in a water diversion tunnel | |
CN113250613B (zh) | 小煤窑采空区顺煤层定向钻探排查方法 | |
CN112377256B (zh) | 一种基于微震监测的冲击危险性大直径钻孔检验方法 | |
CN113434819B (zh) | 工作面采动对采空区矿震活动影响时间和距离的确定方法 | |
Wang et al. | Evolution Laws of Floor Stress and Stability of Floor Roadway Affected by Overhead Mining | |
CN112630826B (zh) | 一种基于微震监测的工作面采动倾向影响范围的确定方法 | |
Li et al. | Movement law and discriminant method of key strata breakage based on microseismic monitoring | |
CN109598029A (zh) | 一种应用于高瓦斯煤巷co2高压气体致裂的设计方法 | |
CN115128670A (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 |