，

(成都市勘察測(cè)繪研究院，四川成都610081)

中國(guó)西南地區(qū)降雨量充沛，但多為丘陵地區(qū)，山河相間，地層巖性導(dǎo)致土壤蓄水能力差，造成西南山區(qū)區(qū)域性缺水。山區(qū)中小型水庫(kù)可有效解決缺水問(wèn)題，同時(shí)發(fā)揮防洪、發(fā)電的作用。然而水庫(kù)的修建也改變了地質(zhì)環(huán)境，蓄水造成大量的岸坡失穩(wěn)，最直接的影響是滑坡失穩(wěn)造成的堵江問(wèn)題，改變水庫(kù)庫(kù)容，形成涌浪，危及人民生命財(cái)產(chǎn)安全，影響山區(qū)中小型水庫(kù)的調(diào)度和使用。研究滑坡堵江的案例也逐年增長(zhǎng)，但研究方法較為單一，通常都是采用經(jīng)驗(yàn)法[1-2]。滑坡堵江數(shù)值模擬的研究常用FLAC3D分析滑坡體的變形特征，或用UDEC模擬滑坡的破壞過(guò)程，預(yù)測(cè)堰塞壩的高度[3-4]。有學(xué)者用ABAQUS/Explicit動(dòng)力有限元分析滑坡的堵江過(guò)程[5]。也有學(xué)者選用PFC顆粒流軟件模擬土質(zhì)滑坡的破壞過(guò)程[6-8]。本文采用經(jīng)驗(yàn)法和PFC2D顆粒流程序預(yù)測(cè)某中型水庫(kù)黃角樹(shù)滑坡的堵江高度，為此類(lèi)問(wèn)題的研究提供參考和借鑒。

1 工程概況[9-11]

黃石盤(pán)水庫(kù)位于四川省南江縣恩陽(yáng)河左岸，河水最大深度約9 m，水庫(kù)正常蓄水位372.0 m，總庫(kù)容11 233×104m3，死水位357.0 m，汛期限制水位357.0 m。黃角樹(shù)滑坡為一古滑坡，覆蓋層主要為第四系殘坡積、崩坡積碎石土，滑坡堆積塊碎石土等，基巖為侏羅系蓬萊鎮(zhèn)組砂泥巖互層。通過(guò)調(diào)查及對(duì)鉆孔揭露的地層分析，古滑帶為基覆交界面，古滑坡處于穩(wěn)定狀態(tài)，但前緣堆積體向臨空方向發(fā)生了變形，鉆孔勘察確定其存在潛在滑帶，堆積體平均厚度約8 m，分布在高程353～385 m處，變形方向?yàn)镾W64°，根據(jù)現(xiàn)場(chǎng)調(diào)查滑坡體地表裂隙分布情況，分析得出變形體的破壞模式為分塊下滑，將變形體分為a、b、c、d區(qū)，分別對(duì)應(yīng)a、b、c、d剖面，見(jiàn)圖1。

圖1 前緣變形體裂隙分布

2 堵江預(yù)測(cè)分析

2.1 經(jīng)驗(yàn)法預(yù)測(cè)

常用的堵江預(yù)測(cè)方法為滑距推算法，選擇河谷寬度相對(duì)較小的b剖面進(jìn)行分析計(jì)算，各參數(shù)取值及計(jì)算結(jié)果見(jiàn)表1。

表1 滑距推算法堵江預(yù)測(cè)

注：V——滑坡體積；L1——滑坡前緣到江對(duì)岸的距離；B——河流的寬度；Lmax——最大滑距；f——?jiǎng)幽Σ料禂?shù)；h1——滑體的平均厚度；H——滑坡后緣最高點(diǎn)和滑程上計(jì)算點(diǎn)的垂直落差；Hd——堵江高度

2.2 顆粒流PFC2D預(yù)測(cè)

