CN107220513B - 一种计算泥石流平均流速的方法及其应用 - Google Patents
一种计算泥石流平均流速的方法及其应用 Download PDFInfo
- Publication number
- CN107220513B CN107220513B CN201710443829.2A CN201710443829A CN107220513B CN 107220513 B CN107220513 B CN 107220513B CN 201710443829 A CN201710443829 A CN 201710443829A CN 107220513 B CN107220513 B CN 107220513B
- Authority
- CN
- China
- Prior art keywords
- debris flow
- flow
- debris
- velocity
- measuring
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Measuring Volume Flow (AREA)
- Excavating Of Shafts Or Tunnels (AREA)
Abstract
本发明公开了一种计算泥石流平均流速的方法及其应用,属于泥石流防治工程领域,其特征在于,包括以下步骤:a、现场调查确定泥石流屈服应力τ0;b、现场测量泥石流容重ρ;c、现场测量泥石流流深h;d、通过式1计算无量纲化的屈服应力τ;e、现场测量泥石流的表面速度Vs;f、现场测量泥石流运动的底床粗糙率n;g、通过式2计算泥石流平均流速Vm。本发明充分考虑了泥石流运动的底床粗糙率、流深与泥石流性质,通过无量纲化将尺寸效应消除,利于推演到实际发生的泥石流中,具有计算方便、计算准确度高的特点,适用于野外大尺度的实际计算,能够为泥石流评估与防治提供关键的参考依据。
Description
技术领域
本发明涉及到泥石流防治工程领域,尤其涉及一种计算泥石流平均流速的方法及其应用。
背景技术
泥石流的运动速度是泥石流的最重要的参数之一,是评估泥石流和防治泥石流灾害的重要参数。泥石流的运动速度一般包括泥石流的平均流速和表面速度。泥石流的平均流速是代表泥石流速度的主要特征速度,评估和防治泥石流的两个重要参数之一的泥石流峰值流量,就是从泥石流的平均流速计算得来的。但野外条件很难观测到泥石流的平均流速;相反的,泥石流运动的表面速度往往比较容易观测测量,获得表面速度值。如果能通过泥石流表面速度,以及表面速度与平均流速的关系获得泥石流平均流速,也能满足泥石流的评估和防治中的要求,从而解决问题泥石流的评估和防治问题。
目前国内外学者对泥石流的平均流速与表面流速的关系研究相对较少且以流速垂向分布规律为主,由Bagnold的膨胀体模型得出的水石流平均流速与表面流速的关系有:Vm/Vs=0.6(Vm为水石流的平均流速,Vs为水石流的表面流速,费祥俊,舒安平,泥石流运动机理与灾害防治,清华大学出版社,2004,104-108);宏观粘性模型的平均流速与表面流速的关系有:Vm/Vs=0.67(费祥俊,舒安平,泥石流运动机理与灾害防治,清华大学出版社,2004,104-108);由Bingham模型得到的接近于粘性泥石流的粘塑性模型的平均流速与表面流速的关系有:Vm/Vs=0.67(无流核)-1(全部为流核)(仅适用于层流运动,费祥俊,舒安平,泥石流运动机理与灾害防治,清华大学出版社,2004,104-108);由固体颗粒摩擦与碰撞混合模型的平均流速与表面流速的关系有:Vm/Vs=0.6(无流核)-1(全部为流核),(费祥俊,舒安平,泥石流运动机理与灾害防治,清华大学出版社,2004,104-108)。由此可见泥石流的平均流速与表面流速的关系比较复杂,不仅范围很大,从0.6到1;而且随泥石流的性质不同,差别也很大。此外,这些研究中没有直接考虑泥石流的底部粗糙率、流深及泥石流性质这3个关键因素,仅从泥石流运动的机理上可以看出这3个因素的作用:底部粗糙率越大,底部速度越小,底部速度的影响区越大,平均流速与表面流速的比值越小;泥石流流深越大,底部较小速度的影响范围相对越小,平均流速与表面流速的比值越大;泥石流性质的影响主要体现在泥石流的底部阻力:从清水到稀性泥石流,再到粘性泥石流,粘滞阻力随粘性的增加逐渐增大,底部阻力也逐渐增大,底部速度也逐渐变小,平均流速与表面流速的比值逐渐减小。
公开号为CN 104794362A,公开日为2015年07月22日的中国专利文献公开了一种泥石流断面平均流速的测算方法及应用,所述测算方法通过现场调查测量、实际取样实测容重和室内流变实验等手段,分别确定泥石流流体容重、泥石流固体物质容重、泥石流沟道或排导槽坡度、泥石流过流断面水力半径、泥石流沟道或排导槽外部糙率、和泥石流粘滞系数等参数,将获取的参数代入泥石流断面平均流速计算公式,从而获得较为准确的泥石流断面平均流速。
该专利文献公开的泥石流断面平均流速的测算方法及应用,综合考虑了泥石流沟道或排导槽形状变化引起的外部糙率变化和泥石流流体性质产生的内部糙率,能够合理确定泥石流断面平均流速,但是,整个测算方法涉及到的参数较多,测算难度较大,不适用于野外大尺度的实际计算。
发明内容
本发明为了克服上述现有技术的缺陷,提供一种计算泥石流平均流速的方法及其应用,本发明充分考虑了泥石流运动的底部粗糙率、流深与泥石流性质,通过无量纲化将尺寸效应消除,利于推演到实际发生的泥石流中,具有计算方便、计算准确度高的特点,适用于野外大尺度的实际计算,能够为泥石流评估与防治提供关键的参考依据。
本发明通过下述技术方案实现:
一种计算泥石流平均流速的方法,其特征在于,包括以下步骤:
a、现场调查确定泥石流屈服应力τ0,单位Pa;
b、现场测量泥石流容重ρ,单位kg/m3;
c、现场测量泥石流流深h,单位m;
d、通过式1计算无量纲化的屈服应力τ;
τ=τ0/(ρgh) (式1)
其中,g为重力加速度,g=9.8m/s2;
e、现场测量泥石流的表面速度Vs,单位m/s;
f、现场测量泥石流运动的底床粗糙率n;
g、通过式2计算泥石流的平均流速Vm,单位m/s;
Vm/Vs=(18.63n2-3.53n+0.6)τ-0.07 (式2)。
所述步骤f中,泥石流运动的底床粗糙率n的取值范围为0.02-0.09。
本发明适用于稀性泥石流和粘性泥石流的平均流速计算。
本发明适用于不同屈服应力泥石流的平均流速计算。
本发明的基本原理如下:
泥石流的速度是泥石流评估与防治中最重要的参数之一。往往采用泥石流的平均流速作为特征量,作为泥石流评估与泥石流防治的重要依据。
泥石流的平均流速很难测量,但是泥石流的表面速度能比较容易地测量到。如果能测量泥石流的表面速度,同时有表面速度计算平均流速的方法,就可以得到泥石流的平均流速,可以为泥石流的评估与防治提供重要的依据。
泥石流的平均流速与表面速度之间的关系,受泥石流运动底床的粗糙率和泥石流性质的影响:泥石流运动底床的粗糙率越大,泥石流的下部速度越小,泥石流的平均流速越小,平均流速与表面速度之比越小。
无量纲化的泥石流屈服应力(泥石流屈服应力与泥石流容重、重力加速度、泥石流水深之比)可以表达泥石流的粘性性质:泥石流的屈服应力与泥石流容重越大,泥石流的粘性越强,泥石流的内部阻力越大,泥石流的平均流速与表面速度之比越小;但是泥石流的屈服应力的变化远大于泥石流的容重变化,因此泥石流的屈服应力与泥石流的容重之比越大,泥石流的粘性越强,泥石流的内部阻力越大,泥石流的平均流速与表面速度之比越小;泥石流运动水深越大,泥石流下部受底床粗糙率的影响范围相对越小,平均流速与表面速度之比越大。综合分析可以得出:无量纲化的泥石流屈服应力越大,泥石流的平均流速与表面速度之比越小。
本发明的有益效果主要表现在以下方面:
一、本发明,“a、现场调查确定泥石流屈服应力τ0;b、现场测量泥石流容重ρ;c、现场测量泥石流流深h;d、通过式1计算无量纲化的屈服应力τ;e、现场测量泥石流的表面速度Vs;f、现场测量泥石流运动的底床粗糙率n;g、通过式2计算泥石流平均流速Vm”,作为一个完整的技术方案,本发明充分考虑了泥石流运动的底床粗糙率、流深与泥石流性质,通过无量纲化将尺寸效应消除,利于推演到实际发生的泥石流中,具有计算方便、计算准确度高的特点,适用于野外大尺度的实际计算,能够为泥石流评估与防治提供关键的参考依据。
二、本发明,将泥石流运动的底床粗糙率、流深与泥石流性质综合起来考虑得出的计算式,对于使用者更加方便,普适性好。
三、本发明,以泥石流的流深、泥石流屈服应力和容重、底床糙率等为基本参数,通过实验研究得出泥石流平均流速与表面流速关系,对于泥石流防治及将来的理论研究有更好的指导作用,具有更高的防灾适用性。
具体实施方式
实施例1
一种计算泥石流平均流速的方法,包括以下步骤:
a、现场调查确定泥石流屈服应力τ0,单位Pa;
b、现场测量泥石流容重ρ,单位kg/m3;
c、现场测量泥石流流深h,单位m;
d、通过式1计算无量纲化的屈服应力τ;
τ=τ0/(ρgh) (式1)
其中,g为重力加速度,g=9.8m/s2;
e、现场测量泥石流的表面速度Vs,单位m/s;
f、现场测量泥石流运动的底床粗糙率n;
g、通过式2计算泥石流的平均流速Vm,单位m/s;
Vm/Vs=(18.63n2-3.53n+0.6)τ-0.07 (式2)。
“a、现场调查确定泥石流屈服应力τ0;b、现场测量泥石流容重ρ;c、现场测量泥石流流深h;d、通过式1计算无量纲化的屈服应力τ;e、现场测量泥石流的表面速度Vs;f、现场测量泥石流运动的底床粗糙率n;g、通过式2计算泥石流平均流速Vm”,作为一个完整的技术方案,本发明充分考虑了泥石流运动的底床粗糙率、流深与泥石流性质,通过无量纲化将尺寸效应消除,利于推演到实际发生的泥石流中,具有计算方便、计算准确度高的特点,适用于野外大尺度的实际计算,能够为泥石流评估与防治提供关键的参考依据。
实施例2
一种计算泥石流平均流速的方法,包括以下步骤:
a、现场调查确定泥石流屈服应力τ0,单位Pa;
b、现场测量泥石流容重ρ,单位kg/m3;
c、现场测量泥石流流深h,单位m;
d、通过式1计算无量纲化的屈服应力τ;
τ=τ0/(ρgh) (式1)
其中,g为重力加速度,g=9.8m/s2;
e、现场测量泥石流的表面速度Vs,单位m/s;
f、现场测量泥石流运动的底床粗糙率n;
g、通过式2计算泥石流的平均流速Vm,单位m/s;
Vm/Vs=(18.63n2-3.53n+0.6)τ-0.07 (式2)。
所述步骤f中,泥石流运动的底床粗糙率n的取值范围为0.02。
实施例3
一种计算泥石流平均流速的方法,包括以下步骤:
a、现场调查确定泥石流屈服应力τ0,单位Pa;
b、现场测量泥石流容重ρ,单位kg/m3;
c、现场测量泥石流流深h,单位m;
d、通过式1计算无量纲化的屈服应力τ;
τ=τ0/(ρgh) (式1)
其中,g为重力加速度,g=9.8m/s2;
e、现场测量泥石流的表面速度Vs,单位m/s;
f、现场测量泥石流运动的底床粗糙率n;
g、通过式2计算泥石流的平均流速Vm,单位m/s;
Vm/Vs=(18.63n2-3.53n+0.6)τ-0.07 (式2)。
所述步骤f中,泥石流运动的底床粗糙率n的取值范围为0.04。
将泥石流运动的底床粗糙率、流深与泥石流性质综合起来考虑得出的计算式,对于使用者更加方便,普适性好。
实施例4
一种计算泥石流平均流速的方法,包括以下步骤:
a、现场调查确定泥石流屈服应力τ0,单位Pa;
b、现场测量泥石流容重ρ,单位kg/m3;
c、现场测量泥石流流深h,单位m;
d、通过式1计算无量纲化的屈服应力τ;
τ=τ0/(ρgh) (式1)
其中,g为重力加速度,g=9.8m/s2;
e、现场测量泥石流的表面速度Vs,单位m/s;
f、现场测量泥石流运动的底床粗糙率n;
g、通过式2计算泥石流的平均流速Vm,单位m/s;
Vm/Vs=(18.63n2-3.53n+0.6)τ-0.07 (式2)。
所述步骤f中,泥石流运动的底床粗糙率n的取值范围为0.09。
本发明适用于稀性泥石流和粘性泥石流的平均流速计算。
本发明适用于不同屈服应力泥石流的平均流速计算。
以泥石流的流深、泥石流屈服应力和容重、底床糙率等为基本参数,通过实验研究得出泥石流平均流速与表面流速关系,对于泥石流防治及将来的理论研究有更好的指导作用,具有更高的防灾适用性。
本发明用于泥石流水槽验证,泥石流平均流速与表面流速关系试验数据参见表1。
表1
通过表1可以看出,公式2在泥石流试验的验证中与实测数据吻合度很好,最大误差为-7.8%(仅1个,占6.25%),其余误差在4.92%(绝对值)以下。故该泥石流平均流速与表面流速比值的计算公式可靠度高,计算准确性好。
Claims (3)
1.一种计算泥石流平均流速的方法,其特征在于,包括以下步骤:
a、现场调查确定泥石流屈服应力τ0,单位Pa;
b、现场测量泥石流容重ρ,单位kg/m3;
c、现场测量泥石流流深h,单位m;
d、通过式1计算无量纲化的屈服应力τ;
τ=τ0/(ρgh) (式1)
其中,g为重力加速度,g=9.8m/s2;
e、现场测量泥石流的表面速度Vs,单位m/s;
f、现场测量泥石流运动的底床粗糙率n;
g、通过式2计算泥石流的平均流速Vm,单位m/s;
Vm/Vs=(18.63n2-3.53n+0.6)τ-0.07 (式2)。
2.根据权利要求1所述的一种计算泥石流平均流速的方法,其特征在于:所述步骤f中,泥石流运动的底床粗糙率n的取值范围为0.02-0.09。
3.根据权利要求1所述的一种计算泥石流平均流速的方法的应用,其特征在于:适用于不同屈服应力泥石流的平均流速计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710443829.2A CN107220513B (zh) | 2017-06-13 | 2017-06-13 | 一种计算泥石流平均流速的方法及其应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710443829.2A CN107220513B (zh) | 2017-06-13 | 2017-06-13 | 一种计算泥石流平均流速的方法及其应用 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107220513A CN107220513A (zh) | 2017-09-29 |
CN107220513B true CN107220513B (zh) | 2020-09-08 |
Family
ID=59948465
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710443829.2A Active CN107220513B (zh) | 2017-06-13 | 2017-06-13 | 一种计算泥石流平均流速的方法及其应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107220513B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109766513B (zh) * | 2019-01-10 | 2023-02-28 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种泥石流断面平均流速的测算方法 |
CN112100567B (zh) * | 2020-09-11 | 2023-06-06 | 中国科学院、水利部成都山地灾害与环境研究所 | 稀性泥石流垂向流速分布、表面流速、平均流速测量方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015113567A (ja) * | 2013-12-09 | 2015-06-22 | 一般財団法人砂防・地すべり技術センター | 土石流の移動の解析方法 |
CN104794362A (zh) * | 2015-05-06 | 2015-07-22 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种泥石流断面平均流速的测算方法及应用 |
CN104848825A (zh) * | 2015-05-15 | 2015-08-19 | 中国科学院水利部成都山地灾害与环境研究所 | 泥石流排导槽肋槛后最大冲刷深度的测算方法及应用 |
WO2016137886A1 (en) * | 2015-02-24 | 2016-09-01 | Hayward Industries, Inc. | Pool cleaner with optical out-of-water and debris detection |
CN106096320A (zh) * | 2016-07-19 | 2016-11-09 | 成都理工大学 | 一种泥石流的平衡坡度测算方法及其应用 |
CN106683184A (zh) * | 2017-01-04 | 2017-05-17 | 朱军 | 网络环境下泥石流灾害过程快速模拟与可视化分析方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7564741B2 (en) * | 2004-04-06 | 2009-07-21 | Newsco Directional And Horizontal Drilling Services Inc. | Intelligent efficient servo-actuator for a downhole pulser |
-
2017
- 2017-06-13 CN CN201710443829.2A patent/CN107220513B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015113567A (ja) * | 2013-12-09 | 2015-06-22 | 一般財団法人砂防・地すべり技術センター | 土石流の移動の解析方法 |
WO2016137886A1 (en) * | 2015-02-24 | 2016-09-01 | Hayward Industries, Inc. | Pool cleaner with optical out-of-water and debris detection |
CN104794362A (zh) * | 2015-05-06 | 2015-07-22 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种泥石流断面平均流速的测算方法及应用 |
CN104848825A (zh) * | 2015-05-15 | 2015-08-19 | 中国科学院水利部成都山地灾害与环境研究所 | 泥石流排导槽肋槛后最大冲刷深度的测算方法及应用 |
CN106096320A (zh) * | 2016-07-19 | 2016-11-09 | 成都理工大学 | 一种泥石流的平衡坡度测算方法及其应用 |
CN106683184A (zh) * | 2017-01-04 | 2017-05-17 | 朱军 | 网络环境下泥石流灾害过程快速模拟与可视化分析方法 |
Non-Patent Citations (5)
Title |
---|
Adaptive computer-based spatial-filtering method for more accurate estimation of the surface velocity of debris flow;MS Uddin 等;《appliedoptics》;19991231;第28卷(第32期);第6714-6721页 * |
地表糙率对坡面流流速影响的试验研究;施明新 等;《水力发电学报》;20150630;第34卷(第6期);第117-124页 * |
栾川县泥石流容重与流速参数应用研究;左伟 等;《地球科技》;20140731;第38-40页 * |
泥石流流速公式探讨;洪正修;《中国地质灾害与防治学报》;19960331;第7卷(第1期);第26-32页 * |
粘性泥石流的平均运动速度研究;余斌;《地球科学进展》;20080531;第23卷(第4期);第524-532页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107220513A (zh) | 2017-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Strom et al. | ADV measurements around a cluster microform in a shallow mountain stream | |
Kaushal et al. | Experimental investigation for near-wall lift of coarser particles in slurry pipeline using γ-ray densitometer | |
Righetti et al. | Particle–fluid interactions in a plane near-wall turbulent flow | |
Zhang et al. | Correction factor to dye-measured flow velocity under varying water and sediment discharges | |
Briceno et al. | Self-lubricated transport of aqueous foams in horizontal conduits | |
CN103529237B (zh) | 一种泥沙群体沉速的测量方法及测量装置 | |
CN107220513B (zh) | 一种计算泥石流平均流速的方法及其应用 | |
Schatzmann et al. | Rheometry for large-particulated fluids: analysis of the ball measuring system and comparison to debris flow rheometry | |
Sahin et al. | Observations of suspended sediment stratification from acoustic backscatter in muddy environments | |
Madadi et al. | Investigation of flow characteristics above trapezoidal broad-crested weirs | |
Matoušek et al. | Validation of transport and friction formulae for upper plane bed by experiments in rectangular pipe | |
Qian et al. | Driven cavity flow: from molecular dynamics to continuum hydrodynamics | |
García et al. | Flow, turbulence, and resistance in a flume with simulated vegetation | |
McLean et al. | Double-averaged velocity profiles over fixed dune shapes | |
CN116838327A (zh) | 一种气水油三相流井筒压降预测方法 | |
Pal et al. | Mathematical model on grain-size distribution in suspension over sand-gravel bed | |
Le Bouteiller et al. | Hydraulics and rheology of natural hyperconcentrated flows from Draix-Bleone observatory, French Alps | |
Carlsen et al. | Simultaneous continuous monitoring of the drilling-fluid friction factor and density | |
Staub et al. | Sediment suspension in oscillatory flow: measurements of instantaneous concentration at high shear | |
Czernuszenko | Drift velocity concept for sediment-laden flows | |
Hoque et al. | Air entrainment and associated energy dissipation in steady and unsteady plunging jets at free surface | |
Zordan et al. | Bed shear stress estimation for gravity currents performed in laboratory | |
Shafquet et al. | Study of void fraction measurement in a two phase flow by using differential pressure and electrical capacitance tomography | |
Crone | The basic sediment transport equations made ridiculously simple | |
RU2583875C1 (ru) | Способ измерения расхода жидкости в скважине (трубопроводе) некалиброванным расходомером |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |