随着我国农业的快速发展，以及工业废水和城市生活污水等点源污染得到有效控制，农业面源污染的重要性越来越突出。据生态环境部2010年发布的《全国第一次污染源普查公报》，我国农业面源污染物总氮和总磷的排放量占总排放量的57.2%和67.4%，已成为水体污染的主要来源。如何有效控制农业面源污染物向水体的排放，是当前学者关注和研究的热点和难点问题。

在多年研究的基础上，杨林章等^[1]系统总结提出了农业面源污染防控的4R技术体系“源头减量（Reduce）-过程阻断（Retain）-养分再利用（Reuse）-水体修复（Restore）”，并在全国范围内得到了推广应用。其中过程阻断是指在污染物向水体迁移过程中，通过一些物理的、生物的以及工程的方法等对污染物进行拦截阻断和强化净化，延长其在陆域的停留时间，最大化减少其进入水体的污染物量。目前应用比较普遍的有生态拦截沟渠技术^[2]、生物篱技术^[3]、生态拦截缓冲带技术^[4]等。其中，由传统农田排水沟渠发展而成的生态沟渠应用最为广泛，其由沟渠、基质和植物组成，具有生态性和生物性的特殊优点，并可根据地形特征来构建^[5]。生态沟渠不仅可以用来净化农田径流排水，还具有湿地生态系统的功能，在维持生态系统物种多样性、促进生态系统生产力和养分循环、净化水质等方面起着重要作用^[6-7]。其对径流污染物的去除净化主要通过减缓流速、植物吸收、基质吸附、泥沙沉降、微生物硝化反硝化等发挥作用^[8]。生态沟渠的规格以及内部构造组合等均会影响生态沟渠对污染物的削减效果，加上农业面源污染物浓度具有极大的时空变异性，因此亟需系统考察生态沟渠的规格及其构造组成对污染物削减效果及其稳定性、耐冲击负荷能力等的影响，寻求污染物去除效率高且稳的最佳生态沟渠。

为此，本研究利用太湖流域平原河网区常用的三种不同深度规格的生态沟渠，通过模拟农业面源污水，比较研究了相同表面水力负荷下三种沟渠对动态连续进水条件下铵态氮（NH₄⁺-N）、总氮（TN）、总磷（TP）和悬浮物（SS）的去除效果以及耐冲击负荷能力，并分析了进水浓度较高情况下，生态沟渠出水达到地表水环境质量Ⅴ类标准（GB3838-2002）所需长度，以期为生态拦截沟渠的应用设计（规格以及长度等）提供科学参考，从而提高对农业面源（尤其是种植业面源）污染物的拦截净化效果。

试验在上海市农业科学院庄行综合试验基地（30º53′N，121º23′E）进行。研究区为亚热带季风气候，多年平均降水量、蒸发量、气温分别为1 192 mm、1 237 mm、16.1℃。全年日照时数1 900 h，无霜期224 d。

设计目前生产上较常见的三种规格的生态沟渠，渠深0.80 m（E_0.80）、1.05 m（E_1.05）、1.30 m（E_1.30），长度各90 m，沟渠内均种植相同密度植被，沟壁为带孔预制板，孔内为狗牙根（Cynodon dactylon），沟底为常绿苦草（Vallisneria natans，购自上海海洋大学）（图 1）。试验开始前，生态沟渠稳定一段时间并对水量进行调试。

图 1 生态沟渠试验区 Fig. 1 Ecological ditch experimental site

于2013年7月23日开始动态进水试验，试验持续20 d，期间的气温和降雨情况见图 2。试验进水根据当地农田尾水中的氮磷浓度^[9]，用附近河道水加尿素和KH₂PO₄配置，进水TN 0.86~6.13 mg·L^-1、TP 0.11~0.28 mg·L^-1。沟渠断面规格及采样点断面设置见图 3，分别在进水处（1#）、22.5 m处（2#）、45 m处（3#）、67.5 m处（4#）及出水（5#）设置采样点。试验最初2 d每1 d采集水样一次，之后间隔1 d采样一次直至试验结束。

图 2 试验期间平均日温和降雨量变化 Fig. 2 Mean daily air temperature and rainfall during the experimental period

图 3 生态沟渠断面规格图与取样断面设置 Fig. 3 Sectional of the ecological ditch relative to specification and set-up of sampling section

