第三讲 SWAT水文循环原理

第三讲SWAT水文循环原理

SWAT水文循环提纲1.1SWAT地表水过程1.2SWAT蒸散发过程1.3SWAT土壤水过程1.4SWAT地下水过程1.1SWAT地表水过程水量平衡方程如下：1.1SWAT地表水过程产流计算：产流计算：CN值可针对不同的土壤类型、土地利用和植被覆盖的组合查表获得，CN值是无量纲的反映降雨前期流域特征的一个综合参数，将前期土壤湿度、坡度、土地利用方式和土壤类型状况等因素综合在一起。1.1SWAT地表水过程产流计算：1.1SWAT地表水过程SWAT水库模拟水库的水量平衡方程为：SWAT水库模拟水库出流计算（Reservoiroutflow

）:

实测日出流数据实测月出流数据

没有控制的水库年平均出流

控制水库的目标出流数据SWAT水库模拟SWAT水库模拟水库调控计算（Targetreleaseforcontrolledreservoir）:提纲1.1SWAT地表水过程1.2SWAT蒸散发过程1.3SWAT土壤水过程1.4SWAT地下水过程1.2SWAT蒸散发过程蒸散发包括冠层截留水蒸发、蒸腾和升华及土壤水的蒸发。蒸散发是水分转移出流域的主要途径。准确地评价蒸散发量是估算水资源量的关键，也是研究气候和土地覆被变化对河川径流影响的关键问题。1）潜在蒸散发SWAT提供了Penman-Monteith、Priestley-Taylor和Hargreaves等三种计算潜在蒸散发能力的方法，另外还可以使用实测资料或已经计算好的逐日潜在蒸散发资料。蒸发模拟：1）潜在蒸散发1.2SWAT蒸散发过程1）潜在蒸散发1.2SWAT蒸散发过程1）潜在蒸散发1.2SWAT蒸散发过程2）实际蒸散发1.2SWAT蒸散发过程2）实际蒸散发1.2SWAT蒸散发过程2）实际蒸散发1.2SWAT蒸散发过程2）实际蒸散发1.2SWAT蒸散发过程提纲1.1SWAT地表水过程1.2SWAT蒸散发过程1.3SWAT土壤水过程1.4SWAT地下水过程渗漏percolation浅层地下水Shallowaquifer深层地下水Deepaquifer壤中流lateral

flow基流returnflow入渗infiltrationloss蒸散发Evapotranspiration1.3SWAT土壤水过程1.3SWAT土壤水过程1.3SWAT土壤水过程提纲1.1SWAT地表水过程1.2SWAT蒸散发过程1.3SWAT土壤水过程1.4SWAT地下水过程1.4SWAT地下水过程1.4SWAT地下水过程河流与地下水的交互关系取用水Wateruse入渗补给Recharge深层地下水补给Todeepaquifer基流Returnflow潜水蒸发Revap浅层含水层ShallowAquifer浅层地下水模拟=percolationcoefficient深层地下水模拟Wdeep(mm)onlyifamountofwaterinshallowaquiferexceedstresholdvalue取用水Wateruse渗漏Percolation深层含水层DeepAquifer浅层地下水平衡方程：1.4SWAT地下水过程地下水位模拟（

Groundwaterheight）：深层地下水平衡方程：1.4SWAT地下水过程谢谢！Thanks!附件：SWAT作物生长模拟与农业管理作物生长模块的结构作物生长模块概述有效积温的计算生物量的同化对水分的吸收对N、P的吸收实际生长耕作措施作物生长模块概述

SWAT模型中作物生长模块是一种适合多种作物生长的通用模型，能够模拟上百种大田作物、园艺作物、草原牧草和树木生长，其特点是根据各种作物生理生态过程的共性研制成模型的主体框架，再结合各种作物的生长和田间管理参数分别进行各种作物生长与产量模拟。模型对作物生长与产量形成过程，及其对环境响应进行定量模拟的基本过程如下：作物模型在逐日气候要素的驱动下，首先模拟计算太阳辐射能转化为干物质的数量，涉及到估算光截获数量的叶面积变化动态，主要通过最大叶面积系数、叶面积变化的S形曲线形态参数、叶面积下降速度等作物生长参数来确定。同时，光合产物转化受作物水汽压差和大气CO2浓度的影响。叶面积的增长和变化与温度、水分、养分等作物生长因素密切相关。通过作物生长的最适温度、最低温度与作物生育过程计算湿度对叶面积增长的影响效应。通过计算根系分布层次的土壤水分和养分状况，估计水分和养分胁迫对叶面积增长的影响效应。根据作物生长期生物量增长的、氮、磷和温度等四种胁迫因子，校正干物质的合成数量。计算分配到地上部和根系的干物质数量，估算根系分布深度和根系密度，进而计算作物根系吸收的土壤水分和养分在不同土层的分布情况。而植株水分、氮素、磷素满足程度最终又影响干物质合成数量。通过作物水分和养分临界期水分、养分亏缺状况估计实际收获指数，最后通过地上部生物量和收获指数计算可供收获的作物经济产量。有效积温为日平均温度减去基础温度计算公式如下：在全生育期内将HU中大于零的部分累加，就得到植物生育期有效积温的总和。农事活动可以依据达到有效积温的某个比例进行安排。有效积温的计算生物量的同化植物叶面积截获光能的量，通过Beer法则来计算:

