張 鐸,劉洪亮

(1.深圳粵港工程技術(shù)服務(wù)有限公司,廣東 深圳 518001;2.長(zhǎng)江勘測(cè)規(guī)劃設(shè)計(jì)研究有限責(zé)任公司,武漢 430014)

0 引 言

混凝土面板堆石壩因其具有施工便捷、結(jié)構(gòu)簡(jiǎn)單、壩體穩(wěn)定、工程造價(jià)低以及適應(yīng)性良好等優(yōu)點(diǎn),被廣泛應(yīng)用于水利工程建設(shè)中[1]。隨著水利工程建設(shè)的不斷發(fā)展,混凝土面板堆石壩的相關(guān)施工技術(shù)也隨之更新,但混凝土面板堆石壩仍然存在著滲流破壞導(dǎo)致堆石壩面板開裂的問(wèn)題[2]。

有限元數(shù)值模擬軟件COMSOL Multiphysics是一種功能強(qiáng)大的建模軟件,可以利用該軟件對(duì)堆石壩壩體進(jìn)行模型建模,并通過(guò)數(shù)值模擬,對(duì)堆石壩面板以及防滲墻不同工種狀態(tài)下的滲流場(chǎng)進(jìn)行計(jì)算分析[3]?？紤]到實(shí)際工程中混凝土面板堆石壩處于各種水壓環(huán)境,為了研究水壓力對(duì)裂縫滲透性的影響,通過(guò)數(shù)值計(jì)算和統(tǒng)計(jì)分析,對(duì)水壓下裂縫滲透系數(shù)的計(jì)算方法進(jìn)行分析。同時(shí),為了研究堆石壩面板防滲體的不同特性,對(duì)水壓力下混凝土面板堆石壩進(jìn)行滲流數(shù)值模擬。根據(jù)面板所處位置不同的水壓、不同防滲墻開裂高度以及不同破壞范圍的大小等因素,對(duì)滲流場(chǎng)的滲流量、浸潤(rùn)線以及水頭等數(shù)值進(jìn)行比較分析[4]。通過(guò)對(duì)水壓力下混凝土面板堆石壩滲流的分析研究,以期為水利工程的滲流問(wèn)題提供參考。

1 水壓力下混凝土面板堆石壩滲流數(shù)值模擬分析

1.1 水壓力下混凝土面板堆石壩裂縫滲透系數(shù)計(jì)算方法分析

混凝土面板堆石壩施工技術(shù)從早期的木面板拋填到如今的混凝土面板碾壓,不管是施工技術(shù)還是施工理論都在不斷發(fā)展與完善。當(dāng)前的混凝土面板堆石壩因其具有施工便捷、結(jié)構(gòu)簡(jiǎn)單、壩體安全穩(wěn)定、機(jī)械化程度高、工程造價(jià)低以及適應(yīng)性良好等優(yōu)點(diǎn),已成為水利工程領(lǐng)域中的優(yōu)選壩型[5]。

因面板開裂導(dǎo)致的滲流是混凝土面板堆石壩長(zhǎng)期存在的工程問(wèn)題。根據(jù)以往的混凝土面板堆石壩裂縫統(tǒng)計(jì)分析,可將裂縫按照寬度分為3類:Ⅰ類是指小于0.3mm的裂縫;Ⅱ類是指0.3～1.0mm的裂縫;Ⅲ類是指大于等于1.0mm的裂縫。在實(shí)際工程中,混凝土面板堆石壩最為常見(jiàn)的是Ⅱ類裂縫[6]。

根據(jù)工程實(shí)例研究發(fā)現(xiàn),混凝土面板堆石壩產(chǎn)生裂縫的位置主要集中于面板的左右兩岸和面板的中下部。其中,導(dǎo)致面板左右兩岸產(chǎn)生裂縫的主要原因是壩體兩岸產(chǎn)生了沉降變形,而導(dǎo)致面板中下部產(chǎn)生裂縫的主要原因是壩體發(fā)生沉降、堆石體產(chǎn)生流變以及面板產(chǎn)生變形等[7]?；炷撩姘宥咽瘔萎a(chǎn)生開裂的原因有不均勻沉降變形、施工不規(guī)范、干縮應(yīng)力、溫度應(yīng)力以及地震等不可抗力因素等,其中地震以及施工因素占比較小。因此可以看出,不均勻沉降變形、干縮應(yīng)力和溫度應(yīng)力等因素是導(dǎo)致混凝土面板堆石壩產(chǎn)生開裂的主要原因。

目前,針對(duì)滲流計(jì)算分析的方法有多種,大體可分為試驗(yàn)?zāi)M法和理論分析法兩大類[8]。具體見(jiàn)圖1。

圖1 滲流計(jì)算分析方法

