楊昕光，饒錫保，徐 晗，潘家軍，左永振

（1.中國(guó)長(zhǎng)江三峽集團(tuán)有限公司，北京市 100038；2.長(zhǎng)江科學(xué)院 水利部巖土力學(xué)與工程重點(diǎn)實(shí)驗(yàn)室，湖北省武漢市 430010）

0 引言

在水電工程建設(shè)中，由于壩址地質(zhì)條件的限制，往往需要在深厚覆蓋層等地基條件上修建高土石壩[1，2]，而且受到經(jīng)濟(jì)、環(huán)境條件約束，不宜挖除覆蓋層。如冶勒瀝青混凝土心墻壩，最大壩高124.5m，覆蓋層最大厚度達(dá)420m；瀑布溝心墻堆石壩，最大壩高186m，覆蓋層厚度為78m；長(zhǎng)河壩礫石土心墻堆石壩，最大壩高240m，覆蓋層厚度為65～76.5m。

目前，深厚覆蓋層上心墻堆石壩的防滲體系通常為心墻—廊道—防滲墻結(jié)構(gòu)，但這種傳統(tǒng)的結(jié)構(gòu)型式容易造成防滲墻壓力過(guò)大、施工周期長(zhǎng)等困難[3]?？紤]到縮短工期和方便運(yùn)行期檢修，提出了一種覆蓋層上建心墻堆石壩的新型結(jié)構(gòu)型式，即在壩基面處挖除表層部分覆蓋層，澆筑混凝土蓋板，在蓋板內(nèi)部設(shè)置廊道，以便在壩體填筑的同時(shí)進(jìn)行防滲墻的施工和運(yùn)行期的灌漿固結(jié)工作，可進(jìn)一步提高壩基淺層處理質(zhì)量與縮短工期以方便運(yùn)行期檢修，從而大為減少工程投資[3，4]。

這些建設(shè)在深厚覆蓋層上的高土石壩一般位于我國(guó)水資源豐富的西部地區(qū)，同時(shí)，這一地區(qū)為地震高烈度區(qū)。所提的新型結(jié)構(gòu)在大壩抗震安全設(shè)計(jì)等方面無(wú)成熟經(jīng)驗(yàn)可借鑒，在設(shè)計(jì)、施工和運(yùn)行中均存在極具挑戰(zhàn)性的難題。因此，有必要通過(guò)動(dòng)力有限元數(shù)值分析方法，研究深厚覆蓋層上心墻堆石壩新型結(jié)構(gòu)的動(dòng)力響應(yīng)特性，并分析其抗震安全性，以期為該類結(jié)構(gòu)堆石壩的設(shè)計(jì)與施工提供參考和依據(jù)。

1 計(jì)算模型與參數(shù)

圖1 大壩新型結(jié)構(gòu)示意圖Fig.1 Diagram for new structure of dam

研究依托我國(guó)西部擬建的某水電工程，大壩為礫石土心墻堆石壩，壩頂高程2108m，壩頂長(zhǎng)1527m，壩頂寬18m，最大壩高238m。大壩上游壩坡為1∶2.25，在高程2020m設(shè)3m寬?cǎi)R道，在高程1936m（上游圍堰頂高程）設(shè)10m寬?cǎi)R道，下游壩坡為1∶2，坡面上設(shè)“之”字形馬道，馬道寬12m，綜合坡比1∶2.2。壩體主要由防滲心墻、反濾料、過(guò)渡料和壩殼堆石體組成。心墻頂高程2106m，頂寬7m，上、下游坡1∶0.2，最低建基面高程1885m，最大底寬95.4m。河床覆蓋層為第四系沖積層，厚77～100m，主要為砂礫卵石和漂石，按地質(zhì)年代、成因類型、物質(zhì)組成及土體結(jié)構(gòu)等特征，分為4個(gè)巖組，自上而下依次為Ⅳ、Ⅲ、Ⅱ、Ⅰ巖組。在河床覆蓋層內(nèi)設(shè)兩道混凝土防滲墻，防滲墻厚1.2m，深入基巖1m?；炷辽w板位于覆蓋層上半部，其上與心墻相連，最大高度15m，寬155.4m，由覆蓋層中的防滲墻、混凝土蓋板及以上的壩體形成完整的防滲體系，大壩新型結(jié)構(gòu)示意圖如圖1所示。根據(jù)大壩設(shè)計(jì)剖面圖建立三維有限元計(jì)算模型，采用8節(jié)點(diǎn)六面體單元或者6節(jié)點(diǎn)五面體單元進(jìn)行網(wǎng)格剖分，局部采用4節(jié)點(diǎn)四面體單元，并在混凝土與土之間設(shè)置接觸面，共剖分131547個(gè)單元，有限元計(jì)算網(wǎng)格見(jiàn)圖2。

壩體和覆蓋層的靜力本構(gòu)模型采用鄧肯E-B模型[5，6]，計(jì)算參數(shù)見(jiàn)表1。動(dòng)力本構(gòu)模型采用等效線性黏彈性模型，最大動(dòng)剪切模量Gmax可以表示為：

圖2 有限元計(jì)算模型Fig.2 Finite element model

式中：σ0′為平均有效主應(yīng)力，Pa為大氣壓，Gmax、σ0′、Pa采用同一量綱；km為系數(shù)，m1為指數(shù)，由試驗(yàn)測(cè)定，具體見(jiàn)表1。通過(guò)試驗(yàn)測(cè)得的動(dòng)剪切模量比Gd/Gdmax與動(dòng)剪應(yīng)變?chǔ)茫约皠?dòng)阻尼比λ與動(dòng)剪應(yīng)變?chǔ)玫年P(guān)系見(jiàn)表2。壩體永久變形計(jì)算模型采用沈珠江計(jì)算模型[7]，具體參數(shù)見(jiàn)表3。計(jì)算中混凝土蓋板與防滲墻均采用線彈性模型，密度取2.40g/cm3，彈性模量為30GPa，泊松比為0.18。

表1 E-B模型計(jì)算參數(shù)Tab. 1 Parameters of Duncan-Chang E-B model

表3 材料永久變形模型參數(shù)Tab. 3 Parameters for permanent deformation of rockfill materials

圖3 地震波加速度時(shí)程(a)順河向;(b)豎直向;(c)壩軸向Fig.3 Time history of seismic wave acceleration

靜力計(jì)算中精確模擬了大壩的施工填筑過(guò)程和水庫(kù)水位的蓄水過(guò)程，其中壩基初始應(yīng)力場(chǎng)由自重應(yīng)力構(gòu)造。根據(jù)該大壩地形地質(zhì)條件及地震安全性評(píng)價(jià)報(bào)告，選取P100=0.02的規(guī)范譜人工波，并將其反演為基巖地震輸入波進(jìn)行動(dòng)力時(shí)程計(jì)算。通過(guò)反演后，水平向與壩軸向基巖地震波峰值為0.21g，豎直向基巖地震波峰值為0.14g。

2 大壩抗震安全性分析

2.1 大壩加速度與動(dòng)位移響應(yīng)

在P100=0.02規(guī)范譜人工波作用下，大壩最大加速度及動(dòng)位移響應(yīng)計(jì)算結(jié)果如表4所示，大壩最大橫剖面順河向最大加速度與東位移即分布圖分別如圖4和圖5所示。由計(jì)算結(jié)果可知，大壩最大順河向、壩軸向及豎直向加速度分別為0.34g、0.25g和0.23g，均位于壩頂，對(duì)應(yīng)放大系數(shù)分別為1.2、0.9和1.2。大壩順河向、壩軸向及豎直向動(dòng)位移最大值分別為13.0、8.8cm及4.3cm，且各方向的動(dòng)位移均隨著壩體高程的增加而增大，在壩頂達(dá)到最值。圖6為壩體最大剖面中軸線上絕對(duì)反應(yīng)加速度沿壩高的分布趨勢(shì)圖，其中虛線所示為壩基面與覆蓋層的交界處。由此看出，反應(yīng)加速度在覆蓋層部位有明顯的衰減現(xiàn)象。這是由于本工程壩基覆蓋層較為深厚，阻尼比較大，致使自由場(chǎng)地表（即壩體與壩基交界面處）的反應(yīng)加速度小于基巖處的加速度，因而壩體內(nèi)加速度反應(yīng)較小。

