辽东湾晚第四纪层序地层

陈晓辉，孟祥君，李日辉

1. 青岛海洋地质研究所，青岛 266071

2. 青岛海洋科学与技术国家实验室，海洋矿产资源评价与探测技术功能实验室，青岛 266071

辽东湾是一个西、北、东三面被陆地包围的半封闭海湾，南面与渤海相通，水深普遍小于30 m，坡度较缓，地形自湾顶及东西两侧向中央倾斜。该地区构造上位于中朝准地台华北断坳下辽河断陷带，古新近纪以来一直沉降；每年接受大量辽河、大辽河等周边河流搬运入海的陆源碎屑物质，第四系厚度为300～500 m[1]。受海平面变化的影响，海陆交互沉积是这一地质时期的重要沉积特征。

近年来国内外学者采用沉积物粒度[2-5]、矿物学[6-7]、古生物学[8-10]或多指标结合[11-13]及浅地层剖面与岩芯相结合[14-16]等方法，对辽东湾及周边陆续开展了物源、沉积地层、晚第四纪沉积演化等方面的研究，取得了一系列的研究成果，然而这些研究多集中于表层沉积物，针对沉积地层或环境演化的研究多局限于晚第四纪晚期以来，如孙荣涛等[10]识别出了辽东湾南部钻孔MIS4 期以来的沉积层序。针对辽东湾晚第四纪地层的研究仍然比较薄弱，仅吴建政等[14]识别出了辽东湾北部晚第四纪以来的3 个海相层与陆相层，然而总体来讲对地层表现的不同特征尚缺乏系统的分析和研究，其根本原因是缺乏长尺度和取芯的良好钻孔及高分辨率过孔浅地层剖面的记录。

本文通过高分辨率浅地层剖面资料的分析解译，结合钻孔沉积物的岩性、微体古生物及测年数据，对辽东湾晚第四纪以来的地层格架进行了划分，揭示了辽东湾地层层序，该成果对渤海地区晚第四纪沉积环境的研究具有重要的指示意义。

1 材料与方法

青岛海洋地质研究所于2014 年6—8 月在辽东湾海域采集了2 797 km 浅地层剖面资料（图1）。海上作业采用美国Trimble 公司制造的DSM132 亚米级DGPS 导航定位系统，定位精度的误差小于1 m。浅地层剖面测量使用英国AAE 公司制造的SBP/AAE 浅地层剖面系统，采用的震源激发能量为400～500 J，触发间隔为800 ms，滤波为500～5 000 Hz，作业期间，施工船速控制在5 节以下。本文厚度计算时，声学平均速度统一采用1 600 cm/s。

本研究所分析的辽东湾YKC-2 孔（40.5082°N、121.9074°E，岩芯长70.1 m，水深15.2 m）（图1）是由青岛海洋地质研究所于2016 年8 月利用“勘407”轮所采集。在室内对岩芯进行了详细描述和精细分样。

选取不同层位的岩芯样品173 个（大致40 cm间隔取样）进行微体古生物鉴定。微体古生物鉴定在山东理工大学完成，样品充分浸泡后过250 目（0.063 mm）的标准筛冲洗并烘干，鉴定时过120 目（0.125 mm）标准筛，显微镜下对粗样（＞0.125 mm）中的有孔虫和介形虫进行统计，对于微体化石丰度较大的样品按二分法进行缩分，每样统计底栖有孔虫不少于100 枚、介形类不少于50 枚，丰度较低的样品则统计全样。由于两钻孔中介形虫含量较低，仅在个别层位出现，个数太少而不具有统计意义，本研究仅对底栖有孔虫的数据资料进行了分析。

图1 辽东湾浅地层剖面测线与钻孔位置图Fig.1 Locations of shallow seismic profiles and cores in the Liaodong Bay

表 1 YKC-2 孔光释光年龄Table 1 OSL ages of borehole YKC-2

此外，18 个沉积物样品在中国科学院盐湖研究所完成光释光测年（表1）。样品首先过300 μm 干筛去除粗颗粒；剩余部分依次浸泡10% 稀盐酸和30%双氧水，以去除碳酸盐和有机质；然后再过湿筛提取38～63 μm 组分，将该组分样品浸泡30%氟硅酸约3 周，去除长石；最后滴入少量10%的稀盐酸去除样品与氟硅酸反应过程中产生的氟化物沉淀。在直径为0.97 cm 的不锈钢圆片上均匀涂一层硅胶，将样品均匀黏附在直径约为0.67 cm 的范围内。采用的仪器型号为Risø TL/OSL-DA-20 光释光测年仪，主要通过测量样品中石英介质中放射性物质（U、Th、K 等）等效剂量（De）和年剂量率（Dose rate）来获取[17]。

2 声学地层单元

研究区自下而上，浅地层剖面上识别出了可连续或不连续追踪的6 个声学地层界面（R5、R4、R3、R2、R1、R0），划定了介于上述界面之间具有层序意义的6 个声学地层单元（SU5、SU4、SU3、SU2、SU12、SU11）（图2—4）。其中R3 及其以上的界面全区可连续追踪识别，R3 以下的界面由于二次反射的影响，绝大多数剖面不能连续追踪。

2.1 界面R5 以下的声学地层单元

SU5声学地层单元是研究区识别出的最下部的地层单元，其内部反射特征表现为强振幅、连续性较好、产状近似平行的反射层，局部显示杂乱反射结构。该地层单元的顶界面R5 为一层序界面，该界面横向呈“U”或“V”型结构，对下伏地层存在明显的侵蚀作用。SU5声学单元底界面及下地层单元由于多次波干扰未能识别，其顶界面R5 在研究区西南部及中部海域埋深较大，普遍超过54 m，而在北部、西部及东部埋深明显变浅，大部分埋深小于48 m（图5）。

2.2 界面R5 与R3 之间的声学地层单元

图2 辽东湾NW-SE 向典型浅地层剖面 (Z17 线,剖面位置见图1)Fig.2 NW-SE seismic profile of Line Z17 in Liaodong Bay （see Fig. 1 for location）

图3 辽东湾NW-SE 向典型浅地层剖面 (Z25 线,剖面位置见图1)Fig.3 SE-NW seismic profile of Line Z25 in Liaodong Bay （see Fig. 1 for location）

图4 辽东湾NW-SE 向过LD01 孔[14]典型浅地层剖面 (L8 线,剖面位置见图1)Fig.4 NE-SW seismic profile of Line L8 through core LD01 in Liaodong Bay （see location in Fig. 1）

侵蚀面R5 与R3 之间自下而上共识别出SU4与SU3两个声学地层单元，由于研究区水深普遍较浅，这两个地层单元不同程度地受到二次反射的影响。

SU4声学地层单元内部显示振幅弱的杂乱反射结构，局部似层状反射特征。其底界面为起伏较大的层序界面R5，顶界面为亚水平的海进侵蚀面R4。该地层单元厚度变化较大，大部分厚度为4～9 m，研究区中部厚度较大，最大厚度超过12 m，东部也存在该地层单元的一个沉积中心，厚度普遍超过6 m；研究区西南部该地层单元厚度较小，普遍小于4 m（图5）。

SU3声学地层单元表现为强振幅、连续性较好、近似水平或缓波状的反射层。其底界面为侵蚀面R4，对下伏地层表现为明显的侵蚀削截，顶界面R3 为一起伏不平、局部“U”或“V”型结构的侵蚀间断面，对下伏的SU4及其下的地层单元产生明显的侵蚀作用。该地层单元的厚度南北差异明显，呈北厚南薄的分布特征，其中北部厚度大部分超过8 m，最厚处超过16 m，位于研究区西北部，南部厚度普遍小于4 m（图6）。

2.3 界面R3 以上的声学地层单元

SU2声学地层单元内部反射表现为弱振幅、杂乱，或较强振幅、弱水平成层等反射特征。其底界面为波状起伏的层序界面R3，顶界面R2 对其存在明显的侵蚀削截。研究区SU2声学地层单元区域厚度变化大，下切深度和残存状态的纵、横向变化较大，大部分厚度为8～24 m，最大厚度可达30 m。该地层单元厚度自北向南逐渐增加，与SU3单元分布特征相反，北部厚度普遍小于10 m，至南部大部分厚度超过16 m（图6）。

图5 辽东湾SU4 单元厚度(A)与SU5 单元顶界面埋深(B)图Fig.5 Thickness of unit SU4（A）and depth of unit SU5 top boundary（B）in Liaodong Bay

