基于黏弹塑性本构模型的泥石流数值模拟

泥石流是一个世界范围的地质灾害问题,多年来受到广泛而深入的研究。大多数工作致力于分离地研究岩土材料破坏前的失稳机制和失稳后的快速流动扩展。尝试以一种连续方法模拟完整的泥石流过程,仅用一个本构模型描述泥石流的起始、扩展和停滞。首先讨论了重力增大和基质吸力降低诱发泥石流的力学机制,描述了泥石流过程中岩土材料力学性质的弹塑性–黏性转变现象,对泥石流进行了黏弹塑性的全程模拟,并对挡土墙的冲击效应进行了启发性研究,证明了研究中所使用的本构模型可以准确地描述岩土材料的固流转化现象,数值方法可以很好地对均质连续介质进行大位移模拟。     

Characteristics and remediation of a landslide complex triggered by the 2008 Wenchuan, China earthquake-case from Yingxiu near the earthquake epicenter

During the 2008 Wenchuan earthquake, the river valley from Yingxiu to Wenchuan experienced numerous landslides and became a prominent area of landslide complexes. The present large landslide complex near the earthquake epicenter consisted of Laohuzui slide 1, Laohuzui slide 2 and Douyaping slide. The scale, geology, morphology, sliding process, and failure mechanism of the landslide complex are analyzed by means of field investigation, aerial photograph and stereographic projection technique. Characteristics of these three slides including seismic response of slope, landslide debris, damage and potential failure are discussed: the convex slope and the upslope of fractured granitic rock at high altitude are highly prone to landsliding under earthquake; the high source altitude and long travel path determine grain sizes and the deposit angle of the slide debris; the landslide complex completely buries the G213 roadway and dams up the Minjiang River in these sections; after the earthquake, rainfall, aftershocks and river erosion may retrigger new failures, such as retrogressive slide of weathered fractured rock, colluvial landslide, debris flow, embankment failure and rockfall. The following are presented as suggested remedial measures to protect the roadway and stabilize the slope: the removing and trenching, protective concrete/rock blocks against erosion, retaining structure, rockfall stopping wall, rockfall restraining net, rock bolt, and the planting of vegetation.     

Hydraulic and numerical study on the generation of a subaqueous landslide-induced tsunami along the coast

By carrying out the hydraulic experiments in a one-dimensional open channel and two-dimensional basin, we clarified the process of how a landslide on a uniform slope causes the generation of a tsunami. The effect of the interactive force that occurs between the debris flow layer and the tsunami is significant in the generation of a tsunami. The continuous flow of the debris into the water makes the wave period of the tsunami short. The present experiments apply numerical simulation using the two-layer model with shear stress models on the bottom and interface, and the results are compared. The simulated debris flow shows good agreement with the measured results and ensures the rushing process into the water. We propose that the model use a Manning coefficient of 0.01 for the smooth slope and 0.015 for the rough slope, and a horizontal viscosity of 0.01 m²/s for the landslide; an interactive force of 0.2 for each layer is recommended. The dispersion effect should be included in the numerical model for the propagation from the shore.     

Regional landslide susceptibility following the Mid NIIGATA prefecture earthquake in 2004 with NEWMARK’S sliding block analysis

This study proposes a calculation method for regional earthquake-induced landslide susceptibility that applies the permanent seismic displacement calculated using Newmark’s sliding block analysis with estimated vertical and horizontal seismic motions. The proposed method takes into account the direction of slope failure based on the specified slope azimuth. The study results reveal the importance of predominant slope failure direction using a simple infinite slope model subjected to earthquakes. The target area for the earthquake-induced landslide susceptibility analysis constituted a region of more than 2000 km² surrounding the epicenter of the Mid Niigata prefecture earthquake in 2004. An earthquake-induced landslide susceptibility map was created based on the proposed method with a specific combination of friction angle and cohesion, and the resulting data were compared to the landslide inventory map produced from aerial photographs following the Mid Niigata prefecture earthquake in 2004. To create the susceptibility map, geomaterial cohesion values for the slope were back-calculated to satisfy the minimum safety factor in the static state. This study also proposes a calculation method for the prediction rate and determines the back-calculated strength parameters of geomaterials. The proposed regional landslide susceptibility map will be useful for understanding potential slope failure locations and magnitude of damage, as well as for planning field investigation and preventing secondary disasters immediately after earthquakes.     

