王冠昌，金居平

（浙江禹川勘測(cè)設(shè)計(jì)有限公司,浙江 桐鄉(xiāng) 314500）

0 引言

我國(guó)大部分地區(qū)水電資源豐富,隨著國(guó)家科技、經(jīng)濟(jì)的發(fā)展,大量的水利工程應(yīng)運(yùn)而生。而這些項(xiàng)目往往難度較大,對(duì)設(shè)計(jì)人員和施工人員的技術(shù)水平也是一種巨大的考驗(yàn)。對(duì)于規(guī)模較大的水利項(xiàng)目往往建立在強(qiáng)度高、承載能力強(qiáng)、抗變形能力強(qiáng)、結(jié)構(gòu)完整的基巖上,但實(shí)際上巖體通常都有著裂隙,水分會(huì)在這些裂隙中流通會(huì)對(duì)裂隙形成附加應(yīng)力,從而降低裂隙的有效應(yīng)力,對(duì)基巖的穩(wěn)定性造成負(fù)面影響[1-3]。大型的水利工程通常建設(shè)在變形小,結(jié)構(gòu)完整的基巖上。所以研究巖體滲流對(duì)水利工程有著非常關(guān)鍵的作用,眾多學(xué)者對(duì)此也進(jìn)行了廣泛的研究[4-7]。

基于此,為了研究巖體裂隙滲流對(duì)水利工程的影響,本研究以某大壩泄水閘為例,根據(jù)其基巖的結(jié)構(gòu)特點(diǎn),通過(guò)建立數(shù)值模型分析水閘裂隙巖體的滲流情況。

1 大壩泄水閘概況

實(shí)體重力壩為某水利樞紐大壩類(lèi)型,其壩高92 m。所研究泄水閘位于左岸階地,區(qū)域高程為24 m～44 m。郁江階第D1y1-2 層至第D1y1-4 層和那高嶺組的第D1n13-2 至第D1n13-4 層為泄水閘壩基的主要地層。此區(qū)域郁江階巖性主要是灰?guī)r、泥質(zhì)灰?guī)r,那高嶺組巖性主要是泥巖、泥質(zhì)粉砂巖和含泥細(xì)砂巖。巖層產(chǎn)狀表現(xiàn)為從下游側(cè)走向N15°～30°E,傾角為22°,傾向SE。主壩區(qū)地質(zhì)構(gòu)造類(lèi)型較多,比較復(fù)雜,發(fā)育又裂隙、斷層、褶皺、層面等,但各類(lèi)型構(gòu)造的發(fā)育程度因地貌單元而異。單斜巖層為主壩區(qū)主要類(lèi)型,產(chǎn)狀比較穩(wěn)定。泄水閘閘壩部位有8 條明顯出露的斷層,且都和壩軸線(xiàn)相交。通過(guò)實(shí)地勘測(cè)水電站28#、27#、26#壩段的地質(zhì)情況,發(fā)現(xiàn)場(chǎng)地裂隙大部分產(chǎn)狀離散,且為陡傾角。

2 離散元滲流模擬

2.1 建立滲流模型

因?yàn)榱严稁r體屬于非連續(xù)介質(zhì),所示選擇離散元程序來(lái)對(duì)其進(jìn)行模擬。由于對(duì)水閘壩基巖體的滲流進(jìn)行三維分析比較困難,并且壩基巖體地質(zhì)類(lèi)型較多,結(jié)構(gòu)分布復(fù)雜,很難全面地掌握壩基巖體的滲流情況。所以此次研究選擇二維模型來(lái)對(duì)壩基巖體的滲流進(jìn)行模擬。通過(guò)地質(zhì)勘查報(bào)告能夠得出28#壩段閘室下部存在一發(fā)育的斷層,其會(huì)對(duì)此壩段的滲流狀態(tài)造成嚴(yán)重的影響,所以此次研究將此壩段作為對(duì)象進(jìn)行模擬,通過(guò)離散元軟件來(lái)對(duì)泄水閘和其下部基巖的模型進(jìn)行建立,模型圖見(jiàn)圖1。

圖1 28#壩段模型圖

2.2 選擇模型參數(shù)

對(duì)于巖體參數(shù),主要包括土體抗剪強(qiáng)度、密度等,具體見(jiàn)表1。在進(jìn)行滲流過(guò)程中比較重要的參數(shù)有裂隙隙寬ares、裂隙初始隙寬ao和滲透系數(shù)kj,參數(shù)具體取值見(jiàn)表2,動(dòng)力學(xué)的相關(guān)參數(shù)為地表和基巖的地震動(dòng)峰值加速分別為0.171g和0.136g,特征周期均為0.45 s。

表1 巖體和帷幕相關(guān)參數(shù)

表2 結(jié)構(gòu)面水力學(xué)參數(shù)

2.3 選擇模型參數(shù)

選擇合適的邊界條件能夠使模型模擬得出的結(jié)果更加準(zhǔn)確,對(duì)于此次所用模型,因?yàn)樵O(shè)置的模型尺寸超過(guò)了實(shí)際重點(diǎn)研究區(qū),所以在此次所用模型中固定其左右邊界和底部邊界,將模型的邊界速度設(shè)置成0。再根據(jù)正常蓄水位（上下游水位分別為62 m 和23 m）來(lái)設(shè)置水壓力。離散元軟件能夠根據(jù)實(shí)際地震記錄值來(lái)模擬地震反應(yīng),因?yàn)椴荒苤苯釉O(shè)定邊界加速度,所以能夠?qū)⑵滢D(zhuǎn)變成應(yīng)力,并且考慮到在現(xiàn)實(shí)中整個(gè)項(xiàng)目都會(huì)有地震荷載,所以將邊界兩側(cè)的應(yīng)力值設(shè)置為最大值的一半。

3 分析滲流計(jì)算結(jié)果

3.1 有無(wú)防滲帷幕在無(wú)動(dòng)力作用下的結(jié)果分析

此次選擇28#壩段模型來(lái)計(jì)算其有帷幕和無(wú)帷幕時(shí)的滲流,并根據(jù)結(jié)果對(duì)其防滲帷幕的作用進(jìn)行分析和驗(yàn)證。模擬選擇的工況為大壩處于正常水位（上下游水位分別為62 m 和23 m）,水壓施加荷載時(shí)都垂直于上下游基巖。因?yàn)樵O(shè)置的模型大小要遠(yuǎn)遠(yuǎn)大于實(shí)際重點(diǎn)研究區(qū),所以在計(jì)算時(shí)設(shè)置邊界為不透水邊界。對(duì)壩段滲流模擬計(jì)算后得到裂隙巖體滲流結(jié)果,見(jiàn)圖2 和圖3,圖中線(xiàn)段代表模型中滲流的賦存情況。從圖2和3能夠得出,裂隙巖體滲流大多出現(xiàn)在導(dǎo)水管道中（裂隙、斷層等 。滲流在整體上形成一條閉合回路,與滲流在實(shí)際工程中的情況相接近。從圖3 能夠得到,當(dāng)壩段有帷幕時(shí),與圖2 無(wú)帷幕裂隙巖體相比其滲流量大大減小,特別是壩底部位的流量減小幅度較大,所以能夠得出防滲帷幕有較好的防水效果。

圖2 無(wú)帷幕時(shí)28#壩段滲流模擬圖

圖3 有帷幕時(shí)28#壩段滲流模擬圖

從上述模擬結(jié)果能夠得到,防滲帷幕能夠很好地隔絕裂隙巖體的滲流,為了進(jìn)一步分析防滲帷幕的防滲效果,對(duì)其開(kāi)展定量分析,在壩段布設(shè)了一些監(jiān)測(cè)點(diǎn),因?yàn)?8#壩段內(nèi)有一斷層F215 穿過(guò),所以在斷層左右和內(nèi)部設(shè)置監(jiān)測(cè)點(diǎn),具體見(jiàn)圖4。各監(jiān)測(cè)點(diǎn)中1、4、6、7、8 號(hào)點(diǎn)在同一個(gè)結(jié)構(gòu)面上,其主要功能是對(duì)壩底揚(yáng)壓力進(jìn)行記錄,并通過(guò)測(cè)得的揚(yáng)壓力來(lái)評(píng)價(jià)滲透壓力。同時(shí)有無(wú)帷幕時(shí)的工況相類(lèi)似,所以監(jiān)測(cè)點(diǎn)在有帷幕時(shí)的布置未發(fā)生改變。有無(wú)帷幕時(shí)28# 壩段流量比較圖見(jiàn)圖5（a）。從圖中能夠看出,模型中未設(shè)置防滲帷幕時(shí),監(jiān)測(cè)所得流量最大值、最低值和平均值分別是2.7×103m3/s、1.67×103m3/s 和2.27×103m3/s；模型中設(shè)置防滲帷幕后,監(jiān)測(cè)所得流量最大值、最低值和平均值分別是2.13×103m3/s、1.12×103m3/s 和1.73×103m3/s,結(jié)構(gòu)面內(nèi)的流量平均值縮小至0.54×103m3/s,降低幅度達(dá)23.3%。圖5（b）為28#壩段有無(wú)帷幕時(shí)各監(jiān)測(cè)點(diǎn)測(cè)得模型最終收斂時(shí)的揚(yáng)壓力比較結(jié)果。從圖5（b）中能夠得出,有帷幕時(shí)的揚(yáng)壓力要低于無(wú)帷幕工況,7 號(hào)測(cè)點(diǎn)位于帷幕左側(cè),所以無(wú)論帷幕是否存在,其數(shù)值都相同。同時(shí)有帷幕時(shí)和無(wú)帷幕時(shí)的平均揚(yáng)壓力值分別為0.43 MPa 和0.56 MPa,前者與后者相比降低了近30%。從圖5（b）能夠得出,8 號(hào)監(jiān)測(cè)點(diǎn)的揚(yáng)壓力明顯低于其余監(jiān)測(cè)點(diǎn),對(duì)此現(xiàn)象進(jìn)行分析發(fā)現(xiàn)8 號(hào)監(jiān)測(cè)點(diǎn)位置處于斷層內(nèi)部,繼而得出地下水流過(guò)斷層中的裂隙,且斷層中存在斷層泥,減小了斷層內(nèi)的裂隙寬度,從而導(dǎo)致滲流量和揚(yáng)壓力減小。

