考虑残余强度速率依赖性的皖南残积土滑坡运动特征

缪海波; 徐盛; 崔玉龙

doi:10.12404/j.issn.1671-1815.2502950

科学技术与工程 ›› 2026, Vol. 26 ›› Issue (8) :3202-3212 DOI: 10.12404/j.issn.1671-1815.2502950

天文学、地球科学

考虑残余强度速率依赖性的皖南残积土滑坡运动特征

作者信息 +

Movement Characteristics of Residual Soil Landslide in Southern Anhui Considering the Rate-dependency of Residual Strength

Author information +

文章历史 +

摘要

滑坡破坏后的运动过程模拟对于获取滑坡运动学参数及评估其致灾范围和强度具有重要意义。滑动面土体残余强度从根本上控制着滑坡破坏后的运动学特征,然而目前的数值模拟中较少地考虑残余强度在滑坡运动过程中的复杂变异性(如速率依赖性)。以皖南一处典型残积土滑坡——畈章组滑坡为例,对取自滑动面的土样测定其基本物理力学性质指标,然后开展速率步进环剪试验,揭示残余强度对剪切速率的依赖性规律。在此基础上,利用PFC^2D软件开展考虑滑动面土体残余强度速率依赖性的滑坡运动过程模拟。结果表明:在实际滑坡应力水平下,滑动面土体的残余强度具有显著的速率强化特征,这一规律控制了畈章组滑坡破坏后不会表现出快速运动特征;滑坡破坏后的运动时间约78.7 s,在整个运动过程中滑坡表现出明显的运动分区特征,且沿滑动面的运动速度呈现振荡特征,最终堆积状态与现场调查结果基本一致。研究成果可为皖南山区残积土滑坡防灾减灾和风险评价提供有价值的参考。

Abstract

Numerical simulation to the movement of landslides in the post-failure stage is crucial for determining kinematic parameters and assessing the impacted scope and intensity. The residual strength of the sliding surface soil fundamentally governs the kinematic characteristics of landslides in post-failure stage. However, current numerical simulations rarely consider the complex variability of residual strength during landslide movement, for instance, the shear rare-dependency. A shallow landslide on the residual soil slope named Fanzhangzu landslide in southern Anhui was taken as a case study. The soil samples taken from the sliding surface were tested to determine their basic physical and mechanical properties. Rate-stepped ring shear tests were then conducted to reveal the rate-dependency of the residual strength. On this basis, rate-dependent residual strength of sliding surface soil was considered to simulate the movement process of the Fanzhangzu landslide by using PFC^2D. The results show that the residual strength of the sliding surface soil shows significant rate-strengthening at the landslide stress level in field-scale. This law demonstrates that the Fanzhangzu landslide will not exhibit high-speed movement after initiation. the movement duration after initiation of the Fanzhangzu landslide is approximately 78.7 s. Throughout the movement process, the landslide exhibits distinct zoning characteristics in its movement pattern. The velocity along the sliding surface demonstrates oscillatory features, final deposition state is basically consistent with the on-site investigation. The research results can provide valuable references for the disaster prevention, mitigation, and risk assessment of residual soil landslides in mountainous area in southern Anhui.

Graphical abstract

关键词

残积土滑坡 / 运动特征 / 残余强度 / 速率依赖性 / 数值模拟 / 皖南山区

Key words

residual soil landslide / movement characteristics / residual strength / rate-dependency / numerical simulation / mountainous area in southern Anhui

引用本文

引用格式 ▾

缪海波, 徐盛, 崔玉龙. 考虑残余强度速率依赖性的皖南残积土滑坡运动特征[J]. 科学技术与工程, 2026, 26(8): 3202-3212 DOI:10.12404/j.issn.1671-1815.2502950

MIAO Hai-bo, XU Sheng, CUI Yu-long. Movement Characteristics of Residual Soil Landslide in Southern Anhui Considering the Rate-dependency of Residual Strength[J]. Science Technology and Engineering, 2026, 26(8): 3202-3212 DOI:10.12404/j.issn.1671-1815.2502950

登录浏览全文

4963

注册一个新账户忘记密码

近年来,因气候变化、地震和人类工程活动等影响,滑坡发生的强度和规模日益加剧,已成为山区生态系统中最危险的地质灾害类型之一,对人类生命财产、基础设施和生态环境构成了严重威胁^[1-2]。皖南山区地属亚热带湿润季风气候,雨量充沛,地表岩石风化强烈,形成不同厚度的风化残积层。作为安徽省地质灾害最为发育的地区,每年汛期的强降雨可诱发大范围的群发性地质灾害,且以小型残积土浅层滑动为主^[3-4]。滑坡破坏后致灾范围和强度的准确评估是开展此类地质灾害防灾减灾的前提。因此,开展皖南山区残积土滑坡失稳后的运动过程研究,对保障皖南山区旅游业的持续开发、人民群众的生命财产安全以及工农业基础设施的正常运行等均具有重要意义。

