CN111859288B - Goaf spontaneous combustion risk prediction method - Google Patents

Goaf spontaneous combustion risk prediction method Download PDF

Info

Publication number
CN111859288B
CN111859288B CN202010540009.7A CN202010540009A CN111859288B CN 111859288 B CN111859288 B CN 111859288B CN 202010540009 A CN202010540009 A CN 202010540009A CN 111859288 B CN111859288 B CN 111859288B
Authority
CN
China
Prior art keywords
index
weight
spontaneous combustion
goaf
evaluation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010540009.7A
Other languages
Chinese (zh)
Other versions
CN111859288A (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.)
Shenyang Research Institute Co Ltd of CCTEG
Original Assignee
Shenyang Research Institute Co Ltd of CCTEG
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 Shenyang Research Institute Co Ltd of CCTEG filed Critical Shenyang Research Institute Co Ltd of CCTEG
Priority to CN202010540009.7A priority Critical patent/CN111859288B/en
Publication of CN111859288A publication Critical patent/CN111859288A/en
Application granted granted Critical
Publication of CN111859288B publication Critical patent/CN111859288B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0635Risk analysis of enterprise or organisation activities
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Mining
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Abstract

The invention relates to the technical field of coal mine safety engineering, in particular to a goaf spontaneous combustion risk prediction method, which is characterized in that according to a fuzzy statistics principle, an AHP-entropy combined weight method is improved to simulate a human decision thinking process, and expert scoring evaluation quantification is carried out on uncertainty multi-factors causing spontaneous combustion risk of coal in a working face goaf, so that a numerical value capable of describing spontaneous combustion risk of the goaf is obtained, and further, the spontaneous combustion risk of the goaf is objectively judged, so that scientific and reasonable goaf spontaneous combustion preventive measures are adopted. According to the method, the spontaneous combustion risk of the goaf is researched by using a fuzzy mathematical prediction model, the importance of each index is judged by using an improved combined weighting method, the weight of each index is calculated, and the spontaneous combustion risk of the goaf is comprehensively evaluated by combining a gray theory method, so that the weight ranking of each factor of the spontaneous combustion risk level of the goaf can be rapidly and accurately performed, and the estimated problem of objectively evaluating the spontaneous combustion risk by using less information can be solved.

Description

