基于malab的碳氧化aq活性污泥法模型.docxVIP

基于malab的碳氧化aq活性污泥法模型 在许多生物污污泥数学模型中，国际水公司（iwq，现改为国际水公司iwa）开发的生态污污泥数学模型（asmno.1.3）最为广泛。但是,根据我国目前污水处理厂的设计、运行和水质检测水平,直接应用该套模型还有一定的困难。因此,作者在学习和总结IAWQ活性污泥数学模型的基础上,根据国内典型活性污泥法污水处理厂的工艺状况,对IAWQ活性污泥模型作了适当简化,建立了适合于普通推流式活性污泥法的碳氧化数学模型,并应用MATLAB数学软件在WINDOWS操作平台上开发出一套模拟系统。该模拟系统考虑了8个组分、3个生化过程,为了便于模拟系统在污水处理厂设计和运行中的应用,还开发了模型参数优化模块、常规检测数据与模型数据之间的转换模块等。应用该模拟系统模拟了天津纪庄子污水处理厂的实际运行资料,检验了模型的实用性能。 1 数学模型的构建 1.1 盐废水处理系统—模型的假设条件 数学模型主要针对传统推流式活性污泥法,为了简化模型,仅考虑异养菌好氧条件下的碳氧化过程和颗粒性有机物质的水解过程等,而对其他过程暂时不予考虑。因此,该模型称之为活性污泥工艺模拟系统之碳氧化篇(Activated Sludge Process Simulation-Carbon Oxidation Version),简称为ASPS—CO模拟系统。ASPS—CO模拟系统的具体假定条件为: a.系统在恒定温度下运行; b.pH值恒定且接近中性; c.速率表达式中系数恒定,不随废水成分的变化而变化; d.不考虑氮、磷和其他无机营养物对有机物去除及细胞生长的影响; e.异养微生物均一且种属稳定; f.颗粒性有机物质被微生物吸附过程是瞬间完成的; g.假设二沉池不发生任何生化反应,仅为物理沉淀作用。 1.2 生物化学过程和组成部分矩阵的矩阵 表1 给出了ASPS—CO模拟系统中包含的8个组分和3个生化过程的关系,表2给出了各组分符号定义。 1.3 系统活性污泥浓度及催化物系统泥龄的表示 普通活性污泥法通常采用长方廊道推流式曝气池,进水和污泥回流在反应器的一端汇合。为便于模拟,通常这类反应器可假定为由3～5个连续、混合反应器串联而成。此处选择了5池串联模型,其流程见图1。 其中,F为废水进水流量,Fw为剩余污泥排放流量,Fr为回流污泥流量。反应器有效容积用V表示。其他符号同表2。 ASPS—CO模拟系统中将5个反应器分别作为控制体,对其分别按活性污泥量、易降解基质、颗粒性有机基质和惰性颗粒生物量建立质量守恒方程组,通过数值逼近方法求得方程组的近似解(利用MATLAB数学工具软件来实现)。在稳态条件下,各组分的质量平衡方程式通式为: FCO+F1C′+rV=0(1) 式中F、F1——分别为进、出该级反应器的流量 CO、C′ ——分别为进、出该级反应器的组分浓度 r——反应速率 V——反应器体积 例如,1级反应器的质量平衡方程组为: Fr(F+Fr)FW+FrXBΗ5-(F+Fr)XBΗ1+(?μΗ?SS1ΚS+SS1SΟ1ΚΟΗ+SΟ1XBΗ1-bLΗXBΗ1)V1=0(2)FSSΟ+FrSS5-(F+Fr)SS1+[(-1YΗ)?μΗ?SS1ΚS+SS1SΟΚΟΗ+SΟXBΗ1+khXS1ΚX+XS1/XBΗ1?SΟ1ΚΟΗ+SΟ1]V1=0(3)FXSΟ+Fr(F+Fr)FW+FrXS5-(F+Fr)XS1+[(1-f′D)bLΗXBΗ1-khXS1ΚX+XS1/XBΗ1SΟ1ΚΟΗ+SΟ1]?V1=0(4)FXΙΟ+Fr(F+Fr)FW+FrXD5-(F+Fr)XD1+[f′DbLΗXBΗ1]V1=0(5) 其余4级反应器的质量平衡方程组与以上方程组相似。 系统泥龄的表示方程为: Θc=5∑i=1[Vi(XBΗi+XSi+XDi)]FW(F+Fr)FW+Fr(XBΗ5+XS5+XD5)(6) 式中XBH1、XBH5、XBHi——分别为1级、5级和第i级反应器中活性污泥浓度,mg/L SSO、SS1、SS5——分别为进水、1级和5级反应器中的溶解性易降解有机物浓度,mg/L XSO、XS1、XS5、XSi——分别为进水、1级、5级和第i级反应器中的颗粒性生物可降解有机物浓度,mg/L XIO、XD1、XD5、XDi——分别为进水及1级、5级和第i级反应器中的惰性颗粒性有机物浓度,mg/L Θc——污泥龄,d 其余符号的含义见表2。 将5级反应器的20个质量平衡方程和泥龄联立方程组用MATLAB数学模拟软件进行数值求解。 1.4 如图2所示，asps-co模拟系统的程序图 ASPS—CO模拟系统可通过菜单随时调整模型参数值,并终止程序的运行。 2 模型参数值的修正 污水处理厂目前最普遍的检测数据仍然局限于COD(多为CODCr)、BOD(多为B