辽东湾是一个西、北、东三面被陆地包围的半封闭海湾，南面与渤海相通，水深普遍小于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。

图  1  辽东湾浅地层剖面测线与钻孔位置图

Figure  1.  Locations of shallow seismic profiles and cores in the Liaodong Bay

下载: 全尺寸图片 幻灯片

本研究所分析的辽东湾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枚，丰度较低的样品则统计全样。由于两钻孔中介形虫含量较低，仅在个别层位出现，个数太少而不具有统计意义，本研究仅对底栖有孔虫的数据资料进行了分析。

此外，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]。

表  1  YKC-2孔光释光年龄

Table  1.  OSL ages of borehole YKC-2

样品深度/m	OSL年龄/kaBP
1.96	4.5±0.5*
2.90	14.1±1.4
4.01	15.3±1.5
5.27	12.7±1.3*
5.86	13.4±1.3*
7.00	30.3±3.0
18.47	59.4±5.9
37.15	54.4±5.4*
40.66	76.8±7.7
40.87	70.2±7.0*
44.70	70.6±7.1*
45.30	69.5±7.0*
49.89	＞93.9*
55.27	93.7±9.4*
60.37	91.2±9.1*
61.21	94.4±9.4*
68.05	85.4±8.5*
69.20	＞98.6*
注：*年龄为本文采用年龄。

下载: 导出CSV 

| 显示表格

2.   声学地层单元

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

图  2  辽东湾NW-SE向典型浅地层剖面 （Z17线，剖面位置见图1）

Figure  2.  NW-SE seismic profile of Line Z17 in Liaodong Bay （see Fig. 1 for location）

下载: 全尺寸图片 幻灯片

图  3  辽东湾NW-SE向典型浅地层剖面 （Z25线，剖面位置见图1）

Figure  3.  SE-NW seismic profile of Line Z25 in Liaodong Bay （see Fig. 1 for location）

下载: 全尺寸图片 幻灯片

图  4  辽东湾NW-SE向过LD01孔^[14]典型浅地层剖面 （L8线，剖面位置见图1）

Figure  4.  NE-SW seismic profile of Line L8 through core LD01 in Liaodong Bay （see location in Fig. 1）

下载: 全尺寸图片 幻灯片

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

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

图  5  辽东湾SU₄单元厚度（A）与SU₅单元顶界面埋深（B）图

Figure  5.  Thickness of unit SU₄（A）and depth of unit SU₅ top boundary（B）in Liaodong Bay

下载: 全尺寸图片 幻灯片

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

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

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

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

图  6  辽东湾SU₂（A）与SU₃（B）单元厚度图

Figure  6.  Isopach maps of units SU₂（A）and SU₃（B）in Liaodong Bay

下载: 全尺寸图片 幻灯片

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

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

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

图  7  辽东湾SU₁¹（A）与SU₁²（B）单元厚度图

Figure  7.  Isopach maps of units SU₁¹（A）and SU₁²（B）in Liaodong Bay

下载: 全尺寸图片 幻灯片

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

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

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

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

图  8  辽东湾YKC-2孔沉积地层与过孔浅地层剖面对比 （剖面位置见图1）

Figure  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）

下载: 全尺寸图片 幻灯片

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]，推测这一地质时期构造作用为该钻孔位置沉积环境演化的主控因素，导致这一海平面较高时期的海相性沉积较弱。

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

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

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

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

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

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

图  9  辽东湾晚第四纪层序地层与海平面变化（Chappell et al^[21]，Lea et al^[22]）的对比

Figure  9.  Sequence stratigraphy in Liaodong Bay and its correlation with sea-level fluctuation （Chappell et al^[21],Lea et al^[22]） during Late Quaternary

下载: 全尺寸图片 幻灯片

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

图  10  辽东湾YKC-2孔与周边Bc-1孔^[1]、LDD7孔^[10]地层对比（钻孔位置见图1）

Figure  10.  Stratigraphic correlation of core YKC-2 with cores Bc-1^[1] and LDD7^[10] （see Fig. 1 for location）

下载: 全尺寸图片 幻灯片

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

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

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

SQ1上覆于层序界面R3之上，是研究区晚第四纪以来最上部的层序，包括DU₂、DU₁²、DU₁¹ 3个地层单元。DU₂位于层序界面R3之上，对应末次冰盛期低海面及之后海平面上升期间的河流与河道充填沉积，代表低水位体系域；DU₁²地层单元位于海进侵蚀面R2之上，对应全新世早中期滨海沉积，代表海侵体系域；DU₁¹地层单元上覆于最大海泛面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个声学地层单元（SU₅、SU₄、SU₃、SU₂、SU₁²、SU₁¹），而且，研究区声学地层的划分与钻孔6个沉积地层单元（DU₅、DU₄、DU₃、DU₂、DU₁²、DU₁¹）的划分具有良好的对应关系。

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

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