復(fù)雜采空區(qū)處理方案的模糊推理系統(tǒng)

張東煒， 熊雪強(qiáng)， 趙 奎， 鐘林傳

(1.江西理工大學(xué)資源與環(huán)境工程學(xué)院，江西 贛州 341000；2.武平紫金礦業(yè)有限公司，福建 上杭 364200)

復(fù)雜采空區(qū)處理方案的模糊推理系統(tǒng)

張東煒1， 熊雪強(qiáng)1， 趙 奎1， 鐘林傳2

(1.江西理工大學(xué)資源與環(huán)境工程學(xué)院，江西 贛州 341000；2.武平紫金礦業(yè)有限公司，福建 上杭 364200)

對(duì)于復(fù)雜采空區(qū)的現(xiàn)狀，采空區(qū)處理是亟待解決的問(wèn)題.在選擇充填和封閉聯(lián)合處理采空區(qū)的前提下，通過(guò)使用關(guān)系矩陣與模糊數(shù)學(xué)聯(lián)合的采空區(qū)穩(wěn)定性分析方法，得到了各采空區(qū)的穩(wěn)定性系數(shù)，結(jié)合采空區(qū)穩(wěn)定性系數(shù)及地壓監(jiān)測(cè)結(jié)果，進(jìn)行了采空區(qū)處理方案的模糊推理系統(tǒng)研究，通過(guò)研究能得到與每個(gè)采空區(qū)相適應(yīng)的處理方法，為礦山的地壓處理提供了相應(yīng)的理論依據(jù)，并為礦山的安全生產(chǎn)提供了保障，同時(shí)也對(duì)類似采空區(qū)處理提供了一定的依據(jù).

采空區(qū)處理；模糊推理；復(fù)雜采空區(qū)；采空區(qū)穩(wěn)定性

0 引 言

空?qǐng)龇ㄊ悄壳笆褂幂^廣的一種開(kāi)采方法，但是該方法遺留下許多未處理的采空區(qū)，復(fù)雜的采空區(qū)可能會(huì)引發(fā)大規(guī)模的沖擊性地壓，造成不必要的財(cái)產(chǎn)及人員傷害，目前對(duì)采空區(qū)處理的方法大致分為崩落圍巖處理采空區(qū)、充填料充填處理空區(qū)、留永久礦柱或構(gòu)筑人工石柱處理空區(qū)、封閉和隔離處理采空區(qū)、聯(lián)合法處理采空區(qū)[1].目前隨著對(duì)環(huán)境保護(hù)要求的提高，使用充填法處理采空區(qū)越來(lái)越廣泛.充填體主要有3種作用原理：應(yīng)力吸收與轉(zhuǎn)移，接觸支撐和應(yīng)力隔離[2].文中主要針對(duì)采用充填和封閉聯(lián)合處理采空區(qū)時(shí)的具體區(qū)域進(jìn)行研究.

目前對(duì)于采空區(qū)穩(wěn)定性的分析方法很多，主要有理論分析、模糊分析、數(shù)值模擬分析.其中模糊理論分析運(yùn)用較廣.文中采用的基于關(guān)系矩陣與模糊數(shù)學(xué)的綜合評(píng)價(jià)方法就是建立在模糊理論之上的一種有效的采空區(qū)穩(wěn)定性分析方法，該方法考慮了6種采空區(qū)的主要影響因素，結(jié)合各影響因素的活動(dòng)性指數(shù)與其對(duì)應(yīng)的隸屬度值算出各采空區(qū)穩(wěn)定性系數(shù)，以此系數(shù)進(jìn)行采空區(qū)穩(wěn)定性的初步評(píng)價(jià)，但是，現(xiàn)場(chǎng)的情況往往是非常復(fù)雜的，且不斷變化著.此方法中各主要因素隸屬度的選取也有很大的模糊性，且不能反應(yīng)采空區(qū)的動(dòng)態(tài)變化，變化著的采場(chǎng)其穩(wěn)定性只有通過(guò)地壓監(jiān)測(cè)才能有所反應(yīng)，地壓監(jiān)測(cè)也具有一定的模糊性，為了能得到更加符合實(shí)際的結(jié)果，必須結(jié)合靜態(tài)的采空區(qū)穩(wěn)定性系數(shù)與動(dòng)態(tài)的地壓監(jiān)測(cè)的模糊結(jié)果，而模糊推理系統(tǒng)便能有效的整合各種模糊輸入變量，得到與實(shí)際相符的精確結(jié)果[3].下面詳細(xì)介紹該模糊推理系統(tǒng)的建立.

