陳慶峰 黃文忠 程 斌 陳 濤 韓 流
(1.寶鋼資源控股(上海)有限公司,上海 200080;2.安徽皖寶礦業(yè)股份有限公司,安徽 池州 247100;3.中國(guó)礦業(yè)大學(xué)礦業(yè)工程學(xué)院,江蘇 徐州 221116)
隨著國(guó)民經(jīng)濟(jì)的不斷增長(zhǎng)和道路交通建設(shè)快速發(fā)展,對(duì)石灰?guī)r資源需求量也不斷增加,在世界范圍內(nèi)露天礦巖質(zhì)邊坡滑坡災(zāi)害頻繁發(fā)生,造成了嚴(yán)重的經(jīng)濟(jì)損失和人員傷亡,對(duì)鐵路、公路、露天礦等行業(yè)的安全運(yùn)營(yíng)影響重大[1,2]。未來(lái)一段時(shí)間內(nèi)我國(guó)石灰?guī)r供需問(wèn)題仍將呈不斷緊張的趨勢(shì)[3],使得石灰?guī)r礦山面臨著邊坡高度和邊坡角不斷加大[4,5],才能保證開(kāi)采效率最大化。
秀山礦在施工挖掘過(guò)程中頻頻出現(xiàn)斷層現(xiàn)象,局部石方脫落嚴(yán)重,斷層破碎帶形成楔形體情況明顯,且?guī)r體結(jié)構(gòu)松散,極易發(fā)生邊坡巖體傾倒破壞現(xiàn)象。誘發(fā)邊坡產(chǎn)生傾倒變形的因素眾多,程?hào)|幸等[6]指出,初始水平應(yīng)力、邊坡高度和邊坡角、巖體力學(xué)參數(shù)、結(jié)構(gòu)面發(fā)育情況等是影響巖質(zhì)邊坡傾倒變形破壞的關(guān)鍵性因素。解析分析法、相似模型試驗(yàn)法和數(shù)值分析法是反傾邊坡穩(wěn)定性分析常用的方法。安明旭等[7]基于懸臂梁模型,并考慮坡腳巖層的剪切破壞,認(rèn)為邊坡在傾角較陡、坡角較大時(shí)穩(wěn)定性最差。鄭達(dá)等[8]采用離心模型試驗(yàn)對(duì)古水水電站壩前傾倒變形體進(jìn)行開(kāi)挖,試驗(yàn)表明巖性、結(jié)構(gòu)面和采礦活動(dòng)等是發(fā)生傾倒變形的關(guān)鍵因素。隨著計(jì)算機(jī)技術(shù)的發(fā)展,有限元、有限差分、離散元等數(shù)值模擬方法在邊坡工程中應(yīng)用廣泛,如FLAC3D[9,10]、UDEC[11,12]、PFC[13,14]等。學(xué)者們對(duì)此也開(kāi)展了大量研究,提出了眾多耦合建模方法,但與真實(shí)地質(zhì)體之間仍有較大差異[15,16]。因此,建立復(fù)雜高精度的真實(shí)地質(zhì)模型是數(shù)值模擬分析的關(guān)鍵。
本研究首先對(duì)秀山礦采場(chǎng)進(jìn)行詳細(xì)的地質(zhì)調(diào)查,對(duì)楔形體彎曲傾倒破壞機(jī)制和邊坡變形機(jī)理進(jìn)行分析,然后以整個(gè)采場(chǎng)為例,介紹了Rhino6.0—FLAC3D6.0高精度耦合建模方法,對(duì)采場(chǎng)整體以及楔形體傾倒邊坡進(jìn)行數(shù)值模擬穩(wěn)定性分析。
秀山熔劑石灰?guī)r礦位于池州市南直距約24 km處,礦區(qū)位于皖南山區(qū)北緣,山勢(shì)大致呈北東向,屬丘陵低山區(qū)。山體走向?yàn)榻鼥|西向,山頂渾圓較平坦,以橢圓及長(zhǎng)條形為主,寬一般10~20 m。區(qū)內(nèi)山峰最高海拔354.84 m,一般184~338 m,最低海拔43.69 m,最大相對(duì)高差310.64 m,切割淺。礦山設(shè)計(jì)生產(chǎn)量為250萬(wàn)t/a,由于礦體分布差異,形成東、西2個(gè)獨(dú)立的露天采場(chǎng)(圖1)。東采場(chǎng)礦層?xùn)|西南北受控于 F8、F7、F3、F1斷層,礦體總體傾向 230°,傾角39°。西采場(chǎng)因受F4、F5斷層破壞,礦體沿走向呈階梯狀被切成東、中、西三段,礦層總體走向130°~310°傾向220°,傾角35°。
隨著西采場(chǎng)采剝臺(tái)階的逐步到界,采場(chǎng)邊坡受F1斷層影響,西采場(chǎng)西北幫揭露了一個(gè)寬度為10 m左右的破碎帶,不僅自身破碎、易風(fēng)化崩落,而且會(huì)影響到上覆蓋層的穩(wěn)定性。F1斷層在礦區(qū)內(nèi)呈北西走向,地表出露長(zhǎng)2 200 m,總體走向100°,傾向南,傾角約為70°,為南盤(pán)下降北盤(pán)上升的正斷層,該斷層被北東向斷層錯(cuò)斷,使西礦段呈階梯狀,東礦段落呈齒狀。
由于西采場(chǎng)北幫和西幫邊坡的下覆巖體中存在F1斷層,造成北幫和西幫邊坡形成順傾傾倒型臺(tái)階。這種斷層破碎帶楔形體傾倒破壞現(xiàn)象類(lèi)似于順層傾斜巖質(zhì)邊坡彎曲傾倒破壞,在降段開(kāi)采過(guò)程中,楔形體下部巖體極易發(fā)生折斷,后部巖體會(huì)隨著巖體軟弱結(jié)構(gòu)面發(fā)生滑移—彎曲變形;當(dāng)變形破壞較為嚴(yán)重時(shí),楔形體發(fā)生強(qiáng)烈折剪,在邊坡下部巖層發(fā)生離層反向的傾倒變形,如圖2所示。根據(jù)材料力學(xué)中梁板理論,將邊坡巖層簡(jiǎn)化成梁板力學(xué)模型,分析其垂直于坡面作用力下變形的臨界坡長(zhǎng),從而判斷發(fā)生滑移—彎曲傾倒變形破壞[7]的概率。
2.1.1 基本力學(xué)模型的建立
本研究將秀山礦北幫、西幫邊坡頂部臺(tái)階的順傾楔形體簡(jiǎn)化成圖3所示的基本梁板力學(xué)模型。假設(shè)巖層傾向兩側(cè)不受約束,邊坡的AB段發(fā)生離層且向外鼓起彎曲,BC段巖層向下滑動(dòng)。
將模型進(jìn)行簡(jiǎn)化,在巖層走向上選取單位寬度,可知巖體自重q沿著x軸方向分力q1和沿y軸方向上的分力q2分別可以表示為
式中,γ為巖層的容重,N/m3;t為滑動(dòng)巖層的厚度,m;α為巖層的傾角,(°)。
則BC段滑動(dòng)巖層的自重,分解為沿巖層滑動(dòng)的分力P可以表示為
式中,l0為滑動(dòng)段的巖層長(zhǎng)度,m;φ為層間摩擦角,(°);c為黏聚力,MPa。
2.1.2 力學(xué)機(jī)理分析
采用能量法原理對(duì)圖3中所示的順傾楔形體邊坡模型進(jìn)行計(jì)算分析,則BC段巖層的分力P所做的總功ΔW應(yīng)與巖體內(nèi)部所積蓄的變形能ΔU相互平衡,即滿(mǎn)足ΔW=ΔU。最終可得到能量平衡方程式為
式中,y′是BC段巖層自重產(chǎn)生的變形量y的一階導(dǎo)數(shù);EI(y″)2是BC段巖體的變形能的內(nèi)力,其中y″是BC段巖層自重產(chǎn)生的變形量y的二階導(dǎo)數(shù)。
假設(shè)在AB段取垂直于巖層傾向上的彈性變形方程有:
則對(duì)式(4)進(jìn)行分析,若巖層發(fā)生滑移—彎曲式破壞時(shí),其變形量y將趨于無(wú)窮大,即y→∞,則可以計(jì)算得到變形破壞時(shí)的臨界荷載:
在方程式Pcr=P中,如果令l0=0,Pcr=0,則可以計(jì)算出發(fā)生變形破壞的極限長(zhǎng)度:
當(dāng)坡長(zhǎng)超過(guò)了臨界長(zhǎng)度,即l>lcr時(shí),巖層在自重作用下,將發(fā)生滑移—彎曲式傾倒變形破壞。
由于采場(chǎng)邊坡受斷層破碎帶影響,邊坡上部巖體強(qiáng)度較大,斷層破碎帶巖體強(qiáng)度較低,導(dǎo)致在開(kāi)采過(guò)程中形成高陡邊坡,其邊坡角超過(guò)臺(tái)階下覆巖層傾角,導(dǎo)致頂部臺(tái)階形成“頭重腳輕”的楔形體,在受到降雨、爆破等影響下,極易形成楔形體傾倒現(xiàn)象。
2.2.1 巖石力學(xué)性質(zhì)
斷層破碎帶巖性是影響邊坡穩(wěn)定性的關(guān)鍵因素。通過(guò)現(xiàn)場(chǎng)勘查,楔形體傾倒區(qū)巖層發(fā)生大規(guī)模斷裂和碎石滾落,然后開(kāi)展了斷層破碎帶巖樣(由于斷層破碎帶巖性較為松散破碎和現(xiàn)場(chǎng)不便開(kāi)展原位試驗(yàn),取斷層上下盤(pán)巖樣進(jìn)行試驗(yàn),判定斷層破碎帶與周?chē)鷰r性的完整程度)的室內(nèi)單軸抗壓強(qiáng)度試驗(yàn),對(duì)其抗壓性能進(jìn)行研究。如圖4(a)為WDW-600巖石單軸試驗(yàn)機(jī),圖4(b)為現(xiàn)場(chǎng)巖樣品經(jīng)精加工后的安裝過(guò)程。
試驗(yàn)結(jié)果顯示:斷層破碎帶上下盤(pán)試樣(上盤(pán)試樣為微晶灰?guī)r、下盤(pán)試樣為灰質(zhì)白云巖)在荷載不斷增加的過(guò)程中,試樣有明顯壓密階段,變形從1.5%~2.8%左右開(kāi)始表現(xiàn)出彈性特征,直至試樣破壞。破壞過(guò)程發(fā)出巨響,試樣整體瞬間爆裂呈現(xiàn)一系列碎塊,如圖5所示。試樣的強(qiáng)度較高,彈性模量較大,殘余強(qiáng)度幾乎為零,如圖6所示,試樣呈X狀共軛斜面剪切破壞方式。總體而言,斷層破碎帶上下盤(pán)巖性較好,在采場(chǎng)開(kāi)拓延深過(guò)程,當(dāng)邊坡下部的斷層破碎帶被采動(dòng),形成了以斷層上盤(pán)滑動(dòng)為主的楔形體,當(dāng)超過(guò)一定的臨界長(zhǎng)度,邊坡斷層破碎帶將會(huì)發(fā)生滑移—彎曲式破壞現(xiàn)象。
2.2.2 結(jié)構(gòu)面發(fā)育
斷層、節(jié)理等結(jié)構(gòu)面的存在對(duì)巖質(zhì)邊坡的穩(wěn)定性起關(guān)鍵作用。受地質(zhì)構(gòu)造作用,楔形體彎曲傾倒體區(qū)域發(fā)育大量斷層破碎帶及節(jié)理密集帶。根據(jù)現(xiàn)場(chǎng)勘查,楔形體傾倒區(qū)發(fā)育1條控制性斷層F1,該斷層傾向南,傾角較陡,約70°,為一南盤(pán)下降北盤(pán)上升的正斷層。邊坡開(kāi)挖過(guò)程中,坡體發(fā)生應(yīng)力重分布,巖體處于卸荷狀態(tài),斷層、節(jié)理等結(jié)構(gòu)面逐漸張開(kāi)發(fā)育貫通形成大的滑裂面。此外巖體內(nèi)儲(chǔ)存的彈性應(yīng)變能不斷向邊坡臨空面釋放,導(dǎo)致臨空面巖體具有一定運(yùn)動(dòng)趨勢(shì)或已經(jīng)發(fā)生變形,且采礦過(guò)程中爆破振動(dòng)等對(duì)結(jié)構(gòu)面的發(fā)育貫通以及對(duì)整個(gè)滑坡的穩(wěn)定性都會(huì)產(chǎn)生不利的影響。
大量研究證明,由于算法和維度的差異性,三維模型的安全系數(shù)計(jì)算結(jié)果通常比二維的略大[17]。但是,三維地質(zhì)模型能夠較為真實(shí)地反映實(shí)際情況。因此,對(duì)于露天礦重要邊坡采用三維數(shù)值模擬計(jì)算是有必要的。本次研究真實(shí)地質(zhì)建模、數(shù)值計(jì)算采用Rhino6.0作為建模前處理,F(xiàn)LAC3D6.0作為數(shù)值計(jì)算軟件,其建模計(jì)算流程如圖7所示。
對(duì)采場(chǎng)實(shí)體地質(zhì)模型進(jìn)行邊界條件約束,采用彈塑性模型平衡初始應(yīng)力,基于Mohr-Coulom破壞準(zhǔn)則,采用強(qiáng)度折減法對(duì)采場(chǎng)邊坡進(jìn)行系數(shù)折減,并計(jì)算收斂準(zhǔn)則為不平衡力比率小于10-5的求解要求。本次數(shù)值計(jì)算參數(shù)確定采用“局部—整體—局部”的反演思路,即通過(guò)選取典型的已發(fā)生滑坡,調(diào)整由室內(nèi)試驗(yàn)獲得的巖性參數(shù),對(duì)滑坡進(jìn)行參數(shù)的反算,最終確定適用于秀山礦采場(chǎng)的巖體強(qiáng)度參數(shù)。經(jīng)過(guò)一系列的參數(shù)反演過(guò)程,最終獲得適用于FLAC3D軟件的巖體強(qiáng)度參數(shù)如表1所示。
采場(chǎng)復(fù)雜三維數(shù)值模擬計(jì)算結(jié)果如圖8所示。圖8(a)和8(b)分別為采場(chǎng)現(xiàn)狀位移云圖和應(yīng)變?cè)隽吭茍D。邊坡失穩(wěn)區(qū)域主要發(fā)生于西采場(chǎng)北幫,并且滑坡區(qū)域主要集中于北幫靠西部,其破壞范圍相對(duì)較小,主要為多臺(tái)階的局部破壞,與實(shí)際邊坡斷層破碎帶的勘察變形情況基本一致。圖8(c)和8(d)分別為西采場(chǎng)最終境界的位移云圖和應(yīng)變?cè)隽吭茍D。其破壞范圍相對(duì)于現(xiàn)狀出現(xiàn)了較大擴(kuò)展,北幫西側(cè)、中部整體失穩(wěn),并且北幫西側(cè)也出現(xiàn)大規(guī)模失穩(wěn)現(xiàn)象。因此,需要對(duì)采場(chǎng)最終境界方案做出一定調(diào)整或?qū)吰逻M(jìn)行合理有效的加固,以防止開(kāi)采過(guò)程中滑坡災(zāi)害事故的發(fā)生。
對(duì)西采場(chǎng)開(kāi)挖最終數(shù)值模型選取具有代表性的楔形體傾倒區(qū),分析其傾倒破壞機(jī)理。圖9(a)為邊坡塑性區(qū)云圖,邊坡內(nèi)部破壞區(qū)域主要以剪切破壞為主,在邊坡頂部出現(xiàn)拉張破壞。圖9(b)為邊坡總位移計(jì)算結(jié)果,變形量最大位置發(fā)生在邊坡坡面處,且由坡面向內(nèi)逐漸遞減,最大位移量已經(jīng)達(dá)到了3.37 m。圖9(c)為邊坡剪應(yīng)變?cè)隽吭茍D,其滑動(dòng)帶已經(jīng)發(fā)育貫通,且在坡腳位置剪應(yīng)變?cè)隽肯鄬?duì)集中,表明坡腳處受上覆巖層重力作用已發(fā)生的嚴(yán)重破壞,無(wú)法提供足夠的支撐力,與現(xiàn)場(chǎng)勘查坡腳處發(fā)生擠壓隆起現(xiàn)象一致。由于斷層上下盤(pán)巖性較為完整,斷層破碎帶的巖性松散破碎,巖體強(qiáng)度較差,造成邊坡滑動(dòng)面沿著斷層破碎帶貫通,當(dāng)采場(chǎng)開(kāi)拓延深至一定階段,邊坡整體彎曲傾倒失穩(wěn)破壞發(fā)生。
(1)針對(duì)秀山礦邊坡斷層破碎帶形成的楔形體彎曲傾倒破壞,根據(jù)梁板理論建立滑移—彎曲傾倒力學(xué)模型,計(jì)算分析出發(fā)生變形的臨界坡長(zhǎng),判斷其發(fā)生滑移—彎曲式傾倒變形的概率。
(2)西采場(chǎng)邊坡的上下盤(pán)巖性較為完整,在開(kāi)采過(guò)程易形成沿?cái)鄬悠扑閹缮⒖迓鋮^(qū)域,楔形體傾倒區(qū)斷層破碎帶較為發(fā)育,邊坡開(kāi)挖過(guò)程中,斷層、節(jié)理等破碎結(jié)構(gòu)面逐漸張開(kāi)發(fā)育貫通形成大的滑裂面。
(3)采用Rhino6.0—FLAC3D6.0高精度耦合建模,有效解決了FLAC3D建模與實(shí)際地質(zhì)情況不相符的情況?;诖私⒉蓤?chǎng)現(xiàn)狀和最終模型,對(duì)整個(gè)采場(chǎng)的穩(wěn)定性進(jìn)行分析。數(shù)值模擬結(jié)果表明:隨著采場(chǎng)開(kāi)拓延深,邊坡高度不斷加大,西采場(chǎng)北幫受斷層破碎帶影響,變形區(qū)域由局部變?yōu)檎w失穩(wěn)。