最终篇——SWAT模型的建立

基于以上数据，我们可以开始进行SWAT模型的计算了。（SWAT模型只要把数据处理完，计算基本上就是比较简单了。）

1. 1. 建立SWAT项目

SWAT模型可以模拟的项目主要有：

• 天气发生器（可以通过历史天气统计数据来进行模拟阶段的天气随机生成，但非常不准确）
• 水文循环（包括地表和地下水）
• 氮磷和农药迁移转化
• 土壤侵蚀
• 植物生长（对水量和N,P的消耗，生物量积累）
• 河流水文水质过程（河流水文为稳态模型，只输出流量；水质为QUALE2模型）
• 水体（湿地，水塘，湖泊，水库：但是只是水平衡方程，非常简化）
• 管理措施（施肥，耕作方式，灌溉等）

建立模型的初步需要的文件有：流域文件、数字高程文件、土壤文件与土地利用文件

文件建立时，要新建文件夹，选择SWAT Project Setup  New SWAT Project,双击打开文件夹窗口。文件夹储存位置，千万要选择英文名称，否则文件将会创建失败（所有件夹不要出现中文名称）。SWAT文件库会自动从安装文件中读取，然后建立在所储存的文件中，找到这个文件，替换成我们之前修改过的文件。

1. 1. SWAT模型子流域划分

SWAT文件建立，并且将数据加载进来后，点击Automatic Watershed Delineation,出现窗口。

选取文件夹，点击Select from map,即从我们之前导入的DEM中选择。也可以选择从文件中读取，但是这样读取文件非常容易报错，不建议采用。选择图层后，点击OK即可。

 

 

文件读取成功后，显示如下。点击DEM project ，打开DEM project窗口，将Z Unit转化成米，即meter。

随后，选择Burn，用同样的方式将流域导入。导入成功的界面如下，点击Flow direction and accumulation 进行数据计算。计算出来的数值为Area区域，显示数字为195.2172，这里的数字表示每195.2172Ha为一个小计算单元，共有2169个计算单元。Area 计算的195.2172数值可以进行手动修改，我们这里将这个数值填写为100.

计算完后点击Stream network 下方的Create streams and outlects根据DEM进行计算。

在Outlet and Inlet Definition下方选择Point source input ,点击下方ADD，添加点源。随后再Watershed Outlects(s) Selection下方选择添加流域出口，点击

进行框选，流域出口选择完毕后，点击Delineating watershed 进行流域划分自动计算，计算时间较长，耐心等待计算结果。点击OK计算结束，发现研究区域自动划分了流域与径流，左侧也出现相应的图层清单。

计算完后，回到工具选项卡，点击Calculate subbasin in parameters，进行流域计算，计算结束后，点击Exit退出。保存的数据储存在我们一开始建立的文件夹里。

1. 1. HRU划分

子流域划分结束后，我们对划分的子流域模块建立HRU模块，选择SWAT工具条的HRU Analysis Land Use/Soils/Slope Definition，定义HRU模型。

打开Land Use/Soils/Slope Definition后出现窗口，我们在Land Use Grid 选择我们下载的栅格数据（切割过的），通过图层选择（把栅格图拖至Arcgis图层后，从图层选择），Choose Grid Field 下面选择Value，点击OK。点击LookUp Table，找到我们之前建立的土地利用数据的索引表，（我之前用的Arcgis10.2中文版，导致索引表一直无法成功导入，将中文文件删除或关闭解决了这个问题），也可以手动选择输入数据，数据导入成功后，点击Reclassify进行重分类。

 

同理，Soil Data也是同样的方法导入，在Soil Database Options的选项中，选择UserSiol，用同样的方法导入土壤类型，进行重分类。坡度Slope的定义，我们可以选用单一坡度与多坡度，这里我们的研究区域坡降较小，选择单一坡度，即Single Slope。若选择多坡降，可以在Number of Slope选择分级量，在Current Slope Class选择编辑的级别（选择1、2，取决于你设定的分级量），选择高程范围。坡度设置完成后，点击Reclassify进行重分类，最后勾选Create HRU Feature Class，点击 overlay（下方图片勾选错误，运行的时候注意一下，不是点击cancel，而是点击overlay）生成HRU单元。

 

将土地利用数据、土壤数据和坡度设置完后，点击选择HRU Definnition，定义HRU的计算方式。打开窗口后，只定义HRU Thresholdes，其余两个不用设置，设定各数据单元的计算率，这里的10%的意思是计算精度，即分成级HRU计算单元太小，小于平均计算单元的10%，我们就不对这HRU单元进行计算（也就是说太破碎的计算单元不进行计算），若选择0%意味着全部进行计算，需要耗费大量的内存和时间，对计算机要求较高，不建议。

定义完后，选择Create HRUs，这个过程较长，耐心等待。

1. 1. 气象数据及其他数据输入

HRU单元划分结束后，我们进入天气发生器的编辑。点击Write Input Table      Weather Stations打开天气数据编辑窗口，选择Weather Generator Data     WGEN_user，其他模块的选择，导入我们之前做的索引表，点击OK。

 

索引表的对应如下：Rainfall Data     pcp，Solar Radiation Data       sor，Wind Speed Data     wd，Relative Humidity   rh，Temperature Data      tm。

天气数据导入成功后，选择Weather Input Table Write SWAT Input Tables，打开天气发生器，未写入信息时，数据显示红色的Incomplete字体。点击Select      Create Tables，会发现表格中一个个Incomplete会生成绿色的complete，点击OK，天气发生器建立完成。

1. 1. 输入数据的编辑

我们在之前建立SWAT项目中加入了一个点源，即烔炀镇污水处理厂。点击Edit SWAT Input      Point Source Discharges进行编辑数据，假如我们在子流域划分阶段添加了水库或管道等信息，则在Edit SWAT Input     Reservoirs或Inlet Discharges编辑水库或城市管道信息。

以点源信息输入为例，打开窗口，选择编辑的点源编号点击OK开始编辑。打开编辑窗口，点击Edit Values，选择点源输入的数据类型，保存编辑，点击EXIT退出。

点源数据的输入包括流速、有机氮、有机磷、硝态氮和亚硝酸盐、溶解氧等。可以输入连续性数据，也可以输入日变化数据、年变化数据、月变化数据等。输入随时间变化的点源污染物排放数据类型，需要在Observed Loadings Input Files下方选择数据文件夹。

若点源导入时间序列所创立的文件夹形式如下，这里显示每一天的数据。

ARCSWAT21中文手册.pdf-原创力文档这是网站信息。

1. 1. SWAT运行及结果读取

SWAT模型的运行选择SWAT选项卡SWAT Simulation Run SWAT，设置模拟的时间，输出的数据类型，设置完后，点击Setup SWAT   Run SWAT，运行完毕。

模型运行成功后，将数据进行导出。选择我们需要导出的数据，通常选择前三个，分别为河道数据（.rch）、子流域数据（.sub）、水文响应单元数据（.hru），选择完毕后，点击Import Files to Database      OK，填写储存的名称后，点击Save simulation。最后，数据的查看在我们之前建立SWAT创建的文件夹中。（例如，我创立的文件夹为EFDC.mlb）

我保存的文件名字为simulate20240104，我们在之前保存的文件夹中找到Scenarios这个文件夹，打开，找到我们保存的simulate20240104     TablesOut打开Access，如下所示：

HRU输出文件包含流域中每个水文响应单元的摘要信息。该文件以电子表格格式编写。以下是对HRU输出文件中的输出变量的简要描述。（在SWAT操作手册中有描述swat-io-documentation-2012.pdf，在操作手册中的第32章，详细描述了各字段意思）