DOI： 10.11817/j.issn.1672-7207.2021.05.013

渤海湾盆地沾化凹陷罗69井沙三下亚段旋回地层学研究

张若琳^{1, 2}，金思丁^{1, 2}

(1. 成都理工大学 能源学院，四川 成都，610059；

2. 油气藏地质及开发工程国家重点实验室，四川 成都，610059)

摘要：对渤海湾盆地沾化凹陷罗69井的自然伽马测井数据进行频谱分析和滑动窗口频谱分析。通过对0.405 Ma和0.100 Ma的周期进行高斯带通滤波，并以沙三下亚段顶和底的磁性地层年龄为锚点，建立罗69井沙三下亚段的天文年代标尺。沙三下亚段总共记录了6个长偏心率周期和24个短偏心率周期，将0.405 Ma的滤波曲线作为四级层序沉积旋回的参考曲线，在沙三下亚段地层划分出6个四级层序。研究结果表明：罗69井沙三下亚段记录了稳定的天文旋回；由长偏心率、短偏心率、斜率和岁差周期产生的地层旋回厚度分别为39.4，9.74，3.94~2.54和1.92~1.70 m。罗69井沙三下亚段的最优沉积速率为8.2×10³ cm/Ma，沉积时限为2.38 Ma。在沙三下亚段早期沉积地层中，长偏心率周期控制了氧化还原条件从而促进了有机质的富集；而在沙三下亚段晚期沉积地层中，长偏心率周期控制了初级生产力，是有机质富集的主要影响因素。

关键词：旋回地层学；罗69井；湖相细粒沉积；沾化凹陷；沙三下亚段

中图分类号：P539.2          文献标志码：A               开放科学(资源服务)标识码(OSID)

文章编号：1672-7207（2021）05-1516-16

Cyclostratigraphy research on lower Member 3 of Shahejie Formation in Well Luo 69 in Zhanhua Sag Bohai Bay Basin

ZHANG　Ruolin^{1, 2}, JIN　Siding^{1, 2}

(1. School of Energy, Chengdu University of Technology, Chengdu, 610059, China;

2. State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Chengdu University of Technology, Chengdu 610059, China.)

Abstract: The nature gamma logging data of Well Luo 69 were chosen for a proxy, and the detailed spectra analysis and evolutional spectra analysis were conducted in the Zhanhua Sag, Bohai Bay Basin. The astronomical time scale was constructed for the Well Luo 69 in the Es³l based on the Gaussian filtering on 0.405 Ma and 0.100 Ma cycles and the anchored age was based on the magnetostratigraphic data of the lower part of Member 3 of the Shahejie Formation. The Es³l was found to have 6 long and 24 short eccentricity cycles; the 0.405 Ma long eccentricity-filtering curve could be regarded as reference curve for fourth order sequences division. It could be divided into 6 fourth-order sequences in the Es³l. The results show that the Well Luo 69 in the Es³l preserves stable Milankovitch cycles. The major cycles of 39.4, 9.74, 3.94-2.54 and 1.92-1.70 m are controlled by long eccentricity, short eccentricity, obliquity and precession, respectively. The optimum sedimentation accumulation rate(SAR) is 8.2×10³ cm/Ma, and the accumulation time of the Es³l of the Well Luo 69 can be estimated at 2.38 Ma. The preservation condition controlled by long eccentricity period in the lower Es³l, and the productivity controlled by long eccentricity period in the upper Es³l can be demonstrated, which are the main factors of organic accumulation.

Key words: Cyclostratigraphy; Well Luo 69; lacustrine fine-gained sediments; Zhanhua Sag; Es³l Formation

在地学研究领域拓宽的过程中，很多地质现象和过程无法仅从地球系统本身得到有效解释的答案，这促使了旋回地层学即米兰科维奇理论^[1]的诞生：由地球轨道驱动力造成沉积地层的旋回性^[2-7]。细粒沉积岩同时具有相对不间断、分辨率高、古环境信息丰富且持续时间长的特征，是保存天文信号的理想载体，也是开展旋回地层研究的合适对象。JIN等^[8-12]在细粒沉积岩中发现了米兰科维奇旋回。曾有学者认为湖泊沉积环境相对更易受构造作用、陆源输入的影响，不利于天文信号的保存，国际上细粒沉积岩的旋回地层研究多集中于海相沉积中^[13-15]。自20世纪90年代以来，在多处湖相沉积中发现了典型的米兰科维奇旋回^[16-18]。随非常规陆相细粒沉积岩油气勘探的开展，WU等^[19]在松辽盆地晚白垩地层、鄂尔多斯盆地晚三叠延长组地层^[20]、渤海湾盆地古近系地层^[21-23]等湖相细粒沉积物中发现了米氏旋回。由于天文轨道周期从根本上控制着气候变化，气候变化不仅影响着湖泊的机械沉积作用，而且影响着营养物质供给、古生产力条件、氧化还原条件等^[24-25]，利用旋回地层学理论不仅可以尝试以轨道周期控制沉积作用的角度进行解释，从成因上解决地层划分问题，亦能从天文角度还原湖泊沉积的古气候、古环境变化特征，有机质富集因素，对湖相细粒沉积岩的勘探开发具有重要科学意义。本次研究拟对渤海湾盆地济阳坳陷沾化凹陷罗69井的自然伽马测井序列进行旋回地层学研究，运用普通频谱(MTM)、滑动窗口频谱(LAH)和相关系数法(“COCO”)分析技术，匹配最优沉积速率，排除误差信号，提取沙三下亚段(Es³l)细粒沉积岩中的米兰科维奇旋回信息，通过滤波与调谐，建立天文年代标尺，计算不同时期的沉积速率，并将长偏心率周期(0.405 Ma)的滤波曲线作为沙三下亚段内部四级层序划分的参考曲线，定量划分细粒沉积岩高频层序地层。最后以GR曲线为主，并结合不同的能反映沉积环境、古气候的地球化学替代指标，验证天文驱动对有机质富集控制作用的模型。

