0.   引言

红层是指岩石颜色以红色系为主，岩性主要为砾岩、砂岩、黏土质岩等中生代-新生代的湖相、河湖交替相或是山麓洪积相等陆相碎屑岩^[1]。四川是我国西南地区红层分布极为广泛的省份之一，主要分布在盆地内部及边缘地区^[2]，因其软硬相间的岩性组合，地质灾害频发、灾害严重，其中滑坡、崩塌是主要的灾害类型。在灾害发育分布规律方面，彭华等^[3]初步探讨了我国红层特点及分布规律。曹先康等^[4]将四川巴中的红层滑坡划分为浅表层土质滑坡、深厚层土质滑坡和缓倾顺层岩质滑坡等3类，研究了其分布规律和发育特征。在成灾机理与变形过程方面，马贤杰等^[5]基于滑坡的地层岩性和滑坡倾角，提出了红层地区滑坡的分类方案，总结了滑坡的发育模式。张涛等^[6]通过对川东地区典型红层缓倾岩质滑坡的调查研究，提出了该类滑坡的变形演化阶段及野外识别标志。王家柱等^[7]以马边县红层地区的黄子树滑坡为例，在工程地质调查的基础上研究了滑坡的发育特征和形成机制，并通过数值模拟还原了滑坡的运动特征。徐伟等^[8]以云南彝良县红层地区崩塌、滑坡为研究对象，总结了不同斜坡结构条件边坡的成灾模式和破坏过程。王军朝等^[9]在川东缓倾角红层地区划分出贯通型倾倒式、贯通型滑移式、非贯通型错断式、非贯通型拉裂式等4类崩塌，并对其岩性组合、斜坡结构特征、演变过程进行了系统总结。胡斌等^[10]以西南红层地区软硬互层型公路边坡为例，在边坡地质结构特征调查的基础上，研究了砂泥岩互层边坡崩塌形成机理、变形破坏模式及其稳定性。在灾害影响因素方面，郭永春等^[11]在分析我国红层分布特征的基础上，提出红层中存在砂岩、泥岩等软硬互层的结构，在工程活动等因素的影响下，沿软弱结构面易发生滑坡。李江等^[12]、张群等^[13]以四川省南江县50个红层缓倾角土质滑坡为研究对象，深入研究了该类型滑坡的发育环境、分布规律和影像因素，在降雨条件下，浸润峰深度是影响该类滑坡失稳的主要原因。连续性降雨可导致坡体地下水位逐步升高，稳定性逐步降低^[14]。近年来，红层软岩吸水膨胀引起的工程问题越来越多，受到学者的广泛关注^{[15 − 16]}。值得一提的是，红层后期受构造影响相对较弱，地层产状平缓，常形成典型的软硬互层缓倾角岩质边坡，当雨水入渗至层理面及优势结构面后，岩体抗剪强度大幅度下降，尤其是软弱夹层遇水膨胀、脱水收缩，与非红层地区的崩塌相比，其崩塌机理具有一定的特殊性，红层区缓倾角岩质崩塌受后缘拉裂缝、软硬岩体差异性风化形成凹岩腔的控制明显，降雨入渗和临空条件是诱发崩塌的主要因素。

2021年4月5日，四川省眉山市洪雅县柳江镇赵河村13组铁匠湾发生高位崩塌地质灾害，源区红层的特点是岩层产状较为平缓，崩塌发生时危岩体大面积垮塌，最远移动至2 km外的赵河河道，损毁洪瓦路约200 m、居民房屋5处，造成3人失踪、68人生产生活受到不同程度影响。本文在室内解译、野外调查的基础上，详细划分崩塌的分区特征，初步探讨崩落体的动力学特征和危岩体的变形历史，开展崩塌的成因机理研究，分析崩塌邻区潜在崩塌隐患的发育情况和发展趋势，以期为区域上类似灾害的防范提供资料支撑。

1.   地质环境条件

