沉水植物型生态净化系统处理农田退水的总磷去除动力学研究

　　摘要：【目的】探讨以沉水植物为先锋物种的生态净化系统处理农田退水磷污染的动力学性能。【方法】通过批式试验研究苦草Vallisnerianatans型和金鱼藻Ceratophyllumdemersum型2种沉水植物型净化系统对农田退水总磷的净化效率及其动力学特征。【结果】至试验结束时(第49天)，2种植物净化系统总磷去除效率无显著差异(P>0.05)，总磷去除率分别为82.8%(苦草型)和84.0%(金鱼藻型)，但两者总磷去除效率的差异在时间尺度上存在不确定性;动力学模拟分析发现：除Grau二级动力学外，一级动力学、Monod 动力学和修正的Gompertz模型均可以描述试验条件下总磷的去除过程[判定系数(R2)>0.930，相对均方根误差(RRMSE)<0.200]，其中Monod动力学和Gompertz模型具有更高的拟合度(R2>0.970)，预测值与实验观测值之间吻合程度更好(RRMSE<0.110)。【结论】3种有效模型拟合获得的动力学常数在植物种类水平上均存在显著差异(P<0.05)，其中指示总磷去除效率的动力学常数值一致表现为苦草型>金鱼藻型，表明苦草型净化系统除磷效率好于金鱼藻型。

　　关键词：沉水植物;生态净化系统;农田退水污染;总磷去除效率;动力学模拟

　　农业面源污染已成为全球许多国家流域水污染和水体富营养化的首要污染来源[1−2]。近年来，国内诸多流域尺度的污染源调查解析表明：中国许多河湖、水库等地表水体中的氮磷输入主要来自农田种植、畜禽水产养殖等农业生产以及农村生活所产生的非点源污染，其中种植业所引起的农田退水污染正逐步成为污染的主要贡献者[3−5]。农田退水污染具有排放时间瞬时多变，迁移途径复杂多样，污染负荷时空变化幅度大等特点，治理难度相对更大[1−2]。目前基于生态工程原理的生态沟渠、人工湿地等技术受到格外重视[6−8]。

　　水资源生态论文： 水生植物凤眼莲入侵对湖泊湿地甲烷排放的影响

　　此类通常建立在培植大型水生植物基础上的生态拦截工程具有建设成本低且可控、技术难度低、生态友好等优点，近年来在中国许多地区特别是生态敏感区域得到了较大规模的推广应用。2018年以来，浙江省在全省范围内实施了300项农田面源污染氮磷生态拦截沟渠示范推广项目[9]。研究已表明：水生植物对包括氮磷营养物质在内的众多污染物的净化效率及其稳定性具有重要的影响[10−11]。沉水植物占据了根系-基质(根际)、水-基质和植物茎叶-水3个关键微界面，由其介导形成的微环境结构更为复杂，往往能够表现出更加稳定和持续的水质净化潜力[12−14]。

　　近几年以沉水植物为先锋物种的“水下植被”技术拓展应用到农业面源污染的拦截处理[15−17]。目前，国内外对挺水植物型人工湿地去污动力学已开展了比较深入的研究[18−20]，但对“水下植被”型净化系统污染物去除动力学特征还十分欠缺，特别是对拦截净化农田退水、池塘养殖尾水等具有瞬间排放特点的非点源污染的除磷动力学研究至今还未见报道。

　　在动力学模型应用研究方面，一级动力学、Grau二级动力学和Monod动力学因具有模型构式可线性化处理、参数求解过程简单、适用范围广等优点而被广泛采用[19−22];近年来，修正的Gompertz模型作为一种多参数非线性模型在污染物去除动力学研究领域逐步被重视，并且已被成功应用于新型生物处理系统的工艺设计和污染物去除效率的预测[23−25]。

　　本研究以苦草Vallisnerianatans和金鱼藻Ceratophyllumdemersum等2种典型沉水植物为供试材料，研究苦草型和金鱼藻型2种“水下植被”型净化系统对农田退水中磷污染的净化效率。在此基础上，采用一级动力学、Grau二级动力学、Monod动力学和修正的Gompertz模型模拟总磷的去除过程并进行统计验证，获取相关动力学参数，进而探讨2种沉水植物净化系统除磷效率的差异，以期为“水下植被”型生态净化系统的工艺调控及其设计应用提供理论依据和技术参考。

　　1 材料与方法

　　1.1 供试材料

　　选择不同生长型的沉水植物苦草和金鱼藻作为供试植物，这2种植物均为长江中下游流域常见的优势沉水植物，偏好于富营养化水体，在溪沟、河流、池塘、湖泊等水域中均有分布。2种供试植物均取自浙江农林大学乌镇大学生实践基地。采集时选取生长状态良好且株高相近的植株，运至实验基地培养10 d，统一修剪至株高25 cm，去离子水清洗干净后用作试验材料。受降雨强度、降雨历时、施肥量、土壤类型、土地利用类型等因素影响，农田退水中总磷质量浓度波动极大，低值时可低于地表水Ⅳ类标准，而峰值时高达5.00 mg·L−1以上[26−27]。

