养殖尾水达标排放是建设渔业强国和水产养殖绿色发展的内在要求。目前，我国海水工厂化养殖水体已达3.37×10⁷ m³ (2019年中国渔业统计年鉴)，未经处理的养殖尾水排放会影响近海生态环境，增加水体富营养化程度(温志良等, 2000; 舒廷飞等, 2002)。

由于海水盐度效应，单纯针对海水工厂化养殖尾水处理的技术较少，而常规的设施设备等处理技术则存在投资运行费用高、能耗高等不足(曲克明等, 2018; 乔卫龙等, 2019)。人工湿地通过基质、植物和微生物的物理、化学和生物协同作用来去除水体中的氮、磷和有机物，具有投资建设成本低、运行费用少、净化率高和易维护等优点，已引起国内外学者的广泛关注(夏汉平, 2002; 贺锋等, 2005; Ji et al, 2007)。

国内外有许多人工湿地应用在污水处理中的报道，但是大部分研究侧重于植物选取、基质改良、微生物筛选等方面(丁怡等, 2019; Zhang et al, 2016; Gary et al, 2017)，而对于水力负荷等因素考察较少，特别是在高盐度背景下，对人工湿地的净化效果及其影响因素(许永辉等, 2018)。研究发现，水力负荷、水深、水力停留时间是影响人工湿地运行的三大要素；适宜的水力负荷能减少人工湿地占地面积，提升处理效率(王世和等, 2003; 胡小芳等, 2008)。因此，本研究探讨不同水力负荷状态下复合垂直流人工湿地系统对牙鲆养殖尾水中主要污染物质的处理效果，分析系统内部化学需氧量(COD)、活性磷酸盐(PO₄^3–-P)、总氮(TN)、硝酸盐氮(NO₃^–-N)、亚硝酸盐氮(NO₂^–-N)和氨氮(NH₄⁺-N)的浓度变化特征，探究水力负荷对人工湿地处理海水养殖尾水的影响机制，为研究人工湿地处理海水养殖尾水提供数据支持。

1 材料与方法 1.1 人工湿地构建

构建了一套实验室规模的海水养殖尾水处理系统(图 1)，主要由牙鲆(Paralichthys olivaceus)养殖池、沉淀池、复合垂直流人工湿地和储水池组成。养殖池为圆柱体，直径×高为0.6 m×0.8 m，有效容积为130 L [π× (0.3 m) ²×0.46 m]；人工湿地长×宽×高为0.8 m × 0.6 m× 0.8 m，有效容积为300 L。人工湿地由面积均为0.1 m²的下行池和上行池2部分组成，有效表面积为0.1 m²，由下往上依次填入碎石(粒径5~15 cm)、煤渣(粒径3~10 cm)和细砂(粒径1~5 mm)；外置采样管，在各层填料中部取样，按照水流方向设置7个采样点(图 1)；滤砂层种植互花米草(Spartina alterniflora Loisel)，种植密度为60株/m²。系统各个单元通过水管和水泵连接。

1.2 实验内容

2019年8—10月开展实验。将人工湿地系统置于室外，处理对象为牙鲆养殖的尾水，每日喂食牙鲆2 h后，将尾水排入复合垂直人工湿地系统。尾水沉淀过滤之后，由水泵提升进入人工湿地，依次经过下行池和上行池，净化后进入储水池中。每次实验使用新鲜的养殖尾水，处理后的尾水不进行再循环。经过1个月的试运行，微生物群落形成，出水水质稳定，开展水力负荷对净化效果影响的实验。

参考《人工湿地污水处理工程技术规范》(HJ 2005—2010)，设置3组不同水力负荷，V₁=0.50、V₂=0.19、V₃=0.10 m/d，系统在V₃ (0.10 m/d)下运行时间为1.5 h，将其设定为所有实验组的运行时间。设定水力负荷之后，运行4 d，处理效果稳定后再进行取样分析，每组水力负荷共运行7 d，实验总时长为30 d。实验期间，天气良好、无雨天，避免了雨水进入系统造成的干扰。

在人工湿地各填料层设置取样点，通过水质指标监测系统运行状态。每组水力负荷运行的第7天，系统开始运行1.5 h后进行采样。测定的主要水质指标包括COD、PO₄^3–-P、TN、NO₃^–-N、NO₂^–-N、NH₄⁺-N、温度、pH和盐度。

1.3 分析方法

