干旱发生频率及重现期理论与应用研究初步报告课件.ppt

1、水利部水利水电规划设计总院水利部水利水电规划设计总院合肥工业大学合肥工业大学20122012年年2 2月月汇汇 报报 提提 纲纲一、研究背景一、研究背景一、研究背景一、研究背景（一）研究目的与意义（一）研究目的与意义 指标值的演变态势，利用干旱频率分析判指标值的演变态势，利用干旱频率分析判别面临干旱的重现期，据此判别抗旱预案别面临干旱的重现期，据此判别抗旱预案的响应级别的响应级别l在旱灾风险管理的作用在旱灾风险管理的作用：干旱频率分析是：干旱频率分析是干旱灾害危险性分析的核心内容，是干旱干旱灾害危险性分析的核心内容，是干旱灾害风险管理的重要组成部分灾害风险管理的重要组成部分一、研究背景一、研究

2、背景（二）国内外研究进展（二）国内外研究进展 一、研究背景一、研究背景（二）国内外研究进展（二）国内外研究进展 社会经济干旱社会经济干旱是指由自然系统与人类社会是指由自然系统与人类社会经济系统中水资源供需不平衡造成的异常经济系统中水资源供需不平衡造成的异常水分短缺现象，如由于用水管理的实际操水分短缺现象，如由于用水管理的实际操作或设施的破坏引起的缺水作或设施的破坏引起的缺水一、研究背景一、研究背景（二）国内外研究进展（二）国内外研究进展 方法对这些一般干旱进行合并。目前提出了方法对这些一般干旱进行合并。目前提出了干旱时间间隔准则、水量准则、综合考虑时干旱时间间隔准则、水量准则、综合考虑时间间隔

3、和水量两个准则三类干旱合并方法间间隔和水量两个准则三类干旱合并方法一、研究背景一、研究背景（二）国内外研究进展（二）国内外研究进展 并依据重现期公式来设计洪水、干旱。点并依据重现期公式来设计洪水、干旱。点干旱频率分析是通过以上步骤来描述各站干旱频率分析是通过以上步骤来描述各站点干旱特征的概率特征，计算干旱发生概点干旱特征的概率特征，计算干旱发生概率、重现期、得到设计干旱成果率、重现期、得到设计干旱成果一、研究背景一、研究背景（二）国内外研究进展（二）国内外研究进展 分布的尾部信息，而且参数估计也较为方便分布的尾部信息，而且参数估计也较为方便等特点，目前它已成为多变量分析中最为广等特点，目前它已

4、成为多变量分析中最为广泛应用的方法，也是多变量干旱频率分析中泛应用的方法，也是多变量干旱频率分析中使用最多的方法使用最多的方法一、研究背景一、研究背景（二）国内外研究进展（二）国内外研究进展 分成多个一致区，对各一致区分别进行区域分成多个一致区，对各一致区分别进行区域频率分析频率分析。主要包括一致区划分和区域频率。主要包括一致区划分和区域频率分析两个步骤。一致区划分主要采用基于各分析两个步骤。一致区划分主要采用基于各种水文、气象、地形、地貌、干旱特征指标种水文、气象、地形、地貌、干旱特征指标的聚类方法。该途径能大致反映干旱的空间的聚类方法。该途径能大致反映干旱的空间分布和特性，也能较好地反映干

5、旱的区域特分布和特性，也能较好地反映干旱的区域特性，若能进一步配合有关干旱特性的数量指性，若能进一步配合有关干旱特性的数量指标的等值线图等，则可望能很好地描述区域标的等值线图等，则可望能很好地描述区域干旱特征干旱特征一、研究背景一、研究背景（三）研究目标与工作内容（三）研究目标与工作内容 及演变特征，为干旱区划和旱灾风险及演变特征，为干旱区划和旱灾风险管理提供相关指标管理提供相关指标一、研究背景一、研究背景（三）研究目标与工作内容（三）研究目标与工作内容 将干旱频率分析的研究成果，应用到试点将干旱频率分析的研究成果，应用到试点区域、干旱区划和旱灾风险分析中区域、干旱区划和旱灾风险分析中一、研究

6、背景一、研究背景经济经济社会资料社会资料自然地理资料自然地理资料水文气象资料水文气象资料旱情旱灾资料旱情旱灾资料干旱指标选择和干旱特征变量提取（游程理论）干旱指标选择和干旱特征变量提取（游程理论）干旱历时干旱历时干旱烈度干旱烈度基于基于copula函数函数的单个干旱指标的干旱频率分析的单个干旱指标的干旱频率分析多个干旱指多个干旱指标的干旱频标的干旱频率分析率分析基于小区域干旱频率分基于小区域干旱频率分析结果的大区域干旱重析结果的大区域干旱重现期计算现期计算频率分析计算结果频率分析计算结果与区域实际旱情和与区域实际旱情和灾情的对比分析灾情的对比分析历史旱情信息历史旱情信息干旱发生频率及重现期理论

