Development Feature and Risk Assessment of Geological Hazards Caused by the Ms 6.2 Earthquake in Jishishan County, Linxia, Gansu Province
-
摘要:
2023年12月18日23时59分，甘肃省临夏州积石山县发生Ms6.2级地震，诱发了大量崩滑地质灾害，严重威胁人民生命财产安全。基于震区现场实地调查成果，开展地震诱发加剧地质灾害的发育特征、危险性评价研究，提出防灾减灾措施建议。根据实地调查统计，此次地震后新增地质灾害隐患点64处，加剧的地质灾害隐患点有63处。新增和变形加剧的隐患点以崩塌为主，滑坡次之；规模等级以小型为主，中型次之。崩滑地质灾害主要集中在黄土丘陵区，以黄土陡坎地带切坡建房、切坡修路为主。利用GIS技术的加权信息量法评价积石山县域地质灾害易发性。结果显示，高易发区、中易发区、低易发区、非易发区占比分别为5.45%、9.83%、32.70%和52.02%，其中高易发区主要分布在积石山东部黄土丘陵区的山梁地带。基于区域地质灾害易发性评价，开展地震活动断裂、地震动分布以及不同降雨工况条件（10年、20年、50年、100年一遇）下积石山县域地质灾害危险性评价。结果显示：100年一遇极高危险区较10年一遇极高危险区的增幅最大为18.26%，说明未来区内遭遇极端降雨会显著提高区内地质灾害危险程度。研究认为，积石山地震地质灾害后效应将增强，未来崩塌、滑坡发生频次升高，地质灾害易在降水、冻融条件下形成，需针对性地采取防控措施，有效降低其威胁程度。
Abstract:On December 18, 2023 at 23:59, an Ms6.2 earthquake occurred in Jishishan County, Linxia Prefecture, Gansu Province, triggering a large number of landslides and geological disasters, seriously threatening the safety of people's lives and property. This article is based on the results of on-site investigations in earthquake prone areas, conducting research on the development characteristics and risk assessment of geological disasters induced by earthquakes, and proposing suggestions for disaster prevention and reduction measures. According to on-site investigations and statistics, there are 64 new geological hazard points added after this earthquake. There are 63 potential geological hazards that have intensified. The main hidden danger points of newly added and intensified deformation are landslides, followed by landslides; The scale level is mainly small, followed by medium. Landslide geological disasters are mainly concentrated in loess hilly areas, with the main focus on building houses and roads by cutting slopes in steep loess slopes. The weighted information method based on GIS technology was used to evaluate the susceptibility of geological disasters in Jishishan County. The evaluation results showed that the proportions of high susceptibility areas, medium susceptibility areas, low susceptibility areas, and non susceptibility areas were 5.45%, 9.83%, 32.70%, and 52.02%, respectively. The high susceptibility areas were mainly distributed in the mountain ridge area of the loess hilly area in the eastern part of Jishi Mountain. On the basis of regional geological hazard susceptibility assessment, a geological hazard risk assessment of Jishishan County was conducted based on seismic activity faults, seismic motion distribution, and different rainfall conditions (10 year, 20 year, 50 year, 100 year return period). The evaluation results showed that the maximum increase in the 100 year return period extremely dangerous area compared to the 10 year return period extremely dangerous area was 18.26%, In the future, encountering extreme rainfall in the area will significantly increase the risk of geological disasters in the area. Research suggests that the post earthquake effects of geological disasters in Jishishan will be enhanced, and the frequency of landslides and collapses will increase in the future. Geological disasters are prone to form under precipitation and freeze-thaw conditions, and targeted prevention and control measures need to be taken to effectively reduce their threat level.
-
Keywords:
- earthquake /
- geological hazards /
- developmental feature /
- risk assessment /
- Jishishan in Gansu
-
雅鲁藏布江与尼洋河位于青藏高原东南部，盆地内山脉纵横起伏，形成大量冲沟、峡谷和河流。内部动力作用非常活跃，地壳中初始高压应力释放，盆地岩石结构松弛。崩塌、滑坡和泥石流等自然灾害频繁发生（苏立彬，2020；武辰爽，2021）。滑坡是自然和人类活动引起的对土壤的破坏（Taalab et al.，2018）。它是一种以大量岩石、碎屑或泥土向坡面移动为特征的自然灾害。无论是由自然还是人类活动造成的滑坡，每年都会造成重大的经济损失（Tien et al.，2018）。因此，利用高效稳定的滑坡灾害评估技术，针对滑坡易发区，快速准确地识别高易发区的灾害，预测滑坡灾害的发生，可以有效地提高灾害预测的效率，减少滑坡灾害造成的损失，为防灾减灾提供参考（张琪等，2023；周硼焜等，2023）。
滑坡易发性划区是通过滑坡发生后的影响因子属性来预测滑坡发生的概率，是滑坡预测的有效方法（沈玲玲等，2016；孟晓捷等，2022）。滑坡易发性评价通常采用传统的定性方法和定量方法（贾俊等，2023）。定性方法依赖于专家在历史资料和滑坡清单的经验和意见，如加权线性组合与层次分析法（Rehman et al.，2022），但计算结果受人为因素影响。定量方法包括数据模型和确定性模型。确定性模型可以提供精确的分析结果，但需要大量的数据，尤其是在大尺度地区实践中难以获得（杨创奇等，2022）。近年来，包括机器学习和统计学的数据驱动模型在地质灾害研究方面取得了显著进展，如证据权模型（WoE）（Batar et al.，2021）、频率比（FR）（ Khan et al.，2019）和确定性系数法（CF） （乔德京等，2020）等。这些算法计算简便，甚至在一些大型区域也能适用，但是过分依赖样本质量且无法有效处理复杂的滑坡及其影响因子之间的关系。机器学习中的随机森林（Arabameri et al.，2019）、决策树（Hong et al.，2018）、BP神经网络（康孟羽等，2022；张林梵等，2022）、梯度提升等也被广泛地运用在滑坡识别中（张文龙等，2023），较好地解决了非线性关系表达的问题，提高了滑坡识别的精度。然而，这些模型通常依赖于单一的学习器，滑坡易发性所涉及的影响因子众多，通常很难获得理想的预测结果，容易发生过拟合现象。因此，笔者利用集成学习将多个单学习器组合起来进行区域滑坡易发性评估，以比较其与传统方法更具有优越性和高效性。
近年来，大量基于机器学习的方法被成功应用于地质灾害研究，而较新的梯度提升（Boosting）模型，包括XGBoost和LightGBM模型，在滑坡易发性方面很少被研究与比较，且不平衡类分布可能会影响特征选择的假设。在此基础上，笔者以雅鲁藏布江与尼洋河两岸为例，首次引入了基于基尼系数的加权随机森林作为特征选择过程过程，并与基于Boosting算法的XGBoost和LightGBM模型对研究区滑坡易发性进行分析和比较。
1. 研究区与数据
1.1 研究区
笔者选取雅鲁藏布江下游与尼洋河两岸为研究对象（图1）。研究区位于西藏自治区林芝市西部，E 92°09′～95°51′，N 27°55′～30°36′，总面积约为68 000 km2，包括工布江达县、波密县、米林县、朗县、墨脱县。研究区属于典型的高原丘陵、高山峡谷地貌，是世界陆地垂直地貌落差最大的地带，区内地形起伏大，呈现北高南低走势，山脉多为东西走向，绝大多数为高海拔大起伏山地，其次为高海拔极大起伏山地与中高海拔极大起伏山地，最高海拔7 782 m，地处米林县与墨脱县的交界地带。研究区位于高原温带湿润半湿润季风区气候带寒带跨越到热带。地区水汽含量高,雨季开始得早，结束晚，持续时间长，年平均降水量约为650 mm，年平均气温为9.1 ℃。研究区内有日土–青丁断裂、达机翁–朗县断裂、贾桑断裂、札达–邛多断裂等断裂带，主要出露底层有盆地相上三叠统的砂岩、夹板岩、火山岩以及海相下—中三叠统的千枚岩、砂岩、含砾状灰岩等。由于高降雨量以及土壤和板块内动力活跃，该区域极易发生滑坡。
![]()
图 1 研究区地理位置及滑坡分布Figure 1. Geographical location and landslide distribution of the study area1.2 数据来源与处理
主要数据来源包括：①地理空间数据云的ASTER GDEM 30 m分辨率数字高程数据，基于ArcGIS软件对坡度信息进行了提取。②1∶5万地质图来源于中国地质调查局，用来提取地层岩性性质。③Landsat8影像来源于地理国情普查，用于土地利用数据的提取。④滑坡数据出自中国科学院资源环境科学数据中心。⑤断层带从地震活动断层探察数据中心获取。
笔者在已有的研究方法上将30 m×30 m栅格大小设定为基础的评价单元（Tanyas et al.，2019），研究区域划分为123 156 296个网格。同时为了解决样本不均衡问题，笔者采用下采样方式从非滑坡区选取等量滑坡点组成188个样本点（Polykretis et al.，2018），滑坡单元设为1，非滑坡单元设为0，从中随机抽取70%（131）数据作为训练样本，剩余30%（57）作为测试样本。滑坡点具体流程见图2。
2. 评价因子选取与独立性检验
2.1 评价因子选取
已有对雅鲁藏布江流域的研究结果和现场勘查表明：河水对河谷的不断侵蚀作用加上高海拔高寒区冻融加剧滑坡区岩石的风化，使得雅鲁藏布江流域极易孕育滑坡（赵永辉，2019）；地层岩性是滑坡产生的重要因素（赵永辉，2021）；坡度为滑坡发生的主控因素（王瑞琪等，2019）。再根据对研究区的地质灾害形成条件与地质环境背景研究分析，选取高程、坡度、断裂带与断层、河流、道路、地层岩性、土地利用7个评价因子。利用ArcGIS软件，将高程、坡度、地层岩性、土地利用4个连续型因子结合分布规范，采用自然间断法将研究区分为5个等级（图3a～图3d），对于离散型因子例如断裂带与断层、河流、道路利用多环缓冲区工具建立0～200、200～400、400～600、600～800、＞800 m共5个等级范围（图3e～图3g）。
统计各评价因子分级范围内滑坡点数量并绘制蔟类柱状图（图4）。结果表明：当高程处在32～1 544 m时，滑坡发生的最多，占总数的30.9%，其次是出现在2 722～3 752 m范围内。其原因是在海拔低于1 544 m时，开挖坡脚等人类活动频繁，随着海拔的提升，坡度增大，加剧了滑坡的发生；随着坡度上升，滑坡数也增加，直到坡度上升达到阈值40°,发生灾害的概率降低，由原来的41.5%逐渐降低到16.0%；当地层岩性为雅鲁藏布江带闪片岩时，相比其他岩性，滑坡发生最频繁；草地土壤侵蚀严重，是浅层滑坡的重要原因。本研究中大量滑坡点分布在坡度为10°～20°的草地上；断裂带与断层会降低岩层的强度和完整性，是滑坡易发性增大的关键，在距断层带200 m以内容易发生滑坡，滑坡点占总数的41.5%，离断裂带与断层越远滑坡灾害越少；河岸受水流不断冲刷，土石在地下水及重力作用下越发失稳，因此越靠近河流越容易发生滑坡，滑坡在距河流200 m以内，发生次数最多，达到40.4%；修建铁路、公路时因大力爆破、强行开挖，常使坡体下部失去支撑而发生下滑，距离道路200 m以内的滑坡数占了总数一半以上达到52.1%，距离道路越远，滑坡活动减少。文中结论与相关研究均吻合（Kouhartsiouk et al.，2021；Zweifel et al.，2021）。
![]()
图 4 各评价因子与灾害点的关系Figure 4. Relationship between assessment factors and disaster points2.2 评价因子独立性检验
为了研究各评价因素的相对独立性以及评价模型的准确性和可靠性，笔者采用皮尔逊相关系数计算影响评价因子的相关性。皮尔逊相关系数是用于度量两个变量之间的线性关系，利用两个变量间的协方差和变量的标准差进行计算而来（Lee et al.，2020）。
$$ {\rho _{X,Y}} = \frac{{ \displaystyle \sum {XY} - \frac{{ \displaystyle \sum X \displaystyle \sum Y }}{N}}}{{\sqrt {\left( { \displaystyle \sum {{X^2} - \frac{{{{\left( { \displaystyle \sum X } \right)}^2}}}{N}} } \right)\left( { \displaystyle \sum {{Y^2} - \frac{{{{\left( { \displaystyle \sum Y } \right)}^2}}}{N}} } \right)} }} $$ (1) 式中：X，Y表示变量，N表示取值个数。
变量间呈现极弱相关时，相关系数为0.0～0.2；0.2～0.4表示变量之间弱相关性。将评价因子的7个属性值代入式（1）计算，结果见表1，发现相关性最高为坡度与道路（R=0.349 3），其他变量间相关关系均小于0.4。总体而言，变量的共线性不强。
表 1 因子间皮尔逊相关系数表Table 1. Pearson correlation coefficient between factors因子 高程 道路 河流 坡度 断裂带与断层 地层岩性 土地利用类型 高程 1.000 0 −0.162 4 0.155 4 −0.170 8 0.231 7 −0.256 4 −0.029 8 道路 −0.162 4 1.000 0 0.140 5 0.349 3 −0.207 6 −0.093 0 0.002 5 河流 0.155 4 0.140 5 1.000 0 0.126 9 −0.067 2 0.301 1 0.012 2 坡度 −0.170 8 0.349 3 0.126 9 1.000 0 −0.237 1 −0.051 0 −0.064 9 断裂带与断层 0.231 7 −0.207 6 −0.067 2 −0.237 1 1.000 0 −0.196 0 −0.265 4 地层岩性 −0.256 4 −0.093 0 0.301 1 −0.051 0 −0.196 0 1.000 0 0.072 5 土地利用类型 −0.029 8 0.002 5 0.012 2 −0.064 9 −0.265 4 0.072 5 1.000 0 3. 雅鲁藏布江与尼洋河两岸滑坡易发性评价
3.1 基于Gini–RF的滑坡易发性评价
随机森林（Random Forest）是一种基于决策树模型的Bagging（Bootstrap AGgregation）的优化版，由于其具有对特征鲁棒性强、适用于高维稠密性数据、并行集成、对不平衡的数据集可自动调整误差、微调超参数等优势，可以获得准确结果，常被用于各种分类和回归任务（Alsahaf et al.，2018）。它的基本单元是决策树，但其本质是集成学习方法，是机器学习的一个分支，其核心思想始终为Bagging。然而，已经做了一些特有的改进，随机森林使用CART决策树作为基学习器。
基于Gini系数的随机森林建立在许多决策树上并支持各种特征权重度量。其中之一为特征与不平衡数据输出的相关性，一旦分类器测量了Gini系数，这种特征选择技术就在 RF 中采用了权重调整技术。Gini指数具有在特定节点中划分二进制类的能力（Disha et al.，2022）。对于具有两个以上不同值的属性，考虑属性子集，通过调整不平衡类分布的随机森林算法中的权重，使用Gini系数标准来分裂树，计算特征重要性得分。GI值越高，特征对模型预测的平均贡献越大，模型的解释能力越好，所有GI特性之和为1。
$$ G I_{\mathrm{m}}=\sum\nolimits_{k=1}^{|K|} \sum\nolimits_{k^{\prime} \neq \mathrm{k}} p_{m k} p_{m k^{\prime}}=1-\sum\nolimits_{k=1}^{|K|} p_{m k}^2 $$ (2) $$V I M_{i j}^{(\text {Gini })}=\sum\nolimits_{m \in M} V I M_{j m}^{(\text {Gini })}$$ (3) $$ V I M_j=\frac{V I M_j}{\displaystyle \sum\nolimits_{i=1}^c V I M_i}$$ (4) 公式（2）：
$ GI_{\text{m}}^{} $ 为基尼指数，K代表k个类别，$ {p_{mk}} $ 表示节点m中k的比例；公式（3）：$ VIM_{ij}^{(Gini)} $ 表示特征i在第j颗树的权重；公式（4）表示对所求出的所有重要度得分进行归一化处理。笔者把94个滑坡点记为‘1’，等量非滑坡点记为‘0’，将7个评价指标因子的属性提取至训练集，构造随机森林二分类模型，并从sklearn库中调用Random Forest Classifier方法，将训练集代入RF模型进行训练。同时为了确保结果的可靠性和准确性，在原本的参数设定基础上，采用贝叶斯优化算法搜索最优参数值。优化结果中，当每次迭代完成后更新权重时的步长取0.1，max_depth取4，num round取30时，效果最佳。用测试集对RF模型进行预测，结合公式（3），将得到各评价因子的权重归一化后导入ArcGIS中的栅格计算器生成滑坡易发性图，采用自然间断法将分区图划分为极高、高、中、低、极低5个等级（图5），易发性越高代表越容易发生滑坡。
3.2 XGBoost易发性评价
XGBoost是一种基于决策树模型和梯度提升的集成机器学习算法，为了控制模型的复杂度，它将正则化项添加到损失函数中，正则项包括每个叶子节点权重的平方和与节点个数。XGBoost处理缺失值并通过学习模型选取缺失值最佳的默认分割方向（Inan et al.，2021）。
描述的数据在预处理过程之后，基于Python3.6与R语言，采用Scikit-learn构建XGBoost多分裂滑坡易发性模型（Alsahaf et al.，2018）。同时为了在独立的验证数据集上对子序列进行测试降低偶然性，选取最优子树，通过贝叶斯算法优化，利用五折交叉验证获得每个模型评价度量的平均值，所有测试集的平均指标被认为是最终结果。将预测结果导入ArcGIS绘制滑坡易发性图（图6）。样本集在所选参数值上的交叉验证准确度结果显示：当进行第5次五折交叉后，训练集和测试集的AUC值达到最大值并趋于稳定（图7）。
3.3 LightGBM易发性评价
Light Gradient Boosting Machine （LightGBM）是一种高性能、开源、快速的分类、回归、排名的方法，同时也是基于决策树算法的梯度提升算法。LightGBM 采用直方图算法来降低内存消耗，使数据分割更简单，将浮点的连续特征离散化为式子中的k个离散值，构造一个宽度为 k的直方图，将数据进行遍历训练，计算直方图中每个离散值的累积统计信息，在特征选择中，只要根据直方图离散值搜索最佳的分割点即可（Zeng et al.，2019）。
在4.2使用方法基础上，将研究区的123 156 296个栅格提取各评价因子的属性值到点，生成123 156 296×7的表格，导入训练好的机器学习模型中，预测每个栅格发生滑坡的概率，利用点转栅格工具将所有的点生成栅格数据，再用自然间断法将研究区的滑坡易发区分为极高、高、中、低、极低5个类别（图8）。图9为LightGBM的学习曲线。
4. 滑坡易发性评价结果验证
4.1 易发性分区结果与对比
基于ArcGIS，分别统计3种不同机器学习模型在每个易发性分区的栅格个数与滑坡点个数（表2），3种模型的滑坡易发性结果呈现出一定的差异，但整体趋同。Gini–RF、XGBoost 和LightGBM模型均在极低类别中的百分比值最高。对于Gini–RF模型，从极高到极低易发性的面积比分别为11.99%、12.63%、19.58%、26.77%和29.03%。XGBoost模型的极高、高、中、低和极低易发性区域分别占12.05%、12.50%、19.62%、26.78%和29.05%。对于LightGBM模型，极低、低、中、高和极高易发性区域分别占12.14%、12.41%、19.43%、26.47%和29.55%。根据滑坡位置的分布可以看出，大多数历史滑坡记录位于高易发性地区，正如Gini–RF、XGBoost和 LightGBM模型所预测的那样。LightGBM模型的性能最高，其次为XGBoost与Gini–RF。
表 2 机器学习模型易发性分区对比Table 2. Comparison of machine learning model vulnerability zones类别 机器学习模型 Gini–RF XGBoost LightGBM 栅格
个数栅格
占比（%）滑坡
点个数滑坡
占比（%）栅格
个数栅格
占比（%）滑坡
点个数滑坡
占比（%）栅格
个数栅格
占比（%）滑坡
点个数滑坡
占比（%）极高 14766439 11.99 44 23.40 14840333 12.05 52 27.66 14951174 12.14 56 29.79 高 15554640 12.63 68 36.17 15394537 12.50 72 38.30 15283696 12.41 75 39.89 中 24114003 19.58 38 20.21 24163265 19.62 40 21.28 23929268 19.43 42 22.34 低 32968940 26.77 22 11.70 32981256 26.78 10 5.32 32599471 26.47 8 4.26 极低 35752274 29.03 16 8.51 35776905 29.05 14 7.45 36392714 29.55 7 3.72 根据评价因子的选取及易发性评价分区图可知，滑坡高和极高易发区多位于墨脱县的达木乡、帮辛乡，林芝县的丹娘、里龙、扎西饶登乡，朗县的陇村，工布江达的江达乡。在这些地区应采取相应的地质灾害防治措施。特别是位于雅鲁藏布江与尼洋河两岸海拔较低、坡度为30°～40°，距河流、道路、断裂带200 m以内的区域。
究其原因，这类地区位于雅鲁藏布江与尼洋河两岸南部与印度板块和亚欧板块交界，地壳运动剧烈，孕育一系列区域性断裂，断裂带与断层降低了岩层的完整性和强度，并且高程多位于200～1 000 m，大多数坡度小于40°，在此范围内人工多进行切坡建房和道路建设等强烈活动，造成大量的裸露斜坡，加上长期的流水作用，使河流两岸遭受严重的侵蚀和冲刷，导致沉积物饱和，从而降低斜坡的完整性，使斜坡运动或质量运动，且距道路越近，道路建设所造成的破坏性会对边坡稳定性产生负面影响，因此滑坡灾害频发。
相反，滑坡低易发区主要分布在工布江达县的错高、朱拉区，林芝市的冲果俄、港阿如，米林县的苏鲁胖地区，其特点是坡度较缓、人类活动较少，远离道路、河流、断裂带。
4.2 模型精度比较
在机器学习中，性能指标通常用于二进制分类中测试集的正确预测数。笔者使用准确度（Accuracy）、精确度（Precision）、召回率（Recall）、F1分数、（ROC）曲线和AUC值6个指标对不同机器学习模型的精度进行了评价。准确度分数是评估模型在二元分类问题中的性能的最常用指标，表示在所有样本中，能被正确识别的概率；精确度是通过计算模型预测为真时实例为正样本的频率来评估模型性能的度量；召回率是模型正确检测真阳性实例的度量；F1分数是召回率和精度之间的权衡指数，同时考虑了FP和FN，使模型整体更具准确性。具体公式如下：
$$ 准确度=\frac{TP+TN}{TP+FP+TN+FN} $$ (5) $$ 精确度=\frac{TP}{TP+FP} $$ (6) $$ 召回率=\frac{TP}{TP+FN} $$ (7) $$ {\text{F}}1 = \frac{{2Precision \times Recall}}{{Precision + Recall}} $$ (8) 式中：TP和TN分别为真阳性和真阴性，代表正确分类的像素数；FP和FN分别是假阳性和假阴性，代表错误分类的像素数。
为了得到不同机器学习算法在测试数据集上的预测准确性，基于上述方法，利用公式（5）～公式（8）计算精确度、精确度、召回率和F1指数，随机抽取30%样本作为测试样本，得出模型的泛化能力和准确率（表3）。可以看出，基于不同框架算法的预测性能不一样。3种机器学习模型中，LightGBM模型在超参数优化下其AUC（0.843 2）、ACC（0.853 1）、F1分数（0.834 5）、Precesion（0.825 1）均高于另外两种机器学习模型。
表 3 各机器学习模型准确率Table 3. Accuracy of each machine learning model机器学习模型 Gini–RF XGBoost LightGBM AUC 0.752 4 0.803 5 0.825 6 5–fold 0.822 5 0.835 8 0.843 2 ACC 0.723 4 0.814 8 0.825 6 5–fold 0.753 4 0.835 9 0.853 1 F1-score 0.775 2 0.786 7 0.802 1 5-fold 0.802 6 0.825 6 0.834 5 Precesion 0.783 4 0.796 8 0.804 5 5–fold 0.802 6 0.813 2 0.825 1 在机器学习中，ROC曲线被广泛应用于二分类问题中来评估分类器的可信度（张玘恺等，2020）。AUC为ROC曲线下面积。AUC=1表示该曲线存在至少一个阈值能得出完美预测。曲线纵轴为真阳率TPR，横轴为假阳率FPR，越靠近左上角，则认为该判断指标预测能力越好。从这条 ROC曲线可以看出，经过网格搜索与5折交叉验证后的蓝色曲线LightGBM模型更接近左上角，AUC值为0.843 2，与Gini–RF模型的0.822 5有较大提升，且准确率高于XGBoost模型的0.935 8（图10）。XGBoost相比Gini–RF而言，对模型的损失函数进行了改进，并加入了模型复杂度的正则项，而LightGBM是在XGBoost基础上，优化了模型的训练速度。因此，LightGBM的泛化能力最好，易发性划区可靠性高。
4.3 典型滑坡验证
对比近几年来雅鲁藏布江与尼洋河两岸发生的滑坡事件（表4），将9个滑坡信息导入生成的滑坡易发性图中，可知3个滑坡点位于中易发区，3个滑坡点位于高易发区，剩余均出现在极高易发区。
表 4 近几年以来滑坡事件Table 4. Landslide events in recent years地区 位置 发生时间 来源 易发性分区 林芝市加拉村 E 94°54′04″，N 29°41′45″ 2018.10.29 新华社 中 林芝市加拉村下游7公里处 E 94°54′24″，N 29°41′27″ 2022.01.22 中国青年网 中 林芝市波密县古乡索通村羌纳自然村 E 95°27′41″，N 30°00′21″ 2017.8.24 中国军视网 中 林芝市朗县辖区560国道K80处 E 92°49′24″，N 29°04′03″ 2022.7.22 朗县公安局 高 林芝市米林县派镇加拉村 E 94°54′04″，N 29°41′45″ 2018.10.17 西藏之声 高 林芝市朗县 E 93°00′48″，N 29°04′42″ 2022.7.23 朗县住建局 高 林芝市墨脱县达木乡 E 95°27′46″，N 29°29′35″ 2021.7.4 中国自然资源报 极高 国道559线波密至墨脱路段 E 97°02′03″，N 29°19′14″ 2019.5.16 西藏自治区交通运输厅 极高 林芝市墨脱县达木珞巴民族乡小学 E 95°27′52″，N 29°29′46″ 2020.8.26 新京报 极高 为了进一步验证本研究分析方法的可靠性，选择羌纳巴嘎滑坡与墨脱县公路滑坡两处滑坡现场调查进行对比验证（图11）。
西藏自治区林芝地区米林县羌纳乡巴嘎村滑坡位于E 94°24′34″，N 29°20′16″；所处地形地貌为高山河谷地貌；下付基为板岩；斜坡结构为岩土复合斜坡，坡度为30°；植被覆盖率一般，土地利用较低；滑坡前缘至斜坡下方公路，后缘至斜坡山脊处，滑坡体主要为碎土石，滑床为板岩。该滑坡变形特征主要为前方公路开挖斜坡坡脚，导致斜坡失稳。
林芝地区墨脱县公路地处E 93°38′10″，N 29°08′28″，滑坡长为30 m，宽为40 m，厚度为2 m，面积为1 200 m²，体积为2 400 m3，坡度为35°，坡向为260°，滑坡侧边界、前缘清晰可辨。该滑坡微地貌为陡坡，地层岩性为泥岩，位于白龙断层附近，斜坡结构类型为土质斜坡，坡形为凸形，滑坡下方人类活动较少，仅有一小段公路，植被覆盖率较低，为低矮灌丛，滑坡位于河流右凸岸。目前状况为不稳定。
两处滑坡均处于滑坡高易发区，再次验证了本研究机器学习模型划区的准确性。研究结果可供区域滑坡防治相关部门参考。
5. 结论
（1）统计各评价因子分级范围内滑坡点数量，表明在高程为32～1 544 m与2 722～3 752 m、坡度为30°～40°、地层岩性为雅鲁藏布江带闪片岩、土地利用为草地、距断裂带、河流与道路200 m以内滑坡发生的次数最多。
（2）采用五折交叉验证后，基于贝叶斯优化算法的Gini–RF模型准确率由原来的0.752 4提升到0.822 5，XGBoost与LightGBM模型准确率也提升了0.032 3与0.017 6。3种模型对研究区的滑坡分区都具有很高的准确性，其中LightGBM模型的性能最好，AUC值、精确度、F1分数、泛化能力、拟合程度、精确率更高。
（3）利用Gini–RF、XGBoost、LightGBM等3种集成机器学习模型对滑坡易发性进行分析，表明滑坡高和极高易发区多位于 墨脱县的达木乡、帮辛乡，林芝县的丹娘、里龙、扎西饶登乡，朗县的陇村，工布江达的江达乡。特别是位于雅鲁藏布江与尼洋河两岸海拔较低、坡度为30°～40°、距河流、道路、断裂带200 m以内的区域。在这些地区应采取相应的地质灾害防治措施。
（4）滑坡极高与高易发性区占比分别为12.14%和12.41%，低和极低易发区分别占26.47%与29.55%，区内一半以上的地区不容易发生滑坡灾害。滑坡易发性分区结果与现场滑坡灾害调查结果吻合较好，同时利用研究区近几年已发生的滑坡点进行验证，表明模型的可靠性高，滑坡分区图可为有关地方部门的防灾减灾活动提供指导。
-
![]()
图 1 甘肃积石山6.2级地震烈度图（据中华人民共和国应急管理部，2023修）
Figure 1. Intensity map of the 6.2 magnitude earthquake in Jishishan, Gansu Province
![]()
图 4 积石山县地震地质灾害规模类型统计图
Figure 4. Statistical chart of the scale and types of seismic geological hazards in Jishishan county
![]()
图 5 积石山县典型地震崩滑地质灾害
Figure 5. Typical geological hazards of earthquake collapse and sliding in Jishishan county
![]()
图 6 积石山县域地质灾害易发性评价指标图
Figure 6. Evaluation index map of geological hazard susceptibility in Jishishan County
![]()
图 7 积石山县域地质灾害易发性评价图
Figure 7. Assessment map of geological hazard susceptibility in Jishishan County
![]()
图 8 积石山县域地质灾害危险性评价分区图
Figure 8. Zoning map for geological hazard risk assessment in Jishishan County
![]()
图 9 积石山县域地质灾害危险性评价分区面积占比统计图
Figure 9. Statistical chart of the proportion of geological hazard risk assessment zoning area in Jishishan County
表 1 易发评价指标权重统计表
Table 1 Weight of susceptibility evaluation indicators
权重 坡高 坡向 坡度 工程地质岩组 灾点分布 地震烈度 矩阵法 0.2 0.08 0.19 0.22 0.2 0.11 专家打分法 0.23 0.09 0.18 0.22 0.2 0.08 综合 0.22 0.08 0.19 0.22 0.2 0.09 
下载: 导出CSV
表 2 各因子图层分类情况及其对应的信息量值统计表
Table 2 Classification of each factor layer and its corresponding information value
评价指标图层 评价因素子类 坡度 分类范围（°） ＜5 5～15 15～25 25～35 35～45 45～55 ＞55 分类值 1 2 3 4 5 6 7 信息量 0.27 0.13 −0.533 −0.87 −0.25 −0.24 −0.68 坡高 分类范围（°） 0～20 20～40 40～60 ＞60 分类值 1 2 3 4 信息量 0.22 −0.07 −1.14 −1.83 坡向 分类范围（°） 0～90 90～180 180～270 270～360 分类值 1 2 3 4 信息量 −1.57 −0.18 0.23 1.10 工程地质岩组 分类范围（°） 土体类 碎屑岩类 岩浆岩类 变质岩类 碳酸盐岩类 混合岩类 分类值 1 2 3 4 5 6 信息量 0.67 −0.27 −1.37 −2.30 −1.390 0.340 灾点分布 分类范围（m） 0～200 200～500 500～1000 >1000 分类值 1 2 3 4 信息量 −10.132 −10.693 −11.180 −11.00 地震烈度 分类范围（°） 8 7 6 分类值 1 2 3 信息量 1.690 1.350 1.16
下载: 导出CSV
表 3 积石山地质灾害易发性分区说明表
Table 3 Explanation table for the susceptibility zoning of geological hazards in Jishi Mountain
易发区划 面积
（km2）百分比
（%）分区说明 高易发区 49.54 5.45 主要分布在积石山东部黄土丘陵区的山梁地带，影响最为严重的乡镇分别为郭干乡、中咀岭乡、寨子沟乡、小关乡、徐扈家乡、刘集乡、胡林家乡 中易发区 89.34 9.83 主要分布在积石山东部黄土丘陵区的山梁地带，影响积石山县域东部大部分乡镇，主要影响大河家镇、石塬镇、柳沟乡、关家川、安吉镇、银川镇、铺川乡、癿藏乡、居集镇、郭干乡、中咀岭乡、寨子沟乡、小关乡、徐扈家乡、刘集乡、胡林家乡 低易发区 297.24 32.70 主要分布在积石山东部侵蚀剥蚀低山丘陵、侵蚀堆积河谷平原区，地质灾害发育较少 非易发区 472.88 52.02 主要分布在积石山以西高、中侵蚀山地区，地质灾害不发育
下载: 导出CSV
表 4 积石山地质灾害危险性分区面积及占比统计表
Table 4 Statistical table of the area and proportion of geological hazard zoning in Jishi Mountain
分区 10年一遇降雨 20年一遇降雨 50年一遇降雨 百年一遇降雨 面积
（km2）占比
（%）面积
（km2）占比
（%）面积
（km2）占比
（%）面积
（km2）占比
（%）极高危险区 42.57 4.68 56.73 6.24 109.45 12.04 208.56 22.94 高危险区 84.08 9.25 98.24 10.80 159.64 17.56 183.56 20.19 中危险区 283.44 31.18 306.89 33.76 349.53 38.45 389.31 42.83 低危险区 498.89 54.89 447.14 49.20 290.38 31.95 127.57 14.04
下载: 导出CSV
-
丛凯, 马宗源, 李瑞冬. 立节北山滑坡地震稳定性及破坏过程三维有限元分析[J]. 西北地质, 2023, 56(2): 283−291. CONG Kai, MA Zongyuan, LI Ruidong. Three–dimensional Analysis of Seismic Stability and Failure Process of Lijie Landslide in Gansu Province, China[J]. Northwestern Geology, 2023, 56(2): 283−291.
戴岚欣, 许强, 范宣梅, 等 . 2017年8月8日四川九寨沟地震诱发地质灾害空间分布规律及易发性评价初步研究[J]. 工程地质学报,2017 ,25 （4 ）:1151 −1164 .DAI Lanxin, XU Qiang, FAN Xuanmei, et al . A preliminary study on spatial distribution patterns of landslides triggered by Jiuzhaigou earthquake in Sichuan on August 8th, 2017 and their susceptibility assessment[J]. Journal of Engineering Geology,2017 ,25 （4 ）:1151 −1164 .冯卫, 唐亚明, 马红娜, 等. 基于层次分析法的咸阳市多灾种自然灾害综合风险评价[J]. 西北地质, 2021, 54(2): 282−288. FENG Wei, TANG Yaming, MA Hongna, et al. Comprehensive Risk Assessment of Multi-hazard Natural Disasters in Xianyang City Based on AHP[J]. Northwestern Geology, 2021, 54(2): 282−288.
黄润秋, 李为乐 . “5.12”汶川大地震触发地质灾害的发育分布规律研究[J]. 岩石力学与工程学报,2008 ,27 （12 ）:2585 −2592 .HUANG Runqiu, LI Weile . Research on development and distribution rules of geohazards induced by Wenchuan earthquake on 12th May, 2008[J]. Chinese Journal of Rock Mechanics and Engineering,2008 ,27 （12 ）:2585 −2592 .黄润秋 . 汶川地震地质灾害后效应分析[J]. 工程地质学报,2011 ,19 （2 ）:145 −151 . doi: 10.3969/j.issn.1004-9665.2011.02.001HUANG Runqiu . After effect of geohazards induced by the Wenchuan earthquake[J]. Journal of Engineering Geology,2011 ,19 （2 ）:145 −151 . doi: 10.3969/j.issn.1004-9665.2011.02.001贾丽娜, 陈世昌. 基于AHP和GIS的舟曲地质灾害易发性评价[J]. 西北地质, 2024, 57(1): 23−33. JIA Lina, CHEN Shichang. Geological Hazard Susceptibility Evaluation Based on AHP and GIS in Zhouqu County, Gansu[J]. Northwestern Geology, 2024, 57(1): 23−33.
李秀珍, 孔纪名 . 芦山和汶川地震诱发次生地质灾害的规律及特征对比分析[J]. 自然灾害学报,2014 ,23 （5 ）:11 −18 .LI Xiuzhen, KONG Jiming . Comparative analysis of development law and characteristics of secondary geo-hazard induced by Lushan and Wenchuan earthquakes[J]. Journal of Natural Disasters,2014 ,23 （5 ）:11 −18 .李为乐, 许强, 李雨森, 等. 2023年积石山Ms6.2级地震同震地质灾害初步分析[J]. 成都理工大学学报(自然科学版), 2024, 51（1）: 33−45+90. LI Weile, XU Qiang, LI Yusen, et al. Preliminary Analysis of the Coseismic Geohazards Induced by the 2023 Jishishan Ms 6.2 Earthquake[J]. Journal of Chengdu University of Technology (Science & Technology Edition), 2024, 51(1): 33−45+90.
马红娜, 刘江, 冯卫, 等. 地质灾害风险评估在国土空间规划中的应用以陕北榆林高西沟为例[J]. 西北地质, 2023, 56(3): 223−231. MA Hongna, LIU Jiang, FENG Wei, et al. Application of Geological Hazard Risk Assessment in Territorial Space Planning: A Case Study of Gaoxigou Village in Yulin City of Northern Shaanxi Province[J]. Northwestern Geology, 2023, 56(3): 223−231.
裴向军, 黄润秋 . “4·20”芦山地震地质灾害特征分析[J]. 成都理工大学学报(自然科学版),2013 ,40 （3 ）:257 −263 .PEI Xiangjun, HUANG Runqiu . Characteristics analysis of the "April 20"Lushan earthquake geological disaster[J]. Journal of Chengdu University of Technology (Science & Technology Edition),2013 ,40 （3 ）:257 −263 .铁永波, 张宪政, 卢佳燕, 等 . 四川省泸定县Ms6.8级地震地质灾害发育规律与减灾对策[J]. 水文地质工程地质,2022 ,49 （6 ）:1 −12 .TIE Yongbo, ZHANG Xianzheng, LU Jiayan, et al . Characteristics of geological hazards and it's mitigations of the Ms 6.8 earthquake in Luding County, Sichuan Province[J]. Hydrogeology Engineering Geology,2022 ,49 （6 ）:1 −12 .万飞鹏, 杨为民, 邱占林, 等. 甘肃岷县纳古呢沟滑坡-泥石流灾害链成灾机制及其演化[J]. 中国地质, 2023, 50(3): 911-925. WAN Feipeng, YANG Weimin, QIU Zhanlin, et al. Disaster mechanism and evolution of Nagune Gully landslide-debris flow disaster chain in Minxian County, Gansu Province[J]. Geology in China, 2023, 50(3): 911-925.
王东坡, 何思明, 葛胜锦, 等 . “9·07”彝良地震诱发次生山地灾害调查及减灾建议[J]. 山地学报,2013 ,31 （1 ）:101 −107 . doi: 10.3969/j.issn.1008-2786.2013.01.014WANG Dongpo, HE Siming, GE Shengjin, et al . Mountain Hazards Induced by the Earthquake of Sep 07, 2012 in Yiliang and the Suggestions of Disaster Reduction[J]. Mountain Research,2013 ,31 （1 ）:101 −107 . doi: 10.3969/j.issn.1008-2786.2013.01.014王兰民, 柴少峰, 薄景山, 等 . 黄土地震滑坡的触发类型、特征与成灾机制[J]. 岩土工程学报,2023 ,45 （8 ）:1543 −1554 . doi: 10.11779/CJGE20220531WANG Lanmin, CHAI Shaofeng, BO Jingshan, et al . Triggering types, characteristics, and disaster mechanism of seismic loess landslides[J]. Chinese Journal of Geotechnical Engineering,2023 ,45 （8 ）:1543 −1554 . doi: 10.11779/CJGE20220531王运生, 赵波, 吉锋, 等. 2023年甘肃积石山Ms6.2级地震震害异常的启示[J]. 成都理工大学学报(自然科学版), 2024, 51（1）: 1−8. WANG Yunsheng, ZHAO Bo, JI Feng, et al. Preliminary insights into hazards triggered by the 2023 Jishishan Ms 6.2 earthquake, Gansu Province[J]. Journal of Chengdu University of Technology (Science & Technology Edition), 2024, 51(1): 1−8.
吴常润, 赵冬梅, 刘澄静, 等. 基于GIS和信息量模型的陇川县滑坡易发性评价[J]. 西北地质, 2020, 53(2): 308−320. WU Changrun, ZHAO Dongmei, LIU Chengjing, et al. Landslide Susceptibility Assessment of Longchuan County Based on GIS and Information Value Model[J]. Northwestern Geology, 2020, 53(2): 308−320.
许强, 李为乐 . 汶川地震诱发大型滑坡分布规律研究[J]. 工程地质学报,2010 ,18 （6 ）:818 −826 . doi: 10.3969/j.issn.1004-9665.2010.06.002XU Qiang, LI Weile . Distribution of large-scale landslides induced by the Wenchuan earthquake[J]. Journal of Engineering Geology,2010 ,18 （6 ）:818 −826 . doi: 10.3969/j.issn.1004-9665.2010.06.002许强, 彭大雷, 范宣梅, 等. 甘肃积石山 6.2 级地震触发青海中川乡液化型滑坡-泥流特征与成因机理[J/OL]. 武汉大学学报(信息科学版), 2024: 1−19. XU Qiang, PENG Dalei, FAN Xuanmei, et al.Preliminary Study on the Characteristics and Initiation Mechanism of Zhongchuan Town Flowslide Triggered by Jishishan Ms 6.2 Earthquake in Gansu Province[J/OL]. Geomatics and Information Science of Wuhan University, 2024: 1−19.
徐敏, 裴向军, 张晓超 . 强震触发石碑塬滑坡黄土动力特性[J]. 成都理工大学学报(自然科学版),2014 ,41 （4 ）:492 −498 .XU Min, PEI Xiangjun, ZHANG Xiaochao . Dynamic characteristics of loess from Ningxia Shibeiyuan landslide triggered by Hai yuan earthquak[J]. Journal of Chengdu University of Technology (Science & Technology Edition),2014 ,41 （4 ）:492 −498 .殷跃平, 张永双, 马寅生, 等 . 青海玉树Ms7.1级地震地质灾害主要特征[J]. 工程地质学报,2010 ,18 （3 ）:289 −296 .YIN Yueping, ZHANG Yongshuang, MA Yansheng, et al . Research on major characteristic of geohazards induced by the Yushu Ms7.1 earthquake[J]. Journal of Engineering Geology,2010 ,18 （3 ）:289 −296 .殷志强, 徐永强, 陈红旗, 等 . 2013年甘肃岷县—漳县Ms6.6级地震地质灾害展布特征及主控因素研究[J]. 第四纪研究,2015 ,35 （1 ）:88 −99 .YIN Zhiqiang, XU Yongqiang, CHEN Hongsheng, et al . Study on the distribution characteristics of geohazards and the causative tectonic of the Minxian-Zhangxian Ms6.6 Earthquake on 22 July, 2013, Gansu, China[J]. Quaternary Sciences,2015 ,35 （1 ）:88 −99 .杨博, 田文通, 孙军杰, 等 . 海原大地震诱发石碑塬黄土滑坡机制探讨[J]. 地震工程学报,2020 ,42 （5 ）:1165 −1172 . doi: 10.3969/j.issn.1000-0844.2020.05.1165YANG Bo, TIAN Wentong, SUN Junjie, et al . Mechanism of SHIbeiyuan Loess Landslide Induced by the Haiyuan Earthquake[J]. China Earthquake Engineerring Journal,2020 ,42 （5 ）:1165 −1172 . doi: 10.3969/j.issn.1000-0844.2020.05.1165张宪政, 铁永波, 李光辉, 等 . 四川泸定Ms6.8级地震区湾东河流域泥石流活动性预测[J]. 地质力学学报,2022 ,28 （6 ）:1035 −1045 . doi: 10.12090/j.issn.1006-6616.20222827ZHANG Xianzheng, TIE Yongbo, LI Guanghui, et al . Characteristics and risk assessment of debris flows in the Wandong catchment after the Ms 6.8 Luding earthquake[J]. Journal of Geomechanics,2022 ,28 （6 ）:1035 −1045 . doi: 10.12090/j.issn.1006-6616.20222827中华人民共和国应急管理部. 应急管理部发布甘肃积石山6.2级地震烈度图[OL]. 中华人民共和国应急管理部. 2023-12-22. https://www.mem.gov.cn/xw/yjglbgzdt/202312/t20231222_472849.shtml 周洪福, 方甜, 韦玉婷 . 国内外地震滑坡研究: 现状、问题与展望[J]. 沉积与特提斯地质,2023 ,43 （3 ）:615 −628 .ZHOU Hongfu, FANG Tian, WEI Yuting . Research situation sand suggestions on earthquake-induced landslides[J]. Sedimentary Geology and Tethyan Geology,2023 ,43 （3 ）:615 −628 .Dai F C, Xu C, Yao X, et al . Spatial distribution of land-slide striggered by the 2008 MS 8.0 Wenchuan earthquake, China[J]. Journal of Asian EarthSciences,2011 ,40 :883 −895 .Tang C, Zhu J, Qi X, et al . Landslides induced by the Wen-chuan earthquake and the subsequent strong rainfall event: A case study in the Beichuan area of China[J]. Engineering Geology,2011 ,122 :22 −33 . doi: 10.1016/j.enggeo.2011.03.013 -
期刊类型引用(5)
1. 孟勇，张欣，陈奋宁，李佐臣，冯益民，祁尧刚，郑若蕾. 古亚洲洋构造域中西段新元古代晚期至晚泥盆世的构造运动. 地质学报. 2025(01): 58-77 . 
百度学术
2. 张永玲，张治国，刘希军，田昊，李得超，肖让，肖尧. 内蒙朝克山辉长岩中单斜辉石矿物化学特征及地质意义. 西北地质. 2024(01): 122-138 . 本站查看
3. 李平，朱涛，吕鹏瑞，王洪亮，陈隽璐. 西天山早寒武世夏特辉长岩：南天山洋早期俯冲的岩浆记录. 西北地质. 2024(03): 44-58 . 本站查看
4. 杨高学，刘晓宇，朱钊，李海，佟丽莉. 中亚造山带西准噶尔地区达尔布特蛇绿岩研究进展与展望. 西北地质. 2024(03): 1-10+301 . 本站查看
5. 刘嘉情，钟世华，李三忠，丰成友，戴黎明，索艳慧，郭广慧，牛警徽，薛梓萌，黄宇. 基于机器学习和全岩成分识别东昆仑祁漫塔格斑岩–矽卡岩矿床成矿岩体和贫矿岩体. 西北地质. 2023(06): 41-56 . 本站查看
其他类型引用(0)
计量
- 文章访问数: 209
- HTML全文浏览量: 34
- PDF下载量: 70
- 被引次数: 5
