基于SWAT滇池流域农业非点源污染空间分布特征研究.docVIP

基于SWAT滇池流域农业非点源污染空间分布特征研究 　　摘 要：该研究应用SWAT模型对昆明市滇池流域进行了流域农业非点源污染负荷的模拟计算，利用昆明水文站干海子2008―2014年的实测日径流数据进行了模型的率定与验证，验证结果显示，SWAT模型适用于滇池流域。并且利用验证后的模型分析了2008―2014年滇池河流域年均的非点源污染空间分布特征，结果表明，有机氮和有机磷负荷输出的地区大体一致。总体上，滇池河流域农业非点源污染的较大负荷产生量集中在中部地区，少量分布在南北部角落其他地区输出负荷相对较少。 　　关键词：农业非点源污染；滇池；空间分布；SWAT模型 　　中图分类号 X32 文献标识码 A 文章编号 1007-7731（2016）24-0080-04 　　水污染从源头上划分可分为点源和非点源污染2种类型，与点源污染污染相比，非点源污染由于涉及范围广、控制难度大，目前已成为影响水体环境质量的重要污染源[1]。国外对农业非点源污染的研究和认识始于20世纪60年代，主要是美国、英国、日本等发达国家开展了对农业非点源污染的分析和研究。我国的农业非点源污染研究始于20世纪80年代初，国内众多流域已有SWAT模型研究的案例，如长江流域[2]、黄河流域[3]以及黑河流域[4]。国内学者郝芳华[5]以官厅水库作为研究区，对水文年的非点源污染负荷进行了模拟，结果表明流域非点源污染负荷与降水量成正相关。万超[6]通过对3个不同的水平年进行非点源污染负荷计算，得到了潘家口水库流域的负荷时空特征及主要影响因素。 　　1 SWAT模型概述 　　SWAT模型是一种具有较强物理机制的长期分布式水文模型，可应用于大型流域的研究。SWAT模型集成了遥感（RS）和DEM（数字高程模型），可以模拟和预测长期不同管理对流域的非点源污染负荷输出影响，是农业非点源污染建模的有效工具，目前已广泛应用于世界各地。 　　2 研究区概况 　　2.1 自然地理状况 滇池流域地处云贵高原中部，隶属金沙江水系。滇池流域面积2 920km2 为南北长东西窄的湖盆地，地形主要以山地丘陵居多，面积2 030km2，约占69.5%；湖滨平原面积590km2，占20.2%；滇池流域水域面积约300km2，占10.3%。滇池水域分为草海、外海2个部分。草海外海分别位于滇池北部和南部外海为滇池主体，草海、外海均有人工控制出口，分别为西北部的西园隧洞和西南端的海口中滩闸，滇池出水经螳螂川、普渡河流入金沙江。 　　2.2 流域农业非点源状况及原因 滇池流域的非点源污染方式主要包括生活污水排放、水土流失、地表径流氮磷流失、农田农药和化肥的不合理施用等。研究表明，水体污染和富营养化的主要来源为农药和化肥的使用，美国每年生产5亿t的农药，其中70%用于农田，由此引发的水污染事件不计其数，USEPA把农业列为全美河流和湖泊污染的第一污染源[8]。据统计，滇池流域的水田面积最大，占38%，旱地35%，菜地19%，花卉地占8%[9]，随着农业产业结构的调整和经济发展的需求，蔬菜花卉等经济作物在滇池流域的种植规模将越来越大，作物施肥量也不断递增，化肥的不合理使用导致氮、磷元素大量流失，随着降水和地表径流进入水体，从而导致滇池水体污染和富营养化。 　　3 滇池流域非点源研究的SWAT模型应用 　　3.1 模型基础数据库建立 　　3.1.1 空间数据库 SWAT空间数据库主要包括数字高程模型，土地利用数据和土壤分类数据。数字高程模式DEM和土地利用遥感影像数据来自地理空间数据云，DEM主要用于确定累积流量计算的流向，河网以及特征参数的提取等信息。土地利用遥感数据主要通过ENVI5.1软件预处理，使用监督分类方法提取土地利用分布。土壤分类数据来源世界和谐土壤数据库（HWSD），根据流域边界裁剪后进行重分类。滇池土地利用分类如图1所示。 　　3.1.2 流域属性数据 SWAT模型属性数据库主要包括气象数据、土壤分属性数据，水文数据等，其中气象数据来自由中国水利水电科学研究院、中国气象局及中国科学院寒区寒区环境与工程研究所等多家参与研发的SWAT模型中国大气同化驱动数据集CMADS（自CMADSV1.1版本数据集，水文数据采用昆明市和干海子水文站实测水文数据。 　　3.2 SWAT模型模拟 　　3.2.1 子流域划分 划分子流域时，首先加载流域DEM数据，然后ArcSWAT调用ArcGIS的水文模块进行无洼地DEM生成、汇流累积量计算、水流长度计算河网提取以及流域分割等过程，并最终将整整个流域划分为若干个子流域。在子流域划分过程中，需要设定子流域面积阈值，子流域阈值面积大小决定了子流域划分的个数，子流域面积越小，河网密度越大河网节点越多，不过集水面积阔值并非确定的越小越好，阔值越小更