泥石流动力特征计算.doc

3.3.2泥石流的力学特征(1)容重泥石流静力学特征主要指泥石流体或浆体的容重、含水量、物质组成、流变特征、化学性质及其静力特征等。在一般地区出于泥石流的突发性、冲击力大等条件所限，难以直接测得天然泥石流容重。一般采用现场调查试验法进行泥石流容重的测定，即在现场请当地亲眼看见泥石流暴发的居民多人，在需要测试的沟段，选取有代表性的堆积物搅拌成暴发时泥石流流体状态，进行样品鉴定，然后分别测出样品的总质量和总体积，求出泥石流流体容重。在无法取得代表性样品时，根据规范中泥石流沟易发程度数量化评分标准，对某泥石流沟进行泥石流沟易发程度数量化评分（详见表4-3），按照规范中附表“数量化评分（N）与重度、（1）关系”，可以得到泥石流的容重。本报告采用后一种方法，查表得到泥石流的容重为1.68t/m3。(2)泥石流流速泥石流的流量是泥石流重要的特征值之一。它不仅反映了泥石流的强度，规模和流体性质，而且决定着防治泥石流工程建筑物的类型、结构和尺寸。因此，泥石流的洪峰流量是泥石流研究和防治工程中不可缺少的参数。流速VC 按照铁道部推荐的稀性泥石流的计算公式进行计算:式中： 泥石流中含沙量变化引起的流速修正系数，；R水力半径(m)，2.5m；IC泥石流水力坡度（），用沟床纵坡代替；清水河床糙率系数；泥石流泥沙修正系数，；泥石流容重（t/m3），为1.68t/m3；清水容重（t/m3），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 泥石流堵塞系数（见表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个测流断面，仔细查找泥石流过境后留下的痕迹，然后测量这些断面上的泥石流流面比降（若不能由泥痕确定，则用沟床比降代替）、泥位高度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可通过计算法和实测法确定。实测法精度高，但因往往不具备测量条件，只是一个粗略的概算。计算法根据泥石流历时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；建筑物受力面与泥石流冲击压力方向的夹角，取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)泥石流中石块冲击力泥石流中石块冲击力的计算参照以下公式（吴积善，泥石流及其综合治理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。根据以上的计算，将泥石流各动力参数的计算结果汇总于表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年一遇泥石流洪峰流量（雨洪法）（m3/s）224.43最大冲起高度（m）3.550年一遇一次泥石流可能最大冲出量（雨洪法）（m3）50480(3)、泥石流动力学参数计算 A、流速计算：据勘查所得泥石流流体水力半径、纵坡、沟床糙率及重度等参数计算；也可按泥石流的性质和所在地域，选择适合的地区性经验公式计算。泥石流流速是决定泥石流动力学性质的最重要参数之一。目前泥石流流速计算公式为半经验或经验公式，概括起来一般分为稀性泥石流流速计算公式、粘性泥石流计算公式和根据泥石流中大石块运动速度推算泥石流流速等三种办法。a、 稀性泥石流流速计算公式VC =1/aRIC (221)式中： VC 泥石流断面平均流速（ms）；1/a =1 / (H 1)1/2 -泥石流中由含沙量变化而引起的流速修正系数，查表23-1R 水力半径（m），一般可用平均水深H（m）代替；IC 泥石流水力坡度（），一般可用沟床纵坡代替。 清水河床糙率系数，查当地水文手册或查铁路桥渡勘测设计规范（TBJ17-86）。表22-1； 泥石流河道河床清水糙率表河段特征河床物质组成及形态沟岸状况1 / n河段较顺直，河宽逐渐扩展，断面较规则，水流通畅，没有堵塞。砂、土质河床，河底平顺，沟床纵坡 30平顺的土岸或人工堤防，沟岸山坡坡度 120。崩、坡积物或岩质沟岸，沟岸山坡坡度 320。，坍塌很发育，有明显堵塞痕迹。25（1535）山区峡谷河段、弯曲河段、阻塞的复式断面河段砂、园砾石河床，边滩沙州犬牙交错人工堤防强制弯曲者35（3040）有矾石或丁坝挑流者30（2535）卵石、园石河床，河底起不平，或长有水生植物的沟床。参差不齐的卵石、园石河岸，或长中密灌丛的河岸25（2030）参差不齐的岩岸，或长灌木丛生的河岸20（1525）卵石、块石、大漂石河床，石梁、跌水、孤石交错、或水生植物稠密，阻水严重的沟床0。参差不齐的岩岸，或长灌木丛生的河岸15（1220）两岸时有岩咀突出，很不平顺，形成强烈斜流、回水、死水的河岸12（1015）b、粘性泥石流流速计算公式 粘性泥石流的河床糙率，用内插法由表222查得。表222 粘性泥石流河床糙率序号泥石流体特征沟床状况糙率值1流体呈整体运动；石块粒径大小悬殊，一般在3050cm，25m粒径的石块约占20；龙头由大石块组成，在弯道或河床展宽处易停积，后续流可超越而过，龙头流速小于龙身流速、堆积呈垄岗状河床极粗糙，沟内有巨石和挟带的树木堆积，多弯道和大跌水，沟内不能通行，人迹罕见，沟床流通段纵坡在100150，阻力特征属高阻型平均值0.2702m时，0.4453.572.252流体呈整体运动，石块较大，一般石块粒径2030cm，含少量粒径23m的大石块；流体搅拌较为均匀；龙头紊动强烈，有黑色烟雾及火花；龙头和龙身流速基本一致；停积后呈垄岗状堆积河床比较粗糙，凹凸不平，石块较多，有弯道、跌水；沟床流通段纵坡70100，阻力特征属高阻型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的大石块；龙头和龙身物质组成差别不大；在运动过程中龙头紊动十分强烈，浪花飞溅；停积后浆体与石块不分离，向四周扩散呈叶片状泥石流铺床后原河床粘附层泥浆体，使干而粗糙的河床变得光滑平顺，利于泥石流体运动，阻力特征属低阻型0.10.5m0.0220.02.0m0.0332.04.0m0.050462620c、 泥石流中石块运动速度推算泥石流流速计算公式在缺乏大量实验数据和实测数据阶情况下，为便于以堆积后的泥石流冲出物最大粒径大体推求石块运动速度推算泥石流流速的经验公式： Vskd max.（222）式中：泥石流中大石块的移动速度推算泥石流流速（ms）；泥石流堆积物中最大石块的粒径（m）；k全面考虑的摩擦系数（泥石流重度、石块密度、石块形状系数、沟床比降等因素）。变化范围3.55.5，B、 流量计算：泥石流流量可采用形态调查法（据泥痕勘测所得的过流断面面积乘以流速）或雨洪法（按暴雨洪水流量乘以泥石流修正系数）确定。暴雨小径流的地区性经验公式较多，暴雨洪水流量应采用适合当地的经验公式计算。a、 形态调查法 在泥石流沟道中选择23个有代表性的过流断面。查找泥石流过境后留下的痕迹，然后确定泥位。最后测量这些断面上的泥石流流面比降（若不能由痕迹确定，则用沟床比降代替）、泥位高度HC（或水力半径）和泥石流过流断面面积等参数。用相应的泥石流流速计算公式，求出断面平均流速VC后，即可用下式求泥石流断面峰值流量QC。QCWCVC .（223）式中：WC 泥石流过流断面面积（m2）；VC 泥石流断面平均流速（ms）。 雨洪法假设泥石流与暴雨同频率、且同步发生，计算断面的暴雨洪水设计流量全部转变成泥石流流量的前题下建立的计算方法。其计算步骤是先按水文方法计算出断面不同频率下的小流域暴雨洪峰流量 (计算方法查阅水文手册)，然后选用泥石流重度和堵塞系数，按式I-2计算泥石流流量。 QC =（1十）QPDC.（224）式中：QC 频率为P的泥石流洪峰值流量（m3s）；QP 频率为P的暴雨洪水设计流量（m3s）； 根据各省水文手册中给出的计算公式计算，或用下式 QB=0.278 riF 计算式中：r按小时平均雨强（毫米 / 小时）设计，用实测最大小时雨强校核 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； 表223 泥石流堵塞系数Dc值堵塞程度特 征堵塞系数Dc严重 河槽弯曲，河段宽窄不均，卡口、陡坎多。大部分支沟交汇角度大，形成区集中。物质组成粘性大，稠度高，沟槽堵塞严重，阵流间隔时间长2.5 中等 沟槽较顺直，沟段宽窄较均匀，陡坎、卡口不多。主支沟交角多小于60，形成区不太集中。河床堵塞情况一般，流体多呈稠浆稀粥状1.52.5轻微 沟槽顺直均匀，主支沟交汇角小，基本无卡口、陡坎，形成区分散。物质组成粘度小，阵流的间隔时间短而少1.11.5无1.0一次泥石流过程总量计算一次泥石流总量Q可通过计算法和实测法确定。实测法精度高，但因往往不具备测量条件，只是一个粗略的概算。计算法根据泥石流历时T（s）和最大流量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），取g9.8ms2；建筑物受力面与泥石流冲击力方向的夹角（）；建筑物形状系数： 圆形建筑物1.0，矩形建筑物1.33，方形建筑物：1.47。D、 弯道超高与冲高计算：泥石流流动在弯曲沟道外侧产生的超高值和泥石流正面遇阻的冲起高度。a)、泥石流的弯道超高由于泥石流流速快，惯性大，故在弯道凹岸处有比水流更加显著的弯道超高现象。根据弯道泥面横比降动力平衡条件，推导出计算弯道超高的公式：.（229）式中：h 弯道超高（m）；R2 凹岸曲率半径（m）；R1 凸岸曲率半径（m）；b)、泥石流最大冲起高度 . （2210）c）泥石流在爬高过程中由于受到沟床阻力的影响，其爬高H：.（2211）（4）泥石流的形成区、流通区和堆积区测绘。 a）工程治理区实测剖面至少按一纵三横控制； b）重点区应有1-3个探槽或探坑（井）控制；5)、对各类防治工程提供以下主要设计参数。（1）各类拦挡坝: 覆盖层和基岩的重度、承载力标准值、抗剪强度，基面摩擦系数，泥石流的性质与类型，发生频次，泥石流体的重度和物质组成，泥石流体的流速、流量和设计暴雨洪水频率，泥石流回淤坡度和固体物质颗粒成分，沟床清水冲刷线。（2）其它工程：桩林着重于其锚固段基岩深度、风化程度和力学性质；排导槽、渡槽着重于泥石流运动的最小坡度、冲击力、弯道超高和冲高；导流堤、护岸堤和防冲墩着重于基岩的埋藏深度和性质、泥石流冲击力和弯道超高、墙背摩擦角；停淤场着重于淤积总量、淤积总高度和分期淤积高度。2泥石流流体重度(c)的确定泥石流流体重度(c)的确定，目前测试手段还不能解决，因此暂时的替代办法有：1)、根据数量化评分法确定重度：根据数量化评分（N）建立与重度(c)的线性关系，制成备查表直接查找，参见表23-12)、现场仿制：可靠性很难保证。现场请当地曾亲眼看见过该沟泥石流暴发的老居民多人次，在需要测试的沟段，选取有代表性的堆积物搅拌成暴发时的泥石流流体状态，进行样品仿制，然后分别测出样品的总质量和总体积，按下式求出泥石流流体重度。c .（23-1）式中：c泥石流流体重度，tm3； 样品的总质量，t； 样品的总体积，m3。3)、现场评估在泥石流沟现场，请多人次当地亲自目睹过该沟泥石流暴发或受过灾害的村民，描述泥石流流体特征和流体运动状况。然后按表23-2的特征确定泥石流流体重度。表232 泥石流流体稠度特征表 稠度特征稀浆状稠浆状稀粥状稠粥状重度c (t /m3)1.201.401.401.601.601.801.802.30表231 数量化评分（N）与重度、（1j）关系对照表 (gh=2.65/2.70)评分重度gc(t/m3)1j1/a评分重度gc(t/m3)1j1/a评分重度gc(t/m3)1j1/a441.3001.223/1.2140.794/0.796731.5021.459/1.4190.672/0.6851021.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.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.763811.5581.524/1.4890.6470.6571101.7591.868/1.8070.5500.561531.3631.296/1.2720.7490.759821.5651.532/1.4980.6440.6531111.7661.881/1.8200.5480.558541.3701.304/1.2780.7440.756831.5721.540/1.5070.6410.6501121.7721.894/1.8320.5450.555551.3771.313/1.2850.7390.752841.5791.549/1.5170.6380.6461131.7791.907/1.8460.5420.552561.3841.321/1.2920.7350.748851.5861.557/1.5140.6360.6471141.7861.919/1.8600.5390.549571.3911.329/1.2990.7310.744861.5931.565/1.5360.6330.6391151.7931.932/1.8740.5370.546581.3981.337/1.3060.7270.740871.6001.577/1.5450.6290.6361161.8001.945/1.8890.5340.542591.4051.345/1.3120.7230.737881.6071.586/1.5550.6260.6331171.8432.208/1.9840.4880.523601.4121.353/1.3200.7190.733891.6141.599/1.5650.6220.6291181.8862.471/2.0880.4520.504611.4191.361/1.3270.7150.729901.6211.611/1.5760.6180.6261191.9292.735/2.2040.4230.485621.4261.370/1.3340.7110.725911.6281.624/1.5860.6140.6221201.9712.998/2.3320.3990.466631.4331.378/1.3420.7070.721921.6341.637/1.5950.6100.6201212.0143.216/2.4780.3810.448641.4401.386/1.3490.7040.717931.6411.650/1.6050.6060.6161222.0573.524/2.6440.3610.429651.4471.394/1.3570.6990.714941.6481.663/1.6160.6020.6131232.1003.788/2.8330.3450.410661.4531.402/1.3630.6960.711951.6551.676/1.6270.5990.6091242.1434.051/3.0520.3320.391671.4601.410/1.3710.6920.707961.6621.688/1.6380.5950.6061252.1864.314/3.3070.3200.372681.4671.418/1.3790.6890.703971.6691.701/1.6490.5920.6031262.2294.577/3.6090.3090.353691.4741.426/1.3870.6850.699981.6761.714/1.6600.5880.5991272.2714.840/3.9630.2990.333701.4811.435/1.3950.6820.695991.6831.727/1.6720.5860.5961282.3145.104/4.4040.2900.313711.4881.443/1.4030.6780.6921001.6901.740/1.6830.5810.5931292.3575.367/4.9560.2820.293721.4951.451/1.4110.6750.6881011.6971.753/1.6950.5780.5901302.4005.630/5.6670.2750.271注：j-为泥石流中含沙量的修正好系数。3、泥沙级配取样分析1)、体积法在需要试验的沟段，选择有代表性的试验点，取样坑的边长取23倍最大粒径、深取12个最大粒径，清除表层杂质，取出其全部土、砂、石，从中挑出粒径大于200mm以上的石块单个分别称重，其余按粒径分别按200150mm，150100mm，10050mm，5020mm，20mm以下分成若干级，每级分组称重，计算分组粒径，并绘制颗粒级配曲线。按粒径与重量关系的理论公式：D(cm)8.91 3G(kg) (232)求算泥砂粒径(也可以直接查表)，然后求得所需特征值。2)、方格网法在取样地段，选出代表性沟段画出100个1m1m的小方格，取每个小方格交点上的一石块（剔除个别大孤石）来作统计。量取每个石块的三轴向尺寸（长、宽、高），计算三边尺寸的几何均值或算术平均值，作为该石块的平均直径。然后按粒径大小分成若干个粒径组，称出各粒径组的质量与总质量之比，绘制颗粒级配曲线，求算颗粒级配曲线并求得所需特征值。取泥石流浆体，使用标准粘度计或旋转粘度计和泥浆静切力计测试。1.形态调查法 形态调查法指根据沟床内以往发生过的泥石流痕迹，测量泥位高和过流断面面积，计算平均流速和流量，再根据其发生的日期，确定其经验频率，推算设计流量的方法。这种方法在无或少资料区经常采用，在有资料区也是作为相互效核验证的主要方法。其流量计算式和经验频率计算式如下: Qc=Uc式中: Qc为调查断面的泥石流流量(米3/秒); Uc为通过查断面的泥石流平均流速(米/秒);。为调查断面的过程面积(米2)。 P=ni/(n+1)*100% (这个式子在这里起什么作用？)式中:P为本次泥石流出现的经验频率(%) ;n为调查期的总年数;ni为n年中按泥石流流量大小排列所得到的本次泥石流的序号。在有比较清晰泥痕和比较好断面的情况下，这种方法的计算结果是最可靠的。 2.配方法 配方法是在清水流量中配上土体含量，使其成为泥石流流量的方法。但需假定土体是在水体产流和汇流过程中通过冲刷同步加入，泥石流的频率和洪水的频率一致，洪水全部变成泥石流。这种方法的计算步骤是，先按水文方法计算出断面的不同频率的小流域暴雨洪峰流量，再根据不同情况计算泥石流流量。 (1)不考虑泥石流土体的天然含水量，其计算公式如下: QC = (1+c)Qwc=（c-w）（s-c）式中:Qw为某一频率的暴雨洪峰设计流量(米3/秒)。Qc为Qw同频率的泥石流流量(米3/秒);为泥石流流量增加系数;rc为泥石流体容重(克/厘米3;吨/米3);Yw为清水容重(克/厘米3;吨/米3);Ys为土体实体容重(克/厘米3;吨/米3(2)考虑泥石流土体的天然含水量，其计算公式为: QC = (1+c)Qwc=（c-1）/s(1+Pw)- c(1+s Pw)式中: c为考虑泥石流土体的天然含水量的流量增加系数，Pw为泥石流土体的天然含水量(%):其它符号同前。(3)考虑堵塞，其计算公式为: QC = (1+c)QwDw QC = (1+c)QwDw式中:D，为泥石流堵塞系数，其它符号同前。 3，地区性经验公式 前人根据对某一地区泥石流流量观测资料的统计分析，建立了一些地区经验性公式，这些公式大多也是在配方法的基础上建立起来的，适合于该地区同类型泥石流沟的泥石流流量计算，也可作为类似地区的参考。 4.实测法 在经常暴发泥石流的沟谷，对泥石流流量进行实地观测，同时进行雨量观测，利用流量与降雨的关系，并以降雨的频率作为这次泥石流出现的频率进行统计分析，来求算设计泥石流流量。 但是，这些间接计算方法都不能算是可靠的，也不能推广到各种类型的泥石流和各个泥石流形成流域中，因此，在选用这些公式时，应当考虑不同地区、不同类型、不同性质的泥石流特点，其计算结果与现场调查和历史资料进行比较，然后确定应选用的数据。泥石流及其综合治理1. 流量 泥石流峰值流量是泥石流重要的特征值之一它不仅反映了泥石流的强度、规程和流体的性质，而且决定着防治泥石流工程建筑物的类型、结构和尺寸因此，泥石流峰值流量是泥石流研究和防治工程中不可缺少的参数，故引起国内外泥石流专家所关注本世纪50年代我国大量采用苏联学者的计算方法以适应公路和铁路建设的迫切需要当时一般采用配方法和形态调查法进行峰值流量的估算但这些方法之间误差很大，通过修改参比改进了一些，但其计算结果偏差仍顾大 确定泥石流最大流量Qc 最可靠的方法是直接测量自然界里奔流着的泥石流或泥石流遗留的痕迹可是这样确定泥石流最大流量时，首先，要确定泥石流最高泥位时的断面面积和相应的平均流速断面与流速的测定又受泥石流频率的影响，要取得其参数比较困难配方法：此法根据泥石流体中水和固体物质的比例，用在一定设计标准下可能出现的洪水流量加上按比例所需的团体物质体积配合而成的泥石流流量此法的基木表达式为：Qc=Qb（1+）.1Qc为泥石流流量Qb为清水设计流量为泥石流修正系数，=（C-B）/（H-C）C为泥石流容重B为水的容重H泥石流中泥沙比重配方法是目前泥石流流量计算的基本方法采用配方法进行泥石流流量计算，对于稀性泥石流来说，与实际结果差别不大，对于粘性泥石流来说，特别是塑性泥石流，推求值远小于实测值因此国内外许多学者都在配方法的基础上，再乘上一个系数，来表示某种原因而引起的站性泥石流流量的增值 (1)考虑补给沙石体含水量的配方法，QcQb (1十P)2式中，P为考虑补给固体物质含水量的泥石流流量修正系致P=(1-CVB)/ CVB -(1- CVB)CVB=（H-C）/（H-1）为泥石流补给区中固体物质的体积比原始水量或天然含水量，可以实测也可以采用当地水文气象站资料。该法是假定在给定的设计容重下，水在形成泥石流的过程中起决定性的作用因此，其缺点未能考虑到固体物质的补给规律以及泥石流形成过程中水和固体物质遭遇的情况，所以计算结果偏小但此法比式(1)要精确，因为式(2)不仅考虑了固体的补给量，而旦还考虑了补给沙石体中的水的含量 (2）考虑堵塞条件下的配力法 泥石流在通过卡口、急弯、纵坡突然变缓的情况下的沟段，常常发生泥石流停积堵塞、加积增大而又开始流动，成为泥石流流量增大的重要原因一殷表达式为：Qc=QB（1+）q式中，q为泥石流堵塞系数(其值为l一3.0），其中微弱堵塞11.4，一殷堵塞15一19，较严重堵塞20一25，严重堵塞263.0）据云南东川大白泥沟19601952年观测资料的分析，求得帖性泥石流的堵塞修正系数q与堵塞时间t的经验关系qO87t0.24 式中，t为堵塞时间，即阵性泥石流之间的断流时间根据上述资队求得泥石流堵塞系数q与产生泥石流流量Qc的经验关系为： q58Qc0.21 适用于粘性阵流关于泥石流堵塞现象的成因很多它与泥石流浆体的性质、组成、规模和固体物质来源的方式及数量等因素有关由于泥石流堵塞后又溃决，使泥石流流量值增大，其变化的大小随泥石流堵塞的时间长短而定，其堵塞程度随泥石流流量大小而变泥石流流量与堵塞系数之间的关系反映了泥石流流体堵塞加积，最后溃决的基本规律：即泥石流流量小、堵塞加积大时才可能达到溃决的流量，流量大、堵塞加积小时，即堵塞系数小，泥石流体溃决便开始运动，当流量达到某种程度时，不需要加积增大，就可能在河床中顺畅通过2形态调查法(又称泥鹰调查法) 从形态调查求算泥石流峰值流里的方法，系按照泥石流在河槽岸边上遗留的最高痕迹，测该处河床横断面积，乘上选用适合的泥石