高棟臣，郭 超，姜呈馥，張麗霞，王 暉，史 鵬，陳奕奕

(陜西延長(zhǎng)石油(集團(tuán))有限責(zé)任公司研究院，西安 710075)

生烴熱模擬實(shí)驗(yàn)是認(rèn)識(shí)烴源巖的重要手段，可以再現(xiàn)地質(zhì)體中有機(jī)質(zhì)熱解演化過(guò)程，為認(rèn)識(shí)生烴過(guò)程與機(jī)理、建立成烴模式、評(píng)價(jià)盆地成烴潛力提供依據(jù)[1-2]。熱模擬最早可追溯到18 世紀(jì)末期的魚(yú)油加熱實(shí)驗(yàn)。隨后的發(fā)展過(guò)程中，實(shí)驗(yàn)物質(zhì)逐漸豐富，裝置不斷改進(jìn)。實(shí)驗(yàn)物質(zhì)從簡(jiǎn)單的純有機(jī)化合物到煤，到干酪根、全巖、抽提物、瀝青和原油等各類(lèi)物質(zhì)，實(shí)驗(yàn)裝置發(fā)展出開(kāi)放體系、半開(kāi)放體系和封閉體系等，實(shí)驗(yàn)方法從最初的無(wú)水熱模擬過(guò)渡到加水、加催化劑熱模擬[3-11]。石油地質(zhì)工作者在生烴熱模擬實(shí)驗(yàn)方面進(jìn)行了諸多嘗試[12-27]，取得了豐富的成果，推動(dòng)了實(shí)驗(yàn)技術(shù)的進(jìn)步，使實(shí)驗(yàn)逐漸貼近地質(zhì)條件下的生烴過(guò)程。

鄂爾多斯盆地上古生界地層天然氣資源豐富，是國(guó)內(nèi)天然氣勘探開(kāi)發(fā)的熱點(diǎn)。上古生界烴源巖的認(rèn)識(shí)和評(píng)價(jià)對(duì)于掌握天然氣及頁(yè)巖氣的成藏規(guī)律、有利區(qū)分布和資源儲(chǔ)量等意義重大。然而，該地區(qū)上古生界二疊系山西組低成熟度頁(yè)巖成烴演化研究成果尚少。本次研究選取了上古生界山西組低成熟度頁(yè)巖進(jìn)行生烴熱模擬實(shí)驗(yàn)，旨在揭示山西組烴源巖的生烴過(guò)程，了解生烴產(chǎn)物的變化規(guī)律及影響因素，建立古生界烴源巖生烴模式，為認(rèn)識(shí)山西組頁(yè)巖的生烴潛力及天然氣、頁(yè)巖氣等資源評(píng)價(jià)提供參考。

1 樣品與實(shí)驗(yàn)方法

1.1 樣品準(zhǔn)備

頁(yè)巖全巖樣品代表沉積盆地?zé)N源巖自然組成，較干酪根更貼近地質(zhì)真實(shí)情況，一般選取低成熟度、高有機(jī)碳含量、有機(jī)質(zhì)類(lèi)型典型的樣品進(jìn)行實(shí)驗(yàn)。對(duì)低成熟度頁(yè)巖樣品進(jìn)行封閉體系下的加水高壓熱模擬實(shí)驗(yàn)，可以明確烴源巖的生烴、演化過(guò)程及熱解產(chǎn)物的變化規(guī)律。

鄂爾多斯盆地延長(zhǎng)地區(qū)山西組頁(yè)巖埋深基本超過(guò)3 km，有機(jī)質(zhì)成熟度高，有機(jī)質(zhì)類(lèi)型為Ⅲ型，因此，選取盆地內(nèi)地質(zhì)條件相近的其他地區(qū)低成熟度樣品，開(kāi)展生烴實(shí)驗(yàn)來(lái)模擬整個(gè)成烴演化過(guò)程。在盆地東緣的山西保德縣扒樓溝剖面山西組露頭中采集了新鮮碳質(zhì)頁(yè)巖樣品，開(kāi)展封閉體系下加水高壓熱模擬實(shí)驗(yàn)。實(shí)驗(yàn)樣品屬于濱淺海附近的沉積產(chǎn)物，沉積相類(lèi)型和有機(jī)質(zhì)類(lèi)型與油氣區(qū)內(nèi)部相同，熱演化程度較低，鏡質(zhì)體反射率Ro≈0.4%，適合于生烴熱模擬實(shí)驗(yàn)研究。熱模擬實(shí)驗(yàn)樣品的地球化學(xué)分析測(cè)試表明，山西組頁(yè)巖樣品熱演化程度較低，生烴潛力較好，有機(jī)質(zhì)類(lèi)型為Ⅲ型(表1)。

1.2 實(shí)驗(yàn)方法

頁(yè)巖生烴熱模擬實(shí)驗(yàn)由中科院廣州地化所有機(jī)地化重點(diǎn)實(shí)驗(yàn)室完成，實(shí)驗(yàn)過(guò)程參照參考文獻(xiàn)[28-29]。實(shí)驗(yàn)采用密封金管—高壓釜加水生烴熱模擬實(shí)驗(yàn)裝置，頁(yè)巖樣品用量為20～50 mg，實(shí)驗(yàn)壓力為50 MPa，分別以20 ℃/h和2 ℃/h升溫速率恒壓加熱。利用GC6890氣相色譜儀進(jìn)行氣體產(chǎn)物分析，用GC-IRMS同位素質(zhì)譜儀進(jìn)行氣體碳同位素測(cè)定。

根據(jù)EASY%Ro原理，一定升溫速率下的鏡質(zhì)體反射率與熱模擬溫度存在對(duì)應(yīng)關(guān)系[30]，從而建立有機(jī)質(zhì)成熟度與模擬產(chǎn)物成分以及同位素特征的演化關(guān)系，明確有機(jī)質(zhì)生烴演化過(guò)程。

2 結(jié)果與討論

2.1 氣態(tài)產(chǎn)物演化特征

2.1.1 升溫速率對(duì)氣體產(chǎn)率的影響

生烴熱模擬實(shí)驗(yàn)分別在20 ℃/h和2 ℃/h 2種升溫速率條件下進(jìn)行。實(shí)驗(yàn)結(jié)果顯示，實(shí)驗(yàn)樣品在20 ℃/h條件下生成烴類(lèi)氣體及非烴氣體產(chǎn)率比在2 ℃/h條件下低，說(shuō)明熱模擬過(guò)程中高升溫速率不利于氣體生成(圖1)。

2.1.2 溫度對(duì)氣體產(chǎn)率的影響

實(shí)驗(yàn)顯示，隨著溫度增加，C1、C1-5、CO2、H2S、非烴和總熱解烴產(chǎn)率均呈增大趨勢(shì)。C1-5/總熱解烴隨溫度增加而增加，其變化速率呈現(xiàn)出緩慢增加到快速增加到緩慢增加的趨勢(shì)。烴/非烴隨溫度變化規(guī)律與C1-5/總熱解烴表現(xiàn)出相似特征(圖1)。

2.1.3 熱解烴氣體特征

隨溫度增加，實(shí)驗(yàn)樣品的C1/C1-5比值曲線呈現(xiàn)“S”型特征。認(rèn)為該趨勢(shì)的形成與烴類(lèi)二次裂解有關(guān)。低溫階段隨著溫度增加，熱演化程度升高，多種氣體生成，甲烷所占比例逐漸下降，達(dá)到二次裂解臨界點(diǎn)時(shí)，甲烷所占比例達(dá)到低點(diǎn)。此后隨溫度增加，烴類(lèi)二次裂解形成大量甲烷，甲烷比例快速增加；此外液態(tài)烴產(chǎn)率在高溫階段存在下降曲線也佐證了二次裂解的存在(圖1，2)。

2.1.4 熱解非烴氣體特征

隨溫度升高，CO2、H2S和H2等氣體產(chǎn)率逐漸增大，烴類(lèi)/非烴氣體不斷增大，表明隨溫度增加，烴類(lèi)氣體生成速率高于非烴氣體。隨著溫度升高，CO2/C1-5比值快速降低，從247.64減小到2.32(20 ℃/h速率升溫)，說(shuō)明烴類(lèi)氣體生成速率隨溫度增加明顯高于CO2(圖1)。研究認(rèn)為該現(xiàn)象的形成與CO2活化能較低有關(guān)，以致CO2在先期快速生成，后期生成量較少，造成比值不斷降低[31]。

表1 鄂爾多斯盆地山西組頁(yè)巖實(shí)驗(yàn)樣品地球化學(xué)指標(biāo)Table 1 Geochemical properties of shale sample from Shanxi Formation, Ordos Basin

