一种震区泥石流沟道堆积物起动流深阀值的计算方法技术

本发明专利技术公开了一种震区泥石流沟道堆积物起动流深阀值的计算方法,该方法通过获取震区泥石流沟的沟道坡度α，沟道纵比降J，沟道堆积物的内摩擦角θ，泥石流黏滞系数η，泥石流容重γ，泥石流体中水的容重γW，沟道堆积物饱和容重γsat；将所获得的参数带入计算公式，可以获得泥石流起动过程的预测泥石流流深的阀值。该方法适用于震区泥石流沟道堆积物起动时的流深预测，与现有技术相比，本发明专利技术基于严密的理论推导，能更准确的得到震区泥石流起动流深，为震区泥石流预防提供依据，其运用简单、高效，可满足地质灾害预防需要。

Method for calculating threshold value of incipient flow depth of debris flow channel accumulation in earthquake area

The invention discloses an earthquake debris flow gully debris flow starting deep threshold calculation method, the method by obtaining the earthquake debris flow gully slope alpha, gully slope J, channel deposits of internal friction angle, the viscous coefficient of debris flow, debris flow density, the density of W debris fluid channel deposits, saturation density of sat; into the calculation formula of parameters will be obtained, can predict debris flow process of debris flow depth threshold. This method is applicable to the earthquake debris flow gully debris flow starting depth prediction, compared with the prior art, the invention is based on the strict theoretical derivation, can get more accurate earthquake debris flow starting flow depth, provide the basis for earthquake debris flow prevention, the use of simple, efficient, and can meet the needs of geological disaster prevention.

【技术实现步骤摘要】
一种震区泥石流沟道堆积物起动流深阀值的计算方法
本专利技术属于泥石流预警程
，特别涉及一种震区泥石流沟道堆积物起动流深阀值的计算方法。
技术介绍
“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。其中本文档来自技高网...

【技术保护点】
一种震区泥石流沟道堆积物起动流深阀值的计算方法，其特征在于：将震区泥石流沟道特征、沟道堆积物性质应用于震区泥石流起动的流深预测，其具体计算方法如下：A.通过获取震区泥石流沟道的参数，确定泥石流沟的沟道纵比降J(‰)，沟道坡度α(°)，泥石流沟道堆积物饱和状态时的饱和容重γsat(KN/m

【技术特征摘要】
1.一种震区泥石流沟道堆积物起动流深阀值的计算方法，其特征在于：将震区泥石流沟道特征、沟道堆积物性质应用于震区泥石流起动的流深预测，其具体计算方法如下：A.通过获取震区泥石流沟道的参数，确定泥石流沟的沟道纵比降J(‰)，沟道坡度α(°)，泥石流沟道堆积物饱和状态时的饱和容重γsat(KN/m3)，泥石流容重γ(KN/m3),泥石流体中水的容重γW(KN/m3)，γW＝10KN/m3，泥石流堆积物的内摩擦角θ(°)，泥石流的黏滞系数η(Pa·s)；B.通过以下公式确定震区泥石流起动的流深阀值H，单位m；H＝[(γcosαta...

【专利技术属性】
技术研发人员：屈永平，
申请(专利权)人：四川建筑职业技术学院，
类型：发明
国别省市：四川,51