长江科学院院报 ›› 2015, Vol. 32 ›› Issue (6): 116-119.DOI: 10.3969/j.issn.1001-5485.2015.06.020

• 水循环演变 • 上一篇    下一篇

基于极大似然法的负偏水文序列参数估计方法研究

胡诗松a,陈进a,b,尹正杰a,c   

  1. 长江科学院 a.流域水资源与生态环境科学湖北省重点实验室;b.院长办公室;c.水资源综合利用研究所,武汉 430010
  • 收稿日期:2015-04-07 出版日期:2015-06-20 发布日期:2015-06-04
  • 作者简介:胡诗松(1990-),男,河南信阳人,助理工程师,硕士研究生,主要从事水文水资源研究工作,(电话)027-82828101(电子信箱)hushisong@mail.crsri.cn。
  • 基金资助:
    水利部公益性行业科研专项(201401020)

Maximum Likelihood Estimation of Negative-skewness Hydrological Series

HU Shi-song1,CHEN Jin1,2,YIN Zheng-jie1,3   

  1. 1.Hubei Provincial Key Laboratory of Basin Water Resource and Eco-environmental Sciences, Yangtze River Scientific Research Institute, Wuhan 430010, China;
    2.Administration Office, Yangtze River Scientific Research Institute, Wuhan 430010, China;
    3.Water Resources Department, Yangtze River Scientific Research Institute, Wuhan 430010, China
  • Received:2015-04-07 Online:2015-06-20 Published:2015-06-04

摘要: 针对皮尔逊Ⅲ型曲线参数估计方法中尺度参数β>0,无法适用于潮位、水位等负偏水文序列(β<0)参数估计的问题,从正偏序列与负偏序列的关系出发,提出了一种极大似然估计的新方法。该方法是将负偏序列先后转换为正偏序列、伽马分布序列,并利用负偏序列确定ξ初值,通过伽马分布序列迭代计算得到α,β的极大似然估计值。以某潮位站44 a的最高潮位数据为算例,比较了负偏极大似然法、负偏矩法、正偏查表法3种估计方法对负偏序列拟合的精度,结果发现负偏极大似然估计法的3项拟合指标:PPCC检验值为0.993,OLS检验值为0.010,KS检验值为0.052,均优于负偏矩法、正偏查表法。为此,负偏极大似然法对负偏序列拟合的精度为最优。

关键词: 水文频率, 负偏序列, 极大似然估计, 皮尔逊Ⅲ型曲线, MATLAB软件

Abstract: The commonly used probability density function of Pearson-III distribution is invalid for negative-skewness water level and tide level series (β<0). To solve this problem, a means of maximum likelihood estimation (MLE) based on the relationship between negative-skewness series and positive-skewness series is proposed. In this method negative-skewness Pearson-Ⅲ distribution is translated into positive-skewness Pearson-Ⅲ and gamma distribution in sequence. The initial value of ξ is decided by the negative-skewness distribution and the MLE of α and β are calculated through the gamma distribution. In the iterative process the sampling error with evaluating coefficient of skewness will be avoided. Moreover, the 44-year record of maximum tide level at a tide station is used as calculation example. The fitting accuracy of negative-skewness MLE, negative-skewness moment method and positive-skewness table look-up are compared. The fitting results of negative-skewness MLE (test value of PPCC is 0.993, OLS 0.010, and KS 0.052) are superior to those of the other two methods. In conclusion, negative-skewness MLE has the optimum fitting accuracy for negative-skewness hydrological series.

Key words: hydrological frequency, negative-skewness series, MLE, Pearson-Ⅲ frequency curve, MATLAB software

中图分类号: