崔永強 孟德福 李娜娜

1．大慶油田有限責任公司勘探開發(fā)研究院，黑龍江大慶，163712；2．大慶油田化工有限公司醋酸分公司，黑龍江大慶，163411

C1-C4正構(gòu)烷烴產(chǎn)物碳同位素的分布是研究中最重要的一個方面。一般認為：沉積有機質(zhì)經(jīng)降解作用生成簡單的烷烴氣體，其碳同位素組成受同位素動力學分餾效應(yīng)的制約。12C-12C鍵較12C-13C鍵弱將優(yōu)先斷裂，輕的碳同位素(12C)將優(yōu)先進入烷烴氣體中，從而構(gòu)成隨烷烴分子量(碳數(shù))增加、δ13C值變大的正序分布特征，即δ13C1＜δ13C2＜δ13C3＜δ13C4【1】。相反由簡單含碳分子(CO、CO2、CH4)聚合生成烷烴氣體的非生物過程，受同位素動力學分餾效應(yīng)的制約，12C-12C鍵優(yōu)先于12C-13C鍵斷裂，輕的碳同位素(12C)將優(yōu)先進入聚合形成的烷烴氣體中，而使初始形成的殘留烷烴氣體保留較重的δ13C值，后生成的烷烴氣體具有輕的 δ13C 值，構(gòu)成δ13C1＞δ13C2＞δ13C3＞δ13C4的反序分布特征【2】。

1 幔源天然氣碳同位素特征

關(guān)于天然氣甲烷同系物碳同位素值反序排列特征研究，國內(nèi)外學者（Des Marais， D.J.， and Moore，J.G.，1984；Chang，S.，Des Marais，D.J.， Mack，R.，1983；Szatmari P.，1989；Anderson R.B.，1984；戴金星等，2009；2008；2006；2003；2001；1997；王先彬，2009；2008；2006；2000；1997；郭占謙，1994；1997；2000）開展了大量的工作，認為碳同位素反序是天然氣地幔來源的直接證據(jù)【3】。

由 C1單元聚合形成高分子烴類的過程中將會出現(xiàn) δ13C1＞δ13C2＞δ13C3＞δ13C4的碳同位素反序排列分布模式，這已被火花放電促使甲烷生成乙烷、丙烷的實驗所證實【4，5】。南格林蘭堿性流體包裹體中甲烷、乙烷、丙烷的碳同位素組成具有這種反序排列特征，暗示無機成因甲烷及其同系物的碳同位素是反序排列的。在地球排氣過程中，甲烷流自然可通過費-托反應(yīng)（Fisher-Tropsch）過程形成各種各樣的烴類分子【6，7】。天然油氣中見到的分子的特殊分布與這種過程產(chǎn)生的分布顯著地相似【8】。戴金星等【9】明確將松遼盆地徐家圍子斷陷產(chǎn)自火山巖的深層烴類氣確定為幔來源。

2 費-托合成天然氣碳同位素特征

胡桂興等【10】指出，在以Co為催化劑的費-托合成過程中，C2-C4的同位素組成呈現(xiàn)出隨碳數(shù)增加而逐漸降低的趨勢。即 δ13C2＞δ13C3＞δ13C4，并且隨著碳數(shù)增大，相鄰碳數(shù)的烷烴之間的同位素分餾值變小，但反序排列特征明顯。因此，胡桂興等認為費-托合成過程形成的乙烷、丙烷和丁烷的碳同位素組成與生物成因的烷烴同系物的碳同位素組成相反。即具有反序排列的碳同位素分配模式。不過.最近也有研究發(fā)現(xiàn).費-托合成形成的甲烷及其同系物之間的碳同位素分配模式更加趨向于正序排列，即 δ13C1＜δ13C2＜δ13C3＜δ13C4【11】。

Taran等【12】發(fā)現(xiàn)同位素效應(yīng)依賴于 CO的轉(zhuǎn)化率。在穩(wěn)態(tài)條件(CO轉(zhuǎn)化率大于 90%)下，CO與CO2、CO與烴類之間發(fā)生了強烈的動力學分餾效應(yīng).分餾值分別大約為-10×10-3和 38×10-3，而且甲烷、乙烷和高分子烴類之間碳同位素組成沒有顯著的差異；在CO低轉(zhuǎn)化率的條件下C1-C4烴類之間存在明顯的負碳同位素序列排列特征?？傮w上，甲烷與乙烷之間的同位素分餾與CO的轉(zhuǎn)化率有良好的相關(guān)性，即在低轉(zhuǎn)化率下乙烷的碳同位素組成比甲烷輕大約6×10-3，而在穩(wěn)態(tài)高轉(zhuǎn)化率下乙烷的碳同位素組成比甲烷重大約2×10-3。

楊春等【13】認為，40多年來的費-托合成模擬實驗表明，只有部分結(jié)果出現(xiàn)了烷烴氣碳同位素反序現(xiàn)象。通過較為完整的分析發(fā)現(xiàn)，其原因可能是：隨著模擬實驗時間的增加與溫度的升高并隨著產(chǎn)物轉(zhuǎn)化率的增高，控制烷烴氣碳同位素分餾的因素逐漸由動力學機制轉(zhuǎn)變?yōu)闊崃W平衡，烷烴氣碳同位素序列將由反序轉(zhuǎn)變?yōu)椴糠值罐D(zhuǎn)再到類似常規(guī)大然氣的正序分布。只有在較短時間(轉(zhuǎn)化率較低)或者開放體系(隨生隨排)條件下才遵從動力學分餾。電火花放電合成實驗僅代表了理想狀態(tài)下動力學分餾過程。

隨著研究的深入和資料的積累，結(jié)合最近研究者【11，12，14】得出的結(jié)果(McCollom 和 Qi Fu 等【11】，nTaran 等【12】， Seewald【14】)，隨著模擬溫度的升高，楊春等【13】發(fā)現(xiàn)δ13C co2-c1漸減小，并逐漸趨向熱力學平衡的分餾數(shù)值范圍(圖1)。

圖1 實驗?zāi)M溫度與δ13Cco2-c1關(guān)系Fig.1 The temperatures and δ13Cco2-c1 in the simulation experiment

Galimov【15】利用 Des Marais等【16】的火花放電實驗數(shù)據(jù)對原始圖進行了修改，電火花實驗中合成的C1-C3烴類與火山巖中C1-C3烴類的同位素分布具有相同的模式(圖2)。圖2中，巖漿來源的 C1-C4烴類的碳同位素值與在 400℃時根據(jù)Galimov【15】提出的同位素動力學方法計算的理論曲線"A”非常相似，因為該計算就是在復雜的瞬間理論框架下完成的。

圖2 電火花放電C1-C3與有機質(zhì)熱分解正構(gòu)烷烴產(chǎn)物碳同位素的分布特征( Galimov【15】）Fig.2 The carbon isotopic distribution patterns of C1-C3 hydrocarbons synthesized in a spark-discharge experiment and C1-C4 hydrocarbon gases produced by cracking (in sedimentary rocks) (Galimov【15】)

將Hu等【10】(圖3a)，Taran等【12】(圖 3b) ，Qi Fu等【11】(圖3c)的實驗結(jié)果按照反應(yīng)時間長短順序放在同一張圖上(圖3)，可以看出如下趨勢：隨著時間的增加，反應(yīng)物轉(zhuǎn)化率的提高，烷烴氣碳同位素序列逐漸由反序逐漸轉(zhuǎn)為倒轉(zhuǎn)然后再轉(zhuǎn)變?yōu)檎颉u等【10】反應(yīng)時間最短的第一個實驗結(jié)果是一個典型的碳同位素反序分布【10】。Taran等【12】的實驗是在開放體系下進行，隨著反應(yīng)時間的增加和轉(zhuǎn)化率的提高， 碳同位素序列從反序逐漸過渡為局部倒轉(zhuǎn)，較好的反映了這一轉(zhuǎn)變的過程。 Qi Fu等【11】實驗反映時間最長并且是在封閉體系下進行，實驗的結(jié)果類似正常天然氣的正序分布。Des Marais等【16】電火花實驗的時間最短，而烷烴氣碳同位素序列也更符合理論上的反序關(guān)系(圖2)，也就是說，隨著反應(yīng)時間的增加，費-托合成無機烷烴氣的合成將逐漸由動力學過程轉(zhuǎn)變?yōu)闊崃W平衡過程。這樣才能更加合理解釋碳同位素分餾理論與實驗現(xiàn)象之間的矛盾。結(jié)合圖1中溫度與δ13Cco2-c1關(guān)系，隨著溫度的提高和反應(yīng)時間的增長，δ13Cco2與δ13C1之間的分餾值逐漸減小，更趨近于熱力學平衡分餾的范圍。因此，在實驗中想預(yù)期得到同位素反序的結(jié)果需要在實驗設(shè)計時盡量提供滿足動力學反應(yīng)的實驗條件，例如在開放體系下隨生隨排的實驗條下。

圖3 費-托合成實驗烷烴氣碳同位素序列Fig.3 δ13c values of n-alkanes distribution patterns in the diferrent Fischer-Tropsch Synthetic experiments

3 芳深9井CO2加H2的費-托合成實驗

大慶油田研究院2009年，利用芳深9井CO2制備實驗材料，以金管體系開展費-托合成實驗和產(chǎn)物的組成、同位素分析研究?大慶油田勘探開發(fā)研究院2009完成內(nèi)部報告《徐家圍子斷陷無機成因烴類氣判別》合同編號：2008-JS-934。在不同溫壓條件下進行了費-托合成實驗。30MPa下分別在350℃、380℃和 390℃的合成實驗；390℃恒溫下，開展了不同壓力（30MPa、100 MPa和200MPa）的實驗；研究了H2對費-托合成的影響。30MPa不同溫度下費-托合成結(jié)果，包括產(chǎn)率和碳同位素數(shù)據(jù)，分別列于表 1至表 3中。390℃不同壓力（30MPa、100MPa和 200MPa）的費-托合成實驗結(jié)果列于表4。

表1 30MPa、350℃下費-托合成實驗結(jié)果Table 1 The test results of F-T synthesis in 30MPa and 350℃

表2 30MPa、380℃下費-托合成實驗結(jié)果Table 2 The test results of F-T synthesis in 30MPa and 380℃

表3 30MPa、390℃下費-托合成實驗結(jié)果Table 3 The test results of F-T synthesis in 30MPa and 390℃

表4 390℃、100MPa和200MPa下費-托合成實驗結(jié)果Table 4 The test results of F-T synthesis in 100MPa, 200MPa and 390℃

3.1 費-托合成氣組分變化

30MPa、350至390℃下，合成天然氣產(chǎn)率隨時間的變化示于圖4。很明顯，隨時間的增加，合成天然氣的產(chǎn)率增加，特別是甲烷產(chǎn)率變化非常明顯。

圖4 30MPa、380℃下合成氣產(chǎn)率隨時間的變化Fig. 4 Change of the Synthesis gas yield with the reaction time in 30MPa and 380℃

在烴氣產(chǎn)物中，甲烷占絕對優(yōu)勢，據(jù)30MPa、350～390℃下費-托合成天然氣數(shù)據(jù)的統(tǒng)計，甲烷在烴氣中占85%～95%之間，天然氣顯示為濕氣特征；高溫下，合成氣產(chǎn)率明顯高于低溫下的產(chǎn)率。

3.2 天然氣碳同位素變化

