1. 引言
三峡库区是我国重要的淡水水源地，是我国南水北调中线工程的重要补给水源。由于库区次级河流氮磷浓度较高，自三峡水库蓄水以来，受水位壅升和回水顶托的影响，“水华”频发。随着库区点源污染逐步控制，非点源污染(又称面源污染)逐渐成为库区水环境污染的主要原因，因此开展对三峡库区非点源氮磷负荷的估算和预测研究，对库区的水环境安全和水污染的治理具有重要的理论和实践意义。
随着国家重大水专项的进展，水污染治理进入新阶段，开展对非点源污染的控制和治理成为重要的课题，而输出系数模型在非点源污染年负荷模拟和分析中发挥着重要的作用。目前，很多研究者利用输出系数模型，对国内外相关流域的非点源污染负荷进行了模拟分析研究，结果均表明输出系数模型在资料缺测地区进行非点源污染负荷模拟研究优势明显 [1] - [12]。同时，总体上可以看出，目前的研究多基于输出系数模型采用定性分析等方法对面源污染负荷进行估算，模型的相关计算参数主要通过理论计算或者查阅相关文献来确定，未在研究区开展相关观测实验。
本文拟针对三峡库区的香溪河典型小流域，结合研究区径流小区实验结果及相关文献，基于输出系数模型估算研究区非点源污染负荷，主要包括总氮(Total nitrogen, TN)负荷、总磷(Total phosphorus, TP)负荷，分析不同污染源对研究区非点源污染负荷的贡献，以期为香溪河流域乃至整个三峡库区非点源污染防控提供决策依据。
2. 研究区概况
香溪河是三峡库首的第一大支流，位于湖北省兴山县境内，流域总面积约为3150 km2，拥有南阳河、古夫河、高岚河3条主要支流(图1)。香溪河流域属亚热带大陆型季风气候，多年平均气温17.1℃，多年平均降雨量为900~1200 mm，年均降雨日为130 d，降雨主要集中在5~9月。土壤类型繁多，有黄壤，黄棕壤，石灰土，紫色土，水稻土，潮土等。香溪河流是三峡库区典型的农林复合小流域，土地利用类型主要分为水田、旱地、林地、果园、草地、水域和城镇用地。土地利用类型及面积见表1。
Table 1. The land use of the Xiangxi River basin
表1. 香溪河流域土地利用情况表
随着三峡水库的建设，库区内的人口进行了大规模的迁移，区域内人口快速的增长，库区的畜牧养殖业主要以养殖牛、猪、羊和家禽为主，根据查阅《兴山县统计年鉴》及香溪河流域相关区县2010年政府工作报告等相关资料，得出2010年香溪河流域农业人口为13.1251万人，香溪河流域畜牧养殖规模为：牛2624头，猪312,519头，羊107,249头，家禽378,200只。
因库区人多地少，实践“高投入，高产出”传统的经营模式，在降雨条件下水土流失、化肥流失及养殖污水排放引起的非点源污染是香溪河水体富营养化的重要原因。
Figure 1. The location of the Xiangxi River basin
图1. 香溪河流域意图
3. 模型模拟计算
经过对香溪河流域面源污染源进行分类调查统计分析，明确了库区氮、磷面源污染负荷主要来源分别为土地利用、畜牧养殖和居民生活污染源。
本文在三峡库区香溪河流域及其他与香溪河流域相似的典型小流域不同土地利用区选取典型观测试验区，对不同土地利用单位面积污染负荷输出进行试验观测。通过实验，确定不同土地利用类型的氮、磷负荷输出系数，然后根据整个小流域面积及输出系数，利用输出系数模型，估算得到整个小流域不同土地利用类型的面源污染负荷输出量。
畜牧养殖及生活污染源的输出系数则通过对小流域的人口、社会经济、养殖规模等指标进行调查，结合相关计算标准，并参考目前已有相关文献，计算得出畜牧养殖及居民生活污染源的氮、磷负荷输出系数，根据畜牧养殖规模及居民人口等，估算出小流域畜牧养殖及居民生活面源污染负荷输出量。
3.1. 输出系统模型
Johns [13] 于1996年建立了输出系数模型，该模型引入了牲畜、人口等因素，建立在土地利用与受纳水体的非点源污染负荷量的基础上。模型大大简化了非点源污染形成过程，提高了土地利用状态发生变化时的灵敏性，降低了对监测数据的依赖性和监测费用，为大流域非点源污染的量化模拟研究提供了一种具有一定精度的模型方法。
输出系数模型计算方程式如下，
式中：L为营养物的输出量；Ei为第i种营养源的输出系数(即单位面积年负荷量)；Ai为第i类土地利用类型的面积或第i种牲畜的数量、人口数量；Ii为第i种营养物输入量。
3.2. 确定输出系数
3.2.1. 土地利用输出系数
在明确流域土地利用类型的基础上，获取输出系数的常用途径是查阅文献和现场监测法。现场监测法需要通过长时间监测流域水质水量，通过计算负荷量，得到相应的输出系数值，该方法所获取的参数有较高的精度，但是耗时长、费用高。另一种途径是利用前人的研究成果，通过查阅文献获取系数值。
通过收集三峡库区香溪河流域兴山监测站12个不同土地利用径流小区共6次产流过程，各不同土地利用径流小区全年观测结果见表2。
Table 2. Runoff and TN and TP load of observation station under different land cover
表2. 不同植被覆盖类型径流小区的径流、TN和TP负荷
由表中数据，可将12个不同植被覆盖类型的径流小区分为3类：旱地(对照、秋葵瓜)、林地(杨树、柏树)、草地(荒草)。这三类土地利用类型的输出系数可通过实验数据计算得出。研究区其他土地利用类型如城镇用地、水田、水域等的输出系数，则查阅相关文献获取 [14] [15] [16]。由此得出香溪河流域不同土地利用类型输出系数见表3。
Table 3. The export coefficients values of different land use in the Xiangxi River basin
表3. 不同土地利用养分输出系数
3.2.2. 畜牧养殖输出系数
在三峡库区香溪河流域及其他与香溪河流域相似的典型小流域开展了对农户畜牧养殖调查，结果见表4。
Table 4. The investigation results of livestock and poultry manure disposal in typical watershed
表4. 香溪河典型小流域畜禽养殖调查
结合小流域典型农户调查结果，参照相关文献资料 [17] [18] [19] [20]，确定香溪河流域畜牧养殖污染输出系数见表5。
Table 5. The export coefficients of livestock and poultry in the Xiangxi River basin
表5. 香溪河流域畜禽养殖氮磷输出系数表
3.2.3. 生活污染源输出系数
在调研过程中发现，居民对生活排泄物、垃圾方式处理不一：在部分居民户中，使用的是浅坑旱厕，建于屋外，浅坑旱厕的使用，极不利于农村环境的改善。部分居民户使用水冲厕，这部分居民普遍认为水冲厕使用方便、干净卫生。此外，在调研过程中，课题组发现，约90%的居民户将粪便作为肥料还田，约有1%的左右的居民，直接将粪便堆放在路边或倒入沟渠，极少数统一收集并集中处理使用。
关于人粪尿的排污系数，目前已有相关研究。马啸 [17] 对三峡库区的农村生活调查表明，每个成年人均年产氮磷量为5.715 kg、1.17 kg。农村人的粪尿，部分当作肥料使用，其余直接排放到环境或水体中，流失率约为21%。而农村生活污水基本上全部直接排放到自然环境中。库区内农村地区生活垃圾和水处理设施建设较为滞后，农村生活污水和粪便很难得到有效的处理。结合库区实际情况，确定农村生活污染氮磷输出系数为：农村生活氮输出系数为：1.872 [kg/(ca∙a)]，农村生活磷输出系数为0.214 [kg/(ca∙a)]。
3.3. 计算结果
根据本文通过实验及查阅文献所得的各种污染源的输出系数，根据研究区基础数据资料，基于输出系数模型计算可得不同污染源的污染负荷结果见表6。
Table 6. TN and TP load of different pollution sources in the Xiangxi River basin
表6. 香溪河流域不同污染源非点源TN、TP负荷
4. 结果分析
4.1. 不同土地利用总氮、总磷负荷特征分析
根据香溪河流域土地利用类型情况及本文所得的土地利用输出系数，所得的不同土地利用类型非点源TN、TP负荷估算值见表7。不同土地利用利用总氮、总磷贡献率见图2。
Table 7. Estimated TN and TP load of land use source in the Xiangxi River basin
表7. 香溪河流域土地利用污染源总氮、总磷负荷估算表
从图中可以看出，从总量上看，林地的总氮、总磷排放量最大，其次是旱地、草地和水田。分析其中原因可知，因林地面积占整个流域面积的比例高达88.08%，所以林地的总氮排放量占流域总排放量的63%，总磷排放量占流域总排放量的61%。旱地的总氮、总磷排放量在整个流域所占比例也较大，分别为16%和24%。因此，在香溪河小流域，林地、旱地氮、林流失的控制应为基于土地利用类型的面源污染防控的重点。
Figure 2. Contribution figures of TN and TP of different land use in the Xiangxi River basin
图2. 香溪河流域不同土地利用总氮、总磷负荷贡献率示意图
4.2. 不同污染源总氮、总磷负荷特征分析
香溪河流域不同污染源总氮、总磷负荷示意图见图3。
从图中可以看出，香溪河流域总氮负荷量由高到低依次为畜牧养殖污染 > 土地利用污染 > 居民生活污染，而总磷负荷量由高到低依次为土地利用污染 > 畜牧养殖污染 > 居民生活污染。因此，控制总氮污染负荷的关键在于控制畜牧养殖污水排放，而控制总磷污染负荷的关键在于合理优化农业种植模式，减少农田施肥量。
Figure 3. TN and TP load of different pollution sources in the Xiangxi River basin
图3. 香溪河流域不同污染源总氮、总磷负荷
5. 结论
本文基于输出系数模型，利用试验观测数据计算和考参考文献查阅的方法，确定了土地利用、畜牧养殖及农村居民生活污染源的输出系数，估算了三峡库区香溪河小流域2010年各污染源的总氮、总磷负荷量，并进行了相关分析，结果表明不同土地利用类型非点源氮、磷排放负荷总量最大的为林地和旱地，合理管理林地及旱地化肥施用是防控土地利用污染的关键。对于不同污染源的总氮、总磷负荷排放量，畜牧养殖污染源排放的总氮负荷最高，而土地利用排放的总磷负荷最高。因此，选择合理的农田管理措施，有效控制林地、旱地施肥量及畜牧养殖污水排放是防控三峡库区香溪河小流域非点源污染的关键。
致谢
感谢国家水体污染控制与治理科技重大专项课题对本文的支持。
基金项目
国家重点研发计划(2017YFC0404504)；中央级公益性科研院所基本科研业务费(CKSF2019410TB)。
NOTES
*通讯作者。
作者简介：杜学才，男，1974年8月出生，云南省昆明市人，本科学历，高级工程师，主要从事工程建设与管理、水环境研究等方面的工作。