郝永剛

（山西省水利水電勘測設計研究院 山西太原 030024）

1 工程概況

北掌水庫位于山西省汾西縣對竹河一支溝內，總庫容465.4萬m3，壩型為均質土壩，壩頂高程973.5 m，最大壩高46 m，水庫正常蓄水位971.02 m，庫水回水長度約2.0 km。

庫區(qū)地貌類型為梁峁狀黃土地貌，庫區(qū)受長期流水侵蝕切割，地形較復雜，溝谷發(fā)育，黃土梁與溝谷交錯排列。在水庫回水范圍內，兩岸多為土質岸坡，僅在庫尾基巖裸露。庫區(qū)左岸壩址至三岔口段，自然邊坡15°～20°，三岔口至庫尾段，自然邊坡 30°～60°。庫區(qū)右岸壩址至三岔口段，自然邊坡50°～70°。岸坡上部為上更新統(tǒng)（Q3pl）淡黃色低液限粉（黏）土，厚3～10 m，下部為中更新統(tǒng)（Q2pl）淡紅色低液限黏土夾鈣質結核層構成。水庫蓄水后兩岸土質岸坡將會產生不同程度的塌岸。

兩岸土層的物理力學指標見表1。上部Q3pl低液限粉（黏）土天然含水率10.6%～14.0%，平均值12.4%，天然密度 1.37～1.58 g/cm3，平均密度 1.50 g/cm3，干密度1.23～1.43 g/cm3，平均值 1.34 g/cm3，孔隙比 0.890～1.199，平均值 1.027，飽和抗剪凝聚力 c=6.6～17.0 kPa，平均值 12.3 kPa，內摩擦角 Φ=20.0°～22.0°，平均值21.0°。下部Q2pl低液限黏土天然含水率7.5%～18.1%，平均值14.6%，天然密度1.39～1.92 g/cm3，平均密度1.70 g/cm3，干密度 1.29～1.67g/cm3，平均值 1.48 g/cm3，孔隙比0.633～1.080，平均值0.837，飽和抗剪凝聚力c=10.4～22.1 kPa，平均值 15.3 kPa，內摩擦角 Φ=19.0～24.0°，平均值 21.5°。

2 水庫塌岸影響因素

水庫塌岸影響因素主要有地形地貌、巖層結構和巖性、水的作用、波浪作用、沿岸流及凍融作用、水流沖刷及水庫淤積等。

水庫所處地貌單元，岸坡形態(tài)，包括岸高、岸的坡度、水下岸形、岸坡的曲率及岸坡切割程度等對塌岸的型式、速度、塌岸量和淺灘的形態(tài)有很大的影響。一般岸形彎曲地段，容易發(fā)生塌岸，凸岸地段更易受到?jīng)_蝕，塌岸速度快。山區(qū)河谷型水庫一般岸坡較高，塌岸較寬，但速度較小。庫區(qū)陡坡岸段塌岸嚴重、寬度范圍也大。

表1 庫區(qū)岸坡原狀土拌物理學性質統(tǒng)計表

庫岸巖層的層厚、上下層序、各層的出露位置及其物質組成和性質，直接控制岸壁坍塌的寬度、速度和型式，還決定著淺灘的形狀、寬度和坡角。一般黃土的粉土粒組含量大，孔隙率高，崩解速度快，易形成快速、強烈的坍塌。且水下淺灘坡角小。

水庫回水影響著庫岸地下水動態(tài)變化，也對庫岸土體的物理、力學、水理性質起著很重要影響。庫水的變化幅度與持續(xù)時間對淺灘寬度、坡角影響強烈，一般庫水變化幅度越大，淺灘越寬。庫水持續(xù)時間越長，淺灘坡度越小。

波浪作用主要表現(xiàn)為擊岸浪對岸壁土體的淘刷和磨蝕，對塌落物質進行搬動，從而加速塌岸過程。

沿岸流主要對凸岸、岸嘴進行沖刷，加速塌岸過程。北方水庫凍融作用破壞岸坡土層結構，使土體發(fā)生裂縫和坍塌。庫岸的物理地質作用，如風化、滑坡現(xiàn)象和地表水沖刷，也會加速塌岸過程。黃土地區(qū)水庫蓄水后泥沙淤積，會導致河床抬高，水深減小，波浪淘刷作用減弱，進行減小塌岸寬度和速度。

3 水庫塌岸預測

3.1 預測方法

采用二段預測法和卡丘金預測法選取典型斷面進行預測計算。

1）兩段法：

