肖屈日，趙國(guó)彥，劉 建，簡(jiǎn) 箏

(1.中南大學(xué)資源與安全工程學(xué)院，湖南 長(zhǎng)沙 410083；2.長(zhǎng)沙礦山研究院有限責(zé)任公司，湖南 長(zhǎng)沙 410012)

0 引言

在地下采礦工程中，為了安全和便利通常會(huì)有意保留部分礦體作為支柱，礦柱的合理留設(shè)不僅有支撐頂板的作用，還能有效地保護(hù)巷道和地面建筑物[1]。隨著采掘工作推進(jìn)、開(kāi)采向深部延伸，頂板暴露面積慢慢變大、礦柱應(yīng)力逐漸增加，礦柱失穩(wěn)的風(fēng)險(xiǎn)也愈發(fā)嚴(yán)峻。作為影響采空區(qū)穩(wěn)定的關(guān)鍵結(jié)構(gòu)單元[2]，礦柱失穩(wěn)可能會(huì)導(dǎo)致柱體不同程度的自由面巖石剝落、剪切破壞、內(nèi)部劈裂和地質(zhì)結(jié)構(gòu)面的滑移與潰曲破壞[3]，更甚者將引發(fā)冒頂和空區(qū)塌陷事故[4]，給礦山企業(yè)(特別是房柱法開(kāi)采[5]礦山企業(yè))人員和財(cái)產(chǎn)安全造成巨大的威脅。因此，開(kāi)展關(guān)于礦柱穩(wěn)定性的研究對(duì)于地下礦山高效開(kāi)采、災(zāi)害防治具有重要現(xiàn)實(shí)意義。

國(guó)內(nèi)外學(xué)者針對(duì)礦柱的設(shè)計(jì)與穩(wěn)定分析問(wèn)題提出了許多理論方法，如：面積承載理論[6-7]、經(jīng)驗(yàn)公式法[7-8]、數(shù)值模擬技術(shù)[4-5,9]等。傳統(tǒng)的礦柱穩(wěn)定性分析方法主要根據(jù)這些理論與手段得到礦柱強(qiáng)度、應(yīng)力等數(shù)據(jù)，以礦柱強(qiáng)度與應(yīng)力之比來(lái)衡量礦柱的安全性。盡管這些方法已經(jīng)發(fā)展的十分成熟，但導(dǎo)致礦柱失穩(wěn)的因素眾多，單一經(jīng)驗(yàn)指標(biāo)評(píng)判存在可靠性低、判別效果差的問(wèn)題。機(jī)器學(xué)習(xí)方法以多特征指標(biāo)為輸入，將比較評(píng)判過(guò)程無(wú)形納入學(xué)習(xí)算法的“黑匣子”中，并借助大量的實(shí)例數(shù)據(jù)進(jìn)行監(jiān)督學(xué)習(xí)，使得分析更為可靠和高效，近年來(lái)被廣泛應(yīng)用于礦柱設(shè)計(jì)和穩(wěn)定性評(píng)估中[10-12]。本質(zhì)上，同經(jīng)驗(yàn)公式法一樣，機(jī)器學(xué)習(xí)算法也是一種基于數(shù)據(jù)的統(tǒng)計(jì)方法，但機(jī)器學(xué)習(xí)算法能夠借助計(jì)算機(jī)的高速計(jì)算能力總結(jié)出人類暫時(shí)無(wú)法歸納出的“經(jīng)驗(yàn)”。通過(guò)機(jī)器學(xué)習(xí)算法逆構(gòu)礦柱穩(wěn)定性影響因素與礦柱狀態(tài)之間的非線性映射關(guān)系來(lái)判定未知礦柱的穩(wěn)定性，對(duì)指導(dǎo)工程實(shí)踐有著積極而深刻的作用。

旋轉(zhuǎn)森林算法是由JUAN J RODRIGUEZ等[13]于2006年提出的新型集成學(xué)習(xí)算法。與隨機(jī)森林等經(jīng)典集成學(xué)習(xí)算法相比，旋轉(zhuǎn)森林算法能夠生成差異度大且精度較高的基分類器，具有更好的泛化能力，近年來(lái)被廣泛應(yīng)用于生物醫(yī)學(xué)[14-15]和模式識(shí)別[16]等諸多領(lǐng)域。目前，尚未有關(guān)旋轉(zhuǎn)森林算法在礦柱穩(wěn)定分析方面的應(yīng)用研究。因此，本文在概述旋轉(zhuǎn)森林相關(guān)理論研究的基礎(chǔ)上，將其引入到礦柱穩(wěn)定性分析中，構(gòu)建出一種基于ICA-RoF算法的硬巖礦柱穩(wěn)定性判別模型，以期為地下礦山的礦柱設(shè)計(jì)與安全評(píng)價(jià)工作提供參考。

