CN111640038A - 水稻作物系数计算方法以及水稻灌溉系统 - Google Patents
水稻作物系数计算方法以及水稻灌溉系统 Download PDFInfo
- Publication number
- CN111640038A CN111640038A CN202010446542.7A CN202010446542A CN111640038A CN 111640038 A CN111640038 A CN 111640038A CN 202010446542 A CN202010446542 A CN 202010446542A CN 111640038 A CN111640038 A CN 111640038A
- Authority
- CN
- China
- Prior art keywords
- crop
- growth
- rice
- module
- coefficient
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 241000209094 Oryza Species 0.000 title claims abstract description 127
- 235000007164 Oryza sativa Nutrition 0.000 title claims abstract description 126
- 235000009566 rice Nutrition 0.000 title claims abstract description 126
- 238000004364 calculation method Methods 0.000 title claims abstract description 65
- 238000003973 irrigation Methods 0.000 title claims abstract description 51
- 230000002262 irrigation Effects 0.000 title claims abstract description 51
- 230000012010 growth Effects 0.000 claims abstract description 171
- 238000012360 testing method Methods 0.000 claims abstract description 36
- 238000005457 optimization Methods 0.000 claims abstract description 32
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 35
- 238000004891 communication Methods 0.000 claims description 22
- 239000003621 irrigation water Substances 0.000 claims description 20
- 230000006870 function Effects 0.000 claims description 17
- 238000009795 derivation Methods 0.000 claims description 16
- 238000000034 method Methods 0.000 claims description 16
- 239000002245 particle Substances 0.000 claims description 14
- 238000002474 experimental method Methods 0.000 claims description 8
- 230000004907 flux Effects 0.000 claims description 6
- 230000005855 radiation Effects 0.000 claims description 6
- 239000002689 soil Substances 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 238000004383 yellowing Methods 0.000 claims description 4
- 238000003306 harvesting Methods 0.000 claims description 2
- 229920006395 saturated elastomer Polymers 0.000 claims description 2
- 238000009331 sowing Methods 0.000 claims description 2
- 230000005068 transpiration Effects 0.000 abstract description 6
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 239000000047 product Substances 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 240000002834 Paulownia tomentosa Species 0.000 description 1
- 235000010678 Paulownia tomentosa Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000003698 anagen phase Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Agronomy & Crop Science (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Husbandry (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了水稻作物系数计算方法以及水稻灌溉系统,能够精确计算得到水稻作物系数,有利于科学合理地实施水稻灌溉。水稻作物系数计算方法包括:收集灌溉试验站多年逐日的气象资料以及试验推导的水稻作物系数;根据水稻生长情况,将水稻生长期划分为生长初期、快速生长期、生长中期、生长末期;根据气象资料计算参考作物腾发量,通过作物系数kc与参考作物腾发量乘积得到实际作物腾发量;修订不同生长期作物系数基础值;进行优化求解,得到生长初期作物系数基础修订值、生长中期作物系数基础修订值和生长后期末作物系数基础修订值;统计水稻生长中期、生长末期的日平均风速、日最低相对湿度;计算得到水稻不同生长阶段的作物系数。
Description
技术领域
本发明属于节水灌溉领域,具体涉及水稻作物系数计算方法以及水稻灌溉系统。
背景技术
灌溉用水定额是农业灌溉节水、水资源优化配置、水资源合理利用的重要依据,作物需水量是灌溉用水定额计算中重要水平衡项,而作物系数是作物需水量计算和预测中重要参数,直接影响其结果的准确性和可靠性。水稻作为灌溉需水量较大的作物,确定水稻作物系数对合理实施灌溉具有重要意义。水稻作物系数推求一般采用试验数据反推法和单作物系数法,其中,试验数据反推法计算方法简单,但需要多年的水稻需水量试验实测数据,而单作物系数法在推求水稻作物系数时,由于受区域、作物品种、气候等因素影响,该方法在实际使用时误差较大,并不适用,因此,如何精确计算推求水稻作物系数成为了亟待解决的问题。
发明内容
本发明是为了解决上述问题而进行的,目的在于提供一种水稻作物系数计算方法以及水稻灌溉系统,能够精确计算得到水稻作物系数,有利于科学合理地实施水稻灌溉。
本发明为了实现上述目的,采用以下方案:
<方法>
本发明提供一种水稻作物系数计算方法,其特征在于,包括如下步骤:
步骤1.收集灌溉试验站多年逐日的日照时数、平均相对湿度、平均气温、最高气温、最低气温、平均风速气象资料以及试验(水稻需水量试验)推导的水稻作物系数kc;
步骤2.根据水稻生长情况,将水稻生长期划分为生长初期、快速生长期、生长中期、生长末期;
步骤3.根据气象资料计算参考作物腾发量ET0,通过作物系数kc与参考作物腾发量乘积得到实际作物腾发量ETc:
式中,Rn为太阳净辐射,MJ/m2·d,可由经验公式计算;G为土壤热通量,MJ/m2·d;γ为湿度表常数,kPa·℃-1;u2为2m处风速,m/s;T为平均气温,℃;ed为实际水汽压,kPa;ea为饱和水汽压,kPa;Δ为温度~饱和水汽压关系曲线T处的切线斜率,kPa·℃-1;
步骤4.修订不同生长期作物系数基础值
假定作物系数试验推导值推求的作物腾发量为实际作物腾发量,目标函数为作物系数基础修订值推求的作物腾发量与作物系数试验推导值推求的作物腾发量相对误差最小;
目标函数F:
约束条件:
ETc,i=Kc,i×ET0,i,
ET′c,i=K′c,i×ET0,i,
a≤K′cini≤b,K′cini≤K′cmid≤c,d≤K′cend≤K′cmid,
式中:ETc,i为作物系数试验推导值推求的日序数为i的作物腾发量,mm;Kc,i为日序数为i的作物系数试验推导值;ET0,i为日序数为i的参考作物腾发量,mm;K′c,i为基于基础修订值的日序数i作物系数理论推导值;K′cini为生长初期作物系数基础修订值,自变量;K'cmid为生长中期作物系数基础修订值,自变量;K'cend为生长后期末作物系数基础修订值,自变量;n0为生长初期起始日序数;n1为生长初期终止日序数;n2为快速生长期终止日序数;n3为生长中期终止日序数;n4为生长后期终止日序数;a和b为生长初期作物系数基础修订值下限和上限;c为生长中期作物系数基础修订值的上限;d为生长后期末作物系数基础修订值的下限;
步骤5.采用粒子群算法进行优化求解,得到生长初期作物系数基础修订值、生长中期作物系数基础修订值和生长后期末作物系数基础修订值;
步骤6.统计水稻生长中期、生长末期的日平均风速u′2、日最低相对湿度RH′min;
步骤7.采用以下公式计算得到水稻不同生长阶段的作物系数:
其中,水稻生长中期、生长后期末水稻作物系数的修订公式为:
优选地,本发明提供的水稻作物系数计算方法,还可以具有如下特征:在步骤2中,生长初期为从水稻播种到作物覆盖率接近10%,快速生长期为从覆盖率10%到大田作物充分覆盖,覆盖率达到70%~80%,生长中期为从充分覆盖到成熟期开始到叶片开始变黄,生长后期为从叶片开始变黄到生理成熟或收获。
优选地,本发明提供的水稻作物系数计算方法,还可以具有如下特征:在步骤5中是采用粒子群算法进行优化求解,
vmin,d,vmax,d分别为最小和最大速度,
其中,1≤d≤D,1≤j≤M,(D=3,M为最大迭代次数)
则粒子在t+1时刻的位置可通过下式更新:
式中,r1,r2为均匀分布在(0,1)区间的随机数;ω为惯性权重,取0.92;c1,c2为学习因子,通常取c1=c2=2。
<系统>
进一步,本发明还提供一种水稻灌溉系统,其特征在于,包括:灌溉装置;和控制装置,基于上文<方法>中所描述的水稻作物系数计算方法得到水稻灌溉用水量,并根据水稻灌溉用水量控制灌溉装置对水稻进行灌溉。
优选地,本发明提供的水稻灌溉系统,其特征在于,控制装置包括:
参数获取模块,获取灌溉试验站多年逐日的日照时数、平均相对湿度、平均气温、最高气温、最低气温、平均风速气象资料以及试验推导的水稻作物系数kc;
生长期划分模块,与参数获取模块通信相连,根据水稻生长情况,将水稻生长期划分为生长初期、快速生长期、生长中期、生长末期;
实际作物腾发量计算模块,与参数获取模块相连,根据气象资料计算参考作物腾发量ET0,通过作物系数kc与参考作物腾发量乘积得到实际作物腾发量ETc:
式中,Rn为太阳净辐射;G为土壤热通量;γ为湿度表常数;u2为2m处风速;T为平均气温;ed为实际水汽压;ea为饱和水汽压;Δ为温度~饱和水汽压关系曲线T处的切线斜率;
作物系数基础值修订模块,假定作物系数试验推导值推求的作物腾发量为实际作物腾发量,目标函数为作物系数基础修订值推求的作物腾发量与作物系数试验推导值推求的作物腾发量相对误差最小;
目标函数F:
约束条件:
ETc,i=Kc,i×ET0,i,
ET′c,i=K′c,i×ET0,i,
a≤K′cini≤b,K′cini≤K′cmid≤c,d≤K′cend≤K′cmid,
式中:ETc,i为作物系数试验推导值推求的日序数为i的作物腾发量;Kc,i为日序数为i的作物系数试验推导值;ET0,i为日序数为i的参考作物腾发量;K′c,i为基于基础修订值的日序数i作物系数理论推导值;K′cini为生长初期作物系数基础修订值,自变量;K'cmid为生长中期作物系数基础修订值,自变量;K'cend为生长后期末作物系数基础修订值,自变量;n0为生长初期起始日序数;n1为生长初期终止日序数;n2为快速生长期终止日序数;n3为生长中期终止日序数;n4为生长后期终止日序数;a和b为生长初期作物系数基础修订值下限和上限;c为生长中期作物系数基础修订值的上限;d为生长后期末作物系数基础修订值的下限;
优化求解模块,与参数获取模块、生长期划分模块、实际作物腾发量计算模块均通信相连,基于作物系数基础值修订模块中的目标函数F和约束条件进行优化求解,得到生长初期作物系数基础修订值、生长中期作物系数基础修订值和生长后期末作物系数基础修订值;
统计模块,与参数获取模块、生长期划分模块均通信相连,统计水稻生长中期、生长末期的日平均风速u′2、日最低相对湿度RH′min;
作物系数计算模块,与参数获取模块、生长期划分模块、实际作物腾发量计算模块、优化求解模块、统计模块均通信相连,根据以下公式计算出水稻不同生长阶段的作物系数:
其中,水稻生长中期、生长后期末水稻作物系数的修订公式为:
灌溉用水量计算模块,与作物系数计算模块通信相连,基于水稻不同生长阶段的作物系数计算得到水稻灌溉用水量;
控制模块,与参数获取模块、生长期划分模块、实际作物腾发量计算模块、作物系数基础值修订模块、优化求解模块、统计模块、作物系数计算模块、灌溉用水量计算模块、灌溉装置均通信相连,控制它们的运行。
优选地,本发明提供的水稻灌溉系统,其特征在于,优化求解模块采用粒子群算法进行优化求解:
vmin,d,vmax,d分别为最小和最大速度,
其中,1≤d≤D,1≤j≤M,(D=3,M为最大迭代次数)
则粒子在t+1时刻的位置可通过下式更新:
式中,r1,r2为均匀分布在(0,1)区间的随机数;ω为惯性权重,取0.92;c1,c2为学习因子,通常取c1=c2=2。
优选地,本发明提供的水稻灌溉系统,其特征在于,还可以包括:输入显示模块,与参数获取模块、生长期划分模块、实际作物腾发量计算模块、作物系数基础值修订模块、优化求解模块、统计模块、作物系数计算模块、灌溉用水量计算模块、控制模块、灌溉装置均通信相连,用于让用户输入操作指令,并进行相应显示。
发明的作用与效果
本发明所提供的水稻作物系数计算方法由于采用了以上步骤,因此能够科学合理地推求水稻作物系数;进一步,本发明所提供的水稻灌溉系统,采用控制装置基于上述水稻作物系数计算方法得到水稻灌溉用水量,并根据水稻灌溉用水量控制灌溉装置对水稻进行灌溉,能够高效、准确、科学、合理地实施水稻灌溉。
附图说明
图1为本发明实施例中涉及的鄂中丘陵区中稻生长期划分情况示意图;
图2为本发明实施例中涉及的迭代过程目标函数变化过程示意图;
图3为本发明实施例中涉及的不同方法确定的kc推求中稻ETC的对比图。
具体实施方式
以下结合附图对本发明涉及的水稻作物系数计算方法以及水稻灌溉系统的具体实施方案进行详细地说明。
<实施例>
本实施例所提供的水稻作物系数计算方法,包括如下步骤:
步骤1.根据湖北省灌溉用水定额分区,选择湖北省鄂中丘陵地区作为参照地区,搜集了长渠站、团林站和徐家河试验站的中稻灌溉试验资料,根据试验观测资料推求了鄂中丘陵地区中稻作物系数,具体见下表1。收集了鄂中丘陵区宜城、荆门、钟祥、京山、应城、安陆、云梦、孝感等主要气象站点的1973~2013年逐日的日照时数、平均相对湿度、平均气温、最高气温、最低气温、平均风速等气象资料。
表1鄂中丘陵区不同生育期中稻作物系数试验推导值
步骤2.如图1所示,根据鄂中丘陵区水稻生长情况,将水稻生长期划分为生长初期、快速生长期、生长中期、生长末期。
步骤3.根据鄂中丘陵区各气象站的气象资料利用彭曼公式计算各气象站1973-2013年逐日参考作物腾发量,确定参照区域1973-2013年逐日参考作物腾发量。
式中,Rn为太阳净辐射,MJ/m2·d,可由经验公式计算;G为土壤热通量,MJ/m2·d;γ为湿度表常数,kPa·℃-1;u2为2m处风速,m/s;T为平均气温,℃;ed为实际水汽压,kPa;ea为饱和水汽压,kPa;Δ为温度~饱和水汽压关系曲线T处的切线斜率,kPa·℃-1。
步骤4.修订不同生长期作物系数基础值。假定作物系数试验推导值推求的作物腾发量为实际作物腾发量,目标函数为作物系数基础修订值推求的作物腾发量与作物系数试验推导值推求的作物腾发量相对误差最小。
目标函数F:
约束条件:
ETc,i=Kc,i×ET0,i
ET′c,i=K′c,i×ET0,i
a≤K′cini≤b
K′cini≤K′cmid≤c
d≤K′cend≤K′cmid
式中:Kc,i为日序数为i的作物系数试验推导值,取表1中对应数值;ET0,i为日序数为i的参考作物腾发量,取步骤3计算值,mm;n0为生长初期起始日序数,本实施例中取153;n1为生长初期终止日序数,本实施例中取175;n2为快速生长期终止日序数,本实施例中取201;n3为生长中期终止日序数,本实施例中取230;n4为生长后期终止日序数,本实施例中取255;a和b为生长初期作物系数基础修订值下限和上限,本实施例中分别取0.7和1.4;c为生长中期作物系数基础修订值的上限,取1.7;d为生长后期末作物系数基础修订值的下限,取0.8。
步骤5.采用粒子群进行优化求解:
其中,1≤d≤D,1≤j≤M,(D=3,M为最大迭代次数。)
则粒子在t+1时刻的位置可通过下式更新:
式中,r1,r2为均匀分布在(0,1)区间的随机数;ω为惯性权重,取0.92;c1,c2为学习因子,通常取c1=c2=2;粒子数取600个,迭代次数为400次。目标函数指标变化过程见图2。中稻作物系数基础修订值见下表2。
表2鄂中丘陵区中稻作物系数基础修订值
参照地区 | K<sub>cini</sub> | K<sub>cmid</sub> | K<sub>cend</sub> |
鄂中丘陵区 | 1.11 | 1.42 | 0.80 |
步骤6.统计水稻生长中期、生长末期的日平均风速u′2、日最低相对湿度RH′min。
步骤7.采用以下公式计算得到水稻不同生长阶段的作物系数:
另外,为验证该计算方法精度,推求荆门市中稻作物系数,重复步骤6,并带入步骤7的计算公式,得到中稻作物系数修正值。同时,直接采用FAO推求荆门市中稻作物系数。具体见下表3。采用作物系数试验值kc、FAO推导值、作物系数修正值分别推导荆门市中稻作物需水量,具体见图3。直接采用FAO推导值推求作物需水量误差为-13.1%,而采用本发明提出的方法推求作物需水量误差仅为0.7%,极大提高了计算精度。
不同方法推求KC | K<sub>cini</sub> | K<sub>cmid</sub> | K<sub>cend</sub> |
作物系数修正值 | 1.11 | 1.42 | 0.80 |
FAO推导值 | 1.05 | 1.14 | 0.85 |
进一步,本实施例还涉及一种用于自动化实施水稻灌溉的水稻灌溉系统,该系统包括灌溉装置和控制装置。
控制装置基于上文所描述的水稻作物系数计算方法得到水稻灌溉用水量,并根据水稻灌溉用水量控制灌溉装置对水稻进行灌溉。具体地,控制装置包括参数获取模块、生长期划分模块、实际作物腾发量计算模块、作物系数基础值修订模块、优化求解模块、统计模块、作物系数计算模块、灌溉用水量计算模块、输入显示模块、控制模块均通信相连。
参数获取模块用于获取灌溉试验站多年逐日的日照时数、平均相对湿度、平均气温、最高气温、最低气温、平均风速气象资料以及试验推导的水稻作物系数kc;
生长期划分模块与参数获取模块通信相连,根据水稻生长情况,将水稻生长期划分为生长初期、快速生长期、生长中期、生长末期;
实际作物腾发量计算模块与参数获取模块相连,根据气象资料计算参考作物腾发量ET0,通过作物系数kc与参考作物腾发量乘积得到实际作物腾发量ETc:
式中,Rn为太阳净辐射;G为土壤热通量;γ为湿度表常数;u2为2m处风速;T为平均气温;ed为实际水汽压;ea为饱和水汽压;Δ为温度~饱和水汽压关系曲线T处的切线斜率。
作物系数基础值修订模块中将作物系数试验推导值推求的作物腾发量假定为实际作物腾发量,设定目标函数为作物系数基础修订值推求的作物腾发量与作物系数试验推导值推求的作物腾发量相对误差最小;
目标函数F:
约束条件:
ETc,i=Kc,i×ET0,i,
ET′c,i=K′c,i×ET0,i,
a≤K′cini≤b,K′cini≤K′cmid≤c,d≤K′cend≤K′cmid,
式中:ETc,i为作物系数试验推导值推求的日序数为i的作物腾发量;Kc,i为日序数为i的作物系数试验推导值;ET0,i为日序数为i的参考作物腾发量;K′c,i为基于基础修订值的日序数i作物系数理论推导值;K′cini为生长初期作物系数基础修订值,自变量;K'cmid为生长中期作物系数基础修订值,自变量;K'cend为生长后期末作物系数基础修订值,自变量;n0为生长初期起始日序数;n1为生长初期终止日序数;n2为快速生长期终止日序数;n3为生长中期终止日序数;n4为生长后期终止日序数;a和b为生长初期作物系数基础修订值下限和上限;c为生长中期作物系数基础修订值的上限;d为生长后期末作物系数基础修订值的下限。
优化求解模块与参数获取模块、生长期划分模块、实际作物腾发量计算模块均通信相连,基于作物系数基础值修订模块中的目标函数F和约束条件进行优化求解,得到生长初期作物系数基础修订值、生长中期作物系数基础修订值和生长后期末作物系数基础修订值;
具体地,优化求解模块采用粒子群算法进行优化求解:
vmin,d,vmax,d分别为最小和最大速度,
其中,1≤d≤D,1≤j≤M,(D=3,M为最大迭代次数),
粒子在t+1时刻的位置通过下式更新:
式中,r1,r2为均匀分布在(0,1)区间的随机数;ω为惯性权重;c1,c2为学习因子。
统计模块与参数获取模块、生长期划分模块均通信相连,统计水稻生长中期、生长末期的日平均风速u′2、日最低相对湿度RH′min。
作物系数计算模块与参数获取模块、生长期划分模块、实际作物腾发量计算模块、优化求解模块、统计模块均通信相连,根据以下公式计算出水稻不同生长阶段的作物系数:
其中,水稻生长中期、生长后期末水稻作物系数的修订公式为:
灌溉用水量计算模块与作物系数计算模块通信相连,基于水稻不同生长阶段的作物系数计算得到水稻灌溉用水量。
输入显示模块与参数获取模块、生长期划分模块、实际作物腾发量计算模块、作物系数基础值修订模块、优化求解模块、统计模块、作物系数计算模块、灌溉用水量计算模块、控制模块、灌溉装置均通信相连,用于让用户输入操作指令,并进行相应显示,例如,显示出各个区域水稻的气象资料信息,并显示出各区域水稻的灌溉用水量计算值,正在执行的灌溉操作,灌溉已用水量等,便于让工作人员监控和了解灌溉情况。
控制模块与参数获取模块、生长期划分模块、实际作物腾发量计算模块、作物系数基础值修订模块、优化求解模块、统计模块、作物系数计算模块、灌溉用水量计算模块、灌溉装置均通信相连,控制它们的运行。
以上实施例仅是对本发明技术方案所做的举例说明。本发明所涉及的水稻作物系数计算方法以及水稻灌溉系统并不仅仅限定于在以上实施例中所描述的内容,而是以权利要求所限定的范围为准。本发明所属领域技术人员在该实施例的基础上所做的任何修改或补充或等效替换,都在本发明的权利要求所要求保护的范围内。
Claims (7)
1.一种水稻作物系数计算方法,其特征在于,包括如下步骤:
步骤1.收集灌溉试验站多年逐日的日照时数、平均相对湿度、平均气温、最高气温、最低气温、平均风速气象资料以及试验推导的水稻作物系数kc;
步骤2.根据水稻生长情况,将水稻生长期划分为生长初期、快速生长期、生长中期、生长末期;
步骤3.根据气象资料计算参考作物腾发量ET0,通过作物系数kc与参考作物腾发量乘积得到实际作物腾发量ETc:
式中,Rn为太阳净辐射;G为土壤热通量;γ为湿度表常数;u2为2m处风速;T为平均气温;ed为实际水汽压;ea为饱和水汽压;Δ为温度~饱和水汽压关系曲线T处的切线斜率;
步骤4.修订不同生长期作物系数基础值
假定作物系数试验推导值推求的作物腾发量为实际作物腾发量,目标函数为作物系数基础修订值推求的作物腾发量与作物系数试验推导值推求的作物腾发量相对误差最小;
目标函数F:
约束条件:
ETc,i=Kc,i×ET0,i,
ET′c,i=K′c,i×ET0,i,
a≤K′cini≤b,K′cini≤K′cmid≤c,d≤K′cend≤K′cmid,
式中:ETc,i为作物系数试验推导值推求的日序数为i的作物腾发量;Kc,i为日序数为i的作物系数试验推导值;ET0,i为日序数为i的参考作物腾发量;K′c,i为基于基础修订值的日序数i作物系数理论推导值;K′cini为生长初期作物系数基础修订值,自变量;K′cmid为生长中期作物系数基础修订值,自变量;K′cend为生长后期末作物系数基础修订值,自变量;n0为生长初期起始日序数;n1为生长初期终止日序数;n2为快速生长期终止日序数;n3为生长中期终止日序数;n4为生长后期终止日序数;a和b为生长初期作物系数基础修订值下限和上限;c为生长中期作物系数基础修订值的上限;d为生长后期末作物系数基础修订值的下限;
步骤5.进行优化求解,得到生长初期作物系数基础修订值、生长中期作物系数基础修订值和生长后期末作物系数基础修订值;
步骤6.统计水稻生长中期、生长末期的日平均风速u′2、日最低相对湿度RH′min;
步骤7.采用以下公式计算得到水稻不同生长阶段的作物系数:
其中,水稻生长中期、生长后期末水稻作物系数的修订公式为:
2.根据权利要求1所述的水稻作物系数计算方法,其特征在于:
其中,在步骤2中,生长初期为从水稻播种到作物覆盖率接近10%,快速生长期为从覆盖率10%到大田作物充分覆盖,覆盖率达到70%~80%,生长中期为从充分覆盖到成熟期开始到叶片开始变黄,生长后期为从叶片开始变黄到生理成熟或收获。
4.一种水稻灌溉系统,其特征在于,包括:
灌溉装置;和
控制装置,基于权利要求1至3中任意一项所述的水稻作物系数计算方法得到水稻灌溉用水量,并根据水稻灌溉用水量控制所述灌溉装置对水稻进行灌溉。
5.根据权利要求4所述的水稻灌溉系统,其特征在于:
其中,所述控制装置包括:
参数获取模块,获取灌溉试验站多年逐日的日照时数、平均相对湿度、平均气温、最高气温、最低气温、平均风速气象资料以及试验推导的水稻作物系数kc;
生长期划分模块,与所述参数获取模块通信相连,根据水稻生长情况,将水稻生长期划分为生长初期、快速生长期、生长中期、生长末期;
实际作物腾发量计算模块,与所述参数获取模块相连,根据气象资料计算参考作物腾发量ET0,通过作物系数kc与参考作物腾发量乘积得到实际作物腾发量ETc:
式中,Rn为太阳净辐射;G为土壤热通量;γ为湿度表常数;u2为2m处风速;T为平均气温;ed为实际水汽压;ea为饱和水汽压;Δ为温度~饱和水汽压关系曲线T处的切线斜率;
作物系数基础值修订模块,假定作物系数试验推导值推求的作物腾发量为实际作物腾发量,目标函数为作物系数基础修订值推求的作物腾发量与作物系数试验推导值推求的作物腾发量相对误差最小;
目标函数F:
约束条件:
ETc,i=Kc,i×ET0,i,
ET′c,i=K′c,i×ET0,i,
a≤K′cini≤b,K′cini≤K′cmid≤c,d≤K′cend≤K′cmid,
式中:ETc,i为作物系数试验推导值推求的日序数为i的作物腾发量;Kc,i为日序数为i的作物系数试验推导值;ET0,i为日序数为i的参考作物腾发量;K′c,i为基于基础修订值的日序数i作物系数理论推导值;K′cini为生长初期作物系数基础修订值,自变量;K′cmid为生长中期作物系数基础修订值,自变量;K′cend为生长后期末作物系数基础修订值,自变量;n0为生长初期起始日序数;n1为生长初期终止日序数;n2为快速生长期终止日序数;n3为生长中期终止日序数;n4为生长后期终止日序数;a和b为生长初期作物系数基础修订值下限和上限;c为生长中期作物系数基础修订值的上限;d为生长后期末作物系数基础修订值的下限;
优化求解模块,与所述参数获取模块、所述生长期划分模块、所述实际作物腾发量计算模块均通信相连,基于作物系数基础值修订模块中的目标函数F和约束条件进行优化求解,得到生长初期作物系数基础修订值、生长中期作物系数基础修订值和生长后期末作物系数基础修订值;
统计模块,与所述参数获取模块、所述生长期划分模块均通信相连,统计水稻生长中期、生长末期的日平均风速u′2、日最低相对湿度RH′min;
作物系数计算模块,与所述参数获取模块、所述生长期划分模块、所述实际作物腾发量计算模块、所述优化求解模块、所述统计模块均通信相连,根据以下公式计算出水稻不同生长阶段的作物系数:
其中,水稻生长中期、生长后期末水稻作物系数的修订公式为:
灌溉用水量计算模块,与所述作物系数计算模块通信相连,基于水稻不同生长阶段的作物系数计算得到水稻灌溉用水量;
控制模块,与所述参数获取模块、所述生长期划分模块、所述实际作物腾发量计算模块、所述作物系数基础值修订模块、所述优化求解模块、所述统计模块、所述作物系数计算模块、所述灌溉用水量计算模块、所述灌溉装置均通信相连,控制它们的运行。
7.根据权利要求4的水稻灌溉系统,其特征在于,还包括:
输入显示模块,与所述参数获取模块、所述生长期划分模块、所述实际作物腾发量计算模块、所述作物系数基础值修订模块、所述优化求解模块、所述统计模块、所述作物系数计算模块、所述灌溉用水量计算模块、所述控制模块、所述灌溉装置均通信相连,用于让用户输入操作指令,并进行相应显示。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010446542.7A CN111640038B (zh) | 2020-05-25 | 2020-05-25 | 水稻作物系数计算方法以及水稻灌溉系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010446542.7A CN111640038B (zh) | 2020-05-25 | 2020-05-25 | 水稻作物系数计算方法以及水稻灌溉系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111640038A true CN111640038A (zh) | 2020-09-08 |
CN111640038B CN111640038B (zh) | 2022-03-04 |
Family
ID=72330910
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010446542.7A Active CN111640038B (zh) | 2020-05-25 | 2020-05-25 | 水稻作物系数计算方法以及水稻灌溉系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111640038B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113378381A (zh) * | 2021-06-09 | 2021-09-10 | 水利部交通运输部国家能源局南京水利科学研究院 | 基于气温分布和地表热量计算冬小麦作物系数的方法 |
CN114002951A (zh) * | 2021-09-16 | 2022-02-01 | 江苏农林职业技术学院 | 一种水稻硬地育秧模糊控制灌溉方法 |
CN115067162A (zh) * | 2022-07-27 | 2022-09-20 | 天津市农业科学院 | 一种提高口感型番茄可溶性固形物的基质栽培管理模式 |
CN118521431A (zh) * | 2024-07-23 | 2024-08-20 | 湖北省水利水电科学研究院 | 灌区多级渠道自适应的实时动态配水方法及系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104521699A (zh) * | 2014-11-18 | 2015-04-22 | 华北水利水电大学 | 田间智能灌溉在线控制管理方法 |
CN106713342A (zh) * | 2017-01-06 | 2017-05-24 | 武汉大学 | 一种基于b/s架构的灌区配水综合管理系统及方法 |
CN107862148A (zh) * | 2017-11-16 | 2018-03-30 | 清华大学 | 一种基于遥感数据的参考作物腾发量计算方法 |
CN107945042A (zh) * | 2017-11-29 | 2018-04-20 | 上海华维节水灌溉股份有限公司 | 一种作物生长灌溉决策控制系统 |
WO2018173045A1 (en) * | 2017-03-20 | 2018-09-27 | Supplant Ltd. | Systems and methods for planning crop irrigation |
WO2019145895A1 (en) * | 2018-01-24 | 2019-08-01 | The State Of Israel, Ministry Of Agriculture & Rural Development Agricultural Research Organization | Method and system for estimating crop coefficient and evapotranspiration of crops based on remote sensing |
CN110209077A (zh) * | 2019-05-05 | 2019-09-06 | 武汉大学 | 基于互联网的远程灌排系统实时动态监测系统 |
CN110754344A (zh) * | 2019-10-08 | 2020-02-07 | 京蓝物联技术(北京)有限公司 | 基于气象预告的灌溉决策方法及装置 |
-
2020
- 2020-05-25 CN CN202010446542.7A patent/CN111640038B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104521699A (zh) * | 2014-11-18 | 2015-04-22 | 华北水利水电大学 | 田间智能灌溉在线控制管理方法 |
CN106713342A (zh) * | 2017-01-06 | 2017-05-24 | 武汉大学 | 一种基于b/s架构的灌区配水综合管理系统及方法 |
WO2018173045A1 (en) * | 2017-03-20 | 2018-09-27 | Supplant Ltd. | Systems and methods for planning crop irrigation |
CN107862148A (zh) * | 2017-11-16 | 2018-03-30 | 清华大学 | 一种基于遥感数据的参考作物腾发量计算方法 |
CN107945042A (zh) * | 2017-11-29 | 2018-04-20 | 上海华维节水灌溉股份有限公司 | 一种作物生长灌溉决策控制系统 |
WO2019145895A1 (en) * | 2018-01-24 | 2019-08-01 | The State Of Israel, Ministry Of Agriculture & Rural Development Agricultural Research Organization | Method and system for estimating crop coefficient and evapotranspiration of crops based on remote sensing |
CN110209077A (zh) * | 2019-05-05 | 2019-09-06 | 武汉大学 | 基于互联网的远程灌排系统实时动态监测系统 |
CN110754344A (zh) * | 2019-10-08 | 2020-02-07 | 京蓝物联技术(北京)有限公司 | 基于气象预告的灌溉决策方法及装置 |
Non-Patent Citations (4)
Title |
---|
ZANOTELLI, DAMIANO等: "Evapotranspiration and crop coefficient patterns of an apple orchard in a sub-humid environment", 《AGRICULTURAL WATER MANAGEMENT》 * |
吴慧蓉等: "作物系数影响因素分析", 《节水灌溉》 * |
汤鹏程: "西藏高海拔地区ET0计算公式试验率定与青稞作物系数推求", 《中国博士学位论文全文数据库(农业科技辑)》 * |
潘少斌等: "湖北省早稻灌溉定额修订方法研究", 《节水灌溉》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113378381A (zh) * | 2021-06-09 | 2021-09-10 | 水利部交通运输部国家能源局南京水利科学研究院 | 基于气温分布和地表热量计算冬小麦作物系数的方法 |
CN113378381B (zh) * | 2021-06-09 | 2022-01-28 | 水利部交通运输部国家能源局南京水利科学研究院 | 基于气温分布和地表热量计算冬小麦作物系数的方法 |
CN114002951A (zh) * | 2021-09-16 | 2022-02-01 | 江苏农林职业技术学院 | 一种水稻硬地育秧模糊控制灌溉方法 |
CN114002951B (zh) * | 2021-09-16 | 2023-12-29 | 江苏农林职业技术学院 | 一种水稻硬地育秧模糊控制灌溉方法 |
CN115067162A (zh) * | 2022-07-27 | 2022-09-20 | 天津市农业科学院 | 一种提高口感型番茄可溶性固形物的基质栽培管理模式 |
CN118521431A (zh) * | 2024-07-23 | 2024-08-20 | 湖北省水利水电科学研究院 | 灌区多级渠道自适应的实时动态配水方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN111640038B (zh) | 2022-03-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111640038B (zh) | 水稻作物系数计算方法以及水稻灌溉系统 | |
CN107945042B (zh) | 一种作物生长灌溉决策控制系统 | |
CN106718695B (zh) | 一种智能节水灌溉物联网控制系统 | |
CN106359005B (zh) | 一种间作农田的自动灌溉装置及自动灌溉方法 | |
JP6609800B2 (ja) | 生育診断方法、生育診断装置、データ測定装置、生育診断プログラム、生育補助方法、生育補助装置及び生育補助プログラム | |
CN109447426B (zh) | 基于作物需水机理的灌溉需水对变化环境的响应分析方法 | |
CN111557158A (zh) | 智能灌溉控制方法及系统 | |
Jia et al. | Optimization of an extreme learning machine model with the sparrow search algorithm to estimate spring maize evapotranspiration with film mulching in the semiarid regions of China | |
CN110896836B (zh) | 一种无土栽培营养液管控方法及系统 | |
CN110432046B (zh) | 一种温室内的智能灌溉系统 | |
CN115310680A (zh) | 一种番茄种苗模型建模与生长预测方法 | |
CN117981663A (zh) | 水稻田自动化节水精准灌溉系统 | |
CN105993720B (zh) | 一种日光温室基质袋培作物灌水量的模拟计算方法 | |
CN108235887A (zh) | 一种温室中宽垄覆膜沟灌的方法 | |
CN112273026B (zh) | 一种潮汐育苗营养液精准决策方法和系统 | |
CN115989763B (zh) | 温室灌溉控制方法、装置、系统及设备 | |
CN109934400B (zh) | 基于改进神经网络的集雨调亏作物需水量预测方法 | |
Van Keulen et al. | Potential wheat yields in Zambia—a simulation approach | |
CN116127705A (zh) | 温室作物生长模拟方法及装置 | |
CN114637353A (zh) | 基于多环境因素分析的农业灌溉控制方法、系统及终端 | |
CN114568239A (zh) | 一种棉花高温热害预测方法 | |
CN110574642B (zh) | 通过调控通风模式提高春季日光温室番茄产量的方法 | |
CN113615545A (zh) | 一种提供水田作物精准灌溉的装置及方法 | |
KR102715861B1 (ko) | 기상데이터를 이용한 작물수분스트레스지수(cwsi) 예측방법 | |
CN117193437B (zh) | 一种林业育苗培育系统及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |