魏小松　嚴德天　龔銀　牛杏　梁萬樂　伏海蛟　劉紫璇　楊向榮　張寶

摘 要 【目的】鄂西—黔南地區(qū)下寒武統(tǒng)富有機質(zhì)頁巖層系是頁巖氣勘探的重要目標(biāo)，查明不同地區(qū)富有機質(zhì)頁巖層段的等時地層關(guān)系及其成因已成為頁巖氣勘探的關(guān)鍵問題?！痉椒ā客ㄟ^對鄂西ND1井及黔南HY1井下寒武統(tǒng)測井自然伽馬數(shù)據(jù)的旋回地層學(xué)分析，識別出下寒武統(tǒng)主要的天文周期，以文獻報道的鋯石年齡作為錨點建立了兩口井的天文年代標(biāo)尺，并基于沉積噪音模型恢復(fù)了兩口井下寒武統(tǒng)時期的相對海平面變化?！窘Y(jié)果】ND1井牛蹄塘組36 m和9 m的沉積旋回及HY1井16.6 m和4.0 m的沉積旋回代表了軌道周期405 kyr和100 kyr偏心率；ND1井牛蹄塘組和HY1井九門沖組持續(xù)時間分別為3.30 Myr和2.64 Myr。給出兩種方案的天文年代標(biāo)尺，方案一以ND1井牛蹄塘組頂作為錨點起點，對應(yīng)的天文年代標(biāo)尺年齡區(qū)間處于535.20±1.70 Ma ～538.50±1.70 Ma；方案二以HY1 井九門沖組頂作為錨點起點，對應(yīng)的天文年代標(biāo)尺年齡區(qū)間處于535.20±1.70 Ma～537.84±1.70 Ma。經(jīng)過對比分析，兩口井海平面變化對應(yīng)于1.20 Myr旋回和兩個三級層序（SQ1和SQ2）。將下寒武統(tǒng)富有機質(zhì)頁巖總有機碳（Total Organic Carbon，TOC）含量和405 kyr旋回濾波及海平面變化曲線對比，結(jié)果顯示TOC高值段對應(yīng)于405 kyr旋回最大值，而ND1井TOC低值段對應(yīng)于兩個405 kyr旋回之間的低值。此外，TOC高值段既對應(yīng)于海平面較高時期，也對應(yīng)于海平面較低時期?！窘Y(jié)論】由于太陽系的混沌行為，無法準確預(yù)測寒武紀早期理論軌道偏心率對富有機質(zhì)頁巖驅(qū)動的機制。假設(shè)偏心率最大值和總有機碳含量最大值對應(yīng)，則更長的時間段內(nèi)強季節(jié)性變化在某種程度上觸發(fā)了黑色頁巖的富集。此外，相對海平面變化和總有機碳含量對比表明海平面變化和有機質(zhì)富集不具有因果關(guān)系。

關(guān)鍵詞 牛蹄塘組；九門沖組；軌道驅(qū)動；偏心率；沉積噪音模型

第一作者簡介 魏小松，男，1992年出生，博士，旋回地層學(xué)，E-mail： weixiaoson_g@163.com

通信作者 嚴德天，男，教授，博士生導(dǎo)師，沉積學(xué)，E-mail： yandetian@cug.edu.cn

中圖分類號 P539 文獻標(biāo)志碼 A

0 引言

新元古代晚期—早古生代早期是華南克拉通地質(zhì)演化的關(guān)鍵階段[1]。這一時期，發(fā)生了包括羅迪尼亞超大陸裂解[2?3]、岡瓦納泛大陸聚合[4]、全球冰川形成與消退[5]、寒武紀生命大爆發(fā)等一系列重大地質(zhì)事件。寒武紀早期，由于寒武紀生命大爆發(fā)為有機質(zhì)富集和烴源巖形成提供了物質(zhì)條件，因此在全球多個地區(qū)廣泛發(fā)育烴源巖[1]。華南上揚子地區(qū)下寒武統(tǒng)普遍沉積了一套黑色頁巖及硅質(zhì)巖互層層系，在空間上呈現(xiàn)北東—南西向狹長型帶狀展布[6]，從云南地區(qū)延伸至浙江一帶。目前，已在川西南筇竹寺組、云南東部麥地坪組、湘鄂西—黔東南地區(qū)牛蹄塘組及黔南黃平九門沖組等地區(qū)獲得工業(yè)頁巖氣流，充分表明上揚子地區(qū)下寒武統(tǒng)具有良好的頁巖氣勘探潛力[7]。

查明不同地區(qū)富有機質(zhì)頁巖層段的等時地層關(guān)系及其成因已成為頁巖氣勘探的關(guān)鍵問題[7]。在頁巖層系中，巖相、凝縮段、沉積旋回及其疊加樣式、典型巖性標(biāo)志層及自然伽馬能譜測井曲線可作為地層劃分的重要依據(jù)[7?12]，從而建立頁巖層系地層格架。然而，海相地層頁巖沉積一般較連續(xù)，在無明顯響應(yīng)標(biāo)志和明顯侵蝕作用面的地層，識別層序界面并建立層序地層格架極其困難。除此之外，華南上揚子不同地區(qū)下寒武統(tǒng)沉積序列時空變化復(fù)雜，巖石地層名稱存在較大差異[13]，不同地區(qū)下寒武統(tǒng)橫向?qū)Ρ冗€未形成統(tǒng)一的認識，因此，基于層序地層及重要標(biāo)志層的等時地層對比存在很大限制。

對火山碎屑層或火山灰層開展放射性鋯石年代學(xué)分析可為地層提供了高精度年齡錨點，但其孤立應(yīng)用具有嚴重的局限性。原因在于，極少研究可以確定巖石所處地層真正的年齡[14]。另外，地層中沉積的火山碎屑層及火山灰層有限，因此難以獲得多個層位的鋯石年齡來對沉積地層進行綜合約束。近年來，放射性同位素年代學(xué)和天文年代學(xué)相結(jié)合的方法在年代地層研究中得到廣泛的應(yīng)用。該方法為地質(zhì)歷史時期的沉積地層提供高分辨率且連續(xù)的年代標(biāo)尺，已成為約束地質(zhì)年代的精確工具[15]。旋回地層學(xué)利用沉積地層的旋回性（沉積韻律）作為地質(zhì)計時器，這些周期性的旋回與地球軌道參數(shù)變化有關(guān)。地球軌道參數(shù)（偏心率、斜率和歲差）的周期性變化引起地球的季節(jié)性變化和固定緯度的太陽輻射量大小，從而影響地球表層沉積過程的變化，并在沉積記錄中體現(xiàn)為周期性的旋回特征[16]。例如，北緯65°冰期旋回是由軌道驅(qū)動引起的太陽輻射量大小控制的[17]。

