王 偉 李云安 王欽剛 胡麗珍 雷 銀 鄒濟(jì)韜

(中國地質(zhì)大學(xué)工程學(xué)院，湖北 武漢 430000)

上向水平分層充填法采礦地表變形研究

王 偉 李云安 王欽剛 胡麗珍 雷 銀 鄒濟(jì)韜

(中國地質(zhì)大學(xué)工程學(xué)院，湖北 武漢 430000)

開采地下礦產(chǎn)會(huì)造成地表開裂、變形，破壞建筑物，對(duì)人們的生產(chǎn)、生活造成嚴(yán)重的影響，礦區(qū)地表變形問題一直是影響礦山安全生產(chǎn)的一個(gè)重要方面。應(yīng)用FLAC3D數(shù)值分析軟件，建立三維地質(zhì)模型，模擬計(jì)算司家營鐵礦分步開采回填后相應(yīng)地表變形情況，并與經(jīng)驗(yàn)公式計(jì)算的地表沉降值進(jìn)行對(duì)比分析。提出了司家營礦區(qū)地表變形的一般規(guī)律及變形機(jī)制，對(duì)地表變形的地質(zhì)災(zāi)害做出分析評(píng)價(jià)。研究結(jié)果表明，上向水平分層填充法采礦能有效的控制地表變形問題，填充體對(duì)減小變形具有重要的作用。對(duì)比分析表明數(shù)值模擬結(jié)果比經(jīng)驗(yàn)公式結(jié)果偏小，但在合理許可范圍內(nèi)，表明FLAC3D是預(yù)測地表變形的一種有效手段,對(duì)其他相似采礦引起的地表變形研究和預(yù)測具有重要的參考價(jià)值。

上向水平分層充填法 地表變形 數(shù)值模擬 經(jīng)驗(yàn)公式

隨著國內(nèi)外資源日趨減少，而國民經(jīng)濟(jì)的迅速發(fā)展對(duì)鋼鐵需求的不斷增長。由于開采深度增加，環(huán)境保護(hù)要求嚴(yán)格，充填采礦法適應(yīng)性強(qiáng)，礦石回采率高，貧化率低，作業(yè)較安全，能利用工業(yè)廢料，保護(hù)地表等，充填采礦法的應(yīng)用比例逐年增加。上向水平分層充填法是自下而上分層回采，每分層先采出礦石，而后填入充填料，以支撐采空區(qū)兩幫和作為工作平臺(tái)。該方法為工作面循環(huán)作業(yè)，鑿巖爆破、出礦、充填和護(hù)頂完成一個(gè)循環(huán)后，進(jìn)行下一分層的循環(huán)；回采空間和范圍可以控制，可以用任何充填材料進(jìn)行充填[1-3]。

礦區(qū)地表變形問題一直是礦山安全生產(chǎn)的一個(gè)重要方面，研究人員運(yùn)用很多方法研究這一問題，常用的研究方法有現(xiàn)場監(jiān)測、理論研究、物理模擬及數(shù)值模擬等方法，每種方法有其各自的特點(diǎn)。礦區(qū)的地表變形是一個(gè)非常復(fù)雜的問題，影響因素有很多，包括礦山的開采方法、地層產(chǎn)狀、巖土體的物理力學(xué)性質(zhì)、結(jié)構(gòu)構(gòu)造、地下水等等。

近年來，數(shù)值模擬方法應(yīng)用非常廣泛，它可以模擬任意尺度的地質(zhì)體，有非常強(qiáng)的可操作性和可重復(fù)性，這是其他研究方法所不具有的特點(diǎn)。對(duì)于本研究的地表變形問題，數(shù)值模擬可以很好的模擬真實(shí)開挖填充的全過程，把影響范圍內(nèi)的巖土體作為一個(gè)整體進(jìn)行分析，獲得普遍適用的規(guī)律。

本研究采用有限差分方法，選擇由美國Itasca公司開發(fā)的大型巖土工程計(jì)算軟件FLAC3D建立三維數(shù)值模擬，研究上向分層水平充填法采礦引起的地表變形規(guī)律及機(jī)制。

