基于Pair-copula函数的文干旱频率分析—以南盘江流域为例_杨茂灵.docVIP

基于Pair-copula函数的文干旱频率分析—以南盘江流域为例_杨茂灵 　 第２３卷第９期２０１４年９月 长江流域资源与环境 ＲｅｓｏｕｒｃｅｓａｎｄＥｎｖｉｒｏｎｍｅｎｔｉｎｔｈｅＹａｎｔｚｅＢａｓｉｎ ｇ Ｖｏｌ．２３Ｎｏ．９ 　 Ｓｅｔ．２０１４ｐ 基于Ｐａｉｒｃｏｕｌａ函数和标准化径流－ｐ 指数的水文干旱频率分析 ———以南盘江流域为例 杨茂灵，王　龙＊，余　航，高　瑞，雷腾云 （）云南农业大学水利水电与建筑学院，云南昆明６５０２０１ （摘　要：选取南盘江流域３个水文站５的逐月径流资料，联合两个相邻站点，分别计算标准化２ａ１９６１～２０１２年）），运用游程理论识别５并用Ｐ径流指数（ＳＳＦＩ２ａ来不同站点干旱特征，ａｉｒｃｏｕｌａ函数计算不同站点的联合重现期。－ｐ（）结果表明：从单个站点看，１５２ａ来３个站点干旱特征值的最大值都出现在２０１１～２０１２年；５２ａ来沾益站的干旱次数、最大干旱历时、最大干旱烈度均大于高谷马站和江边街站，最大干旱峰值出现在江边街站；从联合站点看，中游以上高谷马站和沾益站的干旱特征值均大于中游以下高谷马站和江边街站，３个站点联合识别的干旱特征值大（小于中上游站点联合识别干旱特征值，上游干旱更严重；于中下游站点联合识别干旱特征值，２）ＦｒａｎｋＰａｉｒｃｏｕｌａ　－ｐ（函数的ＲＭ拟合程度最高，适合运用于南盘江流域干旱频率分析中；运用ＦＳＥ、ＡＩＣ和ＢＩＣ值最小，３）ｒａｎｋＰａｉｒ　－，是５ｃｏｕｌａ函数计算得到３个水文站点２００９～２０１２年的连续干旱重现期都在１００～２００ａ２ａ来最严重的一次ｐ干旱。 关键词：游程理论；南盘江流域；重现期；标准化径流指数；干旱频率Ｐａｉｒｃｏｕｌａ函数；－ｐ（）中图分类号：Ｐ３３８＋．６　　文献标识码：Ａ　　文章编号：１００４８２２７２０１４０９１３１５０７－－－：／ＤＯＩ１０．１１８７０ｃｌｚｈ２０１４０９０１９ｊｙｙｙｊ 水循环过程加速，降水时　　在全球变暖的背景下， １］ 。自然降水不空格局及其结构将发生显著变化［ 足，以及降雨量的季节分配比例失调或时空分布不均是干旱形成的直接原因，大气环流异常、地理环境因素、植被生态状况和人类生产活动等都对干旱的形成和发展产生不同程度的影响。云南属于季风气 ２］ ，而且是高原地形；喀斯特地貌，地表水下渗［容候， 模拟，模拟最大负轮次和轮次序列，估计干旱历时的重现期。然而水文干旱的发生与众多因素有关，单一变量存在一定的局限性，不能客观全面的反映整个干旱事件，因此，多变量分析就成为人们的必然选择。总体上多变量水文频率的传统研究方法主要有正态变换法、经验频率法和同边缘分布构ＦＧＭ法、 ６］ ，建法等［但是这些方法都存在一定不足，在实际水 易引起区域性干旱以及极端干旱事件的发生，其中分析研究水文干旱频率水文干旱是重要类型之一， 越来越受到人们的重视。干旱频率的分析方法主要有单变量水文频率分析和多变量水文干旱特征研 ［］ 究。Ｙｅｖｅｖｉｃｈ３在１９６７年最早应用游程理论研究ｊ 文干旱分析中有些假定条件很难得到满足，因此，上述这些方法的有效性和应用范围受到很大的限制。为了克服这些困难和缺点，国内外学者将Ｃｏｕｌａ函ｐ数理论与方法应用于多变量水文干旱分析计算领 ５］ 。域［ 定义了干旱历时、烈度和强度等水文意义干旱特性， ４］ 上的干旱指标；马秀峰［分析了随机序列轮长与轮Ｃｏｕｌａ函数的主要特点在于各个单因子变量ｐ 的边缘分布可以采取任何形式，变量之间可以具有各种相关关系，具有极强的灵活性和适应性等特点。 ［］ 洪量的Ｓｈｉａｕ７根据Ｃｏｕｌａ函数分析了二维洪峰、ｐ 次的统计规律，指出平稳随机序列的轮次服从二阶 ５］ 丁晶等［利用径流和降水随机线性齐次差分方程； ；收稿日期：修回日期：２０１３１０２１２０１３１１２７－－－－ ）；）基金项目：水利部公益性行业专项云南旱灾应急响应系统研究项目（云南省教育厅科学研究基金项目（２０１００１０４４２０１３Ｙ４７７，：作者简介：杨茂灵（女，硕士研究生，主要从事区域水资源开发与利用方向研究．１９８９～　）Ｅ－ｍａｉｌ２６．ｃｏｍｔ８２７９７＠１ｙｙ：ａｉｌｎｎｄｗｌ２６．ｃｏｍ＊通讯作者Ｅ－ｍ＠１ｙ １３１６ 长江流域资源与环境 第２３卷　 联合分布，并求得相应的联合分布概率进行频率分 ［］析；ＳａｌｖａｔｏｒｅＧｒｉｍａ