宁夏引黄灌区典型排水沟道农田退水污染负荷特征分析

陆 阳，仝炳伟，王 乐，韩小龙

（宁夏水利科学研究院，宁夏银川 750021）

当工业点源污染问题得以解决后，面源污染成为导致目前全球水质恶化的主要原因，且其对水体污染的贡献率逐年增加[1]。研究表明，面源污染中以农业面源污染对水体污染贡献率最大，全球有30%～50%的地表水体受到农业面源污染的影响[2-3]。2021年，生态环境部、农业农村部联合印发《农业面源污染治理与监督指导实施方案（试行）》通知，明确指出，农业面源污染治理是生态环境保护的重要内容，应根据农业污染源类型分布、地理气候条件、环境质量状况等，确定农业面源污染优先治理区域，推进重点区域农业面源污染防治；在种植业面源污染突出区域，实施化肥农药减量增效行动，优化生产布局，推进“源头减量-循环利用-过程拦截-末端治理”工程等。

宁夏银北灌区地处青铜峡河西灌区下游，面积3 471 km2。灌区现有西干渠、唐徕渠、汉延渠、惠农渠四大引水干渠，总长726.9 km。灌区排水沟系统自20 世纪50 年代开始修建，形成了目前的主要干沟，并逐年将各级支、斗、农、毛级沟道进行合理配置。主要排水干沟有第二排水沟、银新沟、第三排水沟、第四排水沟、第五排水沟，总长277.6 km。研究排域农业面源污染，对规范肥料、灌溉水的施用量和方式，提高氮磷的利用效率，降低农业氮磷对水环境的荷载，更好地推动“两减一增”行动，促进农业高质量绿色发展具有十分重要的意义[5]。本文以引黄自流灌区第五排水沟排域为典型区域，开展农田退水对干支沟排水氮磷污染物负荷变化特征规律研究。在第五排水沟的干沟和主要支干沟布设监测断面[4]，监测排水沟农田退水水量及水质中总氮（TN）、铵态氮、硝态氮、总磷（TP）4 项农业面源污染特征指标，旨在评估排域农业面源污染程度，揭示第五排水沟农田区域退水氮磷输移量。

1 材料与方法

1.1 第五排水沟基本情况

第五排水沟于1958 年建成，西起贺兰县东南，东至石嘴山市园艺镇入黄河，流经贺兰县、平罗县、惠农区，与第四排水沟相通，全长87.2 km，排水面积8.33 万hm2，现状排水能力56.5 m3·s-1，年排水量1.02 亿m3。第五排水沟排域主要种植作物有水稻、玉米、小麦、露地蔬菜、油葵、枸杞等，地下水埋深1.1～4.0 m，主要作物的种植结构、灌溉及施肥方式在引黄自流灌区具有代表性，区域沟道排水以农田退水为主，农田退水按农-斗-支-干沟进行排放最终汇入黄河。

第五排水沟主要支干沟有8 条，分别为：周城支沟（全长8.15 km），新五一支沟（13.17 km），老五一支沟（3.60 km），五二支沟（17.95 km），第六排水沟（五三支沟）（28.7 km），五四支沟（18.75 km），五五支沟（12.8 km），黑龙沟（14.94 km）。

1.2 监测断面布设

在第五排水沟首部、主要支干沟汇入点后、县界、入黄口等关键节点处，沿排水沟水流方向分别设置监测断面[6]。监测于2021 年进行。自上游到下游依次选取的监测断面为姚伏沙渠、新桥四队、李岗、庙台和熊家庄，支干沟选取的监测断面为周城、新五一、老五一、五二、五三、五四、五五支沟。监测断面布设如图1所示。

图1 第五排水沟监测断面布设示意图

1.3 监测指标及方法

主要监测指标有：排水干沟、支干沟监测断面的沟道水量、水质，排水沟农田退水的污染负荷。

1.3.1 排水沟沟道水量

采用标准断面量水法监测沟道排水量，用ADCP或流速仪率定水位流量关系曲线，用电子水尺和雷达水位计测量水位，经过计算分析得出沟道实时排水流量和不同时段累计水量[7]。

1.3.2 排水沟沟道水质

水样采集方法：在排水沟各断面水面以下1/3 处取样，采集水样位点距离排水沟岸边2/4 和3/4 处，分别各取一个水样，每个水样1 L。

水样采集频次：水样采集时间为灌区头水至冬灌结束；主要作物灌溉期每灌溉1次采集水样1次；作物全生育期采集水样4～5 次；灌溉高峰期加测1～2 次，春灌及冬灌加测1次。

水质指标检测方法：TN采用碱性过硫酸钾消解紫外分光光度法，TP 采用钼酸铵分光光度法采用靛酚蓝法采用酚二磺酸分光光度法[8]。

1.3.3 排水沟农田退水污染负荷的水文估算方法

