间歇性强降雨下基于Green-Ampt入渗模型的边坡稳定性分析.docx

间歇性强降雨下基于Green-Ampt入渗模型的边坡稳定性分析 苏永华 李诚诚摘 要：研究了间歇性强降雨作用下边坡稳定性. 间歇性强降雨下湿润锋以上土体将会经受干湿循环作用，考虑干湿循环对边坡土体渗透特性及土水特性的影响，改进了传统Green-Ampt入渗（GA）模型;考虑干湿循环对土体强度的劣化作用，建立间歇性强降雨下边坡稳定性表达式.对于不考虑间歇性强降雨作用的边坡，改进GA模型可简化为传统GA模型，说明传统GA模型是改进GA模型的一个特例;利用改进GA模型和稳定性评价方法得到的结果与工程实际情况基本吻合，证明了该方法的可靠性;根据降雨入渗速率的变化特征，降雨过程可分为稳定阶段、持减阶段和突变阶段等3个阶段;传统GA模型与改进GA模型相比在失稳深度和时间上有较大延缓;仅改进稳定评价方法的传统GA模型与改进GA模型相比在失稳深度上一致，在失稳时间上有延缓.关键词：间歇性强降雨;Green-Ampt入渗模型;渗透特性;土水特性;强度衰减;稳定性：P642.22 文献标志码：ASlope Stability Analysis Based on Green-Ampt InfiltrationModel under Intermittent Heavy RainfallSU Yonghua？覮，LI Chengcheng（College of Civil Engineering，Hunan University，Changsha 410082，China）Abstract：The stability of slopes under intermittent heavy rainfall is studied. The soil above the wetting front is subjected to dry-wet cycles under intermittent heavy rainfall. Considering the influence of dry-wet cycles on the soil permeability and soil-water characteristics of slope，the traditional Green-Ampt infiltration （GA） model is improved. The expression of slope stability coefficient under intermittent heavy rainfall is established considering the deterioration of dry-wet cycles on soil strength. For slopes without considering intermittent heavy rainfall，the improved GA model can be simplified to the traditional GA model，which indicates that the traditional GA model is a special case of the improved GA model. The results obtained by using the improved GA model and the stability evaluation method are basically consistent with the actual engineering situation，which proves the reliability of the method. According to the variation characteristics of rainfall infiltration rate，the rainfall process can be divided into three stages： stable stage，continuous reduction stage and mutation stage. The slope infiltration depth and instability time of traditional GA are obviously behind that of the improved GA model. Compared with the improved GA model，the traditional GA model which only improves the stability evaluation method is consistent in the instability depth and delayed in the instability time.Key words：intermittent heavy rainfall;Green-Ampt infiltration model;permeability;soil-water characteristics;strength deterioration;stability隨着我国经济的高速发展，高陡型路堑边坡日益增多，加之极端天气多发，边坡稳定性问题日益突出.强降雨是诱发边坡失稳的最主要、最常见因素之一. 据统计，1980年以来，我国大陆所发生的大型灾害性滑坡约50%由强降雨直接触发，造成大量人员伤亡，经济损失巨大[1]. 因此，建立行之有效的强降雨作用下边坡失稳预测预报方法，对防灾减灾具有非常重要的理论意义和应用价值.基于降雨入渗模型的边坡稳定分析是评价降雨时边坡稳定性的有效方法.基于毛细管理论提出的Green-Ampt土壤水分入渗（GA）模型[2]，因其参数少、物理意义明确、精度高、适用性强、求解简单而在降雨型滑坡研究中逐渐得到应用与发展. 但传统GA模型仅适用于均质且初始含水率均匀分布的土壤水分入渗过程，因此，国内外不少学者对该模型的适用条件和范围以及主要参数的确定等问题进行了大量研究与改进，从而使该模型不断得到完善和发展. 常金源等[3]以GA模型为基础，在考虑动水压力的基础上分析了降雨入渗条件下浅层边坡稳定性.石振明等[4]提出适用于多层土边坡降雨入渗的GA模型，并采用強度折减法分析了不同降雨强度和历时下边坡稳定. Sung[5]基于初始含水率非均匀分布的GA模型提出了浅层基岩边坡稳定性分析方法. Loiciga等[6]基于GA模型与径流波动方程建立了入渗-径流耦合模型，并探讨了该模型下的边坡稳定性. Yao等[7]提出了GA模型和分层假定下的边坡稳定性分析方法.由于间歇性强降雨将导致坡内土体处于饱和-非饱和干湿循环过程，因此有必要用发展动态的观点研究间歇性降雨对边坡稳定性的影响.然而目前GA模型应用于边坡稳定性研究时主要考虑持续性降雨下边坡稳定性，鲜有利用GA模型研究间歇性降雨下边坡稳定性问题的报道，这限制了GA模型在边坡稳定性分析中的进一步应用.本文首先提出一种适用于间歇性强降雨的入渗模型，该模型以 GA 模型为基础，考虑了间歇性强降雨引发的干湿循环效应对饱和层土体渗透特性及土水特性的影响;然后考虑干湿循环对土体强度的衰减作用，利用双强度折减法分析边坡稳定性;最后分析了本文模型与GA模型的关系，并通过与实际工程情况对比，验证了本文计算方法的可靠性.1 降雨入渗模型Mein等[8]将GA 模型应用到降雨入渗情况，认为降雨入渗过程应分为降雨强度控制阶段和土体入渗能力控制阶段.如图1所示，假设强降雨强度为q，饱和体积含水率为θs，初始体积含水率为θi，边坡倾角为β. 降雨初期，坡表土体处于非饱和状态，土体入渗能力大于降雨强度，降雨全部渗入土体，此时，垂直于坡面方向的降雨入渗速率为：以每次强降雨湿润锋所达到的最大深度边界线，对土体进行分层.记第1次强降雨时湿润锋所达到的最大深度为hs，m1，第2次强降雨时湿润锋所达到的最大深度为hs，m2，第N次强降雨时湿润锋所达到的最大深度为hs，m N，则土层1的深度范围为h∈[0，hs，m1]，土层2的深度范围为h∈（hs，m1，hs，m2]，土层N的深度范围为h∈（h s，m（N-1），hs，m N]，其余土层为初始土层，其具体分层如图3所示. 由于干湿循环由强降雨引起，故第N次强降雨时，地表土体渗透特性参数与土水特性参数采用土体经历的n次干湿循环之后相应的参数，对土层1：n = N - 1，对于土层2：n = N - 2，其余依次类推，且n ≥ 0.当n = 0时，表示该层土体为未受降雨影响的初始土层，即此次降雨前雨水入渗深度未到达该土层.假设降雨强度控制阶段转为土体入渗能力控制阶段发生在h < hs，m1范围内，令hs，N、hsp，N和tp，N分别为第N次强降雨时湿润锋深度、临界时刻湿润锋深度和临界时刻，故第N次强降雨时土体入渗能力控制阶段的降雨入渗速率式（2）可修改为：式（26）反映了某次强降雨下湿润锋深度随降雨历时增长的动态变化规律，考虑了间接性强降雨对边坡土体渗透特性及土水特性的影响，推广了GA模型的应用范围.2 边坡稳定性分析目前边坡稳定性分析方法主要有极限平衡法[12]和强度折减法[13]. 降雨入渗条件下，非饱和土边坡多发生平行于边坡表面的浅层滑坡，破坏面往往发生在湿润锋或相对隔水层处.本文以无限长边坡体为研究对象，假定湿润锋以上土体处于饱和状态，并以湿润锋处作为潜在滑动面，在考虑干湿循环对土体力学特性的衰减作用下，利用双强度折减法评价边坡稳定性随强降雨湿润锋深度、强降雨次数的动态变化过程，其计算简图如图4所示.3 模型验证与分析3.1 与传统GA模型解析解比较传统GA模型不考虑间歇性降雨对边坡土体渗透特性及土水特性的影响，即令改进GA模型强降雨次数N = 1，则式（26）可化简为：3.2 案例计算与分析广西百隆路某段右侧高10 m、坡率1 ∶ 1.5的新修路堑边坡在经历第1个雨季时发生浅层坍滑.土体设计参数见表1.土体土水特征曲线函数表达式、有效黏聚力c′、有效内摩擦角φ′和孔隙比与干湿循环次数n满足的函数关系分别如式（35）～（39）所示.假设当地平均气温25 ℃，每次强降雨历时为4 h，降雨强度 q = 0.000 8 cm/s.不同干湿循环次数下，土体脱湿曲线和吸湿曲线的偏移系数见表2.图5所示为n = 0和n = 10土水特征曲线.同文献[11]所反映的规律相同，该土体不同干湿循环下土水特征曲线均存在明显的滞后现象，且随着干湿循环次数的增加，滞回效应明显减小，并趋于稳定;若干次干湿循环后土体土水特征曲线中脱湿曲线的偏移量大于吸湿曲线的偏移量.图6所示为不同干湿循环次数下，土体吸湿过程湿润锋处基质吸力与土层饱和导水率变化曲线.可知，土体吸湿时湿润锋处基质吸力随着干湿循环次数的增加而减小，且首次循环降幅最大，多次循环后趋于稳定.饱和导水率随着干湿循环次数的增加而增大，但增大速率逐渐放缓.土体吸湿时湿润锋处基质吸力和饱和导水率的变化皆表明了随着干湿循环次数的增加，土体逐渐松散，土体持水能力逐渐降低，符合相关试验研究和自然条件下土体实际渗流特性[18].根据降雨入渗速率的变化特征，将降雨过程分为稳定阶段、持减阶段和突变阶段等3个阶段，如图7所示. 但当N = 1时，只有稳定阶段、持减阶段2个阶段;不同降雨次数下每个阶段开始时间和持续时间也并不相同.阶段的入渗速率与强降雨次数无关，保持为qcos β.某次降雨下，该阶段湿润锋深度与降雨历时呈正比，边坡稳定性系数随降雨历时的增长快速下降.不同强降雨次数下，该阶段持续时间随强降雨次数的增加而增加，其原因在于干湿循环过程中，饱和。