-
城市内涝是指强降水或连续性降水超过城市排水能力致使城市内产生积水灾害的现象.内涝灾害发生时，城市交通、网络、通信、水、电、气、暖等生命线工程系统瘫痪，社会经济活动中断，其灾害损失已远远大于因建筑物和物资破坏所引起的直接经济损失.随着社会经济的迅猛发展，城市规模不断扩大、人口快速增长，大量的硬质铺装，硬化路面越来越多，使雨水无法下渗透或渗透慢，容易造成积水，城市内涝灾害日趋严重. 2007年7月16-19日，重庆出现了百年以来范围最广、持续时间最长、强度最大的强降雨天气，道路中断积水严重，造成的巨大损失已引起各级政府的高度重视，故迫切需要针对城市内涝的有效预警和区划方法来为政府提供服务.
城市内涝预警预报的关键是确定城市内涝临界雨量，目前有关临界雨量的计算方法主要有实测雨量分析法、降雨灾害频率分析法、产汇流对比分析法及暴雨临界曲线法、分布式水文模型法等[1-12]，也有很多专家学者开发了针对城市内涝的仿真预警模型：如天津市暴雨沥涝仿真模型[13]、城市排水系统非恒定流模拟模型[14]等.应用这些方法可以计算得出内涝发生时的临界雨量，但是由于此类物理模型法涉及城市排水管网数据，而城市管网的数据获取难度大，且随着年代久远很可能由于杂物堵塞等因素造成其排水功能与设计不符.
重庆属山地城市，地形分布复杂，针对城市内涝预警和区划方面的研究较少.本文采用了德国研发的二维水动力模型“FloodArea”(Geomer，2003)建立城市内涝情景模拟，以不同重现期下的暴雨雨型为输入，研究沙坪坝区致灾临界雨量的计算方法和基于GIS栅格叠加运算原理的风险区划方法.
HTML
-
沙坪坝区属于川东平行岭谷低山丘陵区的一部分，呈丘陵、台地和低山组合的地貌结构，西部为缙云山低山区，中部歌乐山海拔高度在550~650 m之间，西面为中梁山.嘉陵江由北往东南流经沙坪坝区19.3 km，国道、省道等高等级公路5条，成渝、渝长、渝遂、上界、绕城高速公路纵横域内.
-
地理信息数据：DEM采用的SRTM空间分辨率为30 m数据，投影方式为ALBERS等面积圆锥投影，WCG1984大地坐标系，中央经线为105°E，双标准纬线分别为25°N和47°N.与DEM相同投影的1：25万沙坪坝区土地利用数据.在选择模拟的范围上，对沙坪坝区行政边界进行了一定范围的扩大，以确保降水在整个沙坪坝区模拟范围内能够下渗和汇流，制作风险区划图时根据沙坪坝区的行政边界对模拟的结果进行了裁剪.气象观测数据：沙坪坝区、渝北区、巴南区、璧山县和北碚区逐分钟降水自记纸观测资料，资料年限为1961-2013年，由重庆市气象信息与技术保障中心提供.
-
根据电子化后的逐分钟降雨量自记纸资料，采用皮尔逊—Ⅲ型分布曲线拟合计算得到沙坪坝5，10，20，30，50和100 a重现期的24 h历时设计暴雨量，见表 1.采用耿贝尔分布曲线拟合得到24 h历时各重现期设计暴雨量分布情况，见表 2.
-
FloodArea模型计算基于水动力方法，同时考虑了一个栅格的周围8个单元.对邻近栅格单元的泻入量由manning-stricker公式计算[15].坡度由单元之间的最低的水位和最高的地形高程之间的差异所决定，对每一个单元都进行计算，相邻单元的流量长度被认为是相等的，位于对角线的单元，以不同的长度算法来计算.
水力糙度指的是流体力学上的粗糙度，是衡量河道或冲沟边壁形状不规则性和粗糙程度影响的一个综合性系数，目前关于水力糙度的研究较多[16].利用地面覆盖类型，根据张洪江等人[7]的研究结论对不同地面覆盖类型的水力糙度赋初值. 图 1使用了区域内土地利用类型及其糙率转变的数值进行赋值，表 3为具体土地利用类型转变赋值.
SCS(Soil Conservation Service)模型应用于无资料地区径流模拟，依据地表覆盖类型来确定SCS模型的参数.模型能够客观反映土壤类型、土地利用方式及前期土壤含水量对降雨径流的影响，其显著特点是模型结构简单、所需输入参数少，是一种较好的小型集水区径流计算方法[17]. SCS模型计算公式如下：
其中：Q为产流量；P为过程降水量；K为产流系数；S为潜在入渗量，CN反映流域前期土壤湿润程度、坡度、土壤类型和土地利用现状的综合特性，可以较好地反映下垫面条件对产汇流过程的影响，CN值由地表覆盖类型确定，可查表获得[18].根据地面覆盖类型，通过公式(1)-(3)计算得出流域地表产流系数值.在计算整个过程降雨量时采用的是沙坪坝及其周边4个气象台站5 a重现期逐小时降雨量资料，利用泰森多边形方法计算了整个降雨过程产流系数空间分布图(图 2).
-
风险区划方法采用矩阵分析法如图 3所示，根据暴雨诱发的城市内涝影响因子等级条件(表 4)进行重采样划分，然后通过矩阵法将影响因子进行叠加分析得到暴雨诱发的城市内涝风险等级区划图.本文采用的城市内涝风险影响因子为积涝水深、DEM和坡度.积涝水深是基于不同重现期的设计暴雨雨型的二维水动力模型的动态模拟结果.重庆市地形起伏大与平原地区不同，地势的高低和坡度将直接影响降雨量的汇流情况.将DEM和地表坡度作为城市内涝风险区划矩阵中的影响因子进行4个等级划分(表 4)，DEM越高处汇水风险较低，反之越低地段越容易积水，风险较高.而坡度越小越容易汇流，反之坡度越大越不容易积水.采用重采样方法将DEM和坡度划分为4个级别(如图 4-图 5)，从图中可以看出嘉陵江、中梁山、缙云山以及山谷地段分布清晰且与坡度分布相互吻合.
1.1. 研究区概况
1.2. 数据源
1.3. 方法
1.3.1. 雨量数据处理
1.3.2. FloodArea模型
1.3.3. 风险区划
-
将24 h历时设计暴雨雨型5 a重现期逐时雨量作为模型输入雨量，采用FloodArea水文模型中的暴雨模块进行动态模拟，并对输出的淹没结果(图 6-图 9)进行分析，从图中可以看出整个汇流和淹没过程非常清晰.设计雨型的峰值区在第3 h和4 h(表 2)，从淹没结果中可以看出在第4~9 h区间，积水也达到此次过程的一个峰值水平，之后积水开始逐渐减缓.再统计沙坪坝区5 a重现期暴雨过程模拟分时淹没面积表(表 5)，在第4 h，超过0.6 m淹没水深的面积陡增.由于降水大部分集中在过程前半段且降水的汇流需要时间，第9 h之后淹没深度在0.6~1.2 m的面积开始逐渐减少，1.2~1.8 m淹没面积在第19 h也开始逐渐减少，而淹没水深大于1.8 m的范围在继续增加.由于城市排水管道数据的可获取性有限，管网参数并未考虑，仅通过调整地表渗透率来弥补排水管网的缺失，可能导致模拟结果与真实情况存在误差.
配合隐患点的经纬度坐标提取整个模拟过程各个隐患点淹没水深数据，建立雨量-淹没水深关系的回归方程.本文分析了此次模拟过程中淹没较深的5个隐患点，结果分析2 h累计降雨量与逐时的淹没深度的拟合相关系数最高，均在0.7以上，根据漫堤淹没预警点0.6 m、1.2 m和1.8 m的致灾临界条件，计算得出不同等级累积雨量(表 6).将本次5个隐患点模拟过程的逐时淹没深度与降水进行了对比(图 9)，可以看出此过程的淹没深度与降雨过程联系紧密，当降雨量在第4 h达到峰值，随着雨水汇流各隐患点淹没深度也达到峰值，随着降雨量逐渐减小，各隐患点淹没深度也开始逐步下降.
-
将5，10，20，30，50，100 a重现期24 h雨型的模拟结果按照表 4中等级进行重采样划分，得到不同重现期的沙坪坝城市内涝淹没因子空间分布图(图 11-图 14)和不同等级淹没面积占整个沙坪坝区面积的比重表(表 7).从空间上可以看出不同重现期的淹没区域空间分布相似，易产生积涝地区分布在缙云山和中梁山之间低洼地，主要为梁滩桥、五福、白鹤、西永及磁器口街道.从重现期尺度来分析，重现期越大，1级和2级淹没深度范围越小，而3级和4级淹没深度范围越大.
-
根据矩阵法将淹没等级、DEM影响因子和坡度影响因子进行叠加计算，得到沙坪坝区城市内涝风险区划图，因为不同重现期的内涝风险区划图在空间分布上大致相似，故仅展示了5 a重现期的沙坪坝城市内涝风险区划图(图 15).从图中可以看出高风险区主要集中在贯穿沙坪坝区的梁滩河附近和磁器口街道附近，这些区域海拔较低且坡度不大，容易产生积涝.
2.1. 沙坪坝区城市内涝致灾临界阈值计算
2.2. 沙坪坝区城市内涝风险分析
2.3. 沙坪坝区城市内涝风险区划分析
-
本文以重庆市沙坪坝区为例，利用FloodArea二维水动力模型对5，10，20，30，50，100 a重现期24 h暴雨雨型进行了动态模拟，得到不同隐患点致灾临界雨量和风险区划图.
结果表明：1) FloodArea在区域性山洪致灾动态模拟方面有大量的应用，在城市内涝的动态模拟上运用较少，由于缺乏城市管网数据，过程的排水量可能存在误差，但是通过增加模型中地表透水交流率参数设定进行相应弥补，模拟结果能够反映出降雨过程的地表汇流情况.
2) 通过提取整个模拟过程各个隐患点淹没水深数据，建立雨量-淹没水深关系的回归方程计算的致灾临界阈值具有一定的物理基础.
3) 从矩阵法计算的沙坪坝城市内涝风险区划图中可以看出，高风险区主要集中在贯穿沙坪坝区的梁滩河附近和磁器口街道附近，这些区域海拔较低且坡度不大，容易产生积涝，与实际情况相符.
由于城市排水管网数据较难获取，FloodArea水文模型计算的淹没深度数值相对较大，制作的致灾临界阈值和风险区划图仅限于技术研究.基于水文模型的风险临界雨量和风险区划的研究属于近年来开展的新兴研究，对于水文模型的选择和优化还需要进一步深入研究.
DownLoad: