P_对数P_和GL分布线型比较研究_刘洁.pdfVIP

收稿日期： 2012－01－29 基金项目：国家重点基础研究发展规划 (973 计划 ) （ 2010CB428405 ）； 国家自然科学基金重大项目课题（） 作者简介：刘洁 （ 1983- ），女，江西景德镇人，讲师，主要从事城市水文与市政工程研究。 E-mail:lj2000liujie@ 通讯作者：陈华 （ 1977- ），女，福建建瓯人，副教授，主要从事水文水资源研究。 E-mail:chua@ 水 文 JOURNAL OF CHINA HYDROLOGY 第 33 卷第 3 期 2013 年 6 月 Vol．33 No．3 Jun．, 2013 1 引言 目前国外设计洪水的计算主要通过两种途径，一 是依据实测流量资料系列直接计算设计洪水 （或称数 理统计法），二是根据暴雨资料间接推算设计洪水（或 称水文气象法）。 我国设计洪水的计算和分析，主要采 用的是第一种途径，从历年实测流量资料中极值选样， 计算其经验频率后，用 P-Ⅲ型适线的方法最后确定设 计洪水参数值。由流量资料系列计算设计洪水，国内已 经在洪水频率及样本参数的估计、 洪水频率线型的选 择、 不连续系列设计洪水计算方法等多方面进行了广 泛的研究 [1-6]，其目的都是为了寻求更佳的概率模型来 拟合样本资料， 以使估计的样本参数最能反映总体的 内在分布特性，使计算的设计洪水更加合理可靠。目前 还无法得到哪一种洪水频率线型最适合国内设计洪水 计算的研究结论。 在样本参数估计上，已有很多学者 认为线性矩法是目前估计样本参数稳健性较好的方 法[7]。 由流量推求设计洪水，目前美国主要采用的是对 数 P-III ， 英 国 主 要 采 用 广 义 逻 辑 分 布 （ Generalised Logistic Distribution ，简记为 GL ） [1] ，本文在收集淮河和 珠江流域长系列水文资料的基础上， 用线性矩法估计 参数， 对中英美三国目前采用的洪水频率计算线型进 行比较分析，为洪水频率计算提供参考。 2 洪水频率理论分布线型 2.1 P-Ⅲ分布曲线 中国《水利水电工程设计洪水规范》规定用 P-Ⅲ 型分布曲线作为水文随机变量的概率分布线型， 其概 率密度函数为[1] f （ x ） = β α Γ （ α ） e -β （ x-a 0 ） （ x-a 0 ） α-1 （ 1 ） 其中： Γ （ α ） = ∞ 0 乙 x α-1 e -x dx a 0 =x 1- 2C v C s 乙 乙 α= 4 C 2 s β= 2 xC v C s 2.2 GL 分布曲线 英国 1999 年的《洪水估算手册》 [8]推荐采用广义 逻辑分布（ GL ）作为频率曲线分布线型， GL 分布是一 个三参数分布，表达式如下： P-Ⅲ、对数 P-Ⅲ和 GL 分布线型比较研究 刘 洁 1 ，陈 华 2 ，王金星 3 ，陈 星 2 ，李天元 2 （ 1. 成都大学城乡建设学院，四川 成都 610106 ； 2. 武汉大学水资源与水电工程科学国家重点实验室， 湖北 武汉 430072 ； 3. 水利部水文情报预报中心，北京 100053 ） 摘 要：以淮河和珠江流域为研究对象，收集了淮河 18 个和珠江 16 个水文站 1951~2010 年长系列日 流量资料。 通过年极大值取样，应用线性矩法估计参数，比较分析了 P-Ⅲ、对数 P-Ⅲ和 GL分布曲线在 两个流域洪水频率分析计算情况。 研究结果表明，对于小概率事件，在这两个流域对数 P-Ⅲ和 GL分布 曲线设计洪水值要大于 P-Ⅲ曲线分布计算的值，而通过误差分析，P-Ⅲ曲线的误差整体上小于对数 P- Ⅲ曲线。 关键词：淮河；珠江； P-Ⅲ；对数 P-Ⅲ；GL 中图分类号： P333 文献标识码： A 文章编号： 1000-0852(2013)03-0001-04 } （ 2 ） 第 33 卷水 文 表 1 流域及测站基本信息表 Table1 The information of the basin and hydrometry stations Q （ F ） =ξ+ α k 1- 1-F F - - k - - k≠0 （ 3 ） 式 中： ξ 为 位 置 参 数； α 为 尺 度 参 数； k 为 形 状 参 数； F 为不超过概率。 由重现期与不超过概率的关系式 T= 1 1-F ，又可将上式改写成如下形式： Q T =ξ+ α k 1- T- - - 1 - - - k （ k≠0 ） （ 4 ） GL 分布的可能值域范围为 -∞Qξ+ α k k0 ξ+ α k 燮Q∞ k0 （ 5 ） 因此，当 k0 时， GL 有上限约束。 当 k0 时，有下限 约束。 2.3 对数 P-Ⅲ分布曲线 美国联邦政府建