图6 辽东湾SU2(A)与SU3(B)单元厚度图Fig.6 Isopach maps of units SU2（A）and SU3（B）in Liaodong Bay

SU12声学地层单元内部呈现强振幅、高频率、近水平的反射特征，对下伏地层产生明显的削蚀作用。该地层单元的顶底界面近似平行，其中底界面R2 为一海进侵蚀面，顶界面R1 为最大海泛面。SU12地层单元厚度变化较小，大部分地层厚度为4～6 m，自西向东地层厚度逐渐增加，其中西部地层厚度大部分小于4 m，东部地层厚度普遍超过5 m，最大厚度超过8 m，位于研究区东南角（图7）。

SU11声学地层单元内部反射振幅较弱，多呈现声学半透明或透明反射层，呈平缓垂向加积或低倾角向海进积，下超于底部侵蚀界面R1。该地层单元沉积厚度变化较大，自东北至西南厚度逐渐变薄，其中研究区西南部厚度普遍小于1 m，局部缺失；东北部该地层单元厚度大部分为2～7 m，自岸向海逐渐变薄，其中现行大辽河、辽河河口为两沉积中心，厚度普遍超过6 m，大辽河河口外地层最大厚度可达10 m（图7）。

3 声学地层与沉积地层的对比

前人对辽东湾及周边海域的地层、沉积环境等做了大量的研究工作，涉及的钻孔以辽东湾海域的LD01 孔[14]与LDD7 孔[10]、渤海东部的BC-1 孔[1]等比较典型。基于钻孔沉积特征、微体古生物变化及光释光年代数据，辽东湾采集的YKC-2 孔的沉积地层自下而上可划分为6 个地层单元或亚单元（DU5、DU4、DU3、DU2、DU12、DU11），而且，依据钻孔沉积地层的划分与浅地层剖面声学地层的划分存在良好的对应关系。需要特别说明的是，钻孔位置处晚第四纪以来的地层单元均存在不同程度的缺失，这主要由于低海面时期河流下切侵蚀与地层发育情况导致。以下对YKC-2 孔沉积地层与过孔剖面进行对比。

（1）70.10～39.60 m（DU5）：对应声学地层单元SU5地层单元，依据岩性及微体古生物组合的变化可 分 为 两 段：70.10～48.10 m 与48.10～39.60 m 段（图8）。

70.10 ～48.10 m 段：岩性主要为灰色中砂，局部夹黏土质条带，微体古生物鉴定未见底栖有孔虫与介 形 虫。55.27 与61.21 m 两 处OSL 年 龄 分 别为93.7±9.4 kaBP 与94.4±9.4 kaBP，显示这一时期对应MIS5 期的沉积。YKC-2 孔位于晚第四纪以来仍在活动的郯庐断裂的边缘地带[18]，推测这一地质时期构造作用为该钻孔位置沉积环境演化的主控因素，导致这一海平面较高时期的海相性沉积较弱。

图7 辽东湾SU11(A)与SU12(B)单元厚度图Fig.7 Isopach maps of units SU11（A）and SU12（B）in Liaodong Bay

图8 辽东湾YKC-2 孔沉积地层与过孔浅地层剖面对比 (剖面位置见图1)Fig.8 Correlation between the depositional units of the core YKC-2 and the seismic units from the profile through the core in Liaodong Bay （see Fig. 1 for location）

48.10 ～39.60 m 段：岩性主要为灰色粉砂质黏土，局部夹灰色中细砂薄层，见贝壳碎片。微体古生物主要为底栖有孔虫，丰度、简单分异度与复杂分异度均较低，优势种以广温广盐浅水种毕克卷转虫Ammonia becarii（Linné）var.与冷水种具瘤先希望虫Protelphidium turberculatum为主。毕克卷转虫Ammonia becarii（Linné）var.是我国东部陆架海区最为常见的滨岸浅水种，广见于潮上带、潮间带、瀉湖与河口湾等半咸水及正常盐度的滨岸浅海，具瘤先希望虫Protelphidium turberculatum是我国渤、黄、东海常见的低温低盐环境代表种，主要集中在沿岸流影响区以及冷水团控制的范围以内[19-20]。40.87、44.70 与45.30 m 3 处OSL 年龄分别为70.2±7.0 kaBP、70.6±7.1 kaBP 与69.5±7.0 kaBP，总体显示该阶段对应MIS5 期晚期滨海相沉积。

