青海同德抽水蓄能電站上水庫瀝青混凝土面板壩三維有限元應(yīng)力變形分析

趙智偉,劉曉光,米占寬,馬秀偉

(1. 國家能源投資集團(tuán)有限責(zé)任公司,北京 100011;2. 中國電建集團(tuán)西北勘測設(shè)計(jì)研究院有限公司,西安 710065;3. 南京水利科學(xué)研究院,南京 210029)

0 前 言

抽水蓄能電站上水庫海拔較高,水源相對匱乏[1],一般需采用全庫盆表面防滲方案[2],其中庫底防滲面積占比較大(約65%以上)、投資占比較高。西北地區(qū)地處高海拔、嚴(yán)寒、強(qiáng)紫外線、高地震烈度等嚴(yán)苛自然條件,抽水蓄能電站建設(shè)中存在混凝土低溫抗裂、材料老化、防震抗震設(shè)計(jì)復(fù)雜等突出問題,青海同德抽水蓄能電站上水庫創(chuàng)新性地采用庫岸瀝青混凝土面板+庫底土工膜防滲型式,有效地降低了工程投資,對西北地區(qū)抽水蓄能電站防滲方案設(shè)計(jì)有引領(lǐng)示范作用。

同德抽蓄上水庫工程區(qū)地層巖性復(fù)雜,地震烈度高,土石方工程量大。筑壩材料中,礫巖軟化系數(shù)低,砂巖夾板巖具有板巖含量高、埋深大、儲(chǔ)量少的特點(diǎn)。充分利用工程區(qū)開挖料,以料定壩,是土石壩設(shè)計(jì)的基本原則。為此,開展了包括靜動(dòng)力特性、濕化特性、長期流變等特性的壩料試驗(yàn)研究工作。本文針對工程區(qū)嚴(yán)寒(極端低溫-37.2℃)、強(qiáng)紫外線的惡劣工作環(huán)境,開展瀝青混凝土配合比試驗(yàn)[3]和全項(xiàng)性能試驗(yàn),針對防滲體系的薄弱環(huán)節(jié),采用三維非線性有限元計(jì)算軟件,對上水庫大壩各種運(yùn)行工況進(jìn)行了三維有限元數(shù)值模擬分析,為強(qiáng)震區(qū)百米級軟巖筑壩及庫盆防滲系統(tǒng)設(shè)計(jì)積累了寶貴經(jīng)驗(yàn)。

1 工程概況

同德抽水蓄能電站位于青海省海南藏族自治州同德縣河北鄉(xiāng)境內(nèi),電站總裝機(jī)2 400 MW,額定水頭378 m,安裝8臺300 MW的可逆式水輪發(fā)電機(jī)組。下水庫利用正在建設(shè)的瑪爾擋水電站水庫,正常蓄水位3 275.00 m,死水位3 240.00 m,調(diào)節(jié)庫容7.06億m3;上水庫為新建水庫,利用黃河右岸岸頂緩坡地形(距瑪爾擋壩址上游約11 km),挖填形成庫盆,正常蓄水位3 665.00 m,死水位3 637.00 m,調(diào)節(jié)庫容1 765萬m3。樞紐建筑物主要由上水庫、輸水系統(tǒng)、地下廠房及開關(guān)站等組成。工程屬Ⅰ等大(1)型工程,主要建筑物為1級。樞紐平面布置見圖1。

上水庫擋水壩采用瀝青混凝土面板堆石壩,壩頂高程 3 670.00 m,最大壩高112 m(壩軸線處),壩頂長2 052 m,壩頂寬10 m。上游壩坡1∶1.7,下游壩坡1∶1.7～1∶1.8,采用混凝土網(wǎng)格梁+干砌石防護(hù)。自上游至下游可分為:墊層區(qū)、過渡區(qū)、排水區(qū)、上游堆石區(qū)、下游堆石區(qū)、下游護(hù)坡、壩后任意料區(qū)。其中墊層區(qū)、過渡區(qū)水平寬度3 m,排水區(qū)水平寬度4～6 m。大壩防滲采用瀝青混凝土面板,與庫底土工膜形成全封閉的防滲體系。瀝青混凝土面板采用簡式結(jié)構(gòu),封閉層厚度2 mm,防滲層厚度10 cm,整平膠結(jié)層厚度10 cm。壩后設(shè)任意料區(qū),頂高程3 630.00 m,頂寬約154 m,下游坡比為1∶2.5。大壩典型設(shè)計(jì)斷面見圖2。

圖2 上水庫大壩典型設(shè)計(jì)斷面 單位:高程,m;其余,cm

工程區(qū)地震基本烈度為Ⅶ度,上水庫大壩抗震設(shè)防類別為甲類[4],設(shè)計(jì)地震、校核地震標(biāo)準(zhǔn)分別為基準(zhǔn)期100年超越概率2%、100年超越概率1%,相應(yīng)的基巖水平地震動(dòng)峰值加速度分別為401gal、493.3gal。參照《建筑物抗震設(shè)計(jì)規(guī)范》[5]和重大工程反應(yīng)譜的確定原則,確定的基巖水平向設(shè)計(jì)規(guī)準(zhǔn)反應(yīng)譜參數(shù)見表1。根據(jù)GB 18306-2015《中國地震動(dòng)參數(shù)區(qū)劃圖》和NB 35047-2015《水電工程水工建筑物抗震設(shè)計(jì)規(guī)范》中5.5.8條規(guī)定,確定的標(biāo)準(zhǔn)設(shè)計(jì)反應(yīng)譜參數(shù)見表2。

表1 場地基巖設(shè)計(jì)規(guī)準(zhǔn)反應(yīng)譜參數(shù)

表2 壩址基巖規(guī)范目標(biāo)反應(yīng)譜參數(shù)

根據(jù)表1、2得到的場地譜和規(guī)范譜比較見圖3,場地譜的動(dòng)力放大系數(shù)總體上大于規(guī)范譜,從偏安全角度出發(fā),抗震動(dòng)力分析采用場地譜進(jìn)行。

圖3 規(guī)范譜和場地譜對比

2 數(shù)值計(jì)算條件和模型

2.1 模型試驗(yàn)情況

2.1.1筑壩材料性能試驗(yàn)

上水庫位于極端低溫(-37.2℃)和強(qiáng)震區(qū)(設(shè)計(jì)地震0.409g、校核地震0.503g)。筑壩材料中,礫巖軟化系數(shù)低,砂巖夾板巖具有板巖含量高的特點(diǎn),筑壩材料具有抗壓強(qiáng)度低、級配寬且變異性大,濕化和流變性明顯等工程特性。為了確保工程安全和大壩料區(qū)布設(shè)合理,開展了筑壩材料性能試驗(yàn),除了常規(guī)的比重及相對密度試驗(yàn)、壓縮試驗(yàn)、滲透和滲透變形試驗(yàn)、反濾試驗(yàn)、靜動(dòng)力三軸壓縮試驗(yàn)外,增加了三軸流變和濕化變形特性試驗(yàn),試驗(yàn)的部分成果見表7、8。

2.1.2瀝青混凝土配合比試驗(yàn)

瀝青混凝土粗細(xì)骨料均為新鮮堅(jiān)硬的石灰?guī)r,填料采用磨制的石灰?guī)r礦粉,瀝青采用遼河SBS改性瀝青,質(zhì)量滿足規(guī)范要求。

(1) 瀝青混凝土防滲層

結(jié)合工程經(jīng)驗(yàn),選取最大骨料粒徑16 mm,不同油石比(7.2%、7.5%、7.8%、8.1%)、不同級配指數(shù)(0.20、0.25、0.30、0.35)和不同填料含量(11%、13%、15%)組成3組11種配合比共計(jì)32個(gè)試件,進(jìn)行瀝青混凝土防滲層配合比試驗(yàn)。根據(jù)試驗(yàn)結(jié)果,從斜坡穩(wěn)定性、防滲性(馬歇爾試件孔隙率≤2.0%)、變形性能考慮,推薦的防滲層配合比參數(shù)見表3。

表3 瀝青混凝土防滲層配合比

(2) 瀝青混凝土整平膠結(jié)層

結(jié)合工程經(jīng)驗(yàn),選取最大骨料粒徑19 mm,不同油石比(3.5%、4.0%、4.5%、5.0%)、不同級配指數(shù)(0.6、0.7、0.8、0.9)和不同填料含量(3%、5%、7%)組成3組11種配合比共計(jì)32個(gè)試件,進(jìn)行瀝青混凝土整平膠結(jié)層配合比試驗(yàn)。根據(jù)試驗(yàn)結(jié)果,從孔隙率、滲透性考慮,推薦的防滲層配合比參數(shù)見表4。

表4 整平膠結(jié)層瀝青混凝土配合比

2.1.3防滲接頭型式試驗(yàn)

土工膜與瀝青混凝土面板連接方式為在兩者之間設(shè)一道混凝土連接板,連接板與瀝青混凝土之間采用工程上已經(jīng)成熟應(yīng)用的滑動(dòng)式接頭;連接板與土工膜之間通常采用機(jī)械錨固方式連接,采用角鋼、螺栓錨固連接,土工膜與混凝土之間通過設(shè)置GB填料墊層來密封止水,土工膜與混凝土錨固連接見圖4。試驗(yàn)證明,通過對接頭部位的HDPE膜設(shè)置一個(gè)50 cm2鼓包,HDPE膜與混凝土連接板的錨固連接接頭可以適應(yīng)沉陷7.2 cm+張開4.2 cm+剪切6.6 cm的位移(接頭三向位移見第4.3節(jié),并考慮2.0的安全系數(shù)),并在0.45 MPa的水壓力下保持防滲性能。

圖4 土工膜與混凝土錨固連接

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

結(jié)合三軸試驗(yàn)結(jié)果并參考同類型工程[6-7],瀝青混凝土面板的靜動(dòng)力模型參數(shù)分別見表5、6,壩體各分區(qū)采用的靜動(dòng)力模型計(jì)算參數(shù)分別見表7、8。

表5 瀝青混凝土本構(gòu)模型參數(shù)

表6 瀝青混凝土動(dòng)力計(jì)算參數(shù)

表7 筑壩材料靜力模型計(jì)算參數(shù)

表8 筑壩材料動(dòng)力模型計(jì)算參數(shù)

2.3 施工過程模擬

根據(jù)表9所示的施工進(jìn)度安排,大壩填筑工期為35.5個(gè)月,每年12月1日至次年2月15日大壩停止填筑2.5個(gè)月,填筑工期為18個(gè)月,庫盆防滲施工工期6個(gè)月,瀝青面板施工工期3個(gè)月,第六年9月底開始下閘蓄水,假定6個(gè)月蓄水至正常蓄水位。

表9 大壩填筑歷時(shí)時(shí)間與填筑高程

表中:ρ為密度,φ0、Δφ為強(qiáng)度指標(biāo),nd、cd、Rd為體變參數(shù)。

2.4 有限元模型

根據(jù)上水庫壩體和庫盆填筑情況,三維有限元網(wǎng)格模型見圖5。三維模型共計(jì)節(jié)點(diǎn)123 268個(gè),單元128 471個(gè)。其中,土工膜為柔性材料,不能承受壓力和彎矩,只能承受拉力,并具有良好的拉伸性能,采用空間無厚度的膜單元進(jìn)行模擬。瀝青混凝土面板與周圍材料的剛度差異較大,在荷載作用下,兩者接觸面因變形不協(xié)調(diào)會(huì)發(fā)生相對位移[8],計(jì)算采用接觸摩擦單元模擬。

圖5 三維網(wǎng)格剖分

3 靜力計(jì)算成果及分析

3.1 壩體應(yīng)力變形

