Study on sequence stratigraphy and sedimentary systems of the Wenchang Formation in the southern Xijiang depression of the Pearl River Mouth Basin
-
摘要: 珠江口盆地西江凹陷南部包括两个相邻洼陷：西江36洼和番禺4洼，钻探结果却揭示了截然不同的油气勘探前景，其中，番禺4洼探明储量已过亿吨，而西江36洼却未有商业发现，导致勘探上对西江36洼烃源潜力有所顾虑。为了揭示这种差异油气地质条件产生的原因，亟待开展两个洼陷的对比研究。以主要烃源层系文昌组为切入点，通过开展精细的层序-沉积研究，进一步系统评估西江36洼烃源潜力。依据地震、钻井及分析化验等资料，在西江36洼与番禺4洼文昌组共识别了6个三级层序。西江36洼和番禺4洼文昌组发育扇三角洲、辫状河三角洲、湖底扇及湖泊沉积等沉积相类型，垂向上沉积相带具有旋回特征。伴随裂陷不同阶段和边界断层活动差异，陡坡带和缓坡带三角洲沉积体系规模呈现差异演化特征，而湖盆和半深湖-深湖相规模总体呈现出先增大后减小的规律。沉积相带在横向上也具有明显的迁移特征，表现为文六段时期先在西江36洼开始沉积，文五段开始再扩展到番禺4洼，同时文一段时期西江36洼湖盆已经消失，只在番禺4洼沉积充填。进一步对比两个洼陷各三级层序反映优质烃源岩条件的半深湖-深湖规模及物源供给量，结果显示在文六段至文四段时期两个洼陷烃源条件基本一致。鉴于该层段作为番禺4洼的主力源岩，其巨大的生烃潜力已经被勘探所证实，因此认为不能低估具有相似烃源条件的西江36洼生烃潜力。Abstract: The Xijiang depression in the Pearl River Mouth Basin consists of two adjacent sags, the Panyu 4 sag and the Xijiang 36 sag, in which quite different oil and gas exploration prospects have been revealed by drilling. The proved reserves of Panyu 4 sag have exceeded 100 million tons, while no commercial discoveries been found in the Xijiang 36 sag so far. In order to reveal the causes of such a difference in oil and gas geological conditions, it is necessary to carry out a comparative studies of the two sags. In this paper, the Wenchang Formation, the main source series, is selected as the key point for breakthrough, and the hydrocarbon source potential of Xijiang 36 sag is systematically evaluated through detailed sequence stratigraphic studies. Based on the data of seismic, drilling and laboratory analysis, six third-order sequences have been identified in the Wenchang Formation in both of the Xijiang 36 and Panyu 4 sags. There are fan delta, braided river delta, sublacustrine fan and lacustrine sedimentary facies in the Wenchang Formation of Xijiang 36 and Panyu 4 sags, and the sedimentary facies belts also show a cyclic patterns vertically. In the different stages of rifting and/or boundary fault activities, the scale of delta depositional system in steep slope zone and gentle slope zone shows different evolutionary features, whereas the scales of Lake Basin and its deep part in the middle lake expanded in the beginning and shrank later on. The sedimentary facies belt also has obvious migration in lateral direction. In the Wen-6 Period, the deposition was initiated in the Xijiang 36 sag, and then extended to Panyu 4 sag with time. In the period of Wen-1, the lake basin occurred only in the Panyu 4 sag and disappeared in the Xijiang 36 sag. The distribution pattern of semi deep and deep lake deposits and study of sediment sources suggest that the high-quality hydrocarbon source potential values of each third-order sequence in the two sags are comparable and the source rock conditions are basically the same in the two sags in the wen-6 and wen-4 periods. As the main source rock of the Panyu 4 sag has been confirmed by exploration, therefore, it is considered that the exploration potential of Xijiang 36 sag with similar hydrocarbon source conditions should not be underestimated too early.
-
20世纪90年代以来，下切河谷（Incised valley）及其沉积充填层序的研究颇受关注[1-3]，这是由于其是油气资源的勘探目标、蕴含丰富的地层演化和古环境变化信息。晚第四纪海平面经历了大幅升降变化，在全球海岸带和陆架发育了多期规模较大的下切河谷，其充填地层保存相对完整，已经成为层序地层学、古气候-海平面变化研究的热点[4-7]。
近三十多年来长江三角洲地区晚第四纪下切河谷的形成与演化得到了学术界的重视[8-10]。晚第四纪以来长江三角洲地区的下切河谷可分为三期（早期、中期和晚期），其分布范围有所差别[11]（图1），推测其沉积充填层序的形成时间分别对应于125～60 ka（日历年龄）、60～25 ka和25 ka以来的时期[12]。然而，对于晚第四纪早期和中期下切谷充填层序年龄的研究尚不多见[13-14]。
![]()
图 1 长江三角洲和晚第四纪下切河谷分布图[11]以及本文研究的JC-1205孔位置I. 三角洲主体，II. 三角洲南翼，III. 三角洲北翼；高程为相对于现今海平面的高程；晚第四纪早期下切河谷的北界尚未确定，未标出。Figure 1. Schematic map showing the distribution of the Yangtze River delta and locations of late Quaternary incised valleys [11], as well as the location of core JC-1205 investigated in this studyThe main body (I), southern (II) and northern (III) flanks of the delta are also shown.2012年我们在实施江苏海岸带地质调查时在江苏南通近岸区施工了一个全取芯钻孔，该孔钻取了晚第四纪地层。本文基于该孔岩芯的综合测试分析，确定该地层的沉积相变化特征和年龄，并在此基础上探讨晚第四纪下切谷及其沉积充填序列的形成与海平面变化的关系。
1. 研究区地质概况
长江三角洲是中全新世以来长江泥沙堆积形成的以镇江-扬州为顶点、杭州湾北岸为南界、大致以弶港辐射沙洲的陆上延续部分为北界的三角洲沉积体系[11, 15]（图1），在地理上包括江苏省南部、上海市和浙江省杭嘉湖地区。该地区为地势平缓的广阔平原，海拔高度一般为3～5 m，总趋势为西高东低。三角洲可分为主体和南、北两翼[11]：主体部分大致与末次冰期（晚第四纪晚期）古河谷位置相当，以存在不同时期雁行排列的河口坝为特征，曾为长江口摆动的地带；两翼部分大致与末次冰期古河间地分布区相当。
长江三角洲地处扬子陆块下扬子地块东段，自元古代以来区域经历了复杂的构造运动，形成了丰富多样的构造样式，其中燕山期火山活动较为强烈，是中国东部火山岩浆活动带的重要组成部分。新近纪以来，由于喜马拉雅运动的影响，中国东部盆地强烈沉降。研究表明，长江三角洲地区的松散沉积物之下是扬子板块元古宙变质岩系与古生代、中生代和新生代岩浆岩和沉积岩；基岩之上的松散沉积物总体上西薄东厚，主要为第四纪堆积物（厚度50～450 m），下部见新近纪沉积物（厚度20～100 m，局部缺失）[16-17]。自上新世以来长江三角洲表现为“间歇式”的沉降，沉降速率主要为1～3 mm/a [17]。新近纪以来该区经历了山间盆地冲积环境（新近纪—早更新世）、辫状河和曲流河环境（早更新世晚期—中更新世）、滨海-浅海-三角洲环境（晚更新世以来）[18]。
晚第四纪全球海平面在冰期-间冰期气候旋回下经历了大幅波动，在长江三角洲地区发育了3期下切河谷[11-12]（图1），不同期次的河谷在分布范围上呈现部分叠置的现象，造成早期和中期的河谷充填沉积物被部分侵蚀[12]。末次盛冰期，海平面下降幅度超过120 m[19]，现今的长江三角洲地区发育古长江的下切河谷及其南北两侧的古河间地（图1），而该下切谷一直延伸至东海外陆架[20-21]。冰后期海平面快速上升，下切河谷逐渐被海侵成为河口湾；至8.0 ka，研究区形成以镇江-扬州为顶点的巨型古河口湾；此后海平面上升趋缓并在约7 ka达到最高水位后保持相对稳定，长江三角洲自约8.0 ka开始发育并逐渐向海推进直至岸线达到现今位置[15]。
2. 材料与方法
2012年9月，青岛海洋地质研究所在江苏省南通市如东县东凌港镇近岸区完成一口全取芯钻孔JC-1205孔的施工，钻孔坐标32°10′47.436″N、121°24′40.207″E（图1）。在去掉钻孔最顶部的人工填土之后，该孔的高程是+1.30 m，进尺79 m。钻探施工使用GXY-1A型岩芯钻机，采用液压回转钻进工艺，取芯管直径为110 mm，所取样品直径为96 mm。钻孔泥质沉积物的取芯率大于85%，砂质沉积物的取芯率大于60%。
在室内首先对岩芯对半分开，然后进行沉积特征描述、拍照并采取沉积物粒度、有孔虫鉴定、AMS 14C测年和光释光（OSL）测年等实验样品。按0.10～0.25 m间距，采取粒度样品596个；按0.10～0.20 m间距，采取有孔虫鉴定样品567个；在0～68 m岩芯中，挑选了13个腹足类和双壳类的壳体样品和1个有孔虫样品开展14C测年；在22～79 m岩芯中，挑选11了个样品开展OSL测年。
沉积物样品的粒度分析在青岛海洋地质研究所实验检测中心完成。对于含砾石的较粗粒样品，采用筛析法计算其粒度组分。对于不含砾石组分的样品，选取沉积物样品10～20 g，经双氧水和稀盐酸浸泡处理，除掉有机质和碳酸盐，然后洗盐，用六偏磷酸钠溶液经超声波分散后，用英国马尔文（MALVERN）公司生产的Mastersizer-2000型激光粒度分析仪进行粒度测试。仪器测量范围为0.02～2000 μm，粒级分辨率为0.01 Φ，重复测量的相对误差＜2%。沉积物命名和粒度参数分别采用文献[22]和[23]提供的方案。
底栖有孔虫样品的鉴定：称50 g干样加入清水与双氧水（H2O2）浸泡分散，经用250目（Φ=0.063 mm）分析筛淘洗，筛上颗粒烘干后，视样品量及有孔虫丰度将其缩分至1/2、1/4、1/8……后，在双目实体显微镜下挑选标本，进行鉴定与数量统计。一般统计个体数在100枚以上，不足100枚的做全部样本统计。底栖有孔虫的丰度表示每50 g干样中有孔虫个体数，简单分异度表示所有属种的种数，重要属种的百分含量按其个体数占底栖有孔虫全群个体总数的百分含量来表示。
样品的AMS 14C测年由美国Beta分析有限公司承担，主要测年物质包括：有孔虫壳体、双壳类和腹足类壳体等。直接测定年龄是以5568年为半衰期，同时测量样品的δ 13C值，并根据分馏效应进行校正，即获得惯用年龄。日历年龄是惯用年龄经过校正程序Calib Rev 8[24]校正所得（表1）。CALIB Rev 8没有提供研究区的海洋储库校正值（△R），对于海洋生物样品日历年龄的校正，本文根据CALIB Rev 8提供的3个黄海地区（包括1个青岛近岸海区和2个韩国西南海岸区）和6个台湾北部近岸样品的△R值，利用CALIB Rev 8计算得到△R值为−128±104 yr。本文所用的14C日历年龄都是从公元1950向前起算的，以“cal yr BP”表示，或用“ka”来代替“cal kyr BP”。本文或前人没有经过日历年龄校正的年龄以14C年龄表示（14C yr BP）。
表 1 JC-1205孔 AMS 14C测年结果Table 1. AMS 14C dating results for core JC-1205深度/m 测年材料 δ 13C/‰ 惯用年龄
/14C yr BP日历年龄/cal yr BP Beta编号 中值 范围 (1σ) 0.29 腹足类 (Semisulcospira gredleri) −2.7 102.3 ± 0.3 − − 371210 0.55 腹足类 (Semisulcospira gredleri) −3 100.1 ± 0.4 − − 371211 13.01 腹足类 −0.6 1010 ± 30 BP 561 463～663 468133 21.42 底栖有孔虫 −2.1 3790 ± 30 BP 3724 3562～3876 468132 60.25 腹足类
(Angulyagra sp.)−2.6 ＞43500 BP − − 371221 64.15 腹足类
(Angulyagra sp.)−5.2 39110 ± 550 BP 42703 42432～42917 371212 64.71 双壳类和腹足类碎片 −8.6 ＞43500 BP − − 371213 65.56 腹足类
(Angulyagra sp.)−8.6 37770 ± 470 BP 42165 41983～42377 371214 65.8 腹足类
(Angulyagra sp.)−5.8 ＞43500 BP − − 371215 65.96 腹足类
(Angulyaga sp.)−4.7 ＞43500 BP − − 371216 66.51 腹足类
(Angulyaga sp.)−9.5 ＞43500 BP − − 371217 67.05 腹足类
(Angulyaga sp.)−7.6 39350 ± 500 BP 42805 42519～43016 371218 67.47 腹足类
(Angulyaga sp.)−5.6 ＞42000 BP − − 371219 67.63 腹足类
(Angulyaga sp.)−7.8 ＞43500 BP − − 371220 样品的OSL测年是在中国科学院青海盐湖研究所光释光实验室完成（表2），样品的前处理和测试方法详见过去的文献[25-26]。从样品中提取38～63 μm的石英颗粒来进行测试，测试仪器为Risø TL/OSL DA-20，等效剂量（De）测试采用单片再生剂量法与标准生长曲线法相结合的方法[27-28]。每个样品共获得18个等效剂量值，采用最小年代模型（MAM）计算样品的最终等效剂量。
表 2 JC-1205孔光释光测年结果Table 2. OSL dating results for core JC-1205孔深/m 粒径/μm U/10-6 Th/10-6 K/10-6 含水率/% 离散值 等效剂量/Gy 剂量率/(Gy/ka) 年代/ka 22.74 38～63 1.36 ± 0.2 8.06 ± 0.5 1.72 ± 0.03 23 ± 5 0.08 ± 0.01 2.53 ± 0.08 2.13 ± 0.09 1.2 ± 0.1 24.4 38～63 1.12 ± 0.2 6.84 ± 0.4 1.34 ± 0.03 18 ± 5 0.09 ± 0.02 7.15 ± 0.25 1.79 ± 0.08 4.0 ± 0.2 31.37 38～63 1.82 ± 0.3 10.04 ± 0.6 1.47 ± 0.03 21 ± 5 0.16 ± 0.03 133.1 ± 10.7 2.16 ± 0.10 61.5 ± 5.7 35.74 38～63 1.71 ± 0.3 9.09 ± 0.5 1.46 ± 0.03 23 ± 5 0.15 ± 0.03 131.3 ± 10.2 2.04 ± 0.09 64.5 ± 5.9 37.56 38～63 1.34 ± 0.2 8.08 ± 0.5 1.41 ± 0.03 14 ± 5 0.13 ± 0.02 138.2 ± 7.3 2.04 ± 0.09 67.7 ± 4.8 41.67 38～63 1.44 ± 0.2 8.07 ± 0.5 1.58 ± 0.03 15 ± 5 0.24 ± 0.04 116.1 ± 15.3 2.18 ± 0.10 53.2 ± 7.4 47.31 38～63 1.07 ± 0.2 5.83 ± 0.4 1.37 ± 0.03 17 ± 5 0.24 ± 0.04 84.4 ± 12.1 1.75 ± 0.08 48.3 ± 7.3 49.78 38～63 1.04 ± 0.2 5.55 ± 0.4 1.74 ± 0.03 20 ± 5 0.20 ± 0.03 125.3 ± 11.7 1.98 ± 0.09 63.3 ± 6.6 68.52 38～63 2.03 ± 0.3 11.32 ± 0.6 1.61 ± 0.03 21 ± 5 – 339 ± 30 2.38 ± 0.11 143 ± 14* 73.11 38～63 1.87 ± 0.3 10.43 ± 0.6 1.75 ± 0.03 22 ± 5 – 357 ± 31 2.39 ± 0.11 149 ± 15* 78.94 38～63 1.78 ± 0.3 10.56 ± 0.6 1.27 ± 0.03 18 ± 5 – 330 ± 42 2.07 ± 0.10 160 ± 22* 注：*测年结果超出了石英颗粒OSL测年的上限，仅供参考。 3. 结果
3.1 地层序列和沉积相
根据沉积特征和底栖有孔虫分布特征，JC-1205孔岩芯沉积物可以划分为4个沉积单元，从上至下编号为DU 1—DU 4（图2）。
![]()
图 2 JC-1205孔岩芯岩性柱状图附加AMS 14C和OSL测年结果。Figure 2. Lithologic log of core JC-1205, showing calibrated 14C and OSL ages3.1.1 DU 4（79.00～68.32 m）
DU 4主要由深灰色黏土质粉砂夹粉砂质线理-透镜体、粉砂夹黏土质线理、粉砂与黏土质条带密集互层（毫米至厘米级）组成（图3a，3b），生物扰动中等，含少量钙质结核，偶见贝壳碎片及碳质斑点。本单元不含有孔虫（图4），中部（73.5～75.3 m）含较多淡水腹足类 (Angulyaga sp.)。
![]()
图 3 JC-1205孔典型岩芯的照片a. DU 4底部（78.48～78.72 m），深灰色粉砂夹较多黏土质线理（毫米级），生物扰动中等；b. DU 4上部（71.50～71.77 m），深灰色黏土质粉砂夹较多粉砂质线理（毫米级）和透镜体，生物扰动中等；c. DU 3下段下部（65.76～66.01 m），深绿灰-灰黄色细砂质中砂（含粗砂，偶见细砾），分选较差，含少量贝壳碎片；d. DU 3上段中部（33.44～33.68 m），深灰黏土质粉砂与粉砂密集互层（毫米级），生物扰动较弱；e. DU 2中部（27.75～27.96 m），棕色斑杂灰绿色黏土质粉砂，含少量粉砂-砂质斑点，生物扰动较强；f. DU 1下段中部（23.60～27.85 m），深黄灰色含细砂-黏土质粉砂，夹较多深灰色黏土质条带以及较多粉砂-细砂线理、透镜体，生物扰动中等—较弱；g. DU 1中段上部（9.86～10.07 m），深黄灰色黏土质粉砂与深灰色粉砂密集互层（毫米级），潮汐层理发育，生物扰动中等—较弱；h. DU 1上段中部（1.97～2.19 m），深灰黄色黏土质条带夹粉砂-细砂质粉砂薄层（毫米至厘米级），生物扰动强烈。箭头指示岩芯的上方。Figure 3. Photographs of representative sedimentary facies in core JC-1205![]()
图 4 JC-1205孔岩芯中底栖有孔虫丰度（a）、简单分异度（b）和主要属种含量（c-k）垂向分布特征Figure 4. Down-core variations in (a) benthic foraminiferal abundance, (b) simple diversity, and (c–k) relative abundance of the main foraminiferal species in core JC-1205粒度分析结果显示（图5），DU 4沉积物以粉砂为主，黏土和和砂次之，平均粒径主要为5～7 Φ，平均值为6.2；分选系数为1.3～1.9（平均为1.66），分选较差；偏态主要为0.1～0.3（平均0.24），为正偏；峰度主要为0.9～1.1（平均1.03），峰度较扁平。
![]()
图 5 JC-1205孔岩芯沉积物的粒度参数（a-d）和粒度组分含量（e-h）垂向分布特征Figure 5. Down-core variations in grain-size parameters (a-d) and grain-size compositions (e-h) in core JC-1205综上所述，DU 4被解释为河流环境之下的洪泛平原沉积。
3.1.2 DU 3（68.32～29.25 m）
根据岩性特征，DU 3可划分为两段：下段（68.32～35.40 m）和上段（35.40～29.25 m）。DU 3与下伏DU 4之间为侵蚀接触。
下段总体上显示正粒序（图2，5），底部0.64 m（67.69～68.32 m）为灰-深灰色砾石层，含中-粗砂及粉砂质条带，含较多贝壳碎片，砾石直径可达2～3 cm，有磨蚀，以次棱角状为主；下部（56.2～67.69 m）为深灰-深灰黄色细砂质中砂（含粗砂，底部偶见细砾）（图3c），向上渐变为深灰色极细砂-粉砂质细砂（含中砂），顶部0.8 m为深灰色极细砂质粉砂。分选一般—较差，板状/槽状交错层理发育，偶夹少量碳质线理、碳屑和贝壳碎片；下段沉积物不含底栖有孔虫（图4），60.0 m以下局部见较多的淡水腹足类 (Angulyaga sp.)。
下段沉积物的粒度分析结果（图5）显示，下段沉积物以砂为主，粉砂和黏土次之，平均粒径从底部的1.6 Φ左右向上逐渐波动变化至4.5 Φ左右，正粒序明显，平均值3.2 Φ；分选系数大多变化于0.5～2.2（平均1.6），分选较好—较差；偏态主要变化于−0.50～0.75（平均0.42），以正偏为主；峰度主要变化于0.5～3.0（平均1.42），下段中下部沉积物的峰度以尖锐为主，上部沉积物以较扁平—扁平为主。
上段为深灰色粉砂-黏土质粉砂、粉砂与黏土质粉砂密集互层（毫米级）夹少量细砂条带（＜1 cm）（图3d），生物扰动较弱；在中部（32.54～34.2 m）夹少量碳屑线理（3～4 mm）。上段沉积物不含底栖有孔虫（图4）。
上段沉积物的粒度分析结果（图5）显示，沉积物以粉砂为主，砂和黏土次之，平均粒径主要变化于4～6.5 Φ，平均值5.43 Φ；分选系数主要变化于1.5～2.5，有向上变大（即分选性变差）的趋势（平均1.86），分选较差—差；偏态主要变化于0.1～0.5（平均0.31），有向上变小的趋势，为正偏—极正偏；峰度主要变化于0.9～1.5（平均1.11），没有明显的趋势，为较扁平—尖锐。
综上所述，DU 3被解释为河流沉积，其中下段为河床-点坝沉积，上段为洪泛平原沉积。
3.1.3 DU 2（29.25～25.93 m）
DU 2与下伏DU 3之间为突变接触，该单元以棕色（混杂灰绿色）黏土质粉砂和粉砂为主（图3e），比较致密，偶见钙质结核，生物扰动较强。下部（28.19～28.60 m）夹深灰色粉砂质细砂薄层，分选一般，见板状/槽状交错层理，偶见贝壳碎片，其底部为侵蚀面与下伏细粒沉积物分隔。本单元不含底栖有孔虫（图4）。
粒度分析结果（图5）显示，DU 2沉积物以粉砂为主，黏土次之，砂又次之（仅在下部的薄夹层（28.19～28.60 m）为主）；除去下部的粉砂质细砂薄层，平均粒径为5.0～7.2 Φ，平均值5.49 Φ，细砂层的平均粒径则为3.4～4.3 Φ（平均值3.65 Φ）；分选系数为1.1～2.5（平均值1.72），分选较差—差；偏态主要变化于0.10～0.37（平均值0.18），以正偏为主；峰度主要变化于0.9～1.4（平均值0.81），为较扁平—尖锐。
从岩性来看，DU 2沉积物可与长江三角洲地区“杂色硬黏土层”[29]对比，可解释为末次盛冰期下切河谷间地上的洪泛平原沉积，且受到了成岩改造；其中的细砂薄层被解释为洪泛平原上的决口扇沉积[30]。
3.1.4 DU 1（25.93～0 m）
DU 1与下伏DU 2之间为侵蚀接触，根据岩性本单元可划分为3段：下段（25.93～21.39 m）、中段（21.39～ 3.08 m）和上段（3.08～0 m）。
下段（25.93～21.39 m）：由两个显示正粒序的旋回构成，两个旋回以23.37 m为界（图2）。下旋回（25.93～23.37 m）与下伏地层侵蚀接触，包括：最底部6 cm厚的深黄灰色中砂，含细砂和粗砂，分选差，含较多贝壳碎片，向上变为深灰色细砂（含中砂）夹较多深灰色黏土质条带（厚5～15 mm）（图3f）、深灰-深黄灰色细砂薄层（厚0.15～0.26 m）与深灰黏土质条带夹细砂-粉砂线理-透镜体的薄层（厚0.15～0.50 m）相交替，潮汐层理发育，细砂层分选一般且显示板状/槽状交错层理，贝壳碎片零星分布，生物扰动较弱。上旋回（23.37～21.39 m）与下伏旋回之间为侵蚀接触，其底部为0.17 m厚的深灰色中砂（含细砂），分选一般，含大量贝壳碎片及较多泥砾；向上变为深灰色细砂夹黏土质条带（3～6 mm）（显示正粒序），细砂分选一般，潮汐层理发育，偶见贝壳碎片。
下段的底部缺失底栖有孔虫，中上部含少量底栖有孔虫，其丰度和简单分异度都较低，其平均值分别为942枚/50g和8.4种。底栖有孔虫的优势种是滨岸种Ammonia beccarii (Linné) vars.、滨岸-内陆架种Elphidium magellanicum Heron-Allen et Earland和冷水种Protelphidium tuberculatum (d'Orbigny)，次要种主要包括近岸陆架种Epistominella naraensis (Kuwano)、内陆架常见种Elphidium advenum (Cushman)、近岸低氧种Ammonia pauciloculata (Phleger et Parker)和沿岸低盐水种Cribrononion vitreum Wang。
中段（21.39～ 3.08 m）：与下段突变接触，9.20 m以下以深灰色黏土质粉砂与细砂-粉砂密集互层（毫米至厘米级）为主（图3g），潮汐层理发育，生物扰动中等，偶见贝壳碎片；9.20～3.08 m粒度有所变粗，以深灰色黏土质粉砂与细砂密集互层（毫米至厘米级）为主，其中7.74～5.20 m为深灰色细砂夹少量（下部）—较多（上部）黏土质条带（毫米级）（与下伏沉积物为侵蚀接触）。生物扰动弱—中等。
中段底栖有孔虫的丰度和简单分异度都比下段有明显增加，其平均值分别为20606枚/50g和25.8种；优势种是近岸陆架种E. naraensis和近岸低氧种A. pauciloculata，次要种主要包括近岸低氧种Ammonia convexidorsa Zheng 和Bolivina cochei Cushman et Adams、滨岸种A. beccarii vars.，滨岸-内陆架种E. magellanicum、冷水种P. tuberculatum、内陆架常见种E. advenum和沿岸低盐水种C. vitreum。
上段（3.08～0 m）：从底部的深灰黄色黏土质粉砂夹较多粉砂-细砂条带（厘米级）、粉砂，向上渐变为深灰-深黄灰色的黏土质粉砂，中上部含较多碳质斑点和少量碳屑，潮汐层理发育，生物扰动中等—强烈。
上段有孔虫的丰度和简单分异度的平均值分别为34934枚/50g和22.8种；优势种是滨岸种A. beccarii vars.，次要种主要包括近岸低氧种A. pauciloculata和B. cochei以及A. convexidorsa、近岸陆架种E. naraensis。
综上，DU 1下段被解释为近岸潮道沉积[31]；中段被解释为浅水潮下带-内陆架沉积，其中上部（7.74～5.20 m）的细砂层为潮道沉积；上段为潮坪沉积。
3.2 测年结果
3.2.1 OSL测年数据可靠性分析
以分别位于JC-1205孔31.37 m和41.67 m的两个代表性样品为例，其OSL信号生长曲线（Growth curve）和衰减曲线（Decay curve）如图6所示，两个样品的OSL信号均衰减较快，大约在2 s内衰减到本底值，表明其OSL信号以快组分为主。此外，沉积物样品OSL年龄的可靠性还与其晒退情况密切相关。前人研究表明，搬运距离短、沉积过程快以及水体浑浊度高等因素限制了河流和三角洲等水成沉积物的曝光时间和几率[32-33]，可能导致样品中只有部分颗粒的释光信号被归零或降低至可忽略水平[34]。该现象在海岸带–陆架区的沉积物（尤其是水成沉积物）光释光测年中尤为显著：样品中若混入了未完全晒退的矿物颗粒，获得的De值将会由于残余剂量的影响而导致高估。图7显示了上述两个样品De值分布的放射图（Radial Plot），较高的离散度（OD）表明样品存在不完全晒退现象，通过中值年代模型（Central age model, CAM）和最小年代模型（Minumin age mode, MAM）分别计算得到的De值是有所差别的。本文采用MAM模式得到的De值，最有可能代表其真实值[35-36]。
![]()
图 6 代表性样品1(31.37 m)和样品2(41.67 m)的光释光信号衰减曲线和生长曲线Figure 6. The decay curves and growth curves of the representative sample 1 (31.37 m) and sample 2 (41.67 m)![]()
图 7 代表性样品1(31.37 m)和样品2(41.67 m)De值分布的放射图OD：离散度，CAM：中值年代模型，MAM：最小年龄模型。Figure 7. Radial plots of De values for the representative sample 1 (31.37 m) and sample 2 (41.67 m)OD: over-dispersion, CAM: central age model, MAM: minimum age model.所测矿物的饱和剂量及其剂量率共同限定了OSL测年的年龄上限。对石英颗粒而言，其饱和剂量在200 Gy左右[37-38]，以2～3 Gy/ka的剂量率估算，其年龄上限约为60～100 ka。DU 3的OSL测年样品，其De在84.4～138.2 Gy之间，并未超出石英的饱和剂量，且样品的剂量率较低（约2 Gy/ka），因此我们认为在DU 3测年样品中所得到的大约 60～70 ka的年龄并未超出OSL年龄上限。
3.2.2 OSL测年结果
如表2和图2所示，11个OSL测年数据分布在DU 4（3个）和DU 3（6个）和DU 1（2个）。DU 4中的3个OSL年龄分布在单元的顶部、中部和下部，从上至下显示递增的年龄序列，年龄范围为143±14～160±22 ka。考虑到石英颗粒OSL测年的上限是约80 ka[37], 这3个年龄数据供参考。
DU 3中的6个OSL年龄数据分布在单元的中部至上部，除了中上部的2个年龄数据（53.2±7.4 ka和48.3±7.3 ka）较低之外，其余4个年龄数据介于61.5±5.7～67.7±4.8 ka。
DU 1中的2个OSL测年数据分布在下段的中部和上部，分别为4.0±0.2 ka和1.2±0.1 ka。
3.2.3 AMS 14C测年结果
在全部14个AMS 14C测年数据中，10个年龄数据位于DU 3的底部（孔深60 m以下）（表2和图2），它们主要是利用淡水腹足类（Angulyagra sp. ）的壳体所测定的，其中7个年龄数据超越了AMS 14C测年的上限（＞43500 yr BP），另外3个年龄数据为42.1～42.8 ka。考虑到本单元OSL的测年结果总体上＞60 ka，这3个较年轻的AMS 14C年龄数据被解释为测年样品在沉积或埋藏过程中受到了轻碳的污染而致使年龄变得显著年轻[39]。
另外4个AMS 14C测年数据位于DU 1，其中2个位于上段的顶部、1个位于中段的中部、1个位于下段的上部；最下部的年龄约为3.7 ka，上面几个年龄基本上是向上递减的。
4. 讨论
由上述的沉积相分析可知，JC-1205孔地层序列中除了最上部的DU 1单元之外，其余沉积单元（DU 2—DU 4）都是河流沉积，它们的成因显然与海平面的变化密切相关。晚第四纪时期，入海河流的侵蚀基准面－海平面经历了多次大幅度（＞100 m）的升降变化[40]，因此在海岸带地区发育的大型河流体系也伴随着海平面波动而经历了多次下切的过程，发育了多期深切河谷及其充填沉积序列[30, 41]。依据JC-1205孔地层序列的测年结果，结合晚第四纪全球海平面变化特征（图8），对该孔所揭示的沉积序列开展地层划分，探讨各沉积单元的形成时代。
![]()
MIS 6（191～130 ka）[43] 时期是中更新世以来大陆冰盖范围仅次于末次冰期最盛期（LGM）（26.0～19.5 ka）[44]的一个大冰期，期间最低海平面出现在MIS 6晚期（约140～135 ka），达到约125 m bpsl（现今海平面之下）[45]。DU 4沉积地层的OSL年龄数据位于MIS 6时期，但其超过了石英颗粒OSL测年的上限，可视为该单元沉积物形成时代的下限。由于MIS 6时期中国东部海岸带和大部分陆架区出露地表，且海平面的大幅度下降势必引起该区域河流体系的强烈下切，加之MIS 6之后为MIS 5e（约130～116 ka）高海平面时期（高于现今海平面6～9 m）[46]，比较有利于MIS 6河流沉积物的保存。因此，河流沉积DU 4被解释为形成于MIS 6。由于孔深的限制，JC-1205孔只揭露了MIS 6时期的洪泛平原沉积，而没有揭露同期的河床沉积。根据以往的研究，位于弶港南侧（JC-1205孔西北方向）的钻孔揭示了晚第四纪早期的下切谷[47]，推测DU 4（洪泛平原沉积）与该下切谷的河床沉积是成因上相关联的同期河流沉积。
全球海平面在MIS 5e之后呈波动下降，至MIS 4（约71.0～60.0 ka）海平面高程显著降低至约80～100 m bpsl[45, 48]；进入MIS 3（约60.0～30.0 ka）[49-50]的早期后海平面明显上升，在中国东部可以达到约20～30 m bpsl[26, 51]，然后海平面波动下降至MIS 3末期的80～90 m bpsl[46, 52]。除去DU 3中部2个明显不可靠的年轻OSL年龄数据，DU 3上部和中部其他OSL年龄数据主要为60～70 ka。考虑到DU 3顶界面和底界面的高程分别为约28 m bpsl和约67 m bpsl，加之晚第四纪以来研究区沉降缓慢，晚第四纪中期下切河谷及其河流沉积DU 3可以合理地解释为形成于MIS 4。
MIS 2（约30.0～14.0 ka）全球海平面再进一步下降，在LGM（26.0～19.5 ka）[44]海平面可达～125 m bpsl[19]；MIS 2期间研究区发育下切河谷。如前所述，DU 2洪泛平原沉积可与长江三角洲地区广泛发育的“杂色硬黏土层”[29]对比。因此河流沉积DU 2被解释为形成于MIS 2，此时钻孔位于MIS 2下切河谷北侧的河间地，下切河谷位于钻孔以南几千米之远（图1）。
自MIS 1（14 ka以来）的早期开始，南黄海西部陆架逐渐被冰后期海侵所覆盖。DU 1的中段和上段为浅水潮下带-内陆架沉积和潮坪沉积，其下段为近岸潮道沉积，其年龄数据表明DU 1形成于约4.0 ka以来近岸海洋环境。我们推测，DU 1在约4.0 ka开始发育的潮道沉积侵蚀了下伏MIS 1早—中期的海洋沉积和部分MIS 2晚期的河流沉积，从而造成了DU 1和DU 2之间超过10 kyr左右的沉积间断；该潮道向下侵入到下伏LGM河流沉积，加之潮道沉积期间较强的水动力条件，导致该潮道沉积（DU 1下段）的底部缺少底栖有孔虫、中上部底栖有孔虫的丰度和简单分异度也比其上的DU 1中段（浅水潮下带-内陆架沉积）和DU 1上段（潮坪沉积）显著降低。
从以上分析可知，在DU 4和DU 3之间缺失了MIS 5时期的沉积物，而在DU 3和DU 2之间缺失了MIS 3时期的沉积物（图8）。DU 3的沉积相特征表明，MIS 4时期JC-1205孔正好位于下切谷之中，这也得到了以往研究的印证（图1）；DU 3和DU 4之间MIS 5海洋沉积的缺失应是该下切谷发育过程中对下伏沉积物的侵蚀所致。DU 2发育于MIS 2河间地的泛滥平原环境，DU 2底界面的高程为约28 m bpsl，如果考虑到晚第四纪研究区一直经历缓慢沉降，推测DU 3形成之后钻孔位置附近因地势较高而缺失了MIS 3的近岸海洋沉积，即钻孔位置附近在MIS 3期间出现了沉积间断。
晚第四纪长江三角洲地区沉降速率较小，下切谷的发育主要受控于海平面变化。在MIS 6、MIS 4和MIS 2低海平面时期，钻孔位置附近区发育河流下切河谷或河间地，后期的下切谷（如MIS 4时期）对前期的下切河谷沉积充填形成了破坏，使得它们之间相对高海平面时期（如MIS 5）的近岸海洋沉积物被侵蚀殆尽，或者因地势较高而产生了沉积间断（如MIS 3），导致这3期河流沉积物直接叠置（图8）。这与以前的研究报道[11-12]是一致的。如果考虑到MIS 1晚期近岸海洋沉积（DU 1）与其下伏河流沉积（DU 2）之间也存在超过10 kyr的沉积间断，研究区晚第四纪沉积序列的4个沉积单元之间都存在重要的沉积缺失。
5. 结论
（1）江苏南通近岸区JC-1205孔揭示了该区晚第四纪沉积序列可划分为4个沉积单元（从下至上DU 4—DU 1）：中—下部3个沉积单元（DU 4—DU 2）为3期直接接触的河流沉积，上部DU 1为近岸海洋沉积（从下至上包括潮道沉积、浅水潮下带-内陆架沉积和潮坪沉积）。
（2）晚第四纪海平面变化是控制钻孔位置附近区沉积序列形成的主要因素。晚第四纪早期和中期的下切河谷及其充填沉积（DU 4和DU 3）分别形成于MIS 6和MIS 4，而与晚期下切谷相关的河间地的洪泛平原沉积（DU 2）形成于MIS 2；MIS 4时期下切谷的发育使得其下伏的MIS 5海洋沉积被侵蚀殆尽，而MIS 3时期相对较高的地势导致了沉积间断的出现。
（3）DU 1形成于MIS 1晚期，其底部的潮道沉积在其发育过程中对下伏沉积物的侵蚀作用形成了DU 1和DU 2之间超过10 kyr的沉积缺失。
致谢：本文作者多年来在开展海岸带地质调查和海洋沉积学研究过程中得到了青岛海洋地质研究所何起祥先生的悉心指导和无私帮助，谨以此文纪念何起祥先生！
-
![]()
![]()
![]()
图 3 各级层序界面地震反射特征及识别标志
Figure 3. Seismic reflection characteristics and identification of various sequence interfaces
![]()
图 4 文昌组控洼断层活动性分析
Figure 4. Analysis on the activity of depression-controlling faults of Wenchang Formation
![]()
图 5 文昌组三级层序地震剖面
Figure 5. Seismic profile showing the third-order sequence of Wenchang Formation
![]()
图 6 不同沉积相类型地震相、测录井、岩心、地化特征
Figure 6. Seismic facies, logging, core and geochemical characteristics of different sedimentary facies
![]()
图 7 西江36洼与番禺4洼文昌组三级层序地震相分布
Figure 7. Distribution of seismic facies of each three-order sequence of Wenchang Formation in Xijiang 36 sag and Panyu 4 sag
![]()
图 8 西江36洼与番禺4洼文昌组三级层序沉积体系展布
Figure 8. Distribution of three-level sequence sedimentary systems of the Wenchang Formation in Xijiang 36 sag and Panyu 4 sag
表 1 西江36洼与番禺4洼烃源岩条件对比
Table 1 Comparison of source rock conditions between Xijiang 36 sag and Panyu 4 sag
西江36洼 番禺4洼 层序 中深湖规模 物源供给量 生烃
潜力值中深湖规模 物源供给量 生烃
潜力值面积
/km2平均
厚度/m面积
/km2平均
厚度/m面积
/km2平均
厚度/m面积
/km2平均
厚度/m上文昌 文一段 0 0 − − 0 0 0 − − 0 文二段 0 0 − − 0 180 400 150 250 1.92 文三段 42 450 60 125 2.52 300 600 147 200 6.12 下文昌 文四段 83 650 80 110 6.13 132 600 110 100 7.2 文五段 45 500 70 120 2.68 50 350 80 120 1.82 文六段 0 0 − − 0 0 0 − − 0 注：生烃潜力值=（中深湖面积×中深湖平均厚度）/（物源面积×物源平均厚度）。 
下载: 导出CSV
-
[1] 朱筱敏, 黄捍东, 代一丁, 等. 珠江口盆地番禺4洼文昌组层序格架与沉积体系研究[J]. 岩性油气藏, 2014, 26(4):1-8 doi: 10.3969/j.issn.1673-8926.2014.04.001 ZHU Xiaomin, HUANG Handong, DAI Yiding, et al. Study on depositional system and sequence framework of Wenchang Formation in Panyu 4 sag of the Pearl River Mouth Basin [J]. Lithologic Reservoirs, 2014, 26(4): 1-8. doi: 10.3969/j.issn.1673-8926.2014.04.001
[2] 张向涛, 朱俊章, 熊万林, 等. 番禺4洼文昌组烃源岩生物标志化合物特征与油源判识[J]. 中国海上油气, 2020, 32(4):12-23 ZHANG Xiangtao, ZHU Junzhang, XIONG Wanlin, et al. Biomarker characteristics and oil-source discrimination of source rocks in Wenchang Formation of Panyu 4 sag [J]. China Offshore Oil and Gas, 2020, 32(4): 12-23.
[3] 代一丁. 珠江口盆地西江南洼古近系构造演化与沉积特征[J]. 中国海上油气, 2013, 25(3):1-7 DAI Yiding. Paleogene tectonic evolution and sedimentation in South Xijiang Sag, Pearl River Mouth Basin [J]. China Offshore Oil and Gas, 2013, 25(3): 1-7.
[4] 江宁, 全志臻, 张向涛, 等. 珠江口盆地番禺4洼古近系层序地层及储层分布预测[J]. 天然气勘探与开发, 2015, 38(4):23-27 doi: 10.3969/j.issn.1673-3177.2015.04.006 JIANG Ning, QUAN Zhizhen, ZHANG Xiangtao, et al. Paleogene sequence stratigraphy and reservoir distribution, Panyu 4 Sag, Pearl River Mouth Basin [J]. Natural Gas Exploration and Development, 2015, 38(4): 23-27. doi: 10.3969/j.issn.1673-3177.2015.04.006
[5] 吴宇翔, 舒誉, 丁琳, 等. 珠江口盆地番禺4洼文昌组基于层序地层格架约束下的优质烃源岩预测[J]. 海洋地质前沿, 2021, 37(3):41-49 WU Yuxiang, SHU Yu, DING Lin, et al. Prediction of high quality source rocks based on sequence stratigraphic framework of Wenchang Formation, Panyu 4 sag, the Pearl River Mouth Basin [J]. Marine Geology Frontiers, 2021, 37(3): 41-49.
[6] 张丽丽, 舒誉, 蔡国富, 等. 珠江口盆地东部始新世渐新世沉积环境演变及对烃源条件的影响[J]. 石油学报, 2019, 40(S1):153-165 doi: 10.7623/syxb2019S1013 ZHANG Lili, SHU Yu, CAI Guofu, et al. Eocene-Oligocene sedimentary environment evolution and its impact on hydrocarbon source conditions in eastern Pearl River Mouth basin [J]. Acta Petrolei Sinica, 2019, 40(S1): 153-165. doi: 10.7623/syxb2019S1013
[7] 施和生, 何敏, 张丽丽, 等. 珠江口盆地(东部)油气地质特征、成藏规律及下一步勘探策略[J]. 中国海上油气, 2014, 26(3):11-22 SHI Hesheng, HE Min, ZHANG Lili, et al. Hydrocarbon Geology, accumulation pattern and the next exploration Stratagy in there eastern Pearl River Mouth basin [J]. China Offshore Oil and Gas, 2014, 26(3): 11-22.
[8] 舒誉, 施和生, 杜家元, 等. 珠一坳陷古近系油气成藏特征及勘探方向[J]. 中国海上油气, 2014, 26(3):37-42 SHU Yu, SHI Hesheng, DU Jiayuan, et al. Paleogene characteristics in hydrocarbon accumulation and exploration direction in Zhu I depression [J]. China Offshore Oil and Gas, 2014, 26(3): 37-42.
[9] 施和生, 舒誉, 杜家元, 等. 珠江口盆地古近系石油地质[M]. 北京: 地质出版社, 2017: 1-4 SHI Hesheng, SHU Yu, DU Jiayuan, et al. Petroleum Geology of the Paleogene in the Pearl River Mouth Basin[M]. Beijing: China Geology Press, 2017: 1-4.
[10] 陈长民. 珠江口盆地东部石油地质及油气藏形成条件初探[J]. 中国海上油气(地质), 2000, 14(2):73-83 CHEN Changmin. Petroleum geology and conditions for hydrocarbon accumulation in the eastern Pearl River mouth basin [J]. China Offshore Oil and Gas (Geology), 2000, 14(2): 73-83.
[11] 米立军, 张向涛, 陈维涛, 等. 珠江口盆地珠一坳陷古近系油气富集规律及下一步勘探策略[J]. 中国海上油气, 2018, 30(6):1-13 MI Lijun, ZHANG Xiangtao, CHEN Weitao, et al. Hydrocarbon enrichment law of Paleogene Zhu 1 depression and its next exploration strategy in Pearl River Mouth basin [J]. China Offshore Oil and Gas, 2018, 30(6): 1-13.
[12] 施和生, 杜家元, 梅廉夫, 等. 珠江口盆地惠州运动及其意义[J]. 石油勘探与开发, 2020, 47(3):447-461 SHI Hesheng, DU Jiayuan, MEI Lianfu, et al. Huizhou Movement and its significance in Pearl River Mouth Basin, China [J]. Petroleum Exploration and Development, 2020, 47(3): 447-461.
[13] Vail P R, Mitchum P M. Seismic stratigraphic and global changes in sea leavel. Parts Ⅰ-Ⅱ [J]. AAPG Memoir, 1977, 26: 51-212.
[14] Mitchum P M, Sangree J B, Vail P R, et al. Recognizing sequences and systems tracts from well logs, seismic data, and biostratigraphy: Examples from the late Cenozoic of the Gulf of Mexico[M]//Weimer P, Posamentier H. Siliciclastic Sequence Stratigraphy: Recent Developments and Applications. California: American Association of Petroleum Geologists, 1993: 163-199.
[15] Mitchum R M, Van Wagoner J C. High-frequency sequences and their stacking patterns: sequence-stratigraphic evidence of high-frequency eustatic cycles [J]. Sedimentary Geology, 1991, 70(2-4): 131-160. doi: 10.1016/0037-0738(91)90139-5
[16] 王英民. 对层序地层学工业化应用中层序分级混乱问题的探讨[J]. 岩性油气藏, 2007, 19(1):9-15 doi: 10.3969/j.issn.1673-8926.2007.01.002 WANG Yingmin. Analysis of the mess in sequence hierarchy applied in the industrialized application of the sequence stratigraphy [J]. Lithologic Reservoirs, 2007, 19(1): 9-15. doi: 10.3969/j.issn.1673-8926.2007.01.002
[17] 卓海腾, 王英民, 徐强, 等. 南海北部陆坡分类及成因分析[J]. 地质学报, 2014, 88(3):327-336 ZHUO Haiteng, WANG Yingmin, XU Qiang, et al. Classification and genesis of continental slopes on the northern South CHINA Sea [J]. Actav Geologica Sinica, 2014, 88(3): 327-336.
[18] 于兴河, 姜辉, 李胜利, 等. 中国东部中、新生代陆相断陷盆地沉积充填模式及其控制因素: 以济阳坳陷东营凹陷为例[J]. 岩性油气藏, 2007, 19(1):39-45 doi: 10.3969/j.issn.1673-8926.2007.01.007 YU Xinghe, JIANG Hui, LI Shengli, et al. Depositional filling models and controlling factors on Mesozoic and Cenozoic fault basins of terrestrial facies in eastern China: A case study of Dongying Sag of Jiyang Depression [J]. Lithologic Reservoirs, 2007, 19(1): 39-45. doi: 10.3969/j.issn.1673-8926.2007.01.007
[19] 冯有良. 大民屯凹陷沙四段—沙三段层序地层格架及岩性油气藏预测[J]. 岩性油气藏, 2008, 20(4):14-19 doi: 10.3969/j.issn.1673-8926.2008.04.003 FENG Youliang. Sequence stratigraphy framework and lithologic reservoirs prediction of the third and fourth members of Shahejie Formation, Damintun Depression [J]. Lithologic Reservoirs, 2008, 20(4): 14-19. doi: 10.3969/j.issn.1673-8926.2008.04.003
[20] 吴嘉鹏, 王英民, 马贵明, 等. 孟加拉扇某区块地层层序划分及典型地震相识别[J]. 岩性油气藏, 2010, 22(2):69-73 doi: 10.3969/j.issn.1673-8926.2010.02.012 WU Jiapeng, WANG Yingmin, MA Guiming, et al. Sequence stratigraphy and typical seismic facies in one exploration block of Bengal Fan [J]. Lithologic Reservoirs, 2010, 22(2): 69-73. doi: 10.3969/j.issn.1673-8926.2010.02.012
[21] 张功成, 梁建设, 徐建永, 等. 中国近海潜在富烃凹陷评价方法与烃源岩识别[J]. 中国海上油气, 2013, 25(1):13-19 ZHANG Gongcheng, LIANG Jianshe, XU Jianyong, et al. An evaluation method of potential hydrocarbon-rich sags and their source rock identification offshore China [J]. China Offshore Oil and Gas, 2013, 25(1): 13-19.
[22] 郭来源. 陆相断陷湖盆富有机质页岩非均质性及其控制因素分析: 以泌阳和沾化凹陷为例[D]. 中国地质大学博士学位论文, 2017 GUO Laiyuan. Heterogeneity and controlling factors of organic-rich shale in continental rift basin: a case study of Biyang and Zhanhua depressions[D]. Doctor Dissertation of China University of Geosciences, 2017.
-
期刊类型引用(1)
1. 邹庆超，石学法，葛淑兰，吴永华，邹建军，贺湘锋，Sergey A.Gorbarenko，刘建兴. 氧同位素3早期以来千年尺度气候事件在日本海北部的磁学记录. 海洋学报. 2024(04): 106-121 . 
百度学术
其他类型引用(1)
计量
- 文章访问数: 2007
- HTML全文浏览量: 374
- PDF下载量: 45
- 被引次数: 2
