灌溉对西北干旱荒漠区排土场边坡稳定的影响—以乌海市新星煤矿为例
Influence of irrigation on slope stability of dump in arid desert area of Northwest China: a case study of Xinxing coal mine, Wuhai
通讯作者:
收稿日期: 2018-08-15
Received: 2018-08-15
作者简介 About authors
邱驰(1992—),男,硕士生,从事岩土工程研究.orcid.org/0000-0002-5444-6265.E-mail:
为了揭示我国西北干旱荒漠化地区的煤矿排土场边坡在灌溉作用下的流固耦合力学机制,采用渗流力学理论与岩土力学理论相结合的方法,引入Darcy定律和质量守恒定律进行公式推导以建立流固耦合分析方法. 以乌海市新星煤矿为例,采用ABAQUS有限元软件构建边坡非饱和土的流固耦合稳定性分析模型,根据植被需水量以及降雨强度选定5、20、45、70 mm/h这4种等效灌溉强度并交叉搭配2、9、15 h这3种持续灌溉时间,形成12套灌溉方案进行数值模拟. 将数值结果与已有试验结果进行对比,结果表明,当入渗水量增加,边坡竖向变形位移、孔隙水压力分布、等效塑性应变以及安全系数均往不利于边坡稳定的方向发展. 当等效灌溉强度为70 mm/h、灌溉持续时间为15 h时,边坡竖向位移出现大幅度增长,安全系数急剧降至1.03. 在满足植被生长需水量的条件下,应控制灌溉强度,采用短持续时间、间隔灌溉,减少坡面径流出现的可能性,同时应做好截排水和坡面防护措施.
关键词:
The seepage mechanics theory was combined with the geomechanics theory to reveal the fluid-solid coupling mechanical mechanism of slopes on the coal mine dump under irrigation in an arid desert region of Northwest China. The Darcy law and the mass conservation law both were introduced to derive the fluid-solid coupling analysis method. Taking Xinxing coal mine of Wuhai city as a sample, the finite element software ABAQUS was used to establish the slope stability model of unsaturated soil under the fluid-solid coupling theory. Twelve types of irrigation schemes were used in the numerical simulation by choosing four types of equivalent irrigation intensity (i.e. 5 mm/h, 20 mm/h, 45 mm/h and 70 mm/h) and three sorts of continuous irrigation time (i.e. 2 h, 9 h and 15 h) according to the water requirement of vegetation and the rainfall intensity, respectively. The numerical results were compared with the existing experimental results. Results showed that the four indicators, including vertical deformation of slope, pore water pressure distribution, equivalent plastic strain, and safety factor, developed in the direction which was not conducive to the slope stability when the infiltration water increased. When the equivalent irrigation intensity was 70 mm/h and the irrigation time was 15 h, the vertical deformation of the slope increased significantly and the safety factor decreased rapidly to 1.03. Under the condition of meeting the water requirement of vegetation growth, the irrigation intensity should be controlled, and short duration and interval irrigation should be adopted to reduce the possibility of slope runoff. Simultaneously, the measures of interception and drainage as well as slope protection should be taken.
Keywords:
本文引用格式
邱驰, 黄建坤, 胡雨村, 郭小平, 胡云卿.
QIU Chi, HUANG Jian-kun, HU Yu-cun, GUO Xiao-ping, HU Yun-qing.
为了对干旱矿区排土场边坡进行生态修复,须对边坡进行灌溉补水. 灌溉水入渗土体会对边坡稳定性产生复杂的影响,水是诱发滑坡最关键的因素之一[6]. 在灌溉初期,水分浸入土体表层,改变表层土体孔隙率,增大土体中水的质量分数[7-9],进而改变边坡土体的物理力学性质;针对公路边坡的人工降雨模拟试验表明,降雨入渗导致边坡土体进行吸湿过程,随着水分浸入土体内部,土的抗剪强度不断降低,对边坡稳定性产生较大影响[10];在下渗过程中,在土体内部,水-气分界面承受着大于水压力的空气压力,这种压力差即基质吸力,是不可忽略的重要因素,既影响土的渗透系数,又影响土体饱和度,可同时作用于应力场和渗流场,构架起应力场和渗流场共同作用的桥梁[11-12]. 因此,基于渗流力学理论和岩土力学理论,应用ABAQUS有限元软件,通过流固耦合分析,针对位移变形、孔隙水压力变化、等效塑性应变区的延伸以及安全系数进行讨论,解析入渗水引起的边坡坍塌现象[13-14].
已有的基于流固耦合理论的坡体稳定性研究主要针对中国中东部地区自然降雨情况下的斜坡. 我国西北干旱矿区排土场是人工堆积体,坡度高坡角大,松散性离散程度大,降雨量极低,人工灌溉是保证边坡生态效益的主要措施之一,针对此类人工堆体在渗流作用下的边坡稳定性的研究较少. 灌溉水的引入,使边坡安全与生态效益之间出现矛盾. 在此类复杂地质边坡中,灌溉水量的安全阈值是多少,如何在保证边坡安全的同时保障此类生态脆弱区的生态效益,目前尚缺乏科学依据.
本研究以内蒙古乌海市新星煤矿第4排土场为研究对象,通过将渗流理论与岩土力学理论相结合,围绕渗流场与应力场的相互影响机理,建立平衡方程,展开力学响应机制研究. 应用ABAQUS有限元软件选定喷灌强度并搭配不同灌溉时间进行数值分析. 为了避免灌溉水量过低导致植被无法正常生长,以及过量灌溉导致坡面侵蚀、滑坡等灾害,讨论灌溉强度以及灌溉时间组合对边坡稳定性的影响,通过数值计算得到关键节点位移、孔隙水压力、等效塑性应变和安全系数的变化规律,进而揭示边坡失稳破坏机理,为干旱荒漠化地区煤矿排土场研建适合的灌溉方案提出理论依据,同时也为降雨对边坡稳定性的影响提供参考.
1. 非饱和土渗流理论
1.1. 渗透系数函数
在非饱和土中,基质吸力是影响渗透系数以及饱和度的重要因素,Genuchtenm[20]提出经典的基质吸力VG模型:
式中:w为水的质量分数;ws为饱和状态下水的质量分数;wr为残余的水的质量分数;
非饱和土渗透系数表达式为
式中:K为非饱和土渗透系数,Ks为饱和土渗透系数.
1.2. 应力场对渗流场的影响机理
在多孔介质中,针对非饱和土体,土体孔隙比n与渗透系数K成正比例关系,K=K(n). 在自然条件下,土体的孔隙比为n0,经过应力场作用,施加荷载,土体产生一定挤压或拉伸,得到体积应变εv. 假定体积应变是由于孔隙度变化所引起的,则应力场作用后的孔隙比n=n0+εv. 应力场作用于土体,可产生一定的土体体积应变,从而改变土体孔隙比,K与应力场σij可建立函数关系,即K=K(σij). 因此,应力场通过影响孔隙比,从而影响土体的渗透系数,进而对渗流场产生影响.
1.3. 渗流场对应力场的影响机理
渗流场分布的变化将改变土体水分分布的变化,渗流场通过改变水荷载,即通过改变应力场的外荷载形式来影响应力场分布. 土体中的水头分布为h(x,y),作用面的渗透水压力p=γw(H−y),y为水头位置,H为设置的边坡初始地下水总深度,
在有限元计算中,将渗透体积力和渗透水压力转化为单元的节点外荷载. 渗流场对应力场的影响可以简述为,通过改变土体单元渗透压力和渗透体积力来影响作用于应力场的外荷载.
1.4. 流固耦合分析控制方程
根据质量守恒定律,水分在土体单元入渗流动的增减速率,等于其在该土体单元中体积流量速率之差. 根据达西定律,假设单元体为常数且为各向同性,x、y方向的渗透系数相同,即Kx=Ky,均质土体和均质水分子不可压缩,压缩体积Ss=0,得到渗流场微分方程表达式如下:
基于水力学基本原理,渗流体积力与水力梯度呈正比例关系. 采用有限元法,将其转化为单元节点外荷载进行叠加计算:
式中:FS为渗流体积力所引起的等效单元节点力矩阵函数,ΔFS为渗流体积力增量所引起的等效单元节点力矩阵函数,Ω为全部有限单元土体,N为引入的形矩阵函数.
设定位移边界和渗流边界条件,联立渗流场和应力场,得到土体在应力场与渗流场共同作用下的流固耦合控制方程,进行有限元计算求解:
式中:K为各向渗透系数(应力场矩阵函数),Δδ为位移增量矩阵,ΔF为节点荷载增量矩阵,
2. 数值分析
2.1. 研究区概况
西北干旱荒漠区位于我国生态安全格局“两屏三带”的北方防沙带,该区域生态环境脆弱,气候干旱、风力较大、植被退化、土地沙化严重. 区域干旱灾害程度为重旱. 受到煤炭资源开发的影响,地表挖损与塌陷、土地荒漠化和植被退化加剧,恢复难度较大. 研究区位于内蒙古乌海市海勃湾区境内的新星煤矿,选择与西北干旱荒漠区其他煤矿排土场地形地貌以及土质状况相近的第4排土场边坡,具有典型性和代表性. 此区域已展开大规模的边坡覆绿工程,并采取节水灌溉等一系列抚育手段进行边坡生态修复. 本研究针对灌溉水入渗对边坡稳定的影响,可以为此类地区的灌溉覆绿措施提供参考和指导.
2.2. 边坡计算模型
采用无人机和遥感技术调研收集研究区的影像资料,对其进行处理和解译,结合实地观测,获取研究区地形地貌的特征信息,如图1(a)所示. 研究区域属于干旱荒漠地区,生态环境脆弱,植被退化严重,年均降雨量仅为159.8 mm,不利于植被自然恢复,须通过灌溉进行人工干预. 通过对排土场边坡的实地调查,在修复区设置随机样地,在样地中随机设置样方应用土钻分别采集0~10 cm、10~20 cm、20~30 cm的土壤样品,通过固结不排水(consolidated undrained,CU)直剪试验测定土壤物理力学性质[21]. 人工堆体的物理参数离散性较高,因此测量物理参数平均值,并结合文献资料确定参数,如表1所示. 表中,ρd为干密度,μ为泊松比,E为弹性模量,c为黏聚力,φ为内摩擦角. 边坡的几何尺寸通过实测取得,坡高为20 m,坡角为39.5°,详细尺寸如图1(b)所示. 图中,v为表面孔隙流速.
图 1
图 1 内蒙古乌海市海勃湾区境内新星煤矿第4排土场
Fig.1 Fourth dump of Xinxing coal mine, Haibowan district, Wuhai, Inner Mongolia
表 1 黄土的物理力学参数
Tab.1
参数 | 数值 | 参数 | 数值 | |
ρd/(g∙cm−3) | 1.3 | c/kPa | 8.5 | |
μ | 0.3 | φ/(°) | 30 | |
E/MPa | 17.5 | K/(mm·h−1) | 18 |
实地调研发现,在过量灌溉及强降雨的情况下,研究区域出现多处侵蚀沟甚至所覆表层黄土出现浅层滑坡. 所覆黄土介质相对均匀,且研究表明当入渗水量较大时,试验尺度的变化对最大入渗深度和土壤水流运动非均匀程度的影响不明显[22]. 因此,统一假定土体为匀质黄土,不考虑地震耦合和植被蒸腾作用,将坡脚处作为初始地下水位处(即地下水位与地表持平),坡脚外表土层顶面的孔压为0 kPa. 在左右两侧水位线以下的边界处起始,随地下水位深度不断增加,应用式(7)设置随深度线性增加的静水孔压,即坡底最大孔压为100 kPa. 在模型底部设置水平和竖直位移约束,即固定约束;在模型左右两侧设置水平位移约束. 以往针对渗流自由面问题的数值模拟,多假定模型边界为不透水边界,试验与数值结果的比较证明,在此假设下数值解有较高的精度[23]. 因此,模型底部和两侧边界设为不排水边界. 利用式(1)、(2),定义渗透系数和饱和度之间的关系以及孔压和饱和度之间的关系,并利用式(7)定义静水孔压和水头位置的关系.
式中:
由于西北干旱荒漠区排土场灌溉缺乏科学参数,根据植被全生育期需水量[26]及当地降雨状况,选择5.0、9.0、17.5 mm/h这3种灌溉强度(分别对应中雨2.6~5.0 mm、大雨8.0~16.0 mm、暴雨16.0~50.0 mm),讨论灌溉强度对边坡稳定性的影响. 考虑到坡面径流,根据水文地质中给排水领域的推理计算方法(rational method),可通过公式Qp=CipA(Qp为汇水总量;C为流失系数,取0.5[25];ip为灌溉强度;A为汇水面积)估算总汇水量[27]. 水分从各级坡面汇流至第1级台阶后,灌溉强度相当于增加到20、45、70 mm/h. 在低灌溉强度下,(雨)水存在完全渗透的可能性,因此将5 mm/h的灌溉强度作为对照组. 根据现代灌溉原理[26],选择2 h作为灌溉时间. 研究区域在2018年发生过百年一遇的暴雨,导致边坡出现侵蚀沟以及浅层滑坡现象. 考虑产流汇流滞后效应,延长灌溉持续时间至9、15 h,讨论灌溉持续时间对边坡稳定的影响. 在坡顶设置5、20、45、70 mm/h这4种灌溉强度,同时交叉搭配2、9、15 h这3种持续灌溉时间,共计进行12组仿真试验. 根据不同灌溉强度,利用式(3)计算相应的渗透压力以及渗透体积力,并作用于应力场.
2.3. 边坡稳定性分析
通过静力分析,对全域模型施加重力荷载并施加孔隙水压边界,边坡土体饱和度分布如图2所示. 图中,SAT为饱和度. 坡脚以下饱和度约为1.00,坡脚以上土体饱和度迅速缩小至约0.08,在坡脚上下两侧出现明显差值. 表明在此处,土壤由饱和逐渐转变为非饱和,此处为过渡的分界线,即为初始水位线. SAT可利用式(2)进行计算.
图 2
2.3.1. 竖向变形位移
针对边坡模型,采用5、20、45、70 mm/h这4种灌溉强度,并交叉搭配2、9、15 h这3种灌溉时间,共计进行12组仿真试验. 在灌溉时,平台以及坡体表层土吸收水分,饱和度增加、土体黏聚力降低,致使土颗粒间摩擦力减小、土壤阻力降低,土体将产生一定的沉降. 灌溉对边坡的水平变形位移影响不大,故只讨论竖向变形位移的变化规律. 选取坡脚处拐点作为关键监测点,监测点的竖向位移随灌溉持续时间的变化如图3所示. 图中,T为持续灌溉时间,U为坡脚处监测点的竖向位移. 当灌溉时间为2 h时,在灌溉强度分别为5、20、45 mm/h的情况下,竖向位移均平缓增长,当灌溉强度提升至70 m/h时,竖向位移出现突变并首次呈现显著增长. 在灌溉时间为9 h或15 h的情况下,当灌溉强度为5、20、45 mm/h时,边坡竖向变形位移均呈缓慢增长的趋势,当等效灌溉强度为70 mm/h且灌溉时长达到9 h时,竖向位移出现突变,且急剧增长,累积竖向位移远高于灌溉强度为5、20、45 mm/h时的竖向位移. 此时,水分不仅局限于表层土,而是开始沿竖直方向向下扩散,导致边坡整体水的质量分数升高,基质吸力减小,土体黏聚力减小,位移急剧增长. 边坡由最初的表层土收缩发展为边坡整体收缩,在坡脚处产生堆积甚至隆起,诱发边坡失稳.
图 3
图 3 坡脚节点处的竖向位移随灌溉时间的变化曲线
Fig.3 Curve of vertical displacement of joints at slope foot with different irrigation time
在图3中,在相同的灌溉强度下,持续灌溉时间由2 h延长至9 h,再延长至15 h,坡角位移总体均呈扩张趋势进行增长. 当灌溉强度为5、20、45 mm/h时,延长灌溉时间,边坡位移均呈现小幅度的缓慢增长;当等效灌溉强度为70 mm/h时,灌溉时间由9 h延长至15 h,边坡竖向位移增幅为104.7%. 灌溉时间的延长导致总灌溉水量的持续增加,边坡土体内部水的质量分数升高,单元土体内部水分过饱和,继而抗剪强度降低,土体内部持续收缩,边坡位移进一步扩大. 在满足植被正常生长的条件下,当灌溉强度为5 mm/h、灌溉持续时间为2 h时,边坡处于稳定状态. 当考虑地表径流时,在强灌溉(降雨)条件下,总入渗水量过多,即当等效灌溉强度为70 mm/h时,灌溉持续时间若达到15 h将会对边坡整体位移变形产生较大的影响.
2.3.2. 孔隙水压力
随着灌溉水浸入边坡土体内部,孔隙水压是判断边坡稳定的重要指标之一[28]. 在水的质量分数较高的土体中,孔隙水压力随时间累积而升高[29]. 孔隙水压力不断升高将导致岩土体内部基质吸力减小,抗剪强度降低,诱发边坡失稳破坏[30]. 研究区属于干旱荒漠环境,模拟发现在相同灌溉强度下延长较短的时间,对孔隙水压力分布的影响不大. 如图4(a)~(d)所示分别为5、20、45、70 mm/h的灌溉强度下产生的孔隙水压力. 图中,a为坡高、b为坡宽. 经过灌溉水的初期入渗,土体中水的质量分数逐渐增大. 产生新的孔隙比,从而出现孔隙水压力. 水的质量分数的分布随边坡土层总入渗水量变化[31]. 当灌溉强度为5 mm/h时,入渗水量较少,表层水的质量分数高于内部水的质量分数,单位时间内对土体孔隙比的改变较小,从而仅在坡面浅层产生一定的孔隙水压力,并且孔隙水压力基本呈线性分布,如图4(a)所示.
图 4
图 4 不同灌溉强度下边坡不同位置处的孔隙水压力
Fig.4 Porewater pressure at different positions of slope under different irrigation intensities
2.3.3. 等效塑性应变及贯通区
如图5(a)、(b)所示分别为5、70 mm/h灌溉强度下9 h持续灌溉所产生的等效塑性应变贯通区;如图5(c)、(d)所示分别为5、70 mm/h灌溉强度下15 h持续灌溉所产生的等效塑性应变贯通区. 图中,PEEQ为等效塑性应变. 由图5可以看出等效塑性应变的发展趋势. 在相同的灌溉时间下,随着喷灌强度的提高,等效塑性应变均先在坡脚出现,并呈流线型,在边坡内部沿坡面逐渐延伸形成贯通区. 提高喷灌强度,等效塑性应变贯通区呈现扩张趋势,即坡脚处等效塑性应变区增大,且沿坡面产生一定程度的延伸. 在相同的喷灌强度下,随着灌溉时间的延长,等效塑性应变同样沿坡面呈现扩张趋势. 相比于改变灌溉强度,在低灌溉强度下,随着时间的延长,等效塑性应变扩张较缓慢,表明降水入渗的渗流体积力带来的单元节点应力增长缓慢,减缓了单元体的应力传递. 在15 h持续灌溉后,坡脚处出现塑性应变,但未出现连续贯穿滑动面,表明边坡相对稳定.
图 5
图 5 不同灌溉强度和时间下的等效塑性应变区
Fig.5 Equivalent plastic strain region under different irrigation intensities and time
2.3.4. 安全系数
在利用强度折减法计算安全系数时,可以通过位移突变和迭代计算不收敛计算终止2种方法取值[34]. 本研究以不收敛点作为安全系数. 如图6所示,在5 mm/h灌溉强度下持续灌溉9、15 h对应的安全系数分别为1.28、1.27;在20 mm/h灌溉强度下持续灌溉9、15 h对应的安全系数分别为1.19、1.16;在45 mm/h灌溉强度下持续灌溉9、15 h对应的安全系数分别为1.15、1.13;在70 mm/h灌溉强度下持续灌溉9、15 h对应的安全系数分别为1.14、1.03. 通过与规范进行对比,70 mm/h灌溉强度下持续灌溉15 h后的安全系数小于规范的最小值,边坡存在失稳破坏的风险. 如图6所示,其余7组均在折减系数为1.10时出现突变,在Fs′=1.10之后出现位移陡增,验证了规范建议安全系数的合理性.
图 6
图 6 折减系数随位移的变化关系
Fig.6 Relationship between reduction coefficient and displacement
对比持续灌溉时间分别为9、15 h的2组曲线可以看出,在相同的持续灌溉时间内,提高喷灌强度,边坡安全系数明显降低. 提高灌溉强度,单位时间内边坡土体内部水的质量分数增大,抗剪强度降低,影响边坡稳定性. 对比灌溉强度分别为5、20、45、70 mm/h的4组曲线可以看出,在相同的灌溉强度下,延长灌溉时间,会导致边坡内部的入渗水增多,边坡土体中水的质量分数增加. 非饱和区岩土体内部孔隙水压力增大导致岩土体基质吸力减小,抗剪强度减小,同时增加边坡的下滑力,减小土颗粒之间的摩擦,降低土壤阻力. 因此,在延长灌溉时间后边坡安全系数明显降低,对边坡稳定性造成不良影响. 当入渗水浸入土体内部,会减少土体内部孔隙,致使土体自身重力增加,从而诱发边坡失稳坍塌破坏.
为了进一步说明安全系数计算的准确性,采用极限平衡理论,在SLOPE/W软件上利用简化BISHOP法进行安全系数的核算,反复迭代直至前后2次计算结果差值接近,即计算结果收敛,最终得到安全系数[35]. 简化BISHOP法和SSR法的计算结果如表2所示. 表中,5 mm/h,9 h是指以5 mm/h的灌溉强度持续灌溉9 h,其他类似. 可以看出,简化BISHOP方法计算得到的安全系数均略低于SSR法的计算结果,平均相对差为2.20%. BISHOP法操作简单,应用较多,但是仅能满足部分平衡条件. SSR法能充分反映复杂实际地质情况,结果更接近实际. 可以看出,增加灌溉强度和延长喷灌时间,位移、等效塑像应变贯通区、孔隙水压力等均向不利于边坡稳定的方向发展,特别是灌溉强度的增幅,会对边坡产生较明显的影响.
表 2 简化BISHOP法和SSR法计算结果的对比
Tab.2
工况 | Fs | |
BISHOP | SSR | |
5 mm/h,9 h | 1.27 | 1.28 |
5 mm/h,15 h | 1.25 | 1.27 |
20 mm/h,9 h | 1.17 | 1.19 |
20 mm/h,15 h | 1.13 | 1.16 |
45 mm/h,9 h | 1.12 | 1.15 |
45 mm/h,15 h | 1.10 | 1.13 |
70 mm/h,9 h | 1.12 | 1.14 |
70 mm/h,15 h | 1.01 | 1.03 |
为了保证植被的需水量以及考虑到节水目的,控制灌溉强度起到至关重要的作用. 因此,建议通过控制流速和持续灌溉时间来限制灌溉水总量. 对比计算结果,增加流速和灌溉时间后,位移、孔隙水压力以及等效塑性应变均呈增大趋势. 特别是在提高灌溉流速后,边坡整体位移产生明显增幅.
在地表渗透性能良好的情况下,在5 mm/h的灌溉强度下持续灌溉2 h,边坡处于稳定状态,不会发生失稳破坏. 然而,研究区域广泛存在土壤表层物理结皮现象,此现象极易引发径流灾害. 在20、45 mm/h的等效灌溉强度下持续灌溉9、15 h,边坡依旧处于稳定状态. 在70 mm/h的等效灌溉强度下持续灌溉15 h,有引起边坡失稳的风险,这是由于强灌溉(或暴雨)导致上级平台出现径流,在最底层的一级平台形成汇流,对坡体稳定产生影响.
在满足植被全生育期生长需水量的条件下,应控制灌溉强度,采用短持续时间、间隔灌溉,减少坡面径流出现的可能性. 当采用高强度灌溉时,在坡面径流汇流后,会有暴雨级别的雨量入渗坡体,对边坡造成极大的影响. 此外,坡顶和坡面须做好截排水措施,及时疏导多余水量,避免强灌溉(降雨)导致的径流在底部坡面形成汇流.
2.4. 试验与数值案例对比
本研究对象为人工堆体,土体物理力学性质离散度较高,考虑到实地灌溉过程会受到自然条件如风力、蒸发等众多因素的干扰,试验操作难度较大. 本研究通过与室内模型试验[36]的对比,说明数值计算结果的准确性. 根据试验测得的土壤物理、力学参数以及模型试验的几何参数建立数值模型,选取灌溉强度为70 mm/h、持续灌溉时间为1 h,验证非饱和土体在水分下渗过程中对边坡的稳定性影响.
在数值模拟过程中,边坡土体孔隙水压上部呈现负值,说明该区域为非饱和区,同时存在饱和渗流和非饱和渗流现象. 通过对比试验和数值模拟的浸润线位置,即孔隙水压为零的等值线位置,验证数值结果的准确性. 如图7所示,试验和数值浸润线的移动趋势相同,均随水分的下渗,由坡顶开始沿斜坡表层下移,直至坡脚处,呈现堆积情况. 试验从侧面直接观测最外层坡面的土壤浸润情况,无法获得内部分层入渗规律. 数值模拟通过孔压分布等值线区分负孔压,即非饱和区域,因此,通过数值模拟能更直观地观测到浸润面位置以及入渗分层情况. 在数值计算中未考虑边坡浅层坡面的水分蒸发,因此,相比于试验结果,数值计算得到的浸润面位置更贴近浅层坡面.
图 7
模拟持续灌溉4 h,等效塑性应变分布如图8所示. 图中,深色阴影区域表示数值计算得到的塑性应变贯通区,该区由坡脚延伸至坡顶,指示滑坡面的产生位置. 将试验得到的破坏面、坍塌体、沉降区和堆积区与数值计算得到的等效塑性应变贯通区进行对比,吻合较好,说明数值结果和试验结果的边坡失稳破坏时间和破坏形式一致. 等效塑性应变从坡脚处开始延伸,验证了边坡失稳破坏是由坡脚处最先开始的.
图 8
3. 结 论
(1)通过数值分析可以看出,入渗水极易存储于边坡浅层区域,水分集中分布在边坡表层,对边坡底部的影响较小;延长灌溉持续时间,入渗水量增多,边坡内部易形成积水,岩土体自身重力增加,从而诱发边坡失稳破坏.
(2)在灌溉过程中,延长灌溉时间,土体中水的质量分数增大,进而影响内摩擦角和抗剪强度. 针对此类浅层滑坡,总灌水量的增加是导致边坡安全系数降低的最关键的外在因素.
(3)在保证边坡稳定、满足植被需水和节水的前提下,应控制灌溉强度,采用短持续时间、间隔灌溉,减少坡面径流出现的可能性. 此外,坡顶和坡面须做好截排水措施,避免强灌溉(降雨)引发的径流在底部坡面形成汇流.
(4)边坡失稳首先从坡脚开始,沿斜坡从内部延伸,直至等效塑性应变区贯通,形成连续滑动面导致失稳破坏. 因此,边坡变形破坏易发生于坡脚,建议在坡脚处加强防护措施.
参考文献
植被护坡技术及其应用
[J].DOI:10.3969/j.issn.1672-3007.2004.04.022 [本文引用: 1]
Technology and application of vegetation slope protection
[J].DOI:10.3969/j.issn.1672-3007.2004.04.022 [本文引用: 1]
排土场空间效应及其稳定性评价方法研究
[J].
Researches on space effect of waste dumps and its stability evaluation method
[J].
矿山废弃地生态修复研究进展
[J].
Review on ecological restoration of mine wasteland
[J].
露天矿高陡边坡稳定性研究现状及发展趋势
[J].DOI:10.3969/j.issn.1000-7598.2011.05.025
Research situation of open-pit mining high and steep slope stability and its developing trend
[J].DOI:10.3969/j.issn.1000-7598.2011.05.025
公路岩质高边坡稳定性的综合评价研究
[J].DOI:10.3321/j.issn:1000-4548.2006.05.003 [本文引用: 1]
Comprehensive evaluation and research on stability of high lithological slope of expressway
[J].DOI:10.3321/j.issn:1000-4548.2006.05.003 [本文引用: 1]
降雨与滑坡灾害相关性分析及预警预报阀值之探讨
[J].DOI:10.3969/j.issn.1000-7598.2007.05.039 [本文引用: 1]
Discuss on the correlations between landslides and rainfall and threshold for landslide early-warning and prediction
[J].DOI:10.3969/j.issn.1000-7598.2007.05.039 [本文引用: 1]
Causative mechanisms of rainfall-induced fill slope failures
[J].DOI:10.1061/(ASCE)1090-0241(2004)130:6(593) [本文引用: 1]
土壤含水率对浅层滑坡体不同层次土壤抗剪强度的影响
[J].DOI:10.3321/j.issn:1009-2242.2009.06.012
Effects of soil water content on soil shearing strength to different soil layer of shallow landslide
[J].DOI:10.3321/j.issn:1009-2242.2009.06.012
边坡土壤含水率对不同植被土壤抗剪强度的影响
[J].DOI:10.3969/j.issn.1001-005X.2009.03.019 [本文引用: 1]
Effects of slope soil water content on soil shear strength of different vegetation
[J].DOI:10.3969/j.issn.1001-005X.2009.03.019 [本文引用: 1]
降雨入渗对泥岩-土混填路堤稳定性的影响
[J].DOI:10.3321/j.issn:1000-6915.2008.11.012 [本文引用: 1]
Effect of rainfall infiltration on stability of mudstone-soil mixture embankment
[J].DOI:10.3321/j.issn:1000-6915.2008.11.012 [本文引用: 1]
A hydro-mechanical elastoplastic model for unsaturated soils under isotropic loading conditions
[J].
Investigation of a pore pressure driven slope failure using a coupled hydro-mechanical model
[J].
Rainfall-induced slope collapse with coupled material point method
[J].
应用ABAQUS程序进行渗流作用下边坡稳定分析
[J].
Analysis of slope stability under seepage by using ABAQUS program
[J].
Modeling regional initiation of rainfall-induced shallow landslides in the eastern Umbria region of central Italy
[J].DOI:10.1007/s10346-006-0037-0 [本文引用: 1]
岩质边坡降雨入渗过程的数值模拟
[J].DOI:10.3321/j.issn:1000-6915.2003.04.023 [本文引用: 1]
Numeric simulation on rainfall infiltration on rock slope
[J].DOI:10.3321/j.issn:1000-6915.2003.04.023 [本文引用: 1]
强降雨作用下浅层滑坡尖点突变模型研究
[J].
Cusp-catastrophic model of shallow landslide under intensive rainfall
[J].
岩质边坡在降雨条件下的稳定性分析
[J].
Stability of rock slope under rainfall
[J].
强度折减有限元法中边坡失稳判据的研究
[J].
Slope failure criteria of shear strength reduction finite element method
[J].
A closed-form equation for predicting the hydraulic conductivity of unsaturated soils
[J].
长武塬边坡不同植被下土壤团聚体及入渗特征
[J].
Soil aggregates and infiltration characteristics under different vegetations in Changwu tableland slope of northwestern China
[J].
入渗水量和试验尺度对土壤水非均匀流动的影响
[J].
Effects of infiltration amount and experimental scale on heterogeneous soil water flow
[J].
Trefftz间接法解自由面渗流问题
[J].DOI:10.3969/j.issn.1007-4708.2005.03.008 [本文引用: 1]
Solution of tree surface seepage problems by Trefftz indirect method
[J].DOI:10.3969/j.issn.1007-4708.2005.03.008 [本文引用: 1]
流动沙漠地区灌溉林地盐结皮层特征的初步研究
[J].DOI:10.3321/j.issn:1000-1522.2007.02.007 [本文引用: 1]
Characteristics of salt crust layers in the forests irrigated with saline water in mobile desert
[J].DOI:10.3321/j.issn:1000-1522.2007.02.007 [本文引用: 1]
降雨型滑坡现场模拟试验研究
[J].
Research on on-site simulation test of landslide by artificial rainfall
[J].
Extending the rational method for assessing and developing sustainable urban drainage systems
[J].DOI:10.1016/j.watres.2018.07.022 [本文引用: 1]
降雨入渗条件下软岩边坡稳定性分析
[J].DOI:10.3969/j.issn.1000-7598.2012.08.019 [本文引用: 1]
Stability analysis of soft rock slope under rainfall infiltration
[J].DOI:10.3969/j.issn.1000-7598.2012.08.019 [本文引用: 1]
饱和土体固结3D比例边界有限元法分析
[J].
Modeling of saturated soil dynamic coupled consolidation problems using three-dimensional scaled boundary finite element method
[J].
Soil suction monitoring for landslides and slopes
[J].DOI:10.1144/1470-9236/09-010 [本文引用: 1]
土壤水分变异对降雨和植被的响应
[J].
Responses of soil moisture variations to rainfall and vegetation
[J].
黄土高原植物根系增强土体抗剪强度的模型与试验研究
[J].
Research on model and test of reinforcing shear strength by vegetation roots in the Loess Plateau of northern China
[J].
有限元强度折减法中边坡失稳的判据探讨
[J].DOI:10.3969/j.issn.1000-7598.2005.02.035 [本文引用: 1]
Study on slope failure criterion in strength reduction finite element method
[J].DOI:10.3969/j.issn.1000-7598.2005.02.035 [本文引用: 1]
地下水渗流作用下土坡稳定性的简化Bishop法解
[J].
Solution for stability of soil slopes subjected to seepage of groundwater via simplified Bishop method
[J].
/
〈 |
|
〉 |