在相同壓力（30MPa）、350～390℃下合成的天然氣，天然氣碳同位素主要呈現(xiàn)正序或者局部反序，沒有發(fā)現(xiàn)碳同位素完全反序現(xiàn)象。在相對較高溫度下，碳同位素局部反序現(xiàn)象有增多的趨勢，且主要為甲烷、乙烷碳同位素反序（圖5）。350℃下合成的天然氣，天然氣碳同位素均為正序，沒有發(fā)生碳同位素局部反序現(xiàn)象。

高壓（100～200MPa）下合成的天然氣，其碳同位素主要表現(xiàn)為正序，特別是甲烷、乙烷間碳同位素為正序；部分C2、C3同位素輕微的局部反序、部分C2、C3同位素為正序（圖5）；個別樣品碳同位素顯示接近負碳系列（反序）。這些實驗結(jié)果表明，壓力不但影響合成天然氣的產(chǎn)率，而且對天然氣碳同位素樣式具有很大影響。低壓下C2、C3同位素分布范圍相對較寬，個別的同位素很重，而高壓下，C2、C3同位素相對集中，分布在-30×10-2～-25×10-3之間（圖 5）。

圖5 30MPa、390℃下合成氣碳同位素的變化Fig.5 Change of the Synthesis gas carbon isotope with the reaction time in 30MPa and 390℃

幔源氣是在地幔高溫高壓條件下形成的，在非常高的壓力下，烴氣的碳同位素是否會完全反序？這取決于高溫高壓下C-H-O平衡體系中，非生物甲烷向 C2+烴轉(zhuǎn)化導致烴鏈增長程度有關(guān)。不能排除幔源氣具有碳同位素完全反序的可能性。這并非是問題的關(guān)鍵，問題的關(guān)鍵在于幔源氣運移倒盆地儲層時，烴氣碳同位素是否為完全反序。根據(jù)最新模擬計算成果（Zhang & Duan，2009），幔源氣從地幔向地殼運移過程中，乙烷很可能被氧化，幔源流體上升到巖石圈-軟流圈邊界時，流體成分由主要為H2O-CH4-H2-C2H6混合物變成主要H2O-CO2-CO混合物組成。

3.3 溫度、壓力對合成氣的影響

前已述及，隨時間的增加，合成烴氣產(chǎn)率增加。圖6直觀展示了30MPa下，不同溫度、時間對合成烴氣的影響。很明顯，合成經(jīng)歷的時間越長、溫度越高，烴氣（C1-C5）產(chǎn)率越高，而且產(chǎn)率隨溫度的增加呈非線性增長。以恒溫時間36h為例，從350℃增加到 380℃時，溫度增加了 30℃，烴氣產(chǎn)率增加不足0.8mL/g；而溫度從380℃增至390℃時，溫度僅增加10℃，烴氣產(chǎn)率增加了1.5mL/g。可見，溫度對合成氣產(chǎn)率具有很重要的影響。溫度越高，合成烴氣越多。換句話說，在相同壓力、時間條件下，溫度越高CO2越容易合成烴氣；反之，溫度越低越不利于合成烴氣。

圖6 30MPa不同溫度/時間下費-托合成氣產(chǎn)率的變化Fig.6 Change of the Synthesis gas yield with the reaction time and temperature in 30MPa

對于大多數(shù)含油氣沉積盆地而言，天然氣儲層溫度不會超過 300℃。這意味著，盡管地質(zhì)時間十分漫長，地幔中合成的烴氣遠高于在地殼或者儲層中合成的非生物氣量。由于在幔源氣在地幔中已達到平衡狀態(tài)，合成反應(yīng)基本完成，在向地殼運移過程中，由于溫度逐漸降低很可能不會或很少有新的合成反應(yīng)發(fā)生。

壓力對合成氣也有重要影響。圖 7顯示了390℃不同壓力下合成甲烷和烴氣（C1-C5）產(chǎn)率的影響。在18～54h的時間范圍內(nèi)，高壓下合成的甲烷和烴氣是低壓下的 3～7倍。而 100MPa和200MPa下的烴氣產(chǎn)率差別不大，暗示200MPa以上的壓力對非生物氣形成影響較??；壓力的影響主要體現(xiàn)在200MPa以下的壓力范圍內(nèi)。

