地質(zhì)缺陷的面板堆石壩滲流特性分析及處理措施研究

徐力群，黃柏云，陸譽(yù)婷，蔣裕豐

(1.河海大學(xué)水利水電學(xué)院，江蘇南京210098；2.江蘇省交通規(guī)劃設(shè)計(jì)院股份有限公司，江蘇南京210014；3.南京河?？萍加邢薰?，江蘇南京210098)

地質(zhì)缺陷的面板堆石壩滲流特性分析及處理措施研究

徐力群1，黃柏云1，陸譽(yù)婷2，蔣裕豐3

(1.河海大學(xué)水利水電學(xué)院，江蘇南京210098；2.江蘇省交通規(guī)劃設(shè)計(jì)院股份有限公司，江蘇南京210014；3.南京河海科技有限公司，江蘇南京210098)

新疆某水電工程地質(zhì)條件復(fù)雜，層面及構(gòu)造發(fā)育，左岸壩基存在斷層構(gòu)造，水平寬度達(dá)0.5 m且上下游貫通，左岸壩肩存在寬5 m的白云巖溶蝕帶，溶穴較為發(fā)育。根據(jù)壩址區(qū)地質(zhì)地形特性，建立了工程壩址區(qū)三維有限元滲流性態(tài)分析模型，模擬了溶蝕帶和大斷層等重要地質(zhì)構(gòu)造，采用改進(jìn)折算滲透系數(shù)法以精確模擬溶蝕帶的滲透性，分析了左岸白云巖溶蝕帶和F4斷層地質(zhì)缺陷對(duì)壩址區(qū)滲流性態(tài)的影響，研究了采用不同深度混凝土洞塞處理溶蝕帶的防滲效果，開展了斷層滲透參數(shù)敏感性分析，提出了白云巖溶蝕帶采用混凝土洞塞截滲及斷層采用混凝土帷幕灌漿的防滲處理方式。

溶蝕帶；改進(jìn)折算滲透系數(shù)法；面板堆石壩；斷層；混凝土洞塞

0 引 言

國(guó)內(nèi)外有許多由于地質(zhì)缺陷而致使水庫(kù)大壩和壩肩產(chǎn)生嚴(yán)重滲漏問題的水利水電工程，有些甚至導(dǎo)致大壩塌陷和破壞失事。例如，1928年美國(guó)圣弗朗西斯(St.Francis)重力壩由于壩基右側(cè)處存在順河向的斷層，在蓄水初期庫(kù)水滲入斷層，從而引發(fā)突然潰壩[1]；1959年法國(guó)馬爾帕賽拱壩(Malpasest arch dam)由于左岸壩肩巖體裂隙發(fā)育，蓄水后庫(kù)水滲入裂隙產(chǎn)生較大較大水壓力，造成大壩失事[2]；1976年美國(guó)Teton心墻壩蓄水后，節(jié)理裂隙中的集中滲流將與之接觸的粉土帶走致使接觸沖刷，致使初期蓄水就引起潰壩事故[3]；貓?zhí)铀募?jí)水電站水庫(kù)蓄水后，由于巖溶滲漏造成電站枯水期六分之一電量損失[4]。

國(guó)內(nèi)外學(xué)者開展地質(zhì)缺陷對(duì)大壩滲流性態(tài)的影響研究，尤其是關(guān)于裂隙等結(jié)構(gòu)面和巖溶引起的滲漏問題，提出了不同的研究方法，并給出了處理地質(zhì)缺陷的不同防滲措施。Snow[5]采用平行板模擬了裂隙張開，提出了二階張量裂隙；Long[6]提出了將復(fù)雜裂隙網(wǎng)格等效為連續(xù)介質(zhì)處理的條件和方法；柴軍瑞[7]提出了基于多層網(wǎng)絡(luò)的裂隙巖體地下水滲流有限元分析方法；陳崇希[8]建立了巖溶三種空隙介質(zhì)地下水統(tǒng)一的控制方程，提出了采用折算系數(shù)法模擬巖溶三重介；沈振中[9]通過巖溶管道和裂隙交叉試驗(yàn)，建立了基于遺傳算法的巖溶管道和裂隙交叉匯流量計(jì)算方法；史光前[10]針對(duì)隔河巖水利工程采用混凝土封堵溶洞，及對(duì)斷層進(jìn)行帷幕灌漿，解決了壩基及兩岸繞滲；Ghobadi[11]針對(duì)伊朗Shahid Abbaspour拱壩指出水庫(kù)主要通過溶洞石灰?guī)r滲漏，對(duì)溶洞進(jìn)行了壩基固結(jié)灌漿和高壓帷幕灌漿，同時(shí)對(duì)右岸斷層進(jìn)行了修建混凝土防滲墻的加固處理方式；紀(jì)偉針[12]對(duì)水布埡面板堆石壩的巖溶化程度較高的問題，提出了優(yōu)化的防滲布置措施。

