CN104764867B - 一种农田面源污染县域地下淋溶总磷排放量的预测方法 - Google Patents

一种农田面源污染县域地下淋溶总磷排放量的预测方法 Download PDF

Info

Publication number
CN104764867B
CN104764867B CN201510142587.4A CN201510142587A CN104764867B CN 104764867 B CN104764867 B CN 104764867B CN 201510142587 A CN201510142587 A CN 201510142587A CN 104764867 B CN104764867 B CN 104764867B
Authority
CN
China
Prior art keywords
total phosphorus
approach
subrosion
county
achieve
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.)
Expired - Fee Related
Application number
CN201510142587.4A
Other languages
English (en)
Other versions
CN104764867A (zh
Inventor
刘申
翟丽梅
刘宏斌
雷秋良
王洪媛
武淑霞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Agricultural Resources and Regional Planning of CAAS
Original Assignee
Institute of Agricultural Resources and Regional Planning of CAAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Agricultural Resources and Regional Planning of CAAS filed Critical Institute of Agricultural Resources and Regional Planning of CAAS
Priority to CN201510142587.4A priority Critical patent/CN104764867B/zh
Publication of CN104764867A publication Critical patent/CN104764867A/zh
Application granted granted Critical
Publication of CN104764867B publication Critical patent/CN104764867B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Fertilizers (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)基于县级区域的生物化学检测,获得第N年和第N-x年全国各县的农田面源污染造成的总磷地下淋溶排放量(SX),x为1~5的正整数;
3)选择影响总磷地下淋溶排放量(SX)的排放量数据:种子带入途径(ZT)、化肥施用途径(FH)、有机肥施用途径(FY)、大气降水途径(SJ)、灌溉带入途径(SG)、秸秆还田途径(JT)、作物收获带出途径(ZS)、地表径流流失途径(SB)、和秸秆移除途径(JY)的磷排放量数据,
建立地下淋溶总磷县域预测模型如下所示:
SX=a+b×ZT+c×FH+d×FY+e×SJ+f×SG+g×JT+h×ZS+i×SB+j×JY (1)
式(1)中,a为常数项,b、c、d、e、f、g、h、i、j为回归系数;等式左边为第N年的地下淋溶总磷排放量,等式右边为影响地下淋溶总磷排放量的磷排放量数据为第N-x年的数据;
选择该农田面源污染发生分区内所有县级行政区划数据,代入所述数学模型中,用多元统计分析法求得a~j的数值;
4)应用马尔科夫分析方法,在所述数学模型中代入步骤2)求得的a~m的数值、第N年的影响地下淋溶总磷排放量的磷排放量数据,求得第N+x年的地下淋溶总磷排放量(SX)。
优选地,所述x为1~3的正整数。
其中,总磷排放量数据中,种子带入途径(ZT)通过测量县域农田不同种子播种量、不同种子中的含磷量求得;化肥施用途径(FH)和有机肥施用途径(FY)根据县域农田不同种植模式面积、化肥与有机肥施用量求得;大气降水途径(SJ)通过测定农田年降水量和降水中含磷量确定;灌溉带入途径(SG)通过测定灌溉量和灌溉用水中含磷量来确定;秸秆还田途径(JT)、秸秆移除途径(JY)和作物收获带出途径(ZS)通过作物播种面积、秸秆还田比例、秸秆移除比例以及测定秸秆和籽粒中的含磷量求得;地下淋溶流失排放量(SX)通过测定地下淋溶水量和采集淋溶水含磷量求得定。
其中,所述选择的县级行政区划属于同一个农田面源污染发生区,所述农田面源污染发生区为南方湿润区、南方湿润丘陵区、北方高原区、东北半湿润平原区、西北干旱区、西北干旱半干旱区中的一种。
本发明的有益效果在于:
针对县域农田面源污染地下淋溶总磷无法预测的难题,在农田面源污染发生区划的基础上,基于农业生态学物质循环理论,提出地下淋溶总磷影响因素,建立地下淋溶总磷排放量预测模型,预测农田地下淋溶总磷排放量,为农业面源污染防治提供支撑。
附图说明
图1农田面源污染途径示意图。
图中,1是种子带入途径、2是化肥施用途径、3是有机肥施用途径、4是大气降水途径、5是灌溉带入途径、6是秸秆还田途径、7是作物收获带出途径、8是地下淋溶流失途径、9是地表径流流失途径和10是秸秆移除途径。
具体实施方式
现以以下实施例来说明本发明,但不用来限制本发明的范围。实施例中使用的手段,如无特别说明,均使用本领域常规的手段。
实施例1:
本实施例选择黄淮海半湿润平原区农田面源污染地下淋溶总磷排放量预测。
黄淮海半湿润平原区包括黄河、淮河、海河流域中下游的北京、天津、河北、山东、河南大部以及苏北、皖北、黄河支流的汾渭盆地和长江流域的南阳盆地,共计636个县,耕地面积2735万ha,土壤类型以潮土、褐土、棕壤为主。该区地形平坦,肥料投入量高,灌溉条件好。主要考虑肥料和农药淋溶污染,特别是集约化蔬菜种植区;在棉区和露地蔬菜种植区地膜残留污染较为普遍。
预测方法包括步骤:
1)基于全国各区县地貌、气候、肥料投入等农田面源污染发生影响因素,聚类各县域农田面源污染地下淋溶总磷排放类型,统计黄淮海半湿润平原区内相同类型区县,统计其中212个县的总磷排放量(SX)和地下淋溶影响因素导致的磷排放量;
其中,总磷排放量数据中,种子带入途径(ZT)通过测量县域农田不同种子播种量、不同种子中的含磷量求得;化肥施用途径(FH)和有机肥施用途径(FY)根据县域农田不同种植模式面积、化肥与有机肥施用量求得;大气降水途径(SJ)通过测定农田年降水量和降水中含磷量确定;灌溉带入途径(SG)通过测定灌溉量和灌溉用水中含磷量来确定;秸秆还田途径(JT)、秸秆移除途径(JY)和作物收获带出途径(ZS)通过作物播种面积、秸秆还田比例、秸秆移除比例以及测定秸秆和籽粒中的含磷量求得;地下淋溶流失排放量(SX)通过测定地下淋溶水量和采集淋溶水含磷量求得定。
其中总磷地下淋溶影响因素包括种子带入途径(ZT)、化肥施用途径(FH)、有机肥施用途径(FY)、大气降水途径(SJ)、灌溉带入途径(SG)、秸秆还田途径(JT)、作物收获带出途径(ZS)、地下淋溶流失途径(SX)和秸秆移除途径(JY),相互关系见图1。
表1 2012年县域地下淋溶总磷排放量及2010年影响因素(吨)
3)应用马尔科夫分析方法,建立地下淋溶总磷县域预测模型如下所示:
SX=a+b×ZT+c×FH+d×FY+e×SJ+f×SG+g×JT+h×ZS+i×SB+j×JY
4)将2012年地下淋溶总磷排放量、2010年份对应的地下淋溶磷排放量各影响因素的数值输入R统计软件,采用多元统计分析法,求得总磷排放量预测模型中常数项a的值、回归系数b、c、d、e、f、g、h、i、j的值。
SX=1606+(0.1313)×ZT+(0.00009661)×FH+(0.001184)×FY+(1.626)×SJ+(0.7193)×SG+(0.1508)×JT+(-0.06694)×ZS+(-0.4051)×SB+(0.1599)×JY
回归分析的拟合度R2=0.7539。
5)基于上式中确定了的常数项a的值、回归系数b、c、d、e、f、g、h、i、j的值,代入2012年地下淋溶磷排放量各影响因素的数值(表2),预测2014年黄淮海半湿润平原区农田面源污染总磷排放量(SB,表3),为农田面源污染防治提供技术支撑。
表2 2012年县域地下淋溶总磷排放量影响因素(吨)
表3 2014年县域地下淋溶总磷预测排放量(吨)
区县代码 2014年SX
210113 1
210114 2
210122 15
210123 12
210124 4
210181 20
210211 2
210212 2
210213 4
以上的实施例仅仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通工程技术人员对本发明的技术方案作出的各种变型和改进,均应落入本发明的权利要求书确定的保护范围内。

