一种智能判定旱作作物缺水状态的方法及系统与流程

本发明涉及农作物监测领域，具体地说，涉及一种智能判定旱作作物缺水状态的方法及系统。

背景技术：

日常应用中，常用“缺水胁迫”来描述旱作作物在干旱条件下由于土壤水分含量过低导致不能提供足量水供应进而胁迫迟滞作物生长发育状态。由于缺乏合适的技术手段，日常农业管理中如何实时、快速评判作物缺水状态一直处于技术空缺状态。

技术实现要素：

本发明所要解决的技术问题是针对现有技术的不足，提供一种智能判定旱作作物缺水状态的方法及系统。

本发明解决上述技术问题的技术方案如下：

一种智能判定旱作作物缺水状态的方法，包括如下步骤：

s1.传感器管将实时采集的土壤水分数据传递通过信号传输装置传递至远程终端中；

s2.远程终端接收信号传输装置传递传输来的土壤水分数据并制作虚拟根系二维图；

s3.远程终端根据土壤水分数据并判别作物的缺水状态。

进一步，所述步骤s2包括：

s2-1.将水分传感器按照时间间隔分别实时采集各个感应点h的土壤水分含量记为sm_h；

s2-2.结合即时(t0)和上一个点(t0-1)的两组数据，在每一个时间点计算土壤水分变化率sd_h；

s2-3.对垂直维度作物根系范围内每一个测量点h的土壤水分时间变化率sd_h进行日积分,获得sd_h_dsum

s2-4.对sd_h_dsum在作物根系深度进行归一化处理，cm_h(h),h＝0,1,…l。结果表示出每一个垂直深度的土壤水分变化相对贡献率；

s2-5.依照cm_h(h)结果构建虚拟根系二维图。

进一步，所述步骤s2-1中的时间间隔为5-15分钟。

进一步，所述步骤s3包括：

测定在“作物枯萎点”时的土壤平均水分含量值sw，即为100％缺水胁迫阀值；在“田间持水量”时的土壤平均水分含量值sf，即为0％缺水胁迫(不胁迫)阀值；以及实测土壤水分含量sm，则作物缺水胁迫度pms通过公式：pms＝(sf-sm)/(sf-sw)*100％计算。

本发明的有益效果在于：利用高分辨率管式土壤水分传感器判断作物缺水状态，包括传感器布置方法、缺水状态判定依据和判定方式；利用“虚拟根系二维图”来定性判定作物缺水状态，包括“虚拟根系二维图”的具体制作方式；利用高分辨率管式土壤水分传感器测量的“土壤水分含量”、“土壤水分在垂直维度的分布”以及“土壤水分在垂直维度的最大变化率”来定量判断“作物缺水状态”，包括以“土壤水分在垂直维度的最大变化率”来定量判断作物缺水状态的具体方式。

需要说明的是，根系深度的获得可以根据经验数值获得，也可以根据下述方法获得，其中以本发明的传感器管为第一传感器管，在距离第一传感器管5-10cm处，设置第二传感器管，所述第二传感器管垂直设置在地面中，所述第二传感器管距离作物第一传感器的距离为5-20cm，所述第二传感器、第一传感器和主杆中心轴在同一平面上：

w1第一传感器管和第二传感器管将实时采集的数据传递通过信号传输装置传递至远程终端中；

w2远程终端根据高度和时间将数据分类和处理；

w3将处理得到的结果与之前的七天的数据结果相比对并进行平滑处理，得到根系深度数据。

所述步骤w2包括如下步骤：

w2-1.将远程终端接收到的数据按照其来源进行分类，由第一传感器管传递来的数据记为a，由第二传感器管传递来的数据记为b,其中每个传感器距离土壤表层的距离记为h(cm)，传感器测量数据记录的时间为主判断依据，将一天的时间分为：日照区段(日升至日落)dbin、黑夜区段(月照或全黑)nbin、过渡(介于白天和黑夜之间)区段gbin；

w2-2.从土壤表面往下，分别针对每个高度h判断：

x：日照区段在高度h上的土壤水分变化率的日平均值(a-dbin-h-dmr-dayaverage)和黑夜区段在高度h上的土壤水分变化率的日平均值(a-nbin-h-dmr-dayaverage)的差异；

y：日照区段在高度h上的土壤温度变化率的日平均值(a-dbin-h-dtr-dayaverage)和黑夜区段在高度h上的土壤温度变化率的日平均值(a-nbin-h-dtr-dayaverage)的差异；