1  区域地质背景

济阳坳陷位于渤海湾盆地东南部，是渤海湾盆地的1个次级构造单元，也是渤海湾盆地油气最丰富的地区之一^[26]。济阳坳陷总面积约为2.65×10⁴ km²，内部发育东营、沾化、惠民、车镇4个凹陷及若干凸起，展现出“多凸多凹、凹凸相间”的构造格局(图1(a))^[27-28]。始新世沙三段是1段厚度为120~700 m的深湖相沉积，湖盆中心厚度可达1 000 m以上，主要岩性为富有机质的细粒沉积岩，含砂岩、灰岩夹层，是古近纪的主要生油岩层段^[29]。沙三段从下到上可划分为3段，分别为沙三下亚段(Es³l)、沙三中亚段(Es³m)、沙三上亚段(Es³u)(图1(b))。

图1　渤海湾盆地济阳坳陷构造简图及地层垂向充填序列^{[23-24, 30]}

Fig. 1　Structural sketch and stratigraphic filling histogram of Jiyang Depression, Bohai Bay Basin^{[23-24, 30]}

沾化凹陷位于济阳坳陷东北部，属于济阳坳陷的1个次级构造单元，其南部为陈家庄凸起，北部为义和庄凸起，是典型的“北断南超”的箕状盆地，面积约为2 800 km^2[30]。从始新世早期开始，沾化凹陷进入断拗阶段，湖盆持续下沉，接受沙河街组沉积，其中沙三下亚段(Es³l)沉积时期处于沾化凹陷的深断陷期，盆地进入非补偿阶段，气候潮湿，处于稳定而持久的深湖，半深湖的沉积环境，发育了厚为100~500 m的暗色泥页岩夹灰岩沉积，为沾化凹陷的主力烃源层。目标井罗69井位于沾化凹陷的罗家鼻状构造区，沙三下亚段(2 930~3 127.5 m)属1套深灰色—灰褐色的深湖—半深湖相富有机质纹层状、层状泥页岩、灰质泥页岩沉积，累计厚度为197.5 m^[31]。这口井在2 911.00~3 140.75 m井段实施了连续密闭取芯，彭丽等^[28]将罗69井沙三下亚段进一步划分为湖扩体系域(EST)、早期高位体系域(EHST)和晚期高位体系域(LHST)。

2  测试方法与数据预处理

本次研究采用的软件包为基于Matlab平台下的Acycle2.2^[32]。对比Kaleidagraph^TM[33]，Analyseries 2.0.8^[34]，SSA-MTM^[35]，基于R软件的“Astrochron”程序包^[36]以及Past3^[37]等，Acycle2.2将旋回地层学研究所要做的所有工作统一编写在同一种操作系统下的同一个程序中，无须来回切换操作系统或软件，操作方便简洁。

本研究中涉及的主要计算方法包括Multi-taper method (MTM)频谱分析、Fast Fourier transform(LAH)滑动窗口频谱分析、滤波分析、相关系数“COCO”及“eCOCO”分析等。罗69井的GR测井序列数值范围为20~100 Ω·m，采样间隔为0.1 m。GR曲线的高值与暗色泥页岩对应，低值与灰岩及粉砂岩相对应，肉眼可见明显的旋回特征(图2)^[28]。由于GR测井曲线是由各种地质因素所引起的周期性变化的综合反映，其中也包括与地质因素无关的干扰信号，且地层记录会出现随时间或深度逐渐增加或减小的趋势。为了更好地提取GR测井序列中存在的天文旋回频段，在进行频谱分析之前，要将这种低频部分的形变影响去除(见图2中去趋势化后的Detrended GR曲线)。去趋势化后GR数据序列中受天文周期影响的信号将更为突出。

图2　罗69井旋回特征

Fig. 2　Cycle characteristic of Well Luo 69

3  旋回地层学研究

3.1　天文周期解决方案

地球轨道3个重要参数为偏心率(eccentricity)、斜率(obliquity)及岁差(precession)，受太阳系其他行星及月球与地球间的引力影响，这3个参数在地史时期呈周期性缓慢变化。LASKAR等^[38]引入月球潮汐力的损耗等因素将米兰科维奇理论模型延伸至过去的250 Ma，提出La2004解决方案，后采用天文历表INPOP08进行计算^[39]，得到更精确的偏心率变化曲线，提出La2010解决方案。目前的天文周期解决方案对50 Ma以来的地球轨道参数的计算非常精确。

LIU等^[22]通过对渤海湾盆地沙河街组四口科学探井的伽马测井数据进行天文旋回计算，得到了22~66 Ma的高精度天文年代标尺，他们将渤海湾盆地济阳坳陷沙三下亚段的顶和底年龄锚定在40.20~42.47 Ma，该年龄也被ZHAO等^[23]应用在渤海湾盆地东营凹陷牛页1井的天文年代调谐的年龄控制点中，并取得了较好的对应关系。本文拟采用La2010天文周期解决方案，计算40.2~42.47 Ma之间的偏心率、斜率与岁差变化(E+O-P)的理论值，采样间隔为0.001 Ma(图3(a))。随后对E+O-P曲线进行频谱分析与滑动窗口频谱分析(图3(b)和(c))，结果显示40.20~42.47 Ma期间的主要天文周期为0.405，0.125，0.095，0.051，0.040，0.023，0.022和0.019 Ma。其中，0.405 Ma为长偏心率(E)周期，0.125 Ma与0.095 Ma为短偏心率(e)周期；0.051 Ma与0.040 Ma属于斜率(O)周期；0.023，0.022和0.019 Ma属于岁差周期。利用计算出的各轨道参数的比值，与伽马测井序列通过频谱分析之后的谱峰进行对比，若存在相等或相近的旋回比率，则可以初步判断地层中记录了米兰科维奇旋回信息^[40]。

