1. 引言
捕鱼作业区域作为近岸海域渔船活动的主要区域,渔业生产过程是造成该区域污染的一个主要来源 [1],对捕鱼作业区水环境质量进行科学、综合、客观的评价十分重要 [2] [3]。近年来,大量的数学方法应用于水环境质量评价中,如水质指数法(WQI) [4]、主成分分析法(PCA) [5] [6]、主因子分析法(PFA) [7] [8]、模糊综合评价(FSE) [9] [10] [11] 等,这些方法在构建水环境评价体系的应用过程中往往忽略了水环境质量指标和其他相关调研因素的相关性。灰色系统理论是20世纪80年代创立的一门基于数学理论的系统工程新兴学科 [12],灰色关联模型主要根据各因素之间发展趋势的几何相似程度,确定各评价因子各指标间强弱、大小和次序和关系,目前己经成为水环境质量评价中应用较广泛的一种量化方法 [13] [14] [15]。但是,传统的灰色关联法存在一些缺陷,例如评价步骤不规范,计算过程极其复杂,评价结果仅限于给出水质污染现状,而对其污染物影响程度难以进行判断。
为了能够综合反应捕鱼作业区域水环境质量,客观说明各评价因子之间的关系和优劣程度,本文将灰色关联法进行改进,根据监测数据和评价内容确定评价矩阵构造判别矩阵,通过对评价矩阵的归一化和标准化优化了计算过程,同时构建母序列并计算各个指标与母序列的灰色关联度,通过计算每个指标的权重确定污染物影响程度,最后对评价对象的得分归一化后进行优劣排名,将定性分析和定量分析相结合,使评价结果更为公正客观。在将该模型应用到天津地区主要近岸捕鱼作业区域水质监测指标的评价过程中,充分考虑了各评价指标的内在联系和相互作用,对区域水环境做出了综合评价。在将该模型与传统的综合指数法对比结果可以看出改进的灰色关联法模型表现出了较好的稳定性和有效性,更加全面反应了水质综合情况。
2. 研究区域概况
天津市位于华北平原东北部,所邻近的渤海湾海域是重要的海洋经济水产物种的繁育区,海岸线全长153.67 km,传统海域面积约3000 km2,其中渔业用海共计18 km2。根据农业农村部《渤海渔港认定标准(试行)》,2019年12月,滨海新区人民政府对渔港重新进行了认定,分别是北塘渔港、东沽渔港、大神堂渔港、蔡家堡渔港和天津中心渔港。渔船产生的含油污水、渔民排放生活废水是作业区水环境污染一个重要来源,在一定程度上也会对近岸海域水质状况造成污染。因此,针对5个渔港所辖的近岸主要作业区进行了综合水质评价研究。
3. 材料与方法
3.1. 样品采样及分析方法
天津开海时间为每年的9月1日,因此在捕鱼期对渔民主要作业区域进行海水水质样品的采集,根据《海水水质标准》(GB3097-1997)和《近岸海域环境监测规范》(HJ442-2008)确定pH、水温、溶解氧、悬浮物质、粪大肠菌群、化学需氧量(COD)、生化需氧量(BOD5)、无机氮、非离子氨(以N计)、活性磷酸盐(以P计)、石油类、铅、镉、铜等14项指标作为监测指标,样品根据《海洋监测规范第3部分:样品的采集、贮存与运输》(GB 17378.3-2007)采集表层水样,2021年10月对每个渔港临近的主要捕鱼作业区选择2个监测点位,主要作业区域位于近岸7~10海里以内(即13~18公里以内)采样点位分布情况如图1所示,监测方法依照《海洋监测规范第四部分海水分析》(GB 17378.4-2007)执行。
Figure 1. Monitoring points of fishing operation areas in Tianjin
图1. 天津主要捕鱼作业区域监测点位
3.2. 基于灰色关联的综合评价步骤
灰色关联度分析是指用灰色关联度来描述因素间关系的强弱、大小和次序。它是一种相对性的排序分析,根据各序列间的关联度判断联系是否紧密,关联度越大,说明被评价受污染地区与区域内受污染程度最高的地区越相似。因此,关联度的大小顺序,就是受污染程度强弱的次序。灰色综合评价方法将定性分析和定量分析相结合的综合评价方法,较好地解决了评价指标难以准确量化和统计的问题,已被广泛应用于水质评价中。
灰色关联度分析基本评价步骤如下:
1) 构造判别矩阵。根据监测数据和评价内容确定评价矩阵
,其中
为n个指标中第m个评价对象。
2) 对指标进行了正向化和标准化,并确定指标母序列。
监测过程评价指标共分为了3类,其中pH和水温属于区间型,溶解氧属于极大型指标,其余监测因子属于及小型监测指标。为了以消除各指标量纲,缩小变量的范围,需要对区间型和极小型指标进行正向化处理,将其转化为极大型指标,再对正向化后的指标进行标预处理,以简化计算过程。
极小型指标正向化:
(1)
区间型指标正向化:
(2)
其中,
是将正向化后的指标与原来极大型指标共同构建成正向化矩阵
,然后对其进行预处理,得到预处理后矩阵
:
(3)
其中
为第j个指标正向化后的平均值。
3) 构建母序列并计算各个指标与母序列的灰色关联度。
首先,求出每一行最大值构成列向量表示母序列,子序列就是预处理后的数据矩阵。
(4)
其中,母序列为
,
,即每一行最大值,子序列为
为预处理后矩阵Z。
然后,利用公式(5)计算各个指标与母序列的灰色关联度
:
(5)
其中,
其中c为两极最小差:
;
d为两级最大差:
;
为分辨系数,一般取
。
4) 计算每个指标的权重。
(6)
其中
为每个指标的权重。
5) 计算评价对象的得分,归一化后进行优劣排名。
,(
) (7)
(8)
其中
为评价对象的得分,
为归一化后得分。
以上5步构成基于灰色关联的水环境质量综合评价方法。
4. 研究区水环境污染状况评价
4.1. 单因子水质评价结果
作业区域水环境质量评价采用单因子评价法,评价与计算方法如式(9),结果见表1:
(9)
式中:
——第i类污染物的污染指数;
——第i类污染物的平均浓度;
——第i类污染物的标准值。
Table 1. Assessment results of fishing operation areas water environment by single factor method
表1. 主要作业区水环境单因子质量评价结果
通过单因子对比结果可知,各主要作业区发现,北塘渔港、中心渔港和蔡家堡渔港作业区水环境质量以劣四类为主,主要超标的污染因子包括化学需氧量、生化需氧量、无机氮、非离子氨和活性磷酸盐。其中化学需氧量浓度范围在2.64~11.8 mg/L之间,生化需氧量为1.0~6.9 mg/L,无机氮为0.093~1.167 mg/L,非离子氨为0.004~0.070 mg/L,活性磷酸盐为0.002~0.074 mg/L。
4.2. 灰色关联评价结果分析
首先根据14类污染物构造判别矩阵X,在根据公式(1)~(3)对各项指标进行正向化和标准化,标准化后的结果如表2所示,各污染指标权重wj值,如图2所示。
Table 2. Standardization results of each evaluation index
表2. 各评价指标标准化结果
构建母序列并计算各个指标与母序列的灰色关联度
,并根据公式
计算权重
,结
果如表3所示。
Table 3. Grey correlation degree y j and weight values w j
表3. 灰色关联度
和权重
Figure 2. The weight values
of each pollution indicators
图2. 各污染指标权重
值
从权重矩阵可以看出作业区主要的污染影响因子为化学需氧量、生化需氧量、无机氮、非离子氨和活性磷酸盐,这几类污染因子均超标,这与单因子评价法超标情况一致,表明方法可靠。最后根据公式(7)、(8)计算评价对象的得分
,归一化后进行优劣排名如表4所示。
Table 4. Evaluation results of improved grey correlation method and composite index method
表4. 灰色关联法和综合指数法评价结果
为了证明方法的有效性和评价结果的准确定,将灰色关联法与综合指数法的结果进行对比,结果如表4和图3所示。图3是对各个监测点位进行排名的对比,从图中的结果可以看出,两种方法的排名结果变化趋势基本一致,即评价结果较为可靠性较高。
根据评价结果可以看出,7、8点位所临近的中心渔港附近的捕鱼作业区水环境质量最好,而5、6和1、2所临近的北塘渔港和大神堂渔港水质较差,造成这种现象的主要原因是,北塘渔港和大神堂渔港船只较多,且以捕鱼作业生产为主,渔民活动十分密集,而中心渔港以旅游为主,船只较少,渔民活动仅限于渔业旅游,且渔船各种设施配套完善,因此,该区域作业区水质相对较好。
Figure 3. Comparison of ranking results of improved grey correlation method and composite index method
图3. 灰色关联法和综合指数法排名结果对比
5. 结论及建议
本文将灰色关联法进行改进,根据监测数据和评价内容确定评价矩阵,通过归一化和标准化优化计算过程,构建母序列计算各个指标与母序列的灰色关联度,通过计算每个指标的权重确定污染物影响程度,最后对评价对象的得分归一化后进行优劣排名,并应用到天津地区主要近岸捕鱼作业区域水质监测指标的评价过程中,将定性分析和定量分析相结合,使评价结果更为公正客观,在将该模型与传统的综合指数法对比结果可以看出改进的灰色关联法模型表现出了较好的稳定性和有效性,更加全面反应了水质综合情况。
通过评价和调研结果,可以看出,由于大部分渔船配有简易厨房和卫生间,但没有生活污水收集器,仅有简单的水桶等,存在生活污水直排入海的情况,生活污水的主要因素和评价结果显示的主要污染影响因子一致,主要包括化学需氧量、生化需氧量、无机氮、非离子氨和活性磷酸盐。
基金项目
天津市科技计划项目(18ZXSZSF00130)。
NOTES
*通讯作者。