7、与应用研究的技术路线图干旱发生频率及重现期理论与应用研究的技术路线图二、干旱发生频率及重现期理论研究二、干旱发生频率及重现期理论研究 （一）干旱发生频率的概念分析（一）干旱发生频率的概念分析（二）干旱过程识别理论（二）干旱过程识别理论 （三）基于（三）基于copula的干旱频率分析理论的干旱频率分析理论 二、干旱发生频率及重现期理论研究二、干旱发生频率及重现期理论研究 （一）干旱发生频率的概念分析（一）干旱发生频率的概念分析 系列，用联合概率分布函数分析这些干旱系列，用联合概率分布函数分析这些干旱特征与出现频率之间的定量关系，据此关特征与出现频率之间的定量关系，据此关系得到某一干旱频率对应的干

8、旱特征值、系得到某一干旱频率对应的干旱特征值、判别典型干旱年的重现期判别典型干旱年的重现期二、干旱发生频率及重现期理论研究二、干旱发生频率及重现期理论研究 （一）干旱发生频率的概念分析（一）干旱发生频率的概念分析 二、干旱发生频率及重现期理论研究二、干旱发生频率及重现期理论研究 二、干旱发生频率及重现期理论研究二、干旱发生频率及重现期理论研究 二、干旱发生频率及重现期理论研究二、干旱发生频率及重现期理论研究 理论依据：理论依据：在给定阈值下运用游程理论从干旱指标过程中识别干旱特在给定阈值下运用游程理论从干旱指标过程中识别干旱特征征干旱历时干旱历时是指一次干旱过程开始到结束的时间是指一次干旱过程

9、开始到结束的时间；干旱烈度；干旱烈度是指阈值指是指阈值指标与干旱指标之差的累积和标与干旱指标之差的累积和难点：难点：指标阈值、干旱特征变量的确定，特征变量的样本满足频率分指标阈值、干旱特征变量的确定，特征变量的样本满足频率分析的基本假设和样本容量的要求析的基本假设和样本容量的要求 二、干旱发生频率及重现期理论研究二、干旱发生频率及重现期理论研究 二、干旱发生频率及重现期理论研究二、干旱发生频率及重现期理论研究 干旱历时干旱历时D和干旱烈度和干旱烈度S一般可分别用指数分布和一般可分别用指数分布和GammaGamma分分布描述布描述当直接用当直接用样本样本数据估计各参数来确定特征变量的分布时，数据

10、估计各参数来确定特征变量的分布时，往往会出现不符合实际的情况，造成计算结果的不合理往往会出现不符合实际的情况，造成计算结果的不合理在基于在基于样本样本数据估计分布函数参数基础上，通过适线法进数据估计分布函数参数基础上，通过适线法进行参数估计的调整行参数估计的调整/()1edDFd10( ) e()ssSsFsds二、干旱发生频率及重现期理论研究二、干旱发生频率及重现期理论研究 相关性结构相关性结构 Copula函数函数C 干旱历时和干旱烈干旱历时和干旱烈度联合分布度联合分布 干旱烈度边际干旱烈度边际分布分布干旱历时边际干旱历时边际分布分布 二、干旱发生频率及重现期理论研究二、干旱发生频率及重现

11、期理论研究 1/,( , )( , ) exp ( ln )( ln ) DSFd sCuvuv 1/,( , )( , )(1)D SFd sC u vuv,1(e1)(e1)(,)(,)ln 1(e1)uvDSFdsCuv 二、干旱发生频率及重现期理论研究二、干旱发生频率及重现期理论研究 ( )0.44( )0.12om iP in211( )( ) ncoiR M S EPiPin二、干旱发生频率及重现期理论研究二、干旱发生频率及重现期理论研究 ,( )( )( , )1( )( )( , )aDSD SE LE LT d sP DdSsFdF sFd s历时与干旱烈度的联合概率分布历时

