1/24/2025, 7:52:33 AM 星期五
网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

四川省地形起伏度与斜坡地质灾害空间分布关系研究  PDF

  • 李强 1,2,3
  • 伍剑波 1,2,3
  • 孙东 1,2,3
  • 杨涛 1,2,3
  • 罗小惠 1,2,3
1. 四川省地质环境调查研究中心,四川 成都 610081; 2. 四川省华地建设工程有限责任公司,四川 成都 610081; 3. 四川省地质灾害防治工程技术研究中心,四川 成都 610081

中图分类号: P642.2

最近更新:2024-07-31

DOI:10.12143/j.ztgc.2024.04.016

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

地形起伏度与斜坡地质灾害的空间分布密切相关,但同时考虑地形起伏度最佳统计单元、灾害类型、灾害规模的研究偏少。以地貌单元跨度巨大的四川省为例,运用均值变点分析法确定地形起伏度最佳统计单元,采用频率比法、信息量法、确定性系数法等3种数据统计模型,研究地形起伏度对斜坡地质灾害空间分布的影响。研究表明:平均地形起伏度与统计窗口大小存在显著正相关关系,最佳统计单元的窗口数量为20×20、面积大小为3.24×106 m2;斜坡地质灾害主要分布于丘陵和小起伏度山地,数量总占比为68.9%,其中滑坡易发优势区间为50~350 m,崩塌易发优势区间为50~200 m;小起伏区域利于发育小型灾害,大起伏度区域更利于发育中型和大型灾害。

0 引 言

崩塌、滑坡是我国山区斜坡地质灾害的两种主要类型,可导致重大人员伤亡和财产损

1。斜坡地质灾害发育分布与地形地貌条件密切相关,尤其地形复杂山区地形作为重要的地质环境因子,直接影响斜坡表层的物质流动与能量转2,深刻影响着区域斜坡物质结构、空间格局,从而控制了崩塌、滑坡的空间分布。

地形起伏度作为反映区域地形高差变

3、地表切割程4和地貌特征的重要指5,已被学者广泛应用于地质灾害的规律认识和评价研究,主要体现在孕灾条件分6-12、敏感性分13-17、易发性评18-20等三个方面。苗朝6、艾国栋7、李乾坤8、聂智强9、徐志11等学者考虑了地形起伏度对地质灾害孕灾条件和成灾机制进行研究;郭芳芳13分析了鄂尔多斯及其周缘地区地形起伏度与滑坡灾害发育分布的相关性;宿星14研究了黄土地区地形起伏度对滑坡的影响;张明媚16以太原市西山地质块体为研究区,分析了斜坡地质灾害敏感性评价中地势起伏度提取最佳尺度。而在易发性评价领域,如黄发明18、叶振南19、田乃满20大量学者把地形起伏度作为区域地质灾害易发性评价的主要指标之一,并以此分析地质灾害发育分布与地形的关系。

地形起伏度是一定统计单元内的地形高差,其值会随着统计单元面积大小改变而改变,直接影响研究区域地形起伏度的提取精

21。因此,确定最佳统计单元大小是探讨起伏度与地质灾害相关性研究的重要前22,并可提升区域地形特征分析的科学准确性。上述大部分研究仅将地形起伏度作为一个指标参与地质灾害评价,而对如何选取最佳统计单元的探讨较少。有鉴于此,本文以四川省为研究区,采用均值变点分析22-24确定地形起伏度的最佳统计单元,基于频率比、信息量、确定性系数等3种统计模型,探讨地形起伏度对不同类型、不同规模斜坡地质灾害空间分布的影响,为区域易发性评价中地形起伏度的应用提供更合理的技术思路,具有重要的实际应用意义。

1 研究区概况

四川省位于中国西南部,地处长江上游,范围介于北纬26°03'58''~34°19'13'',东经97°21'58''~108°31'58''之间,面积约48.6万km2,辖21个市(州)、183个县(市、区)。四川省自西向东地跨青藏高原、横断山脉、四川盆地等三个大地貌单元,地势整体呈西北高东南低的态势,地形起伏非常之大,地貌以山地为主,山区面积占2/3以

25。川西主要为高原地貌和高山峡谷区,西南主要为次高山和中山峡谷区,川东为四川盆地和盆缘山地,川中四川盆地以平原、丘陵地貌为主,盆缘山地以中低山为26。由于地形地貌类型的控制,导致研究区气候具有显著地域性差异、气候垂直变化大。四川盆地及东部周围山地属中亚热带湿润气候区。西南山地属亚热带半湿润气候区,其河谷地区受焚风影响形成典型的干热河谷气候,山地形成显著的立体气候。川西北高山高原高寒气候区,该区海拔高差大,气候立体变化明显。

四川省大地构造单元格局以龙门山—盐源一线为界,西为松潘—甘孜地槽褶皱系,西南为三江褶皱系,北为秦岭褶皱系,东为相对稳定的扬子准地台

11。研究区活动断裂发育,地质构造运动活跃,影响最大的断裂为龙门山断裂、鲜水河断裂和安宁河断裂,形成四川著名的“Y”字型构26图1)。四川省内地层岩性差异明显,东部地台区以红层陆相碎屑岩沉积岩类为主,西部地槽区以海相碎屑岩沉积为主,间夹碳酸盐岩和火山岩、凝灰岩27

图1  研究区地理位置与斜坡地质灾害空间分布

Fig.1  Geographic location and the spatial distribution of slope geohazards in Sichuan Province

2 地形起伏度的最佳统计单元计算

2.1 数据来源与处理

本研究采用的DEM数据为SRTMDEM 90M,栅格单元大小为90 m×90 m,来源于中国科学院计算机网络信息中心地理空间数据云(https://www.gscloud.cn)。通过ArcGIS进行拼接投影等前期处理,并将原始数据坐标系调整为CGCS2000,高斯投影,6度带中央经线105°,利用最新的四川省行政范围矢量图对DEM数据进行裁剪,获得研究区相应的DEM数据。

研究区有记录以来的历史崩塌和滑坡地质灾害数据来源于地勘单位、网络文献、野外调查和遥感解译等多种途径,通过对数据进行清洗剔除数据重复、坐标异常、信息缺失严重等无效点,最终得到有效数据共74844个,其中滑坡56093个、崩塌18751个。

2.2 计算原理:均值变点分析法

2.2.1 地形起伏度计算

地形起伏度(Relief Amplitude)是一定范围内最大高程与最小高程的差值,用于描述地貌形态、划分地貌类型的重要定量指标,其计算表达式如

23

RA=Hmax-Hmin (1)

式中:RA——地形起伏度,m;HmaxHmin——分别为单位面积内最大、最小高程值,m。

地形起伏度采用ArcGIS软件的领域分析工具进行计算,鉴于DEM栅格为方形,本次计算采用矩形统计范围。单位窗格大小为DEM原始栅格大小90 m×90 m,依次对窗格数量为2×2、3×3、4×4、…55×55时的统计单元进行地形起伏度计算,建立起伏度数据表。

2.2.2 均值变点分析法

地形起伏度与统计窗口面积大小呈logarithmic曲线形态特征,依据地形起伏度的定义和地貌发育理论,该曲线上必定存在唯一的一个由陡变缓的拐点,该点对应的统计窗格值即为最佳统计单

10,而均值变点分析法能精准识别出该拐点。因此,本文采用均值变点分析法计算四川省地形起伏度的最佳统计单元,具体计算过程如1422

(1)计算各窗格下单位地势度T

Ti=tisi    (i=2,3,4,,55) (2)

式中:Ti——各分析窗格下的单位起伏度,m;ti——各分析窗格的平均起伏度,m;si——各分析窗格的面积,m2i——矩形邻域的窗格大小。

(2)对单位地势度T取对数lnT,得序列XX为{Xii=2,3,4,…,55}。

(3)计算序列X的算数平均值X¯和离差平方和S

X=i=2nXin-1       (n=55) (3)
S=i=2nXi-X¯      (n=55) (4)

(4)令i=2,3,4,···,55,依次对每个i将上述Xi样本分为两段:X2X3,…,XiXi+1Xi+2,…,X55,分别计算每段样本的算术平均值Xi1¯Xi2¯及两段样本的离差平方和之和Si

Si=t=1i-1(Xt-Xi1¯)2+t=in(Xt-Xi2¯)2 (n=55) (5)

(5)计算SSi的差值ΔS

ΔS=S-Si (i=2,3,4,…,54) (6)

2.3 最佳统计单元分析

按照前文均值变点分析法计算公式,各统计单元对应的参数如表1所示。表1显示平均地形起伏度与窗格大小存在明显的正相关,统计窗格数量越多平均地形起伏度越大,为此对统计单元面积大小和平均起伏度进行了拟合(图2)。指数方程与对数方程两种拟合方式显示,相关系数分别是0.9792和0.9696。从曲线形态可以明显看出,平均起伏度随着统计单元面积增长趋势由陡变缓,明显存在一个变化趋势转折拐点。均值变点分析法就是从统计学的角度定量判断出该拐点,根据分析原理其对应的窗格大小即为地形起伏度最佳统计单元。

表1  统计窗格数量与地形起伏度关系统计
Table 1  Topographic relief amplitude changing with statistical window number
窗格数量平均地形起伏度/m统计单元面积/m2Si/m2ΔS/m2窗格数量平均地形起伏度/m统计单元面积/m2Si/m2ΔS/m2
2×2 41.57 3.24×104 50.17 1.18 29×29 606.35 6.81×106 15.85 35.50
3×3 80.90 7.29×104 44.43 6.92 30×30 618.34 7.29×106 16.60 34.75
4×4 117.30 1.30×105 39.15 12.20 31×31 630.35 7.78×106 17.40 33.95
5×5 150.98 2.03×105 34.59 16.76 32×32 641.15 8.29×106 18.25 33.09
6×6 182.29 2.92×105 30.68 20.67 33×33 653.11 8.82×106 19.15 32.20
7×7 211.62 3.97×105 27.34 24.01 34×34 663.89 9.36×106 20.09 31.25
8×8 239.04 5.18×105 24.49 26.86 35×35 674.44 9.92×106 21.07 30.27
9×9 265.01 6.56×105 22.06 29.28 36×36 683.91 1.05×107 22.09 29.26
10×10 289.56 8.10×105 20.00 31.34 37×37 694.59 1.11×107 23.15 28.20
11×11 312.87 9.80×105 18.26 33.08 38×38 704.93 1.17×107 24.24 27.11
12×12 335.13 1.17×106 16.80 34.54 39×39 714.46 1.23×107 25.36 25.98
13×13 356.21 1.37×106 15.59 35.76 40×40 723.93 1.30×107 26.51 24.83
14×14 376.42 1.59×106 14.60 36.75 41×41 734.56 1.36×107 27.70 23.65
15×15 395.86 1.82×106 13.80 37.54 42×42 742.96 1.43×107 28.90 22.44
16×16 414.25 2.07×106 13.19 38.16 43×43 752.24 1.50×107 30.13 21.21
17×17 432.36 2.34×106 12.74 38.61 44×44 760.22 1.57×107 31.39 19.96
18×18 449.57 2.62×106 12.43 38.91 45×45 769.31 1.64×107 32.67 18.68
19×19 466.05 2.92×106 12.26 39.09 46×46 778.47 1.71×107 33.97 17.37
20×20 482.36 3.24×106 12.21 39.14 47×47 786.68 1.79×107 35.29 16.05
21×21 497.73 3.57×106 12.26 39.08 48×48 793.98 1.87×107 36.63 14.71
22×22 512.59 3.92×106 12.43 38.92 49×49 802.49 1.94×107 37.99 13.35
23×23 527.62 4.28×106 12.68 38.66 50×50 810.12 2.03×107 39.37 11.97
24×24 541.48 4.67×106 13.03 38.32 51×51 817.81 2.11×107 40.77 10.58
25×25 554.69 5.06×106 13.45 37.90 52×52 827.31 2.19×107 42.18 9.17
26×26 568.32 5.48×106 13.95 37.40 53×53 833.29 2.28×107 43.60 7.74
27×27 581.43 5.90×106 14.52 36.83 54×54 841.61 2.36×107 45.04 6.31
28×28 594.14 6.35×106 15.15 36.19 55×55 849.71 2.45×107 46.49 4.85

图2  统计单元面积与平均地形起伏度拟合曲线

Fig.2  Fitting curve of statistical unit area and average topographic relief amplitude

通过对统计单元大小、离差平方和之和Si、差值ΔS进行分析,结果如表1图3所示,SSi的差值ΔS变化曲线形态呈抛物线:窗格数量2~20段,ΔS随统计单元面积增大而增大;窗格数量20~55段,ΔS随统计单元面积增大而减小。窗格数量20×20对应的ΔS值为抛物线顶点,此点ΔS值最大(图3),即为logarithmic曲线的拐点,其对应的窗格大小即为研究区地形起伏度最佳统计单元。因此,本DEM数据源窗格数量20×20所对应的统计单元为研究区提取地形起伏度的最佳统计单元,即统计单元面积大小为3.24×106 m2

图3  ΔS值变化曲线

Fig.3  Change curve of ΔS

3 地形起伏度对斜坡地质灾害空间分布的影响分析

3.1 研究方法

统计模型是研究影响因子与地质灾害空间分布相关性的常用手

14,各种模型和方法得到了很好的应用和推广。本文选取常用的频率比法、信息量法和确定性系数法综合对比分析地形起伏度与四川省斜坡地质灾害发育的相关性。

3.1.1 频率比法(FR)

频率比法是计算影响因子在不同分级区间内地质灾害发生的概率,本文采用此方法分析地形起伏度分级之间与斜坡地质灾害空间分布的相关性。本文频率比计算原理是地形起伏度在某一分级区间内地质灾害发育的数量和研究区总地质灾害数量的比值以及与该分级下的面积和研究区总面积的比值之比。其频率比公式及计算过程如

18

FRi=Ni/NAi/Ai=1,2,3,…,m (7)

式中:FRi——频率比值;Ni——起伏度的第i个分级区间中发生的地质灾害数量;N——研究区的地质灾害总数量;Ai——起伏度的第i个分级区间所占的面积;A——研究区总面积。

频率比值越大,表明发生灾害事件的可能性越大。

3.1.2 信息量法(IV)

信息量法一般认为地质灾害发生的可能性与预测过程中所获取信息的数量和质量有关,可以用信息量来衡量。本文采用单因子的信息量计算模型,具体函数表达

19

Ii= log2Ni/NAi/Ai=1,2,3,…,m (8)

式中:Ii——信息量值;其余同上。

信息量越大,表明发生灾害事件的可能性越大。

3.1.3 确定性系数法(CF)

确定性系数法是一个概率函数,属于一种二元统计方法,假设将来发生地质灾害的条件和过去发生地质灾害的条件相同。函数表达式

17

CFi=PPi-PPsPPi1-PPs,PPiPPsPPi-PPsPPs1-PPi,PPi<PPs  i=1,2,3,…,m (9)

式中:PPi——地质灾害在地形起伏度分类i中发生的条件概率,即地形起伏度分类i中发育的地质灾害点数量与地形起伏度分类i面积的比值;PPs——地质灾害事件发生的先验概率,即整个研究区的地质灾害总数与研究区总面积的比值,是一个定值。

确定性系数CF的值域区间为[-1,1]。其中正值代表事件发生确定性的增长,即该地形起伏度区间易于发生地质灾害;负值代表事件发生的确定性降低,即该类地形起伏度不易于发生地质灾害。

3.2 地形起伏度空间分布

基于均值变点分析结果,以窗格数量20×20作为统计单元,计算出四川省地形起伏度范围为0~1966 m,平均值482.36 m,标准差316.37 m。采用等间距法以50 m为间隔将四川省地形起伏度分为25个等级(表2),其中50~100 m段分布面积最大,面积占比7.49%;其次是≤50 m,面积占比6.26%。四川省地形起伏,总体上呈西高东低形势,地质构造是其重要影响因素,如构造复杂的龙门山断裂带西部山地、鲜水河断裂南段区域起伏度最大,构造相对简单的成都平原起伏度最小(图4)。

表2  起伏度分级情况统计
Table 2  The statistics of the classification of topographic relief amplitude
起伏度分级代号起伏度值/m分级面积/104 km2灾点数量

分级点密度/

(个·km-2)

起伏度分级代号起伏度值/m分级面积/(个·km-2)

灾点

数量

分级点密度/(个·km-2)
R1 ≤50 3.04 3269 0.11 R14 650~700 2.26 1575 0.07
R2 50~100 3.64 9992 0.27 R15 700~750 2.03 1348 0.07
R3 100~150 2.45 7258 0.30 R16 750~800 1.84 1148 0.06
R4 150~200 2.11 5502 0.26 R17 800~850 1.64 943 0.06
R5 200~250 2.22 5541 0.25 R18 850~900 1.43 712 0.05
R6 250~300 2.45 5256 0.21 R19 900~950 1.21 717 0.06
R7 300~350 2.61 4646 0.18 R20 950~1000 0.97 488 0.05
R8 350~400 2.71 3642 0.13 R21 1000~1050 0.81 420 0.05
R9 400~450 2.73 3205 0.12 R22 1050~1100 0.58 286 0.05
R10 450~500 2.72 2804 0.10 R23 1100~1150 0.47 208 0.04
R11 500~550 2.59 2539 0.10 R24 1150~1200 0.34 173 0.05
R12 550~600 2.46 2106 0.09 R25 >1200 0.86 330 0.04
R13 600~650 2.42 1916 0.08

图4  地形起伏度与斜坡地质灾害空间分布关系

Fig.4  The spatial distribution map of topographic relief amplitude and slope geological hazards

3.3 地形起伏度对斜坡地质灾害空间分布的影响

鉴于每个地形起伏度分级区间面积大小不同,其灾害发育数量不能准确衡量地形起伏度对斜坡地质灾害空间分布的影响。本研究基于GIS平台,分别采用频率比法(FR)、信息量法(IV)、确定性系数法(CF)3种统计模型,分析斜坡地质灾害对地形起伏度的影响。

3.3.1 地貌类型与灾害空间分布特征

基于GIS软件空间分析功能,提取地形起伏度各等级中斜坡灾害数量、类型、规模等统计数据。参考文献[

7]结合研究区地形起伏度空间分布特点,将四川省划分为平原台地、丘陵、小起伏度山地、中起伏度山地、大起伏度山地等五种地貌类型。如图5所示,不同地貌类型发育的斜坡地质灾害在空间分布上呈现出明显差异性。平原台地,起伏度≤50 m,主要分布于成都平原及周边、川西高原地区,面积占比6.3%,斜坡地质灾害数量占比4.4%,点密度为0.11个/km2;丘陵,起伏度50~200 m,主要分布于四川盆地腹地地区,面积占比16.9%,斜坡地质灾害数量占比31.5%,0.29个/km2;小起伏度山地,起伏度200~500 m,主要分布于川西、川南、川东北山地,面积占比37.1%,斜坡地质灾害数量占比37.4%,0.18个/km2;中起伏度山地,起伏度500~1000 m,主要分布于川东北大巴山、川西山地支流河谷地区,面积占比38.8%,斜坡地质灾害数量占比23.6%,0.09个/km2;大起伏度山地,起伏度1000~1966 m,主要分布于龙门山断裂带西北部、金沙江、雅砻江中下游段、大渡河、岷江上游段两岸及附近区域,面积占比6.3%,斜坡地质灾害数量占比3.0%,0.07个/km2

图5  地貌类型与斜坡地质灾害空间分布关系

Fig.5  The spatial distribution map of geomorphic type and slope geological hazards

3.3.2 地形起伏度与灾害类型的空间分布关系

图6可以看出,50~100 m区间为转折区间,该区间的斜坡地质灾害、滑坡、崩塌发育数量均高于其他区间,其中滑坡占该类灾害的10.7%,崩塌占该类灾害的22.5%,崩塌滑坡占总斜坡地质灾害的13.6%。大部分灾害分布于起伏度500 m以内,其中滑坡的占比为73.6%,崩塌的占比为72.7%,斜坡地质灾害占比为73.4%。总体上,地形起伏度>100 m的区间,灾害发育数量与地形起伏度大小总体上呈负相关,即起伏度越大灾害发育数量越少,并且斜坡地质灾害、滑坡、崩塌三者的曲线形态基本一致。

图6  地形起伏度分区与斜坡地质灾害数量分布相关性统计

Fig.6  The statistical chart of the correlation between topographic relief amplitude zoning and the

distribution of slope geological hazards

灾害类型方面,FR值、I值和CF值(简称“三值”)的变化趋势几乎一致,峰值段均为小起伏段,不同的灾害类型有一定的差异性:斜坡地质灾害、滑坡的峰值段为50~350 m区间;崩塌为50~200 m。斜坡地质灾害“三值”的最大值区间均为100~150 m。滑坡在区间100~150、150~200、200~250 m的“三值”基本相近,高于其余区间。崩塌区间50~100 m的“三值”明显高于其他区间,其次是区间100~150 m。滑坡总体形态上为单峰,峰值后“三值”表现为地形起伏度越大值越小(图7)。

图7  地形起伏度分级的FR值、IV值、CF值对比

Fig.7  The comparison of FR, I and CF value for classification of topographic relief amplitude

3.3.3 地形起伏度与灾害规模的空间分布关系

对四川省崩塌滑坡灾害数据统计显示,灾害规模以小型为主,全区74844个斜坡地质灾害中小型、中型、大型及以上的数量分别为66024个、7825个、995个,占比分别为88.2%、10.5%、1.3%。由于小型灾害的占比具有绝对优势,小型灾害与斜坡地质灾害数量分布曲线形态基本一致(参见图6b),其中50~100 m区间分布数量最多。中型灾害的分布曲线近似圆弧形,两端的低、高起伏度区间分布的灾害数量少,200~800 m等中间区段分布的灾害数量相对较多。大型及以上灾害分布曲线呈波动形,小起伏度<200 m区间的分布数量较少,中起伏度500~1000区间分布的数量相对较多。因此,地形起伏度对不同规模斜坡地质灾害空间分布的影响程度不同。

灾害规模方面,三种统计模型中仅FR值的变化趋势具有明显的差异性。小型斜坡地质灾害FR值形态呈单峰,峰值段为50~300 m,最大值区间为100~150 m,150 m以后区间表现为地形起伏度越大FR值越小。中型斜坡地质灾害FR值没有明显的峰值段,总体趋势变化较缓,>700 m的区间属于相对高值段,最大值区间为1000~1050 m和1150~1200 m;大型及以上斜坡地质灾害FR值与地形起伏度具明显正相关关系,呈现出随地形起伏度值增大FR值由低到高,区间950~1200 m属于峰值段,其中1150~1200 m为FR最大值区间。

3.4 结果讨论

研究区滑坡灾害的易发优势区间为50~350 m,崩塌灾害易发优势区间为50~200 m,小型灾害的易发优势区间为50~300 m,中型灾害的易发优势区间为700 m以上区域,大型灾害的易发优势区间为950~1200 m。小起伏区域更易发生小型斜坡地质灾害,大起伏度区域易发生中型、大型斜坡地质灾害。其主要原因为地形起伏度仅是影响斜坡地质灾害发育的主要因素之一,研究区内的其他地质环境条件对斜坡地质灾害空间分布也具有重要控制作用。

不同规模崩滑灾害的地形起伏度易发优势区间的差异性,是由多种因素导致。四川省低地形起伏度主要分布于成都平原及周边和川东丘陵区,其大部分区域的地形起伏度<300 m(参见图4图5),据统计地形起伏度平均值为136 m,灾害数量占总数的54.9%,其中小型灾害约有39000处,约占全省灾害总数的52.4%;该区域人类工程活动相对最为强烈、人口密度大,切坡修路建房、乡镇开发等工程活动可能导致边坡失稳,且该区域斜坡高差小,地形切割相对小,灾害体长、宽等几何参数也相对较小,据此灾害体点规模以小型为主。四川省中等起伏度及以上区域主要分布于川西、西南、西北山区,与三大断裂带构造格局空间上高度重叠,该地区多为中高山峡谷区,地形切割深度大,构造活动性强烈,在内外动力因素叠加影响下,易发生规模相对较大的斜坡地质灾害。本研究结果与以往学者的相关研究成果基本一

2528

综上所述,斜坡地质灾害空间分布规律是多种因素耦合的结果,地形起伏度对斜坡地质灾害发育影响是内外地质作用、人工干扰因素综合作用的一个表征。

4 结论

(1)平均地形起伏度与统计单元的窗格数量存在显著的正相关,其值随统计窗格数量的增加而增加;均值变点分析方法计算出四川省地形起伏度最佳统计单元的窗格数量为20×20,该单元面积大小为3.24×106 m2

(2)四川省依据地形起伏度划分出平原台地、丘陵、小起伏度山地、中起伏度山地、大起伏度山地等5种地貌类型,其面积占比分别为6.3%、16.9%、37.1%、38.8%、6.3%,即地貌类型以起伏度200~1000 m的山地为主;斜坡地质灾害主要分布于地形起伏度<500 m的区域,其中丘陵、小起伏山地灾害数量占比分别为37.4%、31.5%。

(3)四川省斜坡地质灾害规模以小型为主、类型以滑坡为主,地形起伏度对不同类型、不同规模斜坡地质灾害空间分布的影响具有明显差异性。滑坡、崩塌的易发优势区间分别为50~350 m、50~200 m;小型、中型、大型易发优势区间分别为50~300 m、>700 m、950~1200 m。小起伏度易于发生小型灾害,大起伏度更易于发生大型灾害,这与前人研究结果基本一致。

(4)斜坡地质灾害空间分布规律是人类工程活动与地质环境因素耦合影响的反映,地形起伏度可作为区域斜坡地质灾害敏感性分析的重要因子,其研究结果可为确定斜坡地质灾害易发靶区提供判断依据。

参考文献(References)

1

李益敏谢亚亚蒋德明.怒江州斜坡地质灾害孕灾环境因素敏感性研究[J].水土保持研究2018255):300-305. [百度学术] 

