SWAT水文模型.doc
《SWAT水文模型.doc》由会员分享,可在线阅读,更多相关《SWAT水文模型.doc(32页珍藏版)》请在咨信网上搜索。
1、(完整版)SWAT水文模型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产汇流模型模型中采用的水量平衡表达式为:式中:为土壤最终含水量,m
5、m;为土壤前期含水量,mm;t为时间步长,d;为第i天降水量,mm;为第i天的地表径流,mm;为第i天的蒸发量,mm;为第i天存在于土壤剖面底层的渗透量和侧流量,mm;为第i天地下水含量,mm。SWAT模型水文循环陆地阶段主要有水文、天气、沉积、土壤温度、作物产量、营养物质和农业管理等部分组成.模型径流量产生计算流程图如下图2-1所示。图21 SWAT模型产流计算流程图2。1.1地表径流当落到地表的降水量多余入渗量时产生地表径流。SWAT模型采用SCS径流曲线法计算。SCS曲线方程自上世纪50年代逐渐得到广泛使用,属于经验模型,是对全美小流域降水与径流关系20多年的研究成果。模型能反应不同土壤
6、类型和土地利用方式及前期土壤含水量对降雨径流的影响,它是基于流域的实际入渗量(F)与实际径流量(Q)之比等于流域该场降雨前的最大可能入渗量(S)与最大可能径流量()之比的假定基础上建立的SCS模型的降雨径流基本关系表达式如下:式中:假定潜在径流量()为降水量(P)与由径流产生前植物截留、初渗和填洼蓄水构成的流域初损()的差值.由此推到上式有:初损受土地利用、耕作方式、灌溉条件、冠层截留、下渗、填洼等因素的影响,它与土壤最大可能入渗量S呈一定的正比关系,美国农业部土壤保持局在分析了大量长期的实验结果基础上,提出了二者最合适的比例系数为0.2,即:由此可得SCS方程为: 流域当时最大可能滞留量在空
7、间上与土地利用方式、土壤类型和坡度等下垫面因素密切相关,模型引入的值可较好地确定,公式如下:CN是一个无量纲参数,CN值是反映降雨前期流域特征的一个综合参数,它是前期土壤湿度、坡度、土地利用方式和土壤类型状况等因素的综合.2.1。2蒸散发模型考虑的蒸散发是指所有地表水转化为水蒸气的过程,包括树冠截留的水分蒸发、蒸腾和升华及土壤水的蒸发。蒸散发是水分转移出流域的主要途径,在许多江河流域,蒸发量都大于径流量.准确地评价蒸散发量是估算水资源量的关键,也是研究气候和土地覆盖变化对河川径流影响的关键。(1)潜在蒸散发模型提供了PenmanMonteith、Priestley-Taylor和Hargrea
8、ves三种计算潜在蒸散发的方法,另外还可以使用实测资料或已经计算好的逐日潜在蒸散发资料.一般采用Penman-Monteith方法来计算流域的潜在蒸散发。(2)实际蒸散发实际蒸散发以潜在蒸散发为计算基础。在计算流域实际蒸散发量的时候,模型首先计算植物冠层截留水分的蒸发,然后计算最大蒸腾量、最大升华量和最大土壤蒸发量,最后计算实际的升华量和土壤水分蒸发量.(3)冠层截留蒸发量模型在计算实际蒸发时假定尽可能蒸发冠层截留的水分,如果潜在蒸发量小于冠层截留的自由水量,则:式中:为某日流域的实际蒸发量,mm;为某日冠层自由水蒸发量,mm;为某日的潜在蒸发量,mm;为某日植被冠层自由水初始含量,mm;为某
9、日植被冠层自由水终止含量,mm。如果潜在蒸发量大于冠层截留的自由水含量则:当植被冠层截留的自由水被全部蒸发掉,继续蒸发所需的水分就会从植被和土壤中得到。(4)植物蒸腾假设植物生长在一个理想的条件下,植物蒸腾可用以下表达式计算:当0LAI3.0时当LAI3.0时式中:为某日最大蒸腾量,mm;为植被冠层自由水蒸发调整后的潜在蒸发,mm;LAI为叶面积指数。因为没有考虑到植物下面图层的含水量问题,由此公式计算处的蒸腾量可能比实际蒸腾量要大一些。(5)土壤水分蒸发在计算土壤水分蒸发时,首先区分出不同深度土壤层所需要的蒸发量,土壤深度层次的划分决定土壤允许的最大蒸发量,可由下式计算:式中:为z深度处蒸发
10、需要的水量,mm;z为地表以下土壤的深度,mm。表达式中的系数是为了满足50%的蒸发所需水分来自土壤表层10mm,以及95的蒸发所需水分来自0100mm土壤深度范围内.土壤水分蒸发所需要的水量是有土壤上层蒸发需水量与土壤下层蒸发需水量决定的:式中:为ly层的蒸发需水量,mm;为土壤下层的蒸发需水量,mm;为土壤上层的蒸发需水量,mm。土壤深度的划分假设50的蒸发需水量由010mm内土壤上层的含水量提供,因此100mm的蒸发需水量中50mm都要由10mm的上层土壤提供,显然上层无法满足需要,这就需要建立一个系数来调整土壤层深度的划分,以满足蒸发需水量,调整后的公式可以表示为:式中:esco为土壤
11、蒸发调节系数,该系数是SWAT为调整土壤因毛细作用和土壤裂隙等因素对不同土层蒸发量二提出的,对于不同的esco值对应着相应的土壤层划分深度。2。1.3土壤水渗入到土壤中的水有多种不同运动方式。土壤水可以被植物吸收或蒸腾而损耗,可以渗透到土壤底层最终补给地下水,也可以在地表形成径流,即壤中流.由于主要考虑径流量的多少,因此对壤中流的计算简要概括.模型采用动力储水方法计算壤中流。相对饱和区厚度计算公式为:式中:为土壤饱和区内可流出的水量,mm;为山坡坡长,m;为土壤可出流的孔隙率;表示土壤层总孔隙度,即与土壤层水分含量达到田间持水量的孔隙度之差。山坡出口断面的净水量为:式中:为出口断面处的流速,m
12、m/h。其表达式为:式中:为土壤饱和导水率,mm/h;为坡度。总结上面表达式,模型中壤中流最终计算公式为:2.1。4地下水模型采用以下表达式来计算流域地下水:式中:为第i天进入河道的地下水补给量,mm;为第(i1)天进入河道的地下水补给量,mm;为时间步长,d;为第i天蓄水层的补给流量,mm;为基流的退水系数。其中补给流量由下式计算:式中:为第i天蓄水层补给量,mm;为补给滞后时间,d;为第i天通过土壤剖面底部进入地下含水层的水分通量,mm/d;为第(i1)天蓄水层补给量,mm。2。2土壤侵蚀模型泥沙生成量用MUSLE方程来预测,计算渠道泥沙输移量的公式为:式中:为输移能力,t/m3;为流速,
13、m/s;a和b是常数.根据天气条件,泥沙输移量可以高于或者低于输移能力,导致沉积过量的泥沙通过渠道侵蚀再悬浮输移泥沙。流速方程为:式中:为流量,m3/s;为渠道宽度,m;为径流深,m.对于低于齐岸深度的径流,径流深使用Manning方程来计算,假定渠道宽度远大于深度:式中:为渠道曼宁系数;为渠道坡度,m/m。由于降水和径流产生的土壤侵蚀是用MUSLE方程来计算的,MUSLE是修正的通用土壤流失方程(USLE).USLE方程是通过降水动能函数预测年均侵蚀量,而在MUSLE中,用径流因子代替降水动能,改善了泥沙产量的预测,这样就不需要泥沙输移系数,并且可以将方程用于单次暴雨事件,因为径流因子是先行
14、湿度和降水动能的函数。USLE中需要输移系数是因为降水动能因子表示的能量只在作用流域内起作用。修正的通用土壤流失方程为;式中:为土壤侵蚀量,t;为地表径流,mm/h;为洪峰流量,m3/s;为水文响应单元的面积,hm2;为土壤侵蚀因子;为植被覆盖和管理因子;为水土保持措施因子;为地形因子;为粗碎屑因子。2。2。1土壤侵蚀因子当其他影响侵蚀的因子不变时,因子反映不同类型土壤抵抗侵蚀力的高低。它与土壤物理性质的影响,如机械组成、有机质含量、土壤结构、土壤渗透性等有关.当土壤颗粒粗、渗透性大时,值就低,反之则高;一般情况下值得变幅在0.02 0.75之间。值得直接测定方法是:在标准小区(坡长为22.1
15、m,宽为1.83m,坡度为9)上没有任何植被,完全休闲,无水土保持措施,降水后收集由于坡面径流而冲蚀到集流槽内的土壤,烘干、称重,由公式计算得到值。试验测算值既费时又费力,1971年一些学者发展了一个通用方程来计算土壤侵蚀因子值,该方程在土壤黏土和壤土组成少于70%时适用。式中:为颗粒尺度参数;为有机物含量百分比,%;为土壤分类中的结构代码;为土壤剖面可渗透性类别。1995年另一位学者提出了一个替换方程:式中:为粗糙砂土质地土壤侵蚀因子;为粘壤土土壤侵蚀因子;为土壤有机质因子;为高砂质土壤侵蚀因子。各因子的计算公式如下:式中:为粒径在0。052。00mm沙粒的百分含量;为粒径在0.0020。0
16、5mm的淤泥、细砂百分含量;为粒径小于0.002mm的粘土百分含量;为各土壤层中有机碳含量,%。2。2。2植被覆盖因子植被覆盖和管理因子表示植物覆盖和作物栽培措施对防止土壤侵蚀的综合效益,其含义是在地形、土壤、降水条件相同的情况下,种植作物或林草地的土地与连续休养地土壤流失量的比值,最大取值为1。0。由于植被覆盖受植物生长期的影响,SWAT模型通过下面的方程调整植被覆盖和管理因子:式中:是最小植被覆盖和管理因子值;是地表植物残留量,kg/hm2。最小C因子可以由已知年平均C值,通过以下方程计算.式中:表示不同植被覆盖的年均C值。2。2.3水土保持措施因子水土保持因子是指有保持措施的地表土壤流失
17、与不采取任何措施的地表土壤流失的比值,这里的水土保持措施包括等高耕作、带状种植和梯田.等高耕作对于中低强度的降水侵蚀具有保护水土流失的作用,但对于高强度的降水其保护作用则很小,等高耕作对坡度为38%之间的土地非常有效。2。2。4地形因子地形因子的计算公式如下:式中:为坡长;为坡长指数;为坡度。坡长指数的计算公式如下:式中:为水文相应单元的坡度,.2.2.5因子因子是通过下面公式计算的:式中:为第一层土壤中砾石的百分比,%.2。3污染负荷模型SWAT模型可以模拟不同形态氮的迁移转化过程,地表径流流失、入渗淋失、化肥输入等物理过程,有机氮矿化、反硝化等化学过程以及作物吸收等生物过程,氮可以分为有机
18、氮、作物氮和硝酸盐氮三种化学状态,氮的生物固定、有机氮向无机氮的转化以及溶解性氮随侧向壤中流的迁移等过程,有机氮又被划分为活泼有机氮和惰性有机氮两种状态,以及铵态氮挥发过程的模拟。2.3.1硝态氮污染负荷模型硝态氮主要随地表径流、侧向流或渗流在水体中迁移,要计算随水体迁移的硝态氮量必须先计算自由水中硝态氮浓度,用这个浓度乘以各个水路流动水的总量,即可得到从土壤中流失的硝态氮总量。自由水部分的硝态氮浓度可用下面公式计算:式中:为自由水中硝态氮浓度,kg/mm;为土壤中硝态氮的量,kg/hm2;为土壤中自由水的量,mm;为孔隙度;为土壤饱和含水量.(1)通过地表径流流失的溶解态氮计算公式:式中:为
19、通过地表径流流失的硝态氮,kg/hm2;为硝态氮渗流系数;为地表径流,mm。(2)通过侧向流流失的溶解态氮的量计算公式:对于地表10mm图层:对于10mm以下的土层:式中:为通过侧向流流失的硝态氮,kg/hm2;为硝态氮渗流系数;为自由水的硝态氮浓度,kg/mm;为侧向流,mm。(3)通过渗流流失的溶解态氮量计算公式:式中:为通过渗流流失的硝态氮,kg/hm2;为自由水的硝态氮浓度,kg/mm;为渗流,mm。2。3.2有机氮污染负荷模型有机氮通常是吸附在土壤颗粒上随径流迁移的,这种形式的氮负荷与土壤流失量密切相关,土壤流失量直接反映了有机氮负荷。有机氮随土壤流失的输移量计算公式为:式中:为有机
20、氮流失量,kg/hm2;为有机氮在表层土壤中的浓度,kg/t;为土壤流失量,t;为水文相应单元的面积,hm2;为氮富集系数,氮富集系数是随土壤流失的有机氮浓度和土壤表层有机氮浓度的比值。计算富集系数的公式如下:式中:为地表径流中泥沙含量。的计算公式如下:式中:为土壤流失量,t;为水文相应单元面积,hm2;为地表径流,mm。2.3。3溶解态磷污染负荷模型溶解态磷在土壤中的迁移主要是通过扩散作用实现的,扩散是指离子在微小尺度下由于浓度梯度而引起的溶质迁移,由于溶解态磷不很活跃,所以由地表径流以溶解态形式带走的土壤表层的磷很少,地表径流输移的溶解态磷可由下面公式计算:式中:为通过地表径流流失的溶解态
21、磷,kg/hm2;为土壤中溶解态磷,kg/hm2;为土壤溶质密度,mg/m3;为表层土壤深度,mm;为土壤磷分配系数,表层土壤中溶解态磷的浓度和地表径流中溶解态磷浓度的比值。2.3.4有机磷和矿物质磷污染负荷模型有机磷和矿物质磷通常是吸附在土壤颗粒上通过径流迁移的,这种形式的磷负荷与土壤流失量密切相关,土壤流失量直接反映了有机磷和矿物质磷负荷,有机磷和矿物质磷随土壤流失输移量计算公式为:式中:为有机磷流失量,kg/hm2;为有机磷在表层土壤中的浓度,kg/t;为土壤流失量,t;为水文相应单元的面积,hm2;为磷富集系数。2。4河道演算模型模型中定义的河道均是明渠流,SWAT用曼宁公式来定义河道
22、糙率和水流流速.水流在河道中演进过程使用变储量演算法或马斯金根法,两种方法都是动力波方程.在模拟中一般使用马斯金根法来计算.(1)马斯金根法马斯金根法假设河道内水体形状是由一个楔形蓄水体和一个棱形蓄水体组成,如下图22所示。图2-2 河道水体示意图当洪水波行进到某个河段槽,入流量大于出流量便形成了楔形蓄水体。当洪水波退去,在河段槽便出现了出流量大于入流量的负楔蓄体.另外对于楔蓄水体,河段槽内始终包含一个体积为流域长度上横截面不变的棱柱状水体。总的蓄水容量为:式中:为蓄水容量,m3;为入流量,m3/s;为出流量,m3/s;为稳定流情况下的河段传播时间;为流量比重因素。该公式可以重新整理为如下形式
23、:流量比重因素的下限为0。0,上限为0.5。这个因子是楔蓄量的函数。对于水库式蓄水,没有楔蓄,;而对于一个完整的楔蓄,;对于河流,落在0.0和0。3之间,其平均值接近0。2。对于蓄水容量的定义可以加入连续公式并简化为:式中:为该时间段开始时的入流量,m3/s;为该时间段结束时的入流量,m3/s;为该时间段结束时的出流量,m3/s;为该时间段结束时的出流量,m3/s。其中,.为了用体积单位表示所有值,在蓄水容量公式两端乘以该时间段得到:为了保持数值稳定和避免出现负出流量的计算,必须满足以下条件:流量比重因素X的值由使用者输入,蓄水时间常数的值估算如下:式中:为稳定流情况下的河段传播时间,s;和为
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- SWAT 水文 模型
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【天****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【天****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。