李 影

(營口市河務(wù)管理處,遼寧 營口 115000)

錦州市市區(qū)供水水源大部分來源于地下水,但是隨著經(jīng)濟建設(shè)的迅猛發(fā)展,水資源超采、污染現(xiàn)象嚴(yán)重。根據(jù)國家、省、市相關(guān)要求,目前錦州市的地下水資源已不適合繼續(xù)開采,因此,建設(shè)地表水資源工程迫在眉睫。但是,小凌河干流上缺少大型控制性水利工程。通過對小凌河進行現(xiàn)場勘測,該地區(qū)地形、水文、地質(zhì)等條件滿足修建大型水利樞紐工程[1]。此外,錦凌水庫建成后,對于解決錦州市市區(qū)缺水問題具有較大的作用。因此,興建錦凌水庫樞紐工程,深入開發(fā)小凌河水資源,對于解決錦州市市區(qū)嚴(yán)重缺水問題至關(guān)重要。

大型土石壩設(shè)計歷史悠久,研究方法較多,在本構(gòu)模型構(gòu)建方面：Duncan.J.M[2]通過運用有限元方法確定了能夠適當(dāng)選取材料的本構(gòu)模型；顧淦臣,郭軍,朱亞林等[3]利用動力彈塑性分析方法,揭示了地震對壩體切應(yīng)變、變形的影響規(guī)律,包承綱[4]等引進了土動力學(xué)理論,研究了長江三峽區(qū)域厚細(xì)砂沉積層的穩(wěn)定性和整體受力問題；陳群等[5]通過引進并應(yīng)用飽和——非飽和滲流理論,分析了青海溝后大壩蓄水后庫水滲透問題,解決了非飽和滲流時土壩的穩(wěn)定問題；劉先珊,劉少煒,余永康[6]研究了包含防滲墻的土石壩體滲流情況,確定汛期壩體滲流場,得到滲透坡降與材料滲透系數(shù)的關(guān)系。在壩體結(jié)構(gòu)形式和強度方面：李星[7]通過分析水泥土攪拌樁擋墻強度和剛度變化,確定了其對基坑變形的單一影響因素。在壩體結(jié)構(gòu)數(shù)值分析方面：劉奉銀,劉國棟等[8]采用數(shù)值分析方法,研究了防滲墻深度與壩體滲流特性的關(guān)系；劉建軍等[9]在巖土飽和——非飽和滲流理論基礎(chǔ)上,對徑流水滲流進行了動態(tài)模擬,揭示了地下水水頭、壓力、水流流速等的變化規(guī)律；彭華[10]在這一研究方法和思路上進行了改進和優(yōu)化,揭示了迭代步間容水度變化曲線關(guān)系。

綜上可知,目前對于水利工程土石壩等的研究較多,對于考慮滲流穩(wěn)定性作用下土石壩體、土石混凝土混合壩體的理論、設(shè)計和施工技術(shù)等的研究理論也較為成熟[13]。但是,對于小凌河復(fù)雜地質(zhì)條件下土石壩施工及其穩(wěn)定性的研究目前尚存在不足和值得商榷之處,特別是壩體整體的應(yīng)力位移變化數(shù)值分析研究較少。

本文基于錦凌水庫小凌河土壩設(shè)計和施工技術(shù)、防滲防護技術(shù)和穩(wěn)定性分析研究以及目前尚存在的不足之處,對錦凌水庫擋水建筑物結(jié)構(gòu)設(shè)計和施工進行了研究,并通過有限元模擬計算對結(jié)構(gòu)內(nèi)部應(yīng)力、變形等進行了分析計算,通過定性和定量相結(jié)合的方法分析驗證了該土壩壩體的穩(wěn)定性。本文的研究結(jié)論可為土壩結(jié)構(gòu)設(shè)計和安全施工,保證壩體穩(wěn)定性和正常運行等提供工程指導(dǎo)和借鑒。

表2 結(jié)構(gòu)設(shè)計參數(shù)確定

1 水庫概況及其防汛標(biāo)準(zhǔn)

錦凌水庫位于小凌河要道,水庫容水量8.08×108m3,攔河壩頂最高為48.29m,具有用水、防洪、蓄水、保證用水安全等作用,水庫下游方向9.5km處為錦州市。根據(jù)水庫等級劃分：工程為大(2)型II等,永不撤除類關(guān)鍵構(gòu)筑物為2級；普通構(gòu)筑物(如下游右岸防沖導(dǎo)墻)為3級；可撤除類構(gòu)筑物則為4級。防汛標(biāo)準(zhǔn)見表1。

