CN114169789A - Coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment - Google Patents
Coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment Download PDFInfo
- Publication number
- CN114169789A CN114169789A CN202111523840.2A CN202111523840A CN114169789A CN 114169789 A CN114169789 A CN 114169789A CN 202111523840 A CN202111523840 A CN 202111523840A CN 114169789 A CN114169789 A CN 114169789A
- Authority
- CN
- China
- Prior art keywords
- factor
- rock burst
- coal
- matrix
- judgment
- 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.)
- Withdrawn
Links
- 239000011435 rock Substances 0.000 title claims abstract description 114
- 239000003245 coal Substances 0.000 title claims abstract description 90
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000004458 analytical method Methods 0.000 title claims abstract description 20
- 239000011159 matrix material Substances 0.000 claims abstract description 62
- 238000011156 evaluation Methods 0.000 claims abstract description 53
- 238000005065 mining Methods 0.000 claims abstract description 25
- 239000013598 vector Substances 0.000 claims description 35
- 239000010410 layer Substances 0.000 claims description 25
- 239000011241 protective layer Substances 0.000 claims description 17
- 230000033001 locomotion Effects 0.000 claims description 9
- 238000013461 design Methods 0.000 claims description 6
- BTCSSZJGUNDROE-UHFFFAOYSA-N gamma-aminobutyric acid Chemical compound NCCCC(O)=O BTCSSZJGUNDROE-UHFFFAOYSA-N 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 238000007405 data analysis Methods 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 3
- 238000007619 statistical method Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 238000011161 development Methods 0.000 claims description 2
- 230000003116 impacting effect Effects 0.000 claims 1
- 230000002265 prevention Effects 0.000 abstract description 5
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- 238000007796 conventional method Methods 0.000 abstract description 2
- 238000012216 screening Methods 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 17
- 238000012544 monitoring process Methods 0.000 description 5
- 230000033228 biological regulation Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 230000006698 induction Effects 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- 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/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Data Mining & Analysis (AREA)
- Operations Research (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Analysis (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Probability & Statistics with Applications (AREA)
- Agronomy & Crop Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Marine Sciences & Fisheries (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Mining & Mineral Resources (AREA)
- Algebra (AREA)
- Animal Husbandry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment, relating to the technical field of mine safety production, by analyzing, screening and considering the influence of a plurality of factors of the geological dynamic environment type, the geology type and the mining type on the impact mine pressure danger, compared with other methods only considering single or a plurality of factors, the method has more comprehensive and reasonable analysis, secondly, each factor is endowed with respective corresponding weight according to different influence degrees, a fuzzy relation matrix is formed by combining the membership degrees of the factors, the risk of the rock burst is predicted by adopting a fuzzy comprehensive evaluation method, the calculability and the quantifiability of the risk degree of the rock burst are realized, meanwhile, the dangers of rock burst are divided into different grades, which is more detailed and reasonable than the conventional method of adopting a two-stage system with impact and without impact, so that the rock burst prevention and treatment work is more targeted and effective.
Description
Technical Field
The invention relates to the technical field of mine safety production, in particular to a coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment.
Background
The rock burst is a dynamic phenomenon of sudden and violent damage caused by instantaneous release of elastic energy of coal and rock mass in a coal mining space, and the coal and rock mass damage process is accompanied by the characteristics of vibration, loud sound, air wave and the like, so that the rock burst has strong destructiveness. In recent years, the safety situation of the coal mine in China is improved day by day, but a plurality of large and serious rock burst accidents occur successively. According to incomplete statistics, as of 2019, the rock burst mine in China reaches more than 180 seats, is mainly distributed in 25 provinces such as Shandong, Henan, Heilongjiang, Shaanxi and Shanxi and is almost spread over all major coal mining areas, and the rock burst causes secondary disasters such as large-area tunnel collapse, personal casualties, equipment damage, gas outburst induction and the like, and becomes one of the main disasters seriously restricting the safety production of the mine. At present, the prediction method of rock burst mainly includes two types: a kind of local detection method mainly using drilling cutting method, including coal-rock body deformation observation method, coal-rock body stress measurement method, flowing ground sound detection method, rock cake method, etc., mainly used for surveying the impact danger degree of excavating the local section; the second type is a system monitoring method, which includes a geophone system monitoring method and a microseismic system monitoring method, and other geophysical methods such as electromagnetic radiation, earth temperature, geomagnetism and the like.
The induction mechanism of rock burst is complex, the influence factor is more, and when the impact risk evaluation and prediction is carried out, if the considered influence factor is single and one-sided, the accuracy of the prediction result is influenced. Therefore, it is necessary to comprehensively analyze the influence factors of rock burst systematically, determine the influence level of each factor on rock burst by a scientific and objective method, and comprehensively determine the risk of rock burst by combining the weight of each factor.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment, which comprises the following steps:
s1: on the basis of statistical analysis of the occurrence of rock burst disasters, determining a judgment factor set U of the rock burst dangerousness, wherein U is (U ═ U-11，u12，u13，u14，u21，u22，u23，u24，u25，u31，u32，u33，u34，u35，u36) Obtaining values of all judgment factors according to field test, geological data analysis and mine design related data;
s2: determining a judgment set of rock burst dangerousness according to the apparent intensity of rock burst;
s3: analyzing the relationship between each factor and rock burst development strength, establishing a fuzzy relationship between each evaluation factor and an evaluation grade in an evaluation set, and calculating by combining the magnitude of each factor to obtain a fuzzy relationship matrix;
s4: determining the importance degree of each factor of the rock burst to the previous level by using an analytic hierarchy process, and calculating to obtain the weight of each rock burst influence factor to obtain the weight vector of the rock burst influence factor;
s5: calculating the membership degree of each impact danger grade according to the fuzzy relation matrix and the weight vector by adopting a fuzzy comprehensive evaluation method, and determining the danger grade of the rock burst according to a maximum membership degree principle;
the judging factor set of the rock burst risk in the S1 comprises a geological dynamic environment type factor U1Geological factor U2And mining factor U3；
The geomechanical environment type factor U1Including making a concave angle of reflection u11Ratio u of structural stress to vertical stress12Vertical motion rate u of broken block13And the range of influence u of fracture structure14；
The geological factor U2Including the depth of mining u21Coefficient of variation u of coal seam thickness22Fault dip angle u23Uniaxial compressive strength u of coal24Distance u between thick and hard rock stratum and coal seam25；
The mining factor U3Including a degree of substantial collapse u of overburden31Protective layer pressure relief degree u32The distance u between the coal bed and the upper coal pillar33The relation u of the working face to the adjacent goaf34Section coal pillar width u35Thickness u of coal left at bottom36；
In S2, the evaluation set V for the risk of rock burst according to the rules for preventing and treating coal mine rock burst pressure includes: without risk of impact v1Risk of weak impact v2Moderate impact risk v3And the danger of strong impact v4；
The S3 establishes a fuzzy relationship between each of the evaluation factors and the evaluation level in the evaluation set, including:
s3.1: determining membership degree vectors of each evaluation factor in the evaluation factor set to each evaluation grade in the evaluation set according to the known relationship between each factor and rock burst showing strength;
s3.2: the membership degrees of all factors jointly form a fuzzy relation matrix R;
the S4: according to the category of rock burst influencing factors, including the geological dynamic environment type factor U1Geological factor U2And mining factor U3Establishing a hierarchical structure model for rock burst prediction and providing a multilevel index system for rock burst prediction, which specifically comprises the following steps:
s4.1: establishing a hierarchical structure model for rock burst prediction, wherein a target layer is a rock burst risk U, and a criterion layer is a geological dynamic environment factor U1Geological factor U2And mining factor U3The index layer is constructed with a concave reflection angle u11Ratio u of structural stress to vertical stress12Vertical motion rate u of broken block13And the range of influence u of fracture structure14Mining depth u21Coefficient of variation u of coal seam thickness22Fault, faultInclination u23Uniaxial compressive strength u of coal24Distance u between thick and hard rock stratum and coal seam25Sufficient collapse degree u of overburden stratum31Protective layer pressure relief degree u32The distance u between the coal bed and the upper coal pillar33The relation u of the working face to the adjacent goaf34Section coal pillar width u35Thickness u of coal left at bottom36；
S4.2: and comparing the influence degrees of the indexes of the same layer in the evaluation factor set on the indexes belonging to the upper layer pairwise to form a judgment matrix:
Z＝(zij)n×n
zij＞0
wherein z isijZ is the importance of the factor i to the factor j relative to the upper level indexjiThe importance of the factor j to the upper-level index is greater than that of the factor i, and n is the order number of the judgment matrix;
s4.3: calculating a weight vector of the judgment matrix;
normalizing each column vector of the judgment matrix to obtain:
b is toijSumming by rows:
to wiNormalization treatment:
wherein the content of the first and second substances,
to determine the normalized arithmetic mean of each column vector in matrix a,
judging the weight vector of the matrix for the solved feature vector;
s4.4: calculating the maximum characteristic root of the judgment matrix according to the product ZW of the judgment matrix and the corresponding weight vector:
s4.5: calculating a consistency index CI:
s4.6: setting average random consistency index RI, and calculating consistency ratio of judgment matrix
I.e. the ratio of the average random consistency index to the average random consistency index, sets a threshold C1When the consistency ratio of the matrix is larger than C1If so, judging that the matrix does not meet the consistency check, and returning to S4.1 to reconstruct the judgment matrix; when the consistency ratio of the matrix is less than C1And then, judging that the matrix meets the consistency check, and taking the characteristic vector of the judgment matrix as a weight vector.
Advantageous technical effects
1. According to the method, the influence of a plurality of factors such as geological dynamic environment, geology and mining on the rock burst danger is analyzed, screened and considered, compared with other methods only considering single or a plurality of factors, the analysis is more comprehensive and reasonable, each factor is given corresponding weight according to different influence degrees, the calculability and quantifiability of the rock burst danger degree are realized, the rock burst danger is divided into different grades, the rock burst danger is more detailed and reasonable compared with the conventional method adopting a two-level system with impact and without impact, and the rock burst prevention and control work is more targeted and effective;
2. the method introduces the geological dynamic environment factors into the rock burst risk evaluation system, evaluates the rock burst from the angles of macroscopic crustal motion, regional structure and the like, and provides a new research idea for evaluating the rock burst;
3. the rock burst evaluation method overcomes the problem that the weight of each influence factor is not considered in the existing evaluation method, so that a rock burst evaluation system is more scientific and reasonable, the method is simple to operate and high in feasibility, and better theoretical support can be provided for the prediction and prevention work of the rock burst.
Drawings
Fig. 1 is a flowchart of a coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive evaluation according to an embodiment of the present invention;
fig. 2 is a model diagram of a hierarchical structure of rock burst prediction according to an embodiment of the present invention.
Detailed Description
The following describes the embodiments of the present invention in further detail with reference to the drawings and examples;
the invention provides a coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment, which comprises the following steps as shown in figure 1:
s1: on the basis of statistical analysis of the occurrence of rock burst disasters, the judgment factor set U for determining the rock burst danger is (U)11， u12，u13，u14，u21，u22，u23，u24，u25，u31，u32，u33，u34，u35，u36) Obtaining values of all judgment factors according to field test, geological data analysis and mine design related data;
according to the judgment factor set influencing the occurrence of rock burst, the factors comprise the geomechanical environment type factor U1Geological factor U2And openSampling factor U3(ii) a Wherein, the geomechanical environment factors U1Including making a concave angle of reflection u11Ratio u of structural stress to vertical stress12Vertical motion rate u of broken block13And the range of influence u of fracture structure14Geological factor U2Including the depth of mining u21Coefficient of variation u of coal seam thickness22Fault drop u23Uniaxial compressive strength u of coal24Distance u between thick and hard rock stratum and coal seam25(ii) a Mining factor U3Including a degree of substantial collapse u of overburden31Protective layer pressure relief degree u32The distance u between the coal bed and the upper coal pillar33The relation u of the working face to the adjacent goaf34Section coal pillar width u35Thickness u of coal left at bottom36；
Meanwhile, values of all judgment factors are obtained according to field test, geological data analysis and mine design related data, and are shown in a table 1;
TABLE 1 evaluation of rock burst of a certain mine
S2: determining an evaluation set of the rock burst dangerousness, namely an evaluation grade set of the rock burst dangerousness;
based on the requirement that the rock burst classification tends to be objective and the judgment standard of multi-index comprehensive judgment, determining the judgment set V-V according to the regulation of the fine rule for preventing and treating coal mine rock burst1，v2，v3，v4The method comprises the following steps: without risk of impact v1Risk of weak impact v2Moderate impact risk v3High impact risk v4See table 2;
TABLE 2 evaluation of rock burst hazard for various influencing factors
Wherein b is the active fracture influence range, and the thick hard rock stratum refers to a rock stratum with uniaxial compressive strength of more than 60MPa and thickness of more than 10 m; the vertical movement rate of the fault block can be read by a Chinese continental structure environment monitoring network, the ratio of the structural stress to the vertical stress can be obtained by mine ground stress measurement, and the mining depth, the coal seam thickness variation coefficient, the coal uniaxial compressive strength, the fault fall and the distance between a thick hard rock stratum and a coal seam can be obtained by mine geological data such as a drilling column diagram, a geological structure outline drawing and the like of a well field; the distance between the coal bed and the upper coal pillar, the relation between the working face and the adjacent goaf, the width of the coal pillar of the section and the thickness of the coal with remained bottom can be comprehensively determined according to the mining engineering plan, the mine design specification, the working face design specification and the actual mining conditions of the mine, and the process for determining the indexes of the structural concave contrast angle, the active fracture influence range, the coal bed thickness variation coefficient, the sufficient caving degree of the overlying rock stratum and the pressure relief degree of the protective layer is as follows:
(1) texture concave contrast angle
One of the remarkable features of the construction of a valley is the relatively high elevation of the raised areas on both sides of the valley; evaluating the geological dynamic environment of the constructed concave land by utilizing the concave land contrast angle; the following calculation formula is adopted:
in the formula, beta is a reflection angle of a concave ground;
Δ h-difference between highest and lowest elevation of the formation pit, km;
Δ l — the distance between the center of the formation pit and the boundary, km;
in general, the smaller the structure contrast angle, the less the risk of impact ground pressure occurring in the concave structure;
(2) range of effect of active fracture
When the straight-line distance between the well field boundary and the active fracture is smaller than the influence range b of the fracture structure, the fracture structure belongs to the dangerous area of rock burst, and the formula (2) shows that:
b＝K·10h (2)
wherein, the K-activity coefficient (K ═ 1, 2, 3), when the fracture activity is strong, K ═ 3, when the fracture activity is medium, K ═ 2, when the fracture activity is weak, K ═ 1; h-fracture vertical drop;
according to the regulation of geotechnical engineering investigation standard, activities are available since the middle and late renewals, the activities of a new era are strong, the average fracture activity rate v is more than 1mm/a, the seismic magnitude M of the historical earthquake is more than or equal to 7, and the fracture belongs to strong activity fracture; the activity of the middle and late updated generations is stronger in the whole new generation, v is more than or equal to 0.1mm/a and less than or equal to 1mm/a, and M is more than or equal to 5 and less than 7, the fracture belongs to the middle activity fracture; the activity of the old and the new generations is stronger, when v is less than 0.1mm/a, M is less than 5, and the activity belongs to weak activity fracture; therefore, under the condition that the vertical fall h of the fracture structure is the same, the fracture structure activity is stronger, the fracture activity K value is larger, the fracture influence range b value is larger, and when the linear distance b between the well field boundary and the active fracture is smaller than or equal to L, the geological dynamic environment of the well field is influenced;
(3) coefficient of variation of coal seam thickness
The coefficient of variation r for coal thickness was calculated as follows:
in the formula, n is the total number of coal points participating in evaluation;
xi-measured coal thickness, m, for each coal-seeing point;
the larger the variation coefficient of the coal seam thickness is, the more unstable the coal seam occurrence is, and the higher the impact risk degree is;
(4) degree of full collapse of overburden
Observing whether a main key layer of the overburden rock is broken, wherein the specific process for judging the key layer is as follows: (1) determining the position of a hard rock layer in the overlying strata from bottom to top, (2) calculating the breaking distance of each hard rock layer, and (3) comparing the breaking distances of each hard rock layer to determine the position of a key layer;
(5) degree of pressure relief of protective layer
According to the national standard' rock burst determination, monitoring and prevention method part 12: a protective layer mining prevention and control method "; when the evaluation area is in the effective pressure relief range and within the effective pressure relief period, the pressure relief degree of the protective layer is 'good'; when the pressure is in an effective pressure relief range and exceeds the effective pressure relief time, the pressure relief degree of the protective layer is 'medium'; the protective layer is not in the effective pressure relief range of the protective layer, and the pressure relief degree of the protective layer is 'normal'; a bearing coal pillar is left when the protective layer is mined, the plane projection of the coal pillar is positioned in an evaluation area, and the pressure relief degree of the protective layer is poor;
s3: establishing a fuzzy relation between each factor in the judging factors and the judging level in the judging set to obtain a fuzzy relation matrix;
for the ith factor u in the factor setiPerforming judgment to judge the jth element v of the setjDegree of membership of rijFactor uiThe result of the evaluation can be represented as Ri＝(ri1，ri2，ri3，ri4)；
The degree of membership may be determined by establishing a membership function for each factor, the selected constructional relief angle u being concave11Ratio u of structural stress to vertical stress12And the range of influence u of fracture structure14Mining depth u21Coefficient of variation u of coal seam thickness22Fault drop u23Uniaxial compressive strength u of coal24Distance u between thick and hard rock stratum and coal seam25The distance u between the coal bed and the upper coal pillar33Thickness u of coal left at bottom36The membership functions are descending (ascending) half trapezoid and trapezoid distribution on the basis of analysis of various factors, and are shown as formulas (4) to (7):
selected vertical motion rate u of a broken block13Sufficient collapse degree u of overburden stratum31Protective layer pressure relief degree u32The relation u of the working face to the adjacent goaf34Section coal pillar width u35The membership functions of (a) are rectangular distributions on the basis of analysis of various factors;
wherein r is1(x)，r2(x)，r3(x)，r4(x) Are respectively a set of judgments v1，v2，v3，v4Membership function of each evaluation factor, x is the value of each factor, x0，x1，x2Are respectively due toCritical value of the evaluation grade of the expression intensity of the mineral pressure;
s3.2: the membership degree of each factor forms a fuzzy relation matrix R of 17 multiplied by 4 order, and the fuzzy relation matrix R is calculated according to the membership degree function
S4: determining the importance degree of each factor to the previous level by using an analytic hierarchy process, and calculating the weight of each influencing factor to obtain the weight vector of each factor; as shown in fig. 2, the method specifically includes the following steps:
s4.1: according to the types of rock burst influence factors, a hierarchical structure model for rock burst prediction is established, the target layer is a rock burst risk U, and the standard layer is a geological dynamic environment factor U1Geological factor U2And mining factor U3The index layer is constructed with a concave reflection angle u11Ratio u of structural stress to vertical stress12Vertical motion rate u of broken block13And the range of influence u of fracture structure14Mining depth u21Coefficient of variation u of coal seam thickness22Fault dip angle u23Uniaxial compressive strength u of coal24Distance u between thick and hard rock stratum and coal seam25Sufficient collapse degree u of overburden stratum31Protective layer pressure relief degree u32The distance u between the coal bed and the upper coal pillar33The relation u of the working face to the adjacent goaf34Section coal pillar width u35Thickness u of coal left at bottom36(ii) a The hierarchical structure is shown in FIG. 2;
s4.2: comparing the importance degree of the index belonging to the previous layer with each other between the indexes of the same layer in the evaluation factor set, and establishing a judgment matrix Z ═ Z (Z)ij)n×n，Zij>0，
Wherein Z isijThe factor i is greater than the importance of the factor j relative to the index of the previous level, ZjiThe factor j is more important than the factor i relative to the index of the previous levelN is the order of the judgment matrix;
the importance judgment matrix of the criterion layer relative to the target layer is as follows:
the importance judgment matrixes of indexes of the index layer to the criterion layer are respectively Z1,Z2And Z3
S4.3: calculating a weight vector of the judgment matrix;
normalizing each column vector of the judgment matrix Z to obtain:
i,j＝1,2,…,n (18)
b is toijSumming by rows:
to wiNormalization process
Wherein the content of the first and second substances,
to determine the normalized arithmetic mean of each column vector in the matrix,
judging the weight vector of the matrix for the solved feature vector;
s4.4: calculating the maximum characteristic root of the judgment matrix according to the product ZW of the judgment matrix and the corresponding weight vector:
s4.5: calculating a consistency index CI, wherein the larger the consistency index CI is, the more serious the inconsistency degree of the matrix is judged; wherein the consistency index is as follows:
s4.6: selecting average random consistency index RI shown in Table 4, and calculating consistency ratio of judgment matrix
When the consistency ratio of the judgment matrix is less than 0.1, the judgment matrix meets consistency test, and takes the characteristic vector of the judgment matrix as a weight vector, otherwise, the judgment matrix does not meet the consistency test, and reconstructs the judgment matrix;
TABLE 4 average random consistency index RI
The weight vector W is calculated as [0.0430,0.1598,0.0232,0.0859,0.1891,0.0328,0.0805, 0.1114,0.0766,0.0160,0.0270,0.0270,0.0455,0.0719,0.0103]Maximum feature root of each decision matrixλmaxThe consistency ratios CR of the judgment matrix are shown in the table 5, and the calculated consistency ratios CR all meet the consistency check requirement;
TABLE 5 maximum feature root λ of decision matrixmaxAnd the consistency ratio CR of the judgment matrix
S5: calculating the membership degree of each impact danger grade according to the evaluation matrix and the weight set by adopting a fuzzy comprehensive evaluation method, and determining the rock burst danger grade according to a maximum membership degree principle;
after determining the weight vector W and the evaluation matrix R, carrying out synthetic operation to obtain a fuzzy comprehensive evaluation set B, namely
To synthesize an operator, bjCalled as fuzzy comprehensive evaluation index; after the fuzzy comprehensive evaluation index is solved, a maximum membership method can be adopted to determine a final evaluation result; the fuzzy relation matrix R and the weight vector W are substituted into a fuzzy comprehensive evaluation formula
Among them, the evaluation vector B can be obtained as (0.31, 0.60, 0.44, 0.30); and determining the grade of the rock burst of the ore as a weak impact danger according to the maximum membership principle, wherein the weak impact danger is consistent with the actual rock burst appearance condition.
Claims (9)
1. The coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment is characterized by comprising the following steps: comprises the following steps:
s1: on the basis of statistical analysis of the occurrence of rock burst disasters, determining a judgment factor set U of the rock burst dangerousness, wherein U is (U ═ U-11，u12，u13，u14，u21，u22，u23，u24，u25，u31，u32，u33，u34，u35，u36) Obtaining values of all judgment factors according to field test, geological data analysis and mine design related data;
s2: determining a judgment set of rock burst dangerousness according to the apparent intensity of rock burst;
s3: analyzing the relationship between each factor and rock burst development strength, establishing a fuzzy relationship between each evaluation factor and an evaluation grade in an evaluation set, and calculating by combining the magnitude of each factor to obtain a fuzzy relationship matrix;
s4: determining the importance degree of each factor of the rock burst to the previous level by using an analytic hierarchy process, and calculating to obtain the weight of each rock burst influence factor to obtain the weight vector of the rock burst influence factor;
s5: and calculating the membership degree of each impact danger grade according to the fuzzy relation matrix and the weight vector by adopting a fuzzy comprehensive evaluation method, and determining the danger grade of the rock burst according to a maximum membership degree principle.
2. The coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment as claimed in claim 1, wherein: the judging factor set of the rock burst risk in the S1 comprises a geological dynamic environment type factor U1Geological factor U2And mining factor U3。
3. The coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment as claimed in claim 2, wherein: the geomechanical environment type factor U1Including making a concave angle of reflection u11Ratio u of structural stress to vertical stress12Vertical motion rate u of broken block13And the range of influence u of fracture structure14。
4. The hierarchical analysis and fuzzy heddle based on claim 2The coal mine rock burst prediction method based on combined judgment is characterized by comprising the following steps: the geological factor U2Including the depth of mining u21Coefficient of variation u of coal seam thickness22Fault dip angle u23Uniaxial compressive strength u of coal24Distance u between thick and hard rock stratum and coal seam25。
5. The coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment as claimed in claim 2, wherein: the mining factor U3Including a degree of substantial collapse u of overburden31Protective layer pressure relief degree u32The distance u between the coal bed and the upper coal pillar33The relation u of the working face to the adjacent goaf34Section coal pillar width u35Thickness u of coal left at bottom36。
6. The coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment as claimed in claim 1, wherein: the specification of the S2 is in accordance with the rules for preventing coal mine from impacting and thinning.
7. The coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment as claimed in claim 1, wherein: the evaluation set V of the risk of rock burst in S2 includes: without risk of impact v1Risk of weak impact v2Moderate impact risk v3And the danger of strong impact v4。
8. The coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment as claimed in claim 1, wherein: the S3 establishes a fuzzy relationship between each of the evaluation factors and the evaluation level in the evaluation set, including:
s3.1: determining membership degree vectors of each evaluation factor in the evaluation factor set to each evaluation grade in the evaluation set according to the known relationship between each factor and rock burst showing strength;
s3.2: the membership degrees of all the factors form a fuzzy relation matrix R together.
9. The coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment as claimed in claim 1, wherein: the S4 specifically includes the following steps: according to the category of rock burst influencing factors, including the geological dynamic environment type factor U1Geological factor U2And mining factor U3Establishing a hierarchical structure model for rock burst prediction and providing a multilevel index system for rock burst prediction, wherein the method comprises the following steps;
s4.1: establishing a hierarchical structure model for rock burst prediction, wherein a target layer is a rock burst risk U, and a criterion layer is a geological dynamic environment factor U1Geological factor U2And mining factor U3The index layer is constructed with a concave reflection angle u11Ratio u of structural stress to vertical stress12Vertical motion rate u of broken block13And the range of influence u of fracture structure14Mining depth u21Coefficient of variation u of coal seam thickness22Fault dip angle u23Uniaxial compressive strength u of coal24Distance u between thick and hard rock stratum and coal seam25Sufficient collapse degree u of overburden stratum31Protective layer pressure relief degree u32The distance u between the coal bed and the upper coal pillar33The relation u of the working face to the adjacent goaf34Section coal pillar width u35Thickness u of coal left at bottom36；
S4.2: and comparing the influence degrees of the indexes of the same layer in the evaluation factor set on the indexes belonging to the upper layer pairwise to form a judgment matrix:
Z＝(zij)n×n
zij＞0
wherein z isijZ is the importance of the factor i to the factor j relative to the upper level indexjiThe factor j is more important than the factor i relative to the upper-level index, and n isJudging the order of the matrix;
s4.3: calculating a weight vector of the judgment matrix;
normalizing each column vector of the judgment matrix to obtain:
b is toijSumming by rows:
to wiNormalization treatment:
wherein the content of the first and second substances,
to determine the normalized arithmetic mean of each column vector in the matrix,
judging the weight vector of the matrix for the solved feature vector;
s4.4: calculating the maximum characteristic root of the judgment matrix according to the product ZW of the judgment matrix and the corresponding weight vector:
s4.5: calculating a consistency index CI:
s4.6: setting average random consistency index RI, and calculating consistency ratio of judgment matrix
I.e. the ratio of the average random consistency index to the average random consistency index, sets a threshold C1When the consistency ratio of the matrix is larger than C1If so, judging that the matrix does not meet the consistency check, and returning to S4.1 to reconstruct the judgment matrix; when the consistency ratio of the matrix is less than C1And then, judging that the matrix meets the consistency check, and taking the characteristic vector of the judgment matrix as a weight vector.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111523840.2A CN114169789A (en) | 2021-12-14 | 2021-12-14 | Coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111523840.2A CN114169789A (en) | 2021-12-14 | 2021-12-14 | Coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114169789A true CN114169789A (en) | 2022-03-11 |
Family
ID=80486236
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111523840.2A Withdrawn CN114169789A (en) | 2021-12-14 | 2021-12-14 | Coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114169789A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115640632A (en) * | 2022-10-14 | 2023-01-24 | 港珠澳大桥管理局 | Method, device, equipment and medium for evaluating service state of steel box girder bridge coating |
| CN115860582A (en) * | 2023-02-28 | 2023-03-28 | 山东科技大学 | Intelligent impact risk early warning method based on self-adaptive lifting algorithm |
| CN117108295A (en) * | 2023-10-20 | 2023-11-24 | 山东能源集团有限公司 | Advanced presplitting area prevention and control method for rock burst of high-position thick hard rock stratum |
-
2021
- 2021-12-14 CN CN202111523840.2A patent/CN114169789A/en not_active Withdrawn
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115640632A (en) * | 2022-10-14 | 2023-01-24 | 港珠澳大桥管理局 | Method, device, equipment and medium for evaluating service state of steel box girder bridge coating |
| CN115860582A (en) * | 2023-02-28 | 2023-03-28 | 山东科技大学 | Intelligent impact risk early warning method based on self-adaptive lifting algorithm |
| CN117108295A (en) * | 2023-10-20 | 2023-11-24 | 山东能源集团有限公司 | Advanced presplitting area prevention and control method for rock burst of high-position thick hard rock stratum |
| CN117108295B (en) * | 2023-10-20 | 2024-03-08 | 山东能源集团有限公司 | Advanced presplitting area prevention and control method for rock burst of high-position thick hard rock stratum |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhu et al. | A prediction method of coal burst based on analytic hierarchy process and fuzzy comprehensive evaluation | |
| CN114169789A (en) | Coal mine rock burst prediction method based on hierarchical analysis and fuzzy comprehensive judgment | |
| CN109447837B (en) | Danger assessment method for rock burst in mining area | |
| US20160070828A1 (en) | Vulnerability Assessment Method of Water Inrush from Aquifer Underlying Coal Seam | |
| CN104533443B (en) | A kind of tight roof Rock Burst disasters danger prediction analysis method | |
| CN109798106B (en) | Method for predicting risk of rock burst and prevention and treatment measures | |
| He et al. | Integrated rockburst early warning model based on fuzzy comprehensive evaluation method | |
| CN109798149B (en) | Method for predicting inflow intensity of mining field of thick coal seam with hard roof in grading manner | |
| CN113958366B (en) | Dynamic quantitative early warning method for impact risk based on vibration-stress double-field monitoring | |
| CN105741029A (en) | Cumulative stress concentration factor based impact risk assessment method | |
| CN111222254A (en) | Working face rock burst danger grade dividing method and system based on stress superposition method | |
| Zhou et al. | Rockburst prediction for hard rock and deep-lying long tunnels based on the entropy weight ideal point method and geostress field inversion: a case study of the Sangzhuling Tunnel | |
| Wang et al. | Risk Assessment of Water Inrush in Karst Tunnels Based on the Efficacy Coefficient Method. | |
| CN109918697B (en) | Method and device for predicting mine rock burst strength | |
| CN110298107A (en) | A kind of working face shock hazard evaluation method based on increment superposition | |
| Xue et al. | An analytical model for assessing soft rock tunnel collapse risk and its engineering application | |
| Zhang et al. | Risk assessment of floor water inrush in coal mines based on MFIM-TOPSIS variable weight model | |
| Liu et al. | Assessment of water inrush risk using the principal component logistic regression model in the Pandao coal mine, China | |
| CN113339073B (en) | Impact risk evaluation method based on roof rock stratum structure | |
| CN111046595A (en) | Typical and atypical rock burst mine type dividing method | |
| Liu et al. | Assessing water and sand inrushes hazard reductions due to backfill mining by combining gis and entropy methods | |
| Chi et al. | Prediction of top-coal caving and drawing characteristics by the analytic hierarchy process-fuzzy discrimination method in extra-thick coal seams | |
| CN109918696A (en) | A kind of classification method and device of bump strength grade | |
| CN115146971A (en) | Rock burst mine distinguishing and grading method based on equivalent depth | |
| CN114692410A (en) | Mine dynamic disaster estimation method based on geostress field inversion numerical range analysis |
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 | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20220311 |
|
| WW01 | Invention patent application withdrawn after publication |