SWAT水文模型.docx

1、SWAT 水文模型介绍1 概述SWAT（Soil and Water Assessment Tool）模型是美国农业部（USDA）农业研究局（ARS）开发的基于流域尺度的一个长时段的分布式流域水文模型。它主要基于 SWRRB 模型，并吸取了 CREAMS、GLEAMS 、EPIC 和 ROTO 的主要特征。SWAT 具有很强的物理基础，能够利用 GIS 和 RS 提供的空间数据信息模拟地表水和地下水的水量和水质，用来协助水资源管理，即预测和评估流域内水、泥沙和农业化学品管理所产生的影响。该模型主要用于长期预测，对单一洪水事件的演算能力不强，模型主要由 8 个部分组成：水文、气象、泥沙、土壤温度

2、、作物生长、营养物、农业管理和杀虫剂。SWAT 模型拥有参数自动率定模块，其采用的是 Q.Y.Duan 等在 1992 年提出的 SCE-UA 算法。模型采用模块化编程，由各水文计算模块实现各水文过程模拟功能，其源代码公开，方便用户对模型的改进和维护。2 模型原理SWAT 模型在进行模拟时，首先根据 DEM 把流域划分为一定数目的子流域，子流域划分的大小可以根据定义形成河流所需要的最小集水区面积来调整，还可以通过增减子流域出口数量进行进一步调整。然后在每一个子流域内再划分为水文响应单元 HRU。HRU 是同一个子流域内有着相同土地利用类型和土壤类型的区域。每一个水文响应单元内的水平衡是基于降水

3、、地表径流、蒸散发、壤中流、渗透、地下水回流和河道运移损失来计算的。地表径流估算一般采用SCS 径流曲线法。渗透模块采用存储演算方法，并结合裂隙流模型来预测通过每一个土壤层的流量，一旦水渗透到根区底层以下则成为地下水或产生回流。在土壤剖面中壤中流的计算与渗透同时进行。每一层土壤中的壤中流采用动力蓄水水库来模拟。河道中流量演算采用变动存储系数法或马斯金根演算法。模型中提供了三种估算潜在蒸散发量的计算方法Hargreaves、Priestley-Taylor 和Penman-Monteith。每一个子流域内侵蚀和泥沙量的估算采用改进的 USLE 方程，河道内泥沙演算采用改进的 Bagnold 泥沙

4、运移方程。植物吸收的氮采用供需方法计算，植物的氮日需求量是植物与生物量中氮浓度的函数。土壤中向植物供给氮，当需求超过供给时，出现营养物压力。地表径流、壤中流和渗透过程运移的硝态氮量由水量和土壤层中的平均硝态氮浓度来估计。泥沙中运移的有机氮采用 McElroy et al.开发的负荷方程，后经进一步改进。该负荷方程基于土壤表层的有机氮浓度、泥沙量和富集率来估计径流中的有机氮损失。植物吸引的磷采用与氮相似的供需方法。径流中带走的可溶解磷采用土壤表层中的不稳定磷、径流量和磷土分离系数来计算。泥沙运移的磷采用与有机氮运移相同的方程。河道中营养物的动态模拟采用 QUAL2E 模型。2.1 产汇流模型模型

5、中采用的水量平衡表达式为：=0+=1()式中： 为土壤最终含水量，mm； 为土壤前期含水量，mm；t 为时间 0步长，d； 为第 i 天降水量，mm； 为第 i 天的地表径流，mm； 为第 i 天的蒸发量，mm ； 为第 i 天存在于土壤剖面底层的渗透量和侧流量，mm； 为第 i 天地下水含量， mm。SWAT 模型水文循环陆地阶段主要有水文、天气、沉积、土壤温度、作物产量、营养物质和农业管理等部分组成。模型径流量产生计算流程图如下图 2-1 所示。图 2-1 SWAT 模型产流计算流程图2.1.1 地表径流当落到地表的降水量多余入渗量时产生地表径流。SWAT 模型采用 SCS 径流曲线法计算

