粘性泥石流运动流速与流量计算.doc

粘性泥石流运动流速与流量计算舒安平1，费祥俊2（1. 中国水利水电科学研究院； 2. 清华大学 ）摘要： 粘性泥石流作为最常见的一种泥石流，普遍存在于固体物质组成松散、降雨持续集中的陡峻 山区作者首先对现有的粘性泥石流运动速度的有关成果进行扼要评述根据曼宁公式的结构形式， 通过对大量泥石流沟的实测资料进行统计分析， 得出涉及参数较为全面、 具有一定普 遍意义的粘性泥石流运动速度公式， 经验表明该公式的可靠度令人满意， 据此提出粘性泥石流的 流量公式及其计算方法，从而为泥石流灾害治理工程规划设计提供了科学依据关键词： 粘性泥石流；曼宁公式；阻力参数；流速；流量1 引言我国是一个泥石流频发的国家， 特别是在西南和西北山区， 每年雨季由于滑 坡泥石流等山地灾害造成的人员伤害及经济损失均十分严重就在 2002年5 月至 8 月中旬不足 4 个月的时间里，仅仅云南省 14 个地市州因连降暴雨，引发受 灾超过 2 000 万人、死亡 230多人、经济损失高达 36 亿元的重大洪水泥石流灾 害即为一个典型的例证 可见，泥石流作为一种破坏性自然灾害， 由于来势迅猛、 影响深远， 一旦成灾，其后果相当严重， 因此对泥石流运动流速与流量进行分析 研究，不仅为工程规划设计所急需， 而且也是工程界和学术界普遍关注的重点课 题。

根据固体物质颗粒组成， 泥石流一般可分为泥流、 粘性泥石流和水石流三种 类型，其中粘性泥石流由于其固体颗粒组成范围广， 并具有粗颗粒多、 细颗粒含 量大、颗粒分布呈 “双峰 ”形态的特点，是自然界中最为常见的一种泥石流，一直 是许多专家学者研究的重点课题 目前盛行的粘性泥石流运动理论模型， 主要基 于两相流体内部阻力特点来求解泥石流运动流速 尽管这种方法理论性较强， 但 由于各种模型应用时存在着一定的局限与不足， 特别是由于粘性泥石流体内部阻 力的复杂性而不得不假定固体颗粒呈均匀分布， 加之模型中存在着一些难以确定 的参数等问题，使得目前人们提出的一些泥石流运动速度模型及流量计算方法尚 难达到实用水平 [1]在这种情况下，只能借助经验方法建立适合某种特定条件下的泥石流运动流速公式，进而求解流量，显然这些公式带有很强的经验性， 其 适用范围有限而难以推广应用有鉴于此，作者试图基于曼宁公式的结构形式， 通过对我国西部大量泥石流 沟的实测资料进行统计分析，提出涉及参数较为全面、具有一定普遍意义的粘性 泥石流运动速度公式及流量的计算方法 本项研究成果不仅有利于丰富河流动力 学学科内容，而且具有重要的实用价值 。

2现有成果的概述目前，有关粘性泥石流运动速度的研究成果较多， 这里受篇幅限制，只将具有代表性的成果进行分析云南省东川蒋家沟是一个典型的泥石流沟，每年洪水季节粘性泥石流频发， 素有天然泥石流博物馆之称，闻名于世康志成等⑵根据60年代以来收集有关 蒋家沟泥石流实测资料，通过分析泥石流内部与外部阻力变化规律， 建立阻力综合参数1/nc与泥深h之间的关系，获得了具有典型代表意义的粘性泥石流流速 公式，即n c=0.035h°.34 (1)或sc — -0.34 z. 2/3 1/2、 sc 一 0.327 .0.5 .、nc= 28.5h (h J )=28.5h J (2)我们根据云南大盈江浑水沟泥石流的观测资料 ：3：，通过分析拟合也得到类似于式(1)的阻力关系式%= 0.051h0.25 (3)此外，基于河流水力学流速公式而得到的粘性泥石流流速经验公式也不少，式（1）和式（3）均表明，泥石流阻力因子 nc值随着泥深的增加而增大一方面可理解为，泥深的增加使得挟带的固体含量也越多，因而阻力也越大 但从另一角度来看， 一定坡降下泥深越大，反映泥石流的流量也大，而流量增加为什么会导致综合阻力因子 nc也增加，令人难以理解。

