Environmental Evaluation of Coal Mines Based on Generalized Linear Model and Nonlinear Fuzzy Analytic Hierarchy

To create a good working environment for workers in coal mines and ensure that they can work safely and comfortably, environmental assessments of coal mines represent an important method used to achieve this goal. In this paper, a comprehensive environmental evaluation model of coal mines was established based on generalized linear theory and fuzzy analytic hierarchy processes. First, the importance degree of each index factor was obtained by analyzing the statistical source data of air, water, soil, ecological compensation, and other indexes that affect the ecological environmental safety of coal mines using generalized linear theory. Based on the importance degree of each index factor, a pairwise comparison matrix was constructed by the logarithmic fuzzy preference programming method, and the influence weight value of each index factor was accurately obtained through MATLAB software. The results indicated that the model can effectively reduce the impact of expert subjectivity on the evaluation results, which is consistent with the reality of coal mines. In addition, in the process of environmental assessment of coal mines, the model is proven to be convenient, precise, and easy to operate.


Introduction
Resource safety is an important part of national strategic safety, and the development of mineral resources is always an important mission [1][2][3]. In particular, in the past 30-40 years, the development of mineral resources has provided strong support and protection for the rapid development of the economy and society. Due to technological reasons and the pursuit of high profits, the strength of mine ecological protection and restoration management has not been as strong as people would prefer, leading to the intensification of the contradiction between the coal mine development and ecological safety, e.g., the high production rate of three wastes of coal mines, the long potential harm of radiation, and the serious pollution to the air, water and soil [2]. In particular, the disordered development in some areas has led to geological disasters, such as surface collapse and debris flow, which have intensified land desertification and led to the degradation of ecological environmental quality [3]. Therefore, it is urgent to explore a new technology of mine geological environmental governance and ecological restoration to prevent the further degradation of ecological health. The current ecological environmental safety assessment of coal mines can not only measure the overall status of the ecological quality of coal mines but also reflect the governance of the ecological environmental pollution. Additionally, it can quickly identify the important factors that cause the deterioration of the mine environment, allowing for the quick adoption of the most effective measures to control the "culprits" that make the mine environment worse.
Currently, many scholars have constructed a series of safety evaluation models, e.g., analytic hierarchy process, fuzzy comprehensive evaluation, grey theory, and neural network method, and applied them to mine ecological environmental safety evaluations. For instance, Zhao et al. [4] established an evaluation index system of the mine geological environment using the analytic hierarchy process (AHP) and applied this method to the mining field of Hainan Island; the authors found that the evaluation value basically conformed to the situation of the mine site. Bascetin et al. [5] constructed a decision support system based on an analytical hierarchy process to evaluate the production of coal mines in Turkey. Samantra et al. [6] created an improved decision method based on fuzzy set theory and distinguished the ranges of evaluated risk ratings, ultimately suggesting an action requirement plan to guide the management of the mine environment. Yang et al. [7] selected four factors that affect the mine environment, including the water environment, atmospheric environment, geological disaster, and ecological environmental characteristics, as the evaluation factors to build a fuzzy comprehensive evaluation model of the mine environment for the quantitative study of the mine environment. Su et al. [8] established an engineering geological evaluation model based on grey mathematics theory and evaluated the engineering geological situation of typical mines by determining the evaluation weight index. Zhao et al. [9] used the principle of BP neural networks and the random interpolation method to build an evaluation model of the BP artificial neural network on mine ecological safety.
Each evaluation model has its own shortcomings. To fully consider the influence of all factors and present the comprehensive evaluation results scientifically and accurately, many scholars have combined various models to carry out comprehensive evaluations according to the characteristics of the different models, namely, the combined comprehensive evaluation model. For instance, based on fuzzy theory, Chen et al. [10] collected the key factors influencing the mine environment by remote sensing technology, determined the influence weight of the environmental evaluation factors by using AHP, and finally built the AHP-fuzzy mine environment evaluation model. Besides, Wei et al. [11] established the comprehensive evaluation of the underground heat source based on the AHP and fuzzy comprehensive evaluation (FCE) to construct the evaluation model, and some suggestions on thermal stress control are put forward. Based on the theories of uncertainty measurement evaluation (WME) and AHP, Dong et al. [12] established the comprehensive risk evaluation of underground goafs using multiple indexes. Dong et al. [13] also established and calibrate the comprehensive models using the random forest, support vector machine, Bayes' classifier, Fisher's classifier, logistic regression, and neural networks. Based on grey system theory, Wei et al. [14] used the analytic hierarchy process to determine the weight of each influencing factor and constructed a grey comprehensive evaluation model of the mine geological environment. To scientifically and accurately obtain the membership degree of each influencing factor, Liu et al. [15] built a fuzzy-grey combination evaluation model based on grey correlation analysis and fuzzy evaluation theory.
However, there are many factors that affect ecological environmental safety, and the relationship between the evaluation factors and the safety level is fuzzy, which makes the shortcomings of the current evaluation model more obvious. For instance, although the neural network model can overcome the subjectivity of expert experience and ensure a higher accuracy of evaluation results, it is easy to fall into the dilemma of a local minimum and slow convergence speed in the calculation process of weight, and this method cannot address the relationships between various index factors clearly. Although the AHP and grey theory evaluation methods can deal with the relationship of each index well, they rely on the experience of experts for evaluation when determining the weight of influencing factors; therefore, they cannot avoid subjectivity. The fuzzy analytic hierarchy process is the best method to address the fuzzy relationship of all factors; however, it also cannot avoid the evaluation of experts when building the comparison matrix.
To overcome the subjectivity of evaluation experts in the evaluation of the influence degree of each index factor and to further improve the uncertainty of data processing, in this paper, a comprehensive evaluation model of the coal mine environment is established based on generalized linear theory and fuzzy level theory. First, the importance of each index in relation to the target index is obtained by the initial statistics of the source data of each lowest level index based on generalized linear theory. Then, according to the statistical value obtained by generalized linear theory, the pairwise comparison matrix is constructed by the logarithmic fuzzy preference programming method (LFPP). Finally, the influence weight value of each index factor is accurately obtained through MATLAB software. The model is applied to a coal mine for an environmental assessment.

