会员体验
专利管家(专利管理)
工作空间(专利管理)
风险监控(情报监控)
数据分析(专利分析)
侵权分析(诉讼无效)
联系我们
交流群
官方交流:
QQ群: 891211   
微信请扫码    >>>
现在联系顾问~
首页 / 专利库 / 污染与排放控制 / 排放量 / 一种农田面源污染县域地表径流总氮排放量预测方法

一种农田面源污染县域地表径流总氮排放量预测方法

阅读:1054发布:2020-07-18

IPRDB可以提供一种农田面源污染县域地表径流总氮排放量预测方法专利检索,专利查询,专利分析的服务。并且本发明提出一种农田面源污染县域地表径流总氮排放量预测方法,包括步骤:1)通过农田面源污染发生分区,获取不同分区中县级行政区划范围,针对一种分区的县级生物化学检测,获得N年县级地表径流总氮排放量;2)选择N‑x年影响地表径流总氮排放量的氮排放量数据,建立与N年地表径流总氮县级预测模型;3)应用地表径流总氮县级预测模型,通过N年影响地表径流总氮排放量的数据,求得第N+x年的地表径流总氮排放量。本发明提出的方法,针对县域农田面源污染地表径流总氮无法预测的难题,在农田面源污染发生区划的基础上,建立地表径流总氮排放量预测模型,预测农田地表径流总氮排放量,为农业面源污染防治提供技术支撑。,下面是一种农田面源污染县域地表径流总氮排放量预测方法专利的具体信息内容。

1.一种农田面源污染县域地表径流总氮排放量预测方法,其特征在于,包括步骤:

1)基于全国各区县地貌、气候、肥料施用自然因素,确定农田面源污染发生区划,确定各农田面源污染地表径流总氮发生分区及各分区涉及县级行政区划范围;

2)基于农田生态养分循环理论,选择影响地表径流总氮排放量SB因素:种子带入途径ZT、化肥施用途径FH、有机肥施用途径FY、生物固氮途径NG、大气降水途径SJ、灌溉带入途径SG、秸秆还田途径JT、作物收获带出途径ZS、地下淋溶流失途径SX、氨挥发途径QH、氧化亚氮排放途径QP和秸秆移除途径JY;

3)在一个农田面源污染发生分区内,基于县级区域的生物化学检测,获得第N年分区内各县的农田面源污染地表径流总氮排放量SB,及第N年和第N-x年地表径流总氮影响因素量;固定步长x为2年;

4)建立县域N年地表径流总氮与N-x年地表径流总氮影响因素预测模型如下所示:

SB=a+b×ZT+c×FH+d×FY+e×NG+f×SJ+g×SG+h×JT+i×ZS+j×SX+k×QH+l×QP+m×JY  (1)

式(1)中,a为常数项,b、c、d、e、f、g、h、i、j、k、l、m为回归系数;影响地表径流总氮排放量的氮排放量数据为第N-x年的数据;

选择该农田面源污染发生分区内所有县级行政区划数据,代入所述数学模型中,用多元统计分析法求得a~m的数值;

5)在所述预测模型中代入步骤4)求得的a~m的数值、第N年的影响地表径流总氮排放量的氮排放量数据,求得第N+x年的地表径流总氮排放量SB。

2.根据权利要求1所述的预测方法,其特征在于,总氮排放量数据中,种子带入途径ZT通过测量县域农田不同种子播种量、不同种子中的含氮量求得;化肥施用途径FH和有机肥施用途径FY根据县域农田不同种植模式面积、化肥与有机肥施用量求得;生物固氮途径NG通过测定豆科植物面积及固氮量测定;大气降水途径SJ通过测定农田年降水量和降水中含氮量确定;灌溉带入途径SG通过测定灌溉量和灌溉用水中含氮量来确定;秸秆还田途径JT、秸秆移除途径JY和作物收获带出途径ZS通过作物播种面积、秸秆还田比例、秸秆移除比例以及测定秸秆和籽粒中的含氮量求得;地表径流总氮排放量SB通过测定地表径流水量和采集地表径流水含磷量求得;地下淋溶流失途径SX通过测定地下淋溶水量和采集淋溶水含氮量求得;氨挥发途径QH与氧化亚氮排放途径QP通过测量不同作物氮肥使用过程中NH3、N2O排放总量确定。

