MBBR生物预处理工艺硝化过程动态模型的建立.docx
- 文档编号:9601857
- 上传时间:2023-05-20
- 格式:DOCX
- 页数:16
- 大小:149.31KB
MBBR生物预处理工艺硝化过程动态模型的建立.docx
《MBBR生物预处理工艺硝化过程动态模型的建立.docx》由会员分享,可在线阅读,更多相关《MBBR生物预处理工艺硝化过程动态模型的建立.docx(16页珍藏版)》请在冰点文库上搜索。
MBBR生物预处理工艺硝化过程动态模型的建立
第38卷 第5期
2006年5月
哈 尔 滨 工 业 大 学 学 报
JOURNALOFHARBININSTITUTEOFTECHNOLOGY
Vol138No15May2006
MBBR生物预处理工艺硝化过程动态模型的建立
徐 斌,夏四清,胡晨燕,高乃云
(同济大学污染控制与资源化研究国家重点实验室,上海200092,E2mail:
tjwenwu@mail.tongji.edu.cn
摘 要:
采用MBBR工艺对微污染黄埔江原水硝化过程动力学和反应器动力学进行了研究.采用考虑最小基质质量浓度的Michaelis-Menten方程,氨氮去除速率方程为N=1879(S-0108/(S+1152,各动力学参数分别为:
最
大氨氮去除速率Nmax为1879mg/(m2
・d;半速率常数Ks为116mg/L;最小基质氨氮质量浓度为0108mg/L;最大细胞比增长速率μmax为1105d-1
.试验表明,MBBR反应器为典型的完全混合反应器,结合硝化反应动力学,建立了
完全混合式原水生物预处理硝化反应器动态模拟模型.通过模型计算与实际中试运行效果比较可以得出:
在低进水氨氮质量浓度条件下,模型计算值与实验数据略有差别;在中高进水氨氮质量浓度条件下,模型计算值与实验结果较为一致.模型较好地反应了工艺硝化的过程,可方便地应用于工艺的控制和管理.关键词:
硝化;动态模拟;生物预处理;MBBR中图分类号:
X506
文献标识码:
A
文章编号:
(2006-05
Dynamicsiafioninbio2pretreatment
gMBBRtechnology
XUBin,XIASi2qing,HUChen2yan,GAONai2yun
(StateKeyLaboratoryofPollutionControlandResourceReuse,TongjiUniversity,Shanghai200092,China,
E2mail:
tjwenwu@mail.tongji.edu.cn
Abstract:
Thenitrifyingandreactorkineticsofrawwaternitrificationinbiologicalpretreatmentofmicro2pol2lutedwaterfromHuangpuRiverbyMBBRtechnologywasstudied.AccordingtotheequationofMichaelis2Menten,takingtheminimumammonianitrogenconcentrationintoconsideration,therelationshipbetweenam2monianitrogenconcentrationsandammonianitrogenremovalratescouldbepresentedasN=1879(S-0108/(S+1152.Fournitrificationkineticsparameterscouldbedirectlyorindirectlyobtainedbythemeas2urementsofthereactorsystem,includingtheminimumammonianitrogenconcentration,thehalfsaturationconstant,themaximumammoniaremovalrateandthemaximumspecificbacterialgrowthrate.Fortheexperi2
mentalconditions,theparametersweredeterminedtobe0108mg/L,116mg/L,1859mg/(m2
・dand1148d
-1
respectively.ItcanbedrawnthattheMBBRbelongstotheCSTRandconsideringthenitrifyingki2
neticsthenitrificationdynamicsimulationmodelofthistypeofreactorwasbuilt.Thevalidityofthismodelwastestifiedbythepilot2scaleexperimentofHuangpuriverrawwaterbio2pretreatmentbyMBBRtechnology.Comparingthecalculationresultsofthemodelwiththeresultsofexperiment,theconclusioncouldbedrawn:
thehighertheconcentrationofammonianitrogeninrawwateris,thelessthedifferencebetweenthemodelandtheexperimentis.Asamatteroffact,theresultsofthismodelareveryclosetotheexperiments;anditcanbeusedtomanageandcontrolthepracticaloperationofthistechnologyeffectively.Keywords:
nitrification;dynamicsimulationmodel;biologicalpretreatment;MBBR
收稿日期:
2004-08-22.
基金项目:
国家高技术研究发展计划资助项目(2002AA601130;
国家科技攻关计划重大资助项目(2003BA808A17.
作者简介:
徐 斌(1976-,男,博士,讲师;
夏四清(1965-,男,博士,教授,博士生导师;高乃云(1949-,女,博士,教授,博士生导师.
原水中的氨氮在水的输送和处理过程中,容易造成管网中亚硝酸盐质量浓度的增高、除锰困难和降低有机物氧化效率等问题;另外原水中氨氮质量浓度过高,需加大氯的投加量,从而引起有
害副产物增加[1]
.因此,如何有效去除原水中氨
氮也是目前预处理工程的主要任务之一.生物硝化是去除水中氨氮的有效方法之一,它具有成本低、效率高等特点.微污染原水中的氨氮质量浓度一般<3mg/L,与研究较多的污水生物硝化动力学有较大的不同[2,3]
.如何正确建立原水生物处理过程中硝化动力学及其动态模拟模型,对于深入研究原水预处理过程和指导工程实践是十分必要的.移动床生物膜反应器(MBBR是在生物滤池和流化床的工艺基础上发展起来的一种水处理新技术.它具有处理能力高、能耗低、不需要反冲洗、水头损失小、不发生堵塞的工艺特点,已广泛应用于污废水处理以及城市污水处理场的改造方面,但应用于微污染原水处理中尚少见报道.本文利用MBBR工艺对黄浦江微污染原水进行了为期1年,规模为60~160m3
/h的生物预处理中试实验.对该工艺硝化过程动态模型进行研究,以期
能为该工艺的实际应用提供参考.
1 实验工艺与设备
MBBR(钢板制尺寸为:
长×宽×高=610m×115m×915m(有效水深910m,池体在
长度方向上平均分为串联式两格;池内按填充率50%投加同济大学专利悬浮填料(填料为球型,
单体直径为100mm,比表面积为100m2/m3
;池底设置48根DN25穿孔曝气管(L=500mm,并均匀分为4组.实验时使用两组曝气管曝气,以保证填料的悬流状态.
2 氨氮的降解规律
试验前MBBR池运行效果良好,各进出水水质指标如表1所示.,,探.
水质指标水温
/U(DO/(mg・L-1
/(mg・L-1
ρ(NH3-N
/(mg・L-1
ρ(NO3
-
-N
/(mg・L-1
ρ(NO2
-
-N
/(mg・L-1
进水2011138355126108113521940114出水
2011
7161
36
918
5192
0135
3185
0105
在MBBR池两格分别投加1L25%工业氨
水,并保持气量为45m3/h,池体内水量为80m3
.氨氮、亚硝酸盐、硝酸盐、pH与DO平均质量浓度随时间的变化曲线如图1、2所示.
图1 氨氮、亚硝酸盐和硝酸盐变化曲线
图2 DO与pH值变化曲线
从图1、2可以看出:
随着时间的变化,氨氮在
开始阶段迅速降解,随后氨氮质量浓度降低导致降解速率逐步降低.硝酸盐的增长基本伴随着氨氮的降解.但同时可以看出,亚硝酸盐质量浓度在氨氮降解过程中,会有一个积累过程,但积累量非常小,随后会逐步降解.在降解过程中,氨氮降解量稍大于亚硝酸盐和硝酸盐的生成量之和,反应最后硝酸盐的生成量小于氨氮的降解量大致1mg/L.pH值在降解初期高氨氮质量浓度条件下,大幅度下降,随着曝气时间的延长,在氨氮基本降解结束后,由于曝气作用,逐步上升.溶解氧的变化规律与pH的变化规律一致,溶解氧质量浓度保持在7mg/L
以上.在进水氨氮质量浓度5~6mg/L条件下,该曝气量可满足硝化过程需氧量.
3 原水生物预处理硝化反应动力学
311 最小基质氨氮质量浓度
生物膜模型中有一个重要概念为最小基质质量浓度,它是保证稳态生物膜存在的必要条件.如果基质质量浓度小于最小基质质量浓度,则生物膜的损失量大于生物细胞生长那个量,从而导致生物膜不断变薄,甚至根本不形成生物膜.Ritt2
man和McCarty提出一般维持稳态生物膜的最小
基质质量浓度可以表示为Smin=Ksb/(μmax-b,
・
637・哈 尔 滨 工 业 大 学 学 报 第38卷
μmax=Yqm
[4]
.Smin为维持稳态生物膜的最小基质质量浓度,mg/L;Ks为半最大比基质去除率时的基质质量浓度,又称半速率常数,mg/L;b为总比
生物膜衰减率,取0105d-1[3]
;Y为细胞消耗单位基质量的实际生长量(mg/mg;qm为最大比基质
去除率,即单位细胞量的最大比基质去除率,mg/
(mg・d;μmax为最大细胞比增长速率,d-1
.
最小基质质量浓度可以从各中试运行工况与降解试验结果得出.中试运行中,当进水氨氮质量浓度接近011mg/L以后,中试基本无去除效果;在降解试验中,随曝气时间的延长,在氨氮质量浓度降低到0108mg/L左右,继续曝气对氨氮去除效果基本无影响.通过数据分析,在水温为20℃左右,最小基质氨氮质量浓度为0108±0104mg/L.312 硝化动力学及其参数确定通过反应器物质的量平衡,氨氮去除速率可利用下式进行计算
Vdt=Q(S0-SerA0.
Nr=
Q(S0-Se
A
.
式中:
V为反应器体积,m3
;Q为进水水量,m3
/d;S0为进水质量浓度,mg/L;Se为出水质量浓度,mg/L;Nr为氨氮去除速率,
mg/(m2
・d;A为填
料的总表面积,m2
.
在降解试验条件下,氨氮去除速率与时间的变化如图3所示;氨氮去除速率随氨氮质量浓度的变化如图4所示.
图3 氨氮去除速率与时间的变化曲线
从图3、4可以看出:
在曝气初期氨氮质量
浓度为6mg/L,其降解速度最大可达到
1500mg/(m2
・d,随曝气时间延长,氨氮质量浓度逐步降低,相应的氨氮降解速度也随之迅速降低,在降解到最后,氨氮质量浓度为011mg/L,
池体内基本不发生降解.
图4 氨氮去除速率随氨氮质量浓度的变化曲线
对于单基质的降解,氨氮去除速率的变化可利用生物处理中所经常采用的Michaelis-Menten
方程来进行描述[6]
:
N=S/s.
-Ment2Nmax(S-Smin/(Ks+S-Smin.Smin可通过试验获得,在本试验中大致为0108
mg/L;为了能正确求出曲线方程,采用MATLAB
软件中优化工具包中任意曲线拟合的函数curve2fit,计算后得出的曲线方程为
N=1879(S-0108/(116+S-0108.
该曲线与试验数据点的分布如图5所示.图5曲线变化与试验数据点基本一致,计算数据与试验数据具有良好的相关性,相关性系数达到01998.通过计算可得出硝化过程最大去除速率为
Nmax=1879mg/(m2
・d,按此值计算氨氮比表面负荷为01078g/(m2
・h,与本中试试验测得
的最大氨氮比表面负荷01072g/(m2
・h较为接近(如图6.半最大比基质去除率时的基质质量浓度(半速率常数Ks
为116mg/L.利用试验计算出最小基质氨氮质量浓度,半速率常数Ks以及总比生物膜衰减率b;依据公式可计算最大细胞
比增长速率μmax=1105d-1
.
图5 氨氮去除速率随氨氮质量浓度变化相关性曲线
・
737・第5期徐 斌,等:
MBBR生物预处理工艺硝化过程动态模型的建立
图6 氨氮填料表面负荷与水温关系图
313 硝化反应控制因子
[7,8]
在微生物生活条件适宜的条件下,硝化反应的速率主要受溶解氧与氨氮两种基质的制约.当
SNH3-N>0127SO2时,氧为限制因素;当SNH3-N<
0127SO2时,氨氮为限制因素
[7,8]
;对于完全混合形
式的反应池如MBBR反应器,如果仅采用单反应池,按照目前处理水质的要求,作为饮用水源的原水氨氮质量浓度最大质量浓度标准为015只要保证生化池内DO11185mg/L,控制的要求.0mg计算,只要第一级池体内能处理到氨氮质量浓度为210mg/L,即DO质量浓度不低于714mg/L,就可满足生
物硝化不受DO质量浓度控制.本中试运行中第
一级池体氨氮出水质量浓度均<2mg/L,而DO质量浓度一般均高于7mg/L,因此,可以认为在黄浦江原水生物预处理过程中,本工艺硝化反应速率主要控制条件为氨氮质量浓度.
4 动态模拟模型的建立
按照中试运行的情况,对于单级池体,由于采
用底部穿孔管曝气带动水流与填料流动,水流紊动性大,达到全池流化效果,该反应器应属于典型的完全混合流反应器.411 反应器动力学的建立
主要针对该典型的完全混合反应器建立动力学方程.两级反应器示意图见图7.Q为进水流量;
A1和A2分别为两级反应池体内填料总表面积;S0、S1和S2分别为原水、第一和第二级池体内氨氮的
质量浓度;V1和V2分别为两级池体的体积
.
图7 两级MBBR反应器示意图
对第一级和第二级反应器建立物料平衡方程,
QS0-N1×A1=QS1,QS1-N2×A2=QS2
.
对于i级反应器其方程为
QSi-1-Ni×Ai=QSi.
依据前述内容,原水进水氨氮质量浓度较低,在较高溶解氧条件下,总反应速率的主要限定因素为氨氮的质量浓度,代入已得的降解速率公式后得
QS0-1879
S1-0108
S1+1152×A1=QS1,
QS1-1879S-0108S2+1152
×A2=QS2
.
(1
求解方程可求出最终出水S2的质量浓度值.方程组(1的第一个方程解为
S1=
8×105
[4×105
QS0-A1-6108×105
Q1×8
Q2
2
×107
QS0A1+108
0+641×106
A2
1+
4×108QA1+51914624×108Q21/2
].
如果采用多级处理时,列出多级方程组,先解第一级值,其余方程的解通过上述解的形式逐级解出.
412 模型效果的验证
模型验证采用中试试验中水温均维持在25℃左右的几个工况,几个工况下原水能够提供
硝化所需要的足够碱度及适合pH.工况1、2和3主要控制工况为:
进水流量80m3
/h,曝气量45m3
/h,水力停留时间HRT为1h,气水比0156.3个工况对氨氮去除效果如表2所示.模型预测
时,可采用中试条件下水温与氨氮填料表面负荷
指数关系式(如图6,将水温20℃时的氨氮最大去除速率转化到相应温度的值.
表2 各工况出水氨氮质量浓度和处理效果与实测值对比
编号工艺控制条件水量
t/℃
HRTρ(氨氮/(mg・L-1
进水
出水
效率/%
工况180********65013005015模型0165011807213工况280261360116013907615模型116013227918工况3802419603105016807716模型3105017007711工况4120
1815
40
2129018006216模型
2129
01940
59.0
中试试验两级池体中,填料的填充率相同,池型也相同.按照模型预测在原水氨氮质量浓度在
015、1、115、2、215、3、315和4mg/L情况下,HRT=60min条件下,出水氨氮质量浓度与去除
效率如图8所示;相同质量浓度范围,水温为20
・837・哈 尔 滨 工 业 大 学 学 报 第38卷
℃,HRT=40min条件下,出水氨氮质量浓度与去
除效率如图9所示.
图8 HRT=60min条件下模型预测值
图9 HRT=40min条件下模型预测值
从模型模拟值与实际工况运行结果分析比较可知:
采用的动力学方程与反应器方程,在停留时
间60min条件下,对于中质量浓度与高质量浓度进水氨氮质量浓度,模型预测值与实际工况运行值基本一致,而对于低进水氨氮质量浓度,模型计算值与实测值有一定差别.停留时间40min条件下,对于中高质量浓度氨氮进水质量浓度,计算值与实际测定值也基本一致.
依据以上结果可知,建立的动力学模型与动态模拟模型基本能正确描述中试过程中氨氮降解规律.该动态模拟模型的建立可方便地应用于考虑不同影响因素(如水量、进水质量浓度、填料比表面积、填料填充率、反应器级数、各级池体体积分配等对处理效果的影响,对于MBBR生物反应器处理微污染黄浦江原水运行具有较好的指导作用.
5 结 论
1建立了微污染原水硝化反应速率方程,并
计算出原水硝化过程Nmax为1879mg/(m2
・d;
Ks为116mg/L;最小基质氨氮质量浓度为
0108mg/L;μmax为1105d-1
.
2采用两级MBBR反应器进行微污染原水
硝化,控制气水比在015以上时,硝化反应控制因素为原水中氨氮的质量浓度.由于曝气作用,池体内悬浮填料不断运动,该类型反应器为典型的完全混合反应器.
3在停留时间60min条件下,对于中质量浓度与高质量浓度进水氨氮质量浓度,模型预测值与实际工况运行值基本一致,而对于低进水氨氮质量浓度,模型计算值与实测值有一定差别.停
留时间40min条件下,对于中高质量浓度氨氮进水质量浓度,计算值与实际测定值也基本一致.该,可方便地应用于MBBR.
[1].给水中的氨氮问题[J].净水技术,2001,20
(2:
12-14.
[2]LIUYu,CAPDEVILLEB.Kineticbehaviorsofnitrif2
yingbiofilmgrowthinwastewaternitrificationprocess
[J].EnvironmentalTechnology,1994,15(11:
1001-1013.
[3]KIMURAK,WATANABEY,OKABES,etal.Kinetic
analysisofnitrifyingbiofilmgrowingontherotatingmem2branedisk[J].WaterSciTechnol,2001,43(4:
111-118.
[4]RITTMANBE,MCARTYPL.Modelofsteady-state
-biofilmkinetcs[J].BiotechnolBioeng,1980(22:
2343-23574.
[5]RITTMANBE,MANEMJA.Developmentandexperi2mentalevaluationofasteady-state,multispeciesbiofilmmodel[J].BiotechnolBioeng,1992,39:
914-922.[6]刘 雨,赵庆良,郑兴灿.生物膜法污水处理技术
[M].北京:
中国建筑工业出版社,2000.
[7]徐 斌,夏四清,胡晨燕,等.微污染原水生物预处
理过程中硝化动力学模型研究[J].中国给水排水,
2003,19(4:
15-18.
[8]徐 斌,夏四清,高廷耀,等.悬浮填料床处理微污
染原水硝化试验研究[J].环境科学学报,2003,23
(6:
742-747.
(编辑 刘 彤
・
937・第5期徐 斌,等:
MBBR生物预处理工艺硝化过程动态模型的建立
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MBBR 生物 预处理 工艺 硝化 过程 动态 模型 建立