1 鐵礦工程地質(zhì)概況

司家營南區(qū)位于華北臺(tái)塊東北端次級(jí)構(gòu)造單元—燕山臺(tái)褶皺帶的古隆起東南部位。礦區(qū)地層主要為前震旦系、震旦系和第四系地層。司家營鐵礦為一大型鞍山式沉積變質(zhì)鐵礦床，以S6勘探線為界，分南北兩區(qū)，總面積約為36 km2，南區(qū)為大賈莊礦段。選取大賈莊礦段Ⅰ礦體為數(shù)值模擬對(duì)象，位于S66線以北，S58線以南。該段礦體均為第四系覆蓋之下的隱伏礦體。礦體主要賦存于黑云變粒巖中，產(chǎn)狀與圍巖一致。礦體走向近南北，傾向西，傾角30°～60°。礦石普氏硬度系數(shù)f=12～16。區(qū)內(nèi)構(gòu)造不發(fā)育。礦區(qū)地表?xiàng)l件較為復(fù)雜，新河位于礦區(qū)東側(cè)并穿過礦區(qū)上部。礦區(qū)地表為村莊和農(nóng)田，村莊比較密集。采礦條件復(fù)雜對(duì)研究采礦引起的地表變形具有非常重要的影響。

綜合礦體賦存條件、開采條件和開發(fā)建設(shè)環(huán)境、要求等因素，設(shè)計(jì)確定采用地下開采，采礦方法采用上向水平分層充填采礦法。上向水平分層充填法的優(yōu)點(diǎn)主要是對(duì)礦體賦存條件適應(yīng)性好，夾石易剔除混入率較低，上盤暴露面積小并且容易控制。本研究采用現(xiàn)場勘查分析、數(shù)值模擬研究和經(jīng)驗(yàn)公式計(jì)算分析等綜合手段研究了地表規(guī)律，為沉降變形預(yù)測及防治提供了重要的理論依據(jù)。

2 建立三維數(shù)值模型

2.1 計(jì)算模型的建立

基于司家營鐵礦大賈莊Ⅰ礦體的開采設(shè)計(jì)，為了實(shí)現(xiàn)多階段開采的數(shù)值模擬，本研究應(yīng)用ANSYS 建立礦山三維數(shù)值模型并劃分網(wǎng)格，然后使用FLAC3D-ANSYS 接口程序?qū)氲紽LAC3D中進(jìn)行數(shù)值模擬。從下到上建立了8個(gè)盤區(qū)，模擬依次開挖回填過程，盤區(qū)內(nèi)沿走向每隔100 m設(shè)4 m厚間柱，在盤區(qū)沿走向平均劃分為4個(gè)采場，每個(gè)采場長48 m。根據(jù)礦體開挖范圍，考慮數(shù)值計(jì)算模型的邊界效應(yīng)，以礦體傾向方向?yàn)閄軸方向，總長800 m; 以礦體走向方向?yàn)閅軸方向，取一個(gè)盤區(qū)的長度，即200 m; 以豎直方向?yàn)閆軸方向，總高450 m，地表高程450 m。整個(gè)三維模型的尺寸定為長×寬×高=800 m×200 m×450 m，具體模型見圖1。模型從上到下依次為地表第四系、強(qiáng)風(fēng)化巖體、弱風(fēng)化巖體、礦體及圍巖，此外模型中間還包括開挖礦體，具體模型見圖2。模型總共劃分51 960個(gè)單元，56 616個(gè)節(jié)點(diǎn)。為了研究地表變形的規(guī)律及機(jī)制分析，通過FLAC3D具有的跟蹤功能，地表中線處每40 m布置2個(gè)監(jiān)測點(diǎn)，分別監(jiān)測垂直位移及水平位移，共布置監(jiān)測點(diǎn)42 個(gè)。對(duì)這些點(diǎn)的位移情況進(jìn)行全過程的追蹤記錄[4-6]。

圖1 ANSYS三維數(shù)值計(jì)算模型

圖2 工程地質(zhì)分組模型

