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

慈竹对小型落石的拦截效果分析

陈海东, 叶四桥, 曾彬, 张立舟, 崔开林, 刘蓓, 辛立平

陈海东,叶四桥,曾彬,等. 慈竹对小型落石的拦截效果分析[J]. 中国地质灾害与防治学报,2023,34(5): 97-106. DOI: 10.16031/j.cnki.issn.1003-8035.202206024
引用本文: 陈海东,叶四桥,曾彬,等. 慈竹对小型落石的拦截效果分析[J]. 中国地质灾害与防治学报,2023,34(5): 97-106. DOI: 10.16031/j.cnki.issn.1003-8035.202206024
CHEN Haidong,YE Siqiao,ZENG Bin,et al. Analysis on the interception effect of neosinocalamus affinis on small falling rocks[J]. The Chinese Journal of Geological Hazard and Control,2023,34(5): 97-106. DOI: 10.16031/j.cnki.issn.1003-8035.202206024
Citation: CHEN Haidong,YE Siqiao,ZENG Bin,et al. Analysis on the interception effect of neosinocalamus affinis on small falling rocks[J]. The Chinese Journal of Geological Hazard and Control,2023,34(5): 97-106. DOI: 10.16031/j.cnki.issn.1003-8035.202206024

慈竹对小型落石的拦截效果分析

基金项目: 重庆市研究生导师团队建设项目(JDDSTD2022009)
详细信息
    作者简介:

    陈海东(1997-),男,重庆石柱人,硕士研究生,主要从事岩土工程防灾减灾方面研究。E-mail:1217430861@qq.com

    通讯作者:

    叶四桥(1978-),男,湖北孝感人,博士(后)、教授,博士生导师,主要从事岩土工程防灾减灾方面研究和教学工作。E-mail:yesiqiao@cqjtu.edu.cn

  • 中图分类号: P642.21

