膜生物反应器“两次分离”数学模型

摘要：膜生物反应的实质是生物降解与膜分离相互影响、共同作用的过程。膜对活性污泥混合液的分离过程是由“两次分离”实现的：一次分离是混合液在压力作用下进行固液分离，二次分离是分离出来的液相透过沉积层到达膜表面，实现有机物大分子及胶体物质与水的分离。固液分离过程可用阻力模型来描述；而膜表面有机物大分子与水的分离可用凝胶极化模型描述。稳态时，固液分离的速率与有机物大分子和水分离的速率相等。

关键词：膜生物反应器 两次分离 数学模型

　　膜生物反应器是由超(微)滤膜组件和生物反应器构成的。传统超(微)滤过程的进料一般为单相流体(即只有液相)，而膜生物反应器中超(微)滤分离的活性污泥混合液却是固—液两相流体，既包括活性污泥悬浮固体，也包括液体以及溶解在其中的有机物质，膜不但要截留大分子有机物，而且要分离固体颗粒。因此，固体颗粒在错流中的迁移规律(包括固体颗粒在剪切流中的运动规律以及固体颗粒在膜面的沉积规律)就必然会对整个膜分离过程产生影响，对此过程进行分析并建立新的数学模型加以描述是十分有意义的。

　　1 固体颗粒在错流过滤中的迁移规律

　　在错流过滤中，固体颗粒的运动受沿膜面平行流动的剪切流和垂直于膜面的过滤渗透流的共同作用。渗透流趋向于将固体颗粒曳向膜面，而剪切流力图保持颗粒悬浮，将其随循环液带出膜面。固体颗粒是沉积于膜面还是随循环液流出，取决于剪切流和渗透流曳力的相互作用结果。
　　固体颗粒迁移运动趋于膜面的速率和流出膜组件的横向迁移速率应该和颗粒直径、流道几何尺寸、剪切流速率和膜渗透过滤速率有关。
　　借用Alter和Belfort固体颗粒在错流微滤过程中的迁移机理^［1］，对于达到稳定状态时，即膜面的渗透过滤速率保持稳定时，颗粒迁移运动趋于膜面的速率为：

　　V_P=V_m［g(β)+R_ewh₂(β)］ (1)

　　式中　VP——颗粒迁移运动趋于膜面的速率
　　　　　Vm——膜渗透过滤速率
　　　　　Rew——渗透系数
　　　　　β——流道无量纲宽度
　　　　　g(β)、h2(β)——β的函数
　　颗粒横向迁移速率为：

　　V_e=R_ep(d_p/2B)²u_mf(β ) (2)

　　式中　Ve ——颗粒横向迁移速率
　　　　　um——流道剪切流速率
　　　　　dp——颗粒直径
　　　　　B——流道宽度
　　　　　Rep——系数
　　　　　f(β)——β的函数
　　由式(1)、(2)可见，颗粒迁移运动趋于膜面的速率与膜渗透过滤速率成正比；颗粒横向迁移 速率与流道剪切流速率成正比。

　　2 固体颗粒沉积阻力模型的建立

　　错流过滤中，固体颗粒在过滤渗透流和剪切流的综合作用下，有向膜面沉积的趋势。随着过滤的进行，会在膜面逐步形成沉积层。沉积层累积阻力Rd可由下式确定：

　　Rd=αd·Md (3)

　　式中　Rd ——固体颗料沉积层累积阻力
　　　　　αd——比阻
　　　　　Md——沉积负荷(单位膜面积的沉积量)
　　比阻αd可由Kozeny—Carman关系计算［2］：

　　αd＝180(1-ε)²/（d_s²ε³） ( 4)

　　式中　ε ——沉积层孔隙率
　　　　　ds——固体颗粒粒径
　　由式(4)可见，沉积层愈密实，固体颗粒粒径愈小，沉积层比阻愈大。
　　比阻受操作压力P的影响，并有：
　　α_d＝αd₀·Pⁿ (5)

　　式中　αd0——系数
　　　　　n——压缩系数
　　由超滤过程的阻力模型［3］得膜通量为：

　　J_v=P/μ(R_m+R_d)=P/μ·R _总 (6)

　　式中　Rm ——膜本身的阻力
　　　　　Rd——累积阻力
　　　　　R总——总阻力
　　　　　P——膜组件的操作压力
　　　　　μ——动力粘滞系数
　　将式(3)、(5)代入式(6)中，有：

　　Jv=1/μ（P/（R_m+αd₀M _dP_n）） (7)

　　当操作压力较大时，Rm＜αd0MdPn，膜通量可近等于：

　　Jv=1/μ（P^1-n/αd₀ M_d ） (8)

　　由于在操作压力较大时已形成凝胶层，在此情况下，Jv已基本与P无关，亦即相当于n＝1，故式(7)可简化为：

　　Jv=1/μ[P/（R_m+αd₀M_dP）] (9)