狭陡型泥石流成灾特征与防治对策研究——以阿坝州金川县刘家沟泥石流为例

崔华丽，杨东旭，杨 栋，李元灵，贾 逸，吴 波*

（1.中国地质科学院探矿工艺研究所，四川 成都 611734；2.成都华建地质工程科技有限公司，四川 成都 611734）

0 引言

2008年汶川8.0级地震诱发了大量的崩塌、滑坡，据粗略估算区内新增松散固体物源量达55.86×108t［1］，为泥石流的活动提供了丰富的物源。受地震影响，大量沟道从山洪沟转化为泥石流沟，针对重震区、（特）大规模的泥石流灾害，众多学者从泥石流的暴发成因、泥石流特征、成灾（堵河）分析、防治经验等方面均进行了大量的研究［2-7］。其中一类“狭陡型泥石流”的特征指标、防治对策以及动力学特性与震前常规泥石流存在显著差异，烧房沟、小岗剑、磨子沟、连山大桥沟、福堂沟、蒋家沟均属于此类。屈永平等［8-9］、杨东旭等［10］在对汶川震区泥石流灾害的调查与分析中得出，狭（急）陡（窄）沟道泥石流流域面积＜5 km2，沟道纵比降＞300‰，沟道横断面以V形为主，流域内单位面积物源量丰富，区域强降雨集中，活动旺盛。韩玫等［11］以都汶高速沿线磨子沟为例，分析了汶川震区窄陡沟道型泥石流动力学过程，并分析其危险范围、评价其危害性，为此类泥石流防治工程措施提供了依据。杨东旭等［12］以小岗剑泥石流灾害为例分析了强震-暴雨叠加区泥石流在形成、活动及危害方面的典型特征，并提出了相应的防治对策。陈明等［13］以汶川瓦窑沟为例，研究其动力特性及启动、侵蚀、堆积演化过程，提出应开展对震区急陡沟道型泥石流风险评估、监测、早期预警等研究，治理方案也需有效组合运用。

狭（急）陡（窄）型泥石流暴发时间短、冲出速度快、冲出距离远、成灾性强，危害大，因此有必要针对单条狭陡型泥石流沟道的特点，深入分析治理难点，提出针对性、适宜性的防治新技术。近年来高频小流域、成灾性与危害性较弱的狭陡型泥石流受关注较少，在灾害防治方面未引起重视，在治理经验方面还存在一些不足。本文以阿坝州金川县刘家沟泥石流为例，分析该泥石流沟的成灾机理、堵江可能性并结合沟道情况提出适宜的治理措施，以期为同类型泥石流防治提供思路。

1 研究区概况

1.1 流域概况

刘家沟位于阿坝州金川县卡撒乡脚姆塘村，为卡撒沟的支沟。流域面积0.36 km2，主沟长1.43 km，沟域内最高点高程2710 m，沟口高程2124 m，相对高差达566 m，沟道平均坡降409.8‰。刘家沟沟域内无支沟发育，且无清水汇流区，遥感图像显示为白色条带状（见图1）。沟域内植被不发育，流域形态指数δ=A L2定量计算［14］为0.168。

图1 刘家沟流域全貌Fig.1 Liujia gully drainage area

1.2 地貌环境

刘家沟沟道切割较深，跌水和陡坎较发育，沟道两侧地形坡度在50°左右，局部可达65°以上，整个沟域由形成流通区、堆积区组成。形成流通区在2161～2710 m，区纵坡432.28‰，沟道较窄，沟底宽约3～5 m（邻脊宽度170～280 m），呈典型V字形深切沟谷，两侧斜坡在30°～50°，局部可达70°，此段滑坡堆积体和残坡积体发育；2161 m至沟口段纵比降111.45‰，地形平缓，沟道较浅，为泥石流的堆积区。

1.3 地层岩性

研究区内地层主要为中生界三叠系中统杂谷脑组（T2z），岩性为浅灰、灰、深灰色中层、厚层及块状变质长石石英砂岩、石英砂岩、凝灰质砂岩、粉砂岩夹灰色粉砂质板岩、炭质板岩。砂、板岩之比为5∶1～7∶1，砂岩以单层厚、分布集中为特征。

1.4 地质构造与地震

刘家沟流域位于滇藏歹字形构造头部外围与金汤弧形构造西半弧复合部位。其构造的显著特点是：一方面区内褶皱构造发育，断裂不发育，以挤压紧密、轴部尖棱为特征之褶皱为主，一个大的褶皱往往由数个至数十个次级褶皱组成，一个次级褶皱，又包含多个揉皱、挠曲，断裂主要有向阳断裂和烧日断裂，皆为推测断裂；另一方面，两构造体系之构造形迹表现形式不一，并在平面上构成明显分区，其中歹字形构造体系为区域性强大的压应力作用的产物，由一系列强烈挤压之线性褶皱组成，构造线呈北西线展布，而金汤弧形构造体系为局部应力场的产物，应力相对较弱，主要由一系列强烈挤压之紧密褶皱组成，由北向南褶皱紧密程度递增，构造轴线呈东西向、北西西向，地层多以塑性变形为主。