铁匠湾崩塌属于浅切割中山地貌，区域上该地貌主山脊近东西向展布，支脊近南北向延伸，沟深坡陡，峰峦叠嶂，瓦屋山一带多形成陡崖，植被发育，林地覆盖率75%～90%。年降雨量约1435.5 mm，年平均气温约16.6 °C，属于中亚热带湿润气候。研究区水系属于青衣江右岸支流花溪河的一部分，以小型溪沟（赵河）为主。

崩塌区及邻区出露夹关组（K₂ j）、蓬莱镇组（J₃p）、遂宁组（J₂sn）、沙溪庙组（J₂s）等中生代陆相沉积地层^[17]。其中，夹关组（K₂ j）以棕红色厚—巨厚层砂岩夹少量粉砂质泥岩为主，底部有数米厚的砾岩；蓬莱镇组（J₃p）以棕红、砖红色中厚层—厚层泥岩、粉砂质泥岩为主；遂宁组（J₂sn）以棕红、砖红色中厚层—厚层泥岩、粉砂质泥岩夹微量钙泥质细砂岩为主；沙溪庙组（J₂s）上部为灰绿、紫灰色中厚层—厚层岩屑砂岩、粉砂岩、粉砂质泥岩、薄层页岩，下部为紫红、绿、灰色薄—中厚层粉砂质泥岩、泥质粉砂岩、钙砂岩互层（图1）。

图  1  铁匠湾崩塌区域地质简图

Figure  1.  Regional geological map of the Tiejiangwan rockfall

下载: 全尺寸图片 幻灯片

崩塌区整体为单斜构造，以崩源区为中心，5 km范围内未见断层发育，断层对崩塌源区无直接影响。研究区及邻区近年来发生的地震活动有2013年4月20日雅安芦山地震，震中距离崩源区约71 km，崩塌发生前无可能影响坡体稳定性的相关地震报道，研究区区域构造背景相对稳定。

2.   数据来源及研究方法

本研究收集的数据主要包括1∶20万区域地质调查资料，崩塌前的多期亚米级光学遥感影像、数字高程模型（DEM）数据，崩塌后的无人机影像（正射、倾斜摄影），崩塌前后的降雨资料等。收集的数据基于不同的存储平台，数据格式、坐标系统多样化，本次采用2000国家大地坐标系、高斯-克吕格投影为匹配基准，进行相应的格式转化、坐标系投影变换，最终将区域地质、地质灾害数据统一到同一平台，建立研究数据库。本次首先通过对区域地质调查资料、多期光学遥感影像的分析，查明铁匠湾崩塌的地质环境条件、崩源区变形历史，其次在获取无人机影像（崩塌后）的基础上探讨崩塌分区特征和崩落体动力特性，最后开展崩塌成因及邻区崩塌灾害隐患的发展趋势分析。

3.   铁匠湾崩塌特征

3.1   崩塌分区特征

崩塌发生后，根据无人机光学影像，结合现场调查，铁匠湾崩塌光学影像色调整体呈暗紫红色，纹理较粗糙；崩源区可见砖红色带状图斑，走向近东西向，与周边墨绿色植被形成对比，为崩塌危岩体，岩层产状近于水平，可见明显的近期崩落面；崩塌堆积区呈带状沿斜坡分布，色调呈紫红色，纹理较粗糙，堆积区下部有雅赵公路和分散农户。崩塌纵长约1810 m，顶宽约181 m，主崩方向178°，平面面积约21.64×10⁴ m²。崩塌海拔高程介于1610（崩源区）～871 m（堆积区前缘），垂直高差约739 m。根据变形和运动特征，可以分为主崩塌区和崩塌影响区两大部分：主崩塌区分为崩源区1处、铲刮区1处、堆积区1处、流水二次搬运堆积区1处；崩塌影响区为崩塌发生时，受牵引作用的影响，局部发育拉裂缝，在崩源区后缘和铲刮堆积区两侧形成欠稳定块体，包括潜在崩源区1处、扰动变形区5处（表1、图2）。

表  1  铁匠湾崩塌光学遥感精细解译统计表

