CN104794362B - 一种泥石流断面平均流速的测算方法及应用 - Google Patents

一种泥石流断面平均流速的测算方法及应用 Download PDF

Info

Publication number
CN104794362B
CN104794362B CN201510223986.3A CN201510223986A CN104794362B CN 104794362 B CN104794362 B CN 104794362B CN 201510223986 A CN201510223986 A CN 201510223986A CN 104794362 B CN104794362 B CN 104794362B
Authority
CN
China
Prior art keywords
mud
rock flow
flow
section
rock
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.)
Expired - Fee Related
Application number
CN201510223986.3A
Other languages
English (en)
Other versions
CN104794362A (zh
Inventor
陈剑刚
陈晓清
游勇
赵万玉
于献彬
李云
李昆
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Mountain Hazards and Environment IMHE of CAS
Original Assignee
Institute of Mountain Hazards and Environment IMHE of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Mountain Hazards and Environment IMHE of CAS filed Critical Institute of Mountain Hazards and Environment IMHE of CAS
Priority to CN201510223986.3A priority Critical patent/CN104794362B/zh
Publication of CN104794362A publication Critical patent/CN104794362A/zh
Application granted granted Critical
Publication of CN104794362B publication Critical patent/CN104794362B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measuring Volume Flow (AREA)
  • Underground Or Underwater Handling Of Building Materials (AREA)

Abstract

本发明公开了一种泥石流断面平均流速的测算方法及应用。所述测算方法通过现场调查测量、实际取样实测容重和室内流变实验等手段,分别确定泥石流流体容重、泥石流固体物质容重、泥石流沟道或排导槽坡度、泥石流过流断面水力半径、泥石流沟道或排导槽外部糙率、和泥石流粘滞系数等参数,将获取的参数代入泥石流断面平均流速计算公式,从而获得较为准确的泥石流断面平均流速。该方法综合考虑了泥石流沟道或排导槽形状变化引起的外部糙率变化和泥石流流体性质产生的内部糙率,能够合理确定泥石流断面平均流速,为泥石流勘查和防治工程设计提供依据,适应泥石流野外勘察和实际工程设计需要。

Description

一种泥石流断面平均流速的测算方法及应用
技术领域
本发明涉及一种泥石流断面平均流速的测算方法,及其在泥石流天然沟道条件下或在泥石流排导槽设计合理性校核中的应用。
背景技术
泥石流是中国山区一种常见的地质灾害,具有发生突然、历时短暂、来势凶猛、大冲淤、破坏力极强的特点。近年来,在趋于增强的地震活动扰动下,加之全球气候变化导致的极端天气显著增加,泥石流活动非常活跃,重大泥石流灾害频繁发生,人民生命财产安全受到严重危害和威胁。为了保障山区经济的可持续发展,泥石流防治治理就显得十分必要。
泥石流流速是泥石流动力学性质中最为重要的参数,也是泥石流防治工程设计不可缺少的参数。泥石流流速关系到泥石流流量、冲击力、运动摩阻、弯道超高等重要工程设计参数的计算。目前确定泥石流断面平均流速的计算方法,经验或半经验的计算公式有数十种之多。例如:云南东川蒋家沟粘性泥石流估算公式,云南东川大白泥沟和蒋家沟粘性泥石流经验公式,西藏波密古乡沟粘性泥石流估算公式等。各家经验公式皆以水力学中的曼宁公式为原型,但均没有考虑粘滞系数对流速的影响。
发明内容
本发明的目的就是针对现有技术的不足,提供一种以曼宁公式为原型的泥石流断面平均流速测算方法,该方法综合考虑了泥石流沟道或排导槽形状变化引起的外部糙率变化和泥石流流体性质产生的内部糙率,计算所需参数能够通过野外调查、取样和实验分析来获得,能够合理确定泥石流断面平均流速,为泥石流勘查和防治工程设计提供依据,适应泥石流野外勘察和实际工程设计需要。
为实现上述目的,本发明的技术方案是:
本发明提出一种泥石流断面平均流速的测算方法,所述测算方法的推导过程如下:泥石流断面平均流速设计取值与过流时波动范围容许值是表达泥石流运动与动力状况之第一要素,它规范了泥石流过流时的不冲、不淤势态暨对导向物、障碍物的受(承)力特性。泥石流归属固液两相流,其流动形态、动力来源和阻力功耗等与水流相似,将两者类比,泥石流流速计算公式与水流流速计算公式具有相关性,结合水力学中曼宁公式,故泥石流断面平均流速表达公式应具有统一性,有
公式1
其中,1/nTotal为综合糙率系数,表示泥石流在运动过程中的综合阻力,考虑到泥石流受到的综合阻力来自于外部沟道和泥石流自身,因此将综合糙率分解为泥石流沟道或排导槽形状变化造成的外部糙率和泥石流本身性质产生的内部糙率,其中内部糙率考虑泥沙特性和泥石流粘性作用两方面,即:
公式2
因此有:
公式3
公式3中,
vC-泥石流断面平均流速,单位m/s;
k-考虑泥石流体粘性作用的修正系数;
γH-泥石流固体物质容重,单位t/m3
ФC-泥石流泥沙修正系数;其中γC为泥石流流体容重、单位t/m3,γH为泥石流固体物质容重、单位t/m3,γW为水的容重、取值1.0t/m3
noutside-泥石流沟道或排导槽外部糙率;如果是天然沟道,可根据野外调查获得泥石流沟道外部糙率;如果是排导槽,需按照排导槽的材料性质及其附属结构特点等确定泥石流排导槽外部糙率;
RC-泥石流过流断面水力半径,单位m;如果是天然沟道,可采用平均泥深代替;如果是排导槽,由于排导槽深宽比要比天然沟道大得多,其水力半径不能采用平均泥深代替,应为过水断面面积除以湿周;
I-泥石流沟道或排导槽坡度;如果是天然沟道,就是沟道坡度(即沟床纵坡);如果是排导槽,就是排导槽坡度(即排导槽纵坡降)。
接下来,确定考虑泥石流体粘性作用的修正系数k。泥石流浆体对泥石流体的运动速度起着至关重要的作用,随着泥石流容重的增加泥石流体粘性的作用也产生不同的变化,进而泥石流的运动流速也会产生相应的变化。高容重的粘性泥石流能够在坡度较小的沟道内运动,其中泥石流体具有的粘性给泥石流运动流速的确定带来了困难。根据上述分析,泥石流沟道或排导槽的外部糙率和泥石流泥沙修正系数易于确定,那么考虑泥石流体粘性作用的修正系数k的取值范围可通过观测数据来确定。根据东川泥石流观测站泥石流取样测试得到的实测粘度系数η和修正系数k进行拟合得到以下关系(泥石流流体容重γC>1.7t/m3):
k=1.7532η0.3941,R2=0.7552 公式4
将公式4代入公式3可得:
公式5
本发明在将综合糙率分解为外部糙率和考虑泥沙特性、泥石流粘性作用的内部糙率的基础上,提出了一种泥石流断面平均流速的测算方法。具体而言,所述泥石流断面平均流速的测算方法步骤如下:
(一)通过现场调查实测,确定泥石流沟道或排导槽坡度I,泥石流过流断面水力半径Rc、单位m,泥石流沟道或排导槽外部糙率noutside;通过对泥石流体的流变实验,确定泥石流粘滞系数η,单位Pa.s。当测算方法用于泥石流天然沟道时,泥石流沟道外部糙率可通过野外调查获得,一般情况下泥石流沟道外部糙率noutside取值0.033;当测算方法用于泥石流排导槽时,需按照排导槽的材料性质及其附属结构特点等确定泥石流排导槽外部糙率noutside
(二)通过实际取样实测容重,确定泥石流流体容重γC和泥石流固体物质容重γH,单位t/m3;将上述参量代入公式确定泥石流泥沙修正系数ФC,式中,γW为水的容重,取值1.0t/m3
(三)通过以下公式确定泥石流断面平均流速vC
式中,vC-泥石流断面平均流速,单位m/s;
η-泥石流粘滞系数,单位Pa.s,由步骤(一)确定;
γH-泥石流固体物质容重,单位t/m3,由步骤(二)确定;
ФC-泥石流泥沙修正系数,由步骤(二)确定;
noutside-泥石流沟道或排导槽外部糙率,由步骤(一)确定;
Rc-泥石流过流断面水力半径,单位m,由步骤(一)确定;
I-泥石流沟道或排导槽坡度,由步骤(一)确定。
所述泥石流断面平均流速测算方法适用于泥石流天然沟道条件下,或用于泥石流排导槽设计合理性的校核。当用于泥石流排导槽设计合理性的校核时,通过现场调查,确定泥石流峰值流量Qc;通过所述泥石流断面平均流速测算方法,确定排导槽中的泥石流断面平均流速vC,然后计算得到排导槽的过流流量Q;将排导槽的过流流量Q与泥石流峰值流量Qc进行比对,如果Q大于Qc,则所设计的排导槽合理。所述泥石流断面平均流速测算方法适用于泥石流流体容重γC大于1.7t/m3、泥石流沟道或排导槽坡度I为0.1-0.4的泥石流断面平均流速测算。
与现有技术相比,本发明的有益效果是:在泥石流断面平均流速的测算中既考虑了泥石流沟道或排导槽形状变化造成的外部糙率,也考虑了泥石流本身性质产生的内部糙率,即综合考虑了外部条件和泥石流自身性质,同时能够通过切实可行的方法获得具体的计算参数,将获取的参数代入泥石流流速计算公式,从而获得较为准确的泥石流断面平均流速。
具体实施方式
下面对本发明的优选实施例作进一步的描述。
实施例一
某泥石流沟位于绵远河左岸,流域面积约为1.36km2,主沟上游山体地形陡峻,最高点位于流域东侧,海拔1987m;最低点位于沟口绵远河附近,海拔810m,相对高差达1177m;该泥石流沟主沟长度2.59km,沟床纵坡比降412‰。该流域平面形态呈桃叶形,有利于泥石流体汇流。地震后发生多次泥石流,沟道下切严重,切割深度约30m,沟道形态呈“V”型,平均宽度约3~5m。陡峻的地形条件为暴雨洪水的汇集提供了良好的条件,同时较好的临空条件为沟域内不良地质现象的发育以及泥石流松散固体物源的汇集提供了有利的条件。加之沟谷纵坡大,为松散固体物质的搬运和泥石流的形成提供了极为有利的地形。
经野外调查计算得知,20年一遇的泥石流峰值流量约为115.3m3/s。为了减轻泥石流灾害对沟口公路安全畅通的影响,拟在泥石流沟口修建排导槽跨越公路。该排导槽过水断面设计为矩形,宽度为6m,高度为3.5m,20年一遇泥石流设计泥深为2.5m,排导槽的排导纵比降为0.35,由于该排导槽的坡度很大,在排导槽内布置消能结构物。接下来,利用本发明的泥石流断面平均流速测算方法进行泥石流排导槽设计合理性的校核,具体步骤如下:
第一步,通过现场调查实测,确定泥石流排导槽坡度I为0.35,泥石流排导槽外部糙率noutside为0.063,泥石流过流断面水力半径Rc为过水断面面积除以湿周、即为通过对泥石流体的流变实验,确定泥石流粘滞系数η为0.72Pa.s。
第二步,通过实际取样实测容重,确定泥石流流体容重γC为1.9t/m3、泥石流固体物质容重γH为2.65t/m3;将上述参量代入公式确定泥石流泥沙修正系数ФC为1.2。
第三步,将第一步和第二步中得到的参量分别代入以下公式,确定泥石流断面平均流速vC
经计算得到排导槽中的泥石流断面平均流速为8.70m/s,过水断面面积为2.5×6=15m2,则泥石流排导槽的过流流量Q为8.70×15=130.5m3/s。在排导20年一遇的泥石流时,泥石流峰值流量Qc为115.3m3/s,泥石流排导槽的过流流量Q大于泥石流峰值流量Qc,表明该泥石流排导槽设计合理,能够满足实际工程要求。
实施例二
某泥石流沟位于四川省绵竹市西北部山区的清平场镇北,属长江流域的沱江水系上游绵远河左岸一支沟,沟口坐标N31.5°,E104.1°。在地貌上属构造侵蚀中切割陡峻低-中山地貌、斜坡冲沟地形。该流域总体东西向伸展,汇水面积7.81km2,主沟全长3.25km,流域内最低点位于沟口海拔883m,最高峰位于东部分水岭九顶山的顶子崖,海拔2402m,相对高差1519m。该沟多次暴发超大规模泥石流,对清平乡场镇构成极大威胁。
利用本发明的泥石流断面平均流速测算方法在泥石流天然沟道条件下进行泥石流流速的测算,具体步骤如下:
第一步,通过现场调查实测,确定泥石流沟道坡度I为0.11,泥石流过流断面水力半径Rc采用平均泥深代替、为2.70m,泥石流沟道外部糙率noutside为0.033;通过对泥石流体的流变实验,确定泥石流粘滞系数η为0.79Pa.s。
第二步,通过实际取样实测容重,确定泥石流流体容重γC为2.25t/m3、泥石流固体物质容重γH为2.7t/m3;将上述参量代入公式确定泥石流泥沙修正系数ФC为2.78。
第三步,将第一步和第二步中得到的参量分别代入以下公式,确定泥石流断面平均流速vC
通过现场调查实测,过水断面面积为112m2,通过进一步计算可得通过该断面的泥石流流量为112×8.81=987m3;而通过原型观测的录像分析获得的泥石流流速为10.19m/s,泥石流流量为1141m3/s,误差约为14%,能够满足泥石流防治工程设计的需要。

Claims (5)

1.一种泥石流断面平均流速的测算方法,其特征在于:针对泥石流流体容重γC大于1.7t/m3、泥石流沟道或排导槽坡度I为0.1-0.4,所述泥石流断面平均流速的测算方法步骤如下:
(一)通过现场调查实测,确定泥石流沟道或排导槽坡度I,泥石流过流断面水力半径Rc、单位m,泥石流沟道或排导槽外部糙率noutside;通过对泥石流体的流变实验,确定泥石流粘滞系数η,单位Pa.s;
(二)通过实际取样实测容重,确定泥石流流体容重γC和泥石流固体物质容重γH,单位t/m3;将上述参量代入公式确定泥石流泥沙修正系数ΦC,式中,γW为水的容重,取值1.0t/m3
(三)通过以下公式确定泥石流断面平均流速vC
式中,vC-泥石流断面平均流速,单位m/s;
η-泥石流粘滞系数,单位Pa.s,由步骤(一)确定;
γH-泥石流固体物质容重,单位t/m3,由步骤(二)确定;
ΦC-泥石流泥沙修正系数,由步骤(二)确定;
noutside-泥石流沟道或排导槽外部糙率,由步骤(一)确定;
Rc-泥石流过流断面水力半径,单位m,由步骤(一)确定;
I-泥石流沟道或排导槽坡度,由步骤(一)确定。
2.根据权利要求1所述泥石流断面平均流速的测算方法,其特征在于:所述泥石流断面平均流速测算方法用于泥石流天然沟道时,泥石流沟道外部糙率noutside取值0.033。
3.根据权利要求1所述泥石流断面平均流速测算方法,其特征在于:适用于泥石流天然沟道条件下,或用于泥石流排导槽设计合理性的校核。
4.根据权利要求3所述泥石流断面平均流速测算方法,其特征在于:当用于泥石流排导槽设计合理性的校核时,通过现场调查,确定泥石流峰值流量Qc;通过所述泥石流断面平均流速测算方法,确定排导槽中的泥石流断面平均流速vC,然后计算得到排导槽的过流流量Q;将排导槽的过流流量Q与泥石流峰值流量Qc进行比对,如果Q大于Qc,则所设计的排导槽合理。
5.根据权利要求1所述泥石流断面平均流速测算方法,其特征在于:适用于泥石流流体容重γC大于1.7t/m3、泥石流沟道或排导槽坡度I为0.1-0.4的泥石流断面平均流速测算。
CN201510223986.3A 2015-05-06 2015-05-06 一种泥石流断面平均流速的测算方法及应用 Expired - Fee Related CN104794362B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510223986.3A CN104794362B (zh) 2015-05-06 2015-05-06 一种泥石流断面平均流速的测算方法及应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510223986.3A CN104794362B (zh) 2015-05-06 2015-05-06 一种泥石流断面平均流速的测算方法及应用

Publications (2)

Publication Number Publication Date
CN104794362A CN104794362A (zh) 2015-07-22
CN104794362B true CN104794362B (zh) 2018-06-05

Family

ID=53559153

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510223986.3A Expired - Fee Related CN104794362B (zh) 2015-05-06 2015-05-06 一种泥石流断面平均流速的测算方法及应用

Country Status (1)

Country Link
CN (1) CN104794362B (zh)

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105115695B (zh) * 2015-08-26 2017-05-31 中交第二航务工程局有限公司 用于测定泥水盾构管路临界流速的方法
CN105256768B (zh) * 2015-09-09 2017-03-22 中国科学院水利部成都山地灾害与环境研究所 箱体消能式泥石流排导槽的箱体消能段设计方法
CN105678062B (zh) * 2015-12-30 2018-11-09 广州地理研究所 泥石流流量和容重沿程变化预测方法
CN105926542B (zh) 2016-05-13 2018-06-12 中国科学院水利部成都山地灾害与环境研究所 一种非对称式泥石流排导槽的设计方法和应用
CN106245573B (zh) * 2016-07-19 2018-08-21 成都理工大学 一种粘性泥石流弯道超高计算方法及应用
CN106446411B (zh) * 2016-09-26 2019-06-28 中国科学院、水利部成都山地灾害与环境研究所 粘性泥石流内部流速、冲击力空间分布测量方法
CN106638454B (zh) * 2017-01-18 2018-12-14 中国地质科学院探矿工艺研究所 粘性泥石流v型排导槽流速的计算方法
CN106968217A (zh) * 2017-04-14 2017-07-21 中国地质科学院探矿工艺研究所 基于降雨频率的泥石流淤积危险范围划分方法
CN106951725A (zh) * 2017-05-12 2017-07-14 成都理工大学 一种小流域沟道堆石体失稳的计算方法
CN107179419B (zh) * 2017-05-14 2019-06-14 中国科学院、水利部成都山地灾害与环境研究所 粘性泥石流断面平均流速测算方法
CN107220513B (zh) * 2017-06-13 2020-09-08 成都理工大学 一种计算泥石流平均流速的方法及其应用
CN107169615B (zh) * 2017-07-19 2020-04-10 四川建筑职业技术学院 一种拦砂坝淤满条件下泥石流流速衰减值的计算方法
CN107357998B (zh) * 2017-07-19 2020-08-25 四川建筑职业技术学院 一种泥石流堆积区流量衰减值的计算方法
CN107169252B (zh) * 2017-07-19 2020-05-19 四川建筑职业技术学院 一种基于泥石流固体颗粒粒径的震区泥石流容重计算方法
CN108221842B (zh) * 2018-01-31 2020-01-14 中国科学院、水利部成都山地灾害与环境研究所 泥石流格子坝过坝流速、流量计算方法
CN108797533B (zh) * 2018-06-14 2020-05-19 中国科学院、水利部成都山地灾害与环境研究所 一种泥石流拦砂坝圆形泄流孔的设计方法
CN109766513B (zh) * 2019-01-10 2023-02-28 中国科学院、水利部成都山地灾害与环境研究所 一种泥石流断面平均流速的测算方法
CN115128591B (zh) * 2022-06-08 2023-06-20 中国地质环境监测院(自然资源部地质灾害技术指导中心) 泥石流监测雷达参数验证方法
CN116519264B (zh) * 2023-06-30 2024-03-12 中国科学院、水利部成都山地灾害与环境研究所 泥石流沟道断面综合糙率系数的动态测定方法及装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
浆体与级配颗粒组合条件下泥石流冲击特性实验研究;何晓英;《中国博士学位论文全文数据库(电子期刊)》;20150228;第2015年卷(第2期);第A011-15页 *
粘性泥石流流速计算方法;朱兴华等;《泥沙研究》;20130630;第2013年卷(第3期);第59-64页 *
粘性泥石流的平均运动速度研究;余斌;《地球科学进展》;20080531;第23卷(第5期);第524-532页 *
粘性泥石流运动流速与流量计算;舒安平等;《泥沙研究》;20030630;第2003年卷(第3期);第7-11页 *

Also Published As

Publication number Publication date
CN104794362A (zh) 2015-07-22

Similar Documents

Publication Publication Date Title
CN104794362B (zh) 一种泥石流断面平均流速的测算方法及应用
Wu et al. Effect of altered flow regime on bankfull area of the Lower Yellow River, China
Hwang et al. The physical processes in the Yellow Sea
CN107727355B (zh) 一种流体挑出水平距离测算方法及其应用
Sousa et al. Influence of the Minho River plume on the Rias Baixas (NW of the Iberian Peninsula)
Shuwei et al. Morphology and mechanism of the very large dunes in the tidal reach of the Yangtze River, China
Al-Ansari et al. Flow of River Tigris and its effect on the bed sediment within Baghdad, Iraq
Zhang et al. Seasonal and interannual variations of flow discharge from Pearl River into sea
Jia et al. Wave-forced sediment erosion and resuspension in the Yellow River Delta
Yu et al. Numerical analysis on the Brisbane River plume in Moreton Bay due to Queensland floods 2010–2011
Ismagilov et al. Hydraulic parameters on the curvilinear section of the river channel in conditions of regulated water flow
Sadeghi et al. Partitioning and analyzing temporal variability of wash and bed material loads in a forest watershed in Iran
CN108563807A (zh) 一种粘性泥石流的冲击力测算方法及其应用
Liang et al. Typhoon-induced suspended sediment dynamics in the mouth-bar region of a river/wave-dominated estuary
Hu Impacts of human activities in the Hanjiang River Basin, China
Wiejaczka et al. Water temperature dynamics in a complex of reservoirs and its effect on the temperature patterns of a mountain river
Barahimi et al. Flow structure and shear stress in the presence of both ice cover on water surface and leafless vegetation in channel bed
Lin et al. Tailings dam break flow and sediment numerical simulation
Jia et al. Erosion Survey of the Modern Yellow River Delta
Waterhouse et al. Tidal asymmetries of velocity and stratification over a bathymetric depression in a tropical inlet
Jia et al. Geo-marine environment and sediment properties of the modern Yellow River Delta
Liu et al. Numerical Study on Influence of Water Exchange for Artificial Island Group
Shu et al. Estimation for the riverbank collapse volume with sandy-riverbank in the desert reach of the upper Yellow River
Umeda et al. VERTICAL STRUCTURE OF VELOCITY DISTRIBUTIONS MEASURED IN A THERMALLY STRATIFIED RESERVOIR
Łajczak Transport of suspended sediment by the Vistula River basin upstream of Kraków, southern Poland, and the human impact during the second half of the 20th century

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180605