CN114063152A - Rock burst main control factor determination method based on mine earthquake statistical characteristics - Google Patents
Rock burst main control factor determination method based on mine earthquake statistical characteristics Download PDFInfo
- Publication number
- CN114063152A CN114063152A CN202111364520.7A CN202111364520A CN114063152A CN 114063152 A CN114063152 A CN 114063152A CN 202111364520 A CN202111364520 A CN 202111364520A CN 114063152 A CN114063152 A CN 114063152A
- Authority
- CN
- China
- Prior art keywords
- earthquake
- rock burst
- elements
- influence
- block
- 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
- 238000000034 method Methods 0.000 title claims abstract description 37
- 239000011159 matrix material Substances 0.000 claims abstract description 36
- 238000005065 mining Methods 0.000 claims abstract description 16
- 238000012216 screening Methods 0.000 claims abstract description 16
- 230000001939 inductive effect Effects 0.000 claims description 40
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 22
- 239000011707 mineral Substances 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 16
- 239000003245 coal Substances 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 7
- 230000035939 shock Effects 0.000 claims description 6
- 238000009933 burial Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 4
- 239000000411 inducer Substances 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 5
- 238000012544 monitoring process Methods 0.000 description 9
- 230000006872 improvement Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 2
- 230000003116 impacting effect Effects 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 230000009528 severe injury Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/288—Event detection in seismic signals, e.g. microseismics
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/903—Querying
- G06F16/9035—Filtering based on additional data, e.g. user or group profiles
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/906—Clustering; Classification
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- Environmental & Geological Engineering (AREA)
- General Engineering & Computer Science (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Geophysics (AREA)
- Mathematical Analysis (AREA)
- Acoustics & Sound (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Linguistics (AREA)
- Computing Systems (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a rock burst main control factor determining method based on ore earthquake statistical characteristics, which comprises the steps of screening and classifying induced impact factors, screening target ore earthquake groups, defining block range, calculating correlation indexes, determining influence indexes and determining rock burst main control factors, wherein through the correlation between high-energy ore earthquakes and impact risk degrees, different block induced impact factor influence indexes are calculated by using an analytic hierarchy process, an influence index matrix is constructed, the weight of each induced impact factor is calculated by using an entropy method, the influence degree of each induced impact factor on rock burst is determined, and the rock burst main control factors are finally determined. The complex mining conditions on site are fully considered, and the accuracy and the credibility of the rock burst master control factor determination are improved.
Description
Technical Field
The invention relates to the technical field of rock burst factor determination, in particular to a rock burst main control factor determination method based on mine earthquake statistical characteristics.
Background
The rock burst is a dynamic phenomenon of severe damage of surrounding rocks, the safety production of a mine is seriously influenced, the rock burst monitoring and early warning method can provide various data related to the occurrence of the rock burst, and the conventional monitoring and early warning means at present comprise a drilling cutting method, a stress monitoring method, an electromagnetic radiation method, a ground sound monitoring method, a micro-seismic monitoring method and the like, wherein the micro-seismic monitoring method can analyze and determine the direction of vibration propagation by recording vibration energy released by coal-rock body fracture induced in the excavation process, accurately position a seismic source and calculate the energy, and has wide field application and better monitoring effect. The mine earthquake data monitored by the method can effectively reflect the stress concentration and the impact danger degree of the mining area;
the main control factor is used as an important index of rock burst and is a basis for revealing, monitoring, early warning and pressure relief measures of a coal mine working face impact mechanism, and the existing method for determining the main control factor of the rock burst mainly comprises methods such as qualitative analysis, numerical simulation and the like, but the methods are theoretical analysis neglecting on-site complex mining conditions, are widely applied to all mine stope working faces, and lack of objective monitoring basis and calculation method, so that the invention provides the method for determining the main control factor of the rock burst based on the mine earthquake statistical characteristics to solve the problems in the prior art.
Disclosure of Invention
In view of the above problems, the present invention aims to provide a method for determining a main factor of rock burst based on statistical characteristics of mineral earthquakes, which specifically determines the main factor of rock burst according to classification of induced impact factors, screening of target mineral earthquakes, and block range definition, as well as correlation indexes and influence indexes of induced impact factors and mineral earthquakes of different block induced impact factors, and calculation of weights of the induced impact factors.
In order to realize the purpose of the invention, the invention is realized by the following technical scheme: a rock burst main control factor determining method based on mine earthquake statistical characteristics comprises the following steps:
step one, screening and classifying the inducing Chongsu
Screening the inducing elements of the target block section, wherein the inducing elements comprise four types of elements of coal bed burial depth, coal bed property, geological structure and mining stress concentration area, and then dividing the inducing elements into a first type of elements and a second type of elements;
step two, screening target ore seismic groups
Preprocessing the target block segment in the step one, namely screening the ore earthquake in the target area and the surrounding 100m range, wherein the energy level of the ore earthquake is more than or equal to 104J;
Step three, defining block segment range
Firstly, defining an influence interval of the inductive elements, and then defining a target block segment, namely if a plurality of second-type element influence intervals are not overlapped, respectively dividing the plurality of intervals into single element block segments, if the second-type element influence intervals are overlapped, independently dividing the repeated part into one block segment, and dividing the part which is not influenced by the second-type element into non-element block segments;
step four, calculating the correlation index
For the first class of elements, the variation range of any element f in one block section is defined as a-b, the mineral earthquake data is averagely divided into four intervals, and the ratio of the frequency of the high-energy mineral earthquake in the four intervals to the total frequency of the block section is respectively counted as cf1、cf2、cf3、cf4For the second kind of elements, defining any element s in one block segment, wherein the variation range of the vertical distance is d-e, averagely dividing the mineral earthquake data into four intervals, and respectively calculating the ratio of the frequency of the high-energy mineral earthquake in the four intervals to the total frequency of the block segment as gs1、gs2、gs3、gs4And then defining a mine earthquake correlation index k and calculating the mine earthquake correlation index k through a first formula.
Step five, determining the influence index
First, two inducer correlation contrast indexes P are definedmnAccording to a second formula, to the index PmnPerforming calculation according to the calculated PmnTaking the corresponding C in the interval of the valuemnThe values are used as elements for constructing a matrix, an analysis matrix of r rows and t columns is constructed, and then the influence index S of each block induced impact element is calculated;
step six, determining the main control factors of rock burst
Taking the calculation result of the impact inducing element influence index S in the step five as an element of a construction matrix, constructing an influence index matrix of each block, then carrying out standardization processing on data, constructing a matrix after the standardization processing after the processing is finished, and calculating the weight W of the impact inducing elementjWherein W isjThe larger the value, the larger the influence degree of the shock element on rock burst, WjThe smaller the value is, the smaller the influence degree of the shock element on the rock burst is, and the weight W isjThe largest inducing element is the main control factor of rock burst.
The further improvement lies in that: in the first step, one or more elements of the coal seam burial depth, the coal seam property, the geological structure and the mining stress concentration area are selected from the inducing element of the target block section.
The further improvement lies in that: in the third step, the influence interval of the inducing elements is defined in a way that a single element influence interval is defined in a rectangular range which is 100m away from the trend and the inclination vertical distance of the second type of elements in the target area.
The further improvement lies in that: in the third step, before determining the rock burst master control factor, defining the block sections containing different inducing factors.
The further improvement lies in that: in the fourth step, the first formula expression is as follows:
in the formula, kfIs the correlation index, k, between the first kind of elements and the mineral earthquakesAnd the correlation index of the second type of elements and the mineral earthquake.
The further improvement lies in that: in the fifth step, the second formula expression is as follows:
Pmn=km/kn
in the formula, kmIs the mine earthquake correlation index, k, of any induced impact element mnAnd the mining earthquake correlation index of any induced impact element n.
The further improvement lies in that: in the sixth step, in each block influence index matrix, when any element in the matrix has no corresponding influence index calculation value, the value of the element is 0.
The invention has the beneficial effects that: the method for determining the rock burst main control factor based on the ore earthquake statistical characteristics calculates the influence indexes of the induced impact factors of different blocks by using an analytic hierarchy process through the correlation between the high-energy ore earthquake and the impact risk degree, constructs an influence index matrix, calculates the weight of each induced impact factor by using an entropy method, determines the influence degree of each induced impact factor on the rock burst, and finally determines the rock burst main control factor. Meanwhile, the method is used for carrying out regional division on the working face, the complicated mining conditions on site are fully considered, the accuracy and the credibility of the rock burst master control factor determination are improved, and a basis is provided for the rock burst master control factor determination and the mine risk decision under similar conditions.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
Fig. 1 is a schematic flow chart according to a first embodiment of the present invention.
Fig. 2 is a schematic diagram of target mine seismic cluster screening according to a second embodiment of the present invention.
FIG. 3 is a diagram illustrating the definition of the range of target segments according to a second embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc., indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings, and are only for convenience of description and simplicity of description, but do not indicate or imply that the device or element being referred to must have a particular orientation, be constructed and operated in a particular orientation, and thus, should not be construed as limiting the present invention. Furthermore, the terms "first," "second," "third," "fourth," and the like are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.
In the description of the present invention, it should be noted that, unless otherwise explicitly specified or limited, the terms "mounted," "connected," and "connected" are to be construed broadly, e.g., as meaning either a fixed connection, a removable connection, or an integral connection; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.
Example one
As shown in fig. 1, the embodiment provides a method for determining a rock burst main control factor based on a mineral earthquake statistical characteristic, which includes the following steps:
step one, screening and classifying the inducing Chongsu
Screening the inducing elements of the target block section, wherein the inducing elements comprise four types of elements of a coal seam buried depth, a coal seam property, a geological structure and a mining stress concentration area, and then, the inducing elements are divided into a first type of element and a second type of element;
step two, screening target ore seismic groups
Preprocessing the target block segment in the step one, namely screening the ore earthquake in the target area and the surrounding 100m range, wherein the energy level of the ore earthquake is more than or equal to 104J, i.e. having an energy level of 10 or more4J, projecting the mine earthquake on a target area of a plane base map of mining engineering, and drawing a large-energy mine earthquake distribution map;
step three, defining block segment range
The method comprises the following steps of defining block ranges of a target area of a mining engineering plane base map by using straight lines, then defining influence intervals of induced impact elements, and then defining target blocks, namely if the influence intervals of a plurality of second-type elements are not overlapped, respectively dividing the plurality of intervals into single-element blocks, if the influence intervals of the second-type elements are overlapped, independently dividing repeated parts into one block, and dividing parts which are not influenced by the second-type elements into non-element blocks, wherein in the third step, the influence intervals of the induced impact elements are defined in such a way that in the target area, a rectangular range which is 100m away from the trend and the inclination vertical distance of the second-type elements is defined as a single-element influence interval, and in the third step, before determining the main control factor of the impact pressure, the blocks containing different induced impact elements are defined;
step four, calculating the correlation index
Defining the block range of a target area of the mining engineering plane base map by using straight lines, defining the variation range of any element f in one or a certain block as a-b for the first type of elements, and averagely dividing the mineral seismic data into The ratio of the large energy ore earthquake frequency to the total frequency of the block segments in the four intervals is respectively calculated as cf1、cf2、cf3、cf4For the second class of elements, any element s in one or a certain block is defined, the variation range of the vertical distance is d-e, and the mine earthquake data are averagely divided into The ratio of the earthquake frequency of the high-energy mine to the total frequency of the block segments in the four intervals is gs1、gs2、gs3、gs4And then defining a mine earthquake correlation index k and calculating the mine earthquake correlation index k through a first formula, wherein in the fourth step, the first formula expression is as follows:
in the formula, kfIs the correlation index, k, between the first kind of elements and the mineral earthquakesAnd (3) a correlation index of the second type of elements and the mineral earthquake, wherein l is a correction coefficient, when the inducing and impacting elements are geological structures, l is 2, and other inducing and impacting elements l are 1.
Step five, determining the influence index
First, two inducer correlation contrast indexes P are definedmnAccording to a second formula, to the index PmnPerforming calculation according to the calculated PmnTaking the corresponding C in the interval of the valuemnAnd taking the value as an element of a construction matrix, constructing an analysis matrix of r rows and t columns, and then calculating the influence index S of the induced impact element of each block, wherein in the step five, a second formula expression is as follows:
Pmn=km/kn;
in the formula, kmIs the mine earthquake correlation index, k, of any induced impact element mnAnd the mining earthquake correlation index of any induced impact element n.
An analysis matrix of r rows and t columns is constructed as follows:
Pmnvalue and CmnThe corresponding relationship of the values is as follows:
table-matrix table for taking value of any element
Calculating the influence index S of the induced impact factors of each block, wherein the calculation formula is as follows:
in the formula, CmnFor the m-th row and n-th column element of the matrix, UmIs the value after matrix normalization, S is the maximum eigenvalue of the analysis matrix, i.e. the influence index of the induced impact element of each block, t is the order of the analysis matrix, CI and RI are both indexes in the analytic hierarchy process, CR is a verification index, and the requirement that CR is less than 0.1 is met, otherwise C is properly adjustedmnTaking values until CR is less than 0.1.
Step six, determining the main control factors of rock burst
Taking the calculation result of the impact inducing element influence index S in the step five as an element of a construction matrix, constructing an influence index matrix of each block, then carrying out standardization processing on data, constructing a matrix after the standardization processing after the processing is finished, and calculating the weight W of the impact inducing elementjWherein W isjThe larger the value, the larger the influence degree of the shock element on rock burst, WjThe smaller the value is, the smaller the influence degree of the shock element on the rock burst is, and the weight W isjIn the sixth step, in each block influence index matrix, when any element in the matrix has no corresponding influence index calculation value, the value of the element is 0;
assuming that i block segments are total, j inducing elements, and constructing an influence index matrix of each block segment as follows:
in the fifth step, the data is standardized, and the calculation formula is as follows:
in the formulaIs the value after standardized processing of any element, SijTo influence any element of the index matrix, Smax、SminIn order to influence the maximum value and the minimum value of the ith row element of the index matrix, a normalized matrix is constructed:
calculating the weight W of the inducing elementsjThe calculation formula is as follows:
in the formula QijThe proportion of the ith block segment in the inducing element under the jth inducing element, m 'and n' represent the number of the inducing elements and the block segments to be evaluated, HjEntropy of the jth inductive element, WjIs the weight of the jth inductive element.
Example two
According to the embodiment shown in fig. 2-3, the major factor of rock burst in a mining area 2501 is determined by selecting microseismic events monitored during stoping of the working face 250104-1 of the northern coal mine of inkstone;
first, the Chongchongsu was screened and classified, and any Chongchongsu element was replaced with a numeric designation, as shown in Table two:
second-class Chongsu classification table
Then selecting the energy level of the target area to be more than or equal to 104J, mineral earthquake, as shown in figure 2, defining a target block range, as shown in figure 3, and calculating a mineral earthquake correlation index k of each block induced impact element, as shown in table three:
table three each block k value calculation result
Then, the calculated impact factor influence index S of each block is shown in table four:
results of calculation of influence index S of each induced impact element in table four
And then constructing the influence index matrix of the induced impact factors of each block and the result after the standardization treatment, wherein the result is shown in a fifth table and a sixth table, wherein the fifth table is as follows:
table five influence index matrix
Wherein the sixth table is:
matrix after six standardization processes of table
Recalculating the weight W of the inducing elementsjAs shown in table seven:
weights of inducing elements of the seven tables
Finally, according to the weight WjThe largest inducing impact element is the main control factor of rock burst, namely the determinable main control factor of rock burst is the declination axis.
The foregoing illustrates and describes the principles, general features, and advantages of the present invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are described in the specification and illustrated only to illustrate the principle of the present invention, but that various changes and modifications may be made therein without departing from the spirit and scope of the present invention, which fall within the scope of the invention as claimed. The scope of the invention is defined by the appended claims and equivalents thereof.
Claims (7)
1. A rock burst main control factor determining method based on mineral earthquake statistical characteristics is characterized by comprising the following steps: the method comprises the following steps:
step one, screening and classifying the inducing Chongsu
Screening the inducing elements of the target block section, wherein the inducing elements comprise four types of elements of coal bed burial depth, coal bed property, geological structure and mining stress concentration area, and then dividing the inducing elements into a first type of elements and a second type of elements;
step two, screening target ore seismic groups
Preprocessing the target block segment in the step one, namely screening the ore earthquake in the target area and the surrounding 100m range, wherein the energy level of the ore earthquake is more than or equal to 104J;
Step three, defining block segment range
Firstly, defining an influence interval of the inductive elements, and then defining a target block segment, namely if a plurality of second-type element influence intervals are not overlapped, respectively dividing the plurality of intervals into single element block segments, if the second-type element influence intervals are overlapped, independently dividing the repeated part into one block segment, and dividing the part which is not influenced by the second-type element into non-element block segments;
step four, calculating the correlation index
For the first class of elements, the variation range of any element f in one block section is defined as a-b, the mineral earthquake data is averagely divided into four intervals, and the ratio of the frequency of the high-energy mineral earthquake in the four intervals to the total frequency of the block section is respectively counted as cf1、cf2、cf3、cf4For the second kind of elements, defining any element s in one block segment, wherein the variation range of the vertical distance is d-e, averagely dividing the mineral earthquake data into four intervals, and respectively calculating the ratio of the frequency of the high-energy mineral earthquake in the four intervals to the total frequency of the block segment as gs1、gs2、gs3、gs4Then defining a mine earthquake correlation index k, and calculating the mine earthquake correlation index k through a first formula;
step five, determining the influence index
First, two inducer correlation contrast indexes P are definedmnAccording to a second formula, to the index PmnPerforming calculation according to the calculated PmnTaking the corresponding C in the interval of the valuemnThe values are used as elements for constructing a matrix, an analysis matrix of r rows and t columns is constructed, and then the influence index S of each block induced impact element is calculated;
step six, determining the main control factors of rock burst
Taking the calculation result of the impact inducing element influence index S in the step five as an element of a construction matrix, constructing an influence index matrix of each block, then carrying out standardization processing on data, constructing a matrix after the standardization processing after the processing is finished, and calculating the weight W of the impact inducing elementjWherein W isjThe larger the value, the larger the influence degree of the shock element on rock burst, WjThe smaller the value is, the smaller the influence degree of the shock element on the rock burst is, and the weight W isjThe largest inducing element is rock burstA master control factor.
2. The rock burst master factor determination method based on the mine earthquake statistical characteristics according to claim 1, characterized in that: in the first step, one or more elements of the coal seam burial depth, the coal seam property, the geological structure and the mining stress concentration area are selected from the inducing element of the target block section.
3. The rock burst master factor determination method based on the mine earthquake statistical characteristics according to claim 1, characterized in that: in the third step, the influence interval of the inducing elements is defined in a way that a single element influence interval is defined in a rectangular range which is 100m away from the trend and the inclination vertical distance of the second type of elements in the target area.
4. The rock burst master factor determination method based on the mine earthquake statistical characteristics according to claim 1, characterized in that: in the third step, before determining the rock burst master control factor, defining the block sections containing different inducing factors.
5. The rock burst master factor determination method based on the mine earthquake statistical characteristics according to claim 1, characterized in that: in the fourth step, the first formula expression is as follows:
in the formula, kfIs the correlation index, k, between the first kind of elements and the mineral earthquakesAnd the correlation index of the second type of elements and the mineral earthquake.
6. The rock burst master factor determination method based on the mine earthquake statistical characteristics according to claim 1, characterized in that: in the fifth step, the second formula expression is as follows:
Pmn=km/kn
in the formula, kmIs the mine earthquake correlation index, k, of any induced impact element mnAnd the mining earthquake correlation index of any induced impact element n.
7. The rock burst master factor determination method based on the mine earthquake statistical characteristics according to claim 1, characterized in that: in the sixth step, in each block influence index matrix, when any element in the matrix has no corresponding influence index calculation value, the value of the element is 0.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111364520.7A CN114063152B (en) | 2021-11-17 | 2021-11-17 | Rock burst main control factor determination method based on mine earthquake statistical characteristics |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111364520.7A CN114063152B (en) | 2021-11-17 | 2021-11-17 | Rock burst main control factor determination method based on mine earthquake statistical characteristics |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114063152A true CN114063152A (en) | 2022-02-18 |
CN114063152B CN114063152B (en) | 2024-01-12 |
Family
ID=80277547
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111364520.7A Active CN114063152B (en) | 2021-11-17 | 2021-11-17 | Rock burst main control factor determination method based on mine earthquake statistical characteristics |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114063152B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115390130A (en) * | 2022-08-29 | 2022-11-25 | 吉林建筑大学 | Coal mining high-energy microseismic event prediction method and device |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160070828A1 (en) * | 2013-04-08 | 2016-03-10 | China University of Mining & Technology, Beijng | Vulnerability Assessment Method of Water Inrush from Aquifer Underlying Coal Seam |
CN109707453A (en) * | 2018-12-12 | 2019-05-03 | 中煤能源研究院有限责任公司 | A kind of permanent chamber impact Safety Argumentation method of coal mine |
CN111898267A (en) * | 2020-07-30 | 2020-11-06 | 中煤能源研究院有限责任公司 | Method for evaluating dangerousness of rock burst based on sedimentary geological conditions |
CN112594008A (en) * | 2020-12-17 | 2021-04-02 | 中国矿业大学 | Rock burst early warning method based on dynamic and static load disturbance superposition |
-
2021
- 2021-11-17 CN CN202111364520.7A patent/CN114063152B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160070828A1 (en) * | 2013-04-08 | 2016-03-10 | China University of Mining & Technology, Beijng | Vulnerability Assessment Method of Water Inrush from Aquifer Underlying Coal Seam |
CN109707453A (en) * | 2018-12-12 | 2019-05-03 | 中煤能源研究院有限责任公司 | A kind of permanent chamber impact Safety Argumentation method of coal mine |
CN111898267A (en) * | 2020-07-30 | 2020-11-06 | 中煤能源研究院有限责任公司 | Method for evaluating dangerousness of rock burst based on sedimentary geological conditions |
CN112594008A (en) * | 2020-12-17 | 2021-04-02 | 中国矿业大学 | Rock burst early warning method based on dynamic and static load disturbance superposition |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115390130A (en) * | 2022-08-29 | 2022-11-25 | 吉林建筑大学 | Coal mining high-energy microseismic event prediction method and device |
Also Published As
Publication number | Publication date |
---|---|
CN114063152B (en) | 2024-01-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109798106B (en) | Method for predicting risk of rock burst and prevention and treatment measures | |
CN108594295B (en) | Method for evaluating coal bed blasting pressure relief effect | |
CN105741029B (en) | A kind of danger of burst evaluation method based on cumulative stress coefficient of concentration | |
CN104267435B (en) | Logging calculation method and device for elastic coefficient of transverse isotropic stratum | |
CN113958366B (en) | Dynamic quantitative early warning method for impact risk based on vibration-stress double-field monitoring | |
CN113222347B (en) | Method for evaluating gray system of safety risk of open-air blasting | |
CN111222254A (en) | Working face rock burst danger grade dividing method and system based on stress superposition method | |
CN114169789A (en) | Coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment | |
CN114063152A (en) | Rock burst main control factor determination method based on mine earthquake statistical characteristics | |
CN115392061B (en) | Method for evaluating static and dynamic coupling of rock burst danger | |
CN112213767A (en) | Method for evaluating advanced presplitting blasting effect of top plate | |
CN114135288B (en) | Method for optimizing high-pressure water jet slotting pressure relief parameters of rock burst coal seam roadway | |
CN116485572B (en) | Mining earthquake key layer identification and trend analysis method based on earthquake source parameter inversion | |
CN113914932B (en) | Method for identifying coal and gas outburst dangerous area by using vibration wave tomography | |
CN116822418B (en) | Mining earthquake mine ground hydraulic fracturing construction horizon identification method and system | |
CN108008464A (en) | Crack anisotropism quantitatively characterizing method and its system | |
Yunisovna | Formation of the earthquake database for evaluation of their influence on the slope stability of deep quarries | |
CN115146971A (en) | Rock burst mine distinguishing and grading method based on equivalent depth | |
CN111046595A (en) | Typical and atypical rock burst mine type dividing method | |
CN116070907A (en) | Karst collapse susceptibility assessment method and system based on analytic hierarchy process | |
CN113434819A (en) | Method for determining influence time and distance of working face mining on mine earthquake activity of goaf | |
El Sgher | The Impacts of the Net Stress and Stress Shadow on the Productivity of Marcellus Shale Horizontal Well | |
CN113496006B (en) | Method for calculating centralization degree of mine earthquake space of rock burst mine | |
CN116663097B (en) | Method for designing advanced stress release hole parameters of rock burst source region of deep-buried hard rock construction tunnel | |
CN110748341A (en) | Method for detecting influence radius of controllable seismic source based on carbon dioxide blasting |
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 |