趙正鵬，江 霄，王 假

（四川省廣元水文水資源勘測(cè)局，四川 廣元 628017）

0 引言

巖土體滲流問(wèn)題作為水力學(xué)、巖土力學(xué)計(jì)算中不可或缺的組成部分，也是水利水電工程十分關(guān)心的工程問(wèn)題[1]。但是滲流計(jì)算的可靠性和滲控措施的合理性仍受控于水文地質(zhì)條件和巖土體滲透特性。由于巖體中裂隙網(wǎng)絡(luò)分布與土體多孔介質(zhì)結(jié)構(gòu)有較大差異，裂隙巖體的滲透特性往往表現(xiàn)出各向異性。若按傳統(tǒng)方式將裂隙巖體視為各向同性的多孔介質(zhì)結(jié)構(gòu)，不僅可能產(chǎn)生較大誤差，甚至?xí)?dǎo)致壩體等水工結(jié)構(gòu)偏于不安全[2～3]。因此，結(jié)合巖體自身結(jié)構(gòu)特點(diǎn)，多角度、全方位的了解其滲流特點(diǎn)顯得極為重要。

西藏地區(qū)某水電站雖然主體工程規(guī)模不大，但是工程地質(zhì)問(wèn)題復(fù)雜。壩址河床除分布“深、厚、寬、廣”的覆蓋層外，右壩肩所在的凸形山梁巖體卸荷強(qiáng)烈，鉆孔揭示壩基壩肩巖體透水率普遍在中等透水范圍，局部甚至達(dá)到強(qiáng)透水；平硐揭示除邊坡淺表部巖體存在較強(qiáng)卸荷外，深部還出現(xiàn)深卸荷問(wèn)題[4～5]。這些問(wèn)題直接關(guān)系到邊坡穩(wěn)定、壩基壩肩滲漏、相應(yīng)的工程處理措施以及工程建成后的安全運(yùn)行等諸多方面[6]。本文采用有限元方法，分析27種工況下的某水電站右壩肩卸荷巖體繞壩滲流情況。

1 概況

1.1 工程概況

某電站位于西藏林芝地區(qū)八一鎮(zhèn)多布村尼洋河干流上，是尼洋河巴河口以下河?段水電規(guī)劃的第3個(gè)梯級(jí)，距拉薩公路里程約375 km，距林芝市政府所在地八一鎮(zhèn)公路里程約25 km，距米易機(jī)場(chǎng)公路里程約75 km，川藏公路從庫(kù)壩區(qū)左岸通過(guò)，對(duì)外交通十分便利。

該電站以發(fā)電為主，設(shè)計(jì)正常蓄水位3076.00 m、死水位3074.00 m，正常蓄水位以下庫(kù)容為0.65×108m3，最大壩高27 m，裝機(jī)容量120 MW。主要建筑物由攔河大壩(壩型為砂礫石復(fù)合壩)、泄水閘和左岸引水發(fā)電系統(tǒng)組成。樞紐屬三等中型工程，主要建筑物為3級(jí)，臨時(shí)建筑物為5級(jí)，設(shè)計(jì)洪水標(biāo)準(zhǔn)為100年一遇，校核洪水標(biāo)準(zhǔn)為500年一遇。為了減少投資和縮短工期，直接利用河床覆蓋層建壩，攔河大壩采取砂礫石土工膜復(fù)合壩方案。

1.2 地質(zhì)條件

該電站所在的尼洋河流域位于南念青唐古拉山脈中段，為典型高山地形，相對(duì)高差一般在2000 m～2500 m。壩址微地貌表現(xiàn)為“峽谷”地形特征，河流流向約NE70°，主流偏右岸，平水期河水位為3055 m～3056 m，水深約3 m～4 m，谷寬158 m～187 m，河道平均坡降0.44%，河谷具典型冰蝕“U”型寬谷特征，最大谷寬約1.5 km。

