随着经济社会的发展，部分水体富营养化情况严重，导致其不再适合饮用或水产养殖。据Liebig最小定律，磷常被视为水体富营养化与藻类爆发的主要诱导因素之一^[1-2]。在水体的生态修复中，作为一种有效替代传统水处理的方法，人工湿地系统由于其出水水质好、建设成本低、能耗低、易于操作及维护简单等优点在国内外得到广泛应用。人工湿地是主要利用基质、植物、微生物三重协同作用有效降低水体磷污染负荷^[3-4]的人工建造并控制运行的系统，其中，基质的吸附作用被认为是除磷的主要机制^[5]。同时，相关研究与工程实践发现，人工湿地基质除磷普遍存在去除效率不断下降的问题。不同基质吸附能力不同，有效运行时间不同，加之微生物作用对不同污染特征的水体除磷能力也各异，因此不同基质不同水质的基质除磷效率很难仅根据理论研究成果直接设计人工湿地基质结构，给基质结构优化设计与人工湿地维护带来不确定性。迄今为止，对人工湿地基质的吸附除磷过程与除磷效率问题，基质的静态吸附实验有大量研究报道^[6-7]，在实验室砂箱或野外实际工程尺度也有不少研究与报道^[8]，但将二者紧密结合，明确人工湿地不同构建基质的吸附除磷机制与主次及整体除磷效率，实现工程设计优化，达到更长期有效运行，这方面的研究却少有报道，尤其是针对缺少碳源的微污染含磷水体的净水除磷则更有必要深入研究。

为有效解决人工湿地基质除磷普遍存在的去除率不断下降问题及基质结构的优化设计与维护问题，需针对不同介质掌握其对具体含磷污水的吸附特性及实际去磷效果的演变趋势与规律。基于适用性、经济性、易得性原则，本文选取人工湿地可选用的除磷性能较好的1~2、2~4、4~8 mm这3种粒径的沸石、麦饭石、磁铁矿基质，首先通过静态实验研究其对磷的吸附动力学过程和等温吸附特征，并考察不同基质和粒径对吸附磷的影响。进而选取三者中磷吸附性能最好的沸石进行去磷的砂箱渗流模拟试验，以模拟获取潜流湿地的基质除磷净化效果。

1 材料与方法 1.1 实验基质

实验所用沸石、麦饭石和磁铁矿均来源于河南省巩义市，粒径分为1~2、2~4和4~8 mm共3种。3种基质组成成分如下：沸石(MgO 1.4%、Fe₂O₃ 1.8%、Al₂O₃ 14%)；麦饭石(MgO 0.65%、Fe₂O₃ 1.6%、Al₂O₃ 17.74%)；磁铁矿(Fe₂O₃ 4.5%、Al₂O₃ 0.32%)。

1.2 实验方法 1.2.1 基质对磷的吸附动力学静态实验

分别称取2~4 mm沸石、麦饭石和磁铁矿各1 g于50 mL离心管中，加入40 mL质量浓度为50 mg/L的KH₂PO₄溶液。分别在每个离心管中加入3滴氯仿防止微生物活动，然后置于温度为25 ℃、转速为175 r/min的恒温振荡器中连续震荡。时隔0.5、1、1.5、2、4、6、8、10、12、14、16、18、20、22、24 h后在5 000 r/min转速下离心5 min，采用钼酸铵分光光度法测定上清液中磷的浓度，所有实验都同时进行2个平行实验，以减少偶然误差。吸附量的计算如下(下同)：

$ Q=\frac{C_{0}-C_{\mathrm{e}}}{m} \times V, , $

(1)

其中：C₀为初始磷质量浓度(mg/L)；C_e为上清液中磷质量浓度(mg/L)；m为基质投加量(kg)；V为水样体积(L)。

1.2.2 基质对磷的等温吸附静态实验

分别称取3种粒径的3种基质(共9种样品)各1 g于50 mL离心管中，在每种样品中分别加入质量浓度梯度为1、5、10、20、40、100、150、200、250、300、350和400 mg/L(以P计，下同)的KH₂PO₄溶液40 mL。分别在每个离心管中加入3滴氯仿防止微生物活动，然后置于恒温摇床中以175 r/min、温度25 ℃震荡24 h。最后在5 000 r/min转速下离心，采用钼酸铵分光光度法测定上清液中磷的浓度，用式(1)计算吸附量。

1.2.3 粒径对基质吸附量影响的静态实验

分别称取3种粒径的3种基质(共9种样品)各1 g于50 mL离心管中，在每种样品中分别加入质量浓度为3 mg/L(以P计，下同)的KH₂PO₄溶液40 mL。分别在每个离心管中加入3滴氯仿防止微生物活动，然后置于温度为25 ℃、转速为175 r/min的恒温振荡器中连续震荡24 h。最后在5 000 r/min转速下离心，采用钼酸铵分光光度法测定上清液中磷的浓度，用式(1)计算吸附量。

1.2.4 基质除磷效果的砂箱物理模型渗流模拟试验

1) 渗流模拟试验概况

利用砂箱物理模型对基质除磷进行小试模拟试验研究，砂箱采用有机玻璃制成，长、宽、高分别为96、56、40 cm，其中，两边各有宽为8 cm的窄型水槽，高度分别为67、57 cm，见图 1。水槽隔板上均匀分布着直径为0.5 cm的布水孔，并固定2层30目的滤网，阻止中间部分的基质进入水槽。含水层中水流流速通过调整两端水槽水位差来实现，并在砂箱一侧设有16个取样口。在砂箱中自下而上依次填充10 cm厚的粒径为4~8、2~4、1~2、0.5~1 mm的天然沸石。

2) 试验步骤

注入自来水(不含磷)，待形成稳定的饱和带后，用蠕动泵以0.06 m³/d速率向左侧布水区连续输入3 mg/L的含磷溶液(实际人工湿地系统进水磷质量浓度一般在1~6 mg/L，该浓度有一定代表性)，砂箱水力梯度保持0.001 9，为稳定渗流场，保持水力停留时间(HRT)4层平均值为24 h，总磷表面负荷为0.4 g/(m²·d)(一般实际人工湿地设计负荷为0.3~0.5 g/(m²·d))。运行后前半个月每3 d取样一次，之后一个半月每周取样一次。

1.3 吸附除磷模型的初选

为了解基质对磷的吸附过程，本研究采用一级动力学方程、双常数方程和Elovich方程3种常用的吸附动力学模型对磷的吸附动力学进行拟合。

一级动力学方程可以描述溶液浓度与吸附时间之间的关系，简单表示为

$ \ln C_{t}=a-k t, $

(2)

适用于由扩散步骤控制的动力学过程。

双常数方程可以描述吸附量与吸附时间之间的关系，简单表示为

$ \ln Q=a+k \ln t, $

(3)

适用于非均相复杂反应的动力学过程。

Elovich方程可以描述吸附量与吸附时间之间的关系，简单表示为

$ Q=a+k \ln t, $

(4)

适用于反应过程中活化能变化较大和一些多界面的动力学过程。

以上各式中：C_t为t时刻溶液的浓度(mg/L)，t为反应时间(h), Q为吸附量(mg/kg)，a和k为动力学常数。

在固定的温度条件下，一定时间后，水中颗粒物对溶质的吸附达到平衡。本文采用Henry方程、Freundlich方程和Langmuir方程表征颗粒物表面的吸附量与溶液中溶质平衡浓度之间的关系。

Henry方程表达式为

$ Q=k c, $

(5)

Freundlich方程表达式为

$ Q=k c^{1 / n}, $

(6)

Langmuir方程表达式为

$ Q=Q^{0} c /(A+c), $

(7)

其中：Q为颗粒物表面的吸附量(mg/kg)；c为溶液中溶质的平衡浓度(mg/L)；Q⁰为单位表面上达到饱和时间的最大吸附量(mg/kg)；k、n、A为常数。

