第9卷第5期 2011年1O月 南水北调与水利科技 South-to-North Water Diversion and Water Science&Technology Vo1.9 No.5 Oct.2011 doi:10.3724/SP.J.1201.2011.05055 基于SWMM的不同降水量 对城市降雨径流TSS的影响分析 姜体胜L ,孙艳伟。,杨忠山 ,黄振芳 ,张 蕾 ,焦中志 (1.北京市水文总站,北京100089;2.北京市水环境监测中心,北京100089; 3.华北水利水电学院水利学院,郑州450011) 摘要:受物理、化学、生物等多重因素的影响,城市降雨径流水质的模拟研究比较模糊。以TSS(总悬浮固体量)为 例,利用SWMM的污染物模块及其土地利用模块、污染物累积模块和污染物冲刷模块进行城市雨水径流水质模 拟;考虑了不同土地利用情况对污染物的累积及冲刷的影响,针对不同的降水量,利用EMC法(平均浓度法)以及 指数方程法进行污染物的冲刷模拟,指数方程法的模拟精度大于EMC法的模拟精度,与监测值更接近。指数方程 法的模拟结果表明:当降水量小于2.5 mm时,其TSS冲刷量增长幅度较缓;而当降水量大于3.5 mm时,其TSS冲 刷量呈显著线性增加关系,从而表明,降水量较大的降水事件对于TSS的冲刷量影响更为显著。 关键词:暴雨水管理模型;降雨径流水质;总悬浮固体量;非点源污染;平均浓度法;指数方程法 中图分类号:TV122 文献标识码:A 文章编号:1672—1683(2011)05—0055—04 Simulation of Urban Runoff Quality Based on Stormwater Management Model JIANG Ti—sheng ,SUN Yan-wei。,YANG Zhong-shan ,HUANG Zhen-fang ,ZHANG Lei ,JIAO Zhong-zhi ’ (1.Beijing Hydrologic Center,Beijing 100089,China;2.Beijing Water Environment Monitoring Center, Beijing 100089,China;3.North China University of Water Resources and Electric Power,Zhengzhou 450011,China) Abstract:Due tO influences of physical,chemical and ecological factors,simulation of the urban runoff quality is a very complex issue in the study of wet-weather flows and their environmental impacts.Taking total suspended solid(TSS)as an example,this paper simulated the building-up and washing-off TSS using the modules of pollutant,land use.build-up and wash-off of storm— water management model(SWMM).Different land uses were taken into consideration for analysis.The event mean concentra— tion(EMC)and ex.mnential function were used to simulate the wash-off effects.Calibration process indicated that the exponen— tial method was easi’er to calibrate than EM[C.Simulation results of the exponential function indicated that the TSS had a slow increase when the rainfall depth is less than 1.25cm;however.there was a significant linear increase when the rainfall depth is greater than 0.35 cm,which showed that greater rainfall had more significant impact on TSS wash-off rate. Key words:stormwater management model;runoff quality;total suspended solid;non-point source pollution;event mean concen— tration;exponential function 城市雨水径流水质受到物理、化学、生物等多方面因素 水体,从而给水体功能带来负面影响[5 ]。由于非点源污染 的产生以及污染物的迁移受到众多环境要素以及人类活动 的影响,加之长期水质监测资料的缺乏,使其水质的模拟研 究比较模糊[1 ]。随着点源污染控制的不断完善和城市的快 速发展,非点源污染的问题日益突出,对城市水系构成严重 的威胁『3 ]。非点源污染是指溶解的或固体污染物(包括泥 一的影响,因此利用模型对非点源污染进行定量模拟研究存在 定的困难。 EPA SWMM(Environment Protection Agency Storm 沙、养分、细菌、需氧物质、重金属、微量有机物等)从非特定 的地域在降水和径流冲刷作用下,通过径流过程汇入到受纳 Water Management Model,简称SWMM)是美国环境保护署 研发的一套针对雨水资源管理的软件,该软件能动态模拟长 收稿日期:2010—10—24 修回日期:2011-08—08 网络出版时间:2011一Og一04 网络出版地址:http://www.cnki.net/kcms/detail/13.1334.T 20110904.0821.035.html 基金项目:北京市优秀人才培养基金(PYZZ090419001989);北京市科技计划项目(D07050601510000) 作者简介:姜体胜(1982一),男,山东成武人,博士,主要从事水环境水资源方面的研究。E-mail:jiangtisheng@126.corn 通讯作者:孙艳伟(1983一),女,山东聊城人,博士,主要从事水资源高效利用、生态水文学以及低影响发展研究。E-mail:stmyanwei@}yahoo.o。m cn 55・ ・第9卷总第56期・南水北调与水利科技・2011年第5期 期和短期的降雨径流过程,并能模拟雨水径流的流量及水 质,因而得以广泛应用_7]。SWMM通过一系列可以对参数 进行校核的经验方程,可以对非点源污染的发生及迁移排放 过程进行动态模拟。本文以TSS(总悬浮固体量)为例,利用 SWMM模拟某一城市小区非点源污染的发生及排放过程, 考虑了不同土地利用情况对污染物的产生及冲刷的影响,针 对不同的降水量,利用降雨径流事件平均浓度(Event Mean Concentration,简称EMC)法以及指数方程法进行污染物的 冲刷模拟。通过对城市小区域系统及其内部发生的复杂污 染过程进行定量描述,可以分析不同降水强度下城市降雨径 流水质时空变化特征,识别其主要来源和迁移路径,预测预 报污染负荷及其对水体的影响,从而为城市水文水资源管理 2.1 降水事件 对长序列降水事件的分析可知,降水量较小的降水事件 占降水序列的绝大部分,而且是城市雨水径流量以及污染物 冲刷量最主要的组成部分[9]。通过对研究区降水系列进行 分析,该区2年及1O年重现期的设计暴雨深度分别为0.25 em和0.65cm。为了分析不同的降水量对污染物的累积及 冲刷的影响,本文采用均匀内插的方法,分别模拟降水量为 0.5 mIn、1.5 mill、2.5 mm、3.5 mm、4.5 mm、5.5 mIn和6.5 mlTl历时2 h的降水事件。对研究区的降水数据进行分析发 现,其降水类型呈现明显的SCSII型降水分布,因此本文采 用SCSII型降水分布进行模拟l】 。 提供决策依据。 l基于SWMM的非点源污染模型 SWMM认为污染物的产生及累积过程主要发生在无降 雨时段,并通过降雨及径流的冲刷作用实现污染物的迁移及 排泄。SWMM通过一系列的模块及经验方程来模拟城市雨 水径流水质。其中,SWMM通过雨量站来模拟不同的降水 事件;利用子区域模块模拟区域的产流、汇流及入渗过程的 模拟;污染物模块主要用来定义拟模拟的污染物,一般来讲, 子区域模块是污染物的最主要来源,但是,SWMM也可以通 过直接添加的形式模拟其他污染物的来源,如降水、地下水、 入渗、汇流源头等。土地利用模块可用来模拟不同土地利用 情况,如居民区、商业区、工业区等对污染物的产生及冲刷的 影响;污染物累积模块用来定义污染物在特定的土地利用情 况下无降水时段的累积速率,SWMM采用幂函数、指数函数 以及饱和函数三种不同的函数来定义不同污染物在不同土 地利用情况下的产生速率;污染物冲刷模块主要用来定义在 降雨及径流冲刷作用下,区域表面污染物的侵蚀、溶解过程, SWMM采用指数方程法、关系曲线法和平均浓度法计算污 染物的冲刷速率。SWMM进行城市雨水径流非点源污染物 的模拟流程如图1所示。 污染物模块和土地利用模块 污染物累积模块J J污染物冲刷模块 和函数法l l指数函数法I l幂函数法l l指数方程法ll关系曲线法ll平均浓度法 图1基于SWMM的城市雨水径流非点源污染物的模拟流程 Fig.1 Simulation flowchart of urban runoff quality based on SWMM 2模型的建立及主要参数 TSS是城市雨水中浓度较高的最常见污染物,雨水径流 中过高的TSS会引起下游河道水环境的劣变,如水生生物 栖息环境的变化、河道浑浊以及河道景观及娱乐性能的丧失 等[8]。本文以TSS为对象,利用SWMM的污染物累积和冲 刷模块,对北京市某小区的雨水径流水质进行模拟研究。其 中TSS浓度单位为mg/L,并假设降水中所含的TSS为 10 mg/L。 .56・ 。2.2模拟区域 在SWMM的基础上,通过分析北京市面积为0.1 km2 的某一小区TSS对不同降水量的响应,讨论城市非点源污 染的模拟,揭示不同降水量对污染物的产生及冲刷的影响。 研究区不透水面积比例为50 ,为典型的城市小区不透水面 积比例。区域各输入参数见表1。 表1 SWMM模拟区域的主要输入参数 Table 1 Major input subcatchment parameters of SWMM 参数名称 数值 面积/km 0.1 不透水面积比例( ) 50 区域宽度/m 600 平均坡度( ) 2 不透水区域的积水深度/mm 1.5 透水区域的积水深度/mm 12.7 2.3土地利用 通过SWMM水质菜单下的土地利用模块添加土地利 用情况,共定义三种不同的土地利用情况:①稀疏居民区; ②密集居民区;③商业区。其中,稀疏居民区的面积占 2o ,密集居民区的面积占5o ,商业区面积占3O 。 2.4污染物累积 选用指数函数法模拟污染物的累积过程。指数函数法是 污染物产生应用最为广泛的一种方法,用来描述污染物产生 速率达到最大值之前,与时间的指数函数关系,如式(1)。 B=C1(1--e一 2 ) (1) 式中:C1一单位面积或单位长度上污染物最大产生量(kg/km2 或kg/m);C2一污染物产生速率参数(1/d)。 根据文献资料El13,商业区和居民区所产生的TSS基本 相同,而人口密度大的居民区要比人口密度小的居民区所产 生的TSS大。利用该小区的监测TSS数据对指数方程法的 各项参数进行校核,使其误差控制在5 以内,得到不同土地 利用情况的c 和C2值见表2。 2.5污染物冲刷 选用EMC法和指数方程法两种方法模拟污染物的冲刷 过程。美国EPA提出“降雨径流事件平均浓度”(EMC),用 姜体胜等・基于SWMM的不同降水量对城市降雨径流TSS的影响分析 表2不同土地利用类型下污染物产生指数方程参数 Table 2 The exponential function parameters of olplutions buildir】g up with different land use types 持不变,而当降水量大于1.5 n'llTl时,随降水量增加呈上升 趋势。因此,降水量较小的降水事件对模拟过程中地表所产 生的TSS的影响不大。 来表示在一场降雨径流全过程中某污染物的平均浓度。 EMC法采用TSS的平均浓度进行模拟,即污染物的冲刷速 率在整个模拟过程中保持不变。由于TSS的平均浓度可以 通过实验获得,因此EMC法是SWMM及其他水质模型中 应用最为广泛的方法[12-13]。本文所采用的TSS的EMC值 见表3。 表3不同土地利用类型TSS的EMC值及冲刷参数 Table 3 The value of EMC and parameters of wash off with different land use types 参数 土地利用类型 EMC法 ∞ 指数方程法参数 9 8 7 6 5 4 3 2 1 0 平均浓度值/(rag・L ) F1 Fz 指数方程法相对EMC法较为复杂,用以描述污染物冲 刷速率与区域剩余污染物总量及流量的指数乘积之间的线 性关系。 W=F1・ ・B (2) 式中:w一污染物冲刷速率(kg/h);F 一污染物冲刷系数; g一单位面积内的径流量(mm/h);F2一冲刷指数;B一区域 污染物剩余总量(kg)。 对于指数方程法,根据泥沙运移理论,冲刷指数F2的值 介于1.1和2.6之间[】 ,考虑到密集居民区及商业区的不透 水性比例较高,因此,经校核,对于密集居民区及商业区,F2 值定为2.2,对于稀疏居民区,其值采用1.8。模型针对不同 土地利用类型所使用的参数见表3。 3模拟结果分析 3.1 TSS累积量分析 根据污染物冲刷模型的不同,采用上文的模型参数,分 别建立EMC法SWMM模型和指数方程法SWMM模型,其 中模拟总时间为12 h。由于降水事件持续时间较短,为了增 加模型的精确性,设定模拟的时间步长为1 rain。根据对降 水序列的分析,设定模拟前的干旱天数为8 d。由于各降水 事件模拟前的干旱天数均为8 d,而模型模拟前所产生的 TSS与降水强度无关,因此各降水事件模拟前区域的TsS 保持不变,均为330.2 。各降水事件模拟过程中所产生的 TSS、降水中携带的TSS、通过入渗损失的TSS见图2,由图 可知,TSS与降水入渗和降水量有直接关系,随着降水量增 加,降水中携带的TSS和入渗损失的TSS均呈直线上升趋 势。当降水量小于1.5 n21Tl时,模拟过程中所产生的TSS保 0 1 2 3 4 5 6 7 降水量/mm 图2不同降水量在模拟过程中的TSS累积量 Fig.2 Accumulation of TSS in different rainfall depth 3.2冲刷量分析 当采用EMC法进行冲刷模拟时,由于污染物的冲刷速 率在整个模拟过程中保持不变,因此,TSS浓度与降水量无 关,在整个降水径流过程中保持不变,为196 mg/L。与 EMC法冲刷模型不同,指数方程法冲刷模型中的TSS浓度 与降水强度及区域剩余的TSS总量有关,见图3。由图3可 知,TSS冲刷浓度随着降水量的增加而增加,其冲刷浓度的分 布与降水历程相符合,即经历了由小到大再减少至0的过程。 5∞ 彻 首 300 器 200 誉 盘 100 0 0 Z 4 6 8 LU I2 时间/h 图3采用指数方程冲刷模型时不同降水量的TSS浓度 Fig.3 Concentration of TSS of exponential function under different rainfall depth 区域TSS在降水径流作用下的冲刷总量,实际上为各 个时段TSS冲刷浓度与径流乘积的总和。为了分析EMC 法和指数方程法两种冲刷模型的应用效果,将其模拟结果与 在降水事件中监测的TSS总量进行分析比较,见图4。由图 4可知,指数方程法的模拟精度大于EMC法的模拟精度,与 监测值更接近。在对模型参数进行反复校核的过程中发现, 当采用EMC法进行冲刷模拟时,针对不同的降水事件采用 不同的平均浓度值时,可以获得较好的模拟精度,而在对降 水序列进行长期的模拟时,由于降水量的年内以及年际分布 的不均匀性,针对每次降水事件均改变EMC的值是不现实 的。在对指数方程法各参数进行校核时,发现指数方程法与 EMC法比较可以更容易达到模型所需精度。采用指数方程 法对不同降水强度下的TSS冲刷量进行分析,结果表明:降 ・57・ … … 第9卷总第56期・南水北调与水利科技・2011年第5期 水量越大,所产生的TSS量越大;当降水量小于2.5 nlm时, 其TSS冲刷量增长幅度较缓;而当降水量大于3.5 inin时, 其TSS冲刷量呈显著线性增加关系,简言之,降水量较大的 ZHU Cheng-li.Water Eco-environmenta[Construction and Pro— tection for Urban River[J].Journal of Irrigation and Drainage, 2004,23(6):18—21.(in Chinese)) 降水事件对于TSS的冲刷量影响更为显著。 80 [4]沈雪娇,杨侃,常蒲婷,等.改进的径流一水质非点源污染预测 模型研究[J].水电能源科学,2009,27(5):40—42.(SHEN Xue- jiao,YANG Kan,CHANG Pu-ting,et a1.State Key Laboratory 70 of Hydrology-Water Resources and Hydraulic Engineering[J]. Water Resources and Power,2009,27(5):40—42.(in Chinese)) 60 的50 E53赵其国,徐梦洁,吴志东.东南红壤丘陵地区农业可持续发展研 咖 蓓4O 量 30 究[J].土壤学报,2000,37(4):433—442.(ZHAO Qi—guo,XU Meng—jie.WU Zhi dong.The Agricultural Sustainable Develop— ment Research on Red Soil Upland Region in Southeast China 20 Agriculture Sustainable Development[J].Acta Pedologica Sini— ca,2000,37(4):433—442.(in Chinese)) 10 [6]常蒲婷,杨侃,侯学勇,等.突发性水污染事件模拟分析研究 0 0 5 1 5 2.5 3 5 4.5 5,5 6.5 EJ].水电能源科学,2009,27(3):44—47.57.(CHANG Pu—ting, YANG Kan,HOU Xue-yong,et a1.State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering,Hohai 降水量/arm 图4降水过程中的TSS冲刷量模拟结果与监测值对比 Fig.4 Comparison between the simulated value and observed value University[J].2009,27(3):44—47.57.(in Chinese)) [7]Lewis A.Rossman.Storm water management model user’s manual Version 5.0.Water Supply and Water Resources Divi— 4结语 本文以TSS为例,阐述了在没有特殊的水质处理措施 下,如何应用SWMM模型进行城市雨水径流水质模拟的问 sion,National Risk Management Research Laboratory,Cincin— nati,OH.(2009) 题。本文分别利用SWMM污染物冲刷模块下的EMC法和 指数方程法模拟了TSS的冲刷效果。结果表明,指数方程法 的模拟精度大于EMC法的模拟精度,与监测值更接近;降水 [81夏青.流域水污染物总量控制[M].北京:中国环境科学出版 社,1996.109—110.(XIA Qing.Total amount control of water pollution[M].Beijing:Chinese Environmental Science Press, 1996.109—11O.(in Chinese)) 量较大的降水事件对于TSS的冲刷量影响更为显著,降水 量越大,TSS累积量越大,TSS冲刷量也越大;当降水量小于 [9]Urban Drainage and Flood Control District(UDFCD)(2001).Ur— ban Storm Drainage Criteria Manual,2007 revision.Denver,0 2.5 Iylln时,其TSS冲刷量增长幅度较缓;而当降水量大于 3.5 ITim时,其TSS冲刷量呈显著线性增加关系。 参考文献(References): [13李贵宝,尹澄清,单宝庆.非点源污染控制与管理研究的概况与 [1o3 http nNw.1mnoeng.com/RainfallMaps/RainfallMaps.htm #Rainfall Distribution Graph [1】]Manning M.J,Sullivan R.H,Kipp T.M.(1977).Nationwide Evaluation of Combined Sewer Overflows and Urban Storm— water Discharges—Vo1.III:Characterization of Discharges. 展望_J].农业环境保护,2001,20(3):190—191.(LI Gui~bao, YIN Chen—qing,SHAN Bao—qing.Current Situation and Pros— EPA-600/2—77—064c(NTIS PD2721 07)。U S Environmental Protection Agency。Cincinnati。OH. pect of Control and Management on Non-Point Source Pollution r12]Field,R.and Tafuri A N.(1973).“Storm Pollution Control in U.S”.Combined Sewer Overflow Seminar Paper ̄EPA 670/2— 73—077,U S Environmental Protection Agency。Eclison,NJ. EJ].Agro—environmental Protection,2001,2O(3):190—191.(in Chinese)) Eel贺缠生,傅伯杰,陈利顶.非点源污染的管理及控制EJ].环境科 学,1998,19(5):87—96.(HE Chan-sheng,FU Bo-jie,CHEN Li- [13] Brezonik P.L and Stadelmann T.H_Analysis and predictive models of stormwater runoff volumes,loads,and pollutant concen— trations from watersheds in the Twin Cities metropolitan area, cling.Contml mid Management on Nowpoint Source Pollution[J]. Environmental Science,1998,19(5):87—96.(in Chinese)) Minnesota,USA[J].Water Research,2002,36(7):1743—1757. [14]Vanoni V A(1975).Sedimentation Engineering,ASCE Manual and Report on Engineerig nPractice,No.54,New York,NY. [33张展羽,卢敏,朱成立.城镇河道的水生态环境建设和保护[J]. 灌溉排水学报,2004,23(6):18—21.(zHANG Zhan-yu,LU Min, ・—・卜”—卜”—・卜..-4-”—卜“+ (上接第3O页) [151张光辉,刘少玉,张翠云,等.黑河流域地下水循环演化规律研 究[J].中国地质,2004,31(3):289—293.(zHANG Guang-hui, I.IU Shao-yu,ZHANG Cui—yun,et a1.Evolution of Groundwa— sheng,QIN Da-yong.A Study on Channel Leakage in the Low~ er Heihe River[J].Journal of Natural Resources,2005,20(2): 195~199.(in Chinese)) ter Circulation in the Heihe River Drainage Area[J].Chinese Geology,2004,20(1):63—66.(in Chinese)) [17]中国科学院计算机网络信息中心科学数据中 6[DB]:http:// datamirror.esdb.cn/landsat les/gap.jsp(International sci— entilte data service platforrn:http://datamirror.csdb.cn/land— sat/files/gap.jsp.(in Chinese)) [16]李云玲,裴源生,秦大庸.黑河下游河道渗漏规律研究EJ].自 然资源学报,2005,20(2):195—199.(LI Yun—ling,PEI Yuan一 ・ 8. 釜H 乌