干旱一般可以用干旱历时(drought duration)、干旱严重程度(drought severity)和干旱强度(drought intensity)三个特征来表示[9],而这可以用游程理论来得到。一般而言,在一个有限取值的序列中,满足一定条件的同一符号的一个连串称之为一个“游程”,一个游程中同一符号出现的次数称之为游程的长度。有研究认为长历时的SPI小于0的事件也会导致严重的水资源供给及与干旱有关的问题,取0作为游程的截断水平[10],本文采用了与其相同的截断水平(图2)。
当SPI≤0时定义为发生了一次干旱事件,游程长度为干旱历时D(单位为月),在游程中累积的SPI值与临界值的缺乏值为干旱严重程度(即为干旱事件中累积的干旱缺水状况)S,即,其中干旱严重程度与干旱历时的比值为干旱强度(即为干旱事件中平均的干旱缺水状况)。
干旱历时、干旱严重程度和干旱强度,任一个可以由另外两个得到,因而可以只用其中两个特征来表示一次干旱,本文选取的是直接得到的干旱历时和干旱严重程度。
一般来说,当将干旱历时看做离散型分布时,干旱历时的分布可通过几何分布来拟合[11],而如果将干旱历时分布看做连续型分布时,常用指数分布来拟合[10],本文用指数分布来拟合干旱历时。干旱严重程度一般用伽马分布函数来拟合[10-11],为使结果更可靠,本文还同时采用在水文领域应用比较广泛的三变量广义极值分布、三变量广义帕累托分布、皮尔逊三型分布、三变量对数正态分布、五变量Wakeby分布来拟合。各分布函数的参数估计采用线性矩估计[12],最后使用Kolmogorov-Smirnov(K-S)[13]方法进行拟合优度检验,确定最适分布函数。