冰碛土体起动泥石流的特征研究

针对冰碛土体起动泥石流机理研究薄弱的现状，本文选取波密县帕隆藏布流域的支流嘎弄沟一冰碛土堆积坡面，通过模拟降水与冰雪融水起动冰川泥石流实验，比较 ...   第四纪研究  2019,Vol.39Issue(5):1235-1245   PDF     doi:10.11928/j.issn.1001-7410.2019.05.15   文章编号：1001-7410(2019)05-1235-11 冰碛土体起动泥石流的特征研究 陈宁生1, 王政1,2, 田树峰1,2, 朱云华3     (1中国科学院、水利部成都山地灾害与环境研究所,中国科学院山地灾害与地表过程重点实验室,四川成都610041;2中国科学院大学,北京100049;3西南石油大学地球科学与技术学院,四川成都610050) 摘要：与广泛分布于干旱河谷的宽级配砾石土体特征不同，冰碛土广泛分布在青藏高原地区，属粗大颗粒多、粘粒含量少、摩擦阻力大、粘滞阻力小的宽级配砾石土体。

在冰川融雪与降雨的共同作用下冰碛土体可失稳并起动泥石流，形成灾害。

针对冰碛土体起动泥石流机理研究薄弱的现状，本文选取波密县帕隆藏布流域的支流嘎弄沟一冰碛土堆积坡面，通过模拟降水与冰雪融水起动冰川泥石流实验，比较不同颗粒组成、不同实验条件下的土体起动泥石流特征，分析其起动成因及力学特性，探讨冰碛土体起动泥石流的机理。

研究发现冰碛土体失稳起动泥石流是粘滞阻力降低、孔隙水压力升高、拖曳力与渗流侵蚀共同作用的结果，起动过程受粘土颗粒含量和径流类型的影响。

当粘粒含量较高时（>3%），土体通过铲蚀与面蚀形成泥石流；粘粒含量中低时（不高于3%），大部分坡面土体主要经掏蚀与坍塌起动泥石流；粘粒含量过低时（<0.32%），土体难以起动泥石流。

在降水作用下土体孔隙水压力迅速增加，易造成土体破坏，起动泥石流；而在冰雪融水的作用下，土体孔隙水压力波动幅度不大时，土体同样可能发生失稳破坏起动泥石流。

关键词：冰碛土体    帕隆藏布流域    冰雪融水试验    力学特性    起动机理     中图分类号     P642.3                     文献标识码    A 第一作者简介: 陈宁生,男,54岁,研究员,从事山地灾害形成机理及其防治技术研究,E-mail:[email protected]. *国家自然科学基金项目（批准号：41671112）、国家自然科学基金国际（地区）合作与交流项目（批准号：41861134008）和四川科技计划项目重点研发计划项目（批准号：18ZDYF0329）共同资助 2019-04-25收稿，2019-06-27收修改稿 0引言 在极端气候和地震活动的影响下，泥石流灾害频繁发生。

目前暴雨泥石流的形成机制的研究较多，包括泥石流形成的土力类机理和水力类机理等[1~3]。

特别地，针对土力类泥石流，人们进一步揭示了其剪切破坏和液化过程与机理[4~5]，而对于冰碛土体起动形成冰川泥石流的研究却相对较少。

目前有关冰川泥石流的形成，针对以西藏林芝县念青唐古拉山脉东南缘培龙沟为代表的冰川泥石流，人们提出了融水溃决机制[6]；针对西藏波密县岗日嘎布山北坡的贡扎冰川末端米堆沟为代表的冰湖溃决类泥石流人们提出了水枕机制[7~10]；针对西藏波密县帕隆藏布流域古乡沟地震引发的特大规模泥石流，人们提出了冰滑坡转化为泥石流的机制[11~14]。

然而无论何种冰川类泥石流，其起动过程都是冰碛土体在降水或冰雪融水作用下失稳转化为泥石流的过程[15~17]。

由于冰川类泥石流发生区域的高寒、高海拔，其可到达性较差，观测研究十分困难，所以有关冰碛土体起动泥石流的力学机制方面的研究目前还十分缺乏。

而随着川藏铁路、川藏高速公路和雅鲁藏布江水能等重大工程的上马，冰川泥石流灾害防治的需要，使得对冰川泥石流形成的机理研究十分的需求和紧迫[18~21]。

冰碛土体最为显著的特征是其粘土颗粒含量普遍较少，仅在气候湿润区的后期风化作用较为强烈的地区，粘土颗粒含量可能超过3%(样品为小于60mm的颗粒)[22~24]。

研究表明在云南省北部金沙江流域蒋家沟等地区不同类型砾石土体其强度特征和起动泥石流的难易程度大相径庭[25~26]。

鉴于冰碛土体的组成和结构不同差异，所以笔者企图通过实验研究不同组分结构的土体起动泥石流的不同模式，并探究其孔隙水压力和含水量的变化过程[27~28]。

本文选择波密县帕隆藏布流域的支流嘎弄沟的冰碛物为代表，选择现代冰碛土体(粘粒含量低)、风化改造的老冰碛土体(粘粒含量中)、混合冰碛土体(粘粒含量中低)，在降水和冰雪融水作用下进行现场起动泥石流的实验，研究流域内不同类型冰碛土体起动泥石流的力学特征。

1研究方法 选择帕隆藏布流域的嘎弄沟不同类型冰碛土体作为实验对象，进行冰碛土体起动泥石流实验，实验对象分别为风化改造的老冰碛土体、现代冰碛土体和混合冰碛土体。

实验分为2组，即降水起动和径流起动泥石流，测量土体含水量变化和孔隙水压力变化特征，分析不同土体起动泥石流的难易特征和起动过程的孔隙水压力和含水量变化特征。

1.1研究区域和样品选择 帕隆藏布的支流嘎弄沟(29°57′57″N，95°43′27″E)流域面积5.8km2，由于嘎龙冰川发育，地质历史上冰碛物分布广[29~31]。

在1万年以来的全新世间冰期帕隆藏布大规模径流量的侵蚀下，残留的冰碛物较少，而由于流域的径流量相对小，而且随着冰川的退缩，嘎弄沟的冰碛物还源源不断地被侵蚀形成滑坡泥石流并危害扎墨公路[32]；导致形成了顶部为低粘粒冰碛土体、中部为中粘粒风化冰碛土体、底部为滑坡泥石流改造形成的中低粘粒冰碛土体的空间分布规律。

冰碛物的颗粒组分为燕山期花岗岩和片麻岩[33~34]，野外采集18组土样颗粒分析显示(最大粒径约60mm)，粘粒含量为0.32%~3.57%，粘土矿物主要为石英、斜长石、钠长石、拉长石等。

实验选择3类土体：1~4号样取于风化改造的冰碛土体，这类土体粘土颗粒含量较高，属中等粘土颗粒含量的土体，粘土普遍高于2.5%；9号样取于现代冰川碛物，其特点为粘土颗粒含量低，为0.32%；5~8号样、10~18号样为混杂冰碛土，粘土颗粒含量为中低，为0.5%~2.5%。

参照1号、10号、12号、5号、8号和9号样进行泥石流起动实验，1号样代表中等粘粒含量的风化改造的冰碛土体，9号样代表现代冰碛土体，其他样品为混合冰碛土体(图1)。

图1(Fig.1) 图1 颗粒分析累计曲线 Fig.1 Cumulativeparticle-sizegradationcurvesforsamples 1.2冰碛土体起动泥石流实验 根据试验条件及土体结构综合进行考虑，本次试验坡面的坡度取35°。

试验区(图2)位于墨脱公路K6+400处嘎弄沟左岸，其坡面组成具有大粒径、无分选、快速堆积等特点，且在坡面有少量残坡积物覆盖。

在坡面分别选取风化改造的冰碛土(中粘土颗粒含量)、混合冰碛土(中低粘土颗粒含量)、现代冰碛物(低粘土颗粒含量)，进行模拟降水及融水作用下泥石流起动对比试验，实验土体均为松裸砾石土体，土体的密度控制在1.68g/cm3左右，其颗粒组成见表1。

图2(Fig.2) 图2 研究区地理位置图 Fig.2 Mapofthestudyarea 表1(Table1) 表1土样颗粒组成 Table1Soilsamplecomposition 室内编号 颗粒成分(%) 砾粒(mm) 砂粒(mm) 粉粒(mm) 粘粒(mm) >60 60~20 20~10 10~2 2~1 1~0.5 0.5~0.25 0.25~0.075 0.075~0.05 0.05~0.01 0.01~0.005 <0.005 1 12.19 15.32 11.85 15.16 5.13 8.6 6.94 2.58 3.5 12.74 2.91 3.06 16 3.9 5.73 9.77 19.47 4.72 7.16 17.54 13.23 5.09 10.15 1.65 1.61 12 8.36 13.09 15.95 15.3 3.65 5.44 12.39 14.12 2.81 5.93 1.49 1.48 5 14.18 17.72 18.08 14.66 2.88 4.21 10.84 9.13 2.26 4.19 0.93 0.91 8 22.33 12.05 11.05 12.18 2.29 2.92 10.39 11.89 1.65 2.92 0.71 0.81 9 11.43 7.08 10.82 20.57 4.31 6.36 14.96 11.83 1.84 3.39 0.29 0.32 表1土样颗粒组成 Table1Soilsamplecomposition 在冰碛土松散堆积坡面，根据试验参数对坡面进行修整，划分一块尺寸约为1.5m×2.3m，厚度约为35cm，坡度约为35°的试验坡面，分两层在试验土体内埋设含水量及孔隙水压力传感器。

1号、2号和3号孔压及含水量传感器埋置在上层，埋设深度约为0.06m；4号、5号和6号孔压及含水量传感器埋在下层，埋设深度约为0.18m。

1号与2号、4号与5号的水平距离为0.75m(传感器的埋设方式及距离见图3)。

模拟降水实验采用人工降水进行，试验前对降水系统的雨量进行了标定，试验的降水强度为35.3mm/h。

与降水试验不同的是，模拟融水试验采用喷头水慢慢下渗进入坡体内部来模拟冰雪融水作用。

试验的潜水泵型号为QDX3-16-0.55，扬程为16m，实验的径流量为2.0m3/h。

图3(Fig.3) 图3 模拟降水与融水起动试验布置示意图 Fig.3 Schematicdiagramofsimulatedprecipitationandmeltwaterexperimentarrangement 2冰碛土体力学特性研究 泥石流的起动有三大必要条件即：水源、物源和坡降条件。

根据统计野外泥石流起动特征发现，在水源和坡面条件相同的情况下，有的土体能起动泥石流，而有的无法起动，可以推断是物源的结构组成对泥石流的起动产生了影响。

为探究土体的结构特征和力学特性，我们还进行了土体粒度分维计算、室内直剪试验和大型三轴试验[35]。

刘松玉和方磊[36]通过研究我国粘性土的粒度特征，发现在双对数坐标下粒径小于r的颗粒累计质量百分比含量(M(r)/M)与粒径r存在直线段，即表明土颗粒的粒度分布具有分形特征，并根据该直线段的斜率b，由公式D=3-b求出相应的粒度分维值D。

颗粒的分维值D反映颗粒在平面上的分布情况，一般来说，D越大，相应地级配越好，颗粒在平面上的密集程度越高，分选性越差；反之亦然。

由小于0.01mm的颗粒质量百分比含量得到拟合公式： Dp=-0.013x2+0.152x+2.244 计算结果分析发现，本文试验土样的分维值范围为2.044~2.683(图4)，与一般粘性泥石流沟源区土体平均分维值2.45~2.78[37]比较起来，数值相对较小，这是由于帕隆藏布流域发育的多为稀性泥石流，与粘性泥石流沟的土样相比，其细颗粒含量相对较低。

分维值的大小用来反映宽级配砾石土的颗粒组成，粗颗粒含量越高，表示其粒度分维值越小；而细颗粒含量越高，其分维值就越大[38]。

