下载:
-
四川盆地筇竹寺组是寒武系页岩气的重要勘探目的层[1⁃2],是继志留系龙马溪组之后的又一页岩气勘探开发的热点层段,具有较好的勘探前景[3⁃4]。前人针对该层段沉积、层序地层方面做过诸多工作[5⁃8],但由于寒武系地层时代老,取心长度有限,区域上也缺乏可靠的古生物地层或年代地层等资料的约束,筇竹寺组地层的沉积旋回与充填过程至今尚不十分明朗,很大程度上制约了该层系页岩气的勘探进程。大量学者尝试对四川盆地及周缘地区的寒武纪地层进行了层序地层划分[9⁃11]。大多数划分方案都以典型露头剖面实测为基础,结合岩石薄片、钻测井、古生物、地球化学和地球物理等资料,通过对层序界面特征和层序内部充填样式的识别,建立年代地层、岩石地层与层序地层之间的联系。然而不同学者的层序划分方案和标准并未达成一致,尤其是没有利用统一的时间区间对层序进行分级。三级层序海平面升降的起源是过去30年悬而未决的问题之一,目前被广泛接受的一种级别划分方式是将三级层序的形成归因为冰川型海平面变化的结果,并将三级层序的持续时限与海平面变化的周期相对应[12⁃13]。因沉积相类型多且变化快,无法通过某一个或几个剖面去追踪区内所有三级层序界面,前人对下寒武统泥页岩沉积地层的三级层序的划分存在一定不合理性。近年来,基于米兰科维奇理论的高频层序识别与划分已成为研究的热点和前沿领域[14⁃15],不仅成功应用于厘定地层沉积持续时间及判别重要地质事件的天文影响因素[16⁃19],并且通过提取天文周期曲线为参考曲线的方式,实现了深水页岩沉积的高频层序定量划分[20⁃22],该理论方法可以从成因上解决高频沉积层序单元的划分问题[23]。在解决页岩地层精细划分和对比的问题基础之上,可进一步实现页岩油气勘探有利层段的预测,对页岩油气勘探也具有重要意义。
国内在20世纪90年代初就开始利用测井资料进行旋回地层学研究[24],目前常规测井曲线（自然伽马、电阻率、密度、声波时差等）是最常使用的测井资料[25⁃27]。其中,自然伽马测井综合反映了沉积地层中放射性元素（钍、铀、钾）的含量,能够较好地反映岩性变化,尤其是地层中黏土矿物和有机质含量的变化,常被用作古气候替代指标[20,26,28]；电阻率若仅受黏土矿物含量或有机质含量其中之一的单因素调控,亦能成为有效的古气候替代指标[29],例如：美国西部内陆盆地白垩纪Niobrara组地层电阻率被认为与天文轨道驱动作用下的碳酸盐岩含量相关[30]；位于佛罗里达海峡附近的钻探（航次166,站位1006）的声波曲线记录了1.4 m的天文强迫周期[31]。天文轨道周期驱动作用下的气候变化能对不同沉积组分产生相应的影响,且不同测井参数对古气候的响应方式不同[29,32],单一测井替代指标及划分标准不能同时适用于所有地层。
本次研究拟对四川盆地井研—犍为地区下寒武统筇竹寺组沉积地层开展旋回地层学研究。在考虑测井响应主控地质因素的前提下,以常规测井数据作为替代指标,并对测井数据进行多窗口频谱分析和相关系数分析,从地层中识别出可靠的米兰科维奇信号。通过长偏心率滤波和天文校准,建立金页1井的“浮动”天文年代标尺。根据沉积噪音模型DYNOT曲线,从相对海平面变化周期的角度探讨筇竹寺组三级层序的划分方案,并根据已建立的天文标尺,估算各个三级层序的持续时间。最后,尝试以邻区梅树村剖面石岩头组底部CA-ID-TIMS U-Pb定年结果（526.86±0.16 Ma）为年龄控制点,估算金页1井时间序列关键界面的年龄,并与2020年国际地质年代表进行对比,验证本次旋回地层研究的准确性。
-
四川盆地位于上扬子板块,是一个四面被造山带环绕的菱形叠合盆地,总面积为18×104 km2（图1a）。震旦纪灯影阶至早寒武世麦地坪阶,上扬子地区发生了三幕桐湾运动[33],区域性隆升和间断性剥蚀作用使得岩相古地理极其复杂,由北西向南东的沉积环境展布为：古陆—滨海相—浅水陆棚相—深水陆棚相[34⁃35]。查明不同地区富有机质页岩段的等时对比关系是制约四川盆地页岩气勘探选区的关键地质问题,刘忠宝等[9]将上扬子下寒武统整体划分为8个三级层序,其中筇竹寺阶发育4个三级层序（图1b）。大量研究表明,四川盆地绵阳—长宁拉张裂陷槽在早寒武世充填了巨厚的深水陆棚相泥页岩[36],受裂陷槽影响,筇竹寺组泥页岩在资阳—宜宾一带地层横向发育稳定,无沉积间断且厚度大,较裂陷槽外部浅水陆棚及台地区更适宜开展旋回地层学工作。
图 1 （a）四川盆地及周缘地区下寒武统筇竹寺组岩相古地理平面图；（b）四川盆地下寒武统层序地层对比图
-
与全球寒武纪年代地层划分方案相对比,国内传统划分的下寒武统下部相当于纽芬兰统和第二统[37]。四川盆地及其周缘地区下寒武统自下而上划分为梅树村阶、筇竹寺阶、沧浪铺阶与龙王庙阶,其中四川盆地麦地坪组（下寒武统底部）大致对应的国际年代地层单位为纽芬兰统幸运阶与第二阶,筇竹寺组大致对应第二阶与第三阶[38⁃39]。四川盆地内该套沉积地层具有同时异相性,地层序列的时空变化极为复杂,岩石地层名称复杂且不统一[40]（表1）。金页1井在寒武系筇竹寺组沉积期时处于裂陷槽沉积中心附近,已经有许多学者对金页1井进行了较为丰富的层序、沉积研究[4⁃6,44⁃45],为我们的研究提供了宝贵的基础资料。通过前人对四川盆地下寒武地层的时间框架的搭建,可以初步判定金页1井筇竹寺组沉积于~529 Ma至~521 Ma之间（表1）。同时可将该钻井与东北部发表了丰富同位素和生物地层的肖滩剖面进行横向对比[41⁃43]（图2a,b）。肖滩剖面石岩头组下段存在的一段黑色岩性标志层,在整个华南板块内部台地区域可以进行对比,金页1井筇竹寺组下段同样也发育这样一套黑色页岩（图2c）,表明筇竹寺组底部与石岩头组底部位于同一沉积时期；金页1井麦地坪组白云岩和磷块岩可与肖滩剖面朱家箐组中谊村段的特征岩性进行对比,麦地坪组与筇竹寺组不整合接触,四川盆地麦地坪组与肖滩剖面寒武系早期朱家箐组处于同一沉积时期[52]。
-
本文的研究数据主要为金页1井的常规测井和自然伽马能谱测井,除此之外,还包括TOC含量和黏土矿物含量的实测数据及测井曲线拟合数据。该钻井以寒武系筇竹寺组页岩层系为目的层,钻至下伏灯影组,本文研究主要选取3 184~3 655 m,采样间隔为0.1 m的测井数据进行旋回地层学分析。在进行天文旋回分析之前,对各类测井数据采用滑动窗口相关分析（evolutionary correlation analysis）,可根据其所显示的平面二维热力图,对测井曲线数值的相关程度以不同深度及不同地层厚度变化的角度进行了解[29],并通过测井响应特征分析,选取最佳测井序列作为旋回地层学的替代指标。
旋回地层学分析方法主要是时间序列分析方法,包括Multitaper method（MTM）频谱分析[53]、Fast Fourier transform（LAH）滑动窗口频谱分析[54⁃55]、滤波分析、相关系数分析（correlation coefficient及evolutionary correlation coefficient）、沉积噪音模型（sedimentary noise model）分析等[56]。由于曲线包含各类地质因素引起的周期性变化,除了一些与地质因素无关的干扰信号之外,还包括一些长周期引起的高频信号,在进行频谱分析之前,要将这种信号通过去趋势化的方法去除[57],去趋势化（detrended）的测井序列中受天文周期影响的信号将更突出。本次研究的所有计算使用的软件包为基于Matlab平台下的Acycle V2.4[58]。
-
测井响应特征分析是利用测井资料进行地质分析的基础。古环境、古气候的变化对不同测井数据影响的差异,决定了不同测井资料是否能够解决相关古气候古环境问题。自然伽马（GR）测井的数值主要取决于地层中自然伽马辐射的强度,主要受地层中Th、U、K元素的分布与聚集的影响。K和Th的含量与黏土矿物含量关系密切,在风化和黏土矿物形成过程中,Th元素具有有限的流动性,比K更不容易溶解,完全存在于碎屑组分内。前人研究表明,Th和K的比值可以较好地反映碎屑颗粒的粒度大小,从而指示当时水体环境的能量大小[59]。Th和U的比值能够反映沉积水体的氧化还原条件[60],随着水体还原程度的逐渐增加,Th/U值逐渐减小。一般而言,Th/U小于4为强还原环境,大于4且小于10指示还原环境,大于10指示氧化环境[61]。在具有很高的U值,而Th、K数值较低的层段,随U值增高,Th、K数值并未相应增大,可能是流体滤液十分活跃所离析的铀矿物附着于裂缝中所引起的。另外,U含量与氧化还原条件密切相关,还原环境下,有机质含量高,对应的U含量高。
金页1井筇竹寺组（3 184~3 566 m）的GR、U元素测井曲线具有相似的变化特征,而Th元素和K元素测井曲线变化趋势则与GR、U元素曲线明显不同（图3）,推测U元素主要形成于后生成岩作用,自生铀在还原条件下的富集作用引起部分层位出现的高U值。与此同时,滑动窗口相关分析图显示（图4）,GR曲线与Th、K元素含量的相关性整体上是正相关,但在较大部分的深度段相关系数都低于0.5,如3 303 m,3 384 m,3 496 m深度附近,滑动窗口大小20 m左右的热力图中相关系数的读数约为0.2,不具有相关性。而GR曲线与U元素含量整体上不具有明显的相关性,相关系数在热力图中的读数为0.2左右,但在3 303 m,3 402 m,3 566 m深度附近,相关系数数值接近1.0,显示出强烈的正相关,推测是地层中放射性元素的非气候成因富集导致的异常高GR和高U值。因此,GR和U曲线并不能很好地反映同沉积时期的沉积环境,不是天文旋回研究的有效古气候替代指标。由于Th和K元素都和碎屑黏土密切相关,所以黏土矿物含量与Th、K元素含量显示出明显的正相关性,从另一个角度证明自然伽马能谱测井曲线中的Th、K元素含量能够很好地反映古气候、古环境变化,可以作为金页1井沉积地层的天文旋回研究的良好替代指标。Th/K与黏土矿物含量的相关性,不如Th、K元素与黏土矿物含量的相关性显著（图4）,部分层段表现为负相关（3 303 m,3 566 m）,这可能是由于K元素含量对黏土矿物含量的影响要高于Th元素。
图 3 金页1井筇竹寺组（3 184~3 566 m）测井曲线特征
图 4 金页1井筇竹寺组伽马能谱测井相关分析
由于Th元素和K元素对古气候的响应机理类似,在偏心率较大的地质历史时期,相对温暖湿润的气候条件增加降水、化学风化强度等,导致沉积颗粒粒径减小,黏土矿物沉积量增加；升高的黏土矿物含量通常对应于的高Th和高K值,高Th和K值被认为与高偏心率的地质历史时期对应,本文采用K元素含量作为旋回地层学的替代指标。在3 303 m,3 384 m,3 496 m等深度段测井数据均出现了较为明显的变化（图3,4）。为了分析不同沉积环境的测井响应特征,分别对这四个深度段进行测井值统计分析（图5）。可以看出在筇竹寺组的四段地层中,各类测井曲线的均值和方差都发生了明显的变化,指示出筇竹寺组页岩沉积环境及沉积速率在这四段沉积地层中发生了明显变化。
图 5 金页1井筇竹寺组不同层段伽马能谱测井响应特征统计
-
由于太阳系内的行星混沌作用贯穿始终,50 Ma以前的地层并没有高精度的天文年代学校正[62]。寒武纪早期的天文轨道周期不同于运用在比较新的地层中20∶5∶2∶1的周期比值,斜率和岁差会更接近,部分叠置难以识别；且短偏心率周期与斜率周期也可呈4∶1的比值关系,与长偏心率与短偏心率4∶1的比值相同,而造成多解性。目前有三种轨道周期方案试图通过建立相对准确的地月系统演化模型来确定较老的地层的天文周期信号[63⁃64]。本次研究采用Waltham[65]解决方案,认为405 kyr长偏心率周期在地质历史时期内保持稳定不变,以同时代古潮汐韵律计地球自转速度建立月球撤退模型计算出520 Ma左右的天文轨道周期为405 kyr,131 kyr,95 kyr,36.4 kyr,28.2 kyr,20.4 kyr,19.4 kyr,16.8 kyr。
-
未做处理的K元素测井曲线存在随深度增大先减小,随后又增加的趋势（图3）,这可能是长周期的构造演化所产生的信号,也有可能是测井曲线在原位地层中受地层温度、压力的影响所导致的与天文周期无关的旋回信号。前人研究表明,该沉积地层的沉积速率在2~5 cm/kyr[66],长偏心率周期（405 kyr）控制的沉积旋回厚度为8~20 m,因此本次研究中去趋势化处理采用的窗口大小为30 m,用于保留天文轨道长偏心率周期。
对趋势化后的K元素测井数据进行滑动窗口频谱分析（图6a）,滑动窗口步长为0.7 m,窗口长度为30 m。图谱中峰值所对应的深度（频率的倒数）表明各个深度段均保存了长偏心率、短偏心率及斜率的信号,岁差信号不明显。然而,谱图中峰值所对应的频率在3 303 m,3 384 m和3 496 m处同样出现了明显的峰值信号不连续现象（图6a）,是由于某些深度段的测井序列中包含更多的噪声,高频信号被具有更高功率的噪声所抑制[67],也可能是由于沉积速率的变化导致沉积周期发生了变化[68]。综合考虑滑动频谱图谱上出现信号不连续的深度（图6a）,及测井响应特征（图4,5）发生变化的深度,将整个筇竹寺组（3 184~3 566 m）划分为四个层段,从下到上依次为：一段（3 496~3 566 m）,二段（3 384~3 496 m）,三段（3 303~3 384 m）和四段（3 184~3 303 m）（图6b~e）,并分别进行频谱分析。在筇竹寺组一段地层中（图6b）,高于90%置信度的6个频率谱峰对应的厚度周期为10.44,3.68,2.53,1.34,1.08和1.01 m；筇竹寺组二段地层（图6c）,高于90%置信度的6个频率谱峰对应的厚度周期为18.86,6.25,4.52,2.43,1.70和1.70 m；筇竹寺组三段地层（图6d）,高于90%置信度的7个频率谱峰对应的厚度周期为11.76,5.71,3.41,2.34,1.21和0.87 m；筇竹寺组四段地层（图6e）,高于90%置信度的7个频率谱峰对应的厚度周期为13.33,4.64,3.21,2.51,1.40和1.01 m。筇竹寺组内划分的四个层段识别的部分显著周期均和405∶131∶95∶36.4∶28.2的比值相接近,证明筇竹寺组四段地层沉积时期均受到米兰科维奇旋回的影响。由于沉积环境在不同层段发生了改变,沉积速率发生了微小变化,此处通过长偏心率周期（405 kyr）对应的沉积厚度,对每段沉积速率进行一个初步估计：一段沉积速率约为2.46 cm/kyr,二段沉积速率约为4.65 cm/kyr,三段沉积速率约为2.90 cm/kyr,四段沉积速率约为3.28 cm/kyr。金页1井筇竹寺组的沉积速率呈现先增大后减小,随后又增大的趋势。
图 6 金页1井筇竹寺组（3 184~3 566 m）K元素曲线频谱分析图
-
通过频谱分析,在90%的置信度之上检测到了周期性变化的功率,但所选取的K元素测井曲线的数据是以深度为单位,在对应的沉积时间范围内并没有进行频率计算。因此,本文通过“COCO（correlation coefficient）”和“eCOCO（evolution correlation coeffi cient）”[57]两种方法分析,将深度序列的替代指标转换成时间序列,然后计算该时间序列的功率谱与天文解决方案功率谱之间的相关系数,从而提高旋回识别的准确性（图7）。本次研究的“COCO”计算得出金页1井筇竹寺组的平均沉积速率在2.5,2.9,3.4和4.8 cm/kyr处的H0显著水平（被错误拒绝的零假设,即非天文轨道驱动信号）都低于0.001,且参与计算的天文轨道参数分量多于6个。本次研究的eCOCO分析选取1~6 cm/kyr的沉积速率区间,滑动步长为0.5 m,滑动窗口为30 m,蒙特卡洛模拟2 000次,得到图7所示结果,其中红色（极值）决定了特定深度的最优沉积速率。
图 7 金页1井筇竹寺组深度序列（3 184~3 566 m）的COCO分析和eCOCO分析结果
根据频谱分析,用比值法对不同深度段的平均沉积速率的初步估算（2.46 cm/kyr,4.65 cm/kyr,2.90 cm/kyr和3.28 cm/kyr）,与用检验法（eCOCO）对最优沉积速率的评价的结果相一致,再次证明本次研究选取的功率谱频段组合是满足天文轨道驱动条件的。
-
基于频谱分析及COCO、eCOCO分析的结果,将代表长偏心率（E,405 kyr）的沉积旋回通过高斯带通滤波提取出来。其中,筇竹寺组一段的滤波频率为0.10±0.02 cycles/m,二段的滤波频率为0.053±0.009 cycles/m,三段的滤波频率为0.085±0.015cycles/m,四段的滤波频率为0.075±0.001 cycles/m（图8）。滤波结果显示金页1井筇竹寺组总共记录了29个长偏心率周期,其中：一段（3 496~3 566 m）记录了约7个长偏心率旋回；二段（3 384~3 496 m）记录了6个长偏心率旋回；三段（3 303~3 384 m）记录了7个长偏心率旋回；四段（3 184~3 303 m）记录了9个长偏心率旋回（图8c）。本次研究选取代表长偏心率的滤波曲线建立深度—时间模型,从而得到时间域数据序列,并结合对应的深度域序列获得金页1井筇竹寺组连续变化的沉积速率曲线（图8d）。在3 303 m,3 384 m和3 496 m处沉积速率出现了明显的变化,整个筇竹寺组沉积时期,由下到上呈现先增大后减小的趋势。
图 8 金页1井筇竹寺组K元素曲线滤波、调谐、时间域滑动窗口频谱分析和浮动天文年代标尺
-
以筇竹寺组底部年龄526.86±0.16 Ma作为天文年代调谐的初始锚点,界定本次研究的时间范围,建立“浮动天文年代标尺”（图8h）。为了验证滤波和调谐结果的准确性,对用天文校准后转化成时间域的数据序列再次进行滑动窗口频谱分析（图8g）,结果显示识别的频率有0.002 5,0.007 6,0.010 5,0.027 8,0.038 5 cycles/kyr,在405 kyr（长偏心率期）,131 kyr（短偏心率周期）,95 kyr（短偏心率周期）,36 kyr（斜率周期）,28 kyr（斜率周期）的位置上有较强的功率谱,验证了天文调谐的结果是可靠的。在岁差频率上由于沉积速率较低和数据分辨率不够等问题,没有获得较明显的天文信号。
-
传统的层序地层学研究多应用于盆地边缘的沉积序列中,基于等时的角度研究地层叠加样式及其变化[69],层序地层学的发展极大地完善了与海平面变化相关的沉积物堆砌样式的解释,并通过建立等时地层格架,来重建地质历史时期的海平面变化[70⁃71]。然而,对于盆地中心的深水沉积物,某些“不整合面”可能是非常细微尺度的,甚至是“整合的”,很难通过传统的层序地层学研究方法进行识别。基于米兰科维奇理论的旋回地层学,以地球轨道周期为“成因”控制,对细粒沉积地层进行定量分析,可以利用天文周期的时间属性约束不同级次的海平面变化旋回,从而进行高频率层序地层划分。
-
Li et al.[56]开发的沉积噪音模型（DYNOT和ρ1）是一种基于轨道调谐后的动态沉积噪音模型,也被证明是一种有效的海平面变化重建方法[56,72⁃73]。沉积岩中不同气候替代指标保存的信息包括“信号”（即轨道周期作用下的产物和“噪声”和无轨道控制作用的影响因素）,如定年误差、差异压实、不稳定沉积和成岩作用等,在固定位置与水深有关的噪声的变化可能与海平面的变化有关。当海平面相对较高时,与水深相关的噪声相对较低,反之亦然。选取DYNOT模型中功率谱中非轨道方差的比值,用于度量气候和海平面替代指标中的噪声；如果与替代指标相关的噪声比较小,噪声的方差可以作为海平面相对变化的指标,在滑动时间窗口中计算出非轨道信号方差和总方差的比值,即DYNOT值；当相对海平面较高的时候,DYNOT较低,反之亦然[56]。DYNOT模型之外的补充模型——独立的lag-1自相关系数模型,即ρ1模型,作为相对海平面变化的第二个独立噪声指标进行检验,噪声的增加会导致ρ1值的降低,反之亦然[56]。基于沉积噪音模型的假设和计算,本文通过对K元素时间序列进行沉积噪声曲线重建,获得两条能够替代同沉积时期相对海平面变化的曲线：DYNOT和ρ1曲线（图9f,g）。金页1井筇竹寺组沉积时DYNOT和ρ1曲线呈现相似的模拟结果,整个沉积时期有9处（图9蓝色条形填充）显著增强的噪声（低信噪比）均出现在eCOCO图谱中相关系数的低值区,即相对浅水、不稳定的沉积环境导致了噪声的增加,亦从另一个角度反映了相对海平面处于下降阶段。
-
三级层序海平面波动的驱动机制长期以来备受争议,但越来越多的地质证据表明天文轨道周期通过调制气候波动来控制海平面的变化[75⁃78],长斜率周期（1.2 Myr）和长偏心率周期（405 kyr）控制着三级层序海平面波动及四级层序海平面波动。对旋回地层学与经典层序地层学的研究,在受海平面变化控制的成因和机理上有相似之处,天文旋回的时间内涵可以标定地质年代,可以作为层序划分的时间标尺。
前人针对页岩的层序地层研究,主要运用T-R旋回理论基础作为层序发育的主控因素,即海退旋回（R）和海进旋回（T）[79⁃80]。本次研究利用沉积噪音模型的计算,根据DYNOT和ρ1模型所指示的相对海平面变化特征,将筇竹寺组划分为4个三级层序和4个T-R旋回（图9j）。值得指出的是,本文通过对K元素深度域测井序列进行滑动窗口频谱分析（图6）,及eCOCO进行最优沉积速率估算的同时（图7）,通过找出的功率图谱上出现变化的深度值,将整个数据区间划分为四个段分别进行滤波,最后调谐得出的时间域曲线在这四段内具有明显不同的沉积速率,故此时的4个三级层序与前文所划分的4段地层是一一对应的。其中,4个层序界面均发育于DYNOT模型中值曲线的高值区,及对应ρ1模型中值曲线的低值区,代表相对海平面下降到最低点,发育层序边界,沉积速率发生改变；4个最大海泛面位于海进旋回T和海退旋回R之间,和DYNOT模型中值曲线的低值区及ρ1模型中值曲线的高值区相对应,代表了相对海平面升高至最大值。结合“浮动”天文年代标尺,SQ1三级层序对应地质年代处于526.86~523.68 Ma,SQ2对应523.68~521.65 Ma,SQ3对应521.65~518.68 Ma,SQ4对应518.68~515.05 Ma,得出每个三级层序的沉积时限分别为3.18,2.03,2.97和3.63 Ma,分别包含8个,5个,7个和9个长偏心率（405 kyr）周期。405 kyr的长偏心率旋回的变化常常与四级层序存在成因联系[81],2.4 Myr周期通过地球和火星的天文共振而形成,也是三级海平面变化旋回的共有周期。图9f显示的DYNOT中值曲线的2.4 Myr滤波曲线呈现出4.5个2.4 Myr变化周期,同时证明相对海平面变化受到2.4 Myr的三级海平面变化旋回的调控作用。
-
寒武纪发生了多次碳同位素漂移事件,碳同位素异常事件在寒武系等的划分和对比中有着广泛的应用[26]。根据所建立的天文年代标尺,将全球无机碳同位素曲线[51]与相对海平面变化曲线[82]进行对比（图9i）。发现ZHUCE正异常对应着相对海平面的高值区,SHICE负异常对应着相对海平面的低值区,SHICE之上的两次明显正异常均对应着相对海平面的极高值。推测可能是因为随着海平面下降,大量有机质被氧化,向海水中释放碳同位素偏轻的无机碳,造成海水中的δ13C变轻[83]。Cremonese et al.[51]在对肖滩剖面的有机碳同位素（δ13Corg）研究中识别了SHICE碳同位素漂移事件（图2）,δ13Corg组成的影响因素众多,不同剖面相同层位之间δ13Corg曲线变化巨大,时常出现δ13Ccarb组成与δ13Corg组成不同步变化的情况,在肖滩剖面的δ13Corg曲线中未能识别其他碳同位素异常事件。
-
（1） 对四川盆地金页1井筇竹寺组地层进行了测井曲线的响应特征分析,选取K元素序列为相对最合适的替代指标,并按测井曲线的不同特征筇竹寺组地层分为四段进行频谱分析和滤波分析。发现目标层段记录的天文轨道控制作用下的沉积周期分别为405 kyr,131 kyr,95 kyr,36 kyr,28 kyr；根据长偏心率周期的天文校准建立了连续约11.8 Myr的 “浮动”天文年代标尺。
（2） 通过轨道调谐后的沉积噪音模型,恢复了沉积时期的相对海平面变化。根据DYNOT和ρ1曲线,将筇竹寺组划分为4个三级层序,识别了对应的4个层序界面和4个最大海泛面,证明了405 kyr的长偏心率旋回与四级层序存在成因上的联系。
（3） 通过相对海平面变化曲线与全球性碳同位素漂移曲线的对比,发现海平面下降往往对应着碳同位素的负异常,推测是由于海水下降导致的有机质氧化,从而轻碳富集。
Cyclostratigraphy Research on Well-logging of the Lower Cambrian Qiongzhusi Formation in Southwestern Sichuan Basin
-
摘要: 目的 四川盆地下寒武统的泥页岩不仅包含重要的环境生命演化信息,同时也是我国当前页岩气勘探开发的重点层位,对其开展精细地层年代学研究具有生产及科研的双重价值。 方法 为了更准确地探明四川盆地下寒武统筇竹寺组黑色页岩层序沉积特征,对川西南地区金页1井筇竹寺组页岩地层开展了旋回地层学研究。 结果与结论 （1）通过分析不同测井曲线测井响应特征,认为钾元素测井序列是最适合古环境、古气候的替代指标。（2）将钾元素深度序列进行频谱分析和滑动窗口频谱分析,识别出偏心率和斜率周期（405 kyr,131 kyr,95 kyr,36.4 kyr,28.2 kyr）；并结合相关系数分析,以不同沉积速率为限,将整段曲线划分为四段,分别进行405 kyr长偏心率周期的带通滤波；以周缘地区的年代地层年龄526.86±0.16 Ma为年龄控制点,建立了金页1井筇竹寺组的“浮动”天文年代标尺。（3）筇竹寺组共记录了29个长偏心率周期,四段地层的最优沉积速率分别为2.5,2.9,3.4和4.8 cm/kyr。（4）以405 kyr的轨道调谐为基础的沉积噪音模型DYNOT及ρ1曲线重建了筇竹寺组沉积时期相对海平面变化曲线；以该曲线为理论参照,在筇竹寺组中划分了4个三级层序,识别出4个最大海泛面及与之对应的T-R旋回；4个三级层序的沉积时限从下到上为3.18,2.03,2.97和3.63 Ma,分别包含了8个,5个,7个和9个长偏心率周期。（5）海平面变化曲线的低值对应了全球碳同位素负偏,推测是由于海平面的下降引起的有机质发生氧化,并向海水释放碳同位素偏轻的无机碳,形成碳同位素负偏。Abstract: Objective The Lower Cambrian in the Sichuan Basin is a key area for understanding the evolution of life and the paleoenvironment,and is a most significant area for future oil and gas exploration. Study of the precise chronostratigraphy has both production and scientific value. Methods To define the sedimentary sequence in the region,cyclostratigraphical studies were conducted for well Jinye-1. [ Results and Conclusions ] （1）Well-logging records were taken to investigate the effects of content variation in strata components. （2）The K logging series was chosen to be the best proxy in the Lower Cambrian Qiongzhusi Formation in the southwestern Sichuan Basin. The K logging series was divided into four segments,and spectra for each segment were analyzed in detail. Evolutional spectrum analysis was conducted for the whole Qiongzhusi Formation. Orbital cycles of 405,131,95,6 and 28 ka were identified in the sedimentary records. A floating astronomical time scale was constructed by applying Gaussian filtering of the 405 ka cycle together with the control age of 526.86 ± 0.16 Ma from a related region. （3）The time scale contained 29 long-eccentricity cycles,and the optimal sedimentation rate of each logging series segment was found to be 2.5,2.9,3.4 and 4.8 cm/ka. （4）The DYNOT sedimentary noise model and ρ1 were built on the 405 ka tuning time series to obtain the relative sea-level change plots in the time domain. Analysis of relative sea-level change recognized four third-order sequences and four maximum flooding surfaces,and the transgressive-regressive （T-R） cycles. The deposition times of each sequence from early to late are 3.18,2.03,2.97 and 3.63 Ma respectively,with each containing 8,5,7 and 9 long eccentricity cycles. （5）The relative sea-level decrease corresponds to the global carbon isotope negative excursion； it is presumed that the decreasing sea-level led to oxidation of exposed organic matter,releasing the lighter carbon isotopes and resulting in the negative carbon isotope excursion.
-
图 6 金页1井筇竹寺组（3 184~3 566 m）K元素曲线频谱分析图
(a) evolutionary spectral analysis of 3 184⁃3 566 m; (b) spectral analysis of 3 184⁃3 303 m; (c) spectral analysis of 3 303⁃3 384 m; (d) spectral analysis of 3 384⁃3 496 m; (e) spectral analysis of 3 496⁃3 566 m
图 8 金页1井筇竹寺组K元素曲线滤波、调谐、时间域滑动窗口频谱分析和浮动天文年代标尺
(a) lithological profile; (b) de⁃trended K series in depth domain; (c) 405⁃kyr filtered K series in depth domain (bandpasses of the four sections: 0.10±0.02, 0.053±0.009, 0.085±0.015 and 0.075±0.001 cycles/m respectively upwards from bottom); (d) sedimentation rate in depth domain; (e) 405⁃kyr tuned K series in time domain; (f) 405⁃kyr filtering curve in time domain; (g) evolutional spectra analysis of tuned K series; (h) floating astronomical time scale
图 9 金页1井筇竹寺组沉积噪音模型与相对海平面变化解释[51,74]
(a) lithological profile; (b) de⁃trended K series in depth domain; (c) ρ1 curve in depth domain; (d) eCOCO and sedimentation rate in depth domain; (e) 405 kyr tuned K series (black line) and 405 kyr filtering curve (red line) in time domain; (f) DYNOT model in time domain; (g) ρ1 model in time domain; (h) global sea⁃level change in GTS2020; (i)
-
[1] 黄金亮,邹才能,李建忠,等. 川南下寒武统筇竹寺组页岩气形成条件及资源潜力[J]. 石油勘探与开发,2012,39（1）：69-75. Huang Jinliang,Zou Caineng,Li Jianzhong,et al. Shale gas generation and potential of the Lower Cambrian Qiongzhusi Formation in southern Sichuan Basin,China[J]. Petroleum Exploration and Development,2012,39（1）： 69-75. [2] 王鹏威,刘忠宝,金之钧,等. 川西南地区下寒武统筇竹寺组页岩气纵向差异富集主控因素[J]. 地球科学,2019,44（11）：3628-3638. Wang Pengwei,Liu Zhongbao,Jin Zhijun,et al. Main control factors of shale gas differential vertical enrichment in Lower Cambrian Qiongzhusi Formation,southwest Sichuan Basin,China[J]. Earth Science,2019,44（11）： 3628-3638. [3] 金之钧,胡宗全,高波,等. 川东南地区五峰组—龙马溪组页岩气富集与高产控制因素[J]. 地学前缘,2016,23（1）：1-10. Jin Zhijun,Hu Zongquan,Gao Bo,et al. Controlling factors on the enrichment and high productivity of shale gas in the Wufeng-Longmaxi Formations,southeastern Sichuan Basin[J]. Earth Science Frontiers,2016,23（1）： 1-10. [4] 刘忠宝,高波,张钰莹,等. 上扬子地区下寒武统页岩沉积相类型及分布特征[J]. 石油勘探与开发,2017,44（1）：21-31. Liu Zhongbao,Gao Bo,Zhang Yuying,et al. Types and distribution of the shale sedimentary facies of the Lower Cambrian in Upper Yangtze area,South China[J]. Petroleum Exploration and Development,2017,44（1）： 21-31. [5] 刘犟,张克银. 井研—犍为地区麦地坪组—筇竹寺组沉积相特征[J]. 科学技术与工程,2018,18（2）：20-25. Liu Jiang,Zhang Keyin. Sedimentary facies characteristics of Maidiping-Qiongzhusi Formation in Jingyan-Qianwei area[J]. Science Technology and Engineering,2018,18（2）： 20-25. [6] 刘犟,张克银. 四川盆地井研地区麦地坪组—筇竹寺组高分辨率层序地层特征[J]. 成都理工大学学报,2018,45（5）：585-593. Liu Jiang,Zhang Keyin. Characteristics of high-resolution stratigraphic sequence of Maidiping Formation and Qiongzhusi Formation in Jingyan area,Sichuan,China[J]. Journal of Chengdu University of Technology （Science & Technology Edition）,2018,45（5）： 585-593. [7] 李延钧,赵圣贤,黄勇斌,等. 四川盆地南部下寒武统筇竹寺组页岩沉积微相研究[J]. 地质学报,2013,87（8）：1136-1148. Li Yanjun,Zhao Shengxian,Huang Yongbin,et al. The sedimentary micro-facies study of the Lower Cambrian Qiongzhusi Formation in southern Sichuan Basin[J]. Acta Geologica Sinica,2013,87（8）： 1136-1148. [8] 范海经,邓虎成,伏美燕,等. 四川盆地下寒武统筇竹寺组沉积特征及其对构造的响应[J]. 沉积学报,2021,39（4）：1004-1019. Fan Haijing,Deng Hucheng,Fu Meiyan,et al. Sedimentary characteristics of the Lower Cambrian Qiongzhusi Formation in the Sichuan Basin and its response to construction[J]. Acta Sedimentologica Sinica,2021,39（4）： 1004-1019. [9] 刘忠宝,杜伟,高波,等. 层序格架中富有机质页岩发育模式及差异分布：以上扬子下寒武统为例[J]. 吉林大学学报（地球科学版）,2018,48（1）：1-14. Liu Zhongbao,Du Wei,Gao Bo,et al. Sedimentary model and distribution of organic-rich shale in the sequence stratigraphic framework： A case study of Lower Cambrian in Upper Yangtze region[J]. Journal of Jilin University （Earth Science Edition）,2018,48（1）： 1-14. [10] 张俊明,王海峰,李国祥. 华北与扬子地台晚寒武世末期至早奥陶世早期层序地层对比及海平面升降事件[J]. 地层学杂志,2000,24（增刊1）：359-369. Zhang Junming,Wang Haifeng,Li Guoxiang. Sequence stratigraphic correlation and eustatic events from latest Cambrian to early Early Ordovician between North China and Yangtze platforms[J]. Journal of Stratigraphy,2000,24（Suppl.1）： 359-369. [11] 梅冥相,马永生,邓军,等. 上扬子区下古生界层序地层格架的初步研究[J]. 现代地质,2005,19（4）：551-562. Mei Mingxiang,Ma Yongsheng,Deng Jun,et al. Sequence-stratigraphic framework for the Early Palaeozoic of the Upper-Yangtze region[J]. Geoscience,2005,19（4）： 551-562. [12] Vail P R. The stratigraphic signatures of tectonics,eustacy and sedimentology： An overview[M]//Eisele G,Rieken W,Seilacher A. Cycles and events in stratigraphy. Berlin： Springer,1991： 617-659. [13] Boulila S,Haq B U,Hara N,et al. Potential encoding of coupling between Milankovitch forcing and Earth's interior processes in the Phanerozoic eustatic sea-level record[J]. Earth-Science Reviews,2021,220： 103727. [14] 龚一鸣,杜远生,童金南,等. 旋回地层学：地层学解读时间的第三里程碑[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： Journal of China University of Geosciences,2008,33（4）： 443-457. [15] 梅冥相. 从旋回的有序叠加形式到层序的识别和划分：层序地层学进展之三[J]. 古地理学报,2011,13（1）：37-54. Mei Mingxiang. From vertical stacking pattern of cycles to discerning and division of sequences： The third advance in sequence stratigraphy[J]. Journal of Palaeogeography,2011,13（1）： 37-54. [16] Lourens L J,Sluijs A,Kroon D,et al. Astronomical pacing of Late Palaeocene to Early Eocene global warming events[J]. Nature,2005,435（7045）： 1083-1087. [17] 吴怀春,张世红,冯庆来,等. 旋回地层学理论基础、研究进展和展望[J]. 地球科学：中国地质大学学报,2011,36（3）：409-428. Wu Huaichun,Zhang Shihong,Feng Qinglai,et al. Theoretical basis,research advancement and prospects of cyclostratigraphy[J]. Earth Science： Journal of China University of Geosciences,2011,36（3）： 409-428. [18] Ruhl M,Hesselbo S P,Hinnov L,et al. Astronomical constraints on the duration of the Early Jurassic Pliensbachian Stage and global climatic fluctuations[J]. Earth and Planetary Science Letters,2016,455： 149-165.] [19] Huang C J,Hesselbo S P. Pacing of the Toarcian Oceanic Anoxic Event （Early Jurassic） from astronomical correlation of marine sections[J]. Gondwana Research,2014,25（4）： 1348-1356. [20] 张若琳,金思丁. 渤海湾盆地沾化凹陷罗69井沙三下亚段旋回地层学研究[J]. 中南大学学报（自然科学版）,2021,52（5）：1516-1531. 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. [21] Guo D,Jin Z J. Astronomical time scale of a Middle Eocene lacustrine sequence from the Dongpu Sag,Bohai Bay Basin,eastern China[J]. Journal of Asian Earth Sciences,2021,214： 104747. [22] Jin S D,Deng H C,Zhu X,et al. Orbital control on cyclical organic matter accumulation in Early Silurian Longmaxi Formation shales[J]. Geoscience Frontiers,2020,11（2）： 533-545. [23] 石巨业,金之钧,刘全有,等. 基于米兰科维奇理论的湖相细粒沉积岩高频层序定量划分[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. [24] 陈中强,杨建国. 米兰柯维奇旋回在我国前第四纪地层之保存[J]. 微体古生物学报,1996,13（1）：65-73. Chen Zhongqiang,Yang Jianguo. Records of Milankovitch Cycles in pre-Quaternary deposits[J]. Acta Micropalaeontologica Sinica,1996,13（1）： 65-73. [25] 邵龙义,汪浩,Large D J. 中国西南地区晚二叠世泥炭地净初级生产力及其控制因素[J]. 古地理学报,2011,13（5）：473-480. Shao Longyi,Wang Hao,Large D J. Net primary productivity and its control of the Late Permian peatlands in southwestern China[J]. Journal of Palaeogeography,2011,13（5）： 473-480. [26] Wu H C,Zhang S H,Jiang G Q,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. [27] Ochoa D,Sierro F J,Hilgen F J,et al. Origin and implications of orbital-induced sedimentary cyclicity in Pliocene well-logs of the western Mediterranean[J]. Marine Geology,2018,403： 150-164. [28] Liu Z H,Huang C J,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. [29] 彭诚. 松辽盆地科探井古环境古气候测井分析方法及其关键问题研究[D]. 北京：中国地质大学（北京）,2019. [30] Ma C,Meyers S R,Sageman B B,et al. Theory of chaotic orbital variations confirmed by Cretaceous geological evidence[J]. Nature,2017,542（7642）： 468-470. [31] Kroon D,Williams T,Pirmez C,et al. Coupled Early Pliocene-Middle Miocene bio-cyclostratigraphy of Site 1006 Reveals orbitally induced cyclicity patterns of Great Bahama Bank carbonate production[J]. Proceedings of the Ocean Drilling Program,Scientific Results,2000,166： 155-166. [32] Li M S,Huang C J,Ogg J,et al. Paleoclimate proxies for cyclostratigraphy： Comparative analysis using a Lower Triassic marine section in South China[J]. Earth-Science Reviews,2019,189： 125-146. [33] 武赛军,魏国齐,杨威,等. 四川盆地桐湾运动及其油气地质意义[J]. 天然气地球科学,2016,27（1）：60-70. Wu Saijun,Wei Guoqi,Yang Wei,et al. Tongwan Movement and its geologic significances in Sichuan Basin[J]. Natural Gas Geoscience,2016,27（1）： 60-70. [34] 张满郎,谢增业,李熙喆,等. 四川盆地寒武纪岩相古地理特征[J]. 沉积学报,2010,28（1）：128-139. Zhang Manlang,Xie Zengye,Li Xizhe,et al. Characteristics of lithofacies paleogeography of Cambrian in Sichuan Basin[J]. Acta Sedimentologica Sinica,2010,28（1）： 128-139. [35] 黄博宇. 四川盆地震旦纪：早寒武世岩相古地理与裂陷槽演化[D]. 北京：中国石油大学（北京）,2018. [36] 刘树根,刘殊,孙玮,等. 绵阳—长宁拉张槽北段构造—沉积特征[J]. 成都理工大学学报（自然科学版）,2018,45（1）：1-13. Liu Shugen,Liu Shu,Sun Wei,et al. Tectonic and sedimentary features of the northern Mianyang-Changning intracratonic sag,Sichuan,China[J]. Journal of Chengdu University of Technology （Science & Technology Edition）,2018,45（1）： 1-13. [37] 朱茂炎,杨爱华,袁金良,等. 中国寒武纪综合地层和时间框架[J]. 中国科学（D辑）：地球科学,2019,49（1）：26-65. Zhu Maoyan,Yang Aihua,Yuan Jinliang,et al. Cambrian integrative stratigraphy and timescale of China[J]. Science China （Seri. D）： Earth Sciences,2019,49（1）： 26-65. [38] 朱茂炎,张智新,杨爱华,等. 中国寒武纪岩石地层划分和对比[J]. 地层学杂志,2021,45（3）： 223-249. Zhu Maoyan,Zhang Zhixin,Yang Aihua,et al. Lithostratigraphic subdivision and correlation of the Cambrian in China[J]. Journal of Stratigraphy,45（3）：223-249. [39] 梁霄,李香华,徐剑良,等. 从优质烃源岩到储层：构造—沉积分异格局下的四川盆地中西部下寒武统页岩气勘探前景[J]. 天然气工业,2021,41（5）：30-41. Liang Xiao,Li Xianghua,Xu Jianliang,et al. Exploration prospects of Lower Cambrian shale gas in the central-western Sichuan Basin under the pattern of tectonic-depositional differentiation： From high-quality source rocks to reservoirs[J]. Natural Gas Industry,2021,41（5）： 30-41. [40] 刘满仓,杨威,李其荣,等. 四川盆地蜀南地区寒武系地层划分及对比研究[J]. 天然气地球科学,2008,19（1）：100-106. Liu Mancang,Yang Wei,Li Qirong,et al. Characteristics and stratigraphic classification and correlation of Cambrian on south Sichuan Basin[J]. Natural Gas Geoscience,2008,19（1）： 100-106. [41] 罗超. 上扬子地区下寒武统牛蹄塘组页岩特征研究[D]. 成都：成都理工大学,2014. [42] 薛耀松,周传明. 扬子区早寒武世早期磷质小壳化石的再沉积和地层对比问题[J]. 地层学杂志,2006,30（1）：64-74. Xue Yaosong,Zhou Chuanming. Resedimentation of the Early Cambrian phosphatized small shell fossils and correlation of the Sinian-Cambrian boundary strata in the Yangtze region,southern China[J]. Journal of Stratigraphy,2006,30（1）： 64-74. [43] 四川省地质矿产局. 四川省区域地质志[M]. 北京：地质出版社,1991. Sichuan Bureau of Geology and Mineral Resources. Regional geology of Sichuan province[M]. Beijing： Geological Publishing House,1991. [44] 潘晓强,董晓霞,代乔坤,等. 四川盆地井研—犍为地区寒武系第二统九老洞组古生物化石及其地层划分对比[J]. 微体古生物学报,2021,38（3）：241-256. Pan Xiaoqiang,Dong Xiao-xia,Dai Qiaokun,et al. Fossil assemblages and stratigraphic correlation of the Jiulaodong Formation （Early Cambrian） in the Jingyan-Qianwei area,Sichuan Basin[J]. Acta Micropalaeontologica Sinica,2021,38（3）： 241-256. [45] 王承红. 井研—犍为地区下寒武统筇竹寺组海相页岩孔隙演化[D]. 成都：成都理工大学,2017. [46] Zhu R X,Li X H,Hou X G,et al. SIMS U-Pb zircon age of a tuff layer in the Meishucun section,Yunnan,southwest China： Constraint on the age of the Precambrian-Cambrian boundary[J]. Science in China Series D： Earth Sciences,2009,52（9）： 1385-1392. [47] Sawaki Y,Nishizawa M,Suo T,et al. Internal structures and U-Pb ages of zircons from a tuff layer in the Meishucunian Formation,Yunnan province,South China[J]. Gondwana Research,2008,14（1/2）： 148-158. [48] Jenkins R J E,Cooper J A,Compston W. Age and biostratigraphy of Early Cambrian tuffs from SE Australia and southern China[J]. Journal of the Geological Society,2002,159（6）： 645-658. [49] Compston W,Zhang Z C,Cooper J A,et al. Further SHRIMP geochronology on the Early Cambrian of South China[J]. American Journal of Science,2008,308（4）： 399-420. [50] 杨传. 华南埃迪卡拉系上部：寒武系下部高精度地质年代学[D]. 北京：中国科学院地质与地球物理研究所,2017. [51] Cremonese L,Shields-Zhou G,Struck U,et al. Marine biogeochemical cycling during the Early Cambrian constrained by a nitrogen and organic carbon isotope study of the Xiaotan section,South China[J]. Precambrian Research,2013,225： 148-165. [52] 赵建华,金之钧,林畅松,等. 上扬子地区下寒武统筇竹寺组页岩沉积环境[J]. 石油与天然气地质,2019,40（4）：701-715. Zhao Jianhua,Jin Zhijun,Lin Changsong,et al. Sedimentary environment of the Lower Cambrian Qiongzhusi Formation shale in the Upper Yangtze region[J]. Oil & Gas Geology,2019,40（4）： 701-715. [53] Thomson D J. Spectrum estimation and harmonic analysis[J]. Proceedings of the IEEE,1982,70（9）： 1055-1096. [54] Weedon G. Time-series analysis and cyclostratigraphy[M]. Cambridge： Cambridge University Press,2003： 276. [55] Kodama K P,Hinnov L A. Rock magnetic cyclostratigraphy： New analytical methods in earth and environmental science series[M]. Chichester： Wiley-Blackwell,2014： 1-147. [56] Li M S,Hinnov L A,Huang C J,et al. Sedimentary noise and sea levels linked to land-ocean water exchange and obliquity forcing[J]. Nature Communications,2018,9（1）： 1004. [57] Cleveland W S. Robust locally weighted regression and smoothing scatterplots[J]. Journal of the American Statistical Association,1979,74（368）： 829-836. [58] Li M S,Hinnov L,Kump L. Acycle： Time-series analysis software for paleoclimate research and education[J]. Computers & Geosciences,2019,127： 12-22. [59] 代大经,唐正松,陈鑫堂,等. 铀的地球化学特征及其测井响应在油气勘探中的应用[J]. 天然气工业,1995,15（5）：21-24. Dai Dajing,Tang Zhengsong,Chen Xintang,et al. U geochemical characteristics and the application of U logging response to oil-gas explobation[J]. Natural Gas Industry,1995,15（5）： 21-24. [60] Wignall P B,Myers K J. Interpreting benthic oxygen levels in mudrocks： A new approach[J]. Geology,1988,16（5）： 452-455. [61] Algeo T J,Rowe H. Paleoceanographic applications of trace mental concentration data[J]. Chemical Geology,2012,324-325： 6-18. [62] 刘光泓,张世红,吴怀春. 前寒武纪旋回地层学研究的进展与挑战[J]. 地层学杂志,2020,44（3）：239-249. Liu Guanghong,Zhang Shihong,Wu Huaichun. Progress and challenges in Precambrian cyclostratigraphy research[J]. Journal of Stratigraphy,2020,44（3）： 239-249. [63] Berger A,Loutre M F. Astronomical forcing through geological time[M]//de Boer P L,Smith D G. Orbital forcing and cyclic sequences. Oxford： Blackwell Scientific Publications,1994： 15-24. [64] 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. [65] Waltham D. Milankovitch Period uncertainties and their impact on cyclostratigraphy[J]. Journal of Sedimentary Research,2015,85（8）： 990-998. [66] 纪永朝. 云南会泽—井寒武纪早期旋回地层学研究[D]. 北京：中国地质大学（北京）,2020. [67] 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. [68] 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. [69] Catuneanu O,Abreu V,Bhattacharya J P,et al. Towards the standardization of sequence stratigraphy[J]. Earth-Science Reviews,2009,92（1/2）： 1-33. [70] Miller K G,Kominz M A,Browning J V,et al. The Phanerozoic record of global sea-level change[J]. Science,2005,310（5752）： 1293-1298. [71] Boulila S,Galbrun B,Miller K,et al. On the origin of Cenozoic and Mesozoic “third-order” eustatic sequences[J]. Earth-Science Reviews,2011,109： 94-112. [72] Wang M,Chen H H,Huang C J,et al. Astronomical forcing and sedimentary noise modeling of lake-level changes in the Paleogene Dongpu Depression of North China[J]. Earth and Planetary Science Letters,2020,535： 116116. [73] Jin S D,Liu S B,Li Z,et al. Astrochronology of a Middle Eocene lacustrine sequence and sedimentary noise modeling of lake-level changes in Dongying Depression,Bohai Bay Basin[J]. Palaeogeography,Palaeoclimatology,Palaeoecology,2022,585： 110740. [74] Mudelsee M. Climate time series analysis[M]. New York： Springer,2010： 454. [75] Zhu M Y,Babcock L E,Peng S C. Advances in Cambrian stratigraphy and paleontology： Integrating correlation techniques,paleobiology,taphonomy and paleoenvironmental reconstruction[J]. Palaeoworld,2006,15（3/4）： 217-222. [76] Zhong Y Y,Wu H C,Zhang Y D,et al. Astronomical calibration of the Middle Ordovician of the Yangtze Block,South China[J]. Palaeogeography,Palaeoclimatology,Palaeoecology,2018,505： 86-89. [77] Boulila S,Galbrun B,Miller K G,et al. On the origin of Cenozoic and Mesozoic “third-order” eustatic sequences[J]. Earth-Science Reviews,2011,109（3/4）： 94-112. [78] Crampton J S,Meyers S R,Cooper R A,et al. Pacing of Paleozoic macroevolutionary rates by Milankovitch grand cycles[J]. Proceedings of the National Academy of Sciences of the United States of America,2018,115（22）： 5686-5691. [79] 张喜,张廷山,赵晓明,等. 天文轨道周期及火山活动对中上扬子区晚奥陶世：早志留世有机碳聚集的影响[J]. 石油勘探与开发,2021,48（4）：732-744. Zhang Xi,Zhang Tingshan,Zhao Xiaoming,et al. Effects of astronomical orbital cycle and volcanic activity on organic carbon accumulation during Late Ordovician： Early Silurian in the Upper Yangtze area,South China[J]. Petroleum Exploration and Development,2021,48（4）： 732-744. [80] Abouelresh M O,Slatt R M. Lithofacies and sequence stratigraphy of the Barnett shale in east-central Fort Worth Basin,Texas[J]. AAPG Bulletin,2012,96（1）： 1-22. [81] 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. [82] Ogg J G,Ogg G,Gradstein F M. A concise geologic time scale 2016[M]. Amsterdam： Elsevier,2016： 234-240. [83] Shields G A,Mills B J W,Zhu M Y,et al. Unique Neoproterozoic carbon isotope excursions sustained by coupled evaporite dissolution and pyrite burial[J]. Nature Geoscience,2019,12（10）： 823-827. -
点击查看大图
图(9) / 表 (1)
计量
- 文章访问数: 252
- HTML全文浏览量: 68
- PDF下载量: 86
- 被引次数: 0
