一种基于作物环境协同反馈的灌溉施肥系统

admin管理员组

文章数量:866648

一种基于作物环境协同反馈的灌溉施肥系统，该系统包括气象站、土壤湿度传感器、叶面积指数测定仪和控制器，气象站、土壤湿度传感器和叶面积指数测定仪获得与农作物生长有关的日照时数、气压、气温、空气湿度、降水量、小气候环境2米高风速、太阳辐射、田间持水量、田间持水率、含水饱和度、凋萎点土壤含水量、凋萎系数、叶面积指数参数后，传输给控制器进行数据处理与计算，得到入渗量、潜在入渗量、潜在作物蒸腾量、土壤表面潜在蒸发量和潜在蒸发量参数，进一步运算建立数学模型，形成灌溉计划。

1.一种基于作物环境协同反馈的灌溉施肥系统，其特征在于，该系统包括气象站、土壤湿度传感器、叶面积指数测定仪和控制器，气象站、土壤湿度传感器和叶面积指数测定仪获得与农作物生长有关的日照时数、气压、气温、空气湿度、降水量、小气候环境2米高风速、太阳辐射、田间持水量、田间持水率、含水饱和度、凋萎点土壤含水量、凋萎系数、叶面积指数参数后，传输给控制器进行数据处理与计算，得到入渗量、潜在入渗量、潜在作物蒸腾量、土壤表面潜在蒸发量和潜在蒸发量参数，进一步运算建立数学模型，形成灌溉计划。 2.如权利要求1所述的基于作物环境协同反馈的灌溉施肥系统，其特征在于，所述的数学模型包括以下处理步骤：A1，将土壤参数的单位由体积单位转换为高度单位，WP＝DP×WP×10.0(2)FC＝DP×FC×10.0(3)ST＝DP×ST×10.0(4)其中，WP、FC和ST分别是凋萎系数、田间持水量和饱和含水量；A2，计算土壤蓄积量(S)，根据土壤入渗和径流曲线得到，S＝254×(100/CN－1)(5)A3，计算土壤含水量阈值，低于该阈值发生干旱胁迫(THE)，该阈值的近似公式为，THE＝WP+0.75×(FC－WP)(6)A4，根据初始土壤含水量计算初始胁迫因子(SWFAC1和SWFAC2)，累计降雨量(TRAIN)，灌溉水量(TIRR)，土壤蒸发量(TESA)，作物蒸腾量(TEPA)，地面径流量(TROF)，垂直向排水量(TDRN)和入渗量(TINF)的累积值；A5，进行作物状态计算，读取灌溉量参数，潜在入渗量(POTINF)是降雨量(TRAIN)和灌水量(TIRR)的总和，根据每日土壤含水量高于田间持水量的数值来计算得出土壤垂直向排水量(DRN，以mm为单位)，公式为，DRN＝(SWC－FC)×DRN(7)如果潜在入渗量(POTINF)大于零，则使用SCS曲线计算每日地表水径流量(ROF)；入渗量(INF)是潜在入渗与径流之间的差，基于Priestly-Taylor方法计算每日潜在的蒸腾量(ETp)，基于LAI的表面反射率(ALB)由土壤和作物反射率的加权平均值来估计，与叶面指数测定仪数据对比进行模型参数修正：ALP＝0.1×EXP(－0.7×LAI)+0.2×(1－EXP(－0.7×LAI))(9)由式(10)和(11)计算白天的平均温度(T)和平衡蒸发(EEQ)，Tmed＝0.6×T+0.4×T(10)EEQ＝SRAD×(4.88E－3－4.37E－3×ALB)×(Tmed+29)(11)通过系数(f)来调节平衡蒸发速率，得到ETp的最终值，接下来，用与反射率相同的加权系数来计算潜在的土壤蒸发量(ESp)和作物腾发速率(EPp)：ES＝ET×EXP(－0.7×LAI)(12)EP＝ET×(1－EXP(－0.7×LAI)(13)根据当前土壤水分可用量大小计算实际的每日土壤蒸发量(ESa)，如果土壤含水量低于凋萎系数，则不会发生水分蒸发，如果土壤含水量大于田间持水量，则会发生蒸发；土壤含水量在凋萎系数和田间持水量之间，实际蒸腾量在0和潜在蒸腾量之间呈线性变化，潜在的作物蒸腾速率(EPp)被最小的土壤水分胁迫因子(SWFAC1或SWFAC2)协导致其数值降低从而获得实际作物蒸腾量(EPa)；A6，进行数据的汇总集成，基于入渗量(INF)、土壤蒸发量(ESa)、作物蒸腾(EPa)和垂直向排水(DRN)的值运算土壤水分平衡模块的汇总集成，更新土壤含水量(SWC)的值，公式为：SWC＝SWC+(INF-EPa-DRN)(15)土壤含水量计算值被限制在饱和含量的最大值到零之间，如果计算土壤含水量超过土壤饱和含水量，则应当调整径流量和土壤含水量。如果计算的土壤含水量小于零，则要考虑额外的调整因子(SWC_ADJ)，最后是累积入渗量、蒸发量、蒸腾量、排水量和径流量的更新；A7，计算土壤水分胁迫因子，干旱胁迫因子(SWFAC1)为1.0时发生最小胁迫，如果土壤含水量大于“初始化”部分(THE)中计算的阈值，在凋萎点之下，发生最大胁迫时(SWFAC1＝0.0)，在这范围之间的胁迫因子是线性分布的；土壤剖面底部地下水的埋深(WTABLE，以mm为单位)从超过田间持水量可用的过量水来计算，然后计算地下水的水层厚度(DWT，以mm为单位)，WTABLE＝(SWC-FC)/(ST-FC)×DP×10.0(17)DWT＝DP×10.0－WTABLE(18)当水位下降时，发生土壤最小多余水分胁迫(SWFAC2＝1.0)为零。当地下水位(DWT)的值大于250mm(STRESS_DEPTH＝250mm)时，发生最大多余水分胁迫(SWFAC2＝0.0)，其余的水分胁迫因子在这些水位之间进行线性条件插值。 3.如权利要求2所述的基于作物环境协同反馈的灌溉施肥系统，其特征在于，输入数学模型的气象参数时，在初始阶段按照干旱年、正常年和湿润年三种典型年的全年月数据生成逐日输入模拟量。在实际运行时，根据真实出现的气象数据修正，同时根据3天、5天天气预报数值进行预测。 4.如权利要求2所述的基于作物环境协同反馈的灌溉施肥系统，其特征在于，在灌水决策后如遇到降雨，降雨大于1mm时，停止灌水，同时记录降雨量，在土壤水分平衡模块中对土壤含水量进行相应的扣除，若作物生长环境为温室、大棚，其降雨量、小气候环境2米高风速均视为0。 5.如权利要求1所述的基于作物环境协同反馈的灌溉施肥系统，其特征在于，连通水源的进水管的管道首先经过主管过滤器，在经过水泵后依次连通主管电磁阀A、主管电磁阀B和主管流量传感器，主管过滤器过滤灌溉水源，水泵提供水流动力，主管电磁阀A和主管电磁阀B控制灌水管路的开启和闭合，主管流量传感器测得主管流量，在进水管的主管电磁阀A和主管电磁阀B之间，进水管连通施肥机支管A，在施肥机支管A上设置有支管电磁阀A，在进水管的主管电磁阀B和主管流量传感器之间，进水管连通施肥机支管B，在施肥机支管B上设置有支管电磁阀B，施肥机支管A和施肥机支管B经过本地操作端装置连接肥料桶，所述控制器被置于本地操作端装置内，本地操作端装置的底座上设有包括施肥机进水口、注肥管、肥液过滤器、施肥流量传感器、计量泵、水泵、压力变送器、EC测量仪、持压阀和出水管，肥液流动顺序为注肥管、肥液过滤器、施肥流量传感器、计量泵，混肥管道、水泵、出水管当控制器打开主管电磁阀A、主管电磁阀B，关闭支管电磁阀A、支管电磁阀B，设备只进行灌水，当控制器打开主管电磁阀A、支管电磁阀A、支管电磁阀B，关闭主管电磁阀B，进行施肥。

技术领域

本发明属于现代农业技术领域，特别涉及一种基于作物环境协同反馈的灌溉施肥系统。

背景技术

公开号为CN102165876A的专利文献公开了“一种智能节水浇灌施肥系统，其特征在于，所述系统包括：管理控制单元、核心处理单元、通信单元、土壤环境检测单元、用水检测单元、施肥控制单元以及浇灌控制单元；其中，所述管理控制单元通过所述通信单元接收管理人员的远程/本地指令，并转发给所述核心处理单元；所述核心处理单元自动或响应所述指令向所述环境检测单元和/或用水检测单元发送状态参数采集指令，根据所述状态参数和预先存储的历史统计数据生成浇灌施肥决策信息；根据管理人员的控制指令或所述决策信息发送施肥/浇灌指令；以及，根据需要将所接收的信息转发给所述管理控制单元；所述土壤环境检测单元实时或响应所述参数采集指令获取与待浇灌场地的土壤环境有关的信息，包括土壤养分信息，并发送给所述核心处理单元；所述用水检测单元实时或响应所述参数采集指令获取与浇灌用水的水质和浇灌管道的水压有关的信息，并发送给所述核心处理单元；所述施肥控制单元响应所述施肥指令，控制溶入浇灌用水中的肥料类型和用量；接受施肥设备的指令反馈信息并发送给所述核心处理单元；所述浇灌控制单元响应所述浇灌指令，对浇灌施肥用水的水质进行过滤处理；对浇灌管道的水压和/或出水量进行调节；接受浇灌设备的指令反馈信息并发送给所述核心处理单元。”上述文件对于作物的灌溉施肥提供了一个较为完整的解决方案，但是，由于现代精细农业对于灌溉施肥的控制要求越来越高，现有技术对于智能化的灌溉施肥要求还是无法满足。

发明内容

本发明提供的一种基于作物环境协同反馈的灌溉施肥系统，目的在于解决现有农作物灌溉施肥系统在控制算法上不能满足越来越精细的要求。

一种基于作物环境协同反馈的灌溉施肥系统，该系统包括气象站、土壤湿度传感器、叶面积指数测定仪和控制器，气象站、土壤湿度传感器和叶面积指数测定仪获得与农作物生长有关的日照时数、气压、气温、空气湿度、降水量、小气候环境2米高风速、太阳辐射、田间持水量、田间持水率、含水饱和度、凋萎点土壤含水量、凋萎系数、叶面积指数参数后，传输给控制器进行数据处理与计算，得到入渗量、潜在入渗量、潜在作物蒸腾量、土壤表面潜在蒸发量和潜在蒸发量参数，进一步运算建立数学模型，形成灌溉计划。

所述的数学模型包括以下处理步骤：

A1，将土壤参数的单位由体积单位转换为高度单位，

WP＝DP×WPP×10.0 (2)

FC＝DP×FCP×10.0 (3)

ST＝DP×STP×10.0 (4)

其中，WP、FC和ST分别是凋萎系数、田间持水量和饱和含水量；

A2，计算土壤蓄积量(S)，根据土壤入渗和径流曲线得到，

S＝254×(100/CN-1) (5)

A3，计算土壤含水量阈值，低于该阈值发生干旱胁迫(THE)，该阈值的近似公式为，

THE＝WP+0.75×(FC-WP) (6)

A4，根据初始土壤含水量计算初始胁迫因子(SWFAC1和SWFAC2)，累计降雨量(TRAIN)，灌溉水量(TIRR)，土壤蒸发量(TESA)，作物蒸腾量(TEPA)，地面径流量(TROF)，垂直向排水量(TDRN)和入渗量(TINF)的累积值。

A5，进行作物状态计算，

读取灌溉量参数，潜在入渗量(POTINF)是降雨量(TRAIN)和灌水量(TIRR)的总和，根据每日土壤含水量高于田间持水量的数值来计算得出土壤垂直向排水量(DRN，以mm为单位)，公式为，

DRN＝(SWC-FC)×DRNP (7)

如果潜在入渗量(POTINF)大于零，则使用SCS曲线计算每日地表水径流量(ROF)；

入渗量(INF)是潜在入渗与径流之间的差，

基于Priestly-Taylor方法计算每日潜在的蒸腾量(ETp)，

基于LAI的表面反射率(ALB)由土壤和作物反射率的加权平均值来估计，与叶面指数测定仪数据对比进行模型参数修正：

ALP＝0.1×EXP(-0.7×LAI)+0.2×(1-EXP(-0.7×LAI)) (9)

由式(10)和(11)计算白天的平均温度(T)和平衡蒸发(EEQ)，

Tmed＝0.6×TMAX+0.4×TMIN (10)

EEQ＝SRAD×(4.88E-3-4.37E-3×ALB)×(Tmed+29) (11)

通过平衡蒸发速率，得到ETp的最终值，

接下来，用与反射率相同的加权系数来计算潜在的土壤蒸发量(ESp)和作物腾发速率(EPp)：

ESP＝ETP×EXP(-0.7×LAI) (12)

EPP＝ETP×(1-EXP(-0.7×LAI) (13)

根据当前土壤水分可用量大小计算实际的每日土壤蒸发量(ESa)，如果土壤含水量低于凋萎系数，则不会发生水分蒸发，如果土壤含水量大于田间持水量，则会发生蒸发；

土壤含水量在凋萎系数和田间持水量之间，实际蒸腾量在0和潜在蒸腾量之间呈线性变化，

潜在的作物蒸腾速率(EPp)被最小的土壤水分胁迫因子(SWFAC1或SWFAC2)协导致其数值降低从而获得实际作物蒸腾量(EPa)；

A6，进行数据的汇总集成，

基于入渗量(INF)、土壤蒸发量(ESa)、作物蒸腾(EPa)和垂直向排水(DRN)的值运算土壤水分平衡模块的汇总集成，更新土壤含水量(SWC)的值，公式为：

SWC＝SWC+(INF-EPa-DRN) (15)

土壤含水量计算值被限制在饱和含量的最大值到零之间，如果计算土壤含水量超过土壤饱和含水量，则应当调整径流量和土壤含水量。如果计算的土壤含水量小于零，则要考虑额外的调整因子(SWC_ADJ)，

最后是累积入渗量、蒸发量、蒸腾量、排水量和径流量的更新；

A7，计算土壤水分胁迫因子，

干旱胁迫因子(SWFAC1)为1.0时发生最小胁迫，如果土壤含水量大于“初始化”部分(THE)中计算的阈值，在凋萎点之下，发生最大胁迫时(SWFAC1＝0.0)，在这范围之间的胁迫因子是线性分布的；

土壤剖面底部地下水的埋深(WTABLE，以mm为单位)从超过田间持水量可用的过量水来计算，然后计算地下水的水层厚度(DWT，以mm为单位)，

WTABLE＝(SWC-FC)/(ST-FC)×DP×10.0 (17)

DWT＝DP×10.0-WTABLE (18)

当水位下降时，发生土壤最小多余水分胁迫(SWFAC2＝1.0)为零。当地下水位(DWT)的值大于250mm(STRESS_DEPTH＝250mm)时，发生最大多余水分胁迫(SWFAC2＝0.0)，其余的水分胁迫因子在这些水位之间进行线性条件插值。

进一步的，在输入数学模型的气象参数时，在初始阶段按照干旱年、正常年和湿润年三种典型年的全年月数据生成逐日输入模拟量。在实际运行时，根据真实出现的气象数据修正，同时根据3天、5天天气预报数值进行预测。

在灌水决策后如遇到降雨，降雨大于1mm时，停止灌水，同时记录降雨量，在土壤水分平衡模块中对土壤含水量进行相应的扣除，

若作物生长环境为温室、大棚，其降雨量、小气候环境2米高风速均视为0。

本发明的基于作物环境协同反馈的灌溉施肥系统包括，连通水源的进水管的管道首先经过主管过滤器，在经过水泵后依次连通主管电磁阀A、主管电磁阀B和主管流量传感器，

主管过滤器过滤灌溉水源，水泵提供水流动力，主管电磁阀A和主管电磁阀B控制灌水管路的开启和闭合，主管流量传感器测得主管流量，

在进水管的主管电磁阀A和主管电磁阀B之间，进水管连通施肥机支管A，在施肥机支管A上设置有支管电磁阀A，

在进水管的主管电磁阀B和主管流量传感器之间，进水管连通施肥机支管B，在施肥机支管B上设置有支管电磁阀B，

施肥机支管A和施肥机支管B经过本地操作端装置连接肥料桶，

所述控制器被置于本地操作端装置内，

本地操作端装置的底座上设有包括施肥机进水口、注肥管、肥液过滤器、施肥流量传感器、计量泵、吸肥泵、压力变送器、EC测量仪、持压阀和出水管，肥液流动顺序为注肥管、肥液过滤器、施肥流量传感器、计量泵、混肥管道、水泵、出水管

当控制器打开主管电磁阀A、主管电磁阀B，关闭支管电磁阀A、支管电磁阀B，该设备不进行施肥作业，只进行进行灌水作业，执行灌溉计划。

当控制器打开主管电磁阀A、支管电磁阀A、支管电磁阀B，关闭主管电磁阀B，主管道的灌溉水进入施肥设备，该设备进行肥料配比和肥料添加，混合好的水肥进入灌溉主管道，流入田间进行灌溉施肥作业。

本发明可以分为信息感知模块、决策模块和执行模块。所述的信息感知模块包括小型气象站、土壤湿度传感器、叶面积指数测定仪、主管流量传感器、施肥流量传感器、EC值测量仪、压力变送器，决策模块包括控制器，执行模块包括主管电磁阀A、主管电磁阀B、支管电磁阀A、支管电磁阀B、计量泵、水泵、吸肥泵。针对灌溉计划所需，通过各种传感器采集数据，传输给控制器进行数据处理与运算，进一步建立数据模型，形成灌溉计划。控制器根据灌溉计划调控执行模块，进行精准灌溉施肥作业，同时信息感知模块实时采集参数反馈至控制器，控制器根据灌溉计划再次调控执行模块，形成感知、决策、执行于一体的协调反馈智慧灌溉施肥循环系统。

本发明的有益效果：

1)本发明是基于作物环境协同反馈的智慧灌溉施肥装置系统，能够检测大量作物环境参数，并传输到本地控制器进行系统地运算，为灌溉计划的数学模型分析提供数据支撑；

2)本发明是基于作物环境协同反馈的智慧灌溉施肥装置系统，灌溉计划包含一种估算田间灌溉水量模型的方法，能够系统地形成每日灌溉用水需求数据库，实现精准灌溉；

