山东半岛南部丁字湾口外海底沉积物粒度时空变化及影响因素

王伟，汤世凯，胡艳萍，王红艳，石洪源，战超

1. 山东省海洋地质勘查院，烟台 2640042. 山东省第三地质矿产勘查院，烟台 2640043. 鲁东大学海岸研究所，烟台 264025

山东半岛南岸全新世最大海侵期间，海水沿河谷（沟谷）侵入半岛陆地内部形成了若干溺谷，实际上是山溪性入海河流的河口段，因河流入海水沙有限，河道口窄肚宽，在海洋动力控制下成为具有潟湖性质的河口湾。经过其上游河流泥沙长期充填，这些海湾最迟在20世纪70年代以前就已演变成大型潮汐汊道系统[1-2]，由狭窄的河口河道、树枝状潮汐汊道、宽平的潮间滩涂及潮上低地组成，低地后缘为埋藏的古海蚀崖，海湾水深除局部潮汐汊道外均不到6.0 m，若按照国际湿地公约的界定均属于湿地范畴。其中，位于青岛及烟台两市间、五龙河尾闾的丁字湾规模最大，形态最典型（图1），该海湾地区近几十年来经济发展快，人类活动影响大，这种宽阔滩涂、潮汐汊道式海湾对人类活动影响反应更为灵敏，其动力地貌演变特征对人类活动影响有着更为显著的表现。同时堆积在海湾的沉积物除了保留物源的信息外，还保留着人类活动对沉积物粒度影响的信息。

图1 山东半岛丁字湾地理位置图a为图2位置，b为图8位置。Fig.1 Location map of Dingzi Bay of Shandong Peninsulaa-Figure.2 location，b-Figure.8 location.

近些年来，前人对丁字湾的入海径流泥沙、动力、地貌特征、滩涂开发活动及其演变已有较多研究[1-6]，取得了一系列的成果。然而由于研究侧重方向不同，在沉积物粒度参数特征、粒径大小变化等方面，研究内容较少，精度较低；在沉积物同位素地质学方面的研究，例如沉积物中210Pb放射性比活度的垂直分布、沉积速率等方面也缺乏相关成果；同时在研究区也缺乏将沉积物粒度与沉积年代学相结合分析沉积环境变化的相关研究。因此通过分析海湾沉积物的粒度参数特征及粒度时空分布，一方面能够直观有效地反演该时期沉积物所处的沉积环境和沉积动力等条件，揭示沉积演变规律[7-10]，同时可以了解在历史沉积时期围海养殖、航道疏浚等人类工程建设活动对沉积地貌、粒度变化的影响[3,11-16]。

本文对丁字湾口外近岸浅海区域11站柱状沉积物样品进行粒度分析，并对其中2站柱状样进行210Pb比活度测试，计算粒度参数、沉积速率，划分沉积物类型，分析粒度参数空间分布，并通过与历史资料相结合，探讨丁字湾口外滨浅海海域海底沉积环境、沉积物侵蚀淤积演变规律以及流域人类活动在沉积物粒度记录中的响应，其结果可深化在全球变化背景下，暖温带山地河流、河口海岸沉积环境与地貌演变研究。

1 研究区概况

丁字湾位于山东半岛南部，是半岛最大最典型的溺谷河口-潮汐汊道海湾，高潮时面积143.75 km2、潮间带面积119.01 km2、湾口宽度6.0 km，呈口小肚大的形状，潮汐类型为正规半日潮，最高潮位4.27 m，最低潮位-0.5 m。潮流的主要运动形式为北西-南东向往复式流，湾内流速大于湾外。海浪浪向南东，夏季受台风影响在湾口外常形成大浪，波高一般大于1.5 m，最高可达3.5 m[2,4,17]。汇入丁字湾的河流有五龙河、白沙河等，沉积物来源除了河流输入以外，还有海岸侵蚀来沙和人工来沙，沉积物类型多样，沉积环境复杂。入该湾的五龙河在各入湾河流中规模也最大（长度128 km，流域面积2806.3 km2），且从1952年开始设有各入湾河流中唯一的水文观测站—团旺水文站，拥有最近60年的河流水文连续观测记录资料。研究区域位于丁字湾的湾口以外，为五龙河河口的口外滨海段，发育了以湾口为顶端的落潮三角洲，其地形以轴部高两翼低、上下游低中部高为特征，最高处高出上游潮汐汊道底部约10 m，实际上为落潮流控制下形成于口门以外的五龙河河口大型拦门沙体。人工地貌类型则主要有养殖池、水田、盐田等。

丁字湾海岸类型、地貌过程和海湾开发利用在各个同类海湾中颇具代表性。丁字湾是著名的古金口港所在，该港在青岛、烟台开埠以前是胶东半岛最重要的商港，直到20世纪上半叶才因为淤积严重而遭废弃。近50年来，丁字湾潮间滩涂和潮上湿地普遍被围垦，湾内5个岛屿中的4个因自然淤积和人工建坝而与陆地相连，海湾潮流运动的边界条件因而相应改变，尤其以1980年麻姑岛连岛大坝修建前后变化最为剧烈。与此同时，入湾河流五龙河流域内气候以及土地利用/土地覆被变化显著，各级河道上建设了大量水库（塘）、拦河坝和人工水源地，大量地面径流及地下水被开发利用，河流下游及河口水文水动力条件也相应改变。以上种种人类活动，必然会深刻影响到河口拦门沙及落潮三角洲的动力、沉积条件，进而改变其动力沉积与动力地貌过程。

2 样品采集与测试

根据丁字湾口外落潮三角洲（拦门沙）轴部高两翼低、上下游低中部高的特征，考虑到不同水深的地貌部位其动力、沉积速率也不同，于2017年5月在丁字湾口外浅海区平行海岸线和垂直海岸线方向上，采用震动活塞式柱状取样器进行取样，共采集11站柱状沉积物样品，相邻样品间距约2.5 km，柱状样长度均大于2 m。丁字湾周边相邻海域的沉积速率数值[18-20]均小于1 cm/a，多分布于4～8 mm/a，为了解研究区的沉积速率以及三角洲不同部位的沉积差异，对分别位于三角洲顶端和前缘的HZ20及HZ24开展210Pb的分布模式及沉积速率研究。对HZ20、HZ24两柱状样，按照4 cm间距分段切割，进行210Pb放射性比活度测试，对11站柱状样，按照4 cm间距分段切割，进行沉积物粒度测试，结合计算获得的HZ20、HZ24两站沉积速率值，研究丁字湾近半个世纪的沉积变化，具体采样站位见图2。

沉积物粒度分析由山东省第四地质矿产勘查院实验测试中心采用济南微纳颗粒仪器股份有限公司生产的Winner2008激光粒度分析仪进行，误差精度为0.01%，分析流程按照海洋地质地球物理调查规范[21]进行。测试结果按照Folk等[22-23]的方法进行沉积物的分类、定名。采用矩法公式[24-26]对沉积物进行粒度参数计算，用粒度参数分级进行定性描述（表1）。

表1 沉积物粒度参数分级[27]Table 1 Classification of grain size parameters

210Pb放射性比活度测试由南京师范大学同位素实验室完成，分析步骤如下：称取样品50～80 g，放入烘箱中，在100℃下恒温24小时。烘干的样品准确称重后装入测量盒中，密封15天，再置入仪器的测量室中测量。测量仪器分别为欧洲Canberra公司生产的BE3830型高纯锗探测器和DSC1000数字分析器组成的γ能谱仪，用来测量210Pb、226Ra。226Ra测量所用射线为214Bi的609.3keV γ射线；210Pb测量所用射线为46.5keV γ射线。流程按照土壤中放射性核素的γ能谱分析方法[28]进行。测试结果采用恒定初始放射性浓度模式（CIC模式），利用放射性衰变方程，通过作图得到深度H和lnN拟合线，其斜率即为S/λ，然后根据斜率获得沉积速率S[18-20,29]。

