不同類型降雨聯(lián)合庫水位驟降黏土心墻壩壩坡滲透穩(wěn)定數(shù)值模擬研究

陳 園

(新疆水利水電勘測設(shè)計研究院，烏魯木齊 311231)

0 引 言

降雨與庫水位變化是邊坡失穩(wěn)的重要誘因[1-4]，降雨與庫水位會使得邊坡內(nèi)部的孔壓升降[5-6]，土體的強(qiáng)度參數(shù)因遇水而降低[7-8]，有效應(yīng)力變化[9-10]，因此在最不利工況下會導(dǎo)致邊坡的失穩(wěn)災(zāi)害。對于降雨與庫水位聯(lián)合下的邊坡滲透穩(wěn)定性研究，國內(nèi)外已經(jīng)有大量的研究成果。如王樂[11]利用Geo-studio軟件對降雨與庫水位聯(lián)合下的邊坡滲透穩(wěn)定性進(jìn)行了數(shù)值模擬研究；陳玉華[12]基于重慶奉節(jié)縣某滑坡的實測資料，對降雨聯(lián)合庫水位條件下的滲透穩(wěn)定性進(jìn)行了有限差分模擬；李炎隆等[13]利用Seep/w對秭歸縣八字門滑坡庫水位聯(lián)合降雨條件下的滲流特性以及穩(wěn)定性進(jìn)行了評價。但是以上研究僅僅考慮降雨發(fā)生在庫水位驟降結(jié)束時刻或者庫水位驟降開始時刻，沒有考慮降雨發(fā)生在庫水位驟降的不同時刻。實際上，庫水位驟降持續(xù)時間較長，而降雨發(fā)生的時間持續(xù)相對較短[14]，目前還未有文獻(xiàn)進(jìn)行實際土石壩壩坡在庫水位驟降聯(lián)合降雨這種情況的相關(guān)研究。

本文利用Geo-studio軟件，以宜春市某黏土心墻壩為例，針對庫水位及不同類型降雨條件下的上下游壩坡滲流特性以及穩(wěn)定性進(jìn)行數(shù)值模擬評價，重點研究不同類型降雨發(fā)生在庫水位驟降的不同時刻，為該工程在不同工況下的滲透穩(wěn)定性規(guī)律提供直觀的認(rèn)識。

1 計算理論

1.1 非飽和理論

非飽和滲流的控制方程形式為[15]：

(1)

式中：kr為相對透水率；kij為飽和滲透張量；hc為壓力水頭；Q為源匯；C(hc)為容水度；θ為壓力水頭函數(shù)；n為孔隙率；Ss為單位貯水量。

1.2 非飽和抗剪強(qiáng)度理論

非飽和抗剪強(qiáng)度理論采Fredlund雙應(yīng)力變量公式[15]：

s=c′+σntanφ′+(ua-uw)tanφb

(2)

式中：c′與φ′為有效強(qiáng)度參數(shù)；σn為法向總應(yīng)力與孔隙氣壓力的差值；ua為孔隙空氣壓力；uw為孔隙水壓力；φb為由負(fù)孔隙水壓力而提高的強(qiáng)度。

2 計算模型與邊界條件

計算模型選擇以宜春市溫湯河四方井水利樞紐工程黏土心墻壩的典型剖面為例，四方井黏土心墻壩是一座以防洪、供水為主，兼顧發(fā)電等綜合效益的大(Ⅱ)型水利樞紐工程。壩址以上控制流域面積約173 km2，工程初擬正常蓄水位152 m(黃海高程)，死水位125 m，壩頂高程156.20 m，最大壩高51.70 m，水庫總庫容約1.189 5×108m3。計算模型及模型網(wǎng)格見圖1。

圖1 計算模型及模型網(wǎng)格Fig.1 Computational model and model grid

在圖1(a)中，為減少邊界的影響，地基厚度選取50 m，長度為300 m，設(shè)置4個監(jiān)測點，即上部1監(jiān)測點(簡稱上1)，上部2監(jiān)測點(簡稱上2)，下部1監(jiān)測點(簡稱下1)，下部2監(jiān)測點(簡稱下2)，用來監(jiān)測各個工況下壩體內(nèi)部的孔壓變化。其中，A點、B點分別位于上游壩坡的上部與下部，C點、D點分別位于下游壩坡的上部與下部。在圖1(b)中，整個模型一共劃分為4 235個單元，4 183個節(jié)點。

邊界條件如下：①bcde為庫水位邊界；②efgh為降雨入滲邊界；③hi為自由滲出邊界；④ab、ak、ijk為不透水邊界。

3 計算參數(shù)

黏土心墻壩的各分區(qū)參數(shù)根據(jù)野外勘測及室內(nèi)試驗進(jìn)行綜合確定，物理力學(xué)參數(shù)見表1，不同材料分區(qū)的土水特征曲線見圖2。

表1 壩體材料力學(xué)參數(shù)Tab.1 Mechanical parameters of dam materials

