CN105740618B - 用于沟床起动型泥石流预测的降雨分割方法 - Google Patents

用于沟床起动型泥石流预测的降雨分割方法 Download PDF

Info

Publication number
CN105740618B
CN105740618B CN201610061328.3A CN201610061328A CN105740618B CN 105740618 B CN105740618 B CN 105740618B CN 201610061328 A CN201610061328 A CN 201610061328A CN 105740618 B CN105740618 B CN 105740618B
Authority
CN
China
Prior art keywords
rainfall
value
less
prophase programming
critical value
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
CN201610061328.3A
Other languages
English (en)
Other versions
CN105740618A (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 CN201610061328.3A priority Critical patent/CN105740618B/zh
Publication of CN105740618A publication Critical patent/CN105740618A/zh
Application granted granted Critical
Publication of CN105740618B publication Critical patent/CN105740618B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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

  • Sewage (AREA)

Abstract

本发明属于泥石流防治工程、水利工程领域，涉及用于沟床起动型泥石流预测的降雨分割方法，根据该区域每小时的降雨数据，确定前1h的降雨值是否达到或超过降雨临界值R0：如果小于降雨临界值R0，则不计入前期降雨；如果大于或等于降雨临界值R0，就将该小时降雨量作为最早的前期降雨值R1，并把之后各小时的降雨值作为前期降雨值R2，R3，……，Rn，并逐小时叠加获得第n小时的前期降雨值和R；该段降雨量Rn+1；根据前期降雨值和R与该时段降雨量Rn+1，计算降雨综合值R*。该方法不受时间控制，计算前期降雨对土壤入渗率的影响更符合实际情况，并且不受地区限制。

Description

用于沟床起动型泥石流预测的降雨分割方法
技术领域
本发明属于泥石流防治工程、水利工程领域，涉及用于沟床起动型泥石流预测的降雨分割方法。
背景技术
沟道是一种坡度陡、遭受间歇性洪水冲刷的集水道。坡面径流产生后，顺坡方向随着降雨径流汇集量的不断加大，形成单宽流量与流速都很大的山洪。沟道侵蚀是指由坡面径流集中冲蚀土壤和母岩并切入地面形成较大沟壑的这种侵蚀形态。沟道内松散物质固体较丰富，当松散固体物质受气候、环境、水文影响产生被掀动或遭受揭底，则导致泥石流启动，此类泥石流称为沟道侵蚀型泥石流。沟道侵蚀型泥石流的特点在于沟道松散物质启动形成泥石流是依靠径流提供推移力进而发生和维持运动，各粒径固体颗粒启动概率远较泥河流泥沙启动概率为大，其启动时候沟道遭受水流强烈侵蚀，沟道松散物质是泥石流固体物质主要来源。
降雨是泥石流最主要的激发形式，沟床起动型泥石流的激发也是由降雨引起的。这种类型的泥石流激发都是由2个降雨因素共同的作用结果：前期降雨和短历时的强降雨。前期降雨使土体入渗率降低，同样的降雨强度就会产生更大的超渗产流，并产生更大的表面径流，以及产生更大的洪水，激发泥石流的可能性也会更大。因此前期降雨对泥石流的准确预报非常重要。
自然界的降雨过程非常复杂，且没有重复的降雨模式。如果把激发泥石流的短期强降雨之前的降雨都算前期降雨，那有可能所有的前期降雨都很大，失去了前期降雨用于判断土体入渗率以及激发泥石流所需要的降雨强度的目的和意义，所以正确地判断前期降雨量对于正确预测泥石流的发生非常重要。目前对前期降雨的分割方法主要有以下2中方法：
前期降雨修正及累加法：即将激发降雨之前的降雨按照时间长短，降雨量的系数随距离激发的时间的增加而减小，将修正了系数的降雨量累加得到的降雨总量作为前期降雨量。该方法有的考虑前期的时间为30天，有的为15天，也有的为7天，或5天，或3天。该方法的主要问题有：(1)随时间增长而衰减的系数没有一个科学的方法，是人为确定的；(2)时间长度也没有科学的依据来确定合适的时间，如果时间很长(如30天)，不仅计算量大，而且很长时间的降雨对后期的激发降雨没有什么作用；如果时间很短(如3天)，又可能忽略3天以前很大的降雨的作用。这种方法对于很复杂的降雨过程，特别是降雨与晴天交替发生时，不能很好的反映降雨对土壤的入渗率的影响。
降雨强度与持续时间法：(1)在台湾地区泥石流研究基础上的方法：1小时雨强大于4mm/h时，可以当作一次降雨的开始，在连续6h雨强都小于4mm/h或12h内累积雨量小于10mm时作为这次降雨的结束。这种方法来源于台湾地区的资料，一方面该区域的降雨量远大于我国西部地区(主要发育沟床起动型泥石流的区域)的降雨量，4mm/h的降雨强度可能偏大；另一方面台湾地区发生的泥石流多以台风引起的，降雨历时一般2-3天，激发的泥石流以浅层滑坡型沟谷泥石流为主，沟床起动型泥石流较少，不适用于本发明涉及沟床起动型泥石流。(2)在汶川地震强烈影响区域的泥石流研究基础上的方法：1小时雨强大于1mm/h时，可以当作一次降雨的开始，连续6h雨强都小于1mm/h可以作为这次降雨的结束。这种方法来源于强震区域有大量的崩塌和滑坡发生，且仅在地震后3年以内的时间段，对于地震影响较小的区域以及地震后较长的时间段，该方法不适用。
因此一个适用的前期降雨分割方法，对于预测泥石流的发生，降低由于泥石流造成的重大损失与人员伤亡是非常有用的。减轻、防止泥石流灾害，提出合理的前期降雨分割方法，需要从降雨的产流过程深入地研究。
发明内容
本发明的目的就是针对现有技术的不足，提供一种激发泥石流的前期降雨计算方法，该方法的核心是以降雨产流为基础，判断前期降雨量的方法。
为实现上述目的，本发明的技术方案如下：
用于沟床起动型泥石流预测的降雨分割方法：
布置传感器实时测降雨量，根据研究区域每小时的降雨数据，确定前1h的降雨值是否达到或超过降雨临界值R0：如果小于降雨临界值R0，则不计入前期降雨，即前期降雨量归0；如果大于或等于降雨临界值R0，就将该小时降雨量作为最早的前期降雨值R1，并把之后各小时的降雨值作为前期降雨值R2，R3，……，Rn，并逐小时叠加获得第n小时的前期降雨值和R：R＝R0+R1+R2+……+Rn；
时段降雨量指最后1小时段的降雨，即第n+1小时的降雨Rn+1；
前一个降雨结束后，之后的降雨重新开始计算降雨量和前期降雨量；
最后，根据前期降雨值和R与该时段降雨量Rn+1，计算降雨综合值R*：R*＝R+13.5Rn+1。
其中，所述的前一个降雨结束或前期降雨量归0的判断依据为：
研究区域的植被覆盖率为x，x为0-1之间的小数，降雨临界值为R0＝(1+x)mm；
(1)对于植被覆盖率x为0时，临界值R0＝1mm，连续6h降雨量都小于1mm/h或连续12h累积降雨量小于8mm；
(2)对于植被覆盖率x为1时，临界值R0＝2mm，连续6h降雨量都小于2mm/h或连续12h累积降雨量小于12mm；
(3)当植被覆盖率为x，此时临界值R0＝(1+x)mm，连续6h降雨量都小于(1+x)mm/h或连续12h累积降雨量小于(8+4x)mm。
发明提供的用于沟床起动型泥石流预测的降雨分割方法，用于沟床起动类型沟谷泥石流的前期降雨计算，该方法不受时间控制，只与降雨过程有关，也没有人为地控制前期降雨的衰减系数；因此计算前期降雨降雨对土壤入渗率的影响更符合实际情况；该方法是通过土壤的降雨入渗实验来的，因此不受地区限制，具有更广泛的应用范围；同时考虑了各种植被的影响，因此对于没有地震影响的区域，或受地震影响的区域在很长的时间范围内，该方法都可以使用。
具体实施方式
下面结合优选实施例对本发明的技术方案作进一步的描述。
表1是2011年9月6日至18日南江县赶场降雨观测站的逐小时降雨监测数据。该区域属于植被覆盖良好区域，x为1，因此临界值R0＝2mm。用本发明方法对南江县赶场实时降雨量计算激发泥石流的前期降雨值。表中R为前期降雨值(mm)，D为小于临界值R0＝2mm的时间(h)。由于实际数据中没有用到连续12h累积降雨量小于12mm的判断，采用连续6h雨强都小于2mm/h已经判断出降雨分割了，因此表中没有给出连续12h的较小降雨累加值。
表1 2011年南江县赶场降雨监测数据级前期降雨值
表1中有几次较大降雨过程：
(1)2011年9月6日5时的1h降雨量25mm，超过了临界降雨量2mm，作为一场降雨的开始，累加计算前期降雨量；2011年9月6日5时至19时的前期降雨量达到了101.3mm。但在6日14时至19时的6小时内，每个小时的降雨量都小雨临界降雨量2mm，因此该场降雨结束，前期降雨量为0。
(2)2011年9月7日1时的1h降雨量2.6mm，超过了临界降雨量2mm，作为一场降雨的开始，累加计算前期降雨量；2011年9月7日1时至7日16时的前期降雨量达到了16.9mm。但在7日11时至16时的6小时内，每个小时的降雨量都小雨临界降雨量2mm，因此该场降雨结束，前期降雨量为0。
(3)2011年9月10日23时的1h降雨量2.4mm，超过了临界降雨量2mm，作为一场降雨的开始，累加计算前期降雨量；2011年9月10日23时至11日16时的前期降雨量达到了45mm。但在11日11时至16时的6小时内，每个小时的降雨量都小雨临界降雨量2mm，因此该场降雨结束，前期降雨量为0。
(4)2011年9月11日21时的1h降雨量3mm，超过了临界降雨量2mm，作为一场降雨的开始，累加计算前期降雨量；2011年9月11日21时至12日14时的前期降雨量达到了34.3mm。但在12日9时至14时的6小时内，每个小时的降雨量都小雨临界降雨量2mm，因此该场降雨结束，前期降雨量为0。
(5)2011年9月13日4时的1h降雨量3.5mm，超过了临界降雨量2mm，作为一场降雨的开始，累加计算前期降雨量；2011年9月13日4时至13日21时的前期降雨量达到了30.1mm。但在13日16时至21时的6小时内，每个小时的降雨量都小雨临界降雨量2mm，因此该场降雨结束，前期降雨量为0。
(6)2011年9月16日23时的1h降雨量3.7mm，超过了临界降雨量2mm，作为一场降雨的开始，累加计算前期降雨量；2011年9月16日23时至18日23时的前期降雨量达到了314mm。但在18日18时至23时的6小时内，每个小时的降雨量都小雨临界降雨量2mm，因此该场降雨结束，前期降雨量为0。

Claims (1)

1.用于沟床起动型泥石流预测的降雨分割方法，其特征在于，包括以下过程：
布置传感器实时测降雨量，根据研究区域每小时的降雨数据，确定前1h的降雨值是否达到或超过降雨临界值R0：如果小于降雨临界值R0，则不计入前期降雨，即前期降雨量归0；如果大于或等于降雨临界值R0，就将该小时降雨量作为最早的前期降雨值R1，并把之后各小时的降雨值作为前期降雨值R2，R3，……，Rn，并逐小时叠加获得第n小时的前期降雨值和R：R＝R0+R1+R2+……+Rn；
时段降雨量指最后1小时段的降雨，即第n+1小时的降雨Rn+1；
前一个降雨结束后，之后的降雨重新开始计算降雨量和前期降雨量；
最后，根据前期降雨值和R与该时段降雨量Rn+1，计算降雨综合值R*：R*＝R+13.5Rn+1；
所述的前一个降雨结束或前期降雨量归0的判断依据为：
研究区域的植被覆盖率为x，x为0-1之间的小数，降雨临界值为R0＝(1+x)mm；
(1)对于植被覆盖率x为0时，临界值R0＝1mm，连续6h降雨量都小于1mm/h或连续12h累积降雨量小于8mm；
(2)对于植被覆盖率x为1时，临界值R0＝2mm，连续6h降雨量都小于2mm/h或连续12h累积降雨量小于12mm；
(3)当植被覆盖率为x，此时临界值R0＝(1+x)mm，连续6h降雨量都小于(1+x)mm/h或连续12h累积降雨量小于(8+4x)mm。
CN201610061328.3A 2016-01-28 2016-01-28 用于沟床起动型泥石流预测的降雨分割方法 Expired - Fee Related CN105740618B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610061328.3A CN105740618B (zh) 2016-01-28 2016-01-28 用于沟床起动型泥石流预测的降雨分割方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610061328.3A CN105740618B (zh) 2016-01-28 2016-01-28 用于沟床起动型泥石流预测的降雨分割方法

Publications (2)

Publication Number Publication Date
CN105740618A CN105740618A (zh) 2016-07-06
CN105740618B true CN105740618B (zh) 2018-10-26

Family

ID=56247931

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610061328.3A Expired - Fee Related CN105740618B (zh) 2016-01-28 2016-01-28 用于沟床起动型泥石流预测的降雨分割方法

Country Status (1)

Country Link
CN (1) CN105740618B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106383373B (zh) * 2016-12-10 2018-10-30 四创科技有限公司 基于实测与预报降雨量的等值线、等值面绘制、预警方法
CN107808357A (zh) * 2017-09-28 2018-03-16 成都理工大学 一种红层地区平推式岩质滑坡降雨分割方法及其应用

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103473892A (zh) * 2013-08-23 2013-12-25 成都理工大学 沟道起动型泥石流发生预警方法及应用
CN103472497A (zh) * 2013-08-23 2013-12-25 成都理工大学 沟道起动型泥石流暴发规模测算方法及应用
CN104318058A (zh) * 2014-09-25 2015-01-28 航天科工惯性技术有限公司 基于雨量监测的泥石流预警方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103473892A (zh) * 2013-08-23 2013-12-25 成都理工大学 沟道起动型泥石流发生预警方法及应用
CN103472497A (zh) * 2013-08-23 2013-12-25 成都理工大学 沟道起动型泥石流暴发规模测算方法及应用
CN104318058A (zh) * 2014-09-25 2015-01-28 航天科工惯性技术有限公司 基于雨量监测的泥石流预警方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Preliminary study of the effect of earthquakes on the rainfall threshold of debris flows;Bin Yu et al.;《Engineering Geology》;20140426;第182卷;第130-135页 *
汶川震区暴雨泥石流发生的降雨阈值;周伟 等;《水科学进展》;20131130;第24卷(第6期);第786-793页 *
沟床起动型泥石流的10min降雨预报模型;余斌 等;《水科学进展》;20150531;第26卷(第3期);第347-355页 *
泥石流预报中前期有效降水量的确定;韦方强 等;《山地学报》;20050731;第23卷(第4期);第453-457页 *

Also Published As

Publication number Publication date
CN105740618A (zh) 2016-07-06

Similar Documents

Publication Publication Date Title
Berti et al. Field observations of a debris flow event in the Dolomites
Tang et al. Triggering conditions and depositional characteristics of a disastrous debris flow event in Zhouqu city, Gansu Province, northwestern China
Turowski et al. The impact of exceptional events on erosion, bedload transport and channel stability in a step‐pool channel
Overduin et al. Coastal changes in the Arctic
Navratil et al. High-frequency monitoring of debris-flow propagation along the Réal Torrent, Southern French Prealps
Kasai et al. Impacts of land use change on patterns of sediment flux in Weraamaia catchment, New Zealand
Burkard et al. Patterns and controls of gully growth along the shoreline of Lake Huron
Evans et al. Post‐mining landform evolution modelling: 2. Effects of vegetation and surface ripping
Bovis et al. The July 29, 1998, debris flow and landslide dam at Capricorn Creek, Mount Meager Volcanic Complex, southern Coast Mountains, British Columbia
Sanders et al. The sediment budget of an alpine cirque
Amiel et al. Characterization of the hydrogeology of the sacred Gihon Spring, Jerusalem: a deteriorating urban karst spring
Nott Tropical cyclones and the evolution of the sedimentary coast of northern Australia
David-Novak et al. Modern extreme storms and the rainfall thresholds for initiating debris flows on the hyperarid western escarpment of the Dead Sea, Israel
Church Mountains and montane channels
Coe Sediment production and delivery from forest roads in the Sierra Nevada, California
CN105740618B (zh) 用于沟床起动型泥石流预测的降雨分割方法
Lavigne Rate of sediment yield following small‐scale volcanic eruptions: a quantitative assessment at the Merapi and Semeru stratovolcanoes, Java, Indonesia
Major et al. Sediment erosion and delivery from Toutle River basin after the 1980 eruption of Mount St. Helens: A 30-year perspective
Waythomas et al. The geomorphology of an Aleutian volcano following a major eruption: the 7–8 August 2008 eruption of Kasatochi Volcano, Alaska, and its aftermath
Pelpola et al. Long-term bed load transport rate based on aerial-photo and ground penetrating radar surveys of fan-delta growth, Coast Mountains, British Columbia
Jueyi et al. Changes in sediment transport in the Kuye River in the Loess Plateau in China
Turcotte et al. Icing and aufeis in cold regions I: the origin of overflow
Mikuś et al. Long-term monitoring of the recruitment and dynamics of large wood in Kamienica Stream, Polish Carpathians
Kehew et al. Reconstructing an extreme flood from boulder transport and rainfall–runoff modelling: Wadi Isla, South Sinai, Egypt
McNamara et al. Case study of a large summer flood on the North Slope of Alaska: bedload transport

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into 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: 20181026

Termination date: 20220128