而且进一步比较式（ 1）与（3）还可发现，在相同泥深 h下，浑水沟nc值比蒋家沟要大，因为同为粘性泥石流，蒋家沟的细颗粒含量比浑水沟大，而上 述经验式的系数又不能正确反映颗粒组成及比降对 nc值的影响，表明经验关系式所反 映的因子并不全面此外，通过分析甘肃省武都地区柳湾沟、 火烧沟及泥湾沟的大量粘性泥石流 实测资料⑷进一步发现，1/nc与h之间不仅并不存在简单的相关关系，而且按 式（2）计算的泥石流流速普遍高出实测值约一倍左右 :1：，充分表明具有式（1）或式（2）结构形式的粘性泥石流流速公式存在着结构形式简单、包涵因子不全 面等弊端，而且经验性很强，因此有必要进一步深入分析研究3粘性泥石流的流速3.1流速公式的建立 一般来说，应用曼宁公式来表达泥石流流速，其中的阻力 因子nc实际上是包含泥石流 内部及外部的阻力综合因子泥石流内部阻力由颗 粒间相互作用及液体粘性引起，与由边壁糙度引起的外部阻力完全不同 对于粘性泥石流的外部糙度，据我国泥石流界估算结果，其平均值 1/n -30，即曼宁糙率n=0.033左右对于泥石流内部阻力综合参数，可以通过引入修正系数a对曼 宁公式进行修正，于是得到泥石流流速公式nc =1/n 1/a h2/3J1/2=1/nc h2/3J1/2 ⑷根据对大量的实测资料分析结果，泥石流内部阻力综合参数 1/n c与泥石流固体浓度、颗粒组成以及泥深与坡降等因素有关。

通过分析包含蒋家沟在内的几条粘性泥石流沟的实测资料，首先发现浓度对 nc的影响是复杂的，一方面粘性泥石流浓度S的增加使得内部阻力提高；另一方面当 S很高时，由于颗粒沉速 与内部阻力因浓度提高迅速减少，因而浓度对泥石流内部阻力有着双重作用，宜 用S（1-Sv）作为影响因子其次，对于粘性泥石流流速，如采用修正的曼宁公式 来表达，则式（4）中水深h及坡降J的方次明显偏高，这说明采用a对曼宁公 式进行修正，应使h及J的方次下降尽管如此，但不同 泥石流沟的阻力关系还是无法统一在一起，进一步对比分析蒋家沟及浑水沟的实测资料，发 现蒋家沟粘性泥石流的细颗粒含量要比浑水沟多，而在相同条件下 1/n c值要大，表明1/n c与细颗粒含量呈反变关系，可选择反映颗粒级配特征的下限粒径 dio作为细颗粒含 量的指标由于各组次泥石流的dio值不完全相同，这里只好取各组次的 平均值，蒋家沟粘性泥石流平均 dio=O.O11mm浑水沟粘性泥石流平均 dio=O.O32mm综上分析，泥石流内部阻一；* 二：关系密切应用蒋家沟及浑水沟的力1/nc与水力泥沙综合因子 实测资料，在对数坐标上点绘1/nc与©的关系图，发现两沟的点据能够很好地集中在一条直线的两侧， 其线性关系令人满意，如图1所示。

