文章用于分析震后泥石流物源演化规律的数据包括等高距为5 m的地形图,通过地形图绘制高精度数字高程模型(DEM),通过将影像与DEM 进行叠加,使得解译结果更加准确。另外从Spot,Goole Earth上获取了2008年、2011年、2013年、2017年4期多源高分辨率光学遥感影像(表2),所获取影像云层活动较少,数据质量较高,能满足物源解译目的。

在遥感影像中,土体、植被、水体在红、绿、蓝3个波段上有相应反射值,其中土体反射值最大,因此全色影像中,土体颜色较明亮^[12]。因此对于地震造成松散堆积体的形状,运动轨迹、变化范围等信息都可以根据影像色调、纹理等提取^[13]。由于不同来源影像,其传感器特征不一致,导致影像数据特征,坐标系等不一致,因此在解译前需要对多源影像进行地理配准。同时由于其他物体遮盖、背景反差等原因导致有些物源难以识别,因此需要对图像做增色处理、图像锐化、波段变换等工作^[12]。如图2所示,图2A为未处理影像,物源单体形状、物源流动轨迹、以及基岩坡等信息不明确,经过图像处理以后结果如图2B所示,物源单体形状、物源流动轨迹、基岩坡等信息较为清晰。在此基础上,通过人机交互方式开展泥石流物源解译。本文主要想探讨震后泥石流物源演化规律和侵蚀方式变化,因此将物源分为沟道物源和坡面物源进行解译,并通过ArcGIS软件强大的空间统计分析功能,对物源数量、面积、物源单体长度等进行计算。

表2 遥感影像相关信息

图2 处理前后影像对比

为探究物源演化规律,文章选取物源数量、物源面积、物源密度、物源长轴长、物源形状系数几个参数作为物源演化量化指标(表3)。地震触发滑坡、崩塌等不良地质现象产生大量物源,在影像上每一处崩塌、滑坡形成物源都有一定形状、范围。震后随着地质环境恢复、降雨和重力作用、物源稳定性提高等必将导致物源数量、面积发生变化。坡面或沟道物源在重力和降雨作用下逐渐向地势低洼处运移,因此不同类型物源在运移过程中会逐渐相连接。因此物源长轴长和物源连接度体现了物源的运动特性,其中物源长轴长是沿解译物源单体形状中线量取的长度,可以描述物源运移强度大小; 物源连接度可以描述坡面物源转换为沟道物源程度,同时也可以说明物源运移通道连通完善状况^[14]。物源形状系数是参照流域水系理论定义,是指物源面积与长轴长平方比值,描述物源形状变化的重要参数,其值越大,说明物源形状越接近方形,反之则说明接近长条形^[15]。通过物源连接度和物源形状系数变化,可以揭示震区泥石流物源演化方式变化规律。通过以上指标选取,将物源分为沟道物源和坡面物源两类,在参考前人建立物源解译标志基础上^[16],利用ArcGIS软件根据物源单体在影像上纹理、色调等变化进行详细解译。并利用软件空间分析功能和统计功能获取所需基础数据。利用软件属性统计功能时,分别建立物源面积、物源连接、物源长轴长等字段。将不同类型物源相互连接时属性记为“1”,反之记为“2”; 物源长轴长是通过量测工具量测并记录; 物源面积通过软件自动获取。结合前期他人解译结果、野外考察数据、野外核查、影像特征等将属于基岩沟床,基岩坡面的解译结果进行剔除,由此统计和建立论文所需基础数据。

表3 物源演化量化指标

通过对研究区多期遥感影像物源解译,利用ArcGIS软件计算功能和统计功能计算出物源面积、物源密度、物源数量、连接度等,整理结果见表4—6。由解译结可知,震初大量崩塌、滑坡造成松散物源堆积在沟道或悬挂在岸坡上。震后2008年,影像共解译出物源1 310处,物源面积7.86 km²,物源密度为0.4。由图3可知,地震导致靠近沟道两侧崩塌、滑坡物源直接进入主沟或支沟,支沟物源经过地表径流冲刷或泥石流挟带作用与主沟相连结,有利于震区泥石流物源输出,经过统计可知,2008年,坡面物源连接度为0.34,沟道物源连接度为0.45。到2011年,坡面物源和沟道物源都有所增加,物源数量增长到1 582处,增长了21%,面积为6.67 km²,物源面积增加了15%,物源密度为0.34。这是由于地震导致山体松弛、地表破碎,松散物源充足,在降雨触发条件下,进一步诱发新的崩塌、滑坡等不良地质现象,同时随着时间发展,一些崩塌、滑坡物源逐渐开始恢复稳定,因此表现出物源数量增多,而面积在逐渐减少趋势。另外在降雨冲刷和重力作用下,物源逐渐向下游运移,经沟道径流冲刷或泥石流挟带作用排出沟口。该时期内坡面物源连接度增长到0.56,沟道物源连接度增长到0.83。物源演化至2013年,部分坡面物源、沟道物源逐渐消失,物源数量减少到1 203处,较上期影像解译结果减少了24%,物源面积减少到4.67 km²,减少了30%,物源密度由0.34减少到0.24,减少了29%。这是由于震后,随着时间发展,震区地质环境逐渐改善,植被恢复、泥石流物源趋于稳定等所致。到2017年,坡面物源和沟道物源持续减少,部分物源逐渐消失。统计结果显示,到2017年,物源数量减少到1 004处,减少了17%,物源面积为3.25 km²,减少了30%,物源密度为0.17,减少了29%。

