韩璐遥, 叶兴成, 王 飞, 朱程亮, 栾 慕, 周 宏, 刘 俊
(1.河海大学 水文水资源学院, 江苏 南京 210098; 2.宿迁市水务局, 江苏 宿迁 223800)
近年来,基于源头控制理念的低影响开发措施(LID)被广泛应用于新老城区的改造中,不同的低影响开发措施减控效果主要通过SWMM模型来模拟分析[1-2]。SWMM模型(雨水管理模型)是由美国环保局在1971年首次提出,经过不断的发展,在最新版本SWMM5.1中提供了8种LID措施的控制模拟[3]。国内关于LID措施的研究多侧重于不同LID措施在不同降雨重现期下水量水质的削减控制效果对比,而对于各LID在不同排水条件下的综合削减控制效果没有过多关注。在应对区域尺度的海绵城市规划时,为提高源头控制整体效果,应针对不同的排水条件寻找适用于各片区的最佳LID布设方案。
基于此目标,本文研究了龙岩市中心城区不同排水片区的各LID措施布设效果差异,同时由于研究区域属于山区城市,各排水片区多依据地势就近排入相应河道,在选取LID布设方案时为最大程度缓解相应河道排涝压力,本文对各LID措施的综合削减控制效果分别进行长、短历时降水模拟。模拟结果包括分析短历时降水情况下LID措施对径流总量、洪峰流量及TSS负荷的削减效率,以及长历时降水情况下各LID措施对3个指标削减效果的降水历时敏感性。在此基础上为各片区选取了最佳LID布设方案,可为宏观尺度的优化LID布设方案提供参考。
龙岩市位于福建省西南部,属亚热带海洋性季风气候。该区域雨季长、雨量多、强度大,多年平均降雨量1 450~2 200 mm。中心城区地处龙岩市盆地,境内山岭与河谷相间,整体呈带状分布。地势总体呈南高北低、西高东低的走向。该区域依据地势进行排水分区,各分区排水可根据管网布局及地形因素就近入河。研究区域总面积134.56 km2,平均不透水面积39.3%。城市土地利用类型主要包括住宅(33.94%)、工业用地(21.74%)、道路(16.09%)、绿地(9.56%)、公共设施(9.43%)以及商业和服务设施(4.74%)。由于该区域管网布设标准偏低,且城市化进程较快,土地利用紧张,故雨水管网改扩建困难,暴雨洪水问题突出[4]。
中心城区包括城区、郊区及山丘区。本文通过GIS软件提取城区的地形地貌及土地利用布局,结合城区实际的管网排水情况来划分子汇水区。郊区位于中心城区外侧,未进行雨水管网布设,故依据地形资料进行分区。中心城区上游的山丘区采用《福建省暴雨洪水查算图表》进行产汇流计算,并通过汛期水库调度模式计算出流,将山丘区汇水概化为SWMM模型的上边界入流。根据研究区的管网布局和地形因素,地面径流可就近入河。最终将研究区划分为553个子汇水区,包括711个节点,698条管段和11个堰。模型概化见图1。
依据研究区实际情况并结合相关文献和SWMM用户手册,对模型参数进行选取。本文主要采用霍顿下渗模型来模拟研究区域的产流入渗过程。最大下渗率范围为25.4~50.8 mm/h,稳定下渗率为2.54 mm/h,衰减系数2~7 h-1。子流域坡度依据龙岩市中心城区实际地形及规划资料选取,范围为0.32%~10.01%。不透水区糙率取0.015,洼蓄量1 mm;透水区糙率取0.2,洼蓄量取5 mm。不透水区无洼蓄百分比在整个研究区范围为25%~80%。汇流模块采用非线性水库模型,排水系统的流量演算采用动力波演算模型[5]。
图1 龙岩市中心城区模型概化图
为分析各LID措施对面源污染减控效果的空间差异性,选取TSS作为污染指标[6-7],根据相关文献及研究区实际情况,在SWMM模型水质模块中率定并输入相应参数,如表1所示:
表1 水质模块模型参数表
SWMM模型参数的率定主要是对地表粗糙系数、洼蓄量等具体数值需要通过模型率定确定或者经多项调查研究加经验取值才能确定的参数进行率定。为检验基于SWMM模型构建的暴雨内涝模型是否符合城市内涝研究的要求,此次参数率定采取2010年6月14日22时至6月15日22时共计24 h实测暴雨,降雨总量为120.50 mm,主要用于积水点率定。为分析模拟的可靠性,模拟得出易涝区与福建省龙岩市城乡规划局提供的实际易涝区的对比情况见表2。
表2 模拟易涝区与实际易涝区对比表
从模拟结果来看,模拟易涝区与实际易涝区相比总体情况比较吻合,说明模型的计算结果较为可靠。但有些与实际情况仍存在差别,其原因主要有两方面:(1)由于地面淹水情况是调查资料,而实际
淹水调查资料不够详实,有些是估计得出,部分地区存在漏报错报现象,所以模拟易涝点比实际易涝点多;(2)由于系统无法模拟一些细小管道,对于管道临时性堵塞的情况也无法全面调查,所以对于一些局部地区的淹水点,无法可靠地模拟和分析。
为测试LID措施的效果并分析其对水量、水质在不同排水片区的减控作用,从上述概化区域中选取3个典型的独立排水片区。根据河流水系的分布、地形地貌的变化,结合现状排水系统的建设等情况,按原《龙岩市中心城区管线综合规划》,选取位于中心城区的上游、中游和下游的排水片区,记为片区A、片区B和片区C(见图1)。
3个排水片区均具有多种用地类型,其下垫面组成情况分别为片区A:屋顶38%、广场0%、道路14%、绿地48%;片区B:屋顶40%、广场0%、道路17%、绿地43%;片区C:屋顶33%、广场0%、道路26%、绿地41%。
3个排水片区占地面积和平均不透水率相差不大,但各自用地功能定位不同。通过GIS提取信息及实际调研情况,对比得出各排水片区的排水条件差异性见表3。
表3 排水条件差异性对比表
2.5.1 短历时设计暴雨 短历时降雨强度计算采用福建省住房和城乡建设厅发布的暴雨强度公式:
(1)
式中:q为降雨强度,L/(s·hm2);T为降雨重现期,a;t为降雨历时,min。
按照管道排水设计要求,短历时设计暴雨采用120min作为统计历时,并采用芝加哥雨型作为短历时设计雨型[8]。通过统计分析雨量站历时120min最大雨量时间分布,确定综合雨峰位置系数r=0.54。结合龙岩市暴雨强度公式,得出1年一遇、2年一遇及5年一遇降雨重现期下短历时设计降雨过程见图2。
图2 短历时设计降雨过程
2.5.2 长历时设计暴雨 为研究各LID措施对长历时降水的敏感性,对长历时设计暴雨采用同频率缩放法计算24h设计暴雨[9]。长历时设计暴雨时程分配按照福建省最大3d设计雨型表得到24h设计降雨过程,见图3。
图3 长历时设计降雨过程
根据研究区域特性选取适用于研究区域的低影响开发措施[10],包括雨水花园、绿色屋顶和渗透铺装。在本次制定方案时,除需结合研究区域实际情况选择布设LID措施及相应布设比例外,还要控制各工程投资总数的大致相同。各LID措施投资成本及维护状况如表4所示。综合考虑以上因素,在A、B、C3个排水片区拟定4种方案:(1)设置10%的雨水花园;(2)设置20%的绿色屋顶;(3)设置15%的渗透铺装;(4)组合措施(包括3%的雨水花园、7%的绿色屋顶和5%的渗透铺装)。为兼顾缓解管道排水与相应河道排涝的要求,选取2h短历时设计暴雨(包括重现期为1年一遇、2年一遇和5年一遇)和24h长历时设计暴雨(包括重现期为10年一遇、20年一遇和50年一遇)共6场,来模拟各LID模块在不同排水片区的径流总量、洪峰流量和污染物负荷的削减效果,以寻求在不同重现期降雨情况下的最佳总体效果。本次研究的总体效果评价包括在短历时降水情况下的减控效率,以及对长历时降雨的敏感性,即选取效率较高且较为稳定的LID措施方案。短历时设计雨量范围为44.1~66.7mm,长历时设计雨量范围为128.6~153.9mm。
表4 低影响开发技术设施投资及维护状况
径流总量反映了区域产流能力。为降低区域洪涝风险,需综合考虑LID措施对区域径流总量及洪峰流量的控制[11-12]。相对于无LID布置,不同LID布置方案下各片区径流总量削减率如表5所示。
表5 各LID布置方案下的径流总量削减率
在不同降雨重现期下,片区A、片区B、片区C的径流总量削减范围依次为3.85%~36.57%、4.60%~35.34%、2.83%~33.73%(表5)。综合长历时和短历时共6个重现期的降雨场次模拟结果,比较各方案对径流总量总体的削减效果(包括短历时降雨情况下的削减效率及长历时降雨情况下的敏感性)可以得出在片区A中以绿色屋顶方案最优,为9.84%~36.57%;片区B中组合措施最优,为10.97%~35.34%;片区C中雨水花园方案最优,为15.10%~29.87%。
随着降雨重现期的增大,各LID布设方案对径流削减控制效果均有所降低,但变化程度各有不同(表5)。其中以渗透铺装较为明显,即渗透铺装在短历时低重现期时对径流削减效果较好,在长历时高重现期时效果较差[14]。绿色屋顶及组合措施对不同降雨历时的敏感性一般,雨水花园对不同降雨历时的敏感性不明显,故在结合长历时降水时,LID措施对径流总量的削减效果以雨水花园最稳定。
相对于无LID方案,不同LID布置方案下各片区的峰值流量削减率如表6所示。
表6 各LID布置方案下的洪峰流量削减率
在不同降雨重现期下,片区A、片区B、片区C的洪峰流量削减范围依次为4.35%~29.76%,3.94%~26.03%,4.16%~24.84%(表6)。各LID布设方案对洪峰流量的总体削减效果,在片区A中绿色屋顶最好,为4.35%~29.76%;片区B中组合措施最好,为7.58~25.68%;片区C中渗透铺装最好,为10.37%~22.56%。
从表6中还可以看出,绿色屋顶在短历时低重现期时对洪峰流量削减控制效果较好;雨水花园和组合措施对洪峰流量的削减效果受降雨历时的影响一般;渗透铺装对于降雨历时的敏感性不明显,在不同降雨历时及不同重现期下对洪峰流量的削减效果较稳定。各LID布设方案对洪峰流量的控制效果均随着降雨重现期的增加有不同程度的降低。
本文主要研究了各LID布设方案在不同排水片区对TSS的控制效果的差异。相对于无LID布置,不同LID布置方案下对各片区的TSS削减率如表7所示。
表7 各LID布置方案下的TSS削减率
在不同降雨重现期下,片区A、片区B、片区C的TSS污染负荷削减范围依次为6.24%~44.86%,7.04%~42.25%,8.42%~45.11%(表7)。综合长、短历时共6个重现期的降雨场次分析,得出结论:片区A、片区B、片区C对TSS负荷的削减率均以雨水花园效果最好,分别为19.33%~38.37%、17.93%~39.29%和23.5%~40.89%。
在结合了长历时降雨的情况下,相比于其他布设方案,雨水花园对TSS负荷的削减效果对不同降雨历时及重现期的敏感程度最低(表7),因为降雨历时的增加延长了雨水花园与污染物的接触时间,从一定程度上减弱了高强度降雨对地表TSS冲刷的影响。绿色屋顶和渗透铺装在短历时低重现期下对TSS削减效果较好,但在长历时高重现期时由于受到雨水的冲刷,导致部分被吸附存储的污染物再次进入径流,从而降低了其对TSS的削减控制效果[15]。
为综合评价各LID措施对径流总量、洪峰流量及面源污染负荷的削减控制效果,根据海绵城市的指导理念并结合研究区的实际情况,设定径流总量削减率指标所占权重为0.6,洪峰流量和污染物负荷削减率指标所占权重均为0.2。
综合长、短降水历时情况下对径流总量、洪峰流量和TSS负荷3个指标的削减控制效果可以得出:
片区A的LID最佳布设方案为绿色屋顶方案。由于片区A以仓储用地及居住用地为主,结构相对简单,故划分的子汇水区数量少且面积大,各子汇水区之间影响较小。LID布设方案中以绿色屋顶方案布设比例最大,因此对于径流总量及洪峰流量的削减率较高。且片区A位于上游,污染物浓度相对较低,综合以上因素,片区A优选绿色屋顶方案。
片区B优选组合措施方案。片区B位于中心城区中游,以商业用地及住宅用地为主,人流量较大,用地布局复杂,故子汇水区数量多且面积相对较小,各子汇水区易受其他子汇水区来水的影响。组合措施中的绿色屋顶能够有效削减各子汇水区自身的降雨径流,而雨水花园和渗透铺装具有一定的处理其他不渗透区域径流的能力,在短历时低重现期下,除了处理自身的降雨径流外,还可以处理部分其他来水从而提高整体效率。片区B污染物浓度相对较高,组合措施中各LID均具有削减TSS负荷的能力,且雨水花园在长历时下对污染物的削减仍可发挥显著作用。
片区C优选雨水花园方案。片区C处于中心城区下游,以工业用地和交通枢纽建设用地为主,该片区污染物浓度最高。在该片区,要求LID措施能够更有效减少TSS污染负荷。雨水花园对TSS负荷的削减控制效果受降雨历时的影响小于其他3个方案,故雨水花园方案效果较稳定。该片区的不渗透区域相对分散,而雨水花园能够处理其他不渗透区域来水的比例最高,且雨水花园对降水历时敏感性相对较小,因此对该片区的径流总量及洪峰流量也具有较好的削减效果。综合以上因素,片区C优选雨水花园方案。
(1)本文通过模拟各LID布设方案在不同排水片区,对于径流总量、洪峰流量及面源污染负荷(TSS)的削减控制效果,可以看出各LID措施在不同的排水片区表现相似,但效果不同。
(2)各LID措施对于不同排水片区水量和水质的削减率均随着降雨历时及重现期的增加而降低,但降低程度略有不同。渗透铺装对径流总量削减在短历时降雨下效果较好,绿色屋顶对于洪峰流量削减在短历时降雨情况下效果较好。对于TSS削减,渗透铺装和绿色屋顶均在短历时降雨下较为有效。
(3)通过权重分配对各LID措施效果进行综合评价,可以得出结论如下:上游排水片区A优选绿色屋顶方案;中游排水片区B优选组合措施方案;下游排水片区C优选雨水花园方案。
[1] 栾 慕,刘 俊.LID措施在老城区内涝风险管理中的应用[J].中国农村水利水电,2017(1):69-74.
[2] 官奕宏,吕 谋,王 灿,等.低影响开发技术的雨洪控制效果及水质影响分析-基于SWMM模型[J].中国农村水利水电,2017(1):84-87+93.[3]ROSSMANLA.Stormwatermanagementmodeluser'smanual(Version5.1)[M].WashingtonDC,US:EPA,2015.
[4] 李 朋.低影响开发措施净流控制效果及空间差异性研究[D].南京:河海大学,2016.
[5] 王文亮,李俊奇,车 伍,等.城市低影响开发雨水控制利用系统设计方法研究[J].中国给水排水,2014,30(24):12-17
[6] 张一平.城市化与城市水环境[J].城市环境与城市生态,1998,11(2):20-22+27.[7] 宫 莹,阮晓红,胡晓东.我国城市地表水环境非点源污染的研究进展[J].中国给水排水,2003,19(3):21-23.[8] 符 锐,罗龙洪,刘 俊,等.SWMM模型中的低影响开发模块在排水防涝系统中的应用[J].水电能源科学,2014,32(9):71-74.
[9] 王 静.基于SWMM模型的山地城市暴雨径流效应及生态化改造措施研究[D].重庆:重庆大学,2012.
[10] 李 霞,石宇亭,李国金.基于SWMM和低影响开发模式的老城区雨水控制模拟研究[J].给水排水,2015,51(5):152-156.
[11] 胡爱兵,任心欣,丁 年,等.基于SWMM的深圳市某区域LID设施布局与优化[J].中国给水排水,2015,31(21):96-100.[12] 常 静,刘 敏,许世远,等.上海城市降雨径流污染时空分布与初始冲刷效应[J].地理研究,2006,25(6):994-1002.[13] 任玉芬,王效科,韩 冰,等.城市不同下垫面的降雨径流污染[J].生态学报,2005,25(12):3225-3230.
[14] 罗红梅,车 伍,李俊奇,等.雨水花园在雨洪控制与利用中的应用[J].中国给水排水,2008,24(6):48-52.
[15] 赵剑强.城市地表径流污染与控制[M].北京:中国环境科学出版社,2002.