一种判断沟道泥石流当前是否发生及其未来会否发生的方法转让专利

申请号 : CN201310300656.0

文献号 : CN103413401B

文献日 :

基本信息:

PDF:

法律信息:

相似专利:

发明人 : 姜元俊乔建平姜震王萌黄栋杨宗佶田宏岭孟华君

申请人 : 中国科学院、水利部成都山地灾害与环境研究所

摘要 :

本发明公开了一种判断沟道泥石流当前是否发生及其未来会否发生的方法,该判断沟道泥石流当前是否发生的方法,包括:S1:选定沟道,利用计算hs;S2:利用h0=hdry-hwet计算h0;S3:比较h0是否大于或等于hs,是,发布泥石流启动警报信息,否,不发布警报信息。判断沟道泥石流未来会否发生的方法,包括:S1:选定上述沟道,利用Q0′=P0A0计算Q0′;S2:利用Q0=b0h0V0计算Q0;S3:利用计算Cq,并利用计算n;S4:利用气象部门发布的预测降雨量信息,分别得到P1、P2、P3……Pi;S5:利用Qi=CqPiA0计算Qi;S6:求解S7:比较Pi是否大于或等于P0,并利用迭代法确定hi;S8:比较hi是否大于或等于hs,是,发布泥石流启动预报信息,否,不发布预报信息。本发明流程简洁、判断结论准确。

权利要求 :

1.一种判断沟道泥石流当前是否发生的方法,其特征在于,包括以下步骤:S1:选定沟道,并利用公式 计算出泥

石流启动时该沟道内的水位高度hs,其中,θ为以角度单位表示的沟道的沟床坡度,dm和分别为沟道堆积物的砂砾平均粒径和内摩擦角,C*为堆积体的体积浓度,ρs为砂砾密度,ρw为水的密度;

S2:利用公式h0=hdry-hwet计算出选定的沟道当前时刻的水位高度h0,其中,hdry为该沟道内无水流时的高度,hwet为该沟道内形成汇流时的高度;

S3:比较h0是否大于或等于hs,是,则判定在当前时刻,泥石流已经发生,并发布泥石流启动的警报信息,否,则不需发布泥石流启动的警报信息。

2.根据权利要求1所述的一种判断沟道泥石流当前是否发生的方法,其特征在于,所述步骤S1中,θ的大小通过罗盘或测绳测得。

3.根据权利要求2所述的一种判断沟道泥石流当前是否发生的方法,其特征在于,所述步骤S1中,dm通过筛分仪测得, 的大小则通过三轴仪测得。

4.根据权利要求1~3任意一条所述的一种判断沟道泥石流当前是否发生的方法,其特征在于,所述步骤S2中,hdry和hwet均通过泥位计测得。

5.一种判断沟道泥石流未来会否发生的方法,其特征在于,包括以下步骤:S1:选定如权利要求1所述的沟道,并利用公式Q0′=P0A0计算出该沟道的断面在当前时刻的理论单位小时水流量Q0′,其中,P0为沟道在当前时刻的降雨量,A0为沟道断面的汇水面积;

S2:利用公式Q0=b0h0V0计算出选定的沟道的断面在当前时刻的实际单位小时水流量Q0,其中,b0为沟道的宽度,h0为权利要求1所述的当前时刻沟道的水位高度,V0为当前时刻沟道断面内的流水速度;

S3:根据步骤S1和S2,利用公式 计算出选定的沟道的断面在当前时刻的单位小时水流量的修正系数Cq,并利用公式 计算出泥石流沟槽的曼宁粗糙系数n,其中,R0为当前时刻沟道泥石流沟槽的水力半径,且R0=b0h0/(b0+2h0),S则为以弧度单位表示的沟道的沟床坡度;

S4:根据选定的沟道,利用气象部门发布的未来时刻预测降雨量信息,分别得到该沟道从当前时刻起未来第一小时的降雨量P1、未来第二小时的降雨量P2、未来第三小时的降雨量P3……未来第i小时的降雨量Pi,其中,i为大于0的自然数;

S5:利用公式Qi=CqPiA0计算出选定的沟道从当前时刻起未来第i小时的水流量Qi;

S6:根据步骤S2与S3得到如下公式:

Qi=hib0Vi (1)

Ri=b0hi/(b0+2hi) (3)

结合公式(1)、(2)、(3)得到 其中,

hi为沟道从当前时刻起未来第i小时的水位高度,Vi为沟道从当前时刻起未来第i小时其沟床上的平均流速,Ri为沟道从当前时刻起未来第i小时其泥石流沟槽的水力半径,Qi的值由步骤S5计算得到;

