地震扰动区内地质构造复杂,断裂活动性强,新构造运动强烈,导致地表岩体裂隙发育,易风化破碎,为泥石流活动提供了大量的物源。近年来,两处典型的泥石流灾害2008年M58.0级汉川地震区震后群发性泥石流和2010年8· 8舟曲泥石流均位于地震扰动区。水不仅是泥石流的组成物质,更是泥石流发生的重要触发因素。不同地区、不同季节、不同降雨量条件下,其泥石流源地堆积土体的侵蚀产沙以及对降雨响应的应力应变特点不同。强震后泥石流灾害的发生发展规律研究无论是对灾区重建,还是泥石流学科的发展都很重要。针对断裂活动区内特别是强震后泥石流的研究,众多学者开展了大量的研究工作。2008年汉川地震后,严重的灾情迫使许多研究者对震后泥石流的形成条件、活动特征等进行了深入的调查和分析,其中崔鹏等学者认为震后泥石流启动的临界雨量明显降低,泥石流形成将由降雨控制型逐步转为松散土体控制型,第一个泥石流活跃期可能会持续巧年左右;陈源井等也得出了相似的结论[司,认为小岗剑沟2010年暴发泥石流的临界雨量约为地震前的21%, 2011年暴发泥石流的临界雨量约为地震前的23% ,震后临界雨量大幅度下降。杨成林等根据汉川地震次生泥石流形成的地形地貌、降水和土源条件特征,认为低密度干燥度较高的土体在降雨作用下易湿陷,体积收缩,有利于土体孔隙水压力的升高和土体强度的降低,从而导致泥石流的启动。2010年8月8日舟曲泥石流发生后,许多研究者相继对该区泥石流的形成条件、活动特点、防治工程等进行了评估。张之贤等认为该地区前期干旱,突遇超历史极值的强降雨是此次灾害的主要原因,确定在不同地形和地质地貌背景下的泥石流降水量值对泥石流灾害发生的预报有着十分重要的意义。因此,本文以两个典型泥石流灾害频发区一都江堰汉川公路沿线和甘肃舟曲作为研究区,开展了区内泥石流流域内土体的结构特性的相关测试,并分析泥石流土体在不同含水率条件下强度参数的变化特征,以期探讨地震扰动区土体破坏形成泥石流的降雨量临界阑值。
1研究区环境地质条件
1. 1都汉公路沿线
都汉公路南起都江堰,北到汉川县城威州镇,地处青藏高原东部边缘崛江上游高山河谷区,区内山高坡陡,河流深切,地表起伏巨大,相对高差达3 000m以上。公路沿线地层除奥陶系、志留系大部缺失外,其它时代地层出露较齐全。岩性主要是花岗岩等岩浆岩和片麻岩、千枚岩、炭质板岩等变质岩。区内新构造运动活跃,断裂构造发育,龙门山中央断裂一映秀北川断裂穿越该区,断裂现今活动性强,龙门山前山断裂和后山断裂活动对该区亦有较大影响。龙门山断裂带白公元1169年以来,共发生破坏性地震25次,其中里氏6级以上地震20次。2 008年大地震前最近一次大地震为193 3年8月25日的茂县北部叠溪地震,震级7. 5级。2013年4月20日8时02分,该断裂再次活动,在四川省雅安市芦山县发生7. 0级地震,震源深度13 km。
5· 12汉川地震后,地处极震区的都汉公路沿线,崩塌、滑坡、滚石及不稳定斜坡广布,沟谷和山坡上形成大量的固体松散物质,为泥石流的发育提供了丰富的物质基础,沿线沟谷绝大多数具备了暴发泥石流的可能。根据震后多次调查,公路沿线多次因暴雨泥石流的破坏造成交通中断。
1. 2舟曲地区
舟曲县位于西秦岭构造带西延部分,地层岩性为中泥盆统古道岭组上段炭质板岩、千枚岩夹薄层灰岩和砂岩,下二叠统上段中厚层灰岩,上二叠统中厚层含硅质条带灰岩。受印支、燕山和喜玛拉雅山等多期造山运动的影响,境内褶皱强烈,断裂构造发育,光盖山一迭山断裂带经两条泥石流流域通过。光盖山一迭山断裂带及影响带宽度大,断裂近时期活动性较强,岩体裂隙发育,十分破碎。区内新构造运动十分活跃,表现为山地强烈隆升、沟谷深切,形成高山峡谷地貌。2010年8月7日晚,舟曲三眼峪和罗家峪沟同时发生特大泥石流灾害。两条泥石流沟流域内支沟发育,水系平面上均C7<‘树枝”状;三眼峪沟沟内有常流水,而罗家峪沟沟内无常流水。由于沟谷强烈侵蚀下切,横断面呈“V”字型或窄深的"U”字型。舟曲地区属地震强烈扰动区,多次古地震在沟内形成了体积巨大的堆石坝,是此次泥石流主要物源。5·12汉川地震时舟曲县处于6度烈度地区,地震扰动致使三眼峪沟和罗家峪沟流域内的岸坡岩体松动,坡体稳定性降低。
2松散堆积体工程特性
2. 1堆积土体物理性质
室内对两个研究区的堆积土样开展了颗分、含水率、液塑限、比重等测试。试验土样取白2008年汉川震后地震灾区典型的泥石流沟谷源区崩滑堆积松散土体和舟曲罗家峪沟泥石流汇水区内堆积体,都汉(以下简称DW)公路土体取于老虎嘴大滑坡附近的坡面泥石流灾害点,粗粒部分都汉公路沿线主要为花岗岩风化形成,舟曲(以下简称ZQ)罗家峪堆积体主要为沉积岩灰岩岩粉、岩屑、角砾,x粒含量较少。
(1)堆积体颗粒组成。
土颗粒的含量大于50%,可归为粗粒类土;其中粗粒土中大部分砾石含量》s。%,且细粒土中粉粒含量》50%,属于宽级配粉土质砾石土。依据图1计算土体的不均匀系数。和Cc,其中DW土样。为133 33,。为5. 33; ZQ土样。为46 67,。为a 9s。可见,两组土样均属于级配不良土体,粗粒土体含量明显高于细粒含量。比较两组堆积土体,都汉土体中勃粒(<a 。。sm m}含量为a 74%,舟曲土体中x粒含量为4 4s%,舟曲土体x粒含量高于都坟土体。
(2)堆积土体的物理性质。
两组泥石流土体的比重、液塑限、液性指数、塑性指数等物性指标的测试结果见表1。由表1可知,zQ土样比重较大,其液塑限、塑性指数均大于DW土样。一般,花岗岩抗风化能力强,风化速度较慢,但由于花岗岩节理裂隙发育,易发生崩解作用,加上水的广泛渗入,可形成深厚的风化带,且勃粒含量较高。但都汉公路地处龙门山中央断裂带影响区,2008年汉川地震释放的巨大能量使多个山体迅速崩解成岩块,新堆积的碎裂岩块尚未风化,因此该区域松散土体中的x粒含量较一般花岗岩地区偏少。
舟曲地区岩性多为石灰岩,主要由碳酸盐刁’一物组成,风化时主要是与水和二氧化碳作用产生的溶解和溶蚀,在干旱少雨的地区一般风化壳较薄}勃粒含量较低。但舟曲地区由于受多次历史地震运动的影响,多期次的泥石流堆积体堆积在沟道中,经过长时问的风化淋滤,松散土体x粒含量较一般石灰岩地区风化壳偏高。另外,该地区地层中尚有板岩、炭质板岩、千枚岩等变质岩。因此,ZQ土样细颗粒成分较高,土样粘粒(d动。。S mm)土颗粒的含量为4 04%;而以花岗岩风化碎屑为主的都汉土样x粒含量仅为a 74%。所以,从两种堆积土体基本物理力学性质分析结果可初步推测,舟曲地区堆积土体启动形成泥石流的含水率应高于都汉公路沿线堆积土体。
2. 2堆积土体强度特性
为研究土体变形破坏时的界限含水率, 采用常规直剪仪开展不同含水率条件下土体的不固结快剪试验, 分别测试其在不同含水率下强度指标的变化趋势。试验土样采用粒径 [2 mm 的泥石流源地土体, 按其天然密度制作不同含水率的重塑土样, 土样采用分层击实, 以保证土样与剪切盒间隙充分填充,试验过程中土样受力均匀。
( 1) DW 土样不同含水率的抗剪强度。
都汶公路沿线松散堆积土按干密度 2. 0 g/ cm3分别配置含水率为 7%, 12%, 15%, 18% 的重塑土样, 进行土样不排水剪切试验, 试验结果见图 2、图 3及表 2。
由图2可知,正应力6与剪应力T呈线性正相关关系,施加给土样的6越大其发挥的T也越大。
图3显示了土的r}x聚力与含水率关系,当含水率由7%增加到12%时,勃聚力。增加;含水率增至巧%时,。值急剧减小。一般,随着含水率的增加,土体x聚力应该逐渐减小,直到完全丧失强度。其根本原因在于水在较粗颗粒之问起润滑作用,使粒问摩阻力降低;另一方面x土颗粒表面结合水膜增厚使原始x聚力减小。但在此试验初始阶段,当含水率由7%增至12%时,勃聚力反而由15. 14 kPa增加至34. 67 kPa,这可能是由于其土体级配、岩性因素引起。该土样属于不均匀宽级配粉土质砾石土,x粒含量少,接近砂土的性质。干燥的砂土颗粒问无x聚力,但当砂性土在较低含水率条件下会表现出假凝聚力,假凝聚力是由于土颗粒接触面上一些水的毛细压力所形成。当松散土体完全干燥或完全饱和时,毛细压力消失。因此,当试样含水率由12%增加到15%时,。值变化非常明显,勃聚力从34 67 kPa急剧下降至2 91 kPa。暗示了都汉泥石流源地土体在降雨作用下失稳的界限含水率可能就在12%和15%之问。因此可以判定DW堆积土体遭遇降雨时,当雨水入渗土体至此含水率区问,土体强度快速降低,在合适的地形坡度、水动力条件下易失稳启动形成泥石流。
与c值不同,含水率从7%增加到15%的过程中,内摩擦角φ在260 - 270之问波动,含水率增加至18%时,φ下降至20.15。总体来说,试样内摩擦角φ随含水率变化无明显下降趋势,变化不大,呈波动趋势,表现出接近砂性土的特性。因此,含水率的变化对都汉堆积土体的φ影响不大;再者,在试验过程中,土样为重复试验,制作的试样土粒不均匀,对测得的φ值也有一定影响,但主要原因还是试样x粒很少。可见在降雨作用下,松散土体颗粒之问x聚力决定坡面泥石流土体的强度,当遭遇强降雨作用时,极易爆发泥石流灾害。
(2)Z Q土样不同含水率的抗剪强度。
ZQ土样按干密度1. 9 g/ cm;分别配置含水率为7% , 12% , 15% , 18% , 21%的重塑土样,进行土样不排水剪切试验,试验结果见图4、图5及表3。
试验结果表明,ZQ土样6与T呈线性正相关关系厂图5),且ZQ土体各含水率条件下。值变化与DW土样具有相似的规律。当土样含水率由7%增至12%时,受土粒问毛细压力的影响,勃聚力、:增高;当含水率由15%增加到18%时,勃聚力从26. 89 kPa急剧下降至8 41kPa。据此可以判定ZQ土体在此含水率区问,强度易快速降低。说明舟曲泥石流源地土体在降雨作用下失稳的界限含水率可能就在15%和18%之问。与 DW土样类似,ZQ土样各试样内摩擦角φ随含水率变化无明显趋势,变化不大,呈波动趋势,主要原因是该松散土体主要由石灰岩岩粉、岩屑、角砾等组成,勃粒含量较少。
3泥石流启动机制
3. 1地震扰动区降雨诱发泥石流机制
众多学者对降雨诱发泥石流的启动机制进行了大量深入的研究。崔鹏等[川通过大量水槽实验对泥石流起动机理进行了详尽的研究,阐述了泥石流起动的加速机理,分离机理,和连接机理,强调了水分饱和度、细颗粒、坡度三者在泥石流形成过程中的重要作用。日本学者高桥堡将泥石流启动的原因归纳为水流影响沟床堆积体强度、坡面崩塌滑坡土体与水作用演化为泥石流以及崩滑体堵塞沟道受上游汇水影响形成泥石流三种情况,汉川地震区震后泥石流研究结果表明,震后流域下垫面条件巨变以后,流域的土体入渗条件、坡面径流条件、流域汇流条件和沟道洪水过程的变化,形成有利于起动、汇流和规模发展的水文条件。
对地震扰动区而言,地震、断裂活动使得震区松散堆积岩土体方量大大增加,且堆积体结构松散,孔隙率大,力学强度低,多处于失稳的临界状态。地震扰动区坡面上或沟道里的松散岩土体,受短时问强降雨的作用后,其强度丧失进而形成泥石流的机制主要体现在两个方面。首先,大量雨水迅速入渗并不能及时排出,导致岩土体中孔隙水压力迅速升高,土体抗剪强度降低,并增大土体的白重力。土体局部出现剪切破坏,出现张裂隙和剪切裂隙,并使土体中原有的裂隙、孔隙等扩大。随着孔隙水压的进一步增大,岩土体有效应力逐渐减小,剪切破坏区扩展连通形成剪切面,沿剪切面发生渗流产生渗透力,进一步促进土体强度弱化。当有效应力降低为零时,局部地段松散土体呈悬浮状态。当剪切面上的剪切力大于抵抗力时,斜坡岩土体启动下滑,下滑的土体碰撞、剥离解体,下泻造浆形成泥石流。另一方面,若松散土体覆盖于坚硬基岩上,由于土体松散,下伏基岩坚硬完整,松散岩土体与基岩之问接触面为力学性质差异明显的界面,由暴雨产生的强烈坡面水流沿此接触面发生渗流,产生一个向下的渗透力,这大大增加松散岩土体的下滑力,从而在合适的坡度下启动形成泥石流。另外,松散岩土体在坡脚沟道大量堵塞,影响沟道的过流能力和径流过程,一旦堵塞体溃决,直接冲击下游堵塞坝,逐级破坏,其级联效应使得物质和能量不断累积,形成大规模泥石流。
3. 2临界雨量与松叙土体临界含水率的关系
对于泥石流发生时激发雨量的研究是近年来泥石流灾害预警研究的一个热点。众多学者通过前期降雨或当时降雨的有效降雨指标与泥石流发生关系的试验和统计,建立了一系列基于雨量和雨强的泥石流预测预报模型。近年来,对震后泥石流启动的临界雨量研究也取得了丰富的成果,特别是临界雨量降低程度的定量化分析方面。胡凯衡等采用逐步回归分析方法,建立了临界土体含水量与土体渗透系数、孔隙度和颗粒曲率系数的经验关系,进而提出一种基于临界土体含水量和实时降雨的泥石流预警方法。一般认为,由于受到地表径流转化、土壤入渗、水分蒸发、植物吸收等诸多因素的影响,前期降雨对泥石流发生影响存在一个衰减问题,距泥石流发生时问越久的降雨,对泥石流发生的影响越小。据相关研究,我国西南地区泥石流在发生过程中,前期降雨的贡献一般小于当时雨量的贡献。
研究结果表明,泥石流启动与灾害发生前1h降雨量密切相关。郭小军等根据都汉公路南段200 8年一2011年问的22次泥石流事件及其对应的降雨过程,确定该区域平均1h激发雨强为14. 8mm/h。相关研究表明舟曲地区沟谷型的泥石流1h临界雨量仅为12 mm}z}}。而根据本次试验结果,两组土样x聚力急剧降低的含水率界限分别为12%一15%和15%一18%,为了便于统计分析,取含水率区问的低值作为土体失稳的含水率值,即都汉公路土样临界含水率取12%,舟曲土样取15%。
另外一些研究者也提出了相似的认识。如四川省德阳市绵竹县清平乡文家沟在2011年8月13日泥石流的1h激发降雨量约为70. 6 mm,文家沟泥石流土体的试验分析结果建议将土体含水率15%一18%作为土体界限含水率的参考} zz}陕西宁陕地区坡面泥石流的研究表明当土体含水率超过28%一30%,勃聚力、内摩擦角与含水率关系曲线都出现了明显的拐点。而据全国各地平均降雨量分区的限界值,四川西部1h临界降雨量为15mm,陕西南部泥石流1h降雨量约为20 m m。而汉川地震后,北川县泥石流临界小时雨强降低了25.4%一31. 6%}ZS}。因此绵竹清平区域1h临界降雨量取15 mm减少25. 4%即11. 2 mm。另外,云南东川蒋家沟的泥石流砾石土体内摩擦角φ值在含水率达到7. 07%时具有临界特征,当含水率、由5. 35%增加到7. 07%时,φ值略有上升,而后φ值随含水率增加而减小。根据云南小江流域蒋家沟的降雨资料与29次泥石流过程的对应关系,取泥石流1h激发降雨量约为6 mm 。
综合以往的研究成果,以上述五个区域实际发生泥石流灾害的1h降雨强度和土体含水率进行统计(表4),结合都汉和舟曲松散土体。值急剧降低的含水率试验结果,对表4进行相关分析,得到临界含水率和临界降雨量的经验公式(式(1))如下:
式中:x为松散土体x聚力。急剧降低的界限含水率`为泥石流暴发1h临界雨强(m m} 。
实际开展泥石流灾害防治时,可依据式(1)对断裂带其他影响地区松散土体进行取样分析,得出其x聚力急剧降低的含水率后,进一步推求其泥石流启动的1h临界雨量,为相关区域的泥石流预警预报提供依据。
4结论与讨论
(1)本研究选取的两处地震扰动区泥石流堆积土体为级配不均粉土质砾石土,勃质含量低。舟曲土样因其母岩为石灰岩、千枚岩、板岩,较都汉公路沿线花岗岩风化程度更深,导致其松散堆积体重细粒含量相对较高,因而其比重、液塑限、液性指数、塑性指数均较大。
(2)泥石流源地土体由于其细颗粒含量少,导致土体在低含水率时,受土颗粒问毛细压力的影响,勃聚力。值先随含水率增加而增加,随后在一定含水率区问。值急剧减小。都汉公路沿线坡积土在12%一15%含水率区问急剧降低,舟曲泥石流堆积土在15%一18%含水率区问降幅最大。暗示了当堆积土体含水率达到一定值(临界值)时土体强度丧失,堆积土体会启动形成泥石流。因此泥石流灾害监测预警过程中应密切关注土体含水率变化,以便及时采取相应措施,避免泥石流灾害。
(3)统计结果表明,泥石流发生前1h的临界雨量与坡面土体强度聚降的临界含水率之问存在密切相关性。因此,对泥石流易发地区或特定泥石流沟坡面土体强度与含水率关系的研究,可为泥石流灾害预警提供科学依据。
( 4)采用快剪试验得到两种堆积土体的强度参数,进一步可采用三轴试验进行结果对比,以提高其测试精度。
(5)文中拟合的泥石流启动前1h雨强经验公式的数据样本只有5个地区,公式尚需进一步修正。都汉公路沿线土体仅选用坡积土为实验样本,可继续选用沟谷泥石流的源地土地进行对比分析。因此,开展不同地区泥石流堆积土体强度与含水率的研究,以提高经验公式的样本数和精度是很有必要的。