梁希金，王丽萍，张验科^*，张佳新，叶凯华

华北电力大学可再生能源学院，北京

收稿日期：2019年12月7日；录用日期：2019年12月20日；发布日期：2019年12月27日

当某一条河流有多条支流汇入时，干流洪水与支流洪水的遭遇是威胁下游防洪控制点安全的主要因素，干支流洪水遭遇风险估计对于下游防洪安全具有重要意义。基于Copula函数易于将多个边缘分布构造成联合分布的优点，将其与频率组合法相结合，推导出了n维组合洪水分布函数表达式，并建立了相应的洪水遭遇风险估计模型。以雅砻江、鲜水河和庆大河洪水遭遇为例，采用双参数三维Copula函数对下游的洪水遭遇风险进行了估计，得到了下游组合洪水分布函数及防洪风险率。结果表明，该模型可以在一定程度上量化分析干支流洪水遭遇的随机性，为下游的防洪安全提供了一定的参考。

关键词 :Copula函数，频率组合法，组合洪水，防洪风险率

Copyright © 2020 by author(s) and Wuhan University.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

随着我国各大流域联合防洪调度工作的逐渐深入，流域性洪水引发的多条河流的洪水遭遇分析工作越来越迫切。洪水遭遇风险研究对下游的防洪安全具有重要意义。传统的多变量水文分析计算方法主要包括多元正态分布 [ 1 ] [ 2 ] 、特定边缘分布构成的联合分布 [ 3 ] [ 4 ] 、非参数方法 [ 5 ] 、将多维联合分布转换成一维分布的方法 [ 6 ] 等。传统方法往往存在着一些缺陷，比如多元正态分布在数据转换时容易造成信息失真，非参数方法构造的联合分布的边缘分布类型未知等 [ 7 ]。目前，国内外针对上游有支流汇入的防洪控制点的风险估算，多是针对一条支流汇入的情况，对于多条河流汇入的情况研究较少；并且多是将干流洪水的随机性作为研究对象，区间洪水设为特定值进行计算，这样显然与实际有所偏差。因此，研究新的方法，充分考虑干支流洪水的随机性，对于上游有多条支流汇入的防洪控制点的防洪安全具有重要意义。

频率组合法是研究洪水遭遇的传统方法，目前多是针对二维水文随机变量进行实例研究，对于多维水文随机变量的研究相对较少，当涉及的维数较多时，用传统的解析方法来确定联合概率密度函数的解析式以及对其积分求解是非常困难的，因此，通常采用离散叠加法进行数值求解。但是，干支流洪水变量间往往存在相关性，在采用离散叠加法进行求解时，需要对变量进行独立性处理，此过程难免会出现数据失真。

Copula函数能够通过变量的边缘分布和相关性结构两部分来构造多维联合分布，形式灵活多样，是构建联合分布的一种非常有效的方法，近些年在水文分析计算领域中得到了广泛应用。熊立华等 [ 8 ] 采用Copula函数构造了两站年最大洪水联合分布函数，并且对经验概率和理论概率进行了拟合。李天元等 [ 9 ] 通过双参数Copula函数构造了两江洪峰之间的联合分布函数，估计了洪峰遭遇风险的条件概率。刘章君等 [ 10 ] 根据Copula函数对设计洪水地区组成进行区间估计，计算出地区组成的点估计值，并对估计的不确定性进行评价。

为此，本文结合Copula函数易于构造联合分布的特点，将Copula函数在构造联合分布函数的优势与频率组合法在计算多维联合分布函数的功能相结合，以雅砻江、鲜水河和庆大河洪水遭遇为例，进行了验证分析。相对于传统的频率组合法，此方法可以充分考虑干支流的随机性及相关性，且计算过程中不会出现信息失真；对变量边缘分布类型没有要求，适用性强；解决了频率组合法在多维问题上的求解问题，使其不再局限于求解二维问题。

Copula函数由Sklar在1959年提出，是定义在[0,1]上均匀分布的联合分布函数。其中Archimedean Copula函数结构简单，适应性强，在水文多变量联合分布中的应用较为广泛。n维Archimedean Copula函数表达式如下：

C ( u , φ ) = φ { φ − 1 ( u 1 ) + φ − 1 ( u 2 ) + ⋯ + φ − 1 ( u n ) } , u ∈ [ 0 , 1 ] (1)

其中，算子 φ 决定了Copula函数的类型。双参数三维Archimedean Copula函数，表达式如下：

C ( u 1 , u 2 , u 3 ) = φ 1 ( φ 1 − 1 ( u 1 ) + φ 1 − 1 φ 2 ( φ 2 − 1 ( u 2 ) + φ 2 − 1 ( u 3 ) ) ) (2)

双参数三维Archimedean Copula函数拥有两个发生元，因此可以较好地刻画出变量之间的相依结构。三维GH Copula、Clayton Copula、Frank Copula函数表达式分别如下：

C ( u 1 , u 2 , u 3 ; θ 1 , θ 2 ) = exp { − ( [ ( − log u 3 ) θ 2 + ( − log u 2 ) θ 2 ] θ 1 / θ 2 + ( − log u 1 ) θ 1 ) 1 / θ 1 } (3)

C ( u 1 , u 2 , u 3 ; θ 1 , θ 2 ) = [ ( u 3 − θ 2 + u 2 − θ 2 − 1 ) θ 1 / θ 2 + u 1 − θ 1 − 1 ] − 1 / θ 1 (4)

C ( u 1 , u 2 , u 3 ; θ 1 , θ 2 ) = − θ 1 − 1 log { 1 − c 1 − 1 ( 1 − [ 1 − c 2 − 1 ( 1 − e − θ 2 u 3 ) ( 1 − e − θ 2 u 2 ) ] θ 1 / θ 2 ) ( 1 − e − θ 1 u 1 ) } (5)

其中， c 1 = 1 − e − θ 1 ， c 2 = 1 − e − θ 2 ，参数 θ 1 用来刻画 ( u 1 , u 2 ) 、 ( u 1 , u 3 ) 之间的相依性，参数 θ 2 用来刻画 ( u 2 , u 3 ) 之间的相依性。

当防洪控制点上游有多条区间洪水时，由于区间洪水这一不确定性因素，必须考虑区间洪水与干流洪水组合后带来的风险，其组合形式有很多种，为下游防洪控制点的防洪安全带来了极大的威胁 [ 11 ] [ 12 ] [ 13 ]。

图1. 洪水遭遇组合示意图

如图1所示，干流洪水为 X 1 ，区间洪水为 X 2 , X 3 , ⋯ , X n ，防洪控制点洪水为 ，则 W = X 1 + X 2 + ⋯ + X n 。 X 1 , X 2 , ⋯ , X n 之间可以是相关关系或相互独立。 n − 1 条区间洪水在防洪控制点上游相继汇入，为下游防洪控制点的防洪安全带来了一定的威胁，为保证下游防洪控制点的防洪安全，需要求出下游防洪控制点组合洪水的分布函数，根据频率组合法，组合洪水洪峰流量超越下游安全泄量的概率分布函数可以表示为：

P ( w ≥ w 0 ) = 1 − ∫ ⋯ ∫ D h ( x 1 , x 2 , ⋯ , x n ) d x 1 d x 2 ⋯ d x n (6)

式中，D为 x 1 + x 2 + ⋯ + x n ≤ w 0 ， h ( x 1 , x 2 , ⋯ , x n ) 是随机变量 x 1 , x 2 , ⋯ , x n 的联合概率密度函数。当 w 0 为下游防洪控制点的安全泄量时，防洪控制点发生洪水的概率 P ( w ≥ w 0 ) 即为防洪风险率。

公式(6)中的联合概率密度函数 h ( x 1 , x 2 , ⋯ , x n ) 可以用Copula函数来表示，公式如下：

H ( x 1 , x 2 , ⋯ , x n ) = C [ F ( x 1 ) , F ( x 2 ) , ⋯ , F ( x n ) ] (7)

h ( x 1 , x 2 , ⋯ , x n ) = ∂ n H ( x 1 , x 2 , ⋯ , x n ) ∂ x 1 ∂ x 2 ⋯ ∂ x n = ∂ n C ( u 1 , u 2 , ⋯ , u n ) ∂ x 1 ∂ x 2 ⋯ ∂ x n (8)

式中， H ( x 1 , x 2 , ⋯ , x n ) 为随机变量 x 1 , x 2 , ⋯ , x n 的联合分布函数， C [ F ( x 1 ) , F ( x 2 ) , ⋯ , F ( x n ) ] 为随机变量 x 1 , x 2 , ⋯ , x n 的Copula函数， u 1 = F ( x 1 ) ,   u 2 = F ( x 2 ) ,   ⋯ , u n = F ( x n ) 分别为 x 1 , x 2 , ⋯ , x n 的边缘分布函数。其中，

∂ n C ( u 1 , u 2 , ⋯ , u n ) ∂ x 1 ∂ x 2 ⋯ ∂ x n = ∂ n C ( u 1 , u 2 , ⋯ , u n ) ∂ u 1 ∂ u 2 ⋯ ∂ u n ⋅ ∂ u 1 ∂ x 1 ⋅ ∂ u 2 ∂ x 2 ⋯ ∂ u n ∂ x n = c ( u 1 , u 2 , ⋯ , u n ) ⋅ f x 1 ( x 1 ) ⋅ f x 2 ( x 2 ) ⋯ f x n ( x n ) (9)

式中， c ( u 1 , u 2 , ⋯ , u n ) 为 C ( u 1 , u 2 , ⋯ , u n ) 的密度函数， f x 1 ( x 1 ) , f x 2 ( x 2 ) , ⋯ , f x n ( x n ) 分别为 x 1 , x 2 , ⋯ , x n 的边缘密度函数。

由上述公式可知，下游防洪控制点组合洪水的分布函数可由Copula函数密度函数以及边缘分布密度函数表示，因此利用Copula函数易于将多个边缘分布构造成多维联合分布的特点，对公式(6)进行改进可以得到新的n维组合洪水的分布函数表达式，如式(10)所示：

P ( w ≥ w 0 ) = 1 − ∫ ⋯ ∫ D c ( u 1 , u 2 , ⋯ , u n ) ⋅ f x 1 ( x 1 ) ⋅ f x 2 ( x 2 ) ⋯ f x n ( x n ) d x 1 d x 2 ⋯ d x n (10)

式(10)具有以下优点：对于组合洪水问题推广至n维，适用性较强；随机变量的联合概率密度函数用Copula函数表示，可以充分考虑变量间的相关性；对于边缘分布类型没有限制，任意边缘分布都可以构造联合分布。

函数求解步骤如下：

①确定边缘分布，即洪水流量的分布函数。

②在分析边缘分布相关性的基础上，选择合适的Copula函数。

③对Copula函数密度函数以及边缘分布密度函数的乘积求n重积分。

式(10)具有以下优点：对于组合洪水问题推广至n维，适用性较强；随机变量的联合概率密度函数用边缘分布密度函数以及Copula函数密度函数表示，可以充分考虑变量间的相关性；对于边缘分布类型没有限制，任意边缘分布都可以构造联合分布。

雅砻江流域位于青藏高原东部，整个流域呈南北向条带状，流域平均长度约950 km，平均宽度约137 km，河系呈羽状发育。支流鲜水河与庆大河相继汇入雅砻江，其流域图如图2所示，对下游防洪控制点的防洪安全造成了一定的威胁，本文采用雅砻江新龙站、鲜水河道孚站、庆大河扎巴站历史洪水资料，选取同一时刻三个测站的历史流量组成三元序列对作为样本，充分考虑洪水遭遇的可能性，对下游组合洪水分布函数以及防洪风险率进行计算。

图2. 雅砻江流域图

本文基于Copula函数和频率组合法，推导出n维组合洪水分布函数表达式，并构造出洪水遭遇风险估计模型。以雅砻江洪峰、鲜水河洪峰和庆大河洪峰为例，在分析其相关性基础上，采用双参数三维GH Copula、Clayton Copula、Frank Copula函数对其洪峰联合分布进行计算，通过对比分析，GH Copula函数对变量的拟合效果最优。建立了基于双参数三维GH Copula函数的三河流洪峰联合分布，然后计算得到组合洪水分布曲线，通过组合洪水分布曲线可以得到下游防洪控制点的防洪风险率以及各组合洪水发生的概率。通过实例分析证明本文所构造的模型具有可行性，且推导出的n维组合洪水分布函数对于易发生洪水遭遇地区的防洪安全提供了一定的参考。

国家自然科学基金(51709105)；中央高校基本科研业务费专项资金资助(2019MS031)。

梁希金,王丽萍,张验科,张佳新,叶凯华. 基于三维Copula函数的组合洪水研究Combined Flood Research Based on Three-Dimensional Copula Function[J]. 水资源研究, 2020, 09(01): 33-41. https://doi.org/10.12677/JWRR.2020.91004