然而，在50 Ma前的地質(zhì)歷史時期，太陽系的混沌行為導(dǎo)致軌道參數(shù)配置（即天文周期）表現(xiàn)出很大的不確定性[18]，這對斜率和歲差具有較大影響。由于長偏心率（405 kyr）是木星和金星之間相互作用引起的（g2-g5），且木星的質(zhì)量足夠大，因此，長偏心率（405 kyr）在地質(zhì)歷史時期保持相對穩(wěn)定的狀態(tài)[18]。在中生代及古生代地層中，長偏心率旋回已被頻繁用作天文調(diào)諧的目標(biāo)曲線（即調(diào)諧到405 kyr偏心率周期）。例如，Zhao et al.[19]在寒武紀地層識別出明顯的天文軌道周期，并建立了16 Myr長的天文年代標(biāo)尺。周楊等[20]和Zhang et al.[21]對川南下寒武統(tǒng)筇竹寺組富有機質(zhì)頁巖層系開展了旋回地層學(xué)分析及等時地層對比的研究。然而，目前對鄂西—黔南地區(qū)下寒武統(tǒng)富有機質(zhì)頁巖層系的旋回地層學(xué)及等時地層對比分析缺乏研究。因此，查明該地區(qū)頁巖層系等時地層關(guān)系，開展旋回地層學(xué)分析具有重要意義。

本研究擬對鄂西下寒武統(tǒng)牛蹄塘組及黔南下寒武統(tǒng)九門沖組開展旋回地層學(xué)研究。分別以鄂西及黔南地區(qū)兩口工業(yè)鉆井測井伽馬曲線作為替代指標(biāo)，識別出主要的天文周期，建立“浮動”天文年代標(biāo)尺，基于層序地層對比及已經(jīng)報道的絕對年齡，標(biāo)定建立等時地層格架。其次，基于沉積噪音模型方法重建下寒武統(tǒng)目的層段海平面變化并進行對比分析，最終揭示牛蹄塘組和九門沖組富有機質(zhì)黑色頁巖的成因。

1 地質(zhì)背景

華南板塊在寒武紀位于接近赤道的位置（圖1a），由東南的華夏板塊和西北的揚子板塊組成（圖1b）。揚子克拉通形成于華夏地塊和揚子地塊的碰撞作用，在埃迪卡拉—早寒武世形成了被動大陸邊緣盆地[26?27]。由于差異沉降作用，揚子克拉通內(nèi)發(fā)育多個洼地，東南方向為大規(guī)模開放盆地[27]。由西北向東南依次為淺水碳酸鹽臺地、局限潟湖盆地、臺地邊緣和斜坡—深水盆地區(qū)[1，6]。早寒武世臺—盆轉(zhuǎn)換帶海底熱液活動頻發(fā)，大規(guī)模海侵導(dǎo)致水體進一步缺氧，湘鄂西—黔北—黔南地區(qū)發(fā)育富有機質(zhì)炭質(zhì)泥巖、頁巖和硅質(zhì)巖[25，28?29]。隨著海平面下降以及碳酸鹽巖臺地的形成，區(qū)域上開始發(fā)育大規(guī)模碳酸鹽巖沉積[28]。根據(jù)野外露頭剖面及巖性特征，下寒武統(tǒng)頂部發(fā)育的碳酸鹽巖段及震旦系頂部發(fā)育的灰?guī)r段已被當(dāng)作兩套區(qū)域性標(biāo)志層，用于富有機質(zhì)頁巖層系區(qū)域性的對比。

本研究選取的ND1井位于湖北省宜昌市聶佳河鎮(zhèn)，屬于碳酸鹽臺地相（圖1b）。沉積地層自下而上分別為震旦紀燈影組，下寒武統(tǒng)牛蹄塘組、石牌組、天河板組、石龍洞組（圖1c、表1）。下寒武統(tǒng)牛蹄塘組下段沉積為黑色及灰色含碳頁巖，與下伏燈影組白云巖成平行不整合接觸[25]。HY1井位于貴州省貴陽市東部地區(qū)，屬于斜坡相（圖1b）。沉積地層自下而上分別為上震旦統(tǒng)燈影組、下寒武統(tǒng)九門沖組、杷榔—變馬沖組、清虛洞組，中寒武統(tǒng)高臺組，上寒武統(tǒng)婁山關(guān)組（圖1c、表1）。下寒武統(tǒng)九門沖組沉積以深灰—黑灰炭質(zhì)泥巖、硅質(zhì)泥巖、硅質(zhì)巖為主（圖1）。

劉忠寶等[7]對上揚子地區(qū)下寒武統(tǒng)地層進行了分區(qū)和對比，將上揚子地區(qū)下寒武統(tǒng)劃分為8個三級層序，其中HY1井下寒武統(tǒng)劃分為兩個層序（SQ1和SQ2）（圖2），但SQ1和SQ2存在穿時性，不同地區(qū)缺失程度和沉積特征存在明顯差異。在層序地層對比中，HY1井SQ2層序和JY1井筇竹寺組底部（或麥地坪組）地層（SQ2層序）對應(yīng)，但JY1井底部缺失SQ1 層序（圖2）。周楊等[20]對JY1 井及肖灘剖面巖石地層、古生物地層及年代地層進行了對比。JY1井麥地坪組白云巖可與肖灘剖面朱家箐組含磷灰?guī)r對比。

ND1井下寒武統(tǒng)牛蹄塘組頂部發(fā)育一套微晶灰?guī)r和深灰色薄層黏土巖互層，其上覆地層石牌組以泥質(zhì)粉砂巖沉積為主，而HY1井九門沖組頂部同樣發(fā)育一套灰色—深灰色灰?guī)r，與JY1井麥地坪組頂部灰?guī)r對應(yīng)，其上覆地層杷榔—變馬沖組以灰色泥巖、粉砂質(zhì)泥巖及泥質(zhì)粉砂巖沉積為主（圖1，2）。因此，九門沖組頂部灰?guī)r和牛蹄塘組頂部灰?guī)r層可以作為橫向?qū)Ρ葮?biāo)志層。而且，ND1井牛蹄塘組和HY1井九門沖組頁巖和底部燈影組白云巖接觸，因此ND1井牛蹄塘組可與HY1井九門沖組進行層序?qū)Ρ?。依?jù)巖性變化在ND1井牛蹄塘組內(nèi)部可劃分出兩個沉積旋回，對應(yīng)于層序SQ1和SQ2（圖2）。

2 數(shù)據(jù)及方法

2.1 數(shù)據(jù)