LI YiminXIE YayaJIANG Deminget al. Study on sensitivity in disaster‑pregnant environmental factors of slope geological hazards in Nujiang prefecture[J]. Research of Soil and Water Conservation2018255):300-305. [百度学术] 

2

白希选严翼.武陵山区地形起伏度特征及其与土地利用的变化关系[J].武汉工程大学学报2020421):73-78. [百度学术] 

BAI XixuanYAN Yi. Characteristics of topographic relief amplitude and its relationship with land use change in Wuling mountainous area[J]. Journal of Wuhan Institute of Technology2020421):73-78. [百度学术] 

3

Willett S DSlingerland RHovius N. Uplift, shortening, and steady state topography in active mountain belts[J]. American Journal of Science20013014-5):455-485. [百度学术] 

4

王志恒胡卓玮赵文吉.应用累积和分析算法的地形起伏度最佳统计单元定[J].测绘科学2014396):59-64. [百度学术] 

WANG ZhihengHU ZhuoweiZHAO Wenjiet al. Extracting optimum statistical unit for relief degree of land surface with CUSUM algorithm[J]. Science of Surveying and Mapping2014396):59-64. [百度学术] 

5

李炳元潘保田韩嘉福.中国陆地基本地貌类型及其划分指标探讨[J].第四纪研究2008,(4):535-543. [百度学术] 

