CN113673119B - Coal mine rock burst danger dynamic and static coupling evaluation method based on Bayesian method - Google Patents
Coal mine rock burst danger dynamic and static coupling evaluation method based on Bayesian method Download PDFInfo
- Publication number
- CN113673119B CN113673119B CN202111044168.9A CN202111044168A CN113673119B CN 113673119 B CN113673119 B CN 113673119B CN 202111044168 A CN202111044168 A CN 202111044168A CN 113673119 B CN113673119 B CN 113673119B
- 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.)
- Active
Links
- 239000011435 rock Substances 0.000 title claims abstract description 49
- 230000003068 static effect Effects 0.000 title claims abstract description 45
- 239000003245 coal Substances 0.000 title claims abstract description 20
- 238000011156 evaluation Methods 0.000 title claims abstract description 20
- 238000013398 bayesian method Methods 0.000 title claims abstract description 16
- 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
- 238000012544 monitoring process Methods 0.000 claims abstract description 40
- 238000004364 calculation method Methods 0.000 claims abstract description 37
- 230000002159 abnormal effect Effects 0.000 claims abstract description 14
- 238000006243 chemical reaction Methods 0.000 claims abstract description 9
- 238000005553 drilling Methods 0.000 claims abstract description 9
- 238000012545 processing Methods 0.000 claims abstract description 5
- 238000010606 normalization Methods 0.000 claims abstract description 4
- 238000013210 evaluation model Methods 0.000 claims description 12
- 238000005065 mining Methods 0.000 claims description 6
- 238000000034 method Methods 0.000 claims description 5
- 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 6
- 230000000694 effects Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000010354 integration Effects 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 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
- 238000011160 research Methods 0.000 description 1
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)
- Economics (AREA)
- Strategic Management (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Entrepreneurship & Innovation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Agronomy & Crop Science (AREA)
- Geology (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Geometry (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Computer Hardware Design (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Quality & Reliability (AREA)
- Probability & Statistics with Applications (AREA)
- Algebra (AREA)
- Artificial Intelligence (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
Abstract
The application provides a coal mine rock burst risk dynamic and static coupling evaluation method based on a Bayesian method, which comprises the following steps: s1, obtaining dynamic index parameter data; s2, static index parameter data are obtained; 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 dangerous calculation on the drilling stress index value obtained by the drilling stress sensor; s6, carrying out fusion calculation on each dynamic and static index through a Bayesian probability combination model; and S7, performing rock burst grading on the calculation result of the combined model to realize intelligent grading early warning. According to the application, the data indexes of all monitoring systems are effectively fused, dynamic and static indexes are comprehensively considered, dynamic weight calculation of all indexes based on time sequences is realized, the rock burst safety early warning capability is improved, and the problems that the existing multi-index safety early warning threshold is difficult to determine and the data fusion degree of all monitoring systems is low are solved.
Description
Technical Field
The application relates to the technical field of coal mine underground rock burst dynamic disaster safety monitoring and early warning, in particular to a coal mine rock burst danger dynamic and static coupling evaluation method based on a Bayesian method.
Background
Along with the increase of the mining depth of the coal mine, the stress of the deep coal seam is increased, and the possibility of disasters is also higher, so that the underground monitoring and early warning of underground common dynamic disasters are particularly important. In deep mining areas, rock burst is taken as a typical common coal-rock dynamic disaster, and has important engineering value for underground monitoring and early warning.
Good results are obtained in the work of rock burst monitoring, early warning and research. The monitoring system for microseismic, ground sound, electromagnetic radiation, drilling stress and the like is widely popularized and applied, and good application effects are achieved in many mine enterprises. At present, rock burst monitoring and early warning mainly adopts dynamic and static comprehensive analysis, wherein the static comprehensive analysis is carried out according to two major indexes of geological conditions, mining condition comprehensive indexes and possibility indexes, and the dynamic indexes are data obtained through microseismic events, energy, stress and other monitoring systems to study the rock burst. The inventor of the application finds that the present rock burst safety monitoring and early warning system has the following problems:
(1) Monitoring data mining inadequacy
In the process of preventing rock burst, a large amount of geological data and monitoring data are accumulated in China at present, the data fusion degree of each monitoring system is low, and a Bayesian method is applied to discover beneficial data capable of pre-warning the safety of the impact disaster through the large data, so that a reliable and effective calculation method is lacking at present.
(2) Determination of safety precaution threshold
When the constant analysis of each index of the rock burst is carried out, most of safety early warning thresholds influencing the rock burst risk are obtained according to the actual conditions of the site, the safety early warning thresholds are biased to engineering experience, have certain uniqueness, are difficult to determine and have randomness, and cannot be suitable for all projects.
(3) Linkage lacking dynamic and static indexes
For a large number of monitoring indexes, one or two independent indexes are often adopted as early warning thresholds at present, and a plurality of constant and multi-monitoring data indexes cannot be uniformly characterized, which is obviously insufficient.
Disclosure of Invention
Aiming at the technical problems that the existing rock burst safety monitoring and early warning system has low data fusion degree of each monitoring system, each index safety early warning threshold has difficulty in determining and randomness, and dynamic indexes and static indexes lack linkage, the application provides a coal mine rock burst risk dynamic and static coupling evaluation method based on a Bayesian method.
In order to solve the technical problems, the application adopts the following technical scheme:
a coal mine rock burst risk dynamic and static coupling evaluation method based on a Bayesian method comprises the following steps:
s1, monitoring received data by a microseismic sensor and a borehole stress sensor which are installed underground, transmitting and storing the data to a monitoring system through a data transmission device, and obtaining dynamic index parameter data with time sequence characteristics;
s2, performing comprehensive index evaluation calculation from two aspects of geological factors and mining technologies, performing possibility 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 possibility index calculation;
s3, carrying out normalization processing on the dynamic index parameter data and the static index parameter data preliminarily by introducing an abnormal conversion function to obtain multi-parameter abnormal indexes uniformly constrained in non-dimensional quantities of 0-1 closed interval;
s4, carrying out abnormal index conversion on index values of each microseismic index obtained through monitoring by the microseismic sensor at different moments according to the following formula:
wherein lambda is ij The value range is 0-1 for the membership degree of the corresponding index in the statistical time window t;
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 risk levels at any moment:
P I =k 1 I 1 +k 2 I 2 +k 3 I 3
wherein I is 1 The stress value is a dangerous index of the stress value of the measuring point; i 2 The stress increment value is a dangerous index, 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 3 Is a dangerous index of stress acceleration value; k (k) 1 、k 2 、k 3 Respectively is I 1 、I 2 、I 3 Weight ratio of (2);
s6, analyzing all dynamic and static indexes after comprehensive risk calculation by using a Bayesian probability model, and comprehensively evaluating the risk by using a combined model; let each dynamic and static index be the independent evaluation model, assume that the kth evaluation model weight is P k And coexisting n mutually independent evaluation models, wherein the posterior probability of the occurrence of the combined model is as follows:
wherein O (C) represents the posterior probability of event C; p (C) represents the prior probability of event C, and is generally set to 1/n;is that the event does not occur; the weights of the evaluation models are obtained in accordance with steps S2, S3, S5The result is taken as probability;
and S7, calculating the combined model to obtain a dimensionless quantity P belonging to a closed zone of 0-1, grading the dimensionless quantity P according to a preset rock burst risk grading 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 coal mine rock burst risk dynamic and static coupling evaluation method based on the Bayesian method provided by the application has the advantages that firstly, regional dynamic index parameter data and regional static index parameter data are collected, then, abnormal conversion functions are introduced to convert all dynamic and static indexes to obtain multi-constant abnormal indexes, so that all dynamic and static indexes are constrained in dimensionless quantities of a 0-1 closed zone, then, all dynamic and static indexes are subjected to fusion calculation through a mathematical probability statistical method, the comprehensive risk degree of all index weights and rock burst is obtained, and safety early warning is carried out according to a preset rock burst risk grade division standard. According to the application, based on a Bayesian method, the multi-index integrated rock burst disaster safety monitoring and early warning is established, the data indexes of all monitoring systems are effectively integrated, meanwhile, the dynamic indexes and the static indexes are comprehensively considered, the dynamic weight calculation of all indexes based on a time sequence is realized, the uncertainty, the randomness and the ambiguity in the calculation process are reduced, and the rock burst safety early warning capability is improved, so that the problem that the integration unification effect of all indexes and modules of the existing online monitoring systems is poor 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 the microseismic sensor monitoring in step S4 includes 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.
Further, the parameter lambda in step S4 ij The calculation method of (2) is as follows:
for the forward anomaly index:
λ ij =(S ij -S min )/(S max -S min )
for negative anomaly index:
λ ij =(S max -S ij )/(S max -S min )
wherein S is ij The value of the microseismic index at a certain moment in a time window is S max For the maximum value of the values of each microseismic index at each moment in the time window, S min And taking the minimum value of the values of each microseismic index at each moment in the time window.
Further, the preset rock burst risk grade classification standard in the step S7 specifically classifies the possibility of the dimensionless quantity P into four grades, namely, no grade, weak grade, medium grade and strong grade, wherein the classification standard values of the no grade, the weak grade, the medium grade and the strong grade 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 early-warned by adopting green, yellow, orange and red.
Drawings
FIG. 1 is a schematic flow chart of a coal mine rock burst risk dynamic and static coupling evaluation method based on a Bayesian method.
Detailed Description
The application is further described with reference to the following detailed drawings in order to make the technical means, the creation characteristics, the achievement of the purpose and the effect of the implementation of the application easy to understand.
Referring to fig. 1, the method for evaluating dynamic and static coupling of coal mine rock burst dangers based on a bayesian method provided by the application comprises the following steps:
s1, monitoring received data by a microseismic sensor and a borehole stress sensor which are installed underground, transmitting and storing the data to a monitoring system through a data transmission device, and obtaining dynamic index parameter data with time sequence characteristics; and the data of the multiple monitoring systems are fused and uniformly analyzed by performing interpretation coding on the database types of the monitoring systems and then writing an interface.
S2, performing comprehensive index evaluation calculation from two aspects of geological factors and mining technology, 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 dangerous static indexes from the comprehensive index evaluation calculation and the probability index calculation, wherein specific calculation processing of the comprehensive index evaluation calculation and the probability index calculation is a prior art well known to a person skilled in the art, and is not repeated herein.
S3, carrying out normalization processing on the dynamic index parameter data and the static index parameter data preliminarily by introducing an abnormal conversion function to obtain multi-parameter abnormal indexes uniformly constrained in the dimensionless range of 0-1 closed zone.
S4, monitoring each obtained microseism index by the microseism sensor: the abnormal index conversion is carried out on the index values at different moments according to the following formula:
wherein lambda is ij The value range is 0-1 for the membership degree of the corresponding index in the statistical time window t; parameter lambda ij The calculation method of (2) is as follows:
for the forward anomaly index:
λ ij =(S ij -S min )/(S max -S min )
for negative anomaly index:
λ ij =(S max -S ij )/(S max -S min )
wherein S is ij The value of the microseismic index at a certain moment in a time window is S max For the maximum value of the values of each microseismic index at each moment in the time window, S min And 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 risk levels at any moment:
P I =k 1 I 1 +k 2 I 2 +k 3 I 3
wherein I is 1 The stress value is a dangerous index of the stress value of the measuring point; i 2 The stress increment value is a dangerous index, 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 3 Is a dangerous index of stress acceleration value; k (k) 1 、k 2 、k 3 Respectively is I 1 、I 2 、I 3 Weight ratio of (c).
S6, after dimensionless quantities of dynamic and static indexes of rock burst are calculated, analyzing the dynamic and static indexes after comprehensive risk calculation by using a Bayesian probability model, and comprehensively evaluating the risk by using a combined model; let each dynamic and static index be the independent evaluation model, assume that the kth evaluation model weight is P k And coexisting n mutually independent evaluation models, wherein the posterior probability of the occurrence of the combined model is as follows:
wherein O (C) represents the posterior probability of event C; p (C) represents the prior probability of event C, and is generally set to 1/n;is that the event does not occur; the weight of each evaluation model is taken as probability according to the results obtained in the steps S2, S3 and S5;
and S7, calculating the combined model to obtain a dimensionless quantity P belonging to a closed zone of 0-1, grading the dimensionless quantity P according to a preset rock burst risk grading 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 coal mine rock burst risk dynamic and static coupling evaluation method based on the Bayesian method provided by the application has the advantages that firstly, regional dynamic index parameter data and regional static index parameter data are collected, then, abnormal conversion functions are introduced to convert all dynamic and static indexes to obtain multi-constant abnormal indexes, so that all dynamic and static indexes are constrained in dimensionless quantities of a 0-1 closed zone, then, all dynamic and static indexes are subjected to fusion calculation through a mathematical probability statistical method, the comprehensive risk degree of all index weights and rock burst is obtained, and safety early warning is carried out according to a preset rock burst risk grade division standard. According to the application, based on a Bayesian method, the multi-index integrated rock burst disaster safety monitoring and early warning is established, the data indexes of all monitoring systems are effectively integrated, meanwhile, the dynamic indexes and the static indexes are comprehensively considered, the dynamic weight calculation of all indexes based on a time sequence is realized, the uncertainty, the randomness and the ambiguity in the calculation process are reduced, and the rock burst safety early warning capability is improved, so that the problem that the integration unification effect of all indexes and modules of the existing online monitoring systems is poor 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 rock burst risk classification criteria in step S7 specifically classifies the possibility of the dimensionless number P into four classes of none, weak, medium and strong, and the classification criteria of each class of none, weak, medium and strong are respectively 0.ltoreq.p <0.25, 0.25.ltoreq.p <0.5, 0.5.ltoreq.p <0.75, 0.75.ltoreq.p <1, as shown in table 1 below:
TABLE 1 rock burst hazard classification
The none, weak, medium and strong grades respectively adopt green, yellow, orange and red to perform early warning, so that intelligent grading safety early warning is realized.
Finally, it is noted that the above embodiments are only for illustrating the technical solution of the present application and not for limiting the same, and although the present application has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications and equivalents may be made thereto without departing from the spirit and scope of the technical solution of the present application, which is intended to be covered by the scope of the claims of the present application.
Claims (4)
1. The coal mine rock burst risk dynamic and static coupling evaluation method based on the Bayesian method is characterized by comprising the following steps of:
s1, monitoring received data by a microseismic sensor and a borehole stress sensor which are installed underground, transmitting and storing the data to a monitoring system through a data transmission device, and obtaining dynamic index parameter data with time sequence characteristics;
s2, performing comprehensive index evaluation calculation from two aspects of geological factors and mining technologies, performing possibility 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 possibility index calculation;
s3, carrying out normalization processing on the dynamic index parameter data and the static index parameter data preliminarily by introducing an abnormal conversion function to obtain multi-parameter abnormal indexes uniformly constrained in non-dimensional quantities of 0-1 closed interval;
s4, carrying out abnormal index conversion on index values of each microseismic index obtained through monitoring by the microseismic sensor at different moments according to the following formula:
wherein lambda is ij The value range is 0-1 for the membership degree of the corresponding index in the statistical time window t;
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 risk levels at any moment:
P I =k 1 I 1 +k 2 I 2 +k 3 I 3
wherein I is 1 The stress value is a dangerous index of the stress value of the measuring point; i 2 The stress increment value is a dangerous index, 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 3 Is a dangerous index of stress acceleration value; k (k) 1 、k 2 、k 3 Respectively is I 1 、I 2 、I 3 Weight ratio of (2);
s6, analyzing all dynamic and static indexes after comprehensive risk calculation by using a Bayesian probability model, and comprehensively evaluating the risk by using a combined model; let each dynamic and static index be the independent evaluation model, assume that the kth evaluation model weight is P k And coexisting n mutually independent evaluation models, wherein the posterior probability of the occurrence of the combined model is as follows:
wherein O (C) represents the posterior probability of event C; p (C) represents the prior probability of the event C and is set to be 1/n;is that the event does not occur; the weight of each evaluation model is taken as probability according to the results obtained in the steps S2, S3 and S5;
and S7, calculating the combined model to obtain a dimensionless quantity P belonging to a closed zone of 0-1, grading the dimensionless quantity P according to a preset rock burst risk grading standard to obtain rock burst grades at all times, and realizing intelligent grading early warning according to the rock burst grades.
2. The method for evaluating the risk dynamic and static coupling of the rock burst of the coal mine based on the Bayesian method as claimed in claim 1, wherein the microseismic indexes obtained by monitoring the microseismic sensor in the step S4 comprise microseismic intensity factors, microseismic class parameters, microseismic time sequence factors, microseismic b values and microseismic A (b) values.
3. The bayesian-based coal mine rock burst risk dynamic-static coupling evaluation method according to claim 1, wherein the parameter λ in the step S4 is ij The calculation method of (2) is as follows:
for the forward anomaly index:
λ ij =(S ij -S min )/(S max -S min )
for negative anomaly index:
λ ij =(S max -S ij )/(S max -S min )
wherein S is ij The value of the microseismic index at a certain moment in a time window is S max For the maximum value of the values of each microseismic index at each moment in the time window, S min And 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 and static coupling evaluation method according to claim 1, wherein the preset rock burst risk classification standard in the step S7 specifically classifies the possibility of the dimensionless quantity P into four classes of none, weak, medium and strong, and the classification standard values of the respective classes of none, weak, medium and strong 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 classes of none, weak, medium and strong are respectively pre-warned by green, yellow, orange and red.
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 CN113673119A (en) | 2021-11-19 |
CN113673119B true 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) |
Families Citing this family (4)
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 |
CN114352355B (en) * | 2022-03-15 | 2022-06-24 | 中国矿业大学(北京) | Deep mine rock burst control method |
CN115511379B (en) * | 2022-10-28 | 2023-03-24 | 北京科技大学 | Dynamic dividing method and device for rock burst dangerous area |
CN115860582B (en) * | 2023-02-28 | 2023-05-26 | 山东科技大学 | Impact risk intelligent early warning method based on self-adaptive lifting algorithm |
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 (1)
Title |
---|
露天矿边坡稳定性智能评价研究现状、存在问题及对策;杨天鸿;王赫;董鑫;刘飞跃;张鹏海;邓文学;;煤炭学报(第06期);375-393 * |
Also Published As
Publication number | Publication date |
---|---|
ZA202204188B (en) | 2022-06-29 |
CN113673119A (en) | 2021-11-19 |
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) | A kind of comprehensive sequential method for early warning of microseism multidimensional information of impulsion pressure | |
CN108388921A (en) | A kind of overflow leakage real-time identification method based on random forest | |
CN110132218B (en) | Multi-level early warning determination method based on slope stability monitoring | |
CN104360415A (en) | Method for recognizing tight sandstone reservoir cracks | |
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 | |
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. | |
Kabiesz et al. | Application of rule-based models for seismic hazard prediction in coal mines. | |
CN103615236A (en) | Method for real-time monitoring of formation pressure by means of remote mud logging information | |
Liang et al. | Prediction model for coal-gas outburst using the genetic projection pursuit method | |
CN112907130A (en) | Construction risk grade determination method, device and medium based on shield construction parameters | |
CN110837701A (en) | Mining area bottom plate fault water inrush quantitative evaluation method based on full-coupling analysis | |
CN107942383A (en) | Roof sandstone watery grade prediction technique | |
CN110555277A (en) | method for evaluating sand burst risk of mining water under loose aquifer | |
CN109325470A (en) | Working face in the pit homework type intelligent identification Method based on gas density parameter | |
CN110043317B (en) | Mine disaster multi-parameter local danger area judgment and early warning method | |
CN107422373A (en) | A kind of coal mine rock burst disaster on-line early warning method of feature based drift | |
Qiu et al. | Rockburst prediction based on distance discrimination method and optimization technology-based weight calculation method | |
CN112241835B (en) | Multi-source information evaluation method for flood disaster in deep shaft engineering | |
CN113685166B (en) | Drilling accident early warning method and system | |
CN104268785A (en) | Analysis method for deep well floor water inrush risk | |
CN112598061B (en) | Clustering and grading method for tunnel surrounding rock | |
CN108319915B (en) | Multi-time-window simplified form identification method for dynamically adjusting rock burst signal threshold |
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 |