1、处理气象数据,留下2005-2009年气象数据, ArcSwat模拟2005-2009年液圆
2、SWAT-CUP填入2007-2009年的径流与产沙,率定参数。
3、率定出的参数带回ArcSwat中的数据库,增加气象数山迹据为2005-2012年,再次模拟
4、对比2010-2012年的模拟径流产沙与2010-2012年的实测值。
疑问1,2012年没有雨量站数据,是否全部舍弃2012年?
疑问2, 如果先上游后下游顺序率定?
关于验证,因为验证结果的评价就是R2和NS,而SWAT-CUP每次率定结果中都会有这两个值。所以在验证时,只需要扩展气象数据到2012年,然后按照2007-2009率定出来的新参数再次率定,根据得到的R2和NS判断率定的效果。
关于率定:一次模拟500次,太慢,可以先模拟10次,将推荐的参数带回再模拟10次,如此往复。ps.经过测试,此方法不靠谱,推荐的参数范围太大,不准确。
关于敏感性:可以用SWAT自带的敏感性分析工具,也可以用SWAT-CUP的敏感性分析工具。这里有一处不明,如何选择参数的修改方式,即赋值、乘以某值、加上某值,三选一。ps.参照杨军军硕士毕业论文中的修改方式吧。
查一下研究产沙的文章其中用的指标是什么,单位是t还是kg/m3或者其闹唯塌他的?大部分泥沙的文章都是用sed_out指标,单位为t。但是我拿到手的泥沙数据只有含沙量kg/m3,可以用.rch中的sedconc指标(单位时长内单位体积的产沙量)替代率定(没找到确切证据可以这样做)。研究SWAT模拟的输出文件,flowout*sedconc=sedout(考虑单位转化的话,flowout(m3/s)*1day*sedconc(mg/L)=86.4sedout(t)),发现部分数据符合此规律,部分不符合,不明何故。
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)