隨著科技的發(fā)展,試驗(yàn)?zāi)M法具有試驗(yàn)和解題等雙重功能,可以直觀觀察到滲流現(xiàn)象。因此,基于試驗(yàn)?zāi)M法的計(jì)算機(jī)數(shù)值模擬技術(shù)逐漸成為當(dāng)前的主流技術(shù)。

在實(shí)際工程中產(chǎn)生的裂縫往往位于水下,而根據(jù)裂縫所處的位置不同,裂縫內(nèi)部滲流受到所產(chǎn)生的水壓力影響也各不相同。不同水壓會(huì)產(chǎn)生流態(tài)、流速以及滲流量等的變化,從而改變整個(gè)混凝土面板堆石壩的滲流場(chǎng)分布。

Comsol Multiphysics是一種常被應(yīng)用于裂縫巖體滲流計(jì)算的多物理場(chǎng)耦合模擬仿真軟件,利用該軟件對(duì)水壓力下混凝土面板堆石壩進(jìn)行滲流分析[9]。設(shè)定水壓力為P,水壓力范圍取值設(shè)定為P∈[0.2,0.4,…,1.0Mpa],裂縫寬度為w,裂縫寬度的取值范圍設(shè)定為w∈[0.3,0.5,…,1.1mm],設(shè)定裂縫深度d=200mm,根據(jù)軟件模擬,可以得到裂縫流量以及滲透系數(shù)的數(shù)學(xué)表達(dá)式:

(1)

式中:F為通過(guò)裂縫的流量;K為滲透系數(shù);α、β均為擬合參數(shù);A為過(guò)流面積;G為水力梯度。

將裂縫寬度與擬合參數(shù)在Comsol Multiphysics軟件中進(jìn)行擬合,可得到裂縫寬度與擬合參數(shù)的關(guān)系,公式如下:

(2)

由式(2)可以看出,擬合參數(shù)α與裂縫寬度e呈非線性關(guān)系。將式(2)代入式(1),可以得到在水壓力下滲透系數(shù)與裂縫寬度的關(guān)系,公式如下:

(3)

綜上所述,混凝土面板堆石壩的壩體上產(chǎn)生裂縫會(huì)引起滲流問(wèn)題,從而影響堆石壩的安全性和穩(wěn)定性。為了提高混凝土面板堆石壩的使用性能,對(duì)壩體進(jìn)行滲流控制至關(guān)重要。面板開裂滲流控制的要求、方法及措施見(jiàn)圖2。

圖2 滲流控制的要求、方法及措施

混凝土面板作為面板堆石壩的關(guān)鍵性防滲結(jié)構(gòu),一旦產(chǎn)生裂縫,需要及時(shí)采取相關(guān)措施來(lái)進(jìn)行滲流控制[10]。針對(duì)已出現(xiàn)裂縫的面板,可通過(guò)充填法、注入法以及表面覆蓋法等方式進(jìn)行滲流控制[11]。此外,還可以從增強(qiáng)面板本身的抗裂性能方面來(lái)提高防滲效果。增強(qiáng)面板自身的抗裂性能可以采取的措施有:對(duì)面板堆石壩的堆石進(jìn)行合理的優(yōu)化和分區(qū);采用防滲性以及變形特性均較好的瀝青混凝土面板;在面板自身的防滲基礎(chǔ)上鋪筑土工薄膜來(lái)加固防滲結(jié)構(gòu)等。

1.2 水壓力下混凝土面板堆石壩滲流數(shù)值模擬

為了對(duì)水壓力下混凝土面板堆石壩進(jìn)行滲流數(shù)值模擬,首先利用軟件構(gòu)建堆石壩的三維算例模型。設(shè)置堆石壩的總長(zhǎng)度540m,水流方向設(shè)置為沿Y軸正方向,上游和下游的河床長(zhǎng)度分別設(shè)置為106和118m。設(shè)置堆石壩模型的總寬度545m,其中壩體寬度345m,兩岸巖體寬度各為100m。設(shè)置堆石壩位于Z坐標(biāo)平面,Z軸高度107m則為大壩高度。堆石壩的三維算例模型見(jiàn)圖3。

圖3 堆石壩的三維算例模型圖

由圖3可知,堆石壩的三維算例模型中含有主次堆石體、防滲墻、混凝土面板、覆蓋層、墊層以及兩岸基巖等區(qū)域[12]。邊界條件設(shè)定基巖的底部、上下游邊界以及兩岸巖體的外側(cè)邊界均為不透水邊界。堆石壩壩體以及壩體兩岸巖體在下游水位以上與大氣接觸的逸出邊界,稱之為出水邊界[13]。該模型已知在正常蓄水位的情況下,上游施加200m水頭,下游施加120m水頭。由于只考慮非降雨條件下的穩(wěn)定水位,且不考慮地下水,因此上下游水庫(kù)以下的水位邊界就是已知的水位邊界。通過(guò)不斷校正下游潛在的逸出邊界,就可以確定下游壩體產(chǎn)生的滲流面,從而計(jì)算出大壩的自由面[14]。

根據(jù)一般工程資料,確定對(duì)應(yīng)的材料滲透系數(shù),見(jiàn)表1。

表1 材料滲透系數(shù)

在實(shí)際工程中,混凝土面板堆石壩的面板上,往往會(huì)存在大量的、分布不均勻的裂縫。以河床為基礎(chǔ)面,裂縫的位置高度大部分在0～101m之間,寬度一般在0.1～1.4mm之間,裂縫的粗糙度一般在0.1～0.45之間[15]。將表1中的滲透系數(shù)代入裂縫參數(shù),即可計(jì)算出各裂縫的滲透系數(shù)。根據(jù)面板的上游水位以及面板的長(zhǎng)度,可以選取4種水壓強(qiáng)度范圍。根據(jù)水位從高到低,水壓范圍可分為:范圍一為0～0.15MPa;范圍二為0.15～0.48MPa;范圍三為0.48～0.80MPa;范圍四為0.80～1.01MPa。將對(duì)應(yīng)位置高度的裂縫在以上范圍一至范圍四的4種水壓力強(qiáng)度下,進(jìn)行滲透系數(shù)的計(jì)算,得到裂縫滲透系數(shù)分別為2.48×10-6、7.58×10-6、4.42×10-5以及5.29×10-5m/s。

利用軟件對(duì)以上數(shù)據(jù)進(jìn)行數(shù)值模擬,可以得出堆石壩面板在不同水壓力下的壩內(nèi)流速矢量。當(dāng)范圍一發(fā)生破裂時(shí),流速范圍為4.0×10-5～5.5×10-5m/s,其中最小流速位于壩體上游坡腳內(nèi)部;當(dāng)范圍二發(fā)生破裂時(shí),流速范圍為1.5×10-5～2.3×10-4m/s;當(dāng)范圍三發(fā)生破裂時(shí),流速范圍為6.0×10-5～8.1×10-4m/s;當(dāng)范圍四發(fā)生破裂時(shí),上游壩腳處的流速最大為4.5×10-4m/s。

綜合可以看出,混凝土面板堆石壩產(chǎn)生面板破裂時(shí),滲流場(chǎng)流速中,高流速部位的分布主要集中在面板破裂處和下游壩腳處;低流速部位主要分布在上游壩腳內(nèi)部。堆石壩壩體上的裂縫位置由頂部向底部移動(dòng)時(shí),流速增大,流量也隨之增大;而當(dāng)面板上的裂縫位于底部時(shí),滲流量變化不明顯。

2 水壓力下混凝土面板堆石壩滲流數(shù)值模擬驗(yàn)證

為了對(duì)水壓力下混凝土面板堆石壩滲流特性進(jìn)行驗(yàn)證,使用有限元模型進(jìn)行數(shù)據(jù)模擬。設(shè)置的模型總長(zhǎng)度為840m,模型整體沿Y軸方向,壩高從上游坡腳開始延伸4倍、從下游坡腳開始延伸1倍。整體模型只設(shè)置堆石壩主體結(jié)構(gòu),未設(shè)置左右兩岸的巖體,沿X軸方向的壩體寬度346.85m。設(shè)置堆石壩位于Z坐標(biāo)平面,Z軸高度為107m,則大壩主體高度為107m,其中覆蓋層43m,基巖57m。

為了計(jì)算方便,設(shè)置的邊界條件不考慮地下水和降雨情況,設(shè)置壩體左右兩岸的側(cè)面邊界以及壩體基巖的底部邊界均為不透水邊界,設(shè)定在正常水位下的上游和下游水位分別為200和120m,初始水頭205m。為了分析不同水壓力下產(chǎn)生開裂情況的混凝土面板堆石壩的滲流控制,設(shè)置3種滲流控制方案進(jìn)行對(duì)比分析,即無(wú)防滲墻的滲流控制方案設(shè)置為方案一;防滲墻深18m的滲流控制方案設(shè)置為方案二;防滲墻深30m的滲流控制方案設(shè)置為方案三,水平鋪蓋長(zhǎng)度設(shè)置范圍為[50,100,...,350]。將以上3種方案的滲流場(chǎng)計(jì)算結(jié)果進(jìn)行對(duì)比分析,得到不同方案與壩前水頭值的關(guān)系,見(jiàn)圖4。

圖4 不同方案與壩前水頭值的關(guān)系

