庫(kù)岸斜坡失穩(wěn)及其涌浪災(zāi)害風(fēng)險(xiǎn)分析

羅冠枝，高文偉，王國(guó)衛(wèi)，蔣明光，陳全明

(1.湖南省地質(zhì)環(huán)境監(jiān)測(cè)總站，湖南 長(zhǎng)沙 410007；2.中國(guó)科學(xué)院地質(zhì)與地球物理研究所，北京 100029)

0 引言

庫(kù)岸斜坡失穩(wěn)[1-2]及其引發(fā)涌浪問(wèn)題是庫(kù)岸斜坡災(zāi)害的一項(xiàng)重要研究?jī)?nèi)容，庫(kù)岸斜坡失穩(wěn)會(huì)導(dǎo)致巖土體快速入水激發(fā)涌浪，不僅危害水庫(kù)及大壩，而且對(duì)庫(kù)岸附近的建筑物及居民生命安全構(gòu)成威脅。國(guó)內(nèi)外已發(fā)生多起重大庫(kù)岸斜坡失穩(wěn)引發(fā)滑坡災(zāi)害的案例，造成了嚴(yán)重人員傷亡和財(cái)產(chǎn)損失，因此，開(kāi)展庫(kù)岸斜坡失穩(wěn)及其涌浪災(zāi)害風(fēng)險(xiǎn)分析研究對(duì)指導(dǎo)庫(kù)岸斜坡防災(zāi)減災(zāi)具有重要意義[3-7]。本文以湖南省麻陽(yáng)縣大水沖水庫(kù)滑坡為例，探討庫(kù)岸斜坡失穩(wěn)以及涌浪災(zāi)害風(fēng)險(xiǎn)分析預(yù)測(cè)，對(duì)于庫(kù)岸斜坡失穩(wěn)涌浪的預(yù)測(cè)、庫(kù)岸邊坡治理工程的實(shí)施、水庫(kù)運(yùn)營(yíng)、下游居民的選址等提供參考依據(jù)，具有非常重要的現(xiàn)實(shí)意義。

1 大水沖水庫(kù)滑坡概況

1.1 基本情況

滑坡位于湖南省麻陽(yáng)縣文昌閣鄉(xiāng)羅家沖村3組大水沖水庫(kù)旁(圖1)，地理坐標(biāo)為： 東經(jīng)109°35′24″(N)，北緯27°37′30″(E)。2011年2月，因連降大雨，斜坡發(fā)生明顯變形，坡頂至前緣范圍內(nèi)出現(xiàn)大量鼓脹裂縫，滑坡?lián)p毀房屋1棟，直接經(jīng)濟(jì)損失15萬(wàn)余元。至2013年4月，裂縫又進(jìn)一步加寬加深趨勢(shì)。目前，中部和前緣居民已搬遷，只剩下后緣西晃山林場(chǎng)辦公樓和宿舍。

圖1 大水沖水庫(kù)滑坡全貌Fig.1 Landslide panorama of Dashuichong reservoir

1.2 基本特征

圖2 大水沖水庫(kù)滑坡地質(zhì)縱剖面示意圖Fig.2 Schematic diagram geological profile of landslide on dashuichong reservoir

滑坡區(qū)屬侵蝕剝蝕丘陵地貌，滑坡后緣高程406.6～409.7 m，前緣高程351.0～350.5 m，相對(duì)高差53.5～56.2 m?；缕矫嫘螒B(tài)呈舌狀，縱向剖面形態(tài)為前緣陡、中部稍緩、后緣陡，橫向剖面形態(tài)為“V”形，地形坡度為25°～40°，局部坡度達(dá)55°(圖2)?；滦遍L(zhǎng)約150 m，寬約130 m，滑坡體厚8.8～15.7 m，平均厚12.3 m，面積約1.95×104m2，體積約2.3985×105m3，規(guī)模等級(jí)為中型?；轮骰较?5°。

滑坡物質(zhì)組成是第四系(Qel+dl)黃褐色殘坡積碎石土，滑體表層為6.90～9.00 m土石混合體層，其中碎石含量占45%～70%，粒徑一般為10～50 mm，最大為125 mm，充填約30%的黏性土，碎石呈棱角狀，成分為板巖，該層分布于整個(gè)坡面(圖3)?；鶐r為青白口系五強(qiáng)溪組(Qb2w)灰褐色強(qiáng)風(fēng)化板巖(圖4)，其下為強(qiáng)風(fēng)化板巖，主要礦物成分已強(qiáng)烈風(fēng)化，節(jié)理裂隙極發(fā)育。根據(jù)勘察實(shí)驗(yàn)，各地層巖土參數(shù)見(jiàn)表1。

表1 巖土的物理力學(xué)參數(shù)Table 1 Physical and mechanical parameters of rock and soil

圖3 土石混合體Fig.3 Soil-rock mixture

圖4 灰褐色強(qiáng)風(fēng)化板巖Fig.4 Grey-brown strong weathered slate

根據(jù)調(diào)查資料，2011年2月，因連續(xù)突降暴雨，該斜坡發(fā)生明顯變形破壞，前緣水庫(kù)邊竹林內(nèi)出現(xiàn)大量鼓脹，中部居民房屋出現(xiàn)大量變形破壞，后緣坡頂西晃山林場(chǎng)院內(nèi)出現(xiàn)明顯拉張裂縫，后緣裂縫最初寬10～30 mm，可見(jiàn)深度0.1～0.3 m，最長(zhǎng)達(dá)70余米，滑坡變形見(jiàn)圖5～圖10。

圖5 前緣竹林歪斜Fig.5 The bamboo forest is skewed

圖6 前緣出現(xiàn)大量鼓脹變形Fig.6 A lot of bulging deformation appeared at front

圖7 中部房屋墻體出現(xiàn)開(kāi)裂Fig.7 Cracks in the wall of building at the middle part

圖9 院內(nèi)出現(xiàn)裂縫Fig.9 Cracks appeared in the yard

圖10 辦公樓墻體變形破壞Fig.10 Deformation and destruction of office building wall

2 滑坡危險(xiǎn)性分析

根據(jù)坡體的結(jié)構(gòu)特征，含碎石土的飽和滲透系數(shù)建議采用6.9×10-6m/s，飽和度為84.3%，采用Van Genuchten模型評(píng)價(jià)不同基質(zhì)吸力下含碎塊石土的滲透系數(shù)(圖11)。由于基巖透水性弱，這里假設(shè)滑床為均一性的不透水材料。在此基礎(chǔ)上，采用Geostudio軟件進(jìn)行10年和50年一遇5日累積降雨的滑坡穩(wěn)定性模擬分析，圖12～圖15為不同工況滑坡穩(wěn)定性狀態(tài)分布情況。

圖11 滑體水土特征曲線(xiàn)(體積含水量隨孔隙水壓力的變化、滲透系數(shù)隨孔隙水壓力的變化)Fig.11 Soil and water characteristic curve of sliding body (volume water content changes with pore water pressure, permeability coefficient changes with pore water pressure)

圖12 工況Ⅰ(自重+水庫(kù)349 m)滑坡穩(wěn)定性Fig.12 Condition I (self-weight & 349 m of reservoir) slope stability

圖13 工況Ⅱ(自重+水庫(kù)353 m)滑坡穩(wěn)定性Fig.13 Condition II (self-weight & reservoir 353 m) slope stability

圖14 工況Ⅲ(自重+水庫(kù)353 m+10年一遇暴雨)滑坡穩(wěn)定性Fig.14 Condition Ⅲ (self-weight & 353 m of the reservoir+a rainstorm in 10 years) slope stability

圖15 工況Ⅳ(自重+水庫(kù)353 m+50年一遇暴雨)滑坡穩(wěn)定性Fig.15 Condition Ⅳ (self-weight & 353 m of reservoir+once rainstorm in 50 years) slope stability

