热干旱地区市政污水处理厂进水变异性的概率化设计包络线构建——以科威特Kabd污水处理厂八年运行记录为例

《Journal of Engineering Research》：Influent Variability and Probabilistic Design Envelopes for Municipal Wastewater in Hot Arid Regions

　　

市政污水处理厂进水特征因人均用水量、气候、工业贡献及下水道系统入渗等差异而呈现显著的区域性变异。在热干旱地区，高温、高蒸发率及社会经济模式进一步加剧进水不确定性，常导致处理设计与实际运行条件不匹配。海湾阿拉伯国家合作委员会（GCC）国家合计年产废水约4×109 m3，人均产率100–200 m3/（人·年），由高用水量（400–550 L/（d·人））驱动。尽管基础设施投入可观，有效的污水处理厂（WWTP）设计仍需可靠的进水表征，尤其在水温常超过30°C、加速生物动力学但同时对氧溶解度和硝化菌生物量稳定性构成挑战的热干旱气候下。国内废水进水表征构成WWTP设计基础，指导水力能力、生物反应器尺寸、曝气需求及污泥处理系统。然而，由入流特征、商业活动及工业贡献驱动的进水变异性仍是WWTP性能不确定性的主要来源，因而需要概率化设计方法而非带有任意安全系数的确定性估算。已发表的中东北非（MENA）研究提供了若干市政污水处理厂的实测进水特征，但多侧重均值、范围、季节变化或跨厂比较；能够联合量化完整水力记录、浓度与负荷分布、上尾量级、日历日持续性、季节温度及源匹配国际比较的长期分析仍然有限。本研究以Kabd污水处理厂为科威特案例，目标为：（i）在一致界定的研究期内表征长期进水流量、有机物、固体、总凯氏氮（TKN）、总磷酸盐（T-PO4）及温度；（ii）量化经验P50、P90和P95浓度与负荷包络线，以及上尾流量和负荷事件的重现性与持续时间；（iii）采用匹配分析物、源感知描述符和透明相对偏差指标，将Kabd观测值与MENA及其他气候背景下实测市政进水数据集进行比较。其新颖贡献在于在长期记录内联合评估负荷量级、超限持续性、季节温度及分析物匹配的国际比较，为场地特异性进水设计包络线提供集成依据。

研究背景与问题　市政污水处理厂（WWTP）进水特征受人均用水量、气候、工业贡献及下水道入渗等因素影响而存在广泛区域差异，进而造成水力工况、有机负荷剖面和营养物浓度的分化，给基于模型的设计带来不确定性。在热干旱地区，高温、高蒸发率和社会经济模式进一步放大进水不确定性，常使处理设计偏离实际运行条件，导致生物反应器、曝气系统和脱氮除磷单元在全球范围内出现尺寸过大或不足。海湾阿拉伯国家合作委员会（GCC）地区水资源稀缺，使废水作为非常规资源具有战略意义，但系统性长期进水表征仍十分有限。已发表的中东北非（MENA）研究多侧重均值、范围、季节变化或跨厂比较，尚缺乏在长期记录中联合量化完整水力记录、浓度与负荷分布、上尾量级、日历日持续性、季节温度及源匹配国际比较的分析。为此，研究人员以科威特Kabd污水处理厂为案例开展研究。

研究概况与结论　研究人员基于2015年1月1日至2023年3月31日的连续日流量记录（覆盖全部3,012个日历日）及约每周6天的进水抓样水质数据（2,474个日期；TKN与T-PO₄分别为830和829个日期），构建了经验P50、P90、P95浓度与同日负荷包络线，分析了上尾流量与负荷事件的重现性和持续时间，并开展季节温度表征及源感知的国际基准比较。结论表明：Kabd设计相关进水条件应视为量级与持续时间的分布而非单一流量或浓度；流量变异系数（CV）为0.089，P90和P95分别达名义平均设计流量的99.6%和101.5%，近半数超限事件持续至少2天，而实测最大值仍低于报告峰值水力能力；P95的化学需氧量（COD）、五日生化需氧量（BOD₅）和总悬浮固体（TSS）负荷较中位数高29.3%–33.4%，持续性可区分持续有机需氧与较间歇的固体负荷；季节流量与浓度变化并不同步，夏季中位负荷仅略低于其他季节，温度则定义了冬季（保守生物速率与硝化校核）和夏季（氧传递与溶解氧管理）两个不同工艺边界；国际比较不支持Kabd废水普遍更强或更弱的简单判断，其COD与若干MENA及国际均值接近，但在BOD₅、TSS和TKN上位置发生变化，气候标签或单一COD系数不足以跨流域转移设计假设。其重要意义在于提出从确定性进水规格转向不确定性感知、本地推导的设计基础，数值包络线限于Kabd本地，但工作流程可迁移至其他WWTP以构建区域包络线。

