基于过程模型CROBAS的全局灵敏度分析方法比较

过程机理模型在开发过程中常受限于生理学参数无法直接或准确测量。全局灵敏度分析可以评估模型预测结果对于生理学参数变化的响应,为模型结构改进、数据收集和参数校准提供参考。本研究基于过程模型CROBAS,以华山松为例,选取模型中描述树木结构关系的10个参数,以树高和各器官生物量的Nash-Sutcliffe效率(NSE)为目标函数,比较了3种应用较广泛的全局灵敏度分析方法,即Morris筛选法、基于方差的Sobol指数法和扩展的傅里叶幅度检验法(EFAST)。结果表明: 参数灵敏度排序在不同方法中仅略微有所变化,但对于不同目标函数则区别明显。对算法耗时和收敛效率而言,Morris和EFAST性能较高,Sobol效率相对较低。所有模型输出变量均对单位面积年最大光合速率、比叶面积、消光系数敏感,林冠光截留状态对于林木生长量有着关键性影响,意味着光合固碳量是CROBAS在模型校正和林木生长动态模拟中需要优先进行数据收集、验证与测试的模块。灵敏度分析同时表明,碳平衡理论在林木生物量模拟中最为核心部分是树叶生物量模块的计算与验证。对于复杂过程模型的参数灵敏度分析,如需定性研究可选Morris,而量化评估采用EFAST更适合。