LI BingyuanPAN BaotianHAN Jiafu. Basic terrestrial geomorphological types in China and their circumscriptions[J]. Quaternary Sciences20084):535-543. [百度学术] 

6

苗朝张勇李乾坤.基于PFC2D的缓倾顺层滑移-溃曲型岩质滑坡破坏机制分析[J].钻探工程2023505):11-17. [百度学术] 

MIAO ZhaoZHANG YongLI Qiankunet al. The deformation mechanism analysis of gently inclined and consequent sliding‑buckling rock landslides based on PFC2D[J]. Drilling Engineering2023505):11-17. [百度学术] 

7

艾国栋李云青资丽君.天然和暴雨条件下的物质点滑坡变形机制分析[J].钻探工程2023505):23-28. [百度学术] 

AI GuodongLI YunqingZI Lijunet al. Analysis of deformation mechanism of landslide under natural and rainstorm conditions based on material point method[J]. Drilling Engineering2023505):23-28. [百度学术] 

8

李乾坤苗朝程英建.川东巴中市恩阳区近水平土质滑坡破坏模式及专业监测预警[J].钻探工程2023505):43-51. [百度学术] 

LI QiankunMIAO ZhaoCHENG Yingjianet al. Damage patterns of near‑horizontal soil landslides and professional monitoring and warning in Enyang District, Bazhong City, East Sichuan[J]. Drilling Engineering2023505):43-51. [百度学术] 