采用水质多参数测定仪(YSI 556，美国)对水温(T)、盐度(S)、pH进行现场测定，实验期间，进水温度和盐度分别为(23.78±0.17)℃和(30.11±0.09)，进水和出水的pH分别为7.76±0.08和8.57±0.05。TN使用总有机碳分析仪(TOC-VCPH)进行测定。其他水质指标按照《海洋监测规范》(GB17378.4-2007)进行测定，其中，COD采用碱性高锰酸钾法，PO₄^3–-P采用磷钼蓝分光光度法，NH₄⁺-N采用靛酚蓝分光光度法，NO₂^–-N采用N-1-萘基-乙二胺分光光度法，NO₃^–-N采用锌镉还原法测定。

采用Excel软件进行数据处理，对水力负荷与污染物去除率相关性分析，利用Origin 8.0软件作图。

2 结果与分析 2.1 COD的处理效果

COD浓度变化情况见图 2。如图 2所示，进水中V₁、V₂和V₃组COD的浓度分别为(5.75±0.30)、(3.70±0.45)和(2.79±0.59) mg/L，对应的去除率分别为28.23%、36.40%和36.25%。下行池所有实验组COD浓度均呈下降趋势，且在底部到达最低值。V₃组能保持稳定的COD浓度水平，而V₁组中的COD浓度起伏明显。其中，V₁组在最底部浓度为1.10 mg/L，去除率达到68.50%；而进入上行池之后，COD浓度出现一定程度的升高。

尾水中COD浓度较低，水力负荷对COD去除率的影响并不显著(P > 0.05)，初步沉淀过滤之后，基本上能满足《海水养殖水排放要求》(SC/T9103-2007) (COD≤10 mg/L)。COD进水浓度偏低，水力负荷对COD的影响并不显著。海水鱼类养殖尾水中COD的浓度较低，而水力负荷对于COD浓度较高的虾类养殖尾水等的去除率影响有待进一步研究。

2.2 PO₄^3–-P的处理效果

不同水力负荷下，PO₄^3–-P的浓度变化趋势一致(图 3)，水力负荷状态对PO₄^3–-P去除率的影响不显著(P > 0.05)。V₁、V₂和V₃组进水中PO₄^3–-P的浓度分别为(0.08±0.01)、(0.51±0.06)和(0.49±0.06) mg/L，对应的去除率分别为77.44%、79.08%和88.19%。降低水力负荷对系统的PO₄^3–-P吸附降解能力有一定提升。

经过下行池滤砂层之后，PO₄^3–-P浓度迅速下降，并在上行池依然保持稳定的水平，说明，在经过填料过滤吸附之后，系统对PO₄^3–-P的处理能力达到下限。在水力负荷适宜的情况下，PO₄^3–-P的吸附转化主要发生在下行池的中上层，即滤砂层。V₃组下行池的滤砂层对PO₄^3–-P的去除贡献率在62.03%，但V₁组在碎石层的PO₄^3–-P去除率才达到最大(82.41%)。水力负荷影响填料层对PO₄^3–-P的吸附降解能力，水力负荷越大，对PO₄^3–-P起主要吸附转化作用的填料层就越靠后，最终影响出水水质。

2.3 TN的处理效果

系统对尾水中TN的去除效果与COD/ρ (TN)见图 4。V₁、V₂和V₃组进水中的TN浓度分别为(3.68±0.38)、(3.47±0.39)和(3.16±0.42) mg/L，对应的去除率分别为49.50%、84.02%和85.90%。水力负荷越低，系统对TN的去除率越高，出水水质也越稳定(P < 0.05)。在下行池，TN的去除率均约为80%，去除效果明显；而在上行池，滤砂层出现了TN浓度升高的情况。水力负荷为0.19和0.50 m/d时，出水水质波动较大；而水力负荷为0.10 m/d时，水质保持稳定。

进水中COD/ρ(TN)在1~1.5之间，而经过处理之后，能达到4左右。下行池之中，V₂、V₃组COD/ρ(TN)在滤砂层之后才开始上升，说明，低水力负荷下微生物能利用大量被基质与植物根系拦截吸附的有机污染物，而系统下层对TN的消耗速率更快；而V₁组在煤渣层之后COD/ρ(TN)却开始下降，有机质的消耗速率要高于TN，这说明高水力负荷不利于系统中微生物获取溶解性污染物，使得脱氮效率变差。上行池中，V₃组COD/ρ(TN)上升点要比V₁、V₂组更靠前，适宜的水力负荷有利于脱氮菌群的活动，提升人工湿地系统整体的脱氮能力。

2.4 NO₃^–-N、NH₄⁺-N、NO₃^–-N的处理效果

V₁、V₂和V₃组进水中NO₃^–-N浓度分别为(2.47±0.28)、(2.49±0.34)和(2.60±0.36) mg/L。随着水力负荷的改变，NO₃^–-N、NH₄⁺-N和NO₃^–-N的浓度趋势类似(图 5)。V₁组条件下，出水中NO₃^–-N的浓度为0.88 mg/L，去除率为64.30%，NH₄⁺-N的去除率只有48.20%。V₂、V₃组对NO₃^–-N的去除率在92%以上，并且NH₄⁺-N和NO₂^–-N去除率在98%以上。

