第一作者：黄圣散，男,1981年生，硕士研究生，研究方向为水污染控制工程。
膜生物反应器次临界通量运行的膜污染特性研究
黄圣散 吴志超
(同济大学污染控制与资源化研究国家重点实验室，上海200092)
摘要 膜生物反应器(M BR)是将膜分离与生物反应相结合的污水处理新工艺，近年来已引起广泛的关注，但不可避免的膜污
染限制其更广泛的应用。临界通量在膜污染控制中是个非常重要的概念。本试验研究平板膜生物反应器在次临界通量运行下的膜污染状况，并结合膜污染模型进一步表征膜表面的污染特性。试验结果表明，该平板膜生物反应器在次临界通量运行的情况下，膜污染可分为膜污染缓慢发展阶段(第Ⅰ阶段)和膜污染迅速发展阶段(第Ⅱ阶段)，可分别用膜孔堵塞模型和泥饼阻力模型表征膜阻力与时间的变化关系。同时，对运行后的膜阻力分布进行分析，表明泥饼阻力和孔道吸附堵塞阻力是膜污染的主要组成部分，分别占到总阻力的73%和24%，而膜本身阻力仅占3%。
关键词膜生物反应器临界通量膜污染
Study on membrane bioreactor fouling in sub -critical flux operation H uang Shengs an,W u Zhichao.(State K ey L abo-r ator y of Pol lution Contr ol and Resour ce Reu ,T ongj i Univer sity ,Shanghai 200092)
Abstract: M embrane bior eacto rs (M BRs),in w hich membranes ar e applied to biolog ical wastew ater tr eat ment,hav e been gr eat ly dev elo ped r ecently,but membrane fouling in M BRs r est ricts their widespread u.Critical flux is an import ant parameter to contr ol membrane fouling.T his st udy investig ates the chang e in membrane fo uling in a flat membr ane bio reactor operat ing at a sub -crit ical flux.T he prev ious theo retica l mo dels a re ud to simulate the mem -br ane fouling characteristics.Ex per imenta l results show that membrane fouling can be divided into tw o stages:the first period during w hich fo uling g radually occur s and the cond period in w hich fouling ris rapidly.T he po re bloc -king resistance and cake resistance models can be ud to simulate resistance w ith time in differ ent perio ds.T he analy -sis o f membr ane r esist ance distributio n show s that resistance ar ising from cake fo rmatio n is 73%o f to tal resistance,resistance due t o pore blocking is 24%,and clean membr ane r esistance is only 3%.
Keywords: M embr ane bio reactor (M BR) Cr itical flux M embr ane fo uling
膜生物反应器(M BR)是将膜分离与生物反应相结合的污水处理新工艺，近年来已引起广泛的关注，并逐步应用于城市污水和工业废水的处理[1~3]。由于传统生物处理技术二沉池固液分离性能的限制，使得反应器的微生物浓度不能太高，而M BR 能成功代替二沉池进行高效固液分离，大大提高反应器的污泥浓度，同时具有出水水质好、耐冲击负荷强和污泥产率低等优点。尽管MBR 有许多其他污水处理技术所不具备的优点，但不可避免的膜污染限制其更广泛的应用[4,5]。教育价值
临界通量的概念最早由Field 在1995年提出，此后国内外很多学者对此开展广泛研究[6~8]。MBR 临界通量即保持一个恒定的膜通量低于此值运行时，污泥颗粒不能在膜表面聚集沉积，操作压力保持稳定，没有显著的膜污染产生；而当控制膜通量高于此值运行时，操作压力急剧升高，膜污染随之迅速上升，导致MBR 系统的崩溃。所谓次临界通量
运行，就是在一定的操作条件下，使膜通量维持在临界通量以下以获得M BR 的稳定运行。本文针对膜生物反应器在次临界通量运行下的膜污染状况进行试验研究，并结合膜污染模型进一步表征膜表面的污染特性。1
试验装置和方法
1.1试验工艺流程
试验用水为某城市污水处理厂沉砂池出水，为强化脱氮效果，采用缺氧池-好氧膜生物反应器工艺，工艺流程如图1所示。贮液罐内的污水经过进水泵提升进入反应器，通过进水槽的溢流和止回阀保证反应器内的水位恒定，溢流废水回流至贮液罐。进水经缺氧池重力自流进入好氧膜生物反应器内，通过出水泵对膜组件抽真空运行，真空表显示膜过滤操作的压力，同时由回流泵进行内循环。 缺氧池和好氧膜生物反应器的有效容积分别为
#
512# 环境污染与防治 第27卷 第7期 2005年10月
6、12L，试验中控制内回流比为2，反应器内污泥质量浓度7~8g /L，除正常取样分析外，不排泥。试验采用平板膜式组件，膜材质为聚偏氟乙烯(PVDF)，上海应用物理研究所提供，膜的有效过滤面积为0.1m 2
，截留分子量100ku 。次临界通量运行时采取间歇抽停出水的方式(8m in 开,2m in 停)，并控制曝气量为0.6m 3/h
。
图1 试验工艺流程
1-贮液罐;2-缺氧池;3-好氧M BR;4-进水泵;5-止回阀;6-回流泵;7-搅拌器;8-真空表;9-出水泵;10-膜;11-空压机
1.2
临界通量的测定方法
临界通量采用通量阶梯式递增法进行测量[7,8]
。即在一定的操作条件下，控制MBR 在恒定的低通量下连续运行，观测操作压力$P 在一定时间段$T (至少15min)内的变化，若$P 保持稳定，可认为此时的通量低于临界通量，再使膜通量增加一个阶量，重复上述试验。通过逐步提高膜通量，当MBR 刚好高于某个通量运行时,$P 急剧变化，不能稳定运行，可认为此时的通量已高于临界通量。1.3离别祝福的诗句
膜污染表征
膜通量和操作压力之间的关系可用式(1)表示[9,10]:
J =$P L R =
$P L (R m +R c +R p )
(1)
式中:J 为膜通量,L/(m 2
#h);$P 为膜操作压力,Pa;L 为渗滤液动力粘度系数,Pa #s;R 为过滤总阻力,m -1;R m 为膜本身的阻力,m -1;R c 为膜表面泥饼层所产生的阻力,m -1;R p 为膜孔堵塞与吸附所产生的阻力,m -1。
膜过滤操作中的膜孔堵塞模型和泥饼阻力模型的膜通量与时间的变化关系[11]: 膜孔堵塞模型:J =J 0ex p (-K p t )
(2) 泥饼阻力模型:J 2
=
J 21+J 20K c t
(3)
式中:J 0为初始时刻的膜通量,L/(m 2#h);K p ,K c 分别为膜孔堵塞模型和泥饼阻力模型的常数值。将
式(1)分别代入式(2)、(3)变形后，得到膜孔堵塞模型和泥饼阻力模型的膜阻力与时间的变化关系: 膜孔堵塞模型:ln R p =ln R 0+K p t
(4)
泥饼阻力模型:R 2c =R 2
期中考试总结发言稿>海螵蛸的功效与作用
0+K t
(5)
式中:K =K c @$P 2L
2
;R 0为初始时刻的膜阻力,m -1。 根据上述公式可以对不同形式的膜污染进行评价和表征。2 结果和讨论2.1 临界通量的测定
在曝气量为0.6m 3/h，污泥质量浓度为8g /L 时，膜通量为4L/(m 2
#h)，以后每20m in 将通量提高3L/(m 2
#h)，直至出现操作压力急剧升高，膜污染迅速加剧为止。临界通量测定结果如图2所示，可知当运行通量不高于16L/(m 2#h)时，在20
m in 的连续抽吸间隔内操作压力能稳定为一个数值；而当通量高于19L/(m 2#h)时，操作压力随抽吸
时间变化非常显著，尤其控制通量为28L/(m 2#h)时，操作压力急剧升高，导致系统迅速崩溃。因此，可认为该平板膜在上述操作条件下的临界通量为16~19L/(m 2#h)
。
图2 临界通量的测定
2.2 次临界通量运行情况
根据临界通量的测量结果，在同样的操作条件下，当膜通量低于16L/(m 2#h)时，可视为次临界通量运行。本试验控制通量为13L/(m 2
#h)，考察M BR 在次临界通量运行时的压力和阻力变化情况，结果如图3
所示。
图3 次临界通量运行时压力和阻力随时间变化
由图3可知，膜污染可分为2个阶段:
(1)膜污染缓慢发展阶段(第Ⅰ阶段)。由于
#
513#黄圣散等 膜生物反应器次临界通量运行的膜污染特性研究
MBR 以次临界通量运行时，在开始阶段操作压力较低，溶液中的大分子物质和颗粒物质在低压力的驱动作用下缓慢在膜表面吸附、沉积，逐渐造成膜孔径变小或堵塞。该阶段大致为MBR 运行的前336h，膜阻力以平均4.55@109
m -1
/h 的速度增长，导致操作压力也随之逐渐变大，但压力的增长率($P /$t )非常小，为0.02kPa/h 。
(2)膜污染迅速发展阶段(第Ⅱ阶段)。随着前一阶段膜孔的不断吸附和堵塞，膜对混合液的过滤性能也不断发生改变，导致膜的新临界通量比投入长期运行前膜的临界通量小[12]。当污染积累到一定程度后，
原先加在膜上的通量超过此时膜实际的临界通量，因而膜污染加速，操作压力迅速变大，直至污泥絮体在膜表面不断沉积形成一层致密的泥饼层，此刻即使控制再大的压力也不能维持膜的恒通量运行。该阶段大致持续了144h，膜阻力以平均3.56@1010m -1/h 的速度增长，增长率大约是第Ⅰ阶段的8倍，同时操作压力的增长率($P /$t )也很大，为0.16kPa/h 。
根据膜污染的阻力表征模型(4)和(5)，分别以ln R -t 和R 2-t 对第Ⅰ阶段和第Ⅱ阶段作图，结果如图4和图5
所示。
图4 第Ⅰ阶段ln R 与时间的变化关系
图5 第Ⅱ阶段R 2与时间的变化关系
由图4、图5可知，第Ⅰ阶段较好地符合膜孔堵塞模型，膜污染的产生主要是由于溶液中的大分子
物质和颗粒物质吸附到膜表面，进而堵塞膜孔，发展比较缓慢；第Ⅱ阶段较好地符合泥饼阻力模型，污泥絮体在膜表面的不断沉积最终导致污泥层的形成，膜污染急剧发展。
2.3 膜阻力分布情况
各部分阻力的测定方法:R m 为膜投入运行前在清水中测得;R p 为膜以次临界通量运行时第Ⅰ阶段的总阻力减去R m ;R c 为第Ⅱ阶段的总阻力减去第Ⅰ阶段的总阻力。测定结果如图6所示，可知泥饼阻力和孔道吸附堵塞阻力是膜污染的主要组成部分，分别占到总阻力的73%和24%，而膜本身阻力仅占3%,
相比可忽略不计。
图6 膜阻力分布比例
3结论
(1)该聚偏氟乙烯(PVDF)平板膜的操作条件为曝气量0.6m 3/h，污泥质量浓度8g /L 时，临界通量为16~19L/(m 2
#h)。
(2)平板膜生物反应器在次临界通量运行的情况下，膜污染可分为2个阶段：膜污染缓慢发展阶段(第Ⅰ阶段)和膜污染迅速发展阶段(第Ⅱ阶段)。 (3)第Ⅰ阶段较好地符合膜孔堵塞模型，污染阻力模型方程为ln R =26.47535+0.0045t ；第Ⅱ阶段较好地符合泥饼阻力模型，污染阻力模型方程为R 2
=-1.304@1026
+3.655@1023
t 。
(4)泥饼阻力和孔道吸附堵塞阻力是膜污染的主要组成部分，分别占到总阻力的73%和24%，相
比而言，膜本身阻力仅占3%可忽略不计。参考文献
1 T an l B,Bao W Y,T an l I N.Characteriz ation of fouling kinet -ics in ultrafiltration s ystem s by resistances in ries model.De -salination,2000,129:7~14
2 Tatsuki U,Ken ji H.Domestic w astew ater treatment by a sub -merged membrane bior eactor w ith gravitational filtration.Wat.
Res.,1999,33:2888~2892
3 Silva D G V,Urbain V,Abeysinghe D H ,et al.Advanced an aly -sis of m embrane -bior eactor performan ce w ith aerobic -anoxic cyc -ling.W at.Sci.T ech.,1998,38:505~512
4 Defrance L,Jaffrin M Y,Bh arat G,et al.Contrib ution of various
cons titu ents of activated sludge to membrane bioreactor foulin g.Bioresource T echnology,2000,73:105~112
5 W ontae L,Seoktae K,H angsik S.Sludge characteris tics and their
田园风光contrib ution to microfiltration in submerged membrane bioreac -tor.Journal of M em brane Science,2003,216:217~227
(下转第522页)
#
514#
日上午10B30，西北风向，平均风速取8m/s(风力4~5级)，总云量为5，低云量取1，大气稳定度为D 类
条件下分别对闵行区的吴泾化工、宝山区的吴淞工业区、杨浦老工业集中地区、普陀区的桃浦工业区所在地，污染源叠加情况下计算所得污染物的浓度贡献，得到的地面污染浓度分布图示例，计算网络大小为1km@1km。
表31998年上海市各区县工业污染源情况
行政区面积/km2*煤燃烧量/(t#a-1)工业点源个数/个黄浦125893839
卢湾81263611
徐汇5532242257
长宁38192819.846
静安83604818
普陀55308806.974
闸北2910063246
虹口2363314.432
杨浦611176814105
闵行3728064713126
宝山41513287708122
嘉定459446851.9124
浦东新区5234884530188金山586264661260
松江605132433.455
青浦676148256.150
南汇688195188.776
分发用英语怎么说奉贤687339145.873
崇明1041697499.448
全市634133645814.41320
注:*资料来自中国行政网站
6结论
通常，大气扩散模型要求输入多项参数，对每一像元手动输入参数是件非常繁琐的事，本文将高斯扩散与相关数据整合，模拟工业点源污染对区域大气质量的影响，有效地解决了模型参数输入、模拟结果可视化等问题，实现了模式与GIS功能的完全集合，真正实现了人机交互。结果显示，高架源排放的大气污染物在一定的气象条件下完全可以影响到城市的空气质量。从全市1300多个工业点源的模拟扩散结果可以分析上海市SO2污染的地域分布特征：污染浓度较高的区域由分散的若干个高浓度中心组成，这些中心分别位于宝山区的吴淞工业区、杨浦老工业集中地区、普陀区的桃浦工业区、闵行区的吴泾化工业区等。这是由于这些工业区多大型钢铁厂、发电厂、水泥厂，污染排放量大、排放高度高。
本文只是一次初步尝试，该模型有待进一步改进，在今后的工作中笔者将会扩充子模块(如排放模型、线型/面型污染扩散模型和风险评价模型等等)，最终开发一个城市大气质量管理系统。
参考文献
1陈国海.上海环境空气污染类型的定位.上海环境科学,2003, 22(4):230~233
2Lu R,Tu rco R P。An integrated air p ollution modeling system for u rban and regional scales, 1.Structure and performance.
Journ al of Geoph ysics Res earch1997,102D:6063~6079
3Ku mar A,Bellm an N K,Sud A.Perform ance of an indu strial sour ce complex model:predicting long-term concentration s in an urban ar ea.En vir on mental Progr ess,1999,18(2):93~100
4Sivacoumar R,Bhanarkar A D.Air p ollution modeling for an in-dustrial complex and model perform an ce evaluation.Envir on-men tal Pollution,2001,111:471~477
5刘文野，马建武.安宁市工业大气污染物迁移及分布规律研究.
西南林学院学报,2002,22(2):29~31
6M ehdizadeh F,Rifai H S.M odeling point s ou rce plum es at high altitu des using a modified Gauss ian model.Atm os pheric Envir on-men t,2004,38:821~831
7赵碧云，贺彬，朱发庆，等.大气污染扩散空间信息系统.环境科学研究,1999,12(6):10~12
耳洞发炎8Dai J M,Rocks D M.A GIS-bad app roach to s patial allocation of area s ou rce solvent emission s.En vironmen tal M odelin g& Software,2000,15:293~302
9孙晓梅.工业污染源地理信息系统研究.环境保护科学,2003, 28(111):42~48
10王远飞，张超，何洪林.GIS支持下的上海市工业源SO2污染研究.上海环境科学,2003,22(11):801~804
11牛文胜，孙振海.大气扩散模式的简要回顾.气象科技,2000,28
(2):1~4
12张宝宁，程水源，任毅斌，等.城市工业区大气污染物浓度场确定方法的研究.化工环保,1997,17(2):114~119
13王伟平，苏高利，杨海鹏.城北高架点源对西湖风景区影响的研究.气象,2002,28(1):13~17
14国家环境保护局.中华人民共和国环境保护行业标准之环境影响评价技术导则)))大气环境.北京：中国环境科学出版社, 1993.24~26
15H J/T2.1~2.3)93，中华人民共和国环境影响评价法实施手册.北京：银声音像出版社,2002.277~278
16钱杰，俞立中.上海市化石燃料排放二氧化碳贡献量的研究.
上海环境科学,2003,22(11):836~839
17陆德超，池钊伟.上海地区开发风力发电的必要性和南汇风力发电示范场场址选择.风电场开发,1997,1(3):25~29责任编辑：闵怀(修改稿收到日期:2005-05-21)
(上接第514页)
6Field R W,W u D,H ow ell J A,et al.Critical flux concept for m-i crofiltration foulin g.Journal of M emb rane S cien ce,1995,100: 259~272
六一儿童节简短作文7Pierre L C,Bru ce J,Chang I S,et al.Critical flux determin ation b y th e flu x-step meth od in a subm erged mem bran e bioreactor.Jour-nal of M embrane Science,2003,227:81~93
8Defrance L,J affrin M Y.Comparison between filtrations at fixed transm emb rane pres sure and fix ed perm eate flux:application to a m emb rane bioreactor us ed for w astew ater treatment.Journal of M emb rane S cience,1999,152:203~210
9Song L F,Pierre L C.M embrane bireactors for w astewater treatment.
Journal of Environment Engi neeri ng,2002,128:1018~1029
10L ianfa S.Flux decline in cros sflow microfiltration and ultrafiltra-tion:mechanisms and m odeling of membrane fouling.J ou rnal of M emb ran e Science,1998,139:183~200
11L im A L,Bai R B.M embrane fouling and cleaning in microfiltra-tion of activated sludge w astew ater.J ou rnal of M embrane Sc-i ence,2003,216:279~290
12Ognier S,Wisniew s ki C,Grasmick A.M embrane fouling in s ub-critical filtration conditions:a local critical flux con cept.J ou rnal of M em bran e Science,2004,229:171~290
责任编辑：贺锋萍(修改稿收到日期:2005-05-08)
# 522 #