2.2 計(jì)算參數(shù)的選取

計(jì)算將采用摩爾-庫侖(Mohr-Coulomb)屈服準(zhǔn)則，在查明了礦區(qū)工程地質(zhì)條件的基礎(chǔ)上，結(jié)合現(xiàn)場勘探資料和以往地質(zhì)資料，巖土體和充填體力學(xué)參數(shù)采用表1中所確定的物理力學(xué)參數(shù)。

表1 巖土體的物理力學(xué)參數(shù)

該三維數(shù)值計(jì)算模型的幾點(diǎn)假設(shè)：①采用彈塑性物理模型，選取摩爾-庫侖(Mohr-Coulomb)屈服準(zhǔn)則；②假設(shè)巖層為均質(zhì)、各向同性連續(xù)介質(zhì)，不考慮動(dòng)荷載及復(fù)雜應(yīng)力的作用；③根據(jù)水壓致裂地應(yīng)力測量原位試驗(yàn)得到的應(yīng)力值，能較好的模擬真實(shí)應(yīng)力情況，前后、左右邊界施加梯度應(yīng)力邊界條件，下邊界約束水平X、Y向和豎直Z向位移，上邊界(地面)為自由面；④模擬計(jì)算首先讓模型在應(yīng)力、位移邊界條件作用下穩(wěn)定,建立初始自重應(yīng)力場,然后清除初始應(yīng)力產(chǎn)生的位移,逐步模擬開采回填過程，記錄每一次開挖回填后監(jiān)測點(diǎn)的數(shù)據(jù)。

3 計(jì)算結(jié)果分析

根據(jù)采礦方案設(shè)計(jì)，對(duì)該模型礦體進(jìn)行8步開挖回填，通過模擬得到應(yīng)力場、位移場，研究開挖回填過程和填充體對(duì)地表變形的影響。

通過地表設(shè)置的42個(gè)監(jiān)測點(diǎn)得到的數(shù)據(jù)，繪制各監(jiān)測點(diǎn)垂直位移及水平位移變化曲線，如圖3、圖4所示。約定垂直位移以向上拱起為正，以向下沉降為負(fù)；水平位移以向右為負(fù)，以向左為正。

圖3 地表監(jiān)測點(diǎn)垂直位移變化曲線

圖4 地表監(jiān)測點(diǎn)水平位移變化曲線

由圖3可得，垂直變形曲線為V字形。隨著開采回填的進(jìn)行，沉降范圍逐漸擴(kuò)大，沉降量也逐步增大，最大沉降量為325.3 mm，位于地表水平距離320 m處；而拱起變形范圍逐漸減小，沉降量也慢慢減小，最終變形值為39 mm。因礦體傾斜的緣故，沉降中心位于上盤靠近礦體的地方，隨著開采的進(jìn)行，其沉降中心也發(fā)生變化，逐漸向右移動(dòng)。上盤的變形值和變形范圍大于下盤。

由圖4可得，水平變形曲線為正弦曲線形狀。地表水平位移上盤明顯大于下盤，因?yàn)榈V體傾角的因素，上盤位移大約是下盤位移的2倍，位移值隨開挖深度的減小而增大。向左變形最大值為261 mm，位于地表220 m處；向右變形最大值為138 mm，位于地表520 m處。開采深度對(duì)地表變形特征有很大的影響，隨著開采深度的減少，地表的各種變形隨之增大。

4 經(jīng)驗(yàn)公式計(jì)算分析

依據(jù)該礦區(qū)的覆巖巖性、地質(zhì)構(gòu)造以及司家營鐵礦開采方案等情況，參考《建筑物、水體、鐵路及主要井巷煤柱留設(shè)與壓煤開采規(guī)程》附錄五中類似礦區(qū)的地表移動(dòng)參數(shù)，選取合適的計(jì)算參數(shù)，按下列經(jīng)驗(yàn)公式估算：

式中,q為下沉系數(shù),mm/m，取0.08；Mc為計(jì)算開采厚度,m，

hm為充填前頂?shù)装宓囊平浚琺，按類似礦區(qū)經(jīng)驗(yàn)，取0.15；hf為充填的不充分量,m，按類似礦區(qū)經(jīng)驗(yàn)，取0.1；B為沉縮率，按全部充填后考慮充填材料的沉縮率，取0.03；M為累計(jì)采厚，m，按主要開采礦體厚度累加，取185 m；α為礦層傾角,(°)，取40°。

代入相關(guān)參數(shù)，計(jì)算得最大下沉值為

三維數(shù)值計(jì)算結(jié)果值與經(jīng)驗(yàn)公式計(jì)算值結(jié)果相比后，數(shù)值模擬的結(jié)果小于理論計(jì)算的結(jié)果。經(jīng)分析原因有以下幾點(diǎn)原因：①只模擬了開采過程中的彈塑性變形過程，沒有考慮殘余變形及時(shí)間效應(yīng)；②假設(shè)巖層為均質(zhì)、各向同性連續(xù)介質(zhì)，忽略地質(zhì)構(gòu)造的影響；③在計(jì)算模型的建立中，存在一定的邊界效應(yīng)，會(huì)產(chǎn)生一定的誤差；④使用的是巖石力學(xué)參數(shù)，而不是巖體力學(xué)參數(shù)。

礦山地表變形是一個(gè)復(fù)雜的過程，它與開采、地下水、充填質(zhì)量及爆破擾動(dòng)等因素有關(guān)，它是各種影響因素的外在表現(xiàn)。數(shù)值模擬忽略一些不重要的因素，模擬值與理論值有一定的差距，但兩者差距不大，是比較合理的。取兩者的平均值作為最終的最大沉降值：

地表的破壞及產(chǎn)生地裂縫很大程度上取決于地表傾斜(icm)、地表變形曲率(kcm)、最大水平移動(dòng)值(Ucm)、地表水平變形值(εcm)的大小。

最大傾斜值：

最大曲率值：

最大水平移動(dòng)值：

最大水平變形值：

式中,r為地表影響區(qū)半徑(m)，按式r=H/tanβ計(jì)算；H為開采深度(m)，取425；Β為移動(dòng)角(°)，取70°；b為水平移動(dòng)系數(shù)，取0.33；

計(jì)算結(jié)果如下：

最大傾斜值：

icm=2.19mm/m.

最大曲率值：

最大水平移動(dòng)值：

Ucm=113.3mm.

最大水平變形值：

從上述計(jì)算結(jié)果分析，地面沉降、最大傾斜值、最大曲率值、最大水平移動(dòng)值、最大水平變形值均小于允許值，對(duì)地面建筑影響較小。

充填法采礦具有其不同于其他采礦方向的特征，由于充填了采礦區(qū)，限制了圍巖的變形，調(diào)整圍巖應(yīng)力狀態(tài)，使圍巖應(yīng)力集中現(xiàn)象消失，具有一定的強(qiáng)度和抵抗變形的能力，起到一定轉(zhuǎn)移和傳遞變形的作用，從而減小了地表變形[7-8]。

5 結(jié) 論

(1)研究結(jié)果表明，司家營大賈莊礦段采用上向水平分層充填法采礦能較好的控制地表位移，效果顯著，對(duì)地面建筑物和農(nóng)田影響較小。

(2)采用填充法開采的傾斜礦體，沉降中心位于上盤靠近礦體的地方，隨著開采的進(jìn)行，其沉降中心也發(fā)生變化。

(3)填充體的存在調(diào)整圍巖應(yīng)力狀態(tài)，使圍巖應(yīng)力集中現(xiàn)象消失，具有一定的強(qiáng)度和抵抗變形的能力，起到一定轉(zhuǎn)移和傳遞變形的作用。

(4)數(shù)值模擬沉降值比經(jīng)驗(yàn)公式計(jì)算值偏小，但在合理的范圍內(nèi)，表明FLAC3D數(shù)值模擬能較好的預(yù)測地表變形的趨勢(shì)。