3 结果

3.1 沉积物类型与层序

根据研究区内11站柱状沉积物样品粒度分析结果，分别沿垂直海岸线方向（A-B剖面）和平行海岸线方向（C-D剖面）（图2）绘制柱状沉积物类型剖面分布图（图3、图4）。从图中可以看出海底柱状沉积物类型可划分为含砾砂、砂、粉砂质砂、砂质粉砂、粉砂、黏土和泥等7种类型。

图2 丁字湾口外海域柱状样站位分布图Fig.2 Location of columnar sediment sampling sites

图3 A-B剖面各站位沉积物类型组合Fig.3 Sediment type assemblages of each station in A-B section

图4 C-D剖面各站位沉积物类型组合Fig.4 Sediment type assemblages of each station in C-D section

丁字湾口外滨海不同的海域受控于不同的沉积水动力环境，各站位沉积物类型不同，从A-B剖面图中可以看出，除HZ20以外，由HZ21至HZ26在远离海岸水平方向上沉积物颗粒逐渐变细。同样在C-D剖面图中，以中央潮汐汊道为界，沿两侧水平方向上沉积物颗粒逐渐变细。可以看出在平面上沉积物向四周边缘粒度逐渐变细，这一特征与此处发育的水下扇形三角洲海底地貌相吻合，指示了由湾口向四周水动力环境逐渐减弱的现象。

同样，在不同的历史沉积时期，各柱状样在垂向上亦形成不同的沉积组合（韵律），A-B剖面中，HZ20柱状样下部为砂，上部为粉砂质砂-黏土（泥）沉积韵律组合；HZ21柱状样下部为含砾砂，上部为砂；HZ22柱状样下部为含砾砂，上部为砂夹粉砂质砂；HZ23柱状样沉积物类型单一，以砂为主，夹砂质粉砂；HZ24柱状样主要为粉砂-黏土（泥）的沉积韵律组合，顶部发育少量砂；HZ25柱状样上部为粉砂，下部为黏土（泥）。HZ26柱状样上部为粉砂夹黏土，下部为黏土（泥）。C-D剖面中，HZ27柱状样上部为粉砂-泥（黏土）沉积韵律组合，下部为粉砂质砂；HZ22柱状样下部为含砾砂，上部为砂夹粉砂质砂；HZ12柱状样沉积物类型单一，主要为含砾砂；HZ5和HZ8柱状样沉积物类型一致，为黏土。

由以上沉积组合（韵律）可以发现，垂直海岸方向上，HZ23所处的海域沉积环境动力相对稳定，并以此为界，两侧海域具有不同的沉积动力环境，近岸的HZ20-HZ22海域早期沉积环境动力强于晚期，远岸的HZ24-HZ26所处海域则相反，早期沉积环境动力弱于晚期。平行海岸方向上，潮汐汊道西南海域处于早期强、晚期弱的水动力环境；潮汐汊道东北沉积动力环境变化不大，靠近潮汐汊道处于较强的动力环境，远离潮汐汊道处于较弱的沉积环境。

3.2 沉积物粒度参数

粒度参数分析结果见图5、图6。平均粒径表示粒度分布的集中趋势，能指示源区物质的远近程度和搬运介质的动力大小[23-24,30-31]。A-B剖面各站位平均粒径变化范围不同，其中，HZ21、HZ22和HZ23平均粒径及变化范围幅度均较小，分别为0.71～1.89Φ（平均值1.23Φ）、1.30～3.24Φ（平均值2.32Φ）、2.55～3.48Φ（平均值2.92Φ）；HZ25、HZ26平均粒径Φ值较大，但变化范围幅度小，分别为6.57～8.61Φ（平均值7.38Φ）、7.02～8.74Φ（平均值7.98Φ）；HZ24平均粒径及变化幅度大，为2.54～8.37Φ（平均值6.35Φ），由下至上整体表现为逐渐减小；HZ20下部平均粒径及变化幅度均较小，为1.53～2.22Φ（平均值1.74Φ），上部变化范围大，为2.22～8.35Φ（平均值5.85Φ），表现为粒径大小交替分布的特征。水平方向上除了HZ20异常外，由HZ21至HZ26方向平均粒径平均值逐渐增大，颗粒逐渐变细，反映了远离海岸方向搬运介质动力逐渐降低。垂直方向上HZ20早期沉积动力环境相对较强且稳定，后期在相对强和弱的沉积动力环境中交替变化，HZ24沉积动力环境整体呈现强弱交替变化特征。其余各站位，HZ21-HZ23沉积动力环境相对较强且较稳定，HZ25-HZ26沉积动力环境整体相对较弱，后期较前期略强。

