基于VIP模型评估不同全局敏感性分析方法有效性及效率  

Effectiveness and efficiency of different global sensitivity analysis methods based on the VIP model

作　　者：何力鸿 莫兴国[1,3] 刘苏峡 He Lihong;Mo Xingguo;Liu Suxia(Key Laboratory of Water Cycle and Related Land Surface Processes,Institute of Geographical Sciences and Natural Resources Research,Chinese Academy of Sciences,Beijing 100101,China;University of Chinese Academy of Sciences,Beijing 100049,China;School of Natural Resources and Environment/Sino-Danish College,University of Chinese Academy of Sciences,Beijing 100049,China)

机构地区：[1]中国科学院地理科学与资源研究所陆地水循环及地表过程重点实验室,北京100101 [2]中国科学院大学,北京100049 [3]中国科学院大学资环学院/中丹学院,北京100049

出　　处：《农业工程学报》2020年第16期78-84,共7页Transactions of the Chinese Society of Agricultural Engineering

基　　金：中国科学院战略先导科技专项(A类)(XDA20040301);国家自然科学基金(41471026,41790424)。

摘　　要：参数敏感性分析可提高复杂过程模型参数优化及模型应用效率。为评估不同全局敏感性分析方法在筛选农田生态系统模型敏感参数时的有效性及其效率,该研究以VIP(Vegetation Interface Processes)模型模拟华北平原土壤硝态氮为例,对比分析PSUADE(Problem Solving Environment for Uncertainty Analysis and Design Exploration)提供的8种敏感性分析方法在筛选华北平原农田土壤硝态氮敏感参数时的有效性及其效率。结果表明:在验证敏感性分析方法有效性时,Spearman秩相关系数(Spearman’s correlation coefficient,SPEA)法和Gaussian Process(GP)法与其他方法的敏感参数筛选的结果差异较大;多元自适应回归样条(Multivariate Adaptive Regression Splines,MARS)、Delta Test(DT)、Sum of Trees(SOT)法、McKay法、Morris法和Sobol’法能有效筛选出敏感参数。在分析敏感性方法效率时发现,正交阵列(Orthogonal Array,OA)和基于拉丁超立方的正交阵列(Orthogonal Array-Based Latin Hypercubes,OALH)抽样方式最适于多元自适应回归样条法(Multivariate Adaptive Regression Splines,MARS),所需样本量为361。蒙特卡罗(Monte Carlo,MC)抽样方式最适合DT和SOT敏感性分析方法,所需最小样本量分别为425和510;与其他抽样方式相比,OALH抽样更适合McKay敏感性分析方法,需要样本量最少;Morris法和Sobol’法所需最小样本量分别为340和810。对于复杂过程模型,可先选用定性敏感性分析方法以较低的计算成本选择敏感的参数,再进行模型参数优化和不确定性分析。Global sensitivity analysis generally refers to a sort of numerical approach for the evaluation of the variation in the output of a model,as the input parameter varies.In ecosystem model based on dynamic processes,the parameter verification can be difficult in the model application,due to a complex structure,multiple input parameters,and strong spatial variability.In this case,the parameter sensitivity analysis can effectively identify the factors of main effect,and thereby to improve significantly optimization of parameters and models.Given a fixed value to an insensitive parameter,it is expected to enhance the predict accuracy,while improve computing efficiency of model calibration,verification,and simplification.In this study,a Vegetation Interface Processes(VIP)model was used to simulate the nitrate content of soil in the North China Plain,in order to evaluate the efficiency of different global sensitivity analysis methods in screening sensitive parameters of agro-ecosystem models.Eight sensitivity analysis methods were selected,including six qualitative and two quantitative sensitivity analysis methods,provided by Problem Solving environment for Uncertainty Analysis and Design Exploration(PSUADE).The results show that:1)The screening information of sensitivity parameter differed greatly in two method of Spearman's correlation coefficient(SPEA)and the Gaussian Process(GP)from others,when verifying the effectiveness of six qualitative sensitivity analysis methods.The SPEA method cannot effectively identify the decomposition rate of microbial pool,whereas,the GP method cannot identify the decomposition rate of structural litter pool.The similar sensitive parameter screening can be found in the Morris method,multiple adaptive regression spline(MARS),Delta Test(DT),and Sum of Tree(SOT)methods;2)Sensitive parameters of soil nitrogen cycle that screened by DT,MARS,Morris,and SOT methods were:the potential nitrification rate,urea hydrolysis rate,Michaelis constant,microbial nitrogen-carbon ratio,slow humus nitrogen-

关 键 词：模型 敏感性分析 土壤硝态氮 PSUADE 抽样方式 

分 类 号：S152[农业科学—土壤学]