一寺ln·芷{掣口ER-人民黄河.pdf

第32卷第11期 人民黄河 V01．32，No．11 2010年11月 YELLOWRIVER Nov．．2010 【水文·泥沙】 基于Copula函数的洪峰洪量联合分布研究 侯芸芸，宋松柏 (西北农林科技大学水利与建筑工程学院，陕西杨凌712100) 摘要：采用Archimedean a的洪水资料为例进行了洪峰和洪量的联 Copula函数族中的4种函数，以陕北地区神木站48 合分布及两变量的联合重现期计算。结果表明：通过相依性度量，神木站洪峰和洪量有较高的相依性；Fank Copula函数 拟合洪峰和洪量的联合分布效果较好；运用Copula函数进行洪峰洪量联合分布的计算相对简单、灵活，可应用于多变量 水文频率分析计算。 关键词：Copula函数；洪水频率分析；联合分布；洪峰；洪量 中图分类号：P333 文献标识码：A doi：10．3969／j．issn．1000．1379．2010．11．017 洪水特性包括洪峰流量、洪水总量、峰现时间、洪水历时 (2)ClaytonCopula函数。 等。现行的洪水频率分析中，一是局限于采用单变量分布进行 C(u，口)=(／／,一’+／3一。一1)-1朋0∈(0，∞)(3) 洪水特征量频率分析，忽略了洪水总量、峰现时间、洪水历时、 洪峰流量等因素之间的相关关系．不能全面地描述洪水特 c(Ⅱ，口)=uv／[1—0(1—11,)(1一F)]0E[一l，1) 征一1；二是应用多元正态分布或者采用具有相同的边缘分布描 (4) 述多变量联合分布拉J。这两种方法都不能全面地描述洪水事 (4)FrankCopula函数。 件的内在规律和分析各个特征属性之间的相互关系。 R Copula函数。。是构建多变量联合分布的一种有效方法， 可以构造边缘分布为任意分布的多变量联合分布函数，描述变 (5) 量之间的相关性结构，具有极强的适应性和灵活性。在多变量 水文分析计算中，Copula函数理论主要应用于洪水频率、降水，，(Y)；一为Copula函数的参数。 频率、干旱特征分析、洪水或降水遭遇问题以及水文随机模拟 在洪水事件中，人们通常关注的是洪峰和洪量超过某一特 等方面旧一…。笔者以陕北神木站为例，应用Archimedean Copu· 定值的概率。因此，洪水频率计算主要是用来确定一定频率下 la函数构建边缘分布均服从P一Ⅲ型分布的最大洪峰流量以 的设计洪水值或者是一定洪水设计值下的设计标准，即重现期。设 及洪量之间的联合分布，进行洪水事件的联合分析，推求洪峰、 洪水事件中的两个特征变量为x和y，单变量洪水重现期为 洪量两变量的联合分布函数以及联合重现期，并且将两变茸联 r(茗)=1／[1一Fx(戈)] (6) 合分布结果与单变量的洪水频率分析结果进行比较，以期为陕 r(y)=1／[1一只(y)] (7) 北地区水利工程规划和防洪抗灾提供依据。 式中：r(髫)为单变量x的重现期；T(，，)为单变量y的重现期。 对于两变量重现期，一般有两变量联合重现期和两变量同 Copula函数和重现期