Simulation and prediction of the potential gravity erosion during rainfall
LI Tong, WANG Yunqi, QI Zihan, LI Yaoming, HE Xiangchang, LUO Pizhao
1. Three-Gorges Reservoir Area (Chongqing)Forest Ecosystem Research Station, Ministry of Education School of Soil and Water Conservation, Beijing Forestry University, 100083, Beijing, China; 2. Three-Gorges Reservoir Area (Chongqing) Forest Ecosystem Research Station, School of Soil and Water Conservation, Beijing Forestry University, 100083, Beijing, China
Abstract:[Background] Rainfall induced landslide is a common kind of mountain hazard in China, and its mass in near failure state makes it more susceptible to further erosion and may become the main source of debris flow. The traditional slope stability analysis models mainly focus on slope safety assessment, but they are not up to the calculation of the Potential Gravity Erosion (PGE) induced by rainfall triggered landslide. For the absence of such methods, simulation and prediction of PGE is still an urgent work in Soil and Water Conservation. [Methods] An idealized 3D slope model was built in COMSOL to quantify the PGE. The transient pore water pressure and solid deformation development of slope soil under rainfall scenario were simulated in the fluid-solid coupling framework. Furtherly, both the global and local stability of slope were analyzed by using the strength reduction method (SRD) and defining a Local Factor of Safety (LFS) scalar field, and the volume integration of PGE amount quantified based on the LFS. Finally, the practicability of this work was evaluated and confirmed by comparing its slide surface geometries with SRD. [Results] 1) The pore water pressure (PWP) was linearly distributed with the elevation from +100 kPa to -100 kPa in initial state, and decreased significantly with rainfall at surface layer. After 5 d of rainfall, the infiltration depth was about 3.5 m with a transient saturation depth in 2 m. The effective stress increased with depth and parallel distributed to the slope surface in the initial state, and decreased significantly in surface layer and increased slightly in the deep layer during rainfall. The range of average effective stress in 5 d of rainfall changed from 204-16.9 kPa to 206-1.86 kPa. 2) In the early stage of rainfall, the LFS of slope was generally greater than 1, and then gradually decreased on the slope surface with a trend spreading from shallow to deep layer. The shallow straight failure surface (depth<1 m) appeared in 2 d and subsequently evolved into a deep arc-shaped sliding surface. The simulated data of PWP, effective principal stress and LFS were consistent in time and space dimensions, confirming the objective knowledge that rainfall affects slope stability by changing effective stress. 3) In the calculation case, the PGE appeared at 1 d, increasing slowly at first and then rapidly. After 5 d of rainfall, it reached to 230 m3. [Conclusions] The LFS method can achieve the sliding bed profile that is more consistent with SRD with more advantages in computational efficiency and conservatism results. The quantification of PGE can be well realized as illustrated in this work by following steps: 1) field terrain and soil hydrology survey; 2) establishment of seepage-stress coupling model; 3) extraction of LFS equivalent profile; 4) geometric cutting and volume integral.
李通, 王云琦, 祁子寒, 李耀明, 何相昌, 骆丕昭. 降雨型浅层滑坡潜在重力侵蚀量模拟与计算[J]. 中国水土保持科学, 2024, 22(5): 41-52.
LI Tong, WANG Yunqi, QI Zihan, LI Yaoming, HE Xiangchang, LUO Pizhao. Simulation and prediction of the potential gravity erosion during rainfall. SSWC, 2024, 22(5): 41-52.
LU Ning, GONATHON W G. Hillslope hydrology and stability [M]. Cambridge, UK:Cambridge University Press,2013:369.
[2]
NG C W W, SHI Q. A numerical investigation of the stability of unsaturated soil slopes subjected to transient seepage[J]. Computers and Geotechnics, 1998, 22(1): 1.
[3]
XU Qianjun, ZHANG Limin. The mechanism of a railway landslide caused by rainfall[J]. Landslides, 2010,7:149.
[4]
WANG Xinhao, MA Chao, WANG Yunqi, et al. Effect of root architecture on rainfall threshold for slope stability: Variabilities in saturated hydraulic conductivity and strength of root-soil composite[J]. Landslides, 2020,17(8):1965.
[5]
DAI Zhisheng, MA Chao, MIAO Lyu, et al. Initiation conditions of shallow landslides in two man-made forests and back estimation of the possible rainfall threshold[J]. Landslides, 2022,19(5):1031.
[6]
MORGENSTERN N R, PRICE V E. The analysis of the stability of general slip surfaces[J]. Geotechnique, 1965,15(1):79.
[7]
ZHU D Y, LEE C F, QIAN Q H, et al. A concise algorithm for computing the factor of safety using the Morgenstern Price method[J]. Canadian Geotechnical Journal, 2005,42(1):272.
[8]
GRIFFITHS D V, LANE P A. Slope stability analysis by finite elements[J]. Geotechnique, 1999,49(3):387.
[9]
郑颖人, 赵尚毅. 有限元强度折减法在土坡与岩坡中的应用[J]. 岩石力学与工程学报, 2004,23(19):3381. ZHENG Yingren, ZHAO Shangyi. Application of strength reduction FEM in soil and rock slope[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(19):3381.
[10]
孔郁斐, 宋二祥, 杨军, 等. 降雨入渗对非饱和土边坡稳定性的影响[J]. 土木建筑与环境工程, 2013,35(6):16. KONG Yufei, SONG Erxiang, YANG Jun, et al. Rainfall's effect on the stability of unsaturated slopes[J]. Journal of Civil, Architectural & Environmental Engineering, 2013,35(6):16.
[11]
杨涛, 孙立娟, 成启航, 等. 基于位移等值面的边坡点安全系数分析[J]. 铁道学报, 2018,40(6):115. YANG Tao, SUN Lijuan, CHENG Qihang, et al. Analysis of slope point safety factor based on displacement isosurface[J]. Journal of The China Raillway Society, 2018,40(6):115.
[12]
闵江涛, 赵强. 基于点安全系数法的边坡稳定性评价[J]. 杨凌职业技术学院学报, 2020,19(1):10. MIN Jiangtao, ZHAO Qiang. The slope stability evaluation based on point safety factor method[J]. Journal of Yangling Vocational & Technical College, 2020,19(1):10.
[13]
宋林波, 丁增志, 罗智勇, 等. 基坑转角处三维空间效应对稳定性影响分析[J]. 路基工程, 2021,217(4):120. SONG Linbo, DING Zengzhi, LUO Zhiyong, et al. Study on influence of 3D spatial effect at foundation pit corner on stability[J]. Subgrade Engineering, 2021,4:120.
[14]
齐兵, 杨官润, 陈治炜, 等. 边坡点安全系数计算方法比较及应用研究[J]. 路基工程, 2021(6):88. QI Bing, YANG Guanrun, CHEN Zhiwei, et al. Comparison and application research on calculation methods of side slope point safety factoer [J]. Subgrade Engineering, 2021,219(6):88.
[15]
LU Ning, SENER-KAYA B, WAYLLACE A, et al. Analysis of rainfall-induced slope instability using a field of local factor of safety[J]. Water Resources Research, 2012,48(9): W9524.
[16]
SHAO Wei, BOGAARD T A, BAKKER M, et al. Quantification of the influence of preferential flow on slope stability using a numerical modelling approach[J]. Hydrology and Earth System Sciences, 2015,19(5):2197.
[17]
MORADI S, HUISIMAN J A, CLASS H, et al. The effect of bedrock topography on timing and location of landslide initiation using the local factor of safety concept[J]. Water, 2018, 10(10): 1290.
[18]
DOU Zhi, LIU Yimin, ZHANG Xueyi, et al. Influence of layer transition zone on rainfall-induced instability of multilayered slope[J]. Lithosphere, 2021,2021:2277284.
[19]
BISHOP A W. The principle of effective stress[J]. Tecnisk Ukeblad, 1960,39:859.
[20]
RICHARDS L A. Capillary conduction of liquids through porous mediums[J]. Physics, 1931,1(5):318.
[21]
VAN GENUCHTEN M T. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils[J]. Soil Science Society of America Journal, 1980, 44(5): 892.
[22]
TEMGOUAA G T, KOKUTSE N K, KAVAZOVIC Z. Influence of forest stands and root morphologies on hillslope stability[J]. Ecological Engineering, 2016, 95: 622.
[23]
许强, 汤明高, 徐开祥, 等. 滑坡时空演化规律及预警预报研究[J]. 岩石力学与工程学报, 2008,27(6):1104. XU Qiang, TANG Minggao, XU Kaixiang, et al. Research on space-time evolution laws and early warning-prediction of landslides[J], Chinese Journal of Rock Mechanics and Engineering, 2008, 27(6):1104.
[24]
许强, 郑光, 李为乐, 等. 2018年10月和11月金沙江白格两次滑坡-堰塞堵江事件分析研究[J]. 工程地质学报, 2018,26(6):1534. XU Qiang, ZHENG Guang, LI Weile, et al. Study on successive landslides damming events on Jinsha River Baige Village on October 11 and November 3,2018[J]. Journal of Engineering Geology,2018,26(6):1534.
[25]
周礼, 范宣梅, 许强, 等. 金沙江白格滑坡运动过程特征数值模拟与危险性预测研究[J]. 工程地质学报, 2019,27(6):1395. ZHOU Li, FAN Xuanmei, XU Qiang, et al. Numerical simulation and hazard prediction on movement process characteristics of Baige landslide in Jinsha River[J]. Journal of Engineering Geology, 2019,27(6):1395.
[26]
YANG Wentao, WANG Yujie, WANG Yunqi, et al. Retrospective deformation of the Baige landslide using optical remote sensing images[J]. Landslides, 2019,17:659.
[27]
沈冬成. 降雨诱发根土复合坡体失稳特征研究[D]. 成都:成都理工大学, 2020:24. SHEN Dongcheng. Research on instability characteristics of root soil composite slope caused by rainfall[D]. Chengdu:Chengdu University of Technology, 2020:24.
[28]
ZHANG Ke, CAO Ping, LIU Ziyao, et al. Simulation analysis on three-dimensional slope failure under different conditions[J]. Transactions of Nonferrous Metals Society of China, 2011,21(11):2490.
[29]
CHUI T F M, FREYBERG D L. Implementing hydrologic boundary conditions in a multiphysics model[J]. Journal of Hydrologic Engineering, 2009,14(12):1374.
[30]
窦智, 刘一民, 周志芳, 等. 基于单、双重渗透介质降雨边界处理的改进[J]. 岩土力学, 2022,43(3):789. DOU Zhi, LIU Yimin, ZHOU Zhifang, et al. Improvement of rainfall boundary treatment based on single and double permeable media[J]. Rock and Soil Mechanics, 2022,43(3):789.