湖相泥頁巖地層米氏旋回測井識別及環(huán)境響應(yīng)特征 以渤海灣盆地濟陽坳陷東營凹陷樊頁1井Es4scs為例

自Milankovitch運用地球軌道旋回理論解釋第四紀冰期和間冰期，進而證實第四系中存在米蘭科維奇旋回（米氏旋回）以來［1］，旋回地層學便廣泛應(yīng)用于地質(zhì)年代學研究［2］。旋回地層學的基本理論是地球軌道變化引起全球氣候變化，其將由地球軌道驅(qū)動力造成的旋回性地層記錄稱為米氏旋回，旋回中所蘊含的周期信息在地層中表現(xiàn)為沉積結(jié)構(gòu)、沉積構(gòu)造、巖性和巖相的變化，這些變化會在古氣候替代性指標中得到體現(xiàn)［3］，而運用數(shù)學方法便能提取替代指標中的周期信息［4］。同時，旋回地層學提供了一種全新的定年方式，通過研究沉積旋回中與天文軌道有關(guān)的周期來確定地層單位的延續(xù)時間和年齡。前人研究認為，細粒沉積地層的高頻旋回受米氏旋回的影響而形成［5-7］，但涉及軌道周期環(huán)境響應(yīng)特征的研究卻相對較少，因此，亟需開展湖相細粒沉積地層米氏旋回研究，建立較為精確的年代地層格架，在此基礎(chǔ)上進一步探討軌道周期對細粒沉積過程以及泥頁巖油氣富集規(guī)律的控制。

近年來，諸多學者嘗試運用各種數(shù)據(jù)和方法［8-15］來提取地層中的米氏旋回：運用自然伽馬能譜數(shù)據(jù)開展滑動平均濾波和小波變換分析，對碳酸鹽巖高頻層序進行了研究［6］；運用野外露頭和磁化率等數(shù)據(jù)研究了華南二疊紀地層中記錄的米氏旋回［16］；運用自然伽馬測井數(shù)據(jù)開展旋回分析，建立了松遼盆地白堊系青山口組的天文年代標尺［17］；通過研究早侏羅世普林斯巴期碳循環(huán)，發(fā)現(xiàn)Fe元素濃度時間序列具有明顯的米氏旋回［18］；運用意大利北部La Vedova剖面磁化率和地球化學測試數(shù)據(jù)，校正了地中海地區(qū)中新世的天文年代格架［14］。綜上所述，目前國內(nèi)外學者多采用測井、磁化率和地球化學測試分析數(shù)據(jù)，運用時頻分析方法對海陸相沉積地層的旋回性進行研究，進而計算旋回的持續(xù)時間和平均沉積速率。鑒于此，為了更好地研究湖相泥頁巖天文旋回周期以及環(huán)境響應(yīng)特征，本文以渤海灣盆地濟陽坳陷東營凹陷樊頁1井古近系沙河街組四段上亞段純上次亞段（Es4scs）為例，在前人研究基礎(chǔ)上，從天文地層學理論出發(fā)，綜合測井數(shù)據(jù)和地球化學資料，探索基于測井信息的細粒沉積地層米氏旋回識別，并首次運用有機質(zhì)豐度和表征古氧化-還原性的指標來探討軌道周期的響應(yīng)特征。旨在能為細粒沉積地層的精細劃分與對比、精細地質(zhì)定年、建立高精度的層序地層格架提供研究思路和方法，為研究天文軌道周期控制下的環(huán)境響應(yīng)特征提供借鑒。

1 地質(zhì)背景

渤海灣盆地位于中國東部，是中朝準地臺經(jīng)過古生代沉積并在印支、燕山期運動的基礎(chǔ)上發(fā)展起來的中-新生代斷陷盆地［19-20］。東營凹陷是渤海灣盆地濟陽坳陷內(nèi)的典型中-新生代陸相斷陷湖盆，其東部為青坨子凸起，南部為魯西隆起與廣饒凸起，西部為青城凸起，北部為陳家莊凸起和濱縣凸起（圖1a）［21-22］。

圖1

圖1   東營凹陷構(gòu)造地質(zhì)簡圖和地層綜合柱狀圖［21-25］

a.東營凹陷構(gòu)造帶示意圖； b.古近系綜合柱狀圖

Fig.1   Geological sketch and composite stratigraphic column of the Dongying Sag ［21-25］