图5 A-B剖面中沉积物粒度参数变化Fig.5 Variation of sediment grain size parameters in A-B section

图6 C-D剖面中沉积物粒度参数变化Fig.6 Variation of sediment grain size parameters in C-D section

C-D剖面中各站位平均粒径Φ值及变化范围也各不相同，分别为HZ27平均粒径3.55～8.36Φ（平均值5.99Φ），HZ22平均粒径1.30～3.24Φ（平均值2.32Φ），HZ12平均粒径1.00～2.30Φ（平均值1.76Φ），HZ8平均粒径8.83～9.70Φ（平均值9.41Φ），HZ5平均粒径9.79～10.48Φ（平均值10.18Φ），水平方向表现为由潮汐汊道向两侧平均粒径Φ值逐渐增大的趋势，指示了搬运介质动力向两侧逐渐降低。垂直方向上表现为HZ27水动力环境明显减弱外，其余站位水动力环境变化不大，HZ12、HZ22沉积动力环境相对较强且稳定，HZ5、HZ8沉积动力环境相对较弱且稳定。

沉积物的分选系数代表着沉积物颗粒大小的均匀程度[23-24,30-31]，反映了水环境动力对颗粒的分选程度，A-B剖面中HZ20-HZ24分选系数范围值较大，分别为0.17～1.75（平均值0.87）、0.19～1.33（平均值1.00）、0.13～1.72（平均值0.78）、0.21～1.40（平均值0.51）、0.21～1.26（平均值0.79）。HZ25和HZ26分选系数范围较小，分别为0.44～0.91（平均值0.73）、0.45～0.78（平均值0.61）。根据表1中分选系数分级及定性描述，A-B剖面分选性属于较好—中等级别，表明从物源区到沉积区，搬运路径的增大导致分选系数变小、分选性变好。

C-D剖面中各站位分选系数变化也不尽相同，HZ5和HZ8分选系数范围小，分别为0.14～0.50（平均值0.39）、0.32～0.51（平均值0.40）；HZ12、HZ22和HZ27分选系数范围较大，分别为0.55～1.10（平均值0.85）、0.13～1.72（平均值0.78）和0.22～1.15（平均值0.68）。根据表1中分选系数分级及定性描述可知，C-D剖面分选性属于好—中等级别，并表现为由剖面中间向两侧分选系数逐渐变小、分选性变好的趋势。

偏态能够判别粒度曲线的对称程度，反映颗粒被改造程度。偏态为正时，沉积物粒径偏向粗粒，偏态为负时，沉积物粒径偏向细粒[23-24,30-31]。A-B剖面中除HZ20下部和HZ22上部偏态异常外，各柱状样整体偏态值上下较为一致，范围为-0.90～4.46，剔除异常后总体平均值为0.79，偏态定性描述为正偏，颗粒粒径集中在粗粒的一侧。C-D剖面中偏态值范围为-6.86～6.80，平均值0.11，偏态定性描述为0值两侧近对称分布，粗、细颗粒含量近似一致。

