成 俊 ，李少榮 ，陳朝剛 ，張健強 ，宋科雄 ，鄭永華

頁巖氣作為一種非常規(guī)天然氣，在涪陵礁石壩和長寧–威遠地區(qū)均已取得巨大產(chǎn)能。在地質演化過程中，頁巖氣含量受頁巖的有機質含量、熱演化程度、儲層的溫度和壓力影響[1–2]。含油氣系統(tǒng)模擬技術基于盆地動力學，通過開展構造演化史、埋藏史、熱演化史、生烴史等模擬，分析油氣藏的保存條件[3–4]，近年來在頁巖氣動態(tài)研究和資源評價中有較好的運用[5–10]。

渝東南地區(qū)位于川中隆起與黔中隆起之間的凹陷褶皺帶，屬于上揚子準地臺黔南川北地層分區(qū)，下古生界發(fā)育下寒武統(tǒng)牛蹄塘組和上奧陶統(tǒng)五峰組–下志留統(tǒng)龍馬溪組兩套頁巖氣層，儲層靜態(tài)指標較好，但地質構造復雜，頁巖氣勘探程度較低，頁巖氣資源勘探也未取得突破。本文運用PetroMod軟件對渝東南地區(qū)下古生界頁巖氣層開展含油氣系統(tǒng)模擬，分析成藏主控因素，為頁巖氣的勘探潛力評價提供依據(jù)。

研究區(qū)位于渝東南秀山穹褶帶，面積1 002 km2，發(fā)育除志留系上統(tǒng)、泥盆系下統(tǒng)和石炭系之外的震旦系到白堊系地層（圖1）。該區(qū)已開展不同程度的地質調查、二維地震勘探和頁巖氣井鉆探工作，已完成2口地質井（YD2、YC6）和2口預探井（YY1、YY2），其中YD2、YC6和YY2井針對五峰組–龍馬溪組，YY1井針對牛蹄塘組兼探五峰–龍馬溪組，該區(qū)頁巖氣勘探取得了一定的成果資料（表1）。

圖1 研究區(qū)構造地質簡圖

表1 研究區(qū)目的層頁巖氣地質特征

1 模擬方法

1.1 埋藏史模擬

埋藏史模擬是含油氣系統(tǒng)模擬的基礎，古地層模型受控于地層的沉積和剝蝕，基于現(xiàn)今的地層沉積序列，根據(jù)不同地層孔隙度演化去除壓力效應，建立不同地史時期古構造模型[11–12]。研究區(qū)地層屬于正常壓實帶沉積，宜采用回剝法模擬地層埋藏史?；貏兎ㄊ歉鶕?jù)孔隙度–深度的關系來恢復古厚度，在正常壓實情況下，二者服從指數(shù)分布[13–14]：

式中， (!)是深度為z時的孔隙度，%$是地表孔隙度，*$是壓實系數(shù)，#$是地層各巖性含量。

研究區(qū)地層總體分為砂巖、泥巖、灰?guī)r和白云巖四種巖性，通過實測孔隙度校正各巖性的壓實系數(shù)，結合現(xiàn)今的地層厚度，恢復古地層沉積厚度。渝東南地區(qū)震旦紀之后受華力西期、燕山期和喜馬拉雅期構造抬升運動的影響，分別在泥盆紀–石炭紀、二疊紀、白堊紀–第四紀發(fā)生三期剝蝕事件，最后一期剝蝕厚度較大，剝蝕厚度為1 700～4 000 m[15]。

1.2 熱演化史模擬

熱演化史模擬采用地球熱力學與地球化學相結合，通過確定巖石熱參數(shù)、大地熱流、古地表溫度和古水深等參數(shù)，模擬地溫史和成熟度演化，并通過實測的鏡質體反射率Ro校正[16–17]。結合四川盆地砂巖、泥巖、碳酸鹽巖的熱導率研究成果[18]，根據(jù)研究區(qū)地層巖性含量，確定各地層的巖石熱參數(shù)；根據(jù)相對海平面變化曲線，結合沉積相研究成果，恢復了各時期水深變化曲線；根據(jù)前人研究成果，建立了550 Ma前至今的古地表溫度和熱流演化曲線。

1.3 生烴史模擬

根據(jù)烴源巖地化和吸附參數(shù)，選用合適的生烴動力學模型，對生烴時間、油氣轉化率、生烴量等進行生烴史模擬。油氣轉化率為干酪根轉化為油氣的比例；針對頁巖氣，還需要模擬吸附氣和游離氣含量；烴源巖參數(shù)主要包括有機碳（TOC）和氫指數(shù)（IH）。為確保模擬結果的準確性，對目的層始TOC和IH進行恢復，選擇原始IH為500 mg/g，恢復牛蹄塘組和五峰組–龍馬溪組頁巖TOC分別為8.74%和3.68%；吸附參數(shù)根據(jù)Langmuir等溫吸附實驗確定；根據(jù)研究區(qū)目的層和烴源巖特征，生烴動力學模型選用PetroMod軟件中與之類似的巴黎盆地土阿辛階（對應于下侏羅統(tǒng)）頁巖的生烴動力學模型。

