北黄海大鹿岛附近海域海洋牧场本底地质环境分析*

杜梦迪，张建兴，宋永东，栾振东,3**，阎 军

(1.中国科学院海洋研究所，中国科学院海洋地质与环境重点实验室，山东青岛 266071；2.中国科学院大学，北京 100049；3.中国科学院海洋大科学研究中心，山东青岛 266071)

我国是海洋大国，海域面积十分辽阔，海洋资源得天独厚，渔业产量常年位居世界首位[1,2]。然而，当前我国海域渔业正面临着渔业资源衰退和海洋生态环境污染等问题，传统渔业方式已难以满足渔业资源可持续发展和海洋环境现状的要求。现代海洋牧场是传统渔业转型的新动力，既可以实现海洋环境保护，又能保证渔业资源高效产出[3]，是应对我国渔业资源衰退和海洋生态环境恶化的一个重要举措，是当前乃至今后较长时期内海洋渔业发展的一个重要方向。

在目前海洋牧场快速发展的背景下，我国浅海近岸海洋牧场建设项目不断增加。丹东是国家级海洋牧场示范区之一，近年来，在丹东大鹿岛附近海域海洋牧场等海上工程建设项目逐渐增多，因此有必要对研究海区进行海洋地质环境综合调查。声学探测技术已成为海洋调查最常用和最成熟的手段，获取高分辨率的海底地形地貌资料和海底底质信息资料是进行海底地质环境综合调查的首要内容，也是保障海洋牧场建设顺利进行的重要前提。此外，表层沉积物调查是了解海底底质类型和沉积环境的重要方法。

虽然前人已经对北黄海海底地形地貌[4-6]和沉积特征[7-11]进行了相关方面的调查和研究，但先前的研究多集中在北黄海西部、山东半岛和辽东半岛等近岸海域，而对北黄海北部海域研究较少，且前人的研究多基于大范围海区，对北黄海大鹿岛附近小目标区域精细测量的研究欠缺，无法为研究区内海洋牧场建设提供可靠的指导。因此，本研究对丹东大鹿岛附近海域海底地形地貌、海洋水质环境及沉积特征开展了详细研究，有助于对海域内海洋牧场资源环境进行系统评估，对丹东地区新型海洋牧场的构建具有一定的参考意义。

1 材料与方法

1.1 研究区概况

北黄海位于中国东部陆架最北端，是典型的半封闭陆缘海，环绕山东半岛、辽东半岛和朝鲜半岛，西北通过渤海海峡与渤海相邻，南以成山角与朝鲜长山串的连线与南黄海相接[5]。北黄海海底整体呈“簸箕”状，海底相对平缓，平均水深为38 m，最大水深可达86 m，黄海海槽为水深最深处，由中央海槽向外水深逐渐变浅[12]。研究区位于北黄海北部大鹿岛附近海域(图1)，水深较浅，最大水深为29 m，地形较为平坦。在构造上，研究区属于海洋岛隆起构造单元，北部与辽东隆起相邻，南部与北黄海盆地相接。

据统计发现，北黄海沉积物的来源有两类:一类来自沿岸河流，另一类来自海洋。径流在陆相沉积物的输送中起着重要的作用，辽东半岛近岸河流汇入北黄海的输沙量可达7×106t·a-1，其中鸭绿江和大洋河占输沙量的80%[12,13]。海洋的沉积物主要来自渤海和南黄海，通过海洋环流进行物质输送和物质循环。黄河虽然不直接汇入北黄海，但是因其巨大的入海水沙量，被认为是北黄海陆架沉积的主要物源，占北黄海总沉积物的66%-80%[14]。

北黄海的潮汐类型以规则半日潮为主，辽东半岛自旅顺至鸭绿江河口的北黄海沿岸均为规则半日潮；潮差分布呈由东北向西南逐渐减小的趋势，鸭绿江口大潮平均潮差为5 m，至成山角潮差不到1 m[15]。研究区属于正规半日潮，潮差为4-5 m，潮流流速较大，流速为0.6-1.0 m·s-1。北黄海的潮流强盛，海流较弱，海区内环流主要由辽南沿岸流和黄海暖流两支组成[10]，其中本研究区域受辽南沿岸流的影响较大。

图1 研究区位置与调查站位

1.2 数据来源

利用HY1601单波束测深系统、SBP/AAE浅地层剖面系统、声学多普勒流速剖面仪(Acoustic Doppler Current Profiler，ADCP)、温盐深(Conductivity Temperature Depth，CTD)测量系统和表层沉积物取样等调查手段对研究区开展综合地质环境调查，获得研究区水深数据、浅地层剖面数据、海水理化数据和表层沉积物样品等资料。

1.3 方法

1.3.1 单波束地形测量

调查前需对研究区进行测线布设，单波束测量布设按照垂直于等深线的原则，测区共包含12条主测线和2条联络线，主测线布设间隔为3.3 km，联络线垂直于主测线进行布设，如图1所示。获取的水深数据利用CARIS HIPS & SIPS软件进行数据处理，进行姿态校正、声速校正和潮位校正，并剔除噪声和干扰等异常数据，得到高分辨率水深数据，后经Surfer等成图软件绘制研究区的水深图及地形剖面图。浅地层剖面数据通过SonarWIZ软件进行成图和分析，最终导出较高分辨率的浅地层剖面图像。

1.3.2 海水理化性质测定

海水理化环境调查和沉积物取样站点见图1。为掌握研究海域的海洋水质环境，于2020年秋季11月在北黄海北部大鹿岛附近海域共设置13个站点，在各调查站点利用CTD测量系统和HQ4300多参数水质分析仪(美国HACH公司)进行水温、盐度、溶解氧含量、浊度、叶绿素含量、硝酸盐、亚硝酸盐、铵盐、活性磷酸盐等理化参数的测量。海水样品的采集、保存与检测分析均严格按照《海洋监测规范》(GB 17378.4-2007)[16]的要求进行。根据测定的海水理化参数以及近年来北黄海海水环境状况，选择溶解氧、无机氮(硝酸盐、亚硝酸盐和铵盐中N的总量，单位：mg/L)、活性磷酸盐这3个项目对海域内海水水质状况进行系统评价，水质标准采用《海水水质标准》(GB 3097-1997)[17]，详见表1。

表1 水质参数评价标准

1.3.3 沉积物粒度测定

为查明研究区的沉积物特征，在调查区进行表层沉积物采集，首先利用锚式采泥器抓取各站点底质样本，然后进行室内实验分析。测量前先加入六偏磷酸钠分散，然后经超声波振荡分散后上机测试，沉积物粒度组成采用筛析法和CILAS-1190激光粒度仪(法国CILAS公司)联合进行测量。粒度实验结果在基于MATLAB的粒度分析软件AnalySize中进行分析计算，并按照Shepard沉积物分类方法对沉积物进行定名。粒级标准采用尤登-温德华氏等比制Φ值粒级标准，沉积物的命名采用Shepard的沉积物粒度三角图解法[18]。沉积物的平均粒径(Mz)、分选系数(δ)、偏态(Sk)和峰态(Ku)通过Folk和Ward粒度参数公式计算，具体如下：

Mz=(Φ16+Φ50+Φ84)/3，

δ=(Φ84-Φ16)/4+(Φ95-Φ5)/6.6，

Sk=[(Φ84+Φ16-2Φ50)/2(Φ84-Φ16)]+

[(Φ95+Φ5-2Φ50)/2(Φ95-Φ5)]，

Ku=(Φ95-Φ5)/2.44(Φ75-Φ25)，

式中，Φ=-log2D，D为沉积物粒径(mm)。

2 结果与分析

2.1 地形地层特征

2.1.1 海底地形地貌特征

