基于降雨频率的泥石流淤积危险范围划分方法与流程

本申请属于泥石流灾害防治
技术领域：
，具体地说，涉及一种基于降雨频率的泥石流淤积危险范围划分方法。
背景技术：
：泥石流是山区常见的地质灾害，其暴发突然、来势凶猛、历时短暂，在造成严重灾害的同时，由于其高含沙量，往往在流域下游淤积形成大面积的泥石流滩地。泥石流滩地往往地势平缓、水源充足、热量丰富，是较好的土地资源，具有广阔的开发前景，是山区发展经济和改善生态环境的重要途径。泥石流滩地多分布于泥石流沟口及平坦的沟谷的两侧，由于平坦开阔的地形条件以及相对较好的水热条件，山区居民常常在稳定的滩地上进行各种农业生产、建设村落、修筑道路等活动。但是，泥石流滩地利用面临潜在泥石流灾害的威胁，因此划分泥石流淤积危险范围是防御泥石流风险的重点。目前，有关泥石流淤积危险范围的划分方法主要有两类：①利用海量调查、实验数据统计分析拟合方程，如泥石流冲出距离与一次泥石流冲出量的关系式；②利用数学、物理模型开展数值模拟，通过模拟得到的流速、泥深确定泥石流淤积淤埋危险范围。然而，这两种方法主要存在如下问题：一是基于室内实验数据的数理统计拟合方程得到泥石流最大冲出距离和最大冲出宽度确定的泥石流范围，很难反映泥石流淤积实际情况，同时不能获得淤埋范围的深度、流速，且该类方法很难应用到实际，尤其区域不同，泥石流特征差异较大，导致数理统计的淤埋危险范围差异很大；二是利用数学、物理模型开展的数值模拟参数选取难度较大，且方程求解和模拟过程复杂，而且没有考虑淤积淤埋范围内是否涉及地表固体物质冲刷、侵蚀等加大危险范围的因素。因此，现有技术的泥石流淤积危险范围的划分方法比较复杂，且可靠性不高。技术实现要素：有鉴于此，本申请所要解决的技术问题是提供了一种基于降雨频率的泥石流淤积危险范围划分方法，以提升划分出的泥石流淤积危险范围的准确性，并降低泥石流淤积危险范围划分方法的复杂性。为了解决上述技术问题，所述的方法包括：根据被检测流域的地形图确定泥石流主沟道位置、泥石流初始过流断面m0n0及其初始过流断面的宽度，并确定初始过流断面的水力坡降；根据历史数据资料，设定降雨频率，计算对应降雨暴发频率下泥石流的峰值流量；根据泥石流的峰值流量、初始过流断面的宽度和水力坡降，计算初始过流断面平均流深h0；在泥石流主沟道上确定n个平均流深为h0的过流断面mini，i＝1,2,3,…n，n≥2；依次连接mimi+1，及nini+1，i＝0,1,2,3,…n，n≥2，获得的范围即为对应暴发频率下的泥石流淤积危险范围。如上所述的方法，其中，所述计算初始过流断面平均流深h0，具体为可以为：将过流断面概化呈倒立的等腰三角形，通过以下公式计算求得h0，q＝l0h0v0，其中，q为泥石流峰值流量，单位为m3；l0为沟口过流断面m0n0的宽度，单位为m；h0为沟口对应降雨频率下峰值流量的平均流深，单位为m；v0为泥石流在沟口过流断面m0n0上的平均流速，单位为m/s；m为粘性泥石流流速系数，根据《泥石流灾害防治工程勘察规范da-t0220-2006》内插确定；r为沟口过流断面m0n0的水力半径，单位为m；i为沟口过流断面m0n0处泥石流水力坡度，单位为‰。如上所述的方法，具体的，所述在泥石流主沟道上确定n个平均流深为h0的过流断面mini，具体包括：确定自初始过流断面开始沿泥石流流动方向的等高线ai，泥石流主沟道与等高线ai相交于ci，在ci和ci+1之间确定平均泥深为h0的点oi，在oi点做垂直主沟道的过流断面，与等高线ai相交于mi、ni，从而得过流断面mini，其中，i＝1,2,3,…n，n≥2。如上所述的方法，其中，所述危险范围在同一横断面上，越靠近主沟，流深越大，越靠近主沟两侧，流深越小。如上所述的方法，进一步的，根据过流断面的泥深和单宽流量，确定所述危险范围的危险等级，所述危险等级包括：高、中、低。与现有技术相比，本申请可以获得包括以下技术效果：1、本申请的泥石流淤积危险范围划分方法计算简单、简捷易行、易于推广实施，针对山区缺资料地区小流域具有较高的应用价值；2、基于降雨频率，可实现各种爆发频率条件下，快速获取泥石流堆积区灾损土地淤埋风险范围，可广泛应用于泥石流监测预警、乡村规划等方面；3、本申请的泥石流淤积危险范围划分方法，针对小流域，考虑了淤积淤埋范围内是否涉及地表固体物质冲刷、侵蚀等加大危险范围的因素，提升了划分出的泥石流淤积危险范围的准确性。附图说明此处所说明的附图用来提供对本申请的进一步理解，构成本申请的一部分，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：图1是本申请基于降雨频率的泥石流淤积危险范围划分方法流程图；图2是本申请实施例的泥石流主沟道及沟口过流断面m0n0的示意图；图3是本申请实施例的泥石流主沟道及过流断面mini的示意图；图4是本申请实施例的泥石流主沟道淤积范围的示意图。具体实施方式以下将配合附图及实施例来详细说明本申请的实施方式，藉此对本申请如何应用技术手段来解决技术问题并达成技术功效的实现过程能充分理解并据以实施。本申请提供的泥石流淤积危险范围划分方法特别适用于无大比例尺地形图及观测数据的无资料地区的泥石流堆积区，该方法主要以流域所在区域泥石流爆发频率以及地形参数为依据，快速完成由泥石流堆积导致的灾损土地淤埋范围及区划。图1是本申请基于降雨频率的泥石流淤积危险范围划分方法流程图。如图1示，本实施例的基于降雨频率的泥石流淤积危险范围划分方法包括以下内容。s101：根据被检测流域的地形图确定泥石流主沟道位置、泥石流初始过流断面m0n0及其初始过流断面的宽度，并确定初始过流断面的水力坡降。图2是本申请实施例的泥石流主沟道及沟口过流断面m0n0的示意图。具体的，可以通过现场调查、测绘、勘查的图形资料以及dem(digitalelevationmodel,数字高程模型)确定泥石流主沟道的地理位置及沟口初始断面位置。s102：根据历史数据资料，设定降雨频率，计算对应降雨频率下泥石流的峰值流量。s103：根据泥石流的峰值流量、初始过流断面的宽度和水力坡降，计算初始过流断面平均流深h0。初始过流断面平均流深h0是指对应频率和泥石流峰值流量q条件下流过初始过流断面m0n0的平均流深。具体的，可以将过流断面概化呈倒立的等腰三角形，通过以下公式计算求得h0，q＝l0h0v0，其中，q为泥石流峰值流量，单位为m3；l0为沟口过流断面m0n0的宽度，单位为m；h0为沟口对应降雨频率下峰值流量的平均流深，单位为m；v0为泥石流在沟口过流断面m0n0上的平均流速，单位为m/s；m为粘性泥石流流速系数，根据《泥石流灾害防治工程勘察规范da-t0220-2006》内插确定；r为沟口过流断面m0n0的水力半径，单位为m；i为沟口过流断面m0n0处泥石流水力坡度，单位为‰。s104：在泥石流主沟道上确定n个平均流深为h0的过流断面mini，i＝1,2,3,…n，n≥2。s105：依次连接mimi+1，及nini+1，i＝0,1,2,3,…n，n≥2，获得的范围即为对应暴发频率下的泥石流淤积危险范围。图3是本申请实施例的泥石流主沟道及过流断面mini的示意图。参照图3所示，在如上所述的方法中，所述在泥石流主沟道上确定n个平均流深为h0的过流断面mini，具体可以包括：确定自初始过流断面开始沿泥石流流动方向的等高线ai，泥石流主沟道与等高线ai相交于ci，在ci和ci+1之间确定平均泥深为h0的点oi，在oi点做垂直主沟道的过流断面，与等高线ai相交于mi、ni，从而得过流断面mini，其中，i＝1,2,3,…n，n≥2。依次连接mimi+1，及nini+1，获得的范围即为对应暴发频率下的泥石流淤积危险范围。图4是本申请实施例的泥石流主沟道淤积范围的示意图。根据已有地形图或dem生成小流域等高线，或者现场调查、测绘、勘查确定泥石流流域主沟道基本地理地形数据，包括主沟道地理位置、泥石流过流断面m0n0及其宽度l0，并确定初始过流断面m0n0处泥石流水力坡度i。如上所述的方法，其中，所述危险范围在同一横断面上，越靠近主沟，流深越大，越靠近主沟两侧，流深越小。具体的，可以根据过流断面的泥深和单宽流量，确定危险范围的危险等级，危险等级包括：高、中、低。对泥石流堆积区的危险等级的划分往往以泥石流强度为划分依据，常用的特征参数有泥石流流速、泥深、冲击力等。例如，可以在泥石流堆积区淤积范围内，进一步划分过流断面，间隔以等高线间距的1/2～1/5，并确定指定过流断面的水力坡度，从而以泥深和单宽流量为基本参数，将淤积范围进行危险性高、中、低划分，参照表1。表1泥石流淤积危险范围的危险等级表泥石流危险等级说明高h≥0.5m或q≥1m2/s中0.1m＜h＜0.5m或0.2m2/s＜q＜1m2/s低h≤0.1m或q≤0.2m2/s具体应用案例：南坝镇为龙门山中央断裂北段途径位置，在该处可见到多处龙门山断裂带的地表破裂区域。窖子沟介于南坝镇和建康村之间，距离南坝镇约5.5km，窖子沟汇入涪江支流，由于断裂带横穿窖子沟流域中下游，造成流域中下部山体破碎，崩塌滑坡现象明显，通过测量，穿越窖子沟的龙门山断裂带距离主河道约750m。目前，窖子沟中下游左岸一滑坡发育程度较高，该滑坡体长约550m，宽350m，由于松散滑坡体挤占主沟道，堵塞较严重，在未来强降雨条件下，暴发泥石流的概率较大，可能造成主河堵塞回水，危及河道上方建康村村落及河道两侧农耕地。利用开源dem数据，获取流域等高线，划定泥石流淤积范围。本应用案例以窖子沟流域泥石流堆积区灾损土地为例，对该流域淤埋范围及区划进行案例分析。(1)确定被检测小流域地理参数通过多次野外实地调查、测量获得流域地形、地物，结合地理空间数据云下载高精度dem数据，利用arcgis等软件进行矢量，得到窖子沟流域地形图(主要是堆积区)。从而确定窖子沟主沟道位置、泥石流初始过流断面m0n0及其宽度2l0，可得2l0＝52.78m，确定此过流断面水力坡降i＝252‰。根据走访调查以及窖子沟泥痕调查，得到窖子沟在2010年四川省群发性泥石流暴发中泥石流峰值流量155m3/s，大致相当于泥石流暴发频率为50年一遇的泥石流流量。(2)确定被检测小流域不同降雨频率下泥石流峰值流量q南坝镇窖子沟流域位于四川省中部，根据四川省中小流域水文手册，泥石流暴发频率分别以20年一遇(5％)、50年一遇(2％)、100年一遇(1％)进行设置，泥石流过程为30min，得到窖子沟泥石流流量情况计算表，见表2。表2平武县南坝窖子沟流域泥石流总量计算表(3)确定初始过流断面平均流深h0得到窖子沟不同频率条件下泥石流峰值流量，从而根据初始过流断面l0＝52.78m和水力坡降i＝252‰，进一步获得初始过流断面不同频率条件下泥石流峰值流量对应的平均泥深，见表3。表3平武县南坝窖子沟流域初始过流断面平均泥深(4)确定过流断面mini的位置由前所知，依据设定频率条件下的平均泥深为标准，假定自初始过流断面以下所有过流断面位置平均泥深均为此值，即20年一遇平均泥深为1.8m，50年一遇平均泥深为2.1m，100年一遇平均泥深为2.7m。从而确定自初始过流断面开始，沿泥石流流动方向的等高线ai，i＝1,2,3，…n。其中，以a1和a2等高线为例，其高差h1-2为10m，泥石流主沟道与等高线分别相交于c1、c2，测量c1和c2间图上距离为c1-2，从c1开始沿泥石流流动方向c1-2段定点o1(该点平均泥深为h0)，其c1至o1间图上距离为c1o1。确定o1点做垂直主沟道过流断面，与等高线a1相交于m1、n1，从而得到第一个过流断面m1n1。进而利用a2和a3等高线，a3和a4等高线，…，确定泥石流主沟道与等高线交汇点c1、c2，…，以及第2、3条过流断面m2n2、m3n3，…，可以参考图3和图4所示。(5)确定泥石流堆积区淤积范围依次连接mimi+1，及nini+1，i＝1,2,3,…n，所得范围即为南坝窖子沟设计频率条件下泥石流堆积区淤积范围，可以参考图4所示。通过对比，在泥石流暴发频率为50年一遇情况下，圈定的泥石流堆积区淤积范围与实际窖子沟堆积扇比较一致；当泥石流暴发频率为20年一遇，圈定的泥石流堆积区淤积范围仅分布在堆积扇主沟两侧；当泥石流暴发频率为100年一遇，圈定的泥石流堆积区淤积范围基本全部覆盖堆积扇。(6)淤积危险性划分确定泥石流堆积区淤积范围，进一步以淤积范围为依据，可以确定范围内淤积危险等级，在同一横断面上，越靠近主沟，流深越大，越靠近主沟两侧，流深越小。仍然以南坝窖子沟流域50年一遇泥石流暴发频率为例，通过以流深为依据进行说明，窖子沟流域堆积区初始过流断面至主河之间高程范围介于900-750m，以10m高差间隔可确定指定过流断面16个。通过测量各横断面宽度以及计算水力坡度，可求得16个过流断面平均泥深hi，i＝1,2,3…16。计算结果见表4。表4南坝窖子沟50年一遇泥石流暴发频率过流断面特征参数计算表以每一个横断面最大泥深为基本参数，根据前面表3确定每一个横断面的泥石流淤积危害程度。根据表4窖子沟流域泥石流暴发频率在50年一遇的情况下，堆积区主沟道最大泥深介于0.36～0.73m，基本上属于中度和高度危险区。综上所述，本申请提供的基于降雨频率的泥石流淤积危险范围划分方法，根据被检测流域的地形图确定泥石流主沟道位置、泥石流初始过流断面m0n0及其初始过流断面的宽度，并确定初始过流断面的水力坡降；根据泥石流的暴发频率，通过查阅历史数据资料获得对应暴发频率下泥石流的峰值流量；根据泥石流的峰值流量、初始过流断面的宽度和水力坡降，计算初始过流断面平均流深h0；在泥石流主沟道上确定n个平均流深为h0的过流断面mini，i＝1,2,3,…n，n≥2；依次连接mimi+1，及nini+1，获得的范围即为对应暴发频率下的泥石流淤积危险范围。本申请的泥石流淤积危险范围划分方法具有简捷易行、易于推广实施的特点，针对山区缺资料地区小流域具有较高的应用价值，且划分出的泥石流淤积危险范围准确性高，可广泛应用于泥石流监测预警、乡村规划等方面。上述说明示出并描述了本申请的若干优选实施例，但如前所述，应当理解本申请并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本申请构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本申请的精神和范围，则都应在本申请所附权利要求的保护范围内。当前第1页12