考慮非飽和滲流的谷幅變形對(duì)高拱壩影響分析

辛長虹，趙 引

(河海大學(xué) 力學(xué)與材料學(xué)院，江蘇 南京 210098)

近年來，小灣、溪洛渡、錦屏一級(jí)、烏東德等一批300 m 級(jí)特高拱壩在我國西南地區(qū)建成并先后進(jìn)入蓄水、運(yùn)行階段。我國西南地區(qū)地形地質(zhì)條件復(fù)雜，水庫蓄水后大多出現(xiàn)谷幅異常變形現(xiàn)象，即蓄水后兩岸邊坡向河谷中心收縮。相關(guān)文獻(xiàn)資料[1-3]顯示，在蓄水初期或蓄水完成后，溪洛渡拱壩谷幅收縮高達(dá)87 mm，錦屏一級(jí)拱壩約為30 mm，李家峽拱壩為15～30 mm，二灘拱壩約為3.5 mm。加拿大Oldman 拱壩、意大利Vajont 拱壩、日本黑部拱壩等在蓄水后也出現(xiàn)了不同程度的谷幅收縮現(xiàn)象。

拱壩是高次超靜定結(jié)構(gòu)，對(duì)地基變形尤其是不均勻變形極為敏感[4]。在谷幅收縮作用下，壩體可能產(chǎn)生局部開裂，形成滲漏通道，影響拱壩正常運(yùn)行，甚至威脅大壩安全[5]。如瑞士Zeuzier 拱壩上游壩面開裂、意大利Beauregard 拱壩下游壩面開裂、我國溪洛渡拱壩下游水墊塘開裂、二灘拱壩下游壩面開裂等都是蓄水期谷幅收縮引起壩體破壞的典型案例。因此，有必要對(duì)高拱壩蓄水期產(chǎn)生的谷幅變形及影響進(jìn)行深入研究。

針對(duì)蓄水誘發(fā)的谷幅變形規(guī)律及其相關(guān)機(jī)制，許多學(xué)者已開展了相關(guān)研究。楊學(xué)超等[3,6-7]依據(jù)邊坡監(jiān)測數(shù)據(jù)對(duì)谷幅長期變形進(jìn)行預(yù)測。楊杰等[1,8-9]針對(duì)谷幅收縮現(xiàn)象進(jìn)行了分析。張國新等[10]基于溪洛渡拱壩邊坡監(jiān)測數(shù)據(jù)的回歸分析，指出谷幅變形主要是由近壩區(qū)水文地質(zhì)條件改變所觸發(fā)的岸坡蠕變變形形成。湯雪娟等[11]采用增量荷載法研究了穩(wěn)定滲流場對(duì)溪洛渡拱壩地基變形的影響，表明滲流場是引起谷幅收縮等現(xiàn)象的重要影響因素之一。楊強(qiáng)等[12]提出了非飽和裂隙巖體的有效應(yīng)力原理，并以Drucker-Prage準(zhǔn)則的形式應(yīng)用于錦屏一級(jí)拱壩中，發(fā)現(xiàn)庫盆變形的計(jì)算結(jié)果與實(shí)際監(jiān)測值基本一致。任青文等[13]推導(dǎo)了巖體在非飽和滲流-應(yīng)力耦合作用下的控制方程，并應(yīng)用于圓柱試件和邊坡入滲過程的變形分析，認(rèn)為飽和度的改變是引起巖體變形的重要因素之一。但目前的研究對(duì)谷幅變形的產(chǎn)生機(jī)制還沒有定論，也缺乏定量的分析；谷幅變形對(duì)特高拱壩的應(yīng)力變形影響也研究得較少。

相較于飽和滲流分析，非飽和滲流更適合裂隙巖體入滲過程，更能反映蓄水期壩體及山體異常變形現(xiàn)象。本文結(jié)合錦屏一級(jí)拱壩谷幅變形的監(jiān)測資料，利用有限元分析方法，進(jìn)行非飽和滲流過程中谷幅變形情況及對(duì)高拱壩位移和應(yīng)力狀態(tài)影響的研究。