图3　目的层段天文周期解决方案频谱分析图

Fig. 3　Spectral analysis of target astronomical solution

3.2　深度域频谱分析

对趋势化后的伽马测井数据进行MTM(multitaper method of spectral analysis)频谱分析(图4(b))，图4(b)中，横坐标显示为频率，纵坐标为功率，功率越大表示该功率对应的旋回越明显。每个高置信度峰值(置信度大于95%)对应的厚度旋回分别为39.40，9.74，3.94，2.78，2.54，1.92和1.70 m，显示出的频谱结构具有能够明显区分开的4个优势频率带，分别为39.40~27.77，9.74~6.32，3.94~2.54和1.92~1.70 m(图4(b))。分别将4个频率谱峰值对应的周期比20.44:5.05:2.09:1.00与理论周期比进行比较，发现天文轨道周期20:5:2:1的比例关系与之非常接近。

图4　罗69井沙三下亚段GR曲线频谱分析图

Fig. 4　Spectral analysis of in Well Luo 69

为了识别旋回地层中主频率随深度的变化，观察整体数据序列是否具有稳定性，继而对罗69井沙三下亚段伽马测井曲线进行滑动窗口频谱分析。本次滑动窗口频谱分析中的窗口步长为0.3 m，滑动窗口长度为40 m(图4(c))。谱图中峰值所对应的频率大致与MTM功率谱中高置信度峰值相对应，表明各个深度段均保存了偏心率、斜率和岁差周期的信号。然而，由图4(c)可见在2 988 m及3 040 m处附近出现了明显的峰值信号不连续现象，可能意味着某些深度段序列中包含更多噪声，高频信号被具有更高功率的噪声所抑制^[41]，但更有可能的是陆源碎屑物质的输入导致沉积速率发生变化，从而使频谱的周期性发生变化^[42]。

据彭丽等^[28]对罗69的层序地层学划分结果，本次研究将GR数据序列进一步划分为3小段(2 930~2 988，2 988~3 040和3 040~3 127.5 m)进行频谱分析(图4(d)~(f))。此处只是初步估计沉积速率发生变化的区间段，对罗69井的沉积速率变化区间精确计算。在2 930~2 988 m深度段，高于95%置信度的4个频率谱峰对应的周期为44.42，10.42，4.06和1.75 m，其比值与405:95:40:22相近；在2 988~3 040 m段，高于95%置信度的4个频率谱峰对应的周期为40.25，12.48，3.89和1.78 m，其比值与405:125:40:19相近；在3 040~3 127.5 m段，高于95%置信度的4个频率谱峰对应的周期为45.02，10.31，4.52和2.38 m，其比值与405:95:40:19相近。同时，图4(d)~(f)的频谱结果显示了在不同深度段，由于沉积环境的改变，即使沉积速率发生微小变化，也会导致沉积物中保存的主要天文轨道信号发生变化。

3.3　最优沉积速率的估算

由于选用的GR测井数据是以深度为单位，尽管在95%的置信度之上检测到了功率的周期性，但并未对空间周期所代表的时间范围进行估计，因此，仅通过传统的“比值法”对天文旋回进行识别是比较粗略的。本次研究进一步采用LI等^[32]提出的“COCO(correlation coefficient)”和“eCOCO(evolutionary correlation coefficient)”分析来提高旋回识别的准确性，通过自回归模型(autoregressive-1，AR1)，利用蒙特卡洛模拟(Monte Carlo)对非天文轨道驱动的信号(噪声)进行零假设检验，增大天文轨道信号识别的精度。检验中的显著性水平(H₀，significance level)表示地层记录中零假设被错误拒绝的概率^[43]，例如0.1%的显著性水平表示有99.9%的频谱被检测到，被错误拒绝的零假设(非天文轨道驱动的信号)仅为0.1%。因此H₀越小，分析的结果越可靠。“COCO”分析同时还可以计算研究层段的沉积速率，通过“COCO”计算得罗69井沙三下亚段的平均沉积速率为7.2×10³ cm/Ma或8.2×10³ cm/Ma(图5(a))，两者都超过了要求的显著性水平(图5(b))，并且沉积速率为8.2×10³ cm/Ma，且7个天文轨道参数分量都参与的置信水平低于0.05%，可以拒绝没有轨道信号参与的零假设(图5(b)~(c))。根据前人对研究区对应层段的沉积速率的估算，本次研究中“eCOCO”分析选取(5~50)×10³ cm/Ma的速率区间，滑动步长为0.7 m，滑动窗口为35 m，蒙特卡洛模拟2 000次，得到如图6所示结果。由图6可知：数据全段与图5有良好的对应关系，且(7.2~8.2)×10³ cm/Ma的沉积速率在沙三下亚段整段保持稳定(图6(a))。因此，“eCOCO”分析可用于对罗69井进行的旋回地层学分析，其在沙三下亚段保存有完整的米兰科维奇旋回。

3.4　滤波和调谐

基于频谱分析的结果，将代表长偏心率(E：0.0405 Ma)的沉积旋回(39.40 m)通过高斯带通滤波提取出来，其滤波频率为(0.027±0.003 3) 周/m，滤波结果显示罗69井沙三下亚段约记录了6个长偏心率旋回(图2)；提取代表短偏心率e(0.100 Ma)的沉积旋回(9.74 m)通过高斯带通滤波，其滤波频率为(0.114±0.035) 周/m，滤波结果显示沙三下亚段约记录了24个短偏心率旋回(图2)；提取代表斜率O(0.040 Ma)的沉积旋回(3.94 m)通过高斯带通滤波，其滤波频率为(0.300±0.020) 周/m，滤波结果显示沙三下亚段约记录了59个斜率旋回(图2)。

