CN112213767A - 一种用于顶板超前预裂***效果的评价方法 - Google Patents

一种用于顶板超前预裂***效果的评价方法 Download PDF

Info

Publication number
CN112213767A
CN112213767A CN202010952999.5A CN202010952999A CN112213767A CN 112213767 A CN112213767 A CN 112213767A CN 202010952999 A CN202010952999 A CN 202010952999A CN 112213767 A CN112213767 A CN 112213767A
Authority
CN
China
Prior art keywords
blasting
mine
spontaneous
wave velocity
inversion
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.)
Pending
Application number
CN202010952999.5A
Other languages
English (en)
Inventor
阚吉亮
王绪俊
窦林名
贾晓东
巩思园
程文武
刘强
马志锋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xuzhou Hongyi Science And Technology Development Co ltd
China University of Mining and Technology CUMT
Yanzhou Coal Mining Co Ltd
Original Assignee
Xuzhou Hongyi Science And Technology Development Co ltd
China University of Mining and Technology CUMT
Yanzhou Coal Mining Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xuzhou Hongyi Science And Technology Development Co ltd, China University of Mining and Technology CUMT, Yanzhou Coal Mining Co Ltd filed Critical Xuzhou Hongyi Science And Technology Development Co ltd
Priority to CN202010952999.5A priority Critical patent/CN112213767A/zh
Publication of CN112213767A publication Critical patent/CN112213767A/zh
Pending legal-status Critical Current

Links

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
    • 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/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/10—Aspects of acoustic signal generation or detection
    • G01V2210/12—Signal generation
    • G01V2210/121—Active source
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/10—Aspects of acoustic signal generation or detection
    • G01V2210/14—Signal detection
    • G01V2210/142—Receiver location
    • G01V2210/1429—Subsurface, e.g. in borehole or below weathering layer or mud line
    • 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/64—Geostructures, e.g. in 3D data cubes

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Acoustics & Sound (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)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种用于顶板超前预裂***效果的评价方法，先确定波速反演区域，在实施顶板超前预裂***措施前，矿井微震监测***实时接收产生的多次自发矿震，采用已知的波速反演方法对自发矿震产生的震动波进行被动CT反演，获得顶板岩体在***前的波速分布；接着在设定位置实施顶板超前预裂***，完成后重复***前的波速反演过程，获得顶板岩体在***后的波速分布；根据***前后被动CT反演结果，对顶板不同区域在***前后的波速进行对比分析，计算获得***区域的波速变化系数，对各个区域进行梯度分类，获得每个梯度范围内的区域面积；最后结合区域面积和波速变化系数，通过公式综合评价***致裂程度和致裂范围，保证致裂效果评价的精确性。

Description