1 模糊推理系統(tǒng)的建立

模糊推理系統(tǒng)主要由輸入輸出變量的模糊化、模糊規(guī)則庫(kù)建立、模糊推理方法選擇及去模糊化幾部分組成.

1.1 關(guān)系矩陣及模糊數(shù)學(xué)聯(lián)合的采空區(qū)穩(wěn)定性分析

1.1.1 采空區(qū)穩(wěn)定性的關(guān)系矩陣

由于影響采空區(qū)穩(wěn)定性方面的因素很多，文中采取文獻(xiàn)[4]和[5]所提出的關(guān)系矩陣及模糊數(shù)學(xué)的聯(lián)合采空區(qū)穩(wěn)定性評(píng)價(jià)方法，其選擇的影響采空區(qū)穩(wěn)定性的主要因素有：①巖石單軸抗壓強(qiáng)度R；②頂板暴露面積A；③原巖應(yīng)力狀況；④結(jié)構(gòu)面狀況；⑤地下水W；⑥礦柱狀況P.其中結(jié)構(gòu)面狀況又可以分為結(jié)構(gòu)面間距D和結(jié)構(gòu)面產(chǎn)狀與頂板的方位關(guān)系O2個(gè)因素，原巖應(yīng)力狀況分為最大初始應(yīng)力和圍巖單軸抗壓強(qiáng)度比值ST及最大主應(yīng)力方向與采場(chǎng)頂板方位之間的關(guān)系SL 2個(gè)因素，各因素的相對(duì)重要程度使用活動(dòng)性系數(shù)來(lái)衡量，系數(shù)越大則說(shuō)明其對(duì)整個(gè)評(píng)價(jià)系統(tǒng)的貢獻(xiàn)越大.主要因素aii的活動(dòng)性系數(shù)由式（1）確定：

式中FA（i）對(duì)應(yīng)的活動(dòng)性系數(shù)，n為主要因素個(gè)數(shù)；aij為因素aii對(duì)因素aij的影響程度.

總結(jié)上述各主要因素之間的相互作用，依據(jù)其相互的影響強(qiáng)度，結(jié)合實(shí)踐經(jīng)驗(yàn)，給關(guān)系矩陣賦值，從而得出與其對(duì)應(yīng)的關(guān)系矩陣，如式（2）所示：

從而由活動(dòng)性計(jì)算式便可得到各活動(dòng)性系數(shù)為：FA（1）=0.118， FA（2）=0.191， FA（3）=0.127，F(xiàn)A（4）=0.064， FA（5）=0.127， FA（6）=0.082， FA（7）=0.127，F(xiàn)A（8）=0.164.

1.1.2 采空區(qū)穩(wěn)定性系數(shù)

由上述建立的關(guān)系矩陣，根據(jù)模糊數(shù)學(xué)理論選取與各影響因素合適的隸屬度函數(shù)，其每個(gè)影響因素的隸屬度函數(shù)的構(gòu)造由于篇幅原因暫不進(jìn)行詳細(xì)概述，這里只介紹巖石單軸抗壓強(qiáng)度與最大初始應(yīng)力比值ST的隸屬函數(shù)構(gòu)造，實(shí)測(cè)資料顯示，當(dāng)ST為3～6時(shí)就會(huì)發(fā)生巖爆和巖心餅化，當(dāng)在5～10和10～200時(shí)分別對(duì)應(yīng)高、低應(yīng)力，故由此選擇隸屬函數(shù)的特征點(diǎn)為3、6、10、200；各主要因素的隸屬度函數(shù)曲線如下圖1所示.

通過(guò)具體采空區(qū)的現(xiàn)狀結(jié)合隸屬函數(shù)選擇對(duì)應(yīng)的隸屬度值，最終建立起采空區(qū)穩(wěn)定性的評(píng)價(jià)指標(biāo)，即采空區(qū)穩(wěn)定性系數(shù)：各主要因素活動(dòng)性評(píng)價(jià)指標(biāo)FA（i）與其對(duì)應(yīng)的隸屬度 W（i）進(jìn)行乘積，并將所有因素的乘積結(jié)果進(jìn)行累加，其累加值就定義為采空區(qū)穩(wěn)定性系數(shù)[6]，記為 MS，即：