由于50 Ma以来的天文模型已经十分完善且可靠，可尝试将滤波曲线与天文目标曲线对应。长偏心率0.405 Ma在地质历史时期一直显示比较稳定的状态，本次研究将选择代表长偏心率的滤波曲线建立深度-时间模型，并选取长偏心率、短偏心率、斜率3个天文旋回参数作为理论调谐曲线(图7)。由图7可见：以沙三下亚段底界面年龄42.47 Ma作为天文年代调谐的初始锚点，沙三下亚段顶界面年龄40.20 Ma为控制点，用来界定本次研究的时间范围并作为深度与时间转换的依据。当获得时间域数据列后，结合对应的深度域序列可获得罗69井沙三下亚段连续变化的沉积速率，在2 963.3 m(40.552 4 Ma)，2 993.6 m(40.957 4 Ma)，3 031.6 m(41.362 4 Ma)，3 066.7 m(41.767 4 Ma)这4处出现了沉积速率的明显变化(图7)，沉积速率从下到上存在明显的减小趋势。

图5　罗69井沙三下亚段GR曲线的“COCO”分析

Fig. 5　“COCO” analysis of detrended GR curve in Es³l of Well Luo 69

图6　罗69井沙三下亚段GR曲线的“eCOCO”分析

Fig. 6　“eCOCO” analysis of detrended GR curve in Es³l of Well Luo 69

图7　罗69井沙三下亚段天文调谐结果

Fig. 7　Astronomical tuning of GR series of Es³l in Well Luo 69

3.5　调谐后时间域频谱分析

为了检验滤波和调谐结果的准确性，对调谐后得到的时间域GR数据序列再次进行频谱分析及滑动窗口频谱分析(图8)，GR曲线频谱分析图中高于95%置信度的波峰对应的频率将被用于对应轨道参数比值的计算，被识别的频率有2.5，9.4，21.8，24.9，32.3，42.6，47.8和57.8 周/Ma，其中频率2.5，9.4，24.9，42.6和47.8 周/Ma的倒数为其所对应的天文轨道周期，分别为0.405，0.100，0.040，0.023和0.020 Ma。以40.2~42.47 Ma时期的各天文周期范围为参考(图3)，结果表明对罗69井进行的旋回地层学分析是可靠的。对时间域GR曲线进行滑动窗口频谱分析，发现各峰值对应的频率大致与MTM功率谱(图8(c))中高置信度峰值相对应。

4  天文旋回控制下的湖相细粒沉积

4.1　自然伽马测井曲线的古气候和古环境意义

自然伽马测井(nature gamma-ray)是测量岩层中放射性元素(主要是⁴⁰K，²³²Th和²³⁸U)衰变后放射出来的伽马射线强度，相当于测量沉积地层中钾、钍、铀总含量的变化。不同种类岩石的放射性元素的含量存在明显差异，例如，泥岩中的黏土矿物和有机质对U的吸附能力较强^[44]，而由钙质生物壳体构成的碳酸盐岩，具有较低的放射性元素含量，促使生物有机体发育和生长^[45]。伽马曲线不仅能应用于区分基本岩性，目前已被广泛应用于反演古气候和古环境研究^[46-47]。在地层中GR值的升高在一定程度上能够表征泥质含量的增高，沉积环境水动力相对较弱，水体变深，海(湖)扩张等，从而反映相对降水多、湿润的古气候特征^[48]。刘泽纯等^[49]将碳酸盐岩含量、孢粉和δ¹⁸O等与GR曲线相结合作为古气候研究的替代指标，GR曲线和不同指标之间显示出相对一致的波动趋势及较好的可比性。

图8　罗69井时间域GR的频谱分析

Fig. 8　Spectral analysis of Well Luo 69 in time domain

自然伽马曲线在旋回地层学中被频繁使用^[50-52]。除去自然伽马对气候、环境变化相对敏感之外，GR数据序列最大的优势在于：油田的探井资料较容易获得，且采样间隔(0.125 m)具有高分辨率。本次研究以GR曲线为主进行旋回地层学分析，并通过分析GR数据与其他气候指标之间的相关性，找出影响气候变化的主要天文周期因素及其控制机理，为反演中国东部始新世气候变化提供参考。

4.2　基于米兰科维奇理论的细粒沉积物层序地层划分

米兰科维奇旋回与层序地层的定量划分有着紧密的联系。石巨业等^[53]研究表明，三级层序的持续时间一般为1~10 Ma，与构造因素相关；0.405 Ma的长偏心率旋回已经被证明与四级层序存在直接成因关系，具有高幅水深变化，且彼此有成因联系的水进-水退序列；五级层序受0.100 Ma的短偏心率控制，具有低幅水深特征；六级层序受0.040 Ma的斜率周期控制^[54]。

考虑到自然伽马曲线的意义，GR的极小值代表着湖平面当时的最低值，可作为四级层序界面的划分依据(图9)。综上所述，本次研究可将罗69井沙三下亚段地层划分为7个四级层序，其中Es³l¹和Es³l⁷发育不完整。

图9　时间域GR、TOC、黏土矿物质量分数、陆源碎屑指数、碳酸盐岩质量分数、化学风化强度指标数据序列与0.405 Ma理论轨道周期曲线的相位关系

Fig. 9　Time domained GR, TOC, clay mineral mass fraction, detrital classic index, carbonate mass fraction, chemical weathering index series and their relationships with astronomical theory period (0.405 Ma)

4.3　沙三下亚段细粒沉积物中有机质富集作用的天文轨道影响

图9中使用的其他古气候替代指标数据见文献[55]，包括罗69井沙三下亚段2 930~3 127 m之间的沉积地层所保存的古气候古环境信息记录。以本文所建立的罗69井的天文年代标尺为基础，将各类替代指标的深度数据对应到相对的地质时间上，得到时间域的数据序列。其中GR曲线的0.405 Ma滤波曲线(图9中红色实线)与TOC、黏土矿物质量分数、陆源碎屑指数、碳酸盐岩质量分数(图9蓝色实线为红色实线的镜像)均具有良好的对应关系，与化学风化强度不具有明显对应关系。