1.4 分析方法

吸附动力学过程与粒径对吸附量的影响采用Origin 8.0进行分析，渗流模拟试验中磷平均浓度的沿程变化采用ArcGIS 10.1中的kriging插值方法获得。

2 结果与讨论 2.1 吸附动力学模型分析 2.1.1 吸附动力学过程描述

根据各基质对磷的动态吸附实验结果绘制3种基质的吸附动力学曲线，见图 2。由图 2可以看出，各基质对磷的吸附均包括快速吸附、缓慢吸附和逐渐平衡3个过程。磁铁矿、麦饭石、沸石依次在10、12、16 h时达到平衡，之后随时间的吸附量波动较小，并且最大吸附量越大所需的饱和时间越长。在达到吸附饱和的时间内，前半段属于快速吸附阶段，后半段属于缓慢吸附阶段，这表明基质的表面都存在高中低能量的吸附位点^[9]。并且在相同时间内，3种基质的吸附量从大到小均表现出：沸石、麦饭石、磁铁矿。这表明，吸附量越大的基质与溶质的传质力越强。

2.1.2 基质吸附磷的动力学模型拟合情况

吸附动力学模型能够在一定程度上反映基质的吸附除磷机理，为研究吸附过程和分析吸附原理提供依据。吸附质在吸附剂多孔表面被吸附的过程分为4步：1)污染物从主体溶液扩散到由于水化作用在吸附剂表面形成的水膜表面；2)污染物克服液膜阻力并穿过液膜到达吸附剂表面；3)污染物从颗粒外表面扩散到颗粒内表面吸附位；4)污染物在活性位上发生吸附反应。

应用3种吸附动力学模型，对3种基质吸附除磷动力学过程进行分析，结果如表 1所示。从相关系数来看，双常数方程、Elovich方程均能较好地描述基质吸附磷的动力学过程，相关系数均达0.85以上；一级动力学方程拟合效果较差。这表明，基质对磷的吸附动力学并不是一个简单的一级反应，而是由反应速率和扩散因子综合控制的过程^[10]。比较而言，Elovich方程模型对3种基质吸附动力学过程的描述最为准确。在Elovich方程中，a表示初始吸附速率常数，k是与吸附剂表面覆盖程度及化学吸附活化能有关的参数，结合图 1来看，饱和吸附量越大，其初始吸附速率越大，基质表面影响化学吸附的活化能越大。

2.2 基质磷素吸附等温线拟合情况

吸附等温线在一定程度上反映吸附剂与吸附物的特性。从表 2可以看出，Langmuir方程能够更好地描述人工湿地基质对磷的吸附过程，拟合系数均达0.99以上。这也说明，基质吸附磷的过程更符合单分子层吸附理论，且颗粒表面较为均匀，即表面上各个晶格位置的吸附能力相同，这与相关的研究^[11]结论相一致。Freundlich方程能够较好地描述此吸附过程，拟合系数均达0.95以上。Henry方程对此吸附过程的刻画拟合相对较差，拟合系数相对于前两者来看相对较低。

Langmuir方程中，A为吸附量达到Q⁰/2时溶液的平衡溶度。并且可知，3种基质磷素理论饱和吸附量大小顺序依次为沸石(5 122.0 mg/kg)、麦饭石(4 169.1 mg/kg)、磁铁矿(3 126.3 mg/kg)。

在Henry型等温式中，k表示分配系数，即假定溶质在吸附剂与溶液之间按固定比值进行吸附。由此可知，同种粒径的情况下，沸石对磷的吸附比例最大，麦饭石次之，磁铁矿最小。

在Freundlich型等温式中，k值大致表示吸附能力的强弱，k值越大，对磷的吸附能力就越强^[12]。相同粒径时，3种基质k值大小依次为沸石、麦饭石、磁铁矿，这与吸附等温模型的结论一致。

2.3 基质粒径及基质成分组成对除磷效果的影响