图4(Fig.4) 图4 粒度分维(D)与颗粒含量的关系曲线 Fig.4 Relationshipbetweengrainsizedistributionandparticlecontent 直剪试验中的土样取自野外模拟降水和融水试验中采集的土样，其中采用2mm的颗粒按粗颗粒土的土样制备(SL237-053-1999)方法，对土样中>2mm的颗粒进行等量替代。

试验设备采用四联直剪仪2J(DSJ-2)进行不固结快剪试验，剪切速率为0.8mm/min，每组试样分别在垂直荷载为50kPa、100kPa、150kPa和200kPa作用下进行水平剪切，剪切位移达至6mm时停止试验。

试验得出，在粒度分维值为2.396~2.658的区间内，冰碛土的粘聚力c随着粒度分维值的增大而增加，这是由于细颗粒占比越大，粗颗粒间骨架孔隙被填充的越多，从而增大了颗粒间的相互接触面积，故土体的粘聚力与粒度分维值是呈现正相关关系的。

而在这个区间，冰碛土的内摩擦角ϕ是随着粒度分维值的增大而减小的，这是由于细颗粒占比越大，粗颗粒逐渐被细颗粒包裹，粗颗粒间的摩擦力和咬合力减小，在剪切力作用下，颗粒间比细颗粒含量低时更容易发生剪切错动，使得土体的内摩擦角呈现出现降低的趋势。

本文采用四川大学水利水电学院华西岩土设备研究所研制的大型三轴仪(SZ30-2c)对不同干密度和不同级配的冰碛土进行固结不排水(CU)剪切试验，试验土样取自上述泥石流起动试验。

CU剪切速率为1.2mm/min。

试样直径为300mm，高600mm，最大粒径为60mm。

试样制样干密度分别为1.68g/cm3、1.73g/cm3、1.78g/cm3和1.83g/cm3。

以表1中12号试样(细粒含量为11.71%)为例，如图5a所示应力-应变关系曲线，取轴向应变εa为0.15时的应力为破坏应力，显示应力-应变为非线性关系，应力随着应变的增大而不断提高，曲线无明显峰值，表现为应变硬化型。

由图5a可知，围压作用越低，其轴向应变增长越明显；土的孔隙水压力-应变关系曲线(图5b)则表示为，在剪切初期，孔压随轴向应变增加而增加，当轴向应变达到某一值后，不同围压下的土样孔压均出现下降趋势，说明土体发生严重的剪胀现象。

这是由于在细颗粒含量较高的情况下，土颗粒在低围压下更容易发生错动和旋转从而产生体积膨胀，而在高围压下颗粒被压力所束缚，反而不易发生错动变形。

图5(Fig.5) 图5 12#土样的应力-应变(a)和孔压-应变(b)关系曲线 Fig.5 Stress-straincurveandporepressure-strainrelationshipcurveof12#sample 3不同土体的起动过程特征 模拟降水(A)与融水试验(B)显示不同类型的冰碛土体起动的过程和难易程度有很大的不同(表2)。

对于风化改造的冰碛土体(中等粘土颗粒旱后)，这类土体粘土颗粒含量较高，大于2.5%。

