- 日期:2020-12-30
- 作者:郭财秀 张鹃 李姣 杨开斌 子万玖
- 来源:云南省水文水资源局昆明分局 昆明龙慧工程设计咨询有限公司 云南省水文水资源局红河分局 中国电建集团昆明勘测设计研究院有限公司
本文为原创,可在署名“郭财秀 张鹃 李姣 杨开斌 子万玖”和转自“昆明市水利学会”的前提下转载,否则视为侵权!
隆阳区水库型水源地非点源风险空间分布及评价
隆阳区地处西南边陲,地形起伏大,地貌类型多样,水资源分布不均,人均水资源量较少,属于严重缺水地区。受区域自然条件影响,供水水源类型有水库型、地下水型和河流型,尤以水库型居多。随着隆阳区社会经济的发展,居民生活生产活动带来的水环境恶化、植被退化、水土流失等问题,威胁着饮用水水源安全。饮水安全事关国计民生,保障饮水安全是隆阳区可持续发展的基础保障。
饮用水水源地风险源主要有点源、非点源二类[1,2]。隆阳区水源地多处于山区河谷、沟箐地带,区域内以农业经济为主,生活污染、农业种植、水产养殖等非点源污染是区内水源地水质恶化的主要诱因。目前,水源地调查评价大多是进行水量和水质两方面的安全评价[3-5],对各类风险源尤其是非点源风险源的评价很少。非点源风险源呈零散的面状分布,难以通过传统的污染源调查统计方法进行评价[6]。随着3S技术的发展,RS和GIS技术在区域空间信息的提取、分析及可视化等方面应用越来越广泛,为水源地调查评价提供可视化技术支撑,并逐步运用于水源地水质、污染源空间分布的调查等[7-9]。胡娟[10]基于GIS的空间分析功能实现了汉江、湘江水源地水质的空间可视化分析。殷绪华等[11]将ArcGIS技术运用到污染物入库量的研究。姚延鹃[12]等通过对高空分辨率遥感影像的人工解译,识别水源地保护区内的非点源风险源,揭露了全国25个典型湖库型水源地保护区的非点源风险空间分布。本文以隆阳区内12个水库型水源地为研究对象,基于水源地保护区内的土地利用方式,运用ArcGIS获取水源地保护区内非点源风险的空间分布,并通过建立一套科学的评价方法进行非点源风险定量和定性评价,旨在为隆阳区水源地风险管理提供数据支撑。
1 风险评价方法
1.1 非点源风险源提取
结合实地勘测资料和近期Landsat TM遥感影像,基于ENVI软件,采用人机交互作业方法,通过识别影像光谱、纹理、色调等目视解译得出水源地保护区内的土地利用方式。根据水源地现场走访调查,水源地保护区范围内存在的非点源污染主要由居民生活污水、固体废弃物、农田化肥农药施用、矿山开采和水产养殖等引起的。因此,本次构建的非点源风险源体系由农业用地、居民点、开采用地、水产养殖四类组成,其中农业用地、居民点、开采用地属陆域类风险,水产养殖属水域类风险。基于研究区土地利用方式解译成果,提取农业用地、居民点、开采用地、水产养殖四类土地利用信息,即可统计出非点源风险源的位置和面积。
1.2 非点源风险定量评估
采用非点源风险指数(non-point source pollution risks index, NPPRI)[12]进行水源地非点源风险的定量评估。
式中,NPPRI为非点源风险指数;
wm为m类非点源风险源的权重;
sm为m类非点源风险源的风险值;
非点源风险源的权重由层次分析法进行确定,非点源风险源的风险值根据陆域或水域非点源风险源面积占其所在保护区的对应陆域或水域面积的比例进行赋分。
层次分析法[13]是通过构建判断矩阵,采用专家赋分法确定相互对比的两个因素的相对重要性分值,由此计算出同一层次因子相对上一层因子的相对重要性权重。根据非点源风险源对水源地水质安全的影响程度,结合国内相关研究及水源地相关法律规定,针对本次建立的非点源风险源体系,构建判断矩阵见表1,经计算,矩阵一致性检验参数CR为0.024,小于0.1,具有较好的一致性。通过对判断矩阵(表1)的最大特征根对应的特征向量进行归一化处理得到各非电源风险源的权重见表2。考虑水源地不同级别保护区对水源地保护程度的不同,各级保护区风险源的权重不同,根据一级保护区较二级保护区重要,确定相应权重统计见表2。
表1 非点源风险源判断矩阵
非点源风险源 | 农业用地 | 居民点 | 开采用地 | 水产养殖 |
农业用地 | 1 | 1/2 | 2 | 1/5 |
居民点 | 2 | 1 | 3 | 1/4 |
开采用地 | 1/2 | 1/3 | 1 | 1/6 |
水产养殖 | 5 | 4 | 6 | 1 |
表2 非点源风险源权重
非点源风险源 | 农业用地 | 居民点 | 开采用地 | 水产养殖 |
权重 | 0.1207 | 0.2007 | 0.0744 | 0.6042 |
一级保护区权重 | 0.0805 | 0.1338 | 0.0496 | 0.4028 |
二级保护区权重 | 0.0402 | 0.0669 | 0.0248 | 0.2014 |
1.3 非点源风险定性评估
为了更直观的反映水源地保护区内非点源风险隐患大小,基于非点源风险源对水源地的影响程度及国内相关风险源评价研究[14],采用低、较低、中等、较高、高5个级别对隆阳区水源地保护区内非点源风险进行定性评价。风险的级别越高,水源地遭遇的非点源风险越大。以水源地非点源风险指数为基础进行非点源风险级别划分,计算非点源风险指数的均值(μ)和标准差(σ),按照表3设定范围确定各水源地非点源风险级别[15]。
表3 非点源风险定性分级
风险等级 | 低 | 较低 | 中等 | 较高 | 高 |
设定范围 | <μ-σ | ≥μ-σ 且 <μ-0.5σ | ≥μ-0.5σ 且 <μ+0.5σ | ≥μ+0.5σ 且 <μ+σ | ≥μ+σ |
2 风险评价与分析
2.1 水源地概况
以隆阳区12个水库型水源地为研究对象,分析水源地保护区内非点源风险源的空间分布规律,并进行非点源风险评价。水源地基本信息见表4,水源地保护区面积范围由1.81km2至13.22km2不等,涵盖小(1)型、小(2)型和中型水库,具有较好的代表性。
表4 水源地基本信息表
水源地编号 | 保护区面积/km2 | 保护区中心位置 | |
经度 | 纬度 | ||
1 | 13.22 | 99°20′38″ | 25°9′45″ |
2 | 7.03 | 99°8′50″ | 24°58′39″ |
3 | 5.28 | 99°6′53″ | 25°1′57″ |
4 | 8.01 | 99°10′49″ | 24°57′51″ |
5 | 1.81 | 99°6′1″ | 25°1′15″ |
6 | 7.93 | 99°16′21″ | 25°17′15″ |
7 | 2.21 | 99°17′34″ | 25°9′56″ |
8 | 10.7 | 99°2′24″ | 25°13′42″ |
9 | 14.2 | 99°2′35″ | 25°16′20″ |
10 | 4.19 | 99°0′40″ | 24°55′25″ |
11 | 8.26 | 99°17′28″ | 24°49′15″ |
12 | 3.1 | 98°51′13″ | 25°27′29″ |
2.2 水源地非点源风险空间分布规律
基于遥感影像解译的土地利用数据,运用ArcGIS提取12个水源地非点源风险源的位置和面积,统计得出研究区内非点源风险源信息见表5,非点源风险源分布特点如下:
水源地保护区内均有非点源风险源分布,主要有农业用地、居民点、开采用地和水产养殖4类,水源地非点源风险源总面积为22.9km2,占水源地保护区总面积的26.8%。按非点源风险源类型统计,以农业用地面积占比最大,占非点源风险源总面积的90.0%,其次是居民点、开采用地、水产养殖,分别占9.1%、0.9%、0.003%。按水源地统计,除第5水源地外,其余11个水源地保护区内均有农业用地和居民点的分布,仅第3水源地保护区内有开采用地,仅第12水源地保护区内有水产养殖。按保护区级别统计,除第5水源地一级保护区内无非点源风险源分布外,其余水源地的一级、二级保护区内均有非点源风险源的分布;一级、二级保护区内非点源风险源总面积分别为2.4、20.6km2,分别占一级保护区和二级保护区总面积的16.3%、28.9%。
表5 水源地非点源风险源信息统计表
水源地 编号 | 非点源风险源面积/km2 | ||||||||
一级保护区 | 二级保护区 | 合计 | |||||||
农业用地 | 居民点 | 开采用地 | 水产 养殖 | 农业用地 | 居民点 | 开采用地 | 水产 养殖 | ||
1 | 11.33 | 0.87 | 0.00 | 0.00 | 85.86 | 17.76 | 0.00 | 0.00 | 115.82 |
2 | 21.13 | 1.12 | 0.00 | 0.00 | 235.06 | 26.44 | 0.00 | 0.00 | 283.75 |
3 | 13.28 | 0.00 | 0.00 | 0.00 | 222.70 | 36.25 | 20.79 | 0.00 | 293.02 |
4 | 54.45 | 18.88 | 0.00 | 0.00 | 347.13 | 40.60 | 0.00 | 0.00 | 461.06 |
5 | 0.00 | 0.00 | 0.00 | 0.00 | 10.98 | 0.00 | 0.00 | 0.00 | 10.98 |
6 | 6.98 | 0.26 | 0.00 | 0.00 | 104.46 | 12.56 | 0.00 | 0.00 | 124.26 |
7 | 3.11 | 1.81 | 0.00 | 0.00 | 21.71 | 1.82 | 0.00 | 0.00 | 28.45 |
8 | 31.41 | 0.00 | 0.00 | 0.00 | 193.34 | 14.03 | 0.00 | 0.00 | 238.78 |
9 | 24.52 | 0.00 | 0.00 | 0.00 | 281.17 | 9.51 | 0.00 | 0.00 | 315.20 |
10 | 24.94 | 0.00 | 0.00 | 0.00 | 96.90 | 4.72 | 0.00 | 0.00 | 126.56 |
11 | 18.68 | 0.56 | 0.00 | 0.00 | 111.34 | 19.27 | 0.50 | 0.00 | 150.35 |
12 | 3.08 | 0.57 | 0.00 | 0.06 | 139.68 | 2.00 | 0.00 | 0.00 | 145.39 |
2.3 水源地非点源风险评价结果
2.3.1 定量评价结果
统计12个水源地保护区内各类非点源风险源面积占所在保护区的陆域或水域非点源风险源面积的比例,由公式(1)计算非点源风险指数,成果见表6。绘制非点源风险指数及各级保护区内非点源风险源分布情况图,见图1。
由表6可知,12个水源地中非点源风险指数最小的是第5号水源地,为0.35,最大的是第4号水源地,为9.60。由图1可知,各水源地中均以农业用地面积居首,说明12个水源地中以农业非点源风险隐患较大。由于第5号水源地仅二级保护区内分布有小面积的农业用地,周边环境保护良好,非点源风险指数小;而第4号水源地一级、二级保护区内均有较大面积的农业用地和居民点,非点源风险指数大。
表6 水源地非点源风险值及定量评价结果
水源地 编号 | 非点源风险值 | 非点源风险 指数 | |||||||
一级保护区 | 二级保护区 | ||||||||
农业 用地 | 居民点 | 开采 用地 | 水产 养殖 | 农业 用地 | 居民点 | 开采 用地 | 水产 养殖 | ||
1 | 14.49 | 1.11 | 0.00 | 0.00 | 6.95 | 1.44 | 0.00 | 0.00 | 1.69 |
2 | 35.83 | 1.90 | 0.00 | 0.00 | 36.95 | 4.16 | 0.00 | 0.00 | 4.90 |
3 | 23.69 | 0.00 | 0.00 | 0.00 | 43.31 | 7.05 | 4.04 | 0.00 | 4.22 |
4 | 56.10 | 19.45 | 0.00 | 0.00 | 51.66 | 6.04 | 0.00 | 0.00 | 9.60 |
5 | 0.00 | 0.00 | 0.00 | 0.00 | 8.69 | 0.00 | 0.00 | 0.00 | 0.35 |
6 | 5.22 | 0.19 | 0.00 | 0.00 | 16.32 | 1.96 | 0.00 | 0.00 | 1.23 |
7 | 10.53 | 6.13 | 0.00 | 0.00 | 11.64 | 0.98 | 0.00 | 0.00 | 2.20 |
8 | 23.01 | 0.00 | 0.00 | 0.00 | 23.87 | 1.73 | 0.00 | 0.00 | 2.93 |
9 | 10.39 | 0.00 | 0.00 | 0.00 | 26.48 | 0.90 | 0.00 | 0.00 | 1.96 |
10 | 31.56 | 0.00 | 0.00 | 0.00 | 29.53 | 1.44 | 0.00 | 0.00 | 3.82 |
11 | 17.08 | 0.51 | 0.00 | 0.00 | 15.98 | 2.77 | 0.07 | 0.00 | 2.27 |
12 | 18.78 | 3.48 | 0.00 | 1.08 | 48.55 | 0.70 | 0.00 | 0.00 | 4.41 |
图1 水源地保护区非点源风险源统计与评价
为进一步分析各类因子对非点源风险指数的影响,选取非点源风险指数、一级保护区内非点源风险源面积比例、二级保护区内非点源风险源面积比例、保护区内非点源风险源面积比例、非点源风险源面积、保护区面积数据系列绘制散点图,建立非点源风险指数与各因子的相关关系见图2。
由图2可知,非点源风险指数与各非点源风险源因子的相关系数均在0.76以上,相关关系较好,与水源地保护区面积相关系数仅为0.04,没有相关性,说明水源地遭遇非点源风险隐患的大小与保护区面积大小无关,但与非点源风险的面积及其在各级保护区内占比有关,尤与一级保护区内非点源风险的面积比例相关关系最好。另一方面,由于水源地保护区面积是非点源风险源绝对面积的基数,保护区面积越大,非点源风险的基数越大,治理难度越大。因此进行水源地风险管理时,既要着重一级保护区内风险源的保护与管理,还要考虑整个水源地保护区的污染治理与保护。
a)非点源风险指数与非点源风险源因子相关关系图
b)非点源风险指数与保护区面积相关关系图
图2 水源地保护区非点源风险指数与各因子相关关系图
2.3.2 定性评价结果
12个水源地非点源风险指数的均值为3.30,标准差为2.30,按照表3风险等级的设定范围,确定非点源风险值范围见表7。12个水源地非点源风险定性等级分类结果见图3。由图2与表7可知,12个水源地中,1个水源地(第4水源地)风险定性等级为高,占8.3%;1个水源地(第2水源地)风险定性等级为较高,占8.3%;6个水源地(第3、7、8、10、11、12水源地)风险定性等级为高,占50%;3个水源地(第1、6、9水源地)风险定性等级为较低,占25%;1个水源地(第5水源地)风险定性等级为低,占8.3%。第2、4水源地因一级、二级保护区内均分布有大面积的农业用地和居民点,评价等级分别为较高风险、高风险。风险定性等级较高以上的占16.7%,说明本研究中水源地遭遇非点源风险隐患总体不大。
表7 水源地非点源风险定性等级结果
风险等级 | 低 | 较低 | 中等 | 较高 | 高 |
风险值范围 | <1 | ≥1且<2.15 | ≥2.15且<4.45 | ≥4.45且<5.60 | ≥5.60 |
水源地编号 | 5 | 1、6、9 | 3、7、8、10、11、12 | 2 | 4 |
图3 水源地非点源风险定性等级统计结果
3 结论
(1)本次选取了隆阳区12个水库型水源地,基于遥感影像资料和ENVI软件提取保护区内非点源风险源信息,发现非点源风险源类别中:农业用地分布于12个水源地保护区内,居民点分布于11个水源地保护区内,而开采用地、水产养殖仅分布个别水源地保护区中。说明本研究中水源地保护区内非点源风险主要来自农业面源和居民生活污染,可根据非点源风险的分布位置和面积制定配套的污染治理措施。
(2)根据水源地非点源风险评价结果得出:本研究中的水源地遭遇非点源风险隐患总体不大,但12个水源地保护区内非点源风险现状不同。其中第2、4水源地分别处于较高风险、高风险状态,应优先对其进行保护与管理。
(3)基于不同时期遥感影像即可识别不同时期的非点源风险源信息,以此进行非点源风险源的监测,可为水源地的非点源风险评价提供数据支撑。本文通过对水源地非点源风险的定量与定性评价,直观的反映了保护区内风险状况,明确风险治理类别、位置和数量,并确立了风险管理优先次序,便于制定有效的污染防治措施,对保证区域饮用水水源供水安全具有重要意义。
[2]郑丙辉,付青,刘琰.中国城市饮用水源地环境问题与对策[J].环境保护,2007,(19):59-61.
[3]龙必能.支持向量机在饮用水水源地水质评价中的应用[J].人民珠江,2016,37(5):85-89.
[4]李文攀,朱擎,李东一,刘廷良等.集中式饮用水水源地水质评价方法研究[J].中国环境监测,2015,31(1):24-27.
[5]陈青山.达州市主城区饮用水源评估及保护对策[J].陕西水利,2018,(5):95-96.
[6]管桂玲,卢发周,李萍,张力.长江南京段饮用水水源地风险评估[J].人民珠江,2018,39(08):20-24.
[7]王晓燕,王振刚,王晓峰.GIS支持下密云水库石匣小流域非点源污染[J].城市环境与城市生态,2003,16卷(增刊):26-28.
[8]周慧平,葛小平,许有鹏等.GIS在非点源污染评价中的应用[J].水科学进展,2004,15(4):441-444.
[9]温兆飞,吴胜军,陈吉龙等.三峡库区农田面源污染典型区域制图及其研究现状评价[J].长江流域资源与环境,2014,23(12):1684-1692.
[10]胡娟.石泉县城区水源地水质调查与评价[J].陕西水利,2018,(5):115-116,119.
[11]殷绪华,朱亮,陈琳等.ArcGIS在水源地污染源强度空间分析中的应用[J].河海大学学报(自然科学版),2018,46(5):395-401.
[12]姚延娟,王雪蕾,吴传庆等.饮用水源地非点源风险遥感提取及定量评估[J].环境科学研究,2013,26(12):1349-1355.
[13]郭先华,崔胜辉,赵千钧.城市水源地生态风险评价[J].环境科学研究,2009,22(6):688-694.
[14]刘琰,郑丙辉,万峻,付青.C市饮用水源风险评价实例分析[J].环境科学研究,2009,22(1):52-59.
[15]Colorado Department of Public Health and Envi ronment Water Quality Control Division. Source water assessment methodology for surface water sources and ground water sources under the direct influence of surface water[R].Colorado: EnecoTech Inc.,2004.