基于Copula函数的嘉陵江流域水沙丰枯遭遇关系探究

为了深入了解嘉陵江水沙关系的变化规律。基于嘉陵江出口北碚水文站1965-2022年的年径流量和年输沙量资料,利用Mann-Kendall检验法、双累积曲线法并结合物理成因法诊断水沙关系的变异点,通过AIC、BIC、RMSE对径流量和输沙量的边缘分布函数和Copula函数进行优选,确定最优Copula函数并建立嘉陵江流域水沙联合分布模型,运用模型对流域内水沙丰枯遭遇概率、联合和同现重现期进行计算,进而对比分析水沙单变量和联合变量设计值。结果表明:①水沙关系在1989年发生突变,变异前的水沙相关性强于变异后的水沙相关性;②1965-1989年阶段水沙序列最优边缘分布都为Logn,最优水沙联合分布模型为Gumbel Copula模型。1990-2022年阶段径流量和输沙量的最优边缘分布分别为Gamma、GP,最优水沙联合分布模型为Frank Copula模型;③2个阶段水沙序列同步概率分别是59.77%、61.14%,均大于异步概率,且出现极端情况“水丰沙枯”、“水枯沙丰”类型的概率极低,表明水沙概率存在较大的相关性,来水来沙条件相对稳定;④联合变量设计值大于单变量设计值,通过两变量联合分布计算的径流量和输沙量设计值更加可靠。水沙联合重现期集中在2年附近,同现重现期绝大部分小于50年。通过探讨嘉陵江流域来水来沙联合变化特征,可为流域的水沙调控、防灾减灾等工作提供参考。