CN108090670A - 一种基于水流侵蚀能力的堰塞坝形成评估方法 - Google Patents
一种基于水流侵蚀能力的堰塞坝形成评估方法 Download PDFInfo
- Publication number
- CN108090670A CN108090670A CN201711345347.XA CN201711345347A CN108090670A CN 108090670 A CN108090670 A CN 108090670A CN 201711345347 A CN201711345347 A CN 201711345347A CN 108090670 A CN108090670 A CN 108090670A
- Authority
- CN
- China
- Prior art keywords
- landslide
- potential
- potential landslide
- unit
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 53
- 230000003628 erosive effect Effects 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 26
- 239000002689 soil Substances 0.000 claims abstract description 78
- 238000009825 accumulation Methods 0.000 claims abstract description 17
- 238000011156 evaluation Methods 0.000 claims description 17
- 230000015572 biosynthetic process Effects 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000006467 substitution reaction Methods 0.000 claims description 9
- 238000011835 investigation Methods 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000002474 experimental method Methods 0.000 claims description 4
- 238000009533 lab test Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 238000011157 data evaluation Methods 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 238000003921 particle size analysis Methods 0.000 claims description 2
- 238000001556 precipitation Methods 0.000 claims description 2
- 239000011435 rock Substances 0.000 description 8
- 238000012360 testing method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 235000021384 green leafy vegetables Nutrition 0.000 description 2
- 241000251468 Actinopterygii Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000034994 death Effects 0.000 description 1
- 231100000517 death Toxicity 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000004172 quinoline yellow Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Development Economics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
本发明公开了一种基于水流侵蚀能力的堰塞坝形成评估方法。所述评估方法首先通过遥感及地形数据，确定潜在滑坡位置的面积、坡度、高度、坡面长度、宽度以及河道宽度和坡度，然后分别计算潜在滑坡土体堆积速率和河道水流侵蚀速率，最后通过无量纲速率指标判定潜在滑坡土体是否会堵塞主流河道形成堰塞坝。与现有技术相比，本发明方法综合考虑滑坡土体参数和水流参数，提出无量纲速率指标，能够合理评估潜在滑坡土体是否会堵塞主流河道并形成堰塞坝，实现对堰塞坝形成可能性的有效预测，且本发明方法计算参数易于获得，评估方法简便，可为潜在滑坡土体是否形成堰塞坝的评估提供科学依据。
Description
技术领域
本发明涉及一种基于水流侵蚀能力的堰塞坝形成评估方法，特别是适用于潜在滑坡或滑坡形成堰塞坝的评估方法。
背景技术
堰塞坝是由火山熔岩流，冰碛物或由地质灾害、地震活动使山体岩石崩塌下来等原因引起的山崩滑坡体等堵塞主流河道，形成类似大坝的挡水体。其中，滑坡形成的堰塞坝同时受到滑坡基本特征与河流水文条件的控制。对于山坡坡面坡度较陡且河道宽度较小的山谷，滑坡发生后，土体的运动量足够移动至河道或抵达对岸，当河道水流无法即时将正在堆积于河道中的滑坡土体冲刷至下游，且水流深度未高过滑坡土体最小堆积高度时，容易阻塞河道形成堰塞坝。
滑坡土体进入到河道是否形成堰塞坝主要可分为三种情况：1.滑坡土体进入到河道，土体受河道水流挟带至下游，土体无法堆积于河道；2.滑坡土体进入到河道，部份土体受河道水流挟带至下游，部份土体堆积在河道，逐渐形成堰塞坝，但是受到河道水流持续冲刷，在坝体尚未完全形成时，土体即被河道水流挟带至下游；3.滑坡土体进入到河道，河道水流无法即时将正在堆积于河道中的滑坡土体冲刷至下游，滑坡土体逐渐形成堰塞坝。通过滑坡土体进入到河道形成堰塞坝的状况可以知道，如果水流无法即时将滑坡土体挟带至下游，滑坡土体才能不受水流作用影响，顺利堆积在河道并且形成堰塞坝，因此，水流对于土体的侵蚀能力是影响滑坡土体进入河道后是否形成堰塞坝的关键因素，但目前尚没有探讨水流对于土体的侵蚀能力与堰塞坝形成之间的定量关系，因此也就无法有效地评估滑坡土体是否会形成堰塞坝。
发明内容
本发明的目的就是针对现有技术的不足，提供一种基于水流侵蚀能力计算的堰塞坝形成评估方法，该方法分别获取河道水流侵蚀速率及潜在滑坡土体堆积速率，通过模型实验将河道水流侵蚀速率及潜在滑坡土体堆积速率的定量关系无量纲化，能够合理地评估潜在滑坡土体能否形成堰塞坝，实现对堰塞坝形成可能性的有效预测。
为实现上述目的，本发明的技术方案是：
本发明提出一种基于水流侵蚀能力的堰塞坝形成评估方法，针对潜在滑坡，所述评估方法步骤如下：
(一)通过滑坡早期识别图得到潜在滑坡的位置。通过遥感及地形数据，确定潜在滑坡坡度θs、单位°，潜在滑坡高度HL、单位m，潜在滑坡坡面长度SL、单位m，以及潜在滑坡宽度W、单位m；通过潜在滑坡面积与潜在滑坡体积之间的经验关系式、或滑坡体积与前期降水量关系法、或地震震级与滑坡体积统计关系法、或滑坡前后期遥感数据配合高精度DEM数据评估法，确定潜在滑坡体积VL、单位m3。所述潜在滑坡面积与潜在滑坡体积之间的经验关系式为VL＝0.074AL 1.45；式中，VL为潜在滑坡体积、单位m3，AL为潜在滑坡面积、单位m2、通过遥感及地形数据确定。
(二)通过地形数据或现场测量，确定河道宽度B、单位m，及河道坡度θw、单位°。
(三)通过公式As＝W×B，确定潜在滑坡土体堆积于河道的面积As、单位m2；
式中，W为潜在滑坡宽度、单位m、由步骤(一)确定，B为河道宽度、单位m、由步骤(二)确定。
(四)通过公式确定平均滑坡速度单位m/s；式中，g为重力加速度、取值9.8m/s2，HL为潜在滑坡高度、单位m、由步骤(一)确定，θs为潜在滑坡坡度、单位°、由步骤(一)确定，f为等价摩擦系数，通过统计经验关系式、或理论分析式、或数值计算方法、或室内实验方法确定；将确定的平均滑坡速度代入公式计算得到潜在滑坡总时间t、单位s，式中SL为潜在滑坡坡面长度、单位m、由步骤(一)确定。所述确定等价摩擦系数f的统计经验关系式为logf＝-0.15666logVL+0.62419；式中，VL为潜在滑坡体积、单位m3、由步骤(一)确定。
(五)通过以下公式确定潜在滑坡土体堆积速率vd
式中，vd—潜在滑坡土体堆积速率，单位cm/s；
VL—潜在滑坡体积，单位m3，由步骤(一)确定；
As—潜在滑坡土体堆积于河道的面积，单位m2，由步骤(三)确定；
t—潜在滑坡总时间，单位s，由步骤(四)确定。
(六)通过以下公式确定河道水流侵蚀速率ve
式中，ve—河道水流侵蚀速率，单位cm/s；
B—河道宽度，单位m，由步骤(二)确定；
hw—河道水深，单位m，通过现场调查、或流量测量计算、或曼宁公式确定；
θw—河道坡度，单位°，由步骤(二)确定；
ρs—潜在滑坡土体的密度，单位g/cm3，通过现场调查、或室内实验、或历史文献资料确定；
ρ—水的密度，取值1g/cm3；
d50—潜在滑坡土体平均粒径，单位m，通过现场调查、或筛分法、或图像粒径分析方法确定；
g—重力加速度、取值9.8m/s2。
上述河道水流侵蚀速率计算公式是通过室内实验结果得到的对Nakagawa andTsujimoto(1975)提出公式的修正。
(七)将步骤(五)中得到的潜在滑坡土体堆积速率vd和步骤(六)中得到的河道水流侵蚀速率ve代入公式计算得到无量纲速率指标vde；当vde>54时，判定潜在滑坡土体能堵塞主流河道并形成堰塞坝；当vde<47时，判定潜在滑坡土体不会堵塞主流河道，不会形成堰塞坝；当47≤vde≤54时，判定潜在滑坡土体有可能会堵塞主流河道而形成堰塞坝(即介于47～54为不确定区)。无量纲速率指标评估方法是基于46组模型试验，通过对试验结果进行数据判别分析而得(参见附图1)。
上述基于水流侵蚀能力的堰塞坝形成评估方法也可针对滑坡进行，即上述评估方法中涉及的与潜在滑坡有关参数均可替换为与滑坡有关参数(如：潜在滑坡坡度θs替换为滑坡坡度θs，潜在滑坡高度HL替换为滑坡高度HL，潜在滑坡坡面长度SL替换为坡面长度SL，潜在滑坡宽度W替换为滑坡宽度W，潜在滑坡体积VL替换为滑坡体积VL，潜在滑坡土体堆积于河道的面积As替换为滑坡土体堆积于河道的面积As，潜在滑坡总时间t替换为滑坡总时间t，潜在滑坡土体堆积速率vd替换为滑坡土体堆积速率vd，等等)，从而能够合理地评估滑坡土体能否形成堰塞坝。
与现有技术相比，本发明的有益效果是：通过修正河道水流侵蚀速率计算公式，同时综合考虑了滑坡土体参数(堆积速率)和水流参数(侵蚀速率)，基于室内模型实验结果得到河道水流侵蚀速率及潜在滑坡土体堆积速率之间的定量关系，并将定量关系无量纲化，提出无量纲速率指标，能够合理地评估潜在滑坡土体是否会堵塞主流河道并形成堰塞坝，实现对堰塞坝形成可能性的有效预测。本发明方法计算参数易于获得，评估方法简便，可为潜在滑坡土体是否形成堰塞坝的评估提供科学依据。
附图说明
图1是基于46组模型试验结果进行数据判别分析的无量纲速率指标示意图。
图2是实施例一、二、三的无量纲速率指标示意图。
图中符号如下：
vde无量纲速率指标
vd潜在滑坡土体堆积速率
ve河道水流侵蚀速率
具体实施方式
下面对本发明的优选实施例作进一步的描述。
实施例一
如图2所示。老鹰岩滑坡是“5.12”汶川地震触发山体滑坡形成堰塞坝的其中一处案例，位于四川绵阳市安县高川乡泉水村西北侧，地理坐标N31°37′20.4″，E104°08′25.92″，老鹰岩滑坡发生后，形成高速碎屑流冲入黄洞子沟并完全堵塞，形成坝高105m、坝顶宽150m的堰塞坝，汶川地震发生后10天，老鹰岩堰塞坝的蓄水量已达2200000m3，滑坡不仅形成巨大的堰塞坝，还将老鹰岩电站整体掩埋，造成20多人死亡与近千万元的经济损失。
采用本发明的基于水流侵蚀能力的堰塞坝形成评估方法，对老鹰岩滑坡形成堰塞坝进行评估验证，具体步骤如下：
第一步，首先确定滑坡位置，然后通过遥感及地形数据，确定滑坡面积AL为320000m2，滑坡坡度θs为36.7°，滑坡高度HL为820m，滑坡坡面长度SL为1372m，滑坡宽度W为170m；通过滑坡面积与滑坡体积之间的经验关系式VL＝0.074AL 1.45＝0.074×3200001.45，计算得到滑坡体积VL为7107211m3。
第二步，通过地形数据，确定河道宽度B为90m，河道坡度θw为1.6°。
第三步，通过公式As＝W×B＝170×90＝15300，确定滑坡土体堆积于河道的面积As为15300m2。
第四步，通过统计经验关系式logf＝-0.15666logVL+0.62419＝-0.15666log7107211+0.62419，确定等价摩擦系数f为0.356；通过公式确定平均滑坡速度为45.87m/s；将确定的平均滑坡速度代入公式计算得到滑坡总时间t为30s。
第五步，通过公式确定滑坡土体堆积速率vd为
第六步，通过现场调查，确定河道水深hw为0.2m，滑坡土体的密度ρs为2.6g/cm3，滑坡土体平均粒径d50为0.01m。通过公式
确定河道水流侵蚀速率ve为0.53cm/s。
第七步，将第五步中得到的滑坡土体堆积速率vd和第六步中得到的河道水流侵蚀速率ve代入公式计算得到无量纲速率指标vde为2921.53；vde>54，判定老鹰岩滑坡土体能堵塞黄洞子沟并形成堰塞坝，评估结果参见图2所示。判定结果与实际情况相符。
实施例二
如图2所示。S303线映秀至卧龙公路K24+500滑坡发生于2008年5月13日，为滑坡发生未形成堰塞坝的案例，S303线映秀至卧龙公路K24+500滑坡位于汶川县耿达镇，滑坡位置位于公路左侧，渔子溪河右岸。滑坡发生后，滑坡土体冲入渔子溪河道，由于渔子溪水流湍急，大部分滑坡土体被河水迅速冲蚀搬运，并未堵塞河道形成堰塞坝，但是破坏了过河索桥，并掩埋了滑坡位置对面的映秀至卧龙公路，造成公路中断。
采用本发明的基于水流侵蚀能力的堰塞坝形成评估方法，对S303线映秀至卧龙公路K24+500滑坡未形成堰塞坝进行评估验证，具体步骤如下：
第一步，首先确定滑坡位置，然后通过遥感及地形数据，确定滑坡面积AL为6678m2，滑坡坡度θs为42.1°，滑坡高度HL为113m，滑坡坡面长度SL为169m，滑坡宽度W为47m；通过滑坡面积与滑坡体积之间的经验关系式VL＝0.074AL 1.45＝0.074×66781.45，计算得到滑坡体积VL为26000m3。
第二步，通过地形数据，确定河道宽度B为40m，河道坡度θw为1.4°。
第三步，通过公式As＝W×B＝47×40＝1880，确定滑坡土体堆积于河道的面积As为1880m2。
第四步，通过统计经验关系式logf＝-0.15666logVL+0.62419＝-0.15666log26000+0.62419，确定等价摩擦系数f为0.856；通过公式确定平均滑坡速度为5.42m/s；将确定的平均滑坡速度代入公式计算得到滑坡总时间t为31s。
第五步，通过公式确定滑坡土体堆积速率vd为44.61cm/s。
第六步，通过曼宁公式，确定河道水深hw为1m；通过历史文献资料，确定滑坡土体的密度ρs为2.6g/cm3；通过筛分法，确定滑坡土体平均粒径d50为0.02m。通过公式
确定河道水流侵蚀速率ve为1.024cm/s。
第七步，将第五步中得到的滑坡土体堆积速率vd和第六步中得到的河道水流侵蚀速率ve代入公式计算得到无量纲速率指标vde为43.56；vde<47时，判定S303线映秀至卧龙公路K24+500滑坡土体不会堵塞渔子溪，不会形成堰塞坝，评估结果参见图2所示。判定结果与实际情况相符。
实施例三
如图2所示。1999年9月21日，中国台湾地区发生7.3级地震，地震诱发云林县古坑乡草岭村附近山区边坡发生滑坡，滑坡土体堵塞清水溪而形成草岭堰塞湖，堰塞湖形成后，形成50～100m的坝体高度，堰塞湖并未造成严重二次灾害，在后续多次的台风豪雨事件中，堰塞湖湖区逐渐淤积，目前已被土砂淤满。
采用本发明的基于水流侵蚀能力的堰塞坝形成评估方法，针对草岭潜在滑坡形成堰塞坝进行评估验证，具体步骤如下：
第一步，通过地震前期遥感数据得到草岭潜在滑坡的位置。通过遥感及地形数据，确定潜在滑坡面积AL为4000000m2，潜在滑坡坡度θs为31.4°，潜在滑坡高度HL为1100m，潜在滑坡坡面长度SL为2111m，潜在滑坡宽度W为1900m；通过潜在滑坡面积与潜在滑坡体积之间的经验关系式计算得到潜在滑坡体积VL为276833541m3。
第二步，通过地形数据，确定河道宽度B为30m，河道坡度θw为1.4°。
第三步，通过公式As＝W×B＝1900×30＝57000，确定潜在滑坡土体堆积于河道的面积As为57000m2。
第四步，通过统计经验关系式logf＝-0.15666logVL+0.62419＝-0.15666log276833541+0.62419，确定等价摩擦系数f为0.200；通过公式确定平均滑坡速度为60.21m/s；将确定的平均滑坡速度代入公式计算得到潜在滑坡总时间t为35s。
第五步，通过公式确定潜在滑坡土体堆积速率vd为13876.37cm/s。
第六步，通过历史文献资料，确定河道水深hw为0.56m；通过历史文献资料，确定潜在滑坡土体的密度ρs为2.65g/cm3；通过历史文献资料，确定潜在滑坡土体平均粒径d50为0.07m。通过公式
确定河道水流侵蚀速率ve为0.420cm/s。
第七步，将第五步中得到的潜在滑坡土体堆积速率vd和第六步中得到的河道水流侵蚀速率ve代入公式计算得到无量纲速率指标vde为33038.98；vde>54，判定草岭潜在滑坡土体能堵塞清水溪并形成堰塞坝，评估结果参见图2所示。判定结果与实际情况相符。
Claims (4)
1.一种基于水流侵蚀能力的堰塞坝形成评估方法，其特征在于：所述基于水流侵蚀能力的堰塞坝形成评估方法步骤如下：
(一)通过遥感及地形数据，确定潜在滑坡坡度θs、单位°，潜在滑坡高度HL、单位m，潜在滑坡坡面长度SL、单位m，以及潜在滑坡宽度W、单位m；通过潜在滑坡面积与潜在滑坡体积之间的经验关系式、或滑坡体积与前期降水量关系法、或地震震级与滑坡体积统计关系法、或滑坡前后期遥感数据配合高精度DEM数据评估法，确定潜在滑坡体积VL、单位m3；
(二)通过地形数据或现场测量，确定河道宽度B、单位m，及河道坡度θw、单位°；
(三)通过公式As＝W×B，确定潜在滑坡土体堆积于河道的面积As、单位m2；式中，W为潜在滑坡宽度、单位m、由步骤(一)确定，B为河道宽度、单位m、由步骤(二)确定；
(四)通过公式确定平均滑坡速度单位m/s；式中，g为重力加速度、取值9.8m/s2，HL为潜在滑坡高度、单位m、由步骤(一)确定，θs为潜在滑坡坡度、单位°、由步骤(一)确定，f为等价摩擦系数，通过统计经验关系式、或理论分析式、或数值计算方法、或室内实验方法确定；将确定的平均滑坡速度代入公式计算得到潜在滑坡总时间t、单位s，式中SL为潜在滑坡坡面长度、单位m、由步骤(一)确定；
(五)通过以下公式确定潜在滑坡土体堆积速率vd
<mrow>
<msub>
<mi>v</mi>
<mi>d</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>V</mi>
<mi>L</mi>
</msub>
<mrow>
<msub>
<mi>A</mi>
<mi>s</mi>
</msub>
<mo>&times;</mo>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
式中，vd—潜在滑坡土体堆积速率，单位cm/s；
VL—潜在滑坡体积，单位m3，由步骤(一)确定；
As—潜在滑坡土体堆积于河道的面积，单位m2，由步骤(三)确定；
t—潜在滑坡总时间，单位s，由步骤(四)确定；
(六)通过以下公式确定河道水流侵蚀速率ve
<mrow>
<msub>
<mi>v</mi>
<mi>e</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>0.02</mn>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>Bh</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<mi>B</mi>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>h</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>tan&theta;</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
<mi>&rho;</mi>
</mfrac>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>d</mi>
<mn>50</mn>
</msub>
</mrow>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<mn>0.035</mn>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
<mi>&rho;</mi>
</mfrac>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>d</mi>
<mn>50</mn>
</msub>
</mrow>
<mrow>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>Bh</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<mi>B</mi>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>h</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>tan&theta;</mi>
<mi>w</mi>
</msub>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
<msqrt>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
<mi>&rho;</mi>
</mfrac>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>gd</mi>
<mn>50</mn>
</msub>
</mrow>
</msqrt>
<mo>+</mo>
<mn>0.335</mn>
</mrow>
式中，ve—河道水流侵蚀速率，单位cm/s；
B—河道宽度，单位m，由步骤(二)确定；
hw—河道水深，单位m，通过现场调查、或流量测量计算、或曼宁公式确定；
θw—河道坡度，单位°，由步骤(二)确定；
ρs—潜在滑坡土体的密度，单位g/cm3，通过现场调查、或室内实验、或历史文献资料确定；
ρ—水的密度，取值1g/cm3；
d50—潜在滑坡土体平均粒径，单位m，通过现场调查、或筛分法、或图像粒径分析方法确定；
g—重力加速度、取值9.8m/s2；
(七)将步骤(五)中得到的潜在滑坡土体堆积速率vd和步骤(六)中得到的河道水流侵蚀速率ve代入公式计算得到无量纲速率指标vde；当vde>54时，判定潜在滑坡土体能堵塞主流河道并形成堰塞坝；当vde<47时，判定潜在滑坡土体不会堵塞主流河道，不会形成堰塞坝；当47≤vde≤54时，判定潜在滑坡土体有可能会堵塞主流河道而形成堰塞坝。
2.根据权利要求1所述基于水流侵蚀能力的堰塞坝形成评估方法，其特征在于：步骤(一)中潜在滑坡面积与潜在滑坡体积之间的经验关系式为VL＝0.074AL 1.45；式中，VL为潜在滑坡体积、单位m3，AL为潜在滑坡面积、单位m2、通过遥感及地形数据确定。
3.根据权利要求1所述基于水流侵蚀能力的堰塞坝形成评估方法，其特征在于：步骤(四)中确定等价摩擦系数f的统计经验关系式为logf＝-0.15666logVL+0.62419；式中，VL为潜在滑坡体积、单位m3、由步骤(一)确定。
4.根据权利要求1所述基于水流侵蚀能力的堰塞坝形成评估方法，其特征在于：所述潜在滑坡的位置通过滑坡早期识别图得到。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711345347.XA CN108090670B (zh) | 2017-12-15 | 2017-12-15 | 一种基于水流侵蚀能力的堰塞坝形成评估方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711345347.XA CN108090670B (zh) | 2017-12-15 | 2017-12-15 | 一种基于水流侵蚀能力的堰塞坝形成评估方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108090670A true CN108090670A (zh) | 2018-05-29 |
| CN108090670B CN108090670B (zh) | 2021-08-31 |
Family
ID=62176252
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711345347.XA Expired - Fee Related CN108090670B (zh) | 2017-12-15 | 2017-12-15 | 一种基于水流侵蚀能力的堰塞坝形成评估方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108090670B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109657281A (zh) * | 2018-11-26 | 2019-04-19 | 西安理工大学 | 一种引流槽开挖措施的堰塞坝溃决机制分析方法 |
| CN110688773A (zh) * | 2019-10-14 | 2020-01-14 | 中国电建集团成都勘测设计研究院有限公司 | 一种流域堵塞快速定位的系统及方法 |
| CN111401529A (zh) * | 2020-04-21 | 2020-07-10 | 中国水利水电科学研究院 | 一种基于遥感技术的河流生态流量监管方法 |
| CN111814321A (zh) * | 2020-06-24 | 2020-10-23 | 河海大学 | 一种再现堰塞坝堆积特性的数值模拟方法 |
| CN111984914A (zh) * | 2020-08-17 | 2020-11-24 | 中国科学院、水利部成都山地灾害与环境研究所 | 崩滑型堰塞坝坝体高度测算方法、应用 |
| CN112182993A (zh) * | 2020-09-29 | 2021-01-05 | 中国科学院、水利部成都山地灾害与环境研究所 | 滑坡堰塞坝溃口深度测算方法、滑坡堰塞坝溃决洪水峰值流量测算方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104915495A (zh) * | 2015-06-05 | 2015-09-16 | 中国科学院水利部成都山地灾害与环境研究所 | 一种泥石流灾情评估方法及应用 |
| US20150331143A1 (en) * | 2014-05-14 | 2015-11-19 | Foundation of Soongsil University-lndustry Cooperation | Method of recognizing slope condition, system using the same, and recording medium for performing the same |
| JP5917226B2 (ja) * | 2012-03-28 | 2016-05-11 | 日鐵住金建材株式会社 | 透過型砂防堰堤 |
| CN105679167A (zh) * | 2016-01-29 | 2016-06-15 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种新型泥石流物源溃坝实验水槽装置及其实现方法 |
| CN106290790A (zh) * | 2016-08-17 | 2017-01-04 | 西南交通大学 | 可调多功能滑坡堵江及堰塞坝溃坝机制试验装置及其试验方法 |
| KR101707302B1 (ko) * | 2015-03-20 | 2017-02-27 | 영남대학교 산학협력단 | 생태 및 경관 기능을 갖는 흙댐 리모델링 공법 |
-
2017
- 2017-12-15 CN CN201711345347.XA patent/CN108090670B/zh not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5917226B2 (ja) * | 2012-03-28 | 2016-05-11 | 日鐵住金建材株式会社 | 透過型砂防堰堤 |
| US20150331143A1 (en) * | 2014-05-14 | 2015-11-19 | Foundation of Soongsil University-lndustry Cooperation | Method of recognizing slope condition, system using the same, and recording medium for performing the same |
| KR101707302B1 (ko) * | 2015-03-20 | 2017-02-27 | 영남대학교 산학협력단 | 생태 및 경관 기능을 갖는 흙댐 리모델링 공법 |
| CN104915495A (zh) * | 2015-06-05 | 2015-09-16 | 中国科学院水利部成都山地灾害与环境研究所 | 一种泥石流灾情评估方法及应用 |
| CN105679167A (zh) * | 2016-01-29 | 2016-06-15 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种新型泥石流物源溃坝实验水槽装置及其实现方法 |
| CN106290790A (zh) * | 2016-08-17 | 2017-01-04 | 西南交通大学 | 可调多功能滑坡堵江及堰塞坝溃坝机制试验装置及其试验方法 |
Non-Patent Citations (6)
| Title |
|---|
| JIA-JYUN DONG ET AL.: "Deriving landslide dam geometry from remote sensing images for the rapid assessment of critical parameters related to dam-breach hazards", 《LANDSLIDES》 * |
| KUN-TING CHEN ET AL.: "Rapid geometry analysis for earthquake-induced and rainfall-induced landslide dams in Taiwan", 《JOURNAL OF MOUNTAIN SCIENCE》 * |
| XIANBIN YU ET AL.: "The Influence of an EPS Concrete Buffer Layer Thickness on Debris Dams Impacted by Massive Stones in the Debris Flow", 《SHOCK AND VIBRATION》 * |
| 张玉柱 等: "黄河上游积石峡史前滑坡堰塞湖形成年代与发展演变研究", 《中国科学：地球科学》 * |
| 王道正 等: "不同颗粒级配条件下堰塞坝溃决特征试验研究", 《防灾减灾工程学报》 * |
| 胡桂胜 等: "茂县宗渠堰塞湖漫顶溃决洪水演进规律模型试验研究", 《水土保持研究》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109657281A (zh) * | 2018-11-26 | 2019-04-19 | 西安理工大学 | 一种引流槽开挖措施的堰塞坝溃决机制分析方法 |
| CN110688773A (zh) * | 2019-10-14 | 2020-01-14 | 中国电建集团成都勘测设计研究院有限公司 | 一种流域堵塞快速定位的系统及方法 |
| CN111401529A (zh) * | 2020-04-21 | 2020-07-10 | 中国水利水电科学研究院 | 一种基于遥感技术的河流生态流量监管方法 |
| CN111814321A (zh) * | 2020-06-24 | 2020-10-23 | 河海大学 | 一种再现堰塞坝堆积特性的数值模拟方法 |
| CN111814321B (zh) * | 2020-06-24 | 2023-05-23 | 河海大学 | 一种再现堰塞坝堆积特性的数值模拟方法 |
| CN111984914A (zh) * | 2020-08-17 | 2020-11-24 | 中国科学院、水利部成都山地灾害与环境研究所 | 崩滑型堰塞坝坝体高度测算方法、应用 |
| CN111984914B (zh) * | 2020-08-17 | 2023-07-28 | 中国科学院、水利部成都山地灾害与环境研究所 | 崩滑型堰塞坝坝体高度测算方法、应用 |
| CN112182993A (zh) * | 2020-09-29 | 2021-01-05 | 中国科学院、水利部成都山地灾害与环境研究所 | 滑坡堰塞坝溃口深度测算方法、滑坡堰塞坝溃决洪水峰值流量测算方法 |
| CN112182993B (zh) * | 2020-09-29 | 2023-08-11 | 中国科学院、水利部成都山地灾害与环境研究所 | 滑坡堰塞坝溃决洪水峰值流量测算方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108090670B (zh) | 2021-08-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108090670B (zh) | 一种基于水流侵蚀能力的堰塞坝形成评估方法 | |
| Seenipandi et al. | Seawater intrusion vulnerability in the coastal aquifers of southern India—an appraisal of the GALDIT model, parameters’ sensitivity, and hydrochemical indicators | |
| Dong et al. | The formation and breach of a short-lived landslide dam at Hsiaolin village, Taiwan—part I: post-event reconstruction of dam geometry | |
| Schmidt et al. | Hydrologic conditions and terrestrial laser scanning of post-firedebris flows in the San Gabriel Mountains, CA, USA | |
| CN110955952B (zh) | 一种多尺度泥石流危险性评价方法 | |
| Pareta et al. | Hydromorphogeological study of Karawan watershed using GIS and remote sensing techniques | |
| CN112735098B (zh) | 一种震后冰湖溃决预警方法、装置及存储介质 | |
| CN106157541A (zh) | 一种沟谷泥石流预警方法及其应用 | |
| Liu et al. | Effects of loose deposits on debris flow processes in the Aizi Valley, southwest China | |
| Al-Ansari et al. | Flow of River Tigris and its effect on the bed sediment within Baghdad, Iraq | |
| CN109766513A (zh) | 一种泥石流断面平均流速的测算方法 | |
| Huang et al. | The long-term geologic hazards and consequent risk after the Wenchuan earthquake | |
| CN106968217A (zh) | 基于降雨频率的泥石流淤积危险范围划分方法 | |
| Kandekar et al. | Morphometric analysis of Agadgaon watershed using remote sensing and geographic information system | |
| George et al. | Predicting piping erosion susceptibility by statistical and artificial intelligence approaches-A review | |
| CN106192864A (zh) | 一种沟谷泥石流的危险性划分方法及其应用 | |
| Hultén et al. | Case studies of landslide risk due to climate change in Sweden | |
| Seenipandi et al. | Modeling of coastal vulnerability to sea-level rise and shoreline erosion using modified CVI model | |
| Chilikova-Lubomirova et al. | River Hydraulics during Flood Events: The Balkan Experiences | |
| Mirdarsoltany | Investigation and analysis of debris flow phenomenon | |
| Ochiai et al. | Mechanism of large-scale collapse at Tue Valley in the Shikoku mountainous region, Japan | |
| Anusha et al. | Drainage morphometry of Harangi Watershed, Madikeri District, Karnataka State and its influence on the areas sensitive to soil erosion | |
| Ganguly et al. | Spatio-temporal analysis of flood and identification of flood hazard zone of west Tripura district, Tripura, India using integrated geospatial technique, hill Geographer | |
| Rouhi | Landslide damming: implications for natural hazards and risks | |
| MR et al. | Drainage Morphometry of Harangi Watershed, Madikeri District, Karnataka State and Its Influence on the Areas Sensitive to Soil Erosion. |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210831 |