唐山市曹妃甸区位于河北省东北部，唐山市东南部55 km，北靠滦南，东接乐亭，西邻丰南，南临渤海，地理坐标为：39°07′43″N~39°27′23″N，118°12′12″E~118°43′16″E，海岸线长约80 km，海域面积约2000 km²。曹妃甸原系滦南县南部海域一带状小岛，沙岛高潮时面积约4 km²，低潮时约20 km²，沙岛与大陆岸线之间是大片浅没海滩；甸头前500 m水深即达25 m，深槽处则为36 m^[10]，因其得天独厚的自然条件，加之近年来我国沿海经济的高速发展，唐山曹妃甸港区应运而生，并逐渐发展壮大(如图 1所示)。但环境生态问题也随之而来，2014年《河北省海洋环境状况公报》^[11]显示，曹妃甸海域水环境状况不容乐观：主要污染物为无机氮、活性磷酸盐，存在Pb超二类水质标准的情况；沉积物质量状况较好；底栖生物群落与多样性一般。因此研究并分析曹妃甸近岸海域大范围、高强度开发条件下的环境状况，对于指导近岸海域开发与管理极具现实意义。

图  1  研究海域位置及监测站位

Figure 1.  The research marine district location and stations

图 1是研究海域位置及监测站位，其中A为研究区域位于中国海域的相对位置；B为具体位置和监测站位，B中1、3、5、7、8、11、13和15号站位监测项目包括海水水质和沉积物，其余各站位仅监测海水水质。

文中采用2005年、2007年、2009—2012年6 a间的实测数据(数据来源于国家海洋局秦皇岛海洋环境监测中心站)，研究悬浮物、水相中的Pb含量和沉积相中的Pb含量的变化趋势以及重金属Pb的多相介质富集特性，采样时间均为每年的8月份丰水期，具体的监测站位如图 1所示。采样层次按水深＜10 m，取表层水样；10 m≤水深＜25 m，取表层和底层水样的要求处理。

所有水质和表层沉积物样品的采集、贮存、运输与分析等均按照GB 17378-2007《海洋监测规范》^[12]和《海洋调查规范》^[13]要求进行(其中2005年监测按照GB 17378-1998《海洋监测规范》^[14]要求进行)。样品采集^[12]以走航方式，用采水器对近岸表层海水进行水样采集，使用抓斗式采泥器采集海洋中表层沉积物样品；悬浮物样品应尽快放出采样器，装瓶同时摇动采样器，并去除杂质；重金属样品应立即用0.45μm滤膜过滤后，酸化至pH小于2，塞上塞子存于洁净环境；并采用冷冻或酸化后低温冷藏固定；每个站位采集3个平行样。悬浮物分析方法采用重量法(GB17378.4-2007)^[12]，经过滤、烘干、恒温、称重等过程后，测定悬浮物含量，悬浮物检测过程中采取平行样的方法对其进行质量控制；海水中Pb的分析方法为无火焰原子吸收分光光度法(GB17378.4-2007)^[12]，经pH调整、萃取分离、反萃取等过程，进行分析并加标回收，回收率在95%以上，相对标准偏差为1.2%~2.82%。沉积物样品自然风干后去除杂物，充分混匀、研磨后过100目尼龙筛。取0.2 g干样，采用无火焰原子吸收分光光度法(GB17378.5-2007)^[12]对沉积物中的Pb进行分析并加标回收，经分析得出回收率在95%以上，相对标准偏差为1.5%~3.8%。

趋势，特别是长期趋势是某一现象在较长时期内发展变化的总态势，是时间序列最基本的构成要素。通过趋势分析可以掌握影响时间序列的根本性因素。线性趋势是判断时间序列变化态势最简单最直观的分析方法^[15]，如式(1)所示。

式中：x_i是指时间序列x₁, x₂, …, x_i, …, x_n在时刻t_i对应的第i个值；n为变量的个数；A和B分别表示趋势系数和常数，可采用最小二乘法计算而得，如式(2)所示。

式中：和分别代表均值。

相关性分析是研究多元随机变量之间相关关系的一种统计学方法。两个变量之间的相关性为简单相关关系；而对于多元变量之间的相关性，则采用偏相关^[16]关系予以分析，其能够剔除其他相关变量的影响，正确地表示两个变量间的线性相关的性质与程度，真实地反映两个变量的本质联系。简单相关系数和N阶偏相关系数，分别如式(3)和(4)所示。

式中：X和Y分别为两个变量；n为变量中数据的个数。

式中：N阶偏相关系数可以利用N+2个变量组成的C_N+2²组简单相关系数计算而得^[17]，y, x, x₁, x₂, …, x_i, …, x_n代表n+2个变量。

本研究重点关注悬浮物及重金属Pb，采用前述趋势分析方法，研究了6 a间悬浮物、水相中的Pb含量和沉积相中的Pb含量的变化趋势，如图 2所示。

图  2  悬浮物、水相中的Pb含量和沉积相中的Pb含量的变化趋势

Figure 2.  The changetrend among suspended matter, water and sedimentary facies of Pb

图 2是悬浮物、水相中的Pb含量和沉积相中的Pb含量的变化趋势，图中：横坐标为监测年份，左侧纵坐标为悬浮物含量，右侧纵坐标为水相中的Pb含量和沉积相中的Pb含量。由3个变量的时间过程与趋势可见：在研究年份内，随着悬浮物的增加与减少，水相中的Pb含量呈现出基本相同的增减规律，而沉积相中的Pb含量则呈现出相反的趋势，其主要原因在于2007—2009年是本研究海域海洋工程大量建设的时期^[18]，吹填、溢流等施工作业引起海域悬浮物的大量增加，施工扰动亦会造成海域大量的底质泥沙起悬，而重金属物质极易附着于泥沙颗粒表面，也由此引起了水体中重金属含量的增高，同时底质的大量起悬也势必引起底质沉积物中重金属含量的减少。因此可见，除陆源污染之外，底质泥沙沉积物的再悬浮与再释放也是本区域重金属Pb污染的一个重要来源。

研究中考虑了悬浮物、水相中的Pb含量和沉积相中的Pb含量3个变量，因此式(4)可分别简化如下：

式(5)—式(7)是悬浮物、水相中的Pb含量和沉积相中的Pb含量3个变量的1阶偏相关系数计算公式，式中相关系数的下标悬、水、沉分别代表悬浮物、水相中的Pb含量和沉积相中的Pb含量3个变量；r_悬水、r_悬沉、r_水沉则为3个变量两两之间的简单相关系数，可采用式(3)求得。

将悬浮物、水相中的Pb含量和沉积相中的Pb含量3个变量近6 a的实测数据带入到式(3)和式(5)—式(7)中求得各相关系数(其值越大表明相关性越强，反之则弱；负值则代表负相关性)，如表 1所示。

表 1  悬浮物、水相中的Pb含量和沉积相中的Pb含量间的简单、偏相关系数

Table 1.  Simple and partial correlation coefficient among suspended matter, water and sedimentary phases of Pb

表 1是悬浮物、水相中的Pb含量和沉积相中的Pb含量间的简单、偏相关系数计算结果，由表 1中的简单相关可知：悬浮物与水相之间呈较弱的正相关；悬浮物和沉积相之间呈较强的负相关；水相和沉积相之间呈极弱的负相关。由表 1中的偏相关可知：同简单相关对比，水相和沉积相的变化最大，由极弱的负相关-0.00044变为较弱的正相关0.11716974，说明水相和沉积相受悬浮物影响较大。其次为悬浮物与水相之间的相关性，悬浮物与沉积相之间的相关性则变化不大，说明水相对于两者的影响较小。相关性研究结果表明，悬浮物是重金属污染因子Pb的重要载体^[1]，是其水相和沉积相变化最关键的影响因素。因此，采取必要的手段与措施，控制海域施工悬浮物质的大量产生是减轻重金属Pb污染的又一有效途径。

(1) 曹妃甸近岸工程海域近6 a悬浮物、水相中的Pb含量和沉积相中的Pb含量的富集趋势研究结果表明：水相中的Pb含量呈现出与悬浮物基本相同的增减规律，而沉积相中的Pb含量则反之，且悬浮物含量较高的年份恰恰是本研究海域海洋工程大量建设的时期，此时水体中的重金属Pb含量属典型的工程主导型。

(2) 简单相关与偏相关多元统计分析结果显示：水相和沉积相受悬浮物影响较大，由此可见，悬浮物是影响水相中的Pb含量和沉积相中的Pb含量变化的核心控制因素。

(3) 对于典型的重大工程建设海域，除控制陆源外，应采取必要的手段与措施，控制人为因素对底质沉积物的扰动，避免海域悬浮物质的大量产生，其有利于减轻重金属Pb的二次污染。