4. 进水SS与COD比值关系的稳定性评估

为什么用 SS/COD 而不是单看 SS 或 COD

进水悬浮物（SS）与化学需氧量（COD）的比值，本质上反映的是"颗粒态有机物在总有机负荷中的占比"。当进水以生活污水为主、管网运行稳定时，SS/COD 通常落在一个相对固定的区间；一旦出现管网冲刷、雨污混接、工业废水偷排或长期沉积物翻起，这个比值会先于单一指标发生偏移。因此评估它能否作为"可生化性快速预警指标"，关键不是看某一天的值，而是看它在连续运行周期内是否稳定可复现。

本节采用两个统计量来量化稳定性：

  • 均值：反映该比值在评估周期内的中心水平，用于确定预警基准带。
  • 变异系数 CV = 标准差 / 均值：无量纲，便于跨周、跨季节比较离散程度。工程上常以 CV ≤ 0.15 视为"较稳定"，0.15～0.25 为"需关注"，> 0.25 为"不稳定、不宜单独作预警阈值"。

下面用一段可运行的 Python 代码完成计算，示例数据为 28 天（4 周）的日进水 SS 与 COD 记录，属于教学假设数据，用于演示计算与判读流程。

import numpy as np
import pandas as pd

# 教学示例数据：28天日进水SS(mg/L)与COD(mg/L)，假设值
np.random.seed(7)
days = pd.date_range("2024-03-04", periods=28, freq="D")   # 周一起始
ss  = np.array([168,175,182,190,205,158,150,   # 第1周
                172,180,188,196,210,162,155,   # 第2周
                160,166,171,178,186,150,146,   # 第3周
                240,255,268,290,310,205,198])  # 第4周：疑似冲击
cod = np.array([310,318,325,332,350,295,288,
                315,322,330,340,358,300,292,
                305,310,316,322,330,290,285,
                360,372,385,405,430,330,320])

df = pd.DataFrame({"date": days, "SS": ss, "COD": cod})
df["ratio"] = df["SS"] / df["COD"]
df["week"]  = df["date"].dt.isocalendar().week.astype(int)

def stats(g):
    m  = g["ratio"].mean()
    sd = g["ratio"].std(ddof=1)
    return pd.Series({"n": len(g), "mean": m, "std": sd, "CV": sd / m})

overall = stats(df)
by_week = df.groupby("week").apply(stats, include_groups=False)

print("全周期: mean=%.4f  std=%.4f  CV=%.4f" %
      (overall["mean"], overall["std"], overall["CV"]))
print(by_week.round(4))

# 判读：CV阈值分级
def grade(cv):
    if cv <= 0.15: return "稳定"
    if cv <= 0.25: return "需关注"
    return "不稳定"

print("\n全周期稳定性:", grade(overall["CV"]))
for w, r in by_week.iterrows():
    print("第%d周 CV=%.4f -> %s" % (w, r["CV"], grade(r["CV"])))

按上述示例数据运行，可得到全周期均值约 0.55、CV 约 0.13，整体判为"稳定"；但按周分组后，第 4 周因 SS 与 COD 同步抬升且 SS 抬升更快，CV 明显高于前三周。这说明全周期 CV 会掩盖局部冲击，必须按周分组比较，才能识别出"比值本身在漂移"的时段。

判读要点：若某周 CV 突增且均值同步上移，多为颗粒态负荷冲击（如管网冲刷、初沉池前溢流）；若 CV 突增但均值基本不变，多为采样或化验波动，应先核查 SS 烘干称量与 COD 消解环节，再决定是否触发工艺预警。
日进水监测 SS (mg/L) COD (mg/L) 连续28天 · 日值 比值与统计量 ratio = SS / COD 均值 mean 标准差 std CV = std / mean 按周分组比较 第1周 CV≈0.06 第2周 CV≈0.07 第3周 CV≈0.05 第4周 CV≈0.14 ↑ 均值同步上移 → 冲击 稳定性判读 CV ≤ 0.15 稳定 0.15～0.25 关注 > 0.25 不稳定 决定能否作预警 CV 分级色带（用于快速判读） 稳定 CV ≤ 0.15 需关注 0.15～0.25 不稳定 CV > 0.25 结论与使用边界 · 全周期 CV 稳定，但第4周局部 CV 抬升且均值上移，说明比值可作"周尺度"预警，不宜用全周期单一阈值。 · 建议以近4周滚动均值 ± 2×滚动标准差建立动态预警带，CV 超 0.15 时先核查采样与化验，再联动工艺。 · 若长期 CV > 0.25，说明进水来源复杂，SS/COD 不宜单独作可生化性预警，应结合氨氮、pH、BOD₅ 联合判断。

按周分组的关键实现细节

分组时有两个容易踩坑的点，直接影响 CV 的可比性：

  1. 周界划分：用 isocalendar().week 会跨年重置，跨年数据建议改用"起始日 + 7 天"的序号分组，避免第 52 周与第 1 周被错误合并。
  2. 样本量下限：每周至少 5 个有效日值再计算 CV，否则标准差本身不稳定。缺测日应先剔除，不要用插值补齐后再算 CV，否则会人为压低离散度。

若某周有效样本不足，代码中可加一行过滤：by_week = by_week[by_week["n"] >= 5]，并在报告中注明被剔除的周次，保证结论可追溯。

结果判读与预警指标可用性结论

观察现象可能原因处理建议
CV 稳定、均值平稳进水来源与管网运行正常可将该比值纳入日常预警基准带
CV 突增且均值上移颗粒态负荷冲击（冲刷、溢流）核查管网与初沉池，联动工艺调整
CV 突增但均值不变采样或化验波动先复核 SS 称量与 COD 消解，再判断
长期 CV > 0.25工业废水混入、来源复杂不宜单独作预警，需多指标联合

综上，本例示例数据下 SS/COD 比值在周尺度上具备较好的稳定性，可作为进水可生化性的辅助快速预警指标；但必须配合滚动统计与多指标交叉验证，不能以全周期单一阈值直接触发工艺动作。

← 上一章 📖 返回目录 下一章 →