天然和暴雨条件下的物质点滑坡变形机制分析

2023-11-05艾国栋李云青资丽君梁昌奇

钻探工程 2023年5期

艾国栋，李云青，资丽君，梁昌奇

（1.湖南省矿产资源调查所，湖南郴州 423000；2.中南大学有色金属成矿预测与地质环境监测教育部重点实验室，湖南长沙 410083；3.湖南省有色资源与地质灾害探查湖南省重点实验室，湖南长沙 410083；4.中南大学地球科学与信息物理学院，湖南长沙 410083）

0 引言

边坡失稳产生滑坡是一个典型的岩土工程大变形问题［1］，其造成的破坏后果通常取决于边坡失稳后滑坡的运动速度、滑动距离和堆积形态等［2-5］。滑坡通常具有突发性、滑速快和动能大等特点［6］，严重威胁着其影响范围内的建筑物和人民生命财产安全，如发生在我国的洒勒山滑坡［7］和深圳渣土场滑坡［8-9］，造成了严重的经济损失和不良社会影响。特别是对于体积超过105万方的滑坡，在重力势能的作用下，其滑动距离通常可达其坡高的10 倍以上［10-11］，超过数公里，因而影响范围巨大。因此，为准确定量评估边坡风险，科学合理地分析边坡失稳后的状态，即预测滑体的运动范围和堆积过程，对实际滑坡防灾减灾，具有十分重要的理论和现实意义［12-13］。

传统的边坡稳定性分析方法如极限平衡法和有限元法只能求解边坡破坏前的状态，不能有效处理滑坡大变形问题。物质点法在每个计算时间步都采用规则的背景网格，计算时网格不会扭曲变形，是求解与材料变形有关问题的有效方法，近年来在岩土工程领域中得到了广泛的应用。例如，王双等［14］采用物质点强度折减法研究了边坡失稳后对滑动路径上建筑物的影响程度，并通过算例验证了MPM 预测失稳后滑坡堆积形态和滑动距离的能力；孙玉进等［15］利用物质点法计算大变形的优势模拟了某高边坡从开始变形滑动到稳定的全过程，证明该方法可以较好地再现边坡滑动全过程；Zhao 等［16］将物质点法应用于流体力学模型，模拟了不同纵横比条件下的溃坝流动问题，验证了物质点法的准确性和稳定性；Troncone 等［17］用物质点法分析了1986 年发生在意大利南部的森尼斯滑坡的破坏过程；Yerro 等［18］采用物质点法模拟了华盛顿奥索镇滑坡的滑动过程；Xu 等［19］使用物质点法、DP 本构模型和流变关系，结合ArcGIS 软件对中国红石岩滑坡进行了数值模拟；Conte 等［20］用物质点法对1986 年发生在意大利南部Senise 滑坡的径流过程进行了二维分析；Coelho等［21］和Dong 等［22］使用物质点法研究了刚性基础上的滑坡流动问题，并用物质点法和深度平均法对海底滑坡的实际滑动过程进行了反演分析，再现了从陡坡到中等坡度的滑坡滑动过程；Li 等［23］使用物质点法对王家岩滑坡破坏后的滑动过程进行了数值模拟；Huang 等［24］使采用物质点法和DP 本构模型对粘性土坡和无粘性土坡的大变形破坏过程进行了数值分析。

上述研究证明了物质点法在模拟边坡滑坡大变形方面具有较强可行性，但鲜有研究从工程实际出发，研究不同工况下的边坡大变形机制。为此，本文以湖南省郴州市某滑坡为例，基于物质点法分析该边坡开挖前后在天然工况和暴雨工况下的大变形特征，并与传统边坡稳定性分析方法作对比。

1 物质点法

物质点法是1994 年Sulsky 等［25］提出的一种用于计算大变形问题的数值计算方法。其基本原理是将连续体离散成一系列质点，每个质点代表一块区域并携带了该区域的所有物质信息，如：质量、密度、速度、加速度、应力、应变等，而网格在空间中固定，不携带任何参数，仅用于动量方程的求解和空间导数的计算。该方法具有拉格朗日质点和欧拉网格的双重优势，通过背景网格与物质点间的交替映射完成数值计算，在每个时间步都采用未变形的背景网格，避免了因网格畸变和单元扭曲带来的数值求解困难，可以很好的解决与边坡大变形有关的问题。

物质点法在每个时间步内的计算可分为3 个阶段，即初始化阶段、拉格朗日计算阶段和映射阶段［26］（如图1 所示）。在初始化阶段，物质点携带的所有物质信息通过线性形函数被映射到背景网格上，为拉格朗日计算阶段做准备；在拉格朗日计算阶段，在背景网格上计算动量方程和导数，并通过计算获得新的节点速度、加速度和位移等信息；最后在映射阶段，更新所有物质点状态信息，如物质店的位移、速度等物理量在该时间步计算结束后，背景网格恢复至初始化阶段时的形态，并在下一个时间步重复上述三个阶段的计算。

图1 物质点法计算流程Fig.1 Calculation process of material point method

2 滑坡基本情况介绍

该滑坡位于郴州市第七完全小学北侧，地貌为剥蚀残丘地貌，总体地势北高南低，地势起伏变化较大。场地内及其周边未见有明显活动的断裂构造通过，且新构造运动不明显，处于相对稳定状态，无地表水系，周边未发现崩塌、泥石流、地下采空区等不良地质作用。地层主要为第四系坡残积成因粉质粘土，二叠系下统当冲组全风化泥页岩以及全-强风化炭质页岩，节理裂隙极发育，具体计算参数见表1。

表1 边坡各土层土体参数Table 1 Soil parameters of each soil layer of slope

