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
mine
main control
rock burst
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.)
Huating Coal Industry Group Co ltd
China University of Mining and Technology CUMT
Original Assignee
Huating Coal Industry Group Co ltd
China University of Mining and Technology CUMT
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 Huating Coal Industry Group Co ltd, China University of Mining and Technology CUMT filed Critical Huating Coal Industry Group Co ltd
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.一种基于矿震统计特征的冲击地压主控因素确定方法，其特征在于：包括以下步骤：1. a method for determining the main control factor of ground pressure based on the statistical characteristics of mine earthquake, it is characterized in that: comprise the following steps: 步骤一、诱冲要素筛选及分类Step 1. Screening and classification of decoy elements 对目标块段的诱冲要素进行筛选，其中，诱冲要素包括煤层埋深、煤层属性、地质构造和开采应力集中区四类要素，之后，将诱冲要素分为第一类要素和第二类要素；The inducing factors of the target block are screened. Among them, the inducing factors include four types of factors: the depth of the coal seam, the properties of the coal seam, the geological structure and the mining stress concentration area. class element; 步骤二、目标矿震群筛选Step 2. Screening of the target mine earthquake swarm 对步骤一中目标块段进行预处理，即筛选目标区域及周围100m范围内的矿震，其中矿震的能级大于或等于104J；Preprocess the target block in step 1, that is, screen the mine shock within the target area and the surrounding 100m range, where the energy level of the mine shock is greater than or equal to 10 4 J; 步骤三、块段范围界定Step 3. Block Scope Definition 首先定义诱冲要素的影响区间，然后界定目标块段，即若多个第二类要素影响区间不重合，则将多个区间分别划分为单一要素块段，若第二类要素影响区间相互重合，则将重复的部分单独划分为一个块段，而不受第二类要素影响的部分划分为无要素块段；First define the influence interval of the decoy element, and then define the target segment, that is, if the influence intervals of multiple second-type elements do not overlap, then divide the multiple intervals into a single element block respectively. , the repeated part is divided into a block segment, and the part that is not affected by the second type of elements is divided into a non-element block segment; 步骤四、相关性指数计算Step 4. Calculation of correlation index 对于第一类要素，定义其中一个块段内任意要素f的变化范围为a～b，将矿震数据平均划分为四个区间，该四个区间内大能量矿震频次与块段总频次的比值分别计为cf1、cf2、cf3、cf4，对于第二类要素，定义其中一个块段内任意要素s，垂直距离的变化范围为d～e，将矿震数据平均划分为四个区间，该四个区间大能量矿震频次与块段总频次的比值分别计为gs1、gs2、gs3、gs4，之后定义矿震相关性指数k，并通过第一公式对其进行计算；For the first type of elements, define the variation range of any element f in one of the blocks as a ~ b, and divide the mine shock data into four intervals. The ratios are calculated as c f1 , c f2 , c f3 , and c f4 respectively. For the second type of elements, define any element s in one of the blocks, and the variation range of the vertical distance is d～e. Each interval, the ratio of the high-energy mine shock frequency to the total block frequency in the four intervals is calculated as g s1 , g s2 , g s3 , g s4 , and then the mine shock correlation index k is defined, and the first formula is used to calculate the correlation index k. Calculation; 步骤五、影响指数确定Step 5. Determine the impact index 首先，定义两个诱冲要素相关性对比指数Pmn，根据第二公式对指数Pmn进行计算，根据计算得到的Pmn值所在区间，取对应的Cmn值并作为构建矩阵的元素，构建r行t列的分析矩阵，之后再计算各块段诱冲要素影响指数S；First, define the correlation index P mn of two hedging factors, and calculate the index P mn according to the second formula. According to the interval of the calculated P mn value, take the corresponding C mn value and use it as the element of the construction matrix to construct The analysis matrix with r rows and t columns, and then calculate the impact index S of the decoy elements of each block; 步骤六、冲击地压主控因素的确定Step 6. Determination of main control factors of rock burst 将步骤五中诱冲要素影响指数S的计算结果作为构建矩阵的元素，构建各块段影响指数矩阵，之后对数据进行标准化处理，处理完成后，构建标准化处理后的矩阵，并计算诱冲要素的权重Wj，其中，Wj值越大，诱冲要素对冲击地压的影响程度越大，Wj值越小，诱冲要素对冲击地压的影响程度越小，则权重Wj最大的诱冲要素即为冲击地压主控因素。Using the calculation result of the impact index S of the decoy factor in step 5 as the element of the construction matrix, construct the impact index matrix of each block segment, and then standardize the data. After the processing is completed, construct the standardized matrix and calculate the decoy factor. The weight W j , among which, the larger the value of W j , the greater the influence of the hedging element on the rockburst, the smaller the value of W j , the smaller the influence of the hedging element on the rockburst, the greater the weight W j The lubricating factor is the main controlling factor of rockburst. 2.根据权利要求1所述的一种基于矿震统计特征的冲击地压主控因素确定方法，其特征在于：所述步骤一中，目标块段的诱冲要素选取煤层埋深、煤层属性、地质构造以及开采应力集中区中的一种要素或多种要素。2 . The method for determining the main control factor of rockburst based on statistical characteristics of mine earthquakes according to claim 1 , wherein in the step 1, the inducing factors of the target block are selected from coal seam burial depth and coal seam attributes. 3 . , one or more elements of geological structure and mining stress concentration areas. 3.根据权利要求1所述的一种基于矿震统计特征的冲击地压主控因素确定方法，其特征在于：所述步骤三中，定义诱冲要素的影响区间的方式为在目标区域内，距离第二类要素走向、倾向垂直距离100m的矩形范围内定义为一个单一要素影响区间。3 . The method for determining the main control factor of rock burst based on statistical characteristics of mine earthquakes according to claim 1 , wherein in the step 3, the method of defining the influence interval of the bouncing element is within the target area. 4 . , and a rectangular area with a vertical distance of 100m from the direction and inclination of the second type of elements is defined as a single element influence interval. 4.根据权利要求1所述的一种基于矿震统计特征的冲击地压主控因素确定方法，其特征在于：所述步骤三中，在确定冲击地压主控因素前，将包含不同诱冲要素的块段进行界定。4. A method for determining the main control factor of rock burst based on statistical characteristics of mine earthquakes according to claim 1, characterized in that: in the step 3, before determining the main control factor of rock burst, different inducements will be included. The block segment of the flush element is defined. 5.根据权利要求1所述的一种基于矿震统计特征的冲击地压主控因素确定方法，其特征在于：所述步骤四中，第一公式表达式为：5. The method for determining the main control factor of ground burst based on statistical characteristics of mine earthquakes according to claim 1, wherein: in the step 4, the first formula is:
Figure FDA0003360134820000031
Figure FDA0003360134820000031
 式中，kf为第一类要素与矿震相关性指数，ks为第二类要素与矿震相关性指数。In the formula, k f is the correlation index between the first type of element and mine earthquake, and k s is the correlation index between the second type of element and mine earthquake. 6.根据权利要求1所述的一种基于矿震统计特征的冲击地压主控因素确定方法，其特征在于：所述步骤五中，第二公式表达式为：6. The method for determining the main control factor of rock burst based on statistical characteristics of mine earthquakes according to claim 1, wherein in the step 5, the second formula expression is: Pmn＝km/kn P mn = km /k n 式中，km为任意诱冲要素m的矿震相关性指数，kn为任意诱冲要素n的矿震相关性指数。In the formula, km is the mine-shock correlation index of any decoy element m , and k n is the mine-shock correlation index of any decoy element n. 7.根据权利要求1所述的一种基于矿震统计特征的冲击地压主控因素确定方法，其特征在于：所述步骤六中，在各块段影响指数矩阵内，当矩阵中任意元素没有对应的影响指数计算值时，则该元素取值为0。7. The method for determining the main control factor of rock burst based on statistical characteristics of mine earthquakes according to claim 1, wherein in the step 6, in the influence index matrix of each block, when any element in the matrix When there is no corresponding calculated value of the influence index, this element takes the value of 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 吉林建筑大学 A method and device for predicting large-energy microseismic events in coal mining

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 吉林建筑大学 A method and device for predicting large-energy microseismic events in coal mining

