張鯤華，陳玉明，張紹周

（昆明理工大學(xué)國土資源工程學(xué)院，云南 昆明 650093）

【其 他】

某露天礦排土場邊坡穩(wěn)定性分析

張鯤華，陳玉明，張紹周

（昆明理工大學(xué)國土資源工程學(xué)院，云南 昆明 650093）

排土場邊坡穩(wěn)定性問題成為阻礙我國露天礦發(fā)展的重大因素而得到國家和社會(huì)的廣泛重視。本文總結(jié)排土場的發(fā)展趨勢，結(jié)合某露天礦排土場的工勘地質(zhì)特征，采用極限平衡法和FLAC3D強(qiáng)度折減法對排土場穩(wěn)定性進(jìn)行計(jì)算，計(jì)算結(jié)果相差不超過0.04，計(jì)算結(jié)果顯示FLAC3D強(qiáng)度折減法計(jì)算安全系數(shù)具有很好的可靠性，計(jì)算不收斂判據(jù)與塑性貫通判據(jù)在本例中等效。安全系數(shù)平均值為0.480，反分析排土料的物理力學(xué)參數(shù)，認(rèn)為達(dá)不到穩(wěn)定要求的強(qiáng)度，故認(rèn)為該排土場存在重大安全隱患，應(yīng)及時(shí)治理。通過數(shù)值模擬研究邊坡的失穩(wěn)機(jī)制對于邊坡穩(wěn)定將有更為積極的意義。實(shí)踐過程中，應(yīng)當(dāng)重視計(jì)算參數(shù)獲取的準(zhǔn)確性。

排土場；邊坡穩(wěn)定性；極限平衡法；FLAC3D強(qiáng)度折減法

1 概述

露天開采在采掘工業(yè)中占有舉足輕重的地位。據(jù)估計(jì)，世界上礦產(chǎn)資源的露天開采比重已經(jīng)超過2/3。我國鐵礦石露天開采產(chǎn)量占總產(chǎn)量的86%～91%，有色金屬礦露天開采產(chǎn)量占總產(chǎn)量的46%～49%[1]。露天礦邊坡工程是一項(xiàng)關(guān)系礦山安全生產(chǎn)、礦山可持續(xù)發(fā)展和取得礦業(yè)開發(fā)重大盈利效益的關(guān)鍵課題[2]。隨著露天礦開采進(jìn)入中后期，露天礦排土場滑坡呈遞增趨勢，排土場滑坡不僅影響到礦山自身的安全生產(chǎn)，還會(huì)造成重大人員傷亡和環(huán)境破壞，并產(chǎn)生重大的經(jīng)濟(jì)損失和嚴(yán)重的社會(huì)影響[3]。

我國冶金礦山90%屬于露天開采，每年剝離的巖土量約4億t，占地10km2，其他有色金屬礦每年排棄的巖土約0.5億t，全國煤礦排棄的煤矸石每年約1億t，露天排土場占地面積約為全礦的30%～50%[4]。隨著邊坡穩(wěn)定性評價(jià)的發(fā)展，排土場的排土潛力得以較大的挖掘，以排土場增容研究為題的研究成果不斷涌現(xiàn)，改進(jìn)排土工藝為排土場的增容實(shí)現(xiàn)提供重要保證?；诠?jié)約用地和穩(wěn)定性的排土場優(yōu)化設(shè)計(jì)將成為排土場發(fā)展的必然趨勢。

2 某露天礦排土場工程概況

2.1 水文地質(zhì)

圖1 勘查區(qū)平面圖