Claims (4)

1.一种农田面源污染县域地下淋溶总磷排放量的预测方法,其特征在于,包括步骤:
1)基于全国各区县地貌、气候、肥料施用的自然因素,确定农田面源污染发生区划,确定各农田面源污染地下淋溶总磷发生分区及各分区涉及县级行政区划范围;
2)基于县级区域的生物化学检测,获得第N年和第N-x年全国各县的农田面源污染造成的地下淋溶总磷排放量SX,x为1~5的正整数;
3)选择影响地下淋溶总磷排放量SX的排放量数据:种子带入途径ZT、化肥施用途径FH、有机肥施用途径FY、大气降水途径SJ、灌溉带入途径SG、秸秆还田途径JT、作物收获带出途径ZS、地表径流流失途径SB、和秸秆移除途径JY的磷排放量数据,
建立地下淋溶总磷县域预测模型如下所示:
SX=a+b×ZT+c×FH+d×FY +e×SJ+f×SG+g×JT+h×ZS+i×SB+j×JY (1)
式(1)中,a为常数项,b、c、d、e、f、g、h、i、j为回归系数;等式左边为第N年的地下淋溶总磷排放量,等式右边影响地下淋溶总磷排放量的磷排放量数据为第N-x年的数据;
选择该农田面源污染发生分区内所有县级行政区划数据,代入所述预测模型中,用多元统计分析法求得a~j的数值;
4)应用马尔科夫分析方法,在所述预测模型中代入步骤2)求得的a~j的数值、第N年的影响地下淋溶总磷排放量的磷排放量数据,求得第N+x年的地下淋溶总磷排放量SX。
2.根据权利要求1所述的预测方法,其特征在于,所述x为1~3的正整数。
3.根据权利要求1所述的预测方法,其特征在于,总磷排放量数据中,种子带入途径ZT通过测量县域农田不同种子播种量、不同种子中的含磷量求得;化肥施用途径FH和有机肥施用途径FY根据县域农田不同种植模式面积、化肥与有机肥施用量求得;大气降水途径SJ通过测定农田年降水量和降水中含磷量确定;灌溉带入途径SG通过测定灌溉量和灌溉用水中含磷量来确定;秸秆还田途径JT、秸秆移除途径JY和作物收获带出途径ZS通过作物播种面积、秸秆还田比例、秸秆移除比例以及测定秸秆和籽粒中的含磷量求得;地表径流流失途径SB通过测定地表径流水量和采集地表径流水含磷量求得;地下淋溶总磷排放量SX通过测定地下淋溶水量和采集淋溶水含磷量求得。
4.根据权利要求1所述的预测方法,其特征在于,所述选择的县级行政区划属于同一个农田面源污染发生区,所述农田面源污染发生区为南方湿润区、北方高原区、东北半湿润平原区、西北干旱区中的一种。
CN201510142587.4A 2015-03-27 2015-03-27 一种农田面源污染县域地下淋溶总磷排放量的预测方法 Expired - Fee Related CN104764867B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510142587.4A CN104764867B (zh) 2015-03-27 2015-03-27 一种农田面源污染县域地下淋溶总磷排放量的预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510142587.4A CN104764867B (zh) 2015-03-27 2015-03-27 一种农田面源污染县域地下淋溶总磷排放量的预测方法

