低溶氧下硝化生物膜中菌群的竞争增殖模型

1、低溶氧下硝化生物膜中菌群的竞争增殖模论文作者：王志盈彭党聪袁林江李久义周游 摘要：根据活性生物膜扩张原理，建立了硝化生物膜中微生物相互作用的非稳态解析数学模型，并可预测生物膜厚 度的变化与微生物菌群的空间分布。模型很好地解释了硝化 生物膜中氨氧化菌与亚硝酸氧化菌对空间和氧的竞争，以及 它们在生物膜中的演变关系。关键词：生物膜微生物相互作用竞争模型表面增殖asurfacegrow thcompetitio nmodelofmicr obialspecies innitrifyingbiof ilmatlow dissolvedoxy genactivebiofilabstra ct: based

2、onth eprincipleofmexpansion,a nonsteadystateanalyticmathematicalmodelofmicrobialinteractionsinnitrifyingbiofilmwasdeveloped. themodelcanpredictbiofilm thicknessvariationandspacedistributionofmicrobialspecies itcanalsoexplaincompetitionsforspaceandoxygenbetweenammonia-oxidizersandnitrite-oxidizersand

3、theirevolutioninnit rifyingbiofi lm.tit ionmodel；surfacegrowt h在生物硝化过程中，氨首先由氨氧化菌氧化为亚硝酸盐, 然后由亚硝酸氧化菌氧化为硝酸盐，整个硝化反应是两种菌 群协同作用的结果。因此，菌群之间的生长和基质转化平衡 至关重要，否则将造成中间产物亚硝酸盐的积累。另一方面， 通过对硝化反应的两大菌群进行合理的调控，使硝化过程产 生持久稳定的亚硝酸盐积累并直接进行反硝化，形成短程硝 化一反硝化，则可以达到】1、2：在好氧阶段可减少25% 的需氧量；在缺氧阶段可减少40%的有机碳源；亚硝酸 盐的反硝化速率通常比硝酸盐高倍3 。这对大

4、量的高氨 氮低碳源废水处理具有十分重要的意义。实现亚硝酸盐积累 的方法有游离氨抑制、纯种培养、温度选择等，但它们都只 适用于一定的范围。通过控制反应体系中溶解氧浓度可在生 物膜反应器中实现持久稳定的亚硝酸化，其原因是在硝化生 物膜中存在着氨氧化菌和亚硝酸氧化菌之间对溶解氧的竞 争利用，由此造成生物膜中菌群之间的竞争生长。两种菌群 在硝化生物膜中的空间竞争特点，对反应器的设计以及工程 操作都有很重要的意义。本文目的是在低溶解氧条件下，建立生物流化床内硝化 生物膜中两种菌群的竞争增殖模型，用以解释硝化反应生物 膜中氨氧化菌和亚硝酸氧化菌之间的空间分布和对基质的竞争利用。现有的生物膜模型大多为稳态模

5、型4,认为生物膜的 性质（如生物膜的厚度、活性）、基质和微生物种群在生物膜 中的分布都是恒定的。这些模型及所得结论极大地丰富了生 物膜理论，但实际上即使反应器处在稳态条件下，生物膜的 组成、厚度及种群分布也在不断变化。据此warmer等人建 立了非稳态生物膜模型57,它可以预测生物膜厚度的 变化和描述微生物种群及基质在生物膜中的分布，目前被广 泛采用，但是w anner模型有以下几方面的缺陷：©wanner 模型是平板模型，即认为支持介质可以看作一个平板，而在 生物流化床内采用微粒作为支持介质，其尺度与生物膜为同 一数量级，这容易造成很大的误差；wan ner模型认为整 个生物膜都具

6、有活性，但是随着微型传感技术的发展，证明 在生物膜中02的扩散深度通常为100 um左右。因此，好氧 生物膜中具有活性的部分也为相同的数量级；w anner模 型结构复杂，因素过多，计算过程繁琐，需要专门的软件计 算求解，很难实际应用。1表面活性膜扩张模型的建立基本思路在生物膜体系中，基质、电子受体02通过生物膜表面 的液膜边界层向生物膜内部扩散传递。由于扩散阻力和生物 反应造成基质浓度随生物膜深度变化而变化，与此相应生物 膜中微生物菌群随深度呈非均相分布。当生物膜增加到一定厚度时，在某一深度处，电子受体 的浓度为零。以此为界，外层的微生物可同时得到基质和电 子受体，进行正常的生物化学反应，为

7、活性层；内层的微生 物则由于得不到电子受体，其好氧生物化学反应将终止，为 非活性层（见图1）。一般来说，活性生物层就是距生物膜一水界面100 11 01 左右的一层微生物。生物膜的增长实际上是活性层随时间在 垂直于生物膜一水界面方向上向外扩张，非活性层逐渐变厚 的过程。在活性层中，由于不同菌群对基质的亲和力不同以 及增殖速率的差异，使生物膜中的组分在该层的扩张过程中 不断变化。本模型的主要目的就是描述在硝化生物膜扩张过 程中微生物组分的变化规律。在生物膜反应器中，当溶解氧浓度在2mg/l以下时，硝 化反应的速率受溶解氧水平限制。依此，根据质量守恒原理， 首先建立活性层中溶解氧的扩散一反应方程，

8、求解出生物膜 中溶解氧的分布，然后根据溶解氧分布规律，通过生物膜的 生长方程推导出微生物组分以及生物膜的扩张随时间的变 化情况。基本假定在模型的推导过程中假定： 微生物颗粒为球形，生物膜在同一半径的球面上均匀分布。 忽略外部传质，即不考虑液膜传质阻力。 基质在膜内传质为分子扩散，忽略紊流等其他形式传质。 生物膜的空隙率、密度以及氧的扩散系数8等不 随深度发生变化。 在微生物的增殖过程中不考虑细胞的内源代谢，微生 物对基质的利用符合monod模式。 不考虑微生物膜由于受到流体剪切和生物颗粒之间 的磨擦而发生的脱落。 生物膜中由于硝化反应只发生在距离生物膜表面很 近的薄层中8、9, 一般认为该层厚

9、度约50100 pm,甚至 有的认为仅在152 oum内硝化微生物的活性较强10。 因此假定在距离生物膜一水界面§距离内的生物膜有活性。 在假定中的召薄层中，生物膜的结构、微生物的成 分等保持不变，即活性生物膜的成分只是时间的函数。 温度保持不变，即各种动力学参数为常数。溶解氧扩散一反应模型的建立由于基质扩散所需要的反应时间比微生物生长的特征 时间要小好几个数量级，故相对于微生物增殖来说基质在生 物膜中的分布可认为处于稳定状态。另外，理论和实验均已 证明，当氧氨比11。因此，当液相中溶解氧的浓度1。所以，在生物膜中溶解氧的扩散一反应方程为:式中deff基质在生物膜中的有效扩散系数，m

10、2/sci基质i的浓度，mg/lr生物颗粒半径，mmri基质i的消耗速率,mg/()氧的比消耗速率可以认为是溶解氧浓度的一级反应，即:式中qmax,nh4+,0 2氨氧化反应中氧的最大比反应 速率,mgo2/qmax,nh2-,02亚硝酸氧化反应中氧的最大比反应速率，mgo2 /()ko2,l氨氮化反应中氧气的半饱和常数ko2,2亚硝酸氧化反应中氧气的半饱和常pj生物膜中第j种微生物的体积质量fj生物膜中微生物菌群j所占的体积分率，对于硝化生物膜下标1指氨氮化菌，下标2指亚硝酸氧 化菌£生物膜的空隙率式中p、fj都只是时间的函数，与在薄层中的位置无关。 根据假定qmax,n04+,0

11、2、qmax,n02 -, 02、ko2,l 和 ko2,2 均为常数。因此，式可改为：ro2=kc 02 (3)式中：式可以写为：l微生物膜一一水界面处的半径，mmh生物膜中微生物的种类数8活性层厚度，mmc02, p液相主体溶解氧浓度，mg/l生物膜增殖扩张模型的建立考虑在t时刻活性生物层中第j种微生物的体积分率为 fj,生物颗粒的半径为l,而在t+dt时刻，活性生物层中第 j种微生物的体积分率为fj+df,生物颗粒的半径变为l+dlo 在dt时间内活性生物层中第j种微生物的净增殖为：下面根据质量守恒建立生物膜扩张方程：在dt时间内 活性生物第j种微生物的质量变化二在dt时间内活性生物层

12、中第j种微生物的净增殖：式中yj,o2生物膜中第j种微生物以氧气表示的细 胞产率，mgvss/mg 02式即为生物膜扩张方程，将其简化得：式中ul表面生物膜扩张速率，m/s可以表示为：同时，生物腊怀液相界面处的半径l可以表示为： 式中ro在时间为零时生物颗粒的半径，mm对于硝化生物膜，考虑其生物种群仅为氨氧化菌和亚硝 氧化菌两种，如果某一时刻氨氮化菌在微生物膜中所占的体 积分率fl=f,则亚硝酸氧化菌在微生物膜中所占的体积分率 f2=,将式、联立可以写成以下形式：联立式(5)、(10)、(1 1),根据f和l初始值以及相应 的参数值(表1),可以求出生物膜中不同种类微生物所占比 率和生物膜扩张