　　本研究人工配置中，农田退水总磷质量浓度处于相对较高水平，试验原水取自浙江农林大学平山实验基地内的池塘，原水中总磷质量浓度为 0.17 mg·L−1，NH3-N质量浓度为0.39 mg·L−1，重铬酸盐指数(CODcr)为8.56 mg·L−1。经配置后的试验用水总磷质量浓度为4.17 mg·L−1。

　　1.2 试验方案

　　考虑到农田退水排放的瞬间性、间歇性等特征以及池塘、相对封闭河道等受纳水体的水动力特点，本研究采用批式试验考察2种供试植物净化系统对农田退水总磷的去除效果及其动力学特征。试验共设置3组，分别为苦草组、金鱼藻组和无植物对照组，每组均重复3次。试验于浙江农林大学平山实验基地温室大棚内进行，试验装置采用白色聚乙烯塑料箱(长×宽×高为40 cm×80 cm×60 cm)，各塑料箱内注入试验用水，试验水深50 cm，计为160 L。

　　每个塑料箱底部铺设1层厚度为10 cm火山石基质(过10目筛)。为避免植物直接定植在基质上导致基质分布不均，先将经预培养的供试植物(鲜质量300 g)定植于渗排水网垫网上(厚度为5 mm)，然后再铺设于基质层上。试验周期共计49 d(2018年7−8月)，试验期间最高平均气温为34.0 ℃，最低平均气温为25.0 ℃，隔1周采样1次。采用虹吸管抽取方式采集水样，分别于水面以下15和25 cm处采集等体积水样组成混合样(100 mL)。每次取样后添加去离子水以弥补蒸发、蒸腾以及取样所带来的水损耗。水样及试验用水总磷采用过硫酸钾氧化-钼锑抗分光光度法测定。

　　2 结果与分析

　　2.1 总磷去除效率的时间变化

　　随着处理时间的延长，2组植物净化系统总磷累计去除负荷呈逐渐增加趋势。至试验结束时(第49 天)，苦草组和金鱼藻组净化系统的累计去除负荷分别为1 151.17和1 167.50 mg·m−2，对应总磷去除率分别为82.8%和84.0%。方差分析表明：试验结束时2组植物系统总磷去除效率无显著差异(P>0.05)，但均显著好于无植物的对照组(P<0.05)。进一步对各周总磷去除负荷分析发现：2组植物系统各周之间的总磷去除负荷均存在较大波动，并且均在第3周出现峰值，单周贡献了整个试验期间(共7周)总磷累计去负荷的28.5%(苦草组)和16.7%(金鱼藻组)。无植物的对照组总磷去除负荷在第1周出现峰值。

　　在不同时间尺度上比较2组植物净化系统总磷去除效率发现：尽管在试验结束时(第49 天)和试验中期(第21 天和第28 天)两者差异不显著，但在试验初期(第7 天和第14 天)和中后期(第35 天和第42 天)2组植物净化系统总磷去除效率存在差异显著(P<0.05)。其中试验初期，金鱼藻组净化系统总磷去除效率显著高于苦草组净化系统;而当处理时间增至35和42 d时，苦草组又显著高于金鱼藻组。

　　这种不同沉水植物系统除磷效率的差异在不同时间尺度上出现偏离甚至逆转的现象可能具有一定的普通性，如文献[29]、文献[30]和文献[31]的研究中均存在这种现象。值得注意的是，忽视植物效应差异的不确性，可能会导致水生植物净化系统工艺优化，特别是水生植物的筛选结果出现偏离甚至失真，因此，有必要借助污染物去除动力学模拟分析获取用于指示和预测污染物去除效率的相关特征参数。

　　2.2 一级动力学和Grau二级动力学模拟与验证

　　分别运用一级动力学[公式(2)]和Grau二级动力学[公式(4)]模型对2组沉水植物净化系统试验数据进行拟合，并进行模型统计验证分析。一级动力学模型具有较高的拟合度(判定系数R2>0.930)，模型预测值与实测值之间吻合程度较好(相对均方根误差RRMSE<0.200)，表明一级动力学模型可用来描述沉水植物净化系统批式运行周期内总磷的去除过程。

　　2组净化系统回归求得的总磷去除一级动力学常数(Fa)分别为0.039 d−1(苦草组)和0.035 d−1(金鱼藻组)。不同于一级动力学，Grau二级动力学模型拟合度较差(R2<0.600)，表明该模型不适用于模拟本试验沉水植物净化系统总磷的去除过程，获取的动力学参数(n和m)也不具有现实意义。 一级动力学模型已被广泛应用于挺水植物型人工湿地系统中有机物、营养物、悬浮物等污染物去除效果的预测[19−22]，但该模型是否适用于模拟沉水植物型净化系统的污染物去除动力学至今还未得到证实。本研究发现：一级动力学模型可用于描述沉水植物型净化系统在批式模式下的总磷去除过程。

