張 航

(浙江水利水電學院，浙江 杭州 310018)

陜西省榆林市靖邊縣蔣家窯則水庫工程是國家水利改革發(fā)展“十三五”規(guī)劃確定實施的大型骨干攔沙水庫工程，水庫壩址位于靖邊縣中山澗鎮(zhèn)五道溝村紅柳河的干流上，上游涉及定邊縣的楊伏井和營盤山兩座水庫。該水庫的主要任務是攔截區(qū)間流域泥沙，減輕下游水庫的泥沙壓力，保障紅柳河庫壩群安全，進一步減少入黃泥沙，改善下游河道生態(tài)環(huán)境。

為了對蔣家窯則水庫的地質(zhì)環(huán)境進行準確把握，本文基于蔣家窯則水庫的實地調(diào)查資料，對庫區(qū)工程的基本地質(zhì)問題進行了討論，并對壩址工程進行了對比分析，對水庫的建設提供了一些參考建議。

1 蔣家窯則水庫工程概況

蔣家窯則水庫為工程等級Ⅱ等的大型水利工程，工程位置距離靖邊縣城約52km，距離梁鎮(zhèn)約14km。工程的主要建筑物包括大壩、溢洪道和放水洞等,均為2級建筑物，其他次要建筑物為3級，壩頂?shù)母叱虨?408～1410m，設定的防洪標準為100年一遇設計、2000年一遇校核。

工程區(qū)地處黃土高原北部，跨入毛烏素沙漠邊緣，地勢西南高、東北低，高程為1123～1823m，高差為300～500m，整體地貌較為復雜，大致可分為北部的風沙灘地區(qū)、中部的梁峁?jié)镜貐^(qū)和南部的丘陵溝壑區(qū)。

地質(zhì)上，工程區(qū)內(nèi)地層以第四系地層為主，第四系前地層僅在紅柳河上游溝底處有零星出露。在工程區(qū)內(nèi)沒有歷史地震活動的記錄，綜合構造特征以及地震動參數(shù)分析，認為工程區(qū)的區(qū)域構造穩(wěn)定性較好[1]。

2 水庫區(qū)工程及主要地質(zhì)問題

水庫區(qū)工程所在地的上游溝谷發(fā)育，地形較為破碎，而中、下游的地形簡單，河谷較為平直。紅柳河在水庫區(qū)的走向為近南北向，由于紅柳河的形成時間較晚，所以河流的堆積地貌較為不發(fā)育，形成的深切河谷深度約30～50m。水庫區(qū)存在的不良地質(zhì)問題主要為崩塌和滑坡。除此之外，地表沖溝、落水洞及岸邊裂縫也較為多見，本節(jié)選取了水庫的常見地質(zhì)問題對蔣家窯則庫區(qū)進行分析。

2.1 水庫滲漏

水庫滲漏是指水庫蓄水后，水庫中的蓄水沿巖石的孔隙、裂隙斷層、溶洞等向庫岸分水嶺外的溝谷低地滲漏的現(xiàn)象[2]。

經(jīng)勘查發(fā)現(xiàn)，蔣家窯則水庫區(qū)的地形封閉，河谷兩岸山體雄厚，庫盆周邊地形較為完整，庫區(qū)上游段的地下水位高于河水的水位，而庫區(qū)下游段的河水低于地下水，庫區(qū)在蓄水初期可能存在一定的滲漏，但隨著土體的逐漸飽和，滲漏現(xiàn)象也會逐漸消失，因此可以認為庫區(qū)不存在永久滲漏，加之水庫以攔沙為主，因此可以不考慮水庫的滲漏問題。

2.2 庫岸穩(wěn)定

水庫在達到正常蓄水位1403m時，所形成的庫周長度約為144km，庫岸兩側均為第四系松散堆積層，兩岸斜坡自然坡角35°～45°，主要巖性為粉質(zhì)壤土，因此在水庫正常蓄水時，庫岸受到?jīng)_刷易發(fā)生塌岸及滑坡[3]。

在對庫岸穩(wěn)定性進行分析時，首先明確蔣家窯則水庫的主要功能為攔沙，當洪水來臨時，攔截洪水，待泥沙沉淀后，將清水放空，因此水庫常年保持低水位運行。圖1為預測的水庫水位和淤泥面隨時間的變化曲線。

圖1 水庫水位和淤泥面變化曲線

水庫的水深可由水庫水位和淤泥面計算得出：

H=H1-H2

(1)

式中H——水庫水深,m；

H1——水庫水位,m；

H2——淤泥面,m。

圖2為由圖1計算出的水深隨時間的變化情況。由圖2可知，水庫的蓄水位隨著時間逐漸抬高，淤泥面每年也隨之抬高。由于水位上升速度較緩，水深較淺，所以塌岸呈逐年坍塌，塌岸量較小。

圖2 水庫水深變化曲線

2.3 塌岸寬度預測

塌岸寬度根據(jù)類比圖解法進行預測，蔣家窯則水庫部分庫區(qū)前身為新橋水庫庫區(qū)，兩岸較為穩(wěn)定[4]。通過調(diào)查新橋水庫庫區(qū)兩岸水上和水下的天然休止角后，選用水下天然休止角19°、水上天然休止角50°進行圖解分析。

將庫區(qū)斷面位置進行劃分，為了表達清晰，此處采用簡易編號，各斷面編號及代表長度見表1，預測的各斷面坍塌寬度繪制于圖3中。

表1 庫區(qū)斷面編號及代表長度

續(xù)表

圖3 庫岸坍塌寬度

水庫蓄水至1405.6m高程時，主河道紅柳河距上壩址24.1km范圍內(nèi)，左岸塌岸寬度為5～70m，右岸塌岸寬度為6～65m。支溝邊墻渠溝口以上12.5km，左岸塌岸寬度范圍為27～75m，右岸塌岸寬度范圍為13.0～65.0m。支溝周灣溝口以上12km，左岸塌岸寬度范圍為40～66m，右岸塌岸寬度范圍為8.0～70.0m。老莊溝溝口以上2.0km，左岸塌岸寬度范圍為12.4m，右岸塌岸寬度范圍為13.3m。楊伏井、高溝畔、王河嘴溝不存在塌岸。

預測時采用坡度變化較小且具有代表性的剖面，用其蓄水位上升中期塌岸方量作為平均塌岸量，乘以水庫運行年限30年，確定最終塌岸方量。根據(jù)計算，水庫蓄水至興利水位1405.6m時，上壩址預估塌岸方量2688萬m3，下壩址預估塌岸方量2856萬m3。

該水庫塌岸方量較大，塌岸寬度局部較寬，且塌岸主要集中在近壩區(qū)的兩岸陡坡段，但考慮到水庫每年的蓄水時間僅1個月左右，水深僅2.3～5.8m，平均每年的塌岸量不大，不會產(chǎn)生大體積的速滑。

2.4 庫區(qū)滑坡穩(wěn)定性分析

根據(jù)實地調(diào)查資料，在水庫區(qū)發(fā)育滑坡中選取了三處滑坡進行穩(wěn)定性分析，其他滑坡分析方法相同，三處滑坡的參數(shù)和穩(wěn)定性初判見表2。

表2 庫區(qū)部分滑坡參數(shù)

在對滑坡進行穩(wěn)定性計算時，采用了剩余推力法，將滑坡體劃分為若干土條塊，然后建立作用在土條上的力或力矩的平衡方程，進而確定穩(wěn)定性系數(shù)[5]，理論計算公式為

(2)

式中ψj——第i塊段的剩余下滑力傳遞至第i+1塊段時的傳遞系數(shù)(j=i)；

Wi——第i條塊的重量，kN/m；

Ci——第i條塊黏聚力,kPa；

φi——第i條塊內(nèi)摩擦角,(°)；

Li——第i條塊滑面長度,m；

αi——第i條塊水平面與滑面的夾角,(°)；

A——地震加速度(重力加速度g)；

K——穩(wěn)定系數(shù)。

圖4～圖6為三個滑坡的分條剖面圖[6]，計算得出的穩(wěn)定性系數(shù)K<1時，認為滑坡為不穩(wěn)定，當1≤K<1.05時，認為滑坡欠穩(wěn)定，當1.05≤K

表3 庫區(qū)部分滑坡穩(wěn)定性計算

圖4 滑坡1穩(wěn)定性計算分條剖面

圖5 滑坡2穩(wěn)定性計算分條剖面

圖6 滑坡3穩(wěn)定性計算分條剖面

根據(jù)分析，庫區(qū)滑坡天然狀態(tài)基本穩(wěn)定。水庫的正常蓄水位以逐年遞增的方式抬高，上升速度較緩，庫水位隨淤積面抬升而升高，滑坡體下滑空間相對變小，滑坡2和滑坡3復活可能性不大。水庫建成蓄水后，滑坡1方為一中型滑坡，距上壩址1.8km，土體飽和會使滑坡1的穩(wěn)定性受到影響，滑坡失穩(wěn)后對大壩安全有一定的影響，但由于水庫蓄水位是逐年漸進式抬高，加之水位較淺，滑坡破壞形式以逐年小規(guī)模坍塌為主，不會產(chǎn)生大的涌浪，因此不會對樞紐安全有太大影響。

2.5 水庫淹沒及浸沒

