一种双季稻田土壤锌供应年份的推算方法及系统

1.本发明涉及农业土壤施肥预测领域，特别是涉及一种双季稻田土壤锌供应年份的推算方法及系统。


背景技术：

2.水稻土壤中锌的含量对水稻锌含量具有明显影响，当土壤中全锌含量小于30mg/kg时，表示土壤为缺锌状态，在土壤缺锌状态下，水稻产量以及品质都会降低。因此，在稻田土壤中的锌消耗到极限之前及时补充锌肥，能够避免土壤因缺锌对水稻产量和品质造成的影响。因此，及时、准确地预测出土壤锌供应极限年限(即土壤中当前锌含量消耗到极限所需的时间)，显得很有必要。
3.目前，普遍是根据种植人员的个人经验预测出土壤锌供应极限年限，准确性较低，一旦预测错误，就容易发生在土壤锌消耗到极限后仍没有及时补充锌肥的情况。况且，种植人员没有有效快速简便方法和手段测定土壤中锌含量，更谈不上预测土壤锌消耗情况，只有当作物生长不佳，产量下降才会考虑是否是病害或者缺素。要确定缺素的门槛较高，测定土壤中锌的仪器较为精密，一般种植户根本不具备诊断条件，更谈不上预测了。因此，在有限的条件下，通过一次或少数几次土壤中锌养分的测定，及时、准确地推算出土壤锌供应极限年限，以避免发生在土壤锌消耗到极限后仍未及时补充锌肥导致水稻产量和品质降低的情况，是目前亟待解决的一个技术问题。


技术实现要素：