2 模擬結果

2.1 埋藏史模擬結果

早寒武紀–早志留紀時期，隨著水體加深，研究區(qū)由淺海碳酸鹽巖沉積轉變?yōu)樯钏鄮r、頁巖以及硅質頁巖等沉積為主；中志留紀后過渡為緩慢沉降，以碎屑巖沉積為主，泥盆紀晚期地層整體輕微抬升剝蝕；石炭紀–二疊紀末期，早期近海湖沼演變?yōu)楹Ｋ笠?guī)模進泛，晚期處于海陸交互環(huán)境，沉積淺海碳酸鹽巖建造，早二疊紀末期有輕微剝蝕；三疊紀海侵加大，主要沉積蒸發(fā)碳酸鹽巖建造；侏羅紀–白堊紀早期，研究區(qū)逐漸演變?yōu)榇笮蛢汝懓枷莺?，發(fā)育河流湖泊相沉積；白堊紀中期至今，燕山運動和喜馬拉雅運動致研究區(qū)大面積抬升剝蝕，幾乎無古近系、新近系沉積，發(fā)育大量逆沖斷層，形成現(xiàn)今的北東–北北東向褶皺構造形態(tài)。

根據(jù)埋藏史模擬結果，研究區(qū)經(jīng)歷了寒武紀、志留紀和侏羅紀三次快速沉積時期，沉積速率分別為51 m/Ma、39～76 m/Ma、36 m/Ma；并經(jīng)歷了泥盆紀–石炭紀（400～285 Ma）、二疊紀（270～260 Ma）、白堊紀–第四紀（90～0 Ma）三次剝蝕時期，前兩次剝蝕厚度較小（<100 m），最后一期累積剝蝕厚度較大（YY1井剝蝕厚度3 499 m、YY2井剝蝕厚度2 144 m）（圖 2）。

2.2 熱演化史模擬結果

YY1井牛蹄塘組烴源巖約在 495 Ma前進入生油窗，411 Ma前進入早期生氣階段，275 Ma前進入晚期生氣階段，目前已不處于生氣階段；五峰組–龍馬溪組烴源巖約在393 Ma前進入生油窗，132 Ma前進入早期生氣階段，目前處于晚期生氣階段。YY2井五峰組–龍馬溪組烴源巖約在423 Ma前進入生油窗，284 Ma前進入早期生氣階段，164 Ma前進入晚期生氣階段，目前處于晚期生氣階段。YY2井比YY1井的五峰組–龍馬溪組演化程度更高（圖2）。

圖2 單井熱演化模擬結果

2.3 生烴史模擬結果

YY1井牛蹄塘組烴源巖約在521.5 Ma前開始生烴，至429.8 Ma前轉化率達到50%，到238.6 Ma前接近100%；牛蹄塘組最佳生烴期在晚志留世，三疊紀已停止生烴。由于烴源巖演化成熟度較高，隨著沉積壓實作用，孔隙度降低，且后期強烈的構造運動導致現(xiàn)今保存下來的含氣量較低。YY2井五峰組–龍馬溪組烴源巖約在427 Ma前開始生烴，414 Ma前轉化率達到50%，目前轉化率為99%；五峰組–龍馬溪組從志留紀晚期開始生烴，三疊紀末含氣量達到最大，目前在構造高部位還有繼續(xù)生烴的能力，含氣量以吸附氣為主（吸附氣與游離氣比值約為5.3:1）（圖 3）。

圖3 單井生烴史模擬結果

3 頁巖氣成藏過程分析

現(xiàn)今頁巖氣含量與歷史事件中生烴演化程度和構造運動有關，生烴演化為頁巖氣提供物質基礎，構造運動則影響頁巖氣的保存；四川盆地內與周緣地區(qū)保存條件的差異，主要是由于燕山運動中期所控制的構造變形和沉降、頁巖生排烴歷史、抬升剝蝕作用等差異決定。

劇烈的構造運動對于烴類保存條件的破壞，同構造運動本身一樣也是幕式的，如果在生烴階段有較長時間的構造穩(wěn)定期，則隨著斷裂、裂縫等疏導體系的愈合，烴類又能重新聚集，新生的烴類對侵入烴源巖的水和空氣（氧氣難以在還原條件下長期存在，絕大部分為氮氣）可以產(chǎn)生反置換作用，逐漸排擠出水和氮氣，重新成藏。只要生烴作用仍在繼續(xù)，這種“成藏—構造破壞—再成藏”的過程就會一直延續(xù)。五峰組–龍馬溪組烴源巖就是這種情況，至今仍有一定勘探潛力；而牛蹄塘組烴源巖在燕山期強烈構造運動之前就已經(jīng)停止生烴，構造運動造成的大型斷裂發(fā)育并切穿頁巖層系，牛蹄塘組頁巖氣被破壞后無法重新成藏，并與氮氣發(fā)生置換作用，因此出現(xiàn)了目前氮氣成藏的現(xiàn)象（圖4）。根據(jù)渝黔湘鄂地區(qū)大量鉆井顯示，牛蹄塘組頁巖儲層中氮氣含量普遍很高。