宝鸡市黄土工程边坡抗滑挡墙防治失效问题研究

论文基于对陕西省宝鸡市黄土工程边坡防护工程的系统调查测绘,针对普遍存在的抗滑挡墙防治失效问题进行了研究,提出了抗滑挡墙3种失效模式:黄土泥流越顶、墙体剪切破坏和挡墙“坐船随行”。以岭南滑坡抗滑挡墙防治失效为例,开展了黄土土力学测试,运用极限平衡法对滑坡推力进行了复核,分析了不同失效模式对应的成因机制。研究表明:（1）黄土工程边坡失稳多受工程切坡和降雨的双重影响。（2）坡面排水系统不完善、表层土体松散、降雨及坡体局部存在积水洼地是导致黄土泥流发生的主要影响因素;挡墙有效埋深不足是发生“坐船随行”的主因。（3）最后对解决黄土工程边坡抗滑挡墙失效问题进行了对策探讨。     

施工诱发边坡失稳变形机制及应急治理对策

通过对香溪大道失稳边坡地质灾害调查,对边坡的地形地貌、地层岩性、岩体结构、原设计方案与施工过程进行研究,在此基础上通过FLAC3D数值模拟,结合工程地质条件分析,对原格构挡墙逆作法施工诱发边坡牵引式滑动过程进行了分析,并对其变形破坏机制进行深入探讨。研究结果表明:边坡变形受工程地质条件控制作用明显,边坡开挖改变了坡体原来的力学平衡条件,为边坡变形失稳创造了客观条件;由于超开挖施工,受格构挡墙自重作用,下部支撑岩体受压破坏,导致格构挡墙失稳下沉,带动其后坡体变形,形成牵引式滑动;随着坡体滑移的发展,边坡将逐渐拉裂破坏并形成贯通的剪切滑动面,从已施工的格构挡墙下部贯通剪出,最终整体失稳破坏;基于破坏机制分析的应急治理措施将重点放在为原格构挡墙下部提供有效支撑和控制潜在滑动面变形上;监测结果表明,应急治理后边坡达到稳定性要求。      

An analysis of the supply process of loose materials to mountainous rivers and gullies as a result of dry debris avalanches

A dry debris avalanche will produce different volumes of colluviums or depositions (loose materials), which can have a significant impact on mountainous rivers or gullies. The loose material supply process caused by a debris avalanche is an important issue for understanding secondary disasters that form via the coupling of water flow and loose materials. Two flumes were designed for laboratory tests of the loose materials supply process to rivers/gullies, and the related impact factors were analyzed. Experimental results show that the supply of loose materials is a continuous process that directly relates to the avalanche’s mass movement processes. The sliding masses with smaller particle sizes are more sensitive to the flume slope and exhibited a longer supply time. The time-consuming for the debris avalanche travel in the flume decreased with the increasing particle size (such as flume B, time-consuming is decreased 0.2 s when the particle size increased from <1.0 to 20–60 mm), landslide volume and flume slope (flume A, consuming 1.6–2.1 s when flume slope is 29° decreased to consuming 1.3–1.5 s when flume slope is 41°), which means the increasing mobility of loose materials. The total supply time increased with the increasing landslide volume or decreasing particle size and flume slope. An empirical model for the process is presented based on numerous laboratory tests and numerical simulations, which can successfully describe the supply process for loose materials to a river/gully. The supply process of loose materials to mountainous gully from a dry debris avalanche is controlled by the material compositions of sliding masses, topographical conditions, landslide volume and bed friction, where large-volume debris avalanches that occur in mountainous river regions are more likely to obstruct the river flow and form a landslide-dammed lake.     

Spatial risk analysis of Li-shan landslide in Taiwan

By coupling limit equilibrium analysis and Monte Carlo analysis with a geography information system (GIS), this study implements a method that can evaluate the risk (corresponding to probability of failure in this study) of landslide with consideration of spatial uncertainties. The GIS can adopt the three-dimensional information including surface topography, underground geomaterial distribution and groundwater level to determine slope profiles for analysis. Then the safety of defined slope can be evaluated by limit equilibrium analysis. In this study, the mechanical properties of geomaterial were considered as random variables instead of single values. The slope and groundwater profiles are also randomly adopted. Through a Monte Carlo sampling process, a distribution of safety factor and probability of failure can be determined. This probabilistic risk analysis approach was applied to Li-shan landslide in Central Taiwan.

