裴文明，张 慧,2，鞠昌华，姚素平，任永乐，马孟枭，仇宽彪

煤炭资源大量开发产生地表塌陷、地质环境破坏等问题，在一定程度上阻碍矿井安全生产或地区发展。近年来，采煤沉陷区评价及治理技术研究受到政府职能部门、科研院所等相关机构重视，成为矿山绿色开采的研究热点。为推广该领域最新科技成果，促进学术交流，我刊设立专题，选登4篇论文集中报道相关科研成果。

基于韦伯–费希纳定律的淮南采煤沉陷水域水环境综合预警评价

裴文明1，张 慧1,2，鞠昌华1，姚素平3，任永乐3，马孟枭1，仇宽彪1

(1. 生态环境部南京环境科学研究所，江苏 南京 210042；2. 南京信息工程大学 江苏省大气环境与装备技术协同创新中心，江苏 南京 210044；3. 南京大学 地球科学与工程学院，江苏 南京 210023)

在我国东部高潜水位地区，长期的地下煤炭开采造成地表形成大面积的沉陷水域，受周围人类活动的影响，沉陷水域水环境污染风险逐渐增大。以淮南潘集采煤沉陷水域为研究对象，根据韦伯–费希纳定律，选取与沉陷水域水环境关系密切的水质评价指标，分析不同时间段沉陷水域的水环境状况，通过图层叠加融合定量反映沉陷水域水环境质量时空变化及潜在水环境风险的发生时间和区域。研究结果显示，沉陷水域水质不容乐观，存在水质恶化的风险。从5月到11月，随着时间的变化，沉陷水域都呈现出中警以上的警情，在不同月份沉陷水域出现的警情区域也有所不同，其中沉陷水域的西部和中部警情出现的概率较大。在今后的水环境管理中需针对不同时间不同区域采取针对性的治理措施，避免水环境污染风险的发生。

韦伯–费希纳定律；沉陷水域；水质评价；综合预警；淮南潘集矿区

我国东部平原矿区有其特殊的地理位置、地质环境和水文条件，长期的大规模地下煤炭开采活动造成地表形成了大面积的沉陷水域。据统计，全国采煤沉陷区面积已超过2万km2[1]。随着地下采煤活动的持续进行，沉陷水域范围将持续扩大，对区域地表水系结构和生态系统产生的影响也将越来越大[2-7]。

近年来，社会经济水平不断提高及生态文明建设不断深入，沉陷水域的关注度不断提高，关于沉陷水域的水质监测、水环境治理等研究也不断增多。吴建宇[8]和张维翔[9]利用常规水质监测手段，对比分析了沉陷水域、河流等水体水质指标的差异，研究不同类型水体水环境质量变化的影响因素；叶圆圆[10]和陈晓晴[11]利用遥感和GIS手段，通过建立水质指标模型，反演沉陷水域水质状况，有效弥补了常规水样采集的缺失。然而，目前的研究多针对沉陷水域单次采样的水质指标进行水环境方面的评价研究，尚缺乏对沉陷水域不同时间段，不同季节条件下系统性的水环境监测，监测指标与数据积累性不够，沉陷水域水环境管理所需数据不足，导致沉陷水域水环境安全受到威胁时，不能及时准确掌握沉陷水域水环境状况及其发展趋势。

为此，本文以淮南潘集沉陷水域为例，根据心理学理论韦伯–费希纳定律，利用多期水质采样数据，选取与沉陷水域水环境关系较为密切的水质参数指标，分析沉陷水域水环境的变化过程，评价沉陷水域潜在的水环境污染风险，以期为沉陷水域的水环境安全管理和水质改善提供科学依据和技术支持。

1 研究区概况与样品采集

本研究选择淮南市潘集区潘一矿采煤沉陷水域为研究对象。潘一矿于1984年投产，截至2019年底，原煤产量累计已超过一亿吨。沉陷水域形成于20世纪80年代，主要由农田、村庄沉陷形成，面积约3.50 km2，平均水深4.5 m，是淮南市因地下采煤活动较早形成的沉陷水域之一，现已用于水产养殖，是当地一种特殊的地表水资源[12]。沉陷水域为半封闭水域，与泥河相连，最终汇入淮河。由于沉陷水域周边人类活动频繁，工业和生活污水排放，农业活动产生的污染物等通过地表径流进入水体，对沉陷水域水环境造成严重的影响。

