郭生茂 劉 濤 程 剛 嚴(yán)文炳

(1.西北礦冶研究院，甘肅 白銀 730900；2.新疆工程學(xué)院安全工程系，新疆 烏魯木齊 830091)

基于RG法的空區(qū)群礦柱臨界破壞條件研究

郭生茂1劉 濤1程 剛2嚴(yán)文炳1

(1.西北礦冶研究院，甘肅 白銀 730900；2.新疆工程學(xué)院安全工程系，新疆 烏魯木齊 830091)

采空區(qū)群的穩(wěn)定性對礦山的安全生產(chǎn)有重大影響，針對礦山采空區(qū)賦存形態(tài)，建立了對采空區(qū)群RG一維和二維模型，考慮應(yīng)力重分布對系統(tǒng)穩(wěn)定性的影響，將采空區(qū)群系統(tǒng)失穩(wěn)破壞的復(fù)雜問題轉(zhuǎn)換為了單元破壞失穩(wěn)概率與臨界失穩(wěn)概率之間比較的問題，并給出了礦山采空區(qū)群一維和二維狀態(tài)下的臨界失穩(wěn)概率，運用該方法對七角井鐵礦采空區(qū)I和采空區(qū)II進(jìn)行了計算驗證，計算結(jié)果顯示與實際情況較為符合，說明該方法對評價礦山采空區(qū)穩(wěn)定性有一定的實用性。

RG法 采空區(qū)群 臨界概率 評價

隨著淺部資源的不斷開采、消耗，礦山逐漸向深部發(fā)展，而淺部開挖形成了大量的采空區(qū)，空區(qū)賦存條件復(fù)雜，且大部分的礦山未進(jìn)行充填，為地表安全、深部開采造成了重大安全影響[1-3]。在20世紀(jì)70—90年代，我國大部分礦山通常采用留礦法、空場法等采礦方法，礦山存在大量的采空區(qū)群，大部分的采空區(qū)都未進(jìn)行充填處理，從而帶來安全隱患[4]。從以往的采空區(qū)失穩(wěn)破壞規(guī)律來看，采空區(qū)塌陷主要是因為支撐空區(qū)穩(wěn)定性的礦柱發(fā)生失穩(wěn)造成的，礦柱作為支撐空區(qū)穩(wěn)定性的基本單元，其穩(wěn)定直接影響到整個空區(qū)系統(tǒng)的整體穩(wěn)定性[5-7]?？諈^(qū)群的破壞常伴有多米諾骨牌效應(yīng)，常常造成巨大的人員傷亡和經(jīng)濟損失[8]。本研究基于重整化群(RG)法，建立了采空區(qū)群RG一維和二維模型，對七角井鐵礦兩大采空區(qū)群進(jìn)行了計算分析，計算結(jié)果與現(xiàn)場實際較為吻合，從整體分析的角度為礦山采空區(qū)群安全評價及穩(wěn)定性問題提供了一個新的思路。

1 基于RG法的礦柱臨界破壞分析

RG法最早應(yīng)用于量子場論中，后來被美國著名科學(xué)家Wilson將該方法應(yīng)用到臨界現(xiàn)象的研究中，并取得了巨大成功，Wilson也因此獲得1982年的諾貝爾物理學(xué)獎。重整化群法通過改變觀測尺寸定量獲取物理量的變化，該方法通過對系統(tǒng)內(nèi)部的基本組成單元進(jìn)行一系列的自相似變換，從而獲取宏觀系統(tǒng)的特性描述。通過計算重整化群不動點與臨界相中的關(guān)聯(lián)長度無窮大相聯(lián)系，從而計算系統(tǒng)的臨界條件，該方法已成為計算和分析系統(tǒng)臨界行為強有力的分析方法。

1.1 采空區(qū)群礦柱穩(wěn)定性演化一維RG模型

礦山采空區(qū)群是一個復(fù)雜的時空演變系統(tǒng)，空區(qū)間通過預(yù)留礦柱來支撐采空區(qū)的穩(wěn)定性，采空區(qū)群內(nèi)單個礦柱的失穩(wěn)破壞必定會影響到整個系統(tǒng)的應(yīng)力重新分布，應(yīng)力重新分布的結(jié)果會使某些未發(fā)生破壞失穩(wěn)的礦柱應(yīng)力集中加劇，從而導(dǎo)致次一級的礦柱發(fā)生失穩(wěn)破壞，這樣的連鎖反應(yīng)依次發(fā)生，最終導(dǎo)致整個采空區(qū)群發(fā)生坍塌[9]。

根據(jù)重整化群理論，首先建立礦山采空區(qū)群原包模型，如圖1所示，在一維情況下，一個一級原包有兩個基本單元組成，一個礦柱和相鄰采空區(qū)組成一個基本單元。

圖1 空區(qū)群一維重整化群簡化模型

圖2 2個單元組成一級原包的所有可能情況

如果考慮一個礦柱破壞后對臨近礦柱應(yīng)力轉(zhuǎn)移重新分布的影響，假設(shè)Pab為一個礦柱破壞后向鄰近礦柱應(yīng)力轉(zhuǎn)移的條件概率，且應(yīng)力向鄰近未破壞的礦柱轉(zhuǎn)移概率是相等的，那么，可得一級原包發(fā)生失穩(wěn)破壞的概率為

(1)

其中，應(yīng)力轉(zhuǎn)移條件概率

(2)

式中，Pa為礦柱未破壞前鄰近礦柱的失穩(wěn)概率；Pb為礦柱破壞后鄰礦柱的失穩(wěn)概率。

假設(shè)基本單元的破壞概率由Weibull分布給出，即：

(3)

式中，F(xiàn)0為基本單元的平均承載力；m為反應(yīng)原包強度的相似性的相關(guān)參數(shù)。

可得：

(4)

當(dāng)m=2時，可得一級原包的失穩(wěn)破壞概率為

(5)

根據(jù)logistic一維映射遞推關(guān)系，可得到n級原包的破壞失穩(wěn)概率為

(6)

式(6)即為一維條件下的重整化群變換關(guān)系方程，為了便于分析，將式(6)寫成函數(shù)形式，即

(7)

在0≤x≤1的范圍內(nèi)求解式(6)，圖3為Pn與Pn-1的重整化映射關(guān)系曲線，從圖3中可以看出，在0～1的范圍內(nèi)，存在3個不動點，即P*=0、0.206、1這3個不動點。

圖3 Pn與Pn-1重整化映射曲線

根據(jù)判斷準(zhǔn)則，當(dāng)

那么x為不穩(wěn)定的不動點，反之，當(dāng)

x為穩(wěn)定的不動點，分別將x=0、0.206、1代入判據(jù)(df(x))/(dx)分別為0、1.62、0，則說明0和1為穩(wěn)定不動點，0.206為不穩(wěn)定不動點。當(dāng)P*=0時，表示系統(tǒng)內(nèi)沒有原包破壞，系統(tǒng)處于穩(wěn)定狀態(tài)，當(dāng)P*=1，系統(tǒng)內(nèi)所有原包都被破壞，那么P*=0.206可以將系統(tǒng)分為2種狀態(tài)，當(dāng)P*<0.206，即系統(tǒng)趨近于P*=0，當(dāng)P*>0.206，系統(tǒng)將向P*=1的趨勢發(fā)展，因此，P*=0.206為系統(tǒng)的臨界點。

1.2 采空區(qū)群礦柱穩(wěn)定性演化二維RG模型

首先建立二維重整化群礦山采空區(qū)群重整化原包模型，如圖4所示。

圖4 空區(qū)群二維重整化群簡化模型

圖5 4個單元組成一級原包的所有可能情況

考慮破壞失穩(wěn)單元對鄰近單元應(yīng)力重分布的影響，運用一維情況下同樣的求解方法，可得一級原包失穩(wěn)概率為

(8)

根據(jù)logistic一維映射遞推關(guān)系，可得到n級原包的破壞失穩(wěn)概率為

(9)

結(jié)合matlab求解上述方程，可得P*=0.170 7。

2 工程應(yīng)用

