CN109977601B - 分流制区域控制降雨径流和混接污染的截流厚度确定方法 - Google Patents

分流制区域控制降雨径流和混接污染的截流厚度确定方法 Download PDF

Info

Publication number
CN109977601B
CN109977601B CN201910283826.6A CN201910283826A CN109977601B CN 109977601 B CN109977601 B CN 109977601B CN 201910283826 A CN201910283826 A CN 201910283826A CN 109977601 B CN109977601 B CN 109977601B
Authority
CN
China
Prior art keywords
interception
thickness
emission
discharge
cod
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
Application number
CN201910283826.6A
Other languages
English (en)
Other versions
CN109977601A (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.)
Fuzhou Planning And Design Institute Group Co ltd
Original Assignee
Fuzhou Planning And Design Institute Group Co ltd
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 Fuzhou Planning And Design Institute Group Co ltd filed Critical Fuzhou Planning And Design Institute Group Co ltd
Priority to CN201910283826.6A priority Critical patent/CN109977601B/zh
Publication of CN109977601A publication Critical patent/CN109977601A/zh
Application granted granted Critical
Publication of CN109977601B publication Critical patent/CN109977601B/zh
Active 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

Abstract

本发明涉及一种分流制区域控制降雨径流污染和混接污染的截流厚度确定方法,首先,收集排水管网和下垫面数据,现场监测不同降雨强度下排口的出水量和出水水质;然后,建立排水管网模型,并对模型进行率定和验证,确定COD、NH3‑N和TP三种水质指标的排放要求,利用数学函数预测排口COD、NH3‑N和TP三种污染物随时间推移的变化过程,确定污染物排放达标时刻,利用排水管网模型模拟排口出流过程,计算截流量;最后计算不同降雨强度下满足三种污染物排放要求的截流厚度,并计算不同降雨强度下满足三种污染物排放要求的截流厚度的平均值。本发明能够针对特地地区特定条件下,确定与水环境治理目标匹配、经济上合理的截流厚度。

Description

分流制区域控制降雨径流和混接污染的截流厚度确定方法
技术领域
本发明涉及黑臭水体治理、水环境提升技术领域,特别是一种分流制区域控制降雨径流污染和混接污染的截流厚度确定方法。
背景技术
分流制排水系统是我国主推的排水系统建设方式,降雨径流带来的面源污染、雨污混接产生污水排河是造成分流制区域水环境污染和水体黑臭的主要原因。沿河两岸建设截污系统(含初期雨水调蓄池)是解决分流制区域水环境污染的主要措施。截流厚度是规划设计截流系统的重要参数,在数值上等于截流量与对应集水区的总面积比值,单位为mm。
截流厚度受地区降雨特征、地表污染状况、服务范围和人口、雨污混接情况、受纳水体水质要求、受纳水体水环境容量等因素影响,不同地区存在显著的差异性。截流厚度取值过大,将加重污水处理厂的运行负荷,截流系统工程建设的投资会相应增大,造成不必要浪费;截流厚度取值偏小,达不到水环境治理目标。
目前,规划设计截污系统时,截流厚度一般参照《室外排水设计规范》(GB 50014-2006),取值4-8mm,截流厚度取值较主观,且无法准确定量说明与水环境治理目标的关系。因此,如何根据分流制区域的实际特征(降雨条件、服务人口、受纳水体水环境要求和水环境容量等),科学准确的确定截流厚度,对于黑臭水体治理和水环境提升具有重要的意义。
发明内容
有鉴于此,本发明的目的是提出一种分流制区域控制降雨径流污染和混接污染的截流厚度确定方法,能够针对特地地区特定条件下,确定与水环境治理目标匹配、经济上合理的截流厚度。
本发明采用以下方案实现:一种分流制区域控制降雨径流污染和混接污染的截流厚度确定方法,包括以下步骤:
步骤S1:获取分流制区域的排水管网和下垫面数据以及分流制区域面积,并分别获取小雨、中雨、大雨和暴雨强度下分流制区域排口的出水量和出水水质监测数据;
步骤S2:建立分流制区域排水管网模型,分别利用小雨、中雨、大雨和暴雨强度下的出水量确定和验证模型;
步骤S3:根据分流制区域内受纳水体水质要求,确定COD、NH3-N和TP三种污染物形成的三种水质指标的排放要求;
步骤S4:获取降雨排口水质监测数据,所述降雨排口水质监测数据包括每场降雨每种污染物的浓度以及每场降雨每种污染物的浓度对应的时间;利用所述降雨排口水质监测数据,采用最小二乘法拟合求解数学函数;改变上述数学函数中的每场降雨每种污染物的浓度对应的时间,分别得到降雨排口COD、NH3-N和TP三种污染物的排放浓度,当所述三种污染物的排放浓度达到各自设定的排放要求时,分别得到三种污染物排放达标时刻,记为化学需氧量排放达标时刻T(COD)i、氨氮排放达标时刻T(NH3-N)i和总磷排放达标时刻T(TP)i,i表示第i场降雨;
步骤S5:利用所述分流制区域排水管网模型,模拟分流制区域排口的流量过程,得到每个时刻的出流量,分别将所述化学需氧量排放达标时刻T(COD)i、所述氨氮排放达标时刻T(NH3-N)i和所述总磷排放达标时刻T(TP)i三个时刻之前的出流量累加,得到三种污染物排放的截流量,分别为化学需氧量排放截流量Q(COD)i、氨氮排放截流量Q(NH3-N)i和总磷排放截流量Q(TP)i,单位均为m3
步骤S6:将所述化学需氧量排放截流量Q(COD)i、所述氨氮排放截流量Q(NH3-N)i和所述总磷排放截流量Q(TP)i分别除以所述分流制区域面积,得到三种污染物排放截流厚度,分别为化学需氧量排放截流厚度H(COD)i、氨氮排放截流厚度H(NH3-N)i和总磷排放截流厚度H(TP)i,单位均为mm。
步骤S7:分别计算所述化学需氧量排放截流厚度H(COD)i、所述氨氮排放截流厚度H(NH3-N)i和所述总磷排放截流厚度H(TP)i的平均值,得到化学需氧量排放截流厚度平均值H(COD)、氨氮排放截流厚度平均值H(NH3-N)和总磷排放截流厚度平均值H(TP)。
进一步地,步骤S7所述计算化学需氧量排放截流厚度平均值H(COD)、氨氮排放截流厚度平均值H(NH3-N)和总磷排放截流厚度平均值H(TP)的具体计算公式如下:
Figure GDA0003619106850000031
Figure GDA0003619106850000032
Figure GDA0003619106850000033
式中:m为降雨场数。
进一步地,步骤S1中所述的排水管网包括但不限于雨水井、雨水排口、雨水管段。
进一步地,步骤S2中所述的建立分流制区域排水管网模型具体为:提供所述排水管网的拓扑结构,将所述雨水管段连接至所述雨水排口,并划分集水区;设置各个集水区的产流表面类型,选择非线性水库模型进行地表汇流演算,得到分流制区域排水管网模型。
进一步地,所述设置各个集水区的产流表面类型具体为:将所述各个集水区的产流表面类型分为两种,分别为非渗透表面和渗透表面两种类型。
与现有技术相比,本发明有以下有益效果:
本发明能够针对特地地区特定条件下,确定与水环境治理目标匹配、经济上合理的截流厚度。
附图说明
图1为本发明实施例的流程图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
如图1所示,本实施例提供了一种分流制区域控制降雨径流污染和混接污染的截流厚度确定方法,包括以下步骤:
步骤S1:获取分流制区域的排水管网和下垫面数据以及分流制区域面积,并分别获取(通过现场监测的方式获取)小雨、中雨、大雨和暴雨强度下分流制区域排口的出水量和出水水质监测数据;
步骤S2:建立分流制区域排水管网模型,分别利用小雨、中雨、大雨和暴雨强度下的出水量确定和验证模型;
步骤S3:根据分流制区域内受纳水体水质要求,通过查阅《地表水环境质量标准》中对应的地表水环境质量标准基本项目标准限值,确定COD、NH3-N和TP三种污染物形成的三种水质指标的排放要求;
步骤S4:获取降雨排口水质监测数据,所述降雨排口水质监测数据包括每场降雨每种污染物的浓度以及每场降雨每种污染物的浓度对应的时间;利用所述降雨排口水质监测数据,采用最小二乘法拟合求解数学函数;(以COD为例,利用COD浓度y1、y2…yi,以及对应的时刻t1、t2…ti,通过最小二乘法求解出数学函数中的参数,得到用于预测的数学函数,利用上述数学函数,通过改变时间t,得到降雨排口COD、NH3-N和TP三种污染物的排放浓度)改变上述数学函数中的每场降雨每种污染物的浓度对应的时间,分别得到降雨排口COD、NH3-N和TP三种污染物的排放浓度,当所述三种污染物的排放浓度达到各自设定的排放要求时,分别得到三种污染物排放达标时刻,记为化学需氧量排放达标时刻T(COD)i、氨氮排放达标时刻T(NH3-N)i和总磷排放达标时刻T(TP)i,i表示第i场降雨;
步骤S5:利用所述分流制区域排水管网模型,模拟分流制区域排口的流量过程,得到每个时刻的出流量,分别将所述化学需氧量排放达标时刻T(COD)i、所述氨氮排放达标时刻T(NH3-N)i和所述总磷排放达标时刻T(TP)i三个时刻之前的出流量累加,得到三种污染物排放的截流量,分别为化学需氧量排放截流量Q(COD)i、氨氮排放截流量Q(NH3-N)i和总磷排放截流量Q(TP)i,单位均为m3。
步骤S6:将所述化学需氧量排放截流量Q(COD)i、所述氨氮排放截流量Q(NH3-N)i和所述总磷排放截流量Q(TP)i分别除以所述分流制区域面积,得到三种污染物排放截流厚度,分别为化学需氧量排放截流厚度H(COD)i、氨氮排放截流厚度H(NH3-N)i和总磷排放截流厚度H(TP)i,单位均为mm。
步骤S7:分别计算所述化学需氧量排放截流厚度H(COD)i、所述氨氮排放截流厚度H(NH3-N)i和所述总磷排放截流厚度H(TP)i的平均值,得到化学需氧量排放截流厚度平均值H(COD)、氨氮排放截流厚度平均值H(NH3-N)和总磷排放截流厚度平均值H(TP)。
在本实施例中,步骤S7所述计算化学需氧量排放截流厚度平均值H(COD)、氨氮排放截流厚度平均值H(NH3-N)和总磷排放截流厚度平均值H(TP)的具体计算公式如下:
Figure GDA0003619106850000061
Figure GDA0003619106850000062
Figure GDA0003619106850000063
式中:m为降雨场数。
在本实施例中,步骤S1中所述的排水管网包括但不限于雨水井、雨水排口、雨水管段。
在本实施例中,步骤S2中所述的建立分流制区域排水管网模型具体为:提供所述排水管网的拓扑结构,确保雨水管线上下游关系正确,将所述雨水管段连接至所述雨水排口,并划分集水区;设置各个集水区的产流表面类型,选择非线性水库模型进行地表汇流演算,得到分流制区域排水管网模型。
在本实施例中,所述设置各个集水区的产流表面类型具体为:将所述各个集水区的产流表面类型分为两种,分别为非渗透表面和渗透表面两种类型。
以上所述仅为本发明的较佳实施例,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本发明的涵盖范围。

