泥石流动力特征计算.docx

1、泥石流动力特征计算3.3.2泥石流的力学特征(1)容重泥石流静力学特征主要指泥石流体或浆体的容重、含水量、物质组成、流变特征、化学性质及其静力特征等。在一般地区出于泥石流的突发性、冲击力大等条件所限，难以直接测得天然泥石流容重。一般采用现场调查试验法进行泥石流容重的测定，即在现场请当地亲眼看见泥石流暴发的居民多人，在需要测试的沟段，选取有代表性的堆积物搅拌成暴发时泥石流流体状态，进行样品鉴定，然后分别测出样品的总质量和总体积，求出泥石流流体容重。在无法取得代表性样品时，根据规范中泥石流沟易发程度数量化评分标准，对某泥石流沟进行泥石流沟易发程度数量化评分（详见表4-3），按照规范中附表“数量化评

2、分（N）与重度、（1）关系”，可以得到泥石流的容重。本报告采用后一种方法，查表得到泥石流的容重为1.68t/m3。(2)泥石流流速泥石流的流量是泥石流重要的特征值之一。它不仅反映了泥石流的强度，规模和流体性质，而且决定着防治泥石流工程建筑物的类型、结构和尺寸。因此，泥石流的洪峰流量是泥石流研究和防治工程中不可缺少的参数。流速VC 按照铁道部推荐的稀性泥石流的计算公式进行计算:式中： 泥石流中含沙量变化引起的流速修正系数，；R水力半径(m)，2.5m；IC泥石流水力坡度（），用沟床纵坡代替；清水河床糙率系数；泥石流泥沙修正系数，；泥石流容重（t/m3），为1.68t/m3；清水容重（t/m3），

3、1.0t/m3；泥石流中固体物质重度（t/m3）。根据以上计算公式，泥石流的平均流速为8.28m/s。(3)泥石流流量泥石流流量计算，目前主要有两种方法，一是雨洪法；二是形态调查法。雨洪法假设泥石流与暴雨同频率、且同步发生，先按水文方法计算出断面不同频率下的小流域暴雨洪峰流量 (计算方法查阅四川省水文手册)，然后选用堵塞系数，按下式计算泥石流流量：式中： 频率为P的泥石流洪峰值流量（m3/s）；频率为P的暴雨洪水设计流量（m3/s）；泥石流泥沙修正系数，查规范附表，值为0.71；泥石流容重（t/m3），为1.68t/m3；清水的重度（t/m3），为1.0； 泥石流中固体物质重度；DC 泥石流堵

4、塞系数（见表3-8），可查经验表为1.5。表3-8 泥石流阵流堵塞系数DC值表堵塞程度严重堵塞中等严重堵塞轻微堵塞无堵塞Dc值2.52.51.51.51.11.0按照雨洪法，利用泥石流流量公式计算所得的泥石流最大流量如下表（表3-9）所示：表3-9 雨洪法计算最大流量 沟名项目某泥石流沟流域面积F (km2)35沟长L(km)11.25平均坡降J()222泥沙修正系数0.71设计频率p(%)25暴雨洪峰流量QP(m3/s)87.49 72.54泥石流峰值流量Qc（m3/s）224.43186.06形态调查法在泥石流沟道中选择2-3个测流断面，仔细查找泥石流过境后留下的痕迹，然后测量这些断面上的

5、泥石流流面比降（若不能由泥痕确定，则用沟床比降代替）、泥位高度Hc（或水力半径）和泥石流过流断面面积等参数。用相应的泥石流流速计算公式，求出平均流速Vc后，即用下式求泥石流断面峰值流量Qc：式中： Qc泥石流断面峰值流量（m3/s）；Wc泥石流过流断面面积（m2），取流通区平均沟道宽8 m，水力半径为2.5 m，过流断面面积为20 m2；Vc泥石流断面平均流速（m/s），为8.28m/s。用上式计算的泥石流断面峰值流量Qc为165.66 m3/s。(4)一次泥石流过程总量计算一次泥石流总量Q可通过计算法和实测法确定。实测法精度高，但因往往不具备测量条件，只是一个粗略的概算。计算法根据泥石流历时