竣工期和水位升降期瀝青混凝土面板壩壩體最大沉降分別為46.8 cm和49.7 cm,壩體最大沉降量值約占壩高的0.44%,發(fā)生在壩軸線高程3 615.00 m處,沉降率與同類壩相比處于中等偏下水平。受壩基面總體向下游傾斜影響,竣工期壩體以向下游變形為主,下游向水平位移最大值為11.1 cm,水位升降期下游向水平位移最大值增至13.7 cm,上水庫大壩水位升降期變形等值線如圖6所示。

圖6 上水庫大壩水位升降期變形等值線 單位:cm

3.2 瀝青混凝土面板應(yīng)變與變形

瀝青混凝土面板在壩體填筑完成后一次性碾壓形成,竣工期面板變形和應(yīng)變很小。水位升降期瀝青混凝土面板撓度最大值為16.0 cm,位于反弧段高程3 636.00 m處,瀝青混凝土面板水位升降期撓度等值線如圖7所示。

圖7 瀝青混凝土面板水位升降期撓度等值線 單位:cm

水位升降期軸向和順坡向壓應(yīng)變最大值分別為0.04%和0.03%,軸向和順坡向拉應(yīng)變最大值分別為0.09%和0.10%,瀝青混凝土面板水位升降期應(yīng)變等值線如圖8所示。壓、拉應(yīng)變均在瀝青混凝土材料應(yīng)變允許范圍內(nèi)(1.0%),瀝青混凝土面板是安全的。

圖8 瀝青混凝土面板水位升降期應(yīng)變等值線 單位:%

3.3 土工膜應(yīng)變與變形

水位升降期土工膜最大主拉應(yīng)變?yōu)?.55%,發(fā)生填方區(qū)挖填分界線邊緣,拉應(yīng)變在土工膜材料應(yīng)變允許范圍內(nèi)(伸長率處于30%～100%),土工膜不會(huì)發(fā)生拉伸破壞,如圖9所示。土工膜與面板的相對錯(cuò)動(dòng)、沉陷和張開最大值分別為2.4、1.7 cm和1.6 cm,相對滑移均較小,土工膜與面板接頭部位變形是安全的。

圖9 水位升降期土工膜主拉應(yīng)變分布 單位:%

3.4 考慮流變和濕化

壩體變形主要包括加載引起的瞬時(shí)壓縮變形和堆石應(yīng)力釋放、調(diào)整和轉(zhuǎn)移導(dǎo)致的流變變形[9],同時(shí),堆石料在一定應(yīng)力條件下會(huì)浸水濕化[10],發(fā)生軟化、棱角破碎、相互滑移、填充,繼而產(chǎn)生體積縮小與剪切變形。假設(shè)面板與土工膜連接處發(fā)生破損,流變變形、濕化變形采用初應(yīng)變法計(jì)算,計(jì)算結(jié)果見圖10。蓄水運(yùn)行5年后大壩變形趨于穩(wěn)定,壩頂最終沉降29.9 cm,約為壩高的0.27%,壩體沉降最大值由49.7 cm增加至68.2 cm,為壩高的0.61%?？紤]流變后與濕化作用后,大壩壩體和防滲體系的應(yīng)變變形也有較明顯的增大,其中,瀝青面板撓度最大值由16.0 cm增加至26.5 cm,順坡向壓、拉應(yīng)變最大值分別由0.05%、0.10%增大至0.07%、0.11%;土工膜最大拉應(yīng)變由0.62%增大至0.75%,土工膜與面板的相對錯(cuò)動(dòng)、沉陷和張開最大值分別由2.4、1.7 cm和1.6 cm增加至3.0、2.2 cm和2.0 cm,但防滲體系的安全仍是有保障的。

圖10 上水庫庫1+340.85m剖面和庫1+840.00m剖面壩頂沉降過程線

4 動(dòng)力計(jì)算成果及分析

4.1 壩體動(dòng)力反應(yīng)和震陷

根據(jù)計(jì)算結(jié)果,隨著壩高的增加,大壩動(dòng)力反應(yīng)越大,壩頂15～20 m壩體的加速度放大系數(shù)明顯增加,鞭梢效應(yīng)顯著。大壩動(dòng)位移也隨著壩高的增加而增大,最大動(dòng)位移發(fā)生在壩后任意料區(qū)頂部。正常蓄水位+校核地震壩體最大震陷量為51.7 cm,震陷比約為0.46%,在土石壩震陷比統(tǒng)計(jì)范圍內(nèi)。

4.2 瀝青混凝土面板應(yīng)變與變形

正常蓄水位發(fā)生校核地震時(shí),靜動(dòng)應(yīng)變疊加后,瀝青混凝土面板軸向壓應(yīng)變和拉應(yīng)變最大值均為0.25%,順坡向壓應(yīng)變和拉應(yīng)變最大值分別為0.14%、0.25%,如圖12所示。壓、拉應(yīng)變均在瀝青混凝土材料應(yīng)變允許范圍內(nèi)(1.0%),地震作用下瀝青混凝土面板是安全的。

圖12 瀝青混凝土面板P100=1%地震工況應(yīng)變疊加等值線 單位:%

4.3 土工膜應(yīng)變與變形

校核地震過程中土工膜最大動(dòng)拉應(yīng)變?yōu)?.30%,震后主拉應(yīng)變最大值為0.70%,均發(fā)生在挖填分界線邊緣,如圖13所示。拉應(yīng)變安全系數(shù)Fab=30%/0.7%=42.9,遠(yuǎn)大于規(guī)范允許值5.0,地震工況下土工膜是安全的;土工膜與面板的切向錯(cuò)動(dòng)、沉陷和張開最大值分別為3.3、3.6 cm和2.1 cm,相對滑移均較小,土工膜與面板接頭部位變形是安全的。

圖13 震后土工膜主拉應(yīng)變分布 單位:%

在設(shè)計(jì)和校核地震下,上水庫瀝青混凝土面板堆石壩設(shè)計(jì)方案總體滿足抗震安全性要求,設(shè)計(jì)方案是可行的,滿足大壩正常運(yùn)行和抗震安全性要求。

5 結(jié) 論

基于青海同德抽水蓄能電站上水庫瀝青混凝土面板壩的三維有限元應(yīng)力變形結(jié)果分析,結(jié)合壩料試驗(yàn)、瀝青混凝土配合比試驗(yàn)和防滲接頭型式科研試驗(yàn)成果,確定了大壩及庫盆防滲體系設(shè)計(jì)的合理性。形成結(jié)論如下:

(1) 上水庫最大壩高為112 m,壩體的總體的沉降變形在允許范圍內(nèi),靜力條件下最大為49.7 cm,約為壩高的0.44%,沉降率與同類壩相比處于中等偏下水平。受壩基面總體向下游傾斜影響,壩體以向下游變形為主,下游向水平位移最大值為13.7 cm?？紤]流變與濕化作用時(shí),蓄水運(yùn)行5 a后大壩變形趨于穩(wěn)定,壩頂最終沉降29.9 cm,約為壩高的0.27%。

(2) 上水庫擋水壩應(yīng)力分布規(guī)律性較好,應(yīng)力值較小,應(yīng)力水平均不超過0.45,各工況下均沒有出現(xiàn)明顯的剪切破壞區(qū),壩體發(fā)生剪切破壞的可能性較小。

(3) 靜力條件下瀝青混凝土面板最大壓應(yīng)變和最大拉應(yīng)變分別為0.04%、0.10%,土工膜最大主拉應(yīng)變?yōu)?.55%,地震后最大拉壓應(yīng)變略有增加,遠(yuǎn)小于瀝青混凝土和土工膜允許的應(yīng)變值,應(yīng)變狀態(tài)總體良好,出現(xiàn)拉裂縫的可能性較小。

(4) 地震工況,土工膜與面板接頭的三向變位分別為3.3、3.6 cm和2.1 cm,相對滑移均較小,土工膜與面板接頭部位變形是安全的。

(5) 應(yīng)力、變形計(jì)算及防滲接頭試驗(yàn)成果說明,大壩在各工況下的變形和受力是安全的,大壩分區(qū)合理,瀝青混凝土面板和土工膜應(yīng)變在材料應(yīng)變允許范圍內(nèi),防滲體系是安全的。