摘要：近年来，珠海市的黑臭水体治理取得了一定成效，然而城中村合流制溢流污染问题成为城市水环境进一步改善的难点，迫切需要对其控制策略进行深入研究。为此，采用管道实时水质水量监测和暴雨洪水管理模型（SWMM）模拟相结合的方法，提出了通过理论截流倍数判别，低影响开发（LID）源头管控与截流-调蓄相结合的城中村合流制溢流污染控制策略，并精准确定相关设计参数，旨在为珠海市城中村合流制排水系统的改造和优化运行，以及黑臭水体治理提供理论支撑。
合流制溢流污染是导致珠海市城市水体黑臭的一个重要原因。目前，珠海市的合流制排水系统主要位于成建制的连片城中村。珠海市的农村60%以上以城中村的形式存在，且将与城市长期共存。加之，珠海市属于典型的南方滨海城市，雨季暴雨强度大，降雨场次多，合流制溢流量大且持续时间长，雨季合流制溢流污染问题尤为突出。虽然部分城中村在前期进行了分流制改造的规划设计，但由于施工难度大，且城中村的监管特性造成建成后难以维系完全分流制，故珠海市城中村分流制改造无法全面落地。因此，迫切需要对珠海市城中村合流制溢流污染控制策略展开研究。
笔者采用在线监测与SWMM相结合的方法，充分发挥实时监测数据的定量反演作用和对模型的率定验证功能，评估不同工程措施下珠海市典型城中村的合流制溢流污染负荷，以合流制排水系统年溢流污染负荷不高于分流制排水系统年雨水径流污染负荷作为溢流污染控制的最低标准，以就地化工程措施控制合流制溢流污染为原则，精准确定珠海市城中村合流制溢流污染控制策略的理论截流倍数、各种实际场景下的LID措施和截流-调蓄的相关设计参数，并以此为理论支撑，提出了城中村合流制排水系统的改造、优化运行和黑臭水体治理的具体措施。
1 研究方法
1.1 代表性区域的选取
目前，珠海市共有城中村108 个，主要集中在香洲区和斗门区，大多城中村为成建制且密集型分布。界涌村位于香洲区前山河流域上游的前山拱北105国道排洪渠流域北端，合流制排水系统服务面积约为0.23km2，是典型的成建制城中村，土地利用以混凝土道路、绿地和密集住宅为主，拆建难度大，合流制排水体系可能长期存在；而且其合流制排水系统从前山河流域分流制污水主干管起端接入，雨季溢流污染严重，导致受纳水体（105国道排洪渠）黑臭严重。为此，选取界涌村合流制排水系统为典型代表进行研究。该村原有合流制排水系统的汇水区域边界、合流制干管、市政分流制污水主干管、溢流口和溢流受纳水体的位置见图1。
根据资料收集和现场踏勘，界涌村合流制干管（DN500）埋设于界涌村排洪渠底，该干管末端接入珠海市105国道分流制市政污水主干管，进入前山河流域排水系统。界涌村内各合流制支管接入干管前均设有溢流井，雨季合流制溢流污水直接排入界涌村排洪渠，并最终汇入溢流受纳水体105国道排洪渠。
1.2 监测方法
本研究前期对界涌村合流制排水系统的流量及液位进行了一年的实时在线监测，并对水质进行了一年的人工采样检测，分别获取了合流制系统旱流流量和水质的逐时变化曲线，并作为基于SWMM的界涌村现状合流制模型的节点旱流输入。
为了获取更精确的逐分钟降雨数据，用于分析合流制溢流污染规律，本研究同步对界涌村2020年4月—2021年3月的降雨进行了实时在线监测。统计数据显示，在一年的监测期内，该村的年总降雨量为1 822.7mm，与珠海市九州港近15年的年均降雨量（1 704.3±214.5）mm相比，具有一定的代表性，亦可作为SWMM模型的年降雨数据输入。
为了对SWMM模型中透水区与不透水区的洼蓄量和曼宁糙率系数，入渗模型、污染物累积模型和污染物冲刷模型的相关参数进行准确率定与验证，避免合流制旱流污水的短时水质水量的波动性和不可预测性对此带来的干扰，提高率定精度，本研究选取界涌村一处已进行过分流制改造的地块进行SWMM模型的率定与验证，并对三场降雨条件下，该区域分流制雨水管道排口处的水质和水量进行了实时监测。
1.3 SWMM模型的建立、率定及验证
根据界涌村原有的合流制排水系统摸排图、规划的及部分已实施的分流制雨水管道系统设计图，以及区域用地性质图、区域高精度地形图，分别建立基于SWMM的原有合流制排水系统和规划中的完全分流制排水系统水质水量动态模型（简称：合流制模型和分流制模型），如图2所示。
本研究利用Morris法对模型参数进行了敏感性分析和率定，采用两场实际降雨事件（降雨量分别为51.5、10.5mm，重现期P分别为0.7、0.2a）率定区域雨水管排口处的水量及水质实测数据，对模型的参数进行了逐一率定。此外，选用另外一场实际降雨事件（降雨量为27.5mm，P=0.3a）的实测流量和水质数据，对率定后的模型进行了验证。
2 结果与讨论
2.1合流制理论截流倍数的取值
考虑到城中村合流制溢流污染控制工程的实施难度，本研究仅以合流制排水系统年溢流污染负荷不高于对应分流制排水系统年雨水径流污染负荷作为城中村溢流污染控制的最低标准（以下简称：合流制溢流污染最低控制标准），确定城中村合流制理论截流倍数n0（即不考虑下游处理系统的能力）。采用该理论截流倍数作为珠海市城中村合流制系统是否需要进行溢流污染控制的参考标准，当城中村合流制系统实际的截流倍数低于该理论截流倍数时，则必须要考虑采取源头控制或过程调蓄等就地化的溢流污染控制工程措施。
首先将2020年4月—2021年3月的逐分钟降雨数据输入分流制模型，模拟结果显示，研究区域年雨水径流污染负荷如下：SS负荷为21867.0kg/a，COD负荷为4423.7kg/a，TN负荷为885.9kg/a，NH3-N负荷为648.2kg/a，TP负荷为120.9kg/a。
对于合流制系统，其截流倍数决定了溢流污染总量。因此，本研究在合流制干管末端不同过流能力条件下（即不同理论截流倍数n0条件下），将2020年4月—2021年3月的实测分钟降雨数据输入合流制模型中，模拟计算年溢流次数和年溢流污染负荷，如表1所示。值得注意的是，当n0=0.2～4.8时，合流制系统的年溢流SS污染负荷均低于对应的分流制系统的年雨水径流SS污染负荷，这可能缘于合流制系统对高SS的初期雨水的有效收纳。
通过非线性拟合发现，截流倍数n0与各项污染物的年溢流污染负荷均满足二次函数关系，如图3所示。因此，可通过这些函数关系计算出不同n0对应的各项污染物年溢流污染负荷。由图3可知，当合流制系统的COD、TN、NH3-N、TP年溢流污染负荷分别等于分流制系统年雨水径流污染负荷时，n0分别为4.4、1.1、1.6和0.3。故按最不利原则，当界涌村合流制系统的n0≥4.4时，才能确保达到合流制溢流污染最低控制标准，即城中村合流制理论截流倍数为4.4，参考《室外排水设计标准》（GB 50014—2021）中的n0推荐值2～5，该取值在推荐范围内。但根据界涌村合流制干管末端的流量监测数据，雨季该合流制系统受下游市政主干管壅水的影响，实际截流倍数约为1.2，低于设计值2.0，且远小于4.4。为就地控制合流制溢流污染，不将溢流转嫁至系统下游，应当考虑采用LID源头管控和截流-调蓄相结合的措施对原合流制系统进行改造。
图3 截流倍数n0与各项污染物年溢流污染负荷的关系
2.2 LID控制策略
雨水径流的源头减量和净化是削减合流制系统雨天溢流污染的重要措施，故本研究提出通过海绵措施削减雨水径流，减轻降雨时合流制排水系统的压力，以有效控制溢流污染。研究借助合流制模型计算得到2020年4月—2021年3月界涌村的年径流总量控制率仅为38%，远未达到《珠海市海绵城市规划设计标准与导则（试行）》（修订版）中提出的
居住用地年径流总量控制率（改建、扩建项目）需达到60%的要求。因此，根据界涌村以居住用地为主、建筑物密度大、屋面和道路硬化、径流系数偏大的特点，本研究设置了植草沟、雨水花园及透水沥青路面三种可行的非入户LID措施（具体布置见图4），将年径流总量控制率控制在60%。
将2020年4月—2021年3月的降雨数据输入加载了LID措施的合流制模型中进行计算。当n0分别为2.0、1.2和0.2时，采用LID措施后合流制溢流污染的削减情况见表2。采用LID措施后，年径流总量控制率达到了60%，年溢流水量和污染负荷都有明显的削减，年溢流水量削减率均达到了60%左右，年溢流次数削减了21%～33%。当截流倍数为0.2时，各污染物的年溢流总量削减率在45%～55%之间，当截流倍数为1.2和2.0时，各污染物的年溢流总量削减率在55%～60%之间。因此，采用LID措施后，即使合流制系统雨季排水严重受阻，截流倍数低至0.2，仍能满足合流制溢流污染控制最低标准。
2.3 调蓄池控制策略
除采取源头控制的LID措施以外，修建调蓄池也是合流制系统控制溢流污染最常用的措施之一。调蓄设施可以用于削减峰值流量、控制地表径流污染，并且当雨水水质较差时，调蓄池能对雨水径流进行高效收集和分质截流。借助合流制模型模拟计算发现：当截流倍数n0>1.6时，除COD外，其余污染物的年溢流污染负荷均低于对应的分流制系统雨水径流污染负荷。因此，调蓄池的理论容积只需满足合流制年溢流COD负荷≤分流制年雨水径流COD负荷。本研究进一步借助合流制模型计算出不同截流倍数和不同调蓄池理论容积下的年溢流COD负荷。通过非线性拟合发现，调蓄池理论容积与年溢流COD负荷之间呈显著指数函数关系（R2>0.999），如图5所示。通过该函数关系，可以准确计算出不同截流倍数条件下调蓄池的理论容积，如表3所示。根据《室外排水设计标准》（GB 50014—2021）中推荐的合流制溢流调蓄池调蓄量的计算方法，可以根据调蓄池容积计算调蓄池建成后的截流倍数。
调蓄池的最大排空周期与相邻两场降雨的间隔时间有关。对珠海市2020年5月—9月（雨季）112场降雨的相邻间隔时间进行累积频率分析，如图6所示。可知，雨季珠海市相邻两场降雨的平均间隔时间为33h，50%和70%累计频率对应的间隔时间分别为14h和24h。
综上，从设计的角度进一步考察了调蓄池最大排空周期在12～48 h之间变化时对调蓄池容积的影响，结果如表4所示。可知，在上述区间范围内，调蓄池容积的变化系数为1.3，可作为选择安全系数时的参考。虽然前山河流域污水处理系统旱季处理能力的冗余量与调蓄池的运行密切相关，但鉴于整个排水系统的上下游各子系统调蓄池的复杂性，难以量化其对局部单一调蓄池运行的影响。故对调蓄池的实际运行而言，可依据上游泵站站前液位或污水厂厂前液位的实时反馈来控制。
3 结论与建议
①以合流制系统年溢流污染负荷等于分流制系统年雨水径流污染负荷作为溢流污染控制的最低标准，从污染负荷角度分析，珠海市城中村合流制系统的实际截流倍数不应低于4.4，当低于此值时，则应采用溢流污染控制的工程措施。
②根据溢流污染最低控制标准，在城中村合流制干管降雨时排水严重受阻的情况下，采用LID措施进行合流制溢流污染控制，需将年径流总量控制在60%以上。
③根据溢流污染最低控制标准，调蓄径流深度与合流制系统实际截流倍数满足二次函数关系。故不同实际截流倍数下，设计调蓄径流深度可按该函数关系进行取值，并可按规范考虑1.1～1.5的安全系数以计算调蓄池容积。
④为了提高整个城市污水系统的雨季处理能力的可操作性，以及降低上下游子系统与之匹配运行的控制难度，对汇入主干管的各子系统（如城中村合流制排水系统）应进行科学合理地限流，分散控制溢流污染，避免降雨时下游子流域混合污水无限制、无序地涌入排水主干管，占据上游子系统正常的排水下行通道，导致上游子系统雨季溢流污染加重。
