Trans. Nonferrous Met. Soc. China 22(2012) 425-431
Bayes discriminant analysis method to identify risky of complicated goaf in mines and its application
HU Yu-xi1, LI Xi-bing1, 2
1. School of Resources and Safety Engineering, Central South University, Changsha 410083, China;
2. Hunan Key Laboratory of Resources Exploitation and Hazard Control for Deep Metal Mines,Changsha 410083, China
Received 10 May 2011; accepted 20 October 2011
Abstract: A Bayes discriminant analysis method to identify the risky of complicated goaf in mines was presented. Nine factors influencing the stability of goaf risky, including uniaxial compressive strength of rock, elastic modulus of rock, rock quality designation (RQD), area ratio of pillar, ratio of width to height of pillar, depth of ore body, volume of goaf, dip of ore body and area of goaf, were selected as discriminant indexes in the stability analysis of goaf. The actual data of 40 goafs were used as training samples to establish a discriminant analysis model to identify the stability of goaf. The results show that this discriminant analysis model has high precision and misdiscriminant ratio is 0.025 in re-substitution process. The instability identification of a metal mine was distinguished by using this model and the identification result is identical with that of practical situation.
Key words: goaf; risky identification; Bayes discriminant analysis; metal mines
1 Introduction
In mining process of metal mines, a large number of underground goafs were brought about by using room and pillar method, a comprehensive law method and shrinkage method [1,2]. On the other hand, for most metal mines, such as Dachang mining, Luanchuan molybdenum, Changba zinc, Qinling gold and Kaiyang phosphate rock, there are also some underground goafs because of more than ten years non-governmental predatory exploitation. Statistics show that, in most mine accidents, goaf collapse is a very common incentive. And goaf has become one of main sources of harm, and is also risky in production [3,4]. Therefore, the risky identification of goaf is very important to ensure safety production of open and underground mining operations and to avoid the occurrence of major accidents in mines [5].
Many scholars have done a lot of work about risky identification. CHEN [6] has presented a monograph on the hazard identification, control and evaluation of the system discussed. And fuzzy comprehensive evaluation method (FCEM) [7], grey clustering analysis method (GCAM) [8,9] and artificial neural network (ANN) [10,11] theory have been used and got good effectiveness. Distance discriminant analysis method is a statistical analysis method based on observed characteristics (discriminant factor) of a certain number of samples and discriminant criterion, which had been used in mining and safety engineering recently [12-14]. However, for distance discriminant analysis method, the prior probability of each collectivity (Classification grade) is not taken into account and the difference of loss produced by mistake-discrimination is also neglected [11]. To solve the faults of distance discriminant analysis method, a Bayes discriminant analysis (BDA) model was presented to predict the stability of open pit slope in metal mines.
In the present work, in combination of the Bayes discriminant analysis theory and actual situation and stability factors of goaf, a Bayes discriminant analysis model is built and used in a practical engineering.
2 Bayes discriminant analysis method
2.1 Basic ideology of Bayes discriminant
Bayes discriminant is a probability discriminant analysis and various types of distribution density functions should be obtained before proceeding. The priori distribution was used to descript the level of awareness of the study collectivity before extracting samples, and then the posterior distribution was obtained by modifying the priori distribution based on extracted samples, followed by a variety of statistical inference. Following is a brief introduction and a variety of types of general overall Bayes discriminant [15].
2.2 Bayes discriminant of two normal collectivities
Suppose is a collectivity with p member indexes(Considering p indexes), and there are two collectivities, G1, G2, with distribution density functions, f1(X), f2(X). Suppose the priori distributions of G1 and G2 are
, (1)
with p1+p2=1. c(2|1) is the loss caused by misjudging G1 to G2, and c(2|1) is the loss caused by misjudging G2 to G1.
When for two normal collectivities, G1, G2, with c(2|1)= c(1|2), Bayes discriminant function can be expressed as:
, j=1, 2 (2)
Then, generalized squared distance function can be obtained as follows:
(3)
with μ1 and μ2 as mean vectors of G1, G2. And then posteriori probability function can be obtained:
(4)
Because ∝, then
(5)
Normally, μ1, μ2 and are unknown and their estimation values , and can be obtained from training samples, then
(6)
The estimation of posteriori probability function is
(7)
Bayes discriminant criterion can be expressed as:
(8)
2.3 Bayes discriminant of multi-normal collectivities
Suppose (k>2), and . Bayes discriminant function is expressed as Function (2).
Generalized squared distance function can be obtained as follows:
(9)
Posteriori probability function is
(10)
Then optimal division can be obtained as:
= (11)
μ1, μ2, …, μk and can be replaced with the expected values, and .
2.4 Evaluation of discriminant criterion
The prior probability pα is allocated by the proportion of training samples of collectivity Gα to all samples, i.e.,
, (12)
where pα is the prior probability of collectivity Gα, and nα is the number training samples belonged to collectivity Gα.
2.5 Evaluation of discriminant criterion
To estimate the reliability of discriminant criterion above, the re-substitution method was used to calculate the mis-discrimination rate [16]. All the training samples were regarded as the new samples and re-substituted into the discriminant criterions. The rate of misjudgment can be evaluated as the value of the number of mis-discrimination samples divided by the number of all samples.
3 BDA model for identification of goaf risky
3.1 Flow chart of building model
The process to build the Bayes discriminant model can be divided into five steps: 1) determining the impact factors influencing the goaf risky; 2) dividing risky levels of goafs; 3) building the BDA model by using training samples; 4) testing of BDA model; 5) application of BDA model. The flow chart of building model is shown in Fig. 1.
Fig. 1 Flow chart of building BDA model
3.2 Risky level dipartition of goaf
The characteristics of goaf disaster area and related studies show that the rock movement will not develop normally to the surface when the ore body buried deeper than 4 to 5 times the mining area. When the development of rock movement reached the surface, caving zone, fracture zone and bending zone (referred to as the "three zones") will form in upper rock mass of goaf (shown in Fig. 2). With the expansion of the scope of mining, the fracture zone and range of three zones will gradually change. When only the fracture zone occurs, the risk is low. When the fracture zone gradually develops, bending zone may occur, but the risk is also low; with the further development, the roof rock at footwall of fracture zone will collapse, which may cause a underground harm and can be regarded as a greater hazard. When caving zone appears and gradually expands, there are significant hazards. The serious risk of disaster is affected by many factors. Generally, caving zone and fracture zone have been extremely developed and goaf spans the space, resulting in complete collapse of bending zone, which will bring a great harm to the underground and surface of goaf.
Therefore, according to the severity of the dangerous consequences scale (Table 1), the instability of the goaf is divided into four risk levels: the first grade I (normal risk hazard), grade II (greater harm hazard), grade III (major hazard) and grade IV (large damage hazard), which means that the collectivity number of Bayesian discriminant analysis model has four.
Fig. 2 Schematic diagrams of goaf initial state (a) and goaf hazard state (b)
Table 1 Scale of risky of serious consequences
3.3 Determine of impact factors influencing goaf risky
The risky of goaf is influenced by many factors, which can be divided into several aspects, including rock properties of goaf, geological conditions, exploitation technical factors and treatment methods. After comprehensive analysis, nine specific indexes influencing the stability of goaf risky (shown in Table 2) [17], uniaxial compressive strength of rock (X1), elastic modulus of rock (X2), rock quality designation (X3), area ratio of pillar (X4), ratio of width to height of pillar (X5), depth of ore body (X6), volume of goaf (X7), dip of ore body (X8) and area of goaf (X9), were selected as discriminant indexes.
3.4 Training of learning samples and modeling
From Ref. [17], the actual data of 40 goafs were used as training samples to establish a Bayes analysis model to identify the risky of goaf. The BDA model is shown in Fig. 3.
Table 2 Original data of 40 training goaf samples
Fig. 3 Bayes discriminant model diagram
3.5 Data normalized
In the process of building BDA model, in order to make the model training more effective, the original sample data were normalized and the model input data will be in [0, 1] interval. For quantitative data, using the following formula:
(14)
where is normalized sample data; x is original sample data; xmin and xmax are the minimum and maximum values of original data, respectively. In addition, there is no conversion for pillar area value and rock quality indicators can be directly divided by 100. Normalized sample data are shown in Table 3.
3.6 Test of BDA model
The prior probability is allocated by the proportion of training samples, and then p1=10/40=0.250, p2=13/40=0.325, p3=10/40=0.250 and p4=7/40=0.175. The normalized data are input into BDA model and discriminant functions can be obtained. It can be seen from Table 1 that the forecasting risky level of goaf is the same as the actual status except sample No. 30, and the ratio of mis-distinguish is 1/40=0.025. It can be concluded that the BDA model can be applied to identifying the goaf risky.
Table 3 Input data of 40 training goaf samples and discriminant results
4 Actual engineering application of BDA model
Yaogangxian tungsten in Hunan province, China, is an old metal mine with more than mining 90 years, and there are a large number of goafs. The risky of goaf in Yaogangxian Tungsten Mine [4] was identified by using this BDA model above. The original data of mine goaf and discriminant results are shown in Table 4. The result of BDA model is identical with actual level, which also is the same with the result of ANN method [17]. It can be concluded that the BDA model can be applied in practical mine engineering to identifying the risky level. In fact, some goafs in the upper part of the ore has collapsed, making surface crack and deform. However, there is not catastrophic phenomenon. In the ore stope and tunnel of lower part, the phenomenon of partial roof collapse has occurred accidentally.
Therefore, the actual degree of danger is the level II, which is identical with the BDA model identification results. Table 5 also shows the results of 16 artificial neural networks (Artificial neural networks, referred to as ANN), which is also level II. Compared with artificial neural network method, the overall prior probability of every collectivity was fully considered in Bayesian discriminant analysis model. The Bayesian discriminant model has a fixed structure and the training process is simple and training is quick.
Table 4 Original data of Yaogangxian Tungsten Mine goaf
Table 5 Discriminant results of Yaogangxian Tungsten Mine goaf
5 Conclusions
1) Based on the Bayes discriminant analysis theory and actual characteristics of goaf risky, a Bayes discriminant analysis (BDA) model for instability identification of goaf risky was presented.
2) The results show that this discriminant analysis model has high precision and can be used in practical engineering. Compared with the other prediction methods, BDA model has a stable structure and the discriminant process is very simple and convenient.
3) It is the preliminary attempt that Bayes discriminant analysis theory is applied to analysis of the identification of goaf risky in mines. In the future work, it is necessary to conduct depth-study in selecting the study sample and discriminant genes, and then enhance the practicality of BDA model.
References
[1] ZHOU Hong-wei, XIE He-ping, ZUO Jian-ping. Developments in researches on mechanical behaviours of rocks under the condition of high ground pressure in the depths [J]. Advances in Mechanics, 2005, 35(1): 91-99. (in Chinese)
[2] FENG Chang-gen, LI Jun-ping, YU Wen-yuan, XUE Hua, LI Bao-dong. Mechanism study of abandoned stope disposal in Dongtongyu Gold Mine [J]. Gold, 2002, 23(10): 11-15. (in Chinese)
[3] State technological development programming of work safety-non- coal mine domain research bulletinl(2004-2010) [R]. Beijing: State Administration of Work Safety, State Administration of Coal Mine Safety, 2003. (in Chinese)
[4] LI Xi-bing, LI Di-yuan, ZHAO Guo-yan, ZHOU Zi-long, GONG Feng-qiang. Detecting and dealing approach of the underground goaf and safety evaluation in metal mines [J]. Journal of Mining and Safety Engineering, 2006, 23(1): 24-29. (in Chinese)
[5] LI Shan-cun. Study of hazard identification and risk assessment of underground goaf in metal mines [D]. Changsha: Central South University, 2007. (in Chinese)
[6] CHEN Bao-zhi. Identification, control and evaluation of hazards [M]. Chengdu: Sichuan Publishing House of Science & Technology, 1996. (in Chinese)
[7] SHI Bi-ming. A fuzzy cluster analysis on the danger of coal and gas burst from coal seams [J]. Journal of Huainan Mining Institute, 1994, 14(2): 38-43. (in Chinese)
[8] MUSEE N, LORENZEN L, ALDRICH C. New methodology for hazardous waste classification using fuzzy set theory: Part I. Knowledge acquisition [J]. Journal of Hazardous Materials, 2008, 154(1-3): 1040-1051.
[9] MUSEE N, ALDRICH C, LORENZEN L. New methodology for hazardous waste classification using fuzzy set theory: Part II. Intelligent decision support system [J]. Journal of Hazardous Materials, 2008, 157(1-3): 94-105.
[10] ZHONG Mao-hua, CHEN Bao-zhi. Study on gradation of major hazards based on neural network [J]. China Safety Science Journal, 1996, 6(S): 143-147. (in Chinese)
[11] ZHONG Mao-hua, CHEN Bao-zhi. Study in dynamic risk classification of major hazards based on neural network [J]. China Safety Science Journal, 1997, 7(2): 6-9. (in Chinese)
[12] GONG Feng-qiang, LI Xi-bing. Distance discriminant analysis method to the classification of engineering quality of rock masses [J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(1): 190-194. (in Chinese)
[13] LING Tong-hua, LIAO Yan-cheng, ZHANG Sheng. Safety criterion method of rock blasting vibration damnification based on multivariant discriminant [J]. Journal of Central South University: Science and Technology, 2010, 41(1): 322-327. (in Chinese)
[14] JIN Zhi-ren. Prediction of sand liquefaction based on distance discriminant analysis and its application [J]. Chinese Journal of Geotechnical Engineering, 2008, 30(5): 776-780. (in Chinese)
[15] FAN Jin-cheng, MEI Chang-lin. Data analysis [M]. Beijing: Science Press, 2002. (in Chinese)
[16] GAO Hui-xuan. Multivariate statistical analysis [M]. Beijing: Peking University Press, 2005. (in Chinese)
[17] LUO Yi-zhong. Major hazard source identification of widespread mined-out area instablity [D]. Changsha: Central South University, 2005. (in Chinese)
复杂采空区危险辨识的贝叶斯判别方法及应用
胡玉玺1, 李夕兵1, 2
1. 中南大学 资源与安全工程学院,长沙 410083;
2. 深部金属矿产开发与灾害控制湖南省重点试验室,长沙 410083
摘 要:提出了复杂采空区危险程度辨识的贝叶斯判别方法。基于多元判别分析理论,将贝叶斯判别方法应用于金属矿山采空区危险程度的预测判别问题中,建立了相应的贝叶斯判别分析模型。该模型选用岩石单轴抗压强度、岩石弹性模量、岩石质量指标、矿柱面积比率、矿柱宽高比、矿体埋藏深度、采空区体积、矿体倾角和采空区面积9项指标作为判别因子,将采空区的危险性等级分为4级;以40个采空区实测数据作为学习样本进行训练,建立相应判别函数对待判样本进行分类。研究结果表明,贝叶斯判别模型的学习精度很高,回判估计的误判率为0.025。利用学习后的模型对某金属矿山采空区实例进行了稳定性判别,判别结果和实际情况相符。
关键词:采空区;危险辨识;贝叶斯判别分析;金属矿山
(Edited by YANG Hua)
Foundation item: Project (2010CB732004) supported by the National Basic Research Program of China
Corresponding author: HU Yu-xi; Tel: +86-13975891760; E-mail: huyuxi@csu.edu.cn
DOI: 10.1016/S1003-6326(11)61194-1