6、T（s）和最大流量Qc（m3/s），按泥石流暴涨暴落的特点，将其过程线概化成五角形，按下式计算Q（m3），泥石流历时T取1200 s： Q = 0.264TQc=KTQc一次泥石流冲出的固体物质总量QH（m3）：根据该公式，在计算频率为20年一遇的情况下，某泥石流沟一次泥石流冲出的最大固体物质总量为41850m3；在计算频率为50年一遇的情况下，一次泥石流冲出的最大固体物质总量为50480m3。(5)泥石流整体冲压力泥石流冲击力是泥石流防治工程设计的重要参数。采用规范推荐的泥石流体整体冲压力计算公式：式中： 泥石流整体冲压力；泥石流容重，1.68t/m3；泥石流流速，8.28m/s；建筑物受力

7、面与泥石流冲击压力方向的夹角，取90o；重力加速度；建筑物形状系数，圆形建筑物取1.0，矩形建筑物取1.33，方形建筑物取1.47。本次建筑物主要以矩形为主，取1.33。通过该公式的计算，泥石流的整体冲压力约15.64104Pa。(6)泥石流体中大石块的最大流速按照规范，用以下经验公式计算：式中： Vs泥石流中大石块的移动速度（m/s）；全面考虑摩擦系数、泥石流容重、石块比重、石块形状系数、沟床比降等因素的参数。3.54.5，平均=4.0；dmax泥石流堆积物中最大石块粒径（m），为2.5m。通过计算，泥石流大石块的最大速度为6.32m/s。(7)泥石流中石块冲击力泥石流中石块冲击力的计算参照

8、以下公式（吴积善，泥石流及其综合治理1993）：式中： Pd为泥石流中石块的冲击力（Pa）；动能折减系数，对于圆端正面撞击，采用=0.3；Vc泥石流平均流速（m/s），为8.28m/s；Q石块质量（kg），按最大粒径2.5m来算；受力面与泥石流冲击力方向的夹角，取90；C1、C2巨石及拦挡圬工的弹性变形系数，C1+C2=0.005。根据该公式，泥石流中石块的冲击力为14.35Kpa。(8)泥石流最大冲起高度与爬高泥石流最大冲起高度根据规范推荐公式为：由于泥石流在爬高过程中受到沟床阻力的影响其爬高为：式中：为迎面坡的函数。通过计算，泥石流的最大冲起高度为3.5m，爬高为2.8m。根据以上的计算，

9、将泥石流各动力参数的计算结果汇总于表3-11。表3-11 泥石流动力学参数计算结果参数结果参数结果泥石流平均流速（m/s）8.28泥石流断面设计流量（形态调查法）(m3/s)229.46大石块运动速度(m/s)6.3220年一遇清水洪峰流量（雨洪法）（m3/s）72.54容重（t/m3）1.6820年一遇泥石流洪峰流量（雨洪法）（m3/s）186.06泥石流中石块冲击力（Kpa）14.3520年一遇一次泥石流可能最大冲出量（雨洪法）（m3）41850泥石流整体冲压力（Pa）15.6410450年一遇清水洪峰流量（雨洪法）（m3/s）87.49爬高（m）2.850年一遇泥石流洪峰流量（雨洪法）（

10、m3/s）224.43最大冲起高度（m）3.550年一遇一次泥石流可能最大冲出量（雨洪法）（m3）50480(3)、泥石流动力学参数计算 A、流速计算：据勘查所得泥石流流体水力半径、纵坡、沟床糙率及重度等参数计算；也可按泥石流的性质和所在地域，选择适合的地区性经验公式计算。泥石流流速是决定泥石流动力学性质的最重要参数之一。目前泥石流流速计算公式为半经验或经验公式，概括起来一般分为稀性泥石流流速计算公式、粘性泥石流计算公式和根据泥石流中大石块运动速度推算泥石流流速等三种办法。a、 稀性泥石流流速计算公式VC =1/aRIC (221)式中： VC 泥石流断面平均流速（ms）；1/a =1 / (

11、H 1)1/2 -泥石流中由含沙量变化而引起的流速修正系数，查表23-1R 水力半径（m），一般可用平均水深H（m）代替；IC 泥石流水力坡度（），一般可用沟床纵坡代替。 清水河床糙率系数，查当地水文手册或查铁路桥渡勘测设计规范（TBJ17-86）。表22-1； 泥石流河道河床清水糙率表河段特征河床物质组成及形态沟岸状况1 / n河段较顺直，河宽逐渐扩展，断面较规则，水流通畅，没有堵塞。砂、土质河床，河底平顺，沟床纵坡 30平顺的土岸或人工堤防，沟岸山坡坡度 120。崩、坡积物或岩质沟岸，沟岸山坡坡度 320。，坍塌很发育，有明显堵塞痕迹。25（1535）山区峡谷河段、弯曲河段、阻塞的复式断面

12、河段砂、园砾石河床，边滩沙州犬牙交错人工堤防强制弯曲者35（3040）有矾石或丁坝挑流者30（2535）卵石、园石河床，河底起不平，或长有水生植物的沟床。参差不齐的卵石、园石河岸，或长中密灌丛的河岸25（2030）参差不齐的岩岸，或长灌木丛生的河岸20（1525）卵石、块石、大漂石河床，石梁、跌水、孤石交错、或水生植物稠密，阻水严重的沟床0。参差不齐的岩岸，或长灌木丛生的河岸15（1220）两岸时有岩咀突出，很不平顺，形成强烈斜流、回水、死水的河岸12（1015）b、粘性泥石流流速计算公式 粘性泥石流的河床糙率，用内插法由表222查得。表222 粘性泥石流河床糙率序号泥石流体特征沟床状况糙率值

13、1流体呈整体运动；石块粒径大小悬殊，一般在3050cm，25m粒径的石块约占20；龙头由大石块组成，在弯道或河床展宽处易停积，后续流可超越而过，龙头流速小于龙身流速、堆积呈垄岗状河床极粗糙，沟内有巨石和挟带的树木堆积，多弯道和大跌水，沟内不能通行，人迹罕见，沟床流通段纵坡在100150，阻力特征属高阻型平均值0.2702m时，0.4453.572.252流体呈整体运动，石块较大，一般石块粒径2030cm，含少量粒径23m的大石块；流体搅拌较为均匀；龙头紊动强烈，有黑色烟雾及火花；龙头和龙身流速基本一致；停积后呈垄岗状堆积河床比较粗糙，凹凸不平，石块较多，有弯道、跌水；沟床流通段纵坡70100，

14、阻力特征属高阻型1.5m时，平均0.0401.5m时，0.0500.100平均0.0672030251020153流体搅拌十分均匀；石块粒径一般在10cm左右，挟有个别23m的大石块；龙头和龙身物质组成差别不大；在运动过程中龙头紊动十分强烈，浪花飞溅；停积后浆体与石块不分离，向四周扩散呈叶片状沟床较稳定，河床物质较均匀，粒径10cm左右；受洪水冲刷沟底不平而且粗糙，流水沟两侧较平顺，但干而粗糙；流通段沟底纵坡5570，阻力特征属中阻型或高阻型0.1m0.5m0.0430.5m2.0m0.0772.0m4.0m0.1002313104流体搅拌十分均匀；石块粒径一般在10cm左右，挟有个别23m的

15、大石块；龙头和龙身物质组成差别不大；在运动过程中龙头紊动十分强烈，浪花飞溅；停积后浆体与石块不分离，向四周扩散呈叶片状泥石流铺床后原河床粘附层泥浆体，使干而粗糙的河床变得光滑平顺，利于泥石流体运动，阻力特征属低阻型0.10.5m0.0220.02.0m0.0332.04.0m0.050462620c、 泥石流中石块运动速度推算泥石流流速计算公式在缺乏大量实验数据和实测数据阶情况下，为便于以堆积后的泥石流冲出物最大粒径大体推求石块运动速度推算泥石流流速的经验公式： Vskd max.（222）式中：泥石流中大石块的移动速度推算泥石流流速（ms）；泥石流堆积物中最大石块的粒径（m）；k全面考虑的摩

