泥石流动力特征计算

3.3.2泥石流的力学特征(1)容重泥石流静力学特征主要指泥石流体或浆体的容重、含水量、物质组成、流变特征、化学性质及其静力特征等。

在一般地区出于泥石流的突发性、冲击力大等条件所限，难以直接测得天然泥石流容重。

一般采用现场调查试验法进行泥石流容重的测定，即在现场请当地亲眼看见泥石流暴发的居民多人，在需要测试的沟段，选取有代表性的堆积物搅拌成暴发时泥石流流体状态，进行样品鉴定，然后分别测出样品的总质量和总体积，求出泥石流流体容重。

在无法取得代表性样品时，根据《规范》中泥石流沟易发程度数量化评分标准，对某泥石流沟进行泥石流沟易发程度数量化评分（详见表4-3），按照《规范》中附表“数量化评分（N ）与重度、（1＋Φ）关系”，可以得到泥石流的容重。

本报告采用后一种方法，查表得到泥石流的容重为1.68t/m3。

(2)泥石流流速泥石流的流量是泥石流重要的特征值之一。

它不仅反映了泥石流的强度，规模和流体性质，而且决定着防治泥石流工程建筑物的类型、结构和尺寸。

因此，泥石流的洪峰流量是泥石流研究和防治工程中不可缺少的参数。

流速VC 按照铁道部推荐的稀性泥石流的计算公式进行计算: 式中：a 1——泥石流中含沙量变化引起的流速修正系数，()5.0111+Φ=H a γ；R ——水力半径(m)，2.5m ；IC ——泥石流水力坡度（‰），用沟床纵坡代替；n 1——清水河床糙率系数；Φ——泥石流泥沙修正系数，()()cHcγγγγω--=Φ；cγ——泥石流容重（t/m3），为1.68t/m3；wγ——清水容重（t/m3），1.0t/m3；Hγ——泥石流中固体物质重度（t/m3）。

根据以上计算公式，泥石流的平均流速为8.28m/s。

(3)泥石流流量泥石流流量计算，目前主要有两种方法，一是雨洪法；二是形态调查法。

①雨洪法假设泥石流与暴雨同频率、且同步发生，先按水文方法计算出断面不同频率下的小流域暴雨洪峰流量(计算方法查阅四川省水文手册)，然后选用堵塞系数，按下式计算泥石流流量：式中：c Q——频率为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值表按照雨洪法，利用泥石流流量公式计算所得的泥石流最大流量如下表（表3-9）所示：表3-9雨洪法计算最大流量②形态调查法在泥石流沟道中选择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.66m3/s。

(4)一次泥石流过程总量计算一次泥石流总量Q可通过计算法和实测法确定。

实测法精度高，但因往往不具备测量条件，只是一个粗略的概算。

计算法根据泥石流历时T（s）和最大流量Qc（m3/s），按泥石流暴涨暴落的特点，将其过程线概化成五角形，按下式计算Q（m3），泥石流历时T取1200s：Q=0.264TQc=KTQc一次泥石流冲出的固体物质总量QH（m3）：根据该公式，在计算频率为20年一遇的情况下，某泥石流沟一次泥石流冲出的最大固体物质总量为41850m3；在计算频率为50年一遇的情况下，一次泥石流冲出的最大固体物质总量为50480m3。

(5)泥石流整体冲压力泥石流冲击力是泥石流防治工程设计的重要参数。

采用规范推荐的泥石流体整体冲压力计算公式：式中：δ——泥石流整体冲压力；γ——泥石流容重，1.68t/m3；cV——泥石流流速，8.28m/s；ca——建筑物受力面与泥石流冲击压力方向的夹角，取90o；G——重力加速度；λ——建筑物形状系数，圆形建筑物取1.0，矩形建筑物取1.33，方形建筑物取1.47。

本次建筑物主要以矩形为主，λ取1.33。

通过该公式的计算，泥石流的整体冲压力约15.64×104Pa。

(6)泥石流体中大石块的最大流速按照规范，用以下经验公式计算：式中：Vs——泥石流中大石块的移动速度（m/s）；α——全面考虑摩擦系数、泥石流容重、石块比重、石块形状系数、沟床比降等因素的参数。

3.5≤α≤4.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)泥石流最大冲起高度与爬高∆根据规范推荐公式为：泥石流最大冲起高度H∆为：由于泥石流在爬高过程中受到沟床阻力的影响其爬高H式中：b——为迎面坡的函数。

通过计算，泥石流的最大冲起高度为3.5m，爬高为2.8m。

根据以上的计算，将泥石流各动力参数的计算结果汇总于表3-11。

表3-11泥石流动力学参数计算结果(3)、泥石流动力学参数计算A 、流速计算：据勘查所得泥石流流体水力半径、纵坡、沟床糙率及重度等参数计算；也可按泥石流的性质和所在地域，选择适合的地区性经验公式计算。

泥石流流速是决定泥石流动力学性质的最重要参数之一。

目前泥石流流速计算公式为半经验或经验公式，概括起来一般分为稀性泥石流流速计算公式、粘性泥石流计算公式和根据泥石流中大石块运动速度推算泥石流流速等三种办法。

a 、稀性泥石流流速计算公式VC=1/a ·n 1R 32·IC 21…………………………………(2—2—1) 式中：VC ——泥石流断面平均流速（m ／s ）；1/a=1/(γH ·Φ＋1)1/2-----泥石流中由含沙量变化而引起的流速修正系数，查表2—3--1…R ——水力半径（m ），一般可用平均水深H （m ）代替；IC ——泥石流水力坡度（‰），一般可用沟床纵坡代替。

n 1——清水河床糙率系数，查当地水文手册或查铁路桥渡勘测设计规范（TBJ17--86）。

表2—2--1；泥石流河道河床清水糙率表b、粘性泥石流流速计算公式n——粘性泥石流的河床糙率，用内插法由表2—2—2查得。

c表2—2—2粘性泥石流河床糙率c nc、泥石流中石块运动速度推算泥石流流速计算公式在缺乏大量实验数据和实测数据阶情况下，为便于以堆积后的泥石流冲出物最大粒径大体推求石块运动速度推算泥石流流速的经验公式：Vs＝k·√dmax……………………………….…………..（2—2—2）式中：V——泥石流中大石块的移动速度推算泥石流流速（m／s）；sd——泥石流堆积物中最大石块的粒径（m）；m axk——全面考虑的摩擦系数（泥石流重度、石块密度、石块形状系数、沟床比降等因素）。

变化范围3.5～5.5，B、流量计算：泥石流流量可采用形态调查法（据泥痕勘测所得的过流断面面积乘以流速）或雨洪法（按暴雨洪水流量乘以泥石流修正系数）确定。

暴雨小径流的地区性经验公式较多，暴雨洪水流量应采用适合当地的经验公式计算。

a、形态调查法在泥石流沟道中选择2—3个有代表性的过流断面。

查找泥石流过境后留下的痕迹，然后确定泥位。

最后测量这些断面上的泥石流流面比降（若不能由痕迹确定，则用沟床比降代替）、泥位高度HC（或水力半径）和泥石流过流断面面积等参数。

用相应的泥石流流速计算公式，求出断面平均流速VC后，即可用下式求泥石流断面峰值流量QC。

QC＝WC·VC…………………………………………………….（2—2—3）式中：WC——泥石流过流断面面积（m2）；VC——泥石流断面平均流速（m／s）。

雨洪法假设泥石流与暴雨同频率、且同步发生，计算断面的暴雨洪水设计流量全部转变成泥石流流量的前题下建立的计算方法。

其计算步骤是先按水文方法计算出断面不同频率下的小流域暴雨洪峰流量(计算方法查阅水文手册)，然后选用泥石流重度和堵塞系数，按式I-2计算泥石流流量。

QC=（1十Φ）QP·DC……………………………………………..（2—2—4）式中：QC——频率为P的泥石流洪峰值流量（m3／s）；QP——频率为P的暴雨洪水设计流量（m3／s）；根据各省水文手册中给出的计算公式计算，或用下式QB=0.278r×i×F计算式中：r——按小时平均雨强（毫米/小时）设计，用实测最大小时雨强校核iB——产水系数。

植被具有保水功能，降雨后入渗少，产水系数与雨强的大小和植被的多少呈正变，与松散土层的厚度呈反变，一般产水系数i=0.3——0.8；泥百流地区多在0.6以下为宜。

产沙系数则相反，它与植被的多少呈反变，与雨强的大小和松散土层的厚度呈正变，一般产沙系数is=0.3——0.8；F——流域面积（平方公里）（1十Φ）——泥石流中由含沙量变化而引起的流量修正系数，查表2—3--1。

式中：Φ＝（γC-γW）/（γH–γc）…………………………………（2—2—5）γW——清水的重度（t／m3）；γH——泥石流中固体物质比重（t／m3）；DC——泥石流堵塞系数，可查经验表2—2—3；表2—2—3泥石流堵塞系数Dc值一次泥石流过程总量计算一次泥石流总量Q可通过计算法和实测法确定。

实测法精度高，但因往往不具备测量条件，只是一个粗略的概算。

计算法根据泥石流历时T（s）和最大流量Qc（m3/s），按泥石流暴涨暴落的特点，将其过程线概化成五角形，按式I.6计Q=0.264TQc=KTQc………………………………（2—2—6）式中K值的变化随流域面积(F)的大小而变化：当F<5(km2)时：K=0.202；F=10～100(km2)时：K=0.0378； 泥石流体中的固体物质流量QH （m3）：QH ＝Q （γC -γW）/（γH–γW）…………………………………(2—2—7)C 、冲击力计算：可用附录I 中公式计算泥石流整体冲击力、泥石流中大石块冲击力。