李亨 郭代泉 谢力欢

摘要：重庆市垫江实验中学男生宿舍西侧，危岩的稳定性直接影响校内师生的生命财产安全，需要对其稳定性进行评估。文章通过对地质资料的分析，采用有限元法对危岩所在边坡的稳定性进行评估。采用稳定性系数法对危岩稳定性进行评估，结果表明：在天然工况下，危岩处于基本稳定状态;在暴雨工况下，危岩处于欠稳定—不稳定状态。采用基于运动学参数的Rockfall数值模拟方法对危岩下落后的运动轨迹进行评估，结果表明：危岩失稳下落后会运动停止在坡脚男生宿舍区域。建议在对边坡及危岩采取相应的防护措施。

[作者简介]李亨（1993—），男，硕士，主要从事岩土工程设计工作;郭代泉（1989—），男，本科，工程师，从事岩土工程施工管理工作。

重庆市垫江县实验中学边男生宿舍西侧坡表面基岩裸露，风化严重，裂隙发育强烈，松动垮塌明显。目前处于稳定—欠稳定状态，在长期内外因素作用下，可能发生崩塌，直接威胁到坡脚在校师生6 000余人及坡顶12户居民的生命财产及房屋安全，危害严重。因此，按DB 50/5029《地质灾害防治工程设计规范》[6]对该边坡的稳定性进行研究是十分有必要的。

本文在对重庆市垫江实验中学边坡及危岩进行工程地质资料分析的基础上，采用有限元法对危岩所在边坡的稳定性和应力分布特征进行研究，采用稳定性系数法对危岩稳定性进行研究，采用基于运动学参数的数值分析方法对危岩下落后的运动特征及运动轨迹进行研究。

1 边坡特征

边坡高度约15.5 m，分布高程425.2～440.3 m，总体坡向250～280°，坡顶为缓坡0～8°，边坡下部为陡崖。陡崖高约10 m，长约70 m，倾角约65～85°。边坡坡面多数裸露，局部崩塌岩块仍堆积在坡脚，陡崖中部生长一棵黄桷树，根系发育（图1）。

2 危岩判别

2.1 模型的建立

结合地形图、地勘资料及现场调查结果，建立边坡计算模型（图2），边坡上部及表层为素填土和强风化砂岩，重度取20 kN/m3，弹性模量取1 000 kPa，泊松比取0.15，粘聚力取5 kPa，内摩擦角取28°，下部基岩为中风化砂岩，重度取24.6 kN/m3，弹性模量取2 967 kPa，泊松比取0.23，粘聚力取1 110 kPa，内摩擦角取36°。

2.2 计算结果及分析

利用有限元数值计算软件对边坡进行网格划分见图3（a），对模型进行边界约束：侧面X向不变形，地面Y向不变形。边坡在天然和地震两种工况下剪应力分布特征见图3（b）、图3（c），场区地震动峰值加速度值0.05g。

应力分析：在自然及地震状态下，边坡坡面中部临空区的岩体产生较明显的剪应力集中现象，且相比天然状态，地震状态下剪应力集中区的岩体剪应力增大20 %。剪应力集中区的岩体处于潜在失稳破状态，存在剪切破坏的可能。结合现场调查结果，该剪应力集中区为危岩发育区域。

3 危岩稳定性分析

3.1 危岩特征

危岩体失稳方式受多方面因素的影响，通过对危岩体形态、发育规模、底界特征、空间分布特征、节理裂隙发育特征、水文地质条件分布规律等进行调查和研究，现状危岩W1、W2（基本特征见表1）尚在倾覆点之内，且随时间的增长，危岩重心将向外位移，危巖的失稳破坏方式主要为滑移式。

3.2 危岩稳定性分析

危岩稳定性系数法的理论计算模型见图4，危岩稳定性系数F的计算见式（1）～式（3）。

式中：U为滑面水压力（kN/m）;hw为裂隙充水高度（m），天然状工况下裂隙充水高度取裂隙深度的1/5，暴雨工况下裂隙充水高度取裂隙深度的2/3;γw为取10 kN/m;Q为地震力（kN/m）;K为危岩稳定性系数;c为后缘裂隙黏聚力标准值（kPa），当裂隙未贯通时，取贯通段和未贯通段粘聚力标准值按长度加权和加权平均值，未贯通段粘聚力标准值取岩石粘聚力标准值的0.4倍;φ为后缘裂隙内摩擦角标准值（°），当裂隙未贯通时，取贯通段和未贯通段内摩擦角标准值按长度加权和加权平均值，未贯通段内摩擦角标准值取岩石内摩擦角标准值的0.95倍;θ为滑面倾角（°）;W为危岩体自重（kN/m）。

