CN108717453A - 一种基于gis平台的平原河网区污染负荷计算方法 - Google Patents
一种基于gis平台的平原河网区污染负荷计算方法 Download PDFInfo
- Publication number
- CN108717453A CN108717453A CN201810486632.1A CN201810486632A CN108717453A CN 108717453 A CN108717453 A CN 108717453A CN 201810486632 A CN201810486632 A CN 201810486632A CN 108717453 A CN108717453 A CN 108717453A
- Authority
- CN
- China
- Prior art keywords
- pollutant
- runoff
- pollution
- formula
- rainfall
- 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Processing Of Solid Wastes (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明开发了一种基于GIS平台的平原河网区污染负荷计算方法,该方法针对平原河网区污染源产生和迁移的特点,在GIS平台支撑下,计算各类污染源和各种污染物的产生量及入河量,采用4种计算模式分别计算包括城镇生活污染、农村生活污染、城镇降雨径流污染、旱地降雨径流污染、稻田降雨径流污染、畜禽养殖和渔业养殖7种污染源在内的点源和面源污染负荷,最终确定整个平原河网区的污染负荷。该方法既能反映平原河网区污染物迁移特征,又能在较大的时空尺度上加以推广应用,提高了污染物从源头发生向河网迁移过程的模拟精度。
Description
技术领域
本发明属于环境领域,尤其涉及一种基于GIS平台的平原河网区污染负荷计算方法。
背景技术
污染负荷是指通过各种途径进入受纳水体的污染物数量,包括点源和面源污染物,即地表水体接纳的污染物量,对污染负荷的定量化研究是流域污染治理的重要基础性工作。对于污染负荷中的点源污染,可以通过调查和监测等手段获得废水和污染物排放量;但是对于面源污染,由于其发生的随机性、机理过程的复杂性、排放途径及排放污染物的不确定性以及时空分布的差异性,难以在较大的时空尺度上通过调查和监测获得污染负荷量。尤其是对于地势相对平坦的平原河网地区,由于其汇流区界限很难确定,水文和水动力特征十分复杂,使得试图通过野外监测获得面源污染负荷的方法存在较大难度。
国外早在20世纪60年代就开展了面源污染模型的研究,提出了大量面源污染模型,但这些模型大多是针对机械化大农场开发的,在中国以人工耕作为主的耕种方式下,一些参数及模块需要调整以适应我国的实际情况。此外,国外有关面源污染的定量化研究大多从污染物的迁移和转化出发,建立面向过程的机理模型,模型参数众多,操作十分复杂,这对于污染物的迁移特征比较特殊的平原河网区,在缺乏充足的基础资料和试验数据情况下,难以直接将国外的面源模型在较大的空间尺度上加以运用。因此。如何建立既能反映我国平原河网区污染物迁移特点,又能在较大的空间尺度上加以推广的污染负荷计算方法成为关键。
近年来,GIS技术在污染负荷定量化研究中的应用成为一种发展方向。GIS技术的发展使得污染负荷研究所需空间信息数据的数量和质量都大大提高,极大地推进了污染负荷计算方法。因此建立平原河网区的污染负荷模型需要在吸收国外污染源定量化研究成果的基础上,结合平原河网区污染物产生和迁移的特点,充分利用GIS技术手段,建立具有一定理论基础,形式更为简单,应用更加方便的污染负荷计算方法。
发明内容
发明目的:为了克服现有污染负荷估算方法存在的不足,本发明提供一种基于GIS平台的平原河网区污染负荷计算方法。提出一种既能反映平原河网区污染物迁移特征,又能在较大的时空尺度上加以推广应用的污染负荷计算方法,提高了污染物从源头发生向河网迁移过程的模拟精度。
技术方案:为实现本发明的目的,本发明所采用的技术方案是:一种基于GIS平台的平原河网区污染负荷计算方法,包括以下步骤:
(1)将污染负荷计算分为产生量计算和入河量计算两部分;
(2)对不同的污染源使用不同的模式计算污染负荷产生量;
(3)基于GIS统计各种土地利用类型面积;
(4)计算平原河网区污染负荷入河量。
其中,所述污染源分为城镇生活污染、农村生活污染、城镇降雨径流污染、旱地降雨径流污染、稻田降雨径流污染、畜禽养殖和渔业养殖。
其中,在步骤(2)中,城镇生活、农村生活、畜禽养殖和渔业养殖与降雨-径流无关的污染负荷采用PROD模式计算;城镇降雨径流迁移的污染负荷采用UNPS模式计算;随旱地和稻田降雨径流迁移的污染负荷分别采用DNPS和PNPS模式计算。
其中,在步骤(3)中,统计各种土地利用类型面积方法如下:将平原河网区的土地利用类型分为城镇,包括工业区、商业区和生活区、旱地、稻田和水面,以某一分区作为污染负荷的计算单元,采用GIS平台统计各个计算单元各种土地利用类型的面积。
其中,PROD模式计算过程:
PROD模式也可称为排污系数法,用于计算与降雨-径流无关污染源的污染负荷产生量,按公式(1)计算:
式中:为第i种污染源第j种污染物的污染物产生量;Ni为第i种污染源的数量;为第i种污染源第j种污染物的污染负荷当量;计算城镇和农村居民的污染物产生量时,Ni为城镇和农村居民数量,为城镇和农村居民的排污系数;计算畜禽养殖产污量时,Ni为畜禽数量,为畜禽的排污系数;计算水产养殖产污量时,Ni为水产养殖产量,为水产养殖品种的排污系数。
其中,所述UNPS模式计算过程如下:
①污染物累积模型
将城镇下垫面分为工业区、商业区和生活区3种类型,按公式(2)计算3种城镇土地类型单位面积的地表污染物累积通量:
Xi=αiFiγiRcl/0.9 (2)
式中:Xi为第i种土地类型单位面积的污染物累积通量,kg/(km2·d);αi为城市污染物浓度参数,mg/L;γi为地面清扫频率参数;Rcl为地表污染冲刷降水量,mm/d;Fi为人口密度参数;
其中,γi=Ni/20,清扫间隔Ni<20h
γi=1,清扫间隔Ni≥20h,h单位为小时。
按公式(2)分别计算城镇各种土地利用类型的污染物累积通量后,再按公式(3)计算城镇地表污染物的总累积量:
式中:P为城镇地表污染物的累积速率,kg/d;Pi为第i种土地类型的污染物累积速率,kg/d;Xi为第i种土地类型单位面积的污染物累积速率,kg/(km2·d1);Ai为第i种土地类型的面积,km2;n为土地利用类型个数,包括工业区、商业区和生活区;
若某日的降雨量小于污染物降雨阈值,则地表污染物的累积量按公式(2)和(3)计算;若某日降雨量大于该阈值,则假设该日地表污染物的累积量为0;
②降雨径流冲刷模型
城镇降雨径流的冲刷速率按公式(4)计:
Pt=P(1-e-kRt)(4)
式中:Pt为降雨历时t的地表污染物冲刷速率,kg/d;P为城镇地表污染物的累积速率,kg/d;k为降雨径流对地表污染物的冲刷系数,1/mm;R为城镇的降雨强度,mm/h;经过降雨径流冲刷后的地表污染物剩余量作为后续地表污染物的累积量计算。
其中,所述DNPS模式计算过程:
①建立单位面积农田肥料年流失量与年流失率和施肥量的经验关系,计算得到年流失量:
Wf=mfη+W0 (5)
式中:η为肥料年流失率,%;Wf为预设施肥水平下单位面积肥料年流失量,g/hm2;W0为零施肥条件下单位面积肥料年流失量,g/hm2;mf为单位面积年施肥量,kg/hm2;
②根据农田单位面积年径流量,使用净雨深代表单位面积的年径流量,计算出径流中各种污染物的年平均浓度;
③根据农田逐日净雨深,计算旱地污染物随降雨径流的流失过程;
若Rd=0,即旱地产流量为零,则污染物流失量Wd=0;
若Rd>0,即旱地产流量不为零,相应污染物日流失量按下式计算:
式中:Wd为旱地污染物日流失量,kg;Hs为旱地标准年净雨深,mm;Rd为旱地日净雨深,mm;Ad为计算单元内的旱地面积,hm2。
其中,所述PNPS模式计算过程:
①稻田径流氮素流失模型
田面水中TN和NH3-N浓度变化过程计算如下:
式中:和为前一时刻和后一时刻的田面水深度,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。
②稻田径流磷素流失模型
田面水TP浓度变化过程计算如下:
式中:和为前一时刻和后一时刻田面水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;
③稻田径流耗氧有机物流失模型
式中:和为前一时刻和后一时刻田面水有机物的质量浓度,mg/L;Ri为灌溉速率,mm/d;Ci4为灌溉水中有机物的质量浓度,mg/L;Cr4为降水中有机物的质量浓度,mg/L;Cmax为田面水有机物浓度上限,mg/L;T为田面水有机物释放周期,d。
根据(8)~(10)式计算田面水污染物浓度随时间的变化过程后,根据稻田的排水量按式(11)计算随径流流失的污染物负荷:
若Rd≤0,即水田产流量为零,则产污量Wp=0;
若Rd>0,即水田产流,产污量按下式计算:
Wp=0.01Ca×Rp×Ap(11)
式中:Wp为稻田日产污量,kg;Ca为田面水污染物浓度,mg/L;Rp为稻田日净雨深,mm;Ap为计算单元内的稻田面积,hm2。
其中,在步骤(4)中,平原河网区污染负荷入河量计算方法如下:
式中:Wei为第i种污染源的污染物入河量,kg/d;Wpi为第i种污染源的污染物产生量,kg/d,包括PROD模式的UNPS模式的Pt、DNPS模式的Wd和PNPS模式的Wp;pij为第i种污染源第j条入河路径的比例系数;m为第i种污染源入河路径的数量;fk为第k种处理单元的处理效率,处理单元包括对化粪池、雨污水管网、农村生活污水处理、畜禽养殖污染物处理、湖荡支浜和土壤6种;n为第i种污染源第j条入河路径对应的处理单元数量。
有益效果:与现有技术相比,本发明具有以下优点:
(1)该方法计算的污染源包括城镇生活污染、农村生活污染、城镇降雨径流污染、旱地降雨径流污染、稻田降雨径流污染、畜禽养殖和渔业养殖等7种,覆盖了平原河网区绝大部分污染源类型。
(2)利用GIS平台的空间运算功能,快速统计各个计算单元各种土地利用类型的面积,为污染负荷计算提供数据支撑。
(3)根据不同污染源的产污特征,采用PROD、UNPS、DNPS和PNPS等4种模式分别计算相应污染源的污染物产生量。
(4)根据平原河网区各类污染源入河过程及排放特征,采用入河路径比例系数以及处理单元的处理效率计算污染物入河量。
附图说明
图1是本发明分布式污染负荷模型路径框图。
具体实施方式
下面结合附图对本发明作更进一步的说明。如图1所示,本发明提出了一种基于GIS平台的平原河网区污染负荷计算方法如下:
(1)计算方法基本结构
该计算方法可以分为污染负荷产生量计算和处理量计算两部分。前者用于计算各种污染源各类污染物的产生量,包括4种计算模式。后者计算各种污染源经过各条入河途径和各个处理单元处理后的污染物入河量,包括6种处理单元。该计算方法的结构框图如图1所示。
(2)污染负荷产生量计算
该方法计算的污染源包括城镇生活污染、农村生活污染、城镇降雨径流污染、旱地降雨径流污染、稻田降雨径流污染、畜禽养殖和渔业养殖等7种。
其中城镇生活、农村生活、畜禽养殖和渔业养殖等与降雨-径流无关的污染负荷采用PROD模式计算;随城镇降雨径流迁移的污染负荷采用UNPS模式计算;随旱地和稻田降雨径流迁移的污染负荷分别采用DNPS和PNPS模式计算。根据污染源类型选用相应的计算模式。
(3)统计各种土地利用类型面积
将平原河网区的土地利用类型分为城镇,城镇包括工业区、商业区和生活区、旱地、稻田和水面。以某一分区作为污染负荷的计算单元,采用GIS平台统计各个计算单元各种土地利用类型的面积。其中,分区可以是行政分区、也可以是水利分区和水资源分区,这个根据需要选择。
例如,以县级行政区作为计算单元,采用GIS平台对土地利用类型图层与县级行政区图层进行空间交集运算,得到每个县级行政区城镇,城镇包括工业区、商业区和生活区、旱地、稻田和水面的面积。
(4)PROD模式计算过程
PROD模式也可称为排污系数法,用于计算与降雨-径流无关污染源的污染负荷产生量,按公式(1)计算:
其中,为第i种污染源第j种污染物的污染物产生量;Ni为第i种污染源的数量;为第i种污染源第j种污染物的污染负荷当量;计算城镇和农村居民的污染物产生量时,Ni为城镇和农村居民数量,为城镇和农村居民的排污系数;计算畜禽养殖产污量时,Ni为畜禽数量,为畜禽的排污系数;计算水产养殖产污量时,Ni为水产养殖产量,为水产养殖品种的排污系数。
对于不同的污染源,公式中变量的具体含义有所差别。例如,计算城镇居民的污染物产生量时,Ni为城镇居民数量,可通过查找统计年鉴获取;为城镇居民的排污系数,可通过全国污染源普查城镇生活源产排污系数手册获取。
(5)UNPS模式计算过程
城市降雨径流的产污过程可用地表污染物累积和降雨径流冲刷两个阶段加以描述,地表污染物累积量和冲刷量分别采用污染物累积模型和降雨径流冲刷模型进行计算。
①污染物累积模型
将城镇下垫面分为工业区、商业区和生活区等3种类型。按公式(2)计算3种城镇土地类型单位面积的地表污染物累积通量:
Xi=αiFiγiRcl/0.9 (2)
式中:Xi为第i种土地类型单位面积的污染物累积通量,kg/(km2·d);αi为城市污染物浓度参数,mg/L;γi为地面清扫频率参数;Rcl为地表污染冲刷降水量,mm/d;Fi为人口密度参数。
其中,γi=Ni/20,清扫间隔Ni<20h
γi=1,清扫间隔Ni≥20h,h单位为小时。
按公式(2)分别计算城镇各种土地利用类型的污染物累积通量后,再按公式(3)计算城镇地表污染物的总累积量:
式中:P为城镇地表污染物的累积速率,kg/d;Pi为第i种土地类型的污染物累积速率,kg/d;Xi为第i种土地类型单位面积的污染物累积速率,kg/(km2·d1);Ai为第i种土地类型的面积,km2;n为土地类型个数。
若某日的降雨量小于污染物降雨阈值,则地表污染物的累积量按公式(2)和(3)计算;若某日降雨量大于该阈值,则假设该日地表污染物的累积量为0。
②降雨径流冲刷模型
城镇降雨径流的冲刷速率按公式(4)计。
Pt=P(1-e-kRt)(4)
式中:Pt为降雨历时t的地表污染物冲刷速率,kg/d;P为城镇地表污染物的累积速率,kg/d;k为降雨径流对地表污染物的冲刷系数,1/mm,城市地区取0.14~0.19;R为城镇的降雨强度,mm/h。
经过降雨径流冲刷后的地表污染物剩余量作为后续地表污染物的累积量计算。
(6)DNPS模式计算过程
DNPS用于计算随旱地降雨径流迁移的污染负荷产生量。考虑不同计算单元施肥量的差异对随降雨径流流失的旱地污染负荷的影响,具体步骤如下:
①建立单位面积农田肥料年流失量与年流失率和施肥量的经验关系,计算得到年流失量。
Wf=mfη+W0(5)
式中:η为肥料年流失率,%;Wf为某一施肥水平下单位面积肥料年流失量,g/hm2;W0为零施肥条件下单位面积肥料年流失量,g/hm2;mf为单位面积年施肥量,kg/hm2。
②根据农田单位面积年径流量(净雨深),计算出径流中各种污染物的年平均浓度。
③根据农田逐日净雨深,计算旱地污染物随降雨径流的流失过程。
若Rd=0,即旱地产流量为零,则污染物流失量Wd=0;
若Rd>0,即旱地产流,相应污染物日流失量按下式计算:
式中:Wd为旱地污染物日流失量,kg;Hs为旱地标准年净雨深,mm;Rd为旱地日净雨深,mm;Ad为计算单元内的旱地面积,hm2。
例如,根据研究区所处地域、作物种植类型和模式、施肥量等,查阅相关文献,分析田间试验成果,总结肥料年流失率和零施肥条件下单位面积肥料年流失量。单位面积年施肥量可以通过查阅当地统计年鉴获取。将相关数据带入公式(5),根据水文模型计算旱地日净雨深,最后带入公式(6)即可求出旱地污染负荷产生量。
(7)PNPS模式计算过程
PNPS模式用于计算随稻田降雨径流流失的污染负荷。根据稻田田面水浓度随施肥量的变化特征,从质量守恒原理出发,考虑影响田面水浓度变化的各种因素,尤其是稻季不同阶段施肥量对田面水浓度的影响,建立稻田营养盐运移转化模型,预测稻田营养盐的径流损失量。
①稻田径流氮素流失模型
田面水中TN和NH3-N浓度变化过程计算如下:
式中:和为前一时刻和后一时刻的田面水深度,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。
②稻田径流磷素流失模型
田面水TP浓度变化过程计算如下:
式中:和为前一时刻和后一时刻田面水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。
③稻田径流耗氧有机物(COD、BOD)流失模型
式中:和为前一时刻和后一时刻田面水有机物的质量浓度,mg/L;Ri为灌溉速率,mm/d;Ci4为灌溉水中有机物的质量浓度,mg/L;Cr4为降水中有机物的质量浓度,mg/L;Cmax为田面水有机物浓度上限,mg/L;T为田面水有机物释放周期,d。
④稻田径流污染物流失量
根据(8)~(10)式计算田面水污染物浓度随时间的变化过程后,根据稻田的排水量按式(11)计算随径流流失的污染物负荷:
若Rd≤0,即水田产流量为零,则产污量Wp=0;
若Rd>0,即水田产流量不为零,则产污量按下式计算:
Wp=0.01Ca×Rp×Ap (11)
式中:Wp为稻田日产污量,kg;Ca为田面水污染物浓度,mg/L;Rp为稻田日净雨深,mm;Ap为计算单元内的稻田面积,hm2。
(8)污染负荷入河量计算
该部分用于计算各类污染源的污染物入河量。入河量根据污染负荷产生量、各条污染路径的比例系数以及各种处理单元的处理效率计算得到,公式如下:
式中:Wei为第i种污染源的污染物入河量,kg/d;Wpi为第i种污染源的污染物产生量,kg/d;pij为第i种污染源第j条入河路径的比例系数;m为第i种污染源入河路径的数量;fk为第k种处理单元的处理效率,处理单元包括化粪池、雨污水管网、农村生活污水处理、畜禽养殖污染物处理、湖荡支浜和土壤等6种,处理单元代表各种水污染处理设施;n为第i种污染源第j条入河路径对应的处理单元数量。
以城镇生活污染为例,如图1所示,城镇生活总共有如下6种入河途径:
①城镇生活污染——化粪池——污水管网——污水处理厂——湖荡
②城镇生活污染——化粪池——污水管网——污水处理厂
③城镇生活污染——化粪池——污水管网——湖荡
④城镇生活污染——化粪池——污水管网
⑤城镇生活污染——湖荡
⑥城镇生活污染——直排水体
因此,m=6,以其中第1条路径为例,该路径包括化粪池、污水管网、污水处理厂、湖荡等4种处理单元,则n=4。将城镇生活污染产生量Wpi、每条路径的污水量占比pij和每个处理单元对污染物的处理率fk带入公式(12),可以计算出城镇生活污染的入河量Wei。
Claims (10)
1.一种基于GIS平台的平原河网区污染负荷计算方法,其特征在于,包括以下步骤:
(1)将污染负荷计算分为污染产生量计算和入河量计算两部分;
(2)对不同的污染源使用不同的模式计算污染产生量;
(3)基于GIS统计各种土地利用类型面积;
(4)计算平原河网区污染入河量。
2.根据权利要求2所述的一种基于GIS平台的平原河网区污染负荷计算方法,其特征在于,所述污染源分为城镇生活污染、农村生活污染、城镇降雨径流污染、旱地降雨径流污染、稻田降雨径流污染、畜禽养殖和渔业养殖。
3.根据权利要求1所述的一种基于GIS平台的平原河网区污染负荷计算方法,其特征在于,在步骤(2)中,城镇生活、农村生活、畜禽养殖和渔业养殖与降雨-径流无关的污染产生量采用PROD模式计算;城镇降雨径流迁移的污染产生量采用UNPS模式计算;随旱地和稻田降雨径流迁移的污染产生量分别采用DNPS和PNPS模式计算。
4.根据权利要求1所述的一种基于GIS平台的平原河网区污染负荷计算方法,其特征在于,在步骤(3)中,统计各种土地利用类型面积方法如下:将平原河网区的土地利用类型分为城镇,包括工业区、商业区和生活区、旱地、稻田和水面,以某一分区作为污染负荷的计算单元,采用GIS平台统计各个计算单元各种土地利用类型的面积。
5.根据权利要求3所述的一种基于GIS平台的平原河网区污染负荷计算方法,其特征在于,PROD模式计算过程:
PROD模式也可称为排污系数法,用于计算与降雨-径流无关污染源的污染负荷产生量,按公式(1)计算:
式中:为第i种污染源第j种污染物的污染物产生量;Ni为第i种污染源的数量;为第i种污染源产污系数;计算城镇和农村居民的污染物产生量时,Ni为城镇和农村居民数量,为城镇和农村居民的产污系数;计算畜禽养殖产污量时,Ni为畜禽数量,为畜禽的产污系数;计算水产养殖产污量时,Ni为水产养殖产量,为水产养殖品种的产污系数。
6.根据权利要求3所述的一种基于GIS平台的平原河网区污染负荷计算方法,其特征在于,所述UNPS模式计算过程如下:
①污染物累积模型
将城镇下垫面分为工业区、商业区和生活区3种类型,按公式(2)计算3种城镇土地类型单位面积的地表污染物累积通量:
Xi=αiFiγiRcl/0.9 (2)
式中:Xi为第i种土地类型单位面积的污染物累积通量,kg/(km2·d);αi为城市污染物浓度参数,mg/L;γi为地面清扫频率参数;Rcl为地表污染冲刷降水量,mm/d;Fi为人口密度参数;
其中,当清扫间隔Ni<20h,γi=Ni/20;当清扫间隔Ni≥20h,γi=1,h单位为小时;
按公式(2)分别计算城镇各种土地利用类型的污染物累积通量后,再按公式(3)计算城镇地表污染物的总累积量:
式中:P为城镇地表污染物的累积速率,kg/d;Pi为第i种土地类型的污染物累积速率,kg/d;Xi为第i种土地类型单位面积的污染物累积速率,kg/(km2·d1);Ai为第i种土地类型的面积,km2;n为土地利用类型个数,包括工业区、商业区和生活区;
若某日的降雨量小于污染物降雨阈值,则地表污染物的累积量按公式(2)和(3)计算;若某日降雨量大于该阈值,则假设该日地表污染物的累积量为0;
②降雨径流冲刷模型
城镇降雨径流的冲刷速率按公式(4)计:
Pt=P(1-e-kRt)(4)
式中:Pt为降雨历时t的地表污染物冲刷速率,kg/d,即降雨污染物产生量;P为城镇地表污染物的累积速率,kg/d;k为降雨径流对地表污染物的冲刷系数,1/mm;R为城镇的降雨强度,mm/h;经过降雨径流冲刷后的地表污染物剩余量作为后续地表污染物的累积量计算。
7.根据权利要求3所述的一种基于GIS平台的平原河网区污染负荷计算方法,其特征在于,所述DNPS模式计算过程:
①建立单位面积农田肥料年流失量与年流失率和施肥量的经验关系,计算得到年流失量:
Wf=mfη+W0(5)
式中:η为肥料年流失率,%;Wf为预设施肥水平下单位面积肥料年流失量,g/hm2;W0为零施肥条件下单位面积肥料年流失量,g/hm2;mf为单位面积年施肥量,kg/hm2;
②根据农田单位面积年径流量,使用净雨深代表单位面积的年径流量,计算出径流中各种污染物的年平均浓度;
③根据农田逐日净雨深,计算旱地污染物随降雨径流的流失过程;
若Rd=0,即旱地产流量为零,则污染物流失量Wd=0;
若Rd>0,即旱地产流量不为零,相应污染物日流失量按下式计算:
式中:Wd为旱地污染物日流失量,kg,即旱地污染物产生量;Hs为旱地标准年净雨深,mm;Rd为旱地日净雨深,mm;Ad为计算单元内的旱地面积,hm2。
8.根据权利要求3所述的一种基于GIS平台的平原河网区污染负荷计算方法,其特征在于,所述PNPS模式计算过程:
①稻田径流氮素流失模型
田面水中TN和NH3-N浓度变化过程计算如下:
式中:和为前一时刻和后一时刻的田面水深度,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,Δt代表计算的时间步长。
②稻田径流磷素流失模型
田面水TP浓度变化过程计算如下:
式中:和为前一时刻和后一时刻田面水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;
③稻田径流耗氧有机物流失模型
式中:和为前一时刻和后一时刻田面水有机物的质量浓度,mg/L;Ri为灌溉速率,mm/d;Ci4为灌溉水中有机物的质量浓度,mg/L;Cr4为降水中有机物的质量浓度,mg/L;Cmax为田面水有机物浓度上限,mg/L;T为田面水有机物释放周期,d。
9.根据权利要求8所述的一种基于GIS平台的平原河网区污染负荷计算方法,其特征在于,根据(8)~(10)式计算田面水污染物浓度随时间的变化过程后,根据稻田的排水量按式(11)计算随径流流失的污染物负荷:
若Rd≤0,即水田产流量为零,则产污量Wp=0;
若Rd>0,即水田产流不为零,产污量按下式计算:
Wp=0.01Ca×Rp×Ap(11)
式中:Wp为稻田日产污量,kg,即稻田污染产生量;Ca为田面水污染物浓度,mg/L;Rp为稻田日净雨深,mm;Ap为计算单元内的稻田面积,hm2。
10.根据权利要求3所述的一种基于GIS平台的平原河网区污染负荷计算方法,其特征在于,在步骤(4)中,平原河网区污染入河量计算方法如下:
式中:Wei为第i种污染源的污染物入河量,kg/d;Wpi为第i种污染源的污染物产生量,kg/d,包括PROD模式的UNPS模式的Pt、DNPS模式的Wd和PNPS模式的Wp;pij为第i种污染源第j条入河路径的比例系数;m为第i种污染源入河路径的数量;fk为第k种处理单元的处理效率,处理单元包括对化粪池、雨污水管网、农村生活污水处理、畜禽养殖污染物处理、湖荡支浜和土壤6种;n为第i种污染源第j条入河路径对应的处理单元数量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810486632.1A CN108717453B (zh) | 2018-05-21 | 2018-05-21 | 一种基于gis平台的平原河网区污染负荷计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810486632.1A CN108717453B (zh) | 2018-05-21 | 2018-05-21 | 一种基于gis平台的平原河网区污染负荷计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108717453A true CN108717453A (zh) | 2018-10-30 |
CN108717453B CN108717453B (zh) | 2020-05-05 |
Family
ID=63900097
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810486632.1A Active CN108717453B (zh) | 2018-05-21 | 2018-05-21 | 一种基于gis平台的平原河网区污染负荷计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108717453B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110287615A (zh) * | 2019-06-28 | 2019-09-27 | 南京大学 | 一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法 |
CN110826920A (zh) * | 2019-11-08 | 2020-02-21 | 河北省南运河河务管理处 | 海岸带污染物负荷估算方法 |
CN111077280A (zh) * | 2020-01-14 | 2020-04-28 | 浙江清华长三角研究院 | 一种基于河网的农村污水处理设施与水质监测站点间的溯源分析方法 |
CN111310124A (zh) * | 2020-02-13 | 2020-06-19 | 北京市环境保护科学研究院 | 一种城市径流数据处理方法 |
CN111523088A (zh) * | 2020-04-13 | 2020-08-11 | 杭州领见数据科技有限公司 | 一种基于dpsir模型的生态环境评价方法 |
CN111639707A (zh) * | 2020-05-29 | 2020-09-08 | 河南大学 | 一种基于gis技术的土地污染控制方法 |
CN113111530A (zh) * | 2021-04-23 | 2021-07-13 | 中国水利水电科学研究院 | 基于分布式水文模型的矿山污染物扩散入河量估算方法 |
CN113256469A (zh) * | 2021-05-14 | 2021-08-13 | 郑州大学环境技术咨询工程有限公司 | 一种城镇地表径流面源污染物入河量核算方法 |
CN113361114A (zh) * | 2021-06-11 | 2021-09-07 | 中国科学院精密测量科学与技术创新研究院 | 一种基于径流路径的多尺度面源污染物入河系数测算方法 |
CN113552305A (zh) * | 2021-07-21 | 2021-10-26 | 中国水利水电科学研究院 | 一种多孔纤维材料埋设对产流水质影响的机理识别方法 |
CN115629190A (zh) * | 2022-12-08 | 2023-01-20 | 中农创达(北京)环保科技有限公司 | 农业面源污染数据管理方法、装置、设备及介质 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008210362A (ja) * | 2007-02-27 | 2008-09-11 | Inha-Industry Partnership Inst | Gisを利用した水質汚染負荷算定用の意思決定支援システム及びその運営方法 |
KR20120023259A (ko) * | 2010-09-01 | 2012-03-13 | 한국건설기술연구원 | 기상정보 및 모의 산출 기반의 저류형 비점오염 저감 시스템 및 저류형 비점오염 저감 방법 |
CN103020424A (zh) * | 2012-11-22 | 2013-04-03 | 北京师范大学 | 一种基于人工降雨实验的北方平原农田区非点源污染负荷估算方法 |
CN103544550A (zh) * | 2013-11-08 | 2014-01-29 | 湖南科技大学 | 一种金属矿区土-水界面重金属污染负荷的预测方法 |
CN103810537A (zh) * | 2014-02-12 | 2014-05-21 | 南京大学 | 一种基于水质模型的区域环境风险评估方法 |
CN104361523A (zh) * | 2014-11-06 | 2015-02-18 | 浙江大学 | 一种基于gis的分布式稻田氮素径流流失负荷估算方法 |
CN104732069A (zh) * | 2015-02-27 | 2015-06-24 | 中国水利水电科学研究院 | 一种变入河系数计算污染物入河量的方法 |
CN107066808A (zh) * | 2017-02-28 | 2017-08-18 | 西北农林科技大学 | 一种丘陵区非点源氮磷流失形态构成分布式模拟方法 |
CN107368108A (zh) * | 2017-07-24 | 2017-11-21 | 中国科学院测量与地球物理研究所 | 田沟塘一体化联合调控稻田面源污染的方法 |
KR101846438B1 (ko) * | 2017-08-08 | 2018-04-06 | 대한민국 | 유역모델을 이용한 오염물질 배출경로 재현방법 |
-
2018
- 2018-05-21 CN CN201810486632.1A patent/CN108717453B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008210362A (ja) * | 2007-02-27 | 2008-09-11 | Inha-Industry Partnership Inst | Gisを利用した水質汚染負荷算定用の意思決定支援システム及びその運営方法 |
KR20120023259A (ko) * | 2010-09-01 | 2012-03-13 | 한국건설기술연구원 | 기상정보 및 모의 산출 기반의 저류형 비점오염 저감 시스템 및 저류형 비점오염 저감 방법 |
CN103020424A (zh) * | 2012-11-22 | 2013-04-03 | 北京师范大学 | 一种基于人工降雨实验的北方平原农田区非点源污染负荷估算方法 |
CN103544550A (zh) * | 2013-11-08 | 2014-01-29 | 湖南科技大学 | 一种金属矿区土-水界面重金属污染负荷的预测方法 |
CN103810537A (zh) * | 2014-02-12 | 2014-05-21 | 南京大学 | 一种基于水质模型的区域环境风险评估方法 |
CN104361523A (zh) * | 2014-11-06 | 2015-02-18 | 浙江大学 | 一种基于gis的分布式稻田氮素径流流失负荷估算方法 |
CN104732069A (zh) * | 2015-02-27 | 2015-06-24 | 中国水利水电科学研究院 | 一种变入河系数计算污染物入河量的方法 |
CN107066808A (zh) * | 2017-02-28 | 2017-08-18 | 西北农林科技大学 | 一种丘陵区非点源氮磷流失形态构成分布式模拟方法 |
CN107368108A (zh) * | 2017-07-24 | 2017-11-21 | 中国科学院测量与地球物理研究所 | 田沟塘一体化联合调控稻田面源污染的方法 |
KR101846438B1 (ko) * | 2017-08-08 | 2018-04-06 | 대한민국 | 유역모델을 이용한 오염물질 배출경로 재현방법 |
Non-Patent Citations (9)
Title |
---|
LI FANG ET AL.: ""Pollution load simulation of Dongting Lake basin based on SWAT and GIS"", 《 2015 23RD INTERNATIONAL CONFERENCE ON GEOINFORMATICS》 * |
SHIGEYA MAEDA ET AL.: ""A GIS-aided two-phase grey fuzzy optimization model for nonpoint source pollution control in a small watershed"", 《PADDY AND WATER ENVIRONMENT》 * |
YITING CHEN ET AL.: ""Estimate model of non-point source pollution load in plain river-net area: A case study in Dafeng city"", 《2011 INTERNATIONAL CONFERENCE ON ELECTRICAL AND CONTROL ENGINEERING》 * |
左一鸣 等: ""基于GIS的平原河网非恒定流计算模型"", 《水利水运工程学报》 * |
张荣保: ""典型平原河网地区污染负荷模型研究"", 《中国优秀硕士学位论文全文数据库》 * |
李卉 等: ""平原河网地区农业非点源污染负荷估算方法综述"", 《北京师范大学学报(自然科学版)》 * |
王婷婷 等: ""洞庭湖平原堤垸区非点源污染模拟与分析"", 《水力发电学报》 * |
管仪庆 等: ""平原河网地区水环境模拟及污染负荷计算"", 《水资源保护》 * |
蒋艳 等: ""雁栖河流域点源氮磷污染负荷量的计算与分析"", 《中国水利水电科学研究院学报》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110287615A (zh) * | 2019-06-28 | 2019-09-27 | 南京大学 | 一种基于遥感解译和降雨实验的雨水径流污染负荷测算方法 |
CN110826920A (zh) * | 2019-11-08 | 2020-02-21 | 河北省南运河河务管理处 | 海岸带污染物负荷估算方法 |
CN111077280A (zh) * | 2020-01-14 | 2020-04-28 | 浙江清华长三角研究院 | 一种基于河网的农村污水处理设施与水质监测站点间的溯源分析方法 |
CN111310124A (zh) * | 2020-02-13 | 2020-06-19 | 北京市环境保护科学研究院 | 一种城市径流数据处理方法 |
CN111523088A (zh) * | 2020-04-13 | 2020-08-11 | 杭州领见数据科技有限公司 | 一种基于dpsir模型的生态环境评价方法 |
CN111639707A (zh) * | 2020-05-29 | 2020-09-08 | 河南大学 | 一种基于gis技术的土地污染控制方法 |
CN113111530A (zh) * | 2021-04-23 | 2021-07-13 | 中国水利水电科学研究院 | 基于分布式水文模型的矿山污染物扩散入河量估算方法 |
CN113256469A (zh) * | 2021-05-14 | 2021-08-13 | 郑州大学环境技术咨询工程有限公司 | 一种城镇地表径流面源污染物入河量核算方法 |
CN113361114A (zh) * | 2021-06-11 | 2021-09-07 | 中国科学院精密测量科学与技术创新研究院 | 一种基于径流路径的多尺度面源污染物入河系数测算方法 |
CN113361114B (zh) * | 2021-06-11 | 2022-05-17 | 中国科学院精密测量科学与技术创新研究院 | 一种基于径流路径的多尺度面源污染物入河系数测算方法 |
CN113552305A (zh) * | 2021-07-21 | 2021-10-26 | 中国水利水电科学研究院 | 一种多孔纤维材料埋设对产流水质影响的机理识别方法 |
CN115629190A (zh) * | 2022-12-08 | 2023-01-20 | 中农创达(北京)环保科技有限公司 | 农业面源污染数据管理方法、装置、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN108717453B (zh) | 2020-05-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108717453A (zh) | 一种基于gis平台的平原河网区污染负荷计算方法 | |
CN108664647B (zh) | 一种集成水环境模型的流域精细化管理*** | |
Schaffner et al. | Modeling the contribution of point sources and non-point sources to Thachin River water pollution | |
CN108733915B (zh) | 基于栅格化处理的平原河网区河道面源污染负荷确定方法 | |
CN105973207B (zh) | 大型浅水湖泊水华集聚水域的补水控藻方法 | |
Pai et al. | Identifying priority subwatersheds in the Illinois River drainage area in Arkansas watershed using a distributed modeling approach | |
CN106326405A (zh) | 一种水库流域生态保护方法 | |
CN104143048B (zh) | 一种农田面源污染物入河量的计算方法及装置 | |
Cebe et al. | Monitoring and modeling land-based marine pollution | |
Darracq et al. | Nutrient transport scenarios in a changing Stockholm and Mälaren valley region, Sweden | |
CN112765531A (zh) | 一种流域非点源排放量自动计算方法和*** | |
CN111291305B (zh) | 基于一维水质模型计算流域污染源入河系数的方法 | |
Yang et al. | Pollution load estimation and control countermeasures of Zhangze reservoir | |
CN101858107A (zh) | 雨水处理和利用设施的设计方法 | |
CN116562051B (zh) | 一种陆海氮磷负荷归趋估算方法 | |
Suarez et al. | Evaluation of a coupled hydrodynamic-closed ecological cycle approach for modelling dissolved oxygen in surface waters | |
Cau et al. | Assessment of alternative land management practices using hydrological simulation and a decision support tool: Arborea agricultural region, Sardinia | |
Balcerzak | The Protection of Reservoir Water against the Eutrophication Process. | |
Liu et al. | Dynamic phosphorus budget for lake-watershed ecosystems | |
Faksomboon | Development of a Hydrodynamic Model for Regulating Water Drainage of Reservoir and Water Resources Management, Lamtakong Watershed of Thailand | |
CN109516569A (zh) | 一种黑臭河综合治理方法 | |
Wu et al. | Water environment planning for the Xianghe Segment of China’s Grand Canal | |
Ruskule | Water quality model description | |
Chinh et al. | Estimation of nitrogen transformation in cau river basin, Vietnam using a material flow analysis model | |
CN116542171A (zh) | 基于机器学习水文模型的叶绿素a含量预测方法及装置 |
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 | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20181030 Assignee: Suzhou Sucheng Environmental Technology Co.,Ltd. Assignor: HOHAI University Contract record no.: X2021320000076 Denomination of invention: A calculation method of pollution load in plain river network area based on GIS platform Granted publication date: 20200505 License type: Common License Record date: 20210831 |