肖 杰 宋 杰 牛志飞

(中兵勘察设计研究院有限公司，北京 100053)

1 工程概况

据勘察资料，某工程场地原始地形为东高西低，按设计要求对场地进行整平，整平后场地形成了弓形边坡，边坡总长约1 500 km，坡高范围约5.0 m～28.0 m，整平时按1∶1进行放坡，坡顶无建筑及荷载。由于放坡后长时间未处理，目前边坡局部地段已发生严重滑坡。如图1所示，北侧的边坡为1号滑坡，南侧的边坡为2号滑坡。

2 1号、2号滑坡强度参数反演

2.1 1号滑坡的强度参数反演

1号滑坡滑动前，整平时按1∶1进行了放坡，结合钻孔资料，将1号滑坡分为2层，上部5 m为黏土层，下部为强风化砂岩、泥岩层。利用FLAC 3D建立3D模型，计算采用摩尔库仑模型，选取不同的强度参数组合，计算得到各种组合下的安全系数见表1。1号滑坡3D模型及各强度组合下的计算结果见图2。据此反算出塌滑体在天然状态下，黏土层力学参数c=50 kPa，φ=10°。强风化砂岩、泥岩层的力学参数为：c=15 kPa，φ=15°。

表1 1号滑坡在各强度组合的安全系数

2.2 2号滑坡的强度参数反演

2号滑坡地层全部为强风化砂岩、泥岩层。同1号滑坡，2号滑坡按照相同算法得到各种组合下的安全系数见表2。2号滑坡3D模型及各强度组合下的计算结果见图2。据此反算出塌滑体在天然状态下强风化砂岩、泥岩层的力学参数为：c=13 kPa，φ=13°。

3 边坡(滑后)稳定性数值模拟分析

3.1 模型网格及计算参数

模型以边坡的走向SN向作为Z轴，EW向作X轴，竖直方向为Y轴。竖直方向坐标采用实际高程坐标，建立直角坐标系，如图3所示。根据地质资料，1号滑坡的模型选取了5个剖面，分别为Z=0 m，75 m，150 m，210 m，250 m。……