z：日照区段在高度h上的水平水分变化率的日平均值(ba-dbin-h-hmr-dayaverage)和黑夜区段在高度h上的水平水分变化率的日平均值(ba-nbin-h-hmr-dayaverage)的差异；

w2-3.针对步骤w2-2中的判断结果，

若x有差异时：若z和y均有差异，此时记录h为新的根系深度，并替代上一个根系深度；若z和y中有一个有差异，则h为新的根系深度的几率为67％；若z和y均没有差异，则h为新的根系深度的几率为33％；对下个h的数据进行步骤w2-1和w2-2；

若x没有差异：判断h是否大于根系深度，否，对下个h的数据进行步骤w2-1和w2-2；是，表示已经超越根系深度，结束判定。

所述步骤w2-2还包括包括如下步骤：

w2-2-1.每个水分传感器采集的土壤水分数据记为sw，每个温度传感器采集的土壤温度数据记为st；

w2-2-2.分别构建实时的“土壤水分垂直分布”图形yhw和“土壤温度垂直分布”图形yht；

w2-2-3.针对每一个时间点采样收集的数据，计算：垂直水分变化率vmr(vmr＝dsw/dh，h为传感器感应环间隔距离)，水平水分变化率hmr(hmr＝dsw/dl，l为传两个感器感间隔距离)

w2-2-4.结合即时(t0)和上一个点(t0-1)的两组数据，在每一个时间点分别计算：土壤水分变化率dmr(dmr＝dsw/dt)、土壤温度变化率dtr(dtr＝dt/dt)。

w2-2-5.传感器测量数据记录的时间为主判断依据，将一天的时间分为：日照区段(日升至日落)dbin、黑夜区段(月照或全黑)nbin、过渡(介于白天和黑夜之间)区段gbin，把每一个时间点测量和计算数据(vmr、hmr、dmr、dtr)整合为一个数据子集(dataset)；

w2-2-6.以天为单位对每一个bin中累积的数据进行统计计算，计算“日平均值”dayaverage、“日均方差值”daystd。

本发明还提供一种智能判定旱作作物缺水状态的系统，包括：

水分传感器管，所述水分传感器管垂直设置在地面中，用于实时采集根系附近土壤中的水分数据信息；

信号传输装置，所述信号传输装置与水分传感器管相连接，所述信号传输装置用于将水分传感器管采集到的数据信息通过无线网络传送到远程终端；

远程终端，所述远程终端接收信号传输来的数据信息，并进行数据处理。

本发明的有益效果是：能够实现对根系周围土壤的水分的实时监控，从而为实现生长在非饱和土壤水分环境中各种旱作缺水性能提供数据支持。

进一步，所述水分传感器管距离作物主杆中心轴的距离为10-30cm。

进一步，所述水分传感器管包括多个垂直分布的水分传感器，相邻的两个水分传感器之间的距离为5cm。

进一步，所述无线网路为lorawan网络。

附图说明

图1为本发明的水分传感器管的位置及结构示意图；

图2为本发明的水分传感器管的位置俯视图；

图3为本发明的非缺水状态下的“虚拟根系”二维图形；

图4为本发明的缺水状态下的“虚拟根系”二维图形；

图5为特定状态下，作物枯萎线、实测土壤水分以及田间持水量三条曲线的关系图。

附图标记说明：

1.水分传感器管，2.水分传感器。

具体实施方式

以下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。

作物缺水状态判定原理

在不考虑降雨或灌溉补水导致的短期跃升情况下，土壤中水分垂直分布以及变化主要受蒸发(e，evaporation，土壤主导)和蒸腾(t，transpiration，植物主导)两个过程控制，在传统测量技术中，e和t的测量很难被区分，一般手段只能捆绑测量，以et统一表征。

蒸发过程e以及蒸发量主要受太阳辐射或大风等气象因素影响，而作物蒸腾t则主要受作物从根系、主杆到叶片的吸水、输水和挥发主导，受植物叶面气孔闭合影响。e和t都是以天(24小时)作为时间周期在整个作物生长发育期内变动重复，e对土壤表面的水分变化影响显著，而t则是对植物根系、特别是处于根系末端的须根周围的土壤水分变化影响明显。

因此，只要解决对土壤水分在垂直维度分布、及其在作物根系周围的精准测量技术障碍，区分获得原位e和t的状态成为可能，这里一个高分辨率的、多感应点传感的、垂直布置的土壤水分传感器就成为技术实现的关键。进一步地，通过对e和t的有效区分，测量结果也为评判作物缺水状态提供了有效技术保障。