在滑坡滲流場(chǎng)分析的基礎(chǔ)上，利用Slope模塊對(duì)不同工況下滑坡的穩(wěn)定系數(shù)和破壞概率進(jìn)行計(jì)算，結(jié)果見(jiàn)表2。當(dāng)滑坡處于低水位天然狀態(tài)下時(shí)，穩(wěn)定系數(shù)最高為1.703；在庫(kù)水位升高時(shí)，穩(wěn)定性隨之降低為1.540，但也處于穩(wěn)定狀態(tài)，在降雨工況下，穩(wěn)定系數(shù)隨著降雨量的不斷增大而降低?；略?0年一遇暴雨工況下，滑坡的穩(wěn)定系數(shù)降低為1.048，而在50年一遇暴雨工況下，滑坡穩(wěn)定性系數(shù)為0.976，處于不穩(wěn)定狀態(tài)。

表2 不同工況下的穩(wěn)定性系數(shù)Table 2 Stability coefficient under different working conditions

3 災(zāi)害影響范圍內(nèi)承災(zāi)體易損性分析

根據(jù)現(xiàn)場(chǎng)調(diào)查，大水沖水庫(kù)滑坡范圍內(nèi)承災(zāi)體主要考慮為：經(jīng)濟(jì)承災(zāi)體(建筑、大壩和土地)和室內(nèi)人員。主要威脅坡體上西晃山林場(chǎng)建筑及員工，對(duì)滑坡影響范圍內(nèi)的使用建筑物進(jìn)行編號(hào)，按照辦公室及住宅分布計(jì)算，共影響65人，詳細(xì)調(diào)查結(jié)果見(jiàn)表3。將滑坡影響范圍內(nèi)的經(jīng)濟(jì)體進(jìn)行估算，主要涉及辦公樓、室內(nèi)財(cái)產(chǎn)以及滑坡影響范圍內(nèi)的土地價(jià)值，本次調(diào)查以調(diào)查人口及辦公場(chǎng)所為主，房屋建筑面積根據(jù)勘察資料給出。通過(guò)對(duì)房屋居住的調(diào)查顯示，在其中一間房子中每天呆18的小時(shí)，人員的時(shí)空概率為P(S∶T)=0.75?；麦w上的建筑物由于滑坡的整體破壞易于整體倒塌破壞，故易損性取1。室內(nèi)人員的易損性隨建筑物的易損性增加而增加，在滑坡體上由于滑坡的整體下滑，人員的易損性最高，即人員的易損性也取1。

表3 滑坡影響范圍內(nèi)承災(zāi)體調(diào)查表Table 3 Investigation of disaster-bearing bodies within the scope of landslide influence

4 滑坡風(fēng)險(xiǎn)分析

4.1 經(jīng)濟(jì)風(fēng)險(xiǎn)計(jì)算

地質(zhì)災(zāi)害的定量風(fēng)險(xiǎn)評(píng)價(jià)涉及地質(zhì)災(zāi)害的頻率分析及承災(zāi)體的易損性計(jì)算[8]。對(duì)于財(cái)產(chǎn)風(fēng)險(xiǎn)，可用下列計(jì)算公式計(jì)算：

R(prop)=P(H)×P(S∶H)×P(S∶T)×V(props)×E

(1)

式中：R(prop)——財(cái)產(chǎn)風(fēng)險(xiǎn)；

P(H)——滑坡失穩(wěn)概率；

P(S∶H)——滑坡到達(dá)承災(zāi)體的概率；

P(S∶T)——承災(zāi)體時(shí)空分布概率；

V(props)——承災(zāi)體針對(duì)某一滑動(dòng)事件的易損值；

E——承災(zāi)體價(jià)值。

4.2 人口風(fēng)險(xiǎn)計(jì)算

對(duì)于人的生命損失，單人風(fēng)險(xiǎn)可采用下式計(jì)算：

R(LOL)=P(H)×P(S∶H)×P(S∶T)×V(D∶T)

(2)

式中：R(LOL)——個(gè)體的風(fēng)險(xiǎn)值；

P(H)——滑坡失穩(wěn)概率；

P(S∶H)——滑坡到達(dá)承災(zāi)體的概率；

P(S∶T)——承災(zāi)體時(shí)空分布概率；

V(D∶T)——人員的易損性，在一定的災(zāi)害事件下人員的易損程度。

4.3 經(jīng)濟(jì)和人口風(fēng)險(xiǎn)各概率值計(jì)算說(shuō)明

根據(jù)張茂省等[9]，唐亞明等[10]提出的概率值分析方法來(lái)進(jìn)行分析計(jì)算。

滑坡的失穩(wěn)概率P(H)：即滑坡發(fā)生的概率，滑坡發(fā)生概率是一個(gè)比較復(fù)雜的問(wèn)題，一般可以用以下三種情況進(jìn)行表述：①是指在研究區(qū)內(nèi)具備一定特征的滑坡在一年之內(nèi)的發(fā)生數(shù)量；②在給定期限內(nèi)(一般為一年)特定的邊坡發(fā)生滑動(dòng)的概率；③下滑力超過(guò)抗滑力的概率或可能性。本文通過(guò)下滑力超過(guò)抗滑力的概率或可能性的方式，利用Slope/w模塊計(jì)算滑坡破壞概率作為滑坡的失穩(wěn)概率(表2)。

滑坡到達(dá)承災(zāi)體的概率P(S∶H)：其取決于滑坡體與承災(zāi)體的相對(duì)位置及滑坡可能的運(yùn)動(dòng)路徑。它是一個(gè)取值在0～1的條件概率。本文對(duì)位于滑坡體上的建筑物和人員來(lái)說(shuō)取P(S∶H)=1。

承災(zāi)體時(shí)空分布概率P(S∶T)：也是一個(gè)條件概率，這個(gè)概率值也在0～1。建筑物位于在滑坡體上其時(shí)空概率P(S∶T)為1；人員的時(shí)空概率為0.75。

承災(zāi)體的易損性(V)和承災(zāi)體價(jià)值(E)：對(duì)財(cái)產(chǎn)來(lái)說(shuō)，用0(沒(méi)有破壞)到1(完全破壞)來(lái)表示其易損性；對(duì)人員來(lái)說(shuō)，也用0(沒(méi)有傷害)到1(死亡)來(lái)表示其易損性。承災(zāi)體價(jià)值可根據(jù)當(dāng)?shù)氐氖袌?chǎng)價(jià)格予以估算，利用建筑物的概算價(jià)格確定不同樓體的經(jīng)濟(jì)價(jià)值。

4.4 室內(nèi)人員風(fēng)險(xiǎn)

該滑坡影響范圍內(nèi)共有西晃山林場(chǎng)職工65人，根據(jù)以上公式可分別計(jì)算出不同工況下室內(nèi)人口風(fēng)險(xiǎn)值，詳細(xì)計(jì)算情況見(jiàn)表4。

表中給出了四種工況下所有房屋室內(nèi)人員風(fēng)險(xiǎn)值，水庫(kù)水位升高及降雨都導(dǎo)致了室內(nèi)人員風(fēng)險(xiǎn)值的增大，在工況自重+水庫(kù)353 m+50年一遇暴雨工況下室內(nèi)人員的風(fēng)險(xiǎn)值最大，達(dá)到27.8。

表4 不同工況滑坡室內(nèi)人員風(fēng)險(xiǎn)值計(jì)算表Table 4 Calculation of risk value of indoor personnel under different working conditions

4.5 滑坡經(jīng)濟(jì)風(fēng)險(xiǎn)

