CN116794725A - Method for correcting impact dangerous area division during stoping based on tunneling microseismic data - Google Patents
Method for correcting impact dangerous area division during stoping based on tunneling microseismic data Download PDFInfo
- Publication number
- CN116794725A CN116794725A CN202310758859.8A CN202310758859A CN116794725A CN 116794725 A CN116794725 A CN 116794725A CN 202310758859 A CN202310758859 A CN 202310758859A CN 116794725 A CN116794725 A CN 116794725A
- Authority
- CN
- China
- Prior art keywords
- impact
- tunneling
- microseismic
- coordinate system
- during
- 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
Links
- 230000005641 tunneling Effects 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000011435 rock Substances 0.000 claims abstract description 24
- 238000005065 mining Methods 0.000 claims abstract description 7
- 238000012544 monitoring process Methods 0.000 claims abstract description 7
- 238000011084 recovery Methods 0.000 claims description 10
- 238000000605 extraction Methods 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 4
- 238000010168 coupling process Methods 0.000 claims description 4
- 238000005192 partition Methods 0.000 claims description 4
- 238000007796 conventional method Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 2
- 239000003245 coal Substances 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The method for correcting the impact dangerous area division during stoping based on tunneling microseismic data comprises the following steps: the method comprises the steps of collecting microseismic data generated during roadway tunneling by using a mine microseismic monitoring system, converting a coordinate system of a working face and microseismic data in a mining engineering plane base map, enabling the X-axis direction of the coordinate system to be parallel to the trend of the working face, enabling the Y-axis direction of the coordinate system to be parallel to the trend of the working face, dividing a statistical area along the trend direction of the working face, enabling the size of the statistical area to be equal to the average footage during tunneling, counting the total microseismic energy of each area, dividing impact risk levels during tunneling according to rock burst early warning indexes, dividing impact risk levels during stoping by adopting a traditional method, integrating the impact risk levels during tunneling and stoping, and taking the maximum risk level of the two as the corrected impact risk level in the same area. The invention provides field data support for the impact dangerous area division and improves the accuracy and reliability of the impact dangerous area division.
Description
Technical Field
The invention relates to a method for dividing impact dangerous areas during correction and stoping based on tunneling microseismic data, and belongs to the technical field of coal mine safety exploitation.
Background
Rock burst is a mine dynamic phenomenon which causes great damage due to sudden release of elastic energy accumulated in coal rock mass, and the instantaneous burst of the rock burst can bring serious threat to mine safety. When mining a working surface with rock burst risk, impact risk evaluation and impact risk classification are required to be carried out on the working surface so as to prevent rock burst accidents and ensure mine safety production.
At present, most of methods for dividing impact dangerous areas are based on theoretical geology and mining conditions, and lack of guidance of on-site data, so that the accuracy and reliability of the division of the impact dangerous areas are required to be improved.
Disclosure of Invention
The invention provides a method for dividing an impact dangerous area during correction and extraction based on tunneling microseismic data, which can provide on-site data support for dividing the impact dangerous area, improve the accuracy and reliability of dividing the impact dangerous area and provide powerful guarantee for mine safety production.
In order to achieve the above purpose, the invention provides a method for correcting the impact danger area division during stoping based on tunneling microseismic data, which comprises the following steps:
(1) Acquiring microseismic data generated during roadway tunneling by using a mine microseismic monitoring system;
(2) Converting a coordinate system of working surfaces and microseismic data in a mining engineering plane base map, so that the X-axis direction of the coordinate system is parallel to the trend of the working surfaces, and the Y-axis direction of the coordinate system is parallel to the trend of the working surfaces;
(3) Dividing a statistical region along the trend direction of the working face, wherein the dimension of the statistical region is equal to the average footage during tunneling, and counting the total microseismic energy of each partition;
(4) Dividing a working surface into two parts along a trend direction, respectively counting the energy falling into each grid of the two parts, and dividing impact danger levels during tunneling according to impact ground pressure early warning indexes;
(5) Dividing impact danger levels during stoping by adopting a traditional method;
(6) And integrating impact risk levels during tunneling and extraction, and taking the maximum risk level of the two as the corrected impact risk level in the same area.
Further, the formula for converting coordinates in the step (2) is as follows:
x′=(x-x o )cos(α)+(y-y 0 )sin(α)
y′=(y-y 0 )cos(α)-(x-x o )sin(α)
wherein, x and y represent x and y coordinates of the microseismic event under an original coordinate system; alpha represents a clockwise included angle between the new coordinate and the original coordinate system; x is x 0 、y 0 The coordinate origin point A is the x and y coordinates of the new coordinate system; x ', y' represent the x, y coordinates of the microseismic event in the new coordinate system.
Further, in the step (4), the method for counting the energy within each grid is as follows: the working face is divided into two parts along the inclined direction, and the mine earthquake near the side of the transportation roadway is regarded as falling into the grid of the transportation roadway, and the mine earthquake near the side of the return roadway is regarded as falling into the grid of the return roadway.
Further, in the step (4), the rock burst early warning index is the total energy of micro-vibration occurring near the tunneling working face every 24 hours, and the impact danger level dividing method comprises the following steps: the total microseismic energy in each region is more than 0.75 times of the rock burst early warning index, and the risk of strong impact is considered; the total micro-vibration energy in each area is regarded as medium impact danger between 0.5 and 0.75 times of the rock burst early warning index; the total micro-vibration energy in each area is regarded as weak impact danger between 0.25 and 0.5 times of the rock burst early warning index; the total energy of micro-vibration in each area is smaller than 0.25 times of the rock burst early warning index, and no impact danger exists.
Further, in the step (5), the conventional method for classifying the impact risk level during the recovery is one of a stress superposition method and a multi-factor coupling method.
The invention collects microseismic data generated during roadway tunneling by using a mine microseismic monitoring system, converts a coordinate system of a working face and microseismic data in a mining engineering plane base map, enables the X-axis direction of the coordinate system to be parallel to the trend of the working face, the Y-axis direction of the coordinate system to be parallel to the trend of the working face, divides a statistical area along the trend direction of the working face, the size of the statistical area is equal to the average footage during tunneling, counts the microseismic total energy of each area, divides the impact risk level during tunneling according to rock burst early warning indexes, divides the impact risk level during stoping by adopting a traditional method, synthesizes the impact risk level during tunneling and stoping, and takes the maximum risk level of the two as the corrected impact risk level in the same area. The method provides field data support for the division of the impact dangerous areas, improves the accuracy and reliability of the division of the impact dangerous areas, and provides powerful guarantee for the safety production of mines, and the method is simple, convenient and high in operability.
Drawings
FIG. 1 is a graph showing microseismic profiles during a 250107-1 face tunneling in an embodiment of the present invention;
FIG. 2 is a graph showing microseismic profiles during tunneling after a 250107-1 face coordinate system transformation in an embodiment of the present invention;
FIG. 3 is a 250107-1 face impact risk classification based on heading microseismic data in an embodiment of the present invention;
FIG. 4 is a 250107-1 face extraction time impact hazard classification in accordance with an embodiment of the invention;
FIG. 5 is a 250107-1 face corrected during recovery impact hazard classification in accordance with an embodiment of the invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
A method for correcting impact dangerous area division during stoping based on tunneling microseismic data comprises the following steps:
(1) Acquiring microseismic data generated during roadway tunneling by using a mine microseismic monitoring system;
(2) Converting a coordinate system of working surfaces and microseismic data in a mining engineering plane base map, so that the X-axis direction of the coordinate system is parallel to the trend of the working surfaces, and the Y-axis direction of the coordinate system is parallel to the trend of the working surfaces;
(3) Dividing a statistical region along the trend direction of the working face, wherein the dimension of the statistical region is equal to the average footage during tunneling, and counting the total microseismic energy of each partition;
(4) Dividing a working surface into two parts along a trend direction, respectively counting the energy falling into each grid of the two parts, and dividing impact danger levels during tunneling according to impact ground pressure early warning indexes;
(5) Dividing impact danger levels during stoping by adopting a traditional method;
(6) And integrating impact risk levels during tunneling and extraction, and taking the maximum risk level of the two as the corrected impact risk level in the same area.
Further, the formula for converting coordinates in the step (2) is as follows:
x′=(x-x o )cos(α)+(y-y 0 )sin(α)
y′=(y-y 0 )cos(α)-(x-x o )sin(α)
wherein, x and y represent x and y coordinates of the microseismic event under an original coordinate system; alpha represents a clockwise included angle between the new coordinate and the original coordinate system; x is x 0 、y 0 The coordinate origin point A is the x and y coordinates of the new coordinate system; x ', y' represent the x, y coordinates of the microseismic event in the new coordinate system.
Further, in the step (4), the method for counting the energy within each grid is as follows: the working face is divided into two parts along the inclined direction, and the mine earthquake near the side of the transportation roadway is regarded as falling into the grid of the transportation roadway, and the mine earthquake near the side of the return roadway is regarded as falling into the grid of the return roadway.
Further, in the step (4), the rock burst early warning index is the total energy of micro-vibration occurring near the tunneling working face every 24 hours, and the impact danger level dividing method comprises the following steps: the total microseismic energy in each region is more than 0.75 times of the rock burst early warning index, and the risk of strong impact is considered; the total micro-vibration energy in each area is regarded as medium impact danger between 0.5 and 0.75 times of the rock burst early warning index; the total micro-vibration energy in each area is regarded as weak impact danger between 0.25 and 0.5 times of the rock burst early warning index; the total energy of micro-vibration in each area is smaller than 0.25 times of the rock burst early warning index, and no impact danger exists.
As a preferred embodiment, in the step (5), the conventional method for classifying the impact risk level during the recovery is one of a stress superposition method and a multi-factor coupling method.
Examples:
the invention divides the impact danger grade of a mine 250107-1 working face during stoping, and comprises the following specific steps:
(1) The method comprises the steps that a microseismic monitoring system is installed in a coal mine, microseismic data generated during tunneling of a working face of 250107-1 is collected, and microseismic distribution is shown in figure 1;
(2) Converting the coordinate system of the working surface and the microseismic data in the base map, so that the X-axis direction of the coordinate system is parallel to the trend of the working surface, the Y-axis direction of the coordinate system is parallel to the trend of the working surface, and the converted working surface and microseismic distribution are shown in figure 2;
(3) Dividing a statistical region along the trend direction of the working face, wherein the dimension of the statistical region is equal to the average footage during tunneling, and counting the total microseismic energy of each partition;
(4) Obtaining 250107-1 working face rock burst vibration total energy monitoring and early warning index of 8.72e 4 J, dividing impact risks during tunneling into four grades according to rock burst early warning indexes, wherein the four grades are respectively strong impact risks:>6.54e 4 J. medium impact risk: 4.36e 4 J~6.54e 4 J. Weak impact hazard: 2.18e 4 J~4.36e 4 J. No risk of impact:<2.18e 4 J;
(5) Dividing a working surface into two parts along the inclined direction, wherein the mine earthquake near the side of the transportation roadway is regarded as falling into the grid of the transportation roadway, and the mine earthquake near the side of the return roadway is regarded as falling into the grid of the return roadway; counting the total energy falling into each grid with the width of 5m and the length of 8m, and dividing the impact danger areas during tunneling according to the impact danger indexes determined in the step (4), wherein the dividing results are shown in the table 1 and the figure 3:
TABLE 1 impact hazard classification based on tunneling microseismic data
(5) The impact danger zone during recovery is partitioned according to the multi-factor coupling method, and the partitioning results are shown in table 2 and fig. 4:
TABLE 2 impact hazard classification during recovery
(6) The impact risk level during the tunneling and the extraction is synthesized, the maximum risk level of the tunneling and the extraction is taken as the corrected impact risk level in the same area, and the corrected impact risk level is shown in table 3 and fig. 5:
TABLE 3 impact hazard classification after correction
Claims (5)
1. The method for correcting the impact dangerous area division during stoping based on tunneling microseismic data is characterized by comprising the following steps:
(1) Acquiring microseismic data generated during roadway tunneling by using a mine microseismic monitoring system;
(2) Converting a coordinate system of working surfaces and microseismic data in a mining engineering plane base map, so that the X-axis direction of the coordinate system is parallel to the trend of the working surfaces, and the Y-axis direction of the coordinate system is parallel to the trend of the working surfaces;
(3) Dividing a statistical region along the trend direction of the working face, wherein the dimension of the statistical region is equal to the average footage during tunneling, and counting the total microseismic energy of each partition;
(4) Dividing a working surface into two parts along a trend direction, respectively counting the energy falling into each grid of the two parts, and dividing impact danger levels during tunneling according to impact ground pressure early warning indexes;
(5) Dividing impact danger levels during stoping by adopting a traditional method;
(6) And integrating impact risk levels during tunneling and extraction, and taking the maximum risk level of the two as the corrected impact risk level in the same area.
2. The method for classifying the impact danger zone during the correction and recovery based on the tunneling microseismic data according to claim 1, wherein the formula of the transformation coordinates in the step (2) is:
x′=(x-x o )cos(α)+(y-y 0 )sin(α)
y′=(y-y 0 )cos(α)-(x-x o )sin(α)
wherein, x and y represent x and y coordinates of the microseismic event under an original coordinate system; alpha represents a clockwise included angle between the new coordinate and the original coordinate system; x is x 0 、y 0 The coordinate origin point A is the x and y coordinates of the new coordinate system; x ', y' represent the x, y coordinates of the microseismic event in the new coordinate system.
3. The method for classifying the impact danger zone during the recovery based on the tunneling microseismic data according to claim 1 or 2, wherein in the step (4), the method for counting the energy falling into each grid is as follows: the working face is divided into two parts along the inclined direction, and the mine earthquake near the side of the transportation roadway is regarded as falling into the grid of the transportation roadway, and the mine earthquake near the side of the return roadway is regarded as falling into the grid of the return roadway.
4. The method for classifying impact danger areas during recovery based on tunneling microseismic data according to claim 3, wherein in the step (4), the rock burst early warning index is total microseismic energy occurring near the tunneling working face every 24 hours, and the classifying method of the impact danger level is as follows: the total microseismic energy in each region is more than 0.75 times of the rock burst early warning index, and the risk of strong impact is considered; the total micro-vibration energy in each area is regarded as medium impact danger between 0.5 and 0.75 times of the rock burst early warning index; the total micro-vibration energy in each area is regarded as weak impact danger between 0.25 and 0.5 times of the rock burst early warning index; the total energy of micro-vibration in each area is smaller than 0.25 times of the rock burst early warning index, and no impact danger exists.
5. The method for classifying an impact risk area during recovery based on tunneling microseismic data according to claim 4, wherein in the step (5), the conventional method for classifying the impact risk level during recovery is one of a stress superposition method and a multi-factor coupling method.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310758859.8A CN116794725A (en) | 2023-06-26 | 2023-06-26 | Method for correcting impact dangerous area division during stoping based on tunneling microseismic data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310758859.8A CN116794725A (en) | 2023-06-26 | 2023-06-26 | Method for correcting impact dangerous area division during stoping based on tunneling microseismic data |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116794725A true CN116794725A (en) | 2023-09-22 |
Family
ID=88034058
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310758859.8A Pending CN116794725A (en) | 2023-06-26 | 2023-06-26 | Method for correcting impact dangerous area division during stoping based on tunneling microseismic data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116794725A (en) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1293364A1 (en) * | 1985-09-13 | 1987-02-28 | Кузбасский Политехнический Институт | Method of mining shock-hazardous coal seams |
WO2007114734A1 (en) * | 2006-03-31 | 2007-10-11 | Yury Alexeevich Baurov | Method for the short-range forecast of the risk seismic hazard of a seismically hazardous region |
CN103306722A (en) * | 2013-06-21 | 2013-09-18 | 中国矿业大学 | Micro-seismic multi-dimensional information comprehensive region detection and evaluation method for impact danger region |
CN111271128A (en) * | 2020-01-20 | 2020-06-12 | 天地科技股份有限公司 | Coal rock mass breaking danger monitoring method, device and system |
CN112483174A (en) * | 2020-11-20 | 2021-03-12 | 中国矿业大学 | Arrangement method of tunneling working face impact dangerous vibration wave CT inversion system |
CN115220092A (en) * | 2022-07-15 | 2022-10-21 | 陕西正通煤业有限责任公司 | Microseismic statistical method for determining advanced impact danger range of working face |
CN115392061A (en) * | 2022-10-28 | 2022-11-25 | 北京科技大学 | Method for evaluating static and dynamic coupling of rock burst danger |
CN115511379A (en) * | 2022-10-28 | 2022-12-23 | 北京科技大学 | Dynamic dividing method and device for rock burst dangerous area |
-
2023
- 2023-06-26 CN CN202310758859.8A patent/CN116794725A/en active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1293364A1 (en) * | 1985-09-13 | 1987-02-28 | Кузбасский Политехнический Институт | Method of mining shock-hazardous coal seams |
WO2007114734A1 (en) * | 2006-03-31 | 2007-10-11 | Yury Alexeevich Baurov | Method for the short-range forecast of the risk seismic hazard of a seismically hazardous region |
CN103306722A (en) * | 2013-06-21 | 2013-09-18 | 中国矿业大学 | Micro-seismic multi-dimensional information comprehensive region detection and evaluation method for impact danger region |
CN111271128A (en) * | 2020-01-20 | 2020-06-12 | 天地科技股份有限公司 | Coal rock mass breaking danger monitoring method, device and system |
CN112483174A (en) * | 2020-11-20 | 2021-03-12 | 中国矿业大学 | Arrangement method of tunneling working face impact dangerous vibration wave CT inversion system |
CN115220092A (en) * | 2022-07-15 | 2022-10-21 | 陕西正通煤业有限责任公司 | Microseismic statistical method for determining advanced impact danger range of working face |
CN115392061A (en) * | 2022-10-28 | 2022-11-25 | 北京科技大学 | Method for evaluating static and dynamic coupling of rock burst danger |
CN115511379A (en) * | 2022-10-28 | 2022-12-23 | 北京科技大学 | Dynamic dividing method and device for rock burst dangerous area |
Non-Patent Citations (2)
Title |
---|
刘旭东;张玉良;田向辉;: "基于微震监测的近直立特厚煤层冲击危险区域划分", 煤炭工程, no. 01, 20 January 2020 (2020-01-20), pages 154 - 157 * |
林飞;王萍;陈娜;王娟娟;: "微震分析软件在冲击地压监测预警中的应用", 陕西煤炭, no. 04, 12 August 2020 (2020-08-12), pages 85 - 89 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105741029B (en) | A kind of danger of burst evaluation method based on cumulative stress coefficient of concentration | |
CN104653226A (en) | Stress-gradient-based method for dividing coal impact ground pressure danger area | |
CN109798106B (en) | Method for predicting risk of rock burst and prevention and treatment measures | |
CN109595036A (en) | A kind of method for early warning of Mine Ground Pressure Disaster | |
Gao et al. | The dynamic failure mechanism of coal and gas outbursts and response mechanism of support structure | |
CN108594295A (en) | A kind of coal seam unloading pressure by blasting effect evaluation method | |
CN113962466B (en) | Dynamic prediction method for impact danger in bending region of rock burst mine based on microseism | |
CN112377257B (en) | Working face mining advance influence range determining method based on microseismic monitoring | |
CN105156110A (en) | Dynamic disaster mine district roadway structure and method for developing same | |
CN115220092A (en) | Microseismic statistical method for determining advanced impact danger range of working face | |
CN115169948A (en) | Coal face overlying rock separation layer water inrush risk prediction method and safe mining method | |
CN111768308A (en) | Mining strength evaluation method based on rock burst monitoring data analysis | |
CN116794725A (en) | Method for correcting impact dangerous area division during stoping based on tunneling microseismic data | |
CN104832163A (en) | Method for monitoring rock burst dangerousness in coal underground mining process | |
Ellenberger et al. | Three dimensional microseismic monitoring of a Utah longwall | |
CN107990798A (en) | The definite method of mountain area civil buildings blasting vibration allowable safety distance | |
CN112377258A (en) | Method for evaluating mining influence degree between adjacent mining areas and mining working faces | |
CN110043317B (en) | Mine disaster multi-parameter local danger area judgment and early warning method | |
Prusek et al. | Powered support selection for longwall workings in dynamic load conditions | |
AU2023202324A1 (en) | Method of steeply-inclined horizontal section top-coal caving for pressure relief | |
CN105952455A (en) | Safety assessing method after strip mine caving process processing of complicated cavity | |
Iannacchione et al. | Safer mine layouts for underground stone mines subjected to excessive levels of horizontal stress | |
CN110807269B (en) | Fault activation tendency analysis method based on critical angle | |
CN113434819A (en) | Method for determining influence time and distance of working face mining on mine earthquake activity of goaf | |
Xu et al. | Study on the reasonable cover thickness of a subsea tunnel with the numerical calculation criterion method |
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 |