圖4 大壩最大橫剖面順河向最大加速度分布圖（單位：g）Fig.4 The maximum acceleration in transverse（unit：g）

圖5 大壩最大橫剖面順河向最大動(dòng)位移分布圖（單位：cm）Fig.5 The maximum dynamic displacement in transverse（unit：cm）

圖6 壩體中軸線處最大加速度沿壩高分布趨勢(shì)圖Fig.6 The maximum acceleration along axial line of dam

表4 大壩最大加速度與動(dòng)位移響應(yīng)Tab.4 The maximum acceleration and dynamic displacement of the dam

2.2 混凝土防滲墻與蓋板的動(dòng)應(yīng)力響應(yīng)

作為大壩新型結(jié)構(gòu)防滲體系的一部分，混凝土防滲墻與蓋板如在地震動(dòng)的作用下，其應(yīng)力超過(guò)強(qiáng)度允許值，將引起大壩防滲體系失效，危及大壩的整體安全性。因此，有必要分析在地震工況下防滲墻與蓋板的動(dòng)應(yīng)力響應(yīng)。圖7為防滲墻豎直向最大動(dòng)應(yīng)力變化云圖。由此可知，防滲墻豎直向最大動(dòng)壓應(yīng)力為1.42MPa，最大動(dòng)拉應(yīng)力為1.41MPa。從防滲墻應(yīng)力分布情況來(lái)看，地震過(guò)程中，豎直向最大瞬時(shí)動(dòng)拉應(yīng)力基本出現(xiàn)在其中部位置，該部位在震前為壓應(yīng)力區(qū)域，且混凝土動(dòng)態(tài)抗拉強(qiáng)度的標(biāo)準(zhǔn)值可較其靜態(tài)標(biāo)準(zhǔn)值提高20%[8]，因此動(dòng)靜疊加后，防滲墻拉應(yīng)力不會(huì)超過(guò)其強(qiáng)度允許值。蓋板豎直向最大動(dòng)壓應(yīng)力為2.62MPa，最大動(dòng)拉應(yīng)力為1.20MPa，限于篇幅，本文沒(méi)有給出蓋板最大動(dòng)應(yīng)力的分布圖。同理，動(dòng)靜疊加后，蓋板應(yīng)力在強(qiáng)度允許值范圍之內(nèi)，因此，分析認(rèn)為混凝土防滲墻和蓋板在地震工況下均具有足夠的安全性。

圖7 防滲墻豎直向動(dòng)應(yīng)力變化云圖(a)最大動(dòng)壓應(yīng)力;(b)最大動(dòng)拉應(yīng)力Fig.7 Contour plot for vertical dynamic stresses of diaphragm wall

圖9 壩坡抗滑穩(wěn)定安全系數(shù)時(shí)程曲線(a)上游壩坡;(b)下游壩坡Fig.9 Time history of stability safety factors of dam slopes

2.3 大壩地震永久變形

