用于盐碱地的灌溉控制方法与系统

摘要

本发明揭示了一种用于盐碱地的灌溉控制方法与系统。该方法包括：设置参数，监测根区各土层土壤含水量，基于相对根长密度分布计算根区加权平均土壤含水量和根区加权平均土壤水渗透势，并计算作物水分胁迫指数CWSI，当CWSI计算值大于预定的CWSI临界值时，开始灌水；还可以根据计划湿润层内的土壤含水量分布计算理论灌水量，当实际灌水量达到理论灌水量时，停止灌水。本发明综合考虑了土壤水分、盐分以及根系分布对作物蒸腾与生长的影响，能更加准确、便捷地估算盐碱地作物所受到的水分胁迫程度，可为盐碱地农业节水、增产目标的实现提供有效工具。

权利要求

1. 一种用于盐碱地的灌溉控制方法，其特征在于，该方法包括： 将灌溉区域的土壤从表层至最大扎根深度L处依次分为若干层，层数记为k，用每层土 壤的平均深度Zi除以Lp得到每层土壤的相对深度，用每层土壤的厚度hi除以Lp得到 每层土壤的相对厚度Λζ#设定作物在每层土壤中的相对根长密度L md(zJ、萎焉系数θ"、 饱和土壤含水量9S、最适宜作物生长的土壤含水量上限ΘΗ与下限V、土壤水渗透势上限 外/与下限攸； 测量每层土壤的含水量、温度以及土体电导率； 计算根区加权平均土壤含水量，且当ΘΗ< 03或 9"时，令等于θ";当彡011时，令等于根据每层土壤的含水 量Θ i、温度和土体电导率获得各土层的土壤水渗透势，，并基于相对根长密度分布计算根 区加权平均土壤水渗透势f:，且当|贫| < w |时，令物等于炉/#; 当|#ι| > I#/, |时，令外·等于，； 计算作物水分胁迫指数CWSI: 〇fS/ = l-+ + 其中取/9 时，令3等于取 UL < C7 S C7¥ 0 且当Λ z 或时，令=等于θ«，当 ΘΗ < Θ < u s Θ s θ¥· 0等于当问之|外丨时'令；等于料； 当计算得到的CWSI大于预定的CWSI临界值时，开始灌水。

2. 如权利要求1所述的用于盐碱地的灌溉控制方法，其特征在于， 所述测量每层土壤的含水量Θ i、温度和土体电导率，是将土壤水分-温度-电导率探 针垂直插入根区土壤，并在土壤水分-温度-电导率探针上对应的各土层深度处设有一个 土壤水分-温度-电导率传感器，从而测得每层土壤的含水量Θ i、温度和土体电导率。

3. 如权利要求2所述的用于盐碱地的灌溉控制方法，其特征在于， 所述测得每层土壤的含水量Θ i、温度和土体电导率，是将分布在灌溉区域内的多根土 壤水分-温度-电导率探针上相同土层深度处的土壤水分-温度-电导率传感器的测量值 取平均而获得。

4. 如权利要求1所述的用于盐碱地的灌溉控制方法，其特征在于，该方法还包括： 在所述最大扎根深度范围内设定计划湿润层深度，记录从表层到计划湿润层深度的土 壤层数n，设定田间持水量9f、土壤灌溉湿润比p、田间水分有效利用系数η，淋洗系数R和 灌溉面积A，按下式计算单位面积的灌水定额Μ :然后用Μ乘 以Α，得到理论灌水量； 当实际灌水量达到理论灌水量时，停止灌水。

5. 如权利要求1至4其中任一项所述的用于盐碱地的灌溉控制方法，其特征在于，该方 法还包括： 当计算得到的CWSI大于预定的临界值时，若预报在设定的时间区间内会有降雨，则不 灌水，否则开始灌水。