RACHOLD等^[56]总结了轨道强迫气候过程及其对沉积环境的影响，其中长偏心率周期(0.405 Ma)对长期季风周期的调节作用有重要影响^[57]，由于岁差旋回的最强和最弱振幅变化对应最大和最小长偏心率，从而形成了2个截然不同的季风气候^[58]。JIN等^[8]对古生代页岩地层中有机质富集规律建立2种假设，确定有机质富集作用主要受控于长偏心率轨道周期。由于太阳系的混沌行为，天文学家还无法获得超过250 Ma的精确天文轨道周期方案，因此，对影响古生代细粒沉积岩的季风变化控制沉积过程、沉积环境的探讨仅仅停留在2种模型的假设：第一种假设认为初级生产力是细粒沉积物中有机质富集(高TOC，高GR)的主要因素。偏心率极大值对应岁差幅度变化最剧烈，从而形成了强烈的季风系统，频繁的干湿气候变化及较强的季风降水作用。剧烈的风化作用使向湖盆中输送碎屑物质含有更多的养分，有利于湖盆内生命体的发育，即偏心率极大值控制有机制的富集，对应TOC、GR的极大值(图10(a))；第二种假设认为在长偏心率周期控制作用下，有机质的保存条件是细粒沉积物中有机质富集(高TOC，高GR)的主要因素。在偏心率极大值时，尽管大量的养分输入促进了季节性的高生产力，但在偏心率极小的干旱条件下，有机质又会被氧化分解，相当于并未真正富集。相反，在偏心率极小值时，季节的交替减少，提供了一个稳定的沉积环境，有机物不受保存条件的影响^[59]，可以大量富集(图10(b))。

本次研究获得的GR数据序列的滤波信号与理论轨道周期之间的相位关系在沙三下亚段早晚期呈不同状态(图9)：Es³l¹，Es³l²和Es³l³这3个四级层序(Es³l早期)的调谐GR曲线中，0.405 Ma滤波曲线的高值与E+O-P-0.405 Ma滤波曲线的低值相对应；而Es³l⁴，Es³l⁵，Es³l⁶和Es³l⁷这4个四级层序(Es³l晚期)的调谐GR曲线中，0.405 Ma滤波曲线的高值与E+O-P-0.405 Ma滤波曲线的高值相对应。即提出的模型二在罗69井沙三下亚段早期的沉积记录更可信，而模型一更适用于沙三下亚段晚期的沉积记录；沙三下亚段早期(Es³l¹，Es³l²，Es³l³)长偏心率极小值控制作用下的季节交替减少，稳定的沉积环境降低了有机质的氧化分解效率，在沙三下亚段晚期(Es³l⁴，Es³l⁵，Es³l⁶，Es³l⁷)长偏心率极大值控制作用下，强烈的季风作用导致碎屑物质提供更充足的养分，促进了有机生命体的生长。从图9所示的化学风化指数与偏心率滤波曲线并无明显相关性也可推测，模型一中剧烈风化作用为有机质富集提供养分的假设不成立。

图10　2种有机质富集因素模型

Fig. 10　Two factor models of organic accumulation

5  结论

1) 根据GR测井曲线频谱分析结果，渤海湾盆地沾化凹陷罗69井沙三下亚段地层中沉积旋回周期分别为0.405，0.095，0.040，0.023，0.022和0.019 Ma。其中长偏心率周期0.405 Ma控制的地层厚度为39.4 m；短偏心率周期控制的地层厚度为9.74 m；地轴斜率周期控制的地层厚度为3.94~2.54 m；岁差周期控制的地层厚度为1.92~1.70 m。罗69井沙三下亚段最优沉积速率为8.2×10³ cm/Ma。

2) 基于米兰科维奇旋回周期的识别结果，对去趋势化的自然伽马测井数据进行高通滤波处理，据0.405 Ma长偏心率旋回建立“浮动”天文年代标尺，并以沙三下亚段底部磁性地层年代数据42.47 Ma为锚点，沙三下亚段顶部磁性地层数据40.20 Ma为控制点，建立罗69井的天文年代标尺。

3) 以长偏心率0.405 Ma周期滤波曲线作为四级层序划分的参考曲线，将沙三下亚段划分为7个四级层序，实现湖相细粒沉积岩的高频地层单元划分，并精确计算沙三下亚段沉积速率的变化。

4) 从天文学角度验证了天文轨道周期变化和有机质富集因素之间的联系，并建立了2个有机质富集因素的对比模式：长偏心率控制初级生产力为主的有机质富集模式；长偏心率控制氧化还原条件为主的有机质富集模式。本研究所获得的GR数据序列的滤波信号与理论轨道周期的滤波信号比对结果，证实了沙三下亚段早期GR数据的极大值与长偏心率之间的极小值具有良好的对应关系，沙三下亚段晚期GR数据的极大值与长偏心率之间的极大值具有良好的对应关系，表明沙三下亚段中有机质的富集存在保存条件和生产力条件这2个不同的主控因素。

参考文献：

[1] MILANKOVITCH M M. Kanon der Erdbestrahlung und seine Anwendung auf das Eiszeitenproblem[M]. Académie royale serbe: Editions speciales, 1941: 133.

[2] EINSELE G. Cycles and events in stratigraphy[M]. Berlin: springer-verlag, 1991: 7-20.

[3] SCHWARZACHER W. Cyclostratigraphy and the Milankovitch theory[M]. Amsterdam: Elsevier Science, 1993: 12-25.

[4] KODAMA K P, HINNOV L A. Rock magnetic cyclostratigraphy[M]. Oxford: John Wiley & Sons, Ltd, 2014: 157-160.

[5] 黄春菊. 旋回地层学和天文年代学及其在中生代的研究现状[J]. 地学前缘, 2014, 21(2): 48-66.

HUANG Chunju. The current status of cyclostratigraphy and astrochronology in the Mesozoic[J]. Earth Science Frontiers, 2014, 21(2): 48-66.

[6] STRASSER A H, HECKEL P H. Cyclostratigraphy concepts, definitions, and applications[J]. Newsletters on Stratigraphy, 2006, 42(2): 75-114.

[7] WEEDON G P. Time-series analysis and cyclostratigraphy[M]. Cambridge: Cambridge University Press, 2003: 77-78.

[8] JIN Siding, DENG Hucheng, ZHU Xing, et al. Orbital control on cyclical organic matter accumulation in Early Silurian Longmaxi Formation shales[J]. Geoscience Frontiers, 2020, 11(2): 533-545.

[9] 房吉闯. 湖北宜昌晚奥陶世—早志留世海相地层旋回地层学研究[D]. 北京: 中国地质大学(北京)地球科学学院, 2017: 14-18.

FANG Jichuang. Cyclostratigraphy research of Late Ordovician: Early Silurian in Hubei Province[D]. Beijing: China University of Geosciences(Beijing). School of Earth Sciences, 2017: 14-18.

[10] 李明松. 华南早—中三叠世天文年代学与古环境变化的天文驱动力[D]. 武汉: 中国地质大学(武汉)资源学院, 2016: 22-26.

LI Mingsong. Astrochronology of the Early-Middle Triassic in South China and astronomical forced paleoenvironment change[D]. Wuhan: China University of Geosciences(Wuhan). College of Resources, 2016: 22-26.

[11] 孙善勇, 刘惠民, 操应长, 等. 湖相深水细粒沉积岩米兰科维奇旋回及其页岩油勘探意义: 以东营凹陷牛页1井沙四上亚段为例[J]. 中国矿业大学学报, 2017, 46(4): 846-858.

SUN Shanyong, LIU Huimin, CAO Yingchang, et al. Milankovitch cycle of lacustrine deepwater fine-grained sedimentary rocks and its significance to shale oil: a case study of the upper Es4 member of well NY₁ in Dongying sag[J]. Journal of China University of Mining & Technology, 2017, 46(4): 846-858.

[12] CAO Haiyang, JIN Siding, SUN Ming, et al. Astronomical forcing of sedimentary cycles of late Eocene Liushagang formation in the Bailian sag, Fushan depression, Beibuwan basin, South China sea[J]. Journal of Central South University, 2016, 23(6): 1427-1438.

[13] LOUCKS R G, RUPPEL S C. Mississippian Barnett Shale: lithofacies and depositional setting of a deep-water shale-gas succession in the Fort Worth Basin, Texas[J]. AAPG Bulletin, 2007, 91(4): 579-601.