Goaf spontaneous combustion risk prediction method
Technical Field
The invention relates to the technical field of coal mine safety engineering, in particular to a goaf spontaneous combustion risk prediction method.
Background
Spontaneous combustion in goaf is an internal factor fire influenced by multiple factors, and due to the fact that mining depth is increased in recent years, surrounding rock temperature is increased and other factors influence, the danger of spontaneous combustion in goaf is more and more serious, spontaneous combustion in goaf is usually caused by deep coal, only smoke is generated without naked flame, and the position of a fire source is difficult to determine. The spontaneous combustion fire of the goaf not only causes the waste of a large number of coal resources and the damage of material equipment, but also often accompanies the secondary damage of dust and gas explosion, the spontaneous combustion of the coal of the goaf at present causes the fire accounting for about 90% of the total fire, 46% -49% of key coal mines in China have spontaneous combustion dangers of the coal of the goaf, so that the prediction and evaluation of whether the goaf is influenced by the spontaneous combustion dangers of the coal is a urgent task, how to efficiently predict the spontaneous combustion dangers of the goaf is already the working key and the difficulty of scientific researchers, and the research on the spontaneous combustion dangers of the goaf has important significance.
Since spontaneous combustion fire in a goaf of a coal mine is a dynamic complex system with burstiness and uncertainty, the probability of occurrence of the spontaneous combustion fire is difficult to evaluate by linear change. When the system is evaluated, the primary and secondary relations of different factors need to be comprehensively considered, and the weight coefficient of the accident occurrence is determined according to the importance ranking of the factors. The existing research has important significance for prediction and prevention of spontaneous combustion danger of goaf, but the built model has the problems of difficult weight determination, incomplete index factor consideration, ambiguous primary and secondary relationship of influencing factors and the like in actual operation.
Disclosure of Invention
In order to solve the problems, the invention provides a goaf spontaneous combustion risk prediction method, which is characterized in that according to a fuzzy statistics principle, an AHP-entropy combined weight method is improved to simulate a human decision thinking process, and expert scoring evaluation quantification is carried out on uncertainty multi-factors causing spontaneous combustion risk of coal in a working face goaf, so that a numerical value capable of describing spontaneous combustion risk of the goaf is obtained, and further, the spontaneous combustion risk of the goaf is objectively judged, so that scientific and reasonable goaf spontaneous combustion preventive measures are adopted.
In order to achieve the above purpose, the invention adopts the following technical scheme:
a goaf spontaneous combustion risk prediction method comprises the following steps,
step 1: determining an evaluation index and determining an evaluation index weight: taking the oxygen uptake, CO temperature rise rate, coal sample ignition temperature difference, air leakage rate, coal layer thickness, coal layer inclination angle, air leakage duration, air leakage intensity, coal loss thickness, surrounding rock temperature and coal layer burial depth of the test point spontaneous combustion dangerous single index coal as evaluation indexes, and determining each index weight of the goaf spontaneous combustion danger by utilizing an improved AHP-entropy combination weight-gray theory method;
step 2: establishing a hierarchical analysis model for spontaneous combustion risk evaluation of the goaf: analyzing the relation and influence of each basic element in the goaf spontaneous combustion risk evaluation system, and dividing each index into an index layer, a criterion layer and a target layer three-level hierarchical structure; comparing the importance of each layer of index factors relative to a certain criterion of the previous layer according to a pairwise comparison scale method to construct a judgment matrix A= (a) ij ) n×n
Figure BDA0002537204070000021
Wherein a is ij Representing the ratio of index i to index j importance;
determining index weight by using an improved AHP weighting method, solving n times of square roots for each row of element products of the matrix, and solving a weight formula by using a geometric average method;
Figure BDA0002537204070000031
β={β 1 ,β 2 …β n } T
wherein: n is the index number;
Figure BDA0002537204070000032
representing the product of each row of elements of the matrix; beta i Representing the i-th index weight; beta is an index weight vector;
normalization processing
Figure BDA0002537204070000033
w A ={w A1 ,w A2 …w An } T To approximate the normalized eigenvector, where w Ai The weight corresponding to the fourth index factor is subjective weight, and
Figure BDA0002537204070000034
step 3: consistency test of the judgment matrix;
step 4: constructing decision matrix B 0 : according to the principle of qualitative and quantitative combination, evaluating m samples of spontaneous combustion risk levels, and an original data matrix B= (B) of n spontaneous combustion risk level evaluation indexes ij ) m×n Constructing decision matrix B by dimensionless and standardized processing 0 =(r ij ) m×n
Figure BDA0002537204070000035
Wherein b is ij Is an initial value; r is (r) ij Is a standardized index; b j max ,b j min Respectively representing the maximum value and the minimum value in index values of different samples under the index j;
step 5: determining an evaluation index entropy value: in order to avoid zero index feature specific gravity value, ensure that the evaluation index entropy has mathematical meaning, and r is ij ·ln r ij The influence on the index entropy is controlled in a reasonable range, and then the evaluation index entropy value is determined:
Figure BDA0002537204070000041
Figure BDA0002537204070000042
wherein: x is x ij Is a corrected evaluation index; psi is a correction coefficient; e, e j The entropy value is the j index entropy value; m is the number of samples;
step 6: calculating the information weight of the j-th index:
Figure BDA0002537204070000043
w E =(w E1 ,w E2 …w En )
wherein: w (w) Ej The weight corresponding to the j index in the entropy method is the objective weight, and
Figure BDA0002537204070000044
step 7: obtaining a comprehensive weight vector: coupling subjective weight and objective weight of each index factor, namely:
Figure BDA0002537204070000045
Q=(q 1 ,q 2 …q n )
wherein: q i Coupling weight for the ith index factor; q is a comprehensive weight vector;
step 8: establishing a comprehensive evaluation index system: screening all indexes by combining a Delphi method to construct a 3-level hierarchical structure comprehensive evaluation index system;
step 9: establishing a comment matrix M;
step 10: according to a combination weighting principle, determining the weight of each index, and constructing a judgment matrix R based on an improved AHP-entropy combination weight-gray theory method;
step 11: according to the comprehensive weight vector Q, the judgment matrix R and the judgment matrix M, the spontaneous combustion risk evaluation result of the goaf can be obtained, namely:
L=Q·R·M
wherein: l is an evaluation result; q is a comprehensive weight vector coupling subjective weight and objective weight; r is a comment matrix, and M is a comment matrix.
Preferably, in step 3, consistency test is performed on the judgment matrix to determine whether the weight distribution of each layer is reasonable;
Figure BDA0002537204070000051
Figure BDA0002537204070000052
Figure BDA0002537204070000053
wherein lambda is max Judging the maximum eigenvalue of the matrix; CI is a consistency index; IR is a random consistency index; CR is a consistency ratio.
Preferably, for the CR test: when CR is less than 0.1, the constructed judgment matrix meets the consistency requirement, otherwise, the judgment matrix needs to be corrected.
Preferably, in step 9, the spontaneous combustion risk level of the goaf is established, and the spontaneous combustion risk of the goaf is divided into a plurality of levels of comment matrixes M according to the safety production characteristics of the coal mine.
Preferably, the spontaneous combustion risk of the goaf is divided into a multi-level arithmetic progression comment matrix M.
Preferably, the goaf spontaneous combustion risk is classified into 5 grades of comment matrix M.
Preferably, in step 10, a whitening weight function of multiple evaluation levels of spontaneous combustion risk of the goaf is constructed according to a gray system theory;
evaluation index U ij Having a plurality of scores d ijt Indicating the t-th pair index U ij Giving a score, substituting the score into a whitening weight function, f k (d ijt ) Represents the t expert pair index U in the k-th gray ij And then sum to obtain U ij The kth gray evaluation coefficient X of (2) ijk Further obtain the total gray evaluation coefficient X of each index ij The method comprises the following steps:
Figure BDA0002537204070000061
Figure BDA0002537204070000062
X ijk and X is ij The ratio of the two is used to obtain an index U ij Belongs to the kth gray evaluation weight vector eta ijk Thereby obtaining the second level index U ij Gray evaluation weight vector η ij Comprehensive weight W obtained by combining improved AHP method and entropy value method ij A gray weight judging matrix R of each index is obtained, and a first-level index U is calculated according to the gray weight judging matrix R i Is the evaluation matrix of (1), namely:
Figure BDA0002537204070000063
η ij =(η ij1 ,η ij2 …η ijn )
Figure BDA0002537204070000064
to obtain a judgment matrix R.
The beneficial effects of using the invention are as follows:
according to the method, the spontaneous combustion risk of the goaf is researched by using a fuzzy mathematical prediction model, the importance of each index is judged by using an improved combined weighting method, the weight of each index is calculated, and the spontaneous combustion risk of the goaf is comprehensively evaluated by combining a gray theory method, so that the weight ranking of each factor of the spontaneous combustion risk level of the goaf can be rapidly and accurately performed, and the estimated problem of objectively evaluating the spontaneous combustion risk by using less information can be solved.
Drawings
FIG. 1 is a goaf spontaneous combustion risk prediction evaluation system hierarchical model;
fig. 2 is a modified AHP-entropy combination weight-gray theory method evaluation model flow.
Detailed Description
In order to make the objects, technical solutions and advantages of the present technical solution more apparent, the present technical solution is further described in detail below in conjunction with the specific embodiments. It should be understood that the description is only illustrative and is not intended to limit the scope of the present technical solution.
In combination with the summary of the invention, as shown in fig. 1, the goaf spontaneous combustion risk prediction method provided in this embodiment includes the following steps:
step 1: taking the oxygen uptake, CO temperature rise rate, coal sample ignition temperature difference, air leakage rate, coal layer thickness, coal layer inclination angle, air leakage duration, air leakage intensity, coal loss thickness, surrounding rock temperature and coal layer burial depth of the test point spontaneous combustion dangerous single index coal as evaluation indexes, and determining each index weight of the goaf spontaneous combustion danger by utilizing an improved AHP-entropy combination weight-gray theory method;
step 2: and establishing a hierarchical analysis model for spontaneous combustion risk evaluation of the goaf. And analyzing the relation and influence of each basic element in the goaf spontaneous combustion risk evaluation system, and dividing each index into an index layer, a criterion layer and a target layer three-level hierarchical structure by combining expert opinion. Comparing the importance of each layer of index factors relative to a certain criterion of the previous layer according to a pairwise comparison scale method to construct a judgment matrix A= (a) ij ) n×n Wherein a is ij The ratio of index i to index j importance is represented.
Figure BDA0002537204070000071
Determining index weight by using an improved AHP weighting method, solving n times of square roots for each row of element products of the matrix, and solving a weight formula by using a geometric average method;
Figure BDA0002537204070000081
β={β 1 ,β 2 …β n } T
wherein: n is the index number;
Figure BDA0002537204070000082
representing the product of each row of elements of the matrix; beta i Representing the i-th index weight; beta is an index weight vector.
Normalization processing
Figure BDA0002537204070000087
Thereby obtaining w A ={w A1 ,w A2 …w An } T I.e. an approximation of the normalized feature vector, where w Ai Is the weight corresponding to the ith index factor, i.e. subjective weight, and
Figure BDA0002537204070000083
step 3: and (5) consistency test of the judgment matrix. The judgment matrix which is usually constructed is difficult to meet the requirement of complete consistency, a certain standard needs to be established, and when the judgment matrix meets the standard, the judgment matrix is regarded as having complete consistency, so that consistency test needs to be carried out on the judgment matrix to determine whether the weight distribution of each layer is reasonable or not.
Figure BDA0002537204070000084
Figure BDA0002537204070000085
Figure BDA0002537204070000086
Wherein lambda is max Judging the maximum eigenvalue of the matrix; CI is a consistency index; IR is a random consistency index, which can be determined by looking up a table. CR is a consistency ratio, and when CR is less than 0.1, the constructed judgment matrix is considered to meet the consistency requirement, otherwise, the judgment matrix needs to be corrected.
Step 4: constructing decision matrix B 0 Raw data matrix b= (B) of n spontaneous combustion risk level evaluation indexes of m samples for spontaneous combustion risk level evaluation according to qualitative and quantitative combination principle ij ) m×n Constructing decision matrix B by dimensionless and standardized processing 0 =(r ij ) m×n
Figure BDA0002537204070000091
Wherein b is ij Is an initial value; r is (r) ij Is a standardized index; b j max ,b j min Respectively represent the maximum value and the minimum value in the index values of different samples under the index j.
Step 5: determining an evaluation index entropy value, ensuring that the evaluation index entropy has mathematical significance in order to avoid zero index characteristic specific gravity value, and adding r ij ·ln r ij The influence on the index entropy is controlled in a reasonable range, and partial factors in the decision matrix need to be corrected, so that the evaluation index entropy value is determined:
Figure BDA0002537204070000092
Figure BDA0002537204070000093
wherein: x is x ij Is a corrected evaluation index; psi is a correction coefficient; e, e j The entropy value is the j index entropy value; m is the number of samples.
Step 6: calculating the information weight of the j-th index:
Figure BDA0002537204070000094
w E =(w E1 ,w E2 …w En )
wherein: w (w) Ej The weight corresponding to the j index in the entropy method is the objective weight, and
Figure BDA0002537204070000101
step 7: in order to ensure that subjective knowledge and objective investigation data of decision makers on spontaneous combustion influence factors of the goaf can truly reflect rules, subjective weights and objective weights of all index factors are coupled to obtain a comprehensive weight vector, and the comprehensive weight vector is marked as Q.
Figure BDA0002537204070000102
Q=(q 1 ,q 2 …q n )
Wherein: q i Coupling weight for the ith index factor; q is the integrated weight vector.
Step 8: due to complexity and uncertainty of spontaneous combustion danger of a coal mine goaf, various indexes are screened by combining a Delphi method through field investigation and reference of related documents, and a 3-level hierarchical structure comprehensive evaluation index system is constructed, as shown in table 1.
Table 1 goaf spontaneous combustion risk evaluation index and expert scoring average value
Figure BDA0002537204070000103
Figure BDA0002537204070000111
Step 9: the spontaneous combustion danger grade of the goaf is established, the spontaneous combustion danger of the goaf is divided into 5 grades reasonably according to the safety production characteristics of the coal mine, and the goaf is respectively endowed with grade I (safe-basically impossible), grade II (safer-less likely), grade III (medium safe-likely), grade IV (less safe-more likely) and grade V (unsafe-highly likely) to form a comment matrix M, wherein the higher the index is, the smaller the spontaneous combustion danger of the goaf is caused by the index, and the smaller the spontaneous combustion probability of the goaf is caused by the index is.
Step 10: and constructing whitening weight functions of 5 evaluation grades of spontaneous combustion risks of the goaf according to a gray system theory.
Figure BDA0002537204070000112
Figure BDA0002537204070000113
Figure BDA0002537204070000114
Figure BDA0002537204070000121
Figure BDA0002537204070000122
Evaluation index U ij There are n experts giving scores, d ijt Representing the t expert pair index U ij Giving a score, substituting the score into a whitening weight function, f k (d ijt ) Represents the t expert pair index U in the k-th gray ij And then sum to obtain U ij The kth gray evaluation coefficient X of (2) ijk Further obtain the total gray evaluation coefficient X of each index ij
Figure BDA0002537204070000123
Figure BDA0002537204070000124
X ijk And X is ij The ratio of the two is used to obtain an index U ij Belongs to the kth gray evaluation weight vector eta ijk Thereby obtaining the second level index U ij Gray evaluation weight vector η ij Comprehensive weight W obtained by combining the improved AHP method and entropy value method ij A gray weight judging matrix R of each index is obtained, and a first-level index U is calculated according to the gray weight judging matrix R i Is a matrix of evaluation of (a).
Figure BDA0002537204070000125
η ij =(η ij1 ,η ij2 …η ijn )
Figure BDA0002537204070000126
Step 11: according to the combination weighting principle, determining the weight of each index, and constructing a judging matrix R based on an improved AHP-entropy combination weight-gray theory method to obtain the spontaneous combustion risk evaluation result of the goaf.
L=Q·R·M
Wherein: l is an evaluation result; q is a comprehensive weight vector coupling subjective weight and objective weight; r is a comment matrix, and M is a comment matrix.
The foregoing is merely exemplary of the present invention, and those skilled in the art can make many variations in the specific embodiments and application scope according to the spirit of the present invention, as long as the variations do not depart from the spirit of the invention.