滑坡破坏后的运动特征参数如速度、运动距离、堆积状态等是揭示滑坡运动演化特征、预测滑坡运动范围和致灾强度以及评估类似工程地质条件下滑坡风险至关重要^[5]。目前常用的方法有经验统计模型、力学模型、物理模型试验和数值模拟等。经验模型是根据大量已发生滑坡样本的运动距离与几何因素之间的统计规律来评估滑坡的运动特征参数,但不涉及滑坡动力学参数^[6]。力学模型方面,主要是基于滑块模型建立动力学或运动学方程得出滑动速度的解析解或数值解,模型中考虑了滑动面土体的力学参数或流变学参数等,相比经验统计模型更具物理意义^[7-8]。物理模型试验能直观地观测坡体运动速度、运动距离和堆积过程等特征参数^[9],但耗时长、成本高且易受模型尺寸效应的影响;因此,近年来基于离散元(discrete element method,DEM)的数值模拟方法因其应用方便和可重复性等特点受到学者们的青睐。Wang等^[10]将岩体节理形态变化引起的抗剪强度退化引入到不连续变形分析(discontinuous deformation analysis,DDA)中模拟了东河口滑坡的运动过程。Li等^[11]基于弹塑性和非牛顿流变学Drucker-Prager模型,利用三维光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)分析方法,模拟了白格滑坡的运动过程。宰德志等^[12]基于物质点法(material point method,MPM)开展了随机地震作用下边坡滑动距离概率预测。陈建强等^[13] 利用MatDEM软件对四川省汉源县中海村滑坡运动过程进行了数值模拟分析。值得注意的是,类似研究中或将滑动面的摩擦强度视为常数,或考虑抗剪强度的弱化效应,较少地考虑滑动面土体抗剪强度随滑坡运动过程的复杂变异性,如残余强度(或摩擦系数)的速率效应等^[14-15]。

事实上,在滑坡失稳后的运动过程中,因滑动面土体经历大位移剪切,此时残余强度将控制着滑坡的运动学特征^[16]。大量的室内试验结果表明,土体的残余强度往往具有一定的剪切速率效应,表现出速率强化、速率弱化或两者兼具的复杂特征^[17-21]。因此,滑坡破坏后的运动学数值模拟中需引入滑动面土体抗剪强度的速率依赖性才能得出更符合实际情况的模拟结果。

综上分析,以皖南山区一处降雨型残积土滑坡为例,开展现场地质调查并取滑动面土样,首先测定试样的基本物理力学指标和矿物成分,然后开展速率步进环剪试验获得试样残余强度的剪切速率依赖性规律,最后采用PFC(particle flow code)离散元软件开展考虑残余强度速率依耐性的滑坡运动过程数值模拟,获得其运动演化特征参数。研究成果可为皖南山区类似工程地质条件下的滑坡防灾减灾和风险防控提供有价值的参考。

1 滑坡案例

1.1 滑坡工程地质概况

畈章组滑坡(117°40'02.28″E, 30°40'06.24″N)位于安徽省池州市梅龙街道新龙村畈章组X009县道一侧。该滑坡于2016年7月3日发生整体滑动。滑坡区为皖南低山丘陵地貌,自然坡角13°~40°,坡向110°。该滑坡平面上呈长舌形,面积1.24×10⁴ m²,主滑方向100°~108°,体积4.85×10⁴ m³,平均厚度3.9 m(厚度3~8 m)。滑坡冠部之下的壁高为1~4 m,滑坡堆积物堵塞了坡脚附近的县道,但没有到达前方的九华河。该滑坡失稳破坏后的卫星影像如图1所示,沿主滑方向的工程地质剖面如图2^[22-23]所示。

滑体物质由第四系残积土(

${Q}_{4}^{el}$

)组成,含风化岩屑和少量粒径为0.5~20 cm的砾石。下伏基岩为中风化中粒石英二长岩(η

${o}_{5}^{3-1}$

)和上泥盆统五通组(D₃w)砂岩。滑动面为基岩和上覆残积土的分界面,滑带特征不明显,但可见富含灰白色高岭土条带的粉质黏土(如图1所示的现场采样照片)。

根据安徽省专业气象台池州站的气象监测资料,绘制2016年6月1日—7月31日(共2个月)的日降雨量和累积降雨量,如图3所示。可以看出,2016年6月28日、7月2日和7月3日的降雨量分别为173.3、216.4、123.9 mm,均达到中国气象台界定的大暴雨级别。此外,2016年6月28日—7月3日共6 d的累积降雨量达606.5 mm。在经历前两次大暴雨后,畈章组滑坡在第三次大暴雨时启动破坏,表明该滑坡为皖南山区典型的降雨型滑坡。