圖8為大壩地震永久變形示意圖。由圖可見(jiàn)，大壩以沉陷變形為主，且壩體輪廓是向內(nèi)部收縮，這主要是因?yàn)槎咽显诟邍鷫合录暗卣鸷奢d聯(lián)合作用下發(fā)生顆粒破碎及重組，產(chǎn)生了塑性體積變形，這符合土石壩震后變形的實(shí)測(cè)規(guī)律。根據(jù)計(jì)算結(jié)果，壩體沉陷在壩頂部靠上游面達(dá)到最大，最大值為81.2cm，約為壩高的0.34%。同時(shí)，大壩向下游水平永久變形也很顯著，最大值為49.7cm。根據(jù)土石壩實(shí)際震害調(diào)查[9，10]，一般認(rèn)為最大震陷超過(guò)壩高的1%時(shí)可產(chǎn)生明顯震害，并可能導(dǎo)致嚴(yán)重后果。因此，分析認(rèn)為上述地震不足以產(chǎn)生使壩體難以承受的變形。

圖8 大壩地震永久變形示意圖（變形放大50倍）Fig.8 The permanent deformation of the dam

2.4 大壩壩坡地震穩(wěn)定性

采用動(dòng)力有限元法[11，12]對(duì)大壩壩坡的地震穩(wěn)定性進(jìn)行了計(jì)算分析，圖9給出了上、下游壩坡穩(wěn)定系數(shù)的時(shí)程曲線。由計(jì)算結(jié)果可見(jiàn)，在地震過(guò)程中，上、下游壩坡的最小安全系數(shù)分別為1.16與1.10，最小安全系數(shù)均大于1.0?？梢?jiàn)上、下游壩坡在地震過(guò)程中是穩(wěn)定的。

3 結(jié)論

依托我國(guó)西部擬建的某水電工程，對(duì)提出的深厚覆蓋層上心墻堆石壩新型結(jié)構(gòu)進(jìn)行了有限元?jiǎng)恿Ψ治?，并?duì)其抗震安全性進(jìn)行了評(píng)價(jià)，主要得到以下幾點(diǎn)結(jié)論：

（1）大壩順河向最大順河向、壩軸向及豎直向加速度分別為0.34g、0.25g和0.23g，均位于壩頂，對(duì)應(yīng)放大系數(shù)分別為1.2、0.9和1.2；最大動(dòng)位移分別為13.0cm、8.8cm及4.3cm。

（2）地震過(guò)程中，防滲墻豎直向最大動(dòng)壓應(yīng)力為1.42MPa，最大動(dòng)拉應(yīng)力為1.41MPa；蓋板豎直向最大動(dòng)壓應(yīng)力為2.62MPa，最大動(dòng)拉應(yīng)力為1.20MPa，由于混凝土動(dòng)態(tài)強(qiáng)度較靜態(tài)強(qiáng)度有一定的提高，因此分析認(rèn)為，動(dòng)靜疊加后，防滲墻和蓋板應(yīng)力在強(qiáng)度允許值范圍之內(nèi)，具有足夠的抗震安全性。

（3）地震作用后，大壩整體發(fā)生了震陷與向下游的水平永久變形，且壩體輪廓是向內(nèi)部收縮的。由計(jì)算結(jié)果可知，壩頂最大震陷量為81.2cm，約占?jí)胃叩?.34%，設(shè)計(jì)時(shí)應(yīng)預(yù)留足夠的地震附加沉降。

（4）大壩壩坡抗震穩(wěn)定性分析表明，在地震過(guò)程中，上下游壩坡的最小安全系數(shù)均大于1.0，可認(rèn)為壩坡在地震過(guò)程中是處于穩(wěn)定狀態(tài)的。

（5）綜上，可認(rèn)為深厚覆蓋層上高心墻堆石壩新型結(jié)構(gòu)在P100=0.02的規(guī)范譜人工地震波作用下具有足夠的抗震安全性，其在受力性能、施工組織上是較為積極的處理方式，具有一定的應(yīng)用前景。

[1] 張建華，嚴(yán)軍，吳基昌，等.瀑布溝水電站樞紐工程關(guān)鍵技術(shù)綜述 [J]. 四川水力發(fā)電，2006，25（3）：8-11.ZHANG Jianhua，YAN Jun，WU Jichang，et al. Summary of optimization for Pubugou hydropower project[J]. Sichuan Water Power，2006，25（3）：8-11.

