孔佑龍

（貴陽市水利水電勘測設(shè)計研究院有限公司，貴州 貴陽 550081）

土石壩作為重要的擋水建筑物，其穩(wěn)定性受建壩材料及建筑時間長等因素影響[1]。合理準(zhǔn)確地分析壩體滲流穩(wěn)定性能夠保障壩體正常運行及下游居民生命和財產(chǎn)安全。目前，已有許多學(xué)者對壩體穩(wěn)定性分析進行了大量的研究，提出了很多不同的方法。阿布都沙拉木等[2]依據(jù)北疆某水庫大壩多年滲流資料進行了系統(tǒng)分析，明確了該水庫的滲透穩(wěn)定性。殷躍平[3]對三峽庫區(qū)滑坡防治工程設(shè)計中的滲透壓力問題進行了較為系統(tǒng)的研究。李輝等[4]建立土石壩在不同水位時的滲透特性模型。姚浩等[5]以邢臺市某土石壩為例，計算分析了該粘土心墻壩滲流和穩(wěn)定均滿足設(shè)計要求。本文以始建于1958 年的某老水庫為例，在壩體土工取樣試驗和水利樞紐現(xiàn)場檢查結(jié)果的基礎(chǔ)之上，運用GeoStudio 軟件結(jié)合Simplified Bishop 圓弧法對該庫土石壩滲流穩(wěn)定性進行計算分析，并提出了相應(yīng)的防滲措施建議。

1 工程概況

某水庫壩型為均質(zhì)土壩，始建于1958 年，于1965年擴建成小（1）型水庫，水庫總庫容192 萬m3，主要功能為灌溉，設(shè)計灌溉面積6920 畝，實際灌溉面積為5800 畝。壩址以上集水面積為9.3km2，其中主河道長6.44km，河道平均比降為21.1‰。樞紐工程由大壩、溢洪道、放水設(shè)施等組成。壩頂高程為1317.20m，最大壩高20.2m，壩軸線長為120m，壩頂寬2.1m，大壩上游迎水面坡比為1∶3.15；大壩下游背水面壩坡分為三段：上壩坡坡比1∶1.5，中壩坡坡比1∶2.7，下壩坡坡比1∶3.35，壩腳排水棱體坡比為1∶1.5。

大壩填筑材料為黃褐色粘土，呈可塑至軟塑狀，含巖屑顆粒，局部巖屑顆粒含量較多，粘性差，較松散。為了復(fù)核壩體填筑土料質(zhì)量，對鉆孔取芯的壩體填筑土料進行了室內(nèi)物理力學(xué)試驗。試驗結(jié)果表明：大壩壩中段壩體填筑土料局部滲透系數(shù)大于1×10-4cm/s，最大值為2.14×10-4cm/s；同時結(jié)合現(xiàn)場檢查發(fā)現(xiàn)大壩下游壩坡一級馬道出現(xiàn)的兩處滲漏帶。鑒于此，有必要對該水庫大壩進行滲流安全計算分析，了解大壩運行狀態(tài)，保障水庫各項功能的正常運行。

2 工程地質(zhì)條件

水庫屬高原侵蝕、溶蝕型中山～中低山地貌，地勢較平緩，河谷較寬。碳酸鹽巖區(qū)域地質(zhì)營力以溶蝕作用為主，基巖多為裸露及半裸露；碎屑巖區(qū)域地質(zhì)營力以侵蝕作用為主，一般覆蓋層較厚。庫區(qū)地形起伏較大，高程于1280m～1500m 之間，地形開闊，地形平均坡度15°～30°，局部存在陡坡、陡坎。

壩區(qū)河谷為順向谷，巖層傾向左岸，壩址區(qū)為淺切低中山河谷地貌。兩岸坡覆蓋層為褐色、黃褐色的粘土、砂土和碎石塊，厚0.5～3.5m 不等；左壩肩基巖為灰、深灰色厚層塊狀含硅質(zhì)巖團塊生物碎屑石灰?guī)r；壩基及右壩肩主要為灰、深灰色粘土巖、泥晶石灰?guī)r。根據(jù)鉆孔取芯揭示，壩基基巖為深灰色粘土巖夾灰?guī)r，基巖軟硬相間，差異風(fēng)化強烈，強風(fēng)化厚度一般3～5m，且強風(fēng)化層及弱風(fēng)化層節(jié)理裂隙發(fā)育，巖芯主要呈碎塊狀、短柱狀，透水性強。鉆孔壓水試驗顯示，壩基、壩肩的接觸帶基巖風(fēng)化強烈，透水性強，不起壓。弱風(fēng)化巖體透水率4.7～5.1Lu，屬弱透水層。

3 壩體滲流穩(wěn)定分析

3.1 計算條件

根據(jù)《防洪標(biāo)準(zhǔn)》（GB50201-2014）及《水利水電工程等級劃分及洪水標(biāo)準(zhǔn)》（SL252-2017）的規(guī)定，本工程水庫屬?。?）型水庫，工程等級為Ⅳ等，主要建筑物大壩、溢洪道、取水建物為4 級建筑物，次要建筑物（含臨時建筑物）為5 級。根據(jù)SL274-2020《碾壓土石壩設(shè)計規(guī)范》的規(guī)定，該大壩屬均質(zhì)土壩，壩基無軟弱夾層或軟弱帶等情況，故采用計及條塊間作用力的簡化畢肖普（Simplified Bishop）圓弧法進行壩坡抗滑穩(wěn)定計算，相應(yīng)的壩坡抗滑穩(wěn)定最小允許安全系數(shù)見表1。

表1 壩坡抗滑穩(wěn)定最小安全系數(shù)表

同時，土體的抗剪強度指標(biāo)選用是否合理，將直接影響壩坡抗滑穩(wěn)定計算成果的可靠性及合理性。因此，壩體填筑土料、排水棱體的堆石料的各物理力學(xué)參數(shù)指標(biāo)采用土工試驗成果中的計算參數(shù)建議值，見表2。

表2 物理力學(xué)參數(shù)指標(biāo)表

3.2 計算工況

經(jīng)計算，該水庫的校核洪水位和設(shè)計洪水位均低于現(xiàn)狀壩頂高程，且現(xiàn)狀壩頂高程滿足規(guī)范要求的壩頂超高要求。

水庫地區(qū)歷史上無中、強地震記錄，區(qū)內(nèi)無區(qū)域性斷裂和活動斷裂通過，地震基本烈度為Ⅵ度。故根據(jù)SL274-2020 的規(guī)定，結(jié)合工程實際情況，設(shè)置計算工況如下。

3.2.1 上游壩坡

工況一：正常蓄水位1314.00m 的穩(wěn)定滲流期；

工況二：上游設(shè)計洪水位1315.85m 形成的滲流情況；

工況三：上游校核洪水位1316.57m 形成的滲流情況。

3.2.2 下游壩坡

工況一：正常蓄水位1314.00m 的穩(wěn)定滲流期；

工況二：上游設(shè)計洪水位1315.85m 的穩(wěn)定滲流期；

工況三：上游校核洪水位1316.57m 的穩(wěn)定滲流期。

3.3 計算結(jié)果及分析

壩坡穩(wěn)定分析采用GeoStudio 11.4.0.18 軟件進行計算，按剛體極限平衡法中計及條塊間作用力的簡化畢肖普（Simplified Bishop）圓弧法進行壩坡抗滑穩(wěn)定計算，根據(jù)計算結(jié)果所示，各工況壩坡穩(wěn)定性均滿足要求，詳見表3。

表3 壩體上下游壩坡抗滑穩(wěn)定計算成果表

根據(jù)模擬分析結(jié)果，上下游壩坡條塊受力如圖1、圖4 所示，正常蓄水位、設(shè)計洪水位以及校核洪水位不同水位下壩體浸潤線的模擬分布情況如圖2、圖3 所示，其中上游水位處于設(shè)計洪水位、校核洪水位時，浸潤線從下游壩坡一級馬道上方逸出，正常蓄水位時浸潤線有從下游壩坡一級馬道逸出的跡象，結(jié)合現(xiàn)場檢查時發(fā)現(xiàn)的下游壩坡一級馬道出現(xiàn)的滲漏帶（上游水位1313.75m，接近正常蓄水位1314.00m），判斷出壩體滲流性態(tài)不安全。

圖1 上游壩坡條塊受力示意圖（單位：KN）

圖2 上游壩坡最危險滑移面位置圖

圖3 下游壩坡條塊受力示意圖（單位：KN）

圖4 上游壩坡最危險滑移面位置圖

同時根據(jù)壩坡抗滑穩(wěn)定計算成果顯示，該水庫大壩各工況最危險滑移面對應(yīng)的上、下游壩坡抗滑穩(wěn)定計算成果均大于SL274-2020《碾壓土石壩設(shè)計規(guī)范》中的允許最小安全系數(shù)，即大壩壩坡處于穩(wěn)定狀態(tài)。據(jù)此分析可知，由于大壩自1958 年建成后，經(jīng)歷過擴建以及數(shù)次除險加固，大壩壩體不斷加厚，大壩壩坡坡比也在歷次改造中不斷變緩，使得大壩壩坡在滲流性態(tài)不安全的情況下可以保持穩(wěn)定狀態(tài)。

4 防滲措施

通過對壩體、壩基鉆孔取芯、壓水試驗的成果及水庫樞紐各部位的現(xiàn)場檢查，結(jié)合上述的壩體滲流穩(wěn)定性計算，分別對各個工況的浸潤線、最危險滑移面進行了分析。綜上所述，該水庫大壩存在影響壩體穩(wěn)定的多個安全隱患問題，建議對大壩壩基、壩肩進行帷幕灌漿，壩體充填灌漿，形成連續(xù)的防滲帷幕；防止因壩體滲流穩(wěn)定性降低，改變壩體填筑粘土的物質(zhì)成分或結(jié)構(gòu)，進而改變壩體填筑粘土的內(nèi)摩擦角和凝聚力系數(shù)，減小顆粒間的有效應(yīng)力，導(dǎo)致壩體的抗剪強度降低。

5 結(jié)論

（1）土石壩滲流穩(wěn)定性分析是土石壩安全校核的重點之一，本文通過對某庫土石壩上下游壩坡在正常蓄水位、設(shè)計洪水位和校核洪水位的滲流穩(wěn)定分析，得出該土石壩在不同工況下滲流和穩(wěn)定均滿足設(shè)計要求。綜上所述，該方法可用于土石壩的滲流穩(wěn)定分析，為今后土石壩滲流分析的相關(guān)問題提供一定的借鑒。

（2）根據(jù)實際情況及計算分析結(jié)果，建議對大壩壩基和壩肩進行帷幕灌漿，壩體充填灌漿，形成連續(xù)的防滲帷幕；以此改變壩體填筑粘土的內(nèi)摩擦角和凝聚力系 數(shù)，提高壩體的抗剪強度，加強顆粒間的有效應(yīng)力，進一步地提高壩體的穩(wěn)定性，從而使壩體達(dá)到滲流穩(wěn)定的要求。