专利名称：一种泥石流预警方法
技术领域：
本发明涉及一种泥石流预警方法，特别是涉及一种通过对引发泥石流的降雨量阈 值的预置与监测而实现的泥石流预警方法，属于水利工程领域。
背景技术：
泥石流是一种由分选性很差的岩石、泥土、大量的水以及少量气体所构成的流体 沿着斜坡和沟谷向下流动的自然现象，其中最为常见的泥石流类型是由降雨触发的泥石 流，称为暴雨泥石流。对于暴雨型泥石流而言，泥石流沟流域内的降雨量与降雨强度的大小 是泥石流激发的决定性因素。当流域内降雨条件达到或超过某一量级的雨量与强度时，该 流域形成暴雨泥石流，这时的降雨量与降雨强度称为该条泥石流沟的雨量阈值。因此，雨量 阈值是暴雨型泥石流预报非常重要的报警指标，合理的雨量阈值指标是保障泥石流预警报 准确性的关键。目前，国内外利用降雨量与降雨强度作为泥石流雨量阈值的通行方法是通过对 预警区实际的降雨和泥石流灾害资料进行统计分析，分别计算出泥石流发生时的前期有效 降雨量与特征雨量的临界值，以临界值作为泥石流预警的雨量阈值。应用该方法对预警区 的数据资料要求较高，需要预警区有较丰富的、长期的雨量序列资料与灾害资料，特别是需 要有短历时降雨量(Ihr降雨量，30min降雨量，或IOmin降雨量等)历史数据才能分别计 算确定泥石流发生时的前期有效降雨量与特征雨量的临界值，由此才能绘制出雨量阈值曲 线。因此这种通行的做法只能适用于具有长期观测历史的地区。但在实际中，绝大多数的 泥石流流域没有长期的监测系统，非常缺乏降雨与灾害资料的长期观测数据。因而对这类 地区泥石流发生的预警工作由于缺乏基础数据而无法适用通行的方法，其泥石流预警的雨 量阈值计算确定问题也一直没有可以适用的方法。
发明内容
本发明的目的就是针对现有技术的不足，提供根据泥石流预警区的降雨和水文特 征及泥石流流域内地形、松散固体物质等数据资料，计算确定泥石流预警雨量阈值的泥石 流预警方法，该方法能够实现对灾害与雨量资料缺乏地区的泥石流发生的预警。为实现上述目的，本发明的技术方案如下 一种泥石流预警方法，通过在泥石流预警系统中预置预警雨量阈值，并将实时雨 量数据与预警雨量阈值实时比对结果判断是否启动系统预警报来实现，其特征在于所述 预警雨量阈值是雨量监测站点降雨量数据的16(1+&值，其中I6tl为Ihr雨强，Pa为前期影响 雨量。上述方法是一种基于水力类泥石流起动机理计算确定预警区泥石流预警雨量阈 值的方法，其关键技术在于对预警雨量阈值的选定。预警雨量阈值确定的基本原理在于由 于暴雨泥石流的暴发与雨季发生时间相同，因此在进行泥石流预报时，泥石流产汇流方式 可认为符合“蓄满产流”原理。根据“蓄满产流”原理，泥石流水量平衡方程表达式则为R =P-U = P-(Im-Pa),式中R为总径流深(mm) ；P为一次降雨量(mm) ;Pa为前期影响雨量(mm)； Im为降雨结束时流域达到的最大蓄水量(mm，对一特定流域，Im为常数)；U为一次降雨的损 失量(mm)。本发明以小时雨量为特征雨量进行泥石流预警，则R = P-U = P-(Iffl-Pa)式中 的一次降雨量P换为小时雨量I6tl，该式可以改写得到公式1 I60+Pa = R+Im 公式 1 由公式1可知，由于对一特定流域，Iffl为常数。所以，只要计算出径流深R，I60+Pa 为一定值。也即是说，当流域的I6(l+pa达到！？+^时，就表明该流域即将发生泥石流。本发明 即以I6(l+Pa值作为对预警区泥石流灾害预警的雨量阈值。在实际工程计算中，公式1中使 用到的各项参数都是整个流域的平均值，因此通过公式1计算得出的也是该流域I6(l+Pa的 平均值，可以表示为7^；，则公式1改写得到公式2 I60+Pa=R + lJ^:2式中，—流域平均雨量阈值(mm)。基于上述原理，本发明在预警区内，尤其是在泥石流形成区布设至少一个雨量监 测站点获取预警区降雨的实时观测资料，并依经验公式计算确定预警区泥石流流域总径流 深R，则能依据公式1计算确定预警区内泥石流灾害的流域平均雨量阈值。计算所得的流 域平均雨量阈值可以依据预警区的地形数据与建立的预警区特征雨量与海拔变化关系式 (即公式3)换算得出各雨量监测站点对应的预警雨量阈值。H = AXIa公式 3式中，H——海拔(m)；I——特征降雨量(讓)；A——系数；a——指数。为了使预警雨量阈值的计算更准确，本方法可以做如下优化(1)对预警区雨量 监测站点的布置位置加以优化；(2)在预警区有至少一次泥石流历史资料数据记录的情况 下，可以根据预警区历史资料数据计算出预警区泥石流灾害预警的雨量阈值的值与Pa κ值，并将该计算值与基于水力类泥石流起动机理计算法所得雨量阈值加以比对与综合分 析。与现有技术相比，本发明的有益效果是方法能够解决大多数泥石流发生地区由 于缺乏长期观测数据而无法通过分别计算确定前期有效降雨量与特征雨量的临界值来确 定泥石流沟流域内预警雨量阈值的技术问题，为该类地区通过预置降雨量阈值、监测降雨 过程的方法进行泥石流预警提供了有效的技术手段。
具体实施例方式下面结合附图，对本发明的优选实施例作进一步的描述。实施例一如
图1所示，用本发明方法对宁南县史家沟进行泥石流预警。宁南县位于四川省凉山彝族自治州东南部，地处北纬26° 50' 12〃 27° 18' 34"，东经102° 27' 44〃 102° 55' 09〃。县城坐落在泥石流堆积扇群上 部，西北、东北、东南三面环山而开口向西南。史家沟位于县城正后方，沟内堆积了大量的松散固体物质，遇到一定强度的降雨极易引发泥石流，对县城安全构成严重威胁。对史家沟(预警区)进行泥石流预警按照如下方案实施——步骤Si、预警区现场调查与雨量监测站点布置步骤S11、一般流域调查对史家沟进行一般流域调查，通过现场对地形的考察(必要时结合地形图)确定 如下数据资料通过调查测量确定流域面积A、流域最大高差ΔΗ、形成区界限、流通区界限；根 据调查测量数据依常规方法计算确定主沟平均纵坡降、形成区面积、清水区面积、流通区 面积、土壤颗粒组成。步骤S12、泥石流预警系统布置安装步骤S121、预警系统雨量监测站点布置在史家沟内布置预警系统雨量监测站点，监测站点布置方法为(1)所有监测站 点均布置在史家沟泥石流沟流域范围内；(2)泥石流沟的形成区内布设监测站点；(3)重要 集雨区布设监测站点；(4)泥石流沟口布置监测站点。本实施例中史家沟内雨量监测站点分别布置在流域内水井湾、官村子、预警站3 个地点。步骤S122、泥石流预警系统安装泥石流预警系统主要由控制中心与分布在预警区内的雨量监测站点组成。各雨量 监测站点获取雨量观测值并实时传输到控制中心，由控制中心根据各雨量监测站点的实时 雨量数据与预置的预警雨量阈值比对的结果判断是否启动泥石流灾害预警报。步骤S2、确定泥石流预警系统预置预警雨量阈值I6Q+Pa步骤S21、起动法计算确定预警区流域平均雨量阈值采用泥石流起动机理计算预警区流域平均雨量阈值7^；。步骤S211、计算预警区泥石流启动临界水深hQ依如下高桥保经验公式计算史家沟泥石流启动临界水深hQ
权利要求
1.一种泥石流预警方法，通过在泥石流预警系统中预置预警雨量阈值，并将实时雨量 数据与预警雨量阈值实时比对结果判断是否启动系统预警报来实现，其特征在于所述预 警雨量阈值是雨量监测站点降雨量数据的16(1+&值，其中I6tl为Ihr雨强，Pa前期影响雨量。
2.根据权利要求1所述的方法，其特征在于按如下步骤实施 步骤Si、预警区现场调查与雨量监测站点布置步骤S11、一般流域调查对预警区进行一般流域调查，确定基本地形数据、松散固体物质特征数据； 步骤S12、泥石流预警系统布置安装 步骤S121、预警系统雨量监测站点布置在预警区内布置预警系统雨量监测站点，实时获取的雨量观测值； 步骤S122、泥石流预警系统安装在预警区内安装泥石流预警系统，预警系统依据步骤S121所获取的雨量观测值与系 统预置的预警雨量阈值的比对结果判断是否启动泥石流灾害预警报； 步骤S2、确定泥石流预警系统预置预警雨量阈值I6Q+Pa 步骤S21、起动法计算确定预警区流域平均雨量阈值
3.根据权利要求2所述的方法，其特征在于所述步骤S121中预警系统雨量监测站点 布置方法为所有监测站点均布置在预警区范围内；泥石流沟的形成区内布设监测站点； 重要集雨区布设监测站点；泥石流沟口布置监测站点。
4.根据权利要求2或3所述的方法，其特征在于所述步骤S21按如下方案实施 步骤S21、起动法计算确定预警区流域平均雨量阈值
5.根据权利要求4所述的方法，其特征在于所述步骤S22中各监测站点预警 雨量阈值I6(l+Pa的计算按如下方案实施步骤S22、确定预警区各监测站点预警雨量阈值I6Q+Pa 步骤S221、预警区降雨特征分析根据预警区历史资料数据分析得到预警区特征雨量与海拔关系表达式 H = AXIa式中，H——海拔(m)，由实际测量得到； I——特征降雨量(mm)，由实测降雨资料得到； A——系数； a——指数；步骤S222、计算预警区各监测站点预警雨量阈值I6Q+Pa根据步骤S21所得流域平均雨量阈值7^；、步骤Sll调查计算所得形成区面积、汇流 区面积数据并参照步骤S221所得关系式计算确定各监测站点预警雨量阈值I6(l+Pa，具体计 算依下式联立进行。
6.根据权利要求2或3或5所述的方法，其特征在于在预警区有泥石流历史记录的情 况下，在步骤S22中根据预警区历史泥石流数据资料计算泥石流发生时的I6cik值与PaK值， 再将该结果用于验正或修正依式计算所得的各监测站点预警雨量阈值I6(l+Pa。
7.根据权利要求6所述方法，其特征在于所述验正或修正计算所得的各监测站点预 警雨量阈值I6(l+Pa的方法是首先将依式计算所得各监测站点预警雨量阈值I6(l+Pa绘制在 预警区各监测站点雨量阈值预报图上，预报图横坐标为I6tl，纵坐标为Pa，I6(l+pa值即为一条 阈值曲线；然后将根据预警区泥石流历史数据资料分析计算所得值与PaK值点汇在预 报图上并分析点与阈值曲线的位置关系，必要时对曲线加以修正。
全文摘要
本发明公开了一种泥石流预警方法。针对现有技术中无法解决缺乏基础数据资料地区泥石流预警系统的预警雨量阈值计算确定的技术问题，本发明提供了一种可以适用于缺乏基础数据资料地区泥石流预警的方法。本方法将雨量监测站点布置在泥石流流域的合理地点；以雨量监测站点降雨量数据的1hr雨强值与前期影响雨量值之和作为预警雨量阈值；再基于泥石流起动机理，利用泥石流预警区的降雨、水文特征、流域内地形、松散固体物质等数据计算确定各站点的预警雨量阈值，并将各阈值预置在预警系统中，通过实时雨量数据与阈值实时比对结果判断是否启动系统预警报。本发明方法原理可靠，计算过程科学简便，实施设备简单、安装方便，能够普遍适用于各类泥石流发生地区泥石流预警的需要。
文档编号G01W1/14GK102103787SQ201110046589
公开日2011年6月22日 申请日期2011年2月27日 优先权日2011年2月27日
发明者吕娟, 季宪军, 杨顺, 柳金峰, 欧国强, 潘华利, 王钧, 黄江成 申请人:中国科学院水利部成都山地灾害与环境研究所