四川盆地中二疊統(tǒng)茅口組地層缺失量計(jì)算及成因探討

江青春 胡素云 姜 華 翟秀芬 任夢(mèng)怡 陳曉月 李秋芬 張運(yùn)波

1. 中國(guó)石油勘探開(kāi)發(fā)研究院 2. 數(shù)巖科技股份有限公司

巖溶儲(chǔ)層是全球含油氣盆地最重要的碳酸鹽巖油氣儲(chǔ)層之一?？碧綄?shí)踐表明，中國(guó)發(fā)現(xiàn)的大型海相油氣田大都與古巖溶作用有關(guān)[1-6]。古巖溶儲(chǔ)層的發(fā)育規(guī)模既受控于古巖溶地貌特征，又受控于地層缺失特征。目前，對(duì)于四川盆地中二疊統(tǒng)茅口組頂部存在地層缺失這一認(rèn)識(shí)得到了大多數(shù)學(xué)者的認(rèn)可[7-11]，但對(duì)其缺失特征僅進(jìn)行了簡(jiǎn)單的描述[12-13]。少數(shù)學(xué)者對(duì)局部的地層缺失量進(jìn)行了定量分析[14]，然而對(duì)于盆地內(nèi)部地層缺失量均未開(kāi)展過(guò)深入的分析，特別是對(duì)于地層缺失的成因研究更為薄弱，多數(shù)學(xué)者認(rèn)為是由于東吳運(yùn)動(dòng)的抬升造成的地層整體隆升，形成地層缺失帶，僅有李旭兵等少數(shù)學(xué)者認(rèn)為是由于茅口晚期東吳運(yùn)動(dòng)Ⅰ幕發(fā)生大規(guī)模的海退導(dǎo)致茅口組頂部碳酸鹽巖暴露，遭受廣泛的淡水淋濾風(fēng)化作用，形成地層缺失，但未給出足夠的證據(jù)[15]。為了解決這些科學(xué)問(wèn)題，本文在定性分析和描述茅口組地層缺失特征的基礎(chǔ)上，利用米蘭科維奇高頻旋回分析方法定量計(jì)算了盆地內(nèi)茅口組的頂部剝蝕量，并探討了其成因，為精確預(yù)測(cè)四川盆地茅口組巖溶儲(chǔ)層分布及下一步勘探選區(qū)奠定了基礎(chǔ)。

1 茅口組頂部地層缺失定性識(shí)別

筆者在2012年發(fā)表的《四川盆地茅口組風(fēng)化殼巖溶古地貌及勘探選區(qū)》[16]文章中提到：茅口組由川西南向川東北及川西北方向地層缺失程度逐漸增加；東吳運(yùn)動(dòng)的構(gòu)造抬升作用致使茅口組頂部地層缺失；地層的風(fēng)化及淋濾作用使大氣淡水沿著斷裂及裂縫形成差異性溶蝕，形成巖溶儲(chǔ)層。隨著后期研究工作的逐步深入，筆者利用四川盆地11條格架剖面對(duì)地層進(jìn)行了精細(xì)對(duì)比，并利用盆地內(nèi)二維地震測(cè)線對(duì)茅口組地層厚度變化進(jìn)行了定量計(jì)算，然后通過(guò)鉆井與地震結(jié)合將茅四段和茅三段的地層尖滅線進(jìn)行了精細(xì)刻畫，分析結(jié)果表明，茅口組殘余厚度由雅安—宜賓地區(qū)向?yàn)o州—南充地區(qū)再到川北宣漢—旺蒼地區(qū)逐漸減薄。茅四段僅在川西南部地區(qū)和川東北部發(fā)育，川西南部主要發(fā)育于宜賓—雅安—江油一帶，川東北部?jī)H石柱地區(qū)殘存；茅三段在盆地內(nèi)分布較廣，僅在萬(wàn)州—宣漢—旺蒼一帶普遍缺失，但該區(qū)域內(nèi)的渡口及旺蒼大兩鄉(xiāng)有零星殘余分布（圖1）。川北地區(qū)缺失茅三段的區(qū)域，茅口組殘余厚度很薄，筆者在文初提到的文章 中認(rèn)為該地區(qū)在茅口組末期為巖溶高地，巖溶地貌較高，但該區(qū)在吳家坪期確為深水盆地相沉積。因此為了解決這一矛盾，筆者利用米蘭科維奇旋回分析計(jì)算了該區(qū)地層剝蝕量，通過(guò)四個(gè)方面的論證分析，認(rèn)為茅口末期受冰川海平面下降影響，具有由西南向東北逐漸降低的巖溶地貌格局，基于地層厚度的差異將其劃分為侵蝕上斜坡和侵蝕下斜坡（圖2），并在本文第4部分中對(duì)其成因進(jìn)行了論述和分析。

圖1 茅口組地層尖滅線分布及侵蝕巖溶地貌圖

圖2 四川盆地茅口組巖溶地貌發(fā)育模式圖

2 地層剝蝕量定量計(jì)算方法

碎屑巖地層剝蝕量的恢復(fù)方法較多，也較為成熟，常用的方法包括趨勢(shì)地層對(duì)比法、鏡質(zhì)組反射率法、聲波時(shí)差法、沉積波動(dòng)過(guò)程分析法等[17-19]。通過(guò)對(duì)這些方法的梳理發(fā)現(xiàn)，因其原理或基于砂泥巖地質(zhì)模型或獨(dú)特的適用條件，很難用于碳酸鹽巖地層剝蝕量的恢復(fù)。為了解決這一難題，不斷有學(xué)者探索碳酸鹽巖地層缺失量的恢復(fù)方法。多名學(xué)者的研究證實(shí)，米蘭科維奇旋回分析能夠有效記錄地層的沉積旋回特征，對(duì)于地層部分缺失的地區(qū)可以嘗試?yán)玫貙映练e速率差異及剝蝕量差異的相互對(duì)比與驗(yàn)證，進(jìn)行碳酸鹽巖地層剝蝕量恢復(fù)[20-23]，為定量計(jì)算海相碳酸鹽巖地層剝蝕量提供新的思路。

米蘭科維奇旋回分析理論已經(jīng)在旋回地層學(xué)研究中得到較為普遍的應(yīng)用[24]，其表現(xiàn)形式為沉積地層具有一定程度的旋回性和韻律性。國(guó)內(nèi)外學(xué)者一般運(yùn)用釷鉀比（Th/K）曲線分析碳酸鹽巖沉積環(huán)境下的相對(duì)海平面的變化，但由于自然伽馬能譜測(cè)井系列價(jià)格較高，很多鉆井不對(duì)此類曲線開(kāi)展測(cè)井工作。通過(guò)對(duì)比分析同一鉆井釷鉀比曲線和自然伽馬曲線的米蘭科維奇旋回特征，實(shí)踐分析表明在釷鉀比曲線缺乏的情況下，可以采用自然伽馬曲線代替分析米蘭科維奇旋回，其誤差在可以接受的范圍內(nèi)。米蘭科維奇旋回分析的主要工作是開(kāi)展頻譜分析，其本質(zhì)就是認(rèn)為測(cè)井曲線是由各種地質(zhì)要素的綜合作用在時(shí)間尺度（或深度尺度）上形成的一個(gè)地層規(guī)律性變化的綜合疊加信號(hào)，運(yùn)用傅里葉變換將這一疊加信號(hào)（測(cè)井曲線）從深度尺度（或時(shí)間尺度）轉(zhuǎn)換到頻率尺度而形成頻譜曲線。其分析流程如下：①原始曲線分析，去除異常值；②利用小波變換對(duì)固定采樣間隔的自然伽馬能譜曲線或者自然伽馬曲線進(jìn)行分解，去除曲線的低頻背景和高頻噪音信號(hào)，然后對(duì)包含規(guī)律性低頻成分的測(cè)井?dāng)?shù)據(jù)進(jìn)行重構(gòu)；③運(yùn)用傅里葉變化方法將重構(gòu)后的數(shù)據(jù)從時(shí)間域轉(zhuǎn)換到頻率域，形成一維的頻譜曲線；④根據(jù)頻譜曲線計(jì)算出來(lái)的旋回周期參數(shù)與米蘭科維奇天文參數(shù)進(jìn)行對(duì)比分析；⑤根據(jù)對(duì)比參數(shù)及沉積速率求取地層剝蝕量。張運(yùn)波及Berger等學(xué)者對(duì)米蘭科維奇旋回分析法的具體原理及詳細(xì)過(guò)程在相關(guān)文獻(xiàn)中進(jìn)行了詳細(xì)介紹和論述[25-27]。

基于許國(guó)明、江青春等人的研究成果認(rèn)為，四川盆地茅口組沉積相為緩坡模式?？傮w而言，除盆地西部及北部茅口組相帶發(fā)育深緩坡外，盆地主體范圍內(nèi)茅口組主要為淺緩坡與中緩坡沉積，即茅口組沉積時(shí)期，四川盆地整體水體較淺，盆地范圍內(nèi)的沉積水深變化不大，其結(jié)果是沉積物沉積速率變化很小，沉積地層厚度在全盆地范圍內(nèi)厚度變化相對(duì)較小，在區(qū)域地震剖面上也表現(xiàn)出橫向穩(wěn)定，厚度變化小。

3 茅口組頂部地層缺失量定量計(jì)算

在上述茅口組沉積背景及沉積速率分析的基礎(chǔ)上，結(jié)合前期區(qū)域地層對(duì)比的認(rèn)識(shí)，分別在地層保存完整區(qū)、弱缺失區(qū)、中等缺失區(qū)和較強(qiáng)缺失區(qū)選擇了4口典型井——宮深1井、板東12井、高石1井和河壩1井。其中宮深1井茅口組地層發(fā)育完整，主要有3方面證據(jù)：①牙形石帶發(fā)育較為完整；②發(fā)育時(shí)限與地層沉積時(shí)限接近（見(jiàn)后文）；③區(qū)域地層對(duì)比揭示其茅四段基本完整。因此可以作為參照標(biāo)準(zhǔn)井計(jì)算其他井茅口組頂部地層剝蝕量。通過(guò)對(duì)4口井茅口組的自然伽馬能譜曲線和自然伽馬曲線按照每米8個(gè)點(diǎn)的采樣間隔，在數(shù)據(jù)前處理的基礎(chǔ)上按照上述方法原理與步驟開(kāi)展時(shí)頻變換，得到4口井對(duì)應(yīng)曲線的一維頻譜圖（圖3），結(jié)合茅口組的地層沉積時(shí)限，開(kāi)展旋回對(duì)比分析。結(jié)果表明，宮深1井頻譜中點(diǎn)c和點(diǎn)i分別對(duì)應(yīng)的平均旋回厚度比值為1∶0.28，分別對(duì)應(yīng)于長(zhǎng)偏心率旋回（125 ka）和斜度旋回（34.94 ka）；河壩1井頻譜中點(diǎn)c和點(diǎn)d分別對(duì)應(yīng)125 ka周期的長(zhǎng)偏心率旋回和95 ka周期的短偏心率旋回；高石1井頻譜中的點(diǎn)b和點(diǎn)d別對(duì)應(yīng)125 ka周期的長(zhǎng)偏心率旋回和95 ka周期的短偏心率旋回；板東12井頻譜中的點(diǎn)b和c點(diǎn)分別對(duì)應(yīng)125 ka周期的長(zhǎng)偏心率旋回和95 ka周期的短偏心率旋回。

