基于APSIM模型旱地小麦叶面积指数相关参数的优化

doi:10.3864/j.issn.0578-1752.2019.12.004

摘要/Abstract

摘要：

【目的】模型的有效应用依赖于参数的快速、准确估算。本研究拟解决作物生长模型参数本土化率定过程中运算量大、耗时长、精度低、效率低的问题。【方法】依据甘肃省定西市安定区2个试验点（李家堡镇麻子川村和凤翔镇安家沟村）,多年（2002—2005年和2015—2017年）大田试验数据以及定西市安定区1971—2017年气象资料,利用混合蛙跳算法智能的迭代搜索原理,对APSIM模型旱地小麦叶面积指数相关参数进行了优化,并采用相关性分析方法对模型校正结果进行检验。【结果】利用青蛙群体即相对独立又合作协调的子群内局部深度搜索与子群间全局信息交流生物进化学习策略,有效提高了运算的速度,实现了对APSIM模型中与旱地小麦叶面积指数相关参数的快速、准确估算。相关参数主要包括：主茎上节出现所需的热时间间隔、小麦出苗后初始化的节数、小麦出苗后初始化的叶片数、小麦出苗后初始化的叶面积指数、某日正在生长的节数和最大比叶面积。分别使用穷举试错法所得参数值和混合蛙跳算法所得参数值模拟叶面积指数,参数优化后,叶面积指数模拟值和实测值之间的RMSE（root mean square error）平均值由0.069降低到0.027,NRMSE（normalized root mean square error）平均值由8.09%降低到4.56%,M_E（model effective index）平均值由0.979提高到0.993。【结论】相对于参数率定常用穷举试错法,混合蛙跳算法具有自发学习特征的智能迭代行为,实现了参数的自动率定,提高了效率。基于该算法进行APSIM模型旱地小麦叶面积指数相关参数的优化,使得模型对叶面积指数的模拟精度显著提高,证明该算法的使用对作物生理生态系统复杂模型的校正效果良好,为改善模型参数率定过程存在的运算量大、耗时长、精度低、效率低的缺点提供了一种行之有效的方法。

关键词: APSIM, 小麦, 叶面积指数, 参数优化, 混合蛙跳算法

Abstract:

【Objective】The effective application of the model depends on the fast and accurate estimation of parameters. The current problems in the calibration of crop growth model parameters include large data volume, long time consumption, lack of precision, and low efficiency. The study tried to solve the problems. 【Method】Based on the field experimental data of two experimental sites (Mazichuan village, Lijiabaotown and Anjiagou village, Fengxiangtown) in Andingdistrict, Dingxi city in multiple years (2002-2005 and 2015-2017) and the meteorological data in Andingdistrict, Dingxicity from 1971 to 2017, the parameters related to dryland wheat leaf area index (LAI) in the APSIM (agricultural production systems simulator) model were optimized with the intelligent iteration search principle of shuffled frog leaping algorithm (SFLA) and tested by the correlation analysis method. 【Result】The biological evolution learning strategy of local depth search within sub-group and global information communication between sub-group in frog population, which was relatively independent and coordinated, was used to effectively improve the speed of calculation and realize the fast and accurate estimation of the parameters related to dryland wheat LAI in the APSIM model. The related parameters mainly included: The required thermal time interval for node appearance on the main stem, the initial node number at emergence, the initial leaf number at emergence, the initial leaf area index at emergence, the growing node number, and the maximum specific leaf area. LAI was respectively simulated by using the parameters based on the trial and error method and based on SFLA. After parameter optimization, the root mean square error (RMSE) between simulated and measured wheat LAI reduced from 0.069 to 0.027, the normalized root mean square error (NRMSE) decreased from 8.09% to 4.56%, and the model effective index (M_E) increased from 0.979 to 0.993. 【Conclusion】Compared with the trial and error method, which was usually used in the calibration of APSIM model, the intelligent iterative behavior with spontaneous learning characteristics based on the SFLA could realize automatic calibration of the parameters and improve the efficiency. The parameters estimated based on the SFLA could remarkably improve the simulation accuracy of wheat LAI. The application of SFLA was effective in calibrating crop models involving complex eco-physiological processes, and it could provide an effective parameter optimization method for improving the disadvantages in the model parameter calibration process include large data volume, long time consumption, lack of precision, and low efficiency.

Key words: APSIM, wheat, leaf area index, parameter optimization, shuffled frog leaping algorithm

聂志刚,李广,王钧,马维伟,雒翠萍,董莉霞,逯玉兰. 基于APSIM模型旱地小麦叶面积指数相关参数的优化[J]. 中国农业科学, 2019, 52(12): 2056-2068.

NIE ZhiGang,LI Guang,WANG Jun,MA WeiWei,LUO CuiPing,DONG LiXia,LU YuLan. Parameter Optimization for the Simulation of Leaf Area Index of Dryland Wheat with the APSIM Model[J]. Scientia Agricultura Sinica, 2019, 52(12): 2056-2068.

图/表 8

图1

图2

图3

表1

表2

图4

图5

表3

参考文献 41

[1]	曹卫星, 朱艳, 田永超, 姚霞, 刘小军 . 数字农作技术研究的若干进展与发展方向. 中国农业科学, 2006,39(2):281-288.
	CAO W X, ZHU Y, TIAN Y C, YAO X, LIU X J . Research progress and prospect of digital farming techniques. ScientiaAgriculturaSinica, 2006,39(2):281-288. (in Chinese)
[2]	沈禹颖, 南志标 , BILL B, MICHAEL R, 陈文, 邵新庆. APSIM模型的发展与应用. 应用生态学报, 2002,13(8):1027-1032.
	SHEN Y Y, NAN Z B, BILL B, MICHAEL R, CHEN W, SHAO X Q . Development of APSIM (Agricultural Production Systems Simulator) and its application. Chinese Journal of Applied Ecology, 2002,13(8):1027-1032. (in Chinese)
[3]	JONES J W, HOOGENBOOM G, PORTER C H, BOOTE K J, BATCHELOR W D . The DSSAT cropping system model. European Journal of Agronomy, 2003,18(3):235-265. doi: 10.1016/S1161-0301(02)00107-7
[4]	WILLIAMS J, JONES C, KINIRY J, SPANEL D . The EPIC crop growth model. Transactions of the ASAE, 1989,32(2):497-511. doi: 10.13031/2013.31032
[5]	DIEPEN C, WOLF J, KEULEN H, RAPPOLDT C . Wofost: A simulation model of crop production. Soil Use Manage, 1989,5(1):16-24. doi: 10.1111/sum.1989.5.issue-1
[6]	庄嘉祥, 姜海燕, 刘蕾蕾, 王芳芳, 汤亮, 朱艳, 曹卫星 . 基于个体优势遗传算法的水稻生育期模型参数优化. 中国农业科学, 2013,46(11):2220-2231. doi: 10.3864/j.issn.0578-1752.2013.11.005
	ZHUANG J X, JIANG H Y, LIU L L, WANG F F, TANG L, ZHU Y, CAO W X . Parameters optimization of rice development stages model based on individual advantages genetic algorithm. ScientiaAgriculturaSinica, 2013,46(11):2220-2231.(in Chinese) doi: 10.3864/j.issn.0578-1752.2013.11.005
[7]	房全孝 . 根系水质模型中土壤与作物参数优化及其不确定性评价. 农业工程学报, 2012,28(10):118-123.
	FANG Q X . Optimizing and uncertainty evaluation of soil and crop parameters in root zone water quality model. Transactions of the Chinese Society of AgriculturalEngineering, 2012,28(10):118-123.(in Chinese)
[8]	刘志娟, 杨晓光, 王静, 吕硕, 李克南, 荀欣, 王恩利 . APSIM玉米模型在东北地区的适应性. 作物学报, 2012,38(4):740-746. doi: 10.3724/SP.J.1006.2012.00740
	LIU Z J, YANG X G, WANG J, LÜ S, LI K N, XUN X, WANG E L . Adaptability of APSIM maize model in Northeast China. ActaAgronomicaSinica, 2012,38(4):740-746.(in Chinese) doi: 10.3724/SP.J.1006.2012.00740
[9]	聂志刚, 李广, 雒翠萍, 马维伟, 代永强 . 利用混合蛙跳算法优化基于APSIM的旱地小麦产量形成模型参数. 作物学报, 2018,44(8):1229-1236.
	NIE Z G, LI G, LUO C P, MA W W, DAI Y Q . Parameter optimization in APSIM-based simulation model for yield formation of dryland wheat using shuffled frog leaping algorithm. ActaAgronomicaSinica, 2018,44(8):1229-1236.(in Chinese)
