ArcGIS雨涝风险模拟

所谓雨涝模拟分析,

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

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

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

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

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

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

方法也比较复杂,

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

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

01 基本原理

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

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

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

该图来源于:ArcGIS帮助文档

技术路线示意图

02 数据与模型

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

降雨量

DEM

集水单元

土地利用

地表径流

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

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

粤港澳气象站点分布图

第二步,对数据进行分析检验------

数据分析检验示意图

趋势分析与Voronoi图

第三步,模拟预测降水量分布------

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

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

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

广州市集水单元

**3. 地表径流:**地表径流量需要根据降雨量和径流系数(不同地表覆盖的径流系数)来计算,参考的计算模型为"双评价指南"的模型:

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

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

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

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

集水单元流量模拟示意图

03 雨涝风险斑块识别

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

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

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

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

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

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

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

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

相关推荐
林泽毅8 分钟前
SwanLab x EasyR1:多模态LLM强化学习后训练组合拳,让模型进化更高效
算法·llm·强化学习
小林熬夜学编程9 分钟前
【高并发内存池】第八弹---脱离new的定长内存池与多线程malloc测试
c语言·开发语言·数据结构·c++·算法·哈希算法
刚入门的大一新生16 分钟前
归并排序延伸-非递归版本
算法·排序算法
独好紫罗兰21 分钟前
洛谷题单3-P1980 [NOIP 2013 普及组] 计数问题-python-流程图重构
开发语言·python·算法
独好紫罗兰26 分钟前
洛谷题单3-P1009 [NOIP 1998 普及组] 阶乘之和-python-流程图重构
开发语言·python·算法
曦月逸霜37 分钟前
蓝桥杯高频考点——高精度(含C++源码)
c++·算法·蓝桥杯
ゞ 正在缓冲99%…1 小时前
leetcode152.乘积最大子数组
数据结构·算法·leetcode
闯闯爱编程1 小时前
数组与特殊压缩矩阵
数据结构·算法·矩阵
秋风战士2 小时前
通信算法之255:无人机频谱探测设备技术详解
算法·无人机
laimaxgg2 小时前
数据结构B树的实现
开发语言·数据结构·c++·b树·算法