1 ICA-RoF理論與算法

1.1 旋轉(zhuǎn)森林算法

旋轉(zhuǎn)森林算法[13]通過(guò)對(duì)原始樣本特征進(jìn)行處理，并采用一定的特征變換方法獲取訓(xùn)練基分類器所需的不同樣本集，以實(shí)現(xiàn)基分類器個(gè)體間的差異性，從而達(dá)到提高分類準(zhǔn)確率的目的，其算法流程如圖1所示。首先，定義訓(xùn)練數(shù)據(jù)集D={X,Y}={(x1,y1), (x2,y2),…,(xn,yn)}，其中：xi∈RP是訓(xùn)練集的輸入；yi∈{C1, C2,…, Cm}是訓(xùn)練集的輸出；F為訓(xùn)練集D的特征集。其次，旋轉(zhuǎn)森林算法需要預(yù)先給定參數(shù)：子集特征數(shù)s和集成規(guī)模L，它們分別代表每一個(gè)特征子集所包含的特征指標(biāo)個(gè)數(shù)和集成系統(tǒng)所包含的基分類器個(gè)數(shù)。整個(gè)算法共分成基分類器生成與基分類器合成兩個(gè)階段。

圖1 旋轉(zhuǎn)森林算法流程圖Fig.1 Procedure of rotation forest

(1)基分類器生成階段

①將原始特征集F進(jìn)行隨機(jī)分割，每個(gè)子集包含s個(gè)特征指標(biāo)，任意兩個(gè)子集間互不相交，無(wú)法整除者，將余數(shù)部分歸為一個(gè)子集，共得到k個(gè)子集，第i個(gè)基分類器的第j個(gè)特征子集記為Fij。

②從數(shù)據(jù)集X提取Fij全部樣本，采用bootstrap方法隨機(jī)抽取75%樣本生成樣本子集Xij。

③按照特征變換算法計(jì)算Xij的變換矩陣Wij。

④將變換矩陣Wi1,Wi2,…,Wik按照式(1)排列構(gòu)造稀疏矩陣Wi。

(1)

(2)基分類器合成階段

2）鉆孔沖洗。灌漿前，要進(jìn)行鉆孔孔壁沖洗和裂隙沖洗，沖洗過(guò)程中同步進(jìn)行抬動(dòng)觀測(cè)。鉆孔沖洗后，孔內(nèi)殘存的沉積物厚度不得超過(guò)20cm，需進(jìn)行鉆孔電視（鉆孔全景成像）的鉆孔需達(dá)到孔內(nèi)水清凈。

⑧利用基分類器Gi對(duì)x′i進(jìn)行預(yù)測(cè)，得到分屬C1, C2,…, Cm的概率P1(x′i),P2(x′i),…,Pm(x′i)。

⑨重復(fù)步驟7～8，得到全部基分類器的分類結(jié)果，按照式(2)對(duì)分類結(jié)果進(jìn)行集成。

(2)

⑩將樣本x劃分到概率最大的類別中，得到最終的集成分類結(jié)果。

1.2 特征變換算法

決策樹(shù)是敏感分類器，數(shù)據(jù)的任何微小改變都可能使訓(xùn)練得到的決策樹(shù)構(gòu)造完全不同[14]。因此，特征變換方法對(duì)于構(gòu)造有差異的數(shù)據(jù)分量起著至關(guān)重要的作用，并直接影響最終的分類結(jié)果。常用的特征變換方法包括主成分分析(PCA)、非參數(shù)線性判別(NDA)、稀疏隨機(jī)映射(SR)和隨機(jī)映射(R)等，大量數(shù)據(jù)實(shí)驗(yàn)結(jié)果表明：基于PCA特征變換的旋轉(zhuǎn)森林算法(PCA-RoF)性能最好[17]。

