0 前言

国务院《水污染防治行动计划》《甘肃省水污染防治工作方案(2015—2050年)》中明确提出渭河水系考核断面水质优良比例达到100%的目标。“十三五”以来，为改善渭河水质，甘肃省连续开展了多项渭河流域水污染防治行动，通过现场督查、断面考核、媒体曝光等措施，狠抓污染源头、实施工业企业废水深度治理、城镇污水处理厂提标改造和新扩建工程建设等专项工作，桦林断面水质与往年相比有很大的改善，枯水期水质恶化严重的问题得到有效缓解，水污染防治工作发挥了积极的作用，但长期来看渭河桦林断面依然存在水质波动明显、月度水质不能稳定达标的风险。面源污染作为影响水质安全的重要污染来源，引起了管理部门的广泛关注，其具有的随机性、广泛性、滞后性、模糊性、潜伏性等特点，加大了相应的研究、治理和管理政策制定的难度。开展渭河桦林断面汇水区面源污染负荷的量化分析，掌握桦林断面汇水范围面源污染现状以及点源污染特征，可为断面所在流域污染防治提供有效参考，对控制流域面源污染、改善流域水环境和促进水质稳定达标具有重要意义。

模型评估是全面掌握面源污染情况的有效途径之一，国内外学者采用的面源模型众多^[1-16]，根据模拟参数是否考虑物理机制过程和时空尺度转换，可以归纳为经验模型和机理模型2大类。近年来，我国生态环境部加快推进农业面源污染监测评估技术攻关、试点示范、文件发布和技术帮扶，以改善水环境质量为核心，以“天地协同监测、模型评估核算”为基本思路，探索建立了全国农业面源污染监测评估技术体系。“十四五”期间我国国家层面农业面源污染监测评估任务已出台技术文件，未来将基于我国自主研发的遥感分布式污染估算(DPeRS)模型系统开展监测评估^[17]。DPeRS模型系统以遥感像元为基本模型运算单位，可实现污染负荷空间精细化、可视化表征。与传统总量减排核算方法相比，该模型引入遥感数据，有效降低了对地面数据的要求，并实现了从“点”到“面”的空间突破，可为因地制宜地制定面源污染防控方案提供技术支撑。目前该模型算法已在长三角地区4个省份完成了典型试点示范应用，未来将作为主要技术方案支撑全国层面的农业面源污染监测任务。

现采用DPeRS模型对渭河桦林断面汇水区2018年总氮(TN)、总磷(TP)、氨氮(NH₄⁺-N)和重铬酸盐指数(COD_Cr)面源污染负荷进行空间估算，其中污染类型覆盖农田径流、畜禽养殖、农村生活、城镇生活和水土流失等，基于模型结果和大数据开展流域面源污染源解析及污染成因分析，为该区域面源污染综合防治、科学管控，促进断面持续稳定达标提供技术支撑。

1 研究方法 1.1 研究区概况

渭河发源于甘肃省渭源县鸟鼠山，流经陇东高原、天水盆地、关中平原，至潼关入黄河，是黄河最大支流，地形地貌特征属于黄土丘陵沟壑区，丘陵沟壑纵横，地形起伏较大，植被条件差，降水较少。桦林断面(103°51′—104°55′E，34°13—35°28′N)是渭河上游天水市入境断面，受地形特征影响, 桦林断面汇水区区域气候较为敏感，在我国植被区划中属于温带草原区，在我国生态地理区划中，研究区北部为暖温带半干旱黄土高原中北部草原区，南部为暖温带半湿润汾渭盆地落叶阔叶林、人工植被区。渭河桦林断面所在汇水区平均海拔约2 000 m，地势呈西高东低走势，西部地区生态环境脆弱, 是水土流失较为严重的地区。汇水区总面积约为7 845 km²，地跨甘肃省定西市和天水市，涉及6个县级行政区，其中包括陇西县、漳县全县，以及岷县、渭源县、武山县、临洮县4个不完整县段，其中武山县属天水市辖区，其他5个县域/县段均属定西市辖区，渭河桦林断面所在汇水区地理位置见图 1。

1.2 DPeRS模型方法