[10]	代永强 . 混合蛙跳算法的改进与应用[D]. 兰州: 甘肃农业大学, 2011.
	DAI Y Q . The improvement and application of shuffled frog leaping algorithm[D]. Lanzhou: Gansu Agricultural University, 2011. (in Chinese)
[11]	MANSOURI M, DESTAIN M F . An improved particle filtering for time-varying nonlinear prediction of biomass and grain protein content. Computers and Electronics in Agriculture, 2015,114:145-153. doi: 10.1016/j.compag.2015.04.006
[12]	CALMON M A, JONES J W, SHINDE D, SPECHT J E . Estimating parameters for soil water balance models using adaptive simulated annealing. Applied Engineering in Agriculture, 1999,15(6):703-713. doi: 10.13031/2013.5841
[13]	DAI C N, YAO M, XIE Z J, CHEN C H, LIU J G . Parameter optimization for growth model of green house crop using genetic algorithms. Applied Soft Computing, 2009,9(1):13-19. doi: 10.1016/j.asoc.2008.02.002
[14]	刘铁梅, 王燕, 邹薇, 孙东发, 汤亮, 曹卫星 . 大麦叶面积指数模拟模型. 应用生态学报, 2010,21(1):121-128.
	LIU T M, WANG Y, ZUO W, SUN D F, TANG L, CAO W X . Simulation model of barley leaf area index. Chinese Journal of Applied Ecology, 2010,21(1):121-128. (in Chinese)
[15]	汪定伟 . 智能优化算法. 北京: 高等教育出版社, 2007: 136-148.
	WANG D W. Intelligent Optimization Methods. Beijing: Higher Education Press, 2007: 136-148. (in Chinese)
[16]	孙建平, 闫蕾, 李妍, 张婧 . 基于改进遗传算法的模糊PID控制器设计. 仪器仪表学报, 2006,27(6):1991-1992.
	SUN J P, YAN L, LI Y, ZHANG J . Design of fuzzy PID controllers based on improved genetic algorithms. Chinese Journal of Scientific Instrument., 2006,27(6):1991-1992. (in Chinese)
[17]	EUSUFF M M, LANSEY K E . Optimization of water distribution network design using the shuffled frog leaping algorithm. Journal of Water Resources Planning and Management, 2003,129(3):210-225. doi: 10.1061/(ASCE)0733-9496(2003)129:3(210)
[18]	李广, 黄高宝 , WILLIAM B, 陈文. APSIM模型在黄土丘陵沟壑区不同耕作措施中的适用性. 生态学报, 2009,29(5):2655-2663.
	LI G, HUANG G B, WILLIAM B, CHEN W . Adaptation research of APSIM model under different tillage systems in the Loess hill-gullied region. ActaEcologicaSinica, 2009,29(5):2655-2663. (in Chinese)
[19]	李广, 黄高宝 . 基于APSIM 模型的降水量分配对旱地小麦和豌豆产量影响的研究. 中国生态农业学报, 2010,18(2):342-347.
	LI G, HUANG G B . Determination of the effect of precipitation distribution on yield of wheat and pea in dryland using APSIM. Chinese Journal of Eco-Agriculture, 2010,18(2):342-347. (in Chinese)
[20]	李广, 李玥, 黄高宝, 罗珠珠, 王琦, 刘强, 燕振刚, 赵有益 . 基于APSIM模型旱地春小麦产量对温度和CO₂浓度升高的响应. 中国生态农业学报, 2012,20(8):1088-1095.
	LI G, LI Y, HUANG G B, LUO Z Z, WANG Q, LIU Q, YAN Z G, ZHAO Y Y . Response of dryland spring wheat yield to elevated CO₂ concentration and temperature by APSIM model. Chinese Journal of Eco-Agriculture, 2012,20(8):1088-1095.(in Chinese)
[21]	李玥, 牛俊义, 谢亚萍, 吴兵, 高珍妮, 刘栋, 剡斌 . 基于APSIM的油用亚麻叶面积指数模型构建. 中国油料作物学报, 2015,37(3):329-335. doi: 10.7505/j.issn.1007-9084.2015.03.011
	LI Y, NIU J Y, XIE Y P, WU B, GAO Z N, LIU D, YAN B . Simulation of oilseed flax leaf area index based on APSIM. Chinese Journal of Oil Crop Sciences, 2015,37(3):329-335. (in Chinese) doi: 10.7505/j.issn.1007-9084.2015.03.011
