【SWAT水文模型】SWAT水文模型建立及应用第一期：DEM数据利用/子流域的划分

WW、forever

已于 2023-08-09 11:44:58 修改

阅读量8.5k

点赞数 24

分类专栏： # 水文模型文章标签：数据库 SWAT水文模型

于 2023-03-12 11:01:06 首次发布

本文链接：https://blog.csdn.net/qq_44246618/article/details/129231698

版权

水文模型专栏收录该内容

35 篇文章 181 订阅

订阅专栏

SWAT水文模型原理及数据库简介可参见另一博客文章-【水文模型】SWAT水文模型原理及数据库简介。
SWAT建模过程中出现的各种错误及解决方案见另一博客-【SWAT水文模型】ArcSWAT各种报错总结。

本博客主要介绍SWAT水文模型建立及应用。以下SWAT水文模型建立过程均以洮河流域为例，区域地理位置/水系图如下：
在这里插入图片描述

SWAT模型建立

SWAT可以用来模拟单个流域或有水力联系的流域系统。每一个流域先被划分为子流域，然后根据土地利用和土壤类型来划分水文响应单元(HRUs)。

模型建立流程如下所示：子流域数据准备→子流域划分→土壤、土地数据准备→气象数据的准备→水文相应单元的划分→率定

步骤1：划分流域边界使用数据：DEM、河流GIS图层
步骤2：HRU子流域划分使用数据：土地利用库、土壤库、气象库
步骤3：率定（SWAT CUP）使用数据：水文站（径流数据）

常见疑问

1.研究区域并非完整流域，如何处理？
我的做法是直接将包含该区域的三级流域拿来分析。

1 DEM数据的利用

DEM (Digital Elevation Model) 是描述区域地貌形态空间分布的GIS数据。DEM用一组有序数值阵列形式表示地面高程，是一定范围内规则格网点的平面坐标(X，Y)及其高程(Z)的数据集。在利用SWAT模型建模的过程中，将利用DEM数据进行流域 水系生成 以及 子流域划分 等工作。

1.1 高分辨率DEM数据的获取

目前使用较多的高分辨率DEM数据包括 SRTM（Shuttle Radar Topography Mission） 90米分辨率数据以及ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) GDEM（Global DEM）30米分辨率数据两类。这两类数据都可通过中国科学院国际科学数据服务平台进行下载。
在这里插入图片描述
此外，亦有ALOS_PALSAR_12.5m分辨率DEM数据。

1.2 数据处理

1）DEM数据的下载与拼接处理

本次处理以ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) GDEM（Global DEM）30米分辨率数据为例，根据经纬度1°×1°分片，拼接裁剪得到研究区域DEM数据，具体拼接裁剪过程见另一博客-【ArcGIS】拼接/裁剪得到研究区域DEM高程图。
研究区域内30m分辨率DEM数据如下所示：
在这里插入图片描述

2）数据的投影转换

GIS数据以经纬度或者平面坐标（X,Y坐标系）保存位置信息。ArcSWAT要求所有输入GIS数据具有平面坐标。对于具有经纬度坐标系统的GIS数据，需要将经纬度坐标转化为平面坐标，即进行投影转换。
如若导入投影不对，会弹出以下界面：
在这里插入图片描述
在实际应用中需要根据GIS数据的具体情况选择合适平面坐标系。在选择投影坐标系后点击应用，即将数据采用选定的投影进行显示。本次研究流域采用WGS_1984_UTM_Zone_46N 投影坐标系。

步骤1：设置图层投影坐标

将DEM数据在ArcMap中打开，查看图层（Layer）的属性（Properties）中的坐标系（Coordinate System）信息，如下图所示。该对话框上部显示当前该图层的坐标系统为WGS_1984_UTM_Zone_46N 投影坐标。
在这里插入图片描述

步骤2：更改DEM数据投影

但以上操作并未实际改变DEM数据的投影。需要将数据按照选定投影进行输出。右键点击图层列表中的DEM数据，选择Export Raster Data，如下图所示。
在这里插入图片描述
左上部的输出范围（Extent）选项选择 Raster Dataset，按DEM原始数据范围进行输出。空间参照系（Spatial Reference）选择 Data Frame，即按当前所使用的坐标系统进行输出。在输入生成文件的路径信息后点击保存(Save)，则生成的数据即为转化为平面坐标系的DEM数据。

2 子流域的划分

子流域划分的关键步骤 ，包括：

加载DEM
(可选)添加掩膜指定研究区域(加载或手绘Mask)
(可选)加载河网水系
DEM处理
指定最小子流域面积(临界阈值)
编辑河网节点
计算子流域参数
(可选) 设置水库

2.1 准备工作：启动ArcSWAT分区模块

启动ArcMap，打开一个空文档，在【Customize】菜单中，点击Extensions，激活至少以下四个扩展模块：Spatial Analyst、SAWT Project Manager、SWAT Watershed Delineator、Data Interoperability。其他模块根据所需进行激活。
在这里插入图片描述
在View菜单中，单击Toolbars并选择ArcSWAT工具条，ArcSWAT工具条出现在ArcMap窗口中。如下所示：

2.2 定义SWAT工程