工程區(qū)位于紅河水系左岸，由于區(qū)內(nèi)大部分為碳酸鹽巖分布，溶洞、溶隙發(fā)育，花崗巖出露區(qū)風(fēng)化強(qiáng)烈，節(jié)理裂隙，雖然降雨量充沛，但由于大氣降雨后多以垂直入滲的方式滲入地下，形成地下徑流，從ZK1～ZK6(圖1)反應(yīng)出勘察區(qū)淺層無地下水位，雨季僅有部分形成地表徑流，向南西匯入紅河。地下水類型主要以裸露型巖溶水為主，其次為碎屑巖風(fēng)化裂隙水、第四系為松散層孔隙水。本工程區(qū)基巖裸露，以大氣降水為地下水主要補(bǔ)給源，含水層大面積分布于該單元內(nèi)，巖溶裂隙發(fā)育，形成巖溶水徑流帶，地下水由東向西徑流。碎屑巖風(fēng)化裂隙水中的風(fēng)化裂隙水埋藏較淺，分布受地形、地貌控制，普遍為潛水性質(zhì)，風(fēng)化裂隙的發(fā)育程度與巖性、地形、地貌、植被等有密切關(guān)系，而且隨著深度的增加，風(fēng)化裂隙逐漸減弱，故一般淺，富集條件差。第四系松散層孔隙水含水層巖性單一，孔隙水分布零星，分散，富水性不均一，盆(谷)地內(nèi)地下水接受大氣降水、地表水的補(bǔ)給，由于地勢平坦，水力坡度小，富水性和滲透性較差，降雨后多形成地表徑流，向附近低凹處排泄，具有就地補(bǔ)給就地排泄的特點(diǎn)。

2.2 地形地質(zhì)

邊坡區(qū)屬構(gòu)造侵蝕中低山—低山區(qū)，為卡房斷裂中等切割河谷地貌，地勢陡峻，地形坡度26～50°，巖石風(fēng)化強(qiáng)烈，裂隙發(fā)育，雨水對邊坡不斷沖刷，降低了邊坡的安全等級，并且加速了邊坡的滑坡和崩塌速度。

根據(jù)地面工程地質(zhì)測繪和鉆探揭露，治理區(qū)主要出露地層為中生界三疊系灰?guī)r、印支期花崗巖、殘坡積層及人工堆積物。中生界三疊系灰?guī)r：上部厚度約4m為淺灰、青灰色，中—厚狀的強(qiáng)—弱風(fēng)化石灰?guī)r，溶洞和溶隙發(fā)育，是地表水的滲漏通道；4m以下為微—未風(fēng)化石灰?guī)r，青灰色，中—厚狀，結(jié)構(gòu)致密，脆性大，強(qiáng)度高。印支期花崗巖：強(qiáng)風(fēng)化花崗巖，灰白色，風(fēng)化裂隙發(fā)育，易碎；弱風(fēng)化花崗巖，灰白色，結(jié)構(gòu)致密，硬度高，強(qiáng)度大。殘坡積層：碎石土為淺黃、暗紅色塊石夾粉土，松散，塊石形態(tài)不規(guī)則，棱角明顯，混雜堆積，以石灰?guī)r、花崗巖強(qiáng)風(fēng)化碎塊為主，塊石粒徑100～300mm，粗礫含量40%～60%。人工堆積物(碎石土)：淺黃、暗紅色塊石夾粉土，松散，塊石形態(tài)不規(guī)則，棱角明顯，混雜堆積，以花崗巖強(qiáng)風(fēng)化碎塊為主。

2.3 物理力學(xué)性質(zhì)

采用不同深度鉆孔巖心做直剪試驗(yàn)測定排土的物理力學(xué)參數(shù)如表1。

表1 直剪試驗(yàn)測定排土料的物理力學(xué)參數(shù)

2.4 地震

工程區(qū)附近較為活躍的斷裂是小江斷裂，該斷裂屬殼斷裂，活動(dòng)時(shí)代Q3-4，活動(dòng)方式左行走滑，速率1.0～8.6mm/a，北強(qiáng)南弱，發(fā)生5.0～5.9級地震18次，6.0～6.9級地震9次，≥7.0級地震4次，最大地震8級。

3 穩(wěn)定性分析

3.1 方法選取

(1) 極限平衡法。