关键技术方法　研究采用单一污水处理厂长期运行档案的回顾性分析设计，样本来源为科威特公共工程部（MPW）Kabd污水处理厂监测记录：日流量取自SCADA系统连续日聚合数据，水质为原进水单次抓样，按APHA标准方法测定温度、可沉降固体、TSS、挥发性悬浮固体（VSS）、COD、BOD₅、油和脂、TKN及T-PO₄。数据分析采用经验分位数（P50、P90、P95及最大值）构建浓度与同日负荷包络线，以经验超限曲线描述超限概率；以连续日历日“游程”定义超限事件，统计事件数、单日与多日事件及持续时间中位数、P75和最大值，并对孤立单日采样缺口进行桥接敏感性分析；按气候季节（DJF、MAM、JJA、SON）开展流量、浓度、负荷与温度的分布分析；国际比较采用源感知证据矩阵，保留各来源统计量、采样支撑、流域背景与分析物身份，并以相对偏差指标Δ=100[（Kabd中心统计量/C_j,a）?1]及Kabd观测超过外部中心值的比例进行描述性比较，未开展跨研究推断性检验。

研究结果

进水水力模式　通过对3,012天日流量记录统计得出，平均与中位流量分别为162,349和162,839 m³/d，CV为0.089，低于若干受降雨影响的城市系统；偏度为?1.337，反映偶发偏低流量。P75、P90、P95分别为172,894、179,202和182,643 m³/d，处于围绕中心工况的较窄区间。上尾分析表明，P90对应302个超限日、157次事件，其中74次（47.1%）持续至少2天、6次至少5天，最长16天；P95对应151个超限日、86次事件，39次（45.3%）持续至少2天，最长7天。两个阈值均未接近270,000 m³/d的峰值水力能力。季节中位流量秋季最高（167,518 m³/d）、夏季最低（155,088 m³/d），相差7.4%；年均流量2016年最低（149,183 m³/d）、2021年最高（174,092 m³/d），2022年降至166,101 m³/d；名义平均设计流量在2016–2017年未被超过，2020年、2021年和2022年分别有31.7%、23.8%和3.3%的观测日超过，说明合并百分位数汇总了研究期内不同运行水平。

进水有机物与固体负荷模式　通过浓度与同日负荷分位数统计得出，COD和BOD₅浓度中位数分别为569和296 mg/L，TSS中位浓度为168 mg/L；其P95同日负荷分别较中位数高31.7%、29.3%和33.4%。油和脂中位浓度39 mg/L、P95为68 mg/L。低频营养物序列亦具明显上尾：TKN由中位51.7 mg N/L升至P95的81.6 mg N/L，P95负荷14,347 kg N/d；T-PO₄由17.0 mg PO₄/L升至24.9 mg PO₄/L，P95负荷4,240 kg PO₄/d。同日配对的BOD₅/COD和TSS/COD中位比值分别为0.522和0.297。负荷–持续时间曲线在最低超限概率处急剧上升，表明少数观测远超经常性P95包络线，这是同时报告百分位数与最大值的实践理由。持续性分析显示，P90下持续至少2个连续观测日的事件占COD事件的35.2%、BOD₅事件的50.0%、TSS事件的21.2%；P95下相应比例为26.3%、39.7%和19.3%。最大事件持续时间对采样缺口敏感：严格日历最大值为6 d，而单日桥接假设下最大跨度达11–66 d，其中最长者为COD的P90事件，覆盖2020年3月19日至5月23日，含55个观测超限日和11个桥接未采样日；桥接情形下各参数–阈值对的中位跨度均为1 d。

负荷与浓度的季节调制　通过季节分布比较得出，浓度与负荷的季节变化并不同步。夏季中位负荷最低：COD 91,358 kg/d、BOD₅ 46,870 kg/d、TSS 27,196 kg/d，较其他季节分别低1.3%–7.7%、3.0%–8.0%和0.2%–4.9%，季节差异总体温和，TSS夏季与秋季中位数几乎相同。各成分浓度并非普遍在夏季最低，该模式与夏季流量较低一致，提示流量或居住人口相关影响而非均匀稀释。该组合支持在设计与运行裕度校核中采用季节特异的进水负荷与温度边界。

