丹江口水库水源区农业面源污染物流失量估算

《丹江口水库水源区农业面源污染物流失量估算》由会员分享，可在线阅读，更多相关《丹江口水库水源区农业面源污染物流失量估算（11页珍藏版）》请在装配图网上搜索。

1、丹江口水库水源区农业面源污染物流失量估算王国重;李中原;左其亭;屈建钢;李晓宇【摘 要】丹江口水库是南水北调中线工程的水源地,也是我国规模最大的饮用水水 源保护区,农业面源污染是影响其水质的首要因素,合理估算其水源区农业面源污染 物的流失量对确保该水库水质有着重大意义.在实地调研的基础上,采用分形理论估 算了丹江口水库水源区河南省所在区域2013年TN、TP的流失量,并与输出系数 法的结果进行比较.结果表明:两种方法都显示该区域污染物的流失以氮为主,其中分 形理论结果显示TN流失量是TP流失量的7.156倍（输出系数法为6.773）;禽畜养 殖产生的污染物流失量最多,占总流失量的69.93（输

2、出系数法为70.28）,农田化 肥产生的污染物流失量占总流失量的 21.99（输出系数法为 21.74）;无论是 TN TP的流失总量，还是二者的总流失量,两种算法配对t检验的P值（双侧）均小于0.05, 具有统计学意义;分形方法考虑了下垫面对污染物流失的影响,而输出系数法只是简 单地将各区域的流失率平均化,故分形方法所得结果比输出系数法更为可靠.研究显 示,运用分形方法来研究区域面源污染是适宜的,同时也拓展了流域尺度面源污染的 计算方法.The Danjiangkou Reservoir is the water source of the South-to- North Water Tra

3、nsfer Project,and is also the largest protected zone for drinking water sources in China.Agricultural non-point source pollution is the primary factor that negatively affects its water quality.Therefore,it is of great importance to reasonably estimate the pollution load of agricultural non-point sou

4、rces in the region to ensure the reservoir water quality.Fractal theory was applied to estimate total nitrogen （TN） and total phosphorus （TP） load in the Danjiangkou Reservoir catchment areas in Henan Province in 2013 by onsite investigation.The results from the fractal theory model were compared to

5、 those from the coefficient method.Both methods showed that nitrogen dominates the pollution load in the area,with TN load being 7.156 (6.773 by output coefficient method)times that of TP.The pollution load from livestock farming was the most important contributor,accounting for 69.93 (70.28 by outp

6、ut coefficient method)of the total;the.pollution load from farmland accounted for 21.99 (21.74 by output coefficient method).The P values (double side)from the paired t-test by the two algorithms were less than 0.05 regardless of the total load of TN,TP,or total load of the two,which indicated the a

7、nalysis is statistically significant.The results of the fractal method,which considered the effects of the underlying surface,were better than those of the output coefficient method,which simply averaged the loss rate of the same pollutant from different areas.Thus,it is appropriate to use the fract

8、al method to study regional non-point source pollution,but also extends the calculation method of non-point source pollution in watershed scale.期刊名称】环境科学研究年(卷),期】2017(030)003【总页数】8页(P415-422)【关键词】分形理论;水系分维;农业面源污染;TN、TP流失作 者】 王国重;李中原;左其亭;屈建钢;李晓宇【作者单位】黄河水文水资源科学研究院,河南郑州 450004;河南省水文水资源局, 河南郑州 450003;郑州

9、大学水利与环境学院,河南郑州 450001;河南省水土保持监 督监测总站,河南郑州 450008;黄河水文水资源科学研究院,河南郑州450004【正文语种】中文【中图分类】X524水系是由许多弯弯曲曲的河流交汇而成，如何描述和研究这种弯曲性及其网状结构直困扰着学术界.1977年Mandelbrot 1 首开先河将分形的概念引入了地理水文学，随后他又发表了专著自然界的分形几何，标志着分形理论的诞生 2分形理论为揭示自然界复杂事物的客观规律及其内在联系提供了新的思路 和方法3农业面源污染难以监测，因为其随机性、分散性、隐蔽性、广泛性和不确定性的特 点，般是借助模型来定量分析，如 SWAT、AGNP

10、S、SWMM 等4-6这些 模型多来自国外，需要大量的参数，而且参数的率定要求也很严格，因此其模拟精 度也有限7输出系数法不考虑污染产生与发展的过程，而且参数少、操作便 捷，也有相当的精度，因此也被研究者大量采用8分形理论是人类解决复杂系统的有力工具，其优势在于能够研究传统方法所不能研 究的支离破碎、处处不可微、不光滑，貌似混乱、无规则、随机现象的内部规律， 能够准确、完整地界定实物的主流趋势性质9农业面源污染带有非线性的特 征10，这就为运用分形理论研究农业面源污染提供了可能但是这方面的文 献还十分罕见，则是掌握和运用分形理论的人比较少，再则是能够将分形理论与 面源污染相结合的研究更少河流水