獨(dú)立成分分析(ICA)是基于高階統(tǒng)計(jì)量的多特征數(shù)據(jù)處理方法，它將數(shù)據(jù)集特征轉(zhuǎn)換為若干個(gè)統(tǒng)計(jì)獨(dú)立的特征組合，從而保證數(shù)據(jù)集的高階統(tǒng)計(jì)特性。ICA被視為PCA的一種重要擴(kuò)展，同樣適合作為旋轉(zhuǎn)森林的特征變化方法，且與其它特征變換算法相比，ICA具有以下優(yōu)勢(shì)：提供更符合實(shí)際的統(tǒng)計(jì)模型，可以更好地確定數(shù)據(jù)在高維空間的位置； 提供了一種非正交基空間，當(dāng)原始數(shù)據(jù)集存在離群點(diǎn)時(shí)，能夠更好地對(duì)數(shù)據(jù)信息進(jìn)行重建；對(duì)高階統(tǒng)計(jì)信息比較敏感，而PCA只對(duì)二階信息敏感。目前，ICA的估計(jì)算法主要有4種：FastICA、InfoMax、JADE、Radical-ICA[18]。

2 工程實(shí)例

2.1 確定特征指標(biāo)

關(guān)于礦柱的穩(wěn)定性問(wèn)題，國(guó)內(nèi)外學(xué)者有一個(gè)普遍共識(shí)，即礦柱對(duì)采礦所引起荷載的整體響應(yīng)取決于礦柱的形狀特征、礦柱巖體的地質(zhì)構(gòu)造和圍巖對(duì)礦柱所施加的表面約束特性。當(dāng)圍巖所施加的約束大于礦柱該形狀特征下巖體構(gòu)造所能承受的極限時(shí)，礦柱就會(huì)發(fā)生破壞[19]。因此，礦柱穩(wěn)定性研究可總結(jié)為如圖2所示的評(píng)價(jià)體系：礦柱的穩(wěn)定性分析主要分為宏觀的強(qiáng)度評(píng)價(jià)和微觀的應(yīng)力評(píng)價(jià)兩種途徑，理論上來(lái)說(shuō)宏觀強(qiáng)度與微觀應(yīng)力存在著關(guān)于巖質(zhì)構(gòu)造、形狀特征的一一對(duì)應(yīng)關(guān)系，但由于巖體非均質(zhì)且構(gòu)造受多種因素影響，并無(wú)有效的度量方法。因此，在尋找其對(duì)應(yīng)關(guān)系的問(wèn)題上存在無(wú)法攻克的難點(diǎn)，從而催生出一系列逼近關(guān)系的研究方法。表1整理了近年部分文獻(xiàn)的指標(biāo)選取情況，所有研究方法的考量因素均包含形狀特征量、力學(xué)狀態(tài)量和力學(xué)極限量三類，三者缺一不可。

(3)

圖2 礦柱穩(wěn)定性評(píng)價(jià)體系Fig.2 Framework of pillar stability evaluation

表1 相關(guān)文獻(xiàn)的指標(biāo)選取情況

(4)

PS=0.44·UCS·(0.68+0.52κ)

(5)

其中，κ為礦柱摩擦系數(shù)，計(jì)算公式為：

(6)

為更全面反應(yīng)礦柱穩(wěn)定情況、提高識(shí)別準(zhǔn)確率，本文除選取4個(gè)直接指標(biāo)外，同時(shí)將礦柱高寬比K、礦柱約束Cpav、礦柱摩擦系數(shù)κ以及礦柱強(qiáng)度PS納入數(shù)據(jù)庫(kù)的指標(biāo)考量范圍。盡管衍生指標(biāo)不是直接來(lái)源于實(shí)際工程，礦柱應(yīng)力δP也非基于同一方法得出，但更多特征輸入能夠給予模型必要的信息余度和特征選擇的可能性。最重要的是上述指標(biāo)根據(jù)統(tǒng)計(jì)經(jīng)驗(yàn)得出，是Per John Lunder[7]大量工作的經(jīng)驗(yàn)成果，經(jīng)過(guò)許多研究方法所檢驗(yàn)，其結(jié)果能夠較為準(zhǔn)確地反映礦柱的實(shí)際狀況(圖3)。

2.2 數(shù)據(jù)采集與處理

