[发明专利]一种干旱检测方法

权利要求书

1.一种干旱检测方法，其特征在于，所述方法包括以下步骤：

步骤一，获取某地区内降水量与蒸发量的差值数据，形成在一时间段内的时间序列数据D；

步骤二，利用平滑样条拟合的趋势值来代替Log-logistic分布的位置参数，得到时变的位置参数

步骤三，基于时变Dt时间序列的Log-logistic分布，用基于经验频率的概率加权距法PWMs来估计参数，进行干旱检测指标标准化；

所述步骤二具体为：

步骤2.1，使用平滑样条函数Smoothing Splines确定所述间序列数据D的线性或非线性趋势拟合；

步骤2.2，计算非平稳标准化降水指数NSPEI：

D是随着时间变化的降水与蒸散发的差值，假设Dt为时间序列D存在一个随着时间变化的自变量，

D_t＝P_t-PET_t (5)

利用平滑样条函数对Dt进行拟合：

其中，f为时间序列Dt的线性拟合函数；h为一个平滑参数；t＝1,…,n，为时间；PET采用Penman-Monteith蒸散发公式计算的潜在蒸发量，使用的气象指标分别为降水P、最高气温Tmax、最低气温Tmin、湿度H、风速W、太阳辐射S、海拔、纬度；

求解得到一个时变的位置参数如下：

2.根据权利要求1所述的一种干旱检测方法，其特征在于，所述步骤三具体为：基于时变Dt时间序列的Log-logistic分布为：

其中，α、β、分别是Log-logistic分布函数的尺度、形状、位置参数，用基于经验频率的概率加权距法PWMs来估计参数，

正的NSPEI值表示湿润，负的NSPEI表示干旱。

3.根据权利要求1所述的一种干旱检测方法，其特征在于，所述步骤二中，利用Kolmogorov-Smirnov K-S检验来评估残差是否符合Log-logistic分布，若K-S检验通过，说明残差服从Log-logistic分布，则说明可以对时间序列进行均值变化建模。