3.根据权利要求1所述的预测方法,其特征在于,所述选择的县级行政区划属于同一个农田面源污染发生分区,所述农田面源污染发生分区为南方湿润区、北方高原区、东北半湿润平原区、西北干旱区中的一种。

说明书全文

一种农田面源污染县域地表径流总氮排放量预测方法

技术领域

[0001] 本发明属于测量领域,具体涉及一种总氮排放量的预测方法。

背景技术

[0002] 随着我国农业和农村经济的快速发展,农业面源污染问题日益突出,特别是在集约化农区,农业面源污染已成为农田生态环境质量下降、地表水富营养化、地下水硝酸盐污染、农田生产能力降低的重要原因。根据2007年开展的第一次全国污染源普查结果,农业源氮、磷年排放量已达270.46万吨和28.47万吨,分别占全国氮磷排放总量的57.2%和67.4%。
[0003] 目前用于农田面源污染氮磷排放预测研究方法主要有面源污染物理模型法。面源污染物理模型是基于对污染物产生、迁移过程进行模拟的机理模型,目前应用较多的有SWAT、AGNPS、HSPF、ANSWERS等。面源污染物理模型法不仅可以估算污染物输出负荷,模拟污染物产生、迁移、转化过程,并可对面源污染的发展进行预测,但此类模型研究范围以流域为主,结构复杂,基础数据、参数数据需求量大,计算效率低(Tripathi et al.,2003;Beher et al.,2006;Panagopoulos et al.,2011;Shang et al.,2012;Niraula et al.,2012;Niraula et al.,2013;吴春玲,2013;邓欧平等,2013)。基于含氮量数据预测总氮排放量通常是以田块尺度或流域尺度为地理基准,建立生物地球化学循环模型,存在获取数据困难、计算时间太长的缺陷。
[0004] 区县是我国最为稳定的行政区划,准确预测县域尺度的总氮排放量对肥料的使用管理、作物种植均有重要意义。目前,我国还没有县域尺度上农田面源污染预测模型。如何全面、实时、准确地预测县域农田面源污染发展趋势将有利于今后科学防治农业面源污染。

发明内容