一种用于顶板超前预裂***效果的评价方法
技术领域
本发明涉及一种用于顶板超前预裂***效果的评价方法，属于煤矿开采和煤矿安全技术领域。
背景技术
近年来，随着开采深度和开采强度的增大，煤矿冲击矿压问题越来越严重，其中坚硬顶板型冲击矿压主要是由于采空区上方的坚硬岩层无法垮落而形成大范围悬顶。一般情况下，为促进顶板能及时垮落，可采用超前预裂***破坏顶板岩层的完整性和连续性，使顶板能够及时垮落，降低冲击矿压发生的可能性。现有顶板预裂***效果的评价过程主要通过钻孔窥视仪观察***区域岩层裂隙内部发育情况，进而进行评价；但是钻孔窥视需要钻设大量的探测钻孔，从而导致施工工程量和成本较大且探测精度较低。另外，目前大部分的钻孔窥视钻孔长度一般在50m范围内，因此不适用于深孔预裂***(钻孔长度大于50m)致裂效果的探测。
发明内容
针对上述现有技术存在的问题，本发明提供一种用于顶板超前预裂***效果的评价方法，通过自发矿震进行波速反演评价顶板超前预裂***效果，不仅无需额外钻设探测钻孔，降低施工成本及施工量，而且能有效对深孔预裂***的致裂效果进行探测，保证致裂效果评价的精确性。
为了实现上述目的，本发明采用的技术方案是：一种用于顶板超前预裂***效果的评价方法，具体步骤为：
(1)根据实施顶板超前预裂***的设定位置，以预裂***孔为中心、边长为100m的立方体范围确定为波速反演区域；
(2)在工作面两侧巷道围绕反演区域交错均匀布置数量不少于8个拾震器，同一巷道内相邻布置的拾震器水平间距在150～200m范围内，进而组成矿井微震监测***；这种布置方式使得拾震器尽可能的有利于对反演区域形成立体全方位的监测；
(3)在实施顶板超前预裂***措施前，矿井微震监测***实时接收工作面回采过程产生的多次自发矿震，每次自发矿震均会产生震动波，然后对接收到的震动波进行筛选及被动CT反演可行性条件判别，直至获取符合要求的自发矿震作为震源，最后采用已知的波速反演方法对符合要求的自发矿震产生的震动波进行被动CT反演，获得顶板岩体在***前的波速分布；
(4)在设定位置实施顶板超前预裂***，完成后矿井微震监测***实时接收工作面回采过程产生的多次自发矿震，每次自发矿震均会产生震动波，然后对接收到的震动波进行筛选及被动CT反演可行性条件判别，直至获取符合要求的自发矿震作为震源，最后采用已知的波速反演方法对符合要求的自发矿震产生的震动波进行被动CT反演，获得顶板岩体在***后的波速分布；
(5)根据***前后被动CT反演结果，对顶板不同区域在***前后的波速进行对比分析，计算获得***区域的波速变化系数A，评估***预裂的致裂程度；
Figure BDA0002677656100000021
式中：A为波速变化系数；V1为***前的波速值；V2为***后的波速值；
(6)根据顶板岩体不同区域波速变化系数A，以A＝5％为梯度值进行梯度分类，即A＝0～5％为第一波速变化梯度；A＝5％～10％为第二波速变化梯度；依次类推，完成所有区域的梯度分类，根据分类结果获得每个梯度范围内的区域面积S；
(7)计算***预裂效果评价指标W:
Figure BDA0002677656100000022
其中：i＝1～n表示不同波速变化梯度，
Figure BDA0002677656100000023
表示当前梯度下的平均波速变化系数，Si表示当前梯度下的区域面积；SA表示所有波速变化梯度下的区域总面积；
最终计算获得W值，根据W值综合评价***致裂程度和致裂范围；并且W值越大，表明***预裂措施的致裂效果越好。
进一步，所述对接收到的震动波进行筛选及被动CT反演可行性条件判别，具体步骤为：
①自发矿震的初筛选
在工作面回采过程中会产生多次自发矿震，多个拾震器实时接收数据，提取自发矿震产生震动波的波形信号特征，筛选出波形信号接收通道数量N大于6的自发矿震，基于波形频谱信息计算筛选出的各个自发矿震产生震动波的波形信噪比(S/N)，筛选出波形信噪比(S/N)大于14dB的自发矿震；
②自发矿震的再筛选
根据矿井微震监测***中拾震器的空间位置，对步骤①初筛选出的自发矿震进行震源空间定位，进而绘制各自发矿震及其接收拾震器的传播路径；其中传播路径穿过波速反演区域的自发矿震确定为能用于震动波被动CT反演的有效矿震；
③震动波被动CT反演可行性条件判别
将所有的有效矿震与各个拾震器之间的传播路径形成射线网，射线网的疏密程度决定震动波被动CT反演的结果准确性与效率，设定可行性条件为：筛选后的有效矿震数量要大于100个；波速反演区域内进行网格化划分后，单位网格内的射线数量要大于60条，若满足上述两个条件，则确定当前有效矿震能够达到波速反演的条件。
与现有技术相比，本发明先根据顶板超前预裂***的设定位置，确定波速反演区域，通过矿井微震监测***在实施顶板超前预裂***措施前，实时接收工作面回采过程产生的多次自发矿震，然后对接收到的震动波进行筛选及被动CT反演可行性条件判别，直至获取符合要求的自发矿震作为震源，最后采用已知的波速反演方法对符合要求的自发矿震产生的震动波进行被动CT反演，获得顶板岩体在***前的波速分布；接着在设定位置实施顶板超前预裂***，完成后重复***前的波速反演过程，获得顶板岩体在***后的波速分布；根据***前后被动CT反演结果，对顶板不同区域在***前后的波速进行对比分析，计算获得***区域的波速变化系数，然后设定梯度值，对各个区域进行梯度分类，获得每个梯度范围内的区域面积；最后结合区域面积和波速变化系数，通过公式计算得出***预裂效果评价指标W；根据W值综合评价***致裂程度和致裂范围；并且W值越大，表明***预裂措施的致裂效果越好。因此本发明能利用自发矿震进行波速反演评价顶板超前预裂***效果，不仅无需额外钻设探测钻孔，降低施工成本及施工量，而且能有效对深孔预裂***的致裂效果进行探测，保证致裂效果评价的精确性，为后续的开采工作提供理论指导。
附图说明
图1是本发明中波速反演区域位置选择示意图；
图2是本发明中自发矿震产生震动波的波形信息图；
图3是本发明中震动波的矿震频谱特征图；
图4是本发明中自发矿震传播路径以及有效反演区域示意图；
图5是本发明中***前后震动波被动CT反演结果对比分析示意图。
具体实施方式
下面将对本发明作进一步说明。
如图1至图5所示，本发明的具体步骤为：
(1)根据实施顶板超前预裂***的设定位置，以预裂***孔为中心、边长为100m的立方体范围确定为波速反演区域；
(2)在工作面两侧巷道围绕反演区域交错均匀布置数量不少于8个拾震器，同一巷道内相邻布置的拾震器水平间距在150～200m范围内，进而组成矿井微震监测***；这种布置方式使得拾震器尽可能的有利于对反演区域形成立体全方位的监测；
(3)在实施顶板超前预裂***措施前，矿井微震监测***实时接收工作面回采过程产生的多次自发矿震，每次自发矿震均会产生震动波，然后对接收到的震动波进行筛选及被动CT反演可行性条件判别，直至获取符合要求的自发矿震作为震源，最后采用已知的波速反演方法对符合要求的自发矿震产生的震动波进行被动CT反演，获得顶板岩体在***前的波速分布；所述对接收到的震动波进行筛选及被动CT反演可行性条件判别，具体步骤为：
①自发矿震的初筛选
在工作面回采过程中会产生多次自发矿震，多个拾震器实时接收数据，提取自发矿震产生震动波的波形信号特征，筛选出波形信号接收通道数量N大于6的自发矿震，基于波形频谱信息计算筛选出的各个自发矿震产生震动波的波形信噪比(S/N)，筛选出波形信噪比(S/N)大于14dB的自发矿震；
②自发矿震的再筛选
根据矿井微震监测***中拾震器的空间位置，对步骤①初筛选出的自发矿震进行震源空间定位，进而绘制各自发矿震及其接收拾震器的传播路径；其中传播路径穿过波速反演区域的自发矿震确定为能用于震动波被动CT反演的有效矿震；
③震动波被动CT反演可行性条件判别
将所有的有效矿震与各个拾震器之间的传播路径形成射线网，射线网的疏密程度决定震动波被动CT反演的结果准确性与效率，设定可行性条件为：筛选后的有效矿震数量要大于100个；波速反演区域内进行网格化划分后，单位网格内的射线数量要大于60条，若满足上述两个条件，则确定当前有效矿震能够达到波速反演的条件。
表1：震动波筛选及被动CT反演可行性条件判别：
Figure BDA0002677656100000051
(4)在设定位置实施顶板超前预裂***，完成后矿井微震监测***实时接收工作面回采过程产生的多次自发矿震，每次自发矿震均会产生震动波，然后对接收到的震动波进行筛选及被动CT反演可行性条件判别，直至获取符合要求的自发矿震作为震源，最后采用已知的波速反演方法对符合要求的自发矿震产生的震动波进行被动CT反演，获得顶板岩体在***后的波速分布；
(5)根据***前后被动CT反演结果，对顶板不同区域在***前后的波速进行对比分析，计算获得***区域的波速变化系数A，评估***预裂的致裂程度；
Figure BDA0002677656100000052
式中：A为波速变化系数；V1为***前的波速值；V2为***后的波速值；
(6)根据顶板岩体不同区域波速变化系数A，以A＝5％为梯度值进行梯度分类，即A＝0～5％为第一波速变化梯度；A＝5％～10％为第二波速变化梯度；依次类推，完成所有区域的梯度分类，根据分类结果获得每个梯度范围内的区域面积S；
(7)计算***预裂效果评价指标W:
Figure BDA0002677656100000053
其中：i＝1～n表示不同波速变化梯度，
Figure BDA0002677656100000061
表示当前梯度下的平均波速变化系数，Si表示当前梯度下的区域面积；SA表示所有波速变化梯度下的区域总面积；
最终计算获得W值，根据W值综合评价***致裂程度和致裂范围；并且W值越大，表明***预裂措施的致裂效果越好。

