环境科学  2014, Vol. 35 Issue (4): 1311-1317   PDF    
循环流廊道湿地中氮归趋过程模拟研究
彭剑峰1,2, 宋永会1,2, 袁鹏1,2, 张雪妍3, 胡小明3    
1. 中国环境科学研究院城市水环境科技创新基地,北京 100012;
2. 中国环境科学研究院环境基准与风险评估国家重点实验室,北京 100012;
3. 沈阳环境保护局东陵分局,沈阳 110179
摘要:针对新型高效脱氮循环流廊道(CFC)湿地,构建了涵盖6种氮形态、3类介质、10种代谢途径的N循环模型,探索了湿地内N迁移转化模式.结果表明,沸石吸附(53.3%)、植物吸收NH4+-N(27.6%)、NO3--N反硝化(10.2%)、植物吸收NO3--N(2.9%)和NO2--N短程反硝化(1.5%)对TN去除贡献依次降低.NH4+-N去除机制存在季节差异,其中1月NH4+-N主要通过沸石吸附去除(84.5%);4~6月通过植物吸收去除(76.4%~85.3%);7月通过沸石吸附(36.1%)、亚硝化(45.8%)及植物吸收(21.4%)共同去除.此外,定期收获植物、按期再生沸石及种植水生植物可分别提升TN去除率1.7%~7.7%、43.1%~72.2%和19.8%~36.2%.综之,CFC湿地去除途径多样性保障了TN的长期高效去除.
关键词循环流廊道湿地          迁移转化模型     沸石吸附     去除贡献    
Modeling Nitrogen Transformation in a Novel Circular-Flow Corridor Wetland
PENG Jian-feng1,2, SONG Yong-hui1,2 , YUAN Peng1,2, ZHANG Xue-yan3, HU Xiao-ming3    
1. Department of Urban Water Environmental Research, Chinese Research Academy of Environmental Sciences, Beijing 100012, China;

2. State Key Laboratory of Environmental Criteria and Risk Assessment, Chinese Research Academy of Environmental Sciences, Beijing 100012, China;

3. Dongling Branchof Shenyang Environmental Protection Bureau, Shenyang 110179, China

Abstract: Nitrogen transformation and removal in a novel circular-flow corridor (CFC) wetland was simulated. A computer model covering 6 nitrogen forms, 3 media and 10 transferring pathways, was developed. The results show that the TN removal is attributed to zeolite adsorption (53.3%), plants NH4+-N uptake (27.6%), denitrification (10.2%), plants NO3--N uptake (2.9%) and short-cut denitrification (1.5%). The major removal pathway for NH4+-N varies from season to season. In January, the zeolite adsorption played the most important role in NH4+-N removal with a contribution of 84.5%. From April to June, the plants uptake accounted for 76.4%-85.3% NH4+-N removal. While in July, there were three removal pathways, namely zeolite absorption (36.1%), nirtosation (45.8%) and plants uptake (21.4%), playing an important role in nitrogen removals. Additionally, regular plants harvesting, periodical zeolite regenerating and aquatic plants growing can improve the TN removal rate by 1.7%-7.7%, 43.1%-72.2% and 19.8%-36.2%, respectively. In short, it is by the diverse removal pathways in CFC wetland that nitrogen can be continuously removed.
Key words: circular-flow corridor wetland     nitrogen     transformation and conversion model     zeolite adsorption     removal contribution    

在人工湿地生态处理系统中,N的迁移转化过程复杂.它不仅涉及N在水体、 植物、 底泥(或填料)多介质间的迁移,还涉及到在NH+4-N、 NO-2-N、 NO-3-N和ON等不同形态间的转化,而湿地对N的最终去除效果取决于各代谢途径的协同作用[1,2,3,4].因而,如何通过优化运行环境,促进湿地内N代谢转向有利于TN去除的方向成为湿地技术研究的热点[5,6,7,8].

循环流廊道(CFC)湿地是一种新型人工湿地,它通过水体在不同类型湿地单元内的持续循环流动,实现对N和有机物的持续高效去除.该类湿地具有运行灵活、 抗冲击能力强,且造价及运行费用低的优点,可实现高浓度氨氮废水的高效低成本达标排放,因而具有较好的发展及推广前景[9].然而,目前循环流廊道湿地仅处于研发阶段,对N的迁移转化过程尚未系统研究,这一定程度上制约了该类湿地的发展及推广.

