Risk Identication of Coal Spontaneous Fire Based on COWA Modied G1 Combination Weighting Cloud Model

: T o realize the scientific judgment of spontaneous fire risk in the coal mine, the spontaneous fire influence factors were analyzed from the three aspects of coal spontaneous combustion tendency, air leakage, and oxygen supply, heat storage and heat dissipation. And the basis for the evaluation of t spontaneous fire grade was constructed. Combination ordered weighted averaging (COWA) calculation was introduced to optimizes G1 subjective weighting, and a COWA modified G1 combined weighting cloud model was proposed to identify the spontaneous fire risk in the coal mine. Finally, the rationality of the model was verified with actual cases. The research results show that the spontaneous fire risk level in the Lingquan coal mine is relatively safe, which is consistent with the actual situation. And the spontaneous combustion tendency of coal is the leading factor affecting spontaneous fire risk.


Introduction
Internal-caused fire (coal spontaneous combustion) occupies the main body of mine fire. Therefore, it is of great significance to identify the main influencing factors of internal fire and scientifically identify the risk of internal fire to ensure the intrinsic safety of the coal mine system [1][2]. Many scholars have in-depth research on the danger of fire in mines, and the application methods include the CW-TOPSIS model [3], entropy weight matter-element extension theory [4][5], gray theory [6], set pair analysis [7][8][9], and so on. Although the above methods have achieved certain results, there are still the following problems: first, the evaluation index weights are subject to large subjective factors. Secondly, the primary and secondary relationship between internal factors of fire is not clear. Thirdly, the extreme value of the evaluation index affects comprehensive weighting.
Based on the above considerations, the causes of spontaneous fire are analyzed comprehensively in the mine, and a combination ordered weighted averaging (COWA) modified G1 comprehensive weighting method is proposed based on combination numbers. Combined with cloud model theory, the mine fire risk identification model is constructed, which provides a new method for the scientific identification of mine fire risk.

Theoretical Analysis of Causes of Fire in Mine
The mine internal fire is the result of the comprehensive effect of internal and external factors of the system, in which the internal factor is the spontaneous combustion tendency of coal. And the external factor includes the air leakage and oxygen supply conditions and heat storage and heat dissipation conditions [10]. Based on the above factors, the author comprehensively analyzes the causes of mine internal fire.

Conditions of spontaneous combustion propensity of coal
The degree of coal metamorphism reflects the physical and chemical properties of coal. The higher the degree of coal metamorphism, the higher the volatile content of coal, the easier the coal spontaneous combustion. The higher the moisture content in coal, the higher the degree of coal looseness, and the higher the oxidation rate. When sulfur in coal is oxidized at low temperature, the expansion of coal is loose, the oxidation area of coal is increased, and its decomposition products enhance the oxygen absorption of coal. Therefore, the higher the sulfur content of coal, the easier coal spontaneous combustion. The higher the ash content in coal, the less easily broken coal, the smaller the oxidation area, the lower the possibility of spontaneous combustion. Due to gas adsorption in coal, the contact between coal and oxygen is isolated, and the oxidation time of coal is prolonged. Therefore, the spontaneous combustion tendency of coal is determined by five factors, including the degree of carbonization and metamorphism of coal, water content, sulfur content, ash content, and coal seam gas content.

Air leakage oxygen supply conditions
The roof is harder, the filling quality is worse, and the greater the air is leakage. When the thickness of residual coal is large, the contact between residual coal and oxygen increases, and the risk of spontaneous combustion of coal is higher. The faster the working face advancing, the shorter the oxidation zone retention time, and the smaller the risk of spontaneous combustion. Ventilation management defects affect the air leakage in goaf. When the air leakage intensity increases, the risk of coal spontaneous combustion increases. Therefore, the oxygen supply conditions of air leakage are determined by four factors: roof lithology, residual coal thickness, advancing speed, and ventilation management.

Heat storage and heat dissipation conditions
Surrounding rock temperature is an important influencing factor of heat storage in goaf. The higher the surrounding rock temperature is, the higher the risk of coal spontaneous combustion is. When the depth of the coal seam is deep, the higher the initial temperature of coal is, the shorter the spontaneous combustion period of coal is. The complex geological structure in coal seam affects the mining speed and prolongs the contact time between coal and air. Therefore, the heat storage and dissipation conditions are determined by the surrounding rock temperature, coal seam depth, and coal seam geological structure.