12、与干旱烈度的联合概率分布二、干旱发生频率及重现期理论研究二、干旱发生频率及重现期理论研究 根据研究区域干旱的类型、成因特点，以及降水量、蒸发量、根据研究区域干旱的类型、成因特点，以及降水量、蒸发量、径流量、土壤含水量、地下水等干旱指标影响研究区域干旱的径流量、土壤含水量、地下水等干旱指标影响研究区域干旱的重要性以及专家咨询等信息，确定降水量和蒸发量、径流量、重要性以及专家咨询等信息，确定降水量和蒸发量、径流量、土壤含水量等干旱指标的权重土壤含水量等干旱指标的权重用这些干旱指标的干旱频率分析的重现期计算结果和这些干旱用这些干旱指标的干旱频率分析的重现期计算结果和这些干旱指标的权重进行加权和，得到

13、基于多个干旱指标的研究区域干指标的权重进行加权和，得到基于多个干旱指标的研究区域干旱频率分析的干旱合成重现期计算结果旱频率分析的干旱合成重现期计算结果二、干旱发生频率及重现期理论研究二、干旱发生频率及重现期理论研究 综合分析综合分析单个干旱单个干旱指标频率指标频率分析结果分析结果Copula多个干旱多个干旱指标频率指标频率分析结果分析结果降水量降水量蒸发量蒸发量径流量径流量土壤含水量土壤含水量 根据各小区域（如县）的面积、人口、根据各小区域（如县）的面积、人口、GDPGDP、政治社会影响等指、政治社会影响等指标以及专家咨询等信息，确定各小区域的权重标以及专家咨询等信息，确定各小区域的权重用这些

14、小区域的干旱合成重现期计算结果和这些小区域的权重用这些小区域的干旱合成重现期计算结果和这些小区域的权重进行加权和，得到大区域（如省）的干旱综合重现期计算结果进行加权和，得到大区域（如省）的干旱综合重现期计算结果根据各级区域的干旱重现期计算结果，进行相应区域的抗旱管根据各级区域的干旱重现期计算结果，进行相应区域的抗旱管理理二、干旱发生频率及重现期理论研究二、干旱发生频率及重现期理论研究 专家咨询专家咨询 二、干旱发生频率及重现期理论研究二、干旱发生频率及重现期理论研究 南方水文干旱试点区域（例如湖北邵阳、江西鄱阳湖）的南方水文干旱试点区域（例如湖北邵阳、江西鄱阳湖）的降水量和蒸发量、径流量资料系

15、列降水量和蒸发量、径流量资料系列中部农业干旱试点区域（例如安徽）的降水量和蒸发量、中部农业干旱试点区域（例如安徽）的降水量和蒸发量、土壤含水量、地下水资料系列土壤含水量、地下水资料系列北方气象干旱试点区域（例如辽宁）的降水量和蒸发量资北方气象干旱试点区域（例如辽宁）的降水量和蒸发量资料系列料系列二、干旱发生频率及重现期理论研究二、干旱发生频率及重现期理论研究 三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （一）基于降水指标的干旱频率分析（一）基于降水指标的干旱频率分析 （二）基于降水和蒸发指标的干旱频率分析（二）基于降水和蒸发指标的干旱频率分析 （三）基于径流指

16、标的干旱频率分析（三）基于径流指标的干旱频率分析 （四）基于地下水指标的干旱频率分析（四）基于地下水指标的干旱频率分析 （五）基于土壤含水量的干旱频率分析（五）基于土壤含水量的干旱频率分析 （六）基于干旱综合（六）基于干旱综合Z Z指数的干旱频率分析指数的干旱频率分析 三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 两个典型研究区域两个典型研究区域 邵阳市是中国南方山地丘陵地区，属于湖南省邵阳市是中国南方山地丘陵地区，属于湖南省“衡邵干旱走廊衡邵干旱走廊”的一部分，由于邵阳地区降雨时空分布不均，干旱灾害频繁的一部分，由于邵阳地区降雨时空分布不均，干旱灾害频繁 三、

17、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 两个典型研究区域两个典型研究区域 淮北平原位于安徽省北部，主要包括界首、临泉、太和、阜南、阜淮北平原位于安徽省北部，主要包括界首、临泉、太和、阜南、阜阳、亳州、涡阳、利辛、蒙城、淮北、濉溪（包括淮北）、宿州、砀阳、亳州、涡阳、利辛、蒙城、淮北、濉溪（包括淮北）、宿州、砀山、萧县、灵璧和泗县等山、萧县、灵璧和泗县等15个市县。由于地处我国南北气候的过渡带，个市县。由于地处我国南北气候的过渡带，此区域内的降水量年际变化大、年内分配不均，每年旱涝灾害频发此区域内的降水量年际变化大、年内分配不均，每年旱涝灾害频发宿州濉溪亳州蒙城