S7:比较Pi是否大于或等于P0,并令hi=h0,根据P0与Pi的大小关系,结合步骤S6,并利用迭代法最终确定沟道从当前时刻起未来第i小时的水位高度hi;

S8:比较hi是否大于或等于hs,是,则判定从当前时刻起未来第i小时,泥石流将会发生,并发布泥石流启动的预报信息,否,则不需发布泥石流启动的预报信息。

6.根据权利要求5所述的一种判断沟道泥石流未来会否发生的方法,其特征在于,所述步骤S1中,P0的大小通过当地的雨量站实时测得。

7.根据权利要求6所述的一种判断沟道泥石流未来会否发生的方法,其特征在于,所述步骤S2中,b0通过测距仪测得,V0则通过雷达水面波速仪测得。

8.根据权利要求5~7任意一条所述的一种判断沟道泥石流未来会否发生的方法,其特征在于,所述步骤S4包括以下步骤:S401:根据选定的沟道,利用气象部门发布的未来时刻预测降雨量信息,分别得到该沟道从当前时刻起未来一小时的降雨量P1s、未来两小时的降雨量P2s、未来三小时的降雨量P3s……未来i小时的降雨量Pis;

S402:分别利用公式P1=P1s、P2=P2s-P1s、P3=P3s-P2s……Pi-1=P(i-1)s-P(i-2)s和Pi=Pis-P(i-1)s计算出沟道从当前时刻起未来第一小时的降雨量P1、未来第二小时的降雨量P2、未来第三小时的降雨量P3……未来第i小时的降雨量Pi。

9.根据权利要求8所述的一种判断沟道泥石流未来会否发生的方法,其特征在于,所述步骤S4中,i=3。

10.根据权利要求9所述的一种判断沟道泥石流未来会否发生的方法,其特征在于,所述步骤S7的具体过程为:令hi=h0,并比较Pi是否大于或等于P0,是,则做出循环:hi=hi+0.00001h0,直到 此时得到的hi即确定为沟道从当前时刻起未来第i小时的水位高度;否,则做出循环:hi=hi-0.00001h0,直到此时得到的hi即确定为沟道从当前时刻起未来第i小时的水位高度。

说明书 :

一种判断沟道泥石流当前是否发生及其未来会否发生的方

技术领域

[0001] 本发明涉及一种判断沟道泥石流发生情况的方法,具体地说,是涉及一种判断沟道泥石流当前是否发生及其未来会否发生的方法。

背景技术

[0002] 泥石流作为一种重大的自然灾害,严重影响了山区经济的发展,并威胁着人民的生命和财产安全,同时,它也严重威胁着社会的可持续发展。
[0003] 目前,泥石流预警作为一项重要的减灾手段,已逐步受到减灾科技工作者和政府相关管理部门的广泛关注。现有的泥石流预警多采用降雨量作为判断泥石流是否发生的指标,然而,对于降雨量与泥石流启动之间的直接关系,目前还没有完善的理论和方法进行确定,大多依赖于个人经验来确定其所研究的泥石流沟道对应的泥石流启动降雨量,然后根据该泥石流启动降雨量数据对泥石流当前是否已经发生和未来是否会发生进行判断。显而易见,这些判断方式并不可靠,其与实际情况相比,普遍存在着较大的误差,很多时候都未能准确判断出泥石流的启动情况,这不仅严重影响了泥石流启动预判的准确性,而且也造成了减灾成本的浪费。

发明内容

[0004] 本发明的目的在于提供一种判断沟道泥石流当前是否发生及其未来会否发生的方法,主要解决现有的技术手段对沟道泥石流当前是否发生及未来会否发生的判定结论与实际情况之间存在误差大的问题。
[0005] 为了实现上述目的,本发明采用的技术方案如下:
[0006] 一种判断沟道泥石流当前是否发生的方法,包括以下步骤:
[0007] S1:选定沟道,并利用公式 计算出泥石流启动时该沟道内的水位高度hs,其中,θ为以角度单位表示的沟道的沟床坡度,dm和 分别为沟道堆积物的砂砾平均粒径和内摩擦角,C*为堆积体的体积浓度,ρs为砂砾密度,ρw为水的密度;
[0008] S2:利用公式h0=hdry-hwet计算出选定的沟道当前时刻的水位高度h0,其中,hdry为该沟道内无水流时的高度,hwet为该沟道内形成汇流时的高度;
[0009] S3:比较h0是否大于或等于hs,是,则判定在当前时刻,泥石流已经发生,并发布泥石流启动的警报信息,否,则不需发布泥石流启动的警报信息。
[0010] 作为优选,所述步骤S1中,θ的大小通过罗盘或测绳测得。
[0011] 作为优选,所述步骤S1中,dm通过筛分仪测得, 的大小则通过三轴仪测得。
[0012] 作为优选,所述步骤S2中,hdry和hwet均通过泥位计测得。
[0013] 基于上述基础,本发明还提供了一种判断沟道泥石流未来会否发生的方法,包括以下步骤:
[0014] S1:选定如上所述的沟道,并利用公式Q0′=P0A0计算出该沟道的断面在当前时刻的理论单位小时水流量Q0′,其中,P0为沟道在当前时刻的降雨量,A0为沟道断面的汇水面积;
[0015] S2:利用公式Q0=b0h0V0计算出选定的沟道的断面在当前时刻的实际单位小时水流量Q0,其中,b0为沟道的宽度,h0为当前时刻沟道的水位高度,V0为当前时刻沟道断面内的流水速度;
[0016] S3:根据步骤S1和S2,利用公式 计算出选定的沟道的断面在当前时刻的单位小时水流量的修正系数Cq,并利用公式 计算出泥石流沟槽的
曼宁粗糙系数n,其中,R0为当前时刻沟道泥石流沟槽的水力半径,且R0=b0h0/(b0+2h0),S则为以弧度单位表示的沟道的沟床坡度;
[0017] S4:根据选定的沟道,利用气象部门发布的未来时刻预测降雨量信息,分别得到该沟道从当前时刻起未来第一小时的降雨量P1、未来第二小时的降雨量P2、未来第三小时的降雨量P3……未来第i小时的降雨量Pi,其中,i为大于0的自然数;
[0018] S5:利用公式Qi=CqPiA0计算出选定的沟道从当前时刻起未来第i小时的水流量Qi;
[0019] S6:根据步骤S2与S3得到如下公式:
[0020] Qi=hib0Vi (1)
[0021]
[0022] Ri=b0hi/(b0+2hi) (3)
[0023] 结合公式(1)、(2)、(3)得到其中,hi为沟道从当前时刻起未来第i小时的水位高度,Vi为沟道从当前时刻起未来第i小时其沟床上的平均流速,Ri为沟道从当前时刻起未来第i小时其泥石流沟槽的水力半径,Qi的值由步骤S5计算得到;
[0024] S7:比较Pi是否大于或等于P0,并令hi=h0,根据P0与Pi的大小关系,结合步骤S6,并利用迭代法最终确定沟道从当前时刻起未来第i小时的水位高度hi;
[0025] S8:比较hi是否大于或等于hs,是,则判定从当前时刻起未来第i小时,泥石流将会发生,并发布泥石流启动的预报信息,否,则不需发布泥石流启动的预报信息。
[0026] 进一步地,所述步骤S1中,P0的大小通过当地的雨量站实时测得。
[0027] 作为优选,所述步骤S2中,b0通过测距仪测得,V0则通过雷达水面波速仪测得。
[0028] 具体地说,所述步骤S4包括以下步骤:
[0029] S401:根据选定的沟道,利用气象部门发布的未来时刻预测降雨量信息,分别得到该沟道从当前时刻起未来一小时的降雨量P1s、未来两小时的降雨量P2s、未来三小时的降雨量P3s……未来i小时的降雨量Pis;
[0030] S402:分别利用公式P1=P1s、P2=P2s-P1s、P3=P3s-P2s……Pi-1=P(i-1)s-P(i-2)s和Pi=Pis-P(i-1)s计算出沟道从当前时刻起未来第一小时的降雨量P1、未来第二小时的降雨量P2、未来第三小时的降雨量P3……未来第i小时的降雨量Pi。
[0031] 作为优选,所述步骤S4中,i=3。
[0032] 再进一步地,所述步骤S7的具体过程为:令hi=h0,并比较Pi是否大于或等于P0,是,则做出循环:hi=hi+0.00001h0,直到 此时得到的hi即确定为沟道从当前时刻起未来第i小时的水位高度;否,则做出循环:hi=hi-0.00001h0,直到 此时得到的hi即确定为沟道从当前时刻起未来第i小时
的水位高度。
[0033] 与现有技术相比,本发明具有以下有益效果:
[0034] (1)本发明构思巧妙,符合实际,流程清晰明了,且操作便捷。
[0035] (2)本发明在大量的计算和实验的基础上推导出了公式并利用该公式与选定的沟道当前时刻的
实际水位高度h0进行比较,判断出该沟道当前时刻是否已经发生泥石流,进而决定是否需要发布相应的泥石流启动警报;本专利申请的发明人通过研究、计算及合理的推导后,提出了将沟道泥石流启动时的水位高度与沟道断面当前的实际水位高度结合作为判断泥石流当前是否启动的标准,该标准不仅具有设计合理、判断准确的特点,而且改变了常规的利用降雨量来判断泥石流启动的方式,克服了技术偏见,因此,本发明的技术方案是一个富有创造性的方案,无论是根据现有技术的单独对比原则,还是根据组合对比原则,甚至是组合对比并结合常规技术手段的原则,均不能显而易见地得出本发明的技术方案。
[0036] (3)由于是在野外环境进行测量和计算,其地形较为复杂,计算结果误差较大,因而本发明首先计算出选定的沟道的断面在当前时刻的理论单位小时水流量,然后再结合实际单位小时水流量得到修正系数Cq,利用该修正系数对理论单位小时水流量进行校正,从而便可得到符合实际情况且精度较高的水流量Q,然后根据气象部门发布的未来时刻预测降雨量信息得到的Pi,并结合公式 从而利用迭代法最终确定沟道从当前时刻起未来第i小时的水位高度hi,最后再将确定的hi与hs进行比较,即可准确判断出选定的沟道从当前时刻起未来第i小时时泥石流是否会发生,并确定是否需要发布相应的泥石流预报信息;本发明在基于事实的基础上准确预判了未来泥石流的发生情况,其预测方式新颖、切合实际、数据可靠、预测精度高,可大幅降低减灾成本和预算。
[0037] (4)本发明中所有的参数均可通过现有的测量设备、实验手段和合理的计算推导得到,因此,其具有数据获取方便、逻辑紧密的特点。
[0038] (5)本发明由从事地质灾害研究工作多年的高级工程师研究和设计,其在大量的实验基础上总结出了符合泥石流启动的规律,并使该规律完美地应用于实际,因此,本发明突破现有技术的瓶颈,实现了重大的创新,其具有很高的实用价值和推广价值。