采用DPeRS模型对渭河桦林断面所在汇水区面源污染进行分类型分指标估算。DPeRS模型是基于二元结构原理构建的遥感分布模型，该模型综合考虑了自然要素(降水、植被盖度、地形和地貌等)和社会经济要素(施肥利用效率、人口、牲畜和家禽等)对面源污染量的影响。模型算法以遥感数据为驱动数据，耦合定量遥感模型和生态水文过程模型对流域尺度面源污染负荷的时空动态进行定量分析。模型采用交互式数据语言-可视化图像环境(IDL-ENVI)平台进行开发调试。模型空间计算单元为影像栅格，模型参数考虑流域水文特征、土地利用和土壤特征等确定。具体技术路线见图 2。

DPeRS模型具体包括农田氮磷平衡核算、植被覆盖度定量遥感反演、溶解态污染负荷估算、颗粒态污染负荷估算和入河估算5大模块。模拟类型包括农田径流、农村生活、畜禽养殖、城镇生活和水土流失5个类型。污染物模拟形态分为溶解态和颗粒态2种，模拟指标包括TN、TP、NH₄⁺-N和COD_Cr。模型耦合的遥感参数包括植被覆盖度和土地利用，DPeRS模型核心算法详见文献[17-25]。

1.3 数据来源与预处理

模型输入数据主要包括渭河桦林断面所在汇水区2018年的月度降水量、土地利用、月度植被覆盖度、坡度坡长、土壤数据和农田氮磷表观平衡量等，所有数据均统一转换为30 m空间分辨率的栅格数据。月降水量、植被覆盖度、农田氮磷表观平衡量数据空间分布见图 3(a)(b)(c)(d)。数据来源与处理方法见表 1。

2 结果与讨论 2.1 土地利用分布特征

采用2~8 m高分系列卫星遥感影像数据，利用人机交互式解译和面向对象信息提取方法获取2018年桦林断面所在汇水区的土地利用信息。渭河桦林断面所在汇水区8 m分辨率土地利用类型空间分布见图 4。由图 4可见，渭河桦林断面北部上游的临洮县县段、渭源县段、陇西县和武山县段内土地利用类型主要为旱地，林地和草地主要分布在西部和南部山区的漳县和岷县县段。

渭河桦林断面所在汇水区土地利用类型面积及占比统计信息见表 2。由表 2可见，土地利用类型中面积占比较多的依次为耕地、草地、林地、建设用地，其中耕地包括水田和旱地。耕地面积合计为4 936.39 km²，占比62.16%；林地、草地面积合计为2 596.05 km²，占比32.69%；建设用地为370.87 km²，占比4.67%；水域、未利用地面积较少，占比均不足1%。

2.2 面源污染物排放分布特征

渭河桦林断面所在汇水区面源污染排放负荷的空间分布见图 5(a)(b)(c)(d)。各行政区面源污染排放负荷及排放量见表 3。由图 5可见，TN和NH₄⁺-N面源污染排放负荷空间分布规律相似，氮磷型面源污染(TN、NH₄⁺-N和TP)高负荷区主要分布于陇西县、武山县县段和岷县县段；不同县级行政区之间面源TP排放量差异较大，武山县TP排放负荷较为突出；COD_Cr面源污染高负荷区主要分布在陇西县、渭源县县段和武山县县段。

由表 3可见，渭河桦林断面所在汇水区TN、TP、NH₄⁺-N和COD_Cr面源污染排放量分别为11 591，2 697，7 141和1 458 t。

2.3 面源污染入河分布特征

基于研究区2018年各指标面源污染排放负荷，结合各水文分区的入河系数，估算了渭河桦林断面所在汇水区面源污染入河量。汇水区面源污染入河排放负荷空间分布差异明显，见图 6(a)(b)(c)(d)。各行政区面源污染入河量见表 4。由图 6可见，氮磷型面源污染(TN、NH₄⁺-N和TP)高负荷区主要分布在武山县县段、陇西县、临洮县县段的大部分区域，岷县县段的局部区域；COD_Cr面源污染高负荷区除了汇水区内漳县分布较少，其他县域污染负荷区呈分散分布。由表 4可见，2018年甘肃省渭河桦林断面所在汇水区TN、TP、NH₄⁺-N和COD_Cr面源污染入河量分别为2 184，512，1 347和263 t。

2.4 水土流失量空间解析

基于DPeRS模型得到2018年渭河桦林断面所在汇水区水土流失量为74.7万t，产沙量为381.8万t。渭河桦林断面所在汇水区水土流失量及泥沙负荷空间分布见图 7(a)(b)。

