陳 園

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

0 引 言

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

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

1 計(jì)算理論

1.1 非飽和理論

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

(1)

式中：kr為相對(duì)透水率；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 計(jì)算模型與邊界條件

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

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

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

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

3 計(jì)算參數(shù)

黏土心墻壩的各分區(qū)參數(shù)根據(jù)野外勘測(cè)及室內(nèi)試驗(yàn)進(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 計(jì)算工況

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

表2 計(jì)算工況Tab.2 Calculation conditions

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

5 計(jì)算結(jié)果分析

5.1 孔壓變化規(guī)律

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

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

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

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

6 結(jié) 論

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

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

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

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