NH₄⁺-NN与NO₂^–-N的曲线变化类似于NO₃^–-N，水力负荷对于无机氮(DIN)各组分的影响具有一致性。上行池滤砂层，V₁组的NO₃^–-N的浓度增加了0.31 mg/L，NH₄⁺-N的浓度减少了0.09 mg/L，pH值降为8.38。硝化反应中，NH₄⁺+2O₂→NO₃^–+H₂O+2H⁺，消耗NH₄⁺-N会产生等质量NO₃^–-N，并产生一定量H⁺，但实际上却有0.22 mg/L多余的NO₃^–-N累积，水力负荷过大导致了NO₃^–-N的累积。

系统内部的pH变化情况如图 6所示，3组实验中，进水pH值为7.76±0.08，进入下行池后pH都迅速上升，系统底部则达到了8.57±0.05，出水pH都在8以上。在系统底层，V₃组NO₃^–-N的浓度减少了2.43 mg/L，pH值从7.73增加到8.74，反硝化反应消耗NO₃^–-N并积累了大量的OH^–离子，使水质呈碱性。进入上行池后，V₃组pH曲线保持稳定的下降趋势，其他2组曲线略有起伏。

2.5 不同水力负荷下DIN的占比

DIN由NH₄⁺-N、NO₂^–-N和NO₃^–-N组成。进水中，NH₄⁺-N和NO₃^–-N占DIN主要成分，经过处理后，NO₃^–-N则占主要成分。3组实验中，DIN在系统内部的占比如图 7所示。V₁组进水中，NO₃^–-N占DIN的75.30%，NH₄⁺-N和NO₂^–-N则占据25.70%；出水中，NH₄⁺-N和NO₂^–-N总占比达到29.30%，而NO₃^–-N则只有70.70%。而在V₂与V₃组进水中，NO₃^–-N分别占DIN的54.70%、86.00%，且出水中，V₂、V₃组的NO₃^–-N占比都在86%以上。

水力负荷对出水中DIN组成成分影响明显，水力负荷越大，出水中NH₄⁺-N和NO₂^–-N的占比就越大。低水力负荷下，各个形态下的氮污染物都能被微生物充分利用；高水力负荷下，氮污染物与微生物接触时间短，部分有机氮被转化为NH₄⁺-N后未被充分吸收转化，导致出水NO₃^–-N的占比降低。

V₁组出水中，NH₄⁺-N的占比要比进水中的高33%左右，而V₃组NH₄⁺-N的占比则降低了74.89%。V₁组中，NH₄⁺-N占比在20%左右，且全程变化不大；在V₂组中，NH₄⁺-N在下行池中占DIN的40%左右，而进入上行池NH₄⁺-N则开始减少；V₃组实验中，经过下行池煤渣层处理之后，NH₄⁺-N的占比开始提升，在上行池中，NH₄⁺-N的占比接近50%。

3 讨论 3.1 水力负荷对氮污染物去除的影响

V₁与V₂组的TN在上行池中浓度升高，这可能与颗粒态污染物的分解有关(赵联芳等, 2006)。高水力负荷下，污水在上行池中流速变缓，使颗粒态污染物得以被分解，导致了TN指标的升高。系统对氮污染物的去除主要在系统中下层，尤其是煤渣层以及碎石层。Rousseau等(2004)研究显示，湿地系统对NO₃^–-N主要通过厌氧条件下微生物的反硝化反应去除，而系统内部pH及DIN浓度的变化趋势都说明反硝化活动激烈。

各个形态氮污染物的浓度及其占比变化都表明，高水力负荷容易影响脱氮微生物的活动，不利于有机氮的转化和无机氮的吸收。这主要是由于高水力负荷下，部分脱氮细菌随水流带出系统，且污染物在系统内停留的时间过短，与微生物接触不充分。而系统中上层由于溶解氧的补充，使得硝化细菌更活跃，造成了系统末端NO₃^–-N或者NO₂^–-N的累积，影响出水水质，从V₁和V₂组中NO₃^–-N浓度的上升可以说明这一点。