[1] 胡 華，孫恒虎.礦山充填工藝技術(shù)的發(fā)展及似膏體充填新技術(shù)[J] .中國礦業(yè)，2001(6)：47-50. Hu Hua，Sun Henghu.Development of backfill technology and the new backfill progress using paste like material [J].China Mining Magazine，2001(6)：47-50.

[2] 解世俊.金屬礦床地下開采[M] .北京：冶金工業(yè)出版社,1986. Xie Shijun.Underground Mining on Metal Deposits [M] .Beijing：Metallurgical Industry Press，1986.

[3] 袁世倫，胡國斌，楊承祥.金屬礦山充填技術(shù)的回顧與展望[J] .江西有色金屬，2004(3):11-15. Yuan Shilun，Hu Guobin，Yang Chengxiang.Review and exhibition of metal mine filling technology hope [J].Jiangxi Nonferrous Metal，2004(3):11-15.

[4] 陳育民，徐鼎平.FLAC /FLAC3D基礎(chǔ)與工程實(shí)例[M].北京:中國水利水電出版社,2009. Chen Yumin，Xu Dingping.FLAC / FLAC3DFoundation and Project Examples[M].Beijing：China Water Power Press，2009.

[5] 鄭天涯，周建軍.基于FLAC方法的煤礦開采沉降分析[J].三峽大學(xué)學(xué)報(bào):自然科學(xué)版,2010,32(4)：63-65. Zheng Tianya，Zhou Jianjun.Settlement analysis of coal well mining based on FLAC [J].Journal of China Three Gorges University：Natural Science，2010,32 (4)：63-65.

[6] 廉 海，魏 秀，甘德清.地下開采引起地表沉陷的數(shù)值模擬[J].礦業(yè)快報(bào),2006(1):29-32. Lian Hai，Wei Xiu，Gan Deqing.Numerical simulation of ground surface subsidence caused by underground mining [J].Express Information of Mining Industry，2006 (1)：29-32.

[7] 謝曉龍.采空老區(qū)地表變形及沉陷機(jī)理研究[D].合肥:合肥工業(yè)大學(xué),2010. Xie Xiaolong.Ground Deformation and Subsidence Mechanism of Old Mined Area[D].Hefei：Hefei University of Technology,2010.

[8] 黃路平，陳從新，肖國峰，等.復(fù)雜地質(zhì)條件下礦山地下開采地表變形規(guī)律的研究[J].巖土力學(xué)，2009，30(10)：3020-3024. Huang Luping，Chen Congxin，Xiao Guofeng，et al.Study of rock movement caused by underground mining in mines with complicated geological conditions [J].Rock and Soil Mechanics，2009，30(10)：3020-3024.

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

Analysis of Surface Deformation Caused by Upward Horizontal Slicing and Filling Stoping Method

Wang Wei Li Yun′an Wang Qingang Hu Lizhen Lei Yin Zou Jitao

(EngineeringCollege，ChinaUniversityofGeoscience，Wuhan430000，China)

Exploration of underground minerals usually results in surface cracking，surface deformation，and damage to buildings，which has a serious impact on people′s production and life.Surface deformation in mining area is always an important aspect affecting mining safety.FLAC3Dsoftware was applied to build a three-dimensional geological model for simulating the surface deformation caused by the upward horizontal slicing and filling stoping in Sijiaying Iron Mine and a contrast with the surface subsidence value calculated by the empirical formula was made.Thus，the general rule and the deformation mechanism concerning surface deformation in Sijiaying Iron Mine were put forwards，and the geological hazards of surface deformation were assessed.The results showed that the upward horizontal slicing and filling stoping method effectively controlled surface deformation and the filling body played an important role in reducing deformation.The contrastive analysis showed that the numerical simulation result is smaller than that by the empirical formulas，but still within a reasonable limit permitted，indicating that FLAC3Dis an effective way of predicting surface deformation.The research has important reference value for research and prediction of surface deformation caused by other similar mining method.

Upward horizontal slicing and filling stoping method，Surface deformation，Numerical simulation，Empirical formula

2014-03-09

王 偉(1989—)，男，碩士研究生。

TD325

A

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