水样铵态氮（NH₄⁺-N）采用纳式试剂分光光度法测定，总氮采用过硫酸钾氧化法测定，总磷采用钼锑抗分光光度法测定，悬浮物经过滤后，105℃烘干至恒质量后测定^[10]。

试验过程中，维持三种生态沟渠水位处于同一水平，通过控制不同规格沟渠的进水流量，控制三种规格生态沟渠的表面水力负荷保持一致。综合考虑进水水质分布特征、处理效率及经济性，借鉴表流湿地水力负荷后，确定沟渠表面水力负荷（N）为0.60 m³·m^-2·d^-1，在保证沟渠水位及相同的表面水力负荷下，计算三种规格沟渠的进水流量（Q）。计算方法如下：

沟渠有效容积V=S×L

沟渠水表面积S_A=N×L

进水流量Q=(N×S_A)/24 h

水力停留时间HRT=V/Q

式中，S为有效断面面积，L为沟渠长度，N为沟渠表面水力负荷。具体试验参数如表 1所示。

表 1 Table 1

表 1 不同规格生态沟渠实验基本设计参数表 Table 1 Parameters of basic design of the ditch relative to specifications 沟渠规格 Specification of ditches 表面水力负荷Surface hydraulic loading/ (m3·m-2·d-1) 水深 Water depth/m 长度 Length /m 有效断面面积Effective sectional area/m2 有效容积 Effective volume/m3 水表面积Water surface area/m2 进水流量 Inflow flux /(m3·h-1) 水力停留时间Hydraulic retention time/h E0.80 0.60 0.40 90 0.20 18.0 54.0 1.35 13.3 E1.05 0.60 0.55 90 0.37 32.9 65.3 1.63 20.2 E1.30 0.60 0.80 90 0.56 50.6 76.5 1.91 26.5

表 1 不同规格生态沟渠实验基本设计参数表 Table 1 Parameters of basic design of the ditch relative to specifications

采用污染物去除率的变异系数来表征污染物去除效果的稳定性：

变异系数=标准偏差/平均值×100%

采用污染物表面负荷与污染物面积去除量相关系数来表征耐冲击负荷能力，其中污染物表面负荷表示沟渠单位面积每天接受的进水中污染物的量：

污染物表面负荷（g·m^-2·d^-1）=C_in×N

式中，C_in为进水浓度，N为沟渠表面水力负荷。

污染物面积去除量（g·m^-2·d^-1）=（C_in-C_out）×Q×24 h/S_A

式中，C_in为进水浓度，C_out为出水浓度，Q为进水流量，S_A为沟渠水表面积。

采用SPSS 13.0软件进行统计分析，采用Sigmaplot 12.0软件进行制图。

试验期间生态沟渠进水NH₄⁺-N浓度为0.17~1.23 mg·L^-1，平均0.67 mg·L^-1；E_0.80、E_1.05、E_1.30三种沟渠的最终出水NH₄⁺-N浓度分别为0.02~0.49 mg·L^-1，0.05~1.00 mg·L^-1和0.07~0.33 mg·L^-1，平均总去除率分别为61.4%、54.7%、64.8%，可见E_1.30处理效率相对最高，而且出水浓度比较稳定（图 4，表 2）。三种规格生态沟渠对TN的去除效果也以E_1.30的最好，出水TN浓度平均0.66 mg·L^-1，低于E_0.80（平均0.74 mg·L^-1）和E_1.05（平均0.85 mg·L^-1），平均去除率63.1%（44.7~87.1%），也高于E_0.80（58.7%）和E_1.05（62.5%），而E_1.05和E_1.30的TN去除效果均比较稳定，去除率的变异系数分别为21.9%和23.1%。

图 4 铵态氮（NH4+-N）、总氮（TN）、总磷（TP）和悬浮物（SS）浓度变化情况 Fig. 4 Dynamic change in concentration of NH4+-N, TN, TP and SS in ecological ditch relative to specification

表 2 Table 2

表 2 不同规格生态沟渠对面源污水中污染物的去除效率和变异系数 Table 2 Removal efficiency of pollutants form non-point source pollution and coefficient variation relative to specification of the ditch 沟渠规格 Specification of ditches NH4+-N TN TP SS 去除效率 Removal efficiency/% 变异系数 Coefficient variation/% 去除效率 Removal efficiency/% 变异系数 Coefficient variation/% 去除效率 Removal efficiency/% 变异系数 Coefficient variation/% 去除效率 Removal efficiency/% 变异系数 Coefficient variation/% E0.80 61.4±29.0 47.2 58.7±18.6 31.7 64.0±16.3 25.5 58.5±26.7 45.7 E1.05 54.7±31.2 57.1 62.5±13.7 21.9 70.2±13.8 19.7 50.3±36.6 73.7 E1.30 64.9±14.2 21.9 63.1±14.6 23.1 71.8±8.6 11.9 60.8±21.8 35.8

表 2 不同规格生态沟渠对面源污水中污染物的去除效率和变异系数 Table 2 Removal efficiency of pollutants form non-point source pollution and coefficient variation relative to specification of the ditch

三种规格生态沟渠对TP的去除效果要普遍高于NH₄⁺-N和TN，平均去除率分别为64.0%（E_0.80）、70.2%（E_1.05）和71.8%（E_1.30），以E_1.30最高，且去除效果相对最为稳定，去除率的变异系数最小，为11.9%，出水浓度变化在0.02~0.10 mg·L^-1之间。

与NH₄⁺-N和TN结果类似，三种规格生态沟渠对SS的去除效果整体较好，SS总去除率分别为E_0.80平均58.5%，E_1.05平均50.3%，E_1.30平均60.8%；E_0.80、E_1.05和E_1.30最终出水的SS平均浓度分别为20.9 mg·L^-1（2.0~52.0 mg·L^-1）、23.1 mg·L^-1（6.0~56.0 mg·L^-1）和18.1 mg·L^-1（2.0~40.0 mg·L^-1），远低于进水SS浓度的49.1 mg·L^-1（16.0~70.0 mg·L^-1）。整体上，E_1.30对SS的处理效率相对较高且较为稳定。

生态沟渠主要拦截农田尾水中的氮、磷，因此，本研究仅对进水浓度较高情况下（TN≥2mg·L^-1，TP≥0.2mg·L^-1），TN、TP沿生态沟渠长度的浓度变化情况进行分析，以期探索适宜建设长度。对不同规格生态沟渠长度与TN、TP浓度的变化进行拟合可以看出，随着长度的增加，三种规格沟渠中TN、TP浓度均逐渐降低（图 5）。E_0.80和E_1.05的TN浓度与沟渠长度之间的递减规律较好，拟合方程相关系数分别为0.981 1和0.991 6。三种规格沟渠的TP浓度与长度之间的递减规律均较好，E_0.80、E_1.05和E_1.30的相关系数分别到达1、0.983 8和0.967 8。虽然E_1.30沿程上的TN、TP浓度波动相对较大，但是在进水浓度一致的情况下，该沟渠在22.5 m和45 m时的出水TN、TP平均浓度最低，说明短距离的E_1.30生态沟渠对N、P就有较好的处理效果。

图 5 不同规格生态沟渠中总氮（TN）和总磷（TP）沿程浓度变化 Fig. 5 Dynamic change in concentration of TN and TP along the ditch relative to specifications

将TN、TP出水浓度统一设置为地表水环境质量Ⅴ类标准（GB3838-2002），即：TN为2 mg·L^-1，TP为0.2 mg·L^-1，根据拟合方程可以计算出达标所需沟渠的适宜长度（表 3），综合而言，在进水浓度统一的情况下，E_1.30处理TN、TP所需长度与E_0.80和E_1.05相比较短，分别为27.4 m和4.9 m；其次为E_0.80，处理TN、TP所需长度分别为54.5 m和9.8 m。

表 3 Table 3

表 3 不同规格生态沟渠TN、TP达到地表水Ⅴ类标准的适宜长度 Table 3 Optimal length of the ditch to reduce TN and TP to Class Ⅴ standard for surface water 污染物 Pollutant 达到Ⅴ类标准所需长度Optimum length/m E0.80 E1.05 E1.30 TN 54.5 60.8 27.4 TP 9.8 10.5 4.9

表 3 不同规格生态沟渠TN、TP达到地表水Ⅴ类标准的适宜长度 Table 3 Optimal length of the ditch to reduce TN and TP to Class Ⅴ standard for surface water