1 非飽和滲流基本理論

1.1 非飽和滲流定律

在非飽和滲流場分析過程中，孔隙介質(zhì)的滲透系數(shù)與飽和度相關(guān)，而飽和度通常被定義成基質(zhì)吸力的函數(shù)。基質(zhì)吸力即為孔隙介質(zhì)的毛細(xì)管力，可由負(fù)孔壓來表示，其與飽和度的函數(shù)關(guān)系稱為吸濕曲線，表達(dá)式[13]如下：

飽和度對(duì)滲透系數(shù)的影響可通過折減系數(shù)ks來表示，折減系數(shù)和飽和度的關(guān)系[13]可表示為：

式中：s為飽和度；pc為基質(zhì)吸力；a、m為曲線參數(shù)，對(duì)裂隙巖體而言，可取a=1.0 MPa，0

圖1 折減系數(shù)與飽和度的關(guān)系Fig.1 Relationship between reduction coefficient and saturation

圖2 飽和度與基質(zhì)吸力的關(guān)系Fig.2 Relationship between saturation and matric suction

因此，孔隙介質(zhì)非飽和狀態(tài)下的滲透系數(shù)和達(dá)西定律可表述為：

式中：k、分別為飽和、非飽和狀態(tài)下的滲透系數(shù)；→v為矢量滲流速度；J為水力梯度；H為孔隙介質(zhì)水頭；?為拉普拉斯算子。

1.2 非飽和滲流數(shù)學(xué)模型

孔隙介質(zhì)中流體運(yùn)動(dòng)的微分方程可根據(jù)達(dá)西定律和質(zhì)量守恒定律進(jìn)行推導(dǎo)。根據(jù)質(zhì)量守恒定律，非飽和孔隙介質(zhì)的連續(xù)性方程[14]可表示為：

式中：vx、vy、vz分別為x、y、z方向的滲流速度；ρ為水的密度；n為孔隙比。

當(dāng)不考慮孔隙介質(zhì)中水的壓縮性時(shí)，體積變化量可表示為：

聯(lián)合達(dá)西定律，并忽略位置水頭的影響，可得到非飽和滲流微分方程：

式中：?v為孔隙介質(zhì)的體積應(yīng)變率；uw為孔隙水壓。

利用非飽和滲流微分方程，結(jié)合給定的初始條件和邊界條件，即可確定非飽和滲流場的分布。

1.3 Drucker-Prage 準(zhǔn)則

在進(jìn)行彈塑性計(jì)算時(shí)，脆性材料的屈服準(zhǔn)則采用Drucker-Prage 系列準(zhǔn)則中的DP1 準(zhǔn)則（M-C 外角點(diǎn)外接圓），表達(dá)式[15]如下：

式中：α、K為關(guān)于黏聚力c和內(nèi)摩擦角φ的參數(shù)；I1為應(yīng)力張量的第一不變量；J2為應(yīng)力偏張量的第二不變量。

2 錦屏一級(jí)高拱壩有限元模型

2.1 三維有限元網(wǎng)格

錦屏一級(jí)拱壩有限元模型的坐標(biāo)系為：x軸正向指向左岸，y軸正向指向上游，z軸正向豎直向上，頂拱上游面與拱冠梁的交點(diǎn)為坐標(biāo)原點(diǎn)O，如圖3 所示。模型的橫河向范圍：自拱壩中心線向左右兩岸各延伸2 倍壩高，共寬1 220 m；順河向范圍：自壩軸線向上游延伸1 倍壩高，向下游延伸2 倍壩高，共長915 m；鉛直向范圍：▽1 450 m～▽2 100 m，其中，壩底▽1 580 m，壩頂▽1 885 m。壩體沿厚度剖分8 層，沿高度剖分26 層。有限元網(wǎng)格全部采用等參八節(jié)點(diǎn)六面體和六節(jié)點(diǎn)五面體單元，節(jié)點(diǎn)總數(shù)66 689，單元總數(shù)61 620，其中壩體單元數(shù)6 005?；灸M了大壩的結(jié)構(gòu)特征、地形地貌、巖體構(gòu)造（包括左岸斷層F5、F8、F42-9、煌斑巖脈FX 和右岸斷層F13、F14）、壩肩墊座及防滲帷幕。斷層及巖脈的分布位置如圖4 所示，防滲帷幕布置情況如圖5 所示。