表4 物源解译统计

表5 坡面物源解译统计

表6 沟道物源解译统计

图3 物源解译结果

汶川地震后,对于地震导致崩塌、滑坡、泥石流等地质灾害后效应问题一直备受广大专家学者关注。由于数据限制,对地质灾害演化趋势预测的问题一直难以解决。震后有专家学者结合台湾集集地震、日本关东地震演化规律指出汶川震后滑坡和泥石流活跃性保持高强度水平,且将持续10 a以上^[4],还有学者通过分析地震前后地质灾害数据,指出地震灾害后效应将持续20～25 a^[3]结 论。震区泥石流物源不同演化阶段对泥石流暴发频率、规模等影响也不相同,为探究震区泥石流物源演化阶段,本文在前人研究成果基础上(图4),结合物源连接度、物源密度、物源数量、物源长轴长4个指标,将物源演化阶段划分为幼年期、青年期、壮年期、老年期4个阶段。

幼年阶段:该阶段约为地震后至2009年,该时期特征主要为地震触发大量崩塌、滑坡等不良地质现象,形成大量松散堆积体堆积在沟道和坡面上(图3)^[3],同时由于地震导致山体松弛、破裂,出现较多震裂和松动的山体,在降雨入渗增大孔隙水压和重力作用下,破碎山体逐渐趋向于失稳。

青年阶段:该阶段表现为由于前期降雨入渗增大孔隙水压,使得崩塌、滑坡等不良地质现象异常活跃,快速增加至最大值(图4,5)。同时在降雨形成地表径流冲刷和重力作用下,物源逐渐向下游运移,部分物源经支沟排入主沟。因此该时期物源连接度、物源数量、物源长轴长不断增长,而随着部分物源逐渐恢复稳定及大量物源输出流域,物源密度在逐渐减少(图5)。随着物源数量、物源连接度增加,使得该时期泥石流进入高暴发期,泥石流异常活跃。据图4,图5演化特征可知,该时期大致为2009—2011年前后。

壮年期:随着前期坡面物源连接度增加,坡面物源大量输出补给沟道物源,以致于泥石流暴发频率高,泥石流松散物源被大量输出流域。另外随着震区地质环境改善、植被恢复、崩塌和滑坡等造成松散固体物质稳定性提高,物源数量、物源密度、物源面积、物源连通度将逐渐降低(图5),但其所处水平仍高于震前最大暴雨年(图4)。结合唐川、黄润秋等人研究成果^[3-4],该时期大致为2011—2025年前后。

老年期:该时期是指随着崩塌、滑坡、泥石流等不良地质现象恢复稳定,地质灾害活跃性、数量等都低于震前最大暴雨年水平(图4)。随着崩塌、滑坡等不良地质现象活跃性持续降低,泥石流物源连接度、物源数量、物源密度将不断减少,震区泥石流活跃性将逐渐降低。

图4 汶川震区地质灾害演化趋势[3]

地震形成大量松散堆积体,Tang C等^[5]将震区泥石流物源总结为悬挂物源和沟道物源两类。在强降雨条件坡面物源或产生新的滑坡物源将向下运移输入沟道成为沟道物源^[17],据统计分析,震区泥石流物源约有50%以上来源于坡面物源^[18]。震后一段时间内,坡面物源稳定性较差,且物源量充足,在降雨条件下多以坡面流、碎屑流等方式补给沟道物源^[19],沟道物源常由于消防管效应、坡体液态化、以及地表径汇流冲刷等形式启动形成泥石流^[18-19]。泥石流在运动过程中不断发生揭底侵蚀以及沟道两岸物源不断坍塌补给泥石流,输出流域出山口。由坡面物源和沟道物源解译结果表5,表6可知,该时期坡面物源和支沟物源数量和连接度不断在增加。据Zhang S等^[20]研究发现,随着震区泥石流物源内细颗粒物质被冲走,物源抗侵蚀能力提高,使得物源供给能力下降,泥石流类型将由土力类泥石流转换为水力类泥石流。由表5,表6可知,随着时间发展,坡面物源连接度快速降低,而沟道物源连接度则将长期处于较高水平。另外物源形状系数越小,说明物源形状呈长条形,图6显示坡面物源和沟道物源形状系数整体上呈下降趋势,说明物源形状上逐渐向长条形演变。这主要是由于随着后期坡面物源植被恢复,物源稳定性提高,坡面物源侵蚀方式转变为细沟、冲沟侵蚀,由此导致物源形状呈现长条形。而沟道物源侵蚀由于坡面物源补给量减少,沟道物源细粒物质被冲走,因此沟道物源多以径流下切侵蚀为主,沟道两岸物源较少出现崩塌或坍塌补给现象。

图5 物源演化趋势