针对循环流廊道湿地中不同形态N在填料、 水体和植物多介质间复杂的循环转化过程,本研究构建了多形态、 多过程的CFC湿地氮归趋模型,表征了吸附、 生物同化吸收、 氨化、 亚硝化、 硝化和反硝化等作用机制对各形态氮代谢过程的贡献,以期为优化CFC湿地脱氮效果提供理论及参数支撑. 1 材料与方法 1.1 工艺构型

选择“回”字形布局的CFC湿地开展研究(如图1).该CFC湿地包括稀释区、 湿地区和储水区,待处理水在各功能区内循环流动.首先,高浓度污水由进水口直接进入稀释区,高浓度的COD和NH+4-N可以得到稀释; 经稀释和沉淀后由提升装置定时、 定量提升至廊道湿地区,污水以间歇流形式依次流经廊道湿地区1~5单元; 最终处理出水进入储水区,部分回用于农业用水、 养殖场洗涤用水或直接达标排放,剩余部分自流进稀释区,重新进入循环[9].

图1 循环流廊道湿地平面示意 Fig.1 Plane figure of the CFC wetland

该CFC湿地有效容积330 L,水力停留时间7~13 d,进水COD表面负荷为5.3~16.8 g ·(m2 ·d)-1,NH+4-N为9.5~18.3 g ·(m2 ·d)-1.为强化湿地脱氮能力,湿地廊道区前端2个单元(单元1、 2)采用大孔隙率陶粒(粒径8~12 cm)作为填料,后端3个单元(单元3、 4、 5)采用改性沸石功能填料(粒径5~8 cm),填料层厚40 cm,砾石承托层厚10 cm.湿地表面种植芦苇(Phragmites australis)、 蒲草(Acorus calamus Linn)等水生植物,种植密度分别为80~100株 ·m-2和5~6株 ·m-2[9]. 1.2 模型构建

考虑到CFC湿地内N的赋存形态及循环途径,N的迁移转化模型不仅考虑了NH+4-N、 NO-2-N、 NO-3-N、 ON、 沸石吸附N和植物所含氮等6种N形态,还涵盖了硝化/反硝化、 植物吸收/腐败、 挥发、 氨化和沸石吸附/解吸等10种作用过程,详见图2.

图1 循环流廊道湿地中氮迁移模式 Fig.2 Nitrogen transformation in CFC wetland
参照CFC湿地内N迁移模式,根据物质流代谢平衡关系,建立CFC湿地中氮迁移转化的动力学模型[10,11,12,13,14]:
式中,Q为湿地流量,m3 ·d-1V为湿地有效容积,m3c(NH+4-N)ic(NH+4-N)分别为进水及湿地内氨氮浓度,mg ·L-1c(NO-2-N)为湿地内NO-2-N浓度,mg ·L-1c(NO-3-N)为湿地内NO-3-N浓度,mg ·L-1; WPN为湿地内植物的N储量, mg ·g-1; ON为进水及湿地内有机氮浓度,mg ·L-1; ZEN为湿地内沸石吸附的N含量, mg ·g-1r为氮反应速率,各速率的计算方法见表1.
表1 循环流廊道湿地中各形态N反应速率 [12,13,14,15,16,17,18,19,20,21] Table 1 Reaction rate of various nitrogen forms in CFC wetland
1.3 分析方法

CFC湿地进水采用模拟高浓度有机废水-畜禽养殖废水,其中进水的NH+4-N为870 (600~1100) mg ·L-1,PO3-4-P为160 (75~230) mg ·L-1,COD为570 (320~1010) mg ·L-1,pH值5.5~8.0.该CFC湿地置于具有温度调节装置的温室内,年温度变化范围12~32℃.

CFC湿地中TN、 NH+4-N、 NO-2-N和NO-3-N采用标准方法分析[22]; WPN采用凯式定氮法测定; ZEN采用沸石吸附/解吸平衡实验计算.ON主要来自于进水有机物携带或水生植物腐败的释放物,其计算方法为[18, 9]

2 模型求解及有效性分析

微分方程组的求解采用4阶龙格-库塔法[9, 10].模型参数部分采用经验数值,部分通过最小二乘法最优化确定.模型校正及验证所需数据分别采用相邻两年的运行结果.图3为CFC湿地中NH+4-N、 NO-2-N、 NO-3-N、 ON、 WPN和ZEN的检测值和模型预测值对比.从中可知,检测值和模型预测值之间的相关系数依次为0.9439、 0.9087、 0.9014、 0.7822、 0.8780和0.8029,其中NH+4-N、 NO-2-N、 NO-3-N和WPN的相关系数均大于或接近0.90; ON和ZEN的模拟效果略差,但相关系数也接近0.80.由于整个湿地系统测定值和模拟值之间具有较好的相关性,且实测值在模拟值两侧均匀分布,因而该模型具有较好的模拟效果,能够有效模拟湿地中氮的迁移转化过程.

