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

    沉積條件下白云石與伴生方解石碳氧同位素組成的差別及影響因素

    2016-04-26 07:49:11李小寧黃思靜黃可可鐘怡江胡作維
    石油實驗地質(zhì) 2016年6期
    關(guān)鍵詞:碳氧白云石方解石

    李小寧,黃思靜,黃可可,鐘怡江,胡作維

    (成都理工大學油氣藏地質(zhì)及開發(fā)工程國家重點實驗室,成都 610059)

    沉積條件下白云石與伴生方解石碳氧同位素組成的差別及影響因素

    李小寧,黃思靜,黃可可,鐘怡江,胡作維

    (成都理工大學油氣藏地質(zhì)及開發(fā)工程國家重點實驗室,成都 610059)

    以巖石學為基礎(chǔ),運用元素、X射線衍射和碳氧同位素分析等方法,在手標本中觀尺度上探討同一樣品中白云石與伴生方解石的碳氧同位素組成及影響因素。研究表明,白云石比伴生方解石具有更高的δ13C和δ18O值。白云石的δ18O值高于伴生方解石可能與如下2種因素有關(guān):(1)在白云石—水和方解石—水系統(tǒng)中,前者的氧同位素分餾系數(shù)大于后者,因此在白云石—方解石系統(tǒng)中,氧同位素分餾系數(shù)大于1;(2)白云石沉淀于比方解石鹽度更高、更富18O的水體中。白云石的δ13C值高于伴生方解石可能與如下因素有關(guān):在白云石—CO2和方解石—CO2系統(tǒng)中,前者的碳同位素分餾系數(shù)大于后者,因而對交代過程的白云石—方解石系統(tǒng)而言,碳同位素分餾系數(shù)大于1。準同生白云石與伴生方解石的氧同位素組成差值較大,而熱液/深埋藏白云石與伴生方解石的氧同位素組成差值較小。

    分餾系數(shù);碳氧同位素;伴生方解石;白云石

    白云石自從1791年命名以來,由于近地表常溫常壓條件下實驗室合成白云石的困難性和沉積巖溫壓環(huán)境中白云石的廣泛分布,使得白云石形成機制的研究一直是碳酸鹽沉積學中最困難也是最有魅力的領(lǐng)域之一。這也使得一些和白云石有關(guān)的相應(yīng)研究陷入困境,如白云石—水—CO2(或白云石/方解石)系統(tǒng)碳氧同位素分餾系數(shù)難以通過實驗獲取,造成了人們難以解釋同一流體中形成(沉淀或交代方式)的白云石和伴生方解石之間碳氧同位素的差別。前人利用磷酸與方解石/白云石反應(yīng)速率差異、X射線衍射等方法研究表明,古代白云石與伴生方解石之間存在碳氧同位素差異[1-4],但誰對古海水碳氧同位素組成更具有代表性仍然是亟待解決的問題。如Horacek等[5]在研究南阿爾卑斯下三疊統(tǒng)和二疊系/三疊系界線的碳同位素記錄時認為,白云巖的碳同位素沒有偏離海水,只有氧同位素偏離了海水,顯示白云石和伴生方解石碳氧同位素組成的差別。即使在同一地層中,白云石也比相鄰的方解石具有較正的δ13C和δ18O值,尤其是氧同位素更為明顯。美國猶他州瓦薩奇高原西部邊緣暴露的古新統(tǒng)Flagstaff組下部的白云石比伴生方解石的氧同位素值偏正約7‰[4],從高溫實驗獲得的同位素數(shù)據(jù)外推到25℃時白云石和方解石在平衡狀態(tài)下白云石比方解石富集18O約4‰~7‰[6-8],變質(zhì)巖組合中白云石和方解石在平衡狀態(tài)下白云石與方解石的Δ18O(即白云石與方解石的δ18O值之差)值約為5‰~9‰[2,9-10],白云石與方解石的Δ18O值約為3.2‰[11]。新近的白云石—水—CO2系統(tǒng)與方解石—水—CO2系統(tǒng)碳氧同位素分餾實驗表明,兩者存在碳氧同位素分餾系數(shù)的差別[12-16],說明白云石的氧同位素組成與方解石的差別是分餾系數(shù)的差別所致,它們所代表的仍然是同一流體或類似流體。

    近期我們研究了四川盆地合川鹽井溪下三疊統(tǒng)嘉陵江組的白云巖和石灰?guī)r的碳氧同位素組成,結(jié)果表明白云巖/不完全白云巖化地層比相鄰的灰?guī)r地層具較高的δ13C和δ18O值[17]。但是白云石和方解石分布于不同的地層間隔,無法排除年代效應(yīng)對這2種礦物碳氧同位素差別的影響。基于這種考慮,我們嘗試性地選擇了3件同時存在白云石和方解石的樣品,在手標本中觀尺度上進行方解石和白云石碳氧同位素組成差別的研究,以進一步理解白云石與伴生方解石在碳氧同位素組成上的差別及可能的影響因素。

    1 地質(zhì)背景

    研究區(qū)位于四川盆地東北部宣漢縣渡口鎮(zhèn)(圖1),構(gòu)造上隸屬于大巴山弧形斷褶帶與川東弧形褶皺帶的交會部位[18]。該地區(qū)在下三疊統(tǒng)飛仙關(guān)組沉積時期是典型的盆地—臺地沉積體系[19],開江—梁平海槽東北側(cè)是孤立的碳酸鹽臺地,開江—梁平海槽西南側(cè)是與大陸毗鄰的碳酸鹽臺地,該海槽在飛仙關(guān)組沉積末期完全關(guān)閉。研究區(qū)位于開江—梁平海槽東北側(cè)的孤立碳酸鹽臺地內(nèi)部,臺地周邊為邊緣淺灘環(huán)境,臺地內(nèi)部主要為受限制的潟湖—潮坪環(huán)境。飛仙關(guān)組主要以碳酸鹽巖沉積為主,少量的紫紅色鈣質(zhì)泥巖,飛仙關(guān)組第一段和第四段因局限環(huán)境發(fā)育白云巖或云質(zhì)灰?guī)r等過渡巖石類型。在嘉陵江組沉積時期,研究區(qū)是開闊臺地—局限臺地沉積體系,嘉一、嘉三段以開闊臺地相為主,發(fā)育泥微晶灰?guī)r,嘉二段和嘉四段以局限臺地相為主,發(fā)育白云巖和鹽溶角礫巖,各段之間的巖性地層界線容易識別。樣品采自川東北宣漢渡口剖面下三疊統(tǒng)飛仙關(guān)組第四段(簡稱飛四段)和嘉陵江組第二段(簡稱嘉二段)地層(圖1)。

    圖1 川東北下三疊統(tǒng)飛仙關(guān)組沉積期古地理圖及采樣位置Fig.1 Paleogeographic map of northeastern Sichuan Basin during the deposition of Feixianguan Formation(Early Triassic),also showing sampling locations

    2 巖石學特征

    選取的3組樣品具有不同的巖石學組構(gòu)特征及代表意義。樣品BL30采自飛四段,具泥微晶結(jié)構(gòu)(圖2a),缺乏生物化石和交代結(jié)構(gòu)。研究區(qū)飛四段同時存在代表蒸發(fā)環(huán)境的石膏沉積,表明這些白云石可能是準同生期高鹽度水交代先驅(qū)CaCO3礦物形成的[20]。

    樣品BL103采自嘉二段,具泥微晶藻粒結(jié)構(gòu)(圖2b),顯微鏡下幾乎看不到泥微晶白云石的晶體形態(tài),部分藻粒因富含有機質(zhì)而呈黑色,有機質(zhì)往往富集在白云石分布的區(qū)域,顯示白云石形成機制與藻粒之間的有機聯(lián)系。另一個重要特征是白云石僅僅分布在藻粒中,呈不均勻斑塊狀分布,明顯不同于第一類白云石,此類白云石與Hips等[21]報道的具有微生物結(jié)構(gòu)的原白云石類似;巖石中存在少量孔隙分布在藻粒的內(nèi)部,其形成機制可能與藻尸體腐爛有關(guān),提示巖石經(jīng)歷了大氣水成巖作用。方解石主要分布于藻粒之間,具微亮晶結(jié)構(gòu),薄片下難以判斷其形成的成巖環(huán)境。依據(jù)白云石僅分布在藻粒中而粒間膠結(jié)物沒有白云巖化這一事實,推測膠結(jié)物的初始礦物是具有更好熱力學穩(wěn)定性的低鎂方解石,同時其沉淀作用也可能是在近地表的低Mg/Ca比值的大氣水環(huán)境中發(fā)生的。

    圖2 四川盆地宣漢渡口剖面實驗樣品的顯微照片F(xiàn)ig.2 Microphotographs of experiment samples from Dukou section in Xuanhan county,Sichuan Basin

    樣品BL98A采自嘉二段,具殘余鮞粒結(jié)構(gòu)(圖2c),大小約200~250 μm的鮞粒是主要的粒屑類型。膠結(jié)物以亮晶方解石為主,白云巖化局限在鮞粒的內(nèi)部(圖2d),白云石呈平直晶面自形晶,粒徑分布在50~70 μm。由于地質(zhì)歷史中的大多數(shù)鮞粒都是典型的文石質(zhì)非生物顆粒,顯示白云巖化作用是選擇文石質(zhì)的鮞粒進行的,同時基于鮞粒在巖石中呈漂浮狀這一特征,可以認為膠結(jié)作用發(fā)生的時間很早。結(jié)合樣品非常低的錳(<10×10-6)、鐵(<500×10-6)含量(表1),說明巖石的膠結(jié)作用主要是在準同生期海水潛流環(huán)境中發(fā)生的。

    3 地球化學特征

    3.1 X射線衍射及元素分析

    對3組樣品任意7個區(qū)域進行了X射線衍射分析,結(jié)果表明巖石基本上由白云石和方解石2種礦物組成,21個樣品(3組樣品任意7個區(qū)域)白云石和方解石含量之和變化在96.84%~100%之間,平均值為98.53%,其他礦物主要為少量的石英。基于X射線衍射的分析結(jié)果,我們可以通過CaO、MgO含量來較好地定量計算巖石中方解石和白云石的相對含量,而可以不采用可靠性較差的、僅有半定量意義的X射線衍射分析結(jié)果。同時,為了與化學分析計算獲得的白云石和方解石相對含量對應(yīng),表1中我們也采用由Füchtbauer和Goldschmidt[22]提供的方法計算了白云石和方解石的相對含量,21個樣品的X射線衍射結(jié)果中白云石的104面網(wǎng)間距值(d104)變化在2.877 3~2.891 8之間,顯示這些白云石接近理想組成。因而通過化學分析CaO、MgO含量計算巖石中白云石和方解石的相對含量時,假定白云石具有理想組成(即白云石中CaCO3和MgCO3的摩爾分數(shù)都為50%),不會對計算結(jié)果造成實質(zhì)性影響。

    總體上來講,全部樣品的MgO含量和X射線衍射分析獲得的白云石含量呈正相關(guān)關(guān)系,兩者之間的相關(guān)系數(shù)高達0.96(圖3a)。在3個樣品中,BL103的MgO含量和白云石含量也具有較好的相關(guān)性,相關(guān)系數(shù)達到0.93(圖3b);其余2個樣品相關(guān)性較差,尤其是BL98A樣品,這與樣品中白云石含量較低,X射線衍射半定量分析的誤差相對較大有關(guān)?;瘜W分析中CaO和MgO含量分析的檢測限為0.1%,相對誤差達到了2%,而由X射線衍射分析獲得的礦物含量通常是一種誤差較大的半定量分析?;诖?,后面有關(guān)樣品中白云石和方解石相對含量的討論中,我們采用了由化學分析結(jié)果換算的白云石和方解石相對含量(表1)。

    表1 四川盆地宣漢渡口剖面不同樣品的元素分析、碳氧同位素和X射線衍射分析結(jié)果Table 1 Element,carbon and oxygen isotope values and XRD data of samples from Dukou section in Xuanhan county,Sichuan Basin

    圖3 四川盆地宣漢渡口剖面樣品X射線衍射白云石含量與MgO含量交會圖Fig.3 Cross-plots of dolomite content from XRD and MgO content of samples from Dukou section in Xuanhan county,Sichuan Basin

    3.2 碳氧同位素

    3組樣品的碳氧同位素分析結(jié)果(均為PDB標準,全文相同)見表1,樣品BL30的δ13C值變化在0.82‰~2.26‰,平均值為1.49‰,δ18O值變化在-6.16‰~-4.91‰,平均值為-5.46‰。樣品BL103的δ13C值變化在0.58‰~1.32‰,平均值為0.96‰,δ18O值變化在-9.80‰~-6.58‰,平均值為-7.93‰。樣品BL98A的δ13C值變化在5.19‰~5.47‰,平均值為5.32‰,δ18O值變化在-5.60‰~-5.36‰,平均值為-5.52‰。除了從樣品BL98A獲得的7個數(shù)據(jù)外,其余2個樣品各自7個數(shù)據(jù)的δ13C和δ18O值都有較大的變化范圍,并與元素含量和相應(yīng)的白云石、方解石含量變化具有較好的對應(yīng)變化關(guān)系。

    4 討論

    4.1 白云石含量和碳氧同位素組成的關(guān)系

    樣品BL30和BL103的碳氧同位素組成和白云石含量(或相應(yīng)MgO含量)之間都顯示出良好的線性關(guān)系(圖4),表現(xiàn)為樣品的碳氧同位素組成隨白云石含量變化而變化的特征。按線性函數(shù)關(guān)系,二者之間的相關(guān)系數(shù)在0.7~0.86之間,顯示在手標本中觀尺度上,白云石與伴生方解石具有不同的碳氧同位素組成這一基本事實是存在的,表現(xiàn)為隨著樣品中白云石含量的增加,δ13C和δ18O值增加,白云石比伴生方解石具有更高的δ13C和δ18O值。樣品BL98A的碳氧同位素組成未顯示出隨樣品中白云石含量變化而變化的特征(圖4),我們推測與如下因素有關(guān):

    圖4 四川盆地宣漢渡口剖面不同樣品白云石含量和δ13C、δ18O值的對應(yīng)關(guān)系Fig.4 Correlation of dolomite content and corresponding δ13C,δ18O values of samples from Dukou section in Xuanhan county,Sichuan Basin

    (1)白云石含量較低且具有較小的變化范圍。由MgO換算的白云石含量變化在0.87%~11.96%之間,較小的白云石含量變化導致較小的δ13C和δ18O值變化。分析結(jié)果顯示,樣品MgO含量變化在0.19%~2.60%之間,δ13C值變化在5.19‰~5.47‰之間,δ18O值變化在-5.60‰~-5.36‰之間,碳氧同位素這種較小的變化范圍可能被分析誤差疊加或掩蓋,尤其是該樣品的氧同位素組成,7個分析數(shù)據(jù)的標準偏差為0.081,已接近儀器0.08標準偏差。因而巖石中較低的白云石含量和較小的碳氧同位素變化范圍是這7個樣品白云石含量和同位素組成之間缺乏系統(tǒng)相關(guān)關(guān)系的主要原因。

    (2)白云石形成機制。上述巖石學研究表明,樣品BL98A的白云石具有和BL30、BL103不同的形成機制,后2個樣品中的白云石都是泥微晶結(jié)構(gòu),具有準同生環(huán)境交代白云石或原白云石的特征。BL30樣品的白云巖化流體是高鹽度的海水,BL103樣品與藻類活動有關(guān)的白云石,藻類新陳代謝活動產(chǎn)生的胞外聚合物對白云石的沉淀有著不可忽視的作用。這兩類白云石都是在近地表常溫常壓環(huán)境下形成的,Sun等[24]從磷酸鹽氧同位素反演的南盤江盆地早三疊世古海水的溫度變化在30~40℃,表明當時的古海水溫度較高。而BL98A中的白云石顯然是成巖過程中交代鮞粒形成的,晶體相對較大,溫度可能略高于前兩類白云石形成的溫度,白云巖化作用的流體在時間上和空間上都相對遠離先驅(qū)方解石或文石的沉淀流體,因而白云石含量和碳氧同位素變化之間缺乏協(xié)同變化關(guān)系。不過,在分析誤差和白云石成因2個因素中,分析誤差是主要的,成巖交代白云石和先驅(qū)方解石或文石碳氧同位素之間的關(guān)系仍然需要進一步研究。

    4.2 碳同位素組成和白云石含量的關(guān)系

    對樣品BL30而言,當白云石含量從61.87%變化到77.19%(MgO從13.45%變化到16.78%)時,其δ13C值從0.82‰增加至2.26‰,增加幅度超過1.4‰(圖4a)。如果按圖4a中白云石含量和相應(yīng)δ13C值的函數(shù)關(guān)系,外延到方解石和白云石的端元礦物,則方解石和白云石的δ13C值分別為-3.2‰和3.3‰,白云石和方解石之間碳同位素的差值超過6‰。對樣品 BL103而言,當白云石含量從56.30%變化到89.56%(MgO從12.24%變化到19.47%)時,其δ13C值從0.58‰增加至1.32‰,增加幅度超過0.7‰(圖4a)。如果按圖4a中白云石含量和相應(yīng)δ13C值的函數(shù)關(guān)系,外延到方解石和白云石的端元礦物,則方解石和白云石的δ13C值分別為-0.5‰和1.5‰,白云石和方解石之間碳同位素的差值為2‰。2個樣品不同白云石含量和相應(yīng)δ13C值的變化關(guān)系說明如下問題:

    (1)對于手標本中觀尺度而言,白云石和伴生方解石具有不同的碳同位素組成,白云石的δ13C值高于伴生方解石是可以肯定的事實;

    (2)上述巖石學研究表明2個樣品都不具有交代組構(gòu),可以推測白云石和方解石可能都是從其母流體中沉淀的。假定是從同一母流體中沉淀的,則白云石—CO2的碳同位素分餾系數(shù)應(yīng)大于方解石—CO2系統(tǒng),但二者都是從同一母流體中沉淀的這一假定可能并不成立,因為白云石沉淀時,海水已演化到鹽度更高、Mg/Ca比值更高的階段,無論是按傳統(tǒng)的滲透回流模式還是蒸發(fā)泵模式,二者沉淀時間都存在時間差,由于白云石形成需要較長時間[25],這種時間差已足以造成海水碳同位素組成的差別(后邊將進一步討論這一問題);

    (3)雖然上述巖石學研究表明2個樣品都不具有交代結(jié)構(gòu),但這并不能完全說明這些白云石不是交代成因。傳統(tǒng)的滲透回流模式和蒸發(fā)泵模式中的白云石都是準同生階段交代先驅(qū)方解石或文石形成的,如果是這樣,則說明白云巖化過程中存在碳同位素的分餾,并不是簡單的鎂和鈣之間的置換,可能是一個溶解—沉淀過程,因而也與白云石—CO2和方解石—CO2的碳同位素分餾系數(shù)的差別有關(guān);

    (4)2個樣品的δ13C值和白云石含量(或相應(yīng)MgO含量)的直線函數(shù)具有不同的斜率,分別為0.065和0.02(圖4a),因而樣品BL30中的方解石和白云石具有比BL103中的方解石和白云石更大的δ13C差值,前者為6.5‰,后者為2‰(表2)。我們認為這里疊加了年代效應(yīng),在圖5所示的宣漢渡口剖面下三疊統(tǒng)碳同位素演化曲線中,樣品BL30正處于δ13C快速上升的起點,就海水演化而言,白云石是形成于(沉淀或交代)伴生方解石或文石之后,海水此時具有更高的鹽度和Mg/Ca比值,白云石形成時的海水顯然已具有更高的δ13C值。由于早三疊世Induan期(可大致與上揚子地區(qū)的飛仙關(guān)組、大冶組或夜郎組對應(yīng)[26])持續(xù)時間只有約1 Ma時間段[27],因而BL30樣品所在的Induan期末期的海水δ13C值快速上升持續(xù)的時間可能達到萬年的尺度范圍。

    表2 四川盆地宣漢渡口剖面樣品的白云石含量與碳氧同位素之間的函數(shù)關(guān)系和相應(yīng)的同位素值Table 2 Function expression of dolomite content and carbon and oxygen isotopes and corresponding isotope values of BL30 and BL103 samples from Dukou section in Xuanhan county,Sichuan Basin

    圖5 四川盆地宣漢渡口剖面下三疊統(tǒng)碳同位素演化曲線Fig.5 Carbon isotope evolution curves for Early Triassic,Dukou section in Xuanhan county,Sichuan Basin

    4.3 氧同位素組成和白云石含量的關(guān)系及與碳同位素的對比

    對樣品BL30而言,白云石含量從61.87%變化到77.19%(MgO從13.45%變化到16.78%)時,其δ18O值從-6.16‰增加至-4.91‰,增加幅度超過1.2‰(圖4b)。如果按圖4b中白云石含量和相應(yīng)δ18O值的函數(shù)關(guān)系,外延到方解石和白云石的端元礦物,則其方解石和白云石的δ18O值分別為-10.54‰和-3.54‰,白云石和方解石之間氧同位素的差值為7‰。對樣品BL103而言,白云石含量從56.30%變化到89.56%(MgO從12.24%變化到19.47%)時,其δ18O值從-9.80‰增加至-6.58‰,增加幅度超過3‰(圖4b)。如果按圖4b中白云石含量和相應(yīng)δ18O值的函數(shù)關(guān)系,外延到方解石和白云石的端元礦物,則其方解石和白云石的δ18O值分別為-15.39‰和-5.39‰,白云石和方解石之間氧同位素的差值為 10‰。按照實驗得出的 BL30和BL103樣品端元方解石的氧同位素值,根據(jù)已有的白云石—方解石系統(tǒng)中氧同位素分餾方程[12],假設(shè)白云石的形成溫度是30℃,計算獲得的白云石的氧同位素分別為-7‰和-11.87‰,與兩者各自的方解石氧同位素差值分別為3.54‰和3.52‰,可見白云石—方解石氧同位素分餾方程計算的白云石與方解石氧同位素之間的差值要小于我們實驗方法獲得的。這與計算過程中參數(shù)選取的合理性、白云石—方解石氧同位素分餾方程的由來以及地質(zhì)事實中白云石形成的復雜性相關(guān)(如與白云石伴生的方解石可能在一定程度上受到大氣水的影響而其δ18O值更低)。

    從另一方面也說明,我們需要加強地質(zhì)事實中白云石—方解石—水系統(tǒng)的氧同位素分餾方程的研究。2個樣品不同白云石含量和相應(yīng)δ18O值的變化關(guān)系說明如下問題:

    (1)對于手標本中觀尺度而言,白云石和伴生方解石具有不同的氧同位素組成,白云石的δ18O值高于伴生方解石是可以肯定的事實;

    (2)如果白云石和方解石可能都是從同一母流體中沉淀的,則說明白云石—水系統(tǒng)的氧同位素分餾系數(shù)大于方解石—水系統(tǒng),如果白云石和方解石不是從同一母流體中沉淀的,除了說明白云石—水系統(tǒng)的氧同位素分餾系數(shù)大于方解石—水系統(tǒng)以外,還可能說明白云石形成(沉淀或交代先驅(qū)文石或方解石)的流體的鹽度高于伴生方解石的沉淀流體,當然相應(yīng)的Mg/Ca比值也會更高。

    (3)2個樣品的δ18O值和白云石含量(或MgO相應(yīng)含量)的直線函數(shù)具有不同的斜率,分別為0.07和0.1,因而樣品BL103中的方解石和白云石具有比BL30中的方解石和白云石更大的δ18O差值,前者為10‰,后者為7‰,同時BL103外延至方解石后,δ18O值低至-15.39‰(表2)。結(jié)合巖石學研究說明該樣品藻粒之間的方解石可能是大氣水成因,白云石和方解石分別沉淀于具不同鹽度的水體中是造成這種現(xiàn)象的主要原因。

    (4)對比氧同位素和碳同位素組成與白云石含量之間的變化關(guān)系,可以看出,白云石與伴生方解石氧同位素組成的差值都大于碳同位素組成,尤其是樣品BL103,白云石和方解石δ13C和δ18O的差值分別為2‰和10‰,這主要與白云石或方解石—水的氧同位素分餾系數(shù)大于白云石或方解石—CO2的碳同位素分餾系數(shù)有關(guān)。

    4.4 溫度對氧同位素組成差別的影響

    雖然我們研究的樣品中白云石和伴生方解石都形成于沉積溫度條件下,但我們可以利用前人實驗與理論研究成果,來進行溫度對白云巖化過程中方解石和白云石氧同位素組成差別影響的討論。較新的研究成果是由Horita完成的[12],在綜合方解石—水[28-29]和白云石—水氧同位素分餾方程的基礎(chǔ)上,(原)白云石—方解石的18O/16O分餾系數(shù)可以用如下方程表達:

    式中:α為分餾系數(shù),T為開氏溫度。

    雖然該氧同位素分餾方程是在80~350℃溫度范圍下獲得的,但氧同位素分餾系數(shù)與溫度1/T2呈線性關(guān)系,可以外延至其他溫度范圍。顯然該研究有助于我們理解白云巖化過程中方解石和白云石氧同位素組成的差別,以及溫度對這種差別的影響。我們?nèi)?個溫度點來討論這種差別,分別是30,80,150,350℃,假定先驅(qū)方解石的δ18O值為0‰,則在這4種溫度條件下形成的白云石的δ18O值分別為3.6‰,2.6‰,1.7‰,0.7‰,這4個數(shù)值顯然也就是在這4種溫度條件下白云石和先驅(qū)方解石的δ18O的差值,顯示溫度越高,白云巖化過程中方解石和白云石δ18O值的差別越小。因此,準同生白云石與伴生方解石的氧同位素組成差值較大,而熱液/深埋藏白云石與伴生方解石的氧同位素組成差值較小。

    5 結(jié)論

    (1)在手標本中觀尺度上,同一樣品中白云石與伴生方解石具有不同的碳氧同位素組成,因而隨著樣品中白云石含量的增加,δ13C和δ18O值增加,顯示白云石比伴生方解石具有更高的δ13C和δ18O值。

    (2)白云石的δ18O值高于伴生方解石可能與如下2種因素有關(guān):①白云石—水系統(tǒng)中和方解石—水系統(tǒng)中,前者的氧同位素分餾系數(shù)應(yīng)大于后者,對交代過程而言,白云石—方解石系統(tǒng)中,氧同位素分餾系數(shù)大于1;②白云石沉淀于比方解石鹽度更高、更富18O的水體中,因而白云石也更富18O。

    (3)白云石的δ13C值高于伴生方解石可能與如下因素有關(guān):白云石—CO2系統(tǒng)中和方解石—CO2系統(tǒng)中,前者的碳同位素分餾系數(shù)大于后者,對交代過程而言,白云石—方解石系統(tǒng)中,碳同位素分餾系數(shù)大于1?;诖?,我們不能簡單地用白云石或方解石的碳同位素去代替古海水,因而需要進一步加強沉積溫壓條件下白云石—方解石—CO2系統(tǒng)碳同位素分餾系數(shù)的定量研究。

    (4)由于白云巖化作用需要時間,即使是萬年尺度的年齡變化,對于海水碳同位素組成大幅度、高頻率變化的時間間隔來說,已足以造成手標本中觀尺度白云石和伴生方解石碳同位素組成的差別,這在早三疊世時間間隔顯得尤其明顯,因此,手標本中觀尺度白云石和方解石碳同位素組成的差別也可能是時間造成的。

    (5)相對低溫條件的白云巖化過程中,白云石與先驅(qū)方解石氧同位素的差值較大,隨著白云巖化過程溫度的升高,白云石與先驅(qū)方解石之間氧同位素組成的差值減小,因而準同生白云石與伴生方解石的氧同位素組成差值較大,而熱液/深埋藏白云石與伴生方解石的氧同位素組成差值較小。

    [1]Degens E T,Epstein S.Oxygen and carbon isotope ratios in coexisting calcites and dolomites from recent and ancient sediments[J].Geochimica et Cosmochimica Acta,1964,28(1):23-44.

    [2]Sheppard S M F,Schwarcz H P.Fractionation of carbon and oxygen isotopes and magnesium between coexisting metamorphic calcite and dolomite[J].Contributions to Mineralogy and Petrology,1970,26(3):161-198.

    [3]Clayton R N,Jones B F.Isotope studies of dolomite formation under sedimentary conditions[J].Geochimica et Cosmochimica Acta,1968,32(4):415-432.

    [4]Weber J N.Oxygen isotope fractionation between coexisting calcite and dolomite[J].Science,1964,145(3638):1303-1305.

    [5]Horacek M,Brandner R,Abart R.Carbon isotope record of the P/T boundary and the Lower Triassic in the southern Alps:Evidence for rapid changes in storage of organic carbon[J].Palaeogeography,Palaeoclimatology,Palaeoecology,2007,252(1/2):347-354.

    [6]Epstein S,Graf D L,Degens E T.Oxygen isotope studies on the origin of dolomites[C]//Isotopic and cosmic chemistry.Amsterdam:North-Holland Publishing Company,1964:169-180.

    [7]Northrop D A,Clayton R N.Oxygen-isotope fractionation in systems containing dolomite[J].The Journal of Geology,1966,74(2):174-196.

    [8]O’Neil J R,Epstein S.Oxygen isotope fractionation in the system dolomite-calcite-carbon dioxide[J].Science,1966,152 (3719):198-201.

    [9]Clayton R A,Epstein S.The relationship between O18/O16ratios in coexisting quartz,carbonate,and iron oxides from various geological deposits[J].The Journal of Geology,1958,66(4):352-371.

    [10]Engel A E J,Clayton R N,Epstein S.Variations in the isotopic composition of oxygen in the Leadville limestone(Mississippian,Colorado)and in its hydrothermal and metamorphic phases[J].The Journal of Geology,1958,66(4):374-393.

    [11]McKenzie J A.Holocene dolomitization of calcium carbonate sediments from the coastal Sabkhas of Abu Dhabi,U.A.E.:A stable isotope study[J].The Journal of Geology,1981,89(2):185-198.

    [12]Horita J.Oxygen and carbon isotope fractionation in the system dolomite-water-CO2to elevated temperatures[J].Geochimica et Cosmochimica Acta,2014,129:111-124.

    [13]Deines P.Carbon isotope effects in carbonate systems[J].Geochimica et Cosmochimica Acta,2004,68(12):2659-2679.

    [14]Jiménez-López C,Romanek C S,Huertas F J,et al.Oxygen isotope fractionation in synthetic magnesian calcite[J].Geochimica et Cosmochimica Acta,2004,68(16):3367-3377.

    [15]Mavromatis V,Schmidt M,Botz R,et al.Experimental quantification of the effect of Mg on calcite–aqueous fluid oxygen isotope fractionation[J].Chemical Geology,2012,310/311:97-105.

    [16]Chacko T,Deines P.Theoretical calculation of oxygen isotope fractionation factors in carbonate systems[J].Geochimica et Cosmochimica Acta,2008,72(15):3642-3660.

    [17]李小寧,黃思靜,胡博,等.合川鹽井溪嘉陵江組第二段鹽溶角礫巖—次生灰?guī)r的碳氧同位素組成[J].成都理工大學學報(自然科學版),2015,42(6):734-745.

    Li Xiaoning,Huang Sijing,Hu Bo,et al.Carbon and oxygen isotope composition of evaporite-solution breccia-secondary limestone of Member 2 of Triassic Jialingjiang Formation in Hechuan,eastern Sichuan Basin,China[J].Journal of Chengdu University of Technology(Science& Technology Edition),2015,42(6):734-745.

    [18]曾云賢,劉微,李西寧.川東北地區(qū)飛仙關(guān)組構(gòu)造演化與油氣成藏研究[J].天然氣工業(yè),2015,25(S1):21-24.

    Zeng Yunxian,Liu Wei,Li Xining.Research on structural evolution and oil/gas reservoiring of Feixianguan Formation in northwest Sichuan Basin[J].Natural Gas Industry,2015,25(S1): 21-24.

    [19]王興志,張帆,馬青,等.四川盆地東部晚二疊世—早三疊世飛仙關(guān)組礁、灘特征與海平面變化[J].沉積學報,2002,20(2):249-254.

    Wang Xingzhi,Zhang Fan,Ma Qing,et al.The characteristics of reef and bank and the fluctuation of sea-level in Feixianguan period of Late Permian-Early Triassic,east Sichuan Basin[J].Acta Sedimentologia Sinica,2002,20(2):249-254.

    [20]陳琪,胡文瑄,李慶,等.川東北盤龍洞長興組—飛仙關(guān)組白云巖化特征及成因[J].石油與天然氣地質(zhì),2014,33(1): 84-93.

    Chen Qi,Hu Wenxuan,Li Qing,et al.Characteristics and genesis of dolomitization in Changxing and Feixianguan formations in Panlongdong,northeastern Sichuan Basin[J].Oil&Gas Geology,2014,33(1):84-93.

    [21]Hips K,Haas J,Poros Z,et al.Dolomitization of Triassic microbial mat deposits(Hungary):Origin of microcrystalline dolomite[J].Sedimentary Geology,2015,318:113-129.

    [22]Füchtbauer H,Goldschmidt H.Beziehungen zwischen calciumgehalt und bildungsbedingungen der dolomite[J].Geologische Rundschau,1966,55(1):29-40.

    [23]黃可可,黃思靜,胡作維,等.四川盆地宣漢渡口和重慶北碚下三疊統(tǒng)海相碳酸鹽碳同位素組成與演化[J].古地理學報,2016,18(1):101-114.

    Huang Keke,Huang Sijing,Hu Zuowei,et al.Carbon isotopic composition and evolution of the Lower Triassic marine carbonates from Dukou of Xuanhan and Beibei of Chongqing,Sichuan Basin[J].Journal of Palaeogeography,2016,18(1):101-114.

    [24]Sun Yadong,Joachimski M M,Wignall P B,et al.Lethally hot temperatures during the Early Triassic greenhouse[J].Science,2012,338(6105):366-370.

    [25]Warren J.Dolomite:Occurrence,evolution and economically importantassociations[J].Earth-Science Reviews,2000,52(1/3):1-81.

    [26]祝幼華,羅輝,蔡華偉,等.四川廣安謝家槽海相早、中三疊世地層劃分[J].地層學雜志,2012,36(4):784-791.

    Zhu Youhua,Luo Hui,Cai Huawei,et al.Stratigraphic division of the Early and Middle Triassic at the Xiejiacao section in Guang’an,Sichuan[J].Journal of Stratigraphy,2012,36(4): 784-791.

    [27]International Commission on Stratigraphy.International chronostratigraphic chart[EB/OL].[2015-11-04].www.stratigraphy.org.

    [28]O’Neil J R,Clayton R N,Mayeda T K.Oxygen isotope fractionation in divalent metal carbonates[J].The Journal of Chemical Physics,1969,51(12):5547-5558.

    [29]Friedman I,O’Neil J R.Compilation of stable isotope fractionation factors of geochemical interest[R].Washington:United States Government Printing Office,1977:KK1-KK12.

    (編輯 徐文明)

    Differences and controls of carbon and oxygen isotope composition in dolomite and coexisting calcite under deposition conditions

    Li Xiaoning,Huang Sijing,Huang Keke,Zhong Yijiang,Hu Zuowei
    (State Key Laboratory of Oil/Gas Reservoir Geology and Exploitation,Chengdu University of Technology,Chengdu,Sichuan 610059,China)

    We discussed the composition and factors influencing carbon and oxygen isotopes in dolomites and coexisting calcites in hand specimen mesoscale based on petrography and combined with elemental,X-ray diffraction and carbon and oxygen isotope analyses.Results indicated that dolomites have higher δ13C and δ18O values than coexisting calcites.Two factors may account for the higher δ18O values of dolomites.(1)The oxygen isotope fractionation factors of a dolomite-water system are higher than those of a calcite-water system.Therefore,oxygen isotope fractionation factors are greater than 1 in a dolomite-calcite system.(2)The fluids of dolomite precipitation have higher salinity and higher18O abundance compared to coexisting calcite.In addition,the δ13C values of dolomites show the same trends with δ18O values.That is to say,the carbon isotope fractionation factors in a dolomite-CO2system are higher than those in a calcite-CO2system.Therefore,carbon isotope fractionation factors are greater than 1 in a dolomitecalcite system during replacement process.Oxygen isotope values varied obviously between syndepositional dolomites and coexisting calcites,while they were similar between hydrothermal/deeply burial dolomites and coexisting calcites.

    fractionation factors;carbon and oxygen isotope;coexisting calcite;dolomite

    TE12<2.2 class="emphasis_bold">2.2 文獻標識碼:A2.2

    A

    1001-6112(2016)06-0828-08

    10.11781/sysydz201606828

    2016-02-06;

    2016-09-28。

    李小寧(1986—),女,博士生,從事沉積地球化學研究。E-mail:l-x-n@139.com。

    黃思靜(1949—),男,教授,博士生導師,從事沉積學的教學與科研工作。E-mail:hsj@cdut.edu.cn。

    國家自然科學基金項目(41272130)“川渝地區(qū)早三疊世海相碳酸鹽巖的碳同位素研究”資助。

    猜你喜歡
    碳氧白云石方解石
    水泥回轉(zhuǎn)窯用白云石質(zhì)耐火材料的研究進展
    工業(yè)爐(2023年6期)2023-03-29 05:29:18
    不只是詞典筆 碳氧全科學習筆大屏版
    F-在方解石表面的吸附及其對方解石表面性質(zhì)的影響
    硅酸鹽通報(2022年8期)2022-09-08 04:25:42
    遵化市魏家井白云石礦冶金用白云巖礦地質(zhì)特征
    碳氧血紅蛋白在新生兒ABO溶血性黃疸中的臨床意義
    氯化鈣和碳酸鈉對方解石浮選的影響及其機理研究
    貴州重晶石與方解石常溫浮選分離試驗研究
    白云石有序度與流變特征的研究進展
    川西甘孜-理塘結(jié)合帶碳酸鹽巖碳氧同位素特征
    碳氧血紅蛋白與急性一氧化碳中毒分級診治的關(guān)系
    久久国产乱子免费精品| 国产在视频线精品| 久久人人爽人人片av| 精品人妻一区二区三区麻豆| 色5月婷婷丁香| 亚洲电影在线观看av| 亚洲欧美精品专区久久| 大片电影免费在线观看免费| 在线 av 中文字幕| av在线app专区| 国产精品一区二区三区四区免费观看| 精华霜和精华液先用哪个| 欧美另类一区| 国产精品嫩草影院av在线观看| 久热久热在线精品观看| 中国美白少妇内射xxxbb| 22中文网久久字幕| 中文精品一卡2卡3卡4更新| 国模一区二区三区四区视频| 尤物成人国产欧美一区二区三区| av福利片在线观看| 深夜a级毛片| 嫩草影院精品99| 亚洲成人久久爱视频| 亚洲三级黄色毛片| 欧美三级亚洲精品| 美女国产视频在线观看| 毛片一级片免费看久久久久| 新久久久久国产一级毛片| 成人免费观看视频高清| 亚洲精品日韩av片在线观看| 91精品伊人久久大香线蕉| 极品教师在线视频| 久久久久久久久久成人| 亚洲精品久久久久久婷婷小说| 精品一区二区免费观看| 国产伦在线观看视频一区| 国产永久视频网站| 国产淫片久久久久久久久| 91aial.com中文字幕在线观看| 国产黄片美女视频| 久久久欧美国产精品| 日本三级黄在线观看| 青春草国产在线视频| 色婷婷久久久亚洲欧美| 插阴视频在线观看视频| 久久久久网色| 少妇的逼水好多| 国产成人精品久久久久久| 亚洲av男天堂| 在线 av 中文字幕| 一级爰片在线观看| 熟女电影av网| 又粗又硬又长又爽又黄的视频| 毛片一级片免费看久久久久| videos熟女内射| 国产欧美亚洲国产| 日韩中字成人| 亚洲aⅴ乱码一区二区在线播放| 嘟嘟电影网在线观看| 97超碰精品成人国产| 久久精品人妻少妇| 国国产精品蜜臀av免费| 亚洲国产色片| 久久亚洲国产成人精品v| 亚洲精品自拍成人| 久久久久久久午夜电影| 国产成人aa在线观看| 亚洲丝袜综合中文字幕| 边亲边吃奶的免费视频| 精品久久久久久久久av| 人妻少妇偷人精品九色| 秋霞伦理黄片| 亚洲真实伦在线观看| 青春草亚洲视频在线观看| 乱系列少妇在线播放| 精品视频人人做人人爽| 亚洲av欧美aⅴ国产| 2018国产大陆天天弄谢| 国产老妇伦熟女老妇高清| 亚洲人成网站高清观看| 身体一侧抽搐| 久久影院123| 少妇熟女欧美另类| 亚洲精华国产精华液的使用体验| 国产一区二区三区av在线| 国产精品女同一区二区软件| 内地一区二区视频在线| 一个人看的www免费观看视频| 亚洲欧美成人综合另类久久久| 国产成人a∨麻豆精品| 免费黄网站久久成人精品| 伦理电影大哥的女人| 最新中文字幕久久久久| 国产成人a∨麻豆精品| 久久精品国产鲁丝片午夜精品| 亚洲成人av在线免费| 最近最新中文字幕免费大全7| 午夜老司机福利剧场| 亚洲天堂国产精品一区在线| 午夜福利高清视频| 国产精品国产三级国产av玫瑰| 亚洲综合精品二区| 国产免费又黄又爽又色| 两个人的视频大全免费| 亚洲欧洲日产国产| 一级黄片播放器| 欧美老熟妇乱子伦牲交| 亚洲欧美精品专区久久| 99九九线精品视频在线观看视频| 国产爽快片一区二区三区| 国产精品一区二区在线观看99| 99热网站在线观看| 国产亚洲精品久久久com| 免费人成在线观看视频色| 精品人妻偷拍中文字幕| 91久久精品国产一区二区三区| 国产乱来视频区| 色5月婷婷丁香| 日本午夜av视频| 国产精品久久久久久av不卡| 国产色婷婷99| 日本与韩国留学比较| 久久99热这里只频精品6学生| 精品99又大又爽又粗少妇毛片| .国产精品久久| 日韩大片免费观看网站| 亚洲精品日本国产第一区| 波多野结衣巨乳人妻| 欧美亚洲 丝袜 人妻 在线| 国产综合精华液| 国产爽快片一区二区三区| 九九在线视频观看精品| 777米奇影视久久| 狂野欧美激情性bbbbbb| 色5月婷婷丁香| 欧美bdsm另类| 国产淫片久久久久久久久| 成人综合一区亚洲| 精品国产一区二区三区久久久樱花 | 久久久久国产精品人妻一区二区| 人妻制服诱惑在线中文字幕| 大码成人一级视频| 狂野欧美白嫩少妇大欣赏| 亚洲性久久影院| 99热国产这里只有精品6| 97热精品久久久久久| 一级爰片在线观看| 亚洲精品日韩在线中文字幕| 91在线精品国自产拍蜜月| 亚洲精品一区蜜桃| 三级国产精品欧美在线观看| 久久久久久久精品精品| 99热这里只有精品一区| 成年女人看的毛片在线观看| 欧美一区二区亚洲| 国产欧美日韩一区二区三区在线 | av女优亚洲男人天堂| 日韩电影二区| 国内揄拍国产精品人妻在线| 水蜜桃什么品种好| 亚洲国产日韩一区二区| 亚洲av.av天堂| 亚洲国产精品成人久久小说| 免费大片黄手机在线观看| 高清日韩中文字幕在线| 黄色视频在线播放观看不卡| 国产精品.久久久| 色吧在线观看| 成年免费大片在线观看| 黄色一级大片看看| 亚洲精品色激情综合| 亚洲欧美日韩另类电影网站 | 亚洲国产高清在线一区二区三| 久久久久久国产a免费观看| 成年免费大片在线观看| 色网站视频免费| 男女边吃奶边做爰视频| 欧美成人精品欧美一级黄| 国产av码专区亚洲av| 亚洲一级一片aⅴ在线观看| 欧美zozozo另类| 国模一区二区三区四区视频| 91aial.com中文字幕在线观看| 国产 一区精品| av又黄又爽大尺度在线免费看| 欧美一区二区亚洲| 国产在视频线精品| 干丝袜人妻中文字幕| 中文字幕久久专区| h日本视频在线播放| 亚洲欧美日韩无卡精品| 97超视频在线观看视频| 精品少妇久久久久久888优播| 亚洲四区av| 日本三级黄在线观看| 少妇的逼水好多| .国产精品久久| 欧美日韩在线观看h| 激情五月婷婷亚洲| 日韩亚洲欧美综合| 日本一本二区三区精品| 在线a可以看的网站| 国产精品久久久久久精品电影| 精品人妻熟女av久视频| 日韩国内少妇激情av| 亚洲综合色惰| 久久久久久久大尺度免费视频| 欧美日韩一区二区视频在线观看视频在线 | 欧美性猛交╳xxx乱大交人| 欧美成人午夜免费资源| 日韩一本色道免费dvd| 亚洲国产精品成人综合色| 精品久久久噜噜| 街头女战士在线观看网站| 日韩欧美精品免费久久| 欧美高清成人免费视频www| 亚洲国产av新网站| 国产精品99久久99久久久不卡 | 麻豆成人午夜福利视频| 免费看光身美女| 大码成人一级视频| 如何舔出高潮| 日韩av不卡免费在线播放| 少妇猛男粗大的猛烈进出视频 | 校园人妻丝袜中文字幕| 国产色婷婷99| 久久久久久久久大av| 麻豆乱淫一区二区| 九九在线视频观看精品| 少妇裸体淫交视频免费看高清| 久久久a久久爽久久v久久| 日韩欧美精品v在线| 高清欧美精品videossex| 人人妻人人爽人人添夜夜欢视频 | 欧美精品人与动牲交sv欧美| 国产在视频线精品| 国产探花极品一区二区| 午夜激情久久久久久久| 成人亚洲精品av一区二区| 国产女主播在线喷水免费视频网站| 亚洲欧美日韩卡通动漫| 久久这里有精品视频免费| 九九久久精品国产亚洲av麻豆| 最近2019中文字幕mv第一页| 在线看a的网站| 丝袜美腿在线中文| 成人午夜精彩视频在线观看| 丰满人妻一区二区三区视频av| 婷婷色av中文字幕| 成人毛片60女人毛片免费| 久久99热6这里只有精品| 97超碰精品成人国产| 各种免费的搞黄视频| 婷婷色av中文字幕| 边亲边吃奶的免费视频| 亚洲高清免费不卡视频| 久久亚洲国产成人精品v| 亚洲av福利一区| 狠狠精品人妻久久久久久综合| 亚洲精品视频女| 国模一区二区三区四区视频| 久久久久久久久久成人| 麻豆国产97在线/欧美| 亚洲国产精品国产精品| 汤姆久久久久久久影院中文字幕| 久久6这里有精品| 国产一区二区在线观看日韩| 中文资源天堂在线| 91午夜精品亚洲一区二区三区| 成人特级av手机在线观看| 毛片女人毛片| 欧美高清性xxxxhd video| 国产精品不卡视频一区二区| 永久网站在线| 99热网站在线观看| 毛片一级片免费看久久久久| 日本三级黄在线观看| 五月开心婷婷网| 精华霜和精华液先用哪个| 国产高清有码在线观看视频| 精品亚洲乱码少妇综合久久| 成人亚洲欧美一区二区av| 精品久久久久久电影网| 国产高清有码在线观看视频| 午夜爱爱视频在线播放| 欧美3d第一页| 91精品一卡2卡3卡4卡| 永久免费av网站大全| 国产成人精品婷婷| 午夜爱爱视频在线播放| 午夜视频国产福利| av在线观看视频网站免费| 久久影院123| 久久精品久久久久久久性| 黄色一级大片看看| 七月丁香在线播放| 国产精品久久久久久精品电影小说 | 国产精品久久久久久久电影| 国产爱豆传媒在线观看| 国产国拍精品亚洲av在线观看| 97热精品久久久久久| 久久久久久久久久久丰满| 国产av码专区亚洲av| 国产成年人精品一区二区| av女优亚洲男人天堂| 成人午夜精彩视频在线观看| 免费看a级黄色片| 日本爱情动作片www.在线观看| 欧美成人午夜免费资源| 国产成人freesex在线| 久久久亚洲精品成人影院| 黑人高潮一二区| 久久久久九九精品影院| 国产一区有黄有色的免费视频| 看十八女毛片水多多多| 亚洲欧洲国产日韩| 欧美日韩国产mv在线观看视频 | 美女视频免费永久观看网站| 亚洲色图av天堂| 日韩三级伦理在线观看| 国产精品嫩草影院av在线观看| 夜夜爽夜夜爽视频| 久久热精品热| 久久99热这里只频精品6学生| 熟妇人妻不卡中文字幕| 大话2 男鬼变身卡| 夫妻午夜视频| 麻豆精品久久久久久蜜桃| 夫妻性生交免费视频一级片| 少妇的逼好多水| 啦啦啦啦在线视频资源| 国产有黄有色有爽视频| 久久久久久国产a免费观看| 97精品久久久久久久久久精品| 国产在线一区二区三区精| 成人亚洲欧美一区二区av| 26uuu在线亚洲综合色| 一区二区三区乱码不卡18| av黄色大香蕉| 亚洲第一区二区三区不卡| 精品久久久久久久久av| 午夜激情福利司机影院| 亚洲国产日韩一区二区| 丰满人妻一区二区三区视频av| 国产精品久久久久久精品电影小说 | 久久人人爽人人片av| 欧美激情久久久久久爽电影| 日韩中字成人| 王馨瑶露胸无遮挡在线观看| 欧美极品一区二区三区四区| www.av在线官网国产| 天天躁夜夜躁狠狠久久av| av在线天堂中文字幕| 亚洲不卡免费看| 国产精品一区二区性色av| 久久久国产一区二区| 极品教师在线视频| 亚洲成人av在线免费| 国产精品国产三级国产av玫瑰| 日本黄色片子视频| 天堂中文最新版在线下载 | 啦啦啦中文免费视频观看日本| 51国产日韩欧美| 亚洲国产精品999| 亚洲欧美日韩无卡精品| 18禁在线无遮挡免费观看视频| 一级毛片aaaaaa免费看小| 国产午夜精品久久久久久一区二区三区| 噜噜噜噜噜久久久久久91| 晚上一个人看的免费电影| 亚洲欧美一区二区三区黑人 | 九九爱精品视频在线观看| 人妻夜夜爽99麻豆av| 国产高清有码在线观看视频| 国产一区有黄有色的免费视频| 免费高清在线观看视频在线观看| 亚洲高清免费不卡视频| a级一级毛片免费在线观看| 婷婷色综合大香蕉| 夜夜看夜夜爽夜夜摸| 噜噜噜噜噜久久久久久91| 高清在线视频一区二区三区| 男人舔奶头视频| 搡老乐熟女国产| 涩涩av久久男人的天堂| 亚洲精品aⅴ在线观看| 亚洲最大成人av| 免费看日本二区| 亚洲国产日韩一区二区| 18禁裸乳无遮挡动漫免费视频 | 亚洲三级黄色毛片| 日日啪夜夜撸| 春色校园在线视频观看| 高清视频免费观看一区二区| 久久99热这里只有精品18| 人妻少妇偷人精品九色| 麻豆成人av视频| 色吧在线观看| 又爽又黄无遮挡网站| 日本黄大片高清| 舔av片在线| 国产久久久一区二区三区| 欧美精品人与动牲交sv欧美| videossex国产| 国产毛片在线视频| 国产精品秋霞免费鲁丝片| 少妇人妻精品综合一区二区| 黄色配什么色好看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日本欧美国产在线视频| 亚洲久久久久久中文字幕| 成人毛片a级毛片在线播放| 搡老乐熟女国产| 亚洲欧美一区二区三区国产| 亚洲国产日韩一区二区| 久久久成人免费电影| 免费在线观看成人毛片| 国产免费福利视频在线观看| 亚洲高清免费不卡视频| 国产伦精品一区二区三区四那| 国产精品福利在线免费观看| 亚洲av一区综合| 丝袜喷水一区| 成人亚洲精品一区在线观看 | 国产成人一区二区在线| www.色视频.com| 简卡轻食公司| 久久鲁丝午夜福利片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久欧美国产精品| 99久久中文字幕三级久久日本| 国产v大片淫在线免费观看| 精品熟女少妇av免费看| 国产成人精品婷婷| 久久精品国产亚洲网站| 欧美xxⅹ黑人| 成人漫画全彩无遮挡| 婷婷色av中文字幕| 全区人妻精品视频| 欧美激情在线99| 小蜜桃在线观看免费完整版高清| 国产午夜精品久久久久久一区二区三区| 亚洲欧美清纯卡通| 99re6热这里在线精品视频| 中文精品一卡2卡3卡4更新| 韩国av在线不卡| 两个人的视频大全免费| 最近手机中文字幕大全| 日韩成人伦理影院| 亚洲欧美成人综合另类久久久| 春色校园在线视频观看| 伊人久久精品亚洲午夜| 久久久国产一区二区| 在线亚洲精品国产二区图片欧美 | 成年女人看的毛片在线观看| 成人亚洲精品av一区二区| 热re99久久精品国产66热6| 国产熟女欧美一区二区| 久久人人爽av亚洲精品天堂 | 极品教师在线视频| 欧美精品人与动牲交sv欧美| 日韩欧美精品免费久久| 人妻制服诱惑在线中文字幕| 舔av片在线| a级一级毛片免费在线观看| 大片电影免费在线观看免费| tube8黄色片| 久久久久久久久大av| 午夜精品国产一区二区电影 | 深爱激情五月婷婷| 秋霞伦理黄片| 亚洲av福利一区| 国产大屁股一区二区在线视频| 国产极品天堂在线| 色婷婷久久久亚洲欧美| 国产毛片a区久久久久| 国产精品爽爽va在线观看网站| 99久久人妻综合| 国产精品无大码| 人妻系列 视频| 午夜福利在线观看免费完整高清在| 国产av国产精品国产| 亚洲无线观看免费| 男女国产视频网站| 精品久久久久久久人妻蜜臀av| 精品久久久久久久久av| 在线亚洲精品国产二区图片欧美 | 亚洲最大成人中文| 嫩草影院精品99| 国产免费又黄又爽又色| 精品一区二区免费观看| 极品教师在线视频| 欧美亚洲 丝袜 人妻 在线| 久久精品国产亚洲网站| 亚洲av中文av极速乱| 黄色视频在线播放观看不卡| 亚洲精品aⅴ在线观看| 在线免费十八禁| 麻豆久久精品国产亚洲av| 久久久久网色| 成人毛片60女人毛片免费| 毛片女人毛片| 国产老妇女一区| 18禁在线无遮挡免费观看视频| 国产成人精品福利久久| 国产高清国产精品国产三级 | 狠狠精品人妻久久久久久综合| 边亲边吃奶的免费视频| 国语对白做爰xxxⅹ性视频网站| 日韩,欧美,国产一区二区三区| 欧美性猛交╳xxx乱大交人| 99久久精品一区二区三区| 又粗又硬又长又爽又黄的视频| 免费黄网站久久成人精品| 少妇高潮的动态图| 欧美另类一区| 免费大片黄手机在线观看| 在现免费观看毛片| 亚洲精品成人久久久久久| 国产成人a区在线观看| 超碰av人人做人人爽久久| 国产熟女欧美一区二区| 久久精品久久久久久久性| 成人国产av品久久久| 日本免费在线观看一区| 亚洲av一区综合| 人妻一区二区av| 日韩强制内射视频| 看黄色毛片网站| 少妇的逼好多水| 成人一区二区视频在线观看| 精品午夜福利在线看| 天堂网av新在线| 亚洲内射少妇av| 国产一区二区三区综合在线观看 | 水蜜桃什么品种好| av在线蜜桃| 毛片女人毛片| 永久网站在线| 午夜福利视频1000在线观看| 人妻制服诱惑在线中文字幕| 午夜免费鲁丝| 免费av观看视频| 欧美国产精品一级二级三级 | 我的女老师完整版在线观看| 乱码一卡2卡4卡精品| 亚洲怡红院男人天堂| 老师上课跳d突然被开到最大视频| 噜噜噜噜噜久久久久久91| 美女脱内裤让男人舔精品视频| 99热这里只有是精品在线观看| 最近最新中文字幕免费大全7| 亚洲性久久影院| 大片免费播放器 马上看| 人体艺术视频欧美日本| 高清欧美精品videossex| 亚洲电影在线观看av| 毛片一级片免费看久久久久| 成人亚洲精品一区在线观看 | 久久久久久九九精品二区国产| 最近手机中文字幕大全| av.在线天堂| 国产老妇伦熟女老妇高清| 亚洲欧美精品自产自拍| 日本欧美国产在线视频| 中文天堂在线官网| 亚洲欧美日韩另类电影网站 | 中文精品一卡2卡3卡4更新| 亚洲欧美日韩东京热| 黄色视频在线播放观看不卡| 菩萨蛮人人尽说江南好唐韦庄| 国产又色又爽无遮挡免| 午夜福利视频1000在线观看| 亚洲真实伦在线观看| 国产欧美日韩一区二区三区在线 | 日日撸夜夜添| 涩涩av久久男人的天堂| 爱豆传媒免费全集在线观看| 国产色爽女视频免费观看| 免费不卡的大黄色大毛片视频在线观看| 国产在线一区二区三区精| 日韩三级伦理在线观看| 黄色欧美视频在线观看| 亚洲内射少妇av| 日日摸夜夜添夜夜添av毛片| 最近中文字幕高清免费大全6| 老司机影院毛片| 日韩亚洲欧美综合| 国内揄拍国产精品人妻在线| 成人综合一区亚洲| 中国国产av一级| 网址你懂的国产日韩在线| 男插女下体视频免费在线播放| 永久网站在线| 欧美xxxx性猛交bbbb| 99精国产麻豆久久婷婷| 一本久久精品| 免费观看在线日韩| 国产老妇伦熟女老妇高清| 成人毛片60女人毛片免费| 丝袜脚勾引网站| av网站免费在线观看视频| 色婷婷久久久亚洲欧美| 免费高清在线观看视频在线观看| 色视频www国产| 久久精品久久久久久噜噜老黄| av在线老鸭窝| 丰满乱子伦码专区| 亚洲欧美日韩东京热| 亚洲综合精品二区| 亚洲欧美精品自产自拍| 亚洲精品成人av观看孕妇| 日本色播在线视频| 高清毛片免费看|