swat使用记录教学提纲.docx

1、swat使用记录教学提纲swat使用记录C:UsersAdministratorDocumentsArcGISDefault.gdbReclass_19902C:UsersAdministratorDocumentsArcGISDefault1.gdbReclass_Extr3预处理：要先建立一个模型库一、 DEM数据下载地理空间数据云下载二、土地利用数据准备：土地利用分布图 & 土地利用类型索引表（矢量、栅格）土地利用栅格图（已准备好） 首先进行投影转换Arctoolbox-data management tools-projection and transformation-raster-

2、define projection提取研究区那一部分的土地利用Arctoolbox-extraction-extract by mask重分类土地利用（分成10以内），首先要查询刚刚剪裁完的LUCC数据库，查看其具体分类attributes of 然后再进行重分类，要将剪裁完成的土地利用数据库分类成swat可以识别的数据（查询数据库中的crop）spatial analyst tools-reclass-reclassify（重分类的输出路径可以更改，但是名字好像要取默认的）重分类完成后，填写txt文本索引表。 三、土壤数据库建立（麻烦）C:UsersAd ministratorDocumen

3、tsArcGISDefault1.gdbExtract_chin1HWSD中国土壤数据集土壤数据也需要进行投影转换（同土地利用数据） Data management tools -projections and transformation-raster-define projection提取研究区的土壤 Spaial anylyst tools -extraction-extract by mask打开研究区土壤属性表 Open attribute table查看value打开HWSD.mdb-HWSD-DATA表，按照对应的value,栅格数据的value值与HWSD-DATA表中的MU-

4、GLOBAL字段的数值相对应，寻找出每一种类型的参数。HWSD土壤数据库分为两层，上层用T来表示，下层用S来表示。由于FAO-90 的土壤分级采用是 的土壤分级采用是的USDA分级制，因此不必在对土壤粒径含量进行转换。 打开spaw软件的swc模块，根据上面查出的变量，输入软件运算（以value值为11222为例）输出结果为：TEXTURE、土壤层结构Available water(SQL-AWC)即土壤有效持水量=田间持水量Field Capacity-饱和导水率Wilting pointSat Hydraulic Cond(SQL-K)、饱和导水率/饱和水力传导系数Matric Bulk

5、Density(SQL-BD) 土壤湿密度（mg/m3或g/cm3）另外，在SPAW模型中单位要选择Metric 国际单位制，在options 下拉菜单中选择units下的Metric 即可USLE-K的计算（USLE_K）USLE 方程中土壤侵蚀力因子在利用Williams 等在EPIC模型中发展起来的土壤可蚀性因子K值得估算方法，只需要土壤的有机碳和颗粒组成资料即可计算。其中，fcsand为粗糙沙土质地土壤侵蚀因子；fcl-si为粘壤土土壤侵蚀因子；forgc为土壤有机质因子；fhisand为高沙质土壤侵蚀因子。Sd (sand)为砂粒含量百分数；si(silt) 为粉粒含量百分数；cl(

6、clay) 为粘粒含量百分数；c(T-oc)为有机碳含量百分数。e，高等数学里以自然常数e为底的指数函数.得到土壤分层的大部分参数参数参数含义获取来源SOL-Z1土壤第一层深度HWSD默认为30cm,单位要转化成mmSOL-BD1土壤湿密度运用Spaw软件得出SOL-AWC1土壤有效持水量SOL-K1饱和导水率/饱和水力传导系数SOL-CBN1土壤中有机碳含量（有机质含量乘以0.58）HWSD.mdb-HWSD-DATA表参数查询CLAY1黏土含量SILT1壤土含量SAND1砂土含量ROCK1砾石含量SOL-ALB1地表反射率默认为0.01USLE-K1USLE方程中土壤侵蚀力因子通过计算得出

7、SOL-EC1土壤电导率默认为0土壤数据库（user soil）的其他字段，层数、最大扎根深度、质地可以得到，水文分组可以根据前面算出的SOL-K来判定。MUID、SEQN、S5ID、CMPPCT字段是美国数据特有的，其他地方不需要这些字段的内容，可以参照例子随意填写？ANION-EXCL、SOL-CRK为可选参数，可以参照例子MUID填?SEQN不填SNAM土壤名称在D_SYMBOL90 表中找到土壤数据库的名字SNAMS5ID不填CMPPCT土壤最大可伸缩量默认为1NLAYERS土壤分层数HYDGRP土壤水文学分组（A/B/C/D）选择上层和下层水力传导洗漱中较小的一个，有图判断属于哪个分

8、组SOL-ZMX土壤剖面最大根系深度（mm）对应REF_DEPTH，HWSD为100cm，上层30cm，下层70cm；(注意：此处注意单位为cm，转换为SWAT土壤数据库需要设置为mmANION-EXCL阴离子交换孔隙度模型默认0.5可选SOL-CRK土壤最大可压缩量，以所占总土壤体积的分数来表示模型默认0.5TEXTURE土壤层结构SPAW中可看出T-ECE上层盐度：默认为0或对应T-ECE为0.1注意：1 HYDGRP的确定计算土壤最小下渗率，采用文献试论土壤渗透系数的经验公式和曲线图介绍的计算公式X=(20Y)18其中：X为土壤渗透系数；Y为土壤平均颗粒直径Y=(砂粒含量的百分数/10)

9、*0.03+0.002在D_SYMBOL90 表中找到土壤数据库的名字SNAM,SYMBOL对应 HWSD.mdb-HWSD-DATA表格中SU-SYM90最后将计算结果录入SWAT土壤数据库在Swat-ArcSWAT-Datebases-SWAT2009.mdb 目录中打开usersoil，将做好的土壤数据库中的数据依次填到每一个的后面，结果保存两位有效数字，无参数列一律填0。土壤数据库索引表 solic.txt 的建立。土壤数据库通过土壤索引表来读取数据。inputdata-soil-solic.txt,VALUE 与code 一定要对应.索引表索引了栅格图，Arcgis 读取栅格图，就会

10、读到 name，索引表就会读取到数据库，数据库就把所有参数带入模型进行运算1、 水文站经纬度坐标添加（便于找到流域出口）添加所要选择的水文站坐标右侧工具栏中右键打开的文件：display x y data-Edit-Beijing1954-确定2新建swat project.3子流域划分Watershed delineator -automatic watershed delineatorDem set up加载准备好的dem数据： open dem raster单击 Dem projection setup打开DEM properties 对话框，设置Z unit 为“meter”Burn

11、in-load from disk-选择要加载的河网Stream 定义区被激活1 选择Dem based（基于加载的DEM，使用加载的DEM 自动划分河网和流域2 Pre-defined streams andwatersheds(预定义河网和流域)，需要用户提供河网和子流域数据，并输入到ArcSWAT 中 点击Flow direction and accumulation对DEM 进行填洼，接着计算流向和水流累积量，流向和水流累积量被用来定义河网和计算流域边界。 Stream network-点击create stream and outlets 按钮 Outlet and inlet de

12、finition定义区被激活Add按钮，选择完出口后点stop editingWatershed outlets selection and definition激活Calculation of subbasin parameters激活单击 Calculate subbasin parameters 按钮，计算子流域和河道参数。 在子流域参数计算完成之后，弹出一个提示框。单击确定。流域划分完成，点exit.4、HRU分析 HRU Analysis- Land Use/Soils/Slope New 2005reclassify205(文件夹里的每个文件怎么得到的) Choose grid f

13、ieldvalue 点击ok Lookup table -user table -ok 选择luc.dbf（SWAT土地利用分类会显示在SWAT Land Use Classification Table 中。一旦LandUseSwat代码与所有图中的分类都对应起来了，Reclassify 按钮就被激活了。）点击reclassify2 Soil data3 Slope先在arcgis里面进行坡度分析后再进行该操作单击Multiple Slope 选项，将坡度栅格分为多个坡度分类。Number of Slope Classes 组合框的下方，选择自定义个坡度分类5、HRU定义选择 HRU Ana

14、lysis 菜单下的HRUdefination命令。将低于百分之十的归于其他类中HRU 创建过程中，会产生一个报表。选择HRU Analysis 菜单下的HRU Analysis Reports 可以打开查看这个报表。在列表中，选择Final HRU Distribution，单击OK。流域中创建的HRU 总数以粗体列在报表的顶部。报表的其余部分列出了每一个子流域模拟的土地利用、土壤和坡度，和子流域占流域的百分比、HRU 占子流域的百分比。查看之后，关闭报表。6、气象站单击Write Input Table 下的Weather Stations ，打开Weather Data Locations Table，填入降水测站位置表风站等测站如上单击 OK，以生成这些气象站点的空间图层，并加载这些监测气象数据到SWAT 气象文件中。界面将会自动为流域中的每一个子流域分配不同的气象站。