CN112417657B - 基于不同下垫面污染物贡献率的海绵城市优化设计方法 - Google Patents
基于不同下垫面污染物贡献率的海绵城市优化设计方法 Download PDFInfo
- Publication number
- CN112417657B CN112417657B CN202011247357.1A CN202011247357A CN112417657B CN 112417657 B CN112417657 B CN 112417657B CN 202011247357 A CN202011247357 A CN 202011247357A CN 112417657 B CN112417657 B CN 112417657B
- Authority
- CN
- China
- Prior art keywords
- research area
- target research
- point source
- underlying surfaces
- various
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000003344 environmental pollutant Substances 0.000 title claims abstract description 83
- 231100000719 pollutant Toxicity 0.000 title claims abstract description 83
- 238000013461 design Methods 0.000 title claims abstract description 49
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000005457 optimization Methods 0.000 title claims abstract description 28
- 238000011160 research Methods 0.000 claims abstract description 112
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 41
- 238000012544 monitoring process Methods 0.000 claims description 24
- 238000004088 simulation Methods 0.000 claims description 15
- 238000005070 sampling Methods 0.000 claims description 14
- 239000000126 substance Substances 0.000 claims description 12
- 239000000356 contaminant Substances 0.000 claims description 10
- 238000011161 development Methods 0.000 claims description 6
- 230000008595 infiltration Effects 0.000 claims description 6
- 238000001764 infiltration Methods 0.000 claims description 6
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 5
- 229910052760 oxygen Inorganic materials 0.000 claims description 5
- 239000001301 oxygen Substances 0.000 claims description 5
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 4
- XKMRRTOUMJRJIA-UHFFFAOYSA-N ammonia nh3 Chemical compound N.N XKMRRTOUMJRJIA-UHFFFAOYSA-N 0.000 claims description 4
- 230000001186 cumulative effect Effects 0.000 claims description 4
- 229910052698 phosphorus Inorganic materials 0.000 claims description 4
- 239000011574 phosphorus Substances 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 229920006395 saturated elastomer Polymers 0.000 claims 1
- 238000010276 construction Methods 0.000 abstract description 18
- 230000000694 effects Effects 0.000 abstract description 5
- 238000004364 calculation method Methods 0.000 description 11
- 238000009825 accumulation Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000009991 scouring Methods 0.000 description 6
- 230000009466 transformation Effects 0.000 description 4
- 238000007726 management method Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000000746 purification Methods 0.000 description 2
- KHOITXIGCFIULA-UHFFFAOYSA-N Alophen Chemical compound C1=CC(OC(=O)C)=CC=C1C(C=1N=CC=CC=1)C1=CC=C(OC(C)=O)C=C1 KHOITXIGCFIULA-UHFFFAOYSA-N 0.000 description 1
- 244000025254 Cannabis sativa Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 239000010426 asphalt Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000007716 flux method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000012847 principal component analysis method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000005211 surface analysis Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Development Economics (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了基于不同下垫面污染物贡献率的海绵城市优化设计方法,包括:建立目标研究区的基础资料数据库;根据基础资料数据库构建目标研究区的SWMM模型;基于构建的SWMM模型,模拟计算目标研究区内不同下垫面的各类非点源污染物负荷贡献率,确定各类下垫面的非点源污染特征;对目标研究区进行海绵城市优化设计。本发明基于SWMM模型对目标研究区的不同下垫面的径流污染物进行模拟,获得各类下垫面的非点源污染物负荷贡献率,根据贡献率针对性的进行海绵城市建设优化设计,达到提高海绵城市改造项目的控制径流污染物源头的效果。
Description
技术领域
本发明属于海绵城市建设技术领域,具体涉及基于不同下垫面污染物贡献率的海绵城市优化设计方法。
背景技术
近年来,随着城市化建设不断推进,导致下垫面发生变化,尤其是路面和屋面等不透水下垫面面积大幅增加,绿地面积比例减小,发生降雨时不仅易造成城市内涝问题,而且降雨径流会挟带大量污染物,诱发严重的城市面源污染(也称为非点源污染)问题。
研究表明在一场降雨事件中,不同下垫面对于不同污染物质的冲刷作用不同,下垫面主要分屋面、路面和绿地三种,由于屋面材料多为沥青,因此屋面COD(化学需氧量)和SS浓度(悬浮物浓度)较高;车辆行驶和道路材料的缘故,造成路面TN浓度(总氮浓度)和COD浓度较高;绿地SS浓度较高。所以下垫面的占比会直接影响城市的径流污染程度。现在的城市规划和海绵城市改造中都存在着不透水路面、屋面和绿化面积之间设置比例不合理,盲目开发建设的问题。通过掌握不同下垫面对污染物输送的贡献率,利用贡献率合理规划各下垫面的占比,可为城市规划和海绵城市改造提供有效的技术指导,所以急需一种计算快速,结果精准的污染物贡献率计算方法。
目前对于不同下垫面污染物贡献率的计算方法,主要有以下几种:一是通量法,收集监测点的断面污染物通量,通过物质通量公式计算贡献率;二是利用多元统计法计算污染物贡献率;三是利用主成分分析法分析污染物的贡献率。
然而以上方法存在的问题是:忽视了降雨冲刷的影响、计算结果误差大、不能准确地反映不同下垫面污染物的贡献率;因为降雨强度、降雨历时的不同,需要不断的更新监测数据,导致工作量大。
发明内容
本发明实施方式的目的在于提供基于不同下垫面污染物贡献率的海绵城市优化设计方法,以克服上述技术缺陷。
为解决上述技术问题,本发明提供了基于不同下垫面污染物贡献率的海绵城市优化设计方法,包括:
建立目标研究区的基础资料数据库;
根据基础资料数据库构建目标研究区的SWMM模型;
基于构建的SWMM模型,获取目标研究区内不同下垫面的各类非点源污染物负荷贡献率,确定各类下垫面的非点源污染特征;
对目标研究区进行海绵城市优化设计。
进一步地,建立目标研究区的基础资料数据库,具体包括:
获取目标研究区的平面设计图和卫星影像图;
根据平面设计图和卫星影像图,对目标研究区内的不同下垫面类型进行划分。
进一步地,基于构建的模型,获取目标研究区内不同下垫面的各类非点源污染物负荷贡献率,确定各类下垫面的各类非点源污染物特征之前,还包括:
获取所述目标研究区不同下垫面的实测流量数据和实测水质数据;
根据实测流量数据对SWMM模型的水文参数进行率定;
根据实测水质数据对SWMM模型的水质参数进行率定。
优选地,获取目标研究区的实测流量数据,具体包括:
监测所述目标研究区的各类下垫面的地表径流流量和关键排水管段流量;
同时对所述目标研究区的实际降雨场次进行监测;
其中对所述关键排水管段进行连续监测,监测频率≥15min/次,监测的降雨场数不少于三场。
进一步地,SWMM模型的水文参数包括管道糙率、透水区曼宁系数、不透水区曼宁系数、最大下渗系数、稳定下渗系数、透水区洼蓄量以及不透水区洼蓄量。
优选地,获取目标研究区的实测水质数据,具体包括:
在所述目标研究区对已划分的各类下垫面设置地表径流采样点;
同时收集关键排水管段流量监测点的水样;
利用自动采样器或手动收集方式,在所述地表径流采样点和关键管段流量监测点以预设的时间点取样,水样数量≥6,对每个水样分别检测悬浮物浓度、化学需氧量浓度、氨氮浓度以及总磷浓度。
进一步地,SWMM模型的水质参数包括最大累积量、半饱和累积时间、冲刷系数以及冲刷指数。
进一步地,基于构建的SWMM模型,获取目标研究区内不同下垫面的各类非点源污染物负荷贡献率,确定各类下垫面的各类非点源污染物特征,具体包括:
S1,基于目标研究区建立SWMM模型,并通过收集的不同下垫面的场次雨量、流量及水质数据对目标研究区总体及各类下垫面的水文水质参数进行率定及验证;
S2,基于已建立的SWMM模型,依次利用目标研究区的设计降雨条件进行降雨模拟,获得在不同重现期的设计降雨条件下目标研究区的所有下垫面的单类非点源污染物的总负荷P 0 ;
S3,依次取消所有下垫面的其中同一类非点源污染物设置参数,在相同设计降雨条件下重新进行降雨模拟,获得各类下垫面的其他类非点源污染物的总负荷P i ;
S4,按照以下公式获取目标研究区内不同下垫面的各类非点源污染物负荷贡献率:
式中:a i 为第i类下垫面的单类非点源污染物负荷贡献率,P 0 为目标研究区的单类非点源污染物的总负荷,P i 为第i类下垫面的其他类非点源污染物的总负荷;
S5,根据各类非点源污染物负荷贡献率,确定各类下垫面的各类非点源污染物特征,具体包括:
确定目标研究区不同下垫面的各类污染物的浓度变化规律及不同重现期下目标研究区的不同下垫面的各类污染负荷的排放量。
进一步地,对所述目标研究区进行海绵城市优化设计,具体包括:
在年径流总量控制率的条件约束下,针对不同下垫面的污染负荷贡献率针对性地设置海绵城市开发措施。
优选地,基础资料数据库包括目标研究区平面图、DEM高程图、管网管线以及降雨资料。
本发明的有益效果如下:
本发明基于SWMM模型对目标研究区的各类下垫面污染物模拟,获得各类下垫面非点源污染物负荷贡献率,根据贡献率针对性的优化布设海绵城市控制措施,达到提高海绵城市改造项目径流污染物源头控制效果的目的。
为让本发明的上述内容能更明显易懂,下文特举优选实施例,并结合附图,作详细说明如下。
附图说明
图1是基于不同下垫面污染物贡献率的海绵城市优化设计方法的流程图。
具体实施方式
以下由特定的具体实施例说明本发明的实施方式,本领域技术人员可由本说明书所揭示的内容轻易地了解本发明的其他优点及功效。
需说明的是,在本发明中,图中的上、下、左、右即视为本说明书中所述的基于不同下垫面污染物贡献率的海绵城市优化设计方法的上、下、左、右。
现参考附图介绍本发明的示例性实施方式,然而,本发明可以用许多不同的形式来实施,并且不局限于此处描述的实施例,提供这些实施例是为了详尽地且完全地公开本发明,并且向所属技术领域的技术人员充分传达本发明的范围。对于表示在附图中的示例性实施方式中的术语并不是对本发明的限定。在附图中,相同的单元/元件使用相同的附图标记。
除非另有说明,此处使用的术语(包括科技术语)对所属技术领域的技术人员具有通常的理解含义。另外,可以理解的是,以通常使用的词典限定的术语,应当被理解为与其相关领域的语境具有一致的含义,而不应该被理解为理想化的或过于正式的意义。
第一实施方式:
本实施方式涉及基于不同下垫面污染物贡献率的海绵城市优化设计方法,参见图1,包括:
建立目标研究区的基础资料数据库;
根据所述基础资料数据库构建所述目标研究区的SWMM模型;
基于构建的SWMM模型,获取所述目标研究区内不同下垫面的各类非点源污染物负荷贡献率,确定各类下垫面的非点源污染特征;
对所述目标研究区进行海绵城市优化设计。
其中,可以在ArcGIS软件内建立目标研究区的基础资料数据库。
建设海绵城市,即构建低影响开发雨水系统,主要是指通过“渗、滞、蓄、净、用、排”等多种技术途径,实现城市良性水文循环,提高对径流雨水的渗透、调蓄、净化、利用和排放能力,维持或恢复城市的海绵功能,目前的海绵城市建设往往根据项目地形进行设计,忽略了各类下垫面的非点源污染负荷排放贡献率,导致海绵城市对城市非点源污染物去除率效率不高,而目前的城市非点源污染情况日益严重,必须对于区域非点源污染情况进行更为深入的了解,以采取有效的污染控制措施。
针对以上问题,本发明采取SWMM模型作为工具模拟目标研究区的降雨条件,需要说明的是,SWMM(storm water management model,暴雨洪水管理模型)是由美国环保部(EPA)开发用于模拟动态降水-径流过程的一种软件,在国内外降雨径流管理方面都得到了广泛应用。
本发明利用SWMM模型计算不同下垫面的污染物贡献率,具体为:可以跟踪模拟不同降雨强度下、不同降雨历时下目标研究区下垫面的污染物贡献率,灵活性高;SWMM模型包含了径流模块、输送模块、扩展输送模块和贮存处理模块等4个计算模块,计算结果准确度高;为城市规划和海绵城市改造提供可靠的技术指导,使不同下垫面的比例规划更加合理;模型建立成功后,输入降雨重现期,即可快速计算出不同下垫面的污染物贡献率,操作简单,易于推广使用。
基础资料数据库:包括目标研究区平面图、DEM高程图、管网管线以及降雨资料,但不仅限于此,可以根据模拟目的自行调整所需数据,此外以上数据的获取都是本领域的常规技术,可查表获得或监测获得等等,在此不作详细说明。
在构建SWMM模型时,将取得的基础资料数据按照需求输入SWMM模型软件,然后通过收集的不同下垫面的场次雨量、流量及水质数据对目标研究区总体及各类下垫面的水文水质参数进行率定及验证。
基于建立的SWMM模型,通过当地暴雨强度公式下的不同重现期设计降雨对目标研究区模型进行降雨模拟,计算目标研究区不同降雨重现期情况下的各类污染物负荷总量,再进一步计算目标研究区内各类下垫面的各类非点源污染物负荷贡献率,确定各类下垫面的非点源污染特征,针对性地布设目标研究区的海绵城市建设方案,为城市非点源污染控制和海绵城市建设提供科学依据。
第二实施方式:
本实施方式涉及基于不同下垫面污染物贡献率的海绵城市优化设计方法,包括:
建立目标研究区的基础资料数据库;
获取目标研究区的平面设计图和卫星影像图;
根据平面设计图和卫星影像图,对目标研究区内的不同下垫面类型进行划分;
根据所述基础资料数据库构建所述目标研究区的SWMM模型;
获取所述目标研究区不同下垫面的实测流量数据和实测水质数据;
根据实测流量数据对SWMM模型的水文参数进行率定;
根据实测水质数据对SWMM模型的水质参数进行率定;
基于构建的SWMM模型,获取所述目标研究区内不同下垫面的各类非点源污染物负荷贡献率,确定各类下垫面的非点源污染特征;
对所述目标研究区进行海绵城市优化设计。
其中对目标研究区内的不同下垫面类型进行划分,可以将目标研究区内的下垫面划分为四类,分别是屋顶、路面、绿地以及未开发土地;
相比第一实施方式,本实施方式有两处不同:
(1)对目标研究区的下垫面进行了分类,根据目标研究区的平面图纸以及卫星影像图,将下垫面划分为屋顶、路面、绿地、及未开发土地四类,通过ArcGIS地理信息软件(计算机制图应用)对不同下垫面进行分割和分类,形成不同矢量图层数据,作为水力模型分析城市降雨产流和汇流机制的基础,此外在雨水分区划分中,需对各下垫面采样点地块做单独划分。
(2)在构建了SWMM模型后,进行参数率定,具体方法如下:
对目标研究区不同下垫面的实际降雨场次径流过程和污染过程及总排口流量及水质监测,进行模型参数的率定与检验,即依据取得的实测数据,运用手动试错法对模型参数进行率定和检验,其中主要率定参数分为水文参数和水质参数,水文参数包括管道糙率、透水区曼宁系数、不透水区曼宁系数、最大下渗系数、稳定下渗系数、透水区洼蓄量以及不透水区洼蓄量等敏感参数取值范围;水质参数包括最大累积量、半饱和累积时间、冲刷系数以及冲刷指数。
需要说明的是,参数率定的实质就是先假定一组参数,代入模型得到计算结果,然后把计算结果与实测数据进行比较,若计算值与实测值相差不大,则把此时的参数作为模型的参数;若计算值与实测值相差较大,则调整参数代入模型重新计算,再进行比较,直到计算值与实测值的误差满足一定的范围。
在本申请中,率定参数是模型水文参数和模型水质参数。
实测流量数据具体包括:
通过对所述目标研究区的各类下垫面的地表径流和关键排水管段的流量进行监测;
同时对所述目标研究区的实际降雨场次进行监测;
对管网关键节点均进行连续监测,监测频率≥15min/次,监测的降雨场数不少于三场;
在实测流量数据时,可以结合目标研究区的相关设计图纸与现场勘测,对目标研究区内雨水设施的衔接关系、汇水面积,判断设施的设计降雨量标准及对应的年径流总量控制率。
实测水质参数具体包括:
在所述目标研究区对已划分的各类下垫面设置地表径流采样点;
同时收集关键排水管段流量监测点的水样;
利用自动采样器或手动收集方式,在所述地表径流采样点和关键管段流量监测点以预设的时间点取样,水样数量≥6,对每个水样分别检测悬浮物浓度、化学需氧量浓度、氨氮浓度以及总磷浓度。
为了清楚地表述实测水质参数的获取过程,以下将以具体数值举例说明:
针对不同下垫面及关键节点布设地表径流采样点,通过自动采样器或手动收集方式,在出现径流后以0min、5min、10min、20min、30min、60min、90min、120min的时间间隔取样,每个样品体积约500ml;若降雨历时较长,后期采样数量可适当增加、取样间隔可适当延长;若降雨历时较短,可适当减少采样数量,但一般不少于6个水样,水质检测项目包括悬浮物(SS)、化学需氧量(COD)、氨氮(NH3-N)、总磷(TP)等污染物指标。
建模过程:
在子汇水区划分及管网概化过程中:基于所收集的数据对目标研究区进行下垫面分类,运用ArcGIS软件划分汇水区并计算汇水区面积、特征宽度、坡度、不透水面百分比及不同下垫面面积等参数;并将研究区管网、节点内底标高、管段长度等数据导入ArcGIS软件中,输出为inp文件;在SWMM软件中,依据手册及设计资料将其他参数设置完成后,设置对应的SS、COD、NH3-N、TN、TP等主要污染因子,依据各下垫面监测各类污染物数据,输入地表物累积与冲刷所需的参数,地表物累积参数包括最大累积量和累积时间,冲刷参数包括冲刷系数和冲刷指数。
本实施方式运用SWMM模型对目标研究区地表径流进行模拟评估,结合现有的资料和下垫面解析结果,将ArcGIS计算的各类雨水分区面积和特征宽度输入SWMM模型;模型的管道传输演算模型选取运动波模型;渗入模型选取Horton模型;管道数据按照研究地区管线资料进行设置;污染物模型选用饱和函数来作为污染物非雨期和雨期的累积过程。
第三实施方式:
本实施方式涉及基于不同下垫面污染物贡献率的海绵城市优化设计方法,参见图1,包括:
建立目标研究区的基础资料数据库;
根据所述基础资料数据库构建所述目标研究区的SWMM模型;
基于构建的SWMM模型,获取所述目标研究区内不同下垫面的各类非点源污染物负荷贡献率,确定各类下垫面的非点源污染特征;
对所述目标研究区进行海绵城市优化设计。
其中,基于构建的SWMM模型,获取目标研究区内不同下垫面的各类非点源污染物负荷贡献率,确定各类下垫面的各类非点源污染物特征,具体包括:
S1,基于目标研究区建立SWMM模型,并通过收集的不同下垫面的场次雨量、流量及水质数据对目标研究区总体及各类下垫面的水文水质参数进行率定及验证;
S2,基于已建立的SWMM模型,依次利用目标研究区的设计降雨条件进行降雨模拟,获得在不同重现期的设计降雨条件下目标研究区的所有下垫面的单类非点源污染物的总负荷P 0 ;
S3,依次取消所有下垫面的其中同一类非点源污染物设置参数,在相同设计降雨条件下重新进行降雨模拟,获得各类下垫面的其他类非点源污染物的总负荷P i ;
S4,按照以下公式获取目标研究区内不同下垫面的各类非点源污染物负荷贡献率:
式中:a i 为第i类下垫面的单类非点源污染物负荷贡献率,P 0 为目标研究区的单类非点源污染物的总负荷,P i 为第i类下垫面的其他类非点源污染物的总负荷;
S5,根据各类非点源污染物负荷贡献率,确定各类下垫面的各类非点源污染物特征,具体包括:
确定目标研究区不同下垫面的各类污染物的浓度变化规律及不同重现期下目标研究区的不同下垫面的各类污染负荷的排放量。
在上述步骤中,基于已建立好的模型,依次运用当地典型雨型对目标研究区进行降雨模拟,计算在典型降雨状况下目标研究区的各类污染物负荷总量;然后依次取消各类下垫面的污染物设置参数,在相同降雨重现期条件下重新对模型进行模拟,计算各类污染物负荷总量,前后差值作为某类下垫面对目标研究区总污染负荷的贡献。
而布设目标研究区的海绵城市建设方案,具体包括:
在年径流总量控制率的条件约束下,针对不同下垫面的污染负荷贡献率针对性地设置海绵城市开发措施。
例如道路污染物贡献率较高,在建立好的SWMM模型中针对性增加道路海绵城市措施的数量,在相同的设计重现期下进行模拟,通过模拟结果判断所设立的海绵城市措施是否有效达到有效控制道路非点源污染的目的。
通过上述计算方法得出的各类下垫面各类污染物贡献率,可提供直观的目标研究区域各类下垫面径流的非点源污染特征,通过对不同下垫面各类污染物贡献率的分析,并结合当地实际状况可针对性的设计相应的LID措施,更加直接有效的缓解当地降雨径流的非点源污染状况。
本发明通过对研究区地表径流的准确监测及数据模拟,通过计算不同下垫面中各类污染物对研究区整体降雨径流的污染物贡献率,确定该地区的非点源污染特征,来针对性设计海绵城市建设方案。对城市建设而言,传统的海绵城市建设往往过于强调用径流控制率作为海绵城市建设的依据,而忽略了研究区内降雨径流的非点源污染对城市整体水环境的威胁,因而导致现行的海绵城市建设方案对降雨径流非点源污染的处理效果不好。而目前的城市地区非点源污染情况日益严重,要求对于区域非点源污染情况进行更为深入的了解,以采取有效的污染控制措施。以研究区域不同下垫面污染物贡献率的海绵城市建设方法可以较好地满足区域非点源污染控制的需要。本发明基于城市区域非点源污染控制的现实需求和非点源污染模拟研究发展的需要,开展的对于城市海绵城市建设的研究具有多重意义。模拟的结果可为城市非点源污染控制和海绵城市建设提供科学依据。
第四实施方式:
某研究区根据本发明提供的基于不同下垫面污染物贡献率的海绵城市优化设计方法,进行了以下实施:
以西安市某小区为研究对象,基于该研究区用地类型,将该研究区下垫面划分为屋面、道路、绿地三类;
对不同下垫面监测点的流量水质数据进行监测收集,并通过监测的场次降雨数据对模型进行率定验证;
结果表明,模型的参数设定在误差范围内,可以用于建立模型。
利用西安市暴雨强度公式,设定降雨重现期分别为P=0.5a、P=1a、P=3a和P=5a,其中p是重现期,a是年;
以上述不同设计重现期降雨对已建立的SWMM模型进行模拟,通过计算最终得到不同重现期下的各类下垫面污染贡献率;
结果表明:
研究区屋面的污染贡献率最大,其次为道路和绿地,屋面道路和绿地对研究区的SS平均污染贡献率分别为62.74%、32.47%和4.79%;COD的平均污染贡献率分别为52.56%、39.82%、7.62%;NH3-N的平均污染贡献率分别为65.95%、21.78%、12.27%;TP的平均污染贡献率分别为59.58%、28.75%、11.67%;
通过模拟结果,针对性的增加研究区绿色屋顶及雨水桶等措施,用以有效削减屋顶污染贡献率,对于道路可以增设透水铺装和植草沟等措施用以有效削减道路污染贡献率。
本领域的普通技术人员可以理解,上述各实施方式是实现本发明的具体实施例,而在实际应用中,可以在形式上和细节上对其作各种改变,而不偏离本发明的精神和范围。
Claims (8)
1.基于不同下垫面污染物贡献率的海绵城市优化设计方法,其特征在于,包括:
步骤S1,建立目标研究区的基础资料数据库;
步骤S2,根据所述基础资料数据库构建所述目标研究区的SWMM模型;
步骤S3,基于构建的SWMM模型,获取所述目标研究区内不同下垫面的各类非点源污染物负荷贡献率,确定各类下垫面的非点源污染特征,具体包括:
S301,基于目标研究区建立SWMM模型,并通过收集的不同下垫面的场次雨量、流量及水质数据对目标研究区总体及各类下垫面的水文水质参数进行率定及验证;
S302,基于已建立的SWMM模型,依次利用所述目标研究区的设计降雨条件进行降雨模拟,获得在不同重现期的设计降雨条件下所述目标研究区的所有下垫面的单类非点源污染物的总负荷P 0;
S303,依次取消所述所有下垫面的其中同一类非点源污染物设置参数,在相同设计降雨条件下重新进行降雨模拟,获得所述各类下垫面的其他类非点源污染物的总负荷P i;
S304,按照以下公式获取目标研究区内不同下垫面的各类非点源污染物负荷贡献率:
式中:a i为第i类下垫面的单类非点源污染物负荷贡献率,P 0为目标研究区的单类非点源污染物的总负荷,P i为第i类下垫面的其他类非点源污染物的总负荷;
S305,根据各类非点源污染物负荷贡献率,确定各类下垫面的各类非点源污染物特征,具体包括:
确定目标研究区不同下垫面的各类污染物的浓度变化规律及不同重现期下目标研究区的不同下垫面的各类污染负荷的排放量;
步骤S4,对所述目标研究区进行海绵城市优化设计,具体包括:
在年径流总量控制率的条件约束下,针对不同下垫面的污染负荷贡献率针对性地设置海绵城市开发措施。
2.如权利要求1所述的基于不同下垫面污染物贡献率的海绵城市优化设计方法,其特征在于,建立目标研究区的基础资料数据库,具体包括:
获取所述目标研究区的平面设计图和卫星影像图;
根据所述平面设计图和卫星影像图,对所述目标研究区内的不同下垫面类型进行划分。
3.如权利要求1或2所述的基于不同下垫面污染物贡献率的海绵城市优化设计方法,其特征在于,基于构建的SWMM模型,获取所述目标研究区内各类下垫面的各类非点源污染物负荷贡献率,确定各类下垫面的各类非点源污染物特征之前,还包括:
获取所述目标研究区不同下垫面的实测流量数据和实测水质数据;
根据所述实测流量数据对所述SWMM模型的水文参数进行率定;
根据所述实测水质数据对所述SWMM模型的水质参数进行率定。
4.如权利要求3所述的基于不同下垫面污染物贡献率的海绵城市优化设计方法,其特征在于,获取目标研究区的实测流量数据,具体包括:
监测所述目标研究区的各类下垫面的地表径流流量和关键排水管段流量;
同时对所述目标研究区的实际降雨场次进行监测;
其中对所述关键排水管段进行连续监测,监测频率≥15min/次,监测的降雨场数不少于三场。
5.如权利要求4所述的基于不同下垫面污染物贡献率的海绵城市优化设计方法,其特征在于,所述SWMM模型的水文参数包括管道糙率、透水区曼宁系数、不透水区曼宁系数、最大下渗系数、稳定下渗系数、透水区洼蓄量以及不透水区洼蓄量。
6.如权利要求3所述的基于不同下垫面污染物贡献率的海绵城市优化设计方法,其特征在于,获取目标研究区的实测水质数据,具体包括:
在所述目标研究区对已划分的各类下垫面设置地表径流采样点;
同时收集关键排水管段流量监测点的水样;
利用自动采样器或手动收集方式,在所述地表径流采样点和关键管段流量监测点以预设的时间点取样,水样数量≥6,对每个水样分别检测悬浮物浓度、化学需氧量浓度、氨氮浓度以及总磷浓度。
7.如权利要求6所述的基于不同下垫面污染物贡献率的海绵城市优化设计方法,其特征在于,所述SWMM模型的水质参数包括最大累积量、半饱和累积时间、冲刷系数以及冲刷指数。
8.如权利要求1所述的基于不同下垫面污染物贡献率的海绵城市优化设计方法,其特征在于,所述基础资料数据库包括所述目标研究区平面图、DEM高程图、管网管线以及降雨资料。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011247357.1A CN112417657B (zh) | 2020-11-10 | 2020-11-10 | 基于不同下垫面污染物贡献率的海绵城市优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011247357.1A CN112417657B (zh) | 2020-11-10 | 2020-11-10 | 基于不同下垫面污染物贡献率的海绵城市优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112417657A CN112417657A (zh) | 2021-02-26 |
CN112417657B true CN112417657B (zh) | 2024-01-26 |
Family
ID=74780507
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011247357.1A Active CN112417657B (zh) | 2020-11-10 | 2020-11-10 | 基于不同下垫面污染物贡献率的海绵城市优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112417657B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113240260A (zh) * | 2021-05-07 | 2021-08-10 | 中冶南方城市建设工程技术有限公司 | 一种公共建筑海绵设施运行效果评价方法 |
CN113379299B (zh) * | 2021-06-28 | 2022-11-11 | 哈尔滨工业大学 | 一种基于模型的合流制溢流污染情况的评价方法 |
CN113657744B (zh) * | 2021-08-11 | 2024-04-12 | 厦门市城市规划设计研究院有限公司 | 一种海绵城市建设面源污染控制指标的确定方法 |
CN113723854B (zh) * | 2021-09-08 | 2023-10-31 | 北京市水科学技术研究院 | 一种合流制溢流的调控方法、装置及电子设备 |
CN113836710A (zh) * | 2021-09-15 | 2021-12-24 | 西安理工大学 | 一种排水管网外水入侵的快速追踪方法 |
CN117077566A (zh) * | 2023-08-24 | 2023-11-17 | 哈尔滨师范大学 | 基于swmm的典型流域非点源污染模拟与控制决策系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100828968B1 (ko) * | 2007-11-27 | 2008-06-10 | (주)웹솔루스 | 지리정보시스템과 연계된 하수관거 유지관리방법 및 그기능을 탑재한 하수관거 유지관리시스템 |
CN109948866A (zh) * | 2019-04-11 | 2019-06-28 | 沈阳建筑大学 | 一种基于swwm模拟的北方滨海地区海绵城市规划方法 |
CN110728062A (zh) * | 2019-10-16 | 2020-01-24 | 中建水务环保有限公司 | 一种基于swmm的农村非点源污染模拟方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050273300A1 (en) * | 2003-09-29 | 2005-12-08 | Patwardhan Avinash S | Method and system for water flow analysis |
-
2020
- 2020-11-10 CN CN202011247357.1A patent/CN112417657B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100828968B1 (ko) * | 2007-11-27 | 2008-06-10 | (주)웹솔루스 | 지리정보시스템과 연계된 하수관거 유지관리방법 및 그기능을 탑재한 하수관거 유지관리시스템 |
CN109948866A (zh) * | 2019-04-11 | 2019-06-28 | 沈阳建筑大学 | 一种基于swwm模拟的北方滨海地区海绵城市规划方法 |
CN110728062A (zh) * | 2019-10-16 | 2020-01-24 | 中建水务环保有限公司 | 一种基于swmm的农村非点源污染模拟方法 |
Non-Patent Citations (3)
Title |
---|
SWMM的校园区降雨径流污染及LID控制模拟;周明来;潘璐;李泽实;刘瑞芬;;湖北工业大学学报(第01期);第79-85页 * |
基于GIS的辽河中游城市非点源污染负荷时空估算新方法研究;滕凡全;;水利规划与设计(第06期);第45-50、107页 * |
海绵城市规划之径流污染控制研究分析――以株洲市为例;石勇;;中国市政工程(第02期);第72-78页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112417657A (zh) | 2021-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112417657B (zh) | 基于不同下垫面污染物贡献率的海绵城市优化设计方法 | |
Bisht et al. | Modeling urban floods and drainage using SWMM and MIKE URBAN: a case study | |
Jeong et al. | Development and integration of sub-hourly rainfall–runoff modeling capability within a watershed model | |
Yang et al. | Quantification of groundwater recharge in the city of Nottingham, UK | |
Li et al. | Modeling the Quality and Quantity of Runoff in a Highly Urbanized Catchment Using Storm Water Management Model. | |
Li et al. | Simulation of the hydrological and environmental effects of a sponge city based on MIKE FLOOD | |
CN111062125B (zh) | 海绵型综合管廊水文效应评估方法 | |
CN102201034A (zh) | 一种临河村落非点源污染负荷的估算方法 | |
CN109559098B (zh) | 一种海绵城市试点区低影响开发设施模拟方法 | |
Alamdari et al. | Evaluating the impact of climate change on water quality and quantity in an urban watershed using an ensemble approach | |
Paithankar et al. | Investigating the hydrological performance of green roofs using storm water management model | |
Li et al. | Performance assessment of stormwater GI practices using artificial neural networks | |
CN106650044B (zh) | 一种确定swmm软件中地表径流污染物参数w的方法 | |
Clary et al. | Integration of low-impact development into the international stormwater BMP database | |
Saghafian et al. | Application of the WEPP model to determine sources of run‐off and sediment in a forested watershed | |
CN110544192A (zh) | 集镇区全年面源污染入河负荷的估算方法 | |
Mohammed et al. | Modeling the quality of sewage during the leaking of stormwater surface runoff to the sanitary sewer system using SWMM: a case study | |
CN113792367B (zh) | 基于PySWMM的排水系统多来源入流入渗和出渗量动态估算方法 | |
Shaneyfelt et al. | Hydrologic modeling of distributed stormwater control measure retrofit and examination of impact of subcatchment discretization in PCSWMM | |
Dougaheh et al. | Urban stormwater management using low-impact development control measures considering climate change | |
Zhang et al. | Restoration of urban water environment landscape system and SWMM technology integration in water-scarce cities | |
Zhang et al. | Regulation characteristics of underlying surface on runoff regime metrics and their spatial differences in typical urban communities across China | |
Li et al. | Scale transfer and simulation of the infiltration in chestnut soil in a semi-arid grassland basin | |
CN115712970A (zh) | 一种管网模型概化程度对参数敏感性与模型精度影响分析方法、系统及存储介质 | |
Arvand et al. | Evaluating the impact of LID-BMPs on urban runoff reduction in an urban sub-catchment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |