经济快速增长背景下，人们对生产生活的空间需求日益增加，长江中下游的巢湖、太湖、洞庭湖和鄱阳湖等滨湖区域开展了大规模的围湖造田活动^[1]。围垦区是滨湖低地生态系统的一种主要地理单元，在荷兰和中国，围垦区是滨湖区主要农业分布区，XIE等^[2]调查发现长江中下游围垦导致湖泊面积减少13.80%，鄱阳湖滨湖围垦区面积达到4 180 km^2[3]。

围垦区是一个较为封闭的生产生活区域，通常采用圈圩筑堤的方式建设，其水文管理需要通过电排灌站调控。围垦区中密集的农业生产活动产生大量的氮磷等营养物质，通过人工设施(如涵洞、电排灌站等)被排放至受纳水体，导致周边水生态环境质量的改变^[4]。因此，研究人员认为围垦区的排水量及排水形式对受纳水体的影响值得关注^[5]。YAN等^[6]模拟了太湖流域低地圩田与泵站的水文过程与排水磷负荷的关系，发现与自然排水相比，电排灌站管控下的年均排水量减少了8.6%。电排灌站排水运行初时会引起沟渠表层沉积物再悬浮，导致沟渠水中TP和浊度骤然上升^[7]。HUANG等^[8]利用氮动态模型对太湖流域低地圩田和非圩田进行了模拟，发现圩田排放的氮负荷高于非圩田。汛期电排灌站排水是引起下游国控断面TP超标的主要原因^[9]，其已成为受纳水体汛期水质超标的潜在风险源^[10]。此外，围垦区由于人类活动密集，过量的营养物质输出导致受纳水体富营养化^[11]。因此，持续关注低地圩田\围垦区内人为控制电排灌站排水水量水质特征，有利于精细刻画其对受纳水体的影响并提出相对有效的管控措施。

鄱阳湖是长江流域最大的“吞吐型”通江湖泊，入湖面源污染及其控制日益得到关注。已有研究^[12]表明，赣江、抚河等五河入湖水质中的TN和TP是引起鄱阳湖水质下降的主要因子，模型模拟估算结果表明农业源和城镇生活源的TP入湖负荷贡献率分别为56.4%和30.6%^[13]。本研究以鄱阳湖流域南昌湖区典型围垦区——蒋巷联圩为主要研究区，该区域地势低洼平坦，是四面环水的独立岛镇，全域排灌水均通过联圩内电排灌站调控。收集与整理电排灌站2021年、2023年全年运维数据，并于2023年开展现场水质监测调研，明确围垦区电排灌站排水的水量水质时空分布特征，并估算全域排水污染负荷，以期为鄱阳湖围垦区农业面源污染防控提供参考。

1.   材料与方法

1.1   研究区概况

蒋巷联圩地处赣江三角洲下游，是赣江南支和中支所夹区域，属于典型赣江-鄱阳湖河湖交错带围垦区^[14](图1)。全域总面积241 km²，根据土地利用类型及人口分布，可将全域分为上、中、下3段，其中上段面积为76.71 km²，人口约7.49×10⁴人，农田面积为23.14 km²；中段面积为67.33 km²，人口约1.84×10⁴人，农田面积为24.78 km²；下段为黄湖蓄滞洪区，面积为96.96 km²，人口约0.14×10⁴人，农田面积为55.61 km²。蒋巷联圩内共建有16座电排灌站，排水去向是赣江南支(9座)、赣江中支(5座)，其中高梧和山尾仅灌不排。上段共有高梧、叶娄、后河、九成和山尾5个电排灌站；中段共有联圩、三斗角、胜利新和胜利4个电排灌站；下段共有五丰、翻身、西舍、澎泽湖、玉丰、石头河和东沙湖7个电排灌站。

图 1  蒋巷联圩电排灌站位置分布图

Figure 1.  Map of electric drainage and irrigation stations in Jiangxiang Lianwei

下载: 全尺寸图片 幻灯片

1.2   数据来源

电排灌站排水量数据来源于2021年和2023年完整的电排灌站的运维资料，2022年因特大干旱，记录不完整。根据式(1)核算该电排灌站的年度总排放量。

式中：Q为年度总排放量，m³；t运行时间，h；n该电排灌站开泵台数；q结合泵站流量，m³·h⁻¹。

降雨量数据来源于中国气象科学数据共享服务网(气象数据共享网(mlogcn.com))南昌气象站2020—2023年逐日数据，并按照《降雨量等级(GB/T 28592—2012)》划分为小雨(0.1~9.9 mm)、中雨(10~24.9 mm)、大雨(25~49.9 mm)和暴雨(>50 mm)4种类型。前期研究发现，退水期、枯水期是鄱阳湖水质超标的主要时期，其中典型国控断面——南矶山断面在枯水期TP超标达2.28倍^[15]。蒋巷联圩围垦区是南矶山断面的主要汇水区域，电排灌站排水水质直接影响南矶山断面水环境质量，因此，本研究中电排灌站排水水量及水质分析将按照鄱阳湖水期进行划分，即涨水期(4—6月)、洪水期(7—9月)、退水期(10—11月)和枯水期(12—翌年3月)。

1.3   样品采集与水质指标分析

基于运维记录整理和分析发现，蒋巷联圩电排灌站排水发生在4—9月份，即涨水期和洪水期。为明确降雨和人类活动强度对电排灌站排水水量和水质的影响，根据2021年总排水量分析，分别在蒋巷联圩上、中、下段选择排水量最大、排水去向分别为赣江南支、中支的2个电排灌站进行排水水质样品采集，其中上段为九成和后河；中段选择三斗角和联圩；下段选择玉丰和西舍(图1)。采样时间为2023年4—9月，每月采集1次排水水样，采用干净的聚乙烯塑料瓶(1 L)于电排灌站排口处各采集1瓶，共采集36个样品。采集现场使用便携式浊度计(WZB-175，中国)测定浊度后，样品当天带回实验室且于4 ℃保存，并在7 d内完成相关指标测定。水质测定指标包括化学需氧量(COD)、总氮(TN)、氨氮(NH₄⁺-N)、硝态氮(${\text{NO}}_{\text{3}}^{{-}}$-N)、总磷(TP)、溶解性总磷(DTP)和颗粒态磷(PP)；有机氮由总氮减去硝态氮、氨氮和亚硝酸盐氮得到，其中亚硝酸盐氮含量低且不稳，易转化为其他形态的氮，故忽略不计^[16]；PP由TP减去DTP得到。COD采用分光光度法(DR2800，美国HACH)测定，其余指标均采用国标法进行测定^[17]；溶解态水质指标测定前需将水样经0.45 μm滤膜(Millipore，USA)过滤后测定。以超纯水为空白，每个采样点取3个平行样品进行分析，测得结果进行平均值计算，误差限在5%以内。

采用SPSS 23.0(IBM，USA)进行数据统计分析，Origin 2023(OriginLab，USA)和ArcMap10.8(ESRI，USA)进行可视化表达。

1.4   污染物排放负荷估算

经调研发现，围垦区地表径流主要通过电排灌站进行内外水力交换，且排放时间集中在涨水期和洪水期。基于区域独特的水文特性及文献调研^[18-19]，围垦区污染负荷排放主要受排水的排放量及污染物质量浓度影响。故本研究的电排灌站排放污染负荷根据式(2)估算。

式中：W为围垦区电排灌站灌排的污染物负荷量，t·a⁻¹；Q为排水量，m³·a⁻¹；C_i为第i种污染物质量浓度，mg·L⁻¹。

2.   结果与讨论

2.1   电排灌站的排水量时空分布特征

1)降雨量变化特征。土壤N、P流失特征与降雨量显著相关，较大降雨量条件下对土壤的冲刷力度加大，易携带较多N、P^[20]。平原圩区水网密布，水系联通，降雨是影响电排灌站排水量和水质的重要因子^[21]。由图2(a)可知，2020—2023年的年降雨量分别为2 096.2、1 856.7、1 559.0和1 658.4 mm，其中2022年降雨量最低，长江中下游区域遭遇罕见的大旱，鄱阳湖流域出现了1949年以来最严重的干旱^[22]。由图2(b)可知，小雨频次最高，年降雨量的平均贡献为13.98%；而暴雨频次最低，但年降雨量的平均贡献占33.65%，其中2022年暴雨贡献占比最低(24.21%)，2023年最高(41.84%)。在各个水期中，小雨虽然发生频次最高，但其对各个水期的降雨量影响不显著(P>0.05)。暴雨是各水期降雨量增加的主要降雨类型，尤其在涨水期和洪水期，例如，2020年洪水期暴雨发生5次，降雨量占比为51.10%；2021年涨水期暴雨发生了6次，占该水期降雨量的54.58%；2022年涨水期暴雨发生了4次，降雨量占比为41.67%；2023年涨水期暴雨发生了5次，降雨量占比为48.92%

图 2  不同水期下南昌市2020—2023年降雨量变化

Figure 2.  Changes in rainfall of Nanchang City from 2020 to 2023 under different water periods of Poyang Lake

下载: 全尺寸图片 幻灯片

2)排水量时空分布特征。由图3(a)可知，2021和2023年蒋巷联圩电排灌站总排水量分别为1.08×10⁸ m³和2.69×10⁷ m³，其中2023年排水量降低了75.04%，这可能与2022年极端干旱天气影响下，地下水储量下降有关^[23]。曹思佳等^[24]研究发现2022年鄱阳湖极端干旱导致地下水迅速外泄，是正常年份地下水排泄量的14.5倍。2023年总降雨量虽与2021年相差不大，但基于2022年的极端干旱，2023年降雨在较大程度上用于地下水补给，导致2023年电排灌站排水量大幅下降。由图3(b)可知，涨水期电排灌站排水量大于洪水期，占总排水量的63.36%(2021年)和89.08%(2023年)。鄱阳湖汛期气候变化对径流深度变化的影响相较于非汛期显著增加^[25]，同时长江中下游地区洪涝集中在5—7月，受涝次数占全年的80%^[26]，因此，受汛期降雨影响，涨水期电排灌站排水量占比增加。

图 3  电排灌站排水量时空分布

Figure 3.  Spatiotemporal distribution of drainage water from electric drainage and irrigation stations

下载: 全尺寸图片 幻灯片

在排水方向上，由于赣江南支的电排灌站布设数量大于赣江中支，因此赣江南支排水方向的排水量大于赣江中支，其中南支排水量占总排水量的54.82%(2021年)和53.81%(2023年)。在排水区域分布上，电排灌站排水量为中段>下段>上段。这与实际调研发现蒋巷联圩中段地势最低且分布有大面积水域有关（图1）。Pearson分析结果表明，电排灌站的排水量与降雨量和降雨类型呈显著正相关(P<0.05)，暴雨雨型下电排灌站日均排水量最大，为6.33×10⁵ m³。有研究^[27]表明，在下垫面类型不变的条件下，降雨强度越大，则地表径流形成时间越短，径流强度越大。本研究结果表明暴雨后的7 d内，围垦区内电排灌站排水总量占年总排水量的61.99%。

2.2   电排灌站排水水质的时空分布特征

1)排水水质的时间分布特征。如图4所示，2023年电排灌站排水的TN、TP质量浓度均值和COD值分别为(3.22±2.73)、(0.22±0.13)和(16.92±8.98) mg·L⁻¹，超过地表水环境质量标准(GB 3838-2002)III类限值，超标倍数均值分别为2.23、0.29和0.12。HUANG等^[28]发现太湖低地圩区排水TP质量浓度均值0.14 mg·L⁻¹；储茵等^[29]调研巢湖低地圩区排水TN、TP质量浓度分别为0.8~6.5 mg·L⁻¹、0.05~0.6 mg·L⁻¹；HUA等^[30]调研发现汉江平原农业地表径流TP质量浓度均值为0.14 mg·L⁻¹。太湖、巢湖圩区、汉江平原和鄱阳湖流域位于长江中下游地区，属于亚热带季风气候，降雨量差异较小；土地利用类型以农田为主，主要作物类型为水稻，研究发现地表径流氮磷浓度相似。

图 4  电排灌站排水的主要水质参数

Figure 4.  Drainage water quality of electric drainage and irrigation stations

下载: 全尺寸图片 幻灯片

不同水期下电排灌站排水水质分析结果表明(图4)，涨水期排水TN、TP质量浓度高于洪水期；涨水期排水TN、TP质量浓度均值和COD值分别为(4.03±1.61)、(0.23±0.15)和(16.36±9.92) mg·L⁻¹，其超标倍数均值分别是3.03、0.31和0.09；洪水期排水TN、TP质量浓度均值和COD值分别为(2.42±1.61)、(0.20±0.15)和(17.47±9.92) mg·L⁻¹，其超标倍数均值分别是1.43、0.26和0.15，这与围垦区农事活动有关^[31]。现场调研结果表明，围垦区沟渠、坑塘密集，在4—6月(涨水期)的早稻生长季，稻田分别施加尿素和有机氮肥作为基肥和分蘖期追肥，随着高频次暴雨和大雨(图2b)导致溢流出来的田面水汇入沟渠和坑塘中，水体氮以有机氮为主^[32]，洪水期围垦区内沟渠和坑塘的水体停留时间长，沟渠植被硝化速率高^[33]，排水水体以${\text{NO}}_{\text{3}}^{-}$-N为主。因此，涨水期排水TN质量浓度高于洪水期，且涨水期以有机氮为主，而洪水期以无机氮为主，其中无机氮以${\text{NO}}_{\text{3}}^{-}$-N为主(59.6%)。不同土地利用类型和水体污染特征在氮形态上的差异，农业低地集水区，电排灌站排水以${\text{NO}}_{\text{3}}^{-}$-N主要形式^[34]，而城市以NH₄⁺-N为主^[35]。受降雨驱动影响下的地表径流水体磷形态以PP为主^[36]，电排灌站排水PP占TP的60.67%，且与浊度呈显著正相关性(P<0.05)。涨水期电排灌站排涝运行时间长，引起水体扰动大，排水沟渠中底泥再悬浮，导致涨水期浊度高于洪水期(图4c)，高浊度与高TP浓度密切相关^[37]；鄱阳湖南昌湖区南矶山断面汇入区域研究发现TP浓度与浊度呈显著正相关（P<0.05），且PP占比超60%^[38]。

2)排水水质的空间分布特征。如图5所示，排水TN、TP平均质量浓度和浊度为上段>下段>中段，排水${\text{NO}}_{\text{3}}^{-}$-N、NH₄⁺-N平均质量浓度为上段>中段>下段，排水COD值、DTP平均质量浓度为下段>中段>上段。上段排水TN、TP平均质量浓度最高，分别为(3.45±2.09)、(0.26±0.17) mg·L⁻¹，分别超标2.3倍、0.45倍；下段排水COD值最高(19.64±8.9) mg·L⁻¹，超标倍数均值为0.18。人类活动和农田种植面积是影响排水水质空间分布的重要因素，人类活动对鄱阳湖水质产生的影响不容忽视^[39-40]；上段、中段和下段的人口分别占78.73%、19.78%和1.48%，其农田面积占各自区域的30.17%、36.81%和57.34%。蒋巷联圩上段以集镇为主，人口稠密，生活污水处理设施有待完善，未经处理的污水直排入中心沟渠，最终通过电排灌站排放至受纳水体，因此，导致排放水体中TN、TP浓度最高。

图 5  电排灌站排水主要污染物的空间变化特征

Figure 5.  Spatial variation characteristics of drainage pollutants from electric drainage and irrigation stations

下载: 全尺寸图片 幻灯片

在排水方向上，排往南支的排水TN、TP质量浓度均值和COD值高于排往中支，其值分别为(3.41±2.9)、(19.18±9.29)、(0.21±0.15) mg·L⁻¹，超标倍数均值分别为2.03、0.17、0.28。王朔月等^[41]的研究表明，2017—2018年赣江上游流向鄱阳湖湖区的氮磷浓度呈阶梯状递增，其中赣江(下游)污染物浓度显著高于香溪河和架竹河(上游)。刘文强等^[42]研究赣江中支周坊断面TP超标规律，发现其主要发生在涨水期，而电排灌站在涨水期的大量排水加剧了断面TP超标的风险。

