The Comprehensive Identiﬁcation of Roof Risk in a Fully Mechanized Working Face Using the Cloud Model

: Roof accidents seriously affect the safe and efﬁcient mining of the working faces. Therefore, it is necessary to assess and identify the possible and inﬂuencing factors on the occurrence of roof risk in a fully mechanized mining workface. In this study, based on the analytic hierarchy process and fuzzy comprehensive evaluation, a comprehensive standard cloud model was established through constructing a quantitative grade interval and calculating the weight of each index to achieve the aim of a roof risk assessment and identiﬁcation. The accuracy of risk assessment was ensured by using the comprehensive analyses of various aspects, such as cloud digital features, risk assessment cloud image and standard cloud image. This showed that the main inﬂuencing factors on the occurrence of roof accidents were roof separation distance, weighting intensity and rib spalling followed by the coal body stress concentration, initial support force and geological conditions. Taking 42,115 fully mechanized working faces in the Yushen coal mining area as an engineering background, this model was adopted to assess and identify the risk of roof accidents through generating comprehensive assessment cloud images and introducing the Dice coefﬁcient to calculate the similarity degree. The results showed that the overall risk of roof accidents in 42,115 working faces was regarded as grade II (general risk) through the overall index of comprehensive risk evaluation and a similarity degree of 0.8606. The impact of roof condition was mainly inﬂuenced by the risk of roof accidents, while the support status, personal working status and coal body condition had a limited effect on the risk of roof accidents. The comprehensive standard cloud model proposed in this study had strong visibility and discovered the key parts of risk indexes easily to solve the problems of ambiguity and quantitative identiﬁcation in traditional roof risk evaluation methods. Therefore, this model was worth promoting, because it laid the foundation for the intelligent identiﬁcation and early warning system of roof accident risk in a fully mechanized mining workface.


Introduction
The main disasters from underground coal mining activities in China involve coal dust, gas, water, fire, and roof accidents. With the increase in mining intensity, various disasters in coal mines occur frequently. From January 2020 to July 2020, a total of 63 coal mine accidents occurred in China causing 88 deaths. Among them, the numbers of coal mine roof accidents and casualties were relatively high. The accidents were mainly caused by factors such as the instability of the surrounding rock structure of the working faces, the failure of supports, rib spalling, and the development of roof cracks [1][2][3][4][5][6][7], which resulted in roof collapse or end-face roof collapse. More serious is the fact that such roof accidents occurred in fully mechanized working faces, directly affecting the safety of mining activities, causing damage to machinery and equipment and even threatening the safety of several people [8][9][10][11][12][13]. Table 1 describes the specific conditions of five roof accidents which occurred in China in recent years. It can be shown that the occurrence of roof accidents mainly depended on support condition, the mechanical properties of the roof, the quality of the personnel, and the stability of the coal body. Obviously, the occurrences of roof accidents not only caused huge economic losses and more casualties, but they also had serious effects on the coal mining community. Therefore, it was necessary to strengthen the roof risk assessment and identify the fully mechanized mining workface to ensure the safety and productivity of coal mines. Due to the influence of the geological structure, stress superposition and other factors, the hydraulic support was compressed during fault crossing.
Thus, the broken rock above the supports flowed into the working faces from the front of the top beam under the action of water. In the area of the supports 57#-58# close to the coal face, operators did not obey the relevant operator regulations, and, over time, the move supports resulted in the roof control distance exceeding the requirements of the operating regulations and further roof collapse. 23 June 2020 1221 air-return roadway, Ping'an coal mine, Hebei Province 1 Operators did not implement a temporary support, and the lumber was removed by hand causing coal and rock to cave. 19 August 2020 91,803 working faces, Xingyun coal mine, Yunnan Province 1 The roof collapsed due to the influence of geological structures (e.g., faults, a broken roof and stress concentration), and operators failed to reinforce the supports in time.
In recent years, various methods have been adopted to analyze the factors causing roof accidents and to accurately assess roof risk. Normally, the neural network and neuro-fuzzy models have been widely adopted in the risk assessment of engineering [14][15][16][17][18][19]. Scholars have proposed relevant models to assess the risk probability and to use decision-making principles to manage risks based on the analysis of most coal mine roof accidents [20]. In order to avoid the structural problem of the method, the triangular fuzzy theory and the traditional fault tree could be combined to quantitatively analyze the reasons for roof collapse accidents [21]. Furthermore, the dynamic fuzzy evaluation method was used to construct the risk matrix to determine the mine roof risk level, which further effectively predicted the occurrence of risk possibility [22][23][24]. In addition, roof separation was an obvious precursor for the occurrence of roof accidents in a fully mechanized working face. Therefore, the gray theory analysis method was used to predict the roof separation distance. Combined with the gray algebraic curve model, the early warning system could be established to evaluate the occurrence of roof accidents [25]. The classification of roof grade was closely related to the geological conditions (e.g., faults) and roof strength (e.g., roof caving), which further affected the roof risk assessment [26]. In order to obtain the roof risk assessment of the soft rock roadway, the analytic hierarchy process and the entropy weight method were used to determine the comprehensive weight of each evaluation index. Combining MATLAB software, a risk assessment model was constructed to comprehensively evaluate its stability. It showed that the lithology was the main factor affecting roof stability [27]. Moreover, according to the US National Safety Commission, more than 80% of coal mine roof accidents were caused by human factors (e.g., personnel attention and irregular operation) [28,29]. Scholars also proposed a model based on the rock engineering system (RES) to assess the risk of roof collapse. Additionally, the tendency of roof collapse could be divided into four categories of safe, moderate, dangerous and critical regions, respectively [30]. All these results could effectively identify the damaged and safety areas of the roof to avoid the occurrence of roof accidents.
Throughout the literature review, although the methods mentioned above could qualitatively evaluate and predict roof risk, they were difficult to use to identify the main factors causing roof accidents. Additionally, the specific quantitative characteristics of each factor and the possibility of risk accidents could not be illustrated. In recent years, the cloud evaluation model with the characteristics of qualitative, quantitative and strong visibility has been successfully used in the assessment of the ecological environment of coal mines [31], gas drainage [32], mine ventilation [33], slope stability analysis [34,35] and other aspects [36][37][38][39][40]. Therefore, in this study, the analytic hierarchy process (AHP) was adopted to identify the causing factors of roof accidents and calculate the index weight of each factor. Additionally, a cloud model was then introduced to assess the possibility of the occurrence of roof accidents by using cloud images and the quantitative grade assessment method. Obviously, it provided theoretical guidance for identification, prevention and controlling of roof accidents to ensure the safe and efficient production of coal mines.

Risk Evaluation Index System of Roof Accidents
Through summarizing and analyzing the reasons of loads of roof accident cases, the causing factors of risk accidents could mainly be divided into four levels of support conditions (e.g., support distances, the touch degree of the supports and roof, initial support force, the efficiency of a single support), roof conditions (e.g., hanging roof areas before the first weighting, hanging roof areas after the first weighting, hanging roof area during the movement of working faces equipment, roof separation distances, geological conditions, weighting intensity) [41,42], personnel working status (e.g., illegally moving the supports, moving the support in time, dealing with the hanging roof in time, the integrity of a single support) [43,44] and coal body status (e.g., rib spalling, soft coal seam and the stress concentration of coal body), including 17 indexes total. Table 2 describes the characteristics of each factor, and Figure 1 illustrates the comprehensive evaluation index system of the risk of roof accidents in a fully mechanized working face.
According to the risk evaluation index system (4 levels with 17 indexes total) of roof accidents, the risk grade could be divided into four categories of Grade I (low risk), Grade II (general risk), Grade III (larger risk) and Grade IV (major risk), which indicated the potential occurrence of roof accidents under the action of various causing factors. The standard of each risk grade corresponding to the value range of each index is illustrated in Table 3. Table 2. Description of each factor and the characteristics of the risk of roof accidents.

Support distances
If the hydraulic supports were overhauled, the support distances were too large, and the gangue slipped from the supports The touch degree of the supports and roof If roof and hydraulic supports were not tightly connected, the roof easily collapsed Initial support force If the initial support force was not up to the standard, the roof easily collapsed The efficiency of single support If the pillar of hydraulic supports was not intact and could not be effectively utilized, it easily caused roof collapse and hurt people Hanging roof areas before first weighting If the hanging areas of the main roof were too large before the first weighting, the roof collapsed suddenly, and it caused the occurrence of risk accidents If the hanging roof areas between the side of the working faces and the side of the goaf exceeded the regulations, and the relevant measures were not taken in time, the roof easily collapsed The integrity of single support If the conditions of the pillars and the upper roof were not checked before the pillars were moved, the pillar caps could fall off, and the high-pressure liquids of hydraulic supports were ejected to cause the occurrence of roof collapse and people being hurt According to the risk evaluation index system (4 levels with 17 indexes total) of roof accidents, the risk grade could be divided into four categories of Grade I (low risk), Grade II (general risk), Grade III (larger risk) and Grade IV (major risk), which indicated the potential occurrence of roof accidents under the action of various causing factors. The standard of each risk grade corresponding to the value range of each index is illustrated in Table 3.

