孫世國(guó)，苗子臻，張玉娟，解禹

基于FLAC-3D的地下轉(zhuǎn)露天復(fù)雜采空區(qū)安全頂板厚度確定

孫世國(guó)，苗子臻，張玉娟，解禹

（北方工業(yè)大學(xué)土木工程學(xué)院，北京，100144）

為了保證礦山安全生產(chǎn)，采空區(qū)安全頂板厚度對(duì)保證安全有著重要作用。以紫金山金銅礦地下轉(zhuǎn)露天采區(qū)為例，基于FLAC-3D數(shù)值模擬方法，以不同安全頂板厚度數(shù)據(jù)，對(duì)采區(qū)不同頂板厚度穩(wěn)定性開展了研究。研究結(jié)果表明：紫金山金銅礦采區(qū)最小安全頂板厚度為45 m，開采方案頂板設(shè)計(jì)厚度為56 m時(shí)，采區(qū)頂板區(qū)域塑性區(qū)未連通，不會(huì)發(fā)生破壞，但礦柱區(qū)域塑性區(qū)有連通現(xiàn)象，需作相應(yīng)處理。

FLAC-3D；地下轉(zhuǎn)露天；采空區(qū)；安全頂板厚度

引言

近年來，隨著經(jīng)濟(jì)的快速發(fā)展，社會(huì)對(duì)礦產(chǎn)資源的需求急劇增加，礦產(chǎn)開采力度加大，許多礦山由地下開采轉(zhuǎn)為露天開采。地下轉(zhuǎn)露天開采使得礦山企業(yè)效益顯著提升的同時(shí)，也帶來了許多安全問題[1-4]。其中，地下采空區(qū)上方預(yù)留安全頂板作為控制地下開采與露天開采之間的相互影響，解決采場(chǎng)塌陷、突水等事故的重要措施，其厚度的確定十分關(guān)鍵。頂板厚度過小，容易引起頂板產(chǎn)生離層垮冒、折斷垮落，嚴(yán)重影響井下設(shè)備及工作人員的安全；反之，頂板過厚，因其回采率較低，開采難度較大，開采效益較低，易造成資源的浪費(fèi)。目前，國(guó)內(nèi)外確定安全頂板厚度的方法主要包括K.B.魯別捏依塔法、B.N.波哥留波夫法、平板梁理論推導(dǎo)法、工程類比法、厚跨比法、荷載傳遞線交匯法、數(shù)值分析法等。其中數(shù)值分析方法，因其能夠從整體上模擬復(fù)雜采空區(qū)的應(yīng)力分布規(guī)律及變形趨勢(shì)，計(jì)算效果良好，從而被廣泛應(yīng)用[5-8]。因此，本文基于FLAC-3D數(shù)值模擬方法，對(duì)紫金山金銅礦復(fù)雜采空區(qū)不同頂板厚度穩(wěn)定性進(jìn)行研究，確定了其安全頂板厚度，為礦山安全開采提供了技術(shù)支持及科學(xué)依據(jù)。

1　工程概況

紫金山金銅礦區(qū)位于福建省上杭縣城以北約26 km處。礦區(qū)初期采用露天與地下聯(lián)合開采，而后轉(zhuǎn)為露天開采，礦區(qū)露天采場(chǎng)底部設(shè)計(jì)標(biāo)高為+616 m，采場(chǎng)上部最高標(biāo)高為+1 096 m，上口面積1 020×700 m2，下口面積184×85 m2，到界邊坡最大高差480 m，屬高陡邊坡。根據(jù)礦區(qū)工程地質(zhì)條件，結(jié)合礦區(qū)多為低品位礦石的特點(diǎn)，礦區(qū)前期采用大直徑深孔階段空?qǐng)龇ㄟM(jìn)行井采，然后采用崩落法和房柱法進(jìn)行回采，應(yīng)用“采一留一”或“采一留二”回采方式，先采礦房，保留礦柱，礦柱高度50 m[9,10]。

目前已在標(biāo)高+520 m中段和+460 m中段形成多個(gè)礦房，礦房平均跨度15 m，長(zhǎng)度分為30 m、60 m、80 m、120 m四種。后期經(jīng)過露天開采，露天礦區(qū)現(xiàn)已形成單臺(tái)階高度為12 m的采場(chǎng)邊坡地形，臺(tái)階坡面角60°～70°。

2　安全頂板厚度確定的數(shù)值模擬研究

2.1模型的建立

根據(jù)礦區(qū)開采設(shè)計(jì)方案，結(jié)合現(xiàn)場(chǎng)開采情況，選取具有代表性的1-3剖面采區(qū)，依次建立了露天開采底標(biāo)高為+652 m、+640 m、+628 m、+616 m、+604 m時(shí)的地下轉(zhuǎn)露天開采模型，相應(yīng)的頂板厚度分別為92 m、80 m、68 m、56 m、44 m。模型長(zhǎng)度為2 000 m，寬度為135 m，高度為950 m，模型及網(wǎng)格劃分如圖1所示。