表1 主要建筑物的防訊標(biāo)準(zhǔn)

2 土壩壩體結(jié)構(gòu)設(shè)計及其穩(wěn)定性計算

2.1 壩體結(jié)構(gòu)設(shè)計

土壩分左、右岸兩部分,左岸土壩樁號為0+059.00m～0+558.00m,長度499.0m；右岸土壩樁號為0+855.50m～1+207.00m,長度351.5m。左、右岸土壩總長850.5m。具體壩體結(jié)構(gòu)設(shè)計見表2。

2.2 土壩穩(wěn)定性計算

2.2.1 滲流穩(wěn)定計算

根據(jù)設(shè)計規(guī)范選擇粘土心墻結(jié)構(gòu)作為土壩的防滲體結(jié)構(gòu),壩體下游選擇表層排水。按水庫洪水調(diào)節(jié)運行方式,滲流計算水位組合見表3。

表3 滲流計算水位組合表 單位:m

根據(jù)拉普拉斯方程：

(1)

式中，H—滲流平面域內(nèi)水頭；Kx、Ky—土體X、Y向滲透系數(shù)。通過計算,土壩滲流穩(wěn)定分析計算成果見表4。計算成果如圖1～3所示。

由表3可以得出：壩體滲透量小,心墻滲透坡降小于允許滲透坡降,符合規(guī)范要求。但是,壩腳滲透坡降大于允許值,不符合要求,需要加設(shè)排水體。排水體采用表面式貼坡排水,排水體尺寸高×寬為39.00m×2.0m。

圖4 典型剖面計算簡圖及滑弧位置

圖1 正常蓄水位滲流

圖2 設(shè)計水位滲流

圖3 校核水位滲流

2.2.2 壩體邊坡穩(wěn)定計算

壩體邊坡的穩(wěn)定應(yīng)根據(jù)SL274- 2001《碾壓式土石壩設(shè)計規(guī)范》的規(guī)定進行復(fù)核驗算,正常運用條件：K=1.35；非正常運用條件I∶K=1.25。水庫水位降落期的安全系數(shù)計算公式為：

式中，R—圓弧半徑;u—孔隙水壓力(土條底面)；γw—水的容量；W—條塊重量；α—通過條塊底面中心點的直線與此條塊的重力線之間的夾角;b—條塊寬度；Z—條塊底面中心與水庫存的距離(水庫位于條塊底面上方)；Q—水平地震慣性力(上負(fù)，下正)；c′、φ′—有效強度指標(biāo)；Mc—力矩(水平向地震慣性力對圓心)；V—垂直地震慣性力(上負(fù)，下正)。

壩體穩(wěn)定計算參數(shù)見表5。

表5 壩體穩(wěn)定計算參數(shù)表

壩坡抗滑穩(wěn)定分析計算結(jié)果見表6,典型剖面計算簡圖及滑弧位置如圖4所示。

由表5土壩壩坡穩(wěn)定分析結(jié)果可知：在各組合荷載作用時,土壩穩(wěn)定性滿足規(guī)范要求。

3 應(yīng)力及變形計算

3.1 沉降計算

粘土心墻竣工及最終的沉降量計算：

(3)

表6 土壩壩坡穩(wěn)定分析結(jié)果

壩基最終沉降量計算：

(4)

式中，S∞—最終沉降量(壩基),m；hi—第i層土層厚度,m;n—土層分層數(shù)；Ei—第i計算土層變形模量,kN；ei0—第i層的起始孔隙比；Si—沉降量,m;Pi—第i計算土層由壩體荷載產(chǎn)生的豎向應(yīng)力,kN；eit—孔隙比(第i層相應(yīng)于竣工時或最終豎向的的有效應(yīng)力作用下)。

對壩體部分的沉降量計算時,將壩體沿高度分成幾個土層,每一土層厚度以不大于1/5～1/10壩高為宜。本設(shè)計心墻高度為34.9m,每層厚度均為5.0m。在計算壩基部分的沉降量時,對于均質(zhì)壩基,壩底寬度的1/4作為每一層厚度的最大值,本設(shè)計壩底寬度約190m,每一層厚度約47.5m,但實際壩基最大覆蓋層厚度為13.5m,采用實際厚度13.5m為分層厚度。經(jīng)計算,壩體總沉降量為904.1mm,施工期沉降量為673.5mm,壩基沉降量為165.5mm,則壩體最終沉降量為904.1-673.5+165.5=396.1mm。