16、擦系数（泥石流重度、石块密度、石块形状系数、沟床比降等因素）。变化范围3.55.5，B、 流量计算：泥石流流量可采用形态调查法（据泥痕勘测所得的过流断面面积乘以流速）或雨洪法（按暴雨洪水流量乘以泥石流修正系数）确定。暴雨小径流的地区性经验公式较多，暴雨洪水流量应采用适合当地的经验公式计算。a、 形态调查法 在泥石流沟道中选择23个有代表性的过流断面。查找泥石流过境后留下的痕迹，然后确定泥位。最后测量这些断面上的泥石流流面比降（若不能由痕迹确定，则用沟床比降代替）、泥位高度HC（或水力半径）和泥石流过流断面面积等参数。用相应的泥石流流速计算公式，求出断面平均流速VC后，即可用下式求泥石流断面峰值

17、流量QC。QCWCVC .（223）式中：WC 泥石流过流断面面积（m2）；VC 泥石流断面平均流速（ms）。 雨洪法假设泥石流与暴雨同频率、且同步发生，计算断面的暴雨洪水设计流量全部转变成泥石流流量的前题下建立的计算方法。其计算步骤是先按水文方法计算出断面不同频率下的小流域暴雨洪峰流量 (计算方法查阅水文手册)，然后选用泥石流重度和堵塞系数，按式I-2计算泥石流流量。 QC =（1十）QPDC.（224）式中：QC 频率为P的泥石流洪峰值流量（m3s）；QP 频率为P的暴雨洪水设计流量（m3s）； 根据各省水文手册中给出的计算公式计算，或用下式 QB=0.278 riF 计算式中：r按小时平

18、均雨强（毫米 / 小时）设计，用实测最大小时雨强校核 iB产水系数。植被具有保水功能，降雨后入渗少，产水系数与雨强的大小和植被的多少呈正变，与松散土层的厚度呈反变，一般产水系数i=0.30.8；泥百流地区多在0.6以下为宜。产沙系数则相反，它与植被的多少呈反变，与雨强的大小和松散土层的厚度呈正变，一般产沙系数is=0.30.8； F流域面积（平方公里）（1十）泥石流中由含沙量变化而引起的流量修正系数，查表23-1。 式中 ：（C -W）/（H c）（225） C 泥石流重度（tm3）；W 清水的重度（tm3）；H 泥石流中固体物质比重（tm3）；DC 泥石流堵塞系数，可查经验表223； 表22

19、3 泥石流堵塞系数Dc值堵塞程度特 征堵塞系数Dc严重 河槽弯曲，河段宽窄不均，卡口、陡坎多。大部分支沟交汇角度大，形成区集中。物质组成粘性大，稠度高，沟槽堵塞严重，阵流间隔时间长2.5 中等 沟槽较顺直，沟段宽窄较均匀，陡坎、卡口不多。主支沟交角多小于60，形成区不太集中。河床堵塞情况一般，流体多呈稠浆稀粥状1.52.5轻微 沟槽顺直均匀，主支沟交汇角小，基本无卡口、陡坎，形成区分散。物质组成粘度小，阵流的间隔时间短而少1.11.5无1.0一次泥石流过程总量计算一次泥石流总量Q可通过计算法和实测法确定。实测法精度高，但因往往不具备测量条件，只是一个粗略的概算。计算法根据泥石流历时T（s）和最

20、大流量Qc（m3/s），按泥石流暴涨暴落的特点，将其过程线概化成五角形，按式I.6计Q = 0.264TQc= K TQc （226）式中 K 值的变化随流域面积 (F) 的大小而变化：当 F 100 (km2 )时： K = 0.0.0252；泥石流体中的固体物质流量QH（m3）： QH Q（C -W）/（H W）(227)C、 冲击力计算：可用附录I中公式计算泥石流整体冲击力、泥石流中大石块冲击力。泥石流中大石块冲击力的计算方法较多，计算结果可信度较低。a) 泥石流体整体冲击压力计算公式F ( c /g Vc2 sin) （228）式中：F泥石流体整体冲击压力（Pa）；g重力加速度（ms2

