基于作物模型、遥感数据和气候预测信息的作物产量早期预报技术的制作方法

本发明涉及农业信息技术领域，具体而言涉及一种基于作物模型、遥感数据和气候预测信息的作物产量早期预报技术。

背景技术：

作物生长监测和产量早期预报信息对及时优化农业管理、减灾防灾、预测农业产品贸易行情及价格至关重要。因此，从上世纪80年代以来，该研究领域就受到广泛关注，试图通过遥感、作物模型和各种气候预测信息发展作物生长监测和粮食安全早期预警技术及系统。其中，作物生长监测和产量早期预报技术是发展粮食安全预警系统的核心技术。近年来，国内外学者结合作物模型、气候预测、遥感数据融合等技术对作物产量预报的可行性进行了探索，并在站点尺度取得了一定的成果。相比之下，目前针对大面积作物生长监测和产量早期预报技术的探索仍处于初级阶段，主要依靠遥感数据产品(如遥感植被指数、遥感叶面积指数)和简单的统计学模型进行预报。预报技术的机理性不强，难以广泛推广到实际应用中。同时，准确度也有待提升，对实际生产的指导意义有所欠缺。另外，现有方法没有系统地对预报时效进行阐述，这使得“作物产量可以提前多久进行预报”这一问题始终没有明确的答案。因此，目前大面积作物产量早期预报仍面临准确度、机理性、时效性等方面的限制。因此，发展一套机理性强，易于推广，且可以尽早地、足够准确地预报区域作物产量的技术及系统，对于客观、及时地发布预报信息至关重要，也是该领域的研究热点以及亟待解决的关键问题。

技术实现要素：

为了至少部分解决上述现有技术中存在的问题，本发明将作物机理模型、遥感观测信息、数据融合新方法(即空间格局同化方法)、气候早期预报以及集合预报、滚动预报技术等结合在一起，由此提出了一种作物产量早期预报技术。

根据本发明一方面，提供了一种基于作物模型、遥感数据和气候预测信息的作物产量早期预报技术，包括如下步骤：

s1:获取研究区域内的预报日期d之前的观测气象数据以及遥感数据；

s2：采用所述观测气象数据来驱动作物机理模型，从作物播种日期d0开始，逐日模拟作物生长发育，直至预报日期d；在该步骤过程中，利用从日期d0到d之间的遥感数据，对作物机理模型进行同化，所述遥感数据包括叶面积指数lai；

s3：获取预报日期d之后的气候预测数据，并利用该气候预测数据驱动作物模型继续运行，模拟至作物成熟日期d1，然后输出模拟产量y；

其中，所述气候预测数据包括多套预测数据，利用所述多套预测数据分别进行模拟，获得多个模拟产量，采用所述多个模拟产量的平均值作为预报产量；以及

s4：调整预报日期d，以在日期d0到d1之间获得多个不同的预报日期，然后针对所述多个不同的预报日期，重复步骤s1-s3，由此获得针对不同预报日期的不同预报产量。

根据本发明的实施方案，所述作物机理模型可以选自mcwla系列模型、ceres系列模型、apsim模型、wofost模型。

根据本发明的实施方案，所述作物可以选自玉米、水稻、小麦和大豆等。

根据本发明的实施方案，所述气象数据可以包括日太阳辐射、日最高温度、日最低气温、日降雨量、日相对湿度和日平均风速。

根据本发明的实施方案，各个预报日期之间的间隔可以为7-10天，例如8，9天。

根据本发明的实施方案，所述基于作物模型、遥感数据和气候预测信息的作物产量早期预报技术还包括步骤s5：将各预报产量与实际统计产量相比较，得到不同的预报误差；然后根据预报误差的大小，确定研究区域内所述作物的对应预报的时长。

根据本发明的实施方案，可以采用平均绝对百分百误差(mape)来表示误差。

本发明将作物机理模型、遥感观测信息、数据融合(即空间格局同化方法)、气候早期预报以及集合预报、滚动预报技术等结合在一起，发展了一套农作物产量早期预报技术，并且能够实现有益的技术效果：

(1)可以及时地、准确地实现大面积作物产量预报。该技术准确度高、机理性强、易于推广；

(2)利用滚动预报技术对“作物产量可以提前多久进行预报”这一问题进行了回答，阐明了有效预报时间；

(3)本发明显著提升了对产量预报潜力的认识，并将极大地推动产量预报技术的实用化。

附图说明

附图中相同的附图标记标示了相同或类似的部件或部分。本发明的目标及特征考虑到如下结合附图的描述将更加明显，附图中：

图1是根据本发明的一个实施方案的研究区及基于遥感叶面积指数提取的冬小麦种植区。

图2是根据本发明的一个实施方案的基于作物模型、遥感数据和气候预测信息的作物产量早期预报技术的流程示意图。

图3是根据本发明一个实施方案的使用不同气象预测数据时的预报精度结果图。

图4是根据本发明一个实施方案的对比三种预报方案实现不同预报精度的概率随预报时间的变化结果图。

具体实施方式

