• <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| 美女xxoo啪啪120秒动态图| 长腿黑丝高跟| 99热网站在线观看| 成人三级黄色视频| 国内精品美女久久久久久| 波野结衣二区三区在线| 亚洲色图av天堂| 尾随美女入室| 国产在线男女| 国产伦精品一区二区三区四那| 国产精品一区二区免费欧美| 久久99热6这里只有精品| 男人的好看免费观看在线视频| 欧美黑人巨大hd| 午夜老司机福利剧场| 91久久精品国产一区二区成人| 国产一区二区激情短视频| 国产色爽女视频免费观看| 久久久久九九精品影院| 成人午夜高清在线视频| 日韩av在线大香蕉| 精品免费久久久久久久清纯| АⅤ资源中文在线天堂| 丰满人妻一区二区三区视频av| 亚洲美女视频黄频| 男女做爰动态图高潮gif福利片| 久久欧美精品欧美久久欧美| 动漫黄色视频在线观看| 国产亚洲精品久久久久久毛片| 欧美不卡视频在线免费观看| 欧美激情久久久久久爽电影| 亚洲,欧美,日韩| 黄色日韩在线| av在线蜜桃| 夜夜看夜夜爽夜夜摸| 国产极品精品免费视频能看的| 999久久久精品免费观看国产| 日本熟妇午夜| 日日摸夜夜添夜夜添小说| 久久久精品大字幕| 国内久久婷婷六月综合欲色啪| 级片在线观看| 久久精品国产亚洲av天美| 国产综合懂色| 欧美国产日韩亚洲一区| 尤物成人国产欧美一区二区三区| 国产高清三级在线| 精品一区二区三区av网在线观看| 精品99又大又爽又粗少妇毛片 | 黄色配什么色好看| 黄色一级大片看看| 国产精品国产高清国产av| 日韩在线高清观看一区二区三区 | 亚洲中文字幕日韩| 国产亚洲精品久久久久久毛片| 国产精品久久久久久久电影| 人妻制服诱惑在线中文字幕| 一区二区三区激情视频| 一区二区三区高清视频在线| 国产爱豆传媒在线观看| 亚洲自拍偷在线| 18禁黄网站禁片免费观看直播| 黄色女人牲交| 国产高清视频在线观看网站| 九色成人免费人妻av| av在线天堂中文字幕| 一区二区三区四区激情视频 | 我的老师免费观看完整版| 久久亚洲真实| 非洲黑人性xxxx精品又粗又长| 中文资源天堂在线| 一区福利在线观看| 久久久国产成人精品二区| 国产精品综合久久久久久久免费| 18禁裸乳无遮挡免费网站照片| 精品久久国产蜜桃| 中文字幕av成人在线电影| 国产精品一区二区免费欧美| 欧美黑人巨大hd| 欧美又色又爽又黄视频| 最近最新中文字幕大全电影3| 我要看日韩黄色一级片| netflix在线观看网站| 国产国拍精品亚洲av在线观看| 淫妇啪啪啪对白视频| 在线免费观看的www视频| 九色国产91popny在线| 色噜噜av男人的天堂激情| 色av中文字幕| 国产一区二区在线观看日韩| 一级av片app| 人人妻人人澡欧美一区二区| 日韩欧美国产在线观看| ponron亚洲| 狠狠狠狠99中文字幕| 97超级碰碰碰精品色视频在线观看| 亚洲精品亚洲一区二区| 国产麻豆成人av免费视频| 日韩精品青青久久久久久| 无人区码免费观看不卡| 我的老师免费观看完整版| 亚洲成人久久爱视频| 国产精品久久久久久亚洲av鲁大| 美女 人体艺术 gogo| av中文乱码字幕在线| 精品人妻熟女av久视频| 国产高清视频在线播放一区| 成人av在线播放网站| 91久久精品电影网| 成人av在线播放网站| 国产午夜精品久久久久久一区二区三区 | 国产 一区 欧美 日韩| 日韩欧美三级三区| av专区在线播放| av在线天堂中文字幕| 特级一级黄色大片| 69人妻影院| 欧美日韩综合久久久久久 | 男女视频在线观看网站免费| av在线亚洲专区| а√天堂www在线а√下载| 网址你懂的国产日韩在线| av在线天堂中文字幕| 搞女人的毛片| 狠狠狠狠99中文字幕| 国产白丝娇喘喷水9色精品| 久久午夜福利片| 我的老师免费观看完整版| 夜夜夜夜夜久久久久| 无遮挡黄片免费观看| 精品久久久久久成人av| 精品人妻视频免费看| 国内久久婷婷六月综合欲色啪| 深夜精品福利| 亚洲成av人片在线播放无| 免费大片18禁| 熟女电影av网| 熟女人妻精品中文字幕| 日本五十路高清| 欧美日韩亚洲国产一区二区在线观看| 日韩欧美一区二区三区在线观看| 精品久久久久久久久久免费视频| 色综合站精品国产| 大又大粗又爽又黄少妇毛片口| 深爱激情五月婷婷| 又紧又爽又黄一区二区| 一本一本综合久久| 2021天堂中文幕一二区在线观| 国产精品久久电影中文字幕| 亚洲精品色激情综合| 欧美日本视频| 国产色爽女视频免费观看| 九九爱精品视频在线观看| 久久久久国内视频| 久久中文看片网| 国产精品综合久久久久久久免费| 亚洲中文字幕日韩| 色哟哟哟哟哟哟| 免费观看精品视频网站| 国产精华一区二区三区| 日韩精品有码人妻一区| 国产精品国产高清国产av| 99久国产av精品| 69人妻影院| 在线观看一区二区三区| 少妇熟女aⅴ在线视频| 久久精品国产99精品国产亚洲性色| 我要看日韩黄色一级片| 久久精品国产自在天天线| 国国产精品蜜臀av免费| 69av精品久久久久久| 人妻少妇偷人精品九色| 久久精品国产99精品国产亚洲性色| 黄色欧美视频在线观看| 欧美精品啪啪一区二区三区| 日本成人三级电影网站| ponron亚洲| 国产精品福利在线免费观看| 波多野结衣高清作品| 一级毛片久久久久久久久女| 99精品久久久久人妻精品| 久久久久久久久久久丰满 | 精品日产1卡2卡| 三级毛片av免费| 在线播放国产精品三级| 精品久久久久久久人妻蜜臀av| 亚洲精华国产精华液的使用体验 | 黄色女人牲交| 成人午夜高清在线视频| 国产日本99.免费观看| 欧美一级a爱片免费观看看| 内射极品少妇av片p| 亚洲中文字幕日韩| 少妇被粗大猛烈的视频| 精品久久久久久久久av| 亚洲aⅴ乱码一区二区在线播放| 久久久久性生活片| 婷婷丁香在线五月| 免费一级毛片在线播放高清视频| 在线观看免费视频日本深夜| eeuss影院久久| 男女边吃奶边做爰视频| 嫁个100分男人电影在线观看| 日本撒尿小便嘘嘘汇集6| 精品国内亚洲2022精品成人| 午夜精品一区二区三区免费看| 亚洲av二区三区四区| 黄片wwwwww| 色综合亚洲欧美另类图片| 亚洲一级一片aⅴ在线观看| 噜噜噜噜噜久久久久久91| 春色校园在线视频观看| 99热这里只有精品一区| 波多野结衣巨乳人妻| 国产亚洲欧美98| 亚洲avbb在线观看| 欧美日韩亚洲国产一区二区在线观看| 欧美日本视频| 日韩人妻高清精品专区| 国产精品久久久久久精品电影| 亚洲性夜色夜夜综合| 俄罗斯特黄特色一大片| 日韩欧美精品免费久久| 女人被狂操c到高潮| 一本精品99久久精品77| 免费在线观看日本一区| 国产精品久久电影中文字幕| 欧美xxxx性猛交bbbb| 亚洲人成伊人成综合网2020| 亚洲图色成人| x7x7x7水蜜桃| 麻豆国产97在线/欧美| 日本精品一区二区三区蜜桃| 欧美在线一区亚洲| 在线观看美女被高潮喷水网站| 99riav亚洲国产免费| x7x7x7水蜜桃| 最近最新中文字幕大全电影3| 少妇熟女aⅴ在线视频| 中文字幕av成人在线电影| 国内揄拍国产精品人妻在线| 欧美性猛交╳xxx乱大交人| 熟女人妻精品中文字幕| 99热这里只有是精品50| 看十八女毛片水多多多| 国内少妇人妻偷人精品xxx网站| 国产精品,欧美在线| 精品免费久久久久久久清纯| 男女边吃奶边做爰视频| 成人午夜高清在线视频| 麻豆国产av国片精品| 亚洲欧美精品综合久久99| 啦啦啦观看免费观看视频高清| 国产精品亚洲一级av第二区| 少妇人妻一区二区三区视频| 淫秽高清视频在线观看| 狂野欧美白嫩少妇大欣赏| 女的被弄到高潮叫床怎么办 | 国产精品一区www在线观看 | 网址你懂的国产日韩在线| 少妇丰满av| 日本 欧美在线| 真实男女啪啪啪动态图| 亚洲成人久久性| 欧美一区二区国产精品久久精品| bbb黄色大片| 国产av一区在线观看免费| 亚洲成人免费电影在线观看| 一区二区三区激情视频| 神马国产精品三级电影在线观看| 日韩精品有码人妻一区| 国产黄色小视频在线观看| 97热精品久久久久久| 成年人黄色毛片网站| 丰满人妻一区二区三区视频av| 亚洲av一区综合| 午夜精品久久久久久毛片777| 久久精品国产亚洲av天美| 国产精品99久久久久久久久| 午夜福利高清视频| 亚洲乱码一区二区免费版| 亚洲美女黄片视频| 午夜免费成人在线视频| 免费av毛片视频| 日本a在线网址| 中国美女看黄片| 久久久久国内视频| 亚洲天堂国产精品一区在线| 亚洲经典国产精华液单| 亚洲乱码一区二区免费版| 成人特级av手机在线观看| 99精品久久久久人妻精品| 国产在视频线在精品| 一本一本综合久久| 欧美一级a爱片免费观看看| 成年版毛片免费区| 精华霜和精华液先用哪个| a级毛片免费高清观看在线播放| 国产精品女同一区二区软件 | 内射极品少妇av片p| 日日啪夜夜撸| 综合色av麻豆| 精品久久久久久久久久久久久| 亚洲 国产 在线| 精品免费久久久久久久清纯| 国产伦在线观看视频一区| 五月伊人婷婷丁香| 亚洲欧美日韩高清专用| 黄色配什么色好看| 天堂√8在线中文| 国产精品一区www在线观看 | 毛片一级片免费看久久久久 | 国产精品一区二区性色av| 黄色视频,在线免费观看| 亚洲第一电影网av| 亚洲美女黄片视频| 禁无遮挡网站| 人妻制服诱惑在线中文字幕| 亚洲18禁久久av| 国产亚洲欧美98| 全区人妻精品视频| 好男人在线观看高清免费视频| 日韩欧美免费精品| 午夜激情福利司机影院| 色综合婷婷激情| 少妇人妻一区二区三区视频| 国产乱人视频| 99精品久久久久人妻精品| 午夜激情欧美在线| 日韩精品青青久久久久久| 美女cb高潮喷水在线观看| 成年版毛片免费区| 日日干狠狠操夜夜爽| 亚洲无线观看免费| 免费看美女性在线毛片视频| 赤兔流量卡办理| 国产精品人妻久久久影院| 九九爱精品视频在线观看| 五月伊人婷婷丁香| 亚洲av不卡在线观看| 国产av在哪里看| 国产成人aa在线观看| 久久精品91蜜桃| 国模一区二区三区四区视频| 在线观看66精品国产| 69av精品久久久久久| 伦精品一区二区三区| 亚洲精品国产成人久久av| 国产av麻豆久久久久久久| 88av欧美| 日本一本二区三区精品| 波多野结衣高清作品| 国产精品一区二区免费欧美| 国产精品一区二区三区四区久久| 99热这里只有是精品50| 亚洲av成人av| 日本欧美国产在线视频| 又黄又爽又免费观看的视频| 国产一区二区亚洲精品在线观看| 精品99又大又爽又粗少妇毛片 | 色综合站精品国产| 欧美3d第一页| 成年女人看的毛片在线观看| 舔av片在线| 日韩欧美在线二视频| 国产高潮美女av| 一个人免费在线观看电影| 日本撒尿小便嘘嘘汇集6| 日日摸夜夜添夜夜添小说| 熟女电影av网| 中文字幕av在线有码专区| 1000部很黄的大片| 精品乱码久久久久久99久播| 最后的刺客免费高清国语| 国产av一区在线观看免费| 亚洲专区国产一区二区| 久久精品国产亚洲av香蕉五月| 亚洲黑人精品在线| 国产v大片淫在线免费观看| 欧美日韩国产亚洲二区| av在线观看视频网站免费| 女生性感内裤真人,穿戴方法视频| 亚洲五月天丁香| 欧美极品一区二区三区四区| 成年女人看的毛片在线观看| 欧美激情在线99| 国产亚洲精品综合一区在线观看| 国产精品久久久久久精品电影| 国产美女午夜福利| 波野结衣二区三区在线| 久久精品综合一区二区三区| 欧美极品一区二区三区四区| 欧美+亚洲+日韩+国产| 少妇猛男粗大的猛烈进出视频 | 国产在线精品亚洲第一网站| 国产黄a三级三级三级人| 99精品在免费线老司机午夜| 狂野欧美激情性xxxx在线观看| 麻豆一二三区av精品| 九九热线精品视视频播放| 18+在线观看网站| 综合色av麻豆| 91在线观看av| 亚州av有码| 麻豆精品久久久久久蜜桃| 久久香蕉精品热| 一夜夜www| 中文在线观看免费www的网站| 亚洲成人中文字幕在线播放| 国产淫片久久久久久久久| av在线亚洲专区| 国产黄a三级三级三级人| 尤物成人国产欧美一区二区三区| 亚洲av中文av极速乱 | 亚洲中文字幕日韩| 男女视频在线观看网站免费| 在线国产一区二区在线| 亚洲熟妇中文字幕五十中出| 亚洲男人的天堂狠狠| 一区二区三区激情视频| 最近中文字幕高清免费大全6 | 色哟哟哟哟哟哟| 国模一区二区三区四区视频| 一个人看视频在线观看www免费| 精品久久久久久久久av| 亚洲人成网站在线播| 国产精品久久久久久亚洲av鲁大| 好男人在线观看高清免费视频| 国产视频一区二区在线看| 一级毛片久久久久久久久女| 日韩一区二区视频免费看| 日日撸夜夜添| av黄色大香蕉| 国产国拍精品亚洲av在线观看| 日韩欧美在线乱码| 69人妻影院| .国产精品久久| 日本欧美国产在线视频| av天堂在线播放| 婷婷丁香在线五月| 亚洲一级一片aⅴ在线观看| 国内精品久久久久精免费| 亚洲无线在线观看| 免费看光身美女| 一区二区三区免费毛片| 精品久久久久久成人av| 久久精品国产清高在天天线| 大又大粗又爽又黄少妇毛片口| 九九热线精品视视频播放| 国产成人一区二区在线| 日本 欧美在线| 亚洲中文字幕一区二区三区有码在线看| 51国产日韩欧美| 国产午夜精品久久久久久一区二区三区 | 国产精品免费一区二区三区在线| 99精品在免费线老司机午夜| 少妇裸体淫交视频免费看高清| 亚洲av二区三区四区| 一本一本综合久久| 国产三级在线视频| 亚洲国产精品久久男人天堂| 日本精品一区二区三区蜜桃| 亚洲国产日韩欧美精品在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产精品久久久久久久电影| 热99re8久久精品国产| 少妇人妻一区二区三区视频| 国产综合懂色| 国国产精品蜜臀av免费| 亚洲国产精品合色在线| 亚洲四区av| 精品人妻1区二区| 午夜福利成人在线免费观看| 美女高潮喷水抽搐中文字幕| 18+在线观看网站| 69av精品久久久久久| 亚洲第一区二区三区不卡| 99热6这里只有精品| 亚洲欧美日韩无卡精品| 免费看av在线观看网站| 无遮挡黄片免费观看| 最新在线观看一区二区三区| 一进一出抽搐gif免费好疼| 午夜精品在线福利| 亚洲人成网站在线播放欧美日韩| 久久亚洲真实| 精品福利观看| 色播亚洲综合网| 成人国产综合亚洲| 少妇猛男粗大的猛烈进出视频 | 国产一区二区在线观看日韩| 成人二区视频| 午夜影院日韩av| 桃色一区二区三区在线观看| 老熟妇仑乱视频hdxx| av女优亚洲男人天堂| 男女啪啪激烈高潮av片| 亚洲久久久久久中文字幕| 免费高清视频大片| 亚洲性久久影院| 一区二区三区免费毛片| 欧美+亚洲+日韩+国产| 亚洲人成网站在线播| 久久精品影院6| 一个人看视频在线观看www免费| 高清毛片免费观看视频网站| 国产av在哪里看| 久久久久精品国产欧美久久久| 久久久久久久午夜电影| 成年版毛片免费区| 可以在线观看的亚洲视频| 搡老岳熟女国产| 淫妇啪啪啪对白视频| 制服丝袜大香蕉在线| 国产老妇女一区| 韩国av在线不卡| 日韩欧美一区二区三区在线观看| 窝窝影院91人妻| 亚洲欧美日韩卡通动漫| av在线蜜桃| 男人狂女人下面高潮的视频| 免费不卡的大黄色大毛片视频在线观看 | 国产精品久久电影中文字幕| 午夜a级毛片| 女的被弄到高潮叫床怎么办 | 女的被弄到高潮叫床怎么办 | 午夜视频国产福利| 国产av不卡久久| 亚洲在线自拍视频| 88av欧美| 中文字幕久久专区| 亚洲av日韩精品久久久久久密| 国产亚洲精品av在线| 国产黄色小视频在线观看| 91狼人影院| 国产精品女同一区二区软件 | 日本免费a在线| 日韩国内少妇激情av| 亚洲av熟女| 在线a可以看的网站| 国产精华一区二区三区| 午夜福利在线观看免费完整高清在 | 亚洲熟妇熟女久久| 少妇高潮的动态图| 一进一出抽搐动态| 久久婷婷人人爽人人干人人爱| 中文字幕高清在线视频| 亚洲五月天丁香| 亚洲精品在线观看二区| 欧美极品一区二区三区四区| 亚洲av美国av| 国产精品99久久久久久久久| 好男人在线观看高清免费视频| 精品久久久久久久久亚洲 | 两个人视频免费观看高清| 色哟哟哟哟哟哟| 成人亚洲精品av一区二区| 国产一区二区亚洲精品在线观看| 亚洲专区国产一区二区| 桃色一区二区三区在线观看| 亚洲精品色激情综合| 内射极品少妇av片p| 最近在线观看免费完整版| 一级a爱片免费观看的视频| 欧美一级a爱片免费观看看| 久久6这里有精品| ponron亚洲| 亚洲精品乱码久久久v下载方式| 狂野欧美激情性xxxx在线观看| 国产熟女欧美一区二区| 成年女人毛片免费观看观看9| 少妇的逼好多水| 欧美人与善性xxx| 国产精品久久视频播放| 免费在线观看影片大全网站| 国内精品久久久久久久电影| 亚洲国产精品sss在线观看| 人妻久久中文字幕网| 久久婷婷人人爽人人干人人爱| 亚洲人成网站高清观看| 欧美性猛交黑人性爽| 亚洲最大成人中文| 午夜福利在线观看吧| 国产精品一区二区性色av| 久久精品国产亚洲网站| 亚洲精品日韩av片在线观看| 九色国产91popny在线| 91在线精品国自产拍蜜月| 22中文网久久字幕| 久久精品国产鲁丝片午夜精品 | av天堂中文字幕网| 桃色一区二区三区在线观看| 国产国拍精品亚洲av在线观看| h日本视频在线播放| 桃色一区二区三区在线观看| 国产av不卡久久| 在线观看66精品国产| 国产单亲对白刺激| 亚洲精品粉嫩美女一区| 国产老妇女一区| 亚洲欧美精品综合久久99| 国产一区二区亚洲精品在线观看| 国产亚洲精品久久久久久毛片| 婷婷丁香在线五月| 露出奶头的视频| 国产午夜福利久久久久久| 九九爱精品视频在线观看| 国产精品女同一区二区软件 | 亚洲欧美精品综合久久99| a级毛片a级免费在线| 日本-黄色视频高清免费观看|