图3 CFC湿地模拟值与实测值相关性 Fig.3 Comparison between the simulated data and the measured data in CFC wetland
3 N迁移转化过程分析 3.1 NH+4-N的循环途径

在传统人工湿地中,亚硝化及植物吸收是NH+4-N的主要去除机制[23, 24],而在CFC湿地内NH+4-N的主导去除机制存在季节差异(如图4).由图4可知,随水温升高,1~7月沸石吸附、 亚硝化及植物吸收等去除途径对NH+4-N的去除速率明显增加,其中亚硝化增速最大,从0.09 mg ·(L ·d)-1上升至51.21 mg ·(L ·d)-1; 其次为沸石吸附,从0.53 mg ·(L ·d)-1上升至40.45 mg ·(L ·d)-1.尽管反应速率均有增加,但其对NH+4-N的去除贡献却存在明显季节差异.1月沸石吸附对NH+4-N的去除贡献高达84.5%; 在4~6月的生长期,植物吸收对NH+4-N去除贡献高达76.4%~85.3%; 而在7月沸石吸附、 亚硝化及植物吸收对NH+4-N的去除贡献依次为36.1%、 45.8%和21.4%.

图4 CFC湿地中NH+4-N的增长速率分布 Fig.4 Growth rate variation of NH+4-N in CFC wetland
这说明在CFC湿地内NH+4-N的主导去除机制存在明显季节差异,但多样化的去除机制确保了不同季节CFC湿地均具有较高的NH+4-N去除率.其中低温期沸石吸附是氨氮的主要去除途径; 4~6月植物吸收对NH+4-N去除贡献最大; 而高温期沸石吸附、 亚硝化及植物吸收共同实现NH+4-N去除,其去除贡献依次为亚硝化>沸石吸附>植物吸收 . 3.2 NOx的循环途径

CFC湿地中NO-2-N和NO-3-N的增长速率季节变化如图5所示.从中可知,1~5月NH+4-N亚硝化、 NO-2-N短程反硝化和硝化等途径对NO-2-N的平均增长速率为1.50、 -0.01和-0.66 mg ·(L ·d)-1,远低于NH+4-N的增长速率,因而低温期水体中NO-2-N不会大量积累.而6~7月,随着水温及DO增长,NH+4-N亚硝化和NO-2-N硝化速率分别增长至51.2 mg ·(L ·d)-1和8.16 mg ·(L ·d)-1,两者差值高达35.6~43.1 mg ·(L ·d)-1,导致NO-2-N在CFC湿地内快速积累.此外,模拟结果也表明高温期CFC湿地内存在显著的NO-2-N短程反硝化反应,其速率最高可达4.1 mg ·(L ·d)-1,对NO-2-N的去除贡献可达10.5%.综上,高温期亚硝化是CFC湿地NO-2-N变化的决定性因素.

图5 循环流廊道湿地中NO-2-N和NO-3-N的增长速率分布 Fig.5 Growth rate variation of NO-2-N and NO-3-N in CFC wetland
与NO-2-N的变化规律类似,1~5月CFC湿地内NO-3-N浓度低,导致植物吸收与反硝化对NO-3-N的平均去除速率仅有0.14 mg ·(L ·d)-1和0.49 mg ·(L ·d)-1,难以实现对水体中TN有效去除.而6~7月随着DO和NO-3-N的快速增长,植物吸收与反硝化速率增加到3.10 mg ·(L ·d)-1和9.71 mg ·(L ·d)-1,成为NO-3-N的主导去除机制.对比1~7月反硝化及植物吸收速率可知,反硝化对NO-3-N去除贡献是植物吸收的3.1倍,是CFC湿地中NO-3-N的主导去除机制. 3.3 ON和WPN的循环途径

CFC湿地中ON和WPN的增长速率季节变化如图6所示.从中可知,CFC湿地内ON变化主要由ON氨化及水生植物腐败共同决定,其中1~4月ON氨化速率和水生植物腐败速率的比值在0.01~0.19之间,湿地内ON的变化主要受植物腐败速率影响; 而5~7月两者比值快速增长到0.37~0.49,氨化对ON去除贡献增大,但植物腐败仍是ON的主导变化机制.这与模拟养殖废水中ON浓度持续偏低相一致.

