基于运动扩散模型的含零值降水序列频率计算.pdf

水 利 学 报 2018年3月 SHUILI XUEBAO 第49卷 第3期 文章编号：0559-9350（2018）03-0387-09 基于运动扩散模型的含零值降水序列频率计算 李丹丹，宋松柏，李运平 （西北农林科技大学 水利与建筑工程学院，陕西 杨凌 712100） 摘要：应用运动扩散（KD）模型，研究含零值降水序列频率计算。以陕西省6个测站2月份降水序列为例，在 Strupczewski推导矩法、极大似然法参数计算公式的基础上，利用数学变换原理和数值计算原理，推导了KD模型 含零值序列分布概率权重矩法参数计算公式。采用AIC准则、OLS准则和残差平方和（RSS）最小准则，进行KD模 型、频率比例法和II型乘法分布模型拟合效果评价。结果表明：KD模型可以应用于含零值降水序列频率计算， 且该模型拟合效果优于频率比例法和II型乘法分布法；3种KD模型参数估计方法中，概率权重矩法估计KD模型 参数拟合度最优。文中方法以期为我国含零值水文序列频率计算提供一种新的计算方法。 关键词：含零值降水序列；频率分析；KD模型；概率权重矩法；频率比例法；II型乘法分布 中图分类号：P333.9 文献标志码：A doi ：10.13243/ki.slxb 1 研究背景 水文频率计算是根据水文现象的统计特性，利用现有水文资料，分析水文变量设计值与出现频 率（或重现期）之间的定量关系的工作过程，对于水利工程规划、设计、管理以及水资源配置等具有 ［1］ ［2］ 至关重要的作用 。降水设计值是进行用水管理的主要依据 。通常情况下，水文频率计算采用的水 文序列不含有零值，但在干旱和半干旱地区的水文频率分析中，经常遇到包含若干个零值的水文序 ［3-6］ 列，而零是变异范围的下限 。这些水文序列可能是干旱地区的月降水量，水库设计中的防洪库 ［7］ ［8］ 容系列以及河道断流等 。因此，常规的水文频率计算方法无法进行这类水文序列频率计算 ，而 ［9］ ［4］ ［9］ ［5］ 且以往研究中对于这类水文序列的频率分析关注较少 。Jennings等 、Woo等 、Wang 等 提出 了增加小水文值、忽略零值、全概率法以及应用次序统计量构造全序列分布等方法计算含零值水 文序列的频率。1980年代以来，我国一些学者也进行了含零值水文序列的计算方法研究，如间接 法、中值适线法和Ⅱ型乘法分布等方法［7，11-14］。《水利水电工程设计洪水计算手册》推荐使用频率比 ［11，15-17］ ［18］ 例法和 Ⅱ型乘法分布法 。张良等 利用这两种方法对安固里河洪峰流量进行了分析。2003 ［8］ 年，Strupczewski等 提出了基于运动扩散模型（KinematicDiffusionModel ，简称KD模型）的含零值序 列频率计算模型，这是一种新型的含零值水文序列频率计算方法，该模型已成功地应用于美国干旱 地区洪水频率计算，取得了较好的拟合效果。 目前，Strupczewski等推导了KD模型矩法、极大似然法参数计算公式。众所周知，概率权重矩 法是水文频率分布参数计算精度较高的方法［19-21］。但是，KD模型目前没有概率权重矩法参数计算公 式，且这类模型尚未用于我国含零值水文序列频率计算。本文从概率权重矩定义出发，应用数学变 换原理和数值计算原理，严格地推导了概率权重矩法推求KD模型参数计算公式，给出了相应的数值 计算方法，并与矩法、极大似然法估计KD模型参数进行比较，为KD模型参数估计方法选择提供依 收稿日期：2017-09-18；网络出版日期：2018-02-26 网络