基于贝叶斯分析的水文组合预报模型.pdfVIP

第l6卷2期 应用基础与工程科学学报 Vo1．16，No．2 2008年4月 JOURNAL OF BASIC SCIENCE AND ENGINEERING April 2O08 文章编号：1005-0930(2008)02-0287-09 中图分类号：TVI24 文献标识码：A 基于贝叶斯分析的水文组合预报模型 张 弛， 李 伟， 周惠成 (大连理工大学土木水利学院，辽宁大连I16o24) 摘要：水文组合预报方法是对多种预报模型的预报结果进行组合分析的一种预 报方法．针对历史洪水数据不丰富的流域难于利用多种水文组合预报模型进行 洪水预报的现状，构建了以贝叶斯分析为基础，同时结合专家经验、马尔科夫蒙 特卡洛模拟、Gibbs抽样法，并引入实时校正的组合预报模型．以嫩江流域为实 例，对基于贝叶斯分析的组合预报模型的精度进行了验证．验证结果表明，该模 型可行而且实用，预报精度明显高于单个模型的预报精度． 关键词：组合预报；贝叶斯分析；马尔科夫蒙特卡洛模拟；Gibbs抽样 水文预报是防洪决策支持系统的核心，及时准确的洪水预报对于防洪决策至关重要． 在现阶段，没有一个水文预报模型是通用的，每个模型都有各自的优点与局限性，针对某 个特定流域如果仅依靠一个模型进行预报往往在特定条件下存在潜在的危险．因此，目前 许多研究集中在利用多种水文预报模型进行组合预报．一般的组合预报模型针对多个预 报模型，分析它们在不同条件下的适用情况，利用智能算法训练组合预报模型，这势必需 要更多的数据，而某些情况下，往往缺乏数据，针对这个问题，本文提出一种基于贝叶斯分 析的组合预报模型，这种预报模型需要较少的历史数据，利用贝叶斯模型、专家经验、马尔 科夫蒙特卡洛模拟等方法求出预报值．当水文水资源系统信息贫乏时，特别在发展中国 家，其有关水文和气象等的历史纪录较短，不可避免地面临着对不确定性和风险的认识和 处理问题，贝叶斯分析提供了一条可行而有效的途径…．贝叶斯定理起源于统计分析，是 在数据挖掘领域充分发展的算法，具体形式如下： L(~ ／O)fl (O) ㈩ 式中 (0)表示的先验密度， (x／O)表示给定0时样本信息 的似然函数， (O／x)表示 汲取样本后所校正的0的后验密度． 为0的论域．显然，贝叶斯分析将参数看作随机变 量来分析，这与水文预报模型中的参数的不确定性十分类似．目前有许多基于贝叶斯分析 的水文学方面的研究，主要集中在水文预报中不确定性的处理方面，代表方法是 Krzysztofwicz提出的贝叶斯概率水文预报理论 剖(Bayesian forecasting system，BFS)，其主 收稿日期；2006—10—10；修订日期：2008-01-09 基金项目：国家自然科学基金委员会、二滩水电开发有限公司雅砻江水电联合研究基金资助项目；大连 市科技计划项目(2007E21SFI65) 作者简介：张弛(1975一)，男，博士，讲师．E—mail：czhang@dlut．edu．cn 应用基础与工程科学学报 要思路是：根据其特征和对水文预报的影响大小，将总不确定度分为输入不确定度和水文 不确定度分别计算．针对降雨径流模型来说，输入不确定度主要指预见期内的降雨量，而 水文不确定度指的是模型误差、测量误差等．将两种不确定度分别计算后综合为总不确定 度积分后得到预报结果．而本文是针对各种模型预报结果，从组合预报的角度利用贝叶斯 理论． 1 模型的关键技术 1．1 贝叶斯分析 假定有n个模型进行水文预报，得到n个结果{ ， ，…， }，而 为水文真实值． 利用贝叶斯分析，问题就变为将n个确定水文预报模型的预报结果作为条件，预报在该条 件下水文值的问题，基本公式如下： P(Xo邶 …X2… = 式中，由于n个预报模型相互独立，则有 P({ ， ，…，X ~／Xo)=Hp(x ／Xo)