按照承災(zāi)體經(jīng)濟(jì)價(jià)值調(diào)查結(jié)果，結(jié)合財(cái)產(chǎn)計(jì)算公式可計(jì)算滑坡不同區(qū)域、不同工況下的財(cái)產(chǎn)損失，詳細(xì)計(jì)算結(jié)果見(jiàn)表5。

結(jié)果表明，水庫(kù)水位升高及暴雨都會(huì)導(dǎo)致滑坡經(jīng)濟(jì)風(fēng)險(xiǎn)值增大，水位升至353 m+50年一遇暴雨工況下經(jīng)濟(jì)風(fēng)險(xiǎn)達(dá)到249.3萬(wàn)元，其中土地風(fēng)險(xiǎn)值為0.6萬(wàn)元，房屋及其室內(nèi)財(cái)產(chǎn)風(fēng)險(xiǎn)值為236.9萬(wàn)元，水庫(kù)風(fēng)險(xiǎn)值為11.8萬(wàn)元。

表5 四種工況滑坡經(jīng)濟(jì)風(fēng)險(xiǎn)值計(jì)算表(萬(wàn)元)Table 5 Calculation table of economic risk value of landslide under four working conditions (ten thousand yuan)

為了使得滑坡風(fēng)險(xiǎn)值直觀(guān)展現(xiàn)且便于決策部門(mén)做出最優(yōu)決策方案，將滑坡體本身影響范圍內(nèi)的建筑物及室內(nèi)人員風(fēng)險(xiǎn)值展示見(jiàn)圖16和圖17。

圖16 不同工況下建筑物風(fēng)險(xiǎn)變化圖Fig.16 Risk variation charts of buildings under different working conditions

圖17 不同工況下室內(nèi)人員風(fēng)險(xiǎn)變化圖Fig.17 Risk variation diagram of indoor personnel under different working conditions

4.6 滑坡涌浪風(fēng)險(xiǎn)分析

4.6.1滑坡涌浪計(jì)算方法

(1)美國(guó)土木工程協(xié)會(huì)推薦的滑坡速度計(jì)算方法(修正)：

(3)

根據(jù)對(duì)水阻力的計(jì)算，將美國(guó)土木工程協(xié)會(huì)推薦公式進(jìn)行修正，設(shè)滑坡的初始速度為零，則滑坡入水時(shí)的速度為：

a=g(sinα-cosαtanφ)-CL/m-Cwρfv2S

(4)

(5)

其中：W——滑體單寬重量；

C，φ——滑面巖土體的內(nèi)聚力(kPa)和內(nèi)摩擦角/(°)；

α——滑面平均傾角/(°)；

CW——粘滯阻力系數(shù)，取0.15～0.18；

v——水下運(yùn)動(dòng)速度/(m·s-1)；

S——迎水面積/m2；

H——入水部分滑體重心距離水面的高度/m；

m——滑體質(zhì)量/kg；

L——滑面長(zhǎng)/m。

(2)滑坡入江體積計(jì)算

庫(kù)岸斜坡失穩(wěn)引發(fā)滑坡后，滑坡進(jìn)江條塊運(yùn)動(dòng)過(guò)程中受到水體阻力與條塊速度的平方成正比，使得滑坡各條塊入水時(shí)速度總體趨勢(shì)呈現(xiàn)先增大后減小。在已取得滑坡各條塊運(yùn)動(dòng)速度的基礎(chǔ)上，取滑坡靜止時(shí)入水條塊的總體積為滑坡入江體積。

不同金屬離子對(duì)殼聚糖酶酶活性的影響，Mg2+，Ca2+，Zn2 +，Mn2+，K+，Na+和Cu2+等離子是該酶的激活劑，F(xiàn)e2+，Hg2+，Co2+和Ag+等離子對(duì)該酶有強(qiáng)烈的抑制作用。例如Wang等在研究Serratia marcescens TKU011菌株中得出Mn2+，Cu2 +也可作為該酶的抑制劑。

(3)涌浪計(jì)算