研究区新构造运动十分强烈，第四纪以来主要表现为强烈升降运动和断裂活动。川西北断块的3条边界——托索湖-玛沁-文县断裂带、龙门山断裂带和鲜水河断裂带均为活动断裂。晚近期以来地震活动十分强烈且频繁，这些地震主要分布于川西北断块的边界断裂带上。流域地处鲜水河地震带和松潘地震带之间，抗震设防烈度为7度，地震加速度值为0.10g。刘家沟沟道纵、横断面及物源如图2、图3所示。

图2 刘家沟沟道纵断面及物源Fig.2 Longitudinal-section and source of Liujia gully

图3 刘家沟沟道典型横断面Fig.3 Cross-section of Liujia gully

1.5 气象气候

刘家沟流域属明显的大陆性高原季风气候，流域内为高山峡谷区，地形复杂，相对高差大，气候呈立体分布，随海拔高度、地形不同而有差异。区内大致可分为2种气候类型区：（1）暖温带季风气候，在海拔2400 m以下，活动积温3400～4500℃，年降雨量606.8～654.4 mm；（2）在海拔2400～2800 m中温带季风气候，活动积温3400～1600℃，年降雨量654～745 mm。金川县泥石流发育地区的特征为1 h降雨量≥10 mm，24 h降雨量≥30 mm时就有发生泥石流的可能。

2 泥石流发育特征及致灾机制

2.1 物源与势能

泥石流发生的势能条件与泥石流流域地形地貌关系密切，金川县陡峻的地形为泥石流发生提供了良好的势能和动能条件，主要表现在沟谷流域高差一般较大，沟谷为典型的深V形谷，使得沟道内或两岸松散堆积物具有较大势能，有利于向动能转化进而引发泥石流。

刘家沟形成流通区沟道两侧发育有5处滑坡以及5处坡面倾蚀物源，上游泥石流启动区（海拔2710～2335 m），沟 道 长801 m，平 均 纵 坡 为468.16‰，此段沟道在泥石流侧蚀作用影响下，概算大约有17.52×104m3物源储量，动储量约6.03×104m3；中游启动区（海拔2335～2180 m），沟道长362 m，平均纵坡为375.23‰，概算大约有17.8×104m3物源储量，动储量约5.97×104m3。形成流通区单位面积上的物源量达98.11×104m3，该段沟道在下蚀和侧蚀作用下，大量残坡积物以垮塌和滑坡的形式持续补给泥石流，增大了泥石流的规模和破坏力。

2.2 泥石流暴发历史及危害

刘家沟泥石流威胁对象是沿沟两岸农田、沟口两岸的居民以及道路。自1997年起，每年雨季均会暴发不同规模的泥石流，受泥石流灾害的影响，目前沟口乡道已经改道。在近两年（2020年、2021年）发生的泥石流灾害中，泥石流漫出沟道，堵塞卡撒沟，形成2.0 m高的堰塞坝，造成上游淤积，回水淹没了沟道堆积区左侧平坝区。由于大颗粒物质残留在主河道内，很难被主河携带走，卡撒沟在此处过流能力严重减弱，沟道被压缩变形。参见图4。

图4 2021年6月泥石流淤积情况Fig.4 Debris flow event on June,2021

2.3 灾害链机制

根据刘家沟泥石流灾害暴发历史分析，其暴发时间与连续降水或短时强降雨具有明显一致性，表明降水是该泥石流发生的触发因素。

降雨及其形成的径流直接导致沟谷两岸的滑坡和坡面倾蚀在流水侧蚀作用下解体，松散土体被席卷进入狭陡的沟道，形成堵塞体，在上游暴雨洪峰流量作用下，堵塞体溃决形成泥石流，这种堵溃并非单一的，而是沿沟道发育，多次堵溃作用下，对泥石流峰值流量产生极大的放大作用［15］。泥石流进入主河，堵塞主河，形成堰塞坝，造成危害。在整个泥石流形成和发生过程中，呈现出的是一个降雨—侵蚀下切—崩滑—堵塞体—堵塞体溃决—席卷沟床物质—进一步侵蚀下切—进一步崩滑—进一步堵溃—增大的泥石流规模—堰塞湖—堰塞坝溃决—卡撒沟中、下游损害的成灾过程。

2.4 泥石流流量特征值计算

2.4.1 不同降雨频率下泥石流的流量