6、。SCS 曲线方程自上世纪 50 年代逐渐得到广泛使用，属于经验模型，是对全美小流域降水与径流关系 20 多年的研究成果。模型能反应不同土壤类型和土地利用方式及前期土壤含水量对降雨径流的影响，它是基于流域的实际入渗量（F）与实际径流量（Q）之比等于流域该场降雨前的最大可能入渗量（S ）与最大可能径流量（ ）之比的假定基础上建立的SCS 模型的降雨- 径流基本关系表达式如下：=式中：假定潜在径流量（ ）为降水量（P）与由径流产生前植物截留、初渗和填洼蓄水构成的流域初损（ ）的差值。由此推到上式有：=()2+初损 受土地利用、耕作方式、灌溉条件、冠层截留、下渗、填洼等因素的影响，它与土壤最大可能入

7、渗量 S 呈一定的正比关系，美国农业部土壤保持局在分析了大量长期的实验结果基础上，提出了二者最合适的比例系数为 0.2，即： =0.2由此可得 SCS 方程为：=(0.2)2+0.8 0.2=0 0.2流域当时最大可能滞留量 在空间上与土地利用方式、土壤类型和坡度等下垫面因素密切相关，模型引入的 值可较好地确定 ，公式如下： =25400254CN 是一个无量纲参数， CN 值是反映降雨前期流域特征的一个综合参数，它是前期土壤湿度、坡度、土地利用方式和土壤类型状况等因素的综合。2.1.2 蒸散发模型考虑的蒸散发是指所有地表水转化为水蒸气的过程，包括树冠截留的水分蒸发、蒸腾和升华及土壤水的蒸发。

8、蒸散发是水分转移出流域的主要途径，在许多江河流域，蒸发量都大于径流量。准确地评价蒸散发量是估算水资源量的关键，也是研究气候和土地覆盖变化对河川径流影响的关键。（1）潜在蒸散发模型提供了 Penman-Monteith、Priestley-Taylor 和 Hargreaves 三种计算潜在蒸散发的方法，另外还可以使用实测资料或已经计算好的逐日潜在蒸散发资料。一般采用 Penman-Monteith 方法来计算流域的潜在蒸散发。（2）实际蒸散发实际蒸散发以潜在蒸散发为计算基础。在计算流域实际蒸散发量的时候，模型首先计算植物冠层截留水分的蒸发，然后计算最大蒸腾量、最大升华量和最大土壤蒸发量，最后计

9、算实际的升华量和土壤水分蒸发量。（3）冠层截留蒸发量模型在计算实际蒸发时假定尽可能蒸发冠层截留的水分，如果潜在蒸发量小于冠层截留的自由水量 ，则：0 =0()=()式中： 为某日流域的实际蒸发量，mm； 为某日冠层自由水蒸发量， mm； 为某日的潜在蒸发量，mm ； 为某日植被冠层自由水初始含量，0 ()mm； 为某日植被冠层自由水终止含量，mm。()如果潜在蒸发量 大于冠层截留的自由水含量 则：0 =()()=0当植被冠层截留的自由水被全部蒸发掉，继续蒸发所需的水分就会从植被和土壤中得到。（4）植物蒸腾假设植物生长在一个理想的条件下，植物蒸腾可用以下表达式计算：当 0LAI3.0 时=03.

10、0当 LAI3.0 时=0式中： 为某日最大蒸腾量，mm； 为植被冠层自由水蒸发调整后的潜在蒸发 0，mm ；LAI 为叶面积指数。因为没有考虑到植物下面图层的含水0=0量问题，由此公式计算处的蒸腾量可能比实际蒸腾量要大一些。（5）土壤水分蒸发在计算土壤水分蒸发时，首先区分出不同深度土壤层所需要的蒸发量，土壤深度层次的划分决定土壤允许的最大蒸发量，可由下式计算：， = +(2.3470.00713)式中： 为 z 深度处蒸发需要的水量，mm ；z 为地表以下土壤的深度， mm。表达式中的系数是为了满足 50%的蒸发所需水分来自土壤表层 10mm，以及 95%的蒸发所需水分来自 0100mm 土

11、壤深度范围内。土壤水分蒸发所需要的水量是有土壤上层蒸发需水量与土壤下层蒸发需水量决定的： ， =， ， 式中： 为 ly 层的蒸发需水量，mm； 为土壤下层的蒸发需水量， ， mm； 为土壤上层的蒸发需水量，mm。， 土壤深度的划分假设 50%的蒸发需水量由 010mm 内土壤上层的含水量提供，因此 100mm 的蒸发需水量中 50mm 都要由 10mm 的上层土壤提供，显然上层无法满足需要，这就需要建立一个系数来调整土壤层深度的划分，以满足蒸发需水量，调整后的公式可以表示为： ， =， ， 式中：esco 为土壤蒸发调节系数，该系数是 SWAT 为调整土壤因毛细作用和土壤裂隙等因素对不同土层

12、蒸发量二提出的，对于不同的 esco 值对应着相应的土壤层划分深度。2.1.3 土壤水渗入到土壤中的水有多种不同运动方式。土壤水可以被植物吸收或蒸腾而损耗，可以渗透到土壤底层最终补给地下水，也可以在地表形成径流，即壤中流。由于主要考虑径流量的多少，因此对壤中流的计算简要概括。模型采用动力储水方法计算壤中流。相对饱和区厚度 计算公式为：00=2， 1000式中： 为土壤饱和区内可流出的水量，mm； 为山坡坡长，m；， 为土壤可出流的孔隙率； 表示土壤层总孔隙度，即 与土壤层水分含 soil量达到田间持水量的孔隙度 之差。fcd=soil-fc山坡出口断面的净水量为：Qlat=24H0vlat式中

13、： 为出口断面处的流速，mm/h。其表达式为：vlatvlat=Ksatslp式中： 为土壤饱和导水率，mm/h； 为坡度。Ksat slp总结上面表达式，模型中壤中流最终计算公式为：Qlat=0.0242SWly， excessKsatslpLhill2.1.4 地下水模型采用以下表达式来计算流域地下水：Qgw， i=Qgw， i-1exp(-gwt)+wrchrg1-exp(-gwt)式中： 为第 i 天进入河道的地下水补给量，mm； 为第（i-1）天Qgw， i Qgw， i-1进入河道的地下水补给量，mm； 为时间步长，d； 为第 i 天蓄水层的补t wrchrg给流量，mm； 为基流

14、的退水系数。其中补给流量由下式计算：gwwrchrg， i=1-exp(-1/gw)Wseep+exp(-1/gw)wrchrg， i-1式中： 为第 i 天蓄水层补给量，mm； 为补给滞后时间，d； 为wrchrg， i gw Wseep第 i 天通过土壤剖面底部进入地下含水层的水分通量，mm/d； 为第wrchrg， i-1（i-1）天蓄水层补给量，mm 。2.2 土壤侵蚀模型泥沙生成量用 MUSLE 方程来预测，计算渠道泥沙输移量的公式为：T=aVb式中： 为输移能力，t/m 3； 为流速，m/s；a 和 b 是常数。T V根据天气条件，泥沙输移量可以高于或者低于输移能力，导致沉积过量的

15、泥沙通过渠道侵蚀再悬浮输移泥沙。流速方程为：V= Fwd式中： 为流量，m 3/s； 为渠道宽度，m； 为径流深，m。F w d对于低于齐岸深度的径流，径流深使用 Manning 方程来计算，假定渠道宽度远大于深度：d=( Fnwcs0.5)0.6式中： 为渠道曼宁系数； 为渠道坡度，m/m。n cs由于降水和径流产生的土壤侵蚀是用 MUSLE 方程来计算的，MUSLE 是修正的通用土壤流失方程（USLE） 。USLE 方程是通过降水动能函数预测年均侵蚀量，而在 MUSLE 中，用径流因子代替降水动能，改善了泥沙产量的预测，这样就不需要泥沙输移系数，并且可以将方程用于单次暴雨事件，因为径流因子

16、是先行湿度和降水动能的函数。USLE 中需要输移系数是因为降水动能因子表示的能量只在作用流域内起作用。修正的通用土壤流失方程为;msed=11.8(QsurfqpeakAhru)0.56KUSLECUSLEPUSLELSUSLEFRG式中： 为土壤侵蚀量，t； 为地表径流，mm/h； 为洪峰流量，msed Qsurf qpeakm3/s； 为水文响应单元的面积，hm 2； 为土壤侵蚀因子； 为植被覆Ahru KUSLE CUSLE盖和管理因子； 为水土保持措施因子； 为地形因子； 为粗碎屑PUSLE LSUSLE FRG因子。2.2.1 土壤侵蚀因子 KUSLE当其他影响侵蚀的因子不变时， 因

