CN107180150A - 一种震区泥石流沟道堆积物起动流深阀值的计算方法 - Google Patents
一种震区泥石流沟道堆积物起动流深阀值的计算方法 Download PDFInfo
- Publication number
- CN107180150A CN107180150A CN201710588586.1A CN201710588586A CN107180150A CN 107180150 A CN107180150 A CN 107180150A CN 201710588586 A CN201710588586 A CN 201710588586A CN 107180150 A CN107180150 A CN 107180150A
- Authority
- CN
- China
- Prior art keywords
- flow
- mud
- raceway groove
- rock
- rock flow
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Abstract
本发明公开了一种震区泥石流沟道堆积物起动流深阀值的计算方法,该方法通过获取震区泥石流沟的沟道坡度α,沟道纵比降J,沟道堆积物的内摩擦角θ,泥石流黏滞系数η,泥石流容重γ,泥石流体中水的容重γW,沟道堆积物饱和容重γsat;将所获得的参数带入计算公式,可以获得泥石流起动过程的预测泥石流流深的阀值。该方法适用于震区泥石流沟道堆积物起动时的流深预测,与现有技术相比,本发明基于严密的理论推导,能更准确的得到震区泥石流起动流深,为震区泥石流预防提供依据,其运用简单、高效,可满足地质灾害预防需要。
Description
技术领域
本发明属于泥石流预警程技术领域,特别涉及一种震区泥石流沟道堆积物起动流深阀值的计算方法。
背景技术
“5.12”地震以后,强震区内诱发了大量的地震次生灾害,其中强震区内每年都有不同规模的泥石流事件暴发,给强震区带来了严重的人员伤害和巨大的经济损失,且强震区地质灾害的活跃期将持续于地震后的20年。泥石流的预警是泥石流防治的主要方法之一,现有的泥石流预警主要基于泥石流流域的降雨条件和山洪流深,且主要的方法是基于统计方法和室内实验。根据汶川震区泥石流沟的沟道纵比降大、泥石流流域面积小等特征,使得震区泥石流的汇水速度快,泥石流冲击力强。为此,震区泥石流的预警和防治是泥石流研究的重点。
一定流深的地表径流形成一定的动水压力,对泥石流沟道堆积物形成冲击力,当泥石流沟道堆积物的抗剪强度小于等于地表径流的冲击力时,泥石流沟道堆积物处于临界起动状态,随着地表径流流深的增加,冲蚀泥石流沟道堆积物进而诱发泥石流过程。泥石流流深是泥石流暴发的最直接影响因素,也是泥石流预警预报最主要的参数之一。地表径流流深与其所在的沟道冲刷断面、沟道纵比降,以及沟道糙率系数等因素相关。由于汶川震区泥石流沟道纵比降大、泥石流流域面积小,使得在强降雨汇流条件下泥石流流域内的汇水速度快,形成地表径流冲蚀泥石流物源进而诱发降雨型泥石流。例如,小岗剑泥石流沟灾害点在2009年的“9.8”、2011年的“9.5”以及2016年的“1.21”均发生地质灾害事件,严重堵断公路,多次堵塞绵远河。
然而,目前基于泥石流流深的预警预报主要是根据不同地区的经验值和统计模型,目前泥石流预警预报主要还是基于降雨强度对泥石流进行。由于降雨强度的分布的不均一性,使得泥石流预警存在一定误差,且相对于雨量资料较少的地区,根据降雨强度的泥石流预警存在风险,不能高效、便捷的推广到设计单位和生产单位。
发明内容
本发明的目的在于:针对上述存在的问题,提供一种震区泥石流沟道堆积物起动流深阀值的计算方法,解决强震区泥石流沟道堆积物起动过程中临界泥石流流深,构建基于泥石流沟道坡度、沟道堆积物性质的泥石流流深阀值的预测模型,并进行实例运用,为泥石流预警的提供新方法,适用于泥石流预警实际工程的需要。
本发明技术的技术方案实现方式:一种震区泥石流沟道堆积物起动流深阀值的计算方法,其特征在于:将震区泥石流沟道特征、沟道堆积物性质应用于震区泥石流起动的流深预测,其具体计算方法如下:
A.通过获取震区泥石流沟道的参数,确定泥石流沟的沟道纵比降J(‰),沟道坡度α(°),泥石流沟道堆积物饱和状态时的饱和容重γsat(KN/m3),泥石流容重γ(KN/m3),泥石流体中水的容重γW(KN/m3),γW=10KN/m3,泥石流堆积物的内摩擦角θ(°),泥石流的黏滞系数η(Pa·s);
B.通过以下公式确定震区泥石流起动的流深阀值H,单位m;
H=[(γcosαtanθ-(γ-γW)sinα)/(ηJ1/2)]-3/4
式中:J为沟道纵比降(‰),α为沟道坡度(°),γ为泥石流容重(KN/m3),γW为泥石流体中水的容重(KN/m3),γW=10KN/m3,θ为泥石流堆积物的内摩擦角(°),η为泥石流的黏滞系数(Pa·s);均由步骤A所确定。
本发明所述的震区泥石流沟道堆积物起动流深阀值的计算方法,其所述步骤B中的确定震区泥石流起动的流深阀值公式适用于震区泥石流沟道堆积物在地表径流冲蚀条件下的泥石流起动的流深阀值的预测,将得到的泥石流流深将作为震区泥石流暴发的最小流深。
与现有技术相比,本发明的有益效果是:本发明突破传统经验的限制,从泥石流沟道纵比降、泥石流沟道宽度等多方面入手,构建了基于泥石流沟道参数的泥石流流深阀值计算模型。首先,在泥石流沟道坡度方面,利用泥石流沟道纵比降变化特征定量反映了震区泥石流沟道堆积物起动时的流深特征;其次,根据泥石流沟道堆积物的内摩擦角反映了汶川震区泥石流流深与沟道堆积物物理性质间的相互关系,计算理论完善,计算结果符合泥石流现场调查,能够为震区泥石流预警参数提供技术支持。
附图说明
图1是泥石流沟道堆积物受力分析示意图。
图中标号如下:△h为沟道堆积物任意侵蚀厚度,△l为堆积物任意侵蚀长度,h为泥石流堆积厚度,α为泥石流沟道坡度,γsat为泥石流沟道堆积物饱和状态时的饱和容重,dN为堆积物所受支撑应力,dGY为泥石流固体颗粒的自重应力,dτ为堆积物所受剪应力。
具体实施方式
下面结合附图,对本发明作详细的说明。
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明技术进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定发明。
如图1所示,一种震区泥石流沟道堆积物起动流深阀值的计算方法,其主要思路是:将震区泥石流沟道特征、沟道堆积物性质应用于震区泥石流起动的流深预测。首选根据土力学推导泥石流沟道堆积物起动时的受力条件;再通过震区典型泥石流沟的调查以及泥石流沟道特征参数的测量等手段确定沟道坡度、沟道纵比降、沟道糙率系数、泥石流容重、泥石流粘滞系数;将所得参数带入震区泥石流沟道堆积物起动流深阀值的计算模型中,得到泥石流流起动的临界流深预测值。
具体计算方法及步骤如下:
首先,根据震区泥石流沟道堆积物的受力特征绘制泥石流沟道堆积物的受力示意图,分析泥石流沟道堆积物起动过程中的驱动力条件和阻力条件。假设震区泥石流沟道坡度为α,泥石流沟道纵比降为J,泥石流容重为γ,泥石流沟道的糙率系数为n,泥石流的粘滞系数为η,如图1所示,泥石流沟道堆积物起动时,堆积物受地表径流的受剪应力(τ2)公式表示为:
式中,τ2为泥石流剪应力(Pa),τB为泥石流浆体的屈服剪应力(Pa),η为泥石流的黏滞系数(Pa·s),在常温20℃时,粘滞系数为1.005Pa·s,v为泥石流的流速(m/s),α为泥石流沟道坡度(°),γ为泥石流容重(KN/m3),γW为泥石流体中水的容重(KN/m3),γW=10KN/m3。
则泥石流沟道堆积物受到的剪应力(τ2)公式表示为:
然后,泥石流沟道物源堆积体在上覆地表径流自重压力和沟道堆积物的自重应力下的正应力(δ)和剪应力(τ1)公式表示为:
式中,H为地表径流的流深(m),△h为沟道堆积物的任意侵蚀厚度(m),γsat为泥石流沟道堆积物饱和状态时的饱和容重(KN/m3),其他符号同上所示。
然后,根据土力学理论得到泥石流沟道堆积物的抗剪强度τf公式表示为:
式中,c为泥石流沟道堆积物的内聚力(Pa),其他符合同上。
由于汶川震区泥石流沟道堆积物沉积时间短,堆积体粘粒含量低,即泥石流沟道堆积物的内聚力c忽略不计,当泥石流沟道堆积物的抗剪强度τf与泥石流的剪应力τ相等时,泥石流沟道堆积物处于极限稳定状态。则泥石流沟道堆积物侵蚀深度△h=0的泥石流起动流深阀值公式表示为:
H=[(γcosαtanθ-(γ-γW)sinα)/(ηJ1/2)]-3/4…………………(5)
式中:J为沟道纵比降(‰),α为沟道坡度(°),θ为泥石流堆积物的内摩擦角(°),η为泥石流的黏滞系数(Pa·s),在常温20℃时,粘滞系数为1.005Pa·s,γ为泥石流容重(KN/m3),γW为泥石流体中水的容重(KN/m3),γW=10KN/m3。
其中,所述步骤B中的确定震区泥石流起动的流深阀值的公式适用于震区泥石流沟道堆积物在地表径流冲蚀条件下泥石流起动的流深阀值的预测,将得到的泥石流流深将作为震区泥石流暴发的最小流深。
以下为本发明的具体实施例:
红椿沟流域形态为扇形,流域面积5.35km2,流域最高点位于新店子沟的沟源处(望乡石),高程2168.4m,沟口与岷江交汇处高程为880m,相对高差1288.4m。红椿沟发育有甘溪铺沟、大水沟和新店子沟等3条较大的支沟,其中主沟长约3.6km,上陡下缓,沟谷平均纵比降358‰。
为了有效的防治和治理泥石流,防止岷江被堵塞威胁下游安全,根据泥石流沟道特征对泥石流起动的流深阀值进行有效预测,具体计算方法及步骤如下:
A.根据红椿沟泥石流沟的野外调查和泥石流沟道参数测试,确定了红椿沟泥石流沟道堆积物的内摩擦角为33°~35°,沟谷平均纵比降358‰,粘性泥石流沟床糙率为0.08,泥石流重度为17.10KN/m3,粘滞系数为1.005Pa·s。
B.将A步骤确定的各参数代入下公式,
H=[(γcosαtanθ-(γ-γW)sinα)/(ηJ1/2)]-3/4
C.通过公式计算获取红椿沟泥石流起动时的流深阀值约为15cm,则红椿沟泥石流起动的最小流深约为15cm。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种震区泥石流沟道堆积物起动流深阀值的计算方法,其特征在于:将震区泥石流沟道特征、沟道堆积物性质应用于震区泥石流起动的流深预测,其具体计算方法如下:
A.通过获取震区泥石流沟道的参数,确定泥石流沟的沟道纵比降J(‰),沟道坡度α(°),泥石流沟道堆积物饱和状态时的饱和容重γsat(KN/m3),泥石流容重γ(KN/m3),泥石流体中水的容重γW(KN/m3),γW=10KN/m3,泥石流堆积物的内摩擦角θ(°),泥石流的黏滞系数η(Pa·s);
B.通过以下公式确定震区泥石流起动的流深阀值H,单位m;
H=[(γcosαtanθ-(γ-γW)sinα)/(ηJ1/2)]-3/4
式中:J为沟道纵比降(‰),α为沟道坡度(°),泥石流容重γ(KN/m3),泥石流体中水的容重γW(KN/m3),γW=10KN/m3,θ为泥石流堆积物的内摩擦角(°),η为泥石流的黏滞系数(Pa·s);均由步骤A所确定。
2.根据权利要求1所述的震区泥石流沟道堆积物起动流深阀值的计算方法,其特征在于:所述步骤B中的确定震区泥石流起动的流深阀值的公式适用于震区泥石流沟道堆积物在地表径流冲蚀条件下的泥石流起动的流深阀值的预测,将得到的泥石流流深将作为震区泥石流暴发的最小流深。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710588586.1A CN107180150B (zh) | 2017-07-19 | 2017-07-19 | 一种震区泥石流沟道堆积物起动流深阀值的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710588586.1A CN107180150B (zh) | 2017-07-19 | 2017-07-19 | 一种震区泥石流沟道堆积物起动流深阀值的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107180150A true CN107180150A (zh) | 2017-09-19 |
CN107180150B CN107180150B (zh) | 2020-05-19 |
Family
ID=59838432
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710588586.1A Active CN107180150B (zh) | 2017-07-19 | 2017-07-19 | 一种震区泥石流沟道堆积物起动流深阀值的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107180150B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108021764A (zh) * | 2017-12-15 | 2018-05-11 | 四川建筑职业技术学院 | 一种冰川泥石流起动流深的计算方法 |
CN108335034A (zh) * | 2018-01-31 | 2018-07-27 | 中国科学院、水利部成都山地灾害与环境研究所 | 格子坝拦挡粘性泥石流闭塞度评价方法 |
CN108446525A (zh) * | 2018-05-16 | 2018-08-24 | 四川大学 | 山区流域山洪泥沙灾害易灾区识别方法 |
CN110147609A (zh) * | 2019-05-20 | 2019-08-20 | 成都理工大学 | 一种基于流深的泥石流沟道物源侵蚀深度的计算方法 |
CN110147394A (zh) * | 2019-05-25 | 2019-08-20 | 成都理工大学 | 一种基于统计学方法的泥石流百年淤积深度预测方法 |
CN110162866A (zh) * | 2019-05-16 | 2019-08-23 | 成都理工大学 | 一种基于钻孔数据的泥石流沟道淤积深度计算方法 |
CN110532683A (zh) * | 2019-08-28 | 2019-12-03 | 四川建筑职业技术学院 | 一种基于沟道流深、压应力的泥石流容重计算方法 |
CN115290578A (zh) * | 2022-08-15 | 2022-11-04 | 中国科学院空天信息创新研究院 | 一种泥石流沟细粒堆积物粘聚力高光谱探测方法和系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6478954B1 (en) * | 2001-06-06 | 2002-11-12 | Fresh Creek Technologies, Inc. | Debris collecting apparatus |
US6519856B1 (en) * | 1998-06-22 | 2003-02-18 | Delphi Oracle Corp | Safety razor head with intrinsic fencing and lateral skin tensioning |
CN103453944A (zh) * | 2013-09-06 | 2013-12-18 | 中国科学院、水利部成都山地灾害与环境研究所 | 泥石流起动临界清水流量测算方法及起动监测方法 |
CN103556602A (zh) * | 2013-11-21 | 2014-02-05 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种防治沟道泥石流起动的结构体及其设计方法 |
CN104809345A (zh) * | 2015-04-24 | 2015-07-29 | 成都理工大学 | 一种泥石流屈服应力和最大淤积厚度计算方法 |
CN106223259A (zh) * | 2016-09-28 | 2016-12-14 | 中国科学院、水利部成都山地灾害与环境研究所 | 沟道泥石流模拟实验系统及确定沟道泥石流起动流量阈值的方法 |
-
2017
- 2017-07-19 CN CN201710588586.1A patent/CN107180150B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6519856B1 (en) * | 1998-06-22 | 2003-02-18 | Delphi Oracle Corp | Safety razor head with intrinsic fencing and lateral skin tensioning |
US6478954B1 (en) * | 2001-06-06 | 2002-11-12 | Fresh Creek Technologies, Inc. | Debris collecting apparatus |
CN103453944A (zh) * | 2013-09-06 | 2013-12-18 | 中国科学院、水利部成都山地灾害与环境研究所 | 泥石流起动临界清水流量测算方法及起动监测方法 |
CN103556602A (zh) * | 2013-11-21 | 2014-02-05 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种防治沟道泥石流起动的结构体及其设计方法 |
CN104809345A (zh) * | 2015-04-24 | 2015-07-29 | 成都理工大学 | 一种泥石流屈服应力和最大淤积厚度计算方法 |
CN106223259A (zh) * | 2016-09-28 | 2016-12-14 | 中国科学院、水利部成都山地灾害与环境研究所 | 沟道泥石流模拟实验系统及确定沟道泥石流起动流量阈值的方法 |
Non-Patent Citations (2)
Title |
---|
屈永平等: "石棉县熊家沟"7·04"泥石流堵江调查与分析", 《水利学报》 * |
屈永平等: "西藏林芝地区冰川降雨型泥石流起动实验初步研究", 《岩石力学与工程学报》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108021764B (zh) * | 2017-12-15 | 2021-05-07 | 四川建筑职业技术学院 | 一种冰川泥石流起动流深的计算方法 |
CN108021764A (zh) * | 2017-12-15 | 2018-05-11 | 四川建筑职业技术学院 | 一种冰川泥石流起动流深的计算方法 |
CN108335034A (zh) * | 2018-01-31 | 2018-07-27 | 中国科学院、水利部成都山地灾害与环境研究所 | 格子坝拦挡粘性泥石流闭塞度评价方法 |
CN108335034B (zh) * | 2018-01-31 | 2022-11-18 | 中国科学院、水利部成都山地灾害与环境研究所 | 格子坝拦挡粘性泥石流闭塞度评价方法 |
CN108446525A (zh) * | 2018-05-16 | 2018-08-24 | 四川大学 | 山区流域山洪泥沙灾害易灾区识别方法 |
CN108446525B (zh) * | 2018-05-16 | 2020-01-21 | 四川大学 | 山区流域山洪泥沙灾害易灾区识别方法 |
CN110162866A (zh) * | 2019-05-16 | 2019-08-23 | 成都理工大学 | 一种基于钻孔数据的泥石流沟道淤积深度计算方法 |
CN110162866B (zh) * | 2019-05-16 | 2022-04-22 | 成都理工大学 | 一种基于钻孔数据的泥石流沟道淤积深度计算方法 |
CN110147609A (zh) * | 2019-05-20 | 2019-08-20 | 成都理工大学 | 一种基于流深的泥石流沟道物源侵蚀深度的计算方法 |
CN110147609B (zh) * | 2019-05-20 | 2022-11-15 | 成都理工大学 | 一种基于流深的泥石流沟道物源侵蚀深度的计算方法 |
CN110147394B (zh) * | 2019-05-25 | 2022-06-21 | 成都理工大学 | 一种基于统计学方法的泥石流百年淤积深度预测方法 |
CN110147394A (zh) * | 2019-05-25 | 2019-08-20 | 成都理工大学 | 一种基于统计学方法的泥石流百年淤积深度预测方法 |
CN110532683A (zh) * | 2019-08-28 | 2019-12-03 | 四川建筑职业技术学院 | 一种基于沟道流深、压应力的泥石流容重计算方法 |
CN110532683B (zh) * | 2019-08-28 | 2023-03-10 | 四川建筑职业技术学院 | 一种基于沟道流深、压应力的泥石流容重计算方法 |
CN115290578A (zh) * | 2022-08-15 | 2022-11-04 | 中国科学院空天信息创新研究院 | 一种泥石流沟细粒堆积物粘聚力高光谱探测方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107180150B (zh) | 2020-05-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107180150A (zh) | 一种震区泥石流沟道堆积物起动流深阀值的计算方法 | |
CN105369768B (zh) | 一种粘性泥石流最大冲起高度的测算方法及其应用 | |
CN103453944B (zh) | 泥石流起动临界清水流量测算方法及起动监测方法 | |
CN107169615B (zh) | 一种拦砂坝淤满条件下泥石流流速衰减值的计算方法 | |
CN107133434A (zh) | 一种泥石流临界体积浓度的计算方法 | |
CN104631378B (zh) | 一种粘性泥石流矩形排导槽深度测算方法及应用 | |
CN107169252A (zh) | 一种基于泥石流固体颗粒粒径的震区泥石流容重计算方法 | |
CN104848825B (zh) | 泥石流排导槽肋槛后沟床最大冲刷深度的测算方法及应用 | |
CN107194134A (zh) | 一种基于泥石流固体颗粒粒径的震区泥石流容重计算方法 | |
CN107391924A (zh) | 一种基于泥石流物源粒径参数的震区泥石流容重计算方法 | |
Yangiev et al. | Results and analysis of field research in flood reservoirs in Kashkadarya region | |
Chen et al. | Simulation of runoff and suspended sediment transport rate in a basin with multiple watersheds | |
Marangu et al. | Modeling open channel fluid flow with trapezoidal cross section and a segment base | |
CN107357998A (zh) | 一种泥石流堆积区流量衰减值的计算方法 | |
Kashiwai et al. | Hydraulic examination of Koshibu dam’s intake facilities for sediment bypass | |
CN107180287A (zh) | 一种基于汶川强震区的泥石流物源转换为泥石流的转换率计算方法 | |
Wang et al. | Total sediment transport from an urbanizing watershed in the upper Yellow River, China | |
Abbs et al. | A model study of the hydraulics related to fish passage through backwatered culverts | |
CN110532683A (zh) | 一种基于沟道流深、压应力的泥石流容重计算方法 | |
CN107330565A (zh) | 一种饱水状态下沟道堆积体破坏角的预测方法 | |
Ikeda et al. | Study of prediction methods of debris-flow peak discharge | |
Jiang et al. | Combined impacts of sand mining activities: The Nandu River downstream segment | |
CN110147394A (zh) | 一种基于统计学方法的泥石流百年淤积深度预测方法 | |
Bong et al. | Self-cleansing design of rectangular open storm sewer | |
Zhang et al. | Analysis of a debris flow after Wenchuan Earthquake and discussion on preventive measures |
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 |