反映非一致性干旱特征变量的干旱预警方法和装置与流程
本发明涉及干旱预警技术领域，尤其涉及一种反映非一致性干旱特征变量的干旱预警方法和装置。
背景技术：
在相关技术中，采用气象干旱指数和游程理论确定干旱特征变量，并根据干旱变量频率分析方法确定干旱重现期等。
其中，气象干旱指数包括降水量距平百分率、相对湿润度指数、标准化降水指数、土壤相对湿度干旱指数、帕默尔干旱指数、气象干旱指数等干旱指标值；另外，游程理论中，预先设置干旱指标阈值，若干旱指标值小于该干旱指标阈值则认为出现干旱，大于该干旱指标阈值则认为是正常或者洪涝的情况。进而根据干旱指标值与干旱指标阈值的大小关系，可以确定干旱事件及该干旱事件的干旱历时、干旱烈度、干旱强度、干旱影响面积等干旱特征变量，以便利用干旱变量频率分析方法揭示干旱的时空变化特征。
在相关技术中，各个干旱特征变量之间是独立同分布的，必须满足独立性、平稳性、一致性等要求，但是，在现实环境中，各个干旱特征变量将受到气候变化的影响而跟随时间发生非一致性的改变，例如：z国在1961-2013年中，平均年雨日呈显著减少趋势，每10年减少3.9天。2000-2013年，z国极端高温范围平均达626县，占z国总面积的27.4％，是常年的2.2倍等。由此可知，相关技术中仅对一致性干旱特征变量进行干旱分析的干旱频率分析方法存在与非一致性的干旱特征变量不匹配的缺陷，从而造成该方法得出的干旱重现期不准确。
由此可知，相关技术中采用干旱频率分析方法对非一致性干旱特征变量进行干旱分析方法存在分析结果不准确的缺陷。
技术实现要素：
本本发明实施例提供一种反映非一致性干旱特征变量的干旱预警方法和装置，以解决相关技术中的采用干旱频率分析方法对非一致性干旱特征变量进行干旱分析方法存在的分析结果不准确的问题。
为解决以上技术问题，本发明采用如下技术方案：
第一方面，本发明实施例提供了一种反映非一致性干旱特征变量的干旱预警方法，包括：
获取目标区域在预设历史时间段内的气象干旱指数；
根据所述气象干旱指数，确定干旱特征变量序列，其中，所述干旱特征变量序列包括：所述预设历史时间段内每一年的干旱特征变量；
根据目标概率分布函数，确定所述目标区域与预估干旱特征变量对应的干旱重现期，其中，所述目标概率分布函数与所述干旱特征变量序列匹配；
在所述干旱重现期大于或者等于预设重现期的情况下，输出干旱预警信息。
可选的，所述干旱特征变量序列包括：最长干旱历时序列或者最大干旱强度绝对值序列，其中，所述最长干旱历时序列包括所述预设历史时间段内每一年的最长干旱历时，所述最大干旱强度绝对值序列包括所述预设历史时间段内每一年的最大干旱强度的绝对值。
可选的，所述最长干旱历时序列通过以下公式计算得出：
其中，所述
所述i为1至n中的任一整数，所述n为所述预设历史时间段包括的年度的数量；
所述k为1至si中的任一整数，所述si为所述第i年的干旱事件的总数；
所述
所述
所述最大干旱强度绝对值序列通过以下公式计算得出：
其中，所述
所述
可选的，在所述根据目标概率分布函数，确定所述目标区域与预估干旱特征变量对应的干旱重现期之前，所述方法还包括：
根据所述干旱特征变量序列，确定预设备选概率分布集中每一个备选概率分布函数的时变参数，其中，所述时变参数包括以下至少一项：位置参数、尺度参数和形状参数，且所述时变参数为时间的三次样条函数；
从所述预设备选概率分布集中选取准确率最高的一个备选概率分布函数作为所述目标概率分布函数。
可选的，所述根据所述干旱特征变量序列，确定预设备选概率分布集中每一个备选概率分布函数的时变参数的步骤，包括：
将所述干旱特征变量序列中的干旱特征变量分别代入所述预设备选概率分布集中的每一个备选概率分布函数对应的对数似然函数
其中，所述f(xtθj(t))为所述预设备选概率分布集中的第j个备选概率分布函数在第t年内的概率密度函数；
所述θj(t)为所述预设备选概率分布集中的第j个备选概率分布函数在第t年内的时变参数；
所述j为小于或者等于所述预设备选概率分布集中备选概率分布函数的数量的任一正整数；
所述xt为所述第t年内的干旱特征变量；
所述t为1至n中的任一整数；
所述n为所述预设历史时间段包括的年度的数量。
可选的，所述从所述预设备选概率分布集中选取准确率最高的一个备选概率分布函数作为所述目标概率分布函数的步骤，包括：
根据公式：
确定与数值最小的评价指标对应的一个备选概率分布函数为所述目标概率分布函数。
可选的，所述根据目标概率分布函数，确定所述目标区域与预估干旱特征变量对应的干旱重现期的步骤，包括：
在所述干旱特征变量序列包括所述最长干旱历时序列的情况下，采用以下公式，计算得出所述目标区域与预估最长干旱历时对应的干旱重现期：
其中，所述预估干旱特征变量包括所述d，其表示所述预估最长干旱历时，所述预估最长干旱历时由预估的气象干旱指数确定；所述
或者，
在所述干旱特征变量序列包括所述最大干旱强度绝对值序列的情况下，采用以下公式，计算得出所述目标区域与预估最大干旱强度的绝对值对应的干旱重现期：
其中，所述预估干旱特征变量包括所述di，其表示所述预估最大干旱强度的绝对值，所述预估最大干旱强度的绝对值由所述预估的气象干旱指数确定；所述
第二方面，本发明实施例还提供一种反映非一致性干旱特征变量的干旱预警装置，包括：
获取模块，用于获取目标区域在预设历史时间段内的气象干旱指数；
第一确定模块，用于根据所述气象干旱指数，确定干旱特征变量序列，其中，所述干旱特征变量序列包括：所述预设历史时间段内每一年的干旱特征变量；
第二确定模块，用于根据目标概率分布函数，确定所述目标区域与预估干旱特征变量对应的干旱重现期，其中，所述目标概率分布函数与所述干旱特征变量序列匹配；
输出模块，用于在所述干旱重现期大于或者等于预设重现期的情况下，输出干旱预警信息。
第三方面，本发明实施例还提供一种电子设备，包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现本发明实施例提供的所述反映非一致性干旱特征变量的干旱预警方法中的步骤。
第四方面，本发明实施例还提供一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现本发明实施例提供的所述反映非一致性干旱特征变量的干旱预警方法中的步骤。
在本发明实施例中，获取目标区域在预设历史时间段内的气象干旱指数；根据所述气象干旱指数，确定干旱特征变量序列，其中，所述干旱特征变量序列包括：所述预设历史时间段内每一年的干旱特征变量；根据目标概率分布函数，确定所述目标区域与预估干旱特征变量对应的干旱重现期，其中，所述目标概率分布函数与所述干旱特征变量序列匹配；在所述干旱重现期大于或者等于预设重现期的情况下，输出干旱预警信息。其中，所述干旱特征变量序列中包括随时间变化的干旱特征变量，所述目标概率分布函数与所述干旱特征变量序列匹配使所述目标概率分布函数也跟随干旱特征变量的变化趋势而发生变化，从而使得所述目标概率分布函数与非一致性的干旱特征变量适配，进而提升了所述目标概率分布函数计算出的干旱重现期的准确率。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种反映非一致性干旱特征变量的干旱预警方法的流程图；
图2是本发明实施例提供的一种反映非一致性干旱特征变量的干旱预警装置的结构图之一；
图3是本发明实施例提供的一种反映非一致性干旱特征变量的干旱预警装置的结构图之二；
图4是本发明实施例提供的一种反映非一致性干旱特征变量的干旱预警装置的结构图之三；
图5是本发明实施例提供的一种电子设备的结构图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
请参阅图1，本发明实施例提供一种反映非一致性干旱特征变量的干旱预警方法，该方法可以包括以下步骤：
步骤101、获取目标区域在预设历史时间段内的气象干旱指数。
在具体实施中，所述气象干旱指数还可以称之为“综合气象干旱指数”、“复合气象干旱指数”或者气象干旱ci(compositeindex，综合)指数等。在实际应用中，所述预设历史时间段内的气象干旱指数可以包括：所述目标区域在预设历史时间段内的逐日气象干旱指数，例如：以cii,j表示所述目标区域在预设历史时间段内的逐日气象干旱指数序列，则所述i可以是年份号，所述j可以是第i年内的第j日，例如：ci2016,200表示2016年中第200天的气象干旱指数。
在具体实施中，若所述cii,j≤-0.6，则可以判断为目标区域在该天内出现了气象干旱旱情，若连续10天的cii,j均小于或者等于﹣0.6，则可以判断该目标区域发生了干旱事件，另外，若在发生干旱事件之后连续10天之内所述气象干旱ci指数均大于﹣0.6，则表示所述干旱事件结束。其中，所述目标区域和所述预设历史时间段可以根据用户的需求设置，例如：获取城市a在2000年至2016年内的逐日气象干旱指数。在此并不限定所述目标区域的大小和所述预设历史时间段的具体时间。
步骤102、根据所述气象干旱指数，确定干旱特征变量序列，其中，所述干旱特征变量序列包括：所述预设历史时间段内每一年的干旱特征变量。
在具体实施中，所述干旱特征变量序列可以包括：最长干旱历时序列或者最大干旱强度绝对值序列，其中，所述最长干旱历时序列包括所述预设历史时间段内每一年内的最长干旱历时，所述最大干旱强度绝对值序列包括所述预设历史时间段内每一年内的最大干旱强度的绝对值。
需要说明的是，所述干旱特征变量还可以包括除了所述最长干旱历时和所述最大干旱强度的绝对值之外任意其他干旱特征变量，在此不作具体限定。
在具体实施中，根据所述气象干旱指数，确定年最长干旱历时序列和年最大干旱强度绝对值序列，其中，以所述年最长干旱历时序列为例，其包括所述历史时间段内按照时间顺序排列的每一年的最长干旱历时。
步骤103、根据目标概率分布函数，确定所述目标区域与预估干旱特征变量对应的干旱重现期，其中，所述目标概率分布函数与所述干旱特征变量序列匹配。
在具体实施中，可以对目标区域在未来时间段内的气象干旱指数进行预测或者估计，以得到该目标区域的预估气象干旱指数，从而根据该预估气象干旱指数确定预估干旱特征变量。
需要说明的是，所述目标概率分布函数与所述干旱特征变量序列匹配，且所述干旱特征变量序列为跟随时间发生变化的非一致性序列，从而使得所述目标概率分布函数也跟随时间发生改变，在具体应用中，可以通过在所述目标概率分布函数中设置时变参数，以实现所述目标概率分布函数跟随时间发生改变。其中，所述目标概率分布函数的时变参数可以包括以下至少一项：形状参数、尺度参数和位置参数，且所述时变参数均为时间的三次样条函数。在具体实施中，可以将所述干旱特征变量序列依次输入至所述目标概率分布函数，并根据其输出结果与实际干旱情况的对比确定所述时变参数在各年内的取值。
另外，所述预估干旱特征变量对应的干旱重现期与所述干旱特征变量序列的种类、所述预估干旱特征变量的取值相关，且所述干旱特征变量序列的种类与所述预估干旱特征变量的种类相同，例如：所述干旱特征变量序列的种类为年最大干旱强度绝对值序列，则所述预估干旱特征变量也是一个预估年最大干旱强度的绝对值，向与所述年最大干旱强度绝对值序列匹配的目标概率分布函数输入所述预估年最大干旱强度的绝对值后计算得出的结果便为与所述预估年最大干旱强度的绝对值对应的干旱重现期。
步骤104、在所述干旱重现期大于或者等于预设重现期的情况下，输出干旱预警信息。
在具体实施中，所述预设重现期与所述干旱特征变量序列的种类、所述预估干旱特征变量的取值相关，例如：若所述干旱特征变量序列的种类为年最长干旱历时序列，所述预估干旱特征变量为年最长干旱历时等于40天，则对应的预设重现期可以设置为10年，表示将所述年最长干旱历时等于40天的预估干旱特征变量输入至所述目标概率分布函数时，若计算出的干旱重现期小于或者等于10年则输出干旱预警信息；若所述干旱特征变量序列的种类为年最长干旱历时序列，所述预估干旱特征变量为年最长干旱历时等于60天，则对应的预设重现期可以设置为15年，表示将所述年最长干旱历时等于60天的预估干旱特征变量输入至所述目标概率分布函数时，若计算出的干旱重现期小于或者等于15年则输出干旱预警信息。
当然，在实际应用中，还可以根据用户的需求、目标区域内蓄水设施的需水量、目标区域内土壤的特征等确定所述预设重现期。
在具体实施中，可以将所述目标概率分布函数对应的程序存储与计算机中，并根据输入的预估干旱特征变量计算出相应的干旱重现期，在实际应用中，所述输出干旱预警信息可以是发出声音告警信息或者输出干旱风险报告等，其中，所述声音告警信息或者干旱风险报告中可以包括所述目标地区发生预估干旱特征变量对应的干旱重现期小于预设重现期等信息。另外，还可以根据计算出的干旱重现期提供预防措施，以将所述目标区域内发生所述预估干旱特征变量对应的干旱事件的干旱重现期延长至大于所述预设重现期。
在本发明实施例中，获取目标区域在预设历史时间段内的气象干旱指数；根据所述气象干旱指数，确定干旱特征变量序列，其中，所述干旱特征变量序列包括：所述预设历史时间段内每一年的干旱特征变量；根据目标概率分布函数，确定所述目标区域与预估干旱特征变量对应的干旱重现期，其中，所述目标概率分布函数与所述干旱特征变量序列匹配；在所述干旱重现期大于或者等于预设重现期的情况下，输出干旱预警信息。其中，所述干旱特征变量序列中包括随时间变化的干旱特征变量，所述目标概率分布函数与所述干旱特征变量序列匹配使所述目标概率分布函数也跟随干旱特征变量的变化趋势而发生变化，从而使得所述目标概率分布函数与非一致性的干旱特征变量适配，进而提升了所述目标概率分布函数计算出的干旱重现期的准确率。
作为一种可选的实施方式，所述最长干旱历时序列通过以下公式计算得出：
其中，所述
所述i为1至n中的任一整数，所述n为所述预设历史时间段包括的年度的数量；
所述k为1至si中的任一整数，所述si为所述第i年内的干旱事件的总数；
所述
所述
所述最大干旱强度绝对值序列通过以下公式计算得出：
其中，所述
所述
在实际应用中，可以取干旱事件发生的日期作为所述
在实际应用中，所述
需要说明的是，如果某一次干旱事件跨越两个连续的日历年，且其干旱历时或干旱强度较强，则依据其主要发生的时间和前后两段的干旱强度将此干旱过程的指标归入前一年或后一年中。例如：假设一次干旱事件发生在12月15日至1月5日，12月15日至12月31日期间的干旱强度为-20,1月1日至1月5日期间的干旱强度为-4，由于干旱时间主要发生在前一年且前一年的干旱过程较强，所以计算时可将此次干旱事件对应的干旱特征变量纳入前一年中考虑。
另外，在计算所述最大干旱强度绝对值序列和所述最长干旱历时序列的过程中，需要分别针对所述预设历史时间段内的每一年分别进行计算，以得出所述预设历史时间段内的每一年的最大干旱强度的绝对值或者最长干旱历时。
作为一种可选的实施方式，在所述根据目标概率分布函数，确定所述目标区域与预估干旱特征变量对应的干旱重现期之前，所述方法还包括：
根据所述干旱特征变量序列，确定预设备选概率分布集中每一个备选概率分布函数的时变参数，其中，所述时变参数包括以下至少一项：位置参数、尺度参数和形状参数，且所述时变参数为时间的三次样条函数；
从所述预设备选概率分布集中选取准确率最高的一个备选概率分布函数作为所述目标概率分布函数。
在具体实施中，所述预设备选概率分布集中的备选概率分布函数可以包括韦伯分布函数、耿贝尔分布函数、弗雷歇分布函数、对数正态分布函数、伽玛分布函数等，在此并不穷举。在具体实施中，可以用mj(j＝1,2,…,l)表示所述预设备选概率分布集，用xt(t＝1,2,…,n)表示所述干旱特征变量序列，其中，所述l为所述预设备选概率分布集中包含的备选概率分布函数的数量，所述j可以取1至l之间的任一整数，所述n表示所述干旱特征变量序列中包含的干旱特征变量的数量，即所述预设历史时间段内包含的年度的数量，且所述t可以取1至n之间的任一整数，则与所述xt(t＝1,2,…,n)和第j个备选概率分布函数对应的概率密度函数可以表示为f(xtθj(t))，所述θj(t)为所述第j个备选概率分布函数的时变参数，其可以是时间t的三次样条函数，且所述时间t的三次样条函数中的可以包括待确定的参数。
需要说明的是，本实施方式中，所述i与所述t的取值相同，且均表示所述预设历史时间段内的第i年或者第t年。
在具体实施中，所述从所述预设备选概率分布集中选取准确率最高的一个备选概率分布函数作为所述目标概率分布函数，可以理解为：从所述预设备选概率分布集中选取能够良好的拟合所述干旱特征变量序列的一个备选概率分布函数作为所述目标概率分布函数。
本实施方式中，可以根据干旱特征变量序列确定预设备选概率分布集中每一个备选概率分布函数的时变参数，并在确定时变参数的备选概率分布函数中选取一个准确率最高的作为所述目标概率分布函数，从而使得所述目标概率分布函数能够更加准确的描述非一致性的干旱特征变量序列，从而使计算出的结果更加准确。
进一步的，所述根据所述干旱特征变量序列，确定预设备选概率分布集中每一个备选概率分布函数的时变参数的步骤，包括：
将所述干旱特征变量序列中的干旱特征变量分别代入所述预设备选概率分布集中的每一个备选概率分布函数对应的对数似然函数
其中，所述f(xtθj(t))为所述预设备选概率分布集中的第j个备选概率分布函数在第t年内的概率密度函数；
所述θj(t)为所述预设备选概率分布集中的第j个备选概率分布函数在第t年内的时变参数；
所述j为小于或者等于所述预设备选概率分布集中备选概率分布函数的数量的任一正整数；
所述xt为所述第t年内的干旱特征变量；
所述t为1至n中的任一整数；
所述n为所述预设历史时间段包括的年度的数量。
在具体实施中，所述θj(t)与所述
本实施方式中，提供计算时变参数的公式，并在
更进一步的，所述从所述预设备选概率分布集中选取准确率最高的一个备选概率分布函数作为所述目标概率分布函数的步骤，包括：
根据公式：
确定与数值最小的评价指标对应的一个备选概率分布函数为所述目标概率分布函数。
在具体实施中，所述p可以根据用户的需求预设值，所述dfj与所述预设备选概率分布集中第j个备选概率分布函数关联，从而分别将所述预设备选概率分布集中各个备选概率分布函数对应的参数代入公式
需要说明的是，在所述xt(t＝1,2,…,n)表示最大干旱强度绝对值序列的情况下，由于最大干旱强度的值小于0，在从预设备选概率分布函数中确定所述目标概率分布函数的过程中可以取所述最大干旱强度绝对值序列中的最大干旱强度的绝对值xt等于最大干旱强度的值的负值，以确保所述xt大于0。
本实施方式中，分别计算出所述预设备选概率分布集中每一个备选概率分布函数的评价指标eij，并通过所述评价指标eij来评价各个备选概率分布函数的误差率，从而选取评价指标eij最小的一个备选概率分布函数作为所述目标概率分布函数，可以确保所述目标概率分布函数的准确度为所述预设备选概率分布集中最高的一个，从而提升计算出的干旱重现期的准确性。
作为一种可选的实施方式，所述根据目标概率分布函数，确定所述目标区域与预估干旱特征变量对应的干旱重现期的步骤，包括：
在所述干旱特征变量序列包括所述最长干旱历时序列的情况下，采用以下公式，计算得出所述目标区域与预估最长干旱历时对应的干旱重现期：
其中，所述预估干旱特征变量包括所述d，其表示所述预估最长干旱历时，所述预估最长干旱历时由预估的气象干旱指数确定；所述
或者，
在所述干旱特征变量序列包括所述最大干旱强度绝对值序列的情况下，采用以下公式，计算得出所述目标区域与预估最大干旱强度的绝对值对应的干旱重现期：
其中，所述预估干旱特征变量包括所述di，其表示所述预估最大干旱强度的绝对值，所述预估最大干旱强度的绝对值由所述预估的气象干旱指数确定；所述
在具体实施中，在所述干旱特征变量序列为最大干旱强度绝对值序列的情况下，所述预设干旱重现期与预估最大干旱强度的绝对值对应，从而可以将所述预估最大干旱强度的绝对值输入与所述最大干旱强度绝对值序列匹配的目标概率分布函数时，根据输出结果确定干旱重现期，并将与预设干旱重现期对应的预设干旱重现期与该干旱重现期进行比较，已根据比较结果确定干旱重现期是否大于预设的预设干旱重现期，从而确定是否输出干旱预警信息。
在实现时，预估最长干旱历时和预估最大干旱的强度的绝对值可以通过以下过程确定：
1)根据当前干旱过程的实际日监测气象资料和未来时期的日预估气象资料，计算当前干旱过程的逐日综合气象干旱指数，其中，所述当前干旱过程的逐日综合气象干旱指数包括已经监测的时间段内的逐日综合气象干旱指数和未来时间段内预估的逐日综合气象干旱指数；
2)根据当前干旱过程的逐日综合气象干旱指数，确定当前干旱过程的开始日，并预估当前干旱过程的结束日期，计算当前干旱过程的预估持续时间和预估干旱强度；
3)将当前干旱过程的预估持续时间作为所述d，并将当前干旱过程的预估干旱强度的绝对值分别作为所述di。
另外，在所述干旱特征变量序列为最长干旱历时序列的情况下，所述预设干旱重现期与预估最长干旱历时对应，其原理与所述干旱特征变量序列为最大干旱强度绝对值序列的情况相同，在此不再赘述。
本实施方式中，分别针对所述干旱特征变量序列为最大干旱强度绝对值序列和最长干旱历时序列两种情况，分别提供计算干旱重现期的方法，并在计算出干旱重现期之后，分别与预估最大干旱强度的绝对值和预估最长干旱历时对应的预设干旱重现期进行比较，以确定计算出的干旱重现期是否大于预设的预设干旱重现期，从而确定是否输出干旱预警信息。
请参阅图2，本发明实施例还提供一种反映非一致性干旱特征变量的干旱预警装置200，该装置200包括：
获取模块201，用于获取目标区域在预设历史时间段内的气象干旱指数。
第一确定模块202，用于根据所述气象干旱指数，确定干旱特征变量序列，其中，所述干旱特征变量序列包括：所述预设历史时间段内每一年的干旱特征变量。
第二确定模块203，用于根据目标概率分布函数，确定所述目标区域与预估干旱特征变量对应的干旱重现期，其中，所述目标概率分布函数与所述干旱特征变量序列匹配。
输出模块204，用于在所述干旱重现期大于或者等于预设重现期的情况下，输出干旱预警信息。
可选的，所述干旱特征变量序列包括：最长干旱历时序列或者最大干旱强度绝对值序列，其中，所述最长干旱历时序列包括所述预设历史时间段内每一年内的最长干旱历时，所述最大干旱强度绝对值序列包括所述预设历史时间段内每一年内的最大干旱强度的绝对值。
可选的，所述最长干旱历时序列通过以下公式计算得出：
其中，所述
所述i为1至n中的任一整数，所述n为所述预设历史时间段包括的年度的数量；
所述k为1至si中的任一整数，所述si为所述第i年内的干旱事件的总数；
所述
所述
所述最大干旱强度绝对值序列通过以下公式计算得出：
其中，所述
所述
可选的，如图3所示，装置200还包括：
第三确定模块205，用于在所述根据目标概率分布函数，确定所述目标区域与预估干旱特征变量对应的干旱重现期之前，根据所述干旱特征变量序列，确定预设备选概率分布集中每一个备选概率分布函数的时变参数，其中，所述时变参数包括以下至少一项：位置参数、尺度参数和形状参数，且所述时变参数为时间的三次样条函数；
选取模块206，用于从所述预设备选概率分布集中选取准确率最高的一个备选概率分布函数作为所述目标概率分布函数。
可选的，所述选取模块206具体用于：
将所述干旱特征变量序列中的干旱特征变量分别代入所述预设备选概率分布集中的每一个备选概率分布函数对应的对数似然函数
其中，所述f(xtθj(t))为所述预设备选概率分布集中的第j个备选概率分布函数在第t年内的概率密度函数；
所述θj(t)为所述预设备选概率分布集中的第j个备选概率分布函数在第t年内的时变参数；
所述j为小于或者等于所述预设备选概率分布集中备选概率分布函数的数量的任一正整数；
所述xt为所述第t年内的干旱特征变量；
所述t为1至n中的任一整数；
所述n为所述预设历史时间段包括的年度的数量。
可选的，如图4所示，所述选取模块206，包括：
计算单元2061，用于根据公式：
确定单元2062，用于确定与数值最小的评价指标对应的一个备选概率分布函数为所述目标概率分布函数。
可选的，所述根第二确定模块203具体用于：
在所述干旱特征变量序列包括所述最长干旱历时序列的情况下，采用以下公式，计算得出所述目标区域与预估最长干旱历时对应的干旱重现期：
其中，所述预估干旱特征变量包括所述d，其表示所述预估最长干旱历时，所述预估最长干旱历时由预估的气象干旱指数确定；所述
或者，
在所述干旱特征变量序列包括所述最大干旱强度绝对值序列的情况下，采用以下公式，计算得出所述目标区域与预估最大干旱强度的绝对值对应的干旱重现期：
其中，所述预估干旱特征变量包括所述di，其表示所述预估最大干旱强度的绝对值，所述预估最大干旱强度的绝对值由所述预估的气象干旱指数确定；所述
本发明实施例提供的反映非一致性干旱特征变量的干旱预警装置能够实现图1所示的方法实施例中反映非一致性干旱特征变量的干旱预警装置实现的各个过程，且能够取得相同的有益效果，为避免重复，这里不再赘述。
请参阅图5，本发明实施例还提供一种电子设备，包括：收发机501、处理器502、存储器503、总线接口以及存储在存储器503上并可在处理器501上运行的计算机程序。
收发机501用于获取目标区域在预设历史时间段内的气象干旱指数。
处理器402执行所述计算机程序，并根据所述气象干旱指数，确定干旱特征变量序列，其中，所述干旱特征变量序列包括：所述预设历史时间段内每一年的干旱特征变量。
处理器402还用于根据目标概率分布函数，确定所述目标区域与预估干旱特征变量对应的干旱重现期，其中，所述目标概率分布函数与所述干旱特征变量序列匹配；并在所述干旱重现期大于或者等于预设重现期的情况下，输出干旱预警信息。
可选的，所述干旱特征变量序列包括：最长干旱历时序列或者最大干旱强度绝对值序列，其中，所述最长干旱历时序列包括所述预设历史时间段内每一年内的最长干旱历时，所述最大干旱强度绝对值序列包括所述预设历史时间段内每一年内的最大干旱强度的绝对值。
可选的，所述最长干旱历时序列通过以下公式计算得出：
其中，所述
所述i为1至n中的任一整数，所述n为所述预设历史时间段包括的年度的数量；
所述k为1至si中的任一整数，所述si为所述第i年内的干旱事件的总数；
所述
所述
所述最大干旱强度绝对值序列通过以下公式计算得出：
其中，所述
所述
可选的，处理器402在所述根据目标概率分布函数，确定所述目标区域与预估干旱特征变量对应的干旱重现期之前，还用于：
根据所述干旱特征变量序列，确定预设备选概率分布集中每一个备选概率分布函数的时变参数，其中，所述时变参数包括以下至少一项：位置参数、尺度参数和形状参数，且所述时变参数为时间的三次样条函数；
从所述预设备选概率分布集中选取准确率最高的一个备选概率分布函数作为所述目标概率分布函数。
可选的，处理器402执行的所述根据所述干旱特征变量序列，确定预设备选概率分布集中每一个备选概率分布函数的时变参数的步骤，包括：
将所述干旱特征变量序列中的干旱特征变量分别代入所述预设备选概率分布集中的每一个备选概率分布函数对应的对数似然函数
其中，所述f(xtθj(t))为所述预设备选概率分布集中的第j个备选概率分布函数在第t年内的概率密度函数；
所述θj(t)为所述预设备选概率分布集中的第j个备选概率分布函数在第t年内的时变参数；
所述j为小于或者等于所述预设备选概率分布集中备选概率分布函数的数量的任一正整数；
所述xt为所述第t年内的干旱特征变量；
所述t为1至n中的任一整数；
所述n为所述预设历史时间段包括的年度的数量。
可选的，处理器402执行的所述从所述预设备选概率分布集中选取准确率最高的一个备选概率分布函数作为所述目标概率分布函数的步骤，包括：
根据公式：
确定与数值最小的评价指标对应的一个备选概率分布函数为所述目标概率分布函数。
可选的，处理器402所述根据目标概率分布函数，确定所述目标区域与预估干旱特征变量对应的干旱重现期的步骤，包括：
在所述干旱特征变量序列包括所述最长干旱历时序列的情况下，采用以下公式，计算得出所述目标区域与预估最长干旱历时对应的干旱重现期：
其中，所述预估干旱特征变量包括所述d，其表示所述预估最长干旱历时，所述预估最长干旱历时由预估的气象干旱指数确定；所述
或者，
在所述干旱特征变量序列包括所述最大干旱强度绝对值序列的情况下，采用以下公式，计算得出所述目标区域与预估最大干旱强度的绝对值对应的干旱重现期：
其中，所述预估干旱特征变量包括所述di，其表示所述预估最大干旱强度的绝对值，所述预估最大干旱强度的绝对值由所述预估的气象干旱指数确定；所述
本发明实施例中，处理器执行计算机程序时，实现如上所述的反映非一致性干旱特征变量的干旱预警方法的步骤。取得与上述方法实施例中相同的有益效果，为避免重复，在此不再赘述。
本发明实施例还提供一种计算机存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如上所述的反映非一致性干旱特征变量的干旱预警方法中的步骤，且能达到相同的技术效果，为避免重复，这里不再赘述。
在本申请所提供的几个实施例中，应该理解到，所揭露方法和装置，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理包括，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。
以上，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。
- 还没有人留言评论。精彩留言会获得点赞！