3.2 應(yīng)力計算

根據(jù)規(guī)范要求,壩體豎向應(yīng)力為計算點以上土柱重量。通過計算,蓄水期壩體豎向應(yīng)力最大值為0.334MPa,竣工期應(yīng)力最大值為0.571MPa,蓄水期壩體豎向應(yīng)力最大值為2.3MPa,竣工期最大值為0.833MPa。

4 有限元計算及其結(jié)果分析

選取剖面0+400.00為典型剖面進行計算,計算采用總應(yīng)力法,將心墻、壩殼、防滲墻和砂礫地基視為整體結(jié)構(gòu),分別對壩基、壩體進行應(yīng)變及應(yīng)力計算。

單元體內(nèi)任一點的應(yīng)變和應(yīng)力均由三個分量來確定：

(5)

式中，εx—沿x坐標(biāo)方向的應(yīng)變；εy—沿x坐標(biāo)方向的應(yīng)變；μ—泊松比;γxy—沿y坐標(biāo)方向的應(yīng)變;E—彈性模量,MPa。

4.1 計算條件

材料的計算參數(shù)見表7。

表7 計算參數(shù)表

4.2 位移和應(yīng)力計算結(jié)果

由于受到壩體自重的影響,壩體沿著心墻及防滲墻為分界線,向上、下游產(chǎn)生較大變形,圖5和圖6分別為壩體竣工期間的豎向、水平位移云圖,圖7、圖8分別為壩體蓄水期間的豎向、水平位移云圖。

圖5 土壩豎向位移云圖

圖6 土壩水平位移云圖

圖7 土壩豎向位移云圖

圖8 土壩水平位移云圖

4.2.1 壩體位移和應(yīng)力

從以上圖中可看出上、下游水平位移最大值為14.0mm、13.0mm。心墻部位產(chǎn)生豎向最大位移,為214.0mm,壩高的0.65%為垂直沉降量。研究壩體應(yīng)力分布規(guī)律,可以得出其應(yīng)力分布從壩基向壩坡逐漸減小,從心墻向兩側(cè)逐漸減小,其等值線與壩坡平行。

竣工期,壩體水平方向大應(yīng)力為0.875MPa、小應(yīng)力別0.93MPa。壩體豎直方向應(yīng)力為0.597MPa。蓄水期,壩體水平方向大、小應(yīng)力分別為0.633MPa和0.082MPa,豎直方向應(yīng)力較竣工期小,為0.356MPa。

通過對比發(fā)現(xiàn)：壩體竣工期大主應(yīng)力與垂直應(yīng)力基本相等,同時也適用于蓄水期壩體盈利變化,并與壩體的上覆自重較為接近。此外,小主應(yīng)力逐漸成水平方向,這主要是因為壩體結(jié)構(gòu)受力均衡的結(jié)果。

4.2.2 粘土心墻位移

粘土心墻竣工期的水平位移較小,最大值為13.0mm。蓄水期,粘土心墻水平位移最大值增至30.0mm,位于心墻頂部,其值在通常范圍內(nèi)是允許的?？⒐て谛膲Υ怪蔽灰谱畲笾禐?12.0mm,蓄水期為218.0mm,均發(fā)生在心墻的中上部。粘土心墻水平位移與基礎(chǔ)防滲墻水平位移在沿高度方向上能夠表現(xiàn)光滑連續(xù),通過這一點可以得出接頭部位不會產(chǎn)生剪切破壞和錯動。

4.2.3 混凝土防滲墻應(yīng)力

竣工期和蓄水期混凝土防滲墻豎直方向應(yīng)力最大值、水平方向應(yīng)力最大值為0.8MPa、0.6MPa。通過對兩種狀態(tài)下的應(yīng)力最大值比較分析,蓄水期的心墻各點靜水壓力小于該處的垂直應(yīng)力,該位置將不會產(chǎn)生水力劈裂、破壞現(xiàn)象。

壩體工程施工期間,必須對結(jié)構(gòu)的應(yīng)力應(yīng)變、較大變形及其位移變化、原型觀測及數(shù)據(jù)等深入分析并及時修正,并根據(jù)計算結(jié)果和參數(shù)調(diào)整不斷進行校核,以確保成果的準(zhǔn)確性和合理性。同時,可以作為工程修正設(shè)計和制定工程措施的重要依據(jù),保障工程設(shè)計、施工和運行管理的順利進行。

5 結(jié)語