Claims (7)

1. A goaf spontaneous combustion risk prediction method is characterized by comprising the following steps: comprises the steps of,
step 1: determining an evaluation index and determining an evaluation index weight: taking the oxygen uptake, CO temperature rise rate, coal sample ignition temperature difference, air leakage rate, coal layer thickness, coal layer inclination angle, air leakage duration, air leakage intensity, coal loss thickness, surrounding rock temperature and coal layer burial depth of the test point spontaneous combustion dangerous single index coal as evaluation indexes, and determining each index weight of the goaf spontaneous combustion danger by utilizing an improved AHP-entropy combination weight-gray theory method;
step 2: establishing a hierarchical analysis model for spontaneous combustion risk evaluation of the goaf: analyzing the relation and influence of each basic element in the goaf spontaneous combustion risk evaluation system, and dividing each index into an index layer, a criterion layer and a target layer three-level hierarchical structure; comparing the importance of each layer of index factors relative to a certain criterion of the previous layer according to a pairwise comparison scale method to construct a judgment matrix A= (a) ij ) n×n
Figure FDA0002537204060000011
Wherein a is ij Representing the ratio of index i to index j importance;
determining index weight by using an improved AHP weighting method, solving n times of square roots for each row of element products of the matrix, and solving a weight formula by using a geometric average method;
Figure FDA0002537204060000012
β={β 1 ,β 2 …β n } T
wherein: n is the index number;
Figure FDA0002537204060000013
representing the product of each row of elements of the matrix; beta i Representing the i-th index weight; beta is an index weight vector;
normalization processing
Figure FDA0002537204060000014
w A ={w A1 ,w A2 …w An } T To approximate the normalized eigenvector, where w Ai Is the weight corresponding to the ith index factor, i.e. subjective weight, and
Figure FDA0002537204060000021
step 3: consistency test of the judgment matrix;
step 4: constructing decision matrix B 0 : according to the principle of qualitative and quantitative combination, evaluating m samples of spontaneous combustion risk levels, and an original data matrix B= (B) of n spontaneous combustion risk level evaluation indexes ij ) m×n Constructing decision matrix B by dimensionless and standardized processing 0 =(r ij ) m×n
Figure FDA0002537204060000022
Wherein b is ij Is an initial value; r is (r) ij Is a standardized index; b jmax ,b jmin Respectively representing the maximum value and the minimum value in index values of different samples under the index j;
step 5: determining an evaluation index entropy value: in order to avoid zero index feature specific gravity value, ensure that the evaluation index entropy has mathematical meaning, and r is ij ·lnr ij The influence on the index entropy is controlled in a reasonable range, and then the evaluation index entropy value is determined:
Figure FDA0002537204060000023
Figure FDA0002537204060000024
wherein: x is x ij Is a corrected evaluation index; psi is a correction coefficient; e, e j The entropy value is the j index entropy value; m is the number of samples;
step 6: calculating the information weight of the j-th index:
Figure FDA0002537204060000031
w E =(w E1 ,w E2 …w En )
wherein: w (w) Ej The weight corresponding to the j index in the entropy method is the objective weight, and
Figure FDA0002537204060000032
step 7: obtaining a comprehensive weight vector: coupling subjective weight and objective weight of each index factor, namely:
Figure FDA0002537204060000033
Q=(q 1 ,q 2 …q n )
wherein: q i Coupling weight for the ith index factor; q is a comprehensive weight vector;
step 8: establishing a comprehensive evaluation index system: screening all indexes by combining a Delphi method to construct a 3-level hierarchical structure comprehensive evaluation index system;
step 9: establishing a comment matrix M;
step 10: according to a combination weighting principle, determining the weight of each index, and constructing a judgment matrix R based on an improved AHP-entropy combination weight-gray theory method;
step 11: according to the comprehensive weight vector Q, the judgment matrix R and the judgment matrix M, the spontaneous combustion risk evaluation result of the goaf can be obtained, namely:
L=Q·R·M
wherein: l is an evaluation result; q is a comprehensive weight vector coupling subjective weight and objective weight; r is a comment matrix, and M is a comment matrix.
2. The goaf spontaneous combustion risk prediction method according to claim 1, wherein:
in step 3, consistency test is carried out on the judgment matrix to determine whether the weight distribution of each layer is reasonable or not;
Figure FDA0002537204060000041
Figure FDA0002537204060000042
Figure FDA0002537204060000043
wherein lambda is max Judging the maximum eigenvalue of the matrix; CI is a consistency index; IR is a random consistency index; CR is a consistency ratio.
3. The goaf spontaneous combustion risk prediction method according to claim 2, wherein:
for CR test: when CR is less than 0.1, the constructed judgment matrix meets the consistency requirement, otherwise, the judgment matrix needs to be corrected.
4. The goaf spontaneous combustion risk prediction method according to claim 1, wherein:
in step 9, the spontaneous combustion risk level of the goaf is established, and the spontaneous combustion risk of the goaf is divided into a multi-level comment matrix M according to the safety production characteristics of the coal mine.
5. The goaf spontaneous combustion risk prediction method according to claim 4, wherein: the spontaneous combustion danger of the goaf is divided into a multi-level arithmetic array comment matrix M.
6. The goaf spontaneous combustion risk prediction method according to claim 4, wherein: the spontaneous combustion risk of the goaf is classified into 5 grades of comment matrixes M.
7. The goaf spontaneous combustion risk prediction method according to claim 4, wherein:
in step 10, according to gray system theory, constructing whitening weight functions of multiple evaluation grades of spontaneous combustion risks of the goaf;
evaluation index U ij Having a plurality of scores d ijt Indicating the t-th pair index U ij Giving a score, substituting the score into a whitening weight function, f k (d ijt ) Represents the t expert pair index U in the k-th gray ij And then sum to obtain U ij The kth gray evaluation coefficient X of (2) ijk Further obtain the total gray evaluation coefficient X of each index ij The method comprises the following steps:
Figure FDA0002537204060000051
Figure FDA0002537204060000052
X ijk and X is ij The ratio of the two is used to obtain an index U ij Belongs to the kth gray evaluation weight vector eta ijk Thereby obtaining the second level index U ij Gray evaluation weight vector η ij Comprehensive weight w obtained by combining improved AHP method and entropy value method ij A gray weight judging matrix R of each index is obtained, and a first-level index U is calculated according to the gray weight judging matrix R i Is the evaluation matrix of (1), namely:
Figure FDA0002537204060000053
η ij =(η ij1 ,η ij2 …η ijn )
Figure FDA0002537204060000054
to obtain a judgment matrix R.
CN202010540009.7A 2020-06-12 2020-06-12 Goaf spontaneous combustion risk prediction method Active CN111859288B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010540009.7A CN111859288B (en) 2020-06-12 2020-06-12 Goaf spontaneous combustion risk prediction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010540009.7A CN111859288B (en) 2020-06-12 2020-06-12 Goaf spontaneous combustion risk prediction method

Publications (2)

Publication Number Publication Date
CN111859288A CN111859288A (en) 2020-10-30
CN111859288B true CN111859288B (en) 2023-06-23

Family

ID=72987044

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010540009.7A Active CN111859288B (en) 2020-06-12 2020-06-12 Goaf spontaneous combustion risk prediction method

Country Status (1)

Country Link
CN (1) CN111859288B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112308427B (en) * 2020-11-02 2024-03-19 江苏省电力试验研究院有限公司 New energy consumption constraint factor evaluation method and system based on combined weighting-ash correlation
CN112785292B (en) * 2021-01-29 2024-02-27 湖南大学 Prepaid card management method and system based on blockchain, equipment and storage medium
CN113408878A (en) * 2021-06-07 2021-09-17 山东电工电气集团有限公司 Cable tunnel running state evaluation method based on multi-source information and area control system
CN113344403B (en) * 2021-06-18 2023-06-27 安徽理工大学 Stability evaluation method for goaf construction site
CN113723817A (en) * 2021-08-31 2021-11-30 北京石油化工学院 Enterprise dust explosion risk assessment method, device and equipment
CN116151674A (en) * 2023-02-22 2023-05-23 辽宁艾特斯智能交通技术有限公司 Decision-making method, device, equipment and medium for maintenance of road electromechanical equipment
CN117291426B (en) * 2023-10-11 2024-03-12 中铁十九局集团矿业投资有限公司 Goaf safety risk level prediction method and goaf safety risk level prediction system based on Fisher discriminant method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110578555A (en) * 2019-09-12 2019-12-17 辽宁工程技术大学 method for predicting spontaneous combustion danger of gas extraction drill hole
CN110599033A (en) * 2019-09-12 2019-12-20 辽宁工程技术大学 Dynamic prediction method for goaf spontaneous combustion danger by introducing update factor

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20190339688A1 (en) * 2016-05-09 2019-11-07 Strong Force Iot Portfolio 2016, Llc Methods and systems for data collection, learning, and streaming of machine signals for analytics and maintenance using the industrial internet of things
US20180284758A1 (en) * 2016-05-09 2018-10-04 StrongForce IoT Portfolio 2016, LLC Methods and systems for industrial internet of things data collection for equipment analysis in an upstream oil and gas environment

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110578555A (en) * 2019-09-12 2019-12-17 辽宁工程技术大学 method for predicting spontaneous combustion danger of gas extraction drill hole
CN110599033A (en) * 2019-09-12 2019-12-20 辽宁工程技术大学 Dynamic prediction method for goaf spontaneous combustion danger by introducing update factor

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
"AHP+熵权法"的CW-TOPSIS煤矿内因火灾评价模型;秦忠诚;陈光波;李谭;孙伟;付彪;;西安科技大学学报(02);全文 *
基于AHP优化的采空区自燃危险性预测研究;郇华清;祁云;汪伟;齐庆杰;;能源技术与管理(02);全文 *
基于熵权物元可拓的采空区自燃危险性评价;左云飞;;煤炭技术(07);全文 *
基于物元可拓灰聚类模型的采空区火灾风险预测;吕辰;薛生;刘辉;徐晓萌;;煤炭技术(04);全文 *
采空区遗煤自燃危险性评价的熵权可拓方法;王金宝;侯金玲;张浪;舒龙勇;;矿业安全与环保(02);全文 *