[22]	聂志刚, 李广 . 基于APSIM的旱地小麦叶面积指数模拟模型构建. 干旱地区农业研究, 2013,31(4):94-98.
	NIE Z G, LI G . Modeling of APSIM-based simulation of leaf area index of wheat in dryland. Agricultural Research in the Arid Areas, 2013,31(4):94-98.(in Chinese)
[23]	GOUDRIAAN J , VAN LAAR H H. Modelling Potential Crop Growth Processes: Textbook with Exercises. The Netherlands: Kluwer Academic Publishers, 1994: 7-28.
[24]	柏军华, 王克如, 初振东, 陈兵, 李少昆 . 叶面积测定方法的比较研究. 石河子大学学报(自然科学版), 2005,23(2):216-218.
	BEI J H, WANG K R, CHU Z D, CHEN B, LI S K . Comparative study on the measure methods of the leaf area. Journal of ShiheziUniversity(Natural Science), 2005,23(2):216-218. (in Chinese)
[25]	鲍士旦 . 土壤农化分析. 第三版. 北京: 中国农业出版社, 2000: 263-268.
	BAO S D. Soil Agricultural Chemistry Analysis.3rd edition. Beijing: China Agriculture Press, 2000: 263-268. (in Chinese)
[26]	ZHENG B Y, CHENU K, DOHERTY A, CHAPMAN S . The APSIM-Wheat Module (7.5 R3008) Documentation. Toowoomba: APSRU, 2014: 20-25.
[27]	RITCHIE J T, GODWIN D C, OTTER-NACKE S . CERES-wheat: A user-oriented wheat yield model.MichiganStateUniversity, East Lansing: AGRISTARS Publication, 1985: 252.
[28]	SINCLAIR T R . Water and nitrogen limitations in soybean grain productionI. Model development. Field Crops Research, 1986,15(2):125-141. doi: 10.1016/0378-4290(86)90082-1
[29]	MONTEITH J L . How do crops manipulate water supply and demand. Philosophical Transactions of the Royal Society B-Biological Sciences, 1986,316:245-259. doi: 10.1098/rsta.1986.0007
[30]	陈国庆, 朱艳, 曹卫星 . 小麦叶鞘和节间生长过程的模拟研究. 麦类作物学报, 2005,25(1):71-74. doi: 10.7606/j.issn.1009-1041.2005.01.018
	CHEN G Q, ZHU Y, CAO W X . Modeling leaf sheath and internode growth dynamics in wheat. Journal of Triticeae Crops, 2005,25(1):71-74.(in Chinese) doi: 10.7606/j.issn.1009-1041.2005.01.018
[31]	杨文雄 . 中国西北春小麦. 北京: 中国农业出版社, 2016: 155-161.
	YANG W X. Spring Wheat in Northwest China. Beijing: China Agriculture Press, 2016: 155-161. (in Chinese)
[32]	刘铁梅, 曹卫星, 罗卫红, 郭文善 . 小麦叶面积指数的模拟模型研究. 麦类作物学报, 2001,21(2):38-41. doi: 10.7606/j.issn.1009-1041.2001.02.039
	LIU T M, CAO W X, LUO W H, GUO W S . Simulation on leaf area index in wheat. Journal of Triticeae, 2001,21(2):38-41. (in Chinese) doi: 10.7606/j.issn.1009-1041.2001.02.039
[33]	ZHANG X C . Calibration, refinement, and application of the WEPP model for simulation climatic impact on wheat production. Transactions of the ASAE, 2004,47(4):1075-1085. doi: 10.13031/2013.16580
[34]	何亮, 赵刚, 靳宁, 庄伟, 于强 . 不同气候区和不同产量水平下APSIM-Wheat模型的参数全局敏感性分析. 农业工程学报, 2015,31(14):148-157.
	HE L, ZHAO G, JIN N, ZHUANG W, YU Q . Global sensitivity analysis of APSIM-Wheat parameters in different climate zones and yield levels. Transactions of the Chinese Society of Agricultural Engineering, 2015,31(14):148-157.(in Chinese)
[35]	董朝阳, 刘志娟, 杨晓光 . 北方地区不同等级干旱对春玉米产量影响. 农业工程学报, 2015,31(11):157-164.
	DONG C Y, LIU Z J, YANG X G . Effect of different grade drought on grain yield of spring maize in Northern China. Transactions of the Chinese Society of Agricultural Engineering, 2015,31(11):157-164.(in Chinese)