圖3 典型井茅口組FFT頻譜分析曲線圖

宮深1井茅口組自然伽馬能譜測(cè)井?dāng)?shù)據(jù)的頻譜顯示，四川盆地中二疊統(tǒng)沉積地層記錄的米蘭科維奇旋回信息較為明顯，長(zhǎng)偏心率旋回（125 ka）和短偏心率旋回（95 ka）為高頻層序發(fā)育的主控因素。根據(jù)前期研究成果結(jié)合地質(zhì)年代表分析茅口組沉積時(shí)間為6.29 Ma，結(jié)合宮深1井的頻譜分析，四川盆地茅口組大約發(fā)育50.3 ka短偏心率旋回。宮深1井茅口組地層測(cè)井曲線的頻譜圖顯示，125 ka對(duì)應(yīng)的長(zhǎng)偏心率旋回信號(hào)最強(qiáng)，最適用于計(jì)算其他相關(guān)的米氏旋回參數(shù)（表1）。與長(zhǎng)偏心率旋回周期對(duì)應(yīng)的平均旋回厚度在宮深1井、河壩1井、高石1井和板東12井均為6.65 m。河壩1井的測(cè)井綜合解釋及區(qū)域地層對(duì)比表明，河壩1井茅口組頂部缺失茅四段和茅三段，驗(yàn)證了頻譜分析結(jié)果的可靠性。

茅口組頂部地層缺失現(xiàn)象的定性分析結(jié)果表明，宮深1井茅口組地層近似無(wú)剝蝕。據(jù)此，根據(jù)宮深1井茅口組測(cè)井曲線頻譜中的長(zhǎng)偏心率旋回周期，計(jì)算出四川盆地茅口組的沉積時(shí)限為（334.5 m÷6.65 m）×125 ka＝6.29 Ma。再根據(jù)宮深1井茅口組地層的發(fā)育時(shí)限和各井的平均旋回沉積厚度，計(jì)算各井茅口組頂部剝蝕量。河壩1井茅口組長(zhǎng)偏心率旋回（125 ka）對(duì)應(yīng)的平均沉積厚度為6.65 m，其頂面缺失量為（6.29 Ma/125 ka）×6.65 m－164.2 m＝170.43 m；高石1井茅口組長(zhǎng)偏心率旋回（125 ka）對(duì)應(yīng)的平均沉積厚度為6.65 m，其頂面剝蝕量為（6.29 Ma/125 ka）×6.65 m－192.6 m = 142.03 m；板東12井茅口組長(zhǎng)偏心率旋回（125 ka）對(duì)應(yīng)的平均沉積厚度為6.65 m，其頂面缺失量為（6.29 Ma/125 ka）×6.65 m－ 252.5 m＝82.13 m。

表1 中二疊世米蘭科維奇旋回周期間關(guān)系分析表

研究發(fā)現(xiàn)，宮深1井位于火山巖覆蓋區(qū)，頂部受火山巖沉積保護(hù)，茅口組無(wú)明顯缺失；高石1井地層缺失較大，頂部缺失量達(dá)142.03 m；河壩1井缺失量最大，達(dá)170.43 m；板東12井地層頂部缺失量?jī)H 82.13 m。位于中國(guó)石化礦權(quán)區(qū)內(nèi)的河壩場(chǎng)構(gòu)造的河壩1井頂部剝蝕現(xiàn)象明顯，其茅口組大約發(fā)育40.43 ka短偏心率旋回，頂部剝蝕現(xiàn)象明顯。綜合錄井報(bào)告也表明，中二疊世時(shí)期，河壩1井位于川北地區(qū)地層殘余厚度較小區(qū)，區(qū)域地質(zhì)背景揭示其為較深緩坡沉積，具有較低的沉積速率，中二疊世末期的冰期海平面下降導(dǎo)致河壩1井茅口組遭受剝蝕。根據(jù)河壩1井茅口組的平均沉積速率計(jì)算，其頂部剝蝕厚度大約為169.92 m。板東12井茅口組發(fā)育有37.98 ka短偏心率旋回，頂部缺失量達(dá)82.13 m。