为清楚地说明本发明中的方案，下面给出优选的实施例并结合附图详细说明。以下的说明本质上仅仅是示例性的而并不是为了限制本公开的应用或用途

应该理解的是，本发明所引用的作物模型、遥感模块以及同化方法等本身是已知的，例如模型的各个子模块、各种参数、运行原理机制等等，因此本发明重点阐述如何作物机理模型、遥感观测信息、气象数据、气候早期预报等结合在一起，以设计农作物产量早期预报技术。

本发明的基于作物模型、遥感数据和气候预测信息的作物产量早期预报技术基于作物机理模型模拟，可以针对年y的作物生长季产量预报，假设当前日期为d，即需要在日期d预测年份y的作物产量，从播种期日期d0开始，采用观测气象数据驱动作物机理模型逐日(也即，模拟步长为日尺度)模拟作物生长发育直至模拟至日期d，d也被称为预报日期。

所述气象数据例如可以包括日太阳辐射、日最高温度、日最低气温、日降雨量、日相对湿度和日平均风速等。观测数据例如可以来观测站和气象局的发布数据。另外，作物机理模型无具体模型限制，例如可以为mcwla系列模型、ceres系列模型、apsim模型、wofost模型等。本发明可应用于多种作物，例如玉米、水稻、小麦和大豆等。

在模拟过程中，如遇到有遥感数据(例如叶面积指数lai等)的日期，则可以进行数据同化操作(本地化)，例如利用空间格局同化方法，对作物机理模型模拟的作物生长特征参数(如lai)和遥感数据进行同化操作。本操作在模拟至预报日期d前实施，一般存在多个日期可以进行同化操作。应该理解的是，上述同化为本领域中所熟知，本文中不再赘述。本发明方法对于多种遥感数据具有普适性，同化运算的空间分辨率及时间频率可以与遥感数据的时空分辨率保持一致。例如lai可以采用copernicuslai数据。

当模型模拟至日期d之后，进入产量预报阶段。使用预测的气象数据序列驱动作物模型继续运行，直至模型模拟作物成熟日期d1，输出模拟产量y。

产量预报阶段所使用的预测的气象数据可以不是单一数据，而是采用多数据序列集合模拟技术进行操作。针对每一个预测的气象数据序列，均可得到一个对应的模拟产量。当使用n个预测气象数据序列(也即，多个气象研究中心发布的、针对同一个时段某个地区的多套气候预测数据)时，可得到y1，y2……yn共n组模拟产量。这些模拟产量的均值，即(y1，y2……yn)/n作为预报产量。所述气候预测数据可以包含多种来源，例如中国气象局(cma)气候预报产品，欧洲中期天气预报中心(ecmwf)气候预报产品、英国气象局(ukmo)气候预报产品等。亦可以利用历史观测数据来代表气候预测数据。

本发明中，所述产量可以匹配实际统计产量的空间尺度，例如作物模型运算所得产量为格网尺度，可以聚合至统计产量的尺度，例如县级。可以将各预报产量与实际统计产量相比较，得到不同的预报误差，可以采用平均绝对百分百误差(mape)、均方根误差(rmse)、或者其他合适的参数来表示误差；然后根据预报误差的大小，确定研究区域内所述作物的对应预报的时长。

本发明中，所述预报日期d为可变变量，不同的d值会影响预报准确度，d值越接近作物成熟期，预报准确度越高。通过滚动变化d值，例如d+1,d+2……d+m…d-1…可获得不同日期下的产量预报结果。对不同预报日期下产量预报准确度进行评价可以获得预报准确度随d值的滚动变化情况，并可以据此回答“作物产量可以提前多久进行预报”这一问题。

下面结合具体的实施例对本发明的技术进行进一步的详细说明。

实施例

下面以中国华北平原中部(图1)为研究区，研究对象为小麦，示例性说明本发明技术方法的具体应用。实施例用于说明本发明，但不用来限制本发明的范围。

更具体地，在本例中，产量预报技术被应用于2008-2015年进行产量预报回测以检验预报技术的性能。本示例采用mcwla-wheat模型作为示例作物机理模型，使用的遥感lai数据为copernicuslai(空间分辨率1km×1km，时间频率为每10天一期)；使用了多种气象数据作为预测数据，包括中国气象局(cma)气象预测产品，欧洲中期天气预报中心(ecmwf)气象预测产品、英国气象局(ukmo)气象预测产品，同时，历史气象观测数据也被用作可能的预测数据。

针对其中的一个年份y，构建作物机理模型mcwla-wheat，应该理解的是，作物机理模型mcwla-wheat为本领域技术人员所熟知，该模型的构建/本地化也是已知的，例如可以利用研究区内以及研究时段的土壤、气候、农田管理等数据来输入模型进行操作，在此就不具体赘述；然后从播种期d0起开始模拟，以年y实测气象数据(国家气象局气象站观测共享数据，包括日太阳辐射、日最高温度、日最低气温、日降雨量、日相对湿度和日平均风速)驱动所构建的mcwla-wheat。同时获取copernicuslai数据，当遇到遥感数据存在的日期，对模型和遥感数据进行空间格局同化，实时修正模型模拟过程；模型与遥感数据的同化没有特别的限制，并且该技术为本领域所熟知。