[2] 何蘊(yùn)龍，曹學(xué)興，馮蕊，等. 長(zhǎng)河壩高土質(zhì)心墻堆石壩抗震安全性研究 [J]. 武漢大學(xué)學(xué)報(bào)（工學(xué)版），2014，47（1）：1-7.HE Yunlong，CAO Xingxue，F(xiàn)ENG Rui，et al. Study of seismic safety of Changheba High rockfill dam with earth core[J].Engineering Journal of Wuhan University，2014，47（1）：1-7.

[3] 饒錫保，李會(huì)中，束一鳴，等. 超百米級(jí)覆蓋層上高土石壩壩基變形控制技術(shù)[R]. 長(zhǎng)江科學(xué)院水利部巖土力學(xué)與重點(diǎn)實(shí)驗(yàn)室，武漢，2014.

[4] 潘家軍，饒錫保，周欣華，等. 深厚覆蓋層上面板堆石壩壩新型結(jié)構(gòu)應(yīng)力變形性狀影響因素研究[J]. 水利學(xué)報(bào)，2015，46（增刊 1）：163-167.PAN Jiajun，RAO Xibao，ZHOU Xinhua，et al. Research on influence factors of stress-deformation behavior of new structure for concrete face rockfill dam on deep overburnden[J]. Journal of Hydraulic Engineering，2015，46（S1）：163-167.

[5] Duncan J M，Chang C Y. Nonlinear analysis of stress and strain in soils[J]. Journal of Soil Mechanics and Foundation Division，ASCE，1970，96（5）：1629-1653.

[6] Duncan J. M.，Byrne P. M.，Wong K. S. Strength，stress-strain and bulk modulus parameters for finite element analysis of stresses and movements in soil masses[R]. Berkeley ：Berkeley University of Calufornia，1980.

[7] 沈珠江，徐剛. 堆石料的動(dòng)力變形特性[J]. 水利水運(yùn)科學(xué)研究，1996，2：143-150.SHEN Zhujiang，XU Gang. Deformation behavior of rock materials under cyclic loading[J]. Hydro-Science and Engineering，1996，2：143-150.

[8] NB 35047—2015水電工程水工建筑物抗震設(shè)計(jì)規(guī)范[S]. 北京：中國(guó)電力出版社，2015.

[9] 顧淦臣，沈長(zhǎng)松，岑威鈞. 土石壩地震工程學(xué)[M]. 北京：中國(guó)水利水電出版社，2009.

[10] 陳生水，李國(guó)英，傅中志. 高土石壩地震安全控制標(biāo)準(zhǔn)與極限抗震能力研究 [J]. 巖土工程學(xué)報(bào)，2013，35（1）：59-65.CHEN Shengshui，LI Guoying，F(xiàn)U Zhongzhi. Safety criteria and limit resistance capacity of high earth-rock dams subjected to earthquakes[J]. Chinese Journal of Geotechnical Engineering，2013，35（1）：59-65.

[11] 劉漢龍，費(fèi)康，高玉峰. 邊坡地震穩(wěn)定性時(shí)程分析方法[J]. 巖土力學(xué)，2003，24（4）：553-556.LIU Hanlong，F(xiàn)EI Kang，GAO Yufeng. Time history analysis method of slope seismic stability[J]. Roch and Soil Mechanics，2003，24（4）：553-556.

[12] 趙劍明，常亞屏，陳寧. 高心墻堆石壩地震變形與穩(wěn)定分析[J].巖土力學(xué)，2004，25（S2）：423-428.ZHAO Jianming，CHANG Yaping，CHEN Ning. Study on earthquake-induced permanent deformation and dynamic stability of high core rockfill dam[J]. Roch and Soil Mechanics，2004，25（S2）：423-428.