圖1 鄂爾多斯盆地山西組頁(yè)巖樣品熱模擬氣態(tài)產(chǎn)物產(chǎn)率隨溫度變化特征Fig.1 Characteristics of gas product yield with experimental temperature changes of shale samples from Shanxi Formation, Ordos Basin

圖2 鄂爾多斯盆地山西組頁(yè)巖樣品生烴熱模擬液態(tài)產(chǎn)物產(chǎn)率隨溫度變化特征Fig.2 Characteristics of liquid product yield with experimental temperature changes of shale samples from Shanxi Formation, Ordos Basin

2.2 液態(tài)產(chǎn)物演化特征

頁(yè)巖生烴熱模擬液態(tài)產(chǎn)物隨溫度變化特征如圖2所示，隨熱解溫度升高，總熱解油的產(chǎn)率總體呈現(xiàn)出先增加后降低的特征。在不同升溫速率的實(shí)驗(yàn)中，液態(tài)烴最大產(chǎn)率時(shí)的溫度不同。

以20 ℃/h速率升溫時(shí)，樣品總熱解油產(chǎn)率在7.29～57.31 mg/g之間。在實(shí)驗(yàn)最低溫度335.9 ℃(對(duì)應(yīng)的鏡質(zhì)體反射率Ro=0.57%)時(shí)，山西組頁(yè)巖樣品總熱解油產(chǎn)率為16.46 mg/g，隨溫度升高總熱解油產(chǎn)率增大，在432 ℃時(shí)(對(duì)應(yīng)的Ro=1.20%)總熱解油產(chǎn)率達(dá)到峰值57.31 mg/g，此后產(chǎn)率隨溫度升高反而降低，產(chǎn)率隨溫度變化曲線呈單峰。以2 ℃/h速率升溫時(shí)，在實(shí)驗(yàn)最低溫度337.2 ℃(對(duì)應(yīng)的Ro=0.74%)時(shí)，總熱解油產(chǎn)率為36.10 mg/g，隨溫度升高總熱解油產(chǎn)率增大，在385.3 ℃時(shí)(對(duì)應(yīng)的Ro=1.08%)總熱解油產(chǎn)率達(dá)到峰值51.73 mg/g，隨后產(chǎn)率隨溫度升高而降低，產(chǎn)率隨溫度變化曲線也呈單峰。

不同升溫速率條件下，總熱解油產(chǎn)率曲線變化趨勢(shì)相同。2 ℃/h與20 ℃/h升溫速率條件下產(chǎn)率曲線存在交叉點(diǎn)，交叉點(diǎn)的意義是兩者在某一相同溫度處產(chǎn)率相同。從實(shí)驗(yàn)開(kāi)始到2 ℃/h的產(chǎn)率曲線達(dá)到峰值之前，2 ℃/h條件下產(chǎn)率高于20 ℃/h條件下產(chǎn)率；此后2 ℃/h條件下產(chǎn)率逐漸下降，20 ℃/h條件下產(chǎn)率逐漸上升，直至交叉點(diǎn)處兩者產(chǎn)率相同。交叉點(diǎn)后20 ℃/h條件下產(chǎn)率始終高于2 ℃/h條件下產(chǎn)率。

C6-14和C14+兩者產(chǎn)率曲線與總熱解油表現(xiàn)出相似特征。與總熱解油相比，不同升溫速率下輕質(zhì)油(C6-14)和重質(zhì)油(C14+)最大產(chǎn)率形成的溫度與總熱解油最大產(chǎn)率形成的溫度稍有不同(圖2)。升溫速率對(duì)液態(tài)烴最大產(chǎn)率(峰值)有明顯影響，無(wú)論輕烴(C6-14)、重?zé)N(C14+)還是總熱解油，低升溫速率(2 ℃/h)條件下其產(chǎn)率峰值形成的溫度低于高的升溫速率(20 ℃/h)條件下產(chǎn)率峰值形成的溫度(圖2)，二者相差60～80 ℃。

20 ℃/h升溫速率時(shí)最大總熱解油產(chǎn)率為57.31 mg/g，最大重質(zhì)油(C14+)產(chǎn)率為41.90 mg/g，最大輕質(zhì)油(C6-14)產(chǎn)率為24.47 mg/g；在2 ℃/h升溫速率時(shí)最大總熱解油產(chǎn)率為51.73 mg/g，最大重質(zhì)油(C14+)產(chǎn)率為37.51 mg/g，最大輕質(zhì)油(C6-14)產(chǎn)率為24.77 mg/g。

