新安江模型参数全局优化研究-中国科技论文在线.PDFVIP

第 32 卷第 4 期 河 海 大 学 学 报 ( 自 然 科 学 版 ) Vol. 32 No. 4 2004 年 7 月 ( ) Journal of Hohai University Natural Sciences Jul. 2004 新安江模型参数全局优化研究 李致家 ,周　轶 ,哈布·哈其 (河海大学水资源环境学院 ,江苏 南京 　210098) 摘要 :采用单纯多边形进化算法对 3 个气候与流域条件各不相同的流域的新安江模型计算参数优 化问题进行了研究. 结果表明 :对于没有误差的水文资料 ,采用单纯多边形进化算法 ,可以优化出新 安江模型的全部参数并使参数收敛到真值 ;对于实测的水文资料 , 由于模型参数之间的共线性, 固 定参数 B 与 EX ,可同时率定模型其他参数 ,随着实测资料系列的增加 ,率定参数值将趋于稳定. 关键词 :新安江模型参数 ;单纯多边形进化算法 ;全局优化 ;共线性 ( ) 中图分类号:P338 　　文献标识码 :A 　　文章编号 2004 040376 04 新安江模型是一个半分散性概念性水文模型. 该模型的参数率定方法有人工试错法和传统的数学优化 方法. 人工试错法的缺点是主观性较强且缺乏客观标准 , 而传统的数学优化方法 , 如 Rosenbrock 方法和 Simplex 方法等 ,大多为局部优化方法. 因此 ,如何寻求概念性水文模型参数率定的全局最优解 ,成为 20 世纪 90 年代以来水文模型研究的热点问题. 美国亚利桑那大学于 20 世纪 90 年代开始研究萨克拉门特模型参数优化问题 ,后来提出了单纯多边形 ( ) [1 ] ( ) 进化算法 SCEUA 并得到了全局最优解 . SCEUA 结合了现有算法 包括基因算法等 中的一些优点 ,可以 解决高维参数的全局优化问题 ,且不需要显式的目标函数或目标函数的偏导数. 因此 ,本文采用单纯多边形 进化算法对 3 个气候与流域条件各不相同的流域的新安江模型计算参数优化问题进行了研究. 1 　新安江模型 ( 应用新安江模型进行水文模拟时 ,首先需按下垫面和地形等情况对流域进行分块 称为分单元或子流 ) 域 ,然后将每个子流域的出口流量演算至整个流域的出口断面. 这样单元流域的新安江模型有 14 个参数. 如果假定各个单元流域的模型参数一致 ,河道汇流采用马斯京干法 ,每个子流域的汇流河段数可以通过河道 汇流的长度来确定 ,那么又增加了汇流方面的 2 个参数. 因此对一个流域来说 ,需要率定的参数有16 个. 2 　SCEUA SCEUA 本身有多个参数 ,这些参数可分为不确定性参数和确定性参数 2 个部分, 只有选择好算法中的 参数, 才能使算法达到全局最优. 实际上, 参数最优值取决于问题的类型, 但通常可用 m = 2 n + 1 , q = n + 1 , y ( ) q 为每个子复合体的顶点数, y = 1 , z = 2 n + 1 作为缺省值. 其中, m 为每个复合体 complex 多边形的顶点数, 为每一个子复合体连续繁殖后代的个数, z 为每个子复合体进化的步数, n 为水文模型中需要优化的参数个 数. 理论上 m 可取比 1 大的任何