彭亞敏， 沈振中， 甘 磊

(河海大學(xué) 水利水電學(xué)院， 江蘇 南京 210098)

我國(guó)西部大型水電站裝機(jī)容量大、水頭高，且位于高山峽谷地區(qū)，地下水位高，因此，水工洞室常深埋于地下水位以下，需要考慮洞室襯砌的外水壓力作用和防滲排水措施，即如何選擇合適的防滲排水處理方案，高效且經(jīng)濟(jì)地減小襯砌所受滲透壓力，提高水工洞室的整體穩(wěn)定性[1]。關(guān)于水工洞室防滲排水的研究已有許多成果。班宏泰等[2]研究了在高壓隧洞中水平段和下水平段布置排水洞和排水孔幕的滲流場(chǎng)分布特征，結(jié)果表明該方案有利于降低外水壓力。李巖松等[3]對(duì)固結(jié)灌漿前后的引水隧洞襯砌所受的壓力進(jìn)行分析，認(rèn)為通過不同灌漿施工方法可有效減小襯砌所受壓力。朱伯芳等[4]通過三維有限元計(jì)算，分析結(jié)果表明孔徑對(duì)排水效果影響很小，而孔距、孔深在一定范圍內(nèi)影響較大。簡(jiǎn)文彬等[5]利用波速計(jì)算研究了固結(jié)灌漿在不同方向上的作用效果。根據(jù)《水工隧洞設(shè)計(jì)規(guī)范》(SL 279—2016)，當(dāng)需要降低地下水變化對(duì)環(huán)境影響和水土保持的要求時(shí)，水工洞室需進(jìn)行固結(jié)灌漿。固結(jié)灌漿作為水工隧洞提高圍巖穩(wěn)定性的有效措施之一，有利于加固巖體、封閉洞室周圍巖體裂隙，大幅度減小地下水外滲，從而增強(qiáng)洞室整體的抗變形能力并減小地下水位變化[6]。當(dāng)水工洞室位于地下水位以下，無地下水環(huán)境影響要求時(shí)，可以設(shè)置排水措施，降低外水壓力強(qiáng)度，如許多地下廠房采用防滲帷幕、排水孔幕以及廠房頂?shù)娜俗中闻潘椎却胧?。但是，?dāng)水工洞室位于地下水位以下，需要考慮地下水環(huán)境影響時(shí)，既要減小襯砌所受外水壓力，又要盡可能減小排水對(duì)地下水位的影響，這時(shí)需要考慮其他防滲排水措施[7]。

本文結(jié)合某水電站工程泄洪洞設(shè)計(jì)，考慮排水孔和固結(jié)灌漿組合方案，建立了泄洪洞圍巖三維有限元模型，精細(xì)模擬了泄洪洞襯砌和排水孔，圍巖材料分區(qū)及固結(jié)灌漿，計(jì)算分析地下泄洪洞滲流場(chǎng)分布和襯砌所受外水壓力，研究了洞室設(shè)置淺排水孔、圍巖進(jìn)行深固結(jié)灌漿的作用效果，提出防滲排水處理建議。

1 某水電站泄洪洞布置

某水電站工程攔河壩為礫石土心墻堆石壩，泄水建筑物布置于左岸，其中洞式溢洪道軸線穿越大壩帷幕，最大壩高為293 m。壩區(qū)河道下切十分強(qiáng)烈，河谷狹窄、谷坡陡峻，為典型的高山峽谷。該區(qū)域出露地層以三疊系上統(tǒng)為主，屬淺變質(zhì)巖系，具有巖性較為單一、巖層延展變化大、分布較廣的特點(diǎn)。深孔泄洪洞、放空洞、豎井泄洪洞位于大壩現(xiàn)有防滲帷幕線以外。深孔泄洪洞全長(zhǎng)1 640 m,最大垂直埋深350 m,洞段周圍巖體質(zhì)量為Ⅲ2～Ⅴ類，具有足夠的承載及抗變形能力，可滿足洞段應(yīng)力及變形要求。深孔泄洪洞典型斷面如圖1所示。

深孔泄洪洞頂拱部位布置排水孔，排水孔間距為3.0 m。排水孔直徑為50 mm，孔深入巖1 m，2個(gè)排水孔軸線夾角為20°，相鄰兩排奇偶數(shù)交錯(cuò)布置。沿縱軸線深孔泄洪洞頂部奇偶孔布置見圖2。

圖1 深孔泄洪洞典型斷面(單位：cm)Fig.1 Typical section of deep hole spillway tunnel(unit: cm)

