李全明 張 紅 董紹艷

(1.中國安全生產(chǎn)科學(xué)研究院：2.華北科技學(xué)院)

?

·安全·環(huán)?！?/p>

尾礦壩壩坡穩(wěn)定分析的新思路及實例驗證*

李全明1,2張 紅2董紹艷1

(1.中國安全生產(chǎn)科學(xué)研究院：2.華北科技學(xué)院)

尾礦庫是具有高勢能的重大危險源，一旦發(fā)生事故不僅造成財產(chǎn)損失，還會導(dǎo)致周邊環(huán)境污染，尤其是“三邊庫”和“頭頂庫”，更是時刻威脅著人民群眾的生命財產(chǎn)安全?，F(xiàn)行法規(guī)標(biāo)準(zhǔn)均采用基于極限平衡分析的邊坡穩(wěn)定方法，雖然方便操作，但難以準(zhǔn)確把握尾礦壩壩坡的實際狀況。為此提出在顆粒流理論基礎(chǔ)上的強(qiáng)度折減法，定量計算尾礦壩壩坡的穩(wěn)定性，模擬潰壩過程中壩體內(nèi)部微觀粒子的運動特性及潰壩演化規(guī)律，克服了極限平衡法各種假設(shè)的缺點。通過某一尾礦壩壩坡實例，計算了尾礦壩壩坡強(qiáng)度折減過程中微觀粒子的相互作用過程，動態(tài)反映壩坡失穩(wěn)狀態(tài)，驗證了該方法的適用性。

顆粒流 強(qiáng)度折減法 尾礦壩 安全系數(shù)

尾礦庫是指筑壩攔截谷口或圍地構(gòu)成的用以貯存金屬非金屬礦山礦石選別后排除尾礦的場所[1]。我國自2007年開始逐步加強(qiáng)含赤泥庫在內(nèi)的非煤礦山尾礦庫隱患清查及綜合治理工作，截止到2014年底，我國共有尾礦庫11 359座，其中危庫、險庫數(shù)量達(dá)到31座，病庫達(dá)到772座[2]。尾礦庫是具有高勢能的重大危險源，一旦發(fā)生事故，不僅危及到人民群眾的生存安全，更會對環(huán)境造成無法挽回的破壞。而赤泥庫中強(qiáng)堿性的化學(xué)物質(zhì)對生物和金屬、硅質(zhì)材料具有強(qiáng)腐蝕性，一旦泄露會造成土地堿化及嚴(yán)重的水污染。2010年10月4日，匈牙利一鋁礦發(fā)生潰壩事故，約70萬m3有毒赤泥外泄，流入附近河流引發(fā)多國恐慌，該事故共造成至少10人死亡，120人受傷，災(zāi)后治理18個月。

據(jù)統(tǒng)計，尾礦庫事故中潰壩事故約占事故總數(shù)的60%以上[3],因此，對尾礦壩壩體進(jìn)行穩(wěn)定性分析，是保證尾礦庫安全的重點。目前常用的尾礦庫邊坡穩(wěn)定性分析方法有：有限元法[4]，極限平衡法[5]及快速拉格朗日法[6]，從宏觀角度分析尾礦壩的穩(wěn)定性，將壩體模型進(jìn)行細(xì)部劃分，建立穩(wěn)定性計算方程，再通過強(qiáng)度折減計算壩體安全系數(shù)。但由于對散體材料的變形及受力分析有一定的局限性，無法模擬壩體潰壩的過程及微觀粒子之間的運動關(guān)系。1979年，cundall和strack[7]提出了適用于土力學(xué)的離散元法(DEM)，在隨后的幾年中，在巖土工程、水電工程和采礦工程中應(yīng)用廣泛，該方法為分析尾礦庫壩體穩(wěn)定性提供了新思路。離散元法分析尾礦庫壩體穩(wěn)定性的優(yōu)勢在于：不用假定滑動面，處理非連續(xù)介質(zhì)問題中不受變形量的限制，從微觀上模擬壩體粒子間的相對運動、受力變化及坡體滑移趨勢，有效的模擬壩體開裂及分離現(xiàn)象。本文以某尾礦庫堆積壩為研究對象，在PFC2D顆粒流基礎(chǔ)上，結(jié)合強(qiáng)度折減法研究壩體失穩(wěn)趨勢及受力變化，計算尾礦庫安全系數(shù)。

1 基本原理

1.1 顆粒流理論原理

PFC2D(Particle Follow Code 2 Dimension)是通過離散元法(DEM)模擬圓形顆粒運動及其相互作用的軟件[8]，其主要特點為：①將分析對象分解為微觀顆粒集合體，從微觀角度研究其力學(xué)行為；② 顆粒運動遵循力-位移定律和牛頓第二定律，在迭代計算過程中兩種定律交互作用，直至模型達(dá)到平衡或產(chǎn)生破壞無法保持穩(wěn)定時止[9]。PFC中主要有以下幾個假設(shè)：①模型中所有顆粒為剛性的圓盤或球體；②顆粒間的接觸近似“點接觸”；③顆粒間接觸行為具有柔性特性，并允許顆粒間存在一定重疊；④顆粒間的接觸部位可以建立黏性特性；⑤接觸力通過力-位移定律與顆粒間重疊量建立聯(lián)系。

松散材料堆積體可以看做是無數(shù)微觀圓形顆粒組成的集合體，通過彼此間的約束力保持穩(wěn)定。PFC2D針對此種情況提供了兩種黏結(jié)模型模擬其間的接觸行為，即：接觸黏結(jié)和平行黏結(jié)。其中接觸黏結(jié)如圖1(a)所示，顆粒間僅在接觸點周圍極小范圍內(nèi)產(chǎn)生黏結(jié)作用，接觸鍵僅承受拉力和剪力作用，不能承受彎矩作用力。平行黏結(jié)如圖1(b)所示，顆粒間在接觸點一定范圍內(nèi)存在一個圓形或矩形的膠結(jié)物，這種材料特性的有效剛度與接觸點剛度平行，可以承受拉力、剪力和彎矩三者共同的作用力。平行黏結(jié)受力性質(zhì)與壩體顆粒間的受力形式一致，故本次模型中的顆粒選用平行黏結(jié)。

圖1 接觸鍵連接模型

1.2 強(qiáng)度折減方法

1975年，Zienkiewicz、Humpheson C[10]等人提出強(qiáng)度折減的概念，基本原理[11]就是將材料的抗剪強(qiáng)度參數(shù)c、tanφ，同時除以系數(shù)F(F>1)，得到一個全新的c′、 tanφ′,利用新的強(qiáng)度參數(shù)測試材料的穩(wěn)定性并不斷進(jìn)行調(diào)整，直到材料達(dá)到臨界破壞，此時的折減系數(shù)就是材料的安全系數(shù)。具體方程為：

(1)

(2)

1996年，Duncan[12]提出安全系數(shù)定義，為土體的實際剪切強(qiáng)度與臨界破壞剪切強(qiáng)度的比值，同時也可以解釋為土壤剪切強(qiáng)度經(jīng)折減后達(dá)到臨界破壞時的系數(shù)值。

2 模型建立及分析

2.1 尾礦壩實例及參數(shù)選取

PFC主要模擬微觀粒子之間的運動關(guān)系，其微觀參數(shù)值主要為顆粒法向和切向剛度、顆粒間平行黏結(jié)剛度、顆粒間平行黏結(jié)強(qiáng)度等。宏觀參數(shù)如泊松比、彈性模量等數(shù)值不能直接用于該軟件，需要根據(jù)軟件進(jìn)行一系列的模擬,找到與之相對應(yīng)的微觀粒子參數(shù)值。為了弄清楚兩者之間的關(guān)系，眾多學(xué)者進(jìn)行了大量研究[13-15]，發(fā)現(xiàn)宏觀參數(shù)與微觀參數(shù)之間存在一定的線性相關(guān)性，但尚未發(fā)現(xiàn)彼此之間明確的數(shù)量級關(guān)系。目前獲取粒子微觀參數(shù)的通用方法是利用PFC模擬雙軸壓縮實驗和巴西劈裂實驗，通過不斷調(diào)整各個微觀參數(shù)值，使其模擬結(jié)果與實驗結(jié)果相一致，此時各個參數(shù)值就是模擬對象的微觀參數(shù)值。本文中同樣利用PFC2D模擬雙軸試驗進(jìn)行微觀參數(shù)選取。