Table  1.  Statistical table of optical remote sensing interpretation in Tiejiangwan rockfall

名称	编号	面积/m2	厚度/m	估算体积/（104 m3）
潜在崩源区	Ⅰ	2311	110	25.42
崩源区	Ⅱ	9578	110	105.36
铲刮区	Ⅲ	30573	15	45.86
堆积区	Ⅳ	146296	10	146.30
流水二次搬运堆积区	Ⅴ	8852	3	2.66
扰动变形区	A	4836	15	7.25
	B	2874	10	2.87
	C	9190	5	4.60
	D	3104	5	1.55
	E	2004	5	1.00

下载: 导出CSV 

| 显示表格

图  2  铁匠湾崩塌分区特征及野外照片

Figure  2.  Subarea features and field photos of the Tiejiangwan rockfall

下载: 全尺寸图片 幻灯片

3.1.1   主崩塌区

（1）崩源区（Ⅱ）

崩塌源区位于铁匠湾斜坡上部的陡崖处，斜坡整体坡度约86°，顶部高程约1610～1600 m，下部剪出口高程约1500 m。崩源区岩层产状近于水平，产状1°∠5°，为岩层面缓倾坡内的岩质坡，崩落岩体平均长度约181 m，宽度约65 m，平均厚度（高度）约110 m，崩塌体积约105.36×10⁴ m³，见图2（a）—（c）。崩塌物质主要为夹关组棕红色砂岩、砾岩，岩体发育两组近于垂直的节理，产状分别为9°∠80°、97°∠73°。

（2）铲刮区（Ⅲ）

铲刮区主要分布在高程1500～1310 m，斜坡整体坡度约42°，当上部砂岩、砾岩岩体突然失稳向下崩落，以巨大的动能和势能作用于下部斜坡表面，被铲刮的物质包括残坡积物、老崩塌堆积及其下部基岩，铲刮作用使得下部基岩出露，坡体后退约15 m。铲刮区平均长度约287 m，宽度约238 m，平均铲刮厚度约15 m，铲刮体积约45.86×10⁴ m³，见图2（a）、（d）。

（3）堆积区（Ⅳ）

崩塌堆积区斜坡坡度相对较缓，堆积区顶部向下至洪瓦路区间段斜坡整体坡度约21°，洪瓦路向下至赵河区间段整体坡度约13°，具有上陡下缓的坡度特征。斜坡体表层以松散的残坡积层为主，下部基岩以易风化的泥岩、页岩为主，在洪瓦路两侧分布有分散农户，当途经铲刮区的崩落体继续高速向下运动，各崩落块体间相互碰撞解体，巨大的动能使得堆积区后部坡表的残坡积及上部植被、建筑物、公路等被高速铲刮和推挤，并将残坡积物推挤到堆积区中前部，根据洪瓦路铲刮错动出露的剖面推测被铲刮残坡积平均厚度约5 m。堆积区顶部向下至洪瓦路区间段位于堆积区中后部，分布有大量岩石碎块，成份主要为砂岩、砾岩，块径最大者可达20 m，以5～8 m居多，向下块径呈逐渐减小的趋势；洪瓦路向下至赵河区间段位于堆积区中前部，斜坡坡度进一步变缓，堆积物粒度亦相对较小，堆积体厚度逐渐加大。堆积区平均长度约1500 m，宽度约130 m，呈倒梯形，平均堆积厚度约10 m，堆积体积约146.30×10⁴ m³，见图2（a）、（e）。

（4）流水二次搬运堆积区（Ⅴ）

崩塌堆积区前部堵塞赵河河道，成分以泥沙为主，夹少量碎石块，在赵河流水二次搬运作用下，沿河道继续向下游运动，最远运动至双溪村，运动距离约580 m，按平均厚度3 m估算，流水二次搬运堆积体积约2.66×10⁴ m³，见图2（a）。

3.1.2   崩塌影响区

崩塌影响区为崩塌发生时，受牵引作用的影响，局部发育拉裂缝，在崩源区后缘和铲刮堆积区两侧形成欠稳定块体，分为潜在崩源区和扰动变形区两种类型。