研究選取HY1井和ND1井自然伽馬能譜測井（GR）數(shù)據(jù)進行旋回地層學(xué)分析。HY1井九門沖組選取2 310～2 420 m的深度段，其分辨率（采樣間距）為0.125 m。ND1 井牛蹄塘組選取481～750 m 的深度段，其分辨率為0.125 m。此外，HY1 井總有機碳（Total Organic Carbon，TOC）數(shù)據(jù)和ND1井地球化學(xué)數(shù)據(jù)（V、Cr和Ni）用于綜合對比分析。HY1井TOC數(shù)據(jù)來源于文獻[24]。ND1井TOC數(shù)據(jù)在中國地質(zhì)大學(xué)（武漢）教育部重點實驗室完成。步驟如下：（1）將巖心樣品研磨成200目，加入稀鹽酸溶液進行前處理；（2）殘余樣品中稱取1 g粉末樣置于5 mL離心管，加入2 mL濃度為4 mol/L的HCl溶液搖勻混合，靜置12 h直至無氣泡；（3）第一步用鹽酸去除無機碳酸鹽后，多次加入去離子水進行離心分離直至取上清液的pH變?yōu)橹行?。然后將處理的樣品用冷凍干燥機干燥；（4）稱取20 mg樣品用錫紙包好，用鑷子將通過工具壓實折疊包好的樣品放入元素分析儀進行測試。

ND1井地球化學(xué)數(shù)據(jù)（V、Cr和Ni微量元素）測試在武漢上譜科技有限責(zé)任公司利用Agilent 7700eICP-MS分析完成。處理步驟如下：（1）將巖心樣品磨成200目粉末狀置于105 ℃烘干箱12 h；（2）取50 mg粉末樣品放置于高壓密閉的Teflon溶樣器，分別先后依次緩慢加入1.5 mL 的高純氫氟酸和高純硝酸；（3）將其密封后置于烘干箱，擰緊后置于190 ℃恒溫烘干24 h以上；（4）待溶樣冷卻，開蓋后置于140 ℃電熱板蒸干，然后加入1 mL HNO3并再次蒸干；（5）加入1 mL高純氫氟酸、1 mL MQ水和1 mL內(nèi)標(biāo)In（濃度為1×10-6），再次將Teflon溶樣氣的放入鋼套，擰緊后置于190 ℃烘箱中加熱12 h以上；（6）將溶液轉(zhuǎn)入聚乙烯料瓶中，并用2% HNO3稀釋至100 g以備ICPMS測試。

自然伽馬測井?dāng)?shù)據(jù)作為一種理想的古氣候替代指標(biāo)，已被廣泛用于旋回地層學(xué)的研究[35?36]。自然伽馬能譜測量巖石中放射性物質(zhì)（Th、U、K）在衰變過程中放射出的自然伽馬輻射強度。放射性元素易富集在黏土礦物和有機質(zhì)中。泥頁巖中黏土礦物和有機質(zhì)含量高，砂礫巖及碳酸鹽巖中黏土礦物含量較少，因此，泥頁巖一般具有較高的GR值[37?38]。GR值變化反映了巖性的變化，巖性變化則反映了風(fēng)化強度、降雨及河流輸入能力[37?38]，這些均與古氣候變化關(guān)系密切。

2.2 時頻分析方法

本研究采用多個旋回地層學(xué)方法，包括多窗譜估計（multitaper method，MTM）頻譜分析、快速傅里葉變換（Fast Fourier Transform）滑動窗口頻譜分析、高斯帶通濾波、相關(guān)系數(shù)（correlation coefficient，COCO）和演化的相關(guān)系數(shù)（evolutionary correlation coefficient，eCOCO）分析。在頻譜分析之前，首先對等間距的測井曲線進行去趨勢化（detrended）處理，以去除長期趨勢[39]。對去趨勢后的測井GR曲線開展MTM頻譜分析、快速傅里葉變換及沉積速率估算，以確定潛在的天文周期?；诿商乜_模擬的COCO 和eCOCO分析，可以給出最優(yōu)的沉積速率及其隨深度的變化。

在邊緣海環(huán)境中，氣候替代指標(biāo)記錄的動態(tài)非軌道信號（即“噪音”）影響了氣候和相對海平面變化，這種沉積噪音可歸因于風(fēng)暴、潮汐、生物擾動等因素[40]。基于此，Li et al.[40]開發(fā)了沉積噪音模型方法（dynamic noise after orbital tuning，DYNOT 和lag-1autocorrelation coefficient，ρ1）分析，用于模擬相對海平面變化。一般情況下，高海平面對應(yīng)于DYNOT低值和ρ1高值，反之亦然。以上所有方法都在Acyclev2.4.1[39]開源軟件上運行。

3 結(jié)果

3.1 深度域頻譜分析

HY1井頻譜分析顯示超過99%置信度的峰值有16.6 m，4 m，2.32 m，1.78 m和1.43 m，且16.6 m和4 m的沉積旋回比值約為4：1（圖3a）?？焖俑道锶~變化圖譜顯示在2 390～2 400 m和2 310～2 365 m深度段16.6 m的波長信號較強，而在2 365～2 390 m深度段4 m的波長信號強（圖3b）。ND1井頻譜分析顯示超過99%置信度的峰值有36 m，9 m，3.3 m，1.9 m，1.5 m和1.25 m，且36 m和9 m的沉積旋回比值為4∶1，表明36 m 和9 m 的沉積旋回可能為405 千年（kyr）和100 kyr偏心率（圖3c）。快速傅里葉變換圖譜顯示36 m的波長在480～750 m間隔信號都很強，而9 m的波長僅在480～550 m間隔具有強信號（圖3d）。

對HY1井開展COCO和eCOCO分析，結(jié)果顯示了三個最佳沉積速率：4.1 cm/kyr，13.5 cm/kyr 和16.5 cm/kyr（圖4）。由于4.1 cm/kyr對應(yīng)的天文周期個數(shù)為7，且eCOCO圖譜顯示整個層段約4.1 cm/kyr的強信號，因此，HY1井最佳沉積速率為4.1 cm/kyr。對應(yīng)地，16.6 m，4 m，2.32 m，1.78 m和1.43 m的沉積旋回的持續(xù)時間大致為405 kyr，98 kyr，57 kyr，43 kyr和35 kyr。同樣地，對ND1井開展COCO和eCOCO分析，結(jié)果顯示最佳沉積速率為9 cm/kyr和12 cm/kyr（圖5）。由于eCOCO圖譜中顯示的沉積速率大部分均值為9 cm/kyr，因此，ND1井36 m，9 m，3.3 m，1.9 m和1.5 m 的沉積旋回持續(xù)時間大致為405 kyr，100 kyr，36 kyr，21 kyr和16.6 kyr。

由于太陽系行星之間的混沌行為，寒武紀早期的天文軌道周期斜率和歲差更接近，部分周期可能出現(xiàn)疊置難以識別[20]，而405 kyr長偏心率周期在地質(zhì)歷史時期保持穩(wěn)定。因此，HY1井16.6 m和ND1井36 m的沉積旋回代表了405千年長偏心率周期。按照Waltham[41]的天文周期方案，520 Ma左右的天文周期為405 kyr，131 kyr，95 kyr，36.4 kyr，28.2 kyr，20.4 kyr，19.4 kyr 及16.8 kyr。因此，HY1 井4 m 和1.43 m的沉積旋回分別代表了短偏心率周期和斜率周期；ND1井9 m和3.3 m的沉積旋回代表了短偏心率周期和斜率周期，1.9 m和1.5 m的沉積旋回代表了歲差周期。

3.2 時間域頻譜分析

以HY1井16.6 m和ND1井36 m沉積旋回對應(yīng)405 kyr偏心率周期建立年代模型，基于該年代模型分別將兩口井深度測井序列調(diào)諧至?xí)r間域。對調(diào)諧后的時間域數(shù)據(jù)進行插值和去趨勢處理，然后進行頻譜分析和快速傅里葉變換分析。HY1井時間域頻譜分析顯示了明顯的405 kyr，100 kyr及35 kyr的長偏心率、短偏心率及斜率周期，快速傅里葉變換圖譜顯示了對應(yīng)的米蘭科維奇周期信號（圖3e，f）。ND1井時間域頻譜分析顯示了明顯的405 kyr，100 kyr及36 kyr的長偏心率、短偏心率及斜率周期以及22 kyr和20 kyr的歲差（圖3g，h）。

4 討論

4.1 鄂西—黔南地區(qū)下寒武統(tǒng)富有機質(zhì)頁巖天文年代標(biāo)尺

肖灘剖面年代地層學(xué)研究表明，朱家箐組頂界面年齡約為526.50±1.10 Ma～526.86±0.16 Ma[30?31]。根據(jù)層序地層對比，HY1 井九門沖組頂部灰?guī)r及ND1井牛蹄塘組灰?guī)r頂界面可能為526.86±0.16 Ma，但灰?guī)r底界面年齡尚不清楚。由于肖灘剖面朱家箐組上部第一套硅質(zhì)白云巖年齡為535.20±1.70 Ma～538.20 Ma±1.50 Ma[32?34]，因此推測HY1井和ND1井灰?guī)r底界面年齡（2 300 m 和480 m 處）年齡可能為535.20±1.70 Ma～538.20 Ma±1.50 Ma 或更年輕。因此，本次研究將535.20±1.70 Ma作為天文年代標(biāo)尺的錨點。

旋回地層學(xué)分析表明鄂西—黔南下寒武統(tǒng)沉積旋回受控于軌道周期。通過將HY1井和 ND1井深度域序列調(diào)諧到時間域序列，建立了兩口井“浮動”天文年代標(biāo)尺。研究結(jié)果表明，HY1井九門沖組富有機質(zhì)頁巖持續(xù)時間為2.64 Myr（圖3f、圖6），ND1井下寒武統(tǒng)含碳頁巖持續(xù)時間僅為3.30 Myr（圖3h、圖6）。由于二者持續(xù)時間不一致，因此本研究給出兩種絕對天文年代標(biāo)尺方案。方案一將535.20±1.70 Ma 作為ND1 井牛蹄塘組頂部錨點，方案二以535.20±1.70 Ma作為HY1井九門沖組頂部錨點。由于層序地層對比并不能建立等時地層，但同一區(qū)域的海平面變化是同步的。因此，對比兩口井下寒武統(tǒng)時期的相對海平面變化，可以在時間上建立對應(yīng)關(guān)系。

ND1井位于臺地相，HY1井位于深水斜坡相，二者所處位置反映的海平面變化應(yīng)具有一致性。Li etal.[40]開發(fā)了沉積噪音模型（DYNOT和ρ1）重建大陸邊緣海平面變化，該方法已被證明是有效的海（或湖）平面重建方法，被廣泛用于全球大洋及陸地古湖平面變化的重建中。本研究對調(diào)諧后的自然伽馬序列進行沉積噪音模擬，分別獲取了HY1井和ND1井的DYNOT和ρ1曲線（圖6）。HY1井相對海平面變化反映了兩次升降旋回，這和劉忠寶等[7]劃分的兩個三級層序?qū)?yīng)。此外，兩次海平面升降旋回大致對應(yīng)兩個1.20 Myr長周期旋回，表明了該時期三級海平面層序可能與1.20 Myr斜率周期存在聯(lián)系。相比之下，ND1井總體上以牛蹄塘組下段高海平面和上段低海平面為特征，相對海平面變化反映了一個下降半旋回和兩次上升—下降旋回。根據(jù)其變化的幅度，將兩次上升—下降旋回和HY1井九門沖組海平面變化對應(yīng)，結(jié)果表明4條海平面變化曲線吻合度較高，且HY1井和ND1井405 kyr長偏心率旋回及1.20 Myr長周期旋回均對應(yīng)良好（圖6）。

盡管本研究建立了兩種方案的天文年代標(biāo)尺，且根據(jù)相對海平面變化、偏心率及斜率調(diào)幅旋回及沉積層序?qū)Ρ却笾陆⒘藘煽诰牡葧r關(guān)系，將這兩種方案的誤差控制在0.10 Ma，但由于控制錨點年齡本身具有誤差和限制性，故本研究所建立的天文年代標(biāo)尺需要進一步檢驗和驗證。對比中國寒武紀最新年代地層和時間框架[42]，表明ND1井牛蹄塘組及HY1井九門沖組處于寒武紀幸運階。

4.2 富有機質(zhì)頁巖的驅(qū)動因素

貴州黔南地區(qū)下寒武統(tǒng)九門沖組巖性主要為硅質(zhì)、碳質(zhì)頁巖及泥巖，具有厚度大、有機質(zhì)豐度高的特征，總有機碳含量介于2%～10%，自然伽馬測井值高，其變化趨勢和TOC含量具有相似性（圖7）。研究表明，下寒武統(tǒng)九門沖組為深水陸棚沉積環(huán)境[43]。前期地球化學(xué)分析表明九門沖組黑色頁巖沉積過程受熱液活動的影響[44]，深部物質(zhì)通過深大斷裂上涌進入海水，提供了豐富的物質(zhì)來源。

