基于Copula函数的丰枯遭遇概率分析

      将年径流丰枯等级划分为丰水、平水、枯水3级。根据遭遇分析方法，定义累积概率 Pf＝62.5%为丰水与平水的分割频率，Pk=37.5%为平水与枯水的分割频率。丰枯遭遇情况总体概括为以下6种：丰丰遭遇、丰枯遭遇、平丰遭遇、平枯遭遇、平平遭遇、枯枯遭遇。

具体计算概念图如下：

相应的遭遇概率公式如下：

根据遭遇概率，结合联合分布和生存分布，可以由R语言代码计算如下：

#-------丰丰遭遇P(x>0.625,y>0.625)
para1<-3.5797
ff<-surCOP(0.375,0.375,cop=GHcop,para=para1)
1-0.625-0.625+GHcop(0.625,0.625,para = para1)

#-------丰枯遭遇P(x>0.625,y<0.375)+P(x<0.375,y>0.625)
fk<-0.375+0.375-GHcop(0.375,0.625,para = para1)-GHcop(0.625,0.375,para = para1)

#-------平丰遭遇P(0.375<x<0.625,y>0.625)+P(x>0.625,0.375<y<0.625)
pf<-1-GHcop(0.625,0.625,para = para1)-ff-fk


#-------平枯遭遇P(0.375<x<0.625,y<0.375)+P(x<0.375,0.375<y<0.625)
pk<-2*(GHcop(0.625,0.375,para = para1)-GHcop(0.375,0.375,para = para1))

#-------平平遭遇P(0.375<x<0.625,0.375<y<0.625)
pp<-GHcop(0.625,0.625,para = para1)-GHcop(0.375,0.375,para = para1)-pk

#------枯枯遭遇P(x<0.375,y<0.375)
kk<-GHcop(0.375,0.375,para = para1)

#-----检验总概率是否为1#####
ff+fk+pf+pk+pp+kk