11、系分形是目前分形研究的热点在大量研究的 基础上，国内外学者发现水系分维数与地貌的发育、土壤、土地的利用等因素密切 相关，水系各支流的分维数不仅反映了水系的发育程度，还代表它所处流域的地貌 侵蚀发育的阶段11-14 但许多研究只是通过GIS软件计算水系的分形维数， 而没有作更进一步的研究15-17，如将其与环境污染、流域生态相联系 丹江口水库是南水北调工程的水源地，也是我国规模最大的饮用水水源保护区18 农业面源污染是影响水体水质的首要污染源19,TN、TP又是导致 农业面源污染的主要因素20.该研究以丹江口水库水源区河南省所在的区域 为研究对象，根据ArcGIS软件计算水系分维值，同时结合对该

12、区域农业面源污染 情况的调查结果，运用分形理论来估算该区域TN、TP的污染情况，并与输出系 数法的结果相比较，以证明分形计算结果的可靠性，同时为拓展流域尺度面源污染 的估算方法提供参考. 丹江口水库由丹库和汉库组成，汉库位于湖北省境内，丹库是该研究所在区域，位 于河南省西南部，南接湖北省，西邻陕西省，包含南阳市的西峡县、淅川县，以及 内乡县、邓州市、三门峡卢氏县、洛阳栾川县的部分乡镇，流域总面积为8 047 km2，其中南阳市6 668 km2，三门峡市1 072 km2，洛阳市307 km2 .地理位 置为 11052E112E、3254N34N,海拔为 121-2 212 . 5 m.生物

13、资源 调查结果显示，该水源区为寡污型水体，处于中营养状态18.区域水系主要为丹江及其支流老灌河、淇河和滔河（见图1）.丹江全长378.6 km， 源于陕西凤凰山，经商南县月亮湾入淅川，淅川境内长116 . 6 km,自滔河镇申 明铺以下为丹江口库区;老灌河发源于栾川县小庙岭，经卢氏、西峡，在淅川县注 入丹江，上集镇奕子营以下为库区;淇河发源于卢氏县花园寺西，至淅川县寺湾乡 老君洞附近汇入丹江，全长150 km;滔河发源于陕西商南县白鲁础乡白龙洞，流 经商南县、郧县，在淅川县滔河乡注入丹江，河流全长155 . 2 km .上述水系信 息由中国科学院计算机网络信息中心提供的DEM数据获得，采用Ar

14、cGIS软件提 取水系并按照Strahler分级方法，将上述水系划分成11个子流域根据2014年 河南省统计年鉴，结合走访调查，确定各水系的人口、耕地面积、畜禽养殖情况(见表1)21 计盒维数的计算分形维数是分形理论的主要参数，有多种定义和计算方法，如相似维数(Ds)、豪斯 道夫维数(DH)、计盒维数(Db)、信息维数(Di)、关联维数(Dg)、容量维数(DC) 等其中，计盒维数(Db)因计算简便而最为常用，反映的是分形体对空间的占据 程度，用边长为r的小方盒子去覆盖分形体，非空小方盒的数量记为N(r)，不断缩 小r的取值，相应地得到一系列的N(r)值，当r-0时，则计盒维数21 定义为 实际

15、工作中，一般借助地理信息系统(GIS)软件利用网格法来计算计盒维数22. GIS技术可以获得研究区的总体数据，比传统的人工计算或计算机编程 计算更为客观、准确、经济、高效常用的GIS软件主要有ArcGIS、MapGIS、 GeoStar等，这里采用ArcGIS 10.0,其水系分维计算的具体步骤:先由 ArcGIS具箱中的conversion tools工具，设置需要输出的网格尺寸，同时将矢量 河网数据转换成栅格文件;打开相应栅格文件的属性表，可以看到各子流域所占 的网格数量;由Excel软件将网格尺寸和各子流域相应的网格数量点绘在双对数 坐标图上，相应直线的斜率即为各子流域水系的分维数. 该