本研究共在沉陷水域中进行了5次水体样品采集，时间分别为2016年3月11日、5月4日、7月17日、9月23日和11月23日，采样点布设如图1所示。水样采集、样品保存和处理均按照国家相关标准和规定完成，样品采集后及时送往实验室进行水质指标的分析测试。

图1 淮南潘集矿区地理位置及采样点位分布

2 评价指标与方法

2.1 预警评价方法

韦伯–费希纳(Weber-Fechner，记作W-F)定律是定量描述人类感觉强度与外界环境刺激强度关系的心理物理学公式[13]。近些年，由于其特别的数学思想和基本含义，逐渐被一些学者应用到环境评价预警领域[14-17]。本研究将沉陷水域水质指标的含量作为外界环境刺激强度，水环境质量作为人体反应量，对沉陷水域水环境进行定量评价预警。

W-F定律基本表达式为：

=lg(1)

式中：为人体反应量；为外界刺激量；为韦伯常数。

W-F定律应用于水环境评价中，基于以下3点假设：①将外界刺激量视为水环境中某项指标的浓度大小；②将人体反应量视为该污染物对人体的影响程度[18]；③韦伯常数视为环境评价中各污染物的权重[19]。

基于以上假设，W-F定律函数关系式表示为：

k=alg(c+1) (2)

式中：k为第个采样点第个水质参数对人体的影响程度；a为第个采样点第个水质参数的权重；c为第个采样点第个水质参数的标准化值；c的目的是为了让lg(c+1)＞0，经数学推导证明此处理方法不影响最终评价结果[20]。

综合影响指数k计算公式为：

利用式(3)计算得出5次采样时间沉陷水域水质的综合影响指数(k)，并依次计算5次采样每个采样点的综合影响指数。考虑到各项指标的高低区别，对水体环境影响不同，需要对其进行权重计算。某项指标的变异系数较大，就可以认为此指标样点差异显著，可以作为环境评价的主要依据[17]。本文采用变异系数法计算5次采样水体指标的权重值。

2.2 预警评价指标体系

评价体系分为目标层、准则层和指标层3层。其中，目标层反映本次评价体系建立的主要目标任务；准则层是将一个大目标分解为多个目标或准则，进而分解成多个指标层；指标层是详细描述各准则层的独立因子，这些因子是构建整个评价体系的核心内容。本文从两个方面选取具有综合性、代表性和可操作性的评价指标。一方面是对地表水水质影响较大的6种常规水质参数指标(溶解氧DO、悬浮物SS、总有机碳TOC、叶绿素a、总磷TP、总氮TN)，另一方面是水体中对人体健康影响较大且采煤活动极易释放的8种重金属元素指标(铁Fe、锰Mn、锌Zn、铜Cu、铬Cr、镉Cd、铅Pb、镍Ni)，共计14个指标组成预警评价指标体系(图2)。

一般情况下，预警级别可分为无警、轻警、中警、重警和巨警5个等级[20-21]。各水质指标相应预警级别的分级标准主要采用GB 3838—2002《地表水环境质量标准》和文献[22]，各等级标准见表1。

图2 预警评价指标体系框架

表1 预警评价指标体系分级标准

3 结果与讨论

3.1 权重分析

由表2可知，常规水质参数5次采样准则层的权重值分别为0.40、0.41、0.28、0.23、0.09。其中，3月份常规指标中叶绿素a对水体影响最大，权重值达到0.49，其次是悬浮物和总氮，权重值分别为0.17和0.11。5月份常规指标中对水体影响最大的指标是悬浮物和总磷，权重值分别为0.34和0.21。7月份和9月份常规指标中悬浮物、叶绿素a、总磷和总氮的权重值相差不大，基本在0.2左右。11月份常规指标的悬浮物、叶绿素a、总氮权重值基本在0.23左右，说明此月份对水体影响较大的指标为悬浮物、叶绿素和总氮。总体上看，5次不同采样时间6种水质常规参数指标中，叶绿素a和悬浮物的平均权重值分别为0.25和0.24，说明在沉陷水域，浮游植物和水中悬浮物的含量对水质影响较大，是影响水环境最主要的两项指标，与之前的研究结果基本一致[12]。

从5组样品的8项重金属元素指标权重值分析，5次采样的平均权重值分别为0.60、0.59、0.72、0.77、0.91，说明重金属是影响沉陷水域水质警级的主要因素。3月份，锌、镍和镉3种重金属元素对水体警级的划分影响较大，权重值分别为0.24、0.18和0.15。5月份，镉、锌和铅3种重金属元素对水体警级的划分影响较大，权重值分别为0.28、0.19和0.19。7月份，铅和铜重金属元素对水体警级的划分影响较大，权重值为0.17。9月份，镉元素对水体警级的划分影响较大，权重值达到0.41。11月份，铜、铅和锌3种重金属元素对水体警级的划分影响较大，权重值分别为0.26、0.25和0.23。从单要素来看，水体中重金属镉、锌、铅对沉陷水域水质警级影响较大，与王兴明[23]研究成果基本一致。其中，镉元素5次采样平均权重值为0.21，9月份权重值达到0.41。其次为锌元素和铅元素，5次采样平均权重值分别为0.21和0.17。最小为铬和铁元素，5次采样平均权重值分别为0.056和0.066。

表2 水环境预警评价指标的权重值

3.2 预警结果

为综合反映沉陷水域的水质状况，本文利用韦伯–费希纳定律，将6种常规水质参数指标和8种重金属元素指标按照权重占比进行综合计算。通过采用GIS的栅格计算器，获得5次采样各项指标的空间分布图。叠加融合各指标图层后得到不同采样时间的沉陷水域水质警级状态(图3)，并统计5次采样各警级水域的占比情况(表3)。

从图3可以看出，沉陷水域水环境状况存在明显的时空变化特征。时间上，3月份水体整体呈现无警状态。5月份，沉陷水域中警及以上占比达到45%，其中巨警范围所占比例居全年最高，约5%，主要集中在沉陷水域西部，泥河入水口处。这可能与农业生产有关，上游河流两岸农田开始播种，农业化肥使用增加，污染物流入泥河后汇入沉陷水域。据水质监测数据显示，5月总磷和总氮指标值均达到或超过Ⅴ类标准；加上春季入梅之前，降雨较少，对污染物的排泄能力较差。7月份，警情有所减缓，沉陷水域中警及以上占比达到30%左右，这可能是因为春末夏初，降雨增多，水体流动性增强，对污染物存在一定的稀释作用。重警、巨警水域主要位于水体东西部交汇处和水域东部出水口处。这与特殊的区域位置存在一定关系，中部和东部出水口处通道较为狭窄，水流较缓，水中营养物质容易在此处富集。9月份，警情稍微有所增加，沉陷水域中警及以上占比超过30%，主要位于沉陷水域西部。这可能与夏秋季的暴雨和台风有关，位于沉陷水域西部的煤矸石堆受雨水冲刷淋溶出来的重金属污染元素流入沉陷水域[24]。11月份，水体状态主要以无警为主，占比超过60%，重警、巨警水域仅占8%，主要集中在沉陷水域中部。

总体上看，沉陷水域存在水质恶化的风险，这与已有研究结果基本一致，淮南矿区沉陷水域水质已呈现富营养状态[25]。沉陷水域周围人口密集，人类活动频繁，农业生产过程中使用的化肥、农药，含有磷、氮化合物，在地表和土壤中累积后，随地表径流进入沉陷水域。同时，周边的工矿企业以及居民生活排放的各种污染物、营养物质进到沉陷水域内，给沉陷水域带来不同程度的影响。时间上分析，从5月份开始，随着温度升高，降雨较少，水体中的藻类、浮游物质繁殖旺盛；再者是当地农业生产播种施肥的季节，造成水体藻类、浮游植物过度繁殖的磷、氮元素会被排入沉陷水域内，加重水体的污染程度，沉陷水域出现预警。到11月份，随着气温下降，沉陷水域内藻类、浮游植物等大量死亡，水体的污染情况减弱，沉陷水域预警情况逐渐消失。空间上分析，由于沉陷水域有泥河沟通，水体流动性较强，在一定程度上对沉陷水域水质造成影响，沉陷水域不同时间段水环境呈现出明显的区域变化。沉陷水域的西部和中部出现重警、巨警的预警级别比例最高，成为容易出现水质恶化的水域，也是沉陷水域水环境管理中需要关注的区域。在今后的水环境管理中，可以根据沉陷水域水环境预警时空变化规律，提出有针对性的管理措施，并进行因地制宜的水环境治理。

