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

    土耳其西部新生代斑巖成礦巖漿氧逸度-揮發(fā)份元素組成演化特征*
    ——來(lái)自鋯石和磷灰石組成約束

    2021-09-24 06:11:08趙俊興蘇本勛UYSALIbrahimAYDINFaruk肖燕SENCuneyt回凱旋秦克章
    巖石學(xué)報(bào) 2021年8期
    關(guān)鍵詞:逸度磷灰石斑巖

    趙俊興 蘇本勛 UYSAL Ibrahim AYDIN Faruk 肖燕 SEN Cuneyt回凱旋, 3 秦克章, 3

    1. 中國(guó)科學(xué)院礦產(chǎn)資源研究重點(diǎn)實(shí)驗(yàn)室,中國(guó)科學(xué)院地質(zhì)與地球物理研究所,北京 100029

    2. 中國(guó)科學(xué)院地球科學(xué)研究院,北京 100029

    3. 中國(guó)科學(xué)院大學(xué)地球與行星科學(xué)學(xué)院,北京 100049

    4. Department of Geological Engineering, Karadeniz Technical University, Trabzon 61080

    5. 巖石圈演化國(guó)家重點(diǎn)實(shí)驗(yàn)室,中國(guó)科學(xué)院地質(zhì)與地球物理研究所,北京 100029

    圖1 特提斯土耳其段主要斑巖型礦床分布圖(據(jù)修改)其中土耳其西部主要為四個(gè)成礦階段:1-俯沖階段;2-碰撞階段(Tavanl帶);3-晚碰撞-后碰撞階段(Biga半島成礦帶);4-后碰撞伸展階段(Afyon-Konya成礦帶);成礦帶劃分來(lái)自Kucu et al. (2019). 本次研究選擇2-4區(qū)域內(nèi)典型礦床開展研究Fig.1 Distribution map of main porphyry deposits of Tethys Turkish section (modified after 2009)The western Turkish is mainly composed of four ore-forming stages: 1-subduction stage; 2-collision stage (Tavanl belt); 3-from late-collision to post-collision stage (Biga Peninsula metallogenic belt); 4-post-collision extension stage (Afyon-Konya metallogenic belt); the division of metallogenic belt refers to Kucu et al. (2019). The typical deposits of 2~4 region were selected to research in this study

    1 土耳其西部斑巖礦床構(gòu)造背景

    2 土耳其西部新生代典型斑巖成礦系統(tǒng)基本特征

    2.1 Tavanl帶始新世斑巖成礦系統(tǒng)

    2.2 Biga半島始新世至漸新世斑巖成礦系統(tǒng)

    2.3 Afyon-Konya中新世-上新世斑巖金-鉬-銅礦及淺成低溫成礦系統(tǒng)

    3 研究對(duì)象和測(cè)試方法

    鋯石和磷灰石在分選后制靶,并拍攝透反射光照片和選擇無(wú)裂縫和礦物-流體包裹體顆粒區(qū)域,以便后續(xù)開展礦物主微量元素測(cè)試工作。鋯石和磷灰石的陰極發(fā)光圖像主要在中國(guó)科學(xué)院地質(zhì)與地球物理研究所使用配置MonoCL4的Nova NanoSEM 450場(chǎng)發(fā)射掃描電子顯微鏡(FSEM)完成。鋯石微量元素含量中國(guó)科學(xué)院地質(zhì)與地球物理研究所巖石圈演化重點(diǎn)實(shí)驗(yàn)室的Agilent 7500a四級(jí)桿電感耦合等離子體質(zhì)譜儀上完成(與193nm ArF準(zhǔn)分子激光剝蝕儀聯(lián)機(jī))。每個(gè)時(shí)間分辨分析數(shù)據(jù)包括大約20~30s的空白信號(hào)和50s的樣品信號(hào)。對(duì)分析數(shù)據(jù)的離線處理(包括對(duì)樣品和空白信號(hào)的選擇、儀器靈敏度漂移校正、元素含量及U-Th-Pb同位素比值和年齡計(jì)算)采用軟件ICPMSDataCal(Liuetal., 2008)完成。詳細(xì)的儀器操作條件和數(shù)據(jù)處理方法同Liuetal.(2008, 2010)。鋯石微量元素選點(diǎn)主要為先前鋯石U-Pb定年測(cè)試點(diǎn),激光束斑選擇為44~60μm,略大于之前的年齡測(cè)定束斑。鋯石微量元素含量利用NIST610、USGS參考玻璃BCR-2G與BIR-1G作為多外標(biāo)、Zr作內(nèi)標(biāo)的方法進(jìn)行定量計(jì)算(Liuetal., 2010)。鋯石微量元素結(jié)果見于表2。

    表2 土耳其西部新生代斑巖礦床中巖漿巖鋯石LA-ICPMS微量元素?cái)?shù)據(jù)(×10-6)Table 2 The analytical data of LA-ICPMS zircon trace elements from magmatic rocks in the Cenozoic porphyry ore deposits, western Turkey (×10-6)

    續(xù)表2Continued Table 2

    磷灰石主量元素在中國(guó)科學(xué)院地質(zhì)與地球物理研究所(IGGCAS)使用JEOL-JXA8100電子探針完成。本次研究工作條件為加速電壓15kV,束流10nA,探針束斑大小為3~5μm,Si、Fe、Mn、Mg、Na、S和Cl的計(jì)數(shù)時(shí)間為20s,F(xiàn)的計(jì)數(shù)時(shí)間為40s,Ca和P的計(jì)數(shù)時(shí)間為10s。選點(diǎn)過(guò)程中我們只將感興趣的磷灰石顆粒在電子束下進(jìn)行移動(dòng),以防止F、Cl、Ca和P成分的遷移。Stormeretal.(1993)指出這些元素在EPMA分析中的Kα-X射線計(jì)數(shù)率與電子束的照射時(shí)間密切相關(guān)。所有數(shù)據(jù)均采用ZAF進(jìn)行校正,所有結(jié)果見表3。

    表3 土耳其西部新生代斑巖礦床中巖漿巖磷灰石EPMA主量元素?cái)?shù)據(jù)(wt%)Table 3 The analytical data of EPMA apatite major elements from magmatic rocks in the Cenozoic porphyry ore deposits, western Turkey (wt%)

    續(xù)表3Continued Table 3

    續(xù)表3Continued Table 3

    4 測(cè)試結(jié)果

    4.1 鋯石微量元素組成

    本次研究共測(cè)定333顆鋯石微量元素?cái)?shù)據(jù),結(jié)合全巖數(shù)據(jù)(趙俊興等,未發(fā)表數(shù)據(jù))和Lietal.(2019)的Geo-fO2來(lái)計(jì)算鋯石鈦飽和溫度和氧逸度(表2)。結(jié)果顯示土耳其西部新生代斑巖礦床的成礦巖漿具有如下的特征:(1)鋯石微量元素表明土耳其西部新生代斑巖礦床成礦巖漿大部分落于高水含量-高氧逸度區(qū)間(圖2);(2)巖漿巖鋯石具有相似的稀土元素配分曲線(圖3),且暗色包體鋯石較同期巖體有較高的稀土總量(圖3a, d);(3)絕大部分鋯石的EuN/EuN*位于0.4~0.7之間,斑巖Mo和斑巖Mo-Cu礦相關(guān)鋯石在結(jié)晶時(shí)熔體受斜長(zhǎng)石結(jié)晶影響更大,具有更低EuN/EuN*比值(圖4a);(4)利用Ferry and Watson (2007)的鋯石Ti飽和溫度計(jì)(假設(shè)TiO2活度為0.7)估算土耳其西部斑巖成礦及其相關(guān)巖體的形成溫度在650~900℃之間,鋯石Ti溫度顯示巖體結(jié)晶溫度表現(xiàn)為:斑巖Au礦>斑巖Cu-Mo=斑巖Mo-Cu礦>斑巖Mo礦,呈現(xiàn)遞減趨勢(shì) (圖4b);(5)鋯石微量元素顯示(Yb/Gd和Hf/Y-Th/U,圖4c, d)鋯石在巖漿房結(jié)晶過(guò)程中,其微量元素組成大部分受角閃石±榍石±磷灰石分離結(jié)晶控制,該分離結(jié)晶控制趨勢(shì)與全巖微量元素指示演化受角閃石結(jié)晶相控制一致(圖5);(6)斑巖Au和斑巖Cu-Mo相關(guān)巖體鋯石形成時(shí)氧逸度高于斑巖Mo和斑巖Mo-Cu礦巖體,結(jié)合各巖體形成年齡,鋯石氧逸度Ce4+/Ce3+比值呈現(xiàn)出先降低(斑巖Cu-Mo到斑巖Mo)后升高(斑巖Mo-Cu到斑巖Au)的趨勢(shì)(圖4d)。

    圖2 鋯石(Ce/Nd)/Y-10000×(Eu/Eu*)/Y圖解其中鋯石微量元素組成圈定高氧逸度-高含水量成礦巖體(△FMQ>2和H2O<12%)與低氧逸度-低含水量(△FMQ<2和H2O<5%)非成礦巖體范圍來(lái)自Lu et al. (2016)Fig.2 Zircon (Ce/Nd)/Y vs. 10000×(Eu/Eu*)/Y diagramThe zircon trace element ranges to define high oxygen fugacity-high water content fertile magma (△FMQ>2 and H2O<12%) and low oxygen fugacity-low water content (△FMQ<2 and H2O<5%) infertile suite refer to Lu et al. (2016)

    圖3 土耳其西部新生代斑巖礦床巖漿巖鋯石球粒隕石標(biāo)準(zhǔn)化稀土元素配分圖解(標(biāo)準(zhǔn)化值據(jù)Sun and McDonough, 1989)Fig.3 Chondrite-normalization REE patterns of western Turkish Neozoic porphyry deposit magmatic rock zircon (normalization values after Sun and McDonough, 1989)

    圖4 鋯石Hf-EuN/EuN*(a)、Hf-鋯石Ti飽和溫度(b)、Th/U-Yb/Gd(c)和Th/U-Hf/Y(d)之間演化關(guān)系圖a中氧化條件演化趨勢(shì)來(lái)自Dilles et al. (2015),圖c和d中分離結(jié)晶趨勢(shì)來(lái)自Lee et al. (2021). 圖例見圖2Fig.4 The evolution relationship of zircon Hf vs. EuN/EuN* (a), Hf vs. zircon Ti saturation temperature (b), Th/U vs. Yb/Gd (c) and Th/U vs. Hf/Y (d)The evolution trend of oxidizing condition in Fig.4a refers to Dilles et al. (2015), the fractional crystallization trend in Fig.4c, d refers to Lee et al. (2021). Legends refer to Fig.2

    圖5 土耳其西部斑巖礦床中鋯石年齡和鋯石Th/U、Ti溫度、EuN/EuN*和Ce4+/Ce3+比值之間演化關(guān)系Fig.5 The evolution relationship of zircon age and zircon Th/U, Ti-in-zircon temperature, EuN/EuN* and Ce4+/Ce3+ ratio in the porphyry ore deposits, western Turkey

    4.2 磷灰石主量元素組成

    圖6 斑巖礦床中磷灰石中F含量-F/Cl比值(a)和Cl含量-F/Cl比值(b)之間的關(guān)系Fig.6 The relationship of apatite F content vs. F/Cl ratios (a) and Cl content vs. F/Cl ratios (b) in porphyry deposits

    5 討論

    5.1 土耳其西部新生代斑巖系統(tǒng)巖漿氧逸度演化趨勢(shì)

    通過(guò)斑巖成礦帶的構(gòu)造演化歷史和巖漿巖微量元素-同位素研究,能夠獲得俯沖-碰撞過(guò)程中地殼屬性(Sillitoe, 2018)、地殼厚度演化(如, (La/Yb)N比值, Profetaetal., 2015; Zhuetal., 2017; Y/Sr比值, Dengetal., 2018; 巖漿金屬元素含量和俯沖帶地殼厚度關(guān)系, Chiaradia, 2014)、斑巖成礦帶的巖石圈結(jié)構(gòu)(如, 同位素填圖, Houetal., 2015b; 地?;蛳碌貧ぐw研究, Wangetal., 2016, Xuetal., 2017)、俯沖板塊邊界信息(如, 配合板塊匯聚速度和地殼縮短估算, Wangetal., 2015)等宏觀要素。

    本次研究得到土耳其西部新生代斑巖系統(tǒng)的巖漿氧逸度自晚古新世-早始新世碰撞階段到后碰撞環(huán)境經(jīng)歷先降低后升高的過(guò)程(圖5d,Ce4+/Ce3+比值的變化關(guān)系)。該特征與西藏岡底斯碰撞環(huán)境斑巖礦床(Wangetal., 2014a, b; Yangetal., 2016)和安第斯型俯沖環(huán)境斑巖礦床(Caoetal., 2016; Lietal., 2017; Angereretal., 2018)的氧逸度演化趨勢(shì)不同。西藏岡底斯碰撞環(huán)境和安第斯俯沖型成礦巖漿演化趨勢(shì)表明為成礦巖漿含水量和氧逸度均顯著升高,地殼加厚過(guò)程可能是控制其氧逸度升高的重要因素(Tangetal., 2020)。因此,估算土耳其西部斑巖礦床成礦巖體及其同期巖漿巖形成時(shí)地殼厚度,并結(jié)合鋯石氧逸度評(píng)價(jià)地殼厚度演化過(guò)程對(duì)巖漿氧逸度的控制過(guò)程,是解釋土耳其西部新生代斑巖系統(tǒng)巖漿氧逸度演化趨勢(shì)的重要研究?jī)?nèi)容。

    本次研究使用Huetal.(2017)在碰撞造山帶巖漿巖研究中提出的地殼厚度估算公式。Huetal.(2017)利用全巖Sr/Y和(La/Yb)N比值以及Crust 1.0模型獲得碰撞帶地殼厚度(H)計(jì)算公式估算地殼厚度:H=0.67×Sr/Y+28.21和H=27.78×ln[0.34×(La/Yb)N],并指出與前人在俯沖島弧帶巖漿巖中獲得的地殼厚度估算模型相比(Chiaradia, 2014; Profetaetal., 2015),在相同地殼厚度下碰撞環(huán)境的中酸性巖漿巖(SiO2含量為55%~72%)具有較低的Sr/Y和較高的 (La/Yb)N比值(Huetal., 2017)。我們結(jié)合前人在土耳其西部地區(qū)發(fā)表的新生代巖漿巖數(shù)據(jù)(選擇SiO2含量在55%~72%),使用該模型對(duì)土耳其西部50Ma以來(lái)的斑巖礦床及同期巖漿巖形成時(shí)地殼厚度進(jìn)行估算(圖7)。估算結(jié)果表明三條成礦帶內(nèi)成礦與非成礦巖體盡管估算值略有不同,但隨時(shí)間演化趨勢(shì)基本一致,即Tavanl帶和Biga半島成礦帶在始新世到漸新世巖漿巖形成時(shí)地殼厚度估算值略微升高,Afyon-Konya成礦帶早中新世地殼依然有與Biga半島成礦帶相似的估算厚度。但到中晚中新世Afyon-Konya成礦帶的地殼厚度明顯增大。而針對(duì)本次研究所選擇的樣品,全巖的(La/Yb)N和估算地殼厚度(圖8a, b)顯示從始新世到早中新世有著地殼逐漸加厚的趨勢(shì)。而在漸新世-中新世的Tepeoba斑巖Mo礦、Pnarba斑巖Mo-Cu礦和斑巖Au礦中,從成礦前到成礦期巖漿演化過(guò)程明顯經(jīng)歷地殼減薄事件(圖8a, b,隨時(shí)間明顯降低的(La/Yb)N)。在該過(guò)程中Pnarba斑巖Mo-Cu礦和斑巖Au礦的巖漿氧逸度明顯上升(圖8c),說(shuō)明與地殼減薄相關(guān)的深部過(guò)程可能是控制土耳其西部中新世高氧逸度成礦巖漿形成的重要機(jī)制。同時(shí),全巖微量元素(圖9a,b)和鋯石微量元素(圖4d)指示,土耳其西部斑巖礦床成礦巖漿地球化學(xué)特征可能主要受到以角閃石為主的分離結(jié)晶作用控制。因此,在該地殼減薄過(guò)程中可能發(fā)生基性巖漿注入過(guò)程,使得成礦期巖漿氧逸度升高。

    圖7 土耳其西部斑巖礦床巖體及同時(shí)期巖漿巖地殼厚度估算與形成時(shí)代演化圖圖中地殼厚度估算公式來(lái)自Hu et al. (2017),并分別利用全巖(La/Yb)N(a)和Sr/Y比值(b)進(jìn)行估算,其中Tavanl帶巖體數(shù)據(jù)來(lái)自Altunkaynak (2007),Kucu et al. (2019)和Güraslan and Altunkaynak (2019);Biga半島成礦帶數(shù)據(jù)來(lái)自Kucu et al. (2019);Afyon-Konya成礦帶數(shù)據(jù)來(lái)自 et al. (2015), Deliba et al. (2017)和Kucu et al.(2019)Fig.7 The evolution diagrams of western Turkish porphyry deposits rock and contemporary magmatic rock crustal thickness estimation and formation ageThe computational formula of crustal thickness refers to Hu et al. (2017), whole rock (La/Yb)N (a) and Sr/Y (b) ratio was used to calculate respectively, the Tavanl belt rock values refer to Altunkaynak (2007), Kucu et al. (2019) and Güraslan and Altunkaynak (2019); Biga Peninsula metallogenic belt values refer to Kucu et al. (2019); Afyon-Konya metallogenic belt values refer to et al. (2015), Deliba et al. (2017) and Kucu et al. (2019)

    圖8 斑巖礦床中各類巖體鋯石年齡-(La/Yb)N(a)、鋯石年齡-地殼估算厚度(b)和巖漿氧逸度-地殼估算厚度(c)圖解△FMQ通過(guò)鋯石計(jì)算獲得,以平均值和1StDev表示Fig.8 Magmatic zircon age vs. whole-rock (La/Yb)N (a), zircon age vs. crustal estimation thickness (b), magmatic oxygen fugacity (vs. crustal estimation thickness diagrams△FMQ was obtained by zircon calculation, and showed in terms of average and 1StDev

    圖9 本次研究的斑巖礦床中各類巖體SiO2-Dy/Yb(a)和SiO2-(La/Yb)N(b)圖解其中角閃石和石榴石演化趨勢(shì)來(lái)自Davidson et al. (2007),全巖數(shù)據(jù)來(lái)自趙俊興等(未發(fā)表數(shù)據(jù))Fig.9 Different types of magmatic rock SiO2 vs. Dy/Yb (a) and SiO2 vs. (La/Yb)N (b) diagrams in this studyThe evolution trends of amphibole and garnet refer to Davidson et al. (2007). The whole-rock geochemical compositions are cited from Zhao et al. (unpublished data)

    5.2 巖漿中揮發(fā)份F-Cl-S演化特征——來(lái)自磷灰石組成的約束

    斑巖型礦床的金屬元素主要以Cl、S絡(luò)合物的形式搬運(yùn),以金屬硫化物的形式沉淀。如Cu與Au在熔體中低氧逸度條件下以氯絡(luò)合物形式遷移,金在高氧逸度條件下以硫氫絡(luò)合物形式遷移,巖漿系統(tǒng)高氯含量有助于提高銅和金的溶解度(Zajacz and Halter, 2009)。而從斑巖銅礦到斑巖鉬礦,巖體酸性程度常與磷灰石F/Cl比值成正比(Mathez and Webster, 2005; 陳雷等, 2017; Zhaoetal., 2020)。因此,研究斑巖礦床成礦物質(zhì)富集-演化過(guò)程,有助于理解揮發(fā)分F、Cl、S以及成礦元素本身的地球化學(xué)行為(Richards, 2015b)。鑒于巖漿磷灰石和共存熔體處于平衡狀態(tài),可使用磷灰石和平衡熔體實(shí)驗(yàn)分配系數(shù)從磷灰石成分估算在給定P-T條件下熔體的F和Cl含量(Stocketal., 2016)。近期研究表明,在磷灰石和熔體之間存在F和Cl的非能斯特分配,在確定熔體中F、Cl和/或H2O的濃度時(shí),各種實(shí)驗(yàn)結(jié)果和熱力學(xué)模型利用磷灰石和熔體之間不同的F和Cl交換系數(shù)(KD)進(jìn)行估算(Pan and Fleet, 2002; Rikeretal., 2018; Li and Costa, 2020)。然而,該方法要求已知熔體中的三種成分之一(F或Cl或H2O含量),以便使用交換系數(shù)確定巖漿的揮發(fā)分及水含量。但由于熔體中的F、Cl、H2O含量常需要利用較難找尋的熔體包裹體或火山玻璃進(jìn)行測(cè)定,該類研究對(duì)象較難找尋。因此,在本研究中,我們使用鹵族元素在磷灰石與熔體間的分配系數(shù)值來(lái)估算巖漿的F和Cl含量??紤]到巖漿成分(全巖成分熔體的摩爾A/CNK和N/NK)和估計(jì)的溫度-壓力條件,我們使用磷灰石與英安質(zhì)-流紋質(zhì)熔體之間的Cl和F分配系數(shù)(A/CNK為0.90~1.05,N/NK為0.60~0.77,Websteretal., 2009; Dohertyetal., 2014),其所代表的壓力和溫度范圍是50~200MPa和900~924℃。長(zhǎng)英質(zhì)硅酸鹽熔體的F濃度方程式為(Websteretal., 2009):CF(%)=[(磷灰石中的XF)-0.12]/3.02,其中XF表示磷灰石中的摩爾F含量。對(duì)于熔體中Cl含量的計(jì)算,我們使用Baoetal.(2016)的公式,該公式利用了Websteretal.(2009)和Dohertyetal.(2014)數(shù)據(jù)提供了一個(gè)簡(jiǎn)單的基于權(quán)重的能斯特分配系數(shù)方法:CCl(熔體,%)=CCl(磷灰石,%)×0.16052。我們利用該公式估算熔體Cl含量。對(duì)于巖漿的硫含量,巖漿磷灰石中SO3含量主要受溫度、氧化狀態(tài)和硅酸鹽熔體中硫含量的控制(Pengetal., 1997; Parat and Holtz, 2005)??紤]到熔體中硫含量的變化,早期Pengetal.(1997)利用與溫度相關(guān)的磷灰石-熔體分配系數(shù),開展巖漿相對(duì)硫含量的估算,公式為ln(S磷灰石/熔體(%))=21130/T(單位為開氏度)-16.2。Paratetal.(2011)總結(jié)前人發(fā)表的實(shí)驗(yàn)數(shù)據(jù),給出了SO3磷灰石(%)=0.157×lnSO3熔體(%)+0.9834(r2=0.68)。我們假設(shè)各巖漿磷灰石從硫酸鹽飽和熔體結(jié)晶,利用Paratetal.(2011)的方法估算巖漿硫含量。雖然兩種計(jì)算方法會(huì)獲得不同的巖漿硫含量,但仍可判斷其總體演化趨勢(shì)。

    圖10 土耳其西部斑巖礦床中磷灰石的F-Cl含量(a),及鋯石年齡和對(duì)應(yīng)巖體磷灰石估算的巖漿F(b)、Cl(c)和S(d)含量的演化關(guān)系其中巖漿中氟含量計(jì)算來(lái)自Webster et al. (2009),巖漿中氯含量計(jì)算來(lái)自Webster et al. (2009)和Doherty et al. (2014),巖漿中硫含量計(jì)算來(lái)自Parat et al. (2011),圖例見圖7Fig.10 The evolution relationship of apatite F-Cl content (a), zircon age and magmatic F (b), Cl (c), S (d) content calculated from apatite in the corresponding rock in western Turkish porphyry depositThe calculation of F content in magmatic refers to Webster et al. (2009), Cl content calculation refers to Webster et al. (2009) and Doherty et al. (2014), S content calculation refers to Parat et al. (2011). Legends refer to Fig.7

    斑巖成礦巖漿是否富硫,目前主要利用巖體中的巖漿礦物(巖漿硬石膏和磷灰石)和同期火山巖中熔融包裹體開展研究。如在成礦相關(guān)的斑巖中發(fā)現(xiàn)巖漿硬石膏,表明其形成時(shí)巖漿的氧逸度至少大于NNO+0.5,巖漿中硫含量為0.05%~1%(Carrol and Rutherford, 1987; Paratetal., 2011; Xiaoetal., 2012)。而利用高硫含量的磷灰石和硫在磷灰石和熔體間的分配系數(shù),計(jì)算平衡巖漿中的硫含量可達(dá)1000×10-6以上(Streck and Dilles, 1998; Liaghat and Tosdal, 2008),而通過(guò)礦床內(nèi)未發(fā)生硫化物飽和的基性端元中熔體包裹體測(cè)定硫含量為100×10-6~900×10-6(Seoetal., 2011; Grondahl and Zajacz, 2017; Zhang and Audétat, 2017),屬于正常的巖漿硫含量范圍。由于高氧逸度使得巖漿中的硫以+6價(jià)形式賦存,并記錄在結(jié)晶的磷灰石中(Paratetal., 2011),因此前人在俯沖與碰撞環(huán)境斑巖成礦帶區(qū)域?qū)Ρ妊芯恐邪l(fā)現(xiàn)成礦期磷灰石常具有較高的SO3含量。如中亞造山帶哈薩克斯坦地區(qū)斑巖銅礦研究(Caoetal., 2016)指出,俯沖作用可以加厚新生下地殼和巖石圈地幔,并使楔形地幔氧逸度不斷升高,較高的氧逸度使得成礦期巖石磷灰石較成礦前具有明顯較高的SO3含量。西藏班公湖-怒江帶俯沖期斑巖銅金礦(Lietal., 2021)和岡底斯后碰撞環(huán)境斑巖銅鉬礦(Tangetal., 2020)中,巖漿的高氧逸度使得成礦巖體中巖漿磷灰石具有較高的硫含量。上述成礦帶的研究與我們?cè)赥avanl帶、Biga半島斑巖成礦帶和Afyon-Konya成礦帶斑巖礦床中觀察到高氧逸度(圖5d)條件下的磷灰石高硫組成相一致。表明在俯沖、碰撞-后碰撞和伸展環(huán)境的斑巖礦床形成過(guò)程中,成礦期高氧逸度巖漿中的磷灰石常具有較高的SO3含量。

    5.3 觸發(fā)土耳其西部斑巖礦床形成的動(dòng)力學(xué)機(jī)制——與新特提斯單向匯聚過(guò)程有關(guān)

    根據(jù)斑巖礦床的形成時(shí)代、成礦巖漿的性質(zhì)及與區(qū)內(nèi)(或重大)地質(zhì)事件的時(shí)空耦合關(guān)系,可以表征斑巖礦床形成時(shí)的動(dòng)力學(xué)機(jī)制,進(jìn)而約束板塊之間俯沖、碰撞和陸內(nèi)伸展等一系列構(gòu)造過(guò)程及細(xì)節(jié)(Tosdal and Richards, 2001)。例如在洋-陸或洋-洋俯沖環(huán)境中,若成礦有關(guān)的埃達(dá)克質(zhì)巖石,與高鎂安山巖/富鈮玄武巖、A型花崗巖具有密切的時(shí)空關(guān)系,則其形成可能與洋脊俯沖有關(guān)(Cookeetal., 2005)。依據(jù)安第斯位于平板俯沖正上方或緩俯沖向陡俯沖過(guò)渡區(qū)域的眾多超大型斑巖銅礦建立的“平板俯沖模型”(Kay and Mpodozis, 2001; Reichetal., 2003),斑巖銅礦的形成可能與俯沖帶結(jié)構(gòu)和俯沖板片的幾何形態(tài)有關(guān),如俯沖板片上存在洋底高原(Gutscheretal., 2000; Rosenbaumetal., 2005),或俯沖板片上的轉(zhuǎn)換斷層(Richards and Holm, 2013)均可觸發(fā)斑巖銅礦的形成。俯沖帶極性反轉(zhuǎn)(Solomon, 1990)或者碰撞使得俯沖作用中止的階段(Sillitoe, 1997),能夠在陸緣弧后背景下形成富金斑巖銅礦。在增生和碰撞造山帶內(nèi),斑巖礦床形成的時(shí)空分布特征可以指示造山帶增生及小陸塊拼合過(guò)程(如中亞造山帶西段,秦克章, 2000; Wanetal., 2011, 2017)。通過(guò)斑巖礦床的形成時(shí)代和區(qū)域上成巖-成礦時(shí)代朝縫合帶方向變年輕的現(xiàn)象,可能說(shuō)明深部發(fā)生俯沖板片回轉(zhuǎn)過(guò)程(slab roll-back,如保加利亞地區(qū)俯沖相關(guān)礦床, Von Quadtetal., 2005; 西藏岡底斯主碰撞階段斑巖礦床, Zhaoetal., 2014)。在碰撞型斑巖銅礦則集中分布在后碰撞階段(Houetal., 2004)或擠壓向伸展轉(zhuǎn)換階段(Qinetal., 2005),可能與俯沖大陸板片斷離作用有關(guān)(slab break-off, De Boorderetal., 1998)。陸內(nèi)環(huán)境的斑巖銅金礦或者斑巖鉬礦,則可根據(jù)其形成時(shí)代判斷與陸內(nèi)伸展背景或者走滑轉(zhuǎn)換背景相關(guān)(Wangetal., 2006; Richards, 2009)。更為重要的是,不同的構(gòu)造動(dòng)力學(xué)背景之間可以相互疊置衍生出復(fù)合模型,用來(lái)解釋某些大型-超大型斑巖礦床及礦集區(qū)形成機(jī)制,如Gaoetal.(2018)利用“洋脊俯沖和俯沖改造的基性下地殼”模型解釋中亞造山帶的部分斑巖礦床形成機(jī)制。Dengetal.(2021)提出三江地區(qū)新生代大尺度巖石圈構(gòu)造變形——地殼差異性旋轉(zhuǎn)造成殼幔圈層非耦合變形,活化早期預(yù)富集的成礦物質(zhì),形成北衙和玉龍巨型斑巖-矽卡巖礦集區(qū)。由此可見,當(dāng)斑巖礦床形成于俯沖-碰撞-后碰撞一系列構(gòu)造演化過(guò)程時(shí),觸發(fā)斑巖礦床形成的動(dòng)力學(xué)機(jī)制往往是較為復(fù)雜的。

    鑒于特提斯構(gòu)造域鮮明的地質(zhì)演化特征是南側(cè)岡瓦納大陸不斷裂解、塊體向北漂移并與勞亞大陸持續(xù)聚合(吳福元等, 2020),該特征被形象稱為“特提斯號(hào)”單向列車(Wanetal., 2019)。而這種板塊單向裂解-聚合過(guò)程在土耳其-伊朗地區(qū)尤為突出(eng?r and Yilmaz, 1981; Stampfli and Borel, 2004)。數(shù)值模擬研究(Menantetal., 2016)顯示板塊動(dòng)力學(xué)和相關(guān)的軟流圈影響著該地區(qū)弧與弧后區(qū)域巖漿活動(dòng)的形成、運(yùn)移和排布,而形成的熔體將弱化巖石圈并促進(jìn)巖石圈的變形及其與軟流圈的耦合。該三維板塊動(dòng)力學(xué)模型強(qiáng)調(diào)俯沖動(dòng)力學(xué)過(guò)程與相關(guān)的軟流圈流動(dòng)對(duì)礦床在空間和時(shí)間上分布的影響,以及在伸展的弧后環(huán)境中變質(zhì)穹窿對(duì)含礦流體循環(huán)-集中的重要性(Menantetal., 2018)。由此可見,這種俯沖體制是控制包括土耳其在內(nèi)的白堊紀(jì)晚期以來(lái)地中海東部俯沖帶的地球動(dòng)力學(xué)和成礦演化重要?jiǎng)恿W(xué)因素。在該種向北單向俯沖體制控制下,土耳其西部斑巖型礦床形成時(shí)代具有向南逐漸變年輕的特點(diǎn)。這反映其形成與分布受控于碰撞后俯沖的新特提斯洋板片(Vardar洋)后撤-回轉(zhuǎn)過(guò)程(slab roll-back),即地幔楔熱結(jié)構(gòu)在板片后撤-板片回轉(zhuǎn)過(guò)程中發(fā)生改變,熱的軟流圈地?;亓饔|發(fā)地幔楔部分熔融,形成Tavanl帶與Biga半島斑巖成礦帶的斑巖型礦床。這種俯沖的洋殼發(fā)生板片回轉(zhuǎn)和斷裂過(guò)程,在土耳其以西的Apuseni-Banat-Timok-Srednogorie成礦帶上形成時(shí)代向南遷移島弧性質(zhì)的火山巖和花崗巖,von Quadtetal.(2005)認(rèn)為這些與鈣堿性巖漿有關(guān)的斑巖型Cu-Au礦和淺成低溫礦床可能與碰撞后俯沖板片的回轉(zhuǎn)過(guò)程有關(guān)。而在特提斯東段西藏岡底斯斑巖礦床中,僅在古新世-始新世巖漿活動(dòng)有關(guān)的斑巖礦床具有相似時(shí)代分布特征(Zhaoetal., 2014, 2016),并被認(rèn)為是與俯沖的新特提斯洋殼板片后撤的構(gòu)造過(guò)程相關(guān)。這一時(shí)期,陸陸碰撞造成俯沖的新特提斯洋板片發(fā)生回轉(zhuǎn),引發(fā)的軟流圈回流改變了地幔楔的熱結(jié)構(gòu),并誘發(fā)同時(shí)期藏南林子宗組島弧火山巖向南遷移(Leeetal., 2009)。由此可見,這種新特提斯洋向北俯沖過(guò)程中由于陸-陸碰撞形成板片后撤-板片回轉(zhuǎn)過(guò)程是特提斯中-東段形成的斑巖礦床的一類重要觸發(fā)機(jī)制。

    此后土耳其西部地區(qū)經(jīng)歷了地殼隆升過(guò)程(圖7和圖8),在中新世Afyon-Konya帶形成了以PnarbaMo-Cu礦和Au礦為代表的斑巖型礦床。新近研究指出(Rabayrol and Hart, 2021),在西安納托利亞中新世成礦作用主要集中在晚中新世(15~9Ma),之前是成礦較差的早中新世巖漿(21~18Ma),之后為貧礦的上新世至第四紀(jì)堿性巖漿。這種成礦集中期與中-東安納托利亞地區(qū)巖漿巖有明顯差別,即中東部早期為成礦巖漿,隨后巖漿向含水少、氧化程度低和以地幔組分為主的源區(qū)特征演化。這種整體上巖漿活動(dòng)特征與中新世俯沖的愛琴海與塞浦路斯板塊回轉(zhuǎn)-撕裂,及阿拉伯板塊斷裂過(guò)程有關(guān),該過(guò)程往往造成貧水的軟流圈地幔上涌和回流,進(jìn)而減薄了安納托利亞地區(qū)深部巖石圈地幔(Rabayrol and Hart, 2021)。此時(shí)控制該區(qū)域巖漿演化和斑巖礦床形成的動(dòng)力學(xué)條件已由北側(cè)Vardar洋俯沖消減造成的Tauride-Anatolide板塊碰撞,轉(zhuǎn)為由土耳其以南愛琴海板片俯沖控制的伸展區(qū)域。我們的研究指出,在這一中新世深部過(guò)程響應(yīng)下同時(shí)期土耳其西部地區(qū)淺部地殼經(jīng)歷了明顯減薄過(guò)程(圖8a, b),深部富集的巖石圈地幔減壓發(fā)生部分熔融(Richards, 2009; Bakeretal., 2016),形成的巖漿向著更高氧化程度演化(圖5d與圖8c中綠色與黃色),這類高氧逸度的巖漿條件使得硫元素更容易保存在熔體中(如最高的熔體硫估計(jì)量),為后期成礦提供有利條件。該動(dòng)力學(xué)推論與前人利用地震數(shù)據(jù)獲得土耳其西部地殼厚度(Zhuetal., 2006)與Moho面深度研究結(jié)果相一致(Karabulutetal., 2013),即從Izmir-Ankara-Erzincan縫合帶到Menders Massif再到地中海沿岸,地殼厚度從36km到28~30km減薄至25km,Moho面深度從32km到25km到20km。

    因此,我們認(rèn)為觸發(fā)土耳其西部新生代斑巖礦床形成的動(dòng)力學(xué)機(jī)制是:在新特提斯洋向北單向匯聚的背景下,北部始新世-漸新世斑巖礦床受控于碰撞后俯沖的新特提斯洋板片(Vardar洋)后撤-回轉(zhuǎn)-斷離過(guò)程;南部中新世斑巖礦床的形成受控于愛琴海板片俯沖控制的地殼伸展-減薄過(guò)程。

    5.4 土耳其西部斑巖型礦床形成的深部過(guò)程

    當(dāng)關(guān)注到斑巖礦床形成的深部殼幔相互作用時(shí),研究的著力點(diǎn)便集中到兩個(gè)概念模型:一個(gè)是地幔楔部分熔融形成的玄武質(zhì)巖漿MASH過(guò)程(熔融-同化-儲(chǔ)存-均一過(guò)程,Hildreth and Moorbath, 1988; Richards, 2003; Lietal., 2016),普遍用來(lái)解釋俯沖環(huán)境斑巖系統(tǒng)中酸性巖漿的形成;另一個(gè)是深部地殼熱區(qū)過(guò)程(Deep Crustal Hot Zone過(guò)程,簡(jiǎn)稱DCHZ過(guò)程,Annenetal., 2006; Kempetal., 2007),該模型強(qiáng)調(diào)巖漿成分取決于玄武質(zhì)巖漿的侵位深度及與其平衡的地溫梯度,其中富水玄武質(zhì)母巖漿形成于玄武巖分離結(jié)晶和地殼圍巖部分熔融過(guò)程(由冷卻玄武巖的熱和水造成, Annenetal., 2006),目前有學(xué)者利用該模式解釋特提斯中段碰撞環(huán)境斑巖礦床(Dengetal., 2018; Wanetal., 2018) 和我國(guó)華南晚中生代與花崗巖有關(guān)礦床(Wangetal., 2021) 形成的深部過(guò)程。但應(yīng)用這兩個(gè)概念模型有一定的適用條件:MASH模式需要有持續(xù)的俯沖作用(俯沖流體造成地幔楔部分熔融)并提供穩(wěn)定的熱和基性巖漿,以便在殼幔邊界附近停留和均一化。而DCHZ模式需要考慮到斑巖礦床的形成過(guò)程中發(fā)生基性巖漿注入過(guò)程(Hattori and Keith, 2001)。 因此,MASH模式更容易解釋北側(cè)Izmir-Ankara-Erzincan縫合帶附近的Tavanl與Biga半島斑巖成礦帶斑巖礦床的形成,因?yàn)榘讏准o(jì)-古新世的Pontides斑巖銅礦帶是由Vardar洋向北持續(xù)俯沖形成的,該時(shí)期巖漿巖的弧型巖漿特征2009; Kucuetal., 2019)表明至少在古新世該區(qū)域深部仍存在持續(xù)的俯沖作用。加之板片回轉(zhuǎn)-斷裂形成廣泛的中酸性巖漿作用和混合作用,指示深部MASH模式可能是控制該區(qū)域斑巖礦床形成的重要深部過(guò)程。而對(duì)于Afyon-Konya帶中Pnarba斑巖Mo-Cu礦(18~17Ma)和斑巖Au礦(16~12Ma),本研究指出形成該時(shí)期其高氧逸度成礦巖漿可能直接繼承于富集的巖石圈地幔部分熔融。伸展過(guò)程中地殼減薄造成富集地幔部分熔融形成玄武質(zhì)巖漿,并且該時(shí)期熔體Cl含量明顯高于Tavanl與Biga半島斑巖礦床中熔體Cl含量(圖10c),表明巖漿磷灰石結(jié)晶的熔體環(huán)境存在明顯差異。由于深部DCHZ帶內(nèi)底侵的玄武質(zhì)巖漿能夠同時(shí)為淺部巖漿房提供S和Cl(Zhuetal., 2018),本文更傾向于利用DCHZ模式解釋中新世Afyon-Konya帶斑巖礦床形成的深部過(guò)程。但在DCHZ過(guò)程中底侵玄武巖注入巖漿房的深度和時(shí)限則需要下一步更為精細(xì)的研究確定。

    6 結(jié)論

    (1)鋯石微量元素表明土耳其西部新生代斑巖礦床成礦巖漿大部分落于高水含量-高氧逸度區(qū)間;鋯石Ti溫度計(jì)估算土耳其西部斑巖成礦及其相關(guān)巖體的形成溫度在650~900℃之間,結(jié)晶溫度表現(xiàn)為斑巖Au礦>斑巖Cu-Mo=斑巖Mo-Cu礦>斑巖Mo礦,呈現(xiàn)遞減趨勢(shì);鋯石微量元素顯示(Yb/Gd和Hf/Y-Th/U)其組成大部分受角閃石±榍石±磷灰石分離結(jié)晶控制;從始新世到中新世,鋯石氧逸度Ce4+/Ce3+比值呈現(xiàn)出先降低(斑巖Cu-Mo到斑巖Mo)后升高(斑巖Mo-Cu到斑巖Au)的趨勢(shì),且斑巖Au和斑巖Cu-Mo相關(guān)巖體鋯石形成時(shí)氧逸度高于斑巖Mo和斑巖Mo-Cu礦巖體。

    (2)始新世-漸新世斑巖型Cu-Mo與Mo礦床中巖漿相對(duì)具有高F和低Cl組成,中新世伸展環(huán)境形成的斑巖Mo-Cu礦和斑巖Au礦中巖漿Cl含量普遍較高,其中與成礦巖體同期的暗色包體或基性巖脈中磷灰石計(jì)算獲得的巖漿硫含量均大于侵入體對(duì)應(yīng)巖漿的硫含量,表明基性巖漿注入可能為巖漿房提供硫。

    (3)通過(guò)區(qū)域動(dòng)力學(xué)分析和地殼厚度估算,本文認(rèn)為觸發(fā)土耳其西部新生代斑巖礦床形成的動(dòng)力學(xué)機(jī)制是:在新特提斯洋向北單向匯聚的背景下,北部始新世-漸新世斑巖礦床受控于碰撞后俯沖的新特提斯洋板片(Vardar洋)后撤-回轉(zhuǎn)-斷離過(guò)程;南部中新世斑巖礦床的形成受控于愛琴海板片俯沖控制的地殼伸展-減薄過(guò)程。北側(cè)Izmir-Ankara-Erzincan縫合帶附近的Tavanl與Biga半島斑巖成礦帶始新世-漸新世斑巖礦床的形成與熔融-同化-儲(chǔ)存-均一過(guò)程(MASH)有關(guān);深部地殼熱區(qū)過(guò)程(DCHZ)與中新世Afyon-Konya帶斑巖礦床有關(guān)。

    致謝野外工作得到中國(guó)科學(xué)院地質(zhì)與地球物理研究所李繼磊副研究員、褚楊副研究員、馮振天博士、劉霞博士、景揭俊博士的幫助;測(cè)試工作得到中國(guó)科學(xué)院地質(zhì)與地球物理研究所余可龍博士和賈立輝工程師的協(xié)助;成文過(guò)程受益于“特提斯地球動(dòng)力系統(tǒng)重大研究計(jì)劃”研究團(tuán)隊(duì)的學(xué)術(shù)交流,并得到了中國(guó)科學(xué)院地質(zhì)與地球物理研究所萬(wàn)博研究員、曹明堅(jiān)研究員、胡方泱副研究員和鄒心宇博士,中國(guó)科學(xué)院青藏高原研究所李金祥研究員的幫助、指導(dǎo)和啟發(fā);審稿人的寶貴修改意見和建議,以及期刊編輯的精心修改,使文章得以完善。在此筆者及其團(tuán)隊(duì)一并致以誠(chéng)摯的謝意。

    在筆者學(xué)生時(shí)期對(duì)斑巖系統(tǒng)構(gòu)造背景研究過(guò)程中,深得李繼亮老師的幫助和教導(dǎo),謹(jǐn)以此文紀(jì)念李繼亮老師。

    猜你喜歡
    逸度磷灰石斑巖
    地幔氧逸度的時(shí)空變化
    不同溫度、壓強(qiáng)、氧逸度條件下斜方輝石含水性的實(shí)驗(yàn)研究
    羥基磷灰石在鈾富集成礦中的作用
    濕法冶金(2019年5期)2019-10-18 09:00:00
    斑巖型礦床含礦斑巖與非含礦斑巖鑒定特征綜述
    幔源巖漿氧化還原狀態(tài)及對(duì)巖漿礦床成礦的制約*
    巖型礦床含礦斑巖與非含礦斑巖鑒定特征綜述
    煌斑巖的研究進(jìn)展
    氣體混合爐中氧逸度控制
    PHBV膜與珊瑚羥基磷灰石聯(lián)合修復(fù)頜骨缺損的研究
    濕化學(xué)法合成羥基磷灰石晶體及其表征
    国产亚洲一区二区精品| 哪个播放器可以免费观看大片| .国产精品久久| 这个男人来自地球电影免费观看 | 久久国产精品男人的天堂亚洲 | 五月开心婷婷网| 免费人妻精品一区二区三区视频| 97在线视频观看| 国产在视频线精品| 晚上一个人看的免费电影| 在线观看免费视频网站a站| 99视频精品全部免费 在线| 亚洲欧美日韩无卡精品| 美女国产视频在线观看| 久久6这里有精品| 99热这里只有是精品在线观看| 国产色婷婷99| 日韩欧美 国产精品| 日本wwww免费看| 天美传媒精品一区二区| 亚洲精品456在线播放app| 成年女人在线观看亚洲视频| 午夜福利网站1000一区二区三区| 爱豆传媒免费全集在线观看| 亚洲va在线va天堂va国产| 亚洲精品国产av蜜桃| 久久av网站| 精品人妻偷拍中文字幕| 日本猛色少妇xxxxx猛交久久| 一级毛片我不卡| 精品99又大又爽又粗少妇毛片| 日日摸夜夜添夜夜添av毛片| 日产精品乱码卡一卡2卡三| 亚洲国产欧美人成| 在线 av 中文字幕| 日韩电影二区| 日韩人妻高清精品专区| 国产精品欧美亚洲77777| 在现免费观看毛片| 国产色爽女视频免费观看| 女人久久www免费人成看片| 国产午夜精品久久久久久一区二区三区| 国国产精品蜜臀av免费| 97精品久久久久久久久久精品| 国产成人a区在线观看| 国产精品一区二区在线观看99| 哪个播放器可以免费观看大片| 久久综合国产亚洲精品| 国产视频内射| 99re6热这里在线精品视频| 国产亚洲5aaaaa淫片| 国产精品三级大全| 亚洲精品成人av观看孕妇| 精品一品国产午夜福利视频| 亚洲精品亚洲一区二区| a级毛片免费高清观看在线播放| 成年人午夜在线观看视频| 国产高清三级在线| 身体一侧抽搐| 麻豆国产97在线/欧美| 在线 av 中文字幕| 成人午夜精彩视频在线观看| 亚洲av国产av综合av卡| 中国国产av一级| 最近的中文字幕免费完整| 久热这里只有精品99| 1000部很黄的大片| 国国产精品蜜臀av免费| 韩国高清视频一区二区三区| 亚洲在久久综合| 国产精品偷伦视频观看了| 国产亚洲最大av| 亚洲欧美一区二区三区国产| 91精品一卡2卡3卡4卡| 亚洲精品日韩在线中文字幕| 中国国产av一级| 99热这里只有是精品在线观看| 欧美日本视频| 成人黄色视频免费在线看| 欧美性感艳星| 干丝袜人妻中文字幕| 91狼人影院| 中文字幕av成人在线电影| 日本猛色少妇xxxxx猛交久久| 国产精品久久久久久精品电影小说 | 丝袜喷水一区| 1000部很黄的大片| 亚洲精品国产色婷婷电影| 少妇 在线观看| 国产又色又爽无遮挡免| 日本色播在线视频| 久久 成人 亚洲| 精品99又大又爽又粗少妇毛片| 性高湖久久久久久久久免费观看| 99热网站在线观看| 热99国产精品久久久久久7| 亚洲第一av免费看| 久久热精品热| 亚洲国产精品一区三区| 亚洲国产成人一精品久久久| 黑丝袜美女国产一区| 亚洲综合色惰| 五月天丁香电影| 国产精品麻豆人妻色哟哟久久| 亚洲国产高清在线一区二区三| 亚洲精品色激情综合| 成年人午夜在线观看视频| 午夜免费观看性视频| 久久久亚洲精品成人影院| 男女国产视频网站| 男人狂女人下面高潮的视频| 欧美最新免费一区二区三区| 18禁裸乳无遮挡动漫免费视频| 免费观看的影片在线观看| 夜夜爽夜夜爽视频| 午夜福利在线在线| 国产av国产精品国产| 国产精品久久久久久久久免| 国产男女超爽视频在线观看| 大片电影免费在线观看免费| 麻豆乱淫一区二区| 欧美激情极品国产一区二区三区 | 日韩欧美一区视频在线观看 | 99九九线精品视频在线观看视频| 97在线人人人人妻| 国产黄色免费在线视频| 超碰av人人做人人爽久久| 嘟嘟电影网在线观看| 欧美最新免费一区二区三区| 欧美国产精品一级二级三级 | 欧美激情国产日韩精品一区| 高清日韩中文字幕在线| 午夜激情福利司机影院| 国产精品一区二区在线观看99| 亚州av有码| a 毛片基地| 91精品一卡2卡3卡4卡| 婷婷色综合大香蕉| 夜夜骑夜夜射夜夜干| h视频一区二区三区| 免费人妻精品一区二区三区视频| av卡一久久| tube8黄色片| 91午夜精品亚洲一区二区三区| 久久人人爽av亚洲精品天堂 | 国产成人精品婷婷| 国产无遮挡羞羞视频在线观看| 看十八女毛片水多多多| 日韩av免费高清视频| 久久人人爽人人片av| 国产成人freesex在线| 日本vs欧美在线观看视频 | 久久精品熟女亚洲av麻豆精品| 国产永久视频网站| 欧美日韩视频高清一区二区三区二| 国产爽快片一区二区三区| 在线天堂最新版资源| 日本猛色少妇xxxxx猛交久久| 久久精品久久久久久久性| 久久久久久久久久久免费av| 精品久久久久久久久av| 妹子高潮喷水视频| 在线天堂最新版资源| 黑人猛操日本美女一级片| 国产精品久久久久久久久免| 中文天堂在线官网| 少妇猛男粗大的猛烈进出视频| 亚洲人成网站高清观看| 中文欧美无线码| 嫩草影院入口| 蜜桃在线观看..| 日韩精品有码人妻一区| 麻豆成人午夜福利视频| 少妇高潮的动态图| 91久久精品国产一区二区三区| 嘟嘟电影网在线观看| 你懂的网址亚洲精品在线观看| 免费av中文字幕在线| 建设人人有责人人尽责人人享有的 | 精品久久久久久久久亚洲| 美女cb高潮喷水在线观看| 亚洲精品一二三| 欧美bdsm另类| 国产精品99久久99久久久不卡 | 国产精品人妻久久久久久| 黄片wwwwww| 毛片女人毛片| 免费大片黄手机在线观看| 各种免费的搞黄视频| 韩国av在线不卡| h日本视频在线播放| 色哟哟·www| 极品教师在线视频| 久久热精品热| 国产在视频线精品| 精品少妇黑人巨大在线播放| 高清av免费在线| 久久久久久久大尺度免费视频| 国产成人精品久久久久久| 国产精品精品国产色婷婷| 久久久久久久久久久丰满| 国产永久视频网站| 欧美成人a在线观看| 免费人成在线观看视频色| 精品久久久久久久久亚洲| 国产精品熟女久久久久浪| 欧美人与善性xxx| 黑人高潮一二区| 成年免费大片在线观看| 身体一侧抽搐| 99久久综合免费| 国产伦精品一区二区三区四那| 少妇丰满av| 日韩 亚洲 欧美在线| 国产黄片美女视频| 成人特级av手机在线观看| 欧美高清性xxxxhd video| 国产爱豆传媒在线观看| 亚洲欧美清纯卡通| 99国产精品免费福利视频| 免费不卡的大黄色大毛片视频在线观看| 国产乱来视频区| 国产高清不卡午夜福利| 视频中文字幕在线观看| 夜夜爽夜夜爽视频| 国产成人免费观看mmmm| 欧美日韩亚洲高清精品| 久久精品人妻少妇| videossex国产| 亚洲欧美成人精品一区二区| www.色视频.com| 亚洲精品成人av观看孕妇| 国产av精品麻豆| 偷拍熟女少妇极品色| 一个人免费看片子| 亚洲欧美成人综合另类久久久| 国产黄色视频一区二区在线观看| 高清午夜精品一区二区三区| 一二三四中文在线观看免费高清| 久久久久久久久久成人| 91久久精品国产一区二区三区| 又爽又黄a免费视频| av在线观看视频网站免费| 一本久久精品| www.av在线官网国产| 国产亚洲午夜精品一区二区久久| 国产成人a区在线观看| 美女主播在线视频| 国产精品一区www在线观看| 人人妻人人澡人人爽人人夜夜| 舔av片在线| 99热这里只有是精品50| 精品一区二区免费观看| 嘟嘟电影网在线观看| 在线观看av片永久免费下载| 2021少妇久久久久久久久久久| 不卡视频在线观看欧美| 麻豆国产97在线/欧美| 久久精品国产鲁丝片午夜精品| 国产精品偷伦视频观看了| 女性被躁到高潮视频| 大香蕉久久网| 成人漫画全彩无遮挡| 日本欧美视频一区| 视频区图区小说| a级毛色黄片| 亚洲av日韩在线播放| 国产成人a∨麻豆精品| 久久av网站| 免费久久久久久久精品成人欧美视频 | 久久鲁丝午夜福利片| 久久久精品94久久精品| 在线观看免费视频网站a站| 亚洲综合色惰| 亚洲国产精品国产精品| 街头女战士在线观看网站| 一级爰片在线观看| 菩萨蛮人人尽说江南好唐韦庄| 日本黄色片子视频| 亚洲av国产av综合av卡| 狂野欧美激情性bbbbbb| 国产爱豆传媒在线观看| 在线免费十八禁| 国产高清不卡午夜福利| 嫩草影院入口| h日本视频在线播放| 国产综合精华液| 日韩一区二区视频免费看| 丰满人妻一区二区三区视频av| 免费看日本二区| 午夜日本视频在线| 欧美精品亚洲一区二区| 永久免费av网站大全| 最后的刺客免费高清国语| 亚洲最大成人中文| 欧美日韩视频高清一区二区三区二| 久久这里有精品视频免费| 成人国产av品久久久| 校园人妻丝袜中文字幕| 欧美人与善性xxx| 激情五月婷婷亚洲| 国产av一区二区精品久久 | 中文字幕免费在线视频6| 久久99蜜桃精品久久| 精品一区二区三区视频在线| 九九在线视频观看精品| 嫩草影院入口| 免费av不卡在线播放| 日本爱情动作片www.在线观看| 亚洲人成网站在线观看播放| 亚洲欧洲日产国产| 秋霞伦理黄片| 国产精品三级大全| 人妻 亚洲 视频| 国产一区亚洲一区在线观看| 国产成人精品一,二区| 亚洲国产av新网站| 丝瓜视频免费看黄片| 亚洲成人中文字幕在线播放| 欧美 日韩 精品 国产| 免费人妻精品一区二区三区视频| 熟女人妻精品中文字幕| 91午夜精品亚洲一区二区三区| 亚洲精品日本国产第一区| 大陆偷拍与自拍| 亚洲aⅴ乱码一区二区在线播放| 亚洲性久久影院| 男女国产视频网站| h视频一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 国产成人精品久久久久久| 亚洲综合精品二区| 插逼视频在线观看| 我要看黄色一级片免费的| av福利片在线观看| 成人无遮挡网站| 黄色日韩在线| 18禁裸乳无遮挡免费网站照片| 五月开心婷婷网| 久久97久久精品| 国产成人freesex在线| 超碰av人人做人人爽久久| 国产精品人妻久久久久久| 多毛熟女@视频| 成人一区二区视频在线观看| 少妇的逼好多水| 精品一区二区三卡| 夫妻午夜视频| 蜜桃在线观看..| 91久久精品国产一区二区三区| 久久这里有精品视频免费| 九色成人免费人妻av| 中文字幕久久专区| 少妇的逼好多水| 精品视频人人做人人爽| 日本猛色少妇xxxxx猛交久久| 中文字幕久久专区| 夫妻午夜视频| 国产精品不卡视频一区二区| 久久综合国产亚洲精品| 国产成人freesex在线| 视频中文字幕在线观看| 岛国毛片在线播放| 精品久久久久久久末码| 亚洲国产高清在线一区二区三| 婷婷色综合www| 你懂的网址亚洲精品在线观看| 日韩免费高清中文字幕av| 国产av精品麻豆| 亚洲婷婷狠狠爱综合网| 大陆偷拍与自拍| 欧美日韩国产mv在线观看视频 | 少妇人妻精品综合一区二区| 只有这里有精品99| 国产v大片淫在线免费观看| 成人黄色视频免费在线看| 免费看日本二区| 哪个播放器可以免费观看大片| 国产精品久久久久久久久免| 亚洲欧美成人综合另类久久久| 国产久久久一区二区三区| 亚洲欧洲国产日韩| 精品人妻一区二区三区麻豆| 亚州av有码| 夜夜爽夜夜爽视频| 免费av中文字幕在线| 人人妻人人看人人澡| 精品一区在线观看国产| 久久久久久久久久成人| 国产免费视频播放在线视频| 免费av不卡在线播放| 日韩欧美精品免费久久| av在线蜜桃| av女优亚洲男人天堂| 国产爱豆传媒在线观看| 天堂俺去俺来也www色官网| 欧美日韩在线观看h| 一级av片app| 哪个播放器可以免费观看大片| 亚洲国产成人一精品久久久| 日本vs欧美在线观看视频 | 在线看a的网站| 国产人妻一区二区三区在| 欧美精品人与动牲交sv欧美| 久久精品国产亚洲网站| 精品久久久久久久久亚洲| 欧美成人午夜免费资源| 国产欧美另类精品又又久久亚洲欧美| 在线观看免费视频网站a站| 一级毛片电影观看| 内射极品少妇av片p| 日日摸夜夜添夜夜爱| 欧美最新免费一区二区三区| 在现免费观看毛片| 国产黄色免费在线视频| 如何舔出高潮| 亚洲综合色惰| 久久婷婷青草| 在线观看美女被高潮喷水网站| 久久鲁丝午夜福利片| 欧美成人a在线观看| 国产有黄有色有爽视频| 国产精品一区www在线观看| 午夜福利在线在线| 亚洲av不卡在线观看| 少妇人妻 视频| 涩涩av久久男人的天堂| 亚洲内射少妇av| 欧美变态另类bdsm刘玥| 色视频在线一区二区三区| 国产成人精品福利久久| 美女cb高潮喷水在线观看| 亚洲精品aⅴ在线观看| 亚洲精品成人av观看孕妇| 亚洲欧美日韩无卡精品| 亚洲av在线观看美女高潮| 亚洲一级一片aⅴ在线观看| 国产淫语在线视频| 色哟哟·www| 亚洲无线观看免费| 91精品国产国语对白视频| 久久久久久久久久人人人人人人| 一级毛片我不卡| 熟女电影av网| 卡戴珊不雅视频在线播放| 91狼人影院| av在线老鸭窝| 欧美日韩在线观看h| 亚洲国产色片| h日本视频在线播放| 少妇裸体淫交视频免费看高清| h视频一区二区三区| 中国三级夫妇交换| 中文字幕亚洲精品专区| 国产精品国产三级专区第一集| 日韩亚洲欧美综合| 成人亚洲欧美一区二区av| 久久青草综合色| 国产精品不卡视频一区二区| 九九久久精品国产亚洲av麻豆| 欧美日韩亚洲高清精品| 我的女老师完整版在线观看| 人妻 亚洲 视频| 久久久久久久精品精品| 99精国产麻豆久久婷婷| 久久久久国产精品人妻一区二区| 婷婷色综合大香蕉| 亚洲人成网站高清观看| 国产精品精品国产色婷婷| 成人免费观看视频高清| 热99国产精品久久久久久7| freevideosex欧美| 国产欧美日韩一区二区三区在线 | 1000部很黄的大片| 亚洲成人一二三区av| 久久久成人免费电影| 国产淫语在线视频| 日本爱情动作片www.在线观看| 99精国产麻豆久久婷婷| 少妇人妻精品综合一区二区| 视频区图区小说| 夜夜爽夜夜爽视频| 午夜免费鲁丝| 免费观看av网站的网址| 国产乱来视频区| 搡女人真爽免费视频火全软件| 精品一区二区三卡| 天美传媒精品一区二区| 日本免费在线观看一区| 久久99热6这里只有精品| 搡女人真爽免费视频火全软件| 联通29元200g的流量卡| 亚洲国产精品一区三区| 大陆偷拍与自拍| 亚洲人与动物交配视频| 99热国产这里只有精品6| 欧美一级a爱片免费观看看| 国产精品精品国产色婷婷| 国产精品一区二区在线观看99| 久久久久精品性色| 日韩中字成人| 亚洲第一av免费看| 日本av手机在线免费观看| 亚洲成人av在线免费| 下体分泌物呈黄色| av国产免费在线观看| 在线观看一区二区三区激情| 丰满迷人的少妇在线观看| 成人亚洲欧美一区二区av| 亚洲人与动物交配视频| 国产精品av视频在线免费观看| 精品一品国产午夜福利视频| 欧美zozozo另类| 我的老师免费观看完整版| 嘟嘟电影网在线观看| 日韩中文字幕视频在线看片 | 97热精品久久久久久| 天天躁夜夜躁狠狠久久av| 欧美极品一区二区三区四区| 欧美精品国产亚洲| 亚洲在久久综合| 2021少妇久久久久久久久久久| 亚洲欧美精品专区久久| 国产一区有黄有色的免费视频| 美女xxoo啪啪120秒动态图| 精品国产露脸久久av麻豆| 麻豆精品久久久久久蜜桃| 亚洲综合色惰| 久久亚洲国产成人精品v| 99久久综合免费| 亚洲一区二区三区欧美精品| 久久久久久久久久久丰满| 亚洲精品国产色婷婷电影| 免费av不卡在线播放| 久久久久国产精品人妻一区二区| 亚洲怡红院男人天堂| 人妻 亚洲 视频| 国产精品久久久久久av不卡| 男人舔奶头视频| 亚洲av在线观看美女高潮| 精品一区在线观看国产| 国产av精品麻豆| 国产精品一二三区在线看| 中国三级夫妇交换| 内射极品少妇av片p| 久久久色成人| 一区在线观看完整版| 亚洲三级黄色毛片| 亚洲色图综合在线观看| 国产永久视频网站| 久久精品国产a三级三级三级| 国产视频首页在线观看| 成人亚洲欧美一区二区av| 青春草国产在线视频| 亚洲综合色惰| 人人妻人人爽人人添夜夜欢视频 | 久久精品久久久久久久性| 久久久久视频综合| 久久鲁丝午夜福利片| 国产精品一区www在线观看| 一区二区av电影网| 欧美老熟妇乱子伦牲交| 2021少妇久久久久久久久久久| 亚洲成人中文字幕在线播放| h视频一区二区三区| 大片免费播放器 马上看| 日产精品乱码卡一卡2卡三| 老女人水多毛片| 一级毛片久久久久久久久女| 国产精品一区二区性色av| 狠狠精品人妻久久久久久综合| 丰满乱子伦码专区| 国产熟女欧美一区二区| 美女视频免费永久观看网站| 国产精品三级大全| 99精国产麻豆久久婷婷| 国产精品久久久久久精品古装| 日韩av免费高清视频| 国产精品一区二区在线观看99| 建设人人有责人人尽责人人享有的 | 97热精品久久久久久| 亚洲欧美精品专区久久| 大又大粗又爽又黄少妇毛片口| 女人十人毛片免费观看3o分钟| 亚洲精品乱码久久久久久按摩| 精品人妻一区二区三区麻豆| 久久精品国产鲁丝片午夜精品| 丰满迷人的少妇在线观看| 寂寞人妻少妇视频99o| 亚洲av不卡在线观看| 97在线视频观看| 国产熟女欧美一区二区| 欧美成人精品欧美一级黄| 一级毛片电影观看| av.在线天堂| 欧美极品一区二区三区四区| 成人影院久久| 亚洲第一区二区三区不卡| 纯流量卡能插随身wifi吗| 亚洲精品国产成人久久av| 国产熟女欧美一区二区| 18禁动态无遮挡网站| 激情五月婷婷亚洲| 精品熟女少妇av免费看| 国产爽快片一区二区三区| 成人漫画全彩无遮挡| 纵有疾风起免费观看全集完整版| 久久精品国产鲁丝片午夜精品| 日韩成人av中文字幕在线观看| 91久久精品电影网| 中文在线观看免费www的网站| 国产av国产精品国产| 国产精品国产三级专区第一集| 精品一区二区三卡| 欧美日韩国产mv在线观看视频 |