鄂西地區(qū)下寒武統(tǒng)牛蹄塘組巖性主要為含碳頁巖，具有厚度大、有機質(zhì)豐度高的特征。TOC含量介于2%～20%，且牛蹄塘組下段TOC含量明顯高于上段（圖8）。地球化學(xué)元素比V/Cr和V/（V+Ni）顯示和TOC含量相似的變化，表明高TOC段水體處于缺氧的環(huán)境[28，45]。前期研究已經(jīng)證實，ND1井牛蹄塘組下段頁巖處于缺氧—硫化環(huán)境，上段頁巖下部處于缺氧環(huán)境，上部處于次氧化—缺氧環(huán)境[26]。熱液活動同樣為黑色頁巖的形成提供了物質(zhì)條件，且早期熱液活動明顯強于晚期，這主要與寒武紀早期揚子板塊的快速拉張裂解有關(guān)[28]。上涌的熱液與碎屑物質(zhì)混合后，則形成了大陸邊緣背景下具有洋盆特征的黑色巖系。

生產(chǎn)力升高，大洋內(nèi)部氧氣含量下降、適宜的沉積速率、有機碳埋藏增加及火山活動通常被認為是有機質(zhì)富集的主要因素[14，46?49]，而軌道驅(qū)動則被認為是這些因素的誘因[50?51]。研究表明，地球軌道周期對黑色頁巖的形成具有控制作用[51]。偏心率調(diào)幅的歲差氣候可以影響大陸風(fēng)化、區(qū)域徑流、海洋富營養(yǎng)化以及生產(chǎn)力和氧化還原條件變化[52?55]。其次，歲差可以影響不同緯度季風(fēng)強度、降雨和溫度的差異。很多研究提出了偏心率和歲差氣候?qū)Ω挥袡C質(zhì)頁巖形成的解釋模型。如Wei et al.[45]和Zhang et al.[56]提出在新生代陸相湖盆中，歲差處于夏季的極小值時期更有利于有機質(zhì)的富集和保存，湖泊生產(chǎn)力水平也在此時處于高值。當(dāng)夏季歲差處于極大值時期，湖平面下降導(dǎo)致有機質(zhì)被動處于氧化的環(huán)境，因此不利于有機質(zhì)的保存。

鄂西—黔南地區(qū)處于被動大陸邊緣背景下的淺水—深水陸棚環(huán)境，軌道的405 kyr偏心率旋回對頁巖層系沉積具有明顯的控制作用（圖3，7，8）。研究結(jié)果顯示，頁巖沉積層段有機質(zhì)含量和405 kyr旋回在相位上具有同步幅度的變化，這在ND1井牛蹄塘組下段尤為明顯。ND1井牛蹄塘組下段兩個TOC高值段對應(yīng)于405 kyr旋回最大值，其持續(xù)時間大約為405 kyr（圖8灰色區(qū)域），而兩個TOC含量高值之間的低值段正好對應(yīng)了兩個405 kyr旋回之間的低值（圖8，約720 m）。此外，高TOC值段明顯對應(yīng)的405 kyr旋回幅度變化大，而相對低的TOC值段對應(yīng)的405 kyr旋回幅度變化?。▓D7和圖8灰色區(qū)域）。由于歲差旋回受偏心率的調(diào)制作用影響。因此，推測405 kyr旋回幅度大的時期對應(yīng)的歲差周期幅度波動也較大。然而，由于太陽系混沌行為，無法得到精確的偏心率結(jié)論[18]。在早寒武世，軌道對有機質(zhì)富集的影響只能通過假設(shè)，無法確定偏心率最大值對應(yīng)TOC的最大值或最小值，從而能夠解釋偏心率驅(qū)動的生產(chǎn)力或偏心率控制的保存條件對有機質(zhì)富集的影響[57]。Wichern et al.[58]提出了海相沉積的偏心率假設(shè)模型，當(dāng)偏心率最小值處于2.40 Myr長周期節(jié)點，季節(jié)性反差是可以避免的，這將導(dǎo)致大陸地區(qū)風(fēng)化層建造的形成。在第一次劇烈的歲差最大值期間，這些累計的風(fēng)化層及其所含的養(yǎng)分將被侵蝕，導(dǎo)致水體富營養(yǎng)化。隨后的歲差最小值期間，水文循環(huán)減弱，但營養(yǎng)元素的再循環(huán)也能夠?qū)е氯毖鯒l件下黑色頁巖的形成。本次研究援引這一假設(shè)模型對本文結(jié)果進行解釋（圖9）。由于歲差引起的太陽輻射變化在南、北兩個半球是相反的，因此，兩個半球低緯季風(fēng)氣候的歲差尺度響應(yīng)也具有相反的變化趨勢[59]。早寒武世上揚子地區(qū)處于北半球低緯度，因此在歲差最小值期間，水文循環(huán)和亞洲夏季風(fēng)降雨增強，使得風(fēng)化層的養(yǎng)分被侵蝕，導(dǎo)致水體富營養(yǎng)化。假如總有機碳含量最大值和理論軌道偏心率最大值同步，推測在更長的時間段內(nèi)強季節(jié)性變化在某種程度上觸發(fā)了黑色頁巖的富集[14]，而偏心率最大值也可以導(dǎo)致氣候變暖，更有利于碳的吸收[60]。

從總有機碳含量和相對海平面變化的對比可以看出，海平面下降時期（2 320～2 330 m 和630～660 m）對應(yīng)于較低TOC段（圖7，8），而在TOC高值段，海平面總體較高，這表明富有機質(zhì)頁巖的形成和海平面高低具有較大的聯(lián)系。然而，HY1 井2350～2 370 m沉積間隔，相對下降的海平面對應(yīng)于相對高的總有機碳含量；ND1井660～690 m間隔，高海平面段對應(yīng)于低總有機碳含量。以上結(jié)果表明，海平面變化和有機質(zhì)富集并不是因果關(guān)系，可能存在其他的驅(qū)動因素。由于海平面變化受到1.20 Myr旋回的調(diào)制，而高有機質(zhì)可能受到偏心率最大值的調(diào)制。因此，推測地球軌道參數(shù)變化是驅(qū)動海平面變化和富有機質(zhì)黑色頁巖富集的主要因素。ND1 井高TOC段持續(xù)了約810 kyr，TOC突然升高表明這一時期可能穿過了某個閾值，使得大洋缺氧[14]。前期研究表明深部熱液為有機質(zhì)富集提供了豐富的營養(yǎng)物質(zhì)來源，推測熱液活動可能形成于有機質(zhì)富集之前，而在高總有機碳含量期間死亡的生物得以快速沉降富集。

5 結(jié)論

（1） 旋回地層學(xué)分析表明軌道周期信號記錄于下寒武統(tǒng)沉積地層。兩口井36 m和16.6 m的沉積旋回代表了長偏心率405 kyr，9 m 和 4m 的沉積旋回代表了短偏心率 100 kyr。以長偏心率作為天文調(diào)諧基準，以文獻報道的絕對年齡535.20±1.70 Ma作為HY1井九門沖組和ND1井牛蹄塘組頂界錨點，建立了兩種方案的天文年代標(biāo)尺。ND1 井年齡區(qū)間為535.20±1.70 Ma～538.50±1.70 Ma；HY1井年齡區(qū)間為535.20±1.70 Ma～537.84±1.70 Ma。

（2） 基于沉積噪音模型建立了兩口井下寒武統(tǒng)沉積時期的相對海平面變化。研究結(jié)果表明，兩口井的相對海平面變化具有同步性，對應(yīng)于兩個沉積層序和1.20 Myr的長周期斜率波動旋回。

（3） 下寒武統(tǒng)富有機質(zhì)頁巖TOC高值段對應(yīng)于405 kyr最大值，TOC低值段對應(yīng)于405 kyr最小值，這表明黑色頁巖的富集可能與更長時間段內(nèi)強季節(jié)變化有關(guān)。相比之下，TOC高值段既對應(yīng)相對海平面較高時期，也對應(yīng)較低時期，暗示了相對海平面變化與有機質(zhì)富集不具有因果聯(lián)系。

致謝 審稿專家、編輯部老師提供了建設(shè)性意見，使得論文質(zhì)量顯著提升，在此一并感謝。同時感謝中國石油化工股份有限公司提供的測井?dāng)?shù)據(jù)。

參考文獻（References）

［1］ 祝慶敏，盧龍飛，潘安陽，等. 湘西地區(qū)下寒武統(tǒng)牛蹄塘組頁巖沉積環(huán)境與有機質(zhì)富集[J]. 石油實驗地質(zhì)，2021，43（5）：797-809， 854.［Zhu Qingmin， Lu Longfei， Pan Anyang， et al. Sedimentaryenvironment and organic matter enrichment of the LowerCambrian Niutitang Formation shale， western Hunan province，China[J]. Petroleum Geology & Experiment， 2021， 43（5）： 797-809， 854.］

［2］ Ernst R E， Wingate M T D， Buchan K L， et al. Global record of1600-700 Ma Large Igneous Provinces （LIPs）： Implications forthe reconstruction of the proposed Nuna （Columbia） and Rodiniasupercontinents[J]. Precambrian Research， 2008， 160（1/2）：159-178.

［3］ Ernst R E， Hamilton M A， S?derlund U， et al. Long-lived connectionbetween southern Siberia and northern Laurentia in the Proterozoic[J]. Nature Geoscience， 2016， 9（6）： 464-469.

［4］ Zhao G C， Wang Y J， Huang B C， et al. Geological reconstructionsof the East Asian blocks： From the breakup of Rodinia to theassembly of Pangea[J]. Earth-Science Reviews， 2018， 186：262-286.

［5］ 趙彥彥，鄭永飛. 全球新元古代冰期的記錄和時限[J]. 巖石學(xué)報，2011，27（2）：545-565.［Zhao Yanyan， Zheng Yongfei. Recordand time of Neoproterozoic glaciations on earth[J]. ActaPetrologica Sinica， 2011， 27（2）： 545-565.］

［6］ 周航兵. 湘西北下寒武統(tǒng)底部黑色巖系地球化學(xué)特征及地質(zhì)意義[D]. 南昌：東華理工大學(xué)，2019.［Zhou Hangbing. Geochemicalcharacteristics and geological significance of the LowerCambrian black rock series in northwestern Hunan[D]. Nanchang：East China University of Technology， 2019.］

［7］ 劉忠寶，杜偉，高波，等. 層序格架中富有機質(zhì)頁巖發(fā)育模式及差異分布：以上揚子下寒武統(tǒng)為例[J]. 吉林大學(xué)學(xué)報（地球科學(xué)版），2018，48（1）：1-14.［Liu Zhongbao， Du Wei， Gao Bo， et al.Sedimentary model and distribution of organic-rich shale in the sequencestratigraphic framework： A case study of Lower Cambrianin Upper Yangtze region[J]. Journal of Jilin University （EarthScience Edition）， 2018， 48（1）： 1-14.］

［8］ Hemmesch N T， Harris N B， Mnich C A， et al. A sequencestratigraphicframework for the Upper Devonian Woodford shale，Permian basin， West Texas[J]. AAPG Bulletin， 2014， 98（1）：23-47.

［9］ Brett C E， Baird G C， Bartholomew A J， et al. Sequence stratigraphyand a revised sea-level curve for the Middle Devonian of easternNorth America[J]. Palaeogeography， Palaeoclimatology，Palaeoecology， 2011， 304（1/2）： 21-53.

［10］ Abouelresh M O， Slatt R M. Lithofacies and sequence stratigraphyof the Barnett shale in East-Central Fort Worth Basin， Texas[J]. AAPG Bulletin， 2012， 96（1）： 1-22.

［11］ Halgedahl S L， Jarrard R D， Brett C E， et al. Geophysical andgeological signatures of relative sea level change in the UpperWheeler Formation， Drum Mountains， west-central Utah： A perspectiveinto exceptional preservation of fossils[J]. Palaeogeography，Palaeoclimatology， Palaeoecology， 2009， 277（1/2）：34-56.

［12］ 吳靖，姜在興，吳明昊. 細粒巖層序地層學(xué)研究方法綜述[J].地質(zhì)科技情報，2015，34（5）：16-20.［Wu Jing， Jiang Zaixing，Wu Minghao. Summary of research methods about the sequencestratigraphy of the fine-grained rocks[J]. Geological Science andTechnology Information， 2015， 34（5）： 16-20.］

［13］ 梅冥相，張海，孟曉慶，等. 上揚子區(qū)下寒武統(tǒng)的層序地層劃分和層序地層格架的建立[J]. 中國地質(zhì)，2006，33（6）：1292-1304.［Mei Mingxiang， Zhang Hai， Meng Xiaoqing， et al. Sequencestratigraphic division and framework of the Lower Cambrianin the Upper Yangtze region[J]. Geology in China， 2006，33（6）： 1292-1304.］

［14］ Silva A F， Santos T P， Xavier P L A， et al. Constraining the durationof the southern Gondwana Irati-Whitehill Sea throughcyclostratigraphy and its relation with deep-time astronomical solutions[J]. Palaeogeography， Palaeoclimatology， Palaeoecology，2023， 629： 111791.

［15］ Zeebe R E， Lourens L J. Solar System chaos and the Paleocene–Eocene boundary age constrained by geology and astronomy[J].Science， 2019， 365（6456）： 926-929.

［16］ Hinnov L A. Cyclostratigraphy and its revolutionizing applicationsin the earth and planetary sciences[J]. GSA Bulletin， 2013，125（11/12）： 1703-1734.

