CN106596347A - 基于等流时线法不透水地表径流污染物冲刷效率测定方法 - Google Patents
基于等流时线法不透水地表径流污染物冲刷效率测定方法 Download PDFInfo
- Publication number
- CN106596347A CN106596347A CN201611101150.7A CN201611101150A CN106596347A CN 106596347 A CN106596347 A CN 106596347A CN 201611101150 A CN201611101150 A CN 201611101150A CN 106596347 A CN106596347 A CN 106596347A
- Authority
- CN
- China
- Prior art keywords
- rainfall
- pollutant
- sigma
- runoff
- concentration
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
Abstract
基于等流时线法不透水地表径流污染物冲刷效率测定方法,属于市政工程和环境监测的交叉领域。该发明在地表径流污染物监测的基础上通过对研究区域进行概化,确定以单位时间作为分割点确定径流速度、划分等流时线、确定等流时面积;其次根据等流时线绘出出水口流量过程线,将每一时间段内监测的地表径流污染物浓度Ct减去天然降雨的浓度C0后与对应的流量Qt相乘获得这一时段内污染物负荷量Mt;接着将各时段径流污染负荷的进行累加得到总负荷量∑M,将∑M与雨前地表清扫物的污染负荷量M进行比较得到δ值;最后通过多场降雨得到δ的取值范围,再将δ反推估算未来某场降雨形成的地表径流污染物负荷∑M的大小范围。该方法简单易懂,并得到应用和验证。
Description
技术领域
本发明涉及基于等流时线法不透水地表径流污染物冲刷效率测定方法,属于市政工程技术和环境工程监测技术的交叉领域。
背景技术
随着城市化水平的不断加快,导致城市不透水地表面积的增加,城市热岛效应的形成,加之全球气候变暖,使得短历时强降雨频率的增加,容易导致城市内涝的发生,如2012年北京7·21事件、2015年南方大面积内涝,伴随城市内涝的发生致使地表径流污染物的形成,尤其是当不透水地表污染物经冲刷流落到城市的各个角落时给城市环境卫生造成了严重的危害,当形成径流通过城市管渠进入河湖容易造成水体的富营养化,据研究发现地表降雨径流污染已经是除工业污染、农业污染以外的第三大污染源。
地表径流污染物除了大气干湿沉降外主要来源于地表的累积,由于人类活动造成地表累积污染物的指标复杂,造成水体污染严重的常见污染物的指标有SS、COD、TN、TP和NH4-N。目前我国对地表径流污染物的研究已从传统的定性研究开始向模型的概化研究转变,主要通过计算机模型软件模拟地表污染物的迁移过程,从而计算出污染物在一次降雨中的冲刷量,但由于目前大部分计算机模型如SWMM水质模型是根据国外的环境特点、气候特点设计的,造成参数上与我国具有差异性,无法真实地反映污染的冲刷过程,同时不同降雨强度对地表的冲刷能力也不同,得到的污染物的浓度也就有差异,从而使得冲刷效率不同,对地表冲刷能力的核算也具有影响
发明内容
针对上述问题,本发明涉及基于等流时线法不透水地表径流污染物冲刷效率测定方法,属于市政工程技术和环境工程监测技术的交叉领域。该发明在地表径流污染物监测的基础上通过对研究区域进行概化,运用数学理论模型进行分布研究,首先确定研究区域下垫面类型(N)、汇水面积(A)、平均坡度(S)和地面流行距离(L)等参数,通过经验公式计算出汇流时间tc,根据地形特点以单位时间作为分割点确定径流速度、划分等流时线、确定等流时面积;其次根据等流时线绘出出水口流量过程线,将每一时间段内监测的地表径流污染物浓度Ct减去天然降雨的浓度C0后与对应的流量Qt相乘获得这一时段内污染物负荷量Mt;接着将各时段径流污染负荷的进行累加得到总负荷量∑M,将∑M与雨前清扫物的污染负荷量M比较得到δ值;最后通过多场降雨得到δ的取值范围,再将δ反推估算未来某场降雨形成的地表径流污染物负荷∑M的大小范围。
1.基于等流时线法不透水地表径流污染物冲刷效率测定方法,其特征在于步骤如下:
(1)根据等流时线法得到各时刻下出水口流量,将出水口流量和对应浓度进行相乘得到径流污染物负荷总量∑M:
式中Qt表示第t时刻出水口的流量;Ct表示t时刻时段出水口污染物的浓度;Co表示天然降雨水质浓度;∑M表示研究区域径流污染物负荷总量,n表示降雨场次;
(2)接着将雨前清扫样品的污染物负荷量M同研究区域污染物的冲刷量(∑M)进行比较,得到冲刷效率δ值,即:
其中δ表示一场降雨地表污染物的冲刷效率;
(3)最后通过n场降雨次数获得n个冲刷效率数值δi,其中i=1、2、3、…、n,,δ=(minδi,maxδj);i,j∈(1,n)
通过降雨前地表清扫物的含量预测每场降雨形成的地表径流污染负荷量∑M的取值范围:
∑M=(minδi,maxδj)×M (3)
式中M表示未来某场降雨前地表清扫物的污染负荷量;minδi表示前n场降雨中冲刷效率最小的值;maxδj表示前n场降雨中冲刷效率最大的值;i,j∈(1,n)。
附图说明
图1流域等流时线示意图
图2不同净雨历时情况下的径流过程;
图3技术流程图。
具体实施方式
(1)针对不透水地表污染物冲刷负荷效率的测定方面研究,目前还没有专门的研究,该方法简单易懂,适合城市大小流域的研究;
(2)等流时线法属于水文学中常用的数学模型,应用广泛,具有科学性和实用性。本发明属于首次运用到水质模型理论中,为定量研究地表径流污染物负荷量具有重要意义。
基于等流时线法不透水地表径流污染物冲刷效率测定方法,其特征在于:所述方法具体步骤如下:
雨前工作:
(1)确定研究区域的范围,统计出研究区域各不透水地表的面积大小以及其所占总面积比例;
(2)根据研究区域下垫面类型N、汇水面积A、平均坡度S和地面流行距离L计算出汇流时间tc,接着根据下垫面坡度以及单位时间△t绘制等流时线,并确定等流时面积(f1、f2、f3、f4…fm),见图1;
(3)降雨前将不同下垫面按其面积比例清扫出不同比例面积的地面沉积物样品,通常采用长和宽为1m的正方形的塑料采集器,将在同一下垫面下雨前清扫出的样品一起溶解于1L的容器中加蒸馏水浸泡24小时后进行检测,将检测的各污染物指标转化成质量化即“mg/m2”,最后将各不同垫面下采集到的各污染物的质量通过加权平均法求出单位面积的平均质量,接着乘以研究区域的面积,从而知道雨前总研究区域的各污染物的质量M。
雨中工作:
(4)降雨时在研究区域出水口收集地表径流样品,在前一小时以内每隔5分钟收集一次,超出一小时之后以每隔20分钟收集一次,直到降雨结束,在收集完的径流水样应该在24小时之内完成检测,同时在研究区域安装一台雨量计,记录降雨强度,获取每分钟降雨量。当研究区域有两个或者多个出水口时,应对研究区域的汇流面积按照实际坡度和水流方向进行分解,再重复(2)(3)步骤;
雨后工作:
(5)根据记录的每分钟降雨量采用等流时线法绘制出水口的流量Qt(t=1、2、3、4、5…),将每一时间段内监测的地表径流污染物浓度Ct(t=1、2、3、4、5…)减去天然降雨的浓度C0后与对应的流量Qt相乘获得这一时段内污染物负荷量Mt,再将各时段污染物的量进行累加,即为研究区域径流污染物的负荷总量
根据图2中a、b和c三种不同情况下的降雨情况
1)当tc<τ时(tc表示净雨历时,τ表示流域汇流时间)
以图2(a)所示,假设一次均匀净雨的历史tc=△t(h);净雨深为h(mm);净雨强度i=h/△t(mm/h)。
式中K——单位换算系数,当流量Q以m3/s计,h以mm计,f以km2计,如△t以小时计时,则K=0.278;如△t以min计时,则K=16.7;f1表示第1个等流时面积;h表示降雨强度;Q1表示第1时刻出水口的流量;Qt表示降雨时间t时刻下的出水口流量;Ct表示t时刻时段出水口污染物的浓度;CO表示天然降雨水质浓度;∑M表示研究区域径流污染物负荷总量;
2)当tc=τ时(tc表示净雨历时,τ表示流域汇流时间)
以图2(b)所示,当tc=4△t
3)当tc>τ时(tc表示净雨历时,τ表示流域汇流时间)
如图2(c)所示,假设一次均匀净雨历时tc=5△t,洪峰流量Qm由全部流域面积F上的部分净雨汇集而成,Qm的数值与tc=τ时求得的洪峰流量相同,只是它多延续了一个tc-τ=△t的时段,从而使图2(c)所示的径流过程线呈梯形形状。
综合1)、2)和3)三种根据实际降雨情况对照选择,最后将流量和对应浓度进行相乘得到径流污染物负荷总量∑M:
式中Qt表示第t时刻出水口的流量;Ct表示t时刻时段出水口污染物的浓度;CO表示天然降雨水质浓度;∑M表示研究区域径流污染物负荷总量;
(6)最后将雨前清扫样品的污染物负荷量M同研究区域污染物的冲刷量(∑M)进行比较,得到冲刷效率δ值,即:
其中δ表示一场降雨地表污染物的冲刷效率。
(7)通过n场降雨次数获得n个冲刷效率数值δi(i=1、2、3、…、n,这里n建议最少取到10,即有十场及以上降雨场次可保证冲刷效率取值范围的准确性),δ=(minδi,maxδj);i,j∈(1,n)
通过降雨前后地表清扫物的含量可预测每场降雨形成的地表径流污染负荷量∑M的取值范围:
∑M=(minδi,maxδj)×M (9)
式中M表示未来某场降雨前地表清扫物的污染负荷量minδi表示前n场降雨中冲刷效率最小的值;maxδj表示前n场降雨中冲刷效率最大的值;i,j∈(1,n)。
该方法是基于水文模型等流时线法提出的,模型简单易懂,在2016年北京工业大学实验室某区域污染物研究方面得到应用和验证,效果良好。
Claims (1)
1.基于等流时线法不透水地表径流污染物冲刷效率测定方法,其特征在于步骤如下:
(1)根据等流时线法得到各时刻下出水口流量,将出水口流量和对应浓度进行相乘得到径流污染物负荷总量∑M:
式中Qt表示第t时刻出水口的流量;Ct表示t时刻时段出水口污染物的浓度;Co表示天然降雨水质浓度;∑M表示研究区域径流污染物负荷总量,n表示降雨场次;
(2)接着将雨前清扫样品的污染物负荷量M同研究区域污染物的冲刷量(∑M)进行比较,得到冲刷效率δ值,即:
其中δ表示一场降雨地表污染物的冲刷效率;
(3)最后通过n场降雨次数获得n个冲刷效率数值δi,其中i=1、2、3、…、n,,δ=(minδi,maxδj);i,j∈(1,n)
通过降雨前地表清扫物的含量预测每场降雨形成的地表径流污染负荷量∑M的取值范围:
∑M=(minδi,maxδj)×M (3)
式中M表示未来某场降雨前地表清扫物的污染负荷量;minδi表示前n场降雨中冲刷效率最小的值;maxδj表示前n场降雨中冲刷效率最大的值;i,j∈(1,n)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611101150.7A CN106596347A (zh) | 2016-12-05 | 2016-12-05 | 基于等流时线法不透水地表径流污染物冲刷效率测定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611101150.7A CN106596347A (zh) | 2016-12-05 | 2016-12-05 | 基于等流时线法不透水地表径流污染物冲刷效率测定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106596347A true CN106596347A (zh) | 2017-04-26 |
Family
ID=58595567
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611101150.7A Pending CN106596347A (zh) | 2016-12-05 | 2016-12-05 | 基于等流时线法不透水地表径流污染物冲刷效率测定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106596347A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108009399A (zh) * | 2017-12-27 | 2018-05-08 | 贵州东方世纪科技股份有限公司 | 一种流域降雨的等流时线计算方法 |
CN108827405A (zh) * | 2018-06-07 | 2018-11-16 | 中国农业科学院农业资源与农业区划研究所 | 农业面源污染径流监测装置与方法 |
CN109272143A (zh) * | 2018-08-15 | 2019-01-25 | 青岛理工大学 | 一种基于未来天气预报的山区汇水量预测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102930153A (zh) * | 2012-10-26 | 2013-02-13 | 中国科学院生态环境研究中心 | 一种地表街尘及其负载的重金属降雨冲刷污染潜势的估算方法 |
CN103020424A (zh) * | 2012-11-22 | 2013-04-03 | 北京师范大学 | 一种基于人工降雨实验的北方平原农田区非点源污染负荷估算方法 |
-
2016
- 2016-12-05 CN CN201611101150.7A patent/CN106596347A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102930153A (zh) * | 2012-10-26 | 2013-02-13 | 中国科学院生态环境研究中心 | 一种地表街尘及其负载的重金属降雨冲刷污染潜势的估算方法 |
CN103020424A (zh) * | 2012-11-22 | 2013-04-03 | 北京师范大学 | 一种基于人工降雨实验的北方平原农田区非点源污染负荷估算方法 |
Non-Patent Citations (5)
Title |
---|
VASSILIOS A.ETC: "Modeling and Management of Urban Stormwater Runoff Quality: A Review", <WATER RESOURCE MANAGEMENT> * |
李立青,尹澄清,孔玲莉,何庆慈: "2次降雨间隔时间对城市地表径流污染负荷的影响", 《环境科学》 * |
潘安君,张书函,陈建: "《城市雨水综合利用技术研究与应用》", 30 August 2010, 中国水利水电出版社 * |
王小梅: "北京地区街尘-径流污染特征及潜在污染负荷估算", 《中国优秀硕士学位论文全文数据库工程科技1辑》 * |
韩龙,余麟: "深圳市建成区地表污染物累积特征研究", 《广东化工》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108009399A (zh) * | 2017-12-27 | 2018-05-08 | 贵州东方世纪科技股份有限公司 | 一种流域降雨的等流时线计算方法 |
CN108009399B (zh) * | 2017-12-27 | 2021-09-17 | 贵州东方世纪科技股份有限公司 | 一种流域降雨的等流时线计算方法 |
CN108827405A (zh) * | 2018-06-07 | 2018-11-16 | 中国农业科学院农业资源与农业区划研究所 | 农业面源污染径流监测装置与方法 |
CN109272143A (zh) * | 2018-08-15 | 2019-01-25 | 青岛理工大学 | 一种基于未来天气预报的山区汇水量预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111784052B (zh) | 一种城市面源污染入河负荷量预测方法 | |
Goodwin et al. | The temporal and spatial variability of sediment transport and yields within the Bradford Beck catchment, West Yorkshire | |
Jalón-Rojas et al. | Turbidity in the fluvial Gironde Estuary (southwest France) based on 10-year continuous monitoring: sensitivity to hydrological conditions | |
CN102201034B (zh) | 一种临河村落非点源污染负荷的估算方法 | |
CN106596347A (zh) | 基于等流时线法不透水地表径流污染物冲刷效率测定方法 | |
CN106874625B (zh) | 一种前置库容积的确定方法及确定系统 | |
CN112084671B (zh) | 城市时变增益降雨-径流过程模拟计算方法 | |
Barbé et al. | Modeling the buildup and washoff of pollutants on urban watersheds 1 | |
Horner et al. | Hydrologic monitoring of the Seattle ultra-urban stormwater management projects: Summary of the 2000–2003 water years | |
CN109583035A (zh) | 基于元胞自动机的城市面源污染物累积冲刷过程计算方法 | |
Wang et al. | Hydrological modeling and field validation of a bioretention basin | |
CN109979172A (zh) | 一种基于新安江模型的动态山洪临界雨量预报方法 | |
CN106250659A (zh) | 一种基于降雨预测合流制溢流水力学行为的方法及其应用 | |
van den Elsen et al. | Discharge and sediment measurements at the outlet of a watershed on the Loess plateau of China | |
Ryan et al. | Changes in streambed sediment characteristics and solute transport in the headwaters of Valley Creek, an urbanizing watershed | |
CN111291305A (zh) | 基于一维水质模型计算流域污染源入河系数的方法 | |
Chang et al. | Use of intercepted runoff depth for stormwater runoff management in industrial parks in Taiwan | |
Ahmad et al. | Nonpoint source pollution | |
Hixon et al. | State of the practice: Assessing water quality benefits from street sweeping | |
Shaikh et al. | Impacts of urbanization on watershed hydrology and runoff water quality of a watershed: a review | |
Afzalimehr et al. | Field investigation on friction factor in mountainous cobble-bed and boulder-bed rivers | |
Schaffner | Numerical investigations on the function of flush waves in a reservoir sewer | |
Sulistyono et al. | An analysis of the influence of variability rainfall on flow rate based on the watershed characteristics | |
Rim | Analyzing runoff dynamics of paved soil surface using weighable lysimeters | |
Sutherland et al. | Sophisticated Stormwater Quality Modeling is Worth the Effort. |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170426 |