17、子反映不同类型土壤抵抗侵蚀力的高低。K它与土壤物理性质的影响，如机械组成、有机质含量、土壤结构、土壤渗透性等有关。当土壤颗粒粗、渗透性大时， 值就低，反之则高；一般情况下 值得K K变幅在 0.02 0.75 之间。值得直接测定方法是：在标准小区（坡长为 22.1m，宽为 1.83m，坡度为 9%）K上没有任何植被，完全休闲，无水土保持措施，降水后收集由于坡面径流而冲蚀到集流槽内的土壤，烘干、称重，由公式计算得到 值。K试验测算 值既费时又费力，1971 年一些学者发展了一个通用方程来计算K土壤侵蚀因子 值，该方程在土壤黏土和壤土组成少于 70%时适用。KKUSLE=0.00021M1.14(

18、12-OM)+3.25(csoilstr-2)+2.5(cperm3)100式中： 为颗粒尺度参数； 为有机物含量百分比，%； 为土壤分类中的M OM csoilstr结构代码； 为土壤剖面可渗透性类别。cperm1995 年另一位学者提出了一个替换方程： =式中： 为粗糙砂土质地土壤侵蚀因子； 为粘壤土土壤侵蚀因子； 为土壤有机质因子； 为高砂质土壤侵蚀因子。各因子的计算公式如 下：=0.2+0.30.256(1100)=( +)0.3=1 0.25+(3.722.95)=1 0.7(1100)(1100)+5.51+22.9(1100)式中： 为粒径在 0.052.00mm 沙粒的百分含量

19、； 为粒径在 0.002 0.05mm 的淤泥、细砂百分含量； 为粒径小于 0.002mm 的粘土百分含量； 为各土壤层中有机碳含量，% 。2.2.2 植被覆盖因子 CUSLE植被覆盖和管理因子 表示植物覆盖和作物栽培措施对防止土壤侵蚀的CUSLE综合效益，其含义是在地形、土壤、降水条件相同的情况下，种植作物或林草地的土地与连续休养地土壤流失量的比值，最大取值为 1.0.由于植被覆盖受植物生长期的影响，SWAT 模型通过下面的方程调整植被覆盖和管理因子 ：CUSLECUSLE=0.8(CUSLE， )(0.00115)+(CUSLE， )式中： 是最小植被覆盖和管理因子值； 是地表植物残留量，

20、CUSLE， kg/hm2。最小 C 因子可以由已知年平均 C 值，通过以下方程计算。CUSLE， =1.463(CUSLE， )+0.1034式中： 表示不同植被覆盖的年均 C 值。CUSLE， 2.2.3 水土保持措施因子 PUSLE水土保持因子 是指有保持措施的地表土壤流失与不采取任何措施的地PUSLE表土壤流失的比值，这里的水土保持措施包括等高耕作、带状种植和梯田。等高耕作对于中低强度的降水侵蚀具有保护水土流失的作用，但对于高强度的降水其保护作用则很小，等高耕作对坡度为 3%8%之间的土地非常有效。2.2.4 地形因子地形因子 的计算公式如下：LSUSLELSUSLE=(22.1)(6

21、5.412+4.56sin+0.065)式中： 为坡长； 为坡长指数； 为坡度。 坡长指数 的计算公式如下：=0.61(35.835)式中： 为水文相应单元的坡度， 。 =tan2.2.5 因子 FRG因子是通过下面公式计算的：FRGFRG=(0.053)式中： 为第一层土壤中砾石的百分比，% 。2.3 污染负荷模型SWAT 模型可以模拟不同形态氮的迁移转化过程，地表径流流失、入渗淋失、化肥输入等物理过程，有机氮矿化、反硝化等化学过程以及作物吸收等生物过程，氮可以分为有机氮、作物氮和硝酸盐氮三种化学状态，氮的生物固定、有机氮向无机氮的转化以及溶解性氮随侧向壤中流的迁移等过程，有机氮又被划分为活

