李景娟，謝羅峰， 段祥寶

(1.南京水利科學(xué)研究院，南京 210029；2. 河海大學(xué)水利水電學(xué)院，南京 210098；3. 水文水資源與水利工程科學(xué)國家重點試驗室，南京 210029)

榆樹溝灰場是山西省樟澤發(fā)電廠用于貯灰并利用粉煤灰淤地造田的山谷型灰場，機組除灰均采用水力除灰方式?；覉龅墓こ桃?guī)模為小(一)型灰場，工程等別為Ⅳ等，灰場屬于永久性建筑物中的主要建筑物，工程等級為4級。初期壩于2007年年初投入運行，灰場初期庫容105萬m3，壩底高程895 m，壩頂高程920 m，壩長300 m，壩頂寬5 m。初期壩外坡坡度1∶2.75，內(nèi)坡坡度1∶2.5，壩腳設(shè)堆石棱體，壩基設(shè)排水褥墊，上游壩面設(shè)貼坡反濾層，下游壩面設(shè)干砌石護(hù)面。初期壩壩體距下游村莊約320 m。自灰場運行以來，基本每年都加筑加高子壩，至2013年11月已加高至五級子壩，各級子壩均高為4 m，壩頂高程為935 m，現(xiàn)仍在運行中。運行7年時間，從未進(jìn)行任何 有關(guān)安全性評估工作。近年隨著農(nóng)村的擴建，壩體距民宅已不足300 m，為了下游居民的安全，有必要對灰壩進(jìn)行應(yīng)力變形和滲流安全分析，從而預(yù)測灰面高程增加的可行性。

1 地質(zhì)概況與地震

灰場范圍谷坡地層大部分為第四系上更新統(tǒng)(Q3)黃土(粉土)、中更新統(tǒng)(Q2)棕紅色黃土(粉質(zhì)黏土)，局部出露石炭系泥頁巖與砂巖地層。泥頁巖地層滲透系數(shù)較小，預(yù)測水力貯灰場運行后整體不會產(chǎn)生太大的滲漏。灰場范圍內(nèi)地下水為石炭系砂巖裂隙水和灰?guī)r裂隙溶隙水，局部出露石炭系的灰?guī)r，滲透系數(shù)較大，因出露面積局部分布于壩址區(qū)，為防滲可以用黏土鋪墊或用水泥漿噴抹。地質(zhì)剖面如圖1所示。

該灰場的工程抗震設(shè)防類別為丁類。根據(jù)《中國地震動峰值加速度區(qū)劃圖》(GB18306-2001)，灰場地震動峰值加速度為0.10 g(對應(yīng)地震烈度為Ⅶ度)，地震動反應(yīng)譜特征周期0.40 s。

2 安全監(jiān)測儀器布設(shè)及觀測結(jié)果

根據(jù)規(guī)范要求結(jié)合灰場實際地形地質(zhì)條件，在該灰場東側(cè)穩(wěn)定區(qū)域設(shè)有3個沉降基準(zhǔn)點，5個位移基準(zhǔn)點，8個位移觀測點，沉降觀測點均勻布設(shè)在位移觀測點基礎(chǔ)上，共計8個，壩頂部5個，下部馬道3個，與東側(cè)沉降基準(zhǔn)點組成水準(zhǔn)網(wǎng)形，2009年在子壩處又加了7個監(jiān)測點。

圖1 壩體最大斷面0+150Fig.1 The most dangerous section 0+150 of the dam

從灰壩的7 a沉降觀測可知，測點H4的總沉降量最大，累積沉降量為299.4 mm，其中年沉降量最大值為79.7 mm，是由于該年多雨季節(jié)導(dǎo)致整個壩體沉降量較之前幾年明顯加大，但是量值小，沉降趨于穩(wěn)定。馬道處的水平位移和沉降基本不變，位于壩頂?shù)臏y點H5水平位移最大為16 mm，指向壩外，而H6、H7、H8水平位移最大為19 mm，指向庫內(nèi)，說明初期壩壩頂可能存在拉裂縫[1]。

3 滲流有限元計算分析

灰壩自運行以來，壩體壩基無明顯滲漏。壩體前壩坡設(shè)混凝土護(hù)坡，同時由于風(fēng)浪的沖刷，作為壩體的防滲體面板有個別斷裂，表層的橡塑油膏脫落，嵌縫止水個別發(fā)生分離?，F(xiàn)場查看壩后段情況可知，壩體滲漏量較小，在可控范圍內(nèi)。

由于在滲流計算時，壩體內(nèi)自由面或浸潤面未知，因此計算按非穩(wěn)定滲流考慮，計算時將浸潤面以上的非飽和區(qū)以及自由面以下的飽和區(qū)同時予以考慮，計算所得的零孔隙壓力線即為浸潤線。計算時對壩底排水褥墊及壩腳排水棱體的滲透系數(shù)根據(jù)以往經(jīng)驗選取，庫內(nèi)灰渣由于沉積容易形成滲透的各向異性，一般是水平向滲透系數(shù)大于垂直向滲透系數(shù)，參考試驗結(jié)果及地質(zhì)勘查資料選取(見表1)。

