劉余欣，易飛，徐裕華，樊莉

（1.南昌大學第一附屬醫(yī)院，江西 南昌 330006；2.江西工業(yè)貿(mào)易職業(yè)技術(shù)學院，江西 南昌 330038）

1 引言

尾礦壩作為一種礦石尾料堆積形成的工業(yè)構(gòu)筑物，一旦發(fā)生潰壩破壞力極強。相關(guān)統(tǒng)計結(jié)果表明[1]，尾礦壩潰壩造成的損失在世界范圍內(nèi)的各類災(zāi)害中排名第18位。例如，發(fā)生在1997年圭亞那的尾礦壩潰壩事件造成了20多人死亡的嚴重后果[2]；發(fā)生在2008年的陜西臨汾尾礦壩潰壩事故累積造成128人死亡，35人受傷[3]。由此可見，尾礦壩潰壩造成的危害性極大。以往的研究表明[4-5]，尾礦壩的失穩(wěn)潰壩往往發(fā)生于連續(xù)降雨期間，連續(xù)降雨會增加壩體的飽和度，降低土體的基質(zhì)吸力，使得土體抗剪強度下降，最終導(dǎo)致尾礦壩失穩(wěn)破壞。鑒于上述原因，對連續(xù)降雨期間尾礦壩失穩(wěn)的機制進行分析并采取有針對性的措施進行預(yù)防是極為有必要的。

本文從某尾礦壩實際工程案例出發(fā)，利用GeoStudio軟件分析了連續(xù)降雨條件下該尾礦壩總水頭、浸潤線以及壩坡穩(wěn)定性的變化情況，為降雨入滲條件下的尾礦壩滲流及穩(wěn)定性分析提供了參考，結(jié)論可供類似工程參考應(yīng)用。

2 滲流及穩(wěn)定性計算原理

2.1 滲流計算原理

尾礦壩的滲流問題屬于連續(xù)介質(zhì)的飽和-非飽和非恒定滲流問題，連續(xù)介質(zhì)飽和-非飽和非恒定滲流理論問題的描述及其求解概要敘述如下。

2.1.1 滲流控制方程

假設(shè)水在非飽和土體中仍舊能滿足達西定律，水流在滿足連續(xù)性原理的同時孔隙氣壓力不隨時間變化，在忽略不同流體流動同土結(jié)構(gòu)平衡條件之間的相互作用條件下，飽和-非飽和區(qū)土地地下水的非穩(wěn)定滲流控制方程可以寫成：

當（1）式中總水頭H恒定不變時，地下水穩(wěn)態(tài)流控制方程可以寫成：

上式中，H為總水頭，h為壓力水頭，y為位置水頭，總水頭為壓力水頭和位置水頭之和，即H=h+y。當處于飽和區(qū)時，滲透系數(shù)為飽和值，與h大小無關(guān)；當處于非飽和區(qū)時，滲透系數(shù)是壓力水頭的函數(shù)；kx、ky分別代表x與y方向的滲透系數(shù)；Q為微元體邊界流量；mw代表體積含水量變化系數(shù)，其值為土水特征曲線的斜率；γw為水的容重。

2.1.2 邊界條件

在水位變動的條件下，尾礦壩滲流問題的邊界條件可以寫成：

上式中，H（x,y）表示在二維條件下尾礦壩的某一邊界S在t時刻的水頭，根據(jù)尾礦壩的庫水位變動條件，可確定邊界S的水頭邊界值H1(x,y,t)。

2.1.3 有限元計算方法

由式（1）和（3）構(gòu)成滲流問題的定解條件，由該定解條件，通過適當?shù)臄?shù)學方法可求出尾礦壩滿足邊界條件（3）的唯一水頭分布的解。

2.2 穩(wěn)定性計算原理

本次邊坡穩(wěn)定性計算采用的方法為剛體極限平衡法，以簡化Bishop法為主。為了便于對結(jié)果對比分析，同時采用瑞典圓弧法以及摩根斯坦-普萊斯法對邊坡穩(wěn)定性進行計算。簡化Bishop法的基本原理為條分法，其假設(shè)條塊間的作用力方向只為水平方向，以此建立整體力矩平衡方程并進行求解。簡化Bishop法在忽略條塊間的豎向剪力并簡化安全系數(shù)的求解過程同時，精度卻幾乎沒有降低，在工程中被廣泛運用，該法已被寫入各國的相關(guān)規(guī)程中。

3 工程實例

3.1 數(shù)值模型的建立

本次研究的尾礦壩壩體材料主要包括尾細砂、尾粉砂、粉質(zhì)黏土、填土以及塊石棱體等，根據(jù)其典型剖面建立有限元網(wǎng)格模型，見圖1，該網(wǎng)格模型共有7486個節(jié)點，7293個單元。

圖1 計算模型網(wǎng)格圖

3.2 材料的參數(shù)

通過一系列室內(nèi)土工試驗并參考相關(guān)勘察報告，最終確定了尾礦壩壩體中各個材料計算所需的力學參數(shù)數(shù)值，如表1所示。

3.3 計算工況

本次計算選取的計算工況為尾礦壩在上游正常水位運行條件的前提下，同時以50mm/d的降雨強度（暴雨）連續(xù)降水3天，計算此時尾礦壩的滲流及壩坡穩(wěn)定性。由于壩體材料中尾砂的滲透系數(shù)較低，雨水不可能完全滲入到壩體內(nèi)部材料，大量雨水集結(jié)于坡體表面形成徑流被排出，因此本次計算不考慮徑流對計算結(jié)果的影響。

材料物理力學參數(shù)表 表1

穩(wěn)定性計算成果匯總表 表2

3.4 計算結(jié)果

3.4.1 滲流計算結(jié)果

利用Geostudio軟件的滲流模塊對上述建立的模型進行滲流計算，獲得尾礦壩在上述工況下的滲流計算結(jié)果，如圖2所示。

由圖2可知，連續(xù)降雨3天后壩體的總水頭相較于降雨一天后的總水頭有明顯提升。隨著雨水不斷滲入到壩體內(nèi)部，壩體的飽和度持續(xù)增加，壩體材料的滲透系數(shù)也隨之增加。降雨3天后壩體內(nèi)部浸潤線的高度相較于第一天有明顯增加，壩體失穩(wěn)的可能性也隨之增加。

3.4.2 穩(wěn)定性計算結(jié)果

圖2 連續(xù)降雨條件下尾礦壩總水頭圖

將上述滲流計算結(jié)果導(dǎo)入到GeoStudio軟件的穩(wěn)定性計算模塊，利用軟件內(nèi)置的簡化bishop法對壩體的穩(wěn)定性進行計算，結(jié)果如圖3所示。

圖3 連續(xù)降雨一天后尾礦壩壩體穩(wěn)定分析圖

由圖3可知，在連續(xù)降雨一天后壩體的安全系數(shù)為0.919，明顯小于《尾礦設(shè)施設(shè)計規(guī)范》（GB 50863-2013）[6]規(guī)定的最低安全系數(shù)1.25，尾礦壩壩體會發(fā)生失穩(wěn)。為了便于比較分析，同時給出了連續(xù)降雨2天、3天以及3種不同算法計算得出的安全系數(shù)值，見表2。從表2可以看出，隨著降雨時間的增加，尾礦壩壩體的穩(wěn)定系數(shù)進一步下降，連續(xù)降雨3天后，由簡化Bishop法計算得出的安全系數(shù)僅為0.766，明顯小于規(guī)范規(guī)定的最低安全系數(shù)，極易發(fā)生潰壩。

4 結(jié)論

本文運用GeoStudio軟件內(nèi)置的滲流及穩(wěn)定性計算模塊對連續(xù)降雨條件下的尾礦壩滲流及穩(wěn)定性情況進行分析，得出的主要結(jié)論如下：

①連續(xù)降雨條件下，壩體的總水頭明顯上升，壩體材料逐漸趨于飽和，浸潤線明顯抬升，進而導(dǎo)致壩體失穩(wěn)，且隨著降雨時間的增加，壩體的安全系數(shù)進一步下降；

②在連續(xù)降雨條件下，應(yīng)做到對尾礦壩位移的24小時監(jiān)測，及時采取措施應(yīng)對發(fā)生的突發(fā)情況，避免潰壩現(xiàn)象的發(fā)生。