Also Published As

Publication number Publication date
CN114063152B (en) 2024-01-12

Similar Documents

Publication Publication Date Title
CN103291364B (en) Microseismic multidimensional information comprehensive time sequence early warning method for rock burst
CN104267435B (en) Logging calculation method and device for elastic coefficient of transverse isotropic stratum
CN113222347B (en) Method for evaluating gray system of safety risk of open-air blasting
CN108594295B (en) Method for evaluating coal bed blasting pressure relief effect
CN112377257B (en) A method for determining the influence range of mining advance in working face based on microseismic monitoring
CN115929304B (en) Method for preventing impact of artificial relief layer of stope face
CN111222254A (en) Method and system for classifying rockburst hazard level of working face based on stress superposition method
CN114087021A (en) Rock burst multi-parameter dynamic trend early warning method
CN114135288A (en) An optimization method for pressure relief parameters of high-pressure water jet slits in rock burst coal seam roadways
CN114063152A (en) Rock burst main control factor determination method based on mine earthquake statistical characteristics
CN116485572A (en) A Method for Key Layer Identification and Trend Analysis of Mining Earthquakes Based on Inversion of Seismic Source Parameters
CN109558976A (en) Method for identifying dangerousness of rock burst based on multi-dimensional information
CN108008464A (en) Crack anisotropism quantitatively characterizing method and its system
CN111458742B (en) A method for excitation of explosive source single well dose combination
CN113914932B (en) A method for identifying dangerous areas of coal and gas outburst by using shock wave tomography
Hudaverdi et al. An alternative approach to predict human response to blast induced ground vibration
CN116822418B (en) Mining earthquake mine ground hydraulic fracturing construction horizon identification method and system
CN118346365A (en) A three-dimensional representation method for tunnel rockburst warning information
CN115146971A (en) Rock burst mine distinguishing and grading method based on equivalent depth
CN113434819B (en) Method for Determining the Time and Distance of the Impact of Working Face Mining on Mining Seismic Activities in Goaf
CN115099584A (en) An Evaluation Method for Roof Blasting Pressure Relief Effect Based on Analytic Hierarchy Process
CN116577820B (en) A method for identifying dangerous mine earthquakes in coal mines
CN112241835A (en) Multi-source information evaluation method for water inrush disaster in deep shaft engineering
CN116359013B (en) Tunnel surrounding rock mechanical parameter value-taking method based on multi-data source analysis
CN108008467A (en) Fractue spacing quantitatively characterizing method and its system

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