CN104966128B - 一种农田面源污染县域地表径流总氮排放量预测方法 - Google Patents
一种农田面源污染县域地表径流总氮排放量预测方法 Download PDFInfo
- Publication number
- CN104966128B CN104966128B CN201510303181.XA CN201510303181A CN104966128B CN 104966128 B CN104966128 B CN 104966128B CN 201510303181 A CN201510303181 A CN 201510303181A CN 104966128 B CN104966128 B CN 104966128B
- Authority
- CN
- China
- Prior art keywords
- approach
- total nitrogen
- rainwash
- nitrogen
- agricultural non
- 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.)
- Active
Links
Landscapes
- Fertilizers (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Hydroponics (AREA)
- Processing Of Solid Wastes (AREA)
Abstract
本发明提出一种农田面源污染县域地表径流总氮排放量预测方法,包括步骤:1)通过农田面源污染发生分区,获取不同分区中县级行政区划范围,针对一种分区的县级生物化学检测,获得N年县级地表径流总氮排放量;2)选择N‑x年影响地表径流总氮排放量的氮排放量数据,建立与N年地表径流总氮县级预测模型;3)应用地表径流总氮县级预测模型,通过N年影响地表径流总氮排放量的数据,求得第N+x年的地表径流总氮排放量。本发明提出的方法,针对县域农田面源污染地表径流总氮无法预测的难题,在农田面源污染发生区划的基础上,建立地表径流总氮排放量预测模型,预测农田地表径流总氮排放量,为农业面源污染防治提供技术支撑。
Description
技术领域
本发明属于测量领域,具体涉及一种总氮排放量的预测方法。
背景技术
随着我国农业和农村经济的快速发展,农业面源污染问题日益突出,特别是在集约化农区,农业面源污染已成为农田生态环境质量下降、地表水富营养化、地下水硝酸盐污染、农田生产能力降低的重要原因。根据2007年开展的第一次全国污染源普查结果,农业源氮、磷年排放量已达270.46万吨和28.47万吨,分别占全国氮磷排放总量的57.2%和67.4%。
目前用于农田面源污染氮磷排放预测研究方法主要有面源污染物理模型法。面源污染物理模型是基于对污染物产生、迁移过程进行模拟的机理模型,目前应用较多的有SWAT、AGNPS、HSPF、ANSWERS等。面源污染物理模型法不仅可以估算污染物输出负荷,模拟污染物产生、迁移、转化过程,并可对面源污染的发展进行预测,但此类模型研究范围以流域为主,结构复杂,基础数据、参数数据需求量大,计算效率低(Tripathi et al.,2003;Beheret al.,2006;Panagopoulos et al.,2011;Shang et al.,2012;Niraula et al.,2012;Niraula et al.,2013;吴春玲,2013;邓欧平等,2013)。基于含氮量数据预测总氮排放量通常是以田块尺度或流域尺度为地理基准,建立生物地球化学循环模型,存在获取数据困难、计算时间太长的缺陷。
区县是我国最为稳定的行政区划,准确预测县域尺度的总氮排放量对肥料的使用管理、作物种植均有重要意义。目前,我国还没有县域尺度上农田面源污染预测模型。如何全面、实时、准确地预测县域农田面源污染发展趋势将有利于今后科学防治农业面源污染。
发明内容
针对本领域的不足之处,本发明的目的是提供一种农田面源污染县域地表径流总氮排放量预测方法。
实现本发明目的的技术方案为:
一种农田面源污染县域地表径流总氮排放量预测方法,包括步骤:
1)基于全国各区县地貌、气候、肥料施用等自然因素,确定农田面源污染发生区划,确定各农田面源污染地表径流总氮发生分区及各分区涉及县级行政区划范围;
2)基于农田生态养分循环理论,选择影响地表径流总氮排放量(SB)因素:种子带入途径(ZT)、化肥施用途径(FH)、有机肥施用途径(FY)、生物固氮途径(NG)、大气降水途径(SJ)、灌溉带入途径(SG)、秸秆还田途径(JT)、作物收获带出途径(ZS)、地下淋溶流失途径(SX)、氨挥发途径(QH)、氧化亚氮排放途径(QP)和秸秆移除途径(JY);
3)在一个农田面源污染地表径流排放发生分区内,基于县级区域的生物化学检测,获得第N年分区内各县的农田面源污染造成的地表径流总氮排放量(SB),及第N年和第N-x年地表径流总氮影响因素量;
其中,x为1~5的正整数;
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)
式(1)中,a为常数项,b、c、d、e、f、g、h、i、j、k、l、m为回归系数;影响地表径流总氮排放量的氮排放量数据为第N-x年的数据;
选择该农田面源污染发生分区内所有县级行政区划数据,代入所述数学模型中,用多元统计分析法求得a~m的数值;
5)在所述数学模型中代入步骤2)求得的a~m的数值、第N年的影响地表径流总氮排放量的氮排放量数据,求得第N+x年的地表径流总氮排放量(SB)。
预测数学模型的步长x可以是1年到5年甚至更长,但超过5年,种植模式、农药种类等会有较大变化导致预测准确度下降。进一步优选地,所述x为1~3的正整数。
其中,所述的总氮排放量数据中,种子带入途径(ZT)通过测量县域农田不同种子播种量、不同种子中的含氮量求得;化肥施用途径(FH)和有机肥施用途径(FY)根据县域农田不同种植模式面积、化肥与有机肥施用量求得;生物固氮途径(NG)通过测定豆科植物面积及固氮量测定;大气降水途径(SJ)通过测定农田年降水量和降水中含氮量确定;灌溉带入途径(SG)通过测定灌溉量和灌溉用水中含氮量来确定;秸秆还田途径(JT)、秸秆移除途径(JY)和作物收获带出途径(ZS)通过作物播种面积、秸秆还田比例、秸秆移除比例以及测定秸秆和籽粒中的含氮量求得;地表径流流失途径(SB)通过测定地表径流水量和采集地表径流水含磷量求得;地下淋溶流失途径(SX)通过测定地下淋溶水量和采集淋溶水含氮量求得;氨挥发途径(QH)与氧化亚氮排放途径(QP)通过测量不同作物氮肥使用过程中NH3、N2O排放总量确定。
测定地下淋溶水量可采用本申请单位提出的地下淋溶原位收集装置经测定。
其中,所述选择的县级行政区划应属于同一个农田面源污染发生分区,所述农田面源污染发生分区应为南方湿润区、南方湿润丘陵区、北方高原区、东北半湿润平原区、西北干旱区、西北干旱半干旱区中的一种。
本发明的有益效果在于:
针对县域农田面源污染地表径流总氮排放量无法预测的难题,在农田面源污染发生区划的基础上,基于农业生态学物质循环理论,应用马尔科夫方法,提出地表径流总氮影响因素,建立地表径流总氮排放量预测模型,预测农田地表径流总氮排放量,为农业面源污染防治提供支撑。
附图说明
图1农田面源污染途径示意图。
图中,1是种子带入途径、2是化肥施用途径、3是有机肥施用途径、4是生物固氮途径、5是大气降水途径、6是灌溉带入途径、7是秸秆还田途径、8是作物收获带出途径、9是地下淋溶流失途径、10是地表径流流失途径、11是氨挥发途径、12是氧化亚氮排放途径和13是秸秆移除途径。
具体实施方式
现以以下实施例来说明本发明,但不用来限制本发明的范围。实施例中使用的手段,如无特别说明,均使用本领域常规的手段。
实施例1:
本实施例选择黄淮海半湿润平原区,进行农田面源污染地表径流总氮排放量预测。
黄淮海半湿润平原区包括黄河、淮河、海河流域中下游的北京、天津、河北、山东、河南大部以及苏北、皖北、黄河支流的汾渭盆地和长江流域的南阳盆地,共计636个县,耕地面积2735万公顷,土壤类型以潮土、褐土、棕壤为主。该区地形平坦,肥料投入量高,灌溉条件好。主要考虑肥料和农药淋溶污染,特别是集约化蔬菜种植区;在棉区和露地蔬菜种植区地膜残留污染较为普遍。
(1)基于全国各区县地貌、气候、肥料投入等农田面源污染发生影响因素,聚类各县域农田面源污染地表径流总氮排放类型,统计黄淮海半湿润平原区内相同类型区县;
(2)结合统计年鉴数据、文献发表数据、现场调查数据和实地监测数据,核算黄淮海半湿润平原区内区县农田面源污染2012年地表径流总氮排放量(SB)和2010年影响因素氮通量。
其中地表径流总氮影响因素包括种子带入途径(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排放总量确定。
表1 2012年县域地表径流总氮排放量及2010年影响因素(吨)
(3)以黄淮海半湿润平原区内所辖区县为统计单元,建立2012年地表径流总氮排放量与2010年影响因素总氮通量数据集。限于篇幅,表1为部分数据,实际输入R统计计算语言软件的区县有212个区县,所有区县均属于黄淮海半湿润平原区。
(4)应用马尔科夫分析方法,建立地表径流总氮县域预测模型如下:
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)
(5)根据各2012年地表径流总氮排放量、各历史2010年份对应的地表径流总氮排放量影响因素(表2),采用多元统计分析法,将表1数据输入R统计计算语言软件,获取总氮排放量预测模型中常数项a的值、回归系数b、c、d、e、f、g、h、i、j、k、l、m的值因此式(1)成为:
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
回归分析的拟合度R2=0.9219。
(6)根据已知a的值、回归系数b、c、d、e、f、g、h、i、j、k、l、m的值和式(1)预测得到2014年黄淮海半湿润平原区农田面源污染总氮排放量(表3),为农田面源污染防治提供技术支撑。
表2 2012年县域地表径流总氮排放量影响因素(吨)
表3 2014年县域地表径流总氮预测排放量(吨)
区县代码 | 2014年SB |
210113 | 113 |
210114 | 69 |
210122 | 221 |
210123 | 170 |
210124 | 211 |
210181 | 408 |
210211 | 28 |
210212 | 34 |
210213 | 91 |
以上的实施例仅仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通工程技术人员对本发明的技术方案作出的各种变型和改进,均应落入本发明的权利要求书确定的保护范围内。
Claims (3)
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所述的预测方法,其特征在于,所述选择的县级行政区划属于同一个农田面源污染发生分区,所述农田面源污染发生分区为南方湿润区、北方高原区、东北半湿润平原区、西北干旱区中的一种。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510303181.XA CN104966128B (zh) | 2015-03-27 | 2015-06-04 | 一种农田面源污染县域地表径流总氮排放量预测方法 |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510142408 | 2015-03-27 | ||
CN2015101424087 | 2015-03-27 | ||
CN201510303181.XA CN104966128B (zh) | 2015-03-27 | 2015-06-04 | 一种农田面源污染县域地表径流总氮排放量预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104966128A CN104966128A (zh) | 2015-10-07 |
CN104966128B true CN104966128B (zh) | 2018-04-03 |
Family
ID=54220162
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510303181.XA Active CN104966128B (zh) | 2015-03-27 | 2015-06-04 | 一种农田面源污染县域地表径流总氮排放量预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104966128B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105279324B (zh) * | 2015-10-11 | 2018-07-24 | 中国农业科学院农业资源与农业区划研究所 | 一种基于区域尺度旱地作物氮投入环境阈值的确定方法 |
CN105844366A (zh) * | 2016-04-08 | 2016-08-10 | 中国科学院南京土壤研究所 | 一种区域尺度种植业面源氮损失量统计方法 |
CN108564246B (zh) * | 2018-03-01 | 2021-07-06 | 北京师范大学 | 一种农田土壤氧化亚氮温室气体排放估算方法及装置 |
CN109063904B (zh) * | 2018-07-19 | 2021-11-30 | 中国地质大学(武汉) | 一种流域农业区地下水向地表水输运氮污染量估算方法 |
CN110646252A (zh) * | 2019-09-27 | 2020-01-03 | 沈阳农业大学 | 一种干湿交替稻田用便携式定时多点氨挥发采样装置 |
CN111008789B (zh) * | 2019-12-20 | 2022-11-25 | 生态环境部华南环境科学研究所 | 一种核算种植业氮肥施用后氨排放量的方法 |
CN112232680B (zh) * | 2020-10-19 | 2021-11-09 | 三峡大学 | 基于Copula函数的流域氮代谢环境风险评价方法 |
CN112348241B (zh) * | 2020-10-28 | 2024-01-05 | 中国农业科学院农业资源与农业区划研究所 | 一种保护地农用地面源污染排放量预测方法及系统 |
CN112464434B (zh) * | 2020-10-28 | 2021-12-07 | 中国农业科学院农业资源与农业区划研究所 | 一种灌溉农用地面源污染排放量预测方法及系统 |
CN114936721B (zh) * | 2022-07-15 | 2022-10-28 | 北京师范大学 | 河流氧化亚氮排放量的确定方法、装置和电子设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102087681A (zh) * | 2010-11-30 | 2011-06-08 | 浙江大学 | 预测稻田化肥尿素氮多维通量迁移转化的模型 |
CN104268787A (zh) * | 2014-10-21 | 2015-01-07 | 何文 | 农业氮面源污染量的一种确定技术 |
-
2015
- 2015-06-04 CN CN201510303181.XA patent/CN104966128B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102087681A (zh) * | 2010-11-30 | 2011-06-08 | 浙江大学 | 预测稻田化肥尿素氮多维通量迁移转化的模型 |
CN104268787A (zh) * | 2014-10-21 | 2015-01-07 | 何文 | 农业氮面源污染量的一种确定技术 |
Non-Patent Citations (4)
Title |
---|
三峡库区重庆段农业面源污染负荷的区域分异与预测;肖新成 等;《应用基础与工程科学学报》;20140815;第22卷(第4期);第634-646页 * |
模糊线性回归模型在河流水体总氮浓度预测中的应用;周九州 等;《生态学杂志》;20091215;第28卷(第12期);第2628-2632页 * |
浙江省农业面源污染的影响因素研究;徐建芬;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20140115(第1期);第B027-385页 * |
长三角地区典型稻作流域农田氮平衡及氮负荷研究;杜伟;《中国优秀硕士学位论文全文数据库 农业科技辑》;20120615(第6期);第D043-74页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104966128A (zh) | 2015-10-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104966128B (zh) | 一种农田面源污染县域地表径流总氮排放量预测方法 | |
Yang et al. | Optimising nitrogen fertilisation: A key to improving nitrogen-use efficiency and minimising nitrate leaching losses in an intensive wheat/maize rotation (2008–2014) | |
Khan et al. | Phosphorus efficiency, soil phosphorus dynamics and critical phosphorus level under long-term fertilization for single and double cropping systems | |
Gaydon et al. | Evaluation of the APSIM model in cropping systems of Asia | |
CN104965057B (zh) | 一种农田面源污染县域地表径流总磷排放量的预测方法 | |
Hou et al. | Annual net input fluxes of heavy metals of the agro-ecosystem in the Yangtze River delta, China | |
Ma et al. | Using stable isotopes to determine seasonal variations in water uptake of summer maize under different fertilization treatments | |
Zhao et al. | Maintaining rice yield and reducing N pollution by substituting winter legume for wheat in a heavily-fertilized rice-based cropping system of southeast China | |
Sun et al. | Simulating the fate of nitrogen and optimizing water and nitrogen management of greenhouse tomato in North China using the EU-Rotate_N model | |
Chen et al. | Recovery efficiency and loss of 15N-labelled urea in a rice-soil system under water saving irrigation in the Songnen Plain of Northeast China | |
Li et al. | Regional simulation of nitrate leaching potential from winter wheat-summer maize rotation croplands on the North China Plain using the NLEAP-GIS model | |
CN107220967A (zh) | 一种草地土壤退化评价方法 | |
Yang et al. | Estimation of groundwater use by crop production simulated by DSSAT‐wheat and DSSAT‐maize models in the piedmont region of the North China Plain | |
Krysanova et al. | Assessment of nitrogen leaching from arable land in large river basins: Part I. Simulation experiments using a process-based model | |
Liu et al. | Using the DSSAT model to simulate wheat yield and soil organic carbon under a wheat-maize cropping system in the North China Plain | |
CN104778505B (zh) | 一种农田面源污染县域地下淋溶总氮排放量的预测方法 | |
Bi et al. | Soil phosphorus loss increases under drought-flood abrupt alternation in summer maize planting area | |
Chen et al. | Spatial distribution and dynamics of cotton fine root under film-mulched drip irrigation | |
Ali et al. | Management of salt-affected soils | |
Zhang et al. | Multiple sources characteristics of root water uptake of crop under oasis farmlands in hyper-arid regions | |
Guo et al. | Sustainability evaluation of protected vegetables production in China based on emergy analysis | |
CN104764867B (zh) | 一种农田面源污染县域地下淋溶总磷排放量的预测方法 | |
Ye et al. | Acidification sensitivity and critical loads of acid deposition for surface waters in China | |
Cai et al. | Study on physicochemical properties of high-density sludge sediment in copper mine after phytoremediation | |
Zhang et al. | Soil erosion and its impacts on environment in Yixing tea plantation of Jiangsu Province |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |