《SWAT模型培训综合》由会员分享，可在线阅读，更多相关《SWAT模型培训综合（238页珍藏版）》请在金锄头文库上搜索。
1、SWATSWAT模型培训讲义模型培训讲义单位：中国水科院水资源所单位：中国水科院水资源所 讲解人：讲解人： 陆垂裕陆垂裕日期：日期：2012-2-222012-2-22自我介绍1976年10月 陆垂裕，出生于湖北省黄石市19951999 武汉水利电力大学，本科， 农田水利19982001 武汉大学，硕士， 水利水电工程20012004 武汉大学，博士， 水利水电工程20042006 中国水科院，博士后，水文水资源2006至今 中国水科院，水资源所，高工研究领域： 土壤水/地下水，地表/地下水质， 行业软件开发，水文模型等讲 义 提 纲水文科学发展简述SWAT模型简介SWAT主要水循环模拟方法详
2、解SWAT的主要物质循环模拟简介SWAT的输入输出SWAT的软件界面SWAT的参数敏感性分析、自动率定等SWAT应用研究示例讲解基于SWAT原理改进的水循环模拟模型MODCYCLE自我介绍联系方式：中国水利水电科学研究院水资源所北京市复兴路甲1号A座920室QQ：63711279手机：13601081649水文水文/ /水循环模型研水循环模型研究简述究简述水文科学发展简述水文科学发展简述 国际上，特别是西方认为以国际上，特别是西方认为以1717世纪进世纪进行的有史以来首批重要水文量测为标志，行的有史以来首批重要水文量测为标志，作为科学水文学的开始。这些量测是：作为科学水文学的开始。这些量测是：
3、Pierre PerraultPierre Perrault的塞纳河流域降雨、蒸发的塞纳河流域降雨、蒸发和毛管势量测；和毛管势量测；Edme MariotteEdme Mariotte的塞纳河流量量测；的塞纳河流量量测；Edmond HalleyEdmond Halley的地中海蒸发和流量量测等的地中海蒸发和流量量测等等。等。 1818、1919世纪主要在水力学世纪主要在水力学世纪主要在水力学世纪主要在水力学试验研究试验研究试验研究试验研究方面发展方面发展方面发展方面发展18561856年年 Henry DarcyHenry Darcy通过渗透试验得出水在多孔介质中的渗透定律，通过渗透试验得出
4、水在多孔介质中的渗透定律，即著名的达西公式；即著名的达西公式；18631863年年J.DupuitJ.Dupuit以以DarcyDarcy公式为基础提出了公式为基础提出了DupuitDupuit假设，为开采井假设，为开采井的稳定流计算奠定了基础的稳定流计算奠定了基础18801880年雷诺根据试验提出了针对流体力学的雷诺公式；年雷诺根据试验提出了针对流体力学的雷诺公式；18891889年曼宁提出了明渠稳定流公式；年曼宁提出了明渠稳定流公式； 进入进入进入进入2020世纪初水文学科开始快速发展世纪初水文学科开始快速发展世纪初水文学科开始快速发展世纪初水文学科开始快速发展4040年代以前：水文要素实
5、验和理论研究年代以前：水文要素实验和理论研究入渗：入渗：入渗：入渗：19111911年年, Green-Ampt , Green-Ampt 建立了饱和入渗理论公式，至今应用建立了饱和入渗理论公式，至今应用仍然很广泛；仍然很广泛；蒸发：蒸发：蒸发：蒸发：19261926年年, Bowen, Bowen提出第一个理论公式，其中将蒸发通量的潜提出第一个理论公式，其中将蒸发通量的潜热和能量通量的感热通过一个比值联系起来，这就是著名的波温比。热和能量通量的感热通过一个比值联系起来，这就是著名的波温比。这是脱离经验公式的显著进步。这是脱离经验公式的显著进步。非饱和水分运动：非饱和水分运动：非饱和水分运动：
6、非饱和水分运动： 1907 1907 年，年， BurkinghamBurkingham引进毛管势概念，即土壤水基质势概念，引进毛管势概念，即土壤水基质势概念，建立建立非饱和水分运动的理论基础；非饱和水分运动的理论基础； 1931 1931 年，年，Richards Richards 建立控制建立控制非饱和土壤水的等温输送的基本方程，非饱和土壤水的等温输送的基本方程，提出一种量测土壤水基质势的方法，并且研究了土壤水运动的滞后提出一种量测土壤水基质势的方法，并且研究了土壤水运动的滞后现象。现象。饱和地下水分运动：饱和地下水分运动：饱和地下水分运动：饱和地下水分运动： 19351935年年C.V.
7、TheisC.V.Theis（泰斯）提出了地下水向承压水井运动的非稳定（泰斯）提出了地下水向承压水井运动的非稳定 流公式，开创了现代地下水运动理论的新纪元；流公式，开创了现代地下水运动理论的新纪元；河川径流：河川径流：河川径流：河川径流：1907 1907 年，年， Hoyt Hoyt 发表了流速仪测流的经典专著，总结了用发表了流速仪测流的经典专著，总结了用流速仪测流的实践经验。流速仪测流的实践经验。降雨产流：降雨产流：降雨产流：降雨产流： 19331933年，年，HortonHorton提出降雨产流的入渗理论，开始了对陆面通量的从提出降雨产流的入渗理论，开始了对陆面通量的从经验向推理的过渡；
8、经验向推理的过渡； 19321932年，年，ShermanSherman提出流域响应的单位线方法。提出流域响应的单位线方法。 两者结合形成完整的流域水文推理和计算体系，开创了流域水文的两者结合形成完整的流域水文推理和计算体系，开创了流域水文的定量途径。定量途径。4040404060606060年代：应用水文学全面发展年代：应用水文学全面发展年代：应用水文学全面发展年代：应用水文学全面发展河道洪水演算：河道洪水演算：河道洪水演算：河道洪水演算： 19391939年，年，McCarthyMcCarthy提出提出MuskingumMuskingum洪水演算方法洪水演算方法（水文学法）（水文学法）；
9、19411941年，年，MeyerMeyer提出提出LagLag和和routeroute方法方法（槽蓄方程，水文学法）（槽蓄方程，水文学法）； 19511951年，年，HayamiHayami提出扩散相似法提出扩散相似法（扩散波方程，水力学法）（扩散波方程，水力学法）流域汇流：流域汇流：流域汇流：流域汇流： 19451945年，年，ClarkClark提出在等流时线基础上，通过一个概念性蓄水单元提出在等流时线基础上，通过一个概念性蓄水单元调蓄，建立单位线的推理方法，开创了单位线的理论研究时期。调蓄，建立单位线的推理方法，开创了单位线的理论研究时期。 蒸发：蒸发：蒸发：蒸发： 蒸发研究更系统，更
10、具理论价值。蒸发研究更系统，更具理论价值。 19351935，19481948年，年，ThornthweiteThornthweite对蒸发估算的改进；对蒸发估算的改进； 19481948年，年，PenmanPenman提出潜在腾发概念及潜在腾发估算的混合公式；提出潜在腾发概念及潜在腾发估算的混合公式； 19481948年开始，年开始，BudykoBudyko发表关于蒸发的系列文章。发表关于蒸发的系列文章。6060606080808080年代：流域产流试验研究年代：流域产流试验研究年代：流域产流试验研究年代：流域产流试验研究19631963年，年，HewlittHewlitt和和Hilbert
11、Hilbert的排水期斜坡土壤水分和能量实验的排水期斜坡土壤水分和能量实验19691969年，年，BetsonBetson的暴雨径流产流面积实验的暴雨径流产流面积实验19701970年，年，DunneDunne和和BlackBlack的小流域透水土壤产流实验的小流域透水土壤产流实验 这些研究大大丰富了对流域产流过程的认识，提出了若干不同于这些研究大大丰富了对流域产流过程的认识，提出了若干不同于HortonHorton产流入渗理论的降雨产流机制，并突出认识到流域部分面积产产流入渗理论的降雨产流机制，并突出认识到流域部分面积产流概念和山坡非饱和水分侧向运动的重要作用流概念和山坡非饱和水分侧向运动的
12、重要作用. . 70707070年代以后：水文学进入模型计算时代年代以后：水文学进入模型计算时代年代以后：水文学进入模型计算时代年代以后：水文学进入模型计算时代地下水数值模型地下水数值模型19881988年，年，MODFLOWMODFLOW（美国）（美国）流域水文模型流域水文模型19771977年，新安江模型（中国）年，新安江模型（中国）19791979年，年，TOPMODELTOPMODEL（纽芬兰）（纽芬兰）19861986年，年，SHESHE模型（英国、法国、荷兰联合开发）模型（英国、法国、荷兰联合开发）本世纪前半叶：本世纪前半叶：水文测验水文测验水文测验水文测验积累水文资料积累水文资料
13、积累水文资料积累水文资料工程实践工程实践工程实践工程实践积累经验和问题积累经验和问题积累经验和问题积累经验和问题本世纪后半叶：本世纪后半叶：水文研究水文研究水文研究水文研究山坡水文实验山坡水文实验山坡水文实验山坡水文实验流域水文模拟流域水文模拟流域水文模拟流域水文模拟计算机出现计算机出现计算机出现计算机出现水循环的各个子系统大气环流系统大气环流系统地表水系统地表水系统土壤水系统土壤水系统地下水系统地下水系统气象学领域地表水文学领域农田水利学领域地下水文学领域大气环流系统地表水系统土壤水系统地下水系统各个学科领域对应的水循环子系统水循环研究的目的生产实践中与水循环有关的各种问题：生产实践中与水循
14、环有关的各种问题：天气预报天气预报降雨时间？降雨量？降雨范围？气候变化气候变化干旱化？湿润化？防洪调度防洪调度降雨径流关系？洪峰流量？洪峰移动速度？泥沙淤积泥沙淤积冲淤平衡流量？河道输沙量？水力发电水力发电上游来水量？来水过程？气象学地表水文学农田灌溉农田灌溉灌溉方式？灌溉制度？土壤盐渍化土壤盐渍化埋深控制？排水沟间距？大坝、基坑防渗大坝、基坑防渗/排水排水渗透压力？减压控制？地下水水源地开发地下水水源地开发开采井间距？开采总量控制？海水入侵海水入侵咸/淡界面位置？平衡关系?农田水利学地下水文学流域水文流域水文/水循环模型水循环模型优势：优势： 各种水文过程统一模拟的平台，使人们可以站在流域各
15、种水文过程统一模拟的平台，使人们可以站在流域整体视角研究各种水文过程的综合影响整体视角研究各种水文过程的综合影响 降雨/蒸发、产流、河网汇流、土壤水运动、地下水运动流域水文流域水文/水循环模型水循环模型中国中国10大一级流域分区大一级流域分区流域：由分水线所流域：由分水线所包围的河流集水区包围的河流集水区流域水文流域水文/水循环模型水循环模型淮河流域淮河流域淮南市淮河干流土壤水模拟土壤水模拟产流、入渗、渗漏、蒸发问题产流、入渗、渗漏、蒸发问题地下水模拟地下水模拟河道基流问题地表水模拟地表水模拟数字高程数字高程流域划分、汇流路径问题河道洪水演进问题地下水运动土壤水运动降水入渗蒸散水文过程地表径流
16、汇流流域出流流流域域层层面面上上的的统统一一水资源研究方面的需求水资源评价水资源开发利用评价重点：“四水转化”规律反演及水平衡分析水文研究方面的需求：防洪、径流预测重点：产汇流计算，洪水过程模拟 流域水文模型的应用层面流域水文模型的应用层面概念性概念性概念性概念性：斯坦福模型，斯坦福模型，斯坦福模型，斯坦福模型，1966196619661966 水箱模型，水箱模型，水箱模型，水箱模型，1967196719671967 新安江模型，新安江模型，新安江模型，新安江模型，1977197719771977 物理性物理性物理性物理性：FreezeFreezeFreezeFreeze模型，模型，模型，模型
17、，1971197119711971 SHE SHE SHE SHE模型，模型，模型，模型，1986198619861986 TOPMODEL TOPMODEL TOPMODEL TOPMODEL，1979197919791979 WEP WEP WEP WEP，2001200120012001 流域水文模型流域水文模型 流域水文模型流域水文模型水箱模型水箱模型水箱模型水箱模型(TANK)(TANK)顶层水箱：模拟地表径流，设置个出流口顶层水箱：模拟地表径流，设置个出流口第二层水箱：模拟壤中流第二层水箱：模拟壤中流第三层水箱：模拟地下径流第三层水箱：模拟地下径流、水箱结构、水箱结构 流域水文模型
18、流域水文模型水箱模型水箱模型水箱模型水箱模型(TANK)(TANK)以水箱作为蓄水容器，将降雨径流过程模拟为若干水箱的调蓄作用；模型结构简单，而且不定，水箱个数和孔数等均可以改变，参数值也不受物理概念的约束，所以适应性好，在湿润地区易于取得成功；模型的物理概念不是直接的，没有明确的土壤含水量的概念，参数多且没有一个客观的数学方法来描述，主要依靠试错法来确定参数；由于线性水库的出流没有洪峰滞时，所以计算的出流过程需要作一定时间的滞后才能与实测出流过程相符，这就相当于河网汇流。 流域水文模型流域水文模型新安江（三水源）模型新安江（三水源）模型新安江（三水源）模型新安江（三水源）模型 河海大学赵人俊
19、提出的一个水文模型，是中国少有的一个具有世界影响力的水文模型。三水源指按“山坡水文学”产流理论用一个具有有限容积和测孔、底孔的自由水蓄水库把总径流划分成饱和地面径流、壤中水径流和地下水径流。适合于湿润地区与半湿润地区水文模拟。 流域水文模型流域水文模型26 流域水文模型流域水文模型MIKE-SHEMIKE-SHE圣圣-维南方程组维南方程组 分布式分布式物理基础物理基础模型主要依据物理定律，其参数都是明确模型主要依据物理定律，其参数都是明确的物理量，可直接测定。但对模型输入信息质量和数量有较的物理量，可直接测定。但对模型输入信息质量和数量有较高要求，计算成本较高，一般较难满足。同时对复杂的应用高
20、要求，计算成本较高，一般较难满足。同时对复杂的应用条件条件( (如人类活动如人类活动) )难以适应，计算效率不高；难以适应，计算效率不高； 概念性模型计算原理简单，相对灵活，但主要参数只有抽概念性模型计算原理简单，相对灵活，但主要参数只有抽象意义，物理意义不强，多数需要反复调算和率定。另外象意义，物理意义不强，多数需要反复调算和率定。另外对具体水循环过程刻画不足；对具体水循环过程刻画不足；概念概念- -物理相结合的半分布式模型？物理相结合的半分布式模型？概念性模型和分布式模型优缺点概念性模型和分布式模型优缺点SWATSWAT模型简介模型简介（SOIL AND WATER ASSESSMENT
21、TOOL）SWAT简介开发者开发者: : Dr. Jeff Arnold for the USDA Agricultural Dr. Jeff Arnold for the USDA Agricultural Research Service (ARS)-Research Service (ARS)-美国农业部农业研究局美国农业部农业研究局应用方向应用方向: :评价土地利用管理等人类活动对流域水循环、泥评价土地利用管理等人类活动对流域水循环、泥沙、农业污染物质迁移的长期影响和作用沙、农业污染物质迁移的长期影响和作用特点：特点：1.1.属于物理概念结合的模型，具有很强的物理基础，能属于物理概念结
22、合的模型，具有很强的物理基础，能够考虑天气、土壤性质、地形、植被、人类土地管理的够考虑天气、土壤性质、地形、植被、人类土地管理的综合作用，同时能够灵活处理各种复杂应用条件；综合作用，同时能够灵活处理各种复杂应用条件；2.2.适合于长时间尺度的水文循环和物质循环研究，而非短适合于长时间尺度的水文循环和物质循环研究，而非短时期水文预报；时期水文预报；3.3.适合于宏观尺度的模拟；适合于宏观尺度的模拟；SWAT简介4.4.不仅模拟水循环过程，还能以水循环为载体，研究水土不仅模拟水循环过程，还能以水循环为载体，研究水土流失、营养物质输移、农药、病原菌等物质循环过程；流失、营养物质输移、农药、病原菌等物
23、质循环过程；5.5.能够灵活处理资料缺失问题。具有强大的模型数据库，能够灵活处理资料缺失问题。具有强大的模型数据库，除地形和土地利用等少量基本数据资料外，很多参数，除地形和土地利用等少量基本数据资料外，很多参数，如作物相关参数、土壤参数等可直接选用备用数据；如作物相关参数、土壤参数等可直接选用备用数据；6.6.分布式计算，先将流域分成子流域，子流域内继续细分分布式计算，先将流域分成子流域，子流域内继续细分水文响应单元（水文响应单元（HRUHRU），单独研究每个水文响应单元的），单独研究每个水文响应单元的内部循环，并在子流域进行累计汇总，再通过河网对子内部循环，并在子流域进行累计汇总，再通过河网
24、对子流域进行有机连接以模拟地表汇流，计算效率很高流域进行有机连接以模拟地表汇流，计算效率很高SWAT开发历史前身前身：SWRRB model (Simulator for Water Resources in Rural Basins) (Williams et al., 1985; Arnold et al., 1990)该模型为美国农业部的三个模型的集成成果该模型为美国农业部的三个模型的集成成果：CREAMS2 (Chemicals, Runoff, and Erosion from Agricultural Management Systems) (Knisel, 1980)；(PS:农
25、业管农业管理系统中的化学作用、径流和土壤侵蚀）理系统中的化学作用、径流和土壤侵蚀）EPIC4 (Erosion-Productivity Impact Calculator) (Williams et al., 1984). (ps:土壤侵蚀和生产力影响估算模型）土壤侵蚀和生产力影响估算模型）GLEAMS3 (Groundwater Loading Effects on AgriculturalManagement Systems) (Leonard et al., 1987) (ps:农业管理系统对农业管理系统对地下水的负荷影响）地下水的负荷影响）SWAT开发历史SWRRBSWRRB（Sim
26、ulator for Water Resources in Rural BasinsSimulator for Water Resources in Rural Basins） 初始时直接源于初始时直接源于CREAMSCREAMS模型的修正，扩展以下功能：模型的修正，扩展以下功能： 1) 1) 多个子流域同时演算；多个子流域同时演算； 2) 2) 地下水基流计算；地下水基流计算； 3) 3) 水库蓄滞计算水库蓄滞计算; 4) ; 4) 气象模拟与插值；气象模拟与插值； 5) 5) 提高洪峰模拟精度；提高洪峰模拟精度； 6) EPIC6) EPIC植物生长模型；植物生长模型； 7) 7) 简单的
27、洪水演算模块；简单的洪水演算模块； 8) 8) 流域产沙及运移；流域产沙及运移； 9) 9) 考虑地表径流传输损失考虑地表径流传输损失 8080年代末，又针对污染物运移模拟进行了增强，包括：年代末，又针对污染物运移模拟进行了增强，包括： 1)1)借鉴借鉴GLEAMSGLEAMS模型的农药循环模型的农药循环; 2); 2)增加增加SCSSCS曲线法模拟产流过程；曲线法模拟产流过程； 3)3)更新产沙计算公式；更新产沙计算公式；SWAT开发历史SWRRBSWRRB的主要问题：的主要问题：1 1、SWRRBSWRRB最多只能模拟最多只能模拟1010个子流域，适合于几百平方公里面积的模拟，个子流域，适
28、合于几百平方公里面积的模拟，如面积太大如面积太大( (几千平方公里以上几千平方公里以上) )会影响模拟精度；会影响模拟精度；2 2、各子流域的产出直接输出到流域出口，无河道循环过程、各子流域的产出直接输出到流域出口，无河道循环过程ROTO (Routing Outputs to Outlet) 在在8080年代末期美国印第安事务局研究亚利桑那州和新墨西哥州的印年代末期美国印第安事务局研究亚利桑那州和新墨西哥州的印第安保留区水管理对下游影响评估时开发的，是个单独的连接程序。初第安保留区水管理对下游影响评估时开发的，是个单独的连接程序。初始开发时的目的是把多个始开发时的目的是把多个SWRRBSWR
29、RB的输出连接起来，并进行河道循环演算。的输出连接起来，并进行河道循环演算。后来发现虽然效果很好，但后来发现虽然效果很好，但SWRRBSWRRB和和ROTOROTO分别独立运行很不方便，于是分别独立运行很不方便，于是进行整合。进行整合。SWAT = SWRRB + ROTOSWAT开发历史 版本进展：版本进展： SWAT90: SWAT90: 第第1 1版，融合版，融合SWRRBSWRRB和和ROTOROTO； SWAT94.2SWAT94.2：子流域划分多个水文响应单元；：子流域划分多个水文响应单元； SWAT96.2SWAT96.2：自动施肥、自动灌溉、作物模型考虑：自动施肥、自动灌溉、作
30、物模型考虑CO2CO2作用、作用、 Penman-MonteithPenman-Monteith、河水演进中考虑水质、农药传输过程；、河水演进中考虑水质、农药传输过程； SWAT98.1SWAT98.1：放牧活动、有机肥料、暗管排水、：放牧活动、有机肥料、暗管排水、 SWAT99.2SWAT99.2：重金属、城市尘埃的沉积、冲洗；：重金属、城市尘埃的沉积、冲洗； SWAT2000SWAT2000：微生物运移、：微生物运移、Green-AmptGreen-Ampt模型、模型、气象发生器气象发生器、 无限水库模拟、马斯京根法无限水库模拟、马斯京根法 SWAT2005SWAT2005：气象情景预测、
31、日降雨细化分布：气象情景预测、日降雨细化分布 SWAT2009SWAT2009：源代码和正式文档近期已经发布源代码和正式文档近期已经发布SWATSWAT流域模拟结构流域模拟结构可选可选水库也在河水库也在河网系统内网系统内SWATSWAT子流域划分及河道系统子流域划分及河道系统SWATSWAT的概念性和物理性的概念性和物理性概念性：概念性：1 1、模拟结构的概念性、模拟结构的概念性 全流域/区域子流域水文响应单元，松散耦合 物理模型：网格、地形单元等，强耦合2 2、模拟方法的概念性、模拟方法的概念性 降雨/灌溉入渗：SCS、Green-Ampt 土壤水模拟：概念性的分层下渗模型，田间持水度 地下
32、水：均衡模式，不考虑侧向径流 地表水：河槽蓄量法、马斯京根法 物理模型：Richard方程、运动波/动力波方程、地下水动力学方程物理性：物理性： 水循环各过程刻画比较全面和贴合实际，几乎所有的实际物理水文过程均有模型计算过程对应水文响应单元？水文响应单元？1 1、代表了流域的下垫面、代表了流域的下垫面2 2、是集合体的概念、是集合体的概念 土地利用 + 土壤类型 + 管理方式草地(30%)麦田(50%)林地(20%)灰潮土(70%）红壤土(30%）林地-红壤土(20%)草地-灰潮土(30%)麦田-红壤土(30%)麦田-灰潮土(20%)子子流流域域雨养灌2水灌3水林地-红壤土-雨养(20%)草地
33、-灰潮土-雨养(30%)麦田-红壤土-灌2水(30%)麦田-灰潮土-灌2水(8%)麦田-灰潮土-灌3水(12%）划分结果：划分结果：全子流域全子流域 5 个水文响应单元个水文响应单元340340个子流域，个子流域，1616种种土地利用，土地利用，2828种土种土壤，共分出壤，共分出27112711个个水文响应单元水文响应单元同一HRU子流域子流域1 1、单个水文响应单、单个水文响应单元用一维土柱模拟元用一维土柱模拟土壤水及作物过程土壤水及作物过程2 2、水文响应单元相、水文响应单元相互独立，无空间联系，互独立，无空间联系，不发生水分交换不发生水分交换与全分布式模型的区别基于网格的全分布式水文模
34、型基于网格的全分布式水文模型 - MIKE-SHE- MIKE-SHE单元格一维非饱和渗流模型饱和流动模型地表水模型两阶段的水文过程模拟 第第1 1阶段：阶段：陆面过程，控制每个子流域水、泥沙、营养物陆面过程，控制每个子流域水、泥沙、营养物微生物等的内部循环转化过程，并产出到主河道微生物等的内部循环转化过程，并产出到主河道两阶段的水文过程模拟 第第2 2阶段：阶段：河道演进过程，模拟水、泥沙、营养物等通过河道演进过程，模拟水、泥沙、营养物等通过河道演进直到流域出口河道演进直到流域出口 SWATSWAT水循环模拟路径水循环模拟路径陆面过程陆面过程河道过程河道过程水文响水文响应单元应单元主河道主河
35、道池塘池塘/湿地湿地浅层地下水浅层地下水深层地下水深层地下水SWATSWAT主要水循环模拟方主要水循环模拟方法详解法详解SWATSWAT主要过程模拟主要过程模拟1、大气水过程、大气水过程2、土壤水过程、土壤水过程3、地下水过程、地下水过程4、地表水过程、地表水过程5、植物生长过程、植物生长过程6、土地与水分管理、土地与水分管理7、泥沙过程、泥沙过程8、营养物迁移过程、营养物迁移过程9、农药迁移、农药迁移10、重金属迁移、重金属迁移11、病原菌迁移、病原菌迁移水循环过程原理水循环过程原理物质循环过程原理物质循环过程原理+一、SWAT的大气水过程大气水主过程大气水主过程2、气象站点与子流域的关系、
36、气象站点与子流域的关系1、SWAT对气象数据的管理对气象数据的管理 五大气象要素分开管理：降水五大气象要素分开管理：降水站、气温站、湿度站、辐射站、站、气温站、湿度站、辐射站、风速站风速站 采用类似泰森多边形的方法，采用类似泰森多边形的方法，子流域所用的气象数据来自与子流域所用的气象数据来自与其形心距离最近的气象站其形心距离最近的气象站一、SWAT的大气水过程气象发生器（气象发生器（WEATHER GENERATORWEATHER GENERATOR）？）？作用：作用：1 1、实际日序列气象数据有缺失时气象发生器可填补之、实际日序列气象数据有缺失时气象发生器可填补之2 2、气象资料缺乏时可用气
37、象发生器模拟的气象数据替代、气象资料缺乏时可用气象发生器模拟的气象数据替代3 3、用于情景方案计算时的气象模拟预测、用于情景方案计算时的气象模拟预测4 4、降雨过程的日内分布预测，用于、降雨过程的日内分布预测，用于Green-AmptGreen-Ampt法计算产流法计算产流- -入渗入渗 根据长测站气象数据的统计规律按照某种随机模拟方法产生日气根据长测站气象数据的统计规律按照某种随机模拟方法产生日气象数据的计算模块，五大气象要素：降雨、气温、日平均相对湿度、象数据的计算模块，五大气象要素：降雨、气温、日平均相对湿度、日辐射、日平均风速都能随机生成日辐射、日平均风速都能随机生成局限：必须提供气象
38、站气象数据的统计特征值局限：必须提供气象站气象数据的统计特征值一、SWAT的大气水过程单个气象发生器为单个气象发生器为1 1个长测气象站，需给出的信息包括：个长测气象站，需给出的信息包括：1 1、站点位置、站点位置：经：经/ /纬度，高程纬度，高程2 2、站点气象数据统计特征值：、站点气象数据统计特征值： 1 1）月最大半小时降雨的统计年数。）月最大半小时降雨的统计年数。1 1个数据个数据 2 2）多年平均每月的日最高气温。）多年平均每月的日最高气温。1212个数据个数据 3 3）多年平均每月的日最低气温。）多年平均每月的日最低气温。1212个数据个数据 4 4）多年每月的日最高气温的标准方差
39、。）多年每月的日最高气温的标准方差。1212个数据个数据 5 5）多年每月的日最低气温的标准方差。）多年每月的日最低气温的标准方差。1212个数据个数据 6 6）多年平均每月的降雨量。）多年平均每月的降雨量。1212个数据个数据 7 7）多年平均每月的日降雨量的标准方差。）多年平均每月的日降雨量的标准方差。1212个数据个数据 8 8）多年平均每月的日降雨量的偏差系数。）多年平均每月的日降雨量的偏差系数。1212个数据个数据 9 9）每月单日降雨之后第二天放晴的概率。）每月单日降雨之后第二天放晴的概率。1212个数据个数据 1010）每月单日降雨之后第二天继续降雨的概率。）每月单日降雨之后第二
40、天继续降雨的概率。1212个数据个数据 1111）多年平均每月的降雨天数。）多年平均每月的降雨天数。1212个数据个数据 1212）史上每月的最大半小时降雨量，）史上每月的最大半小时降雨量，1212个数据个数据 1313）多年平均每月的日平均辐射。）多年平均每月的日平均辐射。1212个数据个数据 1414）多年平均每月的露点温度。）多年平均每月的露点温度。1212个数据个数据 1515）多年平均每月的日风速。）多年平均每月的日风速。1212个数据个数据1 1个气象发生器个气象发生器共需给出共需给出172172个个数据信息数据信息一、SWAT的大气水过程大气水主过程大气水主过程 1 1、降雨、降
41、雪、降雨、降雪 2 2、冠层截留过程、冠层截留过程 3 3、积雪、积雪/ /融雪融雪/ /升华升华1、降雨降雪的分离，决定于当日平均气温：、降雨降雪的分离，决定于当日平均气温：降雨降雨降雪降雪T TBase?T -增加有机氮和有机磷；增加有机氮和有机磷；2 2、藻类的生长、藻类的生长-氮磷固持（包括铵氮、硝氮、无机磷）；氮磷固持（包括铵氮、硝氮、无机磷）；3 3、有机氮、有机氮/ /有机磷的底泥沉降；有机磷的底泥沉降；3 3、有机氮的矿化及底泥吸附态铵氮的释放、有机氮的矿化及底泥吸附态铵氮的释放-铵氮；铵氮；4 4、铵氮的硝化、铵氮的硝化-亚硝氮亚硝氮-硝氮；硝氮；5 5、有机磷的矿化及底泥吸
42、附态无机磷的释放、有机磷的矿化及底泥吸附态无机磷的释放-无机磷无机磷滞蓄水体中（池塘、湿地、水库滞蓄水体中（池塘、湿地、水库, ,不包括洼地）：不包括洼地）：1 1、随地表径流进入滞蓄水体，并出流进入主河道、随地表径流进入滞蓄水体，并出流进入主河道2 2、仅考虑营养元素的沉降，不模拟转化过程；、仅考虑营养元素的沉降，不模拟转化过程；农药循环叶面施药叶面施药挥发挥发分解分解冲刷冲刷入渗入渗淋溶淋溶挥发挥发产流带走产流带走分解分解地表施药地表施药农药的地表水体循环主河道中（仅模拟一种农药）：主河道中（仅模拟一种农药）：1 1、水中：降解、挥发、底泥沉降、出流；、水中：降解、挥发、底泥沉降、出流；2
43、 2、底泥中：降解、再悬浮、解吸附、掩埋；、底泥中：降解、再悬浮、解吸附、掩埋；滞蓄水体中（仅水库，过程与主河道基本一致）滞蓄水体中（仅水库，过程与主河道基本一致）病原菌运移病原菌的来源：病原菌的来源： SWAT SWAT仅概念性地模拟两种菌群，一种为可长久生存的，一种为相仅概念性地模拟两种菌群，一种为可长久生存的，一种为相对短命的。病原菌来自于绿肥对短命的。病原菌来自于绿肥( (动物粪便），动物粪便），SWATSWAT的肥料类型数据库的肥料类型数据库中有各种绿肥的含菌量参数。在施肥过程中，病原菌也一同施入，中有各种绿肥的含菌量参数。在施肥过程中，病原菌也一同施入，并在植物叶面和表层并在植物叶
44、面和表层10mm10mm的土层中进行分配。另外可通过点源输入。的土层中进行分配。另外可通过点源输入。病原菌的运移模拟：病原菌的运移模拟：HRU的病原菌分布：的病原菌分布：1 1、植物叶面；、植物叶面；2 2、土壤溶液；、土壤溶液；3 3、吸附在土壤颗粒中、吸附在土壤颗粒中1 1、植物叶面的雨水冲刷；、植物叶面的雨水冲刷；2 2、菌群的死亡、菌群的死亡- -生长生长(3(3种介质）；种介质）；3 3、土壤、土壤溶液中细菌的淋溶（淋溶的细菌认为死亡）；溶液中细菌的淋溶（淋溶的细菌认为死亡）；4 4、地表产流带走；、地表产流带走；5 5、在地表水体中（主河道、仅水库）随水流演进，并仅考虑死亡过程在地
45、表水体中（主河道、仅水库）随水流演进，并仅考虑死亡过程生化需氧量及溶解氧生化需氧量生化需氧量(CBOD) ： 指水体中有机物质分解所需要的氧气的总量。可由地表径流和点指水体中有机物质分解所需要的氧气的总量。可由地表径流和点源进入到河网系统，在主河道和水库中循环。源进入到河网系统，在主河道和水库中循环。溶解氧溶解氧(DO)：1 1、复氧过程。大气复氧、藻类的光合作用、水体扰动复氧、复氧过程。大气复氧、藻类的光合作用、水体扰动复氧( (坝体坝体) )2 2、消耗过程。藻类的呼吸作用、消耗过程。藻类的呼吸作用、CBODCBOD耗氧、底泥耗氧、铵氮硝化、耗氧、底泥耗氧、铵氮硝化、亚硝氮硝化亚硝氮硝化.
46、 .1 1、CBODCBOD的底泥沉降；的底泥沉降；2 2、CBODCBOD的氧化分解（耗氧）；的氧化分解（耗氧）；重金属的运移重金属（允许模拟重金属（允许模拟3种）：种）： SWAT SWAT的重金属来源于点源输入。重金属运移是的重金属来源于点源输入。重金属运移是SWATSWAT的物质运移模的物质运移模拟中唯一一个不考虑拟中唯一一个不考虑HRUHRU因素的物质运移，仅是简单地计算随河网系因素的物质运移，仅是简单地计算随河网系统水流演进过程和物质量平衡。统水流演进过程和物质量平衡。SWAT的优缺点总结优点：优点：1 1、开源模型。模型相关研究文档丰富、源代码容易获得；、开源模型。模型相关研究文
47、档丰富、源代码容易获得；2 2、开发时间较长，比较成熟。从、开发时间较长，比较成熟。从19901990至今至今2020余年；余年；3 3、半分布式模型。比全分布式模型易于理解和掌握，同时计算速度较快；比、半分布式模型。比全分布式模型易于理解和掌握，同时计算速度较快；比概念性模型更能体现水循环的物理机制，同时精度要好；概念性模型更能体现水循环的物理机制，同时精度要好；4 4、物质循环考虑比较全面，能模拟的物质迁移种类多，应用面广；、物质循环考虑比较全面，能模拟的物质迁移种类多，应用面广；5 5、能够充分利用土地利用等遥感信息；、能够充分利用土地利用等遥感信息；6 6、积累了较丰富的基础数据库，如
48、作物、农药、化肥等；、积累了较丰富的基础数据库，如作物、农药、化肥等；缺点：缺点：1 1、某些水循环机理刻画不尽人意，如灌溉、潜水蒸发、地下水循环等；、某些水循环机理刻画不尽人意，如灌溉、潜水蒸发、地下水循环等；2 2、涉及的模拟方向太多，初学容易找不着北；、涉及的模拟方向太多，初学容易找不着北；3 3、没有水量平衡误差分析和物质平衡误差分析；、没有水量平衡误差分析和物质平衡误差分析；4 4、模块化程度不高，源程序不易学习和修改（、模块化程度不高，源程序不易学习和修改（300300多个多个FORTRANFORTRAN源程序）；源程序）；5 5、输入文件数量巨大，不易管理。、输入文件数量巨大，不
49、易管理。SWATSWAT的输入输出的输入输出SWAT的输入文件TXTTXT文件文件水文响应单水文响应单元个数元个数* *5 + 5 + 子流域个数子流域个数* *6 + 6 + 若干点若干点源文件和气源文件和气象文件象文件模拟控制模拟控制气象数据气象数据子流域子流域HRU水库水库点源点源参数库参数库SWAT的输出文件SWATSWAT软件界面操作软件界面操作SWAT应用界面系统1 1、AVSWAT - AVSWAT - 需需 ARCVIEW ARCVIEW 平台支持；平台支持；2 2、ARCSWAT - ARCSWAT - 需需 ARCGIS ARCGIS 平台支持；平台支持；3 3、MAPWI
50、NDOW - MAPWINDOW - 需需 MAPWINDOW MAPWINDOW 平台支持平台支持http:/www.brc.tamus.edu/swat/SWAT官方网站：SWAT应用界面系统1 1、ARCVIEW3.3ARCVIEW3.3；2 2、AVSWAT-2005AVSWAT-2005；3 3、SWATSWAT相关文档。原理说明、软件操作、输入相关文档。原理说明、软件操作、输入输出说明；输出说明；4 4、SWAT2005SWAT2005源码源码软件培训材料软件培训材料 ：SWAT应用界面系统1 1、DEMDEM（数值高程数据）；（数值高程数据）；2 2、土地利用、土地利用GISGI
51、S图件；图件；3 3、土壤分布、土壤分布GISGIS图件；图件；4 4、研究区的气象数据站点。、研究区的气象数据站点。模拟必备基础数据：模拟必备基础数据：SWATSWAT的参数敏感性的参数敏感性分析、自动率定等分析、自动率定等SWAT的参数敏感性分析问题：问题： SWAT是一个涉及多个参数的综合水文是一个涉及多个参数的综合水文/水循环模拟水循环模拟模型，针对具体的一个应用过程，各个参数对模型影模型，针对具体的一个应用过程，各个参数对模型影响程度如何？哪些参数对模拟结果的影响最大？响程度如何？哪些参数对模拟结果的影响最大？敏感性分析的作用：敏感性分析的作用： 通过合理的抽样多次试算并统计，对各参
52、数的影响通过合理的抽样多次试算并统计，对各参数的影响等级进行评价并排序，使用户清楚地知道影响模型模等级进行评价并排序，使用户清楚地知道影响模型模拟结果的主要参数，有利于进一步的模型调算工作。拟结果的主要参数，有利于进一步的模型调算工作。SWAT的参数敏感性分析SWAT2005自带的敏感性分析工具：自带的敏感性分析工具：The LH-OAT sensitivity analysisLH：指指Latin-Hypercube随机采样方法，与蒙特卡洛法类随机采样方法，与蒙特卡洛法类似。将每个模型参数按其取值范围等分为似。将每个模型参数按其取值范围等分为N个取值区个取值区间，每次模拟时该参数按平均概率随
53、机采样。以等分间，每次模拟时该参数按平均概率随机采样。以等分为为10个取值区间计，模型运行个取值区间计，模型运行10次，参数在每个取值次，参数在每个取值区间内的采样概率为区间内的采样概率为1/10.如对于土壤如对于土壤Ks值，其取值范围为值，其取值范围为0100cm/hr. Latin-Hypercube随机采样过程中随机采样过程中Ks值将被分为（值将被分为（0-10）、）、（10-20）、（）、（20-30)、(90-100)共共10个取值区间个取值区间SWAT的参数敏感性分析SWAT2005自带的敏感性分析工具：自带的敏感性分析工具：The LH-OAT sensitivity analy
54、sisOAT：One-factor-At-a-Time。指每次模拟仅改变。指每次模拟仅改变1个参数个参数的值，这样便于将模拟结果的变化无偏差地归因到该的值，这样便于将模拟结果的变化无偏差地归因到该参数的变化上参数的变化上.LH-OAT结合：结合：每次仅对每次仅对1个参数进行敏感性分析，将其按固定比例调个参数进行敏感性分析，将其按固定比例调整大小整大小(如取值范围的如取值范围的5%），其他参数按），其他参数按Latin-ypercube采样方法随机变化，模拟采样方法随机变化，模拟m次并统计分析对模次并统计分析对模拟结果的影响，拟结果的影响，m为取值区间的个数。如对为取值区间的个数。如对n个参数进
55、个参数进行敏感性分析，共需行敏感性分析，共需m*(n+1)次模拟。次模拟。SWAT的参数敏感性分析共共61个参数可供率定个参数可供率定三个主要的输入控制文件：三个主要的输入控制文件：1、Senin.dat2、changepar.dat3、Responsmet.datAVSWAT界面自动选取界面自动选取27个固定参数进行敏感性分析个固定参数进行敏感性分析（280次模拟），如有特殊次模拟），如有特殊要求，需自己手动修改这三要求，需自己手动修改这三个文件。个文件。SWAT的参数敏感性分析LH-OAT基本参数设置基本参数设置Sensin.datSWAT的参数敏感性分析LH-OAT基本参数设置基本参数设
56、置changepar.dat针对随机变化参数的变动方法SWAT的参数敏感性分析LH-OAT基本参数设置基本参数设置responsmet.dat敏感性分析对象及评价方法控制。可多个对象分析。敏感性分析对象及评价方法控制。可多个对象分析。每个一行数据控制。每个一行数据控制。对象编号评价方法评价方法2阈值临时文件号水质选项SWAT的参数敏感性分析LH-OAT基本参数设置基本参数设置敏感性分析对象为进行敏感性分析对象为进行敏感性评判的数据对象：敏感性评判的数据对象：流量流量(1)？ 产沙量产沙量(2)?有机氮有机氮(3)？有机磷有机磷(4)？SWAT的参数敏感性分析LH-OAT敏感性分析主要敏感性分析
57、主要输出文件：输出文件：SWAT的参数自动率定/不确定性分析参数自动率定（参数寻优）参数自动率定（参数寻优）目的：在模型各参数可行的取值范围内，寻找一套确目的：在模型各参数可行的取值范围内，寻找一套确定的参数，使得模拟值与实测值的整体差别最小。定的参数，使得模拟值与实测值的整体差别最小。SWAT自带的参数寻优/不确定性分析方法：SCE-UA(Shuffled Complex Evolution)算法目标函数：目标函数： 该算法结合了单纯形法、受控随机搜索、生物竞争进化和种群交叉等方法的优点。算法第一步(第零个循环)先运用随机抽样在所有可行参数空间中选择一个初始种群。该可行空间为由changep
58、ar.dat文件指定的参数及其变化范围。随后，该种群被分割成多个个体(complexs)，对应P个参数来说就得到2p+1个点。每个个体运用单纯形法进行独立进化，个体之间定期进行交叉形成新的个体，从而可以获得更多的优化信息。该算法可以搜索全部参数的可行空间，找到全局最优参数的成功率是100%。SCE-UA算法的特点：算法的特点：SWAT的参数自动率定/不确定性分析SWAT的参数自动率定/不确定性分析不确定性分析不确定性分析目的：由于参数有不确定性因素，因此给出某一置信目的：由于参数有不确定性因素，因此给出某一置信区间内（如区间内（如90%、95%、97.5%）模拟结果的分布范）模拟结果的分布范围
59、，供用户参考。围，供用户参考。SCE-UA主要输入文件：主要输入文件：SWAT的参数自动率定/不确定性分析SCE-UA主要输出文件：主要输出文件：SWAT的参数自动率定/不确定性分析SWAT-CUPSWAT-CUP 是一个用于SWAT率定验证的计算机程序，为开源程序，可以自由使用和拷贝。该程序将GLUE、ParaSol，SUFI2，MCMC、PSO等处理过程与SWAT关联起来，用于参数敏感性分析（sensitivity analysis） 、参数率定（calibration）、模型验证（validation）、不确定性分析（uncertainty analysis）SWAT-CUP 特色: 并
60、行处理可视化界面SWATSWAT应用研究示例讲解应用研究示例讲解SWAT应用研究示例1、径流预测2、ET模拟3、土地利用变化4、气候变化5、水土保持6、水资源评价7、融雪产流8、水资源管理9、水量平衡分析10、模型耦合研究11、模型改进研究12、参数敏感性分析研究13、流域面源污染SWAT应用类1-径流预测1袁军营等，基于SWAT模型的柴河水库流域径流模拟研究J,北京师范大学学报(自然科学版),2006. 柴河水库流域位于辽宁省铁岭市，是辽河中游左侧支流柴河上的一座以防洪、灌溉、工业和城市供水为主，兼顾发电、养鱼等综合利用的多年调节大()型水利枢纽工程。总库容6.36亿m3，其控制流域面积13
61、55 km2，流域多年平均降水量737 mm.SWAT应用类1-径流预测1袁军营等，基于SWAT模型的柴河水库流域径流模拟研究J,北京师范大学学报(自然科学版),2006. 利用20022003年的实测水量、降水数据进行参数率定,应用20042005年的实测数据进行模型验证。模拟结果表明，SWAT模型在柴河水库流域径流模拟中的适用性和可靠性良好，可为流域水文过程的预测提供可靠的模型基础SWAT应用类1-径流预测2宋轩等，SWAT模型在淅川县丹江口库区的应用研究J,郑州大学学报(工学版),2010. 淅川县位于河南省西南边陲,豫、鄂、陕三省的结合部，面积2820 km2，是南水北调中线工程渠首所
62、在地。属北亚热带向暖温带过渡的季风性气候，多年平均降水量804 mm。其主要河流有丹江、鹳河、淇河、滔河、刁河。丹江及其支流域面积占全县面积的93. 5%.SWAT应用类1-径流预测2宋轩等，SWAT模型在淅川县丹江口库区的应用研究J,郑州大学学报(工学版),2010. 通过表1对模拟结果与实测值进行比较，可以看出模型预测相对误差控制在15%的范围内，平均精度达到91%，纳西效率系数也高达0. 88，表明年径流的预测值与观测值变化基本一致，说明模型能够较好地模拟丹江口库区年径流的变化SWAT应用类2-ET模拟钱坤等，基于SWAT模型的房山区不同情景方案下的蒸发蒸腾模拟J. 农业工程学报，201
63、1. 房山区位于北京市西南部，全区的总面积2019 km2，山地和丘陵面积占65.3%，平原面积占34.7%。地处华北平原与太行山交界地带，处于暖温带半湿润地区多年平均降雨量为587.6 mm。SWAT应用类2-ET模拟钱坤等，基于SWAT模型的房山区不同情景方案下的蒸发蒸腾模拟J. 农业工程学报，2011. 耕地园林、灌木林/绿地、荒草地为三个占主导地位的土地利用，占总面积的87%。SWAT应用类2-ET模拟钱坤等，基于SWAT模型的房山区不同情景方案下的蒸发蒸腾模拟J. 农业工程学报，2011. 模型模拟的房山区20022007年蒸腾蒸发值除2003年与遥感监测值相差5%以上外，其余年分均
64、在5%以内。平均值相差为1.9%。不同土地利用情况下的ET值与遥感监测值对比平均相差误差为6.01%。表明模型具有较好的适用性。不同土地利用不同年份区域整体SWAT应用类2-ET模拟钱坤等，基于SWAT模型的房山区不同情景方案下的蒸发蒸腾模拟J. 农业工程学报，2011.以以2020年和年和2030年为水平年，情景方案设置年为水平年，情景方案设置1、考虑节水灌溉措施下的农业灌溉定额管理；2、考虑农业种植结构优化；3、考虑逐步减少地下水超采，2020年零超采；4、考虑南水北调来水(1.2亿方）种植结构调整数据农业节水灌溉设置8个情景SWAT应用类2-ET模拟钱坤等，基于SWAT模型的房山区不同情
65、景方案下的蒸发蒸腾模拟J. 农业工程学报，2011. 以2004年为现状年，根据未来2020、2030年房山区水平衡，设定目标ET为546.1 mm。将按照8种情景模拟的全年ET结果与目标ET进行比较，从结果看，预测2020、2030年，在情景3（即中度节水、种植结构为优化及地下水超采逐步减少）时模拟的ET值与目标ET值最接近评论：似乎各个情景区别不大？是不是与耕地面积比重较小有关？SWAT应用类3-土地利用变化1庞靖鹏，刘昌明等，密云水库流域土地利用变化对产流和产沙的影响J. 北京师范大学学报（自然科学版），2010. 研究区为密云水库上游流域。密云水库流域是指潮白河流域密云水库所控制的部分
66、， 流域面积为15369 km2，其中白河流域面积约为9197 km2，潮河流域面积为6172 km2。流域内多年平均降水量为488.9 mm.SWAT应用类3-土地利用变化1庞靖鹏，刘昌明等，密云水库流域土地利用变化对产流和产沙的影响J. 北京师范大学学报（自然科学版），2010. 先做现状模拟。以19861991年共5年进行连续模拟。率定期(19861988)和验证期(19891991)下会和张家坟2个站点的模拟效果都非常好，纳西效率系数达到0.870.95SWAT应用类3-土地利用变化1庞靖鹏，刘昌明等，密云水库流域土地利用变化对产流和产沙的影响J. 北京师范大学学报（自然科学版），20
67、10. 收集密云水库流域内90年代初、90年代中期、90年代末期三期土地利用。通过分析，地类转移主要在耕地、林地和草地间进行。以草地和耕地向林地、林地和耕地向草地转变最为剧烈。SWAT应用类3-土地利用变化1庞靖鹏，刘昌明等，密云水库流域土地利用变化对产流和产沙的影响J. 北京师范大学学报（自然科学版），2010. 将90年代初、中和末期3期的土地利用图分别输入SWAT模型中,对密云水库流域土地利用变化对密云水库径流和产沙的影响作出定量的分析.（降水气象数据沿用19861991年系列）与1990年相比，1995年的土地利用条件下径流量略有下降，然而泥沙量却有大幅度的减少。说明采取退耕还林、还草
68、等措施将极大地减少流域内的水土流失潮河白河宋艳华等，SWAT模型辅助下的生态恢复水文响应-以陇西黄土高原华家岭南河流域为例J. 生态学报，2008.SWAT应用类3-土地利用变化2华家岭南河流域位于陇西黄土高原的中部地区，为葫芦河一级支流。面积约为1125km2， 流域年平均气温38。多年平均降水量在400550mm之间。宋艳华等，SWAT模型辅助下的生态恢复水文响应-以陇西黄土高原华家岭南河流域为例J. 生态学报，2008.SWAT应用类3-土地利用变化2对模型进行校准与验证，将19811983年和19851989年分别作为模型的校准期与验证期。年均径流深的模拟值与实测值的相对误差在5%以内
69、, Nash效率系数均在0. 6以上。年过程月过程宋艳华等，SWAT模型辅助下的生态恢复水文响应-以陇西黄土高原华家岭南河流域为例J. 生态学报，2008.SWAT应用类3-土地利用变化2情景分析方案宋艳华等，SWAT模型辅助下的生态恢复水文响应-以陇西黄土高原华家岭南河流域为例J. 生态学报，2008.SWAT应用类3-土地利用变化2结论：结论：1、模拟年均径流深裸地耕地林地草地，草地和森林植被更涵养水源，同时草地比森林植被涵养水源的作用更强； 2、模拟年均蒸散发草地林地裸地耕地, 草地的蒸散发量最大,比林地高出22%。（注：耕地最小？）不同情景径流深度不同情景平均蒸发SWAT应用类3-土地
70、利用变化3张圣威等，土地覆被和气候变化对拉萨河流域径流量的影响J. 水资源保护，2010.流域面积2.6万km2，年均降水450mm（半湿润区），年均气温7.5度。张圣威等，土地覆被和气候变化对拉萨河流域径流量的影响J. 水资源保护，2010.SWAT应用类3-土地利用变化3先做现状率定验证张圣威等，土地覆被和气候变化对拉萨河流域径流量的影响J. 水资源保护，2010.SWAT应用类3-土地利用变化3设置4个土地变化情景进行分析不同情景下的径流量张圣威等，土地覆被和气候变化对拉萨河流域径流量的影响J. 水资源保护，2010.主要结论：1、气候变化对年径流量的影响显著，若温度升高2度，降水增加2
71、0%（年均降水450mm)，径流将增加89%；2、土地利用变化最多影响径流10%；SWAT应用类3-土地利用变化3再考虑气候变化朱利等，基于径流模拟的汉江上游区水资源对气候变化响应的研究J. 资源科学，2005.SWAT应用类4-气候变化1汉江上游两河流域，积水面积2816km2，降水730-1230mm（湿润区），年均气温12度朱利等，基于径流模拟的汉江上游区水资源对气候变化响应的研究J. 资源科学，2005.SWAT应用类4-气候变化1先做现状率定验证朱利等，基于径流模拟的汉江上游区水资源对气候变化响应的研究J. 资源科学，2005.SWAT应用类4-气候变化1现状平均量各情景平均量进行比
72、较朱利等，基于径流模拟的汉江上游区水资源对气候变化响应的研究J. 资源科学，2005.SWAT应用类4-气候变化1主要结论：1、径流：温度增加4度，径流减少4%；降水增加10%，径流增加12%；2、蒸发：温度增加4度，蒸发变化7%；降水增加10%，蒸发变化2%。张利平等，南水北调中线工程水源区水文循环过程对气候变化的响应J. 水利学报，2010.SWAT应用类4-气候变化2摘要：根据联合国政府间气候变化专业委员会第四次评估报告中大摘要：根据联合国政府间气候变化专业委员会第四次评估报告中大气环流模型多模式输出结果，分析了特别情境下气环流模型多模式输出结果，分析了特别情境下2121世纪降水、气温、
73、世纪降水、气温、径流、蒸发的响应过程径流、蒸发的响应过程流域面积9.5万km2，平均降水820mm，平均气温14.6度。张利平等，南水北调中线工程水源区水文循环过程对气候变化的响应J. 水利学报，2010.SWAT应用类4-气候变化2GCMsGCMs模式数据：模式数据：采用世界气候研究计划（采用世界气候研究计划（World Climate World Climate Research Programme, WCRP)Research Programme, WCRP)的耦合模式比较计划的耦合模式比较计划- -阶段阶段3 3的多模式数的多模式数据。未来气候情景预估采用据。未来气候情景预估采用IPC
74、CIPCC第四次评估报告中假定的第四次评估报告中假定的3 3个情景：个情景：1 1、A2A2情景，高排放，注重经济增长的区域发展情景；情景，高排放，注重经济增长的区域发展情景；2 2、A1BA1B情景，中等排放，注重经济增长的全球共同发展情景；情景，中等排放，注重经济增长的全球共同发展情景；3 3、B1B1情景，低排放，注重全球可持续发展情景；情景，低排放，注重全球可持续发展情景；气候数据降尺度方法：气候数据降尺度方法：将未来将未来2011-20992011-2099年划分为年划分为3 3个时期，个时期，2011-2040,2041-2070,2071-2099.2011-2040,2041-
75、2070,2071-2099.利用利用SWATSWAT的气象发生器的气象发生器将月尺度将月尺度资料降解为日尺度资料，并输入资料降解为日尺度资料，并输入SWATSWAT模型中。模型中。张利平等，南水北调中线工程水源区水文循环过程对气候变化的响应J. 水利学报，2010.SWAT应用类4-气候变化2降水状况：降水状况：A2A2、A1BA1B、B1B1情景年降水变化率分别为情景年降水变化率分别为1.59%1.59%，1.41%1.41%，0.65%0.65%；气温状况：气温状况：A2A2、A1BA1B、B1B1情景分别增长情景分别增长4.334.33、3.783.78、2.192.19度；度；201
76、1-20992011-2099年间气候预测：年间气候预测：（基准期：（基准期：1961-19901961-1990）降水变化气温变化张利平等，南水北调中线工程水源区水文循环过程对气候变化的响应J. 水利学报，2010.SWAT应用类4-气候变化2径流预测：先减小后增大。径流预测：先减小后增大。A2A2、A1BA1B、B1B1情景在情景在2080S2080S期间分别增大期间分别增大13%13%、7%7%、10%10%蒸发预测：蒸发预测：A2A2、A1BA1B、B1B1情景比较接近，约增加情景比较接近，约增加4-5%,4-5%,且增加幅度有减小趋势；且增加幅度有减小趋势；2011-20992011
77、-2099年间相对基准期年间相对基准期径流变化蒸发变化王林等，退化山地生态系统植被恢复水文效应的SWAT模拟J. 山地学报，2008.SWAT应用类5-水土保持1晋江流域安溪水文站上游，集水面积晋江流域安溪水文站上游，集水面积2466km22466km2，年降水，年降水1715mm1715mm王林等，退化山地生态系统植被恢复水文效应的SWAT模拟J. 山地学报，2008.SWAT应用类5-水土保持1情景设定主要针对农用地和林草地，进行流域内植被恢复的产流产沙模拟王林等，退化山地生态系统植被恢复水文效应的SWAT模拟J. 山地学报，2008.SWAT应用类5-水土保持1结论：各情景模拟结果都为减
78、水减沙。产流的相对变化系数在1%之下，产沙相对变化系数为2%-8%。还林的减沙效果优于还草。影响不大？是否本身农用地面积较小？彭顺风等，基于数字流域的地表水资源模拟J. 水资源保护，2009.SWAT应用类6-水资源评价1研究区为史灌河流域，面积5930km2，年降水量1100mm彭顺风等，基于数字流域的地表水资源模拟J. 水资源保护，2009.SWAT应用类6-水资源评价1评价思路：评价思路：通过逐项还原法将水文站的实测值还原为水资源的真值，然后模型仅模拟自然过程下的径流量与之对比甄婷婷等，蓝水绿水资源量估算方法及时空分布规律研究-以卢氏流域为例J. 资源科学，2010.SWAT应用类6-水
79、资源评价2蓝水：传统意义上的水资源量，为地表水资源量与地下水资源量之 和减去重复量绿水：源于降水，存储于土壤并被植被蒸散发消耗的水资源位于黄河中下游，流域面积4655km2，多年平均气温14度，降水600-1100mm甄婷婷等，蓝水绿水资源量估算方法及时空分布规律研究-以卢氏流域为例J. 资源科学，2010.SWAT应用类6-水资源评价2利用利用SWATSWAT模型进行蓝绿水资源评价：模型进行蓝绿水资源评价：蓝水资源：等于产水量与深层地下水补给之和绿水资源：绿水流等于总蒸散量，绿水储量等于给定时间内的土壤含水量甄婷婷等，蓝水绿水资源量估算方法及时空分布规律研究-以卢氏流域为例J. 资源科学，2
80、010.SWAT应用类6-水资源评价2将计算结果反映到GIS图上展示结果李慧等，基于SWAT模型的山区冰雪融水河流的日径流模拟研究J. 灌溉排水学报，2010.SWAT应用类7-融雪产流结论：冬季径流极小，春季过后径流明显增大，表明冬季气温低，冰雪融水极少；春季气温回升，冰雪融水对玛纳斯河径流产生了很大影响，说明在终年积雪的玛纳斯河山区，冰雪融水对其1年的径流起着关键的调控作用杨威等，基于改进SWAT模型的县级综合管理工具在项目县IWEMP中的应用J，海河水利，2011.SWAT应用类8-水资源管理摘要摘要：基于SWAT2000源代码，面向县级管理人员开发了县级水资源与水环境综合管理工具。能够
81、基于项目县DEM、土壤、土地利用等基本信息，构建项目县分布式水循环模型，实现对项目县地表水、地下水以及水质时空变化过程的模拟以及对不同规划方案包括节水措施（调整作物种植结构、优化灌溉水量、调整灌溉水源等）和减污措施（点源、面源）的情景预测，模拟分析结果为县级水资源与水环境综合管理规划的制定与选取提供数据支撑。相当于为SWAT重新开发针对特定应用的界面杨威等，基于改进SWAT模型的县级综合管理工具在项目县IWEMP中的应用J，海河水利，2011.SWAT应用类8-水资源管理情景模拟情景模拟1：作物种植结构调整：作物种植结构调整以成安县为例以成安县为例将现状成安县小麦玉米轮作（高耗水）种植面积的2
82、0%、40%、60%和80%改为棉花种植结论：将小麦玉米轮作面积的20%、40%、60%和80%改为棉花种植后，全县ET平均下降了1.96%、3.92%、5.89%和7.85%。杨威等，基于改进SWAT模型的县级综合管理工具在项目县IWEMP中的应用J，海河水利，2011.SWAT应用类8-水资源管理情景模拟情景模拟2：优化灌溉水量：优化灌溉水量以馆陶县为例以馆陶县为例结论：减少灌溉量对于降低ET作用显著，在减少20%和50%的灌溉量后，ET可以降低的百分比分别为2.3%和12.6%。原馆陶县现状灌水定额为70mm，在此基础上分别将灌溉水量减少20%和50%杨威等，基于改进SWAT模型的县级综
83、合管理工具在项目县IWEMP中的应用J，海河水利，2011.SWAT应用类8-水资源管理情景模拟情景模拟3：调整灌溉水源：调整灌溉水源以肥乡县为例以肥乡县为例结论：外流域调水主要用来补充有咸水区深层地下水开采，因此以旧店乡、东漳堡乡为代表的乡镇深层地下水位在外流域调水后得到明显回升。将肥乡县现状条件下部分开采深层地下水进行农业灌溉的子流域（主要是有咸水区）改为外流域调水灌溉，外流域调水量约3 000万m3。朱新军等，基于分布式模拟的流域水平衡分析研究以海河流域为例J，地理科学进展，2008.SWAT应用类9-流域水量平衡分析将海河流域划分为283个子流域和2010个水文响应单元(HRU)，概化
84、入海河口17个，考虑大、中型水库48座海河流域总面积31.8万km2，其中山地和高原面积18.9万km2，约占60%；平原面积12.9万km2，约占40%。属于温带半湿润、半干旱大陆性季风气候区。多年平均降水量一般为500650mm朱新军等，基于分布式模拟的流域水平衡分析研究以海河流域为例J，地理科学进展，2008.SWAT应用类9-流域水量平衡分析在子流域水量模拟基础上，统计出15个水资源三级分区的降水、蒸发及径流的空间分布，给出了2003、2004、2005三个年份对比变化。朱新军等，基于分布式模拟的流域水平衡分析研究以海河流域为例J，地理科学进展，2008.SWAT应用类9-流域水量平衡
85、分析V=P-E-R+Qin-Qout流域总的水量平衡流域总的水量平衡王宏等，SWAT/GMS联合模型在华北平原地下水库研究中的应用J，世界地质，2005.SWAT应用类10-模型耦合研究研究区位于河北平原子牙河流域。总面积约3.9万km2。根据1979年DEM数据，提取平原区子流域71个、山区子流域4个，其中山区子流域作为平原区侧向补给量的分区基础，平原区71个子流域作为计算地下水入渗量的单元划分王宏等，SWAT/GMS联合模型在华北平原地下水库研究中的应用J，世界地质，2005.SWAT应用类10-模型耦合研究1SWAT：以流域为单元，由SWAT根据水均衡方程，分别算出本区1979年71个子
86、流域水资源量，将计算出的流域内的地下水渗入量直接输入GMS模型中；王宏等，SWAT/GMS联合模型在华北平原地下水库研究中的应用J，世界地质，2005.SWAT应用类10-模型耦合研究1GMS：将地表小流域再次按岩性分区，降雨入渗量在此分区中进行分配，其分配结果在GMS中直接输入。其他汇项则根据统测的数据，结合实际用水情况分配到每个区中。地下水边界：地下水边界：山前边界认为是定流量边界，侧向补给量为19.51亿m3(怎么来的？)；入海处为定水头边界，水头值为0；他界因考虑到与地下水流向大致平行，处理为零流量边界；王宏等，SWAT/GMS联合模型在华北平原地下水库研究中的应用J，世界地质，200
87、5.SWAT应用类10-模型耦合研究1应用全区49个观测孔1979年112月的实际观测水位，对模型进行了调参、校核、验证。模型运行稳定，残差范围合理。所建模型参数能真实地表达地下水系统的地质状况张雪刚等，SWAT模型与MODFLOW模型的耦合计算及应用J，水资源保护，2010.SWAT应用类10-模型耦合研究2研究区：张集地区，位于江苏省徐州市东南,距徐州市中心25 km，总面积360 km2。张雪刚等，SWAT模型与MODFLOW模型的耦合计算及应用J，水资源保护，2010.SWAT应用类10-模型耦合研究2方法：构建SWAT中HRU和MODFLOW中有限差分网格的HRU-cell的交互界面
88、见下参考文献，将SWAT模型计算的模型要素的空间分布赋值给MODFLOW，从而实现SWAT与MODFLOW的耦合计算。 NAM W K,MOON C, YOO S W,et al.Development and application of the integrated SWAT-MODFLOW modelJ. Journal of Hydrology, 2008, (356):1-16.张雪刚等，SWAT模型与MODFLOW模型的耦合计算及应用J，水资源保护，2010.SWAT应用类10-模型耦合研究2第第1步：步：根据地区土地利用和土壤类型，共划分12个子流域。由于区域面积较小，且土地利用
89、以农田为主，土地开发利用方式及土壤覆盖情况比较单一，因此HRU划分采用优势土被占有法，每个子流域对应1个HRU，共分为12个HRU(见前)；第第2步：步：选用地区2005年的实测降雨和蒸发资料，进行SWAT月模拟计算，得到区域内每个HRU的地下水补给量地下水补给量和潜水蒸发量潜水蒸发量值；第第3步：步：应用ARCGIS软件，将SWAT模型计算的地下水补给量和潜水蒸发量结果引入MODFLOW模型边界条件中的地下水补给子模块RCH和潜水蒸发子模块EVT中，进行地下水模拟计算。 根据张集地区的水文地质条件根据张集地区的水文地质条件,建立区域的建立区域的MODFLOW地下水模型。地下水模型。模型分为模
90、型分为3层，模型网格划分为每层层，模型网格划分为每层46行行60列，每个单元格列，每个单元格(cell)的大的大小为小为500m500m。子流域形状与网格不兼容，如果和处理的？张雪刚等，SWAT模型与MODFLOW模型的耦合计算及应用J，水资源保护，2010.SWAT应用类10-模型耦合研究2区域内35口地下水观测井2005年4月5日的实测地下水位与模型计算地下水位的散点图；确定性系数为078，相关系数R2为081第7个HRU内的孔隙含水层一观测井实测与计算地下水水位过程线。计算得到的地下水水位随地下水补给量的变化趋势与实测地下水位的变化基本吻合(?)张永勇等，基于SWAT模型的闸坝水量水质优
91、化调度模式研究J，水力发电学报，2010.SWAT应用类10-模型耦合研究3 以北京市温榆河流域为例，建立了以北京市温榆河流域为例，建立了SWAT与闸坝调度模型的耦合模型与闸坝调度模型的耦合模型(SWAT-QCmode)，从流域的角度探讨了闸坝群的水量水质联合调度模，从流域的角度探讨了闸坝群的水量水质联合调度模式，为温榆河流域水污染整治，生态修复提供技术支撑式，为温榆河流域水污染整治，生态修复提供技术支撑SWAT模型调度模型张永勇等，基于SWAT模型的闸坝水量水质优化调度模式研究J，水力发电学报，2010.SWAT应用类10-模型耦合研究3闸坝调度模型闸坝调度模型(QCmode)：目标函数：目
92、标函数：1、河流水资源可利用量最大；2、河流受污染程度最小约束条件：约束条件：1、闸坝下泄流量和防洪约束；2、库容约束耦合方法：耦合方法：采用Fortran语言编译QCmode模型，然后将其嵌入SWAT源程序中。二者之间以文本读入读出形式进行软连接，QCmode模型输出为所需优化闸坝的月出流过程，作为SWAT模型的输入；SWAT模型演算得到的目标闸坝的水量水质过程作为QCmode模型的输入读入，对所需优化闸坝的月出流过程进行筛选张永勇等，基于SWAT模型的闸坝水量水质优化调度模式研究J，水力发电学报，2010.SWAT应用类10-模型耦合研究3温榆河流域面积2478 km。共分36个子流域，1
93、26个水文响应单元，河道上修建闸坝10个。张永勇等，基于SWAT模型的闸坝水量水质优化调度模式研究J，水力发电学报，2010.SWAT应用类10-模型耦合研究3现状年的率定验证。2004年为率定期，2005年为验证期铵氮COD流量张永勇等，基于SWAT模型的闸坝水量水质优化调度模式研究J，水力发电学报，2010.SWAT应用类10-模型耦合研究3鲁曈闸新堡闸苇沟闸北关闸结论：通过调度，温榆河干流水量水质状况有明显改善SWAT应用类11-模型改进研究1仕玉治等，SWAT模型在水稻灌区的改进及应用研究J,水电能源科学,2010.主要改进原因：1、原模型仅能进行单一水源取水灌溉,不能考虑流域内多水源
94、灌溉情况即地下水和地表水联合取水灌溉；2、原模型中自动灌溉是以田间含水量判断最大灌溉水量,仅适用于旱田,而水稻不同生长阶段具有不同蓄水层深度研究区：努敏河流域长阁灌区，流域面积为3310km2, 年平均降雨量为534mmSWAT应用类11-模型改进研究1仕玉治等，SWAT模型在水稻灌区的改进及应用研究J,水电能源科学,2010.SWAT应用类11-模型改进研究1仕玉治等，SWAT模型在水稻灌区的改进及应用研究J,水电能源科学,2010.方案方案0：利用原模型中单一水源和自动灌溉进行计算的方案;方案方案1：依据取水规则进行计算的方案，为地表水和地下水混合灌溉，该方案为改进模型用于计算灌溉水量和原
95、模型性能比较;方案方案2：改进模型中不考虑灌溉的方案，为河道水源单一灌溉，与方案1比较分析灌溉取用水对河道径流的影响。结论：方案1较方案0的模拟结果精度有所提高,考虑农业灌溉的方案1和方案0均比不考虑灌溉的方案2模拟结果好,表明农业灌溉用水对水循环存在一定的影响SWAT应用类11-模型改进研究2郑捷等，改进的SWAT模型在平原灌区的应用J,水利学报, 2011.主要改进原因：1、平原区的地形高差使模型通常无法提取完整的自然河网；2、灌区是以输配水渠系与排水沟网覆盖的灌域为单元进行用水管理，而非集水子流域；3、灌溉水量只能以田持为上限，多余退回；4、模型中的叶面积指数是根据美国作物的生长情况给出
96、的一系列默认值，这不适合于中国灌区作物实际腾发量的计算研究；研究区：山西汾河灌区，面积为1366km2, 年平均降雨量为455mm。SWAT应用类11-模型改进研究2郑捷等，改进的SWAT模型在平原灌区的应用J,水利学报, 2011.改进方法：1、根据研究区自然河道与人工排水沟网的空间分布，运用“burn-in”算法，生成新的流域河网；2、对于人工输配水渠道，借鉴高程增量迭加算法，增加渠道上栅格单元的高程，将数字渠系叠加到经河网修正过的DEM上。以骨干渠道与排水干沟覆盖的区域为划分标准，设置合理的最小水道集水区面积阈值划分子流域，干支渠道作为子流域划分的边界处理(不模拟输配水损失)3、灌溉水量
97、从土壤表层开始，使各层含水量分别达到田持，当最底层土壤达到田持后，剩余的灌水量认为补充地下水；4、将SWAT模型中实际腾发量计算程序进行了改进，采用国内常用的作物系数法计算农作物实际腾发量，计算方法为：SWAT应用类11-模型改进研究2郑捷等，改进的SWAT模型在平原灌区的应用J,水利学报, 2011.根据前述子流域与水文响应单元的划分原则，灌区共划分73个子流域和324个水文响应单元，子流域面积从231016hm2，水文响应单元面积为1.57191hm2SWAT应用类11-模型改进研究2郑捷等，改进的SWAT模型在平原灌区的应用J,水利学报, 2011.逐月与逐日模拟的结果满足模型评价要求，
98、纳西效率系数分别为0.72、0.65，表明构建的模型适应于平原型灌区的水量平衡模拟SWAT应用类12-参数敏感性分析研究1雷晓辉等，北京市Swat模型参数自动率定J, GEF项目国际研讨会中文论文集, 2009.问题：问题：SWAT2005版本中新增了LH-OAT灵敏度分析模块和SCE-UA自动校准模块，计算耗费的时间随着子流域个数的增加而增加。以北京市为例，共划分了138个子流域，连续5年的日过程进行模拟大概仅需要一分钟时间。但对10个参数进行率定，模型要重复运行5000次，大概需要4天的时问。因此，这种参数率定方法对大流域来说是不合适的，不能令人满意的。解决办法：解决办法：只对部分子流域进
99、行模拟，如只模拟目标水文站控制的上游的几个子流域，参数自动率定需要的时间就会大大减少。在把该方法应用到北京市的SWAT模型参数率定时，对选定的前辛庄站的S年的同径流过程进行率定时总共用时不超过30分钟，而且率定效果较好，说明了该方法的高效性。SWAT应用类12-参数敏感性分析研究1雷晓辉等，北京市Swat模型参数自动率定J, GEF项目国际研讨会中文论文集, 2009.SWAT应用类12-参数敏感性分析研究1雷晓辉等，北京市Swat模型参数自动率定J, GEF项目国际研讨会中文论文集, 2009.问题：如果水文站为流域出口，就没法偷懒了基于基于SWATSWAT原理的改进原理的改进水循环模拟模型
100、水循环模拟模型MODCYCLEMODCYCLEMODCYCLE模型主要特点借鉴SWAT模拟思想开发的半分布式水文/水循环模型基于面向对象的c+语言，完全模块式开发(2.5万行代码)基于ACCESS数据库平台开发，输入输出均基于数据库具有较完善的水平衡校验机制支持并行运算MODCYCLE: An Object Oriented Modularized Model for Basin Scale Water Cycle Simulation 开发历时开发历时3 3年多，目前水循环已经开发完成并稳定，在邯郸市、天津年多，目前水循环已经开发完成并稳定，在邯郸市、天津市、浙江铜山源灌区、内蒙古通辽市等进
101、行了应用。市、浙江铜山源灌区、内蒙古通辽市等进行了应用。MODCYCLE模型MODCYCLE模型的模型的26 C+ 模块模块MODCYCLE模型层次化的水平衡校验机制层次化的水平衡校验机制1 1、基本单元模拟结果输出表、基本单元模拟结果输出表(31(31项输出项输出) )2 2、水库模拟结果输出表水库模拟结果输出表(10(10项输出项输出) )3 3、湿地模拟结果输出表、湿地模拟结果输出表(7(7项输出项输出) )4 4、池塘模拟结果输出表池塘模拟结果输出表(8(8项输出项输出) )5 5、地下水模拟结果输出表、地下水模拟结果输出表(15(15项输出项输出) )6 6、主河道模拟结果输出表、主
102、河道模拟结果输出表(14(14项输出项输出) )7 7、子流域模拟结果输出表、子流域模拟结果输出表(43(43项输出项输出) )8 8、全流域模拟结果输出表、全流域模拟结果输出表(33(33项输出项输出) )水循环模拟实体层水循环模拟实体层次独立水量校核次独立水量校核子流域层次综合水量校核子流域层次综合水量校核全流域层次综合水量校核全流域层次综合水量校核MODCYCLE模型MODCYCLE模型双核心并行计算双核心并行计算MODCYCLE模型4核八线程并行计算（核八线程并行计算（Core i7 920）MODCYCLE模型常规（非并行）常规（非并行）运算耗时运算耗时并行运算并行运算耗时耗时4核心
103、八线程下速度提升核心八线程下速度提升3.3倍（倍（Core i7 920）MODCYCLE模型模拟功能的增强模拟功能的增强1.1.地表滞蓄功能，解决平原区产流地表滞蓄功能，解决平原区产流/ /入渗模拟问题；入渗模拟问题；2.2.无限量灌溉，解决无限量灌溉，解决SWATSWAT最多只灌溉到田间持水度的问题；最多只灌溉到田间持水度的问题；3.3.多水源灌溉，解决多水源灌溉，解决SWATSWAT水文响应单元只能单一水源灌溉的问题；水文响应单元只能单一水源灌溉的问题；4.4.增加灌溉损失处理；增加灌溉损失处理；5.5.结合灌溉制度的动态灌溉；结合灌溉制度的动态灌溉；6.6.地下水模拟修正，解决地下水模拟修正，解决SWATSWAT地下水位无物理意义，潜水蒸发不通过地下水位无物理意义，潜水蒸发不通过土壤层的问题，并可进行数值模拟；土壤层的问题，并可进行数值模拟；7.7.湿地补水功能；湿地补水功能；MODCYCLE模型水资源形成转化定量评价(天津，973)水量调控下水循环响应分析(邯郸市，十一五支撑)干旱区水-生态耦合模拟（通辽，横向)灌区水平衡分析(铜山源，水利部公益项目)MODCYCLE模型MODCYCLE模型黑：实测黑：实测红：模拟MODCYCLE模型MODCYCLE模型MODCYCLE模型MODCYCLE模型讲解结束讲解结束
