第九讲 SWAT模型参数率定和验证.ppt

相关网站、资料 SWAT官网： / 论坛： SWAT软件： /software/ AVSWAT、ARCSWAT、SWAT-CUP、 AVSWATX 下载地址： /pub/swat/pc/swatav/avswatx.zip 或者：/swat/downloads/pc/swatav/avswatx.zip SWAT-cup下载： http://www.eawag.ch/forschung/siam/software/swat/index * * 参数敏感性分析 对于径流：CN2、sol_awc、ESCO最敏感 对于泥沙：CN2、SPCON、BIOMIX最敏感 对于氮磷等污染物： CN2、 ESCO、Alpha_Bf、surlag等参数比较敏感 Out 1 Out 2 Out 3 Out 4 Out 5 Out 6 Out 7 Out 8 Out 9 参数 6 11 3 3 4 4 3 4 5 ALPHA_BF 3 5 6 6 9 5 6 5 6 ESCO 1 1 2 1 1 2 2 1 1 CN2 2 4 5 5 2 3 5 6 3 SOL_AWC 16 6 7 7 7 6 7 7 7 surlag 4 9 10 11 6 10 10 10 12 Sol_z 17 3 4 4 23 8 4 3 4 BIOMIX 13 7 1 9 5 15 1 2 10 SOL_ORGN 自动校核 单击Tools菜单下的auto-calibration and uncertainty，弹出下图所示对话框 自动校核 根据敏感性分析结果选择参数 选择子流域14进行自动校核 自动校核 操作过程 自动校核 查看结果 不确定性分析 自动校核后出现校核结果 可进行不确定性分析 关键步骤 观测数据文件的制作 注意参数界限 ArcSWAT2005中有手动调参的工具 SWAT-CUP参数率定 SWAT-CUP (SWAT Calibration Uncertainty Procedures) 功能： 敏感性分析 校核 验证 模型不确定性分析 SWAT-CUP参数率定 特性： 并行处理 输出结果的可视化 通过参数提取实现多目标函数公式化 参数95ppu提取和可视化 无观测数据时： 温度数据和降水数据在不确定性中考虑 SUFI2中目标函数极值可获得 SWAT-CUP运行过程 安装过程 SWAT-CUP运行过程 建立文档 SWAT-CUP运行过程 选择一次模拟结果 选择校准模拟方法 SWAT-CUP运行过程 校准时的模拟方法： SUFI2（ Sequential Uncertainty Fitting version 2 ） PSO GLUE(Generalized Likelihood Uncertainty Estimation ParaSol MCMC procedures 选择一种运行 SWAT-CUP运行过程 参数个数以及模拟次数，参数 个数要和下面的实际个数一样 设置校准参数： 更改参数并为每个参数设定取值范围 参数的修改方法有v__, a__, r__三种，分别是赋值，加值，乘以某值（百分比浮动） 每个参数后面的两个数字别是在下一次模拟中该参数取值的下界和上界，可以根据需要随便赋值，前提是不能超过理论上的最大范围 SWAT-CUP运行过程 开始模拟数及模拟次数，注意和par_inf对应 SWAT-CUP运行过程 File.cio文件为swatcup运行的基础参数，可不修改 SWAT-CUP运行过程 Absolute_swat_values:接受修改的参数的上下界，参考，一般不修改 SWAT-CUP运行过程 观测文件： 1是指变量个数，如果要同时校准多个出口，这里可以更改数字，但此处写几下面就要写几组观测数据，格式一样 Flow_out_29是命名用的，29是模拟的出口所在subbasin 84是观测值个数，需要和模拟年对应好，一年365或366天 SWAT-CUP运行过程 观测数据格式： 序号（空格）FLOW_OUT_月（日）_年（空格）观测值 多个出口制作多个文件 参看说明 1修改为出口个数，同时将29修改为出口subbasins的序号，例如，3出口可以写成：27 28 29 SWAT-CUP运行过程 将观测数据再输入一遍，修改部分参数：第二行的Objective function type, 1=mult,2=sum,3=r2,4=chi2,5=NS,6=br2,7=ssqr SWAT-CUP运行过程 执行文件设置 SWAT-CUP运行过程 前两个不用动，第三个和第四个如下 Calibration按钮 SWAT-CUP运行过程 执行 Y，回车 确定 率定结果 查看结果 95ppu是模拟和实际图 率定结果 参数点状图分析 率定结果