東營凹陷內(nèi)發(fā)育3個正向構(gòu)造帶（北部陡坡帶、中央背斜帶和南部緩坡帶）和4個負向構(gòu)造帶（利津洼陷、牛莊洼陷、博興洼陷和民豐洼陷），凹陷東西長約90 km，南北寬約65 km，總面積約5 700 km2。古近系沙河街組四段上亞段為半深湖-深湖沉積，形成了大套細粒沉積，主要巖性為泥質(zhì)灰?guī)r、灰質(zhì)泥巖和泥巖等（圖1b）［23-26］。該段地層完整地記錄了湖泊環(huán)境和氣候變化的信息，是開展米氏旋回分析的理想載體。

2 米氏旋回識別

本次米氏旋回研究的總體思路為運用數(shù)學方法（頻譜分析、小波變換、功率譜估計和濾波分析）分析時間序列（自然伽馬測井數(shù)據(jù)）中蘊含的周期性旋回信息，再通過米氏周期計算時間跨度，并進一步計算出沉積速率。頻譜分析是將時域信號轉(zhuǎn)換至頻域信號進而對周期性現(xiàn)象進行分析，主要用來初步判斷測井信號中的優(yōu)勢頻率成分是否受軌道周期的影響；再通過小波變換進一步分析這些頻率信息是否受某一具體軌道周期的驅(qū)動，識別地層中保留的米氏旋回類型；最后通過功率譜估計準確讀取優(yōu)勢旋回對應(yīng)的尺度值，得出不同尺度下的功率譜，進而提取不同級別的旋回。運用這些方法能夠提取與地球軌道參數(shù)相關(guān)的地層響應(yīng)信息，進而能夠確定百萬年及以下級別的地層年代［27-30］。

2.1 古近紀地球軌道要素理論周期及特征

米蘭科維奇理論是描述天文因素變化導致地球軌道三要素（偏心率、地軸傾斜度和歲差）周期性變化的方法［31］，而地球軌道參數(shù)的變化是迄今為止地質(zhì)尺度上唯一可以精確定量計算的周期現(xiàn)象，可以用作地質(zhì)計時標尺［32］。前人運用La 2004方案計算了北緯38°附近50～30 Ma期間的偏心率、斜率與歲差變化的理論值［33］，采樣間隔為l kyr。據(jù)此計算出濟陽坳陷50～30 Ma期間的偏心率、斜率和歲差周期之間比例關(guān)系（表1）。天體運動規(guī)律的固定性使得各級天文周期之間的比例關(guān)系相對穩(wěn)定，因此，在穩(wěn)定的沉積環(huán)境中各級旋回厚度比值與相應(yīng)軌道周期之間的比值是對應(yīng)的，故可將這些比例關(guān)系作為米氏旋回識別的依據(jù)。

表1   50～30 Ma地球軌道參數(shù)周期及其比值關(guān)系

Table 1  Plot showing the earth orbital parameter periodicities in the past 50 to 30 Ma and their ratios

注：周期比值來源于參考文獻［33］。

新窗口打開| 下載CSV

2.2 米氏旋回分析

2.2.1 頻譜分析

本文借助了Matlab提供的一維連續(xù)小波工具箱，對樊頁1井Es4scs（深度為3 251~3 441 m）自然伽馬測井數(shù)據(jù)開展消噪和去趨勢化的預(yù)處理。

運用Past 3.0統(tǒng)計分析軟件，選擇Spectral analysis模塊中的REDFIT程序，對預(yù)處理后的測井數(shù)據(jù)開展頻譜分析［28］，從而得出自然伽馬測井信號頻譜分析圖（圖2）。提取圖中的優(yōu)勢頻率成分，運用頻率與其對應(yīng)的周期旋回厚度呈倒數(shù)關(guān)系求出相應(yīng)的地層厚度，再通過地層記錄中的米蘭科維奇性質(zhì)判斷所研究地層是否存在受軌道周期驅(qū)動的沉積旋回。

圖2

圖2   東營凹陷樊頁1井Es4scsGR曲線Redfit頻譜分析

Fig.2   Redfit spectral analysis of GR curves of Es4scs in Well Fanye 1 in the Dongying Sag

