CN110761841B - Mine group working face mining mutual interference range calculation method based on microseismic monitoring - Google Patents

Mine group working face mining mutual interference range calculation method based on microseismic monitoring Download PDF

Info

Publication number
CN110761841B
CN110761841B CN201911027646.8A CN201911027646A CN110761841B CN 110761841 B CN110761841 B CN 110761841B CN 201911027646 A CN201911027646 A CN 201911027646A CN 110761841 B CN110761841 B CN 110761841B
Authority
CN
China
Prior art keywords
microseismic
mining
rock stratum
working
microseismic event
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
Application number
CN201911027646.8A
Other languages
Chinese (zh)
Other versions
CN110761841A (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.)
China Coal Research Institute CCRI
Original Assignee
China Coal Research Institute CCRI
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 Coal Research Institute CCRI filed Critical China Coal Research Institute CCRI
Priority to CN201911027646.8A priority Critical patent/CN110761841B/en
Publication of CN110761841A publication Critical patent/CN110761841A/en
Application granted granted Critical
Publication of CN110761841B publication Critical patent/CN110761841B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21FSAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
    • E21F17/00Methods or devices for use in mines or tunnels, not covered elsewhere
    • E21F17/18Special adaptations of signalling or alarm devices
    • 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)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Business, Economics & Management (AREA)
  • Emergency Management (AREA)
  • Acoustics & Sound (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention provides a mine group working face mining mutual interference range calculation method based on microseismic monitoring, which relates to the technical field of mine monitoring and data processing, and specifically comprises the following steps: arranging microseismic monitoring between two working faces of a mine group; constructing and drilling holes in the top plate and the bottom plate, constructing and blasting in each rock stratum, and measuring the wave velocity of vibration wave propagation of each rock stratum; monitoring the two working faces in the mining process to obtain the coordinates and time of the microseismic events, and calculating and solving the maximum value of the interval between the seismic events of the same seismic source when the two working faces monitor the microseismic events generated by the same seismic source; and step four, sequentially counting the ratio of the microseismic events in each time period range, and calculating and determining the space distance of the two microseismic events in the time period group, wherein the maximum value is the mining mutual interference range of the working face. By the method, the mining mutual interference range of the two working faces of the mine group can be determined quickly and reasonably, and a basis is provided for reasonable optimal arrangement of the two working faces of the area to be mined in the mining area.

Description

Mine group working face mining mutual interference range calculation method based on microseismic monitoring
Technical Field
The invention relates to the technical field of mine monitoring and data processing, in particular to a mine group working face mining mutual interference range calculation method based on microseismic monitoring.
Background
For a particular geological condition, the overburden typically forms a caving zone, a fracture zone and a bending subsidence zone with the coal being mined, wherein the movement of the bending subsidence zone rock layer has a macroscopic control effect on the stress environment of the coal body of the working face covered by the bending subsidence zone. Practice shows that when two working faces are produced simultaneously in a local area of a mine group, the mining activity of one working face causes the high-level roof of a goaf to bend and sink, and the overall action of a thick hard rock stratum causes the stress of the coal body in the area to be redistributed, so that stronger stress disturbance is generated on the coal body of the other working face.
The construction of future coal bases is developed towards a new sustainable development mode with multiple mining areas, oversize mines, high yield, high efficiency and high benefit, with the large centralized and high-efficiency intensive mining of coal mines, the simultaneous mining of multiple working faces in local areas becomes a new normal state, the phenomenon of strong stress disturbance caused by mining is more and more obvious, and the possibility of causing coal rock dynamic disasters of the working faces by the strong disturbance is greatly increased.
At present, a direct and effective method for accurately determining the mining mutual disturbance range of two working faces does not exist, and a method is needed to be provided, wherein the mining disturbance range of the two working faces which are mined under different geological conditions, different exploitation layouts and different stoping design conditions can be calculated, so that a basis is provided for reasonable and optimal arrangement of the two working faces of a to-be-mined area of a mining area in the future, and safe mining of coal is realized.
Disclosure of Invention
In order to quickly and reasonably determine the mining mutual interference range of the two working faces of the mine group and provide basis for reasonable optimal arrangement of the two working faces of the to-be-mined area of the mining area, the invention provides a mine group working face mining mutual interference range calculation method based on microseismic monitoring, and the specific technical scheme is as follows.
A mine group working face mining mutual interference range calculation method based on microseismic monitoring specifically comprises the following steps:
step A, arranging microseismic monitoring between two working faces of a mine group;
b, drilling holes in the top plate and the bottom plate, performing construction blasting on each rock stratum, and respectively measuring wave velocity of each rock stratum vibration wave transmitted to the coal bed and the ground surface;
c, monitoring in the mining process of the two working faces to obtain coordinates and time of the microseismic event, monitoring microseismic event data generated by the same seismic source according to the two working faces, and calculating and solving the maximum value of the seismic interval of the same seismic source;
and D, taking the microseismic event data of any two working faces as a microseismic event group, sequentially counting the proportion of the microseismic event group in each time period range, and calculating and determining the space distance of any two microseismic events in the time period group to obtain the mining mutual interference range of the two working faces of the mine group.
Preferably, step a comprises:
A1. arranging microseismic monitoring systems in the mining ranges of the two working faces of the mine group, and adjusting the time service of the microseismic monitoring systems of the two working faces to be the same;
A2. arranging microseismic detectors on a coal seam roadway and the ground surface of the working face, and withdrawing and installing the microseismic detectors along with the mining advancing direction of the working face in sequence to continuously receive rock stratum fracture microseismic signals in the mining process of the working face;
A3. collecting microseismic signals generated by fracture instability of the coal bed and the overlying rock stratum in the working face mining process, processing and inverting the spatial position of the microseismic event through a microseismic system, and determining the position of the rock stratum of the microseismic event by combining a stratum comprehensive histogram of the working face position.
Preferably, in the step a, the information of the rock stratum at the occurrence position of the microseismic event is monitored and recorded, and the highest position to the lowest position of the rock stratum at which the microseismic event occurs are marked as the rock stratum 1 to the rock stratum m, and the rock stratum m to the rock stratum n respectively; wherein n is a positive integer, and m is any positive integer between 1 and n.
Further preferably, step B comprises:
B1. b, respectively constructing and drilling holes in the crossheading of the two working surfaces or in the earth surface above the working surfaces towards the positions of a top plate and a bottom plate of the coal seam, wherein the hole bottom of each constructed and drilled hole is the middle part of the rock stratum marked in the step A and subjected to the microseismic event, 200-2000 g of explosive is buried in each hole bottom, and the spatial position coordinates of the explosive are (a, b and c);
B2. determining the three-dimensional coordinates of all microseismic detectors in the coal seam in the same coordinate system and recording the three-dimensional coordinates as (d)1,e1,f1),(d2,e2,f2),···(dn,en,fn) (ii) a The three-dimensional coordinates of all microseismic detectors on the earth surface are marked as (d)1’,e1’,f1’),(d2’,e2’,f2’),···(dn’,en’,fn’);
B3. Determining the distances between the explosive space coordinate and each microseismic detector in the coal bed as follows:
Figure GDA0002897868050000021
and the distances between the explosive space coordinate and each microseismic detector on the earth surface are respectively as follows:
Figure GDA0002897868050000022
Figure GDA0002897868050000023
B4. controlling the explosive to detonate, controlling the rock stratum marked to have the microseismic event in the step A to detonate once layer by layer, and recording the detonation moment as t for the rock stratum 10Simultaneously recording the shock initiation time t of the microseismic event caused by the detonation measured by the microseismic detector in the coal seam1~tnCalculating the average wave velocity from the explosive space coordinate to each microseismic detector in the coal bed as
Figure GDA0002897868050000024
V is to be1~vnTaking the mean value
Figure GDA0002897868050000025
As the average wave velocity of the formation 1 shock wave propagating into the coal seam, wherein
Figure GDA0002897868050000026
Further preferably, step B further comprises: repeating the calculation process in the step B4 to respectively obtain the average wave velocity of the rock stratum 1 vibration wave transmitted to the surface
Figure GDA0002897868050000027
And the average wave velocity of the rock stratum 2-rock stratum n vibration waves transmitted to the coal bed
Figure GDA0002897868050000028
Average wave velocity of rock stratum 2-rock stratum n vibration waves transmitted to surface
Figure GDA0002897868050000029
It is also preferred that step C specifically comprises:
C1. marking the positions of the two working faces and all position coordinates of the microseismic detectors in the recovery process during the simultaneous recovery period of the two working faces of the mine group, and calculating the maximum space distance l between any two microseismic detectors in the two working facesmax
C2. Selecting the average wave velocity of the vibration waves from the rock stratum 1 to the rock stratum n to be transmitted to the coal bed
Figure GDA0002897868050000031
And the average wave velocity of the vibration waves from the rock stratum 1 to the rock stratum n transmitted to the surface
Figure GDA0002897868050000032
Minimum value v ofmin
C3. When two working surfaces monitor the microseismic events generated by the same seismic source, the maximum value of the microseismic event seismic time interval is calculated
Figure GDA0002897868050000033
Still preferably, step D specifically includes:
D1. recording the three-dimensional coordinates and the earthquake initiating time of the microseismic event information monitored during the simultaneous recovery of the two working surfaces, wherein the microseismic event information of one working surface is respectively recorded as (x)1,y1,z1,t1)、(x2,y2,z2,t2)、…、(xk,yk,zk,tk) And the microseismic event information of the other working surface is respectively marked as (X)1,Y1,Z1,T1)、(X2,Y2,Z2,T2)、…、(Xk,Yk,Zk,Tk);
D2. Selecting one from the microseismic events monitored by the two working surfaces to form a microseismic event group, and calculating any oneTime interval of microseismic events within a group of microseismic events, Δ T ═ Tk-tkStatistics of time periods (. DELTA.t)'max,10s]、(10s,20s]、(20s,30s]、…、((n-1)×10s,n×10s]The number of microseismic event groups in a time period of (n × 10s, (n +1) × 10s]If no microseismic event group exists, stopping counting, wherein n is an integer greater than 1;
D3. calculating the total number C of the microseismic event groups as C1+C2+…+CnCounting the number of the microseismic event groups in each time period in the step D2 to be C1、C2、…Cm…、Cn
D4. Respectively calculating the number of the microseismic event groups in the first n time periods:
Figure GDA0002897868050000034
Figure GDA0002897868050000035
when P is satisfiedm-1Less than or equal to 70 percent and PmSelecting a time period (delta t ') when the concentration is more than 70 percent'max,m×10s]All microseismic event groups are arranged in the seismic reservoir, and m is any positive integer between 1 and n; calculating the space distance l of two microseismic events in all the microseismic event groups in the selected time period, and taking the maximum value l of lmaxNamely the mining mutual interference range of the two working faces of the mine group.
The invention provides a mine group working face mining mutual disturbance range calculation method based on microseismic monitoring, wherein the disturbance range in the simultaneous mining process of two working faces of a mine group is determined by using microseismic monitoring, firstly, the propagation speed of a shock wave in a coal rock layer under the geological condition is measured by controlling blasting, then, the coordinate and time of a microseismic event monitored by a microseismic detector in the mining process of the working faces are determined according to the propagation speed, the method determines the condition which is more consistent with the actual condition, the maximum time interval of the origin time is determined according to microseismic event monitoring data with the same seismic source, and finally, the mining mutual disturbance range is determined by the statistics of a microseismic event group and the space distance of the microseismic event, so that the blank that the mining mutual disturbance range of the two working faces of the mine group cannot be determined is.
The method has the advantages of simple steps, reasonable design, flexible monitoring and calculation and strong applicability, can effectively calculate the mutual disturbance range in the mining process of the two working faces under various working conditions, can reasonably optimize the arrangement of the working faces according to the range after determining the mutual disturbance range, and provides a basis for the design of the two working faces of the to-be-mined area of the mining area in the future.
Drawings
FIG. 1 is a schematic diagram of the location of adjacent work surfaces and the placement of microseismic detectors;
FIG. 2 is a schematic diagram of microseismic event monitoring;
in the figure: 1-explosives, 2-blastholes, 3-microseismic detectors, 4-various rock formations that generate microseismic events.
Detailed Description
Referring to fig. 1 and fig. 2, a specific embodiment of a method for calculating a mining mutual disturbance range of a mine group working face based on microseismic monitoring provided by the invention is as follows.
Because a method for determining mining mutual disturbance ranges of two working faces of a mine group does not exist at present, the mutual disturbance ranges cannot be determined when the two working faces are mined simultaneously, and the working faces cannot be reasonably arranged under the condition that the mining disturbance ranges are uncertain. The invention provides a mine group working face mining mutual disturbance range calculation method based on microseismic monitoring, wherein the disturbance range in the simultaneous mining process of two working faces of a mine group is determined by using microseismic monitoring, firstly, the propagation speed of a shock wave in a coal rock layer under the geological condition is measured by controlling blasting, then, the coordinate and time of a microseismic event monitored by a microseismic detector in the mining process of the working faces are determined according to the propagation speed, the method determines the condition which is more consistent with the actual condition, the maximum time interval of the origin time is determined according to microseismic event monitoring data with the same seismic source, and finally, the mining mutual disturbance range is determined by the statistics of a microseismic event group and the space distance of the microseismic event, so that the blank that the mining mutual disturbance range of the two working faces of the mine group cannot be determined is.
A mine group working face mining mutual interference range calculation method based on microseismic monitoring specifically comprises the following steps:
and step A, arranging microseismic monitoring between two working faces of the mine group.
Specifically, the step A comprises the following steps:
A1. and micro-seismic monitoring systems are arranged in the mining ranges of the two working faces of the mine group, the time service of the micro-seismic monitoring systems of the two working faces is adjusted to be the same, and the accuracy of monitoring and recording at each moment is ensured.
A2. Microseismic detectors are arranged on a coal seam roadway and the ground surface of the working face, and are sequentially withdrawn and mounted along with the mining advancing direction of the working face, so that a rock stratum fracture microseismic signal in the mining process of the working face is continuously received.
A3. Collecting microseismic signals generated by fracture instability of the coal bed and the overlying rock stratum in the working face mining process, processing and inverting the spatial position of the microseismic event through a microseismic system, and determining the position of the rock stratum of the microseismic event by combining a stratum comprehensive histogram of the working face position.
Specifically, the method comprises the steps of monitoring vibration signals generated by fracture instability of a coal bed and an overlying rock stratum in the mining process of two working faces in real time through a microseismic sensor, performing space positions of microseismic events by using microseismic system processing software in an inverted mode, importing comprehensive columnar information of strata near the two working faces into software, and displaying the positions of the rock strata where the microseismic events occur by using the software.
In addition, monitoring and recording the rock stratum information of the occurrence position of the microseismic event in the step A, and marking the highest layer position to the lowest layer position of the rock stratum with the microseismic event as a rock stratum 1 to a rock stratum m and marking the rock stratum m to a rock stratum n (rock stratum 1 to n); wherein n is a positive integer, and m is any positive integer between 1 and n.
And B, drilling holes in the top plate and the bottom plate, performing construction blasting on each rock stratum, and respectively measuring wave velocity of each rock stratum vibration wave transmitted to the coal bed and the earth surface.
Specifically, the step B comprises the following steps:
B1. and B, respectively constructing and drilling holes in the crossheading of the two working surfaces or in the earth surface above the working surfaces towards the positions of the top plate and the bottom plate of the coal seam, wherein the hole bottom of each constructed and drilled hole is the middle part of the rock stratum marked in the step A and subjected to the microseismic event, 200-2000 g of explosive is buried at each hole bottom, and the spatial position coordinates of the explosive are (a, b and c).
B2. Determining the three-dimensional coordinates of all microseismic detectors in the coal seam in the same coordinate system and recording the three-dimensional coordinates as (d)1,e1,f1),(d2,e2,f2),···(dn,en,fn) (ii) a The three-dimensional coordinates of all microseismic detectors on the earth surface are marked as (d)1’,e1’,f1’),(d2’,e2’,f2’),···(dn’,en’,fn’)。
B3. Determining the distances between the explosive space coordinate and each microseismic detector in the coal bed as follows:
Figure GDA0002897868050000051
and the distances between the explosive space coordinate and each microseismic detector on the earth surface are respectively as follows:
Figure GDA0002897868050000052
Figure GDA0002897868050000053
B4. artificially controlling the explosive to detonate, controlling the rock stratum marked to have the microseismic event in the step A to detonate once layer by layer, and recording the detonation moment t as t for the rock stratum 10Simultaneously recording the shock initiation time t of the microseismic event caused by the detonation measured by the microseismic detector in the coal seam1~tnCalculating the average wave velocity from the explosive space coordinate to each microseismic detector in the coal bed as
Figure GDA0002897868050000054
Figure GDA0002897868050000055
V is to be1~vnTaking the mean value
Figure GDA0002897868050000056
As the average wave velocity of the formation 1 shock wave propagating into the coal seam, wherein
Figure GDA0002897868050000057
In addition, step B further comprises: repeating the calculation process in the step B4 to respectively obtain the average wave velocity of the rock stratum 1 vibration wave transmitted to the surface
Figure GDA0002897868050000058
And the average wave velocity of the rock stratum 2-rock stratum n vibration waves transmitted to the coal bed
Figure GDA0002897868050000059
Average wave velocity of rock stratum 2-rock stratum n vibration waves transmitted to surface
Figure GDA00028978680500000510
And C, monitoring in the mining process of the two working faces to obtain the coordinates and time of the microseismic event, monitoring microseismic event data generated by the same seismic source according to the two working faces, and calculating and solving the maximum value of the seismic interval of the same seismic source.
Specifically, the step C includes:
C1. marking the positions of the two working faces and all position coordinates of the microseismic detectors in the recovery process during the simultaneous recovery period of the two working faces of the mine group, and calculating the maximum space distance l between any two microseismic detectors in the two working facesmax
C2. Selecting the average wave velocity of the vibration waves from the rock stratum 1 to the rock stratum n to be transmitted to the coal bed
Figure GDA00028978680500000511
And the average wave velocity of the vibration waves from the rock stratum 1 to the rock stratum n transmitted to the surface
Figure GDA0002897868050000061
Minimum value v ofmin
C3. When two working surfaces monitor the microseismic events generated by the same seismic source, the maximum value of the microseismic event seismic time interval is calculated
Figure GDA0002897868050000062
And D, taking the microseismic event data of any two working faces as a microseismic event group, sequentially counting the proportion of the microseismic event group in each time period range, and calculating and determining the space distance of any two microseismic events in the time period group to obtain the mining mutual interference range of the two working faces of the mine group.
Specifically, the step D includes:
D1. recording the three-dimensional coordinates and the earthquake initiating time of the microseismic event information monitored during the simultaneous recovery of the two working surfaces, wherein the microseismic event information of one working surface is respectively recorded as (x)1,y1,z1,t1)、(x2,y2,z2,t2)、…、(xk,yk,zk,tk) And the microseismic event information of the other working surface is respectively marked as (X)1,Y1,Z1,T1)、(X2,Y2,Z2,T2)、…、(Xk,Yk,Zk,Tk)。
D2. Selecting one from the microseismic events monitored by the two working surfaces to form a microseismic event group, and calculating the time interval delta T of the microseismic events in any group of microseismic event groups as Tk-tkStatistics of time periods (. DELTA.t)'max,10s]、(10s,20s]、(20s,30s]、…、((n-1)×10s,n×10s]The number of microseismic event groups in a time period of (n × 10s, (n +1) × 10s]If no microseismic event group exists, stopping counting, wherein n is an integer greater than 1.
D3. Calculating the total number C of the microseismic event groups as C1+C2+…+CnCounting the number of the microseismic event groups in each time period in the step D2 to be C1、C2、…Cm…、Cn
D4. Respectively calculating the number of the microseismic event groups in the first n time periods:
Figure GDA0002897868050000063
Figure GDA0002897868050000064
when P is satisfiedm-1Less than or equal to 70 percent and PmSelecting a time period (delta t ') when the concentration is more than 70 percent'max,m×10s]All microseismic event groups are arranged in the seismic reservoir, and m is any positive integer between 1 and n; calculating the space distance l of two microseismic events in all the microseismic event groups in the selected time period, and taking the maximum value l of lmaxNamely the mining mutual interference range of the two working faces of the mine group.
The method has the advantages of simple steps, reasonable design, flexible monitoring and calculation and strong applicability, can effectively calculate the mutual disturbance range in the mining process of the two working faces under various working conditions, can reasonably optimize the arrangement of the working faces according to the range after determining the mutual disturbance range, and provides a basis for the design of the two working faces of the same coal seam in the to-be-mined area of the mining area in the future.
It is to be understood that the above description is not intended to limit the present invention, and the present invention is not limited to the above examples, and those skilled in the art may make modifications, alterations, additions or substitutions within the spirit and scope of the present invention.

Claims (5)

1. A mine group working face mining mutual disturbance range calculation method based on microseismic monitoring is characterized by specifically comprising the following steps:
step A, arranging microseismic monitoring between two working faces of a mine group;
b, drilling holes in the top plate and the bottom plate, performing construction blasting on each rock stratum, and respectively measuring wave velocity of each rock stratum vibration wave transmitted to the coal bed and the ground surface;
c, monitoring in the mining process of the two working faces to obtain coordinates and time of the microseismic event, monitoring microseismic event data generated by the same seismic source according to the two working faces, and calculating and solving the maximum value of the seismic interval of the same seismic source;
the method specifically comprises the following steps:
C1. marking the positions of the two working faces and all position coordinates of the microseismic detectors in the recovery process during the simultaneous recovery period of the two working faces of the mine group, and calculating the maximum space distance l between any two microseismic detectors in the two working facesmax
C2. Selecting rock stratum 1 to rock stratumAverage wave velocity of n vibration waves transmitted to coal seam
Figure FDA0002981996450000011
And the average wave velocity of the vibration waves from the rock stratum 1 to the rock stratum n transmitted to the surface
Figure FDA0002981996450000012
Minimum value v ofmin
C3. When two working surfaces monitor the microseismic events generated by the same seismic source, the maximum value of the microseismic event seismic time interval is calculated
Figure FDA0002981996450000013
D, taking the microseismic event data of any two working faces as a microseismic event group, sequentially counting the proportion of the microseismic event group in each time period range, and calculating and determining the space distance of any two microseismic events in the time period group to obtain the mining mutual interference range of the two working faces of the mine group;
the method specifically comprises the following steps:
D1. recording the three-dimensional coordinates and the earthquake initiating time of the microseismic event information monitored during the simultaneous recovery of the two working surfaces, wherein the microseismic event information of one working surface is respectively recorded as (x)1,y1,z1,t1)、(x2,y2,z2,t2)、…、(xk,yk,zk,tk) And the microseismic event information of the other working surface is respectively marked as (X)1,Y1,Z1,T1)、(X2,Y2,Z2,T2)、…、(Xk,Yk,Zk,Tk);
D2. Selecting one from the microseismic events monitored by the two working surfaces to form a microseismic event group, and calculating the time interval delta T of the microseismic events in any group of microseismic event groups as Tk-tkCounting each time period (Δ t)max,10s]、(10s,20s]、(20s,30s]、…、((n-1)×10s,n×10s]The number of microseismic event groups in a time period of (n × 10s, (n +1) × 10s]No micro-shock in the interiorIf the group is grouped, stopping counting, wherein n is an integer greater than 1;
D3. calculating the total number C of the microseismic event groups as C1+C2+…+CnCounting the number of the microseismic event groups in each time period in the step D2 to be C1、C2、…Cm…、Cn
D4. Respectively calculating the number of the microseismic event groups in the first n time periods:
Figure FDA0002981996450000014
Figure FDA0002981996450000021
when P is satisfiedm-1Less than or equal to 70 percent and PmWhen the time is more than 70 percent, selecting a time period (delta t)max,m×10s]All microseismic event groups are arranged in the seismic reservoir, and m is any positive integer between 1 and n; and calculating the spatial distance l of the two microseismic events in all the microseismic event groups in the selected time period, and taking the maximum value of l as the mining mutual interference range of the two working faces of the mine group.
2. The method for calculating the mining mutual disturbance range of the mine group working face based on the microseismic monitoring as claimed in claim 1, wherein the step A comprises the following steps:
A1. arranging microseismic monitoring systems in the mining ranges of the two working faces of the mine group, and adjusting the time service of the microseismic monitoring systems of the two working faces to be the same;
A2. arranging microseismic detectors on a coal seam roadway and the ground surface of the working face, and withdrawing and installing the microseismic detectors along with the mining advancing direction of the working face in sequence to continuously receive rock stratum fracture microseismic signals in the mining process of the working face;
A3. collecting microseismic signals generated by fracture instability of the coal bed and the overlying rock stratum in the working face mining process, processing and inverting the spatial position of the microseismic event through a microseismic system, and determining the position of the rock stratum of the microseismic event by combining a stratum comprehensive histogram of the working face position.
3. The method for calculating the mining mutual interference range of the mine group working face based on the microseismic monitoring as claimed in claim 2, wherein in the step A, the information of rock strata at the position where the microseismic event occurs is monitored and recorded, the rock strata where the microseismic event occurs is marked, and the positions from the highest position to the lowest position of the rock strata are rock strata 1 to rock stratum n; wherein n is a positive integer.
4. The method for calculating the mining mutual disturbance range of the mine group working face based on the microseismic monitoring as claimed in claim 1 or 3, wherein the step B comprises the following steps:
B1. b, respectively constructing and drilling holes in the crossheading of the two working surfaces or in the earth surface above the working surfaces towards the positions of a top plate and a bottom plate of the coal seam, wherein the hole bottom of each constructed and drilled hole is the middle part of the rock stratum marked in the step A and subjected to the microseismic event, 200-2000 g of explosive is buried in each hole bottom, and the spatial position coordinates of the explosive are (a, b and c);
B2. determining the three-dimensional coordinates of all microseismic detectors in the coal seam in the same coordinate system and recording the three-dimensional coordinates as (d)1,e1,f1),(d2,e2,f2),···(dn,en,fn) (ii) a The three-dimensional coordinates of all microseismic detectors on the earth surface are marked as (d)1’,e1’,f1’),(d2’,e2’,f2’),···(dn’,en’,fn’);
B3. Determining the distances between the explosive space coordinate and each microseismic detector in the coal bed as follows:
Figure FDA0002981996450000022
and the distances between the explosive space coordinate and each microseismic detector on the earth surface are respectively as follows:
Figure FDA0002981996450000023
Figure FDA0002981996450000024
B4. controlling the explosive to detonate, controlling the detonation once layer by layer for the rock stratum marked to have the microseismic event in the step A, and controlling the detonation once for the rock stratum marked to have the microseismic event in the step AFormation 1, recording the detonation time t0Simultaneously recording the shock initiation time t of the microseismic event caused by the detonation measured by the microseismic detector in the coal seam1~tnCalculating the average wave velocity from the explosive space coordinate to each microseismic detector in the coal bed as
Figure FDA0002981996450000025
V is to be1~vnTaking the mean value
Figure FDA0002981996450000031
As the average wave velocity of the formation 1 shock wave propagating into the coal seam, wherein
Figure FDA0002981996450000032
5. The method for calculating the mining mutual disturbance range of the mine group working face based on the microseismic monitoring as claimed in claim 4, wherein the step B further comprises the following steps: repeating the calculation process in the step B4 to respectively obtain the average wave velocity of the rock stratum 1 vibration wave transmitted to the surface
Figure FDA0002981996450000033
And the average wave velocity of the rock stratum 2-rock stratum n vibration waves transmitted to the coal bed
Figure FDA0002981996450000034
Average wave velocity of rock stratum 2-rock stratum n vibration waves transmitted to surface
Figure FDA0002981996450000035
CN201911027646.8A 2019-10-28 2019-10-28 Mine group working face mining mutual interference range calculation method based on microseismic monitoring Active CN110761841B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911027646.8A CN110761841B (en) 2019-10-28 2019-10-28 Mine group working face mining mutual interference range calculation method based on microseismic monitoring

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911027646.8A CN110761841B (en) 2019-10-28 2019-10-28 Mine group working face mining mutual interference range calculation method based on microseismic monitoring