為了探討滑坡初始涌浪高度，潘家錚[11]提出初始浪高的計(jì)算方法，按庫(kù)岸斜坡變形分為水平運(yùn)動(dòng)和垂直運(yùn)動(dòng)兩種，當(dāng)庫(kù)岸斜坡發(fā)生水平運(yùn)動(dòng)時(shí)，引發(fā)的初始浪高可表示為：

(6)

當(dāng)庫(kù)岸斜坡發(fā)生垂直運(yùn)動(dòng)時(shí)，引發(fā)的初始浪高可用下列函數(shù)表示：

(7)

式中：ξ0——引發(fā)的初始涌浪高度/m；

H——水庫(kù)平均水深/m；

v——庫(kù)岸斜坡水平運(yùn)動(dòng)速度/(m·s-1)；

v′——庫(kù)岸斜坡垂直運(yùn)動(dòng)速度/(m·s-1)；

(4)滑坡涌浪在爬坡高度計(jì)算

在計(jì)算涌浪爬坡高度時(shí)考慮斜坡坡角及爬坡方位角兩方面[12-13]。對(duì)于斜坡坡度的影響而言，斜坡坡度越大，涌浪爬坡高度越小，其關(guān)系可表示為

(8)

式中：Δh——爬坡高度(該高度從靜止水面算起，已包括涌浪高度)；

hc——庫(kù)岸邊的涌浪高度；

α——斜坡坡腳(rad)。

對(duì)于爬坡方位角的影響而言，爬坡方位角越大，涌浪爬坡高度越小，其關(guān)系可表示為

(9)

式中：Δh1——爬坡高度(該高度從涌浪表面算起，不包括涌浪高度)；

β——爬坡方位角。

根據(jù)式(8)和式(9)，可以得到涌浪爬坡高度綜合估算式：

(10)

4.6.2滑坡涌浪計(jì)算結(jié)果

根據(jù)穩(wěn)定性計(jì)算結(jié)果，該滑坡僅在自重+地下水+水庫(kù)353 m+50年一遇暴雨工況下發(fā)生滑動(dòng)，為最不利工況，按照美國(guó)土木工程協(xié)會(huì)推薦的速度計(jì)算公式對(duì)滑坡入江速度進(jìn)行了計(jì)算，按運(yùn)動(dòng)方程分條塊的公式計(jì)算了滑坡入江體積，按潘家錚中垂直運(yùn)動(dòng)時(shí)的公式計(jì)算了初始涌浪。計(jì)算參數(shù)見(jiàn)表6，結(jié)果見(jiàn)表7。

表6 涌浪計(jì)算表Table 6 Calculation table of surge

根據(jù)勘察資料，主壩校核洪水位353.75 m，設(shè)計(jì)洪水位350.5 m，水庫(kù)壩體迎水面的坡度為22.5°(rad=0.392 7)，根據(jù)從水庫(kù)壩體與滑坡方向的關(guān)系，取β以最不利計(jì)算，取值為0度，根據(jù)上述公式可計(jì)算得到，在工況Ⅲ、Ⅳ下涌浪水庫(kù)處爬坡高度分別為10.2 m和15.0 m，因此，涌浪會(huì)產(chǎn)生翻壩，涌浪的影響范圍包括壩體下游溝谷及溝口處的建筑。

表7 滑坡入江速度、入江體積與初始涌浪計(jì)算表Table 7 Calculation table of landslide entry velocity, volume and initial surge

4.6.3滑坡涌浪風(fēng)險(xiǎn)初步分析

根據(jù)涌浪危險(xiǎn)性和承災(zāi)體易損性，可計(jì)算滑坡產(chǎn)生涌浪的風(fēng)險(xiǎn)情況[14-15]。涌浪風(fēng)險(xiǎn)計(jì)算公式一般表述為：R=HV(R為風(fēng)險(xiǎn)值；H為危險(xiǎn)性值；V為承災(zāi)體易損性值)。

