PDF《第39 卷第 6 期》

3 个 函数进行优化,且当选择 UCB 作为采集函数时,取κ= 2.5. 2.3 研究对象 本研究以异龙湖为案例地,以IWIND?LR 模型 为待进行参数估值的复杂水质模型.研究中采用的 Sobol 参数敏感性分析运用 Matlab 进行数值计算, 贝叶斯优化基于 Python 软件的 BayesianOptimization 模块.异龙湖(102°30′~102°38′E,23°39′~ 23°42′N) 位于云南省红河哈尼族彝族自治州石屏县境内,湖 泊总面积 28.4 km2 ,最大体积 114.9 万m3 ,平均海拔

1414 m,平均水深 3.9 m,最大水深 5.7 m,全年平均 水温约

20 ℃,由于风的扰动作用与湖水深度较浅, 水体处于完全混合状态. 本研究采用的复杂水质模型为 IWIND?LR,已 被广泛地应用于支撑水质污染防治决策(王冰等, 2016;

张月霞等, 2018). IWIND?LR 模型以国际上 广泛应用的高级水动力水质模型―环境流体力学 代码(EFDC) 作为计算内核,同时针对当前面临的 环境问题,在EFDC 原始版本的基础上开发了一些 高级模块.该模型包含三维水动力、温度动力学和内 部耦合的水质模型模块,能够实现湖泊内部水质过 程的时空精确模拟.该模型可模拟

26 个水质变量,

7 2

0 2 环境科学学报39 卷 模拟底泥与水体之间的营养盐交换及底泥对溶解 氧的动态作用;

基于物质守恒定律,IWIND?LR 模型 能够模拟物理输送、大气交换、吸附解析、藻类吸 收、底泥?水界面交换、硝化和反硝化、沉积成岩等过 程(邹锐等, 2018). 本研 究将异龙湖的地形和边界条件作为IWIND?LR 的输入,模型采用笛卡尔直角坐标网格, 共划分

813 个水平网格,垂向采用 σ 坐标,平均分 为2层.模型计算步长为

15 s,模拟时间为

365 d.模 型所用的输入条件和边界条件以异龙湖

2008 年9月―2009 年8月的........