18、灵璧涡阳利辛阜南萧县泗县太和临泉阜阳砀山界首降水量距平百分率反映的是一段时期内的降水量与常年同期降水量距平百分率反映的是一段时期内的降水量与常年同期平均降水量之间的关系，其表示为平均降水量之间的关系，其表示为三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （一）基于降水指标的干旱频率分析（一）基于降水指标的干旱频率分析 () /1 0 0 %aPPPP根据邵阳市根据邵阳市15个雨量站的逐月降水资料，采用降水量距平百分率个雨量站的逐月降水资料，采用降水量距平百分率法对气象干旱过程进行识别，依据法对气象干旱过程进行识别，依据气象干旱等级气象干旱等级确定干旱指确定干旱指

19、标阈值分别为标阈值分别为R0取取0，R1取取- -30%，R2取取- -40%，下表为相应的气，下表为相应的气象干旱统计结果象干旱统计结果三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （一）基于降水指标的干旱频率分析（一）基于降水指标的干旱频率分析 邵阳市各雨量站干旱特征计算结果邵阳市各雨量站干旱特征计算结果站点站点干旱次干旱次数数干旱间隔干旱间隔均值均值(月月)干旱历时干旱历时均值均值(月月)干旱烈度干旱烈度均值均值(mm)历时最历时最大值大值(月月)烈度最大值烈度最大值(mm)干旱历时与烈度干旱历时与烈度的的Kendall相关相关系数系数黄桥黄桥916.20

20、2.0148.028244.010.54高沙高沙876.431.9451.278268.720.45邵东邵东787.172.0951.858185.270.45诸甲亭诸甲亭806.992.1048.878189.840.53茅坪茅坪936.021.7548.135195.580.46罗家庙罗家庙906.221.9047.499219.330.54白毛坪白毛坪866.561.8147.355132.670.41武冈武冈866.501.8549.0110208.280.49新宁新宁876.431.8353.039276.240.45小沙江小沙江806.991.9053.477209.360.49大

21、河滩大河滩886.361.8951.969186.050.48回龙寺回龙寺916.151.7845.109216.820.47六都寨六都寨906.221.9251.926216.410.49党坪党坪975.771.7847.056180.470.42三门江三门江836.742.0747.959212.660.45三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （一）基于降水指标的干旱频率分析（一）基于降水指标的干旱频率分析 上表的邵阳市各雨量站干旱特征计算结果说明：上表的邵阳市各雨量站干旱特征计算结果说明：各站干旱历时、干旱烈度均值及其各站干旱历时、干旱烈度均值及

22、其Kendall相关系数较为接近，说明各站的相关系数较为接近，说明各站的干旱平均水平大体相当干旱平均水平大体相当党坪站干旱次数最多，而邵东站干旱次数最少，两者相差党坪站干旱次数最多，而邵东站干旱次数最少，两者相差19次次各站的干旱历时、烈度最大值有较大差异。茅坪站、白茅坪站的历时最大值各站的干旱历时、烈度最大值有较大差异。茅坪站、白茅坪站的历时最大值仅为仅为5个月，而武冈站的历时最大值长达个月，而武冈站的历时最大值长达10个月；新宁站的烈度最大值为个月；新宁站的烈度最大值为276.24mm，而白茅坪站的烈度最大值仅为，而白茅坪站的烈度最大值仅为132.67mm可见，各站的干旱特征具有显著的空间

23、差异性，研究它们之间的关系对区域可见，各站的干旱特征具有显著的空间差异性，研究它们之间的关系对区域干旱频率分析具有重要意义干旱频率分析具有重要意义 三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （一）基于降水指标的干旱频率分析（一）基于降水指标的干旱频率分析 通过适线法确定各站中单个干旱特征变量的分布曲线，根据通过适线法确定各站中单个干旱特征变量的分布曲线，根据Copula函数的估计函数的估计参数，选用参数，选用Gumbel-Hougaard（GH） Copula函数对邵阳市各雨量站典型干旱函数对邵阳市各雨量站典型干旱年份的重现期进行计算，结果见下表：年份的重现

24、期进行计算，结果见下表：站名站名1956重现期重现期/a1960重现期重现期/a1963重现期重现期/a1985重现期重现期/a1998重现期重现期/a黄桥黄桥11.936.7866.173.1852.56高沙高沙23.022.0883.9621.9252.28邵东邵东15.3241.762.6919.4566.17诸甲亭诸甲亭10.987.105.866.5059.62茅坪茅坪31.8252.722.9924.0123.63罗家庙罗家庙12.618.8225.596.0068.22白毛坪白毛坪18.512.555.893.1845.23武冈武冈21.532.5645.702.7362.61新

