本发明公开泥石流灾害险情动态预警方法。在流域调查与监测系统实时数据基础上，采用二维潜水波方程数值模拟泥石流沿程运动，在DEM栅格水平上预测泥石流淹没范围与淹没深度，通过叠加承灾体数据分析出承灾体遭受危害的程度、安全疏散线路等信息，及时向流域内人群发布。本发明还公开泥石流灾害精细化分级监测预警方法，通过全方位获取分析影响泥石流灾害要素的监测信息，将整个泥石流形成、发生、演进过程视为高度动态过程加以预测监测，通过灾前预测、临灾预警、发灾预警、险情动态预警四级逐步完成泥石流灾害的精细预警。与现有的一步、一指标监测预警技术相比，本发明方法充分结合泥石流灾害的动态性特征与规律，更为科学，减灾实务性更高。
1.泥石流灾害精细化分级监测预警方法，在实施泥石流灾害险情动态预警前先完成泥石流灾害灾前预测，其特征在于：首先进行流域调查并安装监测系统，在流域调查中，确定泥石流沟道内流通区、危险区范围，在流通区布置监测断面M；现场调查确定监测断面M处坡降IC、监测断面M宽度B、泥石流沟床糙率系数nC、监测断面M至泥石流沟道危险区的距离L、危险区内承灾体分布；通过GIS软件进行流域内各断面的基本地形地貌参数调查，提取流域内DEM数据，根据DEM数据确定计算栅格数量、各栅格位置坐标、各栅格面积，确定数值模拟计算采用的网格的面积a，确定形成区径流系数C，在形成区的末端设置监测断面N，调查确定监测断面N以上的流域面积A1、监测断面N宽度B1、监测断面N所在处坡降Ic1，根据现场调查取样与现场土力学实验确定形成区堆积物颗粒平均粒径dm；
其次安装监测系统，监测系统包括以无线通信联接的控制中心与监测传感器，监测传感器安装在监测断面M处，包括流速传感器、雨量传感器、监测断面M上方的超声波泥位传感器；
步骤A1、确定流域内泥石流形成的临界水量Q依式8、式9、式10计算泥石流灾害临界水量Q，Q＝qB1 式8
θ＝arctan (Ic1) 式10式中，Q—泥石流灾害形成临界水量，单位m3/s，q—泥石流起动单宽流量，单位m3/s，
B1—监测断面N宽度，单位m，流域调查确定，dm—形成区堆积物颗粒平均粒径，单位m，流域调查确定，θ—监测断面N纵坡坡度，单位°，
步骤A2、确定流域内泥石流形成临界降雨阈值iL依式11计算泥石流形成临界降雨阈值iL，式中，iL—泥石流形成临界降雨阈值，mm/h，Q—泥石流灾害形成临界水量，单位m3/s，步骤A1确定，C—形成区径流系数，流域调查确定，
A1—监测断面N以上流域面积，单位ha，流域调查确定；
从气象卫星云图数据中读取流域内未来最大可降雨量ik，若未来最大可降雨量ik≥临界雨量阈值iL，判断流域内出现泥石流发生可能性，发出预警信息，进入步骤D；
采用二维潜水波方程进行数值模拟，得到泥石流沿程运动方程；其中，运动平衡动量方程依式1、式2，质量守恒方程依式3：式中，g—重力加速度常数，
t—泥石流持续时间，单位s，根据各地《中小流域暴雨洪水计算手册》方法框算确定，u—泥石流在x方向的平均流动速度，单位m/s，流速传感器实时监测确定，h—泥石流流深监测值，单位m，超声波泥位传感器实时监测确定，v—泥石流在y方向的平均流动速度，单位m/s，流速传感器实时监测确定，i—一次降雨强度，单位mm/h，雨量传感器实时监测确定，根据泥石流沿程运动方程确定泥石流最终淤埋区域A、泥石流最终不会淤埋的区域A′、每一栅格内的颗粒数Ni,j、泥石流颗粒体积ΔV、每一栅格中泥石流在x、y方向的速度u、v、每一栅格的动态泥深h'、泥石流密度ρ；
步骤D2、确定每个栅格泥石流淤埋深度Dh依式4计算每个栅格(i,j)的泥石流淤埋深度Dh式中，Dh—每一栅格的泥石流淤埋深度，单位m，Ni,j—每一栅格内的颗粒数，步骤D1数值模拟计算确定，ΔV—泥石流颗粒体积，单位m3，步骤D1数值模拟计算确定，A—泥石流最终淤埋区域，单位m2，步骤D1数值模拟计算确定，步骤D3、确定最终淤埋区域A内每个栅格的泥石流最大动能De依式5计算最终淤埋区域A内每个栅格的泥石流最大动能DeDe＝a·maxt>0[(u2+v2)h′ρ] 式5式中，De—每个栅格的泥石流最大动能，单位J，a—步骤D1数值模拟计算采用的网格的面积，单位m2，流域调查确定，t—数值模拟时间，单位s，
u、v—每一栅格中泥石流在x、y方向的速度，单位m/s，步骤D1数值模拟计算确定，h'—每一栅格的动态泥深，单位m，步骤D1数值模拟计算确定，ρ—泥石流密度，单位kg/m3，步骤D1数值模拟计算确定；
步骤D4、确定最终淤埋区域A内承灾体危险等级系数D依式6计算最终淤埋区域A内承灾体危险等级系数DD＝De+Dh 式6
D—最终淤埋区域A内承灾体危险等级系数，De—每一栅格的泥石流最大动能，单位J，步骤D3确定，Dh—每一栅格的泥石流淤埋深度，单位m，步骤D2确定，根据D判断各承灾体危险程度：若0＜D≤0.4，判断承灾体处于低危险性，0.4＜D≤0.8，判断承灾体处于中危险性，若0.8＜D，判断判断承灾体处于高危险性；
确定灾害警示信息，包括疏散时间t′、承灾体危险等级D、潜在危害区域、逃生路线，具体：依式7计算疏散时间t′，
L—监测断面M至泥石流沟道危险区的距离，单位m，流域调查确定，Vc—泥石平均流流速，单位m/s，流速传感器实时监测确定，潜在危害区域是泥石流最终淤埋区域A，由步骤D1数值模拟计算确定，逃生路线是从泥石流最终淤埋区域A内各承灾体至泥石流最终不会淤埋的区域A′的最短路径，采用深度优先搜索方法确定；
2.根据权利要求1所述的泥石流灾害精细化分级监测预警方法，其特征在于：在步骤A3中判断流域内出现泥石流发生可能性后进入步骤B、泥石流临灾预警，首先进行流域调查并安装监测系统，在安装的监测系统中还包括在流域内设置海拔梯度雨量站；然后依如下步骤实施：步骤B、泥石流临灾预警
采集流域内不同海拔雨量站雨量观测值，组成流域梯度降雨观测数据，计算降雨观测平均值；
将降雨观测平均值与临界降雨阈值iL对比，若降雨观测平均值≥临界雨量阈值iL，判定流域内泥石流灾害临近，发出预警信息；进入泥石流灾害精细化分级监测预警方案。
3.根据权利要求2所述的泥石流灾害精细化分级监测预警方法，其特征在于：所述海拔梯度雨量站是在每200m海拔差上设置雨量站。
4.根据权利要求2或3所述的泥石流灾害精细化分级监测预警方法，其特征在于：在步骤B2中判断流域内泥石流灾害临近后进入步骤C、泥石流发灾预警，首先进行流域调查并安装监测系统，在流域调查中，还包括确定监测断面M处外阻力系数m、监测断面M下游沟道最大安全排泄流量Qp，在安装监测系统中，监测传感器还包括埋设在监测断面M沟床内的总应力传感器，总应力传感器与超声波泥位传感器在同一竖直轴上；然后依如下步骤实施：步骤C、泥石流发灾预警
疑似泥石流发生时，监测传感器实时采集并向控制中心发送泥石流流深监测值h、总应力监测值P；
步骤C11、控制中心依式12计算监测断面M上的泥石流重度γC：式中，h—泥石流流深监测值，单位m，由超声波泥位传感器确定，IC—监测断面M处坡降，根据流域调查确定，P—总应力监测值，kpa，由总应力传感器确定，g—重力加速度常量；
步骤C12、控制中心根据泥石流重度γC判断泥石流发灾信息若γc＜1.3，判断为非泥石流发灾，进入步骤C11；反之判断为泥石流发灾，继续步骤C2；
控制中心实时测算泥石流峰值流量QC，依据QC与下游沟道最大安全排泄流量Qp间流量判别条件确定泥石流发生规模及预警级别；
式中，VC—监测断面M上的平均流速，单位m·s-1，流速传感器实时监测确定，h—泥石流流深监测值，单位m，超声波泥位传感器实时确定，B—监测断面M宽度，单位m，流域调查确定，所述流量判别条件及预警级别是：若0.6Qp≤QC，判断泥石流规模达到黄色预警级别，发出预警信息；进入步骤D、泥石流灾害险情动态预警。
5.根据权利要求4所述的泥石流灾害精细化分级监测预警方法，其特征在于：在步骤C12中，进一步根据泥石流重度γC判断泥石流类型：若1.3≤γc＜1.6，判断为稀性泥石流，若1.6≤γc＜1.9，判断为过渡性泥石流，若γc≥1.9，判断为粘性泥石流。
6.根据权利要求1或2或3或5所述的泥石流灾害精细化分级监测预警方法，其特征在于：在进入步骤D的同时依判别条件判断泥石流规模预警级别并发出相应的预警信号；所述判别条件是：若0.6Qp≤QC＜0.8，预警级别为黄色预警，若0.8Qp≤QC＜Qp，预警级别为橙色预警，若Qp≤QC，预警级别为红色预警。
泥石流灾害险情动态预警方法、精细化分级监测预警方法
技术领域
[0001] 本发明涉及一种灾害预警方法，特别是涉及一种基于演进过程监测的泥石流精细化分级监测预警方法，属于水利工程技术领域、地质灾害防治技术领域。
背景技术
[0002] 地质灾害监测预警是指通过一定的技术手段，监测潜在灾害体的形成条件、激发条件与灾害体的物理特征，在灾害尚未发生或到达危险区之前，预先向受到威胁的地区、人群发出警报。泥石流灾害预警的目的是在泥石流即将发生或者已经发生尚未到达危险区的时刻，发出预警信息，为人员疏散撤离赢得时间，减免人员伤亡。更为理想的泥石流灾害预警技术是能够在灾害过程中，根据监测信息能够动态划分潜在危险区，计算危险度，为不同空间位置的潜在受灾对象提供科学的避灾指导信息。其核心是做到过程监测、精细化分级预警、个性化预警信息发布。
[0003] 现有的监测预警方法主要有通过CCD传感器(照相或录像)对泥石流流域进行监测，直观判断泥石流是否发生以及发生的规模的方法，以及以降雨作为控制参数的监测预警方法，是利用市售雨量计测得降雨数据，通过灾害发生事件与相应的降雨指标统计分析建立的模型进行预警的方法。两类方法都有其缺陷：对于前者，CCD影像信息数据量大，对数据传输要求高，一般需要在野外铺设专用的通讯光纤，建设成本高，且在野外不易实现。而基于GPRS技术的图像抓拍传感器只能定时拍摄泥石流发生区域照片，提供泥石流是否发生的影像，很难给出定量的预警报信息，且往往在雨夜难以保障正常运行和获取信息，容易造成灾害漏报。另外，根据影像估计的灾害规模，不能定量，提供的预警信息比较模糊。对于后者，基于统计规律的模型需要较大的样本才能保证预报的精度，加之泥石流流域地质和地貌条件的区域差异，不同区域的降雨历时—强度曲线具有较大差异，模型很难在建模样本区域以外应用，其推广和应用受到较大限制。
[0004] 现有技术整体而言比较笼统，大多是通过对单一参数的监测，建立统计模型，难以反映灾害的动态过程。这些方法还有一显著缺陷在于不能提供灾害危险性的详细信息，使得疏散撤离实施中难以采取有针对性的措施，防灾减灾效果不理想。并且，已有方法无法做到基于过程监测信息的精细化的预警，只能给出灾害是否发生的判断与发生时间预警信息，不能提供详细的有针对性的临灾信息(如灾害的空间、强度、过程、潜在危险区及其变化等信息)。整体上，预警信息对灾害影响区域内的人群较为抽象，预警结果的可靠性与减灾实务性相对比较弱。
发明内容
[0005] 本发明的目的就是针对现有技术的不足，提供一种泥石流灾害监测预警方法。该方法能够对泥石流灾害的发生演进过程进行精细分级监测，提高预警精度。
[0006] 为实现上述目的，本发明首先提供一种泥石流灾害演进过程中的险情动态预警，用以在判识泥石流灾害形成发生后对险情过程加以精细监测，实现动态警示。其技术方案如下：
[0007] 一种泥石流灾害险情动态预警方法，用以在泥石流灾害发生后监测险情动态过程，实现分级预警；其特征在于：
[0008] 首先进行流域调查，流域调查中，确定泥石流沟道内流通区、危险区范围，在流通区布置监测断面M；现场调查确定监测断面M处坡降IC、监测断面M宽度B、泥石流沟床糙率系数nC、监测断面M至泥石流沟道危险区的距离L、危险区内承灾体分布；通过GIS软件进行流域内各断面的基本地形地貌参数调查，提取流域内DEM数据，根据DEM数据确定计算栅格数量、各栅格位置坐标、各栅格面积，确定数值模拟计算采用的网格的面积a；
[0009] 其次安装监测系统，监测系统包括以无线通信联接的控制中心与监测传感器，监测传感器安装在监测断面M处，包括流速传感器、雨量传感器、监测断面M上方的超声波泥位传感器；
[0010] 其后依如下步骤实施：
[0011] 步骤D1、模拟泥石流沿程运动过程
[0012] 采用二维潜水波方程进行数值模拟，得到泥石流沿程运动方程；其中，运动平衡动量方程依式1、式2，质量守恒方程依式3：
[0013]
[0014]
[0015]
[0016] 式中，g—重力加速度常数，
[0017] Sfx、Sfy—沟道摩擦坡降，
[0018] Sox、Soy—沟道纵坡降，
[0019] t—泥石流持续时间，单位s，根据各地《中小流域暴雨洪水计算手册》方法框算确定，
[0020] u—泥石流在x方向的平均流动速度，单位m/s，流速传感器实时监测确定，[0021] h—泥石流流深监测值，单位m，超声波水位传感器实时监测确定，[0022] v—泥石流在y方向的平均流动速度，单位m/s，流速传感器实时监测确定，[0023] i—一次降雨强度，单位mm/h，雨量传感器实时监测确定，
[0024] 根据泥石流沿程运动方程确定泥石流最终淤埋区域A、泥石流最终不会淤埋的区域A′、每一栅格内的颗粒数Ni,j、泥石流颗粒体积ΔV、每一栅格中泥石流在x、y方向的速度u、v、每一栅格的动态泥深h、泥石流密度ρ；
[0025] 步骤D2、确定每个栅格泥石流淤埋深度Dh
[0026] 依式4计算每个栅格(i,j)的泥石流淤埋深度Dh
[0027]
[0028] 式中，Dh—每一栅格的泥石流淤埋深度，单位m，
[0029] Ni,j—每一栅格内的颗粒数，步骤D1数值模拟计算确定，
[0030] ΔV—泥石流颗粒体积，单位m3，步骤D1数值模拟计算确定，
[0031] A—泥石流最终淤埋区域，单位m2，步骤D1数值模拟计算确定，
[0032] 步骤D3、确定最终淤埋区域A内每个栅格的泥石流最大动能De
[0033] 依式5计算最终淤埋区域A内每个栅格的泥石流最大动能De
[0034] De＝a·maxt＞0[(u2+v2)hρ] 式5
[0035] 式中，De—每个栅格的泥石流最大动能，单位J，
[0036] a—步骤D1数值模拟计算采用的网格的面积，单位m2，流域调查确定，[0037] u、v—每一栅格中泥石流在x、y方向的速度，单位m/s，步骤D1数值模拟计算确定，[0038] h—每一栅格的动态泥深，单位m，步骤D1数值模拟计算确定，
[0039] ρ—泥石流密度，单位kg/m3，步骤D1数值模拟计算确定；
[0040] 步骤D4、确定最终淤埋区域A内承灾体危险等级系数D
[0041] 依式6计算最终淤埋区域A内承灾体危险等级系数D
[0042] D＝De+Dh 式6
[0043] 式中，D—最终淤埋区域A内承灾体危险等级系数，
[0044] De—每一栅格的泥石流最大动能，单位J，步骤D3确定，
[0045] Dh—每一栅格的泥石流淤埋深度，单位m，步骤D2确定，
[0046] 根据D判断各承灾体危险程度：若0＜D≤0.4，判断承灾体处于低危险性，0.4＜D≤0.8，判断承灾体处于中危险性，若0.8＜D，判断判断承灾体处于高危险性；
[0047] 步骤D5、确定灾害警示信息
[0048] 确定灾害警示信息，包括疏散时间t′、承灾体危险等级D、潜在危害区域、逃生路线，具体：
[0049] 依式7计算疏散时间t′，
[0050] t′＝L/Vc 式7
[0051] 式中，t′—疏散时间，单位s，
[0052] L—监测断面M至泥石流沟道危险区的距离，单位m，流域调查确定，[0053] Vc—泥石平均流流速，单位m/s，流速传感器实时监测确定，
[0054] 潜在危害区域是泥石流最终淤埋区域A′，由步骤D1数值模拟计算确定，[0055] 逃生路线是从最终淤埋面积A内各承灾体至泥石流最终不会淤埋的区域A′的最短路径，采用深度优先搜索方法确定；
[0056] 步骤D6、发送灾害警示信息
[0057] 向流域内人群发送灾害警示信息。
[0058] 上述泥石流灾害险情动态预警方法，其技术原理的关键在于数值运算模拟泥石流沿程运动过程。采用二维潜水波方程进行数值模拟，基于非牛顿流体模式及中央有限差分的数值方法来计算洪水、泥石流运动控制方程的集成算法。该算法在洪水灾害管理、基本工程设计、城市淹没分析、泥石流灾害危险性划分等方面应用广泛。为使用该算法适应于本发明解决的技术问题，本发明构思为模型设置了必要假设条件与限制条件。二维模型基本方程式包含连续方程和运动方程，它可以求解泥石流在水平坐标系X方向和Y方向上基本的流速和泥深。在进行数值模拟时，其数据的稳定性是通过在极小的时间步长内进行计算来保证，即在没有调查数据的前提下，通过各种不同泥石流流动路线作为关键因素，用数据稳定的准则限制了时间步长,从而达到避免数据调查情况下又能够有允许时间步长，尽最大努力在合理的时间范围进行泥石流的模拟。
[0059] 以上述泥石流灾害险情动态预警方法为基础，本发明还提供一种泥石流灾害精细化分级监测预警方法，该方法在实施泥石流灾害险情动态预警前先完成泥石流灾害灾前预测，其技术方案如下：
[0060] 利用上述泥石流灾害险情动态预警方法实现的泥石流灾害精细化分级监测预警方法，在实施泥石流灾害险情动态预警前先完成泥石流灾害灾前预测，其特征在于：
[0061] 首先进行流域调查并安装监测系统，在流域调查中，还包括确定形成区径流系数C，在形成区的末端设置监测断面N，调查确定监测断面N以上的流域面积A1、监测断面N宽度B1、监测断面N所在处坡降Ic1，根据现场调查取样与现场土力学实验确定形成区堆积物颗粒平均粒径dm；然后依如下步骤实施：
[0062] 步骤A、泥石流灾害灾前预测
[0063] 步骤A1、确定流域内泥石流形成的临界水量Q
[0064] 依式8、式9、式10计算泥石流灾害临界水量Q，
[0065] Q＝qB1 式8
[0066]
[0067] θ＝arctan(Ic1) 式10
[0068] 式中，Q—泥石流灾害形成临界水量，单位m3/s，
[0069] q—泥石流起动单宽流量，单位m3/s，
[0070] B1—监测断面N宽度，单位m，流域调查确定，
[0071] dm—形成区堆积物颗粒平均粒径，单位m，流域调查确定，
[0072] θ—监测断面N纵坡坡度，单位°，
[0073] Ic1—监测断面N处坡降，单位°，流域调查确定；
[0074] 步骤A2、确定流域内泥石流形成临界降雨阈值iL
[0075] 依式11计算泥石流形成临界降雨阈值iL，
[0076]
[0077] 式中，L—泥石流形成临界降雨阈值，mm/h，
[0078] Q—泥石流灾害形成临界水量，单位m3/s，步骤A1确定，
[0079] C—形成区径流系数，流域调查确定，
[0080] A1—监测断面N以上流域面积，单位ha，流域调查确定；
[0081] 步骤A3、判断泥石流形成条件
[0082] 从气象卫星云图数据中读取流域内未来最大可降雨量ik，若未来最大可降雨量ik≥临界雨量阈值iL，判断流域内出现泥石流发生可能性，发出预警信息；进入泥石流灾害精细化分级监测预警方案。
[0083] 为进一步提高灾害预警的精细程度，上述泥石流灾害精细化分级监测预警方法可进一步优化，是在完成泥石流灾害灾前预测后，即在步骤A3中判断流域内出现泥石流发生可能性后，进行泥石流灾害临灾预警，具体技术方案是：
[0084] 在步骤A3中判断流域内出现泥石流发生可能性后进入步骤B、泥石流临灾预警，首先进行流域调查并安装监测系统，在安装的监测系统中还包括在流域内设置海拔梯度雨量站；然后依如下步骤实施：
[0085] 步骤B、泥石流临灾预警
[0086] 步骤B1、获取梯度降雨观测数据
[0087] 采集流域内不同海拔雨量站雨量观测值，组成流域梯度降雨观测数据，计算降雨观测平均值；
[0088] 步骤B2、判断泥石流临灾条件
[0089] 将降雨观测平均值与临界降雨阈值iL对比，若降雨观测平均值≥临界雨量阈值iL，判定流域内泥石流灾害临近，发出预警信息；进入泥石流灾害精细化分级监测预警方案。
[0090] 为进一步提高灾害预警的精细程度，上述泥石流灾害精细化分级监测预警方法可进一步优化，是在完成泥石流灾害临灾预警后，即在步骤B2中判断流域内泥石流灾害临近后，进行泥石流灾害发灾预警，具体技术方案是：
[0091] 在步骤B2中判断流域内泥石流灾害临近后进入步骤C、泥石流发灾预警，首先进行流域调查并安装监测系统，在流域调查中，还包括确定监测断面M处外阻力系数m、监测断面M下游沟道最大安全排泄流量Qp，在安装监测系统中，监测传感器还包括埋设在监测断面M沟床内的总应力传感器，总应力传感器与超声波水位传感器在同一竖直轴上；然后依如下步骤实施：
[0092] 步骤C、泥石流发灾预警
[0093] 疑似泥石流发生时，监测传感器实时采集并向控制中心发送泥石流流深监测值H、总应力监测值P；
[0094] 步骤C1、泥石流力学参数监测发灾信息
[0095] 步骤C11、控制中心依式12计算监测断面M上的泥石流重度γC：
[0096]
[0097] 式中，H—泥石流流深监测值，单位m，由超声波水位传感器确定，[0098] IC—监测断面M处坡降，根据流域调查确定，
[0099] nC—泥石流沟床糙率系数，根据流域调查确定，
[0100] P—总应力监测值，kpa，由总应力传感器确定，
[0101] g—重力加速度常量；
[0102] 步骤C12、控制中心根据泥石流重度γC判断泥石流发灾信息
[0103] 若γc＜1.3，判断为非泥石流发灾，进入步骤C11；反之判断为泥石流发灾，继续步骤C2；
[0104] 步骤C2、泥石流发生规模预警
[0105] 控制中心实时测算泥石流峰值流量Qc，依据Qc与下游沟道最大安全排泄流量Qp间流量判别条件确定泥石流发生规模及预警级别；
[0106] 泥石流峰值流量Qc依式13计算确定：
[0107] Qc＝VcHB 式13
[0108] 式中，VC—监测断面M上的平均流速，单位m·s-1，流速传感器实时监测确定，[0109] H—泥石流流深监测值，单位m，超声波水位传感器实时确定，
[0110] B—监测断面M宽度，单位m，流域调查确定，
[0111] 所述流量判别条件及预警级别是：若0.6Qp≤QC，判断泥石流规模达到黄色预警级别，发出预警信息；进入步骤D、泥石流灾害险情动态预警。
[0112] 上述泥石流灾害精细化分级监测预警方法可进一步优化，具体优化技术方案分别如下：
[0113] 优化一：在步骤C12中，进一步根据泥石流重度γC判断泥石流类型：若1.3≤γc＜1.6，判断为稀性泥石流，若1.6≤γc＜1.9，判断为过渡性泥石流，若γc≥1.9，判断为粘性泥石流。
[0114] 优化二：在步骤C2中，若判断泥石流规模进入黄色预警级别，发出预警信息，进入步骤D，同时依表1条件进一步判断泥石流规模预警级别并发出相应的预警信号：
[0115] 表1泥石流流量判别条件及预警级别
[0116]
[0117] 与现有技术相比，本发明的有益效果是：(1)本发明提供的一种泥石流灾害险情动态预警方法，该方法能够在对泥石流运动过程数值模拟的基础上，实时预测出泥石流淤埋范围与淤埋深度，再通过叠加流域内承灾体分布数据与物理特征数据分析出承灾体遭受危害的程度、安全疏散线路等信息，并及时向流域内人群发布。方法实现了泥石流灾害监测中对险性的高度动态化与精细化监测，提高监测预警信息的有效性与实际价值。(2)本发明提供的泥石流灾害精细化分级监测预警方法是在流域地形特征数据调查的基础上，在流域内不同位点合理布置监测传感器，通过全方位获取监测影响泥石流灾害的要素信息，将整个泥石流形成、发生、演进过程视为高度动态过程加以预测监测。尤其是从泥石流可能形成的源头起，依据泥石流灾害不同阶段特征规律，分级分类监测不同变量指标。前期集中利用对雨量条件的数值监测与条件判别，后期集中利用对监测断面的数值监测与条件判别，方法整体经由泥石流灾前预测、临灾预警、发灾预警、险情动态预警四级逐步完成泥石流灾害的精细预警，从而逐步提高对泥石流灾害全面预警的效率。与现有的一步、一指标监测预警技术相比，本发明方法充分结合泥石流灾害过程的动态性特征，采用的过程监测有效利用了灾害过程不同阶段不同指标重要性不同自然规律，更为科学；且能够根据监测信息能够动态划分潜在危险区，计算危险度，为不同空间位置的潜在受灾对象提供科学的避灾指导信息，减灾实务性更高。
附图说明
[0118] 图1是都江堰区某流域示意图。
[0119] 图2是泥石流灾害精细化分级监测预警流程图。
[0120] 图3A、图3B、图3C、图3D是泥石流运动数值模拟计算结果图(局部)。
具体实施方式
[0121] 下面结合附图，对本发明的优选实施例作进一步的描述。
[0122] 实施例一
[0123] 如图1～图3所示，用本发明方法对都江堰区某流域中某一子流域(沟道)实施泥石流灾害精细化分级监测预警。
[0124] 图1是都江堰区某流域示意图(图中图标为各类传感器)。该流域居龙门山断裂构造带，陡峻的峡谷地貌使该流域具有以下地形特征：一、流域相对高差大：全流域内最高峰为北端的龙池岗山顶，海拔3280m，最低点位于南端紫平铺水库边，海拔790m，相对高差2490m。各泥石流支沟相对高差723～1605m，最大的是八一沟1605m，最小的是椿牙树沟
723m。二、沟床纵坡陡：整个流域平均纵比降126‰，泥石流支沟纵比降376‰～573‰。三、流域山坡坡度大：流域山坡坡度在30°～70°间；在流域内龙池湖附近及其上游地区的山坡坡度最为陡峻，与之相反的是在流域的东南部分的晏家坪下游到紫平铺水库左岸山坡坡度较缓，这两个区域与流域内别的区域有明显的不同。四、大部分泥石流流域都＜1.0km2，仅个
2
别泥石流沟流域面积＞3.0km 。由此可见，该流域具有山高、坡陡、沟床比降大、支沟面积小的特征。
[0125] 为便于描述，本实施例中以下文字所称“流域”实际均为“子流域”，即该流域中的某一泥石流沟道。
[0126] 图2是泥石流灾害精细化分级监测预警流程图。
[0127] 首先进行流域调查，流域调查中，确定泥石流沟道内流通区、危险区范围，在流通区布置监测断面M；现场调查确定监测断面M处坡降IC＝14.1％、监测断面M宽度B＝30m、泥石流沟床糙率系数nC＝0.14、监测断面M至泥石流沟道危险区的距离L＝3000m、危险区内承灾体分布及；确定形成区径流系数C＝0.8；在形成区的末端设置监测断面N，调查确定监测断面N以上的流域面积A1＝2.14km2、监测断面N宽度B1＝14.2m、监测断面N所在处坡降Ic1＝29.1％，根据现场调查取样与现场土力学实验确定形成区堆积物颗粒平均粒径dm＝0.14m。
[0128] 通过GIS软件进行流域内各断面的基本地形地貌参数调查，提取流域内DEM数据，根据DEM数据确定计算栅格数量(共计1000个)、各栅格位置坐标、各栅格面积，部分数据如表2所示；确定数值模拟计算采用的网格的面积a＝12m2。
[0129] 表2部分栅格基础数据
[0130]
[0131] 其次安装监测系统，监测系统包括以无线通信联接的控制中心与监测传感器，监测传感器安装在监测断面M处，包括流速传感器、雨量传感器、监测断面M上方的超声波水位传感器、监测断面M沟床内的总应力传感器，总应力传感器与超声波水位传感器在同一竖直轴上。在流域内设置海拔梯度雨量站(本实施例中在每200m海拔差上设置雨量站)。
[0132] 其后依如下步骤实施：
[0133] 步骤A、泥石流灾害灾前预测
[0134] 步骤A1、确定流域内泥石流形成的临界水量Q
[0135] 将数据B1＝14.2m、dm＝0.14m m、Ic1＝29.1％代入依式8、式9、式10计算泥石流灾害形成临界水量Q＝12.6m3/s。
[0136] 步骤A2、确定流域内泥石流形成临界降雨阈值iL
[0137] 将数据Q＝12.6m3/s、C＝0.8、A1＝2.14km2代入式11计算泥石流形成临界降雨阈值iL＝26.4mm/h。
[0138] 步骤A3、判断泥石流形成条件
[0139] 在2017年7月18日从气象卫星云图数据中读取流域内未来24h最大可降雨量ik＝49mm/h，因ik≥iL＝26.4mm/h，判断流域内出现泥石流发生可能性，发出预警信息；进入步骤B。
[0140] 步骤B、泥石流临灾预警
[0141] 步骤B1、获取梯度降雨观测数据
[0142] 采集流域内不同海拔雨量站雨量观测值，组成流域梯度降雨观测数据，计算降雨观测平均值＝39.1mm/h。
[0143] 步骤B2、判断泥石流临灾条件
[0144] 因降雨观测平均值≥iL＝26.4mm/h，判定流域内泥石流灾害临近，发出预警信息；进入步骤C。
[0145] 步骤C、泥石流发灾预警
[0146] 此时疑似泥石流发生，监测传感器实时采集并向控制中心发送泥石流流深监测值H、总应力监测值P；
[0147] 步骤C1、泥石流力学参数监测发灾信息
[0148] 将数据IC＝14.1％、nC＝0.14、P＝117kpa代入式12计算监测断面M上的泥石流重度γC＝1.95g/cm，且属于1.9<γC，故判断为泥石流发灾，且泥石流类型为粘性泥石流，继续步骤C2。
[0149] 步骤C2、泥石流发生规模预警
[0150] 将数据VC＝8.9m/s、H＝6m、B＝30m代入式13计算确定泥石流峰值流量Qc＝1602m3/s。因0.6Qp＝1440≤QC，判断泥石流规模达到黄色预警级别，发出预警信息，进入步骤D。同时进一步判断泥石流规模预警级别：依据表1判断泥石流规模预警级别为黄色。
[0151] 步骤D、泥石流灾害险情动态预警
[0152] 步骤D1、模拟泥石流沿程运动过程
[0153] 依式1、2、3建立泥石流数值模拟的二维潜水波方程，得到泥石流沿程运动方程。其中根据《四川省中小流域暴雨洪水计算手册》方法框算泥石流持续时间t＝3.1h。
[0154] 根据泥石流沿程运动方程确定：泥石流最终淤埋区域A、泥石流最终不会淤埋的区域A′、每一栅格内的颗粒数Ni,j、每一栅格中泥石流在x、y方向的速度u、v(部分结果数据见表3)、泥石流颗粒体积ΔV＝0.9m2、泥石流密度ρ＝1.9kg/m3。
[0155] 表3部分栅格计算结果数据
[0156]
[0157] 步骤D2、确定监测断面M的每个栅格泥石流淤埋深度Dh
[0158] 将数据A＝1.8km2、ΔV＝0.9m2、Ni,j代入式4计算得监测断面M的每个栅格(i,j)的泥石流淤埋深度Dh。部分结果如表4所示。
[0159] 表4部分栅格淤埋深度数据
[0160]
[0161] 步骤D3、确定最终淤埋区域A内每个栅格的泥石流最大动能De
[0162] 将数据a＝12m2、ρ＝1.9kg/m3、每一栅格的u、v、h代入式5计算最终淤埋区域A内每个栅格的泥石流最大动能De。部分计算判断结果如表5所示。
[0163] 表5部分栅格最大动能数据
[0164]
[0165] 步骤D4、确定最终淤埋区域A内承灾体危险等级系数D
[0166] 将数据De、Dh代入式6计算最终淤埋区域A内承灾体危险等级系数D，并判断各承灾体危险程度。
[0167] 以某承灾体为例，流域调查确定该承灾体是四层框架结构、坐标31°26'16.00"北，103°34'14.29"东，步骤D3确定该坐标位置处(257号栅格)De＝0.3、步骤D2确定该坐标位置处Dh＝0.6。依式6计算确定D＝0.9，属于高危险程度。
[0168] 图3A、图3B、图3C、图3D所示为泥石流预警监测计算结果(局部)，其中图3A中部的填色区域为泥石流最终淤埋区域A，图3B是分色显示的不同淤埋深度区域，图3C流速分布图，图3D是泥石流危险分区结果(在主沟道两侧根据承灾体危险程度分布区域显示了不同的泥石流危险分区。深灰色区域内承灾体处于高危险性(包含承灾体已被泥石流淤埋的区域)，浅灰色区域内承灾体处于中危险性，近白色区域内承灾体处于低危险性。)[0169] 步骤D5、确定灾害警示信息
[0170] 将数据L＝3000m、Vc＝8.9m/s代入式7计算疏散时间t′＝337s。
[0171] 潜在危害区域是泥石流最终淤埋区域A(如图3所示)。
[0172] 采用深度优先搜索方法计算最终淤埋区域A内各承灾体至泥石流最终不会淤埋的区域A′的最短路径确定逃生路线。具体以淤埋区域到不会被淤埋的区域的最短距离为半径画圆，用决策树分类的方法遍历决定最佳逃生路线。
[0173] 步骤D6、发送灾害警示信息
[0174] 通过公共通讯手段、媒体、网络向流域内人群发送灾害警示信息，包括泥石流规模预警级别、疏散时间t′、承灾体危险等级D、潜在危害区域、逃生路线。
