多变量情况下的水文频率分析方法及其应用.pdfVIP

水 利 学 报 2009年1月 SHUILIXUEBAO 第40卷第1期 文章编号：0559—9350(2009)01-0033·05 多变量情况下的水文频率分析方法及其应用 冯 平1，毛慧慧1，王 勇2 (1．天津大学建筑工程学院，天津300072；2．天津市水文水资源勘测管理中心，天津300061) 摘要：本文探讨了以Gumbel分布为边际分布的二维联合分布模型，给出了其经验频率和理论频率的计算方法，建立 并探讨了联合重现期和条件重现期分布模型的实用功能。由此，可以计算不同程度降雨和入境水遭遇组合的频 率，也可以在一定的降雨量条件下计算不同年入境水的发生频率，或特定的年入境水条件下计算不同降雨量的发 生频率，还可以给出某一重现期对应的某一安全指标下的风险区间，方便于解决与风险有关的多因素影响下的水 文计算问题。以天津市北三河系的降雨和入境水为例，分析了该模型的实用性，结果表明，采用Clayton作为联结函 数来描述水文随机变量的联合分布是有效的。 关键词：Gumbel分布；二维联合分布；频率；重现期 中图分类号：P333 文献标识码：A 人们逐渐意识到复杂的水文事件，例如洪水和暴雨，是一种多变量随机事件，并且其中的随机变量 是相互关联的。Cunnaneuo认为传统的单变量分析方法只能揭示有限的流域水文特性，这样，许多复杂 的水文问题要求对流域水文特性有更加深入和全面的了解。如果想真正掌握水文事件的特性，就必须 了解构成这些水文事件的两个或两个以上相关随机变量的联合概率分布特性，有关学者已经在这方面 做了一定的工作。国内学者郑红星等旧1和闫宝伟等¨1分别对不同水文区域降水的丰枯遭遇进行了分 析。冯平等H1认为暴雨或径流等水文现象总是表现为多变量事件，并采用风险分析的方法探讨了多因 素条件下的海河流域的防洪计算问题。国外学者S．Yue等b’61采用二维Gumbel模型对洪峰和洪量、洪 量和历时进行了分析，采用Gumbel 用Copula函数探讨了双变量降雨频率分布问题。综合看来，现阶段的研究多数是将联合分布用于同一 场次的洪水、降雨等相关因素的分析中，如洪峰和洪量，洪量和历时等。此外，实际研究表明许多水文极 位¨引、土壤含水量【11。、降雨极值u21等都拟合较好。因此，本文采用Gumbel分布作为降雨量和入境水量 的边际分布，并尝试将二维Gumbel模型、Gumbel 水间的频率分析，重点对重现期及其与风险区间的相关关系做了探讨。 1 Gumbel模型 F(戈)和F(Y)分别是随机变量X和Y的边际分布函数，其表达式为： 眦≤善，=唧--exp(一等)】；哪≤，，，一p[- 唧(一Y_-，ur)】 其中，Gumbel分布的参数口和址可由矩法求得，即： 收稿日期：2008．04．28 基金项目：国家科技支撑计划(2006BABl4805)；天津市科技发展计划项目(OYFSYSF602400) 一33— 口：曼s： “：肘一o．577口 (2) ，r 式中：吖和5是样本的均值和均方差。 2二维联合分布模型 2．1二维Gumbel模型和Gumbel Logistic模型 Gumbel 式分别为 F(聋，Y)=F(菇)F(y)exp{一，。【1／lnF(z)+1／lnF(Y)】-1)(0≤rI≤1)(3) (4) 在式(3)、式(4)中，r。、r：都为描述随机变量工和y的相关系数，其表达式分别为‘13’141： rl=2