该滑坡目前处于蠕滑变形阶段，局部已经处于变形加速阶段。地表水下渗后富集在粉质粘土内部滑动面，水量不断富集、运移，该粉质粘土的抗剪强度不断降低，形成软弱滑带，使粉质粘土层内部产生滑动。该滑坡在天然状态下处于基本稳定状态，但在长期降雨或暴雨时，地表水不断入渗滑体，造成滑体自重增加，滑面抗剪强度降低，可能发生更大的位移，严重威胁滑坡体下方居民生命、财产安全。

因此本文用有限元强度折减法［27-28］计算了该边坡开挖前后在天然工况和暴雨工况下的稳定性，并用物质点法［29］计算了开挖前后在天然工况和暴雨工况下的滑动距离，以分析边坡失稳后对滑动路径上建筑物的威胁程度。

3 稳定性分析

利用有限元分析软件MIDAS GTS NX，分析该边坡开挖前后在天然工况和暴雨工况下的稳定性。在构建有限元模型时，根据现场实际情况考虑了不同土层材料。同时，边坡模型左右边界在水平方向进行了固定，底边界在水平方向和垂直方向进行了固定。

其中，局部计算结果通过在有限元计算过程中对边坡模型进行全局搜索滑面得到，而整体计算结果通过在搜索滑面时排除前述安全系数小的局部区域得到。开挖前边坡在天然工况下的计算得到的等效塑性应变云图见图2。由图2 可知，该边坡首先在坡脚处形成较小的等效塑性应变贯通区，安全系数为1.08，随后，在边坡上部产生较大的等效塑性应变贯通区，安全系数为1.61。表明：边坡整体处于稳定状态，但在坡脚处可能会发生局部破坏。

图2 天然工况下开挖前边坡等效塑性应变云图Fig.2 Equivalent plastic strain nephogram of slope before excavation under natural condition

开挖前边坡在暴雨工况下的塑性应变云图见图3。与天然工况下的等效塑性应变云图类似，边坡首先在坡脚处形成等效塑性应变贯通区，安全系数为0.80，随后，在边坡上部形成较大的等效塑性应变贯通区，安全系数为1.22。表明：在暴雨工况下，边坡处于不稳定状态。

图3 暴雨工况下开挖前边坡等效塑性应变云图Fig.3 Equivalent plastic strain nephogram of slope before excavation under rainstorm condition

开挖后边坡在天然工况和暴雨工况下的等效塑性应变云图见图4，安全系数分别为0.60 和0.50，处于不稳定状态。这是因为开挖使得边坡前缘形成临空面，边坡抗滑力减小，下滑力不变，进而导致边坡失稳［30］。

图4 开挖暴雨工况边坡安全系数Fig.4 Slope safety factor after excavation under rainstorm condition

4 大变形分析

由前面的分析可知，开挖前边坡在暴雨工况下处于基本稳定状态，开挖后处于不稳定状态，因此本节使用物质点法［31］计算开挖前后边坡在天然工况和暴雨工况下的滑动距离，以分析其对滑动路径上建筑物的威胁程度。

4.1 计算参数

本研究中，边坡长260 m，高87 m，物质点大小为2.0 m，共离散成6668 个物质点，如图5 所示。为了简化问题，本研究仅考虑平面应变情况。为了更好的还原边坡现场的真实状态，边坡底部选用固定边界条件以模拟基岩，两侧和前后使用对称边坡边界条件以模拟平面应变状态下土体的各种状态。本次计算选用8 节点六面体单元，单元尺寸4 m，采用非关联流动准则和服从Drucker-Prager（DP）屈服准则的弹塑性本构模型。计算时间步长为0.465×10-3s，总计算时间20 s。

图5 边坡离散示意Fig.5 Schematic diagram of slope dispersion

4.2 大变形模拟

开挖前边坡在天然和暴雨工况下的滑动距离计算云图如图6 所示，开挖后边坡在天然和暴雨工况下的滑动距离计算云图如图7 所示。由图6 可知，开挖前边坡在天然工况和暴雨工况下的滑动距离分别为1.37 和3.75 m。由图7 可知，开挖后边坡在天然工况和暴雨工况下的滑动距离分别为20.11 和24.12 m。由此可知，开挖会使边坡形成不稳定临空面，使边坡抗滑力减小，稳定性降低；在雨水的长期作用或暴雨情况下，地表水向边坡土体渗透，滑体自重增加，滑面抗剪强度降低，其稳定性进一步降低，失稳时滑动的距离也进一步增大［32-34］；并且，边坡对于在滑动距离范围内的建筑物威胁程度较大，计算结果可为提高边坡安全性防护提供理论参考。

图6 开挖前边坡滑动距离云图Fig.6 Contour plot of slope sliding distance before excavation

图7 开挖后边坡滑动距离云图Fig.7 Contour plot of slope sliding distance after excavation

由稳定性分析可知，在天然状态下，开挖前边坡的整体稳定性系数为1.61，坡脚的局部稳定性系数为1.08，表明：边坡不会滑动，坡脚处于基本稳定状态。由大变形分析结果可知，天然状态下边坡滑动距离为1.37 m，这是因为在本次计算中视土体在空间上为各项同性，且土体是经过室内试验测试获得的，在测试前土体性质已经扰动，原有性质发生了改变，所以安全系数计算结果存在一定的误差。开挖前暴雨工况和开挖后边坡安全系数均＜1，所以边坡会发生滑动。需要说明的是，由于本文进行的有限元模拟假定土体发生小变形，因此该方法不用于计算滑坡的大变形过程，而用于进行边坡的初始破坏阶段安全系数的计算。在物质点法计算过程中，物质点法可以计算土体的大变形过程，从而用于本文边坡滑动过程的模拟。