基于PFC3D的鱼鳅坡滑坡运动过程分析

张家勇，邹银先，杨大山

（贵州省地质环境监测院，贵州 贵阳 550001）

0 引言

贵州省地处西南腹地，是无平原支撑的省份，地形切割剧烈，地质环境脆弱，地质灾害频发[1−3]。2016年5月23日，贵州省开阳县龙岗镇大石板村鱼鳅坡组突发滑坡灾害，造成下方1 户居民2 人死亡，5 间房屋毁坏，10 余亩耕地损毁，直接经济损失90 万元（图1）。鱼鳅坡滑坡为堆积层失稳，是贵州省内典型突发性灾害类型之一。贵州及我国西部山区常见堆积体，在自然、人为因素影响下，容易诱发滑坡地质灾害，给人民的生命财产安全造成了严重的威胁。黄润秋[4]通过大量的研究表明：查清滑坡变形破坏的地质力学模式是滑坡地质灾害防治的基础所在。

图1 滑坡发生前后影像图Fig.1 Image before and after landslide

当前，针对滑坡机理分析主要是借助地质力学方法开展定性分析，判断破坏模式。近年来，随着计算机技术的发展，数值模拟方法在滑坡领域得到了广泛应用。目前的数值模拟方法中，主要包括确定性分析和非确定性分析两类方法，而确定性分析又可分为连续介质分析、非连续介质分析方法以及近年来李世海等[5]提出的基于连续介质力学的离散元方法(CDEM)。其中连续介质数值分析方法有有限单元法、边界元法、有限差分法等，非连续介质分析方法有块体离散元法、颗粒离散元法、关键块体理论、不连续变形分析(DDA 法) 等。CDEM 方法将有限元与离散元进行耦合，在块体内部进行有限元计算，在块体边界进行离散元计算，通过块体内部及块体边界的断裂，不仅可以模拟材料在连续状态下及非连续状态下的变形、运动特性，更可以实现材料由连续体到非连续体的渐进破坏过程。……