CN106709160B - 一种确定swmm软件中累积参数的方法 - Google Patents

一种确定swmm软件中累积参数的方法 Download PDF

Info

Publication number
CN106709160B
CN106709160B CN201611101142.2A CN201611101142A CN106709160B CN 106709160 B CN106709160 B CN 106709160B CN 201611101142 A CN201611101142 A CN 201611101142A CN 106709160 B CN106709160 B CN 106709160B
Authority
CN
China
Prior art keywords
days
formula
determining
swmm
pollutants
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.)
Expired - Fee Related
Application number
CN201611101142.2A
Other languages
English (en)
Other versions
CN106709160A (zh
Inventor
周玉文
张静
赵树旗
卢兴超
刘原
娄富豪
吴献平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nantong Saiyasi Technology Co ltd
Original Assignee
Beijing University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201611101142.2A priority Critical patent/CN106709160B/zh
Publication of CN106709160A publication Critical patent/CN106709160A/zh
Application granted granted Critical
Publication of CN106709160B publication Critical patent/CN106709160B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00Water conservation; Efficient water supply; Efficient water use
    • Y02A20/152Water filtration

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

一种确定SWMM软件中累积参数确定的方法,属于市政工程、环境工程、海绵城市建设和计算机模型软件参数确定的交叉领域。本发明发明了SWMM软件中累积参数C1和C2的确定方法。采用模拟实验和现场采样相结合,确定前期晴天天数和地表沉积物累积量变化关系,从而建立前期晴天数和地表沉积物累积量的确定方法,并确定不同污染物C1的取值范围,通过相关计算确定累积参数C2的取值范围。该发明可为SWMM软件在降雨径流水质模拟中的实际应用提供数据支持和技术帮助。

Description

一种确定SWMM软件中累积参数的方法
技术领域
本发明涉及一种确定SWMM软件中累积参数确定的方法,属于市政工程、环境工程、海绵城市建设和计算机模型软件参数确定的交叉领域。
背景技术
SWMM(Storm Water Management Model)又称暴雨洪水管理模型是美国环境保护署开发的一款动态的降水-径流模拟模型,除了用于模拟城市某一单一降水事件或长期的水量也可以用于水质模拟。其径流模块部分综合处理各子流域所发生的降水,径流和污染负荷。其汇流模块部分则通过管网、渠道、蓄水和处理设施、水泵、调节闸等进行水量传输。该模型可以跟踪模拟不同时间步长任意时刻每个子流域所产生径流的水质和水量,以及每个管道和河道中水的流量、水深及水质等情况。
SWMM模型可以动态的模拟降雨径流,但在模拟水流水质的单一或者连续方面较为常用,即城市面源污染的研究方面,该模型是模拟评估使得面源污染负荷有效降低的最佳管理措施(BMPs)。在美国、加拿大和欧洲等国家,该模型都被用来解决当地的排水的水量和水质问题,在分析有关防止雨污合流制管道的溢流情况、城市大降水管理规划以及降低污染负荷方面时,该模型都能得到较好的模拟效果。在我国随着海绵城市建设的推广,城市生态系统的规划,借助计算机模型可以反映城市地表径流水质的运移情况,在添加LID措施的作用下通过对SWMM模型水质参数的设定预测海绵城市建设后的水质情况,从而反映LID措施在海绵城市规划建设中的效果。
SWMM模型参数是美国环保署根据美国的土地利用类型、地表污染状况以及环境质量等因素确定的,其使用上具有局限性。由地形和城市生活环境的差异性,其参数所具有的局限性导致无法直接在我国城市地表径流水质中进行模拟和使用。目前在国内SWMM水质模型的使用上,多数照搬国外SWMM模型手册上的参数,没有适合我国国情的一套参数,同时也没有一套关于SWMM水质参数的率定方法。本专利旨在提供一套根据我国国情确定SWMM水质冲刷模型参数W为污染物浓度(mg/L);B为单位面积污染物沉积量(kg/104m2);参数S1为冲刷系数、S2为冲刷指数;参数q为单位面积径流量(mm/h);累积模型中C1为最大增长可能和C2为增长速率常数等参数的方法。在原有的SWMM水质模型的基础上对其进行继承和改进,并提出一套完整的测定方法,分别为:《一种直接获得SWMM水质模型中参数W的方法》、《一种确定SWMM软件中地表径流污染物参数W的方法》、《一种确定SWMM中含不同下垫面地表径流污染物参数W的方法》、《一种受大气污染影响自然降雨水质参数W1的确定方法》、《一种无大气影响下SWMM中地表径流水质参数W0的确定方法》、《一种直接获得SWMM水质模型参数B的方法》、《一种确定SWMM软件中累积参数确定的方法》、《一种确定SWMM水质冲刷模型参数的方法》,共计八个测定方法。
由于地表沉积污染物受大气污染物累积的影响,大气中累积污染物在雨天经雨水的冲刷作用下进入地表径流,与地表沉积污染物共同成为城市面源污染的重要来源。在SWMM软件的累积参数C1和C2,考虑大气污染累积因素,可提高该软件对降雨径流水质模拟的准确度,有利于SWMM软件在降雨径流水质模拟中的应用。
本发明建立了SWMM软件的累积参数C1、C2的确定方法,将大气污染物和地表沉积污染物合并为城市降雨径流的污染源,建立大气污染物和地表沉积物最大累积量的确定方法,为SWMM软件在降雨径流水质模拟中的应用提供帮助。同时,不同的污染物由于其沉积性能不同,在大气和地表累积规律也存在显著差异,因此本发明针对大气和地表沉积的不同的污染物进行了累积参数取值范围的确定,提高SWMM软件在降雨径流水质模拟中的精度。
发明内容
本发明为提高SWMM软件在降雨径流水质模拟中的精度,提供一种水质参数的确定方法。
SWMM软件中累积参数的确定是制约SWMM软件在降雨径流水质模拟中的关键因素。SWMM软件中累积模型包括幂函数模型、指数函数模型、饱和浸润方程和外部时间序列。其中指数函数模型(见式(1))受到大家普遍关注,该模型中最大累积量C1和累积速率常数C2的确定是该模型应用的关键所在。
Figure GDA0002542584710000031
式中:C1——最大累积量,g/m2
C2——累积速率常数,d-1
SWMM软件中累积参数的确定是制约SWMM软件在降雨径流水质模拟中的关键因素。SWMM软件中累积模型包括幂函数模型、指数函数模型、饱和浸润方程和外部时间序列。其中指数函数模型(见式(1))受到大家普遍关注,该模型中最大累积量C1和累积速率常数C2的确定是该模型是否能够准确模拟的关键所在。
Figure GDA0002542584710000032
式中:C1——最大累积量,g/m2
C2——累积速率常数,d-1
一种确定SWMM软件中累积参数确定的方法,应用以下模型:
Figure GDA0002542584710000033
式中:C1——最大累积量,g/m2
C2——累积速率常数,d-1
B为单位面积污染物沉积量,kg/104m2
其特征在于其步骤如下:
①分别在研究区域的不同功能区选择采样点;
②在每个采样点铺设8块地表沉积物采样装置,采样装置由厚度为2.5mm油毡及边框组成,大小为1.2m×0.6m;
③分别经过1、3、5、7、15、30、45、60天后取回采样装置,每次取回1块,采用冲刷装置对采样装置进行冲刷;
④冲刷过程:布水器对下垫面冲刷,下垫面为上述采样装置;下垫面较低一端设有收集槽,用于收集冲刷液;
⑤冲刷液搅拌促溶24h后,测定冲刷液中各类污染物浓度;
⑥计算不同功能区不同前期干旱天数下各类污染物累积总量,污染物累积总量计算公式见式(2);
Figure GDA0002542584710000041
式中:L——各采样点地表沉积物总量,g/m2
Xi——冲刷液中污染物浓度,mg/L;
V——冲刷液体积,m3
⑦建立不同功能区地表沉积污染物累积总量与前期晴天天数曲线,确定最大累积量C1的值;
⑧计算累积参数C2的值;
通过转化,式(1)书写为:
Figure GDA0002542584710000042
将式(3)两遍取自然对数,获得关系式(4):
Figure GDA0002542584710000043
通过建立
Figure GDA0002542584710000044
与t之间的直线关系,获得各类污染物累积速率常数C2的取值范围。
本发明收集器的优势:
(1)解决SWMM软件中水质参数确定方法欠缺问题;
(2)将大气污染物和地表沉积污染物合并为城市降雨径流的污染源,建立大气污染物和地表沉积物最大累积量的确定方法;
(3)在较大尺度范围内采用自然降雨污染物总量间接反应大气污染物总量,建立大气污染物累积量与前期晴天天数的相关关系;
(4)建立了地表沉积污染物累积量与前期晴天天数相关关系确定方法,为SWMM模型中相关参数的确定提供技术支持;
(5)确定不同污染物水质参数取值范围,为SWMM软件准确模拟提供支持。
附图说明
图1为地表沉积物采样装置图;
图2为地表沉积物冲刷装置图。
具体实施方式
目前对该模型中C1和C2的确定较为困难,因此使该模型的使用受到限制。本研究发明了一种实验方法确定SWMM软件模拟降雨径流水质指数函数模型所需要的参数的确定方法,其步骤如下:
①分别在研究区域的不同功能区(如生活区、工业区、商业区、文教区等)选择采样点,采样点一般位于距离路边1m处;
②在每个采样点铺设8块地表沉积物采样装置(见附图1),采样装置由厚度为2.5mm油毡及边框组成,大小为1.2m×0.6m;
③分别经过1、3、5、7、15、30、45、60天后取回采样装置,每次取回1块,采用冲刷装置对采样装置进行冲刷。
④冲刷过程中采用自来水代替雨水,布水器的一端接自来水水龙头,带有压力的自来水通过两个圆形(直径为20cm)布水器对下垫面冲刷2h,冲刷强度采用600mL/min,下垫面长1.2m,宽0.6m,坡度为1‰。下垫面较低一端开直径为10cm圆形收集槽,用于收集冲刷液。
⑤冲刷液搅拌促溶24h后,测定冲刷液中各类污染物浓度;
⑥计算不同功能区不同前期干旱天数下各类污染物累积总量,污染物累积总量计算公式见式(2);
Figure GDA0002542584710000051
式中:L——各采样点地表沉积物总量,g/m2
Xi——冲刷液中各类污染物平均浓度,mg/L;
V——冲刷液体积,m3
⑦建立不同功能区地表沉积污染物累积总量与前期晴天天数曲线,确定最大累积量C1的值;
⑧计算累积参数C2的值。
通过转化,式(1)可书写为:
Figure GDA0002542584710000061
将式(3)两遍取自然对数,获得关系式(4):
Figure GDA0002542584710000062
通过建立
Figure GDA0002542584710000063
与t之间的直线关系,可获得各类污染物累积速率常数C2的取值范围。
具体效果
课题组在河北科技大学校园及周边选定的文教区、生活区、商业区、工业区的采样点分别铺设八块1.2m×0.6m油毡,并标注“实验区域,请勿清扫”字样,分别经过1、3、5、7、15、30、45、60天后取回油毡,采用冲刷装置对油毡进行冲刷,收集冲刷液并搅拌促溶24h;而后测定冲刷液中SS、COD和氨氮浓度,计算不同功能区污染物不同累积天数下污染物累积总量,建立不同功能区地表沉积污染物累积总量与前期晴天天数曲线,确定SS、COD、氨氮最大累积量C1取值范围为:6.37-12.86g/m2、1.403-2.166g/m2和8.89-12.88mg/m2,并通过式(4)计算出,SS、COD、氨氮最大累积量C2取值范围为:0.055-0.100d-1、0.061-0.085d-1、0.133-0.167d-1

Claims (1)

1.一种确定SWMM软件中累积参数的方法,应用以下模型:
Figure FDA0002542584700000011
式中:C1——最大累积量,g/m2
C2——累积速率常数,d-1
B为单位面积污染物沉积量,kg/104m2
其特征在于其步骤如下:
①分别在研究区域的不同功能区选择采样点;
②在每个采样点铺设8块地表沉积物采样装置,采样装置由厚度为2.5mm油毡及边框组成,大小为1.2m×0.6m;
③分别经过1、3、5、7、15、30、45、60天后取回采样装置,每次取回1块,采用冲刷装置对采样装置进行冲刷;
④冲刷过程:布水器对下垫面冲刷,下垫面为上述采样装置;下垫面较低一端设有收集槽,用于收集冲刷液;
⑤冲刷液搅拌促溶24h后,测定冲刷液中各类污染物浓度;
⑥计算不同功能区不同前期干旱天数下各类污染物累积总量,污染物累积总量计算公式见式(2);
Figure FDA0002542584700000012
式中:L——各采样点地表沉积物总量,g/m2
Xi——冲刷液中污染物浓度,mg/L;
V——冲刷液体积,m3
⑦建立不同功能区地表沉积污染物累积总量与前期晴天天数曲线,确定最大累积量C1的值;
⑧计算累积参数C2的值;
通过转化,式(1)书写为:
Figure FDA0002542584700000013
将式(3)两遍取自然对数,获得关系式(4):
Figure FDA0002542584700000014
通过建立
Figure FDA0002542584700000021
与t之间的直线关系,获得各类污染物累积速率常数C2的取值范围。
CN201611101142.2A 2016-12-05 2016-12-05 一种确定swmm软件中累积参数的方法 Expired - Fee Related CN106709160B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611101142.2A CN106709160B (zh) 2016-12-05 2016-12-05 一种确定swmm软件中累积参数的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611101142.2A CN106709160B (zh) 2016-12-05 2016-12-05 一种确定swmm软件中累积参数的方法

Publications (2)

Publication Number Publication Date
CN106709160A CN106709160A (zh) 2017-05-24
CN106709160B true CN106709160B (zh) 2020-11-03

Family

ID=58935557

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611101142.2A Expired - Fee Related CN106709160B (zh) 2016-12-05 2016-12-05 一种确定swmm软件中累积参数的方法

Country Status (1)

Country Link
CN (1) CN106709160B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110867130A (zh) * 2018-08-28 2020-03-06 广州八爪鱼教育科技有限公司 海绵城市课堂教学模型

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201780296U (zh) * 2010-07-08 2011-03-30 常州市环境监测中心站 一种用于城市降雨径流污染研究的物理模型试验装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101858107B (zh) * 2010-06-09 2011-12-21 北京建筑工程学院 雨水处理和利用设施的设计方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201780296U (zh) * 2010-07-08 2011-03-30 常州市环境监测中心站 一种用于城市降雨径流污染研究的物理模型试验装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"SWMM中典型水质参数值确定方法的研究";李海燕等;《给水排水》;20110731;第37卷(第S1期);正文第159-162页第1-3节 *
"北方城市天然降雨中氨氮污染特征研究";张静等;《环境科学与技术》;20160630;第39卷(第6期);正文第129-134页摘要,第2节 *
"城市地表污染物累积模型研究";边博;《土木建筑与环境工程》;20101231;第32卷(第6期);正文第137-141页第1-2节 *

Also Published As

Publication number Publication date
CN106709160A (zh) 2017-05-24

Similar Documents

Publication Publication Date Title
Liu et al. Enhancing a rainfall-runoff model to assess the impacts of BMPs and LID practices on storm runoff
Abi Aad et al. Modeling techniques of best management practices: Rain barrels and rain gardens using EPA SWMM-5
Nie et al. Impacts of climate change on urban drainage systems–a case study in Fredrikstad, Norway
Arjenaki et al. Modeling and investigating the effect of the LID methods on collection network of urban runoff using the SWMM model (case study: Shahrekord City)
Lee et al. Comparative evaluation of runoff and water quality using HSPF and SWMM
Ren et al. Evaluating the stormwater management model to improve urban water allocation system in drought conditions
CN108984823B (zh) 一种合流制溢流调蓄池规模的确定方法
Li et al. Modeling the Quality and Quantity of Runoff in a Highly Urbanized Catchment Using Storm Water Management Model.
CN106650044B (zh) 一种确定swmm软件中地表径流污染物参数w的方法
CN109522645B (zh) 用于河流流域的水质模拟方法及系统
CN110544192B (zh) 集镇区全年面源污染入河负荷的估算方法
Ahamed et al. Urban flood modeling and management using SWMM for new RR Pet region, Vijayawada, India
Shuster et al. Front-loading urban stormwater management for success–a perspective incorporating current studies on the implementation of retrofit low-impact development
Taji et al. LID coupled design of drainage model using GIS and SWMM
CN101858107B (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
Thorndahl et al. Assessment of runoff contributing catchment areas in rainfall runoff modelling
Obaid et al. Modeling sewer overflow of a city with a large floating population
Zhang et al. Restoration of urban water environment landscape system and SWMM technology integration in water-scarce cities
CN106709160B (zh) 一种确定swmm软件中累积参数的方法
Zeleňáková et al. Impact of land use changes on surface runoff in urban areas-Case study of Myslavsky Creek Basin in Slovakia.
Shon et al. Analysis of the characteristics of non-point pollutant runoff applied LID techniques in industrial area
Greenberg Urban Hydrological Modeling of the Malden River Using the Storm Water Management Model (SWMM)
CN106771026A (zh) 一种无大气影响下swmm中地表径流水质参数w0的确定方法
CN106682271B (zh) 一种确定swmm水质冲刷模型参数的方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210218

Address after: Room 1424, Floor 4, Peony Pioneer Building, No. 2 Garden Road, Haidian District, Beijing, 100191

Patentee after: Beijing Zhonglian Technology Service Co.,Ltd.

Address before: 100124 No. 100 Chaoyang District Ping Tian Park, Beijing

Patentee before: Beijing University of Technology

TR01 Transfer of patent right

Effective date of registration: 20210319

Address after: 1205-30, 12 / F, building 3, yard 8, East Road of Automobile Museum, Fengtai District, Beijing

Patentee after: Li'an (Beijing) Technology Service Group Co.,Ltd.

Address before: Room 1424, Floor 4, Peony Pioneer Building, No. 2 Garden Road, Haidian District, Beijing, 100191

Patentee before: Beijing Zhonglian Technology Service Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210425

Address after: 226600 No. 288 Changjiang West Road, Haian Town, Haian City, Nantong City, Jiangsu Province

Patentee after: NANTONG HAIJIA ITELLIGENT TECHNOLOGY Co.,Ltd.

Address before: 1205-30, 12 / F, building 3, yard 8, East Road of Automobile Museum, Fengtai District, Beijing

Patentee before: Li'an (Beijing) Technology Service Group Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220324

Address after: 226600 No. 300, Baichuan Road, hi tech Zone, Hai'an City, Nantong City, Jiangsu Province

Patentee after: Nantong saiyasi Technology Co.,Ltd.

Address before: 226600 No.288, Changjiang West Road, Hai'an town, Hai'an City, Nantong City, Jiangsu Province

Patentee before: NANTONG HAIJIA ITELLIGENT TECHNOLOGY Co.,Ltd.

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20201103