一种估算缺资料地区参考作物蒸散量的方法

摘要

本发明公开了一种估算缺资料地区参考作物蒸散量的方法，包括资料收集，参数优化，地理因子选择，基于机器学习的参数区域化模型，计算得到的Hargreaves‑Samani模型参数，采用气温数据进一步计算缺资料地区参考作物蒸散量。本发明充分借助研究区内资料齐全地区的气象数据和地理因子信息，结合机器学习算法构建了Hargreaves‑Samani模型参数的缺资料区优化估计方法，相对于原始的Hargreaves‑Samani模型能显著提高缺资料地区参考作物蒸散量计算精度，同时相比于以往直接利用气象数据和机器学习模型建模的方式，本发明以Hargreaves‑Samani模型为基准，具备良好的物理基础，为缺资料地区精准的参考作物蒸散量估算提供了技术手段。

说明书

技术领域

本发明涉及农业水管理及水文水资源应用领域，具体涉及一种估算缺资料地区参考作物蒸散量的方法。

背景技术

蒸散量是气候学、农业科学和水文水资源领域的基本要素，也是农田精准灌溉决策和水资源管理的重要标准。然而在全球大多数地方仍然无法直接观测蒸散过程，蒸散量的计算精度主要取决于参考作物蒸散量，参考作物蒸散量代表大气蒸发能力的大小，通常指在生长茂盛，水分充足条件下，高约0.12m完全覆盖地面的草地(如苜蓿)的蒸散量。参考作物蒸散量是作物蒸散量计算的关键参照标准，其计算精度直接决定了实际蒸散量的估算结果，对于陆地水热循环模拟和精准的农业水管理工作至关重要。

目前农业领域常采用联合国粮农组织FAO推荐的Penman-Monteith(PM)模型计算参考作物蒸散量，但PM模型要求全面的最高和最低气温、湿度、日照时数(或辐射)、风速和大气压强众多气象数据，而全球多数地区难以提供PM模型所需的全部气象信息。另一方面，仅基于温度信息的参考作物蒸散量计算公式如Hargreaves-Samani模型被FAO 推荐为唯一的替补PM模型的方法，可用于仅有温度信息的区域参考作物蒸散量估算。但在全球范围内的应用表明，Hargreaves-Samani模型与PM模型方法计算结果相比差异较大，无法给出缺资料地区精准的参考作物蒸散量计算结果。

申请号为202010675336.3的中国专利提出了利用BP神经网络和降水信息预测参考蒸散发的方法，但该方法训练中所需要的参考值是以PM模型计算结果为准，这就限制了这种方法在缺资料地区的有效应用，因为大多数缺资料地区根本没有足够气象资料满足PM模型的计算，自然无法实现机器学习模型的训练。另外，这种方法完全依赖机器学习模型的训练，脱离了物理模型基础，这导致该方法在复杂地球物理模型的子过程建模中适用性不高。因此，在多数情况下，需要一种能够以物理模型为基础，可精准估计缺资料地区参考作物蒸散量的方法。

发明内容

本发明的目的在于以现有物理模型为基础，针对缺资料地区参考作物蒸散量计算难度大、缺乏物理模型支撑、精度不足的缺点，提供一种以物理模型为基础，可以更精确估算缺资料地区参考作物蒸散量的方法，借助机器学习充分挖掘有资料地区信息，以提高缺资料地区的参考作物蒸散量估算的准确性与合理性。

本发明提供了一种估算缺资料地区参考作物蒸散量的方法，包括以下步骤：

(1)资料收集：收集区域内有资料站点的最高和最低气温、湿度、日照时数(或辐射)、风速和大气压强气象要素资料，以及缺资料站点的最高和最低气温资料。

有资料区指具备长系列的最高和最低气温、湿度、日照时数(或辐射)、风速和大气压强气象观测信息；所述无资料区指缺乏长系列的气象观测信息，但具备多年平均降水量、多年平均最高气温、多年平均最低气温和多年平均参考作物蒸散量气候态信息，以及具备高程和站点纬度位置信息。

(2)参数优化：以Penman-Monteith模型计算的ET0作为参考值，以KGE系数作为目标函数，以KGE系数最大为优化目标，采用SCE-UA优化算法对有资料区站点的 Hargreaves-Samani模型参数进行优化，计算得到每个有资料区站点的Hargreaves-Samani 模型优化参数a

(3)地理因子选择：以方便获取为基本原则，选取每个气象站的地理因子，包括高程、多年平均降水量、多年平均最高气温、多年平均最低气温、多年平均潜在蒸散量和干旱指数(多年平均潜在蒸散量/多年平均降水量)。

(4)基于机器学习的参数区域化模型：根据有资料地区的气象站点地理因子和优化后的Hargreaves-Samani模型优化参数a

(5)基于步骤(4)中参数区域化方法模拟的Hargreaves-Samani模型参数，采用气温数据进一步计算缺资料地区参考作物蒸散量。

本发明的有益效果在于：

本发明充分借助研究区内资料齐全地区的气象数据和地理因子信息，结合机器学习算法构建了Hargreaves-Samani模型参数的缺资料区优化估计方法，相对于原始的Hargreaves-Samani模型能显著提高缺资料地区参考作物蒸散量计算精度，同时相比于以往直接利用气象数据和机器学习模型建模的方式，本发明以Hargreaves-Samani模型为基准，具备良好的物理基础，为缺资料地区精准的参考作物蒸散量估算提供了技术手段。

附图说明

图1为本发明方法的流程图；

图2为随机森林机器学习参数区域化模拟参数值和优化参数值的散点图，其中，2a为参数a

图3为本发明计算的某一站点区域化模拟得到的参考蒸散量与PM模型计算的参考蒸散量值对比；