2.3 碳同位素組成

頁(yè)巖樣品生烴熱解氣體產(chǎn)物碳同位素特征分析表明，δ13C1值隨溫度升高先變輕后變重，而δ13C2和δ13C3總體上隨溫度升高而變重，δ13CCO2隨溫度變化不明顯(圖3)。不同升溫速率下，δ13C1、δ13C2、δ13C3和δ13CCO2變化趨勢(shì)基本相同，同一溫度下δ13C值相近。

樣品在熱解溫度低于384.5 ℃(Ro<0.8%)時(shí)，氣態(tài)烴的δ13C呈現(xiàn)局部倒轉(zhuǎn)現(xiàn)象，即δ13C1<δ13C2>δ13C3。通常認(rèn)為有機(jī)成因原生烷烴氣的碳同位素系列是正碳同位素系列(δ13C1<δ13C2<δ13C3)，造成烴類(lèi)氣體中碳同位素倒轉(zhuǎn)現(xiàn)象可能由于有機(jī)成因和無(wú)機(jī)成因氣混合、煤成氣和油型氣混合、同型不同源或同源不同期氣體的混合以及細(xì)菌的氧化作用等幾種原因。

圖3 鄂爾多斯盆地山西組頁(yè)巖熱解氣態(tài)產(chǎn)物δ13C隨熱演化程度變化特征Fig.3 δ13C vs. Ro of gaseous product from pyrolysis simulation of shale samples from Shanxi Formation, Ordos Basin

氣態(tài)烴的同位素分餾效應(yīng)由12C-12C和12C-13C化學(xué)鍵能所決定。在熱解初期碳同位素較輕是由2個(gè)原因造成的：(1)12C-12C鍵能弱，易斷裂，使得初次裂解產(chǎn)物12C更易富集，造成前期生成的氣態(tài)烴碳同位素變輕，而后期生成的氣態(tài)烴及其二次裂解產(chǎn)物碳同位素值變重[32]。(2)隨著濕度的增加，12C可能更易富集，致使烴類(lèi)氣體碳同位素變輕[33-34]。液態(tài)烴二次裂解時(shí)的濕度比生成液態(tài)烴的濕度大，更易形成較低δ13C值的現(xiàn)象。熱模擬過(guò)程中生油和油裂解階段，由于濕度的不斷增加，δ13C值逐漸減小，與殘余有機(jī)質(zhì)熱解形成的重δ13C烴類(lèi)氣體混合，形成碳同位素倒轉(zhuǎn)系列。綜合以上分析認(rèn)為，該樣品碳同位素局部倒轉(zhuǎn)的現(xiàn)象是由δ13C3變輕引起的，即熱模擬實(shí)驗(yàn)中，來(lái)自有機(jī)質(zhì)早期熱演化階段的丙烷和液態(tài)烴二次裂解的丙烷相互混合，導(dǎo)致碳同位素變輕。

2.4 生烴演化模式

通過(guò)開(kāi)展頁(yè)巖樣品封閉體系下加水高壓熱模擬實(shí)驗(yàn)，統(tǒng)計(jì)分析隨溫度升高多種烴類(lèi)產(chǎn)物產(chǎn)率的變化特征，將頁(yè)巖熱演化過(guò)程劃分為3個(gè)階段。

2.4.1 熱催化生油氣階段

該階段Ro<1.20%，熱模擬溫度不高于440 ℃，屬于有機(jī)質(zhì)熱演化的初始階段。此階段烴類(lèi)產(chǎn)物以液態(tài)烴占絕對(duì)優(yōu)勢(shì)，氣態(tài)烴產(chǎn)物生成較少。Ro=1.20%時(shí)達(dá)到生油高峰。在達(dá)到峰值之前，氣態(tài)烴產(chǎn)率較低，總烴、氣態(tài)烴和液態(tài)烴等產(chǎn)物隨溫度增加，產(chǎn)率持續(xù)增加，峰值之后液態(tài)烴逐漸降低，總烴和氣態(tài)烴不斷增加。

2.4.2 熱裂解生濕氣階段

該階段Ro=1.20%～3.50%，熱模擬實(shí)驗(yàn)溫度在440～570 ℃。本階段典型特征是實(shí)驗(yàn)產(chǎn)物中氣態(tài)烴大量生成，液態(tài)烴產(chǎn)物因二次裂解不斷消耗，其產(chǎn)率持續(xù)降低，氣態(tài)烴與液態(tài)烴的產(chǎn)率比顯著增大。

圖4 鄂爾多斯盆地山西組頁(yè)巖生烴模式Fig.4 Hydrocarbon generation model of shale from Shanxi Formation, Ordos Basin

以20 ℃/h速率升溫，當(dāng)熱解溫度達(dá)到575.7 ℃時(shí)，氣態(tài)烴產(chǎn)率達(dá)到89.15 mg/g，而液態(tài)烴產(chǎn)率降為10.23 mg/g，氣態(tài)烴/液態(tài)烴比值為8.72。以2 ℃/h速率升溫，當(dāng)熱解溫度達(dá)到528.6 ℃時(shí)，氣態(tài)烴產(chǎn)率達(dá)到90.23 mg/g，液態(tài)烴產(chǎn)率降為7.69 mg/g，氣態(tài)烴/液態(tài)烴比值為11.73(圖4)。

2.4.3 熱裂解生干氣階段

該階段Ro>3.50%，熱模擬溫度超過(guò)570 ℃。本階段有機(jī)質(zhì)熱演化產(chǎn)物主要為甲烷氣體，液態(tài)烴產(chǎn)物由于二次裂解消耗殆盡，氣態(tài)產(chǎn)物產(chǎn)率不斷增加，重?zé)N氣在高溫條件下裂解生成大量甲烷。

從山西組頁(yè)巖生烴模式圖(圖4)可以看出，有機(jī)質(zhì)熱演化過(guò)程中，生成的產(chǎn)物類(lèi)型及產(chǎn)物間的轉(zhuǎn)化過(guò)程與溫度密切相關(guān)。隨溫度增加，總烴、氣態(tài)烴和液態(tài)烴持續(xù)增加，液態(tài)烴的生成存在峰值；當(dāng)液態(tài)烴產(chǎn)率達(dá)到峰值后，隨溫度增加，液態(tài)烴二次裂解形成氣態(tài)產(chǎn)物，其自身不斷消耗，產(chǎn)率持續(xù)降低。

升溫速率對(duì)熱模擬產(chǎn)物產(chǎn)率存在影響。低升溫速率條件下，達(dá)到同樣的反應(yīng)溫度需要更長(zhǎng)的熱模擬反應(yīng)時(shí)間。與高升溫速率相比，同一溫度下，低升溫速率生成的氣態(tài)烴和總烴產(chǎn)率較高。對(duì)于液態(tài)烴，生油高峰前，同一溫度下，低升溫速率生成液態(tài)烴產(chǎn)率較高；生油高峰后，同一溫度下，由于低升溫速率增加了液態(tài)烴二次裂解時(shí)間，造成液態(tài)烴大量消耗，產(chǎn)率較低。

3 結(jié)論

(1)升溫速率是氣態(tài)烴產(chǎn)率和液態(tài)烴產(chǎn)率重要影響因素，然而對(duì)熱解氣碳同位素影響不明顯。相對(duì)低的升溫速率，需要持續(xù)較長(zhǎng)的升溫時(shí)間，因此烴類(lèi)產(chǎn)率較高，該實(shí)驗(yàn)結(jié)果充分證實(shí)了生烴過(guò)程中時(shí)間對(duì)溫度的補(bǔ)償作用。溫度對(duì)熱解氣碳同位素有明顯影響，δ13C1值隨溫度升高先變輕后變重，而δ13C2和δ13C3總體上隨溫度升高而變重。

(2)根據(jù)生烴熱模擬實(shí)驗(yàn)建立了頁(yè)巖生烴模式，將有機(jī)質(zhì)生烴過(guò)程分為3個(gè)演化階段。熱催化生油氣階段(Ro<1.20%)以形成液態(tài)烴為主，伴生少量氣，當(dāng)Ro=1.20%時(shí)達(dá)到生液態(tài)烴高峰。熱裂解生濕氣階段(Ro=1.20%～3.50%)以形成濕氣為主，伴生有少量輕質(zhì)油，氣態(tài)烴產(chǎn)率明顯增大，該階段形成的氣可能包括了干酪根熱裂解氣、液態(tài)烴裂解氣和重?zé)N氣二次裂解氣。熱裂解生干氣階段(Ro>3.50%)基本全為氣態(tài)產(chǎn)物，氣態(tài)烴產(chǎn)物成分主要為甲烷氣體，以重?zé)N氣的二次裂解為主，氣體為干氣。

(3)根據(jù)生烴熱模擬結(jié)果，在559.7 ℃時(shí)，甲烷產(chǎn)率可以達(dá)到119.79 mg/g，C1-5產(chǎn)率達(dá)到120.35 mg/g，表明山西組頁(yè)巖生烴潛力較好。

參考文獻(xiàn):

[1] 鄒艷榮,帥燕華,孔楓,等.油氣生成過(guò)程實(shí)驗(yàn)研究的思考與展望[J].石油實(shí)驗(yàn)地質(zhì),2004,26(4):375-382.

ZOU Yanrong,SHUAI Yanhua,KONG Feng,et al.Experiments on petroleum generation：considerations and outlook[J].Petro-leum Geology & Experiment,2004,26(4):375-382.

[2] 湯慶艷,張銘杰,張同偉,等.生烴熱模擬實(shí)驗(yàn)方法述評(píng)[J].西南石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,35(1):52-62.

TANG Qingyan,ZHANG Mingjie,ZHANG Tongwei,et al.A review on pyrolysis experimentation on hydrocarbon generation[J].Journal of Southwest Petroleum University (Science & Technology Edition),2013,35(1):52-62.

[3] 王云鵬,耿安松,劉德漢,等.頁(yè)巖、煤、瀝青和原油的生氣實(shí)驗(yàn)研究[J].沉積學(xué)報(bào),2004,22(S1):106-109.

WANG Yunpeng,GENG Ansong,LIU Dehan,et al.Comparison of experimental results between two simulation methods[J].Acta Sedimentologica Sinica,2004,22(S1):106-109.

[4] 劉德漢,張惠之,戴金星,等.煤巖顯微組分的成烴實(shí)驗(yàn)研究與評(píng)價(jià)[J].科學(xué)通報(bào),2000,45(4):346-352.

LIU Dehan,ZHANG Huizhi,DAI Jinxing,et al.Experimental study and evaluation on hydrocarbon generation of macerals[J].Chinese Science Bulletin,2000,45(14):1270-1276.

[5] 賈望魯,帥艷華,彭平安,等.塔北輪南地區(qū)原油瀝青質(zhì)生烴動(dòng)力學(xué)研究[J].科學(xué)通報(bào),2004,49(S1):76-80.

JIA Wanglu,SHUAI Yanhua,PENG Ping’an,et al.Kinetic study of hydrocarbon generation of oil asphaltene from Lunnan area,Tabei uplift[J].Chinese Science Bulletin,2004,49(S1):83-88.

[6] LEWAN M D,WINTERS J C,MCDONALD J H.Generation of oil-like pyrolyzates from organic-rich shales[J].Science,1979,203(4383):897-899.

[7] 侯讀杰.熱模擬實(shí)驗(yàn)技術(shù)簡(jiǎn)介[J].國(guó)外地質(zhì)勘探技術(shù),1990(10):39-40.

HOU Dujie.Brief introduction of thermal simulation experiment technology[J].Foreign Geoexploration Technology,1990(10):39-40.

[8] 姜峰,張友聯(lián),杜建國(guó).油氣生成熱模擬實(shí)驗(yàn)研究進(jìn)展[J].地球科學(xué)進(jìn)展,1996,11(5):453-459.

JIANG Feng,ZHANG Youlian,DU Jianguo.Advance of pyrolysis experimentation on hydrocarbon genesis[J].Advance in Earth Sciences,1996,11(5):453-459.

[9] BEHAR F,KRESSMANN S,RUDKIEWICZ J L,et al.Experimental simulation in a confined system and kinetic modeling of kerogen and oil cracking[J].Organic Geochemistry,1992,19(1/3):173-189.

[10] LEWAN M D.Laboratory simulation of petroleum formation[M]//ENGEL M H,MACKO S A.Organic geochemistry.Boston,MA:Springer,1993:419-442.

[11] 高喜龍,肖賢明,劉中云,等.用開(kāi)放體系的熱解方法對(duì)烴源巖生烴動(dòng)力學(xué)的研究:以東營(yíng)凹陷某生油巖為例[J].地球化學(xué),2003,32(5) :485-490.