［17］ Milankovitch M. Kanon der erdbestrahlungen und seine anwendungauf das eiszeitenproblem[M]. Belgrade： Royal SerbianAcademy， 1941.

［18］ Laskar J， Robutel P， Joutel F， et al. A long-term numerical solutionfor the insolation quantities of the Earth[J]. Astronomy &Astrophysics， 2004， 428（1）： 261-285.

［19］ Zhao Z F， Thibault N R， Dahl T W， et al. Synchronizing rockclocks in the Late Cambrian[J]. Nature Communications， 2022，13（1）： 1990.

［20］ 周楊，金思丁，劉巖，等. 川西南下寒武統(tǒng)筇竹寺組頁巖旋回地層學(xué)研究[J]. 沉積學(xué)報，2024，42（1）：142-157.［Zhou Yang，Jin Siding， Liu Yan， et al. Cyclostratigraphy research on wellloggingof the Lower Cambrian Qiongzhusi Formation in southwesternSichuan Basin[J]. Acta Sedimentological Sinica， 2024，42（1）： 142-157.］

［21］ Zhang T， Li Y F， Fan T L， et al. Orbitally-paced climate changein the Early Cambrian and its implications for the history of theSolar System[J]. Earth and Planetary Science Letters， 2022，583： 117420.

［22］ Scotese C R. Atlas of Cambrian and Early Ordovician PaleogeographicMaps （Mollweide Projection）， maps 81-88， volumes5. The Early Paleozoic， PALEOMAP Atlas for ArcGIS[M].PALEOMAP Project， Evanston， IL， 2014， doi： 10. 13140/2. 1.4756. 7369.

［23］ Ye Y T， Wang H J， Wang X M， et al. Tracking the evolution ofseawater Mo isotopes through the Ediacaran–Cambrian transition[J]. Precambrian Research， 2020， 350： 105929.

［24］ Liu Z B， Gao B， Hu Z Q， et al. Pore characteristics and formationmechanism of high-maturity organic-rich shale in LowerCambrian Jiumenchong Formation， southern Guizhou[J].Petroleum Research， 2018， 3（1）： 57-65.

［25］ 牛杏. 湘鄂西下寒武統(tǒng)牛蹄塘組石英成因類型及其對頁巖儲層物性的影響[D]. 武漢：中國地質(zhì)大學(xué)（武漢），2019.［NiuXing. Genetic types of quartz and its influence on shale reservoirphysical properties of the Lower Cambrian Niutitang Formationin western Hunan and Hubei[D]. Wuhan： China Universityof Geosciences （Wuhan）， 2019.］

［26］ Shu L S， Faure M， Yu J H， et al. Geochronological and geochemicalfeatures of the Cathaysia block （South China）： Newevidence for the Neoproterozoic breakup of Rodinia[J]. Precam‐brian Research， 2011， 187（3/4）： 263-276.

［27］ 付勇，周文喜，王華建，等. 黔北下寒武統(tǒng)黑色巖系的沉積環(huán)境與地球化學(xué)響應(yīng)[J]. 地質(zhì)學(xué)報，2021，95（2）：536-548.［FuYong， Zhou Wenxi， Wang Huajian， et al. The relationship betweenenvironment and geochemical characteristics of blackrock series of Lower Cambrian in northern Guizhou[J]. ActaGeologica Sinica， 2021， 95（2）： 536-548.］

［28］ 梁小聰，牛杏，胡明毅，等. 湘鄂西下寒武統(tǒng)牛蹄塘組黑色頁巖發(fā)育特征及沉積環(huán)境[J]. 巖性油氣藏，2023，35（4）：102-114.［Liang Xiaocong， Niu Xing， Hu Mingyi， et al. Developmentcharacteristics and sedimentary environment of black shaleof Lower Cambrian Niutitang Formation in western Hunan andHubei[J]. Lithologic Reservoirs， 2023， 35（4）： 102-114.］

［29］ 李海，陳孝紅，彭中勤，等. 湘鄂西地區(qū)下寒武統(tǒng)牛蹄塘組頁巖氣儲層特征[J]. 地質(zhì)學(xué)報，2022，96（4）：1421-1433.［Li Hai，Chen Xiaohong， Peng Zhongqin， et al. Shale gas reservoircharacteristics of the Lower Cambrian Niutitang Formation inthe western Hunan and Hubei areas[J]. Acta Geologica Sinica，2022， 96（4）： 1421-1433.］

［30］ 楊傳. 華南埃迪卡拉系上部：寒武系下部高精度地質(zhì)年代學(xué)[D]. 北京：中國科學(xué)院地質(zhì)與地球物理研究所，2017.［YangChuan. High precision geochronology of the Upper Ediacaran-Lower Cambrian in South China[D]. Beijing： Institute of Geologyand Geophysics， 2017.］

［31］ Compston W， Zhang Z C， Cooper J A， et al. Further SHRIMPgeochronology on the Early Cambrian of South China[J]. AmericanJournal of Science， 2008， 308（4）： 399-420.

［32］ Zhu R X， Li X H， Hou X G， et al. SIMS U-Pb zircon age of atuff 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.

［33］ Sawaki Y， Nishizawa M， Suo T， et al. Internal structures and UPbages of zircons from a tuff layer in the Meishucunian Formation，Yunnan province， South China[J]. Gondwana Research，2008， 14（1/2）： 148-158.

［34］ Jenkins R J F， Cooper J A， Compston W. Age and biostratigraphyof Early Cambrian tuffs from SE Australia and southern China[J]. Journal of the Geological Society， 2002， 159（6）： 645-658.

［35］ 吳懷春，張世紅，黃清華. 中國東北松遼盆地晚白堊世青山口組浮動天文年代標(biāo)尺的建立[J]. 地學(xué)前緣，2008，15（4）：159-169.［Wu Huaichun， Zhang Shihong， Huang Qinghua. Establishmentof floating astronomical time scale for the terrestrial LateCretaceous Qingshankou Formation in the Songliao Basin ofnortheast China[J]. Earth Science Frontiers， 2008， 15（4）：159-169.］

［36］ Huang C J， Hinnov L. Astronomically forced climate evolutionin a saline lake record of the Middle Eocene to Oligocene， JianghanBasin， China[J]. Earth and Planetary Science Letters， 2019，528： 115846.

［37］ Schnyder J， Ruffell A， Deconinck J F， et al. Conjunctive use ofspectral gamma-ray logs and clay mineralogy in defining LateJurassic-Early Cretaceous palaeoclimate change （Dorset， U. K. ）[J]. Palaeogeography， Palaeoclimatology， Palaeoecology， 2006，229（4）： 303-320.