圖3 錦屏一級(jí)拱壩有限元網(wǎng)格Fig.3 Finite element mesh of Jinping Ⅰ arch dam

圖4 錦屏一級(jí)拱壩斷層及巖脈分布示意Fig.4 Distribution of faults and dikes of dam

圖5 錦屏一級(jí)拱壩防滲帷幕布置Fig.5 Distribution of anti-seepage curtain of dam

2.2 材料參數(shù)

錦屏一級(jí)拱壩彈塑性計(jì)算使用的各類材料參數(shù)如表1 所示。在進(jìn)行滲流場分析時(shí)，不考慮壩體、墊座、建基面和防滲帷幕等混凝土材料的滲透性。

表1 材料的物理力學(xué)參數(shù)Tab.1 Physical and mechanical parameters of the material

2.3 計(jì)算工況

考慮到谷幅變形的時(shí)效特點(diǎn)，按不考慮滲流場（工況1）、考慮非飽和滲流場（工況2）及考慮飽和滲流場（工況3）共3 種工況進(jìn)行研究。同時(shí)考慮到蓄水期水位的變化，結(jié)合錦屏的4 個(gè)蓄水階段，工況2 又分為2-H1、2-H2、2-H3 和2-H4 共4 種情況。具體計(jì)算工況及荷載組合見表2。

表2 計(jì)算工況及荷載組合Tab.2 Calculation cases and load combinations

3 非飽和滲流場作用的谷幅變形分析

3.1 吸濕曲線對(duì)谷幅變形的影響

為能較好反映非飽和滲流場引起的谷幅變形情況，采用最高滲流水位1 880 m，借助圖6 中的3 個(gè)剖面（y軸坐標(biāo)依次為50、?150、?300 m），取各剖面位置河谷面上從右岸到左岸的34 個(gè)節(jié)點(diǎn)，來分析谷幅變形對(duì)吸濕曲線的敏感性。河谷面節(jié)點(diǎn)分布如圖7 所示。

圖6 錦屏一級(jí)拱壩邊坡剖面位置分布Fig.6 Distribution of slope profile of dam

圖7 河谷面節(jié)點(diǎn)分布Fig.7 Distribution of nodes on valley surface

圖8 為采用不同吸濕曲線時(shí)非飽和滲流場引起的各剖面谷幅變形。由圖8 可知，在橫河向?90～30 m的范圍內(nèi)，參數(shù)m的取值對(duì)谷幅變形的影響很小，可忽略不計(jì)；但該范圍以外，同一位置的谷幅變形隨參數(shù)m的增大明顯減小。不同曲線參數(shù)下，谷幅變形基本呈收縮狀態(tài)，但當(dāng)參數(shù)m超過0.6 時(shí)，在坡頂出現(xiàn)一定程度的谷幅擴(kuò)張變形，且參數(shù)m越大擴(kuò)張現(xiàn)象越明顯。由此可見，吸濕曲線的參數(shù)對(duì)谷幅變形影響較大，應(yīng)該慎重選擇。

圖8 吸濕曲線對(duì)谷幅變形的影響Fig.8 Influence of absorption curve on valley width deformation

3.2 蓄水期谷幅變形規(guī)律

根據(jù)錦屏一級(jí)拱壩壩址區(qū)的邊坡地質(zhì)條件與溫度場變化情況，滲流場是影響其谷幅變形的主要因素[12]。在3.1 節(jié)谷幅變形對(duì)吸濕曲線敏感性分析的基礎(chǔ)上，將錦屏一級(jí)拱壩典型谷幅測線的橫河向位移計(jì)算值與蓄水期監(jiān)測數(shù)據(jù)[16]進(jìn)行對(duì)比，如表3 所示?？梢姡琺=0.5 的吸濕曲線更適合錦屏一級(jí)拱壩邊坡巖體的非飽和滲流場分析。

表3 典型谷幅測線位移監(jiān)測值與計(jì)算值對(duì)比Tab.3 Monitoring and calculated transverse displacement of valley width measuring lines 單位：cm

在m=0.5 的基礎(chǔ)上，分別計(jì)算工況2-H1、工況2-H2、工況2-H3、工況2-H4 的滲流場，研究非飽和滲流過程中谷幅變形規(guī)律。圖9 為非飽和滲流過程中不同水位下各剖面的谷幅變形。由圖9 可知，在非飽和滲流場作用下，兩岸邊坡向河谷變形，且上游谷幅變形相對(duì)下游偏大，壩基處谷幅變形介于上下游變形之間。在各剖面，各谷幅變形曲線的變化趨勢基本一致，極值點(diǎn)大致相同；在同一位置，谷幅變形值隨滲流水位的升高而增大；在極值點(diǎn)處，不同水位間的谷幅變形差值最大，向兩岸延伸差值逐漸減小。最大谷幅變形發(fā)生在Ⅰ-Ⅰ剖面的▽1 700 m 高程，兩岸相對(duì)收縮值為6.79 cm。此外，Ⅱ-Ⅱ剖面穿過左岸壩肩混凝土墊座，墊座對(duì)谷幅變形起到一定的抑制作用，且水位越高影響效果越明顯。

圖9 滲流水位對(duì)谷幅變形的影響Fig.9 Influence of seepage water level on valley width deformation

4 谷幅變形對(duì)拱壩的影響

4.1 對(duì)壩體位移的影響

根據(jù)谷幅變形計(jì)算結(jié)果，飽和滲流場引起的谷幅變形大于非飽和滲流場的變形。而在非飽和滲流過程中，壩體位移和應(yīng)力的分布規(guī)律基本保持不變，但量值略有差別。隨著滲流水位的升高，壩體最大順河向位移和最大主壓應(yīng)力逐漸減小，最大主拉應(yīng)力逐漸增大。因此，以工況1、2-H4 和3 為例，對(duì)比分析滲流場引起的谷幅變形對(duì)拱壩位移和應(yīng)力的影響。如圖10～11 所示，考慮谷幅變形后，拱壩的位移分布規(guī)律與正常工況（工況1）基本一致，但橫河向位移值略有增加，順河向變形位移值略有減小。隨谷幅變形的增加，壩體最大橫河向位移逐漸增大，且位置逐漸向兩岸壩肩移動(dòng)；最大順河向位移仍然出現(xiàn)在拱冠梁下游面，但量值從正常工況的7.88 cm 減小到考慮非飽和滲流場的6.96 cm 和考慮飽和滲流場的5.39 cm。由此可見，谷幅變形對(duì)壩體產(chǎn)生擠壓作用，導(dǎo)致壩體向上游變形，且谷幅變形越大，壩體向上游變形的趨勢越明顯。

圖10 壩體橫河向位移u1Fig.10 Comparison of displacement of dam in x direction

圖11 壩體順河向位移u2Fig.11 Comparison of displacement of dam in y direction

此外，谷幅變形對(duì)壩基位移的影響也比較明顯?，F(xiàn)以下游壩面為例，統(tǒng)計(jì)壩基位移變化規(guī)律。由圖12可知，考慮谷幅變形后，壩基橫河向變形依舊成反對(duì)稱分布，但變形方向與正常工況相反，且位移峰值從壩中移到了壩頂；壩基順河向位移均有所增加，且左岸變化更為明顯，這可能與左岸不良地質(zhì)條件有關(guān)。相比于非飽和滲流場，飽和滲流場引起的壩基橫河向位移偏大，且順河向變形偏小。