[0005] 针对本领域的不足之处,本发明的目的是提供一种农田面源污染县域地表径流总氮排放量预测方法。
[0006] 实现本发明目的的技术方案为:
[0007] 一种农田面源污染县域地表径流总氮排放量预测方法,包括步骤:
[0008] 1)基于全国各区县地貌、气候、肥料施用等自然因素,确定农田面源污染发生区划,确定各农田面源污染地表径流总氮发生分区及各分区涉及县级行政区划范围;
[0009] 2)基于农田生态养分循环理论,选择影响地表径流总氮排放量(SB)因素:种子带入途径(ZT)、化肥施用途径(FH)、有机肥施用途径(FY)、生物固氮途径(NG)、大气降水途径(SJ)、灌溉带入途径(SG)、秸秆还田途径(JT)、作物收获带出途径(ZS)、地下淋溶流失途径(SX)、氨挥发途径(QH)、氧化亚氮排放途径(QP)和秸秆移除途径(JY);
[0010] 3)在一个农田面源污染地表径流排放发生分区内,基于县级区域的生物化学检测,获得第N年分区内各县的农田面源污染造成的地表径流总氮排放量(SB),及第N年和第N-x年地表径流总氮影响因素量;
[0011] 其中,x为1~5的正整数;
[0012] 4)建立县域N年地表径流总氮与N-x年地表径流总氮影响因素预测模型如下所示:SB=a+b×ZT+c×FH+d×FY+e×NG+f×SJ+g×SG+h×JT+i×ZS+j×SX+k×Q H+l×QP+m×
JY            (1)
[0013] 式(1)中,a为常数项,b、c、d、e、f、g、h、i、j、k、l、m为回归系数;影响地表径流总氮排放量的氮排放量数据为第N-x年的数据;
[0014] 选择该农田面源污染发生分区内所有县级行政区划数据,代入所述数学模型中,用多元统计分析法求得a~m的数值;
[0015] 5)在所述数学模型中代入步骤2)求得的a~m的数值、第N年的影响地表径流总氮排放量的氮排放量数据,求得第N+x年的地表径流总氮排放量(SB)。
[0016] 预测数学模型的步长x可以是1年到5年甚至更长,但超过5年,种植模式、农药种类等会有较大变化导致预测准确度下降。进一步优选地,所述x为1~3的正整数。
[0017] 其中,所述的总氮排放量数据中,种子带入途径(ZT)通过测量县域农田不同种子播种量、不同种子中的含氮量求得;化肥施用途径(FH)和有机肥施用途径(FY)根据县域农田不同种植模式面积、化肥与有机肥施用量求得;生物固氮途径(NG)通过测定豆科植物面积及固氮量测定;大气降水途径(SJ)通过测定农田年降水量和降水中含氮量确定;灌溉带入途径(SG)通过测定灌溉量和灌溉用水中含氮量来确定;秸秆还田途径(JT)、秸秆移除途径(JY)和作物收获带出途径(ZS)通过作物播种面积、秸秆还田比例、秸秆移除比例以及测定秸秆和籽粒中的含氮量求得;地表径流流失途径(SB)通过测定地表径流水量和采集地表径流水含磷量求得;地下淋溶流失途径(SX)通过测定地下淋溶水量和采集淋溶水含氮量求得;氨挥发途径(QH)与氧化亚氮排放途径(QP)通过测量不同作物氮肥使用过程中NH3、N2O排放总量确定。
[0018] 测定地下淋溶水量可采用本申请单位提出的地下淋溶原位收集装置经测定。
[0019] 其中,所述选择的县级行政区划应属于同一个农田面源污染发生分区,所述农田面源污染发生分区应为南方湿润区、南方湿润丘陵区、北方高原区、东北半湿润平原区、西北干旱区、西北干旱半干旱区中的一种。
[0020] 本发明的有益效果在于:
[0021] 针对县域农田面源污染地表径流总氮排放量无法预测的难题,在农田面源污染发生区划的基础上,基于农业生态学物质循环理论,应用马尔科夫方法,提出地表径流总氮影响因素,建立地表径流总氮排放量预测模型,预测农田地表径流总氮排放量,为农业面源污染防治提供支撑。

附图说明

[0022] 图1农田面源污染途径示意图。
[0023] 图中,1是种子带入途径、2是化肥施用途径、3是有机肥施用途径、4是生物固氮途径、5是大气降水途径、6是灌溉带入途径、7是秸秆还田途径、8是作物收获带出途径、9是地下淋溶流失途径、10是地表径流流失途径、11是氨挥发途径、12是氧化亚氮排放途径和13是秸秆移除途径。

具体实施方式