Determining Each Risk Index Weight of Roof Accidents
The analytic hierarchy process (AHP) was a qualitative and quantitative analysis and evaluation method, which decomposed the elements of the decision-making problem into three levels of objectives, criteria and programs [45][46][47][48]. The evaluation processes of using AHP mainly included the establishment of the hierarchical structure model, the construction of the judgment matrix groups and the calculation of each index weight. Moreover, consistency testing was conducted to verify the credibility of calculated weight results. Subsequently, the relevant elements were compared in pairs to assess the risk grade. In this study, according to comparison results, the importance grade from 1 to 9 was adopted to build two judgment matrices from program level to criterion level and from criterion level to objective level, respectively. Additionally, the weight vector was calculated through the maximum eigenvalue (λ m ) of the judgment matrices. The corresponding eigenvector U j was expressed as follows: According to Equation (1), index weight (ω j ) and weight vector were calculated as follows: In order to ensure the credibility of the weight results, a consistency test needs to be conducted based on the maximum eigenvalue (λ m ) of judgment matrices. Additionally, the consistency index (CI) and consistency ratio (CR) were expressed as follows: where RI was random index. And it was 0.58, 0.9, 1.12 and 1.24 when the n was 3, 4, 5 and 6, respectively [45].
If CR was less than 0.1, the test was regarded as pass [49,50]. Otherwise, the test failed. Additionally, the judgment matrices then needed to be further adjusted until the consistency test passed.

Establishment of Judgment Matrices Groups and Calculation of Index Weight
Based on the risk evaluation analysis of roof accidents and the hierarchy index system, the Delphi method was adopted to construct the judgment matrix of each index at the criteria level and program level. Moreover, through an expert scoring system and the loads of the relevant references, the values from 1 to 9 were assigned to represent the importance grade of each risk index of roof accidents in the working faces as shown in Table 4. Table 4. Importance grade and its assignment value.

Importance Grade Assignment
Equally important 1 Slightly important 3 Obviously important 5 Strongly important 7 Extremely important 9 Note: 2, 4, 6 and 8 represent the assignments whose importance was between the corresponding values of 1, 3, 5, 7 and 9, respectively.
Firstly, the judgment matrix from criteria level to objectives level was shown as follows: Obviously, the eigenvalue and eigenvector of this judgment matrix were calculated as 4.0458 and (0.1564, 0.46213, 0.087717, 0.29376), respectively. According to Equation (1), the vector of index weight (ω A (B i )) was expressed as follows:  (4) and (5), the consistency (CI) index and consistency ratio (CR) were 0.015273 and 0.01697, respectively. Due to a CR of less than 0.1, the consistency test passed.
Secondly, the judgment matrices from the program level to each criterion level (support conditions (B 1 ), roof conditions (B 2 ), personnel working status (B 3 ) and coal body status (B 4 )) were shown as follows: Support conditions (B 1 ). .
Personnel working status (B 3 ). At the same time, based on Equations (1), (4) and (5), the vector of index weight, consistency index (CI) and consistency ratio (CR) of each criterion level and the corresponding results of the consistency test are also illustrated in Table 5. The eigenvalue and eigenvector of each judgment matrix are summarized in Table 6.  Table 6. Total ranking results of each risk index weight.

Calculation of Total Weight of Each Risk Index
Assuming that ω A (C j ) was the weight vector of program level (C) in objective level (A), the index weight was regarded as 0 if the index had a limited effect in program level (C). Based on Equation (6), the rank results of each risk index weight are shown in Table 6. Table 6 illustrates that the weight of three indexes (C 8 , C 15 and C 10 ) was larger than the others, followed by C 17 , C 3 and C 9 . This demonstrates that the most important factors causing roof accidents were roof separation distances, weighting intensity and rib spalling followed by the stress concentration of coal body, initial support force and geological conditions.
In view of the traditional coal mine roof risk assessment methods (e.g., fuzzy mathematics, AHP and other evaluation methods), although the factors causing accidents and their risk weights were obtained, it was difficult to intuitively express the specific importance level of risk factors and quantify the risk evaluation grade due to the ambiguity characteristics of these traditional methods. Therefore, the cloud theory was introduced to construct a comprehensive cloud model of roof risk assessment in the working faces and then combined with qualitative and quantitative measures to effectively identify the specific roof risk grade.