9

聂智强崔杰孔梁钰.江西省永新县地质灾害孕灾环境分析[J].钻探工程2023505):67-72. [百度学术] 

NIE ZhiqiangCUI JieKONG Liangyu. Analysis of geological disaster‑breeding environment in Yongxin County, Jiangxi Province[J]. Drilling Engineering2023505):67-72. [百度学术] 

10

伍剑波孙强张泰丽.地形起伏度与滑坡发育的相关性——以丽水市滑坡为例[J].华东地质2022432):235-244. [百度学术] 

WU JianboSUN QiangZHANG Tailiet al. Research for the correlation between relief amplitude and landslides: a case study of Lishui City[J]. East China Geology2022432):235-244. [百度学术] 

11

徐志文.四川省地质环境状况及地质灾害发育特征研究[J].地质与勘探2006424):97-102 [百度学术] 

XU Zhiwen. Proposing the frame and measures for prevention and cure system of geological disasters in Sichuan Province[J].Geology and Prospecting2006424):97-102 [百度学术] 

12

龚诚黄海陈龙.藏东地区冻错曲流域崩塌发育特征与链式成灾模式研究[J].钻探工程2023505):1-10. [百度学术] 

GONG ChengHUANG HaiCHEN Longet al. Characteristics of collapse development and geohazard chain model in the Dongcuoqu Basin, eastern Tibet[J]. Drilling Engineering2023505):1-10. [百度学术] 