極限平衡法是邊坡穩(wěn)定性分析計(jì)算的主要方法，也是工程實(shí)踐中運(yùn)用最多的一種方法。極限平衡法是根據(jù)邊坡上的滑體或滑體分塊的靜力平衡原理分析邊坡各種破壞模式下的受力狀態(tài)，以邊坡滑體上的抗滑力和下滑力之間的比值來評價(jià)邊坡的穩(wěn)定性。在大多數(shù)情況下，邊坡穩(wěn)定性問題是靜不定的，邊坡受力的平衡方程中未知數(shù)的個(gè)數(shù)大于方程個(gè)數(shù)，為了使方程有唯一解，通常引入假設(shè)將靜不定轉(zhuǎn)化為靜可解，這種處理使解決問題的嚴(yán)密性受到破壞。但由于極限平衡法計(jì)算模型簡單、計(jì)算參數(shù)量化準(zhǔn)確、計(jì)算結(jié)果直接適用而被廣泛應(yīng)用。極限平衡法理論體系的發(fā)展完善形成了以不同適用范圍和假設(shè)條件的計(jì)算方法，李小春等[5]認(rèn)為不應(yīng)指定某一種方法作為驗(yàn)算的唯一方法，可以選擇多種精度較高的方法(如M-P、Spencer、Bishop、L-K法)的均值作為驗(yàn)算方法。當(dāng)條件具備時(shí)，應(yīng)盡量采用滿足所有平衡條件的方法，其具有較高的精度與較低的標(biāo)準(zhǔn)差，計(jì)算結(jié)果穩(wěn)定可靠。Bishop法具有良好的替代性，計(jì)算精度較高。

(2) FLAC3D強(qiáng)度折減方法。

FLAC3D中采用強(qiáng)度折減系數(shù)來定義安全系數(shù)[6-7]。若邊坡采用Mohr-Coulomb準(zhǔn)則，影響其穩(wěn)定性的強(qiáng)度參數(shù)粘結(jié)力C 和內(nèi)摩擦角φ，將坡體原始粘結(jié)力C0

和內(nèi)摩擦角φ0同時(shí)除以一個(gè)折減系數(shù)K，然后進(jìn)行數(shù)值分析。通過不斷增大K，反復(fù)分析直至邊坡達(dá)到臨界破壞狀態(tài)。假設(shè)此時(shí)的粘結(jié)力和內(nèi)摩擦角為Ccr和φcr，由于邊坡處于臨界狀態(tài)，所對應(yīng)的安全系數(shù)Kcr=1，可得原始邊坡所對應(yīng)的安全系數(shù)為：

本文采用極限平衡法與有限差分方法相結(jié)合計(jì)算排土場邊坡安全系數(shù)，對排土場邊坡的穩(wěn)定性作出合理評價(jià)。在有限差分方法中采用計(jì)算不收斂判據(jù)，破壞之前計(jì)算收斂，破壞之后滑面上巖土體無限流動(dòng)計(jì)算不收斂。計(jì)算過程中記錄最后一個(gè)循環(huán)最大不平衡力及最后一步塑性貫通狀態(tài)。

3.2 計(jì)算過程

(1) 模型建立。

排土料在山谷中部堆置較厚，山谷兩側(cè)堆置變薄。由此認(rèn)為沿原地形山谷線剖面為排土場存在安全隱患最突出的剖面，以此剖面作為排土場安全系數(shù)計(jì)算的剖面(如圖2)，圖2中人工堆積物水平總長度約310m，沿坡向下游堆積，厚度逐漸變厚，坡頂ZK4位置堆積存5.4m，中部ZK3位置堆積厚33.6m，下游邊緣距ZK2約50m，堆積物結(jié)構(gòu)松散，易產(chǎn)生滑坡和泥石流，從坡腳向坡頂堆積物粒徑逐漸變小，分選明顯，其下覆地層為強(qiáng)風(fēng)化易破碎結(jié)構(gòu)松散花崗巖和強(qiáng)—弱風(fēng)化石灰?guī)r。

圖2 排土場最危險(xiǎn)剖面

由于殘、坡積層厚度和范圍較小，為了計(jì)算簡化，將其物理力學(xué)參數(shù)用印支期花崗巖的代替。計(jì)算模型既要包括整個(gè)潛在滑坡體，又要有一定的外延，盡可能減少邊界條件對計(jì)算結(jié)果的影響。當(dāng)坡腳到左端邊界的距離為坡高的1.5倍，坡頂?shù)接叶诉吔绲木嚯x為坡高的2.5倍，且上下邊界總高不低于2倍坡高時(shí)，計(jì)算精度最為理想[8]。建立如圖3的計(jì)算模型。

圖3 計(jì)算模型

(2) 計(jì)算結(jié)果。

將圖3中簡化的計(jì)算模型導(dǎo)入到理正軟件中，應(yīng)用復(fù)雜邊坡計(jì)算模塊，輸入物理力學(xué)參數(shù)和控制參數(shù)計(jì)算排土場邊坡的安全系數(shù)。假設(shè)滑坡為圓弧形，采用通用方法自動(dòng)搜索最危險(xiǎn)滑裂面，土條寬0.5m，搜索時(shí)的圓心步長和半徑步長都為5m。土柱位置誤差為1m。

圖4 FLAC3D計(jì)算模型

在FLAC3D軟件中建立如圖4的計(jì)算模型，自編強(qiáng)度折減程序計(jì)算排土場邊坡的安全系數(shù)。折減過程只對排土料進(jìn)行。模型建立了45 312個(gè)單元共57 970個(gè)節(jié)點(diǎn)的有限差分網(wǎng)格，采用Mohr-Coulomb模型并賦相應(yīng)材料參數(shù)，模型受重力影響，不考慮孔隙水壓，模型四周和底面固定。通過試算，模型在計(jì)算30 000步以內(nèi)已經(jīng)達(dá)到平衡狀態(tài)，同時(shí)將計(jì)算最大不平衡力與典型內(nèi)力的比率不超過1.0e-5為收斂準(zhǔn)則，計(jì)算精度設(shè)置為0.001，計(jì)算結(jié)果如圖5、圖6和圖7所示。

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

圖5 最大不平衡力

表2 安全系數(shù)計(jì)算結(jié)果

安全系數(shù)計(jì)算結(jié)果見表2，由表2可以看出，排土場邊坡安全系數(shù)都小于1，排土場邊坡已經(jīng)破壞，實(shí)地踏勘排土場邊坡現(xiàn)狀并未出現(xiàn)坡體拉伸破壞?？赡艽嬖谌缦略颍孩偃訋r心的級配組成與排土場的實(shí)際的級配組成不一致。取樣巖心來自ZK3和ZK4，取了6組樣進(jìn)行直剪試驗(yàn)；②直剪試驗(yàn)存在缺陷，剪切面限定在上下盒之間的平面，而不是沿土樣最薄弱的面剪切破壞。剪切面上的剪應(yīng)力分布不均勻，且豎向荷載會(huì)發(fā)生偏轉(zhuǎn)，主應(yīng)力的大小和方向都是變化的，對試驗(yàn)結(jié)果產(chǎn)生影響；③粘性土具有觸變性，其強(qiáng)度會(huì)因?yàn)槭艿綌_動(dòng)而削弱。在運(yùn)輸過程中，巖土體的結(jié)構(gòu)可能因?yàn)檎饎?dòng)發(fā)生改變。

圖6 塑性貫通區(qū)

圖7 X 方向的位移云圖

雖然巖土體的物理力學(xué)參數(shù)存在一定的誤差，但是采用不同計(jì)算方法計(jì)算的安全系數(shù)之間差別不大，最大值與最小值之間相差為0.04，由此證明計(jì)算的正確，該排土場邊坡穩(wěn)定性系數(shù)為0.480。排土場邊坡計(jì)算過程未考慮孔隙水壓及地震的影響，孔隙水壓和地震因素都是降低邊坡穩(wěn)定性的因素，即使試驗(yàn)得到的巖土體物理力學(xué)參數(shù)偏小，不考慮孔隙水壓和地震將使計(jì)算結(jié)果誤差減小。另外，當(dāng)排土場邊坡達(dá)到臨界狀態(tài)時(shí)，需要材料的粘結(jié)力為40.899kPa和內(nèi)摩擦角為27.937°，礫砂、角礫和粉質(zhì)粘土是在天然狀態(tài)下難以達(dá)到此強(qiáng)度[9]。所以排土場邊坡存在重大安全隱患，應(yīng)當(dāng)及時(shí)采取治理滑坡措施。

4 結(jié)論

工程案例安全系數(shù)計(jì)算說明FLAC3D強(qiáng)度折減系數(shù)法具有很好的可靠性。FLAC3D強(qiáng)度折減不需要假設(shè)危險(xiǎn)滑移面形狀，采用計(jì)算不收斂作為邊坡破壞判據(jù)是可行的。邊坡滑動(dòng)失穩(wěn)是在內(nèi)外綜合因素的作用下，系統(tǒng)從量變到質(zhì)變，從劣化因素的不斷積累到最終發(fā)生突變失穩(wěn)的一個(gè)變化過程[10]，通過數(shù)值模擬研究邊坡的失穩(wěn)機(jī)制對于邊坡穩(wěn)定將有更為積極的意義。考察FLAC3D計(jì)算的最大不平衡力(圖5)、塑性區(qū)貫通狀況(圖6)和X方向的位移云圖(圖7)，圖5顯示最大不平衡力是收斂的，模型在30 000步內(nèi)能達(dá)到平衡狀態(tài)。圖6顯示在計(jì)算不收斂判據(jù)成立條件下塑性區(qū)貫通仍然成立，說明這兩個(gè)判據(jù)在判斷邊坡破壞時(shí)具有相同效應(yīng)。圖7顯示邊坡的破壞為圓弧形(X方向位移為0)，由此證明在極限平衡法中的假設(shè)成立。

在計(jì)算過程中，參數(shù)的準(zhǔn)確與否將決定計(jì)算結(jié)果的可靠性。在生產(chǎn)實(shí)踐中，應(yīng)當(dāng)重視巖土體物理力學(xué)參數(shù)的獲取，采取有效措施減小參數(shù)誤差。

[1]梁克鈞.國內(nèi)外露天采礦技術(shù)水平與動(dòng)向[J].河北礦冶學(xué)院學(xué)報(bào)，1983(1)：45-59.

[2]張四維.我國露天礦邊坡工程研究成果與進(jìn)展[J].地質(zhì)災(zāi)害與環(huán)境保護(hù)，2000，11(2)：98-101.

[3]孫世國，楊宏.典型排土場邊坡穩(wěn)定性控制技術(shù)[M].北京：冶金工業(yè)出版社，2011.

[4]薄建芬，候豆豆.露天礦山排土場安全管理措施研究與應(yīng)用[C]∥魯冀晉瓊粵川遼七省金屬(冶金)學(xué)會(huì)第十九屆礦山學(xué)術(shù)交流會(huì)，2009：982-988.

[5]李小春，任偉，王少泉，等.論金屬礦山排土場設(shè)計(jì)規(guī)范中邊坡極限平衡計(jì)算方法的選取[J].巖石力學(xué)與工程學(xué)報(bào)，2011，30(Z2)：4136-4142.

[6]趙尚毅，鄭穎人，時(shí)衛(wèi)民，等.用有限元強(qiáng)度折減法求邊坡穩(wěn)定安全系數(shù)[J].巖土工程學(xué)報(bào)，2002，24(3)：343-346.

[7]張魯渝，鄭穎人，趙尚毅，等.有限元強(qiáng)度折減系數(shù)法計(jì)算土坡穩(wěn)定安全系數(shù)的精度研究[J].水利學(xué)報(bào)，2003，34(1)：21-27.

[8]鄭穎人，趙尚毅.有限元強(qiáng)度折減法在土坡與巖坡中的應(yīng)用[J].巖石力學(xué)與工程學(xué)報(bào)，2004，23(19)：3381-3388.

[9]林宗元.巖土工程勘察設(shè)計(jì)手冊[M].沈陽：遼寧科學(xué)技術(shù)出版社，1996.

[10]孫書偉，林杭，任連偉.FLAC3D在巖土工程中的應(yīng)用[M].北京：中國水利水電出版社，2011.

TD824.71；TD824.8

A

1007-9386(2014)05-0059-04

2014-02-28