樊頁1井Es4scs自然伽馬測井曲線的頻譜分析結(jié)果表明，研究層段具有明顯的周期性，優(yōu)勢頻率0.031 574，0.107 350，0.246 270，0.315 740，0.606 210和0.681 990 m-1所對應(yīng)的主要旋回厚度分別為31.672，9.315，4.061，3.167，1.650和1.466 m（圖2），其旋回厚度（m）之間的比值為31.672∶9.315∶4.061∶3.167∶1.650∶1.466=21.600∶6.353∶2.769∶2.160∶1.125∶1，該比值與古近紀米氏旋回理論軌道周期（kyr）之間的比例405.00∶124.22∶51.68∶39.76∶220∶18.82十分接近，即31.672 m和9.315 m的旋回厚度分別對應(yīng)405.00 kyr偏心率長周期（E1）和124.22 kyr偏心率短周期（E2），4.061 m和3.167 m的旋回厚度對應(yīng)于51.68 kyr斜率周期（O1）和39.76 kyr斜率周期（O2），而1.650 m和1.466 m的旋回厚度對應(yīng)22.00 kyr歲差周期（P1）和18.82 kyr歲差周期（P2）。因此，樊頁1井Es4scs極有可能受米氏旋回控制。

2.2.2 小波變換

由于軌道信號進入沉積系統(tǒng)并不是一個線性過程，可能伴隨很多其他雜信號，因此通過頻譜分析不能完全判斷研究層段是否存在明顯的米氏旋回，還需要進一步對其進行小波分析來判斷，小波分析能夠識別出測井曲線中不同頻率的曲線旋回［34-35］。本次運用Matlab提供的一維連續(xù)小波工具箱中的Morlet小波作為子波對研究層段的自然伽馬測井數(shù)據(jù)信號開展小波分析。圖3a—d為樊頁1井Es4scs自然伽馬測井信號小波時頻能譜圖，從圖中可以看出，不同成分的頻率被分解到不同的尺度范圍上，高頻成分對應(yīng)于小尺度值，低頻成分則對應(yīng)大尺度值，總的來說不同成分的頻率對應(yīng)于某幾個特定的尺度值，即說明頻率的成分相對單一穩(wěn)定，進一步證實了研究層段受明顯的米氏旋回控制。

圖3

圖3   東營凹陷樊頁1井Es4scs一維連續(xù)小波處理結(jié)果（尺度值=256）

Fig.3   Results of one?dimensional continuous wavelet transform for Es4scs in Well Fanye 1 in the Dongying Sag （scale value： 256）

2.2.3 功率譜估計

為了進一步判斷這些相對單一穩(wěn)定的頻率成分是否受某一具體軌道周期控制，在小波分析的基礎(chǔ)上開展功率譜估計，將不同尺度的能譜值累加做出相應(yīng)的模極值圖。在Office Excel表格中自行定義公式y=｛SUM［ABS（1∶n）］｝/n，式中n為第n列處的數(shù)據(jù)值，運用公式得出3個明顯的模極大值點，其對應(yīng)的小波尺度值分別為30，51和163，其中，尺度值163和51的比值3.20與E1和E2的比值3.26非常相近。結(jié)合頻譜分析結(jié)果，可以確定該套地層受明顯的偏心率周期的控制，即尺度值163和51所對應(yīng)的優(yōu)勢頻率分別受E2和E1的驅(qū)動。故兩個尺度值所對應(yīng)的小波曲線可理解為E1和E2的旋回曲線，用以代表研究層段地層的周期旋回曲線。通過旋回曲線圖（圖4）可以得到研究層段保存的E1和E2旋回個數(shù)分別為6個和22個，兩者比值3.67∶1也與405.00 kyr∶124.22 kyr比值3.26∶1相近，進一步證實了研究層段主要受E1和E2的驅(qū)動，同時也受到O2和P1的控制。

圖4

圖4   東營凹陷樊頁1井Es4scs米氏旋回地層格架

Fig.4   Milankovitch cyclic stratigraphic framework of Es4scs in Well Fanye 1 in the Dongying Sag

2.2.4 米氏旋回周期

根據(jù)頻譜分析、小波變換和功率譜估計的結(jié)果，樊頁1井Es4scs的沉積過程受E1，E2，O2和P1的控制。對頻譜分析結(jié)果中的不同主峰進行帶通濾波處理，濾掉時間序列中與米氏旋回無關(guān)的高頻和低頻信息，提取與偏心率、斜率和歲差相關(guān)的天文軌道信號［36］。本次通過帶通濾波，將主要峰值頻率信息提取出來，用于建立不同級別的地層格架，同時獲得全井段受不同軌道周期控制的地層沉積時間，用于地層沉積速率的計算。根據(jù)濾波結(jié)果，研究得出樊頁1井Es4scs發(fā)育約6個E1，22個E2，65個O2及110個P1（圖4）。