以周喻[16]等的尾礦庫數(shù)據(jù)資料進(jìn)行模擬計算，尾礦庫長225 m，高90 m，壩坡角26°，壩體安全平臺寬7.5 m，堆積壩下方土層分別為黏土層和巖石層，其宏觀參數(shù)如表1所示。利用PFC2D進(jìn)行雙軸試驗和巴西劈裂實驗計算出其微觀參數(shù)，如表2所示。

表1 尾礦庫宏觀參數(shù)

表2 尾礦庫微觀參數(shù)

2.2 模型建立

嚴(yán)格按照工程實際尺寸計算PFC模型尺寸，尾礦庫簡圖如圖2所示。

圖2 尾礦庫簡圖(單位：m)

巖層邊界及模型周邊約束均按照墻體繪制，墻體剛度設(shè)置為1×1010N/m，根據(jù)每個巖層面積及孔隙率要求生成相應(yīng)數(shù)量的粒子，粒徑0.75～1.25 m隨機(jī)分布，經(jīng)初步迭代達(dá)到平衡后刪除接觸少于3個的浮動粒子及巖層間的墻體，并將微觀參數(shù)賦予各個巖層，在自重作用下再次迭代達(dá)到平衡，最終生成粒子數(shù)4 673個。迭代過程中模型兩側(cè)約束X方向位移，底部約束X、Y方向位移。堆積壩最終效果如圖3所示。

圖3 堆積壩數(shù)值模擬效果

2.3 強(qiáng)度折減計算安全系數(shù)

堆積壩變形穩(wěn)定后，依照強(qiáng)度折減法計算原理，將堆積壩抗剪強(qiáng)度參數(shù)及內(nèi)聚力均按照同比例折減，即在PFC中將影響強(qiáng)度的平行黏結(jié)強(qiáng)度、黏結(jié)半徑及摩擦系數(shù)按照同比例折減。每次折減計算時步為20 000，折減系數(shù)以1.0開始，以0.1遞增，為了避免折減系數(shù)遞增過程中造成堆積壩累積應(yīng)變而影響分析結(jié)果，每次折減計算均取壩體穩(wěn)定后的模型為基礎(chǔ)進(jìn)行平行計算。

取壩體底端、中端和頂端及黏土層邊緣4個位置布置監(jiān)測點，監(jiān)測并記錄壩體在各折減系數(shù)下總位移及監(jiān)測點位移變化量。監(jiān)測點位置如圖4所示，壩體總位移及監(jiān)測點位移變化如圖5、圖6所示。從圖5中可以看出，折減系數(shù)為1.0～1.7的壩體位移基本無變化，折減系數(shù)為1.8～2.3壩體出現(xiàn)較小位移，但并未對其穩(wěn)定性造成影響；折減系數(shù)自2.4以后，壩體位移值增大，隨著強(qiáng)度進(jìn)一步減弱，整體邊坡失穩(wěn)，出現(xiàn)大變形。由此可以判斷，折減系數(shù)為2.3時，壩體位于邊坡失穩(wěn)臨界值，故其安全系數(shù)為2.3。

圖4 監(jiān)測點位置

圖5 壩體總位移變化量

監(jiān)測點位移變化中出現(xiàn)X方向位移遠(yuǎn)遠(yuǎn)大于Y方向位移，與實際壩體初始滑坡時巖體Y方向位移大于X方向位移存在差異，其原因是壩體堆積角度不同，在垮塌過程中粒子在X、Y方向位移變化量不同，本文堆積壩坡度為26°，理論上壩體在滑坡運動過程中，顆粒水平方向每運動1 m，垂直方向降低0.48 m，但由于具體運動過程中粒子間的相互作用，使顆粒兩個方向的位移值不是絕對的線性相關(guān)，與整體滑坡運動過程中出現(xiàn)壩體水平位移大于垂直位移且位移值無明顯線性相關(guān)的情況是合理的。

圖6 監(jiān)測點位移變化量

2.4 壩體失穩(wěn)分析

