我国村镇生活污水处理率较低，处理设施比较落后，污水厂尾水氮磷等污染物浓度较高。目前，城镇污水厂执行GB 18918-2002《城镇污水处理厂污染物排放标准》一级A标准，其ρ(TN)(15.0 mg·L-1)和ρ(TP)(1.0 mg·L-1)比GB 3838-2002《地表水环境质量标准》V类水标准分别高6.5倍和4倍，污水处理厂尾水已经成为水体富营养化的重要污染源。对污水处理厂尾水进行深度净化，源头削减N和P含量，对湖泊或入湖河道水质改善具有重要意义。目前，传统尾水深度处理技术包括混凝一沉淀一过滤、膜过滤、曝气生物滤池和高级氧化等。但传统处理技术具有费用高、容易造成二次污染等缺点。水生植物修复技术具有效果好、投资少、运行成本低、易管理和景观效果好等优点。目前，水生植物修复技术已被广泛应用于湖泊、生活污水尾水等污染水体的生态处理中。
水生植物修复污染水体过程中，因有机物的降解及氮磷的去除，释放温室气体，对环境造成二次污染。温室效应造成的气候变化引起了人们的广泛关注，温室气体浓度的增加是引起温室效应的主要原因，因此温室气体的“源一汇”受到了广泛的关注。大气中CO2,CH4、和N2O的浓度增加对温室效应增强的总贡献率占了将近80 %，是温室效应的主要贡献者，且其大气浓度仍分别以年均0.5% ,0.8%和0.3%的速率在增长。目前，对于温室气体排放的研究多集中于农田、水库、湖泊及天然湿地等方面，对于污水处理过程中温室气体(CO2,CH4、和N2O)排放研究很少，而水生植物修复污染水体过程中温室气体排放的研究鲜见。
依托生态治理工程，采用江苏省农业科学院自主研发的原位收集和释放气体装置，监测凤眼莲( Eachhornaa crassapes)深度净化污水厂尾水过程中温室气体(CO2,CH4、和N2O)排放通量的季节变化特征和沿程变化特征，并探讨温室气体排放通量的相关环境因素，为凤眼莲深度净化污水厂尾水生态工程提供理论支撑。
1 材料与方法
1.1 污水处理厂与深度净化塘概况
南京市高淳区东坝污水处理厂(31。17'28.0" N ,119。02'29.3" E)，主要污水来源于东坝镇及附近的生活污水，采用A20工艺处理污水，日接纳污水能力为2 000 t，实验期间日均处理生活污水1 024 t o未构建尾水深度净化生态工程前，生活污水经污水厂处理后直接排入连通太湖的青河。
如图1所示，深度净化塘采用三级串联方式组成。深度净化塘各级长度均为105 m，深1.2 m，其中第一级深度净化塘宽为25 m，第二、三级深度净化塘为27.5 m，总有效容积为7 500 m3，之间采用土夯方式隔开，深度净化塘底部和岸堤均铺设防水布防止底部渗漏至地下水。进水口和出水口均设置流量计监测污水净化量。出水口设置溢流堰保持深度净化塘水深为1 m。污水厂尾水全部进入深度净化塘，其水力负荷为((0.13 ± 0.03) m3·m-2·d-1，TN负荷为((1.21 ± 0.10) g·m-2·d-1, CODMn负荷为（0.57 ± 0.02) g·m-2·d-1, TP负荷为(0.05 ± 0.00)mg ·m-2·d-1。 2015年5月底凤眼莲种苗投放完毕，种苗投放量为0.6 kg·m-2。在进水口、一级、二级及三级净化塘出水口沿程设置4个监测点(图1)，将采气装置放置在监测点连续采气，并在附近设置水质监测点采集水样。〕
1.2 进水情况
该尾水深度净化生态工程进水为高淳县东坝污水处理厂尾水，尾水水质执行GB 18918-2002一级A标准，水质因季节和时节不同有所差别。工程运行期间，污水处理厂尾水ρ（TN）为(9.27±3.31)mg·L-1 ,ρ( TP)为(0.39±0.O5)mg·L-1 }P } NHa+-N)为(0.4910.07) mg·L-'，CODn, 为(4.3810.65)mg·L-'，水体p( DO)为(5.4012.21)mg·L-' , pH值为7.3610.28 。
1.3采样及分析方法
采用江苏省农业科学院自主研发的气体收集装置(图2)采集气体，综合考虑凤眼莲的生长特征、温度和产气量变化等因素，在8-9月，一次采气过程持续7d，连续采气，采集3次;10-11月，一次采气过程约持续15 d，连续采气，采集2次。为减少误差，统一在上午8 ; 00-11 ; 00采集气体，气体的采集和测定方法参考文献[21}。每个采样点设置3套采气装置。当集气罩内气体积聚形成气泡时，根据排水集气法原理自动将气体吸入集气瓶，通过集气瓶的质量变化来计算产气量。采用气相色谱仪测定各气体组分浓度，采用峰面积外标法定量各气体浓度，各组分气体释放通量的计算方法为
En=Cn} X Pn}, X E，(1)
E=(oielTj X 2}3.} s/(2}3.} s+t)，(2)
V=(W，一巩)/D。(3)式(1)一(3)中，乓a为气体释放通量，即单位面积水体单位时间释放气体的量，g·m_zm·h-'-乓a为气体组分浓度，%;pn}为标准状态下被测气体密度，g·L-' ; E为标准温度标准压力下水体释放气体的速率，mL·m_zm·h-' ; V为收集的气体体积，L;S为集气罩覆盖水体的面积，mz ; T为收集气体所用时间，h;，为收集气体过程中的平均温度，℃;W，为试验开始前装满水的集气瓶质量，g } }z为收集气体结束后集气瓶质量，g;D为室温(O}t}50℃)下水的密度，g·mL-'。
采用德国SEAL AA3连续流动分析仪测定进水及各级出水总氮(TN),铰态氮(NHQ'-N) ,硝态氮}N03--N}和总磷ATP)浓度，采用酸性高锰酸盐滴定法测定高锰酸盐指数(CODM )，采用多功能水质测定仪(YSI Pro Plus, USA)现场测定水温(c) , DO浓度和pH值。每隔15 d采集凤眼莲植株，采用重量法现场测定生物量。〕
1.4数据分析
采用Excel 2007和Sigmaplot 12.5软件进行数据整理和相关性分析，用Origin 8.5软件作图。统计检验显著性水平为a = 0.OS〕2结果与分析2.1试验期间水体主要理化指标变化
2015年6-11月，深度净化塘凤眼莲单位面积生物量和总生物量分别由(0.6010.09) kg " m-Z和(4.5010.64) t增至(22.73 1 2.82 ) kg " m-Z和(170.50121.17 ) t。由表1可知，水体温度变化范围为13 } 27 0C , 8月水温最高。DO浓度变化维持在3.0 7.0 mg " L-‘之间，属好氧状态，10-11月进水DO浓度大幅增高，各级出水DO浓度也呈递增趋势。水体pH值基本维持在7.07.6左右，属于微生物硝化反硝化的最佳pH值范围，随月份推移变化的幅度高于沿水流方向上的变化幅度。由上述结果可知，凤眼莲三级净化生态工程水体主要理化指标季节变化较明显，沿程变化较小，基本维持在一个较稳定的生态系统中。〕
水体氮磷污染物指标如图3所示，水体主要污染物TN , NHQ'-N , N03--N及TP都得到有效降解。监测周期内，进水p(TN) ,p(NHQ'-N) ,p(N03 -N)及p(TP)平均值为9.27,0.49,7.63和0.39 mg·L-'，三级净化出水平均值为2.96,0.21,2.20和0.14 mg·L-'，其中TN浓度接近地表V类水标准，TP浓度优于地表V类水标准，三级净化去除率达68.07% ,71.14% , 57.28%和64.21 %，凤眼莲深度净化生态工程对污水厂尾水具有明显的氮磷去除及水质改善效果。监测周期内，进水CODn, 均值为4.38 mg } L-'，三级净化出水均值为4.75 mg·L-'略高于进水，原因可能是污水厂尾水CODn, 处于较低水平，深度净化塘对尾水有机物的进一步去除效率不高，且水生植物根系的分泌物会在一定程度上增加CODn, 。三级净化出水CODn, 低于111类水标准〕
2.2温室气体排放特征
2.2.1排放通量
2015年8-11月，对凤眼莲深度净化生态工程中温室气体(COZ,CH、和NZO)排放进行监测，根据每月实际采样分析结果，计算凤眼莲深度净化尾水系统中COz,CH、和Nz0的月平均排放通量(表2) o
表2显示，凤眼莲深度净化塘COZ,CH、和NZO排放通量范围分别为。}0.136,0}0.263和0.6082.561 mg·m_Zm·h-'，平均排放通量为0.05 8 , 0.076和1.539 mg } m-Z } h-'。在整个试验周期内，凤眼莲深度净化塘累积排放1.273 kg C0z,1.685 kg CHQ及33.590 kg NzO。
2.2.2月份变化特征
如表2所示，随着月份变化，COZ , CH、排放通量呈现明显降低趋势，8月排放通量达最大值，排放通量分别为0.136和0.608 mg } m-Z } h-' ,10月和11月排放通量接近零，这可能与冬天水温降低及DO浓度、pH值升高有关。由表3可知，COZ和CHQ排放通量与水温的相关系数分别为O.s67(P<0.Os)和0.s24(P<0.Os)，呈显著正相关关系;COZ排放通量与 DO浓度、pH值的相关系数分别为-o.sss ( P<o.os ) , -o.606 ( P< o.os，呈显著负相关关系;CH、排放通量与DO浓度、pH值的相关系数分别为一0.3s4和一0.471，呈负相关关系，但相关性不显著。
NZO排放通量没有明显的季节变化趋势，排放通量从大到小依次为9,11,10和8月。9月排放通量达最大值，为2.s61 mg·m_zm·h-' } NZO是硝化过程中的副产物，反硝化过程的中间产物，是不完全硝化或不完全反硝化的产物。研究表明，NZO的生成及排放与水温、DO浓度、pH值、底物浓度及植物覆盖度等因素密切相关。该研究中Nz0排放通量与水温、DO浓度及pH值相关系数分别为-0.130,-0.217和一0.178，均未表现出相关性。
2.2.3沿程变化特征
三级净化生态工程温室气体排放通量沿程变化特征如图4所示。在沿程方向上，温室气体排放通量呈现出先升高后降低趋势，呈现明显的沿程变化特征，总体上进水端高于出水端。COZ排放通量在二级净化塘出水口达到最大值，排放通量为0.092 mg·m_Zm·h-' , C H、和Nz0在一级净化塘出水口达到最大值，排放通量分别为0.178和3.657mg " m_Z " h_'。由表1可知，沿程方向上水温没有明显变化，DO浓度维持在好氧状态，且pH值维持在在最佳范围，NZO产生量与碳氮浓度密切相关，排放量与水生植物覆盖度有关，TN和N03--N呈递减趋势。相关性分析结果(表3)表明，NZO排放通量与TN和N03--N相关系数分别为0.477和0.428 ,呈正相关关系。
3讨论
3.1凤眼莲三级净化生态工程温室气体排放通量
与相关研究相比，该研究中COZ和CH、排放通量较小，NZO排放通量较大。沙晨燕等[z3}运用静态箱一气相色谱法对Olentangy河湿地4种不同类型河滨湿地的CH、和COZ排放通量进行研究，发现不同类型河滨湿地CH、和COZ排放通量从大到小依次为自然湿地(( 0.33 } 85.7 mg } m-Z } h-' )、人工湿地( 0.02 20.5 mg·m_zm·h-')和半人工湿地(-0.040.09 mg } m-Z } h-' ) , COZ排放通量由大到小依次为自然湿地(13.1 } 53.5 mg } m-Z } h-' )、半人工湿地(一0.7一132.9 mg·m_zm·h一‘)和人工湿地(一13.3-51.6 mg·m_zm·h-' )。黄国宏等应用封闭箱法对辽河三角洲芦苇湿地CH、释放通量的研究结果表明，在5-11月，其释放通量为一968 } 2 734 },g·m_2m·h-' } WU等[251利用人工湿地系统处理污水的研究表明，潜流和表面流人工湿地系统N20平均通量为296.5和28.2 },g·m_Zm·h-'，远低于笔者研究结果。根据KHALIL等对全球N20产生源的估计，污水处理过程N20年释放量为0.3x10'2一3.Ox10'2 kg，占全球N20总释放量的2.5%一25 % } KA-MPSCHREUR等综合分析相关文献得到:在实验室规模的生物脱氮过程中可能有。一90%的氮会转化为N20;在大规模城镇污水厂的污水生物脱氮过程中可能有。一14.6%的氮转化为N20}
3.2 COZ和CHq排放通量影响因素
尾水深度净化生态工程系统内，C02和CH、主要通过植物传输由水体进入大气，植物传输受水生植物种类、覆盖度及植物传输机制的影响。水温不仅可以通过影响气体分子的扩散速度及其在水体中的溶解度来直接影响气体交换通量，还可以通过影响微生物活性间接影响温室气体产生的地球化学过程[2A1。监测周期内，C02和CH、释放通量与水温呈显著正相关关系，这与以往的研究结果[zy-3z}相一致。pH值直接影响水体碳酸盐体系(C02 , C032和HC03-)的动态平衡及分布，控制水体C02浓度,，水一气界面C02交换通量与pH值通常表现为负相关关系。笔者研究结果表明:COZ释放通量与pH值呈显著负相关关系，CH、释放通量与pH值呈负相关关系，与以往研究结果相同。但COZ和CH、排放通量与凤眼莲生物量呈显著负相关关系，与以往研究结果不一致。这可能是因为水温是控制COZ和CH、排放的关键因素，11月凤眼莲生物量增加，但生长缓慢，水温下降幅度很大。
TREMBLAY等[351的研究显示:DO浓度与水库中COZ,CH、释放通量呈显著负相关关系。沉积物中产生的甲烷不完全进入气泡中，一部分通过扩散上升到水面。上升过程中，由于DO浓度逐渐升高，产生的大部分甲烷被有氧一缺氧临界面的甲烷氧化菌消耗。笔者研究发现，COZ释放通量与DO浓度呈显著负相关关系，CH、释放通量与DO浓度呈负相关关系。对碳循环而言，有机物在有氧状态下产生COZ和CHQ，在缺氧状态下主要产生CHQ，因此，COZ和CH、排放通量与水体有机物浓度有关。笔者研究中COZ,CH、与CODn, 无相关性，可能是因为进水有机物浓度过低，基本不降解，因此由有机物降解产生的COZ和CH、量很少。
3.3 NZO排放通量影响因素
水温直接影响微生物活性及酶活性，笔者研究结果表明，NZO释放通量与水温没有相关性，这与以往研究结果不符，但目前对于水生植物修复技术及人工湿地处理系统中水温与NZO释放的相关关系没有明确结论。可能是由于水生植物的存在造成了复杂的硝化一反硝化微生物环境，不是简单的水温影响微生物活性进而影响NZO产生的过程。有研究表明在植物生长季，由于植物组织向根系传输了更多氧气，改变了根际溶氧微环境，从而促进人工湿地系统释放出较多。但也有研究表明人工湿地系统的最高释放量发生在植物枯萎衰败的秋季。笔者研究结果显示:11月，凤眼莲开始腐败脱落，NZO释放通量开始增加，此与上述研究结果相符。植物可通过吸收作用除氮，植物生物量越多，吸收的氮也越多，NZO的排放就越少该研究结果显示NZO排放通量与凤眼莲生物量呈正相关(P>0.OS)，与其他文献结果不一致。
pH值通过影响微生物的活性间接影响NZO释放通量，微生物活性一般在中性或弱碱性环境下最高，pH值越低，NZO释放通量越大，两者之间呈负相关关系[ao。笔者研究中，NZO释放通量与pH值没有相关关系，可能是pH值变化范围较小，基本维持在最佳的反应条件，pH值不是控制CH、和NZO产生的关键因素，而是其他因素造成Nz0释放通量的变化。NZO是硝化过程中的副产物，反硝化过程的中间产物。硝化过程中DO浓度过低是造成NZO产生的最主要原因;反硝化过程中DO浓度过高可导致NZO还原酶活性降低或失活进而造成NZO积累。
4 结论
(1)通过凤眼莲生态工程深度净化污水厂尾水，出水水质得到较大改善。出水p STN)和P }Tp}分别为(2.9611.77)和(0.1410.08) mg·L-'，远低于GB 3838-2002一级A标准。
(2)工程运行期间，凤眼莲三级净化污水厂尾水生态工程COZ,CH、和NZO平均排放通量分别为0.058 ,0.076和1.539 mg·m_Zm·h-'，COZ和CHQ排放具有明显的月份和沿程变化特征，NZO排放具有明显的沿程变化特征，但没有明显的月份变化特征。
(3)通过对监测结果的相关性分析可知，COZ和CH、排放通量与水温具有显著正相关关系;COZ排放通量与pH值和DO具有显著负相关关系;CHQ排放通量与pH值和DO呈现负相关关系，但相关性不显著;NZO排放通量与水体TN和N03--N浓度呈正相关关系。〕