圖12 壩基位移分布Fig.12 Comparison of displacements of dam foundation

4.2 對(duì)壩體應(yīng)力的影響

圖13 給出了工況1、工況2-H4 和工況3 的壩體最大主拉應(yīng)力分布。由圖13 可知，考慮谷幅變形后，拱壩最大主拉應(yīng)力仍然出現(xiàn)在上游面，但位置由壩踵轉(zhuǎn)移到左岸壩肩，且量值也從正常工況的2.15 MPa 減小到考慮非飽和滲流場的1.50 MPa 和考慮飽和滲流場的1.54 MPa。可見飽和滲流場產(chǎn)生的最大主拉應(yīng)力略大于非飽和滲流場的。

圖13 壩體最大主拉應(yīng)力σ1Fig.13 Comparison of major principal stress of dam

圖14 對(duì)比了不考慮和考慮谷幅變形時(shí)，壩體最大主壓應(yīng)力分布情況。由圖14 可知，考慮滲流場引起的谷幅變形后，壩體最大主壓應(yīng)力的位置基本固定在左岸壩肩的下游面，但量值從正常工況的18.10 MPa減小到考慮非飽和滲流場的17.64 MPa 和考慮飽和滲流場的17.02 MPa；下游壩面的高壓應(yīng)力區(qū)開始向拱冠梁中部轉(zhuǎn)移，且谷幅變形越大，主壓應(yīng)力減小和轉(zhuǎn)移的效果越明顯。表4 為各工況的壩體應(yīng)力計(jì)算結(jié)果，扣除建基面應(yīng)力集中，其大小基本滿足規(guī)范要求。

表4 壩體應(yīng)力計(jì)算結(jié)果Tab.4 Calculation results of stress of dam

圖14 壩體最大主壓應(yīng)力σ3Fig.14 Comparison of minor principal stress of dam

5 結(jié)語

針對(duì)特高拱壩在蓄水初期出現(xiàn)的谷幅異常變形情況，以錦屏一級(jí)拱壩為例，通過對(duì)大壩進(jìn)行精細(xì)化建模，從滲流場的角度研究了蓄水期谷幅變形規(guī)律，分析了滲流場產(chǎn)生的谷幅收縮對(duì)壩體位移和應(yīng)力的影響，得到以下結(jié)論：

（1）谷幅變形對(duì)吸濕曲線比較敏感，吸濕曲線關(guān)鍵參數(shù)m越大，谷幅變形值越小。當(dāng)m=0.5 時(shí)，非飽和滲流場產(chǎn)生的谷幅變形更接近實(shí)際監(jiān)測值。

（2）在非飽和滲流場作用下，兩岸邊坡向河谷中心變形，且上游比下游的谷幅變形值大。隨著水位升高，谷幅變形值不斷增大，且飽和滲流場作用下的谷幅變形大于非飽和滲流場的。

（3）在非飽和滲流過程中壩體位移和應(yīng)力的分布規(guī)律基本保持一致，但隨滲流水位的升高，壩體最大順河向位移和最大主壓應(yīng)力逐漸減小，最大主拉應(yīng)力逐漸增加。

（4）滲流場引起的谷幅收縮對(duì)壩體產(chǎn)生擠壓作用，導(dǎo)致壩體最大順河向位移減小。谷幅變形越大，壩體向上游變形的量值越大。

（5）滲流場引起的谷幅變形導(dǎo)致上游壩面高拉應(yīng)力區(qū)向頂拱的左右拱端轉(zhuǎn)移，下游壩面高壓應(yīng)力區(qū)向拱冠梁中部轉(zhuǎn)移，但對(duì)應(yīng)力大小影響有限，不會(huì)威脅拱壩安全。