1.2 试样物理力学性质和矿物成分测定

试样取自滑坡前部堆积体中的岩土分界面处(图1)。考虑到试验用土量较大,除少部分用于测定试样孔隙比和干密度的原状环刀样外,其余均为扰动样。将自然风干的扰动样碾散后,过筛去除粒径大于5 mm的风化岩屑得到重塑样,标记为FZZ,其中细粒组(粒径<0.075 mm)、砂粒组(0.075~2 mm)和砾粒组(2~5 mm)的含量分别为2.47%,63.93%和33.60%。根据《土工试验方法标准》(GB/T 50123—2019)^[24],利用环刀样测定试样FZZ的基本物理性质指标,结果如表1所示。

针对重塑样FZZ,按照1.50 g/cm³的干密度(与取样时的干密度接近,如表1所示)制备三轴试样并在真空饱和器中完全饱和。按照《土工试验方法标准》(GB/T 50123—2019)^[24],在50、100、150和200 kPa的围压(σ₃)下分别开展固结不排水三轴压缩试验。试验过程中的有效应力路径及拟合的破坏主应力线如图4所示。根据破坏主应力线与Mohr-Coulomb强度包络线的关系,可确定试样FZZ的有效峰值黏聚力(c' _p)为21.9 kPa,有效峰值内摩擦角(ϕ' _p)为37.5°。

利用X射线粉晶衍射测得试样FZZ中粒径小于0.075 mm土颗粒(即细颗粒)的矿物成分,其全岩矿物及黏土矿物的衍射图谱如图5所示。可以看出,试样FZZ细颗粒的矿物成分主要为黏土矿物(78.0%)和石英(17.7%),并含少量钾长石、斜长石和黄铁矿,其中黏土矿物主要为高岭石(60.1%)和伊利石(26.7%),以及少量的绿泥石(6.0%)及伊蒙混层(I/S)(7.2%)。

2 环剪试验

2.1 试验设备

环剪试验因其能实现无限位移的剪切而被广泛用于研究土体的残余剪切行为^[25-26]。利用DPRI-5型环剪仪^[27],对取自畈章组滑坡滑动面处的重塑土样FZZ开展速率步进环剪试验并测定其残余强度指标。剪切盒的内、外径分别为12 cm和18 cm,高度为11.5 cm,剪切面积为 141.37 cm²,所能施加的最大法向应力为2 MPa,剪切速率范围为0.01~100 mm/s。该型环剪仪不仅可以测定土体的运动强度,还能模拟剪切带的形成。剪切模式有应力控制和速率控制模式两种。排水条件可由上、下剪切盒上的排水阀控制。为确保高速剪切过程中试样和孔隙流体的损失,在下剪切盒上安装有特制的强抗磨损、低摩阻力和良好密封性的橡皮垫圈。试样的剪切阻力(即抗剪强度)由上剪切盒上的传感器读数减去橡皮垫圈的阻力而得。

2.2 试验程序

将烘干后的试样FZZ采用干落法^[28]装入剪切盒,分层控制装样厚度并轻轻夯实(每层厚度约2 cm),装样高度约7.5 cm,总装样质量为1 563.9 g,干密度约1.48 g/cm³,与取样时的干密度接近(表1)。待试样装填完毕后,首先打开上下排水管阀门通CO₂气体置换试样孔隙中的空气,然后通脱气水使试样饱和。根据Sassa等^[27]提出的饱和度测定方法,测得饱和度系数(B_D)为0.98(>0.95),可认为试样FZZ达到饱和状态。

考虑到畈章组滑坡的坡体厚度(3~8 m)及滑动面倾角(约15°),假定坡体重度为20 kN/m³,从而确定环剪试验的总法向应力为140 kPa。在测定试样饱和度系数(B_D)后将法向应力增加到140 kPa并在排水条件下固结稳定,随后利用速率控制模式开展速率步进环剪试验。剪切速率设为11级,分别为0.01、0.02、0.05、0.10、0.20、0.50、1.0、2.0、5.0、10和20 mm/s。整个试验过程在排水条件下进行,并保持总应力为140 kPa不变。在速率步进环剪试验完成后,在180、220和260 kPa的总应力下,以同一剪切速率(0.01 mm/s)在排水条件下分别剪切至残余状态,获得试样FZZ的残余强度指标。

2.3 试验结果及分析

整个速率步进环剪试验过程中的总应力、抗剪强度、竖向位移变化量和剪切面附近的孔隙水压力随剪切位移的变化如图6所示。需指出的是,竖向位移变化量增加代表试样的高度(或体积)减小,反之则表示试样高度(或体积)的增加。