圖1 壩體標(biāo)準(zhǔn)剖面及分區(qū)(單位：m)

水利水電工程中易造成滲漏問題的地質(zhì)缺陷主要包括巖溶、古河道和斷層、節(jié)理、裂隙、層間錯(cuò)動(dòng)帶等結(jié)構(gòu)面以及卸荷帶、次生夾泥層等。本文針對(duì)存在地質(zhì)缺陷中巖溶滲漏問題的新疆某水電站工程，采用改進(jìn)折算系數(shù)法，開展溶蝕帶和斷層等地質(zhì)缺陷下面板堆石壩的滲流性態(tài)研究，并由此提出最優(yōu)的工程措施，實(shí)現(xiàn)對(duì)該工程的滲流控制，保證大壩滲流安全性。

1 工程概述

新疆開都河上某水電工程壩址處為基本對(duì)稱且較完整的“V”形谷，河床覆蓋層約厚50 m，左岸巖體主要為白云巖，局部存在溶蝕集中現(xiàn)象，右岸基巖主要為英安質(zhì)、安山質(zhì)凝灰?guī)r，層面及構(gòu)造發(fā)育，巖體卸荷風(fēng)化較深，完整性較差，壩址區(qū)斷裂和裂隙構(gòu)造較發(fā)育，共發(fā)現(xiàn)20余條斷層，裂隙多成組發(fā)育，少則4～5條，多則13～20條。攔河壩為趾板建在覆蓋層上的混凝土面板壩，工程擬定正常蓄水位1 820.00 m，壩頂高程1 825.00 m，最大壩高158.0 m，壩體標(biāo)準(zhǔn)剖面及分區(qū)如圖1所示。

本工程地質(zhì)條件復(fù)雜，特別是位于高程1 740 m附近的左岸白云巖溶蝕帶的局部溶蝕地質(zhì)條件，溶蝕帶寬1～2 m，長(zhǎng)約52 m，溶洞直徑最大0.4 m，最深1.5 m，溶洞順層密集發(fā)育，洞內(nèi)有塌落現(xiàn)象，塌落巖塊和溶洞洞壁多附著重結(jié)晶方解石顆粒。左岸F4結(jié)構(gòu)面走向?yàn)镹W285°，傾向NE，傾角67°，水平寬0.5 m，順河向延伸，貫穿上下游。包含白云巖溶蝕帶及F4斷層分布的面板壩壩軸線工程地質(zhì)如圖2所示。

圖2 面板壩壩軸線工程地質(zhì)示意

2 計(jì)算原理

巖溶介質(zhì)的多樣性決定了求解巖溶滲流的復(fù)雜性，根據(jù)本工程巖溶蝕帶的地質(zhì)特性，采用改進(jìn)折算滲透系數(shù)法的三重介質(zhì)模型進(jìn)行壩址區(qū)滲流場(chǎng)性態(tài)分析[12]。改進(jìn)折算滲透系數(shù)法通過建立基于巖溶介質(zhì)的非達(dá)西滲流規(guī)律的統(tǒng)一表達(dá)式，由此反應(yīng)巖溶滲流場(chǎng)中層流和紊流型滲流共存現(xiàn)象的折算系數(shù)。如果忽略毛細(xì)現(xiàn)象、非飽和部分和蒸發(fā)的影響，考慮基于改進(jìn)折算滲透系數(shù)法的巖溶滲流基本微分方程為[13]

(1)

式中，kij為滲透張量，對(duì)于孔隙介質(zhì)為達(dá)西滲流下的滲透張量，對(duì)于巖溶介質(zhì)為折算滲透系數(shù)張量；h為水頭；Ss為貯水系數(shù)；t為時(shí)間。為其達(dá)西滲流下的滲透張量，對(duì)于巖溶介質(zhì)kij為其折算滲透系數(shù)張量