图3 不同月份水环境预警评价各警级分布和占比

表3 不同月份不同预警级别占比

4 结论

a. 基于韦伯–费希纳定律对沉陷水域进行水环境预警，通过选取6种常规水质参数和8种重金属元素指标，并根据不同采样时间赋予不同的权重值，可定量反映出沉陷水域水环境质量变化及潜在的水环境风险发生的时间和区域，未来能够更好地与水环境管理目标相结合，广泛应用于沉陷水域的水环境质量评价。

b.受降水、气温、水文条件的影响，沉陷水域水环境预警结果显示出明显的时空变化特征。3月，沉陷水域基本处于无警状态，5月到11月，沉陷水域都呈现出中警以上的警情，其中，5月、7月、9月出现的警情较高，水域面积较大。同时，在不同月份沉陷水域出现警情的水域有所变化，沉陷水域的西部和中部因特殊的区位条件，造成警情出现的比重较高，在今后的水环境管理中需要在不同时间和不同区域制定针对性的预防措施，避免水环境风险的发生。

请听作者语音介绍创新技术成果等信息，欢迎与作者进行交流

[1] 胡炳南，郭文砚. 我国采煤沉陷区现状、综合治理模式及治理建议[J]. 煤矿开采，2018，23(2)：1–4. HU Bingnan，GUO Wenyan. Mining subsidence area status，syntheses governance model and governance recommendation[J]. Coal Mining Technology，2018，23(2)：1–4.

[2] 段洪涛，罗菊花，曹志刚，等. 流域水环境遥感研究进展与思考[J]. 地理科学进展，2019，38(8)：1182–1195. DUAN Hongtao，LUO Juhua，CAO Zhigang，et al. Progress in remote sensing of aquatic environments at the watershed scale[J]. Progress in Geography，2019，38(8)：1182–1195.

[3] 陈珏. 我国水环境监测存在的问题及对策[J]. 环境与发展，2019，31(4)：181–184. CHEN Jue. Problems and countermeasures of water environment monitoring in China[J]. Environment & Development，2019，31(4)：181–184.

[4] 李茜，张鹏，彭福利，等. 国家水环境质量预报预警研究进展及业务发展思路[J]. 中国环境监测，2019，35(1)：8–16. LI Qian，ZHANG Peng，PENG Fuli，et al. Research progress and preliminary plan of national water quality forecasting and alarming system[J]. Environmental Monitoring in China，2019，35(1)：8–16.

[5] 李佳洺，余建辉，张文忠. 中国采煤沉陷区空间格局与治理模式[J]. 自然资源学报，2019，34(4)：867–880. LI Jiaming，YU Jianhui，ZHANG Wenzhong. Spatial distribution and governance of coal-mine subsidence in China[J]. Journal of Natural Resources，2019，34(4)：867–880.

[6] 潘仁飞. 煤矿开采生态环境综合评价及生态补偿费研究[D]. 北京：中国矿业大学，2010. PAN Renfei. Research on the synthetical assessment to eco-environment and ecological compensation of mining[D]. Beijing：China University of Mining and Technology，2010.

[7] 陈晨，李兵，徐燕飞. 平原高潜水位采煤沉陷区规划与综合利用策略：以淮南矿区为例[J]. 安徽农学通报，2019，25(16)：125–126. CHEN Chen，LI Bing，XU Yanfei. Discussion on planning and comprehensive utilization strategy of coal mining subsidence area in high phreatic level in plain：Taking Huainan mining area as an example[J]. Anhui Agicultural Science Bulletin，2019，25(16)：125–126.

[8] 吴建宇. 封闭式采煤沉陷积水区水环境特征及水质评价研究[D]. 淮南：安徽理工大学，2018. WU Jianyu. Study on water environmental characteristics and water quality evaluation of closed coal mining subsidence area[D]. Huainan：Anhui University of Science & Technology，2018.

[9] 张维翔. 淮南高潜水位采煤沉陷区水质特征及变化趋势[D]. 合肥：安徽大学，2019. ZHANG Weixiang. Water quality characteristics and changing trend of coal mining subsidence area with high underground water in Huainan[D]. Hefei：Anhui University，2019.

[10] 叶圆圆. 基于RS淮南采煤沉陷水域水质实时监测技术研究[D].淮南：安徽理工大学，2014. YE Yuanyuan. Research on real-time monitoring of water quality based on remote sensing technology[D]. Huainan：Anhui University of Science & Technology，2014.

[11] 陈晓晴. 基于GIS潘谢矿区氮磷流失特征研究：以潘集、谢桥塌陷区为例[D]. 淮南：安徽理工大学，2013. CHEN Xiaoqing. Study on the losses characteristic of nitrogen and phosphorus based on GIS in Panxie mine[D]. Huainan：Anhui University of Science & Technology，2013.

[12] PEI Wenming，YAO Suping，DONG Shaochun，et al. Using field spectral measurements to estimate chlorophyll-a in waterlogged areas of Huainan，China[J]. GIScience & Remote Sensing，2015，52(6)：660–679.

[13] DRÖSLER J. An N-dimensional Weber law and the corresponding Fechner law[J]. Journal of Mathematical Psychology，2000，44(2)：330–335.

[14] 巩如英，王飞，刘雅莉，等. 韦伯–费希纳定律评价模型在景观环境质量评价中的应用[J]. 西北林学院学报，2006，21(1)：131–135. GONG Ruying，WANG Fei，LIU Yali，et al. Assessment of landscape environmental quality based on Weber-Fechner’s law[J]. Journal of Northwest Forestry University，2006，21(1)：131–135.

[15] 张宝，刘静玲，陈秋颖，等. 基于韦伯–费希纳定律的海河流域水库水环境预警评价[J]. 环境科学学报，2010，30(2)：268–274. ZHANG Bao，LIU Jingling，CHEN Qiuying，et al. Research of water environment forewarning for reservoirs in Haihe River Basin[J]. Acta Scientiae Cirumstantiae，2010，30(2)：268–274.

[16] 李小燕，王菲凤，张江山. 基于韦伯–费希纳定律的湖泊富营养化评价[J]. 水电能源科学，2011，29(3)：37–39. LI Xiaoyan，WANG Feifeng，ZHANG Jiangshan. Lake eutrophication assessment based on Weber-Fechner law[J]. Water Resources and Power，2011，29(3)：37–39.

[17] 林秀珠，饶清华，欧晓敏. 韦伯–费希纳模型在湖泊(水库)预警评价中的应用[J]. 水土保持通报，2018，38(2)：285–291.LIN Xiuzhu，RAO Qinghua，OU Xiaomin. Application of Weber-Fechner model to water environment forewarning of lakes(reservoirs)[J]. Bulletin of Soil and Water Conservation，2018，38(2)：285–291.

[18] 钟龙芳，王菲凤，张江山. 基于韦伯–费希纳定律的地下水环境质量评价[J]. 环境科学与管理，2012，37(12)：189–192. ZHONG Longfang，WANG Feifeng，ZHANG Jiangshan. Groundwater quality assessment based on Weber-Fechner law[J]. Environmental Science and Management，2012，37(12)：189–192.

[19] 李因果，李新春. 综合评价模型权重确定方法研究[J]. 辽东学院学报(社会科学版)，2007，9(2)：92–97. LI Yinguo LI Xinchun. Weight determination of comprehensive evaluation model[J]. Journal of Eastern Liaodong University(Social Sciences)，2007，9(2)：92–97.

[20] 李祚泳，彭荔红. 基于韦伯–费希纳拓广定律的环境空气质量标准[J]. 中国环境监测，2003，19(4)：17–18. LI Zuoyong，PENG Lihong. Environmental air quality standard based on Weber-Fischna’s law[J]. Environmental Monitoring in China，2003，19(4)：17–18.

[21] 周爱仙. 煤矿区生态环境现状评价及预警研究以南屯煤矿区为例[D]. 济南：山东师范大学，2006. ZHOU Aixian. The actuality evaluation and prewarning study of the environment of coal mining area[D]. Jinan：Shandong Normal University，2006.

