所谓雨涝模拟分析，

就是模拟降雨量达到一定强度，

城市的哪些区域容易被淹没形成内涝。

雨涝模拟更重要的是提前预测，

可在预测结果的基础上进行实地勘察，

为项目规划、风险防控等工作提供指导作用。

雨涝模拟的原理和思想多种多样，

方法也比较复杂，

本文以模拟广州市雨涝风险区域为例，

介绍一种“体积比较法”的原理和分析思路。

01基本原理

首先搞明白什么是地表径流：地表径流指的是降水或融雪强度一旦超过下渗强度，超过的水量暂时留于地表，当地表贮留量达到一定限度时，即向低处流动，成为地表水而汇入溪流的水量。

本文的原理是基于体积比较法，通俗简明的解释下：

以某种降水量（年均降水量、极端降水量等）为基础，以集水流域为分析单元（简称“集水单元”），计算每个集水单元的径流量，从而得到集水单元的体积，当集水单元的体积（高度）>该集水单元内栅格表面的体积（高度），即形成雨涝淹没区。

该图来源于：ArcGIS帮助文档

技术路线示意图

02数据与模型

原理虽然清楚，但是涉及的内容可不简单，需要制备的数据包括——

降雨量

DEM

集水单元

土地利用

地表径流

1. 降雨量：降雨量的数据可以到中国气象数据网去获取，本次我们模拟了粤港澳30年长时间序列的平均降雨量，从中提取出广州市降水量。（你也可以获取极端降雨数据等等，不同的降水量模拟的结果肯定也是不同的哦）。

第一步，得到带有降水量的气象站点数据——

粤港澳气象站点分布图

第二步，对数据进行分析检验——

数据分析检验示意图

趋势分析与Voronoi图

第三步，模拟预测降水量分布——

克里金插值半变异函数拟合分析

粤港澳年均降水量的插值模拟

2. 集水单元：区域的模拟分析，都要有一个分析单元，一般而言，和水有关的分析，常常结合DEM数据，利用水文分析，提取出集水流域，作为分析单元，本例也不例外——

广州市集水单元

3. 地表径流：地表径流量需要根据降雨量和径流系数（不同地表覆盖的径流系数）来计算，参考的计算模型为“双评价指南”的模型：

计算地表径流基本的步骤是：

第一步：为每种地类赋予对应的径流系数；

第二步：集水单元与地类叠加，计算出每个集水单元的加权平均径流系数；

第三步：将计算出加权平均径流系数的集水单元与降水量的数据叠加计算，得到每个集水单元的径流量。

集水单元流量模拟示意图

03雨涝风险斑块识别

下面我们开始进行雨涝风险斑块的识别。

第一，计算每个集水单元的集水体积：这一步倒是不难，就是用每个集水单元径流量与每个集水单元面积相乘即可；

第二，体积比较法识别淹没区域：如果只是分析一个集水单元的雨涝淹没情况，可以利用“表面体积”工具，不断的试验集水高度，直到接近该集水单元的体积，则该高度以下即为淹没区，当然，这个方法工作量很大，而且只适用于一个集水单元，很多集水单元的话，这个方法就扑街了；

下面，就需要介绍下体积比较法，简明介绍下基本思想——

每个集水单元的体积是一定的，某区域集水单元的面积<该集水单元区域地表的面积，则，该区域集水高度>该区域的地表高度，即为淹没区！

这基本思想看起来简单，但是涉及到的操作可真是不少！通过集水单元栅格面积的计算与比较，基于上述思想，得到了集水高度>地表高度的区域，也就是潜在的雨涝淹没风险区——

第三，风险集中度分级：这一步又简单了，就是将淹没区进行核密度的分析，根据集中程度进行分级就可以啦——

最后再说明一下，城市雨涝淹没分析原理多样，方法也比较复杂，本文主要基于“体积比较法”的一个思路框架。

​