Also Published As

Publication number Publication date
CN111859288A (en) 2020-10-30

Similar Documents

Publication Publication Date Title
CN111859288B (en) Goaf spontaneous combustion risk prediction method
CN110163500B (en) Fuzzy fault tree-based oil storage tank area fire explosion risk assessment method
Lin et al. Application of cloud model in rock burst prediction and performance comparison with three machine learning algorithms
CN107886235A (en) A kind of Fire risk assessment method for coupling certainty and uncertainty analysis
CN107967571A (en) Index system generation method based on entropy assessment
CN110889440A (en) Rockburst grade prediction method and system based on principal component analysis and BP neural network
Li et al. Real-time warning and risk assessment of tailings dam disaster status based on dynamic hierarchy-grey relation analysis
Zhou et al. A new approach for selecting best development face ventilation mode based on G1-coefficient of variation method
CN114065223A (en) Multi-dimensional software security risk assessment method based on CVSS
CN114049026A (en) Industrial internet identification analysis risk assessment method based on analytic hierarchy process
CN112508416A (en) Oil and gas storage and transportation station safety grade evaluation method based on cloud fuzzy analytic hierarchy process
CN115471097A (en) Data-driven underground local area safety state evaluation method
CN108776866B (en) Adaptability evaluation method for crude oil storage and transportation system of oil field and mine field
CN110020758A (en) Natural gas pipeline inherent peril appraisal procedure
CN116415836A (en) Security evaluation method for intelligent power grid information system
CN113722230B (en) Integrated evaluation method and device for vulnerability mining capability of fuzzy test tool
CN115358454A (en) Coal and gas outburst prediction method based on extension-fuzzy hierarchical analysis theory
Lin et al. Decision Tree analysis of the relationship between defects and construction inspection grades
Su et al. Research on evaluation of coal mine safety risk based on AHP-MCS coupling TOPSIS
CN111461558B (en) Key data safety destruction evaluation method
CN117196313A (en) Tunnel construction collapse accident coupling risk source identification method
CN108446459B (en) Fuzzy semantic reasoning-based coking process heat consumption influence factor optimization method
CN117787760A (en) Coal blending and burning fuzzy evaluation method and system based on game theory combined weighting
CN117094571A (en) Regional development status monitoring index determination method and device and electronic equipment
CN117196286A (en) Building construction safety risk pressure evaluation method based on toughness city theory

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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 113122 No.11 Binhe Road, Shenfu demonstration area, Fushun City, Liaoning Province

Applicant after: Shenyang Research Institute Co.,Ltd.

Address before: 113122 No. 11 Binhe Road, Fushun Economic Development Zone, Fushun City, Liaoning Province

Applicant before: CCTEG SHENYANG Research Institute

GR01 Patent grant
GR01 Patent grant