赵冬泉*1，佟庆远2 ，杜鹏飞1，陈吉宁1

（1.清华大学环境科学与工程系，北京100084 ;

2.北京清华城市规划设计研究院环境与市政所，北京 100084）

摘要：近年来利用GIS技术为手段进行水淹分析一直是一个研究热点，各种文献中报道的分析方法都需要编写程序或通过复杂操作实现。利用ArcGIS 9提供的先进的可视化建模工具，可以根据GIS操作流程快速建立复杂的GIS分析模型。本文利用这种技术建立了水淹分析模型，并在滇池湿地生态规划中进行防浪堤拆除后水淹范围的分析。利用这种方法建立的模型具有使用简单、扩充性强、便于和数据同时发布等特点。该方法适用于各种包含复杂GIS操作的模型构建。

关键字：GIS，空间分析，可视化建模，水淹分析

1. 引言

20世纪90年代以来，利用GIS技术为手段进行水淹分析一直是一个研究热点，许多学者在这方面做过研究，并取得一定的进展。如葛小平等人的GIS支持下的洪水淹没范围模拟[1]，丁志雄等人的基于GIS格网模型的洪水淹没分析方法[2]。这些方法的实现，都需要在GIS环境中编写相应的程序。刘小生等利用ArcMap的自带功能，无需编程实现了低于一特定高程的区域计算[3]。但是该方法的实施需要手工执行多个GIS操作步骤，只有对GIS熟悉的人员才能使用。而在ArcGIS 9提供的统一的地理处理（Geoprocessing）框架下，可以利用可视化建模工具，将数据和相应的空间处理工具按照操作流程快速组织建立复杂的GIS 分析模型。本文将以ArcGIS 9中数据叠加和3D分析处理工具为基础，快速构建一个水淹分析模型，并在滇池湿地生态规划中应用，分析防浪堤拆除后的水淹范围。

2. ArcGIS9中的可视化建模工具介绍

ArcGIS 9是美国环境系统研究所公司（Environmental Systems Research Institute，Inc. 简称ESRI）公司最新推出的GIS平台[4]。ArcGIS 9提供了400多种空间处理工具，包括数据管理、数据转换、矢量分析、栅格分析、3D分析等[5]。同时，在ArcGIS 9中还提供了可视化建模工具（ModelBuilder）。在ModelBuilder中，可以使用多种类型的GIS数据和各种空间处理工具可视化地构建GIS分析模型，从而快速设计和实现复杂地理处理模型的过程，并建立相应的地理处理工作流。另外这个模型可以保存到Geodatabase中，最终将数据和处理模型同时发布，用户可以修改和重新执行这个模型，从而实现了GIS数据和操作方法的统一发布。

*作者简介：赵冬泉，1978.10，博士研究生，汉，甘肃兰州，Email：zdq01@mails.tsinghua.edu.cn，研究方向：GIS在环境规划中应用。

3. 水淹分析模型建立

3.1. 高程数据准备

要进行水淹分析，相关区域的高程数据是必不可少的。在根据地面高程计算给定水位条件下的淹没区时, 应当区分一下两种情形：1) 凡是高程值低于给定水位的点, 皆计入淹没区；2) 需考虑淹没区域的“连通性”，即洪水只淹没它能流到的地方。例如对于环形山（一种中间低洼、四周环形隆起的地形）需要考虑如果水位低于山顶标高, 则只能在山环外形成淹没区。

在ArcGIS中，高程数据可以采用栅格数据结构（GRID）或者不规则三角网（TIN）表达。但是考虑的ArcGIS提供了多种栅格分析的工具，而且栅格数据具有数据结构简单，进行地图代数运算快等特点[6]，因此本文采用的高程数据为栅格数据结构，并将研究区域涉及的地面标高点数据、等高线数据通过TIN模型插值生成GRID文件。

3.2. 模型构建

利用GIS可视化建模工具的关键是制定切实可行的GIS操作流程。在本文中采用了如下所示的GIS操作流程：

