一种基于泥石流固体颗粒粒径的震区泥石流容重计算方法技术

本发明专利技术公开了一种基于泥石流固体颗粒粒径的震区泥石流容重计算方法，该方法通过汶川震区泥石流沟的野外调查、测量、取样以及泥石流沟特征参数提取等手段，确定泥石流沟的沟道纵比降J，泥石流流深h，沟道坡度α，泥石流固体颗粒容重γs，泥石流的粘滞系数η，泥石流固体颗粒粒径di；将所得到的参数代入计算公式以得到震区泥石流容重预测值。与现有技术相比，本发明专利技术基于严密的理论推导，能够合理准确的反映汶川震区低黏粒含量的泥石流容重，为震区灾后重建中的泥石流防治设计提供理论依据，计算结果精度高，符合实际工程需要。

【技术实现步骤摘要】
一种基于泥石流固体颗粒粒径的震区泥石流容重计算方法
本专利技术属于泥石流防治工程设计应用
，特别涉及一种基于泥石流固体颗粒粒径的震区泥石流容重计算方法。
技术介绍
“5.12”汶川地震以后，强震区内诱发了大量的地震次生灾害，强震区内每年都有不同规模的泥石流事件暴发，其中在2008年的“9.26”、2010年的“8.13”、2011年的“7.03”、2013年的“7.10”以及2016年的“7.05”强降雨过程中，强震区暴发了多期大规模的群发性泥石流灾害，给强震区带来了严重的人员伤害和巨大的经济损失，也给震后重建项目带来了巨大的破坏，而在震后的10年至15年内也将是泥石流灾害的高发期。为了保证汶川震区经济的可持续发展，泥石流防治工程的有效防治就显得十分必要。泥石流容重特征是泥石流动力学参数的主要参数之一，也是泥石流防治工程设计的主要参数之一。例如，2013年“7.10”特大暴雨期间，羊店1组后山坡面泥石流3条沟道均发生泥石流，泥石流物源以碎石、块石和砂粒堆积为主，粒径范围为10～30cm，砂及粉粘土充填，通过现场配方法和查表法分别得到泥石流的容重为19.6KN/m3和17.2KN/m3，通过黏粒百分含量和粗颗粒百分含量得到的容重约为16.8KN/m3。2008年“9·26”洪灾和2010年“8·14”洪灾期间牛圈沟发生大规模泥石流，其中侵蚀的物源量分别为17.15×104m3和14.87×104m3，牛圈沟泥石流松散固体物源中以碎块石为主，通过现场配方法得到泥石流的容重为17.8KN/m3，通过黏粒百分含量和粗颗粒百分含量得到的容重约为16.4KN/m3。泥石流的容重反映了泥石流的流体的固体颗粒含量和水分含量之间的相互关系，不同的泥石流容重反映了泥石流的内部受力条件、搬运条件和泥石流暴发的降雨条件等，现行的泥石流容重主要根据泥石流固体颗粒(＞2mm)的累积百分含量和黏粒(＜0.05mm)的累积百分含量。汶川震区泥石流沟物源主要为松散的岩土体颗粒，黏粒含量极低，基于泥石流粗颗粒和黏粒累积百分含量的泥石流容重计算方法导致泥石流的容重相对于其实际容重值偏小，使得泥石流的防治工程设计过程中的动力学参数过低，不能满足震区泥石流防治工程设计的需要。
技术实现思路
本专利技术的目的在于：针对上述存在的问题，提供一种基于泥石流固体颗粒粒径的震区泥石流容重计算方法，解决强震区泥石流暴发过程中固体颗粒与泥石流容重的关系，构建基于泥石流固体颗粒粒径的泥石流容重预测模型，并进行实例运用，为泥石流预警预报和防治提供新方法，适用于泥石流防治实际工程的需要。本专利技术技术的技术方案实现方式：一种基于泥石流固体颗粒粒径的震区泥石流容重计算方法，所述震区的泥石流沟包括沟道和在沟道上运动的泥石流体，其特征在于：将震区泥石流的固体颗粒粒径应用于泥石流容重计算，所述震区泥石流容重预测计算方法如下：A.获取震区泥石流沟的泥石流体流深h(m)，通过泥石流沟大比例尺测绘得到沟道的沟道坡度α(°)，沟道纵比降J(‰)，通过泥石流采样和筛分得到泥石流体固体颗粒粒径di(mm)，通过泥石流沟床状况确定泥石流体糙率系数n，通过为泥石流暴发时的野外温度调查确定泥石流体的粘滞系数η(Pa·s)；B.通过以下公式确定震区泥石流浆体密度ρ(Kg/m3)公式表示为：式中：di为累积百分含量为i的泥石流固体颗粒粒径(mm)，η为泥石流的粘滞系数(Pa·s)，ρ为泥石流浆体的密度(kg/m3)，g为重力加速度(m/s2)，γS为泥石流固体颗粒容重(KN/m3)，H为泥石流沟道堆积物的堆积厚度(m)，n为泥石流沟道的糙率系数；C.通过以下公式确定震区泥石流容重γL(KN/m3)公式表示为：γL＝C(γS-γ)+γ,式中：C＝(ρ-ρW)/(ρS-ρW)式中：C为泥石流体积浓度(％)，ρw为水体的密度(kg/m3)，ρ为泥石流浆体的密度(kg/m3)，γ为泥石流浆体容重(KN/m3)，γ＝ρg；γS为泥石流固体颗粒容重(KN/m3)，γS＝ρSg；γL为泥石流容重(KN/m3)。本专利技术所述的基于泥石流固体颗粒粒径的震区泥石流容重计算方法，其所述步骤C中确定震区泥石流容重的公式适用于震区泥石流沟的防治工程设计，将得到的泥石流容重作为泥石流防治工程设计的最小泥石流容重参数。本专利技术所述的基于泥石流固体颗粒粒径的震区泥石流容重计算方法，其所述步骤C中确定震区泥石流容重的公式适用于汶川震区黏粒小于0.05mm含量，泥石流固体颗粒大于2mm级配宽的泥石流沟的泥石流容重预测。与现有技术相比，本专利技术的有益效果是：本专利技术突破传统经验统计模型的限制，从泥石流颗粒在泥石流体中的受力特征入手，通过泥石流固体颗粒在泥石流体的运移时的静力极限平衡条件构建了沿流向方向和垂直流向方向的受力特征，充分考虑了因为泥石流不同容重所形成的不同粘滞性条件对泥石流颗粒运移影响，因地制宜的反映了汶川震区的低粘粒含量泥石流中的固体颗粒粒径与泥石流容重间的相互关系，计算理论完善，计算结果符合泥石流现场调查，能够为汶川震区泥石流防治工程设计参数提供技术支持，避免因为容重设计值低所导致的泥石流防治工程未能满足实际需要，进而造成泥石流防治工程失效。附图说明图1是泥石流沟道堆积物受力分析示意图。图中标号如下：①为沟道，②为泥石流体，H为堆积层厚度，h为泥石流流深，vZ为沿流向方向的泥石流流速，dN为固体颗粒所受的支持力，dGdh为泥石流流体颗粒的自重应力，dGi为泥石流固体颗粒的自重应力，dτh为沿流向方向的泥石流单位剪应力。具体实施方式下面结合附图，对本专利技术作详细的说明。为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本专利技术技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本专利技术，并不用于限定专利技术。一种基于泥石流固体颗粒粒径的震区泥石流容重计算方法，所述震区的泥石流沟包括沟道①和在沟道上运动的泥石流体②，将震区泥石流的固体颗粒粒径应用于泥石流容重计算，其主要思路是：首先根据震区泥石流固体颗粒在泥石流体中的受力平衡绘制泥石流固体颗粒受力示意图，分析泥石流固体颗粒运动过程中的驱动力条件和阻力条件；然后根据理论力学推导泥石流固体颗粒运动时的泥石流容重；再通过震区已发生泥石流事件的泥石流沟的调查、采样、筛分以及泥石流特征参数的测量等手段确定泥石流流深、沟道坡度、沟道纵比降、泥石流固体颗粒粒径泥石流糙率系数泥石流粘滞系数；将所得参数带入震区泥石流容重计算模型中，得到泥石流容重预测值。所述震区泥石流容重预测计算方法及步骤如下：首先，根据震区泥石流固体颗粒在泥石流体中的受力平衡绘制泥石流固体颗粒受力示意图，分析泥石流固体颗粒运动过程中的驱动力条件和阻力条件。假设震区泥石流沟道纵比降为J，沟道坡度为α，泥石流固体颗粒粒径为di，如图1所示，泥石流固体颗粒受到的支持力dfB、固体颗粒自重应力dGS、泥石流体自重应力dGA、泥石流体剪切应力τ公式表示为：式中：dfB为固体颗粒所受的支持力(Pa)，dGS为泥石流流体颗粒的自重应力(Pa)，dGA为泥石流固体颗粒的自重应力(Pa)，τ为流体剪切应力(Pa)，γ为泥石流浆体容重(KN/m3)，γS为泥石流固体颗粒的容重(KN/m3)，VG为泥石流因固体颗粒挤排本文档来自技高网...

【技术保护点】
一种基于泥石流固体颗粒粒径的震区泥石流容重计算方法，所述震区的泥石流沟包括沟道和在沟道上运动的泥石流体，其特征在于：将震区泥石流的固体颗粒粒径应用于泥石流容重计算，所述震区泥石流容重预测计算方法如下：A.获取震区泥石流沟的泥石流体流深h(m)，通过泥石流沟大比例尺测绘得到沟道的沟道坡度α(°)，沟道纵比降J(‰)，通过泥石流采样和筛分得到泥石流体固体颗粒粒径di(mm)，通过泥石流沟床状况确定泥石流体糙率系数n，通过泥石流暴发时的野外温度确定泥石流体的粘滞系数η(Pa·s)；B.通过以下公式确定震区泥石流浆体密度ρ(KN/m

【技术特征摘要】
1.一种基于泥石流固体颗粒粒径的震区泥石流容重计算方法，所述震区的泥石流沟包括沟道和在沟道上运动的泥石流体，其特征在于：将震区泥石流的固体颗粒粒径应用于泥石流容重计算，所述震区泥石流容重预测计算方法如下：A.获取震区泥石流沟的泥石流体流深h(m)，通过泥石流沟大比例尺测绘得到沟道的沟道坡度α(°)，沟道纵比降J(‰)，通过泥石流采样和筛分得到泥石流体固体颗粒粒径di(mm)，通过泥石流沟床状况确定泥石流体糙率系数n，通过泥石流暴发时的野外温度确定泥石流体的粘滞系数η(Pa·s)；B.通过以下公式确定震区泥石流浆体密度ρ(KN/m3)，式中：di为累积百分含量为i的泥石流固体颗粒粒径(mm)，η为泥石流的粘滞系数(Pa·s)，ρ为泥石流浆体的密度(kg/m3)，g为重力加速度(m/s2),h为泥石流...

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