13、随时间的变化关系。表1计算过程中参数值2、12参数取值单位qma x,nh4+,02mg0 2/()qmax,nh2 -,02mg02/()k o2,2mg02ko2, lmg02yl,o2mg vss/mg02y2,o 2mgvss/mgo2d eff, o2x10-9m 2/s £ p 256g/lr 0mmco2, pmg/l2结果与讨论l与f随时间的变化通过基质降解和生物膜生长两过程的计算，可以得出在 一定条件下生物膜结构和生物颗粒半径随时间的变化关系。 图2给出液相主体溶解氧浓度cp=mg/l,活性生物膜厚度6 =50 pm时，氨氧化菌在生物膜中的比率f和生物颗粒半径l 随

14、时间的变化。从图2可以看出，氨氧化菌在生物膜中的比率随时间的 延长逐渐增大，也就是说随着生物膜的扩张氨氧化菌将占据 硝化生物膜表面的位置，使氨氧化过程的反应速率增大；由 于亚硝酸氧化菌的比率减少，从而使亚硝酸氧化菌逐渐失去 在生物膜中的有利位置，亚硝酸的氧化速率减小，出水中的 n02n浓度逐渐增加，直至完全亚硝酸化，这与大多数的试 验观察完全符合。cp对活性生物膜中菌群分布的影响由式（5）、（10 ）和（11）可以看出，影响生物膜菌群分布 变化的因素有活性生物膜的厚度§和液相主体的氧浓度cpo 图3为cp对氨氧化菌在活性生物膜中的比率f的影响。从图3可以看出，当活性生物膜的厚度保持不

15、变时，液 相主体浓度cp越大越有利于氨氧化菌的竞争生长，氨氧化 菌在较短的时间内就可以达到较高的比率，也就是可以在较 短的时间内完成氨氧化菌与亚硝酸氧化菌之间的选择竞争, 从而可在较短的时间内实现出水亚硝酸盐的高积累。液相主 体的溶解氧浓度越高，通过分子扩散进入生物膜内部特定位 置上的溶解氧浓度也就越高，由于氨氧化菌对溶解氧的亲和 力比亚硝酸氧化菌对溶解氧的亲和力要强，导致的氨氧化菌 生长速率的增加比亚硝酸氧化菌生长速率的增加要大，使竞 争过程更为显著。5对活性生物膜中菌群分布的影响从图4可以看出，当液相主体浓度维持一定时，活性生 物膜厚度§越薄，越容易实现氨氧化菌和亚硝酸氧化菌的选

16、 择竞争，氨氧化菌就可以在较短的时间内占据生物膜表面的 位置，完成对氨氧化菌的选择。这也可以用氨氧化菌和亚硝 酸氧化菌对溶解氧的亲和力不同来解释，在一定的溶解氧浓 度条件下，当活性生物膜较薄时，其中的溶解氧平均浓度较 大，也就容易实现氨氧化菌的竞争生长。3结语通过建立基质扩散一反应及微生物生长的数学模型，从 理论上证实了在硝化生物膜中存在着氨氧化菌与亚硝酸氧 化菌之间对空间和基质的竞争。这种竞争是由于微生物本身 的增殖速率差异造成的。竞争的必然结果是随着时间的延续, 氨氧化菌在活性生物膜中的比率逐渐增大，直至占据整个活 性层，从而实现持久稳定的亚硝酸化。参考文献：1 garridojmet o

17、fdissolvedo xygenconcentrationonnitr iteaccumulat ioninabiofil mairliftsusp ensionreacto r j biotech &bioeng, 1997 , 53： 168-178 2 piciorean ucettheeffec tofoxygencon centrationon nitriteaccum ulationinabi ofilmairlift suspensionre actor j . wat scitech, 1997 , 36： 147-1563 ablingu,s eyfriedctre

18、a tmentofhighs trengthammon iumwastewate rnitrogenrem ovalusingnit rite j . watscitech, 1992,26： 1007-1015.4王世和等生物流化床降解有机物的动力学机理研 究j东南大学学报,1993, 23(1 ) : 62-68. 5 f ruhenmetofsp atialdistrib utionofmicro bialspeciesi nmixedcultur ebiofilms jwatscitech,1990,23：13656 re ichertp,wann erofsolidsin biofilms：sig nif icanceof 1 i quidphase tr anspor t j .wat scitech, 19 97,36:321-32 8.7 wanner , reichertmod elingofmixed -culturebiof ilms j . biot ech&bioeng,