冯 卫,毕银强,唐亚明,张乐中,李政国
(1.长安大学 地质工程与测绘学院,陕西 西安 710064; 2.中国地质调查局西安地质调查中心,陕西 西安 710054;3.西北地质科技创新中心,陕西 西安 710054; 4.机械工业勘察设计研究院有限公司,陕西 西安 710043)
我国是地震频发国家,中国大陆分布有很多的活动断裂带与构造区,无论是板块间还是板块内活动构造都十分强烈,特别是活动断裂带或构造区附近最易于发生破坏性地震[1]。因此,对活动断裂带沿线的地质灾害分布规律和断层效应开展研究具有较强的现实意义。
近些年来,关于活动断裂沿线的地质灾害分布规律及断层效应研究成为了地震研究的热点之一。国内外学者利用统计方法对地震地质灾害的分布与地震震级、地震烈度、距震中或发震断层距离、坡度、岩石类型的关系进行了深入分析和研究,形成了许多地震诱发地质灾害的统计规律[2-12]。尤其是黄润秋等[13]对汶川地震触发地质灾害的断层效应开展了深入分析,张永双等[14]以汶川地震诱发的大型滑坡为例,系统分析了活动断裂带与大型滑坡之间的关系,许冲等[15]对汶川地震滑坡的地层岩性、坡度、滑向、高程、所在烈度区、上下盘位置和距发震断裂的距离等7个参数进行了统计分析,他们的研究成果为查明活动断裂带地震地质灾害形成机制和总结地质灾害发育规律起到了明显的指导作用。
礼县至罗家堡断裂带是青藏高原东缘中段的一条重要断裂带,相关学者认为该断裂为1654年罗家堡8级地震的发震断裂[16-17]。笔者通过资料搜集、遥感解译和野外调查,总结了该断裂带沿线的地质灾害分布规律,并基于GIS技术开展了地震地质灾害断层效应研究,以期为类似断裂带沿线确定地质灾害的多发区域提供参考。
礼县至罗家堡断裂带平面上主要由3条斜列的次级断裂段组成,总长约150km,走向60°。西段从宕昌东到礼县东,长约90 km;中段从礼县南到罗家堡东北,长约40 km;东段从天水镇西到平南镇东,长约20 km[17]。本文以该断裂带的核心区(中段区)为研究区,沿线出露的主要地层有泥盆系、石炭系、三叠系、第三系及第四系,礼县至罗家堡断裂带沿线岩土体结构破碎,在地震、降雨等因素作用下,区内滑坡、崩塌等地质灾害较发育。通过资料收集、遥感解译和野外调查,礼县至罗家堡断裂带两侧发育有崩塌、滑坡等地质灾害点702处(区内泥石流发育程度很低,此处不含泥石流),其中地震地质灾害点396处(占总数的56.4%),地质灾害发育密度很大(图1)。
图1 研究区地理位置图Fig.1 Geographical location map of the study area
将研究区地质灾害发生的年代由老到新划分为5期:(Ⅰ)晚更新世以来的地质灾害(13万年~1万年);(Ⅱ)全新世以来的地质灾害(1万年~1654年);(Ⅲ)罗家堡地震诱发的地质灾害(1654年);(Ⅳ)罗家堡地震以来的地质灾害(1654年~2000年);(Ⅴ)2000年至今的地质灾害。
研究表明,绝大部分地质灾害(占总数的96.3%)发生在罗家堡地震同期及之后至今的时间范围内(由于罗家堡地震发生后研究区内未再次发生破坏性地震,因此该时期之后发生的地质灾害主要为降雨、人工开挖和流水侵蚀所致),而发生在晚更新世和全新世的地质灾害数量很少(仅占总数的3.7%)(图2)。初步分析认为,一是在该地区同时期发生地质灾害的数量的确很少;二是许多古滑坡或老滑坡由于发育在黄土覆盖层内,随着水土侵蚀的加剧而日益消亡,因此年龄越老的滑坡相对数量越少。
图2 不同期次的地质灾害发育数量对比图Fig.2 Comparison graph of the number of geological disasters in different stages
(1)区域分布的不均衡性
首先,不同地貌单元地质灾害发育程度不同。从整体上讲,研究区地质灾害主要分布在西汉水和燕子河流域,即中东部黄土梁峁地区,而礼县城区西南部低山丘陵和缓坡丘陵地带地质灾害不发育。其次,同一地貌单元内地质灾害的发育也不均衡,这既包括灾害数量的不均衡也包括灾害类型的不均衡。数量方面,在黄土梁峁地带灾害数量自西向东逐渐增加,最终形成大部分地质灾害在研究区东北部集中分布的局面。类型方面,小型崩塌和滑坡集中分布在研究区西部,大中型崩塌和滑坡主要分布于研究区东部,形成了大中型地质灾害在黄土梁峁区东多西少的局面。这与黄土地层厚度和沟谷发育情况有关,研究区东部黄土地层厚度大于西部黄土地层厚度,且东部沟谷发育深度一般为50~150 m,西部沟谷深度一般为30~100 m,而厚度大的黄土地层和较深的沟壑都有利于较大型地质灾害的发育。此外,东部区域距离罗家堡地震震中更为接近。
(2)沿沟谷的成带性和流域的集中性
研究区内地质灾害呈现“线—带—片状”的分布特征,具体表现为沿306和208省道及县级公路呈线状密集分布,沿西汉水、燕子河等河流或沟谷以及礼县至罗家堡活动断裂带呈带状分布,沿永兴、城关、祁山、盐官等4个乡镇以及高家村、陆家湾、李集寨、磨石村等村落呈片状分布。
以滑坡等地质灾害的剪出口位置为准,将研究区地质灾害发育部位分为斜坡上部、斜坡中部和斜坡下部。虽然各部位发育的地质灾害数量相差不大,但剪出口位置在斜坡中部和上部的地质灾害多由地震诱发,数量共计420处,占总数的60%,剪出口位置处于斜坡下部的地质灾害共计282处,占总数的40%,多为降雨、人工开挖和流水侵蚀所致(图3)。
图3 不同发育部位的地质灾害发育数量对比图Fig.3 Comparison graph of the number of geological disasters in different development parts
研究区地震地质灾害点共计396处(占灾害总数的56.4%),分布面积较大(占灾害总面积的84.6%),本节对统计的396处地震地质灾害点开展相应的断层效应研究。
地震地质灾害的分布位置与距断裂带远近密切相关(图4和图5),39.1%的地震地质灾害发育于距断裂带2.0 km范围内,28.3%的地震地质灾害发育于距断裂带2.0~4.0 km范围内,13.1%的地震地质灾害发育于距断裂带4.0~6.0 km范围内,5.1%的地震地质灾害发育于距断裂带6.0~8.0 km范围内,4.5%的地震地质灾害发育于距断裂带8.0~10.0 km范围内,9.8%的地震地质灾害发育于距断裂带10 km以外区域。总之,随着距断裂带越来越远,地震地质灾害发育的数量和密度逐渐减小。
图4 礼县至罗家堡断裂带两侧缓冲区地震地质灾害点密度分布图Fig.4 Distribution map of earthquake geological disaster points in the buffer zone on both sides of the Lixian-Luojiabu fault zone
图5 地震地质灾害发育数量与距断裂带远近的关系对比图Fig.5 Comparison graph of relationship between the number of earthquake geological disasters and thedistance from the fault zone
Abrahamson等[18]研究了1994年美国加利福尼亚州Northridge地震的近场强震记录和其他逆断层型地震的强震记录,首先提出了“上下盘效应”概念,认为断层上盘的加速度峰值普遍高于下盘的加速度峰值。黄润秋等[19]在研究汶川地震地质灾害的分布规律时发现,不仅发震断层上盘较下盘地质灾害分布密度大,分布范围更广,而且上盘地质灾害的规模也远较下盘大。
对于礼县至罗家堡断裂而言,绝大多数区域的地震地质灾害分布基本符合上述规律,尤其是在永兴乡境内,地质灾害分布规律更为明显,其正断层的上盘(西汉水右岸)比下盘(西汉水左岸)更为发育(图6)。其上盘灾点密度可达4.0~5.5个/km2,而下盘灾点密度最大为2.0~3.0个/km2。
图6 永兴乡段断裂上下盘地震地质灾害点密度分布图Fig.6 Distribution map of earthquake geological disaster points on the upper and lower walls of the Yongxingxiang segment fault
将研究区内滑坡等地震地质灾害的滑动方向按照45°进行分带,统计每个滑向区域的滑坡数量与滑坡面积(图7和图8)。结果表明,区内滑坡的优势主滑方向为北西—南东(NW-SE)向315°~360°,近似垂直于发震断裂走向,只有个别滑坡滑动方向接近断裂走向。这表明,在强震作用下,起控制作用的是强震地震波在垂直断裂方向的传播,正是由于地震波在垂直断裂方向具有最大的强度,从而将坡体在这个方向上震裂、滑移,出现一系列中大型滑坡、崩塌。这同时也表明,地震波和地面运动的强度具有显著的方向性,在垂直断裂方向最强。
图7 不同滑向的地震滑坡发育数量对比图Fig.7 Comparison graph of the number of earthquake landslides in different sliding directions
图8 不同滑向的地震滑坡发育面积对比图Fig.8 Comparison graph of developed areas of earthquake landslides in different sliding directions
研究区地处地震活动强烈的南北地震带北段,历史上曾发生过多次强震,其中影响最大的地震为1654年罗家堡8级地震,其最大烈度达到Ⅺ度。本节围绕地震地质灾害点密度和面密度与地震烈度的关系开展分析。
3.4.1 地震地质灾害点密度与地震烈度关系
将地震地质灾害点与罗家堡8级地震烈度分区图叠加后可以看出,研究区主要位于Ⅷ~Ⅺ烈度区内(图9)。从数量上来说,共有83处地震地质灾害点(占总数的21.0%)落在VIII度区,共有174处地震地质灾害点(占总数的43.9%)落在IX度区,共有103处(占总数的26.0%)地震地质灾害点落在Ⅹ度区,共有36处地质灾害点(占总数的9.1%)落在XI度区(图10)。仅从数量来看,这并不符合地震烈度越大地质灾害发育数量越多的规律。为了进一步分析两者之间关系,本文引入地震地质灾害点密度的概念,即单位面积内地震地质灾害发育的数量,其量纲为处/km2。以地震烈度值为横轴,地震地质灾害点密度为纵轴,得到地震地质灾害点密度随地震烈度变化的关系曲线(图11)。该曲线斜率为正,表明随着地震烈度的增加,相应地震地质灾害点密度逐渐增加,尤其从VIII度区到IX度区,点密度从0.18处/km2增加到0.57处/km2,增加了3倍有余。但该曲线的递增速率并非直线上升,从IX度区到X度区以及X度区到XI度区,尽管灾点密度持续上升,但上升趋势明显变缓。这表明越靠近极震区,地震地质灾害点密度尽管有所增加,但增加速率明显减小。
图9 地震地质灾害点与罗家堡8级地震烈度分区关系图Fig.9 Map of relationship between the earthquake geological distaster points and the intensity of the Luojiabu Ms8 earthquake
图10 不同地震烈度区内地震地质灾害数量统计图Fig.10 Statistical graph of the number of earthquake geological disasters in different earthquake intensity areas
图11 地震地质灾害点密度与地震烈度区关系曲线图Fig.11 Graph of relationship between point density of earthquake geological disasters and earthquake intensity area
3.4.2 地震地质灾害面密度与地震烈度关系
将每处地震地质灾害的面积与罗家堡8级地震烈度分区图叠加后可以看出(图12),共有3.15 km2(占总面积的5.4%)落在了VIII度区,共有24.97 km2(占总面积的43.1%)落在了IX度区,共有20.40 km2(占总面积的35.2%)落在了Ⅹ度区,共有9.42 km2(占总面积的16.3%)落在了XI度区(图13)。仅从面积来看,这并不符合地震烈度越大地震地质灾害发育面积越大的规律。为了进一步分析两者之间关系,本文引入地震地质灾害面密度的概念,即单位面积内地震地质灾害发育的面积,它是一个无量纲的量。以烈度值为横轴,地震地质灾害面密度为纵轴,得到地震地质灾害面密度随地震烈度变化曲线(图14)。该曲线斜率为正,但它与图11中的地震地质灾害点密度随地震烈度变化曲线不同,其递增速率是明显上升的,表明随着距离极震区越来越近,地震烈度越来越大,尽管单位面积上发育的地震地质灾害点数有增长缓慢的趋势,但地震地质灾害面积增长速率却是直线上升的。这说明距离震中越近的区域,越易发生大型或特大型滑坡,如海头村特大型滑坡就分布在距罗家堡地震震中1 km范围内。
图12 地震地质灾害分布面积与罗家堡8级地震烈度分区关系图Fig.12 Map of relationship between the earthquake geological distaster distribution area and the intensity of the Luojiabu Ms8 earthquake
图13 不同地震烈度区内地震地质灾害面积统计图Fig.13 Statistical graph of the area of earthquake geological disasters in different earthquake intensity areas
图14 地震地质灾害面密度与地震烈度区关系曲线图Fig.14 Graph of the relationship between area density of earthquake geological disasters and earthquake intensity area
(1)礼县至罗家堡断裂带沿线的地质灾害在历史时间上绝大部分发生在1654年罗家堡地震同期及之后至今的时间范围内,在水平空间上具有区域分布的不均衡性、沿沟谷的成带性和流域的集中性等特点,在垂直空间上斜坡各部位发育的地质灾害数量相差不大,但大部分地震滑坡的剪出口位于斜坡中部及上部。
(2)礼县至罗家堡断裂带沿线地震地质灾害发育的断层效应作用较明显,距断裂带越近,地震地质灾害发育数量和密度越大;活动断裂的上盘比下盘地震地质灾害分布密度更大;地震地质灾害的优势主滑方向近似垂直于法阵断裂走向;越靠近极震区,地震地质灾害点密度尽管有所增加,但增加速率明显减小;随着离极震区越来越近,地震烈度越来越大,尽管单位面积发育的地震地质灾害点数有增长缓慢的趋势,但地震地质灾害面积增长速率却是直线上升的。
(3)研究成果既为深入分析礼县至罗家堡断裂带沿线地质灾害的发育特征、形成机理及预测评价提供了较有价值的信息,也为类似活动断裂沿线地质灾害多发区域的划定提供了参考。