在SWAT Project Setup 下拉菜单中选择建立工程，在如下图所示对话框中选择工程存放路径及名称，选择读取数据库的位置。
在这里插入图片描述
在Watershed Delineation下拉菜单中选择Watershed Delineation。进入如下图所示流域划分对话框，该对话框分为5个部分：DEM设置，河网定义，0utlet，Inlet定义，流域总出口选择以及子流域参数的计算。

如图需要注意标注的几个方面：

步骤1：DEM的加载

DEM通过Load from Disk，路径选择前文保存的DEM数据。
在这里插入图片描述

加载DEM数据之后，单击DEM projection setup按钮定义DEM属性，如下图，查看并确认信息是否准。Z Unit选择为meter。
在这里插入图片描述

步骤2：MASK定义（可选）

说明：若加载的DEM大小已剪裁好，则可不使用Mask。

加载DEM之后，就可以选择加载Mask，其主要作用是：准确地确定研究区域，减少处理DEM的数据量。如下图所示，共有三个选项：
（1）从硬盘加载；
（2）从图中选择；
（3）手动绘制。手绘Mask完成之后，点击右键选择Stop Editing即可。
在这里插入图片描述

步骤3：数字化河网添加（可选）

数字化河网通过Burn In进行添加，可对河网及子流域生成进行“引导”。
加载河网的好处在于，可以更好的生成与实际较符合的河网水系，尤其在河流下游的平坦区域。可根据情况，添加“水系”所在文件路径。
在这里插入图片描述

步骤4：河网生成

一种是基于DEM（DEM-based），这种方法需要先进行Flow direction and accumulation
计算完成后下方Area会出现一个范围及默认值，Area指的是Upstream drainage area即汇水面积，实际上是形成河流的最小汇水面积，设置该值会影响河网的详细程度和最终形成的子流域大小，取值越小，河网将会越详细，子流域面积越小。

单击按钮，进行河网的生成和计算。
在这里插入图片描述
河网计算完毕后，会提示【End of DEM grid preprocessing】。

Flow direction and accumulation完成后，对生成的Area阈值可根据与实际河网的匹配进行修改，该值越小，划分的河网就会越详细。本操作中修改阈值Area为15000，点击按钮：生成新的河网。
在这里插入图片描述
运行完成后，可发现，河网已经成功生成，提示【End of DEM grid preprocessing】，如下图所示，

步骤5： Outlet、Inlet定义（可选）

流域内将要用来做径流、泥沙和水质校准的水文站点、水质站点，可以定义为OUTLET，以便更好地定位监测点的位置，以作校准、验证之用。在这里也可以直接指定流域总出口的位置。注意这些站点必须是临近河边的测站。也可以手动编辑所有的Outlets、Inlets，单击ADD按钮可以直接在水系上添加一个节点，也可以单击 DELETE删除节点。编辑完成之后，单击右键选择Stop Editing。

如果选择Add by Table进行流域出口的加载，打开如下所示窗口：
在这里插入图片描述
按手册制作Outlet表：

将此文件保存为.dbf格式后，可成功导入在SWAT中：
在这里插入图片描述
确定后，显示以下界面：（成功导入14个outlet）

如果研究区域有水库需要添加，在后期定位时只能定位在最临近的outlet 上，所以需提前利用Subbasin outlet 为水库提前添加节点（如下图所示），以便后期选择水库定位。
在这里插入图片描述
说明：SWAT2009 增加了自动为子流域添加预留点源的功能，选中“Add pointsource to each subbasin”，即可自动在每段子流域中设定点源排污口位置。如若在SWAT2005 中为每条子流域添加点源排污口，需要预先做完watershed Delineation划分完子流域后，记下位置再返回来添加点源排污口，或通过“Add by table”直接输入文件，确定点源排污口位置。

步骤6：流域总出口指定及子流域划分

单击按钮，选择全流域outlet，可以选一个或多个（按住shift 进行复选），按钮可以取消所作选择。
在这里插入图片描述
本研究区域出口有两个，选择如下：
（如若只需要其中一个子流域，可只选取其中一个出口即可）

确定后按Delineation watershed 按钮，开始生成子流域。

完成后，会提示【Watershed delineation is done.】，如下图：
在这里插入图片描述
需要注意的是，流域边界并不与输入的DEM边界相同，存在些许差异，如下图所示：

步骤7：子流域参数的计算

单击Calculate subbasin parameter 按钮，进行子流域各参数的计算。
在这里插入图片描述
计算完毕出现计算成功字样。

此时，GIS界面如下：

此结果明显子流域过多，可重新设置阈值分类。结果如下：

如若有水库，可通过按钮，可添加流域内的水库，该水库位置会自动吸附在附近的节点处。
在这里插入图片描述
子流域生成完毕，按Exit键，退出并保存工程。

总结： 当流域划分完成之后，单击Exit按钮退出划分过程。ArcSWAT产生的栅格数据集，将从SWAT项目目录Watershed\Grid转移到Project Raster Geodatabase。流域划分完成之前，Watershed\Grid目录中的栅格以ESRI GRID格式存储，以提高执行效率。一旦划分完成，它们将会被转移到Raster Geodatabase，以简化项目的数据存储。