流域非点源污染负荷估算模型的研究
白龟山水库非点源污染负荷计算

12 .9脒 (s P I )C
式中: 为单位面积土壤侵蚀量( h) t a; / 为降雨侵蚀参数; 为土壤侵蚀性参数; 为地形参数; 为植被覆盖因 S C
子; P为管 理参 数。
2 2 参数确定 . 2 2 1 参数 的确定 ..
与降雨量、 降雨强度、 历时、 雨滴大小和下降速度有关。 它反映了降雨对土壤的潜在侵蚀能力。
全流域的总潜在侵蚀量为 :
‘
=
l
式中: 为流域潜在侵蚀量; 玉为 f 小区单位面积潜在侵蚀量;f i A 为 小区面积。 计算得白龟山水库流域潜在侵蚀量为 158 ×l a 3. 6 ot 。 4/
24 溶解态氮磷输 入量 .
陆地输入水库 中的氮磷污染物分为两大类来计算… , 一类是 由泥沙吸附 的吸附态 污染物 ( 固态 污染物 ), 它随泥 沙
lE g i= 15培 (i/ )一0039 . P P 2 .2
式 中: l P 为月 降雨 量( m ; m )P为年 降雨量 (m 。 m )
l 2
全 侵 参 年的 蚀 数为: E=∑ 墨
i ll
பைடு நூலகம்
全年 的侵蚀参数见表 2 。 表 2 白龟山水库 降雨量分布情 况表
根据表 2 可得 : 年降雨量 为 78a, 3rn年侵蚀参数 1 . 。 o 1 6 9 3
l 通用流失方程的计算
11 流域土壤侵蚀■计算模型与参数 .
1 11 流域分 区 .. .
为计算方便 , 将白龟山水库流域分为 3 个小流域, 每个小流域作为 1 个计算小区, 每个小区的情况见表 1 。
表 1 白龟山水库漉域分 区情况表
2 12 土壤 侵蚀量模型 ..
城市非点源污染负荷估算研究

年 降雨量 的相关关 系, 其相关性显著 ; 在不 同降雨强度条件下 , 对典 型示范 区城市排 污 口实施降雨过程监测 , 获得 了有代表性 的 城 市径流污染物平均浓度 , 而估 算出辽 宁省城 市径流污染 负荷总量。 从
关键 词 : 点 源 污染 ; 表 径 流 ; 流 系数 ; 染 物 平 均 浓 度 ; 染 负荷 非 地 径 污 污
S 模 型 是 由美 国 水 土 保 持 局 ( ol o sra ( S i C nev —
式() , 2 中 S=( 54 0 C 一2 4 C 为 径 流 曲线 2 0 / N) 5 , N 数值 , 一个 无量 纲参 数 。 是
t nSri ) 制 的小 流 域设 计 洪 水 模 型 , 算 过 程 i evc 研 o e 计
非 点 源 污染 是 指 大 气 、 面 或 土 壤 中的污 染 物 地
式 ( ) , 为 地表 径流 深 mm)P 为降雨 量 mm) 1中 D , ,
为初 损量 mm) S为 流域 饱和储 水 量 。 , 由于 不 易 准 确 测 定 , C S S建 立 了经 验 公 式 :
/ a=0 2 . S。则 : D = ( 一0 2 ( +0 8 P . S)/ P . S) () 2
地、 工业用地 、 仓储用地 、 对外交通用地 、 道路广场用地 、 市政公用设施用地 、 绿地和特殊用地等九类 。分别选定
1 S S降 雨 径 流 模 拟 C
美 国水 土保 持 局 提 出 s s模 型 , 计 算 一 次 降 c 是
雨过程 径 流深 度 的经 验 公 式 , 构 简 单 , 算 方 便 。 结 计 其计算 公式 :
简 单 , 需资 料简 单易取 , 反 映不 同土壤 和地 表 覆 所 能 盖 条 件 对产 流影 响 的 特 点 , 一 种较 好 的计 算小 流 是
丹江口水库入库非点源污染负荷的计算与讨论

2 0 1 5年 5 月
民 长 江
Ya n g t z e Ri v e r
Vo 1 . 4 6. No . 1 0
Ma y, 2 01 5
文章编号 : 1 0 0 1— 4 1 7 9 ( 2 0 1 5 ) 1 0— 0 0 4 2— 0 6
摘要 : 正 确 估 算 丹 江 口水 库 入 库 非 点 源 污染 负荷 对 于 水 源 区 水环 境 保 护 具 有 重要 意 义 。根 据 丹 江 口库 区 6条
主要 入 库 河 流 汉 江 、 天河、 堵 河、 丹江、 老灌河、 淇河控制水文站 2 0 1 3年 逐 日流 量 数 据 , 采 用 数 字 滤 波 法 对 基 流
流 中较 稳 定 的 部 分 ( 河川基流 ) 与波动 部分 ( 地 面 径 流) 分 割 开 来 的一 种 水 文 学 方 法 。对 于 污 染 负 荷 而 言, 点 源污 染负 荷和 自然 背景 具有稳 定性 特征 , 年 内变
收 稿 日期 : 2 0 1 5—0 2—1 5
大坝 以上 流域面 积 9 . 5 2万 k m , 多年 平 均入 库 流量 为 1 2 3 0 . 4 m / s , 多年 平 均 年 径 流量 3 8 8亿 m 。丹 江 口 库 区有 湖北 郧 西 、 郧县 、 张湾、 茅箭 、 丹 江 口和 河 南 西 峡、 淅 川 7个 县 ( 市、 区) , 分 布着 1 6条 主要 入库 支 流 , 见表 1 , 其 中汉江 、 天河 、 堵河 、 丹江、 淇 河 和老灌 河等
制方 案 和 主要 入 库河 流 水 环 境 治 理 方 案 具 有 重 要 意
义, 相关 学 者 分 别 采 用 输 出 系 数 法 、 A n n A G N P S模 型 、 S WA T模 型 对丹 江 口库 区非 点 源 污 染 负 荷 进
SWAT模型点源输入识别方法研究--以鞍山市南沙河流域为例

SWAT模型点源输入识别方法研究--以鞍山市南沙河流域为例杨武志;苏保林;罗运祥;张倩【期刊名称】《北京师范大学学报:自然科学版》【年(卷),期】2012(48)5【摘要】点源污染的划分和核定是正确估算流域非点源污染负荷的重要工作.我国城市地区往往缺乏可靠的点源入河数据,因此,在SWAT模型的应用中需要根据有限资料进行点源污染输入识别,以确保模拟结果的正确性.本文以鞍山市南沙河流域为例构建了非点源SWAT模型,探讨了点源输入识别中总量估算、时空分配、营养物质形态划分和迁移过程中降解估算等4个方面的问题和处理方法,评估了单位负荷法、现场试验法和统计数据法等3种点源负荷总量估算方法的优劣.结果表明,统计数据法以实际统计资料为依据,能较好反映点源负荷的年际变化,估算结果较好;点源对总氮和总磷的贡献率分别达到了79%和82.5%,是南沙河流域的主要污染源.因此点源控制是改善南沙河流域水质状况的首选措施.【总页数】5页(P515-519)【关键词】SWAT模型;点源;识别;南沙河【作者】杨武志;苏保林;罗运祥;张倩【作者单位】北京师范大学水科学研究院,水沙科学教育部重点实验室,北京100875【正文语种】中文【中图分类】X522【相关文献】1.鞍山市南沙河流域污染情况浅析 [J], 郭富兴2.基于SWAT模型的南沙河流域日流量过程研究 [J], 杨呜婵;贺国平;刘冀宏3.基于 SWAT模型的小流域面源污染负荷时空分异研究--以重庆市万州区陈家沟小流域为例 [J], 穆婧;史明昌;郭宏忠;魏曦;赵锦序4.基于SWAT模型的北京沙河水库流域非点源污染模拟 [J], 刘博;徐宗学5.流域分布式水文模型SWAT空间输入数据的不确定性研究 [J], 姚允龙;吕宪国;王蕾因版权原因,仅展示原文概要,查看原文内容请购买。
基于SWAT模型的流域非点源污染模拟

