CN112182993A - 滑坡堰塞坝溃口深度测算方法、滑坡堰塞坝溃决洪水峰值流量测算方法 - Google Patents

Abstract

本发明公开滑坡堰塞坝溃口深度测算方法、滑坡堰塞坝溃决洪水峰值流量测算方法。针对现有技术缺陷，本发明滑坡堰塞坝坝体水流易侵蚀性判断方法，以坝体物质中值粒径为指标，快速判断坝体水流易侵蚀程度。本发明滑坡堰塞坝溃口深度测算方法、滑坡堰塞坝溃口底宽度测算方法，将与坝体侵蚀特征有关的特定计算系数引入到测算方法中，能够充分考虑堰塞坝的水流易侵蚀性对溃口形成与演化的影响。本发明滑坡堰塞坝溃决洪水峰值流量测算方法，采用半解析公式计算溃决峰值流量，计算结果更接近灾害现场测量值。本发明技术方案测算原理科学、所需计算量小、数据采集与计算硬件要求低，能够充分满足滑坡堰塞坝形成与溃决事件中快速响应的防灾减灾措施需要。

Description

滑坡堰塞坝溃口深度测算方法、滑坡堰塞坝溃决洪水峰值流 量测算方法

技术领域

本发明涉及一种堰塞坝溃决特征指标的测算方法，特别是涉及一种滑坡堰塞坝溃口深度测算方法与一种滑坡堰塞坝溃决洪水峰值流量测算方法，属于地质灾害监测防治技术领域、地质灾害应急响应技术领域。

背景技术

河谷山体在崩塌、滑坡、泥石流等灾害的作用下，横向阻塞山谷、河谷或河床形成堰塞坝，导致上游段壅水形成堰塞湖。堰塞湖在世界各地广泛分布。统计数据显示，堰塞坝从形成至发生溃决往往历时较短，短则1天，长则1周。因而自堰塞坝形成开始，围绕其展开的救援治理工作都必须得到快速响应。堰塞坝溃决时，常形成汹涌洪水，引发较大洪水灾害。这使得快速科学合理计算堰塞坝溃决洪水峰值流量十分重要。

滑坡堰塞坝溃决洪水峰值流量是地质灾害危险性评估和工程防治中的关键参数。在滑坡堰塞坝形成后，需要尽快对自然溃坝的形式做出预测，即对溃口的特征值进行准确计算，主要溃口特征值为溃口深度和溃口底宽。在确定这两个特征值之后，即可依据相关公式计算出溃决洪水峰值流量、水位、流速及其在下游河道的演进过程，进而预测溃坝洪水在下游的危害范围和程度，为制定有效的防灾减灾预案提供科技支撑。

申请号为2019106536250、名称为“堰塞湖危险性分级和溃决型泥石流及山洪险情预警方法”的中国发明专利申请公开了一种堰塞湖溃决峰值流量计算方法，是基于堰塞坝溃决数学模型根据时间步长的迭代计算出不同时刻下的溃决流量，选择出最大值，即峰值流量。该方法计算工作量大、硬件设备要求较高，难以在除特定的灾害监测实验位点以外的灾害现场及时实施应用。除此技术以外，现有的堰塞坝溃决洪水峰值流量计算方法多是通过历史滑坡堰塞坝溃决洪水资料建立，通常堰塞坝高度、堰塞坝溃口深度、堰塞坝库容和堰塞坝溃决水量共4个参数中的1个或2个参数，大多未考虑坝体的易侵蚀性对溃决洪水峰值流量的影响。此外，现有滑坡堰塞坝溃决洪水峰值流量计算公式多存在量纲不和谐问题。

发明内容

本发明的目的就是针对现有技术的不足，提供一种滑坡堰塞坝溃决洪水峰值流量测算方法。

为实现上述目的，本发明首先提供一种滑坡堰塞坝坝体易侵蚀性判断方法，其技术方案如下：

一种滑坡堰塞坝坝体水流易侵蚀性判断方法，用于判断水流对滑坡堰塞坝坝体侵蚀的容易程度，其特征在于：以坝体物质中值粒径d₅₀为指标，若d₅₀<0.50mm，判定为高度易侵蚀性，若0.50mm≤d₅₀≤20mm，判定为中度易侵蚀性，若d₅₀>20mm，判定为低度易侵蚀性。

