CN106570334A - 一种地统计学模型计算城市降雨径流污染负荷的方法 - Google Patents
一种地统计学模型计算城市降雨径流污染负荷的方法 Download PDFInfo
- Publication number
- CN106570334A CN106570334A CN201610986571.6A CN201610986571A CN106570334A CN 106570334 A CN106570334 A CN 106570334A CN 201610986571 A CN201610986571 A CN 201610986571A CN 106570334 A CN106570334 A CN 106570334A
- Authority
- CN
- China
- Prior art keywords
- runoff
- class
- pollution
- smc
- model
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种地统计学模型计算城市降雨径流污染负荷的方法。本方法公开了一种基于地理信息技术的计算城市降雨径流污染负荷的方法,它包括降雨径流污染信息采集;空间信息数据库建立;最优GIS模型构建和模型处理;路面(地面)径流污染负荷计算;区域径流污染负荷计算等步骤。本发明将GIS技术引入对降雨径流污染空间分布特征的分析,减少主观方法带来的失误;同时,采用GIS技术对降雨径流污染进行空间模拟和分析,是一种快速、精确的计算城市降雨径流污染负荷的方法,可供制定有效的降雨径流污染控制策略、实施污染控制措施、开展海绵城市建设提供科学指导。
Description
技术领域
本发明涉及一种对城市降雨径流污染的分析方法,尤其涉及一种基于地统计学模型计算城市降雨径流污染负荷的方法。
背景技术
城市降雨径流污染指降雨过程中,降雨冲刷城镇地面,如商业区、街道、停车场等,形成了汇聚原油、盐分、氮、磷、有毒物质等一系列的污染物的径流,随之流入河流或湖泊,造成城市水体污染的现象。
随着我国各城市逐渐加大对工业废水、城市污水等点源污染的控制,城市降雨径流污染(城市非点源污染的主要类型)对城市水环境污染贡献的比重逐渐凸显。国外的城市水环境治理的相关学者认为即使完全控制住城市点源污染,但忽视非点源污染,特别是城市降雨径流污染的话,城市水体的水质仍然不会改善。因此,研究城市降雨径流污染负荷(污染产生量)对于开展城市降雨径流污染控制、制定污染控制措施,从而达到改善城市水环境的目的有着非常重要的意义。
目前,评估城市降雨径流污染负荷常用的方法为输出系数法,该方法起源于1970s北美地区,经过Johnes等学者的不断修正后,目前输出系数法已经在全世界得到了广泛的应用,并且积累了大量的输出的系数值。具体公式如下:
上式中,L是指降雨径流污染物的产生负荷(t/a),Cj是指第i类土地类型(商业用地、居住用地、绿地等)产生的降雨径流污染物的年平均负荷(t/ha/a),Ai是指第i类土地利用的面积(ha)。
应该指出,尽管输出系数法经过不断的修正,具有所需参数少、数据获取简单、结果精度尚可等优点,但其本质仍沿用城镇用地类型的划分方式(国内通常采用《城市用地分类与规划建设用地标准》GB50137),且将用地类型的径流负荷输出量单一化。这种径流污染负荷的计算方法存在如下问题,而导致区域径流污染负荷的计算结果误差较大:
①应根据下垫面类型来计算径流污染负荷,而非城镇用地类型。降雨径流污染程度主要跟下垫面类型有关,而下垫面的划分并不能简单的采用城镇用地类型的划分,而应分为屋面和路面(地面)径流。典型的例子为商业区的径流污染应分为商业区屋顶和商业区路面(地面),不能仅通过路面(地面)径流代表整个商业区域的径流污染情况;另商业区内的绿地、水体不对外输出径流污染,另在计算整个商业区径流污染时剔除。
②应考虑同一下垫面的异质性和相邻下垫面的关联性。同一块用地类型的下垫面并不是径流污染的均质体,其污染程度受人为活动强度和距离的影响;而相邻的不同用地类型的下垫面并不是相互独立的污染输出板块,其因人为活动强度在空间上的关联性而导致污染输出强度亦具有关联性。
本发明从下垫面的分类和径流污染的空间相关性的角度处理上述两个问题,并形成一种快速、精确的计算城市降雨径流污染负荷的方法。
发明内容
本发明的目的是提供一种地统计学模型计算城市降雨径流污染负荷的方法,采用GIS技术对降雨径流污染进行空间模拟和分析,为相对于输出系数法而更为精确的计算降雨径流污染负荷,从而有效制定降雨径流污染控制策略、实施污染控制措施、开展海绵城市建设提供科学指导。
为实现上述目的,本发明采用如下技术方案:
一种地统计学模型计算城市降雨径流污染负荷的方法,其特征在于它包括以下具体步骤:
A、降雨径流污染信息采集,对地表类不透水面和屋面类不透水面分别进行降雨采样,并获取采样点的GPS定位数据;;
B、建立地表类不透水面的降雨径流污染属性数据库,通过采样点GPS定位数据录入基础地理信息数据库,将地面类不透水面采样点的SMC作为属性数据与相应的空间坐标进行关联,建立地表类不透水面的降雨径流污染属性数据库;
C、最优GIS模型构建,进行空间模拟,对各地面类不透水面采样点的地理位置和SMC进行拟合分析,通过拟合分析选择最优半方差模型以及最优参数,并利用该最优半方差模型以及最优参数对各地面类不透水面采样点之间的位置进行SMC的插值计算;
D、计算地面径流污染负荷LRd,输出根据步骤C得到的降雨径流污染插值图,将插值图栅格化后,根据栅格化的降雨径流污染插值图计算地面径流污染负荷LRd;
E、计算屋面类不透水面的径流污染负荷LRf,将屋面类不透水面的降雨径流污染SMC进行加权平均,作为屋面类不透水面的统一SMC,根据所述统一SMC计算屋面类不透水面的径流污染负荷LRf;
F、计算区域径流污染负荷L,L=LRd+LRf。
优选地,步骤A中,降雨采样分为地表类不透水面和屋面类不透水下垫面的采样,采用人工采样或降雨径流采样器进行降雨采样,采样过程中,前30min内径流采样间隔时间宜为5min,30min~90min内采样间隔时间宜为10min,90min之后的采样间隔时间可适当扩大,采样监测指标为COD、BOD、NH4 +-N、TN、TP、重金属和POPs。
优选地,步骤B包括以下步骤:1)得到获取空间信息数据,空间解译数据;2)删除空间解译数据中的透水下垫面或为自然水体的区域;3)将不透水下垫面进一步细分为地表类不透水下垫面和屋面类不透水下垫面;4)通过采样点GPS定位数据录入基础地理信息数据库,将地面类不透水面采样点的SMC作为属性数据与相应的空间坐标进行关联,建立地表类不透水面的降雨径流污染属性数据库,根据监测数据,加权计算点位降雨径流平均浓度SMC,具体公式如下:
上述公式(1)中,EMC是指事件平均浓度,单位mg/L;M是指降雨事件径流污染物总量,单位mg;V是指降雨事件径流总量,单位L;Qj是指第j次取样时所对应的径流量,单位L/min;Cj指j次取样时的径流污染物浓度,单位mg/L;上述公式(2)中,SMC是指点位降雨径流平均浓度,通过加权计算多场降雨的次降雨径流平均浓度求得,P是指第i次降雨量,单位L。
优选地,步骤C包括以下步骤:
1)采用GS+(Version 9,Gamma Design Software company)拟合半方差函数并优化其参数。将对路面(地面)径流采样点的地理信息和污染因子SMC输入GS+软件中;用GS+对污染因子SMC进行空间模拟,通过拟合分析,选择最优半方差模型以及最优参数;
2)将GS+拟合的半方差函数最优参数输入ArcGIS(Version 9.3,ESRI.company)中的地统计模块(Geostatistical analyst),采用最常用的普通克里金法(OrdinaryKriging)进行最优化无差的空间插值,利用配对t检验(Paired t test)对空间差值结果与实测结果的进行精度检验,最终确定最优GIS分析模型。
所述步骤C中最优半方差模型线性包括基台模型、球状模型、指数模型和高斯模型,其空间模拟的最优标准为决定系数(R2)值最大和残差平方和(RSS)值最小。
优选地,步骤D又包括以下步骤:
1)输出降雨径流污染插值图,将插值图栅格化后,并对插值图进行提取边界和裁剪的优化措施,剔除研究范围内的下凹式绿地和水体的范围;
2)点击插值图层的属性,输出插值图层的属性表;
3)通过如下公式,计算路面(地面)径流污染负荷:
式中,LRd——地表类不透水面年径流污染负荷,单位kg;
D——年均降雨量,单位mm/a;
R——研究区域的年径流系数;
S——插值单元面积,单位m2;
SMCi——第i个插值单元的点位降雨径流平均浓度,单位mg/L。
优选地,步骤E中,计算屋面类不透水面的径流污染负荷LRf的公式为:
LRf=0.001D×R×SMCRf×SRf
式中,LRf——屋面类不透水面年径流污染负荷,单位kg;
R——研究区域的年径流系数;
S——屋面类不透水面的面积,单位m2;
SMCRf——屋面类不透水面的点位平均浓度,单位mg/L。
本发明将GIS技术引入对降雨径流污染空间分布和污染负荷计算,可减少主观方法带来的失误;同时,这是一种快速、精确的计算城市降雨径流污染负荷的方法,可供制定有效的降雨径流污染控制策略、实施污染控制措施提供科学指导。
附图说明
图1为本发明的流程图。
图2为本发明一实施例的监测点位总图。
图3为该实施例EPASS的最优半方差函数模型及拟合结果。
图4为利用普通克立格方法插值获得的地表径流SMCSS分布图。
具体实施方法
一种地统计学模型计算城市降雨径流污染负荷的方法,其特征在于它包括以下具体步骤:
A、降雨径流污染信息采集
在现场调研研究区域用地类型、用地面积等用地信息的基础上,制定降雨径流采样方案,包括采样点位(包括屋面类采样点和地表类采样点)、采样次数、监测指标等。调研研究区域时需现场获取采样点的GPS坐标。
采用人工采样或自动径流采样器进行降雨采样。采样过程中,前30min内径流采样间隔时间宜为5min,30min~90min内采样间隔时间宜为10min,90min之后的采样间隔时间可适当扩大。采样后应尽快进行试验分析,常见分析指标为COD、BOD、NH4 +-N和TP,以及重金属和POPs等。
B、空间信息数据库建立
①获取空间信息数据。可通过研究区域的遥感影像资料解译用地类型和向相关部门索取研究区域用地类型的解译数据。②删除空间解译数据中的绿地、湿地、林地、湖泊等透水下垫面或为自然水体的区域。③将不透水下垫面进一步细分为地表类不透水下垫面(RD)和屋面类不透水下垫面(RF)。④进行各污染物的相关性分析,通过采样点GPS定位数据录入基础地理信息数据库,将地面类不透水面采样点的SMC作为属性数据与相应的空间坐标进行关联,建立地表类不透水面的降雨径流污染属性数据库。
根据监测数据,加权计算SMC(点位降雨径流平均浓度),具体公式如下:
上述公式(1)中,EMC是指事件平均浓度,即次降雨径流平均浓度,mg/L;M是指降雨事件径流污染物总量,mg;V是指降雨事件径流总量,L;Qj是指第j次取样时所对应的径流量,L/min;Cj指j次取样时的径流污染物浓度,mg/L。上述公式(2)中,SMC是指点位降雨径流平均浓度,通过加权计算多场降雨的次降雨径流平均浓度求得,P是指第i次降雨量,L。
C、最优GIS模型构建。
1)采用GS+(Version 9,Gamma Design Software company)拟合半方差函数并优化其参数。将对路面(地面)径流采样点的地理信息和污染因子SMC输入GS+软件中;用GS+对污染因子SMC进行空间模拟,通过拟合分析,选择最优半方差模型以及最优参数;
2)将GS+拟合的半方差函数最优参数输入ArcGIS(Version 9.3,ESRI.company)中的地统计模块(Geostatistical analyst),采用最常用的普通克里金法(OrdinaryKriging)进行最优化无差的空间插值,利用配对t检验(Paired t test)对空间差值结果与实测结果的进行精度检验,最终确定最优GIS分析模型。
D、路面(地面)径流污染负荷计算。
1)输出降雨径流污染插值图,将插值图栅格化后,并对插值图进行提取边界(Extract by mask)和裁剪等优化措施,剔除研究范围内的下凹式绿地和水体的范围;2)点击插值图层的属性(Properties),输出插值图层的属性表。3)通过如下公式,计算路面(地面)径流污染负荷:
式中,LRd——地表类不透水面年径流污染负荷,kg;
D——年均降雨量,mm/a;
R——研究区域的年径流系数;
S——插值单元面积,m2;
SMCi——第i个插值单元的点位降雨径流平均浓度,mg/L;
E、屋面类不透水面的径流污染负荷计算
LRf=0.001D×R×SMCRf×SRf
式中,LRf——屋面类不透水面年径流污染负荷,kg;
R——研究区域的年径流系数;
S——屋面类不透水面的面积,m2;
SMCRf——屋面类不透水面的点位平均浓度,mg/L;
F、区域径流污染负荷计算
L=LRd+LRf
式中,L——研究区域的不透水面年径流污染负荷,kg。
所述步骤A中降雨径流污染信息采集的步骤包括:在现场调研研究区域用地类型、用地面积等用地信息的基础上,制定降雨径流采样方案,包括采样点位、采样次数、监测指标等。调研研究区域时需现场获取采样点的GPS坐标。
所述步骤A中采用人工采样或自动径流采样器进行降雨采样。
所述步骤A中采样点位设计应对研究区域进行网格化均匀布点,并在可接受采样强度范围内对特殊区域(如大片的广场等)进行密集布点,以达到分析其降雨径流污染的空间差异性的目的。
所述步骤A中采样次数设置应满足如下要求:需对研究范围的降雨径流进行长期采样监测;时间跨度至少一年,并保证每季度至少采样监测一次;降雨量为大雨、中雨和小雨的降雨各一次。
所述步骤A采样点包括路面(地面)采样点和屋面采样点。
所述步骤B中解译依据为《城市用地分类与规划建设用地标准》GB50137。
所属步骤B中删除上述用地类型区域的方法为通过现场调研和空间数据分析相结合的方式,以下垫面的渗透性和高程低于周边硬化下垫面为依据进行空间数据预处理。
所述步骤B中降雨径流污染检测数据预处理和数据库建立的步骤包括:根据监测数据,进行描述性统计,统计内容包括相关性分析、方差分析、正态分布检验等,以分析降雨径流的基本特征,并便于选取典型污染物的SMC进行空间模拟;通过计算不同采样点的单次降雨径流污染物的EMC(次降雨径流平均浓度),并进一步加权计算SMC(点位降雨径流平均浓度)。
所述步骤B中GPS定位数据通过采样点录入基础地理信息数据库,将SMC作为属性数据与相应的空间坐标进行关联,建立降雨径流污染属性数据库。
所述步骤B中监测数据的正态分布检验为必须环节,常见检验方法包括Kolmogorov-Smirnov检验和Shapiro-Wilk检验,原始数据不符合正态分布的需进行对数转换、倒数转换等数据处理。
所述步骤C中最优半方差模型线性包括基台模型、球状模型、指数模型和高斯模型,其空间模拟的最优标准为决定系数(R2)值最大和残差平方和(RSS)值最小。
所述步骤D中插值图栅格化的栅格不应大于1m×1m。
为了更好的说明本发明专利的具体实施方案,现以本发明所采用的方法在上海临港新城的应用为例,对本发明内容的实施过程进行实例说明。
上海临港新城位于上海东南角,其陆地区域45%的面积由填海造陆而成,土壤盐碱度较高,移植的植株生长受限临港新城主城区的中心是滴水湖,它是世界上最大的人工景观湖,面积为5.56km2,与沿湖以环状的“三涟”(内涟河,中涟河,外涟河)和“七射”构成了临港水系的主构架。临港主城区划分为滴水湖二级集水区,主城区采用的雨污分流的市政排水系统,暴雨时期,溢出的径流污水将排入临港水系,进入滴水湖对其景观水质造成影响。
(1)降雨径流污染信息采集。
在现场调研研究区域用地类型、用地面积等用地信息的基础上,制定降雨径流采样方案。采样点位如下图下表所示。
表1监测点位信息
采用自动径流采样器共采集有效降雨事件5次,具体见表2。
表2采集降雨事件
分析指标为COD、BOD、NH4 +-N和TP。
表3各污染物描述性统计
(2)空间信息数据库建立。
上海临港新城2008年土地利用数据解译于华东师范大学地理信息科学教育部重点实验室,根据其人为活动的空间关联性,本研究将不透水面分为道路广场类不透水面与屋面类不透水面。
通过采样点GPS定位数据录入基础地理信息数据库,将地面类不透水面采样点的SMC作为属性数据与相应的空间坐标进行关联,建立地表类不透水面的降雨径流污染属性数据库。
根据监测数据,加权计算SS污染的SMC(点位降雨径流平均浓度)。污染物SS的EMC值与CODCr、TP和NH3-N均呈显著相关,在此本研究仅选SS的SMC讨论临港新城地表径流污染的空间差异性。
并通过采样点GPS定位数据录入基础地理信息数据库,将SMCss作为属性数据与相应的空间坐标进行关联,建立降雨径流污染属性数据库。
(3)最优GIS模型构建。
根据公式计算不同监测点位的SMCss,并以SMCss为区域化变量,利用半方差函数,对临港新城道路交通类EIA的SMCss的半方差值随步长变化的理论模型进行拟合,以残差平方和(RSS)最小和决定系数(R2)最大为最优模型的筛选标准,最终确定SMCss在空间模拟的最优模型为高斯模型。
表4临港径流SMCSS的全方位的半方差模型及其参数
(4)路面(地面)径流污染负荷计算。
A.输出降雨径流污染插值图,将插值图栅格化后,并对插值图进行提取边界(Extract by mask)和裁剪等优化措施,剔除研究范围内的下凹式绿地和水体的范围;B.点击插值图层的属性(Properties),输出插值图层的属性表。C.按照下式,计算路面(地面)径流污染负荷:
式中,LRd——地表类不透水面年径流污染负荷,kg;
D——年均降雨量,mm/a,根据上海市水资源公报,2006-2013年上海市年均降雨量为1198.7mm/a;
R——研究区域的年径流系数,鉴于道路、广场等路面(地面)以水泥与大理石为主,此处选取0.95;
S——插值单元面积,m2;
SMCi——第i个插值单元的点位降雨径流平均浓度,mg/L;
(5)屋面类不透水面的径流污染负荷计算
屋面类EIA的计算方法为输出系数法,有别于传统经验系数模型,本研究采用地统计学模型模拟区域的道路广场类EIA的单位面积输出系数来计算区域的非点源污染负荷,具体如下式:
LRf=0.001D×R×SMCRf×SRf
式中,LRf——屋面类不透水面年径流污染负荷,kg;
R——研究区域的年径流系数,鉴于屋面以水泥与大理石为主,此处选取0.95;
S——屋面类不透水面的面积,m2;
SMCRf——屋面类不透水面的点位平均浓度,mg/L;
(6)区域径流污染负荷计算
L=LRd+LRf
式中,L——研究区域的不透水面年径流污染负荷,kg。
结果如下:
表5临港新城主城区非点源SS的产生负荷
*括号内为均值。
Claims (7)
1.一种地统计学模型计算城市降雨径流污染负荷的方法,其特征在于它包括以下具体步骤:
A、降雨径流污染信息采集,对地表类不透水面和屋面类不透水面分别进行降雨采样,并获取采样点的GPS定位数据;;
B、建立地表类不透水面的降雨径流污染属性数据库,通过采样点GPS定位数据录入基础地理信息数据库,将地面类不透水面采样点的SMC作为属性数据与相应的空间坐标进行关联,建立地表类不透水面的降雨径流污染属性数据库;
C、最优GIS模型构建,进行空间模拟,对各地面类不透水面采样点的地理位置和SMC进行拟合分析,通过拟合分析选择最优半方差模型以及最优参数,并利用该最优半方差模型以及最优参数对各地面类不透水面采样点之间的位置进行SMC的插值计算;
D、计算地面径流污染负荷LRd,输出根据步骤C得到的降雨径流污染插值图,将插值图栅格化后,根据栅格化的降雨径流污染插值图计算地面径流污染负荷LRd;
E、计算屋面类不透水面的径流污染负荷LRf,将屋面类不透水面的降雨径流污染SMC进行加权平均,作为屋面类不透水面的统一SMC,根据所述统一SMC计算屋面类不透水面的径流污染负荷LRf;
F、计算区域径流污染负荷L,L=LRd+LRf。
2.如权利要求1所述的一种地统计学模型计算城市降雨径流污染负荷的方法,其特征在于步骤A中,降雨采样分为地表类不透水面和屋面类不透水下垫面的采样,采用人工采样或降雨径流采样器进行降雨采样,采样过程中,前30min内径流采样间隔时间宜为5min,30min~90min内采样间隔时间宜为10min,90min之后的采样间隔时间可适当扩大,采样监测指标为COD、BOD、NH4 +-N、TN、TP、重金属和POPs。
3.如权利要求1所述的一种地统计学模型计算城市降雨径流污染负荷的方法,其特征在于步骤B包括以下步骤:1)得到获取空间信息数据,空间解译数据;2)删除空间解译数据中的透水下垫面或为自然水体的区域;3)将不透水下垫面进一步细分为地表类不透水下垫面和屋面类不透水下垫面;4)通过采样点GPS定位数据录入基础地理信息数据库,将地面类不透水面采样点的SMC作为属性数据与相应的空间坐标进行关联,建立地表类不透水面的降雨径流污染属性数据库;5)根据监测数据,加权计算点位降雨径流平均浓度SMC,具体公式如下:
上述公式(1)中,EMC是指事件平均浓度,单位mg/L;M是指降雨事件径流污染物总量,单位mg;V是指降雨事件径流总量,单位L;Qj是指第j次取样时所对应的径流量,单位L/min;Cj指j次取样时的径流污染物浓度,单位mg/L;上述公式(2)中,SMC是指点位降雨径流平均浓度,通过加权计算多场降雨的次降雨径流平均浓度求得,P是指第i次降雨量,单位L。
4.如权利要求1所述的一种地统计学模型计算城市降雨径流污染负荷的方法,其特征在于步骤C包括以下步骤:
1)采用GS+(Version 9,Gamma Design Software company)拟合半方差函数并优化其参数,将对路面(地面)径流采样点的地理信息和污染因子SMC输入GS+软件中;用GS+对污染因子SMC进行空间模拟,通过拟合分析,选择最优半方差模型以及最优参数;
2)将GS+拟合的半方差函数最优参数输入ArcGIS(Version 9.3,ESRI.company)中的地统计模块(Geostatistical analyst),采用最常用的普通克里金法(Ordinary Kriging)进行最优化无差的空间插值,利用配对t检验(Paired t test)对空间差值结果与实测结果的进行精度检验,最终确定最优GIS分析模型。
5.如权利要求1所述的一种地统计学模型计算城市降雨径流污染负荷的方法,其特征在于所述步骤C中最优半方差模型线性包括基台模型、球状模型、指数模型和高斯模型,其空间模拟的最优标准为决定系数值最大和残差平方和值最小。
6.如权利要求1所述的一种地统计学模型计算城市降雨径流污染负荷的方法,其特征在于步骤D又包括以下步骤:
1)输出降雨径流污染插值图,将插值图栅格化后,并对插值图进行提取边界和裁剪的优化措施,剔除研究范围内的下凹式绿地和水体的范围;
2)点击插值图层的属性,输出插值图层的属性表;
3)通过如下公式,计算路面(地面)径流污染负荷:
式中,LRd——地表类不透水面年径流污染负荷,单位kg;
D——年均降雨量,单位mm/a;
R——研究区域的年径流系数;
S——插值单元面积,单位m2;
SMCi——第i个插值单元的点位降雨径流平均浓度,单位mg/L。
7.如权利要求1所述的一种地统计学模型计算城市降雨径流污染负荷的方法,
其特征在于步骤E中,计算屋面类不透水面的径流污染负荷LRf的公式为:
LRf=0.001D×R×SMCRf×SRf
式中,LRf——屋面类不透水面年径流污染负荷,单位kg;
R——研究区域的年径流系数;
S——屋面类不透水面的面积,单位m2;
SMCRf——屋面类不透水面的点位平均浓度,单位mg/L。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610986571.6A CN106570334B (zh) | 2016-11-09 | 2016-11-09 | 一种地统计学模型计算城市降雨径流污染负荷的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610986571.6A CN106570334B (zh) | 2016-11-09 | 2016-11-09 | 一种地统计学模型计算城市降雨径流污染负荷的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106570334A true CN106570334A (zh) | 2017-04-19 |
CN106570334B CN106570334B (zh) | 2021-12-28 |
Family
ID=58540633
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610986571.6A Active CN106570334B (zh) | 2016-11-09 | 2016-11-09 | 一种地统计学模型计算城市降雨径流污染负荷的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106570334B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108062454A (zh) * | 2018-01-19 | 2018-05-22 | 宁波市镇海规划勘测设计研究院 | 污染物时空分布不确定性特征分析方法、系统及存储介质 |
CN110675091A (zh) * | 2019-10-16 | 2020-01-10 | 中国水利水电科学研究院 | 一种城市地表径流污染物负荷计算方法 |
CN110849415A (zh) * | 2019-11-07 | 2020-02-28 | 中交天航港湾建设工程有限公司 | 一种降雨点位布设及监测系统及方法 |
CN111126702A (zh) * | 2019-12-24 | 2020-05-08 | 武汉新烽光电股份有限公司 | 道路面源污染信息预测方法、装置、设备及存储介质 |
CN111191905A (zh) * | 2019-12-24 | 2020-05-22 | 上海市政工程设计研究总院(集团)有限公司 | 一种基于三级管控体系的城市降雨径流污染控制方法 |
CN111310124A (zh) * | 2020-02-13 | 2020-06-19 | 北京市环境保护科学研究院 | 一种城市径流数据处理方法 |
CN111552763A (zh) * | 2020-05-13 | 2020-08-18 | 水利部交通运输部国家能源局南京水利科学研究院 | 城市面源污染负荷监测方法 |
CN113689099A (zh) * | 2021-08-16 | 2021-11-23 | 河海大学 | 基于水量构成的环状河网水系水质考核断面污染负荷贡献率确定方法 |
CN115687855A (zh) * | 2022-09-26 | 2023-02-03 | 中国科学院生态环境研究中心 | 基于粒径与下垫面特征对传统指数冲刷方程的改进方法 |
CN115712813A (zh) * | 2022-11-23 | 2023-02-24 | 北京中科三清环境技术有限公司 | 确定非点源的污染负荷输出系数的方法、装置及电子设备 |
Citations (6)
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 |
US20060262963A1 (en) * | 2005-05-23 | 2006-11-23 | Digitalglobe, Inc. | Method and Apparatus for Determination of Water Pervious Surfaces |
CN102201034A (zh) * | 2011-04-26 | 2011-09-28 | 北京师范大学 | 一种临河村落非点源污染负荷的估算方法 |
CN102841081A (zh) * | 2012-08-30 | 2012-12-26 | 湖南科技大学 | 有色金属矿区土-水界面污染流中各重金属分布预测方法 |
CN103020424A (zh) * | 2012-11-22 | 2013-04-03 | 北京师范大学 | 一种基于人工降雨实验的北方平原农田区非点源污染负荷估算方法 |
CN103577640A (zh) * | 2013-11-01 | 2014-02-12 | 上海市水务规划设计研究院 | 城市初期雨水截留量标准的制定方法 |
-
2016
- 2016-11-09 CN CN201610986571.6A patent/CN106570334B/zh active Active
Patent Citations (6)
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 |
US20060262963A1 (en) * | 2005-05-23 | 2006-11-23 | Digitalglobe, Inc. | Method and Apparatus for Determination of Water Pervious Surfaces |
CN102201034A (zh) * | 2011-04-26 | 2011-09-28 | 北京师范大学 | 一种临河村落非点源污染负荷的估算方法 |
CN102841081A (zh) * | 2012-08-30 | 2012-12-26 | 湖南科技大学 | 有色金属矿区土-水界面污染流中各重金属分布预测方法 |
CN103020424A (zh) * | 2012-11-22 | 2013-04-03 | 北京师范大学 | 一种基于人工降雨实验的北方平原农田区非点源污染负荷估算方法 |
CN103577640A (zh) * | 2013-11-01 | 2014-02-12 | 上海市水务规划设计研究院 | 城市初期雨水截留量标准的制定方法 |
Non-Patent Citations (4)
Title |
---|
吕永鹏: "平原河网地区城市集水区非点源污染过程模拟与系统调控管理研究", 《中国优秀博硕士学位论文全文数据库(博士) 工程科技I辑》 * |
李春林 等: "沈阳市降雨径流污染物排放特征", 《生态学杂志》 * |
陈莹等: "城市路面径流污染负荷及采样次数对其估算准确性的影响", 《城市环境与城市生态》 * |
黎巍等: "滇池流域城市降雨径流污染负荷定量化研究", 《环境监测管理与技术》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108062454A (zh) * | 2018-01-19 | 2018-05-22 | 宁波市镇海规划勘测设计研究院 | 污染物时空分布不确定性特征分析方法、系统及存储介质 |
CN110675091B (zh) * | 2019-10-16 | 2022-08-09 | 中国水利水电科学研究院 | 一种城市地表径流污染物负荷计算方法 |
CN110675091A (zh) * | 2019-10-16 | 2020-01-10 | 中国水利水电科学研究院 | 一种城市地表径流污染物负荷计算方法 |
CN110849415A (zh) * | 2019-11-07 | 2020-02-28 | 中交天航港湾建设工程有限公司 | 一种降雨点位布设及监测系统及方法 |
CN111126702A (zh) * | 2019-12-24 | 2020-05-08 | 武汉新烽光电股份有限公司 | 道路面源污染信息预测方法、装置、设备及存储介质 |
CN111191905A (zh) * | 2019-12-24 | 2020-05-22 | 上海市政工程设计研究总院(集团)有限公司 | 一种基于三级管控体系的城市降雨径流污染控制方法 |
CN111126702B (zh) * | 2019-12-24 | 2024-02-27 | 武汉新烽光电股份有限公司 | 道路面源污染信息预测方法、装置、设备及存储介质 |
CN111310124A (zh) * | 2020-02-13 | 2020-06-19 | 北京市环境保护科学研究院 | 一种城市径流数据处理方法 |
CN111552763A (zh) * | 2020-05-13 | 2020-08-18 | 水利部交通运输部国家能源局南京水利科学研究院 | 城市面源污染负荷监测方法 |
CN113689099A (zh) * | 2021-08-16 | 2021-11-23 | 河海大学 | 基于水量构成的环状河网水系水质考核断面污染负荷贡献率确定方法 |
CN113689099B (zh) * | 2021-08-16 | 2022-05-06 | 河海大学 | 基于水量构成的环状河网水系水质考核断面污染负荷贡献率确定方法 |
CN115687855A (zh) * | 2022-09-26 | 2023-02-03 | 中国科学院生态环境研究中心 | 基于粒径与下垫面特征对传统指数冲刷方程的改进方法 |
CN115687855B (zh) * | 2022-09-26 | 2023-04-07 | 中国科学院生态环境研究中心 | 基于粒径与下垫面特征对传统指数冲刷方程的改进方法 |
CN115712813A (zh) * | 2022-11-23 | 2023-02-24 | 北京中科三清环境技术有限公司 | 确定非点源的污染负荷输出系数的方法、装置及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN106570334B (zh) | 2021-12-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106570334A (zh) | 一种地统计学模型计算城市降雨径流污染负荷的方法 | |
Shen et al. | Relating landscape characteristics to non-point source pollution in a typical urbanized watershed in the municipality of Beijing | |
Shi et al. | Modelling hydrology and water quality processes in the Pengxi River basin of the Three Gorges Reservoir using the soil and water assessment tool | |
Demesouka et al. | Suitability analysis for siting MSW landfills and its multicriteria spatial decision support system: method, implementation and case study | |
CN109740562B (zh) | 一种生态海绵型城市建设适宜区的系统及方法 | |
Pisciotta et al. | Groundwater nitrate risk assessment using intrinsic vulnerability methods: A comparative study of environmental impact by intensive farming in the Mediterranean region of Sicily, Italy | |
US20180058932A1 (en) | Method for analyzing the types of water sources based on natural geographical features | |
Tong et al. | Modeling the relationship between land use and surface water quality | |
CN107563139A (zh) | 一种点源突发性水污染源事故溯源贡献程度计算方法 | |
CN103544550B (zh) | 一种金属矿区土-水界面重金属污染负荷的预测方法 | |
Gardes et al. | Statistical prediction of the nocturnal urban heat island intensity based on urban morphology and geographical factors-An investigation based on numerical model results for a large ensemble of French cities | |
CN106528788A (zh) | 基于gis技术分析地面雨径流污染空间分布特征的方法 | |
Capri et al. | Assessment of nitrate contamination risk: The Italian experience | |
CN104361523A (zh) | 一种基于gis的分布式稻田氮素径流流失负荷估算方法 | |
CN103294902A (zh) | 一种基于遥感影像与gis空间分析确定自然湿地恢复方案的方法 | |
Hosseini | Effect of land use change on water balance and suspended sediment yield of Taleghan catchment, Iran | |
Zhou et al. | Establishment and validation of an amended phosphorus index: refined phosphorus loss assessment of an agriculture watershed in Northern China | |
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 | |
Richardson et al. | Assessing watershed vulnerability in Bernalillo County, New Mexico using GIS-based fuzzy Inference | |
Bai et al. | Simulating the impacts of land-use changes on non-point source pollution in Lugu Lake watershed | |
Fu et al. | GIS and ANN-based spatial prediction of DOC in river networks: a case study in Dongjiang, Southern China | |
Walters et al. | Examining the relationship between development patterns and total phosphorus in the Galveston Bay Estuary | |
Bellamy et al. | A GIS-based approach for determining potential runoff coefficient and runoff depth for the Indian River Lagoon, Florida, USA | |
Ruibal-Conti et al. | Performance assessment of the Anvil Way Compensation Basin living stream: 2004-2013 | |
Chen et al. | An integrated simulation-monitoring framework for nitrogen assessment: a case study in the Baixi watershed, China |
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 |