一种临河村落非点源污染负荷的估算方法技术

本发明专利技术涉及一种临河村落非点源污染负荷的估算方法，该估算方法主要包括：1)基础资料收集：获取某地小雨、中雨、大到暴雨等三种情况的多年累积平均雨量R1、R2、R3，以及村庄面积A等基本情况；2)野外实地的监测：对有效降雨场次的降雨历时、降雨量和降雨径流的次径流系数、径流历时、产流面积进行监测，收集的水样进行化学需氧量、总氮、总磷的监测；3)非点源污染负荷的模型估算：采用公式计算临河村落非点源污染物化学需氧量、总氮、总磷的年污染负荷。本发明专利技术提出的计算方法更具有指导意义，便于开展相关的控制和管理措施。

【技术实现步骤摘要】

本专利技术属于农业非点源污染负荷量的估算领域，具体涉及到一种确定临河村落非点源污染负荷的方法。
技术介绍
近年来，随着人们对环境问题的关注，一类十分普遍而又不为人们所熟悉的环境污染问题逐渐得到各国政府环境保护部门的高度重视，这就是非点源污染，国内有的文献或媒体中又称之为面源污染污染，通常在比较正式的、学术性较强的文献中多称“非点源污 m^K ο与点源污染相比，非点源污染起源于分散、多样的地区，地理边界和发生位置难以识别和确定，随机性强、成因复杂、潜伏周期长，因而防治十分困难。随着各国政府对点源污染控制的重视，点源污染在包括我国在内的许多国家已经得到较好的控制和治理，而非点源污染，由于涉及范围广、控制难度大，目前已成为影响水体环境质量的重要污染源。目前，非点源污染主要源头是农业非点源污染和城市地表径流污染两部分，以下将具体介绍在此两方面相关的污染负荷的计算方法第一种是农业非点源污染负荷计算方法在整个农业非点源污染负荷估算领域，主要有以下四种研究方法1)断面实测总负荷减去点源负荷的方法该方法原理简单，操作方便；主要问题是点源估算难以准确。2)水文估算法该方法将枯水期的基流代表点源负荷，而把超过基流时的流量则视作同时包含有点源和非点源污染的情况，两者相减得到非点源污染负荷；主要问题是，当基流中污染物含量较高时，将难以区分点源与非点源之间的差别。3)输出系数法该方法利用污染物输出系数来估算流域非点源污染输出负荷，直接建立土地利用与受纳水体非点源负荷之间的关系，简单实用；但也因此割裂了单位源和受纳水体之间的水力联系，对于干旱缺水的北方地区并不适用。4)模型估算法目前，我国应用较多的是SWAT模型和AGNPS模型。SWAT模型是连续的，适用于几千km2的流域，由气候、水文、泥沙、养分、农药等模块组成，用于预测大型流域管理措施对水质、泥沙和化学物质的作用。AGNPS模型用于暴雨事件，1 200km2流域，由水文模块、侵蚀沉积模块和化学物质模块组成，用于模拟流域管理措施对径流和泥沙产生的影响。第二种是城市地表径流负荷估算方法3城市地表径流污染负荷计算方法一般分为以下三类1)浓度法浓度法是根据地表径流的水质水量同步监测数据计算负荷的方法，即年污染负荷 Ly = 0. OlEMCs X RX Cv，式中EMCs为多场降雨的次降雨平均浓度的均值，R为多年平均降雨量，Cv为径流系数；浓度法从概念出发，在计算中未加考虑降雨产流有效场次的影响，直接应用于临河村落的非点源污染负荷估算不够准确，尤其临河坡地多以透水地面为主，雨强的不同会造成径流系数和产流面积较大范围的波动，大大影响了其准确性。2)统计模型统计模型则用于估算城市开发区污染物排放量，采用分析大量的实测资料，直接建立污染负荷与有关影响因素相关关系的方法，模型为Lt =0.01[(CF)(W(P)(C)],式中Lt为计算时段(t)内径流排放污染负荷；CF为用于对不产生地表径流的降雨进行校正的因子 (产生径流的降雨事件占总降雨事件的比例)滑为径流区内平均径流系数；P为计算时段 (t)内降雨量；C为污染物的径流量加权平均浓度；在我国由于缺乏相关的统计数据，应用统计模型进行相关的污染负荷估算不具有可行性。三、概念模型以污染形成的动力学为基础，分别从污染物的累积、冲刷、迁移过程建立数学模型，从而得出污染物的输出。应用较多的是SWMM模型，后者在城市区域内排水区和排水管网的非点源污染负荷计算方面有着较明显的优势。缺点是其对数据和人员的要求较高，整个模型比较复杂；且鉴于我国村落排水管网铺设上的匮乏，应用此类模型的条件极不满足。目前在太湖、滇池、三峡等流域已开展农村生活污染的调查，但也仅限于由调查结果进行潜在污染负荷的估算，在临河村落这样较小的尺度上，尚无系统的计算公式用于非点源污染负荷的核算。无法实现在村落尺度上的非点源污染负荷的核算，这限制了当前的研究。针对现有的问题，本专利技术根据临河村落的临河坡地多以透水地面为主的特性，按照雨强的不同，对降雨产流情况进行分类，综合考虑了径流系数、产流面积的变化，同时考虑了降雨产流有效场次的影响，对现有方法进行改进，最终提出了本专利技术的数学模型。
技术实现思路
本专利技术的目的是提供一种针对临河村落临河坡地多以透水地面为主，不同雨强下径流系数与实际产流面积会有较大的变化的估算临河村落非点源污染负荷的方法。本专利技术提供了，包括以下步骤1)基础资料收集从当地气象部门或相关文献中获取该地小雨、中雨、大到暴雨的多年累积平均雨量Rl、R2、R3，以及村庄面积A等基本情况；2)野外实地的监测通过野外现场考察，掌握临河坡地基本的降雨产流状况，选定2 3个产流小区的出口作为监测点，对有效降雨场次进行降雨径流的监测，统计全年的有效降雨场次，得到CFi ；同时，在选定的监测点监测降雨历时、降雨量、次径流系数、径流历时、产流面积；收集的水样应依照预先制定的保存方案和监测分析方案，完成化学需氧量、总氮和总磷的监测分析，然后利用EMC的计算公式得出次降雨径流的平均浓度EMC(COD)、 EMC(TN)和EMC(TP)；通过对3种降雨情况下，多次的次降雨径流污染物浓度求算术平均值， 从而可得 EMCsi (COD)、EMCsi (TN)、EMCsi (TP)；M TCtQtAt其中EMC的计算公式为EMC = —= ^q ^，式中M表示COD、TN、TP的量，单位为g ；V表示总径流量，单位为m3 ;Ct为t时刻 COD、TN、TP的浓度，单位为mg/L ;Qt为t时刻径流流量，单位为m3/min ；所述Δ t为采样时间间隔，单位为min ；3)非点源污染负荷的模型估算，根据下面的公式进行估算，即可得到临河村落非点源污染物化学需氧量COD、总氮TN、总磷TP的年污染负荷，3J^Ly = OmY4EMCsi χCn χR1 χCFl χ —式中所述i = 1、2、3，分别表示小雨、中雨、大到暴雨三种情况；所述EMCsi表示在i降雨情况下，多次降雨的次降雨径流浓度EMC的平均值；所述Cvi表示在i降雨情况下，多场降雨的次径流系数平均值，即多场降雨中径流深度与降雨深度的比值的平均值；所述氏表示在i种降雨多年平均累积降雨量；所述CFi表示在i降雨情况下，产生径流的降雨事件占总降雨事件的比例；所述Ai表示在i降雨情况下，平均产流面积；所述A为村庄面积。本专利技术所述的方法中所述小雨一般指日04小时)内降水量小于IOmm的降雨；所述中雨一般指日( 小时)内降水量在10-24. 9mm的降雨；所述大雨一般指日( 小时)内降水量为25-49. 9mm的降雨；所述暴雨一般指日04小时)内降水量为50mm或以上的强降雨；所述有效降雨场次能产生地表径流的降雨场次。所述降雨径流是指由降雨所形成的径流；所述降雨历时(duration of rainfall)是指在一次降雨过程中，从开始到结束所持续的时间为降雨历时，一般以min，h或d计；所述降雨量是指从天空降落到地面上的雨水，未经蒸发、渗透、流失而在水面上积聚的水层深度，称为降雨量，以毫米为单位；所述次径流系数在一次降雨过程中，产生的径流深与降雨深的比值，取值范围为 0 1 ；所述径流历时在一次降雨过程中，从开始产生径流到径流消失所持续的时间为径流历时，以m本文档来自技高网...