采用同樣的方法，逐一對(duì)全盆地108口井茅口組層段的自然伽馬測(cè)井曲線或者自然伽馬能譜曲線（Th/K）開(kāi)展頻譜分析，確定米蘭科維奇旋回周期，計(jì)算沉積速率，并最終求取每口井茅口組頂部的地層缺失量。同時(shí)為了確保頻譜分析在深度域的可靠性，研究過(guò)程中根據(jù)最新編制的四川盆地茅口組沉積相圖，進(jìn)行了分相帶的旋回計(jì)算，對(duì)計(jì)算結(jié)果進(jìn)行了沉積速率的相互對(duì)比驗(yàn)證。川西南部地區(qū)的周公1井、漢深1井、大深1井、老龍1井處于淺緩坡相，沉積速率高，計(jì)算的沉積速率也比較接近；川中地區(qū)處于開(kāi)闊臺(tái)地內(nèi)，局部發(fā)育臺(tái)內(nèi)灘，沉積速率中等，高石2井的計(jì)算結(jié)果得到了證實(shí)；川東北地區(qū)張2井為較深水沉積，計(jì)算結(jié)果證實(shí)沉積速率相對(duì)較低，對(duì)比川西北地區(qū)的河壩1井（深水沉積），河壩1井沉積速率比張2井高。通過(guò)以上分析，利用以上原理對(duì)全盆地108井的一維傅里葉變換進(jìn)行分析，最終編制茅口組地層缺失厚度圖（圖4）。四川盆地茅口組存在兩個(gè)地層弱缺失區(qū)：①成都—宜賓一線向西南方向地層缺失相對(duì)較少，缺失厚度在0～80 m范圍內(nèi)，缺少茅四段的中上部；②重慶—萬(wàn)州一帶，地層剝蝕厚度介于40～100 m，地層保存相對(duì)比較完整。蜀南、川西北部及川東北部地區(qū)地層缺失相對(duì)嚴(yán)重，缺失厚度介于140～210 m，其中川東北部和蜀南地區(qū)地層缺失要比川中地區(qū)的地層缺失程度高。

4 茅口組頂部地層缺失成因探討

圖4 四川盆地茅口組剝蝕厚度等值線圖

對(duì)于四川盆地地層缺失的成因，多數(shù)學(xué)者認(rèn)為是由于茅口期末的東吳運(yùn)動(dòng)的構(gòu)造活動(dòng)造成 。然而，對(duì)于其運(yùn)動(dòng)性質(zhì)到底是水平擠壓的造山運(yùn)動(dòng)還是垂直升降運(yùn)動(dòng)，目前存在爭(zhēng)議。筆者的研究團(tuán)隊(duì)前期也曾提出由于東吳運(yùn)動(dòng)的擠壓抬升作用導(dǎo)致四川盆地在中二疊世末期發(fā)育一個(gè)面積可超過(guò)8×104km2的東吳期古隆起，并命名為“瀘州—通江古隆起”，古隆起斜坡帶地層暴露剝蝕是風(fēng)化殼巖溶儲(chǔ)層發(fā)育的有利地區(qū)[30]。然而最新地貌學(xué)研究結(jié)合沉積背景分析認(rèn)為茅口冰期海平面下降是導(dǎo)致地層缺失的主要原因，并明確了四川盆地由茅口早期到吳家坪期的構(gòu)造演化與地層沉積充填過(guò)程，在茅口沉積期早—中期，四川盆地主體為碳酸鹽巖緩坡沉積，由川中地區(qū)向川北地區(qū)逐漸由中緩坡演化為深斜坡相沉積。茅口組沉積后，由于全球冰期的影響，海平面急劇下降，侵蝕基準(zhǔn)面可能下降到茅口組下部，其上覆地層均遭受侵蝕，但盆地內(nèi)西南高、北東低的古構(gòu)造格局沒(méi)有發(fā)生改變。吳家坪期海平面逐漸上升，茅口組侵蝕后的古構(gòu)造格局仍得以保存，沉積環(huán)境及相帶的平面展布和茅口組相似。這一認(rèn)識(shí)合理地解釋了四川盆地北部地區(qū)茅口組地層缺失層位多與吳家坪組深水沉積兩者之間的矛盾。總體看，由盆地西南地區(qū)向川北地區(qū)可能整體表現(xiàn)為一個(gè)大型古侵蝕斜坡，侵蝕高地主要分布在盆地西南地區(qū)及川東石柱地區(qū)，侵蝕上斜坡分布在南充—瀘州一帶，侵蝕下斜坡分布在川北地區(qū)（圖1），其成因證據(jù)主要有以下三方面。

4.1 茅口末期為全球性的冰期事件

全球海平面變化曲線揭示[31]，晚古生代全球冰期始于泥盆紀(jì)末，鼎盛于石炭紀(jì)—早二疊世，可延續(xù)到中晚二疊世瓜德魯普統(tǒng)末期（茅口末期），由于冰期活動(dòng)導(dǎo)致全球海平面下降[32-35]，冰期海平面下降幅度為20～120 m[36]。二疊紀(jì)，揚(yáng)子陸塊呈“孤島”立于古特提斯洋[37]，更易于受海平面升降變化影響，由于冰期的持續(xù)作用，海平面在茅口中后期開(kāi)始下降，茅口末期快速下降到最低點(diǎn)，揚(yáng)子大部分地區(qū)準(zhǔn)平原化；吳家坪—長(zhǎng)興早期海平面緩慢上升，長(zhǎng)興中期升至高點(diǎn)，這在中—上揚(yáng)子地區(qū)二疊紀(jì)海平面變化曲線上也有所印證[38]。

4.2 川北地區(qū)牙形石化石缺失嚴(yán)重，局部存在侵蝕殘丘

受全球海平面下降影響，茅口組頂部地層缺失普遍。在野外剖面的牙形石化石帶上體現(xiàn)非常明顯，在這一時(shí)期歐美等全球大部分地區(qū)茅口組頂部缺失1～2個(gè)化石帶[38]；華南地區(qū)茅口組頂部缺失2個(gè)以上牙形化石帶。為了落實(shí)四川盆地北部地區(qū)茅口組頂部牙形石化石帶的缺失情況，對(duì)川北地區(qū)的廣元烏木沱、廣元上寺、廣元西北鄉(xiāng)、宣漢渡口等6條剖面的牙形石化石帶進(jìn)行了分析，除宣漢渡口和旺蒼大兩剖面露頭缺失兩個(gè)牙形石化石帶外，其他地區(qū)一般缺失4～6個(gè)牙形石化石帶（圖5）。表明川北地區(qū)地層缺失較其他地區(qū)缺失嚴(yán)重，但由于海平面下降的侵蝕不均一性，可能導(dǎo)致川北地區(qū)存在侵蝕殘丘。

4.3 茅口組與吳家坪組碳同位素漂移特征明顯

圖5 四川盆地川北地區(qū)茅口組典型野外剖面牙形石化石帶分布圖

導(dǎo)致海相碳酸鹽巖的δ13C值發(fā)生變化的因素較多，例如有機(jī)碳埋藏量、生物生產(chǎn)率、海洋與大氣交換及陸地物質(zhì)輸入等等，其中最重要的兩個(gè)因素是有機(jī)碳埋藏量和生物生產(chǎn)率。自然界中存在有機(jī)和無(wú)機(jī)兩大碳庫(kù)，當(dāng)其中一個(gè)碳庫(kù)同位素發(fā)生變化，就會(huì)間接導(dǎo)致另外一個(gè)碳庫(kù)同位素的變化。海洋中碳同位素的分餾作用主要由生物活動(dòng)造成，一般情況下，12C優(yōu)先為活的有機(jī)體所吸收利用，當(dāng)海洋生物生產(chǎn)率高且被埋藏時(shí)候，就會(huì)從海洋中帶走輕的碳同位素，導(dǎo)致海洋中溶解的無(wú)機(jī)碳同位素變重，引起海洋沉積物碳酸鹽巖的δ13C值增加。當(dāng)有機(jī)質(zhì)被暴露氧化，富集12C的碳析出，海洋庫(kù)溶解的無(wú)機(jī)碳同位素組成變輕，引起海洋沉積的無(wú)機(jī)碳酸鹽的δC值降低。在廣元上寺、渡口等剖面，均可見(jiàn)茅口末期無(wú)機(jī)碳同位素出現(xiàn)快速下降、吳家坪時(shí)期逐漸快速升高，預(yù)示茅口末期發(fā)生冰封海洋作用，生物生產(chǎn)率迅速降低，而到了吳家坪期時(shí)期冰川逐漸開(kāi)始消融，生物生產(chǎn)率開(kāi)始逐漸增大，導(dǎo)致同位素從－3.8‰增加到4‰。從野外剖面的同位素曲線形態(tài)可以看出，受全球性冰期作用影響，茅口末期碳同位素曲線波動(dòng)較大，吳家坪期碳同位素曲線相對(duì)穩(wěn)定，反映茅口后期沉積環(huán)境不穩(wěn)定及海平面快速下降、吳家坪期沉積環(huán)境逐漸趨于穩(wěn)定[39]（圖6）。

圖6 四川盆地及其周緣茅口組與吳家坪組牙形石化石帶及碳同位變化曲線圖（據(jù)本文參考文獻(xiàn)[40]修改）

5 結(jié)論

1）茅四段僅在宜賓—雅安—江油地區(qū)和石柱地區(qū)地層殘存，其余地區(qū)地層普遍存在缺失，且表現(xiàn)為川南—川中—川北地層缺失強(qiáng)度逐加大。

2）茅口組地層缺失量介于0～200 m，其中川西南部和川東北部地區(qū)缺失厚度介于0～60 m，蜀南、川中和川北地區(qū)地層缺失厚度介于140～200 m。

3）冰期背景、牙形石化石帶分布特征及同位素漂移特征揭示茅口末期的地層剝蝕是由于冰期海平面下降侵蝕所致，川北地區(qū)海平面下降幅度大，導(dǎo)致其缺失茅三段，但牙形石化石帶分析揭示其存在局部殘丘。

4）由于冰期海平面下降，導(dǎo)致下二疊統(tǒng)巖溶地貌繼承了西南部東北低的沉積特征，巖溶地貌由西南到川中至川北地區(qū)逐漸由侵蝕高地逐漸過(guò)渡為巖溶上斜坡和巖溶下斜坡，與吳家坪期西南高、東北低的沉積特征具有較好的一致性。

[1] 顧家裕. 塔里木盆地輪南地區(qū)下奧陶統(tǒng)碳酸鹽巖巖溶儲(chǔ)層特征及形成模式[J]. 古地理學(xué)報(bào), 1999, 1(1): 54-60.Gu Jiayu. Characteristics and evolutional model of karst reservoirs of Lower Ordovician carbonate rocks in Lunnan area of Tarim Basin[J]. Journal of Palaeogeography, 1999, 1(1): 54-60.

[2] 許效松, 杜佰偉. 碳酸鹽巖地區(qū)古風(fēng)化殼巖溶儲(chǔ)層[J]. 沉積與特提斯地質(zhì), 2005, 25(3): 1-7.Xu Xiaosong & Du Baiwei. The palaeoweathering crust-type karst reservoir rocks in carbonate rocks[J]. Sedimentary Geology and Tethyan Geology, 2005, 25(3): 1-7.

[3] 歐陽(yáng)永林, 馬曉明, 郭曉龍, 耿晶, 張秀平, 代春盟. 利用分頻地震屬性進(jìn)行古風(fēng)化殼巖溶儲(chǔ)層預(yù)測(cè)——以千米橋潛山凝析氣田為例[J]. 天然氣地球科學(xué), 2008, 19(3): 381-384.Ouyang Yonglin, Ma Xiaoming, Guo Xiaolong, Geng Jing, Zhang Xiuping & Dai Chunmeng. Prediction of ancient weathering crust karst reservoirs using separate frequency seismic attributes[J].Natural Gas Geoscience, 2008, 19(3): 381-384.