由图 7(a)可见，漳县西部地区水土流失较为严重，陇西县东部和武山县地区的水土流失量表现为轻度水平以上的零散分布。由图 7(b)可见，漳县西部、陇西县、渭源县县段北部地区泥沙负荷量较高，由此携带的颗粒态面源污染物可能较多，应该引起足够重视，此外，临洮县县段也存在较高泥沙负荷量，存在水土流失量较高的风险。

2.5 桦林断面水质、水量及超标规律分析

2011—2015年桦林断面月度水质类别见表 5。有13个监测月实测水质为劣Ⅴ类，其中2013—2015年最为严重；“十三五”期间，2019年3月出现1次Ⅳ类水质的情况。相较“十二五”期间水质波动明显、月度水质不能稳定达标的情况，“十三五”以来，桦林断面水质有很大的改善，枯水期(1—3月)水质恶化问题得到有效缓解。

桦林断面2011—2020年的月平均流量监测数据见表 5。由表 5可见，年均流量呈现先升后降再上升的波动态势，其中2012—2013年出现小波峰后，2018—2020年出现最高点，2018年7月丰水期流量达到145 m³/s。数据显示，桦林断面枯水期(1—3月)最差水质出现的年份流量相对较小，而枯水期多出现水质不达标现象，因此可以看出，与我国东南部汛期(6—9月)污染特点相反，桦林断面汇水区较少存在入河污染物旱季(1—3月)藏污纳垢，汛期“零存整取”的情况，枯水期是出现水质不达标的关键时期。

2.6 面源污染类型解析

渭河桦林断面所在汇水区不同类型面源污染入河量统计信息见表 7。

由表 7可见，农田径流是甘肃省渭河桦林断面所在汇水区氮型(TN和NH₄⁺-N)面源污染的首要污染类型，TN的次要影响类型为水土流失，NH₄⁺-N的次要影响类型为农村生活污水；农田径流是磷型面源污染的首要污染来源，其次为水土流失。畜禽养殖是COD_Cr型面源污染的首要污染类型，其次为农村生活类型。

3 结论与建议 3.1 结论

(1) 渭河桦林断面所在汇水区生态本底脆弱，水资源匮乏，河床萎缩，生态水量难以保障，枯水期污水厂排出的废水成为河流主要的水源，也是“十二五”期间河流水质超标的疑似重要原因之一。桦林断面所在汇水区气候以干流为界，南北两岸差异明显，南岸降水较多；北岸及上游区属气候干旱，降雨量较少，土地裸露，植被稀疏区域易导致水土流失严重。汇水区内渭河上游陇西县、渭源县县段人口密度较大，生产、生活用水量较多，同时受引蓄水工程等影响，河川径流量明显减少。“十二五”期间水资源供需矛盾紧张，生态流量不足，有13个监测月实测水质为劣Ⅴ类水质，枯水期水质超标问题凸显，“十三五”期间通过施行一系列水源涵养与保护修复措施，水量有所提高，仅有1个监测月为Ⅳ类水质，其他月份均达到了Ⅱ、Ⅲ类，因此稳定的生态流量是水质达标的重要因素。

(2) 由于生态流量较小，桦林国控断面所在汇水区很少出现汛期污染强度较大的情况，枯水期多因人为排放而出现水质不达标的情况，而汇水区人口、企业点源、工业比重较大的城镇主要集中在干流，因此枯水期污染治理仍是保障水质稳定达标的关键期。

(3) 桦林国控断面所在汇水区面源污染主要来自种植业、畜禽养殖业，汇水区面源污染防控重点应是控源(源头治理)，主抓种植和畜禽养殖等污染源排放，尤其是农田养分流失导致的面源污染排放。

3.2 建议

从甘肃省渭河桦林断面所在汇水区面源污染排放和入河负荷空间分布来看，面源污染防控应着重从农田管理措施方面下手，推行测土配方施肥技术，提高化肥利用率，有效降低土壤养分盈余；采用滴灌、沟灌等节水灌溉措施，逐步推广水肥一体化技术，有效降低土壤养分流失，进而降低农田生产过程导致的面源污染排放。综上，建议从节水节肥等方面提高农业生产效率，保证农田土壤养分供给的前提下减少养分流失，以降低农田径流导致的面源污染排放。