13

郭芳芳杨农孟晖.地形起伏度和坡度在区域滑坡灾害评价中的应用[J].中国地质2008351):131-143 [百度学术] 

GUO FangfangYANG NongMENG Huiet al. Application of the relief amplitude and slope analysis to regional landslide hazard assessments[J]. Geology in China2008351):131-143 [百度学术] 

14

宿星魏万鸿郭万钦.基于SRTMDEM的地形起伏度对天水市黄土滑坡的影响分析[J].冰川冻土2017393):616-622. [百度学术] 

SU XingWEI WanhongGUO Wanqinet al. Analyzing the impact of relief amplitude to loess landslides based on SRTM DEM in Tianshui prefecture[J]. Journal of Glaciology and Geocryology2017393):616-622. [百度学术] 

15

张峰赵忠国谢大伟.基于确定系数法和信息量模型对新疆特克斯县滑坡敏感性评价的对比研究[J].新疆地质2019374):575-579. [百度学术] 

ZHANG FengZHAO ZhongguoXIE Daweiet al. Comparative study on landslide sensitivity assessment in Turks County of Xinjiang based on determined coefficient method and information quantity model[J]. Xinjiang Geology2019374):575-579. [百度学术] 

16

张明媚薛永安.斜坡地质灾害敏感性评价中地势起伏度提取最佳尺度研究[J].太原理工大学学报2020516):881-888. [百度学术] 