[4] 劉家洪, 楊平, 汪正江, 卓皆文, 杜秋定. 黔北震旦系燈影組頂部古風(fēng)化殼特征及油氣意義[J]. 中國(guó)地質(zhì), 2012, 39(4):931-938.Liu Jiahong, Yang Ping, Wang Zhengjiang, Zhuo Jiewen & Du Qiuding. Paleo-weathering crust at the top of Sinian Dengying Formation in northern Guizhou and its petroleum exploration signiベcance[J]. Geology in China, 2012, 39(4): 931-938.

[5] 趙文智, 沈安江, 潘文慶, 張寶民, 喬占峰, 鄭劍鋒. 碳酸鹽巖巖溶儲(chǔ)層類型研究及對(duì)勘探的指導(dǎo)意義——以塔里木盆地巖溶儲(chǔ)層為例[J]. 巖石學(xué)報(bào), 2013, 29(9): 3213-3222.Zhao Wenzhi, Shen Anjiang, Pan Wenqing, Zhang Baomin, Qiao Zhanfeng & Zheng Jianfeng. A research on carbonate karst reservoirs classiベcation and its implication on hydrocarbon exploration: Cases studies from Tarim Basin[J]. Acta Petrologica Sinica,2013, 29(9): 3213-3222.

[6] 司馬立強(qiáng), 黃丹, 韓世峰, 馮春珍, 焦雨佳, 寧治軍. 鄂爾多斯盆地靖邊氣田南部古風(fēng)化殼巖溶儲(chǔ)層有效性評(píng)價(jià)[J]. 天然氣工業(yè), 2015, 35(4): 7-15.Sima Liqiang, Huang Dan, Han Shifeng, Feng Chunzhen, Jiao Yujia & Ning Zhijun. Effectiveness evaluation of palaeo-weathering crust-type karst reservoirs in the southern Jingbian Gasベeld,Ordos Basin[J]. Natural Gas Industry, 2015, 35(4): 7-15.

[7] 蘇成鵬, 唐浩, 黎虹瑋, 陳虹宇, 趙東方, 萬(wàn)偉超, 等. 四川盆地東部中二疊統(tǒng)茅口組頂部鈣結(jié)殼的發(fā)現(xiàn)及其發(fā)育模式[J].古地理學(xué)報(bào), 2015, 17(2): 229-240.Su Chengpeng, Tang Hao, Li Hongwei, Chen Hongyu, Zhao Dongfang, Wan Weichao, et al. Discovery of caliches at top of the Middle Permian Maokou Formation, eastern Sichuan Basin and their developmental model[J]. Journal of Palageography, 2015,17(2): 229-240.

[8] 胡世忠. 論東吳運(yùn)動(dòng)構(gòu)造事件與二疊系分統(tǒng)界線問(wèn)題[J]. 地層學(xué)雜志, 1994, 18(4): 309-315.Hu Shizhong. On the event of Dongwu Movement and its relation with Permian subdivision[J]. Journal of Stratigraphy, 1994, 18(4):309-315.

[9] 馮少南. 東吳運(yùn)動(dòng)的新認(rèn)識(shí)[J]. 現(xiàn)代地質(zhì), 1991, 5(4): 378-384.Feng Shaonan. New knowledge on Dongwu movement[J]. Geoscience, 1991, 5(4): 378-384.

[10] 陳維濤, 周瑤琪, 馬永生, 楊恒書(shū). 關(guān)于龍門山地區(qū)東吳運(yùn)動(dòng)的存在及其性質(zhì)的認(rèn)識(shí)[J]. 地質(zhì)學(xué)報(bào), 2007, 81(11): 1518-1525.Chen Weitao, Zhou Yaoqi, Ma Yongsheng & Yang Hengshu. The knowledge on the existence and nature of the Dongwu movement in the Longmen mountain area[J]. Acta Geologica Sinica, 2007,81(11): 1518-1525.

[11] 張祖圻. 論東吳運(yùn)動(dòng)[J]. 煤田地質(zhì)與勘探, 1983, 4(3): 14-22.Zhang Zuxin. Disscussing on Dongwu movement[J]. Coal Geology & Exploration, 1983, 4(3): 14-22.

[12] 何斌, 徐義剛, 王雅玫, 肖龍. 東吳運(yùn)動(dòng)性質(zhì)的厘定及其時(shí)空演變規(guī)律[J]. 地球科學(xué)——中國(guó)地質(zhì)大學(xué)學(xué)報(bào), 2005, 30(1):89-96.He Bin, Xu Yigang, Wang Yamei & Xiao Long. Nature of the Dongwu movement and its temporal and spatial evolution[J].Earth Science—Journal of China University of Geoscience, 2005,30(1): 89-96.

[13] 李志宏, 牛志軍, 陳立德, 黃照先. 據(jù)牙形石論鄂西茅口組頂部古巖溶不整合面形成時(shí)代[J]. 地球?qū)W報(bào), 2001, 22(2): 157-158.Li Zhihong, Niu Zhijun, Chen Lide & Huang Zhaoxian. The age of paleokarst unconformity at the top of Maokou Formation in western Hubei inferred from conodonts[J]. Acta Geoscientia Sinica, 2001, 22(2): 157-158.