上述滑坡堰塞坝坝体水流易侵蚀性判断方法中，水流易侵蚀性是指水流对滑坡堰塞坝坝体产生侵蚀作用的容易程度。在滑坡堰塞坝溃决发生与演变过程中，坝体的易侵蚀性高低直接影响溃口形态的变化，进而影响到溃决的严重程度、溃决峰值流量的变化、下游区域受灾程度。

以上述滑坡堰塞坝坝体水流易侵蚀性判断方法为基础，本发明同时提供一种滑坡堰塞坝溃口深度测算方法，其技术方案如下：

一种滑坡堰塞坝溃口深度测算方法，其特征在于：首先，现场调查，获取基本调查数据、判断坝体水流易侵蚀性，所述基本调查数据包括堰塞坝高度h_d、堰塞坝库容V_l；其次，根据坝体水流易侵蚀性依表1确定4个计算系数λ、ξ、α、β的取值，

表1

坝体水流易侵蚀性	λ	ξ	α	β
					高度	0.27	4.554	0.283	0.433
中度	0.12	5.641	0.047	0.532
					低度	0.07	0.7247	0.404	0.612

最后，依式1计算滑坡堰塞坝溃口深度h_b

式中，h_b—滑坡堰塞坝溃口深度，单位m，

H_r—单位长度，取值1.00m，

h_d—堰塞坝高度，单位m，基本调查数据确定，

V_l—堰塞坝库容，单位m³，基本调查数据确定，

λ、ξ、α、β—计算系数，依表1确定。

以上述滑坡堰塞坝坝体水流易侵蚀性判断方法为基础，本发明同时提供一种滑坡堰塞坝溃口底宽度测算方法，其技术方案如下：

一种滑坡堰塞坝溃口底宽度测算方法，其特征在于：首先，现场调查，获取基本调查数据、判断坝体水流易侵蚀性，所述基本调查数据包括堰塞坝库容V_l；其次，根据坝体易侵蚀性依表2确定计算系数λ的取值，

表2

坝体水流易侵蚀性	λ
		高度	0.27
中度	0.12
		低度	0.07

最后，依式2计算滑坡堰塞坝溃口底宽度W_b

式中，W_b—滑坡堰塞坝溃口底宽度，单位m，

V_l—堰塞坝库容，单位m³，基本调查数据确定，

λ—计算系数，依表2确定。

本发明上述滑坡堰塞坝溃口深度测算方法与滑坡堰塞坝溃口底宽度测算方法，是在滑坡灾害发生并形成滑坡堰塞坝后，快速测算可能的溃决事件中溃口形态特征指标的方法。测算方法是在滑坡堰塞坝溃决过程中在考虑了坝体易受水流侵蚀性的前提下对溃口形态特征变化测算方法。借助滑坡堰塞坝溃口底宽度W_b与滑坡堰塞坝溃口深度h_b两项变量，能够很好地将溃口断面形状与溃口水深两项指标特征化，从而为计算溃决中的峰值流量Q提供了简洁而有效的中间量。据此，本发明同时提供一种滑坡堰塞坝溃决洪水峰值流量测算方法，用于测算可能的溃决事件中溃决洪水锋值流量规模，其技术方案如下：

利用上述滑坡堰塞坝溃口深度测算方法、上述滑坡堰塞坝溃口底宽度测算方法实现的滑坡堰塞坝溃决洪水峰值流量测算方法，其特征在于：

首先，现场调查，获取基本调查参数、判断堰塞坝坝体易侵蚀性，所述基本调查参数包括堰塞坝高度h_d、堰塞坝库容V_l、溃口左岸坡比M_l、溃口右岸坡比M_r；

其次，测算滑坡堰塞坝溃口深度h_b；

再次，测算滑坡堰塞坝溃口底宽度W_b；

再次，依式3计算确定溃口断面形状结合参数q

式中，q—溃口断面形状结合参数，无量纲参数，

W_b—滑坡堰塞坝溃口底宽，单位m，

M_l、M_r—溃口左右岸坡比，基本调查数据确定；

再次，依式4、式5计算确定溃口水深特征参数W

式中，W—溃口水深特征参数W，无量纲参数，

