我国污水厂进水普遍具有低C/N的水质特征，限制了污水厂的脱氮除磷效率^[1-2]。为满足出水氮磷的达标排放，通常选择在污水处理流程中投加有机碳源^[3-4]。虽然提升了氮磷脱除效果，但存在碳源药剂费用高昂、污泥产量增加、后续处理处置费用高、大部分有机物转化为CO₂ (不符合碳达峰、碳中和的时代发展需求)^[5]等诸多问题。因此，急需找寻新的方法或策略以摆脱低C/N污水处理对外加碳源的依赖。

污水处理工艺中好氧区会消耗大量有机物，造成碳源的浪费，控制溶解氧能够节约碳源，提高碳源在缺氧段的脱氮利用率。近几年，较多研究提出通过控制污水处理工艺中好氧区溶解氧(DO)可实现低C/N污水的高效脱氮除磷，无需外加碳源的同时可节省曝气能耗。JIN等在低C/N进水条件下，通过2步控制DO<0.5 mg·L⁻¹，实现87.9%的TN去除率^[6]。赵智超等在AAO反应器中通过间歇曝气降低好氧区DO浓度，在无硝化液内回流条件下实现了同步硝化反硝化脱氮除磷，出水TN和TP分别为5.8 mg·L⁻¹和0.3 mg·L^−1[7]。此外，低DO控制过程中微生物群落与功能基因的变化，也是实现低C/N进水条件下各工艺高效脱氮除磷的关键原因之一^[8-9]。池玉蕾等^[10]的研究表明，对于低C/N污水而言，好氧区DO的降低能显著提高系统内反硝化菌和聚磷菌的丰度，从而提高污水处理系统的脱氮除磷效率。

然而，在低C/N进水条件和低DO控制策略下污水处理效能及微生物群落结构的研究大多停留在小试及中试阶段，缺乏针对实际工程的应用研究。对于城镇污水处理厂而言，DO的变化会对整个工艺系统产生较大影响，其准确调控至关重要。因此，有必要针对实际工程研究低DO调控策略下污水处理系统的综合效能。本研究选择山东某城镇污水处理厂(低C/N进水)进行实验，比较在无外加碳源并调控好氧区DO措施实施前后污水厂脱氮除磷效能及微生物种群结构的变化，并综合分析低DO控制策略产生的经济效益。面对碳达峰、碳中和的时代发展目标，本研究通过低DO调控，保障低C/N进水污水厂出水达标的同时，实现药耗(无需外加碳源)和能耗(降低曝气)的同步节省，不仅大大节约了运行费用，也为同类型实际工程的调控及改造提供参考。

1.   材料与方法

1.1   工程项目概况

本研究依托于山东省某城镇污水处理厂，该厂设计处理规模为8×10⁴ t·d⁻¹，实际运行规模8×10⁴ t·d⁻¹，采用A²O工艺。进水C/N平均值为4.2:1，属典型低C/N进水水质，常年大量投加有机碳源，出水水质执行《城镇污水处理厂污染物排放标准》(GB 18918-2002)一级A标准。主要工艺流程见图1，主要构筑物参数见表1。

图 1  污水厂主要工艺流程

Figure 1.  Main process flow of sewage treatment plant

下载: 全尺寸图片 幻灯片

表 1  工艺构筑物参数

Table 1.  Parameters of process structures

构筑物	HRT/h	池体数/个	容积/m3	水深/m
厌氧区	3.6	1	12 000	4.3
缺氧区	2.5	2	8 332	6.0
好氧区	5.0	2	16 664	6.0
二沉池	6.0	4	20 400	4.5

 | Show Table

DownLoad: CSV

1.2   实验方法

本研究的运行周期分为2个阶段，第1阶段(P1)为2020-01-01—2020-12-31，共365 d，每天根据出水TN浓度投加碳源(乙酸钠)；第2阶段(P2)为2021-01-01—2021-05-22，共143 d。其中P2阶段开始逐渐减少碳源投加量，在416 d后停止投加碳源。同时，P2阶段硝化液回流比分4个阶段下调：366~399 d维持在370%；400~425 d逐渐下调至280%；426~455 d逐渐下调至180%；456~484 d逐渐下调至90%，之后不再改变。2个阶段运行参数见表2。

表 2  污水厂2个阶段运行参数

Table 2.  Operating parameters of sewage treatment plant in two phases

阶段	乙酸钠投加量/(mg·L−1)	DO/(mg·L−1)	硝化液回流比/%	污泥回流比/%
P1	73~575	2~4	370	100
P2	0~284	0.5~1.4	370~90	100

 | Show Table

DownLoad: CSV

1.3   水质指标及微生物检测分析方法

1)水质指标测定与分析。水样COD、NH₄⁺-N、NO₃⁻-N、NO₂⁻-N、TN、TP、MLSS等指标均根据国家标准方法测定^[11]。水样DO浓度通过便捷式溶解氧仪(YSI-pro20)测量。

2)微生物样品采集与保存。基于表观实验数据，分别于P1(第140天)和P2(第499天)2个阶段从污水厂采集2组活性污泥样品。采样位置为奥贝尔氧化沟好氧区的前段、中段和后段，每段于液面下1、3和5 m处采集3个样品，将9个样品混合形成1组微生物样品。采集后样品保存于−20 ^oC，用于后续微生物测序。

3)宏基因测序与分析。利用E.Z.N.A.® Soil DNA Kit(Omega Bio-tek，美国)试剂盒对活性污泥微生物样品进行DNA抽提；使用NEXTflex^TM Rapid DNA-Seq(Bioo Scientific，美国)建库；使用Illumina NovaSeq测序平台进行宏基因组测序；使用上海美吉生物医药科技有限公司云平台进行数据分析。

1.4   计算方法

好氧区硝化比例、好氧区亚硝化比例、好氧区同步脱氮比例按式(1)、式(2)、式(3)计算。

式中：R₁为好氧区硝化比例，%；R₂为好氧区亚硝化比例，%；R₃为好氧区同步脱氮比例，%；${C_{{\rm{NH}}_4^ + - {\rm{N}},{\rm{in}}}}$,${C_{{\rm{NH}}_4^ + - {\rm{N}},{\text{ef}}}}$分别为好氧区进出水NH₄⁺-N质量浓度，mg·L⁻¹；${C_{{\rm{NO}}_3^ - - {\rm{N}},{\rm{in}}}}$,${C_{{\rm{NO}}_3^ - - {\rm{N}},{\text{ef}}}}$分别为好氧区进出水NO₃⁻-N质量浓度，mg·L⁻¹；${C_{{\rm{NO}}_2^ - - {\rm{N}},{\rm{in}}}}$,${C_{{\rm{NO}}_2^ - - {\rm{N}},{\rm{ef}}}}$分别为好氧区进水NO₂⁻-N质量浓度，mg·L⁻¹。

2.   结果与讨论

2.1   运行效能

针对P1、P2 2个阶段，考察系统在调控DO前后对有机物、氮、磷等污染物的转化去除效果。实验过程中污水厂进出水COD变化情况见图2。P1和P2阶段污水厂出水COD值无显著差异。P1阶段实际进水COD平均值为215 mg·L⁻¹(C/N=4.4∶1)，外加碳源COD当量平均值为41 mg·L⁻¹(总C/N=5.2∶1)；P2阶段实际进水COD平均值为244 mg·L⁻¹(C/N=4.2∶1)。在P2阶段降低好氧池DO浓度后，污水厂出水COD值基本维持在25 mg·L⁻¹以下，490 d后出水COD略有上升，但仍然满足排放要求。COD平均去除率由P1阶段的91.5%变为P2阶段的92.5%，变化不明显，可见DO的降低对有机物的去除几乎无影响。

图 2  进出水COD、乙酸钠投加量及COD去除率变化

Figure 2.  Changes of influent and effluent COD, sodium acetate dosage and COD removal ratio

下载: 全尺寸图片 幻灯片

2个阶段的NH₄⁺-N、TN去除及NO₂⁻-N积累情况见图3。P1和P2阶段对NH₄⁺-N的去除效果无显著差异。2个阶段污水厂进水NH₄⁺-N平均质量浓度分别为46.6 mg·L⁻¹和49.6 mg·L⁻¹，降低好氧区DO浓度后，P1和P2阶段污水厂出水NH₄⁺-N质量浓度基本都维持在1 mg·L⁻¹以下，平均去除率在97%以上。可见DO浓度的降低对NH₄⁺-N的去除效果无显著影响。

