泥石流动力特征计算 (1)

泥石流的力学特征

(1)容重

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

(2)泥石流流速

泥石流的流量是泥石流重要的特征值之一。它不仅反映了泥石流的强度，规模和流体性质，而且决定着防治泥石流工程建筑物的类型、结构和尺寸。因此，泥石流的洪峰流量是泥石流研究和防治工程中不可缺少的参数。

流速VC 按照铁道部推荐的稀性泥石流的计算公式进行计算:

式中： a 1——泥石流中含沙量变化引起的流速修正系数，()5

.0111+Φ=H a γ；

R ——水力半径(m)，2.5m ；

IC ——泥石流水力坡度（‰），用沟床纵坡代替；

n ——清水河床糙率系数；

Φ——泥石流泥沙修正系数，

()()c H c γγγγω--=

Φ；

c γ——泥石流容重（t/m3）

，为m3；

w γ——清水容重（t/m3）

，m3；

H γ——泥石流中固体物质重度（t/m3）。 根据以上计算公式，泥石流的平均流速为s 。

(3)泥石流流量

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

①雨洪法

假设泥石流与暴雨同频率、且同步发生，先按水文方法计算出断面不同频率下的小流域暴雨洪峰流量 (计算方法查阅四川省水文手册)，然后选用堵塞系数，按下式计算泥石流流量：

式中：

c

Q ——频率为P 的泥石流洪峰值流量（m3/s ）；

p

Q ——频率为P 的暴雨洪水设计流量（m3/s ）；

Φ——泥石流泥沙修正系数，查《规范》附表，值为； c

γ——泥石流容重（t/m3），为m3；

w——清水的重度（t/m3），为；

γ——泥石流中固体物质重度；

H

DC ——泥石流堵塞系数（见表3-8），可查经验表为。

表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为 m3/s。

(4)一次泥石流过程总量计算

一次泥石流总量Q可通过计算法和实测法确定。实测法精度高，但因往往不具备测量条件，只是一个粗略的概算。计算法根据泥石流历时T（s）和最大流量Qc（m3/s），按泥石流暴涨暴落的特点，将其过程线概化成五角形，按下式计算Q（m3），泥石流历时T取1200 s：Q = =KTQc

一次泥石流冲出的固体物质总量QH（m3）：

根据该公式，在计算频率为20年一遇的情况下，某泥石流沟一次泥石流冲出的最大固体物质总量为41850m3；在计算频率为50年一遇的情况下，一次泥石流冲出的最大固体物质总量为50480m3。

(5)泥石流整体冲压力

泥石流冲击力是泥石流防治工程设计的重要参数。

采用规范推荐的泥石流体整体冲压力计算公式：

式中：δ——泥石流整体冲压力；

γ——泥石流容重，m3；

c

V——泥石流流速，8.28m/s；

c

a——建筑物受力面与泥石流冲击压力方向的夹角，取90o；

G——重力加速度；

λ——建筑物形状系数，圆形建筑物取，矩形建筑物取，方形建筑物取。本次建筑物主要以矩形为主，λ取。

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

(6)泥石流体中大石块的最大流速

按照规范，用以下经验公式计算：

式中： Vs——泥石流中大石块的移动速度（m/s）；

α——全面考虑摩擦系数、泥石流容重、石块比重、石块形状系数、沟床比降等因素的参数。≤α≤，平均α=；

dmax——泥石流堆积物中最大石块粒径（m），为2.5m。

通过计算，泥石流大石块的最大速度为s。

(7)泥石流中石块冲击力

泥石流中石块冲击力的计算参照以下公式（吴积善，《泥石流及其综合治理》1993）：式中： Pd——为泥石流中石块的冲击力（Pa）；

γ——动能折减系数，对于圆端正面撞击，采用γ=；

Vc——泥石流平均流速（m/s），为8.28m/s；

Q——石块质量（kg），按最大粒径2.5m来算；

α——受力面与泥石流冲击力方向的夹角，取90°；

C1、C2——巨石及拦挡圬工的弹性变形系数，C1+C2=。

根据该公式，泥石流中石块的冲击力为。

(8)泥石流最大冲起高度与爬高

?根据规范推荐公式为：

泥石流最大冲起高度H

?为：

由于泥石流在爬高过程中受到沟床阻力的影响其爬高H

式中：b——为迎面坡的函数。

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

根据以上的计算，将泥石流各动力参数的计算结果汇总于表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 、粘性泥石流流速计算公式

c

n ——粘性泥石流的河床糙率，用内插法由表2—2—2查得。

表2—2—2 粘性泥石流河床糙率

c

n

c、泥石流中石块运动速度推算泥石流流速计算公式

在缺乏大量实验数据和实测数据阶情况下，为便于以堆积后的泥石流冲出物最大粒径大体推求石块运动速度推算泥石流流速的经验公式：

Vs＝k·√d max……………………………….…………..（2—2—2）

式中：

V——泥石流中大石块的移动速度推算泥石流流速（m／s）；

s

d——泥石流堆积物中最大石块的粒径（m）；

m ax

k——全面考虑的摩擦系数（泥石流重度、石块密度、石块形状系数、沟床比降等因素）。

变化范围～，

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= r×i×F 计算