從上述的分析可以看出，對于礦山采空區(qū)群系統(tǒng)的整體穩(wěn)定性可以歸納為單元體失穩(wěn)概率與臨界失穩(wěn)概率的比較，當(dāng)單元體破壞概率大于臨界概率，則存在整體失穩(wěn)的可能，反之，則無。

單元的破壞失穩(wěn)概率可用礦柱承受的平均應(yīng)力與礦柱的承載強度的比值來近似表示，礦柱承載平均應(yīng)力按照公式(10)計算，礦柱的承載強度如按公式(11)計算。

(10)

(11)

(12)

式中，σS為礦柱承受的平均應(yīng)力；σR為礦柱的承載強度；σc為巖石的單軸抗壓強度；Vc,Vm為巖石彈性波速和巖體波速；a、b、c為空區(qū)跨度、礦柱長度和寬度，具體見圖1所示；H為空區(qū)上覆蓋巖層的厚度；h為礦柱的高度。

七角井鐵礦位于甘肅北部，礦山初期為小規(guī)模露天開采，后轉(zhuǎn)入地下開采，礦山自開采以來，對采空區(qū)未進(jìn)行充填處理，形成了大量的采空區(qū)，對礦山的安全生產(chǎn)帶來了隱患，為防止采空區(qū)坍塌對工人安全和礦山的持續(xù)生產(chǎn)產(chǎn)生影響，需對已形成的采空群進(jìn)行安全評價和治理。表1為七角井鐵礦淺部(采空區(qū)群I)和深部采空區(qū)群(采空區(qū)群II)基本參數(shù)。

表1 兩采空區(qū)形態(tài)和巖石力學(xué)參數(shù)

根據(jù)RG法和式(10)～式(12)，對于采空區(qū)群I分布較廣且數(shù)量大，分布較為不規(guī)則，故采用二維RG對其進(jìn)行分析，采空區(qū)群II由于層間頂柱較厚，故采用一維RG模型對其進(jìn)行分析。

對于采空區(qū)Ⅰ：

說明采空區(qū)I不穩(wěn)定狀態(tài)，存在整體失穩(wěn)的可能，從現(xiàn)場調(diào)查的情況來看，礦區(qū)黃鐵礦受到分化影響，節(jié)理發(fā)育，巖體較為破碎，露天采場下的多處空區(qū)產(chǎn)生坍塌，與理論分析較為符合。

對于采空區(qū)Ⅱ：

說明采空區(qū)Ⅱ是穩(wěn)定的，不存在整體破壞失穩(wěn)的可能，從現(xiàn)場調(diào)查的情況來看，深部采空區(qū)較為穩(wěn)定，沒有出現(xiàn)大范圍的坍塌，整體性良好，與理論分析結(jié)果一致。

3 結(jié) 論

(1)運用重整化群方法，建立了對采空區(qū)群RG一維和二維模型，考慮應(yīng)力重分布對系統(tǒng)穩(wěn)定性的影響，給出了礦山采空區(qū)群一維和二維狀態(tài)下的臨界失穩(wěn)概率。

(2)對七角井鐵礦空區(qū)I和空區(qū)群Ⅱ進(jìn)行了計算驗證，計算結(jié)果顯示與實際情況較為符合，說明該方法對評價礦山采空區(qū)穩(wěn)定性有一定的實用性。

(3)由于礦山采空區(qū)本身的復(fù)雜性及其在礦山生產(chǎn)過程中，采空區(qū)常常受到爆破等動力擾動，式(10)與式(11)是通過大量經(jīng)驗總結(jié)出來的公式，因此，應(yīng)用上述方法對礦山采空區(qū)進(jìn)行評價需結(jié)合工程實際，靈活選用礦柱強度的計算公式。

[1] 孫國權(quán)，李 娟，胡杏保．基于FLAC3D程序的采空區(qū)穩(wěn)定性分析[J]．金屬礦山，2007(2)：29-32． Sun Guoquan，Li Juan，Hu Xingbao.FLAC3D-based stability analysis of mined-out area[J].Metal Mine，2007(2):29-32.

[2] 陳建功，周陶陶，張永興．深部洞室圍巖分區(qū)破裂化的沖擊破壞機制研究[J]．巖土力學(xué)，2011，32(9)：2629-2634． Chen Jiangong，Zhou Taotao，Zhang Yongxing.Shock failure mechanism of zonal disintegration within surrounding rock in deep chamber[J].Rock and Soil Mechanics，2011，32(9):2629-2634.