3种不同粒径的沸石、麦饭石和磁铁矿对磷的吸附效果如图 3所示。3种基质均表现出，粒径为1~2 mm的基质对磷的吸附量最高，2~4 mm的次之，4~8 mm的最低，即随着基质粒径的减小，对磷的吸附量逐渐增大。基质粒径越小，其比表面积越大，与溶液中磷接触的机会越多^[13]，吸附质的传质速率越大^[14]，更有利于基质与磷发生吸附作用。相关研究也表明，基质的粒径与吸附量呈明显的负相关关系^{[2, 15-17]}。

特别需要指出的是，粒径对沸石除磷吸附量的影响较大，随着粒径从4~8 mm降低到1~2 mm，除磷吸附量从26.47 mg/kg提升到46.46 mg/kg；但值得注意的是，不同基质其影响有一定差异，相对而言，粒径的减小对麦饭石和磁铁矿的除磷吸附量提升较小，仅使其除磷吸附量分别提高16.57和9.32 mg/kg。此外，在相同粒径的情况下，3种基质对磷的吸附量大小表现一致，从大到小均为沸石、麦饭石、磁铁矿。这与基质理化性质相关，Mg、Fe和Al等元素是促进磷吸收的重要指标^[18]，通过形成Mg-P、Fe-P和Al-P体系降低磷的有效性^[19]。由3种基质中Mg、Fe和Al的含量可知，沸石中MgO的含量为1.4%，远远高于另外两种基质，沸石和麦饭石中Al₂O₃含量比磁铁矿中高出40多倍，依次为14%和17.74%。比较而言，Mg和Al所形成的结合体系对磷有效性的降低更为显著，Mg和Al含量的增加明显地提升了磷的去除效率。在相关的基质吸附磷的性能研究中^[20]，与本实验的同种基质吸附性能差别较大，但并未给出相关成分具体含量，因此，在选用基质时，对不同基质进行理化参数检测显得至关重要。

2.4 基质除磷砂箱物理模型渗流模拟试验结果分析与讨论

将实验第一个月所测得的每个取样点的磷质量浓度取平均值，并利用ArcGIS所做的砂箱磷平均浓度沿程变化如图 4所示。由图 4可知，在垂向上自下而上随着基质粒径的减小，磷浓度逐渐降低，即更多的磷素被从水中去除，且浓度分布呈较明显的层状分布，这与砂箱层状铺设的方式表现出较大的相关性，也再次证明基质的粒径与吸附量呈明显的负相关关系。水平方向上，磷浓度基本呈现推流式递减。

砂箱中各层(自上而下依次为第1到第4层)实际去除率如图 5所示。根据基质吸附静态实验结果，依据基质Langmuir吸附理论，计算对应进水磷质量浓度为3 mg/L的基质最大理论吸附量，根据各层水流通量及磷监测浓度，可获得各层已有最大吸附量，据此估算针对进水磷质量浓度保持为3 mg/L时底部3层不同粒度基质吸附除磷理论去除率(图中虚线)。从图 5可以看出，不同粒径的沸石基质吸附除磷及总体除磷效率呈现较大差异，但各粒径沸石基质的总体除磷效率具有波动现象，总体呈先上升再明显下降趋势。由于在无植物种植的人工湿地系统中磷的去除只来自于基质吸附与微生物或微生物膜除磷，并且静态实验中所体现的吸附效果为去除微生物干扰的基质本身的除磷特性，总体除磷效率明显高出理论吸附去除效率，以及初期的磷去除率升高及波状变化特征等，这些都表明，在无外加碳源的微污染含磷水体除磷净化中，除明显的吸附除磷机制，微生物或微生物膜除磷作用也可产生相当的去磷效果。即砂箱每层总去除率与基质吸附去除率差值代表微生物膜对去除率的贡献估值。由图 5可知，砂箱中磷的实际去除率总体上自第4层到第1层随着基质粒径的减小而逐渐增大，Freundlich方程中的k值代表对磷吸附能力的强弱，而从前面等温吸附实验结果看，4~8、2~4和1~2 mm沸石的k值依次为61.07、101.63和126.36，小粒径所引起的更强吸附性能在该渗流实验中也得到了较好印证。另外，各层基质在砂箱运行一段时间后，各层去除率均出现明显下降趋势，这一方面由于沸石对磷的吸附过程更倾向于单分子层吸附，废水中的磷不断地占据基质表面的吸附位点，使得可吸附磷的位置减少。

根据前期静态吸附实验结果，可以看出，在各层中随着时间的延长，微生物的作用均经历初期不断增加且不断波动现象，但3层中变化幅度不尽相同。小粒径层微生物作用相对贡献较大，最大颗粒的底层微生物作用很小，3个粒度基质层微生物机制磷去除率大致为5%~20%。这是由于粒径的增大使得微生物可附着的比表面更少，并且砂箱中渗透系数自顶层73 m/d增加至底层的222 m/d，下层流速明显快于上层，底层水力停留时间较短，微生物去除时间也不充分。因此，在人工湿地系统工程设计中最佳磷处理效果应更好发挥较小颗粒的高吸附性能及明显的微生物去除作用，同时也应注意到其渗透速度低带来的污水净化处理量小的不足，以及如何最大可能降低湿地系统易堵塞的风险。在一些实际人工湿地系统中，磷的去除效率大致为22%~67%^[21-24]，这为本研究提供了一定的参考依据。

3 结论

本研究通过静态实验与砂箱渗流模拟试验相结合，揭示不同类型及不同成分组成基质的除磷动力学过程和吸附机理，查明不同粒径的沸石基质渗流状态下实际去磷效果的演变趋势与规律。具体研究成果总结如下：

1) 沸石、麦饭石、磁铁矿等3种基质对磷的吸附包括快速吸附、缓慢吸附和逐渐平衡3个过程。磁铁矿、麦饭石、沸石依次在10、12和16 h时达到吸附平衡，并且最大吸附量越大所需的饱和时间越长。Elovich方程能够更为准确地描述基质吸附磷的动力学过程，3种基质的拟合系数均达到0.93以上。这表明，基质对磷的吸附动力学并不是一个简单的一级反应，而是由反应速率和扩散因子综合控制的过程。

2) Langmuir方程和Freundlich方程均能较好地描述人工湿地基质对磷的吸附过程，拟合系数均达到0.95以上。从Langmuir方程可知，3种基质磷素理论饱和吸附量大小顺序依次为沸石(5 122.0 mg/kg)、麦饭石(4 169.1 mg/kg)、磁铁矿(3 126.3 mg/kg)。

3) 基质粒径与除磷吸附量呈明显的负相关关系，粒径对沸石除磷吸附的影响较大，随着粒径从4~8 mm降低到1~2 mm，除磷吸附量提高约20 mg/kg，麦饭石和磁铁矿仅分别提高16.57和9.32 mg/kg。

4) 比拟潜流人工湿地系统基质部分，利用含磷低碳源微污染进水的沸石砂箱渗流模拟系统，模拟除磷过程变化规律与趋势。可以得到，不同粒径的沸石基质吸附除磷及总体除磷效率呈现较大差异，在一定粒度范围内，粒径小者除磷效率较大；各粒径沸石基质的总体除磷效率具有波动现象，且总体呈先上升再明显下降趋势。总体除磷效率明显高出理论吸附去除效率，同时考虑到模拟初期的磷去除率逐渐升高及波状变化特征等，这些都表明，在无外加碳源的微污染含磷水体除磷净化中，除明显的吸附除磷机制，微生物或微生物膜除磷作用也可产生相当的去磷效果。通过此种综合研究，可有效确定人工湿地基质除磷机制及其主次关系与基质整体除磷效率，有望为实际潜流人工湿地基质选取与结构优化设计及后期基质及时维护提供理论依据与工作方法借鉴。