ZHANG MingmeiXUE Yong’an. Optimal scale for extracting relief amplitude in slope geological hazard sensitivity evaluation[J]. Journal of Taiyuan University of Technology2020516):881-888. [百度学术] 

17

薛永安王玉洁朱婧聪.县域国土空间斜坡地质灾害敏感性评价研究[J].自然灾害学报2022314):219-230. [百度学术] 

XUE Yong’anWANG YujieZHU Jingconget al. Study of slope geological hazard sensitivity evaluation in the county territory of territorial spatial[J]. Journal of Natural Disasters2022314):219-230. [百度学术] 

18

黄发明曹昱范宣梅.不同滑坡边界及其空间形状对滑坡易发性预测不确定性的影响规律[J].岩石力学与工程学报202140S2):3227-3240. [百度学术] 

HUANG FamingCAO YuFAN Xuanmei. et al. Effects of different landslide boundaries and their spatial shapes on the uncertainty of landslide susceptibility prediction[J]. Chinese Journal of Rock Mechanics and Engineering202140S2):3227-3240. [百度学术] 

19

叶振南田运涛陈宗良.西藏芒康县斜坡地质灾害空间分布特征与易发性区划[J].自然灾害学报2021303):199-208. [百度学术] 

YE ZhennanTIAN YuntaoCHEN Zonglianget al. Spatial distribution and susceptibility of geological disasters in Mangkang County[J]. Journal of Natural Disasters2021303):199-208. [百度学术] 