圖4 28#壩段監(jiān)測(cè)點(diǎn)布設(shè)情況

圖5 28#壩段各監(jiān)測(cè)點(diǎn)揚(yáng)壓力、滲流量對(duì)比圖

3.2 分析動(dòng)力作用下滲流變化情況

選擇28#壩段模型來(lái)對(duì)其動(dòng)力作用下滲流的變化進(jìn)行計(jì)算,模擬所選工況為大壩處于正常水位（上下游水位分別為62 m 和23 m）,圖6 和圖7 為裂隙巖體在有動(dòng)力荷載時(shí)的滲流情況。圖中線(xiàn)段代表模型中滲流的賦存情況。從圖6 和圖7 中能夠得到,裂隙巖體滲流大多出現(xiàn)在導(dǎo)水管道中（裂隙、斷層等）。滲流在整體上形成一條閉合回路,與滲流在實(shí)際工程中的情況相接近。但無(wú)論動(dòng)力作用是否存在,滲流場(chǎng)區(qū)域均未發(fā)生明顯變化,表示裂隙巖體滲流情況基本不會(huì)受到動(dòng)力作用的影響。

圖6 無(wú)動(dòng)力作用時(shí)28#壩段滲流模擬圖

圖7 有動(dòng)力作用時(shí)28#壩段滲流模擬圖

從圖6 和圖7 中能夠得出,加上動(dòng)力荷載后使?jié)B流量表現(xiàn)出略微增大,因此為了深入分析裂隙巖體是否會(huì)被動(dòng)力作用影響,選擇圖4 布設(shè)的監(jiān)測(cè)點(diǎn)來(lái)對(duì)其揚(yáng)壓力進(jìn)行分析。

有動(dòng)力作用時(shí)28#壩段流量比較圖見(jiàn)圖8（a）。從圖8（a）中能夠得出,加上動(dòng)力荷載能夠明顯增大滲流量,滲流量在無(wú)動(dòng)力工況和有動(dòng)力工況的平均值分別為2.13×10-3m3/s和2.66×10-3m3/s,提高了近19.9%。圖8（b）為28#壩段有無(wú)帷幕時(shí)各監(jiān)測(cè)點(diǎn)測(cè)得模型最終收斂時(shí)的揚(yáng)壓力比較結(jié)果。從圖8（b）中能夠得出,揚(yáng)壓力在有動(dòng)力荷載時(shí)較大,且揚(yáng)壓力在有帷幕工況和無(wú)帷幕工況下的平均值分別為0.743 MPa 和0.589 MPa,增加近20.7%。綜上所述,能夠得出當(dāng)結(jié)構(gòu)面中孔隙水壓力顯著提高時(shí),表示外部荷載會(huì)對(duì)巖體裂隙造成一定的影響,使其出現(xiàn)體積應(yīng)變,繼而增大了孔隙水壓力。

圖8 有無(wú)動(dòng)力時(shí)28#壩段各監(jiān)測(cè)點(diǎn)滲流量和揚(yáng)壓力對(duì)比圖

4 結(jié)論

本研究以某大壩泄水閘為例,根據(jù)其基巖的結(jié)構(gòu)特點(diǎn),通過(guò)建立數(shù)值模型分析了水閘裂隙巖體的滲流情況,得出以下結(jié)論：

1）未設(shè)置防滲帷幕時(shí),監(jiān)測(cè)所得流量最大值、最低值和平均值分別是2.7×103m3/s、1.67×103m3/s 和2.27×103m3/s；模型中設(shè)置防滲帷幕后,監(jiān)測(cè)所得流量最大值、最低值和平均值分別是2.13×103m2/s、1.12×103m2/s 和1.73×103m3/s,結(jié)構(gòu)面內(nèi)的流量平均值減小了0.54×103m3/s,減小了23.3%。

2）有帷幕時(shí)的揚(yáng)壓力要低于無(wú)帷幕工況,有帷幕時(shí)和無(wú)帷幕時(shí)的平均揚(yáng)壓力值分別為0.43 MPa 和0.56 MPa,前者與后者相比降低了近30%。

3）加上動(dòng)力荷載能夠明顯增大滲流量,滲流量在無(wú)動(dòng)力工況和有動(dòng)力工況的平均值分別為2.13×10-3m3/s 和2.66×10-3m3/s,提高了近19.9%；揚(yáng)壓力在有帷幕工況和無(wú)帷幕工況下的平均值分別為0.743 MPa 和0.589 MPa,增加近20.7%。