Claims (5)

1.一种分流制区域控制降雨径流污染和混接污染的截流厚度确定方法,其特征在于:包括以下步骤:
步骤S1:获取分流制区域的排水管网和下垫面数据以及分流制区域面积,并分别获取小雨、中雨、大雨和暴雨强度下分流制区域排口的出水量和出水水质监测数据;
步骤S2:建立分流制区域排水管网模型,分别利用小雨、中雨、大雨和暴雨强度下的出水量确定和验证模型;
步骤S3:根据分流制区域内受纳水体水质要求,确定COD、NH3-N和TP三种污染物三种水质指标的排放要求;
步骤S4:获取降雨排口水质监测数据,所述降雨排口水质监测数据包括每场降雨每种污染物的浓度以及每场降雨每种污染物的浓度对应的时间;利用所述降雨排口水质监测数据,采用最小二乘法拟合求解数学函数;改变上述数学函数中的每场降雨每种污染物的浓度对应的时间,分别得到降雨排口COD、NH3-N和TP三种污染物的排放浓度,当所述三种污染物的排放浓度达到各自设定的排放要求时,分别得到三种污染物排放达标时刻,记为化学需氧量排放达标时刻T(COD)i、氨氮排放达标时刻T(NH3-N)i和总磷排放达标时刻T(TP)i,i表示第i场降雨;
步骤S5:利用所述分流制区域排水管网模型,模拟分流制区域排口的流量过程,得到每个时刻的出流量,分别将所述化学需氧量排放达标时刻T(COD)i、所述氨氮排放达标时刻T(NH3-N)i和所述总磷排放达标时刻T(TP)i三个时刻之前的出流量累加,得到三种污染物排放的截流量,分别为化学需氧量排放截流量Q(COD)i、氨氮排放截流量Q(NH3-N)i和总磷排放截流量Q(TP)i,单位均为m3
步骤S6:将所述化学需氧量排放截流量Q(COD)i、所述氨氮排放截流量Q(NH3-N)i和所述总磷排放截流量Q(TP)i分别除以所述分流制区域面积,得到三种污染物排放截流厚度,分别为化学需氧量排放截流厚度H(COD)i、氨氮排放截流厚度H(NH3-N)i和总磷排放截流厚度H(TP)i,单位均为mm;
步骤S7:分别计算所述化学需氧量排放截流厚度H(COD)i、所述氨氮排放截流厚度H(NH3-N)i和所述总磷排放截流厚度H(TP)i的平均值,得到化学需氧量排放截流厚度平均值H(COD)、氨氮排放截流厚度平均值H(NH3-N)和总磷排放截流厚度平均值H(TP)。
2.根据权利要求1所述的一种分流制区域控制降雨径流污染和混接污染的截流厚度确定方法,其特征在于:步骤S7计算化学需氧量排放截流厚度平均值H(COD)、氨氮排放截流厚度平均值H(NH3-N)和总磷排放截流厚度平均值H(TP)的具体计算公式如下:
Figure FDA0003619106840000021
Figure FDA0003619106840000022
Figure FDA0003619106840000023
式中:m为降雨场数。
3.根据权利要求1所述的一种分流制区域控制降雨径流污染和混接污染的截流厚度确定方法,其特征在于:步骤S1中所述的排水管网包括雨水井、雨水排口、雨水管段。
4.根据权利要求3所述的一种分流制区域控制降雨径流污染和混接污染的截流厚度确定方法,其特征在于:步骤S2中所述的建立分流制区域排水管网模型具体为:提供所述排水管网的拓扑结构,将所述雨水管段连接至所述雨水排口,并划分集水区;设置各个集水区的产流表面类型,选择非线性水库模型进行地表汇流演算,得到分流制区域排水管网模型。
5.根据权利要求4所述的一种分流制区域控制降雨径流污染和混接污染的截流厚度确定方法,其特征在于:所述设置各个集水区的产流表面类型具体为:将所述各个集水区的产流表面类型分为两种,分别为非渗透表面和渗透表面两种类型。
CN201910283826.6A 2019-04-10 2019-04-10 分流制区域控制降雨径流和混接污染的截流厚度确定方法 Active CN109977601B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910283826.6A CN109977601B (zh) 2019-04-10 2019-04-10 分流制区域控制降雨径流和混接污染的截流厚度确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910283826.6A CN109977601B (zh) 2019-04-10 2019-04-10 分流制区域控制降雨径流和混接污染的截流厚度确定方法

Publications (2)

Publication Number Publication Date
CN109977601A CN109977601A (zh) 2019-07-05
CN109977601B true CN109977601B (zh) 2022-06-14

Family

ID=67083858

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910283826.6A Active CN109977601B (zh) 2019-04-10 2019-04-10 分流制区域控制降雨径流和混接污染的截流厚度确定方法

Country Status (1)

Country Link
CN (1) CN109977601B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103577640A (zh) * 2013-11-01 2014-02-12 上海市水务规划设计研究院 城市初期雨水截留量标准的制定方法
CN107368623A (zh) * 2017-06-06 2017-11-21 西安建筑科技大学 一种基于径流污染控制的调蓄池容积确定方法
RU2667745C1 (ru) * 2017-08-07 2018-09-24 Государственное Унитарное Предприятие "Водоканал Санкт-Петербурга" Способ оптимизации потоков сточных вод

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103577640A (zh) * 2013-11-01 2014-02-12 上海市水务规划设计研究院 城市初期雨水截留量标准的制定方法
CN107368623A (zh) * 2017-06-06 2017-11-21 西安建筑科技大学 一种基于径流污染控制的调蓄池容积确定方法
RU2667745C1 (ru) * 2017-08-07 2018-09-24 Государственное Унитарное Предприятие "Водоканал Санкт-Петербурга" Способ оптимизации потоков сточных вод

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"宜兴典型村落不同下垫面降雨径流污染物排放特征";汪楚乔 等;《生态与农村环境学报》;20160725;全文 *
"常州市湖塘纺织工业园降雨径流污染负荷分析";王宇翔 等;《水资源保护》;20170520;全文 *

