目前,国内外常用的水质评价法主要包括单因子评价法、综合污染指数法、主成分分析法、内梅罗污染指数法、(改进)模糊评价法、人工神经网络评价法等[1-6]。内梅罗指数法因计算简单,能够对污染环境做出定性分析和定量评价,可综合反映出污染的性质和程度等优点而被广泛应用。1974年,美国叙拉古大学内梅罗教授在《河流污染科学分析》中第1次提出内梅罗指数的概念,并采用该法对河流水环境质量进行评价[7]。1980年,中国学者金士博等[8]首次采用内梅罗指数法对漓江水质进行评价。随后在地下水[9-11]、地表水[12-14]、近岸海域[15]、湿地水环境[16-17]等水质评价,以及沉积物[18-19]、土壤[20-22]等环境要素中污染物的环境健康风险评估中,内梅罗指数法被广泛研究应用,也是我国GB/T 14848—1993《地下水质量标准》中规定的评价方法。
传统的内梅罗指数法在实际应用过程中因存在只考虑算术平均值和最大值,未考虑污染因子权重等问题,会导致评价结果不能真实反映环境要素质量[23],因此,国内外学者采用不同修正方法对内梅罗指数进行修正,以便更加可靠应用于环境要素评价。为了探究内梅罗指数法在北戴河国家湿地公园水质评价中是否合理可靠,本文以2018年9月北戴河国家湿地公园22个水样点监测数据为基础,对比多种修正内梅罗指数法的评价结果,并与综合水质标识指数和营养状态指数进行相关性分析。
传统的内梅罗指数是一种包含单因子指数的算术平均值和最大值Fmax,并强调最大值污染因子对环境质量的影响程度。数学表达式为:
P传
(1)
(2)
(3)
式中:Ci为第i种污染因子的实测浓度;Sij为第i种污染因子在第j种标准下的标准值;P传为传统内梅罗指数。
实际水质评价中,存在某污染因子实测浓度较低,但所占权重较大,对水质的评价影响较大。因此,在传统方法上考虑各污染因子权重,对均值进行修正[12]。
(4)
(5)
(6)
(7)
式中:Wi为第i种污染因子的权重值; Ri为第i种评价指标的相关性比值;Smax为第j类标准下污染因子中的最大标准值。
传统的内梅罗指数法忽略了权重最大污染因子对环境质量的影响,有学者对Fmax进行修正,修正后的内梅罗数学表达式为[24]:
(8)
(9)
式中:FW为权重值最大的评价指标值。
针对传统内梅罗指数法的缺陷,有学者引入权重系数,从而增强评价结果的合理性[25-26]。改进后的数学表达式为:
(10)
上述改进方法未同时考虑最值和均值的问题,势必在一定程度上会影响水质评价结果的正确性,故同时对最值和均值做出修正[12]。
(11)
(12)
(13)
式中:F′为为前n项组分的评分均值;m为的项数。
本文以北戴河国家湿地公园2018年9月水质监测数据为基础,以Ⅲ类水质作为评价基准。在园区内设置22个水样取样点,有针对性地选取项水质指标为评价指标,其中Chlorophyll-a仅作为营养状态评价指标。相关项目的检测按照《水和废水监测分析方法》(第四版)[27]操作,监测数据如表1所示。
本文以GB 3838—2002《地表水环境质量标准》中的Ⅲ类水为计算基准,污染因子的权重计算结果见表2。计算结果表明,TP指标对水体环境质量的影响最大,次之,COD的影响最小。
表1 北戴河国家湿地公园2018年9月水质监测数据和评价结果
Table 1 Water quality monitoring data and evaluation results of Beidaihe National Wetland Park
样点ρ(NH+4-N)/(mg·L-1)ρ(TN)/(mg·L-1)ρ(TP)/(mg·L-1)ρ(COD)/(mg·L-1)Chlorophyll-a/(μg·L-1)P传PⅠPⅡPⅢPⅣIwq营养状态指数10.1363.42 0.188 4086.52.677Ⅴ2.559Ⅴ1.922Ⅴ2.508Ⅴ2.324Ⅴ4.621 Ⅳ72.122.1563.86 0.346 84170.63.324Ⅴ3.324Ⅴ2.976Ⅴ3.212Ⅴ2.947Ⅴ6.443劣Ⅴ,不黑臭 78.730.0822.91 0.124 5262.82.332Ⅴ2.150Ⅴ1.662Ⅴ2.104Ⅴ2.098Ⅴ4.321 Ⅳ67.940.0554.23 0.218 203.73.196Ⅴ3.148Ⅴ2.193Ⅴ3.089Ⅴ2.446Ⅴ4.121 Ⅳ59.550.0272.86 0.096 1640.02.152Ⅴ2.092Ⅴ1.392Ⅳ2.051Ⅴ2.092Ⅴ3.410 Ⅲ64.560.1911.50 0.140 4461.41.755Ⅴ1.644Ⅴ1.308Ⅳ1.632Ⅴ1.527Ⅴ4.020 Ⅳ65.270.0271.26 0.156 3659.71.445Ⅳ1.378Ⅳ1.140Ⅳ1.387Ⅳ1.291Ⅳ3.720 Ⅲ64.880.0271.56 0.144 4080.41.606Ⅴ1.510Ⅴ1.227Ⅳ1.503Ⅴ1.438Ⅳ3.920 Ⅲ66.890.0551.50 0.238 5625.22.209Ⅴ2.123Ⅴ1.718Ⅴ2.151Ⅴ1.808Ⅴ4.331 Ⅳ63.710—1.35 0.156 7635.62.884Ⅴ2.742Ⅴ1.929Ⅴ2.743Ⅴ2.319Ⅴ4.121 Ⅳ62.71102.55 0.214 20149.31.979Ⅴ1.972Ⅴ1.518Ⅴ1.955Ⅴ1.651Ⅴ3.820 Ⅲ 73.9120.5191.13 0.140 6011.22.323Ⅴ2.187Ⅴ1.614Ⅴ2.178Ⅴ1.868Ⅴ4.321 Ⅳ 56.1130.1091.29 0.164 2431.71.094Ⅳ1.069Ⅳ0.960Ⅲ1.081Ⅳ1.056Ⅳ3.520 Ⅲ 62.2140.1911.14 0.136 3630.81.440Ⅳ1.362Ⅳ1.106Ⅳ1.361Ⅳ1.253Ⅳ3.820 Ⅲ 60.6150.0271.20 0.148 3222.81.295Ⅳ1.238Ⅳ1.040Ⅳ1.247Ⅳ1.174Ⅳ3.420 Ⅲ59.9160.5732.10 0.152 4052.51.768Ⅴ1.625Ⅴ1.395Ⅳ1.579Ⅴ1.609Ⅴ4.521 Ⅳ 66.5170.2181.62 0.112 4457.11.755Ⅴ1.627Ⅴ1.270Ⅳ1.605Ⅴ1.529Ⅴ4.121 Ⅳ 64.1180.1091.58 0.200 4852.91.921Ⅴ1.829Ⅴ1.501Ⅴ1.838Ⅴ1.590Ⅴ4.221 Ⅳ 66.5190.0820.49 0.036 168.00.629Ⅱ0.586Ⅱ0.442Ⅱ0.580Ⅱ0.321Ⅱ2.400 Ⅱ 44.0200.3820.87 0.064 204.40.841Ⅲ0.765Ⅲ0.652Ⅲ0.742Ⅲ0.765Ⅲ3.200 Ⅲ 46.8210.2180.60 0.076 1610.80.667Ⅲ0.630Ⅱ0.547Ⅱ0.626Ⅱ0.396Ⅱ2.800 Ⅱ49.9220.1640.52 0.052 418.70.448Ⅰ0.419Ⅰ0.376Ⅱ0.411Ⅰ0.273Ⅰ2.100 Ⅱ49.9
以GB 14848—1993中的Ⅲ类水为计算基准,计算传统和修正后的内梅罗指数等级划分,从而确定计算的水质类别。例如,Ⅰ类水环境中待评价污染因子的内梅罗指数计算如下:将Ⅰ类水待评价水质指标的临界值与Ⅲ类水水质指标的临界值相比[24],结果见表3。
表2 污染因子权重值Wi
Table 2 Pollution factors’ weight value
指标Ⅲ类水标准/(mg·L-1)WiNH+4-N1.00.142TP0.20.709TN1.00.142COD200.007
表3 Ⅰ类水质指标临界值与Ⅲ类水质指标临界值的比值
Table 3 Ratios of the critical value of the class Ⅰ
index to the critical value of the class Ⅲ
水质指标NH+4-NTPTNCOD临界值0.150.10.20.75
由表3可计算出:为0.3,最大值Fmax为0.75,由式(1)计算得P传为0.57,即只要待评价水体内梅罗指数<0.57,则水质为Ⅰ类水体。同理可得其他水质类别及修正内梅罗指数,如表4所示。
表4 内梅罗指数类别
Table 4 Categories of Nemero index
水质类别ⅠⅡⅢⅣⅤP传P<0.5710.571≤P<0.6630.663≤P<11≤P<1.51.5≤PPⅠP<0.5380.538≤P<0.6380.638≤P<11≤P<1.51.5≤PPⅡP<0.3680.368≤P<0.5950.595≤P<11≤P<1.51.5≤PPⅢP<0.5350.535≤P<0.6370.637≤P<11≤P<1.51.5≤PPⅣP<0.2800.280≤P<0.4430.443≤P<11≤P<1.51.5≤P
注:Ⅳ和Ⅴ类内梅罗指数相同,是因为的Ⅳ和Ⅴ类临界值比上Ⅲ类临界值的数值完全相等,若增加其他参评水质指标,则可避免出现这种情况。
应用修正前后的内梅罗指数法对北戴河国家湿地公园22个水样点的水质进行评价,评价结果见表1。根据评价结果可知:北戴河国家湿地公园水质受污染较为严重,其水质主要为Ⅳ及Ⅴ类,未达到本湿地水体功能区规划要求;湿地水体主要超标因子为TN、COD和Chlorophyll-a,TP也存在个别样点超标现象;修正前后的5种内梅罗指数法对水质评价结果有所差别,其中,修正内梅罗指数法Ⅱ的评价结果(达标率22.7%)要优于其他4种内梅罗指数评价结果(达标率为18.2%)。
综合水质标识指数法是徐祖信在单因子水质标识指数法的基础上提出的一种全新的水质评价方法,该法不仅可以定性、定量合理评价水质,还可比较劣Ⅴ类水体污染的严重程度[28],评价方法见参考文献[28-29]。
综合水质标识指数法评价结果显示(表1),北戴河国家湿地公园22个样点水质达标率为50.0%,远高于内梅罗指数法的评价结果;综合水质标识指数法可以评价出劣Ⅴ类水质(样点2)。将修正前后的内梅罗指数和综合水质标识指数进行相关行分析(图1),其相关系数分别为0.72(P传)、0.68(PⅠ)、0.84(PⅡ)、0.67(PⅢ)、0.74(PⅣ)。其中,内梅罗指数Ⅱ和综合水质标识指数具有较强的线性相关性,相关系数达到0.84。内梅罗指数法和综合水质标识指数法都具有定性、定量评价功能,都能综合反映水质污染程度,在对同一水体进行评价时,两者评价结果之间理应存在一定程度的相关性,虽鲜有文献报道,但从本文结果分析得,两者具有较强的线性相关性。
P传; PⅠ; PⅡ; PⅢ; PⅣ。
图1 内梅罗指数与综合水质标识指数的相关性
Fig.1 Correlation between Nemerow index and
comprehensive water quality index
Chlorophyll-a是湿地水环境中非常重要的水质指标,北戴河国家湿地公园水质污染也主要是由富营养化引起。根据《中国湖泊环境》(第一册)[30]中的综合营养状态评价方法,对北戴河湿地公园22个水样点的水质进行营养状态评价,评价结果见表1。
营养状态指数表明,大部分样点的水质处于中度富营养状态,样点1、2、11的水质达到重富营养化的状态,说明北戴河国家湿地公园水体富营养化较为严重。采用SPSS.22统计软件对修正前后内梅罗指数和营养状态指数评价结果进行皮尔森相关性分析(表5),结果表明:置信水平都在0.01以上,传统内梅罗指数(R=0.704)、内梅罗指数Ⅰ(R=0.698)、内梅罗指数Ⅱ(R=0.762)、内梅罗指数Ⅲ(R=0.697)、内梅罗指数Ⅳ(R=0.770)与营养状态指数都显著相关,且内梅罗指数Ⅳ的相关性最大,其次为内梅罗指数Ⅱ。根据对北戴河国家湿地公园水质监测,其主要超标指标为TN、Chlorophyll-a和COD,而TN和Chlorophyll-a也刚好是综合营养状态指数评价指标,因此,营养状态指数和内梅罗指数之间理应具有一定的相关性,对比结果则印证了这点。徐彬等[7]的对比结果也表明内梅罗指数和营养状态指数具有很强的相关性。
表5 内梅罗指数与营养状态指数相关性分析
Table 5 Correlation analysis between Nemero index
and nutritional status index
内梅罗指数P传PⅠPⅡPⅢPⅣ相关系数R20.704**0.698**0.762**0.697**0.770**P0.000250.000310.000030.000320.00003
注:**表明(P<0.01),在双侧水平上极显著相关。
1)参评水质指标的权重计算结果表明:北戴河国家湿地公园的水环境质量受TP指标的影响最大,次之,COD影响最小,水质属于氮限制状态。
2)内梅罗指数评价结果表明:不同修正方法修正的内梅罗指数法在评价北戴河国家湿地公园水质结果上存在差异。数据表明,修正内梅罗指数法Ⅱ对水体评价结果要稍优于其他内梅罗指数评价结果,样点水质达标率为22.7%,其余为18.2%;与综合水质标识指数法对比结果表明,修正内梅罗指数法Ⅱ的评价结果和综合水质标识指数的相关性最强,相关系数达到0.84。
3)营养状态指数评价结果表明:北戴河国家湿地公园水质营养化较为严重;营养状态指数与传统内梅罗指数、内梅罗指数Ⅰ、内梅罗指数Ⅱ、内梅罗指数Ⅲ、内梅罗指数Ⅳ的相关性分别为0.704、0.698、0.762、0.697、0.770,均表现出一定的相关性。
4)通过与综合水质标识指数和营养状态指数的相关性分析,修正内梅罗指数法Ⅱ在北戴河国家湿地公园水环境质量评价中可能更为合理可靠,并且能较好地反映各污染因子对水质影响大小,评价结果拟合程度也较高。
[1] Ji X L, Randy A Dahlgren, Zhang M H. Comparison of seven water quality assessment methods for the characterization and management of highly impaired river systems[J]. Environmental Monitoring & Assessment, 2016,188(1):1-16.
[2] Uddameri V, Honnungar V, Hernandez E A. Assessment of groundwater water quality in central and southern Gulf Coast aquifer, TX using principal component analysis[J]. Environmental Earth Sciences, 2014,71(6):2653-2671.
[3] Li D Q, Ji L, Di Y M, et al. The spatial-temporal variations of water quality in controlling points of the main rivers flowing into the Miyun Reservoir from 1991 to 2011[J]. Environmental Monitoring & Assessment, 2016,188(1):42-54.
[4] Tiri Ammar, Belkhiri Lazhar, Mouni Lotfi. Evaluation of surface water quality for drinking purposes using fuzzy inference system[J]. Groundwater for Sustainable Development, 2018,6:235-244.
[5] 卢俊平, 张晓晶, 马太玲, 等. 基于可变模糊数学模型的大河口水库富营养化评价[J]. 环境工程, 2016,34(11):136-140.
[6] Charulatha G,Srinivasalu S, Uma Maheswari O, et al. Evaluation of ground water quality contaminants using linear regression and artificial neural network models[J]. Arabian Journal of Geosciences, 2017,10(6):1-9.
[7] 徐彬, 林灿尧, 毛新伟. 内梅罗水污染指数法在太湖水质评价中的适用性分析[J]. 水资源保护, 2014,30(2):38-40.
[8] 金士博, 张孟威. 内梅罗污染指数公式与漓江水质评价[J]. 环境科学, 1980(2):1-7.
[9] 刘萌, 陈世俭. 基于内梅罗指数与主成分分析的洪湖地区地下水水质评价[J]. 华中师范大学学报(自然科学版), 2016,50(4):633-640.
[10] Zheng Q L, Ma T, Wang Y Y, et al. Hydrochemical characteristics and quality assessment of shallow groundwater in xincai river basin, northern China[J]. Procedia Earth and Planetary Science, 2017,17:368-371.
[11] 马骁, 张永祥, 王昊. 改进的尼梅罗污染指数法在地下水水质评价中的应用[J]. 环境工程, 2017,35(3):158-162.
[12] 李小丽, 黎小东, 敖天其. 改进内梅罗指数法在西充河水质评价中的应用[J]. 人民黄河, 2016,38(8):65-68.
[13] 吴彬, 臧淑英, 那晓东. 灰色关联分析与内梅罗指数法在克钦湖水体重金属评价中的应用[J]. 安全与环境学报, 2012,12(5):134-137.
[14] Al-Othman,Ahmed A. Evaluation of the suitability of surface water from Riyadh Mainstream Saudi Arabia for a variety of uses[J]. Arabian Journal of Chemistry, 2015(1):1-7.
[15] 韩术鑫, 王利红, 栾玲玉, 等. 基于浅层地下水特征因子构建的内梅罗指数法在东营市海水入侵评价中的应用[J]. 海洋环境科学, 2018,37(2):294-298.
[16] 罗嗣卿, 贾子书. 改进内梅罗指数法在东方红湿地水质评价中的应用[J]. 黑龙江大学自然科学学报, 2018,35(1):15-21.
[17] 吕俊芳, 王书锦, 刘云根, 等. 罗时江河口湿地水质季节变化特征及评价[J]. 环境工程, 2016,34(增刊1):845-849,877.
[18] Chi Thanh Vu, Lin Chitsan, Shern Chienchuan, et al. Contamination, ecological risk and source apportionment of heavy metals in sediments and water of a contaminated river in Taiwan[J]. Ecological Indicators, 2017,82:32-42.
[19] Yan N, Liu W B, Xie H T, et al. Distribution and assessment of heavy metals in the surface sediment of Yellow River, China[J]. Journal of Environmental Sciences, 2016,39:45-51.
[20] Gąsiorek M, Kowalska J, Mazurek R, et al. Comprehensive assessment of heavy metal pollution in topsoil of historical urban park on an example of the Planty Park in Krakow (Poland)[J]. Chemosphere, 2017,179:148-158.
[21] Adamu C I, Nganje T N, Edet Aniekan. Heavy metal contamination and health risk assessment associated with abandoned barite mines in Cross River State, southeastern Nigeria[J]. Environmental Nanotechnology, Monitoring & Management, 2015,3:10-21.
[22] 李秉辉, 李海燕. 海绵公园土壤重金属现状调查与评价[J]. 环境工程, 2016,34(增刊1):1016-1020,1024.
[23] 韩术鑫, 王利红, 赵长盛. 内梅罗指数法在环境质量评价中的适用性与修正原则[J]. 农业环境科学学报, 2017,36(10):2153-2160.
[24] 蔡晔, 林怡雯, 李月娥, 等. 利用改进的内梅罗指数法模型评价苏州市内外城河水质[J]. 化学分析计量, 2015,24(2):84-87.
[25] 王宇博, 梁秀娟, 乔雨. 2种方法在地下水水质评价中的应用:以吉林市为例[J]. 节水灌溉, 2014(5):34-37.
[26] 刘璐瑶, 冯民权. 几种水质评价方法在涑水河的应用与比较研究[J]. 黑龙江大学工程学报, 2017,8(3):6-14.
[27] 国家环境保护总局《水和废水监测分析方法》编委会. 水和废水监测分析方法[M]. 4版.北京:中国环境科学出版社,2002.
[28] 徐祖信.我国河流综合水质标识指数评价方法研究[J].同济大学学报(自然科学版),2005,33(4):482-488.
[29] 徐祖信.我国河流单因子水质标识指数评价方法研究[J].同济大学学报(自然科学版),2005,33(3):321-325.
[30] 金相灿. 中国湖泊环境(第一册)[M]. 北京:海洋出版社,1995.