侯恩科 閆會 鄭永飛 楊帆

摘要：風化基巖含水層是對礦井安全生產構成威脅的主要含水層之一，準確預測其富水性是煤礦水害防治的基礎。針對風化基巖富水性預測問題，以張家峁井田為例，在總結風化基巖巖性特征、分布特征和富水性特征的基礎上，以風化基巖厚度、風化程度、巖性組合指數、脆塑性巖厚度比和地形地貌作為判別指標，以研究區(qū)及其外圍29組風化基巖抽水試驗數據作為訓練及檢驗樣本，建立了風化基巖含水層富水性B aves判別模型，可對無抽水試驗資料區(qū)域的風化基巖富水性類型進行預測。研究表明，張家峁井田風化基巖整體富水性較弱，且空間分布不均。在井田中部和西北角局部區(qū)域為強富水性，中部大部分區(qū)域為中等富水性，東部及東北部區(qū)域幾乎全為極弱富水性，其余部分為弱富水性。通過與礦井生產過程中井下涌水量觀測數據對比，顯示所建判別模型具有較高的準確率，表明運用此模型進行風化基巖富水性預測是可行的，為風化基巖含水層的富水性預測提供了一種新方法。

關鍵詞：風化基巖;Bayes判別;富水性;張家峁井田

中圖分類號：P 641.4;TD 741

文獻標志碼：A

文章編號：1672 -9315（2019）06 -0942 -08

DOI：10. 13800/j. cnki. xakjdxxb.2019. 0604

開放科學（資源服務）標識碼（OSID）：

收稿日期：2019 - 04 - 15

責任編輯：李克永

基金項目：國家自然科學基金面上項目（ 41472234）;陜西煤業(yè)化工集團科研項目（2017 SMHKJ-C-23，2018 SMHKI-C-31）

通信作者：侯恩科（1963 -），男，陜西扶風人，教授，博士生導師，E-mail：houek@ xust. edu. cn

O 引言

侏羅紀煤炭是中國的主要煤炭資源，也是中國煤炭開發(fā)的主要對象之一。風化基巖是巖石遭受物理風化、化學風化和生物風化等作用而形成的產物，與正常巖石相比，其在顏色、結構構造、孔隙率和含水率等方面均有顯著差異[1-2]。風化基巖在陜北侏羅紀煤田普遍發(fā)育，近年來在煤礦生產中發(fā)現(xiàn)侏羅系風化基巖水是造成礦井水害的重要水源之一，對井下安全生產構成嚴重威脅[3-4]。侏羅系風化基巖水害具有突水量大、持續(xù)時間長、防治難度較大等特點。但是由風化基巖含水層所引起的水害并不是在所有地方都會發(fā)生，而是在富水性強的地方多有發(fā)生。因此，在侏羅系風化基巖水害防治中，準確預測其富水性空間變化具有重要意義[5]。

目前針對含水層富水性探測評價的方法，主要有地球物理探測方法和數學地質分析方法等，常用的綜合物探法[6]、瞬變電磁法以及直流電法[7]等地球物理探測方法應用廣泛，但存在費用較高、采集數據粗糙、準確性較低等不足[8];數學地質分析有層次分析法[9 -10]、聚類分析法[11]、模糊數學[12]等方法，有些方法考慮的影響因素相對較少，且確定指標權重的主觀性較強，往往受到人為因素的干擾，在不同程度上影響了富水性預測的客觀性。Logistic回歸分析模型和證據權模型對各主控因素的綜合作用考慮相對不足[13];基于多元統(tǒng)計的距離判別[14]和費舍爾判別[15]等方法則未考慮判別方法和各總體出現(xiàn)概率的關系，對判別方法和誤判損失也缺乏考量[16]。Bayes判別分析法具有方法簡單、求解快速和準確性高的特點[19]。因此，文中以張家峁煤礦為例，探討利用Baves判別法進行風化基巖的富水性預測，以期為風化基巖富水性預測提供一種新方法。

1 風化基巖特征

1.1 風化基巖巖性特征

風化基巖是巖石遭受風化作用而形成的產物，與正常巖石相比，風化基巖在巖性方面有獨特的性質。風化后的巖石，裂隙發(fā)育，破碎嚴重，且隨著風化程度的提高，其巖石抗壓強度逐漸下降，孔隙率逐漸增大，煤層采動影響下，內摩擦角也會增大[17]。砂質泥巖、泥質粉砂巖風化后趨于粘土化，泥巖在風化嚴重時容易泥化[18]。區(qū)內風化基巖大部分為延安組砂巖、泥巖風化而成，只有西北部有少量風化基巖是直羅組砂巖、泥巖層的風化產物。風化基巖原巖巖性主要為中粒砂巖、細粒砂巖、粉砂巖以及少量粗砂巖和泥巖。泥巖、粉砂巖等塑性較強的巖石一般風化程度較嚴重，而中粗粒砂巖等硬脆性巖石的風化程度相對較弱。

1.2 風化基巖分布特征

區(qū)內風化基巖廣泛發(fā)育，風化程度一般由上到下逐漸減弱。風化基巖厚度與地形有一定的相關性，溝谷處的風化基巖一般厚度較小，梁峁等地勢較高處的風化基巖厚度較大。區(qū)內風化基巖厚度分布不均，據鉆孔資料分析，風化基巖最厚達44. 91 m，最薄僅1.6 m（圖1）。

1.3 風化基巖富水性特征

風化基巖風化裂隙較發(fā)育，含水性條件相對較好[19]。在上覆風積沙層和薩拉烏蘇組地層存在而新近系保德組隔水層缺失的地區(qū)，由于接受潛水補給，富水性相對較好。風化基巖的富水性與風化基巖的巖性、厚度、風化程度以及補給條件等因素有關，且富水性分布不均。

2 Bayes判別方法

2.1

Bayes判別原理

判別分析是根據樣本的總體特征，按照一定的判別準則對新樣本的歸屬做出判斷[20-22]。與其他判別方式不同，B aves判別對先驗概率十分注重。在對某個總體進行判別之前，對所研究的總體已經有一定的認識，實際中表現(xiàn)為某些總體出現(xiàn)的可能性大于其他總體。在獲得先驗概率之后，通過隨機抽取的樣本對其進行修正，從而得到后驗概率?；诤篁灨怕始纯蛇M行各種統(tǒng)計推斷，并考慮誤判造成的損失。基于Baves統(tǒng)計思想的判別分析即為B aves判別法。

2.1.1 先驗概率

設有k個p維總體G1，G2，…，G，樣品來自各個總體，其概率密度函數

3 風化基巖富水性Bayes判別模型

3.1

Bayes判別指標的確定

風化基巖含水層的富水性規(guī)律和空間分布特征受諸多因素的影響，依據井田及其外圍施工的風化基巖抽水試驗鉆孔資料，通過比較各鉆孔信息，發(fā)現(xiàn)風化基巖的風化程度、巖性等條件不同，單位涌水量有較大差別，參考前人相關研究成果，挑選出風化基巖厚度、風化程度、巖性組合指數、脆塑性巖厚度比和地形地貌5個對富水性影響較大的因素作為Baves判別指標[23 -24]。

3.1.1 風化基巖厚度

風化基巖厚度的大小，決定了含水層儲水空間的大小。一般而言，遭受風化的巖石，其內部孔隙增多，節(jié)理裂隙相對發(fā)育，透水性也明顯增強。風化基巖厚度越大，其富水性也越強。如在均為中等風化中粒砂巖的情況下，風化厚度為26.3 m的8 - HB6孔（圖2（g））單位涌水量為0. 473 9 L/（s·m），而風化厚度為9.1 m的8- HB3孑L（圖2（k）0單位涌水量為0.246 9 L/（s·m）。同樣情況的還有H05孔（圖2（a））和Hl7孑L（圖2（b）），在巖層結構相似的情況下，風化厚度較大的H17孔富水性強于風化基巖厚度相對較小的H05孔。

3.1.2 風化程度

就研究區(qū)而言，巖石遭受風化的程度越強.其巖體結構就越疏松，孔隙也就越發(fā)育，富水性就越強。將研究區(qū)風化基巖的風化程度劃分為強、中等和弱3類，并量化賦值（表1）。

如在均為粗粒砂巖厚度接近的情況下，強風化的HB2 -8孑L（圖2（f））單位涌水量為3. 775 6L/（s·m），而中等風化的8- HB6孔（圖2（g））單位涌水量為0.473 9 L/（s·m）。

3.1.3 巖性組合指數

風化基巖通常由多層巖性不同的巖石組成，巖性不同其富水性也有差異。據抽水試驗數據，以硬脆性中粗粒砂巖為主的風化基巖，孔隙度較大，儲水能力強，滲透性好，其富水性一般較強;以軟塑性泥巖、砂質泥巖為主的風化基巖孔隙度較小，儲水能力較弱，滲透性差，其富水性一般較弱。如在均為中等風化的條件下，風化基巖厚度為9.1m巖性為中粒砂巖的8 - HB3孔（圖2（k））單位涌水量為0.246 9 L/（s.m），大于風化基巖厚度為16. 37 m巖性為泥巖的8- HB9孔（圖2（j））0. 08L/（s.m）的單位涌水量。同為風化砂巖，粒徑較大的粗砂巖富水性強于粒徑較小的細砂巖，如H21孑L（圖2（（c））的單位涌水量0.006 579 L/（s.m）明顯小于HB2 -8（圖2（f））的單位涌水量3. 775 6L/（s.m）。在風化基巖厚度接近的情況下，巖層中砂巖占比多的富水性要強于泥巖占比多的，如風檢3號孔（圖2（d））巖性組合中泥巖占比50. 77%，其單位涌水量為0.000 52 L/（s.m），小于巖性組合全為砂巖的BK1孔（圖2（e））的單位涌水量0.005 591 L/（s.m）。

因此，將風化基巖巖性和其在風化基巖中的厚度占比結合起來，構建風化基巖巖性組合指數F，其數值越大，表明富水性越強（表2）。

3.1.4 脆塑性巖厚度比

在地應力的作用下，地下巖層會發(fā)生變形甚至破壞，不同巖性的巖石在地應力作用下發(fā)生變形破壞的形式不同。硬脆性的中粗粒砂巖在受力時多發(fā)生脆性破壞，以破裂的形式釋放應力，在巖石內部產生大量裂隙，從而具有較強的富水性;而軟塑性的泥巖、粉砂巖等，在受力時往往產生塑性變形，以塑性形變的方式釋放應力，因此巖體內部結構相對完整，裂隙發(fā)育較少，富水性也相應較弱。因此，可用脆性巖與塑性巖厚度之比作為反映風化基巖富水性的一個指標，通常比值越大富水性越強。如在風化程度及厚度均接近的情況下，脆塑性巖厚度比小的風檢3號孔（圖2（d））單位涌水量0. 000 52 11（s.m）小于脆塑性巖厚度比大的BK3孑L（圖2（1））的單位涌水量0. 004 435L/（s.m）。

3.1.5 地形地貌

地形地貌條件影響著風化基巖含水層的補給。研究區(qū)西部多為風沙灘地，地表覆蓋著松散沙層，且地形平坦，因此降雨入滲強烈，補給條件較好，富水性強。中部常家溝、烏蘭不拉溝等溝谷地帶，補給條件相對較好，富水性也相對較好。常家溝東北側主要為黃土丘陵溝壑區(qū)，降雨易于流失，不易入滲，滲透率較低，因此補給條件較差，富水性也相對較弱。據此，對不同地貌條件予以量化賦值（表3）。

3.2

Bayes判別模型的建立及驗證

選擇風化基巖厚度（Xl）、風化程度（X2）、巖性組合指數（X3）、脆塑性巖厚度比（X4）和地形地貌因素（x5）作為Baves判別指標，以張家峁井田及其外圍有風化基巖抽水試驗資料的29組完整鉆孔數據作為訓練及驗證樣本，隨機選取其中21組數據用來訓練樣本，其余8組用來檢驗判別結果的準確性。采用SPSS統(tǒng)計分析軟件進行Baves判別訓練，先驗概率采用“根據組大小計算”的原則確定[12]。

富水性分區(qū)預測時富水性等級通常根據鉆孔單位涌水量的大小，以《煤礦防治水細則》（以下簡稱《細則》）為標準進行確定。《細則》中將含水層富水性分為4級：弱富水性（q≤0.1 L/s.m）、中等富水性（0.15 L/s.m）。對井田內及其外圍有抽水試驗的29個鉆孔分析發(fā)現(xiàn)，其最大單位涌水量為3. 776 L/s.m，屬于強富水性，未見富水性極強的鉆孑L，其中富水性弱的鉆孑L占58.6%.為了使富水性預測結果更便于礦井防治水工作的開展，將富水性弱這一級別進行進一步細分，增加富水性極弱這一級別，并去掉井田內未見的富水性極強這一級別，從而將研究區(qū)內風化基巖的富水性分為極弱、弱、中等和強富水性4個等級（表4）。

以這4個等級作為Baves判別的4個總體，且假定各總體協(xié)方差矩陣相等，建立風化基巖富水性等級預測的Baves判別模型。運用SPSS軟件對隨機選取的21組訓練樣本進行統(tǒng)計分析，得到Baves線性判別函數式[23]。將訓練結果與實際情況對比，21組訓練樣本中有17組判別結果與實際情況相一致，判別準確率為80.95%（表5）。

訓練所得Baves線性判別函數如下

Y1 =0. 031X， -0.238X2 +2. 813X3+8.244X4+7. 314Xs - 19. 759

Y2 =0. 177X， -0.243X2 +3. 810X3+8.402X4+7. 438Xs - 24. 073

Y3 =0. 043X， -0. 314X2 +3. 917X3+ 12. 352X4+ 11. 192Xs - 42. 680

Y4=0. 02IXl -0. 394X2+3.869X3+ 14. 745X4+ 12. 932X5 - 57. 030

將其余8組數據作為驗證樣本，帶入判別函數進一步驗證判別模型的準確性。8組驗證樣本中，除了BKI鉆孔高判了一級外，其余7組樣本的判別結果均與實際情況一致，判別準確率達87. 5%（表6）。從整體結果看，由樣本訓練所得Baves判別函數模型準確率相對較高，可用于對未進行抽水試驗的探煤鉆孑L進行富水性預測。

4 風化基巖富水性分區(qū)預測

對于井田內其他80個未進行抽水試驗但信息完整的探煤鉆孔，按照訓練樣本的量化賦值規(guī)則進行取值，帶人Bayes判別函數得到對應鉆孔的富水性預測類型，再采用插值軟件將所有鉆孔信息進行插值運算，生成井田風化基巖富水性分區(qū)預測圖（圖3）。

從預測圖可以看出，井田內風化基巖整體富水性較弱，且空間分布不均。在井田中部和西北角局部區(qū)域為強富水性，強富水區(qū)分布范圍較小，中部大部分區(qū)域為中等富水性，東部及東北部區(qū)域幾乎全為極弱富水性，其余部分為弱富水性。

結合礦井已開采工作面涌水量生產觀測數據，烏蘭不拉溝北部源頭處布置有22201綜采工作面，該工作面的涌水量在8.6 - 33.6 m3/h，且該涌水量中有一部分來源于2。2煤燒變巖含水層，與該區(qū)風化基巖含水層極弱富水性預測相一致;在上房溝流域附近的23102綜采工作面涌水量則只有4.7 rri3/h，與該區(qū)的極弱一弱富水性預測結果相符;井田東部及東北部5。2煤層15201 - 15208和N15201 - N15203工作面開采后，全井田的井下正常涌水量在105. 97 m3/h，涌水量較小，說明開采區(qū)富水性弱，與此區(qū)域的極弱一弱富水性預測結果吻合。

5 結論

1）通過已有水文地質數據，建立起風化基巖富水性Bayes判別模型，對未進行抽水試驗的區(qū)域進行預測，結果表明，井田中部區(qū)域風化基巖富水性相對較強，東部和北部富水性相對較弱。

2）利用Bayes判別分析法建立的風化基巖富水性判別模型的預測結果與實際觀測結果～致，說明利用Baves判別模型對無風化基巖抽水試驗資料區(qū)域進行富水性預測是可行的。

參考文獻（ References）：

[1]李定龍，周治安，臨渙礦區(qū)煤系風化帶及其水文工程地質意義的初步研究[J].中國煤炭地質，1994（1）： 49 - 54.

LI Ding-long，ZHOU Zhi-an.A preliminary study on theWeathering zone and ifs hydrogeological and engineeringgeological significance in Linhuan Mining Area[J].Coal Geology of China，1994，6（1）：49 - 54.

[2] Woithington S R H，Davies G J，Alexander E C.En-hancement of bedrock permeability by weathering[J].Earth Science Reviews， 2016， 160： 188 - 202.