選擇斜坡失穩(wěn)時(shí)產(chǎn)生涌浪最大影響范圍作為承災(zāi)體的調(diào)查依據(jù)，承災(zāi)體對(duì)象主要為溝口居民和壩體。根據(jù)計(jì)算可知滑坡在工況Ⅳ條件下產(chǎn)生的涌浪影響范圍最大，具體影響范圍見(jiàn)圖18，并且涌浪會(huì)產(chǎn)生翻壩，由德?tīng)柗品ㄈ∥ｋU(xiǎn)性值為0.6，會(huì)對(duì)溝口居民產(chǎn)生影響，但根據(jù)現(xiàn)場(chǎng)調(diào)查，溝口居民距離水庫(kù)較遠(yuǎn)，房屋均為磚混結(jié)構(gòu)，溝口居民處溢洪道修建完整且合理，涌浪翻壩產(chǎn)生洪水的能夠及時(shí)且安全排泄(圖19)，取易損性值為0.1，因此，風(fēng)險(xiǎn)值R=0.06，對(duì)溝口居民安全造成威脅小。

圖18 大水沖水庫(kù)滑坡涌浪范圍預(yù)測(cè)圖Fig.18 Forecast chart of landslide surge area in Dashuichong reservoir

圖19 涌浪翻壩產(chǎn)生洪水的排泄示意圖Fig.19 Schematic diagram of flood discharge caused by dam overturn induced by surge

大水沖水庫(kù)主壩壩高34.79 m，壩長(zhǎng)130 m，為小型均質(zhì)重力壩。根據(jù)現(xiàn)場(chǎng)調(diào)查，大壩位于滑坡左側(cè)，沒(méi)有出現(xiàn)變形現(xiàn)象，結(jié)構(gòu)性質(zhì)良好，能滿(mǎn)足最大涌浪對(duì)它的沖擊，取易損性值為0.15，所以，風(fēng)險(xiǎn)值R=0.6×0.15=0.09，斜坡失穩(wěn)涌浪對(duì)水庫(kù)壩體會(huì)產(chǎn)生影響，但造成破壞的風(fēng)險(xiǎn)低。

因此，根據(jù)初步分析，滑坡產(chǎn)生涌浪對(duì)溝口居民和壩體的風(fēng)險(xiǎn)低。

5 結(jié)論

(1)本文采用Van Genuchten模型評(píng)價(jià)不同基質(zhì)吸力下大水沖水庫(kù)滑坡含碎塊石土的滲透系數(shù)，并利用Geostudio軟件進(jìn)行10年和50年一遇5日累積降雨的滑坡穩(wěn)定性模擬分析。

(2)根據(jù)滑坡風(fēng)險(xiǎn)公式計(jì)算并繪制滑坡影響范圍內(nèi)建筑物風(fēng)險(xiǎn)變化圖和室內(nèi)人員風(fēng)險(xiǎn)變化圖?？偨Y(jié)了大水沖水庫(kù)滑坡最危險(xiǎn)工況為高水位+50年一遇暴雨工況雨，在該工況下，大水沖水庫(kù)滑坡總經(jīng)濟(jì)風(fēng)險(xiǎn)為249.3萬(wàn)元，室內(nèi)人員的風(fēng)險(xiǎn)值達(dá)到27.8。

(3)對(duì)于庫(kù)岸斜坡失穩(wěn)引發(fā)涌浪災(zāi)害是滑坡風(fēng)險(xiǎn)的重要研究部分，因而探討了庫(kù)岸斜坡失穩(wěn)引發(fā)涌浪災(zāi)害的分析方法，大水沖水庫(kù)滑坡變形受庫(kù)水位和降雨影響較大，在工況Ⅳ條件下涌浪影響范圍最大，會(huì)產(chǎn)生翻壩的可能，但滑坡產(chǎn)生涌浪對(duì)溝口居民和水庫(kù)壩體的風(fēng)險(xiǎn)低。