模型所需基礎(chǔ)數(shù)據(jù)取自文獻(xiàn)[7](共含有178例樣本)，剔除掉16例直接指標(biāo)缺失樣本，得到一個(gè)含有162個(gè)完整硬巖礦柱樣本的數(shù)據(jù)庫(kù)。關(guān)于礦柱穩(wěn)定性的分級(jí)標(biāo)準(zhǔn)，各組織和學(xué)者觀點(diǎn)不盡相同，其范圍從簡(jiǎn)單的“穩(wěn)定/破壞”兩級(jí)標(biāo)準(zhǔn)到五級(jí)或六級(jí)的更精細(xì)化標(biāo)準(zhǔn)。為實(shí)現(xiàn)分類標(biāo)準(zhǔn)的簡(jiǎn)單實(shí)用功能和模型對(duì)新樣本數(shù)據(jù)的兼容性，本文選取最簡(jiǎn)單的“穩(wěn)定/破壞”兩級(jí)分類作為模型的分級(jí)標(biāo)準(zhǔn)，同時(shí)考慮到不穩(wěn)定礦柱的破壞漸進(jìn)特性，將文獻(xiàn)[7]中通用分級(jí)方法的不穩(wěn)定礦柱和破壞礦柱合歸為同一類，穩(wěn)定礦柱依然獨(dú)成一類，分別采用“F”和“S”標(biāo)記礦柱的破壞和穩(wěn)定狀態(tài)。圖3為所建數(shù)據(jù)庫(kù)8類特征指標(biāo)的箱線圖。

2.3 模型檢驗(yàn)及應(yīng)用實(shí)例分析

模型輸入量為(W,H,δP,UCS,K,Cpav,κ,PS)，輸出為“穩(wěn)定/破壞”兩級(jí)類別V：“F”和“S”，輸入輸出間存在ICA-RoF算法的映射關(guān)系F(W,H,δP,UCS,K,Cpav,κ,PS)→V，具體檢驗(yàn)過(guò)程如圖4所示。

圖3 礦柱樣本特征箱線圖Fig.3 Boxplots of pillar sample data

圖4 礦柱穩(wěn)定性判別的ICA-RoF模型及其檢驗(yàn)過(guò)程Fig.4 The ICA-RoF model of pillar stability recognition and its testing procedure

模型的構(gòu)建與檢驗(yàn)均在Python2.7環(huán)境下進(jìn)行，相關(guān)程序通過(guò)調(diào)用NumPy、SciPy和sklearn等庫(kù)進(jìn)行編寫(xiě)，且在同一平臺(tái)上運(yùn)行，具體的電腦配置為：操作系統(tǒng)：Windows 7旗艦版64位；CPU：Intel(R) Core(TM) i3-4160 3.60 GHz；RAM：4.00 GB?；诸惼骱吞卣髯儞Q方法均直接引用于sklearn庫(kù)。其中，基分類器采用CART決策樹(shù)(CDT)，按信息增益劃分節(jié)點(diǎn)數(shù)據(jù)集，為防止模型過(guò)擬合，將節(jié)點(diǎn)最少樣本數(shù)量限制為3%(min_samples_split=0.03)，其余采用默認(rèn)設(shè)置；ICA特征變換矩陣基于FastICA函數(shù)獲得，其參數(shù)設(shè)置為：algorithm='parallel'，whiten=True，max_iter=200，tol=0.000 1。

為獲得可靠的ICA-RoF算法映射F→V，用5折交叉驗(yàn)證方法對(duì)訓(xùn)練集數(shù)據(jù)進(jìn)行遍歷尋優(yōu)計(jì)算，并根據(jù)40次運(yùn)行結(jié)果的平均識(shí)別準(zhǔn)確率選取ICA-RoF算法的最優(yōu)參數(shù)，運(yùn)行結(jié)果如圖5所示。尋優(yōu)過(guò)程同時(shí)對(duì)比PCA-RoF算法，可以明顯看到：隨著集成規(guī)模的增加，兩種算法的平均準(zhǔn)確率不斷提高，且在集成規(guī)模達(dá)到10后逐漸趨于穩(wěn)定；當(dāng)子集特征數(shù)為4時(shí)，兩種算法均獲得最高的識(shí)別準(zhǔn)確率。兩種基于不同特征變換矩陣的旋轉(zhuǎn)森林算法的規(guī)律不盡相同，但I(xiàn)CA-RoF算法的平均準(zhǔn)確率高于PCA-RoF算法，且這一結(jié)果不因集成規(guī)模和子集特征數(shù)大小而改變?？紤]到計(jì)算時(shí)間與集成規(guī)模成正比，故選取L=10、s=4作為模型性能和計(jì)算消耗的折中優(yōu)選參數(shù)。

圖5 參數(shù)尋優(yōu)運(yùn)行結(jié)果Fig.5 Results of searching parameter