22、泼有机氮和惰性有机氮两种状态，以及铵态氮挥发过程的模拟。2.3.1 硝态氮污染负荷模型硝态氮主要随地表径流、侧向流或渗流在水体中迁移，要计算随水体迁移的硝态氮量必须先计算自由水中硝态氮浓度，用这个浓度乘以各个水路流动水的总量，即可得到从土壤中流失的硝态氮总量。自由水部分的硝态氮浓度可用下面公式计算：3,=3 (1)式中： 为自由水中硝态氮浓度，kg/mm； 为土壤中硝态氮的量，3, 3kg/hm2； 为土壤中自由水的量，mm； 为孔隙度； 为土壤饱和含 水量。（1）通过地表径流流失的溶解态氮计算公式：3=33， 式中： 为通过地表径流流失的硝态氮， kg/hm2； 为硝态氮渗流系3 3数； 为

23、地表径流，mm。 （2）通过侧向流流失的溶解态氮的量计算公式：对于地表 10mm 图层：3， =33， ， 对于 10mm 以下的土层：3， =3， ， 式中： 为通过侧向流流失的硝态氮， kg/hm2； 为硝态氮渗流系数；3， 3为自由水的硝态氮浓度，kg/mm； 为侧向流，mm。3， ， （3）通过渗流流失的溶解态氮量计算公式：3， =3， ， 式中： 为通过渗流流失的硝态氮， kg/hm2； 为自由水的硝3， 3， 态氮浓度，kg/mm； 为渗流，mm。， 2.3.2 有机氮污染负荷模型有机氮通常是吸附在土壤颗粒上随径流迁移的，这种形式的氮负荷与土壤流失量密切相关，土壤流失量直接反映了有

24、机氮负荷。有机氮随土壤流失的输移量计算公式为：=0.001 式中： 为有机氮流失量，kg/hm 2； 为有机氮在表层土壤中的浓度， kg/t； 为土壤流失量，t； 为水文相应单元的面积，hm 2； 为氮富集系数， 氮富集系数是随土壤流失的有机氮浓度和土壤表层有机氮浓度的比值。计算富集系数的公式如下：=0.78()0.2468式中： 为地表径流中泥沙含量。 的计算公式如下： = 10式中： 为土壤流失量，t ； 为水文相应单元面积，hm 2； 为地表径流， mm。2.3.3 溶解态磷污染负荷模型溶解态磷在土壤中的迁移主要是通过扩散作用实现的，扩散是指离子在微小尺度下由于浓度梯度而引起的溶质迁移，

25、由于溶解态磷不很活跃，所以由地表径流以溶解态形式带走的土壤表层的磷很少，地表径流输移的溶解态磷可由下面公式计算：=， ， 式中： 为通过地表径流流失的溶解态磷，kg/hm 2； 为土壤中 ， 溶解态磷，kg/hm 2； 为土壤溶质密度， mg/m3； 为表层土壤深度，mm； 为土壤磷分配系数，表层土壤中溶解态磷的浓度和地表径流中溶解态磷， 浓度的比值。2.3.4 有机磷和矿物质磷污染负荷模型有机磷和矿物质磷通常是吸附在土壤颗粒上通过径流迁移的，这种形式的磷负荷与土壤流失量密切相关，土壤流失量直接反映了有机磷和矿物质磷负荷，有机磷和矿物质磷随土壤流失输移量计算公式为：=0.001 式中： 为有机

26、磷流失量， kg/hm2； 为有机磷在表层土壤中的浓度， kg/t； 为土壤流失量，t； 为水文相应单元的面积，hm 2； 为磷富集系数。 2.4 河道演算模型模型中定义的河道均是明渠流，SWAT 用曼宁公式来定义河道糙率和水流流速。水流在河道中演进过程使用变储量演算法或马斯金根法，两种方法都是动力波方程。在模拟中一般使用马斯金根法来计算。（1）马斯金根法马斯金根法假设河道内水体形状是由一个楔形蓄水体和一个棱形蓄水体组成，如下图 2-2 所示。图 2-2 河道水体示意图当洪水波行进到某个河段槽，入流量大于出流量便形成了楔形蓄水体。当洪水波退去，在河段槽便出现了出流量大于入流量的负楔蓄体。另外对