式（3）中：W（i）為主要評(píng)價(jià)因素的隸屬度；FA（i）為前述的各影響因素活動(dòng)性指數(shù).

1.2 輸入輸出變量的隸屬度函數(shù)確立

選擇空區(qū)穩(wěn)定性系數(shù)及地壓監(jiān)測(cè)結(jié)果作為輸入變量，采空區(qū)處理方案作為輸出變量，從靜態(tài)和動(dòng)態(tài)兩個(gè)方面進(jìn)行采空區(qū)處理研究，其結(jié)果符合實(shí)際狀況，其中輸入、輸出變量可采用高斯型隸屬度函數(shù)如式（4）所示：輸入、輸出隸屬度函數(shù)具體建立過(guò)程如下.

根據(jù)穩(wěn)定性系數(shù)計(jì)算結(jié)果，將穩(wěn)定性系數(shù)劃分為3個(gè)變量：不穩(wěn)定、中等穩(wěn)定、穩(wěn)定.結(jié)合許多其他開(kāi)采礦山的采空區(qū)賦存狀況進(jìn)行空區(qū)穩(wěn)定性的計(jì)算分析，發(fā)現(xiàn)采空區(qū)穩(wěn)定性系數(shù)一般在0.5附近，故選擇采空區(qū)系數(shù)在0.5時(shí)對(duì)應(yīng)于中等穩(wěn)定的隸屬度為1，穩(wěn)定與不穩(wěn)定隸屬度函數(shù)曲線則隨著穩(wěn)定性系數(shù)的增大而分別增大和減小，函數(shù)曲線如圖2所示.

現(xiàn)場(chǎng)監(jiān)測(cè)使用最廣泛的監(jiān)測(cè)儀器便為光彈應(yīng)力計(jì)，因?yàn)槠涞蛢r(jià)實(shí)用，故文中所有監(jiān)測(cè)數(shù)據(jù)均針對(duì)光應(yīng)力計(jì)結(jié)果，地壓監(jiān)測(cè)結(jié)果主要包括應(yīng)力大小和一段時(shí)期的應(yīng)力增量大小，應(yīng)力增量有正負(fù)之分，其中正說(shuō)明處于應(yīng)力集中狀態(tài)，負(fù)則說(shuō)明處于應(yīng)力釋放狀態(tài).監(jiān)測(cè)應(yīng)力的大小與許多因素有關(guān)，如儀器布置離空區(qū)越近則其所受的力越大，且埋藏越深對(duì)應(yīng)的應(yīng)力也越大，反之亦然，所以僅僅依靠應(yīng)力大小無(wú)法確定應(yīng)力集中區(qū)域，只能作為判斷的因素之一，故同時(shí)選擇應(yīng)力增量作為影響采空區(qū)處理的另一個(gè)重要因素，應(yīng)力增量可以看作采空區(qū)在一定時(shí)期內(nèi)其內(nèi)部的動(dòng)態(tài)變化，應(yīng)力及應(yīng)力增量在排除施工因素的影響后能反應(yīng)監(jiān)測(cè)區(qū)域的應(yīng)力集中狀況.

圖1 各主要因素隸屬度函數(shù)曲線

圖2 空區(qū)穩(wěn)定性隸屬度函數(shù)曲線圖

對(duì)于監(jiān)測(cè)的應(yīng)力增量，如何判斷其危險(xiǎn)性大小，目前只有依據(jù)其與礦柱的強(qiáng)度進(jìn)行比較，而礦柱影響因素也非常復(fù)雜，對(duì)于目前許多對(duì)礦柱強(qiáng)度計(jì)算的公式，選擇最具代表性的公式，如式（5）所示：

其中：σp為礦柱的強(qiáng)度，/MPa，σp1為寬高比為1時(shí)的礦柱強(qiáng)度，/MPa，ω、h分別為礦柱的寬度和高度，/m，其他為不同巖體所對(duì)應(yīng)的常系數(shù).