针对一个特定的作物，在整个作物生长期,从土壤表面垂直向下直至作物根系深度范围，土壤中水分含量的绝对值、在垂直上的分布和相对时间变化率与土壤表面水分蒸发、作物根系的吸水速率、吸水量有一一对应的关系。垂直布置的土壤水分传感器对这几个被观测量进行高分辨率数据收集，通过进一步数学和统计学处理，可以对作物缺水状态进行客观的评价。

对垂直维度作物根系深度l范围内每一个测量点h，依据采集时间划分，计算出

土壤水分含量(sm_h)的时间变化率,sd_h

土壤水分含量点日平均,sm_dave_h

土壤水分含量点日均方差,sm_dstd_h

土壤水分含量点白昼(6:00-18:00)平均,sm_mave_h

土壤水分含量点白昼均方差,sm_mstd_h

土壤水分含量点黑夜(18:00-6:00)平均,sm_nave_h

土壤水分含量点黑夜均方差，sm_nstd_h。

对垂直维度作物根系范围内每一个测量点h的土壤水分时间变化率sd_h进行日积分(0:00-24:00),获得sd_h_dsum

然后在作物根系深度进行归一化处理，结果直观表示出每一个垂直深度的土壤水分变化相对贡献率，cm_h(h),h＝0,1,…l

本发明以“虚拟根系”二维图形式来表示cm_h(h)。

在土壤水分含量充沛状态下(根据土壤水分绝对含量和垂直分布初步判断)，作物生长所需要的水分可以通过根系从土壤中充分获取。虚拟根系图示中，作物根系活跃区(须根所在区域)对作物需水总量的贡献比例会占显著比例，整个虚拟根系二维图示呈现“上瘦下胖”葫芦形状。并且，在根系活跃区呈现昼夜交替变化(由作物气孔昼开夜闭特性导致)的特征，统计获得的点白昼平均sm_mave_h与点黑夜平均sm_nave_h之间存在显著差异。

在土壤干旱状态下(根据土壤水分绝对含量和垂直分布初步判断)，作物生长所需要的水分不能通过根系从土壤中充分获取。虚拟根系图示中，作物根系活跃区(须根所在区域)对作物需水总量的贡献比例很小，土壤表面由太阳照射或大风吹扫导致的蒸发占有比例显著，整个虚拟根系二维图示呈现“上宽下窄”雨伞形状，并且昼夜交替变化(由作物气孔昼开夜闭特性导致)特征被极大抑制，统计获得的点白昼平均sm_mave_h与点黑夜平均sm_nave_h之间不存在显著差异。

发明人定义，针对一个特定土壤、一个特定作物，在垂直维度作物根系范围，在“作物枯萎点”时的土壤平均水分含量值为100％缺水胁迫阀值sw、在“田间持水量”时的土壤平均水分含量值为0％缺水胁迫(不胁迫)阀值sf，介于两个阀值之间土壤水分含量sm可以被用来定量描述“作物缺水胁迫度”pms。

如图5所示，图5表示

在收集大量数据、不断深度学习的过程中萃取和不断迭代判定参数(作物枯萎点sw、田间持水量sf)集，总结出合适的、针对一个特定作物的“作物缺水状态(水胁迫度)”判定模式。

需要根据垂直维度土壤水分绝对含量和相对分布，对收集的数据进行清洗处理(预处理)，滤除调不适合导入到判定算法的数据。

实施例1

土壤水分传感器的采样频率定在5-15分钟/次，分别采集各个感应点h的土壤水分含量sm_h；

结合即时(t0)和上一个点(t0-1)的两组数据，在每一个时间点计算土壤水分变化率sd_h，sd_h是标量；

对垂直维度作物根系范围内每一个测量点h的土壤水分时间变化率sd_h进行日积分,获得sd_h_dsum；

对sd_h_dsum在作物根系深度进行归一化处理，cm_h(h),h＝0,1,…l。结果表示出每一个垂直深度的土壤水分变化相对贡献率；

依照cm_h(h)结果构建“虚拟根系”二维图；

进行“作物缺水状态”定性判断，当构建出的“虚拟根系”二维图为倒锥形时(如图4所示)，表示作物处于缺水状态，当构建出的“虚拟根系”二维图为葫芦形时(如图3所示)，表示作处于非缺水状态；

分别根据“作物枯萎线”计算作物缺水状态时土壤平均水分含量值sw、“田间持水量”时土壤平均水分含量值sf、以及实测土壤水分含量sm，再根据公式“pms＝(sf-sm)/(sf-sw)*100％计算”定量计算“作物缺水胁迫度”pms。

进行“作物缺水状态”定量判断。以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。