（1）潜在崩源区（Ⅰ）

崩塌发生后，在崩源区后部形成环状危岩带，根据无人机光学影像，主崩区正后方发育有一系列走向约90°的拉裂缝，影像上呈凹槽状，岩性为夹关组棕红色砂岩、砾岩，潜在崩源区紧邻已崩落的崩源区，推测岩体仍发育两组近于垂直的节理。环状危岩带长度约180 m，宽度约15 m，平均厚度（高度）参照已崩落的崩源区厚度约110 m，崩塌体积约25.42×10⁴ m³，见图2（a）。

（2）扰动变形区（A—E）

铁匠湾崩塌主崩方向约178°，受南北走向近于平行排列的山脊等微地貌影响，崩落体运动方向大致与南北向微地貌走向一致，沿沟谷向下高速运动的崩落体铲刮、推挤沟道两侧岩土体，对岩土体造成扰动变形破坏，局部形成拉裂缝。扰动变形区沿原有山脊呈断续透镜状排列，铲刮区发育扰动变形区（A），变形区长度约148 m，宽度约64 m，根据邻区铲刮出露的剖面推测厚度约15 m，推测体积约7.25×10⁴ m³；堆积区发育扰动变形区4处，变形区（B）长度约120 m，宽度约33 m，根据邻区铲刮出露的剖面推测厚度约10 m，推测体积约2.87×10⁴ m³；变形区（C）长度约240 m，宽度约45 m，根据邻区铲刮出露的剖面推测厚度约5 m，推测体积约4.60×10⁴ m³；变形区（D）长度约122 m，宽度约30 m，根据邻区铲刮出露的剖面推测厚度约5 m，推测体积约1.55×10⁴ m³；变形区（E）长度约100 m，宽度约17 m，根据邻区铲刮出露的剖面推测厚度约5 m，推测体积约1.0×10⁴ m³，见图2（a）。

3.2   崩落体动力特征分析

铁匠湾崩塌垂直高差H=739 m，水平距离L=1853 m，崩落体的等价摩擦系数（崩塌碎屑流运动特性参数）μ=H/L=0.3988，远小于国际公认的岩质材料摩擦系数经验值0.6，说明铁匠湾崩塌主崩落体失稳后发生了远程运动，为高速远程崩塌碎屑流灾害。为探讨崩塌发生后所形成碎屑流的动力学特征，采用谢得格尔法对对崩落体沿途运动速度进行估算^[18]，计算公式为：v=$\sqrt{2g(H-\mu L)}$，其中，v为崩落体速度，单位m/s；g为重力加速度，取9.8 m/s²；H为崩落体自崩塌壁后缘顶点运动至计算点的垂直高差，单位m；L为崩落体自崩塌壁后缘顶点运动至计算点的水平距离，单位m；μ为崩落体的等价摩擦系数，取μ=H/L=0.3988计算。由此可计算出崩塌-碎屑流从失稳、高速运动到停止各阶段的运动速度（图3），大致在堆积区后缘，崩塌-碎屑流达最大运动速度约60 m/s。

图  3  剖面图和进度分布图

Figure  3.  Geological cross-section diagram of the rockfall and velocity distribution of the debris flow movement

下载: 全尺寸图片 幻灯片

3.3   崩源区危岩体变形历史与形成机理