21、），取g9.8ms2；建筑物受力面与泥石流冲击力方向的夹角（）；建筑物形状系数： 圆形建筑物1.0，矩形建筑物1.33，方形建筑物：1.47。D、 弯道超高与冲高计算：泥石流流动在弯曲沟道外侧产生的超高值和泥石流正面遇阻的冲起高度。a)、泥石流的弯道超高由于泥石流流速快，惯性大，故在弯道凹岸处有比水流更加显著的弯道超高现象。根据弯道泥面横比降动力平衡条件，推导出计算弯道超高的公式：.（229）式中：h 弯道超高（m）；R2 凹岸曲率半径（m）；R1 凸岸曲率半径（m）；b)、泥石流最大冲起高度 . （2210）c）泥石流在爬高过程中由于受到沟床阻力的影响，其爬高H：.（2211）（4）泥石流的

22、形成区、流通区和堆积区测绘。 a）工程治理区实测剖面至少按一纵三横控制； b）重点区应有1-3个探槽或探坑（井）控制；5)、对各类防治工程提供以下主要设计参数。（1）各类拦挡坝: 覆盖层和基岩的重度、承载力标准值、抗剪强度，基面摩擦系数，泥石流的性质与类型，发生频次，泥石流体的重度和物质组成，泥石流体的流速、流量和设计暴雨洪水频率，泥石流回淤坡度和固体物质颗粒成分，沟床清水冲刷线。（2）其它工程：桩林着重于其锚固段基岩深度、风化程度和力学性质；排导槽、渡槽着重于泥石流运动的最小坡度、冲击力、弯道超高和冲高；导流堤、护岸堤和防冲墩着重于基岩的埋藏深度和性质、泥石流冲击力和弯道超高、墙背摩擦角；停

23、淤场着重于淤积总量、淤积总高度和分期淤积高度。2泥石流流体重度(c)的确定泥石流流体重度(c)的确定，目前测试手段还不能解决，因此暂时的替代办法有：1)、根据数量化评分法确定重度：根据数量化评分（N）建立与重度(c)的线性关系，制成备查表直接查找，参见表23-12)、现场仿制：可靠性很难保证。现场请当地曾亲眼看见过该沟泥石流暴发的老居民多人次，在需要测试的沟段，选取有代表性的堆积物搅拌成暴发时的泥石流流体状态，进行样品仿制，然后分别测出样品的总质量和总体积，按下式求出泥石流流体重度。c .（23-1）式中：c泥石流流体重度，tm3； 样品的总质量，t； 样品的总体积，m3。3)、现场评估在泥石

24、流沟现场，请多人次当地亲自目睹过该沟泥石流暴发或受过灾害的村民，描述泥石流流体特征和流体运动状况。然后按表23-2的特征确定泥石流流体重度。表232 泥石流流体稠度特征表 稠度特征稀浆状稠浆状稀粥状稠粥状重度c (t /m3)1.201.401.401.601.601.801.802.30表231 数量化评分（N）与重度、（1 ）关系对照表 ( h=2.65/2.70)评分重度 c(t/m3)1 1/a评分重度 c(t/m3)1 1/a评分重度 c(t/m3)1 1/a441.3001.223/1.2140.794/0.796731.5021.459/1.4190.672/0.6851021.

25、7031.765/1.7050.575/0.587451.3071.231/1.2200.7880.792741.5091.467/1.4270.6690.6811031.7101.778/1.7170.5690.584461.3141.239/1.2260.7820.788751.5161.475/1.4360.6650.6781041.7171.791/1.7290.5670.580471.3211.247/1.2330.7770.784761.5231.483/1.4440.6620.6741051.7241.804/1.7420.5650.577481.3281.256/1.2390.

26、7710.779771.5301.492/1.4530.6590.6711061.7311.817/1.7540.5620.576491.3351.264/1.2450.7670.776781.5371.500/1.4620.6560.6671071.7381.830/1.7670.5590.571501.3421.272/1.252 0.7620.772791.5441.508/1.4710.6530.6641081.7451.842/1.7800.5560.568511.3491.280/1.2580.7580.767801.5511.516/1.4800.6500.6601091.7521.855/1.7930.5530.564521.3561.288/1.2650.7530.763