式中：r——按小时平均雨强（毫米 / 小时）设计，用实测最大小时雨强校核 iB——产水系数。植被具有保水功能，降雨后入渗少，产水系数与雨强的大小和植被的多少呈正变，与松散土层的厚度呈反变，一般产水系数i=——；泥百流地区多在以下为宜。产沙系数则相反，它与植被的多少呈反变，与雨强的大小和松散土层的厚度呈正变，一般产沙系数is=——；

F——流域面积（平方公里）

（1十Φ）——泥石流中由含沙量变化而引起的流量修正系数，查表2—3--1。

式中：Φ＝（γC -γW）/（γH –γc）…………………………………（2—2—5）

γC ——泥石流重度（t／m3）；

γW ——清水的重度（t／m3）；

γH ——泥石流中固体物质比重（t／m3）；

DC ——泥石流堵塞系数，可查经验表2—2—3；

表2—2—3 泥石流堵塞系数Dc值

一次泥石流过程总量计算

一次泥石流总量Q可通过计算法和实测法确定。实测法精度高，但因往往不具备测量条件，只是一个粗略的概算。计算法根据泥石流历时T（s）和最大流量Qc（m3/s），按泥石流暴涨暴落的特点，将其过程线概化成五角形，按式计

Q = = K TQc……………………………… （2—2—6）

式中 K 值的变化随流域面积 (F) 的大小而变化：

当 F < 5 (km2 )时： K = ；

F = 5～10 (km2 )时： K = ；

F = 10～100 (km2 )时： K = ；

泥石流体中的固体物质流量QH（m3）：

QH＝ Q（γC -γW）/（γH –γW）…………………………………(2—2—7)

C、冲击力计算：可用附录I中公式计算泥石流整体冲击力、泥石流中大石块冲击力。泥石流中大石块冲击力的计算方法较多，计算结果可信度较低。

a) 泥石流体整体冲击压力计算公式

F ＝( Υ·γc /g · Vc2 · sinα) · λ………（2—2—8）

式中：

F——泥石流体整体冲击压力（Pa）；

g——重力加速度（m／s2），取g＝9.8m／s2；

α——建筑物受力面与泥石流冲击力方向的夹角（°）；

λ——建筑物形状系数：圆形建筑物λ＝，矩形建筑物λ＝，

方形建筑物：λ＝。

D、弯道超高与冲高计算：泥石流流动在弯曲沟道外侧产生的超高值和泥石流正面遇阻的冲起高度。

a)、泥石流的弯道超高

由于泥石流流速快，惯性大，故在弯道凹岸处有比水流更加显着的弯道超高现象。

根据弯道泥面横比降动力平衡条件，推导出计算弯道超高的公式：

12

2lg

3.2R R g V h c =?……………………………………….（2—2—9） 式中：

Δh —— 弯道超高（m ）；

R2 —— 凹岸曲率半径（m ）；

R1 —— 凸岸曲率半径（m ）；

b)、泥石流最大冲起高度

g V H c 22=

? ……………………………………..…… （2—2—10） c ）泥石流在爬高过程中由于受到沟床阻力的影响，其爬高ΔH ：

g V g bV H c c 2

28

.02≈=?……………………………………..（2—2—11）

（4）泥石流的形成区、流通区和堆积区测绘。

a ）工程治理区实测剖面至少按一纵三横控制；

b ）重点区应有1-3个探槽或探坑（井）控制；

5)、对各类防治工程提供以下主要设计参数。

（1）各类拦挡坝: 覆盖层和基岩的重度、承载力标准值、抗剪强度，基面摩擦系数，泥石流的性质与类型，发生频次，泥石流体的重度和物质组成，泥石流体的流速、流量和设计暴雨洪水频率，泥石流回淤坡度和固体物质颗粒成分，沟床清水冲刷线。

（2）其它工程：桩林着重于其锚固段基岩深度、风化程度和力学性质；排导槽、渡槽着重于泥石流运动的最小坡度、冲击力、弯道超高和冲高；导流堤、护岸堤和防冲墩着重于基岩的埋藏深度和性质、泥石流冲击力和弯道超高、墙背摩擦角；停淤场着重于淤积总量、淤积总高度和分期淤积高度。

2泥石流流体重度(γc)的确定

泥石流流体重度(γc)的确定，目前测试手段还不能解决，因此暂时的替代办法有：

1)、根据数量化评分法确定重度：

根据数量化评分（N ）建立与重度(γc)的线性关系，制成备查表直接查找，参见表2—3--1

2)、现场仿制：可靠性很难保证。

现场请当地曾亲眼看见过该沟泥石流暴发的老居民多人次，在需要测试的沟段，选取有代表性的堆积物搅拌成暴发时的泥石流流体状态，进行样品仿制，然后分别测出样品的总质量和总体积，按下式求出泥石流流体重度。

γc

V G c

……………………………………………………..（2—3--1）

式中：

γc ——泥石流流体重度，t ／m3；

G——样品的总质量，t；

c

V——样品的总体积，m3。

3)、现场评估

在泥石流沟现场，请多人次当地亲自目睹过该沟泥石流暴发或受过灾害的村民，描述泥石流流体特征和流体运动状况。然后按表2—3--2的特征确定泥石流流体重度。

表2—3—2 泥石流流体稠度特征表

表2—3—1 数量化评分（N）与重度、（1＋?）关系对照表 (?h=