近年来,沿海地区建立了越来越多的配有烟气海水脱硫系统(SFGD)的燃煤发电厂[1]。这些电厂利用海水弱碱性吸收烟气中SO2的同时,也会将部分重金属洗脱一起携带入海[2],从而导致附近海域重金属含量增加[3-6]。重金属由于具有生物毒性强、难生物降解以及易在生物体内累积等性质,可对海洋生物造成严重威胁,被认为是海洋环境的重要污染物[7-8]。例如,Zn抑制黑海胆(Arbacia lixula)生长发育的3 d-EC50为21 μg·L-1 [9],导致大凤螺(Lobatus gigas)幼虫死亡的4 d-LC50为8.9 μg·L-1[10];Cu的暴露浓度为1 μg·L-1时就会对微藻产生毒性作用[11],18 μg·L-1时即可造成紫贻贝(Mytilus edulis)DNA损伤[12]。而郭娟等[13]指出脱硫海水排水口附近海域中Zn和Cu的检测范围分别为8.63~123.39 μg·L-1和19.24~72.23 μg·L-1。因此,有必要开展脱硫海水长期排放对排污口附近海域的生态风险评价工作。目前燃煤电厂脱硫海水重金属的研究主要集中于具有挥发性质的Hg在烟气和水环境介质中的迁移转化[14-17],而对于其他重金属,如As、Cd、Cr、Cu、Pb和Zn的长期排放对排污口附近海域的影响的研究尚未见报道。
水环境中的重金属,由于吸附、水解和沉淀的共同作用,只有一小部分的自由金属离子溶解在水中,大部分都沉积于沉积物中[18]。因此,对重金属的生态风险评价主要集中在沉积物上。目前,已有多种污染指数法可用于沉积物中重金属的风险评价[19]。其中,沉积物质量基准(SQGs)[20-21]、潜在生态危害指数(potential ecological risk index, PERI)[22-23]、地累积指数(Igeo)[1,24]、富累积因子(EF)[25]、污染严重指数(CSI)[26]等均是较为常用的方法。Ranjbar Jafarabadi等[8]采用SQGs、PERI、CSI等10种方法评价了波斯湾10个珊瑚岛沿海表层沉积物中13种重金属的生态风险,其中基于EF和PERI的结果显示V、Ni、Hg和Cd为中等污染水平,需要引起关注。然而,上述污染指数法通常仅使用(或联合使用)数据的平均值、保守值(90%或95%分位数)和最大值来获得风险的一个平均或保守的估计[25,27]。这种将监测数据压缩为一个单点值的方法会导致生态风险评价中信息的丢失,产生不确定性,从而造成对风险的高估或者低估。因此,美国国家科学院和环境保护局推荐使用Monte Carlo模拟方法来量化风险评价中的可变性和不确定性[28-29]。这种随机模拟技术通过产生大量的符合一定规律的随机数,将其带入到有效的风险评价模型,进而得到风险的概率分布,为决策者提供更多的信息。
燃煤电厂附近海域重金属浓度一般是通过浓度增量(根据电厂及SFGD的运行参数)叠加环境背景值或对电厂运行后的短期检测获得的[30]。然而,由于估计值或短期检测值的不确定因素较多,导致以其为基础的风险评价结果具有较大的不确定性。因此,本研究连续3年跟踪监测了青岛某燃煤电厂排污口附近海域沉积物中As、Cd、Cr、Cu、Hg、Pb和Zn的浓度。基于重金属的长期暴露数据,采用Monte Carlo模拟技术与SQGs和PERI相结合的概率方法评价了排污口附近海域沉积物中重金属的生态风险。本研究能够为燃煤电厂排口附近海域的污染控制和生态风险减缓提供更多的科学依据。
在青岛胶州湾某燃煤电厂排污口附近海域共设置了12个站位(图1)。利用采泥器分别于2012年12月26日和2013年6月21日的低潮期采集表层2 cm处的沉积物混合样,去除上覆水后,装入密封的聚乙烯袋中带回实验室。样品运回实验室后,晾干,研磨至完全通过160目(96 μm)的尼龙筛,然后存储到预先准备好的干净样品袋中以待进一步分析。此外,本研究还引用了王云鹏等[31]于2011年12月23日和2012年6月5日对该排污口附近表层沉积物的重金属监测数据。其站位设置与采样方法均与本研究相同。2011年12月23日未在4#站位采集到沉积物样本,2013年6月21日未在1#和4#站位采集到沉积物样本。沉积物样品采集回实验室后,置于冷冻干燥机中进行冻干,并研磨过160目筛,之后密封冷藏保存。
图1 调查站位设置
Fig. 1 Map of sampling stations near the outlet of the coal-fired power plant
取沉积物样品0.04 g置于Teflon容器中,加入0.8 mL HNO3(V∶V=1∶1)和0.8 mL HF,再加入HClO4(V∶V=1∶3)后,密封超声震荡60 s,然后在100 ℃下恒温加热蒸干。再加入0.8 mL HNO3(V∶V=1∶1),在100 ℃下恒温加热蒸干。之后加入0.8 mL的HF和HClO4,密封置于高压釜中,在170 ℃烘箱内加热。样品蒸干之后再加入4 mL 4 mol·L-1的HNO3,在170 ℃烘箱内加热4 h。最后用3%(V∶V)HNO3稀释并移至50 mL的容量瓶中,加入Rh-Re内标液,以1%(V∶V)HNO3至稀释至50 mL以备分析使用。
重金属分析采用美国安捷伦公司生产的等离体质谱仪ICP-MS(Agilent 7500CX),标准物质为安捷伦公司多元素混合标准溶液(批号:5183-4688)。通过该方法,金属元素As、Cd、Cr、Cu、Hg、Pb和Zn的最低检出限分别为0.2、0.02、1、0.2、0.01、0.5和2 mg·kg-1;所有金属样品的回收率均在90%~110%范围内,重复样品的相对标准差均<10%。
SQGs方法是基于北美沉积物生物效应数据库构建的多阈值型基准,广泛用于海洋沉积物中污染物的风险评价[20,32-33]。Long等[20]将导致10%和50%有害生物效应的污染物浓度分别称为效应范围低值(effect range low, ERL)和效应范围中值(effect range median, ERM),随后MacDonald等[21]又对该阈值做了修改,给出了一对更为严格的阈值:临界效应水平(threshold effect level, TEL)和可能发生效应水平(probable effect level, PEL)。当污染物的环境浓度EEC(estimated environmental concentration)
表1 7种重金属的不同类型海洋沉积物质量基准值
Table 1 Different types of marine sediment quality guideline values for the seven heavy metals (mg·kg-1)
重金属Heavy metalsTELPELERLERM砷 As7.2441.68.270镉 Cd0.684.211.29.6铬 Cr52.316081370铜 Cu18.710834270汞 Hg0.130.70.150.71铅 Pb30.211246.7218锌 Zn124271150410
注:数据引自文献[21];TEL指临界效应水平;PEL指可能发生效应水平;ERL为效应范围低值;ERM为效应范围中值。
Note: The data come from Reference [21]; TEL means threshold effect level; PEL means probable effect level; ERL means effect range low; ERM means effect range median.
基于水生生态系统的敏感性依赖于其生产力的假设,Hakanson[22]提出了PERI方法。该方法结合监测浓度与生态效应、环境效应和毒性的关系,对沉积物中重金属的污染程度进行评估。其具体计算公式如下:
(1)
(2)
(3)
式中:RI为多种金属的综合潜在生态危害指数,为金属i的潜在生态危害指数,
为金属i的毒性响应系数,
为金属i的污染系数,
为重金属i的实测浓度,
为金属i的沉积物背景值。各类重金属的背景值是该海域未受到影响时的测量值,本研究参考王云鹏等[31]所给出的胶州湾沉积物中重金属背景值和毒性响应系数(表2)。基于上述参考值,重金属潜在生态危害指数对应的危害程度分级标准见表3。
表2 胶州湾沉积物中重金属浓度背景参考值及毒性响应系数
Table 2 The background values and toxic-response factors of heavy metals in sediments of Jiaozhou Bay
砷As镉Cd铬Cr铜Cu汞Hg铅Pb锌ZnCir/(mg·kg-1)7.750.6658.1660.340.02426.2678.28Tir1030254051
注:为金属的沉积物背景值;
为金属的毒性响应系数。
means toxic response factor for heavy metals.
表3 沉积物中重金属潜在生态危害指数对应的危害程度分级标准
Table 3 Criteria for degrees of the ecological risk of heavy metal in sediments
Eir范围Range of Eir单个因子危害程度Risk degree of individual factorRI范围Range of RI多个因子综合危害程度Risk degree of multiple factorsEir<40轻微 LowRI<150轻微 Low40≤Eir<80中等 Moderate150≤RI<300中等 Moderate80≤Eir<160强 Considerable300≤RI<600强 Considerable160≤Eir<320很强 HighRI≥600很强 HighEir>320极强 Extreme
注:为金属的潜在生态危害指数;RI为多种金属的综合潜在生态危害指数。
Note: means potential ecological risk index for a certain heavy metal; RI means the potential ecological risk index of heavy metals.
本研究采用Monte Carlo模拟技术与传统的SQGs和PERI相结合的概率方法评价表层沉积物中重金属的生态风险。Monte Carlo模拟技术是将模型或者方法中需要输入的关键变量(如污染物的浓度)视为符合某种概率分布的随机变量,通过重复抽样,带入到模型或者方法中得到输出结果的概率分布[34]。本研究采用log-normal、log-logistic和Weibull分布分别对每种重金属的监测数据进行分布拟合,根据Kolmogorov-Smirnov检验选择最优模型:P值越接近于1,拟合效果越好[35]。总体上,log-logistic的拟合效果要优于其他2种分布。基于log-logistic,运用软件Matlab v.9.3 (MathWorks Inc., Natick, Massachusetts)对每个数据集进行10万次的随机抽样。然后再结合SQGs和PERI法,得到超过各个SQGs标准值的概率以及和RI的概率分布,进而得到重金属所处各个风险水平的概率等更为详细的信息。
4次监测排污口附近海域表层沉积物中As、Cd、Cr、Cu、Hg、Pb和Zn的浓度范围分别为1.54~7.33、0.08~0.72、12.16~54.04、5.08~33.10、0.02~0.08、8.37~35.91和14.58~85.01 mg·kg-1,平均浓度分别为3.54、0.36、26.22、18.04、0.04、24.37和47.21 mg·kg-1。7种重金属4次监测的平均浓度均符合中国海洋沉积物质量标准(GB 18668—2002)的一类标准值。各个监测时间的详细结果见表4。排污口附近下游区域(2#~11#站位)沉积物中重金属的平均浓度大体上都高于排水口上游的1#站位和胶州湾对照点12#站位。如2012年12月,2#~11#站位As、Cd、Cr、Cu、Pb和Zn的平均浓度分别为4.25、0.34、37.51、22.03、26.90和57.61 mg·kg-1,高于1#站位的2.86、0.10、17.03、15.29、23.75和40.31 mg·kg-1,高于12#站位的3.81、0.08、23.34、8.58、16.01和50.53 mg·kg-1,Hg的浓度变化不明显。
表4 电厂排污口附近海域表层沉积物重金属含量
Table 4 Statistics of concentrations of the heavy metals in the marine surface sediment near the outlet of the power plant (mg·kg-1)
重金属Heavy metal2011.12a2012.06a2012.122013.06范围Range平均值Mean标准差SD范围Range平均值Mean标准差SD范围Range平均值Mean标准差SD范围Range平均值Mean标准差SD砷 As3.23~7.334.991.502.04~3.282.710.363.33~5.044.210.541.54~2.992.250.40镉 Cd0.11~0.550.260.120.25~0.720.520.150.08~0.680.310.190.09~0.610.320.17铬 Cr14.75~42.7824.417.0715.76~38.4123.256.4123.34~54.0436.228.7212.16~45.7820.3010.02铜 Cu5.35~31.4018.936.706.43~18.2513.873.928.58~33.1020.818.615.08~25.7118.716.45汞 Hg0.03~0.050.040.0080.03~0.060.0450.010.02~0.050.030.0080.04~0.080.0540.013铅 Pb8.37~25.0917.984.7717.46~35.9125.765.5315.03~35.6425.916.7013.48~33.7127.516.28锌 Zn24.13~85.0155.1515.8419.84~52.7635.849.7741.66~83.4956.9613.7214.58~63.8141.0317.79
注:a表示引自文献[31];SD为标准差。
Note: a represents that the data come from Reference [31]; SD means standard deviation.
由Monte Carlo模拟与SQGs相结合,可以得到实际监测数据高于各个SQGs标准值的概率。表5列出了7种重金属超越4种SQGs标准值的具体概率。由表5可知,如果以ERL/ERM这对较为宽松的标准来看,4次监测7种重金属均以90%以上的概率处于低风险,因此以下分析均依据TEL/PEL这对较为严格的标准进行。
表5 7种重金属超越ERL和ERM的概率
Table 5 Probabilities of exceeding ERL and ERM for the seven heavy metals
重金属Heavy metal时间Sampling dateTELPELERLERM砷As2011.127.77E-0203.89E-0202012.0600002012.126.00E-0401.20E-0402013.060000镉Cd2011.128.43E-0305.70E-0402012.061.87E-0109.01E-0302012.127.52E-024.10E-041.69E-0202013.066.72E-022.60E-041.28E-020铬Cr2011.128.10E-0402.07E-0502012.061.63E-0305.42E-0502012.125.40E-0202.00E-0302013.065.24E-032.66E-055.80E-040铜Cu2011.125.00E-013.00E-045.45E-0202012.061.54E-013.00E-056.10E-0302012.125.34E-011.00E-039.72E-021.52E-052013.065.08E-014.00E-046.842E-022.00E-05汞Hg2011.121.20E-0405.74E-0502012.063.90E-0401.20E-0402012.126.00E-0501.97E-0502013.065.00E-0401.40E-040铅Pb2011.123.15E-021.00E-051.88E-0302012.061.71E-0105.35E-0302012.122.56E-013.71E-051.87E-0202013.063.58E-013.21E-052.00E-020锌Zn2011.125.31E-0301.56E-0302012.063.30E-0409.25E-0502012.122.63E-031.00E-056.40E-0402013.061.79E-021.21E-039.26E-033.30E-04
注:代号意义同表1。
Note: The meanings of symbols are the same as those in Table 1.
从不同监测时间来看,2011年12月Cu以49.91%的概率为中等风险,49.99%的概率为低风险,其他6种重金属均以90%以上的概率为低风险。2012年6月Cd、Cu和Pb处于低风险的概率分别为81.26%、84.60%和82.94%,处于中等风险的概率分别为18.74%、15.40%和17.06%,其他4种重金属均以99%以上的概率处于低风险。2012年12月Cu和Pb处于低风险的概率分别为46.58%和74.39%,处于中等风险的概率分别为53.32%和25.60%,其他5种重金属均以90%以上的概率处于低风险。2013年6月,Cu和Pb处于低风险的概率分别为49.24%和64.22%,处于中等风险的概率分别为50.72%和35.77%,其他5种重金属均以90%以上的概率处于低风险。
从不同重金属的角度来看,4次监测As、Cr、Hg和Zn均以90%以上的概率处于低风险。Cd除2012年6月以18.74%的概率处于中等风险外,其他3次均以90%以上的概率处于低风险。Cu 4次监测处于中等风险的概率分别为49.91%、15.40%、53.32%和50.72%。Pb 4次监测处于中等风险的概率分别为3.15%、17.06%、25.60%和35.77%。虽然2011年12月Cu的平均浓度为18.93 mg·kg-1,略大于其TEL(18.7 mg·kg-1),处于中等风险,但其仍有49.99%的概率处于低风险。2012年12月Cu的平均浓度为20.81 mg·kg-1,为其TEL(18.7 mg·kg-1)的1.11倍,处于中等风险,但其仍有46.58%的概率为低风险,可见应用Monte Carlo模拟与SQGs相结合能够得到更加详细的结果。
总体来讲,4次监测As、Cr、Hg和Zn均为低风险,Cd以80%以上的概率处于低风险,Pb以65%以上的概率处于低风险,Cu有3次以大约50%的概率处于中等风险。7种重金属对脱硫排污口附近海域的污染程度较小,其大小顺序为Cu>Pb>Cd>As≈Cr≈Hg≈Zn。
采用Monte Carlo模拟与PERI相结合的评价方法,可以得到各重金属生态危害指数的累积概率曲线(图2),以及其所处每个风险级别的概率(表6)。
表6 基于平均值估计和概率估计的单个重金属的生态风险
Table 6 Ecological risk of each heavy metal derived from both average estimation and probabilistic estimation
时间Sampling date重金属Heavy metal平均值估计Estimation of mean valueMonte Carlo模拟估计每个风险级别的概率Probability of each risk degree estimated by Monte Carlo simulation平均Eir值Mean value of Eir风险级别Risk degree轻微/%Low/%中等/%Moderate/%强/%Considerable/%很强/%High/%极强/%Extreme/%2011.12As6.44Low1000000Cd11.89Low99.690.31000Cr0.84Low1000000Cu1.57Low1000000Hg61.23Moderate5.8979.7214.200.180.01Pb3.42Low1000000Zn0.70Low10000002012.06As3.49Low1000000Cd23.76Low94.934.960.1100Cr0.80Low1000000Cu1.15Low1000000Hg75.76Moderate1.0160.3938.260.340Pb4.91Low1000000Zn0.46Low10000002012.12As5.43Low1000000Cd14.25Low96.093.290.530.070.02Cr1.25Low1000000Cu1.72Low1000000Hg57.58Moderate7.2985.916.770.030Pb4.93Low1000000Zn0.73Low10000002013.06As2.91Low1000000Cd14.41Low96.722.840.380.050.01Cr0.70Low1000000Cu1.55Low1000000Hg90.00High0.1534.4164.860.580Pb5.24Low1000000Zn0.52Low1000000
注:As、Cd、Cr、Cu、Hg、Pb和Zn分别指重金属砷、镉、铬、铜、汞、铅和锌;Low、Moderate和High分别表示风险级别为轻微、中等和很强。
Note: As, Cd, Cr, Cu, Hg, Pb and Zn means arsenic, cadmium, chromium, copper, mercury, lead and zinc, respectively; Low, Moderate, and High indicate that risk degree is low, moderate and very high, respectively.
图2 单个重金属潜在生态危害指数的累积概率曲线
Fig. 2 Cumulative probability curves of potential ecological risk factor for individual heavy metal
4次监测As、Cr、Cu、Pb和Zn的潜在生态危害均以100%的概率为轻微程度,4次监测Cd的潜在生态危害均以90%以上的概率为轻微程度。Hg前3次监测的平均潜在生态危害为中等程度,其处于中等程度的概率分别为79.72%、60.39%和85.91%,2013年6月Hg以64.68%的概率为强危害程度。
2011年12月、2012年6月、2012年12月和2013年6月的综合生态危害指数(RI)分别为86.10、110.32、85.89和115.33,均<150(轻微程度的上限值),属于轻微危害程度。Monte Carlo模拟结果显示,4次监测沉积物中重金属的综合生态危害以90%以上的概率为轻微程度(图3)。
图3 7种重金属综合潜在生态危害指数(RI)的累积概率曲线
Fig. 3 Cumulative probability curves of comprehensive potential ecological risk index (RI) of seven heavy metals
总体来讲,4次监测沉积物中重金属污染程度均为轻微。SQGs显示Cu、Cd和Pb污染程度相对较大。PERI方法结果显示Hg的污染程度最大,其中夏季(2012.06和2013.06)重金属的污染程度要略高于冬季(2011.12和2012.12)。
本研究海域表层沉积物重金属平均浓度与胶州湾以及其他受人类活动影响较大海域的对比如表7所示。对比分析发现,本研究4次监测表层沉积物中Cd、Cr、Cu、Pb和Zn的平均浓度均高于2004年5月胶州湾的平均浓度[36],表明脱硫海水的长期大量排放确实引起了附近海域沉积物中重金属的累积。本研究Cd、Hg和Pb的平均浓度高于2009年胶州湾跨海大桥附近沉积物中的平均浓度[37],但大部分重金属含量低于胶州湾东部沉积物中的含量[38]。类似的,李玉等[39]运用平均富集因子发现胶州湾表层沉积物中重金属主要分布在胶州湾东部李村河口、娄山河口以及附近海域。魏璟弢等[40]应用连续提取法研究了2008年11月和2010年11月青岛近海和胶州湾内沉积物的富集状况,发现Pb和Cr存在明显富集,胶州湾东部重金属污染程度最为明显。本研究As、Cd、Cu、Hg、Pb和Zn的浓度比五里河河口附近海域(葫芦岛锌厂排污口邻近海域)低一个数量级[41]。与我国高雄湾[15]、香港海岸带[40]、泉州湾[42]以及国外一些其他海湾(表7)相比,本研究沉积物中重金属含量均低于上述区域。总体上,该燃煤电厂排污口附近海域沉积物中重金属的浓度较低。
表7 电厂排污口附近海域与其他受影响海域中表层沉积物重金属的平均浓度比较
Table 7 Mean concentrations of heavy metals found in sediment near the outlet of the power plant compared to the reported average concentrations for other impacted coastal systems (mg·kg-1)
地区Area砷As镉Cd铬Cr铜Cu汞Hg铅Pb锌Zn参考文献Reference中国胶州湾Jiaozhou Bay, China3.540.3626.2218.040.0424.3747.21本研究This study中国胶州湾Jiaozhou Bay, China11.150.05922.1411.99-10.3131.1[36]中国胶州湾Jiaozhou Bay, China-0.055-19.050.03112.751.14[37]中国胶州湾Jiaozhou Bay, China10.900.5577.4036.23-52.94161.32[38]中国五里河口Wulihe Estuary, China88.254.11-116.600.56104.651 008.75[41]中国泉州湾Quanzhou Bay, China21.70.5982.071.40.4067.70179.6[42]中国香港近岸Hong Kong Coasts, China-0.3348.93118.680.1953.56147.73[47]中国台湾高雄港Kaohsiung Harbor, Taiwan, China-0.93154.15117.600.4655.23318.33[48]意大利那不勒斯湾Gulf of Naples, Italy2.00.572827.20.70221602[49]巴西里贝拉湾Ribeira Bay, Brazil-0.20711824.6-22.9109[50]巴西塞佩提巴湾Sepetiba Bay, Brazil-3.226631.9-40567[50]
目前燃煤电厂脱硫海水重金属的研究主要集中于具有挥发性质的Hg在烟气和水环境介质中的迁移转化以及其不同形态在水环境中的分布特征[1]。Liang等[43]研究发现燃煤电厂排污口附近海域表层沉积物中Hg的浓度范围为0.055~0.201 mg·kg-1,平均值为0.126 mg·kg-1,为本研究Hg浓度的3倍左右。本研究采用PERI方法过程中由于选用了该地区Hg的较低背景值0.024 mg·kg-1[31,44],导致其潜在生态危害指数范围为57.58~90.00,污染程度远大于其他重金属。此外,PERI方法不仅考虑了重金属的环境背景值,还考虑了其毒性系数。Hg由于具有高毒性,尤其是通过生物体产生的有机汞,其毒性系数采用了所有重金属元素中最高的40。Li等[35]采用地累积指数法和PERI法对2009年7月胶州湾表层沉积物8种重金属的污染和潜在生态危害进行了评价,发现重金属总体污染程度较轻,但Hg污染较为严重,与本研究结果一致。此外,由于夏季(2012.06和2013.06)Hg的监测浓度要略高于冬季(2011.12和2012.12),进而导致夏季重金属的综合生态危害指数(RI)也略高于冬季。本研究4次监测重金属的RI均<150(轻微程度的上限值),属于轻微危害程度。类似的,郭军辉等[38]采用潜在生态危害指数法对2009年胶州湾东岸表层沉积物重金属的生态危害进行了评价,发现Cd对沉积物生态环境的潜在危害为中等程度,As、Cr、Cu、Pb和Zn的危害程度均为轻微,多个重金属的RI值均<150,属于轻微危害程度。
目前大部分针对沉积物重金属生态风险评价的方法,包括SQGs和PERI,都采用监测数据的平均值、保守值(90%或95%分位数)或最大值来获得风险的一个平均或保守的估计[25-26]。这种将监测数据压缩为一个单点值的方法会导致生态风险评价中信息的丢失,产生不确定性,从而造成对风险的高估或者低估。本研究中2012年12月Cu的平均浓度为20.81 mg·kg-1,为其TEL(18.7 mg·kg-1)的1.11倍,处于中等风险,但结合Monte Carlo模拟发现其仍有46.58%的概率为低风险。类似的,Sawe等[45]发现坦桑尼亚的瓦米河口中As、Cd、Cr、Pb和Zn的浓度以99%的概率低于SQGs的效应范围低值,风险可忽略,但Cu却有一定的概率对该河口造成风险。此外,本研究2012年6月Hg的平均值为75.76(为中等程度上限值80的0.947倍),生态危害程度为中等,但由Monte Carlo模拟的结果可得其仍有38.26%的概率处于强危害程度。Qu等[34]指出虽然香江和天池湖的平均RI值略高于高危害程度(high risk)的上限值600(分别为其1.08倍和1.17倍),处于高危害程度,但Monte Carlo模拟结果显示其仍有43.3%和47.1%的概率处于较高危害程度(considerable risk)。Li等[35]指出尽管莱州湾西部Hg的平均
为20.24(为轻微程度上限值40的0.51倍),但Monte Carlo模拟显示仍有一定的概率(0.03%)处于中等危害程度。可见只用平均值进行评价,可能会高估或者低估重金属的真实风险。采用Monte Carlo模拟与传统评价方法相结合能够有效降低风险评估的不确定性,为污染水域沉积物重金属风险缓解提供更多的决策支持。
目前对于像燃煤电厂这样的点污染源的暴露评价主要是基于电厂建设前的模型预测研究,缺少电厂运行后的长期跟踪监测[2]。并且模型涉及大量的污染物特有参数(如分子量、蒸气压、吸附系数等)和环境场景参数,其结果往往具有一定的不确定性[46]。本研究将脱硫海水工艺运行10年后,连续3年的跟踪监测数据作为环境暴露数据能更好地表征重金属排放对排污口附近海域的长期生态风险。SQGs和PERI方法均显示4次监测沉积物中重金属污染程度为轻微,表明燃煤电厂脱硫海水的长期排放并未造成排污口附近海域表层沉积物重金属的污染。除重金属外,脱硫海水中的其他特征性污染物pH(H+)、SO42-、温度等也可能对排污口附近海域造成不利影响[2]。为研究脱硫海水排放对海洋生态环境的综合影响,本研究还根据调查了排污口附近海域的生物群落结构(本研究未列出),共发现浮游植物4门93种,其中硅藻78种;浮游动物58种,其中节肢动物类和浮游幼虫类分别有27种和17种;底栖动物8门128种,其中环节动物门69种。依据Shannon-Wiener多样性指数和Pielou均匀度指数[31],排污口附近海域各类群生物多样性较高,群落结构较稳定,表明脱硫海水排放未对生物群落结构和生物多样性造成明显不利影响。
[1] Liu X Y, Sun L M, Yuan D X, et al. Mercury distribution in seawater discharged from a coal-fired power plant equipped with a seawater flue gas desulfurization system [J]. Environmental Science and Pollution Research International, 2011, 18(8): 1324-1332
[2] Córdoba P. Status of flue gas desulphurisation (FGD) systems from coal-fired power plants: Overview of the physic-chemical control processes of wet limestone FGDs [J]. Fuel, 2015, 144: 274-286
[3] Clarke L B. The fate of trace elements during coal combustion and gasification: An overview [J]. Fuel, 1993, 72(6): 731-736
[4] Aunela-Tapola L, Hatanpää E, Hoffren H, et al. A study of trace element behaviour in two modern coal-fired power plants: Ⅱ. Trace element balances in two plants equipped with semi-dry flue gas desulphurisation facilities [J]. Fuel Processing Technology, 1998, 55: 13-34
[5] Radojevi M. The use of seawater for flue gas desulphurization [J]. Environmental Technology, 1989, 10(1): 71-76
[6] Barrero F V, Ollero P, Gutiérrez Ortiz F J, et al. Catalytic seawater flue gas desulfurization process: An experimental pilot plant study [J]. Environmental Science &Technology, 2007, 41(20): 7114-7119
[7] Zhuang W, Gao X L. Integrated assessment of heavy metal pollution in the surface sediments of the Laizhou Bay and the coastal waters of the Zhangzi Island, China: Comparison among typical marine sediment quality indices [J]. PLoS One, 2014, 9(4): e94145
[8] Ranjbar Jafarabadi A, Riyahi Bakhtiyari A, Shadmehri Toosi A, et al. Spatial distribution, ecological and health risk assessment of heavy metals in marine surface sediments and coastal seawaters of fringing coral reefs of the Persian Gulf, Iran [J]. Chemosphere, 2017, 185: 1090-1111
[9] Carballeira C, De Orte M R, Viana I G, et al. Assessing the toxicity of chemical compounds associated with land-based marine fish farms: The sea urchin embryo bioassay with Paracentrotus lividus and Arbacia lixula [J]. Archives of Environmental Contamination and Toxicology, 2012, 63(2): 249-261
[10] Garr A L. Toxicity and availability of copper and zinc to queen conch: Implications for larval recruitment in the Florida Keys [D]. Boca Raton, FL: Florida Atlantic University, 2012: 31-35
[11] Debelius B, Forja J M, DelValls A, et al. Toxicity and bioaccumulation of copper and lead in five marine microalgae [J]. Ecotoxicology and Environmental Safety, 2009, 72(5): 1503-1513
[12] Al-Subiai S N, Moody A J, Mustafa S A, et al. A multiple biomarker approach to investigate the effects of copper on the marine bivalve mollusc, Mytilus edulis [J]. Ecotoxicology and Environmental Safety, 2011, 74(7): 1913-1920
[13] 郭娟, 袁东星, 陈进生, 等. 燃煤电厂海水脱硫工艺的排水对海域环境的影响[J]. 环境工程学报, 2008, 2(5): 707-711
Guo J, Yuan D X, Chen J S, et al. Impact of discharge water from seawater flue gas desulfurization system of coal-fired power plant on the environment of surrounding sea area [J]. Chinese Journal of Environmental Engineering, 2008, 2(5): 707-711 (in Chinese)
[14] lvarez-Ayuso E, Querol X, Tomás A. Environmental impact of a coal combustion-desulphurisation plant: Abatement capacity of desulphurisation process and environmental characterisation of combustion by-products [J]. Chemosphere, 2006, 65(11): 2009-2017
[15] Chen Y W, Liu G J, Wang L, et al. Occurrence and fate of some trace elements during pyrolysis of Yima coal, China [J]. Energy &Fuels, 2008, 22(6): 3877-3882
[16] Liu B L, Hu K, Jiang Z L, et al. Distribution and enrichment of heavy metals in a sediment core from the Pearl River Estuary [J]. Environmental Earth Sciences, 2011, 62(2): 265-275
[17] Sun L M, Feng L F, Yuan D X, et al. The extent of the influence and flux estimation of volatile mercury from the aeration pool in a typical coal-fired power plant equipped with a seawater flue gas desulfurization system [J]. The Science of the Total Environment, 2013, 444: 559-564
[18] Gaur V K, Gupta S K, Pandey S D, et al. Distribution of heavy metals in sediment and water of River Gomti [J]. Environmental Monitoring and Assessment, 2005, 102(1-3): 419-433
[19] Kowalska J B, Mazurek R, Gąsiorek M, et al. Pollution indices as useful tools for the comprehensive evaluation of the degree of soil contamination-A review [J]. Environmental Geochemistry and Health, 2018, 40(6): 2395-2420
[20] Long E R, MacDonald D D, Smith S L, et al. Incidence of adverse biological effects within ranges of chemical concentrations in marine and estuarine sediments [J]. Environmental Management, 1995, 19(1): 81-97
[21] MacDonald D D, Carr R S, Calder F D, et al. Development and evaluation of sediment quality guidelines for Florida coastal waters [J]. Ecotoxicology, 1996, 5(4): 253-278
[22] Hakanson L. An ecological risk index for aquatic pollution control. A sedimentological approach [J]. Water Research, 1980, 14(8): 975-1001
[23] Li D P, Yu R D, Chen J, et al. Ecological risk of heavy metals in lake sediments of China: A national-scale integrated analysis [J]. Journal of Cleaner Production, 2022, 334: 130206
[24] 孔令晖, 罗明科, 余辉, 等. 长江中游湖南段湖泊重金属时空累积因素研究[J]. 中国环境科学, 2023, 43(4): 1812-1823
Kong L H, Luo M K, Yu H, et al. Spatial-temporal accumulation factors of heavy metals in Hunan Lakes along the middle reach of the Yangtze River [J]. China Environmental Science, 2023, 43(4): 1812-1823 (in Chinese)
[25] Su L Y, Liu J L, Christensen P. Spatial distribution and ecological risk assessment of metals in sediments of Baiyangdian wetland ecosystem [J]. Ecotoxicology, 2011, 20(5): 1107-1116
[26] Pejman A, Nabi Bidhendi G, Ardestani M, et al. A new index for assessing heavy metals contamination in sediments: A case study [J]. Ecological Indicators, 2015, 58: 365-373
[27] Bai J H, Cui B S, Chen B, et al. Spatial distribution and ecological risk assessment of heavy metals in surface sediments from a typical plateau lake wetland, China [J]. Ecological Modelling, 2011, 222(2): 301-306
[28] National Research Council (NRC). Science and judgment in risk assessment [R]. Washington DC: National Academy Press, 1994
[29] United States Environmental Protection Agency (US EPA). Guiding principles for Monte Carlo analysis [R]. Washington DC: US EPA, 1997
[30] Taylor M R G, Heaton R, Baty R. The impact of flue-gas desulphurization on the water environment [J]. Water and Environment Journal, 1989, 3(3): 227-234
[31] 王云鹏, 纪良, 时国梁, 等. 脱硫海水排放对附近海域的生态环境影响初步研究[J]. 中国海洋大学学报(自然科学版), 2014, 44(5): 61-68
Wang Y P, Ji L, Shi G L, et al. A preliminary study of the impact of discharge water from seawater flue gas desulphurization on the ecological environment of surrounding sea area [J]. Periodical of Ocean University of China, 2014, 44(5): 61-68 (in Chinese)
[32] McCready S, Birch G F, Long E R. Metallic and organic contaminants in sediments of Sydney Harbour, Australia and vicinity—A chemical dataset for evaluating sediment quality guidelines [J]. Environment International, 2006, 32(4): 455-465
[33] Christophoridis C, Dedepsidis D, Fytianos K. Occurrence and distribution of selected heavy metals in the surface sediments of Thermaikos Gulf, N. Greece. Assessment using pollution indicators [J]. Journal of Hazardous Materials, 2009, 168(2-3): 1082-1091
[34] Qu C S, Li B, Wu H S, et al. Probabilistic ecological risk assessment of heavy metals in sediments from China’s major aquatic bodies [J]. Stochastic Environmental Research and Risk Assessment, 2016, 30(1): 271-282
[35] Li X, Chi W Q, Tian H, et al. Probabilistic ecological risk assessment of heavy metals in western Laizhou Bay, Shandong Province, China [J]. PLoS One, 2019, 14(3): e0213011
[36] Wang X C, Feng H, Ma H Q. Assessment of metal contamination in surface sediments of Jiaozhou Bay, Qingdao, China [J]. CLEAN - Soil, Air, Water, 2007, 35(1): 62-70
[37] 张珂, 王朝晖, 冯杰, 等. 胶州湾表层沉积物重金属分布特征及污染评价[J]. 分析测试学报, 2011, 30(12): 1406-1411
Zhang K, Wang Z H, Feng J, et al. Distribution and pollution assessment of heavy metals in surface sediments of Jiaozhou Bay [J]. Journal of Instrumental Analysis, 2011, 30(12): 1406-1411 (in Chinese)
[38] 郭军辉, 殷月芬, 陈发荣, 等. 胶州湾表层沉积物重金属污染分布特征及其生态风险评价[J]. 环境污染与防治, 2012, 34(3): 13-21
Guo J H, Yin Y F, Chen F R, et al. The distribution characteristics of heavy metals in sediments of Jiaozhou Bay and its potential ecological risk evaluation [J]. Environmental Pollution &Control, 2012, 34(3): 13-21 (in Chinese)
[39] 李玉, 俞志明, 曹西华, 等. 重金属在胶州湾表层沉积物中的分布与富集[J]. 海洋与湖沼, 2005, 36(6): 580-589
Li Y, Yu Z M, Cao X H, et al. Distribution and enrichment of heavy metals in surface sediments of Jiaozhou Bay [J]. Oceanologia et Limnologia Sinica, 2005, 36(6): 580-589 (in Chinese)
[40] 魏璟弢, 张焕玲, 李铁, 等. 胶州湾及青岛近海表层沉积物重金属赋存形态研究[J]. 中国海洋大学学报(自然科学版), 2012, 42(S1): 157-164
Wei J T, Zhang H L, Li T, et al. Study of the geochemical forms of heavy metals in surface sediments of Jiaozhou Bay and Qingdao inshore areas [J]. Periodical of Ocean University of China, 2012, 42(S1): 157-164 (in Chinese)
[41] Feng Y L, Chen Y Z, Wang J, et al. Ecological risk assessment of heavy metals in surface seawater and sediment near the outlet of a zinc factory in Huludao City, Liaoning Province, China [J]. Chinese Journal of Oceanology and Limnology, 2016, 34(6): 1320-1331
[42] Yu R L, Yuan X, Zhao Y H, et al. Heavy metal pollution in intertidal sediments from Quanzhou Bay, China [J]. Journal of Environmental Sciences (China), 2008, 20(6): 664-669
[43] Liang Y, Liu X Y, Yuan D X, et al. Mercury species in seawater and sediment of Xiamen western sea area adjacent to a coal-fired power plant [J]. Water Environment Research: A Research Publication of the Water Environment Federation, 2010, 82(4): 335-341
[44] 迟清华, 鄢明才. 应用地球化学元素丰度数据手册[M]. 北京: 地质出版社, 2007: 67-75
[45] Sawe S F, Shilla D A, Machiwa J F. Assessment of ecological risk of heavy metals using probabilistic risk assessment model (AQUARISK) in surface sediments from Wami Estuary, Tanzania [J]. BioMed Research International, 2021, 2021: 6635903
[46] 周林军, 冯洁, 刘济宁, 等. 污水处理厂中化学品暴露预测模型研究进展[J]. 环境科学与技术, 2015, 38(8): 68-74
Zhou L J, Feng J, Liu J N, et al. Models of exposure prediction for chemicals in the sewage treatment plant: A review [J]. Environmental Science &Technology, 2015, 38(8): 68-74 (in Chinese)
[47] Chen C W, Kao C M, Chen C F, et al. Distribution and accumulation of heavy metals in the sediments of Kaohsiung Harbor, Taiwan [J]. Chemosphere, 2007, 66(8): 1431-1440
[48] Zhou F, Guo H C, Hao Z J. Spatial distribution of heavy metals in Hong Kong’s marine sediments and their human impacts: A GIS-based chemometric approach [J]. Marine Pollution Bulletin, 2007, 54(9): 1372-1384
[49] Romano E, Ausili A, Zharova N, et al. Marine sediment contamination of an industrial site at Port of Bagnoli, Gulf of Naples, Southern Italy [J]. Marine Pollution Bulletin, 2004, 49(5-6): 487-495
[50] de Carvalho Gomes F, Godoy J M, Godoy M L, et al. Metal concentrations, fluxes, inventories and chronologies in sediments from Sepetiba and Ribeira Bays: A comparative study [J]. Marine Pollution Bulletin, 2009, 59(4-7): 123-133