16、研究将栅(网)格尺寸从大到小依次取 5 000、3 000、1 000、500、300、100、 50、10、5、1 mm,则各子流域所占栅格数如表2所示根据上述步骤即可得到 各子流域的分维数. 水系分维数乘以不同污染物的平均流失率，即为考虑下垫面作用的面源污染物的流 失率23. 污染物的平均流失率是在实地监测和调查的基础上计算而来,虽然不能尽善尽美, 但有一定的科学依据，也是目前所能达到的水平肥料的平均流失率是以我国农业 种植区划和优势农产品区划为依据,在主要农作物种植区选取典型种植制度和有代 表性地貌的农田,按照不施肥和常规施肥两种处理分别设置监测小区进行对照试验, 通过对其地表径流和地下

17、淋溶进行为期一年的实地监测，分别计算两种处理下的肥 料流失量，二者之差占肥料使用量的比例即为肥料的平均流失率24畜禽粪 便的平均流失率是指畜禽粪便在堆放、冲洗过程中流失到水体的比率，地形、降水 条件的差异会导致各区域的流失率存在差别农村生活污水是指农村居民在日常生 活或为日常生活提供服务的活动中产生的生活污水，生活污水和畜禽粪便的平均流 失率是在对流域调查、监测分析的基础上获得2522 输出系数法 输出系数法是20世纪70年代初北美国家在评价土地利用和湖泊富营养之间的关 系时提出的，该方法由于所需参数少、操作简便，在大中尺度的流域上也能保证一 定的精度该方法认为不同污染物具有不同的输出系数输出

18、系数乘以污染物的产 生量即为流域的面源污染负荷，其表达式8 一般为式中:Lj为流域内污染物j的 总负荷，t (hm2a);i为流域内土地、禽畜或人口分类，共有n类;Cij为污染物j在 第i类土地、牲畜或人口中的输出系数；Ai为第i类土地的面积(hm2)或牲畜(头)、 人口(人)的数目H为降水输入的污染物数量，t (hm2a),该研究不考虑该项的影 响该式中相关参数的取值根据研究区实际情况同时结合已有研究结果 8确定3 1各水系面源污染物TN、TP的流失率根据对研究区的走访调查以及已有成果26-30 ,确定农田化肥中TN、TP平均 流失率为12 5%和3%，村民粪尿和生活污水中TN、TP的平均流

19、失率为5%和 50%，禽畜粪尿中TN、TP的平均流失率如表3所示.各水系所在的地域和环境条件不同，导致污染物的流失率必然有差异分形作用下 农田化肥、禽畜养殖、农村生活污染物的流失率如表4、5所示.水系分维反映了水系发育的复杂程度，分维数越大说明水系发育的越好31 由表4、5可知，由于引入了分形维数，各子流域的TN、TP流失率随着 分维值的不同而存在差异，子流域3的TN、TP流失率稍大，因其分维值稍大， 相应地子流域9的TN、TP流失率则较低.3 2面源污染物TN、TP的流失量污染物的流失率乘以其产生量即为污染物的流失量据2014年河南省统计年鉴和 相关文献32-35，同时结合河南省农村用水定额

20、和污水排放标准，确定研究区 2013年流域农田化肥中TN、TP的施用量，禽畜养殖中TN、TP的排泄量以及生 活污染物中TN、TP的产生量(见表6)，相应污染物的流失量见表7 由表7可见，研究区2013年共流失污染物6 864.46 t,其中TN流失量为6 022 . 86 t,TP流失量为841.610 t,前者是后者的7 . 156倍，说明研究区产 生的污染物以氮的流失为主;禽畜流失的污染物量为4 800 . 120 t,占研究区污染 物总量的69 93% ,这表明应当加强对禽畜粪尿的管理;农田中化肥产生的污染物 量为1 509 . 612 t,占研究区污染物总量的21 . 99%，说明农田

21、中的养分流失也 是农业面源污染的重要来源.从水系角度来看,子流域11产生的污染物最多,对研究区污染的贡献也最大,达 26 . 33%；其次是子流域5;污染物流失最少的是子流域3，对环境污染的贡献率仅 为1.66%.相应输出系数法计算的污染物流失量如表8所示.3.3两种方法结果的比较该研究中的两种算法都是针对同一个区域进行的,其结果存在着必然的联系,可采 用配对样本均数的t检验进行分析比较这里以两种方法计算的TN流失总量、TP 流失总量和污染物总流失量为例进行分析,结果如表9所示.由表9可知，TN和TP流失总量、污染物总流失量检验的t值分别为2 . 774、2 . 869、2.790,均大于t0

22、 . 05(10)=2 . 228，其双侧检验的P值分别为 0.020、0.017、0.019,均小于0.05,表明检验结果具有统计学意义.分 形方法考虑了下垫面对污染物流失的影响,同一污染物在不同区域的流失率不同; 而输出系数法则是将污染物的流失平均化,认为各区域的流失率相同,由此可知, 分形算法的结果更符合实际，优于输出系数方法a）研究区污染物的流失主要表现为氮的流失，TN流失量是TP流失量的7156（输出系数法为6773）倍;禽畜养殖产生的污染物占污染物总量的6993%（输出系数法为7028%）农田化肥产生的污染物流失量占流失总量的2199%（输出系数法为2174%）表明应当加强对禽畜粪