20

田乃满兰恒星伍宇明.人工神经网络和决策树模型在滑坡易发性分析中的性能对比[J].地球信息科学学报20202212):2304-2316. [百度学术] 

TIAN NaimanLAN HengxingWU Yuminget al. Performance comparison of BP artificial neural network and cart decision tree model in landslide susceptibility prediction[J]. Journal of Geo‑information Science20202212):2304-2316. [百度学术] 

21

钟静卢涛.中国西南地区地形起伏度的最佳分析尺度确定[J].水土保持通报2018381):175-181,186. [百度学术] 

ZHONG JingLU Tao. Optimal statistical unit for relief amplitude in southwestern China[J]. Bulletin of Soil and Water Conservation2018381):175-181,186. [百度学术] 

22

南希李爱农景金城.中国山地起伏度计算中地形自适应滑动窗口获取与验证[J].地理与地理信息科学2017334):34-39,2. [百度学术] 

NAN XiLI AinongJING Jincheng. Calculation and verification of topography adaptive slide windows for the relief amplitude solution in mountain areas of China[J]. Geography and Geo‑Information Science2017334):34-39,2. [百度学术] 

23

张学儒官冬杰牟凤云.基于ASTERGDEM数据的青藏高原东部山区地形起伏度分析[J].地理与地理信息科学2012283):11-14. [百度学术] 

ZHANG XueruGUAN DongjieMOU Fengyunet al. Analysis on the relief amplitude based on ASTER GDEM data in mountain area of eastern Tibetan Plateau[J]. Geography and Geo‍‑Information Science2012283):11-14. [百度学术] 

24

封志明李文君李鹏.青藏高原地形起伏度及其地理意义[J].地理学报2020757):1359-1372. [百度学术] 

FENG ZhimingLI WenjunLI Penget al. Relief degree of land surface and its geographical meanings in the Qinghai-Tibet Plateau, China[J]. Acta Geographica Sinica2020757):1359-1372. [百度学术] 

25

刘乾坤.基于SPARK的四川省滑坡灾害成因分析与临界阈值研究[D].成都电子科技大学2020. [百度学术] 

LIU Qiankun. Analysis of causes and critical threshold of landslide hazards in Sichuan Province based on SPARK[D]. ChengduUniversity of Electronic Science and Technology of China2020. [百度学术] 

26

杨红娟韦方强马振峰.四川省泥石流灾害的时空分布规律和降水特征[J].灾害学2017324):102-107. [百度学术] 

YANG HongjuanWEI FangqiangMA Zhenfenget al. Spatial‑temporal distribution of debris flows in Sichuan Province in China and the corresponding rainfall characteristics[J]. Journal of Catastrophology2017324):102-107. [百度学术] 

27

解征帆赵文华.四川省泥石流滑坡类型分布与防治对策[J].四川地质学报1992121):58-65 [百度学术] 

XIE ZhengfanZHAO Wenhua. Distribution and control measures of mud‑rock flows and landslides in Sichuan[J]. Acta Geologica Sichuan1992121):58-65 [百度学术] 

28

李云君刘志红吕远洋.四川省滑坡灾害气象预警模型建立与验证[J].地球信息科学学报2017197):941-949. [百度学术] 

LI YunjunLIU ZhihongLÜ Yuanyanget al. Establishment and validation of a meteorological warning model for landslide hazards in Sichuan Province[J]. Journal of Geo‑Information Science2017197):941-949. [百度学术]