根据2013—2021年光学遥感数据（图4），对崩源区危岩体宏观变形情况进行动态分析，2013年4月17日WorldView-2光学影像显示，崩源区右侧已开始出现小规模的崩落现象，在危岩带下部斜坡形成长约161 m，宽约16 m的带状堆积体，植被基本崩落体铲刮和覆盖，崩落岩块块径最大者约10 m，运动最远处距崩落区约230 m；到2018年11月—2019年2月，WorldView-2光学影像显示危岩体未出现大面积的失稳崩落迹象，但危岩体下部斜坡溜滑迹象的数量和规模呈逐渐增加的趋势，溜滑作用可在危岩体下部形成岩腔，为危岩体的失稳准备了条件；直至2021年4月5日，危岩体大面积失稳，引发了铁匠湾崩塌-碎屑流高位远程灾害。2013—2021年变形历史表明，铁匠湾崩塌变形特征可大致分为“坡面局部崩落—整体失稳崩塌—撞击铲刮坡面—碎屑流”四个发展阶段，崩塌形成机理上属于拉裂式崩塌，具有如下特征：在岩性方面，崩源区为砂岩、砾岩夹泥岩的软硬互层结构，岩性差异大；在岩体结构方面，崩源区发育两组近于垂直的风化裂隙和重力拉张裂隙，继承了地质历史时期的构造节理；在微地貌与受力状态方面，随着岩体差异性风化，砂岩、砾岩区较泥岩区抗风化能力强，形成突出的悬崖，在野外较易识别上凸下凹的陡崖地貌，应力状态上部受拉、下部受压，上部岩体在拉应力、水压力、冻胀作用、根劈作用下，近于垂直的构造裂隙不断向下扩展，在裂隙末端形成应力集中区，并逐渐与邻区裂隙贯通形成裂缝带，一旦拉应力超过岩石抗拉强度，上部岩石发生失稳，表现为向坡外的拉裂倾倒变形。

图  4  铁匠湾崩塌源区多期遥感影像对比图

Figure  4.  Multi-period remote sensing images comparison of the source area of the Tiejiangwan rockfall

下载: 全尺寸图片 幻灯片

4.   成因及趋势分析

4.1   崩源区岩体结构条件

铁匠湾崩塌崩源区危岩体位于夹关组底部，岩性为棕红色砂岩、砾岩，岩层产状1°∠5°，危岩体下部为蓬莱镇组棕红、砖红色泥岩，形成“上硬下软”的岩石组合，因差异性风化等因素，泥岩易于风化而形成岩腔，风化较慢的砂岩在上部形成陡崖或凸出的剖面形态，陡崖高度约110 m。从应力状态来看，该类型剖面往往上部受拉形成拉张裂隙，下部处于挤压状态，同时受地质时期区域构造运动的影响，砂岩发育两组近于直立的节理，产状分别为：9°∠80°（J₁），97°∠73°（J₂），节理与岩层面共同将砂岩切割成巨大块体，其中J₁节理与斜坡走向大体一致，且为陡倾节理，为岩体的失稳创造了条件。

4.2   气候条件

崩塌发生前一周内洪雅县柳江镇赵河村进入昼夜温差大、持续降雨的过程，本次收集了距离崩塌较近的峨眉山雨量站的气象数据，3月31日研究区开始降雨，一直持续到崩塌发生后，由于崩塌位于斜坡后缘，崩源区砂砾岩裂缝发育，下部泥岩为隔水层，有利于雨水的汇聚，逐渐降低了岩体的抗剪切强度，增加了斜坡岩土体的静水压力，增大了斜坡的不稳定性，加上4月3日以来，昼夜温差逐渐增大，在崩塌发生当天昼夜温差达12.6 °C，在水压、冻胀等作用下裂隙、节理在山体内部不断扩展、贯通，直至拉应力超过岩石的抗剪强度，危岩体失稳引发崩塌地质灾害（图5）。

图  5  崩塌前降雨过程和温差变化柱状图

Figure  5.  Histogram of rainfall process and temperature variation before rockfall

下载: 全尺寸图片 幻灯片

4.3   铁匠湾陡崖区崩塌发展趋势