　　3 讨论

　　3.1 动力学模型的比较

　　污染物去除动力学是指导工程设计的重要基础。许多动力学模型已用于模拟污水生物处理和人工湿地系统污染物的去除过程，其中一级动力学、Grau二级动力学和Monod 动力学模型简化了污染物去除过程中所涉及的生物、化学、物理及其相互影响等复杂过程，模型构式可线性化处理而使得参数求解十分简单，拟合求得的动力学参数又可以较好地代表处理系统污染物的去除能力。

　　这些有利因素促使此类模型越来越多地应用到包括水平潜流[32]、垂直流[33]、潮汐流[34]、静止状态[35]及其组合工艺[19]等多种人工湿地系统中污染物的去除动力学特征的研究。但近年来研究发现：此类建立在众多假定条件基础上的表观动力学模型的适用性，不仅与湿地工艺类型(水平潜流和垂直流)及其水动力模式(推流或完全混合态)有关，还与所描述的污染物类型有关[19−20]。

　　本研究结果表明：一级动力学和Monod 动力学模型均可以用来描述沉水植物型净化系统批式试验周期内的总磷去除过程，而Grau二级动力学却不适用;相较于一级动力学，Monod 动力学模型具有更高的拟合优度。近年来，修正的Gompertz模型作为一种多参数非线性模型，在污染物去除动力学研究领域逐步被重视。

　　众多研究已发现：这种非线性预测模型描述异养硝化菌、厌氧氨氧化等新型好氧或厌氧生物处理系统的脱氮除磷效果具有较好的拟合精确度[23−25]。但该模型能否用于模拟人工湿地、植物塘等生态处理系统中污染物的去除过程至今还未见报道。本研究发现：修正的Gompertz模型可以很好地模拟沉水植物净化系统总磷的去除动力学过程。相较于一级动力学和Monod动力学模型，修正的Gompertz模型的拟合除了可以输出代表着污染物去除能力的关键性能参数(如最大体积去除速率)，还可以输出可用于指示前期阶段污染物去除的迟缓时间。这些参数的获取有助于对污染物去除机理的理解，也可用于协助工艺系统的优化设计和过程调控。

　　4 结论

　　苦草型和金鱼藻型2种植物净化系统总磷去除效率无显著差异，总磷去除率分别为82.8%(苦草组)和84.0%(金鱼藻组)，但两者总磷去除效率在不同时间尺度上的差异存在不确定性。4种动力学模型模拟与验证结果表明：除Grau二级动力学外，一级动力学、Monod 动力学模型和修正的Gompertz模型均可以有效描述2种沉水植物型净化系统试验条件下总磷的去除过程，其中Monod动力学和Gompertz模型的模拟具有更高的拟合优度，模型预测值与实测值之间吻合程度更好。

　　3种有效模型拟合获得的动力学常数在植物种类之间均存在显著差异(P<0.05)，其中指示除磷效率的动力学速率常数一致表现为苦草组大于金鱼藻组，表明苦草型净化系统对总磷的去除效率好于金鱼藻型。但苦草组除磷延迟时间显著长于鱼藻组(P<0.05)，表明其生境适应周期相对较长。

　　5 参考文献

　　[1]杨林章, 施卫明, 薛利红, 等. 农村面源污染治理的“4R”理论与工程实践——总体思路与“4R”治理技术[J].农业环境科学学报, 2013,32(1): 1 − 8.YANG Linzhang, SHI Weiming, XUE Lihong,etal. Reduce-retain-reuse-restore technology for the controlling theagricultural non-point source pollution in countryside in China: general countermeasures and technologies[J].JAgroEnvironSci, 2013,32(1): 1 − 8.

　　[2]武淑霞, 刘宏斌, 刘申, 等. 农业面源污染现状及防控技术[J]. 中国工程科学, 2018,20(5): 23 − 30.WU Shuxia, LIU Hongbin, LIU Shen,etal. Review of current situation of agricultural non-point source pollution and itspreventionand control technologies[J].EngSci, 2018,20(5): 23 − 30.

　　[3]翟敏婷, 辛卓航, 韩建旭, 等. 河流水质模拟及污染源归因分析[J].中国环境科学, 2019,39(8): 3457 − 3464.ZHAI Minting, XIN Zhuohang, HAN Jianxu,etal. Water quality simulation and multi-source attribution analysis[J].ChinaEnvironSci, 2019,39(8): 3457 − 3464.

　　作者：龚苗苗1，蔡飞翔2，姜培坤1，曹玉成1