采用Galerkin逼近方法，有限元支配方程可寫為

(2)

式中，K為總滲透矩陣；P為節(jié)點(diǎn)水頭列陣；S為貯水矩陣；G為給水矩陣；F為不平衡流量矩陣。采用截止負(fù)壓法[14]對(duì)方程進(jìn)行自由面和滲流場(chǎng)解。

3 計(jì)算模型和工況

3.1 模型范圍

模型中三維坐標(biāo)系為以壩頂左壩端中心為模型坐標(biāo)原點(diǎn)，Y方向?yàn)閴屋S線方向，Z方向按實(shí)際高程。根據(jù)邊界截取的基本原則以及本工程的實(shí)際工程特點(diǎn)，確定計(jì)算模型的范圍。三維模型截取邊界范圍為模型右邊界為原點(diǎn)以右500 m(Y=-500 m)，模型左邊界為原點(diǎn)以左969 m(Y=969 m)；X方向?yàn)轫樅恿鞣较?，模型上游邊界為壩軸線上游700 m(X=-700 m)，下游邊界為壩軸線下游700 m(X=700 m)，底高程截至1 417 m，頂高程截至1 826.2 m，低于1 826.2 m的地形按實(shí)際高程考慮，高于1 826.2 m地形按削平處理。

采用自行編制的三維滲流有限元分析程序，建立了三維有限元滲流分析模型，模型較精確的模擬了地質(zhì)條件和壩體結(jié)構(gòu)，例如巖層分區(qū)、斷層、巖溶蝕帶等以及壩體材料分區(qū)、防滲墻、防滲帷幕和導(dǎo)流洞等。該滲流分析模型生成的有限元網(wǎng)格結(jié)點(diǎn)總數(shù)為62 269個(gè)，單元總數(shù)為60 177個(gè)，三維有限

元模型網(wǎng)格如圖3所示。

圖3 滲流分析計(jì)算范圍內(nèi)三維有限元模型網(wǎng)格

3.2 模型邊界

本工程計(jì)算模型邊界類型主要包括：①已知水頭邊界包括壩址區(qū)上、下游水位線以下的水庫(kù)庫(kù)岸和庫(kù)底、壩體上游坡和下游坡、河道，以及給定地下水位的截取邊界；②出滲邊界為壩址區(qū)上下游水位線以上的左、右岸山坡面，以及壩體上、下游坡面和壩頂；③不透水邊界包括不透水邊界包括模型上、下游兩側(cè)和左右岸兩側(cè)截取邊界除給定地下水位以外的部分邊界，以及模型底面。

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

根據(jù)工程地質(zhì)和水文地質(zhì)資料，壩址區(qū)巖體分為4類地層，包括中等透水層、弱透水層、相對(duì)不透水層、不透水層，同時(shí)考慮河床覆蓋層、巖溶蝕帶和F4斷層等，各巖體和地質(zhì)缺陷等介質(zhì)的等效滲透系數(shù)如表1所示。根據(jù)工程設(shè)計(jì)方案和相關(guān)工程經(jīng)驗(yàn)，壩體各料區(qū)分為面板、墊層、過渡層、主堆石區(qū)、次堆石區(qū)等，防滲系統(tǒng)包括防滲帷幕和混凝土防滲墻等，壩體各料區(qū)和防滲系統(tǒng)滲透系數(shù)如表2所示。此外擬定白云巖溶蝕帶的粗糙度為0.03 m，水的粘滯系數(shù)為1.1×10-6m2/s。

表1 巖體和地質(zhì)缺陷地層滲透系數(shù)

地層巖層分類壩基滲透系數(shù)/m·s-1備注含漂石砂卵礫石層—6.0×10-4—含礫中粗砂層—4.0×10-4—中等透水10Lu

表2 面板壩壩體各料區(qū)和防滲系統(tǒng)滲透系數(shù)

壩料分區(qū)滲透系數(shù)/m·s-1面板1.0×10-10墊層1.0×10-5過渡層1.0×10-4主堆石1.0×10-4下游堆石1.0×10-3壩基防滲帷幕3.0×10-7混凝土防滲墻1.0×10-9泄洪洞洞身混凝土襯砌1.0×10-10

3.4 計(jì)算工況

針對(duì)工程存在的左岸白云巖溶蝕帶和F4斷層等地質(zhì)缺陷，擬定7種不同方案(見表3)研究地質(zhì)缺陷對(duì)工程滲流場(chǎng)性態(tài)的影響，并進(jìn)一步分析地質(zhì)缺陷的處理措施以及優(yōu)化方案。