图1 水淹分析GIS操作流程

figure 1 GIS operation process of submerge analysis

按照图1所示的操作流程，在ArcCatlog中选定数据所在的Geodatabase文件，新建一个工具箱（ToolBox），然后新建一个GIS操作工具（Tool），按照操作顺序从ArcToolBox中拖拽相应的工具到模型构建窗口中，并建立这些操作和数据之间的对应关系以及操作之间的先后顺序。如图2所示为在ArcGIS中建立的水淹分析模型。

在模型的建立过程中，考虑了淹没区域的“连通性”。模型将首先利用用户指定的水面高度对地面高程图层进行逻辑比较计算，将水面以上的区域设为0，将水面以下的区域设为1；然后只保留其中水面以下的区域，并将该GRID文件转换为多边形矢量文件。接着将防浪堤破损的多边形区域和防浪堤图层进行裁剪（Clip）操作，获得防浪堤破损的确切位置，并与水面以下的区域进行相交运算，只保留与防浪堤破损位置相互连通的淹没区域，即保证了淹没区域的连通性，消除了环形山和孤岛等因素的影响。最后计算淹没区图层中的面积字段，以获取各个淹没区域的面积大小。如果用户还有该区域的土地利用图层，可以扩充这个模型，继续添加相应的GIS操作方法，以获取淹没区的土地利用情况。

图2 在ArcGIS9中建立的水淹分析模型

figure 2 Submerge analysis model in ArcGIS9

4. 实例应用

在滇池湿地生态规划中，需要了解防浪堤拆除对周围环境的影响。利用以上构建的模型，准备相关数据，就可以在ArcGIS 9中直接运行这个模型分析计算防浪堤拆除后的淹没区。在分析过程中，用户不需要了解模型的内部结构，也不需要掌握很多的GIS知识，只要选

择相应的图层、设定相应的参数就可以运行这个模型，如图3所示为模型执行界面。

图3 水淹模型执行界面

figure 3 Submerge model implement interface

在项目中，计算了多个水位高度和不同地区防浪堤拆除的水淹情景，下图所示为在滇池

水位高度为1887.4米的情景下，防浪堤全部拆除和部分拆除后对淹没区域的覆盖范围。

图4 滇池防浪堤拆除淹没区域分析图

figure 4 Submerge analysis map of Dianchi lake with all breakwater removed

5. 结论

传统的水淹模型的实现都需要编写程序实现特定的算法或者通过一系列的手工操作实现。本文利用ArcGIS 9中提供的先进的可视化建模技术，通过一系列的GIS操作实现了水淹区域的分析模型。该模型可以在ArcGIS9的体系下与数据捆绑在一起通过Geodatabase发布，从而实现数据和GIS操作的同时共享。另外建立的模型应用简单方便，普通用户不需要掌握很多GIS知识就可以通过选择数据图层和设置相关参数执行这个模型；同时这个模型也具有良好的扩充性，可以根据淹没区域后期分析的要求在ArcGIS 9提供的ModelBuilder 中添加相应的GIS操作方法。这种模型建立的方法不仅适用于水淹模型，而且也可应用于其他与GIS操作紧密相关的模型建立。

6. 参考文献

[1] 丁志雄, 李纪人, 李琳. 基于ＧＩＳ格网模型的洪水淹没分析方法. 水利学报, 2004, 6(6): 56-61

[2] 葛小平, 许有鹏. GIS支持下的洪水淹没范围模拟. 水科学进展, 2002, 13(2): 456-460.

[3] 刘小生, 黄玉生. 基于Arc/Info的洪水淹没面积的计算方法. 测绘通报, 2003, (6): 46-48

[4] ESRI. What’s New In ArcGIS Desktop 9.0. 2004

[5] ESRI中国（北京）有限公司. What is ArcGIS. 2004

[6] 赵冬泉, 贾海峰, 郭茹等. 图片资料的矢量和栅格处理方法比较. 测绘通报, 2003, (3):

57-58下载本文