基于全局敏感性分析和贝叶斯方法的WOFOST作物模型参数优化.pdf

《基于全局敏感性分析和贝叶斯方法的WOFOST作物模型参数优化.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第32卷第2期农业工程学报·、，r01．32No．2基于全局敏感性分析和贝叶斯方法的WOFOST作物模型参数优化何亮1，侯英雨1，赵刚2，邬定荣3，于强4,5(1．国家气象中心，北京100081；2．德国波恩大学作物科学研究组，波恩D．53115；3．中国气象科学研究院，北京100081；4．中国科学院地理科学与资源研究所陆地水循环与地表过程重点实验室，北京100101；5．澳大利亚悉尼科技大学生命科学学院，悉尼2007)摘要：作物模型参数的敏感性分析、标定和验证可以提高模型的效率和精准度，进而为模型应用做好准备工作。该研究结

2、合参数全局敏感性分析方法以及贝叶斯后验估计理论的马尔科夫蒙特卡洛(MarkovChainMonteCarlo，MCMC)方法，以华北栾城站三年的冬小麦观测数据(叶面积和地上生物量)为参照，对WOFOST模型的55个品种参数进行了敏感性分析、筛选和优化。发现：1)对叶面积影响较大的参数为：生育期为0、O．5、0．6和O．75时的比叶面积、生育期为1．5时的最大光合速率、叶面积指数最大增长率；对地上干物质影响较大的参数为：生育期为1。5时的最大光合速率、生育期为0时的比叶面积、35℃时叶面积的生命周期、生育期为0时的散射消光系数、

3、生育期为1．8时的最大光合速率、储存器官的同化物转换效率。2)潜在和雨养产量水平下，最大叶面积和地上生物量对参数的敏感性差异不大。3)马尔科夫蒙特卡洛方法(MCMC)可以对WOFOST模型品种参数较好地优化；设计的3种校正．验证方案中，第1种方案(用1998—1999年作为校正年份，1999--2000年，2000--2001年作为验证年份)模拟效果最好。4)优化后的参数，模型对潜在产量水平模拟较好，一致性指数均大于0．9，相对均方根误差小于20％；而对有水分胁迫的雨养情况下比潜在产量水平的模拟结果差，表明模型对水分胁迫的模拟

4、不足。该研究为WOFOST模型区域应用和模型调整优化提供科学理论依据。关键词：模型；作物；优化；WOFOST；全局敏感性分析；MCMC；模型参数优化doi：10．11975／j．issn．1002-6819．2016．02．025中图分类号：S512．1文献标志码：A文章编号：1002—6819(2016)一02—0169—11何亮，侯英雨，赵刚，邬定荣，于强．基于全局敏感性分析和贝叶斯方法的WOFOST作物模型参数优化[J]．农业工程学报，2016，32(2)：169—179．doi：10．11975／j．issn．1002

5、—6819．2016．02．025http：／／www．tcsae．orgHeLiang，HouYingyu，ZhaoGang，WuDingrong，YuQiang．ParametersoptimizationofWOFOSTmodelbyintegrationofglobalsensitivityanalysisandBayesiancalibrationmethod[J】．TransactionsoftheChineseSocietyofAgriculturalEngineering(TransactionsoftheCS

6、AZ)，2016，32(2)：169—179．(inChinesewithEnglishabstract)doi：10．11975，j．issn．1002—6819．2016．02．025http：／／www．tcsae．org0引言模型的标定(calibration)和验证(validation)等准备工作是基于机理作物模型应用的前提。对于基于机理的作物模型，参数众多。作物模型的很多参数(如比叶面积)不能通过田间观测直接获得，而是需要通过可观测的变量(如叶面积指数)进行反演。通过准确的方法获得模型参数是进行模型应用、提高模型

7、可预报性的前提【l】。根据观测值反推参数属于参数估计问题。对于线性方程或者简单的非线性方程可以用最／bz．乘法来解决。基于过程的作物模型，刻画了光合、干物质分配、土壤水分运移和蒸发等众多生物物理过程，往往包含的方程较多，模型的非线性效应很明显。对于一般的最小二乘法和非线性参数估计方法难以获得全局的最优解。一种收稿日期：2015．11—13修订日期：2015．12—21基金项目：公益性行业(气象)科研专项(GⅥ{Y201306052，GYHY201506001)作者简介：何亮，男(汉族)，湖南永兴县人，博士，主要从事作物模型、农

8、业气象和全球变化研究。北京国家气象中心，100081。Emaihheliang_hello@163．com常用的模型参数标定是基于蒙特卡洛试错法。这种方法根据以往的经验或者参数文献参考值重复和随机地选择参数，使得观测值和模拟值拟合程度指标：如决定系数(尺2)、相关系数、均根方