V₃组中，NH₄⁺-N在系统空间内的占比变化可以说明，适宜水力负荷下，系统中氨氧化细菌不断的将有机氮转化为NH₄⁺-N，而反硝化细菌消耗NO₃^–-N速率更快，使得NH₄⁺-N的占比逐渐上升。从而可以认为，下行池煤渣层到上行池煤渣层之间，水力负荷适宜的情况下，氨氧化与反硝化细菌主导氮污染物的转化过程。吴晓磊(1995)认为，复合垂直流人工湿地系统在能够承受的水力负荷范围内出水水质波动不大。如表 1所示，不同类型的人工湿地系统中TN的去除率均与水力负荷负相关。适宜的水力负荷下，TN的去除率为50.30%~88.70%；水力负荷过大，则导致TN的去除率明显降低。康传磊等(2018)研究显示，异养硝化—好氧反硝化菌能够提高对氮污染物的去除效果；而对高水力负荷考察之后，高水力负荷人工湿地中水力停留时间是污染物去除效果的限制因素(高奇英等, 2018)。设定合理的水力负荷与水力停留时间，或者投加适应性更强的脱氮细菌，能够在高水力负荷的条件下，提升氮污染物处理效果。

3.2 水力负荷对pH和DO的影响

在系统中，微生物活动产生的副产物如OH^–、H⁺会改变水体的酸碱性，影响到菌群的生长，反硝化过程的最佳pH一般为7.0，pH小于或大于7.0，反硝化速率都会随之降低(徐乐中, 1996)。而水力负荷状态对pH的改变影响较小，但从V₃和V₁组中pH曲线变化可以看出，在系统后程，低水力负荷能保持pH平稳下降，高水力负荷状态则会导致pH上升。

研究表明，在反硝化过程中会产生大量的碱度(王淑莹等, 2008)，硝化过程则消耗碱度(王建龙, 2000)。实验中DIN与pH值的变化说明，上行池在低水力负荷下，硝化细菌更活跃，有利于NH₄⁺-N的去除；而高水力负荷不利于硝化反应的进行，使得出水中的NH₄⁺-N与NO₂^–-N浓度过高。张欢欢等(2012)发现，湿地进水有机负荷和温度对碱度有显著性影响，但目前相关的研究较少。

由于植物的生长，上中层DO含量更充足，所以系统的上中层硝化作用活跃。脱氮过程集中在系统中下层，水力负荷越低，进水中DO富集到填料中的越多，脱氮效果也就越明显。焦玉恩等(2017)研究证明，随着水力负荷的增加，填料层内各处DO浓度呈现先增加后降低的趋势。适宜的水力负荷能够维持人工湿地系统中的溶氧水平，保证脱氮率的提升。

3.3 水力负荷对COD和PO₄^3–-P的影响

COD作为衡量水体中耗氧有机物含量的重要指标，能够反映出水体的受污染程度(过锋等, 2012; 陈芳等, 2013)。水力负荷对COD的去除率影响较小，虽然进水中COD浓度不一样，但去除率却大致相近。和丽萍等(2014)也认为，人工湿地中，COD的去除率受水力负荷影响程度较小。

本研究中，低水力负荷条件下，系统内部的COD浓度稳定，而高水力负荷可能会造成系统中生物膜的脱落，使得COD浓度波动较大。系统对于低COD养殖尾水的处理效果一般，水力负荷更多的影响到上行池中COD浓度的变化，所以，对高COD浓度的海水养殖尾水在复合垂直流人工湿地中不同水力负荷状态下的处理效果考察十分必要。

PO₄^3–-P的吸附转化主要发生在下行池的中上层，即滤砂层附近，并且吸附转化比较彻底，使得上行池中PO₄^3–-P的浓度保持稳定；同时，植物发达的根系也起到了吸收利用的效果。基质的吸附是磷污染物主要的去除方式(Vymazal, 2007)，从V₃和V₁组在各填料层PO₄^3–-P的去除率可以发现，水力学因素会导致发挥吸附作用的基质层靠后，影响最终处理效果。因此，本研究中高水力负荷条件下PO₄^3–-P的去除率都不如低水力负荷条件，这与凌祯等(2011)的研究结果相似，磷污染物的去除率与水力负荷负相关。

4 结论

复合垂直流人工湿地系统对海水养殖尾水中的氮、磷污染物有良好稳定的去除效果。水力负荷对COD去除率影响不大，去除率最高只有36.25%；水力负荷为0.50 m/d时，PO₄^3–-P的去除率为77.44%；水力负荷降低到0.10 m/d时，PO₄^3–-P的去除率达88.19%。水力负荷越大，PO₄^3–-P的吸附转化越靠系统后程。水力负荷为0.50 m/d时，TN去除率只有49.50%；水力负荷降低到0.10 m/d时，TN去除率达85.90%。水力负荷状态同样影响脱氮能力，不同水力负荷状态下，下行池的氮污染物转化最多，去除率达到80%以上，高水力负荷会导致上行池中NO₂^–-N与NO₃^–-N累积，适宜的水力负荷才能维持系统的稳定。