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 PDF

Info

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
Application number
CN202111364520.7A
Other languages
Chinese (zh)
Other versions
CN114063152B (en
Inventor
田文辉
曹安业
闫林
王崧玮
金建成
上官科峰
薛再君
高利军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology CUMT
Huating Coal Group Co Ltd
Original Assignee
China University of Mining and Technology CUMT
Huating Coal Group Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology CUMT, Huating Coal Group Co Ltd filed Critical China University of Mining and Technology CUMT
Priority to CN202111364520.7A priority Critical patent/CN114063152B/en
Publication of CN114063152A publication Critical patent/CN114063152A/en
Application granted granted Critical
Publication of CN114063152B publication Critical patent/CN114063152B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Geophysics (AREA)
  • Acoustics & Sound (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Geology (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Algebra (AREA)
  • Emergency Management (AREA)
  • Software Systems (AREA)
  • Business, Economics & Management (AREA)
  • Computing Systems (AREA)
  • Computational Linguistics (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

Rock burst main control factor determination method based on mine earthquake statistical characteristics
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:
Figure BDA0003360134830000041
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
Figure BDA0003360134830000071
Figure BDA0003360134830000072
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
Figure BDA0003360134830000073
Figure BDA0003360134830000074
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:
Figure BDA0003360134830000081
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:
Figure BDA0003360134830000082
Pmnvalue and CmnThe corresponding relationship of the values is as follows:
Figure BDA0003360134830000083
Figure BDA0003360134830000091
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:
Figure BDA0003360134830000092
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:
Figure BDA0003360134830000101
in the fifth step, the data is standardized, and the calculation formula is as follows:
and aligning the related elements:
Figure BDA0003360134830000102
for negative correlation elements:
Figure BDA0003360134830000103
in the formula
Figure BDA0003360134830000104
Is 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:
Figure BDA0003360134830000105
calculating the weight W of the inducing elementsjThe calculation formula is as follows:
Figure BDA0003360134830000111
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:
Figure BDA0003360134830000112
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:
Figure BDA0003360134830000113
Figure BDA0003360134830000121
table three each block k value calculation result
Then, the calculated impact factor influence index S of each block is shown in table four:
Figure BDA0003360134830000122
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:
Figure BDA0003360134830000123
table five influence index matrix
Wherein the sixth table is:
Figure BDA0003360134830000124
Figure BDA0003360134830000131
matrix after six standardization processes of table
Recalculating the weight W of the inducing elementsjAs shown in table seven:
Figure BDA0003360134830000132
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:
Figure FDA0003360134820000031
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.
CN202111364520.7A 2021-11-17 2021-11-17 Rock burst main control factor determination method based on mine earthquake statistical characteristics Active CN114063152B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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) A kind of well logging computational methods of transverse isotropy stratum coefficient of elasticity and device
CN113958366B (en) Dynamic quantitative early warning method for impact risk based on vibration-stress double-field monitoring
CN110595887B (en) Quantitative evaluation method for cave forming conditions of shallow hard rock large-scale underground cavern
CN104653226A (en) Stress-gradient-based method for dividing coal impact ground pressure danger area
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
CN112213767A (en) Method for evaluating advanced presplitting blasting effect of top plate
CN112343602A (en) Construction method for upper and lower layer overlapping tunnel
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
CN115392061B (en) Method for evaluating static and dynamic coupling of rock burst danger
CN111046595A (en) Typical and atypical rock burst mine type dividing method
CN114135288B (en) Method for optimizing high-pressure water jet slotting pressure relief parameters of rock burst coal seam roadway
CN116070907A (en) Karst collapse susceptibility assessment method and system based on analytic hierarchy process
CN115983639A (en) Evaluation method for water inrush risk of karst aquifer bottom plate
CN113434819A (en) Method for determining influence time and distance of working face mining on mine earthquake activity of goaf
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

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