圖7 390℃不同壓力下合成甲烷（上）和烴氣（下）產(chǎn)率的變化Fig.7 Yield change of synthesis methane(upper) and hydrocarbon gas(under) with the reaction time and pressure in 390℃

不同壓力下的恒溫合成實驗結(jié)果，揭示了壓力對費-托合成天然氣的影響，意味著在高壓條件下，費-托合成反應(yīng)很可能達到平衡狀態(tài)。暗示在地幔高壓條件下，C-O-H體系中化學熱力學已經(jīng)平衡，在降溫降壓過程中，很難或很少發(fā)生新的費-托合成反應(yīng)。這與不同溫度實驗結(jié)果相呼應(yīng)，是以前研究沒有注意到的。

3.4 H2/CO2摩爾比對碳同位素的控制

雖然不同溫壓條件下，烴氣的碳同位素變化的總體趨勢是清楚的。然而，其同位素變化樣式卻顯得有些復雜。前人的實驗雖然也觀察到相似的結(jié)果，但未進行深入研究，現(xiàn)有文獻亦很少涉及。本項研究中，在文獻調(diào)研的基礎(chǔ)上，采用H2含量控制實驗，有可能進一步揭示控制烴氣碳同位素樣式變化的本質(zhì)因素。

時間、溫度、壓力無疑都是影響費-托合成天然氣同位素的重要因素，但并非本質(zhì)因素。大多數(shù)研究常常把重點放在時、溫、壓的影響上，沒有關(guān)注到隱含在這些影響背后的本質(zhì)因素，這也是合成天然氣碳同位素及樣式變化沒有得到很好解釋的主要原因。

圖8展示了不同壓力390℃下合成甲烷碳同位素隨H2/CO2摩爾比的變化?？梢?，在相同溫度、不同壓力條件下，合成甲烷的碳同位素明顯受反應(yīng)物H2/CO2摩爾比的控制。H2/CO2摩爾比越大，甲烷碳同位素越重。

圖8 390℃不同壓力下合成甲烷碳同位素隨H2/CO2比的變化Fig.8 Change of the synthesis methane carbon isotope with the H2/CO2 ratio and pressure in 390℃

圖 9為不同溫度相同壓力下合成甲烷碳同位素隨H2/CO2摩爾比的變化。同樣，合成甲烷碳同位素也受到 H2/CO2摩爾比的控制：在高H2/CO2摩爾比下，甲烷碳同位素重。H2/CO2摩爾比越大，反應(yīng)達到平衡是消耗的 CO2越多。根據(jù)化學反應(yīng)動力學原理，12C-CO2較13C-CO2優(yōu)先與H2發(fā)生費-托反應(yīng)，生成碳同位素相對輕的甲烷，隨著反應(yīng)的進行，更多的13C-CO2加入到合成反應(yīng)，合成的甲烷碳同位素越來越重。當H2/CO2摩爾比接近4∶1、反應(yīng)達到平衡時，CO2將耗盡，這時近 4∶1、反應(yīng)達到平衡時，CO2將耗盡，這時CO2-CH4間分餾就會很小。

圖9 30MPa不同溫度下合成甲烷碳同位素隨H2/CO2比的變化Fig.9 Change of the synthesis methane carbon isotope with the H2/CO2 ratio and temperature in 30MPa

H2/CO2摩爾比不但控制了合成甲烷的碳同位素，而且控制烴氣的同位素樣式。圖10展示了30MPa不同溫度下合成天然氣碳同位素隨H2/CO2摩爾比的變化；圖11為390℃不同壓力下合成天然氣碳同位素隨 H2/CO2摩爾比的變化規(guī)律。為了直觀反映反應(yīng)物H2/CO2摩爾比對合成烴氣碳同位素洋式的控制，圖中左上角標明了H2/CO2摩爾比，以便比較溫度、壓力對碳同位素樣式的影響。