3)本发明是基于作物环境协同反馈的智慧灌溉施肥装置系统，可以根据已经形成的灌溉计划，自动开始或停止灌溉施肥，灌溉过程无需人为监控；

4)本发明是基于作物环境协同反馈的智慧灌溉施肥装置系统，能够自动采集灌溉施肥后的土壤、小气候环境等数据，动态监测水流压力值和肥液EC值，并自行反馈至本地控制器，重新生成灌溉计划。

附图说明

通过参考附图阅读下文的详细描述，本发明示例性实施方式的上述以及其他目的、特征和优点将变得易于理解。在附图中，以示例性而非限制性的方式示出了本发明的若干实施方式，其中：

图1为发明的智慧灌溉施肥系统组成示意图；

图2为发明的施肥装置立体示意图图；

图3为本发明的施肥装置的侧视图；

图4为本发明智慧灌溉施肥系统控制器外形图。

图5是本发明实施例中EC值与KCl溶液浓度关系图。

1-小型气象站、2-土壤温湿度传感器、3-叶面积指数测定仪、4-进水管、5-主管过滤器、6-水泵、7-主管电磁阀A、8-主管电磁阀B、9-支管、10-支管电磁阀A、11-支管电磁阀B、12-主管流量传感器、13-施肥机进水口、14-注肥管、15-肥液过滤器、16-施肥流量传感器、17-计量泵、18-水泵、19-压力变送器、20-EC测量仪、21-持压阀、22-出水管、23-底座、24-肥料桶、25-控制器、26-施肥装置27-远程操作端。