由图6(a)可知,试样FZZ在0.01 mm/s的初始剪切速率下,抗剪强度随剪切位移先增大至峰值后逐渐降低至一稳定值,表明试样在此剪切速率下达到残余状态。此后,随着剪切速率的增加,试样FZZ的残余强度呈增大趋势,表现出明显的速率强化特征。如图6(b)所示,当剪切速率为0.01~0.10 mm/s时,竖向位移变化量逐渐增加(即剪缩效应),随后在0.20~20 mm/s的范围内有微弱的降低(即弱剪胀效应)。值得注意的是,在0.01 mm/s的初始剪切速率下,试样的剪缩效应最为明显。研究表明^[29-31],环剪试验中土体在剪切过程中的体积减小意味着剪切带的形成,且其形成是粗颗粒破碎与细颗粒运移聚集两者共同作用的结果。因此,在0.01~0.10 mm/s的剪切速率范围内竖向位移变化特征表明试样FZZ的剪切带已经形成。随后在快速剪切过程中(0.20~20 mm/s),可能由于快速剪切扰动了已定向排列的片状或团聚状土颗粒而改变了剪切模式^[32-33],导致了微弱的剪胀效应。

剪切面附近孔隙水压力的变化也从侧面反映了试样的剪胀和剪缩效应。如图6(c)所示,在0.01 mm/s的排水剪切过程中,剪切面附近超孔隙水压力因试样剪缩而快速消散至负值,随后在0.02、0.05和0.10 mm/s三级剪切速率下逐渐上升至0。在快速剪切过程中,因试样有微弱的剪胀效应,剪切面附近的孔隙压力因体积增大而成负压状态,但其绝对值小于3 kPa。如此小的负孔压(弱剪胀效应产生)对试样的残余强度速率强化效应的贡献甚微。因此,在剪切带形成后,快速剪切时残余强度的速率强化可能与剪切面由滑动剪切模式的转变有关。

饱和试样FZZ在0.01 mm/s下的残余强度包络线如图7所示。根据Mohr-Coulomb强度准则,试样FZZ的有效残余黏聚力(c'_r)为28.5 kPa,有效残余内摩擦角(ϕ'_r)为14.5°。与试样FZZ的峰值强度指标相比(c'_p=21.9 kPa,ϕ'_p=37.5°),经历大剪切位移后试样FZZ的黏聚力有较大增长,但内摩擦角则显著降低,降幅约61.3%。这是由于长距离剪切后,土颗粒沿剪切面定向排列,显著降低了摩擦强度,且因体积的压缩,土颗粒连接更为紧密,粒间连结力更强,故而黏聚强度更大^[34]。

环剪试验完成后从剪切盒中取出试样,如图8所示,可以清晰地观察到形成的剪切面以及厚度约23 mm的剪切带,且其土颗粒明显比周围试样要细且颜色较深,表明剪切带的形成是粗颗粒破碎和细颗粒聚集的直接结果。

根据速率与状态摩擦准则(rate-and state-dependent friction law,RSFL)^[35-36],滑带土的稳定摩擦系数(μ_ss,即残余摩擦系数)与滑动(或剪切)速率(v)之间的关系可用对数函数表达,即

(1)

${\mu }_{ss}={\mu }_{0}+(a-b)lg\frac{v}{{v}_{0}}$

式(1)中:v₀为参考速率;μ₀为初始摩擦系数;a-b为RSFL中的摩擦参数,定义为剪切速率效应系数,其数值大小反映土的残余强度剪切速率效应的类型及强弱程度,(a-b)>0为速率强化,反之则为速率弱化。

根据式(1),以0.01 mm/s的剪切速率为参考速率(v₀),将各剪切速率作归一化处理,得到的拟合曲线如图9所示,其表达式为

(2)

${\mu }_{ss}=0.292+0.107lg\frac{v}{{v}_{0}}$

式(2)中:v₀=1.0×10^-5 m/s。

式(2)中剪切速率效应系数(a-b)为0.107,表明畈章组滑坡滑动面土体的残余强度具有显著的运动速率依赖性,且表现出明显的速率强化效应。这一规律将控制着滑坡破坏后的运动特征。具体为:滑坡在破坏后的初始加速度下滑动,因滑动面土体由峰值强度衰减至残余强度,故滑坡加速运动到某一速度,随后因残余强度的速率强化效应导致滑坡逐渐减速,与此同时残余强度也随之减小,此后,滑坡进入到新一轮的加速与减速运动,如此反复直至滑动停止。因此,畈章组滑坡的运动速度随时间表现出振荡特征,且不会表现出快速长距离的运动特征。

3 畈章组滑坡运动过程模拟