【技术保护点】
１．一种临河村落非点源污染负荷的估算方法，包括以下步骤：１）基础资料收集：从当地气象部门或相关文献中获取该地小雨、中雨、大到暴雨三种情况的多年累积平均雨量Ｒ１、Ｒ２、Ｒ３，以及村庄面积Ａ的基本情况；２）野外实地的监测：通过野外现场考察，掌握临河坡地基本的降雨产流状况，然后选定２～３个产流小区的出口作为监测点，对降雨进行降雨径流的监测，然后统计全年的有效降雨场次，得到ＣＦｉ；同时，在选定的监测点监测降雨历时、降雨量、次径流系数、径流历时、产流面积；收集的水样应依照预先制定的保存方案和监测分析方案，完成化学需氧量、总氮、总磷的监测分析，然后利用ＥＭＣ的计算公式得出次降雨径流的平均浓度ＥＭＣ（ＣＯＤ）、ＥＭＣ（ＴＮ）、ＥＭＣ（ＴＰ）；通过对３种降雨情况下，多次的次降雨径流污染物浓度求算术平均值，从而可得ＥＭＣＳｉ（ＣＯＤ）、ＥＭＣＳｉ（ＴＮ）、ＥＭＣＳｉ（ＴＰ）；其中ＥＭＣ的计算公式为（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｉ）ＥＭＣ（／ｍｉ）?（ｍｏ）＝（／ｍｏ）?（ｍｆｒａｃ）?（ｍｉ）Ｍ（／ｍｉ）?（ｍｉ）Ｖ（／ｍｉ）?（／ｍｆｒａｃ）?（ｍｏ）＝（／ｍｏ）?（ｍｆｒａｃ）?（ｍｒｏｗ）?（ｍｉ）＆Ｓｉｇｍａ；（／ｍｉ）?（ｍｓｕｂ）?（ｍｉ）Ｃ（／ｍｉ）?（ｍｉ）ｔ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｓｕｂ）?（ｍｉ）Ｑ（／ｍｉ）?（ｍｉ）ｔ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｉ）＆Ｄｅｌｔａ；ｔ（／ｍｉ）?（／ｍｒｏｗ）?（ｍｒｏｗ）?（ｍｉ）＆Ｓｉｇｍａ；（／ｍｉ）?（ｍｓｕｂ）?（ｍｉ）Ｑ（／ｍｉ）?（ｍｉ）ｔ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｉ）＆Ｄｅｌｔａ；ｔ（／ｍｉ）?（／ｍｒｏｗ）?（／ｍｆｒａｃ）?（ｍｏ），（／ｍｏ）?（／ｍｒｏｗ）?（／ｍａｔｈ）所述Ｍ表示ＣＯＤ、ＴＮ、ＴＰ的量，单位为ｇ；所述Ｖ表示总径流量，单位为ｍ３；所述Ｃｔ为ｔ时刻ＣＯＤ、ＴＮ、ＴＰ的浓度，单位为ｍｇ／Ｌ；所述Ｑｔ为ｔ时刻径流流量，单位为ｍ３／ｍｉｎ；所述Δｔ为采样时间间隔，单位为ｍｉｎ；３）非点源污染负荷的模型估算：根据下面的公式进行估算，即得临河村落非点源污染物化学需氧量、总氮、总磷的年污染负荷，（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｓｕｂ）?（ｍｉ）Ｌ（／ｍｉ）?（ｍｉ）Ｙ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ）＝（／ｍｏ）?（ｍｎ）０．０１（／ｍｎ）?（ｍｕｎｄｅｒｏｖｅｒ）?（ｍｉ）＆Ｓｉｇｍａ；（／ｍｉ）?（ｍｒｏｗ）?（ｍｉ）ｉ（／ｍｉ）?（ｍｏ）＝（／ｍｏ）?（ｍｎ）１（／ｍｎ）?（／ｍｒｏｗ）?（ｍｎ）３（／ｍｎ）?（／ｍｕｎｄｅｒｏｖｅｒ）?（ｍｉ）ＥＭ（／ｍｉ）?（ｍｓｕｂ）?（ｍｉ）Ｃ（／ｍｉ）?（ｍｉ）Ｓｉ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ）＆ｔｉｍｅｓ；（／ｍｏ）?（ｍｓｕｂ）?（ｍｉ）Ｃ（／ｍｉ）?（ｍｉ）Ｖｉ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ）＆ｔｉｍｅｓ；（／ｍｏ）?（ｍｓｕｂ）?（ｍｉ）Ｒ（／ｍｉ）?（ｍｉ）ｉ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ）＆ｔｉｍｅｓ；（／ｍｏ）?（ｍｓｕｂ）?（ｍｉ）Ｃ（／ｍｉ）?（ｍｉ）Ｆｉ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ）＆ｔｉｍｅｓ；（／ｍｏ）?（ｍｆｒａｃ）?（ｍｓｕｂ）?（ｍｉ）Ａ（／ｍｉ）?（ｍｉ）ｉ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｉ）Ａ（／ｍｉ）?（／ｍｆｒａｃ）?（／ｍｒｏｗ）?（／ｍａｔｈ）式中：所述ｉ＝１、２、３，分别表示小雨、中雨、大到暴雨三种情况；所述ＥＭＣＳｉ表示在ｉ降雨情况下，多次降雨的次降雨径流浓度ＥＭＣ的平均值；所述ＣＶｉ表示在ｉ降雨情况下，多场降雨的次径流系数平均值，即多场降雨中径流深度与降雨深度的比值的平均值；所述Ｒｉ表示在ｉ种降雨多年平均累积降雨量；所述ＣＦｉ表示在ｉ降雨情况下，产生径流的降雨事件占总降雨事件的比例；所述Ａｉ表示在ｉ降雨情况下，平均产流面积；所述Ａ为村庄面积。...

【技术特征摘要】
1.一种临河村落非点源污染负荷的估算方法，包括以下步骤1)基础资料收集从当地气象部门或相关文献中获取该地小雨、中雨、大到暴雨三种情况的多年累积平均雨量队、&、R3,以及村庄面积A的基本情况；2)野外实地的监测通过野外现场考察，掌握临河坡地基本的降雨产流状况，然后选定2 3个产流小区的出口作为监测点，对降雨进行降雨径流的监测，然后统计全年的有效降雨场次，得到CFi ；同时，在选定的监测点监测降雨历时、降雨量、次径流系数、径流历时、 产流面积；收集的水样应依照预先制定的保存方案和监测分析方案，完成化学需氧量、总氮、总磷的监测分析，然后利用EMC的计算公式得出次降雨径流的平均浓度EMC(COD)、EMC(TN)、 EMC(TP)；通过对3种降雨情况下，多次的次降雨径流污染物浓度求算术平均值，从而可得 EMCsi (COD)、EMCsi (TN)、EMCsi (TP)；M TCtQtAt其中EMC的计算公式为EMC = -= ^q ^，所述M表示COD、TN、TP的量，单位为g ；所述V表示总径流量，单位为m3 ；所述Ct为t 时刻CO...

【专利技术属性】
技术研发人员：陈家军，孙添伟，王浩，史震天，
申请(专利权)人：北京师范大学，
类型：发明
国别省市：11