將計算結(jié)果與已有的測壓管實測值進(jìn)行對比(圖2中已標(biāo)出測壓管實測值)，擬合情況與實測情況相近，誤差在3%以內(nèi)，說明該計算模型與參數(shù)的選取合理，可用該模型來預(yù)測子壩加高后的滲流安全性態(tài)。計算選取最大斷面0+150(見圖2和圖3)。

表1 各土層物理力學(xué)參數(shù)Tab.1 The physical parameters of each soil layer

圖2 0+150斷面五級子壩等勢線分布圖Fig.2 The equipotential line of the 0+150 fifth stage dam

圖3 0+150斷面六級子壩等勢線分布圖Fig.3 The equipotential line of the 0+150 sixth stage dam

從0+150斷面的流場可以看出，在無干灘運行條件下五級子壩內(nèi)浸潤線位置也不高，現(xiàn)場觀測時也未見到壩體下游有滲水，而且壩腳排水棱體里的出滲量也較小，說明目前為止壩體滲流穩(wěn)定是安全的。初期壩上游貼坡排水體與壩底排水褥墊的功能十分重要，對于降低壩體浸潤線起著主導(dǎo)作用，因此在貯灰場的運行過程中，要確保上游貼坡排水與壩基排水褥墊及坡腳排水棱體有效的工作。由滲流計算結(jié)果可知，壩體下游滲透坡降小于允許滲透比降，五級子壩時單寬滲流量為5.6×10-5m3/(s·m)，加高至六級子壩時滲流量為5.5×10-5m3/(s·m)，差異不大。壩面高程繼續(xù)增加時，可考慮在子壩前鋪設(shè)一定長度的干灘，能有效降低貯灰區(qū)的浸潤線，以防子壩下游坡大面積出逸[2]。

4 滲流作用下應(yīng)力計算分析

將SEEP計算結(jié)果代入SIGMA中來分析體內(nèi)應(yīng)力分布規(guī)律與變形成因，根據(jù)現(xiàn)有的試驗參數(shù)結(jié)果，結(jié)合變形監(jiān)測資料，利用Geo-Studio SIGMA模塊彈塑性模型，對土壩內(nèi)部應(yīng)力進(jìn)行數(shù)值模擬計算[3]。

灰壩壩體內(nèi)部有無裂縫產(chǎn)生可采用傾度法計算，設(shè)壩體某一平面上任意兩點的總沉陷差為Δh，兩點間的距離為ΔL，兩者的比值為傾度d，即：

(1)

傾度d實質(zhì)上是兩點間的不均勻沉陷斜率。當(dāng)傾度d大于A、B間土體的破壞臨界傾度dc時，即表明在兩點之間可能產(chǎn)生裂縫。經(jīng)過國內(nèi)幾十座土石壩的調(diào)查分析，當(dāng)傾度d小于1%時壩體不會產(chǎn)生裂縫，在1%～2%之間時壩體可能會有裂縫產(chǎn)生，大于2%時壩體會有裂縫產(chǎn)生[4]。

圖4和圖5給出了現(xiàn)狀五級子壩及加高至六級子壩后的壩體內(nèi)部變形分布，壩頂處水平變形一側(cè)指向庫內(nèi)，一側(cè)指向壩外，產(chǎn)生裂縫的可能性較大，根據(jù)傾度法可知，斷面0+150目前最大傾度為0.34%，加高至六級子壩時最大傾度為0.79%，壩體不會產(chǎn)生裂縫，說明施工時對壩基的濕陷性黃土處理得比較徹底，上下游濕陷的時間與量值差異小，不均勻沉降小，不至于產(chǎn)生裂縫。

圖4 五級子壩變形云圖Fig.4 The deformation contour of the fifth stage dam

圖5 六級子壩變形云圖Fig.5 The deformation contour of the sixth stage dam

5 滲流作用下邊坡穩(wěn)定計算分析

5.1 靜力穩(wěn)定計算

貯灰場壩坡靜力穩(wěn)定分析采用簡化畢肖普法。本次計算采用Geo-Studio SLOPE模塊的等效線性模型，計算時將SEEP模塊算得的孔隙水壓力及浸潤線等結(jié)果代入SLOPE模塊來進(jìn)行極限平衡抗滑穩(wěn)定計算。

計算結(jié)果如表3所示，五級子壩在無干灘情況下正常運行的最小抗滑穩(wěn)定系數(shù)為2.01，加高至六級子壩時最小抗滑穩(wěn)定系數(shù)稍有降低為2.00，說明壩體目前是穩(wěn)定的，不會出現(xiàn)整體滑動現(xiàn)象。