作为离散元数值模拟的代表性方法,PFC^2D凭借其丰富的接触模型被广泛运用于岩土体的大变形破坏的精细化模拟,如崩塌滑坡的成灾过程等^[37-38]。利用PFC^2D软件对畈章组滑坡破坏后的运动过程进行数值模拟,获取滑坡的运动速度、运动距离和堆积状况等。根据图2所示的主滑剖面工程地质图,建立畈章组滑坡的PFC^2D数值模型,并在滑动面上自后缘至前缘依次设置12个监测点(P1~P12),如图10所示。滑坡数值模型采用Ball-Wall模型,其中滑坡体为Ball模型,共1 076颗粒,颗粒间为平行黏结;滑床及后缘稳定区为Wall模型,其与颗粒间为线性黏结。

3.1 模型细观参数的标定

模型细观参数的合理性对准确模拟滑坡堆积范围和运动距离有直接影响。中外学者通常利用岩土体室内与数值试验获得的强度参数的一致性来标定PFC数值模型的细观参数^[39-40]。通过比对PFC^2D双轴压缩试验与室内三轴压缩试验结果来标定数值模型的细观参数。在表2所示的细观参数下,得到50、100、150和200 kPa围压下试样破坏的应力-应变关系曲线,如图11所示。

根据图11确定4个围压下的峰值偏差应力值(σ₁-σ₃),结合Mohr-Coulomb强度准则及破坏主应力线得到数值试验的峰值黏聚力(c_p)和内摩擦角(ϕ_p)分别为21.6 kPa和34.5°,与试样FZZ的室内常规三轴压缩试验结果(c'_p=21.9 kPa,ϕ'_p=37.5°,如图4所示)基本吻合,表明表2所列的细观参数可用于畈章组滑坡失稳破坏后的运动过程模拟。但需指出的是,表2中颗粒摩擦系数(μ)在滑坡运动过程的模拟中则采用式(2)来赋值,即考虑滑动面土体残余强度的速率依赖性。

3.2 结果分析

本次模拟得到畈章组滑坡总运动时长约78.7 s。不同时刻滑坡体总位移云图和堆积状态如图12所示。由图12可知,畈章组滑坡自启滑至运动停止的整个过程中表现出明显的运动分区特征(图12中虚线分隔出的运动分区A和B),即中前部滑体(A区)的运动距离明显大于中后部(B区)。在滑坡启滑后30 s内,A区滑体最先开始运动,且A区前缘坡脚处的运动距离较大。此后,A区中后部滑体加速运动,并在整个运动过程中A区中后部和前缘滑体的位移几乎相同,表现出整体运动趋势。B区滑体约在30 s后开始运动,且当运动停止后,B区滑体的运动距离最大约在2 m。从整个滑体的位移云图来看,畈章组滑坡没有超过前缘X009县道,其最终堆积状态与现场调查结果基本一致。

根据图12所示的滑坡位移云图的分区特征,将滑动面上布设的监测点划归为后缘(P1~P3)、中后部(P4~P6)、中前部(P7~P9)和前缘(P10~P12)共4个区域,各区域的平均运动速度随时间变化曲线如图13所示。

由图13可知,后缘和中后部滑体沿滑动面的运动速度基本相同(监测点P1~P6),总体上不超过0.5 m/s,而中前部和前缘滑体的运动速度明显增大(监测点P7~P12),但总体上不超过1.8 m/s。上述各区域滑体的运动速度均表现出振荡式变化特征,表明畈章组滑坡失稳后处于加速运动和减速运动的交替循环中。从运动速度的大小来看,滑坡破坏后的整个运动过程较为缓慢,没有表现出快速运动特征。畈章组滑坡的这一运动特征与滑动面土体的残余强度所表现出的速率强化特征有关(图9)。当滑体加速运动时,滑动面摩擦系数随速度增大而增大,继而增加了抗滑力,随后滑坡减速运动,如此反复直至滑动停止。

图14为滑动面上各监测点的平均运动距离随时间变化曲线。从图14可以看出,滑体中前部和前缘的运动距离最大(A区,监测点P7~P12),且该区域内滑坡体几乎同时运动,最大运动距离约为8 m;滑体中后部(监测点P4~P6)的运动距离约1.6 m,稍小于滑体后缘(监测点P1~P3,约2 m)。此外,前缘及中前部两个区域监测点的平均运动距离随时间的发展几乎同步,表明畈章组滑坡在前缘坡脚处剪切破坏后,中前部滑体随之整体运动。与之明显区别的是,在滑坡启滑后约34 s时刻,中后部及后缘滑体才开始加速运动,这是由于原本支撑该区域的滑体向下运动所致。上述模拟结果显示,考虑滑带土残余强度的速率依赖性能真实地反映畈章组滑坡破坏后的运动特征。

4 讨论