25、宁新宁2.371.8548.823.7787.40小沙江小沙江1.543.3231.332.8013.80大河滩大河滩6.496.636.8816.007.06回龙寺回龙寺1.911.7447.509.0375.97六都寨六都寨8.6717.669.8210.1116.04党坪党坪5.021.5911.628.3823.12三门江三门江12.639.8951.504.7882.63三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （一）基于降水指标的干旱频率分析（一）基于降水指标的干旱频率分析 上表的邵阳市各雨量站典型干旱年份的重现期计算结果说明：上表的邵阳市各雨量

26、站典型干旱年份的重现期计算结果说明：1963、1998年邵阳市受旱严重，个别站点的干旱重现期超过了年邵阳市受旱严重，个别站点的干旱重现期超过了80年一年一遇遇对于同一年份，区域内各站干旱程度差异显著。例如：对于同一年份，区域内各站干旱程度差异显著。例如：1963年，干旱年，干旱严重区主要集中在邵阳市中部，西部和东部的干旱程度则相对较轻；而严重区主要集中在邵阳市中部，西部和东部的干旱程度则相对较轻；而1998年，邵阳市大部分地区干旱严重，南部各站的干旱重现期普遍高于年，邵阳市大部分地区干旱严重，南部各站的干旱重现期普遍高于北部的小沙江、六都寨和大河滩站北部的小沙江、六都寨和大河滩站1963(a)

27、、1998(b)年邵阳市干旱重现期空间分布特征 三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （一）基于降水指标的干旱频率分析（一）基于降水指标的干旱频率分析 为了更好地反映区域内各站在平均干旱水平下的干旱程度，由各站的干旱历时为了更好地反映区域内各站在平均干旱水平下的干旱程度，由各站的干旱历时均值、干旱烈度均值得干旱重现期，见下图，结果表明：均值、干旱烈度均值得干旱重现期，见下图，结果表明：邵阳市南部的平均干旱程度相对较重邵阳市南部的平均干旱程度相对较重邵阳市邵阳市中部和东北部相对较轻中部和东北部相对较轻上述干旱空间特征分析结果可为区域干旱灾害风险管理提供有效

28、的科学依据上述干旱空间特征分析结果可为区域干旱灾害风险管理提供有效的科学依据 平均干旱水平下的邵阳市干旱重现期空间分布特征 三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （一）基于降水指标的干旱频率分析（一）基于降水指标的干旱频率分析 取罗家庙站取罗家庙站19562002年的降水量距平百分率样本资料系列直接进行频率分析，见下年的降水量距平百分率样本资料系列直接进行频率分析，见下图和下表，可见：基于图和下表，可见：基于Copula的干旱历时干旱烈度多变量频率分析的计算结果较合理，的干旱历时干旱烈度多变量频率分析的计算结果较合理，降水量距平百分率单变量频率分析的结果

29、显然不合理。例如：邵阳市降水量距平百分率单变量频率分析的结果显然不合理。例如：邵阳市1998年干旱程度年干旱程度显然明显比显然明显比1963年严重，年降水量不能表达因各月降水量分布不均而导致的干旱现象，年严重，年降水量不能表达因各月降水量分布不均而导致的干旱现象，无法识别干旱过程的起止时间，也无法用以判别干旱现象发生的重现期无法识别干旱过程的起止时间，也无法用以判别干旱现象发生的重现期罗家庙站年降水量距平百分率指标下不同频率分析方法的典型干旱年份的重现期罗家庙站年降水量距平百分率指标下不同频率分析方法的典型干旱年份的重现期 年份年份干旱历时干旱烈度多变量频率分析干旱历时干旱烈度多变量频率分析/

30、a降水量距平百分率单变量频率分析降水量距平百分率单变量频率分析/a年降水量年降水量/mm195612.612.361182.719608.824.241092.5196325.59111.11917.419856.005.631062.5199868.226.501049.3罗家庙站年降水量频率曲线 P/E指数考虑了降水量指数考虑了降水量P与蒸发量与蒸发量E的比值，的比值，K指数反映的是降水量相对变指数反映的是降水量相对变率与蒸发量相对变率之间的关系：率与蒸发量相对变率之间的关系：三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （二）基于降水和蒸发指标的干旱频率分

