3. 生物膜过滤的传质与反应基础

生物膜过滤之所以能同时完成“截留”和“降解”,关键在于污染物必须先穿过液膜、进入生物膜内部,才能被微生物酶系捕获并转化。也就是说,反应速率并不只由微生物活性决定,还受“底物能不能及时送到”的限制。本节围绕三个问题展开:底物如何扩散到生物膜、膜内传质阻力如何量化、以及用什么动力学方程把传质与反应耦合起来做工程估算。

一、从主体液相到生物膜:两级传质串联

以一块附着在滤料或填料表面的好氧生物膜为例,底物(以溶解性有机物 BOD 或氨氮计)从主体水流迁移到膜内反应位点,要经过两段阻力:第一段是膜外侧的滞流液膜(厚度通常几十到数百微米,随流速增大而变薄),第二段是生物膜本体(厚度常在 0.1~2 mm 量级)。两段串联,总通量可写成:

J = k_L · (S_b − S_s) = D_e · (dS/dz)|膜表面
其中:J   底物通量,g/(m²·d)
      k_L 液膜传质系数,m/d
      S_b 主体液相底物浓度,g/m³
      S_s 生物膜表面底物浓度,g/m³
      D_e 生物膜内有效扩散系数,m²/d

注意符号约定:本文统一取“由液相指向膜内”为正方向,z 坐标自膜表面向膜底(载体侧)增大。液膜阻力大时,S_s 明显低于 S_b,此时即使膜内酶活性很高,整体去除速率也被“供料”卡住;反之,流速高、液膜薄,S_s≈S_b,过程转入膜内传质或反应控制。工程上判断属于哪种控制,是选型和放大前必须做的第一步。

二、生物膜内的有效扩散系数与传质阻力

生物膜不是清水,而是由胞外聚合物、菌体和高含水孔隙构成的非均相介质。底物在其中的扩散路径曲折、部分孔隙被堵塞,因此有效扩散系数 D_e 通常只有同温度清水中分子扩散系数 D_w 的 0.2~0.8 倍。对 BOD、氨氮这类小分子,可取 D_e ≈ 0.5 D_w 作为初步估算;对难降解大分子或胶体态底物,D_e 会更低,需实测或按经验折减。

膜内传质阻力的大小,用无量纲的蒂勒模数(Thiele modulus)φ 表征,它本质上是“反应速率”与“扩散速率”的比值。对平板状生物膜、一级反应简化情形:

φ = L · sqrt(k₁ / D_e)
其中:L  生物膜厚度,m
      k₁ 一级反应速率常数,1/d
      D_e 有效扩散系数,m²/d
当 φ ≪ 1:反应控制,膜内浓度近似均匀;
当 φ ≫ 1:扩散控制,底物在膜表层即被消耗,深层“饥饿”。

φ 的工程意义很直接:它告诉你生物膜“有效厚度”有多大。扩散控制时,盲目加厚生物膜不会提高去除率,反而增加曝气负担和脱落风险;此时应减薄膜、提高流速或降低有机负荷。下面这张图把液膜、生物膜、载体三层的浓度分布与通量关系画在一起,便于对照理解。

生物膜过滤中的底物传质与浓度分布 主体液相(流动水) S_b 较高,湍流混合 滞流液膜 k_L 控制 生物膜(厚度 L) D_e 有效扩散 酶促反应消耗 载体/滤料 无传质 S_b S_s S≈0 底物通量 J 方向(液相 → 膜内) 控制步骤判断 液膜控制: S_s ≪ S_b,提高流速有效 膜内扩散控制: φ 大,减薄膜厚有效 反应控制: φ 小,提高活性/负荷 工程先判控制步骤, 再定膜厚与流速。

三、反应动力学基本方程与简化假设

生物膜内的底物消耗通常用莫诺方程描述,并与扩散项组成反应—扩散方程。对平板膜、稳态、忽略对流的一维情形:

D_e · d²S/dz² = (k · X / Y) · S / (K_s + S)
边界条件:z = 0 处  D_e·dS/dz = k_L·(S_b − S_s)
          z = L 处  dS/dz = 0(载体不透水)
其中:k 比降解速率,1/d;X 生物膜密度,g/m³;
      Y 产率系数;K_s 半饱和常数,g/m³

该方程一般无解析解,工程计算依赖三条常用简化假设,选型时必须清楚自己用了哪一条:

  • 零级简化:当 S ≫ K_s,反应速率近似恒定,膜内浓度呈抛物线分布。适用于高浓度进水段,如初沉后 BOD 较高的过滤段。
  • 一级简化:当 S ≪ K_s,速率与浓度成正比,方程化为线性,可得解析解,用 φ 判断扩散影响。适用于低浓度深度处理段,如二级出水过滤。
  • 薄膜/薄层假设:当 φ 很小或膜很薄时,可忽略膜内浓度梯度,把整膜视为一个均匀反应器,直接用表面浓度 S_s 计算去除速率。

四、一个可复算的估算示例

设某深度处理生物膜滤池,进水氨氮 S_b = 8 g/m³,水温 20 ℃,氨氮分子扩散系数 D_w ≈ 1.5×10⁻⁴ m²/d,取 D_e = 0.5 D_w = 7.5×10⁻⁵ m²/d;生物膜厚 L = 1.0 mm = 1×10⁻³ m;一级速率常数 k₁ = 30 /d。先算 φ:

φ = L·sqrt(k₁/D_e) = 1×10⁻³ × sqrt(30 / 7.5×10⁻⁵)
  = 1×10⁻³ × sqrt(4×10⁵) ≈ 1×10⁻³ × 632 ≈ 0.63

φ≈0.63 属中等偏反应控制,膜内浓度下降有限,说明该膜厚下扩散尚未成为绝对瓶颈,但已不可忽略。若把膜厚增至 2 mm,φ 翻倍到约 1.26,进入明显扩散控制区,此时有效去除率不会随膜厚线性增长。工程判断依据是:在保持其他条件不变的前提下,实测去除率随膜厚增加而趋于平台,即说明已进入扩散控制,应停止加厚、转而优化水力或减薄膜层。上述数值为教学示例假设,实际设计须以现场小试或中试数据校核 D_e、k₁ 与 k_L。

选型提示:液膜传质系数 k_L 随滤速提高而增大,但过高滤速会带来水头损失和生物膜脱落。建议在 φ 计算后,用“提高滤速看去除率是否上升”做一次现场判别:上升明显说明液膜控制,应优先提流速;几乎不变则说明已进入膜内扩散或反应控制,应调整膜厚或负荷,而不是继续加大流速。