1.一种地下水储量异常反演产品适用性评价方法，其特征在于，包括以下步骤：步骤S1：根据GRACE/GRACE-FO总水储量异常产品及水文模型模拟地表水储量和土壤水储量异常产品反演生成区域地下水储量异常反演产品，计算像元尺度实测水位时间序列和地下水储量异常时间序列之间的互相关系数、相位差变异系数和幅度比变异系数，从趋势、相位和幅度多个视角评价地下水储量异常反演产品适用性；步骤S2：利用评价区域所有像元互相关系数、相位差变异系数和幅度比变异系数组成评价样本集，利用主成分分析法确定互相关系数、相位差变异系数和幅度比变异系数权重，加权计算像元尺度地下水储量异常反演产品综合适用性指数；步骤S3：根据地下水储量异常反演产品综合适用性指数结果，分级划分适用性等级，评价区域地下水储量异常反演产品用于地下水储量时空变化分析的适用性；所述步骤S1包括以下子步骤：子步骤S11：生产基于GRACE/GRACE-FO总水储量异常产品和GLDAS水文模型模拟地表水储量和土壤水储量异常产品反演生成区域地下水储量异常反演产品；子步骤S12：获取区域实测地下水位数据；子步骤S13：多视角评价基于实测地下水位的区域地下水储量异常反演产品适用性；所述子步骤S13的多视角适用性评价指标包括互相关系数、相位差变异系数和幅度比变异系数；所述互相关系数计算公式为：式中：σ x 、σ y 分别为时间序列x(t)和y(t)的均方差； 分别为x(t)和y(t)的均值；k是滞后量，表示两个时间序列的偏移程度；C xy (k)和r xy (k)分别为时间序列x(t)和y(t)在时滞k下的互协方差和互相关系数；n是时间序列的长度；互相关系数值域为-1到1，两个时间序列互相关系数越接近1，说明两个时间序列的变化趋势越相同；所述相位差变异系数的计算公式为：式中：σ PCE 为PCE的标准差；μ PCE 为PCE的平均值；PCV越接近0，说明两个时间序列相位一致性越高；所述幅度比变异系数的计算公式为：式中：σ AR 为AR的标准差；μ AR 为AR的平均值；ACV越接近0，说明两个时间序列幅度一致性越高。
2.根据权利要求1所述的地下水储量异常反演产品适用性评价方法，其特征在于，所述子步骤S11的地下水储量异常的表达式为：GWSA＝TWSA-SWSA-SMSA式中，GWSA表示地下水储量异常，TWSA表示总水储量异常，SWSA表示地表水储量异常，SMSA表示土壤水储量异常。
3.根据权利要求1所述的地下水储量异常反演产品适用性评价方法，其特征在于，所述步骤S2采用利用主成分分析法确定互相关系数、相位差变异系数和幅度比变异系数权重，加权计算像元尺度地下水储量异常反演产品综合适用性指数，所述地下水储量异常反演产品综合适用性指数，其计算公式如下：式中：PCV p 为第p个主成分；α p 为对应主成分的贡献率；对地下水储量异常反演产品综合适用性指数进行归一化处理，计算公式为：式中：GSEI norm 为归一化后的综合适用性指标值；GSEI、GSEI max 和GSEI min 分别表示地下水储量异常反演产品综合适用性指数的原始值、所有像元中的最大值和最小值。
4.根据权利要求1所述的地下水储量异常反演产品适用性评价方法，其特征在于，所述步骤S3的分级划分适用性等级，具体为基于地下水储量异常反演产品综合适用性指数结果，对适用性指数进行分级处理，以0.2为一个区间将划分为五个等级：一级数值区间为[0，0.2]，表示地下水储量异常反演产品的适用性“差”；二级数值区间为(0.2，0.4]，表示适用性“较差”；三级数值区间为(0.4，0.6]，表示适用性“一般”；四级数值区间为(0.6，0.8]，表示适用性“良好”；五级数值区间为(0.8，1]，表示适用性“优”，根据以上分级结果对地下水储量异常反演产品的适用性进行综合评价。
