ISSN 1003-8035 CN 11-2852/P
  • 中国科技核心期刊
  • CSCD收录期刊
  • Caj-cd规范获奖期刊
  • Scopus 收录期刊
  • DOAJ 收录期刊
  • GeoRef收录期刊
欢迎扫码关注“i环境微平台”

降雨对花岗岩风化层路堑边坡滑动模式影响以福建云平高速云霄段为例

蔡荣坤, 戴自航, 徐根连, 胡长江

蔡荣坤, 戴自航, 徐根连, 胡长江. 降雨对花岗岩风化层路堑边坡滑动模式影响——以福建云平高速云霄段为例[J]. 中国地质灾害与防治学报, 2021, 32(2): 27-35. DOI: 10.16031/j.cnki.issn.1003-8035.2021.02.04
引用本文: 蔡荣坤, 戴自航, 徐根连, 胡长江. 降雨对花岗岩风化层路堑边坡滑动模式影响——以福建云平高速云霄段为例[J]. 中国地质灾害与防治学报, 2021, 32(2): 27-35. DOI: 10.16031/j.cnki.issn.1003-8035.2021.02.04
Rongkun CAI, Zihang DAI, Genlian XU, Changjiang HU. Influence of rainfall on sliding modes of cutting slope of weathered granite stratum: Taking Yunxiao section in the Yunping freeway in Fujian for example[J]. The Chinese Journal of Geological Hazard and Control, 2021, 32(2): 27-35. DOI: 10.16031/j.cnki.issn.1003-8035.2021.02.04
Citation: Rongkun CAI, Zihang DAI, Genlian XU, Changjiang HU. Influence of rainfall on sliding modes of cutting slope of weathered granite stratum: Taking Yunxiao section in the Yunping freeway in Fujian for example[J]. The Chinese Journal of Geological Hazard and Control, 2021, 32(2): 27-35. DOI: 10.16031/j.cnki.issn.1003-8035.2021.02.04

降雨对花岗岩风化层路堑边坡滑动模式影响——以福建云平高速云霄段为例

基金项目: 泉州市科技局项目(2017T013)
详细信息
    作者简介:

    蔡荣坤(1994-),男,福建漳州人,硕士研究生,助理工程师,主要从事边坡稳定性方面的研究工作。E-mail:1571131783@qq.com

    通讯作者:

    戴自航(1966-),男,湖南长沙人,博士,教授,主要从事边坡稳定分析、建筑基础、基坑支护等方面的教学研究。E-mail:dzhang@fzu.edu.cn

  • 中图分类号: P642.22