[3]侯恩科，童仁劍，王蘇健，等.陜北侏羅紀煤田風化基巖富水性Fisher模型預測方法[J].煤炭學報，2016，41（9）：2312 -2318.

HOU En-ke， TONG Ren-jian， WANC Su-jian，et al.Pre-diction method for the water enrichment of weatheredbedrock based on Fisher model in Northern Shaaxi Ju-rassic coalfield[J].Journal of China Coal Society， 2016，41（9）：2312 - 2318.

[4]侯恩科，車曉陽，馮潔，等，榆神府礦區(qū)含水層富水特征及保水采煤途徑[J].煤炭學報，2019，44（ 3）： 812 - 819.

HOU En-ke， CHE Xiao-yang， FENG Jie， et al.Abun-dance of aquifers in Yushenfu coal field and the meas-ures for water-preserved coal mining[J].Joumal of Chi-na Coal Society ，2019 ，44.

[5]靳德武，周振方，趙春虎，等，西部淺埋煤層開采頂板含水層水量損失動力學過程特征[J].煤炭學報，2019 ，44（3）：30 -40.

JIN De-wu， ZHOU Zhen-fang， ZHAO Chun-hu， et al.Dynamics process analysis of groundwater quantity lossof roof aquifer with shallow seam mining in Western Chi-na[J].Journal of China Coal Society ，2019 ，44（3）：690- 700.

[6]離俊良，段建華，郭粵蓮.綜合物探技術在探測煤礦采空區(qū)及其富水性中的應用[J].中國煤炭地質，2007 ，19（ A02）：111 - 113.

GAO Jun-liang， DUAN Jian-hua， GUO Yue-Iian. Appli-cation of comprehensive geophysical exploration technol-ogy in detecting goaf and its water enrichment in coalmine[J].Coal Geology of China， 2007， 19（ A02）： 111- 113.

[7]施龍青，翟培合，魏久傳，等，三維高密度電法技術在巖層富水性探測中的應用[J].山東科技大學學報（自然科學版），2008 ，27（6）：1-4.

SHI Long-qing ，ZHAI Pei-he， WEI Jiu-chuan，et al.Ap-plication of three-dimensional high-density electricmethod in detecting water-rich strata[ Jl. Journal ofShandong University of Science and Technology（ Natural Science Edition） ，2008 ，27（6）：1-4.

[8]韉超，泮曉華，李國梁，等.基于GIS多源信息集成的含水層富水性模糊層次分析法[J].水文地質工程地質，2012 ，39（4）：19 - 25.

HAN Chao， PAN Xiao-hua， LI Cuo-Iiang， et al.Fuzzyanalytic hierarchy process for aquifer water enrichmentbased on GIS multi-source Information integration[J].Hydrogeologic Engineering Geology， 2012， 39（ 4）： 19 - 25.

[9] 武強，樊振麗，劉守強，等.基于GIS的信息融合型含水層富水牲評價方法——富水性指數法[J].煤炭學報，2011 ，36（7）：1124 - 1128.

WU Qiang， FAN Zhen-li， LIU Shou-qiang， et al.Waterenrichment index method of information fusion aquiferbased on GIS[J].Journal of Coal，2011 ，36（7）：1124 - 1128.

[10] Zeynep Didem Unutmaz Durmu pglu. Assessment oftechno-entrepreneurship projects by using Analytical Hi-erarchy Process（ AHP）[J].Technology in Society，2018： S016079IX17302749.

[11]李青松，李紹泉，韓真理.基于模糊聚類分析法的近距離突出煤層群突出危險程度評價[J].煤礦安全， （11）.

LI Qing-song， LI Shao-quan， HAN Zhen-li. Coal and gasoutburst hazard degree evaluation of close distance out-burst coal bed group based on fuzzy cluster analysismethod[J].Safety in Coal Mines，（ II）.

[12]劉偉韜，李加祥，張文泉.頂板涌水等級評價的模糊數學方法[J].煤炭學報，2001，26（4）：399 -403.

LIU Wei-tao， LI Jia-xiang， ZHANC Wen-quan. Fuzzymathematical method for evaluation of roof water gushinggrade[J].Journal of Coal，2001 ，26（ 4）：399 - 403.

[13]武強，張波，趙文德，等.煤層底板突水評價的新型實用方法V-基于GIS的ANN型、證據權型、

Logistic回歸型脆弱性指數法的比較[J].煤炭學報，2013 ，38 （1）：21 - 26.

WU Qiang， ZHANG Bo， ZHAO Wen-de， et al.A newpractical method for evaluating water inrush from coalseam floor V：comparison of ANN type， evidence weighttype and Logistic regression vulnerability index methodbased on GIS[J].Journal of Coal， 2013 ，38（1）：21 - 26.

[14]宮鳳強，劉科偉，李志國，礦區(qū)采空塌陷危險性預測的Bayes判別分析法[J].采礦與安全工程學報，2010 ，27 （1）：34 - 38 ，43.

GONC Feng-qiang， LIU Ke-wei， LI Zhi-guo. Bayes dis-criminant analysis method for prediction of mining col-lapse risk in mining area[ Jl. Journal of Mining andSafety Engineering ，2010 ，27（1）：34 - 38 ，43.

[15]鐘世勇.基于費舍爾判別分析的半監(jiān)督故障分類方法研究[D].浙江：浙江大學，2015.

[16]李春萍，郝會兵.煤巷圍巖分類的Bayes判別分析法

[J].煤炭學報，2011（ S2）：304 - 307.

LI Chun-ping， HAO Hui-bing. Bayes discriminant analy-sis of surrounding rock classification in coal roadway[J]. Journal of Coal，2011（ S2）：304 - 307.

[17]姚多喜，成春奇.淮北百善煤礦基巖風化帶的隔水性研究[J].煤田地質與勘探，1997（2）：38 -39.

YAO Duo-xi， CHENG Chun-qi. Study on water repellen-cy of weathering zone of bedrock in baishan coal mine，Huaibei[J].Coalfield Ceology and Exploration， 1997（2）：38 -39.

[18]張楓林，孫云飛，風化基巖斷層破碎帶巷道注漿加固技術研究[J].煤炭科技，2011（2）：8 -9.

ZHANG Feng-Iin， SUN Yun-fei. Study on grouting rein-forcement technology in broken zone of weathered bed-rock[J].Coal Science and Technology ，2011（2）：8-9.

[19] Contributions of water supply from the weathered bed-rock zone to forest soil quality[J].Geoderma ，2003， 114（3 -4）：0 -400.

[20]文暢平.基于Bayes判別分析法的沖擊地壓預測與危險性分級[J].自然災害學報，2015 （5）：229 - 236.

WEN Chang-ping. Prediction and hazard classification ofrock burst based on Bayes discriminant analysis[J].Journal of Natural Disasters ，2015（5）：229 - 236.

[21]史秀志，周 健，鄭 緯，等.邊坡穩(wěn)定性預測的Bayes判別分析方法及應用[J].工程科學與技術，2010 ，42（3）：63 - 68.

SHI Xiu-zhi， ZHOU Jian， ZHENG Wei， et al.Bayes dis-criminant analysis method for slope stability predictionand its application[J].Engineering Science and Tech-nology，2010，42（3）：63 - 68.

[22]譚杰.基于Bayes判別分析的配電網故障選線方法[J].電力學報，2016（3）：217 - 223.

TAN Jie. Fault line selection method for distribution net-work based on Bayes discriminant analysis[J].Journalof Electric Power， 2016（3）：217 - 223.

[23]武強，許珂，張維.再論煤層頂板涌（突）水危險性預測評價的“三圖一雙預測法”[J].煤炭學報，2016 ，41（6）：1341 - 1347.

WU Qiang， XU Ke， ZHANG Wei. Further research on“three maps-two predictions" method for prediction oncoal seam roof water bursting risk[J].Joumal of ChinaCoal Society ，2016 ，41（6）：1341 - 1347.

[24]曾一凡，李哲，宮厚建，等.頂板風化基巖含水層富

水特征與涌（突）水危險性預測[J].煤炭工程，2018.ZENG Yi-fan， LI Zhe， CONG Hou-jian， et al.Water en-

richment characteristics and water inrush risk predictionof roof weathered bedrock aquifer[J].Coal Engineering，2018.

[25]付玉華，董隴軍.巖爆預測的Bayes判別模型及應用[J].中國礦業(yè)大學學報，2009，38（4）：528 -533.

FU Yu-hua， DONG Long-jun. Bayes discriminant modelfor rock burst prediction and its application[J].Journalof China University of Mining and Technology ，2009 ，38（4）：528 -533.