31、析（二）基于降水和蒸发指标的干旱频率分析 根据根据气象干旱等级气象干旱等级标准及干旱分类相关研究标准及干旱分类相关研究:当以降水量距平百分率当以降水量距平百分率Pa为干旱指标时，本专题选用单月距平百分率为为干旱指标时，本专题选用单月距平百分率为- -30%和和 - -40%时所对应的降水量分别作为阈值时所对应的降水量分别作为阈值R1、R2，而，而R0为月同期平均降为月同期平均降水量水量当以当以K指数为干旱指标时，将所有指数为干旱指标时，将所有K指数序列中（从小到大进行排序）百分位指数序列中（从小到大进行排序）百分位数为数为50%、35%和和25%时所对应的指数值作为阈值时所对应的指数值作为阈值

32、R0、R1和和R2当以当以P/E指数为干旱指标时，将单月指数序列中（从小到大进行排序）百分指数为干旱指标时，将单月指数序列中（从小到大进行排序）百分位数为位数为50%、35%和和25%和时所对应的指数值作为阈值和时所对应的指数值作为阈值R0、R1和和R2 EPK/根据邵阳市罗家庙站根据邵阳市罗家庙站19562002年的月降水量、月蒸发量，得年的月降水量、月蒸发量，得Pa指数、指数、 K指数和指数和P/E指数。用干旱过程识别方法，确定在指数。用干旱过程识别方法，确定在Pa指数、指数、K指数和指数和P/E指数指数下的干旱次数及干旱特征变量，见下表下图下的干旱次数及干旱特征变量，见下表下图 ：三、干

33、旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （二）基于降水和蒸发指标的干旱频率分析（二）基于降水和蒸发指标的干旱频率分析 干旱指数干旱指数干旱次数干旱次数干旱间隔均值干旱间隔均值(月月)干旱历时均值干旱历时均值(月月)干旱历时最大值干旱历时最大值(月月)干旱历时与烈度的干旱历时与烈度的Kendall相关系数相关系数Pa906.221.9090.54K876.482.0380.59P/E876.482.0390.51不同干旱指数下各年份的干旱月数不同干旱指数下各年份的干旱月数 由上表和上图可知：由上表和上图可知：在在47年共年共564个月中，以个月中，以Pa指数、指数

34、、K指数和指数和P/E指数为干旱指标识别得到指数为干旱指标识别得到的干旱月数分别为的干旱月数分别为171、177和和177，且，且3种指数下的干旱次数、干旱历时均种指数下的干旱次数、干旱历时均值和干旱历时最大值都较为接近值和干旱历时最大值都较为接近有有483个月为在个月为在3种指数下具有相同的识别结果（同为干旱或无旱），一种指数下具有相同的识别结果（同为干旱或无旱），一致性率达致性率达85.64%Pa指数与指数与K指数、指数、K指数与指数与P/E指数、指数、Pa指数与指数与P/E指数的识别一致性率指数的识别一致性率分别为分别为90.07%、89.36%、91.84%可见，可见，3种干旱指数对研

35、究区域的干旱指示程度具有较高的一致性，综种干旱指数对研究区域的干旱指示程度具有较高的一致性，综合考虑三者的干旱频率分析结果，可增强区域干旱频率分析的合理性和可合考虑三者的干旱频率分析结果，可增强区域干旱频率分析的合理性和可靠性靠性 三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （二）基于降水和蒸发指标的干旱频率分析（二）基于降水和蒸发指标的干旱频率分析 3种干旱指数在不同干旱历时和干旱烈度时的干旱重现期估计种干旱指数在不同干旱历时和干旱烈度时的干旱重现期估计见下图和下表，见下图和下表，在干旱历时分别为在干旱历时分别为2个月、个月、3个月和个月和4个月的平均干旱水

36、平下，个月的平均干旱水平下，K指数和指数和P/E指指数的相应于烈度均值的干旱重现期的计算结果相对接近，数的相应于烈度均值的干旱重现期的计算结果相对接近，Pa指数的相应结果指数的相应结果略微偏大略微偏大三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （二）基于降水和蒸发指标的干旱频率分析（二）基于降水和蒸发指标的干旱频率分析 Pa指数 K指数 P/E指数 干旱历时干旱历时(月月)Pa指数指数K指数指数P/E指数指数烈度均值烈度均值重现期重现期/年年烈度均值烈度均值重现期重现期/年年烈度均值烈度均值重现期重现期/年年245.962.760.602.130.932.46

37、384.484.971.184.481.904.894119.168.751.467.542.317.78以以3种指数为干旱指标计算了罗家庙站典型干旱年份的重现期，见下表。由种指数为干旱指标计算了罗家庙站典型干旱年份的重现期，见下表。由该表可知：该表可知：3种指数对种指数对1956、1960、1985和和1998年的干旱重现期计算结果相近年的干旱重现期计算结果相近，其中，其中，1956年的干旱重现期超过了年的干旱重现期超过了10 a，而，而1998年超过了年超过了60 a，干旱程度十分严重，干旱程度十分严重对于对于1963年，年，K指数和指数和P/E指数的干旱重现期计算结果相近，而指数的干旱重

38、现期计算结果相近，而Pa指数相指数相对偏小。对偏小。通过综合比较可知，通过综合比较可知，K指数和指数和P/E指数具有相对一致的干旱频率分指数具有相对一致的干旱频率分析结果，由于它们包含了降水量和蒸发量两个方面的信息，可进一步提高析结果，由于它们包含了降水量和蒸发量两个方面的信息，可进一步提高区域干旱区域干旱频率分析的可靠性频率分析的可靠性 三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （二）基于降水和蒸发指标的干旱频率分析（二）基于降水和蒸发指标的干旱频率分析 年份年份干旱重现期干旱重现期/年年PaKP/E195612.6114.0110.3019608.826

39、.756.39196325.5936.8140.8919856.005.857.70199868.2264.7965.26根据根据19562002年邵阳市茅坪站逐月径流资料，对该区域水文干旱频率进行年邵阳市茅坪站逐月径流资料，对该区域水文干旱频率进行分析研究。选取径流量距平百分率法（记为分析研究。选取径流量距平百分率法（记为A识别法）和径流量累积频率法识别法）和径流量累积频率法（记为（记为B识别法）两种阈值标准来识别水文干旱过程：识别法）两种阈值标准来识别水文干旱过程：根据根据水文情报预报规范水文情报预报规范中的标准，选取阈值中的标准，选取阈值R0为单月同期平均径流量，为单月同期平均径流量，而

40、将单月径流量距平百分率为而将单月径流量距平百分率为- -20%和和- -40%时所对应的径流量分别作为阈值时所对应的径流量分别作为阈值R1、R2 径流量序列一般认为服从径流量序列一般认为服从P-型曲线。对于多年同一月份的径流量序列，型曲线。对于多年同一月份的径流量序列，当确定当确定p后即可得后即可得xp。在径流量频率分析的基础上，本专题选用径流量的阈。在径流量频率分析的基础上，本专题选用径流量的阈值为值为R0 = x55%，R1 = x65%，R2 = x70% 三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （三）基于径流指标的干旱频率分析（三）基于径流指标的干

41、旱频率分析 A识别法和识别法和B识别法的干旱特征统计结果识别法的干旱特征统计结果三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （三）基于径流指标的干旱频率分析（三）基于径流指标的干旱频率分析 阈阈值值方方法法干旱干旱次数次数干旱间干旱间隔均值隔均值(月月)干旱历干旱历时均值时均值(月月)干旱烈度干旱烈度均值均值(m3/s)月月干旱历时干旱历时最大值最大值(月月)干旱烈度最干旱烈度最大值大值(m3/s)月月干旱历时与烈干旱历时与烈度的度的Kendall相关系数相关系数A717.893.9444.6213217.710.67B698.102.7022.4113139

42、.790.52由上表可知：由上表可知： 在在A、B两种阈值标准下的干旱次数、干旱间隔平均值及两种阈值标准下的干旱次数、干旱间隔平均值及Kendall相关相关系数都较为接近，都出现了最大历时达系数都较为接近，都出现了最大历时达13个月的干旱事件个月的干旱事件 A中平均干旱历时、干旱烈度和干旱烈度最大值，均大于中平均干旱历时、干旱烈度和干旱烈度最大值，均大于B的相应值。的相应值。可见，在不同阈值标准下的水文干旱特性具有一定的差异性，而合理确可见，在不同阈值标准下的水文干旱特性具有一定的差异性，而合理确定干旱发生阈值是区域水文干旱频率研究的关键点定干旱发生阈值是区域水文干旱频率研究的关键点 三、干旱

43、发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （三）基于径流指标的干旱频率分析（三）基于径流指标的干旱频率分析 005.1005.1G HC o p u la 理 论 联 合 概 率经 验 联 合 概 率005.1005.1C la y to nC o p u la 理 论 联 合 概 率经 验 联 合 概 率005.1005.1F ra n kC o p u la 理 论 联 合 概 率经 验 联 合 概 率005.1005.1GH Copula 理 论 联 合 概 率经验联合概率005.1005.1Clayton Copula 理 论 联 合 概 率经验联合概率00

