搜索
1021380135火炎焱燚Swatcup简单使用说明（以SUFI方法为例）1．数据准备：已经进行过至少一次模拟的swat输出文件，一般放于simxx目录下（xx是数字）swatcup需要的是其中的txtinout文件夹。2．Swatcup运行流程：选择new建立工程1021380135火炎焱燚选择txtinout文件夹只有选择正确才能使用选择校准时的模拟方法，具体见SWATCUP说明书1021380135火炎焱燚Finish之后工程就开始自动建立，工程建立完成后是一个“模子”，具体的东西还要我们自己来修改。命名以及选择输出文件夹，这里就随便了这里是设置校准参数的地方，我们需要做的是更改参数并为每个参数设定取值范围，需要注意的是，参数的修改方法有v__,a__,r__三种，分别是赋值，加值，乘以某值（百分比浮动）。具体可见说明书。每个参数后面的两个数字别是在下一次模拟中该参数取值的下界和上界，可以根据需要随便赋值，前提是不能超过理论上的最大范围参数个数以及模拟次数，参数个数要和下面的实际个数一样1021380135火炎焱燚File.cio是一些swatcup运行的基础参数，非专业人士请勿乱动哈（反正我是不会）。我们需要更改的有NBYR,IYR,IDAL（年中的天数，有时会出问题，在365和366之间改改试试）以及NYSKIP（缓冲年）开始模拟数及模拟次数，注意和par_inf对应1021380135火炎焱燚这里是swatcup接受修改的参数的上下界，参考用，一般禁止修改。在par_inf中的参数的范围需要注意不要越界。观测数据的使用，这里仅说一下rch的（其他的没用过）。本图中，1是指变量个数，如果要同时校准多个出口，这里可以更改数字，但此处写几下面就要写几组观测数据，格式一样。Flow_out_29是命名用的，29是模拟的出口所在subbasin。84是观测值个数，需要和你的模拟年对应好，一年12个月哦，365或366天，一定要算准，不然出错。1021380135火炎焱燚格式是：序号（空格）FLOW_OUT_月（日）_年（空格）观测值。建议在excel中做，很方便可以使用CONCATENATE函数。注意这里的月（日）是指这一年的第几个月（日）。如果不模拟sub以及hru的话，observation大类下的其他两个小类就不用管了，里面写什么也无所谓。后面会说为什么。此处就是输出的命名，注意和oberved_rch中一样就行了。同样注意，多个出口要写多个名字1021380135火炎焱燚此处每一行都有详细注释，因此很清晰，图中的32指的是模型中划分的subbasin的数量。需要注意的是如果是多出口校准，需要更改红框部分，1修改为出口个数，同时将29修改为出口subbasins的序号，例如，3出口可以写成：272829。另外注意模拟年份要对应，beginning年是初始年减去缓冲年之后的年份。下面的步长和你的模拟对应好就行了。关于下面那堆绿字，是为了校准不同目标的人们服务的，我们校准径流，所以是output.rch文件中的第7列，如果校准泥沙呢？那就是sed_out,在第11列。对应更改的就是variablecolumnnumber(s)intheswatoutputfile。Swatcup允许同时进行多目标校准，如果我同时想要校准flow_out和sed_out,那就如下图接下来，见下页：1021380135火炎焱燚这里是将观测数据再输入一遍，另外就是改些参数（一般我是默认的）第二行的Objectivefunctiontype,1=mult,2=sum,3=r2,4=chi2,5=NS,6=br2,7=ssqr根据个人需求更改。具体方法含义见说明书。这里之前写过类此的，同var_file_rch1021380135火炎焱燚其他的就不用管了然后是这里，前两个不用动，第三个和第四个如下这里就可以解释一下了没因为我们不使用sub，hru，以及各种带NO的校准过程，所以都点了取消掉，这样前面的那些对应的项目就不会被运行，里面写什么也都无所谓了。1021380135火炎焱燚另外就是不要忘记及时保存更改。然后就运行吧，少年们。直接all，你会发现直接Y，回车，然后一路回车，等着看结果。1021380135火炎焱燚这是结果，95ppu是模拟和实际图Bestpar是在你所给范围一及模拟次数内得到的最好结果对应的参数Bestsim是在以上参数下的模拟结果Newpar是swatcup的建议你下一次使用的参数取值范围。Summary_stat就是统计结果了，r2，ens都在接下来要做的就是1.如果你满意了，那么把bestsim带回到swat中去改参数，模拟，进行接下来的工作去吧。2.如果你不满意，那么根据bestpar和newpars参考着修改你的参数取值范围，或者再加减一些参数（记得改参数个数），个人喜好了，然后重复的calibration吧，直到结果满意了为止。大概就是这么多了，模拟泥沙，污染的步骤都差不多，是传说中的说明书所在。其实每个参数后面都有介绍，大家详细看看。另外注意的就是swat中的参数并不是所有的都1021380135火炎焱燚被swatcup支持，注意筛选。Cup的参数带回swat模拟结果可能有一些误差。另外就是如果模拟结果出错，一般95ppu是看不到的，这时候你就要检查名称，变量个数，模拟个数等等多处对应的地方是不是有些地方没对应。365和366也常常引发问题。大家多多注意。以上是大致步骤以及一些心得，错误纰漏欢迎指正。鸣谢“一叶舟”在我刚弄swatcup时的帮助，最后祝大家都能方便快捷干净利落的校准好。