水产种质资源保护区是指为保护水产种质资源及其生存环境，在具有较高经济价值和遗传育种价值的水产种质资源的主要生长繁育区域(即产卵场、索饵场、越冬场以及洄游通道等)，依法划定并予以特殊保护和管理的水域、滩涂及其毗邻的岛礁、陆域等(农业部, 2011)。近年来，工程建设、环境污染等大量占用、破坏重要渔业生物栖息地，严重影响到渔业可持续发展和国家生态文明建设。自2007年起我国积极推进水产种质资源保护区的建立，目前，全国已划定水产种质资源保护区535处(农业农村部, 2018)。建立水产种质资源保护区是协同推进经济发展和资源环境保护的有效举措，对维护渔业生物资源和水域生态健康具有重要意义。

为了掌握保护区生态环境质量状况，同时，为保护区管理部门制定和调整管理措施提供科学依据，使水产种质资源及其栖息环境得到有效的保护，对保护区生态环境质量状况进行监测和评价是十分重要的。目前，已对国家级水产种质资源保护区开展了一系列生态环境调查、监测和综合质量评价。夏斌等(2010)、崔毅等(2013)和滕瑶等(2018)对靖海湾松江鲈鱼(Trachidermus fasciatus)国家级水产种质资源保护区开展了多年连续监测和环境质量综合评价。结果表明，该海域长期处于富营养化状态且有机污染严重。袁健美等(2013)报道了吕泗渔场小黄鱼(Larimichthys polyactis)、银鲳(Pampus argenteus)国家级水产种质资源保护区春、夏季环境因子水平，监测结果表明，水体中悬浮物、无机氮(DIN)和磷酸盐(PO₄-P)超标严重，与历史资料相比，该海区富营养化程度加剧。徐衡等(2016)监测和评价了乐清湾泥蚶(Tegillarca granosa)国家级水产种质资源保护区水质及沉积物质量状况，表明该保护区综合水质为良，优于乐清湾北部海水养殖区的水质状况，但沉积物中Cu、Zn、Cr和粪大肠菌群超标。

海州湾大竹蛏(Solen grandis)国家级水产种质资源保护区，由农业部于2007年12月12日批准建立，主要保护对象为海珍品——大竹蛏。夏斌等(2014)研究表明，2008~2013年夏季海州湾水产种质资源保护区海水环境质量均符合国家Ⅰ类海水水质标准，属于尚清洁等级，但近年来，对海州湾保护区海水环境质量状况的研究尚未见报道。本研究根据2018年5月和8月对海州湾大竹蛏国家级水产种质资源保护区的生态环境监测数据，运用综合质量指数法对海州湾保护区的海水环境质量进行综合评价，并从水体营养状况和有机污染状况等不同角度对该海域的海水水质状况进行分析，旨在及时掌握海州湾渔业水域的生态环境质量状况，为海州湾水产种质资源保护提供科学依据。

海州湾大竹蛏国家级水产种质资源保护区位于山东省日照市东港区邻近的黄海海域。该保护区总面积为42.88 hm²，其中，核心区面积为15.20 hm²，实验区面积为27.68 hm²。核心区特别保护期为每年4~7月。该海域海水水质清澈，海底底质为泥砂质，生态环境良好，是大竹蛏天然资源集中生长繁殖的重要场所。

在海州湾大竹蛏国家级水产种质资源保护区内共布设5个调查站位进行海洋环境调查。其中，1~4号站位于核心区与实验区交界处，5号站位于核心区(图 1)。采样时间分别为2018年5月(春季)和8月(夏季)，分别为核心区特别保护期中和保护期后。

图 1(Fig. 1)

图 1 海州湾保护区调查站位 Fig.1 Location of sampling sites in Haizhou Bay

海水样品使用Niskin采水器进行采集，采样层次为表层。检测项目包括水温、盐度、pH、溶解氧(DO)、化学需氧量(COD)、DIN、PO₄-P、Cu、Zn、Pb、Cd、Hg、As以及石油类。所有样品的采集、现场处理及检测分析均依据《海洋监测规范》(GB17378.4-2007)中所规定的方法执行。其中，温度、盐度及pH使用YSI多参数分析仪(556MPS, 美国)现场测定；DO采用碘量法测定；COD采用碱性高锰酸钾法测定；亚硝酸盐采用萘乙二胺分光光度法测定；硝酸盐采用锌—镉还原法测定；铵盐采用次溴酸盐氧化法测定；PO₄-P采用磷钼蓝分光光度法测定；Cu、Zn、Pb、Cd采用原子吸收分光光度法测定；Hg、As采用原子荧光分光光度法测定；石油类采用荧光分光光度法测定。

综合质量指数法能很好地反映某一海域环境综合质量(贾晓平等, 2003)，评价公式：

式中，Q为综合质量指数，n为环境因子项数，P_i为环境因子的污染指数，C_i为环境因子实测值，C_s为环境因子评价标准值。

由于海水pH的评价标准是一个范围，其污染指数计算公式(崔毅等, 2013)如下：

式中，S_{i, pH}为第i站pH的质量指数，pH_i为第i站pH实测值，pH_su为pH评价标准上限，pH_sd为pH评价标准下限。

海水中DO含量随污染程度的增大而减少，其污染指数计算公式(崔毅等, 2013)如下：

式中，S_{i, DO}为第i站DO的质量指数，DO_f为相同温度相同盐度条件下，DO的饱和浓度值(mg/L)，DO_i为DO的实测浓度(mg/L)，DO_s为DO的评价标准值(mg/L)。

本研究调查海域位于海州湾大竹蛏国家级水产种质资源保护区内，采用《海水水质标准》(GB 3097-1997)Ⅰ类标准作为评价标准值，综合质量指数与环境分级的关系见表 1(夏斌等, 2014)。

表 1(Tab.1)

表 1 综合质量指数与环境分级关系 Tab.1 The relationship between comprehensive quality index and environmental classification 分级Classification 水质质量综合指数Comprehensive quality index (Q) 清洁 Clean Q < 0.3 尚清洁 Slightly clean 0.3≤Q < 0.7 允许 Permited 0.7≤Q < 1.0 轻污染 Lightly polluted 1.0 < Q < 2.0 污染 Polluted 2.0≤Q < 3.0 重污染Seriously polluted 3.0≤Q < 5.0 恶性污染 More seriously polluted Q≥5.0

表 1 综合质量指数与环境分级关系 Tab.1 The relationship between comprehensive quality index and environmental classification

海水营养水平采用富营养化指数评价，该方法是目前应用较为广泛的一种方法，其计算公式(邹景忠等, 1983)为：

式中，E为富营养化指数，COD_i为第i站化学需氧量(mg/L)，DIN_i和PO₄-P_i分别为第i站DIN和PO₄-P的实测浓度(μg/L)。当0 < E < 0.5时，海水为贫营养；0.5≤E < 1时，海水为中营养；E≥1时，海水呈富营养化，且E值越高，海水富营养化程度越严重。

采用有机污染指数法评价海水有机污染程度，其计算公式(蒋国昌等, 1987)为：

式中，A为有机污染指数，COD_i和DO_i分别为第i站COD和DO的实测值(mg/L)，DIN_i和PO₄-P_i分别为第i站DIN和PO₄-P的实测值(μg/L)，COD_s、DIN_s、PO₄-P_s和DO_s分别为对应因子的Ⅰ类《海水水质标准》(GB3097- 1997)评价值。依据有机污染分级标准(表 2)(蒋国昌等, 1987)评价以上公式计算结果。

表 2(Tab.2)

表 2 有机污染指数与海水有机污染等级关系 Tab.2 The relationship between organic pollution index and organic pollution classification in seawater 分级Classification 有机污染指数Organic pollution index (A) 良好 Good A < 0 较好 Preferable 0 < A < 1 开始受到污染 Begin to be polluted 1 < A < 2 轻度污染 Lightly polluted 2 < A < 3 中度污染 Moderately polluted 3 < A < 4 严重污染 Seriously polluted A > 4

表 2 有机污染指数与海水有机污染等级关系 Tab.2 The relationship between organic pollution index and organic pollution classification in seawater

2018年春季和夏季，海州湾大竹蛏国家级水产种质资源保护区海水中相关环境因子的含量见表 3。从表 3可以看出，2018年春季和夏季表层海水盐度的变化范围分别为30.02~30.23和28.06~28.11，平均值分别为30.13和28.08。根据我国《海水水质标准》(GB 3097-1997)有关规定，海州湾大竹蛏国家级水产种质资源保护区的海水水质应符合第Ⅰ类海水水质标准。实际调查数据显示，2018年春季2号站和夏季5号站的DIN含量略超出Ⅰ类海水水质标准限值(分别超出1.51%和0.16%)，其他监测因子均符合Ⅰ类海水水质标准，表明海州湾大竹蛏国家级水产种质资源保护区的海水水质较好。该结果与2008~2013年夏季海州湾保护区的海水环境质量评价结果较为相似(夏斌等, 2014)，其中，DIN含量有所升高，PO₄-P、Zn、Cd和As含量有所下降。与靖海湾松江鲈鱼国家级水产种质资源保护区近年的海水环境质量(滕瑶等, 2018)相比较，海州湾保护区的pH、DO水平较高，营养盐、重金属含量水平较低。这可能有2个原因，一是相比于靖海湾保护区，海州湾保护区离陆地较远，受陆源排污影响较小；二是靖海湾为半封闭海湾，而海州湾为半开阔海湾，后者的海水交换能力和自净能力更强。

表 3(Tab.3)

表 3 2018年海州湾大竹蛭国家级水产种质资源保护区主要海水环境因子含量 Tab.3 Content of the water quality factors in Haizhou Bay in 2018 时间Time 指标Index 湿度Temperature (℃) 盐度Salinity pH 溶解氧DO(mg/L) 化学需氧量COD(mg/L) 溶解无机氮DIN(ug/L) 活性磷酸盐PO4-P(ug/L) 铜Cu(ug/L) 锌Zn(ug/L) 铅Pb(ug/L) 镉Cd(ug/L) 汞Hg(ug/L) 砷As(ug/L) 石油类Petroleum(ug/L) 5月 May 范围Range 19.29~19.99 30.02~30.23 7.98~8.10 8.12~8.27 0.74~1.07 111.03~203.02 1.81~7.41 1.65~2.14 4.83~5.68 0.26~0.35 0.18~0.26 0.024~0.039 0.42~0.82 10.27~36.38 平均值Average 19.55 30.13 8.06 8.20 0.89 163.67 3.49 1.87 5.19 0.30 0.20 0.030 0.71 17.38 8月 August 范围Range 28.03~28.23 28.06~28.11 8.02~8.11 6.17~6.42 0.91~1.15 193.45~200.32 1.49~4.46 1.33~2.34 5.34~6.04 0.27~0.37 0.16~0.28 0.014~0.022 0.92~1.40 9.83~18.06 平均值Average 28.14 28.08 8.07 6.27 1.05 196.34 2.67 1.73 5.82 0.32 0.22 0.019 1.20 12.81 平均值Average 23.84 29.11 8.06 7.23 0.97 180.00 3.08 1.80 5.51 0.31 0.21 0.024 0.96 15.10

表 3 2018年海州湾大竹蛭国家级水产种质资源保护区主要海水环境因子含量 Tab.3 Content of the water quality factors in Haizhou Bay in 2018

基于海州湾大竹蛏国家级水产种质资源保护区环境因子指标(表 3)，运用综合质量指数法对该保护区的生态环境质量进行综合评价见图 2。从图 2可以看出，2018年春季海州湾调查海域海水综合质量指数范围是0.31~0.39，平均值为0.34；夏季海水综合质量指数范围是0.34~0.38，平均值为0.36。从调查时间分析，海州湾保护区海水综合质量指数春季 < 夏季。从调查站位分析，2018年春季和夏季各站位间海水综合质量指数相差较小。总体而言，2018年海州湾大竹蛏国家级水产种质资源保护区综合质量指数均处于0.31~0.39之间，平均值为0.35，根据综合质量指数与环境分级关系(表 1)，判定调查海域海水质量等级为尚清洁。与历史调查数据比较，2008~2013年海州湾大竹蛏国家级水产种质资源保护区综合质量指数为0.33~0.68，同为尚清洁等级(夏斌等, 2014)，表明近年来海州湾保护区环境质量保持较好。贾晓平等(2003)研究表明，综合评价得出广西北部湾渔场的海水环境综合质量指数为0.45，海州湾保护区综合质量指数与北部湾渔场水平相当，海水环境综合质量较好，有利于保护大竹蛏的生长和繁殖。

图 1(Fig. 1)

图 2 2018年海州湾调查站位的综合质量指数 Fig.2 The comprehensive quality index of each station in Haizhou Bay in 2018

