• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    四川盆地海、陸相烴源巖有機(jī)質(zhì)穩(wěn)定碳同位素組成變化及其地球化學(xué)意義

    2018-01-04 09:00:28朱揚(yáng)明李穎郝芳鄒華耀郭旭升
    沉積學(xué)報(bào) 2017年6期

    朱揚(yáng)明,李穎,郝芳,鄒華耀,郭旭升

    1.浙江大學(xué)地球科學(xué)學(xué)院,杭州 310027 2.中國地質(zhì)大學(xué)(武漢),武漢 430074 3.中國石油大學(xué)(北京)地球科學(xué)學(xué)院,北京 102249 4.中國石化勘探分公司,成都 610000

    四川盆地海、陸相烴源巖有機(jī)質(zhì)穩(wěn)定碳同位素組成變化及其地球化學(xué)意義

    朱揚(yáng)明1,李穎1,郝芳2,鄒華耀3,郭旭升4

    1.浙江大學(xué)地球科學(xué)學(xué)院,杭州 310027 2.中國地質(zhì)大學(xué)(武漢),武漢 430074 3.中國石油大學(xué)(北京)地球科學(xué)學(xué)院,北京 102249 4.中國石化勘探分公司,成都 610000

    采用質(zhì)譜和色譜—同位素質(zhì)譜技術(shù),測定了500余個(gè)干酪根及飽和烴、芳烴組份、正構(gòu)烷烴的13C值,以此揭示了四川盆地海、陸相烴源巖有機(jī)碳同位素組成隨地質(zhì)時(shí)代的變化特點(diǎn)及其在高成熟階段的演化特征,并結(jié)合其他相關(guān)分析資料,應(yīng)用碳同位素剖析了不同類型海、陸相烴源巖的有機(jī)質(zhì)生源及沉積環(huán)境。研究結(jié)果表明,該盆地從震旦系燈影組到中三疊統(tǒng)雷口坡組的海相地層中,干酪根的碳同位素組成隨層位變新呈逐漸變重的趨勢,可能的原因之一是浮游植物的進(jìn)化作用;而上三疊統(tǒng)須家河組至中侏羅統(tǒng)千佛崖組的陸相有機(jī)質(zhì)碳同位素組成則有反向的年代變化,主要與生源構(gòu)成和沉積環(huán)境性質(zhì)有關(guān)。這些海、陸相烴源巖的有機(jī)源難以用干酪根13C值進(jìn)行區(qū)分,而可用飽和烴與芳烴組份的13C及其CV值來區(qū)別。海洋和湖泊不同沉積相帶烴源巖的干酪根具有明顯不同的13C值,結(jié)合其他相關(guān)資料可識別其有機(jī)質(zhì)來源。煤系地層中煤與泥巖在干酪根碳同位素組成上沒有可區(qū)分性,而兩者有不同的正構(gòu)烷烴碳同位素分布曲線。從成熟晚期到高—過成熟階段,海、陸相干酪根的碳同位素均變重1~2‰,煤系泥巖的正構(gòu)烷烴碳同位素分布曲線由負(fù)向線型分布向平直型演變。這些碳同位素的變化特征為油氣源對比和烴源識別提供了依據(jù)。

    烴源巖;碳同位素;干酪根;正構(gòu)烷烴;四川盆地

    0 引言

    沉積有機(jī)質(zhì)的穩(wěn)定碳同位素組成主要繼承于生物體,在成巖作用過程中變化不大,能用于烴源巖有機(jī)質(zhì)類型劃分[1-2]和油氣源對比[3],也可用來重建古環(huán)境[4-5]和區(qū)分海、陸相有機(jī)質(zhì)[6]。四川盆地海、陸相烴源巖有機(jī)質(zhì)的熱演化大都在高—過成熟階段,它們的常規(guī)生源和指相參數(shù)大多失效。其有機(jī)顯微組分的光學(xué)性質(zhì)趨同,難以判別其原始性質(zhì)[2];生物標(biāo)志物的組成和分布出現(xiàn)了異常變化,基本失去了其常規(guī)的地球化學(xué)意義[7]。因而在這樣的高熱演化地區(qū),烴源巖碳同位素組成的地球化學(xué)應(yīng)用顯得尤為重要。前人已在這方面做過不少研究工作。黃藉中[8]依據(jù)干酪根碳同位素值,結(jié)合其他參數(shù)劃分了海相烴源巖的有機(jī)相。黃藉中等[9]和梁狄剛等[2]借鑒陸相烴源巖有機(jī)質(zhì)類型的碳同位素劃分標(biāo)準(zhǔn),利用干酪根13C值確定海相烴源巖的有機(jī)質(zhì)類型。王順玉等[10]闡述了大巴山、米倉山南緣地區(qū)震旦系至上三疊統(tǒng)烴源巖干酪根的碳同位素組成及變化特征。本研究在四川盆地東北部、東南部及東部地區(qū)的主要探井較為系統(tǒng)地采集了各套海、陸相烴源巖樣品,分離出干酪根及飽和烴、芳烴組份,進(jìn)行穩(wěn)定碳同位素分析,擬揭示它們的分布特征和變化規(guī)律,并結(jié)合其他資料探討其生源及沉積環(huán)境意義,為油氣源對比和烴源識別提供依據(jù)。

    1 地質(zhì)背景與樣品

    四川盆地屬揚(yáng)子準(zhǔn)地臺西北隅的一個(gè)呈北東向延展的菱形狀次級構(gòu)造單元,是古生代克拉通盆地與中新生代前陸盆地的復(fù)合型盆地,發(fā)育有從震旦系至第四系厚逾13 000 m的沉積地層。其中,震旦系—中三疊統(tǒng)主要為海相沉積地層,上三疊統(tǒng)—第四系主要為陸相沉積地層[11]。根據(jù)基底性質(zhì)、沉積蓋層、氣藏特征及天然氣類型等,把四川盆地劃分為川東、川南、川西和川中4個(gè)油氣聚集區(qū)(構(gòu)造區(qū)塊,圖1)[12]。該盆地發(fā)育多套海、陸相烴源巖,在縱向上可劃分出下寒武統(tǒng)(∈1)、下志留統(tǒng)(S1)、下二疊統(tǒng)(P1)、上二疊統(tǒng)(P2)、上三疊統(tǒng)(T3)、下侏羅統(tǒng)(J1)等6套主要有效烴源巖。其中,∈1、S1、J1主要為泥質(zhì)烴源巖,P1主要為碳酸鹽巖烴源巖,P2、T3主要為含煤泥質(zhì)烴源巖[13]。此外,局部地區(qū)的中侏羅統(tǒng)千佛崖組、中三疊統(tǒng)雷口坡組和震旦系部分層段也具備烴源條件[8,14-15]。

    本研究在川東、川中、川西及川南地區(qū)的探井及露頭剖面采集了不同層系、不同巖性類型的烴源巖樣品500余個(gè),主要的取樣探井、剖面和地層及樣品數(shù)等見圖1和表1。相關(guān)地層烴源巖的熱演化程度都較高,除川中、川西部分地區(qū)的中、下侏羅統(tǒng)及上三疊統(tǒng)外,基本上均處于高—過成熟階段。這些烴源巖樣品經(jīng)TOC、瀝青A抽提等常規(guī)分析后,用酸解方式分離和提純出干酪根。部分樣品的瀝青A用石油醚沉淀去除瀝青質(zhì)后,在硅膠/氧化鋁層析柱上分別用石油醚和二氯甲烷分離出飽和烴和芳烴組份。一些飽和烴組份用尿素的甲醇飽和溶液進(jìn)行絡(luò)合,分離出正構(gòu)烷烴組份。干酪根的碳同位素分析用DELTA PLUS XL 穩(wěn)定同位素質(zhì)譜儀完成;正構(gòu)烷烴單體烴碳同位素組成的測定儀器為GV Isoprime 色譜—同位素質(zhì)譜儀,所測碳同位素比值的誤差為0.5‰(PDB)。為了作對比研究,還有6個(gè)千佛崖組泥巖樣品進(jìn)行了熱脫附輕烴氣相色譜分析,熱脫附溫度為250℃。此外,部分烴源巖干酪根碳同位素?cái)?shù)據(jù)收集于中石化勘探分公司及相關(guān)文獻(xiàn)。

    圖1 四川盆地構(gòu)造單元與主要取樣探井和剖面位置圖Fig.1 Distribution of tectonic units and some wells and profiles sampled in Sichuan Basin

    地層沉積相d13C范圍/‰平均值/‰樣品個(gè)數(shù)主要取樣探井或剖面侏羅系J2q湖相-23.2~-27.7-25.613普光5、元陸6、川石44等J1z湖相-22.5~-28.7-24.965普光3、元陸4、興隆1等三疊系T3x湖沼相-23.1~-27.0-24.570川科1、雷北1、元壩6、馬201等T2l淺海臺地相-23.6~-26.6-25.359川科1、毛壩3、元陸6、馬201等T1j-T1f淺海臺地-24.5~-31.0-26.121川岳83、元壩2、河壩1、建深1等二疊系P2c(d)臺地—陸棚相-24.9~-29.2-27.031川岳84、元壩2、馬2、龍8等P2l(w)海陸交互相-22.3~-28.7-25.8113見下川東北海灣潟湖相-24.0~-28.7-26.754普光5、元壩6、河壩1、金溪1等川東南近海湖沼相-22.3~-26.9-23.834隆盛1、丁山1、南川、習(xí)水等過渡區(qū)灰泥臺坪相-23.1~-28.2-26.625三星1、福石1、建深1、冷水溪P1m臺地—緩坡相-25.1~-29.6-27.651焦頁66、毛壩3、元壩3、冷水溪P1q臺地相-26.5~-29.5-27.924元壩3、龍8、冷水溪等志留系S1l陸棚相-25.2~-30.8-29.139建深1、焦頁1、丁山1、南江寒武系1q淺?!獫暫?26.6~-32.9-30.020丁山1等震旦系Z1d淺海相-31.6~-36.7-34.98大巴山-米倉山南緣?

    注:數(shù)據(jù)引自文獻(xiàn)[10]

    2 不同地質(zhì)時(shí)代烴源巖干酪根碳同位素組成變化

    2.1 海相地層干酪根

    分析數(shù)據(jù)表明,四川盆地在震旦系燈影組到中三疊統(tǒng)雷口坡組的海相地層中,除海陸交互相的上二疊統(tǒng)龍?zhí)?P2l)/吳家坪(P2w)組外,其干酪根碳同位素組成隨層位變新也呈逐層變重的趨勢(圖2),具有海相沉積有機(jī)質(zhì)碳同位素組成全球性變化的一般特征。其中,震旦系燈影組(Z2d)的干酪根碳同位素最輕,其13C值在-31.6‰~-36.7‰之間,平均值為-34.9‰(表1)。

    圖2 四川盆地地層與其干酪根13C值分布綜合圖Fig.2 Generalized stratigraphy column and distribution of stable carbon isotopes for kerogens from various strata in Sichuan Basin

    中生界中、下三疊統(tǒng)地層的干酪根碳同位素進(jìn)一步變重。其中,飛仙關(guān)組(T1f)—嘉陵江組(T1j)的13C值在-24.5‰~-31.0‰,平均值為-26.1‰;雷口坡組(T2l)的13C值在-23.6‰~-26.6‰,平均值為-25.3‰,在海相地層中最重。

    關(guān)于海相沉積有機(jī)質(zhì)碳同位素隨時(shí)代變重的原因目前尚無定論,存在多種觀點(diǎn)。Petersetal.[3]認(rèn)為,在前寒武紀(jì)和下古生代沉積中干酪根碳同位素之所以輕,是因其含有較多的貧13C且抗生物降解能力強(qiáng)的類脂物,而富含13C的碳水化合物及蛋白質(zhì)則在成巖作用中已被降解掉,因這個(gè)時(shí)期的浮游植物僅具有機(jī)包膜,抗降解性低;到了中、新生代時(shí)期,許多浮游植物有了鈣質(zhì)、硅質(zhì)的介殼保護(hù),使碳水化合物不易被降解,因而干酪根中相對富13C。Galimov[5]則認(rèn)為,引起碳同位素隨地質(zhì)時(shí)代發(fā)生全球性變化的原因,是有機(jī)碳與碳酸鹽碳之間的比例及兩者間的同位素分餾值不同所致。另外,還有人將這種現(xiàn)象解釋為是光合作用加劇[16]和大氣CO2濃度變化[20]的結(jié)果。

    2.2 陸相地層干酪根

    從以上敘述可看出,研究區(qū)連同川東南地區(qū)龍?zhí)督M在內(nèi)的陸相地層有機(jī)碳同位素組成具有隨時(shí)代變新而變輕的分布。陸相地層中干酪根的碳同位素組成受多種因素影響,有機(jī)質(zhì)生源及其沉積環(huán)境應(yīng)是主要的因素,但還有一些需進(jìn)一步探索的問題。如須家河組煤層的干酪根碳同位素(平均值-24.3‰)要比龍?zhí)督M煤(-23.4‰)輕不少,而它們的有機(jī)源均主要為裸子植物,成熟度也相近,似乎煤層有機(jī)碳同位素也存在某種時(shí)代效應(yīng)。

    3 烴源巖有機(jī)質(zhì)生源及沉積環(huán)境區(qū)分

    3.1 海相與陸相烴源巖

    據(jù)Galimov[5]及相關(guān)的實(shí)際分析資料,漸新世及之后的陸地有機(jī)質(zhì)碳同位素比海生有機(jī)質(zhì)輕,而漸新世之前的海源有機(jī)碳比陸源碳輕。因而,同層位的海、陸相烴源巖有機(jī)質(zhì)可借鑒這種觀點(diǎn)加以區(qū)分,而不同時(shí)代的海、陸相烴源巖難以適用這種方法。Sofer[6]依據(jù)大量分析數(shù)據(jù),認(rèn)為飽和烴和芳烴組份的碳同位素組成可區(qū)分海相與陸相有機(jī)質(zhì),并統(tǒng)計(jì)回歸出兩者的分界線:δ13C芳=1.14δ13C飽+5.46(圖3);同時(shí)還提出用CV值(CV= —2.53δ13C飽+2.22δ13C芳—11.65)來表征海、陸相有機(jī)質(zhì),兩者界限為0.47。

    如圖3所示,研究區(qū)下侏羅統(tǒng)自流井組和上三疊統(tǒng)須家河組陸相烴源巖的飽和烴和芳烴組份的δ13C數(shù)據(jù)點(diǎn)均落在此圖版的陸相區(qū),相應(yīng)的CV值都高于0.47,變化在1.12~13.29之間,平均值達(dá)4.46,表明這種方法的有效性。而震旦系至下三疊統(tǒng)海相地層的這兩組份δ13C值大多分布在海相區(qū)域,但有部分上古生界二疊系及下三疊統(tǒng)海相烴源巖的樣點(diǎn)落在陸相區(qū)。這些海相烴源巖樣品的CV值總體較低,在-7.60~4.51之間,其中下古生界樣品基本上都低于0.47,總平均值為-0.23,整體上表征海相有機(jī)質(zhì)性質(zhì)。不過,有部分上古生界樣品高于此界限值。

    圖3 四川盆地海、陸相烴源巖飽和烴和芳烴組份13C值分布圖Fig.3 Plot of stable carbon isotope ratios for saturate vs.aromatic hydrocarbons from marine and terrestrial source rocks in Sichuan Basin

    實(shí)際上,上述方法不是區(qū)分烴源巖的沉積環(huán)境,而是有機(jī)質(zhì)來源[3],即利用飽和烴與芳烴組份之間的δ13C差值來識別海洋水生生物與陸源高等植物有機(jī)質(zhì)。由于海、陸相烴源巖分別主要含這兩類有機(jī)質(zhì),因而可用這種方法加以區(qū)分。本研究所分析的部分上古生界烴源巖樣點(diǎn)之所以落在陸相區(qū),CV值高于0.47,可能原因是其有機(jī)源中有部分陸源輸入。

    3.2 龍?zhí)督M海陸交互相烴源巖

    四川盆地上二疊統(tǒng)龍?zhí)督M(吳家坪)組是一套優(yōu)質(zhì)的海陸交互相烴源巖[21],呈區(qū)域性分布。在該地層沉積時(shí)期,盆地具有西南高、東北低的西陸東海的古地理格局,因而其沉積自西向東出現(xiàn)明顯的由陸到海的相變??傮w來說,成都—南充一線以南及華鎣山以西為陸相、海陸交互相沉積區(qū),巖性主要為煤、碳質(zhì)泥巖和泥巖;而盆地東北部為海灣潟湖環(huán)境,沉積深灰色、灰黑色泥巖及泥灰?guī)r[22]。這兩種相帶之外的區(qū)域則為灰泥臺坪、淺水陸棚環(huán)境,主要為沉積海相碳酸鹽巖為主的同時(shí)異相吳家坪組相區(qū),巖石組合主要為灰色泥晶灰?guī)r、微晶灰?guī)r夾泥巖。

    該套烴源巖的干酪根碳同位素比值呈區(qū)域性變化(圖4)。在川東北海灣潟湖相區(qū)的各探井中,該層位泥質(zhì)巖、泥灰?guī)r樣品的碳同位素較輕,其δ13C值變化在-24.0‰~-28.7‰之間,主要集中在-25.5‰~-28.0‰范圍,平均值為-26.7‰。其中,大普光構(gòu)造帶的普光5和毛壩3井樣品的δ13C值普遍較小,平均值最低,分別為-27.5‰和-27.6‰;通南巴地區(qū)的河壩1、金溪1井δ13C平均值次之,均為-27.0‰;其他探井的此平均值稍大些,在-25.9‰~-26.8‰之間。

    圖4 四川盆地龍?zhí)?吳家坪)組不同沉積相帶烴源巖干酪根13C平均值變化圖Fig.4 Carbon isotope variation of kerogens in Longtan(Wujiaping) source rocks from different sedimentary facies in Sichuan Basin

    相比之下,川東南近海湖沼相區(qū)中這套含煤烴源巖樣品的碳同位素顯得很重,其δ13C值變化在-22.3‰~-26.9‰之間,主要在-23.0‰~-24.5‰范圍,平均值為-23.8‰,比川東北地區(qū)重近3‰。其煤和泥巖的干酪根δ13C值分別為-22.3‰~-24.8‰和-22.4‰~-26.9‰,兩者平均值相近,分別為-23.4‰和-23.9‰。區(qū)內(nèi)丁山1井、隆盛1井和南川、習(xí)水良村露頭剖面樣品的δ13C大都較重,平均值分別為-24.3‰、-24.1‰和-23.6‰、-23.4‰。

    地處上述兩個(gè)區(qū)域之間的吳家坪組灰泥臺坪相烴源巖的干酪根碳同位素組成變化較大,所取樣品的δ13C值分布在-23.1‰~-28.2‰之間,主要在-24.0‰~-28.0‰范圍,平均值為-26.6‰,總體上接近于川東北地區(qū)。該區(qū)域的吳家坪組烴源巖碳同位素組成呈明顯的分區(qū)性。地處渝東涪陵地區(qū)的福石1、建深1和三星1井樣品的δ13C平均值較小,分別為-26.7‰、-26.1‰和-27.0‰,與川東北地區(qū)相近;而位于其南邊的焦頁1井樣品的碳同位素較重,δ13C平均值為-23.5‰,與川東南地區(qū)相近。

    龍?zhí)?吳家坪)組不同沉積相帶烴源巖干酪根δ13C值的上述區(qū)域性變化,表征其有機(jī)質(zhì)生源構(gòu)成存在明顯的差異性?;诠派!㈥懺从袡C(jī)質(zhì)碳同位素的組成關(guān)系和實(shí)際沉積相分析,可認(rèn)為川東北地區(qū)海灣潟湖相烴源巖中貧13C的有機(jī)質(zhì)主要生源為海生生物,陸源輸入較少,有機(jī)質(zhì)類型以ⅡA型為主;而川東南地區(qū)近海湖沼相地層中富13C的有機(jī)質(zhì)生源以陸源高等植物為主,類型主要為Ⅲ型;在兩者的過渡相區(qū),生源構(gòu)成變化較大,其東北部以水生生物占優(yōu)勢,而東南部則以陸源輸入為主。它們的芳烴生源和指相參數(shù)佐證了這種認(rèn)識[7]。

    3.3 不同相帶湖相烴源巖

    通常認(rèn)為半深湖—深湖相烴源巖的有機(jī)質(zhì)生源以水生生物為主,碳同位素組成較輕;而濱淺湖相地層中陸源輸入較多,碳同位素較重,但湖相有機(jī)質(zhì)碳同位素組成的影響因素較多。例如,浮游生物勃發(fā)導(dǎo)致水體中CO2供給度不足,致使碳同位素分餾程度降低,造成有機(jī)質(zhì)碳同位素變重。在國內(nèi)外含油氣盆地都出現(xiàn)過深湖相烴源巖中有機(jī)質(zhì)δ13C異常重的例子[23]。因而,對湖相烴源巖有機(jī)質(zhì)進(jìn)行生源分析時(shí),僅依據(jù)碳同位素?cái)?shù)據(jù)可能會得出錯(cuò)誤的認(rèn)識,需結(jié)合其他地球化學(xué)參數(shù)來綜合考慮。下面以元陸6井千佛崖組湖相烴源巖為例,利用碳同位素和其他分析資料來區(qū)分其不同相帶烴源巖的生源構(gòu)成。

    元陸6井地處川東北地區(qū)元壩構(gòu)造帶,其千佛崖組中上部(3 542 m以上)為薄層灰色泥巖與砂巖、粉砂巖互層(圖5),屬濱淺湖相沉積;中部(3 542~3 558 m)為深湖相厚層狀灰黑色泥頁層;下部(3 564~3 596 m)為深灰色泥巖與粉砂巖、粉砂質(zhì)泥巖互層,應(yīng)為淺湖—半深湖相沉積。

    該井所分析的6個(gè)千佛崖組代表性泥巖樣品的干酪根碳同位素組成變化較大,與它們所在的沉積相帶相一致。上部地層2個(gè)濱淺湖相灰色泥巖樣品的干酪根碳同位素較重,δ13C值分別為-24.4‰和-23.3‰(圖5、表2),指示其有機(jī)質(zhì)生源中高等植物占優(yōu)勢,類型為Ⅲ型。中部的2個(gè)灰黑色深湖相泥巖的碳同位素較輕,分別為-26.5‰和-27.0‰,比上部的2個(gè)樣品輕2‰~3‰,表征生源以水生生物為主,類型為Ⅱ型。下部2個(gè)深灰色泥巖樣品的碳同位素組成介于上述兩處樣品之間,分別為-25.8‰和-25.6‰,應(yīng)具有水生和陸源雙重生源。

    這些烴源巖樣品的瀝青A中,飽和烴/芳烴比值所指示的生源意義與上述碳同位素組成相符。淺湖相烴源巖樣品中飽/芳比值低(<0.5),表征有機(jī)質(zhì)類型較差,生源以陸源輸入為主;而半深湖—深湖相泥巖中此比值較高(>2),指示有機(jī)源以水生生物為主。它們的熱脫附輕烴分子標(biāo)志物則直觀佐證了這些烴源巖有機(jī)質(zhì)生源上的差別。上部淺湖相地層樣品的C6~C8輕烴中富含甲苯和甲基環(huán)已烷等源于高等植物有機(jī)質(zhì)的化合物(圖5a),相應(yīng)的正庚烷/甲基環(huán)已烷比值低于0.4(表2);而中部深湖相泥巖中這些陸源化合物相對較少(圖5b),此輕烴比值在1.9左右;下部半深湖相泥巖的輕烴組成(圖5c)及參數(shù)(1.0上下)介于前兩者之間。另外,這幾個(gè)樣品的TOC與δ13C在數(shù)值上呈負(fù)相關(guān)性,表明這些烴源巖的有機(jī)質(zhì)豐度主要與沉積環(huán)境的氧化還原性條件有關(guān)[24-25],取決于相帶的變化。在深湖環(huán)境中,水體深保存條件好,沉積物中TOC值較高(>1.4%),所含類脂物較多,碳同位素較輕(<-26.5‰);而濱淺湖相中,沉積有機(jī)質(zhì)的生物降解作用強(qiáng),TOC值低(<0.75%),類脂物較少,碳同位素較重(>-24.5‰)。

    圖5 元陸6井千佛崖組烴源巖干酪根13C值和相關(guān)地球化學(xué)參數(shù)變化剖面圖Fig.5 Profile of kerogen 13C values and related geochemical parameters of Qianfoya source rocks in Well YL-6

    井號井段/m巖性干酪根δ13C/‰TOC/%飽/芳甲苯/正庚烷正庚烷/甲基環(huán)已烷Ro/%元陸63497~3498灰色泥巖-24.40.500.29.130.391.36元陸63508~3510灰色泥巖-23.30.740.426.700.191.39元陸63546~3549灰黑色泥巖-26.51.412.51.141.921.41元陸63556~3558灰黑色泥巖-27.01.405.50.221.821.44元陸63585~3586深灰色泥巖-25.80.963.20.291.051.46元陸63587~3588深灰色泥巖-25.60.863.50.621.041.47

    3.4 煤與泥巖

    在煤系地層中,由于泥巖和煤的有機(jī)質(zhì)生源相近,均以高等植物為主,因而兩者在總碳同位素組成上相差不多,在高熱演化階段更是如此。如本區(qū)須家河組、龍?zhí)督M的煤和泥巖干酪根δ13C值就很接近。令人關(guān)注的是,它們的正構(gòu)烷烴單體烴碳同位素存在明顯差別,可作為兩者烴源的識別標(biāo)志。

    如圖6所示,大普光構(gòu)造帶雷北1井的須家河組煤和泥巖樣品的正構(gòu)烷烴單體烴δ13C值分布曲線明顯不同。除在C13~C16范圍內(nèi)有所交疊外,在C17~C29范圍兩者明顯分離,且在C17~C25之間隨碳數(shù)增大其差距呈增大的趨勢。在δ13C值上,泥巖樣品比煤樣輕1~4‰。這兩個(gè)樣品深度相近,具有相近的成熟度(Ro在1.1%~1.2%之間),對它們碳同位素組成的差異沒有影響。川西北地區(qū)川科1井須家河組的煤與泥巖樣品之間,正構(gòu)烷烴δ13C值分布曲線有類似差別,同樣是泥巖的δ13C比煤要輕得多,只是隨碳數(shù)的變化稍有不同。從圖6還可觀察到,這兩井的煤與泥巖之間,正構(gòu)烷烴δ13C值分布曲線形態(tài)也明顯不同。煤樣在C13~C22范圍內(nèi)δ13C值變化幅度不大(1‰上下),其曲線形態(tài)呈平緩型;而在高碳數(shù)部分,隨碳數(shù)升高其δ13C顯著變輕,從C22到C28(C29)變輕超過3‰,曲線呈負(fù)向線型。而泥巖樣品中,C25之前的正構(gòu)烷烴δ13C隨碳數(shù)呈逐漸變輕的負(fù)向線型變化,在C25之后則變化不大,呈平緩狀。

    圖6 須家河組煤與泥巖樣品正構(gòu)烷烴系列碳同位素分布曲線對比圖Fig.6 Carbon isotope profiles of n-alkanes for Xujiahe coals and mudstones

    沉積有機(jī)質(zhì)中正構(gòu)烷烴的來源很多,包括陸地高等植物、水生生物和細(xì)菌及其他微生物。一般認(rèn)為,碳數(shù)小于C22的正構(gòu)烷烴主要來自水生藻類,而高碳數(shù)化合物則主要源于高等植物,但細(xì)菌及其他微生物,甚至一些水生植物也能合成長鏈的正構(gòu)化合物[26]。上述須家河組煤與泥巖的正構(gòu)烷烴碳同位素分布不同,表明它們的有機(jī)質(zhì)生源構(gòu)成和沉積環(huán)境性質(zhì)有所不同,泥巖有機(jī)質(zhì)中可能有更多的細(xì)菌貢獻(xiàn)。

    4 成熟度對碳同位素組成的影響

    4.1 干酪根碳同位素

    模擬實(shí)驗(yàn)表明,在生烴演化過程中干酪根的碳同位素組成變重不會超過2‰~3‰[4]。其中,下古生界海相腐泥型(Ⅰ型)干酪根的δ13C值在從未成熟到過成熟乃至淺變質(zhì)演化過程中發(fā)生的變化很小(<1‰)[27]。而蘇愛國[28]和熊永強(qiáng)等[29]的模擬實(shí)驗(yàn)結(jié)果有所不同,從低成熟到過成熟,Ⅲ、Ⅱ和Ⅰ型干酪根的δ13C分別變重0.8‰、1.5‰和2‰,類型好的干酪根卻變化更大。其原因可能是所用的干酪根地質(zhì)年代不同(非下古生代)所致。目前,關(guān)于干酪根的碳同位素組成在天然熱演化過程中的變化狀況報(bào)道較少。由于烴源巖中有機(jī)質(zhì)類型存在非均質(zhì)性而導(dǎo)致碳同位素組成上的差別,因而研究成熟度對干酪根δ13C值的影響,只能選擇有機(jī)質(zhì)生源及其沉積環(huán)境相對一致的地層。為此,本文以須家河組煤巖和茅口組灰?guī)r的干酪根為例來探討這個(gè)問題。

    如圖7所示,須家河組煤樣的Ⅲ型干酪根δ13C值隨熱演化程度增高有較大的變化。在Ro值為1.2%上下的成熟晚期時(shí),其δ13C值在-25.2‰左右;當(dāng)達(dá)到Ro值高于2.0%的過成熟階段時(shí),它們的δ13C值大多重于-24‰,且隨Ro值的進(jìn)一步增高呈持續(xù)變重趨勢,在Ro值為2.4%的1個(gè)樣品中δ13C變重到-23.5‰。在整個(gè)樣品的成熟度范圍(Ro值1.2%~2.45%),須家河組Ⅲ型干酪根的δ13C變重約1.7‰。

    圖7 須家河組煤和茅口組灰?guī)r的干酪根δ13C隨Ro的變化Fig.7 Plot of δ13C vs. Ro of kerogens from Xujiahe coals and Maokou carbonates

    四川盆地茅口組主要為臺地和緩坡相沉積,巖性主要為泥晶灰?guī)r和含泥灰?guī)r,有機(jī)質(zhì)生源主要為水生生物,類型主要為Ⅱ型。因存在一定的相變,所分析樣品的干酪根δ13C有較大的變化范圍,但隨成熟度增高還是呈現(xiàn)出明顯的變重趨勢(圖7)。在Ro(由固體瀝青反射率換算)為1.6%~2.0%的高成熟中晚期階段,δ13C值主要分布在-28.0‰左右;在Ro高于2.0%的過成熟階段時(shí),δ13C明顯變重,大都重于-28‰;當(dāng)Ro增高到2.5%以上時(shí),δ13C重于-27‰。在樣品的Ro值(1.6%~2.9%)范圍內(nèi),茅口組干酪根的δ13C約變重1.6‰。從圖7中的演化趨勢線看,在相同的成熟度范圍內(nèi),茅口組Ⅱ型干酪根的δ13C變重程度似乎稍低于須家河組Ⅲ型干酪根。

    4.2 正構(gòu)烷烴分子碳同位素

    像宏觀有機(jī)組份一樣,單體烴的碳同位素組成也隨成熟度增高而變重。據(jù)文獻(xiàn)[30],烴源巖的正構(gòu)烷烴等化合物在生油窗范圍內(nèi)一般變重2‰~3‰,而在更高熱演化階段時(shí)的變化未見有報(bào)道。本文利用不同成熟度的須家河組泥巖樣品分析資料,闡明正構(gòu)烷烴碳同位素在高—過成熟階段的變化。

    取自大普光地區(qū)雷北1井和元壩地區(qū)元壩204井須家河組Ro值在1.2%上下的泥巖樣品中,正構(gòu)烷烴的δ13C隨碳數(shù)增高而變輕,從C14的-28‰左右逐漸減小到C29的-34‰上下,整個(gè)系列的變化范圍在6‰左右(圖8)。正構(gòu)烷烴碳同位素隨碳數(shù)增加而變輕的分布表征陸源高等植物有機(jī)質(zhì)生源[31-32],與須家河組含煤烴源巖的有機(jī)源相一致。對于正構(gòu)烷烴δ13C隨碳數(shù)而變輕的分布模式,Murrayetal.[31]認(rèn)為最可能的原因是細(xì)菌和真菌對高等植物有機(jī)質(zhì)的改造作用所致;而有的學(xué)者則認(rèn)為是因不同碳數(shù)的正構(gòu)烷烴在樹葉生長周期內(nèi)形成的比例不同所造成,因植物在不同生長期的水分、營養(yǎng)條件和光照強(qiáng)度不同[33-34]。

    圖8 須家河組不同成熟度泥巖的正構(gòu)烷烴δ13C變化曲線Fig.8 Carbon isotope profiles of n-alkanes in Xujiahe mudstones with different maturities

    令人關(guān)注的是,本區(qū)高—過成熟的須家河組泥巖樣品中,其正構(gòu)烷烴的碳同位素分布完全不同于上述成熟晚期的樣品。如圖8所示,元壩6等井Ro值在2.0%上下的泥巖中,正構(gòu)烷烴單體烴碳同位素的分布曲線總體上呈平直狀,各碳數(shù)化合物的δ13C值主要集中在-30‰~-32‰范圍。與前幾個(gè)成熟度相對較低的樣品相比,它們的中、低碳數(shù)(C15~C23)正構(gòu)烷烴的碳同位素變輕,而高碳數(shù)(C24~C32)化合物變重。這些成熟度不同的烴源巖中正構(gòu)烷烴碳同位素出現(xiàn)這樣的分布變化,可能與在高—過成熟階段長鏈化合物發(fā)生的裂解作用有關(guān)。有研究表明,在有機(jī)質(zhì)熱降解過程中,分子量低的產(chǎn)物有較多的12C[35]。據(jù)此及動(dòng)力學(xué)分餾原理不難推斷:當(dāng)正構(gòu)烷烴在較高地溫條件下發(fā)生裂解時(shí),碳數(shù)高的化合物裂解作用較強(qiáng)烈,失去較多碳同位素輕的分子,導(dǎo)致其殘余部分的碳同位素變重;而低碳數(shù)化合物的裂解強(qiáng)度較低,且有大量由高碳數(shù)化合物裂解而來的富12C化合物的加入,使得其整體碳同位素變輕。這種現(xiàn)象在油(瀝青)—源巖對比時(shí)要加以考慮。

    5 結(jié)論

    (1) 四川盆地震旦系燈影組至中三疊統(tǒng)雷口坡組海相烴源巖的干酪根碳同位素組成隨層位變新呈逐漸變重的趨勢,海洋生物的演化可能是其重要的影響因素;而包括川東南上二疊統(tǒng)龍?zhí)督M在內(nèi)的陸相地層有機(jī)質(zhì)碳同位素組成則具有反向的年代變化,主要與其生源構(gòu)成及其沉積環(huán)境性質(zhì)不同有關(guān)。

    (3) 從成熟晚期到高—過成熟階段,海相Ⅱ型和煤系Ⅲ型干酪根的碳同位素均持續(xù)變重,13C值升高1‰~2‰;煤系泥巖中正構(gòu)烷烴系列的碳同位素分布曲線,由負(fù)向線型分布模式向平直型演變,因而在油氣與烴源巖對比時(shí)要應(yīng)用成熟度相近的樣品。

    )

    [1] 胡見義,黃第藩. 中國陸相石油地質(zhì)理論基礎(chǔ)[M]. 北京:石油工業(yè)出版社,1991:189. [Hu Jianyi, Huang Difan. The bases of nonmarine petroleum geology in China[M]. Beijing: Petroleum Industry Press, 1991: 189.]

    [2] 梁狄剛,郭彤樓,陳建平,等. 中國南方海相生烴成藏研究的若干新進(jìn)展(二):南方四套區(qū)域性海相烴源巖的地球化學(xué)特征[J]. 海相油氣地質(zhì),2009,14(1):1-15. [Liang Digang, Guo Tonglou, Chen Jianping, et al. Some progresses on studies of hydrocarbon generation and accumulation in marine sedimentary regions, southern China (part 2): Geochemical characteristics of four suits of regional marine source rocks, South China[J]. Marine Oil and Gas Geology, 2009, 14(1): 1-15.]

    [3] Peters K E, Walters C C, Moldowan J M. The biomarker guide: I. biomarkers and isotopes in the environment and human history[M]. 2nd Ed. Cambridge: Cambridge University Press, 2005: 136-156.

    [4] Sackett W M. Determination of kerogen maturity by the pyrolysis-carbon isotope method[J]. Organic Geochemistry, 1984, 6: 359-363.

    [5] Galimov E M. Isotope organic geochemistry[J]. Organic Geochemistry, 2006, 37(10): 1200-1262.

    [6] Sofer Z. Stable carbon isotope compositions of crude oils: application to source depositional environments and petroleum alteration[J]. AAPG Bulletin, 1984, 68(1): 31-49.

    [7] 朱揚(yáng)明,顧圣嘯,李穎,等. 四川盆地龍?zhí)督M高熱演化烴源巖有機(jī)質(zhì)生源及沉積環(huán)境探討[J]. 地球化學(xué),2012,41(1):35-44. [Zhu Yangming, Gu Shengxiao, Li Ying, et al. Biological organic source and depositional environment of over-mature source rocks of Longtan Formation in Sichuan Basin[J]. Geochemica, 2012, 41(1): 35-44.]

    [8] 黃籍中. 中上揚(yáng)子區(qū)海相沉積烴源研究(之二)[J]. 天然氣勘探與開發(fā),2000,23(1):9-27. [Huang Jizhong. Study on hydrocarbon source of marine sediments in Middle-Upper Yangtze platform (part 2)[J]. Natural Gas Exploration and Development, 2000, 23(1): 9-27.]

    [9] 黃籍中,冉隆輝. 四川盆地震旦系燈影灰?guī)r黑色瀝青與油氣勘探[J],石油學(xué)報(bào),1989,10(1):27-36. [Huang Jizhong, Ran Longhui. Bitumen and oil-gas exploration in Sinian “Denying limestone” in Sichuan Basin[J]. Acta Petrolei Sinica, 10(1): 27-36.]

    [10] 王順玉,戴鴻鳴,王海清,等. 大巴山、米倉山南緣烴源巖特征研究[J]. 天然氣地球科學(xué),2000,11(4/5):4-16. [Wang Shunyu, Dai Hongming, Wang Haiqing, et al. Source rock feature of the south of the Dabashan and Mi-cangshan area[J]. Natural Gas Geoscience, 2000, 11(4/5): 4-16.]

    [11] 童崇光. 四川盆地油氣地質(zhì)研究[J]. 西安地質(zhì)學(xué)院學(xué)報(bào),1990,12(3):9-17. [Tong Chongguang. New progress of the oil-gas geology research in the Sichuan Basin[J]. Journal of Xi’an College of Geology, 1990, 12(3): 9-17.]

    [12] 四川油氣區(qū)石油地質(zhì)志編寫組. 中國石油地質(zhì)志(卷十) 四川油氣區(qū)[M]. 北京:石油工業(yè)出版社,1989:299. [Editorial Group of Sichuan Oil and Gas Field. Petroleum geology of China (Vol. 10) Sichuan oil and gas field[M]. Beijing: Petroleum Industry Press, 1989: 299.]

    [13] 黃籍中,陳盛吉,宋家榮,等. 四川盆地?zé)N源體系與大中型氣田形成[J]. 中國科學(xué)D輯:地球科學(xué),1996,26(6):504-510. [Huang Jizhou, Chen Shengji, Song Jiarong, et al. Hydrocarbon source systems and formation of gas fields in Sichuan Basin[J]. Science in China Series D: Earth Sciences, 1997, 40(1): 32-42.]

    [14] 黃仁春. 川東北元壩地區(qū)雷口坡組天然氣來源與成藏分析[J]. 現(xiàn)代地質(zhì),2014,28(2):412-418. [Huang Renchun. Source and accumulation of natural gas in Leikoupo Formation, Yuanba area, eastern-northern Sichuan Basin[J]. Geoscience, 2014, 28(2): 412-418.]

    [15] 陳盛吉,萬茂霞,杜敏,等. 川中地區(qū)侏羅系油氣源對比及烴源條件研究[J]. 天然氣勘探與開發(fā),2005,28(2):11-14. [Chen Shengji, Wan Maoxia, Du Min, et al. Jurassic oil-gas source correlation and hydrocarbon source condition in central Sichuan region[J]. Natural Gas Exploration & Development, 2005, 28(2): 11-14.]

    [16] Welte D H, Kalkreuth W. Age-trend in carbon isotopic composition in Paleozoic sediments[J]. Naturwissenschaften, 1975, 62(10): 482-483.

    [17] Galimov E M.13C/12C in kerogen[M]//Durand B. Kerogen-Insoluble Organic Matter from Sedimentary Rocks. Paris: Editions Technip, 1980: 271-299.

    [18] Lewan M D. Stable carbon isotopes of amorphous kerogens from Phanerozoic sedimentary rocks[J]. Geochimica et Cosmochimica Acta, 1986, 50(8): 1583-1591.

    [19] 張大江,王大銳,章聆. 我國不同地質(zhì)時(shí)期沉積有機(jī)質(zhì)中穩(wěn)定碳同位素組成特征及其意義[J]. 科學(xué)通報(bào),1991,36(4):302-305. [Zhang Dajiang, Wang Darui, Zhang Ling. Characteristics of stable isotopic composition of sedimentary organic matter in different geological ages, China [J]. Chinese Science Bulletin, 1991, 36(4): 302-305.]

    [20] Kump L R, Arthur M A. Interpreting carbon-isotope excursions: carbonates and organic matter[J]. Chemical Geology, 1999, 161(1/2/3): 181-198.

    [21] 梁狄剛,郭彤樓,陳建平,等. 中國南方海相生烴成藏研究的若干新進(jìn)展(一):南方四套區(qū)域性海相烴源巖分布[J]. 海相油氣地質(zhì),2008,13(2):1-16. [Liang Digang, Guo Tonglou, Chen Jianping, et al. Some progresses on studies of hydrocarbon generation and accumulation in marine sedimentary regions, southern China (part 1): Distribution of four suits of regional marine source rocks[J]. Marine Origin Petroleum Geology, 2008, 13(2): 1-16.]

    [22] 梁狄剛,郭彤樓,邊立曾,等. 中國南方海相生烴成藏研究的若干新進(jìn)展(三):南方四套區(qū)域性海相烴源巖的沉積相及發(fā)育的控制因素[J]. 海相油氣地質(zhì),2009,14(2):1-19. [Liang Digang, Guo Tonglou, Bian Lizeng, et al. Some progresses on studies of hydrocarbon generation and accumulation in marine sedimentary regions, southern China (part 3): Controlling factors on the sedimentary facies and development of Paleozoic marine source rocks[J]. Marine Oil and Gas Geology, 2009, 14(2): 1-19.]

    [23] 孫玉梅,李友川,黃正吉. 部分近海湖相烴源巖有機(jī)質(zhì)異常碳同位素組成[J]. 石油勘探與開發(fā),2009,36(5):609-616. [Sun Yumei, Li Youchuan, Huang Zhengji. Abnormal carbon isotopic compositions in organic matter of lacustrine source rocks close to sea[J]. Petroleum Exploration and Development, 2009, 36(5): 609-616.]

    [24] Freeman K H, Hayes J M, Trendel J M, et al. Evidence from carbon isotope measurements for diverse origins of sedimentary hydrocarbons[J]. 1990, Nature, 343(6255): 254-256.

    [25] Harris N B, Freeman K H, Pancost R D, et al. The character and origin of lacustrine source rocks in the Lower Cretaceous synrift section, Congo Basin, west Africa[J]. AAPG Bulletin, 2004, 88(8): 1163-1184.

    [26] Chikaraishi Y, Naraoka H.13C andD relationships among threen-alkyl compound classes (n-alkanoic acid,n-alkane andn-alkanol) of terrestrial higher plants[J]. Organic Geochemistry, 2007, 38(2): 198-215.

    [27] Buchardt B, Clausen J, Thomsen E. Carbon isotope composition of Lower Palaeozoic kerogen: effects of maturation[J]. Organic Geochemistry, 1986, 10(1/2/3): 127-134.

    [28] 蘇艾國. 干酪根碳同位素在成熟和風(fēng)化過程中變化規(guī)律初探[J]. 礦物巖石地球化學(xué)通報(bào),1999,18(2):79-84. [Su Aiguo. Discussion on variation of carbon isotope of kerogen during thermal maturation and weathering[J]. Bulletin of Mineralogy, Petrology and Geochemistry, 1999, 18(2): 79-84.]

    [29] 熊永強(qiáng),張海祖,耿安松. 熱演化過程中干酪根碳同位素組成的變化[J]. 石油實(shí)驗(yàn)地質(zhì),2004,26(5):484-487. [Xiong Yongqing, Zhang Haizu, Geng Ansong. Variation of carbon isotopic composition of kerogen during thermal evolution[J]. Petroleum Geology and Experiment, 2004, 26(5): 484-487.]

    [30] Clayton C J, Bjor?y M. Effect of maturity on13C/12C ratios of individual compounds in North Sea oils[J]. Organic Geochemistry, 1994, 21(6/7): 737-750.

    [31] Murray A P, Summons R E, Boreham C J, et al. Biomarker and n-alkane isotope profiles for Tertiary oils: relationship to source rock depositional setting[J]. Organic Geochemistry, 1994, 22(3/4/5): 521-542.

    [32] Schouten S, Hartgers W A, Lòpez J F, et al. A molecular isotopic study of13C-enriched organic matter in evaporitic deposits: recognition of CO2-limited ecosystems[J]. Organic Geochemistry, 2001, 32(2): 277-286.

    [33] Badewien T, Vogts A, Rullk?tter J.n-Alkane distribution and carbon stable isotope composition in leaf waxes of C3and C4plants from Angola[J]. Organic Geochemistry, 2015, 89-90: 71-79.

    [34] Tu T T N, Derenne S, Largeau C, et al. Diagenesis effects on specific carbon isotope composition of plantn-alkanes[J]. Organic Geochemistry, 2004, 35(3): 317-329.

    [35] Marais D J D, Donchin J H, Nehring N L, et al. Molecular carbon isotopic evidence for the origin of geothermal hydrocarbons[J]. Nature, 1981, 292(5826): 826-828.

    CompositionalVariationsandGeochemicalSignificancesofStableCarbonIsotopeforOrganicMattersfromMarineandTerrestrialSourceRocksinSichuanBasin

    ZHU YangMing1, LI Ying1, HAO Fang2, ZOU HuaYao3, GUO XuSheng4

    1.SchoolofEarthSciences,ZhejiangUniversity,Hangzhou310027,China2.ChinaUniversityofGeosciences,Wuhan430074,China3.CollegeofEarthSciences,ChinaUniversityofPetroleum,Beijing102249,China4.SouthernExplorationCompany,SINOPEC,Chengdu610000,China

    Over 500 kerogens as well selected saturate and aromatic fractions and n-alkanes from various age source rocks in the Sichuan basin were analyzed using MS and GC-ir-MS for13C ratios, to characterize the age-trend in organic carbon isotopic composition in marine and terrestrial sediments and their thermal evolution during over-mature stage. Moreover, different kind of marine and terrestrial source rocks with respect to their biological organic sources were distinguished using carbon isotopes in combination with other analytical data. The results indicate the kerogens from marine sediments in Sinian Dengying to middle Triassic Leikoupo Formations in this basin show a trend toward isotopically heavy values with decreasing age, possibly due to biological evolution. Whereas, a reversal isotopic-age trend for the terrigenous organic carbon is observed in upper Triassic Xujiahe to middle Jurassic Qingfoyan Formations and is considered to be dependent of biological source and sedimentary environment. The organic inputs of marine and terrestrial source rocks cannot be differentiated using the carbon isotopic ratios of kerogen, however, it can be distinguished by13C values of saturated and aromatic fractions and CV values. The marine and lacustrine source rocks with different facies are of obviously varying13C values of kerogens, thus the isotopic ratio, combined with additional relate-source data, can be acted as an indicator for their organic source. Coal and mudstone in coal-bearing strata is not distinguishable for their kerogen carbon isotopes, but has different profiles of n-alkane isotopes with more negative values in mudstone related to coal. The marine (type-Ⅱ)and terrigenous (type-Ⅲ) kerogens both continuously become isotopically enriched in13C by 1‰~2‰ with maturation during high-over maturity. The negatively sloping curve of carbon isotope for individual n-alkanes from coaly source rock at mature stage is transferred into a flat one at high maturity. These variations in carbon isotope compositions are helpful for oil & gas-rock correlation and hydrocarbon source identification.

    source rock; carbon isotope; kerogen; n-alkane; Sichuan Basin

    1000-0550(2017)06-1254-11

    10.14027/j.cnki.cjxb.2017.06.016

    2017-03-07;收修改稿日期2017-05-03

    國家科技重大專項(xiàng)(2011ZX05005-03-009HZ);教育部高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金博導(dǎo)類資助課題(20130101110051)[FoundationNational Science and Technology Major Project, No. 2011ZX05005-03-009HZ; Specialized Research Fund for the Doctoral Program of Higher Education, No. 20130101110051]

    朱揚(yáng)明,男,1954年出生,教授,博士,油氣地球化學(xué),E-mail: zyming@zju.edu.cn

    TE112.11

    A

    亚洲av电影在线进入| 精品福利观看| 精品国产乱码久久久久久男人| 色尼玛亚洲综合影院| 久久天堂一区二区三区四区| 亚洲男人的天堂狠狠| 久久久水蜜桃国产精品网| 日本撒尿小便嘘嘘汇集6| av中文乱码字幕在线| 精品欧美国产一区二区三| 成人三级黄色视频| 日日爽夜夜爽网站| 69av精品久久久久久| 99久久久亚洲精品蜜臀av| 成人亚洲精品av一区二区| 国产精品免费一区二区三区在线| 国产日本99.免费观看| 十八禁网站免费在线| 久久人妻福利社区极品人妻图片| 久久精品国产99精品国产亚洲性色| 美女免费视频网站| 国产成人一区二区三区免费视频网站| 精品国产超薄肉色丝袜足j| 亚洲av电影不卡..在线观看| 最近最新中文字幕大全免费视频| 亚洲成av人片在线播放无| 国产成人精品久久二区二区91| 国产单亲对白刺激| 亚洲成人精品中文字幕电影| 两个人看的免费小视频| 国内揄拍国产精品人妻在线| 国产精品野战在线观看| 午夜免费激情av| 国产精品久久久人人做人人爽| 中文资源天堂在线| 制服诱惑二区| 亚洲乱码一区二区免费版| 一本精品99久久精品77| 听说在线观看完整版免费高清| 国产精品免费视频内射| 在线观看免费午夜福利视频| 亚洲av成人不卡在线观看播放网| 欧美性猛交黑人性爽| 啦啦啦观看免费观看视频高清| 久久久水蜜桃国产精品网| 特大巨黑吊av在线直播| 亚洲欧美精品综合一区二区三区| 亚洲精品在线观看二区| 丰满的人妻完整版| 日本免费a在线| 别揉我奶头~嗯~啊~动态视频| 中文资源天堂在线| av欧美777| 欧美成人性av电影在线观看| 久久热在线av| 亚洲真实伦在线观看| 日本撒尿小便嘘嘘汇集6| 在线观看一区二区三区| netflix在线观看网站| 热99re8久久精品国产| 久久久国产精品麻豆| 三级国产精品欧美在线观看 | 丁香六月欧美| 久久久久久国产a免费观看| 欧美成人午夜精品| 正在播放国产对白刺激| 欧美日韩中文字幕国产精品一区二区三区| 欧美不卡视频在线免费观看 | 国产精品 国内视频| 国产精品久久久久久久电影 | 很黄的视频免费| 国产精品久久久久久精品电影| 亚洲专区国产一区二区| 丰满人妻熟妇乱又伦精品不卡| 桃红色精品国产亚洲av| 精品人妻1区二区| 成人18禁在线播放| 久久久久久国产a免费观看| 三级国产精品欧美在线观看 | 亚洲男人天堂网一区| 老司机深夜福利视频在线观看| 亚洲人成电影免费在线| 久久精品夜夜夜夜夜久久蜜豆 | 免费在线观看视频国产中文字幕亚洲| cao死你这个sao货| 日本成人三级电影网站| 黄色丝袜av网址大全| 国产亚洲精品久久久久5区| 久久精品影院6| 又粗又爽又猛毛片免费看| 操出白浆在线播放| АⅤ资源中文在线天堂| 又粗又爽又猛毛片免费看| 又黄又粗又硬又大视频| 国产成人啪精品午夜网站| 久久精品91无色码中文字幕| 亚洲自拍偷在线| 国产99久久九九免费精品| 精品国产乱码久久久久久男人| 亚洲精品一区av在线观看| 成人永久免费在线观看视频| 18禁黄网站禁片午夜丰满| 久久欧美精品欧美久久欧美| 久久久久免费精品人妻一区二区| 国产成人aa在线观看| 女生性感内裤真人,穿戴方法视频| 小说图片视频综合网站| 丁香六月欧美| 欧美人与性动交α欧美精品济南到| 人人妻,人人澡人人爽秒播| 91麻豆av在线| 亚洲欧美一区二区三区黑人| 亚洲美女黄片视频| 久久久久久久久中文| 免费在线观看亚洲国产| 国模一区二区三区四区视频 | 女警被强在线播放| 亚洲av成人av| 免费在线观看影片大全网站| 在线观看舔阴道视频| 国产97色在线日韩免费| 91成年电影在线观看| 午夜激情av网站| 日韩欧美国产在线观看| 亚洲色图 男人天堂 中文字幕| 日韩精品青青久久久久久| 国产高清videossex| 免费看十八禁软件| 国产91精品成人一区二区三区| 亚洲色图av天堂| 欧美成人免费av一区二区三区| 国产精品 国内视频| 精品乱码久久久久久99久播| av福利片在线观看| 中文字幕最新亚洲高清| 一进一出抽搐动态| 哪里可以看免费的av片| 波多野结衣高清无吗| 日韩欧美在线乱码| 亚洲欧美精品综合久久99| 在线观看美女被高潮喷水网站 | 一进一出抽搐gif免费好疼| 亚洲第一电影网av| 丝袜人妻中文字幕| 国产三级黄色录像| 亚洲,欧美精品.| 久久欧美精品欧美久久欧美| 又黄又爽又免费观看的视频| 人妻丰满熟妇av一区二区三区| 欧美极品一区二区三区四区| 免费在线观看黄色视频的| 美女黄网站色视频| 久久人人精品亚洲av| 少妇粗大呻吟视频| 国产精品一区二区三区四区免费观看 | 午夜福利成人在线免费观看| 欧美丝袜亚洲另类 | 妹子高潮喷水视频| 欧美人与性动交α欧美精品济南到| 天堂av国产一区二区熟女人妻 | 国产成+人综合+亚洲专区| 99久久精品热视频| 免费无遮挡裸体视频| 伦理电影免费视频| 日本撒尿小便嘘嘘汇集6| 最近在线观看免费完整版| 久久天躁狠狠躁夜夜2o2o| 国产亚洲欧美在线一区二区| 18禁国产床啪视频网站| 97人妻精品一区二区三区麻豆| 在线a可以看的网站| 看黄色毛片网站| 国产在线精品亚洲第一网站| 久久久久久亚洲精品国产蜜桃av| 久9热在线精品视频| 啪啪无遮挡十八禁网站| 淫秽高清视频在线观看| 久久香蕉国产精品| 一二三四社区在线视频社区8| 少妇被粗大的猛进出69影院| 国产av在哪里看| 色哟哟哟哟哟哟| 久久国产精品影院| 国产成人精品久久二区二区免费| 久久人妻av系列| 国产69精品久久久久777片 | 一区福利在线观看| 久久婷婷成人综合色麻豆| 变态另类成人亚洲欧美熟女| 免费在线观看视频国产中文字幕亚洲| 国产探花在线观看一区二区| 成年人黄色毛片网站| 黄片大片在线免费观看| 亚洲av日韩精品久久久久久密| 国产精品av视频在线免费观看| 18禁黄网站禁片免费观看直播| 每晚都被弄得嗷嗷叫到高潮| 麻豆成人av在线观看| 国产熟女xx| 真人做人爱边吃奶动态| 欧美久久黑人一区二区| 精品午夜福利视频在线观看一区| 国产视频内射| 香蕉国产在线看| 99riav亚洲国产免费| 久久人妻福利社区极品人妻图片| 亚洲激情在线av| 久久欧美精品欧美久久欧美| 久久国产精品影院| 精品福利观看| 亚洲一区二区三区不卡视频| 两性午夜刺激爽爽歪歪视频在线观看 | 天天添夜夜摸| 国产精品免费一区二区三区在线| 精品电影一区二区在线| 国产午夜精品论理片| 欧美一级a爱片免费观看看 | 色av中文字幕| 操出白浆在线播放| 中国美女看黄片| 欧美成人性av电影在线观看| 一个人免费在线观看的高清视频| 国产精品一区二区三区四区久久| 一区福利在线观看| 免费看a级黄色片| 成人国产一区最新在线观看| 在线观看免费视频日本深夜| 黑人操中国人逼视频| 国产又黄又爽又无遮挡在线| 成人三级做爰电影| 亚洲国产欧美一区二区综合| 亚洲熟妇中文字幕五十中出| 日韩精品免费视频一区二区三区| 亚洲乱码一区二区免费版| 国产av在哪里看| 久久婷婷成人综合色麻豆| 不卡一级毛片| 婷婷精品国产亚洲av| 757午夜福利合集在线观看| 99久久99久久久精品蜜桃| 亚洲av成人精品一区久久| 久久久久免费精品人妻一区二区| 亚洲avbb在线观看| 精品欧美一区二区三区在线| 午夜免费激情av| 亚洲欧美日韩高清专用| 后天国语完整版免费观看| 最近最新免费中文字幕在线| 精品久久久久久久末码| 中亚洲国语对白在线视频| 视频区欧美日本亚洲| 琪琪午夜伦伦电影理论片6080| 久久精品aⅴ一区二区三区四区| 国产野战对白在线观看| 欧美黑人精品巨大| 琪琪午夜伦伦电影理论片6080| 亚洲九九香蕉| 一个人免费在线观看的高清视频| 日本免费a在线| 黄频高清免费视频| 欧美性猛交黑人性爽| 我要搜黄色片| 狠狠狠狠99中文字幕| 淫秽高清视频在线观看| 美女黄网站色视频| 成人三级黄色视频| www日本黄色视频网| 人人妻人人澡欧美一区二区| 亚洲真实伦在线观看| 欧美+亚洲+日韩+国产| 黄色 视频免费看| 欧美三级亚洲精品| 999精品在线视频| 国产爱豆传媒在线观看 | 91麻豆av在线| 亚洲av熟女| 国产精品免费一区二区三区在线| 国语自产精品视频在线第100页| 中文资源天堂在线| 观看免费一级毛片| 午夜老司机福利片| 亚洲人成电影免费在线| 一级作爱视频免费观看| 别揉我奶头~嗯~啊~动态视频| 国产成人aa在线观看| 18禁黄网站禁片午夜丰满| www.www免费av| 黄片大片在线免费观看| 久久久国产成人精品二区| 精品欧美一区二区三区在线| 国产v大片淫在线免费观看| 可以在线观看的亚洲视频| 亚洲成人久久爱视频| 亚洲av第一区精品v没综合| 深夜精品福利| 免费av毛片视频| 免费在线观看完整版高清| 男女做爰动态图高潮gif福利片| 亚洲真实伦在线观看| 人人妻,人人澡人人爽秒播| 啦啦啦免费观看视频1| 亚洲精品在线观看二区| 欧美乱色亚洲激情| 中文资源天堂在线| 国产激情偷乱视频一区二区| 国内精品久久久久精免费| 国产精品一区二区精品视频观看| 精品人妻1区二区| 午夜福利免费观看在线| 最新在线观看一区二区三区| 9191精品国产免费久久| 亚洲中文字幕日韩| 成人国产一区最新在线观看| 免费在线观看完整版高清| 国产精品九九99| 久久久久久大精品| 久久精品国产亚洲av高清一级| 日本a在线网址| 久久亚洲真实| 大型av网站在线播放| 十八禁人妻一区二区| 久久久国产成人免费| 国产视频内射| 国产麻豆成人av免费视频| 色播亚洲综合网| 无人区码免费观看不卡| 免费在线观看黄色视频的| 女人被狂操c到高潮| 国产亚洲欧美98| 久久久久久久久免费视频了| 日本熟妇午夜| 日本熟妇午夜| 亚洲国产日韩欧美精品在线观看 | 一a级毛片在线观看| 日韩国内少妇激情av| 久热爱精品视频在线9| 小说图片视频综合网站| 欧美精品啪啪一区二区三区| 亚洲av五月六月丁香网| 亚洲精品在线美女| 亚洲色图av天堂| 日韩国内少妇激情av| 国产亚洲精品一区二区www| 日韩大尺度精品在线看网址| 日本一本二区三区精品| 日本一本二区三区精品| 国产视频内射| 亚洲,欧美精品.| 精品国产美女av久久久久小说| 1024香蕉在线观看| 日韩欧美三级三区| 精品久久久久久久毛片微露脸| 亚洲自拍偷在线| 日韩欧美 国产精品| 老司机深夜福利视频在线观看| 久久久久免费精品人妻一区二区| 久久 成人 亚洲| 日日爽夜夜爽网站| 性欧美人与动物交配| 一卡2卡三卡四卡精品乱码亚洲| 中出人妻视频一区二区| 欧美日韩亚洲国产一区二区在线观看| 久久久久国产一级毛片高清牌| 日韩欧美在线二视频| 观看免费一级毛片| 日本熟妇午夜| 日本黄大片高清| 婷婷六月久久综合丁香| 90打野战视频偷拍视频| 国产精品久久久久久人妻精品电影| 亚洲色图 男人天堂 中文字幕| 在线看三级毛片| 精品第一国产精品| 超碰成人久久| 又粗又爽又猛毛片免费看| 精品久久久久久,| 久久精品国产亚洲av香蕉五月| 欧美黄色片欧美黄色片| 天堂√8在线中文| 俄罗斯特黄特色一大片| 日本免费a在线| 日韩精品免费视频一区二区三区| 亚洲黑人精品在线| 亚洲精品av麻豆狂野| 三级国产精品欧美在线观看 | 日韩免费av在线播放| 午夜福利18| 国产av一区二区精品久久| 波多野结衣高清作品| 久久久久久国产a免费观看| 国产精品久久视频播放| 三级毛片av免费| 日韩欧美在线二视频| 亚洲成av人片在线播放无| 动漫黄色视频在线观看| 亚洲av成人不卡在线观看播放网| 精品国产美女av久久久久小说| 淫秽高清视频在线观看| 麻豆国产av国片精品| 三级毛片av免费| 中文字幕最新亚洲高清| 亚洲最大成人中文| 亚洲精品中文字幕在线视频| av在线播放免费不卡| 国产精品久久久久久精品电影| 欧美日韩亚洲综合一区二区三区_| 日本一二三区视频观看| 特大巨黑吊av在线直播| 看片在线看免费视频| 一区二区三区激情视频| 人人妻人人澡欧美一区二区| 国产精品影院久久| av国产免费在线观看| 国产精品香港三级国产av潘金莲| 黄色视频不卡| 丝袜人妻中文字幕| 免费搜索国产男女视频| 亚洲熟妇熟女久久| 麻豆av在线久日| 欧美一级毛片孕妇| 欧美绝顶高潮抽搐喷水| 精品国产亚洲在线| 嫩草影视91久久| avwww免费| 欧美日韩亚洲综合一区二区三区_| 欧美成人一区二区免费高清观看 | 亚洲第一欧美日韩一区二区三区| 国产精品久久久久久久电影 | 麻豆av在线久日| 一本一本综合久久| 成人欧美大片| 色综合婷婷激情| 国产精品一区二区三区四区免费观看 | 亚洲第一欧美日韩一区二区三区| 一级毛片女人18水好多| 国产黄片美女视频| 日韩欧美精品v在线| 欧美在线一区亚洲| 一本一本综合久久| 中文字幕久久专区| 国产三级在线视频| 亚洲av熟女| 人妻久久中文字幕网| 亚洲最大成人中文| 97碰自拍视频| 男女之事视频高清在线观看| 精品久久久久久久毛片微露脸| 哪里可以看免费的av片| 男女视频在线观看网站免费 | 搡老熟女国产l中国老女人| 亚洲av美国av| 久久人妻福利社区极品人妻图片| 精品欧美一区二区三区在线| 窝窝影院91人妻| 午夜福利在线在线| 亚洲欧美一区二区三区黑人| 亚洲avbb在线观看| 特级一级黄色大片| 男人舔女人的私密视频| 久久精品国产清高在天天线| 后天国语完整版免费观看| 色尼玛亚洲综合影院| 久久 成人 亚洲| 97人妻精品一区二区三区麻豆| 国产精品久久久av美女十八| 亚洲av日韩精品久久久久久密| 亚洲色图 男人天堂 中文字幕| 亚洲av成人精品一区久久| 在线观看www视频免费| 精品乱码久久久久久99久播| 很黄的视频免费| 亚洲午夜精品一区,二区,三区| 99热只有精品国产| 亚洲七黄色美女视频| 欧美3d第一页| 午夜激情福利司机影院| 老司机深夜福利视频在线观看| 波多野结衣高清无吗| 日本精品一区二区三区蜜桃| 18美女黄网站色大片免费观看| 99在线视频只有这里精品首页| 国产高清videossex| 91大片在线观看| 免费在线观看视频国产中文字幕亚洲| 黄色女人牲交| 麻豆国产av国片精品| 超碰成人久久| 国产精品久久久av美女十八| 91成年电影在线观看| 中文字幕最新亚洲高清| 日本免费一区二区三区高清不卡| 国产伦人伦偷精品视频| 一级作爱视频免费观看| 可以在线观看毛片的网站| 亚洲自拍偷在线| 国产伦人伦偷精品视频| 久久久久精品国产欧美久久久| 成人国产一区最新在线观看| 日韩三级视频一区二区三区| 成人手机av| 欧美不卡视频在线免费观看 | 露出奶头的视频| 法律面前人人平等表现在哪些方面| 欧美性长视频在线观看| 日韩欧美精品v在线| 国产精品精品国产色婷婷| 精品久久久久久久久久久久久| 在线a可以看的网站| 88av欧美| 精品久久蜜臀av无| 欧美绝顶高潮抽搐喷水| 国产精品一及| 国产v大片淫在线免费观看| av福利片在线观看| 亚洲九九香蕉| 高清在线国产一区| 少妇粗大呻吟视频| а√天堂www在线а√下载| 性欧美人与动物交配| 精品少妇一区二区三区视频日本电影| av在线天堂中文字幕| 精品久久久久久久久久久久久| 成人av一区二区三区在线看| 久久天躁狠狠躁夜夜2o2o| 亚洲天堂国产精品一区在线| 老汉色∧v一级毛片| 一个人免费在线观看电影 | 色综合站精品国产| 天堂√8在线中文| 特大巨黑吊av在线直播| 老汉色∧v一级毛片| 亚洲九九香蕉| 色尼玛亚洲综合影院| 美女大奶头视频| 久久久国产成人免费| 变态另类成人亚洲欧美熟女| 一级片免费观看大全| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品av视频在线免费观看| 一级片免费观看大全| 欧美极品一区二区三区四区| 伊人久久大香线蕉亚洲五| 十八禁人妻一区二区| 国产精品一区二区精品视频观看| 男女之事视频高清在线观看| 色综合婷婷激情| 亚洲精品粉嫩美女一区| 精品第一国产精品| 少妇被粗大的猛进出69影院| 91麻豆av在线| av中文乱码字幕在线| 99久久国产精品久久久| 日本黄大片高清| 国产精品久久久久久人妻精品电影| 色综合婷婷激情| 国产主播在线观看一区二区| 国模一区二区三区四区视频 | 成人一区二区视频在线观看| 在线永久观看黄色视频| 精品欧美一区二区三区在线| 男女午夜视频在线观看| 99精品久久久久人妻精品| 在线观看66精品国产| 午夜福利在线观看吧| 欧美不卡视频在线免费观看 | 日本三级黄在线观看| 久久久久久人人人人人| 一边摸一边做爽爽视频免费| 91av网站免费观看| 9191精品国产免费久久| 国产午夜精品久久久久久| 国产人伦9x9x在线观看| 啦啦啦免费观看视频1| 老司机午夜福利在线观看视频| 999精品在线视频| or卡值多少钱| 丰满人妻熟妇乱又伦精品不卡| 99re在线观看精品视频| 久久性视频一级片| 看片在线看免费视频| 久久久久亚洲av毛片大全| 免费在线观看黄色视频的| 亚洲av五月六月丁香网| 99精品欧美一区二区三区四区| 天堂影院成人在线观看| 亚洲真实伦在线观看| 国产成人av教育| 在线观看免费视频日本深夜| 色噜噜av男人的天堂激情| 99久久综合精品五月天人人| 亚洲 欧美 日韩 在线 免费| 午夜福利在线观看吧| 窝窝影院91人妻| 色综合亚洲欧美另类图片| 欧美一级毛片孕妇| 日本免费一区二区三区高清不卡| 1024视频免费在线观看| 国产黄片美女视频| 国产一区二区激情短视频| 国产精品电影一区二区三区| 日韩欧美国产一区二区入口| 99riav亚洲国产免费| 亚洲熟女毛片儿| 中文字幕熟女人妻在线| 一本精品99久久精品77| 1024香蕉在线观看| 久久久久国产一级毛片高清牌| 国产一级毛片七仙女欲春2| 中亚洲国语对白在线视频| 成人三级做爰电影| 亚洲精品中文字幕在线视频| 一个人观看的视频www高清免费观看 | 久久精品成人免费网站| 99热这里只有是精品50| 村上凉子中文字幕在线|