[3] 徐水太，饒運章，潘建平．地下開采引起地表移動和不均勻沉降的機理分析[J]．礦業(yè)工程，2004，56(1)：31-33． Xu Shuitai，Rao Yunzhang，Pan Jianping.Mechanism of ground surface displacement & uneven subsidence caused by underground mining[J].Mining Engineering，2004，56(1):31-33.

[4] 周科平，鄧環(huán)宇，楊念哥．復(fù)雜空區(qū)群對周邊地質(zhì)體環(huán)境影響研究[J]．礦冶工程，2011，31(3)：1-4． Zhou Keping，Deng Huanyu，Yang Niange.Study on effect of complicated goaf group on surrounding geologic-body environment[J].Mining and Metallurgical Engineering，2011，31(3):1-4.

[5] 潘 岳，張 勇，吳敏應(yīng)，等．非對稱開采礦柱失穩(wěn)的突變理論分析[J]．巖石力學(xué)與工程學(xué)報，2006(S2)：3694-3702． Pan Yue，Zhang Yong，Wu Minying，et al.Analysis of catastrophe theory for pillar destabilization in Dissymmetric Mining[J].Chinese Journal of Rock Mechanics and Engineering，2006(S2):3694-3702.

[6] 劉洪強，張欽禮，潘常甲，等．空場法礦柱破壞規(guī)律及穩(wěn)定性分析[J]．采礦與安全工程學(xué)報，2011(1)：138-143． Liu Hongqiang，Zhang Qinli，Pan Changjia，et al.Analysis of the failure law and stability of the pillar in open stope mining[J].Journal of Mining & Safety Engineering，2011(1):138-143.

[7] 韓 穎，楊志龍．礦柱失穩(wěn)的突變分析[J]．中國礦業(yè)，2013，22(1)：86-88． Han Ying，Yang Zhilong.Pillar CUSP catastrophe analysis of instability[J].China Mining Magazine，2013，22(1):86-88

[8] 麥尚德，王洛鋒，李迎峰，等．用深孔爆破成井次品礦充填技術(shù)治理采空區(qū)[J]．金屬礦山，2012(5)：27-31． Mai Shangde，Wang Luofeng，Li Yingfeng，et al.Mined-out area treatment with the deep-hole blasting shaft-forming method and filling technique with poor ore[J].Metal Mine，2012(5):27-31.

[9] 周子龍，李夕兵，趙國彥．民窿空區(qū)群級聯(lián)失穩(wěn)評價[J]．自然災(zāi)害學(xué)報，2007(5)：91-95． Zhou Zilong，Li Xibing，Zhao Guoyan.Evaluation of cascading collapse for private mined zone groups[J].Journal of Natural Disasters，2007(5):91-95.

(責(zé)任編輯 石海林)

Pillar Critical Failure Condition of Mined-out Zone Groups Based on RG Method

Guo Shengmao1Liu Tao1Cheng Gang2Yan Wengbing1

(1.NorthwestInstituteofMiningandMetallurgy，Baiyin730900，China;2.DepartmentofSafetyEngineering,XinjiangInstituteofEngineering，Urumqi830091，China)

The mined-out zone groups have become one of the most potential safety factors for mines.In view of the mined-out zone′s occurrence，the RG one-dimensional and two-dimensional models were established.Considering the effect of stress redistribution on the stability of the system，the complex problem of the mined-out zone group system instability was transformed to the comparison between the unit instability probability and the critical instability probability.The critical instability probability of the mined-out zone groups under one-dimensional and two-dimensional conditions are put forwards.Mined-out zone 1 and 2 in Qijiaojing Iron Mine were calculated and verified by this method，showing that the calculation results are consistent with the actual situation.It is indicated that the method has a certain practicality to evaluate stability of mined-out zone.

RG method，Mined-out zone group，Critical probability，Evaluation

2014-03-07

郭生茂(1972—)，男，高級工程師。

TD853.391

A

1001-1250(2014)-08-125-04