壩址區(qū)主要發(fā)育河谷3級(jí)階地，壩址“峽谷”進(jìn)出口為漫灘及Ⅰ、Ⅱ級(jí)堆積階地。漫灘和江心灘頂面高出河水面0.8 m～2.6 m，Ⅰ級(jí)階地前緣高程3059 m～3061 m，高出水面6 m～8 m。Ⅱ級(jí)階地前緣高程3068 m～3071 m，高出水面14 m～17 m。Ⅲ級(jí)階地主要分布在左岸多布村一帶，前緣高程為3104 m～3082 m，高出水面28 m～45 m，而右壩肩部位主要分布在3100 m高程水渠附近。Ⅲ級(jí)階地地面起伏較大，呈壟崗狀地貌，相對(duì)高差27 m～29 m，最大34 m，其間寬緩小沖溝較發(fā)育。壩線下游地勢(shì)較平坦，表部高程3081 m～3069 m（如圖1）。壩址基巖為二長(zhǎng)花崗巖，致密堅(jiān)硬，抗風(fēng)化能力強(qiáng)，巖體呈塊狀～次塊狀結(jié)構(gòu)，完整性較好，壩基壩肩巖體淺表部以弱風(fēng)化為主，河床和岸坡弱風(fēng)化巖體厚度大致在25m以內(nèi)。河床覆蓋層物質(zhì)主體為砂卵礫石，為強(qiáng)透水層?；野咨?，似斑狀結(jié)構(gòu)，塊狀構(gòu)造。除坡頂和山脊表部有薄層殘積土和全風(fēng)化外，右壩肩巖體基本不存在全風(fēng)化。

圖1 壩址地形地貌特征(下游拍攝)右壩肩基巖為燕山晚期～喜馬拉雅期二長(zhǎng)花崗巖(γ5-63)，

圖2 PD04現(xiàn)場(chǎng)照片(a)硐深64 m高傾角順河向裂隙張開(kāi)狀態(tài)；(b)坡外側(cè)端部基/覆交界面

根據(jù)地表測(cè)繪、平硐勘探、鉆孔揭示及物探波速資料，以壩走向剖面PD04平硐為代表，現(xiàn)場(chǎng)照片見(jiàn)圖2?？蓪⒂野缎逗蓭r體水平分帶為：0～30 m(PD04為27 m)為強(qiáng)卸荷(大致對(duì)應(yīng)強(qiáng)風(fēng)化)，該段巖體裂隙開(kāi)度普遍大于1 cm，裂隙內(nèi)普遍充填次生黃砂、巖屑、角礫物質(zhì)，局部有架空且有明顯透風(fēng)特征，硐壁潮濕，降雨后硐壁有滴水；30 m～45 m(PD04為27 m～42 m)為弱卸荷(大致對(duì)應(yīng)弱風(fēng)化)，巖體裂隙開(kāi)度一般僅約幾毫米，裂隙內(nèi)充填次生黃砂、巖粉物質(zhì)，硐壁潮濕；45 m～175 m(PD04為42 m～174 m)深卸荷(大致對(duì)應(yīng)微風(fēng)化)，175m以后為正常新鮮巖體。

2 模型建立

右壩肩以及下伏于河床覆蓋層之下的基巖滲透分區(qū)，主要依據(jù)平硐裂隙統(tǒng)計(jì)和鉆孔壓水及示蹤劑法滲透測(cè)試成果，分為強(qiáng)透水、中等透水和弱透水，各滲透分區(qū)大致對(duì)應(yīng)岸坡正常卸荷帶、深卸荷帶和新鮮巖體[7]。三維滲流計(jì)算時(shí)裂隙巖體的滲透特性和滲透張量按各向異性滲透介質(zhì)考慮。壩基河床覆蓋層各介質(zhì)分區(qū)依據(jù)鉆孔資料，且覆蓋層按各向同性滲透介質(zhì)考慮[8]。由于右岸山體壩頂3080 m以上高程位于庫(kù)水位3078 m以上，對(duì)繞壩肩滲流沒(méi)有任何影響，因此，滲流計(jì)算模型不包括右岸3080 m以上部分。而右岸近壩附近的岸坡覆蓋層為強(qiáng)透水層，從工程安全角度考慮，庫(kù)水透過(guò)岸坡覆蓋層作用在基/覆交界面上的水頭均按庫(kù)水位設(shè)置。

圖3 三維滲流計(jì)算模型網(wǎng)格圖

三維滲流計(jì)算模型中，X軸與壩軸線平行，指向右岸為正；Y軸垂直向上為正；Z軸垂直壩軸線，且指向上游為正。上下游范圍從壩坡坡腳分別向下游和上游取200 m；垂直向下取至2650 m高程；順壩軸線方向從PD4硐口向坡內(nèi)延伸約307 m，而從PD4硐口向左岸延伸約330 m。計(jì)算模型包含單元總數(shù)27782個(gè)、節(jié)點(diǎn)總數(shù)28901個(gè)，劃分網(wǎng)格示意圖見(jiàn)圖3。

該電站壩體壩基防滲采用混凝土防滲墻接土工布，壩肩采用防滲帷幕，且懸掛式防滲墻墻底還設(shè)帷幕灌漿。結(jié)合設(shè)計(jì)文件并參考有關(guān)工程資料，混凝土防滲墻寬度取80 cm、深度60 m，計(jì)算分析時(shí)混凝土防滲墻的滲透系數(shù)取1×10-7cm/s=8.64×10-5m/d。防滲帷幕寬度取2.5 m，其中基巖中防滲帷幕按1 Lu控制，即滲透系數(shù)約1×10-5cm/s=8.64×10-3m/d。

3 敏感性分析

對(duì)位于防滲墻底部的覆蓋層中防滲帷幕，根據(jù)工程經(jīng)驗(yàn)，其防滲效果要達(dá)到1 Lu較為困難，因此，滲流計(jì)算中對(duì)覆蓋層中防滲帷幕按1 Lu～10 Lu控制進(jìn)行敏感性分析。另外，對(duì)防滲帷幕水平深度和垂直深度也進(jìn)行多方案的敏感性分析，防滲帷幕水平深度和垂直深度關(guān)系見(jiàn)圖4；對(duì)無(wú)任何防滲措施下的滲流進(jìn)行對(duì)比分析；對(duì)裂隙巖體按各向同性和各向異性的滲流進(jìn)行對(duì)比分析。

圖4 防滲帷幕和防滲墻空間示意圖

3.1 防滲帷幕水平深度敏感性分析

在敏感性計(jì)算時(shí)，防滲墻墻底的帷幕灌漿深度為60 m且不變，而壩肩防滲帷幕水平深度變化范圍在52.52 m～199.25 m，分別對(duì)應(yīng)平硐PD4強(qiáng)透水和中等透水的底界。根據(jù)三維滲流計(jì)算統(tǒng)計(jì)獲得不同部位的滲流量，見(jiàn)表1。

表1 防滲帷幕水平深度滲流敏感性計(jì)算成果 單位：m3

計(jì)算結(jié)果表明，即使設(shè)計(jì)為懸掛式防滲墻，但是相比無(wú)任何防滲措施下的滲流量來(lái)講，由于防滲墻、防滲帷幕和土工布的阻滲作用，有防滲措施下的總滲流量顯著減小，且無(wú)防滲措施下絕大部分滲流量主要穿過(guò)壩體。當(dāng)水平帷幕深度穿過(guò)岸坡正常卸荷區(qū)后，即使增大水平帷幕深度，對(duì)降低繞壩肩滲流量的效果并不顯著。顯然，水平帷幕防滲深度設(shè)在岸坡正常卸荷底界是合適的(方案7大致對(duì)應(yīng)平硐PD4硐深54 m附近)，方案7的總滲流量大致為無(wú)任何防滲措施下的總滲流的5.45%，防滲效果非常顯著。

由于壩肩水平防滲深度的變化，對(duì)應(yīng)防滲帷幕軸線剖面的中等透水區(qū)(大致相當(dāng)于PD4揭示的岸坡深卸荷帶)和新鮮巖體弱透水區(qū)范圍也相應(yīng)變化，因此，透過(guò)防滲帷幕、中等透水區(qū)和微新巖體弱透水區(qū)的滲流量也隨透過(guò)面積的大小呈規(guī)律性變化。但是，透過(guò)防滲墻和墻低防滲帷幕之下的河床覆蓋層的滲流量變化很小。

由于壩基壩肩強(qiáng)透水區(qū)全部被防滲系統(tǒng) (包括壩基防滲墻、墻底垂直防滲帷幕和壩肩水平防滲帷幕)所封堵，庫(kù)水滲漏實(shí)質(zhì)上已不存在透過(guò)強(qiáng)透水區(qū)巖體，而被透過(guò)防滲系統(tǒng)斷面所替代，方案0～7的強(qiáng)透水區(qū)滲漏量均為0.0。方案8計(jì)算假定為無(wú)任何防滲措施，根據(jù)Darcy定律，庫(kù)水是沿滲流路徑最短的通道向下游滲漏。由于右岸強(qiáng)、弱卸荷巖體和地表覆蓋層為強(qiáng)透水介質(zhì)，大部分滲流均沿覆蓋層和卸荷巖體透過(guò)，因此，右岸側(cè)覆蓋層和卸荷巖體繞壩滲漏損失量約為6萬(wàn)m3/d。

3.2 覆蓋層防滲帷幕滲透性敏感性分析

該電站壩高僅27 m，壩前壩后水頭差僅25 m，設(shè)計(jì)采取防滲墻預(yù)留孔方式對(duì)墻底覆蓋層和透水基巖進(jìn)行帷幕灌漿來(lái)增大滲流路徑，從而減緩懸掛式防滲系統(tǒng)的滲流量。對(duì)于基巖來(lái)講，防滲帷幕效果控制在1 Lu相對(duì)來(lái)說(shuō)較易實(shí)現(xiàn)。但是對(duì)于深厚覆蓋層來(lái)講，防滲帷幕效果要控制在1 Lu卻較為困難。為此，對(duì)河床覆蓋層中防滲帷幕效果進(jìn)行敏感性分析，且河床覆蓋層中防滲帷幕的滲透性在1 Lu～10 Lu范圍內(nèi)。根據(jù)三維滲流計(jì)算統(tǒng)計(jì)獲得不同部位的滲流量見(jiàn)表2。

表2 覆蓋層中防滲帷幕效果滲流敏感性計(jì)算成果 單位：m3

計(jì)算結(jié)果表明，覆蓋層中帷幕效果隨透水性的增加、滲流量也增加，并導(dǎo)致總滲流量也相應(yīng)增加。但是即使覆蓋層中帷幕效果為10 Lu時(shí)(即表2中方案7 Lu10)，總滲流量也僅大致為無(wú)任何防滲措施下的總滲流的6.50%，防滲效果仍較為顯著，且施工中控制覆蓋層中防滲帷幕的效果達(dá)到10 Lu以內(nèi)也是較為容易。隨覆蓋層中帷幕效果隨透水性的變化，庫(kù)水滲流路徑對(duì)應(yīng)的滲流梯度也相應(yīng)變化，從而使得各部位(包括微新、弱卸荷、防滲墻墻底帷幕之下的河床覆蓋層、防滲墻以及防滲帷幕等)的滲流量也相應(yīng)呈規(guī)律性變化。

方案0～7，由于壩基壩肩強(qiáng)透水區(qū)全部被防滲系統(tǒng)(包括壩基防滲墻、墻底垂直防滲帷幕和壩肩水平防滲帷幕)所封堵，庫(kù)水滲漏實(shí)質(zhì)上已不存在透過(guò)強(qiáng)透水區(qū)巖體，而被透過(guò)防滲系統(tǒng)斷面所替代，因此，方案0～7的強(qiáng)透水區(qū)滲漏量均為0.0。方案8為無(wú)任何防滲措施，壩體和壩基河床覆蓋層為強(qiáng)透水介質(zhì)，大部分滲流均沿壩體和壩基河床覆蓋層透過(guò)，因此，方案8中透過(guò)壩基壩肩巖體的滲流量較方案0～7有所減小。