3 “浮動”天文年代標尺的建立

利用Laskar關(guān)于地球軌道參數(shù)解的研究成果生成地球軌道要素的理論曲線［33］，將主峰值的濾波分析結(jié)果與理論曲線進行對比，發(fā)現(xiàn)E2濾波曲線與理論偏心率曲線的頻率相似，且兩者具有很好的對應(yīng)關(guān)系。因此，本次采用旋回信號明顯的E2建立東營凹陷樊頁1井Es4scs的“浮動”天文年代標尺。以理論偏心率周期曲線為目標曲線，以從自然伽馬測井曲線中提取的E2濾波曲線為調(diào)諧曲線，結(jié)合前人依據(jù)古地磁和火山巖同位素測年確定的沙四上亞段頂部的絕對年齡42 Ma［37-41］，在時間點42 Ma的控制下，利用每兩條調(diào)節(jié)線之間的E2，最終劃分出22個周期，以沙四上亞段頂部年齡開始向下計算，計算得出每個周期界面的地層年齡值，建立了樊頁1井Es4scs的“浮動”天文年代標尺（圖4）。研究層段共記錄了22個E2，由此計算得出樊頁1井Es4scs的沉積時間大致為2.73 Myr，與古地磁、火山巖同位素測定和地層對比得到的沙四上亞段（包括純上次亞段和純下次亞段）持續(xù)時間3.00 Myr是相對匹配的，亦即沙四上亞段底部的火山巖同位素年齡為45 Ma［38-39，42］，表明本次研究識別出的米氏旋回是準確可靠的，據(jù)此進一步得出沙四上亞段純上次亞段底部的絕對地質(zhì)年齡為44.73 Ma。

本次研究實現(xiàn)了湖相細粒沉積地層米蘭科維奇周期的識別，為地層精細劃分和對比提供了一種有效的手段，并建立了樊頁1井Es4scs的地質(zhì)年代浮動標尺。該天文標尺的建立可用于精確估計樊頁1井Es4scs中每一深度位置的地質(zhì)年齡，而絕對年齡可以為旋回地層學研究提供初始年代約束，同時能夠檢驗其合理性［42］。在年代標尺建立的基礎(chǔ)上研究了沉積速率的垂向變化，運用E2相鄰波峰間地層厚度與沉積持續(xù)時間（沙四上亞段純上次亞段-124.22 kyr），計算得出樊頁1井Es4scs的平均沉積速率為0.069 m/kyr（圖4）。

為了使研究結(jié)果更加可信，在識別樊頁1井沙四上亞段純上次亞段天文旋回的基礎(chǔ)上，進一步分析了樊頁1井沙三下亞段及牛頁1井沙四上亞段純上次亞段存在的米氏旋回，其中，樊頁1井沙三下亞段共識別出16個124.22 kyr偏心率短周期和4.5個405.00 kyr偏心率長周期，由此進一步計算出沉積時間大致為1.905 Myr，平均沉積速率可估算為0.105 m/kyr［27］；綜合牛頁1井沙四上亞段純上次亞段自然伽馬測井數(shù)據(jù)、深度、巖性、尺度值以及小波能譜圖，繪制該井的旋回地層格架圖（圖5）。由圖可知，沙四上亞段純上次亞段地層保存明顯的偏心率、斜率和歲差周期，研究層段識別出約5.7個405.00 kyr長偏心率周期、24個95.24 kyr短偏心率周期、53個39.76 kyr斜率周期和110.5個18.82 kyr歲差周期，該段地層厚度為151 m（深度3 316～3 467 m），由此計算出牛頁1井沙四上亞段純上次亞段的沉積時間大致為2.195 Myr，平均沉積速率估算為0.069 m/kyr。

圖5

圖5   東營凹陷牛頁1井Es4scs米氏旋回地層格架

Fig.5   Milankovitch cyclic stratigraphic framework of Es4scs in Well Niuye 1 in the Dongying Sag

4 軌道周期的環(huán)境響應(yīng)特征

天體行星影響地球繞太陽公轉(zhuǎn)和其自轉(zhuǎn)運動，即太陽系中地球與太陽、月亮和其他行星之間的萬有引力作用導致地球繞太陽的公轉(zhuǎn)運動和地球的自轉(zhuǎn)軸發(fā)生（準）周期性的攝動，表現(xiàn)為地球軌道參數(shù)偏心率、斜率和歲差發(fā)生近似周期性的變化［43］。由天文軌道控制的日-地距離和地軸的傾斜度影響了地球表面接收日照量的強度，并引起氣候的周期性變化，這些變化在地層中以巖性、地球物理和地球化學參數(shù)的旋回性表現(xiàn)出來［44-46］。因此，可通過反映沉積環(huán)境特征的有機質(zhì)豐度及元素地球化學指標探究軌道周期的環(huán)境響應(yīng)特征。

研究數(shù)據(jù)采用中國石化勝利油田勘探開發(fā)研究院測試的間距約0.1 m的元素地球化學分析數(shù)據(jù)1 894個（深度3 251.67～3 440.97 m），這種連續(xù)的數(shù)據(jù)既保證了數(shù)據(jù)的有效性又可作為最優(yōu)替代性指標開展時頻分析。有機質(zhì)豐度指標主要包括總有機碳含量（TOC）和生烴潛量（S1+S2），數(shù)據(jù)同樣來源于中國石化勝利油田勘探開發(fā)研究院測試的187個TOC和S1+S2數(shù)據(jù)，測試間隔約為1 m（深度3 251.67～3 440.58 m），每個旋回內(nèi)基本能測試兩個數(shù)據(jù)點，滿足了替代性指標時頻分析的最低數(shù)據(jù)間隔標準，這樣既不會造成不必要的時間和資金浪費，也可以保證獲得真實的旋回信息。

4.1 有機質(zhì)豐度的軌道周期響應(yīng)

氣候變化導致沉積環(huán)境的不同往往直接影響了有機化合物的聚集與保存，并在TOC指標上呈高低起伏的變化特征，因此，天文軌道周期的變動往往記錄在了TOC和S1+S2等指示有機質(zhì)豐度指標的變化中。本次研究首次在研究區(qū)探索性地運用有機質(zhì)豐度指標來探討有機質(zhì)對軌道周期的響應(yīng)特征，以TOC和S1+S2作為有機質(zhì)豐度指標參數(shù)對其開展米氏旋回研究，探索軌道周期與有機碳含量和生烴潛量的響應(yīng)關(guān)系。

樊頁1井Es4scs泥頁巖的有機碳含量普遍較高，TOC平均值在2.0 %以上。將該層段有機碳含量數(shù)據(jù)繪制成TOC頻率分布直方圖（圖6），TOC占比最大的區(qū)間分布在1 %～2 %和2 %～3 %，比例分別達到了43.8 %和40.0 %，而含量大于3.0 %的比例為10.8 %，僅有5.4 %的TOC小于1.0 %。同時，根據(jù)陸相烴源巖評價標準（SY/T5735—1995），結(jié)合有機質(zhì)豐度評價的主要指標S1+S2，發(fā)現(xiàn)研究區(qū)泥頁巖的有機質(zhì)豐度較高，主要落于好-最好區(qū)間（圖7），表明在樊頁1井Es4scs沉積時期，TOC高（介于1 %～2 %），有機質(zhì)聚集較為集中。

圖6

圖6   東營凹陷樊頁1井Es4scsTOC分布

Fig.6   Histogram showing TOC contents of Es4scs in Well Fanye 1 in the Dongying Sag

圖7

圖7   東營凹陷樊頁1井Es4scs有機質(zhì)豐度分析判別圖（據(jù)SY/T 5735—1995標準）

Fig.7   Organic matter abundance chart of Es4scs in Well Fanye 1 in the Dongying Sag （according to SY/T 5735-1995）

據(jù)TOC測試數(shù)據(jù)，中、下部地層（E1中的2—3旋回、E2中的5—12旋回、P1中的20—60旋回）的TOC高且達到峰值，向上變化趨于平穩(wěn)。對比TOC旋回分析曲線與偏心率周期變化曲線，可見偏心率周期在研究層段的中、下部地層幅值幅度變化較大，而向上周期曲線變化平緩，說明偏心率周期振幅變化趨勢與TOC變化趨勢具有一致性，即偏心率較高時期沉積層段的平均TOC也相對較高，這主要是由于偏心率增大時地球上接收的日照量增大，氣候變得暖濕，物源供給充足，生物有機質(zhì)豐富，該時期有利于發(fā)育細粒且富有機質(zhì)的沉積地層［47-49］，故偏心率值較高時期沉積層段是潛在頁巖油氣的有利儲集層段。

以上研究可知，研究層段的米氏旋回研究成果與TOC在地層中的縱向分布具有一定的相關(guān)關(guān)系。為了進一步確定這種相關(guān)關(guān)系是否受到天文軌道周期的控制，探索其沉積記錄中蘊含的天文軌道周期信號，本文選取TOC和S1+S2數(shù)據(jù)作為原始數(shù)據(jù)，運用Past 3.0統(tǒng)計分析軟件對數(shù)據(jù)開展頻譜分析。結(jié)果表明，研究層段的TOC和S1+S2數(shù)據(jù)顯示出明顯的周期性（表2；圖8a—d），即TOC和S1+S2頻譜分析得到的優(yōu)勢旋回厚度比值與測井數(shù)據(jù)得出的結(jié)果一致，表明研究層段的沉積記錄中蘊含著軌道周期信息，且主要受控于偏心率周期，也進一步證實了地球軌道參數(shù)的周期變化影響了有機質(zhì)豐度的變化。

表2   東營凹陷樊頁1井Es4scsTOC，S1+S2及微量元素頻譜分析旋回厚度比值與理論軌道周期比值對比

Table 2  Cycle thickness ratios from spectrum analysis of TOC， S1+S2 and trace elements， and theoretical orbital periodicity ratios of Es4scs in Well Fanye 1 in the Dongying Sag

新窗口打開| 下載CSV

圖8

圖8   東營凹陷樊頁1井Es4scs地層TOC，S1+S2及微量元素頻譜分析

E1.405 kyr偏心率長周期； E2.124.22 kyr偏心率短周期

Fig.8   Spectrum analysis of TOC， S1+S2 and trace elements of Es4scs in Well Fanye 1 in the Dongying Sag

a.TOC； b.S1+S2； c.V/（V+Ni）； d.Ni/Co

4.2 古氧化-還原性的軌道周期響應(yīng)

有機質(zhì)的富集和保存受古氧化-還原環(huán)境的影響，而有機質(zhì)在不同的環(huán)境下經(jīng)受著消耗和堆積的循環(huán)變化，為研究這種變化的驅(qū)動力是否來自軌道周期，本次研究首次在研究區(qū)探索性地運用表征古氧化-還原性指標探討軌道周期的響應(yīng)特征，即選取兩組微量元素比值V/（V+Ni）和Ni/Co作為替代性指標來分析目的層段古氧化-還原環(huán)境的變化規(guī)律，并探討其軌道周期的響應(yīng)特征。

前人研究表明［50-53］，當V/（V+Ni）大于0.60時，主要為缺氧的還原環(huán)境，當0.46<V/（V+Ni）<0.60時，認為是貧氧環(huán)境，當V/（V+Ni）<0.46時，表示為富養(yǎng)環(huán)境；Ni/Co<5時，表明是氧化環(huán)境，當5<Ni/Co<7時，表示為貧氧環(huán)境，當Ni/Co>7時，顯示為缺氧環(huán)境。從微量元素測試分析結(jié)果可知，研究層段V/（V+Ni）值介于0.24～1.00，均值為0.82，Ni/Co值介于0.01～18.35。為了使微量元素比值對環(huán)境的指示作用更為直觀，利用有機碳含量和兩種微量元素含量比值建立關(guān)系圖版，繪制出古氧化-還原條件判別圖（圖9），從圖中可以看出，樣品點多落于貧氧的區(qū)域。因此，研究層段整體為缺氧的還原環(huán)境，有機質(zhì)含量高的層段多為缺氧的沉積環(huán)境。

圖9

圖9   東營凹陷樊頁1井Es4scs微量元素氧化-還原條件判識圖版

a.TOC-V/（V+Ni）； b.TOC-Ni/Co

Fig.9   Identification of redox conditions by trace elements of Es4scs in Well Fanye 1 in the Dongying Sag