Claims (2)

1.一种用于顶板超前预裂***效果的评价方法，其特征在于，具体步骤为：
(1)根据实施顶板超前预裂***的设定位置，以预裂***孔为中心、边长为100m的立方体范围确定为波速反演区域；
(2)在工作面两侧巷道围绕反演区域交错均匀布置数量不少于8个拾震器，同一巷道内相邻布置的拾震器水平间距在150～200m范围内，进而组成矿井微震监测***；
(3)在实施顶板超前预裂***措施前，矿井微震监测***实时接收工作面回采过程产生的多次自发矿震，每次自发矿震均会产生震动波，然后对接收到的震动波进行筛选及被动CT反演可行性条件判别，直至获取符合要求的自发矿震作为震源，最后采用已知的波速反演方法对符合要求的自发矿震产生的震动波进行被动CT反演，获得顶板岩体在***前的波速分布；
(4)在设定位置实施顶板超前预裂***，完成后矿井微震监测***实时接收工作面回采过程产生的多次自发矿震，每次自发矿震均会产生震动波，然后对接收到的震动波进行筛选及被动CT反演可行性条件判别，直至获取符合要求的自发矿震作为震源，最后采用已知的波速反演方法对符合要求的自发矿震产生的震动波进行被动CT反演，获得顶板岩体在***后的波速分布；
(5)根据***前后被动CT反演结果，对顶板不同区域在***前后的波速进行对比分析，计算获得***区域的波速变化系数A，评估***预裂的致裂程度；
Figure FDA0002677656090000011
式中：A为波速变化系数；V1为***前的波速值；V2为***后的波速值；
(6)根据顶板岩体不同区域波速变化系数A，以A＝5％为梯度值进行梯度分类，即A＝0～5％为第一波速变化梯度；A＝5％～10％为第二波速变化梯度；依次类推，完成所有区域的梯度分类，根据分类结果获得每个梯度范围内的区域面积S；
(7)计算***预裂效果评价指标W:
Figure FDA0002677656090000021
其中：i＝1～n表示不同波速变化梯度，
Figure FDA0002677656090000022
表示当前梯度下的平均波速变化系数，Si表示当前梯度下的区域面积，SA表示所有波速变化梯度下的区域总面积；
最终计算获得W值，根据W值综合评价***致裂程度和致裂范围；并且W值越大，表明***预裂措施的致裂效果越好。
2.根据权利要求1所述的一种用于顶板超前预裂***效果的评价方法，其特征在于，所述对接收到的震动波进行筛选及被动CT反演可行性条件判别，具体步骤为：
①自发矿震的初筛选
在工作面回采过程中会产生多次自发矿震，多个拾震器实时接收数据，提取自发矿震产生震动波的波形信号特征，筛选出波形信号接收通道数量N大于6的自发矿震，基于波形频谱信息计算筛选出的各个自发矿震产生震动波的波形信噪比，筛选出波形信噪比大于14dB的自发矿震；
②自发矿震的再筛选
根据矿井微震监测***中拾震器的空间位置，对步骤①初筛选出的自发矿震进行震源空间定位，进而绘制各自发矿震及其接收拾震器的传播路径；其中传播路径穿过波速反演区域的自发矿震确定为能用于震动波被动CT反演的有效矿震；
③震动波被动CT反演可行性条件判别
将所有的有效矿震与各个拾震器之间的传播路径形成射线网，射线网的疏密程度决定震动波被动CT反演的结果准确性与效率，设定可行性条件为：筛选后的有效矿震数量要大于100个；波速反演区域内进行网格化划分后，单位网格内的射线数量要大于60条，若满足上述两个条件，则确定当前有效矿震能够达到波速反演的条件。
CN202010952999.5A 2020-09-11 2020-09-11 一种用于顶板超前预裂***效果的评价方法 Pending CN112213767A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010952999.5A CN112213767A (zh) 2020-09-11 2020-09-11 一种用于顶板超前预裂***效果的评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010952999.5A CN112213767A (zh) 2020-09-11 2020-09-11 一种用于顶板超前预裂***效果的评价方法