以下采用模拟降水试验第一组(A-1#)为例，可将其起动过程的特征分为7个阶段(表3)：1)在强降雨作用初期，雨水大量入渗导致土体含水量大幅增加，随即土体表面出现部分裂缝，产生蠕滑趋势，并附带形成少量细小冲沟；2)随着降雨的持续作用，土体入渗速率渐渐小于降雨速率，土体饱和度大幅增加，这时产生超渗透产流现象，形成地表径流层，进一步导致表层土的含水量、下滑力增加，粘滞力、抗滑力降低，稳定性下降；3)此时土体趋于饱和，且密度增大，超渗产流开始侵蚀表层细颗粒物质，在地表形成高含沙水流(图6)；4)由于土体的剪胀作用，使得孔压从降低变为增加，而含水量则仅有微小变化，表面裂缝继续扩张，土体的内摩擦角、摩擦阻力、抗剪强度及稳定性都降低；5)这一阶段初期土体有较小破坏，随着动水压力的作用，土体下滑力增加，之后开始整体滑动；6)当下滑力超过抗滑力时，土体从蠕滑转变为坍滑；7)土体滑动后粘滞力急剧下降，并在动水压力作用下带动更多土体滑动，从而形成大规模泥石流。

且从坍滑发生到泥石流起动，其体积相比初始状态增大了数倍。

这类土体的起动特征与目前研究的干旱河谷区多粘粒的土体起动泥石流的特征类似，为铲蚀+面蚀型起动，第一次蠕滑在本试验工况下为从降水开始300~400s[39~45]。

表2(Table2) 表2野外降水与融水径流条件下泥石流起动试验结果 Table2Theexperimentresultsofdebrisflowtriggeredbyrainfallandmeltwater 坡面条件 组别 粘粒含量(%) 是否起动泥石流 起动类型 中粘粒土体 A-1#/B-1# 3.06/3.01 是 土/水力类起动 A-2#/B-2# 1.61/1.86 是 水力类起动 A-3#/B-3# 1.48/1.03 是 水力类起动 中低粘粒土体低粘粒土体 A-4#/B-4# 0.94/0.81 是 水力类起动 A-5#/B-5# 0.81/0.68 是 水力类起动 A-6#/B-6# 0.32/0.32 否 难起动 表2野外降水与融水径流条件下泥石流起动试验结果 Table2Theexperimentresultsofdebrisflowtriggeredbyrainfallandmeltwater 表3(Table3) 表3中等粘粒含量冰碛土体起动泥石流起动过程 Table3Themediumclaycontentmorainesoilfailureanddebrisflowtriggeringprocess 过程 土体坍滑阶段 起动泥石流 雨水入渗 超渗产流 表面侵蚀 孔压升高 土体蠕滑 坍滑 起动及规模放大 特征描述 表面裂缝 产生超渗产流 坡面侵蚀 裂缝扩展 土体开始滑动 起动泥石流 体积增大 表3中等粘粒含量冰碛土体起动泥石流起动过程 Table3Themediumclaycontentmorainesoilfailureanddebrisflowtriggeringprocess 图6(Fig.6) 图6 降水条件下的试验坡面 Fig.6 Theslopeconditionduringtheartificialrainfallexperiment.(a)Soilfailureandtriggerdebrisflow;(b)Hyperconcentratedflow;(c)Gullyformedduringprecipitationtest 对于中低粘土颗粒含量的土体，其起动过程多表现为类水力类泥石流的起动过程，以掏蚀+坍塌型泥石流为主。

以模拟融水试验第三组(B-3#)为例，可分为6个阶段来阐述起动试验特征：1)在融水试验初期，细颗粒发生侵蚀运动，土体有小幅度蠕滑变形，且变形主要集中于破体的中上部；2)在地表径流的作用下，坡面会首先形成部分细小冲沟(图7)，随着径流的持续冲刷，通常会导致形成凹凸不平的纹沟；3)随着纹沟深度和宽度的不断增加，粗大颗粒逐渐开始侵蚀纹沟；4)径流汇集于此导致流量不断增加，并且径流下切侵蚀形成后缘陡坡；5)随着坡体稳定性越来越低，导致陡坡垮塌形成大量松散体并堆积于沟道内；6)径流的冲刷、携带及揭底作用导致冲沟内的固体物质参与径流运动，从而起动泥石流，在其运动过程中，坡体下部的松散物质不断参与其中，导致泥石流规模愈加庞大(图8)[45~47]。

图7(Fig.7) 图7 模拟融水起动试验坡面及坡面形成的细小冲沟 Fig.7 Theslopeconditionofdebrisflowtriggeredbyglaciermeltwaterandthesmallgulliesformedontheslope 图8(Fig.8) 图8 冰碛土体水力侵蚀起动泥石流过程 Fig.8 Debrisflowprocesstriggeredbythehydraulicerosionofmorainesoil 对于现代冰碛土体，其粘土颗粒含量低，通常小于0.5%(小于60mm的样品)，本文粘粒含量为0.32%，这类土体在实验所取的降水强度和径流强度下，不起动泥石流。

4不同冰碛土体起动泥石流的成因分析 不同类型的冰碛土体起动泥石流有两类成因，对于风化改造的冰碛土体，粘土颗粒含量中等，其成因为土力类；对于混合堆积土体，粘土颗粒含量中低，其形成成因为水力类。

粘粒含量较高的冰碛土体，例如A-1#试验土体的粘粒含量大于3%而小于5%，但实际仍属较低粘粒含量土体，该类土体抗剪强度高，渗透性好，但土中粘粒含量低，不能形成相对不透水面，导致雨水直接渗入坡底，无法在土中形成高孔隙水压力。

图9是其模拟降水过程中含水量与孔隙水压力随时间变化曲线，其显示在降水初期，土体的含水量与孔隙水压力均迅速增长；在随后的过程中，土中大量细颗粒被水流带走，孔隙水压力整体上升幅度并不大，却依然导致土体发生失稳滑动。

说明冰碛土体起动泥石流所需的临界粘粒含量并不高，孔隙水压力在泥石流起动过程中仍有一定作用。

图9(Fig.9) 图9 降水作用下风化改造的冰碛土体含水量及孔压变化曲线 Fig.9 Thechangeofwatercontentandporewaterpressureofmorainesoilindebrisflowprocesstriggeredbyrainfall 为巩固以上推测，我们进行了土体综合强度与受力的计算[48~50](表4)，根据边坡稳定的极限平衡理论对土体稳定性分析，引起土体坍滑的动力有土体表面超渗产流形成的拖曳力F1和雨水入渗过程产生的渗流压力F2，以及由土体重力产生的下滑力F3；而阻力则包括土体的抗剪强度即粘滞力c和摩擦阻力f。

当土体处于极限平衡状态时，有： (1) 表4(Table4) 表4土体综合强度与受力计算结果 Table4Soilcomprehensivestrengthandforcecalculationresults 指标力 估算依据 理论估算取值 估算结果(kN/m2) 拖曳力F1 F1= γw×l×J γw=1×103kg/m3，l=0.005m，J=tan35° 0.035 渗流压力F2 F2= γw×h×J γw=1×103kg/m3，h=0.005m，J=tan35° 0.035 下滑力F3 F3= W×sinθ W=1.02kN/m2，θ=35° 0.54 粘滞力c 0.28 摩擦阻力f f=(σ-p)×tanϕ Φ=32.5° 0.33 表4土体综合强度与受力计算结果 Table4Soilcomprehensivestrengthandforcecalculationresults 根据表4计算可得公式(1)左侧=0.61kPa，为使等式平衡，可知土体若要达到失稳起动泥石流，则土体粘聚力c需要降至0.28kPa。

根据极限平衡理论反演结果可知，在相对较高的粘粒含量条件下，冰碛土体失稳起动泥石流仍是孔隙水压力升高，粘滞阻力降低，拖曳力与渗流侵蚀共同作用的结果，其中粘滞阻力的降低占比45%，是最主要的影响因素。

对于更广泛分布于帕隆藏布流域的粘粒含量更低的冰碛土体，其坡面起动泥石流破坏也更具有代表性。

在低粘粒条件下，土体渗透性更强，地下水渗透速度更快。

如图10是B-5#试验模拟融水过程中含水量与孔隙水压力随时间变化曲线(只分析选取上层1号、2号、3号传感器数据)，在模拟融水试验初期，土中含水量快速增长而孔压变化不明显，随着水流下渗并出现地表径流，由此产生冲沟并导致上部土体蠕滑。

在地表径流和下渗水流的持续作用下，坡面发生坍塌，并呈片状向下推移起动泥石流。

针对土体内颗粒体在渗流力、重力、摩擦阻力及支撑反力作用下保持平衡[51~52](图11)，我们列出如下： (2) 图10(Fig.10) 图10 模拟融水作用下风化改造的冰碛土体含水量与孔压变化曲线 Fig.10 Thechangeofwatercontentandporewaterpressureofmorainesoilundertheweatheringbyglaciermeltwater 图11(Fig.11) 图11 渗流力作用下颗粒单元体受力图 Fig.11 Forcediagramofgranularunitunderseepageforce 公式(2)中，P为单位土体沿渗流方向所受到的渗透力(P=rwJ，其中rw为水的容重，J为水力坡降)，G为颗粒重力，f为摩擦阻力，hi为各力力臂长度。

在降水或融水初期，水流不断下渗，土体的含水量不断增加，土中细颗粒在自重、渗流力和摩擦阻力的作用下沿某一轨迹发生运移，并可能在粗颗粒构成的骨架中形成一条颗粒移动的通道。

随着细颗粒不断掏蚀移动，土体抗剪强度逐渐降低，一旦土体强度低于渗流下滑力时，土体发生坍塌下滑。

引起土体发生液态流动的临界水力坡降为： (3) 公式(3)中，Gs为土颗粒的相对密度，e为孔隙比[53~54]。

以模拟降水试验A-3#组为例，起动土体的下切深度为8.0cm，宽度5cm，长度为20cm，土体起动时需要的临界坡降为0.89。

5结论 本文针对不同粘粒含量下的冰碛土体进行了力学实验以及模拟降水和模拟融水试验，分析了不同冰碛土体泥石流的起动特征和成因，主要得出以下结论： (1)本文试验土样的粒度分维值范围为2.044~2.683，与一般粘性泥石流相比较小，表明冰碛土体的粗颗粒含量较一般粘性泥石流土体更大，土体结构更为松散，孔隙率较大。

直剪试验表明，粒度分维值在2.396~2.658区间时，冰碛土的粘聚力随粒度分维值的增大而增大，内摩擦角则与粒度分维值成负相关。

三轴试验则表明，细颗粒含量较高时，围压越低冰碛土体越容易发生剪胀现象。

(2)冰碛土体起动泥石流的过程是一个复杂的过程，针对不同粘土颗粒组成，泥石流起动的类型不同。

当粘粒含量较高时(>3%)，土体发生铲蚀+面蚀型泥石流起动，主要过程特征可分为7个阶段：包括雨水入渗、超渗产流、表面侵蚀、孔压升高、蠕滑、坍塌和起动放大阶段；当粘粒含量中低时(不高于3%)，大部分坡面泥石流起动以掏蚀+坍塌型为主，其特征可分为6个阶段：即细颗粒侵蚀、纹沟形成、粗颗粒侵蚀纹沟、形成后缘陡坡、垮塌和起动阶段；当粘粒含量过低时(<0.32%)，冰碛土体无法起动泥石流。

(3)对于土力类泥石流起动，在强降水条件下，冰碛土体孔隙水压力迅速增加，易造成土体破坏，从而起动泥石流，其失稳主要是受土体粘滞力降低的影响，本试验中由粘滞阻力降低造成的影响占45%；而对于水力类泥石流起动，在降水或融水条件下，当土体内孔隙水压力变化不明显时，土体同样可能发生失稳破坏起动泥石流，其主要受水力坡降的影响。

当水力坡降超过临界值后，会引发土体发生液态化流动，本试验中临界坡降为0.89时土体起动泥石流。

致谢:感谢审稿专家建设性的修改意见。

参考文献（References） [1] 崔鹏.泥石流起动条件及机理的实验研究[J].科学通报,1991,36(21):1650-1652.CuiPeng.Experimentalstudyontriggeringconditionsandmechanismofdebrisflow[J].ChineseScienceBulletin,1991,36(21):1650-1652. [2] 康志成.粘性泥石流稳定运动的力学分析[J].山地学报,1995,13(2):128-132.KangZhicheng.Dynamicanalysisofstablemovementinviscousdebrisflow[J].MountainResearch,1995,13(2):128-132. [3] 唐川,章书成.水力类泥石流起动机理与预报研究进展与方向[J].地球科学进展,2008,23(8):787-793.TangChuan,ZhangShucheng.Studyprogressandexpectationforinitiationmechanismandpredictionofhydraulic-drivendebrisflows[J].AdvancesinEarthScience,2008,23(8):787-793. [4] IversonRM.Debris-flowmechanics[M]//JakobandM,HungrO.Debris-flowHazardsandRelatedPhenomena.BerlinHeidelberg:Springer,2005:105-134. [5] ChenNingsheng,ZhuYunhua,HuangQi,etal.Mechanismsinvolvedintriggeringdebrisflowswithinacohesivegravelsoilmassonaslope:AcaseinSWChina[J].JournalofMountainScience,2017,14(4):611-620.DOI:10.1007/s11629-016-3882-x [6] 陈宁生,陈瑞.培龙沟泥石流及其堵江可能性探讨[J].山地学报,2002,20(6):738-742.ChenNingsheng,ChenRui.Glacialandrainstormdebris-flowinPeilongGullyandpossibilityinitsblockingmainriver[J].MountainResearch,2002,20(6):738-742. [7] 吕儒仁.泥石流与环境演变[J].山地学报,1997,15(2):91-96.LüRuren.Debrisflowandenvironmentalchange[J].MountainResearch,1997,15(2):91-96. [8] 程尊兰,田金昌,张正波,等.藏东南冰湖溃决泥石流形成的气候因素与发展趋势[J].地学前缘,2009,16(6):207-214.ChengZunlan,TianJinchang,ZhangZhengbo,etal.Debrisflowinducedbyglacial-lakebreakinSoutheastTibet[J].EarthScienceFrontiers,2009,16(6):207-214. [9] 党超,褚娜娜,丁瑜.冰湖溃决泥石流形成的临界条件[J].冰川冻土,2014,36(5):1176-1183.DangChao,ChuNana,DingYu.Thecriticalconditionfortheformationofdebrisflowinducedbyglaciallakeoutburstflood[J].JournalofGlaciologyandGeocryology,2014,36(5):1176-1183. [10] 赵万玉,陈晓清,刘建康,等.冰川终碛湖溃决-再生特征与机理[J].山地学报,2015,33(6):703-712.ZhaoWanyu,ChenXiaoqing,LiuJiankang,etal.Outburst-regenerationcharacteristicandmechanismofglacierlake[J].MountainResearch,2015,33(6):703-712. [11] 王培清,郭进军.藏东南山体滑坡成因分析及防治措施[J].西藏大学学报(汉文版),2006,21(1):33-36.WangPeiqing,GuoJinjun.PreliminarystudyontheanalysisofformationmechanismanddesignofcontrolforlandslideinsoutheasternTibet[J].JournalofTibetUniversity,2006,21(1):33-36. [12] 游勇,程尊兰,胡平华,等.西藏古乡沟泥石流模型试验研究[J].自然灾害学报,1997,6(1):54-60.YouYong,ChengZunlan,HuPinghua,etal.AstudyonmodeltestingofdebrisflowinGuxiangGully,Tibet[J].JournalofNaturalDisasters,1997,6(1):54-60. [13] 崔鹏,马东涛,陈宁生,等.冰湖溃决泥石流的形成、演化与减灾对策[J].第四纪研究,2003,23(6):621-628.CuiPegn,MaDongtao,ChenNingsheng,etal.Theinitiation,motionandmitigationofdebrisflowcausedbyglaciallakeoutburst[J].QuaternarySciences,2003,23(6):621-628. [14] 冯自立,崔鹏,何思明.滑坡转化为泥石流机理研究综述[J].自然灾害学报,2005,14(3):8-14.FengZili,CuiPeng,HeSiming.Mechanismofconversionoflandslidestodebrisflow[J].JournalofNaturalDisasters,2005,14(3):8-14. [15] 邓养鑫.冰碛转化为冰川泥石流堆积过程及其沉积特征[J].沉积学报,1995,13(4):37-48.DengYangxin.Processofaccumulationandcharacteristicsofglacialdebrisflowdepositstransformedbymoraine[J].ActaSedimentologicaSinica,1995,13(4):37-48. [16] 邓养鑫.冰碛与冰川泥石流堆积研究的若干新成果[J].冰川冻土,1996(S1):250-256.DengYangxin.Newachievementsofstudiesonthedepositsofmoraineandglacialdebrisflow[J].JournalofGlaciologyandGeocryology,1996(S1):250-256. [17] 铁永波,李宗亮.冰川泥石流形成机理研究进展[J].水科学进展,2010,21(6):861-866.TieYongbo,LiZongliang.Progressinthestudyofglacialdebrisflowmechanisms[J].AdvancesinWaterScience,2010,21(6):861-866. [18] KumarA,BhambriR,TiwariSK,etal.EvolutionofdebrisflowandmorainefailureintheGangotriGlacierregion,GarhwalHimalaya:Hydro-geomorphologicalaspects[J].Geomorphology,2019,333:152-166.DOI:10.1016/j.geomorph.2019.02.015 [19] 屈智炯,刘开明,肖晓军,等.冰碛土微观结构、应力应变特性及其模型研究[J].岩土工程学报,1992,14(6):19-28.QuZhijiong,LiuKaiming,XiaoXiaojun,etal.Studyofmicrostructures,stress-strainbehaviorandconstitutivemodeloftill[J].ChineseJournalofGeotechnicalEngineering,1992,14(6):19-28. [20] 崔鹏,韦方强,谢洪,等.中国西部泥石流及其减灾对策[J].第四纪研究,2003,23(2):142-151.CuiPeng,WeiFangqiang,XieHong,etal.DebrisflowanddisasterreductionstrategiesinWesternChina[J].QuaternarySciences,2003,23(2):142-151. [21] 方学东,黄润秋.青藏高原典型冰碛土的物理力学特性研究[J].工程地质学报,2013,21(1):123-128.FangXuedong,HuangRunqiu.PhysicalandmechanicalpropertiesoftypicalmorainesoilontheQinghai-TibetPlateau[J].JournalofEngineeringGeology,2013,21(1):123-128. [22] 谢春庆.冰碛土工程性能的研究[J].山地学报,2002,20(S1):129-132.XieChunqing.Theengineeringpropertiesofmoraine[J].MountainResearch,2002,20(S1):129-132. [23] 张斌斌.帕隆藏布流域海洋性冰川区泥石流特征研究[D].成都:西南交通大学硕士论文,2016:78-81.ZhangBinbin.StudyonDebrisFlowCharacteristicsinTemperateGlacierAreaofPallonTsangpo[D].Chengdu:TheMaster'sThesisofSouthwestJiaotongUniversity,2016:78-81. [24] 张振,王敬勇,吉锋,等.冰碛体力学性质的研究及探讨[J].人民黄河,2017,39(1):106-109.ZhangZhen,WangJingyong,JiFeng,etal.Researchandinvestigatethephysicalpropertiesofthemoraine[J].YellowRiver,2017,39(1):106-109. [25] 李树德.滑坡型泥石流形成机理[J].北京大学学报(自然科学版),1998,34(4):107-110.LiShude.Formationmechanismofthelandslide-typedebrisflow[J].ActaScientiarumNaturaliumUniversitatisPekinensis,1998,34(4):107-110. [26] 陈晓清,崔鹏,冯自立,等.滑坡转化泥石流起动的人工降雨试验研究[J].岩石力学与工程学报,2006,25(1):106-116.ChenXiaoqing,CuiPeng,FengZili,etal.Artificialrainfallexperimentalstudyonlandslidetranslationtodebrisflow[J].ChineseJournalofRockMechanicsandEngineering,2006,25(1):106-116. [27] ReidME,LaHusenRG,IversonRM.Debris-flowinitiationexperimentsusingdiversehydrologictriggers[C]//Debris-FlowHazardsMitigation:Mechanics,Prediction&Assessment,1997:1-11. [28] SpringmanSM,JommiC,TeysseireP.Instabilitiesonmoraineslopesinducedbylossofsuction:Acasehistory[J].Geotechnique,2003,53(1):3-10.DOI:10.1680/geot.2003.53.1.3 [29] 蔡祥兴,李椷,李念杰.帕尔提巴尔沟冰川泥石流的成因及其发展趋势[J].冰川冻土,1980,2(1):22-25+36.CaiXiangxing,LiHan,LiNianjie.ThegenesisanddevelopmenttrendofthedebrisflowinthePaltibargouGlacier[J].JournalofGlaciologyandGeocryology,1980,2(1):22-25+36. [30] 杨威,姚檀栋,徐柏青,等.近期藏东南帕隆藏布流域冰川的变化特征[J].科学通报,2010,55(18):1775-1780.YangWei,YaoTandong,XuBaiqing,etal.CharacteristicsofrecenttemperateglacierfluctuationsintheParlungZangboRiverbasin,southeastTibetanPlateau[J].ChineseScienceBulletin,2010,55(18):1775-1780. [31] 刘春玲,祁生文,童立强,等.喜马拉雅山地区泥石流发育特征研究[J].工程地质学报,2016,24(3):435-441.LiuChunling,QiShengwen,TongLiqiang,etal.DevelopmentcharacteristicsofdebrisflowinHimalayasusingremotesensingtechnology[J].JournalofEngineeringGeology,2016,24(3):435-441. [32] ZhangJinshan,ShenXingju.Debris-flowofZelongnongRavineinTibet[J].JournalofMountainScience,2011,8(4):535-543.DOI:10.1007/s11629-011-2137-0 [33] 吴锡浩,李永昭.青藏高原的冰碛层与环境[J].第四纪研究,1990(2):146-158.WuXihao,LiYongzhao.MorainesandenvironmentsinQinghai-XizangPlateau[J].QuaternarySciences,1990(2):46-158. [34] 吴中海,赵希涛,朱大岗,等.念青唐古拉山西布冰川区的冰碛层[J].地球学报,2002,23(4):343-348.WuZhonghai,ZhaoXitao,ZhuDagang,etal.ThemorainesofXibuGlacierareaintheNyainqentanglhaRange[J].ActaGeoscienticaSinica,2002,23(4):343-348. [35] 徐鼎平,汪斌,江龙剑,等.冰碛土三轴数值模拟试验方法探讨[J].岩土力学,2008,29(12):3466-3470.XuDingping,WangBin,JiangLongjian,etal.Studyofmethodsoftriaxialnumericalsimulationtestofglacialtill[J].RockandSoilMechanics,2008,29(12):3466-3470. [36] 刘松玉,方磊.试论粘性土粒度分布的分形结构[J].工程勘察,1992,20(2):1-4.LiuSongyu,FangLei.Discussiononfractalstructureofgrainsizedistributionofcohesivesoil[J].GeotechnicalInvestigation&Surveying,1992,20(2):1-4. [37] 黄祺,陈宁生,朱云华,等.泥石流源区砾石土的粒度分形特征[J].山地学报,2012,30(5):578-584.HuangQi,ChenNingsheng,ZhuYunhua,etal.Particlesizefractalcharacteristicsofthegravelsoilinoriginalareaofdebrisflowvalleys[J].MountainResearch,2012,30(5):578-584. [38] 倪化勇,刘希林.泥石流粒度分维值的初步研究[J].水土保持研究,2006,13(1):89-91.NiHuayong,LiuXilin.Apreliminarystudyondebrisflowgrain-sizefractaldimension[J].ResearchofSoilandWaterConservation,2006,13(1):89-91. [39] FuchuD,LeeCF,WangSJ.Analysisofrainstorm-inducedslide-debrisflowsonnaturalterrainofLantauIsland,HongKong[J].EngineeringGeology,1999,51(4):279-290.DOI:10.1016/S0013-7952(98)00047-7 [40] PasutoA,SilvanoS.Rainfallasatriggerofshallowmassmovements:AcasestudyintheDolomites,Italy[J].EnvironmentalGeology,1998,35(2-3):184-189.DOI:10.1007/s002540050304 [41] 胡明鉴,汪稔,张平仓.斜坡稳定性及降雨条件下激发滑坡的试验研究——以蒋家沟流域滑坡堆积角砾土坡地为例[J].岩土工程学报,2001,23(04):454-457.HuMingjian,WangRen,ZhangPingcang.Primaryresearchontheeffectofrainfallonlandslide-TaketheslopepiledbyoldlandslideinJiangjiagouvalleyasexample[J].ChineseJournalofGeotechnicalEngineering,2001,23(04):454-457. [42] 矫滨田,鲁晓兵,王淑云,等.土体降雨滑坡中细颗粒运移及效应[J].地下空间与工程学报,2005,1(7):1014-1016.JiaoBintian,LuXiaobing,WangShuyun,etal.Themovementoffinegrainsanditseffectsonthelandslideanddebrisflowcausedbyraining[J].ChineseJournalofUndergroundSpaceandEngineering,2005,1(7):1014-1016. [43] 赵彦波,游勇,柳金峰,等.粘性泥石流沟床冲刷深度试验研究[J].水利学报,2012,43(S2):92-97.ZhaoYanbo,YouYong,LiuJinfeng,etal.Experimentalstudyongullybederosiondepthofviscousdebrisflow[J].JournalofHydraulicEngineering,2012,43(S2):92-97. [44] 杨成林,陈宁生,邓明枫,等.粘粒含量对泥石流源区砾石土体强度影响的实验研究[J].成都理工大学学报(自然科学版),2011,38(5):522-528.YangChenglin,ChenNingsheng,DengMingfeng,etal.Experimentalstudyontheinfluenceoftheclaycontentonthegravelsoilmassfromtheupstreamareaofadebrisflow[J].JournalofChengduUniversityofTechnology(Science&TechnologyEdition),2011,38(5):522-528. [45] 王涛,陈宁生,邓明枫,等.沟道侵蚀型泥石流起动临界条件研究进展[J].泥沙研究,2014,1(2):75-80.WangTao,ChenNingsheng,DengMingfeng,etal.Progressinstudyofmechanismandcriticalconditionforinitiationofdebrisflowsduetobedfailure[J].JournalofSedimentResearch,2014,1(2):75-80. [46] MajorJJ,IversonRM.Debris-flowdeposition:Effectsofpore-fluidpressureandfrictionconcentratedatflowmargins[J].GeologicalSocietyofAmericaBulletin,1999,111(10):1424-1434.DOI:10.1130/0016-7606(1999)111<1424:DFDEOP>2.3.CO;2 [47] 曾凡伟,徐刚,李青.坡面泥石流发生的坡度阀值研究[J].地理科学,2005,25(2):244-247.ZengFanwei,XuGang,LiQing.Criticalgradientanddebrisflowonmountainslopes[J].ScientiaGeographicaSinica,2005,25(2):244-247. [48] IversonRM,GeorgeDL,LoganM.Debrisflowrunuponverticalbarriersandadverseslopes[J].JournalofGeophysicalResearch:EarthSurface,2016,121(12):2333-2357.DOI:10.1002/2016JF003933 [49] SassaK,WangG.Mechanismoflandslide-triggereddebrisflows:Liquefactionphenomenaduetotheundrainedloadingoftorrentdeposits[M]//JakobandM,HungrO.Debris-flowHazardsandRelatedPhenomena.BerlinHeidelberg:Springer,2005:81-104. [50] 李广信.高等土力学(第2版)[M].北京:清华大学出版社,2016:377-380.LiGuangxin.AdvancedSoilMechanics(SecondEdition)[M].Beijing:TsinghuaUniversityPress,2016:377-380. [51] 李广信,周晓杰.土的渗透破坏及其工程问题[J].工程勘察,2004,47(5):10-13+52.LiGuangxin,ZhouXiaojie.Seepagefailureofsoilanditsproblemsinengineering[J].GeotechnicalInvestigation&Surveying,2004,47(5):10-13+52. [52] 黄春娥,龚晓南.条分法与有限元法相结合分析渗流作用下的基坑边坡稳定性[J].水利学报,2001,32(3):6-10.HuangChun'e,GongXiaonan.CombinationofslicesmethodandFEMforslopestabilityanalysisoffoundationpitundertheseepagecondition[J].JournalofHydraulicEngineering,2001,32(3):6-10. [53] 黄永江,熊耀湘.渗流作用下土质边坡稳定性分析[J].水利科技与经济,2005,11(8):464-466.HuangYongjiang,XiongYaoxiang.Analysisofearthslopestabilityundertheseepagecondition[J].WaterConservancyScienceandTechnologyandEconomy,2005,11(8):464-466. [54] 钱家欢,殷宗泽.土工原理与计算(第二版)[M].北京:中国水利水电出版,1994:135-140.QianJiahuan,YinZongze.GeotechnicalPrinciplesandCalculations(SecondEdition)[M].Beijing:ChinaWater&PowerPress,1994:135-140. [55] 邓东平,李亮,赵炼恒.极限平衡理论下边坡稳定性抗滑强度参数反演分析[J].长江科学院院报,2017,34(3):67-73+79.DengDongping,LiLiang,ZhaoLianheng.Backanalysisofanti-slidestrengthparametersofslopestabilitybasedonlimitequilibriumtheory[J].JournalofYangtzeRiverScientificResearchInstitute,2017,34(3):67-73+79. Studyondebrisflowprocessinducedbymorainesoilmassfailure ChenNingsheng1, WangZheng1,2, TianShufeng1,2, ZhuYunhua3     (1KeyLaboratoryofMountainHazardsandSurfaceProcess,InstituteofMountainHazardsandEnvironment,ChineseAcademyofSciences,Chengdu610041,Sichuan;2UniversityoftheChineseAcademyofSciences,Beijing100049;3SchoolofGeoscienceandTechnology,SouthwestPetroleumUniversity,Chengdu610050,Sichuan) Abstract ThetriggeringmechanismsofdebrisflowwereexploredinthefieldbymeansofartificialrainfallandglacialmeltwaterexperimentintheGanonggou(29°57'57"N,95°43'27"E),ParlungZangbobasinandTibet,China.Thespatialdistributionofthemoraineinthisareais:thelowclaycontentmorainesoilisatthetop,themediumclaycontentweatheredmorainesoilisinthemiddle,andthemedium-lowclaycontentmorainesoilformedbylandslidesanddebrisflowtransformationisatthebottom.Theresultsofthesoiloffractalcalculation,directsheartestandlarge-scaletriaxialtestshows:Thegrainsizedistributionofmoraineisgood.Whenthegrain-sizefractaldimensioniswithintherangeof2.396~2.658,thecohesiveforceincreaseswiththevalueofgrain-sizefractaldimension,whiletheangleofinternalfrictiondecreaseswithit.Wefoundthattheclaycontentisacriticalfactoraffectingthetriggersofdebrisflowinmorainebyanalyzingtheartificialrainfallandtheglacialmeltwaterexperiment.Whentheclaycontentishigherthan3%,themajortriggersofdebrisflowinmorainecombinewithscouringerosionandsurfaceerosion;whentheclaycontentisnothigherthan3%,thedebrisflowinmoraineismainlytriggeredbyundercuttingerosionandcollapse;whentheclaycontentislowerthan0.32%,itisdifficultformorainetotriggerthedebrisflow.Thetriggersofdebrisflowinmorainecanbedividedintotwocategories:Oneisgravityviscoustrigger,whoseinstabilityismainlyaffectedbythereductionofmoraineviscosityresistance.Inthisexperiment,theinfluenceofthedecreaseoftheviscosityresistanceis45%;Theotherishydraulictrigger,whichismainlyaffectedbythehydraulicslope.Whenthehydraulicslopeexceedsthecriticalvalue,thesoilwillundergoliquidflow.Inthisexperiment,thecriticalslopetriggeringdebrisflowinmoraineis0.89. Keywords: morainesoil    ParlungZangbobasin    glacialmeltwatertest    mechanicalproperties    triggeringmechanism