3.3 防滲帷幕垂直深度敏感性分析

上文已初步論證防滲墻底覆蓋層中防滲帷幕透水性較大時(shí)，總滲流量也不大，說(shuō)明防滲墻和右壩肩防滲帷幕的效果顯著。在此分析墻底防滲帷幕的深度對(duì)滲流的敏感程度，以論證是否可以取消防滲墻墻底的垂直防滲帷幕。同樣根據(jù)三維滲流計(jì)算統(tǒng)計(jì)獲得不同部位的滲流量見(jiàn)表3。

表3 防滲墻底防滲帷幕深度滲流敏感性計(jì)算成果 單位：m3

計(jì)算結(jié)果表明，防滲墻底帷幕深度減小后，總滲流量的趨勢(shì)是逐漸增大。但是，當(dāng)取消墻底防滲帷幕(仍保留右岸壩肩基巖防滲帷幕)，總滲流量大致也僅為無(wú)任何防滲措施下的總滲流的7.42%，防滲效果并未明顯減弱。這說(shuō)明壩基防滲墻顯著地削減了壩基滲流量，而防滲墻墻底的灌漿帷幕對(duì)減少壩基滲流量并不顯著。對(duì)60 m深的防滲墻墻底河床覆蓋層再進(jìn)行深達(dá)60 m的帷幕灌漿，不但施工工期會(huì)延長(zhǎng)、工程投資會(huì)增大，而且防滲效果也難以保障。鑒于未來(lái)蓄水之后，近壩水庫(kù)死水位以下將逐漸沉積湖相紋泥層，其滲透性較河床中的砂卵礫石層、砂層等將明顯要低，會(huì)逐漸起到天然鋪蓋效果，從而進(jìn)一步降低壩基和繞壩肩滲流。因此，若比較分析表3中方案7-7的總滲流量未超過(guò)尼洋河平均徑流量的3%(該值是規(guī)范中控制指標(biāo))，則壩基防滲墻下不采取帷幕灌漿是可行的。

另外，隨防滲墻墻底帷幕深度的變化，庫(kù)水滲流路徑對(duì)應(yīng)的滲流梯度也相應(yīng)變化，從而使得各部位(包括微新、岸坡深卸荷中等透水區(qū)、防滲墻墻底帷幕之下的河床覆蓋層、防滲墻以及防滲帷幕等)的滲流量也呈規(guī)律性變化。

圖5 方案7 Lu5三維滲流水位等勢(shì)圖

圖6 方案7 Lu5三維滲流梯度

方案7 Lu5三維滲流水位等勢(shì)圖見(jiàn)圖5，三維滲流梯度見(jiàn)圖6。各方案的水位等勢(shì)特征符合該電站壩基壩肩滲流規(guī)律。滲流梯度較大的部位均在防滲墻、帷幕周圍和滲透性差異較大的覆蓋層各層交界面部位。壩肩巖體滲流梯度沒(méi)有超過(guò)允許滲流梯度；壩基覆蓋層在防滲墻和防滲帷幕附近局部有超過(guò)允許滲流梯度的情況，但因位于較深部位的圍壓效應(yīng)，故不會(huì)發(fā)生滲透破壞；而防滲墻和防滲帷幕滲流梯度雖然較大，但是防滲墻和防滲帷幕因允許滲流梯度較高，因此，防滲墻和防滲帷幕不會(huì)發(fā)生滲透擊穿破壞。

圖7 方案7 Lu5壩基防滲墻滲流梯度

圖8 方案7 Lu5壩基壩肩防滲帷幕滲流梯度

圖9 方案7-7壩基防滲墻滲流梯度

圖10 方案7-7壩肩防滲帷幕滲流梯度