图6CFC湿地中ON和WPN的增长速率分布 FIg.6 Growth rate variation of ON and WPN in CFC wetland

WPN增长速率也存在明显季节差异.1~4月,CFC湿地内NO-3-N含量低,植物吸收NH+4-N成为WPN的主要来源; 而6~7月,尽管植物吸收NH+4-N仍是WPN增长的主要来源,但是植物吸收NO-3-N和植物腐败释放对WPN的增长贡献增长到20.92%和-35.20%.这表明随着高温期水生植物快速生长以及水体内NO-3-N的大量积累,植物的吸收NH+4-N、 腐败释放及吸收NO-3-N共同决定了WPN的变化特征.

4 TN去除季节差异

受反应速率季节差异影响,CFC湿地内TN的去除也存在明显季节差异. 1月TN物质流速率仅有2.38 mg ·(L ·d)-1,而7月增长到88.24 mg ·(L ·d)-1.这意味着高温期有更多TN参与到物质循环中,相应就有更多TN通过植物吸收、 挥发、 反硝化等途径从CFC湿地内去除.图7为CFC湿地内N物质流季节分布差异.如果仅考察CFC湿地内TN变化,可知水体中TN的去除途径多样,有NH+4-N挥发、 NO-2-N/NO-3-N反硝化、 沸石吸附、 植物吸收NH+4-N/NO-3-N和植物释放,各途径对TN的去除贡献依次为沸石吸附(53.3%)>植物吸收NH+4-N(27.6%)>NO-3-N反硝化(10.2%)>植物吸收NO-3-N(2.9%)>NO-2-N反硝化(1.5%)>NH+4-N挥发(<0.1%).此外,TN主导去除机制也存在季节差异. 1月沸石吸附对TN的去除贡献为72.2%,是TN的主导去除机制; 4月沸石吸附及植物吸收对TN的去除贡献最大,分别为44.7%和40.3%; 7月沸石吸附、 植物吸收及反硝化对TN去除贡献最大,依次为43.1%、 29.6%和19.5%.

图7 湿地内N物质流季节差异 Fig.7 Comparison of N material flow rate in different seasons

内环为1月; 中环为4月; 外环为7月

分析N物质流差异有助于优化各技术措施,提升CFC湿地对TN去除效果[25].由模拟结果可知,定期收获植物可提升CFC湿地TN去除率1.7%~7.7%; 未定期再生功能填料,湿地对TN的去除率1月下降72.2%, 7月下降43.1%; 不种植水生植物湿地对TN的去除率下降19.8%~36.2%.因而,定期维护CFC湿地运行,保障水生植物生长对于提升TN去除效果具有重要作用.

5 结论

(1)CFC湿地内NH+4-N的主导去除机制存在明显季节差异,但多样化的去除机制确保了不同季节CFC湿地均具有较高的NH+4-N去除率.低温期沸石吸附是氨氮的主要去除途径; 4~6月植物吸收对NH+4-N去除贡献最大; 而高温期沸石吸附、 亚硝化及植物吸收共同实现NH+4-N去除.

(2)水体内TN主导去除机制也存在季节差异. 1月沸石吸附是TN的主导去除机制; 4月沸石吸附及植物吸收对TN的去除贡献最大; 7月沸石吸附、 植物吸收及反硝化对TN去除贡献最大.各途径对TN的去除贡献依次为沸石吸附(53.3%)>植物吸收NH+4-N(27.6%)>NO-3-N反硝化(10.2%)>植物吸收NO-3-N(2.9%)>NO-2-N反硝化(1.5%)>NH+4-N挥发(<0.1%).

(3)定期维护CFC湿地运行有助于提升TN去除效果.定期收获植物可提升CFC湿地TN去除率1.7%~7.7%; 未种植植物降低TN去除率19.8%~36.2%; 未定期再生功能填料降低TN的去除率43.1%~72.2%.

参考文献
[1] Itoh M, Takemon Y, Makabe A, et al. Evaluation of wastewater nitrogen transformation in a natural wetland (Ulaanbaatar, Mongolia) using dual-isotope analysis of nitrate[J]. Science of the Total Environment, 2011, 409 (8): 1530-1538.
[2] Saeed T, Sun G Z. The removal of nitrogen and organics in vertical flow wetland reactors: Predictive models[J]. Bioresource Technology, 2011, 102 (2): 1205-1213.
[3] Vymazal J. The use of hybrid constructed wetlands for wastewater treatment with special attention to nitrogen removal: a review of a recent development[J]. Water Research, 2013, 47 (14): 4795-4811.
[4] Hu Y S, Zhao Y Q, Zhao X H, et al. High rate nitrogen removal in an alum sludge-based intermittent aeration constructed wetland[J]. Environmental Science Technology, 2012, 46 (8): 4583-4590.
[5] 范改娜, 祝贵兵, 王雨, 等. 河流湿地氮循环修复过程中的新型功能微生物[J]. 环境科学学报, 2010, 30 (8): 1558-1563.
[6] Hu Z, Lee J W, Chandran K, et al. Nitrogen transformations in intensive aquaculture system and its implication to climate change through nitrous oxide emission[J]. Bioresource Technology, 2013, 130 : 314-320.
[7] 林建伟, 朱志良, 赵建夫. 天然沸石覆盖层控制底泥氮磷释放的影响因素[J]. 环境科学, 2006, 27 (5): 880-884.
[8] Borin M, Salvato M. Effects of five macrophytes on nitrogen remediation and mass balance in wetland mesocosms[J]. Ecological Engineering, 2012, 46 : 34-42.
[9] Peng J F, Song Y H, Liu Z H, et al. Performance of a novel Circular-Flow Corridor wetland toward the treatment of simulated high-strength swine wastewater[J]. Ecological Engineering, 2012, 49 : 1-9.
[10] Peng J F, Wang B Z, Song Y H, et al. Modeling N transformation and removal in a duckweed ponds: model development and calibration[J]. Ecological Modelling, 2007, 206 (1-2): 147-152.
[11] Peng J F, Wang B Z, Song Y H, et al. Modeling N transformation and removal in a duckweed ponds: model application[J]. Ecological Modelling, 2007, 206 (3-4): 294-300.
[12] Edelfeldt S, Fritzson P. Evaluation and comparison of models and modelling tools simulating nitrogen processes in treatment wetlands[J]. Simulation Modelling Practice and Theory, 2008, 16 (1): 26-49
[13] 彭剑峰, 宋永会, 袁鹏, 等. 浮萍塘中氮归趋模式模拟分析[J]. 环境科学, 2006, 27 (10): 1963-1968.
[14] Lv T, Wu S B, Hong H, et al. Dynamics of nitrobenzene degradation and interactions with nitrogen transformations in laboratory-scale constructed wetlands[J]. Bioresource Technology, 2013, 133 : 529-536.
[15] Pai T Y, Shyu G S, Chen L, et al. Modelling transportation and transformation of nitrogen compounds at different influent concentrations in sewer pipe[J]. Applied Mathematical Modelling, 2013, 37 (3): 1553-1563.
[16] Wang Y H, Zhang J X, Kong H N, et al. A simulation model of nitrogen transformation in reed constructed wetlands[J]. Desalination, 2009, 235 (1-3): 93-101.
[17] 韩建刚, 曹雪. 典型滨海湿地干湿交替过程氮素动态的模拟研究[J]. 环境科学, 2013, 34 (6): 2383-2389.
[18] Tang W Z, Zhang W Q, Zhao Y, et al. Nitrogen removal from polluted river water in a novel ditch-wetland-pond system[J]. Ecological Engineering, 2013, 60 : 135-139.
[19] 刘晓娜, 丁爱中, 程莉蓉, 等. 潜流人工湿地除氮的生态动力学模拟[J]. 农业环境科学学报, 2011, 30 (1): 166-170.
[20] 白军红, 欧阳华, 邓伟, 等. 湿地氮素传输过程研究进展[J]. 生态学报, 2005, 25 (2): 326-333.
[21] McGechan M B, Moir S E, Sym G, et al. Estimating inorganic and organic nitrogen transformation rates in a model of a constructed wetland purification system for dilute farm effluents[J]. Biosystems Engineering, 2005, 91 (1): 61-75.
[22] 国家环境保护总局. 水和废水监测分析方法[M]. (第四版). 北京: 中国环境科学出版社.
[23] Wolf K, Ahn C, Noe G B. Microtopography enhances nitrogen cycling and removal in created mitigation wetlands[J]. Ecological Engineering, 2011, 37 (9): 1398-1406.
[24] Gómez-Rey M X, Couto-Vázquez A, González-Prieto S J. Nitrogen transformation rates and nutrient availability under conventional plough and conservation tillage[J]. Soil and Tillage Research, 2012, 124 : 144-152.
[25] Yue F J, Li S L, Liu C Q, et al. Using dual isotopes to evaluate sources and transformation of nitrogen in the Liao River, northeast China[J]. Applied Geochemistry, 2013, 36 : 1-9.