活性污泥工艺(CSTR·Monod型)模拟器 返回
环境工程

活性污泥工艺(CSTR·Monod型)模拟器

以Monod方程和完全混合槽(CSTR)稳态物质衡算为基础,计算污水处理的主役——活性污泥法。通过滑块调整进水BOD、反应槽体积、SRT等参数,出水BOD、MLSS、余剩污泥量、F/M比实时更新,直观掌握曝气槽的设计和运行条件。

参数设置
进水BOD S₀
mg/L
处理厂进口的生化需氧量
反应槽体积 V
曝气槽的有效容积
流量 Q
m³/day
处理厂的进水流量
SRT(污泥滞留时间)
day
活性污泥在系统内停留的平均天数 θ_c
比污泥收率 Y
kg/kg
单位基质产生的微生物质量
最大比增殖速度 μ_max
1/day
基质过剩条件下的微生物增殖速度
半饱和常数 K_s
mg/L
μ=μ_max/2时的基质浓度
衰减系数 k_d
1/day
内生呼吸导致的自分解速度
计算结果
HRT (day)
出水BOD (mg/L)
BOD去除率 (%)
MLSS 污泥浓度 (mg/L)
余剩污泥 (kg/day)
F/M比 (kg-BOD/kg-MLSS/day)
曝气槽截面 — 进水·处理·污泥回流动画

褐色有机物流入,与空气气泡和活性污泥(绿色微生物颗粒)作用分解成蓝色处理水。底部污泥回流线维持MLSS。

出水BOD vs SRT
BOD去除率 vs F/M比
理论·主要公式

$$S_{\text{eff}}=\frac{K_s\,(1/\theta_c+k_d)}{\mu_{\max}-1/\theta_c-k_d},\qquad X=\frac{Y\,\theta_c}{\theta_H\,(1+k_d\,\theta_c)}\,(S_0-S_{\text{eff}})$$

CSTR稳态出水基质浓度S_eff和MLSS(污泥浓度)X。θ_c=SRT(污泥滞留时间)、θ_H=HRT(水力停留时间)、K_s=半饱和常数、k_d=内生呼吸系数、Y=收率、μ_max=最大比增殖速度。

$$P_x=\frac{Y\,Q\,(S_0-S_{\text{eff}})}{1+k_d\,\theta_c},\qquad \text{F/M}=\frac{Q\,S_0}{X\,V}$$

余剩污泥发生量P_x(kg/day)和F/M比(kg-BOD/kg-MLSS/day)。F/M比在标准活性污泥法运行域为0.2~0.5。

活性污泥模型(Activated Sludge)简介

🙋
在污水处理厂参观时听说"活性污泥",那冒泡的池子里到底发生了什么?
🎓
那个曝气槽本质上是微生物养殖场。污水中的有机物(BOD成分)作为细菌和原生动物群落(即活性污泥)的食物,微生物增殖同时氧化分解有机物。通过通风供氧,好氧微生物可以在24小时内处理掉进水BOD的95~99%。沉砂池中污泥沉淀分离,上清液排放,沉淀污泥的一部分回流到曝气槽——这是标准的活性污泥法流程。
🙋
明白了!微生物养殖场。那模拟器的滑块里有个"SRT"是什么意思?和流量不同吗?
🎓
好问题。由流量决定的停留时间叫HRT(水力停留时间),是水分子在槽内的平均时间。而SRT(污泥滞留时间,Sludge Retention Time)是微生物在系统内的平均时间,通过调节回流和抽吸比例可以独立控制。比如HRT=0.5天(12小时),但通过减少污泥排出,可以达到SRT=10、20天。SRT越长,微生物越成熟,能降解越低浓度的BOD,出水水质越好。这是活性污泥工艺设计最关键的参数。
🙋
那SRT越长越好,一直延长不就行了?
🎓
这是个很有意思的问题。SRT太长会导致微生物自分解(内生呼吸k_d),污泥"老化"。老化污泥沉降性变差,沉砂池无法有效分离,污泥流出。而且F/M比(食物对微生物比)降得太低,会导致丝状菌增殖引起污泥膨化。标准活性污泥法中F/M=0.2~0.5是最佳值。你可以用本模拟器试试把SRT调到20、25天,看F/M比是怎么下降的。
🙋
Monod方程中的K_s也想了解一下。这代表什么?
🎓
K_s是"半饱和常数",当微生物的比增殖速度降到μ_max的一半时的基质(BOD)浓度。K_s越小,微生物对低浓度BOD的处理能力越强。污水处理菌K_s≈30~100 mg/L,甲烷菌数百mg/L,硝化菌<1 mg/L差异很大。看公式S_eff=K_s(1/SRT+k_d)/(μ_max−1/SRT−k_d),K_s和出水BOD成正比。可以说K_s体现微生物对基质的亲和力。冷水期K_s往往增大,这解释了冬季出水水质变差的原因。
🙋
最后确认一下,这个模型是不是和我听说的"ASM1"是同一个?
🎓
好眼力。真正的ASM1(IWA活性污泥模型第1号)把有机物分为"溶解/颗粒"和"易降解/难降解"四类,还包括从属营养菌、硝化菌、无活性物质和氮成分,共13个变量。本工具简化到"仅碳除去""单一基质""单一菌群"——这是教学用的简化版。实际工程和工厂设计用GPS-X、WEST等商用软件实现的ASM1~ASM3。但要理解Monod×CSTR的本质,本工具的滑块是最快的办法——试试把进水BOD和SRT翻倍,看出水BOD和MLSS怎么变。

常见问题

在完全混合活性污泥反应槽中,出水基质浓度S由SRT(污泥滞留时间)和Monod方程唯一确定。微生物的比增殖速度μ=μ_max·S/(K_s+S)减去内生呼吸k_d的值与1/SRT相平衡,因此S = K_s·(1/SRT + k_d) / (μ_max − 1/SRT − k_d)。出水BOD不依赖于进水BOD,仅由SRT和动力学参数决定,这是CSTR模型的重要特征。SRT越长,S越小,去除率越高。
稳态MLSS(mg/L)通过X = Y·SRT/(HRT·(1+k_d·SRT))·(S₀−S)计算。Y是微生物收率,HRT是水力停留时间。SRT越长,X越大,但受k_d自分解的限制。余剩污泥发生量(kg/day)为Y·Q·(S₀−S)/(1+k_d·SRT)/1000,SRT延长可减少余剩污泥,这是污泥处置成本削减的关键。但需注意F/M比不能过低。
F/M比(食物对微生物比,kg-BOD/kg-MLSS/day)是活性污泥工艺的典型运行指标,标准活性污泥法一般为0.2~0.5。F/M=0.05~0.15为低负荷运行(长时间曝气、氧化沟),污泥老化导致沉降性变差。F/M>0.5为高负荷过载,易引起丝状菌增殖导致污泥膨化。本工具自动计算F/M比,可观察SRT和反应槽体积变化的影响。
当SRT低于微生物最小增殖时间(冲洗出SRT)时,活性污泥持续流出,槽内微生物消失,BOD去除无法进行。临界值为SRT_min = 1/(μ_max·S₀/(K_s+S₀) − k_d),进水BOD越高,该值越小。实际工程中宜采用SRT_min的2~3倍以上作为安全设计。工具在SRT接近临界值时显示警告,可视化出水BOD急剧恶化。冷水期μ_max降低,冬季需延长SRT进行运行控制。

现实世界中的应用

城市污水处理厂(标准活性污泥法):日本国内2200余座污水处理厂大多采用活性污泥法工艺。进水BOD 150~300 mg/L、HRT 6~8小时、SRT 5~10天是代表性设计值,本工具的计算公式(Monod×CSTR)直接用于处理厂初步设计。用滑块把流量Q翻倍,可观察HRT减半、出水BOD急剧恶化的过程,直观说明雨天越流(CSO)对策的必要性。

食品、化工、造纸等工业废水处理:进水BOD 1000~5000 mg/L的高浓度废水,一般在活性污泥前段设置厌氧处理(UASB等),后段用活性污泥法精处理。把工具中S₀调到2000 mg/L,可见MLSS大幅增加,但实际工厂当MLSS超过4000 mg/L时,沉砂池固液分离会破坏,需要改用膜分离活性污泥法(MBR)或稀释处理。

运行成本优化和节能设计:活性污泥工艺的运行成本中约60%用于曝气动力(供氧),约20%用于污泥处理。延长SRT可减少余剩污泥,但增加耗氧量(曝气动力)和污泥老化风险——这是须每天权衡的因素。用本工具对比SRT=5天和20天,可见余剩污泥减半的同时,MLSS增加、沉砂池负荷增加。

高校和高专环境工程教育:Monod方程和CSTR物质衡算是环境工程、生化工程基础教材(如Metcalf & Eddy《污水工程》)的必教内容。本工具通过可视化这一经典解析解,让学生瞬间体验"8参数灵敏度分析"——这在纸笔计算上很困难。可用作毕业设计学习ASM1前的基础。

常见误解和注意事项

第一个误解:"改变进水BOD时出水BOD不变——计算是不是出了问题?"CSTR稳态模型中,出水基质浓度S_eff只由SRT、μ_max、K_s、k_d决定,与进水BOD S₀无关。这是数学性质,进水增加时微生物量MLSS按比例增加,吸收增加的负荷。实际工厂中MLSS有上限(沉砂池能力)、供氧能力有限,高进水BOD会导致处理破坏,但这些约束是模型范围外的。S₀灵敏度平坦不是计算错误。

第二个误解:"文献中的Monod参数μ_max、K_s、k_d、Y直接用就能重现实际工厂"。这些参数随温度(20℃和10℃相差一倍)、pH值(偏离中性急剧下降)、有机物组成(易降解vs难降解)、有毒物质有无而大幅变化。实际工程设计必须用目标废水进行室内回分试验(OUR、SOUR测定)来实测这些参数。本工具初始值是基于标准城市污水,不适用工业废水。

第三个误解:"本模型只处理碳除去,不计算氮磷去除"。实际高度处理中,好氧硝化(NH₄⁺→NO₃⁻)、缺氧脱氮(NO₃⁻→N₂)、厌氧聚磷菌(PAO)的厌氧聚磷等机制相结合,涉及ASM2/ASM2d/ASM3,变量增至13~19个。硝化菌μ_max很小,SRT需要8天以上,低温敏感,与本工具的简单BOD除去完全不同。涉及氮磷去除的现代设计必须用ASM2d以上的模拟器。

使用指南

  1. 输入进水BOD(mg/L)、反应槽体积(m³)、流量(m³/day)、污泥滞留时间SRT(day),开始模拟
  2. 基于Monod型反应方程μ=μ_max×S/(K_s+S),计算稳态微生物增殖速度
  3. 出水BOD、去除率、MLSS浓度、余剩污泥量、F/M比实时更新

具体计算示例

进水BOD 300mg/L、反应槽20m³、流量10m³/day、SRT 10day的标准活性污泥法运行条件下:采用HRT=2.0day、μ_max=0.5/day、K_s=50mg/L、收率Y=0.6、衰减系数k_d=0.05/day参数,稳态出水BOD=15mg/L(去除率95%)、MLSS=3500mg/L、余剩污泥=18kg/day、F/M比=0.48kg-BOD/kg-MLSS/day。表示有机物负荷在适正范围内。

实务注意事项