Simulation prediction and risk evaluation of debris flow in gullyprone ditches of Lajing Village, Lanping County, Yunnan Province, China
-
摘要:
兰坪县位于云南省西北部，地质环境不良，近年来该地区受降雨影响泥石流灾害频发，但是该地区缺乏较为精确的灾害风险评价。为对该地区泥石流灾害进行防范预警，文章应用FLO-2D软件对兰坪县啦井村1#典型泥石流河灾害进行动力过程反演分析，基于反演所得模拟参数对2#沟在20年、50年和100年三种降雨重现期工况进行预测分析和灾害风险评价。结果表明：啦井村2#沟道泥石流淤积深度和最大流速随降雨强度的增大呈增加趋势，极端降雨条件下沟道底部淤积深度将超过2.5 m，流通区最大流速将大于5 m/s，泥石流携带更多的松散物于堆积区堆积，易造成安全隐患；啦井村2#沟中高风险区主要分布在沟道中流速泥深较大、沟道周围民房较为集中位置以及泥石流堆积区域。研究结果可为该地区防灾减灾工程与泥石流早期预警等研究提供一定参考。
Abstract:Lanping County, located in the northwestern part of Yunnan Province, is characterized by unfavorable geological conditions. In recent years, debris flow disasters have frequently occurred in this area due to the influence of rainfall, yet there is a lack of accurate disaster risk evaluations. To enhance early warning preparedness for debris flow disasters in the region, the FLO-2D model was employed to perform an inverse analysis of a historical debris flow event in Lajing Village, Lanping County. Based on the parameters obtained from this analysis, predictions of the dynamic characteristics and risk assessments were conducted for another debris flow-prone gully under three rainfall recurrence period conditions: 20, 50 and 100 years. The results indicate that in Gully No. 2 of Lajing Village, the siltation depth and maximum flow velocity of debris flow increase with the intensification of rainfall. Under extreme rainfall conditions, the siltation depth at the trench bottom is expected to exceed 2.5 m, and the maximum flow velocity in the transportation zone may exceed 5 m/s, leading to the accumulation of more loose materials in the deposition area and creating significant safety hazards. The high-risk zones in gully No. 2 are primarily located in areas with high flow velocity and mud depth, around densely populated residential areas near the gully, and within the debris flow accumulation zone. The findings of this study provide valuable references for disaster prevention and mitigation engineering, as well as for early warning systems in the region.
-
Keywords:
- debris flow /
- FLO-2D /
- numerical simulation /
- dynamic analysis /
- risk evaluation
-
0. 引言
地基稳定性受多种因素影响[1 − 2]，特别是黄土地区，受震陷影响很大。黄土震陷作为一种严重的地质灾害，一旦发生将会造成非常严重的后果。前人对与黄土震陷的研究基本是以动三轴试验为主，以单个震陷量作为震陷性的主要的评价指标[3 − 6]。然而黄土地基震陷的准确性评价受到土性参数不确定性的影响，其中物理力学指标的概率分布类型和变异性会对评价结果产生影响。针对岩土工程中参数、计算模型、荷载等不确定性问题，其主要的不确定性分析方法包括：概率分析方法[7]、随机有限元法[8]、灰色系统评价法[9]、模糊数学法[10]、神经网络评价方法[11]等。这些不确定性方法可以考虑得更加全面，从而提高评价结果的准确性，使结果更为科学合理。
关于黄土震陷在不确定性方面的研究，李波等[12]采用灰色关联方法和逐步判别方法建立了耦合式软土震陷判别模型从而实现软土震陷等级的判定。在黄土震陷概率评价方面，谭芳[13]将概率地震危险性分析方法引入震陷计算中，采用不确定性分析方法分析了黄土震陷与地震参数之间的关系，从而提出一种考虑地震随机性的黄土震陷预测方法。任森等[14]基于大量动三轴试验，结合场地地震危险性概率分析方法，确定了研究区内的震陷区划。孙军杰等[15]依据残余应变估算模型采用概率统计等非确定性分析方法对地基动力沉降进行了概率分析。可见概率分析方法在土体动力沉降中已经得到了初步研究。
然而到目前为止，还没有研究将可靠度理论应用于黄土地基震陷的分析中。因此本文尝试将可靠度理论应用于黄土地基震陷评价中，使用可靠度方法对黄土地基震陷性进行分析。可靠度方法和传统的确定性分析方法相互补充与结合，有助于进行震陷可靠度分析。
1. 黄土地基震陷可靠度指标
可靠度理论在边坡稳定分析中被广泛应用，而边坡工程的不确定因素主要为土体参数的抗剪强度（内摩擦角、黏聚力）[16 − 18]。其依托蒙特卡洛分析法通过计算出多组安全系数值将小于1的安全系数的次数与计算次数的比值称作失效概率。因此可求得对应的可靠度指标（$ \mathrm{\beta } $）为：
$$ \beta ={\varPhi }^{-1}(1-{P}_{\text{f}}) $$ (1) 式中：$ \varPhi $——标准正态分布的累计分布函数；
${P_{\text{f}}}$——失效概率。
在目前已有的地震灾害评估中，尚未有关于地基震陷量的标准划分依据。根据中国地震局发布的《震害调查及地震损失评定工作指南》，可以采用“地基基础破坏等级标准”来评估震陷灾害的程度[16]，将震陷量小于2 cm的定义为基本完好，震陷量大于40 cm的定义为破坏。
参考可靠度理论在边坡稳定性评价中的应用，结合黄土动力学对震陷等级的划分，本文提出了一种对黄土地基震陷评价的可靠度指标。
（1）最大可靠度
在震陷等级划分中将震陷量大于40 cm的定义为震陷破坏，因此可采用破坏概率（Pd）来判定其震陷破坏状态，即$ {P_{\mathrm{d}}} = P\left( {S \geqslant 40} \right) $。将$ \beta_{\mathrm{d}} $定义为：
$$ {\beta _{\mathrm{d}}} = \frac{{40 - {u_{\mathrm{s}}}}}{{{\sigma _{\mathrm{s}}}}} $$ (2) 式中：$ {u_{\mathrm{s}}} $——震陷量的平均值；
$ {\sigma _{\mathrm{s}}} $——震陷量的标准差。
如果震陷量符合标准正态分布，破坏概率与可靠度的关系可表示为：
$$ {P_{\mathrm{d}}} = P\left( {S > 40} \right) = 1 - \varPhi \left(\frac{{40 - {u_{\mathrm{s}}}}}{{{\sigma _{\mathrm{s}}}}}\right) = 1 - \varPhi ({\beta _{\mathrm{d}}}) $$ (3) 则根据可靠度理论，将破坏概率对应的可靠度称为最大可靠度（${\beta }_{\text{d}} $），也就是指衡量地基震陷最大安全程度的标准，其计算公式为：
$$ {\beta }_{\text{d}}={\varPhi }^{-1}(1-{P}_{\text{d}}) $$ (4) 根据蒙特卡罗可靠度计算方法可定义破坏概率为：
$$ {P_{\text{d}}} = \frac{M}{N} $$ (5) 式中：$ {P}_{{\mathrm{d}}} $——地基震陷破坏概率；
M——抽样次数中震陷量大于40 cm的次数；
N——蒙特卡洛分析法总的抽样次数。
（2）最小可靠度
在震陷等级划分中将震陷量小于2 cm的定义为基本完好，因此可采用完好概率（Ph）来判定其完好状态，即$ {P_{\mathrm{h}}} = P\left( {S \leqslant 2} \right) $。将$ {\beta _{\mathrm{h}}} $定义为：
$$ {\beta _{\mathrm{h}}} = \frac{{2 - {u_{\mathrm{s}}}}}{{{\sigma _{\mathrm{s}}}}} $$ (6) 式中：$ {u_{\mathrm{s}}} $——震陷量的平均值；
$ {\sigma _{\mathrm{s}}} $——震陷量的标准差。
如果震陷量符合标准正态分布，完好概率（Ph）与可靠度的关系可表示为：
$$ {P_{\mathrm{h}}} = P\left( {S < 2} \right) = \varPhi \left(\frac{{2 - {u_{\mathrm{s}}}}}{{{\sigma _{\mathrm{s}}}}}\right) = \varPhi ({\beta _{\mathrm{h}}}) $$ (7) 根据可靠度理论，将完好概率对应的可靠度称为最小可靠度（${\beta _{\mathrm{h}}} $），也就是指衡量地基震陷完好程度的标准，其计算公式为：
$$ {\beta _{\text{h}}} = {\varPhi ^{ - 1}}\left( {{P_{\text{h}}}} \right) $$ (8) 根据蒙特卡罗可靠度计算方法可定义完好概率为：
$$ {P_{\text{h}}} = \frac{M}{N} $$ (9) 式中：$ {P}_{{\mathrm{h}}} $——地基完好概率；
M——抽样次数中震陷量小于2 cm的次数；
N——蒙特卡洛分析法总的抽样次数。
2. 基于Kriging代理模型的震陷可靠度计算
2.1 Kriging代理模型的建立
Kriging模型是一种基于方差最小的无偏估计方法，它的原理是通过已知样本点的线性加权法来预测未知点的函数值，即已知加权系数插值就可求得变量空间内任一点的函数值。Kriging代理模型通常情况下由回归部分和随机误差两部分所组成。将其用输出变量S(x)与输入变量x表示为：
$$ S({\boldsymbol{x}}) = \sum\limits_{i = 1}^m {{\beta _i}} f({\boldsymbol{x}}) + {\textit{z}}({\boldsymbol{x}}) = f{({\boldsymbol{x}})^T}{\boldsymbol{\beta}} + {\textit{z}}({\boldsymbol{x}}) $$ (10) 式中：每个样本点中输入变量的个数为n，$ {\boldsymbol{x}} = [ {x_1},{x_2}, \cdots , {x_n} ]^T$，$ f(x) = {\left[ {{f_1}(x),{f_2}(x), \cdots ,{f_m}(x)} \right]^T} $为回归多项式基函数向量，m为回归多项式的数量；${\boldsymbol{\beta}} = {[{\beta _1},{\beta _2},\cdots,{\beta _m}]^T}$为多项式参数向量，z(x)为一均值为0，方差为$ {\sigma ^2} $的正态随机过程，该随机过程的协方差为：
$$ COV\left[{\textit{z}}(u),{\textit{z}}(v)\right]={\sigma }^{2}R({u,v}) $$ (11) 其中，R表示任意两点u，v的相关函数，通常取为高斯相关方程，即：
$$ R(u,v) = \exp \left\{ { - \theta {{(u - v)}^2}} \right\} $$ (12) 其中，$ \theta $为相关参数，可通过最大似然估计得到：
$$ \theta=\max \left\{ { - \left[ {\frac{1}{2}\ln \left( {\left| R \right|} \right)} \right] + \frac{k}{2}\ln \left( {{\sigma ^2}} \right)} \right\} $$ (13) 其中，R为相关矩阵，元素R=R(u,v)。
基于给定的样本点，多项式参数向量$\beta $与随机过程方差$ {\sigma ^2} $的估计值分别为：
$$ \left\{ {\begin{split} & {\hat {\boldsymbol{\beta}} = {{\left( {{{\boldsymbol{F}}^T}{{\boldsymbol{R}}^{ - 1}}{\boldsymbol{F}}} \right)}^{ - 1}}{{\boldsymbol{F}}^T}{{\boldsymbol{R}}^{ - 1}}{\boldsymbol{Y}}} \\ & {{{\hat {\boldsymbol{\sigma}} }^2} = \frac{1}{k}{{\left( {{\boldsymbol{Y}} - {\boldsymbol{F}}\hat {\boldsymbol{\beta}} } \right)}^T}{{\boldsymbol{R}}^{ - 1}}\left( {{\boldsymbol{Y}} -{\boldsymbol{ F}}\hat {\boldsymbol{\beta}} } \right)} \end{split}} \right. $$ (14) 其中，F为由回归多项式函数值构成的矩阵。
对于任意输入变量$ {x_0} $,假定$ {x_0} $与输入变量样本X之间的相关系数向量为$ {{\boldsymbol{r}}_0} = {\left[ {R\left( {{x_0},{x_1}} \right), \cdots ,R\left( {{x_0},{x_k}} \right)} \right]^T} $利用式（10）构建的Kriging模型表达式，可以得到与$ {x_0} $对应的系统响应值$ S({x_0}) $的预测值$ \hat S({x_0}) $服从正态分布$ N(\hat S({x_0}),{\hat \sigma _{\hat S}}^2({x_0})) $,$ S({x_0}) $的最佳线性无偏估计及均方误差可表示为：
$$ \left\{ {\begin{split} &{\hat S\left( {{x_0}} \right) = {f^T}({x_0})\hat {\boldsymbol{\beta}} + {{\boldsymbol{r}}_0}^T{{\boldsymbol{R}}^{ - 1}}({\boldsymbol{Y}} - {\boldsymbol{F}}\hat {\boldsymbol{\beta }})} \\ & {{{\hat \sigma }^2}_{\hat S}({x_0}) = {{\hat \sigma }^2}\left( {1 + {{\boldsymbol{u}}^T}{{({{\boldsymbol{F}}^T}{{\boldsymbol{R}}^{ - 1}}{\boldsymbol{F}})}^{ - 1}}{\boldsymbol{u}} - {{\boldsymbol{r}}_0}^T{{\boldsymbol{R}}^{ - 1}}{{\boldsymbol{r}}_0}} \right)} \end{split}} \right. $$ (15) 其中，$ {{{\boldsymbol{u}} = }}{{{{\boldsymbol{F}}}}^T}{{\boldsymbol{R}}^{ - 1}}{{\boldsymbol{r}}_0} - {\boldsymbol{f}} $。
本文进行震陷可靠度分析时，式（10）—（15）中应用Kriging模型构建和预测其过程，采用MATLAB中的DACE工具箱实现。
2.2 黄土地基数值模型
采用FLAC3D有限差分软件建立黄土地基数值模型，模型尺寸设置为15 m×15 m×15 m，模型网格数为
3375 ，本构模型选择Mohr-Coulomb弹塑性本构模型。在静力分析时，底面和侧面选用固定边界，限制底面与侧面的位移和速度，顶部为自由边界。动力分析与静力分析不同，由于地震荷载会存在波的反射，因此在模型底面设置静态边界，侧面设置为自由场边界，模拟土体在实际中的状态（图1）。本文选用三条地震波分别是Kobe地震波、Imperial地震波和一条人工合成波。为方便后续说明，这里将选用的三条地震波分别命名为地震波1、地震波2及地震波3。地震波1的持续时间为25 s，时间间隔为0.02 s，最大加速度为0.345g，地震波2的持续时间为40 s，时间间隔为0.01 s，最大加速度为0.313g，地震波3的持续时间为20 s，时间间隔为0.02 s，最大加速度为0.153g。图2为经过滤波以及基线校正后所选地震波的时程曲线。
为了研究不同地震水平作用下震陷变形规律，将其改变幅值的大小来满足不同工况的需要，得到 0.1g、0.2g、0.3g、0.4g、0.5g、0.6g、0.7g不同峰值加速度大小的7种水平。其计算所得震陷量结果如表1所示。
表 1 不同地震波及不同峰值加速度作用下地基震陷量Table 1. Seismic subsidence of foundation under different seismic waves and peak accelerations.峰值加速度/g 震陷量/cm 地震波1 地震波2 地震波3 0.1 0.72 1.95 2.65 0.2 2.37 8.40 13.44 0.3 9.10 16.02 21.62 0.4 17.83 21.44 29.52 0.5 19.78 25.29 37.06 0.6 20.08 30.67 40.98 0.7 20.48 35.84 46.98 2.3 震陷变形影响因素敏感性分析
岩土体材料参数存在着变异性，但土体泊松比的变异系数非常小，在考虑土性参数不确定时可将其看作常量。孙军杰等[17]通过理论分析，给出了非饱和黄土动残余应变与黏聚力、内摩擦角、初始孔隙比、土层重度、土层埋深以及地震动荷载之间的定量关系式。压缩模量是体现土体压缩性的指标之一，作为土体的力学性质，在分析地基土体变形或沉降时起着重要作用，而土体的黏聚力和内摩擦角反映土体的强度特征，在静力分析中，黄土地基静力参数中密度对于计算结果的影响较明显。因此，这里选择密度、弹性模量、黏聚力及内摩擦角四种参数。采用FLAC3D结合灰色关联分析对影响黄土地基震陷变形的因素做敏感性分析，求得其密度、弹性模量、黏聚力及摩擦角四个参数的灰色关联度为$ q = [ {0.656\;8},\;{0.856\;2},\;{0.668\;5},\;{0.728\;1} ] $，敏感性排序为：弹性模量>内摩擦角>黏聚力>密度。
2.4 土性参数分布概率模型
本节选取密度（ρ）、压缩模量（Es）、黏聚力（c）及内摩擦角（φ）四种参数作为黄土震陷性评价的参数指标来进行分布概率模型研究。为了全面分析不同地貌单元4种土性参数的最优概率分布，选用可以包含大部分不确定性因素的概率分布形式包括正态分布、对数正态分布和威布尔分布进行假设检验。从而确定每个参数的最优概率分布形式[18]。
为获取原状黄土地基土性不同参数的分布概率模型，基于引汉济渭工程的渭河二级阶地地貌土体为例，通过对原状黄土的各种物理力学指标进行统计分析，得出各指标服从的概率密度函数及分布参数，图3所示为渭河二级阶地分布概率模型图，其中$ f\left( x \right) $表示概率密度函数。
![]()
图 3 不同参数的概率密度分布图Figure 3. Probability density distribution plots of different parameters由图3可知，从分布形态上来看，$ \rho $值的正态分布和对数正态分布整体较一致，两者均属于左偏态；其余三种土性参数正态分布、对数正态分布及威布尔分布三者分布较离散，差异性较大。
通过K-S检验法对4种土性参数的3种分布概率模型进行检验，求出样本值的理论累计概率分布与经验累计概率分布的最大差值（Dmax），与由显著性水平和样本数确定临界值（$ D_n^\alpha $）进行比较，计算结果如表2所示。四种参数均可接受正态分布、对数正态分布和威布尔分布。$ c $、$ {E_{\mathrm{s}}} $值的最优分布是正态分布，$ \varphi $、$ \rho $值的最优分布是威布尔分布，考虑到正态分布的优良性质，将正态分布作为4种土性参数的分布概率模型用于震陷可靠度计算。
表 2 分布假设检验结果Table 2. Distribution hypothesis test results检验参数 样本量 验证结果$ {D_{\max }} $ 可接受的临界值
$ D_n^\alpha /\left( {\alpha = 0.05} \right) $正态分布 对数正态分布 威布尔分布 结果 $ {D_{\max }} $ 结果 $ {D_{\max }} $ 结果 $ {D_{\max }} $ 内摩擦角 24 √ 0.180 √ 0.223 √ 0.093 0.273 黏聚力 24 √ 0.099 √ 0.125 √ 0.102 0.273 压缩模量 24 √ 0.094 √ 0.146 √ 0.164 0.273 密度 48 √ 0.150 √ 0.185 √ 0.115 0.196 2.5 基于Kriging代理模型的震陷可靠度计算方法
本节基于Kriging代理模型对地基震陷可靠度进行计算，这里进行
10000 次模拟，计算步骤如下：（1）根据前文统计出的土性参数（密度、压缩模量、黏聚力、摩擦角）的概率分布模型，确定地基土体参数的均值、标准差、分布类型等统计特性。
（2）FLAC3D6.0内置有Python编程语言，可以实现蒙特卡洛抽样与连续计算，利用Itasca开发的专门用于Python扩展编程的程序库，通过编程实现随机变量拉丁超立方抽样（LHS）、模型的建立以及计算，即会产生N组计算方案$ \left(\rho_1,E_{\mathrm{s}1},c_1,\varphi_1\right),\left(\rho_2,E_{\mathrm{s}2},c_2,\varphi_2\right),\cdots,(\rho_n,E_{\mathrm{s}n},c_n, \varphi_n) $，然后计算出对应的响应值。将计算模型文件保存为“n.f3sav”，后续模型按顺序编号，并将其作为源文件以便后续提取。
（3）对地震波设计不同的工况组合，计算每组工况下LHS的样本所对应的震陷量，并保存计算模型。
（4）利用Python编程语言从结果文件中批量提取竖向变形。
（5）采用MATLAB软件分别建立不同工况下的Kriging代理模型，以4种土体参数为输入，相应的地基震陷量为输出，对代理模型进行训练，得到满足精度要求的代理模型。经适当训练的代理模型能够较好地逼近实际情况，在此基础上可以十分方便地预测震陷值从而实现震陷可靠度的计算（图4）。
![]()
图 4 基于Kriging代理模型的震陷可靠度计算流程Figure 4. Calculation process of seismic subsidence reliability based on the Kriging surrogate model3. 算例分析
为验证本文提出方法的适用性，本节以引汉济渭工程为背景结合Kriging代理模型的震陷概率计算方法开展黄土地基震陷可靠度分析。
（1）最大可靠度
表3为不同地震波及不同地震峰值加速度作用下地基震陷的破坏概率，以渭河二级阶地为例。
表 3 不同地震波作用下地基震陷破坏概率Table 3. Probability of foundation subsidence failure under different seismic wave effects加速度/g 地震波1 地震波2 地震波3 0.1 0 0 0 0.2 0 0 0 0.3 0 0 0 0.4 0 0.072 0.116 0.5 0 0.17 0.378 0.6 0 0.232 0.704 0.7 0.002 0.264 0.888 由表3可知，在地震波1的作用下，该地基土体在地震波峰值加速度为0.7g时会出现破坏的概率，在地震波2和地震波3的作用下，在峰值加速度达到0.4g时会出现破坏的概率。三条地震波在峰值加速度达到最大时，其产生的破坏概率也有巨大差异，地震波1在峰值加速度为0.7g时破坏概率为0.002，地震波2在峰值加速度为0.7g时破坏概率为0.264，地震波3在峰值加速度为0.7g时破坏概率为0.888。
依据式（4）求得最大可靠度指标，图5为不同地震波作用下地基震陷的可靠度指标随峰值加速度的变化规律，最大可靠度指标的范围为−1.216～4.526。总体上可靠度指标随峰值加速度的增加呈非线性减小。但三条地震波的减小斜率并不相同，在地震波3作用下，可靠度指标随峰值加速度的变化斜率最大，说明在地震波3作用下，地基的可靠度受地震波峰值加速度的影响最大。对于已知地震波峰值加速度，可以快速粗略估计地基震陷的破坏概率，而不是单纯以震陷量来评价，这在风险评估中具有更为重要的意义。
![]()
图 5 可靠度分析获得的地基震陷最大可靠度Figure 5. Maximum reliability index of foundation seismic subsidence obtained from reliability analysis（2）最小可靠度
表4为不同地震波及不同地震峰值加速度作用下地基震陷的破坏概率，以渭河二级阶地为例。由表4可知，在地震波1的作用下，该地基土体在地震波峰值加速度为0.1g～0.7g时均会出现完好概率，在地震波2和地震波3的作用下，在峰值加速度分别达到0.4g及以后和0.3g及以后完好概率均为0。
表 4 不同地震波作用下地基震陷完好概率Table 4. Probability of intact seismic subsidence of foundation under different seismic wave effects加速度/g 地震波1 地震波2 地震波3 0.1 0.946 0.854 0.842 0.2 0.392 0.094 0.832 0.3 0.238 0.004 0 0.4 0.01 0 0 0.5 0.018 0 0 0.6 0.012 0 0 0.7 0.002 0 0 依据式（8）求得最小可靠度指标，图6为在不同地震波作用下地基震陷的最小可靠度随峰值加速度的变化规律，最小可靠度指标的变化范围为−2.878～1.607。从整体来说，随着地震峰值加速度的增大，地基在地震作用下的最小可靠度指标逐渐降低。当地震峰值加速度为0.1g时，震陷最小可靠度较高，都为正值；当地震峰值加速度为0.2g时，地震最小可靠度显著减小，且超过0.2g以后最小可靠度均为负值。
![]()
图 6 可靠度分析获得的地基震陷的最小可靠度指标Figure 6. Minimum reliability index of foundation seismic subsidence obtained from reliability analysis（3）变异系数对震陷可靠度的影响
为了更好地分析随机变量的变异性大小对震陷可靠度的影响，本文计算了各变量在不同变异系数组合情况下，四种土体参数变异系数随地基震陷可靠度的变化情况，工况组合见表5。
表 5 各参数变异系数工况组合表Table 5. Combination table of variation coefficient for each parameter工况 COV 密度 弹性模量 黏聚力 摩擦角 1-1 0.01 0.3 0.2 0.3 1-2 0.05 1-3 0.1 1-4 0.15 1-5 0.2 2-1 0.05 0.05 0.2 0.3 2-2 0.1 2-3 0.2 2-4 0.3 2-5 0.4 3-1 0.05 0.3 0.05 0.3 3-2 0.1 3-3 0.2 3-4 0.3 3-5 0.4 4-1 0.05 0.3 0.2 0.1 4-2 0.2 4-3 0.3 4-4 0.4 4-5 0.5 计算结果如图7所示。由图可知，随各参数变异系数的增加，震陷最大可靠度逐渐减小，而最小可靠度逐渐增大。当$ COV(\rho) $从0.01增加至0.2，相应的$ {\beta _{\text{d}}} $从3.72减小至3.05，最小可靠度（$ {\beta _{\text{h}}} $）从−1.72增加至−1.55，降幅为21.97%，增幅为10.96%；当$ COV(E_{\mathrm{s}}) $从0.05增加至0.4，相应的$ {\beta _{\text{d}}} $从3.89降至3.29，$ {\beta _{\text{h}}} $从−1.89增加至−1.64，降幅为18.24%，增幅为15.24%；当$ COV(c) $从0.05增加至0.4，$ {\beta _{\text{d}}} $从4.26降至2.89，$ {\beta _{\text{h}}} $从−1.81增加至−1.60，降幅为47.40%，增幅为13.13%；当$ COV(\varphi)$从0.1增加至0.5，$ {\beta _{\text{d}}} $从3.72降至2.74，$ {\beta _{\text{h}}} $从−1.99增加至−1.63，降幅为35.77%，增幅为22.09%。对比数据，在同样的增幅下，最大可靠度的变化幅度比最小可靠度的变化幅度大说明最大可靠度对变异系数更为敏感。计算各个参数最大可靠度起末两点斜率的绝对值，分别为3.526、1.714、3.914和2.45，可看出黏聚力变异系数对震陷最大可靠度影响更为显著。
![]()
图 7 不同参数变异系数对震陷可靠度的影响规律Figure 7. The influence of different parameter coefficient of variations on seismic subsidence reliability4. 结论
（1）Kriging代理模型可以在小样本抽样的情况下得到较好的近似拟合效果，达到精度要求，提高了震陷可靠度的计算效率。
（2）采用震陷可靠度方法来评价地基土体的震陷性，能够合理地考虑土性参数的变异性，更加符合实际情况，具有更好的可靠性和适用性，对于已知地震波峰值加速度，可以快速粗略估计地基震陷的破坏概率，而不是单纯的以震陷量来评价，这在风险评估中具有更为重要的意义。
（3）综合4种土体参数的变异系数对震陷可靠度的影响可以发现，最大可靠度对各个参数的变异系数更为敏感，此外四种土体参数中黏聚力的变异系数对震陷最大可靠度的影响更为显著。
-
![]()
图 2 啦井村1#沟泥石流模拟结果
Figure 2. Simulation results of debris flow in gully No. 1, Lajing Village
![]()
图 3 不同降雨工况下啦井村2#沟泥石流淤积深度分布
Figure 3. Debris flow deposition depth distribution in gully No. 2, Lajin Village, under different rainfall conditions
![]()
图 4 不同降雨工况下啦井村2#沟泥石流流速分布
Figure 4. Flow velocity distribution of debris flow in gully No. 2, Lajing Village, under different rainfall conditions
表 1 沟道现场调查数据
Table 1 The site investigation data of each ditch
参数 流域面积/km2
（1#/2#）沟长/ km
（1#/2#）平均坡降/‰
（1#/2#）取值 0.96/0.99 2.11/2.08 301/432 
下载: 导出CSV
表 3 FLO-2D手册建议的层流阻力系数值
Table 3 Recommended laminar flow resistance coefficient values from the FLO-2D manual
地面条件 K值范围 级配土 90～400 被侵蚀黏土 100～500 稀疏植被 1000 ～4000 矮草原 3000 ～10000
下载: 导出CSV
表 4 FLO-2D数值模型主要物理力学参数
Table 4 Main physical and mechanical parameters of the FLO-2D numerical model
参数 曼宁
系数摩擦
系数泥石流重度
/（t·m−3）黏滞系数
（1#/2#）屈服应力kPa
/（1#/2#）取值 0.07 2280 1.65 3644.13 /1860.49 4.53/2.62
下载: 导出CSV
表 5 泥石流强度影响划分表[26]
Table 5 Classification table of debris flow intensity impacts
泥石流强度 泥深/m 关系 泥深与最大流速乘积/（m2·s−1） 高 H≥2.5 OR VH≥2.5 中 0.5≤H<2.5 AND 0.5≤VH<2.5 低 0<H<0.5 AND VH<0.5
下载: 导出CSV
表 6 易损性计算公式
Table 6 Vulnerability calculation formula
易损性指标 社会易损性 建筑易损性 交通易损性 经济易损性 表达式 $F = \displaystyle\sum_{{{i}} = 1}^n {{F_{i,j}}} $
F为社会易损性；Fi,j为主级指标
权重标准化的值乘以次级
指标权重标准化的值$M = \displaystyle\sum_{{{i}} = 1}^n {{A_{i,j}}} $
M为物质易损性值（建筑/交通）；
Ai,j为承灾体主级指标权重标准化
的值乘以次级指标权重标准化的值$J = G/N$
J为经济易损性，无量纲，G为区域内每个
承灾体家庭年收入（万元）；
N为区域人均年收入（万元）
下载: 导出CSV
-
[1] LIN Qigen,WANG Ying. Spatial and temporal analysis of a fatal landslide inventory in China from 1950 to 2016[J]. Landslides,2018,15(12):2357 − 2372. DOI: 10.1007/s10346-018-1037-6
[2] ZHANG Shuai,LI Can,PENG Jingyu,et al. Fatal landslides in China from 1940 to 2020:Occurrences and vulnerabilities[J]. Landslides,2023,20(6):1243 − 1264. DOI: 10.1007/s10346-023-02034-6
[3] HE Songtang,WANG Daojie,ZHAO Peng,et al. A review and prospects of debris flow waste-shoal land use in typical debris flow areas,China[J]. Land Use Policy,2020,99:105064. DOI: 10.1016/j.landusepol.2020.105064
[4] 何晓英,唐红梅,陈洪凯. 浆体黏度和级配颗粒组合条件下泥石流冲击特性模型试验[J]. 岩土工程学报,2014,36(5):977 − 982. [HE Xiaoying,TANG Hongmei,CHEN Hongkai. Experimental study on impacting characteristic of debris flow considering different slurry viscosities,solid phase ratios and grain diameters[J]. Chinese Journal of Geotechnical Engineering,2014,36(5):977 − 982. (in Chinese with English abstract)] HE Xiaoying, TANG Hongmei, CHEN Hongkai. Experimental study on impacting characteristic of debris flow considering different slurry viscosities, solid phase ratios and grain diameters［J］. Chinese Journal of Geotechnical Engineering, 2014, 36（5）: 977 − 982. （in Chinese with English abstract）
[5] KEVIN H, HYUNCHUL J, LAWRENCE L, et al. Geohazard study of mudslide event at a subsea escarpment and structural response of crossing flowlines: Part I – FEA Analysis[C].Proceedings of the ASME 2022 41st International Conference on Ocean, Offshore and Arctic Engineering. Volume 9: Offshore Geotechnics. Hamburg, Germany. 2022.
[6] WANG Weiping,YANG Jiansi,WANG Yanbin. Dynamic processes of 2018 Sedongpu landslide in Namcha Barwa–Gyala Peri massif revealed by broadband seismic records[J]. Landslides,2020,17(2):409 − 418. DOI: 10.1007/s10346-019-01315-3
[7] O’BRIEN J S,JULIEN P Y,FULLERTON W T. Two-dimensional water flood and mudflow simulation[J]. Journal of Hydraulic Engineering,1993,119(2):244 − 261. DOI: 10.1061/(ASCE)0733-9429(1993)119:2(244)
[8] CHRISTEN M,KOWALSKI J,BARTELT P. RAMMS:numerical simulation of dense snow avalanches in three-dimensional terrain[J]. Cold Regions Science and Technology,2010,63(1):1 − 14.
[9] 段学良,马凤山,郭捷,等. 基于Massflow模型的西藏仁布杰仲沟泥石流运动特征分析[J]. 中国地质灾害与防治学报,2019,30(6):25 − 33. [DUAN Xueliang,MA Fengshan,GUO Jie,et al. Movement characteristics of Jiezhonggou debris flow of Renbu,Tibet based on Massflow model[J]. The Chinese Journal of Geological Hazard and Control,2019,30(6):25 − 33. (in Chinese with English abstract)] DUAN Xueliang, MA Fengshan, GUO Jie, et al. Movement characteristics of Jiezhonggou debris flow of Renbu, Tibet based on Massflow model［J］. The Chinese Journal of Geological Hazard and Control, 2019, 30（6）: 25 − 33. （in Chinese with English abstract）
[10] 史继帅,姜亮,翟胜强. 四川甘洛县黑西洛沟“8•31”泥石流动力过程分析[J]. 中国地质灾害与防治学报,2024,35(3):1 − 9. [SHI Jishuai,JIANG Liang,ZHAI Shengqiang. Analysis on dynamic process of the “8•31” debris flow in Luoxi gulley of Ganluo County,Sichuan Province[J]. The Chinese Journal of Geological Hazard and Control,2024,35(3):1 − 9. (in Chinese with English abstract)] SHI Jishuai, JIANG Liang, ZHAI Shengqiang. Analysis on dynamic process of the “8•31” debris flow in Luoxi gulley of Ganluo County, Sichuan Province［J］. The Chinese Journal of Geological Hazard and Control, 2024, 35（3）: 1 − 9. （in Chinese with English abstract）
[11] SHEN Wei,WANG Dongpo,HE Siming,et al. Numerical assessment of the impeding effect of check dams in the Hongchun debris flow gully,Sichuan Province,China[J]. Bulletin of Engineering Geology and the Environment,2020,79(6):2833 − 2845. DOI: 10.1007/s10064-020-01755-5
[12] 刘军友. 基于不同方法的泥石流活动强度对比研究[D]. 北京:中国地质大学(北京),2021. [LIU Junyou. Comparative study on debris flow activity intensity based on different methods[D]. Beijing:China University of Geosciences,2021. (in Chinese with English abstract)] LIU Junyou. Comparative study on debris flow activity intensity based on different methods［D］. Beijing: China University of Geosciences, 2021. （in Chinese with English abstract）
[13] 王一鸣,殷坤龙. 台风暴雨型泥石流单沟危险度研究[J]. 水文地质工程地质,2018,45(3):124 − 130. [WANG Yiming,YIN Kunlong. A study of the typhoon-triggered debris flow hazard degree of a single gully[J]. Hydrogeology & Engineering Geology,2018,45(3):124 − 130. (in Chinese with English abstract)] WANG Yiming, YIN Kunlong. A study of the typhoon-triggered debris flow hazard degree of a single gully［J］. Hydrogeology & Engineering Geology, 2018, 45（3）: 124 − 130. （in Chinese with English abstract）
[14] 孙滨,祝传兵,康晓波,等. 基于信息量模型的云南东川泥石流易发性评价[J]. 中国地质灾害与防治学报,2022,33(5):119 − 127. [SUN Bin,ZHU Chuanbing,KANG Xiaobo,et al. Susceptibility assessment of debris flows based on information model in Dongchuan,Yunnan Province[J]. The Chinese Journal of Geological Hazard and Control,2022,33(5):119 − 127. (in Chinese with English abstract)] SUN Bin, ZHU Chuanbing, KANG Xiaobo, et al. Susceptibility assessment of debris flows based on information model in Dongchuan, Yunnan Province［J］. The Chinese Journal of Geological Hazard and Control, 2022, 33（5）: 119 − 127. （in Chinese with English abstract）
[15] 杨志昊,郭永刚. 基于GIS的藏东南泥石流灾害风险预警系统研发[J]. 农业与技术,2024,44(6):91 − 96. [YANG Zhihao,GUO Yonggang. Research and development of risk early warning system for debris flow disaster in southeast Tibet based on GIS[J]. Agriculture and Technology,2024,44(6):91 − 96. (in Chinese)] YANG Zhihao, GUO Yonggang. Research and development of risk early warning system for debris flow disaster in southeast Tibet based on GIS［J］. Agriculture and Technology, 2024, 44（6）: 91 − 96. （in Chinese）
[16] 高瑞源. 基于机器学习方法的区域泥石流易发性评价研究[D]. 长春:吉林大学,2023. [GAO Ruiyuan. Evaluation of regional debris flow susceptibility based on machine learning method[D]. Changchun:Jilin University,2023. (in Chinese with English abstract)] GAO Ruiyuan. Evaluation of regional debris flow susceptibility based on machine learning method［D］. Changchun: Jilin University, 2023. （in Chinese with English abstract）
[17] 蒋庆丰. 山坡泥石流的风险性分析与评价研究[D]. 重庆:西南师范大学,2002. [JIANG Qingfeng. Risk analysis and evaluation of debris flow on hillside[D]. Chongqing:Southwest University,2002. (in Chinese with English abstract)] JIANG Qingfeng. Risk analysis and evaluation of debris flow on hillside［D］. Chongqing: Southwest University, 2002. （in Chinese with English abstract）
[18] 胡凯衡,张晓鹏,罗鸿,等. 丹巴县梅龙沟 “6•17” 泥石流灾害链调查[J]. 山地学报,2020,38(6):945 − 951. [HU Kaiheng,ZHANG Xiaopeng,LUO Hong,et al. Investigation of the “6•17” debris flow chain at the Meilong catchment of Danba County,China[J]. Mountain Research,2020,38(6):945 − 951. (in Chinese with English abstract)] HU Kaiheng, ZHANG Xiaopeng, LUO Hong, et al. Investigation of the “6•17” debris flow chain at the Meilong catchment of Danba County, China［J］. Mountain Research, 2020, 38（6）: 945 − 951. （in Chinese with English abstract）
[19] 崔鹏,庄建琦,陈兴长,等. 汶川地震区震后泥石流活动特征与防治对策[J]. 自然杂志,2010,32(1):19 − 25. [CUI Peng,ZHUANG Jianqi,CHEN Xingchang,et al. Activity characteristics and prevention countermeasures of debris flow after “5•12” Wenchuan devastating earthquake[J]. Nature Magazine,2010,32(1):19 − 25. (in Chinese with English abstract)] CUI Peng, ZHUANG Jianqi, CHEN Xingchang, et al. Activity characteristics and prevention countermeasures of debris flow after “5•12” Wenchuan devastating earthquake［J］. Nature Magazine, 2010, 32（1）: 19 − 25. （in Chinese with English abstract）
[20] 崔鹏,陈晓清,程尊兰,等. 西藏泥石流滑坡监测与防治[J]. 自然杂志,2010,32(1):19 − 25. [CUI Peng,CHEN Xiaoqing,CHENG Zunlan,et al. Monitoring and prevention of debris-flows and landslides in Tibet[J]. Chinese Journal of Nature,2010,32(1):19 − 25. (in Chinese with English abstract)] CUI Peng, CHEN Xiaoqing, CHENG Zunlan, et al. Monitoring and prevention of debris-flows and landslides in Tibet［J］. Chinese Journal of Nature, 2010, 32（1）: 19 − 25. （in Chinese with English abstract）
[21] 崔鹏,杨坤,陈杰. 前期降雨对泥石流形成的贡献——以蒋家沟泥石流形成为例[J]. 中国水土保持科学,2003,1(1):11 − 15. [CUI Peng,YANG Kun,CHEN Jie. Relationship between occurrence of debris flow and antecedent precipitation:Taking the Jiangjia gully as an example[J]. Science of Soil and Water Conservation,2003,1(1):11 − 15. (in Chinese with English abstract)] CUI Peng, YANG Kun, CHEN Jie. Relationship between occurrence of debris flow and antecedent precipitation: Taking the Jiangjia gully as an example［J］. Science of Soil and Water Conservation, 2003, 1（1）: 11 − 15. （in Chinese with English abstract）
[22] 张浩韦,刘福臻,王军朝,等. 基于FLO-2D数值模拟的工布江达县城泥石流灾害危险性评价[J]. 地质力学学报,2022,28(2):306 − 318. [ZHANG Haowei,LIU Fuzhen,WANG Junchao,et al. Hazard assessment of debris flows in Kongpo Gyamda,Tibet based on FLO-2D numerical simulation[J]. Journal of Geomechanics,2022,28(2):306 − 318. (in Chinese with English abstract)] ZHANG Haowei, LIU Fuzhen, WANG Junchao, et al. Hazard assessment of debris flows in Kongpo Gyamda, Tibet based on FLO-2D numerical simulation［J］. Journal of Geomechanics, 2022, 28（2）: 306 − 318. （in Chinese with English abstract）
[23] 刘德玉,蒋兴波,李保强. 基于FLO-2D的矿山泥石流运动特征及危险性评价 ——以甘肃省岷县簸箕沟矿山泥石流为例[J]. 地质与资源,2022,31(5):693 − 699. [LIU Deyu,JIANG Xingbo,LI Baoqiang. Movement characteristics and risk assessment of mine debris flow based on FLO-2D simulation:A case study of the debris flow in Bojigou mine of Minxian County,Gansu Province[J]. Geology and Resources,2022,31(5):693 − 699. (in Chinese with English abstract)] LIU Deyu, JIANG Xingbo, LI Baoqiang. Movement characteristics and risk assessment of mine debris flow based on FLO-2D simulation: A case study of the debris flow in Bojigou mine of Minxian County, Gansu Province［J］. Geology and Resources, 2022, 31（5）: 693 − 699. （in Chinese with English abstract）
[24] 王高峰,陈宗良,毛佳睿,等. 不同工程情景和降雨频率工况下白龙江流域泥石流危险性评价[J]. 山东科技大学学报(自然科学版),2020,39(5):30 − 40. [WANG Gaofeng,CHEN Zongliang,MAO Jiarui,et al. Debris flow risk assessment in Bailong River Basin under different engineering scenarios and rainfall frequency conditions[J]. Journal of Shandong University of Science and Technology (Natural Science Edition),2020,39(5):30 − 40. (in Chinese with English abstract)] WANG Gaofeng, CHEN Zongliang, MAO Jiarui, et al. Debris flow risk assessment in Bailong River Basin under different engineering scenarios and rainfall frequency conditions［J］. Journal of Shandong University of Science and Technology （Natural Science Edition）, 2020, 39（5）: 30 − 40. （in Chinese with English abstract）
[25] 曹鹏,侯圣山,陈亮,等. 基于数值模拟的群发性泥石流危险性评价——以甘肃岷县麻路河流域为例[J]. 中国地质灾害与防治学报,2021,32(2):100 − 109. [CAO Peng,HOU Shengshan,CHEN Liang,et al. Risk assessment of mass debris flow based on numerical simulation:An example from the Malu River Basin in Min County[J]. The Chinese Journal of Geological Hazard and Control,2021,32(2):100 − 109. (in Chinese with English abstract)] CAO Peng, HOU Shengshan, CHEN Liang, et al. Risk assessment of mass debris flow based on numerical simulation: An example from the Malu River Basin in Min County［J］. The Chinese Journal of Geological Hazard and Control, 2021, 32（2）: 100 − 109. （in Chinese with English abstract）
[26] 常鸣,窦向阳,唐川,等. 降雨驱动泥石流危险性评价[J]. 地球科学,2019,44(8):2794 − 2802. [CHANG Ming,DOU Xiangyang,TANG Chuan,et al. Hazard assessment of typical debris flow induced by rainfall intensity[J]. Earth Science,2019,44(8):2794 − 2802. (in Chinese with English abstract)] CHANG Ming, DOU Xiangyang, TANG Chuan, et al. Hazard assessment of typical debris flow induced by rainfall intensity［J］. Earth Science, 2019, 44（8）: 2794 − 2802. （in Chinese with English abstract）
[27] 王晓娅,陈晓文,焦姗,等. 不同尺度的忻州市东北部泥石流风险评价及风险管控研究[J]. 地质找矿论丛,2023,38(4):541 − 546. [WANG Xiaoya,CHEN Xiaowen,JIAO Shan,et al. Research on different scale risk assessment and management and control of mudflows in the northeast area of Xinzhou City[J]. Contributions to Geology and Mineral Resources Research,2023,38(4):541 − 546. (in Chinese with English abstract)] WANG Xiaoya, CHEN Xiaowen, JIAO Shan, et al. Research on different scale risk assessment and management and control of mudflows in the northeast area of Xinzhou City［J］. Contributions to Geology and Mineral Resources Research, 2023, 38（4）: 541 − 546. （in Chinese with English abstract）
[28] 赵聪,梁京涛,铁永波,等. 西藏雅鲁藏布江峡谷特大巨型泥石流活动与泥沙输移特征研究[J]. 中国地质灾害与防治学报,2024,35(4):45 − 55. [ZHAO Cong,LIANG Jingtao,TIE Yongbo,et al. Study on the activities of the massive debris flows and sediment transport characteristics in the Grand Bend of the Yarlung Zangbo River Gorge, Xizang[J]. The Chinese Journal of Geological Hazard and Control,2024,35(4):45 − 55. (in Chinese with English abstract)] ZHAO Cong, LIANG Jingtao, TIE Yongbo, et al. Study on the activities of the massive debris flows and sediment transport characteristics in the Grand Bend of the Yarlung Zangbo River Gorge, Xizang［J］. The Chinese Journal of Geological Hazard and Control, 2024, 35（4）: 45 − 55. （in Chinese with English abstract）
[29] 渠敬凯,杨为民,申俊峰,等. 基于DEM数据预处理的小流域泥石流防治工程效果分析——以漳县方家沟泥石流为例[J]. 水文地质工程地质,2024,51(4):206 − 219. [QU Jingkai,YANG Weimin,SHEN Junfeng,et al. Analysis of debris flow engineering treatment effect based on DEM data preprocessing in small watershed: A case study of Fangjiagou debris flow, Zhangxian County[J]. Hydrogeology & Engineering Geology,2024,51(4):206 − 219. (in Chinese with English abstract)] QU Jingkai, YANG Weimin, SHEN Junfeng, et al. Analysis of debris flow engineering treatment effect based on DEM data preprocessing in small watershed: A case study of Fangjiagou debris flow, Zhangxian County［J］. Hydrogeology & Engineering Geology, 2024, 51（4）: 206 − 219. （in Chinese with English abstract）
[30] 刘鑫,张文,李根,等. 高位远程崩滑碎屑流-泥石流灾害链的演变过程与影响范围预测——以“4•5” 四川洪雅县铁匠湾地质灾害链为例[J]. 吉林大学学报(地球科学版),2023,53(6):1799 − 1811. [LIU Xin,ZHANG Wen,LI Gen,et al. Research on evolution process and impact range prediction of high level remote collapse and landslide-debris flow disaster chain:Taking the “4•5” tiejiangwan geological disaster chain in Hongya County, Sichuan Province as an example[J]. Journal of Jilin University (Earth Science Edition),2023,53(6):1799 − 1811. (in Chinese with English abstract)] LIU Xin, ZHANG Wen, LI Gen, et al. Research on evolution process and impact range prediction of high level remote collapse and landslide-debris flow disaster chain: Taking the “4•5” tiejiangwan geological disaster chain in Hongya County, Sichuan Province as an example［J］. Journal of Jilin University （Earth Science Edition）, 2023, 53（6）: 1799 − 1811. （in Chinese with English abstract）
[31] 泥石流灾害防治工程勘查规范:DZ/T 0220-2006[S]. 2006. [Specification of geological investigation for debris flow stabilization:DZ/T 0220-2006[S]. 2006. (in Chinese)] Specification of geological investigation for debris flow stabilization: DZ/T 0220-2006［S］. 2006. （in Chinese）
[32] 康晓波,杨迎冬,王宇,等. 云南省地质灾害综合防治体系建设系列专项研究进展[J]. 中国地质灾害与防治学报,2023,34(6):146 − 157. [KANG Xiaobo,YANG Yingdong,WANG Yu,et al. Progress of the special-subjects study on the construction of comprehensive geological disaster prevention and control system in Yunnan Province[J]. The Chinese Journal of Geological Hazard and Control,2023,34(6):146 − 157. (in Chinese with English abstract)] KANG Xiaobo, YANG Yingdong, WANG Yu, et al. Progress of the special-subjects study on the construction of comprehensive geological disaster prevention and control system in Yunnan Province［J］. The Chinese Journal of Geological Hazard and Control, 2023, 34（6）: 146 − 157. （in Chinese with English abstract）
[33] 张宪政,铁永波,宁志杰,等. 四川汶川县板子沟 “6•26” 特大型泥石流成因特征与活动性研究[J]. 水文地质工程地质,2023,50(5):134 − 145. [ZHANG Xianzheng,TIE Yongbo,NING Zhijie,et al. Characteristics and activity analysis of the catastrophic“6•26” debris flow in the Banzi catchment, Wenchuan County of Sichuan Province[J]. Hydrogeology & Engineering Geology,2023,50(5):134 − 145. (in Chinese with English abstract)] ZHANG Xianzheng, TIE Yongbo, NING Zhijie, et al. Characteristics and activity analysis of the catastrophic“6•26” debris flow in the Banzi catchment, Wenchuan County of Sichuan Province［J］. Hydrogeology & Engineering Geology, 2023, 50（5）: 134 − 145. （in Chinese with English abstract）
[34] 万飞鹏,杨为民, 邱占林,等. 甘肃岷县纳古呢沟滑坡-泥石流灾害链成灾机制及其演化[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. (in Chinese with English abstract)] 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. （in Chinese with English abstract）
计量
- 文章访问数: 131
- HTML全文浏览量: 12
- PDF下载量: 68
邮件订阅
RSS