基于SWAT模型的流域非点源污染模拟张皓天;张弛;周惠成;沈必成【摘要】利用SWAT(soil and water assessment tool)模型,在GIS技术支持下,通过构建模型所需的空间数据库和属性数据库,以黑龙江蚂蚁河流域为研究区域进行流域非点源污染的模拟研究.分别在时间和空间尺度上对研究区域非点源污染负荷分布进行分析,识别出非点源污染严重的区域及其影响因素.结果表明:在时间尺度上,非点源污染物负荷与降雨量和径流量有较强的相关关系;在空间尺度上,流域西南部地区非点源污染较为严重;不同土地利用类型的非点源污染负荷不同,耕地的非点源污染单位面积负荷最高,疏林地次之,林地等其他土地利用类型的单位面积负荷较小.【期刊名称】《河海大学学报(自然科学版)》【年(卷),期】2010(038)006【总页数】7页(P644-650)【关键词】SWAT;非点源污染;时空分布;土地利用;蚂蚁河流域【作者】张皓天;张弛;周惠成;沈必成【作者单位】大连理工大学建设工程学部,辽宁,大连,116024;大连理工大学建设工程学部,辽宁,大连,116024;大连理工大学建设工程学部,辽宁,大连,116024;黑龙江省水文局,黑龙江,哈尔滨,150001【正文语种】中文【中图分类】X502东北地区是我国重要的农牧业生产基地和最大的老工业基地,城市化水平位于全国前列,近些年来,随着东北地区农业的发展和振兴东北老工业基地战略的逐步深入,东北地区水环境污染问题日益突出.目前,许多发达国家的研究已经证实,非点源污染是导致水环境恶化的主要原因之一[1].为了对非点源污染进行有效的治理和控制,必须定量研究污染物的排放规律.应用GIS技术支持的非点源污染模型,可以对整个流域的非点源污染进行定量描述,分析其产生的时间和空间特征,并评估土地利用的变化对非点源污染负荷的影响,为流域规划和管理提供决策依据.SWAT(soil and water assessment tool)模型是美国农业部农业研究局开发的流域尺度分布式水文模型,该模型已经广泛地应用到流域的水平衡、河流流量预测和非点源污染控制评价等诸多方面.在国外,早期以Arnold为首的工作组分别从美国的国家尺度、流域尺度以及小流域尺度验证了SWAT模型在径流模拟方面的适用性[2-6];Santhi等[7-10]先后将SWAT模型应用于美国得克萨斯州Bosque流域,对模型模拟非点源污染的适用性进行了验证;基于美国经验,SWAT模型在其他国家也得到了广泛的应用[11-15].在国内,SWAT模型的应用主要包括3个方面:产流/产沙模拟、非点源污染研究及输入参数对模拟结果的影响研究.在非点源污染研究方面,胥彦玲等[16]将SWAT模型应用到陕西黑河流域的非点源污染研究中;张永勇等[17]扩展了SWAT模型的水质模块,并将其成功运用到海河流域中;庞靖鹏等[18]以密云水库为例,研究了土地利用变化对产流和产沙的影响.经过国内的研究可以发现SWAT模型是针对北美的土壤植被和流域水文结构来设计的,模型自带的数据库和中国的实际情况有较大出入,而我国东北地区又缺乏比较完善的流域基础数据库,用来描述流域的很多参数如果靠直接测量来获得是很困难的,因此估算模型的参数和验证模型在相关流域的适用性就显得尤为重要.这些问题限制了SWAT等一系列物理型模型在我国东北地区非点源污染研究中的应用,而目前我国东北地区非点源污染负荷量估算仍较粗略[19].针对上述问题,本文选取黑龙江省蚂蚁河流域为研究区域,通过建立和完善研究流域的水文模拟基础数据库,运用SWAT模型来模拟该流域的非点源污染的发生,通过率定SWAT模型的参数和验证模型的适用性,得到了适合研究流域的一套模型参数值;进而计算研究流域内产生的非点源污染负荷量,研究非点源污染负荷的时空分布规律,分析影响其分布的主要因素,并评估不同土地利用类型非点源污染负荷量,为流域的水质管理提供科学的决策依据.1 研究区域概况蚂蚁河流域位于黑龙江省南部,张广才岭西侧,是松花江干流右岸的一级支流,地理位置为东经127°15′~129°33′,北纬43°57′~46°26′.蚂蚁河源海拔高700m,干流河长285km,流域面积10547km2,流域主要为山区和半山区,地势自东向西北,然后向东北倾斜,海拔高程在104~1400m之间,河谷宽阔平坦,地势低平,地面组成物上部为亚黏土,下部为沙砾石[20].流域内的土壤主要有黑土、暗棕壤、白浆土、草甸土、沼泽土等.流域内的主要种植作物为水稻、玉米和大豆,施肥主要以有机肥、尿素、硫酸钾和二铵为主.流域内多年平均气温2.3~3.4℃,年最高气温35.5℃,最低气温-41.6℃,年降雨量在500~900mm之间.研究区域地理位置如图1所示.图1 蚂蚁河流域示意图Fig.1 Mayihe watershed2 基于SWAT模型的非点源污染模拟与验证在收集资料的基础上,构建研究区域水文模拟基础数据库;之后通过划分流域内的空间单元,对非点源污染负荷量进行连续时段的动态模拟;最后利用流域内的实测资料对模型参数进行率定.2.1 流域水文模拟基础数据库的建立在运行SWAT模型前,需要准备必要的地图和数据库以生成SWAT模型输入数据集.模型所需的数据可以分成空间数据和属性数据.其中,SWAT模型所需要的空间数据的准备、修改和存储可以通过ArcGIS 9等GIS软件完成;而降水、温度、土壤等数据则通过多个输入文件以ASCII或者.dbf格式输入.模型根据土地利用和土壤类型表,将土地利用图、土壤类型图与土地利用和土壤数据库进行链接.本研究采用的气象数据是流域内1966—2001年气象站点的实测数据.经过查阅《黑龙江土壤》[21]和中国土壤数据库以及利用美国华盛顿州立大学开发的土壤水特性软件SPAW计算出土壤的属性数据[22].通过统计《黑龙江统计年鉴》和《哈尔滨统计年鉴》中的资料,得到研究区内3个县市的农业管理数据.2.2 子流域划分及水文响应单元(HRU)的生成SWAT模型通过采用数字地形分析技术,利用栅格数字高程图(DEM)来精确描述流域边界、生成流域河网、进行子流域的划分以及生成水文响应单元.流域河网生成的详细程度是由定义的上游集水区面积阈值大小来控制的;子流域的形成由子流域出口位置来控制;SWAT模型采用了不能确定空间位置的水文响应单元的划分方法,使用不同的土地利用类型、土壤植被和坡度组合来生成水文响应单元[23].本研究将流域划分为70个子流域,在整个流域内共生成了509个水文响应单元. 2.3 SWAT模型参数率定与模型验证通过参数敏感性分析,先得到了对径流影响比较大的高敏感参数,然后对这些关键参数进行调整,将1983—1990年定为模型月模拟的率定期,1991—2001年定为模型的验证期.选择相对误差 Re、相关系数R2和Nash-Suttcliffe系数Ens来评价模型的适用性.模型参数的最终取值见表1,率定期和验证期径流率定的结果见表2.模型输出结果与实测平均值之差占实测值的百分比应小于规定标准,即不大于20%;评价系数(Ens和R2)也应达到规定的精度标准,一般要求R2>0.6,且Ens >0.5[24].从表2可以看出研究流域率定期月径流的Re=5.93%<20%,Ens=0.73,R2=0.76,Ens和R2均大于0.7.验证期月径流的Re=2.85%<20%,Ens=0.82,R2=0.82,Ens和R2均大于0.8.说明模拟值和实测值之间的拟合较好,精度可满足模拟要求.将流域内1957—2001年的年径流量观测资料进行频率分析,把频率大于75%的年份定为枯水年,频率小于25%的年份定为丰水年,其余的年份定为平水年.在模型验证期内,有枯水年3年,平水年6年,丰水年2年,基本上体现出了径流量在年际间的变化情况,说明验证期具有代表性.在验证期内得到了较好的月径流模拟值,这说明模型的模拟是符合实际的,经过参数率定的SWAT模型适用于该流域.表1 SWAT模型参数在研究流域的率定值Table 1 Calibrated values of parameters of SWAT model for study areaESCO 土壤蒸发补偿系数 .hru -0.09 CN2 SCS径流曲线系数 .mgt -4 Gw_Revap 地下水再蒸发系数 .gw +0.15 ALPHA_BF 基流a系数 .gw +0.052 NPERCO 硝基氮下渗系数 .bsn +0.3 BC1 氨氮生物氧化速度常数 .swq -0.35 BC2 亚硝氮生物氧化速度常数 .swq +0.5 BC3 有机氮转化为氨氮速度常数 .swq +0.1由于缺乏长时间序列的泥沙和水质监测数据,本研究只对流域的氨氮模拟过程进行了校准.结果表明,氨氮的模拟误差在10%左右,总体反映了实际情况.表2 率定期及验证期月径流率定结果Table 2 Calibrated results of monthly runoff during calibration and test periods率定期 75.52 71.29 5.93 0.76 0.73验证期 57.78 56.18 2.85 0.82 0.823 SWAT模型模拟结果分析选用率定后的模型参数计算流域内产生的非点源污染负荷量.在时间上对非点源污染负荷的变化趋势进行分析,对不同降水频率下的非点源污染负荷量进行模拟和研究;在空间上对非点源污染负荷的分布进行分析,识别出研究区内非点源污染严重的区域;最后评估不同土地利用类型的非点源污染负荷量.3.1.1 非点源污染负荷年际分布通过对研究流域年降雨资料和年径流资料分析得出研究流域的年降雨量和年径流量的相关系数为0.875,有很强的相关性.图2为1993—2001年研究流域内各年径流量实测值与河道总氮、总磷、输沙量变化,由图2可以看出,各年的河道输沙量、河道总氮负荷量和河道总磷负荷量均与径流量的关系密切.这是因为降雨和产汇流是土壤侵蚀和非点源污染物的主要驱动力,当下垫面的条件不变时,土壤侵蚀和非点源污染就受降雨和径流的影响很大.经过计算分析,研究区域内年径流量与相应河道输沙量、总氮负荷量和总磷负荷量三者的相关系数分别为0.973,0.709和0.839.图2 1993—2001年径流量与河道总氮、总磷、输沙量变化Fig.2 Variations of runoff,total nitrogen,total phosphorus and sediment from 1993 to 2001 为了分析降雨量变化对非点源污染的影响,有必要进行不同降雨量下的非点源污染负荷研究.根据流域内1966—2001年的降雨资料,选用皮尔逊Ⅲ型曲线,运用适线法计算流域年降雨量的经验累积频率,得到曲线的算术平均值Ex=591.8mm、变差系数Cv=0.25和偏态系数Cs=0.64.由理论频率曲线可得到在10%,50%,75%,90%降水频率下的年降雨量分别为788.5mm,576.9mm,485.2mm和414.2mm.选取与不同频率下降雨量相近的特征水文年的降雨资料作为模型降雨量输入数据,研究非点源污染负荷与降雨量之间的相关关系.各降水频率下泥沙和氮磷污染负荷估算见表3.由表3中的数据分析,在不同降水频率下,各种非点源污染负荷量随降雨量的增加有增大的趋势,丰水年(P=10%)的总氮、总磷负荷分别为平水年(P=50%)的1.76倍和1.17倍,为枯水年(P=90%)的4.09倍和4.74倍.因此在丰水年应特别注意研究区非点源污染的控制与防治.1996和1997年各月的非点源污染负荷随逐月径流量的变化过程见图3.从图3可以看出年内的总氮和总磷含量随季节大致呈不规则的“W”形分布,一般在春汛期会出现一个小的峰值,进入汛期以后,特别是8月会出现一个较大的峰值,在冬季随着降雪和融雪的发生,又会出现一个较小的峰值.表3 不同水文年非点源污染负荷Table 3 Non-point source pollution loads in different hydrological years90 1996 411.2 37.05 2.54 2112.5 67.6 75 1970 496.2 32.53 5.64 4224.0 206.9 50 2000 557.8 55.39 3.91 4917.6 274.8 10 1988 764.4 93.86 32.83 8649.8 320.2图3 1996和1997年逐月流量实测值与河道总氮、总磷、输沙量变化Fig.3 Variations of observed monthly runoff,total nitrogen,total phosphorus and sediment from1996 to 1997从图3还可以看出,流域内径流量较大的汛期(6—9月)非点源污染的总氮和总磷负荷量均较大,而径流量少的枯水期非点源污染负荷量也较小.因为非点源污染常常是伴随着降雨径流过程特别是暴雨过程而产生的,所以非点源污染主要集中在汛期产生[25].汛期的非点源污染负荷贡献率见表4.汛期的总氮和总磷负荷量所占的比例较大,分别为85.83%和84.07%,这一时期的径流量和输沙量分别占全年总量的86.74%和95.69%,这与氮磷污染负荷在汛期内所占的比例是相近的,说明年内氮磷污染负荷与径流和产沙有关.通过研究可以发现研究区内的总氮负荷同月径流量呈明显的正相关关系,而总磷负荷同泥沙负荷呈明显的正相关关系,其二次拟合方程分别为表4 1993—2001年年均汛期非点源污染负荷贡献率Table 4 Contribution rates of annual average non-point source pollution during flood season from 1993 to 2001 %6 5.29 2.40 4.63 5.60 7 10.73 6.32 9.10 10.62 8 48.09 51.7167.49 74.86 9 22.63 25.40 2.85 4.61汛期总计 86.74 85.83 84.07 95.69式中:y1——月平均流量,m3/s;x1——总氮负荷量,t;y2——月泥沙负荷量,t;x2——总磷负荷量,t.可以得出结论,年内氮磷污染负荷与径流和泥沙负荷有很强的相关性,其中氮负荷同月径流的相关系数为0.848,磷负荷同泥沙负荷的相关系数为0.965.汛期是研究流域内非点源污染产生和发展的重要时段,这一时期的泥沙、总氮和总磷负荷量占全年的总负荷量的80%以上.因此汛期是防治研究流域非点源污染的主要时期. 3.2 非点源污染负荷的空间分布非点源污染有很强的空间差异性,可以根据研究区各子流域污染物负荷的大小来研究非点源污染的空间分布,进而找到对于非点源污染比较敏感的区域.本研究结合降雨和土地利用类型在流域内的空间分布情况,进行泥沙、有机氮、有机磷负荷空间分布的对比和原因分析.研究流域内降雨、泥沙、有机氮、有机磷负荷空间分布情况如图4所示.从图4可以看出,泥沙流失同降雨的空间分布具有相关性,流域内降雨较大的地区泥沙流失也比较严重.这是因为降雨是土壤侵蚀的主要驱动力,特别是暴雨的冲刷会造成土壤侵蚀加剧,进而造成泥沙的高负荷.有机氮负荷和泥沙负荷的空间分布很相似,相关系数为0.8925,具有很强的相关性.分析原因是颗粒物对有机氮有较强的吸附能力[26],有机氮吸附于泥沙而随泥沙输移.另外,有机氮负荷的空间分布同降雨也有一定关系,在降雨丰富的地区有机氮负荷也较大.有机磷负荷的空间分布同有机氮负荷的空间分布相似,二者的相关系数为0.9945,这说明影响二者空间分布的因素是一致的.此外,通过将流域内非点源污染分布图与流域内的土地利用图进行对比可以发现,河道附近农田较多的地区非点源污染负荷较大.图4 研究流域内降雨与泥沙、有机氮、有机磷负荷分布Fig.4 Distribution of rainfall and sediment,organic nitrogen,organic phosphorus loads in study area3.3 不同土地利用类型的非点源污染负荷通过以上分析可以发现除了降雨外,土地利用类型也是影响研究区非点源污染分布的主要因素,所以有必要对不同土地利用类型的非点源污染负荷量进行研究.研究流域内1993—2001年不同土地利用类型的年均非点源污染负荷量见表5.从表5可看出,研究流域内耕地的非点源污染负荷量最大,泥沙、总氮和总磷负荷的比例均超过总负荷的80%,是非点源污染的主要发生地;疏林地是流域内泥沙和磷负荷的第二大来源地;由于林地面积占整个流域的57.36%,所以林地的氮污染负荷总量也较大.研究流域内1993—2001年不同土地利用类型非点源污染单位面积负荷对比见表6.从表6可看出,不同土地利用类型的非点源污染负荷量不同,耕地的单位面积非点源污染负荷最高,疏林地次之,林地、草地等其他土地利用类型较小.这与李俊然等[27]的研究结论一致,即在单一土地利用类型占主导地位的流域内,土地利用类型以林地和草地为主的小流域地表水水质明显比以耕地为主的小流域地表水质好.因此控制研究区域内非点源污染的关键是控制耕地的氮磷污染和泥沙流失,采取退耕还林、等高耕作等水土保持措施[28],同时还要注意耕地中农药化肥的合理施用. 表5 不同土地利用类型1993—2001年年均非点源污染负荷模拟结果Table 5 Simulated results of annual average non-point source pollution loads of different land uses from 1993 to 2001林地 6049.34 57.36 27.85 10384.78 39.38耕地 3274.53 31.05 3803.44 89595.88 22472.71疏林地 422.30 4.01 77.81 3358.04 52.86草地 399.20 3.79 0.61 374.93 2.91滩地 128.78 1.22 4.81 444.85 20.75总计 10274.15 97.43 3914.52 104158.48 22588.614 结论表6 不同土地利用类型1993—2001年年均非点源污染单位面积负荷模拟结果Table 6 Simulated results of annual average non-point source pollution loads of unit area for different land uses from 1993 to2001林地 0.0046 1.72 0.01耕地 1.1615 27.36 6.86疏林地 0.1843 7.95 0.13草地 0.0015 0.940.01滩地 0.0373 3.45 0.16a.通过收集和统计研究流域内大量资料,建立了研究流域的水文模拟基础数据库.通过率定SWAT模型参数,得到研究流域逐月径流模拟值与实测值的Nash-Suttcliffe系数、相关系数均大于0.8,模型的模拟是符合实际的,经过参数率定的SWAT模型适用于该流域.b.在时间尺度上,研究流域内的农业非点源污染负荷与降雨量和径流量有较强的相关关系.年际间的非点源污染负荷的差异主要是由于降雨量不同造成的.丰水年(P=10%)的总氮、总磷负荷分别为平水年(P=50%)的1.76倍和1.17倍,为枯水年(P=90%)的4.09倍和4.74倍.年内非点源污染负荷总体上随季节呈不规则的“W”形态变化,汛期是研究流域内非点源污染产生和发展的重要时段,这一时期的泥沙、总氮和总磷负荷量占全年的总负荷量的80%以上.c.在空间尺度上,研究流域内降雨丰富的中上游地区的非点源污染负荷要大于流域下游地区的非点源污染负荷;流域内河流附近农田面积较多的地区非点源污染较为严重.降雨、产沙和土地利用类型是影响研究区非点源污染空间分布的主要影响因素.d.不同土地利用类型的非点源污染负荷不同,耕地的非点源污染负荷最高,泥沙、总氮和总磷负荷量占总负荷量的比例均超过80%,是非点源污染的主要发生地;单位面积非点源污染负荷最高的是耕地,疏林地次之,林地、草地等其他土地利用类型较小.参考文献:【相关文献】[1]VAN DER MOLEN D T,PORTIELJE R,DE NOBEL W T,et al.Nitrogen in Dutch freshwater lakes:trends and targets[J]. Environmental Pollution,1998,102(S1):553-557.[2]FONTAINE T A,CRUICKSHANK T S,ARNOLD J G,et al.Development of a snowfall-snowmelt routine for mountainous terrain for the soil water assessmenttool(SWAT)[J].Journal of Hydrology,2002,262:209-233.[3]AR NOLD J G,ALLEN P M.Estimating hydrologic budgets for three Illinoiswatersheds[J].Journal of Hydrology,1996,176:57-77.[4]AR NOLD JG,SRINIVASAN R,MUTTIAH R S,et al.Continental scale simulation of the hydrologic balance[J].Journal of the American Water Resources Association,1999,35(5):1037-1051.[5]ECKHARDT K,ARNOLD J G.Automatic calibration of a distributed catchmentmodel[J].Journal of Hydrology,2001,251(112):103-109.[6]ARNOLD JG,SRINIVASAN R,TSETAL rge area hydrologic modeling and assessment PartⅡ:model application[J].Journalof the American Water Resources Association,1998,43(1):91-101.[7]SALEHA,ARNOLD JG,GASSMAN P W,et al.Application of SWAT for the upper north Bosque Riverwatershed[J].Transactions of the American Society of Agricultural and Biological and Biological Engineers,2000,43(5):1077-1087.[8]SANTHI C,AR NOLD J G,WILLIAMS J R,et al.Validation of the SWAT model on large river basin with point and nonpoint sources[J]. Journal of the American Water Resources Association,2001,37(5):1169-1188.[9]DI LUZIO M,SRINIVASAN R,AR NOLD J G.A GIS-coupled hydrological model system for the watershed assessment of agricultural nonpoint and point sources ofpollution[J].Transactions in GIS,2004,8(1):113-136.[10]BOSCH N S.The infuence of impoundments on riverine nutrient transport:an evaluation using the soil and water assessment tool[J]. Journal of Hydrology,2008,355(1/2/3/4):131-147.[11]KANG M S,PAR K S W,LEE J J,ET Al.Applying SWAT for TMDL programs to a small watershed containing rice paddy fields[J]. Agricultural Water Management,2006,79(1):72-92.[12]BEHERA S,PANDA R K.Evaluation of management alternatives for an agricultural watershed in a sub-humid subtropical region using a physical process basedmodel[J].Agriculture Ecosystems&Environment,2006,113(1/4):62-72.[13]JURGEN S,KARIM CA,RAGHAVAN S,et al.Estimation of freshwater availability inthe WestAfrican sub-continent using the SWAT hydrologic model[J].Journal of Hydrology,2008,352(112):30-49.[14]CHAPONNIEREA,BOULET G,CHEHBOUNi A.Unterstanding hydrological processes with scarce date in a mountain environment[J]. Hydrological Processes,2008,22(12):1908-1921.[15]BARLAND I,KIRKKALA T.Examining a model and assessing its performance in describing nutrient and sediment transport dynamics in a catchment in southwestern Finland[J].Boreal Environment Research,2008,13(3):195-207.[16]胥彦玲,秦耀民,李怀恩,等.SWAT模型在陕西黑河流域非点源污染模拟中的应用[J].水土保持通报,2009,29(4):114-117.(XU Yan-ling,QIN Yao-min,LI Huai-en,et al.Applicationof SWAT model to non-point source pollution simulation in Heihe River basin of Shaanxi Province[J].Bulletin of Soil and Water Conservation,2009,29(4):114-117.(in Chinese)) [17]张永勇,王中根,于磊,等.SWAT水质模块的扩展及其在海河流域典型区的应用[J].资源科学,2009,31(1):94-100. (ZHANG Yong-yong,WANG Zhong-gen,YU Lei,etal.Extendedwater quality module ofSWAT model andits applicationto Hai Riverbasin[J].ResourcesScience,2009,31(1):94-100.(in Chinese))[18]庞靖鹏,刘昌明,徐宗学.密云水库流域土地利用变化对产流和产沙的影响[J].北京师范大学学报:自然科学版,2006,46 (3):290-299.(PANG Jing-peng,LIU Chang-ming,XU Zong-xue.Impact of land use change on runoff and sediment yield in the Miyun Reservoir catchment[J].Journal of Beijing Normal University:Natural Sciences,2006,46(3):290-299.(in Chinese))[19]杨育红,阎百兴.中国东北地区非点源污染研究进展[J].应用生态学报,2010,21(3):777-784.(YANG Yu-hong,YAN Baixing.Nonpoint source pollution in Northeast China:research advance[J].Chinese Journal of Applied Ecology,2010,21(3):777-784. (in Chinese)). [20]任友山.蚂蚁河流域水文特性分析[J].黑龙江水利科技,2008,36(4):102-103.(REN You-shan.Hydrological characters analysis of Mayihe river basin[J].Heilongjiang Science and Technology of Water Conservancy,2008,36(4):102-103.(in Chinese))[21]黑龙江省土地管理局,黑龙江省土壤普查办公室.黑龙江土壤[M].北京:农业出版社,1992.[22]魏怀斌,张占庞,杨金鹏.SWAT模型土壤数据库建立方法[J].水利水电技术,2007,38(6):15-18.(WEI Huai-bin,ZHANG Zhan-pang,YANG Jin-peng.Establishing method for soil database of SWAT model[J].Water Resources and Hydropower Engineering,2007,38(6):15-18.(in Chinese))[23]李硕.GIS和遥感辅助下流域模拟的空间离散化与参数化研究与应用[D].南京:南京师范大学,2002.[24]郝芳华,程红光,杨胜天.非点源污染模型:理论方法与应用[M].北京:中国环境科学出版社,2006:54-55.[25]李国斌,王焰新,程胜高.基于暴雨径流过程监测的非点源污染负荷定量研究[J].环境保护,2002(5):46-48.(LI Guo-bin,WANG Yan-xin,CHENG Sheng-gao.A quantitative study on non-point source pollutant loading based on stream runoff monitoring results[J].Environmental Protection,2002(5):46-48.(in Chinese))[26]张学青,夏星辉,杨志峰.水体颗粒物对有机氮转化的影响[J].环境科学,2007,28(9):1954-1959.(ZHANG Xue-qing,XIA Xing-hui,YANG Zhi-feng.Effect of suspended sediment on the transformation of organic nitrogen[J].Environmental Science,2007,28 (9):1954-1959.(in Chinese))[27]李俊然,陈利顶,郭旭东,等.土地利用结构对非点源污染的影响[J].中国环境科学,2000,20(6):506-510.(LI Jun-ran,CHEN Li-ding,GUO Xu-dong,et al.Effect of land use structure on non-point source pollution[J].China Environmental Science,2000,20(6):506-510.(in Chinese))[28]陈雪,蔡强国,王学强.典型黑土区坡耕地水土保持措施适宜性分析[J].中国水土保持科学,2008,6(5):44-49.(CHEN Xue,CAI Qiang-guo,WANG Xue-qiang.Suitability of soil and water conservation measures on sloping farmland in typical black soil regions of Northeast China[J].Science of Soil andWater Conservation,2008,6(5):44-49.(in Chinese))。
东辽河辽源段非点源污染估算与变化规律分析

东辽河辽源段非点源污染估算与变化规律分析李杰;余麟;许君雨;张利静;王宪恩;段海燕【摘要】控制非点源污染排放成为东辽河辽源段水污染治理的重要环节.采用枯水期日均污染负荷对径流分割法进行了修正,将其应用于东辽河辽源段非点源COD和氨氮负荷的估算,并运用Pearson相关系数分析非点源污染负荷与当地降水量之间的相关性.研究表明,东辽河辽源段非点源污染具有年际变化大、年内分布不均的特点,非点源COD占比为50%~ 65%,非点源氨氮占比约为30%~50%,该河段非点源污染严重.东辽河非点源污染年内月际分布不均、月际变化大,COD主要集中在7、8月,主要由夏季大量降水引起;氨氮则有3月和7月两次高峰,主要是由于污染物随春季融雪入河和夏季汛期与水冲刷共同导致;非点源污染与降水量具有较强的相关性.【期刊名称】《科学技术与工程》【年(卷),期】2013(013)029【总页数】6页(P8691-8696)【关键词】东辽河;非点源污染;变化规律分析;径流分割法;Pearson相关系数【作者】李杰;余麟;许君雨;张利静;王宪恩;段海燕【作者单位】吉林大学地下水资源与环境教育部重点实验室,长春130021;吉林大学地下水资源与环境教育部重点实验室,长春130021;吉林大学地下水资源与环境教育部重点实验室,长春130021;吉林大学地下水资源与环境教育部重点实验室,长春130021;吉林大学地下水资源与环境教育部重点实验室,长春130021;吉林大学地下水资源与环境教育部重点实验室,长春130021【正文语种】中文【中图分类】X522近年来,随着流域内经济社会的不断发展,东辽河水污染较为严重[1,2],已经成为影响吉林省和下游辽宁省经济与社会可持续发展的重大制约性因素。
随着工业污染的治理和城镇污水处理厂的投入使用,非点源污染逐渐成为影响东辽河水环境质量的主要因素,由此,控制非点源污染排放成为东辽河辽源段水污染治理的重要环节。
城市水源地非点源污染研究综述

