CN113673119A - Dynamic and static coupling evaluation method for coal mine rock burst danger based on Bayes method - Google Patents
Dynamic and static coupling evaluation method for coal mine rock burst danger based on Bayes method Download PDFInfo
- Publication number
- CN113673119A CN113673119A CN202111044168.9A CN202111044168A CN113673119A CN 113673119 A CN113673119 A CN 113673119A CN 202111044168 A CN202111044168 A CN 202111044168A CN 113673119 A CN113673119 A CN 113673119A
- Authority
- CN
- China
- Prior art keywords
- index
- dynamic
- microseismic
- rock burst
- static
- 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.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 49
- 230000003068 static effect Effects 0.000 title claims abstract description 43
- 238000011156 evaluation Methods 0.000 title claims abstract description 22
- 239000003245 coal Substances 0.000 title claims abstract description 21
- 230000008878 coupling Effects 0.000 title claims abstract description 14
- 238000010168 coupling process Methods 0.000 title claims abstract description 14
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 14
- 238000000034 method Methods 0.000 title claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims abstract description 39
- 238000012544 monitoring process Methods 0.000 claims abstract description 37
- 230000002159 abnormal effect Effects 0.000 claims abstract description 14
- 238000013398 bayesian method Methods 0.000 claims abstract description 13
- 238000005553 drilling Methods 0.000 claims abstract description 9
- 238000006243 chemical reaction Methods 0.000 claims abstract description 6
- 238000010606 normalization Methods 0.000 claims abstract description 4
- 238000012545 processing Methods 0.000 claims abstract description 4
- 238000013210 evaluation model Methods 0.000 claims description 11
- 238000005065 mining Methods 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 claims description 3
- 230000004927 fusion Effects 0.000 abstract description 8
- 238000012821 model calculation Methods 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/01—Probabilistic graphical models, e.g. probabilistic networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- Evolutionary Computation (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Development Economics (AREA)
- General Engineering & Computer Science (AREA)
- Educational Administration (AREA)
- Entrepreneurship & Innovation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Marine Sciences & Fisheries (AREA)
- Agronomy & Crop Science (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Physics (AREA)
- General Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Mathematical Optimization (AREA)
- Primary Health Care (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Animal Husbandry (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Artificial Intelligence (AREA)
- Algebra (AREA)
- Probability & Statistics with Applications (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
Abstract
The invention provides a Bayesian method-based dynamic and static coupling evaluation method for coal mine rock burst danger, which comprises the following steps of: s1, obtaining dynamic index parameter data; s2, obtaining static index parameter data; s3, carrying out normalization processing on the dynamic and static index parameter data; s4, carrying out abnormal index conversion on each microseismic index obtained by monitoring the microseismic sensor; s5, carrying out comprehensive risk calculation on the drilling stress index value obtained by the drilling stress sensor; s6, performing fusion calculation on each dynamic and static index through a Bayes probability combination model; and S7, carrying out rock burst grade division on the combined model calculation result to realize intelligent grading early warning. According to the method and the device, data indexes of each monitoring system are effectively fused, dynamic and static indexes are comprehensively considered, dynamic weight calculation of each index based on a time sequence is realized, the safety early warning capability of rock burst is improved, and the problems that an existing multi-index safety early warning threshold is difficult to determine and the data fusion degree of each monitoring system is low are solved.
Description
Technical Field
The invention relates to the technical field of coal mine underground rock burst dynamic disaster safety monitoring and early warning, in particular to a Bayesian method-based dynamic and static coupling evaluation method for coal mine rock burst danger.
Background
Along with the increase of the coal mining depth, the stress of a deep coal seam is increased, and the possibility of disasters is higher, so that the underground monitoring and early warning device is particularly important for underground common dynamic disasters. In a deep mining area, rock burst is taken as a typical common coal rock dynamic disaster, and has important engineering value for underground monitoring and early warning.
The research work of monitoring and early warning of rock burst has achieved good results. In particular, monitoring systems for micro-seismic, earth sound, electromagnetic radiation, drilling stress and the like are widely popularized and applied, and have good application effects in many mine enterprises. At present, the monitoring and early warning of the rock burst are mainly comprehensively analyzed in two aspects of dynamic state and static state, the static state is analyzed according to two major indexes of geological conditions, comprehensive indexes of mining conditions and possibility indexes, and the dynamic indexes are used for researching the rock burst through data acquired by monitoring systems of microseismic events, energy, stress and the like. The inventor of the invention finds that the existing rock burst safety monitoring and early warning system has the following problems through research:
(1) monitoring data mining inadequacy
In the process of preventing and controlling rock burst, at present, China accumulates a large amount of geological data and monitoring data, the data fusion degree of each monitoring system is low, and how to use a Bayesian method to excavate beneficial data capable of safely early warning impact disasters through big data is lack of a reliable and effective calculation method at present.
(2) Determination of safety precaution threshold
When the constant analysis of each index of rock burst is carried out, most of safety early warning threshold values influencing the rock burst danger are obtained according to the actual situation on site, are biased to engineering experience and have certain uniqueness, and the safety early warning threshold values of each index are difficult to determine and random and cannot be applied to all projects.
(3) Linkage property lacking dynamic and static indexes
For a large number of monitoring indexes, one or two independent indexes are usually adopted as an early warning threshold value at present, and more than constant and more than monitoring data indexes can not be represented uniformly, which is obviously insufficient.
Disclosure of Invention
The invention provides a dynamic and static coupling evaluation method for coal mine rock burst danger based on a Bayesian method, aiming at the technical problems that the existing rock burst safety monitoring and early warning system is low in data fusion degree of each monitoring system, each index safety early warning threshold value is difficult to determine and random, and dynamic indexes and static indexes lack linkage.
In order to solve the technical problems, the invention adopts the following technical scheme:
the dynamic and static coupling evaluation method for the coal mine rock burst danger based on the Bayesian method comprises the following steps:
s1, transmitting and storing the data received by the underground micro-seismic sensor and the borehole stress sensor to a monitoring system through a data transmission device to obtain dynamic index parameter data with time series characteristics;
s2, performing comprehensive index evaluation calculation from two aspects of geological factors and mining technologies, performing probability index calculation from stress and coal seam impact tendency, and obtaining regional static index parameter data based on the comprehensive index evaluation calculation and the probability index calculation;
s3, preliminarily carrying out normalization processing on each dynamic index parameter data and each static index parameter data by introducing an abnormal conversion function, and obtaining multi-parameter abnormal indexes uniformly constrained in a dimensionless quantity in a 0-1 closed interval;
s4, converting the abnormal indexes of the microseismic indexes obtained by monitoring the microseismic sensors at different moments according to the following formula:
in the formula, λijThe membership degree of the corresponding index in the statistical time window t is in a value range of 0-1;
s5, carrying out comprehensive risk calculation on drilling stress index values obtained by monitoring the drilling stress sensor at different moments according to the following formula to obtain the risk level at any moment:
PI=k1I1+k2I2+k3I3
in the formula I1The risk index is the magnitude of the stress value of the measuring point; i is2The stress amplitude value is a dangerous index of the stress amplitude value, namely a value obtained by calculating the difference between the stress value of the current measuring point and the stress value of the initial setting point; i is3A risk indicator which is a stress acceleration value; k is a radical of1、k2、k3Are respectively I1、I2、I3The weight ratio of (a);
s6, analyzing each dynamic and static index after the comprehensive risk calculation by using a Bayesian probability model, and comprehensively evaluating the risk by using a combined model; setting each dynamic and static index as independent evaluation model, and assuming the weight of the kth evaluation modelIs PkIf n independent evaluation models coexist, the posterior probability of the combined model is:
wherein O (C) represents the posterior probability of event C; p (C) represents the prior probability of event C, typically set to 1/n;no event occurs; here, the weight of each evaluation model is determined as a probability according to the result obtained in steps S2, S3, S5;
s7, obtaining a dimensionless quantity P belonging to a 0-1 closed interval through calculation of the combined model, carrying out grade division on the dimensionless quantity P according to a preset rock burst risk grade division standard to obtain rock burst grades at all times, and realizing intelligent grading early warning according to the rock burst grades.
Compared with the prior art, the dynamic and static coupling evaluation method for the coal mine rock burst danger based on the Bayes method comprises the steps of firstly collecting regional dynamic index parameter data and regional static index parameter data, then converting each dynamic and static index by introducing an abnormal conversion function to obtain a plurality of constant abnormal indexes, constraining each dynamic and static index in a dimensionless quantity of a 0-1 closed interval, then performing fusion calculation on each dynamic and static index by a mathematical probability statistical method to obtain the comprehensive danger degree of each index weight and rock burst, and performing safety early warning according to a preset rock burst danger grade division standard. The method is based on a Bayesian method, establishes multi-index-fused rock burst disaster safety monitoring and early warning, effectively fuses data indexes of each monitoring system, comprehensively considers dynamic indexes and static indexes, realizes dynamic weight calculation of each index based on time series, reduces uncertainty, randomness and fuzziness in the calculation process, and improves the rock burst safety early warning capability, so that the problem of poor unified effect of fusion of each index and a module of the existing online monitoring system is effectively solved, and the technical problem that the existing multi-index safety early warning threshold is difficult to determine is solved.
Further, the microseismic index obtained by monitoring the microseismic sensor in the step S4 includes a microseismic intensity factor, a microseismic equivalent energy level parameter, a microseismic timing factor, a microseismic b value, and a microseismic a (b) value.
Further, the parameter λ in the step S4ijThe calculation method of (2) is as follows:
for the forward anomaly indicators:
λij=(Sij-Smin)/(Smax-Smin)
for negative anomaly indicators:
λij=(Smax-Sij)/(Smax-Smin)
in the formula, SijIs the value of the microseismic index at a certain moment in the time window, SmaxThe maximum value S of each microseismic index is taken at each moment in the time windowminAnd taking the minimum value of the values of each microseismic index at each moment in the time window.
Further, the preset criterions for classifying the dangerousness of rock burst in step S7 specifically classify the possibility of the dimensionless quantity P into four grades, i.e., none, weak, medium and strong, the criterions for classifying the none, weak, medium and strong grades are respectively 0 ≦ P <0.25, 0.25 ≦ P <0.5, 0.5 ≦ P <0.75 and 0.75 ≦ P <1, and the four grades are respectively pre-warned with green, yellow, orange and red.
Drawings
FIG. 1 is a schematic flow chart of a dynamic and static coupling evaluation method for coal mine rock burst hazard based on a Bayesian method provided by the invention.
Detailed Description
In order to make the technical means, the creation characteristics, the achievement purposes and the effects of the invention easy to understand, the invention is further explained below by combining the specific drawings.
Referring to fig. 1, the dynamic and static coupling evaluation method for coal mine rock burst hazard based on the bayesian method provided by the invention comprises the following steps:
s1, transmitting and storing the data received by the underground micro-seismic sensor and the borehole stress sensor to a monitoring system through a data transmission device to obtain dynamic index parameter data with time series characteristics; the data of the database is read in real time by interpreting and coding the types of the monitoring system database and then compiling an interface, and the data of multiple monitoring systems are fused and analyzed uniformly.
S2, performing comprehensive index evaluation calculation from two aspects of geological factors and mining technologies, performing probability index calculation from stress and coal seam impact tendency, and obtaining regional static index parameter data based on the comprehensive index evaluation calculation and the probability index calculation, namely discussing static indexes with dangerousness from the comprehensive index evaluation calculation and the probability index calculation, wherein the specific calculation processes of the comprehensive index evaluation calculation and the probability index calculation are the prior art well known to those skilled in the art, and are not repeated herein.
S3, performing normalization processing on each dynamic index parameter data and each static index parameter data preliminarily by introducing an abnormal conversion function, and obtaining multi-parameter abnormal indexes uniformly constrained within a dimensionless quantity in a 0-1 closed interval.
S4, monitoring the microseismic sensors to obtain each microseismic index: the index values of the microseismic intensity factor, the microseismic equivalent energy level parameter, the microseismic time sequence factor, the microseismic b value and the microseismic A (b) value at different moments are converted into abnormal indexes according to the following formula:
in the formula, λijThe membership degree of the corresponding index in the statistical time window t is in a value range of 0-1; parameter lambdaijThe calculation method of (2) is as follows:
for the forward anomaly indicators:
λij=(Sij-Smin)/(Smax-Smin)
for negative anomaly indicators:
λij=(Smax-Sij)/(Smax-Smin)
in the formula, SijIs the value of the microseismic index at a certain moment in the time window, SmaxThe maximum value S of each microseismic index is taken at each moment in the time windowminAnd taking the minimum value of the values of each microseismic index at each moment in the time window.
S5, carrying out comprehensive risk calculation on drilling stress index values obtained by monitoring the drilling stress sensor at different moments according to the following formula to obtain the risk level at any moment:
PI=k1I1+k2I2+k3I3
in the formula I1The risk index is the magnitude of the stress value of the measuring point; i is2The stress amplitude value is a dangerous index of the stress amplitude value, namely a value obtained by calculating the difference between the stress value of the current measuring point and the stress value of the initial setting point; i is3A risk indicator which is a stress acceleration value; k is a radical of1、k2、k3Are respectively I1、I2、I3The weight ratio of (2).
S6, after dimensionless quantities of dynamic and static indexes of the rock burst are calculated, analyzing the dynamic and static indexes after comprehensive risk calculation by using a Bayes probability model, and comprehensively evaluating risks by using a combined model; setting each dynamic and static index as independent evaluation model, and assuming that the weight of the kth evaluation model is PkIf n independent evaluation models coexist, the posterior probability of the combined model is:
wherein O (C) represents the posterior probability of event C; p (C) represents the prior probability of event C, typically set to 1/n;no event occurs; here, the weight of each evaluation model is determined as a probability according to the result obtained in steps S2, S3, S5;
s7, obtaining a dimensionless quantity P belonging to a 0-1 closed interval through calculation of the combined model, carrying out grade division on the dimensionless quantity P according to a preset rock burst risk grade division standard to obtain rock burst grades at all times, and realizing intelligent grading early warning according to the rock burst grades.
Compared with the prior art, the dynamic and static coupling evaluation method for the coal mine rock burst danger based on the Bayes method comprises the steps of firstly collecting regional dynamic index parameter data and regional static index parameter data, then converting each dynamic and static index by introducing an abnormal conversion function to obtain a plurality of constant abnormal indexes, constraining each dynamic and static index in a dimensionless quantity of a 0-1 closed interval, then performing fusion calculation on each dynamic and static index by a mathematical probability statistical method to obtain the comprehensive danger degree of each index weight and rock burst, and performing safety early warning according to a preset rock burst danger grade division standard. The method is based on a Bayesian method, establishes multi-index-fused rock burst disaster safety monitoring and early warning, effectively fuses data indexes of each monitoring system, comprehensively considers dynamic indexes and static indexes, realizes dynamic weight calculation of each index based on time series, reduces uncertainty, randomness and fuzziness in the calculation process, and improves the rock burst safety early warning capability, so that the problem of poor unified effect of fusion of each index and a module of the existing online monitoring system is effectively solved, and the technical problem that the existing multi-index safety early warning threshold is difficult to determine is solved.
As a specific example, the preset criterions for the risk rating of rock burst at step S7 specifically classify the possibility of the dimensionless quantity P into four grades, i.e., none, weak, medium and strong, and the criterions for classification at each of the grades of none, weak, medium and strong are 0 ≦ P <0.25, 0.25 ≦ P <0.5, 0.5 ≦ P <0.75, 0.75 ≦ P <1, as shown in table 1 below:
TABLE 1 grade of rock burst hazard
The four grades of none, weak, medium and strong adopt green, yellow, orange and red respectively to carry out early warning, thereby realizing intelligent grading safety early warning.
Finally, the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting, although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications or equivalent substitutions may be made to the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, and all of them should be covered in the claims of the present invention.
Claims (4)
1. The method for evaluating the dynamic and static coupling of the coal mine rock burst dangerousness based on the Bayesian method is characterized by comprising the following steps of:
s1, transmitting and storing the data received by the underground micro-seismic sensor and the borehole stress sensor to a monitoring system through a data transmission device to obtain dynamic index parameter data with time series characteristics;
s2, performing comprehensive index evaluation calculation from two aspects of geological factors and mining technologies, performing probability index calculation from stress and coal seam impact tendency, and obtaining regional static index parameter data based on the comprehensive index evaluation calculation and the probability index calculation;
s3, preliminarily carrying out normalization processing on each dynamic index parameter data and each static index parameter data by introducing an abnormal conversion function, and obtaining multi-parameter abnormal indexes uniformly constrained in a dimensionless quantity in a 0-1 closed interval;
s4, converting the abnormal indexes of the microseismic indexes obtained by monitoring the microseismic sensors at different moments according to the following formula:
in the formula, λijThe membership degree of the corresponding index in the statistical time window t is in a value range of 0-1;
s5, carrying out comprehensive risk calculation on drilling stress index values obtained by monitoring the drilling stress sensor at different moments according to the following formula to obtain the risk level at any moment:
PI=k1I1+k2I2+k3I3
in the formula I1The risk index is the magnitude of the stress value of the measuring point; i is2The stress amplitude value is a dangerous index of the stress amplitude value, namely a value obtained by calculating the difference between the stress value of the current measuring point and the stress value of the initial setting point; i is3A risk indicator which is a stress acceleration value; k is a radical of1、k2、k3Are respectively I1、I2、I3The weight ratio of (a);
s6, analyzing each dynamic and static index after the comprehensive risk calculation by using a Bayesian probability model, and comprehensively evaluating the risk by using a combined model; setting each dynamic and static index as independent evaluation model, and assuming that the weight of the kth evaluation model is PkIf n independent evaluation models coexist, the posterior probability of the combined model is:
wherein O (C) represents the posterior probability of event C; p (C) represents the prior probability of event C, typically set to 1/n;no event occurs; here, the weight of each evaluation model is determined as a probability according to the result obtained in steps S2, S3, S5;
s7, obtaining a dimensionless quantity P belonging to a 0-1 closed interval through calculation of the combined model, carrying out grade division on the dimensionless quantity P according to a preset rock burst risk grade division standard to obtain rock burst grades at all times, and realizing intelligent grading early warning according to the rock burst grades.
2. The Bayesian method-based coal mine rock burst hazard dynamic-static coupling evaluation method as recited in claim 1, wherein microseismic indexes obtained by monitoring of the microseismic sensor in step S4 include a microseismic intensity factor, a microseismic equivalent energy level parameter, a microseismic time sequence factor, a microseismic b value and a microseismic A (b) value.
3. The Bayesian method-based coal mine rock burst hazard dynamic-static coupling evaluation method according to claim 1, wherein in the step S4, parameter λ isijThe calculation method of (2) is as follows:
for the forward anomaly indicators:
λij=(Sij-Smin)/(Smax-Smin)
for negative anomaly indicators:
λij=(Smax-Sij)/(Smax-Smin)
in the formula, SijIs the value of the microseismic index at a certain moment in the time window, SmaxThe maximum value S of each microseismic index is taken at each moment in the time windowminAnd taking the minimum value of the values of each microseismic index at each moment in the time window.
4. The Bayesian method-based coal mine rock burst risk dynamic-static coupling evaluation method according to claim 1, wherein in step S7, preset rock burst risk grade division criteria specifically divide the possibility of a dimensionless quantity P into four grades of none, weak, medium and strong, the division criteria values of the none, weak, medium and strong grades are respectively 0. ltoreq.P <0.25, 0.25. ltoreq.P <0.5, 0.5. ltoreq.P <0.75, 0.75. ltoreq.P <1, and the four grades of none, weak, medium and strong adopt green, yellow, orange and red for early warning respectively.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111044168.9A CN113673119B (en) | 2021-09-07 | 2021-09-07 | Coal mine rock burst danger dynamic and static coupling evaluation method based on Bayesian method |
ZA2022/04188A ZA202204188B (en) | 2021-09-07 | 2022-04-13 | Dynamic and static coupling method for evaluating coal mine rock-burst hazards based on bayesian method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111044168.9A CN113673119B (en) | 2021-09-07 | 2021-09-07 | Coal mine rock burst danger dynamic and static coupling evaluation method based on Bayesian method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113673119A true CN113673119A (en) | 2021-11-19 |
CN113673119B CN113673119B (en) | 2023-08-18 |
Family
ID=78548656
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111044168.9A Active CN113673119B (en) | 2021-09-07 | 2021-09-07 | Coal mine rock burst danger dynamic and static coupling evaluation method based on Bayesian method |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN113673119B (en) |
ZA (1) | ZA202204188B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114352355A (en) * | 2022-03-15 | 2022-04-15 | 中国矿业大学(北京) | Deep mine rock burst control method |
CN114595581A (en) * | 2022-03-11 | 2022-06-07 | 重庆地质矿产研究院 | Regional geological disaster danger probability model based on influence factor dynamic weight |
CN115511379A (en) * | 2022-10-28 | 2022-12-23 | 北京科技大学 | Dynamic dividing method and device for rock burst dangerous area |
CN115860582A (en) * | 2023-02-28 | 2023-03-28 | 山东科技大学 | Intelligent impact risk early warning method based on self-adaptive lifting algorithm |
CN118154046A (en) * | 2024-05-10 | 2024-06-07 | 太原向明智控科技有限公司 | Top plate pressure grade dividing method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110020749A (en) * | 2019-03-20 | 2019-07-16 | 中国矿业大学 | A kind of comprehensive intelligent prediction technique based on microseism fractal prediction danger of burst grade |
CN110135079A (en) * | 2019-05-20 | 2019-08-16 | 中国石油大学(华东) | A kind of macroscopical photoelastic evaluation method and system of offshore oil well control equipment |
CN112483176A (en) * | 2020-11-26 | 2021-03-12 | 绍兴文理学院 | Rock burst multi-parameter early warning method based on fuzzy mathematics and microseismic monitoring |
-
2021
- 2021-09-07 CN CN202111044168.9A patent/CN113673119B/en active Active
-
2022
- 2022-04-13 ZA ZA2022/04188A patent/ZA202204188B/en unknown
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110020749A (en) * | 2019-03-20 | 2019-07-16 | 中国矿业大学 | A kind of comprehensive intelligent prediction technique based on microseism fractal prediction danger of burst grade |
CN110135079A (en) * | 2019-05-20 | 2019-08-16 | 中国石油大学(华东) | A kind of macroscopical photoelastic evaluation method and system of offshore oil well control equipment |
CN112483176A (en) * | 2020-11-26 | 2021-03-12 | 绍兴文理学院 | Rock burst multi-parameter early warning method based on fuzzy mathematics and microseismic monitoring |
Non-Patent Citations (4)
Title |
---|
张勇;蔡辉;杨永杰;潘立友;: "基于信息融合技术的冲击地压算法模型研究", 中国安全生产科学技术, no. 03, pages 42 - 47 * |
杨天鸿;王赫;董鑫;刘飞跃;张鹏海;邓文学;: "露天矿边坡稳定性智能评价研究现状、存在问题及对策", 煤炭学报, no. 06, pages 375 - 393 * |
王彦彬;孙韶光;: "基于网格搜索和ELM的冲击地压危险等级预测", 中国安全科学学报, no. 08, pages 101 - 105 * |
石永奎;李鹏瑞;王英;甄恩泽;张云朝;: "基于Weka的冲击地压预测方法分析", 煤炭技术, no. 05, pages 225 - 227 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114595581A (en) * | 2022-03-11 | 2022-06-07 | 重庆地质矿产研究院 | Regional geological disaster danger probability model based on influence factor dynamic weight |
CN114352355A (en) * | 2022-03-15 | 2022-04-15 | 中国矿业大学(北京) | Deep mine rock burst control method |
CN114352355B (en) * | 2022-03-15 | 2022-06-24 | 中国矿业大学(北京) | Deep mine rock burst control method |
CN115511379A (en) * | 2022-10-28 | 2022-12-23 | 北京科技大学 | Dynamic dividing method and device for rock burst dangerous area |
CN115860582A (en) * | 2023-02-28 | 2023-03-28 | 山东科技大学 | Intelligent impact risk early warning method based on self-adaptive lifting algorithm |
CN118154046A (en) * | 2024-05-10 | 2024-06-07 | 太原向明智控科技有限公司 | Top plate pressure grade dividing method |
Also Published As
Publication number | Publication date |
---|---|
CN113673119B (en) | 2023-08-18 |
ZA202204188B (en) | 2022-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113673119B (en) | Coal mine rock burst danger dynamic and static coupling evaluation method based on Bayesian method | |
CN103291364B (en) | Microseismic multidimensional information comprehensive time sequence early warning method for rock burst | |
CN103410568B (en) | Dynamic mine disaster integral early warning method and device | |
CN104732070B (en) | A kind of rock burst grade prediction technique based on information vector machine | |
CN112949202A (en) | Bayesian network-based rockburst probability prediction method | |
CN115423213A (en) | Intelligent analysis and early warning system for safety risks of smart mine based on data analysis | |
Kabiesz et al. | Application of rule-based models for seismic hazard prediction in coal mines. | |
CN110132218B (en) | Multi-level early warning determination method based on slope stability monitoring | |
CN102913285B (en) | Mine roof hazard warning method | |
CN110211338A (en) | A kind of method for early warning of monomer pipeline landslide disaster | |
Wang et al. | Risk Assessment of Water Inrush in Karst Tunnels Based on the Efficacy Coefficient Method. | |
Liang et al. | Prediction model for coal-gas outburst using the genetic projection pursuit method | |
CN112483176B (en) | Rockburst multi-parameter early warning method based on fuzzy mathematics and microseismic monitoring | |
CN113792997A (en) | Method for judging risk level of geological disaster caused by water inrush and mud inrush in fault zone | |
CN115980826A (en) | Rock burst intensity prediction method based on weighted meta-heuristic combined model | |
CN116044501A (en) | Advanced geological forecast dynamic monitoring and early warning system and method | |
CN114294062A (en) | Time-space dynamic comprehensive early warning method for rock burst | |
CN109325470A (en) | Working face in the pit homework type intelligent identification Method based on gas density parameter | |
Qiu et al. | Rockburst prediction based on distance discrimination method and optimization technology-based weight calculation method | |
CN116402959A (en) | Tunnel support design construction method, system and storage medium based on early warning model | |
CN108319915B (en) | Multi-time-window simplified form identification method for dynamically adjusting rock burst signal threshold | |
CN114970921A (en) | Early warning method and system for severity of drilling blockage | |
CN108223010B (en) | Mine dynamic disaster integrated early warning method and device | |
CN113586157B (en) | Extraction working face salient dangerous area rapid division method based on Kriging interpolation | |
CN114233386B (en) | Coal mine disaster early warning method based on multi-parameter risk judgment database |
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 |