為探究軌道周期對古氧化-還原環(huán)境的控制作用，在研究區(qū)探索性地運用兩組微量元素比值V/（V+Ni）和Ni/Co作為原始數(shù)據(jù)開展頻譜分析（圖8）。分析結(jié)果顯示，兩組數(shù)據(jù)中均識別出明顯的米氏旋回周期（表2），主要峰值的旋回厚度均為31.33 m和9.40 m，其旋回厚度比值為3.33∶1，該比值與理論軌道周期405 kyr∶124.22 kyr=3.26∶1比值很接近，因此，V/（V+Ni）和Ni/Co的主要峰值均被認為是由E1和E2引起的。通過分析可知，旋回厚度31.33 m對應(yīng)E1，9.40 m對應(yīng)E2。頻譜分析結(jié)果表明，兩組微量元素均顯示出一致的規(guī)律性和相似的天文周期比值，該比值與測井數(shù)據(jù)得出的結(jié)果一致，從而可知微量元素比值的周期性記錄反應(yīng)了目的層段古氧化-還原環(huán)境的變化受到天文軌道周期的影響。

通過分析偏心率周期的變化可知，暖濕氣候有利于生物大量繁殖，即偏心率值較高的層段平均TOC相對較高，對氣候環(huán)境較敏感的方解石相對含量是古環(huán)境信息的重要來源，氣候偏干旱時水體蒸發(fā)作用較強，碳酸鹽礦物（主要為方解石）含量較高，相反，當氣候相對濕潤時碳酸鹽礦物含量較低，因此可知方解石含量與偏心率值呈負相關(guān)性。上述指標的變化進一步說明了偏心率增大時氣候會向暖濕轉(zhuǎn)變，且隨著偏心率增高，TOC也增高。

綜上所述，偏心率周期通過控制TOC、礦物含量及各種表征沉積環(huán)境介質(zhì)條件的指標而對環(huán)境產(chǎn)生影響。有機質(zhì)含量高的層段多處于缺氧的還原環(huán)境，這主要是由于在缺氧的環(huán)境中，有機質(zhì)不會因與氧氣接觸而發(fā)生降解反應(yīng)，而會因缺氧被保存下來。同時，偏心率周期振幅幅度的變化對有機質(zhì)豐度變化有影響，當偏心率振幅幅度較大且處于最大值時期，TOC的振幅和數(shù)值也相應(yīng)增大，進一步說明偏心率周期振幅變化趨勢與TOC變化趨勢具有一致性，即偏心率較高時期沉積層段的平均TOC也相對較高。這主要是由于偏心率增大時地球上界接收的日照量增大，季節(jié)性增強［47，54］，氣候變得暖濕，陸源碎屑輸入使得物源供給充足，大量營養(yǎng)物質(zhì)提高了湖泊的生產(chǎn)力，水體較深且水體偏還原環(huán)境，從而使有機質(zhì)含量升高，該時期有利于發(fā)育細粒且富有機質(zhì)的沉積地層，故偏心率值較高時期沉積層段是潛在頁巖油氣的有利儲集層段。

5 結(jié)論

1） 樊頁1井Es4scs地層自然伽馬測井數(shù)據(jù)頻譜分析和小波變換分析結(jié)果表明，米氏旋回驅(qū)動控制了東營凹陷古近紀的沉積周期。米氏旋回的自然伽馬測井識別，為渤海灣盆地古近系沙河街組湖相細粒沉積地層的高精度等時旋回地層的劃分與對比提供了參考。

2） 研究層段共識別出約6個E1，22個E2，65個O2和110個P1，建立了樊頁1井Es4scs的“浮動”天文年代標尺，計算得出了樊頁1井Es4scs的沉積時間大致為2.73 Myr，平均沉積速率估算為0.069 m/kyr。米氏旋回的識別解決了細粒沉積地層的精細劃分和等時對比的難題，對于頁巖油氣勘探具有重要的地質(zhì)意義。

3） 地球軌道參數(shù)周期的變化在一定程度上影響了有機質(zhì)豐度及古氧化-還原性的變化。當偏心率振幅較大且處于最大值時期，地球整體接受日照量多，為間冰期，氣候暖濕，有利于有機質(zhì)的聚集和保存，易于形成缺氧的富有機質(zhì)沉積地層。因此，偏心率值較高的暖濕氣候背景下發(fā)育的富有機質(zhì)沉積地層，是潛在的頁巖油氣勘探的有利層段。