Publications (2)

Publication Number Publication Date
CN104764867A CN104764867A (zh) 2015-07-08
CN104764867B true CN104764867B (zh) 2017-05-10

Family

ID=53646828

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510142587.4A Expired - Fee Related CN104764867B (zh) 2015-03-27 2015-03-27 一种农田面源污染县域地下淋溶总磷排放量的预测方法

Country Status (1)

Country Link
CN (1) CN104764867B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107656034B (zh) * 2017-10-11 2020-06-30 中国农业科学院农业环境与可持续发展研究所 农田面源污染对地表径流中总磷排放量影响的预测方法
CN112348241B (zh) * 2020-10-28 2024-01-05 中国农业科学院农业资源与农业区划研究所 一种保护地农用地面源污染排放量预测方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102087681A (zh) * 2010-11-30 2011-06-08 浙江大学 预测稻田化肥尿素氮多维通量迁移转化的模型
CN102628856A (zh) * 2012-02-24 2012-08-08 中国环境科学研究院 一种湖泊营养物基准参照状态确定的组合方法
CN103473463A (zh) * 2013-09-17 2013-12-25 中国环境科学研究院 一种定量确定湖泊流域水体氮磷背景浓度的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102087681A (zh) * 2010-11-30 2011-06-08 浙江大学 预测稻田化肥尿素氮多维通量迁移转化的模型
CN102628856A (zh) * 2012-02-24 2012-08-08 中国环境科学研究院 一种湖泊营养物基准参照状态确定的组合方法
CN103473463A (zh) * 2013-09-17 2013-12-25 中国环境科学研究院 一种定量确定湖泊流域水体氮磷背景浓度的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
东辽河流域土地利用变化对非点源污染的影响研究;孙丽娜 等;《中国环境科学》;20130831;第33卷(第8期);第1459-1467页 *
太湖流域典型蔬菜地磷流失因素的交互效应模型研究;苑韶峰 等;《核农学报》;20130531;第27卷(第5期);第663-672页 *