Cloud Model
The cloud model can transform the uncertainty relationship between qualitative and quantitative values into an intuitive and universal model, which mainly reflects the characteristics of fuzziness and randomness [51,52]. Normally, the digital features of expected value (E x ), entropy value (E n ) and super entropy value (H e ) are adopted to describe the characteristics of the cloud model as shown in Figure 2. These parameters could be obtained based on Equation (7) if the range of qualitative indexes was given: where k is the uncertainly cohesion degree of assessment indexes, which was assigned according to the ambiguity and randomness of each risk index. C max and C min are the maximum and minimum values of the level interval in the cloud model, respectively.

Cloud Model
The cloud model can transform the uncertainty relationship between qualitative and quantitative values into an intuitive and universal model, which mainly reflects the characteristics of fuzziness and randomness [51,52]. Normally, the digital features of expected value (Ex), entropy value (En) and super entropy value (He) are adopted to describe the characteristics of the cloud model as shown in Figure 2. These parameters could be obtained based on Equation (7) if the range of qualitative indexes was given: where k is the uncertainly cohesion degree of assessment indexes, which was assigned according to the ambiguity and randomness of each risk index. Cmax and Cmin are the maximum and minimum values of the level interval in the cloud model, respectively. Based on the current evaluation method, the indexes were qualitatively expressed using the normal distribution of the comprehensive assessment system. Meanwhile, the risk evaluation framework of roof accidents in a fully mechanized working face based on cloud model theory is shown in Figure 3. Based on the current evaluation method, the indexes were qualitatively expressed using the normal distribution of the comprehensive assessment system. Meanwhile, the risk evaluation framework of roof accidents in a fully mechanized working face based on cloud model theory is shown in Figure 3.

Establishment of Comprehensive Standard Cloud Model
The grade interval of each index as shown in Table 3 was dimensionlessly processed to change it to the range of 0 to 1 by using the extreme value method. Based on Equation (7), the standard value features of all 17 causing factors and 4 evaluation grades were obtained, and then the weight of each index was also calculated through the analytic hierar-

Establishment of Comprehensive Standard Cloud Model
The grade interval of each index as shown in Table 3 was dimensionlessly processed to change it to the range of 0 to 1 by using the extreme value method. Based on Equation (7), the standard value features of all 17 causing factors and 4 evaluation grades were obtained, and then the weight of each index was also calculated through the analytic hierarchy process. Furthermore, according to Equation (8), the standard value features (E x , E n and H e ) of four indexes in the criterion level and the standard value features (E x , E n and H e ) of four risk assessment grades (e.g., I, II, III and IV) are listed in Tables 7 and 8, respectively: where ω n is the weight of each index.  The forward cloud generator of the MATLAB program was adopted to set 5000 cloud drops to generate the standard cloud model of four indexes in the criterion level and all causing factors as shown in Figures 4 and 5, respectively.

Identification Criteria of Roof Risk Grade
Normally, the roof risk grade of indexes could be obtained through the comparison of the cloud drops of the roof risk assessment cloud images with the corresponding roof risk grade interval of the standard cloud images. However, the roof risk grade of indexes could not be intuitively determined if the cloud drops of the roof risk assessment cloud images fell into two roof risk grade intervals of the standard cloud images at the same time. Under this condition, the Dice coefficient was further introduced to evaluate the roof risk grade by using the entropy value (E n ) of the cloud images. According to Equation (9), the similarity degree between the roof risk assessment cloud images and the roof risk grade of the standard cloud images was obtained. Based on the maximum similarity principle, the corresponding roof risk assessment grade of the maximum Dice coefficient was regarded as the current roof risk grade in the working faces: where |A| and |B| are the interval length of a certain roof risk grade of the assessment cloud images and the standard cloud images, respectively, which were obtained through the entropy value (E n ). |A ∩ B| is the intersection length of a certain risk grade between the roof risk assessment cloud images and the roof risk standard cloud images.
The forward cloud generator of the MATLAB program was adopted to set 5000 cloud drops to generate the standard cloud model of four indexes in the criterion level and all causing factors as shown in Figures 4 and 5, respectively.

Identification Criteria of Roof Risk Grade
Normally, the roof risk grade of indexes could be obtained through the comparison of the cloud drops of the roof risk assessment cloud images with the corresponding roof risk grade interval of the standard cloud images. However, the roof risk grade of indexes could not be intuitively determined if the cloud drops of the roof risk assessment cloud images fell into two roof risk grade intervals of the standard cloud images at the same time. Under this condition, the Dice coefficient was further introduced to evaluate the roof risk grade by using the entropy value (En) of the cloud images. According to Equation (9), the similarity degree between the roof risk assessment cloud images and the roof risk grade of the standard cloud images was obtained. Based on the maximum similarity principle, the corresponding roof risk assessment grade of the maximum Dice coefficient was regarded as the current roof risk grade in the working faces: The forward cloud generator of the MATLAB program was adopted to set 5000 cloud drops to generate the standard cloud model of four indexes in the criterion level and all causing factors as shown in Figures 4 and 5, respectively.

Identification Criteria of Roof Risk Grade
Normally, the roof risk grade of indexes could be obtained through the comparison of the cloud drops of the roof risk assessment cloud images with the corresponding roof risk grade interval of the standard cloud images. However, the roof risk grade of indexes could not be intuitively determined if the cloud drops of the roof risk assessment cloud images fell into two roof risk grade intervals of the standard cloud images at the same time. Under this condition, the Dice coefficient was further introduced to evaluate the roof risk grade by using the entropy value (En) of the cloud images. According to Equation (9), the similarity degree between the roof risk assessment cloud images and the roof risk grade of the standard cloud images was obtained. Based on the maximum similarity principle, the corresponding roof risk assessment grade of the maximum Dice coefficient was regarded as the current roof risk grade in the working faces:

Engineering Background
Taking the 42,115 fully mechanized working faces in the Yushen mining area as the engineering background, the corresponding 4−2 coal seam of the working faces was relatively stable with an average thickness of 3.25 m. Additionally, the immediate roof of the working faces was siltstone and fine sandstone with moderate stability (Grade II). Meanwhile, the main roof of the working faces was siltstone and fine-grained sandstone with dense, hard and good stability. In addition, it was investigated that there was no fault in the mining area based on the geological exploring conditions and the exposure conditions of the mining roadway in the working faces. Moreover, it was revealed that the beams were pressed, and eight existing points were eroded by sandstone. According to the on-site measurements, the first weighting caving distance of the immediate roof was about 13-16 m. Additionally, the first and periodic weighting caving distances of the main roof were about 35-40 m and 18 m, respectively. The average moving distance between the roof and floor was about 0.1-0.25 m. This illustrated that the obvious weighting intensity was observed. On the other hand, the rated initial support force and the maximum average support strength of hydraulic supports in the working faces were 8728 kN and 994 kN/m 2 , respectively. The support was automatically moved according to the speed of the coal mining. The obtained actual overview conditions of all 17 risk evaluation indexes in the 42,115 fully mechanized working faces are shown in Table 9.

Calculating the Digital Features of Each Roof Risk Index
Firstly, each roof risk index was dimensionlessly processed, and the digital features of all indexes were obtained based on Equation (7) as shown in Table 10. According to Equation (8), combined with the index weight and digital features of each index, the comprehensive digital features of four indexes in the criterion level (e.g., support conditions, roof conditions, personnel working status and coal body status) are illustrated in Table 11.

Generating Comprehensive Risk Evaluation Cloud Images
According to the digital features of all risk evaluation indexes, 5000 cloud drops were set up in the standard cloud model to generate risk assessment cloud images as shown in Figures 6 and 7. Compared with the risk assessment cloud images and the standard cloud images, this directly showed that the cloud drops of the risk assessment cloud images were concentrated in the specific grade interval of standard cloud images. Meanwhile, based on the assessment cloud images and the entropy value (E n ) of the standard cloud images, the Dice coefficient was calculated as shown in Table 12 to effectively identify the risk grade of roof accidents in the 42,115 fully mechanized working faces.      As shown in Figure 6a, the risk assessment cloud image of support conditions intersected with Grade I and Grade II of the standard cloud model. However, the evaluation cloud drops were mainly concentrated in Grade I. On the other hand, through the calculation of the similarity degree, the evaluation was 0.6876 and 0.2789 for comparing the risk  As shown in Figure 6a, the risk assessment cloud image of support conditions intersected with Grade I and Grade II of the standard cloud model. However, the evaluation cloud drops were mainly concentrated in Grade I. On the other hand, through the calculation of the similarity degree, the evaluation was 0.6876 and 0.2789 for comparing the risk assessment cloud images with Grade I and Garde II of the standard cloud model, respectively. Therefore, the comprehensive risk grade of support conditions in the 42,115 fully mechanized working faces was regarded as Grade I (low risk). Figure 6b illustrates the risk assessment of roof conditions in the 42,115 fully mechanized working faces. It was observed that the risk assessment cloud image was intersected with Grade II and Grade III of the standard cloud model. Additionally, the entropy values (E n ) of Grade II and Grade III of the standard cloud model were same. It was only seen that the cloud drops of risk assessment were a little bit denser within the range of Grade III. Moreover, through the calculation of the similarity degree, the evaluation was 0.4539 and 0.4943 for comparing the risk assessment cloud images with Grade II and Grade III of the standard cloud model, respectively. Based on the maximum similarity principle, the risk grade of roof conditions was regarded as Grade III (large risk) in the 42,115 fully mechanized working faces.
As shown in Figure 6c, the cloud drops of risk assessment of personnel working status basically fell into the range of Grade I of the standard cloud model. Through the calculation of the similarity degree, the evaluation was 0 for comparing the risk assessment cloud image with Grade II, Grade III and Grade IV of the standard cloud model. This fully revealed that the risk grade of personnel working status was Grade I (low risk) in the 42,115 fully mechanized working faces. Figure 6d shows that the risk assessment cloud image of coal body status was intersected with Grade I and Grade II of the standard cloud model. Additionally, the evaluation cloud drops were mainly concentrated in Grade I. Through the calculation of the similarity degree, the evaluation was 0.44 and 0.2971 for comparing the risk assessment cloud images with Grade I and Grade II of the standard cloud model, respectively. Therefore, the risk grade of coal body status was Grade I (low risk) in the 42,115 fully mechanized working faces.
In Figure 7, it is intuitively illustrated that the comprehensive risk assessment cloud drops of all 17 indexes basically fell into the range of Grade II of the standard cloud model in the 42,115 fully mechanized working faces. In addition, through comparing the risk evaluation assessment cloud image with the standard could model image, the similarity degree of Grade II was 0.8606, while the corresponding similarity degree of Grade I, Grade III and Grade IV was 0. This explained that the risk grade of roof accidents in the 42,115 intelligent fully mechanized working faces was Grade II (general risk). In summary, support conditions, personnel working status and coal body status had limited effect on the occurrence of roof accidents, while the roof conditions in the working faces had a large risk in terms of the potential occurrence of roof accidents. Therefore, the management of roof conditions should be further strengthened in the whole mining processes of the 42,115 fully mechanized working faces.

Conclusions
A comprehensive risk evaluation cloud model of roof accidents was constructed based on the weight and standard grade level of each roof risk evaluation index. It effectively solved the problems of ambiguity and difficulty in quantitative expression and lacking visible results in the traditional risk evaluation methods of roof accidents. In terms of one specific engineering condition, the evaluation cloud images could be generated based on the newly proposed model, and then the Dice coefficient could be used to calculate the similarity degree between the risk evaluation cloud images and the risk grade interval of the standard cloud model. The qualitative and quantitative roof risk grade of the working faces could be effectively identified. Moreover, this model could also be widely used in the risk assessment process of various other potential accidents in coal mines. It provides theoretical guidance for effectively identifying, preventing and controlling the various accidents in coal mines and for ensuring the safe and efficient production of coal mines. Some conclusions are illustrated as follows.
(1) A comprehensive risk assessment model of roof accidents, including 4 indexes in the criterion level and a total of 17 causing factors, was constructed in a quantitative identification roof risk grade. Additionally, the main causing factors of roof accidents were roof separation distances, weighting intensity and rib spalling followed by the stress concentration of coal body, initial support force and geological conditions through the values of index weight determined by using the analytic hierarchy process. (2) Based on the comprehensive risk assessment cloud images, the overall risk grade in 42,115 intelligent fully mechanized working faces was regarded as Grade II (general risk). Additionally, the support conditions, personnel working status and coal body status had limited effect on the occurrence of roof accidents, while there was a large risk (Grade III) for roof conditions in the potential occurrence of roof accidents. (3) Compared with the AHP evaluation results, we showed that the risk assessment cloud model proposed in this study could identify the risk grade of each index and output comprehensive results by considering various aspects, such as cloud digital features, risk assessment cloud images and standard cloud images. It was effective in solving the problems of the ambiguity and difficulty of quantitative identification of traditional risk assessment methods, and it is worth promoting.