具體原理為：預測塌岸線由水下穩(wěn)定岸坡線和水上穩(wěn)定岸坡線的連線組成時，水下穩(wěn)定岸坡線由原河道多年最高洪水位及水下穩(wěn)定坡腳α確定；水上穩(wěn)定岸坡線由設計洪水位和毛細水上升高度H’及水上穩(wěn)定坡腳β確定（圖1）。以原河道多年最高洪水位與岸坡交點A為起點，以α為傾角繪出水下穩(wěn)定岸坡線，該線延伸至設計洪水位加毛細水上升高度的高程點B，再過B點以β為傾角繪出水上穩(wěn)定岸坡線，與原岸坡交于點C，C點即為水上穩(wěn)定岸坡的終點。水上穩(wěn)定岸坡線起點B的高程所對應的原岸坡D點，與該線終點C之間的水平距離即為預測的塌岸寬度St。

圖1 兩段法塌岸預測示意圖

圖2 卡丘金預測法示意圖

2）卡丘金預測法

卡丘金法適用于松散沉積層，且波浪較小的水庫，如圖2。其計算公式為：

式中：St——塌岸帶最終寬度，m；

N——與土的顆粒大小有關的系數(shù)，取N=1.0；

A——庫水位變化幅度，m；

hp——波浪沖刷深度，m，取hp=1.3 m；

hb——波浪爬高，m，取hb=2.23 m；

H——正常蓄水位以上岸坡高度，m；

α——淺灘沖刷后水下穩(wěn)定坡角，（°）；

β——岸坡水上穩(wěn)定坡角，（°）；

γ——原始岸坡坡角，（°）。

3.2 參數(shù)確定

1）水下穩(wěn)定坡角α

臨汾市曲亭水庫庫區(qū)兩岸土質岸坡下部為中更新統(tǒng)洪積低液限黏土，與本庫區(qū)兩岸土質類似。據(jù)曲亭水庫庫區(qū)兩岸水下穩(wěn)定坡角調查，其左右岸水下穩(wěn)定坡角略有差異，左岸水下穩(wěn)定坡角11.0°～19.6°，平均值 14.6°，右岸水下穩(wěn)定坡角 14.1°～21.1°平均值17.3°，總平均值15.8°。水下穩(wěn)定坡角α取值參考臨汾市曲亭水庫庫區(qū)水下坡角數(shù)據(jù)，綜合考慮，α取15°。

2）水上穩(wěn)定坡角β

水上穩(wěn)定坡角β取北掌水庫庫區(qū)現(xiàn)狀兩岸土質岸天然坡角，Q3pl低液限黏土取45°，Q2pl低液限黏土取55°。

3）低液限黏土層毛細水上升高度H′

取低液限黏土室內毛細水上升高度試驗平均值H′=2.84 m。

3.3 計算成果及評價

在庫區(qū)范圍內選取三條典型剖面采用兩段法及卡丘金法進行庫區(qū)塌岸寬度及塌岸方量預測計算。計算成果見表2。

由表可知，兩段法預測庫區(qū)塌岸寬度范圍30.3～80.4 m，庫區(qū)兩岸長期塌岸總方量215.1×104m3，其中水上塌岸方量 141.3×104m3，水下塌岸方量 73.8×104m3。卡丘金預測法庫區(qū)兩岸塌岸寬度范圍33.5～81.1 m，庫區(qū)兩岸長期塌岸總方量161.1×104m3，其中水上塌岸方量 110.5×104m3，水下塌岸方量 50.6×104m3。兩種方法預測塌岸寬度范圍接近，其中卡丘金預測法估算塌岸方量要小于兩段法。上述兩種方法均假定庫區(qū)岸坡塌方土方被庫水帶至下游，直至岸坡處于穩(wěn)定狀態(tài)下所算塌岸寬度及方量。實際情況，隨著庫區(qū)兩岸塌岸進行，庫容達到一定淤積后，兩岸塌岸最終寬度及方量要小于上述計算值。

表2 庫區(qū)土質岸坡塌岸寬度及塌岸方量計算成果表

4 結論及措施建議

北掌水庫總庫容465.4萬m3，庫水回水長度約2.0 km。北掌水庫兩岸土質岸坡塌岸嚴重，預測庫區(qū)塌岸寬度范圍 30～80 m，塌岸方量 161.1～215.1 萬 m3，對庫容及大壩安全均造成威協(xié)，建議對庫區(qū)土質岸坡塌岸線范圍內的岸坡進行削坡或防護處理，邊坡坡度庫水位以上按 45°～55°考慮，庫水位以下按 13°～15°考慮。此外，大壩壩型為均質土壩，庫區(qū)兩岸第四系中、上更新統(tǒng)低液限黏土均可作為土料使用，庫區(qū)兩岸土質岸坡削坡土料可作為上壩土料使用。