Classification of spontaneous fire hazard
Based on the classification basis of mine spontaneous fire [11][12], the risk level of spontaneous fire and influencing factors are divided into four grades: I, II, III, IV, which represent safety, safer, more dangerous and dangerous respectively, as shown in Table 1.

G1 subjective weight determination
The G1 method is a subjective weighting method that can reflect the importance of indicators [13]. The calculation steps are as follows:

1)Expert sorted the importance of evaluation indicators to determine the sequence relationship of indicators;
2)Determine the importance ratio of adjacent indicators Xk-1 and Xk; Where, rk is the importance ratio of adjacent indicators, dk-1 is the importance of indicator Xk-1, and dk is the importance of indicator Xk.
3)Based on the importance ratio rk of adjacent indexes, the weight of indexes is calculated by G1 method; , −1 are the subjective weights of the kth, k-1th indicators determined by the G1 method, and n is the total number of indicators.

COWA Operator
The COWA operator arranges the indicator data in descending order, and combines the position of the indicators for weighting, reducing the influence of subjective extreme values on indicator weights. It is an objective weighting method [14][15]. The calculation steps are as follows: 1)The index data is processed in descending order to obtain a reconstructed data set: 0 ≥ 2)The position weighted +1 calculation of the data bj: Where −1 is the number of combinations of j data obtained from n-1 data.
3)Calculation of the absolute weight value of the indicator: ̅̅̅ = ∑ • =1 (5) Where ̅̅̅ is the absolute weight of the index and bj is the jth data value.

4)Calculation of relative weight values of indicators:
Where is the relative weight value of the index.

Combination weighting under COWA correction condition
To take into account the subjectivity of decision-makers and the objectivity of data, and reduce the impact of subjective weighting extremum on weight, a COWA modified G1 combination weighting method is proposed. Based on game theory, the optimal combination of subjective weighting and objective weighting is realized by establishing combination weights and minimizing the difference between weights [16]. The calculation steps are as follows: 1)Assuming that the number of index weighting methods is m, the number of weight combinations is: where ui is the linear combination coefficient and is the weight of each assignment method.
2)Combinatorial coefficients ui solving. If the optimal point of the combinatorial assignment method is realized, the optimization model can be constructed as: Then the first-order derivative condition for its optimal condition is: 3)Combined coefficients ui normalization processing: 4)Based on the above analysis, the optimal combination weights are: 4 Model for identifying the risk of Spontaneous fire in mines

Cloud Model Theory
Cloud models enable the uncertain transformation of qualitative concepts and quantitative descriptions by establishing a mapping relationship between quantitative and qualitative concepts [17][18].

1)Cloud model definition and its numerical characteristics
In mine endogenous fire evaluation, assuming that U is the theoretical domain corresponding to the values of endogenous fire indicators, C is the qualitative concept in endogenous fire evaluation indicators, x denotes cloud drops, and u(x) is the affiliation degree of any cloud drop x to C in the theoretical domain U, then The cloud model is represented by expectation (Ex), entropy (En ), and superb entropy (He). Ex is the center of the cloud graph, En characterizes the reliability of Ex, and He characterizes the uncertainty of En. The distribution interval of cloud drops x is [Ex-3En，Ex+3En].

2)Cloud Generator
The cloud generator is divided into forward and inverse cloud generators. The forwarding cloud generator mainly realizes the conversion from qualitative concept to quantitative, and the calculation steps are as follows: ① With expectation En, variance He 2 , generate Gaussian random number ′ = ( , 2 ).
③ Calculation of the determination of the indicator: where is the degree of determination of the index.
④ Cloud droplet interval construction, based on the above steps to form a cloud droplet (xi,ui), and then repeat steps ①~③, until the formation of N cloud droplets.

Cloud numerical characteristics of fire hazard indicators in mines
Referring to the related research results [19], the cloud number characteristics of the fire hazard index in the mine can be calculated as follows, Where Fmax is the upper limit of the index value, Fmin is the lower limit of the index value, and k is the degree of index fuzziness, which is taken as 0.05 here.
Then, based on the graded values of Table 1 indicators, the cloud numerical characteristics of the mine endogenous fire evaluation indicators are calculated by equation (14), and for variables with unilateral boundaries, their numerical characteristics are obtained in the form of boundary parameters [20][21], as shown in Table 2.

Comprehensive discriminative model construction
Based on the cloud numerical characteristics of the evaluation indexes of mine internal fires, Matlab software is used to generate the cloud diagram of evaluation indexes, determine the determinacy of each index under different hazard levels, and then combine the G1 combination weights under COWA correction conditions to obtain the comprehensive rating of mine internal fires.
= ∑ * =1 (15) Where is the single indicator determinant and * is the optimal combination weight.

Determination of index weights
The II 3 coal seam of Lingquan Coal Mine in Inner Mongolia is selected as an application example, the coal seam buried depth 447~470m; coal seam gas content 7.8m 3 /t; the geological structure is simple, the average thickness of coal remains 1.2m, the average speed of advance 9m/d, the average thickness 14.69m, the lithology of the roof is medium sandstone, the measured water content 3.19%, ash content 12.14%, sulfur content 0.16%, the temperature of the surrounding rock is about 23℃, and the natural firing period is 40d, which is easy to the spontaneous combustion coal seam. The values of qualitative and quantitative indicators were determined by expert experience, and based on the principle of weight determination by the G1 method, industry experts were hired to analyze the serial relationship and importance among the indicators, and calculate the subjective weights according to Eqs. (1)-(3), while determining the objective weights based on COWA operator according to Eqs. (4)- (6), and finally combining Eqs. (7)- (11) to determine the comprehensive weights of the indicators, as shown in Table 3.

Determination and analysis of risk grade
Based on the cloud numerical characteristics of mine endogenous fire indicators in Table 2, the risk level determinacy of each indicator was calculated by equation (13), and then combined with the comprehensive weights of each indicator of endogenous fire in Table 3, the comprehensive determinacy of the hazard evaluation model was calculated by equation (15) to determine the hazard level of endogenous fire in Lingquan coal mine. The determinacy of the risk level of endogenous fire in Lingquan coal mine was calculated as P(Ⅰ)=0.046，P(Ⅱ)=0.260，P(Ⅲ)=0.074， P(Ⅳ)=0.095, respectively, and it is known that the risk boundary of endogenous fire in Lingquan coal mine is level II according to the principle of maximum determinacy, which is consistent with the actual situation of the mine. The calculation method and evaluation criterion of reference [2] were selected to analyze the risk level of endogenous fire in the Lingquan coal mine, and the result was a safer level, and the results of the two methods were consistent, which showed the reasonableness of the COWA modified G1 combined assignment cloud model in the evaluation of endogenous fire in the mine.
Meanwhile, the weights of each index in Table 3 were analyzed, and the weights were arranged in descending order, in which the secondary index sequence was spontaneous combustion propensity condition of coal > air leakage and oxygen supply condition > heat storage and heat dissipation condition, and the tertiary main index sequence was coal seam depth > surrounding rock temperature > water content > ash content > sulfur content, and the results were approximately the same as the findings of references [1][2], which verified that the COWA modified G1 combined assignment weighting method was more effective in determining the endogenous fire risk. The scientific feasibility of the COWA modified G1 combination weighting method in determining the weights of fire indicators.

Conclusion
1)Considering the three aspects of coal spontaneous combustion tendency condition, air leakage oxygen supply, and heat storage, the grade evaluation basis of fire risk index in mine is constructed.
2)The COWA modified G1 combination weighting cloud model is proposed to identify the fire risk in the mine, and the Lingquan coal mine is taken as the engineering background for verification. The risk level is relatively safe, which is consistent with the actual scene, and is consistent with the conclusion of the CW-TOPSIS method.
3)Based on the COWA modified G1 combination assignment method, the weights of each index are analyzed in descending order, among which, among the secondary indexes, the propensity of coal to spontaneous combustion has the greatest influence on mine internal fire, and among the tertiary indexes, the depth of coal seam, the temperature of surrounding rock, water content and ash content are the main factors affecting mine internal fire. Based on the COWA modified G1 combination assignment method, the weights of each index are analyzed in descending order, among which, among the secondary indexes, the propensity of coal to spontaneous combustion has the greatest influence on mine internal fire, and among the tertiary indexes, the depth of coal seam, the temperature of surrounding rock, water content and ash content are the main factors affecting mine Spontaneous fire.