［38］ Li M S， Huang C J， Hinnov L， et al. Obliquity-forced climateduring the Early Triassic hothouse in China[J]. Geology， 2016，44（8）： 623-626.

［39］ Li M S， Hinnov L， Kump L. Acycle： Time-series analysis softwarefor paleoclimate research and education[J]. Computers &Geosciences， 2019， 127： 12-22.

［40］ Li M S， Hinnov L A， Huang C J， et al. Sedimentary noise andsea levels linked to land-ocean water exchange and obliquityforcing[J]. Nature Communications， 2018， 9（1）： 1004.

［41］ Waltham D. Milankovitch period uncertainties and their impacton cyclostratigraphy[J]. Journal of Sedimentary Research， 2015，85（8）： 990-998.

［42］ 朱茂炎，楊愛華，袁金良，等. 中國寒武紀綜合地層和時間框架[J]. 中國科學(xué)（D 輯）：地球科學(xué)，2019，49（1）：26-65.［ZhuMaoyan， Yang Aihua， Yuan Jinliang， et al. Cambrian integrativestratigraphy and timescale of China[J]. Science China （Seri. D）：Earth Sciences， 2019， 62（1）： 26-65.］

［43］ 劉忠寶，高波，胡宗全，等. 高演化富有機質(zhì)頁巖儲層特征及孔隙形成演化：以黔南地區(qū)下寒武統(tǒng)九門沖組為例[J]. 石油學(xué)報，2017，38（12）：1381-1389. ［Liu Zhongbao， Gao Bo， HuZongquan， et al. Reservoir characteristics and pores formationand evolution of high maturated organic rich shale： A case studyof Lower Cambrian Jiumenchong Formation， southern Guizhouarea[J]. Acta Petrolei Sinica， 2017， 38（12）： 1381-1389.］

［44］ 侯東壯，吳湘濱，鄧鑫楠. 貴州銅仁地區(qū)九門沖組黑色頁巖地球化學(xué)特征及成巖環(huán)境研究[J]. 地質(zhì)與勘探，2019，55（3）：779-788.［Hou Dongzhuang， Wu Xiangbin， Deng Xinnan. Geochemicalcharacteristics and diagenetic setting of the JiumenchongFormation black shale in the Tongren area of Guizhouprovince[J]. Geology and Exploration， 2019， 55（3）： 779-788.］

［45］ Wei X S， Deng Y， Yan D T， et al. Organic matter enrichment inAsias palaeolake controlled by the early and Middle Eoceneglobal warming and astronomically driven precessional climate[J]. Marine and Petroleum Geology， 2023， 154： 106342.

［46］ Barclay R S， McElwain J C， Sageman B B. Carbon sequestrationactivated by a volcanic CO2 pulse during Ocean AnoxicEvent 2[J]. Nature Geoscience， 2010， 3（3）： 205-208.

［47］ Kuypers M M M， Pancost R D， Nijenhuis I A， et al. Enhancedproductivity led to increased organic carbon burial in the euxinicNorth Atlantic Basin during the Late Cenomanian oceanic anoxicevent[J]. Paleoceanography， 2002， 17（4）： 1051.

［48］ Schoepfer S D， Shen J， Wei H Y， et al. Total organic carbon， organicphosphorus， and biogenic barium fluxes as proxies forpaleomarine productivity[J]. Earth-Science Reviews， 2015， 149：23-52.

［49］ Tyson R V. The genesis and palynofacies characteristics of marinepetroleum source rocks[J]. Geological Society， London，Special Publications， 1987， 26（1）： 47-67.

［50］ Batenburg S J， de Vleeschouwer D， Sprovieri M， et al. Orbitalcontrol on the timing of oceanic anoxia in the Late Cretaceous[J]. Climate of the Past， 2016， 12（10）： 1995-2009.

［51］ Meyers S R， Sageman B B， Arthur M A. Obliquity forcing of organicmatter accumulation during Oceanic Anoxic Event 2[J].Paleoceanography， 2012， 27（3）： PA3212.

［52］ Rachold V， Brumsack H J. Inorganic geochemistry of Albiansediments from the Lower Saxony Basin NW Germany：Palaeoenvironmental constraints and orbital cycles[J]. Palaeogeography，Palaeoclimatology， Palaeoecology， 2001， 174（1/2/3）：121-143.

［53］ Gambacorta G， Menichetti E， Trincianti E， et al. Orbital controlon cyclical primary productivity and benthic anoxia： Astronomicaltuning of the Telychian Stage （Early Silurian） [J]. Palaeogeography，Palaeoclimatology， Palaeoecology， 2018， 495：152-162.

［54］ Lu M， Lu Y H， Ikejiri T， et al. Periodic oceanic euxinia and terrestrialfluxes linked to astronomical forcing during the Late DevonianFrasnian–Famennian mass extinction[J]. Earth and PlanetaryScience Letters， 2021， 562： 116839.

［55］ Ma W T， Tian J， Li Q Y， et al. Simulation of long eccentricity（400-kyr） cycle in ocean carbon reservoir during Mioceneclimate optimum： Weathering and nutrient response to orbitalchange[J]. Geophysical Research Letters， 2011， 38（10）：L10701.

［56］ Zhang J G， Jiang Z X， Liang C， et al. Astronomical forcing ofmeter-scale organic-rich mudstone – limestone cyclicity in theEocene Dongying Sag， China： Implications for shale reservoirexploration[J]. AAPG Bulletin， 2022， 106（8）： 1557-1579.

［57］ Jin S D， Deng H C， Zhu X， et al. Orbital control on cyclical organicmatter accumulation in Early Silurian Longmaxi Formationshales[J]. Geoscience Frontiers， 2020， 11（2）： 533-545.

［58］ Wichern N M A， Bialik O M， Nohl T， et al. Astronomicallypacedclimate and carbon-cycle feedbacks in the lead-up to theLate Devonian Kellwasser Crisis[J]. Climate of the Past， 2023，doi： 10. 5194/cp-2023-58.

［59］ 劉曉東，石正國. 歲差對亞洲夏季風(fēng)氣候變化影響研究進展[J]. 科學(xué)通報，2009，54（20）：3097-3107.［Liu Xiaodong， ShiZhengguo. Effect of precession on the Asian summer monsoonevolution： A systematic review[J]. Chinese Science Bulletin，2009， 54（20）： 3097-3107.］

［60］ Feulner G. Formation of most of our coal brought Earth close toglobal glaciation[J]. Proceedings of the National Academy ofSciences of the United States of America， 2017， 114（43）： 11333-11337.

基金項目：國家自然科學(xué)基金項目（41690131）［Foundation： National Natural Science Foundation of China， No. 41690131］