三水源新安江模型异参同效现象的探究.pdfVIP

·151· 第二部分 国际水文科学协会(IAHS)十年科学计划(2003--2012年) 三水源新安江模型异参同效现象的研究’ 熊立华郭生练 (武汉大学水资源与水电工程科学国家重点实验室武汉430072) 摘要在确定流域水文模型参数最优值时，总是会遇到“异参同效(equifinal— ity)”现象，即我们总是可以搜索到几组或者很多组不同的参数值，使得模型的 效率系数达到一样或者几乎一样的水平。产生这种现象的原因的有：流域水文过 程太复杂，而流域水文模型只是对自然现象的简化描述；水文模型结构本身的缺 陷；模型参数的冗余或相关性太强等。异参同效现象所导致的直接后果是我们最 终选定的一组“最优”参数值具有很大的不确定性。为了能够充分地掌握和评估 水文模型参数率定中的这种不确定性，我们就必须比较全面的分析异参同效现 ChainMonte 象。本文将采用基于马尔科夫链蒙特卡洛(Makov Carlo，简称 MCMC)的取样方法来研究三水源新安江模型中的“异参同效”现象。结果发 现，三水源新安江模型存在非常强烈的异参同效现象，给我们选择所谓的最优参 数值带来了巨大的不确定性。 关键词参数优化 新安江模型 异参同效 不确定性 水文模型研究中一个关键问题是模型参数的优化。只有在确定模型参数的最优值之 后，模型才能够用来进行水文预报或其他用途。当模型的参数比较多时，已经无法采用人 工试错法进行模型参数优化，而必须选用一些高效的自动优选法。自动优选方法有很多， 比如单纯形法、遗传算法、SCE—UA[1J等。在采用这些方法来确定流域水文模型参数最 索到几组或者很多组不同的参数值，使得模型的目标函数(比如效率系数)达到一样或者 几乎一样的水平。产生这种现象的原因的有：流域水文过程太复杂，而流域水文模型只是 对自然现象的简化描述；水文模型结构本身的缺陷；模型参数的冗余或相关性太强等。异 参同效现象的存在使得我们在最终选定一组“最优”参数值时具有很大的不确定性；因为 这组“最优”参数值往往与我们主观设置的优化过程控制条件有关，比如循环次数、收敛 误差范围等。为了能够充分地掌握和评估水文模型参数率定中的这种不确定性，我们就必 须比较全面地分析模型参数优化中的异参同效现象。英国著名水文学家Beven于1992年 Likelihood 提出基于贝叶斯方法的GLUE(GeneralizedUncertainty * (和国家自然科学基金项目。 作者简介：熊立华(1972～)，男，汉族，湖北荆门人，博士，武汉大学副教授。 ·152· 水问题的复杂性与不确定性研究与进展 分析和评价模型参数优化中的异参同效现象。这一方法目前已经在国内外得到了广泛的应 ChainMonte 用[3～5J。由于马尔科夫链蒙特卡洛(Makov 样方法能够较好地模拟参数的后验概率分布，本文采用MCMC方法来研究三水源新安江 模型中的异参同效现象。 1三水源新安江模型及目标函数 三水源新安江模型是我国应用最为广泛的水文模型。三水源是指地面径流、壤中流、 地下径流。研究中对壤中流和地下径流分别采用不同的线性调蓄水库模拟其汇流过程；对 地面径流采用Nash单位线法。可以进行优化选择的参数有15个。这些参数包括：流域平 均蓄水容量WM；流域上层平均蓄水容量WUM与WM之比pu；流域下层平均蓄水容量 WLM与WM之比pl；蒸发能力转换系数KE；流域蓄水容量曲线指数B；流域平均自由 水蓄水容量SM；流域自由水蓄水容量曲线指数EX，壤中流出流系数CI；地下水出流系 数CG；不透水面积系数IMP；深层蒸发系数CE；壤中流退水系数CKI；地下水退水系 数CKG；地表汇流单位线的线性水库个数NN；地表汇流单位线的线性水库退水常数 TT。 水文模型中常用的目标函数是Nash模型效率系数R2，模型参数优化目标可表