[发明专利]一种基于多源信息分步校正作物模型的区域作物估产方法

摘要

本发明公开了一种基于多源信息分步校正作物模型的区域作物估产方法，包括：分解作物模型,获取多源信息以及分步校正作物模型，依次校正生育期模块、LAI模块和产量模块；每个功能模块的参数均经历如下参数优化过程：初始参数空间由数十组参数组成，运行作物模型得到每一组参数对应的模拟结果，建立当前优化功能模块的结果和校正数据之间的代价函数，在当前参数空间的结果中搜索代价函数的全局最优解，若最优解不符合优化终止条件，则向全局最优解和局部最优解共同决定的方向更新参数空间中当前功能模块的参数数值，重新计算模拟结果并进行判断；若最优解符合优化终止条件，校正过程结束，输出对应功能模块的校正后最优化参数和最优化模拟结果。

主权项

1.一种基于多源信息分步校正作物模型的区域作物估产方法，包括如下步骤：S1:分解作物模型,将作物模型的初始参数按照控制功能模块的不同分为生育期参数、LAI参数和产量参数三组；S2：获取分步校正所述作物模型的多源信息，包括基础输入信息、遥感校正信息和产量校正信息，其中基础输入信息包括驱动作物模型运行的日值气象数据、土壤属性信息和待校准的初始参数，遥感校正信息包括由LAI时间序列反演得到的生育期以及校正LAI参数所需要的LAI时间序列；产量校正信息包括用于校正作物模型的模拟产量的地面产量观测数据；S3：分步校正所述作物模型，为每一个最小行政单位依次校正生育期模块、LAI模块和产量模块；每个功能模块的参数均经历如下参数优化过程：初始参数空间由数十组参数组成，运行作物模型得到每一组参数对应的模拟结果，建立当前优化功能模块的结果和校正数据之间的代价函数，在当前参数空间的结果中搜索代价函数的全局最优解，若最优解不符合优化终止条件，则按照如下公式(1)和(2)向全局最优解和局部最优解共同决定的方向更新参数空间中当前功能模块的参数数值，重新计算模拟结果并进行判断；若最优解符合优化终止条件，校正过程结束，输出对应功能模块的校正后最优化参数和最优化模拟结果，其中经过校正后的当前功能模块最优化参数将取代原有参数空间中的对应参数，以保障下一功能模块的校正准确性；v_i(t+1)＝w*v_i(t)+c₁*rand1*(p_l(t)‑x_i(t))+c₂*rand2*(p_g(t)‑x_i(t))     (1)x_i(t+1)＝x_i(t)+v_i(t)        (2)式中，p_l(t)和p_g(t)分别为第t次模拟结束后的局部最优解和全局最优解对应的参数组数值；x_i(t)和x_i(t+1)分别为某一组参数中第i个参数在第t次和第t+1次模拟过程中的数值；v_i(t)和v_i(t+1)分别为某一组参数中第i个参数在第t次和第t+1次模拟过程中的更新幅度，c₁和c₂为加速系数，常取值为2；rand1和rand2为[0,1]之间随机数；w为惯性权重因子；其中，上述代价函数的计算优化终止条件为以下之一：连续5次循环后待优化模型参数值收缩到指定的值域范围；代价函数值在5次循环后无法提高0.01％；以及计算代价函数的次数超过100次。