据《四川省中小流域暴雨洪水计算手册》所附暴雨量等值线图，由变异系数查皮尔逊Ⅲ型曲线，求得P＝5%设防频率下的1 h、6 h雨强值分别为33.4、47.1 mm，利用雨洪法计算泥石流的流量：

式中：Qc——频率为P的泥石流峰值流量，m3/s；Qp——频率为P的暴雨洪水流量，m3/s；Ф——泥石流泥沙修正系数，Ф＝（γc-γw）/（γH-γc）；γc——泥石流重度，kN/m3；γw——清水的重度，kN/m3；γH——泥石流中固体物质重度，kN/m3；Dc——泥石流堵塞系数。下同。

泥石流重度是影响泥石流峰值流量的关键因素，从泥沙携带能力考虑，重现期越长，泥石流沟的水力条件越好，携带大块石能力越强，泥石流重度越高［16］。泥石流的流量计算结果见表1。

表1 不同频率下泥石流峰值流量计算Table 1 Calculated peak volumes of debris flow at different frequencies

2.4.2 一次泥石流活动总量预测

泥石流历时T和泥石流流量Qc是决定泥石流一次活动总量的关键，按泥石流暴涨暴落的特点，将其过程概化“三角形”状［17］，通过断面一次泥石流总量Wc由下式计算：

一次冲出固体物质的总量Ws由下式计算：

据式（2）、（3），计算出该沟一次泥石流总量和相应的固体物质总量（见表2）。

表2 泥石流固体物质总量Table 2 Total solids in debris flow at different frequencies

2.4.3 形态调查法

根据2021年6月20日泥石流泥位、过流面积，采用形态调查法（表3）计算泥石流的流量明显大于根据雨洪法计算的泥石流峰值流量（表2），表明沟道内堵溃作用造成了泥石流流量沿程放大。

表3 形态调查法计算刘家沟泥石流流量Table 3 Calculated volumes of debris flow in Liujia gully

2.4.4 泥石流冲击力计算

刘家沟形成流通区平均纵坡为468.16‰，使得泥石流在运动过程中，沿程流速加大，冲击力增强，本文选取拟设工程部位为计算断面，对单块块石最大冲击力进行计算。计算公式如下，计算结果见表4。

表4 单块块石最大冲击力计算结果Table 4 Calculated maximum impact force of a single stone

式中：F——对墩的冲击力，kN；γ——动能折减系数，巨石冲击时取0.3；α——受力面与泥石撞击面接触角，取90°；W——巨块石重力，kN；C1、C2——巨石与建筑物的弹性变形系数，取C1+C2=0.0005；Vs——泥石流中块石移运速度，m/s。

3 泥石流堵河可能性判别

3.1 泥石流与主河交汇条件

刘家沟沟口于2021年修建了长50.0 m的排导槽，排导槽过流面积23.22 m2，根据雨洪法计算是可以满足20年一遇泥石流过流能力的，但由于排导槽与主河卡撒沟交汇角为87°，近似直角交汇，泥沙沉降速度加快［18］，较易形成堆积体，增加堵塞卡撒沟的可能性。

3.2 堵河可能性计算

根据泥石流沟口的断面宽度、流量及泥石流重度参数、主河在汇口处断面宽度、流量及河水重度参数，泥石流沟与主河交汇的夹角等，将主河与泥石流不同设计洪水频率以及两者不同组合条件下的相关参数带入泥石流堵河流量判别指标［19］公式进行计算：

式中：QM——主河单宽流量，m3/s；QB——支槽泥石流单宽流量，m3/s；θ——泥石流沟与主河交汇的夹角，（°）；γG——支槽泥石流重度，kN/m3；γM——主槽水流重度，kN/m3。

计算得CF＝-8.42，即会产生堵河。

徐永年等［20］通过分析堵河影响因素，选取主要因子，得到泥石流堵江的危险程度表达式：

式 中：γm、γc——主 河 和 支 沟 重 度，kN/m3；Qm、Qc——主河和支沟流量，m3/s；Jc、Jm——主河和支沟比降，%；θ——入汇角，（°）。

据金川县河流统计数据，卡撒沟长度23.5 km，落差2250 m，坡降95.25‰，流域面积246.2 km²，多年平均流量4.31 m3/s，由于沟口修建的排导槽使泥石流归流，目前可见主支沟夹角为87°，卡撒沟水流重度取10 kN/m3，刘家沟泥石流重度取18 kN/m3。计算K＝20.74，根据庄建琦等［19］将K=10作为泥石流堵江的临界指标，K大于或接近10后泥石流便会发生堵江。

3.3 堵河规模判别

根据正交入主河的三角形堰塞体堵河，所需最小泥石流总量，则堵沟所需堆积体的最大体积依以下经验公式估算：

式中：B、h——河面宽度、水宽，m；b——滑坡堆积体顺河顶宽，m；ɑ——沟口原堆积扇或岸坡的坡度，（°）。