具体实施方式

下面结合附图和具体实施例对本发明做进一步的描述。

实施例1：接通电源，小型气象站、土壤湿度传感器、叶面积指数测定仪正常工作，灌溉前，各传感器收集数据反馈至控制器，控制器根据程序运算得出灌溉计划。只需灌水时，控制器根据灌溉计划，打开主管电磁阀A、主管电磁阀B、水泵，关闭支管电磁阀A、支管电磁阀B、施肥泵，水流通过进水管、主管过滤器、水泵、主管电磁阀A、主管电磁阀B、主管流主管流量计，到达田间。灌溉后，各传感器再次收集数据反馈至控制器，控制器根据程序运算更新灌溉计划，指导下一次灌溉。

实施例2：接通电源，小型气象站、土壤湿度传感器、叶面积指数测定仪正常工作，假设土层为相对不可渗透层下层的单个均质土层，进行水平衡计算，更新每天的土壤含水量变化。

在测得日照时数、气压、气温、空气湿度、降水量、小气候环境2米高风速、太阳辐射、田间持水量、田间持水率、含水饱和度、凋萎点土壤含水量、凋萎系数、叶面积指数等参数后，传输给控制器进行数据处理与计算，得到入渗量、潜在入渗量、潜在作物蒸腾量、土壤表面潜在蒸发量、潜在蒸发量等，进一步运算建立数据模型，形成灌溉计划。

灌溉计划是一种利用田间作物耗水数学模型实现估算田间灌溉水量的方法，其包括以下步骤：

叶面积指数影响作物灌溉综合系数kc由叶面指数测定仪确定。作物的ETp：

ETp＝kc×ET0 (1)

程序采用数学模型计算分析后，形成以日为单位的作物需要水量数据库。

进一步的，可在本地操作端或远程操作端人为设置参数指标修改灌溉计划。

本发明包括小型气象站、土壤湿度传感器、叶面积指数测定仪，自动采集灌溉施肥后的作物生长状态数据，反馈到控制器重新计算分析，形成新的灌溉计划。

本发明中的压力变送器、EC测量仪、主管流主管流量计、施肥流量计动态监测水流压力值、肥液EC值，过水流量，数据反馈至控制器计算分析，根据设定参数，对执行模块进行调控。

灌溉计划形成流程分为模型初始化、状态计算和汇总集成。

模型初始化。代码的“初始化”部分主要是输入数据文件SOIL.INP和IRRIG.INP和打开输出文件SW.OUT。读取所需的环境数据输入SOIL.INP文件然后关闭，变量名称写入输出文件SW.OUT中。土壤参数的单位由体积转换为高度(mm)，公式如下：

WP＝DP×WPP×10.0 (2)