针对畈章组滑坡滑动面土样在实际滑坡应力水平下的速率步进环剪试验清晰地揭示了其残余强度具有明显的剪切速率效应,且表现出单调的速率强化特征。这一试验结果表明,当畈章组滑坡失稳启动后,其在运动过程中滑动面土体的残余强度将会随着运动速率的变化而变化。相似的试验结果也出现在安徽省岳西县一残坡积土滑坡上,即滑坡滑带土的残余强度也同样表现出单调的速率强化特征^[41]。以上关于安徽省典型残积土和残坡积土滑坡的研究表明,在进行此类型滑坡破坏后运动过程反演或预测滑坡运动致灾范围时,忽略滑动面土体残余强度的速率依赖性,而将其视为常量时,往往得到不合理的结论,从而不能准确评估滑坡运动学参数和致灾范围等。因此,在利用离散元软件PFC模拟畈章组滑坡破坏后的运动过程时,将与滑床(Wall模型)接触的颗粒摩擦系数用式(2)给出的函数赋值,而非采用恒定摩擦系数,这与传统的颗粒摩擦系数设定有显著差别^[42],因而能获得符合滑坡岩土体材料实际力学行为的模拟结果,为类似地质环境条件下科学地制定滑坡防灾减灾方案提供理论和数据支撑。

5 结论

以皖南山区一处降雨型残积土滑坡为案例,在现场地质调查的基础上,对取自滑动面的土样测定其基本物理力学性质指标和矿物成分,并开展了滑坡应力水平下的速率步进环剪试验,获得了滑动面土样残余强度的剪切速率依赖性规律,在此基础上,建立主滑剖面的离散元数值模型对其破坏后的运动特征进行模拟。得出如下主要结论。

(1)畈章组滑坡滑动面土体的矿物成分主要为黏土矿物和石英,其中黏土矿物以高岭石和伊利石为主。相对于峰值强度指标,在经历大位移剪切后其残余黏聚力有一定程度的增大,但残余内摩擦角则显著降低,降幅约61.3%。

(2)在实际滑坡应力水平下,畈章组滑坡滑动面土体的残余强度具有显著的速率依赖性,即表现出明显的速率强化特征。这一规律控制了滑坡启动后处于加速和减速的交替运动中,直至运动停止,因而不会表现出快速运动的特征。

(3)PFC^2D模拟得到畈章组滑坡破坏后的运动时间约78.7 s。在整个运动过程中滑坡表现出明显的运动分区特征,且沿滑动面的运动速度呈现振荡特征。滑坡中前部及前缘几乎同时运动,最大运动距离约为8 m;滑坡中后部及后缘的运动时间较中前部和前缘滞后约34 s,最大运动距离1.6~2 m。最终堆积状态与现场调查结果基本一致。

参考文献

原文顺序 | 出版日期 | 本文引用

[1]	Sharma A, Sajjad H, Roshani, et al. A systematic review for assessing the impact of climate change on landslides: research gaps and directions for future research[J]. Spatial Information Research, 2024, 32: 165-185.

[2]	Tian M, Lan H, Li L, et al. Human activities are intensifying the spatial variation of landslides in the Yellow River Basin[J]. Science Bulletin, 2025, 70: 263-272.

[3]	潘国林, 洪天求, 林桂香. 皖南山区滑坡分类及形成机制研究[J]. 合肥工业大学学报(自然科学版), 2016, 39(4): 543-547.

[4]	Pan Guolin, Hong Tianqiu, Lin Guixiang. Research on classifications and formation mechanism of landslides in mountainous area of southern Anhui[J]. Journal of Hefei University of Technology, 2016, 39(4): 543-547.

[5]	Wang C, Lin Q, Wang L, et al. The influences of the spatial extent selection for non-landslide samples on statistical-based landslide susceptibility modelling: a case study of Anhui Province in China[J]. Natural Hazards, 2022, 112: 1967-1988.

[6]	Broeckx J, Rossi M, Lijnen K, et al. Landslide mobilization rates: a global analysis and model[J]. Earth-Science Reviews, 2020, 201: 102972.

[7]	Johnson B C, Campbell C S. Drop height and volume control the mobility of long-runout landslides on the Earth and Mars[J]. Geophysical Research Letters, 2017, 44(24): 12091-12097.

[8]	Zhang H, He S, Liu W, et al. Creep-to-runout transition of large landslides controlled by frictional velocity strengthening and weakening (Vajont 1963, Italy)[J]. Rock Mechanics and Rock Engineering, 2023, 56: 8471-8483.

[9]	杨亮, 汪洋, 张全, 等. 开挖卸荷下顺层岩质滑坡运动距离理论模型[J]. 地球科学, 2024, 49(8): 2851-2861.

[10]	Yang Liang, Wang Yang, Zhang Quan, et al. A theoretical model about the runout distance of bedding rock landslide[J]. Earth Science, 2024, 49(8): 2851-2861.

[11]	贺旭荣, 殷跃平, 赵立明, 等. 基于大型物理模型试验的高位岩质滑坡碎屑流解体破碎效应[J]. 地球科学, 2024, 49(7): 2650-2661.

[12]	He Xurong, Yin Yueping, Zhao Liming, et al. Disintegration and fragmentation effect of high position rock landslide debris flow based on large scale physical model test[J]. Earth Science, 2024, 49(7): 2650-2661.

[13]	Wang J M, Zhang Y B, Chen Y L, et al. Back-analysis of Donghekou landslide using improved DDA considering joint roughness degradation[J]. Landslides, 2021, 18: 1925-1935.

[14]	Li S, Tang H, Peng C, et al. Sensitivity and calibration of three-dimensional SPH formulations in large-scale landslide modeling[J]. Journal of Geophysical Research: Solid Earth, 2023, 128(2): e2022JB024583.

[15]	宰德志, 庞锐, 刘俊. 基于物质点法的随机地震作用下边坡滑动距离概率[J]. 科学技术与工程, 2023, 23(17): 7478-7484.

[16]	Zai Dezhi, Pang Rui, Liu Jun. Study on the probability of sliding distance of slopes under the action of random earthquakes based on thematerial point method[J]. Science Technology and Engineering, 2023, 23(17): 7478-7484.

[17]	陈建强, 沈贺, 冯文凯. 四川省汉源县中海村滑坡特征与运动过程分析[J]. 科学技术与工程, 2024, 24(7): 7448-7454.

[18]	Chen Jianqiang, Shen He, Feng Wenkai. Characteristics and movement process of Zhonghaicun landslide in Hanyuan County, Sichuan Province[J]. Science Technology and Engineering, 2024, 24(17): 7448-7454.

[19]	刘伦杰, 富海鹰, 张迎宾, 等. 考虑滑床摩擦弱化的大光包滑坡运动机制DDA数值模拟研究[J]. 工程地质学报, 2020, 28(6): 1221-1232.

[20]	Liu Lunjie, Fu Haiying, Zhang Yingbin, et al. DDA simulation of mobility of Daguangbao landslide with frictionally weakened sliding bed[J]. Chinese Journal of Engineering Geology, 2020, 28(6): 1221-1232.

[21]	Guo Z, Zhou X, Huang D, et al. Dynamic simulation insights into friction weakening effect on rapid long-runout landslides: a case study of the Yigong landslide in the Tibetan Plateau, China[J]. China Geology, 2024, 7(2): 222-236.

[22]	Kang X, Wang S, Wu W, et al. Residual-state mechanical behavior of soil-rock interface within landslide shear zones[J]. Acta Geotechnica, 2025, 20: 2221-2236.

[23]	Wang G, Suemine A, Schulz W H. Shear-rate-dependent strength control on the dynamics of rainfall-triggered landslides, Tokushima Prefecture, Japan[J]. Earth Surface Processes and Landforms, 2010, 35: 407-416.

[24]	Miao H, Wang G, Yin K, et al. Mechanism of the slow-moving landslides in Jurassic red-strata in the Three Gorges Reservoir, China[J]. Engineering Geology, 2014, 171: 59-69.

[25]	Scaringi G, Hu W, Xu Q, et al. Shear-rate-dependent behavior of clayey bimaterial interfaces at landslide stress levels[J]. Geophysical Research Letters, 2018, 45(2): 766-777.

[26]	Miao H, Wang G. Prediction of landslide velocity and displacement from groundwater level changes considering the shear rate-dependent friction of sliding zone soil[J]. Engineering Geology, 2023, 327: 107361.

[27]	Cueva M, Kang X, Wang S, et al. Unveiling the roleof saturation and displacement rate in the transition from slow movement to catastrophic failure in landslides[J]. Engineering Geology, 2025, 352: 108042.

[28]	姬淑丽. 安徽省池州市新龙村畈章组滑坡稳定性分析与趋势预测[J]. 地质灾害与环境保护, 2016, 27(1): 8-12.

[29]	Ji Shuli. Chizhou in Anhui Province for transshipment neosaurus village group of landslide stability analysis and trend prediction[J]. Journal of Geological Hazards and Environment, 2016, 27(1): 8-12.

[30]	缪海波, 王功辉. 风振影响下乔木坡地暴雨型浅层滑坡演化机制[J]. 地质科技通报, 2022, 41(2): 60-70.

[31]	Miao Haibo, Wang Gonghui. Evolution mechanism of rainstorm-induced shallow landslides on slope covered by arbors considering the influence of wind-induced vibration[J]. Bulletin of Geological Science and Technology, 2022, 41(2): 60-70.