危岩W1、W2外倾裂隙少量黏土充填，为软弱结构面，结构面天然黏聚力取30 kPa，天然内摩擦角取17°;饱和黏聚力取25 kPa，饱和内摩擦角取15°;砂岩岩体天然重度取24.6 kN/m3，砂岩岩体饱和重度取24.8 kN/m3。

对于滑移式危岩，当稳定性系数F <1.0时危岩处于不稳定状态;当1.0≤F<1.15时危岩处于不稳定状态;当1.15≤F<1.4时危岩处于不稳定状态;当F≥1.4时危岩处于不稳定状态。危岩W1、W2的稳定性计算结果及稳定性评价见表2。

计算结果表明：在天然工况下，W1处于基本稳定状态，不稳定W2处于基本稳定状态;在暴雨工况下，W1处于欠稳定状态，W2处于不稳定状态。

4 危岩运动特征分析

边坡坡面法相恢复系数取0.3，切向恢复系数取0.6，摩擦角取30°，危岩W1、W2失稳下落后的运动轨迹见图5（a），危岩下落后在边坡上运动时的弹跳高度分布曲线见图5（b）、运动能量分布曲线见图5（c）。

Rockfall模拟结果表明：危岩失稳下落后的运动方式主要是滑动和弹跳，危岩运动停止在坡脚的概率为86%，且集中在距离坡脚1.2～2 m位置处，危岩下落运动过程中的弹跳高度最大值为2.97 m，最大动能为15 054.51 kJ。坡脚为男生宿舍，危岩下落后将会威胁到学校师生的生命安全。

5 结论

（1）边坡坡面临空区的岩体发生剪应力集中，岩体处于潜在失稳破状态，存在剪切破坏的可能，结合现场调查结果，剪应力集中区为危岩发育区域。

（2）在天然工况下，危岩W1、W2处于基本稳定状态;在暴雨工况下，W1处于欠稳定状态，W2处于不稳定状态。在暴雨工况下，随着遇水的下渗和坡面流水的冲刷作用，使得危岩与母岩的摩擦力降低，从而导致危岩的稳定性降低而发生失稳下落。

（3）危岩失稳下落后在边坡的运动方式为滑动和弹跳，且能够运动达到坡脚，运动过程中的弹跳高度最大值为2.97 m，最大动能为15 054.51 kJ。

（4）岩体受裂隙切割，坡面岩体较破碎，建议清除松动危岩体，同时对边坡采用肋柱锚杆和主动防护网的治理措施。

参考文献

[1] 胡厚田. 崩塌与落石[M]. 北京：中国铁道出版社，1989： 78-105.

[2] 何建平， 于远忠. 川中危岩崩塌稳定性分析评价[J]. 四川建筑科学研究.2007，33（5）： 69-72.

[3] Yoichi Okura， Hikaru Kitahara， Toshiaki Sammori， etc. The effeets of rockfall volume on runout distance [J]. Engineering Geology.2000（58）：109-24.

[4] Joachim Schweigl， Carlo Ferretti， Ludwig Nosing. Geotechnical characterization and rockfall simulation of a slope： a practical case study from South Tyrol （Ttaly）[J]. Engineering Geology.2003（67）：281-29.

[5] 楊智翔， 裴向军， 袁进科. 高陡边坡危岩体孤石的稳定性分析[J]. 路基工程.2017，23（9）： 25-29.

[6] 重庆市城乡建设委员会. 地质灾害防治工程设计规范： DB50/5029[S]. 2004.

[7] 唐红梅， 易朋莹. 危岩落石运动路径研究[J]. 重庆建筑大学学报.2003（1）： 17-23.

[8] 蔡冻， 陈东， 吴松山， 蒋楠. 变电站边坡稳定性分析及治理技术[J]. 四川建筑科学研究.2015，41（3）： 70-72.

[9] 亚南， 王兰生， 赵其华， 徐进. 崩塌落石运动学的模拟研究[J]. 地质灾害与环境保护.1996. （2）：25-32.

[10] 韦启珍与雷秀丽， 崩塌落石运动参数的数值模拟[J].中国水运.2008，8（3）： 165-166.

[11] 周文超， 肖财和. 基于ROCKFALL软件计算落石轨迹的分析[J]. 建筑工程技术与设计.2016（3）： 1985-1986.