峰态是衡量粒度分布曲线的宽窄尖锐程度[23-24,30-31]，峰态规律同偏态相近，A-B剖面中HZ20下部和HZ22上部峰态值异常，其余各柱状样峰态变化范围不大。整体峰态范围为1.91～55.76，剔除异常后总体平均值为3.93，峰态定性描述为尖锐。C-D剖面中HZ22和HZ27上部峰态值异常，整体范围为1.03～55.76，剔除异常后总体平均值为4.13，峰态定性描述为尖锐。

3.3 沉积速率

结合丁字湾历史水深资料[17]和2015年10月中国航海图书出版社出版发行的海阳港区及附近海图，研究区等深线位置自1950年以来尚未发生明显变化，可认为选取的两个站位是连续性沉积，选用210Pb测年方法计算沉积年龄是合适的。HZ20柱状样中210Pbex随深度分布变化如图7所示，38 cm以上层位210Pbex随深度有明显的指数衰减趋势，该深度内为210Pb的衰变段，38～54 cm层位210Pbex有随深度增加趋势，54 cm以下层位210Pbex又呈降低趋势。由此可知本次采集的沉积物柱样长度未达到210Pb本底层。因此选取38 cm以浅的10个样品，根据线性拟合结果，计算该柱状样2～38 cm层位的平均沉积速率为0.51 cm/a。这一计算层位沉积历时75年左右，年代跨度为1940—2017年。

HZ24柱状样中210Pbex随深度分布变化如图7所示，210Pbex在表层6 cm内活度异常，表明扰动层发育，可能由于底栖生物扰动作用所致。6～42 cm层位210Pbex随深度有明显的指数衰减趋势，该深度内为210Pb的衰变段，42 cm以下层位210Pbex活度值分布规律性差，同样采集的本沉积物柱样长度未达到210Pb本底层。因此选取6～42 cm层位的10个样品，根据线性拟合结果，计算该柱状样6～42 cm层位的平均沉积速率为1.19 cm/a。这一计算层位沉积历时37年左右，年代跨度为1980—2017年。

图7 钻孔HZ20和HZ24柱状沉积物210Pbex垂直分布曲线Fig.7 Vertical distribution curve of 210Pbex

两站位沉积相同厚度的沉积物，沉积历时不同，与分别处于不同的沉积水动力环境密切相关，HZ20整体处于相对高能水动力环境，易侵蚀冲刷，而HZ24站位相反，整体处于相对较弱水动力环境，易接受沉积。

4 讨论

对丁字湾口外滨海而言，在丁字湾落潮流控制下发育了大型落潮三角洲（拦门沙），落潮流强弱及其携带泥沙数量、沉积物粒度粗细受丁字湾纳潮量和五龙河入湾水沙制约，此外还受人类活动尤其是工程建设活动的深刻影响。

历史资料表明，1980年以前，五龙河流域降雨充沛，且人类工程活动建设（水库）主要集中于五龙河上游，植被破坏，大量山地来源沉积物进入五龙河，入海径流、泥沙量明显增大，尤其是粗颗粒泥沙比例较大。1980年以后，五龙河流域降水量显著减少，受到水库等水利设施的围截，入海径流、泥沙量也相应显著减少[3]。同时，80年代初人类开始对湾内滩涂进行围垦，1970年前后丁字湾沿岸有少量的盐田，1981—2000年，养殖池面积剧增，以新修建养殖池或其他类型地貌向养殖池转换为主[5]，尤其是盐田和潮滩等从80年代开始显著减少。2000年以后，养殖池面积增速放缓，人类活动的影响逐渐减弱，海湾地貌格局趋于稳定，变化不大，这在不同年份的海湾地貌资料上有着很好的反映（图8）。

图8 丁字湾沉积地貌演化示意图1971年和1981年海湾地貌图修改自文献[5,32]，2000年和2016年海湾地貌分别根据美国陆地卫星TM、L8多波段遥感图像绘制。Fig.8 Schematic maps of sedimentary landform evolution of Dingzi BayThe 1971 and 1981 gulf topography maps modified from references [5, 32], the 2000 and 2016 gulf topography maps based on US Landsat TM,L8 multi-band remote sensing imagery.

这些诸多人类活动深刻影响着原有的地貌边界条件、动力过程及其演化趋势，进一步影响着丁字湾沉积过程及沉积物粒度变化情况。1980年以前，丁字湾总体受人类活动影响较小，五龙河径流量、输沙量较大，海湾处于由落潮不对称向涨潮不对称转换[3,5-6]，海湾地貌属于潮汐汊道窄深、两侧潮间滩涂宽广的空间分布格局，纳潮量较大，涨潮历时长于落潮历时，涨潮流速小于落潮流速，沉积物主要表现自然沉积。在强劲落潮流冲刷、五龙河携带的泥沙搬运距离较短以及外海传来的波浪在此变形及破碎、水流能动性强等因素共同作用下，落潮三角洲范围成为高能水动力沉积环境，近岸湾口处三角洲顶端的HZ20下部（图3中该站位24 cm以下）砂、粉砂质砂等较粗粒沉积物是这一沉积时期很好的响应。向三角洲前缘方向外海海域水深增大，五龙河携带的泥沙搬运至此处距离增大，落潮流能量降低，且潮流流速随水深递增而明显减小，HZ24底部（图3中该站位44 cm以下）的泥质沉积物相对较细。

1980年以后，沉积物在自然条件沉积的前提下，受到了人类活动的深刻影响。丁字湾滩涂围垦导致海湾纳潮空间显著减小，细颗粒泥沙的淤积空间变小，海湾涨潮历时变短，落潮历时变长，涨潮流流速加大，落潮流流速变小[32]，水动力作用受限，再加上此时五龙河入湾径流、泥沙量都迅速减少[7]，近岸湾口处三角洲顶端的HZ20上部（图3中该站位24 cm以上）沉积物较下部明显变细，垂直方向上呈粗细交替变化，位于三角洲前缘海域的HZ24则受涨潮流及波浪作用，其中上部（图3中该站位44 cm以上）沉积物较底部略粗，垂直方向上也呈粗细交替变化。

5 结论

（1）丁字湾口外浅海区海底48 cm深度内发育含砾砂、砂、粉砂质砂、砂质粉砂、粉砂、黏土和泥等7种沉积物类型。平均粒径Φ值由三角洲顶端向前缘、潮汐汊道向两侧方向逐渐增大，分选性属于较好—中等级别，偏态为正偏或近对称分布，峰态为尖锐型。沉积环境水动力强度由湾口向湾外、由中央潮汐汊道向两侧海域逐渐降低。历史沉积时期三角洲东北翼及前缘以外的海域为较弱且稳定的沉积动力环境，三角洲内潮汐汊道周边海域为较强且稳定的沉积动力环境，三角洲顶端、前缘及西南翼海域为强弱交替不稳定沉积动力环境。

（2）对位于落潮三角洲顶端及前缘的HZ20和HZ24柱状样进行210Pb放射性比活度分析，HZ20样品最近80年以来的沉积速率为0.51 cm/a，揭示三角洲顶端（潮汐汊道）发育于高能水动力环境，易遭受侵蚀冲刷；HZ24样品最近40年以来的沉积速率为1.19 cm/a，揭示三角洲前缘发育于相对低能水动力环境，易接受沉积。

（3）最近40年以来的流域气候和海湾滩涂围垦深刻影响丁字湾动力地貌、沉积动力及沉积物粒度演化模式。1980年以前，人类活动影响有限，海湾沉积物以自然沉积为主，在落潮流主导下潮汐汊道及周围海域沉积颗粒较粗，三角洲前缘外海域沉积颗粒较细。1980年以后，丁字湾滩涂遭围垦严重，海湾纳潮空间显著减小，五龙河入湾径流、泥沙量也都迅速减少，潮汐汊道及上游其水动力影响降低，潮汐汊道及相邻两侧海域沉积物明显变细，三角洲前缘及外海域受波浪作用叠加影响沉积物略微变粗，垂直方向上均呈粗细交替变化。