CN114595426B - 一种稀性泥石流堵江判断方法 - Google Patents

一种稀性泥石流堵江判断方法 Download PDF

Info

Publication number
CN114595426B
CN114595426B CN202210032435.9A CN202210032435A CN114595426B CN 114595426 B CN114595426 B CN 114595426B CN 202210032435 A CN202210032435 A CN 202210032435A CN 114595426 B CN114595426 B CN 114595426B
Authority
CN
China
Prior art keywords
river
flow
main
branch
debris flow
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
Application number
CN202210032435.9A
Other languages
English (en)
Other versions
CN114595426A (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.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
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 Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN202210032435.9A priority Critical patent/CN114595426B/zh
Publication of CN114595426A publication Critical patent/CN114595426A/zh
Application granted granted Critical
Publication of CN114595426B publication Critical patent/CN114595426B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services
    • G06Q50/26—Government or public services
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00—Water conservation; Efficient water supply; Efficient water use
    • Y02A20/40—Protecting water resources
    • Y02A20/402—River restoration

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Databases & Information Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Algebra (AREA)
  • Operations Research (AREA)
  • Evolutionary Biology (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Probability & Statistics with Applications (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • General Business, Economics & Management (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种稀性泥石流堵江判断方法，属于地质灾害防治工程技术领域，其特征在于，包括以下步骤：a、实地测量主河与支沟交汇处主河的平均宽度B1和主河与支沟交汇处支沟的平均宽度B2，计算支沟与主河的单宽流量比ε；b、计算稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径；c、计算体积因子yt；d、计算入汇角因子

Description

一种稀性泥石流堵江判断方法
技术领域
本发明涉及到地质灾害防治工程技术领域，尤其涉及一种稀性泥石流堵江判断方法。
背景技术
泥石流是指在山区或者其他沟谷深壑，地形险峻的地区，因为暴雨、暴雪或其他自然灾害引发的山体滑坡并携带有大量泥沙以及石块的特殊洪流。泥石流具有突然性以及流速快，流量大，物质容量大和破坏力强的特点。
目前国内外对泥石流堵江的计算研究，主要有三种方法：
1、以室内水槽试验模拟为主，主要方法是模拟在不同交汇角的条件下，主支沟的交汇及河床响应过程，选择不同的评价指标，例如：流量比、容重和规模，分析得出泥石流的堵河判别式。此类方法得到的公式往往由于野外数据尺度远大于室内试验数据，在量纲不和谐的情况下实际计算存在较大的误差。
2、基于统计学的相关性分析，该方法主要以文献查阅和实地调查得到各泥石流沟的基本数据和堵河特征资料，对堵河型泥石流基本特征进行分析总结，得出计算泥石流堵断主河的定性或动量的经验公式。此类方法得到的公式往往只适用于部分区域，缺乏普适性。
3、建立起泥石流沉积及泥石流与主河水流相互作用动力学方程组，以动量与动能平衡为主要切入点，分析得到其各自的理论解，理论分析泥石流堵断主河的条件。这类方法虽然切入点较前沿，但是从成熟到能用于实际生产还需时日。
公开号为CN 105527405A，公开日为2016年04月27日的中国专利文献公开了一种泥石流汇江物理模拟试验装置，其特征在于，包括第一泥石流支沟沟槽、第二泥石流支沟沟槽、泥石流主沟形成区沟槽、泥石流主沟流通区沟槽、泥石流主沟堆积区沟槽、宽度可变河流主河道槽、三角堰、蓄水池、模拟降雨装置、数据采集装置和移动支撑装置，所述第一泥石流支沟沟槽、所述第二泥石流支沟沟槽、所述泥石流主沟形成区沟槽和所述泥石流主沟流通区沟槽均倾斜设置，所述第一泥石流支沟沟槽的水平低端、第二泥石流支沟沟槽的水平低端、所述泥石流主沟形成区沟槽的水平低端和所述泥石流主沟流通区沟槽的水平高端通过可活动橡胶带固定连接，所述第一泥石流支沟沟槽和所述第二泥石流支沟沟槽分别设置在所述泥石流主沟形成区沟槽的两侧，所述第一泥石流支沟沟槽、所述第二泥石流支沟沟槽、所述泥石流主沟形成区沟槽和所述泥石流主沟流通区沟槽的上方均设有所述模拟降雨装置，所述宽度可变河流主河道槽和所述蓄水池均设在所述移动支撑装置上，所述蓄水池的出水口通过可活动橡胶带与所述宽度可变河流主河道槽的一端连接，所述蓄水池的出水口上设有阀门，所述宽度可变河流主河道槽的另一端与所述三角堰的一端连接，所述泥石流主沟流通区沟槽的水平低端通过所述泥石流主沟堆积区沟槽与所述宽度可变河流主河道槽的侧面连通，所述数据采集装置设在所述移动支撑装置的周围。
该专利文献公开的泥石流汇江物理模拟试验装置，虽然能模拟泥石流暴发及汇江过程，研究不同固体物源、坡面角度、江水流速、泥石流堆积区与江河河道交汇角这些因素对泥石流汇江的影响，但是，由于野外数据尺度远大于室内试验数据，在量纲不和谐的情况下实际计算存在较大的误差，影响泥石流堵江判断的准确性，不利于提高泥石流防治效果。
公开号为CN106683019A，公开日为2017年05月17日的中国专利文献公开了一种泥石流堵江计算方法及其应用，包括以下步骤：测得支沟与主河的交汇角θ；测量主河与支沟交汇处主河的平均宽度B1，主河流量Qm，确定主河单宽流量Q1；测量主河与支沟交汇处支沟的平均宽度B2，泥石流流量Q，确定支沟单宽流量Q2，确定支沟与主河的单宽流量比ε，测量主河与支沟交汇处主河的水深h，确定体积因子yt，测量泥石流屈服应力τ1，泥石流密度ρ，确定泥石流屈服应力因子x；根据判断因子y值判断泥石流是否会造成堵江。
该专利文献公开的泥石流堵江计算方法及其应用，具有量纲和谐、计算结果准确度高的特点，能够对泥石流堵江进行准确判断；但是，其泥石流堵江计算方法有限定条件，即必须与主河正交、流量比较大才行，其普适性还是较差。
发明内容
本发明为了克服上述现有技术的缺陷，提供一种稀性泥石流堵江判断方法，本发明结合室内试验和野外数据分析，结合稀性泥石流特点与堵江过程，将稀性泥石流堵江的影响因素进行无量纲化，能适应野外复杂多变的灾害情况，为稀性泥石流堵江防治提供科学指导，泥石流堵江判断准确，普适性好。
本发明通过下述技术方案实现：
一种稀性泥石流堵江判断方法，其特征在于，包括以下步骤：
a、实地测量主河与支沟交汇处主河的平均宽度B1和主河与支沟交汇处支沟的平均宽度B2，通过式1计算主河单宽流量Q1，通过式2计算支沟单宽流量Q2，通过式3计算支沟与主河的单宽流量比ε；
Figure BDA0003467049860000031
Figure BDA0003467049860000032
Figure BDA0003467049860000033
其中，Q1为主河单宽流量，m2/s；Qz为主河流量，m3/s，通过调查获取；B1为主河与支沟交汇处主河的平均宽度，m；Q2为支沟单宽流量m2/s；Qn为稀性泥石流流量，m3/s，通过调查获取；B2为主河与支沟交汇处支沟的平均宽度，m；ε为支沟与主河的单宽流量比；
b、计算稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径Dm，按照式4进行计算；
Dm＝3.8Q1 3/h3.5        式4
其中，Dm为稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径，mm；Q1为主河单宽流量，m2/s；h为主河与支沟交汇处主河的水深，m，实地调查获取；
c、计算体积因子yt，体积因子yt按照式5进行计算；
yt＝V/(hB1B2)  式5
其中，yt为体积因子；V为稀性泥石流冲入河道中不被主河水流冲走的固体物质总量，m3，根据计算的稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径Dm以上部分在泥石流中泥沙的百分比确定；h为主河与支沟交汇处主河的水深，m，实地调查获取；B1为主河与支沟交汇处主河的平均宽度，m；B2为主河与支沟交汇处支沟的平均宽度，m；
d、通过式6计算入汇角因子
Figure BDA0003467049860000041
Figure BDA0003467049860000042
其中，
Figure BDA0003467049860000043
为入汇角因子；π为圆周率；θ为支流与主流的交汇角，度，0°＜θ＜180°，实地调查或卫星影像图获取；
e、通过式7建立稀性泥石流堵江判断指标y；
Figure BDA0003467049860000044
其中，y为稀性泥石流堵江判断指标；yt为体积因子；ε为支沟与主河的单宽流量比；D50为泥石流中值粒径，mm；D0为粗颗粒粒径，D0＝2mm；
Figure BDA0003467049860000045
为入汇角因子；
当y＜7.5时，判断为不发生堵江；
当7.5≤y＜13.7时，判断为发生局部堵江；
当13.7≤y＜24.1时，判断为发生大部堵江；
当y≥24.1时，判断为发生全部堵江。
所述步骤e中，局部堵江具体是指流体入河后出露部分的长度与交汇处河宽的比值为30％-69％。
所述步骤e中，大部堵江具体是指流体入河后出露部分的长度与交汇处河宽的比值为70％-99％。
所述步骤e中，全部堵江具体是指流体入河后出露部分的长度与交汇处河宽的比值为100％，且造成断流至少3分钟。
本发明所述泥石流中值粒径是指泥石流颗粒中重量百分数为50％的颗粒粒径。
本发明所述稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径Dm以上部分在泥石流中泥沙的百分比通过泥石流中的泥沙颗粒级配分布获取。
本发明的基本原理如下：
单宽流量比是稀性泥石流堵江的一个重大影响因素。
主河流量越大，越不易发生堵塞；稀性泥石流流量越大，越容易发生堵塞。单宽流量比很好的将流量、河道宽度、水深这些影响因素联系到了一起，合理简化了复杂的实际情况。式3得到的单宽流量比，反映了稀性泥石流流量与主河流量之间的潜在关系。
稀性泥石流中的粗颗粒也是堵江是否会发生的一大影响因素。如果稀性泥石流中的颗粒太细，很容易被主河水流冲走，不易发生堵江。稀性泥石流颗粒越粗，主河水动力越弱，不被主河水流冲走的最小颗粒粒径以上部分的粗颗粒中值粒径越大，粗颗粒越多，那么同样条件下发生堵江的可能性就越高。式4反映了稀性粗颗粒抵抗河流冲刷的粒径。
稀性泥石流的入河体积量，即不被主河冲走的实际体积和沟道交汇处的宽度、泥石流沟道宽度及主河水深也较大的影响了稀性泥石流堵江的结果。在稀性泥石流颗粒粒径和单宽流量比变化不大的情况下，稀性泥石流入河实际体积越大，沟道越窄，水越浅，就越容易发生全堵。式5很好的联系了稀性泥石流入河体积量与沟道及水深。
最后，支沟与主河的交汇角会较大幅度的影响稀性泥石流流体在沟口的冲淤状况，实验与文献调查可知钝角交汇比锐角交汇更容易发生堵江，式6描述了交汇角对堵江的影响。
本发明的有益效果主要表现在以下方面：
1、本发明，“a、实地测量主河与支沟交汇处主河的平均宽度B1和主河与支沟交汇处支沟的平均宽度B2，通过式1计算主河单宽流量Q1，通过式2计算支沟单宽流量Q2，通过式3计算支沟与主河的单宽流量比ε；b、计算稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径Dm，按照式4进行计算；c、计算体积因子yt，体积因子yt按照式5进行计算；d、通过式6计算入汇角因子
Figure BDA0003467049860000051
e、通过式7建立稀性泥石流堵江判断指标y”，结合室内试验和野外数据分析，结合稀性泥石流特点与堵江过程，将稀性泥石流堵江的影响因素进行无量纲化，能适应野外复杂多变的灾害情况，为稀性泥石流堵江防治提供科学指导，泥石流堵江判断准确，普适性好。
2、本发明，使用的参数均经过无量纲处理，能够较好的避免因为单位数量级差距带来的误差，从而能够较好的适用于野外情况。
3、本发明，具有量纲和谐，计算结果准确度高，普适性好的特点。
4、本发明，充分考虑了稀性泥石流的特点，在稀性泥石流进入主河与主河水流的相互作用过程中，会有一部分泥沙被主河水流带走造成实际堵塞主河的泥沙量减少，引入稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径Dm，使得获得的稀性泥石流冲入河道中不被主河水流冲走的固体物质总量V更准确，从而能够更好的判断堵江程度。
具体实施方式
实施例1
一种稀性泥石流堵江判断方法，包括以下步骤：
a、实地测量主河与支沟交汇处主河的平均宽度B1和主河与支沟交汇处支沟的平均宽度B2，通过式1计算主河单宽流量Q1，通过式2计算支沟单宽流量Q2，通过式3计算支沟与主河的单宽流量比ε；
Figure BDA0003467049860000061
Figure BDA0003467049860000062
Figure BDA0003467049860000063
其中，Q1为主河单宽流量，m2/s；Qz为主河流量，m3/s，通过调查获取；B1为主河与支沟交汇处主河的平均宽度，m；Q2为支沟单宽流量m2/s；Qn为稀性泥石流流量，m3/s，通过调查获取；B2为主河与支沟交汇处支沟的平均宽度，m；ε为支沟与主河的单宽流量比；
b、计算稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径Dm，按照式4进行计算；
Dm＝3.8Q1 3/h3.5        式4
其中，Dm为稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径，mm；Q1为主河单宽流量，m2/s；h为主河与支沟交汇处主河的水深，m，实地调查获取；
c、计算体积因子yt，体积因子yt按照式5进行计算；
yt＝V/(hB1B2)  式5
其中，yt为体积因子；V为稀性泥石流冲入河道中不被主河水流冲走的固体物质总量，m3，根据计算的稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径Dm以上部分在泥石流中泥沙的百分比确定；h为主河与支沟交汇处主河的水深，m，实地调查获取；B1为主河与支沟交汇处主河的平均宽度，m；B2为主河与支沟交汇处支沟的平均宽度，m；
d、通过式6计算入汇角因子
Figure BDA0003467049860000071
Figure BDA0003467049860000072
其中，
Figure BDA0003467049860000073
为入汇角因子；π为圆周率；θ为支流与主流的交汇角，度，0°＜θ＜180°，实地调查或卫星影像图获取；
e、通过式7建立稀性泥石流堵江判断指标y；
Figure BDA0003467049860000074
其中，y为稀性泥石流堵江判断指标；yt为体积因子；ε为支沟与主河的单宽流量比；D50为泥石流中值粒径，mm；D0为粗颗粒粒径，D0＝2mm；
Figure BDA0003467049860000075
为入汇角因子；
当y＜7.5时，判断为不发生堵江；
当7.5≤y＜13.7时，判断为发生局部堵江；
当13.7≤y＜24.1时，判断为发生大部堵江；
当y≥24.1时，判断为发生全部堵江。
所述步骤e中，局部堵江具体是指流体入河后出露部分的长度与交汇处河宽的比值为30％-69％。
所述步骤e中，大部堵江具体是指流体入河后出露部分的长度与交汇处河宽的比值为70％-99％。
所述步骤e中，全部堵江具体是指流体入河后出露部分的长度与交汇处河宽的比值为100％，且造成断流至少3分钟。
本实施例为最基本的实施方式，结合室内试验和野外数据分析，结合稀性泥石流特点与堵江过程，将稀性泥石流堵江的影响因素进行无量纲化，能适应野外复杂多变的灾害情况，为稀性泥石流堵江防治提供科学指导，泥石流堵江判断准确，普适性好。
实施例2
一种稀性泥石流堵江判断方法，包括以下步骤：
a、实地测量主河与支沟交汇处主河的平均宽度B1和主河与支沟交汇处支沟的平均宽度B2，通过式1计算主河单宽流量Q1，通过式2计算支沟单宽流量Q2，通过式3计算支沟与主河的单宽流量比ε；
Figure BDA0003467049860000081
Figure BDA0003467049860000082
Figure BDA0003467049860000083
其中，Q1为主河单宽流量，m2/s；Qz为主河流量，m3/s，通过调查获取；B1为主河与支沟交汇处主河的平均宽度，m；Q2为支沟单宽流量m2/s；Qn为稀性泥石流流量，m3/s，通过调查获取；B2为主河与支沟交汇处支沟的平均宽度，m；ε为支沟与主河的单宽流量比；
b、计算稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径Dm，按照式4进行计算；
Dm＝3.8Q1 3/h3.5        式4
其中，Dm为稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径，mm；Q1为主河单宽流量，m2/s；h为主河与支沟交汇处主河的水深，m，实地调查获取；
c、计算体积因子yt，体积因子yt按照式5进行计算；
yt＝V/(hB1B2)  式5
其中，yt为体积因子；V为稀性泥石流冲入河道中不被主河水流冲走的固体物质总量，m3，根据计算的稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径Dm以上部分在泥石流中泥沙的百分比确定；h为主河与支沟交汇处主河的水深，m，实地调查获取；B1为主河与支沟交汇处主河的平均宽度，m；B2为主河与支沟交汇处支沟的平均宽度，m；
d、通过式6计算入汇角因子
Figure BDA0003467049860000091
Figure BDA0003467049860000092
其中，
Figure BDA0003467049860000093
为入汇角因子；π为圆周率；θ为支流与主流的交汇角，度，0°＜θ＜180°，实地调查或卫星影像图获取；
e、通过式7建立稀性泥石流堵江判断指标y；
Figure BDA0003467049860000094
其中，y为稀性泥石流堵江判断指标；yt为体积因子；ε为支沟与主河的单宽流量比；D50为泥石流中值粒径，mm；D0为粗颗粒粒径，D0＝2mm；
Figure BDA0003467049860000095
为入汇角因子；
当y＜7.5时，判断为不发生堵江；
当7.5≤y＜13.7时，判断为发生局部堵江；
当13.7≤y＜24.1时，判断为发生大部堵江；
当y≥24.1时，判断为发生全部堵江。
所述步骤e中，局部堵江具体是指流体入河后出露部分的长度与交汇处河宽的比值为30％-69％。
本实施例为较佳实施方式，使用的参数均经过无量纲处理，能够较好的避免因为单位数量级差距带来的误差，从而能够较好的适用于野外情况。
实施例3
一种稀性泥石流堵江判断方法，包括以下步骤：
a、实地测量主河与支沟交汇处主河的平均宽度B1和主河与支沟交汇处支沟的平均宽度B2，通过式1计算主河单宽流量Q1，通过式2计算支沟单宽流量Q2，通过式3计算支沟与主河的单宽流量比ε；
Figure BDA0003467049860000101
Figure BDA0003467049860000102
Figure BDA0003467049860000103
其中，Q1为主河单宽流量，m2/s；Qz为主河流量，m3/s，通过调查获取；B1为主河与支沟交汇处主河的平均宽度，m；Q2为支沟单宽流量m2/s；Qn为稀性泥石流流量，m3/s，通过调查获取；B2为主河与支沟交汇处支沟的平均宽度，m；ε为支沟与主河的单宽流量比；
b、计算稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径Dm，按照式4进行计算；
Dm＝3.8Q1 3/h3.5        式4
其中，Dm为稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径，mm；Q1为主河单宽流量，m2/s；h为主河与支沟交汇处主河的水深，m，实地调查获取；
c、计算体积因子yt，体积因子yt按照式5进行计算；
yt＝V/(hB1B2)  式5
其中，yt为体积因子；V为稀性泥石流冲入河道中不被主河水流冲走的固体物质总量，m3，根据计算的稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径Dm以上部分在泥石流中泥沙的百分比确定；h为主河与支沟交汇处主河的水深，m，实地调查获取；B1为主河与支沟交汇处主河的平均宽度，m；B2为主河与支沟交汇处支沟的平均宽度，m；
d、通过式6计算入汇角因子
Figure BDA0003467049860000111
Figure BDA0003467049860000112
其中，
Figure BDA0003467049860000113
为入汇角因子；π为圆周率；θ为支流与主流的交汇角，度，0°＜θ＜180°，实地调查或卫星影像图获取；
e、通过式7建立稀性泥石流堵江判断指标y；
Figure BDA0003467049860000114
其中，y为稀性泥石流堵江判断指标；yt为体积因子；ε为支沟与主河的单宽流量比；D50为泥石流中值粒径，mm；D0为粗颗粒粒径，D0＝2mm；
Figure BDA0003467049860000115
为入汇角因子；
当y＜7.5时，判断为不发生堵江；
当7.5≤y＜13.7时，判断为发生局部堵江；
当13.7≤y＜24.1时，判断为发生大部堵江；
当y≥24.1时，判断为发生全部堵江。
所述步骤e中，局部堵江具体是指流体入河后出露部分的长度与交汇处河宽的比值为30％-69％。
所述步骤e中，大部堵江具体是指流体入河后出露部分的长度与交汇处河宽的比值为70％-99％。
所述步骤e中，全部堵江具体是指流体入河后出露部分的长度与交汇处河宽的比值为100％，且造成断流至少3分钟。
本实施例为最佳实施方式，具有量纲和谐，计算结果准确度高，普适性好的特点。
充分考虑了稀性泥石流的特点，在稀性泥石流进入主河与主河水流的相互作用过程中，会有一部分泥沙被主河水流带走造成实际堵塞主河的泥沙量减少，引入稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径Dm，使得获得的稀性泥石流冲入河道中不被主河水流冲走的固体物质总量V更准确，从而能够更好的判断堵江程度。
下面对本发明的实施方式进行具体说明。
案例一：
1988年7月15日西藏波密县米堆沟上游光谢错发生冰碛湖溃决，溃决洪水侵蚀米堆沟内固体物质形成稀性泥石流，米堆沟冰湖溃决形成泥石流直倾帕隆藏布。
根据调查，米堆沟的稀性泥石流流量Qn为1021m3/s，主河与支沟交汇处支沟的平均宽度B2为20m，主河流量Qz为194m3/s，主河与支沟交汇处主河的平均宽度B1为60m，支流与主流的交汇角θ为113°，冲出量36.8万方。计算得到支主沟单宽流量比为15.79，体积因子yt为7.37，稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径Dm为43.04mm，泥石流中值粒径D50为7.5mm，入汇角因子
Figure BDA0003467049860000121
为1.97，最后算出的稀性泥石流堵江判断指标y为40.3。
因此米堆沟稀性泥石流发生全部堵江形成堰塞湖判断正确。
案例二：
西藏波密米堆沟泥石流堵河模型试验，山地学报，2005，2(3)：288-293)开展了米堆沟稀性泥石流堵江模型实验研究，分别研究了不同规模稀性泥石流进入主河与不同规模洪水遭遇，在交汇口表现出不同的堵溃特征。实验特征参数与结果见表1：
表1
Figure BDA0003467049860000122
从表1可知，采用本发明所述方法判断稀性泥石流堵江情况，完全正确的有3个，分别为实验组次1-3，占全部实验组次的75％；低估判断的有1个，为实验组次4，占25％，低估1个档次。总体上本发明判断这一系列实验基本正确。
综上所述，本发明针对各种情况下稀性泥石流暴发时计算是否发生堵江有较高的准确性。

Claims (4)

1.一种稀性泥石流堵江判断方法，其特征在于，包括以下步骤：
a、实地测量主河与支沟交汇处主河的平均宽度和主河与支沟交汇处支沟的平均宽度，通过式1计算主河单宽流量，通过式2计算支沟单宽流量，通过式3计算支沟与主河的单宽流量比；
Figure FDA0004105707860000011
Figure FDA0004105707860000012
Figure FDA0004105707860000013
其中，Q1为主河单宽流量，单位为m2/s；Qz为主河流量，单位为m3/s，通过调查获取；B1为主河与支沟交汇处主河的平均宽度，单位为m；Q2为支沟单宽流量，单位为m2/s；Qn为稀性泥石流流量，单位为m3/s，通过调查获取；B2为主河与支沟交汇处支沟的平均宽度，单位为m；ε为支沟与主河的单宽流量比；
b、计算稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径，按照式4进行计算；
Dm＝3.8Q1 3/h3.5        式4
其中，Dm为稀性泥石流冲入河道中不被主河水流冲走的泥石流最小颗粒粒径，单位为mm；h为主河与支沟交汇处主河的水深，单位为m，实地调查获取；
c、计算体积因子，体积因子按照式5进行计算；
yt＝V/(hB1B2)  式5
其中，yt为体积因子；V为稀性泥石流冲入河道中不被主河水流冲走的固体物质总量，单位为m3，通过选取大于和等于Dm的泥沙在泥石流中所占的百分比确定；
d、通过式6计算入汇角因子；
Figure FDA0004105707860000021
其中，
Figure FDA0004105707860000022
为入汇角因子；π为圆周率；θ为支流与主流的交汇角，单位为°，0°＜θ＜180°，实地调查或卫星影像图获取；
e、通过式7建立稀性泥石流堵江判断指标；
Figure FDA0004105707860000023
其中，y为稀性泥石流堵江判断指标；D50为泥石流中值粒径，单位为mm；D0为粗颗粒粒径，单位为mm；D0＝2mm；
当y＜7.5时，判断为不发生堵江；
当7.5≤y＜13.7时，判断为发生局部堵江；
当13.7≤y＜24.1时，判断为发生大部堵江；
当y≥24.1时，判断为发生全部堵江。
2.根据权利要求1所述的一种稀性泥石流堵江判断方法，其特征在于：所述步骤e中，局部堵江具体是指流体入河后出露部分的长度与交汇处河宽的比值为30％-69％。
3.根据权利要求1所述的一种稀性泥石流堵江判断方法，其特征在于：所述步骤e中，大部堵江具体是指流体入河后出露部分的长度与交汇处河宽的比值为70％-99％。
4.根据权利要求1所述的一种稀性泥石流堵江判断方法，其特征在于：所述步骤e中，全部堵江具体是指流体入河后出露部分的长度与交汇处河宽的比值为100％，且造成断流至少3分钟。
CN202210032435.9A 2022-01-12 2022-01-12 一种稀性泥石流堵江判断方法 Active CN114595426B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210032435.9A CN114595426B (zh) 2022-01-12 2022-01-12 一种稀性泥石流堵江判断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210032435.9A CN114595426B (zh) 2022-01-12 2022-01-12 一种稀性泥石流堵江判断方法

Publications (2)

Publication Number Publication Date
CN114595426A CN114595426A (zh) 2022-06-07
CN114595426B true CN114595426B (zh) 2023-04-07

Family

ID=81803897

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210032435.9A Active CN114595426B (zh) 2022-01-12 2022-01-12 一种稀性泥石流堵江判断方法

Country Status (1)

Country Link
CN (1) CN114595426B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106202771A (zh) * 2016-07-19 2016-12-07 成都理工大学 一种粘性泥石流堵江计算方法及其应用
CN106683019A (zh) * 2016-07-19 2017-05-17 成都理工大学 一种泥石流堵江计算方法及其应用
CN107169229A (zh) * 2017-06-09 2017-09-15 成都理工大学 一种计算粘性泥石流堵江的临界条件的方法
CN112729768A (zh) * 2021-01-20 2021-04-30 中国科学院、水利部成都山地灾害与环境研究所 一种泥石流入汇主河运动过程实验测量系统及测算方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106202771A (zh) * 2016-07-19 2016-12-07 成都理工大学 一种粘性泥石流堵江计算方法及其应用
CN106683019A (zh) * 2016-07-19 2017-05-17 成都理工大学 一种泥石流堵江计算方法及其应用
CN107169229A (zh) * 2017-06-09 2017-09-15 成都理工大学 一种计算粘性泥石流堵江的临界条件的方法
CN112729768A (zh) * 2021-01-20 2021-04-30 中国科学院、水利部成都山地灾害与环境研究所 一种泥石流入汇主河运动过程实验测量系统及测算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
白龙江流域泥石流危险性评价;李淑松;《中国水土保持科学》;第16卷(第4期);83-88 *

Also Published As

Publication number Publication date
CN114595426A (zh) 2022-06-07

Similar Documents

Publication Publication Date Title
Sumi et al. Reservoir sedimentation management with bypass tunnels in Japan
Chiew et al. Failure behavior of riprap layer at bridge piers under live-bed conditions
Garde et al. Mechanics of sediment transportation and alluvial stream problems
CN108984962B (zh) 溃决型黏性泥石流冲击下大颗粒堰塞坝溃决流量计算方法
CN105527405A (zh) 一种泥石流汇江物理模拟试验装置及试验方法
CN107423484B (zh) 一种泥石流在拦挡坝后的流量计算方法及其应用
CN110008608B (zh) 一种震区溃决型泥石流流量的计算方法及应用
CN102116014B (zh) 一种滑坡泥石流坝溃口特征测算方法及其应用
CN106202771B (zh) 一种粘性泥石流堵江计算方法及其应用
CN103343526B (zh) 一种窗口坝拦截泥石流闭塞类型判别方法及其应用
CN108166435B (zh) 格子坝拦挡粘性泥石流闭塞临界综合判断方法
CN110083898B (zh) 窗口坝调控粘性泥石流颗粒粒径评价方法
CN106683019B (zh) 一种泥石流堵江计算方法及其应用
CN114595426B (zh) 一种稀性泥石流堵江判断方法
CN111719492B (zh) 一种泥石流疏导系统及其疏导方法
CN110766270B (zh) 基于山区河流形态与主支流动量比变化的交汇区山洪泥沙灾害易灾区识别方法
Meng et al. Influence of internal structure on breaking process of short-lived landslide dams
CN107169229B (zh) 一种计算粘性泥石流堵江的临界条件的方法
CN105603940A (zh) 梯级鱼脊型水石分离系统及设计方法
CN112381319B (zh) 基于河床比降与水沙耦合水位的山洪灾害类型预测方法
Talchabhadel et al. Experimental study on suspended sediment transport to represent Tidal Basin Management
CN108797510B (zh) 一种泥石流堵江的早期识别方法
CN108755566B (zh) 泥石流拦砂坝坝基扬压力计算方法
CN110765650A (zh) 一种泥石流体积含沙量的测算方法
Lipiński et al. A study of granulation distribution in a beach of tailings dam

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