李 笛,陳 忠,鄧 彤
(1.云南力合礦山工程設(shè)計有限公司, 云南昆明 650093;2.昆明理工大學(xué)國土資源工程學(xué)院, 云南昆明 650093)
巖體質(zhì)量分類的DDA模型及應(yīng)用
李 笛1,陳 忠2,鄧 彤1
(1.云南力合礦山工程設(shè)計有限公司, 云南昆明 650093;2.昆明理工大學(xué)國土資源工程學(xué)院, 云南昆明 650093)
將距離判別分析法(DDA)應(yīng)用于巖體質(zhì)量等級判別與分類中,選取影響巖體質(zhì)量分類的主要因素作為判別因子,建立了巖體質(zhì)量分類評判的距離判別模型。以實測樣本作為訓(xùn)練樣本,建立相應(yīng)線性判別函數(shù)對待判樣本進(jìn)行分類。研究結(jié)果表明,距離判別分析模型是巖體質(zhì)量分類的一種有效方法,可以在實際工程中應(yīng)用。
巖體質(zhì)量;分類;距離判別分析
地下巖體工程的施工,必須確定巖體的穩(wěn)定性,而判別巖體是否穩(wěn)定必須要對周圍的巖體進(jìn)行分類[1]。目前常用的巖體穩(wěn)定性評價方法有以RMR方法和Q分類體系為代表的回歸分析法以及模糊數(shù)學(xué)方法、灰色理論方法、神經(jīng)網(wǎng)絡(luò)方法、物元層次分析法、支持向量機(jī)、物元可拓方法和可拓模糊評價方法等[2~12]。
本文借鑒馬氏距離判別分析法,選取影響圍巖分類的主要因素作為判別因子,建立馬氏距離判別分析模型。然后利用學(xué)習(xí)樣本構(gòu)建線性判別函數(shù),對巖體質(zhì)量進(jìn)行分級,分級結(jié)果取得了良好的效果,為巖體質(zhì)量分類研究提供了一條新思路。
設(shè)有 g個 p維總體 G1,G2,…,Gg,均值向量分別為 μ1,μ2,…,μg,協(xié)方差矩陣分別為∑1,∑2,…,∑g。類似兩總體的距離判別方法,計算樣本X到總體的距離。定義馬氏距離樣本X與總體G的馬氏距離為:

由(1)可知X到Gj和Gi的馬氏距離的平方差為:


則X到Gi的距離最小等價于對所有的j≠i,有Wij(x)>0,從而判別的標(biāo)準(zhǔn)為:
X∈Gi,若對一切

則 μ1,μ2,…,μg和∑未知時,可利用總體的訓(xùn)練樣本對其進(jìn)行估計、設(shè)為來自總體Gk的訓(xùn)練樣本(k=1,2,…,g),令:

利用Sk對∑的聯(lián)合估計為:


則判別準(zhǔn)則為:
X∈Gi,若對一切 j≠i,Wij(x) >0
為驗證判別預(yù)測準(zhǔn)則的優(yōu)良性[14],訓(xùn)練樣本構(gòu)建的模型,然后對學(xué)習(xí)樣本進(jìn)行回代,以誤判率η為閥值來考察評判的優(yōu)良性。
參考相關(guān)文獻(xiàn),結(jié)合實際情況,本文選取單軸抗壓強(qiáng)度R(x1)、體積節(jié)理數(shù)Jv(x2)、節(jié)理面風(fēng)化變異系數(shù)Ja(x3)、巖體聲波縱波速度Vp(x4)、透水性系數(shù)Wk(x5)、節(jié)理面粗糙度系數(shù)Jr(x6)作為巖體質(zhì)量分類的判別影響因子,建立線性判別模型。
以文獻(xiàn)[15]中工程樣本為例,選取其中的19個樣本數(shù)據(jù)進(jìn)行學(xué)習(xí),其余6個樣本數(shù)據(jù)作為待判樣本檢驗,具體見表1。根據(jù)實際情況,將巖體質(zhì)量分為3類,其判別函數(shù)如下:
W1,2=-13.8721x6-0.8589x5+0.0020303x4-0.006348x3-0.22109x2-0.12638x1+26.1118
W1,3=-18.93x6-11.09x5+0.007618x4+0.051598x3-0.24468x2-0.13602x1+29.4125
W2,3=-5.059x6-10.2322x5+0.005588x4+0.057946x3-0.023589x2-0.0096378x1+3.3008
將19個巖體質(zhì)量樣本作為訓(xùn)練樣本,其余的6個樣本作為學(xué)習(xí)樣本進(jìn)行訓(xùn)練,然后進(jìn)行回代。回代結(jié)果見表1。只有樣本13被誤判,誤判率η=1/25=0.04。由此可見,距離判別分析模型用于巖體質(zhì)量分類是高效、可行的。