[22] 李子成，邓义祥，郑丙辉. 中国湖库营养状态现状调查分析[J].环境科学与技术，2012，35(61)：209–213. LI Zicheng，DENG Yixiang，ZHENG Binghui. Investigation of the eutrophication status of the lakes and reservoirs in China[J]. Environmental Science & Technology，2012，35(61)：209–213.

[23] 王兴明. 淮南煤矸石堆积地重金属元素环境生物地球化学研究[D]. 淮南：安徽理工大学，2013. WANG Xingming. Environmental biogeochemistry of heavy metals in different mediums in dumping sites of coal mine spoil in Huainan[D]. Huainan：Anhui University of Science & Technology，2013.

[24] 任永乐，董少春，姚素平. 淮南塌陷塘重金属空间分布特征研究[J]. 煤田地质与勘探，2018，46(1)：125–134.REN Yongle，DONG Shaochun，YAO Suping. Spatial distribution characteristics of heavy metals in Huainan subsidence pond[J]. Coal Geology & Exploration，2018，46(1)：125–134.

[25] 裴文明，张慧，姚素平，等. 淮南矿区不同类型沉陷水域水质遥感反演和时空变化分析[J]. 煤田地质与勘探，2018，46(3)：85–90. PEI Wenming，ZHANG Hui，YAO Suping，et al. Remote sensing inversion and analysis on spatial-temporal variation of water quality in different types of subsided waterlogged zones in Huainan mining area[J]. Coal Geology & Exploration，2018，46(3)：85–90.

Water environment comprehensive forewarning for waterlogged area in Huainan based on Weber-Fechner law

PEI Wenming1, ZHANG Hui1,2, JU Changhua1, YAO Suping3, REN Yongle3, MA Mengxiao1, QIU Kuanbiao1

(1. Nanjing Institute of Environmental Sciences, Ministry of Ecology and Environment, Nanjing 210042, China; 2. Jiangsu Collaborative Innovation Center of Atmospheric Environment and Equipment Technology, Nanjing University of Information Science & Technology, Nanjing 210044, China; 3. School of Earth Sciences and Engineering, Nanjing University, Nanjing 210023, China)

In area with high groundwater level in the east of China, long-time underground coal mining results in a large subsidence water area on the surface. The risk of water environment pollution in subsidence waters is increasing gradually with the influence of human activities.According to Weber-Fisher law, this paper took Panji waterlogged area of Huainan as the research object, selected the water quality evaluation indexes closely related to the water environment of the waterlogged area, analyzed the water environmental conditions of the waterlogged area in different periods, and quantitatively reflected the time and the area of the water environment quality change and potential water environment risk of the waterlogged area through the overlay and fusion of layers. The results show that the water quality of waterlogged area is not optimistic, and there is a risk of water quality deterioration. The waterlogged area presents warning conditions above middle level with the change of time from May to November. The warning areas changes in different months. Among them, the western and central part of the waterlogged area has a higher probability of occurrence of warning situation. In the future water environment management, it is needed to take targeted preventive measures for different areas during different periods to avoid the risk of water environment pollution.

Weber-Fisher law; waterlogged area; water quality assessment; comprehensive forewarning; Panji coalfield of Huainan

P641.69；X141

A

10.3969/j.issn.1001-1986.2020.03.001

1001-1986(2020)03-0001-07

2019-12-03；

2020-04-29

国家重点研发计划课题(2017YFC0506604)

National Key R&D Program of China(2017YFC0506604)

裴文明，1988年生，男，江苏南京人，博士，助理研究员，从事生态环境监测、环境规划等研究工作. E-mail：438661015@qq.com

张慧，1968年生，女，河南开封人，博士，研究员，从事生态安全、环境规划等研究工作. E-mail：zhanghui@nies.org

裴文明，张慧，鞠昌华，等. 基于韦伯–费希纳定律的淮南采煤沉陷水域水环境综合预警评价[J]. 煤田地质与勘探，2020，48(3)：1–7.

PEI Wenming，ZHANG Hui，JU Changhua，et al. Water environment comprehensive forewarning for waterlogged area in Huainan based on Weber-Fechner law[J]. Coal Geology & Exploration，2020，48(3)：1–7.

(责任编辑 周建军)