本示例进行了多个预报日期d(d1,d2,d3……dn)的滚动预报工作,针对每个预报日期d，作物模型耦合遥感数据在日期d之前进行模拟以及同化工作。在d之后将使用气候预测数据驱动模型运行至成熟期d1并输出模拟产量。图2是根据本发明提出的预报技术的一个实施方案的流程示意图。

如图2所示，本发明所提出的作物产量早期预报技术耦合了作物模型、作物模型与遥感数据同化新技术、以及气候预报信息。基于本技术，当前示例给出了一个针对华北地区冬小麦的先进的作物产量早期预报系统。随着遥感lai数据和气候预报信息的更新，开展滚动预报，在不同的预报日期，输出作物产量预报结果。本示例中，不同预报日期d(d1,d2,d3……dn)之间的间隔为8天。

本示例使用了多种气象预测数据进行产量预报效果展示，其中，中国气象局(cma)气象预测产品，欧洲中期天气预报中心(ecmwf)气象预测产品、英国气象局(ukmo)气象预测产品均提供了多套可能的气象预测。当使用历史气象观测作为气候预测数据时，使用预报年y之前30年(即y-1至y-30)的历史数据作为30套可能的气象预测。在使用上述数据时，每一套气象预测均被用于驱动预报日期d之后的模型运行，使用多套气象预测的产量预报结果的均值作为最终预报产量。

然后，针对每个预报日期输出的作物产量预报结果，产量将聚合至县级尺度，并与当年的实际统计产量记录进行对比。

预报准确度采用平均绝对百分比误差(mape)进行评价。

其中，oi,y为y年i县的产量统计结果，pi,dn,y为在y年i县预报日期为dn时的产量预报结果，n为产量模拟预报的组数，也即使用的气象数据的序列数。

本示例所得的2008-2015年产量预报结果以及的准确度评价结果如表1、图3和图4所示。

下表1为本示例研究区各年总的平均预报产量以及统计产量的结果。

表1表明，随着预报期缩短，使用历史气象数据、cma气象预测、ukmo气象预测的预报产量逐渐接近统计产量，其中历史数据表现最优，在提前44天时绝对误差已经在5％以下，cma表现次之，ukmo和ecmwf数据的产量预报结果表现再次之。

图3和图4为研究区中使用不同气象预测数据时的预报精度结果以及各种方案的比较图，其中：

图3的(a)中，基于过去30年历史气象数据的产量预报误差随随预报日期的变化而减小，提前60天预报的平均误差小于12％；

图3的(b)中，基于中国气象局(cma)气象预测的产量预报误差随预报日期的变化呈现下降趋势，提前60天预报的平均误差小于12％；

图3的(c)中，基于欧洲中期天气预报中心(ecmwf)气象预测的产量预报误差随预报日期的变化呈现下降趋势，提前44天预报的平均误差小于14％(ecmwf仅提供最长44天的气象预测)；

图3的(d)中，基于英国气象局(ukmo)气象预测的产量预报误差随预报日期的变化呈现下降趋势，提前60天预报的产量误差小于15％；

图3的(e)作为对照，基于实测气象数据的产量模拟误差随预报日期的变化逐渐减小，提前60天预报的产量误差接近10％；

图3的(f)中，基于不同气候预测数据的产量预报与基于实测气候数据模拟的差异随预报时间的变化表明，基于历史气候数据的集合预报误差较小，其次是基于中国气象局(cma)的气候预报信息。

基于图3的结果，图4展示了实现不同预报准确度的概率。此处选取了基于中国气象局(cma)气象预测的产量预报结果，和基于历史气候数据的产量预报结果，并使用基于实测气候资料的产量模拟结果作为对照。结果表明利用历史数据开展产量集合预报，提前52天小麦产量预报误差小于12.5％的概率大于80％，提前44天小麦产量预报误差小于12.5％的概率可以达到90-100％。在概率≥80％的前提下，实现预报误差小于10.0％，12.5％和15％的对应预报时长分别为提前36、52和60天。基于中国气象局(cma)气象预测的产量预报,在概率≥80％的前提下，实现预报误差小于10.0％，12.5％和15％的对应预报时长分别为提前28、36和36天。

本示例基于上述操作和结果，展示了本发明提出的产量预报技术的实际操作过程以及预报结果。上述结果表明，基于本发明所提出的产量预报技术，结合历史观测数据，可以实现较高准确度的华北地区冬小麦产量预报。针对本示例中华北地区冬小麦产量预报，可以实现提前5周至两个月对产量进行较为准确的预报。本技术可以广泛的应用于各类作物、区域的推广。需要指出的是，针对不同地区、不同作物，可能需要采用不同的遥感或气象预测数据以实现最优的预报效果。本发明可以及时地、准确地实现大面积作物产量预报。同时，基于本发明所提出的滚动预报技术，可以对“作物产量可以提前多久进行预报”这一问题进行回答，阐明有效预报时间。

本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的装置及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。