• <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é)法合成羥基磷灰石晶體及其表征
    国产精品九九99| 国产精品自产拍在线观看55亚洲 | 精品国产超薄肉色丝袜足j| 日本91视频免费播放| 亚洲精品国产精品久久久不卡| 中国美女看黄片| 婷婷丁香在线五月| 一区二区三区乱码不卡18| 欧美成人午夜精品| 久久精品国产亚洲av香蕉五月 | 啦啦啦 在线观看视频| 人成视频在线观看免费观看| 国产男女内射视频| a 毛片基地| 欧美 日韩 精品 国产| 久久精品熟女亚洲av麻豆精品| 宅男免费午夜| a在线观看视频网站| 午夜久久久在线观看| 高清欧美精品videossex| 91字幕亚洲| 欧美午夜高清在线| 久久99一区二区三区| av福利片在线| 91精品国产国语对白视频| 动漫黄色视频在线观看| 一本一本久久a久久精品综合妖精| 中亚洲国语对白在线视频| 人妻久久中文字幕网| 人妻人人澡人人爽人人| 日本五十路高清| 国产精品.久久久| 在线天堂中文资源库| 久久免费观看电影| 久久香蕉激情| 国产成人影院久久av| 男女免费视频国产| 黄色视频,在线免费观看| 久久久久网色| 亚洲 欧美一区二区三区| 欧美日韩一级在线毛片| 中文字幕人妻熟女乱码| 狠狠婷婷综合久久久久久88av| 99国产综合亚洲精品| videosex国产| 久久毛片免费看一区二区三区| 在线亚洲精品国产二区图片欧美| 国产在线视频一区二区| 日本av免费视频播放| 久久热在线av| 国产高清视频在线播放一区 | 午夜免费鲁丝| 999久久久国产精品视频| 久久久久精品人妻al黑| 午夜老司机福利片| 欧美人与性动交α欧美软件| 欧美人与性动交α欧美软件| 在线观看免费视频网站a站| 国产福利在线免费观看视频| 午夜成年电影在线免费观看| 一边摸一边做爽爽视频免费| 在线精品无人区一区二区三| 女警被强在线播放| 国产1区2区3区精品| 午夜视频精品福利| 香蕉丝袜av| 国产成人免费观看mmmm| 无限看片的www在线观看| 999久久久精品免费观看国产| 久久av网站| 国产色视频综合| 成人国产av品久久久| 一级毛片女人18水好多| 在线观看舔阴道视频| 一本综合久久免费| 国产视频一区二区在线看| 久久精品熟女亚洲av麻豆精品| 老司机影院毛片| 一区在线观看完整版| 成人三级做爰电影| 咕卡用的链子| 国产精品久久久久久精品古装| 另类亚洲欧美激情| 欧美激情高清一区二区三区| 18禁黄网站禁片午夜丰满| 久久久国产一区二区| 午夜福利视频精品| 欧美激情久久久久久爽电影 | 欧美另类亚洲清纯唯美| 高清av免费在线| 男女高潮啪啪啪动态图| 日日摸夜夜添夜夜添小说| 国产亚洲一区二区精品| 男女床上黄色一级片免费看| 中文字幕人妻熟女乱码| 精品国产国语对白av| 99国产极品粉嫩在线观看| 久久女婷五月综合色啪小说| 亚洲国产av新网站| 日本av手机在线免费观看| 亚洲国产精品成人久久小说| 在线观看免费视频网站a站| 亚洲激情五月婷婷啪啪| av在线播放精品| 亚洲av成人不卡在线观看播放网 | 国产精品久久久久久精品电影小说| 男女无遮挡免费网站观看| 国产精品 国内视频| www.熟女人妻精品国产| 中文字幕制服av| 成人亚洲精品一区在线观看| 欧美亚洲 丝袜 人妻 在线| 色视频在线一区二区三区| 飞空精品影院首页| 国产精品.久久久| 热99re8久久精品国产| 国产1区2区3区精品| 精品人妻一区二区三区麻豆| 老鸭窝网址在线观看| 色综合欧美亚洲国产小说| 亚洲专区国产一区二区| 啦啦啦 在线观看视频| 精品久久久久久电影网| 一边摸一边抽搐一进一出视频| 一区二区日韩欧美中文字幕| 久久久国产成人免费| 国产1区2区3区精品| 国产亚洲精品第一综合不卡| 久久免费观看电影| 在线观看免费高清a一片| 十八禁高潮呻吟视频| 男女之事视频高清在线观看| 一级片'在线观看视频| 国产亚洲av高清不卡| 99re6热这里在线精品视频| 精品熟女少妇八av免费久了| 欧美激情久久久久久爽电影 | 日本av免费视频播放| 久久国产亚洲av麻豆专区| 日本五十路高清| 久久99一区二区三区| 国产高清国产精品国产三级| 国产深夜福利视频在线观看| 精品国产乱码久久久久久男人| 高清av免费在线| 亚洲欧美激情在线| 国产av精品麻豆| av欧美777| 老鸭窝网址在线观看| 国产片内射在线| av免费在线观看网站| 欧美在线黄色| 免费女性裸体啪啪无遮挡网站| 法律面前人人平等表现在哪些方面 | 在线观看免费视频网站a站| 国产精品成人在线| 久久国产亚洲av麻豆专区| www.自偷自拍.com| 精品国产一区二区三区久久久樱花| av又黄又爽大尺度在线免费看| 国产精品99久久99久久久不卡| 深夜精品福利| 精品久久蜜臀av无| 国产精品国产三级国产专区5o| 成年人免费黄色播放视频| 免费不卡黄色视频| 亚洲三区欧美一区| 亚洲欧美一区二区三区黑人| 久久久久久久精品精品| 91成人精品电影| 电影成人av| 丰满人妻熟妇乱又伦精品不卡| 黄色 视频免费看| 成在线人永久免费视频| av天堂久久9| 国产野战对白在线观看| 人人妻人人澡人人爽人人夜夜| 亚洲成国产人片在线观看| 国产精品熟女久久久久浪| 久久精品aⅴ一区二区三区四区| 麻豆国产av国片精品| 高清欧美精品videossex| 精品人妻在线不人妻| 国产成人精品无人区| 国产欧美日韩精品亚洲av| 少妇粗大呻吟视频| 亚洲五月色婷婷综合| 亚洲全国av大片| 欧美黑人精品巨大| 少妇被粗大的猛进出69影院| 国产福利在线免费观看视频| 亚洲精品乱久久久久久| 国产人伦9x9x在线观看| 一本综合久久免费| 纯流量卡能插随身wifi吗| 欧美+亚洲+日韩+国产| 久久精品熟女亚洲av麻豆精品| 免费观看av网站的网址| 少妇的丰满在线观看| 少妇猛男粗大的猛烈进出视频| 欧美成人午夜精品| 国产成人系列免费观看| 成人亚洲精品一区在线观看| 午夜精品久久久久久毛片777| 人妻久久中文字幕网| 人妻人人澡人人爽人人| 男女午夜视频在线观看| 国产男女超爽视频在线观看| 99热国产这里只有精品6| 亚洲色图综合在线观看| 国产亚洲一区二区精品| 巨乳人妻的诱惑在线观看| 成在线人永久免费视频| 日韩大片免费观看网站| 一级a爱视频在线免费观看| 亚洲精品粉嫩美女一区| 欧美一级毛片孕妇| 乱人伦中国视频| 久久99一区二区三区| av在线播放精品| avwww免费| 午夜精品久久久久久毛片777| 久久久久网色| 人人妻人人澡人人爽人人夜夜| 男女国产视频网站| 亚洲欧美色中文字幕在线| 两性夫妻黄色片| 午夜福利影视在线免费观看| 日韩熟女老妇一区二区性免费视频| 久久久久国产一级毛片高清牌| 另类精品久久| 亚洲av国产av综合av卡| 国产精品久久久久久精品电影小说| 美女高潮到喷水免费观看| 高清av免费在线| 伊人久久大香线蕉亚洲五| 伊人亚洲综合成人网| 91精品伊人久久大香线蕉| 国产欧美日韩一区二区三 | 黑人操中国人逼视频| 国产精品国产av在线观看| 国产麻豆69| 午夜两性在线视频| 纵有疾风起免费观看全集完整版| 丰满少妇做爰视频| 在线天堂中文资源库| 亚洲欧洲日产国产| 久久久国产欧美日韩av| 精品国产乱码久久久久久小说| 狂野欧美激情性bbbbbb| 亚洲一码二码三码区别大吗| 99久久99久久久精品蜜桃| 久久av网站| 成人三级做爰电影| tube8黄色片| 波多野结衣一区麻豆| 欧美激情久久久久久爽电影 | 极品人妻少妇av视频| 建设人人有责人人尽责人人享有的| 国产有黄有色有爽视频| 日日摸夜夜添夜夜添小说| 国产精品久久久人人做人人爽| 日韩大片免费观看网站| 人人妻,人人澡人人爽秒播| 91国产中文字幕| 美国免费a级毛片| 日本精品一区二区三区蜜桃| 欧美大码av| 国产一区有黄有色的免费视频| 久9热在线精品视频| 老熟妇乱子伦视频在线观看 | 一级毛片女人18水好多| 亚洲男人天堂网一区| 熟女少妇亚洲综合色aaa.| 国产深夜福利视频在线观看| 99热国产这里只有精品6| 午夜两性在线视频| 免费看十八禁软件| a级毛片黄视频| 99精国产麻豆久久婷婷| 美女福利国产在线| 在线天堂中文资源库| 国产人伦9x9x在线观看| 国产免费av片在线观看野外av| 国产区一区二久久| 韩国精品一区二区三区| 一进一出抽搐动态| 妹子高潮喷水视频| 国产色视频综合| 精品福利观看| 国产精品一区二区在线观看99| svipshipincom国产片| 国产亚洲精品第一综合不卡| 色婷婷av一区二区三区视频| 亚洲 国产 在线| 成人手机av| 亚洲午夜精品一区,二区,三区| 亚洲九九香蕉| 亚洲av片天天在线观看| 日韩熟女老妇一区二区性免费视频| 久久人人爽av亚洲精品天堂| 国产一卡二卡三卡精品| 久久ye,这里只有精品| 国产日韩欧美在线精品| 在线av久久热| 极品少妇高潮喷水抽搐| 日本wwww免费看| 欧美日韩黄片免| 久久天躁狠狠躁夜夜2o2o| 亚洲av成人一区二区三| 99精品欧美一区二区三区四区| 国产精品 欧美亚洲| 在线观看一区二区三区激情| 一本色道久久久久久精品综合| 亚洲 国产 在线| 亚洲天堂av无毛| 亚洲成人免费电影在线观看| 母亲3免费完整高清在线观看| 高清欧美精品videossex| 精品久久久久久久毛片微露脸 | 国产欧美日韩一区二区三 | 欧美av亚洲av综合av国产av| 亚洲免费av在线视频| 亚洲精品国产av蜜桃| 国产av又大| 一级黄色大片毛片| 中亚洲国语对白在线视频| 不卡av一区二区三区| 国产一区二区三区在线臀色熟女 | 99国产精品一区二区蜜桃av | 中文字幕高清在线视频| 亚洲中文av在线| 国产精品一区二区免费欧美 | 国产1区2区3区精品| 老汉色av国产亚洲站长工具| 一区二区三区四区激情视频| 五月开心婷婷网| 国产99久久九九免费精品| 99精品欧美一区二区三区四区| 制服人妻中文乱码| 久久精品aⅴ一区二区三区四区| 老司机福利观看| 天天影视国产精品| 777米奇影视久久| av在线app专区| 三级毛片av免费| 18禁观看日本| 少妇人妻久久综合中文| 欧美+亚洲+日韩+国产| 欧美成狂野欧美在线观看| 不卡av一区二区三区| 免费观看人在逋| 午夜福利影视在线免费观看| 精品免费久久久久久久清纯 | 午夜老司机福利片| 精品人妻1区二区| 亚洲国产av影院在线观看| 黄色视频在线播放观看不卡| 色94色欧美一区二区| 正在播放国产对白刺激| 纯流量卡能插随身wifi吗| 成人亚洲精品一区在线观看| 精品欧美一区二区三区在线| a 毛片基地| 国产欧美日韩精品亚洲av| 爱豆传媒免费全集在线观看| 在线观看www视频免费| 亚洲国产看品久久| 久久精品亚洲av国产电影网| 午夜成年电影在线免费观看| 国产高清视频在线播放一区 | 亚洲欧美清纯卡通| 男人添女人高潮全过程视频| 超色免费av| 日日夜夜操网爽| 国产精品影院久久| 最近最新免费中文字幕在线| 99国产精品99久久久久| 国产真人三级小视频在线观看| 男女免费视频国产| 黄色a级毛片大全视频| 91成人精品电影| 免费在线观看黄色视频的| 久久久久久人人人人人| 一进一出抽搐动态| 国产日韩欧美亚洲二区| 一本—道久久a久久精品蜜桃钙片| 乱人伦中国视频| 亚洲av男天堂| 男人添女人高潮全过程视频| 正在播放国产对白刺激| 亚洲久久久国产精品| 脱女人内裤的视频| 午夜福利免费观看在线| 我要看黄色一级片免费的| 精品免费久久久久久久清纯 | 女人爽到高潮嗷嗷叫在线视频| 国产av精品麻豆| 一本大道久久a久久精品| 国产有黄有色有爽视频| 亚洲国产精品成人久久小说| 最新的欧美精品一区二区| 97人妻天天添夜夜摸| 电影成人av| 亚洲精品成人av观看孕妇| 欧美日韩成人在线一区二区| 999精品在线视频| 99热国产这里只有精品6| 日韩制服丝袜自拍偷拍| 久久这里只有精品19| 十八禁高潮呻吟视频| 18禁裸乳无遮挡动漫免费视频| av在线老鸭窝| 午夜久久久在线观看| 国产不卡av网站在线观看| 黄色a级毛片大全视频| 岛国毛片在线播放| netflix在线观看网站| 国产精品av久久久久免费| 亚洲中文日韩欧美视频| 91av网站免费观看| 水蜜桃什么品种好| 香蕉国产在线看| 免费在线观看影片大全网站| 日本vs欧美在线观看视频| 天天影视国产精品| 欧美精品av麻豆av| 欧美午夜高清在线| 亚洲国产中文字幕在线视频| 亚洲欧美日韩另类电影网站| 又紧又爽又黄一区二区| 美女高潮到喷水免费观看| 高清在线国产一区| 91字幕亚洲| 国产av又大| 侵犯人妻中文字幕一二三四区| 精品免费久久久久久久清纯 | 操出白浆在线播放| 老熟妇仑乱视频hdxx| 纵有疾风起免费观看全集完整版| 咕卡用的链子| 女人久久www免费人成看片| 久久热在线av| 日本精品一区二区三区蜜桃| 久久综合国产亚洲精品| 一区福利在线观看| 三级毛片av免费| 成人av一区二区三区在线看 | 亚洲av片天天在线观看| 欧美激情 高清一区二区三区| 国产精品免费大片| 久久人人97超碰香蕉20202| 日韩大码丰满熟妇| 国产精品熟女久久久久浪| 亚洲人成77777在线视频| 国产99久久九九免费精品| 国产免费av片在线观看野外av| 老司机午夜福利在线观看视频 | 一个人免费看片子| av网站免费在线观看视频| 老司机在亚洲福利影院| 后天国语完整版免费观看| 亚洲专区中文字幕在线| 婷婷丁香在线五月| 亚洲男人天堂网一区| 美女福利国产在线| 永久免费av网站大全| 一区二区三区激情视频| 18禁黄网站禁片午夜丰满| 欧美精品高潮呻吟av久久| 亚洲欧洲日产国产| 国产欧美亚洲国产| 热99国产精品久久久久久7| 久久久久国内视频| 国产欧美日韩一区二区三区在线| 老司机午夜福利在线观看视频 | 90打野战视频偷拍视频| 热99国产精品久久久久久7| 国产黄色免费在线视频| 纯流量卡能插随身wifi吗| 各种免费的搞黄视频| 久久天堂一区二区三区四区| 美女福利国产在线| 久久久久久久久久久久大奶| 午夜视频精品福利| 色综合欧美亚洲国产小说| 免费在线观看完整版高清| 亚洲国产欧美网| 亚洲av欧美aⅴ国产| 国产成人欧美| 久久久久国产一级毛片高清牌| 国产一卡二卡三卡精品| 国产免费一区二区三区四区乱码| 久久精品国产亚洲av高清一级| 建设人人有责人人尽责人人享有的| 一区二区av电影网| 日韩欧美一区视频在线观看| 啦啦啦视频在线资源免费观看| 91av网站免费观看| 丰满饥渴人妻一区二区三| 老司机午夜十八禁免费视频| 另类亚洲欧美激情| 亚洲国产欧美网| 一级片免费观看大全| 欧美亚洲日本最大视频资源| 成年女人毛片免费观看观看9 | 黄色 视频免费看| 国产精品亚洲av一区麻豆| 91麻豆精品激情在线观看国产 | 电影成人av| 纯流量卡能插随身wifi吗| 中文精品一卡2卡3卡4更新| 日韩视频在线欧美| 美女扒开内裤让男人捅视频| 亚洲av成人不卡在线观看播放网 | 亚洲欧美一区二区三区黑人| 久久青草综合色| 欧美黑人欧美精品刺激| 9色porny在线观看| kizo精华| av一本久久久久| 国产又色又爽无遮挡免| 国产片内射在线| 午夜免费观看性视频| 免费久久久久久久精品成人欧美视频| 亚洲国产毛片av蜜桃av| 午夜老司机福利片| 日韩一区二区三区影片| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩国产mv在线观看视频| 亚洲欧美色中文字幕在线| 色综合欧美亚洲国产小说| 波多野结衣av一区二区av| 菩萨蛮人人尽说江南好唐韦庄| 成年av动漫网址| 欧美成狂野欧美在线观看| 激情视频va一区二区三区| 制服诱惑二区| 国产亚洲一区二区精品| 欧美精品一区二区免费开放| 999精品在线视频| 国产精品亚洲av一区麻豆| 午夜福利,免费看| 欧美 亚洲 国产 日韩一| 欧美日韩视频精品一区| 青春草视频在线免费观看| 日日摸夜夜添夜夜添小说| 久久久久久久久久久久大奶| 午夜精品国产一区二区电影| 免费av中文字幕在线| 成人免费观看视频高清| 久久国产精品影院| 菩萨蛮人人尽说江南好唐韦庄| 欧美变态另类bdsm刘玥| 天天躁日日躁夜夜躁夜夜| 国产黄频视频在线观看| 久久久久久久大尺度免费视频| 国产精品秋霞免费鲁丝片| 久久精品国产a三级三级三级| 精品久久蜜臀av无| 国产又色又爽无遮挡免| 天堂8中文在线网| 大香蕉久久网| 国产精品久久久久成人av| 美女大奶头黄色视频| 麻豆乱淫一区二区| 欧美成狂野欧美在线观看| 老熟妇仑乱视频hdxx| 精品一区二区三区四区五区乱码| 精品亚洲成国产av| 国产男女内射视频| 亚洲成人国产一区在线观看| 精品一区二区三区av网在线观看 | 国产野战对白在线观看| 2018国产大陆天天弄谢| 宅男免费午夜| 久久 成人 亚洲| 日韩,欧美,国产一区二区三区| 美女福利国产在线| 免费女性裸体啪啪无遮挡网站| 五月天丁香电影| 欧美亚洲 丝袜 人妻 在线| 久久久久久人人人人人| 亚洲伊人久久精品综合| tube8黄色片| 午夜久久久在线观看| 国产一卡二卡三卡精品| 国产成人免费无遮挡视频| 肉色欧美久久久久久久蜜桃| 国产精品久久久久久精品电影小说| 亚洲精品自拍成人| 亚洲国产欧美日韩在线播放| 超色免费av| 久久人人97超碰香蕉20202| 国产精品一区二区在线不卡| 亚洲精品中文字幕一二三四区 | 五月开心婷婷网| 黄网站色视频无遮挡免费观看| 国产男人的电影天堂91| 午夜福利免费观看在线| 亚洲精品在线美女| 欧美中文综合在线视频| 丝袜美腿诱惑在线| 欧美日韩亚洲高清精品| 国产高清videossex| 天堂中文最新版在线下载| 久久人人爽人人片av| 一进一出抽搐动态| 在线观看一区二区三区激情| 国产精品一区二区在线不卡| 女人久久www免费人成看片| 免费在线观看视频国产中文字幕亚洲 | 久久精品国产综合久久久| 国产黄频视频在线观看| 在线观看免费视频网站a站|