雅鲁藏布江色东普沟崩滑-碎屑流过程模拟及运动特征分析

李昆仲; 张明哲; 邢爱国

doi:10.16031/j.cnki.issn.1003-8035.2021.01.03

雅鲁藏布江色东普沟崩滑-碎屑流过程模拟及运动特征分析

1.
西藏自治区地质环境监测总站，西藏拉萨　850001
2.
上海交通大学海洋工程国家重点实验室，上海　200240

基金项目: 国家重点研发计划课题(2018YFC1504804)；国家自然科学基金(41530639)

详细信息

作者简介:
李昆仲(1978-)，男，江苏东海县人，硕士，高级工程师，主要从事地质灾害研究。E-mail：910949179@qq.com

通讯作者:
张明哲(1997-)，男，山东泰安人，博士研究生，主要从事高速远程滑坡灾害研究。E-mail：mingzhezhang@sjtu.edu.cn

中图分类号: P642.21
计量
- 文章访问数: 803
- HTML全文浏览量: 445
- PDF下载量: 456
出版历程
- 收稿日期: 2020-09-17
- 修回日期: 2020-09-24
- 网络出版日期: 2021-02-25
- 刊出日期: 2021-02-25

Numerical runout modeling and dynamic analysis of the ice avalanche-debris flow in Sedongpu Basin along Yarlung Zangbo River in Tibet

1.
Geological Environmental Monitoring Station of Tibet Autonomous Region, Lhasa, Tibet　850001, China
2.
State Key Laboratory of Ocean Engineering, Shanghai Jiao Tong University, Shanghai　200240, China

摘要

摘要: 2018年10月17日，西藏自治区雅鲁藏布江西岸色东普沟发生冰崩岩崩，冲刷侵蚀下部冰积堆积体和沟床松散堆积物，形成碎屑流冲入雅鲁藏布江，堆积后形成堰塞坝堵塞河道，严重影响当地居民安全和社会经济发展。本文基于真实遥感影像建立三维数值模型，利用DAN^3D分析其动力学特征，选取Frictional-Voellmy复合模型反演得到碎屑流堆积特征、滑体运动速度和铲刮深度等动力特征参数，与实际勘察结果基本一致，具有较好的可靠性。在野外勘察基础上，结合高密度电阻率法揭示了色东普沟崩滑-碎屑流堆积特征，结果表明DAN^3D能够较好地模拟滑体的堆积形态和最大厚度，验证了DAN^3D模拟的有效性，研究结论及方法对于此类灾害的机理研究及未来防治工程提供了新的研究思路。
- 色东普流域 /
- 冰崩 /
- 碎屑流 /
- DAN^3D /
- 高密度电阻率法
Abstract: On October 17, 2018, a massive ice avalanche-debris flow occurred in the Sedongpu Basin on the left bank of the Yarlung Zangbo River in Tibet. The sliding mass scoured the glacial deposits and loose deposits in the gully bed and then transformed into a debris flow. Eventually, a barrier dam was formed and blocked the Yarlung Zangbo River, which seriously affected the safety of residents and social and economic development. In this paper, a three-dimensional numerical simulation of this event was performed based on remote sensing. A dynamic analysis (DAN^3D) model was utilized to analyze the dynamic characteristics of sliding mass, and the Frictional-Voellmy composite model was selected to reproduce debris flow accumulation characteristics, movement speed, and scraping depth. The simulation results are consistent with the actual survey results and have good reliability. Based on field investigation, the electrical resistivity tomography method was used to reveal the post-event accumulation characteristics of the ice avalanche-debris flow in Sedongpu gully. The results showed that DAN^3D could well simulate the accumulation shape and maximum thickness of sliding mass, which verified the effectiveness of DAN^3D simulation. The results of the study can help forecast similar events and provide new research ideas for the mechanism research of such disasters for future prevention and control projects.
- Sedongpu Basin /
- iceavalanche /
- debris flow /
- DAN^3D /
- electrical resistivity tomography method

HTML全文

0. 引言

地质灾害监测是地质灾害防治手段中最重要的一环，通过分析不同种类地质灾害的发育情况和影响因素，在现场有针对性的布设监测设备进行高精度实时监测。一旦灾体发生倾斜、振动等形变异常情况，监测数据达到或者超过预警阈值，那么就会触发应急预案，再对现场进行核实调查。然而，很多时候监测设备受到卫星信号等环境因素干扰时，易产生误报和漏报，浪费人力物力，造成不良影响。因此，需要不断有针对性地对设备预警阈值进行修正。

Izadi等^[1]提出了一种预警阈值的方法，其基于操作特性曲线并以减少误报和持续监测为目的。其中一些学者^[2-8]介绍了预警与预警阈值设计之间的关系，目的是使预警阈值能处在一个合适的区间。Sharma等^[9]分析了火灾报警器的最新趋势和发展，该研究指出，大多数误报都是由于安装不良、维护不完善等人为因素造成的。研究对当前和潜在的问题和挑战提出了批判性观点。Al-Dabbagh等^[10]专注于操作指标，可视化图和警报洪水作为警报系统的决策支持工具，引入了多层雷达图，以提高报警指标的可读性。Quiñones-Grueiro等^[11]提出了一种基于报警-洪水临界指数的复排序设计方法数据驱动的多模过程建模与监测技术，回顾了多模式流程的建模和现有的监控方法，还研究了建模方法和监测方法的系统聚类和表征。Hollifield等^[12]研究了提高报警系统性能和减少报警迫害的不同方法，在研究的滤波方法中，可以突出显示延迟定时器和死区，但是，在警报监测期间会发生一些延迟。Kondaveeti等^[13]提出了另一种具有接收器工作基于特性曲线的技术，用于设计死区或阈值以降低误报率（false alarm rate，FAR）和漏报率（missing alarm rate，MAR）。Cheng等^[14]引入了一种优化滤波器设计方法，以增强报警系统的性能，而不是经典的移动式滤波器。Yang等^[15]使用相关色图提出了一种新的可视化方法，用于说明不同警报之间的时间滞后相关性。Han等^[16]提出了一种通过优化多变量报警阈值来最小化多维空间中误报概率（FAR）和漏报概率（MAR）的方法。虽然这种方法可以减少误报的数量，但同时减少误报和漏报是不可能的。Han等^[16]预发送了一种结合FAR、MAR和相关性分析的多元报警阈值优化方法，以解决报警阈值优化问题，重新对滋扰报警进行补偿。Xu等^[17]引入了一种基于证据理论的在线变量阈值方法，在优化报警系统设计的同时，考虑了与过程变量相关的不确定性。Gao等^[18]提出了一种基于相关性一致性的多元报警阈值优化方法。该文采用Pearson相关分析方法获得过程数据的相关系数，并采用粒子群优化（PSO）算法对问题进行优化。为了提高现有基于阈值的报警系统的效率，Bahar-Gogani等^[19]利用均值变化点检测方法提出了一种扩展的自适应阈值报警系统。