图4为Hargreaves-Samani模型原始参数计算的某一站点参考蒸散量与PM模型计算的参考蒸散量过程对比；

图5为本发明相对于Hargreaves-Samani模型原始参数模拟方法的提升效果图。

具体实施方式

下面通过实施例，并结合附图，对本发明技术方案进行详细说明，但是本发明的保护范围不局限于所述实施例。

在本实施例中是以中国600个气象站1961-2019年逐日气象观测数据为基础，进行计算得到，在具体操作中，采用10折交叉验证方法将600站点随机分成互斥的10 组样本集合(每组60个站点)，每次选取其中1组共60个站点视作为缺资料地区，其余9组共540个站点视作为有资料地区，共抽样10次，即每个站点均有一次被视作为缺资料地区。

实施例

如图1，一种估算缺资料地区参考作物蒸散量的方法，包括以下步骤：

(1)基础资料收集：收集区域内有资料站点的最高和最低气温、湿度、日照时数(或辐射)、风速和大气压强气象要素资料以及缺资料站点的最高气温、最低气温和站点纬度数据资料。

有资料区指具备长系列的最高和最低气温、湿度、日照时数(或辐射)、风速和大气压强气象观测信息；所述无资料区指缺乏长系列的气象观测信息，但具备多年平均降水量、多年平均最高气温、多年平均最低气温和多年平均潜在蒸散量气候态信息，以及具备高程和站点纬度位置信息。

(2)参数优化：以Penman-Monteith模型计算的ET0作为参考值，以KGE系数作为目标函数，采用SCE-UA优化算法对有资料区站点的Hargreaves-Samani模型参数进行优化。

Penman-Monteith模型计算参考作物蒸散量ET0公式如下：

其中：Δ为饱和水汽压温度曲线斜率，由观测的最高气温和最低气温计算获得；R

目标函数KGE系数计算公式如下：

其中：Cov

Hargreaves-Samani模型计算公式如下：

其中：T

SCE-UA优化算法是一种全局优化算法，以PM模型计算的ET0为参考值，以KGE 系数最大为优化目标，采用SCE-UA优化算法求解得到每个有资料区站点的 Hargreaves-Samani模型参数a

(3)地理因子选择：以方便获取为基本原则，选取每个气象站的地理因子，包括高程、多年平均降水量、多年平均最高气温、多年平均最低气温、多年平均潜在蒸散量和干旱指数(多年平均潜在蒸散量/多年平均降水量)。

(4)基于机器学习的参数区域化模型：根据有资料地区的气象站点地理因子和优化后的Hargreaves-Samani模型参数a

机器学习模型采用随机森林，在本实施例的机器学习建模中，采用10折交叉验证方法将中国的600站点随机分成互斥的10组样本集合(每组60个站点)，每次训练中机器学习模型抽取其中9组共540个站点作为训练集(本实施例中将其视作为有资料地区)，剩余1组共60个站点作为测试集(本实施例中将其视作为缺资料地区)，共训练10次。相应的，共进行了10次参数模拟测试，600个气象站均参与到了缺资料情况下的参数区域化，就有600个缺资料站点的模拟结果，将600个气象站经过参数区域化模拟后的参数值作为更新后的缺资料地区Hargreaves-Samani模型参数值。图2为随机森林参数区域化模拟的参数值与SCE-UA优化参数值的散点关系图，从图中可以看出，随机森林模拟参数值与优化参数值的散点基本分布在图形对角线周边，拟合的KGE 系数均超过0.6，表明随机森林机器学习参数区域化模型很好的再现了缺资料地区的最优参数值。

(5)基于步骤(4)中参数区域化方法模拟的Hargreaves-Samani模型参数，采用最高气温、最低气温和站点纬度数据，根据步骤(2)公式进一步计算全国600个气象站的参考作物蒸散量，并与PM模型计算的ET0进行对比，计算本方法模拟的KGE精度系数。图3展示了某一站点Hargreaves-Samani模型区域化模拟的参数值计算的参考蒸散量和PM模型计算得到的参考蒸散量的对比。同时，基于Hargreaves-Samani模型的原始参数(即FAO的推荐取值)、最高气温、最低气温和站点纬度数据，根据步骤(2) 公式计算全国600个气象站的Hargreaves-Samani模型参考作物蒸散量，并与PM模型计算的ET0进行对比，图4展示了某一站点Hargreaves-Samani模型原始参数计算的参考蒸散量和PM模型计算的参考蒸散量的对比。计算原始Hargreaves-Samani模型模拟的 KGE精度系数，将本方法的KGE系数减去原始Hargreaves-Samani模型模拟的KGE精度系数，得到本方法的提升效果。图5为本发明相对于Hargreaves-Samani模型原始参数模拟方法的提升效果图，提升效果表示本方法的KGE系数减去原始Hargreaves-Samani 方法的KGE系数。

从图3和图4的对比可以看出，基于参数区域化模拟的参数值，Hargreaves-Samani模型计算得到的参考蒸散量与PM模型计算结果最为接近(图3)，相对于原始参数模拟的参考蒸散量提升效果明显(图4)。如图5所示，本方法相对于Hargreaves-Samani 模型原始方法的提升效果在0～20％之间，能较好的提高缺资料地区的参考作物蒸散量计算精度。提升效果表示本方法计算的参考作物蒸散量KGE系数值相对于 Hargreaves-Samani模型原始参数计算的参考作物蒸散量KGE系数值的增加百分比。

如上，尽管参照特定的优选实施例已经表示和表述了本发明，但其不得解释为对本发明自身的限制。在不脱离所附权利要求定义的本发明的精神和范围前提下，可对其在形式上和细节上作出各种变化。