FC＝DP×FCP×10.0 (3)

ST＝DP×STP×10.0 (4)

其中WP、FC和ST分别是凋萎系数、田间持水量和饱和含水量。RUNOFF子程序计算土壤蓄积量(S)，根据土壤入渗和径流曲线得到：

S＝254×(100/CN-1) (5)

调用子程序STRESS来计算土壤含水量阈值，低于此值发生干旱胁迫(THE)，近似公式为：

THE＝WP+0.75×(FC-WP) (6)

然后根据初始土壤含水量计算初始胁迫因子(SWFAC1和SWFAC2)，累计降雨量(TRAIN)，灌溉水量(TIRR)，土壤蒸发量(TESA)，作物蒸腾量(TEPA)，地面径流量(TROF)，垂直向排水量(TDRN)和入渗量(TINF)的累积值，开始模拟初始时假定为零。

状态计算。灌溉量从文件IRRIG.INP中读取，潜在入渗量(POTINF)是降雨量(TRAIN)和灌水量(TIRR)的总和。运用DRAINE子程序，根据每日土壤含水量高于田间持水量的数值来计算得出土壤垂直向排水量(DRN，以mm为单位)，公式为：

DRN＝(SWC-FC)×DRNP (7)

如果潜在入渗量(POTINF)大于零，则调用子程序RUNOFF，使用SCS曲线计算每日地表水径流量(ROF)。