由圖4可知,方案一中,當(dāng)水平鋪蓋為50～200m范圍內(nèi)時(shí),隨著水平鋪蓋長(zhǎng)度的增加,壩前水頭值對(duì)應(yīng)減少,減少幅度9.7%;當(dāng)水平鋪蓋大于200m時(shí),壩前水頭值維持在140m不變。方案二中,當(dāng)水平鋪蓋為50～250m范圍內(nèi)時(shí),隨著水平鋪蓋長(zhǎng)度的增加,壩前水頭值對(duì)應(yīng)減少,減少幅度16.1%;當(dāng)水平鋪蓋大于250m時(shí),壩前水頭值維持在130m不變。方案三中,隨著水平鋪蓋長(zhǎng)度由50m增加至300m,對(duì)應(yīng)的壩前水頭值從165m降低至130m,降低幅度21.2%;當(dāng)水平鋪蓋大于300m時(shí),壩前水頭值維持在130m不變。綜上可以看出,方案三的壩前水頭值下降幅度最大,與其他兩種方案相比具有更好的防滲效果。

為了驗(yàn)證防滲控制中防滲墻在實(shí)際應(yīng)用中的防滲效果,將3種方案下各部位的滲流量進(jìn)行對(duì)比分析,見(jiàn)圖5。

圖5 3種方案下各部位的滲流量對(duì)比圖

由圖5(a)可知,防滲墻深度分別為0、18和30m的3種方案下通過(guò)防滲墻的滲流量,都是隨著水平鋪蓋長(zhǎng)度的增加而減小。其中,當(dāng)水平鋪蓋為0m時(shí),防滲墻深度為0m的方案一達(dá)到最大滲流量為8.0×10-4m3/s,此時(shí)方案二、方案三的滲流量分別為7.0×10-4和6.0×10-4m3/s。當(dāng)水平鋪蓋為350m時(shí),方案一、方案二、方案三的滲流量分別為4.0×10-4、2.5×10-4和2.5×10-4m3/s。可以看出,通過(guò)防滲墻產(chǎn)生的最大流量差為4.5×10-4m3/s;當(dāng)水平鋪蓋為150m時(shí),3種方案均達(dá)到穩(wěn)定值。

由圖5(b)可知,當(dāng)水平鋪蓋為0m時(shí),方案一、方案二、方案三通過(guò)堆石區(qū)的滲流量分別為1.5×10-4、7.5×10-5和5.0×10-5m3/s。當(dāng)水平鋪蓋為350m時(shí),方案一、方案二、方案三的滲流量分別為5.0×10-5、4.5×10-5和4.5×10-5m3/s。通過(guò)堆石區(qū)產(chǎn)生的最大流量差為1.0×10-4m3/s;當(dāng)水平鋪蓋為150m時(shí),3種方案均趨于穩(wěn)定不變。

由圖5(c)可知,當(dāng)水平鋪蓋為0m時(shí),方案一、方案二、方案三通過(guò)壩基的滲流量分別為4.0×10-5、2.5×10-5和2.2×10-5m3/s。當(dāng)水平鋪蓋為350m時(shí),方案一、方案二、方案三的滲流量分別為2.3×10-5、1.5×10-5和1.5×10-5m3/s。通過(guò)壩基產(chǎn)生的最大流量差為1.7×10-5m3/s;當(dāng)水平鋪蓋為200m時(shí),3種方案均達(dá)到穩(wěn)定值。

綜合可以看出,水平鋪蓋的長(zhǎng)度對(duì)減小滲流量的作用有限,有效的最大長(zhǎng)度為200m。滲流控制中,有無(wú)設(shè)置防滲墻對(duì)滲流場(chǎng)的影響極為明顯,方案一產(chǎn)生的流量差最大;方案三的滲流控制效果最為顯著。

3 結(jié) 論

為了分析水壓力下混凝土面板堆石壩的滲流問(wèn)題,本文建立了一種水壓力下混凝土面板堆石壩裂縫滲透系數(shù)的計(jì)算方法,同時(shí)針對(duì)堆石壩面板防滲體的不同特性,對(duì)水壓力下混凝土面板堆石壩進(jìn)行了滲流數(shù)值模擬。結(jié)果顯示,方案一中,當(dāng)水平鋪蓋為50～200m范圍內(nèi)時(shí),壩前水頭值減少9.7%;方案二中,當(dāng)水平鋪蓋為50～250m范圍內(nèi)時(shí),壩前水頭值減少16.1%;方案三中,當(dāng)水平鋪蓋為50～300m范圍內(nèi)時(shí),壩前水頭值減少21.2%。通過(guò)防滲墻、堆石區(qū)以及壩基產(chǎn)生的最大流量差分別為4.5×10-4、1.0×10-4及1.7×10-5m3/s。研究表明,防滲墻深30m的滲流控制具有更好的防滲效果。