6. -种用于盐碱地的灌溉控制系统，其特征在于，该系统包括参数输入子系统、土壤水 盐动态监测子系统、中央决策子系统和灌水控制子系统： 参数输入子系统，用于设定作物的最大扎根深度k，灌溉区域内从土壤表层至最大扎 根深度L处所分的层数k，每层土壤的平均深度Zi以及厚度lv作物在每层土壤中的相对根 长密度Lmd(zJ、萎焉系数θ"、饱和土壤含水量0 S、最适宜作物生长的土壤含水量上限ΘΗ 与下限θ ρ土壤水渗透势上限与下限抖，以及CWSI临界值； 土壤水盐监测子系统，用于测量每层土壤的含水量Θ i、温度以及土体电导率，并发送 给中央决策子系统； 中央决策子系统，包括计算模块和决策模块； 计算模块，用Zi除以L，得到每层土壤的相对深度Zh ;用比除以b得到每层土壤 的相对厚度Λ ^ ;然后计算根区加权平均土壤含水量且当ΘΗ< 03或0"时，令等于0W，当θ,< 0 11时，令01等 于V;根据每层土壤的含水量、温度和土体电导率获得各土层的土壤水渗透势，， 基于相对根长密度分布计算根区加权平均土壤水渗透势且当Μ <|供"|时，令f I.等于fw;当|終|乏时， 令灼等于ft;计算作物水分胁迫指 数 CWSI : £縦=1 - (M + /))(印 + f/)，其中取，且 当· < 3 A或^ S秦时，令5等于Θ w，当^ < ^ S ##时，令^等于Θ ;取

1=1 = ，且当免< |a|时，令：等于_ .当#ϊ $ |f^|时，令：等 I ! 11 φ Ψη·> 1 1 φ 于fi; 决策模块，比较计算得到的CWSI是否大于预定的CWSI临界值，当大于时，则向灌水控 制子系统发出开始灌水的指令； 灌水控制子系统，包括灌溉控制模块，当收到开始灌水的指令时，控制阀门开启而开始 灌水。

7. 如权利要求6所述的用于盐碱地的灌溉控制系统，其特征在于， 所述土壤水盐动态监测子系统，包括土壤水分-温度-电导率探针，土壤水分-温 度-电导率探针上对应每层土壤设有一个土壤水分-温度-电导率传感器，通过将土壤水 分-温度-电导率探针垂直插入根区土壤，从而测得每层土壤的含水量Θ i、温度和土体电 导率。

8. 如权利要求7所述的用于盐碱地的灌溉控制系统，其特征在于， 在灌溉控制区域内，所述土壤水分-温度-电导率探针为多根，所述测得每层土壤的含 水量Θ i，是将灌溉区域内的多根土壤水分-温度-电导率探针上相同土层深度处的土壤水 分-温度-电导率传感器的测量值取平均而获得。

9. 如权利要求6所述的用于盐碱地的灌溉控制系统，其特征在于， 参数输入子系统，还用于输入从表层到计划湿润层深度的土壤层数n，并设定田间持水 量9f、土壤灌溉湿润比p、田间水分有效利用系数η，淋洗系数R和灌溉面积A ; 中央决策子系统的计算模块还计算单位面积的灌水定额Μ :然后用Μ乘以Α，得到理论灌水量； 中央决策子系统的决策模块，还包括比较实际灌水量与理论灌水量，当实际灌水量达 到理论灌水量时，则向灌水控制子系统发出停止灌水的指令； 灌水控制子系统，还包括灌水计量模块，用于计量实际灌水量并发送给中央决策子系 统；灌溉控制模块还用于在收到停止灌水的指令时，关闭控制阀门而停止灌水。

10. 如权利要求6-9其中任一项所述的用于盐碱地的灌溉控制系统，其特征在于， 中央决策子系统还包括天气预报模块，用于接收天气预报，当预报在设定的时间区间 内有降雨时，则决策模块不发出开始灌水的指令。

说明书

用于盐碱地的灌溉控制方法与系统 技术领域

[〇〇〇1] 本发明涉及自动控制技术，更具体涉及用于盐碱地的灌溉控制方法与系统。 背景技术

[0002] 我国盐碱土面积约5. 2亿亩，开发利用潜力巨大。对于盐碱地作物，除土壤水分 夕卜，土壤盐分也是影响其水分吸收、生长以及产量的重要因素。土壤含盐量越高，土壤溶液 盐分浓度就越高，土壤水渗透势就越低，从而降低土壤水分有效性并对作物造成水分胁迫 (也常被称为盐分胁迫)，甚至对作物造成毒害。由此可见，对盐碱地作物实施灌溉需同时考 虑土壤水分与盐分的影响。灌水不足时，作物可能因根区土壤含水量过低或土壤溶液盐分 浓度过高而减缓生长速率甚至减产。然而，过量灌溉也可能会影响作物正常生长，并且因抬 高地下水位而导致土壤次生盐碱化。所以，如何通过灌溉适时调节盐碱地作物根区土壤水 分与盐分状况，使其既有利于作物生长又能减少蒸发与渗漏等损失，是实现盐碱地可持续 利用并提高水分利用效率过程中亟待解决的难题。显然，在盐碱地推广应用自动灌溉控制 系统是解决上述难题的有效途径，而灌溉控制方法则是该系统的核心所在。

[0003] 迄今为止，已有的灌溉控制方法大概可分为以下两类。第一类灌溉控制方法借助 作物对水分胁迫的生理响应（比如冠层温度、叶水势、气孔导度等生理指标的变化）来估算 作物所受到的水分胁迫程度，并据此判断灌水时间，其中基于作物冠层温度估算作物水分 胁迫指数（CWSI，表示因水分胁迫所导致的作物蒸散速率降低的程度）的方法较为常见。理 论上，由于土壤盐分也是通过影响土壤水分有效性并最终影响作物蒸腾以及冠层温度等生 理指标，所以，可直接将该方法应用于盐碱地。然而，作物冠层温度在光照期间很容易随时 间动态变化，并且受气候环境以及观测视野内土壤与作物枝干的影响较大，所以该方法存 在稳定性与代表性较差的问题，在作物生长早期冠层较为稀疏时表现得尤为突出。其次， 绝大多数情况下单纯依靠 CWSI只能判断灌水时间，为确定灌水定额，仍需借助实测的土壤 含水量剖面，除增加前期成本投入外，还给实际应用带来不便。另外，基于作物生理响应估 算CWSI并确定灌水时间时，作物往往已经受到了一定程度的胁迫，所以很难保证作物一直 在最优土壤水、盐条件下生长。因此，绝大多数方法(第二类）都根据土壤含水量进行灌溉 控制，并基于充分灌溉原理，将根区土壤含水量算术平均值(简称为根区算术平均土壤含水 量）或根区某深度处的土壤含水量(或土壤水基质势）作为灌水控制指标，当其低于最适宜 作物生长的土壤含水量下限时开始灌溉，直到其达到最适宜作物生长的土壤含水量上限。 然而，已有的第二类方法在判断灌水时间时都只考虑了土壤水分对作物生长与蒸腾的影响 而忽略了土壤盐分的影响。因此，在应用第二类方法对盐碱地作物进行灌溉控制之前，仍需 再对其进行改进以便考虑土壤水分与盐分对作物的共同胁迫，否则将给灌溉控制精度带来 较大偏差。 发明内容

[0004] (一）要解决的技术问题

[0005] 本发明要解决的技术问题为：如何根据根区土壤水分、盐分以及根系分布快速、准 确地估算CWSI进而控制灌溉。

[0006] (二）技术方案

[0007] 为了解决该技术问题，根据本发明的一方面，提出了一种用于盐碱地的灌溉控制 方法，其特征在于，该方法包括 ：

[0008] 将灌溉区域的土壤从表层至最大扎根深度L处依次分为若干层，层数记为k，用每 层土壤的平均深度 Zi除以b得到每层土壤的相对深度，用每层土壤的厚度比除以Lp 得到每层土壤的相对厚度，设定作物在每层土壤中的相对根长密度Lmd(zJ、萎焉系 数0W、饱和土壤含水量0 S、最适宜作物生长的土壤含水量上限ΘΗ与下限V、土壤水渗透 势上限fcf与下限f/.;

[0009] 测量每层土壤的含水量Θ i、温度以及土体电导率；[〇〇1〇] 计算根区加权平均土壤含水量# :
*且当θΗ < 彡03或0"时，令等于θ";当彡011时，令等于根据 每层土壤的含水量θ i、温度和土体电导率获得各土层的土壤水渗透势基于相对根长密 度分布计算根区加权平均土壤水渗透势^ :
，且当h.Hkl时, 令ft等于f;当|辦,丨2φ/.|时，令f等于f I;

[0011] 计算作物水分胁迫指数CWSI : 〇»7 = 1 -U/ff + /，）（€.ψ + έ/)，其中取日去 » 成= 时令=禁干Q 去 1 白· < 0<乂 ^ 焉 W，'ff 寸丁 W，白 • ，且当|_

[0012] 当计算得到的CWSI大于预定的CWSI临界值时，开始灌水。

[0013] 优选地，所述测量每层土壤的含水量Θ i、温度和土体电导率，是将土壤水分-温 度-电导率探针垂直插入根区土壤，并在土壤水分-温度-电导率探针上对应的各土层深 度处设有一个土壤水分-温度-电导率传感器，从而测得每层土壤的含水量Θ i、温度和土 体电导率。

[0014] 优选地，所述测得每层土壤的含水量Θ i、温度和土体电导率，是将分布在灌溉区 域内的多根土壤水分-温度-电导率探针上相同土层深度处的土壤水分-温度-电导率传 感器的测量值取平均而获得。

[0015] 优选地，该方法还包括：

[0016] 在所述最大扎根深度范围内设定计划湿润层深度，记录从表层到计划湿润层深度 的土壤层数n，设定田间持水量e f、土壤灌溉湿润比P、田间水分有效利用系数η，淋洗系 数R和灌溉面积Α，按下式计算单位面积的灌水定额Μ然后 用Μ乘以Α，得到理论灌水量；

[0017] 当实际灌水量达到理论灌水量时，停止灌水。

[0018] 优选地，该方法还包括：[〇〇19] 当计算得到的CWSI大于预定的临界值时，若预报在设定的时间区间内会有降雨， 则不灌水，否则开始灌水。

[0020] 根据本发明的另一方面，提供了一种基于作物根区土壤水分与根系分布的灌溉控 制系统，其特征在于，该系统包括参数输入子系统、土壤水盐动态监测子系统、中央决策子 系统和灌水控制子系统：

[0021] 参数输入子系统，用于设定作物的最大扎根深度L，灌溉区域内从土壤表层至最 大扎根深度L处所分的层数k，每层土壤的平均深度 Zi以及厚度hy作物在每层土壤中的 相对根长密度Lmd(zJ、萎焉系数θ"、饱和土壤含水量0 S、最适宜作物生长的土壤含水量 上限θ H与下限土壤水渗透势上限 fir与下限灼，以及CWSI临界值；

[0022] 土壤水盐动态监测子系统，用于测量每层土壤的含水量Θ i、温度以及土体电导 率，并发送给中央决策子系统；

[0023] 中央决策子系统，包括计算模块、天气预报模块和决策模块；

[0024] 计算模块，用Zi除以L，得到每层土壤的相对深度;用比除以Lp得到每层土 壤的相对厚度ΛζΗ;然后计算根区加权平均土壤含水量^ :且当ΘΗ< 03或0"时，令等于θ"，当θ,< 0 [1时，令等于 根据每层土壤的含水量、温度和土体电导率获得各土层的土壤水渗透势妁，基于 相对根长密度分布计算根区加权平均土壤水渗透势^
，且当 时，令λ等于衡；当ΙλΙ>|λΙ时，令•等于计算作物水分胁迫指数CWSI: ClfSl = 1 - {a〇 + b) {εφ + d)
二 ΘΗ < Θ 孓 θ& 或时，令；等于％，当@ <5轉时，令5等于匕；

[0025] 决策模块，比较计算得到的CWSI是否大于预定的CWSI临界值，当大于时，则向灌 水控制子系统发出开始灌水的指令。

[0026] 灌水控制子系统，包括灌溉控制模块，当收到开始灌水的指令时，控制阀门开启而 开始灌水。

[0027] 优选地，所述土壤水盐动态监测子系统，包括土壤水分-温度-电导率探针，土壤 水分-温度-电导率探针上对应每层土壤设有一个土壤水分-温度-电导率传感器，通过 将土壤水分-温度-电导率探针垂直插入根区土壤，从而测得每层土壤的含水量Θ i、温度 和土体电导率。

[0028] 优选地，在灌溉控制区域内，所述土壤水分-温度-电导率探针为多根，所述测得 每层土壤的含水量Θ i，是将灌溉区域内的多根土壤水分-温度-电导率探针上相同土层深 度处的土壤水分-温度-电导率传感器的测量值取平均而获得。

[0029] 优选地，参数输入子系统，还用于输入从表层到计划湿润层深度的土壤层数n，并 设定田间持水量9 f、土壤灌溉湿润比p、田间水分有效利用系数η，淋洗系数R和灌溉面积 A ；

[0030] 中央决策子系统的计算模块还计算单位面积的灌水定额Μ :然后用Μ乘以Α，得到理论灌水量；

[0031] 中央决策子系统的决策模块，还包括比较实际灌水量与理论灌水量，当实际灌水 量达到理论灌水量时，则向灌水控制子系统发出停止灌水的指令；

[0032] 灌水控制子系统，还包括灌水计量模块，用于计量实际灌水量并发送给中央决策 子系统；灌溉控制模块还用于在收到停止灌水的指令时，关闭控制阀门而停止灌水。

[0033] 优选地，中央决策子系统还包括天气预报模块，用于接收天气预报，当预报在设定 的时间区间内有降雨时，则决策模块不发出开始灌水的指令。

[0034] (三）有益效果

[0035] 相对于其它第二类自动灌溉控制方法所采用的根区算术平均土壤含水量，本发明 综合考虑了土壤水分、盐分以及根系分布对作物蒸腾与生长的影响，能更加准确、便捷地估 算盐碱地作物所受到的水分胁迫程度，可为盐碱地农业节水、增产目标的实现提供有效工 具。 附图说明[〇〇36] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以 根据这些附图获得其他的附图。

[0037] 图1是表示土壤水分胁迫修正系数与土壤含水量各临界值之间的关系的函数图。

[0038] 图2是表示土壤盐分胁迫修正系数与土壤水渗透势各临界值之间关系的函数图。

[0039] 图3是根据本发明一个较佳实施例的用于盐碱地的灌溉控制方法与系统的流程 图。

[0040] 图4是土壤水盐动态监测子系统一个优选实施例的结构示意图。[〇〇41] 图5是灌水控制子系统一个优选实施例的结构示意图。

[0042] 图6是各情形下基于土层相对深度（zj的土壤含水量（Θ )分布图及土壤水渗透 势（#)分布图。 具体实施方式[〇〇43] 以下将结合附图及实施例来详细说明本发明的实施方式，借此对本发明如何应用 技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据以实施。需要说明 的是，只要不构成冲突，本发明中的各个实施例以及各实施例中的各个特征可以相互结合， 所形成的技术方案均在本发明的保护范围之内。

[0044] 首先说明本发明所采用的计算盐碱地作物水分胁迫指数CWSI的原理。

[0045] 由定义可知，作物水分胁迫指数（CWSI)可用水分胁迫所导致的作物蒸散速率降低 的程度来表征 ：Γ1 Ί

[0046] LU

[0047] 式中ETa为作物实际蒸散速率，cm cf1 ;ETp为作物潜在蒸散速率，cm cf1 ;Ta为作物 实际蒸腾速率，cm cf1 ;TP为作物潜在蒸腾速率，cm cf1屯为实际土面蒸发速率，cm cf1 ;EP 为潜在土面蒸发速率，cm cf1。一般情况下，相对于旱地作物蒸腾强度而言，农田土面蒸发强 度都比较小，当作物叶片生长茂盛时更是如此，因此常被忽略。于是，式[1]可简化为：

[0048][2]
[〇〇49] 除土壤含水量可直接对作物造成水分胁迫外，土壤盐分通过降低土壤水分有效性 也可影响作物对水分的吸收(蒸腾)，从而使作物遭受一定程度的水分胁迫(有时也被称为 盐分胁迫)。当考虑土壤水分与盐分对作物蒸腾的影响时，根系吸水速率通常可表示为：

[0050][3]

[0051] 式中z为垂直坐标，取地表为原点，向下为正，cm ;S(z)为根系吸水速率， cn^cm^f1 ; Θ为土壤体积含水量，cm3cm3 ; γ ( Θ )为基于土壤含水量计算的土壤水分胁 迫修正系数；f为土壤水渗透势，cm ; /i(f)_为基于土壤水渗透势计算的土壤盐分胁迫 修正系数；Smax(z)为最大根系吸水速率，表示在最优土壤水分条件下的根系吸水速率， cn^cm^f1 A为最大扎根深度，cm & (=ζ/Ι^)为土层相对深度；Lmd(zJ为相对根长密度，其中Ld(zJ为根长密度，cm cnT3。当忽略作物体内水分含量变化时， 作物实际蒸腾速率可按下式估算：

[0052][4]

[0053] 将式[3]代入式[4]，可得：

[0054][5]

[0055] 将式[5]代入式[2]，可得：

[0056][6]

[0057] 1)简化方法：基于根区算术平均土壤含水量与土壤水渗透势估算CWSI

[0058] 假设根区范围内土壤含水量与土壤水渗透势都均匀一致，即各深度处土壤含水量 与土壤水渗透势都按根区算数平均值（^，cm3CnT3;i, cm)来计，则由式[6]可得：

[0059] m

[0060] 2)改进方法：基于根区土壤水分、盐分以及作物根系分布估算CWSI

[0061] 大量相关研究结果表明，土壤水分胁迫修正系数Y ( Θ )与土壤盐分胁迫修正系 数可分别表示为土壤含水量与土壤水渗透势的线性函数，如式[8]与式[9]所示：

[0062][8]

[0063] 式中a与b分别为系数；Θ s为饱和土壤含水量，cm3cnf3 ; Θ w为萎焉系数，cm3cnf3 ; ΘΗ与V分别为最适宜作物生长或根系吸水的土壤含水量上限与下限，cm3cm 3。式[8]表 明（如图1所示):当ΘΗ〈Θ < 0S时，由于土壤含水量太高，土壤通气性太差，作物根系无法 吸水，Y ( Θ )=〇 ;当Θ彡Θ w时，由于土壤含水量太低,作物根系也无法吸水，Y ( Θ )=〇 ;当 θ'θ彡011时，最适宜作物根系吸水，Y(e)=i;i θ¥〈θ彡θ,时，作物根系吸水速率随 着土壤含水量的降低而线性递减，Y ( Θ )由1递减为0。

[0064] 聞

[0065] 式中c为一个基于土壤水渗透势的削减系数，cnT1 ;d为无量纲系数；fe与ft分别 为土壤水渗透势上限和下限，cm。式[9]表明（如图2所示）：当|炉| S |%|时，表明根系吸 水不受土壤盐分影响；当时，表明作物根系吸水速率随着土壤水渗透势的降 低而线性递减；当|炉|^妗|时，表明土壤盐分含量过高，作物已无法吸收水分。

[0066] 将式[8]与式[9]分别代入式[6]可得：

[0067][10]

[0068] 由式[8]与式[9]可知，式[10]中的系数a、b、c、d在根区内可能不是常数，会随 各土层深度处土壤含水量或土壤水渗透势发生改变。为了让系数a、b、c、d在根区范围内 为常数(便于计算)，在保证不改变原有方程结果的基础之上，本发明将式[8]与式[9]分别 改写为：

[0069][11]

[0070][12]

[0071] 将式[11]与式[12]分别代入式[10]可得：

[0072][13]

[0073] 其中：

[0074][14]

[0075] 式中^为基于相对根长密度分布计算的根区土壤含水量加权平均值，简称为根区 加权平均土壤含水量，cm3cnT3 ; $为基于相对根长密度分布计算的根区土壤水渗透势加权 平均值，简称为根区加权平均土壤水渗透势，cm ;i为根区土层编号，从表层至最大扎根深 度处依次为第1层、第2层、第3层……第k层；Θ i为第i层土壤的含水量，cm3cm3 ;Zii为 第i层土壤的相对深度。由此可见，当土壤水分胁迫修正系数与盐分胁迫修正系统分别用 土壤含水量与土壤水渗透势的线性函数表达时，CWSI可表示为|与^的函数。由式[13]可 知：当饵^丨=丨即土壤盐分对作物生长或蒸腾没有影响时，作物水分胁迫程度主要由根区 土壤水分与根系分布决定，〇F57 = l-当即土壤含水量本身对作物蒸腾没有 影响时，作物水分胁迫程度主要由根区土壤盐分与根系分布决定，CWS/ = 1 β

[0076] 上面介绍了本发明盐碱地作物CWSI计算公式的推导过程，下面通过一个较佳实 施例介绍本发明的用于盐碱地的灌溉控制方法与系统。图3是本发明一个较佳实施例的用 于盐碱地的灌溉控制方法与系统的流程图：

[0077] 在步骤S1，设置参数，具体可以包括：将灌溉区域从土壤表层至最大扎根深度L 处依次分为若干层，层数记为k，用每层土壤的深度Zi除以L，得到每层土壤的相对深度 z#用每层土壤的厚度匕除以L，得到每层土壤的相对厚度Λ ，设定作物在每层土壤中 的相对根长密度Lmd(zJ、萎焉系数θ"、饱和土壤含水量0S、最适宜作物生长的土壤含水 量上限Θ Η与下限L、土壤水渗透势上限的ί与下限你；

[0078] 在步骤S2,测量每层土壤的含水量Θ i、温度以及土体电导率。其中一种优选的测 量方式是：将土壤水分-温度-电导率探针垂直插入根区土壤，土壤水分-温度-电导率探 针上对应每层土壤设有一组土壤水分-温度-电导率传感器，从而测得每层土壤的含水量 Θ i、温度以及土体电导率；并且进一步将灌溉区域内的多根土壤水分-温度-电导率探针 上相同土层深度处的土壤水分-温度-电导率传感器的测量值取平均，获得平均值。

[0079] 在步骤S3,基于相对根长密度分布计算根区加权平均土壤含水量；： ，且当ΘΗ< 03或0"时，令等于0W，当 i «iθ^< 0i< ΘΗ时，令0i等于V;根据每层土壤的含水量0i、温度和土体电导率获得 各土层的土壤水渗透势基于相对根长密度分布计算根区加权平均土壤水渗透势^ :
>且当hi < hi时，令約等于御/;当μμ μΐ时，令物等于 φΰ

[0080] 在步骤S4,计算作物水分胁迫指数CWSI : CfSf = 1 - iai + AMe; + 〇0 ,其 中取，且当·< 爲时，令i等于％，当 A 时'令i等于0L，取
'且当时，令炉 等于ph;当b I时，令5等于史厂，

[0081] 在步骤S5 :当计算得到的CWSI大于预定的CWSI临界值，开始灌水。在作物生育 期内，可以根据作物各生长发育阶段以及具体的灌溉要求(充分灌溉或非充分灌溉）设定启 动灌水的CWSI临界值。优选地，还需要根据天气预报得到的未来天气决定是否开始灌溉， 即当预报在设定的时间区间内有降雨时则不开始灌水，以避免水资源浪费，设定的时间区 间例如为一至三天。

[0082] 进一步，还可以确定灌水定额，并计算理论灌水量，从而判断灌水终止时间。

[0083] 根据下式确定灌水定额：

[0084][15]

[0085] 式中Μ为灌水定额，m3亩4 ;n为从表层到计划湿润层深度的土壤层数，计划湿润 层深度在所述最大扎根深度范围内选定为第i层土壤的土层厚度，cm ; Θ i为灌水前测 得的第i层土壤的土壤含水量，cm3cnT3 ; Θ f为田间持水量，对于砂土、壤土、粘土可分别取为 土壤水基质势_l〇〇cm、-200cm、-300cm所对应的土壤含水量，cm 3cnT3 ;p为土壤灌溉湿润比； R为淋洗系数；Π 为田间水分有效利用系数。土壤从表层到最大扎根深度分为k层，因为作 物根系主要汇集在浅层土壤中，为了防止深层渗漏，一般取计划湿润层深度小于最大扎根 深度，即n彡k。

[0086] 理论灌水量计算公式如下：

[0087] Qa= ΜΑ [16]

[0088] 式中QaS理论灌水量，m3 ;Α为控制灌溉面积，亩。当实际灌水量达到理论灌水量 Qa时，则停止灌水。[〇〇89] 下面介绍本发明一个实施例的自动灌溉控制系统，该系统包括参数输入子系统、 土壤水盐动态监测子系统、中央决策子系统和灌水控制子系统。

[0090] 参数输入子系统，用于设定作物的最大扎根深度Lp灌溉区域的土壤从表层至最 大扎根深度L所分的层数k，每层土壤的平均深度 Zi以及厚度hp需要灌溉的作物每层的 相对根长密度Lmd(zJ、萎焉系数θ"、饱和土壤含水量0 S、最适宜作物生长的土壤含水量 上限ΘΗ与下限V、土壤水渗透势上限彻与下限妁，以及CWSI临界值；

[0091] 为了计算理论灌水量，参数输入子系统还包括输入从表层到计划湿润层深度的土 壤层数η、田间持水量e f、土壤灌溉湿润比Ρ、田间水分有效利用系数η、淋洗系数R和灌 溉面积Α。

[0092] 上述参数，可以使用现有的经验数据，也可以输入实测数据，比如：[〇〇93] 需根据作物耗水规律输入最适宜作物生长或根系吸水的土壤含水量上限与下限， 对于不同作物甚至作物不同生育阶段，其值可能存在差异。另外，模块中需输入的萎焉系数 一般取为土壤水基质势-15000cm所对应的土壤含水量。

[0094] 可以通过以下三种方式获取最大扎根深度：1)输入实测数据；2)根据经验关系式 进行估算；3)利用作物扎根深度模型进行模拟。相对根长密度分布资料可以直接输入实 测数据，也可利用一些已公开发表的简化关系式（比如将相对根长密度表达为相对深度的1 次、2次或3次函数关系)或特定作物的统计结果进行估算，比如对于小麦，可通过下式进行 估算：L nrd(zr) = 4. 522(l-zr) 5.228exp(9· 644zr2.426)。

[0095] 对于计划湿润层深度，在作物生长苗期，一般使其等于最大扎根深度；在作物生长 中后期，可让其小于最大扎根深度，一般在0. 4-0. 6m之间取值。[〇〇96] 根据作物耗水规律与实际灌水要求输入启动灌水的CWSI临界值：对于充分灌溉， 可将CWSI临界值设为0. 001，对于非充分灌溉，则需根据每次灌水的实际要求而定。需根据 灌水方式与作物种植株行距等资料输入土壤灌溉湿润比（P)，对于全面灌溉(畦灌、沟灌、喷 灌等)，P=l，对于局部灌溉(微喷灌、滴灌等)，P〈l。需根据具体洗盐要求确定淋洗系数，一 般取R3 1.0。另外，需根据灌水方法以及一些经验资料，输入田间水分有效利用系数，灌水 方法越先进，其值就越高。

[0097] 土壤水盐动态监测子系统，用于测量最大扎根深度范围内每层土壤的含水量θ ρ 温度以及土体电导率，并发送给中央决策子系统，优选为分布在灌溉区域内的多根土壤水 分-温度-电导率探针，土壤水分-温度-电导率探针上对应每层土壤设有一个土壤水 分-温度-电导率传感器，将土壤水分-温度-电导率探针垂直插入根区土壤，从而测得每 层土壤的含水量Θ i，并将相同土层深度处的土壤水分-温度-电导率传感器的测量值取 平均。土壤水分-温度-电导率探针可按实际要求调整自动采集数据的时间间隔（时间步 长，如0. 5h或1. Oh等)。土壤水分-温度-电导率探针内埋设的土壤水分-温度-电导率 传感器数量以及垂直方向间距（空间步长）也可根据实际情况(如作物全生育期内最大扎根 深度）进行调整，一般从地表以下5cm开始，以10-20cm为空间步长向下添加土壤水分-温 度-电导率传感器。为了减少空间变异所带来的影响，可以根据具体要求设定灌溉区域内 土壤水分-温度-电导率探针的数量。

[0098] 图4示出了土壤水盐动态监测子系统一个优选实施例的结构示意图，其中11为土 壤水分-温度-电导率探针，插入土壤中；12为土壤水分-温度-电导率传感器，按照土壤 分层数和土层厚度设置在土壤水分-温度-电导率探针上；13为无线传输装置，用于将土 壤水分-温度-电导率传感器12的测量值发送给中央决策子系统。14为GPS定位装置； 15为太阳能板，用于供电；16为支架，支撑太阳能板14。[〇〇99]中央决策子系统，包括计算模块和决策模块；计算模块，用zi除以L，得到每层 土壤的相对深度^ ;计算根区加权平均土壤含水量# :
且当 ΘΗ< 彡03或彡0"时，令等于θ"，当θ,< 彡0[1时，令等于 根据每层土壤的含水量Θ i、温度和土体电导率获得各土层的土壤水渗透势釣，计算根区 加权平均土壤水渗透势;
且当丨时，令狗等于彻;当 |供.| I时，令的等于* ;计算作物水分胁迫指数CWSI : 0157 = 1 - + /?) (αφ + d) s 其中®
，且当& < 久或is焉时，令；等于％，当 6 心时，令g等于取
，且当时，令; 等于当旷》h/.i时，令$等于攸；

[0100] 决策模块，比较计算得到的CWSI是否大于预定的CWSI临界值时，当大于时，则向 灌水控制子系统发出开始灌水的指令。

[0101] 为了计算灌水定额和理论灌水量，中央决策子系统的计算模块还计算从表层到计 划湿润层深度的土壤层数n，计算单位面积的灌水定额Μ :然后用Μ乘以Α，得到理论灌水量Qa ;此时，中央决策子系统的决策模块向灌水控制子系统 发出的开始灌水的指令包括理论灌水量Qa。

[0102] 中央决策子系统还可包括天气预报模块，用于接收天气预报，当预报设定的时间 区间内比如一至三天内有降雨时，决策模块不发出开始灌水的指令。以避免灌水后短时间 内又出现较大降雨的情况，从而提高水分利用效率。

[0103] 中央决策子系统例如可由计算机构成，并包括用于接收和传送数据的无线传输装 置。对于小型灌区，可将中央决策子系统中所需的相关软件或程序直接集成到灌水控制子 系统中，从而将其置于灌溉控制单元田块内部，无需再配置专用计算机。

[0104] 灌水控制子系统，包括灌溉控制模块，当收到开始灌水的指令时，则控制阀门开启 而灌水。灌水控制子系统还可以包括灌水计量模块，用于计量实际灌水量。灌溉控制模块 还用于在实际灌水量达到理论灌水量时，控制阀门关闭而停止灌水。

[0105] 图5示出了灌水控制子系统一个优选实施例的结构示意图，灌水控制子系统可位 于其控制的灌概控制单元26 (即控制的灌溉区域中的灌溉管道装置)中，包括电磁阀21、智 能水表22、无线传输装置23、GPS定位装置24和太阳能供电装置25,在无线传输装置23接 到中央决策子系统所发送的灌水指令后，控制电磁阀21开启实施灌溉，智能水表22记录动 态灌水量，并将数据通过无线传输装置23传送给中央决策子系统。当实际灌水量达到理论 灌水量时，无线传输装置23接到中央决策子系统所发送的终止灌水指令，控制电磁阀21关 闭，停止灌溉。在这里也可以包括一个与无线传输装置23相连的例如为单片机形式的灌溉 控制器件（图中未示出），接收包括理论灌水量的灌水指令并实时获得动态灌水量，从而控 制开始灌水和停止灌水，此时就不需要将实际灌水量发送到中央决策子系统。GPS定位装置 24主要用于确定、跟踪灌水控制子系统的具体位置，便于管理。电磁阀21、智能水表22以 及无线传输装置23和GPS定位装置24所需电能均由太阳能供电装置25提供。

[0106] 为了验证本发明关键技术，特设置了一个数值试验：保持根区算术平均土壤含水 量与土壤水渗透势、相对根长密度分布等条件不变，对比分析各土壤水分与与盐分分布方 式下两种方法所估算的CWSI及其与理论值之间的相对误差。本发明适合于盐碱地旱作物 的灌溉。

[0107] 在该数值试验中，设土壤为粉质壤土，Θ s=〇. 450cm3cnT3 ; Θ ^0. 〇67cm3cnT3 ;van Genuchten (1980)土壤水分特征曲线中 a =〇· 〇2，n=L 41 ; Θ f=〇. 247cm3cnT3。设最适宜作物 生长的土壤含水量上限取为土壤水基质势-50cm所对应的土壤含水量，Θ H=〇. 380cm3Cm_3 ; 最适宜作物生长的土壤含水量下限取为田持的80%，即Θ f〇. 197cm3CnT3 ;萎焉系数取为土 壤水基质势-15000cm所对应的土壤含水量，Θ w=〇. l〇4cm3cnT3 ;影响根系吸水的土壤水渗 透势临界值-720 cm，歡=-5650 cm;最大扎根深度L=4〇cm ;相对根长密度分布为 相对土层深度的函数Lmd(zJ = 4.522(l-%)5_228eXp(9.644z/_426);作物潜在蒸腾速率为 0.60cm cf1;潜在土面蒸发速率为0.03cm cf1;实际土面蒸发速率根据表层土壤含水量确 定，当表层土壤含水量大于田持时，实际土面蒸发速率等于潜在土面蒸发速率，当表层土壤 含水量介于田持与残余含水量之间时，实际土面蒸发速率由0.03cm cT1线性递减至零。以 垂直方向上土壤水分线性分布以及保持根区算术平均土壤含水量不变为原则，设置3种土 壤水分分布方式：根区土壤含水量均勻分布并且j = 0.182 cm3 cm-3,空间步长为icm (W1);以W1为参考，地表处土壤有效含水量（Θ f- Θ J分别增大（W2)或缩小（W3) 0. 6倍，最 大扎根深度处土壤有效含水量则相应地缩小或增大0. 6倍(见图6左侧)。以垂直方向上土 壤水渗透势线性分布以及保持根区算术平均土壤水渗透势不变为原则，设置3种土壤盐分 分布方式：根区土壤水渗透势均匀分布并且j =-3000 cm,空间步长为lcm (si);以si 为参考，地表处土壤水渗透势分别增大（S2)或缩小（W3) 0. 6倍，最大扎根深度处土壤有效 含水量则相应地缩小或增大〇. 6倍(见图6右侧)。本数值试验共讨论五种土壤水分与盐分 分布情形：W1S1、W2S2、W2S3、W3S2、W3S3。

[0108] 在各种土壤水分、盐分分布条件下，先按式[4]计算作物实际蒸腾速率，然后再按 式[1]计算CWSI理论值；按简化方法(式[7])与改进方法(式[13])分别估算CWSI。各情 形下CWSI理论值与估算值以及其对比关系如表1所示。表1表明：①在根区算术平均土 壤含水量与土壤水渗透势保持不变的情况下，作物水分状况因土壤水分与盐分分布方式存 在很大差异，土壤水分与根系分布(表层多，深层少）越一致、土壤盐分与根系分布越背离， CWSI越小，反之则CWSI越大；②相对于简化方法，采用改进方法更能准确估算CWSI，相对误 差小于10%。

[0109] 采用本发明的自动灌溉方法和自动灌溉系统，基于作物根长密度分布计算根区加 权平均土壤含水量与根区加权平均土壤水渗透势，从而更加准确地估算盐碱地作物所受到 的水分胁迫程度，以此为依据控制灌溉时既考虑了根区土壤水分与盐分状况又考虑了作物 水分状况，可为盐碱地农业节水、增产目标的实现提供有效工具。本发明仅适合对旱作物进 行灌溉控制。

[0110] 以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例 对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施 例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者 替换，并不使相应技术方案的本质脱离本发明权利要求所限定的范围。

[0111] 表1各土壤含水量与土壤水渗透势分布情形下CWSI的理论值与不同方法估算值 之间的对比关系

[0112]