圖2 泄洪洞排水孔布置平面投影Fig.2 Plane projection of drainage holes in spillway tunnel

圖3 泄洪洞有限元模型分區(qū)示意Fig.3 Partitions of finite element model of drainage hole

2 泄洪洞圍巖有限元模型

為方便建立三維有限元模型，按照面積等效原則將圓形排水孔等效為44 mm×44 mm的正方形排水孔，不影響模型整體滲流計(jì)算且便于三維有限元建模。固結(jié)灌漿深度為8 m，根據(jù)固結(jié)灌漿孔布置及其灌漿效果，將圍巖灌漿區(qū)域分區(qū)并按實(shí)際效果調(diào)整材料滲透系數(shù)[8]。泄洪洞有限元模型分區(qū)如圖3所示。

圖4 泄洪洞襯砌及排水孔三維有限元網(wǎng)格Fig.4 3D FEM mesh of spillway tunnel lining and drainage holes

2.1 網(wǎng)格及邊界條件

選取泄洪洞地下水位最高的洞段，以泄洪洞中央為基準(zhǔn)，建立泄洪洞及其圍巖三維有限元模型。該模型模擬范圍沿垂直泄洪洞軸線方向(X方向)長(zhǎng)190 m，順河流方向(Y方向)寬15.176 m，即泄洪洞0+685.00 m～ 0+700.176 m，高程方向(Z方向)約196 m。計(jì)算的邊界模型如下：①模型底部邊界取至固結(jié)灌漿最深區(qū)域往外延伸80 m，約為4倍隧洞高度，高程2 642 m。②模型頂部邊界取至固結(jié)灌漿最深區(qū)域往外延伸80 m,約為4倍隧洞高度，高程為2 838 m。③模型左右邊界沿著固結(jié)灌漿最深區(qū)域往外延伸80 m，約為4倍隧洞高度。根據(jù)泄洪洞、排水孔等的實(shí)際尺寸和布置，選取控制斷面10個(gè)。泄洪洞、排水孔等均按實(shí)際尺寸離散，采取控制斷面超單元自動(dòng)剖分技術(shù)，形成超單元結(jié)點(diǎn)2 680個(gè)，超單元1 991個(gè)[9]。加密細(xì)分后生成有限元網(wǎng)格格點(diǎn)總數(shù)37 380個(gè)，單元總數(shù)為33 768個(gè)。泄洪洞有限元模型襯砌及排水孔三維有限元網(wǎng)格如圖4所示。有限元模型的計(jì)算邊界條件如下：左、右側(cè)截取邊界為已知水頭邊界，沿洞軸線方向模型長(zhǎng)度不大，僅15.176 m，忽略其滲流運(yùn)動(dòng)，即上游以及下游側(cè)截取邊界近似為不透水邊界，模型底邊界近似為不透水邊界。邊界水頭由樞紐區(qū)滲流場(chǎng)計(jì)算分析確定，取最高地下水位2 784.00 m[10]。

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

由于頂拱排水孔呈輻射向布置，越往圍巖深部灌漿孔間距越大，固結(jié)效果也越差，因此考慮到固結(jié)灌漿的實(shí)際效果，將固結(jié)灌漿區(qū)劃分為4個(gè)區(qū)(如圖3)，并取不同的滲透系數(shù)[11]。泄洪洞圍巖各分區(qū)滲透系數(shù)為：混凝土1×10-7cm/s，抗沖磨混凝土1×10-7cm/s， 固結(jié)灌漿1區(qū)1×10-6cm/s，固結(jié)灌漿2區(qū)2×10-6cm/s，固結(jié)灌漿3區(qū)3×10-6cm/s，固結(jié)灌漿4區(qū)1×10-6cm/s，圍巖2.7×10-5cm/s。

3 防滲排水作用分析

3.1 防滲排水方案

根據(jù)泄洪洞防滲排水措施的可能情況，擬定4種防滲排水處理方案(見表1)。模型采用鄧肯-張(E-B)模型，進(jìn)行三維非線性有限元計(jì)算，分析水工隧洞在不同防滲排水處理方案下浸潤(rùn)面以及襯砌所受滲水壓力的變化。

表1 泄洪洞防滲排水處理方案Tab.1 Anti-seepage treatment scheme for spillway tunnel

3.2 地下水位變化分析

選取經(jīng)過偶數(shù)孔中心的橫剖面分析泄洪洞圍巖滲流場(chǎng)，方案3的滲流場(chǎng)位勢(shì)分布如圖5所示。當(dāng)泄洪洞圍巖進(jìn)行深固結(jié)灌漿后，泄洪洞圍巖灌漿區(qū)地下水等勢(shì)線密集，其水力梯度遠(yuǎn)大于外側(cè)圍巖和內(nèi)側(cè)排水孔區(qū)，深部圍巖的固結(jié)灌漿可減少圍巖地下水內(nèi)滲[12]。

