CN108733915A - 基于栅格化处理的平原河网区河道面源污染负荷确定方法 - Google Patents
基于栅格化处理的平原河网区河道面源污染负荷确定方法 Download PDFInfo
- Publication number
- CN108733915A CN108733915A CN201810474581.0A CN201810474581A CN108733915A CN 108733915 A CN108733915 A CN 108733915A CN 201810474581 A CN201810474581 A CN 201810474581A CN 108733915 A CN108733915 A CN 108733915A
- Authority
- CN
- China
- Prior art keywords
- river
- pollution
- area
- formula
- grid
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于栅格化处理的平原河网区河道面源污染负荷确定方法,采用栅格化处理技术和GIS空间叠合分析功能,以河道断面面积与栅格到周边河道最小距离之比为权重,计算汇流权重因子,确定各河道的汇流区范围及对应的各类土地利用面积。具体步骤包括:绘制河网多边形、计算域栅格化及空间运算、计算汇流权重因子、确定面源污染流向、划分面源污染汇流区、统计土地利用类型面积、计算单位面积污染负荷、计算河道面源污染负荷。本发明通过提出反映平原河网区水系特征的面源污染汇流区确定方法,解决了该类地区汇流区界限难以确定的问题,易于在具有同类水系特征地区的污染负荷计算中推广应用。
Description
技术领域
本发明涉及基于栅格化处理的平原河网区河道面源污染负荷确定方法,属于环境领域。
背景技术
污染物按排放方式不同,可分为点源和面源两大类。污染负荷计算是开展流域水环境治理的重要基础性工作。然而,由于面源污染发生的随机性、机理过程的复杂性、排放途径的不确定性以及时空分布的差异性,造成难以在较大的时空尺度上通过调查和监测等手段获取其排放量,因此,通常采用数学模型方法估算面源污染负荷。
目前国际上比较成熟的面源污染模型主要有GWLF、AGNPS、SWAT、HSPF等。这些模型多依据地形特征来进行水系和流域划分,适用于流域界限清晰的山地和丘陵地区。然而,我国长江、淮河、珠江流域的下游平原地区分布有大量河网,以太湖流域为例,该流域总面积3.69万km2,其中高程10m以下的平原区面积约占80%,地面高程差别小,并且该地区河道纵横交错,形成独特的网状结构。特殊的地貌和水系特征决定了平原河网区的面源汇流区界限难以划分,现有面源污染模型无法直接应用于该地区。
发明内容
发明目的:为了克服现有技术中存在的不足,本发明提供一种基于栅格化处理的平原河网区河道面源污染负荷确定方法,充分利用GIS的栅格化技术和空间分析功能,准确划分面源污染汇流区,易于在平原河网区面源污染负荷计算中推广使用。
技术方案:为解决上述技术问题,本发明的一种基于栅格化处理的平原河网区河道面源污染负荷确定方法,包括以下步骤:
(1)绘制河网多边形结构图
根据计算区域的水系图,对河网进行概化处理,绘制绘制由河道包围所形成的网状多边形结构图;
(2)计算域栅格化及空间运算
借助工具,对计算区域的土地利用图层和栅格图层进行空间叠合分析,生成每个网格单元不同土地利用类型的范围,统计相应土地利用的面积为河网多边形第i个网格第j种土地利用的面积,km2;
(3)计算汇流权重因子
按公式(1)进行计算:
式中:为第i个网格到第k条河道上的权重因子,%;为第i个网格到第k条河道的最小距离,km;Ak为第k条河道的断面面积,m2;m为包围该网格的河道数量;
(4)确定面源污染流向
统计每个网格到第k条河道的最大权重因子,按公式(2)计算:
Mi=max{Pi k,k=1…m} (2)
式中:Mi为某个河网多边形的第i个网格汇流入第k条河道的最大权重因子,%;
如果第i个网格到第k条河道的权重因子最大,那么该网格的面源污染负荷全部汇入第k条河道;
(5)面源污染汇流区划分
依次计算河网多边形内每个网格的面源污染汇入周边河道的最大权重因子,所有汇入同一条河道的网格即构成该河道的汇流区;
(6)统计汇流区土地利用类型面积
按公式(3)统计河网多边形不同土地利用类型的面积。
式中:Aj分别为河网多边形第j种土地利用类型的面积,km2;n为某一河道汇流区的栅格数量;依次对所有河网多边形重复上述计算过程,得到计算区域全部河道的汇流区范围及对应的各种土地利用类型面积。
(7)将河道汇流区的各类土地利用面积输入面源污染模型,计算与每条河道对应的面源污染量WLk;
(8)根据面源污染量WLk,结合工业废水和生活污水的点源调查资料,统计每条河道各种污染物的污染负荷总量WLs,将该数据与政府公布的水污染物限制排污总量WLL进行比较:
如果WLs≤WLL,说明河道的实际污染负荷小于或等于限制排污总量,无需对该河道制定相应减排措施。
如果WLs>WLL,说明该河道的实际污染负荷已经超过限制排污总量,需要制定相应的污染治理和减排措施。对于实际污染负荷超过限制排污总量的河道,计算相应的污染物削减量。根据面源和点源污染负荷构成,明确水污染治理方向,制定分阶段的污染源削减方案。
作为优选,所述步骤(2)中的工具为GIS软件平台的“Intersect”工具。
作为优选,所述步骤(7)中面源污染类型包含城市径流污染、旱地径流污染、稻田径流污染。
作为优选,城市降雨径流污染单元负荷计算包含以下两个步骤:
a、计算城市单位面积的地表污染物累积通量:
Ws=αiFiγiRcl/0.9 (4)
式中:Ws为城市单位面积污染物累积通量,kg/(km2·d);αi为城市污染物浓度参数,mg/L;γi为地面清扫频率参数;Rcl为地表污染冲刷降水量,mm/d;Fi为人口密度参数。
其中,γi=Ni/20(清扫间隔Ni<20h)
γi=1(清扫间隔Ni≥20h)
b、降雨径流冲刷模型
城镇降雨径流的流失速率按公式(5)计:
Wu=Ws(1-e-kRt) (5)
式中:Wu为单位面积城市污染负荷,kg/(km2·d);k为降雨径流对地表污染物的冲刷系数,1/mm,城市地区取0.14~0.19;R为降雨强度,mm/h;t为降雨历时,h。
作为优选,旱地径流污染单元负荷利用式(6)计算而得
Wf=mfη+W0 (6)
式中:η为肥料流失率,%;Wf为某一施肥水平下单位面积肥料流失量,kg/(km2·d);W0为零施肥条件下单位面积肥料年流失量,kg/(km2·d);mf为单位面积施肥量,kg/(km2·d)。
作为优选,稻田径流污染单元负荷的计算包含以下步骤,
①稻田径流氮素流失模型
田面水中TN和NH3-N浓度变化过程按公式(7)和(8)计算:
式中:和为前一时刻和后一时刻的田面水深度,mm;和为前一时刻和后一时刻田面水NH3-N浓度,mg·L-1;和为前一时刻和后一时刻田面水TN浓度,mg·L-1;Ri为稻田灌溉速率,mm·d-1;Ci1和Ci2为稻田灌溉水NH3-N和TN浓度,mg·L-1;Rr,Rd,Rl分别为降水强度、实际排水速率及渗漏速率,mm·d-1;Cr1和Cr2为降水中NH3-N和TN浓度,mg·L-1;Φn为氮肥向田面水的释放通量,kg·hm-2·d-1;kv为溶液中NH3-N的挥发速率常数,d-1;kn和kdn为水土界面的硝化和反硝化速率常数,d-1;
②稻田径流磷素流失模型
田面水磷素浓度变化过程按公式(9)计算:
式中:和为前一时刻和后一时刻田面水TP的质量浓度,mg·L-1;Ri为灌溉速率,mm·d-1;Ci3为灌溉水中TP的质量浓度,mg·L-1;Cr3为降水中TP的质量浓度,mg·L-1;ka为土壤对TP的吸附速率常数,d-1;Φp为磷肥向田面水的释放通量kg·hm-2·d-1;
③单位面积稻田污染物流失量
根据(7)~(9)式计算田面水污染物浓度随时间的变化过程后,根据单位面积稻田排水量按式(10)计算随污染负荷:
若Rd≤0,即水田产流量为零,则产污量Wp=0;
若Rd>0,即水田产流,产污量按下式计算:
Wp=0.01Ca×Rp (10)
式中:Wp为单位面积稻田污染负荷,kg/(km2·d);Ca为田面水污染物浓度,mg/L;Rp为稻田净雨深,mm/d。
作为优选,所述步骤(7)中每条河道对应的面源污染量WLk为:
式中:WLk为第k条河道的面源污染负荷,kg/d;Au k为利用公式(3)计算得到的第k条河道对应的城市面积,km2;为第k条河道对应的旱地面积,km2;Ap k为第k条河道对应的稻田面积,km2。
有益效果:本发明的基于栅格化处理的平原河网区河道面源污染负荷确定方法,结合平原河网区的地貌和水系特征,充分利用GIS的栅格化离散技术和空间叠合分析功能,不依赖于地形数据,克服了现有面源污染模型无法直接应用于平原河网区的缺点,可准确计算随降雨径流汇入各条河道的面源污染量,且综合考虑了汇流距离和河道过流能力等因素的影响,易于在平原河网区的河道面源污染汇入量计算中推广应用。
附图说明
图1是我国东部某平原河网区水系图。
图2是我国东部某平原河网区河网多边形结构图。
图3是本发明的河网多边形示例图。
图4是我国东部某平原河网区栅格化示意图。
图5是本发明的图层空间运算示意图。
图6是本发明的网格权重因子计算示意图。
图7是本发明的汇流区划分示意图。
具体实施方式
下面结合附图和我国东部某平原河网区实际算例对本发明作进一步的说明。
(1)绘制河网多边形
图1为我国东部某平原河网区水系图,该地区河道纵横交错,水系构成了典型的网状结构。通过对河网结构进行拓扑分析,生成由相邻河道包围所形成的若干多边形,称为“河网多边形”,如图2所示。
河网多边形内产生面源污染会随降雨-径流汇入周边各条河道。为了计算该地区的面源污染负荷及汇流过程,首先需要判断每个河网多边形内的面源污染流入哪条河道。以其中一个河网多边形S1为例,如图3所示。多边形S1内产生的污染物有可能汇入河道L1、L2、L3和L4,因此要确定每条河道的汇流区范围。
(2)计算域栅格化及空间运算
采用正交矩形网格对计算区域进行空间离散化处理,如图4所示。理论上网格尺度越小,汇流区划分精度越高,但是网格尺度过小,计算效率会显著降低,且对划分精度的改善效果有限。经过反复测试,当网格边长小于所有河网多边形最小边长的1/10时,汇流区划分成果可以满足计算精度要求。
借助GIS软件平台的“Intersect”工具,对计算区域的土地利用图层和栅格图层进行空间叠合分析,如图5所示,生成每个网格单元不同土地利用类型的范围,而统计相应的土地利用面积为河网多边形第i个网格第j种土地利用的面积,km2。
(3)计算汇流权重因子
在确定网格内的面源污染物汇入哪条河道之前,需要计算汇流权重因子。综合考虑河道过流能力以及河网多边形的结构特征对面源污染汇流的影响,以河道断面面积与网格到周边河道最小距离之比为权重,按公式(1)计算网格汇流权重因子。
式中:为第i个网格到第k条河道上的权重因子,%;为第i个栅格到第k条河道的距离,km;Ak为第k条河道的断面面积,m2;m为包围该网格的河道个数。
图6为网格汇流权重因子计算示意图。以其中一个网格单元i为例,假设该网格到周边河道L1、L2、L3和L4的距离分别为d1、d2、d3和d4,周边河道的断面面积分别为A1、A2、A3和A4。将河道断面面积和网格距河道的距离代入式(1),计算该网格i到L1、L2、L3和L4河道的汇流权重因子
(4)确定面源污染流向
以网格到周边河道权重因子最大作为判据,确定网格内面源污染物的流向。按公式(2)计算网格到包围该网格的第k条河道的最大权重因子。
Mi=max{Pi k,k=1…m} (2)
式中:Mi为第i个网格到包围该网格的第k条河道上的最大权重因子,%。
如果第i个网格到包围该网格的第k条河道的权重因子最大,那么该网格的面源污染负荷全部汇入第k条河道。
例如,步骤(3)中求得第i个网格到周边河道的权重因子分别为假设最大,即第i个网格的面源污染全部汇入河道L4。
(5)划分面源污染汇流区
依次计算河网多边形内每个网格的面源污染流向,所有汇入同一条河道的网格即构成该河道的汇流区。如图7所示。
例如,重复按步骤(3)计算河网多边形S1内其余网格到河道L1、L2、L3和L4的汇流权重因子,按步骤(4)统计最大权重因子,判断面源污染流向。所有面源流向为同一条河道的网格形成该河道的汇流区,如图7所示。
依次对图1中的每个河网多边形重复步骤(3)~(5),得到计算区域全部河道的汇流区范围。
(6)统计土地利用类型面积
按公式(3)统计河道汇流区不同土地利用的面积。
式中:Aj为河道汇流区第j种土地利用面积,km2;n为河道汇流区的网格数量。
重复步骤(6),统计各河道汇流范围内的各种土地利用面积。
(7)计算单位面积污染负荷
计算随降雨径流流失的单位面积面源污染负荷,包括城市径流污染、旱地径流污染、稻田径流污染3种面源类型。其中城市降雨径流污染负荷采用UNPS模式计算,旱地和稻田径流污染负荷分别采用DNPS和PNPS模式计算。
a.UNPS模式
城市降雨径流产污过程可用地表污染物累积和降雨径流冲刷两个阶段进行描述。
①污染物累积模型
按公式(4)计算城市单位面积的地表污染物累积通量:
Ws=αiFiγiRcl/0.9 (4)
式中:Ws为城市单位面积污染物累积通量,kg/(km2·d);αi为城市污染物浓度参数,mg/L;γi为地面清扫频率参数;Rcl为地表污染冲刷降水量,mm/d;Fi为人口密度参数。
其中,γi=Ni/20(清扫间隔Ni<20h)
γi=1(清扫间隔Ni≥20h)
②降雨径流冲刷模型
城镇降雨径流的流失速率按公式(5)计。
Wu=Ws(1-e-kRt) (5)
式中:Wu为单位面积城市污染负荷,kg/(km2·d);k为降雨径流对地表污染物的冲刷系数,1/mm,城市地区取0.14~0.19;R为降雨强度,mm/h;t为降雨历时,h。
b.DNPS模式
利用公式(6)计算单位面积旱地污染负荷,考虑施肥量的影响。
Wf=mfη+W0 (6)
式中:η为肥料流失率,%;Wf为某一施肥水平下单位面积肥料流失量,kg/(km2·d);W0为零施肥条件下单位面积肥料年流失量,kg/(km2·d);mf为单位面积施肥量,kg/(km2·d)。
c.PNPS模式
根据稻田田面水浓度随施肥量的变化特征,从质量守恒原理出发,计算单位面积稻田径流污染负荷。
①稻田径流氮素流失模型
田面水中TN和NH3-N浓度变化过程按公式(7)和(8)计算:
式中:和为前一时刻和后一时刻的田面水深度,mm;和为前一时刻和后一时刻田面水NH3-N浓度,mg·L-1;和为前一时刻和后一时刻田面水TN浓度,mg·L-1;Ri为稻田灌溉速率,mm·d-1;Ci1和Ci2为稻田灌溉水NH3-N和TN浓度,mg·L-1;Rr,Rd,Rl分别为降水强度、实际排水速率及渗漏速率,mm·d-1;Cr1和Cr2为降水中NH3-N和TN浓度,mg·L-1;Φn为氮肥向田面水的释放通量,kg·hm-2·d-1;kv为溶液中NH3-N的挥发速率常数,d-1;kn和kdn为水土界面的硝化和反硝化速率常数,d-1。
②稻田径流磷素流失模型
田面水磷素浓度变化过程按公式(9)计算:
式中:和为前一时刻和后一时刻田面水TP的质量浓度,mg·L-1;Ri为灌溉速率,mm·d-1;Ci3为灌溉水中TP的质量浓度,mg·L-1;Cr3为降水中TP的质量浓度,mg·L-1;ka为土壤对TP的吸附速率常数,d-1;Φp为磷肥向田面水的释放通量kg·hm-2·d-1。
③单位面积稻田污染物流失量
根据(7)~(9)式计算田面水污染物浓度随时间的变化过程后,根据单位面积稻田排水量按式(10)计算随污染负荷:
若Rd≤0,即水田产流量为零,则产污量Wp=0;
若Rd>0,即水田产流,产污量按下式计算:
Wp=0.01Ca×Rp (10)
式中:Wp为单位面积稻田污染负荷,kg/(km2·d);Ca为田面水污染物浓度,mg/L;Rp为稻田净雨深,mm/d。
(8)计算河道面源污染负荷
将公式(3)以及公式(5)、(6)和(10)的计算结果带入公式(11),计算汇入每条河道的面源污染负荷。
式中:WLk为第k条河道的面源污染负荷,kg/d;Au k为利用公式(3)计算得到的第k条河道对应的城市面积,km2;为第k条河道对应的旱地面积,km2;Ap k为第k条河道对应的稻田面积,km2。
(9)根据面源污染量WLk,结合工业废水和生活污水的点源调查资料,统计每条河道各种污染物的污染负荷总量WLs,将该数据与政府公布的水污染物限制排污总量WLL进行比较:
如果WLs≤WLL,说明河道的实际污染负荷小于或等于限制排污总量,无需对该河道制定相应减排措施。
如果WLs>WLL,说明该河道的实际污染负荷已经超过限制排污总量,需要制定相应的污染治理和减排措施。对于实际污染负荷超过限制排污总量的河道,计算相应的污染物削减量。根据面源和点源污染负荷构成,明确水污染治理方向,制定分阶段的污染源削减方案。
本发明提出了一种适合平原河网区水系特征的面源污染汇流区确定方法,利用GIS空间叠合分析技术以及栅格化处理方法,克服了现有面源污染模型无法应用于平原河网区的缺点,可准确划定河道汇流范围,易于在平原河网区面源污染负荷计算中推广使用。
Claims (7)
1.一种基于栅格化处理的平原河网区河道面源污染负荷确定方法,其特征在于,包括以下步骤:
(1)绘制河网多边形结构图
根据计算区域的水系图,对河网进行概化处理,绘制由河道包围所形成的网状多边形结构图;
(2)计算域栅格化及空间运算
借助工具,对计算区域的土地利用图层和栅格图层进行空间叠合分析,生成每个网格单元不同土地利用类型的范围,统计相应土地利用的面积为河网多边形第i个网格第j种土地利用的面积,km2;
(3)计算汇流权重因子
按公式(1)进行计算:
式中:Pi k为第i个网格到第k条河道上的权重因子,%;di k为第i个网格到第k条河道的最小距离,km;Ak为第k条河道的断面面积,m2;m为包围该网格的河道数量;
(4)确定面源污染流向
统计每个网格到第k条河道的最大权重因子,按公式(2)计算:
Mi=max{Pi k,k=1…m} (2)
式中:Mi为某个河网多边形的第i个网格汇流入第k条河道的最大权重因子,%;
如果第i个网格到第k条河道的权重因子最大,那么该网格的面源污染负荷全部汇入第k条河道;
(5)面源污染汇流区划分
依次计算河网多边形内每个网格的面源污染汇入周边河道的最大权重因子,所有汇入同一条河道的网格即构成该河道的汇流区;
(6)统计汇流区土地利用类型面积
按公式(3)统计河网多边形不同土地利用类型的面积;
式中:Aj分别为河网多边形第j种土地利用类型的面积,km2;n为某一河道汇流区的栅格数量;依次对所有河网多边形重复上述计算过程,得到计算区域全部河道的汇流区范围及对应的各种土地利用类型面积;
(7)将河道汇流区的各类土地利用面积输入面源污染模型,计算与每条河道对应的面源污染量WLk;
(8)根据面源污染量WLk,结合工业废水和生活污水的点源调查资料,统计每条河道各种污染物的污染负荷总量WLs,将该数据与政府公布的水污染物限制排污总量WLL进行比较:
如果WLs≤WLL,说明河道的实际污染负荷小于或等于限制排污总量,无需对该河道制定相应减排措施;
如果WLs>WLL,说明该河道的实际污染负荷已经超过限制排污总量,需要制定相应的污染治理和减排措施。
2.根据权利要求1所述的基于栅格化处理的平原河网区河道面源污染负荷确定方法,其特征在于:所述步骤(2)中的工具为GIS软件平台的“Intersect”工具。
3.根据权利要求1所述的基于栅格化处理的平原河网区河道面源污染负荷确定方法,其特征在于:所述步骤(7)中面源污染类型包含城市径流污染、旱地径流污染、稻田径流污染。
4.根据权利要求3所述的基于栅格化处理的平原河网区河道面源污染负荷确定方法,其特征在于:城市降雨径流污染单元负荷计算包含以下两个步骤:
a、计算城市单位面积的地表污染物累积通量:
Ws=αiFiγiRcl/0.9 (4)
式中:Ws为城市单位面积污染物累积通量,kg/(km2·d);αi为城市污染物浓度参数,mg/L;γi为地面清扫频率参数;Rcl为地表污染冲刷降水量,mm/d;Fi为人口密度参数;
其中,γi=Ni/20(清扫间隔Ni<20h)
γi=1(清扫间隔Ni≥20h)
b、降雨径流冲刷模型
城镇降雨径流的流失速率按公式(5)计:
Wu=Ws(1-e-kRt) (5)
式中:Wu为单位面积城市污染负荷,kg/(km2·d);k为降雨径流对地表污染物的冲刷系数,1/mm,城市地区取0.14~0.19;R为降雨强度,mm/h;t为降雨历时,h。
5.根据权利要求4所述的基于栅格化处理的平原河网区河道面源污染负荷确定方法,其特征在于:旱地径流污染单元负荷利用式(6)计算而得
Wf=mfη+W0 (6)
式中:η为肥料流失率,%;Wf为某一施肥水平下单位面积肥料流失量,kg/(km2·d);W0为零施肥条件下单位面积肥料年流失量,kg/(km2·d);mf为单位面积施肥量,kg/(km2·d)。
6.根据权利要求5所述的基于栅格化处理的平原河网区河道面源污染负荷确定方法,其特征在于:稻田径流污染单元负荷的计算包含以下步骤,
①稻田径流氮素流失模型
田面水中TN和NH3-N浓度变化过程按公式(7)和(8)计算:
式中:和为前一时刻和后一时刻的田面水深度,mm;和为前一时刻和后一时刻田面水NH3-N浓度,mg·L-1;和为前一时刻和后一时刻田面水TN浓度,mg·L-1;Ri为稻田灌溉速率,mm·d-1;Ci1和Ci2为稻田灌溉水NH3-N和TN浓度,mg·L-1;Rr,Rd,Rl分别为降水强度、实际排水速率及渗漏速率,mm·d-1;Cr1和Cr2为降水中NH3-N和TN浓度,mg·L-1;Φn为氮肥向田面水的释放通量,kg·hm-2·d-1;kv为溶液中NH3-N的挥发速率常数,d-1;kn和kdn为水土界面的硝化和反硝化速率常数,d-1;
②稻田径流磷素流失模型
田面水磷素浓度变化过程按公式(9)计算:
式中:和为前一时刻和后一时刻田面水TP的质量浓度,mg·L-1;Ri为灌溉速率,mm·d-1;Ci3为灌溉水中TP的质量浓度,mg·L-1;Cr3为降水中TP的质量浓度,mg·L-1;ka为土壤对TP的吸附速率常数,d-1;Φp为磷肥向田面水的释放通量kg·hm-2·d-1;
③单位面积稻田污染物流失量
根据(7)~(9)式计算田面水污染物浓度随时间的变化过程后,根据单位面积稻田排水量按式(10)计算随污染负荷:
若Rd≤0,即水田产流量为零,则产污量Wp=0;
若Rd>0,即水田产流,产污量按下式计算:
Wp=0.01Ca×Rp (10)
式中:Wp为单位面积稻田污染负荷,kg/(km2·d);Ca为田面水污染物浓度,mg/L;Rp为稻田净雨深,mm/d。
7.根据权利要求6所述的基于栅格化处理的平原河网区河道面源污染负荷确定方法,其特征在于:所述步骤(7)中每条河道对应的面源污染量WLk为:
WLk=Wu×Aw k+Wf×Af k+Wp×Ap k
式中:WLk为第k条河道的面源污染负荷,kg/d;Au k为利用公式(3)计算得到的第k条河道对应的城市面积,km2;为第k条河道对应的旱地面积,km2;Ap k为第k条河道对应的稻田面积,km2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810474581.0A CN108733915B (zh) | 2018-05-17 | 2018-05-17 | 基于栅格化处理的平原河网区河道面源污染负荷确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810474581.0A CN108733915B (zh) | 2018-05-17 | 2018-05-17 | 基于栅格化处理的平原河网区河道面源污染负荷确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108733915A true CN108733915A (zh) | 2018-11-02 |
CN108733915B CN108733915B (zh) | 2019-08-20 |
Family
ID=63937548
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810474581.0A Active CN108733915B (zh) | 2018-05-17 | 2018-05-17 | 基于栅格化处理的平原河网区河道面源污染负荷确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108733915B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109598419A (zh) * | 2018-11-19 | 2019-04-09 | 北京英视睿达科技有限公司 | 一种基于水质数据确定水环境热点网格的方法及装置 |
CN110287615A (zh) * | 2019-06-28 | 2019-09-27 | 南京大学 | 一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法 |
CN111046551A (zh) * | 2019-12-10 | 2020-04-21 | 河海大学 | 一种城市群排水过程模拟方法 |
CN111310124A (zh) * | 2020-02-13 | 2020-06-19 | 北京市环境保护科学研究院 | 一种城市径流数据处理方法 |
CN111552763A (zh) * | 2020-05-13 | 2020-08-18 | 水利部交通运输部国家能源局南京水利科学研究院 | 城市面源污染负荷监测方法 |
CN111639707A (zh) * | 2020-05-29 | 2020-09-08 | 河南大学 | 一种基于gis技术的土地污染控制方法 |
CN111784052A (zh) * | 2020-07-03 | 2020-10-16 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种城市面源污染入河负荷量预测方法 |
CN112526100A (zh) * | 2021-02-09 | 2021-03-19 | 浙江省生态环境低碳发展中心 | 一种基于物联网的点面源排污权储备方法和系统 |
CN112766664A (zh) * | 2020-12-31 | 2021-05-07 | 中国科学院生态环境研究中心 | 基于gis平台的城市面源污染风险识别方法及装置 |
CN113111530A (zh) * | 2021-04-23 | 2021-07-13 | 中国水利水电科学研究院 | 基于分布式水文模型的矿山污染物扩散入河量估算方法 |
CN113361114A (zh) * | 2021-06-11 | 2021-09-07 | 中国科学院精密测量科学与技术创新研究院 | 一种基于径流路径的多尺度面源污染物入河系数测算方法 |
CN113689099A (zh) * | 2021-08-16 | 2021-11-23 | 河海大学 | 基于水量构成的环状河网水系水质考核断面污染负荷贡献率确定方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20110117412A (ko) * | 2010-04-21 | 2011-10-27 | 서울대학교산학협력단 | 반응항을 결합한 이송-분산 방정식의 수치모형을 통해 하천에 유입된 열 오염물의 거동을 해석하는 방법 |
CN102663267A (zh) * | 2012-05-15 | 2012-09-12 | 南京大学 | 一种半湿润区流域面源污染负荷的确定方法 |
CN103020424A (zh) * | 2012-11-22 | 2013-04-03 | 北京师范大学 | 一种基于人工降雨实验的北方平原农田区非点源污染负荷估算方法 |
CN104091052A (zh) * | 2014-06-25 | 2014-10-08 | 安徽省环境科学研究院 | 一种非点源污染空间分布式模型 |
CN104217064A (zh) * | 2014-07-23 | 2014-12-17 | 中国水利水电科学研究院 | 一种农村生活类面源负荷的空间展布方法及装置 |
CN104361523A (zh) * | 2014-11-06 | 2015-02-18 | 浙江大学 | 一种基于gis的分布式稻田氮素径流流失负荷估算方法 |
-
2018
- 2018-05-17 CN CN201810474581.0A patent/CN108733915B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20110117412A (ko) * | 2010-04-21 | 2011-10-27 | 서울대학교산학협력단 | 반응항을 결합한 이송-분산 방정식의 수치모형을 통해 하천에 유입된 열 오염물의 거동을 해석하는 방법 |
CN102663267A (zh) * | 2012-05-15 | 2012-09-12 | 南京大学 | 一种半湿润区流域面源污染负荷的确定方法 |
CN103020424A (zh) * | 2012-11-22 | 2013-04-03 | 北京师范大学 | 一种基于人工降雨实验的北方平原农田区非点源污染负荷估算方法 |
CN104091052A (zh) * | 2014-06-25 | 2014-10-08 | 安徽省环境科学研究院 | 一种非点源污染空间分布式模型 |
CN104217064A (zh) * | 2014-07-23 | 2014-12-17 | 中国水利水电科学研究院 | 一种农村生活类面源负荷的空间展布方法及装置 |
CN104361523A (zh) * | 2014-11-06 | 2015-02-18 | 浙江大学 | 一种基于gis的分布式稻田氮素径流流失负荷估算方法 |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109598419A (zh) * | 2018-11-19 | 2019-04-09 | 北京英视睿达科技有限公司 | 一种基于水质数据确定水环境热点网格的方法及装置 |
CN110287615A (zh) * | 2019-06-28 | 2019-09-27 | 南京大学 | 一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法 |
CN111046551A (zh) * | 2019-12-10 | 2020-04-21 | 河海大学 | 一种城市群排水过程模拟方法 |
CN111046551B (zh) * | 2019-12-10 | 2020-11-20 | 河海大学 | 一种城市群排水过程模拟方法 |
CN111310124A (zh) * | 2020-02-13 | 2020-06-19 | 北京市环境保护科学研究院 | 一种城市径流数据处理方法 |
CN111552763A (zh) * | 2020-05-13 | 2020-08-18 | 水利部交通运输部国家能源局南京水利科学研究院 | 城市面源污染负荷监测方法 |
CN111639707A (zh) * | 2020-05-29 | 2020-09-08 | 河南大学 | 一种基于gis技术的土地污染控制方法 |
CN111784052A (zh) * | 2020-07-03 | 2020-10-16 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种城市面源污染入河负荷量预测方法 |
CN111784052B (zh) * | 2020-07-03 | 2023-10-17 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种城市面源污染入河负荷量预测方法 |
CN112766664B (zh) * | 2020-12-31 | 2023-05-23 | 中国科学院生态环境研究中心 | 基于gis平台的城市面源污染风险识别方法及装置 |
CN112766664A (zh) * | 2020-12-31 | 2021-05-07 | 中国科学院生态环境研究中心 | 基于gis平台的城市面源污染风险识别方法及装置 |
CN112526100A (zh) * | 2021-02-09 | 2021-03-19 | 浙江省生态环境低碳发展中心 | 一种基于物联网的点面源排污权储备方法和系统 |
US11382255B1 (en) | 2021-02-09 | 2022-07-12 | Zhejiang Province Eco-Environmental Low-Carbon Development Center | Method and system for storing emission rights for point and nonpoint source pollution based on internet of things |
CN113111530B (zh) * | 2021-04-23 | 2021-11-02 | 中国水利水电科学研究院 | 基于分布式水文模型的矿山污染物扩散入河量估算方法 |
CN113111530A (zh) * | 2021-04-23 | 2021-07-13 | 中国水利水电科学研究院 | 基于分布式水文模型的矿山污染物扩散入河量估算方法 |
CN113361114B (zh) * | 2021-06-11 | 2022-05-17 | 中国科学院精密测量科学与技术创新研究院 | 一种基于径流路径的多尺度面源污染物入河系数测算方法 |
CN113361114A (zh) * | 2021-06-11 | 2021-09-07 | 中国科学院精密测量科学与技术创新研究院 | 一种基于径流路径的多尺度面源污染物入河系数测算方法 |
CN113689099A (zh) * | 2021-08-16 | 2021-11-23 | 河海大学 | 基于水量构成的环状河网水系水质考核断面污染负荷贡献率确定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108733915B (zh) | 2019-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108733915B (zh) | 基于栅格化处理的平原河网区河道面源污染负荷确定方法 | |
CN108717453B (zh) | 一种基于gis平台的平原河网区污染负荷计算方法 | |
CN108664647B (zh) | 一种集成水环境模型的流域精细化管理系统 | |
Bui et al. | Integration of SWAT and QUAL2K for water quality modeling in a data scarce basin of Cau River basin in Vietnam | |
CN112526100B (zh) | 一种基于物联网的点面源排污权储备方法和系统 | |
CN107563139A (zh) | 一种点源突发性水污染源事故溯源贡献程度计算方法 | |
CN109598428B (zh) | 一种基于行政单元和水系的污染物削减分配方法 | |
CN108984823A (zh) | 一种合流制溢流调蓄池规模的确定方法 | |
Hinsby et al. | Threshold values and management options for nutrients in a catchment of a temperate estuary with poor ecological status | |
Chen et al. | Artificial neural network modelling of concentrations of nitrogen, phosphorus and dissolved oxygen in a non‐point source polluted river in Zhejiang Province, southeast China | |
Kalin et al. | Evaluation of sediment transport models and comparative application of two watershed models | |
CN114169266B (zh) | 一种河流污染源解析和污染物运移计算方法及系统 | |
Wang et al. | Accurately early warning to water quality pollutant risk by mobile model system with optimization technology | |
CN110765213A (zh) | 一种地表水流域污染源排放清单(动态清单)的编制方法 | |
CN114462698B (zh) | 一种面向流域汇水区域的磷排放污染负荷预测方法 | |
CN111291305A (zh) | 基于一维水质模型计算流域污染源入河系数的方法 | |
CN112464434B (zh) | 一种灌溉农用地面源污染排放量预测方法及系统 | |
Venohr et al. | Nitrogen retention in a river system and the effects of river morphology and lakes | |
Yang et al. | Pollution load estimation and control countermeasures of Zhangze reservoir | |
Liu et al. | Research on nonpoint source pollution assessment method in data sparse regions: a case study of Xichong River Basin, China | |
CN115935615A (zh) | 场次降雨径流冲刷面源污染负荷估算方法及系统 | |
CN116384268A (zh) | 基于水污染的污染源溯源方法、装置以及设备 | |
Farzin et al. | Developing the performance of modern methods using multi-objective optimization in urban runoff control | |
Tuppad et al. | SWAT discharge response to spatial rainfall variability in a Kansas watershed | |
Tapia et al. | Economic valuation of the natural service of nitrate regulation provided by rivers including dilution effects: Application to a semiarid region, the Ebro basin (Spain) |
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 | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20181102 Assignee: Suzhou Sucheng Environmental Technology Co.,Ltd. Assignor: HOHAI University Contract record no.: X2021320000077 Denomination of invention: Determination method of river non-point source pollution load in plain river network area based on grid processing Granted publication date: 20190820 License type: Common License Record date: 20210831 |
|
EE01 | Entry into force of recordation of patent licensing contract |