圖1　地下轉(zhuǎn)露天開采數(shù)值計(jì)算模型

2.2巖體參數(shù)的選取及邊界條件的確定

根據(jù)礦區(qū)地質(zhì)勘察資料及巖體試驗(yàn)結(jié)果，綜合確定各巖體力學(xué)參數(shù)如表1所示。相應(yīng)邊界條件為：前、后邊界施加Y方向位移約束，底部邊界施加X、Y、Z方向位移約束；露天采場(chǎng)底部作用荷載按1臺(tái)225 t挖掘機(jī)施工重量考慮。

表1　巖石力學(xué)參數(shù)

2.3模擬結(jié)果分析

隨著開采的進(jìn)行，露采境界線下移，采區(qū)周圍巖體在自重及其他荷載作用下產(chǎn)生移動(dòng)及變形。因巖體一旦進(jìn)入塑性狀態(tài)，其強(qiáng)度及承載力將明顯降低，因此，將圍巖塑性區(qū)的大小及貫通性作為采區(qū)穩(wěn)定性的評(píng)判標(biāo)準(zhǔn)[11]。其中頂板厚度為68 m、56 m、44 m時(shí)采區(qū)1-3剖面塑性區(qū)分布如圖2～圖4所示。

圖2　頂板厚68 m時(shí)模型塑性區(qū)分布圖

從圖中可以看出，采區(qū)塑性區(qū)主要分布在頂板中部、礦柱、露天采場(chǎng)底部邊角及采空區(qū)兩幫。當(dāng)頂板厚度為68 m時(shí)，由于拉應(yīng)力作用，頂板中部有零星塑性區(qū)出現(xiàn)；頂板厚度為56 m時(shí)塑性區(qū)分布范圍有所增加，但仍未貫通，頂板尚未發(fā)生破壞；頂板厚度為44 m時(shí)塑性區(qū)范圍顯著增加，出現(xiàn)明顯貫通，頂板發(fā)生破壞。且隨著頂板厚度的減小，頂板傳遞給空區(qū)兩幫的壓應(yīng)力減小，兩幫受壓產(chǎn)生塑性區(qū)明顯減少；隨著頂板厚度的減小，露天采場(chǎng)底部邊角因頂板自重作用而產(chǎn)生的拉應(yīng)力急劇減小，邊角處塑性區(qū)分布顯著減小。

圖3　頂板厚56 m時(shí)模型塑性區(qū)分布圖

圖4　頂板厚44 m時(shí)模型塑性區(qū)分布圖

為使安全頂板厚度精確至1 m，根據(jù)上述分析，應(yīng)用二分法原理對(duì)44 m～56 m范圍內(nèi)不同頂板厚度采區(qū)穩(wěn)定性進(jìn)行模擬分析，最終確定采區(qū)最小安全頂板厚度為45 m，此時(shí)采區(qū)塑性區(qū)仍未貫通，頂板仍未破壞，其塑性區(qū)分布圖如圖5所示。

3　結(jié)束語

通過對(duì)紫金山金銅礦復(fù)雜采空區(qū)不同頂板厚度穩(wěn)定性進(jìn)行數(shù)值模擬研究，可得出以下結(jié)論：

（1）采區(qū)塑性區(qū)主要分布在頂板中部、礦柱、露天采場(chǎng)底部邊角及采空區(qū)兩幫，且隨著頂板厚度的減小，礦柱、露天采場(chǎng)底部邊角及采空區(qū)兩幫塑性區(qū)范圍逐漸減小，頂板中部塑性區(qū)逐漸增加，不利于頂板安全。

圖5　頂板厚45 m時(shí)模型塑性區(qū)分布圖

（2）紫金山金銅礦采區(qū)最小安全頂板的厚度為45 m。露天采場(chǎng)底部設(shè)計(jì)標(biāo)高+616 m，頂板厚度為56 m時(shí)，采區(qū)頂板區(qū)域塑性區(qū)未連通，不會(huì)發(fā)生破壞，但在礦柱區(qū)域塑性區(qū)有連通現(xiàn)象，需作相應(yīng)處理。

國(guó)家自然科學(xué)基金(41172250)；國(guó)家十二五科技支撐項(xiàng)目(2012BAK09B06)；北京市創(chuàng)新團(tuán)隊(duì)提升計(jì)劃項(xiàng)目(IDHT20140501)；北京市科研基地建設(shè)－科研創(chuàng)新平臺(tái)、科研專項(xiàng)－沖擊地壓微震監(jiān)測(cè)與預(yù)警體系的構(gòu)建(XN083)；新型錨桿加固技術(shù)現(xiàn)場(chǎng)試驗(yàn)研究及研究生能力實(shí)訓(xùn)(XN107)