W_u—溃口初始水深特征参数，无量纲参数，

h_b—滑坡堰塞坝溃口深度，单位m；

最后，依式6计算确定滑坡堰塞坝溃决洪水峰值流量Q

式中，Q—滑坡堰塞坝溃决洪水峰值流量，单位m³/s，

g—重力加速度，单位s/m²。

与现有技术相比，本发明的有益效果是：本发明提供了一种快速判断滑坡堰塞坝坝体水流易侵蚀性的方法。该方法只需一个现场采集指标，便可以将滑坡堰塞坝坝体易受水流侵蚀性划分为高度、中度、低度三个水平，能够为测算溃决过程中与溃口形态变化有关的指标提供基础。本发明同时提供一种滑坡堰塞坝溃口深度测算方法与一种滑坡堰塞坝溃口底宽度测算方法，此两项方法均借助坝体水流易侵蚀性指标将特定计算系数引入到测算方法中，从而在测算过程中充分考虑了堰塞坝的水流易侵蚀性对溃口形成与演化的影响。利用溃口底宽度W_b与溃口深度h_b能够很好地将溃口断面形状与溃口水深两项指标特征化的技术有效性，本发明最终提供一种滑坡堰塞坝溃决洪水峰值流量测算方法。该方法充分考虑滑坡堰塞坝坝体的易侵蚀性、库容和坝高对滑坡堰塞坝溃决峰值流量的影响，采用半解析公式计算滑坡堰塞坝溃决峰值流量。与现有技术测算结果相比，采用本方法测量的峰值流量更接近灾害现场测量值。整体上，本发明提供了一系列建立在快速划分滑坡堰塞坝坝体水流易侵蚀性水平基础上的涉及滑坡堰塞坝溃决过程重要指标的测算方案。各测算方案均是针对性地解决在发生滑坡堵江形成滑坡堰塞坝后，如何在充分考虑坝体水流侵蚀性特征的前提下对溃决事件中的重要指标进行快速测算的技术问题。通过引入受水流易侵蚀性水平控制的计算系数，各项技术方案将坝体易受水流侵蚀性特征纳入到计算过程，从而对溃口形态与溃决洪水锋值流量的测算更具科学原理。整套测算方法所需基本参数少、计算量小、计算过程简洁、数据采集与计算硬件要求低，能够充分满足滑坡堰塞坝形成与溃决事件中快速响应的防灾减灾措施需要。

附图说明

图1是滑坡堰塞坝溃口形状示意图。

具体实施方式

下面结合附图，对本发明的优选实施例作进一步的描述。

实施例一

2018年10月，长江上游金沙江某处发生1次滑坡，滑坡体堵江形成堰塞湖。该滑坡堰塞坝溃决可能对下游沿岸居民和交通等基础设施造成危害，为制定有效的防灾减灾预案，需要科学合理测算滑坡堰塞坝溃决洪水峰值流量。

首先，现场调查，获得基本调查数据。现场调查包括了针对工程所在山洪泥石流沟道现场的各种测绘、测量、模拟实验测试，以及历史灾害记录获取，以及有参照借鉴作用的经验数据获取等。

基本调查数据包括坝体物质中值粒径d₅₀＝5.2mm，并据此判断坝体属中度水流易侵蚀性；滑坡堰塞坝高度h_d＝61.00m；依据该滑坡堰塞坝的高度h_d，结合滑坡坝上游数字高程地形，确定滑坡堰塞坝库容V_l＝249000000m³。由于溃口左右岸坡比M_l、M_r不能现场确定，故依经验分别取值1.0。图1是滑坡堰塞坝溃口形状示意图。

其次，依表1确定4个计算系数λ＝0.12、ξ＝5.641、α＝0.047、β＝0.532。

再次，将H_r＝1.00m、h_d、V_l、λ、ξ、α、β代入式1计算确定滑坡堰塞坝溃口深度h_b＝23.68m，将V_l、λ代入式2计算确定滑坡堰塞坝溃口底宽度W_b＝75.49m。

再次，将W_b、M_l、M_r代入式3计算确定q＝37.75。

再次，将h_b、q代入式4计算确定溃口初始水深特征参数W_u＝0.5601，将W_u＝0.5601代入式5计算确定溃口水深特征参数W＝0.3920。