23、尿的管理，采取措施 减少农田中养分的流失b）分形维数会对污染物的流失强度产生影响流失强度小并不代表分形维数就小，如子流域4的流失强度最小但其分维数却不小，子流域9的分维数最小但其流失 强度却最高，这表明农业面源污染受流域地形地貌、土壤、植被、人类活动等因素 的综合影响，分形方法能够反映这种影响，因为水系分维与之息息相关.c）分别对TN和TP流失总量、污染物总流失量进行配对t检验分析，均显示分形 方法的结果优于输出系数法因为输出系数法没有考虑下垫面因素，将之平均化， 认为其流失率相同，故分形法所得结果更符合实际、更为可靠.【相关文献】1 MANDELBROT B B . Fractals:for

24、m , chanee , and dimension M. San Francisco:W H Freeman and Company , 1977:5-30 .2 MANDELBROT B B.The fractal geometry of natureM.New York:W H Freeman and Company，1982:3-20.3 KURAKIN A.The self-organizing fractal theory as a universal discovery method:the phenomenon of lifeJ.Theoretical Biology and

25、Medical Modelling，2011，8（4）:1-66.4 LUO Y，ZHANG M.Management-oriented sensitivity analysis for pesticide transport in watershed-scale water quality modeling using SWATJ.Environmental Pollution，2009，157:3370-3378.5 JAEPIL C，SEUNGWOO P，SANGJUN I.Evaluation of agricultural nonpointsource（AGNPS）model for

26、 small watersheds in Korea applying irregular cell delineationJ.Agricultural Water Management，2008，95:400-408.6 ZAGHLOUL N A，KIEFA M A A.Neural network solution of inverse parameters used in the sensitivity-calibration analyses of the SWMM model simulations JAdvances in Engineering Software，2001，32(

27、7):587-5957 姚瑞华，王东，赵越，等非点源污染负荷模型的研究进展J.水文，2012 , 32(2):11- 15YAO Ruihua ， WANG Dong ， ZHAO Yue ， et alResearch progress of non-point source pollution models in water environment JJournal of China Hydrology ， 2 0 1 2 ， 32(2):11-158 马亚丽，敖天其，张洪波，等基于输出系数模型濑溪河流域泸县段面源分析J.四川 农业大学学报，2013 , 31(1):53-59 . MA

28、Yali , AO Tia nqi , ZHANG Ho ngbo , et al. No n- point source analysis of Laixi River basin in Luxian County based on output coefficient modelJ.Journal of Sichuan agricultural university，2013，31(1):53-59.9 PEITGEN H O，JURGENS H，SAUPE D.Chaos and fractals:new frontiers of science M.2nd ed.Berlin:Spri

29、nger Press，2004: 32-49.10 PARR A，ZOU S，MCENROE B.Effectsofinfiltrationon agricultural runoff contaminationJ.Journal of Environmental Engineering ， 1998 ， 124(9):863-868 .11 ROBERT A，ROY A G.On the fractal interpretation of the mainstream lengthdrainage area relationships J . Water Resources Research

30、 ， 1990 ， 26(5):839-842 .12 VENEZIANO D ， NIEMANN J D . Self-similarity and multifractality of fluvial erosion topography:2.scaling properties J . Water Resources Research，2000，36(7):1937- 1951.13 TSAKIROGLOU C D ， FLEURY M.Pore network analysis of resistivity index for water-wet porous media J . Tr

31、ansport in Porous Media，1999，35:89-128.14 TURCOTTE D L.Self-organized complexity in geomorphology: observations and models J . Geomorphology ， 2007 ， 91(3 4): 302-310 .15 LIU Dong，LUO Mingjie，FU Qiang，et al.Precipitation complexity measurement using multifractal spectra empirical mode decomposition

32、detrended fluctuation analysisJ.Water Resources Management，2016，30(2):505-522.16 KIM J C，JUNG K.Fractal tree analysis of drainage patterns J . Water Resources Management，2015，29(4):1217-1230.17 DOMBRADI E，TIMAR G，BADA G，et al.Fractal dimension estimations of drainage network in the Carpathian-Pannon