[14] GAMBACORTA G, MENICHETTI E, TRINCIANTI E, et al. Orbital control on cyclical primary productivity and benthic anoxia: astronomical tuning of the Telychian Stage (Early Silurian)[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2018, 495: 152-162.

[15] WEEDON G P, COE A L, GALLOIS R W. Cyclostratigraphy, orbital tuning and inferred productivity for the type Kimmeridge Clay(Late Jurassic), Southern England[J]. Journal of the Geological Society, 2004, 161(4): 655-666.

[16] MACHLUS M L, OLSEN P E, CHRISTIE-BLICK N, et al. Spectral analysis of the lower Eocene Wilkins peak member, green river formation, Wyoming: support for milankovitch cyclicity[J]. Earth and Planetary Science Letters, 2008, 268(1/2): 64-75.

[17] ABELS H A, AZIZ H A, VENTRA D, et al. Orbital climate forcing in mudflat to marginal lacustrine deposits in the Miocene Teruel basin (northeast Spain)[J]. Journal of Sedimentary Research, 2009, 79(11): 831-847.

[18] SZURLIES M. Latest Permian to Middle Triassic cyclo-magnetostratigraphy from the Central European Basin, Germany: implications for the geomagnetic polarity timescale[J]. Earth and Planetary Science Letters, 2007, 261(3/4): 602-619.

[19] WU Huaichun, ZHANG Shihong, JIANG Ganqing, et al. Astrochronology of the Early Turonian-Early Campanian terrestrial succession in the Songliao Basin, northeastern China and its implication for long-period behavior of the solar system[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2013, 385: 55-70.

[20] ZHANG Rui, JIN Zhijun, LIU Quanyou, et al. Astronomical constraints on deposition of the Middle Triassic Chang 7 lacustrine shales in the Ordos Basin, Central China[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2019, 528: 87-98.

[21] SHI Juye, JIN Zhijun, LIU Quanyou, et al. Terrestrial sedimentary responses to astronomically forced climate changes during the Early Paleogene in the Bohai Bay Basin, Eastern China[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2018, 502: 1-12.

[22] LIU Zhanhong, HUANG Chunju, ALGEO T J, et al. High-resolution astrochronological record for the Paleocene-Oligocene (66-23 Ma) from the rapidly subsiding Bohai Bay Basin, northeastern China[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2018, 510: 78-92.

[23] ZHAO Ke, DU Xuebin, LU Yongchao, et al. Are light-dark coupled laminae in lacustrine shale seasonally controlled a case study using astronomical tuning from 42.2 to 45.4 Ma in the Dongying Depression, Bohai Bay Basin, Eastern China[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2019, 528: 35-49.

[24] COOPER M C, O'SULLIVAN P E, SHINE A J. Climate and solar variability recorded in Holocene laminated sediments: a preliminary assessment[J]. Quaternary International, 2000, 68/69/70/71: 363-371.

[25] TUENTER E, WEBER S L, HILGEN F J, et al. The response of the African summer monsoon to remote and local forcing due to precession and obliquity[J]. Global and Planetary Change, 2003, 36(4): 219-235.

[26] JIN Siding, CAO Haiyang, WANG Hua, et al. The Paleogene multi-phase tectono-sedimentary evolution of the syn-rift stage in the Nanpu Sag, Bohai Bay Basin, East China[J]. Energy Exploration & Exploitation, 2018, 36(6): 1519-1545.

[27] FENG Youliang, LI Sitian, LU Yongchao. Sequence stratigraphy and architectural variability in Late Eocene lacustrine strata of the Dongying Depression, Bohai Bay Basin, Eastern China[J]. Sedimentary Geology, 2013, 295: 1-26.

[28] 彭丽, 伍轶鸣, 练章贵, 等. 陆相断陷湖盆高频层序特征及其沉积演化: 以渤海湾盆地济阳坳陷沙三下亚段为例[J]. 石油与天然气地质, 2019, 40(4): 789-798.

PENG Li, WU Yiming, LIAN Zhanggui, et al. Features and sedimentary evolution of high-frequency sequence in continental lacustrine rift basin: example of the lower Shahejie member 3 in Jiyang Depression, Bohai Bay Basin[J]. Oil & Gas Geology, 2019, 40(4): 789-798.

[29] 翟光明, 何文渊. 渤海湾盆地资源潜力和进一步勘探方向的探讨[J]. 石油学报, 2002, 23(1): 1-5.

ZHAI Guangming, HE Wenyuan. Potential resources and future exploration orientation in Bohai bay basin[J]. Acta Petrolei Sinica, 2002, 23(1): 1-5.

[30] 贾屾, 姜在兴, 张文昭. 沾化凹陷页岩油储层特征及控制因素[J]. 海洋地质前沿, 2018, 34(12): 29-38.

JIA Shen, JIANG Zaixing, ZHANG Wenzhao. Characteristics of the shale oil reservoir in Zhanhua sag and its controlling factors[J]. Marine Geology Frontiers, 2018, 34(12): 29-38.

[31] 李志明, 陶国亮, 黎茂稳, 等. 渤海湾盆地济阳坳陷沾化凹陷L69井古近系沙三下亚段取心段页岩油勘探有利层段[J]. 石油与天然气地质, 2019, 40(2): 236-247.

LI Zhiming, TAO Guoliang, LI Maowen, et al. Favorable interval for shale oil prospecting in coring well L69 in the Paleogene Es^3L in Zhanhua sag, Jiyang depression, Bohai bay basin[J]. Oil & Gas Geology, 2019, 40(2): 236-247.

[32] LI Mingsong, HINNOV L, KUMP L. Acycle: time-series analysis software for paleoclimate research and education[J]. Computers & Geosciences, 2019, 127: 12-22.

[33] CLEVELAND W S. Robust locally weighted regression and smoothing scatterplots[J]. Journal of the American Statistical Association, 1979, 74(368): 829-836.

[34] PAILLARD D, DONNADIEU Y. A 100 Myr history of the carbon cycle based on the 400 kyr cycle in marine δ¹³C benthic records[J]. Paleoceanography, 2014, 29(12): 1249-1255.

[35] GHIL M, ALLEN M R, DETTINGER M D, et al. Advanced spectral methods for climatic time series[J]. Reviews of Geophysics, 2002, 40(1): 3-1-3-41.

[36] MEYERS S R. Astrochron: an R Package for Astrochronology. Version 0.7. 2014[EB/OL]. [2021-02-10]. https: //cran.r-project.org/web/packages/astrochron.

[37] HAMMER O, HARPER D A T, RYAN P D. Past: paleontological statistics software package for education and data analysis [J]. Palaeontologia Electronica, 2001, 4(1): 1-9.

[38] LASKAR J, ROBUTEL P, JOUTEL F, et al. A long-term numerical solution for the insolation quantities of the Earth[J]. Astronomy & Astrophysics, 2004, 428(1): 261-285.

[39] FIENGA A, LASKAR J, MORLEY T, et al. INPOP08, a 4-D planetary ephemeris: from asteroid and time-scale computations to ESA Mars Express and Venus Express contributions[J]. Astronomy & Astrophysics, 2009, 507(3): 1675-1686.

[40] 吴怀春, 张世红, 黄清华. 中国东北松辽盆地晚白垩世青山口组浮动天文年代标尺的建立[J]. 地学前缘, 2008, 15(4): 159-169.

WU Huaichun, ZHANG Shihong, HUANG Qinghua. Establishment of floating astronomical time scale for the terrestrial Late Cretaceous Qingshankou formation in the Songliao basin of Northeast China[J]. Earth Science Frontiers, 2008, 15(4): 159-169.

[41] CROWLEY T J, KIM K Y, MENGEL J G, et al. Modeling 100 000-year climate fluctuations in pre-Pleistocene time series[J]. Science, 1992, 255(5045): 705-707.

[42] MEYERS S R, SAGEMAN B B, HINNOV L A. Integrated quantitative stratigraphy of the cenomanian-turonian bridge creek limestone member using evolutive harmonic analysis and stratigraphic modeling[J]. Journal of Sedimentary Research, 2001, 71(4): 628-644.

[43] MA C, MEYERS S R, SAGEMAN B B, et al. Testing the astronomical time scale for oceanic anoxic event 2, and its extension into Cenomanian strata of the Western Interior Basin (USA)[J]. Geological Society of America Bulletin, 2014, 126(7/8): 974-989.

[44] 杨平, 陈晔, 刘泽纯. 柴达木盆地自然伽玛曲线在古气候及沉积环境研究中的应用[J]. 古地理学报, 2003, 5(1): 94-102.

YANG Ping, CHEN Ye, LIU Zechun. Application of gamma ray log to study on palaeoclimate and sedimentary environments of the Jurassic in Qaidam basin[J]. Journal of Palaeogeography, 2003, 5(1): 94-102.

[45] 梁文君, 肖传桃, 代城建, 等. 自然伽马曲线应用于古气候、古环境研究: 以柴达木盆地七个泉地区古近—新近纪地层为例[J]. 矿产与地质, 2015, 29(5): 698-702.

LIANG Wenjun, XIAO Chuantao, DAI Chengjian, et al. A study on application of natural gamma ray curve in paleoclimate and paleoenvironment: a case study of Paleogene-Neogene stratain Qigequan area of Qaidam Basin[J]. Mineral Resources and Geology, 2015, 29(5): 698-702.

[46] RUFFELL A, WORDEN R. Palaeoclimate analysis using spectral gamma-ray data from the Aptian(Cretaceous) of southern England and southern France[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2000, 155(3/4): 265-283.

[47] GHASEMI-NEJAD E, ARDAKANI E P, RUFFELL A. Palaeoclimate change recorded in Upper Cretaceous(albian-cenomanian) kazhdumi formation borehole SPECTRAL gamma-ray logs, south pars gas field, Persian gulf[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2010, 291(3/4): 338-347.

[48] BAKER P A, RIGSBY C A, SELTZER G O, et al. Tropical climate changes at millennial and orbital timescales on the Bolivian Altiplano[J]. Nature, 2001, 409(6821): 698-701.

[49] 刘泽纯, 陈晔, 袁林旺, 等. 应用自然伽玛测井曲线反演 2.85 Ma BP来古气候变化[J]. 中国科学(D辑), 2000, 30(6): 609-618.

LIU Zechun, CHEN Ye, YUAN Linwang, et al. Inversion the palaeoclimate change of 2.85 Ma BP by the natural gamma ray log[J]. China Sciences (D Series), 2000, 30(6): 609-618.

[50] TENVEEN J H , POSTMA G. Astronomically forced variations in gamma-ray intensity: Late Miocene hemipelagic successions in the eastern Mediterranean basin as a test case[J]. Geology, 1996, 24(1): 15-18.

[51] 刘占红, 陈荣, 宋成兵, 等. 塔里木盆地石炭系卡拉沙依组旋回地层与层序地层综合研究[J]. 地球科学, 2012, 37(5): 1051-1065.

LIU Zhanhong, CHEN Rong, SONG Chengbing, et al. Cyclostratigraphy and sequence stratigraphy on carboniferous karashayi formation of Tarim basin, China[J]. Earth Science, 2012, 37(5): 1051-1065.

[52] JIN Siding, CAO Haiyang, WANG Hua, et al. Orbital cyclicity in sedimentary sequence and climatic indications of C-O isotopes from Lower Cretaceous in Qingxi Sag, Jiuquan Basin, NW China[J]. Geoscience Frontiers, 2019, 10(2): 467-479.

[53] 石巨业, 金之钧, 刘全有, 等. 基于米兰科维奇理论的湖相细粒沉积岩高频层序定量划分[J]. 石油与天然气地质, 2019, 40(6): 1205-1214.

SHI Juye, JIN Zhijun, LIU Quanyou, et al. Quantitative classification of high-frequency sequences in fine-grained lacustrine sedimentary rocks based on Milankovitch theory[J]. Oil & Gas Geology, 2019, 40(6): 1205-1214.

[54] 龚一鸣, 杜远生, 童金南, 等. 旋回地层学: 地层学解读时间的第三里程碑[J]. 地球科学, 2008, 33(4): 443-457.

GONG Yiming, DU Yuansheng, TONG Jinnan, et al. Cyclostratigraphy: the third milestone of stratigraphy in understanding time[J]. Earth Science, 2008, 33(4): 443-457.

[55] MA Yiquan, FAN Majie, LU Yongchao, et al. Climate-driven paleolimnological change controls lacustrine mudstone depositional process and organic matter accumulation: Constraints from lithofacies and geochemical studies in the Zhanhua Depression, Eastern China[J]. International Journal of Coal Geology, 2016, 167(1): 103-118.

[56] RACHOLD V, BRUMSACK H J. Inorganic geochemistry of Albian sediments from the Lower Saxony Basin NW Germany: palaeoenvironmental constraints and orbital cycles[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2001, 174(1/2/3): 121-143.

[57] WANG Pinxian. Global monsoon in a geological perspective[J]. Chinese Science Bulletin, 2009, 54(7): 1113-1136.

[58] ZACHOS J C, MCCARREN H, MURPHY B, et al. Tempo and scale of late Paleocene and early Eocene carbon isotope cycles: implications for the origin of hyperthermals[J]. Earth and Planetary Science Letters, 2010, 299(1/2): 242-249.

[59] HERBERT T D, FISCHER A G. Milankovitch climatic origin of mid-Cretaceous black shale rhythms in central Italy[J]. Nature, 1986, 321(6072): 739-743.

(编辑  刘锦伟)

收稿日期： 2020 -09 -30; 修回日期： 2020 -11 -09

基金项目(Foundation item)：国家自然科学基金资助项目(41802030); 国家科技重大专项(2017ZX05036004006) (Project (41802030) supported by the National Natural Science Foundation of China; Project(2017ZX05036004006) supported by the National Science and Technology Major Program)

通信作者：金思丁，博士，副教授，从事旋回地层学研究；E-mail：jinsiding@cdut.edu.cn

引用格式： 张若琳, 金思丁. 渤海湾盆地沾化凹陷罗69井沙三下亚段旋回地层学研究[J]. 中南大学学报(自然科学版), 2021, 52(5): 1516-1531.

Citation: ZHANG Ruolin, JIN Siding. Cyclostratigraphy research on lower Member 3 of Shahejie Formation in Well Luo 69 in Zhanhua Sag Bohai Bay Basin[J]. Journal of Central South University(Science and Technology), 2021, 52(5): 1516-1531.