Influence of rainfall on sliding modes of cutting slope of weathered granite stratum: Taking Yunxiao section in the Yunping freeway in Fujian for example

  • 摘要: 降雨引起路堑边坡的浅层局部滑动是工程中常见的现象,但其内在原因的研究尚不充分。为此,以福建省云平高速云霄段花岗岩风化层中的某两级路堑边坡为例,采用饱和和非饱和土抗剪强度理论,按照强度折减有限元方法分析了该边坡在一般工况、地下水渗流工况和降雨入渗工况下的滑动模式。结果表明:一般工况下该边坡滑动模式为深层整体滑动;地下水渗流作用可使边坡由深层滑动模式向浅层滑动模式转化;长时小或中等降雨作用下坡体内出现三道滑动面,且一级边坡内滑动面最先贯通,滑坡模式由深层整体滑动向浅层局部滑动演化。短时强降雨作用下仅在一级边坡内出现贯通的滑动面,即由整体滑动模式完全转变为浅层局部滑动模式。降雨工况越不利表现越为明显,与实际滑坡现象相吻合,且安全系数不满足规范要求,从而从理论上验证了该路堑边坡发生浅层局部滑动的必然性。此外,为避免花岗岩风化层路堑边坡在降雨入渗下发生浅层局部滑动,对假设的不同放坡方案进行了数值模拟分析。结果表明,适当放缓坡率可使一级边坡内不出现滑动面,滑坡模式仍为较深位置的整体滑动,且安全系数可得到有效提高,使其满足规范要求。因此,实践中,采用文中方法对设计的类似地层中的路堑边坡事先进行降雨入渗影响下的模拟分析极为必要。将使放坡坡率更为科学合理,可在确保边坡整体稳定性的同时,极大降低其发生浅层局部滑动的可能性。
    Abstract: It is a common phenomenon that rainfall may cause locally shallow sliding of cutting slopes, but the study on the immanent cause is still insufficient. Take a two-grade cutting slope in the weathered granite stratum in Yunxiao section in the Yunping freeway in Fujian as an example, the shear strength theories of saturated soil and unsaturated soil and the strength reduction method of finite element were adopted to study the sliding modes of cutting slope under general working condition, groundwater seepage, and rainfall infiltration. The results show that the sliding mode of the slope exhibits the globally deep sliding under the general working condition. The sliding mode turns to be the shallow sliding from the deep sliding under groundwater seepage. Three slip surfaces may develop in the slope after a long-term light and/or moderate rainfalls, and the slip surface within the first-grade slope runs through firstly so that the sliding mode evolves into the locally shallow sliding from the globally deep sliding. A run-through slip surface only occurs within the first-grade slope after short-term strong rainfalls, i.e., the sliding mode completely turns to be the locally shallow one from the globally sliding one, and the more unfavorable the working condition of a rainfall is, the more remarkable it exhibits. It also agrees well with the phenomenon of the practical landslide, and the safety factor cannot meet the requirement of the specification, and thereby the inevitability of the locally shallow sliding of the cutting slope was theoretically verified. In addition, to avoid the occurrence of the locally shallow sliding of cutting slopes in the weathered granite stratum under the infiltration of rainfalls, various hypothetical schemes were numerically modeled. The results show that no slip surface will occur within the first-grade slope if it is sloped at an appropriate gentle slope, and as a result, the sliding mode is still a globally sliding in a deeper location and the safety factor may be effectively enhanced to meet the requirement of the specification. Therefore, in practice, it is extremely necessary to carry out the simulation analysis of a designed cutting slope in similar stratum under the infiltration of rainfall beforehand. Thus, the gradient of sloping will be more scientific and reasonable, so probability of the occurrence of the locally shallow sliding may be significantly reduced.
  • 地裂缝是一种在内外地质营力作用下,地壳浅表部土层中发生的破裂以及错断现象[1-7],表现形式多种多样,常见的有地表裂缝、连续陷坑、陡坎、缓坡、隐伏地裂缝形成的近地表破碎带等。地裂缝在国内外分布广泛,对其影响范围内的农田、房屋建筑、公路铁路、油气管道等线性工程以及地下工程可造成巨大破坏,严重影响了耕地以及建设用地的适宜性,给国民经济带来了巨大损失[8-16]

    东非大裂谷是目前地球表面最大的断层陷落带,具有典型的大陆裂谷特征,谷底地势平缓,边缘为相互平行的阶梯状断层群。受地质构造影响,裂谷区历史上火山活动异常活跃,断裂极为发育,地形高低起伏,工程地质条件异常复杂[17]。根据GPS观测数据,东支裂谷现今仍以0.5 cm/a的速率进行扩张,构造活动强烈[18]。岩浆侵入-地震活动-裂谷开裂是目前认为东非裂谷的主要活动模式。肯尼亚裂谷位于东非大裂谷东支的中段,构造活动和火山活动强烈,裂谷区岩层产生大量隐伏破裂,为地裂缝的产生提供了优势构造条件。肯尼亚裂谷地裂缝常在一场大雨后出露地表,与降雨侵蚀存在密切的相关性。裂谷区属于东非高原区,降雨具有短时降雨量大、频率高且极端降雨天气频发的特点。强烈的降雨可为地裂缝的形成提供有利的水力条件。裂谷区浅表部覆盖层多为第四纪火山沉积物,厚度较小,土质疏松,耐水力侵蚀能力较差,为地裂缝的形成提供了良好的物质条件。这种特殊孕育条件,使得该区域地裂缝较为发育,对当地居民的生活以及工程建设造成了不可忽视的影响,然而针对裂谷区地裂缝的基本特征以及成因机制的研究仍然很薄弱。

    目前对于地裂缝成因机理的研究已经形成了一套较为完整的理论,主要有以下三种观点,即构造成因、地下资源开采成因、构造和地下资源开采复合成因[19-24]。彭建兵等[25]通过对西安地裂缝的系统研究,将地裂缝的成因机制进一步细化为: 深部构造孕育地裂缝、盆地伸展萌生地裂缝、黄土介质响应地裂缝、断层活动伴生地裂缝、应力作用群发地裂缝、抽水作用加剧地裂缝、表水渗透开启地裂缝。在形成地裂缝的众多影响因素中,构造因素往往是地裂缝形成的内因,其控制了地裂缝的活动方式,其次不容忽视的是水的作用,其加剧了地裂缝的活动速率,缩短了其出露地表的进程[26]。地裂缝与表水入渗的关系非常密切,已有学者针对该问题进行了探索研究。LU等[27]对陕西三原县双槐树地裂缝的成因进行了研究,并对地裂缝带土的渗透特性进行了大型原位浸水试验。研究认为该地裂缝与隐伏断裂相连,其对地裂缝的形成具有控制作用,强降雨的渗透侵蚀作用使得地裂缝出露地表。乔建伟等[28]将临汾盆地果场地裂缝的成因机理概化为构造孕缝、抽水诱缝和降雨扩缝。AYALEW等[29]对埃塞俄比亚裂谷的Muleti小镇降雨后出现的地裂缝进行了研究,作者认为该地裂缝与含水层的压缩和水平渗流应力没有确切的关系,地裂缝的产生与大雨期间的管涌和塌陷过程有关。NGECU[30]研究了位于Menengai 火山与Nakuru湖之间的小镇出现的带状地面塌陷。作者认为在水的渗透侵蚀作用下松散堆积层形成了连续的地下空洞,最终由于车辆或火车交通引起的振动或者暴雨的情况下,丧失稳定性从而导致地下坍塌,最终到达地表形成塌陷。

    本文以肯尼亚裂谷段地裂缝为研究对象,通过资料收集、野外调查、槽探、钻探和室内土工试验等研究手段,对该区域地裂缝的基本特征以及成因机制进行了研究,研究成果可为该区域的工程建设提供科学依据,并为裂谷火山区地裂缝的研究提供借鉴。

    研究区位于东非大裂谷东支,肯尼亚裂谷的中部,属于东非高原区,地面高程1 650~1 820m,地跨东经36°17′30″~36°35′45″,南纬0°56′~1°09′(图1)。谷底区为两侧隆起区相夹而形成的“地堑”地貌。裂谷底部地势平坦开阔,植被覆盖率高,零星散布着多个火山,其中规模较大的为研究区北侧的Lognot火山和南侧的Suswa火山。Longonot火山为复式火山,火山口现已被茂密森林所覆盖,是一座休眠火山。Suswa火山为一双层盾形火山,里层为一圈深堑,外层形成大大小小的三十余个火山洞穴,属于休眠火山。研究区冲蚀作用强烈,分布有多条深切的冲沟,沟内无常年流水,在降雨后会形成短时汇流。在研究区发现多处地热气体逸出点,说明研究区地球内动力地质作用仍然很活跃。

    图  1  研究区地理位置
    Figure  1.  Geographical location of the study area

    调查区内活动断裂较少,主要分布于裂谷的肩部。区域上晚更新世-全新世断裂发育,沿东西隆起区和裂谷区分界线正断裂发育(图2)。在研究区南部裂谷区顺着裂谷延伸的方向也发育着众多全新世断裂。东部隆起区与裂谷中央区的边界发育一系列近于平行的正断层,台阶式断陷,整体断裂带影响宽度约15 km。北部断裂构造简单,为典型的地堑式正断层,走向NNW,由平行于裂谷方向的若干断裂组成,1928年1月6日曾发生里氏7.0级地震;西部隆起区断裂较少,走向基本和裂谷走向一致,主要由活动性较弱的晚更新世断裂组成,表现为断续延伸的线性。研究区内的活动断裂发育相对较少,因而地裂缝的形成与活动断裂相关性较小。

    图  2  研究区活动断裂分布图(据中国地震局修改)
    Figure  2.  Distribution map of active faults in the study area (modified according to China Earthquake Administration)

    研究区内主要发育新生界地层,特别是第四纪以来的地层,以第四系火山碎屑沉积及冲、洪积层为主,第四系和新近系火山岩亦分布较广。第四系覆盖层厚度分布不均,具体表现为谷底区及两翼沟谷内分布较厚,其他地段则普遍较薄。新近系岩层则主要分布于裂谷及两翼火山熔岩流区。通过研究区布设多个勘探钻孔,揭露出研究区具有特殊的“软硬软”地层结构(图3)。图中①−④为上部软层,主要以粉砂为主,为岩层上部的覆盖层,厚度5−15 m。中部硬层为厚层凝灰岩,厚度约40 m,上部和下部风化程度较高,中部为弱风化,且该层存在较多的岩体裂隙。凝灰岩下部又是一层粉土、粉砂为主的软质岩层。

    图  3  地层柱状图
    Figure  3.  Stratigraphic histogram

    上部覆盖层土体主要为粉砂,其孔隙比较大,密度很小,容易受到水力侵蚀。通过颗粒分析实验发现,土中含有较多的火山浮石颗粒。浮石的容重小,是一种多孔、轻质的玻璃质酸性火山喷出岩,可在水中浮起,且本身强度不高,风化后可形成不同大小的浮石颗粒,吸水能力很强。浮石土在工程上是一种不良的土质材料,其具低密度、多孔隙、结构较松散、弱胶结、吸水率高、承载力较低、难以碾密、抗冲蚀能力差的特性,容易引起地基不均匀沉降,在雨水作用下容易诱发滑坡等不良地质灾害[31]

    研究区浅表部覆盖层较薄且抗水力侵蚀能力较差,其下部岩层中存在较多的贯通裂隙,因而这种地层结构在长期的水力侵蚀作用下易于发生潜蚀破坏。

    研究区的岩层埋深较浅,存在多个天然岩体露头,其主要位于河道、采石场以及邻近裂谷肩部区域。通过岩体出露区节理的统计以及多个勘探钻孔的资料分析,获得研究区岩体破裂的基本发育特征。

    研究区内岩体裂隙较为发育,主要表现为平行以及X型展布的特征。在多个河道中观测到的岩体裂隙大多呈平行展布,节理面近垂直,部分节理尾部出现分叉,形态上弯曲粗糙呈现树枝状,表现为张节理的特征。在地裂缝带内岩体裂隙尤为发育,其发育的优势方向约为320°,与裂谷展布方向一致(图4)。在邻近裂谷肩部区域的天然岩石露头处,观测到岩体裂隙表现为X型共轭剪节理的特征,节理平直且相互切割形成网状形态(图5),两组优势节理走向分别为110°和203°。

    图  4  平行岩体裂隙
    Figure  4.  Parallel rock stratum fissures
    图  5  X型岩体裂隙
    Figure  5.  X-type rock stratum fissure

    通过钻探揭露发现钻探岩芯中多次出现竖向、斜向的裂缝,多为张裂缝,部分裂缝断面含有铁锈,个别张裂缝中有泥质填充。在勘探孔Z-73-120钻进过程中,发现岩芯出现裂缝。该勘探孔深度62.2 m,裂缝出现在39~44.7 m段,呈不连续的竖向裂隙将岩芯剖分为两半,其中39~40.0 m裂缝宽度约10 mm,竖向长度约1 m,无充填(图6);40.8 m~44.7 m裂缝宽度约为25 mm,竖向长度3.9 m,裂缝为黏土充填,黏土手摸有细腻滑感。两处裂缝岩芯断面处均未发现擦痕,应为张裂缝。

    图  6  钻孔Z-73-120深度40~45 m岩芯(岩芯被张裂缝一分为二,泥质充填)
    Figure  6.  Core of borehole z-73-120 with depth of 40−45 m (the core is divided into two parts by tensile fracture and filled with mud)

    研究区的岩体受张性和剪性岩体裂隙切割,岩体较为破碎。浅部岩体裂缝与地表分布的第四系覆盖层贯通性较好,是有利的地表水、地下水下渗、排泄通道。埋深较深的岩体裂缝与地表贯通性较好,特别是裂缝中含有黏土夹层的裂缝。浅表黏土颗粒,随地表水下渗至岩体裂缝处滞留,形成泥质夹层。此类裂缝与地表贯通性好,形成良好的地表水、地下水渗透、潜蚀、排泄通道,易于在上覆土层内形成自下而上的地裂缝。

    研究区位于东非高原,以热带草原气候区为主,每年的3—6月和10—12月是雨季,其余月份是旱季,其中3—6月是大雨季、10—12月是小雨季。年平均降水量600~1400 mm,年均蒸发量为1450~2200 mm[32]

    研究区西北约8.5 km处的纳库鲁气象站是距研究区最近的一个气象站。对其1980年到2019年40年的日降雨量进行统计分析,在1981年4月26日出现了最大日降雨量329.95 mm,该数据为40年来最大日降雨量,且在2007年到2013年,极端降雨出现的频率有所增加(图7)。

    图  7  肯尼亚裂谷日均降雨量(Nakuru台站,1980—2019年)
    Figure  7.  Daily average rainfall in the Rift Valley of Kenya (Nakuru Station, 1980—2019)

    通过野外调查发现,研究区地裂缝出露地表多发生于雨季,地裂缝的产生与降雨的相关性较高。研究区降雨具有短时降雨量大、频率高且极端降雨天气频发的特点。强烈的降雨可为地裂缝的形成提供有利的水力条件。

    研究区地裂缝研究基础较为薄弱,没有地裂缝活动的历史记录,多数地裂缝的形成历史已无从考证。通过对研究区多期Google卫星影像的解译以及地裂缝现场地质调查工作,研究现阶段地裂缝的平面分布规律,通过对典型地裂缝探槽揭露,研究该区域地裂缝的剖面结构特征。

    研究区共发育有21条地裂缝,大多分布在地势平坦开阔的裂谷中心区两个火山之间,靠近谷肩区域的地裂缝较少(图8)。

    图  8  地裂缝平面分布图
    Figure  8.  Plane distribution map of ground fissures

    研究区地裂缝大多呈平行分布且直线延伸,地表可见延伸长度174~5000 m(表1)。对地裂缝的走向进行统计,绘制地裂缝走向玫瑰图(图9)。该区域地裂缝的走向集中于34°~360°和20°~30°这两个区间,具有明显的定向性。肯尼亚裂谷在该段的走向约为353°,Logonot火山与Suswa火山连线走向约为25°,分别与地裂缝的优势走向重合。由此可见,该区域地裂缝的形成与裂谷水平拉张以及火山活动密切相关,且裂谷拉张作用为主导因素。

    表  1  地裂缝分布状况
    Table  1.  Distribution of ground fissures
    地裂缝编号走向/(°)长度/m主要特征
    DL0122758宽8 m,深4.5~5 m,平面呈直线状,底部可见明显水流冲刷痕迹,地层水平,未见明显位错
    DL023552000宽5 m,深6.6 m,地表为直线状出露
    DL033571900宽6 m,深4.8 m,地表为直线状出露
    DL04175306宽1.6~5 m,深2.7 m,地表为直线状出露
    DL05288393最宽处达12 m,深1.2 m,地表为直线状出露
    DL061811200宽8.4 m,深6.8 m,与B3公路相交,地表为直线状出露
    DL0763950宽10~15 m,深5.4 m,线性延伸较好,侧壁垂直,冲刷强烈
    DL08327665直线延伸,在平面上与DL007在平面垂直相交
    DL0935375宽4~6 m,平面形态为直线状出露。
    DL103405000出现时间为2018年3月13日傍晚,最大深度可达15 m,宽度约2~5 m
    DL11255947宽2~4 m,侧壁垂直,地表为直线状出露
    DL12340420宽10~20 m,底部植被覆盖茂盛,地裂缝北端可观测到岩层出露
    DL13239457裂缝最北端存在一近圆形陷坑,宽3.7 m,深0.3~0.6 cm
    DL1430190宽6 m,底部植被覆盖茂盛,地表为直线状出露
    DL15238174沿南东方向延伸,最终消失在田地里,两侧发育有落水洞,最深处1.5 m,植被发育茂盛
    DL16201400雁列状展布,南段宽2~4 m,北段宽2~3 m,最深处1.5 m
    DL17340500直线延伸,最深处1.5 m
    DL18330459直线延伸,深度0.2 m左右,裂缝侧面发育一排小落水洞
    DL193431256冲沟侧壁直立,现状宽2~5 m,深2~4 m
    DL203401000冲沟侧壁直立,现状宽4~5 m,深3~4 m
    DL213501594现状沟宽3~10 m,深2~3 m,直线延伸
    下载: 导出CSV 
    | 显示表格
    图  9  地裂缝走向玫瑰图(蓝线为裂谷走向,红线为火山连线方向)
    Figure  9.  Rose flower diagram of ground fissure strike (blue line is rift strike, red line is the direction of volcano connection)

    地裂缝平面形态为直线型,多表现为带状地面塌陷,部分为串珠状落水洞。平面组合形式为雁列状、平行状分布。直线型地裂缝的形成多受隐伏断裂所控制,而该区域现代构造活动强烈,岩层埋深较浅。根据野外岩体露头观测以及工程钻探资料,该区域岩体破裂较为发育,因此隐伏岩体破裂对该区域地裂缝的形成具有控制作用。地裂缝的平面组合特征,反映了区域地应力的特征。平行分布的地裂缝常为拉张应力环境的产物[25]。肯尼亚裂谷为主动性裂谷,地下热点的岩浆上涌,到岩石圈后被阻挡,向四周水平流动,牵引岩石圈水平拉张破裂,断陷而形成裂谷[33]。因此,该区域水平拉张应力为该区域的主应力状态。雁列式分布的地裂缝常形成于剪切应力状态。该区域分布有多处火山,因火山活动的影响,在某一时期,该区域存在局部剪切应力,从而形成雁列式的岩体破裂形态,进而可形成雁列式的地裂缝。

    该区域地裂缝平面形态主要为直线型,表现为构造地裂缝的特征,受隐伏岩体破裂控制。该区域地裂缝走向的定向性以及平面组合形式显示裂谷的水平拉张作用和火山活动与地裂缝的形成密切相关。

    为研究该区域地裂缝的剖面结构特征,对铁路北侧约2.2 km的地裂缝进行了探槽揭露。通过Google Earth时间轴观测,该地裂缝形成于2010年至2013年期间,该时间段也为该区域极端降雨天气频发的阶段(图7)。该地裂缝现状表现为地表陷坑及冲沟,沟宽2~7 m,深1~2 m,从北向南逐渐变浅,沟两侧侵蚀减弱。全长398 m,走向为197°,中部断开,呈断续出露。探槽开挖位置位于该地裂缝的北端(图10 a和图10 b)。

    图  10  探槽剖面图
    Figure  10.  Sectional view of the trench

    探槽长约22 m,宽12 m,深6.5 m,垂直地裂缝走向开挖。探槽底部为弱风化的火山角砾岩(图10 f),施工机械开挖至该层后已很难向下挖掘,该层上部为砂土和浮石的互层。火山角砾岩存在两条岩层裂缝,裂缝向上延伸直至地表形成两条地裂缝,主裂缝f1和次级裂缝f2。裂缝两侧地层连续性较好,地层水平,无垂直位错,两条裂缝近似呈平行分布。

    f1裂缝剖面形态呈上大下小的楔形体,岩层面以上部分以及岩体裂缝均已被粉砂填充(图10 e)。裂缝顶部最宽处约2.8 m,底部最窄处0.7 m,对应于地表冲沟的位置。裂缝长期受表水入渗的影响,裂缝处土呈黑棕色,有机质含量较高。f2裂缝宽度较小约5~10 cm,走向215°,岩层面以上部分已被填充,岩体裂缝仍处于张开状态,未被填充,且在开挖后有凉气逸出,持续半天后,气体逐渐消失(图10 g和图10 h)。裂缝在剖面呈近直立形态,在地表未形成明显的地表裂缝,表现为串珠状落水洞。

    探槽揭露出的两条地裂缝,下部均与岩体裂缝相连接,地层无垂直位错,裂缝处土体水力侵蚀现象明显,剖面形态上呈现上大下小的楔形体。岩体裂缝的存在是地裂缝形成的前提条件,强烈的水力侵蚀以及松散易侵蚀的土体是裂缝出露地表的决定性因素。

    研究区地裂缝的形成与区域构造应力、降雨、浅表部土体性质以及地层结构均密切关系。

    研究区位于主动型裂谷的谷底区,区域构造应力主要为拉张应力。该区域散布有多个大大小小的火山,在火山活动作用的影响下,存在局部的剪应力区域。在水平拉张以及火山活动作用下,岩体形成不同宽度以及延伸长度的岩体裂缝,在平面上表现为直线型平行和雁列状分布的地裂缝。根据研究区地裂缝走向统计分析,研究区地裂缝的走向与裂谷走向以及两个主要火山的走向具有明显的相关性,且与裂谷走向的相关性更高。

    研究区降雨具有短时降雨量大,频率高、极端降雨天气频发的特点,从而使得研究区水力侵蚀特别强烈。研究区多条地裂缝均是在一场大雨后出露地表的,与降雨存在较高相关性。从多期谷歌影像解译出的地裂缝,其出现时间也多位于2007年到2013年极端降雨频发的时间段。因而该区地裂缝的形成与降雨入渗存在密切关系。当出现强降雨并导致地面大量积水时,由于岩体裂缝形成的过程中破坏了土层结构,加强了上、下土层之间的连通,增强了土体的渗透性,使得雨水较容易沿地裂缝带下渗到土层深部,导致裂缝带土体强度降低并发生渗透或侵蚀变形,同时地表积水还会增大裂缝带土体的上部荷载[24]。在表水下渗的影响下,可能会触发地裂缝重新开裂和活动。

    该区域浅表部土体主要为粉砂以及浮石的互层,土质疏松,孔隙比较大,密度较小,且其在浸水饱和后,抗剪强度衰减幅度很大,因而其抗水力侵蚀能力较差。土体中含有较多的浮石颗粒,具有较强的吸水能力,从而使得土体浸水饱和后,其饱和重度相对于天然重度变化较大,在重力作用下更容易发生坍塌。在岩体裂缝形成的过程中,浅表部土体也被一定程度的扰动,在裂缝上部可形成一个破碎带。表水沿着该破碎带优势入渗,岩体裂缝上部逐渐在渗流的作用下由下而上发生垮塌,最终出露地表形成地裂缝。

    该区域典型的隐伏岩体破裂以及 “软硬软”地层结构为水流裹挟上部软层物质提供了堆积场所。岩体的构造裂缝是上部松散沉积物在水力侵蚀作用下的堆积场所,裂缝的宽度、破裂岩层的厚度以及裂缝的连通性,决定了其储存物质的能力的大小。裂缝越宽,破裂岩层越厚,连通性越好,其储存物质的能力也就越强。当充填物质填满岩体裂缝时,浅表部地裂缝的发展也就结束了。该区域钻孔揭露的岩层厚度为40~50 m左右,上部软层厚度较小5~15 m左右。岩层裂缝处具有足够的储存空间可以容纳上部软层的水力侵蚀物质。软硬岩交界处也是水力侵蚀最强的位置,因而与中部硬层相接的下部软层的水力侵蚀也会比较强烈,可形成潜蚀空腔,扩大了容纳物质的能力。这种地层结构的形成与火山多期次喷发有关,在深部可能存在多套这种结构地层,从而使得上部土体在渗流作用可下进入深部地层。

    综上,裂谷拉张和火山活动等内动力地质作用是该区域地裂缝形成的控制因素;浅表部松散易潜蚀的土体是地裂缝形成的物质基础;强降雨为地裂缝的形成提供了水力条件,是地裂缝出露地表的动力源;具有岩体裂缝的“软硬软”的特殊地层结构为水力侵蚀物质提供了有利的堆积场所。构造作用产生了早期的岩体破裂,其控制了该区域地裂缝的基本格局,是该区域地裂缝形成的主控因素。强降雨是后期地裂缝出露地表的一个重要诱发因素和动力源。

    可将裂谷区地裂缝的形成过程分为以下3个阶段,如图11所示。

    孕育阶段:在东非大裂谷区域拉张应力、深部热运动以及火山的周期性活动的顶托作用下导致松散层下部岩层产生拉裂破坏,形成隐伏岩体破裂。岩体破裂形成的过程中,对裂缝上部一定范围土体产生扰动,使得其土体强度降低,孔隙度增大,易于被水流冲刷。在自重下作用下,土体产生一定程度塌落,最终形成初始状态的塌落空腔,空腔顶部形成塌落拱。

    扩展阶段:降雨入渗作用下,在土体内部形成渗流。塌落空腔顶部与侧面土体发生潜蚀作用,水裹挟土体颗粒通过隐伏岩体裂缝进入岩层下部土体。随着潜蚀作用的持续进行,塌落拱逐渐向上延伸,上覆土层厚度逐渐减小,而其跨度逐渐增大,当跨度刚好等于极限跨度时,则上覆土体处于极限平衡状态。这一阶段是极其漫长的过程,地表并没有明显的破坏迹象。

    成灾阶段:地表则逐渐出现微弱沉降变形或是细小裂缝,一般不容易被观测到。当遇到一场大雨时,上覆土体容重增大的同时强度降低,且作用有渗流的拖曳力,当达到土体极限破坏强度时,上覆土体发生坍塌或开裂,使得地裂缝出露地表。破坏过程一般较快且具有突发性的特点,很难提前监测到。

    图  11  肯尼亚裂谷地裂缝成因机理
    Figure  11.  The genetic mechanism of ground fissures in the Kenya Rift Valley

    (1)肯尼亚裂谷区地裂缝是典型的构造型地裂缝,地裂缝裂面陡直,以水平拉裂和垂直塌陷为主;

    (2)裂谷区地裂缝的走向与裂谷走向以及区域内两个主要火山的连线方向具有明显的相关性,在剖面上表现为上大下小的楔形体,填充物水力侵蚀现象明显,裂缝两侧无垂直位错;

    (3)裂谷拉张和火山活动等内动力地质作用是该区域地裂缝形成的控制因素;浅表部松散易潜蚀的土体是地裂缝形成的物质基础;强降雨为地裂缝的形成提供了水力条件,是地裂缝出露地表的动力源;“软硬软”的地层结构为水力侵蚀物质提供了有利的堆积场所;

    (4)裂谷区地裂缝的形成可分为3个阶段:孕育阶段、扩展阶段和成灾阶段。

  • 图  1   边坡设计断面

    Figure  1.   Design profile of slope

    图  2   浅层滑坡情形

    Figure  2.   View of a shallow landslide

    图  3   有限元几何模型

    Figure  3.   Geometric model of finite element

    图  4   有限元网格划分

    Figure  4.   Finite element meshes

    图  5   一般工况塑性应变等值云图

    Figure  5.   Plastic strain contours under general working condition

    图  6   一般工况水平位移等值云图

    Figure  6.   Horizontal displacement contours under general working condition

    图  7   地下水渗流工况孔隙水压力等值云图

    Figure  7.   Pore water pressure contours under groundwater seepage

    图  8   地下水渗流工况塑性应变等值云图

    Figure  8.   Plastic strain contours under groundwater seepage

    图  9   长时小雨工况塑性应变等值云图

    Figure  9.   Plastic strain contours under long-term light rainfall

    图  10   长时小雨工况水平位移等值云图

    Figure  10.   Horizontal displacement contours under long-term light rainfall

    图  11   工况3孔隙水压力等值云图

    Figure  11.   Pore water pressure contours under Working Condition 3

    图  12   工况3塑性应变等值云图

    Figure  12.   Plastic strain contours under Working Condition 3

    图  13   工况3水平位移等值云图

    Figure  13.   Horizontal displacement contours under Working Condition 3

    图  14   方案1塑性应变等值云图(1∶1.20)

    Figure  14.   Plastic strain contours of Scheme 1(1∶1.20)

    图  15   方案1塑性应变等值云图(1∶1.25)

    Figure  15.   Plastic strain contours of Scheme 1(1∶1.25)

    图  16   方案2塑性应变等值云图(1∶1.20)

    Figure  16.   Plastic strain contours of Scheme 2 (1∶1.20)

    图  17   方案2塑性应变等值云图(1∶1.25)

    Figure  17.   Plastic strain contours of Scheme 2 (1∶1.25)

    表  1   土工参数

    Table  1   Geotechnical properties of soils

    土层编号γ/(kN·m−3)γd/(kN·m−3)E/kPaυc/kPaφ/(°)c'/kPaφ/(°)
    7-21019.515.680000.2824.916.822.420.6
    7-21120.016.510 0000.3028.917.824.622.2
    下载: 导出CSV

    表  2   不同降雨工况边坡稳定安全系数

    Table  2   Safety factors of slope under different rainfall conditions

    工况降雨强度/(mm·h−1降雨历时/h安全系数
    12041.152
    23041.113
    34041.091
    44031.106
    54021.140
    下载: 导出CSV
  • [1]

    ZIENKIEWICZ O C, HUMPHESON C, LEWIS R W. Associated and non-associated visco-plasticity and plasticity in soil mechanics[J]. Geotechnique,1975,25(4):671 − 689. DOI: 10.1680/geot.1975.25.4.671

    [2] 赵尚毅, 郑颖人, 时卫民, 等. 用有限元强度折减法求边坡稳定安全系数[J]. 岩土工程学报,2002,24(3):343 − 346. [ZHAO Shangyi, ZHENG Yingren, SHI Weimin, et al. Analysis on safety factor of slope by strength reduction FEM[J]. Chinese Journal of Geotechnical Engineering,2002,24(3):343 − 346. (in Chinese with English abstract) DOI: 10.3321/j.issn:1000-4548.2002.03.017
    [3] 郑颖人, 赵尚毅, 宋雅坤. 有限元强度折减法研究进展[J]. 后勤工程学院学报,2005,21(3):1 − 6. [ZHENG Yingren, ZHAO Shangyi, SONG Yakun. Advance of study on the strength reduction finite element method[J]. Journal of Logistical Engineering University,2005,21(3):1 − 6. (in Chinese with English abstract) DOI: 10.3969/j.issn.1672-7843.2005.03.001
    [4] 李湛, 栾茂田, 刘占阁, 等. 渗流作用下边坡稳定性分析的强度折减弹塑性有限元法[J]. 水利学报,2006,37(5):554 − 559. [LI Zhan, LUAN Maotian, LIU Zhange, et al. Elasto-plastic FEM for embankment stability analysis under the condition of seepage based on the technique of shear strength reduction[J]. Journal of Hydraulic Engineering,2006,37(5):554 − 559. (in Chinese with English abstract) DOI: 10.3321/j.issn:0559-9350.2006.05.006
    [5] 张昊, 余巍伟, 林杭, 等. 不同安全系数对应的边坡滑动面位置分析[J]. 岩土力学,2012,33(2):449 − 452. [ZHANG Hao, YU Weiwei, LIN Hang, et al. Analysis of location distribution of slope slip plane with different factors of safety[J]. Rock and Soil Mechanics,2012,33(2):449 − 452. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-7598.2012.02.021
    [6] 刘明维, 郑颖人. 基于有限元强度折减法确定滑坡多滑动面方法[J]. 岩石力学与工程学报,2006,25(8):1544 − 1549. [LIU Mingwei, ZHENG Yingren. Determination methods of multi-slip surfaces landslide based on strength reduction fem[J]. Chinese Journal of Rock Mechanics and Engineering,2006,25(8):1544 − 1549. (in Chinese with English abstract) DOI: 10.3321/j.issn:1000-6915.2006.08.005
    [7] 毕港, 陈征宙, 姜玉平, 等. 一种边坡多滑动面搜索的新方法[J]. 水文地质工程地质,2011,38(6):24 − 28. [BI Gang, CHEN Zhengzhou, JIANG Yuping, et al. A new method for searching all possible circular slips of a slope[J]. Hydrogeology & Engineering Geology,2011,38(6):24 − 28. (in Chinese with English abstract)
    [8] 李小春, 袁维, 白冰, 等. 基于局部强度折减法的边坡多滑面分析方法及应用研究[J]. 岩土力学,2014(3):847 − 854. [LI Xiaochun, YUAN Wei, BAI Bing, et al. Analytic approach of slope multi-slip surfaces based on local strength reduction method and its application[J]. Rock and Soil Mechanics,2014(3):847 − 854. (in Chinese with English abstract)
    [9]

    LIN H, CAO P, GONG F Q, et al. Directly searching method for slip plane and its influential factors based on critical state of slope[J]. Journal of Central South University of Technology,2009,16(1):131 − 135. DOI: 10.1007/s11771-009-0022-6

    [10] FREDLUND D G, RAHARDJO H. 非饱和土力学[M]. 陈仲颐, 张在明, 陈愈炯, 等. 译. 北京: 中国建筑工业出版社, 1997.

    FREDLUND D G, RAHARDJO H. Soil mechanics for unsaturated soils[M]. CHEN Zhongyi, ZHANG Zaiming, CHEN Yujiong, et al. translator. Beijing: China Building Industry, 1997. (in Chinese with English abstract)

    [11]

    ALONSO E, GENS A, LLORET A, et al. Effect of rain infiltration on the stability of slopes[J]. Unsaturated Soils J.,1995,1:241 − 249.

    [12]

    CHO S E, LEE S R. Instability of unsaturated soil slopes due to infiltration[J]. Computers and Geosciences,2001,28(3):185 − 208.

    [13] 黄盛锋, 陈志波, 郑道哲. 基于灰色关联度法和强度折减法的边坡稳定性影响因素敏感性分析[J]. 中国地质灾害与防治学报,2020,31(3):35 − 40. [HUANG Shengfeng, CHEN Zhibo, ZHENG Daozhe. Sensitivity analysis of factors influencing slope stability based on grey correlation and strength reduction method[J]. The Chinese Journal of Geological Hazard and Control,2020,31(3):35 − 40. (in Chinese with English abstract)
    [14] 陈子玉, 宋彦辉, 严豪, 等. 双参数强度折减法研究中存在的问题分析[J]. 水文地质工程地质,2019,46(2):125 − 132. [CHEN Ziyu, SONG Yanhui, YAN Hao, et al. Analyses of the existing problems in the double parameters reduction method[J]. Hydrogeology & Engineering Geology,2019,46(2):125 − 132. (in Chinese with English abstract)
    [15] 中华人民共和国交通运输部. 公路路基设计规范: JTG D 30—2015[S]. 北京: 人民交通出版社, 2015.

    Specifications for design of highway subgrades: JTG D 30—2015[S]. Beijing: China Communications Press, 2015. (in Chinese)

    [16] 陈丽刚. 基于ABAQUS渗流与应力耦合作用的边坡稳定性分析[D]. 郑州: 郑州大学, 2010.

    CHEN Ligang. The slope stability analysis of taking into account the seepage and stress coupling effect based on ABAQUS[D]. Zhengzhou: Zhengzhou University, 2010. (in Chinese with English abstract)

    [17] 费康, 彭劼. ABAQUS岩土工程实例详解[M]. 北京: 人民邮电出版社, 2017.

    FEI Kang, PENG Jie. Detailed explanation of geotechnical engineering cases using ABAQUS[M]. Beijing: Posts & Telecom Press, 2017. (in Chinese)

    [18] 徐翔, 王义兴, 方正. 降雨和库水联合作用下边坡稳定性变化规律[J]. 中国地质灾害与防治学报,2019,30(2):51 − 58. [XU Xiang, WANG Yixing, FANG Zheng. Analysis on stability of a slope under rainfall and reservoir water level variation[J]. The Chinese Journal of Geological Hazard and Control,2019,30(2):51 − 58. (in Chinese with English abstract)
  • 期刊类型引用(4)

    1. 刘阳,王飞永,贺鸣,赵玺. 表水侵蚀地裂缝物理模型试验研究. 工程地质学报. 2024(01): 183-193 . 百度学术
    2. 刘振华,何付兵,崔玉斌,刘晓勇,王凯,张悦泽,赵艳英. 平谷地裂缝基本特征与成因分析. 科学技术与工程. 2023(04): 1436-1446 . 百度学术
    3. 孟振江,彭建兵,李超,乔建伟,王飞永,康尘云,赵俊彦,张凡. 耦合型地裂缝活动特征与成因机制模拟研究——以北京宋庄地裂缝为例. 水文地质工程地质. 2023(03): 138-148 . 百度学术
    4. 王晗,邓亚虹,慕焕东,薛捷. 西安典型地裂缝场地地脉动测试及地震响应特征分析. 中国地质灾害与防治学报. 2022(04): 55-64 . 本站查看

    其他类型引用(1)

图(17)  /  表(2)
计量
  • 文章访问数:  501
  • HTML全文浏览量:  321
  • PDF下载量:  954
  • 被引次数: 5
出版历程
  • 收稿日期:  2020-03-15
  • 修回日期:  2020-05-18
  • 网络出版日期:  2021-04-26
  • 刊出日期:  2021-04-26

目录

/

返回文章
返回