5.2 動力計算

動力計算根據(jù)《水工建筑物抗震設(shè)計規(guī)范》的要求采用擬靜力法。擬靜力法是將壩體各點地震動力作用所引起的慣性力作為靜力作用在質(zhì)點上計算抗滑穩(wěn)定性，方法簡單，但是難以反映土壓力的分布和土體的動力特性。本次地震模擬采用Geo-Studio QUAKE模塊，運用等效線性模型模擬出灰壩在地震前的初始應(yīng)力狀態(tài)，然后在動力分析中輸入唐山遷安地震記錄曲線，持續(xù)時間取振幅較大的10 s，計算時將其最大加速度調(diào)整到0.1 g，設(shè)計反應(yīng)譜特征周期為0.4 s，只考慮水平方向的地震波作用[5](見圖6)。

圖6 地震加速度進(jìn)程曲線Fig.6 Horizontal acceleration earthquake record

邊界條件為：壩基底部水平、豎直方向均約束，壩基兩側(cè)豎直方向約束，正常蓄水位932.0 m。

由動力法對貯灰壩地震情況的模擬可知，地震引起的最大水平位移為48 mm，整個初期壩體發(fā)生液化的可能性不大，但是上游堆灰處密實度不足，在地震工況下孔隙水壓力上升，土體有效應(yīng)力降低，土體處于懸浮狀態(tài)以至于變形增大，而地震引起的超靜孔隙壓力不能及時消散導(dǎo)致粉煤灰噴水冒沙產(chǎn)生液化，下游壩體壩基產(chǎn)生液化可能性不大[6](見圖7和圖8)。

圖7 五級子壩液化區(qū)域云圖Fig.7 Liquefaction zone of the fifth stage dam

圖8 最大動剪應(yīng)力Fig.8 Peak dynamic shear stresses

5.3 動力穩(wěn)定計算

由于地震作用，貯灰壩組成材料的物理力學(xué)性質(zhì)較靜力狀態(tài)有很大的變化，尤其是土體的內(nèi)摩擦角和黏聚力有明顯變化，都較地震前有所降低。將地震過程中的動力計算結(jié)果代入SLOPE模塊作為液化后壩坡穩(wěn)定分析的初始條件，液化后穩(wěn)定計算采用的參數(shù)如表2所示。

從表3可以看出，五級子壩在滲流作用下的穩(wěn)定系數(shù)為2.01，在滲流加地震的工況下穩(wěn)定系數(shù)為1.41，可見地震作用對壩坡穩(wěn)定性的影響比較嚴(yán)重[7]。

6 結(jié) 語

基于Geo-Studio分析軟件，結(jié)合原有的大壩監(jiān)測資料，利用流固耦合原理計算了滲流作用下的應(yīng)力變形分布流固耦合原理計算了滲流作用下的應(yīng)力變形分布，從而判斷壩體內(nèi)有無裂縫。模擬了地震作用下的超孔靜壓、液化情況，最終計算出滲流工況和滲流加地震工況下的大壩邊坡最小安全系數(shù)，判斷灰壩的安全性。

表2 動力計算參數(shù)表Tab.2 Dynamic calculation parameter

表3 壩坡抗震穩(wěn)定計算結(jié)果Tab.3 The results of slope stability under earthquake

壩體加高至六級子壩時，不會產(chǎn)生裂縫，壩內(nèi)浸潤線位置不高，初期壩上游貼坡排水及壩底排水褥墊起著重要的作用，因此要確保排水系統(tǒng)的有效運行，壩前可設(shè)置適當(dāng)長度的干灘。地震設(shè)防烈度七度下，上游庫內(nèi)灰渣在一定范圍內(nèi)可能會發(fā)生液化，但是大壩整體是穩(wěn)定的。綜合上述計算分析，子壩加高至六級是可行的。

□

[1] 劉富蟬，黃曉燕. 某電廠陽坡灰場安全評價及其治理措施評述[J].電力學(xué)報，2008,23(2):101-105.

[2] 毛昶熙.滲流計算分析與控制[M].2版.北京：水利電力出版社，2003:510-528.

[3] 池龍哲，吳世勇.常用滲流有限元計算軟件比較分析[J].中國農(nóng)村水利水電，2014,(6):186-188.

[4] 劉富蟬，黃曉燕.某電廠陽坡灰場安全評價及其治理措施評述[J].電力學(xué)報，2008,23(2):101-106.

[5] 王 宇，謝羅峰，段祥寶.于橋水庫大壩地震液化判別和動力分析[J].中國農(nóng)村水利水電，2013,(8):133-137.

[6] 錢家歡，殷宗澤. 土工原理與計算[M]. 2版. 北京：中國水利水電出版社，1996:558-583.

[7] 毛昶熙，段祥寶.滲流數(shù)值計算與程序應(yīng)用[M].南京：河海大學(xué)出版社，1999:75-76.