Also Published As

Publication number Publication date
CN104764867A (zh) 2015-07-08

Similar Documents

Publication Publication Date Title
CN104966128B (zh) 一种农田面源污染县域地表径流总氮排放量预测方法
CN104965057B (zh) 一种农田面源污染县域地表径流总磷排放量的预测方法
Pan et al. Grain yield, water productivity and nitrogen use efficiency of rice under different water management and fertilizer-N inputs in South China
Li et al. Managing irrigation and fertilization for the sustainable cultivation of greenhouse vegetables
Murtaza et al. Irrigation and soil management strategies for using saline-sodic water in a cotton–wheat rotation
Yang et al. Estimation of irrigation requirement for sustainable water resources reallocation in North China
Hou et al. Detection and attribution of nitrogen runoff trend in China's croplands
Chen et al. Green-depressing cropping system: A referential land use practice for fallow to ensure a harmonious human-land relationship in the farming-pastoral ecotone of northern China
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
Zhang et al. Nitrate-nitrogen dynamics and nitrogen budgets in rice-wheat rotations in Taihu Lake region, China
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
CN104756699A (zh) 兼顾环境容量的豫北潮土区夏玉米高产栽培施肥方法
Li et al. Rational biochar application rate for cotton nutrient content, growth, yields, productivity, and economic benefits under film-mulched trickle irrigation
Huang et al. Linking landscape dynamics to the relationship between water purification and soil retention
Zhou et al. Arable land use intensity change in China from 1985 to 2005: evidence from integrated cropping systems and agro economic analysis
Xu et al. Agricultural non-point source pollution management in a reservoir watershed based on ecological network analysis of soil nitrogen cycling
Bi et al. Soil phosphorus loss increases under drought-flood abrupt alternation in summer maize planting area
Ali et al. Management of salt-affected soils
Shi et al. Influence of climate and socio-economic factors on the spatio-temporal variability of soil organic matter: A case study of Central Heilongjiang Province, China
CN104778505B (zh) 一种农田面源污染县域地下淋溶总氮排放量的预测方法
CN104764867B (zh) 一种农田面源污染县域地下淋溶总磷排放量的预测方法
CN104782301A (zh) 一种兼顾区域产量与环境风险的北京地区夏玉米分区施氮的方法
Li et al. Non-film mulching comprehensively improved plant growth and yield of cotton in a deep-drip irrigation system under arid regions
Zhou et al. Historical variability of cotton yield and response to climate and agronomic management in Xinjiang, China
Song et al. Pervasive soil phosphorus losses in terrestrial ecosystems in China

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170510

Termination date: 20180327

CF01 Termination of patent right due to non-payment of annual fee