基于組合賦權和未確知測度理論的采空區(qū)安全性評價

陳志強，王紅心

(北京金誠信礦山技術研究院有限公司)

引 言

目前，中國多數(shù)礦山逐步由淺部開采向深部及超深部開采過渡，礦石挖掘后未經(jīng)充填處理會形成大面積采空區(qū)，這些采空區(qū)在地質(zhì)因素或人為擾動因素作用下，極易發(fā)生失穩(wěn)破壞，給礦山乃至周邊環(huán)境造成極大損害，如沖擊波、山體滑坡、地表沉陷、人員傷亡等[1]，因此定期對采空區(qū)的狀態(tài)進行失穩(wěn)預測就有比較重要的現(xiàn)實意義，可為后期采空區(qū)治理提供理論參考。

長期以來，正是由于無數(shù)專家、學者的刻苦鉆研探索，才使得采空區(qū)穩(wěn)定性分析方法取得長足的進步。目前，對采空區(qū)的失穩(wěn)預測主要有數(shù)值模擬、力學分析、不確定性分析等方法[2]。例如：MORTAZAVI等[3]利用離散元軟件UDEC對采空區(qū)中礦柱的變形及破壞過程進行模擬分析，相關計算結果與現(xiàn)場實測較吻合。何標慶[4]首先對采空區(qū)進行三維掃描測量，然后結合有限元Flac3D軟件對其進行模擬分析，結果顯示采空區(qū)上下盤易發(fā)生拉破壞，存在大規(guī)模冒落危險，因而對采空區(qū)及時進行處理，避免礦山發(fā)生重大災害事故。李玉飛等[5]采用力學分析方法，將采空區(qū)頂板概化為固支梁力學結構模型，基于突變理論獲得失穩(wěn)判據(jù)，之后進行實例驗證，驗證結果與現(xiàn)場情況較相符。LI等[6]在現(xiàn)有力學模型基礎上，利用黏彈性模型，進行分析計算，成功預測了煤柱的變形。陳嬌[7]運用層次分析法對某鉛鋅礦采空區(qū)穩(wěn)定性影響因子計算權重，然后運用系統(tǒng)工程理論進行分析計算，并在此基礎上提出了針對性建議。

以上方法雖然在實際中取得了較好的效果，但是也存在著一定的局限性。例如：數(shù)值模擬及力學分析在計算時，需要理想化很多條件；不確定性方法中，權重計算容易受主觀因素或客觀實測數(shù)據(jù)偏差影響等。鑒于采空區(qū)是一個綜合的復雜系統(tǒng)，其失穩(wěn)狀態(tài)受多種定性、定量因素影響，如地質(zhì)構造、采空區(qū)埋藏深度、礦區(qū)含水條件、采空區(qū)跨高比等，所以對采空區(qū)的失穩(wěn)預測需要采用主客觀相結合的方法，提出基于AHP-信息熵組合賦權的未確知測度理論對采空區(qū)失穩(wěn)狀態(tài)進行預測。

1 采空區(qū)失穩(wěn)預測模型構建

采空區(qū)失穩(wěn)預測是分析和治理采空區(qū)過程中必不可少的重要環(huán)節(jié)，分析時將采空區(qū)失穩(wěn)情形劃分為極穩(wěn)定、穩(wěn)定、不穩(wěn)定、極不穩(wěn)定4種狀態(tài)，采用主客觀賦權相結合的AHP-信息熵組合賦權的方法確定影響因素權重，以及運用專家打分法及線性公式法作為評價因子隸屬度函數(shù)，構建采空區(qū)未確知測度模型，進行分析預測。

1.1 未確知測度原理

王光遠教授于20世紀90年代提出未確知信息理論，其后被廣大科研學者成功應用于多個領域。例如：秦新生[8]運用未確知測度模型對10家跨境電商物流企業(yè)服務能力進行綜合分析，幫助電商企業(yè)優(yōu)化選擇物流服務，從而降低成本。楊清平等[9]運用未確知測度理論分析保安礦柱的穩(wěn)定性，結果顯示，這種評價方法可行且可為礦山安全生產(chǎn)提供借鑒。黃丹等[10]運用未確知測度-集對分析法對巖質(zhì)邊坡穩(wěn)定性進行評價，其分析不僅滿足工程要求，且與實際相符。

未確知測度理論基本思路是，若某預測對象含有m個因素層，其中包括n個評價指標，用X表示，則，X={x1，x2，…，xn}，且任意xi有p個評價等級C1，C2，…，Ck，…，Cp，記U={C1，C2，…，Ck，…，Cp}，且滿足第k個評價等級大于第k+1個評價等級，即Ck>Ck+1。

1)單指標測度。若μijk=μ(xij∈Ck)表示實測值xij屬于第k個評價等級Ck的程度，則μ應滿足：

0≤μ(xij∈Ck)≤1

(1)

μ(xij∈U)=1

(2)

(3)

式中：μ為未確知測度，簡稱測度；μijk為單指標測度評價矩陣。

2)多指標綜合測度評價向量。若μijk=μ(Ri∈Ck)為評價樣本Ri屬于第k個評價類Ck的程度，則有：

(4)

3)置信度識別準則。最終預測結果，選用置信度識別準則[11]：設為置信度(λ≥0.5，通常取λ=0.6或0.7，這里取λ=0.6)，若C1>C2>…>Ck>…>Cp，且有:

(5)

則認為評價方案Ri屬于第k0評價等級。

1.2 權重計算

1)AHP。AHP法是兩兩指標相互比較，并用1～9整數(shù)值標度其相對大小，形成判斷矩陣，具體比例標度見表1。

