CN115220092A - Microseismic statistical method for determining advanced impact danger range of working face - Google Patents

Microseismic statistical method for determining advanced impact danger range of working face Download PDF

Info

Publication number
CN115220092A
CN115220092A CN202210830765.2A CN202210830765A CN115220092A CN 115220092 A CN115220092 A CN 115220092A CN 202210830765 A CN202210830765 A CN 202210830765A CN 115220092 A CN115220092 A CN 115220092A
Authority
CN
China
Prior art keywords
microseismic
determining
energy
working face
coordinate system
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
CN202210830765.2A
Other languages
Chinese (zh)
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.)
China University of Mining and Technology CUMT
Shaanxi Zhengtong Coal Industry Co Ltd
Original Assignee
China University of Mining and Technology CUMT
Shaanxi Zhengtong Coal Industry 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 China University of Mining and Technology CUMT, Shaanxi Zhengtong Coal Industry Co Ltd filed Critical China University of Mining and Technology CUMT
Priority to CN202210830765.2A priority Critical patent/CN115220092A/en
Publication of CN115220092A publication Critical patent/CN115220092A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/288Event detection in seismic signals, e.g. microseismics

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Emergency Management (AREA)
  • Business, Economics & Management (AREA)
  • Acoustics & Sound (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

A microseismic statistical method for determining the advance impact danger range of a working face comprises the following steps: collecting natural microseismic signals in the working face extraction process by using a mine microseismic monitoring system, and determining the position and energy of a seismic source; selecting any coordinate point on a mining engineering plane graph as a coordinate origin of a new coordinate system, drawing the new coordinate system, introducing microseismic data, obtaining coordinates of a microseismic event under the new coordinate system through a coordinate conversion formula, and obtaining an abscissa of a daily pushing position of a working surface under the new coordinate system; determining the distance of each daily microseismic event from the working surface propulsion position; and (3) counting microseismic events in a range before and after the propulsion position, dividing a counting area at equal intervals, counting total microseismic energy and frequency in each area, drawing a curve chart of total microseismic energy and frequency variation, dividing total microseismic energy and frequency danger levels, and determining the range of the advanced impact danger area. The method can determine the area range with the advance impact danger in the working face extraction period, and is high in reliability and easy to operate.

Description

Microseismic statistical method for determining advanced impact danger range of working face
Technical Field
The invention relates to a microseismic statistic method for determining the advanced impact danger range of a working face, belonging to the technical field of coal mine safety mining.
Background
The microseism monitoring method is used as a leading-edge technology in the aspect of monitoring and forecasting dynamic disasters such as impact mine pressure, and determines and analyzes the position and energy of vibration by recording the vibration spectrum of the vibration caused by mining activities. Along with the increase of the mining depth of the mine, the influence of coal rock dynamic disasters on the safe production of the mine is increasingly prominent. At present, microseismic monitoring systems are installed in more mines, particularly in dangerous mines with disasters such as rock burst and the like, so that monitoring of microseismic events of the mines is realized.
Most of the existing micro-seismic monitoring systems are complex in design and high in cost, and cannot accurately monitor the advanced impact dangerous area range during the working face recovery period.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides the microseismic statistical method for determining the advanced impact danger range of the working face, which can determine the area range with the advanced impact danger during the mining period of the working face, and has the advantages of high reliability, easy operation and low cost.
In order to achieve the above object, the present invention provides a microseismic statistical method for determining a working face advanced impact risk range, comprising the following steps:
(1) Acquiring natural micro-seismic signals in the working face extraction process by using a mine micro-seismic monitoring system, and determining a seismic source position and seismic source energy;
(2) Selecting any coordinate point A (x) on the plane drawing of the excavation project 0 ,y 0 ) Drawing a new coordinate system as the coordinate origin of the new coordinate system, and sitting the new coordinate systemThe x axis of the mark system is parallel to the propelling direction of the working surface;
(3) Introducing microseismic data, and processing the microseismic data through a coordinate conversion formula to obtain the coordinate of the microseismic event in a new coordinate system;
(4) Determining a working face propulsion starting point coordinate B, and introducing working face footage data to obtain an x coordinate of a working face daily propulsion position in a new coordinate system;
(5) Determining the moving distance of each day of microseismic event from the advancing position of the working surface;
(6) Dividing statistical regions at equal intervals along the direction of the working surface, and determining total microseismic energy and microseismic frequency in each statistical region by using cumulative addition;
(7) Drawing a variation curve graph of the total microseismic energy and the microseismic frequency by taking the moving distance from the working surface to the propelling position determined in the step (5) as an abscissa and taking the total microseismic energy and the microseismic frequency in each statistical area counted in the step (6) as an ordinate;
(8) And (4) dividing the total microseismic energy and the microseismic frequency danger level according to the total microseismic energy and microseismic frequency change curve chart obtained in the step (7) and determining the advanced impact dangerous area range.
Further, the method for determining the position of the seismic source and the energy of the seismic source in the step (1) comprises the following steps: the SOS micro-seismic monitoring system is installed in a coal mine and comprises an underground acquisition recording device and a plurality of single-component probes installed underground, the single-component probes are installed on an anchor rod perpendicular to a roadway bottom plate, the mine micro-seismic monitoring system is used for acquiring a plurality of natural mine seismic signals in the stoping process of a working face, and the natural mine seismic signals are processed to obtain the corresponding mine seismic source positions and the corresponding mine seismic source energy.
Further, the x axis of the new coordinate system drawn in the step (2) is parallel to the advancing direction of the working surface, and forms a clockwise included angle alpha with the x axis of the old coordinate system.
Further, in the step (3), the coordinate conversion formula is as follows:
x′ ij =(x ij -x o )cos(α)+(y ij -y 0 )sin(α),
y′ ij =(y ij -y 0 )cos(α)-(x ij -x o )sin(α);
in the formula, x ij 、y ij Respectively representing x and y coordinates of the j microseismic event on the ith day in an old coordinate system; alpha represents a clockwise included angle between the new coordinate and the x axis of the old coordinate system; x is the number of o 、y 0 Respectively representing x and y coordinates of a new coordinate system coordinate origin point A; x' ij 、y′ ij Representing the x, y coordinates of the j microseismic event on day i in the new coordinate system.
Further, in the step (4), the determination method of the x coordinate of the daily advancing position of the working surface in the new coordinate system is as follows: selecting a working face lower cutting hole B point as a working face advancing starting point position, introducing working face footage data, and calculating by adopting the following formula:
x′ Bi =(x B -x 0 )cos(α)+(y B -y 0 )sin(α)±L i
x 'in the formula' Bi Representing the advancing position of the working surface on the ith day; x is the number of B 、y B Coordinates of a point B of a lower cutting eye of the selected working surface; l is a radical of an alcohol i And (4) indicating the cumulative ruler feeding of the day i, and selecting the + calculation when the advancing direction of the working surface is consistent with the positive direction of the x axis of the new coordinate, and otherwise selecting the-calculation.
Further, in the step (5), the calculation formula of the strike distance of the daily microseismic event from the working surface propulsion position is as follows:
d ij =x′ Bi -x′ ij
in the formula, d ij Indicating the j microseismic event strike distance from the face advancement position on day i.
Further, in the step (6), the total microseismic energy and the microseismic frequency in each divided statistical region are calculated according to the following formula:
E n =e 1n +e 2n +e 3n +…+e mn
Q n =q 1n +q 2n +q 3n +…+q mn
in the formula, E n Is one by oneEach statistical area is the total microseismic energy in the subarea; q n The microseismic frequency in each statistical region, namely the subarea; e.g. of the type in I =1,2,3, …, m for the energy of the ith microseismic event of the nth zone; q. q.s in For the ith microseismic event for the nth zone, i =1,2,3, …, m; n is each partition number.
Further, in the step (8), the total microseismic energy and microseismic frequency danger level is determined according to the mining condition of the working face and the total microseismic energy and microseismic frequency danger standard value determined by the monitoring and early warning index.
The method comprises the steps of collecting natural microseismic signals in the working face extraction process by using a mine microseismic monitoring system, determining the position and energy of a seismic source, selecting any coordinate point on a microseismic positioning base map as the origin of coordinates of a new coordinate system, drawing the new coordinate system, introducing microseismic data, obtaining the coordinates of microseismic time under the new coordinate system through a coordinate conversion formula, obtaining the abscissa of a daily propulsion position of the working face under the new coordinate system, and determining the trend distance of a daily microseismic event from the propulsion position of the working face; and (3) carrying out statistics on microseismic events within the range before and after the propulsion position, dividing statistical areas at equal intervals, carrying out statistics to determine total microseismic energy and microseismic frequency in each divided statistical area, drawing a curve graph of total microseismic energy and microseismic frequency change, dividing total microseismic energy and microseismic frequency danger levels according to the obtained curve graph of total microseismic energy and microseismic frequency change, and determining the range of the advanced impact dangerous area. The method determines the area range with the advance impact danger during the stoping period of the working face, has high reliability, is easy to operate, and is convenient for computer programming.
Drawings
FIG. 1 is a schematic diagram of the new coordinate system of the present invention;
FIG. 2 is a diagram illustrating the partitioning of statistical regions according to an embodiment of the present invention;
FIG. 3 is a schematic illustration of a work surface propulsion operation in an embodiment of the present invention;
FIG. 4 is a graph of the total microseismic energy and the frequency of the microseisms in an embodiment of the present invention.
Detailed Description
The invention will be further explained with reference to the drawings.
A microseismic statistical method for determining the advance impact danger range of a working face comprises the following steps:
(1) Acquiring natural micro-seismic signals in the working face extraction process by using a mine micro-seismic monitoring system, and determining a seismic source position and seismic source energy;
(2) Selecting any coordinate point A (x) on the plane drawing of the excavation project 0 ,y 0 ) Drawing a new coordinate system as the coordinate origin of the new coordinate system, wherein the x axis of the new coordinate system is parallel to the propelling direction of the working surface;
(3) Introducing microseismic data, and processing the microseismic data through a coordinate conversion formula to obtain the coordinates of the microseismic event in a new coordinate system;
(4) Determining a working face propulsion starting point coordinate B, and introducing working face footage data to obtain an x coordinate of a working face daily propulsion position in a new coordinate system;
(5) Determining the moving distance between each daily microseismic event and the advancing position of the working surface;
(6) Dividing statistical regions at equal intervals along the direction of the working surface, and determining total microseismic energy and microseismic frequency in each statistical region by using cumulative addition;
(7) Drawing a variation curve graph of the total microseismic energy and the microseismic frequency by taking the moving distance from the working surface to the propelling position determined in the step (5) as an abscissa and taking the total microseismic energy and the microseismic frequency in each statistical area counted in the step (6) as an ordinate;
(8) And (4) dividing the total microseismic energy and the microseismic frequency danger level according to the total microseismic energy and microseismic frequency change curve chart obtained in the step (7) and determining the advanced impact dangerous area range.
As a preferred embodiment of the present invention, the method for determining the position and energy of the seismic source in step (1) comprises: the SOS micro-seismic monitoring system is installed in a coal mine and comprises an aboveground acquisition and recording device and a plurality of single-component probes installed underground, the single-component probes are installed on anchor rods perpendicular to a roadway bottom plate, the mine micro-seismic monitoring system is used for acquiring a plurality of natural mine seismic signals in the working face stoping process, and the natural mine seismic signals are processed to obtain the corresponding mine seismic source positions and energy thereof.
In one embodiment, the x axis of the new coordinate system drawn in step (2) is parallel to the working surface advancing direction, and has a clockwise angle α with the x axis of the old coordinate system.
Specifically, in the step (3), the coordinate conversion formula is:
x′ ij =(x ij -x o )cos(α)+(y ij -y 0 )sin(α),
y′ ij =(y ij -y 0 )cos(α)-(x ij -x o )sin(α);
in the formula, x ij 、y ij Respectively representing x and y coordinates of the j microseismic event on the ith day in an old coordinate system; alpha represents a clockwise included angle between the new coordinate and the x axis of the old coordinate system; x is the number of o 、y 0 Respectively representing x and y coordinates of a new coordinate system origin point A; x' ij 、y′ ij Representing the x, y coordinates of the j microseismic event on day i in the new coordinate system.
In a preferred embodiment, in step (4), the x coordinate of the daily advancing position of the working surface in the new coordinate system is determined by: selecting a point B of a lower cutting hole of the working face as a position of a propelling starting point of the working face, introducing the footage data of the working face, and calculating by adopting the following formula:
x′ Bi =(x B -x 0 )cos(α)+(y B -y 0 )sin(α)±L i
in formula (II), x' Bi Representing the working surface propelling position on the ith day; x is the number of B 、y B Coordinates of a point B of a lower cutting eye of the selected working surface; l is i And (4) indicating the cumulative footage of the ith day, and selecting a plus calculation when the advancing direction of the working face is consistent with the positive direction of the x axis of the new coordinate, otherwise selecting a minus calculation.
Specifically, in the step (5), the calculation formula of the strike distance of the daily microseismic event from the working surface propulsion position is as follows:
d ij =x′ Bi -x′ ij
in the formula (d) ij Indicating the j microseismic event strike distance from the face advancement position on day i.
Specifically, in the step (6), the total microseismic energy and the microseismic frequency in each divided statistical region are calculated according to the following formula:
E n =e 1n +e 2n +e 3n +…+e mn
Q n =q 1n +q 2n +q 3n +…+q mn
in the formula, E n The total microseismic energy of each statistical area, namely the subarea; q n The microseismic frequency in each statistical region, namely a subarea; e.g. of a cylinder in For the energy of the ith microseismic event of the nth partition, i =1,2,3, …, m; q. q.s in For the ith microseismic event for the nth zone, i =1,2,3, …, m; n is the number of each partition.
In a preferred embodiment, in step (8), the total microseismic energy and microseismic frequency risk level is determined according to the total microseismic energy and microseismic frequency risk standard value determined by the working face mining condition and the monitoring and early warning index.
Example (b):
(1) By installing an SOS (seismic isolation system), collecting a plurality of natural mine seismic signals in the recovery process of the 7301 working face by using a mine micro-seismic monitoring system, determining the position and energy of each natural mine seismic source by using a known processing method for each natural mine seismic signal, and obtaining 7301 working face micro-seismic data;
(2) As shown in fig. 1, a point a (20401245.1220, 3916327.9279) is selected as an origin of a new coordinate system in the microseismic location base map, an x axis of the new coordinate system is parallel to the propulsion direction of the 7301 working surface, and an included angle alpha clockwise from the horizontal direction is 28.5 degrees;
(3) Introducing '7301 working face microseismic data', and determining coordinates x 'of microseismic event in new coordinate system by using coordinate conversion formula' ij =(x ij -20401245.1220)cos(28.5°)+(y ij -3916327.9279)sin(28.5°),y′ ij =(y ij -3916327.9279)cos(28.5°)-(x ij -20401245.1220)sin(28.5°) Wherein x is ij 、y ij 7301 working surface day i the j microseismic event;
(4) Selecting a working face undercut hole B point (20402687.0159, 3917116.2743) as a working face advancing starting point;
(5) Introducing '7301 working face footage data' to determine the x coordinate of the daily working face propelling position in the new coordinate system, wherein the propelling direction of the working face is opposite to the positive direction of the x in the new coordinate system, so that the following formula is selected as a 'minus' calculation: x' Bi =(x B -x 0 )cos(α)+(y B -y 0 )sin(α)-L i =(20402687.0159-20401245.1220)cos(28.5°)+(3917116.2743-3916327.9279)sin(28.5°)-L i =1267.16+376.16-L i =1643.32-L i Wherein L is i Cumulative footage for the working face day i;
(6) As shown in fig. 2 and fig. 3, according to the distribution of microseisms, preferably counting the mineral earthquakes within 800m range before and after the 7301 working face propulsion position, dividing the strike direction into a plurality of regions at equal intervals of 10m, and counting the total energy and frequency of the mineral earthquakes in each region by using cumulative methods; the calculation formula is as follows:
E n =e 1n +e 2n +e 3n +…+e mn
Q n =q 1n +q 2n +q 3n +…+q mn
E n the unit of the total microseismic energy in each statistical area, namely the subarea, is J; q n The microseismic frequency in each statistical region, namely the subarea; e.g. of the type in The energy of the ith microseismic event of the nth zone is in units of J, i =1,2,3, …, m; q. q.s in I =1,2,3, …, m for the frequency of the ith microseismic event for the nth zone; n is the number of each partition;
(7) Drawing a variation curve graph of the total microseismic energy and the microseismic frequency by taking the distance according to the propulsion position as an abscissa, taking the total microseismic energy and the microseismic frequency in each statistical area, namely each subarea, as an ordinate and the propulsion position as an origin of coordinates, as shown in fig. 4;
(8) Combining with 7301 working face microseismic early warning index, total energy of subareas is 3 multiplied by 10 5 J energy,The frequency is 100 microseismic events as the impact risk division standard, and as can be seen from fig. 4, the 7301 working surface has impact risk in the range of 160m ahead.