結(jié)合圖10和圖11可見，H2/CO2摩爾比是控制合成烴氣的根本因素。目前的實驗結(jié)果，在低壓下：

當H2/CO2摩爾比小于1.2時，合成烴氣的碳同位素隨碳數(shù)增加同位素線性變重。

當H2/CO2摩爾比大于1.2時，合成烴氣甲烷碳同位素與乙烷碳同位素大體相當，甚至有甲烷碳同位素＞乙烷碳同位素的現(xiàn)象；C2以后，隨碳數(shù)增加同位素值線性增加。加拿大前寒武地盾的非生物氣碳同位素表現(xiàn)為V字型分布，很可能與淺層低壓下H2/CO2比對碳同位素樣式的控制有關(guān)。

當H2/CO2摩爾比接近或大于4時，出現(xiàn)碳同位素完全反序。這很可能與H2/CO2完全反應(yīng)合成烴氣以后烴鏈增長有關(guān)。意味著鏈增長反應(yīng)速率小于合成甲烷的反應(yīng)速率。當甲烷合成反應(yīng)完成以后，烴鏈增長反應(yīng)成為優(yōu)勢反應(yīng)，由于12C甲烷優(yōu)先參與鏈增長過程，使得剩余甲烷碳同位素變重而C2+烴碳同位素相對較輕。

圖10 30MPa不同溫度下H2/CO2比對碳同位素樣式的影響Fig.10 Effect of carbon isotopic distribution patterns with the H2/CO2 ratios and temperatures in 30 MPa

圖11 390℃不同壓力下H2/CO2比對碳同位素樣式的影響Fig.11 Effect of carbon isotopic distribution patterns with the H2/CO2 ratios and pressures in 390℃

4 結(jié)論

徐家圍子斷陷芳深9井二氧化碳與氫合成的費-托合成實驗結(jié)果，與 Taran 等【12】、楊春等【13】、Galimov【15】、倪云燕等【18】的研究結(jié)果相印證。費-托合成氣的碳同位素特征控制因素為：

（1）碳同位素效應(yīng)依賴于CO和CO2的轉(zhuǎn)化率，轉(zhuǎn)化率低傾向于形成甲烷及其同系物碳同位素反序排列，轉(zhuǎn)化率高傾向形成正序排列；

（2）只有在較短的反應(yīng)時間內(nèi)形成的產(chǎn)物才能構(gòu)成反序排列。

芳深9井二氧化碳與氫合成的費-托合成實驗進一步得到：

（3）H2/CO2摩爾比是控制合成烴氣的根本因素。在低壓下，當H2/CO2摩爾比小于1.2時，合成烴氣的碳同位素隨碳數(shù)增加同位素線性變重；當H2/CO2摩爾比大于1.2時，合成烴氣甲烷碳同位素與乙烷碳同位素大體相當，甚至有甲烷碳同位素大于乙烷碳同位素的現(xiàn)象；C2以后，隨碳數(shù)增加同位素值線性增加。加拿大前寒武地盾的非生物氣碳同位素表現(xiàn)為V字型分布，很可能與淺層低壓下H2/CO2比對碳同位素樣式的控制有關(guān)。

當H2/CO2摩爾比接近或大于4時，出現(xiàn)碳同位素完全反序。這很可能與H2/CO2完全反應(yīng)合成烴氣以后烴鏈增長有關(guān)。意味著鏈增長反應(yīng)速率小于合成甲烷的反應(yīng)速率。當甲烷合成反應(yīng)完成以后，烴鏈增長反應(yīng)成為優(yōu)勢反應(yīng)，由于12C甲烷優(yōu)先參與鏈增長過程，使得剩余甲烷碳同位素變重而C2+烴碳同位素相對較輕。

