本发明属于泥石流防治工程设计应用技术领域，特别涉及一种基于泥石流固体颗粒粒径的震区泥石流容重计算方法。
背景技术：
“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)、泥石流固体颗粒容重γs(kn/m3)，通过泥石流沟床状况确定泥石流体糙率系数n，通过为泥石流暴发时的野外温度调查确定泥石流体的粘滞系数η(pa·s)；
b.通过以下公式确定震区泥石流容重γ(kn/m3)，
γ＝(4γsπdmax-9πηv)/(4π+6cosα)
式中，dmax为泥石流最大固体颗粒的直径(mm)，v为山洪泥石流的流速(m/s)，η为流体的粘滞系数(pa·s)，α为泥石流沟的沟道坡度(°)。
本发明所述的基于泥石流固体颗粒粒径的震区泥石流容重计算方法，其所述步骤b中确定震区泥石流容重的公式适用于震区泥石流沟的防治工程设计，将得到的泥石流容重作为泥石流防治工程设计的最小泥石流容重参数。
本发明所述的基于泥石流固体颗粒粒径的震区泥石流容重计算方法，其所述步骤b中确定震区泥石流容重的公式适用于汶川震区低黏粒小于0.05mm含量，泥石流固体颗粒大于2mm级配宽的泥石流沟的泥石流容重预测。
与现有技术相比，本发明的有益效果是：本发明突破传统经验统计模型的限制，从泥石流颗粒在泥石流体中的受力特征入手，通过泥石流固体颗粒在泥石流体的运移时的静力极限平衡条件构建了沿流向方向和垂直流向方向的受力特征，充分考虑了因为泥石流不同容重所形成的不同粘滞性条件对泥石流颗粒运移影响，因地制宜的反映了汶川震区的低粘粒含量泥石流中的固体颗粒粒径与泥石流容重间的相互关系，计算理论完善，计算结果符合泥石流现场调查，能够为汶川震区泥石流防治工程设计参数提供技术支持，避免因为容重设计值低所导致的泥石流防治工程未能满足实际需要，进而造成泥石流防治工程失效。
附图说明
图1是泥石流沟道堆积物受力分析示意图。
图中标号如下：①为沟道，②为泥石流体，hy为泥石流固体颗粒所处的流深，h为泥石流的流深，dfy为泥石流的黏滞力，dff为泥石流的浮应力，dp为泥石流的压强差，dgy为泥石流固体颗粒的自重应力。
具体实施方式
下面结合附图，对本发明作详细的说明。
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定发明。
一种基于泥石流固体颗粒粒径的震区泥石流容重计算方法，所述的震区泥石流沟包括沟道①和在沟道上运动的泥石流体②，将震区泥石流的固体颗粒粒径应用于泥石流容重计算，其主要思路是：首先根据震区泥石流固体颗粒在泥石流体中的受力平衡绘制泥石流固体颗粒受力示意图，分析泥石流固体颗粒运动过程中的驱动力条件和阻力条件；然后根据理论力学推导泥石流固体颗粒运动时的泥石流容重；再通过震区已发生泥石流事件的泥石流沟的调查、采样、筛分以及泥石流特征参数的测量等手段确定泥石流流深、沟道坡度、沟道纵比降、泥石流固体颗粒粒径泥石流糙率系数泥石流粘滞系数；将所得参数带入震区泥石流容重计算模型中，得到泥石流容重预测值。
所述震区泥石流容重预测计算方法及步骤如下：
首先，根据震区泥石流固体颗粒在泥石流体中的受力平衡绘制泥石流固体颗粒受力示意图，分析泥石流固体颗粒运动过程中的驱动力条件和阻力条件。假设震区泥石流沟道纵比降为j，沟道坡度为α，泥石流固体颗粒粒径为di，如图1所示，泥石流固体颗粒受到的自重应力dgy、浮应力dff、黏滞力dfy公式表示为：
式中，γs为泥石流固体颗粒容重(kn/m3)，γ为泥石流体容重(kn/m3)，dmax为泥石流体中的最大固体颗粒粒径(mm)，η为泥石流的粘滞系数(pa·s)。
然后，山洪泥石流运动时，由于沿其流向的流速等差异性，导致压差阻力的存在，根据伯努利效应得到山洪泥石流固体颗粒上下压强差异形成的压差阻力公式表示为：
dp＝p2-p1＝ρg(h1-h2)＝2γdmax………………………(2)
式中，ρ为泥石流的密度(kg/m3)，h1和h2为泥石流固体颗粒上覆面和下伏面的流深(m)，g为重力加速度(m/s2)。
然后，假设山洪泥石流固体颗粒沿其垂直流向方向的速度差异为零，其加速度ay为零，则山洪泥石流的容重与固体颗粒的关系公式表示为：
ay＝dgy-dff-dfy-cosαdp＝0………………………(3)
即：
γ＝(4γsπdmax-9πηv)/(4π+6cosα)……………………(4)
式中，dmax为泥石流最大固体颗粒的直径(mm)，v为山洪泥石流的流速(m/s)，α为泥石流沟的沟道坡度(°)，η为流体的粘滞系数(pa·s)，在常温20℃时，η＝1.005(pa·s)；其他符号同上。
其中，所述步骤b中确定震区泥石流容重的公式适用于震区泥石流沟的防治工程设计，将得到的泥石流容重作为泥石流防治工程设计的最小泥石流容重参数。
其中，所述步骤b中确定震区泥石流容重的公式适用于汶川震区低黏粒小于0.05mm含量，泥石流固体颗粒大于2mm级配宽的泥石流沟的泥石流容重预测。
以下为本发明的具体实施例：
实施例一：
牛圈沟位于映秀镇南侧，岷江右岸，往上游至映秀新镇约1km。牛圈沟流域形态为扇形，流域纵向长度5.4km，平均宽度2.0km，流域面积10.7km2。流域最高点位于勘查区西侧，高程约2700m，最低点位于牛圈沟汇入岷江河口，高程858.3m，相对高差约1842m。沟道总体上以深切割“v”型谷为主，主沟纵长6124m，相对高差1629m，主沟平均纵比降266‰。
为了有效的防治和治理泥石流，防止岷江被堵塞威胁下游安全，提前获得泥石流暴发时容重，根据泥石流物源级配特征进行有效预测，具体计算方法及步骤如下：
a.根据牛圈沟沟口位置的野外调查、采样和泥石流参数测试，确定了牛圈沟沟口的平均泥深0.7m，泥石流沟道纵比降130‰，糙率系数11，固体物质重度26.5(kn/m3)，泥石流固体颗粒粒径dmax为1.3m，粘滞系数为1.005(pa·s)，沟口泥石流流速为3.15(m/s)。
b.将a步骤确定的各参数代入下公式，
γ＝(4γsπdmax-9πηv)/(4π+6cosα)
c.通过公式计算获取泥石流暴发时的容重为18.5kn/m3，则牛圈沟泥石流暴发时的泥石流容重为18.5kn/m3。
实施例二：
高家沟位于银杏乡北侧，岷江右岸，主沟纵向长3.26km，横向宽1.3km，沟域面积3.79km2，平面形态呈近似扇形，流域相对高差约1807m。高家沟纵坡总体较陡，全沟域平均纵比降523‰。沟域内共发育4条支沟，支沟长度较短，一般长度在0.5～1.6km之间，纵坡较大，多在500‰以上。
为了有效的防治和治理泥石流，防止岷江被堵塞威胁下游安全，提前获得泥石流暴发时容重，根据泥石流物源级配特征进行有效预测，具体计算方法及步骤如下：
a.根据高家沟口支沟位置的野外调查、采样和泥石流参数测试，确定了沟口的平均泥深1.1m，泥石流沟道纵比降233‰，糙率系数13，固体物质重度26.5(kn/m3)，泥石流固体颗粒粒径dmax为1.6m，粘滞系数为1.005(pa·s)，沟口泥石流流速为6.68(m/s)。
b.将a步骤确定的各参数代入下公式，
γ＝(4γsπdmax-9πηv)/(4π+6cosα)
c.通过公式计算获取泥石流暴发时的容重为19.5kn/m3，则高家沟口的泥石流容重为19.5kn/m3。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。