Publications (1)

Publication Number Publication Date
CN112213767A true CN112213767A (zh) 2021-01-12

Family

ID=74049968

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010952999.5A Pending CN112213767A (zh) 2020-09-11 2020-09-11 一种用于顶板超前预裂***效果的评价方法

Country Status (1)

Country Link
CN (1) CN112213767A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113339072A (zh) * 2021-07-06 2021-09-03 中国矿业大学 一种基于微震信号波形分析的***卸压效果评价方法
CN114924311A (zh) * 2022-05-17 2022-08-19 中国矿业大学 一种基于顶板***诱发震动能量的释能效果定量评估方法
CN114994791A (zh) * 2022-05-27 2022-09-02 中国矿业大学 一种用于评估井地一体微震监测***监测能力的方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105224742A (zh) * 2015-09-29 2016-01-06 鞍钢集团矿业公司 一种分析***作用下台阶边坡稳定性的方法
CN106291702A (zh) * 2016-10-17 2017-01-04 中国矿业大学 一种矿井采区应力集中区的时移地震监测方法
CN107169686A (zh) * 2017-07-06 2017-09-15 中国神华能源股份有限公司 一种对煤层顶板下方的煤柱进行***时的***效果评价方法
CN108594295A (zh) * 2018-05-31 2018-09-28 天地科技股份有限公司 一种煤层***卸压效果评价方法
CN109596020A (zh) * 2018-11-15 2019-04-09 中国矿业大学(北京) 一种定向断裂光面***的分析评价方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105224742A (zh) * 2015-09-29 2016-01-06 鞍钢集团矿业公司 一种分析***作用下台阶边坡稳定性的方法
CN106291702A (zh) * 2016-10-17 2017-01-04 中国矿业大学 一种矿井采区应力集中区的时移地震监测方法
CN107169686A (zh) * 2017-07-06 2017-09-15 中国神华能源股份有限公司 一种对煤层顶板下方的煤柱进行***时的***效果评价方法
CN108594295A (zh) * 2018-05-31 2018-09-28 天地科技股份有限公司 一种煤层***卸压效果评价方法
CN109596020A (zh) * 2018-11-15 2019-04-09 中国矿业大学(北京) 一种定向断裂光面***的分析评价方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
万洪富 等: "《我国区域农业环境问题及其综合治理》", 31 January 2005, 中国环境科学出版社 *
巩思园: "矿震震动波波速层析成像原理及其预测煤矿冲击危险应用实践", 《中国博士学位论文全文数据库•工程科技I辑》 *
李静等: "基于震动波CT探测技术的冲击危险预测 ", 《煤矿安全》 *
窦林名 等: "断顶***防治冲击矿压技术体系及效果评价研究", 《煤炭科学技术》 *
贺可强 等: "《地质资源与环境及其承载力综合分析与评价》", 31 May 2009, 山东大学出版社 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113339072A (zh) * 2021-07-06 2021-09-03 中国矿业大学 一种基于微震信号波形分析的***卸压效果评价方法
CN113339072B (zh) * 2021-07-06 2022-07-01 中国矿业大学 一种基于微震信号波形分析的***卸压效果评价方法
CN114924311A (zh) * 2022-05-17 2022-08-19 中国矿业大学 一种基于顶板***诱发震动能量的释能效果定量评估方法
CN114994791A (zh) * 2022-05-27 2022-09-02 中国矿业大学 一种用于评估井地一体微震监测***监测能力的方法
CN114994791B (zh) * 2022-05-27 2023-03-31 中国矿业大学 一种用于评估井地一体微震监测***监测能力的方法