尽管到目前为止已经形成了大量的研究，但据作者所知，没有一项已发表的作品通过使用优先级和门和半马尔可夫过程（SMP）的概念来考虑一种基于优先级和门以及半马尔可夫过程(variable thresholds and semi-markov processes, VTAS)的性能分析问题。本研究介绍了两种类型的VTAS（即可变阈值和带可变死区带的可变阈值），它们能够减少报警干扰并提高性能。随后，利用优先级和门的功能行为，在不同条件下对VTAS进行FAR和MAR评估，提出了一种新的方法。FAR和MAR的优先级和门模型通过创建它们的半马尔可夫模型来求解。然后通过蒙特卡罗模拟验证结果，并与其他一些现有的研究工作进行比较。为了便于报警系统的优化设计，利用遗传算法获取并分析了所提出的非高斯分布方法。这种系统化方法及其解决方案使报警系统设计人员能够更深入地了解这些系统的性能行为。

1. 问题分析

报警系统的性能主要取决于3个判断指标：FAR、MAR和平均报警延迟（AAD）。为此，将随机离散信号x(t)视为采样时间间隔为h的过程变量测量值及其相关报警触发点x_tp，如图1所示。

图 1 过程变量测量值的随机离散信号x(t)

Figure 1. Random discrete signal x(t) of process variable measured values

下载: 全尺寸图片幻灯片

根据一种广泛使用的延迟方法，如果假设数据x(t)的n个连续样本超过触发点，则会发出警报。假设正常数据是小于x_tp的x(t)样本，异常数据是超过x_tp的x(t)样本，如图2所示，该图可分为四个部分。

图 2 x(t)的正常、异常、假报警和错过报警部分分类

Figure 2. Classification of normal, abnormal, false alarm and missed alarm parts of x(t)

下载: 全尺寸图片幻灯片

对正常和异常数据进行分类后，可以得到它们的概率密度函数，如图3所示。

图 3 x(t)正常部分和异常部分的分离概率密度函数

Figure 3. Separation probability density function of x(t) normal part and abnormal part

下载: 全尺寸图片幻灯片

将测量值分为正常和异常两类后，误报和漏报概率可分别计算如下：

F A R = q_{1} = \int_{x_{t p}}^{+ \infty} q (x) d x

(1)

M A R = p_{2} = \int_{- \infty}^{x_{t p}} p (x) d x

(2)

式中：q(x)——测量正常部分的概率密度函数；

p(x)——测量异常部分的概率密度函数。

现在，假设作为测量或监控信号的x(t)在时间t₀时等于或大于x_tp，在t_a时报警。报警延迟可定义如下：

T_{d} = t_{a} - t_{0}

(3)

由于信号x(t)是一个离散的随机变量，因此报警延迟T_d是一个随机变量。因此，平均报警延迟（AAD）定义为T_d的预望值。

{\overset{―}{T}}_{d} = E (T_{d})

(4)

式中：E——期望函数。

一般来说，我们需要x(t₀), ···, x(t₀+h)的多维联合概率密度函数来计算 ${\overset{―}{T}}_{d}$ 。为简单起见，假设x(t)是独立的且同分布的（IID），那么T_d=ih的概率密度函数可以写成：

P (T_{d} = i h) = p_{2}^{i} p_{1}

(5)

式中：P——T_d = ih的概率密度函数。

p₁=1−p₂，i=[1, 2, 3, ···]，平均报警延迟也是如此：

{\overset{―}{T}}_{d} = E (T_{d}) = \sum_{i = 0}^{\infty} i h p_{2}^{i} p_{1} = h \frac{p_{2}}{p_{1}}

(6)

2. 阈值报警系统背景

本节简要介绍了有死区和无死区的简单和可变阈值报警系统。作为一种定义，当数据的绝对量大于预定阈值的量时，会出现警报，否则会清除警报。式（7）给出了考虑零平均信号数据的报警的基本定义。为了简化计算，假设系统故障以上升而非间歇故障的形式发生，我们只有一个警报生成上限（自适应限值）。然而，通过所提出的方法，可以考虑故障与回落发生及其低限值。

A = {\begin{aligned} 1 & | S | ⩾ T \\ 0 & | S | < T \end{aligned}

(7)

式中：A——报警状态值；

T——阈值；

S——信号。

注意，该信号可以是基于信号的方法中的真实信号，也可以是基于模型的方法中的残余信号。

2.1 固定阈值报警系统

在这种类型的系统中，对信号应用固定阈值来诊断报警状态。可以通过以下方式实现：

T = m \pm φ v

(8)

式中：m——所有数据的平均值；

$φ$ ——可调指数；

v——标记信号中所有数据的方差值。

图4显示了信号x(t)应用的固定阈值。

图 4 信号应用固定阈值和死区

Figure 4. The generated signal has a static threshold and a deadband

下载: 全尺寸图片幻灯片

表1确定了FAR和MAR这两个用于分析报警管理系统运行的指标。必须确定适当的固定阈值，以尽可能同时降低这两个指标。当固定阈值足够大时，故障灵敏度降低。相反，由于真实数据中的不确定性和干扰行为，如果静态阈值过低，错误警报的数量会增加。在这种情况下，报警信号在开启和关闭状态之间传输，称为抖动报警。

表 1 可变阈值报警系统的结果

Table 1. Results of variable threshold alarm system

方法		固定阈值	固定阈值和死区	可变阈值	可变阈值和死区
指标	FAR	0.2668	0.2668	0.2009	0.2009
指标	MAR	0.2036	0.0452	0.1836	0.0448

下载: 导出CSV

| 显示表格

2.2 具有死区的简单阈值报警系统

引入死区方法，在存在干扰数据的情况下提高报警管理系统的性能。事实上，死区可以被确定为清除噪声信号中警报的另一个限制。根据阈值类型（上限阈值或下限阈值），选择合适的死区。死区通常定义为阈值或变量范围的百分比。由于变量的范围通常未知，因此在这项工作中，死区表示为报警限（阈值）的一部分。

d b w h = T + d b

(9)

d b w h = T - d b

(10)

式中：T——阈值；

db——死区值。

图4考虑了L=4和db=3.2。此外，表1中列出了以多种方式计算的FAR和MAR，以便进行比较。可以看出，死区限制的使用使得警报抖动比以前更好。然而，这种方法只是改进了MAR。这意味着发生的漏报更少，但这种方法不足以减少误报。

2.3 可变阈值报警系统

为了减少建模不确定性和测量噪声的影响，必须应用较大的固定阈值。然而事实上，这种解决方案可以大大减少误报。但是灵敏度也会降低。因此，这种方式是不合适的。此外，死区也不能大大减小。因此，为了使FAR最小化，引入了可变阈值。这种阈值的主要思想：干扰和其他非受控影响（如建模不确定性和测量噪声）随时间变化，阈值随时间变化。

在设计可变阈值时，定义了一个足够长的滑动窗口，并在每次滑动中计算窗口内数据的均值和方差。必须注意的是，根据信号的趋势，窗口大小可能会有所不同。此外，动量因子用于防止在可变阈值中忘记之前的数据样本。在以下等式给出了变量的计算过程。

\begin{aligned} \overset{―}{v} = γ v (k - 1) + (1 - γ) v (k) \\ \overset{―}{m} = γ m (k - 1) + (1 - γ) m (k) \\ T (k) = \overset{―}{m} (k) \pm α \overset{―}{v} (k) \end{aligned}

(11)

式中：v——滑动窗口内数据的方差；

γ——动量因子；

m——滑动窗口内数据的平均值；

k——窗口序号，k=1, 2, 3, ···；

α——可调因子；

T(k)——第k个窗口中的可变阈值。

这里的主要问题是适当地选择滑动窗口长度n。如果选择的n太小，阈值会很快适应由任何因素（例如干扰、噪声或故障）引起的残差变化。如果n太大，阈值的作用方式与常数相似，决策的敏感性降低。

此外，还考虑了对信号应用可变阈值的一些假设。在本文中，假设信号的正常状态可用，并且在这种情况下可以对信号应用可变阈值。然后，将确定的可变阈值应用于实际情况下的信号进行决策。让我们考虑上述假设，k=25，γ=0.5。然后，在应用可变阈值后，异常和正常情况下的数据概率密度函数可以如图5所示。此外，表1中也列出了以这种方式计算的FAR和MAR，并将其与固定阈值进行比较。

图 5 可变阈值a(x)、死区d(x)、正常q(x)和异常p(x)信号的估计概率密度函数

Figure 5. Estimated probability density function for variable thresholds, deadbands, normal and abnormal signals

下载: 全尺寸图片幻灯片

通过将该方法计算的FAR和MAR与前面两种方法进行比较，可以看出，与其他方法相比，该方法降低的FAR值较小。然而，与方法2（带死区的固定阈值）相比，MAR增大。为了克服这个缺点，我们引入了带可变死区的可变阈值。

2.4 具有可变死区的可变阈值报警系统

在该方法中，将可变阈值和死区进行合并，以降低FAR和MAR。如2.2节中方法所述，死区假设为报警限值的一部分。因此，由于阈值的变化，死区也应该变化。如图6所示，在信号上应用可变死区和可变阈值。与其他方法相比，该方法可以同时降低FAR和MAR。此功能是此方法相对于其他方法的优点。表1明确确认了这一问题。

图 6 具有可变阈值和死区的生成信号

Figure 6. Generated signal with adaptive threshold and deadband

下载: 全尺寸图片幻灯片

3. VTAS方法

显然，在正常情况下，当第i次测量的值大于阈值（固定阈值中的x_tp）时，会发出警报，称为假警报。类似地，在异常情况下，当第i个测量值小于阈值时，警报恢复正常，称为错过警报。也就是说，在正常情况下，当x_q>x_tp和在异常情况下发生虚警，当x_tp>x_P时发生漏报，其中x_q表示正常情况下的测量，x_p表示异常情况下的测量，x_tp是阈值。为了计算FAR和MAR，需要一种系统的方法。

动态故障树是一种流行的技术，它有助于对概率和随机现象进行系统建模和说明。在动态故障树中，优先级和门用于顺序故障行为的建模。如果两个事件A和B输入至优先级和门中，则当事件A发生在事件B发生之前时，输出将为真（t_A<t_B）。动态故障树中的这种情况如图7（a）所示。在本文中，利用优先级和门的行为进行FAR和MAR评估。在本文的方法中，如果在x域而不是在时域中使用优先级和门，则将有一个系统的概率门，其中输入是正常测量值（q）、异常测量值（p）和可变阈值（a）的累积分布函数（cumulative distribution function，CDF）。为了实现FAR，x_q>x_tp是正常条件下的对象，并且可以在图7（b）中示出。同样地，为了获得MAR，x_tp>x_p是处于异常条件下的对象，并且可以如图7（c）所描绘。通过使用优先级和门，可以得到一种系统的FAR和MAR计算方法。如前所述，假设测量值的显著增加是异常行为的迹象。相反，如果显著下降被视为异常行为，则必须切换MAR和FAR模型中优先级和门的输入。在这种情况下，对于FAR模型，“输入A”必须分别定义为门的第一个输入，“输入Q”必须分别定义为门的第二个输入。对于MAR模型，“输入P”必须分别设置为门的第一个输入，“输入A”必须分别设置为门的第二个输入。

图 7 动态故障树计算图解

Figure 7. Calculation diagram for a priority

下载: 全尺寸图片幻灯片

在相同的情况下，对于具有死区的VTAS，FAR模型可以如图8（a）所示，MAR模型可以如图8（b）所示。在这些模型中，假设警报将随上升故障而上升。因此，对于下降故障，必须按照前面的解释重建这些模型。

图 8 带死区的优先级和门计算图解

Figure 8. Calculation diagrams for a priority and gate with deadband

下载: 全尺寸图片幻灯片

3.1 假设条件

一般来说，报警系统对环境影响、缺陷和运行模式等因素有着很大的依赖性，这使得报警系统的设计任务具有挑战性。在本文中，作出以下假设，以提高方法对实际用途的适用性。

首先将本文中的测量值假设为独立且同分布（IID）；其次所提出的方法可以对下降和上升故障类型及其相关报警进行建模，而且尖峰和短时间歇性故障将被忽略；最后假设在每个VTAS的设计过程中，正常测量和异常测量在开始时由专家进行分类。

3.2 SMP的解析解

序言中汇总并讨论了不同类型的解决方案，所有这些解决方案都可用于评估VTAS中的MAR和FAR。在本节中，将研究SMP的解析解，并结合优先级和门来计算漏报和误报的参数概率。SMP可以用不同的符号来建模。本文使用元组来进行建模，即(p, P(x), F(x))，其中：p是初始分布的向量，P(x)是条件转移概率矩阵，F(x)描述了测量值的分布函数矩阵；考虑到X_i(i = 0, 1, 2)作为一个随机变量，其次SMP由初始状态概率 $p (0) = [1, 0, \dots, 0]$ ，条件转移概率矩阵 $P (x) = [P_{i j} (x)]$ 的向量确定。其中，P(x)由下式计算：

P_{i j} (x) = P {X (x) = j, (0) = i}

(12)

P_{i j} (x) = δ_{i j} [1 - G_{i} (x)] + \sum \int_{0}^{x} P_{i j} (t - x) d Q_{i j} (x)

(13)

G_{i j} (x) = P {S_{1} ⩽ x, X_{0} = i} = \sum_{j = 1}^{i} Q_{i j} (x)

(14)

Q_{i j} (x) = P {X_{1} = j, S_{1} ⩽ x, X_{0} = i}

(15)

其中，Q是每个状态的分布矩阵，G为对角矩阵；如果i=j，则δ_ij=1，否则，则δ_ij=0；X是由x组成的向量，i, j=[0, 1, 2, ···]。S是系统在不同x下的状态。

通过应用拉普拉斯-斯蒂尔杰斯变换可以得到 $Q_{i j}$ 的解。因此可以得出：

c {\tilde{p}}_{i j} (s) = δ_{i j} [1 - {\tilde{g}}_{i j} (s)] + {\tilde{p}}_{s} {\tilde{q}}_{i j} (s) {\tilde{p}}_{i j} (s)

(16)

式中：c——常数；

$\tilde{g}$ ——G在时间s的值；

$\tilde{q}$ 、 $\tilde{p}$ ——与q、p互补。

矩阵形式的等式（13）可以重写如下：

\tilde{p} (s) = (\sum_{n = 0}^{\infty} \tilde{q} (s)^{n}) (I - \tilde{g} (s))

(17)

x域中的无条件状态概率矩阵确定如下：

P (x) = P (0) P (x)

(18)

最后，通过总结某些状态的瞬态概率，可以得出FAR或MAR，这将在下一小节中进行详细讨论。

3.3 指数概率密度函数测量的FAR和MAR评估

为了简化所提方法的演示，本文考虑了异常测量和正常测量的指数概率密度函数。在优先级和门的半马尔可夫模型中，F_A(x)是优先级和门的第二个输入的累积分布函数（可变阈值的累积分布函数），F_Q(x)是优先级和门的第一个输入的累积分布函数。分别以从上到下和从左到右进行编号。该模型可通过上一小节中描述的SMP定理求解。核矩阵可以写成如式（18）的形式，每个状态的分布矩阵是式（19）的形式。

Q (x) = [\begin{array}{ccccc} 0 & Q_{1, 2} & Q_{1, 4} & 0 & 0 \\ 0 & 0 & 0 & Q_{2, 3} & 0 \\ 0 & 0 & 0 & 0 & Q_{4, 5} \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \end{array}]

(19)

\begin{aligned} \boldsymbol G (x) = & d i a g (1 - (1 - F_{A} (x)) (1 - F_{Q} (x)), F_{A} (x) \\ F_{Q} (x), δ (x - \infty), δ (x - \infty)) \end{aligned}

(20)

式中：diag——创建对角矩阵的函数；

$δ$ ——狄拉克函数。

\begin{aligned} Q_{1, 2} & = \int_{0}^{x} [1 - (1 - e^{- λ_{A} τ})] d {1 - e^{- λ_{Q} τ}} \\ = \int_{0}^{x} e^{- (λ_{A} + λ_{Q}) τ} λ_{Q} d τ = \frac{λ_{Q} [1 - e^{- (λ_{A} + λ_{Q}) x}]}{λ_{A} + λ_{Q}} \end{aligned}

(21)

\begin{aligned} Q_{1, 4} & = \int_{0}^{x} [1 - (1 - e^{- λ_{Q} τ})] d {1 - e^{- λ_{A} τ}} \\ = \int_{0}^{x} e^{- (λ_{A} + λ_{Q}) τ} λ_{A} d τ = \frac{λ_{A} [1 - e^{- (λ_{A} + λ_{Q}) x}]}{λ_{A} + λ_{Q}} \end{aligned}

(22)

Q_{2, 3} = \int_{0}^{x} λ_{A} e^{- λ_{A} τ} d τ = 1 - e^{- λ_{A} x}

(23)

Q_{4, 5} = \int_{0}^{x} λ_{Q} e^{- λ_{Q} τ} d τ = 1 - e^{- λ_{Q} x}

(24)

式中：λ——门的输入系数；

τ——输出系数。

可使用式（23）从失效状态的概率获得优先级和门输出的概率。

F A R = q (x) = T^{- 1} {q_{1, 4} (s) (f_{Q}^{*} (s))}

(25)

同样，MAR的概率将得到：

M A R = P (x) = T^{- 1} {q_{1, 4} (s) (f_{P}^{*} (s))}

(26)

式中：T——阈值；

s——时间/s；

$f_{Q}^{*}$ ——0到 $+ \infty$ 的积分；

$f_{P}^{*}$ ——0到 $- \infty$ 的积分。

3.4 蒙特卡罗模拟验证

通过一千万次重复的蒙特卡罗模拟验证了所提出的半马尔可夫方法，并且在每次迭代中，生成的混淆矩阵的假正值和负值都存储在向量中。仿真完成后，将假阳性和假阴性的均值和方差值分别计算为FAR和MAR。表2提供了VTAS解决方案和蒙特卡罗模拟之间的比较。根据表2，验证了VTAS方法的正确性。

表 2 蒙特卡罗模拟与VTAS方法的比较

Table 2. Comparison between Monte Carlo simulation and VTAS

方法		蒙特卡罗模拟	VTAS的解
指标	FAR	0.159852	0.15867
指标	MAR	0.159856	0.15867

下载: 导出CSV

| 显示表格

3.5 遗传算法优化

在设计最优报警系统时，可以考虑多目标进化算法。遗传算法是求解优化问题的著名元启发式算法之一。遗传算法的第一步是生成随机解的初始总体。下一个步骤是选择操作，其中使用选择函数从当前总体中选择两个或多个解决方案。然后，对这些解执行交叉操作以生成新解。然后以给定的变异概率对新解进行变异。

在这种方法中，输入变量（正常和异常测量中的样本警报数（n, m）、样本窗口宽度（ω）和VTAS中的α被假定为基因型，表型将是FAR、MAR和平均警报检测。输入的分布服从一致函数，其中n和m的低值为1，高值为20。α参数的范围为20到30。此外，窗户的长度在20到100之间。选择每代成本最低的顶级基因作为交叉函数的输入，生成下一个群体。同时，突变函数生成新的随机基因并将其添加到下一个群体中。此过程将重复一定次数的迭代，以达到最佳值（最低成本）。

在优化过程中，每代考虑40个种群；突变率为35%，单点交叉率为50%，重组率为15%。报警系统设计、报警延迟和AAD取决于一些报警样本（n−m）增加n或m会导致AAD值增加。此外，无论α和w对MAR和FAR的影响如何，在大多数情况下，增加报警样本数n，m会导致MAR和FAR值降低。为了在AAD、MAR和FAR的优化之间进行权衡，加权和成本函数可用作 $J (α, w, n, m)$ 。

J (α, w, n, m) = ω_{1} \frac{F A R}{I_{R F A R}} + ω_{2} \frac{M A R}{I_{R M A R}} + ω_{3} \frac{A A D}{I_{R A A D}}

(27)

式中：I_RFAR、I_RMAR和I_RAAD——FAR、MAR和AAD的补数；

$ω_{1}$ 、 $ω_{2}$ 、 $ω_{3}$ ——权重系数。

注意，由于nm样本报警延迟的不同延迟，AAD被分为两部分，分别用于n和m样本的延迟。根据实际情况，每个FAR、MAR和AAD的重要性可分别通过ω₁、ω₂、ω₃进行加权。引入的权重可以根据应用程序进行不同的选择。例如，如果在假设的系统中，延迟减少是最重要的问题，且MAR优先于FAR，则成本函数J(α, w, m, n)的权重应视ω₁<ω₂<ω₃。作为遗传算法（genetic algorithm，GA）优化的结果，可以实现“报警样本数（n，m）”、“自适应阈值的α系数”和“自适应阈值的窗口大小w”的最佳值。换句话说，优化问题的目标可以写为：

(α, w, n, m) = \arg min J (α, w, n, m)

(28)

下面，通过一些案例研究，说明了所提出的方法和GA优化的有效性。

4. 仿真

在本节中，提供了2个实例来说明基于半马尔可夫解的能力，并将所提出方法的有效性与其他已发表的研究进行比较。

4.1 实例1：威布尔分布测量

在前面的章节中，我们研究了一个具有分段白高斯随机分布的测量示例，并给出了这个示例来说明所提出的方法在非高斯随机分布中的有效性。因此，威布尔分布采用以下表达式考虑。假设生成的数据长度为200，其中0到100的测量值服从F_normal，101到200的测量值遵循F_abnormal。

F_{normal} (x) = 1 - e^{- {(\frac{x}{A})}^{B}}, A = 1, B = 1.3

(29)

F_{abnormal} (x) = 1 - e^{- {(\frac{x}{A})}^{B}}, A = 5, B = 2

(30)

式中：F_normal——正常的概率密度；

F_abnormal——异常的概率密度。

正常和异常假设测量的概率密度函数曲线如图9所示。表3分别显示了通过蒙特卡罗模拟、拟议方法和简单马尔可夫解计算的FAR和MAR。从表3中可以看出，所提出的方法和蒙特卡罗模拟的结果是相同的。然而，有两种方法可以使用简单的马尔可夫解：①计算正常和异常测量的平均值和方差，并使用分析正态分布表达式；②使用数值概率密度函数估计。由于测量值的非高斯行为，本例中基于马尔可夫的简单解的结果是不精确的。

图 9 测量的概率密度函数（威布尔分布）

Figure 9. Probability density function (Weibull distribution) of measurements

下载: 全尺寸图片幻灯片

表 3 所提出的解、蒙特卡罗结果和简单马尔可夫解之间的比较

Table 3. Comparison between the proposed solution, Monte Carlo results and simple Markov solutions

方法		VTAS的解	蒙特卡罗结果	简单马尔可夫解
指标	FAR	0.1371	0.1371	0.1455
指标	MAR	0.1371	0.1371	0.1070

下载: 导出CSV

| 显示表格

4.2 实例2：高斯分布的测量

假设生成的测量值x(t)为“分段白高斯随机变量”，不精确的平均值和方差为下列公式。随机数生成的采样率为1 s，x(t)如图10所示。

图 10 生成的随机数

Figure 10. The distribution of the generated random number over time

下载: 全尺寸图片幻灯片

{\begin{aligned} x (t) = N (μ_{1}, σ_{1}^{2}) {\begin{matrix} μ_{1} = U (0.2, 0 .3) \\ σ_{1} = U (1.5, 1.6) \end{matrix}, t < 2 000 \\ x (t) = N (μ_{2}, σ_{2}^{2}) {\begin{matrix} μ_{2} = U (1.2, 1.5) \\ σ_{2} = U (1.5, 1.6) \end{matrix}, t ⩾ 2 000 \end{aligned}

(31)

式中：U——均匀分布的随机变量；

t——时间/s；

μ——期望值；

$σ$ ——方差；

N——正态分布。

在本例中，假设测量的静态分布部分已知，且不精确，因为x(t)混合了认知和假设不确定性。如图10所示，这是简单阈值决策的挑战性示例。然而，通过使用一个采样延迟和权重系数为0.05的指数加权移动平均（EWMA）滤波器，可以从图10中获得图11。很明显，EWMA过滤器的使用简化了决策过程，并改善了MAR和FAR。

图 11 使用具有自适应阈值的EWMA过滤器（alpha=0.5，窗口大小=100）后生成的随机数

Figure 11. Generated random numbers after filtering with an EWMA filter （alpha=0.5, window size =100） with adaptive threshold

下载: 全尺寸图片幻灯片

表4提供了可变阈值与其他现有方法之间性能指标（MAR和FAR）的比较结果，如固定阈值（ST）和指数加权移动平均EWMA报警系统、循证报警系统evidence-based alarm system和三阶移动平均滤波器（3OMAF）、三样本报警延迟计时器（3SADT）。根据这些结果，使用EWMA滤波器的VTAS比其他方法具有更好的性能指标。

表 4 可变阈值与其他方法性能指标的比较结果

Table 4. Comparison results of performance metrics of variable threshold method with other methods

方法		ST	VTAS	EWMA	VTAS （+EWMA filter）	Evidence-based alarm system	3OMAF	3SADT
指标	MAR	0.3519	0.3118	0.0437	0.0518	0.0582	0.2470	0.1511
指标	FAR	0.3806	0.2860	0.0447	0.0241	0.0429	0.2926	0.2527

下载: 导出CSV

| 显示表格

5. 应用

将上述方法应用于地质灾害监测预警系统中某实时监测点位上进行实测。

图12展示出了传感器的输出信号，其中曲线的绿色表示异常状态，蓝色表示正常状态。在图12中，设计的自适应阈值用红色表示。假设系统中不存在持续时间很短的间歇性故障。因此，正常信号（蓝色）中观察到的峰值不被视为异常事件或情况。它们可能是间歇性故障、高振幅噪音或干扰。然而，即使这些峰值是异常的，根据我们的假设，我们也必须忽略它们。

图 12 用自适应阈值（α=20）测量传感器（时间（s）-振动测量x(t) mm/s）

Figure 12. Measuring sensors with adaptive threshold (α = 20)

下载: 全尺寸图片幻灯片

在图13中，正常/异常条件下测量的概率密度函数和设计的自适应阈值分别用绿色、蓝色和红色表示。此图显示了可变阈值的估计概率密度函数遵循与正常数据的概率密度函数类似的模式，正如所预期的那样。

图 13 传感器在正常和异常情况下的概率密度函数（测量范围-概率）

Figure 13. Probability density function of sensor under normal and abnormal conditions

下载: 全尺寸图片幻灯片

表5描述了几种报警系统性能评估的比较结果。可以看出，通过使用带死区的固定阈值，FAR增加MAR降低。此外，在可变阈值带死区系统中，FAR增加。然而，可变阈值带死区系统可以比其他系统更好地降低MAR。根据GA优化部分，不同的需求导致不同的优化结果。根据优先级和确定的成本函数系数，将获得VTAS的不同值。如表6所示，当允许使用样本延迟时，AAD的值将增加，MAR和FAR将按预期降低。

表 5 不同方法性能指标的比较结果

Table 5. Comparison results of different methods’performance metrics

方法		固定阈值	可变阈值（带死区）	固定阈值（带死区）
指标	FAR	0.192783	0.289527	0.228910
指标	MAR	0.460713	0.207569	0.364654

下载: 导出CSV

| 显示表格

表 6 用遗传算法优化VTAS

Table 6. Optimizing VTAS with genetic algorithm

遗传算法中损失函数的权重			报警系统参数（α, ω, n, m）	MAR	FAR	AAD
ω₃	ω₂	ω₁	（25.024, 82, 1, 1）	0.2109	0.2282	1.2674
1	1	10
I_RAAD	I_RFAR	I_RMAR
1	1e−05	1e−05
ω₃	ω₂	ω₁	（25.039, 82, 2, 3）	0.0556	0.1681	2.8742
100	0.1	0.1
I_RAAD	I_RFAR	I_RMAR
1	1e−03	1e−03
ω₃	ω₂	ω₁	（25.680, 83, 2, 3）	0.0596	0.1625	2.8979
100	0.5	0.1
I_RAAD	I_RFAR	I_RMAR
1	1e−03	1e−03

下载: 导出CSV

| 显示表格

首先选择传统的固定阈值方法计算得到固定的预警阈值，其结果如图14所示，由于系统数据跳变触发阈值，产生很多的误报警，使得报警数量增多。而采用本文的动态阈值估算方法，通过选择180 d时间序列的监测值,确定滑动窗口的长度为200，如图15所示。通过阈值的不断调整，得到如图16所示的结果。并结合表7和表8可以看出，采用预警数量得到了进一步的减少，并且预警性能指标得到明显改善。

图 14 固定预警阈值

Figure 14. Fixed threshold alert threshold

下载: 全尺寸图片幻灯片

图 15 滑动窗口长度

Figure 15. Sliding window length

下载: 全尺寸图片幻灯片

图 16 可变预警阈值

Figure 16. Dynamic alert threshold

下载: 全尺寸图片幻灯片

表 7 预警次数对比

Table 7. Comparison of warning times

方法	预警次数
固定阈值	301
可变阈值（带死区）	62

下载: 导出CSV

| 显示表格

表 8 预警性能指标对比

Table 8. Comparison of early warning performance indicators

方法	MAR	FAR
固定阈值	0.0830	0.137
可变阈值（带死区）	0.0007	0.013

下载: 导出CSV

| 显示表格

6. 结论

可变阈值报警系统广泛应用于许多行业，以减少报警干扰。在本研究中，首先介绍了报警系统的主要性能指标FAR和MAR。随后，将优先级和门概念用于MAR和FAR的建模。针对MAR和FAR模型的评估，提出了基于SMP的解析解，该解可扩展到有无可变死区的不同类型VTAS，并利用半马尔可夫过程对具有非高斯特性的报警系统进行了性能优化。此外，本研究将所提出的VTAS方法与遗传算法相结合，以获得最佳的参数选择。

本文通过仿真与相关研究工作进行比较，结果表明使用EWMA滤波器的VTAS比其他方法具有更好的性能指标，证明了该方法的有效性。基于此，本文还将提出的方法应用于实际案例中，结果表明，当使用固定阈值时，FAR为0.19，MAR为0.46，当应用本文提出的方法得到可变阈值后，FAR降低到0.29，MAR降低为0.21。以上结果表明，相比于使用固定阈值，可变阈值系统能有效地减少误报及漏报，增加预警的准确性，全面提升相关单位对突发性地质灾害的分析、预警、处置和服务的能力，为政府相关部门进行地质环境与地质灾害决策管理和社会服务提供技术保障。

图 1 色东普沟崩滑-碎屑流地理位置

a-滑体运动路径边界；b-崩滑区；c-流通区；d-堆积区红色箭头表示滑体的滑动方向

Figure 1. Location of the ice avalanche-debris flow in Sedongpu Basin

下载: 全尺寸图片幻灯片

图 2 色东普流域气温、降水变化曲线

Figure 2. Temperature and precipitation variation in Sedongpu Basion

下载: 全尺寸图片幻灯片

图 3 色东普沟主沟剖面图

Figure 3. Longitudinal profile of the main gully in Sedongpu Basin

下载: 全尺寸图片幻灯片

图 4 色东普沟崩滑-碎屑流数字高程模型

Figure 4. Digital elevation model of the ice avalanche-debris flow in Sedongpu Basin

下载: 全尺寸图片幻灯片

图 5 崩滑-碎屑流过程不同时刻滑坡堆积形态

Figure 5. Time-lapse image of deposit distribution during the ice avalanche-debris flow in Sedongpu Basin

下载: 全尺寸图片幻灯片

图 6 崩滑-碎屑流铲刮分布

Figure 6. Scraping distribution of the ice avalanche-debris flow

下载: 全尺寸图片幻灯片

图 7 崩滑-碎屑流最大速度分布图

Figure 7. Maximum velocity contour of theice avalanche-debris flow

下载: 全尺寸图片幻灯片

图 8 崩滑-碎屑流物探测线分布

Figure 8. ERT survey lines distribution of the ice avalanche-debris flow

下载: 全尺寸图片幻灯片

图 9 横向测线(ERT4-ERT7)物探成果图

Figure 9. Transverse ERT profiles along the lines ERT4 to ERT7

下载: 全尺寸图片幻灯片

图 10 纵向测线(ERT1-ERT3)物探成果图

Figure 10. Longitudinal ERT profiles along the lines ERT1 to ERT3

下载: 全尺寸图片幻灯片

图 11 物探解译和DAN^3D模拟堆积物厚度对比

Figure 11. Comparison of the landslide deposits depth from the ERT interpretation and DAN^3D

下载: 全尺寸图片幻灯片

表 1 流变模型和计算参数

Table 1 Rheological model and simulation parameters

	模型	容重/ (kN·m⁻³)	摩擦角/ (°)	摩擦系数	湍流系数/ m	内摩擦角/ (°)	最大铲刮深度/m
崩滑区	Frictional	25	28	−	−	36	0
流通区	Voellmy	21	−	0.04	1200	24	35
堆积区	Voellmy	21	−	0.05	1000	24	17

下载: 导出CSV

表 2 物探解译与DAN^3D模拟堆积物最大深度对比（单位: m）

Table 2 Comparison of the maximum deposition depth between the ERT method and DAN^3D simulation (Unit: m)

	ERT1	ERT2	ERT3	ERT4	ERT5	ERT6	ERT7
物探解译	29	40	19	12	15	20	18
DAN^3D模拟	30	25	20	15	15	10	10

下载: 导出CSV

参考文献(21)

[1]	SHANG YJ, YANG ZF, LI LH, et al. A super-large landslide in Tibet in 2000: background, occurrence, disaster, and origin[J]. Geomorphology,2003,53:225 − 243.
[2]	YIN YP, XING AG. Aerodynamic modeling of the Yigong gigantic rock slide-debris avalanche, Tibet, China[J]. Bulletin of Engineering Geology and the Environment,2012,71:149 − 160.
[3]	CUI P, JIA Y. Mountain hazards in the Tibetan Plateau: research status and prospects[J]. National Science Review,2015,2(4):397 − 402.
[4]	王宗太, 刘潮海. 中国冰川分布的地理特征[J]. 冰川冻土,2001,23(3):231 − 237. [WANG Zongtai, LIU Chaohai. Geographical characteristics of the distribution of glaciers in China[J]. Journal of Glaciology and Geocryology,2001,23(3):231 − 237. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-0240.2001.03.004
[5]	崔鹏, 贾洋, 苏凤环, 等. 青藏高原自然灾害发育现状与未来关注的科学问题[J]. 中国科学院院刊,2017,32(9):985 − 992. [CUI Peng, JIA Yang, SU Fenghuan, et al. Natural hazards in Tibetan Plateau and key issue for feature research[J]. Bulletin of Chinese Academy of Sciences,2017,32(9):985 − 992. (in Chinese with English abstract)
[6]	曾庆利, 杨志法, 尚彦军, 等. 川藏公路102滑坡群地段基覆界面问题研究[J]. 水文地质工程地质,2008,3:1 − 6. [ZENG Qingli, YANG Zhifa, SHANG Yanjun, et al. Study on discontinuity between bedrock and quaternary covers in No.102 landslide group section of Sichuan-Tibet highway[J]. Hydrogeology & Engineering Geology,2008,3:1 − 6. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-3665.2008.04.002
[7]	胡文涛, 姚檀栋, 余武生, 等. 高亚洲地区冰崩灾害的研究进展[J]. 冰川冻土,2018,40(6):1141 − 1152. [HU Wentao, YAO Tandong, YU Wusheng, et al. Advances in the study of glacier avalanches in High Asia[J]. Journal of Glaciology and Geocryology,2018,40(6):1141 − 1152. (in Chinese with English abstract)
[8]	许强, 郑光, 李为乐, 等. 2018年10月和11月金沙江白格两次滑坡-堰塞堵江事件分析研究[J]. 工程地质学报,2018,26(6):1534 − 1551. [XU Qiang, ZHENG Guang, LI Weile, et al. Study on successive landslide damming events of Jinsha River in Baige Village on Octorber 11 and November[J]. Joumal of Engineering Geology,2018,26(6):1534 − 1551. (in Chinese with English abstract)
[9]	刘传正, 吕杰堂, 童立强, 等. 雅鲁藏布江色东普沟崩滑-碎屑流堵江灾害初步研究[J]. 中国地质,2019,46(2):219 − 234. [LIU Chuanzheng, LYU Jietang, TONG Liqiang, et al. Research on glacial/rock fall-landslide-debris flows in Sedongpu basin along Yarlung Zangbo River in Tibet[J]. Geology in China,2019,46(2):219 − 234. (in Chinese with English abstract)
[10]	童立强, 涂杰楠, 裴丽鑫, 等. 雅鲁藏布江加拉白垒峰色东普流域频繁发生碎屑流事件初步探讨[J]. 工程地质学报,2018,26(6):1552 − 1561. [TONG Liqiang, TU Jienan, PEI Lixin, et al. Preliminary discussion of the frequently debris flow events in Sedongpu Basin at Gyalaperi peak, Yarlung Zangho River[J]. Journal of Engineering Geology,2018,26(6):1552 − 1561. (in Chinese with English abstract)
[11]	HU KH, ZHANG XP, YOU Y, et al. Landslides and dammed lakes triggered by the 2017 Ms6.9 Milin earthquake in the Tsangpo gorge[J]. Landslides,2019,16(2/3):1 − 9.
[12]	WANG WP, YANG JS, WANG YB. Dynamic processes of 2018 Sedongpu landslide in Namcha Barwa–Gyala Peri massif revealed by broadband seismic records[J]. Landslides,2020,17:409 − 418.
[13]	HUNGR O. A model for the runout analysis of rapid flow slides, debris flows and avalanches[J]. Canadian Geotechnical Journal,1995,32(4):610 − 623.
[14]	KEITH B. The 2000 Yigong landslide (Tibetan Plateau), rockslide-dammed lake and outburst flood: Review, remote sensing analysis, and process modelling[J]. Geomorphology,2015,246(1):377 − 393.
[15]	戴兴建, 殷跃平, 邢爱国. 易贡滑坡-碎屑流-堰塞坝溃坝链生灾害全过程模拟与动态特征分析[J]. 中国地质灾害与防治学报,2019,30(5):1 − 8. [DAI Xingjian, YIN Yueping, XING Aiguo. Simulation and dynamic analysis of Yigong rockslide-debris avalanche-dam breaking disaster chain[J]. The Chinese Journal of Geological Hazard and Control,2019,30(5):1 − 8. (in Chinese with English abstract)
[16]	王国章, 李滨, 冯振, 等. 重庆武隆鸡冠岭岩质崩滑-碎屑流过程模拟[J]. 水文地质工程地质,2014,41(5):101 − 106. [WANG Guozhang, LI Bin, FENG Zhen, et al. Simulation of the process of the Jiguanling rock avalanche in Wulong of Chongqing[J]. Hydrogeology & Engineering Geology,2014,41(5):101 − 106. (in Chinese with English abstract)
[17]	XING AG, YUAN XY, XU Q, et al. Characteristics and numerical runout modelling of a catastrophic rock avalanche triggered by the Wenchuan earthquake in the Wenjia valley, Mianzhu, Sichuan, China[J]. Landslides,2017,14(1):83 − 98.
[18]	XING AG, WANG GH, LI B, et al. Long-runout mechanism and landsliding behaviour of large catastrophic landslide triggered by heavy rainfall in Guanling, Guizhou, China[J]. Canadian Geotechnical Journal,2015,52(7):971 − 981.
[19]	张远娇, 邢爱国. 汶川地震触发牛圈沟高速远程滑坡-碎屑流动力学特性分析[J]. 上海交通大学学报,2012,46(10):1665 − 1670. [ZHANG Yuanjiao, XING Aiguo. Dynamics analysis of Niujuanhou rockslide-debris flow triggered by Wenchuan earthquake[J]. Journal of Shanghai Jiao Tong University,2012,46(10):1665 − 1670. (in Chinese with English abstract)
[20]	LOKE MH, BARKER RD. Rapid least-squares inversion of apparent resistivity pseudosections using a quasi-Newton method[J]. Geophysical prospecting,1996,44:131 − 152.
[21]	XINGAG, WANG GH, YIN YP, et al. Dynamic analysis and field investigation of a fluidized landslide in Guanling, Guizhou, China[J]. Engineering Geology,2014,181:1 − 14.

施引文献(3)

期刊类型引用(3)

1.	丁文海，王朋伟，安玉科，宋金光. 基于监测数据的抗滑桩变形机理与预警判据研究. 公路. 2024(08): 152-159 . 百度学术
2.	张文政，邢真强. 煤矿热动力灾害智能预警软件开发及系统集成. 山东煤炭科技. 2024(11): 74-78 . 百度学术
3.	肖自为，马源，李金林，龚弦，邓修林，程曦. 基于过程预警模型的矿山边坡监测系统研究与应用. 地质灾害与环境保护. 2024(04): 91-98 . 百度学术

其他类型引用(0)

资源附件(0)

图(12) / 表(2)

计量

文章访问数: 803
HTML全文浏览量: 445
PDF下载量: 456
被引次数: 3

0. 引言
1. 问题分析
2. 阈值报警系统背景
2.1 固定阈值报警系统
2.2 具有死区的简单阈值报警系统
2.3 可变阈值报警系统
2.4 具有可变死区的可变阈值报警系统
3. VTAS方法
3.1 假设条件
3.2 SMP的解析解
3.3 指数概率密度函数测量的FAR和MAR评估
3.4 蒙特卡罗模拟验证
3.5 遗传算法优化
4. 仿真
4.1 实例1：威布尔分布测量
4.2 实例2：高斯分布的测量
5. 应用
6. 结论

0. 引言
1. 问题分析
2. 阈值报警系统背景
2.1 固定阈值报警系统
2.2 具有死区的简单阈值报警系统
2.3 可变阈值报警系统
2.4 具有可变死区的可变阈值报警系统
3. VTAS方法
3.1 假设条件
3.2 SMP的解析解
3.3 指数概率密度函数测量的FAR和MAR评估
3.4 蒙特卡罗模拟验证
3.5 遗传算法优化
4. 仿真
4.1 实例1：威布尔分布测量
4.2 实例2：高斯分布的测量
5. 应用
6. 结论

参考文献(21)

施引文献

资源附件(0)

雅鲁藏布江色东普沟崩滑-碎屑流过程模拟及运动特征分析

作者简介: 李昆仲(1978-)，男，江苏东海县人，硕士，高级工程师，主要从事地质灾害研究。E-mail：910949179@qq.com

通讯作者: 张明哲(1997-)，男，山东泰安人，博士研究生，主要从事高速远程滑坡灾害研究。E-mail：mingzhezhang@sjtu.edu.cn

计量

出版历程

Numerical runout modeling and dynamic analysis of the ice avalanche-debris flow in Sedongpu Basin along Yarlung Zangbo River in Tibet

0. 引言

1. 问题分析

2. 阈值报警系统背景

2.1 固定阈值报警系统

2.2 具有死区的简单阈值报警系统

2.3 可变阈值报警系统

2.4 具有可变死区的可变阈值报警系统

3. VTAS方法

3.1 假设条件

3.2 SMP的解析解

3.3 指数概率密度函数测量的FAR和MAR评估

3.4 蒙特卡罗模拟验证

3.5 遗传算法优化

4. 仿真

4.1 实例1：威布尔分布测量

4.2 实例2：高斯分布的测量

5. 应用

6. 结论

期刊类型引用(3)

其他类型引用(0)

计量

出版历程

目录

0. 引言

1. 问题分析

2. 阈值报警系统背景

2.1 固定阈值报警系统

2.2 具有死区的简单阈值报警系统

2.3 可变阈值报警系统

2.4 具有可变死区的可变阈值报警系统

3. VTAS方法

3.1 假设条件

3.2 SMP的解析解

3.3 指数概率密度函数测量的FAR和MAR评估

3.4 蒙特卡罗模拟验证

3.5 遗传算法优化

4. 仿真

4.1 实例1：威布尔分布测量

4.2 实例2：高斯分布的测量

5. 应用

6. 结论

友情链接：

作者简介:
李昆仲(1978-)，男，江苏东海县人，硕士，高级工程师，主要从事地质灾害研究。E-mail：910949179@qq.com

通讯作者:
张明哲(1997-)，男，山东泰安人，博士研究生，主要从事高速远程滑坡灾害研究。E-mail：mingzhezhang@sjtu.edu.cn