當(dāng)L=10、s=4時(shí)，訓(xùn)練集數(shù)據(jù)200(50×4)次運(yùn)行結(jié)果的準(zhǔn)確率統(tǒng)計(jì)特性如表2所示。為客觀反映算法性能、科學(xué)檢驗(yàn)?zāi)Ｐ托Ч疚耐瑫r(shí)選擇PCA-RoF算法、CDT(基分類器)算法和文獻(xiàn)[23]中所述性能最為優(yōu)越的GPC算法進(jìn)行了相同條件下的對(duì)比研究，結(jié)果見(jiàn)表2。對(duì)各算法準(zhǔn)確率總體的均值之差進(jìn)行統(tǒng)計(jì)學(xué)檢驗(yàn)(獨(dú)立大樣本的Z檢驗(yàn))，經(jīng)檢驗(yàn)可認(rèn)為：在顯著性水平α=0.05的情況下，ICA-RoF算法的平均準(zhǔn)確率明顯高于其余3種算法。

同時(shí)，根據(jù)上文所構(gòu)建的4種模型對(duì)12組測(cè)試集數(shù)據(jù)進(jìn)行判別，判別結(jié)果見(jiàn)表2。從表中的識(shí)別結(jié)果可知：除了CDT算法的識(shí)別準(zhǔn)確率僅為83.33%(10/12)外，ICA-RoF算法和其余兩種算法識(shí)別結(jié)果一致，僅有1個(gè)誤分類，準(zhǔn)確率達(dá)91.66%，與實(shí)際吻合情況較好，從而說(shuō)明ICA-RoF硬巖礦柱破壞識(shí)別模型同樣具有較強(qiáng)的泛化能力。

與其它模型相比，基于ICA-RoF算法的礦柱穩(wěn)定性判別模型具有如下優(yōu)點(diǎn)：①數(shù)據(jù)處理簡(jiǎn)單，對(duì)新數(shù)據(jù)的兼容性好，有利于數(shù)據(jù)庫(kù)的拓展；② 識(shí)別準(zhǔn)確率高，泛化能力好，算法映射更為可靠；③決策樹(shù)作為基分類器，可以處理連續(xù)變量和類型變量，受異常值影響小。同時(shí)，該模型也存在著不可避免的缺點(diǎn)：①算法復(fù)雜度大，將耗費(fèi)更多的計(jì)算資源與時(shí)間；②對(duì)于不平衡數(shù)據(jù)集，ICA-RoF算法的識(shí)別結(jié)果將更傾向于數(shù)量更多的類別?？傮w來(lái)說(shuō)，使用ICA-RoF算法來(lái)識(shí)別硬巖礦柱的破壞情況，是完全可行且可靠的，無(wú)論對(duì)于前期的采礦工程設(shè)計(jì)，還是后期的工程災(zāi)害防治，都有其實(shí)際意義。更廣泛的收集礦柱穩(wěn)定性數(shù)據(jù)，建立更詳細(xì)的數(shù)據(jù)庫(kù)，開(kāi)發(fā)大數(shù)據(jù)平臺(tái)，將使該模型的可靠度更高。

表2 不同算法的性能比較

3 結(jié)論

(1)通過(guò)分析礦柱破壞機(jī)理，總結(jié)歸納礦柱穩(wěn)定性研究的一般方法體系，明確礦柱穩(wěn)定性判別模型應(yīng)包含形狀特征量、力學(xué)狀態(tài)量和力學(xué)極限量3方面指標(biāo)。三者相互聯(lián)系，缺一不可。

(2)礦柱穩(wěn)定性級(jí)別采用最簡(jiǎn)單的“穩(wěn)定/破壞”兩級(jí)分類標(biāo)準(zhǔn)。該分類標(biāo)準(zhǔn)簡(jiǎn)單實(shí)用，對(duì)新樣本數(shù)據(jù)的兼容性好，有利于對(duì)數(shù)據(jù)庫(kù)進(jìn)行拓展。

(3)通過(guò)對(duì)所建模型進(jìn)行網(wǎng)格式參數(shù)尋優(yōu)計(jì)算，發(fā)現(xiàn)：礦柱識(shí)別準(zhǔn)確率會(huì)隨著ICA-RoF的集成規(guī)模的增加而增大，當(dāng)集成規(guī)模達(dá)到10后逐漸趨于穩(wěn)定；當(dāng)子集特征數(shù)為4時(shí)，模型整體識(shí)別率最高。

(4)通過(guò)模型對(duì)比實(shí)驗(yàn)可知：ICA-RoF識(shí)別模型的識(shí)別精度高、泛化能力好，其平均準(zhǔn)確率顯著高于PCA-RoF算法、CDT算法和GPC算法(α=0.05)。