图 3  进出水氨氮、总氮、亚硝氮及碳氮比、氨氮去除率、总氮去除率变化

Figure 3.  Changes of influent and effluent NH₄⁺-N, TN, NO₂^--N and C/N, removal rate of NH₄⁺-N, TN

下载: 全尺寸图片 幻灯片

P2阶段污水厂对总氮的去除效果明显好于P1阶段。P1阶段污水厂进水和出水TN平均质量浓度分别为48.9 mg·L⁻¹和 11.6 mg·L⁻¹，平均去除率为76.3%；P2阶段污水厂进出水TN平均质量浓度分别为57.9 mg·L⁻¹和 9.9 mg·L⁻¹，平均去除率升高至82.9%。P2阶段，本研究在控制DO浓度、降低外加碳源的同时，硝化液回流比从P1阶段的370%逐渐下调，最终调整为P2阶段的90%，再加上100%的污泥回流比，TN的理论去除率应从82.5%降至65.5%。P1阶段的TN平均去除率低于理论去除率；然而，P2阶段的TN去除率却明显高于理论去除率。控制溶解氧浓度及减少碳源投加控制策略下，污水处理厂的脱氮能力非但没有降低，反而显著提升。

P2阶段开始对污水厂出水NO₂⁻-N浓度进行跟踪检测，结果如图3所示。长期低DO运行条件下，出水TN逐渐以NO₂⁻-N为主。第367~406天，出水NO₂⁻-N质量浓度几乎为0，之后逐渐升高，在508 d时达到最高值8.5 mg·L⁻¹，有明显的NO₂⁻-N积累。这与姚丽婷等在小试实验中的研究结果一致。在低DO条件下，会出现显著的NO₂⁻-N积累现象，主要原因为亚硝酸盐氧化菌(NOB)受到抑制，氨氧化菌(AOB)逐渐成为优势菌，最终实现短程硝化^[12]。P2阶段开始逐渐停止外加碳源，C/N从5.2∶1(含外加碳源)降至4.2∶1，而TN去除效果却明显有提升。其原因之一可能是降低DO后，好氧区发生短程硝化，而反硝化去除NO₂⁻-N所需电子供体仅为NO₃⁻-N的60%，大大降低了污水厂脱氮对碳源的需求^[13]，也是保障低碳源条件下脱氮效能的关键因素之一。

污水厂2个阶段TP去除效果见图4。P2阶段污水厂对TP的去除效果优于P1阶段。降低好氧区DO浓度后，TP的去除效果有明显提升，沉淀池出水平均质量浓度从0.64 mg·L⁻¹降至0.34 mg·L⁻¹，去除率由87.0%上升至92.7%，说明DO浓度的降低提升了污水处理系统的除磷效果。池玉蕾等在小试实验研究中得出同样的结果，认为低DO条件促进了系统中除磷功能微生物的生长，从而表现出更好的除磷效能^[10]。

图 4  进出水总磷及去除率变化

Figure 4.  Changes of influent and effluent TP and removal ratio

下载: 全尺寸图片 幻灯片

为探究低DO条件下污水厂脱氮能力提升的原因，分别于P2阶段的第392、475、493、499天监测各生化单元的出水水质情况，结果如图5所示。COD、NH₄⁺-N与TP在各池体单元的去除规律无显著变化。其中，在厌氧池中的COD显著下降。分析主要原因为非溶解性有机物在厌氧池被活性污泥快速吸附及污泥回流(100%)的稀释作用^[14-15]。

图 5  各工艺阶段COD、氮、磷质量浓度

Figure 5.  COD, N, P mass concentration at each process stages

下载: 全尺寸图片 幻灯片

P2阶段降低DO后，短程硝化作用逐渐显现。第392、475、493、499天，NH₄⁺-N在好氧池的转化量基本无差异；然而，随着时间的增加，好氧池末端出水NO₃⁻-N质量浓度逐渐降低，分别为11.1、7.8、1.7和0.1 mg·L⁻¹；好氧池末端出水NO₂⁻-N质量浓度逐渐升高，分别为0、2.6、5.0、5.9 mg·L⁻¹。长期低DO运行条件下，好氧池对NH₄⁺-N的转化逐渐从全程硝化变为短程硝化。此外，第475、493、499天，好氧池出现了明显的TN脱除现象。

P2阶段低DO运行条件下，针对好氧区的氮转化平衡问题进行分析，结果如图6所示。P2阶段好氧区全程硝化逐渐转化为短程硝化，且表现出明显的脱氮效果。392 d时，好氧区主要发生全程硝化，之后好氧区全程硝化比例不断下降，由392 d的94.1%降至499 d的0.7%，而短程硝化的比例从0增加到38.0%。然而，值得注意的是，在长期低DO运行条件下，除了全程硝化和短程硝化外，仍有部分NH₄⁺-N既未转化为NO₃⁻-N也未转化为NO₂⁻-N。该部分NH₄⁺-N的比例由392 d的0逐渐上升至499 d的61.3%。基于好氧区的低DO及NO₂⁻-N持续积累状态，推测本研究在好氧区发生了同步脱氮作用，可能为同步硝化(也可能短程硝化)反硝化或者厌氧氨氧化作用。

图 6  好氧区氨氮硝化率、亚硝化率、同步脱氮率的变化

Figure 6.  Variation of NH₄⁺-N nitrification ratio, nitrosation ratio and simultaneous denitrification ratio in aerobic zone

下载: 全尺寸图片 幻灯片

2.2   微生物种群结构分析

为进一步探索P1和P2阶段污水厂表观差异的内部机制，对2个阶段污水厂好氧区的微生物进行宏基因组学测序并对比分析。微生物在门水平的群落结构如图7(a)所示。长期低DO运行条件造成污水处理系统内微生物在门水平上差异性显著。P1阶段系统内的优势菌为Proteobacteria(36.3%)，其次为Chloroflexi(20.0%)与Actinobacteria(17.7%)；P2阶段系统中丰度最高的仍为Proteobacteria(40.8%)，而Chloroflexi(15.1%)明显下降Actinobacteria(23.0%)显著升高。此外，值得注意的是，Nitrospirae(亚硝酸盐氧化相关菌门)的丰度由P1阶段的2.2%下降至P2阶段的0.6%。

图 7  微生物相对丰度

Figure 7.  Relative abundance of microorganism

下载: 全尺寸图片 幻灯片

微生物在属水平的丰度如图7(b)所示。P1和P2阶段系统内脱氮除磷主要功能微生物在属水平的差异性同样显著。硝化作用相关功能微生物变化为：氨氧化菌(AOB)Nitrosomonas和Nitrosomonadales的丰度分别由P1阶段的0.41%、0.32%变化至P2阶段的0.45%、0.21%，总丰度略有下降；而亚硝酸盐氧化菌(NOB)Nitrospirae的丰度却从P1阶段的0.15%大幅下降至P2阶段的0.06%^[16]。由此可见，低DO运行条件对NOB的明显抑制作用，可能是造成好氧区NO₂⁻-N积累的主要因素。

长期低DO运行条件下，P2阶段系统内反硝化功能微生物Hyphomicrobium、Lautropia、Dechloromonas、Acidovorax、Thauera、Comamonas、Ignavibacterium总体较P1阶段有所增加^[17-19]，相对丰度为原来的1.4倍。此外，P2阶段系统中除磷菌(PAOs)Candidatus Accumulibacter、Rhodobacteraceae和聚糖菌(GAOs)Defluviicoccus、Candidatus Competibacter均较P1阶段有显著提升，总的相对丰度为P1阶段的1.7倍。有研究表明，这4种功能菌能够利用厌氧期储存的内碳源进行反硝化脱氮^[20-21]。由此可见，低DO条件下系统内各反硝化功能微生物的增加同样是保障低碳源脱氮的关键。

长期低DO运行条件下，好氧区氮平衡的分析结果表明，好氧区有明显的氨氮(总氮)缺失，推测其原因可能为好氧区发生同步硝化反硝化或者厌氧氨氧化。ZHANG等的研究表明，在低DO条件下部分GAOs、PAOs在厌氧段储存的内碳源未被完全消耗，微氧条件下仍然能够进行反硝化脱氮^[22]。而对于微生物的分析中，未发现有厌氧氨氧化相关的功能微生物，认为未发生厌氧氨氧化作用。因此，本研究认为长期低DO运行条件促使具备反硝化功能的PAOs与GAOs丰度升高，在好氧区形成了基于内碳源的同步反硝化(包括亚硝氮反硝化)作用，从而在无外加碳源及降低内回流的条件下提高了污水厂对低碳源污水的脱氮效能。

