泥石流动力特征计算

3.3.2泥石流的力学特征(1)容重泥石流静力学特征主要指泥石流体或浆体的容重、含水量、物质组成、流变特征、化学性质及其静力特征等。在一般地区出于泥石流的突发性、冲击力大等条件所限，难以直接测得天然泥石流容重。一般采用现场调查试验法进行泥石流容重的测定，即在现场请当地亲眼看见泥石流暴发的居民多人，在需要测试的沟段，选取有代表性的堆积物搅拌成暴发时泥石流流体状态，进行样品鉴定，然后分别测出样品的总质量和总体积，求出泥石流流体容重。在无法取得代表性样品时，根据《规范》中泥石流沟易发程度数量化评分标准，对某泥石流沟进行泥石流沟易发程度数量化评分（详见表4-3），按照《规范》中附表“数量化评分（N）与重度、（1＋Φ）关系”，可以得到泥石流的容重。本报告采用后一种方法，查表得到泥石流的容重为1.68t/m3。(2)泥石流流速泥石流的流量是泥石流重要的特征值之一。它不仅反映了泥石流的强度，规模和流体性质，而且决定着防治泥石流工程建筑物的类型、结构和尺寸。因此，泥石流的洪峰流量是泥石流研究和防治工程中不可缺少的参数。流速VC按照铁道部推荐的稀性泥石流的计算公式进行计算:5.03211cCIRnaV式中：a1——泥石流中含沙量变化引起的流速修正系数，5.0111Ha；R——水力半径(m)，2.5m；IC——泥石流水力坡度（‰），用沟床纵坡代替；n1——清水河床糙率系数；——泥石流泥沙修正系数，cHc；c——泥石流容重（t/m3），为1.68t/m3；w——清水容重（t/m3），1.0t/m3；H——泥石流中固体物质重度（t/m3）。根据以上计算公式，泥石流的平均流速为8.28m/s。(3)泥石流流量泥石流流量计算，目前主要有两种方法，一是雨洪法；二是形态调查法。①雨洪法假设泥石流与暴雨同频率、且同步发生，先按水文方法计算出断面不同频率下的小流域暴雨洪峰流量(计算方法查阅四川省水文手册)，然后选用堵塞系数，按下式计算泥石流流量：cpcDQQ1式中：cQ——频率为P的泥石流洪峰值流量（m3/s）；pQ——频率为P的暴雨洪水设计流量（m3/s）；——泥石流泥沙修正系数，查《规范》附表，值为0.71；c——泥石流容重（t/m3），为1.68t/m3；w——清水的重度（t/m3），为1.0；H——泥石流中固体物质重度；DC——泥石流堵塞系数（见表3-8），可查经验表为1.5。表3-8泥石流阵流堵塞系数DC值表堵塞程度严重堵塞中等严重堵塞轻微堵塞无堵塞Dc值2.52.5～1.51.5～1.11.0按照雨洪法，利用泥石流流量公式计算所得的泥石流最大流量如下表（表3-9）所示：表3-9雨洪法计算最大流量沟名项目某泥石流沟流域面积F(km2)35沟长L(km)11.25平均坡降J(‰)222泥沙修正系数0.71设计频率p(%)25暴雨洪峰流量QP(m3/s)87.4972.54泥石流峰值流量Qc（m3/s）224.43186.06②形态调查法在泥石流沟道中选择2-3个测流断面，仔细查找泥石流过境后留下的痕迹，然后测量这些断面上的泥石流流面比降（若不能由泥痕确定，则用沟床比降代替）、泥位高度Hc（或水力半径）和泥石流过流断面面积等参数。用相应的泥石流流速计算公式，求出平均流速Vc后，即用下式求泥石流断面峰值流量Qc：CCCQWV式中：Qc——泥石流断面峰值流量（m3/s）；Wc——泥石流过流断面面积（m2），取流通区平均沟道宽8m，水力半径为2.5m，过流断面面积为20m2；Vc——泥石流断面平均流速（m/s），为8.28m/s。用上式计算的泥石流断面峰值流量Qc为165.66m3/s。(4)一次泥石流过程总量计算一次泥石流总量Q可通过计算法和实测法确定。实测法精度高，但因往往不具备测量条件，只是一个粗略的概算。计算法根据泥石流历时T（s）和最大流量Qc（m3/s），按泥石流暴涨暴落的特点，将其过程线概化成五角形，按下式计算Q（m3），泥石流历时T取1200s：Q=0.264TQc=KTQc一次泥石流冲出的固体物质总量QH（m3）：wHwcHQQ根据该公式，在计算频率为20年一遇的情况下，某泥石流沟一次泥石流冲出的最大固体物质总量为41850m3；在计算频率为50年一遇的情况下，一次泥石流冲出的最大固体物质总量为50480m3。(5)泥石流整体冲压力泥石流冲击力是泥石流防治工程设计的重要参数。采用规范推荐的泥石流体整体冲压力计算公式：gaVccsin2式中：——泥石流整体冲压力；c——泥石流容重，1.68t/m3；cV——泥石流流速，8.28m/s；a——建筑物受力面与泥石流冲击压力方向的夹角，取90o；G——重力加速度；λ——建筑物形状系数，圆形建筑物取1.0，矩形建筑物取1.33，方形建筑物取1.47。本次建筑物主要以矩形为主，λ取1.33。通过该公式的计算，泥石流的整体冲压力约15.64×104Pa。(6)泥石流体中大石块的最大流速按照规范，用以下经验公式计算：5.0maxdaVs式中：Vs——泥石流中大石块的移动速度（m/s）；α——全面考虑摩擦系数、泥石流容重、石块比重、石块形状系数、沟床比降等因素的参数。3.5≤α≤4.5，平均α=4.0；dmax——泥石流堆积物中最大石块粒径（m），为2.5m。通过计算，泥石流大石块的最大速度为6.32m/s。(7)泥石流中石块冲击力泥石流中石块冲击力的计算参照以下公式（吴积善，《泥石流及其综合治理》1993）：12γsinαdCQPVCC式中：Pd——为泥石流中石块的冲击力（Pa）；γ——动能折减系数，对于圆端正面撞击，采用γ=0.3；Vc——泥石流平均流速（m/s），为8.28m/s；Q——石块质量（kg），按最大粒径2.5m来算；——受力面与泥石流冲击力方向的夹角，取90°；C1、C2——巨石及拦挡圬工的弹性变形系数，C1+C2=0.005。根据该公式，泥石流中石块的冲击力为14.35Kpa。(8)泥石流最大冲起高度与爬高泥石流最大冲起高度H根据规范推荐公式为：gVHc22由于泥石流在爬高过程中受到沟床阻力的影响其爬高H为：gVgbVHcc228.02式中：b——为迎面坡的函数。通过计算，泥石流的最大冲起高度为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.64×10450年一遇清水洪峰流量（雨洪法）（m3/s）87.49爬高（m）2.850年一遇泥石流洪峰流量（雨洪法）（m3/s）224.43最大冲起高度（m）3.550年一遇一次泥石流可能最大冲出量（雨洪法）（m3）50480(3)、泥石流动力学参数计算A、流速计算：据勘查所得泥石流流体水力半径、纵坡、沟床糙率及重度等参数计算；也可按泥石流的性质和所在地域，选择适合的地区性经验公式计算。泥石流流速是决定泥石流动力学性质的最重要参数之一。目前泥石流流速计算公式为半经验或经验公式，概括起来一般分为稀性泥石流流速计算公式、粘性泥石流计算公式和根据泥石流中大石块运动速度推算泥石流流速等三种办法。a、稀性泥石流流速计算公式VC=1/a·n1R32·IC21…………………………………(2—2—1)式中：VC——泥石流断面平均流速（m／s）；1/a=1/(γH·Φ＋1)1/2-----泥石流中由含沙量变化而引起的流速修正系数，查表2—3--1…R——水力半径（m），一般可用平均水深H（m）代替；IC——泥石流水力坡度（‰），一般可用沟床纵坡代替。n1——清水河床糙率系数，查当地水文手册或查铁路桥渡勘测设计规范（TBJ17--86）。表2—2--1；泥石流河道河床清水糙率表河段特征河床物质组成及形态沟岸状况1/n河段较顺直，河宽逐渐扩展，断面较规则，水流通畅，没有堵塞。砂、土质河床，河底平顺，沟床纵坡30平顺的土岸或人工堤防，沟岸山坡坡度150。≥55（50～75）砂、土质河床，河底不很平顺，沟床纵坡30～60。平顺的土岸，略有坍塌，沟岸山坡坡度15～250。50（40～55）卵、砾石河床，河底较平顺，沟床纵坡30～60。有坍塌的土岸或岩质沟岸，沟岸山坡坡度15～250。45（35～50）卵、砾石河床，河底不平顺，沟床纵坡60～120。有坍塌的土岸或岩质沟岸，坍塌较发育，沟岸山坡坡度25～320。40（30～45）卵、碎、块石河床，河底不平顺，沟床纵坡。60～120。崩、坡积物或岩质沟岸，沟岸山坡坡度25～320。坍塌发育，有堵塞痕迹。35（25～40）卵、碎、块石河床，河底松散堆积，很不平顺，沟床纵坡120。崩、坡积物或岩质沟岸，沟岸山坡坡度320。，坍塌很发育，有明显堵塞痕迹。≤25（15～35）山区峡谷河段、弯曲河段、阻塞的复式断面河段砂、园砾石河床，边滩沙州犬牙交错人工堤防强制弯曲者35（30～40）有矾石或丁坝挑流者30（25～35）卵石、园石河床，河底起不平，或长有水生植物的沟床。参差不齐的卵石、园石河岸，或长中密灌丛的河岸25（20～30）参差不齐的岩岸，或长灌木丛生的河岸20（15～25）卵石、块石、大漂石河床，石梁、跌水、孤石交错、或水生植物稠密，阻水严重的沟床0。参差不齐的岩岸，或长灌木丛生的河岸15（12～20）两岸时有岩咀突出，很不平顺，形成强烈斜流、回水、死水的河岸12（10～15）b、粘性泥石流流速计算公式21321ccccIHnVcn——粘性泥石流的河床糙率，用内插法由表2—2—2查得。表2—2—2粘性泥石流河床糙率cn序号泥石流体特征沟床状况糙率值cncn11流体呈整体运动；石块粒径大小悬殊，一般在30～50cm，2～5m粒径的石块约占20％；龙头由大石块组成，在弯道或河床展宽处易停积，后续流可超越而过，龙头流速小于龙身流速、堆积呈垄岗状河床极粗糙，沟内有巨石和挟带的树木堆积，多弯道和大跌水，沟内不能通行，人迹罕见，沟床流通段纵坡在100‰～150‰，阻力特征属高阻型平均值0.270cH2m时，0.4453.572.252流体呈整体运动，石块较大，一般石块粒径20～30cm，含少量粒径2～3m的大石块；流体搅拌较为均匀；龙头紊动强烈，有黑色烟雾及火花；龙头和龙身流速基本一致；停积后呈垄岗状堆积河床比较粗糙，凹凸不平，石块较多，有弯道、跌水；沟床流通段纵坡70‰～100‰，阻力特征属高阻型cH＜1.5m时，平均0.040cH≥1.5m时，0.050～0.100平均0.06720～302510～20153流体搅拌十分均匀；石块粒径一般在10cm左右，挟有个别2～3m的大石块；龙头和龙身物质组成差别不大；在运动过程中龙头紊动十分强烈，浪花飞溅；停积后浆体与石块不分离，向四周扩散呈叶片状沟床较稳定，河床物质较均匀，粒径10cm左右；受洪水冲刷沟底不平而且粗糙，流水沟两侧较平顺，但干而粗糙；流通段沟底纵坡55‰～70‰，阻力特征属中阻型或高阻型0.1m＜cH＜0.5m0.0430.5m＜cH＜2.0m0.0772.0m＜cH＜4.0m0.1002313104流体搅拌十分均匀；石块粒径一般在10cm左右，挟有个别2～3m的大石块；龙头和龙身物质组成差别不大；在运动过程中龙头紊动十分强烈，浪花飞溅；停积后浆体与石块不分离，向四周扩散呈叶片状泥石流铺床后原河床粘附—层泥浆体，使干而粗糙的河床变得光滑平顺，利于泥石流体运动，阻力特征属低阻型0.1＜cH＜0.5m0.0220.0＜cH＜2.0m0.0332.0＜c