33、ian system J.Global and Planetary Change，2007，58(1 2 3 4):197-213.18 胡兰群，冯精兰，李怡帆，等.南水北调中线工程水源地丹江口水库生物监测试点研究J.河南师范大学学报(自然科学版),2014,42(3):100-104 . HU Lanqun , FENG Jinglan , LI Yifan，et al.Pilot study on biomonitoring in Danjiangkou Reservoir of the water source area in the middle route of Chinas Sou

34、th to North Water Diversion ProjectJ.Journal of Henan Normal University(Natural Science)，2014，42(3):100-104.19 SHEN Zhenyao，QIU Jiali，HONG Qian，et al.Simulation of spatial and temporal distributions of non-point source pollution load in the Three Gorges Reservoir region JScience of the Total Environ

35、ment，2014，493:138-14620 GUO H Y，WANG X R，ZHU J GQuantification and index of nonpoint source pollution in Taihu Lake region with GISJEnvironmental Geochemistry and Healthy ， 2004，26:147-15621 SKUBALSKA-RAFAJLOWICZ EA new method of estimation of the box-counting dimension of multivariate objects using

36、 spacefilling curves JNonlinear Analysis， 2005，63:1281-128722 GE Meiling，LIN QizhongRealizing the box-counting method for calculating fractal dimension of urban form based on remote sensing image JGeo-Spatial Information Science，2009，12(4):265-27023 姚素珍基于分形理论的清河水库流域非点源污染研究D 阜新:辽宁工程技术大学， 2012:40-4324

37、 国务院第一次全国污染源普查领导小组办公室第一次全国污染源普查:农业污染源肥料流 失系数手册R.北京:国务院第一次全国污染源普查领导小组办公室，2009:8-120 .25 尹微琴，王小治，王爱礼，等太湖流域农村生活污水污染物排放系数研究:以昆山为例J.农业环境科学学报，2010 , 29 (7):1369-1373 .YIN Weiqi n , WANG Xiaozhi , WANGAili，et alDischarge index of pollutants from village sewage in Taihu region:a case studyin KunshanJ.Journa

38、l of Ago-Environment Science，2010，29(7): 1369-1373.26 夏颖，汪荣勇，高立，等.湖北省油菜原棉花轮作系统地表径流氮磷流失特征 J .湖北 农业科学， 2014 ， 53(23):5751-5754 . XIA Ying，WANG Rongyong，GAO Li，et al.Runoff losses characteristics of nitrogen and phosphorus from rapeseed-cotton rotation systems in HubeiProvinceJ.HubeiAgricultural Scienc

39、es ， 2014 ， 53(23):5751-5754 .27 李晓光，周其文，胡梅，等.中国畜禽粪便污染现状及防治对策 J .中国农学通报，2008，24(S1):77-80.LI Xiaoguang，ZHOU Qiwen，HU Mei，et al.Environmental contamination situation from animal faeces in China and control countermeasuresJ.Chinese Agricultural Science Bulletin，2008，24(S1):77-80.28 李荣刚，夏源陵，吴安之，等.江苏太湖地

40、区水污染物及其向水体的排放量 J .湖泊科 学，2000 , 12(2) :147-153 . LI Ron gga ng , XIA Yua nling , WU An zhi , et al. Polluta nts sources and their discharging amount in Taihu Lake area of Jiangsu ProvinceJ.Journal of Lake Sciences，2000，12(2):147-153.29 全国规模化畜禽养殖污染调查办公室.全国规模化畜禽养殖业污染情况调查技术报告 R 北京:全国规模化畜禽养殖污染调查办公室，2001:

41、4-55.30 孙瑞敏.湖北省农村生活污水水量水质调查与分析 D .武汉:武汉理工大学， 2010:4-6 .31 NIKORA V.Fractal structures of river plan formsJ.Water Resources Research ， 1991，27(6):1327-1333.32 刘钦普.河南省化肥施用强度地理分布及其环境风险评价 J .河南农业科学， 2014 ，43(6):66-70LIU QinpuRegional distribution and environmental risk assessment of chemical fertilizer

42、application rate in Henan Province JJournal of Henan Agricultural Sciences，2014，43(6):66-7033 国务院第一次全国污染源普查办公室第一次全国污染源普查畜禽养殖业源产排污系数手 册R.北京:国务院第一次全国污染源普查办公室，2009:3-53 .34 赵东美于桥水库库区农业非点源污染及其控制对策研究D 天津:天津大学， 2007:39-4135 JONSSON R ， VINNERAS BAdapting the nutrient content of urine and feces in different countries using FAO and Swedish data in eco-sanclosingtheloop C WEMER C，AVENDANO V ， DEMSAT S，et alEcosan-closing the loopProceedings of the 2nd InternationalSymposium on EcologicalSanitationLuebeck: Germany，2003:623-626