44、5.1005.1Frank Copula 理 论 联 合 概 率经验联合概率三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （三）基于径流指标的干旱频率分析（三）基于径流指标的干旱频率分析 Copula函数函数A识别法识别法B识别法识别法适线前适线前适线后适线后适线前适线前适线后适线后Gumbel-Hougaard（GH）0.05090.05120.07720.0723Clayton0.05710.06810.08000.0890Frank0.04930.05620.07360.0809由上表可知：由上表可知： A中适线调整前后的中适线调整前后的Copula拟合最

45、小误差分别为拟合最小误差分别为0.0493、0.0512，而，而B中中适线调整前的适线调整前的Copula拟合最小误差为拟合最小误差为0.0736，调整后为，调整后为0.0723。可见经适。可见经适线调整前后的线调整前后的Copula函数拟合效果相当，说明适线调整干旱特征变量的分函数拟合效果相当，说明适线调整干旱特征变量的分布曲线保证了布曲线保证了Copula函数计算的准确性函数计算的准确性 3种种Copula函数都能较好反映干旱历时与干旱烈度间的相关关系，其中函数都能较好反映干旱历时与干旱烈度间的相关关系，其中GH及及Frank Copula函数拟合效果相对更优函数拟合效果相对更优三、干旱发

46、生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （三）基于径流指标的干旱频率分析（三）基于径流指标的干旱频率分析 根据根据Copula函数的估计参数，选用函数的估计参数，选用Gumbel-Hougaard (GH) Copula函数，函数，分别对分别对A、B识别方法中各年份的干旱重现期进行求解，见下图和下表：识别方法中各年份的干旱重现期进行求解，见下图和下表： 2种阈值标准下的干旱重现期计算结果大体上接近，基本能反映受旱严种阈值标准下的干旱重现期计算结果大体上接近，基本能反映受旱严重年份的干旱特征重年份的干旱特征 1963、1985年的干旱重现期都超过了年的干旱重现期都超过

47、了40 a，受旱十分严重，属特旱，受旱十分严重，属特旱 年份（年）年份（年）实际受旱率（实际受旱率（%）干旱重现期（干旱重现期（a）A识别法识别法B识别法识别法195633.279.2512.13196029.056.657.61196338.0845.0745.17198548.3448.7449.50199827.1116.5813.82选择淮北平原最北部的砀山县作为研究区域，该区域内有周寨、苇子园、李选择淮北平原最北部的砀山县作为研究区域，该区域内有周寨、苇子园、李庄等共庄等共1111处地下水埋深观测点处地下水埋深观测点三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的

48、应用研究 （四）基于地下水指标的干旱频率分析（四）基于地下水指标的干旱频率分析 三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （四）基于地下水指标的干旱频率分析（四）基于地下水指标的干旱频率分析 自自19751975年以来，每年以来，每5d5d观测观测1 1次、每月观测次、每月观测6 6次地下水埋深，取次地下水埋深，取6 6次均值作为该月次均值作为该月得平均地下水埋深，取得平均地下水埋深，取1111处观测点的地下水埋深平均值作为研究区域地下水处观测点的地下水埋深平均值作为研究区域地下水埋深。埋深。1111处观测点，大多自处观测点，大多自19751975年开始观测

49、（个别站点观测年份稍迟于年开始观测（个别站点观测年份稍迟于19751975年，极年，极少数站点存在部分缺测数据）。这里取少数站点存在部分缺测数据）。这里取1975-20071975-2007年观测数据进行分析，共年观测数据进行分析，共3333年，年，396396（33331212）个月份）个月份地下水平均埋深数据。记第地下水平均埋深数据。记第i年第年第j月地下水埋深为月地下水埋深为hij（i = 133，j = 112）。）。 第第i（i = 1 n）年第）年第j（j = 1 12）月的地下水埋深的变化）月的地下水埋深的变化hi,j记为记为 ,1,1,1,12(111)(11;12)i ji

50、ji jiihhjhhhinj三、干旱发生频率及重现期理论的应用研究三、干旱发生频率及重现期理论的应用研究 （四）基于地下水指标的干旱频率分析（四）基于地下水指标的干旱频率分析 基于地下水埋深变化的干旱识别及干旱特征变量确定的示意图：基于地下水埋深变化的干旱识别及干旱特征变量确定的示意图：当第当第j月地下水埋深的变化值月地下水埋深的变化值hi,j大于大于R1j时，则初步判断此月为干旱，即图中有时，则初步判断此月为干旱，即图中有a、b、c和和d共共4个干旱过程个干旱过程对于历时只有对于历时只有1个时段的干旱（如个时段的干旱（如a，d），若），若hi,j大于大于R2j（如（如a），则此月最终被），