Similar Documents

Publication Publication Date Title
CN112213767A (zh) 一种用于顶板超前预裂***效果的评价方法
Ghosh et al. Application of underground microseismic monitoring for ground failure and secure longwall coal mining operation: a case study in an Indian mine
CN102644482B (zh) 冲击地压预测预警方法
CN109736796B (zh) 一种加深炮孔超前探测机构及其预报方法
CN112731525B (zh) 微震与电磁辐射同步监测的巷道围岩稳定性智能预测方法
RU2467171C1 (ru) Способ диагностики опасных ситуаций при подземной добыче каменного угля и методика прогноза параметров зон трещиноватости, образованной гидроразрывом пласта
CN111222254A (zh) 基于应力叠加法的工作面冲击地压危险等级划分方法和***
CN111852566A (zh) 煤矿井下水力压裂裂缝扩展规律的微震监测方法
Cao et al. Seismic effort of blasting wave transmitted in coal-rock mass associated with mining operation
CN106646640A (zh) 一种对穿越采空区隧道围岩稳定性评价方法
CN116842411A (zh) 基于震动波探测概率的煤矿微震台网监测能力评价方法
CN110632667B (zh) 一种基于冲击波激震条件下的隐伏陷落柱超前探测方法
CN115755185A (zh) 一种基于微震监测判识大能量矿震致灾性的方法
CN105223274A (zh) 利用孔内落锤探测岩体等效弹性模量的方法
CN106032750B (zh) 一种基于钻探能谱的地质编录仪
Feng et al. Microseismic energy and intensity criterion of rockburst in deep TBM tunnels: A case study of the Neelum-Jhelum hydropower project
CN113914932B (zh) 利用震动波断层扫描识别煤与瓦斯突出危险区域的方法
CN110761841B (zh) 一种基于微震监测的矿井群工作面开采互扰范围计算方法
Dong et al. Case study of microseismic tomography and multi-parameter characteristics under mining disturbances
CN113447982B (zh) 独头掘进巷道超前区域的冲击危险探测方法
FANG et al. Intelligent recognition of tunnel stratum based on advanced drilling tests
CN115263427A (zh) 覆岩破坏高度确定方法、装置、电子设备及存储介质
Vinoth et al. Slope stability monitoring by quantification and behavior of microseismic events in an opencast coal mine
RU2309434C1 (ru) Способ контроля разработки залежей углеводородов по микросейсмической эмиссии
Yu-wei et al. Intelligent recognition of tunnel stratum based on advanced drilling tests

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20210112