1、中 国 地 质 大 学研 究 生 课 程 论 文课程名称 地质灾害风险分析与预警预报 教师姓名 殷坤龙 教授、吴益平 教授 学生姓名 陈仁全 学 号 120110351 专 业 地质工程 所在院系 地质调查研究院 类 别 硕 士 日 期 2012 年 4 月 2 日 评 语对课程论文的评语:平时成绩： 课程论文成绩：总 成 绩： 评阅人签名：注：1、无评阅人签名成绩无效；2、必须用钢笔或圆珠笔批阅，用铅笔阅卷无效；3、如有平时成绩，必须在上面评分表中标出，并计算入总成绩。三峡库区下土地岭滑坡稳定性分析与风险预测陈仁全（中国地质大学（武汉）地质调查研究院 湖北武汉 430073）摘要：滑坡灾害对
2、人类生命、财产造成重大损失，而滑坡的稳定性分析是确定滑坡防治措施的重要依据；同时，滑坡灾害人口、经济易损性及人口伤亡、经济损失风险预测是单体滑坡灾害预警预报工作的一个重要环节。在对三峡库区下土地岭滑坡的有关资料进行分析的基础上，归纳总结影响滑坡稳定性主要影响因素，选定四种工况，应用极限平衡理论，基于岩土力学参数的 Monte Carlo模型，采用 Morgenstern-Price 法分析该滑坡稳定性系数及概率分布。最后，确定滑坡灾害强度、体积、人口密集程度、人口年龄结构、居民对滑坡灾害的防范意识、政府对滑坡灾害的重视程度等 6 个方面对人口易损性的影响，评价了研究区人口易损性；同时确定了研究
3、区经济易损性指数值。在此基础上，通过Rt = EHV 风险预测公式（ UNESCO, 1972; Fournier DAlbe, 1979） ，分别计算出土地岭滑坡人口风险和经济风险预测值，并作简要分析。关键词：下土地岭；极限平衡；稳定性；易损性；风险分析1 引言滑坡的稳定性分析是确定滑坡防治措施的重要依据，也是滑坡灾害研究的重要方面之一。目前，滑坡稳定性分析评价的常用方法主要有两类，一是以极限平衡理论为基础的条分法，二是以弹塑性理论为基础的数值计算方法。其中极限平衡法为我国工程建设所最为常用，主要包括：瑞典圆弧法、Bishop 法、Janbu 法、不平衡推力传递系数法等 1。这类方法主要考虑
4、沿滑动面上的抗滑力（矩）与滑动力（矩）的比值稳定性系数 FS 来衡量大于 1 时，滑坡体是稳定的，反之，滑坡体不稳定 2。但这种方法最大的缺陷是没有考虑材料参数的变异性；孔隙水压力及外荷载的波动性；计算模型的不确定性等 3。因此，将滑坡体的稳定性系数 FS 看作是各种不确定因素作用下的随机变量，并通过滑坡的破坏概率分析来评价滑坡的稳定性更为合理。在此基础上，进行滑坡灾害风险性分析，最终得出风险预测值供风险区划制图、预警预报等工作参考。本文以三峡库区湖北秭归县下土地岭滑坡为例，采用 Morgenstern-Price 法首先计算滑坡稳定性系数及破坏概率，确定灾害的发生概率，然后分析滑坡影响范围内
5、承灾体的价值和易损性，最后计算出灾害体风险预测值，并作简要分析。2 下土地岭滑坡概况2.1 滑坡地理位置下土地岭滑坡区位于三峡库区湖北省秭归县水田坝乡新址规划区北部，水河的南岸。距原县城归州镇约 7 公里，对外交通为公路，并可通过老县城归州镇从长江水路与外界联系（图 2.1） 。图 2.1 地理位置交通图2.2 地形地貌特征滑坡区属于河谷片坡地貌，处在袁水河弯凸岸一侧的凹岸处。地形上为中缓坡过渡到前沿缓坡台地。南西面为山体，其图幅范围内最大高程223m 。北东面为袁水河，为当地最低高程，海拔156m，相对高差 67m。全区可分为袁水河、河谷阶地、剥蚀斜坡地貌、滑坡地貌和人工改造地貌等几个地貌单
6、元。2.3 滑坡空间特征下土地岭滑坡是一老滑坡，滑坡前缘宽度 210m，往中后部滑坡变窄，宽度在 80100m，滑坡纵向长度右侧明显大于左侧，平均长度为 170m。滑坡体平面范围约 20000m2，滑坡体的厚度在纵向上也有变化，中上部滑坡体平均厚度约为 14m，沿江大道下部的滑坡体厚度逐渐变薄，平均厚度为 7m。该滑坡的总体积约为 25 万 m3。滑坡体后缘高程为 204m，中学宿舍坎下 198m 高程处，发育滑坡后壁及陡坎，由于拉裂滑脱作用形成 34 道小坎，坎高 11.5m 左右。后壁平面上呈大曲率弧形。连接后坎的是经人工稍作改造的台地，台地较平坦开阔，高程约 190m，为居民宅地（图 2
7、.2） 。台地到沿江公路之间为中缓坡，坡度 30左右。沿江公路以下总体缓坡地形，经改造成阶梯状耕地。前沿为舌状地貌，已深入到河谷中，中部前凸，两侧发育浅切冲沟，前缘高程 156m。滑坡区地貌形态反映出滑坡的明显特点，与两侧天然地貌极不协调。滑动面总体上呈半圆型弧状，受基岩构造形态控制，两侧向中部倾斜，后部倾角较前部稍陡，总体倾向北东，倾角 10左右，整体滑动方向 NE1020左右。2.4 地质背景 42.4.1 地质构造滑坡区发育一个背斜和一个向斜构造，滑坡区正位于向斜的核部部位。滑坡区未见断层，主要有四组节理裂隙发育，裂隙大多张开，局部充填。图 2.2 下土地岭滑坡平面图2.4.2 地层岩性
8、滑坡出露地层为侏罗系蓬莱镇组中段(J 3p2)砂岩、泥岩及不同成因的第四系堆积物（图2.3） 。其中基岩为侏罗系蓬莱镇组中段(J 3p2)砂岩及粉砂质泥岩互层。岩性为紫色长石砂岩和灰白色长石石英砂岩，厚层状，单层厚度40120cm不等，中细粒结构。紫红色泥岩、泥质粉砂岩，巨厚层状，层理不发育，粉细粒及泥质结构，单层厚120cm，各分层总厚度一般小于5m 。第四系堆积物主要为人工堆填土、残积及残坡积物、冲积物和滑坡堆积物。滑坡堆积物主要成分为粘土、亚粘土混杂块石、碎石。滑坡表层物质以粘土、碎块石为主，往下以大块石加粘土为主，并保留有原基岩的结构特征。 高 程 （ m） 距 离 （ m）18 正
9、常 库 水 位K6172.48TJ315.06J3p2176.4del4Q民 房 沿 江 大 道中 学 宿 舍210591085710651045 工 程 地 质 剖 面 图 （ ） J65K890K3206.地 面 高 程编 号探 井卵 砾 石TJ3编 号孔 口 高 程K3206.8钻 孔杂 填 土 15.0滑 坡 体泥 岩滑 动 面砂 岩粘 土地 下 水 位图 2.3 下土地岭滑坡- 工程地质剖面图2.4.3 水文地质袁水河是本区地表水系干流，为当地地下水排泄最低基准面。滑坡区内发育一条浅切冲沟，位于滑坡东边界，沟内无水，雨季有小股溪流汇水袁水河。地下水类型分为基岩裂隙水及第四系沉积物及崩
10、塌堆积物空隙水。滑坡体中的地下水位变动较大，旱季滑坡体中的地下水埋深410m 。2.4 人类工程活动滑坡体分布有水田改造成的桔树地和居民点，以及沿江大道。滑坡体上人工堆砌的陡砍与平台相间分布，滑坡体中后部平台上为居民居住地。沿江大道从滑坡体中部穿过，为该区主要交通要道。2.5 滑坡物质组成及物理力学参数滑体表层为约 12m 厚的粘土或含砾粘土，其下为 56m 厚的泥岩、砂岩块石夹粘土，粘土含量 1540不等，其余为块石和卵砾石等混杂物，块径从几十厘米到一米不等。滑动面为棕红色粘土组成，可塑，有明显的镜面和擦痕，擦痕平直，深 13mm，宽约 1cm。擦痕走向 345。镜面产状为 3428。滑带土
11、的埋深 6.807.10m。根据土体物理力学性质测试指标，滑动带土体抗剪强度参数包括天然峰值强度和天然残余强度及饱和峰值强度和饱和残余强度。其天然峰值强度标准值 p = 15.80， 3.8 0 ；C p= 11.1kPa， 2.2kPa；天然残余强度标准值 r =14.70, 3.4 0 ；C r= 10kPa， 2.1kPa；其饱和峰值强度标准值 p = 12.90， 3.3 0 ；C p= 9.9kPa， 2.0kPa；饱和残余强度标准值 r =10.50, 3.3 0 ；C r= 9.2kPa， 2.0kPa。3 滑坡稳定性系数计算在过去的几十年间，产生了很多计算稳定性系数的方法。目前
12、边坡稳定性的分析方法归结起来可分为两类：即确定性方法和不确定性方法，确定性方法是边坡稳定性研究的基本方法，它包括极限平衡法（表 3.1） 、数值方法、块体理论法、赤平极射投影法等。它们将影响边坡稳定性的各种因素都作为确定的量来考虑。例如，极限平衡法是通过潜在滑体的受力分析，引入莫尔-库仑强度准则，根据滑坡的力(力矩)平衡，建立边坡安全系数表达式进行定量评价，这种方法由于安全系数的直观性至今仍被工程界广泛应用。数值方法在 20 世纪60 年代被引入边坡稳定性分析中，它包括有限元法、边界元法、离散元法及混合法等。数值方法能从较大范围考虑介质的复杂性、全面地分析边坡的应力应变状态 5。根据地质勘探资
13、料，考虑稳定性计算模型与滑坡特点的适应性，本文采用适合于任意形态滑面的 Morgenstern-price 法来进行稳定性系数计算。表 3.1 极限平衡理论边坡稳定性分析方法基本条件的比较满足平衡条件分析方法力的平衡 力矩平衡 条间力的假定 滑面形状瑞典法 部分满足 部分满足 不考虑土条间作用力 圆弧Bishop 法 部分满足 满足 条间力合力方向水平 圆弧郎畏勒法 部分满足 部分满足 条间力合力方向水平 任意Janbu 法 满足 满足 假定条间力作用于土条底 以上 1/3 处 任意Spencer 法 满足 满足 假定各条间的合力方向相 互平行 任意Morgenstern-price 法 满足
14、 满足 法向和切向条间力存在一 个函数关系 任意Saram 法 满足 满足 对土条侧向力大小分布做 出假定 任意不平衡推力法 满足 不满足 条间的合力方向与前一土 条滑动面倾角一致 任意3.1 计算模型Saskatchewan 大学 Fredlund 于 20 世纪 70 年代提出广义平衡法（GLE） ，GLE 法基于两个安全系数方程，并对条间切向力法向力做了一定假设。其一的方程给出了考虑弯矩 平衡的安全系数法，另一方程则考虑了水平力平衡的mF安全系数法。Spencer （1967 ）提出同时考虑水平力和力矩平衡的想法。在 GLE法中条间剪力是通过方程 Morgensterm Price(19
15、65)给出的方程处理，方程如下：()XEfx式中： 为函数， 为函数的使用率（小数形式表示） ，E 为条间作用的法向)(xf力，X 为条间作用剪力。GLE 法中考虑弯矩平衡的安全系数方程如下： ()tanmcRNuFWxfDd考虑力平衡的安全系数方程如下： (os)()tcos)inf a式中： 有效凝聚力c有效摩擦角孔隙水压力条块底部法向力N条块重量W线荷载 D几何参数,dfxR土体底部倾斜角在两个安全系数的方程中都有一重要变量 ，即各土条底部的法向力，该N力为所有竖向力（vertical）的总和，竖向力的平衡也就得到了满足，底部法向力的定义如下： sinta()cosRLWXFNa当 已除
16、力矩安全系数时， 即为 ，当 已除水平力的安全系数时，mN即为 。FfMorgenstern-Price 法对任意曲线形状的滑裂面进行分析，推导出了既满足力平衡又满足力矩平衡条件的微分方程，是国际公认的最严密的边坡稳定性分析方法，计算简图如 3.1。 ihzi+1il iKW集 中 荷 载滑 动 面 推 力 线abi iPNiTfaPiEXa)滑动面上的力和力臂 b)条块上的作用力图 3.1 Morgenstern-Price 法计算简图3.2 参数概率模型根据统计资料（表 3.1） ，将相应力学参数用于稳定性计算中。又下土地岭滑坡为一老滑坡复活，所以计算参数选取天然残余强度与饱和残余强度进行
17、计算，岩土体天然重度取 21.5kN/m3，饱和重度 23.5 kN/m3。由于岩土性质的差异性，各种参数的随机性，测试和取样的误差，地质体的不均匀性，外界因素变化的随机性等特点，影响滑坡稳定性系数的各参数。因此用 Monte Carlo 模型处理内摩擦角与粘聚力值。Monte Carlo 模型基本原理为：滑坡的稳定性系数 Fs 的一般表达式 5为：（3-1）(,)SFfcp式中： 正应力条件；剪应力条件岩土体内聚力；c岩土体内摩擦角；p外荷载条件。由（3-1 ）确定 FS 也是随机变量，具有期望 Fsc 和方差 Fsv。滑坡的破坏概率 Pf （即 Fs90%危险程度 安全 低危险性 中等危险
18、性 高危险性 不安全4 下土地岭滑坡易损性分析4.1 易损性评价内容及影响因素易损性评价的主要对象是承灾体，目的是确定各类承灾体的易损性参数，为滑坡灾害风险预测提供基础。滑坡灾害易损性由承灾体自身条件和社会经济条件所决定，前者主要包括承灾体类型、数量和分布情况等；后者包括人口分布、城镇布局、厂矿企业分布、交通通讯设施等。主要评价内容包括：划分承灾体类型，调查统计各类承灾体数量及其分布情况，确定承灾体价值，分析各种承灾体遭受不同类型、不同强度滑坡灾害危害时的破坏程度及其价值损失率。经济易损性与灾害强度密切相关，同时也与承灾体自身的抗灾能力有关。经济损失包括直接经济损失和间接经济损失两部分。人口易
19、损性大小受滑坡灾害的强度大小、人口密集程度、人口年龄结构、居民对地质灾害的认识程度及防范风险的意识和观念等多种因素的影响。4.2 承灾体调查数据统计下土地岭滑坡处于西集镇西部，影响范围内有多处居民聚集地，沿江大道贯穿滑坡中下部，滑坡后缘为水田坝初级中学，前缘紧邻袁水河。下土地岭滑坡对象调查如表 4-1，滑坡体上有 14 户共 65 人，滑坡体后缘学生宿舍 3 栋，受影响人数共 1800 人，直接经济损失和间接损失共计 810 万元。下土地岭滑坡受威胁对象调查如表 4.1。表 4.1 承灾体调查表承灾体 承灾体位置 数量 单位 直接损失（万元） 间接损失（万元）居民 滑坡中后部 65 人人口学校
20、学生 滑坡体后缘 1800 人居民住房 滑坡中后部 14 户 120建筑 学生宿舍 滑坡体后缘 9000 m2 450交通 沿江大道 滑坡中部 300 m 404.3 人口易损性分析人口易损性是指人在滑坡灾害中最大可能死伤人数占其灾前人口总数的比重或百分比。影响下土地岭滑坡影响范围内人口易损性的因素主要包括坡灾害的强度、滑坡体积、人口年龄结构、居民对滑坡灾害的防范意识、政府对滑坡灾害的重视程度等。其评价方法可采用公式（4-1） 7进行计算：（4-1）)1()(65432161 CWCWCWVmpi 其中：Vpi下土地岭滑坡影响范围内的易损性指数；W1滑坡灾害的强度易损性影响的权重；W2滑坡体积
21、对易损性影响的权重；W3人口密集程度因素的权重；W4人口年龄结构评价因素的权重；W5受教育程度评价因素的权重；W6政府重视程度评价因素的权重；C1滑坡灾害的强度评价系数；C2滑坡体积评价系数；C3 人口密集程度系数；C4人口年龄结构系数；C5教育程度系数；C6政府重视程度系数。4.3.1 滑坡灾害强度高速滑坡易造成严重的人员伤亡，而蠕变型的滑坡对人员伤亡的易损性远低于高速滑坡。代云霞等 8（2008）将美国土木工程师协会推荐的滑坡速度计算公式修正为公式（4-2） ，以此求下土地岭滑坡的滑动速度。（4-2）)sin(cos2SHmgLtgVfw其中：V入水或剪出口处最大滑速，单位 m/s；滑体重
22、心至水面或剪出口处的垂直落差；H重力加速度，取为 9.8；g滑面倾角；滑动面加权平均摩擦角为滑动时滑面抗剪强度参数；c为滑块与滑面接触面长；L粘滞阻力系数，0.150.18；w浮密度（g/m 3） ；f迎水面积(m 2)。S上式中 可由 CAD 底图中命令 massprop 查找面域物理力学性质功能H（4.3）求得。代云霞等（2008）已通过上式针对 175m 库水位条件下，计算出下土地岭滑坡速度 V=4.34m/s。并根据表 4.3，可以判断该滑坡滑动属于很快速滑动，取其评价系数 C1=0.6。表 4.2 下土地岭滑坡体物理力学参数表面积: 1463.8727周长: 365.1764边界框:
23、 X: 52.4751 - 222.4769Y: 155.1745 - 202.4869质心: X: 122.2752Y: 178.2325惯性矩: X: 46639348.1203Y: 24381969.2939惯性积: XY: 31351126.2397旋转半径: X: 178.4944Y: 129.0573表 4.3 滑动速度评价系数滑速等级 速度(最大速度) 人类反应 评价系数极快 5m/s 无 0.9很快 3m/min 因人而异 0.6快 1.8m/hr 撤退 0.5中等 13m/mon 撤退 0.4慢 1.6m/a 维护 0.3很慢 0.016m/a 维护 0.2极慢 10000评
24、价系数 0.1 0.3 0.5 0.7 0.94.3.3 人口密集程度 评价改滑坡影响范围内人口的密集程度，相同规模的滑坡灾害发生在人烟稀少的地方易损性相对较小；发生在人口密度相对集中的城镇或村庄，滑坡可能造成的生命和财产损失将会增大。该滑坡中后部平台分布居民相对集中，西侧紧邻新集镇，滑坡后缘影响区内分布有水田坝初级中学，因此取该滑坡的影响区内 C3=0.65。4.3.4 人口年龄结构评价该滑坡影响范围内人口的年龄分布情况，一般老人和儿童对滑坡灾害的预防能力较低，若比例较大，则表示这一地区人口易损性较高。该滑坡后缘有学生宿舍，滑坡体上有 14 户人家，则滑坡影响区取 C4=0.2，滑坡分布区取
25、C4=0.25。4.3.5 居民对滑坡灾害的防范意识评价该区居民对地质灾害的认识程度及风险防范的意识和观念等。一般用高、中、低等定性指标来衡量。若居民对地质灾害的认识程度越高，其人口易损性就越低，反之则人口易损性越高。滑坡影响区人口主要为中学生，取教育程度系数 C5=0.8，表示基本都接收了中等教育水平；滑坡分布区上住户取教育程度系数 C5= 0.5。4.3.6 政府的重视程度政府对评价区地质灾害知识的宣传力度、投入减灾防灾工作中的人力和物力等，称为政府重视系数。随着政府重视程度的提高，人口易损性会相应降低，变化范围 01，0 表示政府部门漠视评价区的地质灾害防治工作，1 表示政府非常重视（表
26、 4.5） 7。该滑坡是一个老滑坡，从滑坡平面图可以看出，该滑坡变形比较明显，政府对该滑坡也采取了相应的防治措施，所以评价系数C6=0.3。表 4.5 政府重视系数的确定防治程度 重点防治 次重点防治 一般防治 预防评价系数 0.3 0.5 0.7 0.94.3.7 计算人口易损性值根据以上分析，C 3、C 4 相对于其他因素来说是内因，不难看出，外因对易损性影响较大 9。滑坡体本身的性质对易损性影响也较为明显。滑坡评价系数因滑坡影响和分布范围而不同；同时对于滑坡分布区和影响区，各因素权重有所不同。（1） 滑坡分布区W1=0.25，W 2=0.25，W 3=0.2，W 4=0.1，W 5=0.
27、1，W 6=0.1C1=0.6，C 2=0.3，C 3=0.6， C4=0.25，C 5=0.5，C 6=0.3那么： )1()(61Vpi =0.250.6+0.250.3+0.20.6+0.10.25+0.10.5+0.10.7=0.49（2） 滑坡影响区W1=0.15，W 2=0.15，W 3=0.25，W 4=0.1，W 5=0.15，W 6=0.2C1=0.6，C 2=0.3，C 3=0.6， C4=0.2，C 5=0.8，C 6=0.3那么： )1()(654321 CWCWCVpi =0.150.6+0.150.3+0.250.6+0.10.2+0.150.2+0.20.7=0.
28、4754.3 经济易损性分析经济易损性研究的重要内容是确定承灾体价值损失率，可以用承灾体价值损失率来表示，即指承灾体遭受灾害破坏损失的价值与受灾前承灾体价值的比率 10。在灾后评估中，可通过对承灾体的调查，根据其实际损毁程度，评估核算承灾体的价值损失率。但在以期望损失为目标的风险预测中，准确评价研究区各类承灾体的易损性非常困难。但易损性与滑坡发生强度密切相关，可根据承灾体遭受某种强度的滑坡灾害时可能发生的破坏程度，分析预测承灾体的易损性。滑坡活动越强烈，承灾体的易损性相应越高，并且不同承灾体对不同类型和不同活动强度的滑坡灾害的承受能力不一样，可能的损毁程度及灾后的可恢复性也存在着差异。本研究区
29、主要分布的承灾体类型据现有的资料可以归纳为两类：建筑物和交通设施（表 4.6） 。表 4.6 经济易损性评价表承灾体类型 总量 单位价值总价值（万元）易损性居民住房 14 栋 8.57 万元/栋 120 0.85建筑物学生宿舍 9000m2 500 元/m 2 450 0.15交通设施 沿江大道 300m 1333 元/m 40 0.85滑坡分布区上的承灾体在滑坡发生时将全部遭到破坏，且可恢复程度低， 视为严重损坏，价值损失率取 85%；滑坡后缘等影响区内的承灾体为轻微损坏，价值损失率取 15%9。5 滑坡灾害风险分析进行滑坡灾害风险评估的目的是为实施减灾工程提供决策依据，为布置防治工程和地区
30、规划提供依据。通过风险分析在滑坡灾害中的应用，对单体滑坡进行风险评价，可以对该滑坡进行有效的风险管理，以减少由于滑坡失稳而造成的生命及财产的损失 10。下土地岭滑坡灾害风险分析包括人口风险预测与经济损失风险预测两个部分。采用公式（5-1）进行风险预测计算（UNESCO, 1972; Fournier DAlbe, 1979） 。Rt = EHV (5-1)其中：Rt 总风险，指特定的灾害现象对生命、财产、经济活动、等等，所可能造成的损失；H危险性，特定地区范围内某种潜在的滑坡灾害现象在一定时期内发生的概率；V易损性，灾害现象以一定的强度发生而对受威胁对象所造成的损失程度，易损性可以用 01 的
31、区间来表示，0 表示无损失，1 表示完全损失；E受威胁对象，特定区域内受滑坡灾害威胁的各种对象，包括人口、财产、经济活动、公共设施、土地、资源、环境等等。5.1 人口风险预测根据公式（5-1）确定出未来 1 年内滑坡发生时的伤人口伤亡风险。计算结果见表 5.1。表 5.1 人口伤亡风险预测结果工况 发生概率 易损性 受威胁人数 预测伤亡人数 合计0.49 65 3156m 库水位 8.4%0.475 1800 72 750.49 65 3165m 库水位 10.35%0.475 1800 88 910.49 65 4175m 库水位 12.65%0.475 1800 108 1120.49 6
32、5 18175m 库水位+暴雨 57.29% 0.475 1800 490 5085.2 经济风险预测根据公式（5-1）计算出未来 1 年内滑坡发生时的经济损失风险。计算结果见表 5.2。表 5.2 经济损失风险预测结果工况 发生概率 承灾体 易损性 价值（万元） 风险（万元） 合计（万元）沿江大道 0.85 40 2.9居民住房 0.85 120 8.6156m 库水位 8.4%学生宿舍 0.15 450 5.717.2沿江大道 0.85 40 3.5居民住房 0.85 120 10.6165m 库水位 10.35%学生宿舍 0.15 450 7.021.1沿江大道 0.85 40 4.3居
33、民住房 0.85 120 12.9175m 库水位 12.65%学生宿舍 0.15 450 8.525.7沿江大道 0.85 40 19.5居民住房 0.85 120 58.4175m 库水位+暴雨 57.29%学生宿舍 0.15 450 38.7116.6计算结果表明，受库水位变化及降雨的影响，滑坡稳定性受到严重影响。若滑坡发生破坏，预测伤亡人数 75 人，其经济损失至少 17.2 万元；在强降雨高水位情况下，一旦发生滑坡，预测危险地带伤亡人数将达到 508 人，最大造成约 116.6 万的经济损失。因此，应注意提高当地居民的防灾意识，提高抗灾的应变能力，同时加强监测预报，制定防灾计划，进行
34、防灾演习。一旦滑坡灾害发生，可以立即作出反应，将滑坡灾害造成的损失程度降到最低。特别注意在强降雨条件下，加强“群测群防”工作，及时向上级部门汇报，以便上下做好应急反应工作，将损失减低到最小。6 结论（1）当研究区库水位在 156m 时，稳定性系数较大，滑坡处于低危险，破坏概率为 8.4%， ；当库水位达到 165m 时，稳定系数仍较大，破坏概率稍有增大为 10.35%，处于低危险；当库水位在 175m 时，稳定性系数减小为 1.29，破坏概率为 12.65%，低危险；但在 175m 库水位同时有强降雨出现，滑坡不稳定，稳定性系数1，破坏概率急剧增加为 57.29%，为中等危险。（2）研究区滑坡
35、灾害强度、体积、人口密集程度、人口年龄结构、居民对滑坡灾害的防范意识、政府对滑坡灾害的重视程度等 6 个方面对人口易损性的影响。在滑坡分布区和影响区，其权重值不同。（3）研究区主要分布的承灾体类型为两类类：建筑物和交通设施。滑坡分布区上的承灾体在滑坡发生时将全部遭到破坏，且可恢复程度低，为严重损坏；滑坡后缘等影响区内的承灾体为轻微损坏。（4）风险预测计算结果表明：若滑坡在 156m 库水位时发生破坏，预测伤亡人数 75 人，其经济损失达 17.2 万元；在强降雨和高水位情况下，一旦发生滑坡，预测伤亡人数将达到 508 人，最大造成 116.6 万的经济损失。参考文献1 赵明阶, 何光春 , 王
36、多银. 边坡工程处理技术M. 北京: 人民交通出版社, 20062 周春梅, 殷坤龙. 破坏概率分析及其在浙江张村滑坡中的应用J. 安全与环境工程, 2004(11): 1-43 罗文强, 晏同珍. 蒙特卡洛模拟方法在斜坡稳定性评价中的应用J. 地球科学, 1997(11): 609-6714 吴益平, 殷坤龙, 李崧. 下土地岭滑坡稳定性分析及治理工程设计C. 全国岩土与工程学术大会论文集（上）.北京: 人民交通出版社, 2003: 666-6715 谭文辉，王家臣，刘伟. 边坡稳定性分析方法的探讨J. 露天采煤技术, 1998（2）: 21-236 罗文强, 黄润秋, 张倬元. 斜坡稳定性概率分析的理论与应用M . 武汉: 中国地质大学出版社, 20067 张桂荣, 殷坤龙, 陈丽霞. 浙江省永嘉县区域滑坡灾害人口易损性评价和伤亡风险与预测J. 地质科技情报, 2007(4): 70-758 代云霞,殷坤龙 ,汪洋. 滑坡速度计算及涌浪预测方法探讨J. 岩土力学, 2008, 29（增刊）: 407- 4119 王 琰, 殷坤龙, 陈丽霞. 三峡库区下土地岭滑坡风险分析研究 J. 人民黄河, 2010, 32（9）:114-11710 殷坤龙. 滑坡灾害预测预报M. 武汉: 中国地质大学出版社, 2004: 52-53