（2）39.60～30.40 m（DU4）：对应声学地层单元SU4，岩性下部39.60～30.84 m 段为灰色、灰褐色中砂，局部夹粉砂质薄层或有机质条带，上部30.84～30.40 m 段为深灰色黏土夹细砂薄层或透镜体。该层位微体古生物鉴定未见底栖有孔虫与介形虫。结合下伏地层DU5测年数据，推测该层段为MIS4期低海面时期的河流相沉积及之后海平面上升过程中的河道充填相沉积。

（3）30.40～5.40 m（DU2）：对应声学地层单元SU2，岩性主要为灰色、灰黄色细砂—中砂，局部夹粉砂薄层或有机质条带，见植物碎屑。微体古生物鉴定未见底栖有孔虫与介形虫。该段上部5.27 与5.86 m 两处OSL 年龄分别为12.7±1.3 kaBP 与13.4±1.3 kaBP，总体显示该阶段对应末次冰盛期低海面时期的河流相沉积与之后海平面上升阶段中的河道充填相沉积。

（4）5.40～1.45 m（DU12）：对应声学地层单元SU11，岩性主要为灰色细砂，见贝壳碎片，局部夹黏土薄层。该段底栖有孔虫丰度、简单分异度与复合分异度较高，优势种以冷水种Buccella frigida与Protelphidium turberculatum、浅水种Elphidium advenum、Elphidium magellanicum 为 主。1.96 m 处 的OSL 年龄为4.5±0.5 kaBP，总体显示该阶段对应全新世早中期的滨海相沉积。

（5）1.45～0 m（DU11）：对应声学地层单元SU11，岩性主要为灰色粉砂，含水量高。该段底栖有孔虫丰度、简单分异度与复合分异度达到整个钻孔最高，优势种以冷水种Buccella frigida与Protelphidium turberculatum为主。结合下覆地层测年数据，总体显示该阶段对应全新世高海面以来的浅海相沉积。

4 辽东湾晚第四纪以来的层序地层分析

通过调查区浅地层剖面声学地层与YKC-2 孔沉积地层的综合对比分析，结合该钻孔年代学资料及前人有关辽东湾及周边海域的研究成果，揭示了MIS5 期 以 来 研 究 区6 个 地 层 单 元（DU11、DU12、DU2、DU3、DU4、DU5）和介于其间的6 个地层界面（R0、R1、R2、R3、R4、R5）。R3 与R5 两个层序界面将这些地层单元自上而下划分为3 个层序（SQ1、SQ2、SQ3）（图9）。

R3 层序界面介于层序SQ1 与SQ2 之间，为研究区晚第四纪以来年代最新的层序界面。由于缺乏直接的证据，R3 层序界面的精准年龄很难识别，仅根据上下地层推测。辽东湾YKC-2 孔与周边典型钻孔（Bc-1 孔[1]与LDD7 孔[10]）多项指标表明（图10），R3 层序界面对应末次冰盛期低海面时期，与这一地质时期的陆架裸露及河流下切作用密切相关，发育明显的沉积间断。R5 层序界面介于层序SQ2 与SQ3 之间，其界面特征与R3 层序界面相似，与MIS4 期低海面时间的陆架裸露及河流侵蚀相关（图10），下切下部地层，造成层序SQ2 与SQ3 之间明显的侵蚀间断。

强制海退体系域是海平面下降过程中形成的，在声学反射中通常表现为进积沉积[23]。辽东湾晚第四纪以来地层中的强制海退体系域的沉积，基本被随后的MIS4 期与末次冰盛期低海面期间的古河流作用侵蚀殆尽，很难在研究区的地层中识别。因而，本文仅从低水位、海侵与高水位体系域的角度对辽东湾的地层层序进行了分析。

SQ3 位于层序界面R5 之下，是辽东湾晚第四纪以来可识别出的最下部的层序，主要涉及DU5一个地层单元，其下部的单元由于钻孔或浅地层剖面穿透深度的限制未能识别。该地层单元对应MIS5期海平面相对较高时期的滨浅海沉积，代表海侵体系域与高水位体系域，然而，区分两体系域之间的最大海泛面在研究区的浅地层剖面与钻孔地层中均未识别出，主要受限于浅地层剖面分辨率与钻孔缺乏精确定年数据，同时，可能由于MIS5 期海平面波动过程中的再造与混合作用造成。