在這類公式中選擇較為廣用的Bieniawski提出的強(qiáng)度公式,如式（6）所示：

由于上述未考慮節(jié)理裂隙對(duì)礦柱的影響，為此采用Hoek-Brown計(jì)算節(jié)理巖體強(qiáng)度的方法作為σp1的計(jì)算公式，即：

式（7）中 σc為完整巖石的單軸抗壓強(qiáng)度，/MPa；s表示巖體完整性的一個(gè)常數(shù).對(duì)于含有多組節(jié)理且沒(méi)有明顯的各向異性的巖體s值由式（8）確定：

其中RMR為巖體地質(zhì)分類指標(biāo).

對(duì)于礦柱的平均應(yīng)力的計(jì)算也有很多方法，但是其由于模型進(jìn)行了簡(jiǎn)化，計(jì)算結(jié)果與實(shí)際相差很大，因?yàn)樗鼈兙鶠榭紤]原巖應(yīng)力的影響，為此選擇公式（9）作為礦柱平均應(yīng)力的計(jì)算：

其中K為采場(chǎng)開(kāi)挖前的礦柱集中系數(shù)，一般取1～2；σ0為原巖應(yīng)力，/MPa，可以取自重應(yīng)力；σ1為二次回采增量即為我們監(jiān)測(cè)過(guò)程中的應(yīng)力增量，/MPa.

通過(guò)上述對(duì)礦柱強(qiáng)度的計(jì)算，用礦柱強(qiáng)度除以礦柱平均應(yīng)力便能作為衡量礦柱危險(xiǎn)性大小的一個(gè)系數(shù)，稱為安全系數(shù)[7]，具體公式如式（10）所示：

安全系數(shù)S與應(yīng)力增量σ1的關(guān)系可近似作為選取應(yīng)力增量隸屬函數(shù)的依據(jù).通過(guò)大量礦柱現(xiàn)場(chǎng)實(shí)測(cè)資料結(jié)果發(fā)現(xiàn)安全系數(shù)小于1對(duì)應(yīng)的礦柱破壞嚴(yán)重，一般大于2時(shí)礦柱都很穩(wěn)定，故選擇安全系數(shù)范圍在0～3之間，同樣把安全系數(shù)也劃分為3個(gè)變量：不穩(wěn)定、中等穩(wěn)定、穩(wěn)定.其對(duì)應(yīng)的隸屬度函數(shù)如圖3所示.

圖3 安全系數(shù)隸屬度函數(shù)曲線圖

通過(guò)對(duì)許多礦山光應(yīng)力計(jì)的實(shí)際監(jiān)測(cè)結(jié)果分析，發(fā)現(xiàn)光應(yīng)力計(jì)條紋級(jí)數(shù)大的區(qū)域應(yīng)力均比較集中，因?yàn)楦浇脖厝话殡S著一些地壓顯現(xiàn)，在對(duì)銅陵紫金的現(xiàn)場(chǎng)監(jiān)測(cè)發(fā)現(xiàn)光彈應(yīng)力計(jì)壓壞后其礦柱也隨之破壞，室內(nèi)試驗(yàn)得到光應(yīng)力計(jì)所能承受荷載約為7 MPa，故應(yīng)力的選取范圍為0～7 MPa，同時(shí)對(duì)各礦山的監(jiān)測(cè)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，一般選取0～3 MPa范圍內(nèi)所監(jiān)測(cè)區(qū)域?yàn)榉€(wěn)定區(qū)域，3～5 MPa為中等穩(wěn)定區(qū)域，5～7 MPa為不穩(wěn)定區(qū)域.故將應(yīng)力水平也劃分為3個(gè)變量：不穩(wěn)定、中等穩(wěn)定、穩(wěn)定.函數(shù)曲線如圖4所示.

圖4 應(yīng)力隸屬度函數(shù)曲線圖

輸出變量則為采空區(qū)對(duì)應(yīng)的處理方案，分別對(duì)應(yīng)充填、封閉、暫不處理.函數(shù)曲線如圖5所示.

圖5 處理方案隸屬度函數(shù)曲線圖

1.3 模糊推理系統(tǒng)規(guī)則庫(kù)建立

模糊規(guī)則庫(kù)代表著整個(gè)推理系統(tǒng)的思考法則，其產(chǎn)生方式一般有3種方法：①直接由操作員的經(jīng)驗(yàn)技巧建立模糊規(guī)則；②根據(jù)受控體對(duì)輸入與輸出的反應(yīng)去歸納受控行為，以誤差法進(jìn)行設(shè)計(jì)；③控制系統(tǒng)本身進(jìn)行學(xué)習(xí)或修正控制規(guī)則.文中采用第一種方法建立模糊規(guī)則庫(kù)，模糊規(guī)則以if-Then條列式的形式，表示結(jié)合了人類判斷的模糊性控制法則，每個(gè)控制規(guī)則由前后件組成，通式如式（11）所示：

其中Rl即表示第l條規(guī)則；x為輸入及 Bl為對(duì)應(yīng)于變量的模糊集合；y為輸出，用來(lái)控制受控系統(tǒng).依據(jù)實(shí)際情況共建立27條規(guī)則，其結(jié)果如表1所示.其蘊(yùn)含的規(guī)則思想選擇第一條規(guī)則進(jìn)行解釋即：if輸入的空區(qū)穩(wěn)定性系數(shù)對(duì)應(yīng)不穩(wěn)定，安全系數(shù)對(duì)應(yīng)穩(wěn)定且應(yīng)力對(duì)應(yīng)穩(wěn)定，then處理方案選擇封閉.表1考慮了所有可能的情況，因此能滿足實(shí)際情況.

表1 模糊規(guī)則庫(kù)

1.4 模糊推理系統(tǒng)的合成

依據(jù)實(shí)際選擇Mamdani推理算法作為模糊推理系統(tǒng)的合成算法[8]，該算法與Zadeh提出的推理合成算法本質(zhì)相同[9]，是在選擇合理的規(guī)則后，就建立了一種模糊蘊(yùn)含關(guān)系，在該蘊(yùn)含關(guān)系下選擇各因素的計(jì)算或分析結(jié)果作為輸入，對(duì)輸入結(jié)果進(jìn)行模糊化處理后得到各規(guī)則對(duì)應(yīng)的模糊集，將該模糊集作為輸入，經(jīng)合成處理后得到的一個(gè)綜合模糊集作為輸出.最后選擇合適的反模糊化方法將綜合模糊集以確定的數(shù)值形式輸出.一般采用重心法來(lái)求出精確值，文中采用中心法作為反模糊化方法.

由上述輸入、輸出變量及規(guī)則庫(kù)結(jié)果，運(yùn)用Matlab模糊邏輯工具箱建立模糊推理系統(tǒng).其對(duì)應(yīng)的不同變量之間的關(guān)系曲面圖如圖6～圖8所示.

通過(guò)對(duì)輸入輸出曲線、曲面檢驗(yàn)：輸入變量采空區(qū)穩(wěn)定性系數(shù)、安全系數(shù)、應(yīng)力具有等價(jià)性，通過(guò)對(duì)單輸入曲線得到應(yīng)力為遞減函數(shù)，空區(qū)穩(wěn)定性及安全系數(shù)均為遞增函數(shù)，說(shuō)明與實(shí)際相符合，對(duì)應(yīng)的隸屬度函數(shù)及規(guī)則的確定是合理的.

圖6 安全系數(shù)與空區(qū)穩(wěn)定性系數(shù)關(guān)系曲面圖

圖7 安全系數(shù)與應(yīng)力關(guān)系曲面圖

2 應(yīng)用實(shí)例

2.1 輸入變量分析

武平紫金礦業(yè)有限公司悅洋銀多金屬礦是2007年1月由紫金礦業(yè)集團(tuán)股份有限公司出資收購(gòu)的民營(yíng)企業(yè)，礦山主要采用空?qǐng)龇恐ㄟM(jìn)行開(kāi)采，至今形成了超過(guò)約61.5萬(wàn)m3采空區(qū)，通過(guò)對(duì)礦山各中段的現(xiàn)場(chǎng)調(diào)查及點(diǎn)荷載試驗(yàn)獲得了各采空區(qū)的賦存狀況及巖石的力學(xué)參數(shù)，通過(guò)調(diào)查結(jié)果共布置77個(gè)應(yīng)力監(jiān)測(cè)儀器，其中光彈應(yīng)力計(jì)60個(gè).得到數(shù)據(jù)幾千個(gè)，采用上述采空區(qū)穩(wěn)定性分析方法對(duì)每個(gè)采空區(qū)進(jìn)行了穩(wěn)定性系數(shù)計(jì)算，并通過(guò)地壓監(jiān)測(cè)得到各采空區(qū)相應(yīng)測(cè)點(diǎn)的應(yīng)力及應(yīng)力增量值，由應(yīng)力增量算出其對(duì)應(yīng)的安全系數(shù)，由于篇幅有限故只列舉-145 m中段部分空區(qū)數(shù)據(jù)，如表2所示.

