CN111963243B - Rock burst danger monitoring and early warning method based on dynamic and static combined stress analysis - Google Patents
Rock burst danger monitoring and early warning method based on dynamic and static combined stress analysis Download PDFInfo
- Publication number
- CN111963243B CN111963243B CN202010714159.5A CN202010714159A CN111963243B CN 111963243 B CN111963243 B CN 111963243B CN 202010714159 A CN202010714159 A CN 202010714159A CN 111963243 B CN111963243 B CN 111963243B
- Authority
- CN
- China
- Prior art keywords
- dynamic
- static
- stress
- rock burst
- microseismic
- 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
Links
- 230000003068 static effect Effects 0.000 title claims abstract description 113
- 239000011435 rock Substances 0.000 title claims abstract description 75
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000012544 monitoring process Methods 0.000 title claims abstract description 41
- 238000004458 analytical method Methods 0.000 title claims abstract description 26
- 230000000694 effects Effects 0.000 claims abstract description 72
- 230000006378 damage Effects 0.000 claims abstract description 49
- 238000009825 accumulation Methods 0.000 claims abstract description 31
- 230000007246 mechanism Effects 0.000 claims abstract description 11
- 230000002301 combined effect Effects 0.000 claims abstract description 8
- 230000035939 shock Effects 0.000 claims description 17
- 150000001875 compounds Chemical class 0.000 claims description 15
- 239000000523 sample Substances 0.000 claims description 15
- 230000001186 cumulative effect Effects 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 238000013016 damping Methods 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims 1
- 230000008901 benefit Effects 0.000 abstract description 4
- 238000004445 quantitative analysis Methods 0.000 abstract description 4
- 239000003245 coal Substances 0.000 description 45
- 230000009471 action Effects 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 238000005065 mining Methods 0.000 description 7
- 230000001939 inductive effect Effects 0.000 description 5
- 230000002706 hydrostatic effect Effects 0.000 description 4
- 238000009412 basement excavation Methods 0.000 description 3
- 238000005422 blasting Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000006698 induction Effects 0.000 description 3
- 230000001052 transient effect Effects 0.000 description 3
- 230000005670 electromagnetic radiation Effects 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
- E21F17/18—Special adaptations of signalling or alarm devices
Landscapes
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Geophysics And Detection Of Objects (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention discloses a rock burst danger monitoring and early warning method based on dynamic and static combined stress analysis, which carries out rock burst danger monitoring and early warning based on a dynamic load effect and static load effect superposition-induced rock burst mechanism and combined dynamic and static combined effect; the dynamic and static combined effect is characterized in that dynamic load effect parameters reconstructed by vibration wave attenuation accumulation, static load effect parameters reconstructed by vibration wave CT and the weight of each parameter of static load effect parameters reconstructed by micro-seismic damage are determined respectively, dynamic and static combined stress coefficients are calculated by adopting a weighted average method, and a final danger area and the danger degree thereof are determined comprehensively by combining the dynamic load stress coefficients of mine pressure instant dynamic load stress increment. The method has the advantages of definite calculation model, universality and strong operability, can realize quantitative analysis of the rock burst dangerous area and the prediction of the dangerous degree thereof, and has good application feasibility.
Description
Technical Field
The invention relates to a rock burst hazard monitoring and early warning method based on dynamic and static combined stress analysis.
Background
Rock burst is a phenomenon that burst type damage occurs in an adjacent empty rock mass in a deep part or a region with high structural stress of underground mining, and is a dynamic phenomenon that energy accumulated in a roadway and a coal rock mass of a stope is suddenly released under disturbance of mine earthquake dynamic load, the coal rock mass is thrown to the roadway, and strong sound is generated at the same time, so that vibration and damage of the coal rock mass, damage of a support and equipment, casualties, partial roadway landing damage and the like are caused.
According to long-term theoretical research, laboratory tests and field tests, the rock burst generation is consistently considered to meet the dynamic and static load superposition induction principle, and is particularly shown in the coal cutting, frame moving, mechanical vibration, blasting, roof and floor breakage, coal body and roof structure instability, gas outburst, coal cannon, fault slippage and other mine earthquake dynamic loads and the bearing stress (static load) of coal bodies around a stope and a roadway, wherein the coal cutting, frame moving, mechanical vibration, blasting, roof and floor breakage, the coal body and roof structure instability, the gas outburst, the coal cannon, the fault slippage and the like are shown, and once the dynamic and static load superposition exceeds the bearing limit of the coal bodies, the rock burst is easily induced, so that certain damage is caused.
Due to the characteristics of complexity, burstiness, diversity and the like of rock burst, monitoring and early warning of the rock burst as a multi-dimensional spatial information description problem cannot be expressed by a single parameter, and is more remarkable particularly when nonlinear physical change of a coal rock catastrophe fracture process is considered, for example, in the prior art:
chinese patent application No. CN201811357374.3 discloses a method for evaluating risk of rock burst in mining area, comprising the following steps: s001: collecting geological data of a mining area; s002: predetermining the danger level of rock burst and the danger index range of the rock burst; s003: determining parameters, weights and evaluation index value ranges of coal seam mining influence factors of the rock burst; s004: calculating a rock burst danger comprehensive index Y total of the influence factors on the rock burst danger; s005: and comparing the Y total with the rock burst risk index range in the step S002, and determining the rock burst risk grade corresponding to the Y total according to the comparison result.
Chinese patent application No. CN201210155302.7 discloses a rock burst prediction and early warning method, which is used for monitoring mine earthquake in the whole coal field, positioning the earthquake magnitude, and analyzing the mine earthquake occurrence and motion rules in the whole mining area; installing a Polish microseismic system on the upper and lower mountains of a mining area and a transportation main roadway, and implementing online continuous monitoring; installing a microseismic system probe on a working face, installing a microseismic monitoring probe in a key scour prevention area of the working face, implementing deployment and control, drilling a hole in a coal wall or a top plate, placing the probe into the bottom and the hole wall of the hole, utilizing an electromagnetic radiation instrument to monitor electromagnetic radiation of a return airway and a transportation lane of the working face, arranging a monitoring probe in the key scour prevention area for deployment and control, monitoring mine pressure, inputting parameters obtained by monitoring into a computer for comprehensive analysis, and realizing accurate prediction of rock burst.
However, in the prior art, the calculation of the comprehensive indexes of the impact ground pressure dangerousness of the impact ground pressure caused by the influence factors or the input of the monitored parameters into the computer for comprehensive analysis cannot realize the quantitative description of the impact ground pressure dangerous area and the danger degree thereof, so as to achieve the purpose of accurate early warning, and how to utilize the existing monitoring means to specifically quantify and determine the dynamic and static load effects and the superposition thereof in the impact ground pressure induction mechanism is the key point for realizing the accurate early warning of the impact ground pressure and the quantification of the result thereof; in addition, technicians in coal mines are difficult to effectively mine essential and useful information sources by applying scientific means in the face of multi-source massive monitoring data. Therefore, a rock burst hazard monitoring and early warning method based on dynamic and static combined stress analysis is urgently needed to be researched.
Disclosure of Invention
Aiming at the problems, the invention provides a rock burst danger monitoring and early warning method based on dynamic and static combined stress analysis, which realizes comprehensive monitoring and early warning of rock burst danger by quantitatively analyzing dynamic and static combined stress in the underground coal seam excavation process; furthermore, the method has the advantages of clear calculation model, universality and strong operability, can realize quantitative analysis of the rock burst dangerous area and the prediction of the dangerous degree thereof, and has good application feasibility.
In order to achieve the technical purpose and achieve the technical effect, the invention is realized by the following technical scheme:
a rock burst danger monitoring and early warning method based on dynamic and static combined stress analysis is based on a dynamic load effect and static load effect superposition induced rock burst mechanism and combines a dynamic and static combined effect to carry out rock burst danger monitoring and early warning, wherein:
the dynamic load effect is quantitatively described through dynamic load effect parameters reconstructed by vibration wave attenuation accumulation and dynamic load effect parameters of mine pressure instant dynamic load stress increment;
the static load effect is quantitatively described through a static load effect parameter inverted by the vibration wave CT and a static load effect parameter reconstructed by the microseismic damage;
the dynamic and static combined effect is characterized in that dynamic load effect parameters reconstructed by vibration wave attenuation accumulation, static load effect parameters reconstructed by vibration wave CT and the weight of each parameter of static load effect parameters reconstructed by micro-seismic damage are determined respectively, dynamic and static combined stress coefficients are calculated by adopting a weighted average method, and a final danger area and the danger degree thereof are determined comprehensively by combining the dynamic load stress coefficients of mine pressure instant dynamic load stress increment.
Preferably, the dynamic load effect parameter of the seismic wave attenuation accumulation reconstruction is calculated by adopting the following formula:
in the formula (I), the compound is shown in the specification,
dynamic load stress coefficient, σ, reconstructed for seismic wave attenuation accumulationsdsDynamic load stress parameter, sigma, reconstructed for damping accumulation of shock wavessdsaAverage value, sigma, of dynamic load stress parameter values reconstructed for damping accumulation of shock wavessdsiIs σ at the ith spatial position nodesds，EjIs the jth microseismic event energy, RijIs the distance between the jth microseismic event and the ith spatial location node, beta is the attenuation coefficient of the shock wave in the formation, NMSThe microseismic frequency.
Preferably, the dynamic load effect parameter of the transient dynamic load stress increment of the mine pressure is calculated by adopting the following formula:
in the formula (I), the compound is shown in the specification,
dynamic load stress coefficient, P, being the instantaneous dynamic load stress increment of the mine pressureTiThe weighted average of the stent's operational resistance in the ith time window,
is PTiAverage value of (a) ("sigmaPIs PTiStandard deviation of (1), tpj、pjRespectively corresponding time and pressure values of the real-time monitoring sequence of the working resistance of the bracket, and n is the number of monitoring data points in a time window.
Preferably, the static load effect parameter of the seismic wave CT inversion is calculated by using the following formula:
in the formula (I), the compound is shown in the specification,
static load stress coefficient, v, obtained for seismic wave CT inversionPLongitudinal wave velocity, v, obtained for seismic wave CT inversionPaAnd alpha is an experimental test fitting parameter.
Preferably, the static load effect parameter of the microseismic damage reconstruction is calculated by adopting the following formula:
in the formula (I), the compound is shown in the specification,
coefficient of static load stress, σ, for microseismic damage reconstructionsssStatic load stress parameter, sigma, for microseismic damage reconstructionsssaAverage of static load stress parameter values reconstructed for microseismic damage, NMSFor microseismic frequency, the exp function is an exponential function with a natural constant e as the base,
to accumulate Benioff strain, E is the microseismic event energy, max { εEIs the maximum value of the cumulative Benioff strain, DFIs a critical damage value, corresponding to a complete damage state, εFIs the critical cumulative Benioff strain.
Preferably, the standard deviation method is adopted to determine and) the weights of the dynamic load effect parameter reconstructed by the attenuation accumulation of the vibration wave, the static load effect parameter reconstructed by the CT inversion of the vibration wave and the static load effect parameter reconstructed by the microseismic damage, and the calculation formula is as follows:
in the formula (I), the compound is shown in the specification,
the static load stress coefficient obtained for the inversion of the vibration wave CT is weighted by wvss；
The coefficient of static load stress for reconstructing microseismic damage is weighted by wsss；
The dynamic load stress coefficient reconstructed for the attenuation accumulation of the vibration wave is weighted by wsds(ii) a The subscript StDv indicates the standard deviation of each parameter.
Preferably, dynamic and static combined stress coefficient
Calculated using the following formula:
in the formula, subscripts min and max represent the minimum value and the maximum value of each parameter, respectively.
Preferably, the attenuation coefficient beta of the shock wave in the stratum is obtained by fitting and calculating the attenuation change relation of the speed amplitude of the micro-seismic probe with the distance between the micro-seismic probe and the seismic source, wherein the speed amplitude of the micro-seismic probe is not beyond the measuring range in the far field.
Preferably, max { εEThe position is regarded as critical damage state, DFIs 0.95.
Preferably, when the value of the dynamic and static combined stress coefficient in a certain area is larger than a set first threshold value, or the value of the dynamic load stress coefficient of the dynamic load stress increment at the moment of mine pressure is larger than a set second threshold value, the system immediately gives an early warning prompt.
The invention has the beneficial effects that:
the dynamic load effect is quantitatively depicted by two parameters of shock wave attenuation accumulation reconstruction and mine pressure instant dynamic load stress increment, the static load effect is quantitatively described by two parameters of shock wave CT inversion and microseismic damage reconstruction, the four parameters of the static load effect of the shock wave CT inversion, the static load effect of the microseismic damage reconstruction, the dynamic load effect of the shock wave attenuation accumulation reconstruction and the dynamic load effect of the mine pressure instant dynamic load stress increment have clear physical meanings, the dynamic and static and superposition effects of an impact pressure induction mechanism are comprehensively considered, a calculation model is clear, universality and operability are strong, quantitative analysis of an impact pressure dangerous area and the danger degree prediction can be realized, the application feasibility is good, and accurate prediction can be realized. In addition, the rapid updating and adjusting of the weight and the objective judgment of the final comprehensive prediction result have the advantage of high prediction efficiency.
Drawings
FIG. 1 is a schematic diagram of the mechanism of the present invention for inducing the dynamic and static load of impact ground pressure;
FIG. 2 is a schematic view of dynamic and static combined stress analysis during rock burst generation according to the present invention;
FIG. 3 is an impact case and monitoring data of the present invention;
FIG. 4 is a static load stress coefficient result graph obtained by the inversion of the seismic wave CT of the present invention;
FIG. 5 is a static load stress coefficient result plot for microseismic damage reconstruction in accordance with the present invention;
FIG. 6 is a graph of the peak velocity propagation attenuation characteristics of particles in accordance with the present invention;
FIG. 7 is a graph of the dynamic load stress coefficient results of the seismic wave attenuation accumulation reconstruction of the present invention;
FIG. 8 is a dynamic load stress coefficient result plot of the transient dynamic load stress increment of mine pressure of the present invention;
FIG. 9 is a dynamic and static combined stress coefficient result diagram for rock burst hazard prediction in accordance with the present invention.
Detailed Description
The present invention will be better understood and implemented by those skilled in the art by the following detailed description of the technical solution of the present invention with reference to the accompanying drawings and specific examples, which are not intended to limit the present invention.
The underground coal body excavation causes stress field disturbance, and a stress concentration area which is a supporting stress area and exceeds the original rock stress is formed in the coal body transversely in front of the working face and around the roadway. Meanwhile, excavation causes damage and migration of overlying strata, as shown in fig. 1, a collapse zone, a fissure zone and a bending subsidence zone are formed in the longitudinal direction and respectively correspond to a post-peak damage zone DE, a pre-peak plastic zone BD and an elastic zone AB in a transverse bearing stress zone.
Under the high ground stress condition in deep, the former rock is in accurate hydrostatic pressure state, therefore the stress environment change of the coal rock body in front of the deep working face begins in accurate hydrostatic pressure state, along with the promotion of working face, the bearing pressure (vertical stress) in the coal seam rises gradually to the peak value stress from the hydrostatic pressure state of three-dimensional isobaric, then goes into the release state along with the destruction of coal body, and vertical stress reduces gradually until the single pressure residual strength state of coal wall department. On the other hand, the horizontal stress is gradually reduced from the three-way isobaric hydrostatic pressure state to 0, i.e. pressure relief. As can be seen from fig. 1, the dynamic and static load superposition inducing mechanism of the rock burst can be represented by the following formula, wherein the dynamic and static load superposition inducing mechanism of the rock burst is represented by the following formula, namely, the coal cutting, the frame moving, the mechanical vibration, the blasting, the top and bottom plate breaking, the coal body and top plate structure instability, the gas outburst, the coal cannon, the fault slippage and other mine vibration loads are superposed with the supporting stress (static load) of the coal body around the stope and the roadway, and once the bearing limit of the coal body is exceeded, the dynamic and static load superposition inducing mechanism of the rock burst can be represented by the following formula:
σs+σd≥σbmin (1)
in the formula, σsIs the static load stress, sigma, of the coal-rock massdFor mine vibration loading stress, sigmabminIs the critical stress of rock burst.
According to the rock burst generation model shown in fig. 1, a roof-coal seam-floor impact carrier system model shown in fig. 2 and a dynamic and static combined stress analysis schematic diagram thereof can be established. That is, as the coal seam mining activity advances, the stress-strain relationship of the coal body under load in front of the face may be described by the right curve in fig. 2. Considering the roof and floor as a complete wall rock and its rigidity and strength are much greater than those of the coal seam, the stress-strain relationship under load can be represented by the left curve in fig. 2. In fig. 2, U1 is the elastic energy released by the surrounding rock during impact; u2 is the consumption of the coal bed to release the elastic energy to the surrounding rock in the impact process; u3 is the residual elastic energy released by the whole system during impact; u4 is the additional input energy.
Quasi-static load (sigma)s) Under the action of the strain, the coal bed in the post-peak stage (DE) generates strain increment delta epsilon2Corresponding to the strain increase delta epsilon generated in the roof-floor surrounding rock1Comprises the following steps:
in the formula, k1Loading stiffness for the surrounding rock before the peak; k is a radical of2The post peak unloading stiffness of the coal seam. Therefore, the total strain increment Δ ε generated by the whole system from roof to coal seam to floor is:
wherein, the ratio of the coal seam strain increment to the system integral strain increment is as follows:
when k is1+k2When equal to 0, corresponds to S in FIG. 21Point, at this time,. DELTA.. epsilon2And (/ Δ ∈ → ∞ and the impact process starts, and the roof-coal seam-floor system gradually reaches a new balance along with gradual slowing of the coal body impact destruction process, and the impact process is ended corresponding to the point S in fig. 2.
When the roof-coal bed-floor impact carrier system is superposed with dynamic load (sigma)d) When compared with the impact energy U released under quasi-static load3Its energy release will increase U4As shown in FIG. 2, thisTime equivalent to the rigidity of the surrounding rock from k1Down to k1', simultaneous impact start position from S1Advance to S2. More importantly, for accumulated plastic deformation type dynamic load, such as mine earthquake, especially far-field seismic source, the action mode is equivalent to cyclic loading and unloading, and due to the heterogeneous nature of coal rock materials, each loading and unloading caused by the mine earthquake dynamic load can cause permanent deformation of coal bodies, and when the action time of the dynamic load is long enough, the coal bodies are in S state for stress state2The coal body of the' is positioned at S under the action of superimposed dynamic load and the action of similar quasi-static load2Impact conditions under stress conditions; for stress increment type dynamic load, such as instantaneous dynamic load generated by fault slip near working face, roof fracture and the like, the action mode is equivalent to applying an instantaneous stress increment delta sigma, and the current area S is123>S3D4During the process, under the action of superimposed dynamic load, the coal body in the stress state at the position of 1 before the peak can start the impact condition similar to the impact condition in the stress state at the position of 4 after the peak under the action of quasi-static load.
The invention discloses a rock burst danger monitoring and early warning method based on dynamic and static combined stress analysis, which mainly adopts the inversion of vibration wave CT
And microseismic lesion reconstruction
Quantitative description of static load stress sigmasReconstruction of attenuation accumulation of shock waves
And mine pressure instant dynamic load stress increment
The cumulative plastic deformation type dynamic load and the stress increment type dynamic load are quantitatively described, respectively, as described in detail below.
A rock burst danger monitoring and early warning method based on dynamic and static combined stress analysis is based on a dynamic load effect and static load effect superposition induced rock burst mechanism and combines a dynamic and static combined effect to carry out rock burst danger monitoring and early warning, wherein:
the dynamic load effect is quantitatively described through dynamic load effect parameters reconstructed by vibration wave attenuation accumulation and dynamic load effect parameters of mine pressure instant dynamic load stress increment;
the static load effect is quantitatively described through a static load effect parameter inverted by the vibration wave CT and a static load effect parameter reconstructed by the microseismic damage;
the dynamic and static combined effect is characterized in that dynamic load effect parameters reconstructed by vibration wave attenuation accumulation, static load effect parameters reconstructed by vibration wave CT and the weight of each parameter of static load effect parameters reconstructed by micro-seismic damage are determined respectively, dynamic and static combined stress coefficients are calculated by adopting a weighted average method, and a final danger area and the danger degree thereof are determined comprehensively by combining the dynamic load stress coefficients of mine pressure instant dynamic load stress increment.
Preferably, the static load effect parameter of the seismic wave CT inversion is calculated by using the following formula:
in the formula (I), the compound is shown in the specification,
static load stress coefficient, v, obtained for seismic wave CT inversionPLongitudinal wave velocity, v, obtained for seismic wave CT inversionPaAnd alpha is an experimental test fitting parameter.
Preferably, the static load effect parameter of the microseismic damage reconstruction is calculated by adopting the following formula:
in the formula (I), the compound is shown in the specification,
coefficient of static load stress, σ, for microseismic damage reconstructionsssStatic load stress parameter, sigma, for microseismic damage reconstructionsssaAverage of static load stress parameter values reconstructed for microseismic damage, NMSFor microseismic frequency, the exp function is an exponential function with a natural constant e as the base,
to accumulate Benioff strain, E is the microseismic event energy, max { εEIs the maximum value of the cumulative Benioff strain, DFIs a critical damage value, corresponding to a complete damage state, εFIs the critical cumulative Benioff strain. In practice, max [ epsilon ] is calculatedEThe position is regarded as a critical damage state, and the default value is 0.95.
Preferably, the dynamic load effect parameter of the seismic wave attenuation accumulation reconstruction is calculated by adopting the following formula:
in the formula (I), the compound is shown in the specification,
dynamic load stress coefficient, σ, reconstructed for seismic wave attenuation accumulationsdsDynamic load stress parameter, sigma, reconstructed for damping accumulation of shock wavessdsaAverage value, sigma, of dynamic load stress parameter values reconstructed for damping accumulation of shock wavessdsiIs σ at the ith spatial position nodesds，EjIs the jth microseismic event energy, RijBetween the jth microseismic event and the ith spatial location nodeBeta is the attenuation coefficient of the shock wave in the formation, NMSThe microseismic frequency.
The attenuation coefficient beta of the shock wave in the stratum can be obtained by fitting and calculating the attenuation change relation of the speed amplitude of the micro-seismic probe with the distance between the micro-seismic probe and the seismic source, wherein the speed amplitude of the micro-seismic probe is not beyond the measuring range in the far field.
Preferably, the dynamic load effect parameter of the transient dynamic load stress increment of the mine pressure is calculated by adopting the following formula:
in the formula (I), the compound is shown in the specification,
dynamic load stress coefficient, P, being the instantaneous dynamic load stress increment of the mine pressureTiIs the weighted average of the stent working resistance in the ith time window (generally taking 6 hours, 8 hours, 12 hours, 24 hours and the like),
is PTiAverage value of (a) ("sigmaPIs PTiStandard deviation of (1), tpj、pjRespectively corresponding time and pressure values of the real-time monitoring sequence of the working resistance of the bracket, and n is the number of monitoring data points in a time window.
Preferably, a weighted average method is adopted to calculate the dynamic and static combined stress coefficient of the rock burst danger prediction
：
In the formula, subscripts min and max represent the minimum value and the maximum value of each parameter, respectively.
Wherein, the dynamic load effect parameter reconstructed by the attenuation accumulation of the vibration wave, the static load effect parameter reconstructed by the CT of the vibration wave and the weight (w) of each parameter of the static load effect parameter reconstructed by the microseismic damage can be determined by adopting a standard deviation methodvss、wsssAnd wsds) The calculation formula is as follows:
in the formula (I), the compound is shown in the specification,
the static load stress coefficient obtained for the inversion of the vibration wave CT is weighted by wvss；
The coefficient of static load stress for reconstructing microseismic damage is weighted by wsss；
The dynamic load stress coefficient reconstructed for the attenuation accumulation of the vibration wave is weighted by wsds(ii) a The subscript StDv indicates the standard deviation of each parameter.
The following detailed description is provided with reference to specific embodiments, as shown in fig. 3 to 9, wherein an example analysis selects a rock burst occurrence event occurring in a certain mine 2016/10/2722: 21:51, and simultaneously selects microseismic (2016/10/20-2016/10/26) and mine pressure (2016/10/01-2016/10/27) data in a period of time before the impact occurs for analysis, as shown in fig. 3, the steps of performing rock burst risk monitoring and early warning by using the method of the present invention are as follows, wherein the calculation steps of each parameter can be set by itself:
(1) performing seismic wave CT inversion analysis by using the microseismic data and the microseismic probe shown in FIG. 3 to obtain longitudinal wave velocity distribution, and then calculating by adopting a formula (5) to obtain a static load stress coefficient obtained by seismic wave CT inversion shown in FIG. 4
Results graph, in which α value takes 0.2.
(2) The static load stress coefficient of the microseismic damage reconstruction shown in figure 5 is obtained by calculation of equations (6) to (8) by using the microseismic data shown in figure 3
And (5) a result chart.
(3) By using the microseismic data shown in fig. 3, the beta value is 2.998 obtained by fitting and calculating the relation between the velocity amplitude of the microseismic probe with far field not exceeding the measuring range and the attenuation change of the distance between the microseismic probe and the seismic source, and the result is shown in fig. 6; then, the equations (9) to (10) are adopted to calculate and obtain the dynamic load stress coefficient of the seismic wave attenuation accumulation reconstruction shown in FIG. 7
And (5) a result chart.
(4) The dynamic load stress coefficient of the mineral pressure instant dynamic load stress increment shown in the figure 8 is obtained by utilizing the collected mineral pressure (2016/10/01-2016/10/27) data and calculating by adopting the formulas (11) to (12)
And (5) a result chart.
(5) According to the calculation results of fig. 4, 5 and 7, the weights (w) of the parameters are calculated by using the formulas (14) to (16)vss＝0.125、wsds0.304 and wsss0.571), and then calculating by using a formula (13) to obtain the dynamic and static combined stress coefficient for predicting the rock burst danger as shown in fig. 9
And (5) a result chart.
And (3) comprehensively determining the danger degree of the mine pressure instant dynamic load stress increment by combining the dynamic load effect parameters of the mine pressure instant dynamic load stress increment, preferably, when the value of the dynamic and static combined stress coefficient in a certain area is larger than a set first threshold value or the value of the dynamic load stress coefficient of the mine pressure instant dynamic load stress increment is larger than a set second threshold value, the system carries out early warning prompt. For example, the risk degree may be divided into several levels, and a specific level is determined according to the dynamic-static combined stress coefficient and the dynamic load stress coefficient of the mine pressure instant dynamic load stress increment, for example, in fig. 8 and 9, the risk degree is divided into 4 levels: for each stress coefficient
Or
When the value is in the interval of 0-0.25, the danger level of the area is zero;
when the value is in the interval of 0.25-0.5, the danger level of the area is weak;
when the value is in the interval of 0.5-0.75, the danger level of the area is medium;
when its value is in the interval 0.75-1, then the regional risk rating is strong.
When in actual application, the stress coefficient
A value of greater than 0.5, or
When the value of (2) is more than 0.5, the system carries out early warning prompt.
The example shows that the static load stress coefficient (figure 4) obtained by the seismic wave CT inversion, the dynamic load stress coefficient (figure 7) reconstructed by the seismic wave attenuation accumulation and the dynamic load stress coefficient (figure 8) of the dynamic load stress increment at the ore pressure moment well predict the position of the impact seismic source. The location of the impact failure zone is well predicted by the static load stress coefficient (figure 5) of the microseismic damage reconstruction; the dynamic and static combined stress coefficient (figure 9) simultaneously predicts the positions of an impact seismic source and an impact damage area, and the prediction effect is good. The dynamic and static effects and the superposition effect of the rock burst inducing mechanism are comprehensively considered, the calculation model is clear, the universality and the operability are high, the quantitative analysis of the rock burst dangerous area and the prediction of the dangerous degree of the rock burst dangerous area can be realized, the application feasibility is good, and the accurate prediction can be realized. In addition, the rapid updating and adjusting of the weight and the objective judgment of the final comprehensive prediction result have the advantage of high prediction efficiency.
The above description is only a preferred embodiment of the present invention, and not intended to limit the scope of the present invention, and all modifications of equivalent structures and equivalent processes, which are made by using the contents of the present specification and the accompanying drawings, or directly or indirectly applied to other related technical fields, are included in the scope of the present invention.
Claims (10)
1. A rock burst danger monitoring and early warning method based on dynamic and static combined stress analysis is characterized in that a rock burst danger monitoring and early warning method based on a dynamic load effect and static load effect superposition-induced rock burst mechanism and combined with a dynamic and static combined effect is carried out, wherein:
the dynamic load effect is quantitatively described through dynamic load effect parameters reconstructed by vibration wave attenuation accumulation and dynamic load effect parameters of mine pressure instant dynamic load stress increment;
the static load effect is quantitatively described through a static load effect parameter inverted by the vibration wave CT and a static load effect parameter reconstructed by the microseismic damage;
the dynamic and static combined effect is characterized in that dynamic load effect parameters reconstructed by vibration wave attenuation accumulation, static load effect parameters reconstructed by vibration wave CT and the weight of each parameter of static load effect parameters reconstructed by micro-seismic damage are determined respectively, dynamic and static combined stress coefficients are calculated by adopting a weighted average method, and a final danger area and the danger degree thereof are determined comprehensively by combining the dynamic load stress coefficients of mine pressure instant dynamic load stress increment.
2. The rock burst hazard monitoring and early warning method based on dynamic and static combined stress analysis as claimed in claim 1, wherein the dynamic load effect parameter of the seismic wave attenuation accumulation reconstruction is calculated by adopting the following formula:
in the formula (I), the compound is shown in the specification,
dynamic load stress coefficient, σ, reconstructed for seismic wave attenuation accumulationsdsDynamic load stress parameter, sigma, reconstructed for damping accumulation of shock wavessdsaAverage value, sigma, of dynamic load stress parameter values reconstructed for damping accumulation of shock wavessdsiIs σ at the ith spatial position nodesds，EjIs the jth microseismic event energy, RijIs the distance between the jth microseismic event and the ith spatial location node, beta is the attenuation coefficient of the shock wave in the formation, NMSThe microseismic frequency.
3. The rock burst hazard monitoring and early warning method based on dynamic and static combined stress analysis as claimed in claim 1, wherein the dynamic load effect parameter of the mine pressure instant dynamic load stress increment is calculated by adopting the following formula:
in the formula (I), the compound is shown in the specification,
dynamic load stress coefficient, P, being the instantaneous dynamic load stress increment of the mine pressureTiThe weighted average of the stent's operational resistance in the ith time window,
is PTiAverage value of (a) ("sigmaPIs PTiStandard deviation of (1), tpj、pjRespectively corresponding time and pressure values of the real-time monitoring sequence of the working resistance of the bracket, and n is the number of monitoring data points in a time window.
4. The rock burst hazard monitoring and early warning method based on dynamic and static combined stress analysis as claimed in claim 1, wherein the static load effect parameter of the seismic wave CT inversion is calculated by adopting the following formula:
in the formula (I), the compound is shown in the specification,
static load stress coefficient, v, obtained for seismic wave CT inversionPLongitudinal wave velocity, v, obtained for seismic wave CT inversionPaAnd alpha is an experimental test fitting parameter.
5. The rock burst hazard monitoring and early warning method based on dynamic and static combined stress analysis as claimed in claim 1, wherein the static load effect parameter of the microseismic damage reconstruction is calculated by adopting the following formula:
in the formula (I), the compound is shown in the specification,
coefficient of static load stress, σ, for microseismic damage reconstructionsssStatic load stress parameter, sigma, for microseismic damage reconstructionsssaAverage of static load stress parameter values reconstructed for microseismic damage, NMSFor microseismic frequency, the exp function is an exponential function with a natural constant e as the base,
to accumulate Benioff strain, E is the microseismic event energy, max { εEIs the maximum value of the cumulative Benioff strain, DFIs a critical damage value, corresponding to a complete damage state, εFIs the critical cumulative Benioff strain.
6. The rock burst hazard monitoring and early warning method based on dynamic and static combined stress analysis as claimed in claim 1, wherein the weights of the dynamic load effect parameter reconstructed by vibration wave attenuation accumulation, the static load effect parameter reconstructed by vibration wave CT and the static load effect parameter reconstructed by microseismic damage are determined by a standard dispersion method, and the calculation formula is as follows:
in the formula (I), the compound is shown in the specification,
the static load stress coefficient obtained for the inversion of the vibration wave CT is weighted by wvss；
The coefficient of static load stress for reconstructing microseismic damage is weighted by wsss；
The dynamic load stress coefficient reconstructed for the attenuation accumulation of the vibration wave is weighted by wsds(ii) a The subscript StDv indicates the standard deviation of each parameter.
7. The rock burst hazard monitoring and early warning method based on dynamic and static combined stress analysis as claimed in claim 6, wherein the dynamic and static combined stress coefficient NσsdCalculated using the following formula:
in the formula, subscripts min and max represent the minimum value and the maximum value of each parameter, respectively.
8. The method for monitoring and warning the danger of rock burst based on dynamic and static combined stress analysis as claimed in claim 2, wherein the attenuation coefficient β of the shock wave in the stratum is obtained by selecting the attenuation change relationship of the velocity amplitude of the far-field non-overranging microseismic probe along with the distance between the microseismic probe and the seismic source through fitting calculation.
9. The rock burst hazard monitoring and early warning method based on dynamic and static combined stress analysis as claimed in claim 5, wherein max { epsilon [ [ epsilon ] ]EThe position is regarded as critical damage state, DFIs 0.95.
10. The rock burst hazard monitoring and early warning method based on dynamic and static combined stress analysis as claimed in claim 1, wherein when the value of the dynamic and static combined stress coefficient in a certain area is greater than a set first threshold value, or the value of the dynamic load stress coefficient of the dynamic load stress increment at the moment of mine pressure is greater than a set second threshold value, the system performs early warning prompt.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010714159.5A CN111963243B (en) | 2020-07-22 | 2020-07-22 | Rock burst danger monitoring and early warning method based on dynamic and static combined stress analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010714159.5A CN111963243B (en) | 2020-07-22 | 2020-07-22 | Rock burst danger monitoring and early warning method based on dynamic and static combined stress analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111963243A CN111963243A (en) | 2020-11-20 |
| CN111963243B true CN111963243B (en) | 2021-05-07 |
Family
ID=73362633
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010714159.5A Active CN111963243B (en) | 2020-07-22 | 2020-07-22 | Rock burst danger monitoring and early warning method based on dynamic and static combined stress analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111963243B (en) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112630826B (en) * | 2020-12-04 | 2022-01-21 | 中国矿业大学 | Method for determining mining tendency influence range of working face based on microseismic monitoring |
| CN112594008B (en) * | 2020-12-17 | 2022-03-01 | 中国矿业大学 | Rock burst early warning method based on dynamic and static load disturbance superposition |
| CN113343548A (en) * | 2021-05-08 | 2021-09-03 | 北京科技大学 | Rock burst danger evaluation method based on dynamic and static load simulation |
| CN113107601A (en) * | 2021-05-19 | 2021-07-13 | 中国矿业大学 | Stope face pressure relief effect evaluation method based on multi-parameter monitoring parameters |
| CN113339072B (en) * | 2021-07-06 | 2022-07-01 | 中国矿业大学 | Blasting pressure relief effect evaluation method based on microseismic signal waveform analysis |
| CN113882901A (en) * | 2021-09-16 | 2022-01-04 | 平顶山天安煤业股份有限公司八矿 | Comprehensive early warning system based on dynamic and static load rock burst danger superposition |
| CN114087018B (en) * | 2021-11-17 | 2023-03-24 | 中国矿业大学 | Large-diameter pressure relief drilling hole accurate pressure relief method based on stress sensing |
| CN114810211B (en) * | 2021-11-26 | 2023-03-10 | 中国矿业大学 | Rock burst danger prediction method based on mine seismic group shock wave energy attenuation characteristics |
| CN114878051B (en) * | 2022-06-16 | 2023-02-03 | 中国矿业大学 | Roadway side coal impact risk assessment method based on drilling stress distribution |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015002556A2 (en) * | 2014-10-30 | 2015-01-08 | Instytut Technik Innowacy Jnych Emag | A method and system for detecting and reducing methane hazard in vicinity of a longwall |
| CN108590768A (en) * | 2018-05-11 | 2018-09-28 | 山东科技大学 | One kind is non-uniformly to cover pressure water filling coal seam fluid structurecoupling stress monitoring system |
| CN109296399A (en) * | 2018-11-01 | 2019-02-01 | 中国矿业大学(北京) | Underground goaf, which is leaked out, prevents and treats device |
| CN110737023A (en) * | 2019-10-17 | 2020-01-31 | 北京矿冶科技集团有限公司 | mining micro-seismic monitoring signal processing method |
| CN110765983A (en) * | 2019-11-07 | 2020-02-07 | 山东科技大学 | Coal pillar type rock burst risk evaluation method based on microseismic frequency spectrum characteristics |
| CN111271120A (en) * | 2020-01-22 | 2020-06-12 | 天地科技股份有限公司 | Coal mine rock burst dynamic and static load source-based danger relieving method |
-
2020
- 2020-07-22 CN CN202010714159.5A patent/CN111963243B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015002556A2 (en) * | 2014-10-30 | 2015-01-08 | Instytut Technik Innowacy Jnych Emag | A method and system for detecting and reducing methane hazard in vicinity of a longwall |
| CN108590768A (en) * | 2018-05-11 | 2018-09-28 | 山东科技大学 | One kind is non-uniformly to cover pressure water filling coal seam fluid structurecoupling stress monitoring system |
| CN109296399A (en) * | 2018-11-01 | 2019-02-01 | 中国矿业大学(北京) | Underground goaf, which is leaked out, prevents and treats device |
| CN110737023A (en) * | 2019-10-17 | 2020-01-31 | 北京矿冶科技集团有限公司 | mining micro-seismic monitoring signal processing method |
| CN110765983A (en) * | 2019-11-07 | 2020-02-07 | 山东科技大学 | Coal pillar type rock burst risk evaluation method based on microseismic frequency spectrum characteristics |
| CN111271120A (en) * | 2020-01-22 | 2020-06-12 | 天地科技股份有限公司 | Coal mine rock burst dynamic and static load source-based danger relieving method |
Non-Patent Citations (1)
| Title |
|---|
| 煤矿冲击矿压动静载的"应力场一震动波场"监测预警技术;窦林名;《岩石力学与工程学报》;20170430;第803-811页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111963243A (en) | 2020-11-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111963243B (en) | Rock burst danger monitoring and early warning method based on dynamic and static combined stress analysis | |
| Cai | Prediction and prevention of rockburst in metal mines–A case study of Sanshandao gold mine | |
| Eberhardt et al. | Numerical analysis of initiation and progressive failure in natural rock slopes—the 1991 Randa rockslide | |
| He et al. | Microseismic and electromagnetic coupling method for coal bump risk assessment based on dynamic static energy principles | |
| CN113958366B (en) | Dynamic quantitative early warning method for impact risk based on vibration-stress double-field monitoring | |
| CN109798106B (en) | Method for predicting risk of rock burst and prevention and treatment measures | |
| Masoudi et al. | Reinforcement selection for deep and high-stress tunnels at preliminary design stages using ground demand and support capacity approach | |
| Yu et al. | Experimental study of the mechanical behavior of sandstone affected by blasting | |
| Yardimci et al. | A new protective destressing technique in underground hard coal mining | |
| CN114810211A (en) | Rock burst danger prediction method based on mine seismic group vibration wave energy attenuation characteristics | |
| Vennes et al. | Geomechanical effects of stress shadow created by large-scale destress blasting | |
| Feng et al. | Deterministic seismic hazard assessment of coal fractures in underground coal mine: A case study | |
| Pan et al. | Disturbance response instability theory of rock bursts in coal mines and its application | |
| Lai et al. | Structure instability forecasting and analysis of giant rock pillars in steeply dipping thick coal seams | |
| Stacey | Support of excavations subjected to dynamic (rockburst) loading | |
| Du et al. | Experimental study of the influence of structural planes on the mechanical properties of sandstone specimens under cyclic dynamic disturbance | |
| CN109918697B (en) | Method and device for predicting mine rock burst strength | |
| Luo et al. | Early warning of rock mass instability based on multi-field coupling analysis and microseismic monitoring | |
| CN109918696B (en) | Classification method and device for rock burst strength grade | |
| Wu et al. | Reliability analysis and prediction on tunnel roof under blasting disturbance | |
| Wang et al. | Rock burst evaluation using the CRITIC algorithm-based cloud model | |
| Dai et al. | The sliding instability of rock roof block induced by vertical dynamic disturbance: An experimental study | |
| Lv et al. | Effect of multi-factor dynamic loading on gob-side entry driving during longwall face extraction: a case study | |
| Wu et al. | Mechanism of rock bursts induced by the synthetic action of “roof bending and rock pillar prying” in subvertical extra-thick coal seams | |
| Kocharyan et al. | Impact of massive blasts on stability of tectonic faults |
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 |