PFC2D顆粒流是依據(jù)顆粒之間的碰撞和摩擦產(chǎn)生的力和力矩計(jì)算顆粒間的相對(duì)運(yùn)動(dòng)，沒(méi)有固定的本構(gòu)模型，完全按照牛頓第二定律計(jì)算運(yùn)動(dòng)的形態(tài)，可以給不同的區(qū)域賦予不同的參數(shù)，利用此功能可以模擬坡體內(nèi)的穩(wěn)定水位，在水位線以下的區(qū)域顆粒間的連接參數(shù)強(qiáng)度較低，顆粒間可以破裂，相對(duì)于UDEC和3DEC的塊體不能破裂，PFC2D用于模擬土質(zhì)滑坡具有相對(duì)更好的效果，然而模擬巖制滑坡的破壞效果不如前二者。

2.2.1模型建立

根據(jù)四川省水利院提供的恩陽(yáng)河鳳儀鄉(xiāng)段的洪峰資料，按照P=1%的洪峰概率來(lái)計(jì)算，1965年從9月4日14時(shí)至次日14時(shí)，水位由371 m下降至358 m，最快的下降速率為13 m/d。選擇該工況條件計(jì)算各塊體的穩(wěn)定性，通過(guò)Geo-studio軟件計(jì)算分析，從SEEP中導(dǎo)出地下水位線，以模擬過(guò)程中，滑坡穩(wěn)定性系數(shù)最低時(shí)刻的地下水水位作為最危險(xiǎn)水位，在PFC2D建模中以此為界區(qū)分飽和土體和非飽和土體的參數(shù)值，模擬變形體失穩(wěn)堵江的堆積形態(tài)。

本次計(jì)算模型中，碎石土顆粒采用半徑R=0.25～0.75 m，服從高斯分布，平均半徑0.5 m；考慮到變形體是含塊石的碎石土，因此在剖面中隨機(jī)選取位置加入適當(dāng)?shù)膲K石模型，塊石模型用半徑R=0.1～0.3 m的更小顆粒充填，并賦予較高的黏結(jié)強(qiáng)度參數(shù)。

其中a剖面模型顆粒數(shù)為4 417，碎石土顆粒數(shù)為3 439；其中b剖面模型顆粒數(shù)為7 287，碎石土顆粒數(shù)為6 041；其中c剖面模型顆粒數(shù)為6 940，碎石土顆粒數(shù)為6 126；其中d剖面模型顆粒數(shù)為3 471，碎石土顆粒數(shù)為2 266。

根據(jù)導(dǎo)出的地下水水位建立PFC2D模型，見(jiàn)圖2，圖中黃色部分為非飽和帶，青色部分為地下水位以下飽和帶，紅色部分為塊石，d剖面模型包括既有抗滑樁，模擬比例為1∶500。

a) a剖面

b) b剖面

c) c剖面

d) d剖面圖2 滑坡體不同剖面PFC2D模型

2.2.2參數(shù)取值[12-16]

顆粒流軟件的參數(shù)采用的是顆粒之間的微觀參數(shù)，在實(shí)際的巖土體中并不存在這種參數(shù)，但兩者之間存在一定的聯(lián)系，目前尚沒(méi)有成熟的理論和方法確定宏觀參數(shù)和微觀參數(shù)間的關(guān)系，本次模擬參數(shù)采用試算的方法獲取，以SEEP導(dǎo)出的臨界水位線作為參數(shù)分界點(diǎn)，水位線以上為非飽和土體，以下為飽和土體，分別賦予不同的強(qiáng)度參數(shù)，通過(guò)計(jì)算其位移量，找到最合適的參數(shù)組合，當(dāng)模擬結(jié)果為微小位移且不發(fā)生大面積破壞時(shí)，該參數(shù)即認(rèn)定為本次模擬的微觀參數(shù)取值。再用其進(jìn)行最危險(xiǎn)工況下的模擬。本次模擬對(duì)微觀參數(shù)進(jìn)行適當(dāng)?shù)臉?biāo)定見(jiàn)表2。

表2 微觀力學(xué)參數(shù)取值

2.2.3堆積形態(tài)模擬結(jié)果