得到判斷矩陣后，采用最大特征值法求得特征向量，再將特征向量進行絕對值化和歸一化處理，最后進行一致性檢驗得到因子權重向量(φ1,φ2,…,φj)，記為α。

表1 AHP法比例標度

2)信息熵[12]。

(1)數(shù)據(jù)標準化。將實測數(shù)據(jù)組成原始信息矩陣，按式(6)進行標準化：

(6)

(2)計算評價因子的熵：

(7)

(8)

式中：L為評價對象個數(shù)；j為第j個評價因子。若pij=0，lnpij無意義，所以這里規(guī)定lnpij≠0。

(3)確定權重：

(9)

最后，記因子權重向量β=(w1,w2,…,wj)。

3)組合賦權法[13]。

最終權重向量采用乘法集成法將向量α和β融合，計算公式如下：

(10)

1.3 預測模型指標選取

鑒于采空區(qū)失穩(wěn)狀態(tài)影響因素具有模糊、多變等復雜特點，這里著眼于容易收集且易于量化的原則，選取工程地質(zhì)、采空區(qū)形態(tài)、其他因素三大類共12項因子來建立預測模型指標體系，具體因素分配見圖1。

圖1 預測模型因素分配指標體系

1.4 采空區(qū)失穩(wěn)影響因素分級標準

由于礦山所處的環(huán)境及選用開采工藝的不同，形成的采空區(qū)賦存狀態(tài)千差萬別，影響因素分級標準也不同。根據(jù)礦山具體實際，結合國內(nèi)外相似礦山及參考大量相關文獻、礦山規(guī)范等，采用表2所示的參數(shù)作為各因子分級標準。

表2 預測模型因子分級標準

1.5 單指標測度函數(shù)

綜合各種采空區(qū)影響因子分布特征，選用廣泛應用的專家打分法及線性公式法作為其單指標測度函數(shù)，見表3、圖2。

表3 離散型單指標測度函數(shù)

圖2 連續(xù)型指標測度函數(shù)

2 工程應用

內(nèi)蒙古某金礦，已知控制礦體長2 500 m，控制斜深1 100 m。礦體分布于140勘探線—263勘探線，礦體呈層狀、似層狀產(chǎn)出。礦體傾角45°～85°，平均55°，深部有變緩趨勢，礦體厚度0.43～8.36 m。礦山設計采用淺孔留礦采礦法開采，且挖掘形成的采空區(qū)大部分未進行處理，目前開采標高818 m中段以上均已采掘完畢，采礦最低中段為118 m中段，開采深度950 m，礦山面臨嚴峻的安全隱患。現(xiàn)對礦區(qū)部分采空區(qū)失穩(wěn)狀態(tài)進行分析預測，具體調(diào)查數(shù)據(jù)見表4。

表4 采空區(qū)各指標調(diào)查數(shù)據(jù)

下面以403采空區(qū)為例，進行計算分析。首先由圖2、表3可得403采空區(qū)的測度矩陣為：

其次，根據(jù)AHP法計算原理，得出403采空區(qū)權重向量α=(0.108，0.108，0.044，0.061，0.108，0.060，0.060，0.280，0.108，0.013，0.018，0.035)。

再根據(jù)403采空區(qū)各指標測度矩陣及式(7)～(9)，得出403采空區(qū)權重向量β=(0.090，0.155，0.017，0.017，0.087，0.155，0.102，0.099，0.054，0.155，0.017，0.054)。

最后由組合賦權公式(10)得出因素最終權重向量ψ=(0.107，0.184，0.008,0.011，0.104，0.102，0.067，0.305，0.065，0.022，0.003，0.021)。

由式(4)得出403采空區(qū)綜合評價向量μ11=ψ×(μ1jk)12×4=(0.394，0.391，0.097，0.118)。根據(jù)置信度識別準則原理，10=0.394<0.6，而20=0.394+0.391=0.785>0.6，故403采空區(qū)屬于第20級，即Ⅱ級(穩(wěn)定)。

同理，利用以上計算方法可得其他采空區(qū)失穩(wěn)狀態(tài)預測結果，見表5。

表5 采空區(qū)失穩(wěn)狀態(tài)未確知測度預測結果

從以上計算結果可以看出，各采空區(qū)都處于Ⅱ級(穩(wěn)定)狀態(tài)，與目前現(xiàn)場采空區(qū)表現(xiàn)的實際狀態(tài)相符，表明此采空區(qū)失穩(wěn)狀態(tài)預測分析方法具備一定的可行性和有效性。同時，由于礦山在不斷生產(chǎn)，采空區(qū)一直受爆破等因素擾動，故不能掉以輕心，應加強監(jiān)測管理，防止發(fā)生大規(guī)模災害事故。

3 結 論

1)從影響采空區(qū)失穩(wěn)的諸多復雜因素中，選取單軸抗壓強度、完整性系數(shù)、采空區(qū)頂板暴露面積、埋藏深度、跨高比、地質(zhì)構造等12項因素，作為預測采空區(qū)失穩(wěn)狀態(tài)的指標，構建未確知測度模型。

2)針對模型中的定性指標和定量指標，運用主客觀相結合的方法，即AHP-信息熵組合賦權的指標權重計算方法，使得權重計算結果更專業(yè)、更科學、更合理，進而得到的未確知預測結果更具有信服力。

3)將此分析方法運用于內(nèi)蒙古某金礦，其采空區(qū)狀態(tài)預測結果與現(xiàn)場實際情況較為相符，表明此方法的可行性和有效性，豐富了采空區(qū)穩(wěn)定性分析方法。