依据水文学原理，本文利用有限的水量水质观测资料对排水沟的非点源污染负荷进行估算[9-10]。基于排水沟“年内点源排放相对稳定，农田面源集中发生于灌溉期”的特点，点源污染负荷通过排水沟“基流”推求，非点源污染负荷利用排水沟总负荷减去点源污染负荷间接得出。排水沟汇入接纳的主要为农田退水，年内接纳点源污染较少，排水沟农田退水污染物负荷采用下述公式计算：

（1）式中：L为污染物负荷，kg；φ为污染物浓度，mg·L-1；q为流量，m3·s-1；t为过水时间，s。

2 结果与分析

2.1 退水水量和水质变化特征

2.1.1 排水量年内动态

从5个监测断面流量动态（见图2）可知，姚伏沙渠监测断面处于排水沟最上游，农田较少，主要为湖泊和鱼池，该监测断面流量变化比较平稳，周期性变化不明显。姚伏沙渠夏灌期，6—7 月初，监测断面流量较大，流量在0.509～0.549 m3·s-1；冬灌期，监测断面流量较小，流量在0.033～0.05 m3·s-1。下游4 个监测断面的逐日平均流量动态随排域内灌水呈现较明显的周期性上下波动，6—8 月，各监测断面流量均较大，新桥四队以下断面流量为李岗＞庙台＞熊家庄，主要原因是姚伏沙渠到庙台监测断面之间存在多处取水泵站。比较各监测断面流量，李岗变化最明显，夏秋灌期，6—9 月，监测断面流量较大，在3.42～10.60 m3·s-1；冬灌期，监测断面流量较小，在0.56～1.86 m3·s-1。下游入黄河监测断面熊家庄，断面最大过流量时段为5 月中下旬，最大流量为10.60 m3·s-1，最小过流量时段为10月中下旬至11月上旬，平均流量为0.34 m3·s-1。

图2 第五排水沟5个监测断面流量动态

2.1.2 退水污染物日负荷变化特征

图3 第五排水沟5个监测断面污染物日负荷变化过程

上游姚伏沙渠断面排水污染TP日负荷变化过程与其他断面有明显差别，可能由于养殖废水和生活污水对该断面的影响更大。姚伏沙渠断面排水污染负荷变化特征表现为：5—8月氮素以为主，TN 其次，6 月22 日左右日负荷（185.738 kg·d-1）和TN日负荷（158.877 kg·d-1）同时达到峰值日负荷达到峰值持续十几天后迅速下降，在8月1日和9月22 日左右又各出现一次较小的峰值，而TN 日负荷同时达到峰值即刻迅速下降；在6—8 月日负荷高于TN 日负荷，这是由于上游湖泊来水和农田施肥，排水量相对很小，污染物浓度较高，污染负荷相对较大[11]。

2.1.3 退水污染物负荷年内分布规律及负荷估算

表1 熊家庄监测断面排水污染月负荷估算值单位：kg

5—11 月第五排水沟过熊家庄断面污染物量为377 984.458 kg，其中TN 为185 524.720 kg，占污染排放总量的49.08%；，占污染排放总量的33.33%；，占污染排放总量的10.90%；TP 为25 249.499 kg，占污染排放总量的6.68%。可以看出，TN 是排水沟排放量最大的污染物，其次是，第三是TP 最小。

2.2 第五排水沟支干沟退水水量和水质变化特征

2.2.1 支干沟退水流量及污染物浓度

周城支沟流量和污染物浓度变化过程如图4a 所示，在作物夏秋灌期（5—8 月）支沟有径流，其余时段无径流、TN 浓度在6 月22 日左右出现峰值，分别为6.205、9.605 mg·L-1，这与农田施肥关系密切。4—8 月、TN 浓度的平均值分别为2.863、3.256 mg·L-1。根据《地表水环境质量标准》[14]，TN 浓度值超过2 mg·L-1属超Ⅴ类水，表明该排域农田施肥量过多，肥料流失严重；4—8 月，、TP 浓度均值分别为0.715、0.155 mg·L-1，在Ⅲ类水范围。

新五一支沟、老五一支沟流量和污染物浓度动态如图4b、4c所示，二者比较相近，由于它们控制排域面积大，在作物夏秋灌期及冬灌期支沟均有径流。TN浓度在6月12日和7月2日左右出现较大峰值，表明该时段前排域内农田施肥频繁且施肥量超标；4—8 月，新五一支沟、老五一支沟TN 浓度平均值分别为2.611、3.079 mg·L-1，均属于超Ⅴ类水；两支沟浓度在7 月2 日和8 月11 日左右出现峰值，第一次峰值最高，两支沟分别为3.149 mg·L-1（7 月2 日）、5.105 mg·L-1（6月12日）；两支沟和TP浓度均在6月下旬出现 峰值，4—8 月、TP 浓度平均值分别为0.546～0.556、0.164～0.196 mg·L-1，属Ⅲ-Ⅳ类水。冬灌期，农田无施肥，TN 浓度均值1.3 mg·L-1左右，属Ⅲ-Ⅳ类水、TP 浓度均值分别在0.011、0.280 mg·L-1左右，属Ⅱ-Ⅲ类水。