具体实施方式

[0039] 下面结合实施例对本发明作进一步说明,本发明的实施方式包括但不限于下列实施例。
[0040] 实施例
[0041] 本发明主要用于判断沟道当前时刻泥石流是否已经发生及未来时刻会否发生,其中,用于判断沟道当前时刻泥石流是否发生的步骤如下:
[0042] S1:选定沟道,并利用公式 计算出泥石流启动时该沟道内的水位高度hs,其中,θ为以角度单位表示的沟道的沟床坡度,dm和 分别为沟道堆积物的砂砾平均粒径和内摩擦角,C*为堆积体的体积浓度,ρs为砂砾密度,ρw为水的密度;本步骤中的所有参数均可通过现有的设备或实验手段测得,例如沟床坡度θ可通过罗盘或者测绳测出其大小,沟道堆积物的砂砾平均粒径dm可通过筛分仪测量沟道堆积物的砂砾粒径并求出其平均值得到,沟道堆积物的内摩擦角 可通过三轴仪测出其大小,堆积体的体积浓度C*和砂砾密度ρs可以通过现有的室内密度实验测得,而水的密度ρw则为公知常数;
[0043] S2:利用公式h0=hdry-hwet计算出选定的沟道当前时刻的水位高度h0,其中,hdry为该沟道内无水流时的高度,hwet为该沟道内形成汇流时的高度;本步骤中的hdry和hwet均可通过安装在沟道上的泥位计测得;
[0044] S3:比较h0是否大于或等于hs,是,则判定在当前时刻,泥石流已经发生,并发布泥石流启动的警报信息,否,则不需发布泥石流启动的警报信息。
[0045] 上述步骤是用于判断沟道在当前时刻是否已经发生泥石流,并进行警报。下面利用已经算出的沟道当前时刻的水位高度h0来判断未来时刻是否会发生泥石流,并根据判断得出的结论,决定是否需要发布泥石流启动的预报,其具体的判断步骤如下:
[0046] S1:利用公式Q0′=P0A0计算出选定的沟道的断面在当前时刻的理论单位小时水流量Q0′,其中,P0为沟道在当前时刻的降雨量,A0为沟道断面的汇水面积;本步骤中的P0通过当地的雨量站实时记录并测得,A0则根据选定沟道的地形流域,并通过现有的技术手段测出;
[0047] S2:利用公式Q0=b0h0V0计算出选定的沟道的断面在当前时刻的实际单位小时水流量Q0,其中,b0为沟道的宽度,V0为当前时刻沟道断面内的流水速度;本步骤中的h0即为前面所述的当前时刻沟道的水位高度,沟道宽度b0可通过测距仪测得,而当前时刻的沟道断面内的流水速度V0则可通过安装在沟道上的雷达水面波速仪测得;
[0048] S3:根据步骤S1和S2,利用公式 计算出选定的沟道的断面在当前时刻的单位小时水流量的修正系数Cq,并利用公式 计算出泥石流沟槽的
曼宁粗糙系数n,其中,R0为当前时刻沟道泥石流沟槽的水力半径,且R0=b0h0/(b0+2h0),S则为以弧度单位表示的沟道的沟床坡度;本步骤中的公式 由曼宁公式
推导出,S与前面的θ表示的是同一个沟床坡度的大小,只是单位不相同,
一个以弧度单位表示,一个以角度单位表示,并且S也可通过罗盘或测绳测得其大小;
[0049] S4:根据选定的沟道,利用气象部门发布的未来时刻预测降雨量信息,分别得到该沟道从当前时刻起未来一小时的降雨量P1s、未来两小时的降雨量P2s、未来三小时的降雨量P3s……未来i小时的降雨量Pis;
[0050] S5:分别利用公式P1=P1s、P2=P2s-P1s、P3=P3s-P2s……Pi-1=P(i-1)s-P(i-2)s和Pi=Pis-P(i-1)s计算出沟道从当前时刻起未来第一小时的降雨量P1、未来第二小时的降雨量P2、未来第三小时的降雨量P3……未来第i小时的降雨量Pi;本步骤中,由于气象部门对降雨量的预测精度随着时间的增长会出现下降,因此,为了保证预测的精度,本实施例只利用了气象部门发布的未来3个小时的降雨量数据,从而对沟道自当前时刻起未来3个小时进行泥石流是否会启动的预报;
[0051] S6:利用公式Qi=CqPiA0计算出选定的沟道从当前时刻起未来第i小时的水流量Qi,由于前面步骤S5中只利用了未来3个小时的降雨量数据,因此,本步骤中的i=1、2或3;
[0052] S7:根据步骤S2与S3得到如下公式:
[0053] Qi=hib0Vi (1)
[0054]
[0055] Ri=b0hi/(b0+2hi) (3)
[0056] 结合公式(1)、(2)、(3)得到
[0057] 其中,hi为沟道从当前时刻起未来第i小时的水位高度,Vi为沟道从当前时刻起未来第i小时其沟床上的平均流速,Ri为沟道从当前时刻起未来第i小时其泥石流沟槽的水力半径,本步骤中Qi的值由步骤S6计算得到;
[0058] S8:利用迭代法,令hi=h0,并比较Pi是否大于或等于P0,是,则做出循环:hi=hi+0.00001h0,直到 此时得到的hi即确定为沟
道从当前时刻起未来第i小时的水位高度;否,则做出循环:hi=hi-0.00001h0,直到此时得到的hi即确定为沟道从当前时刻起未来第i小时的
水位高度;
[0059] S9:比较hi是否大于或等于hs,是,则判定从当前时刻起未来第i小时,泥石流将会发生,并发布泥石流启动的预报信息,否,则不需发布泥石流启动的预报信息。按照这些步骤,便可很好地判断选定的沟道从当前时刻起未来第i小时时是否会发生泥石流。
[0060] 本发明在申请专利之前以不公开技术方案的方式进行了多次泥石流启动警报和预报的实验,均取得了不俗的效果,为方便理解本发明,本实施例以四川省成都市彭州市玉石沟安装的泥位计在2013年6月19日早上8点实时记录的沟道水位变化数据为例,对本发明的技术方案进行详细说明。
[0061] 选定沟道,然后通过现场调查和室内土工实验,得到下列表1和表2的参数(本实施例中所有参数的单位均为SI国际标准单位):
[0062]
[0063] 表1
[0064]
[0065] 表2
[0066] 利用公式 计算泥石流启动时该沟道内的水位高度hs,将上述表1和表2的数据代入该公式中,并且根据公知常识,可以确定水的密度ρw为1.0g/cm3,数据代入后即可得到:
[0067]
[0068] 根据泥位计记录的数据可知,在2013年6月19日早上8点(即当前时刻)时,选定的沟道内无水流时的高度hdry为1.23m,选定的沟道内形成汇流时的高度hwet为2.41m,则根据公式h0=hdry-hwet就可得到选定的沟道当前的水位高度h0为1.18m。由于h0=1.18m>hs=0.297m,因此,可得出结论:选定的沟道在8点时已经发生泥石流,需要发布泥石流启动的警报信息。
[0069] 接下来判断选定的沟道从8点起未来3个小时(即9点、10点、11点)是否会发生泥石流。
[0070] 根据选定的沟道流域的地形测出该沟道断面的汇水面积A0为4.07×106m2,同时根据玉石沟雨量站的记录数据,在2013年6月19日早上8点时该地的降雨量P0为0.06m/h,然后利用公式Q0′=P0A0计算出选定的沟道的断面在当前时刻的理论单位小时水流6 6 3
量Q0′,将A0和P0的值代入,得到Q0′=P0A0=0.06×4.07×10 =0.2442×10 m/h=
3
67.8m/s。
[0071] 根据安装在沟道上方的雷达水面波速仪可以测得2013年6月19日早上8点时该沟道断面内的流水速度V0=1.79m/s,结合上述计算出的沟道当前时刻的水位高度h0=1.18m,以及测得的沟道宽度b0=37m,代入公式Q0=b0h0V0,即可得到选定的沟道的断面在
3
当前时刻的实际单位小时水流量Q0=b0h0V0=37×1.18×1.79=78.1514m /s。
[0072] 然后根据计算得到的Q0′和Q0,计算出选定的沟道的断面在当前时刻的单位[0073] 小时水流量的修正系数 接着利用公式:
[0074]
[0075] R0=b0h0/(b0+2h0) (2)
[0076] 结合公式(1)、(2),并将前面测得的b0、S、V0以及计算出的h0代入,即可计算出泥石流沟槽的曼宁粗糙系数n的值:
[0077]
[0078] 根据气象部门发布的未来预测降雨量信息,分别得到该沟道从当前时刻(8点)起未来一小时(9点)的降雨量P1s=0.07m、未来两小时(10点)的降雨量P2s=0.12m、未来三小时(11点)的降雨量P3s=0.15m。然后分别利用公式P1=P1s、P2=P2s-P1s和P3=P3s-P2s即可计算出从当前时刻起未来第一小时的降雨量P1=P1s=0.07m、未来第二小时的降雨量P2=P2s-P1s=0.05m和未来第三小时的降雨量P3=P3s-P2s=0.03m。
[0079] 进而可得到从当前时刻起未来第一、二、三小时沟道断面的水流量Q1、Q2和Q3:
[0080] Q1=CqP1A0=1.15×0.07×4.07×106=3.28×105m3/h=91m3/s;
[0081] Q2=CqP2A0=1.15×0.05×4.07×106=2.34×105m3/h=65m3/s;
[0082] Q3=CqP3A0=1.15×0.03×4.07×106=1.4×105m3/h=39m3/s。
[0083] 根据公式 同时令hi=h0=1.18m,并根据上述计算得到的Q1、Q2和Q3,利用迭代法分别计算沟道从当前时刻起未来第一、二、三小时的水位高度h1、h2和h3,其结果如下:
[0084] (1)从当前时刻(8点)起未来第一小时(9点)的水位高度h1:
[0085] 由 于P1=0.07m>P 0= 0.06m,因 此,有 循 环:hi=hi+0.00001h0,直 到得到h1=1.57m>hs=0.297m,所以判定从当前时刻起未来第一小时沟道内将会发生泥石流,需要发布泥石流启动的预报信息;
[0086] (2)从当前时刻(8点)起未来第二小时(10点)的水位高度h2:
[0087] 由 于P2=0.05m

hs=0.297m,所以判定从当前时刻起未来第二小时沟道内将会发生泥石流,需要发布泥石流启动的预报信息;

[0088] (3)从当前时刻(8点)起未来第三小时(11点)的水位高度h3:
[0089] 由 于P3=0.03m

hs=0.297m,所以判定从当前时刻起未来第三小时沟道内将会发生泥石流,需要发布泥石流启动的预报信息。

[0090] 综上所述,本次泥石流的预报结果为:
[0091]时间 是否会发生泥石流
未来第一小时 是
未来第二小时 是
未来第三小时 是
[0092] 表3
[0093] 根据该地雨量站当时的记录可知,在2013年6月19日早上9点、10点和11点时