CN105045969B - 一种地应力型冲击地压危险性多元信息耦合预测方法 - Google Patents
一种地应力型冲击地压危险性多元信息耦合预测方法 Download PDFInfo
- Publication number
- CN105045969B CN105045969B CN201510375710.7A CN201510375710A CN105045969B CN 105045969 B CN105045969 B CN 105045969B CN 201510375710 A CN201510375710 A CN 201510375710A CN 105045969 B CN105045969 B CN 105045969B
- Authority
- CN
- China
- Prior art keywords
- stress
- gradient
- velocity
- mrow
- coal body
- 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 25
- 230000008878 coupling Effects 0.000 title claims abstract description 12
- 238000010168 coupling process Methods 0.000 title claims abstract description 12
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 12
- 239000003245 coal Substances 0.000 claims abstract description 30
- 239000011435 rock Substances 0.000 claims abstract description 11
- 238000004587 chromatography analysis Methods 0.000 claims abstract description 4
- 238000003384 imaging method Methods 0.000 claims abstract description 4
- 238000009795 derivation Methods 0.000 claims description 6
- 238000005520 cutting process Methods 0.000 claims description 4
- 238000005553 drilling Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 description 4
- 230000006378 damage Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 208000033999 Device damage Diseases 0.000 description 1
- 238000006424 Flood reaction Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000002817 coal dust Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Geophysics And Detection Of Objects (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明公开了一种地应力型冲击地压危险性多元信息耦合预测方法，(1)在冲击危险区域测试煤体中的应力值，得出应力梯度的分布情况；(2)在冲击危险区域煤体的一侧放炮或撞击煤壁，通过另一侧拾波器接收，然后进行CT层析成像技术得出波速梯度的分布情况；(3)由于应力梯度与波速梯度之间的耦合性，通过确定煤体发生冲击地压的应力梯度临界值，得出波速梯度临界值VGm；(4)在采煤机回采过程中，通过公式实时得到煤体中的波速梯度VG；(5)将实时检测得出的波速梯度VG与波速梯度临界值VGm进行比较，若VG小于VGm，则预测结果为安全；若VG大于VGm，则预测结果为危险。能实时预测冲击地压危险性，充分得出冲击地压突出危险性信息。
Description
技术领域
本发明涉及一种冲击地压危险性的预测方法，具体是一种地应力型冲击地压危险性多元信息耦合预测方法。
背景技术
近年来，随着经济发展的需求，采矿往深部发展，地质条件复杂，冲击地压发生的频率越来越高；在西北部地区，采深虽然较浅，冲击地压事故仍然多次发生。冲击地压是聚积在巷道和采场周围煤岩体中的应力能量突然释放，将煤岩抛向巷道，同时发出强烈声响，造成煤岩体震动和破坏，支架与设备损坏，人员伤亡，部分巷道垮落破坏等动力现象。冲击地压还会引发或可能引发其他矿井灾害，尤其是瓦斯、煤尘***，火灾以及水灾，干扰通风***，严重时造成地面震动和建筑物破坏等。因此，冲击地压是煤矿重大灾害之一。传统的预测方法有钻屑法、应力值法，这些方法指标单一，适用性差。现场的采掘工作面的地质条件复杂，根据冲击地压发生的机理及CT反演技术，冲击地压的发生与应力梯度有直接关系。传统的监测煤体中应力值的方法，操作复杂，只能监测局部点的应力值，无法对整个煤体中的应力值进行实时监测，而且破坏了煤体的力学结构。
发明内容
针对上述现有技术存在的问题，本发明提供一种地应力型冲击地压危险性多元信息耦合预测方法，能实时预测冲击地压危险性，充分得出冲击地压突出危险性信息。
为了实现上述目的，本发明采用的技术方案是：该种地应力型冲击地压危险性多元信息耦合预测方法，其具体步骤为：
(1)采用钻屑法或应力计法在冲击危险区域测试煤体中的应力值，然后绘制应力分布曲线，通过对距离求导，得出应力梯度的分布情况，具体公式为；
上式中：ΔS为应力点σ1和应力点σ2之间的距离，Δσ为应力梯度；
(2)在冲击危险区域煤体的一侧放炮或撞击煤壁，通过设置在另一侧拾波器接收，然后进行CT层析成像技术得到煤体中的波速分布，通过对距离求导，得出波速梯度的分布情况，具体公式为；
上式中：ΔS是测量的波速点V1和测量的波速点V2之间的距离，VG为波速梯度；
(3)由于应力梯度与波速梯度之间的耦合性，根据下述线性公式，得出随着应力梯度增加，波速梯度也线性增加，通过确定煤体发生冲击地压的应力梯度临界值，得出波速梯度临界值VGm；
Vp＝a×σλ
上式中，Vp表示P波波速；a是常数，是修正值；σ为应力值；λ是指数，常量；
(4)在采煤机回采过程中，在煤壁一侧每半小时撞击一次，在另一侧通过拾波器接收波，进而得到煤体中波速的分布，通过公式实时得到煤体中的波速梯度VG；
(5)将实时检测得出的波速梯度VG与波速梯度临界值VGm进行比较，若VG小于VGm，则预测结果为安全；若VG大于VGm，则预测结果为危险，可及时采取相应的解危措施。
与现有技术相比，本发明利用波速与应力之间的耦合特性，得出波速梯度与应力梯度之间存在线性关系，从而建立了地应力型冲击地压危险性多元信息耦合预测模型，实时预测冲击地压危险性，能充分得出冲击地压突出危险性信息，通过监测波速解决了应力值法只能监测煤体中局部应力的局限性，能实现对整个煤体进行反演分析，提高了冲击地压危险性预测的准确性，而且操作简单，耗费人力物力较少，能实现实时预测，具有较好的适用性。
附图说明
图1是本发明的结构示意图；
具体实施方式
下面结合附图对本发明作进一步说明。
如图1所示，本发明的具体步骤为：
(1)采用钻屑法或应力计法在冲击危险区域测试煤体中的应力值，然后绘制应力分布曲线，通过对距离求导，得出应力梯度的分布情况，具体公式为；
上式中：ΔS为应力点σ1和应力点σ2之间的距离，Δσ为应力梯度；
(2)在冲击危险区域煤体的一侧放炮或撞击煤壁，通过设置在另一侧拾波器接收，然后进行CT层析成像技术得到煤体中的波速分布，通过对距离求导，得出波速梯度的分布情况，具体公式为；
上式中：ΔS是测量的波速点V1和测量的波速点V2之间的距离，VG为波速梯度；
(3)由于应力梯度与波速梯度之间的耦合性，根据下述线性公式，得出随着应力梯度增加，波速梯度也线性增加，通过确定煤体发生冲击地压的应力梯度临界值，得出波速梯度临界值VGm；
Vp＝a×σλ
上式中，Vp表示P波波速；a是常数，是修正值；σ为应力值；λ是指数，常量；
(4)在采煤机回采过程中，在煤壁一侧每半小时撞击一次，在另一侧通过拾波器接收波，进而得到煤体中波速的分布，通过公式实时得到煤体中的波速梯度VG；
(5)将实时检测得出的波速梯度VG与波速梯度临界值VGm进行比较，若VG小于VGm，则预测结果为安全；若VG大于VGm，则预测结果为危险，可及时采取相应的解危措施。
Claims (1)
1.一种地应力型冲击地压危险性多元信息耦合预测方法，其特征在于，其具体步骤为：
(1)采用钻屑法或应力计法在冲击危险区域测试煤体中的应力值，然后绘制应力分布曲线，通过对距离求导，得出应力梯度的分布情况，具体公式为；
<mrow>
<mi>&Delta;</mi>
<mi>&sigma;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&sigma;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>S</mi>
</mrow>
</mfrac>
</mrow>
上式中：ΔS为应力点A和应力点B之间的距离，σ1为应力点A的应力值，σ2为应力点B的应力值，Δσ为应力梯度；
(2)在冲击危险区域煤体的一侧放炮或撞击煤壁，通过设置在另一侧拾波器接收，然后进行CT层析成像技术得到煤体中的波速分布，通过对距离求导，得出波速梯度的分布情况，具体公式为；
<mrow>
<mi>V</mi>
<mi>G</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>S</mi>
</mrow>
</mfrac>
</mrow>
上式中：ΔS是波速点A和波速点B之间的距离，V1为波速点A的波速值，V2为波速点B的波速值，VG为波速梯度；
(3)由于应力梯度与波速梯度之间的耦合性，根据下述线性公式，得出随着应力梯度增加，波速梯度也线性增加，通过确定煤体发生冲击地压的应力梯度临界值，得出波速梯度临界值VGm；
Vp＝a×σλ
上式中，Vp表示P波波速；a是常数，是修正值；σ为应力值；λ是指数，常量；
(4)在采煤机回采过程中，在煤壁一侧每半小时撞击一次，在另一侧通过拾波器接收波，进而得到煤体中波速的分布，通过公式实时得到煤体中的波速梯度VG；
(5)将实时检测得出的波速梯度VG与波速梯度临界值VGm进行比较，若VG小于VGm，则预测结果为安全；若VG大于VGm，则预测结果为危险。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510375710.7A CN105045969B (zh) | 2015-06-30 | 2015-06-30 | 一种地应力型冲击地压危险性多元信息耦合预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510375710.7A CN105045969B (zh) | 2015-06-30 | 2015-06-30 | 一种地应力型冲击地压危险性多元信息耦合预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105045969A CN105045969A (zh) | 2015-11-11 |
| CN105045969B true CN105045969B (zh) | 2017-12-26 |
Family
ID=54452512
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510375710.7A Expired - Fee Related CN105045969B (zh) | 2015-06-30 | 2015-06-30 | 一种地应力型冲击地压危险性多元信息耦合预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105045969B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106501857B (zh) * | 2016-11-02 | 2018-11-27 | 中国矿业大学 | 一种煤矿巷道冲击地压危险性的声学监测方法 |
| CN107633146B (zh) * | 2017-09-28 | 2021-03-30 | 国家电网公司 | 一种不同尺度模型间地应力计算结果高精度转换的方法 |
| CN109798106B (zh) * | 2018-11-13 | 2022-09-20 | 辽宁工程技术大学 | 一种冲击地压危险性的预测方法及防治措施 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104653226A (zh) * | 2014-12-26 | 2015-05-27 | 中国矿业大学 | 一种基于应力梯度的煤矿冲击地压危险区域的划分方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2002951892A0 (en) * | 2002-10-09 | 2002-10-24 | Bhp Billiton Innovation Pty Ltd | Mining process and design |
-
2015
- 2015-06-30 CN CN201510375710.7A patent/CN105045969B/zh not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104653226A (zh) * | 2014-12-26 | 2015-05-27 | 中国矿业大学 | 一种基于应力梯度的煤矿冲击地压危险区域的划分方法 |
Non-Patent Citations (3)
| Title |
|---|
| Prevention and forcast of rock burst hazards in coal mines;DOU Lin-ming等;《Mining Science and Technology》;20091017;全文 * |
| 冲击倾向煤岩纵波波速与应力关系试验研究;巩思园等;《采矿与安全工程学报》;20120115;全文 * |
| 冲击危险性动态预测的震动波CT 技术研究;窦林名等;《煤炭学报》;20140519;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105045969A (zh) | 2015-11-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ghosh et al. | Application of underground microseismic monitoring for ground failure and secure longwall coal mining operation: a case study in an Indian mine | |
| CN105179018B (zh) | 一种煤矿冲击地压检测方法 | |
| CN104653226B (zh) | 一种基于应力梯度的煤矿冲击地压危险区域的划分方法 | |
| Zhou et al. | Seepage channel development in the crown pillar: Insights from induced microseismicity | |
| CN105891874B (zh) | 一种采动煤岩体突水微震监测方法 | |
| CN105741029B (zh) | 一种基于累积应力集中系数的冲击危险评价方法 | |
| Li et al. | Characteristics of microseismic waveforms induced by hydraulic fracturing in coal seam for coal rock dynamic disasters prevention | |
| Lu et al. | Microseismic low-frequency precursor effect of bursting failure of coal and rock | |
| Dey et al. | Prediction of blast-induced overbreak from uncontrolled burn-cut blasting in tunnels driven through medium rock class | |
| Lu et al. | Microseismic and acoustic emission effect on gas outburst hazard triggered by shock wave: a case study | |
| CN104880729B (zh) | 一种基于连续跟踪槽波信号的煤巷超前探测异常构造方法 | |
| CN105045969B (zh) | 一种地应力型冲击地压危险性多元信息耦合预测方法 | |
| CN106501857B (zh) | 一种煤矿巷道冲击地压危险性的声学监测方法 | |
| CN103104294A (zh) | 一种冲击地压的预测方法 | |
| CN204101748U (zh) | 组合式冲击矿压危险区域地震ct探测用液压式自补偿可控震源 | |
| Ma et al. | A method for numerical simulation based on microseismic information and the interpretation of hard rock fracture | |
| CN112012797A (zh) | 一种用于煤矿冲击危险卸压效果的评估方法 | |
| CN103775073B (zh) | 一种采掘工作面地应力分布特征探测方法 | |
| Tan et al. | Quantitative prop support estimation and remote monitor early warning for hard roof weighting at the Muchengjian Mine in China | |
| CN109239770A (zh) | 一种以co2炮为震源的工作面侧帮地震探测方法 | |
| Ishida et al. | Acoustic emission induced by progressive excavation of an underground powerhouse | |
| CN105756709A (zh) | 一种工作面顶板来压及破断的监测方法 | |
| CN109886600A (zh) | 监测预警煤层分叉合并区冲击地压危险性的方法 | |
| CN113126147A (zh) | 一种确定岩巷底板隐伏陷落柱空间形态的探测方法 | |
| Singh et al. | Review of IT enabled technologies in Indian mining industry for improved productivity & safety |
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 | ||
| CP02 | Change in the address of a patent holder | ||
| CP02 | Change in the address of a patent holder |
Address after: 221000 research and development center of mine safety technology and equipment at China University of Mining and Technology, Xuzhou, Xuzhou, Jiangsu Patentee after: China University of Mining & Technology Address before: 221000 Xuzhou high tech Industrial Development Zone,, Xuzhou Road, China University of Mining and Technology, mine safety technology and equipment research and Development Center Patentee before: China University of Mining & Technology |
|
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171226 Termination date: 20190630 |