Analysis on the interception effect of neosinocalamus affinis on small falling rocks

  • 摘要: 落石灾害是我国西南地区常见的一种小型崩塌。落石对下方公路、桥梁、铁路等基础设施安全造成严重威胁。因此,如何对落石进行有效的拦截就显得尤为重要。常见的拦截设施,如拦石沟,拦石网,挡墙等,在布设时极易受地形、施工等因素的限制,有时难以大面积地布置,导致拦挡效果大为减弱。我国的西南地区生长着大量慈竹,慈竹具有生长快、繁殖能力强、抗弯折等特点,对落石具有很好的拦截效果。文章通过模型试验与数值模拟相结合的方法,对慈竹的抗冲击能力进行了研究并提出了慈竹拦截小型落石的方案。结果如下:(1)慈竹破坏所需能量随着直径和壁厚的变大而上升,随着年龄、长重比和长径比的增大而降低;(2) 一丛慈竹在抗冲击过程中最少能消耗3975.55 J能量,最多能消耗10890.88 J能量;(3)结合CRockfall软件进行慈竹落石拦截效应计算,结果显示拦截边坡高度43 m、直径0.5 m的危岩,需栽种2列3行6丛慈竹。
    Abstract: Dangerous rock disaster is prevalent geological hazard in southwest China. Under the influence of gravity, dangerous rock masses can easily trigger rockfalls, posing a significant threat to infrastructure such as roads, bridges, and railways. Therefore, it is crucial to effectively intercept rockfalls. Common interception facilities, including stone-blocking trenches, stone-blocking nets, and retaining walls, are often constrained by terrain and construction limitations, making it difficult to implement them on a large scale and resulting in reduced interception effectiveness. The southwestern region of China is rich in Cizhu bamboo, which exhibits rapid growth, strong reproductive capabilities, and resistance to bending, rendering it highly effective in intercepting falling rocks. This study employs a combination of model testing and numerical simulation to investigate the impact resistance of Cizhu bamboo and proposes a strategy for using Cizhu bamboo to intercept small-scale rockfalls. The research findings are as follows: (1) The energy requirements of Cizhu bamboo increase with diameter and wall thickness, while decreasing with age, length-to-weight ratio, and aspect ratio. (2) A cluster of Cizhu bamboo can consume a minimum of 3975.55 J and a maximum of 10890.88 J of energy during impact resistance. (3) Utilizing the CRockfall software, the interception effect of Cizhu bamboo on rockfalls is calculated, indicating that a slope with a height of 43 m and a dangerous rock of 0.5 m in diameter requires the planting of 2 rows, 3 lines, and 6 clusters of Cizhu bamboo.
  • 流滑型黄土滑坡是黄土地区沿沟道或斜坡发生远程滑动和形成长条状堆积的特殊类型滑坡,一般具有明显的滑源区、流通区和堆积区。流滑型黄土滑坡多为高速远程滑坡,发生过程犹如“蛟龙出沟”,常造成意想不到的严重灾害。如1978年7月12日,天水市麦积区刘家湾沟脑因暴雨诱发体积约2.50×106 m3的黄土滑坡,前部滑体以平均6 m/s的速度沿沟道滑动1550 m,压埋沟口陇海铁路150 m,中断行车15 d[1]。2013年7月21日天水市马跑泉镇大沟发生体积约29.7×104 m3的流滑型滑坡,滑距达850 m,摧毁9户村民137间房屋及村委会[24]

    长期以来,高速远程滑坡运动机理一直是滑坡研究的热点,学者们提出了众多机理和假说。1932年HEIM[5]对瑞士Elm滑坡研究后提出了颗粒流理论,并建立了高速远程滑坡运动的“雪橇”模型;KENT[6]对美国Madison峡谷滑坡研究后提出了空气润滑理论。基于不排水环剪试验,SASSA等[7]提出了“滑动面液化”的高速滑坡机制;汪发武[8]也提出了土颗粒破碎导致超孔隙水压力的高速滑坡形成机制;HUTCHINSON等[9]提出了流滑性滑坡的不排水荷载效应;殷跃平[10]研究了汶川8级地震触发的高速远程滑坡滑动中的抛掷、碰撞、铲刮和气垫效应;刘传正[11]分析了滑坡规模、不同岩性和状态的滑道对高速远程滑坡的影响。胡广韬等[12]按照滑坡启动加速度、滑速两方面特征进行了组合分类和研究。段钊等[13]对陕西泾阳南塬流滑型黄土滑坡的滑动特征和液化机理进行了研究,统计出其滑距为坡高的4倍,属典型的高速远程滑坡。彭建兵等[2]、张帆宇等[3]、翟张辉等[4]对天水大沟滑坡泥石流的运动过程和速度进行了模拟计算。王玉峰等[14]将目前的研究成果总结为摩擦生热减阻、滑带液化减阻、动力破碎减阻、底部裹挟减阻、剪切振动减阻和动量传递减阻等六大滑动机理类型。

    2013年7月22日7时45分,甘肃发生岷县漳县Ms6.6级地震。在距震中4 km处的岷县梅川镇永光村触发了2处黄土滑坡,其中永光1#滑坡体积约23×104 m3,造成12人遇难。滑坡前后缘高差175 m,滑坡总长度1030 m,二者之比值为0.17,小于0.33,属远程滑坡。作为流滑型黄土滑坡研究的典型事例,许多学者开展了成因和滑动特征的研究[1517]。但在永光1#滑坡的滑动过程、历时长短和滑速快慢等方面有不同认识。本文通过现场调查、影像对比和滑动过程观察资料的综合分析,探讨了其滑动过程特征、不同部位的滑速及变化情况,并简要分析了滑动机理。以期为此类滑坡的进一步研究和数值模拟提供参考。

    岷县地处青藏高原东北边缘阶梯地形带,是甘南高原、陇西黄土高原和陇南山地的过渡地带,海拔2040~3872 m,切割深度500~1000 m。南部为迭山山系,北部为西秦岭山地,中部为相对低缓的洮河谷地。

    本区属西秦岭构造带的北支西伸部分,地质构造总体展布方向为NWW向(图1)。其中,临潭—岷县—宕昌活动断裂带全新世以来活动强烈,地震频发。历史上区内曾发生过3次Ms6.0级以上强震[18],2003年以来岷县发生Ms5.0级以上中强地震4次。地震动峰值加速度0.15 g

    图  1  “7•22”岷县漳县地震烈度与构造背景图
    Figure  1.  Seismic intensity and structural background map of the July 22 earthquake in Minxian-Zhangxian earthquake

    岷县地区高寒阴湿,年均降水量560.8 mm,年最大降水量709.3 mm,5—9月降水量占全年的78%以上,随海拔的升高降水量增加明显。日最大降水量94 mm,1 h最大降水量69.2 mm。

    永光1#滑坡位于岷县北部西秦岭南麓黄土覆盖的斜坡区,斜坡坡向南,平均坡度11.5°。斜坡上部为南北长700 m、东西宽约500 m的缓坡,冲沟不甚发育,坡面较为完整;斜坡下部冲沟较发育,切割深度20~150 m,地形较为破碎。斜坡主要由上更新统马兰黄土和古近系泥岩组成,马兰黄土一般厚度5~20 m。黄土底部地下水较丰富,含水层厚度2~3 m,隔水层为泥岩[17]。永光1#滑坡位于黄土斜坡中下部近南北向小型冲沟的沟脑部位,坡脚曾有泉水分布。

    永光1#滑坡位于岷县梅川镇永光村四社北部黄土斜坡的中下部。滑坡总面积4.2×104 m2,平均厚度5 m,总体积23×104 m3,滑体主要由马兰黄土组成,局部夹少量泥岩碎块,为中型黄土滑坡。“7•22”岷县漳县地震触发其剧烈滑动,大量堆积在下部平台村庄区,破坏8户村民房屋,造成12人死亡。前缘约6×104 m3滑体滑入前部近南北向较宽缓的小型冲沟,再转向进入主沟道远程流滑并沿沟底不断堆积,其前舌直达沟口乡村公路小桥处。滑坡总长度1030 m,滑坡后缘高程约2702 m,滑舌处高程约2527 m,前后缘高差175 m,滑坡前后缘高差与滑坡总长度之比(H/L)为0.17,小于0.33,属远程滑坡。

    永光1#滑坡平面形态复杂,总体呈不规则的“L型”(图2图3)。根据滑坡形态、堆积特征和滑动环境,将滑坡分为滑源区—平台堆积区和流通区—沟道堆积区两大部分(图4)。

    图  2  永光滑坡正射影像图
    Figure  2.  Orthophoto image map of Yongguang landslide
    图  3  永光滑坡平面图
    Figure  3.  Plan view of the Yongguang landslide
    图  4  永光1#滑坡工程地质剖面图
    Figure  4.  Engineering geological cross-section of the Yongguang 1# landslide

    滑源区—平台堆积区为滑坡上部(图5)。平面形态近似矩形,南北长290 m,东西宽75~110 m,面积约2.76×104 m2,平均厚度6 m,堆积体积约17×104 m3,主滑方向210°。

    图  5  永光1#滑坡滑源区—平台堆积区
    Figure  5.  Sliding source area of Yongguang 1# landslide:platform accumulation area

    滑坡后缘高程2702 m,剪出口高程约2645 m,高差57 m,平均坡度13°。上部斜坡较陡,坡度18°左右;下部较缓,坡度约6°,为南北宽度约180 m的平台,是永光四社村民院落区。地震首先触发上部黄土陡坡段失稳滑动,为永光1#滑坡的滑源区,主滑动面位于饱和黄土与泥岩接触面附近[17]。滑坡在前部较宽缓的平台区大量堆积。由图2图3可看出,滑动过程中,滑坡向东侧地形较低的村民院落区扩散和堆积,形成长约150 m、宽20~40 m、高3~5 m的鼓丘,展布方向与滑动方向近于一致,压埋了村庄。滑坡西侧中前部为长约130 m、宽50 m、深3~5 m,向南部冲沟敞开的洼地,洼地内横向拉张裂缝发育。鼓丘与洼地之间因差异滑动形成的剪切带明显。

    前缘中下部有部分坡体保留了滑坡前的原地形,宽约15 m、高约12 m的可见范围土体结构和地表植被完整,说明未发生滑动。其顶部有6~8 m厚的滑坡堆积,地面局部反翘。

    永光1#滑坡后壁平面形态呈圈椅状弧形,高10~15 m,坡度约55°,由马兰黄土组成,其上擦痕较为清晰。

    另外,永光1#滑坡西北部还发生一处相对独立的次级黄土滑坡。滑源区前部宽30 m,后部宽55 m,长60 m,面积2800 m2,体积约1.2×104 m3。宽20~30 m,长约190 m的滑体呈长条状叠加堆积在已滑的滑坡体上,滑动方向由160°向南偏转为210°。从滑坡叠加堆积次序分析,该次级黄土滑坡发生的时间稍晚。

    由上可见,永光1#滑坡的滑源区—平台堆积区滑坡堆积量较大,为该滑坡的主体部分,主滑面位于黄土、泥岩界面处,属地震诱发的黄土—泥岩接触面滑坡。此类滑坡在岷县地区发育较广泛。

    为永光1#滑坡前缘滑体沿沟道发生特殊流滑型远程滑动的部分,以南部正对的小型宽浅冲沟作为流通区,呈流滑状向南滑动,并转向东南方向进入狭窄的主沟道继续发生远程流滑,沿主沟道不断堆积形成长条状滑体(图6)。

    图  6  流通区—沟道堆积区及堰塞湖
    Figure  6.  The sliding body flowing along the channel and barrier lake

    流通区沟道地形较宽缓顺直,沟道坡度大。塑流状滑体沿流通区快速向下游滑移,在流通区滑坡堆积物较少,厚度2 m左右。

    滑体主要沿主沟道堆积,长590 m,宽度10~35 m,面积1.5×104 m2,平均厚度4 m,体积约6×104 m3,占永光1#滑坡总体积的1/4。长宽比达24,呈典型的长条状,平面形态及滑动方向主要受沟道控制。主沟道平均坡度6.5°,最陡段9.3°,最缓段5.5°,在流通区与主沟交汇处、沟道较缓地段滑坡堆积厚度较大,最厚处6 m左右。滑坡前舌位于主沟沟口小桥处,厚度2 m左右,表面散布破坏民房的木材,小桥受滑坡推挤而拱起、开裂(图7)。滑坡在主沟道滑动过程中,未出现因刮铲两岸沟坡形成的滑塌和弯道外侧的明显超高。

    图  7  滑坡前舌及挤压变形的小桥
    Figure  7.  Small bridge at the front lobe of the landslide with extrusion deformation

    同时,在流通区与主沟汇合处有部分滑体向上游滑动约40 m,堵塞沟道,形成小型堰塞湖。

    根据滑坡形态特征、滑坡前后遥感影像、地质环境条件和观察到的部分滑动过程资料分析认为:永光1#滑坡在滑源区—平台堆积区和流通区—沟道堆积区经历了2次加速—减速过程,滑动过程复杂多变、特征差异明显。

    在强烈地震、滑坡灾害叠加的情景下,永光1#滑坡最初的滑动情况难以知晓。根据树木、滑坡迹象调查对比分析,滑源区—平台堆积区的滑距多为50~100 m,主滑方向210°。滑坡后部及东侧滑距50 m左右,滑速较低,东侧形成鼓丘,压埋村民院落。西侧中前部滑动相对强烈,最大滑距约130 m,估计最大滑速3~4 m/s,部分滑体滑入前部沟道,形成向南敞开的长条状洼地。受强烈地震作用,永光1#滑坡失稳后快速滑动,但加速阶段持续时间短;在孔隙水压力作用下,后期减速阶段滑速较低而持续时间较长。

    在平台南部与流通区上部衔接处中部有未滑动坡体,顶部滑坡堆积较厚、地面反翘;前缘部分滑体从其东西两侧分两股滑入流通区。据此分析认为,滑坡前缘一带的滑速较低且部分已发生制动,部分滑体以很低速度持续滑入下部冲沟。否则,前缘中部临空条件好、强度低的未滑动土质坡体很难抵抗上部快速滑坡的推力而保留下来,也将有更多滑体来不及制动而滑入前部流通区。

    为滑坡南部小型冲沟,高差40 m,水平长度约150 m,沟道坡度14.5°,相对较陡。冲沟呈上部宽浅、下部窄深的漏斗状,方向195°。漏斗口正对滑坡前缘,成为永光1#滑坡前缘部分滑体进入冲沟继续滑动的通道,是该滑坡滑动的重要加速段。滑体分两股从未滑动坡体两侧进入流通区,以西侧通道为主,又开始加速滑动。到达沟底时部分滑体堆积,并在对岸斜坡区形成垂直高度约4 m的逆冲超高(图89)。

    图  8  流通区下段与主沟道处滑坡堆积
    Figure  8.  Landslide deposits in the lower section of the runout zone and the accumulation of sliding body in the main channel
    图  9  流通区滑速计算剖面简图
    Figure  9.  Simplified profile for calculating slide velocity in the runout zone

    永光1#滑坡前缘部分滑体从流通区上部到对岸逆冲停滑处的高差Δh为34 m,水平滑距Lmax为190 m,根据常用的架空坡理论公式[5]

    $$ {\mathrm{tan}}{\varphi _{\rm{r}}} = \frac{{\Delta h}}{{{L_{\max }}}} = 34/190 = 0.179 $$ (1)

    得出流通区滑带土的滑动摩擦角$ \varphi_r $为10.1°。

    滑坡到达流通区底部时滑速最大。高差取38 m,滑距取180 m,按式(2)[19]求得最大滑速 :

    $$ v = \sqrt {2g\left( {\left. {\Delta h - L \cdot {\mathrm{tan}}{\varphi_r}} \right)} \right.} = 10.6\;{\text{m/s}} $$ (2)

    另外,滑体从流通区底部急剧转弯75°进入主沟道时,弯道外侧超高(d)为4 m,转弯半径(r)为23 m,滑体宽度(b)为16 m。按照滑坡弯道处滑速计算式(3)[20]

    $$ v=(gdr/b{)}^{0.5}=(9.8\times 4\times 23/16{)}^{0.5}=7.5\;{\rm{m/s}} $$ (3)

    用以上两种方法求得永光1#滑坡前缘部分到达流通区底部时的滑速为7.5~10.6 m/s。实际上,本滑坡为黏滞性较高的流滑型黄土滑坡,上述两种方法计算的滑速均可能大于实际滑速。式(2)也未考虑滑动过程中滑体变形等耗能,计算滑速更大。

    图8可知,滑体在流通区下部与主沟交汇处受到对面斜坡区的迎头阻拦,发生逆冲爬高、大角度转弯和部分滑体突然制动堆积,滑坡运动状态发生了显著变化。但爬高堆积物边界较为规整,并未出现抛出、边缘溅起散落等现象,与本文分析和计算结果基本相符。而文献[15]认为永光1#滑坡在流通区底部的滑速达25 m/s。

    受地形条件控制,从流通区进入主沟道的大部分滑体转向75°,滑向下游主沟道。由于滑坡滑入主沟时的滑速较高,且大量堆积使沟道坡降加大,滑坡在主沟道中起初的滑速也较快,估计为4~5 m/s,否则在支沟与主沟交汇处将有更多的滑坡堆积。滑坡沿主沟下滑200 m后再转向为160°滑至沟口小桥处,滑速沿途逐渐减缓并堆积。

    地震当天几位地质灾害应急调查的专业人员在赶往永光1#滑坡灾害点的途中,中午12时观察到滑坡前舌距小桥还有80 m左右(图7),下午3时返回时滑坡前舌到达了小桥处。地震当天航拍的图2中也显示滑坡前舌还未到达小桥处。这为分析永光1#滑坡在主沟道中的滑动状态和过程提供了现场依据。从地震诱发滑坡发生到下午3时,滑坡在主沟道的滑动时间总体经历了约7 h,滑距590 m,平均滑速0.023 m/s。其中,前4 h的平均滑速为0.035 m/s,后3 h的平均滑速为0.0075m/s。可见,滑坡在主沟道中处于较为缓慢的流滑状态。而文献[15]模拟计算的永光1#滑坡的整体滑动过程仅为120 s,平均滑速8.3 m/s,差异甚大。其原因可能与对永光1#滑坡的实际滑动过程与特征调查分析不够,以及模型计算参数取值有关。

    综合分析永光1#滑坡上述滑动过程和滑速变化情况,建立滑动过程的地质概念模型(图10),显示出该滑坡经历了2次加速—减速的滑动过程。较客观地反映了复杂的滑动过程,为认识永光1#滑坡的滑动过程和机理,进一步进行数值模拟提供了基础资料。

    图  10  永光1#滑坡滑速分析图
    Figure  10.  Analysis of the sliding speed of the Yongguang 1# landslide

    永光1#滑坡在前期降水入渗和地震的耦合作用下发生,其远程滑动主要与地震、高含水率滑带土孔隙水压力升高、液化和特殊沟道地形有关。

    2013年岷县降水量高于年平均值,其中7月份高出当月平均值100 mm(图11)。降水增加坡体自重的同时,也有部分入渗补给地下水,对斜坡稳定性产生不利影响[21]

    图  11  2013年岷县降水量曲线
    Figure  11.  The precipitation curve for Minxian County in 2013

    永光1#滑坡距“7•22”岷县漳县Ms6.6级地震震中仅4 km,地震场地效应明显。据甘肃省地震局王谦模拟[22],该黄土斜坡地表加速度PGA达0.291~ 0.355 g。黄土孔隙发育且底部地下水较丰富,强烈地震使黄土结构破坏并产生体积剪缩,产生孔隙水压力,导致斜坡稳定性急剧降低,发生失稳滑动。结合滑源区地质条件进行动三轴试验,振动作用下饱和滑体产生的动孔隙水压力情况如图12[22]。动应力(σd)为20 kPa时,最大动孔隙水压力比(Ud/σ0')为0.53;动应力为 25 kPa时,试样破坏时动孔隙水压力比为0.60;动应力为30 kPa时,试样破坏时的动孔隙水压力比为0.41;动应力为35 kPa时,试样破坏时的动孔隙水压力比仅为0.31。随着循环振次的增长动孔隙水压力增长速率总体加快。

    图  12  滑体土动三轴试验曲线图
    Figure  12.  Dynamic triaxial test curve of the sliding soil

    滑动过程中饱和滑带土孔隙水压力进一步升高或发生液化,摩阻力大幅降低,饱和黄土滑带的孔隙水压力消散缓慢,持速效应明显。另外,滑源区后期发生的小型黄土滑坡叠加在滑坡洼地区,对中前部滑体产生推挤作用,也进一步加剧了其滑动。使永光1#滑坡在宽缓的平台区滑距达130 m,滑速3~4 m/s。也为滑坡前缘部分进入沟道继续滑动创造了条件,提供了动力。

    永光1#滑坡滑源区的滑动面受黄土与泥岩接触面控制,属黄土—泥岩接触面滑坡。甘肃黄土地区由降水诱发的此类滑坡的滑速较低[2324],滑距一般为数十米[2526]。“7•22”岷县漳县Ms6.6级地震加速了滑源区黄土斜坡的失稳和快速滑动。相对区内一般黄土—泥岩接触面滑坡,在降水和地震的耦合作用下,地震触发的永光1#滑坡在平台区的滑速和滑距均明显较大。

    滑坡前部流通区高差40 m、坡度14.5°,上部漏斗状开口正对滑坡前缘,延伸方向与滑动方向基本一致,顺直通畅,有利于滑坡的滑入和快速滑动。在较陡流通区的滑动过程中,高含水率滑带土因快速滑动进一步产生孔隙水压力并部分液化,滑坡再次加速滑动,在下部与主沟交汇处的滑速达7.5~10.6 m/s,是永光1#滑坡滑速最大的部位,也为滑坡沿主沟道继续滑动提供了动力。

    高含水率滑带土中产生的孔隙水压力,在永光1#滑坡的滑动过程中始终起到重要减阻作用。据试验,滑体的含水率普遍在25%~30%,洛阳铲取出岩芯呈泥状,敲击振动时液化严重(图13)。

    图  13  饱和滑体振动液化现象
    Figure  13.  Vibration-induced liquefaction in the saturated sliding mass

    滑坡沿主沟道滑动过程中孔隙水压力的作用更加明显。主沟道沟底平均坡度6.5°,最大为9.3°。而滑体在流通区滑动时,受孔隙水压力作用的影响,滑带土滑动摩擦角为10.1°。按照此滑动摩擦角,滑坡不会沿主沟道滑动,且滑体沿主沟道滑动时主地震的作用过程也已结束。所以,滑坡沿较缓主沟道的滑动存在更高的孔隙水压力。

    7月21日8时至7月22日8时地震发生时,永光村周围雨量站记录的降水量为11~24 mm,使地表土体处于高含水状态[17]。沟底主要由吸水性和膨胀性强、渗透性弱的古近系软弱泥岩组成,地震前期的大量降水和沟道洪水的浸泡,使沟道泥岩表层非常湿滑,摩阻力很低。当滑坡借助在流通区形成的较高滑速进入表部高含水的主沟道时,强烈的剪切和挤压使沟底浸泡软化的泥岩发生剪缩,在滑带土中进一步产生较高孔隙水压力,甚至发生部分液化,摩擦力大幅降低(图14),使其沿主沟道缓慢滑动。主震后频繁发生的余震(包括90 min后的Ms5.6级强余震)也有助于滑带土的液化和滑坡继续滑动。

    图  14  主沟道中的饱和滑带与液化现象
    Figure  14.  Saturated sliding zone and liquefaction phenomenon in the main channel

    降雨入渗条件下,圈闭的沟底地形、滑带土的低渗透性,易在滑带产生孔隙水压力且消散非常缓慢,托浮滑体缓慢向下滑动,持速效应更加明显(图15)。因此,在圈闭沟底中渗透性差且高含水率的滑道和滑体的有利组合产生孔隙水压力,降低有效应力和摩阻力,为滑坡的滑动提供了湿润的滑道,沿主沟道缓慢而持续地发生远程滑动。

    图  15  沟道中滑带孔隙水压力作用模式图
    Figure  15.  Schematic diagram of pore water pressure in sliding zone of the channel

    滑坡在沟道中缓慢滑动,说明其重力产生的下滑力略大于滑带土摩阻力,接近平衡状态,由于滑动状态的滑带土黏聚力近似于零,滑动摩擦角略小于沟道坡度。考虑滑带土孔隙水压力作用时,受力状态简化为:

    $$ W\sin\alpha > (W\cos\alpha-P\mathit{\mathrm{_W}})\mathit{\mathrm{tan\mathrm{ }}}\varphi $$ (4)

    式中:W——单位面积滑带上部滑体重量/kPa;

    Pw—单位面积滑带土孔隙水压力/kPa;

    α——滑道倾角/(°);

    $ \varphi $——滑带土摩擦角/(°)。

    由于Pw/Wcosα=ruru为孔隙水压力比。代入式(4)整理后得:

    $$ r_{\rm{u}}>1-\mathrm{tan}\alpha/\mathrm{tan}\varphi $$ (5)

    根据资料,α取6.5°,$ \varphi $取20°[15]。按照式(5)计算,主沟道段滑带土中的ru的平均值大于0.68。可见,孔隙水压力是永光1#滑坡在主沟道缓慢滑动的重要因素。

    随着滑坡在流通区和主沟道中累计740 m的远程滑动,该滑坡由最初的黄土—泥岩接触面滑坡转化为流滑型黄土滑坡,也是永光1#滑坡的重要滑动特征之一。

    综上所述,永光1#滑坡总滑距870 m,滑动总历时7 h,属远程非高速滑坡。近年来,黄土高原区此类滑坡时有发生,需引起重视和深入研究。

    (1)“7•22”岷县漳县Ms6.6级地震触发的永光1#滑坡总体积23×104 m3,滑坡主要由马兰黄土组成,平面形态呈不规则的“L型”。大部分滑体堆积在前部平台区,前缘6×104 m3滑体沿下部沟道远程滑动并不断堆积。可分为滑源区—平台堆积区和流通区—沟道堆积区两部分。

    (2)永光1#滑坡最大滑距870 m,前后缘高差与长度比为0.17,属远程滑坡;经历了2次加速—减速的复杂滑动过程,不同部位滑动特征差异较大;滑源区—平台堆积区滑距50~130 m,最大滑速3~4 m/s。流通区—沟道堆积区滑坡沿顺直较陡的流通区加速滑动150 m,最大滑速7.5~10.6 m/s;主沟道段滑距590 m,滑动持续时间7 h,平均滑速0.023 m/s。属远程非高速滑坡。

    (3)滑源区底部饱和黄土广泛分布,地震前连续降水和沟道洪水对沟底古近系泥岩的浸泡和软化,高含水率滑带土在地震和滑动过程中产生孔隙水压力及部分发生液化。圈闭的沟底地形和滑带土的低渗透性,产生更高孔隙水压力且消散非常缓慢,使全滑程摩擦力大幅降低,持速效应明显,是永光1#滑坡远程滑动的主要原因。

    (4)在滑源区—平台堆积区滑坡类型为地震和前期降水耦合作用触发的黄土—泥岩接触面滑坡;前缘部分滑体沿沟道发生远程流滑,形成特殊的长条状流通区—沟道堆积区,最终转化为流滑型黄土滑坡。

    致谢:地震当天到达滑坡现场的何文贵、王世宇、张永军教授级高工、苏永奇博士提供了宝贵照片和滑坡滑动情况的资料,在此一并致谢!

  • 图  1   落石试件

    Figure  1.   Falling stone test specimen

    图  2   固定慈竹

    Figure  2.   Diagram of fixed smash bamboo

    图  3   落石冲击慈竹

    Figure  3.   Falling rocks impact on smash bamboo

    图  4   试验全要素示意图

    Figure  4.   Schematic diagram of full-factor experimental design

    图  5   不同直径慈竹被砸折所需能量

    Figure  5.   The energy required to smash bamboo in different diameters

    图  6   不同壁厚慈竹被砸折所需能量

    Figure  6.   The energy required to smash bamboo in different wall thicknesses

    图  7   不同年龄慈竹被砸折所需能量

    Figure  7.   The energy required to smash bamboo at different ages

    图  8   不同长重比慈竹被砸折所需能量

    Figure  8.   The energy required to smash bamboo at different length-to-weight ratios

    图  9   不同长径比慈竹被砸折所需能量

    Figure  9.   The energy required to smash bamboo at different length-to-diameter ratios

    图  10   不同冲击位置慈竹被砸折所需能量

    Figure  10.   The energy required to smash the bamboo at different impact positions

    图  11   不同径厚比慈竹被砸折所需能量

    Figure  11.   The energy required to smash bamboo at different diameter-to-thickness ratios

    图  12   慈竹的破坏模式图

    Figure  12.   Diagram of the bamboo destruction pattern

    图  13   慈竹有限元模型

    Figure  13.   Finite element model of smash bamboo

    图  14   碰撞过程

    Figure  14.   Illustration of collision process

    图  15   碰撞过程应力变化

    Figure  15.   Illustration of changes in stress during the collision process

    图  16   落石能量变化曲线

    Figure  16.   Curve of falling rock energy variation

    图  17   慈竹总能量变化曲线

    Figure  17.   Curve of total energy variation in bamboo

    图  18   慈竹内能变化曲线

    Figure  18.   Curve of internal energy variation in bamboo

    图  19   慈竹动能变化曲线

    Figure  19.   Curve of kinetic energy variation in bamboo

    图  20   落石坡段概貌图

    Figure  20.   General view of the rockfall section

    图  21   落石运动轨迹图

    Figure  21.   Trajectory of rockfall movement

    表  1   慈竹直径 、壁厚试验结果

    Table  1   Experimental results of diameter and wall thickness of bamboo

    试验参数数值最大值/J最小值/J平均值/J样本/根平均值/J标准差变异系数修正系数标准值
    /J
    直径/mm新生竹35~40275.34197.90251.684251.6836.760.14610.8329209.61
    51~45283.95223.71253.1112253.1120.900.08260.9567242.15
    46~50309.76223.71255.4416255.4422.160.08670.9615245.6
    51~55344.18240.92268.0613268.0627.850.10390.948 0254.13
    老竹39~40240.92206.51223.713223.7117.210.07690.8844197.84
    41~45258.13215.11240.9211240.9215.390.06390.9647232.42
    46~51275.34223.71241.5415245.2218.320.07470.9656236.78
    壁厚/mm新生竹4275.34197.90250.608250.6026.220.10460.9293232.89
    5283.95223.71246.4514246.4519.860.08060.9614236.94
    6283.95240.92260.4315260.4313.980.05370.9753253.99
    7~9344.18240.92281.798281.7933.090.11740.9207259.44
    老竹4249.53206.51240.923238.0513.150.05520.917 0218.29
    5258.13223.71240.9215240.3518.240.07590.965 0231.94
    6~7275.34240.92258.1312243.7919.180.07870.9587233.73
      注:竹龄小于等于1 a的为新生竹,竹龄2~3 a的为老竹。
    下载: 导出CSV

    表  2   长重比、年龄和长径比试验结果

    Table  2   Experimental results of length-to-weight ratio, age and length-to-diameter ratio

    试验参数数值最大值/J最小值/J平均值/J样本/根平均值
    /J
    标准差变异系数修正系数标准值
    /J
    年龄/a新生竹1344.18197.90258.1344258.1324.950.09670.9752251.73
    老竹2~3275.34206.51241.5030241.5017.780.07360.9763235.77
    长重比
    /(m·kg−1
    新生竹0.9~1.3344.18240.92269.6015269.6028.660.10630.951 0256.40
    1.4~1.8275.34223.71254.1614254.4418.400.07230.9653245.63
    1.9~2.7283.95197.90250.6016250.6023.700.09460.958 0240.07
    老竹1.2~1.3275.34223.71245.709245.7022.810.09280.9419231.43
    1.5~1.6258.13215.11241.7810241.7814.880.06150.964 0233.07
    1.7~2.2258.13206.51237.7911237.7916.450.06920.9618228.71
    长径比新生竹65~70292.55240.92259.8510259.8518.500.07120.9583249.02
    71~76344.18249.53278.219278.2130.420.10940.9316259.17
    77~85283.95223.71253.2114253.2118.710.07390.9646244.25
    86~103275.34197.90247.3712247.3725.190.10180.9466234.17
    老竹65~70275.34223.71251.449251.4419.610.078 00.9512239.17
    71~79258.13223.71237.7911237.7913.480.05670.9687230.34
    81~97258.13206.51236.6210236.6218.250.07710.9548225.93
      注:竹龄小于等于1 a的为新生竹,竹龄2~3 a的为老竹
    下载: 导出CSV

    表  3   冲击位置和径厚比试验结果

    Table  3   Experimental results of impact position and diameter-to-thickness ratio

    试验参数数值最大值/J最小值/J平均值/J样本/根平均值
    /J
    标准差变异系数修正系数标准值
    /J
    冲击位置/m新生竹0.5~0.9344.18223.71260.8216260.8226.790.10270.9544248.92
    1.0~1.5275.34197.90251.5113251.5123.60.09390.953 0239.70
    1.7~2.6309.76223.71260.8216260.8224.690.09460.958 0249.85
    老竹0.6~0.9275.34215.11246.669246.6620.630.08370.9477233.75
    1.0~1.5266.74223.71243.799243.7915.510.06360.9602234.08
    1.6~1.9258.13206.51235.9012235.9017.000.072 00.9622226.99
    径厚比新生竹6~7309.76240.92268.898268.8920.440.076 00.9486255.08
    8344.18233.71266.7413266.7431.420.11780.9411251.02
    9283.95197.90250.2412250.2423.910.09550.9499237.71
    10~13275.34223.71249.5312249.5315.990.06410.9664241.14
    老竹7275.34223.71246.055246.0818.850.07660.9273228.19
    8275.34206.51234.6611234.6620.030.08530.9529223.60
    9266.74223.71245.847245.8415.600.06340.9531234.30
    10~12258.13223.71244.617244.6115.600.06380.9528233.07
      注:竹龄小于等于1 a的为新生竹,竹龄2~3 a的为老竹
    下载: 导出CSV

    表  4   不同直径不同厚度模拟结果

    Table  4   Simulation results for different diameters and thicknesses

    直径/cm厚度/m最大值/J最小值/J位移/cm总能量/J内能/J动能/J
    50.003240.86110.670−14.51145.53143.601.9327
    0.004247.7910.194−14.58187.17180.4316.878
    0.005269.1752.5460.12224.54211.9914.560
    0.006338.6337.17412.50314.85301.5213.328
    60.003276.50102.050−12.4086.2753.2333.035
    0.004276.5030.800−0.18158.52116.9541.572
    0.005314.6063.0860.1197.0873.51924.421
    0.006372.0618.0300.19155.13116.1938.930
    70.003299.0777.015−24.83107.9775.0032.939
    0.004306.7844.891−46.61103.4285.5917.833
    0.005330.5337.9100.10142.89126.1432.950
    0.006414.250.830−0.11177.09147.7029.392
    下载: 导出CSV

    表  5   物理模型试验与数值模拟对比

    Table  5   Comparison between physical model experiment and numerical simulation

    厚度/m慈竹破坏所需能量/J
    物理模型试验数值模拟
    0.003153.87145.53
    0.004198.55187.17
    0.005234.44224.54
    0.006328.62314.85
    下载: 导出CSV

    表  6   一丛慈竹消耗能量

    Table  6   The energy consumption of a bamboo cluster

    根/丛<4:4~6:6~8cm<44~66~8J/丛
    2621.38:58.15:20.1761555941.27
    1831.14:61.83:6.9261113975.55
    3072.18:26.99:0.8322805899.98
    下载: 导出CSV

    表  7   落石坡段信息

    Table  7   Falling rock section information

    序号水平投影长度/m竖向投影长度/m坡段长/m坡角/(°)
    19.809.800−0
    21.68.28.35578.959
    32.17.67.88574.554
    43.911.011.67170.478
    54.08.19.03463.719
    614.16.615.56825.084
    725.81.725.8563.770
    825.70.325.702−0.669
    下载: 导出CSV

    表  8   数据采集器结果

    Table  8   Data collector results

    参数最大值95%保证率值平均值
    速度/(m·s−110.4058.8484.177
    冲击能量/kJ7.7941.8490.468
    弹跳高度/m0.7660.6360.244
    采集点位置/m27.043
    横向威胁范围/m3.207
    下载: 导出CSV
  • [1] 叶四桥,陈洪凯,唐红梅. 危岩落石防治技术体系及其特点[J]. 公路,2010,55(7):80 − 85. [YE Siqiao,CHEN Hongkai,TANG Hongmei. Rockfall mitigation techniques and its characteristics[J]. Highway,2010,55(7):80 − 85. (in Chinese with English abstract)

    YE Siqiao, CHEN Hongkai, TANG Hongmei. Rockfall mitigation techniques and its characteristics[J]. Highway, 2010, 55(7): 80-85. (in Chinese with English abstract)

    [2] 铁永波,徐伟,向炳霖,等. 西南地区地质灾害风险“点面双控”体系构建与思考[J]. 中国地质灾害与防治学报,2022,33(3):106 − 113. [TIE Yongbo,XU Wei,XIANG Binglin,et al. The thoughts on construction of “double-control of point and zone” system of geological hazard risk in southwest China[J]. The Chinese Journal of Geological Hazard and Control,2022,33(3):106 − 113. (in Chinese with English abstract)

    TIE Yongbo, XU Wei, XIANG Binglin, et al. The thoughts on construction of “double-control of point and zone” system of geological hazard risk in southwest China[J]. The Chinese Journal of Geological Hazard and Control, 2022, 33(3): 106-113. (in Chinese with English abstract)

    [3]

    DORREN L K A,BERGER F. Stem breakage of trees and energy dissipation during rockfall impacts[J]. Tree Physiology,2006,26(1):63 − 71. DOI: 10.1093/treephys/26.1.63

    [4]

    DORREN L K A,BERGER F,IMESON A C,et al. Integrity,stability and management of protection forests in the European Alps[J]. Forest Ecology and Management,2004,195(1/2):165 − 176.

    [5]

    DORREN L K A,BERGER F,LE HIR C,et al. Mechanisms,effects and management implications of rockfall in forests[J]. Forest Ecology and Management,2005,215(1/2/3):183 − 195.

    [6]

    STOKES A,SALIN F,KOKUTSE A D,et al. Mechanical resistance of different tree species to rockfall in the French Alps[J]. Plant and Soil,2005,278(1/2):107 − 117.

    [7] 黄润秋, 刘卫华, 龚满福, 等. 树木对滚石拦挡效应研究[J]. 岩石力学与工程学报, 2010, 29(增刊1): 2895 − 2901

    HUANG Runqiu, LIU Weihua, GONG Manfu, et al. Study of trees resistance effect test on rolling rock blocks[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(Sup 1): 2895 − 2901. (in Chinese with English abstract)

    [8] 韩国刚,陈立新,程志华,等. 慈竹密度和微纤丝角变异规律研究[J]. 木材加工机械,2012,23(4):22 − 24. [HAN Guogang,CHEN Lixin,CHENG Zhihua,et al. Study on the Tsz bamboo density and the variation of microfibril angle[J]. Wood Processing Machinery,2012,23(4):22 − 24. (in Chinese with English abstract)

    HAN Guogang, CHEN Lixin, CHENG Zhihua, et al. Study on the Tsz bamboo density and the variation of microfibril angle[J]. Wood Processing Machinery, 2012, 23(4): 22-24. (in Chinese with English abstract)

    [9] 汪淑芳,梁梓,杨瑶君,等. 不同秆龄慈竹竹秆含水率与生物量的分配特征及其生物量模型的构建[J]. 西部林业科学,2013,42(1):42 − 45. [WANG Shufang,LIANG Zi,YANG Yaojun,et al. Biomass distribution characteristics and model construction of neosinocalamus affinis at different ages[J]. Journal of West China Forestry Science,2013,42(1):42 − 45. (in Chinese with English abstract)

    WANG Shufang, LIANG Zi, YANG Yaojun, et al. Biomass distribution characteristics and model construction of neosinocalamus affinis at different ages[J]. Journal of West China Forestry Science, 2013, 42(1): 42-45. (in Chinese with English abstract)

    [10] 杨喜. 梁山慈竹多尺度力学性能研究[D]. 长沙: 中南林业科技大学, 2014

    YANG Xi. The research on multi-scale mechanical properties of dendrocalamus farinosus[D]. Changsha: Central South University of Forestry & Technology, 2014. (in Chinese with English abstract)

    [11] 谢九龙,齐锦秋,周亚巍,等. 慈竹材物理力学性质研究[J]. 竹子研究汇刊,2011,30(4):30 − 34. [XIE Jiulong,QI Jinqiu,ZHOU Yawei,et al. A study on bamboo physico-mechanical properties of neosinocalamus affinis[J]. Journal of Bamboo Research,2011,30(4):30 − 34. (in Chinese with English abstract)

    XIE Jiulong, QI Jinqiu, ZHOU Yawei, et al. A study on bamboo physico-mechanical properties of neosinocalamus affinis[J]. Journal of Bamboo Research, 2011, 30(4): 30-34. (in Chinese with English abstract)

    [12] 谢九龙,齐锦秋,黄兴彦,等. 生长发育进程中慈竹秆形结构及物理力学性质[J]. 四川农业大学学报,2012,30(1):46 − 49. [XIE Jiulong,QI Jinqiu,HUANG Xingyan,et al. Culm form structure and physico-mechanical properties of neosinocalamus affinisin growth process[J]. Journal of Sichuan Agricultural University,2012,30(1):46 − 49. (in Chinese with English abstract)

    XIE Jiulong, QI Jinqiu, HUANG Xingyan, et al. Culm form structure and physico-mechanical properties of neosinocalamus affinisin growth process[J]. Journal of Sichuan Agricultural University, 2012, 30(1): 46-49. (in Chinese with English abstract)

    [13] 杜文军. 黔北地区主要竹种生长情况及物理力学性能初步研究[D]. 贵阳: 贵州大学, 2010

    DU Wenjun. Preliminary study on growth and physical and mechanical properties of main bamboo species in northern Guizhou[D]. Guiyang: Guizhou University, 2010. (in Chinese with English abstract)

    [14] 黄志良,刘亿,向波. 山区高速公路崩塌落石勘察设计浅析[J]. 公路交通技术,2010,26(1):12 − 15. [HUANG Zhiliang,LIU Yi,XIANG Bo. Analysis of survey and design against collapse and rockfall in mountainous expressways[J]. Technology of Highway and Transport,2010,26(1):12 − 15. (in Chinese with English abstract)

    HUANG Zhiliang, LIU Yi, XIANG Bo. Analysis of survey and design against collapse and rockfall in mountainous expressways[J]. Technology of Highway and Transport, 2010, 26(1): 12-15. (in Chinese with English abstract)

    [15] 陈其针,仲平,张贤,等. 构建中国自然灾害防灾减灾新体系[J]. 水文地质工程地质,2020,47(4):1 − 4. [CHEN Qizhen,ZHONG Ping,ZHANG Xian,et al. Establishment of an innovative system of natural disaster prevention and mitigation in China[J]. Hydrogeology & Engineering Geology,2020,47(4):1 − 4. (in Chinese with English abstract)

    CHEN Qizhen, ZHONG Ping, ZHANG Xian, et al. Establishment of an innovative system of natural disaster prevention and mitigation in China[J]. Hydrogeology & Engineering Geology, 2020, 47(4): 1-4. (in Chinese with English abstract)

    [16] 中华人民共和国建设部. 岩土工程勘察规范: GB 50021—2019[S]. 中国建筑工业出版社, 2002

    Ministry of Construction of the People’s Republic of China. Geotechnical engineering survey specification: GB 50021-2019[S]. China Architecture and Building Press, 2002. (in Chinese)

    [17] 任海青,张东升,潘雁红. 竹材抗压动态破坏过程分析[J]. 南京林业大学学报(自然科学版),2007,31(2):47 − 50. [REN Haiqing,ZHANG Dongsheng,PAN Yanhong. Dynamic compressive mechanical behavior of bamboo[J]. Journal of Nanjing Forestry University (Natural Sciences Edition),2007,31(2):47 − 50. (in Chinese with English abstract)

    REN Haiqing, ZHANG Dongsheng, PAN Yanhong. Dynamic compressive mechanical behavior of bamboo[J]. Journal of Nanjing Forestry University (Natural Sciences Edition), 2007, 31(2): 47-50. (in Chinese with English abstract)

    [18] 曹小军,李呈翔,魏素才,等. 四川慈竹生长现状调查与分析[J]. 世界竹藤通讯,2009,7(6):24 − 28. [CAO Xiaojun,LI Chengxiang,WEI Sucai,et al. Investigation and analysis of Neosinocalamus affinis growth condition in Sichuan[J]. World Bamboo and Rattan,2009,7(6):24 − 28. (in Chinese with English abstract)

    CAO Xiaojun, LI Chengxiang, WEI Sucai, et al. Investigation and analysis ofNeosinocalamus affinis growth condition in Sichuan[J]. World Bamboo and Rattan, 2009, 7(6): 24-28. (in Chinese with English abstract)

图(21)  /  表(8)
计量
  • 文章访问数:  90
  • HTML全文浏览量:  18
  • PDF下载量:  68
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-06-19
  • 修回日期:  2022-10-08
  • 网络出版日期:  2023-10-15
  • 刊出日期:  2023-10-30

目录

/

返回文章
返回