Claims (8)

1. A microseismic statistical method for determining the advance impact risk range of a working face is characterized by comprising the following steps:
(1) Collecting natural microseismic signals in the working face extraction process by using a mine microseismic monitoring system, and determining the seismic source position and the seismic source energy;
(2) Selecting any coordinate point A (x) on the plane drawing of the excavation project 0 ,y 0 ) Drawing a new coordinate system as the coordinate origin of the new coordinate system, wherein the x axis of the new coordinate system is parallel to the propelling direction of the working surface;
(3) Introducing microseismic data, and processing the microseismic data through a coordinate conversion formula to obtain the coordinates of the microseismic event in a new coordinate system;
(4) Determining a working face propulsion starting point coordinate B, and introducing working face footage data to obtain an x coordinate of a working face daily propulsion position in a new coordinate system;
(5) Determining the moving distance between each daily microseismic event and the advancing position of the working surface;
(6) Dividing statistical regions at equal intervals along the direction of the working surface, and determining total microseismic energy and microseismic frequency in each statistical region by using cumulative addition;
(7) Drawing a variation curve graph of the total microseismic energy and the microseismic frequency by taking the moving distance from the working surface to the propelling position determined in the step (5) as an abscissa and taking the total microseismic energy and the microseismic frequency in each statistical area counted in the step (6) as an ordinate;
(8) And (4) dividing the total microseismic energy and the microseismic frequency danger level according to the total microseismic energy and microseismic frequency change curve chart obtained in the step (7) and determining the advanced impact dangerous area range.
2. The microseismic statistic method for determining the risk range of the leading impact of the working face as claimed in claim 1 wherein the method for determining the source position and the source energy in step (1) is: the SOS micro-seismic monitoring system is installed in a coal mine and comprises an aboveground acquisition and recording device and a plurality of single-component probes installed underground, the single-component probes are installed on anchor rods perpendicular to a roadway bottom plate, the mine micro-seismic monitoring system is used for acquiring a plurality of natural mine seismic signals in the working face stoping process, and the natural mine seismic signals are processed to obtain the corresponding mine seismic source positions and energy thereof.
3. The microseismic statistic method for determining the working surface advanced impact risk range according to claim 1 or 2 wherein the x-axis of the new coordinate system plotted in step (2) is parallel to the working surface advancing direction and has a clockwise angle α with the x-axis of the old coordinate system.
4. The microseismic statistic method for determining the risk range of the leading impact of the working surface according to claim 3 wherein in the step (3), the coordinate transformation formula is as follows:
x′ ij =(x ij -x o )cos(α)+(y ij -y 0 )sin(α),
y′ ij =(y ij -y 0 )cos(α)-(x ij -x o )sin(α);
in the formula, x ij 、y ij Respectively representing x and y coordinates of the j microseismic event on the ith day in an old coordinate system; alpha represents a clockwise included angle between the new coordinate and the x axis of the old coordinate system; x is the number of o 、y 0 Respectively representing x and y coordinates of a new coordinate system coordinate origin point A; x' ij 、y′ ij Representing the x, y coordinates of the j microseismic event on day i in the new coordinate system.
5. The microseismic statistic method for determining the working surface advanced impact risk range according to claim 4 wherein in the step (4), the x coordinate of the working surface daily thrust position in the new coordinate system is determined by: selecting a working face lower cutting hole B point as a working face advancing starting point position, introducing working face footage data, and calculating by adopting the following formula:
x′ Bi =(x B -x 0 )cos(α)+(y B -y 0 )sin(α)±L i
in formula (II), x' Bi Representing the advancing position of the working surface on the ith day; x is the number of B 、y B Cutting down coordinates of a point B of the hole for the selected working surface; l is i And (4) indicating the cumulative ruler feeding of the day i, and selecting the + calculation when the advancing direction of the working surface is consistent with the positive direction of the x axis of the new coordinate, and otherwise selecting the-calculation.
6. The microseismic statistic method for determining the risk range of leading impact on the working surface according to claim 5 wherein in step (5), the distance of each microseismic event from the advancing position of the working surface is calculated by the following formula:
d ij =x′ Bi -x′ ij
in the formula (d) ij Indicating the j microseismic event strike distance from the face advancement position on day i.
7. The microseismic statistic method for determining the risk range of the leading impact of the working face as defined in claim 6 wherein, in the step (6), the total microseismic energy and the microseismic frequency in each divided statistic region are calculated by the following formula:
E n =e 1n +e 2n +e 3n +…+e mn
Q n =q 1n +q 2n +q 3n +…+q mn
in the formula, E n The total microseismic energy of each statistical area, namely the subarea; q n The microseismic frequency in each statistical region, namely the subarea; e.g. of the type in I =1,2,3, ·, m, for the energy of the ith microseismic event of the nth zone; q. q.s in For the nth zone, the ith microseismic event, i =1,2,3, ·, m; n is each partition number.
8. The microseismic statistical method for determining the working face advanced impact risk range according to the claim 7 wherein, in the step (8), the total microseismic energy and frequency risk level is determined according to the working face mining condition and the total microseismic energy and frequency risk standard value determined by the monitoring and early warning index.
CN202210830765.2A 2022-07-15 2022-07-15 Microseismic statistical method for determining advanced impact danger range of working face Pending CN115220092A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210830765.2A CN115220092A (en) 2022-07-15 2022-07-15 Microseismic statistical method for determining advanced impact danger range of working face

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210830765.2A CN115220092A (en) 2022-07-15 2022-07-15 Microseismic statistical method for determining advanced impact danger range of working face

Publications (1)

Publication Number Publication Date
CN115220092A true CN115220092A (en) 2022-10-21

Family

ID=83611393

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210830765.2A Pending CN115220092A (en) 2022-07-15 2022-07-15 Microseismic statistical method for determining advanced impact danger range of working face

Country Status (1)

Country Link
CN (1) CN115220092A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116794725A (en) * 2023-06-26 2023-09-22 华亭煤业集团有限责任公司 Method for correcting impact dangerous area division during stoping based on tunneling microseismic data

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116794725A (en) * 2023-06-26 2023-09-22 华亭煤业集团有限责任公司 Method for correcting impact dangerous area division during stoping based on tunneling microseismic data

Similar Documents

Publication Publication Date Title
Senfaute et al. Micro-seismic precursory cracks prior to rock-fall on coastal chalk cliffs: a case study at Mesnil-Val, Normandie, NW France
Cook Seismicity associated with mining
CN109736796B (en) Advanced detection mechanism for deepening blast hole and prediction method thereof
CN103726851B (en) Excavation method capable of lowering rockburst risk of deep tunnel in process of passing through fault area
CN111648826A (en) Coal mine rock burst prediction early warning system and method
Li et al. Characteristics of microseismic b-value associated with rock mass large deformation in underground powerhouse caverns at different stress levels
CN112377257B (en) Working face mining advance influence range determining method based on microseismic monitoring
CN110671093A (en) Comprehensive detection method for overburden rock activity rule based on drilling shooting
CN115220092A (en) Microseismic statistical method for determining advanced impact danger range of working face
Shynar et al. Methodology of assessment and prediction of critical condition of natural-technical systems
CN110761841B (en) Mine group working face mining mutual interference range calculation method based on microseismic monitoring
CN109239775B (en) Method for tracking and positioning stolen mining of mineral resources
CN112377258B (en) Method for evaluating mining influence degree between adjacent mining areas and mining working faces
CN112593914B (en) Device and method for positioning coal rock stratum fracture coordinates through cross-layer drilling
Deng et al. Fracture mechanisms of competent overburden under high stress conditions: A case study
Rasskazov et al. The researches of burst–hazard on mines in Russian Far East
Zvarivadza et al. In-stope pillar scaling and fracturing in Southern African deep level gold mines
Hudyma et al. Monitoring open stope caving at Goldex Mine
CN113434819B (en) Method for determining influence time and distance of working face mining on goaf mining vibration activities
CN111290022A (en) Rock tunnel potential seepage channel identification method based on microseism monitoring
Gangrade et al. Investigating seismicity surrounding an excavation boundary in a highly stressed dipping underground limestone mine
CN114294059A (en) High-precision control method for hydraulic fracturing layer of hard rock stratum
CN112415585A (en) Micro-seismic monitoring-based method for analyzing mining safety distance of steeply inclined extra-thick coal seam
CN112630826B (en) Method for determining mining tendency influence range of working face based on microseismic monitoring
US20230140874A1 (en) Method for monitoring hydraulic fracturing range of surface vertical shaft

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