基于ArcGIS的洪水淹没分析与三维模拟_孙君

《基于ArcGIS的洪水淹没分析与三维模拟_孙君》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、20122012Vo1.7No.3孙君等：基于ArcGIS的洪水淹没分析与三维模拟第7卷第3期基于ArcGIS的洪水淹没分析与三维模拟孙君，奚赛英，尤迪，郑付涛（常州市规划设计院,江苏常州213003）摘要：洪水淹没范围的确定是洪灾损失评估和防洪决策的核心环节。基于TIN数据，运用ArcMap，采用“无源淹没分析”方法对区域天然防洪能力进行划分；实现了在给定水位条件下，对洪水淹没范围提取与统计计算，建立了洪水水位高程和淹没面积关系公式，并用于洪水淹没快速预测；运用ArcScene，对水位抬升的“无源渐进淹没”情况进行了三维模拟。关键词：洪水；

2、GIS；无源淹没；淹没范围；三维模拟中图分类号：TP311.5文献标识码：A文章编号：1007-1903（2012）03-0031-030引言洪水淹没是一个动态而至平衡的过程。确定洪水最终淹没范围，有2种概化模型：(1)基于水位的洪水淹没我国是一个洪水灾害频发的国家，大约2/3的国土面范围计算：给定某一洪水水位H，由此推算出洪水的淹没[1]积有着不同类型和不同危害程度的洪水灾害。因此，快范围。该模型相对简便。(2)基于水量的洪水淹没范围计速、准确地模拟预测洪水淹没范围及面积，对防洪减灾算：即在给定某一洪水水量Q的条件下，计算相应的洪水具有重要

3、意义。淹没范围。在应对非调度洪水灾情分析时，往往需要利[7]防洪减灾工作自20世纪90年代以来，在水动力-水用获取的淹没区洪水水位或水量。因洪水水位数据的文模型基础上，利用ArcGIS强有力的空间分析和可视化获取较为容易，故采用基于水位的淹没分析。功能，模拟显示洪水淹没区，进行非工程措施防洪减灾，基于水位的淹没分析，分两种情形：(1)无源淹没：凡是研究热点。是高程值低于给定水位的点均为淹没区，相当于整个地结合相关学术文献，发现ArcGIS在洪水淹没分析方区大面积均匀降水，所有低洼处都可能积水成灾；(2)有[2-6]面多有应用。但需要编制复杂程

4、序或应用复杂算法，源淹没：考虑“流通”淹没的情况，即洪水只淹没它能流往往过于繁琐。到的地方，相当于高发洪水流域泛滥，例如洪水决堤，或[8]本文提供了无需编程就可以简单迅速地计算出洪水局部暴雨引起的暴涨水向四周扩散。对于洪水源头不淹没面积的方法。以ArcGIS为基本处理方法，以TIN数易确定、地势相对平坦地区，无源淹没模型较为贴切。据为基础数据，运用ArcMap自身功能，对区域天然防洪依据上述分析方法适用范围和本文研究对象的特能力进行划分，求得低于一定高程的洪水淹没范围，实点，采用基于水位的洪水淹没范围计算中的无源淹没模现了洪水淹没面积的计算，

5、进而建立洪水水位与淹没面型进行分析计算与三维模拟。积间的关系公式，并运用公式进行淹没面积快速预测。运用ArcScenen软件，简便建立了地形和洪水模型，并对1.2划分区域天然防洪能力洪水渐进情况下的淹没情况进行了三维模拟。防洪能力涉及因素较多，如地形地貌、地面高程、河道畅通性能、防洪排涝设施建设情况等。这里以某城市1洪水淹没范围分析建成区为研究对象，根据地区防洪除涝水位控制标准，只考虑利用地面高程评价防洪适宜性，对其天然防洪1.1确定洪水淹没分析方法能力作定性划分。研究区域内高程-9.1m～107.6m，利用作者简介：孙君（1983-），男，

6、江苏扬州人，硕士，从事岩土工程与地理信息系统相关研究，email：sunjun_czpad@163.com3120122012第7卷第3期分析探讨Vo1.7No.322ArcMap，结合地区防洪水位要求，以50年防洪设防水位面积为181306884m，即约181.3km，如图2所示。黄海标高3.75m、200年防洪设防水位标高4.05m为界划分系统建设区内各区域天然防洪能力（图1），直观展示发生各种级别的洪灾时，可能淹没的区域，以便有针对性地进行规划、决策。图2淹没面积计算同样的，计算洪水水位高程为1.0m、1.5m、2.5m、图1天然防洪能

7、力分区图3.9m、4.05m、5.0m的淹没面积，详见表1：表1淹没面积统计1.3洪水淹没范围分析与预测洪水水位总面积淹没面积淹没比例(km2)(km2)(%)（1）数据预处理高程(m)1.00446.7014.603.3TINtoRaster：预先生成分析区数字高程模型TIN数1.50446.7024.905.6据，由于基于栅格的计算比较简单，因此需要将TIN转换2.50446.7073.1016.4成Raster:ArcMap下单击【3DAnalyst】，找到【Convert】3.75446.70181.3040.6的【TINtoRast

8、er】项，Attribute选择Elevation；CellSize3.90446.70205.8046.1输入100；在OutputRaster中输入路径，单击O