对图1的关系进一步拟合分析，得到如下综合阻力参数公式⑸« L枫％」式中dio单位为mm将式（5）关系代入式（4），得到粘性泥石流流速公式―丄為•壯U］加 ⑹叫 L % _式（6）作为粘性泥石流流速公式，其形式并不复杂，而且式中所包涵的影响参数比较全面，与同类公式相比具有更加普遍意义和实用价值3.2流速公式的验证应用甘肃武都地区三条粘性泥石流沟的资料对式（6）进行检验，结果表明 按式（6）计算的泥石流流速与实测值相吻合的程度令人满意，如图 2所示此外，图2同时还点绘有别于前文建立式（6）时的1966年蒋家沟实测资料，我们 发现蒋家沟 点群也均集中在45°线的两侧，基本说明式（6）可以用于不同地区 的泥石流流速计算鉴于目前所收集到的实测资料十分有限，以上的验证只是初 步的，今后有待于进一步丰富完善 图1取1/ n c关系图2粘性泥石流流速公式的验证结果Fig.1 Relati on ship betwee nO andFig.2 Verification of velocity of viscous debris1/n cflow4粘性泥石流的流量泥石流流量作为泥石流治理工程规划设计的重要参数，可以基于式( 6)所示的泥石流流速公式进行计算根据泥石流流量连续律，可得粘性泥石流沟单宽流量计算公式““严丄嘉二1瑟空f审 ⑺叫 L % 」对于一般泥石流沟，泥石流流量不仅取决于流速的大小， 而且还与沟道断面 形态密切相关。

如定义泥石流沟断面形态参数或称沟相关系系数 M为断面湿周P与水力半径R之比，即M=P/R=A/F2 (8)根据流量连续方程，可得泥石流流量为q=au=mRuQ/M=RUC (9)5应用计算5.1流速与流量计算为便于应用，现列举体积浓度S=0.36、下限粒径di°=0.032mm勺粘性泥石流 沟进行计算说明作为一种简化处理，假定具有宽浅特征的泥石流沟床宽度 B=20m和纵坡J=0.06保持不变，则将式(6)中的泥深h可以近似用水力半径R 代替来计 算泥石流流速UC首先将泥深取值为h=0.5、1.0、1.5、2.0m,依次 列表计算R、M UC及F2Uc(=Q /M)，然后绘制出R〜UC及R〜Q关系线，分别如表 1和图3所示因此，只要已知该泥石流沟的流量 Q流速UC及水力半径R中的任意一个参数，即可利用图3所示的关系线求得其他两个变量例如，Q=100m/s， 由图3查得泥石流沟的水力半径 R=113m再由R〜UC关系求得流速UC = 3. 94m/s； 如此类推，可求得 Q=50ri/s时水力半径R=0 69m流速U=3.30m/s，这表明，在 相同浓度和纵坡的前提下，尽管泥石流流 量变幅很大，但由于断面水力半径的调整，造成泥石流的运动流速变幅却比较小。

表1粘性泥石流的流速与流量计算结果( Sv=0.36，d10=0.032mm)Table 1 The calculated resultsfor velocityflowand discharge ofviscous debrish/m0.51.01.52.0备注R=Bh/(B+2h)0.480.911.301.67M=Bh/R44.1024.2017.6314.40UC/m s-12.953.664.134.48式(6)Q/M0.673.037.0312.45Q/m3 s-129.5173.22123.88179.245.2流速计算方法的改进我们仍以蒋家沟泥石流为例进行计算分析，提出对粘性泥 石流流速公式（6）进一步完善 的方法通过分析更多的蒋家沟 实测资料，结果发现当泥石流容 重减小时，颗粒组成相应变化 用固定的dio平均值，作为细颗 粒泥沙特征粒径计算流速结果 误差有所加大因此建议当 ym2. 12时，仍采用10 耳畑① 50 100 20010 50 1A0 却0宀1图3 R〜UC及R〜Q关系线dio=O.1653.6O(1O)Fig.3 Relati on ship betwee n RUc and R 〜Qdi0=O.O11mm 当 T< 2.12 时 经资料分析dio有所增加，建议 采用如下关系按上式取d1O值，我们又用蒋家沟较早时期的观测资料对式（6）进行校核计算，验算结果表明， 计算泥石流流速的平均误差为 芳.82%,其中误差在10%以内的 比例达到80%^上，表明修正后 的流速公式计算结果令人满意。

值得说明的是，式（10）只适用于蒋 家沟泥石流流速计算，对于其他泥石流沟尚需具体问题具体分析6结语（1）根据修正的曼宁公式，通过分析大量的实测资料，提出粘性泥石流综合阻力参数公式 （5）及流速公式（6），而且流速公式（6）具有考虑影响因子 全面、结构形式不复。