表1 DDA分類結(jié)果
巖體質(zhì)量的分類受到眾多因素的影響,本文主要考慮影響巖體質(zhì)量的物理力學(xué)參數(shù),選取單軸抗壓強(qiáng)度、體積節(jié)理數(shù)、節(jié)理面風(fēng)化變異系數(shù)、巖體聲波縱波速度、透水性系數(shù)、節(jié)理面粗糙度系數(shù)6個影響因素作為巖體質(zhì)量分類的判別影響因子,以此構(gòu)建巖體質(zhì)量分類模型。將所建立的分類模型應(yīng)用到工程實例,取得了滿意的效果,說明該模型具有較強(qiáng)的判別能力,可應(yīng)用于巖體質(zhì)量分類。
[1] 馮夏庭,林韻梅.巖石力學(xué)與工程專家系統(tǒng)[M].沈陽:遼寧科學(xué)技術(shù)出版社,1992.
[2] 陳守煜.模糊最優(yōu)歸類理論模型及其在圍巖穩(wěn)定性分類與場地土類別評定中的應(yīng)用[J].水利學(xué)報,1993,(12):26-36.
[3] 費為進(jìn),經(jīng)蘇龍.多級模糊分級聚類法在地下工程圍巖穩(wěn)定性分類中的應(yīng)用[J].勘察科學(xué)技術(shù),2002,(1):20-24.
[4] 邵中勇,馮德順.公路隧洞圍巖穩(wěn)定性模糊評判方法[J].武漢理工大學(xué)學(xué)報(交通科學(xué)與工程版),2004,28(5):771-774.
[5] 韓桂武,劉 斌,范 鶴.模糊擇近原則在隧洞圍巖分類中的應(yīng)用[J].東北大學(xué)學(xué)報(自然科學(xué)版),2006,27(7):819-822.
[6] 馮玉國.灰色優(yōu)化理論模型在地下正程圍巖穩(wěn)定性分類中的應(yīng)用[J].巖土工程學(xué)報,1996,18(3):62-66.
[7] 王 彪,陳劍平,李鐘旭,等.人工神經(jīng)網(wǎng)絡(luò)在巖體質(zhì)量分類中的應(yīng)用[J].世界地質(zhì),2004,23(1):64-68.
[8] 周翠英,張亮,黃照藝.基于改進(jìn)BP網(wǎng)絡(luò)算法的隧洞圍巖分類[J].地球科學(xué)一中國地質(zhì)大學(xué)學(xué)報,2005,30(4):480-486.
[9] 朱曉寧,正惠勇.洞室圍巖開挖穩(wěn)定性等級綜合評判實用方法研究[J].中南公路工程,2005,30(2):48-51.
[10] 趙洪波,馮夏庭,尹順德.基于支持向量機(jī)的巖體工程分級[J].巖土力學(xué),2002,23(6):698-701.
[11] 原國紅,陳劍平,馬 琳.可拓評判方法在巖體質(zhì)量分類中的應(yīng)用[J].巖石力學(xué)與工程學(xué)報,2005,24(9):1539-1544.
[12] 王彥武.地下采礦工程巖體質(zhì)量可拓模糊評價方法[J].巖石力學(xué)與工程學(xué)報,2002,21(1):18-22.
[13] 高惠璇.應(yīng)用多元統(tǒng)計分析[M].北京:北京大學(xué)出版社,2005.
[14] 范金城,梅長林.數(shù)據(jù)分析[M].北京:科學(xué)出版社,2002.
[15] 王 彪,陳劍平,李鐘旭,等.人工神經(jīng)網(wǎng)絡(luò)在巖體質(zhì)量分類中的應(yīng)用[J].世界地質(zhì),2004,23(1):64-68.
2011-02-22)
李 迪(1971-),男,工程師,主要從事采礦方面生產(chǎn)和管理工作,Email:lidi71@163.com。