[36]	ASSENG S, BAR-TAL A, BOWDEN J W, KEATING B A, HERWAARDEN A V, PALTA J A, HUTH N I, PROBERT M E . Simulation of grain protein content with APSIM-N wheat. European Journal Agronomy, 2002,16(1):25-42. doi: 10.1016/S1161-0301(01)00116-2
[37]	刘铁梅, 邹薇, 刘铁芳, 谢国生, 吴江生, 曹卫星, 李卫平, 曹凑贵, 傅廷栋 . 不同冬油菜品种比叶面积的多因子分析. 作物学报, 2006,32(7):1083-1089.
	LIU T M, ZOU W, LIU T F, XIE G S, WU J S, CAO W X, LI W P, CAO C G, FU T D . Analyses of multiple agronomic factors relative to specific leaf area in winter rape cultivars. ActaAgronomicaSinica, 2006,32(7):1083-1089. (in Chinese)
[38]	张怀志, 曹卫星, 周治国, 朱艳, 张立桢 . 棉花适宜叶面积指数的动态知识模型. 棉花学报, 2003,15(3):151-154. doi: 10.3969/j.issn.1002-7807.2003.03.004
	ZHANG H Z, CAO W X, ZHOU Z G, ZHU Y, ZHANG L Z . A dynamic knowledge model for optimail LAI in cotton. Cotton Science, 2003,15(3):151-154. (in Chinese) doi: 10.3969/j.issn.1002-7807.2003.03.004
[39]	GARY C, BARCZI J F, BERTIN N, TCHAMITCHIAN M . Simulation of individual organ growth and development on a tomato plant: A model and a user-friendly interface. ActaHorticulturae, 1995,399:199-205.
[40]	李广 . APSIM模型模拟与应用研究—基于旱地小麦、豌豆保护性耕作定位试验[D]. 兰州: 甘肃农业大学, 2006.
	LI G . Study on simulation and application of APSIM model-based on the conservation tillage of dryland wheat and pea[D]. Lanzhou: Gansu Agricultural University, 2006. (in Chinese)
[41]	贺力霞 . 基于动态客流的城市轨道交通列车牵引能耗仿真及优化[D]. 北京: 北京交通大学, 2018.
	HE L X . Simulation and optimization of traction energy consumption of urban rail transit based on dynamic passenger flow[D]. Beijing: BeijingJiaotongUniversity, 2018. (in Chinese)

项目 Item	土层深度 Soil depth
项目 Item	5 cm	10 cm	30 cm	50 cm	80 cm	110 cm	140 cm	170 cm	200 cm
容重 BD (g·cm^-3）	1.29	1.23	1.32	1.20	1.14	1.14	1.13	1.12	1.11
萎蔫系数 WC (mm·mm^-1)	0.08	0.08	0.08	0.08	0.09	0.09	0.11	0.13	0.13
最大持水量 DUL (mm·mm^-1)	0.27	0.27	0.27	0.27	0.26	0.27	0.26	0.26	0.26
饱和水分含量 SM (mm·mm^-1)	0.46	0.49	0.45	0.50	0.52	0.52	0.48	0.53	0.53
风干系数 CA (mm·mm^-1)	0.01	0.01	0.05	0.07	0.07	0.07	0.07	0.07	0.07
土壤导水率 SWC (mm·h^-1)	0.60	0.60	0.60	0.60	0.60	0.60	0.60	0.60	0.60
小麦有效水分下限 LLW (mm·mm^-1)	0.09	0.09	0.09	0.09	0.09	0.10	0.11	0.13	0.15

参数 Parameter	默认值 Default value	优化值 Optimized value
主茎上节出现所需的热时间间隔 The required thermal timeinterval for node appearance on the main stem (°C·d)	95.0	98.5
小麦出苗后初始化的节数The initial node number at emergence	2.00	1.79
小麦出苗后初始化的叶片数The initial leaf number at emergence	2.00	1.66
小麦出苗后初始化的叶面积指数The initial leaf area index at emergence	0.080	0.074
某日正在生长的节数The growing node number	2.00	1.86
最大比叶面积The maximum specific leaf area (mm²·g^-1)	26000	h_SLA(LAI_d)

模型参数 Model parameter	麻子川村Mazichuan			安家沟村Anjiagou
模型参数 Model parameter	RMSE	NRMSE (%)	M_E	RMSE	NRMSE(%)	M_E
默认值 Default value	0.047	7.06	0.983	0.090	9.12	0974
优化值 Optimized value	0.031	4.53	0.993	0.023	4.59	0.992