biowin51的新进展-envirosim.pdf

2016 年 10 月 BioWin 5.1 的新进展 本文将主要介绍新版本在稳态模拟和动态模拟速度上做出的改进，尤其是对于a) 包含很多工艺单 元的复杂工艺流程; b) 包含生物膜系统反应器的工艺流程 （例如MBBR 和IFAS 工艺反应器、浸 没式曝气滤池和滴滤池等等）。BioWin 5.1 中文版会在近期推出！ 模拟速率的改进 – 稳态模拟 稳态模拟的基本数学模型并没有改变，也就是说BioWin 仍利用包含Decoupled Linear Search (DLS) 方法以及 Modified Newton-Raphson (MNR) 方法的综合方法求解。求解程序通 常会先使用DLS 方法进行若干次迭代，然后转为MNR 方法直至求解成功。MNR 方法在求 解过程中效率更高，但是需要进行大量计算。例如在先前的版本中，对于大型复杂的工艺 流程，MNR 方法每次迭代会耗时若干分钟。通常情况下，稳态模拟需要迭代10 次左右找 到稳态解，这样整个求解过程长达20-30 分钟。 BioWin 5.1 版本对 MNR 方法的后台数学模型进行了改进。对于同样的大型复杂工艺流程， 每次迭代只需4 秒左右，而老版本每次迭代则需要200 秒左右，新版本相比原来速度提高 了大约50 倍。 模拟速率的改进 – 动态模拟 BioWin 5.1 采用全新的动态集成法，叫做Backwards Differentiation Formula, 简称BDF 法。 该方法 进行动态模拟的速度要远远快于原来老版本中的方法。EnviroSim 内部BioWin 5.1 测试结果表明： • 对于不包括生物膜反应器的大型复杂工艺流程，动态模拟速度提高 3-6 倍。 • 对于包括生物膜反应器的系统，动态模拟速度提高 10-50 倍。 用户打开数学模型参数界面(Project Current Project Options Numerical parameters), 如 下所示，BDF 方法是默认的动态模拟采用的数学方法，即使在5.1 版本中打开老版本创建 的文件， 这一设置也不会改变。由于新的BDF 法还无法用于SBR 系统的模拟，新版本中 也保留了老版本中的求解方法BW Heun。 1 模型改进 – 生物膜中胶体物质的扩散过程 BioWin 中的生物膜模型包括胶体物质扩散进入生物膜的扩散速率，该参数可以在 Project Parameters Biofilm Effective diffusivities 界面下查看或者修改，如下图所示: 2 在老版本BioWin 中，该参数默认值为 5.0 x 10-12, 大于颗粒态物质的扩算速率 (例如 颗粒态 可降解COD 扩散速率为 5.0 x 10-14), 低于溶解态物质的扩算速率(溶解态可降解COD 扩散速 率为6.9 x 10-10)。该参数直接影响胶体物质在生物膜系统中的分布。以下图所示四个依次 排列 MBBR 反应器系统为例： 若使用默认的扩散参数，溶解态COD 和BOD 浓度在系统中分布如下图所示： 3 需要注意的是溶解态BOD 直到第三个MBBR 反应器才有显著降低。通过观察溶解态COD 分布图可以看出，尽管快速降解COD (rbCOD) 经过前两个MBBR 反应器后就已经基本完全 降解，但是胶体态物质由于扩散速度的限制，还没有被完全降解，所以大量BOD 仍存在 于第二个反应器出水中。 EnviroSim 研究表明类似上述示例中的实际系统中，溶解态BOD 的去除要高于模拟值。因 此在 BioWin 5.1 中胶体物质的扩散速率提高至 5.0 x 10-10, 该值已经接近溶解态物质的扩散 速率 (溶解态可降解COD 扩散速率为 6.9 x 10-10)。这样在经过前两段MBBR 之后，胶体态 物质去除率提高，溶解态BOD 浓度如下图所示： 4 模型改进 – 聚磷菌(PAO)的溶解氧半饱和系数 在老版本中聚磷菌