ArcGIS在水源地污染源强度空间分析中的应用

Vol.46No.5Sep.2018

DOI:10.3876/j.issn.10001980.2018.05.004

ArcGIS在水源地污染源强度空间分析中的应用

殷绪华1,2,朱摇亮1,2,陈摇琳1,2,吴新明3,方长风4,朱月华5,邱云鹏1,2,谭相君1,2

2.河海大学环境学院,江苏南京摇210098;3.江苏省骆运水利工程管理处,江苏宿迁摇223800;4.溧水区水务局,江苏南京摇211200;5.溧水区水务局方便水库管理所,江苏南京摇211212)

(1.河海大学浅水湖泊综合治理与资源开发教育部重点实验室,江苏南京摇210098;

摘要:以南京市溧水区方便水库为研究对象,运用ArcGIS对其汇水区域进行子流域划分,同时结合遥感技术对各子流域污染源强度和空间分布进行模拟,确定方便水库各入库河流主要污染源及汇水区域内污染源强度分布情况,然后利用实际监测值对ArcGIS模拟值进行验证。结果表明:ArcGIS技术可直观反映各入库河流汇水面积大小及各污染源强度的空间分布情况;ArcGIS估算值与实测值计算的污染物入库量的误差范围在10%以内。排放量;方便水库

中图分类号:X524摇摇摇文献标志码:A摇摇摇文章编号:10001980(2018)05039507

关键词:ArcGIS;水文分析;子流域划分;土地利用;汇水面积;污染源强度;污染物空间分布;污染

ApplicationofArcGISinthespatialdistributionanalysisof

pollutionsourceintensityinwatersource

YINXuhua1,2,ZHULiang1,2,CHENLin1,2,WUXinming3,FANGChangfeng4,

ZHUYuehua5,QIUYunpeng1,2,TANXiangjun1,2

(1.KeyLaboratoryofIntegratedRegulationandResourcesDevelopmentofShallowLakes,

HohaiUniversity,Nanjing210098,China;

2.CollegeofEnvironment,HohaiUniversity,Nanjing210098,China;

3.JiangsuLuoyunWaterProjectManagementDepartment,Suqian223800,China;

4.LishuiDistrictWaterBureau,Nanjing211200,China;

5.LishuiFangbianReservoirManagementOffice,Nanjing211212,China)

Abstract:TheFangbianReservoirintheLishuidistrictofNanjingCitywasselectedasthestudyobject.TheArcGIStechnologywasusedinthewatershedsubdivisionofcatchmentareainFangbianReservoir,andthencombinedwiththeremotesensingtechnology,thisstudysimulatedthepollutionintensityandspatialdistributiontodeterminethemainpollutionsourceofinflowstreamandtheintensitydistributionofpollutionsource.Afterthat,thesimulatingresultsofArcGIShasbeenvalidatedbytheactuallymonitoredwaterquality.TheresultsshowedthatthecatchmentareaofinflowstreamandspatialdistributionofpollutionsourceintensitycanbedirectlyreflectedbytheArcGIStechnology.Besides,thedifferenceofpollutionemissionamountcalculatedbetweentheArcGISsimulationandthesurveyeddatawaslessthan10%.

Keywords:ArcGIS;hydrologicalanalysis;watersheddivision;landuse;catchmentarea;pollutionsourceintensity;spatialdistributionofthepollution;pollutionemissionamount;FangbianReservoir

摇摇作者简介:殷绪华(1992—),男,硕士研究生,主要从事水污染控制理论与技术研究。E鄄mail:84551735@qq.com

通信作者:朱亮,教授。E鄄mail:zhulianghhu@163.com

引用本文:殷绪华,朱亮,陈琳,等.ArcGIS在水源地污染源强度空间分析中的应用[J].河海大学学报(自然科学版),2018,46(5):395鄄401.

YINXuhua,ZHULiang,CHENLin,etal.ApplicationofArcGISinthespatialdistributionofpollutionsourceintensityinwatersource[J].JournalofHohaiUniversity(NaturalSciences),2018,46(5):395鄄401.

396

河海大学学报(自然科学版)第46卷

表明,我国湖库型水源地的主要污染物包括化学需氧量(COD)、总氮(TN)和总磷(TP)等,污染来源主要是水库[4]。传统的污染源调查方法主要采用现场采样、区域入户调查和水质监测分析等,无论是数据的时效性还是准确度都存在一定的局限,而且消耗大量的人力物力[5鄄7]。

泛[8鄄11]。Khalil等[12]运用数据库和GIS技术评估废弃矿场的环境影响及预测污染形成过程,开发出基于GIS

ArcGIS技术由于其强大的空间制图、空间可视化、空间分析和数据管理功能,在各类规划中应用广

我国各类型地表水源地中,湖库型水源地的比例约为40%,是我国饮用水源的重要类型之一[1]。研究

点源污染、面源污染和内源污染[2鄄3]。其中,点源污染和面源污染主要是通过入库河流来水携带污染物进入

的环境数据库,为其他废弃矿场的污染研究提供依据。周慧平等[13]运用GIS技术强大的可视化功能识别非点源污染的主要来源、迁移途径,预报污染产生负荷、变化趋势及对水体水质的影响。王晓燕等[14]在密云水库石匣小流域非点源污染研究中,运用GIS技术确定流域内非点源污染的分布规律和重点控制区。温兆飞等[15]利用GIS的空间叠加分析法,对三峡库区农田面源污染的关键影响因子进行筛选,获取农田面源污染典型区域,并发现不同典型区分布差异的主导因素是土壤类型。

目前,运用ArcGIS技术模拟污染源分布情况已有了相应的研究,但用其进行污染物入库量的模拟计算却很少。笔者运用ArcGIS技术结合调查数据和遥感技术对方便水库土地利用方式进行分析,同时对入库污染源源强空间分布进行模拟,并用模拟结果计算各污染物的入库量,然后通过实际监测数据进行验证。将ArcGIS技术应用到污染物入库量的研究中,旨在为河流的污染源源强分布研究提供有效的分析手段。

1摇研究区概况及监测布点

1.1摇研究区概况

方便水库保护区位于南京市溧水区城区东北7km处,库区地处东经119毅05忆07义~119毅07忆02义、北纬31毅40忆~31毅43忆之间。方便水库在整个汇水区域范围内已进行了保护区的划定:一级保护区范围为校核洪水位(29郾15m)以下的库区,即生态红线范围之内的区域,面积约为12郾21km2;二级保护区范围为生态红线以外所有的汇水区域(生态黄线范围),面积1.2摇监测布点约为郾92km2。

77郾13km2,其中主要入库河流为青龙桥河、四古凹河、西南村河、西村河、杨家坝河、谢家棚河、王家山河、后吕河这河流。为了验证ArcGIS模拟结果的可行性和准确性,在方便水库各入库河流上分别设置采样点(图1),根据实际监测的数据计算各污染物的入库量。监测指标有COD、TN、TP和NH3鄄N,水质指标分析按照《水和废水监测分析方法》[16]。

研究范围为整个方便水库汇水区域,总面积为

2摇研究方法

2.1摇流量计算

图1摇方便水库汇水区域及各入库河流监测点布局Fig.1摇Layoutofmonitoringpointsaroundthe

FangbianReservoiranditsinflowrivers

根据方便水库1975—2012年降水资料,采用水文手册

中的方法,对陈墩库区入库水量进行计算。频率分析计算表明,1978年为枯水年(保证率P=90%),1992年为平水年(P=50%),2002年为丰水年(P=20%)。根据上述3个代表年全年各月的降雨分配情况,查《江苏省水文手册》秦淮河山丘区月降雨-径流关系曲线图,推算得到陈墩库区枯水年、平水年、丰水年入库流量分别为3郾46万m3、32郾83万m3、91郾46万m3。2014年总降雨量为1368郾9mm,为丰水年(P=20%),年入库流量按91郾46万m3计算。根据ArcGIS水文分析计算出的各入库河流所占比例,推算得到各条河流的年入库流量。

第5期殷绪华,等摇ArcGIS在水源地污染源强度空间分析中的应用

397

2.2摇子流域划分

条入库河流的汇水面积[17鄄18]。分析过程为:DEM数据生成、DEM栅格数据填洼处理、流向和洼地识别处理、计算汇流累积量、计算河网栅格数据、计算河流连接数据、计算子流域栅格、转化矢量数据后合并[19]。分析2.3摇土地利用类型分析

过程结合流域的实际情况,选择适当的汇流累积量阈值,将整个小流域划分为8个子流域。

通过ArcGIS对数字高程模型(DEM)数据进行水文分析,得到方便水库汇水区域内河流水系分布及各

土地利用数据解译是基于谷歌卫星影像遥感信息,结合实地调查和其他辅助数据,采用全数字人机交互作业方法,主要根据对图像光谱、纹理、色调等的认识结合地形图目视解译而成。同时,根据溧水方便水库汇水区域水田为主、坑塘众多的特点,采用由点到面、由易到难的方式进行数据解译。按照全国土地利用分类体系,将土地利用类型分为耕地、林地、草地、水域、城乡工矿居民用地、未利用土地等6个一级类型和24个2.4摇ArcGIS模拟污染排放量分布

区域内污染物排放量主要是通过现场调查以及相关资料收集统计。结合土地利用图,将各种污染排放量按土地利用类型进行分类及赋值,区域重合的部分数值进行叠加,以此模拟得到整个汇水区域内各入库河流小流域的污染排放量分布情况。然后,按污染排放量属性对其进行分级,分级按照ArcGIS软件中提供的自然间断点分级法,此方法可自动识别、分组数据,使各级之间差异最大化。最后,运用ArcGIS的矢量剪切2.5摇污染排放量分布计算结果验证

功能计算各入库河流汇水面积内的污染排放量,各入库河流小流域污染排放量及贡献率。

由于各个污染源产生点距离库区的距离不能确定,污染物进入库区的途径差别较大,汇入库区过程中自然降解以及削减情况未能有准确、详细的调查数据,因此污染物入库量通过污染物排放量乘以相关入库系数计算得出,最后用实际监测数据计算污染物入库量,作为对计算结果的验证。计算过程中,不考虑入库沿途物理、化学等方面的自然降解作用,从而核算出水源保护区水库最大污染物入库量,以最不利情况考虑水源保护区的水质污染问题。

a.根据ArcGIS模拟估算污染排放量及排放系数计算污染入库量。各类污染物入库量计算公式如下:

W入库量=(W排放量-兹)茁

(1)

二级类型。

式中:W排放量———不同来源的污染物年排放量;兹———污水处理厂处理掉的污染物量,取零;茁———不同来源生活污染物的入库系数。

根据《全国水环境容量核定技术指南》及其他大中型水库的实际调研结果,水库各类污染源入库系数见表1。

b.根据河流入库处的实际监测数据及河

表1摇各类污染源入库系数

Table1摇Storagecoefficientofdifferentpollutionsources

污染源工业、城镇生活农村生活农田径流禽畜养殖水土流失

COD0郾8~1郾00郾1~0郾20郾1~0郾30郾2~0郾40郾1~0郾3

TN0郾8~1郾00郾1~0郾20郾1~0郾30郾2~0郾40郾1~0郾3

TP0郾8~1郾00郾1~0郾20郾1~0郾30郾2~0郾40郾1~0郾4

0郾8~1郾00郾1~0郾20郾1~0郾30郾2~0郾40郾1~0郾5NH3鄄N

流入库流量计算污染物入库量。计算公式如下:

W入库量

(表2)。

CQ琢

=106

(2)

式中:C———入库河流入库处污染物实测质量浓度;Q———年入库总流量;琢———各入库河流的入库流量占比

表2摇入库河道汇水情况

Table2摇Catchmentconditionofinflowrivers

河道名四古凹河西南村河西村河杨家坝河

汇水面积/km2

1郾1776郾5922郾7811郾614

琢/%1郾814郾282郾49

谢家棚河青龙桥河王家山河后吕河河道名

汇水面积/km2

11郾62535郾6173郾6171郾883

17郾9154郾875郾572郾91琢/%

10郾16

398

河海大学学报(自然科学版)第46卷

3摇结果与讨论

3.1摇河流水系分布及汇水面积

通过ArcGIS对水文DEM数据的模拟分析,得到方便水库汇水区域内河流水系分布(图2)及各条入库河流的汇水面积青龙桥河自溧水区共和山区汇流而来,来水量占总入库水量的流时间短,坡面水流很快汇入河道,并挟带大量泥沙,是主要的入库水沙来源。除青龙桥河外,谢家棚河及西南村河相较于其他河3.2摇土地利用类型分析

流来水量较大,分别占总来水量的17郾91%和10郾16%。

方便水库水源地保护区总面积约77郾13km2,其中水域面积

图2摇方便水库汇水区河网分布Fig.2摇RivernetworkdistributionofcatchmentareaintheFangbianreservoir

(表2)。由表2可以看出青龙桥河是汇入方便水库的主要河流,50%以上。该河系山区性河流,降雨后径流汇聚快,来水较猛,汇

10郾25km2,土地面积66郾88km2,各用地类型面积、所占土地面积比例以及在各级保护区内的分布见表3。从表3可知,方便水库水源地一级保护区内大部分为耕地,耕地比例高达72郾4%,其次为林地,以及城乡、工矿、居民用地,所占比例均为11郾76%。二

级保护区内占比最高的仍为耕地(54郾19%),其次为林地(32郾78%)。综合整个汇水区域来看,水源保护区内耕地和城乡、工矿、居民用地总的比例过高,而林地和草地所占比例偏低。

表3摇方便水库汇水区域土地利用状况

Table3摇LandusestatusinthecatchmentareaofFangbianReservoir

土地利用类型

耕摇地林摇地草摇地

城乡、工矿、居民用地

未利用土地

一级保护区

面积/km2

4郾250郾690郾240郾690

占比/%72郾4011郾754郾0911郾760

二级保护区

面积/km233郾0620郾000郾746郾940郾27

占比/%54郾1932郾781郾2111郾380郾44

整个汇水区

面积/km237郾3120郾690郾987郾630郾27

占比/%55郾7930郾941郾4711郾410郾39

3.3摇饮用水源地污染源强度计算与验证3.3.1摇水源地污染排放量

经现场调查以及溧水区水务局提供的相关资料统计,方便水库主要污染源包括工业企业和集约化畜禽

养殖造成的点源污染和农村生活污水、农田径流、农村固体废弃物、水土流失以及分散式禽畜养殖造成的面源污染。考虑到方便水库水面以下底泥面积和底泥释放通量等监测困难,故在研究方便水库污染源的过程中,只考虑外源污染。方便水库保护区污染物产生量统计见表4。综合点源和面源2种污染来源,方便水库的主要污98郾71%。其中,农田径流污染是面源污染中的第一大污染源,这与种植面积及施肥量密切相关。

表4摇方便水库保护区不同来源污染物排放量及贡献率

TN

TP

COD

染源为面源污染,面源污染的COD、TN、TP、NH3鄄N排放量占各自总污染排放量的97郾60%、99郾67%、98郾10%、

Table4摇EmissionamountandcontributionrateofpollutantsfromdifferentsourcesintheFangbianreservoir

污染类别工业污染源集约化禽畜养殖

生活污水农田径流农村固体废弃物

水土流失分散畜禽养殖

排放量/(t·a-1)1郾8136郾16222郾831830郾12165郾203郾84

贡献率/%0郾081郾609郾8680郾987郾310郾17

排放量/(t·a-1)0郾182郾7362郾37694郾3659郾7360郾790郾79

贡献率/%0郾020郾317郾0878郾826郾786郾900郾09

排放量/(t·a-1)0郾053郾036郾67134郾687郾0110郾400郾41

贡献率/%0郾031郾874郾1183郾014郾326郾410郾25

排放量/(t·a-1)0郾145郾8343郾62385郾1427郾470郾32

NH3鄄N

贡献率/%0郾031郾269郾4383郾275郾940郾07

摇摇注:工业污染源、集约化禽畜养殖为点源污染,其他为面源污染。

第5期殷绪华,等摇ArcGIS在水源地污染源强度空间分析中的应用

399

3郾3郾2摇ArcGIS模拟计算汇水区域内污染物排放量分布

结合RS技术对区域内土地利用方式的调查及各类型土地利用方式污染排放系数,同时结合表4,模拟

计算整个汇水区域内各入库河流小流域的污染物排放量分布情况(图3)。从图3可以直观地看出不同入库河流小流域污染物排放量存在明显差异。图3中红色区域分别表示COD排放量为175郾0~325郾4t/a、TN排放量为61郾2~121郾2t/a、TP排放量为12郾6~23郾4t/a和NH3鄄N排放量为35郾4~65郾8t/a,均分布在青龙桥河66郾2t/a、TP排放量为6郾5~12郾6t/a和NH3鄄N排放量为18郾3~34郾5t/a,主要分布在青龙桥河流域、西村河可见这2个流域污染负荷年排放量较小。

流域和谢家棚河流域内;其次为橙红色区域,分别表示COD排放量为90郾6~175郾0t/a、TN排放量为33郾7~

流域、西南村河流域以及四古凹河流域小部分区域内。后吕河流域和王家山河流域内多为深绿色和浅绿色,

图3摇方便水库保护区污染物年排放量分布

Fig.3摇AnnualdischargedistributionofpollutionintheFangbianReservoir

运用ArcGIS的矢量剪切功能计算各入库河流汇水面积内的污染物排放量,各入库河流小流域污染物排放量及贡献率见表5。由表5可知,青龙桥河流域及谢家棚河流域污染负荷产生量较大,主要原因:(a)二者的流域汇水面积较大,青龙桥河、谢家棚河流域汇水面积各占水库总汇水面积的46郾20%和15郾08%;(b)处在这2条河流汇水区域内的农田及人口较多,且苗木种植也多在这2个流域范围内,故农田径流污染及农村生活、农村固废污染等面源污染的排放量大。

表5摇方便水库保护区入库河流小流域污染排放量及贡献率

Table5摇Pollutionemissionamountandcontributionrateofsub鄄watershedintheFangbianReservoir

编号12345678

流域青龙桥河流域谢家棚河流域西南村河流域杨家坝河流域王家山河流域四古凹河流域后吕河流域西村河流域

排放量/(t·a-1)1330郾77393郾21224郾97113郾2070郾5866郾4132郾1628郾67

COD

贡献率/%58郾8817郾409郾955郾013郾122郾941郾421郾28

排放量/(t·a-1)

508郾41159郾08郾1444郾9227郾3126郾1514郾2311郾73

TN

贡献率/%57郾7118郾065郾103郾102郾971郾621郾3210郾12

排放量/(t·a-1)

95郾7728郾498郾114郾944郾702郾372郾1016郾18

TP

贡献率/%58郾8817郾529郾954郾993郾042郾1郾451郾28

排放量/(t·a-1)

271郾2881郾86郾6223郾5214郾6913郾856郾5郾81NH3鄄N

贡献率/%58郾4017郾6210郾045郾063郾162郾981郾481郾26

3.3.3摇污染排放量分布计算结果的验证

按照不同污染物来源及相应的入库系数,计算出不同河流小流域污染入库量及贡献率。结果显示按照排放量及入库系数计算出来的COD、TN、TP、NH3鄄N年入库量分别为:903郾99t、88郾10t、14郾60t、41郾60t。青龙桥河流域对COD、TN、TP、NH3鄄N贡献率最大,分别达到了58郾88%、57郾71%、58郾88%、58郾40%;其次为谢家棚河,其4种污染物的贡献率分别为17郾40%、18郾06%、17郾52%、17郾62%。同时,按照实际监测数据及入库水量计算出来的COD、TN、TP、NH3鄄N年入库量分别为:1000郾33t、95郾63t、14郾15t、45郾41t。此种计算方法下其次为谢家棚河,其4种污染物的贡献率分别为18郾23%、16郾49%、13郾55%、20郾13%。

仍然是青龙桥河流域对COD、TN、TP、NH3鄄N贡献率最大,分别达到了57郾40%、51郾74%、54郾35%、50郾21%;

从各条入库河流的贡献率结果比较中可以看出(图4),与实际监测数据计算结果相比,用污染排放总量计算出的入库河流污染负荷入库量与其较为吻合,只有贡献率较小的几条河流结果稍有不同,但是相差不大,用污染排放量总量计算出的COD、TN、TP、NH3鄄N总入库量分别相差10%、8%、3%和8%。

400

河海大学学报(自然科学版)第46卷

图4摇不同算法下污染入库量比较

Fig.4摇Comparisonofpollutionstorageunderdifferentalgorithms

4摇结摇摇论

库各类污染源中,面源污染所占比例最大,其对COD、TN、TP、NH3鄄N污染负荷的贡献率分别达97郾60%、家棚河流域,其次为西村河流域和西南村河流域。位于水库南边的青龙桥河流域对COD、TN、TP、NH3鄄N贡献率最大,分别达到了58郾88%、57郾71%、58郾88%、58郾40%。污染物入库量的估算值和实测值的误差范围在10%以内,进而说明用ArcGIS模拟计算污染源强空间分布具有一定的可行性和准确度。参考文献:

[1]冉飞.湖库型饮用水水源地水环境生态管理和保护[J].山西建筑,2011,37(34):184鄄186.(RANFei.Thewater

environmentecosystemmanagementandprotectiontothedrinkingwatersourceareasoflakesandreservoirstype[J].ShanxiArchitecture,2011,37(34):184鄄186.(inChinese))

a.在方便水库主要入库河流中,青龙桥河的汇水面积最大,其次为谢家棚河和西南村河。方便水

99郾67%、98郾10%、98郾71%。

b.ArcGIS模拟的污染排放量分布结果显示:污染负荷排放量最大的区域主要集中在青龙桥河流域和谢

[2]边博,朱伟,李冰,等.太湖流域西部地区面源污染特征及其控制技术[J].水资源保护,2015,31(1):48鄄55.(BIANBo,

ZHUWei,LIBing,etal.Characteristicsandcontroltechniquesofnon鄄pointpollutioninthewesternregionofTaihuBasin[J].WaterResourcesProtection,2015,31(1):48鄄55.(inChinese))

[3]郑丙辉,付青,刘琰.中国城市饮用水源地环境问题与对策[J].环境保护,2007(19):59鄄61.(ZHENGBinghui,FU

Qing,LIUYan.EnvironmentalproblemsandcountermeasuresofurbandrinkingwatersourcesinChina[J].EnvironmentalProtection,2007(19):59鄄61.(inChinese))

[4]唐访良,张明,徐建芬,等.千岛湖库区及其主要入库河流水中有机氯农药残留污染特征及健康风险评价[J].环境科

学,2014,35(5):1735鄄1741.(TANGFangliang,ZHANGMing,XUJianfen,etal.Pollutioncharacteristicsandhealthriskassessmentoforganochlorinepesticides(ocps)inthewaterofLakeQiandaoanditsmajorinputrivers[J].EnvironmentalScience,2014,35(5):1735鄄1741.(inChinese))

[5]李晓燕,李壁成,安韶山.基于GIS的小流域坡耕地调查方法探讨[J].水土保持通报,2001,21(5):37鄄40.(LIXiaoyan,

LIBicheng,ANShaoshan.SlopingfarmlandinvestigationmethodinsmallwatershedbasedonGIS[J].BulletinofSoilandWaterConsercation,2001,21(5):37鄄40.(inChinese))

[6]周超,高乃云,赵世嘏,等.上海青草沙水库水质调查与评价[J].同济大学学报(自然科学版),2012,40(6):4鄄9.

(ZHOUChao,GAONaiyun,ZHAOShijia,etal.EvaluationofwaterqualityinQingcaoshaReservoirofShanghai[J].Journal

第5期殷绪华,等摇ArcGIS在水源地污染源强度空间分析中的应用

401

[7]卫丹,万梅,刘锐,等.嘉兴市规模化养猪场沼液水质调查研究[J].环境科学,2014,35(7):2650鄄2657.(WEIDan,

WANMei,LIURui,etal.Studyonthequalityofdigestedpiggerywastewaterinlarge鄄scalefarmsinJiaxing[J].EnvironmentalScience,2014,35(7):2650鄄2657.(inChinese))

[8]邓必平,严恩萍,洪奕丰,等.基于GIS和DEM的东江湖流域水文特征分析[J].湖北农业科学,2013,52(15):3531鄄

3536.(DENGBiping,YANEnping,HONGYifeng,etal.HydrologicalcharacteristicanalysisofDongjiangLakeWatershedBasedonGISandDEM[J].HubeiAgriculturalScience,2013,52(15):3531鄄3536.(inChinese))Chinese))