利用截获的光，根据光能利用率，就得到时间步长内生物量的最大增量，公式如下：

SWAT考虑了CO2浓度对光能利用率的影响。CO2浓度对光能利用率校正的公式如下：冠层高度的计算公式如下：其中，hc是冠层高度，hc,max是冠层最大的高度，frLAIm是植物最大LAI的分数。LAI的计算公式如下：

LAI达到最大值后，这个最大值将会维持一段时间，直到叶片老化的速率超过叶片生长补充的速率

根系的发育：分配到根系的干物质从出苗期的40%变化到成熟期的20%。根系逐日的分配系数依下式计算：其中，frroot是该日分配给根系的比例，frPHU是在达到生长季指定的某天时积累的潜在热量单位的分数。根深的计算，依作物不同而不同，对一年生植物，开始时根系的深度被设置为10.0mm，随后发育阶段为0.4时达到最大，维持在最大值。在发育阶段小于0.4时，根深依下式计算：对水分的吸收在根区任何深度的土层下，利用下式来计算植物的潜在吸水量：其中，wup,z是某一特定深度的土层内吸收的水量，Et是该日的最大蒸腾量，是水分利用分布参数，z是土层深度，zroot是该日根深。通过代入土层顶层和底层，取其差值后再对公式求解，可以得到任意土层的潜在吸水量。随着土壤含水量的降低，土壤水分被土壤颗粒吸附得越来越牢，植物的吸水因此也变得越来越困难。为了反应这个现象，利用下面这个公式来修正植物的潜在吸收水量：whenwhen吸收的水量经过校正后，就可以用来下式来计算从某一土层吸收的水量：作物对N、P的吸收植物对氮素的吸收受植物氮素平衡的控制。植物氮素平衡公式计算植物生物量中的含氮量，在最适生长条件下，这个含氮量是植物生长阶段的函数。计算植物N需求的公式如下：被吸收的氮素在不同深度土层中的分配，按下式计算：地表根系的密度最大，因此地表层氮素吸收的强度也最大。这种随深度不同而异的分布由氮素吸收分布参数来控制，这是一个用户可以自己调整的参数。如果根区的含氮量不足以满足豆类植物的需要，SWAT允许植物通过固氮获得额外的氮素。固氮量是土壤水分，土壤氮素和植物发育阶段的函数：fgr是固氮能力受生育期的影响最大，需要另外进行校正。植物P吸收的计算与N吸收的计算类似，受植物P平衡的控制，在最适生长条件下，植物对P的吸收受植物中含P量的影响，并且是发育阶段的函数：计算植物P需求的公式如下：吸收的P在不同深度的土层内的分布依下式计算：作物实际生长上述计算的是潜在条件下的生长，需要进行水分、温度、N和P胁迫的校正后就是实际条件下的生长。水分胁迫的计算公式如下：其中，wstrs是该日的水分胁迫影响因子，Et是该日的最大蒸腾量，Et,act是该日的实际蒸腾量，wactualup是该日植物总的吸水量。温度胁迫是日平均温度和植物生长最适温度的函数，在最适温度附近，植物不会受到温度的胁迫。然而，如果偏离较大，植物就会开始受到温度的胁迫。计算公式如下：whenwhenwhenwhenN胁迫的计算仅限于对非豆类作物，SWAT模型不允许豆类作物受到N胁迫。N胁迫的计算公式如下：nstrs是氮素胁迫程因子，依下式计算：

其中，bioN,opt是当前发育阶段下植物组织里的最适含氮量，bioN是实际贮存在植物干物质中的氮素。

P素胁迫程度的计算，可通过比较植物实际的与最适的P素含量得到。当植物的实际含P量等于最适含P量时，P素的胁迫等于0.0，当实际含P量低于最适含P量的50%时，P素的胁迫等于1.0，这个P素胁迫随含P量的下降而增加的过程也是非线性的，依下式计算：pstrs是P素胁迫程因子，依下式计算：植被的实际生产，取水分、温度、N和P胁迫中最大的胁迫因子，即依下式计算：计算出的潜在生长时的干物质就被乘以这个值，计算公式如下：植物在该日实际增加的叶面积依下式计算：

对于植物每年生产的生物量的多少，模型允许用户自己定义一个值。当这个值在模型的管理文件(.mgt)文件中被设立后，模型中每年的上述胁迫因素就不再考虑。当这个值被设立后，生物量的计算方法就如下式所示：模型允许用户指定一个固定的收获指数，当这个参数被指定时，这个参数值就取代了公式计算出来的收获指数，此时不考虑生育期的不同，也不计算水分亏缺对其的影响。其中，HIact是实际的收获指数，HI