表3 正常運(yùn)行期滲流場(chǎng)分析計(jì)算方案

方案研究?jī)?nèi)容方案說明1無(wú)工程措施方案溶蝕帶和斷層不采取防滲措施的設(shè)計(jì)方案2溶蝕帶防滲措施1采用混凝土洞塞截滲,截滲深度75m3溶蝕帶工程措施2采用混凝土洞塞截滲,截滲深度100m4溶蝕帶工程措施3采用混凝土洞塞截滲,截滲深度125m5溶蝕帶工程措施4采用混凝土洞塞截滲,截滲深度150m6F4斷層滲透性影響分析1F4斷層滲透性增加2倍7F4斷層滲透性影響分析2F4斷層滲透性增加5倍8F4斷層滲透性影響分析3F4斷層滲透性增加10倍

4 地質(zhì)缺陷影響分析

4.1 白云巖溶蝕帶滲流性態(tài)分析

左岸白云巖溶蝕帶滲透性較強(qiáng)，白云巖溶蝕帶對(duì)滲流場(chǎng)的影響較大，壩址區(qū)地下水位等值線分布如圖4所示，沿白云巖溶蝕帶上、下游方向斷面滲流場(chǎng)位勢(shì)分布如圖5所示(圖中黑粗部分為溶蝕帶分布)。由圖4、圖5和計(jì)算結(jié)果可得，由于存在白云巖溶蝕帶，形成了沿上下游方向的集中滲漏通道，對(duì)滲流場(chǎng)有較大的影響，其中因計(jì)算工況考慮防滲帷幕的阻滲作用，因此實(shí)質(zhì)上已對(duì)白云巖溶蝕帶進(jìn)行部分防滲處理，并取得一定的防滲效果，但由于溶蝕帶超出防滲系統(tǒng)范圍，仍存在集中滲漏問題。因此，根據(jù)計(jì)算結(jié)果，建議采用混凝土洞塞形式對(duì)左岸白云巖溶蝕帶進(jìn)行進(jìn)一步防滲處理。

左岸白云巖溶蝕帶基本位于高程1 740 m附近，為滲控措施優(yōu)化分析需要，擬定4種不同截滲深度方案，對(duì)比分析其對(duì)滲流特性的影響。白云巖溶蝕帶不同截滲深度的滲流性態(tài)計(jì)算成果見表4，混凝土洞塞截滲深度為150 m的沿白云巖溶蝕帶上、下游方向斷面滲流場(chǎng)位勢(shì)圖如圖5中虛線所示。

由表4可知，采用混凝土洞塞后，壩體浸潤(rùn)面變化幅度較小，尤其在采用75 m以上截滲深度的混凝土洞塞措施后，影響幅度更甚微。左岸白云巖溶蝕帶靠近左岸壩肩，對(duì)左岸壩肩處防滲帷幕的滲透坡降有一定的影響，采用混凝土洞塞后，左壩肩防滲帷幕最大滲透坡降有一定幅度增加，并且隨著混凝土洞塞深度的增加的防滲體系的阻滲效果進(jìn)一步加強(qiáng)。左岸壩基的滲透流量受白云巖溶蝕帶的影響明顯，在僅采用防滲帷幕方案下，由于滲漏通道并未完全截?cái)?，左岸壩基滲透流量達(dá)2 096.9 m3/d，而采用75 m截滲深度的混凝土洞塞，滲透流量降為1 936.0 m3/d，滲透流量減少量達(dá)7.7%，但隨著混凝土洞塞截滲深度的增加，滲漏量并未隨截滲深度的增加呈線性減少的趨勢(shì)，阻滲效果不明顯。

表4 白云巖溶蝕帶不同截滲深度的滲流性態(tài)計(jì)算成果

工況面板后浸潤(rùn)面最高點(diǎn)位置/m左壩肩防滲帷幕最大滲透坡降左壩肩滲透流量/m3·d-111681.376.552096.921681.256.591936.031681.246.621865.541681.236.641781.651681.236.651712.4

注：左壩肩滲透流量是指左岸壩肩處向左延伸直至模型邊界，并通過沿壩軸線橫斷面的滲透流量。

圖4 壩址區(qū)地下水位等值線分布(單位：m)

4.2 F4斷層滲流性態(tài)分析

根據(jù)地質(zhì)勘測(cè)資料分析得到左岸F4斷層滲透系數(shù)約為1.0×10-6m/s，該斷層位于左岸壩肩帷幕灌漿范圍以外，需分析該地質(zhì)缺陷是否需要單獨(dú)進(jìn)行防滲處理。正常運(yùn)行工況下壩址區(qū)地下水位等值線分布如圖4，F(xiàn)4斷層上下游方向斷面滲流場(chǎng)位勢(shì)分布如圖6所示(圖中黑粗線為斷層分布，虛線為位勢(shì)分布)。由圖4、圖6和三維滲流計(jì)算結(jié)果可知，在該滲透參數(shù)下，F(xiàn)4斷層對(duì)壩址區(qū)滲流場(chǎng)影響較小。但由于F4斷層存在一定的不可預(yù)知因素，因此擬定三組不同滲透參數(shù)的斷層參數(shù)，對(duì)比分析其對(duì)壩址區(qū)滲流特性的影響，對(duì)比分析計(jì)算成果如表5所示，F(xiàn)4斷層滲透系數(shù)放大10倍后沿F4斷層上下游方向斷面滲流場(chǎng)位勢(shì)分布如圖6中實(shí)線所示。由表5可知，隨著F4斷層滲透系數(shù)的增加，壩體內(nèi)浸潤(rùn)面不斷抬高，滲透系數(shù)增加10倍后，水位高度增加1.69 m；同時(shí)因F4斷層滲透性增強(qiáng)，使得防滲帷幕的阻滲效果降低，左壩肩防滲帷幕最大滲透坡降由6.55減小為4.71；并且左岸壩基處的滲透流量增幅較為明顯，由2 906.7 m3/d增加為3 196.5 m3/d，增加幅度達(dá)9.97%。由此可見，F(xiàn)4斷層滲透性的增加，使得左岸壩基處形成了較大的滲漏通道，影響水庫(kù)的蓄水和發(fā)電能力，同時(shí)壩體內(nèi)浸潤(rùn)面抬升，進(jìn)而可能影響壩體的安全性。

圖5 白云巖溶蝕帶上下游方向斷面滲流場(chǎng)位勢(shì)分布(單位：m)

圖6 F4斷層上下游方向斷面滲流場(chǎng)位勢(shì)分布(單位：m)

表5 F4斷層不同滲透參數(shù)的滲流場(chǎng)計(jì)算成果

注：左壩基滲透流量是指原河床左岸處至左岸壩肩處，通過沿壩軸線橫斷面的滲透流量。

5 結(jié) 論

(1)通過計(jì)算分析表明，本文所采用改進(jìn)折算滲透系數(shù)法可較好的解決溶蝕帶和大斷層等地質(zhì)缺陷構(gòu)造的滲流模擬問題。

(2)左岸白云巖溶蝕帶對(duì)左壩肩滲流性態(tài)影響較大，建議將整個(gè)白云巖溶蝕帶采用混凝土洞塞截滲，截滲深度為75 m的截滲效果明顯，可滿足阻滲目的，且工程量較少，施工難度較小。

(3)F4斷層可暫不進(jìn)行處理，但由于斷層存在不可預(yù)測(cè)性，因此建議對(duì)此斷層進(jìn)行必要的混凝土灌漿形式的防滲處理，以防止左壩基發(fā)生集中滲漏的危險(xiǎn)。

(4)加強(qiáng)壩址區(qū)地質(zhì)缺陷的前期地質(zhì)勘測(cè)工作，根據(jù)不同的地質(zhì)缺陷性質(zhì)及所處位置，采取相應(yīng)的防滲處理措施，本工程防滲處理范圍和分析方法可為類似工程的溶蝕帶模擬及為溶蝕帶工程處理措施提供借鑒。

[1]馬永鋒， 生曉高. 大壩失事原因分析及對(duì)策探討[J]. 人民長(zhǎng)江， 2001， 32(10)： 53- 54， 66．

[2]LONDE P. The Malpasset Dam failure[J]. Engineering Geology， 1987， 24(12)： 295- 329．

[3]顧淦臣. 國(guó)內(nèi)外土石壩重大事故剖析—對(duì)若干土石壩重大事故的再認(rèn)識(shí)[J]. 水利水電科技進(jìn)展， 1997， 17(1)： 13- 20， 65．

[4]張興仁. 全國(guó)巖溶滲漏學(xué)術(shù)討論會(huì)總結(jié)[J]. 水力發(fā)電學(xué)報(bào)， 1985(4)： 77- 82．

[5]SNOW D T. Anisotropic Permeability of Fractured Media[J]. Water Resources Research， 1969， 5(6)： 553- 566．

[6]LONG J C S. Porous Media Equivalents for Networks of Discontinuous Fractures[J]. Water Resources Research， 1982， 18(3)： 236- 247．

[7]CHAI Junrui， LI Shouyi， LI Shouyi. Multi-level fracture network model and FE solution for ground water flow in rock mass[J]. Journal of Hydraulic Research， 2005， 43(2)： 202- 207．

[8]陳崇希. 巖溶管道一裂隙一孔隙三重空隙介質(zhì)地下水流模型及模擬方法研究[J]. 地球科學(xué)—中國(guó)地質(zhì)大學(xué)學(xué)報(bào)， 1995， 20(4)： 361- 366 ．

[9]沈振中， 陳雰， 趙堅(jiān). 巖溶管道與裂隙交叉滲流特性試驗(yàn)研究[J]. 水利學(xué)報(bào). 2008， 39(2)： 137- 145．

[10]史光前. 隔河巖水利樞紐基礎(chǔ)滲流控制[J]. 水力發(fā)電， 1993(5)： 41- 44， 72．

[11]GHOBADI M H， KHANLARI G R， DJALALY H. Seepage problems in the right abutment of the Shahid Abbaspour dam， southern Iran[J]. Engineering Geology， 2005， 82(12)： 119- 126．

[12]紀(jì)偉， 趙堅(jiān)， 沈振中， 等. 水布埡水利樞紐巖溶體內(nèi)防滲帷幕優(yōu)化布置研究[J]. 水電能源科學(xué)， 2005(2)： 36- 39．

[13]趙堅(jiān)， 賴 苗， 沈振中. 適于巖溶地區(qū)滲流場(chǎng)計(jì)算的改進(jìn)折算滲透系數(shù)法和變滲透系數(shù)法[J]. 巖石力學(xué)與工程學(xué)報(bào)， 2005， 24(8)： 1341- 1347．

[14]速寶玉， 沈振中， 趙堅(jiān). 用變分不等式理論求解滲流問題的截止負(fù)壓法[J]. 水利學(xué)報(bào)， 1996(3)： 22- 29， 35．

(責(zé)任編輯王 琪)

Study on Seepage Characteristics and Anti-seepage Measures for CFRD with Geological Defects

XU Liqun1, HUANG Baiyun1, LU Yuting2, JIANG Yufeng3

(1. College of Water Conservancy and Hydropower, Hohai University, Nanjing 210098, Jiangsu, China;2. Jiangsu Province Communications Planning and Design Institute Co., Ltd., Nanjing 210014, Jiangsu, China;3. Nanjing Hohai Technology Co., Ltd., Nanjing 210098, Jiangsu, China)

A hydropower project in Xinjiang has complex geological conditions. A fault structure (F4) passes through the upstream and downstream with 0.5 m width in the left bank of dam foundation, and a developed dolomite karst zone with 5 m width also passes through the left abutment of dam. According to geological and topographical characteristics of dam foundation, a three-dimensional finite element model is established for analyzing the seepage characteristics of dam. In this model, the fault structure and developed dolomite karst zone are simulated, and the improved converting permeability coefficient method is applied to simulate the permeability of dolomite karst zone. The influences of these geological defects on dam seepage characteristics are analyzed. The anti-seepage effects of concrete plug with different depths in dolomite karst zone are studied. The sensitivity of penetration parameter of fault is also analyzed. The anti-seepage methods by using concrete plug in dolomite karst zone and concrete curtain grouting for fault are finally proposed．

karst zone; improved converting permeability coefficient method; CFRD; fault; concrete plug

2015- 05- 05

江蘇省自然科學(xué)基金青年基金(BK2012410);國(guó)家自然科學(xué)基金項(xiàng)目(51179062)

徐力群(1983—)，男，浙江仙居人，講師，碩士生導(dǎo)師，博士，主要研究方向水工結(jié)構(gòu)工程滲流控制與優(yōu)化計(jì)算與分析．

TV641.4

A

0559- 9342(2015)12- 0048- 06