4 結論

（1）本文介紹的含油氣系統(tǒng)模擬技術為頁巖氣成藏的熱演化、生烴和保存的時空關系研究提供了有益的依據(jù)。

（2）渝東南地區(qū)經(jīng)歷了寒武紀、志留紀和侏羅紀三次快速沉積時期，分別沉積了淺海碳酸鹽巖、深水泥頁巖和河流湖泊相沉積；在泥盆紀、早二疊紀和白堊紀–古近紀分別經(jīng)歷了三次剝蝕時期，最后一期剝蝕作用強烈。

（3）現(xiàn)今頁巖氣含量與歷史事件中生烴演化程度和構造運動有關，渝東南五峰組–龍馬溪組從志留紀晚期開始生烴，三疊紀末含氣量達到最大，在頁巖氣“成藏—構造破壞—再成藏”的過程中仍在持續(xù)生烴，至今仍具一定勘探潛力；牛蹄塘組最佳生烴期在晚志留世，燕山期構造運動之前已經(jīng)停止生烴，頁巖氣藏發(fā)生置換作用，導致現(xiàn)今含氣量極低且以氮氣為主。

圖4 渝東南地區(qū)下古生界含油氣系統(tǒng)事件

參考文獻

[1] 畢赫，姜振學，李鵬.頁巖含氣量主控因素及定量預測方法[J].大慶石油地質與開發(fā)，2014，33（1）：160–164.

[2] 張漢榮.川東南地區(qū)志留系頁巖含氣量特征及其影響因素[J].天然氣工業(yè)，2016，36（8）：36–42.

[3] 石廣仁.油氣盆地數(shù)值模擬方法[M].北京：石油工業(yè)出版社，1999：14–39.

[4] 邢衛(wèi)新，費永濤，楊永毅，等.含油氣盆地成藏期分析方法及進展[J].石油地質與工程，2006，20（6）：12–15.

[5] 陶樹，湯達禎，許浩，等.中、上揚子區(qū)寒武–志留系高過成熟烴源巖熱演化史分析[J].自然科學進展，2009，19（10）：1 126–1 133.

[6] 朱傳慶，饒松，袁玉松，等.川東南地區(qū)古生界主要頁巖層系熱演化[J].煤炭學報，2013，38（5）：834–839.

[7] 曹環(huán)宇，朱傳慶，邱楠生.川東地區(qū)下志留統(tǒng)龍馬溪組熱演化[J].地球科學與環(huán)境學報，2015，37（6）：22–32.

[8] 江強，朱傳慶，邱楠生，等.川南地區(qū)熱史及下寒武統(tǒng)筇竹寺組頁巖熱演化特征[J].天然氣地球科學，2015，26（8）：1 563–1 570.

[9] XIULIAN YAO，YUNPENG WANG.Assessing shale gas resources of Wufeng–Longmaxi shale （O–S） in Jiaoshiba area， SE Sichuan（China）[J].Petroleum Science& Technology，2016，34（11–12）：1 008–1 015.

[10] 徐二社，李志明，楊振恒.彭水地區(qū)五峰–龍馬溪組頁巖熱演化史及生烴史研究——以 PY1井為例[J].石油實驗地質，2015，37（4）：494–499.

[11] OPERA A，ALIZADEH B，SARAFDOKHT H，et al.Burial history reconstruction and thermal maturity modeling for the Middle Cretaceous–Early Miocene Petroleum System，southern Dezful Embayment，SW Iran[J].International Journal of Coal Geology，2013，120（6）：1–14.

[12] UNDERDOWN R，REDFERN J.Petroleum generation and migration in the Ghadames Basin，North Africa：A two–dimensional basin–modeling study[J].AAPG Bulletin，2008，92（1）：53–76.

[13] VAN SICKEL W A，KOMINZ M A，MILLER K G，et al.Late Cretaceous and Cenozoic sea-level estimates：backstripping analysis of borehole data，onshore New Jersey[J].Basin Research，2015，16（4）：451–465.

[14] 倪斌，湯良杰，郭穎，等.塔里木盆地玉北地區(qū)埋藏史及熱史分析[J].現(xiàn)代地質，2017，31（2）：357–366.

[15] 曾祥亮，劉樹根，黃文明，等.四川盆地志留系龍馬溪組頁巖與美國Fort Worth盆地石炭系Barnett組頁巖地質特征對比[J].地質通報，2011，30（2～3）：372–384.

[16] GALUSHKIN Y I，SITAR K A，KUNITSYNA A V.Realization of the oil–generation potential in the estimated source rocks of the sedimentary basin in the shelf of northeastern Sakhalin：Numerical simulation[J].Geochemistry International，2011，49（1）：1–12.

[17] GANJAVAR K K，JOHAN K.The thermal evolution of solid bitumens，bitumen reflectance and kinetic modeling of reflectance：application in petroleum and ore prospecting[J].Energy Sources，1993，15（2）：181–204.

[18] 徐明，朱傳慶，田云濤，等.四川盆地鉆孔溫度測量及現(xiàn)今地熱特征[J].地球物理學報，2011，54（4）：1 052–1 060.