根据铁匠湾崩塌的变形历史，坡面岩体局部崩落在崩塌发生前的数年内可观察到其变形特征，铁匠湾陡崖区地处高位，隐蔽性强、难以接近，遥感技术是进行该类崩塌早期识别与持续监测不可或缺的手段之一。铁匠湾崩塌源区所在铁匠湾陡崖走向约96°，陡崖沿走向延伸约9.6 km，铁匠湾崩塌以西约2.4 km，以东约7.2 km，陡崖所在区域具有极其相似的孕灾地质环境条件，为潜在崩塌危岩带。陡崖区后缘至坡脚的相对高差由东至西呈逐渐递减的趋势，东侧相对高差约200 m，西侧相对高差约100 m。陡崖区受地层控制作用明显，分布范围与夹关组棕红色厚—巨厚层状砂岩夹粉砂质泥岩的位置大体一致，地层变形迹象弱，层理产状近于水平，砂岩与粉砂质泥岩接触部位常形成阶梯状的地形，上部砂岩往往临空。受地质历史时期区域构造作用的影响，陡崖区普遍发育与铁匠湾崩塌源区产状大体一致的两组近于直立的优势结构面，切割水平层理将岩体分成菱形块体，见图6（a）—（c）。陡崖区基岩裸露，植被整体覆盖少，东侧植被较西侧多，以低矮灌木为主，陡崖区后缘为山脊，地势相对平坦，植被茂密，以高大乔木为主，植被根系往往沿近于直立的优势结构面生长，利于地表水下渗，促进了冻胀作用、根劈作用等外动力地质作用的发生，陡崖区坡脚处为地形由陡变缓的过渡部位，同时也是夹关组砂岩与蓬莱镇组泥岩的地质界线，界线以下斜坡坡度约40°～50°，残坡积层厚度大，植被茂密，以高大乔木为主。

图  6  铁匠湾陡崖区及潜在崩塌隐患分布图

注：a为陡崖区WorldView-2光学卫星影像；b—d为陡崖区野外照片；e—j为潜在崩塌隐患WorldView-2光学卫星影像。

Figure  6.  Distribution map of the steep cliff area and potential rockfalls in the Tiejiangwan area

下载: 全尺寸图片 幻灯片

本次基于光学遥感技术和野外调查，在铁匠湾陡崖区识别出潜在的崩塌隐患点6处，见图6（e）—（j），其中铁匠湾崩塌以西2处（编号：QB-1、QB-2），以东4处（编号：QB-3、QB-4、QB-5、QB-6），新识别出的6处崩塌隐患点目前处于坡面局部崩落变形阶段，通过WorldView-2光学影像均可识别出小规模的崩落现象，其中QB-1、QB-4崩落迹象最明显，以上6处崩塌隐患点在强降雨等诱发因素下，危岩体很可能失稳引发崩塌-碎屑流地质灾害。同时，铁匠湾崩塌发生后，受牵引作用的影响，崩塌壁后缘局部发育拉裂缝，形成有新的危岩体，存在二次崩塌的风险。建议采用无人机、机载LiDAR等技术手段开展铁匠湾陡崖区崩塌隐患早期识别与持续监测。

5.   结论

（1）铁匠湾崩塌可分为主崩塌区和崩塌影响区2个区域。其中，主崩塌区包括崩源区1处、铲刮区1处、堆积区1处、流水二次搬运堆积区1处；崩塌影响区包括潜在崩源区1处、扰动变形区5处。崩塌源区在2013年已表现出变形迹象，其变形特征可大致分为坡面局部崩落—整体失稳崩塌—撞击铲刮坡面—碎屑流四个发展阶段，碎屑流最大运动速度约60 m/s。

（2）崩塌源区具有“上硬下软”的岩石组合，受内外动力共同作用，内动力地质作用导致岩体发育两组近于垂直的优势结构面；外动力作用主要为降雨和温差，雨水沿近于垂直的裂隙入渗，产生静水压力，加之温度变化引起的冻胀作用，对崩塌的稳定性产生了一定影响，从而导致源区危岩体垮塌，巨大的冲击力作用于危岩体下方的老崩塌堆积体和基岩，引起崩塌-碎屑流灾害链。

（3）铁匠湾陡崖区地处高位、隐蔽性强、难以接近，基于光学遥感技术和野外调查，在铁匠湾陡崖区识别出潜在崩塌隐患点6处，同时铁匠湾崩塌存在二次崩塌的风险，在强降雨等诱发因素下，危岩体很可能失稳引发类似的崩塌-碎屑流地质灾害。