4.本发明的目的是提供一种双季稻田土壤锌供应年份的推算方法及系统，能够准确推算出土壤锌供应极限年限，实现在土壤缺锌前提前预判，提醒农户为土壤补充锌肥，避免因土壤严重缺锌导致的农作物产量和品质降低的问题。
5.为实现上述目的，本发明提供了如下方案：
6.一方面，本发明提出了一种双季稻田土壤锌供应年份的推算方法，所述方法包括：
7.获取待测定地区土壤的全锌数据；所述全锌数据包含所述待测定地区土壤的全锌含量值；
8.获取所述待测定地区的水稻种植数据，并计算所述待测定地区的双季水稻总产量值；所述双季水稻总产量值为同一年内种植的早稻和晚稻的产量总和；
9.根据所述待测定地区土壤的全锌含量值，判断所述待测定地区土壤当前所处的锌供应阶段；所述锌供应阶段包括锌供应极丰富阶段、锌供应丰富阶段和锌供应缺乏阶段；
10.根据所述待测定地区土壤的全锌含量值和双季水稻总产量值，计算得到从当前锌供应阶段到锌供应缺乏阶段所对应的土壤锌供应极限年限。
11.可选地，所述获取待测定地区土壤的全锌数据，具体包括：
12.采用hf-hno
3-hclo4消解-原子吸收分光光度法，在晚稻收获后对待测定地区土壤的全锌含量进行测定，得到所述待测定地区土壤的全锌含量值。
13.可选地，所述获取所述待测定地区的水稻种植数据，并计算所述待测定地区的双季水稻总产量值，具体包括：
14.采集所述待测定地区的水稻种植数据；所述水稻种植数据包括种植面积、水稻总重量、杂质含量、空瘪率和含水率；
15.根据所述待测定地区的水稻种植数据，分别计算所述待测定地区的早稻产量和晚稻产量；
16.对所述待测定地区的早稻产量和晚稻产量进行求和，计算得到所述待测定地区的双季水稻总产量。
17.可选地，利用下式分别计算所述待测定地区的早稻产量和晚稻产量：
18.yield＝(666.7/s)
×w×
(1－i)
×
(1－e)
×
[(1－m)/(1－mo)]；
[0019]
其中，yield表示产量，s表示种植面积，w表示水稻总重量，i表示杂质含量，e表示空瘪率，m表示含水率，mo为标准干重含水率；
[0020]
利用下式计算所述待测定地区的双季水稻总产量：
[0021]
yield_t＝yield_e+yield_l；
[0022]
其中，yield_t表示双季水稻总产量，yield_e表示早稻产量，yield_l表示晚稻产量。
[0023]
可选地，在所述根据所述待测定地区土壤的全锌含量值，确定所述待测定地区土壤当前所处的锌供应阶段的步骤之前，还包括以下步骤：
[0024]
采集多个地区不同年份下的历史水稻产量数据和历史全锌含量数据；
[0025]
根据所述历史全锌含量数据，将土壤全锌含量水平按区间划分为含锌极丰富水平区间、含锌丰富水平区间和锌缺乏水平区间，分别对应锌供应极丰富阶段、锌供应丰富阶段和锌供应缺乏阶段；
[0026]
根据所述历史水稻产量数据，分别计算在不同产量水平下，含锌极丰富水平区间、含锌丰富水平区间和锌缺乏水平区间对应的土壤全锌消耗速率；
[0027]
根据所述历史水稻产量数据、所述历史全锌含量数据及土壤全锌消耗速率数据，分别确定与锌供应极丰富阶段、锌供应丰富阶段和锌供应缺乏阶段对应的土壤锌供应极限年限计算公式。
[0028]
可选地，采用下式计算所述土壤全锌消耗速率：
[0029]
v＝dx/dt；
[0030]
其中，v表示土壤全锌消耗速率，dx表示从第i年份到第j年份之间的土壤全锌含量差值，dt表示从第i年份到第j年份之间的年份差值，i和j为正整数，且i＜j。
[0031]
可选地，所述根据所述待测定地区土壤的全锌含量值，判断所述待测定地区土壤当前所处的锌供应阶段，具体包括：
[0032]
当80≤qzn＜300mg/kg时，判定所述待测定地区土壤当前正处于锌供应极丰富阶段；
[0033]
当50≤qzn＜80mg/kg时，判定所述待测定地区土壤当前正处于锌供应丰富阶段；
[0034]
当30≤qzn＜50mg/kg时，判定所述待测定地区土壤当前正处于锌供应缺乏阶段；
[0035]
其中，qzn表示待测定地区土壤的全锌含量值。
[0036]
可选地，所述根据所述待测定地区土壤的全锌含量值和双季水稻总产量值，计算
得到从当前锌供应阶段到锌供应缺乏阶段所对应的土壤锌供应极限年限，具体包括：
[0037]
将所述待测定地区土壤的全锌含量值和双季水稻总产量值代入到土壤锌供应极限年限计算公式中，分别计算各锌供应阶段对应的土壤锌供应极限年限，包括：
[0038]
当所述当前锌供应阶段为锌供应极丰富阶段时，则分别计算出锌供应极丰富阶段对应的土壤锌供应极限年限、锌供应丰富阶段对应的土壤锌供应极限年限以及锌供应缺乏阶段对应的土壤锌供应极限年限；
[0039]
当所述当前锌供应阶段为锌供应丰富阶段时，则将锌供应极丰富阶段对应的土壤锌供应极限年限取0，并分别计算锌供应丰富阶段对应的土壤锌供应极限年限以及锌供应缺乏阶段对应的土壤锌供应极限年限；
[0040]
当所述当前锌供应阶段为锌供应缺乏阶段时，则将锌供应极丰富阶段对应的土壤锌供应极限年限以及锌供应丰富阶段对应的土壤锌供应极限年限均取0，并计算锌供应缺乏阶段对应的土壤锌供应极限年限；
[0041]
对所述锌供应极丰富阶段对应的土壤锌供应极限年限、所述锌供应丰富阶段对应的土壤锌供应极限年限以及所述锌供应缺乏阶段对应的土壤锌供应极限年限进行求和，得到所述待测定地区从当前锌供应阶段到锌供应缺乏阶段所对应的土壤锌供应极限年限。
[0042]
可选地，所述土壤锌供应极限年限计算公式，具体包括：
[0043]
采用jff公式计算锌供应极丰富阶段对应的土壤锌供应极限年限：
[0044][0045]
采用ff公式计算锌供应丰富阶段对应的土壤锌供应极限年限：
[0046][0047]
采用qf公式计算锌供应缺乏阶段对应的土壤锌供应极限年限：
[0048][0049]
其中，jff表示锌供应极丰富阶段对应的土壤锌供应极限年限，ff表示锌供应丰富阶段对应的土壤锌供应极限年限，qf表示锌供应缺乏阶段对应的土壤锌供应极限年限；qzn表示待测定地区土壤的全锌含量值，单位为mg/kg；yield_t表示双季水稻总产量，单位为kg/公顷。
[0050]
另一方面，本发明还提出了一种双季稻田土壤锌供应年份的推算系统，包括：
[0051]
全锌含量获取模块，用于获取待测定地区土壤的全锌数据；所述全锌数据包含所述待测定地区土壤的全锌含量值；
[0052]
双季水稻产量计算模块，用于获取所述待测定地区的水稻种植数据，并计算所述待测定地区的双季水稻总产量值；所述双季水稻总产量值为同一年内种植的早稻和晚稻的产量总和；
[0053]
土壤锌供应阶段确定模块，用于根据所述待测定地区土壤的全锌含量值，判断所述待测定地区土壤当前所处的锌供应阶段；所述锌供应阶段包括锌供应极丰富阶段、锌供应丰富阶段和锌供应缺乏阶段；
[0054]
土壤锌供应极限年限计算模块，用于根据所述待测定地区土壤的全锌含量值和双季水稻总产量值，计算得到从当前锌供应阶段到锌供应缺乏阶段所对应的土壤锌供应极限年限。
[0055]
根据本发明提供的具体实施例，本发明公开了以下技术效果：
[0056]
本发明提供了一种双季稻田土壤锌供应年份的推算方法及系统，方法包括：获取待测定地区土壤的全锌数据；全锌数据包含待测定地区土壤的全锌含量值；获取待测定地区的水稻种植数据，并计算待测定地区的双季水稻总产量值；双季水稻总产量值为同一年内种植的早稻和晚稻的产量总和；根据待测定地区土壤的全锌含量值，确定待测定地区土壤当前所处的锌供应阶段；锌供应阶段包括锌供应极丰富阶段、锌供应丰富阶段和锌供应缺乏阶段；根据待测定地区土壤的全锌含量值和双季水稻总产量值，计算得到从当前锌供应阶段到锌供应缺乏阶段所对应的土壤锌供应极限年限。本发明仅需要采集土壤全锌含量以及双季水稻总产量，根据土壤全锌含量即可判断出待测定地区土壤当前所处的锌供应阶段，即锌供应极丰富阶段、锌供应丰富阶段或者锌供应缺乏阶段，从而再根据土壤全锌含量以及双季水稻总产量，即可准确、可靠地推算出从当前锌供应阶段到锌供应缺乏阶段对应的土壤锌供应极限年限，从而能够根据推算结果适时补充锌肥，保证待测定地区水稻种植的产量和品质，避免发生在土壤锌消耗到极限后仍未及时补充锌肥导致水稻产量和品质降低的情况，为施肥时机提供参考。
附图说明
[0057]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。以下附图并未刻意按实际尺寸等比例缩放绘制，重点在于示出本发明的主旨。
[0058]
图1为本发明实施例1提供的双季稻田土壤锌供应年份的推算方法的流程图；
[0059]
图2为本发明实施例1提供的双季稻田土壤锌供应年份的推算方法的原理图；
[0060]
图3为本发明实施例2提供的双季稻田土壤锌供应年份的推算系统的结构框图。
具体实施方式
[0061]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完
整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0062]
如本发明和权利要求书中所示，除非上下文明确提示例外情形，“一”、“一个”、“一种”和/或“该”等词并非特指单数，也可包括复数。一般说来，术语“包括”与“包含”仅提示包括已明确标识的步骤和元素，而这些步骤和元素不构成一个排它性的罗列，方法或者设备也可能包含其他的步骤或元素。
[0063]
虽然本发明对根据本发明的实施例的系统中的某些模块做出了各种引用，然而，任何数量的不同模块可以被使用并运行在用户终端和/或服务器上。所述模块仅是说明性的，并且所述系统和方法的不同方面可以使用不同模块。
[0064]
本发明中使用了流程图用来说明根据本发明的实施例的系统所执行的操作。应当理解的是，前面或下面操作不一定按照顺序来精确地执行。相反，根据需要，可以按照倒序或同时处理各种步骤。同时，也可以将其他操作添加到这些过程中，或从这些过程移除某一步或数步操作。
[0065]
本发明的目的是提供一种双季稻田土壤锌供应年份的推算方法及系统，能够准确估算出土壤锌供应极限年限，从而能够在土壤缺锌前及时为土壤补充锌肥，避免因土壤缺锌导致的农作物产量和品质降低的问题。
[0066]
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
[0067]
实施例1
[0068]
如图1所示，本实施例提供了一种双季稻田土壤锌供应年份的推算方法，该方法具体包括：
[0069]
步骤s1、获取待测定地区土壤的全锌数据。所述全锌数据包含所述待测定地区土壤的全锌含量值。所述全锌含量表示的是土壤中包含的水溶态锌、代换态锌、难溶态锌和有机态锌等所有形态下的锌元素的含量总和。
[0070]
步骤s1具体包括以下步骤：
[0071]
采用hf-hno
3-hclo4消解-原子吸收分光光度法，在晚稻收获后对待测定地区土壤的全锌含量进行测定，得到所述待测定地区土壤的全锌含量值。
[0072]
原子吸收分光光度法简称原子吸收法，是利用被测元素基态原子蒸气对其共振辐射线的吸收特性进行元素定量分析的方法。该原子吸收分光光度法具有灵敏度高、精密度好、应用范围广、干扰少、快速简便等优点。可采用以下原子吸收分光光度法：首先在仪器推荐的浓度范围内，制备含锌元素的对照品溶液至少3份，浓度依次递增，并分别加入各品种项下制备供试品溶液的相应试剂，同时以相应试剂制备空白对照溶液。将仪器按规定启动后，依次测定空白对照溶液和各浓度对照品溶液的吸光度，记录读数。以每一浓度3次吸光度读数的平均值为纵坐标、相应浓度为横坐标，绘制标准曲线。按各品种项下的规定制备供试品溶液，使锌元素的估计浓度在标准曲线浓度范围内，测定吸光度，取3次读数的平均值，从标准曲线上查得相应的浓度，计算锌元素的含量，从而实现对待测定地区土壤的全锌含量的精准测定，得到准确、可靠的全锌含量值，从而为准确推算出土壤锌供应极限年限打下基础，保证土壤锌供应极限年限结果的准确性。本实施例采用的是hf-hno
3-hclo4消解-原子
吸收分光光度法，属于现有技术，具体过程在此不再赘述。
[0073]
步骤s2、获取所述待测定地区的水稻种植数据，并计算所述待测定地区的双季水稻总产量值；所述双季水稻总产量值为同一年内种植的早稻和晚稻的产量总和。
[0074]
步骤s2具体包括以下步骤：
[0075]
步骤s2.1、采集所述待测定地区的水稻种植数据；所述水稻种植数据包括种植面积、水稻总重量、杂质含量、空瘪率和含水率。
[0076]
本实施例中采集水稻种植数据，具体包括以下步骤：
[0077]
步骤s2.1.1、测量所述待测定地区的种植面积；
[0078]
步骤s2.1.2、对所述种植面积下同一年内收获的早稻和晚稻分别进行称重，得到早稻总重量和晚稻总重量；
[0079]
步骤s2.1.3、从早稻和晚稻中分别随机抽取1/10的水稻进行去杂、称重，测定杂质含量；
[0080]
步骤s2.1.4、从去杂后的早稻和晚稻中分别随机抽取1kg的水稻进行水分测定，对其烘干至含水量小于20％，并剔出空瘪粒，测定早稻和晚稻的空瘪率；
[0081]
步骤s2.1.5、采用谷物水分速测仪分别测量早稻和晚稻的含水率，重复10次取含水率平均值。
[0082]
步骤s2.2、根据所述待测定地区的水稻种植数据，分别计算所述待测定地区的早稻产量和晚稻产量。
[0083]
本实施例中，利用下式分别计算所述待测定地区的早稻产量和晚稻产量：
[0084]
yield＝(666.7/s)
×w×
(1－i)
×
(1－e)
×
[(1－m)/(1－mo)]；
[0085]
其中，yield表示产量；s表示种植面积，单位为m2；w表示水稻总重量，单位为kg；i表示杂质含量；e表示空瘪率；m表示含水率；mo为标准干重含水率，本实施例中，籼稻的标准干重含水率为13.5％，粳稻的标准干重含水率为14.5％，北方地区以种植粳稻为主，南方地区以种植籼稻为主。
[0086]
步骤s2.3、对所述待测定地区的早稻产量和晚稻产量进行求和，计算得到所述待测定地区的双季水稻总产量。
[0087]
本实施例中，利用下式计算所述待测定地区的双季水稻总产量：
[0088]
yield_t＝yield_e+yield_l；
[0089]
其中，yield_t表示双季水稻总产量，yield_e表示早稻产量，yield_l表示晚稻产量。
[0090]
在步骤s3之前，还包括以下步骤：
[0091]
(1)采集多个地区不同年份下的历史水稻产量数据和历史全锌含量数据。
[0092]
(2)根据所述历史全锌含量数据，将土壤全锌含量水平按区间划分为含锌极丰富水平区间、含锌丰富水平区间和锌缺乏水平区间，分别对应锌供应极丰富阶段、锌供应丰富阶段和锌供应缺乏阶段。
[0093]
(3)根据所述历史水稻产量数据，分别计算在不同产量水平下，含锌极丰富水平区间、含锌丰富水平区间和锌缺乏水平区间对应的土壤全锌消耗速率。
[0094]
本实施例中，采用下式计算所述土壤全锌消耗速率：
[0095]
v＝dx/dt；
[0096]
其中，v表示土壤全锌消耗速率，dx表示从第i年份到第j年份之间的土壤全锌含量差值，dt表示从第i年份到第j年份之间的年份差值，i和j为正整数，且i＜j。
[0097]
进一步的，根据土壤全锌消耗速率计算公式变换得到dt＝dx/v，其中，被除数就是土壤全锌消耗速率，在实际应用时，即可根据第i年份到第j年份之间的土壤全锌含量差值以及土壤全锌消耗速率，即可计算得到从第i年份到第j年份之间的年份差值，利用这一原理，从而在实际应用时可推算出锌供应极限年限，也即可确定锌元素消耗殆尽的具体年份。
[0098]
本实施例在划分锌供应极丰富阶段、锌供应丰富阶段和锌供应缺乏阶段时，根据大量的历史水稻产量数据和历史全锌含量数据，以时间(年份)为横坐标，以土壤全锌含量为纵坐标，绘制曲线图。根据曲线斜率的不同，划分出3个锌供应阶段区间，即锌供应极丰富阶段、锌供应丰富阶段和锌供应缺乏阶段，根据各个区间上限和下限两个土壤全锌含量的差值，以及所在区间的年份差值，计算得到土壤全锌消耗速率。
[0099]
(4)根据所述历史水稻产量数据、所述历史全锌含量数据及土壤全锌消耗速率数据，分别确定与锌供应极丰富阶段、锌供应丰富阶段和锌供应缺乏阶段对应的土壤锌供应极限年限计算公式。
[0100]
本发明中不同锌供应阶段对应的土壤锌供应极限年限计算公式，具体包括：
[0101]
采用jff公式计算锌供应极丰富阶段对应的土壤锌供应极限年限：
[0102][0103]
采用ff公式计算锌供应丰富阶段对应的土壤锌供应极限年限：
[0104][0105]
采用qf公式计算锌供应缺乏阶段对应的土壤锌供应极限年限：
[0106][0107]
其中，jff表示锌供应极丰富阶段对应的土壤锌供应极限年限，ff表示锌供应丰富阶段对应的土壤锌供应极限年限，qf表示锌供应缺乏阶段对应的土壤锌供应极限年限；qzn表示待测定地区土壤的全锌含量值，单位为mg/kg；yield_t表示双季水稻总产量，单位为kg/公顷。
[0108]
步骤s3、根据所述待测定地区土壤的全锌含量值，判断所述待测定地区土壤当前
所处的锌供应阶段；所述锌供应阶段包括锌供应极丰富阶段、锌供应丰富阶段和锌供应缺乏阶段。
[0109]
本发明判断待测定地区土壤的当前锌供应阶段，具体为：
[0110]
当80≤qzn＜300mg/kg时，判定所述待测定地区土壤当前正处于锌供应极丰富阶段。
[0111]
当50≤qzn＜80mg/kg时，判定所述待测定地区土壤当前正处于锌供应丰富阶段。
[0112]
当30≤qzn＜50mg/kg时，判定所述待测定地区土壤当前正处于锌供应缺乏阶段。其中，qzn表示待测定地区土壤的全锌含量值。
[0113]
需要说明的是，本实施例中的锌供应缺乏阶段表示的是当前全锌含量比一般水平稍低，农作物也会正常生长，但长势可能稍差(此时也可以适当补充锌肥)，一旦当土壤全锌含量处于锌供应缺乏阶段以下时，即qzn＜30mg/kg时，此时全锌含量不利于农作物生长，土壤中的锌元素基本消耗殆尽、极度匮乏，此时必须立即补充锌肥，否则水稻等农作物会明显降低产量和质量，甚至会生长缓慢乃至死亡。通常情况下，qzn＜30mg/kg或者qzn≥300mg/kg时，分别处于锌极度匮乏和富含锌的状态，都不利于农作物生长，对农作物的产量和质量都会造成影响。
[0114]
步骤s4、根据所述待测定地区土壤的全锌含量值和双季水稻总产量值，计算得到从当前锌供应阶段到锌供应缺乏阶段所对应的土壤锌供应极限年限。
[0115]
步骤s4具体包括以下步骤：
[0116]
步骤s4.1、将所述待测定地区土壤的全锌含量值和双季水稻总产量值代入到土壤锌供应极限年限计算公式中，分别计算各锌供应阶段对应的土壤锌供应极限年限，如图2所示，包括以下三种情况：
[0117]
1)当所述当前锌供应阶段为锌供应极丰富阶段时，则分别计算出锌供应极丰富阶段对应的土壤锌供应极限年限、锌供应丰富阶段对应的土壤锌供应极限年限以及锌供应缺乏阶段对应的土壤锌供应极限年限。
[0118]
由于当80≤qzn＜300mg/kg时，判定所述待测定地区土壤当前正处于锌供应极丰富阶段；当50≤qzn＜80mg/kg时，判定所述待测定地区土壤当前正处于锌供应丰富阶段；当30≤qzn＜50mg/kg时，判定所述待测定地区土壤当前正处于锌供应缺乏阶段。因此，当所述当前锌供应阶段为锌供应极丰富阶段时，需要分别计算锌供应极丰富阶段对应的土壤锌供应极限年限jff、锌供应丰富阶段对应的土壤锌供应极限年限ff以及锌供应缺乏阶段对应的土壤锌供应极限年限qf。
[0119]
2)当所述当前锌供应阶段为锌供应丰富阶段时，则将锌供应极丰富阶段对应的土壤锌供应极限年限取0，并分别计算锌供应丰富阶段对应的土壤锌供应极限年限以及锌供应缺乏阶段对应的土壤锌供应极限年限。
[0120]
如果qzn＜80mg/kg，则jff＝0，表示当前全锌含量处于锌供应极丰富阶段以下的阶段，所以锌供应极丰富阶段对应的土壤锌供应极限年限jff取0。此时，此时若处于锌供应丰富阶段，则需分别计算锌供应丰富阶段对应的土壤锌供应极限年限ff以及锌供应缺乏阶段对应的土壤锌供应极限年限qf。
[0121]
3)当所述当前锌供应阶段为锌供应缺乏阶段时，则将锌供应极丰富阶段对应的土壤锌供应极限年限以及锌供应丰富阶段对应的土壤锌供应极限年限均取0，并计算锌供应
缺乏阶段对应的土壤锌供应极限年限。
[0122]
如果qzn＜50mg/kg，则jff＝0，ff＝0，表示当前全锌含量处于锌供应丰富阶段以下的阶段，所以锌供应极丰富阶段对应的土壤锌供应极限年限jff以及锌供应丰富阶段对应的土壤锌供应极限年限ff均取0。此时若处于锌供应缺乏阶段，则仅需计算锌供应缺乏阶段对应的土壤锌供应极限年限qf。
[0123]
步骤s4.2、对所述锌供应极丰富阶段对应的土壤锌供应极限年限、所述锌供应丰富阶段对应的土壤锌供应极限年限以及所述锌供应缺乏阶段对应的土壤锌供应极限年限进行累加求和，得到所述待测定地区从当前锌供应阶段到锌供应缺乏阶段所对应的土壤锌供应极限年限。
[0124]
本实施例以双季水稻产量与土壤全锌含量为基础，且全锌含量范围在30-300mg/kg(一般稻田的锌含量范围)，计算土壤锌供应极限年限，jff为锌供应极丰富阶段时对应的供应极限年限，ff为锌供应丰富阶段对应的供应极限年限，qf为锌供应缺乏阶段时对应的供应极限年限。并且，由于随时间锌元素的消耗会经历从极丰富到丰富再到缺乏的变化过程，因此，最终得到的土壤锌供应极限年限实际上是上面各个锌供应阶段的供应极限年限的累加值gjy，表示为gjy＝jff+ff+qf，如图2所示。
[0125]
为了使本发明的推算土壤锌供应极限年限的过程更加详细，下面进行举例说明：
[0126]
例1：假设某地水稻产量水平为9251公斤/公顷，当前土壤全锌含量为159mg/kg，由于80≤159＜300mg/kg时，则判定全锌含量位于极丰富区间内，即处于锌供应极丰富阶段。
[0127]
首先，土壤全锌含量从极丰富区间降低到丰富区间，就会有个消耗时间，即jff，可以通过计算得出；其次，全锌含量进入丰富区间后，在丰富区间降低到缺乏区间以前，也会有个消耗时间，即ff，可以通过计算得出；当全锌含量降低到缺乏区间，而没有进入极度匮乏需要立即补充状态时，这个过程的时间为qf，所有会有三个阶段，三个时间。土壤锌供应极限总年限为这三个时间之和，计算过程如下：
[0128]
(1)jff＝(qzn1-80)/7.0＝(159-80)/7.0＝11.3
[0129]
(2)ff＝(qzn2-50)/1.9＝(80-50)/1.9＝15.8
[0130]
(3)qf＝(qzn3-30)/3.2＝(50-30)/3.2＝6.3
[0131]
(4)gjy＝jff+ff+qf＝11.3+15.8+6.3＝33.6
[0132]
根据本发明的方法，可以推算出土壤锌供应极限年限为33.6年，也就是说，在目前种植制度与常规管理下，大约33.6年后土壤将极度匮乏，需要立即补充锌肥，否则会严重降低农作物的产量和质量。
[0133]
例2：假设某地的水稻产量水平为11200公斤/公顷，当前土壤全锌含量为48mg/kg，由于30≤48＜50mg/kg，则判定全锌含量位于缺乏区间内，即处于锌供应缺乏阶段。计算过程如下：
[0134]
(1)因为48＜80mg/kg，因此，jff＝0
[0135]
(2)因为48＜50mg/kg，因此，ff＝0
[0136]
(3)qf＝(qzn3-30)/2.3＝(48-30)/2.3＝7.8
[0137]
(4)gjy＝jff+ff+qf＝0+0+7.8＝7.8
[0138]
根据本发明的方法，可以推算出土壤锌供应极限年限为7.8年，也就是说，在目前种植制度与常规管理下，大约7.8年后土壤将消耗殆尽、极度匮乏，需要立即补充锌肥，否则
会严重降低农作物的产量和质量。
[0139]
通过本发明可以较准确地预测出双季水稻出现锌匮乏所需要的时间，并评估出土壤全锌含量所处的阶段，能够为长期不施锌肥的农户提供施肥参考，为稻田施入锌肥量提供依据，能够保障农作物的产量以及质量。同时，由于有些复合肥料中随意盲目添加大量锌肥，而农户在不知自家土壤情况下施入添加大量锌肥的复合肥，从而造成的土壤锌过量和水稻毒害。为了避免土壤锌污染的发生，如果农户知道自家的土壤当前锌含量，可以判断出当前是否需要施入锌肥，施肥的时机是何时，因此，本发明还能够为农户购买和使用含锌复合肥，提供依据和指导。
[0140]
实施例2
[0141]
相对于实施例1中的一种双季稻田土壤锌供应年份的推算方法，本实施例提出了一种双季稻田土壤锌供应年份的推算系统，如图3所示，该系统具体包括：
[0142]
全锌含量获取模块m1，用于获取待测定地区土壤的全锌数据；所述全锌数据包含所述待测定地区土壤的全锌含量值；
[0143]
双季水稻产量计算模块m2，用于获取所述待测定地区的水稻种植数据，并计算所述待测定地区的双季水稻总产量值；所述双季水稻总产量值为同一年内种植的早稻和晚稻的产量总和；
[0144]
土壤锌供应阶段确定模块m3，用于根据所述待测定地区土壤的全锌含量值，判断所述待测定地区土壤当前所处的锌供应阶段；所述锌供应阶段包括锌供应极丰富阶段、锌供应丰富阶段和锌供应缺乏阶段；
[0145]
土壤锌供应极限年限计算模块m4，用于根据所述待测定地区土壤的全锌含量值和双季水稻总产量值，计算得到从当前锌供应阶段到锌供应缺乏阶段所对应的土壤锌供应极限年限。
[0146]
除非另有定义，这里使用的所有术语(包括技术和科学术语)具有与本发明所属领域的普通技术人员共同理解的相同含义。还应当理解，诸如在通常字典里定义的那些术语应当被解释为具有与它们在相关技术的上下文中的含义相一致的含义，而不应用理想化或极度形式化的意义来解释，除非这里明确地这样定义。
[0147]
上面是对本发明的说明，而不应被认为是对其的限制。尽管描述了本发明的若干示例性实施例，但本领域技术人员将容易地理解，在不背离本发明的新颖性和优点的前提下可以对示例性实施例进行许多修改。因此，所有这些修改都意图包含在权利要求书所限定的本发明范围内。应当理解，上面是对本发明的说明，而不应被认为是限于所公开的特定实施例，并且对所公开的实施例以及其他实施例的修改意图包含在所附权利要求书的范围内。本发明由权利要求书及其等效物限定。