圖2 土水特征曲線Fig.2 Soil water characteristic curve

4 計算工況

為綜合研究庫水位驟降聯(lián)合降雨工況下四方井黏土心墻壩上下游壩坡的滲透穩(wěn)定特征，根據(jù)水庫多年運行經(jīng)驗，取庫水位驟降速率分別為0.5、1與1.5 m/d，結(jié)合不同類型降雨(平均型降雨、前鋒型降雨、中鋒型降雨和后鋒型降雨)進(jìn)行相應(yīng)計算，庫水位驟降耦合不同時刻不同類型降雨工況見表2，降雨時程曲線見圖3。

表2 計算工況Tab.2 Calculation conditions

圖3 降雨歷程曲線Fig.3 Rainfall history curve

5 計算結(jié)果分析

5.1 孔壓變化規(guī)律

不同監(jiān)測點的孔壓變化規(guī)律見圖4。

由圖4可見，不同監(jiān)測點的孔壓變化具有較大差異，這種差異在上游壩坡處的監(jiān)測點與下游壩坡處的監(jiān)測點的差異最為明顯。對于上游壩坡處的監(jiān)測點來說，監(jiān)測點的孔壓變化對不同類型降雨不敏感，僅上1監(jiān)測點處在不同類型降雨發(fā)生在44～54 d時孔壓出現(xiàn)小幅的上升，總體孔壓仍然呈現(xiàn)隨庫水位持續(xù)下降的趨勢；對于下游壩坡處監(jiān)測點來說，降雨發(fā)生在庫水位下降的不同時刻下孔壓均有一個大幅上升的過程，不同類型降雨使得孔壓達(dá)到最大的先后順序分別是前鋒型≥中鋒型≥平均型≥后鋒型，下游壩坡上部監(jiān)測點(下1)不同類型降雨所達(dá)到的最大孔壓幾乎一致，下部監(jiān)測點不同類型降雨所達(dá)到的最大孔壓值則不同，呈現(xiàn)前鋒型降雨≤平均型≤中鋒型≤后鋒型的規(guī)律。

5.2 安全系數(shù)變化規(guī)律

上下游壩坡安全系數(shù)的變化規(guī)律見圖5。

圖4 孔壓變化規(guī)律Fig.4 Variation of pore pressure

圖5 安全系數(shù)變化規(guī)律Fig.5 Law of variation of safety factor

對于上游壩坡來說，降雨發(fā)生在庫水位驟降時刻越后，最小安全系數(shù)越小；對于不同降雨類型來說，最小安全系數(shù)大小排序為后鋒型≤平均型≤中鋒型≤前鋒型，總體變化規(guī)律則較為類似。

對于下游壩坡來說，庫水位驟降下安全系數(shù)呈現(xiàn)一直上升的趨勢，在降雨時刻則有個突然下降的過程，不同類型降雨的峰值越前，最小安全系數(shù)出現(xiàn)時間越早。對于平均型降雨來說，降雨發(fā)生在44～54 d安全系數(shù)最??；對于前鋒型與后鋒型降雨來說，降雨發(fā)生在22～32 d安全系數(shù)最?。欢鴮τ谥袖h型降雨來說，降雨發(fā)生在0～10 d安全系數(shù)最小。

6 結(jié) 論

1) 上游壩坡處的監(jiān)測點孔壓變化對不同類型降雨不敏感，下游壩坡處監(jiān)測點不同類型降雨下孔壓變化差異較大，降雨發(fā)生在庫水位下降的不同時刻下孔壓均有一個大幅上升的過程。

2) 不同類型降雨使得孔壓達(dá)到最大的先后順序分別是前鋒型≥中鋒型≥平均型≥后鋒型，下游壩坡上部監(jiān)測點(下1)不同類型降雨所達(dá)到的最大孔壓幾乎一致，下部監(jiān)測點不同類型降雨所達(dá)到的最大孔壓值則不同，呈現(xiàn)前鋒型降雨≤平均型≤中鋒型≤后鋒型的規(guī)律。

3) 對于上游壩坡來說，降雨發(fā)生在庫水位驟降時刻越后，最小安全系數(shù)越??；對于不同降雨類型來說，最小安全系數(shù)大小排序為后鋒型≤平均型≤中鋒型≤前鋒型，總體變化規(guī)律則較為類似。

4) 對于下游壩坡來說，庫水位驟降下安全系數(shù)呈現(xiàn)一直上升的趨勢，在降雨時刻則有個突然下降的過程，不同類型降雨的峰值越前，最小安全系數(shù)出現(xiàn)時間越早。對于平均型降雨來說，降雨發(fā)生在44～54 d安全系數(shù)最??；對于前鋒型與后鋒型降雨來說，降雨發(fā)生在22～32 d安全系數(shù)最?。欢鴮τ谥袖h型降雨來說，降雨發(fā)生在0～10 d安全系數(shù)最小。