本次调查利用单波束测深系统获得的水深数据绘制出研究区高分辨率水深地形图，结果如图2所示。研究海域水深较浅，水深为1-29 m，平均水深为16 m，最深处位于研究区南部，最浅处位于西北部大鹿岛近岸区域，在123.7°-123.9° E内发育有NE-SW向沙脊，海底地貌以水下侵蚀堆积地貌为主。区域整体地形较为平坦，表现为典型的近岸缓坡地形，坡向由北向南。从A-A′与B-B′的海底地形剖面图可以看出，总体坡度小于12‰，水深5 m内地形起伏较为平缓，坡度小于10‰；水深大于5 m坡度变陡，坡度为12.5‰(图3)。

图2 研究区地形图

图3 海底地形剖面图

2.1.2 浅地层剖面特征

根据北黄海东北部的JPL9-1高分辨率浅地层剖面的反射特征，结合海底表层沉积物取样结果，可以分析出区域内海底浅地层结构、沉积层厚度和分布范围。由图4可知，在浅地层剖面图像上，以海底反射面R0和声学反射面R1为分界面，可以划分为U1和U2上下两套声学地层。其中U1声学地层振幅强、连续性好、全区可追踪，与下伏地层呈侵蚀不整合接触，厚度约为6 m，其厚度在研究范围内变化不大，层间具有平行反射结构，无明显构造变形，以加积反射层为特征，沉积物以砂、粉砂和黏土为主，为全新世以来的陆架环境海相沉积[19]。U2声学地层与U1声学地层的地震相特征明显不同，U2层声学信号较弱，说明上覆地层含有粗粒物质，阻挡了部分声学信号的传播，且该层信号较复杂，连续性较差，厚度为10-15 m，近岸方向厚度逐渐减小，发育有明显的信号异常区。U1层为本次浅地层划分的主要声学地层，该地层的岩性和分布受全新世以来的海平面变化和潮流因素的控制。

图4 北黄海东北部海域JPL9-1浅地层剖面图

2.2 海水水质状况

调查水域各站位表、中、底3层海水的水文和营养盐情况如表2所示。从调查结果可以看出，温度为3.23-12.94℃，最大值出现在10站位，最小值出现在9站位，各站位海水温度随水深增大而升高；溶解氧浓度为8.13-12.25 mg/L，最大值出现在9站位，最低值出现在13站位，海域溶解氧的平均含量为9.48 mg/L，各站点溶解氧含量随水深增大而逐渐降低；底层水盐度为26.33-29.88，变化不大，平均值为28.05；浊度为1.88-222.72 NTU，浊度变化较大，最大值出现在1站位，最小值出现在12站位；叶绿素含量为0.30-4.62 μg/L，最大值出现在1站位，最小值出现在10站位；无机氮(以N计)含量为0.170-0.439 mg/L，最大值出现在5站位，最小值出现在12，13站位；活性磷酸盐含量为0.004-0.008 mg/L，最大值出现在站位10，最小值出现在3，6,8站位。调查海区海水温度与盐度空间分布规律趋于一致，由大鹿岛沿岸向海方向逐渐增大，溶解氧、浊度和叶绿素的空间分布特征与之相反；无机氮和活性磷酸盐含量基本稳定，变化不大。

表2 海水理化参数和营养盐含量

续表

分析表明，北黄海大鹿岛附近海域的溶解氧含量均符合Ⅰ类海水水质标准；75.0%区域的无机氮符合Ⅱ类水质标准，16.7%区域的无机氮符合Ⅰ类水质标准，8.3%的无机氮符合Ⅳ类水质标准；各站位活性磷酸盐含量均符合Ⅰ类海水水质标准。综合分析，大鹿岛附近海域海水水质良好，主要符合Ⅱ类海水水质标准，符合海水养殖的水质要求。

2.3 沉积特征

2.3.1 表层沉积物类型与分布

由表3可知，研究区表层多为黑色砂质沉积物，底质类型均为砂质沉积，沉积物类型无明显的空间变化。各采样点之间的粒度参数差异不大，其中，砂的组分为89.29%-96.85%，粉砂含量次之，黏土最少。平均粒径的分布为2.65-2.98Φ，平均值为2.83Φ；分选系数为0.42-0.79，平均值为0.53；偏态均为正偏，平均值为0.22；峰态的平均值为1.20。根据Shepard分类命名方法，北黄海表层沉积物均为砂，砂的含量很高，均大于89%，调查区各站位表层沉积物分类和命名详见图5。