GAO Xilong,XIAO Xianming,LIU Zhongyun,et al.Hydrocarbon generation kinetics of source rocks in open system based on pyrolysis simulation experiment:A case study of the Dongying Depression[J].Geochimica,2003,32(5):485-490.

[12] PAN C C,JIANG L L,LIU J Z,et al.The effects of calcite and montmorillonite on oil cracking in confined pyrolysis experiments[J].Organic Geochemistry,2010,41(7):611-626.

[13] MA X X,ZHENG J J,ZHENG G D,et al.Influence of pyrite on hydrocarbon generation during pyrolysis of type-Ⅲ kerogen[J].Fuel,2016,167:329-336.

[14] HAKIMI M H,ABDULLAH W H.Thermal maturity history and petroleum generation modeling for the Upper Jurassic Madbi source rocks in the Marib-Shabowah Basin,western Yemen[J].Marine and Petroleum Geology,2015,59:202-216.

[15] 余曉露,馬中良,鄭倫舉,等.不同熱模擬方式下烴源巖干酪根演化特征紅外光譜分析[J].石油實(shí)驗(yàn)地質(zhì),2017,39(1):134-140.

YU Xiaolu,MA Zhongliang,ZHENG Lunju,et al.FTIR analyses of source rock kerogen from different hydrous pyrolysis experiments[J].Petroleum Geology & Experiment,2017,39(1):134-140.

[16] 徐進(jìn)軍,金強(qiáng),程付啟,等.渤海灣盆地石炭系—二疊系煤系烴源巖二次生烴研究進(jìn)展與關(guān)鍵問(wèn)題[J].油氣地質(zhì)與采收率,2017,24(1):43-49.

XU Jinjun,JIN Qiang,CHENG Fuqi,et al.Advances and crucial issues on secondary hydrocarbon generation of the Carboniferous-Permian coal-measure source rocks in Bohai Bay basin[J].Petro-leum Geology and Recovery Efficiency,2017,24(1):43-49.

[17] 王寧,李榮西,王香增,等.海陸過(guò)渡相頁(yè)巖氣形成熱模擬實(shí)驗(yàn)研究[J].天然氣地球科學(xué),2016,27(1):189-197.

WANG Ning,LI Rongxi,WANG Xiangzeng,et al.Pyrolytic study on the gas-generating process of transitional shale[J].Natural Gas Geoscience,2016,27(1):189-197.

[18] 徐學(xué)敏,汪雙清,孫瑋琳,等.一種頁(yè)巖含氣性熱演化規(guī)律研究的模擬實(shí)驗(yàn)方法[J].巖礦測(cè)試,2016,35(2):186-192.

XU Xuemin,WANG Shuangqing,SUN Weilin,et al.Thermal simulation experiment for evaluating the influence of thermal evolution on gas-bearing properties of shale[J].Rock and Mineral Analysis,2016,35(2):186-192.

[19] 鄭倫舉,關(guān)德范,郭小文,等.影響海相烴源巖熱解生烴過(guò)程的地質(zhì)條件[J].地球科學(xué)(中國(guó)地質(zhì)大學(xué)學(xué)報(bào)),2015,40(5):909-917.

ZHENG Lunju,GUAN Defan,GUO Xiaowen,et al.Key geolo-gical conditions affecting pyrolysis experiments of marine source rocks for hydrocarbon generation[J].Earth Science (Journal of China University of Geosciences),2015,40(5):909-917.

[20] 薛蓮花,楊巍,仲佳愛(ài),等.富有機(jī)質(zhì)頁(yè)巖生烴階段孔隙演化:來(lái)自鄂爾多斯延長(zhǎng)組地質(zhì)條件約束下的熱模擬實(shí)驗(yàn)證據(jù)[J].地質(zhì)學(xué)報(bào),2015,89(5):970-978.

XUE Lianhua,YANG Wei,ZHONG Jia’ai,et al.Porous evolution of the organic-rich shale from simulated experiment with geological constrains,samples from Yanchang Formation in Ordos Basin[J].Acta Geologica Sinica,2015,89(5):970-978.

[21] 陶樹(shù),湯達(dá)禎,王東營(yíng),等.低成熟油頁(yè)巖的生排烴作用實(shí)驗(yàn)?zāi)M[J].地學(xué)前緣,2009,16(3):356-363.

TAO Shu,TANG Dazhen,WANG Dongying,et al.Simulation of hydrocarbon generation and expulsion of low maturity oil shale[J].Earth Science Frontiers,2009,16(3):356-363.

