3D FEM numerical simulation of infiltration of landslide considering runoff
-
摘要: 现有数值模拟方法在模拟边坡降雨入渗时通常不能考虑径流流量补给对滑体渗流的影响,当滑床和滑带的渗透性相对于滑体极低时,会出现低估降雨入渗量、高估滑坡稳定性的现象。以Richards方程、运动波方程和有限元法为基础,提出了一种降雨时滑坡渗流数值模拟方法。该方法忽略滑床和滑带的渗流而只考虑滑体渗流过程,从而减小计算规模,也避免了因滑体与滑床(滑带)渗透性差异过大引起的数值计算困难;依据地形和滑体降雨入渗边界饱和情况,修正降雨入渗边界类型和流量,实现了考虑径流流量补给的滑坡降雨入渗简化数值模拟。对简单土质边坡降雨入渗进行数值模拟,从水量平衡与边坡渗流场两个角度,对比了考虑和不考虑径流流量补给的区别;结果表明,考虑径流流量补给时才能保证计算结果水量平衡,正确反映降雨入渗对基质吸力大幅降低的作用。Abstract: When the permeability of landslide bed and landslide zone is very low relatively to the sliding mass, the influences of the runoff from upper slope would not be considered by existing methods in the numerical simulation of rainfall infiltration, so that the feed of precipitation infiltration is underestimated, hence the stability of landslide is overestimated. In this paper, based on Richards's equation, the kinematic wave equation and FEM, a 3D simplified numerical method considering the runoff under the rainfall condition is proposed. In the new method, the calculation domain only contains the sliding mass, to decrease the computational scale and to avoid computational difficulty caused by a significant difference of permeability between the landslide bed (and slide zone) and sliding mass. And what's more, the rainfall boundaries are corrected according to the terrain and net rainfall including not only fluxes but also types (head boundary and flux boundary). The new method is suitable for numerical simulation of the landslide infiltration when the permeability of the landslide bed and landslide zone is very low relatively to the sliding mass, and it can consider the influences of water supply from the runoff to the infiltration of the sliding mass. In order to show the necessity in considering the water supply and the validation of the new method, the infiltration process of a simple soil slope is simulated by the new method and an existing method. The comparisons are made in the aspects of water balance and slope seepage field. The simulated results show that the calculation results, the water balance and the effects of rainfall infiltration on the matric suction can be guaranteed when the runoff supplement is taken into account.
-
Key words:
- slope stability /
- rainfall infiltration /
- runoff supplement /
- FEM /
- numerical simulation
-
表 1 土体SWCC与渗透性函数
Table 1. SWCC and permeability function of soil mass
θ 0.028 0.062 0.085 0.116 0.178 0.265 0.306 0.350 h/cm 200 100 80 60 40 20 10 0 kr 0.001 0.0085 0.010 0.018 0.050 0.886 0.995 1.000 注:θ为体积含水率;h为基质吸力;kr为相对渗透系数。 表 2 雨强和饱和渗透系数
Table 2. Values of R and ks
序号 雨强R/
(10-5cm·s-1)饱和渗透系数ks/
(10-5cm·s-1)R/ks 1 2.894 72.2 0.04 2 5.787 72.2 0.08 3 11.574 72.2 0.16 4 23.148 72.2 0.32 5 2.894 7.22 0.4 6 5.787 7.22 0.8 7 11.574 7.22 1.6 8 23.148 7.22 3.2 -
[1] RICHARDS L A. Capillary conduction of liquids in porous mediums[J]. Physics, 1931, 1: 318-333. doi: 10.1063/1.1745010 [2] TSAI T L, YANG J C. Modeling of rainfall-triggered shallow landslide[J]. Environ Geol, 2006, 50: 525-534. doi: 10.1007/s00254-006-0229-x [3] 娄一青.降雨条件下边坡渗流及稳定有限元分析[J].水利学报, 2007(增刊1): 346-351. http://youxian.cnki.com.cn/yxdetail.aspx?filename=XAGL201705006&dbname=CJFDPREP LOU Yiqing. Finite element analysis of slope seepage and stability due to rainfall infiltration[J]. Journal of Hydraulic Engineering, 2007(Suppl1): 346-351. (in Chinese) http://youxian.cnki.com.cn/yxdetail.aspx?filename=XAGL201705006&dbname=CJFDPREP [4] 李萍, 张琦, 陈小念.降雨条件下饱和-非饱和土坡的渗流分析[J].兰州理工大学学报, 2007, 33(6): 115-118. http://www.cnki.com.cn/Article/CJFDTOTAL-GSGY200706030.htm LI Ping, ZHANG Qi, CHEN Xiaonian. Seepage analysis of saturated-unsaturated slope under rainfall[J]. Journal of Lanzhou University of Technology, 2007, 33(6): 115-118. (in Chinese) http://www.cnki.com.cn/Article/CJFDTOTAL-GSGY200706030.htm [5] 沈银斌, 朱大勇, 蒋泽锋, 等.降雨过程中边坡临界滑动场[J].岩土力学, 2013, 34(增刊1): 60-66. http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2013S1010.htm SHEN Yinbin, ZHU Dayong, JIANG Zefeng, et al. Critical slip field of slope in process of rainfall infiltration[J]. Rock and Soil Mechanics, 2013, 34(Suppl1): 60-66. (in Chinese) http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2013S1010.htm [6] 唐栋, 李典庆, 周创兵, 等.考虑前期降雨过程的边坡稳定性分析[J].岩土力学, 2013, 34(11): 3239-3248. http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201311032.htm TANG Dong, LI Dianqing, ZHOU Chuangbing, et al. Slope stability analysis considering antecedent rainfall process[J]. Rock and Soil Mechanics, 2013, 34(11): 3239-3248. (in Chinese) http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201311032.htm [7] 张培文, 刘德富, 郑宏, 等.降雨条件下坡面径流和入渗耦合的数值模拟[J].岩土力学, 2004, 25(1): 109-113. http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200401025.htm ZHANG Peiwen, LIU Defu, ZHENG Hong, et al. Coupling numerical simulation of slope runoff and infiltration under rainfall conditions[J]. Rock and Soil Mechanics, 2004, 25(1): 109-113. (in Chinese) http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200401025.htm [8] MORITA M, YEN B C, Modelling of conjunctive two-dimensional surface three-dimensional subsurface flows[J]. Journal of Hydraulic Engineering, 2002, 128(2): 184-200. doi: 10.1061/(ASCE)0733-9429(2002)128:2(184) [9] HUSSEIN M, SCHWARTZ F W. Modelling of flow and contaminant transport in coupled stream-aquifer systems[J]. Journal of Contaminant Hydrology, 2003, 65 (1/2): 41-64. [10] PRUDIC D E, KONIKOW L F, BANTA E R. A new stream-flow routing (SFR1) package to simulate stream-aquifer interaction with Modflow-2000[R]. United States Geological Survey, Open File Report 2004-1042. [11] LIANG D, FALCONER R A, LIN B. Coupling surface and subsurface flows in a depth averaged flood wave model[J]. Journal of Hydrology, 2007, 337(1-2): 147-158. doi: 10.1016/j.jhydrol.2007.01.045 [12] SPANOUDAKI K, STAMOU A I, NANOU-GIANNAROU A. Development and verification of a 3-D integrated surface water-groundwater model[J]. Journal of Hydrology, 2009, 375(3-4): 410-427. doi: 10.1016/j.jhydrol.2009.06.041 [13] KOLLET S J, MAXWELL R M. Integrated surface-groundwater flow modeling: a free-surface overland flow boundary condition in a parallel groundwater flow model[J]. Advances in Water Resources, 2006, 29(7): 945-958. doi: 10.1016/j.advwatres.2005.08.006 [14] 童富果, 田斌, 刘德富.改进的斜坡降雨入渗与坡面径流耦合算法研究[J].岩土力学, 2008, 29(4): 1035-1040. http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200804036.htm TONG Fuguo, TIAN Bin, LIU Defu. A coupling analysis of slope runoff and infiltration under rainfall[J]. Rock and Soil Mechanics, 2008, 29(4): 1035-1040. (in Chinese) http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200804036.htm [15] TIAN D, LIU D. A new integrated surface and subsurface flows model and its verification[J]. Applied Mathematical Modelling, 2011, 35(7): 3574-3586. doi: 10.1016/j.apm.2011.01.035 [16] 田东方, 郑宏, 刘德富.考虑径流影响的滑坡降雨入渗二维有限元模拟及应用[J].岩土力学, 2016, 37(4): 1179-1186. http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201604034.htm TIAN Dongfang, ZHENG Hong, LIU Defu. 2D FEM numerical simulation of rainfall infiltration for landslide with considering runoff effect and its application[J]. Rock and Soil Mechanics, 2016, 37(4): 1179-1186. (in Chinese) http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201604034.htm [17] WOOLHISER D A, LIGGET J A. Unsteady, one-dimensional flow over a plane-the rising hydrograph[J]. Water Resources Research, 1967, 3(3): 753-771. doi: 10.1029/WR003i003p00753 [18] 李占斌, 鲁克新.透水坡面降雨径流过程的运动波近似解析解[J].水利学报, 2003(6): 8-13, 21. http://www.cnki.com.cn/Article/CJFDTOTAL-SLXB200306001.htm LI Zhanbin, LU Kexin. Approximate analytical solution of rainfall runoff process on permeable slope[J]. Journal of Hydraulic Engineering, 2003(6): 8-13, 21. (in Chinese) http://www.cnki.com.cn/Article/CJFDTOTAL-SLXB200306001.htm -