當(dāng)水庫(kù)蓄水之后死水位為357 m，高于原河道的353 m水位，因此，若變形體失穩(wěn)造成的堰塞壩最低高程超過(guò)357 m，則會(huì)影響到水庫(kù)的正常運(yùn)行，屬于完全堵江，若最低高程小于357 m，則屬于不完全堵江。庫(kù)水位13 m/d下降工況下的模擬結(jié)果見(jiàn)圖3。

a) a剖面

b) b剖面

c) c剖面

d) d剖面圖3 庫(kù)水位13 m/d下降工況下各斷面模擬結(jié)果

在該工況下，模擬結(jié)果顯示4個(gè)區(qū)域均形成SW岸低NE高的堰塞壩，以b區(qū)最為突出，b區(qū)的模擬結(jié)果則是NE岸高，SW岸堆積形態(tài)較平緩的堰塞壩，其他3個(gè)區(qū)域均為坡度較陡的堰塞壩。變形體各區(qū)根據(jù)PFC2D的顆粒位移量計(jì)算結(jié)果見(jiàn)表3。

表3 顆粒流模擬結(jié)果

模擬分析的結(jié)果屬于不完全堵江，一次形成了4座暗壩，形成順序?yàn)閏、a、b、d區(qū)，見(jiàn)圖10。

圖4 庫(kù)水位13 m/d下降工況堵江預(yù)測(cè)

2.3 綜合評(píng)價(jià)

選用經(jīng)驗(yàn)法預(yù)測(cè)滑坡堵江具有較高的效率，僅需要獲取各項(xiàng)計(jì)算數(shù)值即可，利用公式計(jì)算，用時(shí)不超過(guò)20 min，且具有大量的工程實(shí)例作為參考和佐證；而PFC2D顆粒流模擬則需要先建立模型，通過(guò)反演試算參數(shù)的方法確定其參數(shù)，最后模擬計(jì)算其破壞過(guò)程，前后所需的時(shí)間較多，本次模擬使用內(nèi)存4 G的筆記本電腦，僅模擬一次堵江的全過(guò)程尚需2 h，因此時(shí)間效率上不如經(jīng)驗(yàn)法。經(jīng)驗(yàn)法的預(yù)測(cè)堵江高度結(jié)果為351.5 m，顆粒流模擬結(jié)果為353 m，相差不大，且顆粒流模擬結(jié)果更加保守，從保護(hù)人民生命財(cái)產(chǎn)角度考慮，其安全度更高，說(shuō)明顆粒流PFC2D用于模擬山區(qū)中小型滑坡堵江，具有一定的合理性及適用性，為解決相關(guān)問(wèn)題提供參考。

3 結(jié)語(yǔ)

本文研究黃角樹(shù)滑坡前緣變形體的破壞模式，通過(guò)滑距推算法和顆粒流PFC2D法預(yù)測(cè)前緣變形體在水庫(kù)運(yùn)行工況下失穩(wěn)破壞造成的堵江高度，得出以下幾點(diǎn)結(jié)論。

a) 經(jīng)驗(yàn)法和顆粒流模擬法預(yù)測(cè)黃角樹(shù)滑坡前緣變形體失穩(wěn)堵江形成堰塞壩的最低高度分別為351.5、353 m，結(jié)果較為相近，均屬于不完全堵江，使黃石盤(pán)水庫(kù)造成庫(kù)容減少。

b) 顆粒流PFC2D模擬預(yù)測(cè)山區(qū)中小型水庫(kù)滑坡堵江高度，具有一定的適用性，模擬效果較經(jīng)驗(yàn)法更為直觀，對(duì)于巖質(zhì)邊坡則不宜使用，且無(wú)法模擬動(dòng)態(tài)水位是其不足之處，但該方法對(duì)于解決同類(lèi)問(wèn)題具有一定的參考和借鑒意義。

c) 顆粒流模型對(duì)計(jì)算機(jī)要求高，效率不如經(jīng)驗(yàn)法，模擬結(jié)果較經(jīng)驗(yàn)法偏保守，對(duì)治理方案設(shè)計(jì)要求高，從保護(hù)人民生命財(cái)產(chǎn)安全的角度考慮，其較為保守的模擬結(jié)果有一定的益處。