IF(POTINF.GT.0.2＊S)THEN

ROF＝((POTINF-0.2＊S)＊＊2)/(POTINF+0.8＊S)

ELSE

ROF＝0

ENDIF (8)

入渗量(INF)是潜在入渗与径流之间的差。子程序ETpS基于Priestly-Taylor方法计算每日潜在的蒸腾量(ETp)。基于LAI的表面反射率(ALB)由土壤(0.1)和作物(0.2)反射率的加权平均值来估计，与叶面指数测定仪数据对比进行模型参数修正：

ALP＝0.1×EXP(-0.7×LAI)+0.2×(1-EXP(-0.7×LAI)) (9)

由下式计算白天的平均温度(T)和平衡蒸发(EEQ)。

Tmed＝0.6×TMAx+0.4×TMIN (10)

EEQ＝SRAD×(4.88E-3-4.37E-3×ALB)×(Tmed+29) (11)

通过系数(f)来调节平衡蒸发速率，得到ETp的最终值。

接下来，用与反射率相同的加权系数来计算潜在的土壤蒸发量(ESp)和作物腾发速率(EPp)：

ESP＝ETP×EXP(-0.7×LAI) (12)

EPP＝ETP×(1-EXP(-0.7×LAI) (13)

子程序ESAS根据当前土壤水分可用量大小计算实际的每日土壤蒸发量(ESa)。如果土壤含水量低于凋萎系数，则不会发生水分蒸发，如果土壤含水量大于田间持水量，则会发生蒸发。土壤含水量在凋萎系数和田间持水量之间，实际蒸腾量在0和潜在蒸腾量之间呈线性变化。

IF(SWC.LT.WP)THEN

a＝0

ELSEIF(SWC.GT.FC)THEN

a＝1

ELSE

a＝(SWC-WP)/(FC-WP)

ENDIF

ESa＝ESp＊a (14)

潜在的作物蒸腾速率(EPp)被最小的土壤水分胁迫因子(SWFAC1或SWFAC2)协导致其数值降低从而获得实际作物蒸腾量(EPa)。

汇总集成。基于入渗量(INF)、土壤蒸发量(ESa)、作物蒸腾(EPa)和垂直向排水(DRN)的值运算土壤水分平衡模块的汇总集成，更新土壤含水量(SWC)的值，公式为：

SWC＝SWC+(INF-EPa-DRN) (15)

土壤含水量计算值被限制在饱和含量的最大值到零之间。如果计算土壤含水量超过土壤饱和含水量，则应当调整径流量和土壤含水量。如果计算的土壤含水量小于零，则要考虑额外的调整因子(SWC_ADJ)。

最后是累积入渗量、蒸发量、蒸腾量、排水量和径流量的更新。

根据更新的土壤含水量值应用子程序STRESS计算土壤水分胁迫因子。干旱胁迫因子(SWFAC1)为1.0(最小胁迫)，如果土壤含水量大于“初始化”部分(THE)中计算的阈值。在凋萎点之下，发生最大胁迫时(SWFAC1＝0.0)。在这范围之间的胁迫因子是线性分布的。

IF(SWC.LT.WP)THEN

SWFAC1＝0.0

ELSEIF(SWC.GT.THE)THEN

SWFAC1＝1.0

ELSE

SWFAC1＝(SWC--WP)/(THE--WP)

SWFAC1＝MAX(MIN(SWFAC1，1.0)，0.0)

ENDIF (16)

土壤剖面底部地下水的埋深(WTABLE，以mm为单位)从超过田间持水量可用的过量水来计算。然后计算地下水的水层厚度(DWT，以mm为单位)。

WTABLE＝(SWC-FC)/(ST-FC)×DP×10.0 (17)

DWT＝DP×10.0-WTABLE (18)

当水位下降时，发生土壤最小多余水分胁迫(SWFAC2＝1.0)为零。当地下水位(DWT)的值大于250mm(STRESS_DEPTH＝250mm)时，发生最大多余水分胁迫(SWFAC2＝0.0)，其余的水分胁迫因子在这些水位之间进行线性条件插值。

IF(DWT.GE.STRESS_DEPTH)THEN

SWFAC2＝1.0

ELSE

SWFAC2＝DWT/STRESS_DEPTH

ENDIF (19)

以上是在作物灌水计划模拟分析模型。

输入数学模型的气象参数时，在初始阶段按照干旱年、正常年和湿润年三种典型年的全年月数据生成逐日输入模拟量。在实际运行时，根据真实出现的气象数据修正，同时根据3天、5天天气预报数值进行预测。

在灌水决策后如遇到降雨，降雨大于1mm时，停止灌水，同时记录降雨量，在土壤水分平衡模块中对土壤含水量进行相应的扣除。

若作物生长环境为温室、大棚，其降雨量、小气候环境2米高风速均视为0。

实施例3：需要施肥时，控制器形成灌溉计划后，自动关闭主管电磁阀B，打开主管电磁阀A、支管电磁阀A、支管电磁阀B、水泵、施肥泵，调控灌水量和施肥量。

流量计测得主管与输支管的流量Q1和Q2，Q2与Q1的比值与肥液浓度存在必然关系。以氯化钾(KCL)水溶肥为例，假设温度为25℃，EC值与KCl溶液浓度关系如图5所示。参考图5，控制器根据EC测量仪反馈的EC值，在保证Q1不变的情况下，通过调控计量泵来调节Q2的大小，或者在保证Q2不变的情况下，调节Q1的大小，即调节肥液浓度大小，从而达到调节灌溉肥液EC值的目的。

当EC值过大时，保持主管流量Q2不变，在本地操作端或远程操作端重新设置参数，控制器很具设定参数自动调控计量泵，减小Q1，肥液浓度降低，直至EC值符合灌溉计划设定值。

灌溉施肥结束时，本地控制器自动切断电源，停止水泵工作，关闭主管电磁阀和施肥电磁阀。

值得说明的是，虽然前述内容已经参考若干具体实施方式描述了本发明创造的精神和原理，但是应该理解，本发明并不限于所公开的具体实施方式，对各方面的划分也不意味着这些方面中的特征不能组合，这种划分仅是为了表述的方便。本发明旨在涵盖所附权利要求的精神和范围内所包括的各种修改和等同布置。

本文标签： 作物反馈环境系统