[14] 李儒峰, 郭彤樓, 湯良杰, 陳國(guó)飛. 海相碳酸鹽巖層系不整合量化研究及其意義——以四川盆地北部二疊系為例[J]. 地質(zhì)學(xué)報(bào), 2008, 82(3): 407-412.Li Rufeng, Guo Tonglou, Tang Liangjie & Chen Guofei. The imperfection quantitative study of the unconformity record in Sea facies carbonate rock stratum and the signiベcation: A case of the north of Sichuan Basin in Permian[J]. Acta Petrolei Sinica, 2008,82(3): 407-412.

[15] 李旭兵, 曾雄偉, 王傳尚, 劉安, 白云山. 東吳運(yùn)動(dòng)的沉積學(xué)響應(yīng)——以湘鄂西及鄰區(qū)二疊系茅口組頂部不整合面為例[J]. 地層學(xué)雜志, 2011, 35(3): 299-304.Li Xubing, Zeng Xiongwei, Wang Chuanshang, Liu An & Bai Yunshan. Sedimentary response to the Dongwu movement: A case of the unconformity on top of the Permian Maokou Formation in the western Hubei-Hunan and neighboring areas[J]. Journal of Stratigraphy, 2011, 35(3): 299-304.

[16] 江青春, 胡素云, 汪澤成, 池英柳, 楊雨, 魯衛(wèi)華, 等. 四川盆地茅口組風(fēng)化殼巖溶古地貌及勘探選區(qū)[J]. 石油學(xué)報(bào), 2012,33(6): 949-960.Jiang Qingchun, Hu Suyun, Wang Zecheng, Chi Yingliu, Yang Yu, Lu Weihua, et al. Paleokarst landform of the weathering crust of Middle Permian Maokou Formation in Sichuan Basin and selection of exploration regions[J]. Acta Petrolei Sinica, 2012,33(6): 949-960.

[17] 付曉飛, 李兆影, 盧雙舫, 付廣. 利用聲波時(shí)差資料恢復(fù)剝蝕量方法研究與應(yīng)用[J]. 大慶石油地質(zhì)與開(kāi)發(fā), 2004, 23(1):9-11.Fu Xiaofei, Li Zhaoying, Lu Shuangfang & Fu Guang. Recovering denudation thickness by interval transit-time[J]. Petroleum Geology & Oilベeld Development in Daqing, 2004, 23(1): 9-11.

[18] 袁玉松, 鄭和榮, 涂偉. 沉積盆地剝蝕量恢復(fù)方法[J]. 石油實(shí)驗(yàn)地質(zhì), 2008, 30(6): 636-642.Yuan Yusong, Zheng Herong & Tu Wei. Methods of eroded strata thickness restoration in sedimentary basins[J]. Petroleum Geology and Experiment, 2008, 30(6): 636-642.

[19] 牟中海, 唐勇, 崔炳富, 肖又軍, 王國(guó)林. 塔西南地區(qū)地層剝蝕厚度恢復(fù)研究[J]. 石油學(xué)報(bào), 2002, 23(1): 40-44.Mou Zhonghai, Tang Yong, Cui Bingfu, Xiao Youjun & Wang Guolin. Erosion thickness restoration in southwest Tarim Basin[J]. Acta Petrolei Sinica, 2002, 23(1): 40-44.

[20] 梅冥相. 碳酸鹽旋回與層序[M]. 貴陽(yáng): 貴州科技出版社,1995.Mei Mingxiang. Carbonate cycle and sequence[M]. Guiyang:Guizhou Science and Technology Publishing House, 1995.

[21] 趙宗舉, 陳軒, 潘懋, 吳興寧, 鄭興平, 潘文慶. 塔里木盆地塔中—巴楚地區(qū)上奧陶統(tǒng)良里塔格組米蘭科維奇旋回性沉積記錄研究[J]. 地質(zhì)學(xué)報(bào), 2010, 84(4): 518-536.Zhao Zongju, Chen Xuan, Pan Mao, Wu Xingning, Zheng Xingping & Pan Wenqing. Milankovitch cycles in the Upper Ordovician Lianglitage Formation in the Tazhong-Bachu Area, Tarim Basin[J]. Acta Geologica Sinica, 2010, 84(4): 518-536.

[22] 郭穎, 湯良杰, 岳勇, 陳剛, 龍禹, 謝欣彤. 旋回分析法在地層剝蝕量估算中的應(yīng)用: 以塔里木盆地玉北地區(qū)東部中下奧陶統(tǒng)鷹山組為例[J]. 中國(guó)礦業(yè)大學(xué)學(xué)報(bào), 2015, 44(4): 664-672.Guo Ying, Tang Liangjie, Yue Yong, Chen Gang, Long Yu &Xie Xintong. Application of cycle analysis method to estimate the denuded strata thickness: A case study of Middle-Lower Ordovician Yingshan Formation of the eastern Yubei area, Tarim Basin[J]. Journal of China University of Mining & Technology,2015, 44(4): 664-672.

[23] Weedon GP. Time-series analysis and cyclostratigraphy[M]. London: Cambridge University Press, 2003.

[24] Zhang YB, Jia CZ, Zhao ZJ, Jiang ZX & Xu J. Characteristics of Milankovitch Cycles in the Mid-Permian Liangshan and Qixia Formations of the Sichuan Basin—Examples from Well-Long17 and Well-Wujia1[J]. Acta Geologica Sinica, 2012, 86(6): 1045-1059.

