多组分给水管网水质模型的开发与应用.docxVIP

多组分给水管网水质模型的开发与应用 数学模型是研究给排水网络的内部动态过程、给排水网络的规划、设计、评估和优化供水网络的运营，以及供水网络设施和水质安全的重要方法。由于研究人员在网格内部过程中对检测到的单元的理解和计算操作能力的局限性，原始给排水网络水质模型主要为分组模型。换句话说，认为网格中某些组件的反应速率与自身浓度有关。这种模型即使应用于大型管网系统，也可以快速求解。根据大量研究，这些模型可以模拟实际管网的消毒、跟踪、污染等方面的变化。然而，在实际的给水管网络中，许多组件之间的相互作用往往不容忽视。客观上，研究人员需要进行多组分散水质模拟。 EPANET是美国环境保护署(environmental protection agency, EPA)开发的一套源代码公开的给水管网水力学和水质模拟软件,其现行版本为2000年发布的2.0版.EPANET 2.0的水质模拟模块是单组分模型,为了增强其水质模拟功能,美国EPA于2007年8月发布了支持多组分水质模拟的扩展版,即EPANET-MSX(Multi-Species eXtension).用户通过联合使用EPANET-MSX和EPANET 2.0的工具包可以自定义给水管网的反应动力学方程,从而模拟管网中更多的水质组分和反应过程. 本研究基于EPANET-MSX工具包开发机制性的多组分给水管网水质模型,利用管段模拟实验数据进行模型验证,并在案例研究中应用该模型模拟给水管网的水质动态变化及其超标风险. 1 学习方法 1.1 分离产物的代谢特性表7 模型将给水管网的管段视为推流式反应器,并且忽略轴向的扩散过程,则对于如图1所示的任意Δx的微元可以建立如式(1)所示的物料平衡方程. ?[ci]?t=?U?[ci]?x+r(ci)(1)?[ci]?t=-U?[ci]?x+r(ci)(1) 式中,[ci]为微元中模拟变量ci的浓度(mg/L),U为管段流速(m/min),r(ci)为变量ci在微元中的总体反应速率[mg·(L·min)-1],t为时间(min),i为模拟变量的编号. 每个微元包括管壁、生物膜、液膜和液相4个组成部分,其中发生的物质传递和转化过程如图1所示.模型的模拟变量共15个,包括生物膜和液相中各7个变量,即溶解性有机物、不溶性有机物、氨氮、余氯、异养菌、自养菌和惰性颗粒,以及管壁生物膜厚度.为简化结构,模型做出以下基本假设: (1)微元内生物膜和液相均为完全混合式反应器,二者通过液膜发生传质; (2)生物膜均匀分布在管壁且性质均一,传质仅发生在与管壁垂直的方向,而不考虑与管壁平行方向的传质; (3)液相中溶解性物质通过分子扩散穿过液膜进入生物膜,其扩散速率符合Fick定律,且溶解性物质的分子扩散系数均相等; (4)液相中不溶性物质通过吸附过程进入生物膜,其吸附速率符合一级反应,且不溶性物质的吸附速率常数均相等; (5)生物膜中的不溶性物质脱落后进入液相,其脱落速率符合一级反应,且不溶性物质的脱落速率常数均相等; (6)液相和生物膜中的异养菌或自养菌的代谢特性相同,其生长速率均符合Monod方程,内源呼吸引起的自然衰减速率符合一级反应,而余氯灭活导致的死亡速率符合二级反应,微生物衰减或被余氯灭活之后均转化为惰性颗粒; (7)管网中余氯继续与有机物发生氧化和卤代反应,其反应速率符合二级反应; (8)液膜传质速率和生物膜中颗粒脱落速率均与管段流速成正比; (9)管段进出水中不溶性组分的质量浓度与浊度线性相关. 基于上述假设可以推导得到各模拟变量的反应速率r(ci)的表达式,进而建立物料平衡方程.限于篇幅,这里仅给出液相中7个模拟变量的物料平衡方程,即式(2)～(8),式中变量和参数的含义如表1所示.模型共有20个待率定参数,其它参数及方程参见文献. d[SCODL]dt=?μHYH[SCODL]kHC+[SCODL][NH3L]kHN+[NH3L][HBL]?DFAVU([SCODL]?[SCODF])?kchr[SCODL][Cl2L]71×1?000?eC/Cl(2)d[NCODL]dt=kdetAVU?LF[NCODL]?kadpAV[NCODL]?kchr[NCODL][Cl2L]71×1?000?eC/Cl(3)d[NH3L]dt=?μHYH[SCODL]kHC+[SCODL][NH3L]kHN+[NH3L][HBL]?eN/C?μAYA[NH3L]kAN+[NH3L][ABL]?DFAVU([NH3L]?[NH3F])(4)d[HBL]dt=μH[SCODL]kHC+[SCODL][NH3L]kHN+[NH3L][HBL]?kdH[HBL]?kadpAV[HBL]?kint[HBL][Cl2L]+kdetAVU?LF[HBF](5)d[ABL]dt=μA[NH