Due to heavy rains, the sites near the highway 7A (mileage 73 k + 150) and the highway 8 (mileage 82 k) in the Li-shan Township began to subside in mid April 1990. Topography, geology, and groundwater condition of this area were first reviewed. Based on this review, together with field investigations and a series of limit equilibrium back analyses, a general hypothetic model was established to illustrate the failure mechanism of this landslide area. Then the developed probabilistic risk analysis model is applied to spatially evaluate the risk of this landslide area as well as the performance of the remediation treatment.     

低山丘陵区典型滑坡-泥石流链生灾害特征与成灾机理

2019年6月10—13日,龙川县发生持续强降雨,导致全县境内发生大量滑坡、泥石流灾害,贝岭镇米贝村是三个重灾区之一。本文以贝岭镇米贝村6号沟发生的滑坡-泥石流链生灾害为研究对象,在野外精细化调查测量基础上,结合数值模拟分析与计算,对链生灾害特征与成灾机理展开研究。研究发现:①6号沟内共发育7处浅层土质小型滑坡,仅3号滑坡体与部分6号滑坡体转化为泥石流,构成泥石流主要物源,其余滑坡未构成持续性影响;②持续降雨下渗,坡体由非饱和向饱和状态转变,坡表形成连续饱和区,孔隙水压力的增加与孔隙水的软化促使土体强度降低,加之坡体饱和自重的增大,斜坡发生浅表层失稳破坏;③降雨的持续下渗与支沟沟源“漏斗状”地形下的地表汇水快速增大滑坡松散堆积体内的含水率,促使其物理性质发生变化,在重力势能下呈流态状启动、运动转化为泥石流。降雨结构影响滑坡-泥石流链生过程,前期降雨引发滑坡、后期降雨启动形成泥石流,滑坡与泥石流的发生表现出阶段性特征。研究成果有助于指导当地政府进一步开展滑坡-泥石流链生灾害的防灾减灾工作,也为该地区未来区域预警研究工作提供理论支撑。     

某挡土墙边坡滑坡形成机理及治理措施

本文对持续暴雨引起的某一复杂挡土墙边坡产生滑坡的形成机理进行了科学分析和预测,计算了边坡稳定性。结合地形条件、规划道路以及经济性,提出在滑坡坡脚修建锚拉桩板墙的治理方案。本研究可为类似滑坡的治理提供参考。     

颗粒流离散元在滑坡运动过程模拟中的应用

介绍了颗粒流离散元的分析原理,以某岩质边坡为实例,构建了计算模型。根据实际资料,在计算模型中加了节理,节理的力学参数弱于岩体。在下方的抗滑挡土墙完全失去作用的情况下,用此模型对该滑坡的运动过程进行了计算预测,分析了整体的变形、不同部位的位移和应力等运动特点,上部粉质粘土区系拉裂引起的下滑,下部碎石土区因挡墙损坏而产生沿节理面的滑动,但拉应力变化不大。指出了在该滑坡存在2个滑动面,并针对2个滑动面提供了相应的治理建议。     

Probabilistic analysis of landslide potential of an inclined uniform soil layer of infinite length: application

An extended probabilistic model that is a modification of the Chen et al. (2007,) model for evaluating the failure probability of an inclined soil layer with an infinite length was developed in the present paper, and then applied to evaluate the occurrence probability of landslide-related debris flow in Tungmen gully located in the eastern Taiwan, which occurred a devastating debris flow in 1990. The statistical properties of hydrogeological parameters were collected and summarized, and then used to evaluate the landslide-related debris-flow probabilities at various relative water depths for Tungmen gully by using the probabilistic model. Under the assumption that the soil is saturated, the soil’s cohesion is negligible and the specific gravity of the solid particles of soils is a constant, a simplified probabilistic critical slope equation for the stability of an infinite slope of soils was also developed, and used to estimate the occurrence probability of debris flow. The result shows that probabilistic landslide analysis for an infinite slope could provide a suitable approximation for the risk analysis of debris flow mobilization at a given gully.     

地震与强降雨作用下堆积体滑坡变形破坏机理及防治方案分析

位于西南山地堆积体滑坡常受到地震和强降雨的双重作用,查明此类滑坡变形破坏机理是地质灾害防治和风险防控的基础。文章的研究对象是鲜水河断裂带附近的炉霍县马居滑坡。研究表明,地震作用对位于斜坡地带堆积体滑坡体结构损伤明显,不但使滑坡整体稳定性下降,还促使坡体内裂隙大量发育,利于降雨入渗,进一步恶化滑坡的水文地质条件。强降雨形成的大规模洪水和泥石流下切坡脚沟道,牵引滑坡体整体向下。长历时强降雨入渗影响坡体稳定性,且在降雨结束后较长时间持续影响坡体稳定性。因此,对此类滑坡防治的对策应考虑坡脚防护和抗滑支挡设置。在对防治方案的有效性分析后,表明防护方案在极端条件下仍然能保障安全性,达防治和风险管控的目的。     

Deterministic Landslide Risk Assessment at a Past Landslide Site

The slope failure and landslide hazard will possess the same properties within a range including the same engineering geological conditions. To assess the landslide risk of a mountainous area, the study of landslides previously having occurred is very important to evaluate the landslide risk around the area in which they took place. Based on the study of the mechanism of a previous landslide recorded in Kumamoto, Japan, this study initially proposes mechanical parameters for evaluating the landslide hazard using a 3D slope stability method. For each slope unit in the study area, the critical slip surface, which reveals the minimum safety factor of a slope, can be obtained. The affected streams and range of possible landslide masses are analyzed based on the debris flow simulation. This is initially applied to simulate the past landslide event and the result shows the landslide-deduced debris flow effectively re-displayed. Overlayered with layers of infrastructure in Geographic Information Systems (GIS), this risk map indicates which houses and road sections remain in dangerous areas.     

四川都江堰三溪村710高位山体滑坡研究

2013年7月10日上午10时,四川都江堰市中兴镇三溪村受极端暴雨影响发生高位山体滑坡灾害,滑坡-碎屑堆积体方量超过150104m3,其中1#滑坡-碎屑堆积体长度1.26km,造成三溪村一组重大人员伤亡。笔者在野外实地调查和室内研究分析的基础上,总结了都江堰三溪村滑坡的基本特征,研究了其启动运动机制和滑动速度,主要认识如下:（1）该滑坡为一处高位山体滑坡,后缘白垩系砂砾岩地层高速滑动后剧烈撞击-铲刮-偏转后铲动坡体上的松散堆积层而形成高位山体滑坡-碎屑流灾害。（2）根据滑坡的运动及堆积特征,将1#滑坡划分为砂砾岩滑动区、碰撞铲刮区和碎屑流堆积覆盖区3部分。（3）7月8日8时至10日8时,中兴镇三溪村的持续强降雨天气过程（都江堰市3d的降雨量相当于该地区年降雨总量的44.1%）,直接触发了滑坡的发生。（4）三溪村滑坡的发生受2008年汶川地震、特殊的岩土体性质、地形地貌条件以及极端暴雨事件的综合影响,地震、地形为其发育提供了基础条件,极端暴雨事件为其直接诱发因素。（5）建议加强高位山体滑坡的研究,尤其是远程滑坡-碎屑流的早期识别和预警。     

Effects of soil heterogeneity on susceptibility of shallow landslides

Uncertainties in parameters of landslide susceptibility models often hinder them from providing accurate spatial and temporal predictions of landslide occurrences. Substantial contribution to the uncertainties in landslide assessment originates from spatially variable geotechnical and hydrological parameters. These input parameters may often vary significantly through space, even within the same geological deposit, and there is a need to quantify the effects of the uncertainties in these parameters. This study addresses this issue with a new three-dimensional probabilistic landslide susceptibility model. The spatial variability of the model parameters is modeled with the random field approach and coupled with the Monte Carlo method to propagate uncertainties from the model parameters to landslide predictions (i.e., factor of safety). The resulting uncertainties in landslide predictions allow the effects of spatial variability in the input parameters to be quantified. The performance of the proposed model in capturing the effect of spatial variability and predicting landslide occurrence has been compared with a conventional physical-based landslide susceptibility model that does not account for three-dimensional effects on slope stability. The results indicate that the proposed model has better performance in landslide prediction with higher accuracy and precision than the conventional model. The novelty of this study is illustrating the effects of the soil heterogeneity on the susceptibility of shallow landslides, which was made possible by the development of a three-dimensional slope stability model that was coupled with random field model and the Monte Carlo method.

     

