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

    礦物微量元素組成用于火成巖構(gòu)造背景判別

    2016-12-12 08:38:58趙振華中國(guó)科學(xué)院廣州地球化學(xué)研究所廣東廣州510640
    大地構(gòu)造與成礦學(xué) 2016年5期
    關(guān)鍵詞:單斜火成巖輝石

    趙振華(中國(guó)科學(xué)院 廣州地球化學(xué)研究所, 廣東 廣州 510640)

    礦物微量元素組成用于火成巖構(gòu)造背景判別

    趙振華
    (中國(guó)科學(xué)院 廣州地球化學(xué)研究所, 廣東 廣州 510640)

    火成巖中的礦物, 特別是穩(wěn)定副礦物的微量元素組成可用于火成巖構(gòu)造–巖漿背景的判別。產(chǎn)于不同構(gòu)造背景火成巖中單顆粒鋯石原位微量元素含量的準(zhǔn)確測(cè)定, 積累了大量資料。本文介紹了用鋯石微量元素含量的二維核密度分布投影, 構(gòu)建了 U/Yb-Nb/Yb、Sc/Yb-Nb/Yb和 U/Yb-Sc/Yb構(gòu)造背景判別圖解, 可以區(qū)分洋中脊、地幔柱影響的洋島及俯沖有關(guān)的弧(大陸弧)等不同構(gòu)造背景形成的火成巖。單斜輝石、尖晶石、黑云母和角閃石的微量元素或主量元素組成同樣可用于構(gòu)造背景識(shí)別。

    構(gòu)造背景判別; 微量元素; 鋯石; 單斜輝石; 尖晶石; 黑云母; 角閃石

    盡管用巖石微量元素識(shí)別用于構(gòu)造背景的判別圖存在較多爭(zhēng)議, 但用微量元素判別基性巖漿巖構(gòu)造背景仍取得了較明顯的進(jìn)展(Vermeesch, 2006a, 2006b; 趙振華, 2007; Pearce 2014; Li et al., 2015)。主要包括以下三個(gè)方面:

    第一是大量、精確的主、微量元素?cái)?shù)據(jù)的積累。隨著電子探針?lè)治黾夹g(shù), 特別是激光等離子體質(zhì)譜(LA-ICP-MS)技術(shù)的快速發(fā)展, 基性巖全巖(包括單顆粒礦物, 特別是鋯石原位)微量元素含量的準(zhǔn)確測(cè)定積累了大量數(shù)據(jù), 建立了典型巖石類型的數(shù)據(jù)庫(kù),如 PETDB和 GEOROC數(shù)據(jù)庫(kù)(網(wǎng)址: http: //www. earthchem.org/petdb和 http: //www.rankinsider.com/ georoc.mpch-mainz.gwdg.de), 大大拓寬了全巖微量元素組成判別火成巖形成構(gòu)造背景方法的適應(yīng)范圍, 特別是為遭受蝕變巖石的構(gòu)造背景恢復(fù)提供了新的工具(Grimes et al., 2007, 2015; Gonzalez-Jimenez et al., 2012)。

    第二是更有效地運(yùn)用多變量統(tǒng)計(jì)分析方法。在發(fā)表的有關(guān)玄武巖構(gòu)造背景的傳統(tǒng)判別圖解中, 不同構(gòu)造背景分類之間界線的確定往往是用兩個(gè)或三個(gè)變量數(shù)據(jù), 不是數(shù)值化的, 在統(tǒng)計(jì)上不嚴(yán)格, 由此劃分的界線常常是不準(zhǔn)確的、或有較大重疊。例如, 為了將玄武巖構(gòu)造背景的判別分析建立在更嚴(yán)格的基礎(chǔ)上, Vermeesch (2006a, 2006b)收集了45種共756件已知構(gòu)造背景的、有全巖主、微量元素定量分析數(shù)據(jù)的海洋玄武巖樣品(洋島玄武巖259件、洋中脊玄武巖241件、島弧玄武巖256件), 建立了分類樹。分類和回歸樹(classification and regression tree, CART)是用一個(gè)分段常值函數(shù)采用了更嚴(yán)格的統(tǒng)計(jì)方法——線性判別分析和二次判別分析。所選擇的數(shù)據(jù)在統(tǒng)計(jì)分析前均轉(zhuǎn)換為對(duì)數(shù)比值(不是直接用該成分的含量)。判別分析是假定多變量正態(tài)性(multivariate normality), 如果所有的構(gòu)造背景分類具有相同的協(xié)方差結(jié)構(gòu), 則所確定的分類界線是線性的, 稱為線性判別分析(linear discriminant analysis, LDA), 與此相反則為二次判別分析(quadratic discriminant analysis, QDA)。用上述方法產(chǎn)生了 14190個(gè)三元判別圖解, 對(duì)它們的窮舉分析(exhaustive exploration)獲得了最好的線性和二次判別圖解, 分別是不活動(dòng)元素 Ti-V-Sc和 Ti-V-Sm三元圖解。上述方法減少了判別分析的錯(cuò)判, 判別函數(shù)是線性結(jié)合, 它使不同構(gòu)造背景分類之間相對(duì)于每一類型內(nèi)部的變化最大化。Belousova et al. (2002)根據(jù)鋯石的稀土和微量元素含量和組成, 采用遞歸分割法(recursive partitioning, 簡(jiǎn)稱RP法), 該方法關(guān)鍵是選擇節(jié)點(diǎn), 或稱根節(jié)點(diǎn)、終端節(jié)點(diǎn)、葉節(jié)點(diǎn), 進(jìn)而建立決策樹或決定樹。構(gòu)建判別不同類型巖漿巖的分類樹(classification and regression trees)。本文介紹的不同構(gòu)造?巖漿背景鋯石微量元素判別圖則是用其微量元素的二維核密度分布投影(詳見(jiàn)下述)。與上述分類樹類似, 對(duì)于一些特殊類型的火成巖建立了一組圖解, 使用這些判別圖解時(shí)需按一定流程(flow-chart)進(jìn)行, 如對(duì)于鉀質(zhì)火成巖(Müller and Groves, 1997)。

    第三是用礦物, 特別是穩(wěn)定副礦物的、而不是全巖的微量元素組成構(gòu)建構(gòu)造?巖漿背景判別圖解。

    本文重點(diǎn)介紹火成巖中一些造巖礦物和副礦物(特別是鋯石)微量元素組成在構(gòu)造背景判別中的應(yīng)用。

    1 鋯 石

    鋯石是倍受關(guān)注的遍在性礦物, 從地殼巖石到地幔包體、月巖及隕石中都有鋯石存在。鋯石還是化學(xué)上很穩(wěn)定的礦物, 在風(fēng)化、搬運(yùn)過(guò)程, 甚至在變質(zhì)和深熔過(guò)程都可保持化學(xué)上的穩(wěn)定。因此, 火成巖中鋯石U-Pb被廣泛用于定年, 特別是激光等離子體質(zhì)譜的單顆粒礦物定年技術(shù)的發(fā)展, 使鋯石成為巖漿巖U-Pb定年的最重要對(duì)象, 積累了大量年齡數(shù)據(jù)。而在鋯石用于U-Pb定年的同時(shí), 其單顆粒原位微量元素含量的準(zhǔn)確測(cè)定更提供了探討晶出鋯石的巖漿源區(qū)和成巖過(guò)程的重要地球化學(xué)信息。

    1.1鋯石的稀土及微量元素地球化學(xué)特點(diǎn)

    強(qiáng)烈富集重稀土和明顯的 Ce正異常是巖漿鋯石的典型特征, 這種特點(diǎn)主要受控于晶出鋯石的巖漿巖成分以及鋯石的稀土和微量元素分配系數(shù)。

    不同類型巖漿巖中鋯石的稀土和微量元素含量和組成有明顯區(qū)別, 如: 金伯利巖中的鋯石稀土總量低(<50×10-6), 重稀土富集程度也較低, (Yb/Sm)N= 3~30; 但花崗巖中鋯石稀土含量可達(dá)百分之幾, 重稀土富集程度很高, (Yb/Sm)N>100。因此可根據(jù)鋯石稀土元素組成特點(diǎn)區(qū)分源區(qū)的特征(表1)。

    鋯石稀土分配系數(shù)明顯與形成鋯石的巖漿成分有關(guān), 硅不飽和堿性巖中鋯石的稀土分配系數(shù)最高,而基性巖最低, 中性巖、酸性巖居中(Hanchar and Westrenen, 2007)。

    表1 不同類型火成巖鋯石的稀土元素組成Table 1 REE composition of zircon from different type of igneous rocks

    稀土元素在鋯石(ZrSiO4)中的分配系數(shù)還受控于其晶體化學(xué)(離子半徑及電荷)特點(diǎn)。鋯石中稀土元素的進(jìn)入是基于與鋯石等結(jié)構(gòu)的“磷釔礦式”(YPO4)置換:

    稀土元素中重稀土元素Lu3+離子半徑為0.0977 nm, Y3+為 0.1019 nm; 輕稀土元素 La的離子半徑為0.116 nm, 而8配位的Zr4+離子半徑為0.084 nm, 因此, 重稀土元素離子與 Zr4+離子半徑更接近, 更容易進(jìn)入鋯石晶格, 導(dǎo)致鋯石中強(qiáng)烈富集重稀土, 這明顯反映在分配系數(shù)上。Hanchar and Westrenen (2007) 綜合了自然界火山巖和熔體包裹體樣品及由高溫高壓實(shí)驗(yàn)研究所獲得的鋯石的稀土元素分配系數(shù)(圖1), 其特點(diǎn)是: 重稀土(HREE)分配系數(shù)明顯高于輕稀土(LREE), 在一個(gè)樣品中 DLu比 DLa高 1~5個(gè)數(shù)量級(jí); 不同類型巖石之間鋯石稀土分配系數(shù)差別大, 對(duì)于輕稀土變化達(dá) 4個(gè)數(shù)量級(jí), 而重稀土變化相對(duì)較小, 為2個(gè)數(shù)量級(jí)。對(duì)于其他微量元素, 重要的是U、Th、Nb和Sc。U4+、Th4+與Zr4+電價(jià)相同, 八配位Th4+離子半徑為1.0 ?, 八配位的U更大,離子半徑為1.04 ?, 因此鋯石對(duì)它們的相容性比Hf4+低。Hf、U和 Th與鋯石可形成等結(jié)構(gòu)的固溶體: HfSiO4、USiO4和ThSiO4, 但這種固溶體是非理想的、有較大不混溶間斷(Finch and Hanchar, 2003; Geisler et al., 2007)。鋯石對(duì)U和Yb的分配系數(shù)相近, 分別為254和278。Sc和Nb進(jìn)入鋯石的置換方式是:

    上述不同稀土元素分配系數(shù)之間的差異可從其離子半徑及相關(guān)晶格變化得到解釋。在鋯石中稀土元素分配系數(shù)對(duì)離子半徑的 Onuma圖解中, 從 Lu到La構(gòu)成簡(jiǎn)單拋物線(圖 1), 該曲線頂部是最佳被置換的 Zr4+的位置, 從曲線上稀土元素的排列順序可以看出, 重稀土特別是Lu最先置換Zr4+。從圖上還可以看出Ce明顯偏離拋物線, 反映在鋯石Ce的分配系數(shù)DCe明顯高于其相鄰輕稀土元素La和Pr或Nd, 形成顯著Ce正異常。這是由于在氧化條件下 Ce3+氧化成 Ce4+,它與 Zr4+電價(jià)相同, 在置換時(shí)不需電價(jià)補(bǔ)償, 而且, Ce4+的半徑(6配位為 0.087 nm, 8配位為 0.097 nm)比Ce3+更近于Zr4+, 因而比其他LREE更容易進(jìn)入鋯石晶格而形成 Ce正異常。此外, DCe偏離拋物線的程度可作為鋯石生長(zhǎng)過(guò)程中氧化狀態(tài)的指標(biāo)。

    圖1 鋯石–熔體稀土分配系數(shù)與離子半徑和熔體成分關(guān)系的Onuma圖解(引自Hanchar and Westrenen, 2007)Fig.1 Onuma diagrams showing zircon-melt REE partition coefficients as a function of REE ion radius

    1.2不同源區(qū)巖漿巖鋯石的稀土與微量元素組成

    由于鋯石對(duì)重稀土元素的強(qiáng)烈選擇性, 很多巖石中鋯石的稀土含量發(fā)生明顯重疊, 因此難以區(qū)分開(kāi)不同源區(qū)的鋯石。由上述, U與Yb在鋯石中分配系數(shù)相近, 因此, 鋯石的 U/Yb可反映鋯石晶出時(shí)熔體的特點(diǎn)。但是, Nb和Sc在鋯石中的分配系數(shù)變化范圍很大(DNb: 0.2~21; DSc: 2~70), 可達(dá)≥10倍。因此, Grimes et al. (2007, 2015)用鋯石中微量元素U/Nb比值及它們與稀土元素的比值(U/Yb, U/Nb或 Nb/Yb)進(jìn)行不同類型鋯石比較, 這可以使由溫度造成的分配系數(shù)變化導(dǎo)致的濃度改變最小化。對(duì)采自緩慢擴(kuò)張中的大西洋和南印度洋洋脊的 36個(gè)輝長(zhǎng)巖中的300個(gè)鋯石顆粒, 以及1700多個(gè)采自大陸太古宙到顯生宙花崗巖類、金伯利巖和阿拉斯加Talkeetna 島弧石英閃長(zhǎng)巖和英云閃長(zhǎng)巖巖墻的鋯石進(jìn)行了微量元素和 REE元素分析(用 SHRIMPRG離子探針?lè)治?。用 U-Yb、U/Yb-Hf或Y; Th-Yb; Th/Yb-Hf 或Y、P對(duì)上述鋯石作圖(圖 2), 可以將洋殼輝長(zhǎng)巖的鋯石與大陸巖石中的鋯石明顯區(qū)分開(kāi)。雖然對(duì)于同一樣品中的U、Th、Yb的絕對(duì)含量可以有數(shù)量級(jí)的變化, 但洋殼巖石的鋯石 U/Yb比值平均為 0.18, 大陸巖石鋯石平均為 1.07, 金伯利巖為 2.1。這種變化與含這些鋯石的全巖 U/Yb比值一致, 洋中脊玄武巖 MORB平均為 0.01~0.1 (Klein, 2003), 大陸殼平均為0.7(Rudnick and Gao, 2003), 金伯利巖為4~6(Farmer, 2003)。阿拉斯加島弧巖石的鋯石與洋殼巖石的鋯石明顯不同, 但與大陸巖石鋯石重疊, 一些海洋島弧熔巖全巖與MORB明顯重疊, 由此推測(cè)大陸溢流玄武巖和大洋島弧玄武巖中的鋯石與洋殼巖石鋯石相似。

    由上述可見(jiàn), 在構(gòu)筑判別圖解的參數(shù)選擇中, U/Yb是很重要的參數(shù), MORB中鋯石低的U/Yb比值(<0.10)反映了它們的母巖漿及MORB地幔長(zhǎng)期的不相容元素的虧損, 它明顯不同于其他構(gòu)造背景的鋯石。但U/Yb>0.1則不是唯一的, 它包括了多種構(gòu)造背景的鋯石。因此, 需采用包括 U-Nb-Sc-Yb-Gd-Ce等多種陽(yáng)離子的雙變量圖解(Grimes et al., 2015)。例如, 在 U/Yb-Nb/Yb圖解中(實(shí)質(zhì)是鋯石的 U/Nb比值), 包含了多數(shù)MORB、洋島和金伯利巖鋯石的陰影區(qū), 形成了明顯的地幔鋯石排列, 其上部邊界的Nb/Yb和U/Yb值分別為0.004, 0.02和1.0, 10。而巖漿弧和后碰撞大陸鋯石則形成了排列之外的、高U/Yb比值的與之平行的巖漿弧排列(圖3)。

    圖2 區(qū)分大陸和海洋鋯石的U/Yb-Hf (a)和U/Yb-Y (b)圖解(Grimes et al., 2007)Fig.2 Plots of U/Yb vs. Hf (a), and U/Yb vs. Y (b) for discrimination of continental and ocean crust zircons

    圖3 用U/Nb代表火成巖鋯石的構(gòu)造?巖漿源(Grimes et al., 2015)Fig.3 U/Nb proxy for tectono-magmatic source of igneous zircon

    受地幔柱影響的洋島鋯石與MORB的明顯區(qū)別是U/Yb、U/Nb和Nb/Sc比值高; 對(duì)于U/Yb>0.1的鋯石, 高 Sc/Yb可以將弧背景與低 Sc/Yb的受地幔柱影響的洋島區(qū)分開(kāi)。鋯石中Sc的富集可能反映了熔體分異過(guò)程中角閃石的卷入, 進(jìn)而反映了拉斑質(zhì)與鈣堿質(zhì)巖漿系統(tǒng)源區(qū)成分的差異。

    1.3鋯石稀土與微量元素組成對(duì)不同構(gòu)造?巖漿背景的判別

    基于上述火成巖中鋯石微量元素與其源區(qū)及成巖過(guò)程的密切和復(fù)雜關(guān)系, Grimes et al. (2015)對(duì)鋯石微量元素組成開(kāi)展了更系統(tǒng)的研究, 收集了5438個(gè)普通構(gòu)造?巖漿背景產(chǎn)出的火成巖中鋯石SHRIMP-RG分析的微量元素?cái)?shù)據(jù), 包括洋中脊 MOR(大西洋中脊MAR; 西南印度洋脊SWIR; 東太平洋洋隆EPRS)、巖漿弧、大陸碰撞帶和地幔柱影響的洋島; 時(shí)代主要為中生代?新生代。用數(shù)理統(tǒng)計(jì)方法對(duì)這些資料進(jìn)行了計(jì)算, 統(tǒng)計(jì)參數(shù)包括平均值、中值、均方差(1σ)、偏態(tài)、峰態(tài)、第一四分位數(shù)(quartile)和第三四分位數(shù)。表2中給出了不同構(gòu)造背景火成巖中鋯石微量元素含量及U、Th、Ce、Gd、Nb和Sc與Yb的比值及U/Nb比值的平均值、中值、方差和樣品數(shù)。

    根據(jù)上述統(tǒng)計(jì)資料, Grimes et al. (2015)構(gòu)筑了不同構(gòu)造?巖漿背景的鋯石微量元素判別圖解, 這些圖解包括U/Yb-Nb/Yb、Sc/Yb-Nb/Yb和 U/Yb-Sc/Yb (圖4), 圖中坐標(biāo)均為對(duì)數(shù)值(lg)。圖解中陰影區(qū)是所收集的不同構(gòu)造?巖漿背景鋯石微量元素的二維核密度分布投影。等值線分別為 50%、80%、90%和95%, 它們代表了該等值線內(nèi)樣品所占的比例。該密度投影是用二維核密度函數(shù)(kde2d)的 R程序, 即用軸平行(axis-aligned)二元正態(tài)核(bivariated normal kernel)計(jì)算, 用方格網(wǎng)(square grid)評(píng)價(jià)的二維核密度。所產(chǎn)生的構(gòu)造背景圖解中, 洋中脊、地幔柱影響的洋島及俯沖有關(guān)的弧背景明顯分布在不同區(qū)域,每一種背景中的近 80%~90%數(shù)據(jù)都投影于圖中所限定的相應(yīng)區(qū)域中。

    圖4 不同構(gòu)造?巖漿背景鋯石微量元素地球化學(xué)標(biāo)志的密度分布圖(Grimes et al., 2015)Fig.4 Density distribution plots of tectono-magmatic setting based on trace elements in igneous zircon

    在使用上述圖解時(shí)應(yīng)注意的是, Grimes et al. (2015)用鋯石微量元素成分判別構(gòu)造?巖漿背景的圖解是基于現(xiàn)代到中生代的火成巖中的鋯石, 而且對(duì)古生代火成巖的鋯石也做了檢驗(yàn), 其結(jié)果基本一致,但如果擴(kuò)展到前寒武紀(jì)應(yīng)當(dāng)非常小心。此外, 用于構(gòu)建構(gòu)造背景判別圖解的鋯石微量元素資料均來(lái)自SHRIMP-RG分析方法, 不同實(shí)驗(yàn)室之間的數(shù)據(jù)應(yīng)進(jìn)行所分析標(biāo)樣的校正。判別圖中不同區(qū)域之間界線的絕對(duì)值可能有些移動(dòng), 但其空間關(guān)系不會(huì)變化。此外, 應(yīng)綜合采用有 U、Yb、Nb、Sc、Ce和Gd建立的多變量圖解, 不要僅依據(jù)單一圖解。在鋯石微量元素分析中, 應(yīng)特別注意Nb、Sc的分析必須準(zhǔn)確定量測(cè)定, 因?yàn)檫@兩個(gè)元素容易受占鋯石 50%多的Zr的同量異位效應(yīng)的干擾。

    2 單斜輝石

    單斜輝石的化學(xué)成分隨著其寄生火成巖巖石的化學(xué)成分而變化, 對(duì)于單斜輝石的斑晶更為明顯。不同類型玄武巖的單斜輝石斑晶成分之間的差別比相應(yīng)的基質(zhì)成分之間的差別更大。該特征被用作識(shí)別不同構(gòu)造背景的判別要素, 同時(shí)也提供了判別蝕變巖石構(gòu)造背景的可靠手段, 因?yàn)槲g變巖石的單斜輝石斑晶核心的化學(xué)成分可能沒(méi)有變化(Rollinson, 1993)。Leterrier et al. (1982)對(duì) Nisbet and Pearce (1977)提出的單斜輝石MnO-TiO2-Na2O含量判別圖解進(jìn)行了修改完善, 構(gòu)筑了新的 Ti-Cr-Ca-Al-Na投影圖解, 其圖解是根據(jù)一系列更大數(shù)據(jù)庫(kù)為基礎(chǔ)的Ti、Cr、Ca、Al和Na資料而完成的(圖5)。Leterrier et al. (1982)提出的三個(gè)圖解可以區(qū)別堿性玄武巖、擴(kuò)張中心拉斑玄武巖和島弧玄武巖。單斜輝石分析數(shù)據(jù)按六個(gè)氧進(jìn)行晶體化學(xué)計(jì)算, 換算出陽(yáng)離子數(shù),并且只用那些Ca原子數(shù)大于0.5(每個(gè)單斜輝石分子中)的分析數(shù)據(jù)。為了提供計(jì)量比例, 進(jìn)行了 Fe3+和Fe2+的計(jì)算(見(jiàn)Cameron and Papike, 1981)。第一個(gè)圖解(圖 5a)利用 Ti-(Ca+Na)投影, 區(qū)分堿性玄武巖(A: 大洋島和大陸堿性玄武巖)和拉斑玄武巖及鈣堿性玄武巖(T)。第二個(gè)圖解(圖 5b)是(Ti+Cr)-Ca的投影, 區(qū)分非堿性玄武巖及非造山玄武巖(D: MORB, 大洋島拉斑玄武巖和弧后盆地拉斑玄武巖)與火山弧玄武巖(O)。第三個(gè)圖解(圖5c)是Ti-Al總的投影圖解, 區(qū)分火山弧玄武巖和鈣堿性玄武巖(C)與島弧拉斑玄武巖(I)。

    這些圖解不能只用一個(gè)單斜輝石的分析數(shù)據(jù); Leterrier et al. (1982)推薦的最少不得少于10個(gè)分析數(shù)據(jù)。如用20個(gè)分析數(shù)據(jù)投影到圖解上將會(huì)獲得更可靠的結(jié)果。此方法甚至可以應(yīng)用于綠片巖相的變基性巖,其中的單斜輝石成分可能經(jīng)受了變質(zhì)反應(yīng)的調(diào)整。

    Le Bas (1962)曾用單斜輝石的四面體中Al所占比例Al總與Ti含量關(guān)系劃分鎂鐵巖漿系列。Loucks (1990)認(rèn)為, 單斜輝石的 Al總與Ti關(guān)系可用于區(qū)分不同構(gòu)造環(huán)境的鎂鐵–超鎂鐵巖石, 其依據(jù)是在非造山區(qū)(如洋中脊, 弧后盆地、大陸裂谷、熱點(diǎn))的拉斑玄武質(zhì)和硅不飽和巖石中, 單斜輝石晶格結(jié)合的Al大部分是通過(guò)電荷補(bǔ)償VIMgIVSi2=VITiIVAl2, 以CaTiAl2O6分子形式進(jìn)入單斜輝石。這種置換方式由于在熔體中 Ti的活度比Si高, 在硅不飽和鎂鐵-超鎂鐵巖石中更為典型。在匯聚板塊邊緣(拉斑玄武巖、鈣堿性輝長(zhǎng)巖和玄武巖), 鎂鐵質(zhì)火成巖的特點(diǎn)是大部分Al以CaFe3+AlSiO6分子結(jié)合在單斜輝石中,置換方式是VIMgIVSi2?VIFe3+IVAl, 這是由于在島弧軸部水和氧的逸度高。將輝長(zhǎng)巖和超鎂鐵巖堆積巖的輝石中四面體Al對(duì)八面體中的Ti作圖, 可見(jiàn)與島弧有關(guān)的堆積巖中的單斜輝石 Al/Ti比值比與裂谷有關(guān)的拉斑玄武巖中的單斜輝石高, 約為其兩倍。因此,單斜輝石 Al/Ti比值可用于區(qū)分造山帶中蛇綠巖質(zhì)和非蛇綠巖質(zhì)的鎂鐵–超鎂鐵外來(lái)體。我們用這種方法對(duì)新疆克拉瑪依地區(qū)晚石炭世哈圖玄武巖的形成構(gòu)造背景進(jìn)行了討論(Tang et al., 2012, 圖6), 該玄武巖具有與大洋高原玄武巖相似的地球化學(xué)特點(diǎn), 其輝石的Al總對(duì)TiO2含量作圖分布于島弧和裂谷之間區(qū)域, 并從島弧趨勢(shì)偏向裂谷。結(jié)合該區(qū)早石炭世典型島弧環(huán)境及與哈圖玄武巖同時(shí)期的埃達(dá)克巖的存在, 認(rèn)為哈圖玄武巖形成于島弧和伸展的雙重背景, 與洋內(nèi)弧中洋脊俯沖形成的板片窗有關(guān)。

    圖5 玄武巖的單斜輝石斑晶對(duì)構(gòu)造背景的Ti-(Ca+Na), (Ti+Cr)-Ca和Ti-Al總判別圖解(據(jù)Leterrier et al., 1982)Fig.5 Discrimination diagrams for clinopyroxene phenocrysts in basalts

    3 尖 晶 石

    尖晶石的化學(xué)成分對(duì)于其母巖漿的地球化學(xué)非常靈敏, 是其母巖漿形成的構(gòu)造背景的函數(shù), 因此,尖晶石的Cr、Al、Fe、Mg和Ti的比值廣泛用于識(shí)別其母巖漿性質(zhì)和鉻鐵巖形成的構(gòu)造背景。然而,不同構(gòu)造背景形成的鉻鐵巖中尖晶石在 Cr-Al-Fe2+, (Cr/(Cr+Al))/(Mg/(Mg+Fe2+)), TiO2-Al2O3, TiO2-Cr2O3等圖解中出現(xiàn)大量重疊。Gonzalez-Jimenez et al. (2012)用LA-ICP-MS對(duì)已知構(gòu)造背景的鉻鐵巖中的尖晶石進(jìn)行微量元素分析, 并以MORB中的尖晶石為標(biāo)準(zhǔn)進(jìn)行作圖, 發(fā)現(xiàn)可以清楚地將科馬提巖、層狀鎂鐵質(zhì)侵入體及蛇綠巖區(qū)分開(kāi)。蛇綠巖中的尖晶石以富Ni、Mg而與虧損Ni、Mg的科馬提巖、層狀鎂鐵質(zhì)侵入體中的尖晶石區(qū)分開(kāi), 層狀鎂鐵質(zhì)侵入體中的尖晶石則明顯富集Sc、V而與科馬提巖蛇綠巖中的尖晶石相區(qū)分(圖7)。

    圖6 克拉瑪依哈圖玄武巖單斜輝石Al總-TiO2與構(gòu)造背景關(guān)系(據(jù)Tang et al., 2012修改自Loucks, 1990)Fig.6 Plot of Al總vs. TiO2in clinopyroxene for tectonic settings

    4 黑 云 母

    在Debon and Le Fort (1983)提出的侵入體化學(xué)?礦物分類中, 提出了特征礦物(characteristic mineral diagram)圖解, 圖解由參數(shù) A=Al?(K+Na+2Ca); B= (Fe+Mg+Ti)和礦物分類組成。該圖解包含了對(duì)火成巖含Al特點(diǎn)和對(duì)其礦物學(xué)的解釋, 它將火成巖劃分為鈣鐵鎂質(zhì)(cafemic, 主要或全部源于地幔)、Al質(zhì)(主要或全部源于陸殼深熔)及鈣鐵鎂質(zhì)?Al質(zhì)過(guò)渡三種類型。這種特征礦物圖解包含了巖石所固有的構(gòu)造背景信息。據(jù)此, 對(duì)堿性、鈣堿性和過(guò)鋁質(zhì)巖漿中的325個(gè)黑云母典型成分進(jìn)行了分析(Abdel-Rahman, 1994), 其典型特點(diǎn)是堿性非造山系列巖漿巖中黑云母富鐵, 為鐵云母, FeO*/MgO平均值為7.04; 過(guò)鋁質(zhì)系列巖漿巖的黑云母親鐵, FeO*/MgO平均值為 3.48;鈣堿系列黑云母富鎂, FeO*/MgO平均值為1.76。

    根據(jù)上述不同構(gòu)造背景巖漿巖中黑云母成分的顯著差異, 用FeO*-MgO, FeO*-Al2O3, Al2O3-MgO和Al2O3-MgO-FeO*分別作圖(式中 FeO*=[FeO+(Fe2O3×0.89981)]), 在這些圖解中可劃分出伸展(堿性)、俯沖有關(guān)(鈣堿性)和碰撞(過(guò)鋁質(zhì))等不同構(gòu)造背景的分布區(qū)。對(duì)這些圖解中堿性、鈣堿性和過(guò)鋁質(zhì)巖漿巖分布區(qū)進(jìn)行線性統(tǒng)計(jì)分析, 獲得了兩個(gè)線性判別函數(shù), 分別為:

    式中變量為對(duì)數(shù)比轉(zhuǎn)換(log-ratio transformation), 不是重量百分(有關(guān)對(duì)數(shù)比轉(zhuǎn)換方法原理可參考文獻(xiàn):周蒂, 1988, 1998)。用上述判別函數(shù)F1對(duì)F2作圖(圖8), 圖中劃分出A、C和P三部分, 分別代表伸展背景(堿性)、俯沖有關(guān)背景(鈣堿性)和碰撞背景(過(guò)鋁質(zhì)),相應(yīng)的判別準(zhǔn)確率分別為89%、91%和95%。該圖解與火成巖全巖主、微量元素構(gòu)造背景判別圖相結(jié)合,可以更準(zhǔn)確地判斷火成巖形成的構(gòu)造背景。

    圖7 以MORB中尖晶石為標(biāo)準(zhǔn)的不同構(gòu)造背景鉻鐵巖中的尖晶石主、微量元素圖解(據(jù)Gonzalez-Jimenez et al., 2012)Fig.7 Profiles of major and trace elements of spinel in eh chromitites from different tectonic settings normalized to composition of spinel of MORB

    5 角 閃 石

    不同構(gòu)造背景形成的角閃石對(duì)Nb、Ta的富集程度有明顯差異。Coltorfi et al. (2007)研究了來(lái)自巖石圈不同地幔區(qū)地幔包體的角閃石, 發(fā)現(xiàn)它們的微量元素,特別是高場(chǎng)強(qiáng)元素含量明顯不同: 產(chǎn)于俯沖帶上方角閃石(S 型)相對(duì)虧損 Nb(含量<10×10?6), Ti/Nb 和Zr/Nb比值高于球粒隕石, 而產(chǎn)于板內(nèi)的角閃石(I型)富集 Nb(含量在10×10?6~100×10?6), Ti/Nb, Zr/Nb低于球粒隕石。因此, 角閃石是了解俯沖帶上方流體、熔體交代特點(diǎn)的重要工具, 角閃石 Zr/Nb比值以及Zr/Nb-Ti/Nb, Zr/Nb-Ti/Zr作圖可以區(qū)分與不同構(gòu)造背景有關(guān)的、受多期交代的地幔區(qū)很有效(圖9)。

    圖 8 不同構(gòu)造背景火成巖中黑云母成分判別圖解(Abdel-Rahman, 1994)Fig.8 Discrimination diagram of two linear discriminant functions of biotite in different tectonic settings

    圖 9 區(qū)分板內(nèi)和俯沖帶上方角閃石的 Zr/Nb-Ti/Nb(a)和Zr/Nb-Ti/Zr(b)圖解(據(jù)Coltorfi et al., 2007)Fig.9 Zr/Nb vs. Ti/Nb (a) and Zr/Nb vs. Ti/Zr (b) diagrams for intraplate and suprasubduction amphiboles

    6 結(jié) 語(yǔ)

    火成巖中的礦物, 特別是穩(wěn)定副礦物微量元素地球化學(xué)的研究, 為探討火成巖(特別是蝕變的火成巖)構(gòu)造–巖漿背景提供了新的途徑, 可以預(yù)料, 隨著單顆粒礦物原位微量元素準(zhǔn)確測(cè)定技術(shù)的發(fā)展和資料的大量積累, 火成巖形成構(gòu)造背景的地球化學(xué)研究水平將迅速提高。

    趙振華. 2007. 關(guān)于巖石微量元素構(gòu)造環(huán)境判別圖解使用的有關(guān)問(wèn)題. 大地構(gòu)造與成礦學(xué), 31(1): 92–103.

    周蒂. 1988. 成分?jǐn)?shù)據(jù)的對(duì)數(shù)比統(tǒng)計(jì)方法簡(jiǎn)介. 地質(zhì)科技情報(bào), 7(2): 107–114.

    周蒂. 1998. 地質(zhì)成分?jǐn)?shù)據(jù)統(tǒng)計(jì)分析——困難與探索. 地球科學(xué), 23(2): 147–152.

    Abdel-Rahman A M. 1994. Nature of biotite from alkaline, calk-alkaline and peraluminous magmas. Journal of Petrology, 35: 525–541.

    Belousova E A, Griffin W L, O’Relly S Y and Fisher N J. 2002. Igneous zircon: Trace element composition as an indicator of source rock type. Contributions to Mineralogy and Petrology, 143: 602–623.

    Cameron M and Papike J J. 1981. Structual and chemical variations in pyroxene. American Mineralogist, 66: 1–50.

    Coltorti M, Bonadiman C, Faccini B, Grégoire M, O’Reilly S Y and Powell W. 2007. Amphiboles from suprasubduction and intraplate lithospheric mantle. Lithos, 99: 68–84.

    Debon F and Le Fort P. 1983. A chemical-mineralogical classification of common plutonic rocks and association. Transactions of the Royal Society Edinburgh: Earth Science, 73: 135–149.

    Farmer G L. 2003. Continental basaltic rocks // Holland H D and Turekian K K. Treties on Geochemistry, Vol. 3: The Crust. Oxford: Elsevier-Pergamon: 239–263

    Finch R J and Hanchar J M. 2003. Structure and chemistry of zircon and zircon-group minerals. Mineral Society America Reviews in Mineralogy and Geochemistry, 53: 1–25.

    Geisler T, Schaltergger U and Tomascher F. 2007. Re-equilibration of zircon in aqueous fluids and melts. Elements, 3: 43–50.

    Gonzalez-Jimenez J M G, Griffin W L, Locmelis M, O’Reilly S Y and Pearson N J. 2012. Contrasted minor-and trace-element compositions of spinel in chromitites of different tectonic settings. 22nd VM Goldshmidt Conference, Post.

    Grimes C B, John B E and Kelemen P B. 2007. Trace element chemistry of zircon from oceanic crust: A method for distinguishing detrital zircon provenance. Geology, 35(7): 643–646.

    Grimes C B, Wooden J L, Cheadle M J and John B E. 2015.“Fingerprinting” tectono-magmatic provenance using trace elements in igneous zircon. Contributions to Mineralogy and Petrology, 170: 46. DOI: 10.1007/ s00410-015-1199-3

    Hanchar J M and Westrenen W V. 2007. Rare earth element behavior in zircon-melt systems. Elements, 3: 37–42.

    Klein E M. 2003. Geochemistry of the oceanic crust // Holland H D and Turekian K K. Treties on geochemistry, Vol 3: The Crust. Amsterdam: Elsevier-Pogmon: 433–464.

    Le Bas M J. 1962. The role of aluminum in igneous clinopyroxenes with relation to their parentage. American Journal of Science, 260: 267–288.

    Leterrier J, Maury R C, Thonon P, Griard D and Marchal M. 1982. Clinopyroxene composition as a method of identification of the magmatic affinities of paleo-volcanic series. Earth and Planetary Science Letters, 59: 139–154.

    Li C S, Arndt N T, Tang Q Y and Ripley E M. 2015. Trace element indiscrimination diagrams. Lithos, 232: 76–83.

    Loucks R R. 1990. Discrimination of ophiolitic from nonophiolitic ultramafic-mafic allochthons in orogenic belts by the Al/Ti ratio in clinopyroxene. Geology, 18: 346–349.

    Müller D and Groves D I. 1997. Potassic Igneous Rocks and Associated Gold-copper Mineralization. Berlin: Springer: 11–40.

    Nisbet E G and Pearce J A. 1977. Clinopyroxene composition in mafic lavas from different tectonic setting. Contributions to Mineralogy and Petrology, 63: 149–160.

    Pearce J A. 2014. Immoble element fingerpringting of opgiolites. Elements, 10: 101–108.

    Rollinson H R. 1993. Using geochemical data: Evaluation, presentation, interpretation. New York: Longman Scientific Technical: 160–250.

    Rudnick R L and Gao S. 2003. Composition of the Continental Crust. Amsterdam: Elsevier: 3: 1–64.

    Tang G J, Wyman D A, Wang Q, Li J, Li Z X, Zhao Z H and Sun W D. 2012. Arthenosphere-lithosphere interaction triggered by a slab window during ridge subduction: Trace element and Sr-Nd-Hf-Os isotopeic evidence from late Carboniferous tholiites in the western Junggar area(NW China). Earth and Planetary Science Letters, 329-330: 84–96.

    Vermeesch P. 2006a. Tectonic discrimination of basalts with classification trees. Geochimica et Cosmochimica Acta, 70: 1839–1848.

    Vermeesch P. 2006b. Tectonic discrimination diagrams revisited. Geochemistry Geophysics Geosystems, 7: 1–55.

    Discrimination of Tectonic Settings Based on Trace Elements in Igneous Minerals

    ZHAO Zhenhua
    (Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou 510640, Guangdong, China)

    Tectonic settings of igneous rocks can be discriminated based on trace element geochemistry of their constitute minerals, particular the stable accessory minerals. In-situ analyses of trace element contents of minerals in igneous rocks have been compiled and applied to set up a group of discrimination diagrams for tectonic settings. Density distribution plots were created with two-dimensional kernel of trace elements in igneous zircon which present U/Yb-Nb/Yb, Sc/Yb-Nb/Yb and U/Yb-Sc/Yb diagrams. These diagrams are effective for distinguishing the mid-ocean ridge (MOR), magmatic arc (continental arc) and ocean island (and other plume-influenced) settings. Tectonic setting discrimination using clinopyroxene, spinel, biotite and amphibole are also introduced.

    tectonic setting discrimination; trace elements; zircon; clinopyroxene; spinel; biotite; amphibole

    P595

    A

    1001-1552(2016)05-0986-010

    10.16539/j.ddgzyckx.2016.05.008

    2016-03-11; 改回日期: 2016-03-28

    項(xiàng)目資助: NSFC-新疆聯(lián)合基金重點(diǎn)項(xiàng)目(U1203291)和國(guó)家自然科學(xué)基金(41273056)聯(lián)合資助。

    趙振華(1942–), 男, 研究員, 地球化學(xué)專業(yè)。Email: zhzhao@gig.ac.cn

    猜你喜歡
    單斜火成巖輝石
    新疆和巴西薔薇輝石玉的對(duì)比研究
    粉色薔薇輝石的寶石學(xué)及礦物學(xué)特征
    火成巖研磨性試驗(yàn)研究
    不同溫度、壓強(qiáng)、氧逸度條件下斜方輝石含水性的實(shí)驗(yàn)研究
    火成巖巖脈(墻)侵蝕對(duì)工作面的影響
    準(zhǔn)噶爾盆地西緣石炭系火成巖錄井綜合評(píng)價(jià)技術(shù)
    錄井工程(2017年3期)2018-01-22 08:40:24
    低溫單斜相BiPO4∶Tb3+納米晶的合成及其發(fā)光性能
    10位高速分級(jí)比較型單斜模數(shù)轉(zhuǎn)換器
    雙層厚硬火成巖破斷的力學(xué)分析
    南海新生代玄武巖中單斜輝石地球化學(xué)特征及其地質(zhì)意義
    免费女性裸体啪啪无遮挡网站| av片东京热男人的天堂| 国产一区二区三区av在线| 亚洲av欧美aⅴ国产| freevideosex欧美| 欧美激情国产日韩精品一区| 天天躁夜夜躁狠狠久久av| 亚洲美女搞黄在线观看| 99久久精品国产国产毛片| 熟女电影av网| 日韩电影二区| 啦啦啦视频在线资源免费观看| 免费在线观看黄色视频的| 国产成人精品无人区| 人妻一区二区av| 人妻系列 视频| 97在线人人人人妻| 久久青草综合色| 亚洲中文av在线| 丝袜美足系列| 国产麻豆69| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲av电影在线进入| 宅男免费午夜| 午夜视频国产福利| 欧美成人精品欧美一级黄| 亚洲精品日本国产第一区| 黑人欧美特级aaaaaa片| 久久人人爽人人爽人人片va| 亚洲综合色惰| 精品福利永久在线观看| 亚洲精品美女久久av网站| 久久狼人影院| 三上悠亚av全集在线观看| 亚洲人与动物交配视频| 久久女婷五月综合色啪小说| 免费观看性生交大片5| 天天躁夜夜躁狠狠久久av| 久久av网站| 久久久久久人人人人人| 久久人人爽av亚洲精品天堂| 男女边摸边吃奶| 成年女人在线观看亚洲视频| 精品一区二区三区视频在线| 丰满迷人的少妇在线观看| 1024视频免费在线观看| 99国产精品免费福利视频| 国产一区有黄有色的免费视频| 99久久人妻综合| 一个人免费看片子| 国产精品不卡视频一区二区| 熟女人妻精品中文字幕| 日韩一区二区三区影片| 建设人人有责人人尽责人人享有的| 欧美成人精品欧美一级黄| 午夜影院在线不卡| 99re6热这里在线精品视频| 久久毛片免费看一区二区三区| 午夜影院在线不卡| 少妇高潮的动态图| 亚洲精品aⅴ在线观看| 中文欧美无线码| 国产熟女午夜一区二区三区| 欧美少妇被猛烈插入视频| 日韩成人av中文字幕在线观看| 最近中文字幕高清免费大全6| 69精品国产乱码久久久| 午夜老司机福利剧场| 永久免费av网站大全| a 毛片基地| 日韩不卡一区二区三区视频在线| 久久狼人影院| 日本欧美国产在线视频| 最后的刺客免费高清国语| 久久久久久久久久久免费av| 午夜视频国产福利| 好男人视频免费观看在线| 丰满迷人的少妇在线观看| 啦啦啦在线观看免费高清www| 日本午夜av视频| 亚洲欧美一区二区三区黑人 | 丝袜美足系列| 91久久精品国产一区二区三区| 桃花免费在线播放| 视频中文字幕在线观看| 欧美日韩av久久| 三上悠亚av全集在线观看| 久久久精品免费免费高清| 国产精品.久久久| 22中文网久久字幕| 免费观看a级毛片全部| 国产色婷婷99| 在线看a的网站| 飞空精品影院首页| 女性被躁到高潮视频| 欧美日韩视频高清一区二区三区二| 国产欧美另类精品又又久久亚洲欧美| 国产亚洲精品久久久com| 青春草亚洲视频在线观看| 成人国产av品久久久| 亚洲欧美成人综合另类久久久| 国产精品一区二区在线不卡| 最近最新中文字幕大全免费视频 | 男男h啪啪无遮挡| 精品久久国产蜜桃| 成人手机av| 汤姆久久久久久久影院中文字幕| 男女啪啪激烈高潮av片| a级片在线免费高清观看视频| 久久久a久久爽久久v久久| 国产精品熟女久久久久浪| 国产成人a∨麻豆精品| 少妇的丰满在线观看| 亚洲欧美一区二区三区黑人 | 国产欧美日韩一区二区三区在线| 欧美日韩av久久| 女性被躁到高潮视频| 亚洲人与动物交配视频| 久久久久国产精品人妻一区二区| 亚洲国产精品一区三区| 亚洲欧洲日产国产| 波多野结衣一区麻豆| 2022亚洲国产成人精品| 亚洲成人av在线免费| 日本黄色日本黄色录像| 国产乱人偷精品视频| 18禁观看日本| 下体分泌物呈黄色| 在线观看美女被高潮喷水网站| 两个人看的免费小视频| 亚洲情色 制服丝袜| 91久久精品国产一区二区三区| 男女高潮啪啪啪动态图| 国产综合精华液| 亚洲av电影在线观看一区二区三区| www.熟女人妻精品国产 | 国产精品女同一区二区软件| 精品久久久久久电影网| 少妇被粗大猛烈的视频| 汤姆久久久久久久影院中文字幕| 18在线观看网站| 午夜福利视频在线观看免费| 美女视频免费永久观看网站| 丰满迷人的少妇在线观看| 黑人高潮一二区| 国产成人一区二区在线| 中文字幕制服av| 国产乱来视频区| 97精品久久久久久久久久精品| 日韩免费高清中文字幕av| 国产老妇伦熟女老妇高清| 青青草视频在线视频观看| 在线观看国产h片| 桃花免费在线播放| www.色视频.com| 亚洲国产av新网站| 国产精品一国产av| 搡老乐熟女国产| 欧美国产精品va在线观看不卡| 亚洲国产精品国产精品| 免费看不卡的av| 日韩电影二区| 在线精品无人区一区二区三| 久久久久久久久久久免费av| 欧美xxⅹ黑人| 国产成人精品久久久久久| 久久免费观看电影| 一级片'在线观看视频| 国产精品三级大全| 9热在线视频观看99| 国产乱人偷精品视频| 大香蕉久久网| 巨乳人妻的诱惑在线观看| 午夜免费男女啪啪视频观看| 边亲边吃奶的免费视频| 日韩不卡一区二区三区视频在线| 天堂8中文在线网| 五月天丁香电影| 亚洲av在线观看美女高潮| 高清在线视频一区二区三区| 在线天堂中文资源库| 亚洲精品乱码久久久久久按摩| 美女国产视频在线观看| videossex国产| 在线观看国产h片| 美国免费a级毛片| 18禁动态无遮挡网站| 色5月婷婷丁香| 国产精品国产三级专区第一集| 少妇的逼好多水| 久久精品久久久久久久性| 国产男女超爽视频在线观看| 久久热在线av| 欧美精品高潮呻吟av久久| 国产极品天堂在线| av电影中文网址| 国产成人午夜福利电影在线观看| 亚洲伊人久久精品综合| 99热全是精品| 制服丝袜香蕉在线| 亚洲av电影在线观看一区二区三区| 韩国精品一区二区三区 | 在线观看国产h片| 亚洲精品一区蜜桃| 亚洲精品久久久久久婷婷小说| 欧美少妇被猛烈插入视频| 国产老妇伦熟女老妇高清| 国产成人精品一,二区| www.av在线官网国产| 日本免费在线观看一区| 精品视频人人做人人爽| 欧美丝袜亚洲另类| 日本91视频免费播放| 亚洲精品一二三| 插逼视频在线观看| av在线观看视频网站免费| 日韩精品有码人妻一区| 成人国语在线视频| 人人妻人人添人人爽欧美一区卜| 国产片特级美女逼逼视频| 老熟女久久久| 久久精品熟女亚洲av麻豆精品| 亚洲欧美成人精品一区二区| 久久女婷五月综合色啪小说| 在线看a的网站| 大香蕉久久网| 亚洲中文av在线| 日韩熟女老妇一区二区性免费视频| 国产一区二区激情短视频 | 少妇的逼水好多| 国产在线视频一区二区| 国产男女超爽视频在线观看| 精品熟女少妇av免费看| 看非洲黑人一级黄片| 久久久国产一区二区| 视频中文字幕在线观看| 欧美日韩视频高清一区二区三区二| 视频中文字幕在线观看| 国语对白做爰xxxⅹ性视频网站| 人妻少妇偷人精品九色| 久久精品国产a三级三级三级| 99久久人妻综合| 久久精品人人爽人人爽视色| 欧美精品国产亚洲| 久久精品国产综合久久久 | 欧美性感艳星| 两个人看的免费小视频| 午夜福利影视在线免费观看| 国产乱来视频区| 国产乱人偷精品视频| 在线精品无人区一区二区三| 人人妻人人澡人人看| 考比视频在线观看| 日产精品乱码卡一卡2卡三| 日韩免费高清中文字幕av| 国产成人精品一,二区| 亚洲av欧美aⅴ国产| 搡老乐熟女国产| 国产69精品久久久久777片| 日韩,欧美,国产一区二区三区| 视频区图区小说| 一二三四中文在线观看免费高清| 亚洲精品国产av蜜桃| 久久精品人人爽人人爽视色| 日韩中文字幕视频在线看片| a级毛片黄视频| 热99久久久久精品小说推荐| 国产欧美日韩综合在线一区二区| 亚洲av免费高清在线观看| 午夜福利乱码中文字幕| 大香蕉久久成人网| 久久午夜综合久久蜜桃| 在线亚洲精品国产二区图片欧美| av.在线天堂| 女人精品久久久久毛片| 国产成人免费无遮挡视频| 欧美 日韩 精品 国产| 在线观看免费视频网站a站| 色网站视频免费| 国产成人精品婷婷| 9热在线视频观看99| 蜜桃在线观看..| 交换朋友夫妻互换小说| 亚洲国产av影院在线观看| 久久综合国产亚洲精品| 一级黄片播放器| 国产日韩一区二区三区精品不卡| 国产欧美另类精品又又久久亚洲欧美| 制服人妻中文乱码| 国产成人av激情在线播放| 在线 av 中文字幕| 免费黄网站久久成人精品| 亚洲精华国产精华液的使用体验| 90打野战视频偷拍视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品国产av成人精品| 午夜福利乱码中文字幕| 国产欧美日韩综合在线一区二区| 精品一区二区三区视频在线| 性高湖久久久久久久久免费观看| 欧美精品国产亚洲| 亚洲国产欧美日韩在线播放| 最近中文字幕2019免费版| 成人国产av品久久久| 国产精品一区www在线观看| av在线老鸭窝| 国产一区二区三区综合在线观看 | 美女国产视频在线观看| 国产黄色视频一区二区在线观看| 伦理电影大哥的女人| 色吧在线观看| 视频中文字幕在线观看| 精品亚洲成国产av| 纵有疾风起免费观看全集完整版| 中文字幕av电影在线播放| 五月天丁香电影| 国产av码专区亚洲av| 91精品伊人久久大香线蕉| 亚洲久久久国产精品| 夫妻午夜视频| 久久精品aⅴ一区二区三区四区 | 久久久久久伊人网av| 精品99又大又爽又粗少妇毛片| 欧美亚洲 丝袜 人妻 在线| 国产片特级美女逼逼视频| 少妇高潮的动态图| 免费观看无遮挡的男女| 国产午夜精品一二区理论片| 久久这里只有精品19| 国产乱来视频区| 国产福利在线免费观看视频| 国产亚洲精品久久久com| videosex国产| 两个人看的免费小视频| 午夜福利视频在线观看免费| 五月伊人婷婷丁香| 秋霞在线观看毛片| 99视频精品全部免费 在线| 国产视频首页在线观看| 宅男免费午夜| 亚洲情色 制服丝袜| 国产精品国产av在线观看| 国产成人精品婷婷| 日韩精品免费视频一区二区三区 | 熟妇人妻不卡中文字幕| 国产一区二区在线观看日韩| 国产欧美日韩一区二区三区在线| 国产精品一国产av| 免费看av在线观看网站| av免费观看日本| 国产无遮挡羞羞视频在线观看| 91午夜精品亚洲一区二区三区| 中文乱码字字幕精品一区二区三区| 香蕉丝袜av| 精品少妇内射三级| 国产高清不卡午夜福利| 亚洲精品国产av成人精品| 国产男女超爽视频在线观看| 性色avwww在线观看| 免费av不卡在线播放| 中文字幕av电影在线播放| 捣出白浆h1v1| 91精品国产国语对白视频| 九草在线视频观看| 美女福利国产在线| 久久久国产精品麻豆| 91精品国产国语对白视频| 九草在线视频观看| 精品少妇内射三级| 九九在线视频观看精品| 亚洲精品国产av成人精品| 美女脱内裤让男人舔精品视频| 99九九在线精品视频| 欧美最新免费一区二区三区| 精品国产国语对白av| 欧美最新免费一区二区三区| av电影中文网址| 九九在线视频观看精品| 亚洲av免费高清在线观看| 成年美女黄网站色视频大全免费| 国内精品宾馆在线| 亚洲av日韩在线播放| 老司机影院成人| 高清黄色对白视频在线免费看| 色婷婷久久久亚洲欧美| 精品人妻偷拍中文字幕| 国产乱来视频区| 大香蕉久久成人网| 欧美成人精品欧美一级黄| 哪个播放器可以免费观看大片| 国产一级毛片在线| 日韩不卡一区二区三区视频在线| av在线app专区| 观看av在线不卡| 狂野欧美激情性bbbbbb| 九九在线视频观看精品| 观看美女的网站| 丁香六月天网| 久久99热6这里只有精品| 欧美激情 高清一区二区三区| 日本vs欧美在线观看视频| 欧美成人午夜免费资源| 99久久中文字幕三级久久日本| 高清在线视频一区二区三区| 亚洲精品色激情综合| 国产精品麻豆人妻色哟哟久久| 五月玫瑰六月丁香| 成年女人在线观看亚洲视频| 成人影院久久| 国产毛片在线视频| 久久精品久久精品一区二区三区| 波野结衣二区三区在线| 一区二区三区乱码不卡18| 免费观看av网站的网址| 午夜激情av网站| 亚洲av国产av综合av卡| 国产精品一区www在线观看| 免费女性裸体啪啪无遮挡网站| 国产精品嫩草影院av在线观看| 亚洲欧美一区二区三区黑人 | 新久久久久国产一级毛片| 婷婷色综合大香蕉| 久久精品国产鲁丝片午夜精品| 26uuu在线亚洲综合色| 免费黄网站久久成人精品| 久久久久久伊人网av| 国产高清三级在线| 2018国产大陆天天弄谢| 国产精品不卡视频一区二区| 国精品久久久久久国模美| 国产成人精品在线电影| 一本大道久久a久久精品| 国产成人精品无人区| 久久精品夜色国产| 男人操女人黄网站| 免费黄网站久久成人精品| av免费在线看不卡| 成人综合一区亚洲| 国产片内射在线| 久久久久久久久久成人| 国产成人免费无遮挡视频| 90打野战视频偷拍视频| 国产精品熟女久久久久浪| 久久久久久久久久久久大奶| 久久热在线av| 国产精品不卡视频一区二区| 看免费av毛片| 90打野战视频偷拍视频| 午夜福利影视在线免费观看| 22中文网久久字幕| a级毛色黄片| 五月伊人婷婷丁香| 亚洲精品一区蜜桃| 2018国产大陆天天弄谢| 欧美激情 高清一区二区三区| 校园人妻丝袜中文字幕| 精品国产一区二区三区久久久樱花| 一本久久精品| 久久免费观看电影| 美女视频免费永久观看网站| 精品一区二区免费观看| 国产淫语在线视频| 久久ye,这里只有精品| 久久亚洲国产成人精品v| 人人妻人人添人人爽欧美一区卜| 夫妻性生交免费视频一级片| 国产精品熟女久久久久浪| 欧美精品人与动牲交sv欧美| 亚洲精品中文字幕在线视频| 搡女人真爽免费视频火全软件| 涩涩av久久男人的天堂| 一边亲一边摸免费视频| 国产av一区二区精品久久| 激情五月婷婷亚洲| 国产探花极品一区二区| 人妻人人澡人人爽人人| 男女啪啪激烈高潮av片| 欧美成人午夜精品| 波多野结衣一区麻豆| 精品亚洲乱码少妇综合久久| 欧美精品国产亚洲| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 美女xxoo啪啪120秒动态图| 欧美激情国产日韩精品一区| 国精品久久久久久国模美| 黑丝袜美女国产一区| 国产精品久久久久久精品电影小说| 国产永久视频网站| 国产精品一二三区在线看| 成人国产麻豆网| 久久久a久久爽久久v久久| 久久久精品94久久精品| 免费看光身美女| 人妻人人澡人人爽人人| 午夜日本视频在线| 18禁动态无遮挡网站| 一本色道久久久久久精品综合| 80岁老熟妇乱子伦牲交| 免费观看无遮挡的男女| 大话2 男鬼变身卡| 99热这里只有是精品在线观看| 视频区图区小说| 欧美97在线视频| 成人影院久久| 波野结衣二区三区在线| 国产精品99久久99久久久不卡 | 亚洲国产精品一区二区三区在线| 秋霞在线观看毛片| 久久这里只有精品19| 免费人成在线观看视频色| 黑人欧美特级aaaaaa片| 欧美最新免费一区二区三区| 日韩 亚洲 欧美在线| 午夜av观看不卡| 亚洲精品一区蜜桃| 久久久精品94久久精品| videosex国产| 欧美精品亚洲一区二区| 国产欧美日韩综合在线一区二区| 男人添女人高潮全过程视频| xxxhd国产人妻xxx| 欧美日韩一区二区视频在线观看视频在线| 欧美成人午夜精品| 国产成人91sexporn| 成人午夜精彩视频在线观看| 国产成人免费无遮挡视频| 大话2 男鬼变身卡| 亚洲精品456在线播放app| 国产麻豆69| 久热这里只有精品99| xxxhd国产人妻xxx| 日本黄大片高清| 成人二区视频| 精品午夜福利在线看| 伦精品一区二区三区| 丰满少妇做爰视频| 伊人亚洲综合成人网| 纯流量卡能插随身wifi吗| 国产成人精品婷婷| 国产欧美亚洲国产| 曰老女人黄片| 久久精品久久久久久噜噜老黄| 夜夜爽夜夜爽视频| 人人澡人人妻人| 狠狠婷婷综合久久久久久88av| 亚洲欧美日韩另类电影网站| 精品一区二区免费观看| 日韩在线高清观看一区二区三区| 亚洲国产最新在线播放| 欧美亚洲 丝袜 人妻 在线| 精品人妻一区二区三区麻豆| 中国国产av一级| 美女脱内裤让男人舔精品视频| 一边摸一边做爽爽视频免费| 国产成人免费观看mmmm| 午夜免费男女啪啪视频观看| av国产久精品久网站免费入址| a级毛片在线看网站| 三上悠亚av全集在线观看| 日本猛色少妇xxxxx猛交久久| 国产欧美亚洲国产| 国产 精品1| 亚洲欧美成人综合另类久久久| 成年女人在线观看亚洲视频| 卡戴珊不雅视频在线播放| a级毛片在线看网站| av线在线观看网站| 激情五月婷婷亚洲| 精品国产一区二区三区久久久樱花| a级毛色黄片| 十分钟在线观看高清视频www| 国产成人午夜福利电影在线观看| 欧美日韩国产mv在线观看视频| 天天躁夜夜躁狠狠久久av| 久久久久久久国产电影| 亚洲av.av天堂| 男人添女人高潮全过程视频| 中文字幕最新亚洲高清| 十八禁网站网址无遮挡| 美女内射精品一级片tv| 日本黄大片高清| 亚洲五月色婷婷综合| 亚洲熟女精品中文字幕| 亚洲精品一区蜜桃| 啦啦啦在线观看免费高清www| 热99久久久久精品小说推荐| 亚洲精品久久午夜乱码| 成人免费观看视频高清| 2021少妇久久久久久久久久久| 日韩一本色道免费dvd| 欧美3d第一页| 最新中文字幕久久久久| 亚洲国产毛片av蜜桃av| 丰满乱子伦码专区| 日本黄色日本黄色录像| 性色avwww在线观看| 伊人亚洲综合成人网| 国产69精品久久久久777片| 人人妻人人澡人人看| 国产成人精品一,二区| 久久ye,这里只有精品| 国产精品久久久久成人av| 国产福利在线免费观看视频| 国产日韩欧美视频二区| 成人国产麻豆网| 成人综合一区亚洲| 黄片播放在线免费| 色94色欧美一区二区| av女优亚洲男人天堂| 成人黄色视频免费在线看| 91成人精品电影| 狂野欧美激情性xxxx在线观看| 成人免费观看视频高清| 精品国产国语对白av| 一级,二级,三级黄色视频|