摘
要:随着城市水源地点源污染控制水平的不断提高 ,非点源污染也越来越受到人们的重视 。文中在阐述非点源
污染 概念 的基 础 上 ,概述 了国 内外城 市 水源 地 非点源 污 染 的研 究现 状 ,分 析 了 国内外 研 究存在 的问题 ,并讨 论 了城 市水 源地 非点 源污 染 研究 的未 来 发展 趋势 。 关键 词 :非点 源污 染 ;城 市水 源地 ;发展 趋势
荇染。可见非点源污染已经成为城市水源地的重要污染源,
对 非点 源 污 染 研 究 也 已成 为 国 际 上环 境 问题 研 究 的 活 跃领
域。
一
7 0年代我国城市水源地水环境研究的重点是 由点源、 非 点源造成的湖泊富营养化 ,但我国真正意义上的非点源污染
研 究 是从 北 京城 市 径流 污染 研 究 开始 的 ,之 后相 继 在 上海 、
第 1 卷 第 6期 1
21 0 1年
中 国
水
运
V .1 o1 1 Jn ue
N 6 o. 2 1 01
6月
C n Wa er hi a t Tr n or a sp t
城市水源地 非点源污染研 究综 述
徐 洪 庆 ,宋 晓 冉 ,杜 子 龙 。
( 1山东 省胶 东调 水 局 青 岛分 局 , 山东 青 岛 2 60 ; 600 2青 岛水 文水 资 源勘 测局 ,山东 青 岛 26 1 ;3山东 省水 文 水 资源勘 测 局 , 山东 济 南 20 1 ) 60 7 50 4
2 国 内研 究现 状 .
以至我国城市水源地非点源污染问题 日益突出。研究表明 ,
即使 点源 污 染得 到全 部 控制 ,水源 地 水环 境 也得 不 到根 本 的 改善 ,水质 达标 率 仅为 4 % ̄6 %i 2 5 ,北 京密 云水 库 、天津 于 桥 水库 、 巢湖 和 太湖 等水 域 的 非点源 污 染 比例 已超 过 点源
非点源污染负荷的估算方法在TMDL计划中的应用