五二支沟流量和污染物浓度动态如图4d所示，作物夏秋灌期支沟有径流，其余时段无径流TN 和浓度在7月2日左右出现相对较大峰值，较上游三个支沟峰值出现时间晚，可能与该排域内作物灌水时间滞后有关；污染物浓度平均值、TN、、TP分别为3.045、4.756、0.639、0.166 mg·L-1，属Ⅳ-Ⅴ类水。

五四支沟、五五支沟流量和污染物浓度动态如图4e、4f所示，作物夏秋灌期支沟有径流，其余时段无径流浓度平均值2.5 mg·L-1左右，TN 浓度平均值3.5 mg·L-1左右，属超Ⅴ类水、TP浓度平均值分别为0.64、0.23 mg·L-1左右，属Ⅱ-Ⅲ类水。

图4 第五排水沟支干沟退水流量及污染物浓度动态

五三支沟（第六排水沟）流量和污染物浓度动态如图4g 所示，5—6 月初，沟水流量逐渐增大，6 月12 日左右受上游沟道取水影响流量减小，而后上升到最大值并一直持续到8 月初，夏秋灌期结束流量逐渐减小到零，此后随着冬灌期来临，沟水流量又逐渐增大。TN 浓度在7 月2 日和冬灌期10 月22 日左右各出现一次峰值，最大峰值7.566 mg·L-1出现在7 月2 日左右；5—8 月期间，TN 浓度平均值2.754 mg·L-1左右，属超Ⅴ类水浓度虽然也在7 月2 日和冬灌期10 月1 日左右各出现一次峰值，但浓度较小，分别为1.594、1.250 mg·L-1，远低于其他支沟；5—8月时段的浓度平均值为0.397、0.404 mg·L-1，属于Ⅱ或Ⅴ类水。冬灌期，排水中、TN 和TP的浓度均显著升高，可能由于田间灌水大量流入排水沟，冲击沉积在沟底的废水引起TN 和TP 浓度的升高。

2.2.2 支干沟污染物日负荷

计算支干沟各监测断面主要污染物TP 和TN 日负荷量，绘制日负荷动态图。由图5 可以看出，各支沟排水污染负荷动态差别明显，灌溉和排水引起了TN和TP 较大的差异。在灌水期，排水沟TN 主要来自农田，TN 虽然浓度相对较低，但由于排水量较大，TN的污染负荷相对较大，随着排域灌水停止，排水量减少，TN 浓度升高，TN 的污染负荷却相对较小。TP 主要来源于养殖废水和生活污水，周城支沟、五二支沟和五三支沟排放比较稳定，TP负荷变化不大；五一支沟和五四支沟，在8 月中下旬TP 升高，出现了峰值，主要是退水冲洗沟道使沉积在沟底淤泥中的磷释放出来迁移到出水口的缘故[15]。

图5 第五排水沟支干沟监测断面TN、TP负荷动态

2.2.3 支干沟污染物负荷年内分布规律及负荷估算

表2 第五排水沟各支沟监测断面的各类面源污染负荷估算值 单位：kg

3 结论与讨论

1）第五排水沟干沟的5 个监测断面流量变化特征：呈现周期性上下波动趋势，年内峰值出现在6—8月（夏灌期），冬灌历时较短且以漫灌为主，所以退水过程峰形较低。上下游监测断面流量呈现周期性的上下波动但差异性很大，上游姚伏沙渠监测断面以上农田较少，流量变化比较平稳，周期性的变化不明显；下游4 个监测断面的平均流量呈现周期性的上下波动，尤其以李岗断面波动幅度最为显著，下游熊家庄入黄河处流量低于新桥四队、李岗和庙台监测断面流量，主要原因为姚伏沙渠到庙台监测断面之间存在多处取水泵站。入黄河监测断面熊家庄处，监测断面最大过流量时段为5月中下旬，最大流量为10.60 m3·s-1，最小过流量时段为10 月中下旬至11 月上旬，平均流量为0.34 m3·s-1。

4）第五排水沟支干沟流量变化特征：各支干排水沟汇入接纳的水源主要来自农田退水，部分支干沟接纳了少量的养殖废水和生活污水。各支沟的流量呈现周期性的上下波动，这种变化规律与各支沟排域的灌溉制度密切相关。

5）第五排水沟支干沟污染负荷年内分布规律和水质变化特征：支沟退水污染负荷主要发生在5—8月农田灌水期，说明支沟退水污染主要来源为农田退水；TN是排水沟排放量最大的污染物，其次是第三是，TP最小。