因為不同作者、不同條件下費-托合成反應(yīng)的天然氣碳同位素值主要是正序排列，有學者質(zhì)疑碳同位素反序不作為無機合成烴的特征。上述研究已經(jīng)表明，費-托合成氣的正序和反序排列即受控于參加反應(yīng)氣體組分，又受控于熱動力學條件，正序和反序排列天然氣均是合成的產(chǎn)物。這個認識對石油天然氣成因的認識將產(chǎn)生巨大沖擊。因為費-托合成氣主要是正序排列的實驗結(jié)果與大多數(shù)天然氣藏表現(xiàn)出正序排列的地質(zhì)事實相符合。

1 Sherwood Lollar B， Frape S K，Weise S M， et al. Abiogenic methanogenesis in crystalline rocks[J]. Geochimica et Cosmochimica Acta， 1993， 57(23/24)：5087～5097

2 Lancet MS， Anders E. Carbon isotope fractionation in the Fischer-Tropsch synthesis and in meteorites[J].Science，1970，170(3961)：980～982

3 戴金星. 各類烷烴氣的鑒別［J］.中國科學（B輯），1992，（2）：185～193

4 Des Marais， D.J.， and Moore，J，G. Carbon and its istopes in mid-ocean basaltic glass［J］. Earth planet.Sci.Lett.，1984，69：43～57

5 Chang， S.， Des Marais， D.J.， Mack， R.， Miller， S.L. and Strathearn， G.E. Prebio-tic organic synthesis and the origin of life， In： Earth’s Earliest Biosphere： Its Origin & Evolution(ed.J.W.Schopf) ［M］. Princeton University Press， 1983.54～92

6 Szatmari P. Petroleum formation by Fischer-Tropsch synthesis in plate tectonics[J]. AAPG，1989，73(8)：989～998

7 Anderson R.B. The Fisher-Tropsch Systhesis[M]. New York：Academic Press，1984：176

8 Friedel， R.A.， and Sharkey， A.G. Alkanes in natural and synthetic petroleum：comparison of calculated and actural compositions [J］. Science，1963，124： 1203～1205

9 戴金星，胡國藝，倪云燕，等. 中國東部天然氣分布特征天然氣[M]. 地球科學，2009，(20 )4： 471～487

10 Hu G X，Ouyang Z Y. Carbon isotopic fractionation in the process of Fischer～Tropsch reaction in primitive solar nebula[J]. Science in China：Series D，1998，41(2)：202～207

11 Fu Q，Lollar B S， Horita J， et a1. Abiotic formation of hydrocarbon under hydrothermal conditions; Constraints from chemical and isotope data[J]. Geochimica et Cosmochimica Acta，2007，71(8)：1982～1998

12 T aran Yu ri A， George A Kl iger， Vyacheslav S Sevast ianov. Carbon i sot op e eff ect s in th e open- syst em Fischer～T ropsch s ynthes is[ J ] . Geochimica et C osmochim ica Act a，2007，71：4474～4487

13 楊春，王京紅，米敬奎，等. 費-=托合成實驗中烷烴氣碳同位素分餾機理[J]. 天然氣地球科學，2011，22（2）：206～210，223

14 McCollom T M ，Seewald J S. Carbon isotope composition of organic compounds produced by abiotic synthesis under hydrothermal conditions [ J ] . Earth and Planetary Science Letters，2006，243( 1/ 2) ：74～ 84

15 Galimov E M. Isotope organic geochemistry[ J] . Organic Geochemistry，2006，37：1200～1262

16 Des Marais D J， Donchim J H，Nehring N L， et al. Molecular carbon isotopic evidence for the origin of geothermal hydrocarbons[ J]. Nature，1981，292：826～828

17 Davis， B. H.， "The two～alpha value for iron Fischer-Tropsch catalysts： fact or fiction?，" ACS[J].， Division of Fuel Chemistry， Preprints，1992，37(1)：172～83.

18 倪云燕，靳永斌. 費-托合成反應(yīng)中的碳同位素分餾[J]. 石油勘探與開發(fā)，2011，38 (2)：229～256