SQ2 层组介于层序界面R5 与R3 之间，由DU4与 DU3上下两个地层单元组成。DU4地层单元上覆于层序界面R5 之上，对应MIS4 期海平面较低时期的河流与河道充填沉积，代表低水位体系域,该地层 单 元 在LDD7 孔50.10～26.80 m 段 也 有 显 示[10]（图10）；DU3地层单元上覆于海进侵蚀面R4 之上，对应MIS3 期早中期滨海沉积，代表海侵体系域，YKC-2 孔该期地层由于末次冰盛期低海面的河流下切侵蚀殆尽，LDD7 孔对应层位埋深26.80～11.90 m，海侵持续时间为45.84～27.66 cal.kaBP[10]（图10）。该层序无高水位体系域，主要由于后期末次冰盛期低海面时期的研究区裸露于河流下切侵蚀殆尽导致。

图9 辽东湾晚第四纪层序地层与海平面变化(Chappell et al[21],Lea et al[22])的对比Fig.9 Sequence stratigraphy in Liaodong Bay and its correlation with sea-level fluctuation（Chappell et al[21],Lea et al[22]） during Late Quaternary

图10 辽东湾YKC-2 孔与周边Bc-1 孔[1]、LDD7 孔[10]地层对比(钻孔位置见图1)Fig.10 Stratigraphic correlation of core YKC-2 with cores Bc-1[1] and LDD7[10] （see Fig. 1 for location）

SQ1 上覆于层序界面R3 之上，是研究区晚第四纪以来最上部的层序，包括DU2、DU12、DU113 个地层单元。DU2位于层序界面R3 之上，对应末次冰盛期低海面及之后海平面上升期间的河流与河道充填沉积，代表低水位体系域；DU12地层单元位于海进侵蚀面R2 之上，对应全新世早中期滨海沉积，代表海侵体系域；DU11地层单元上覆于最大海泛面R1 之上，对应全新世高海面以来的浅海沉积，代表高水位体系域。全新世海侵开始之后，渤海范围内不同钻孔记录中陆续出现代表海侵体系域与高水位体系域的海侵层，其中Bc-1 孔对应8.60～0 m 段[1]，LDD7 孔对应1.55～0 m 段（10.64 cal.kaBP以来）[10]（图10）。

研究区的海侵体系域，受辽东湾平缓地形的影响，主要表现为厚度较薄且变化较小、分布广泛的特征。然而，低水位体系域则表现为厚度与横向分布均变化较大，这主要归因于低海面时期纵横交错分布的高能河流的侵蚀作用，类似的分布模式在黄海[24-25]、东海[26]均有过报道。

5 结论

（1）辽东湾晚第四纪以来的地层自下而上划定了6 个声学地层界面和介于上述界面之间具有层序意义的6 个声学地层单元（SU5、SU4、SU3、SU2、SU12、SU11），而且，研究区声学地层的划分与钻孔6 个沉积地层单元（DU5、DU4、DU3、DU2、DU12、DU11）的划分具有良好的对应关系。

（2）辽东湾晚第四纪MIS5 期以来的地层层序自下而上被层序界面R5 与R3 划分为3 个层序（SQ3、SQ2、SQ1），其中SQ3 位于层序界面R5 之上，仅识别出该层序上部的MIS5 期海平面相对较高时期形成的DU5地层单元，代表海侵体系域与高水位体系域；SQ2 介于层序界面R5 与R3 之间，自下而上由低水位体系域（DU4地层单元，对应MIS4 期海平面较低时期的河流与河道充填沉积）与海侵体系域（DU3地层单元，对应MIS3 期早中期滨海沉积）组成。SQ1 上覆于层序界面R3，自下而上包括低水位体系域（DU2地层单元，对应末次冰盛期低海面及之后海平面上升期间的河流与河道充填沉积）、海侵体系域（DU12地层单元，全新世早中期滨海沉积）和高水位体系域（DU11地层单元，全新世高海面以来的浅海沉积）。

（3）研究区的海侵体系域，厚度较薄且变化较小，分布广泛，而低水位体系域的厚度与横向分布均变化较大。