最后，将g＝9.81m/s²、q、M_l、M_r、W代入式6计算确定Q＝8304m³/s。

对比测算：采用现有技术(Pierce MW,Thornton CI,Abt SR(2010).Predictionpeak outflow from breached embankment dams.Colorado State University DarylB.Simons Building at the Engineering Research Center Fort Collins,CO)对本滑坡堰塞湖的溃决事件进行测算。计算滑坡堰塞坝溃决洪水峰值流量Q＝16556m³/s

结果验证：该滑坡堰塞坝在形成当月发生溃决，现场观测到的溃决洪水峰值流量约10000m³/s。对比可见本研究计算相对误差小于现有技术。

Claims (5)

1.滑坡堰塞坝坝体水流易侵蚀性判断方法，用于判断水流对滑坡堰塞坝坝体侵蚀的容易程度，其特征在于：以坝体物质中值粒径d₅₀为指标，若d₅₀<0.50mm，判定为高度易侵蚀性，若0.50mm≤d₅₀≤20mm，判定为中度易侵蚀性，若d₅₀>20mm，判定为低度易侵蚀性。

2.滑坡堰塞坝溃口深度测算方法，其特征在于：首先，现场调查，获取基本调查数据、判断坝体水流易侵蚀性，所述基本调查数据包括堰塞坝高度h_d、堰塞坝库容V_l；其次，根据坝体水流易侵蚀性依表1确定4个计算系数λ、ξ、α、β的取值，

表1

坝体水流易侵蚀性 λ ξ α β 高度 0.27 4.554 0.283 0.433 中度 0.12 5.641 0.047 0.532 低度 0.07 0.7247 0.404 0.612

最后，依式1计算滑坡堰塞坝溃口深度h_b

式中，h_b—滑坡堰塞坝溃口深度，单位m，

H_r—单位长度，取值1.00m，

h_d—堰塞坝高度，单位m，基本调查数据确定，

V_l—堰塞坝库容，单位m³，基本调查数据确定，

λ、ξ、α、β—计算系数，依表1确定。

3.滑坡堰塞坝溃口底宽度测算方法，其特征在于：首先，现场调查，获取基本调查数据、判断坝体水流易侵蚀性，所述基本调查数据包括堰塞坝库容V_l；其次，根据坝体易侵蚀性依表2确定计算系数λ的取值，

表2

坝体水流易侵蚀性 λ 高度 0.27 中度 0.12 低度 0.07

最后，依式2计算滑坡堰塞坝溃口底宽度W_b

W_b＝λV_l ^1/3 式2

式中，W_b—滑坡堰塞坝溃口底宽度，单位m，

V_l—堰塞坝库容，单位m³，基本调查数据确定，

λ—计算系数，依表2确定。

4.利用权利要求2所述的滑坡堰塞坝溃口深度测算方法、权利要求3所述的滑坡堰塞坝溃口底宽度测算方法实现的滑坡堰塞坝溃决洪水峰值流量测算方法，其特征在于：

首先，现场调查，获取基本调查参数、判断坝体易侵蚀性，所述基本调查参数包括堰塞坝高度h_d、堰塞坝库容V_l、溃口左岸坡比M_l、溃口右岸坡比M_r；

其次，测算滑坡堰塞坝溃口深度h_b；

再次，测算滑坡堰塞坝溃口底宽度W_b；

再次，依式3计算确定溃口断面形状结合参数q

式中，q—溃口断面形状结合参数，无量纲参数，

W_b—滑坡堰塞坝溃口底宽，单位m，

M_l、M_r—溃口左右岸坡比，基本调查数据确定；

再次，依式4、式5计算确定溃口水深特征参数W

式中，W—溃口水深特征参数W，无量纲参数，

W_u—溃口初始水深特征参数，无量纲参数，

h_b—滑坡堰塞坝溃口深度，单位m；

最后，依式6计算确定滑坡堰塞坝溃决洪水峰值流量Q

式中，Q—滑坡堰塞坝溃决洪水峰值流量，单位m³/s，

g—重力加速度，单位s/m²。

5.根据权利要求4所述的滑坡堰塞坝溃决洪水峰值流量测算，其特征在于：溃口左右岸坡比M_l、M_r分别取值1.0。

Images