[25] 張運(yùn)波, 趙宗舉, 袁圣強(qiáng), 鄭民. 頻譜分析法在識(shí)別米蘭科維奇旋回及高頻層序中的應(yīng)用——以塔里木盆地塔中—巴楚地區(qū)下奧陶統(tǒng)鷹山組為例[J]. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版),2011, 41(2): 400-410.Zhang Yunbo, Zhao Zongju, Yuan Shengqiang & Zheng Min. Application of spectral analysis to identify Milankovitch cycles and high-frequency sequences—Take the Lower Ordovician Yingshan Formation of mid-Tarim Basin as an example[J]. Journal of Jilin Univcrsity (Earth Science Edition), 2011, 41(2): 400-410.

[26] Berger A & Loutre MF. Astronomical forcing through geological time[M]//De Boer PL, Smith DG. Orbital Forcing and Cyclic Sequences. Oxford: Blackwell Scientiベc Publications, 1994: 15-24.

[27] Berger A, Loutre MF & Laskar J. Stability of the astronomical frequencies over the earth's history for paleoclimate studies[J].Science, 1992, 255(5044): 560-566.

[28] 李四光. 中國(guó)東南部古生代后期之造山運(yùn)動(dòng)[J]. 中國(guó)地質(zhì)學(xué)會(huì)志, 1931, 11(2): 209-217.Lee JS. Variskian or Hercynian movement in south-eastern China[J]. Bulletin of the Geological Society of China, 1931, 11(2):209-217.

[29] 陳顯群, 劉應(yīng)楷, 童鵬. 東吳運(yùn)動(dòng)質(zhì)疑及川黔運(yùn)動(dòng)之新見(jiàn)[J].石油與天然氣地質(zhì), 1987, 8(4): 412-423.Chen Xianqun, Liu Yingkai & Tong Peng. Query about Dongwu movement and new opinion on Chuanqian movement[J]. Oil &Gas Geology, 1987, 8(4): 412-423.

[30] 汪澤成, 趙文智, 胡素云, 徐安娜, 江青春, 姜華, 等. 克拉通盆地構(gòu)造分異對(duì)大油氣田形成的控制作用——以四川盆地震旦系—三疊系為例[J]. 天然氣工業(yè), 2017, 37(1): 9-23.Wang Zecheng, Zhao Wenzhi, Hu Suyun, Xu Anna, Jiang Qingchun, Jiang Hua, et al. Control of tectonic differentiation on the formation of large oil and gas ベelds in craton basins: A case study of Sinian–Triassic of the Sichuan Basin[J]. Natural Gas Industry,2017, 37(1): 9-23.

[31] Haq BU & Schutter SR. A chronology of Paleozoic sea-level changes[J]. Science, 2008, 64(5898): 64-68.

[32] 王洪浩, 李江海, 李維波, 周肖貝, 馬麗亞, 劉仲蘭. 岡瓦納大陸古生代冰蓋分布研究[J]. 中國(guó)地質(zhì), 2014, 41(6): 2132-2143.Wang Honghao, Li Jianghai, Li Weibo, Zhou Xiaobei, Ma Liya &Liu Zhonglan. A study of the distribution of ice sheets on Gondwana Supercontinent in Paleozoic[J]. Geology in China, 2014,41(6): 2132-2143.

[33] Van Den Heuvel EPJ & Buurman P. Possible causes of glaciations[M]//Herman Y. Marine geology and oceanography of the Arctic Seas. Berlin Heidelberg: Springer, 1974.

[34] Smith LB Jr & Read JF. Rapid onset of Late Paleozoic glaciation on Gondwana: Evidence from Upper Mississippian strata of the Midcontinent, United States[J]. Geology, 2000, 28(3): 279-282.

[35] Smith LB & Read JF. Discrimination of local and global effects on Upper Mississippian stratigraphy, Illinois Basin, USA[J].Journal of Sedimentary Research, 2001, 71(6): 985-1002.

[36] Golonka J & Ford D. Pangean (Late Carboniferous-Middle Jurassic) paleoenvironment and lithofacies[J]. Palaeogeography,Palaeoclimatology, Palaeoecology, 2000, 161(1/2): 1-34.

[37] 李江海, 楊靜懿, 馬麗亞, 王洪浩. 顯生宙烴源巖分布的古板塊再造研究[J]. 中國(guó)地質(zhì), 2013, 40(6): 1683-1698.Li Jianghai, Yang Jingyi, Ma Liya & Wang Honghao. A study of the distribution of source rocks in Phanerozoic based on paleoplate reconstruction[J]. Geology in China, 2013, 40(6): 1683-1698.

[38] 王成善, 李祥輝, 陳洪德, 覃建雄. 中國(guó)南方二疊紀(jì)海平面變化及升降事件[J]. 沉積學(xué)報(bào), 1999, 17(4): 536-541.Wang Chengshan, Li Xianghui, Chen Hongde & Qin Jianxiong.Permian sea-level changes and rising-falling events in South China[J]. Acta Sedimentologica Sinica, 1999, 17(4): 536-541.

[39] Veizer J, Ala D, Azmy K, Bruckschen P, Buhl D, Bruhn F, et al.87Sr/86Sr, δ13C and δ18O evolution of Phanerozoic seawater[J].Chemical Geology, 1999, 161(1/3): 59-88.

[40] 章明圓. 湖南桑植仁村坪二疊紀(jì)瓜德魯普統(tǒng)—樂(lè)平統(tǒng)生物地層和碳同位素變化[D].北京: 中國(guó)科學(xué)院大學(xué), 2013.Zhang Mingyuan. The Guandalupian-Lopingian biostratigraphy and C-isotopic excursions in Rencunping Section of Sangzhi,Hunan[D]. Beijing: University of Chinese Academy of Sciences,2013.