基于gis的松花湖流域非点源污染物空间分布研究.docxVIP

基于gis的松花湖流域非点源污染物空间分布研究 非点源污染是指影响广泛非特定区域的污染物的污染，以及降水和地表径流引起的污染。近年来，由于非点污染的发生机制复杂、影响严重、监督管理和管理难度，非点源污染引起了国内外学者的关注。其中，污染负荷量化是非点源污染研究的重点和难点。 非点源污染物的流失是与水文早期循环阶段有关的现象,受到区域内降雨强度、土地利用类型、地形坡度及植被覆盖等多种因素影响.因此,污染负荷的量化研究首先应模拟区域内的水文循环和土壤流失,应用相关模型估算非点源污染物的侵蚀量.在众多的估算模型中,美国通用土壤流失方程(USLE)是学术界广泛应用的模型之一.USLE与地理信息系统(GIS)相结合,用于估算小流域集水区内土壤流失、泥沙及污染物的输出量,为非点源污染发生区识别、污染负荷量计算提供了快捷、便利、精度较高的方法及手段.然而GIS与USLE结合应用于集水面积200 km2以上流域的实例在国内外尚未见报道.本文拟在这一方面进行尝试.在地理信息系统Arc/info与Citystar软件支持下,结合通用土壤流失方程和污染物负荷模型,估算了吉林松花湖流域的土壤侵蚀模数及氮磷负荷量.旨在研究GIS技术在大流域污染负荷定量计算中的适用性及污染物流失与相关的环境因子之间的关系. 1 研究领域和方法 1.1 入湖流域流域及河流类型 松花湖流域位于吉林省中部,地理位置为东经125°05′~128°48′、北纬41°43′~44°00′,集水区面积43 370.8 km2.流域内地形复杂,地貌类型多样,70%为海拔大于500 m的山地,平川、谷地海拔多在200~400 m. 松花湖流域处于温带大陆性季风气候区内,年均气温1.9℃~4.4℃,年均降水量650~850 mm,6,7,8三个月降水量占全年降水量的90%左右. 汇入松花湖的较大入湖河流有蛟河-漂河、辉发河、头道江、二道江等.根据几大入湖河流的分水岭位置,将流域分为蛟河-漂河区、吉林松花湖区、辉发河区和松花江源头区4个次级流域,见图1. 1.2 土壤条件的拟合 借助GIS技术对空间地理信息的快速采集、分析,以及表达信息间相互作用关系的功能,将松花湖流域的地形、土壤、土地利用类型等各种地理因素通过资料收集、整理后进行信息处理,建立空间与属性数据的统一框架,实现流域内复杂地理环境的拟合,结合模型估算不同地理条件下土壤的流失量和氮磷元素的负荷量. 1.2.1 模型选择 (1) 土壤流失量更新模型 A=0.224·R·K·L·S·C·P. 本文在预测松花湖流域土壤流失量时,选用了修改模型. A=R·St·Ls·C. 式中:A为土壤流失模数(t/(km2·a));R为降雨侵蚀力因子;St为土壤侵蚀因子;Ls为坡度坡长因子;C为植被类型因子. (2) 降雨径流风速 Lj=α∑j=1mCnj?Xn?Sn.Lj=α∑j=1mCnj?Xn?Sn. 式中:Lj为降雨径流中j类污染物负荷量(kg/a);Cnj为n类土地利用类型降雨径流中j类污染物浓度(mg/L);Xn为n类土地利用类型中土壤的流失模数(t/(km2·a));Sn为n类土地利用类型的面积(m2);n为流域土地利用类型编号,1~10分别代表17.5°旱地、17.5°旱地、林地、草地、水田等;α为换算系数(10-3). 1.2.2 模型参数的选择 (1) 年平均降雨p 采用Wischmeier经验公式计算: R=∑{1.735 exp[1.5 log10(Pi2/P)-0.818 8]}. 式中:P为流域内年平均降雨量(mm);Pi为流域内各月平均降雨量(mm). 根据松花湖流域各分区雨量站1989~1998年降雨资料,整理得出各月各站平均降雨量及年平均降水量,代入公式计算出各分区的R值,结果见表1. (2) 松花湖流域土类及土类转换 St与土壤的机械组成(砂、粉砂、粘粒的含量)、有机质含量等密切相关.根据松花湖流域土壤类型的性质,按国际土壤分类制转换图的土类转换,并据土壤侵蚀性因子的一般指标,获得相应的St值(见表2). (3) 坡度m的确定 松花湖流域地貌类型复杂,Ls的确定是根据流域的地形地貌特征,综合其他地区进行坡地土壤流失试验的经验形成的非线性方程: Ls=(λ22.13)m(0.085+0.045S+0.00025S2).Ls=(λ22.13)m(0.085+0.045S+0.00025S2). 式中:λ为平均坡长;S为平均坡度;m为系数,m值的确定与坡度相关. 坡度≥22.5°时,m=0.30;坡度为17.5°~22.5°时,m=0.25;坡度为12.5°~17.5°时,m=0.20;坡度为7.5°~12.5°时,m=0.15;坡度7.5°时,m=0.10. (4) 不同覆盖度的参数设定 根据松花湖流域的具体情况,将植