《基于copula函数的区域干旱分析方法 Regional drought analysis approach based on copula function》.pdfVIP

第21卷第2期 水科学进展 V01．21，No．2 2010年3月 ADVANCESINWATERSCIENCE Mar．，2010 基于copula函数的区域干旱分析方法 陆桂华，闫桂霞，吴志勇，康燕霞 河海大学水文水资源与水利工程科学国家重点实验室，江苏南京210098 摘要：提出不同网格干旱等级影响下的区域干旱指数。采用copula函数，建立了区域干旱历时和干旱强度的联合 分布，计算联合分布的重现期，并对实际重现期作区间估计。所提出的区域干旱指数能够反映研究区域内发生干旱 的严重程度和该区内不同干旱等级所影响面积的大小，且联合分布同时考虑了干旱历时和干旱强度，更全面地反 映区域的干旱状况。重庆市2006年干旱实证分析表明，区域干旱指数能够较好地描述本次干旱的发生发展过程， 并计算出联合分布的重现期为118年。 关键词：区域干旱；copula函数；重现期；重庆 中图分类号：P467 文献标志码：A 文章编号：1001-679l 2010 02-0188-06 通常干旱的特征分析，仅仅局限于单变量分析。包括干旱历时和干旱强度等多个相关变量的极值水文事 件，单变量分析难以反映多变量特征之间的联合分布关系。并且由于干旱历时和干旱强度是显著相关且服从 不同分布的两个主要干旱特征，已有的具有相同边缘分布的两变量分布模型如两变量对数分布模型等将不再 适用，而具有灵活边缘分布的两变量分布模型往往具有有限的适用范围，如FGM模型只适用于弱相关的随 一种新的途径。copula函数的主要特点在于各单因子变量的边缘分布可以采用任何形式，变量之间可以具有 各种相关关系，具有极强的灵活性和适应性¨J。Shiau、闫宝伟等”刮分别利用此方法对黄河流域和汉江上游 地区的干旱特征进行了分析。但对于区域的干旱研究往往是以区域内单个站点或多个站点的平均值为基础 km×50 的，没有考虑不同站点不同干旱等级对区域干旱的影响。本文基于分辨率为50 km网格的综合气象 干旱指数1，提出区域干旱指数；然后，统计研究区域内的干旱场次，计算干旱历时和干旱强度，并利用 Gumbel—Hougaardcopula函数建立干旱历时和强度的联合分布，借助重现期对研究区域的干旱特征进行分析； 最后，采用本文提出的区域干旱分析方法对重庆市2006年干旱进行实证分析。 1 区域干旱指数 基于帕尔默干旱指数 Palmer DroughtSeverity Precipitation 实际水分供应持续少于气候适宜水分供应的水分亏缺。它能够较好地反映区域干旱受旱／成灾范围，以及河 道径流的丰枯状况，其计算式为”1 ， aPl+6P2+cP3 1 收稿日期：2009-04-29 目 座机电话号码9：座机电话号码1 作者简介：陆桂华 1962一 ，男，浙江东阳人，博士生导师，教授，主要从事水文水资源研究。 E—mail：lugh@hhu．edu．ell 万方数据 第2期 陆桂华，等：基于copula函数的区域干旱分析方法 189 除以历史上出现最小指数值的平均值。，的变化范围为0—1，划分的干旱标准是，无旱：J 0．35；轻微干旱： 0．35≤， 0．25；中等干旱：0．25≤， 0．15；严重干旱：0．15≤， 0．10；极端干旱：11 0．10。 对于综合气象干旱指数，，0．35为发生干旱的临界值，， 值越小干旱越严重。对于某一研究区域 如图1 ，区域面积为 A，区域内任意网格i的干旱指数，i小于0．35，则该网格处于 干旱状态，其在区域内的面积为A；，设网格i对区域干旱的 贡献为 0．35一，； A；，若区域内有n个网格处于干旱状态，则 n 图中阴影部分的体积y ’ 0．35一，； A；即为区域的干旱程 篇 度。为了对比分析不同区域的干旱强度，按下式定义研究区域 第．『时刻的区域干旱指数： 图l 区域干旱指数定义图 VOOI‘A iiI一53．o ∑001 Definitionofthe index 2 Fig．1 regionaldrought B T —鱼j二 式中 Rj．为研究区第J月份的区域干旱指数；，Ⅱ为研究区域内处于干旱状态的第i个网格第．『月份的综合气象 干旱指数值。由于干旱指数，在0—1之间变化，因此，区域干旱指数尺的变化范围为0～35，区域干旱指数 越大，区域干旱越严重。 根据游程理论，设R。为阈值，当尺大于或等于R。时即发生干旱，正的游程长度为干旱历时D，游程总 量为干旱强度S。本文假设如果两次干旱过程 干旱历时和干旱强度分别为d，，d：和s。，s： 中间有且只有1 个月份的干旱指数值小于R。且大于R。 R。小于R。 ， ￡ D D ● 。吐 则这两次干旱是从属干旱¨o，可合并成一次干旱过程， ^d。