河流 DO/BOD 水质模拟器 返回
环境工程

河流 DO/BOD 水质模拟器 — Streeter-Phelps 方程

当污水或工业废水排入河流,微生物分解有机物会大量消耗溶解氧(DO),使下游 DO 出现明显下垂。基于 Streeter-Phelps 方程,平衡脱氧与复氧两个一阶反应,定位下游临界点,并立刻判断排放负荷是否满足水质标准。

参数设置
极限 BOD L₀
mg/L
排放口下游侧的最终(极限)BOD,反映污水负荷大小
饱和 DO D_s
mg/L
由水温与海拔决定的平衡 DO(20℃ 时约 9.1 mg/L)
初始 DO 亏量 D₀
mg/L
混合后 D_s 与实测 DO 的差,即起点值
脱氧速率 k_d
1/day
微生物分解有机物的一阶反应速率,20℃ 通常为 0.1—0.5
复氧速率 k_r
1/day
大气向水面溶氧的一阶速率,由流速、水深与紊动决定
流速 v
km/day
平均流速,将时间 t 换算为下游距离 x = v·t
计算结果
临界时间 t_c (day)
临界距离 x_c (km)
临界亏氧 D_c (mg/L)
临界 DO (mg/L)
1 天后 DO (mg/L)
水质判定
河流断面示意 — 下游 DO 下垂

上游在左,下游在右。色彩深浅表示 DO 浓度(深蓝=高 / 浅=低);临界点位置标有鱼形标志。

DO 下垂曲线 — 下游距离 vs DO / 亏量
负荷敏感性 — 临界 DO 对极限 BOD L₀
理论与主要公式

$$D(t) = \frac{k_d L_0}{k_r - k_d}\left(e^{-k_d t} - e^{-k_r t}\right) + D_0\,e^{-k_r t}$$

下游 DO 亏量 D(t)[mg/L]。L₀:极限 BOD,D₀:初始亏氧,k_d:脱氧速率,k_r:复氧速率。实际 DO = D_s − D。

$$t_c = \frac{1}{k_r - k_d}\ln\!\left[\frac{k_r}{k_d}\left(1 - \frac{D_0(k_r - k_d)}{k_d L_0}\right)\right], \qquad x_c = v\,t_c$$

临界时间 t_c 与下游临界距离 x_c。v:平均流速 [km/day]。k_d 与 k_r 的平衡决定临界点位置。

$$D_c = \frac{k_d}{k_r}\,L_0\,e^{-k_d t_c}, \qquad \text{DO}_c = D_s - D_c$$

临界亏氧 D_c 与临界 DO,由 dD/dt = 0 推得。保证 DO_c 高于环境标准(暖水鱼类约 5 mg/L)是排放许可的核心条件。

河川的 DO/BOD 水质模型(Streeter-Phelps)

🙋
污水排入河流后鱼为什么会死?是被「脏」直接毒死的吗?
🎓
绝大多数情况下不是直接中毒,而是「缺氧」。污水里的有机物(蛋白质、糖、油脂等)被水中的微生物吃掉时,它们会用呼吸消耗大量溶解氧,这一过程叫脱氧。同时大气中的氧气也会缓慢溶进水里——叫复氧。只要脱氧速度超过复氧速度,DO 就会持续下降,到某一点达到最低,这就是「临界点」。低于约 2 mg/L 时大多数鱼就无法呼吸,会大量死亡。
🙋
所以是 Streeter-Phelps 方程预测这条下垂曲线?1925 年的式子现在还够用吗?
🎓
没错。Streeter 和 Phelps 当年从俄亥俄河的数据拟合出这个公式,是水质工程的「祖师」公式。但今天在污水处理厂排放许可、河流环境标准评估里依然是第一个用上的工具。原因有三:它有闭式解析解,能一眼看出哪个参数最敏感;计算瞬时完成,适合做敏感性分析;并且能很好地解释脱氧与复氧的对抗物理。正式设计会用 QUAL2K、WASP 这样的数值模型,但工程师还是会用 Streeter-Phelps 来核对结果是否合常理。
🙋
我把左边的极限 BOD L₀ 从 30 调到 100,临界 DO 变成负值了。这该怎么解读?
🎓
观察得很好。Streeter-Phelps 在数学上不会区分「DO 不可能小于 0」这一物理事实——只要负荷过大,公式就会算出负的 DO。实际 DO 不会真的负值,但负值预测是一个明确信号:这条河会进入厌氧状态。在厌氧条件下,另一类微生物会还原硫酸盐释放硫化氢 H2S(那种臭水沟「臭鸡蛋味」),还会产生甲烷。对工程师来说这就是「设计不通过、必须降低负荷或升级处理」的红灯。本工具在临界 DO 低于 2 mg/L 时会给出「严重污染/厌氧风险」的警告。
🙋
我调流速 v 时,临界「时间」不变但临界「距离」会变。这说明急流山溪和缓流大江的结果差别很大吗?
🎓
对。流速 v 只把时间 t 缩放成下游距离 x = v·t,对反应速率本身没有直接影响。但现实中 v 大的河流通常水面紊动强、复氧 k_r 也大;O'Connor-Dobbins 经验公式给出 k_r ≈ 3.93·v^0.5·H^(-1.5),由流速和水深预测复氧。所以急流河虽然临界距离远,但因 k_r 大、下沉浅;缓流大江临界距离近,但 k_r 小、下沉深、恢复慢。这正是日本上游的山溪自净能力强、下游城市河流(如多摩川下游)最痛苦的典型规律。
🙋
那设计污水处理厂时,是用这个工具来反推 L₀ 必须降到多少吗?
🎓
实务上其实多一步:先做物料平衡,把「污水流量 Q_w・BOD C_w」和「河流流量 Q_r・BOD C_r」按质量混合,得到排口下游的 L₀ 与 D₀;再代入本工具算临界 DO,检验是否高于环境标准(中国地表水 II 类约 6 mg/L,III 类 5 mg/L,IV 类 3 mg/L)。如果不达标,就在「升级处理(活性污泥→A2O→三级处理)」「迁移排口到流速更快的河段」「合流制下水道改分流制」等方案中比较成本。有时候最便宜的方案竟然是「把排口移到有跌水的急流段」——靠复氧的「免费红利」,工艺等级反而可以降一级。

常见问题

由 Streeter 与 Phelps 于 1925 年根据美国俄亥俄河数据建立的解析解,描述含有机物的污水排入河流后下游溶解氧(DO)随时间的变化。它把微生物分解有机物造成的脱氧和大气向水面溶氧的复氧都视为一阶反应,得到闭式解 D(t) = (k_d·L₀)/(k_r−k_d)·(exp(−k_d·t)−exp(−k_r·t)) + D₀·exp(−k_r·t)。模型虽简,但在污水排放许可、纳污容量核定、河流水质达标评估中至今仍是首选工具,并构成 QUAL2K 等专业水质模型的氧平衡核心。
DO 亏量 D(t) 起初因脱氧占优而增大,待复氧追上后转为下降,最大值就是「临界点」,由 t_c = ln[(k_r/k_d)·{1 − D₀(k_r−k_d)/(k_d·L₀)}]/(k_r−k_d) 给出。该点 DO 是全河段最低的氧值,决定了「鱼类能否生存」「会不会厌氧化释放硫化氢和甲烷」等设计性结论。污水处理厂的排放许可往往就以预测临界 DO 是否高于环境标准(如 5 mg/L)作为判定依据。
20 ℃ 下脱氧速率 k_d 大约为 0.1—0.5 day⁻¹(以生活污水为主的城市河流多为 0.2—0.4)。复氧速率 k_r 与水力条件密切相关:缓流大河 0.2—0.5;有坡度的中型河流 0.5—1.5;有急流或瀑布的陡峭河段 2—10 day⁻¹。一般只要 k_r > k_d,下游 DO 一定会恢复;若 k_r < k_d 则复氧跟不上,河流将陷入长期缺氧。温度升高会同时增大 k_d(θ=1.047)和 k_r(θ=1.024),但饱和 DO 反而降低,故夏季是最严格的设计工况。
真实河流还涉及底泥耗氧(SOD)、藻类光合与夜间呼吸、硝化反应的额外耗氧、随流量变化的复氧系数(O'Connor-Dobbins 公式 k_r=3.93·v^0.5·H^(-1.5))、支流稀释等因素,因此单纯的二项 Streeter-Phelps 在精细设计中精度不足。但它仍是 (1) 排放口近场快速筛查、(2) 复杂数值模型的合理性自检、(3) 教学场景下的首选工具。正式设计往往使用 QUAL2K、WASP 或 MIKE11,而工程师仍用 Streeter-Phelps 来核对结果是否「物理上合理」。

实际应用场景

污水处理厂排放许可与提标改造:各国水污染防治法规通常都要求「混合后的临界 DO 不低于受纳水体的水质标准」。中国地表水环境质量标准 GB3838 中,I—II 类要求 DO ≥ 6 mg/L,III 类 ≥ 5 mg/L,IV 类 ≥ 3 mg/L,V 类 ≥ 2 mg/L;日本 A 类型 7.5 mg/L、B 类型 5 mg/L 也是类似设定。若 Streeter-Phelps 预测临界 DO 不达标,处理厂就需要从「活性污泥」升级到 A2O 厌氧—缺氧—好氧组合、膜生物反应器 MBR、或臭氧氧化等高级工艺,进一步削减 BOD、氮和磷。

合流制下水道溢流(CSO)的影响评估:老旧城市采用雨污合流管网,暴雨时未处理污水会直接溢流入河。CSO 事件中 L₀ 短时间内可飙升至 100 mg/L 以上,Streeter-Phelps 模拟会给出临界 DO 趋近于零的结果,这正是东京湾入海河、巴黎塞纳河、伦敦泰晤士河旧城段历史性鱼类大规模死亡事件的机理。如今城市普遍建设深埋 CSO 调蓄隧道(东京的环七地下调节池、伦敦的 Thames Tideway)并逐步把合流制改为分流制,DO 计算就是这些巨额投资的科学依据。

工业排放纳污总量核定:食品、造纸、酿造、畜牧等高 BOD 排放行业,需要根据受纳河流的流量和现状负荷反推「企业被允许排放的最大 BOD」。用 Streeter-Phelps 模拟「现状临界 DO + 本企业新增排放后的临界 DO」,再反算保证达标的最大 L₀。新厂选址、扩建论证、ISO 14001 环境影响评价都会用这套方法。

城市河流复兴与自净能力可视化:韩国首尔清溪川、日本东京隅田川、英国伦敦泰晤士河等城市河流复兴项目,往往通过疏浚、增加流量、建设人工浅滩等方式有意提高复氧系数 k_r。仅调整 k_r 一项,用 Streeter-Phelps 展示下游 10 km 处 DO 提高 2 mg/L 的结果,对公众与决策者具有极强的说服力。一座桥下添置小跌水使 k_r 翻倍,从而带来下游显著改善——这种定量诉求是项目立项的关键论据。

常见误区与注意事项

第一个误区是把「BOD 达标」等同于「水质达标」。BOD 只是「这些有机物最终会消耗多少氧」的预测指标,并不是当前的 DO。即使 BOD 不高,在流速慢、复氧弱的河段,或底泥厚、SOD 大的河段,DO 也可能跌破标准。BOD 是负荷指标,DO 才是鱼类生存的最终判据。本工具想强调的核心信息正是:设计判定必须基于 DO,而不仅仅是 BOD。

第二个陷阱是忽视季节与温度的影响。本工具允许独立输入 k_d、k_r 和饱和 DO,但现实中它们随温度联动变化。盛夏 25 ℃ 时 k_d 约为 20 ℃ 的 1.3 倍,k_r 略增,但饱和 DO 从 9.1 降至 8.4 mg/L。结果临界 DO 大幅降低,这就是每年 7—8 月鱼类大量死亡「夏季危机」的机理。设计时应取最高水温下的最不利工况,而非全年平均工况。西班牙埃布罗河、美国南部河流在夏季干旱时段也年复一年地出现类似事件。

第三个常见误解是把「Streeter-Phelps 给出负 DO」当作「公式坏掉」。实际 DO 不可能为负,但 Streeter-Phelps 是线性模型,本质上回答的是「假定氧无限供给的话还会差多少」。负 DO 预测是「真实情况下河流会厌氧化、改由完全不同的化学过程主导」的明确警告。一旦出现负 DO,应立即判为「设计不通过、存在厌氧化风险」,去思考削减负荷或改善复氧的方案。本工具在临界 DO 低于 2 mg/L 时已自动给出「严重污染/厌氧风险」警告。对负 DO 工况,更严谨的做法是把 DO 截断在 0、并改用厌氧分解模型。

使用指南

  1. 输入废水BOD浓度(mg/L,典型范围10-200)和河流饱和溶解氧(mg/L,20°C时约8.6),模拟器基于Streeter-Phelps方程计算脱氧速率k₁与复氧速率k₂
  2. 设置初始亏氧D₀(mg/L)与脱氧化速率常数k₁(/day,典型值0.1-0.5),复氧化速率k₂由河流流速与深度自动关联
  3. 点击计算后获得临界时间t_c、临界距离x_c、最小溶解氧浓度及对应的水质等级判定,判断排污口下游是否满足GB3838二级或三级标准(DO≥5或≥3 mg/L)

具体计算示例

某制药厂排污口:BOD=80 mg/L,河流饱和DO=8.5 mg/L,初始亏氧D₀=6.2 mg/L,脱氧速率k₁=0.35/day,复氧速率k₂=0.55/day。按Streeter-Phelps方程D(t)=k₁L₀/(k₂-k₁)×[exp(-k₁t)-exp(-k₂t)]+D₀×exp(-k₂t),计算得临界时间t_c≈2.1天,临界距离x_c≈42 km(河流流速20 km/day),最小DO≈3.8 mg/L,判定为三级水体。1天后DO恢复至4.6 mg/L。

实务注意事项

  1. 脱氧速率k₁随温度变化显著,20°C标准值需用θ^(T-20)修正(θ=1.047),冬季低温延长临界时间,夏季高温加剧DO枯竭
  2. 复氧速率k₂与河道比降、糙率相关,浅滩急流k₂可达0.8-1.2/day,平缓河段仅0.2-0.4/day,需采集水文参数准确输入
  3. 实际排污需控制BOD≤150 mg/L并预处理至50 mg/L以下,避免临界点下移至取水口;多污染源叠加时按等效BOD累加计算