Also Published As

Publication number Publication date
CN109977601A (zh) 2019-07-05

Similar Documents

Publication Publication Date Title
CN111428972B (zh) 基于耦合模型的调蓄工程环境效应评估方法及装置
CN111784052B (zh) 一种城市面源污染入河负荷量预测方法
Imteaz et al. Modelling stormwater treatment systems using MUSIC: Accuracy
CN110334468B (zh) 一种城市排水管网雨水汇入量与溢流量的定量方法
Li et al. Modeling the Quality and Quantity of Runoff in a Highly Urbanized Catchment Using Storm Water Management Model.
CN108984823B (zh) 一种合流制溢流调蓄池规模的确定方法
CN112883644B (zh) 一种动态水环境管理方法
CN114169266B (zh) 一种河流污染源解析和污染物运移计算方法及系统
CN102201034A (zh) 一种临河村落非点源污染负荷的估算方法
CN109598428B (zh) 一种基于行政单元和水系的污染物削减分配方法
CN111444595B (zh) 基于swmm与efdc耦合模型的截污工程环境效应评估方法及装置
CN112766669A (zh) 基于水质模型的流域污染限值排放核定方法
CN113902344A (zh) 一种基于模型的水环境水质目标管控方法
RU2592611C2 (ru) Система оценки баланса подачи и отведения воды мегаполиса
Jung et al. The economic effect of green roofs on non-point pollutant sources management using the replacement cost approach
CN109977601B (zh) 分流制区域控制降雨径流和混接污染的截流厚度确定方法
Manache et al. Calibration of a continuous simulation fecal coliform model based on historical data analysis
Bach et al. Analysis of intensively used catchments based on integrated modelling
CN111723976B (zh) 一种雨水调蓄池空间分布规划方法
Kipyego et al. Analysis of nonpoint source pollution loading on water quality in an urban-rural river catchment using GIS-PLOAD model: Case study of Sosiani River watershed
CN114818211A (zh) 一种基于mike+构建污水收集系统量化模拟技术方法
KR100599848B1 (ko) 비점오염원 모니터링 방법 및 시스템
Kim et al. Analysis of the effect of bio-retention cells to improve water cycle and water quality in urban streams
CN113919625A (zh) 一种雨天溢流污染评价与管控方法
Matar et al. Modelling the benefits of urine source separation scenarios on wastewater treatment plants within an urban water basin

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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 350108 No.1, Gaoxin Avenue, high tech Zone, Minhou County, Fuzhou City, Fujian Province

Applicant after: Fuzhou planning and Design Institute Group Co.,Ltd.

Address before: 350108 No.1, Gaoxin Avenue, high tech Zone, Minhou County, Fuzhou City, Fujian Province

Applicant before: FUZHOU PLANNING DESIGN & Research Institute

GR01 Patent grant
GR01 Patent grant