一种大田玉米蒸发蒸腾量计算方法及系统与流程

本发明属于土壤水分参数及降水量检测分析技术领域，具体涉及一种大田玉米蒸发蒸腾量计算方法及系统。

背景技术：

目前，在研发基于玉米日耗水量的自动灌溉过程中，土壤水分状况和大田积累日蒸发蒸腾量是衡量是否启动灌溉的重要参数，测量大田玉米群体蒸发蒸腾量对研究大田作物需水量过程及特征具有重大意义。目前估算作物蒸发蒸腾量的方法主要有基于土壤贮水量变化(δw)的水量平衡法、彭曼公式法和涡度相关法。所述彭曼公式法通过将气象因子代入，计算出参考作物蒸发蒸腾量(et0)，然后乘以玉米不同时期的作物系数(kc)，估算出玉米蒸发蒸腾量；彭曼公式法所取的气象因子较多，获取较难，且作物系数随着品种、灌溉技术、种植区域的变化而变化，较难准确估算农田尺度玉米耗水量，该方法主要用监测农田小气候条件下的作物蒸发蒸腾量。

技术实现要素：

本发明解决的技术问题是提供了一种操作过程简单且测量结果精确的大田玉米蒸发蒸腾量计算方法及系统，该方法基于智墒产品0-30cm水分数据及降水数据实时计算黄淮南部大田玉米的蒸发蒸腾量。

本发明为解决上述技术问题采用如下技术方案，一种大田玉米蒸发蒸腾量计算方法，其特征在于具体过程为：

步骤s1：根据晴天分段0-10cm、10-20cm、20-30cm土壤体积含水率与大田玉米蒸发蒸腾量函数计算出各分段土层对大田玉米蒸发蒸腾量权重；

分段0-10cm土壤体积含水率与大田玉米蒸发蒸腾量函数关系式为：

etc1＝exp(0.6330+0.0123θ1+0.0011θ1²)r1²＝0.528^**

分段10-20cm土壤体积含水率与大田玉米蒸发蒸腾量函数关系式为：

etc2＝exp(4.5292－0.3271θ2+0.0077θ2²)r2²＝0.413^**

分段20-30cm土壤体积含水率与大田玉米蒸发蒸腾量函数关系式为：

etc3＝exp(4.6705－0.3151θ3+0.0069θ3²)r3²＝0.323^**

其中，etc1、etc2、etc3分别为分段0-10cm、10-20cm、20-30cm大田玉米的蒸发蒸腾量，θ1、θ2、θ3分别为分段0-10cm、10-20cm、20-30cm土壤体积含水率；

θ1、θ2、θ3分别对etc1、etc2、etc3贡献所占权重的算法为：

分段0-10cm土壤体积含水率θ1对etc1的权重为：

wt1＝r1²/(r1²+r2²+r3²)

分段10-20cm土壤体积含水率θ2对etc2的权重为：

wt2＝r2²/(r1²+r2²+r3²)

分段20-30cm土壤体积含水率θ3对etc3的权重为：

wt3＝r3²/(r1²+r2²+r3²)

其中，wt1、wt2、wt3分别为分段0-10cm、10-20cm、20-30cm土壤体积含水率对etc1、etc2、etc3的权重；r¹、r²、r³分别为决定系数，则大田玉米晴天的蒸发蒸腾量函数关系式为：.

etc＝etc1×wt1+etc2×wt2+etc3×wt3

根据分段0-10cm、10-20cm、20-30cm土壤体积含水率对大田玉米蒸发蒸腾量的权重计算出基于根层0-30cm土壤体积含水率的大田玉米蒸发蒸腾量；

步骤s2：根据大田玉米蒸发蒸腾量与雨天降水量函数计算出雨天大田玉米蒸发蒸腾量，大田玉米蒸发蒸腾量与雨天降水量函数式为：

etc＝5.720×pre^-0.606r²＝0.878^**

其中，pre为天圻所测量当日00:00至24:00实际测量的降水量，单位为mm；

步骤s3：将基于0-30cm土壤水分数据计算的晴天大田玉米蒸发蒸腾量etc值与基于降雨的雨天大田玉米蒸发蒸腾量etc值进行日期合并，计算出全生育期玉米的蒸发蒸腾量。

优选的，步骤s1中根据大田玉米根部0-10cm、10-20cm、20-30cm土壤体积含水率分别拟合出其与当天etc值的函数关系，根据spss软件levenberg–marquardt算法拟合出最优的拟合方程为指数方程与二次方程的结合方程，根据分段土壤体积含水率与etc的决定系数计算出其对大田玉米蒸发蒸腾量etc的权重。

优选的，步骤s1中晴天分段0-10cm、10-20cm、20-30cm土壤体积含水率通过智墒土壤水分传感器进行测量，该智墒土壤水分传感器安装于玉米行同行中间，测量半径15cm范围内外进行隔离处理，函数计算的大田玉米田间etc值长度均为1d。

本发明所述的大田玉米蒸发蒸腾量计算系统，包括测量模块和计算模块，其中测量模块用于获取大田玉米0-30cm各分段土壤体积含水率和雨天降水量；计算模块用于根据大田玉米各分段土壤体积含水率和etc的函数、土壤体积含水率对etc的权重及雨天降水量与etc的函数关系计算出大田玉米全生育期各种天气情况的蒸发蒸腾量；测量模块通过智墒水分传感器、天圻气象站分别获取0-30cm各分段体积含水率数据、雨天有效降水量数据，计算模块通过权重累加等方法对测量的数据进行运算分析，最终实现了对大田玉米在各种天气状况下的蒸发蒸腾量的有效计算，能够大规模应用在实际生产中。

本发明与现有技术相比具有以下有益效果：本发明提供的一种黄淮南部基于智墒产品0-30cm水分数据及降水数据的大田玉米蒸发蒸腾量计算方法，通过对大田玉米0-30cm分段土壤体积含水率和雨天降水量进行实时无损检测，实现了对大田玉米蒸发蒸腾量的有效估计，测量结果精确，并且能够在实际生产中大规模应用，并弥补智墒产品无法估算雨天etc的系统缺陷。

附图说明

图1是本发明大田玉米蒸发蒸腾量计算方法的流程图；

图2是本发明大田玉米蒸发蒸腾量计算方法及系统流程图。

具体实施方式

下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。

实施例

图1为根据本发明的大田玉米蒸发蒸腾量计算方法的流程图，如图1所示，一种大田玉米蒸发蒸腾量计算方法，包括：s1，根据晴天或多云天气0-30cm土壤各分段体积含水率，得到0-10cm、10-20cm和20-30cm含水量与etc函数，通过函数计算出基于0-30cm的平均日etc值。本研究选取0-30cm土层土壤水分是基于相关试验的前期拟合得出的(表1)；s2，基于天圻所测当日降水量，拟合出降水量与蒸渗仪etc的函数关系，用于计算雨天大田玉米的etc值；s3，根据所求的晴天与雨天的etc，将其合并累计，计算大田玉米全生育期的蒸发蒸腾量，弥补智墒产品在雨天没有etc值的缺陷。

表1称重式蒸渗仪etc值与对应的智墒水分仪分段土层水分的相关系数

注：**表示p<0.01水平显著相关。

上述实施例提供的一种大田玉米蒸发蒸腾量计算方法，通过对大田0-30cm各分段体积含水率和雨天降水量进行实时无损检测，实现了对黄淮南部大田玉米蒸发蒸腾量的有效估计，测量结果精确，并且能够在实际生产中大规模应用。

基于上述实施例的内容，大田玉米根部土壤体积含水量根据智墒水分传感器的测量值计算得到。所述s1具体包括：根据大田玉米根部0-10cm、10-20cm、20-30cm土壤体积含水率分别拟合出其与当天etc值的函数关系。根据spss软件levenberg–marquardt算法拟合出最优的拟合方程为指数方程与二次方程的结合方程。根据分段体积含水率与etc的决定系数r2计算出其对大田玉米etc的权重。所述s2具体包括：根据雨天的降水量与etc函数关系，计算得到所述大田玉米雨天的etc。所述s3具体包括：将晴天和雨天的etc合并得出大田玉米全生育期的etc。