Publications (2)

Publication Number Publication Date
CN110761841A CN110761841A (en) 2020-02-07
CN110761841B true CN110761841B (en) 2021-06-04

Family

ID=69333943

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911027646.8A Active CN110761841B (en) 2019-10-28 2019-10-28 Mine group working face mining mutual interference range calculation method based on microseismic monitoring

Country Status (1)

Country Link
CN (1) CN110761841B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112464143B (en) * 2020-10-23 2024-07-26 中国石油天然气集团有限公司 Method and device for identifying in-situ gasification boundary of underground coal
CN113496006B (en) * 2021-08-02 2023-08-01 中国矿业大学 Method for calculating centralization degree of mine earthquake space of rock burst mine
CN114185081B (en) * 2021-11-11 2023-10-27 北京科技大学 Ground vertical well hydraulic fracturing range monitoring method

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005006020A1 (en) * 2003-07-05 2005-01-20 Westerngeco Seismic Holdings Limited Passive seismic event detection
CN202330733U (en) * 2011-11-21 2012-07-11 大同煤矿集团有限责任公司 Micro-seismic monitoring system for working surface for thick seam full-mechanized caving mining
CN103697999A (en) * 2013-12-30 2014-04-02 中国科学院武汉岩土力学研究所 Method for acquiring micro seismic wave speed of high-stress hard rock TBM (Tunnel Boring Machine) construction tunnel in real time
CN104018790A (en) * 2014-06-04 2014-09-03 天地科技股份有限公司 Roadway rock burst early-warning method based on rock noise monitoring
CN104656123A (en) * 2015-01-23 2015-05-27 四川大学 Method for measuring equivalent wave velocity of regional rock mass
CN105467436A (en) * 2015-12-17 2016-04-06 东北大学 Method for arranging micro-shock sensors in construction of ultra-deep vertical shaft
CN108930554A (en) * 2018-08-20 2018-12-04 煤炭科学技术研究院有限公司 Well-ground-air combined monitoring method of the fracture of coal mine overlying strata and surface strata movement
CN110159347A (en) * 2019-05-05 2019-08-23 北京科技大学 A kind of deep high stress tight roof stope dynamic disaster monitoring and pre-alarming method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10054702B2 (en) * 2014-10-24 2018-08-21 Schlumberger Technology Corporation Method to enhance the resolvability of moment tensor inversion for III conditioned receiver coverage

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005006020A1 (en) * 2003-07-05 2005-01-20 Westerngeco Seismic Holdings Limited Passive seismic event detection
CN202330733U (en) * 2011-11-21 2012-07-11 大同煤矿集团有限责任公司 Micro-seismic monitoring system for working surface for thick seam full-mechanized caving mining
CN103697999A (en) * 2013-12-30 2014-04-02 中国科学院武汉岩土力学研究所 Method for acquiring micro seismic wave speed of high-stress hard rock TBM (Tunnel Boring Machine) construction tunnel in real time
CN104018790A (en) * 2014-06-04 2014-09-03 天地科技股份有限公司 Roadway rock burst early-warning method based on rock noise monitoring
CN104656123A (en) * 2015-01-23 2015-05-27 四川大学 Method for measuring equivalent wave velocity of regional rock mass
CN105467436A (en) * 2015-12-17 2016-04-06 东北大学 Method for arranging micro-shock sensors in construction of ultra-deep vertical shaft
CN108930554A (en) * 2018-08-20 2018-12-04 煤炭科学技术研究院有限公司 Well-ground-air combined monitoring method of the fracture of coal mine overlying strata and surface strata movement
CN110159347A (en) * 2019-05-05 2019-08-23 北京科技大学 A kind of deep high stress tight roof stope dynamic disaster monitoring and pre-alarming method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《矿井群冲击地压发生机理与控制技术探讨》;齐庆新等;《煤炭学报》;20190630;第44卷(第1期);第141-150页 *

Also Published As

Publication number Publication date
CN110761841A (en) 2020-02-07

Similar Documents

Publication Publication Date Title
Konicek et al. Stress changes and seismicity monitoring of hard coal longwall mining in high rockburst risk areas
Ghosh et al. Application of underground microseismic monitoring for ground failure and secure longwall coal mining operation: a case study in an Indian mine
CN110761841B (en) Mine group working face mining mutual interference range calculation method based on microseismic monitoring
Zhao et al. The analysis of rock damage process based on the microseismic monitoring and numerical simulations
CN112965136B (en) Multi-means advanced detection method for water-rich karst tunnel
CN110703320A (en) Up-down combined microseismic monitoring system and method
Fan et al. Advanced stability analysis of the tunnels in jointed rock mass based on TSP and DEM
CN110671093A (en) Comprehensive detection method for overburden rock activity rule based on drilling shooting
Hanson et al. Advanced techniques in site characterization and mining hazard detection for the underground coal industry
CN116049964B (en) Anti-impact method for artificial liberation layer of newly-built mine
Eremenko et al. Effect of mining geometry on natural stress field in underground ore mining with conventional and nature-like technologies
Vinoth et al. Applying real time seismic monitoring technology for slope stability assessment—an Indian opencast coal mine perspective
Potvin et al. Seismic monitoring in highly mechanized hardrock mines in Canada and Australia
Wang et al. Microseismicity evolution related to two extremely intense rockbursts in a water diversion tunnel
Dahnér et al. Transition from non-seismic mine to a seismically active mine: Kiirunavaara Mine
CN112965139B (en) Advanced geological comprehensive forecasting method for tunnel with complex geological condition
Rasskazov et al. The researches of burst–hazard on mines in Russian Far East
CN115045660B (en) Deep hole blasting fracturing method and system for top plate of large-dip angle stope face
Hudyma et al. Monitoring open stope caving at Goldex Mine
Mutke et al. Seismic monitoring and rock burst hazard assessment in Deep Polish Coal Mines–Case study of rock burst on April 16, 2008 in Wujek-Slask Coal Mine
Westley-Hauta et al. Characterisation of seismic activity at a kimberlite block caving operation in a complex geological setting in Quebec, Canada
Vinoth et al. Slope stability monitoring by quantification and behavior of microseismic events in an opencast coal mine
Zubkov et al. Rock pressure control methods based on detected regularities of stress formation in mining structures
Holub et al. Some parameters of rockbursts derived from underground seismological measurements
Wang et al. A polarization migration velocity model building method for geological prediction ahead of the tunnel face

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