滑坡碎屑堆积体形成泥石流的实验

地震或强降雨诱发滑坡,滑坡体碰撞解体形成碎屑物质堆积在沟道内,在后期降雨作用下形成泥石流,这是泥石流形成的一种方式,可称为滑坡碎屑堆积体泥石流。笔者分析了影响碎屑堆积体泥石流起动的特征参数,通过实验研究了碎屑堆积体泥石流形成的过程,分析了堆积体表面坡度、黏粒质量分数、中值粒径（d₅₀）以及不均匀系数（Cu）对泥石流形成的影响。结果表明：碎屑堆积体表面坡度对形成泥石流所需单宽流量无明显影响;黏粒质量分数在不大于5%时仅影响碎屑堆积体侵蚀,对碎屑堆积体揭底所需单宽流量无明显影响;影响碎屑堆积体形成泥石流所需单宽流量的因素主要为中值粒径和不均匀系数--随堆积体中值粒径、不均匀系数的增大而增大。通过实验数据拟合得出了中值粒径和不均匀系数与泥石流形成和揭底所需单宽流量的公式;由于公式中只考虑了影响碎屑堆积体形成泥石流的两个主要因素（d₅₀和Cu）,因此对比水文计算结果偏小,但整体趋势基本一致。在实际应用到这类泥石流沟时,可通过修正进行预测。     

山东邹城雨山村西南滑坡灾害特征与防治方案

山东邹城雨山村西南滑坡在20世纪50年代发生过滑动,导致部分农田被毁。该滑坡为浅层土质滑坡,滑体成分主要为含碎石粉质粘土,滑带位于含碎石粉质粘土层与基岩接触带,由于人工切坡,滑带在滑坡前缘局部已揭露,滑床为下伏强风化闪长岩。近年来,频繁的人类工程活动导致该滑坡存在重新滑动的可能,为消除该地质灾害隐患,采取了悬臂式挡墙支挡和截排水结合的防治方案,达到了治理效果。     

大连市中山区小龙街滑坡治理工程

大连小龙街边坡原始地貌陡峭，高差45.28m，由于上覆残积土及后期在建筑施工中堆积了大量的杂填土，在地表水与地下水浸泡作用下，土体失稳，产生类似泥石流性质的滑坡，毁坏原山坡走梯及下部挡墙，造成人员、住宅及财产损失。本文在大连小龙街滑坡治理工程设计时，重点考虑了滑坡的发生机理、治理工程场地施工条件、建筑美学及原有下水管线、走梯等修复工作，采用土石方清理，滑坡前缘一、二级锚拉钢筋混凝土板墙，滑坡中部花坛及走梯，滑坡后缘一、二级钢筋混凝土柱板墙，排水工程和绿化工程六部分进行施工，效果较好。可供相似城区地质灾害治理工程参考。     

九寨沟震区“6·21”泥石流成因与致灾机制研究

2019年6月20日晚九寨沟普降大到暴雨,导致九寨沟景区内卓追沟、下季节海子沟、则查洼沟等多处发生泥石流。3条泥石流沟分别冲出固体物质2.79×10⁴ m³、3.7×10⁴ m³, 2.67×10⁴ m³,造成多处停淤挡墙及拦砂坝损坏,并淤埋了景区内的部分道路,造成了交通的短时间中断。根据灾后地面调查和遥感解译,初步查明了本次泥石流灾害的成因和致灾过程机制。灾害成因：地震引发的崩塌滑坡产生大量的新增松散固体物源与原有的沟坡堆积物(包括老泥石流堆积物)在强降水的作用下形成泥石流,因此本次泥石流是地震和降雨共同作用的结果。通过泥石流发生前后各流域沟道纵剖面和典型横断面调查分析,发现泥石流致灾过程机制主要包括：(1)泥石流冲刷沟道使老泥石流堆积物参与泥石流运动、同时诱发沟岸失稳崩塌,因而增强泥石流规模和冲击力等动力学参数;(2)沟道卡口和崩滑体滑入沟道后导致泥石流堵塞,在后续泥石流的作用下发生级联溃决,进而增大泥石流峰值流量和冲击力等动力学参数。在两种机制的共同...