由于表面水力负荷相同，因此在同一进水浓度的情况下，三种规格生态沟渠污染物表面负荷也相同。试验期间，NH₄⁺-N进水表面负荷为0.10~0.71 g·m^-2·d^-1；TN为0.52~3.68 g·m^-2·d^-1；TP为0.07~0.17 g·m^-2·d^-1；SS为9.60~42.00g·m^-2·d^-1。

从各污染物表面负荷与三种规格生态沟渠的污染物面积去除量关系来看（图 6），二者均显著相关（P < 0.05），三种规格生态沟渠对NH₄⁺-N和TN的处理能力均较好，E_1.30的耐冲击负荷能力相对较强，NH₄⁺-N和TN的相关系数分别为0.973和0.997，E_1.05的NH₄⁺-N和TN相关系数最小，耐冲击负荷能力相对较弱。而三种规格生态沟渠对TP的处理能力也均较好，E_1.05对TP的耐冲击负荷能力最强（r=0.960），E_1.30其次（r=0.876）。相比较而言，E_1.30对SS的处理能力较好，耐冲击负荷能力也最强（r=0.740）。

图 6 不同规格生态沟渠中铵态氮（a）、总氮（b）、总磷（c）和悬浮物（d）表面负荷与面积去除量的关系 Fig. 6 Relationships between surface loading and per unit area removalrate of NH4+-N (a), TN (b), TP (c) and SS (d) in the ditch relative to specification

由于氮、磷等营养物质的输入，农业面源污染已成为水体富营养化的主要污染源，生态沟渠技术是当前研究的控制技术之一^[11]。研究表明，生态沟渠对农田排水中氮磷拦截效果明显优于传统的混凝土沟渠和土质沟渠^[12]。但是生态沟渠对氮磷的去除效果具有较大的时空变异性，有学者利用生态沟渠对稻田径流水中污染物进行拦截，研究表明其对NH₄⁺-N、TN和TP的净化效率为31.81%、54.18%和58.21%^[13]；还有学者在灌溉和降雨条件下对生态沟渠的氮磷输出特征进行了研究，其对TN、TP去除率达64%和70%^[14]，造成氮磷去除效果不同的原因可能是由于生态沟渠规格不同引起的，沟渠的大小、边坡形式、断面尺寸、水力半径、纵坡等几何尺度均影响着水体中氮磷的转化和去除^[15]，如上述研究中的生态沟渠前者深为1.20 m，上下宽度分别为1.5 m和0.5 m，后者深为1.00 m，上下部宽为4 m和2.8 m。此外，研究表明，进水浓度对沟渠氮磷去除率也有不同程度的影响^[12]。因此，加强对同一条件下、不同规格生态沟渠处理效率的研究和评价，能为生态沟渠在农业面源污染控制上的应用提供技术支持。本研究在进水浓度一致、表面水力负荷相同条件下，比较研究了三种不同规格生态沟渠对农业面源中主要污染物的去除效果，研究结果表明，在动态条件下，三种规格生态沟渠对污染物的去除效果均较好，E_0.80对NH₄⁺-N、TN、TP和SS的平均去除率分别为61.4%、58.7%、64.0%和58.5%；E_1.05分别为54.7%、62.5%、70.2%和51.1%；E_1.30的平均去除效率与前两种规格相比较高，分别达到64.8%、63.1%、71.8%和60.8%。

生态沟渠作为一种特殊的湿地系统，对氮的去除主要是通过微生物的氨化、硝化与反硝化途径^[16]；此外，污染负荷、水力负荷和水力停留时间等也会影响生态沟渠的出水水质，污染负荷和水力负荷越小，水力停留时间越长，则出水水质越好^[17]。在本研究中，E_1.30对NH₄⁺-N和TN的去除效率均相对较高，其出水浓度也较为稳定，与其他两种规格生态沟渠相比，在相同表面水力负荷下，E_1.30的水力停留时间较长（26.5 h，表 1），这可能是其氮素去除率较高的原因之一^[18]。而虽然水力停留时间较短，E_0.80的NH₄⁺-N去除效率较E_1.05高，可能是由于植物等其他原因的影响，由于本文结果没有涉及植物的作用，故在此不做讨论。基质和底泥的吸附是生态沟渠中除磷的主要途径^[19]。在本研究中，E_1.30和E_1.05对TP的去除效率均相对较高，可能是相对于窄浅型生态沟渠，宽深型沟渠内的基质和底泥与水体接触面积更大，有利于其对磷的吸收。泥沙的截控固持能有效降低农田排水中的颗粒物浓度，也是控制农业面源污染的关键之一。在动态试验中，水力停留时间作为最直观的因素，揭示了生态沟渠对农田排水中颗粒物高效沉降的原因，它既是颗粒物的沉降时间，又是水体中污染物与植物、基质、底泥的接触时间^[8]。在本研究中，E_1.30对水体SS的去除效率相对较高，可能是较长的水力停留时间降低了SS输出能力，促进了SS沉淀。而E_0.80的SS去除效率较E_1.05高，可能是该类型沟渠中植物的作用超过了停留时间的影响。总体而言，与其他两种规格生态沟渠相比，E_1.30较长的水力停留时间、与基质和底泥较大的接触面积和较强的耐冲击负荷能力，促进了污染物的去除，使得该规格生态沟渠对农田排水中污染物的处理能力较好。

从主要污染物沿生态沟渠长度的浓度变化可以看出，随着沟渠长度增加，污染物浓度逐渐降低，说明污染物在生态沟渠中的迁移具有规律性，这与余红兵等^[14]的研究结果相似。也有学者对长沟渠的研究结果表明，在铺设盘培牧草的600 m沟渠中，各种污染物指标在靠近进水口处的300 m沟渠内降解速度较快，而在远离进水口处的后300 m沟渠内变化污染物变化幅度相对比较平缓^[20]，说明沟渠前段对污染物的处理效果较好，后段由于污染物浓度已经较低，处理效果并不明显。本研究通过2阶多项式方程计算出在进水浓度较高情况下，出水TN、TP浓度达到地表水Ⅴ类标准时，不同规格生态沟渠处理污染物的适宜长度，综合TN、TP的处理效果来看，处理TN所需沟渠长度与TP相比较长，但E_1.30所需长度与E_0.80和E_1.05相比较短，为27.4 m，而E_0.80和E_1.05分别为54.5 m和60.8 m；E_1.30处理TP所需沟渠长度也最短（4.9 m）。在长江中下游地区（如太湖流域）水网密布，农村经济发达，化肥施用量高，因此，在这一地区适宜利用生态沟渠来强化净化农业面源污染物，减轻水体污染负荷。本文研究的三种规格生态沟渠中，虽然E_1.30边坡建设所需材料较多，但处理污染物所需长度最短，若按27.4 m进行建造，占地面积为28.8 m²；而E_0.80虽然本身占地面积小，但处理污染物所需长度较长，若按54.5 m进行建造，则占地面积为43.6 m²。因此E_1.30总占地面积少，可以在一定程度上节约土地资源。从环境效益上看，按照水稻每667 m²耗水1500 m³计算^[21]，建设E_1.30生态沟渠工程从农田径流中减少污染物排放氮2.12 kg、磷0.20 kg（TN平均进水浓度2.24 mg·L^-1，去除率63.1%；TP平均进水浓度0.19 mg·L^-1，去除率71.8%），可有效减少对周边水体环境的危害。本研究主要关注的是相同表面水力负荷下不同规格生态沟渠对排水中污染物的去除效果，对于其脱氮除磷的内在机理及各去除途径的研究涉及较少，且在1.30 m深的生态沟渠中，如何在不增加规模的情况下，通过改变沟渠线型空间结构（间隔设置小型跌水堰等）、间隔设置基质（沸石基质等）、优化植物配置（沟壁和沟底植物群落）等方式，进一步提高生态沟渠处理效率、节约占地等，还有待于进一步研究。

在相同表面水力负荷下，三种深度规格生态沟渠对农业面源污染物均有较好的去除效果。与E_0.80和E_1.05相比，E_1.30对农业面源污染物NH₄⁺-N、TN、TP和SS的去除效率均相对较高，其出水浓度也较为稳定，且耐冲击负荷能力相对较强。在进水浓度较高条件下，出水TN、TP浓度达到地表水Ⅴ类标准时，E_1.30所需长度最短。综合而言来看，E_1.30对农田排水中主要污染物的去除效率、耐冲击负荷能力最好，在实践应用中应因地制宜根据需要应用或进一步提升该规格生态沟渠的设计，以达到对农业面源污染控制的最大化。