(1)通過對土壩壩體結(jié)構(gòu)設(shè)計進行研究,確定了主要建筑物12項設(shè)計項目研究內(nèi)容及其防汛標(biāo)準(zhǔn),設(shè)計了土壩擋水建筑物的設(shè)計參數(shù),通過構(gòu)建正常蓄水位滲流、設(shè)計水位滲流和校核水位滲流計算模型,對水位降落期的安全系數(shù)計算公式進行了研究,確定并計算了壩體的穩(wěn)定性計算參數(shù)及其施工條件下的安全系數(shù)公式,通過各組合荷載相互作用,土壩穩(wěn)定性滿足設(shè)計要求。

(2)通過有限元分析,采用總應(yīng)力法分別對壩基、壩體進行應(yīng)變及應(yīng)力計算。壩體的應(yīng)力分布從壩基開始逐漸向壩坡減小,切等值線與壩坡基本呈平行狀態(tài)。壩體在竣工期和蓄水期條件下,粘土心墻與混凝土防滲墻結(jié)構(gòu)的應(yīng)力基本保持一致,垂直應(yīng)力與水平方向應(yīng)力基本一致,但是位移出現(xiàn)較大變化,通過位移變化可以反饋確定接頭部位的破壞情況。

(3)通過理論計算和有限元計算對土壩壩體的穩(wěn)定性進行定量計算,兩計算結(jié)果總體趨勢基本保持一致,且計算結(jié)果之間的差異較小,可以作為進行工程設(shè)計和施工的科學(xué)依據(jù),并且可以作為工程修正設(shè)計和制定工程措施的重要依據(jù),保障工程設(shè)計、施工和運行管理的順利進行。

[1] 鐘祖良, 劉新榮, 吳強. 流固耦合作用下馬埡口隧道塌方機理及其災(zāi)害治理研究[J]. 工程勘察, 2011(07): 18- 22.

[2] Duncan J M, et al. Strength stress- strain and bu1k modulus parameters of Finite Element Analysis of stresses and movements in soil masses[R]. Report No. VCB/GT/78- 02. Univ. of California Berkeley, 1978: 47- 61.

[3] 顧淦臣. 土石壩的現(xiàn)狀水平和我國土石壩的前景[M]. 南京: 河海大學(xué)出版社, 1984: 1- 5.

[4] 包承綱, 錢勝國, 馬時冬. 三峽二期圍堰下淤砂層的動力特性及有關(guān)工程問題的研究[J]. 巖土工程學(xué)報, 2000 (04): 402- 407.

[5] 徐江, 趙其興, 袁曉華, 等. ANSYS在混凝土拱壩溫度場仿真計算中的應(yīng)用[J]. 水利規(guī)劃與設(shè)計, 2008(04): 75- 78+98.

[6] 劉先珊, 劉少煒, 余永康. 土石圍堰滲流場分析及穩(wěn)定性評價[J]. 遼寧工程技術(shù)大學(xué)學(xué)報(自然科學(xué)版), 2009(S1): 212- 215.

[7] 夏建軍. 膠結(jié)材料筑壩技術(shù)在白土嶺水電站的應(yīng)用[J]. 中國水能及電氣化, 2014(11): 14- 16.

[8] 劉奉銀, 劉國棟, 戚長軍. 某水電站防滲墻深度對圍堰滲流特性及地震穩(wěn)定性影響的計算研究[J]. 巖土力學(xué), 2009 (S2): 297- 301+311.

[9] 袁劍軍. 土石壩防滲墻黏土混凝土材料的工程應(yīng)用初探[J]. 水利建設(shè)與管理, 2013(02): 33- 35.

[10] 彭華, 陳勝宏. 飽和-非飽和巖土非穩(wěn)定滲流有限元分析研究[J]. 水動力學(xué)研究與進展, 2002(02): 253- 259.

[11] Duncan J M, Zhang CY. 1970. Nonlinear analysis of stress and strain in soils[J]. Journal of Soil Mechanics and Foundation Division, 96(SM5): 1629- 1653.

[12] Zhang J Y, Li Y, Xuan G X, et al . Overtopping breaching of cohesive homogeneous earth dam with different cohesive strength[J]. Science in China, Series E, 2009, 52(10): 3024- 3029.

[13] 李淑貞, 周艷麗, 王霞, 等. 水環(huán)境監(jiān)測實驗室質(zhì)量管理體系的構(gòu)建[J]. 水利技術(shù)監(jiān)督, 2015(06): 1- 2+14.