[1]劉獻(xiàn)華. 地下轉(zhuǎn)露天開采的主要安全隱患與技術(shù)對(duì)策[J]. 金屬礦山, 2002(7): 52-54, doi: 10.3321/j.issn:1001-1250.2002.07.017.

[2]譚卓英, 鐘文, 胡天壽, 等. 地下轉(zhuǎn)露天礦山復(fù)雜采空區(qū)探測(cè)工程實(shí)踐[J]. 金屬礦山, 2014, 43(2): 1-5.

[3]燕恩科, 姚國(guó)棟, 萬忠明. 露天礦采空區(qū)頂板最小安全厚度研究[J]. 化工礦物與加工, 2011, 40(5): 23-25, doi: 10.3969/ j.issn.1008-7524.2011.05.008.

[4]林連寶. 紫金山金礦地下轉(zhuǎn)露天開采采空區(qū)圍巖工程地質(zhì)特征及危害控制[J]. 中國(guó)礦山工程, 2003, 32(2): 10-14, doi: 10.3969/j.issn.1672-609X.2003.02.003.

[5]李廷春, 李術(shù)才, 邱祥波, 等. 三維快速拉格朗日法在安全頂板厚度研究中的應(yīng)用[J]. 巖土力學(xué), 2004, 25(6): 935-939, doi: 10.3969/j.issn.1000-7598.2004.06.020.

[6]孫國(guó)權(quán), 李娟, 胡杏保. 基于FLAC3D程序的采空區(qū)穩(wěn)定性分析[J]. 金屬礦山, 2007(2): 29-32, doi: 10.3321/j.issn:1001-1250.2007.02.009.

[7]魏青軍, 于澤森. 井下開采引發(fā)地面塌陷的定量評(píng)估方法——以新疆伊犁州煤礦為例[J]. 中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào), 2006, 17(2): 168-170, doi: 10.3969/j.issn.1003-8035.2006.02.039.

[8]艾顯恒, 謝建兵, 朱和玲. 地下礦山采場(chǎng)頂板厚度數(shù)值模擬研究[J]. 采礦技術(shù), 2009, 25(2): 30-32, doi: 10.3969/j.issn.1674-6082.2009.04.013.

[9]劉維東. 地下轉(zhuǎn)露天采空區(qū)頂板安全厚度及其采空區(qū)治理方案的研究[D]. 北京: 北方工業(yè)大學(xué), 2016.

[10]關(guān)宇. 復(fù)合開挖擾動(dòng)下邊坡穩(wěn)定性評(píng)價(jià)方法的研究[D]. 北京:北方工業(yè)大學(xué), 2013.

孫世國(guó)(1959-)，男，教授，博士后，博士生導(dǎo)師，北方工業(yè)大學(xué)礦山安全與巖土工程研究中心主任，土木工程一級(jí)學(xué)科責(zé)任教授，國(guó)家和北京市政府安全生產(chǎn)專家組成員，中國(guó)巖石力學(xué)與工程學(xué)會(huì)理事，國(guó)際工程地質(zhì)與環(huán)境協(xié)會(huì)中國(guó)專家組成員。

苗子臻，男，碩士研究生，主要從事巖土工程方面研究工作。

E-mail: 406818699@qq.com

Study on Security Coping Thickness of Complicated Goaf of Mining Converted from Underground to Surface based on FLAC-3D

SUN Shi-guo, MIAO Zi-zhen, ZHANG Yu-juan, XIE Yu
(College of Civil Engineering, North China University of Technology, Beijing, 100144, China)

The determination of security coping thickness of complicated goaf is of great significance to the mine safety production. In order to determine the security coping thicknes of the complicated goaf of ZiJinShan gold copper mine converted from underground to surface based on the FLAC-3D simulation method, the stability of goaf with different coping thickness was researched. The study results show that the minimum security coping thickness in ZiJinShan gold copper mine is 45 m. When the design of coping thickness is 56 m, the plastic zone of coping area in goaf is not connected, so that the goaf is not damaged. However, the plastic zone of pillar area is connected, so corresponding measures should be taken to ensure the safety of mining.

FLAC-3D; Mining Converted from Underground to Surface; Goaf; Security Coping Thickness

TU473

A

2095-8412 (2016) 05-837-04工業(yè)技術(shù)創(chuàng)新 URL: http://www.china-iti.com

10.14103/j.issn.2095-8412.2016.05.001

[11]萬文. 地下空區(qū)對(duì)邊坡穩(wěn)定性的影響研究[D]. 長(zhǎng)沙: 中南大學(xué), 2006, 10.7666/d.y1190362.