其中采空區(qū)穩(wěn)定性采用前述計(jì)算方法，-145 m中段主要回采礦體為Ⅰ號(hào)礦體.礦體受花崗巖“舌狀體”與樓子壩群變質(zhì)巖接觸帶控制，巖石普遍綠泥石化，圍巖單軸抗壓強(qiáng)度在30～95.9 MPa.該中段與表2對(duì)應(yīng)的采空區(qū)暴露面積分別為 1076.7、182.7、1259.4、140、196.2、829.0 m2；結(jié)構(gòu)面間距主要集中在10～20 cm，傾角在 30～80 °；自重應(yīng)力 rH=9.1 MPa；中段留礦柱大部分已破壞，空區(qū)內(nèi)見(jiàn)少量滴水.選用相應(yīng)隸屬度后計(jì)算空區(qū)穩(wěn)定性系數(shù)結(jié)果見(jiàn)表2.

表2 -145 m中段部分空區(qū)穩(wěn)定性系數(shù)表

通過(guò)前期調(diào)查及試驗(yàn)可知礦體賦存巖石為中細(xì)?；◢弾r，其抗壓強(qiáng)度經(jīng)點(diǎn)荷載試驗(yàn)轉(zhuǎn)換約為135 MPa，礦柱寬高比選1，RMR經(jīng)計(jì)算約為70，K取1， 自重應(yīng)力取 10.4 MPa，A、B按 Bieniawski提出的強(qiáng)度取值分別取0.64、0.36[10-11].應(yīng)力增量取前后2次監(jiān)測(cè)應(yīng)力的差值，通過(guò)計(jì)算其對(duì)應(yīng)結(jié)果如表2所示.

應(yīng)力值有光彈應(yīng)力計(jì)監(jiān)測(cè)條紋經(jīng)轉(zhuǎn)換后的應(yīng)力值[12]，同時(shí)也有鉆孔應(yīng)力計(jì)監(jiān)測(cè)的直接應(yīng)力值.其中采空區(qū)附近多組監(jiān)測(cè)結(jié)果選取其中應(yīng)力最大值作為代表.如83-2采空區(qū)選取最大6.651 MPa代表該區(qū)域應(yīng)力.

2.2 模糊推理系統(tǒng)計(jì)算結(jié)果

通過(guò)建立起來(lái)的模糊推理系統(tǒng)使用evalfis函數(shù)根據(jù)已有的各組計(jì)算結(jié)果計(jì)算出與其相對(duì)應(yīng)的采空區(qū)處理系數(shù)，-145 m中段處理系數(shù)分別為：0.3404、0.5498、0.4847、0.5250、0.2949、0.4697 通過(guò)系數(shù)確定不同的采空區(qū)的處理方法，通過(guò)對(duì)輸出變量的隸屬度函數(shù)觀察，3條隸屬函數(shù)曲線的交點(diǎn)為0.35和0.65，即這2點(diǎn)可以理解為最模糊點(diǎn)，故以此作為其處理方案的依據(jù)：系數(shù)在0～0.35之間的選用充填處理采空區(qū)；系數(shù)在0.35～0.65選擇封閉處理采空區(qū)；0.65～1選擇暫不處理采空區(qū).故-145 m中段79-1、83-2 采空區(qū)采取充填處理；83-1、85-1、83-2、1-2 采區(qū)封閉處理.

3 結(jié)束語(yǔ)