不同防滲排水處理方案下圍巖地下水位線見圖6。由圖6可知，當(dāng)泄洪洞圍巖進(jìn)行深部固結(jié)灌漿后，深部圍巖的固結(jié)灌漿可有效減少圍巖地下水內(nèi)滲，地下水位線變化不明顯，泄洪洞拱頂上方地下水位距左右邊界水位下降2.35 m。當(dāng)頂拱設(shè)置排水孔之后，泄洪洞上方地下水位距左右邊界水頭下降27.17 m。因此對(duì)于深埋水工隧洞，設(shè)置排水措施會(huì)引起圍巖地下水位不可恢復(fù)的下降，破壞地下水環(huán)境，對(duì)生態(tài)環(huán)境造成一定的影響[13]。圍巖采取深固結(jié)灌漿且設(shè)置排水孔時(shí)，泄洪洞洞室上方的地下水頭下降了6.48 m。當(dāng)深部圍巖灌漿深度改為5 m，泄洪洞上方距左右邊界水頭下降4.92 m。固結(jié)灌漿深度減小，地下水位變化變小，對(duì)地下水環(huán)境影響也更小?？芍?dāng)需考慮地下水環(huán)境影響時(shí)，可采取固結(jié)灌漿防滲排水措施，可有效減小地下水位的降低，保護(hù)生態(tài)環(huán)境[14]。

圖5 方案3泄洪洞圍巖滲流場(chǎng)位勢(shì)分布(單位:m)Fig.5 Seepage equipotential line distribution of surrounding rock in scheme 3 (unit:m)

圖6 泄洪洞偶數(shù)孔中心橫剖面地下水位線變化Fig.6 Changes of groundwater level at center section of even numbered holes

3.3 襯砌上所受外水壓力分析

對(duì)于深埋地下的泄洪洞，混凝土襯砌材料本身具有一定的滲透性，襯砌有利于減小糙率和保護(hù)圍巖[15]。泄洪洞橫剖面頂拱襯砌上的滲壓水頭分布見圖7，x坐標(biāo)值為襯砌沿拱圈展開坐標(biāo)，x=0表示隧洞拱頂，縱坐標(biāo)表示該位置上襯砌所受滲水壓力值。當(dāng)排水孔失效時(shí)，襯砌上滲水壓力值沿著x方向呈先減小后增大的趨勢(shì)。排水孔正常作用時(shí)排水孔處滲水壓力值為0，相鄰排水孔之間滲水壓力值呈拋物線趨勢(shì)。方案1滲壓水頭最大值為23.35 m，位于模型左、右邊界兩側(cè)。圍巖進(jìn)行深固結(jié)灌漿可減小對(duì)地下水環(huán)境的影響，但可對(duì)襯砌造成額外的附加力。方案2滲壓水頭在左、右兩側(cè)相鄰排水孔之間達(dá)到最大值，最大值為3.43 m，降低了85.31%。當(dāng)襯砌設(shè)置排水孔可有效發(fā)揮排水降壓的功效，大幅減小作用于襯砌上的滲透壓力。方案3滲壓水頭也在左、右兩側(cè)相鄰排水孔之間達(dá)到最大值，最大值為9.83 m，降低了57.9%。采用排水孔和固結(jié)灌漿聯(lián)合防滲加固措施可減小其洞室襯砌所受滲透壓力，同時(shí)提高巖體整體抗變形能力。方案4滲壓水頭最大值為13.74 m，相較于方案3增加了39.78%。固結(jié)灌漿深度減小會(huì)造成襯砌滲壓水頭的增加。

取隧洞橫剖面為研究面，沿襯砌拱頂向圍巖深度方向的滲壓水頭分布見8，x坐標(biāo)值表示襯砌沿圍巖深度方向深度大小，縱坐標(biāo)表示該點(diǎn)所受滲水壓力值。方案1襯砌厚度方向的最大滲壓水頭為17.35 m；方案2為1.68 m，降低了90.32%；方案3的最大滲壓水頭為3.91 m，降低了78.33%。方案4的最大滲壓水頭為6.73 m，比方案3增加了72.12%。分析表明,隧洞圍巖設(shè)置排水孔可以有效減小襯砌所受的滲透壓力。

圖7 偶數(shù)孔中心橫剖面頂拱襯砌上的滲壓水頭分布(沿拱圈展開)Fig.7 Seepage pressure distribution of top arch tunnel lining at section of even numbered holes

圖8 襯砌洞軸線剖面沿圍巖深度方向的滲壓水頭分布Fig.8 Seepage pressure distribution of lining structure along depth of surrounding rock

4 結(jié) 語

本文針對(duì)深埋水工隧洞保護(hù)地下水環(huán)境要求和減小襯砌承受的滲透壓力的情況，結(jié)合實(shí)際工程，建立了泄洪洞三維有限元模型，詳細(xì)模擬了襯砌、排水孔和固結(jié)灌漿區(qū)，分析研究了不同處理方案下泄洪洞圍巖地下水滲流場(chǎng)的變化和襯砌上所受滲透壓力的分布規(guī)律，提出了相應(yīng)的防滲排水處理方案。

(1)當(dāng)需要考慮地下水環(huán)境影響時(shí)，可對(duì)圍巖進(jìn)行深部固結(jié)灌漿，范圍為5～8 m。分析表明固結(jié)灌漿可將深埋水工隧洞對(duì)周圍環(huán)境地下水位的影響降至最低，封閉泄洪洞圍巖裂隙，減少圍巖地下水內(nèi)滲，改善隧洞襯砌的受力狀態(tài)。

(2)當(dāng)不需考慮地下水環(huán)境影響，水工隧洞低于地下水位線且需要降低襯砌外水壓力時(shí)，可在襯砌頂拱位置設(shè)置淺排水孔，排水孔深度為1～2 m。淺排水可以大幅度降低襯砌所受的外水壓力，提高其滲透穩(wěn)定性。

(3)當(dāng)?shù)叵滤h(huán)境影響和襯砌外水壓力兩者均需要考慮時(shí)，可采用深部固結(jié)灌漿聯(lián)合淺排水的防滲排水措施，可減小襯砌所受的外水壓力，并降低對(duì)地下水環(huán)境的影響。

[1] 信春雷. 不同防排水模式對(duì)山嶺隧道襯砌水壓力影響關(guān)系研究[D]. 成都： 西南交通大學(xué)， 2011. (XIN Chunlei. Study on the impact of different waterproof and drainage patterns upon water pressure load on mountain tunnel lining[D]. Chengdu: Southwest Jiaotong University, 2011. (in Chinese))

[2] 班宏泰， 劉昌軍， 丁留謙， 等. 深埋地下管道及地下洞室群三維滲流場(chǎng)的有限元精細(xì)模擬[J]. 水利水電技術(shù), 2012, 43(9): 33- 38. (BAN Hongtai， LIU Changjun， DING Liuqian， et al. Precise finite element simulation on 3-D seepage field of deep-embedded underground pipelines and underground cavern group[J]. Water Resources and Hydropower Engineering, 2012, 43(9): 33- 38. (in Chinese))

[3] 李巖松， 陳壽根， 周澤林. 固結(jié)灌漿前、后引水隧洞襯砌力學(xué)行為研究[J]. 水利與建筑工程學(xué)報(bào), 2014, 12(2): 193- 198. (LI Yansong, CHEN Shougen, ZHOU Zelin. Mechanical behaviors of water diversion tunnel lining before and after consolidation grouting[J]. Journal of Water Resources and Architectural Engineering, 2014, 12(2): 193- 198. (in Chinese))

[4] 朱伯芳， 李玥， 張國(guó)新. 滲流場(chǎng)中排水孔直徑、間距及深度對(duì)排水效果的影響[J]. 水利水電技術(shù), 2008， 39(3)： 27- 29. (ZHU Bofang, LI Yue, ZHANG Guoxin. The influences of the diameter, spacing and depth of drain holes on the effect of draining in the seepage field[J]. Water Resources and Hydropower Engineering, 2008， 39(3)： 27- 29. (in Chinese))

[5] 簡(jiǎn)文彬， 張登， 許旭堂. 基于波速測(cè)試的裂隙巖體固結(jié)灌漿效果分析[J]. 巖土力學(xué), 2014, 35(7): 1945- 1949. (JIAN Wenbin, ZHANG Deng, XU Xutang. Analysis of consolidation grouting effect of fractured rock mass based on wave velocity test[J]. Rock and Soil Mechanics, 35(7): 1945- 1949. (in Chinese))

[6] 張巍， 陳云長(zhǎng)， 黃立財(cái)， 等. 高壓隧洞內(nèi)水外滲三維有限元分析與滲透穩(wěn)定性研究[J]. 水利與建筑工程學(xué)報(bào), 2015, 13(5): 212- 217. (ZHANG Wei, CHEN Yunchang, HUANG Licai, et al. 3D finite element analysis of seepage from a high pressure tunnel and study on its permeation stability[J]. Journal of Water Resources and Architectural Engineering, 2015, 13(5): 212- 217. (in Chinese))

[7] 沈振中, 倪治斌, 趙堅(jiān). 水工洞室圍巖穩(wěn)定性的非連續(xù)變形分析[J]. 巖石力學(xué)與工程學(xué)報(bào), 2003, 22(增刊1): 2299- 2303. (SHEN Zhenzhong, NI Zhibin, ZHAO Jian. Stability analysis of rock mass for underground hydraulic powerhouse by DDA[J]. Chinese Journal of Rock Mechanics and Engineering, 2003, 22(Suppl1): 2299- 2303. (in Chinese))

[8] 沈振中， 張?chǎng)危?陸希， 等. 西藏老虎嘴水電站左岸滲流控制優(yōu)化[J]. 水利學(xué)報(bào), 2006, 37(10): 1230- 1234. (SHEN Zhenzhong, ZHANG Xin, LU Xi, et al. Seepage control optimization of left bank of Laohuzui hydropower station[J]. Journal of Hydraulic Engineering, 2006, 37(10): 1230- 1234. (in Chinese))

[9] 周人杰， 沈振中， 徐力群， 等. 基于三維非穩(wěn)定滲流分析的隧洞開挖地下水環(huán)境影響評(píng)價(jià)[J]. 南水北調(diào)與水利科技, 2016, 14(6): 135- 140. (ZHOU Renjie, SHEN Zhenzhong, XU Liqun, et al. Effect of tunnel excavation on groundwater environment based on three-dimensional unsteady seepage flow analysis[J]. South-to-North Water Transfers and Water Science & Technology, 2016, 14(6): 135- 140. (in Chinese))

[10] 王連廣， 劉陽. 引水隧洞三維滲流場(chǎng)分析研究[J]. 水利規(guī)劃與設(shè)計(jì), 2014(4): 52- 54. (WANG Lianguang, LIU Yang. Analysis of three-dimensional seepage field at diversion tunnel[J]. Water Resources Planning and Design, 2014(4): 52- 54. (in Chinese))

[11] 任華春， 沈振中， 張耀祖. 地下廠房圍巖穩(wěn)定性分析及支護(hù)方案優(yōu)化[J]. 水電能源科學(xué), 2010, 28(4): 96- 99. (REN Huachun, SHEN Zhenzhong, ZHANG Yaozu. Stability analysis and reinforcement scheme optimization of surrounding rock for underground powerhouse[J]. Water Resources and Power, 2010, 28(4): 96- 99. (in Chinese))

[12] 王東， 沈振中， 陶小虎. 尾礦壩滲流場(chǎng)三維有限元分析與安全評(píng)價(jià)[J]. 河海大學(xué)學(xué)報(bào)(自然科學(xué)版), 2012, 40(3): 307- 312. (WANG Dong, SHEN Zhenzhong, TAO Xiaohu. Three-dimensional finite element analysis and safety assessment for seepage field of a tailing dam[J]. Journal of Hohai University(Natural Science), 2012, 40(3): 307- 312. (in Chinese))

[13] 劉干斌， 謝康和， 施祖元， 等. 壓力隧洞襯砌-圍巖(土)相互作用研究[J]. 巖石力學(xué)與工程學(xué)報(bào), 2005， 24(14)： 2449- 2455. (LIU Ganbin, XIE Kanghe, SHI Zuyuan, et al. Interaction of surrounding rock or soil and lining of a pressure tunnel[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(14)： 2449- 2455. (in Chinese))

[14] 白明洲， 陳云， 師海. 山嶺隧道施工誘發(fā)地下水位下降環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)[J]. 鐵道工程學(xué)報(bào), 2016(1)： 5- 10, 15. (BAI Mingzhou, CHEN Yun, SHI Hai. Environmental risk assessment of underground water level falling induced by construction in mountain tunnel[J]. Journal of Railway Engineering Society, 2016(1): 5- 10, 15. (in Chinese))

[15] 蘇凱， 伍鶴皋. 水工隧洞內(nèi)水外滲耦合分析[J]. 巖土力學(xué), 2009, 30(4): 1147- 1152. (SU Kai, WU Hegao. Analysis of hydro-mechanical interaction in hydraulic tunnel with inner water exosmosis[J]. Rock and Soil Mechanics, 2009, 30(4): 1147- 1152. (in Chinese))