[0024] 现以以下实施例来说明本发明,但不用来限制本发明的范围。实施例中使用的手段,如无特别说明,均使用本领域常规的手段。
[0025] 实施例1:
[0026] 本实施例选择黄淮海半湿润平原区,进行农田面源污染地表径流总氮排放量预测。
[0027] 黄淮海半湿润平原区包括黄河、淮河、海河流域中下游的北京、天津、河北、山东、河南大部以及苏北、皖北、黄河支流的汾渭盆地和长江流域的南阳盆地,共计636个县,耕地面积2735万公顷,土壤类型以潮土、褐土、棕壤为主。该区地形平坦,肥料投入量高,灌溉条件好。主要考虑肥料和农药淋溶污染,特别是集约化蔬菜种植区;在棉区和露地蔬菜种植区地膜残留污染较为普遍。
[0028] (1)基于全国各区县地貌、气候、肥料投入等农田面源污染发生影响因素,聚类各县域农田面源污染地表径流总氮排放类型,统计黄淮海半湿润平原区内相同类型区县;
[0029] (2)结合统计年鉴数据、文献发表数据、现场调查数据和实地监测数据,核算黄淮海半湿润平原区内区县农田面源污染2012年地表径流总氮排放量(SB)和2010年影响因素氮通量。
[0030] 其中地表径流总氮影响因素包括种子带入途径(ZT)、化肥施用途径(FH)、有机肥施用途径(FY)、生物固氮途径(NG)、大气降水途径(SJ)、灌溉带入途径(SG)、秸秆还田途径(JT)、作物收获带出途径(ZS)、地下淋溶流失途径(SX)、氨挥发途径(QH)、氧化亚氮排放途径(QP)和秸秆移除途径(JY)。其中,所述的氮排放量数据中,种子带入途径(ZT)通过测量县域农田不同种子播种量、不同种子中的含氮量求得;化肥施用途径(FH)和有机肥施用途径(FY)根据县域农田不同种植模式面积、化肥与有机肥施用量求得;生物固氮途径(NG)通过测定豆科植物面积及固氮量测定;大气降水途径(SJ)通过测定农田年降水量和降水中含氮量确定;灌溉带入途径(SG)通过测定灌溉量和灌溉用水中含氮量来确定;秸秆还田途径(JT)、秸秆移除途径(JY)和作物收获带出途径(ZS)通过作物播种面积、秸秆还田比例、秸秆移除比例以及测定秸秆和籽粒中的含氮量求得;地表径流流失途径(SB)通过测定地表径流水量和采集地表径流水含磷量求得;地下淋溶流失途径(SX)通过测定地下淋溶水量和采集淋溶水含氮量求得;氨挥发途径(QH)与氧化亚氮排放途径(QP)通过测量不同作物氮肥使用过程中NH3、N2O排放总量确定。
[0031] 表1 2012年县域地表径流总氮排放量及2010年影响因素(吨)
[0032]
[0033]
[0034] (3)以黄淮海半湿润平原区内所辖区县为统计单元,建立2012年地表径流总氮排放量与2010年影响因素总氮通量数据集。限于篇幅,表1为部分数据,实际输入R统计计算语言软件的区县有212个区县,所有区县均属于黄淮海半湿润平原区。
[0035] (4)应用马尔科夫分析方法,建立地表径流总氮县域预测模型如下:
[0036] SB=a+b×ZT+c×FH+d×FY+e×NG+f×SJ+g×SG+h×JT+i×ZS+j×SX+k×QH+l×QP+m×JY              (1)
[0037] (5)根据各2012年地表径流总氮排放量、各历史2010年份对应的地表径流总氮排放量影响因素(表2),采用多元统计分析法,将表1数据输入R统计计算语言软件,获取总氮排放量预测模型中常数项a的值、回归系数b、c、d、e、f、g、h、i、j、k、l、m的值因此式(1)成为:
[0038] SB=21650+(-0.3775)×ZT+(0.0006705)×FH+(-0.001449)×FY+(0.04089)×NG+(-0.009531)×SJ+(0.04395)×SG+(0.007778)×JT+(0.0005345)×ZS+(0.08889)×SX+(-0.02315)×QH+(0.05624)×QP+(0.02278)×JY
[0039] 回归分析的拟合度R2=0.9219。
[0040] (6)根据已知a的值、回归系数b、c、d、e、f、g、h、i、j、k、l、m的值和式(1)预测得到2014年黄淮海半湿润平原区农田面源污染总氮排放量(表3),为农田面源污染防治提供技术支撑。
[0041] 表2 2012年县域地表径流总氮排放量影响因素(吨)
[0042]
[0043]
[0044] 表3 2014年县域地表径流总氮预测排放量(吨)
[0045]区县代码 2014年SB
210113 113
210114 69
210122 221
210123 170
210124 211
210181 408
210211 28
210212 34
210213 91
[0046] 以上的实施例仅仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通工程技术人员对本发明的技术方案作出的各种变型和改进,均应落入本发明的权利要求书确定的保护范围内。
高效检索全球专利

IPRDB是专利检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,专利查询、专利分析

电话:13651749426

侵权分析

IPRDB的侵权分析产品是IPRDB结合多位一线专利维权律师和专利侵权分析师的智慧,开发出来的一款特色产品,也是市面上唯一一款帮助企业研发人员、科研工作者、专利律师、专利分析师快速定位侵权分析的产品,极大的减少了用户重复工作量,提升工作效率,降低无效或侵权分析的准入门槛。

立即试用