此外，上节研究结果表明在长期低DO运行条件下，系统除磷效果有明显提升。基于此，对系统内除磷相关功能微生物进行分析。结果表明，低DO运行条件下，除具备反硝化除磷功能的Candidatus Accumulibacter、Rhodobacteraceae(能利用DO或NO₃⁻作为电子受体的PAOs)外，系统内另外一种除磷菌Tetrasphaera(仅能利用DO作为电子受体的PAOs)的丰度也有显著提升，P2阶段3种菌属在系统中的总丰度为P1阶段的2.3倍。分析其原因可能为，低DO运行条件下，提高了有机物作为内碳源用于生物除磷的利用率，促进了除磷相关功能微生物的生长，从而提升了系统的除磷效能^[10]。

基于功能微生物数据，对P1、P2 2个阶段污水处理系统内氮代谢相关功能基因进行分析，结果如图8所示。2个阶段微生物脱氮相关的主要功能基因丰度有显著性差异。P2阶段降低DO后，氨氧化相关功能基因AmoCAB总丰度几乎未变，亚硝酸盐氧化还原酶功能基因NxrAB和硝酸盐还原酶功能基因NarGHI/NapAB总丰度分别下降了6.7%和4.2%，亚硝酸盐还原酶功能基因NirKS总丰度升高了13.5%^[23-24]。结合表观效能及功能微生物丰度数据，推测系统中存在的一个内部机制为：长期低DO运行条件可抑制NOB，导致Nxr基因相对丰度下降，造成NO₃⁻-N生成量减少、NO₂⁻-N积累；而长期低NO₃⁻-N、高NO₂⁻-N条件下，Nar/Nap基因相对丰度降低，Nir基因相对丰度升高，也是保障无外加碳源条件下低C/N进水城镇污水厂脱氮效能的关键因素之一。

图 8  部分脱氮基因相对丰度

Figure 8.  Relative abundance of partial denitrification genes

下载: 全尺寸图片 幻灯片

2.3   经济效益分析

长期低DO运行条件下，为污水厂带来的经济效益主要体现在3个方面：1)节省外加碳源，药剂费用大大降低；2)曝气量降低，电耗费用减少；3)低DO、低/无外加碳源条件下，污泥产量降低，污泥处理处置成本降低。污水采用低DO运行策略所带来的主要经济效益计算见表3，其中乙酸钠(含量20%)750 元·t⁻¹，电费0.52 元·(kWh)⁻¹，污泥处理费(80%含水率)270 元·t⁻¹，P2阶段的污水厂运行费用比P1阶段降低2.44 万元·d⁻¹，吨水处理费用减少0.3元。

表 3  污水厂经济效益表

Table 3.  The economic benefit calculation of sewage treatment plant

项目	乙酸钠用量/( t·d−1)	曝气电耗/(kWh·d−1)	污泥产量/( t·d−1)	费用/(万元·d−1)
P1阶段	20	9 600	76	4.05
P2阶段	0	3 960	52	1.61

 | Show Table

DownLoad: CSV

3.   结论

1)长期低DO运行策略下，低C/N进水城镇污水厂(A²O工艺)逐渐停止外加碳源，在硝化液回流比由370%降低到90%的过程中，TN去除效果不降反升，平均去除率由76.3%提升至82.9%。这主要由于好氧区短程硝化和同步反硝化作用比例的升高所致。

2)短程硝化作用的内在原因为长期低DO条件对AOB抑制作用小，对NOB抑制作用明显，NOB的丰度从0.15%下降至0.06%，亚硝酸盐氧化还原酶功能基因NxrAB丰度显著下降，造成NO₂⁻-N积累。

3)好氧区同步反硝化作用的内在原因为低DO条件下具备反硝化功能的PAOs与GAOs比例升高，总丰度从2.06%提升至3.57%，好氧区低DO条件下利用厌氧期储存的内碳源完成反硝化脱氮；而内碳源利用率的提升又促进了PAOs的生长，提升了除磷效能。

4)低DO运行策略下，污水厂外加碳源、曝气电耗、污泥处理等费用大大降低，吨水处理费用节省0.3元。