方案7 Lu5壩基壩肩防滲墻、防滲帷幕滲流梯度分別見(jiàn)圖7、圖8；方案7-7壩肩防滲墻、防滲帷幕滲流梯度分別見(jiàn)圖9、圖10。由圖可知，方案7-7與方案7 Lu5差異較大的部位主要在防滲墻。由于方案7 Lu5墻底防滲帷幕深度大，阻滲效果更顯著，因此，防滲墻和防滲帷幕附近的壩基覆蓋層和壩肩卸荷巖體的滲流路徑更長(zhǎng)，滲流梯度相對(duì)小些。但是總體來(lái)講，除右壩肩部位的靠防滲墻和防滲帷幕附近滲流梯度較大外，其它部位均較小。因此，兩種方案壩基覆蓋層和壩肩巖體、防滲墻和防滲帷幕發(fā)生滲透的區(qū)域較小，水庫(kù)蓄水后壩基壩肩防滲體系對(duì)應(yīng)發(fā)生滲透變形破壞的程度較小。但是需注意，這些局部靠近地表附近區(qū)域存在發(fā)生淺層滲透破壞跡象。

3.4 壩肩巖體各向同性和各向異性滲流對(duì)比分析

為了說(shuō)明各向同性和各向異性滲流的差異，論證多布電站右壩肩裂隙巖體按各向異性滲流問(wèn)題處理的合理性，對(duì)前面的方案7和方案7-7中的基巖進(jìn)行各向同性滲流計(jì)算。為便于比較各向同性和各向異性滲流的差異，將計(jì)算成果列于同一表中，見(jiàn)表4。

表4 壩肩巖體各向同性和各向異性滲流計(jì)算成果對(duì)比 單位：m3

計(jì)算結(jié)果表明，就各向同性和各向異性滲流對(duì)比方案來(lái)講，如方案7和7TX，雖然總滲流量差異并不顯著，但是各滲透區(qū)尤其是基巖(對(duì)應(yīng)的微新巖體微透水區(qū)、岸坡深卸荷弱透水區(qū)和岸坡正常卸荷強(qiáng)透水區(qū)基巖)中的滲流量差異還是較為顯著的。結(jié)合壩址基巖中裂隙發(fā)育特征和平硐揭露的卸荷裂隙充填河湖相物質(zhì)的規(guī)律可知，用裂隙巖體各向異性滲流模型更能反映右壩肩巖體垂直向和順?biāo)肯驖B透性強(qiáng)、垂河向滲透性弱的特征。因此，對(duì)右壩肩巖體按各向異性滲流問(wèn)題能更好地反映未來(lái)蓄水后繞壩肩滲流特征。

4 結(jié)論

該水電站右壩肩卸荷巖體繞壩滲流是本項(xiàng)目的關(guān)鍵問(wèn)題，一定程度決定了本工程的順利開(kāi)展與否。通過(guò)本文分析，得出以下結(jié)論：

（1）盡管右壩肩巖體完整性較好，但是巖體卸荷較強(qiáng)，尤其是NE向和NW向高傾角裂隙普遍存在卸荷拉張?zhí)卣?，是控制水?kù)蓄水后繞壩肩滲流的主要通道。

（2）通過(guò)滲流計(jì)算，右岸卸荷巖體繞壩滲漏損失量約為6萬(wàn)m3/d。右壩肩水平防滲深度大致布置在距PD4硐口約52 m為相對(duì)不透水層區(qū)段 (大致對(duì)應(yīng)弱卸荷巖體以外10 m左右的微新巖體區(qū))，將有效阻止庫(kù)水繞壩肩滲流。數(shù)值計(jì)算成果與地質(zhì)上分析推薦、定性界定的水平防滲帷幕深度幾乎一致。

（3）基于滲流敏感性計(jì)算成果，右壩肩卸荷巖體水平帷幕灌漿深度可大致控制在距PD4硐口約50 m～60 m的深度。垂直防滲帷幕深度應(yīng)在剖面弱卸荷帶以下10 m左右。