卡撒沟在此段河宽15.0 m，按形成3.0 m高的堰塞坝，水深1.5 m，沟口堆积扇坡度65°计算堵河所需堆积体体积为2419.30 m3，大于20年一遇泥石流固体物质含量。

4 泥石流防治的特殊性

（1）形成流通区沟道两岸松散残坡积物层较厚，岸坡垮塌发育，未见明显基岩出露，中部纵坡达432.28‰，沟道横向平均宽度约6.0 m，最窄处仅3.0 m，沟底跌水。若在此段修建拦挡工程，回淤库容和反压稳固方量极为有限，坝肩、坝基稳定性堪忧，且沟道狭窄、纵坡较陡，设备转运以及布置的施工难度大。

（2）刘家沟泥石流在形成过程中，沿程铲刮沟床质堆积物、席卷坡脚残坡积或崩坡积等松散堆积物以及堵塞体溃决导致流速、流量在向下游运动的过程中都呈现出显著的沿程放大效应，用传统的雨洪法推求泥石流流量，与现场泥痕调查法得出的流量指标相比偏小，防治工程的设计参数不能满足实际流量的要求。

（3）出山口距离主河入口处有600 m，此段纵坡达111.45‰，沟道仅有1.0 m深，沟型不明显，泥石流暴发时极易漫出沟道，因此此段需修建归流工程，避免漫流与堵河成灾。

5 防治思路与措施

5.1 泥石流形成流通区设置潜坎群

防冲潜坎可以起到稳固岸坡，防治物源揭底启动，潜坎应急防冲以格宾石笼为宜。格宾石笼具有施工快速，就地取材，搬运轻便，少扰动的施工工艺。主体联结强，整体稳定，耐冲刷，成本较低，因其柔性强，可抵抗高强度压力，可承受大范围的变形而不破裂，可以依地势灵活布置，整体结构稳定，可有效提高地基承载力和冲刷淘蚀，且具有很强的抵御自然破坏及耐磨蚀和抗恶劣气候影响的能力，使用寿命长。

5.2 出山口设置拦挡坝

刘家沟泥石流为粘性泥石流，并含众多大石块，纵坡较陡，在出山口处流速较大的石块将产生巨大的冲击力，严重威胁下游建筑物及卡撒沟的行洪。基于此，经现场比选在出山口处岸坡覆盖层较薄，且沟道变缓，具有较好的拦挡条件处修建一座拦砂坝，以拦挡沟道中的巨石及较大固体物质，减小泥石流对下游的破坏［21-22］。

5.3 堆积区设置排导工程

针对堆积区沟道纵坡可以排泄泥石流，与出山口拦挡坝衔接做成八字堤导流，修建双边防护堤，对泥石流进行归流，在入主河处修建双边排导槽，与主河呈小角度交汇，最大限度利用主河的输移能力，对泥石流进行排导［22］。

5.4 建立监测预警系统

降雨是刘家沟泥石流暴发的激发因素，流域内局部降雨明显，且流域附近无雨量站。当监测到可能发生泥石流的临界雨量时，及时进行预警预报，减少灾害发生的可能性。另由于沟道内物源比较丰富，在降雨激发情况下，物源一旦启动，极易形成堵溃性泥石流，因此建议加强群测群防，在沟道上游、中游及下游建立雨量监测系统，尤其在沟道中游安装裂缝计，加强对滑坡物源的监测，发现滑坡物源大规模启动时，及时采用声光报警器进行预警［23］。

防治工程规划示意见图5。

图5 防治工程规划示意Fig.5 Mitigation planning for Liujia gully

6 结论与讨论

通过现场调查、理论计算与分析，探讨了刘家沟泥石流的成灾特征并提出可行的防治措施，得出结果如下：

（1）刘家沟泥石流属于高频泥石流，由于沟道狭陡，沟口流速、流量出现沿程放大，使用雨洪法计算数据难以满足工程设计要求，设计理念、方法与设防标准均应有所调整。

（2）刘家沟泥石流流域面积仅0.36 km2，平均纵坡为409.8‰，流域完整性系数为0.168，沟道横断面以典型V形，单位面积上物源量达98.11×104m3。沟口已建排导槽与主河呈大角度相交，加重了堵河可能性，泥石流堵河定量分析以及堵河所需体积计算均表明在发生20年一遇泥石流灾害情况下均会发生堵河。

（3）受限于较窄的沟道以及极陡的纵坡降，施工设备及材料转运、基础开挖难度较大，在施工阶段，需地质灾害专业人员全面科学评估泥石流灾害的链式效应，科学有效组织施工，避免次生灾害发生。

（4）需加强对沟道物源、沟域内降雨量及泥石流活动的动态监测，尤其是对松散物源稳定性和失稳情况的监测，建议安装专业仪器在汛期实行实时监测。