[22] 董春梅,馬存飛,欒國(guó)強(qiáng),等.泥頁(yè)巖熱模擬實(shí)驗(yàn)及成巖演化模式[J].沉積學(xué)報(bào),2015,33(5):1053-1061.

DONG Chunmei,MA Cunfei,LUAN Guoqiang,et al.Pyrolysis simulation experiment and diagenesis evolution pattern of shale[J].Acta Sedimentologica Sinica,2015,33(5):1053-1061.

[23] 劉大錳,金奎勵(lì),姚素平.沉積有機(jī)質(zhì)的成烴熱模擬實(shí)驗(yàn)研究[J].地質(zhì)論評(píng),1995,41(6):544-552.

LIU Dameng,JIN Kuili,YAO Suping.Hydrocarbon-generating thermal simulation experimental study of sedimentary organic matter[J].Geological Review,1995,41(6):544-552.

[24] 李斌,孟自芳,張百濤,等.含鈾物質(zhì)對(duì)有機(jī)質(zhì)生烴的實(shí)驗(yàn)研究[J].礦物學(xué)報(bào),2008,28(4):421-425.

LI Bin,MENG Zifang,ZHANG Baitao,et al.Experimental research on the influence of U-bearing minerals on hydrocarbon-generation of organic matter[J].Acta Mineralogica Sinica,2008,28(4):421-425.

[25] 姜蘭蘭,潘長(zhǎng)春,劉金鐘.礦物對(duì)原油裂解影響的實(shí)驗(yàn)研究[J].地球化學(xué),2009,38(2):165-173.

JIANG Lanlan,PAN Changchun,LIU Jinzhong.Experimental study on effects of minerals on oil cracking[J].Geochimica,2009,38(2):165-173.

[26] 董鵬,李賢慶,仰云峰,等.原油裂解成氣動(dòng)力學(xué)研究進(jìn)展[J].礦物巖石地球化學(xué)通報(bào),2009,28(2):201-208.

DONG Peng,LI Xianqing,YANG Yunfeng,et al.The advancement of kinetic study on crude oil-cracking gas[J].Bulletin of Mineralogy,Petrology and Geochemistry,2009,28(2):201-208.

[27] 肖芝華,胡國(guó)藝,李志生.從烴源巖熱模擬實(shí)驗(yàn)討論其生烴特征[J].天然氣地球科學(xué),2008,19(4):544-547.

XIAO Zhihua,HU Guoyi,LI Zhisheng.An analysis of characteri-stics of hydrocarbon generation from pyrolysis experiment of source rock[J].Natural Gas Geoscience,2008,19(4):544-547.

[28] 劉金鐘,唐永春.用干酪根生烴動(dòng)力學(xué)方法預(yù)測(cè)甲烷生成量之一例[J].科學(xué)通報(bào),1998,43(11):1187-1191.

LIU Jinzhong,TANG Yongchun.One example of predicting methane generation yield by hydrocarbon generating kinetics[J].Chinese Science Bulletin,1998,43(11):1187-1191.

[29] SHUAI Yanhua,PENG Ping’an,ZOU Yanrong,et al.Kinetic modeling of individual gaseous component formed from coal in a confined system[J].Organic Geochemistry,2006,37(8):932-943.

[30] SWEENEY J J,BURNHAM A K.Evaluation of a simple model of vitrinite reflectance based on chemical kinetics (1)[J].AAPG Bulletin,1990,74(10):1559-1570.

[31] SWEENEY J J,BURNHAM A K,BRAUN R L.A model of hydrocarbon generation from type Ⅰ kerogen:Application to Uinta Basin,Utah[J].AAPG Bulletin,1987,71(8):967-985.

[32] ZUMBERGE J,FERWORN K,BROWN S.Isotopic reversal (‘rollover’) in shale gases produced from the Mississippian Barnett and Fayetteville formations[J].Marine and Petroleum Geology,2012,31(1):43-52.

[33] TANG Y,PERRY J K,JENDEN P D,et al.Mathematical modeling of stable carbon isotope ratios in natural gases[J].Geochimica et Cosmochimica Acta,2000,64(15):2673-2687.

[34] TANG Y,XIA X.Kinetics and mechanism of shale gas formation:A quantitative interpretation of gas isotope “rollover” for shale gas formation[C]//AAPG Hedberg Conference.Austin,Texas:AAPG,2010:5-10.