陳慶峰 黃文忠 程 斌 陳 濤 韓 流

（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)定性分析。

1 工程概況

秀山熔劑石灰?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ò)斷，使西礦段呈階梯狀，東礦段落呈齒狀。

2 斷層破碎帶邊坡變形機(jī)制分析

2.1 彎曲傾倒破壞機(jī)理分析

由于西采場(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 邊坡變形機(jī)理分析

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)生不利的影響。

3 邊坡數(shù)值模擬穩(wě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ā)生。

4 結(jié) 論

（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)。