表3 研究区表层沉积物组分含量及粒度分析结果

图5 研究区表层沉积物分类与定名

2.3.2 沉积物粒度参数空间分布

从图6可以看出，研究区沉积物组分的空间变化不大，3，6，9号站位黏土与粉砂占比略高，但整体上砂的组分占绝对优势，沉积物粒径较粗。其他粒度参数中，分选系数、偏态和峰态的空间分布关系较为类似，这三者与平均粒径的空间变化规律相反。其中，测区南北两侧平均粒径较大，而测区的西北角和东部9号站位附近平均粒径偏小，整体变化趋势为由3，6，9号站位的较小值向南北两侧逐渐增大。相反地，测区内3，6，9号站位附近粒度参数表现为分选系数偏大，偏态偏正，峰态偏大；而测区南北两侧分选系数较小，偏态偏负，峰态偏小。

图6 研究区表层沉积物组分含量及粒度参数空间分布

3 讨论

3.1 研究区沉积物粒度参数特征

沉积物的粒度参数及其空间分布是沉积物的基本属性之一。沉积物粒度参数是以一定数值定量地表示沉积物的粒径频率分布[20]，是用来研究沉积物的沉积特征及沉积环境的统计学方法。粒度参数可以指示沉积物的物源条件、水动力条件等信息，其空间分布特征是沉积物物质循环和外动力输运的结果。本研究发现，沉积物的来源主要与平均粒径和分选系数密切相关，而沉积环境主要反映在偏态和峰态上[21]。结合研究区沉积物粒度参数分析，研究区沉积物分选较好，呈正偏态和中等尖锐峰态，峰态的数值变化较大。

3.2 沉积物类型及成因分析

已有的研究表明北黄海北部沉积物类型主要为砂、粉砂质砂和砂质粉砂，且呈自鸭绿江口向西粒径变细的趋势[22-25]，其中砂质沉积物占绝大比例，主要分布在123.3° E以东的海域[10]。由此可见，本研究对北黄海大鹿岛附近海域表层沉积物类型的研究结果与前人的研究结果能较好地吻合。

沉积物粒度参数和分布特征与沉积物的物源条件、水动力条件等多种因素有关。北黄海的沉积物主要为陆源沉积，其沉积物源主要来自黄河、鸭绿江等沿岸河流[14,26,27]，另外还包括近岸和海底侵蚀。北黄海西部沉积物主要来源以黄河入海物质为主，东部沉积物以鸭绿江入海物质为主。研究区属于北黄海北部区域，根据其粒度分析结果和沉积物类型，推测其砂质沉积的主要物源为鸭绿江等沿岸河流的入海物质。

研究区砂质沉积形成与辽南沿岸流及区域内强潮流环境密切相关，海底地形地貌、古环境及海平面变化对其也具有一定影响。辽南沿岸流由鸭绿江和沿岸河流冲淡水组成，一直沿辽东半岛东南岸由鸭绿江口流向渤海，方向终年稳定不变。夏季降雨导致鸭绿江径流量增加，辽南沿岸流流速超过10 cm/s，冬季受东北风影响，长山列岛海域流速也接近10 cm/s[28]，主要为北黄海北部近岸海域分布的砂质沉积提供了物质来源。区域内强的潮流环境对砂质沉积的形成也起到重要作用，辽东半岛东南岸高速的潮流可以为河流携带入海的悬移质物质的抬升和搬运提供动力，推动粗粒物质的分选富集，促进砂质沉积的形成。北黄海北部近岸区海底地形较为平坦，有利于沉积物的稳定输运与沉积，故研究区的沉积组分及空间分布特征也较为稳定。古气候及海平面变动对沉积特征的形成也具有一定的影响，北黄海长期处于滨海相环境，末次盛冰期以来的沉积被持续地侵蚀和改造，细颗粒物质被搬运，粗颗粒物质被留在原地[26]，因此在北黄海地区形成了大量的砂质沉积。

本研究依据北黄海大鹿岛附近表层沉积物类型、粒度参数及空间分布特征，分析了沉积物的物质来源和水动力条件，并总结出了研究区沉积特征的形成原因，分析认为研究区内沉积物整体粒度较粗，以砂质沉积为主，其物质来源主要为鸭绿江等沿岸河流的入海物质；沉积特征还受到辽南沿岸流、海洋潮流、海底地形地貌及古气候等多种因素的影响。

3.3 海洋牧场本底地质环境分析与评估

为促进海洋牧场的选址与建设，从水深地形、水质环境、海底底质类型等3个方面对研究区海洋牧场本底地质环境进行分析与评估。第一，从水深地形条件来看，一般而言，水深适宜、地形平坦的区域有利于海洋牧场工程的顺利开展，考虑到光照和温度对人工鱼礁的影响以及波浪对海洋牧场构建物的冲刷作用，海洋牧场的海区一般位于水深5-30 m的海域[29]。北黄海大鹿岛附近海域平均水深为16 m，地形为典型的近岸缓坡地形，适于海洋牧场的构建和布局。第二，从水质环境来看，海水水质是海洋牧场选址的一个重要环境参数，养殖水域的海水环境应满足不受污染，透明度好，不浑浊，水质应符合渔业水质标准[30]。大鹿岛附近海域海水主要处于国家Ⅱ类海水水质标准，透明度较好，水温适宜，理化环境优越，营养盐含量丰富，具有沿岸河流携带的丰富有机质，适于多种海洋生物的生长，尤其是丹东海域的优势种中国蛤蜊(Mactrachinensis)的活动受水温影响显著，在一定温度范围内，中国蛤蜊的跳跃频率和跳跃高度随水温的升高呈上升趋势。调查海域海水的温度等环境参数适合中国蛤蜊的活动和生长，该海域具备建设丹东特色海洋牧场的海水水质环境。第三，从海底底质类型来看，在选择海底底质时一般选取海底表层含有薄层泥沙的较硬质底质，而避免人工鱼礁投放后因底质太软而沉入海底，大鹿岛附近海域表层沉积类型均为硬质砂质沉积，淤泥含量较少，该海域的海底底质类型有利于人工鱼礁的投放和海洋牧场构建物的布设。

综合多种环境因素对大鹿岛附近海域的海底地质环境进行综合评估，最终可为丹东特色滩涂型海洋牧场的选址、布局和建设提供科学指导。

4 结论

本研究对北黄海大鹿岛附近海域海底地形地貌、海水水质状况、沉积物类型及分布等方面进行了地质环境综合调查和分析，得出了以下结论：(1)北黄海大鹿岛附近海域水深较浅，平均水深为16 m。区域整体地形较为平坦，表现为典型的近岸缓坡地形，地貌以水下侵蚀堆积地貌为主；(2)浅地层剖面资料显示，研究区海底面以下沉积物为松散的砂质沉积物，U1声学地层信号清晰连续，全区可追踪，沉积层厚度约6 m且分布稳定，为全新世以来的陆架环境海相沉积；(3)北黄海大鹿岛南部海域海水水质良好，主要处于Ⅱ类海水水质标准，该区域水质环境符合海水养殖的水质要求；(4)研究区表层沉积物类型均为砂质沉积，大部分海域沉积物分选较好。研究区沉积物类型及分布是沉积物物源和水动力环境等多种因素共同作用的结果。研究区沉积物主要来源于大洋河和鸭绿江等沿岸河流的入海物质，沿岸流、海洋潮流及古气候对沉积物分布具有明显控制作用；(5)研究海域地形地貌和沉积物反映出北黄海大鹿岛附近海域为地形平坦的近岸砂质浅滩沉积且海水水质较好，整体工程环境良好，适于进行丹东滩涂型海洋牧场的构建。