Chenghu,SHAOQuanqin.OceanGISdatamodelandstructure[J].JournalofGeo鄄InformationScience,2003,5(4):25鄄29.(in[9]仉天宇,周成虎,邵全琴.海洋GIS数据模型与结构[J].地球信息科学,2003,5(4):25鄄29.(ZHANGTianyu,ZHOU[10]王云鹏,闵育顺,傅家谟,等.水体污染的遥感方法及在珠江广州河段水污染监测中的应用[J].遥感学报,2001,5(6):

460鄄465.(WANGYunpeng,MINYushun,FUJiamo,etal.RemotesensingmethodofwaterpollutionandapplicationonwaterpollutionmonitoringinGuangzhouSectionofPearlRiver[J].JournalofRemoteSensing,2001,5(6):460鄄465.(inChinese))hydrologyandwaterresourcesanddevelopmenttrend[J].JilinWaterResources,2009(6):1鄄5.(inChinese))ofanabandonedminesiteinMorocco[J].JournalofGeochemicalExploration,2014,144:468鄄477.

(WANGGuangming,LIANGXiujuan,XIAOChanglai,etal.ThestatusoftheapplicationofGIStechnologyinthefieldof[11]王光明,梁秀娟,肖长来,等.GIS技术在水文水资源领域中的应用现状与发展趋势[J].吉林水利,2009(6):1鄄5.

ofTongjiUniversity(NaturalScience),2012,40(6):4鄄9.(inChinese))

[12]KHALILA,HANICHL,HAKKOUR,etal.GIS鄄basedenvironmentaldatabaseforassessingtheminepollution:acasestudy[13]周慧平,葛小平,许有鹏,等.GIS在非点源污染评价中的应用[J].水科学进展,2004,15(4):441鄄444.(ZHOUHuiping,[14]王晓燕,王振刚,王晓峰.GIS支持下密云水库石匣小流域非点源污染[J].城市环境与城市生态,2003(增刊1):29鄄31.

SmallWatershed,MiyunReservoir[J].UrbanEnvironment&UrbanEcology,2003(Sup1):29鄄31.(inChinese))2004,15(4):441鄄444.(inChinese))