27、于楔蓄水体，河段槽内始终包含一个体积为流域长度上横截面不变的棱柱状水体。总的蓄水容量为：Vstored=Kqout+KX(qin-qout)式中： 为蓄水容量，m 3； 为入流量，m 3/s； 为出流量，m 3/s； 为Vstored qin qout K稳定流情况下的河段传播时间； 为流量比重因素。X该公式可以重新整理为如下形式：Vstored=KXqin+(1-X)qout流量比重因素的下限为 0.0，上限为 0.5。这个因子是楔蓄量的函数。对于水库式蓄水，没有楔蓄， ；而对于一个完整的楔蓄， ；对于河X=0.0 X=0.5流， 落在 0.0 和 0.3 之间，其平均值接近 0.2。X对于

28、蓄水容量的定义可以加入连续公式并简化为：qout， 2=C1qin， 2+C2qin， 1+C3qout， 1式中： 为该时间段开始时的入流量，m 3/s； 为该时间段结束时的入qin， 1 qin， 2流量，m 3/s； 为该时间段结束时的出流量，m 3/s； 为该时间段结束qout， 1 qout， 2时的出流量，m 3/s。C1= t-2KX2K(1-X)+tC2= t+2KX2K(1-X)+tC3=2K(1-X)-t2K(1-X)+t其中， 。为了用体积单位表示所有值，在蓄水容量公式C1+C2+C3=1两端乘以该时间段得到：Vout， 2=C1Vin， 2+C2Vin， 1+C3Vqo

29、ut， 1为了保持数值稳定和避免出现负出流量的计算，必须满足以下条件：2KXt2K(1-X)流量比重因素 X 的值由使用者输入，蓄水时间常数的值估算如下：K=coef1Kbnkfull+coef2K0.1bnkfull式中： 为稳定流情况下的河段传播时间，s； 和 为权重系数，由使K coef1 coef2用者输入； 为稳定刘情况下渠道蓄满水的河段传播时间，s；Kbnkfull为渠道蓄满 1/10 水量时河段传播时间， s。K0.1bnkfull要计算 和 ，Cunge 于 1969 年提出一个公式：KbnkfullK0.1bnkfullK=1000Lchck式中： 为渠道长度，km； 为指定

30、深度处的波速， m/s。Lch ck波速的计算公式为：ck= ddAch(qch)其中流速 由曼宁公式求解，将曼宁公式代入上式得：qchck=53(R23chslp12chn )=53vc式中： 为湿周，m； 为河段坡度，%； 为曼宁系数； 为流速，m/s。Rch slpch n vc（2）变储量演算法对于一个给定的河段，储量演算基于连续方程 ，可写Vin-Vout=Vstored为：t(qin， 1+qin， 22 )-t(qout， 1+qout， 22 )=Vstored， 2-Vstored， 1式中： 为稳定流情况下的河段传播时间，s ； 为该时间段开始时的入流t qin， 1量，m

31、 3/s； 为该时间段结束时的入流量，m 3/s； 为该时间段开始时qin， 2 qout， 1的出流量，m 3/s； 为该时间段结束时的出流量，m 3/s； 为该时间qout， 2 Vstored， 1段开始时的蓄水容量，m 3； 为该时间段结束时的蓄水容量，m 3。Vstored， 2演进时间是由渠道中的水容量除以水流流量：TT=Vstoredqout=Vstored， 1qout， 1 =Vstored， 2qout， 2式中： 为蓄水容量，m 3； 为出流量，m 3/s。Vstored qout联合以上两式，简化为：qout， 2=SCqin， ave+(1-SC)qout， 1其中，

32、 为蓄水系数， ； 为入流平均流量， SC SC=2t2TT+t qin， ave，m 3/s。qin， ave=qin， 1+qin， 22整理上式得：qout， 2=SC(qin， ave+Vstored， 1t )为用体积单位来表达所有变量，公式两边都乘以时间段，则有：Vout， 2=SC(Vin+Vstored， 1)3 软件操作3.1 输入数据要创建 SWAT 数据集，ArcSWAT 需要访问提供流域信息的准确类型的ArcGIS 的栅格、矢量和数据文件。这些必需的地图主题和数据文件必须在运行SWAT 之前准备好。3.1.1 空间数据空间数据包括必需的 ARCSWAT 空间数据和可选的

33、 ARCSWAT 空间数据。SWAT 空间数据集可以以任何投影类型来创建，但所有的地图必须用同一投影。（1）必需的 ARCSWAT 空间数据 DEM，ESRI GRID Format。DEM 的高程值可以用整型或者实数型。决定地图分辨率的单位不必与高程的单位保持一致。例如地图分辨率可以是米（m） ，而高程可以是英尺（feet） 。地图分辨率单位必须定义为以下几种：米（m） ，公里（km） ，英尺（feet） ，码（yards） ，英里（miles） 。高程单位必须定义为以下几种：米（m） ，厘米（cm） ，码（yards） ，英尺（feet ） ，英寸（ inches） 。 Land Cove

34、r/Land Use，ESRI GRID，Shapefile ，or Feature Class Format。土地利用图中的类别需要重新分类为 SWAT 需要的土地利用类型。用户可以用三种方法重新分类土地利用，第一种是创建地图时用美国地质勘探局的分类代码；第二种是加载土地利用图时为每一分类选定SWAT 土地利用类型；第三种是为土地利用图的不同分类创建一个能够识别 4 位 SWAT 代码的 look up 表。 Soil，ESRI GRID，Shapefile ，or Feature Class Format。用户需要输入SWAT 土壤文件，在创建工程之前把图中每一类的土壤数据输入 User

35、Soil database 中。要对土壤图的土壤类型进行重分类，信息必须在界面中手动输入，列出土壤信息 look up 表（2）可选的 ARCSWAT 空间数据集 DEM Mask，ESRI GRID，Shapefile or Feature Class Format。模型可以加载一个 mask 在 DEM 之上。 STREAMS，Shapefile or Feature Class Format。在那些地势很低的地方，DEM 地图网格不能精确地判断河流的位置，就需要有河网描绘的shapefile 线性文件。 User-Defined Watersheds，Shapefile or Featu

36、re Class Format。用户可以自定义子流域划分加载到模型中去，如果选择自定义子流域也必须自定义河网。 User-Defined Streams，Shapefile or Feature Class Format。用户自定义子流域后，模型将不会进行河网提取，需要用户自行定义。3.1.2 属性数据属性数据包括：土壤属性数据、气象资料数据和水文数据。SWAT 模型通过三个个数据库来存储属性数据，其中气象资料数据通过模型自带的数据库进行存储，另外自定义的土壤属性数据和气象站参数数据通过两个附加的数据库进行存储，水文数据用于模拟结果率定。（1）土壤属性数据模型需要的土壤属性数据包括两大类：土壤

37、物理属性和土壤化学属性。土壤的物理属性对土壤剖面中水和气的运动状况起决定作用，并显著影响水文响应单元（HRU）中水文过程的模拟。土壤物理属性参数主要包括土层厚度、土壤层数、密度、土壤颗粒大小分布和土壤饱和水力传导率等。由于获得的土壤属性数据难以完全满足构建 SWAT 模型土壤库的要求，因此，部分无法直接获取的参数如土壤可利用有效水、饱和水力传到率等采用美国农业部开发的土壤水特性程序 SPAW 进行估算；部分参数通过对前人的研究成果和经验公式获得。估算结果仅简略反映土壤参数特征，在模型参数校准的过程中还需做进一步调整。（2）气象、水文资料数据SWAT 模型所需要的气象数据包括气象站点的位置以及各

38、气象站点的实测数据。气象站点包括雨量站、湿度站、温度站、风速站、太阳辐射站，其实测数据分别为日降水量、相对湿度、日最高/最低气温、风速和太阳辐射。本文采用北京站的气象资料来构建模型“天气发生器” 所需的数据。 SWAT 模型自带的天气发生器 WXGEN，其功能有二，一是用于生成气候数据，通过大量基础气象数据构建完成之后，可以用来生成任意年份的气象数据，二是对缺失的数据进行填补。3.2 输出数据SWAT 模型每次模拟都会生成五个输出文件：输入汇总文件（input.std ） 、输出汇总文件（output.std） 、HRU 输出文件（output.hru） 、子流域输出文件（output.sub

39、） 、河道输出文件（output.rch ） 。3.2.1 output.hruHRU 输出文件包含了流域中每一个水文响应单元的信息，文件中输出数据介绍如下：MON：时间步长，根据用户选择可以为日、月、年；AREA：HRU 的面积，km 2；PRECIP：时间步长内 HRU 的降雨量，mm；SNOFALL：时间步长内 HRU 的降雪量，mm ；SNOMELT：时间步长内 HRU 的冰雪融化量，mm；IRR：灌溉用水量，mm ；PET：潜在蒸散发量，mm；ET：实际蒸散发量，mm；SW-INIT：土壤初始含水量，mm；SW-END：时间步长后土壤含水量，mm；PERC：植物蒸腾水量，mm；GW-

40、RCHG：入渗量，mm；DA-RCHG：深层地下水补给量，mm；REVAP：浅层地下水向表层土壤回流量，mm；SA-IRR：浅层地下水灌溉量，mm；DA-IRR：深层地下水灌溉量，mm ；SA-ST：浅层地下水储水量，mm；DA-ST：深层地下水储水量，mm ；SURQ-GEN：地表径流量， mm；SURQ-CNT：地表径流进入河道水量，mm ；TLOSS：输移损失水量，mm；LATQ：侧向流水量，mm；GW-Q：地下水补充河道水量，mm；WYLD：总产流量，包括地表径流、地下水、侧向流，扣除损失水量，mm；DAILYCN：日径流曲线值；TMP-AV：日平均温度，摄氏度；TMP-MX：日最高温

41、度，摄氏度；TMP-MN：日最低温度，摄氏度；SOL-TMP： 土壤温度，摄氏度；SOLAR：日辐射量，MJ/m 2；SYLD：产沙量，t/ha；USLE：土壤流失量， t/ha；N-APP：氮肥施用量，kg/ha ；P-APP：磷肥施用量，kg/ha；NRAIN：氮素随降雨进入土壤量，kg/ha；NFIX：植物固氮量，kg/ha ；ORGN：有机氮产生量，kg/ha；ORGP：有机磷产生量， kg/ha；SEDP：随泥沙磷流失量，kg/ha；NSURQ：地表径流中硝态氮量，kg/ha；NLATQ：侧向流中硝态氮含量，kg/ha；NO3L：土壤中被冲刷出的硝态氮，kg/ha；NO3GW：地下水

42、流失硝态氮量，kg/ha；SOLP：溶解态磷产生量，kg/ha；P-GW：地下水溶解态磷流失量，kg/ha。3.2.2 output.subMON：时间步长，据用户选择可以为日、月、年；AREA：子流域面积，km 2；PRECIP：降水量，mm ；SNOMELT：冰雪融化量，mm ；PET：潜在蒸发量，mm；ET：实际蒸发量，mm；SW：土壤含水量， mm；PERC：入渗量，mm ；SURQ：地表径流进入河道水量，mm ；GW-Q：地下水补给河道水量，mm；WYLD：流域产水量， mm；SYLD：流域产沙量，t/ha；ORGN：流域有机氮产生量， kg/ha；ORGP：流域有机磷产生量，kg/ha；NSURQ：地表径流中硝态氮含量，kg/ha；SOLP：流域溶解态磷产生量，kg/ha；SEDP：流域矿化磷产生量，kg/ha；3.2.3 output.rchMON：时间步长，据用户选择可以为日、月、年；AREA：河道集水面积，km 2；FLOW-IN：日平均入流量，m 3/s；FLOW-OUT：日平均出流量， m3/s；EVAP：河道蒸发量， m3/s；TLOSS：河道输移损失水量，m 3/s；SED-IN：泥沙入流量， t；SED-OUT：泥沙出流量，tSEDCONC：泥沙淤积量，mg/L；ORGN-IN：有机氮输入量， kg；ORGN-OUT：有机