水庫區(qū)位于溝谷內(nèi)，兩岸地形較為封閉，在正常蓄水位1405.6m時，不會超越兩岸的耕地，因此兩岸的耕地不存在淹沒現(xiàn)象。但上游河道內(nèi)的耕地在水庫蓄水后將會出現(xiàn)淹沒現(xiàn)象。而在庫尾河道處存在浸沒現(xiàn)象，河床粉土的毛管水上升高度Hk為1.5m，安全超高值ΔH為0.5m，因此浸沒地下水埋深臨界值Hcr為2.0m。正常蓄水位為1405.6m，所以庫水位在1407.6m時不會發(fā)生浸沒，以此預測的水庫區(qū)浸沒面積見表4。

表4 水庫區(qū)浸沒預測

3 壩址工程

3.1 壩址區(qū)基本地質(zhì)條件

上、下壩址位于蔣家窯則村紅柳河河段，壩區(qū)地勢南高北低，地形高差約120m。河谷呈“U”形，深度約45m，寬度約160～210m，兩側岸坡自然坡度約35°～45°，灘面高程為1377.4～1377.8m。壩區(qū)無基巖出露，地下水以第四系松散堆積層孔隙潛水為主，類型較為簡單。

在壩線下游約1.5km處，跨河段公路旁的粉質(zhì)黏土層上預先設置抽水試驗井兩眼，深度都為35m，兩井的間距為8m，井徑為300mm。兩井分別進行抽水試驗，并互為觀測孔，單井抽水按3個降深控制。按《水文地質(zhì)手冊》中非完整井計算公式如下[7]：

(3)

式中Q——出水量,m3/d；

r1、r2——觀測孔1與觀測孔2之間的距離,m；

s1、s2——觀測孔1與觀測孔2的水位下降值,m；

ξ1、ξ2——觀測孔1與觀測孔2的滲透阻力系數(shù)；

H——抽水前潛水層的厚度,m。

水位恢復法計算公式如下[8]：

(4)

式中s′——剩余降深，m；

t——抽水累計時間，min；

t′——停止抽水后恢復水位時間，min。

抽水試驗計算成果見表5[9]。由表5可知，河漫灘粉質(zhì)黏土和粉質(zhì)壤土的滲透系數(shù)K介于1.29×10-5～1.02×10-4cm/s之間，滲透系數(shù)K平均值為3.4×10-5cm/s。

表5 抽水試驗滲透系數(shù)成果

3.2 上、下壩址比選

3.2.1 壩址工程地質(zhì)條件對比

上壩址壩軸線處河床高程1377.8m，河谷寬約166m，環(huán)境水對普通水泥具硫酸鹽型強腐蝕性，對混凝土結構中的鋼筋具中等腐蝕性，對鋼結構具中等腐蝕性。設計的壩頂長度為310.95m，壩高為37.8m，壩頂?shù)母叱虨?413.8m，正常蓄水位為1405.6m。

下壩址位于上壩址下游900m處，設計壩頂長度372.9m，壩高37.5m，壩頂高程1413.5m，正常蓄水位1405.2m。環(huán)境水對普通水泥具硫酸鹽型強腐蝕性，對混凝土結構中的鋼筋具中等腐蝕性，對鋼結構具中等腐蝕性。

上、下壩址的左壩肩均不存在邊坡穩(wěn)定和繞壩滲漏的問題，都具非自重濕陷性，濕陷等級為Ⅱ級，上壩址左壩肩濕陷下限13.5m，下壩址為12.5m。

上壩址的右壩肩不存在邊坡穩(wěn)定、繞壩滲漏問題，具非自重濕陷性，濕陷等級為Ⅱ級，濕陷下限12.5m。而下壩址的右壩濕陷等級為Ⅰ級，濕陷下限深度3.5m，存在繞壩滲漏，應采取截滲墻等方式進行防滲處理。

上、下壩址的壩基滲透破壞形式均為流土，不存在滲透破壞問題及壩基滲漏。

3.2.2 壩址工程地質(zhì)問題對比

a.岸坡穩(wěn)定。上、下壩址兩岸壩肩均為土質(zhì)邊坡，邊坡巖性主要為粉質(zhì)壤土，坡面基本平整，有少量植被。上壩址處邊坡現(xiàn)狀未發(fā)現(xiàn)崩塌、滑坡存在，整體穩(wěn)定性較好。下壩址左岸邊坡未發(fā)現(xiàn)崩塌、滑坡存在，而在右岸邊坡壩肩處存在一小型滑坡，但對整體穩(wěn)定性影響不大。

b.砂土液化。上、下壩址工程區(qū)的地震基本烈度均為Ⅵ度，因此可以不考慮砂土地震液化的影響。

c.壩基及繞壩滲漏。上壩址在水庫蓄水位高程以下不存在強透水的粗粒土層，除粉質(zhì)壤土外，兩岸岸坡從壩頂高程以下都屬弱透水層，結合壩肩開挖清除粉質(zhì)壤土后，一般不存在壩基及繞壩滲漏問題。

下壩址的左岸同上壩址處情況相同，不存在壩基滲漏及左壩肩繞壩滲漏問題，而右岸岸坡存在中等透水層，應進行防滲處理。

d.滲透穩(wěn)定。我國規(guī)范《碾壓式土石壩設計規(guī)范》(SL 274—2020)中分別給出了流土與管涌的臨界水力坡降計算公式，可以有效地進行滲透穩(wěn)定計算。流土型宜采用下式計算[10]：

Jcr=(Gs-1)(1-n)

(5)

式中Jcr——土的臨界水力坡降；

Gs——土粒密度與水的密度之比；

n——土的孔隙率。

管涌型或過渡型宜采用下列計算：

(6)

式中d5、d20——占總土重5%和20%的土粒粒徑,mm。

上壩址處河谷粉質(zhì)黏土不均勻系數(shù)Cu為24.8，細粒含量大于35%，依據(jù)《水利水電工程地質(zhì)勘察規(guī)范》(GB 50487—2008)的附錄G，判定滲透破壞形式為流土。因此采用式(5)進行計算，得到的結果為0.884，安全系數(shù)取2.5，因此允許水力坡降J允=0.35。

上壩壩基的寬度為267m，河床高程為1377.5m，正常蓄水位為1405.6m，28.1m水頭時壩體上、下游實際水力坡降J實=0.105,根據(jù)土的顆分成果、物性成果，兩岸岸坡土及河道淤積土允許水力坡降J允=0.35，J允>J實，故不存在滲透破壞問題。

下壩壩基寬265m，河床高程1377.5m，正常蓄水位1405.2m，27.7m水頭時壩體上、下游實際水力坡降為J實=0.105，根據(jù)土的顆分成果、物性成果，兩岸岸坡土及河道土允許水力坡降J允=0.35，J允>J實，故不存在滲透破壞問題。

e.變形穩(wěn)定。上、下壩址處壩基的粉質(zhì)黏土呈軟塑—可塑狀，壓縮模量及承載力低，屬中壓縮性土，存在變形穩(wěn)定問題，不能直接作為天然地基，建議對壩基進行排水固結處理，可采用振沖碎石樁進行處理，根據(jù)設計計算沉降深度和允許沉降深度，處理到適宜深度。

3.2.3 壩址比選

除主要壩體外，壩址處的建筑物還有放水洞和放水塔等，將這些建筑物的勘察結果對比與主要壩體勘查結果對比一并列于表6中。

表6 壩址工程地質(zhì)條件比較

兩壩址地貌單元相同，層位相近，均具備修建當?shù)夭牧蠅蔚牡刭|(zhì)條件。從工程地質(zhì)問題進行分析，主要的地質(zhì)問題來自于岸坡的穩(wěn)定性以及壩基的滲漏。上壩址庫區(qū)淹沒較少，壩頂長度較短，且不存在繞壩滲漏問題，地質(zhì)條件優(yōu)于下壩址，因此推薦上壩址為選定壩址。

4 結 論

a.蔣家窯則水庫的主要地質(zhì)問題為岸坡的滑坡和崩塌，預測庫區(qū)塌岸總方量上壩址為2688萬m3，下壩址為2856萬m3?；缕茐男问揭灾鹉晷∫?guī)模坍塌為主，對樞紐的安全性影響較小。處理類似工程問題時，可以采用加強庫區(qū)岸坡綠化進行防護。

b.水庫建成蓄水后隨著庫水位的不斷抬高，其穩(wěn)定性將發(fā)生變化。蓄水后，滑坡相比于天然狀態(tài)下更不穩(wěn)定，因此在分析滑坡時，應對蓄水后的滑坡穩(wěn)定性多加考慮。

c.壩基的滲透穩(wěn)定性可以采用臨界水力坡降計算公式進行計算分析。上、下壩址的壩基不存在滲透破壞問題及壩基滲漏，但實際的工程壩基土層呈軟塑狀，土層欠固結，存在變形穩(wěn)定問題，地下水位埋藏淺，排水固結慢，不能直接作為天然壩基。對于此類壩基應進行排水固結處理。

d.在進行類似的壩址選擇時，若壩址的地貌單元相同，則應從工程地質(zhì)條件進行對比，其中主要的地質(zhì)問題來自于岸坡的穩(wěn)定性以及壩基的滲漏。