2.3   电排灌站排水的主要污染物负荷量

由表1可知，2023年蒋巷联圩电排灌站年排放的TN、TP和COD负荷分别为111.07、6.14和501.95 t·a⁻¹，其中涨水期排水TN、TP和COD负荷分别占全年93.01%、90.20%和86.57%。从区域分布看，中段排放的TN、TP和COD负荷量最高，分别占全区域47.50%、52.67%和49.13%；虽然中段排水的氮磷质量浓度低(图5(a)和图5(b))，但其排水量最大(图3(b))，占比超过年总排水量50%以上。从电排灌站排水去向看，排往南支方向的TN、TP和COD负荷大于排往中支方向，分别占总负荷的51.57%、57.69%和56.89%。围垦区TN负荷年排放量是TP负荷的18.07倍，这与高田田等^[43]对巢湖流域农业面源氮磷负荷输出结果一致；在中低强度的降雨下，排放的农田尾水会留存在沟渠和池塘中^[44]。沟渠和池塘能够有效阻断农田氮磷输出过程，增加水力停留时间，且沟渠坑塘内大量水生植物能够充分吸收氮磷^[45-46]。蒋巷联圩土地利用以农田为主，区内河网密布，沟渠坑塘众多，研究发现暴雨7 d内的排放水量占总降水量的61.99%，因此，建议将电排灌站作为鄱阳湖流域围垦区内农业面源污染控制的重要节点，统筹电排灌站的运行管理，结合天气预报，充分利用围垦区内沟渠池塘缓滞作用，减少农业面源污染物排放，降低对受纳水体水环境质量的影响。

表 1  蒋巷联圩电排灌站排水的主要污染物负荷估算表

Table 1.  Estimated loading of drainage main pollutants from electric drainage and irrigation stations in Jiangxiang Lianwei

水期	位置	排水去向(赣江)	排放负荷/(t·a−1)				负荷占比/%
			TN	TP	COD		TN	TP	COD
涨水期	上段	中支	2.39	0.19	9.09		2.15	3.05	1.81
	上段	南支	9.58	0.56	30.63		8.62	9.11	6.10
	中段	中支	33.30	1.53	112.89		29.98	24.84	22.49
	中段	南支	16.16	1.49	106.12		14.55	24.22	21.14
	下段	中支	16.05	0.70	79.11		14.45	11.33	15.76
	下段	南支	25.83	1.08	96.70		23.26	17.64	19.27
洪水期	上段	中支	0.04	0.01	0.15		0.03	0.16	0.03
	上段	南支	0.23	0.03	1.23		0.21	0.45	0.25
	中段	中支	1.41	0.11	10.62		1.27	1.78	2.12
	中段	南支	1.89	0.11	16.98		1.70	1.82	3.38
	下段	中支	0.61	0.07	4.51		0.55	1.13	0.90
	下段	南支	3.59	0.27	33.91		3.23	4.45	6.76
合计			111.07	6.14	501.95		100	100	100

 | Show Table

DownLoad: CSV

3.   结论

1) 2023年排水量(2.69×10⁷ m³)较2021年减少了75.04%，涨水期因暴雨频率增加引起排水量占比增加，由2021年的63.36%增加到2023年的89.08%。年总排水量为中段>下段>上段；赣江南支排水方向的排水量大于赣江中支。

2) 2023年电排灌站排水中TN、TP的质量浓度年均值超过地表III水，其中涨水期TN和TP质量浓度最高，平均值分别为(4.03±1.61) mg·L⁻¹和(0.23±0.15) mg·L⁻¹，其平均超标倍数分别为3.03、0.31；蒋巷联圩上段电排灌站排水中TN、TP质量浓度最高。

3) 2023年电排灌站排水中TN、TP、COD负荷分别为111.07、6.14和501.95 t·a⁻¹，涨水期其对应的占比分别为93.01%、90.20%和86.57%。中段TN、TP、COD负荷量占比最高，分别为47.50%、52.67%和49.13%；赣江南支排水方向排水TN、TP、COD负荷大于赣江中支。