一种干旱重现期计算中确定干旱阈值的方法技术

本发明专利技术公开了一种干旱重现期计算中确定干旱阈值的方法，特别涉及一种以干旱事件最长调查期为约束条件，使频率分析计算的历史干旱事件对应的最大重现期等于历史干旱最大调查期来确定干旱阈值的方法；本发明专利技术克服了目前在干旱重现期计算过程中确定干旱阈值存在的主观性及其导致的计算的历史干旱事件重现期明显大于历史调查期这一不足；在干旱重现期计算方面，具有较强的应用前景。

【技术实现步骤摘要】

本专利技术属于干旱重现期计算
，特别涉及一种干旱重现期计算中确定干旱阈值的方法。
技术介绍
干旱是我国主要的自然灾害，其具有影响面积大、影响范围广、持续时间长、连片发生及连季/连年发生等特点，对农业生产、城乡饮用水安全、社会经济发展以及生态环境造成严重影响。干旱事件严重程度的度量通常采用干旱重现期指标。目前常采用Copula函数构建干旱历时和干旱烈度两变量多维联合分布函数计算干旱事件重现期。其首要步骤是干旱事件的识别，即首先需要确定干旱阈值，然后基于干旱阈值，从历史事件中筛选哪些事件是干旱事件，进而获得干旱历时和干旱烈度样本系列，以供推求干旱重现期之用。目前关于干旱阈值的确定方法，如供-需差异程度或过程模拟等，都存在一定的主观性，且计算的历史干旱事件的最大重现期与实际干旱调查期通常存在较大差异。合理地确定干旱阈值大小是保证干旱事件识别的合理性及干旱重现期计算的可靠性的关键。
技术实现思路
专利技术目的：为了克服目前在干旱重现期计算过程中确定干旱阈值存在的主观性及其导致的计算的历史干旱事件重现期明显大于历史调查期这一不合理现象，本专利技术提供一种干旱重现期计算中确定干旱阈值的方法。技术方案：为实现上述目的，本专利技术采用的技术方案为：一种干旱重现期计算中确定干旱阈值的方法，包括以下步骤：(1)最大调查期确定：收集研究对象与干旱相关的资料，根据历史干旱事件资料，确定干旱事件的最大调查期；(2)干旱指标计算：采用相对距平指标，对降雨数据资料进行标准化处理，获得用于评估干旱事件的干旱指标样本系列；(3)干旱阈值设定：针对研究对象特征，初步设定一组不同的干旱阈值水平；(4)干旱事件识别：根据给定的干旱阈值，识别对应的干旱事件序列，并计算每场干旱的干旱历时和干旱烈度，获得干旱历时和干旱烈度样本系列；(5)边缘分布函数拟合：分别选取合适类型的分布函数对干旱历时和干旱烈度进行拟合，并评估拟合精度；(6)联合分布函数构建：采用多维联合分布函数，构建干旱历时和干旱烈度的两维联合分布函数对干旱事件进行描述；(7)根据步骤(6)构建好的分布函数，计算历史干旱事件中每一场干旱事件的重现期，并确定最大重现期值；(8)重复执行步骤(4)～(7)，计算不同干旱阈值水平对应的最大干旱重现期，选取计算的最大干旱重现期等于或接近步骤(1)中干旱最大调查期的干旱阈值，作为识别干旱事件的最终阈值。进一步的，所述步骤(5)中采用马秀峰提出的游程理论构建干旱历时分布函数。进一步的，所述步骤(5)采用2参数的伽马分布函数对干旱烈度分布函数进行拟合。进一步的，所述步骤(6)中采用G-H Copula函数，构建干旱历时分布函数和干旱烈度分布函数的两维联合分布函数对干旱事件进行描述。进一步的，所述步骤(8)中，选取计算的最大干旱重现期等于或接近历史干旱最大调查期的干旱阈值，作为识别干旱事件的最终阈值。有益效果：与现有技术相比，本专利技术具有以下有益效果：本专利技术提供了一种干旱重现期计算中确定干旱阈值的方法，克服了目前在干旱重现期计算过程中确定干旱阈值存在的主观性及其导致的计算的历史干旱事件重现期明显大于历史调查期这一不足；在干旱重现期计算方面，具有较强的应用前景。具体实施方式下面结合实施例对本专利技术作更进一步的说明。一种干旱重现期计算中确定干旱阈值的方法，包括以下步骤：(1)最大调查期确定：收集研究对象与干旱相关的资料，根据历史干旱事件资料，确定干旱事件的最大调查期；(2)干旱指标计算：采用相对距平指标，对降雨数据资料进行标准化处理，获得用于评估干旱事件的干旱指标样本系列；(3)干旱阈值设定：针对研究对象特征，初步设定一组不同的干旱阈值水平；(4)干旱事件识别：根据给定的干旱阈值，识别对应的干旱事件序列，并计算每场干旱的干旱历时和干旱烈度，获得干旱历时和干旱烈度样本系列；(5)边缘分布函数拟合：分别选取合适类型的分布函数对干旱历时和干旱烈度进行拟合，并评估拟合精度；(6)联合分布函数构建：采用多维联合分布函数，构建干旱历时和干旱烈度的两维联合分布函数对干旱事件进行描述；(7)根据步骤(6)构建好的分布函数，计算历史干旱事件中每一场干旱事件的重现期，并确定最大重现期值；(8)重复执行步骤(4)～(7)，计算不同干旱阈值水平对应的最大干旱重现期，选取计算的最大干旱重现期等于或接近步骤(1)中干旱最大调查期的干旱阈值，作为识别干旱事件的最终阈值。所述步骤(5)中采用马秀峰提出的游程理论构建构建干旱历时样本系列的概率分布函数。所述步骤(5)采用2参数的伽马分布函数对干旱烈度样本系列的分布函数进行拟合。所述步骤(6)中采用G-H Copula函数，构建干旱历时分布函数和干旱烈度分布函数的两维联合分布函数对干旱事件进行描述。进一步的，所述步骤(8)中，选取计算的最大干旱重现期等于或接近历史干旱最大调查期的干旱阈值，作为识别干旱事件的最终阈值。实施例现有某一流域50年的月降雨资料；依据本专利技术方法，干旱重现期计算中确定干旱阈值的过程为：(1)最大限度地收集该研究流域与干旱相关的资料(如水文气象、历史干旱记录等)，若历史资料丰富详实，可判断最远的历史干旱事件发生年份时，将最远年份距离目前的年数作为最大历史调查期；如历史资料缺乏，不足以对实测干旱事件作历史干旱调查时，可将实测的干旱年数作为干旱调查期；(2)采用降雨相对距平法，对50年的月降雨资料p1,p2,…,p600进行统一化处理，获得干旱指标样本系列x1,x2,…,x600；(3)计算干旱指标系列的均值：(4)初步设定一组干旱阈值：k1＝0.9E,k2＝0.85E,…,k16＝0.15E,k17＝0.1E，共17个不同的干旱阈值；(5)根据给定的干旱阈值ki,i＝1,2,…17，识别其对应的干旱事件系列，即当干旱指标xi,i＝1,2,…,600小于干旱阈值ki时，即可判定此月为干旱月份；(6)根据干旱识别的游程理论，识别出Mi,i＝1,2,…17场干旱事件，计算每场干旱的干旱历时D和干旱烈度S，获得干旱历时和干旱烈度样本系列，记为D1,D2,…,DMi和S1,S2,…,SMi；(7)采用马秀峰提出的游程分布函数对干旱历时样本系列D1,D2,…,Dm的分布函数进行拟合；(8)采用2参数的伽马分布函数对干旱烈度样本系列S1,S2,…,Sm的分布函数进行拟合；(9)采用G-H Copula函数，构建干旱历时分布函数(步骤(7))和干旱烈度分布函数(步骤(8))的两维联合分布函数对干旱事件进行描述；(10)根据步骤(9)构建好的分布函数，计算Mi,i＝1,2,…17场历史干旱事件中每一场干旱事件的重现期并确定最大重现期值i＝1,2,…17；(11)通过步骤(5)～(10)，可获得不同干旱阈值水平ki,i＝1,2,…17下，计算的最大干旱重现期i＝1,2,…17，选取计算的最大干旱重现期等于或接近步骤(1)中干旱最大调查期的干旱阈值，作为识别干旱事件的最终阈值。本专利技术是一种干旱重现期计算中确定干旱阈值的方法，特别涉及一种以干旱事件最长调查期为约束条件，使频率分析计算的历史干旱事件对应的最大重现期等于历史干旱的最大调查期来确定干旱阈值的方法。本专利技术提供了，克服了目前在干旱重现期计算过程中确定干旱阈值存在的主观性及其导致的计算的历史干旱本文档来自技高网...

【技术保护点】
一种干旱重现期计算中确定干旱阈值的方法，其特征在于，包括以下步骤：（1）最大调查期确定：收集研究对象与干旱相关的资料，根据历史干旱事件资料，确定干旱事件的最大调查期；（2）干旱指标计算：采用相对距平指标，对降雨数据资料进行标准化处理，获得用于评估干旱事件的干旱指标样本系列；（3）干旱阈值设定：针对研究对象特征，初步设定一组不同的干旱阈值水平；（4）干旱事件识别：根据给定的干旱阈值，识别对应的干旱事件序列，并计算每场干旱的干旱历时和干旱烈度，获得干旱历时和干旱烈度样本系列；（5）边缘分布函数拟合：分别选取合适类型的分布函数对干旱历时和干旱烈度进行拟合，并评估拟合精度；（6）联合分布函数构建：采用多维联合分布函数，构建干旱历时和干旱烈度的两维联合分布函数对干旱事件进行描述；（7）根据步骤（6）构建好的分布函数，计算历史干旱事件中每一场干旱事件的重现期，并确定最大重现期值；（8）重复执行步骤（4）~（7），计算不同干旱阈值水平对应的最大干旱重现期，选取计算的最大干旱重现期等于或接近步骤（1）中干旱最大调查期的干旱阈值，作为识别干旱事件的最终阈值。

【技术特征摘要】
1.一种干旱重现期计算中确定干旱阈值的方法，其特征在于，包括以下步骤：（1）最大调查期确定：收集研究对象与干旱相关的资料，根据历史干旱事件资料，确定干旱事件的最大调查期；（2）干旱指标计算：采用相对距平指标，对降雨数据资料进行标准化处理，获得用于评估干旱事件的干旱指标样本系列；（3）干旱阈值设定：针对研究对象特征，初步设定一组不同的干旱阈值水平；（4）干旱事件识别：根据给定的干旱阈值，识别对应的干旱事件序列，并计算每场干旱的干旱历时和干旱烈度，获得干旱历时和干旱烈度样本系列；（5）边缘分布函数拟合：分别选取合适类型的分布函数对干旱历时和干旱烈度进行拟合，并评估拟合精度；（6）联合分布函数构建：采用多维联合分布函数，构建干旱历时和干旱烈度的两维联合分布函数对干旱事件进行描述；（7）根据步骤（6）构建好的分布函数，计算历史干旱事件中每一场干旱事件的重现期，并确定最大重现期值；（8）重复执行步骤（...

【专利技术属性】
技术研发人员：梁忠民，胡义明，王军，李彬权，王栋，蒋晓蕾，刘甜，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏;32