水文频率线型参数不确定性分析.pdfVIP

水文额率线型参数不确定性分析 水文频率线型参数不确定性分析木 桑燕芳王栋 (南京大学水科学系南京210093) 摘要水文频率计算过程中的线型参数估计结果存在不确定性，相应水文设计结果存在风险。为此，本文 ChainMonte 应用贝叶斯理论和AM—MCMC(AdaptiveMetropolis—MaAOVCarlo)采样方法，定量分析水文线型 参数不确定性问题。选用一实测水文序列，取一(c。)、二(C。、C．)、三(x、C。、C。)个参数三种情况分 别讨论P一砸线型中的参数不确定性。分析结果表明：水文线型参数不确定性是客观存在；应用贝叶斯方法可 以求得各参教的后验分布；通过对各参数后验分布的特征统计量和不同置信区间分别进行计算分析，可最终 实现对线型参数不确定性的定量描述。在参数不确定性分析的基础上，可提高相关水文设计结果的合理性和 可靠性。 关键词 水文时问序列 水文频率分析 贝叶斯方法 不确定性风险 1研究背景 由于水文循环过程中许多复杂因素的共同影响¨-3J，水文模型预报结果的不确定性是客观存在。影响水 o。水文 文预报量不确定性的因素主要分为输入不确定性和水文(包括模型结构、模型参数等)不确定性【4J 频率计算是进行“概率”预报，即求解水文变量发生大于等于某特定值时的概率，但“概率”结果同样存 在不确定性。针对水文频率计算问题，模型的输入(一般为某频率值)确定；我国结合实际情况一般选用P 一Ⅲ型线型。因此，研究水文线型参数的分析和求解方法成为重点和关键。 对于水文频率线型参数分析问题，参数优化求解方法已有多种，有研究对不同方法的优劣进行对比后得 到如下结论¨o】：POME与极大似然法效果相当，优于其他方法。但目前关于水文线型参数不确定性分析方 法的研究很少，因此充分挖掘已知水文序列中包含的信息并分析线型参数的不确定性，以提高水文频率分析 结果的可靠性，对实际生产应用具有重要的意义。 贝叶斯方法能利用已知信息(包括客观信息和经验、直觉和判断等主观知识)，最大程度地描述研究对 象的不确定性，并最终作出概率预报Hjl。引入贝叶斯方法，可以很好地研究水文频率计算中的线型参数不 确定性问题。在本文中，首先探讨引起水文线型参数不确定性的原因，然后提出基于贝叶斯方法的水文线型 参数不确定性分析方法，最后结合实测水文序列进行定量分析，以深入认识和理解水文线型参数不确定性 问题。 2基于贝叶斯理论的水文线型参数不确定性分析方法 2．1水文频率计算中的参数不确定性 影响水文频率计算结果精度和可靠性的因索主要有两方面：一是水文资料；二是水文频率计算方法．。相 比于较为成熟的水文频率计算方法，水文资料问题成为主要影响冈素。实际中，由于我国大部分站点实测的 水文资料长度有限，有时针对某些具体同题时难以获得有效的资料，样本序列特性不能很好地反映水文序列 的总体变化特性，代表性差；加之气候变化和人类活动对下垫面条件的强烈干扰，序列的一致性变差，序列 还原时也存在许多困难和问题。因此在分析实际水文问题的过程中，对水文序列进行代表性、一致性和可靠 性审查十分重要。 水文频率计算时由于掌握的水文资料有限，存在较大的抽样误差和不确定性，而参数的先验分布未知； 且不同特性水文序列的参数先验分布类型未必相同，因此欲准确描述线型参数不确定性十分困难。目前关于 · 基金项目：国家自然科学基金40730635)。 第一作者简介：桑燕芳(1983一)，男(汉)，山西长治人，南京大学水科学系硕士研究生，主要从事工程水文和随 机水文研究。E—mail：sunsangyf@sintom 第四篇 水文水资源斯方法新技术 不确定性。研究结果认为：仅掌搬短序列时抽样误差大；偏态性明显殿求解某稀遇频率对应的水文设计德时 不确定性凝大；提出通过延长序列长度、分析序列代表性、弓I入较为可靠的历史洪水、对成果作合理性分析 等手段霹竣夺这释苓确定缝。 但上述的各种手段并不能完全消除水文线捌参数不确定性，因此实际中需要碍求合理的方法