文中通過(guò)地壓監(jiān)測(cè)結(jié)果、以及采空區(qū)穩(wěn)定性系數(shù)建立了采空區(qū)處理的模糊推理系統(tǒng)，得到了復(fù)雜采空區(qū)群下各個(gè)采空區(qū)的處理系數(shù)，為每個(gè)采空區(qū)的處理提出了具體的處理方案，能起到緩解地應(yīng)力集中的狀況，保障作業(yè)安全，對(duì)類似礦山采空區(qū)處理也提供了一定的依據(jù).

由于輸入、輸出變量的隸屬度函數(shù)的選取有一定人為因素的影響，同時(shí)規(guī)則庫(kù)的建立也包含了許多模糊因素，因此其處理結(jié)果還不能排除人為因素的影響，只能在此基礎(chǔ)上逐漸逼近實(shí)際結(jié)果，故在這些方面還有待于進(jìn)一步的研究.

[1]羅毅莎，姚振鞏.采空區(qū)安全處理的對(duì)策措施研究[J].采礦技術(shù)，2010，10(4)：65-67.

[2]趙傳卿，胡乃聯(lián).膠結(jié)充填對(duì)采場(chǎng)穩(wěn)定性的影響[J].遼寧工程技術(shù)大學(xué)學(xué)報(bào)，2008，27(1)：13-16.

[3]趙 奎，萬(wàn)林海，饒運(yùn)章，等.基于聲波測(cè)試的礦柱穩(wěn)定性模糊推理系統(tǒng)及其應(yīng)用[J].巖石力學(xué)與工程學(xué)報(bào)，2004，23(11)：1804-1809.

[4]趙 奎.巖金礦山采空區(qū)及殘留礦柱回采穩(wěn)定性研究[D].贛州：南方冶金學(xué)院，2003.

[5]鐘春暉，趙 奎，郭秋根，等.基于關(guān)系矩陣和模糊理論的采空區(qū)穩(wěn)定性綜合評(píng)價(jià)[J].采礦科學(xué)技術(shù)前沿論壇論文集，2006，44-46.

[6]馬尚權(quán)，王立兵.基于關(guān)系矩陣和模糊集合的煤與瓦斯突出綜合評(píng)價(jià)[J].現(xiàn)代礦業(yè)，2001，15(10)：49-51.

[7]趙 奎.礦山巖石力學(xué)若干測(cè)試技術(shù)及其分析方法[M].北京：冶金工業(yè)出版社，2009：31-35.

[8]呂大江，石志寒，雷英杰，等.基于Mamdani蘊(yùn)含關(guān)系的直覺(jué)模糊推理仿真[J].計(jì)算機(jī)工程與設(shè)計(jì)，2010，31(11)：2539-2542.

[9]聞 新，周 露，李東江，等.MATLAB模糊邏輯工具箱的分析與應(yīng)用[M].北京：科學(xué)出版社，2011：41-45.

[10]趙 奎，胡慧明，王曉軍，等.某金礦房柱法采場(chǎng)人工礦柱參數(shù)選取[J].采礦技術(shù)，2011，11(2)：15-17.

[11]Z.T.Bieniawski.Estimating the strength of rock materials[J].J.S.Afr.Inst.Min.Metall，1974（74）：312-320.

[12]郭建平.光應(yīng)力計(jì)監(jiān)測(cè)采空區(qū)穩(wěn)定性在遂昌金礦的應(yīng)用[J].有色金屬，2004，56(3)：115-117.

Complicated goaf treatment scheme of fuzzy inference system

ZHANG Dong-wei1，XIONG Xue-qiang1，ZHAO Kui1，ZHONG Lin-chuan2

(1.School of Resource and Environmental Engineering,Jiangxi University of Science and Technology,Ganzhou 341000,China；2.Zijin Mining Co.Ltd.，Shanghang 364200,China.)

This paper studies the stability coefficients of the mined out area by using the fuzzy relation matrix and joint stability analysis goaf under the premise of using the fuzzy relation matrix associated with the goaf stability analysis method.Combined with the stability of mined-out area and ground pressure,goaf treatment options of fuzzy inference systems are obtained.

treatment of mined goaf;fuzzy reasoning；complicated goaf；stability of goaf

TF065.1；TG111.3

A

1674-9669（2012）01-0089-06

2011-11-12

國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃資助項(xiàng)目（2011cb411913）

張東煒（1985- ），男，碩士研究生，主要從事巖石力學(xué)方面研究，E-mail：zhangdongdelove@yahoo.com.cn.