基于极大似然法的负偏水文序列参数估计方法研究-长江科学院.PDFVIP

　 第３２卷第６期 长　江　科　学　院　院　报 Ｖｏｌ．３２　Ｎｏ．６　　 　　２０１５年 ６月 ＪｏｕｒｎａｌｏｆＹａｎｇｔｚｅＲｉｖｅｒＳｃｉｅｎｔｉｆｉｃＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅ Ｊｕｎ．２０ １５　 　 　　ＤＯＩ：１０．３９６９／ｊ．ｉｓｓｎ．１００１－５４８５．２０１５．０６．０２０ ２０１５，３２（０６）：１１６－１１９，１２６ 基于极大似然法的负偏水文序列 参数估计方法研究 ａ ａ，ｂ ａ，ｃ 胡诗松 ，陈　进 ，尹正杰 （长江科学院ａ．流域水资源与生态环境科学湖北省重点实验室；ｂ．院长办公室； ｃ．水资源综合利用研究所，武汉　４３００１０） 摘　要：针对皮尔逊 型曲线参数估计方法中尺度参数 ＞０，无法适用于潮位、水位等负偏水文序列（ ＜０）参数 Ⅲ β β 估计的问题，从正偏序列与负偏序列的关系出发，提出了一种极大似然估计的新方法。该方法是将负偏序列先后 转换为正偏序列、伽马分布序列，并利用负偏序列确定 初值，通过伽马分布序列迭代计算得到 ，的极大似然估 ξ αβ 计值。以某潮位站４４ａ的最高潮位数据为算例，比较了负偏极大似然法、负偏矩法、正偏查表法３种估计方法对负 偏序列拟合的精度，结果发现负偏极大似然估计法的３项拟合指标：ＰＰＣＣ检验值为０．９９３，ＯＬＳ检验值为０．０１０，ＫＳ 检验值为０．０５２，均优于负偏矩法、正偏查表法。为此，负偏极大似然法对负偏序列拟合的精度为最优。 关键词：水文频率；负偏序列；极大似然估计；皮尔逊 型曲线；ＭＡＴＬＡＢ软件 Ⅲ 中图分类号：ＴＶ１２　　　文献标志码：Ａ　　　文章编号：１００１－５４８５（２０１５）０６－０１１６－０４ 线参数估计的极大似然估计方法（后文统一简称为 负偏极大似然法），并以某潮位站４４ａ的最高潮位 １　研究背景 数据为算例，比较负偏极大似然法、负偏矩法、正偏 在水文计算中常需推求 １００ａ一遇、甚至 查表３种方法对负偏序列拟合的精度，同时，对本文 １０００ａ一遇的设计值，水文频率曲线是合理解决这 提出的方法加以了分析和验证。 一问题的有效“工具”。根据《水利水电工程设计洪 水计算规范》（ＳＬ４４—２００６），我国水文频率计算线 型采用皮尔逊 型曲线。目前常用的皮尔逊 型曲 ２　计算原理 Ⅲ Ⅲ 线尺度参数 ＞０，适用于 Ｃ ＞０的正偏水文序 β 负偏水文序列的皮尔逊 型曲线计算主要有 ｓ Ⅲ ［１］ 列 ，但分析水库水位、潮位等水文序列时经常出 ２种思路： 寻找正偏负偏序列的关系，把负偏序列 ① 现负偏的情况，此时如仍沿用上述曲线配线，解得的 计算转化为正偏，然后查表或者用矩法等方法估计 理论频率曲线和经验点偏离较大。为解决该问题， 参数； 直接根据负偏皮尔逊 型曲线密度函数求 ② Ⅲ