[计算机]新安江模型参数有效优化及不确定性评估.pdfVIP

新安江模型参数有效优化及不确定性评估 王文川 ，程春 田 ，邱 林 ，杨斌斌 (1．大连理工大学水电与水信息研究所，辽宁大连 116085；2．华北水利水电学院，郑州 450011) [摘要] 应用新安江模型进行水文模拟时，由于模型本身的不足及参数多、信息量少等原因，会出现率定的 最优参数组不唯一、不稳定等问题。考虑到以往的参数优选 ，都只得出一个参数组，不能反映出其不确定性 状况。提出应用基于马尔可夫链蒙特卡罗(MCMC)理论的SCEM—UA算法，通过双牌流域以1h为时段间隔 的36场典型洪水数据对新安江模型参数进行优选和不确定性评估。结果表明，该算法能很好地推出新安江 模型参数的后验概率分布 ；率定和检验结果分析也表明，应用 SCEM—UA算法对新安江模型进行优选和不确 定评估是有效和可行的。 [关键词] 新安江模型；参数率定；不确定性评估 ；SCEM—UA [中图分类号] P338 ．9 [文献标识码] A [文章编号] 1009—1742(2010)03—0100—08 上建立起来的，加上率定模型的数据包含的测量误 1 前言 差，因此模型不可能真正模拟现实的系统。Kuczera 概念流域水文模型是 以流域为系统，模拟流域 和Parent指出 j：没有对参数不确定性的评估，无 上降雨径流形成过程，在进行水文规律研究和解决 法评定模型参数值和流域特征之间的相互关系及对 生产实际问题时起着重要作用，除要求其结构合理 将来水文事件响应的置信度。即使对模型应用有熟 外，模型模拟效果的好坏很大程度上依赖于模型参 练的技巧和想象力，也不应仅依靠唯一确定的模型 数的选择。由于水文现象具有高度的复杂性和随机 参数。 性，水文模型在使用过程中，参数多，信息量少，就会 近年来，水文模型参数的不确定性评估和水文 产生参数的相关性和不稳定性等问题。正是因为模 预报的不确定性分析成为国际上一个活跃的研究方 型各参数之问可能存在相依性，以及所构造的 目标 向。研究方法主要包括 GLUE (generalizedlikeli— 函数的非单峰(谷)性，按最优化方法求得的最优参 hooduncertaintyestimation)方法 ，BaRE(bayesian 数组可能不是唯一的，令参数优选成为十分复杂和 recursiveestimation)方法 ，马尔可夫链蒙特卡罗 困难的问题 。新安江模型作为一个通过长期实践和 (MCMC，markOVchainmontecarlo)方法… ，SCEM— 对水文规律认识基础上建立起来的具有分散性参数 UA(shuffledcomplexevolutionmetropolisalgorithm) 的概念性水文模型 ¨J，目前对其参数的优选取得了 方法 ¨’”等。相对来说，对我国应用得非常广泛的 许多可喜进展，主要有遗传算法 (GAs)L2,3j，SCE— 新安江模型，其模型参数不确定性分析方面研究还 UA l4j粒子群优化算法 。等。这些方法的目的都 较少见，还没有开展其参数不确定性范围的有效性 ， 是为了找到一个最好的参数组，无法对模型参数的 评估研究。 不确定性进行评估。事实上，模型参数通常会有误 近些年来，随着水情测报系统的广泛深入应用， 差倾向，主要由于模型是在假设、概化和判断的基础 特别是 1998年大洪水后国家防总进行全国水库洪 [收稿 Et期] 2008—09—17；修回日期 2008—12—03 [基金项 目] 国家 自然科学基金；水利部公益性行业科研专项 (200801015)；华北水利水电学院高层次人才科研启动资助项 目 (200821) [作者简介]