具体地，将0-30cm分段的长度均为10cm，水分监测深度依次为10cm、20cm、30cm，土壤体积含水率依次为θ1、θ2、θ3，则大田玉米蒸发蒸腾量etc与各分段土壤含水率函数为：

则分段0-10cm土壤体积含水率与etc的关系函数式为：

etc1＝exp(0.6330+0.0123θ1+0.0011θ1²)r1²＝0.528^**

则分段10-20cm土壤体积含水率与etc的关系函数式为：

etc2＝exp(4.5292－0.3271θ2+0.0077θ2²)r2²＝0.413^**

则分段20-30cm土壤体积含水率与etc的关系函数式为：

etc3＝exp(4.6705－0.3151θ3+0.0069θ3²)r3²＝0.323^**

其中，etc为大田玉米的蒸发蒸腾量，θ1、θ2、θ3分别分段0-10cm、10-20cm、20-30cm土壤体积含水率；θ1、θ2、θ3对etc贡献所占权重的算法为：

则分段0-10cm土壤体积含水率对etc的权重为：

wt1＝r1²/(r1²+r2²+r3²)

则分段10-20cm土壤体积含水率对etc的权重为：

wt2＝r2²/(r1²+r2²+r3²)

则分段20-30cm土壤体积含水率对etc的权重为：

wt3＝r3²/(r1²+r2²+r3²)

其中，wt为土壤体积含水率对etc的权重；r为决定系数。

则大田玉米晴天的蒸发蒸腾量函数关系式为：

etc＝etc1×wt1+etc2×wt2+etc3×wt3；

大田玉米雨天蒸发蒸腾量与降水量pre的函数关系式为：

etc＝5.720×pre^-0.606r²＝0.878^**

其中，pre为实际测量00:00-24:00时间段的降水量(mm)。

具体地，所述大田玉米各分段土壤体积含水率数据为大田玉米种植区布置的智墒水分传感仪测得的数据，所述雨天降水量为天圻多要素自动气象站测定所得，用于计算晴天和雨天的etc。所述大田玉米全生育期的etc值为玉米整个生育期晴天和雨天所得etc的总和。本发明所述晴天包括晴天、多云、阴天等所有无降水天数；雨天指天圻气象站所能获取的有有效降水量的天数。

本发明还提供一种大田玉米蒸发蒸腾量计算系统，包括：测量模块301，用于获取所述大田玉米0-30cm各分段体积含水率和所述雨天的降水量；计算模块302，用于根据大田玉米各分段体积含水率和etc的函数、体积含水率对etc的权重、雨天降水量与etc的函数关系等计算出大田玉米全生育期各种天气情况的蒸发蒸腾速率(图2)。具体的大田玉米蒸发蒸腾量计算方法可参见上述实施例，本发明实施例在此不再赘述。

具体地，测量模块通过智墒水分传感器、天圻气象站分别获取0-30cm各分段土壤体积含水率数据、雨天有效降水量数据，计算模块通过权重累加等方法对测量的数据进行运算分析，最终实现了对大田玉米在各种天气状况下的蒸发蒸腾量的有效计算，能够大规模应用在实际生产中。

本发明提供的一种适用于黄淮南部的大田玉米蒸发蒸腾量计算方法及系统，通过对大田玉米0-30cm分段土壤体积含水率和雨天降水量进行实时无损检测，实现了对大田玉米蒸发蒸腾量的有效估计，测量结果精确，并且能够在实际生产中大规模应用。

以上显示和描述了本发明的基本原理，主要特征和优点，在不脱离本发明精神和范围的前提下，本发明还有各种变化和改进，这些变化和改进都落入要求保护的本发明的范围。