[32]	中华人民共和国住房和城乡建设部, 国家市场监督管理总局土. 工试验方法标准: GB/T 50123—2019[S]. 北京: 中国计划出版社, 2019.

[33]	Ministry of Housing and Urban Rural Development of the People’s Republic of China, State Administration for Market Regulation. Standard for geotechnical testing method: GB/T 50123—2019[S]. Beijing: China Planning Press, 2019.

[34]	Kang X, Wang S, Wu W, et al. Residual state rate effects of shear-zone soil regulating slow-to-fast transition of catastrophic landslides[J]. Engineering Geology, 2022, 304: 106692.

[35]	王顺, 项伟, 崔德山, 等. 不同环剪方式下滑带土残余强度试验研究[J]. 岩土力学, 2012, 33(10): 2967-2972.

[36]	Wang Shun, Xiang Wei, Cui Deshan, et al. Study of residual strength of slide zone soil under different ring-shear tests[J]. Rock and Soil Mechanics, 2012, 33(10): 2967-2972.

[37]	Sassa K, Fukuoka H, Wang G, et al. Undrained dynamic-loading ring-shear apparatusand its application to landslide dynamics[J]. Landslides, 2004, 1(1): 7-19.

[38]	Ishihara K. Liquefaction and flow failure during earthquakes[J]. Géotechnique, 1993, 43(3): 351-415.

[39]	Siman-Tov S, Brodsky E E. Gravity-independent grain size segregation in experimental granular shear flows as a mechanism of layer formation[J]. Geophysical Research Letters, 2018, 45(16): 8136-8144.

[40]	Miao H, Wang G. Shear rate effect on the residual strength of saturated clayey and granular soils under low-to high-rate continuousshearing[J]. Engineering Geology, 2022, 308: 106821.

[41]	Hu W, Zheng Y, McSaveney M, et al. Evolution of the strain localization and shear-zone internal structure in the granular material: Insights from ring-shear experiments[J]. Engineering Geology, 2023, 325: 107283.

[42]	Tika T E, Vaughan P R, Lemos L. Fast shearing of pre-existing shear zones in soil[J]. Géotechnique, 1996, 46(2): 197-233.

[43]	Wang G, Watanabe N, Hoshikawa G, et al. Diverse shear behaviors of clayey materials: implications for differing landsliding behaviors within the same area in Niigata, Japan[J]. Engineering Geology, 2023, 312: 106932.

[44]	Miao H, Wang G. Effects of clay content on the shear behaviors of sliding zone soil originating from muddy interlayers in the Three Gorges Reservoir, China[J]. Engineering Geology, 2021, 294: 106380.

[45]	Dieterich J H. Modeling of rock friction: 1.Experimental results and constitutive equations[J]. Journal of Geophysical Research: Solid Earth, 1979, 84(B5): 2161-2168.

[46]	Ruina A. Slip instability and state variable friction laws[J]. Journal of Geophysical Research: Solid Earth, 1983, 88: 10359-10370.

[47]	孟桓羽, 占洁伟, 卢全中, 等. 陕西山阳“8·12”大型山体滑坡运动特征及数值模拟分析[J]. 工程地质学报, 2023, 31(6): 1910-1928.

[48]	Meng Huanyu, Zhan Jiewei, Lu Quanzhong, et al. Kinematic characteristics and numerical simulation analysis of “8·12” giant landslide in Shanyang County, Shaanxi Province[J]. Journal of Engineering Geology, 2023, 31(6): 1910-1928.

[49]	Nam K, Wang F, Dai Z, et al. Kinetic characteristics and runout behavior of the rainfall-induced Hoengseong landslide at a solar power plant on 9 August 2022 in Gangwon Province, Korea[J]. Landslides, 2023, 10: 1905-1923.

[50]	Li Z, Rao Q H. Quantitative determination of PFC^3D microscopic parameters[J]. Journal of Central South University, 2021, 28(3): 911-925.

[51]	Xu Z H, Wang Z Y, Wang W Y, et al. An integrated parameter calibration method and sensitivity analysis of microparameters on mechanical behavior of transversely isotropic rocks[J]. Computers and Geotechnics, 2022, 142: 104573.

[52]	Xiao Y H, Wei L, Liu X H. Failure mechanism and movement process inversion of rainfall-induced landslide in Yuexi Country[J]. Sustainability, 2025, 17: 5639.

[53]	惠航, 刘锋, 冯兵, 等. 基于PFC^2D的柿子树坪滑坡变形特征模拟研究[J]. 科学技术与工程, 2023, 23(4): 973-981.

[54]	Hui Hang, Liu Feng, Feng Bing, et al. Deformation destroys simulation of Shizishuping landslide based on PFC^2D[J]. Science Technology and Engineering, 2023, 23(4): 973-981.