进水温度表征　通过2,474个观测的温度分布统计得出，温度范围12.3–37.7°C，总体中位数28.9°C；中位数由冬季22.9°C升至夏季34.0°C，季节差11.1°C，春秋为过渡工况。两个季节极端具有不同工程含义：冬季包络线为温度敏感的生物速率与硝化校核提供更保守输入，夏季中位接近34°C则标示溶解氧控制需重点复核的时期。实测记录由此界定了冬季生物动力学边界与夏季氧管理边界。

源感知基准比较与场地特异含义　通过纳入12个同行评议的实测未处理进水原始数据集（9个来自MENA国家，3个来自巴西、希腊及波黑）的证据矩阵比较得出，Kabd平均COD处于实测文献的中间位置，与阿曼Al-Ansab、阿联酋Al Ain、埃及Sharkia 21厂聚合及希腊Evros的COD均值相差在6.3%以内；但COD相近并不意味着组成相近，相对Al-Ansab，Kabd的COD低3.5%而BOD₅和TSS分别高32.6%和25.2%；相对Evros，COD低6.3%而BOD₅和TSS分别低16.2%和44.5%。阈值分析显示，仅26.2%的Kabd COD观测超过Al-Ansab均值，而96.7%和93.5%分别超过其BOD₅和TSS均值；Kabd COD超过Al Ain、Sharkia和Evros均值的观测比例分别为42.6%、45.8%和19.8%，而99.7%超过Doha-South历史COD均值，仅0.9%和0.2%超过较高的M'zar和Fuhais均值。营养物比较因分析物身份而范围更窄：Kabd平均TKN 53.1 mg N/L较Al Ain的29.5 mg N/L高80.1%，较M'zar的142.2 mg N/L低62.6%；其他来源报告的TN、氨氮、NH₄和NNH₃不能替代TKN，TP或不同定义磷酸盐亦不能替代T-PO₄。中位同日VSS/TSS比为0.42，对应固定固体分数约58%，提示悬浮固体中无机贡献显著。上述结果表明表7识别的是不同用水、流域、下水道与采样条件下形成的多种强度与组成组合，而非单一热干旱进水特征；设计验证应对COD、BOD₅和TSS分别采用浓度与负荷包络线而非相互推断。

讨论与结论　讨论部分指出，Kabd的数值包络线为单一厂记录且计划性水质缺口主要出现在周五，故其数值仍属场地特异，主要成分游程长度代表最低确认持续时间，在提出科威特全境或区域数值包络线前需开展多厂验证。比较中较高浓度可源于有限的生活污水稀释、工业或商业负荷、流域增长或其组合，较低浓度则可能反映水力稀释、流域条件或测量支撑差异，因此气候标签或单一COD系数不足以跨流域转移设计假设。可辩护的设计用法是将典型运行锚定于本地P50，以P90和P95校核经常性与较低频应力，加入事件持续时间以区分孤立与持续负荷，并将这些负荷包络线与冬季和夏季温度边界耦合。研究结论为：Kabd设计相关进水条件是一个量级与持续时间的分布而非单一流量或浓度；合并流量记录CV为0.089，P90和P95分别达名义平均设计流量的99.6%和101.5%，相应超限游程近半数持续至少2天，该组合区分了经常性近能力运行与较低频应力，而实测最大值仍低于报告峰值水力能力。仅凭浓度同样不能反映工艺负荷：P95的COD、BOD₅和TSS负荷较中位数高29.3%–33.4%，持续性可区分持续有机需氧与较间歇的固体负荷模式。流量与浓度的季节变化并不同步，夏季中位负荷仅略低于其他季节，而温度界定了两个不同工艺边界：冬季用于保守生物速率与硝化校核，夏季用于氧传递与溶解氧管理。因此水力、质量负荷、持续性与热包络线需联合评估。国际比较不支持Kabd废水普遍更强或更弱的简单判断，其COD与若干MENA及国际均值接近，但考虑BOD₅、TSS和TKN后位置发生变化；部分厂较高值关联于有限生活污水稀释或有记录的工业与合流制输入，较低值则出现在被描述为稀释或采样支撑有限的数据集中。故气候标签或单一COD系数不足以跨流域转移设计假设，Kabd数值包络线在经其他科威特与MENA厂验证前应保持场地特异。实践结论是从确定性进水规格转向不确定性感知、本地推导的设计基础：以P50作为典型负荷参照，P90和P95作为两个上尾条件，以事件持续时间描述持续性，以配对浓度与质量负荷分布指导工艺尺寸确定，以季节温度限值用于生物与曝气校核；该框架揭示单一引入设计值可能掩盖水力、需氧量、澄清、固体处理与营养物处理的不同风险。数值限于Kabd本地，但工作流程可迁移，在其他WWTP一致应用将可基于可比证据而非假定相似性构建未来区域包络线。