海水中的N、P等营养盐是浮游植物生长、繁殖所必需的物质，也是限制初级生产力的关键因子。2018年春季，海州湾调查海域N/P值为27.40~88.18，平均值为58.97；夏季海州湾调查海域N/P值为43.80~132.05，平均值为87.97。2018年海州湾大竹蛏国家级水产种质资源保护区N/P值均大于Redfield系数(Redfield, 1958)，平均值为73.47，表明该海域营养盐结构为磷限制。调查数据显示，该海域海水中PO₄-P含量(表 3)远低于浮游植物正常生长所必需的含量(17.98 μg/L)(郭卫东等, 1998)，可能成为影响浮游植物生长繁殖的限制因子。夏斌等(2014)研究表明，2008~2013年海州湾保护区的N/P值为16.44~91.62，平均值为51.67。与本研究2018年对海州湾保护区的调查结果相似，海州湾保护区N/P值远大于Redfield比值且PO₄-P含量偏低，可能会影响浮游植物种群结构特征，进而干扰整个生态系统(赵玉庭等, 2016)，对大竹蛏的生长产生潜在的不利影响。总体来看，海州湾大竹蛏保护区近10年来营养盐结构失衡状况日趋严重，这与我国近岸海域N、P营养盐结构失衡严重结论相一致(孙伟等, 2016)。

2018年海州湾保护区海水中DIN出现超标，而PO₄-P含量缺乏，进一步使用富营养化指数对该海域的营养状况进行综合评价。结果表明，2018年春季海州湾调查海域富营养化指数的变化范围为0.05~0.27，平均值为0.12；夏季海州湾调查海域富营养化指数为0.07~0.22，平均值为0.12。总体来看，2018年所有调查站位富营养化指数均 > 0.5，海水为贫营养状况。贾晓平等(2003)评价了广西北部湾渔场海水富营养化指数范围为0.003~0.083，平均值为0.022，海水处于贫营养状态，与本研究调查结果相似。胡益峰等(2015)研究了浙江嵊泗马鞍列岛国家级海洋特别保护区的富营养化指数范围为1.29~4.65，该海区在不同季节均呈富营养化，且富营养化程度较严重。海州湾大竹蛏保护区海水营养水平与天然渔场相似，水质较好，主要得益于海州湾大竹蛏国家级水产种质资源保护区的建立，更好地保护了该海域的水质，为大竹蛏提供了良好的栖息环境。

根据有机污染评价方法计算得出海州湾保护区的有机污染指数。2018年春季海州湾调查海域有机污染指数为—0.24~0.57，平均值为0.13；夏季海州湾调查海域有机污染指数为0.54~0.81，平均值为0.64。总体而言，2018年海州湾大竹蛏国家级水产种质资源保护区有机污染指数的平均值为0.39。从调查时间分析，海州湾保护区有机污染指数春季 < 夏季。从调查站位分析，春季1号站和5号站A < 0，属于良好等级，春季其他站位和夏季各站位0 < A < 1，属于较好等级。有机污染指数表明，海州湾保护区有机污染程度较轻。对比富营养化指数和有机污染指数2种评价方法。结果表明，海州湾保护区海水水质较好，富营养化程度和有机污染程度轻微。与其他海域相比，2012~2013年浙江嵊泗马鞍列岛保护区有机污染指数为3.31~4.53，属中度及重度污染(胡益峰等, 2015)。2013年莱州湾海域(包括单环刺螠近江牡蛎国家级水产种质资源保护区)有机污染指数平均值为1.61，开始受到有机污染(赵玉庭等, 2016)。2014~2016年靖海湾松江鲈鱼保护区有机污染指数的变化范围为4.66~9.23，处于严重污染状态(滕瑶等, 2018)。

2018年海州湾大竹蛏保护区海水综合质量指数范围为0.31~0.39，海水质量等级为尚清洁。核心区和实验区内各调查站位间海水综合质量指数变化较小。与2008~2013年调查数据相比，海水综合质量指数有所下降，表明近年来海州湾保护区环境质量保持较好。

从营养结构分析，2018年海州湾保护区营养盐结构为磷限制，海水中较低的PO₄-P含量成为影响浮游植物生长繁殖的限制因子。根据富营养化评价，2018年海州湾保护区内各调查站位富营养化指数均小于0.5，海水为贫营养状况。根据有机污染评价，2018年海州湾保护区有机污染指数为—0.24~0.81，表明该水域有机污染程度较轻。