為進(jìn)一步了解壩體失穩(wěn)過程，分別取折減系數(shù)2.4、2.9、3.2時的位移矢量圖、接觸力圖進(jìn)行潰壩分析。位移矢量圖如圖7所示。

圖7 不同折減系數(shù)位移矢量圖

圖7中箭頭表示顆粒位移方向，長度表示位移量，長度越大顆粒發(fā)生的位移越大。堆積壩變形穩(wěn)定后，壩坡在重力的作用下具有向下運動的趨勢。隨著強(qiáng)度值開始逐步減小，壩體對顆粒運動的約束力越來越小，顆粒開始出現(xiàn)緩慢下滑且滑動越來越明顯，并在壩體中顯現(xiàn)出部分裂隙；強(qiáng)度值進(jìn)一步減小，裂隙越來越大逐漸形成貫通裂隙，坡體發(fā)生破壞，此時顆粒間的位移出現(xiàn)較大的變化；強(qiáng)度再次降低后，坡體破壞程度加大，顆粒之間的約束力急劇減小，壩體最終垮塌形成滑坡。從圖中可以看出，坡體滑坡形成后，影響長度非常長，經(jīng)計算，滑坡位移約58 m。

顆粒間接觸力用線條表示，線條顏色越深代表顆粒間的接觸力越大，如圖8所示。根據(jù)接觸力顯示，壩體初期平衡階段由于顆粒受到重力的作用，因此在堆積壩的底部所受接觸力最大，壩體強(qiáng)度值減小后，壩體開始出現(xiàn)下滑趨勢，邊坡底部顆粒最先出現(xiàn)滑移趨勢，其拉應(yīng)力逐步增大，隨著整個坡體開始出現(xiàn)向下滑移，坡面應(yīng)力開始集中，接觸力逐漸增大，向下滑移過程中，坡頂?shù)牧Ｗ訑D壓底部的粒子促使其滑移，坡頂范圍內(nèi)壓應(yīng)力逐漸增大，為了阻止坡面運動，壩體底部的作用力逐漸由壓力轉(zhuǎn)化為拉力并逐漸增大；隨著剪力進(jìn)一步降低，顆粒間的接觸鍵開始破壞，壩體出現(xiàn)裂隙，此時坡面開始逐步向下滑移，坡面以拉力為主；隨著顆?；萍敖佑|鍵破壞量增多，堆積壩底部重力作用減小拉力逐漸降低，直至壩體失穩(wěn)，此時壩體在滑移過程中，顆粒間的接觸鍵幾乎完全破壞，壩底顆粒作用力接近于零，不再對上體顆粒有約束作用，整個壩體完全破壞。

3 結(jié) 論

(1)通過對顆粒流程序運用強(qiáng)度折減原理，計算出的安全系數(shù)為2.3。

(2)利用PFC與強(qiáng)度折減法結(jié)合分析尾礦壩穩(wěn)定性方式，為尾礦庫的穩(wěn)定性判斷提供了量化方法。該方法可有效地模擬包括赤泥庫在內(nèi)的尾礦庫潰壩過程，可以作為衡量尾礦庫壩體穩(wěn)定性的方法之一，同時也證明了PFC軟件在尾礦庫穩(wěn)定性評價中的可行性。

(3)通過進(jìn)行潰壩模擬，可以直觀的看到尾礦庫潰壩的影響范圍，進(jìn)而可以確定潰壩安全邊界，可以在邊界周圍設(shè)置圍擋或警戒標(biāo)志，減小潰壩可能造成的損失。

(4)算例顯示實際尾礦庫潰壩的滑裂面是從軟弱面發(fā)展演化而成的，并不一定遵守圓弧滑動的形態(tài)。

圖8 不同折減系數(shù)接觸力

[1] 謝旭陽,王云海,張興凱,等.我國尾礦庫數(shù)據(jù)庫的建立[J].中國安全生產(chǎn)科學(xué)技術(shù),2008,4(2):53-56.

[2] 國家安全監(jiān)管總局等七部門.關(guān)于印發(fā)深入開展尾礦庫綜合治理行動方案的通知[2013]58號[Z].北京：國家安全監(jiān)管總局，2013.

[3] 楊麗紅,李全明,程五一,等.國內(nèi)外尾礦壩事故主要危險因素的分析研究[J].中國安全生產(chǎn)科學(xué)技術(shù),2008,4(5):28-31.

[4] 連鎮(zhèn)營,韓國城,孔憲京.強(qiáng)度折減有限元法研究開挖邊坡的穩(wěn)定性[J].巖土工程學(xué)報,2001,23(4):407-411.

[5] 鄒廣電,魏汝龍.土坡穩(wěn)定分析普遍極限平衡法數(shù)值解的理論及方法研究[J].巖石力學(xué)與工程學(xué)報,2006,25(2):363-370.

[6] 鄒 棟,鄭 宏.快速拉格朗日法及其在邊坡穩(wěn)定性分析中的應(yīng)用[J].礦業(yè)研究與開發(fā),2005,25(5):80-83.

[7] Cundall P A,Strack O D L.A discrete numerical model for granular assemblies[J].Geotechnique,1979,29(1):47-65.

[8] Itasca Consulting Group.PFC2D(Particle Flow Code in 2 Dimensions)User's Guide[R].Minnesota:Itasca Consulting Group,2002.

[9] 石 崇,徐衛(wèi)亞.顆粒流(PFC2D/3D)數(shù)值模擬技巧與實踐[M].北京:中國建筑工業(yè)出版社,2015.

[10] Zienkiewicz O C,Humpheson C,Lewis R W.Associated and non-associated visco-plasticity and plasticity in soil mechanics[J].Geotechnique,1975,25(4):671-689.

[11] 周先奇,徐衛(wèi)亞,鄭文棠,等.基于強(qiáng)度折減的離散單元法在巖質(zhì)邊坡未定型分析中的應(yīng)用[J].長江科學(xué)院院報,2008,25(5):107-114.

[12] Duncan J M.State of the art:Limit equilibrium and finite-element analysis of slopes[J].Journal of Geotechnical Engineering,AS-CE,1996,122(7):577-596.

[13] 彭 博,陳玉明,袁利偉.基于PFC2D馬鞍山鐵礦2號露天采場邊坡穩(wěn)定性分析[J].黃金,2014,35(12):37-40.

[14] 周 健,王家全,曾遠(yuǎn),等.顆粒流強(qiáng)度折減法和重力增加法的邊坡安全系數(shù)研究[J].巖土力學(xué),2009,30(6):1549-1554.

[15] 楊嘯銘.基于PFC2D的邱山鐵礦尾礦庫壩體穩(wěn)定性和潰壩模擬分析研究[D].昆明:昆明理工大學(xué),2015.

[16] 周 喻,吳順川,馬 聰,等.基于顆粒流理論的露天礦排土場穩(wěn)定性分析[J].中國礦業(yè),2010,19(7):94-97,101.

New Ideas and Case Analysis of Slope Stability of Tailings Dam

Li Quanming1,2Zhang Hong2Dong Shaoyan1

(1.China Academy of Safety and Science & Technology;2.North China Institute of Science & Technology)

Tailings ponds is a major danger source of high potential energy,if the accidents is happened,not only caused the damage to property,but also lead to the pollution of the surrounding environment,especially the the safety of the lives and property of the masses of the people is threatened by the"trilateral library"and"the head of the library"all the time.The current slope stability analysis method is the limit equilibrium analysis method in some relative codes and standards.The method is easy to operate,but it is difficult to grasp the actual situation of the tailings dam slope accurately.Based on the particle flow theory,the strength reduction method is proposed to analyze the slope stability of tailings dam quantitatively.The movement characteristics of the micro particle inside the dam body and dam break evolution regularity can be simulated effectively by the method,therefore,the shortcomings of the assumptions of limit equilibrium method is overcomed.Taking a tailings dam as the study example,the interaction process of the micro particle in the strength reduction process of the tailings dam slope is calculated,which can reflect the dynamic process of instability of tailings dam slope.The method proposed in this paper is verified by the method proposed in this paper.

Particle flow,Strength reduction method,Tailings dams,Safety coefficient

*國家自然科學(xué)基金項目(編號：71373245)。

2016-09-07)

李全明(1979—)，男，教授級高級工程師，博士，100012 北京市朝陽區(qū)北苑路32號甲1號樓。