[发明专利]一种二维干旱风险估计方法及其系统

权利要求书

1.一种二维干旱风险估计方法，其特征在于，包括：

步骤S1：提取不同时期的干旱历时特征值和干旱烈度特征值；

步骤S2：分别基于干旱历时特征值和干旱烈度特征值拟合出干旱历时概率分布函数F_D(d)和干旱烈度概率分布函数F_S(s)，并构建干旱历时和干旱烈度的联合分布函数F(d,s)＝C(F_D(d),F_S(s))，通过最可能组合法推算重现期T_or的干旱历时d*和烈度s*；

步骤S3：基于联合分布函数F(d,s)，采用随机模拟技术得到B组新样本(d_ij,s_ij)，其中，i＝1,...,n，j＝1,...,B，d_ij和s_ij分别为第j次抽样得到的第i个样本的干旱历时和干旱烈度；

步骤S4：针对每组新样本，分别对干旱历时序列和干旱烈度序列进行拟合，得到第j组样本的干旱历时概率分布函数F_Dj(d)以及干旱烈度概率分布函数F_Sj(s)，并构建出第j组样本的干旱历时和干旱烈度的联合分布函数F_j(d,s)＝C(F_Dj(d),F_Sj(s))；

步骤S5：基于每组联合分布函数F_j(d,s)，通过最可能组合法推算重现期T_or的干旱历时d_Tj和烈度s_Tj；

步骤S6：基于B组干旱历时d_Tj和干旱烈度s_Tj，评估干旱历时d*和烈度s*的不确定性。

2.如权利要求1所述的二维干旱风险估计方法，其特征在于，步骤S6包括：采用置信区域的面积、横向平均偏移幅度、纵向平均偏移幅度、平均欧式距离中的任意一种或几种方式评估干旱历时d*和烈度s*的不确定性。

3.如权利要求2所述的二维干旱风险估计方法，其特征在于，横向平均偏移幅度D_Q、纵向平均偏移幅度D_W和平均欧式距离h的计算公式如下：

4.如权利要求1所述的二维干旱风险估计方法，其特征在于，B大于等于8000。

5.如权利要求1所述的二维干旱风险估计方法，其特征在于，在步骤S1中，获取干旱历时特征值和干旱烈度特征值的步骤包括：

步骤S1-1：获取长系列的栅格定量降水数据；

步骤S1-2：基于长系列的栅格定量降水数据和大气数据集，计算各栅格的降水蒸散指数；

步骤S1-3：基于降水蒸散指数和干旱阈值，提取各栅格长系列的干旱历时特征值和干旱烈度特征值。