5科技资讯科技资讯S I N &T NOLOGY I NFORM TI O N 2008N O .11SC I ENC E &TEC HN OLO GY I NFO RM ATI O N 资源与环境与点源污染相比,非点源污染的发生受到大气、土壤、地形、植被、水文等环境因素和人类生产、生活等行为的多重影响,具有在不确定时间内通过不确定途径排放不确定数量污染物的特性,这一特性决定了非点源污染负荷定量化估算的较大难度,以及非点源污染控制管理任务的长期性和艰巨性。
1TM D L 概况美国环保局(EP A )于1972年提出了TM DL (t ot a l m axi m um da i l y l oa d ,最大日负荷总量)计划。
这比排放许可证制度更科学合理,因为它把水资源保护由过去单纯的污染源控制逐渐转变为根据生态健康和生态功能来决定控制污染。
美国已实施的一些案例表明T M D L 计划能够改善水体质量,尤其在点源和非点源污染综合控制方面成效显著。
这一水污染控制模式已在全球多个国家和地区推广实施,也将成为今后国内借鉴运用的主要方法。
所谓T MD L ,即水体在满足当前适用的水质标准前提下所能接受的某种污染物的最大日负荷量。
T M D L 计划的目标之一就是将可分配的污染负荷分配到各个污染源,包括点源和非点源,同时要考虑安全临界值和季节性的变化,进而采取适当的污染控制措施。
一般地,T M DL =∑W L As +∑LAs +M OS 式中:WL A :包括现存和未来点源的污染负荷;L A :包括非点源污染负荷以及自然本底部分;M O S :安全因数,亦称安全临界值,指关于污染物质负荷与受纳水体水质之间关系的不确定数量,有确定和不确定关系两种说法,可通过假设分析提供一个不确定的数量比例关系,或把它当作负荷容量的储备部分而明确地给出,如负荷容量的5%或10%,等等。
2国外T M D L 中常用的非点源污染负荷方法T MDL 方法中最关键的便是对非点源污染负荷的定量化估算。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
流域非点源污染负荷估算模型的研究
2
———————————————————————————————— 作者: ———————————————————————————————— 日期:
3
流域非点源污染负荷估算模型的研究-工程论文 流域非点源污染负荷估算模型的研究 杜娟① DU Juan;赵湘璧① ZHAO Xiang-bi;李怀恩② LI Huai-en (①榆林学院建筑工程系,榆林 719000;②西安理工大学西北水资源与环境生态教育部重点实验室,西安 710048) (①Department of Architectural Engineering,Yulin University,Yulin 719000,China; ②Key Lab of Northwest Water Resource and Environment Ecology of MOE,Xi′an University of Technology,Xi′an 710048,China) 摘要:非点源污染是目前影响水环境质量的主要污染源之一,本文对国内外水环境非点源污染负荷的估算方法进行了总结和分析,介绍了污染分割法、输出系数法、平均浓度法等经验统计方法和ANSWERS、BASINS、SWAT等分布式模型估算方法的特点。在此基础上,提出了如何选用合适的非点源污染负荷估算方法。 Abstract: Non-point source pollution was a major source that influenced water quality. Based on the comprehensive literature review at home and abroad, a review of recent research advancement in non-point source pollution loads estimation methods was presented. The estimation methods could group into two types of experiential statistical method and model estimation method. The principles and applications of experiential statistical methods, such as Pollution Division, export coefficient model, mean concentration method, were
4
introduced, and the characteristics and applications of major estimation models, such as Areal Non-point Source Watershed Environment Response Simulation(ANSWERS), Better Assessment Science Integrating Point and Non-point Sources(BASINS) and Soil and Water Assessment Tools(SWAT) were summarized and analyzed. Finally, the method of estimating the load of non-point source pollution was put forward. 关键词 : 非点源污染;负荷估算;模型;研究 Key words: non-point pollution;load estimation;model;research 中图分类号:X820.3;F301 文献标识码:A文章编号:1006-4311(2015)25-0247-02 基金项目:榆林学院高层次人才启动基金项目14GK31。 作者简介:杜娟(通讯作者)(1987-),女,陕西榆林人,硕士研究生,研究方向为非点源污染与水资源保护。 0 引言 水环境污染可分为点源污染(PSP)和非点源污染(NSP),点源污染易于识别和治理,且随着PSP控制能力的提高,NSP的严重性逐渐显现出来。研究表明:即使PSP得到全面控制,水环境也得不到根本改善,水质达标率仅达42%~65%[1]。 我国属农业大国,随着人口的持续增加以及人类不合理活动的增多,更加剧了水体NSP,导致可用水资源量愈加减少。因此,加强对NSP的研究对解决人类面临的水危机具有重要意义。
5
1 非点源污染负荷估算模型 目前针对流域NSP的估算、模拟有多种模型和方法。总体来说有两种途径:①通过分析水文、水质数据,用径流分割法来估算流域NSP的输出量;②分析污染源与控制断面之间的响应关系建立NSP模拟模型,进而估算NSP负荷。 常用的NSP模型分为两大类:集总式模型和分布式模型。 有代表性的集总式模型有:降雨量差值法[2]、相关关系法[3]、平均浓度法[4]、污染分割法[5]、水质水量法[6]、输出系数法[7]等。 降雨量差值法[2]认为:NSP的产生受降雨量和降雨径流过程的影响,其负荷与降雨量的大小密切相关。任意两场洪水(或任意两年)产生的污染负荷(包括点源和非点源)之差应为这两场(或这两年)降雨量之差引起的NSP负荷,即: LA-LB=f(RA)-f(RB)=f(RA-RB)=LnA-LnB(1) 在暴雨径流过程中,N、P等颗粒态污染物的迁移与泥沙输移有密切关系,可以建立NSP负荷-泥沙相关关系[3]。 平均浓度法[4]根据各次降雨径流过程的水量、水质同步监测资料,计算各污染物NSP的平均浓度,以各次暴雨的径流量为权重,求加权平均浓度作为地表径流的平均浓度,与地表径流之积为NSP负荷量。 污染分割法[5]以流域出口断面的径流、水质及径流在年内的分配为基础,进行流域的PSP和NSP负荷的分割,属水文学方法,机理明确,资料易得。 水质水量法[6]根据典型流域各次暴雨水质水量同步监测资料,忽略地表径流迁移过程,利用统计分析方法,建立水质水量之间的相关关系。 输出系数法的一般表达式[7]如下: 避开了NSP发生和发展的复杂过程,且所需参数少,操作简便易行,又具有一定的精度,适合水文水质资料缺乏的地区
6
以及大尺度流域,因此得到了广泛的应用。 集总式水文模型(包括经验类模型)主要是通过回归分析的手段建立污染源和监测断面之间的关系,得出污染源合理的相关输出系数,估算整个流域的污染物输出量。 这类模型的精度较低,优点是对数据的需求也比较低,能够简便计算出流域出口或区域单元总污染负荷,表现出较强的实用性和一定的准确性,但由于难以描述污染物迁移的路径与机理,使模型应用受限,且所取经验系数受限于区域,阻碍了系数的通用性和模型的可转移性。因此,这类模型的关键在于如何合理地确定相关经验系数。 在近年NSP负荷的研究中,蔡明[8]对输出系数法进行了改进,利用长年的资料加入了降雨因子和流域损失对于经验系数的影响。 杜娟[9]利用连续多年的实测资料以及两种方法估算结果的比对中,考虑到输移损失的影响,对于已建模型进行修正系数的重新确定,使其更符合流域实际情况。 分布式模型常见的有ANSWERS[10]、 BASINS[11]、 SWAT[12]、GLEAMS[13]、AnnAGNPS等。主要是对流域产汇流过程、水土流失过程和污染物迁移转化过程进行机理层次上的模拟,估算污染源的强度和输出负荷量。NSP模型一般是由四个子模型系统构成的,分别是:降雨径流模型、侵蚀和泥沙输移模型、污染物转化模型、受纳水体水质模型。 分布式模型虽然有较好的耦合作用,模拟精度较高,也可在后期加入LID以及BMPs措施进行具体的流域模拟,但模型前期所需数据庞大(如长期的气象、地
7
形、地表资料,及实测的长期水文、水质监测数据),且建模时要求拟合程度高(需进行反复多次的率定、验证),才有模拟流域(区域)的可行性,估算成本较高,但估算结果的实际意义较大。故而在应用分布式模型模拟流域状况时,关键是采用长期的动态跟踪以获得上述的各类参数。 2 结论 综上所述,在对流域(区域)NSP负荷估算时,可同时采用集总式水文模型(包括经验类模型)和分布式模型,互相检验,既可修正经验类模型的相关系数,又有利于模型率定和验证时参数的调整。但无论选取哪一种模型对于NSP负荷进行估算,都必须有实测的水质水量资料作为依据进行检验和校正,以得到更为精确的NSP负荷。 参考文献: [1]宫莹,阮晓红,胡晓东.我国城市地表水环境非点源污染的研究进展[J].中国给水排水,2003(19). [2]蔡明,李怀恩,庄咏涛.估算流域非点源污染负荷的降雨量差值法[J].西北农林科技大学学报:自然科学版,2005,33(4). [3]李怀恩,蔡明.非点源营养负荷-泥沙关系的建立及其应用[J].地理科学,2003,23(4). [4]李怀恩.估算非点源污染负荷的平均浓度法及其应用[J].环境科学学报,2000,20(4):397-400. [5]施为光,凌文州.用实测资料计算流域非点源污染负荷[J].长江流域资源与环境,1996,5(3). [6]李国斌,王焰新,程胜高.基于暴雨径流过程监测的非点源污染负荷定量研