GEXiaoping,XUYoupeng,etal.ApplicationofGIStononpointsourcepollutionevaluation[J].AdvancesinWaterScience,(WANGXiaoyan,WANGZhengang,WANGXiaofeng.AGIS鄄basedstudyonlossofnon鄄pointsourcepollutionsfromShixia

[15]温兆飞,吴胜军,陈吉龙,等.三峡库区农田面源污染典型区域制图及其研究现状评价[J].长江流域资源与环境,2014,

23(12):1684鄄1692.(WENZhaofei,WUShengjun,CHENJilong,etal.Typicalregionmappingofagriculturalnon鄄pointTheYangtzeBasin,2014,23(12):1684鄄1692.(inChinese))

sourcepollutionintheThreeGorgesReservoirArea:patternanalysisandstatusevaluation[J].ResourcesandEnvironmentin

[16]国家环境保护总局《水和废水监测分析方法》编委会.水和废水监测分析方法[M].4版.北京:中国环境科学出版社,[17]李晶,张征,朱建刚,等.基于DEM的太湖流域水文特征提取[J].环境科学与管理,2009,34(5):138鄄142.(LIJing,

ZHANGZheng,ZHUJiangang,etal.HydrologicalfeatureextractionofTaihuLakeBasinbasedonDEM[J].EnvironmentalScienceandManagement,2009,34(5):138鄄142.(inChinese))

[18]顾用红,舒振文,张火青.DEM在流域水文特征分析中的应用[J].人民珠江,2001,22(4):5鄄6.(GUYonghong,SHU[19]OMRANA,DIETRICHS,ABOUELMAGDA,etal.NewArcGIStoolsdevelopedforstreamnetworkextractionandbasin

delineationsusingPythonandjavascript[J].Computers&Geosciences,2016,94:140鄄149.

(收稿日期:20170901摇编辑:刘晓艳)

2002.

2001,22(4):5鄄6.(inChinese))

Zhenwen,ZHANGHuoqing.ApplicationofDEMmethodtoanalysisofhydrologicalfeaturesinthewaterbasin[J].PearlRiver,