Methodology
2.1. Triangular Fuzzy Number Theory. There are many kinds of distribution functions used to imitate probability theory, such as the normal distribution function, trapezoid distribution function, triangular fuzzy number, k-th parabola distribution, Cauchy distribution, and S-type distribution. The triangular fuzzy numbers (TFN) is widely adopted to express the decision perception of alternatives' performances with respect to each criterion [16]. In order to easy the judgment  3), and it is a kind of variable whose values are not numbers but words and sentences [17]. Therefore, the triangular fuzzy theory is mainly used to build a model to evaluate the ecological environment of coal mines in this paper. In a given region Z, for any z ∈ Z, the fuzzy set A can be defined as where ε A is the membership degree of z to A and ε A ðxÞ ∈ ð0, 1Þ. Let l, m, and u be the minimum possible value, the intermediate value, and the maximum possible value of a certain Ecological environmental risk assessment of uranium mine   Then, the three numbers (l, m, u) constitute a triangular fuzzy number,Ã. LetÃ = ðl, m, uÞ and l ≤ m ≤ u, and the function diagram of the triangular fuzzy number is shown in Figure 1. The membership function [18] is where the parameter m denotes the maximum possible value describing a fuzzy event, and l and u represent the upper and lower bounds of the boundary, respectively.
where μ i is the mathematical expectation; Y i is the dependent variable vector; gðxÞ is the connection function; X ij is the independent variable matrix; β j is the model parameter vector to be estimated; and ξ i is the interference term. Combined with equation (3), the observed value of the dependent variable can be connected with the expected value of the explanatory variable by linear addition; therefore, the relationship between the mathematical expectation μ i and the linear factors can be expressed as where g(μ i ) is a strictly monotonic function and is differentiable, which is called a continuous function. The expected value of the observed value of the dependent variable is related to the explanatory variable by equation (4-1). It is assumed that there are N sets of observation data, i.e., ðyjX i Þ = ðyjx 1 , x 2 ,⋯,x k Þ and ði = 1, 2,⋯,nÞ. y i ∈ ð0, 1, ⋯, jÞ is a multicategory dependent variable with j categories, indicating that the conclusion of the safety evaluation model has j grades. X = ðx 1 , x 2 , ⋯, x k Þ is the explanatory variable, i.e., the influencing factor affecting the conclusion grade of the model. It is assumed that ω i is the cumulative influence weight of a certain level y i under the condition of the index factor variable where ω i is the influence weight of an evaluation level. To better describe the relationship between the cumulative influence weight of the level and the index factor variables, in this paper, the relationship is transformed by logistic transformation, that is where the relationship between the logit (ω) transformed and the independent variable is linear, and its value range is (-∞, + ∞), which is conducive to its estimation. The regression equation can be obtained by substituting the estimated value of parameter β into equation (6) ω 2.3. The Fuzzy Analytic Hierarchy Process. First, the weight value of each index influencing the target value is determined by generalized linear theory, and then the importance of each pair of factors is judged by the weight of each index factor. According to the linguistic variable in Table 1, the relative importance of each pair of factors can be obtained. Then, the comparison matrixŨ can be created.
To facilitate the research, the comparison matrixŨ is processed by the logarithmic transformation [19,20] lnũ ij ≈ ln l ij , ln m ij , ln u ij À Á : By substituting equation (9) into equation (2), the following can be obtained

Geofluids
In this paper, to solve the problem that a specific value λ results in an inconsistency between weights and fuzzy judgements in U, the nonnegative deviation parameters of γ ij are introduced [17,18].
where G is a specific constant, which is used to determine the weights in the fuzzy judgements. The normalized priority calculation of the fuzzy pairwise comparison matrixŨ can be expressed by 3. Application

Determine the Weights of the Responding Index Factors.
The decision problem in the hierarchical structure is decomposed, as shown in Figure 2. The hierarchy consists of three levels, including goals, criteria and alternatives. This goal    6 Geofluids layer is mainly used to assess the ecological environmental risk of coal mines. The criteria and alternatives are located at the second level and the third level, respectively. The ecological environment monitoring data of several typical coal mines in China are selected to evaluate the ecological environmental safety. The ecological compensation index is generally calculated by comprehensively considering the benefits of the mine, the income expectations of leaders and miners, and the laws and regulations. The data of influencing factors of air, soil, and water indexes are mainly collected from the concentration of some harmful substances specified in relevant national standards. Based on the analysis of the evaluation indexes of the coal mine according to references [21][22][23][24][25][26], the ecological safety compensation index value and the index factor index value of air, soil, and water indexes of the coal mine are obtained, as shown in Tables 2  and 3.
The cumulative influence weight of the corresponding grade of the first level indicator and the influence weight of each grade under the second level indicator factor are obtained through the analysis of the generalized linear model (the influence weight of each grade indicator is mainly obtained through the Statistical Analysis System (SAS), as shown in Table 4), and then the maximum value of the influence weight of each grade is taken as the influence weight of the corresponding category. The weight of each index factor is used to judge the relative importance of each index factor. The maximum likelihood ratio of each category is 42.684, 35.033, 43.556, and 41.022, respectively, all of which are within the allowable error range, indicating that the generalized linear model has a good fit.
The influence weight of the grade corresponding to the second level index and its corresponding first level index is analyzed through the generalized linear model, and the comparison matrix of the first level index and the second level index can be obtained according to Table 5, as shown in Tables 6-10.  Table 11, which shows that the soil quality has the most important index influence on the environmental risk of coal mines among the four influencing factors, followed by water environmental quality, atmospheric environment, and ecological compensation indicators. Therefore, it is necessary to monitor the dangerous elements in soil and water on a regular basis and take timely measures to prevent exceeding the standard to control the dangerous sources during the early risk stage.   (Table 7) 0.8200 Water environmental quality (Table 8) 0.6230 Atmospheric environment (Table 9) 0.4819 Ecological compensation indicators (Table 10) 0.6017 7 Geofluids In addition, Table 11 shows that the concentration of coal (C1-1), concentration of Pb (C2-1), concentration of cadmium (C1-3), and concentration of lead (C1-2) have a great influence on the environmental risk of coal mines. However, the perfection of an ecological early warning mechanism (C6-2), proportion of environmental protection input to mine income (C6-3), concentration of SO 2 (C3-3), and concentration of radon (C3-4) have little influence on the environmental risk of coal mines. Table 11 also reveals that the environmental investment management of coal mines in the later stage has little impact on the environmental assessment. This result is mainly because in the later stage, the environment of coal mines has caused many irreversible losses and has had many human and financial costs. From a sustainable point of view, the compensation and management of the mine environment in the later period are definitely not a good long-term strategy.

Verification of the model's Validity.
According to the literature, the validity of a model is usually verified by calculating the value of consistency and the nonnegative deviation variable of the fuzzy pairwise comparison matrices, i.e., the consistency λ and nonnegative deviation variable γ are the key parameters to determine the validity of a model. There are two cases in which λ = 0 and λ ≠ 0. When λ is not equal to zero, the greater the value of parameter λ is, the better the consistency of the model is. When λ is equal to zero, if γ = 0, the consistency of the model is still good. If γ ≠ 0, it indicates that the model is inconsistent. The greater the nonnegative deviation variable γ is, the worse the consistency of the model is. With the help of MATLAB software, the value of λ in every matrix is obtained, as shown in Table 12. From the consistency index aspect, the proposed model is confirmed to be acceptable.
In the second perspective, we compare the results of the present method we used in this paper with the results of the representative comprehensive evaluation method [27]. The results of the representative comprehensive evaluation method are shown in Table 13. We have obtained similar results on some factors. In particular, the impact of the first level index factors on the safety environment of uranium mines is basically the same. From all above, the proposed model was proven to be valid.

Advantages and Disadvantages of the Present Model.
Based on generalized linear theory, the influence value of index factors on the environmental evaluation of coal mines was obtained using SAS software to count the original data of the index factors. Then, the fuzzy pairwise comparison matrices were established according to the influence value of index factors on the environmental evaluation, which greatly reduced the influence of human subjectivity on the construction of the comparison matrix. In addition, the linguistic variables and the TFN were introduced in the process of solving the factor weight, which enabled managers to quickly determine the relative importance of each pair of factors through the linguistic variables.
However, in the present model, although the importance degree of each factor was obtained using a linear generalized model to count the original data of each factor, the relative importance degree of each factor was judged by experts. Therefore, it is urgent to establish an evaluation formula and basis for the relative importance degree of each factor. In addition, according to the literature, although the

Conclusions
This paper develops a comprehensive environmental evaluation model of coal mines based on generalized linear theory and fuzzy analytic hierarchy processes.
(1) The model is consistent with the results of literature [27] and let the evaluation process more efficient (2) In this paper, the linguistic variables are a kind of variable whose values are not numbers but words or sentences, and they are closer to the feelings of the decision makers. Therefore, the present model can effectively solve the problem of multiclassification in safety evaluations (3) In addition, in the environmental evaluation model, based on prior knowledge of the sample statistical source, the importance degree of each index factor on the environmental factor index is preliminarily obtained through generalized linear theory, and the influence weight value of each index factor was accurately obtained based on the logarithmic fuzzy preference programming method, which overcomes the defect of artificial empowerment and has more accurate evaluation results

Data Availability
The data used to support the findings of this study are included within the article.

Conflicts of Interest
The authors declare that they have no conflicts of interest.