基于作物生长机理模型的关键因子产量贡献计算方法
【技术领域】
[0001] 本发明涉及一种计算决定作物产量的关键因子对产量贡献的方法，具体涉及利用 APS頂作物生长机理模型分离和评价作物品种变化、管理措施改变以及气候变化(30年及以 上)对作物产量影响的计算方法，属于农业可持续发展适应性评估技术领域。
【背景技术】
[0002] 在制约农业生产的自然资源中，气候条件是最为重要的组成部分。气候变化对农 业生产的影响已经引起各国政府和科学家们的广泛关注和高度重视。伴随着全球地表温度 的持续升高，一些地区地表辐射降低和降雨量减少等气候因子的变化均对作物生产和产量 形成产生了重要的影响。因此，准确量化作物产量对过去气候变化的响应及其敏感性是理 解和预测未来气候变化对农业生产影响的前提和基础。
[0003] 过去几十年里，作物育种技术迅速发展，人们为了获得更高的产量，平均每3-5会 变换新的作物品种；同时，人们为了获得最大的经济效益，往往采取更好的田间管理措施， 主要表现为加大化肥投入和农业灌溉。因此，品种更替和管理措施改变对作物生产和最终 作物产量起了关键性的作用。
[0004]总之，过去几十年作物生产受到了气候变化、品种变换和管理措施改变等诸多因 子的共同影响。目前研究关键因子产量贡献的主要方法有统计回归和作物模型模拟方法。 统计模型分析简便易行，分析结果的可信度也比较高，因而得到广泛应用。但影响作物生产 的各个因子相互之间往往不是独立的，因此，统计模型不能揭示众多影响因子之间的相互 关系，因此难以利用其分析结果进一步提出针对性较强的适应性对策;然而基于过程的作 物生长机理模型可以克服这一缺点，作物模型可以模拟不同情景下的作物生长，因此可以 摒弃因子之间的相互作用，准确分离和量化各个因子对作物产量的影响。
[0005] APSIM(Agricultural Production Systems Simulator，农业生产系统模拟模型） 由隶属澳大利亚联邦科工组织(CSIR0)和昆士兰州政府的农业生产系统研究组(APSRU)在 过去20多年内开发的，是能够模拟农业系统各主要组分的机理模型。该模型可用于模拟农 业系统中作物生长过程及土壤水氮动态，特别适用于评价农作系统生产潜力及耕作措施的 产量效益受气候波动和环境变化的影响。APSIM模型可以让用户容易地通过选择一系列的 作物、土壤以及其它子模块来配置自己的作物模型。模块之间的逻辑关系可以非常简单地 通过模块的"插拔"功能来规定。由于它的灵活性、可操作性，APSIM模型被认为更应该是一 个模型系统的灵活软件环境，而不是针对某特定作物系统的单个模型。
【发明内容】
[0006] 本发明要解决的技术问题是提供一种基于作物生长机理模型的关键因子产量贡 献计算方法，能够直观的分离品种、管理措施以及气候变化等关键因子对作物产量的贡献。
[0007] 为解决上述技术问题，本发明所采取的技术方案是：
[0008] -种基于作物生长机理模型的关键因子产量贡献计算方法，该方法包括如下步 骤：
[0009] (1)选择研究站点不同时期的多个作物品种，根据实际田间观测试验的作物物候 和产量数据确定各自的品种参数范围；
[0010] (2)基于步骤（1)中确定的品种参数范围，利用研究站点的气象数据和作物管理措 施数据驱动APS頂作物生长机理模型，采用试错法确定具体的作物品种参数值；
[0011] (3)根据步骤(2)中确定的作物品种参数值，对长时间序列的不同品种和不同管理 措施下的作物生长机理模型模拟产量值进行对比分析，分别计算品种变化和管理措施改变 对作物产量的贡献;并对作物品种在长时间序列上的作物产量模拟值进行线性拟合，计算 气候变化对作物产量的贡献。
[0012]进一步的，所述步骤(3)中：
[0013] (3-1)品种变化对作物产量的贡献Cc的计算公式为：
[0014]
[0015]其中，Y(C2) i是品柙Π 的作物模拟产量，Y(Cl)i是品种I的作物模拟产量，i是年份 (1,2,......,η)；
[0016] (3-2)管理措施改变对作物产量的贡献Cm的计算公式为：
[0017]
[0018]其中，Y(M2h是管理措施Π 条件下的作物模拟产量，Y(Mlh是管理措施I条件下的 作物模拟产量；
[0019] (3-3)气候变化对作物产量的贡献的计算方法为：
[0020] (3-3-1)首先对连续多年一致的品种和管理措施下作物模拟产量Y(C)i进行线性 拟合：
[0021] Y(C)i = aXi+b
[0022] 其中，Xi模拟年份(i = 1，2，……，n)，a是线性拟合斜率，b是线性拟合截距；
[0023] (3-3-2)计算长时间序列气候变化对作物产量的贡献Cel:
[0024]
[0025] 其中，η是模拟年份值。
[0026] 采用上述技术方案所产生的有益效果在于：
[0027] 本发明公开了一种基于作物生长机理模型的关键因子产量贡献计算方法，通过作 物田间实验观测资料调整APS頂作物生长机理模型品种参数，利用作物生长机理模型模拟 不同情景下作物生长过程，能够直观的分离出品种、管理措施以及气候变化等关键因子对 作物产量的贡献，对于指导品种选育和改善田间管理等农业生产措施来应对和适应气候变 化具有重要意义。
【附图说明】
[0028]图1是本发明实施的流程图；
[0029]图2是本发明夏玉米产量的实测值与模拟值验证；
[0030]图3是本发明关键因子产量贡献。
【具体实施方式】
[0031]下面结合附图和【具体实施方式】对本发明作进一步详细的说明。
[0032]如图1所示，本发明公开了一种基于APS頂作物生长机理模型的关键因子产量贡献 计算方法。
[0033]研究对象:过去30年（1981-2010年)华北地区气候变化、作物品种变换和管理措施 改变(主要考虑施肥量的变化)对夏玉米产量的贡献。
[0034] 1、确定研究站点不同时期典型作物品种和各自的品种参数范围
[0035] 根据国家农业气象站过去30年的试验观测数据，选择1981-1985和2006-2010两个 时期作为研究对象，确定4个代表性站点两个时期种植的典型作物品种，如表1所示。根据不 同品种的实际