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

    新疆西準(zhǔn)噶爾蘇云河斑巖型鉬礦的剝露和保存條件:來(lái)自裂變徑跡和(U-Th)/He熱年代學(xué)的約束*

    2021-09-24 05:42:22王雅美尹繼元袁超肖文交陳文蔡克大楊雪葉陶再禮
    巖石學(xué)報(bào) 2021年8期
    關(guān)鍵詞:徑跡準(zhǔn)噶爾鉬礦

    王雅美 尹繼元 袁超 肖文交 陳文 蔡克大 楊雪葉 陶再禮

    1. 中國(guó)地質(zhì)大學(xué)(北京)地球科學(xué)與資源學(xué)院,北京 100083

    2. 中國(guó)地質(zhì)科學(xué)院地質(zhì)研究所,自然資源部深地動(dòng)力學(xué)重點(diǎn)實(shí)驗(yàn)室,北京 100037

    3. 中國(guó)科學(xué)院廣州地球化學(xué)研究所,同位素地球化學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,廣州 510640

    4. 中國(guó)科學(xué)院新疆生態(tài)與地理研究所新疆礦產(chǎn)資源研究中心,烏魯木齊 830011

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

    斑巖型鉬礦床是全球Mo金屬最主要的來(lái)源之一(Hart, 2007; Sillitoe, 2010),通常形成于俯沖帶上方的巖漿弧中或俯沖后碰撞帶中(Richards, 2009; Wilkinson, 2013)。已有研究表明,成礦巖體通常就位于3~5km深度(Sillitoe, 2010; 陳華勇和吳超, 2020),礦床平均厚度為2km(Aliyarietal., 2020)。世界范圍內(nèi)保存完好的斑巖型礦床大多形成于中-新生代(Singeretal., 2008),而多數(shù)古生代及其以前的斑巖型礦床常常受到長(zhǎng)期的剝蝕而消失殆盡(Kesler and Wilkinson, 2006; Wilkinson and Kesler, 2007)。然而,還是有少量古老的斑巖型礦床得以保存下來(lái)(如內(nèi)蒙古泥盆紀(jì)奧尤陶勒蓋銅金礦床,新疆天山晚古生代的土屋-延?xùn)|斑巖銅礦床和西準(zhǔn)噶爾晚古生代蘇云河斑巖型鉬礦床等),這些礦床的儲(chǔ)量巨大,具備一定的經(jīng)濟(jì)價(jià)值。前中生代斑巖型礦床經(jīng)過(guò)長(zhǎng)期的構(gòu)造演化而保存良好,必定存在特殊的保存條件,研究古老斑巖型礦床成礦后的剝露和去頂過(guò)程是理解其保存條件的關(guān)鍵。

    近二十年來(lái),低溫?zé)崮甏鷮W(xué)(裂變徑跡和(U-Th)/He定年)分析技術(shù)已經(jīng)逐步發(fā)展成熟,并被廣泛應(yīng)用到礦床地質(zhì)領(lǐng)域(McInnesetal., 2005; 袁萬(wàn)明, 2016; Yinetal., 2019)。磷灰石(AFT)和鋯石(ZFT)裂變徑跡和鋯石 (U-Th)/He(ZHe)方法能有效提取礦床成礦期后的不同時(shí)間和溫度區(qū)間內(nèi)的熱演化信息(Reinersetal., 2004; McInnesetal., 2005; Enkelmann and Garver, 2016)。傳統(tǒng)的鋯石U-Pb同位素體系具有較高的封閉溫度,主要記錄火成巖的侵位時(shí)代,含鉀礦物的40Ar/39Ar測(cè)年(如:角閃石、黑云母和鉀長(zhǎng)石40Ar/39Ar)主要記錄熱液蝕變和成礦時(shí)代,結(jié)合含礦巖體或圍巖中鋯石和磷灰石等礦物的 (U-Th)/He和裂變徑跡測(cè)年,可以獲取金屬礦床礦化、蝕變以及后期剝露過(guò)程完整的時(shí)間-溫度信息。

    本文以新疆西準(zhǔn)噶爾晚古生代的蘇云河斑巖型鉬礦為例,通過(guò)該礦床的鋯石 (U-Th)/He以及鋯石和磷灰石裂變徑跡熱年代學(xué)的研究,結(jié)合前人已發(fā)表的輝鉬礦Re-Os和鋯石U-Pb年齡來(lái)約束該礦床形成和熱液演化過(guò)程,旨在揭示礦床埋藏與剝露過(guò)程,探討其特殊的礦床保存條件。

    1 區(qū)域地質(zhì)背景與礦床地質(zhì)特征

    1.1 區(qū)域地質(zhì)背景

    西準(zhǔn)噶爾地區(qū)位于中國(guó)新疆維吾爾自治區(qū)境內(nèi),處于準(zhǔn)噶爾盆地、天山和阿爾泰山、哈薩克斯坦地體所包圍的一個(gè)三角形區(qū)域(圖1),是巴爾喀什-西準(zhǔn)噶爾成礦帶的重要組成部分(朱永峰, 2014; Shenetal., 2015; 申萍等, 2015)。以NEE向謝米斯臺(tái)斷裂為界,西準(zhǔn)噶爾劃分為北部和南部(Xuetal., 2012)。西準(zhǔn)噶爾北部的謝米斯臺(tái)和薩吾爾兩條斷裂控制著該地區(qū)巖漿巖的分布(Chenetal., 2010)。

    圖1 中亞造山帶構(gòu)造簡(jiǎn)圖(a, 據(jù)Jahn et al., 2000修改)和西準(zhǔn)噶爾地區(qū)地質(zhì)簡(jiǎn)圖(b, 據(jù)Yin et al., 2017修改)TC-塔里木克拉通;NCC-華北克拉通Fig.1 Simplified tectonic divisions of the Central Asian Orogenic Belt (a, modified after Jahn et al., 2000) and simplified geological map of the West Junggar (b, modified after Yin et al., 2017)TC-Tarim Craton; NCC-North China Craton

    西準(zhǔn)噶爾南部NE向的大斷裂(如達(dá)拉布特?cái)嗔?、安齊斷裂、瑪伊勒斷裂和巴爾魯克斷裂等)十分發(fā)育(圖1b),表現(xiàn)為多組、多期次的時(shí)空分布特征,對(duì)區(qū)內(nèi)地層、巖漿巖、構(gòu)造形態(tài)及成礦等起著重要的控制作用。區(qū)域內(nèi)廣泛分布泥盆-石炭系地層,其中巴爾魯克山位于西準(zhǔn)噶爾西南部,主要出露中泥盆統(tǒng)巴爾魯克組、上泥盆統(tǒng)鐵列克組、下石炭統(tǒng)姜巴斯套組等地層。西準(zhǔn)噶爾地區(qū)晚古生代花崗巖廣泛分布,這些花崗巖的形成時(shí)代主要集中在早石炭世(ca. 340~320Ma)和晚石炭世至早二疊世(ca. 310~290Ma)(韓寶福等, 2006; Gengetal., 2009; Chenetal., 2010; 尹繼元等, 2013)。

    晚古生代以來(lái),西準(zhǔn)噶爾地區(qū)主要經(jīng)歷了俯沖增生、洋盆關(guān)閉、走滑和逆沖推覆等多個(gè)階段(Allenetal., 1995)。達(dá)拉布特?cái)嗔言诙B紀(jì)經(jīng)歷了左旋走滑,使克拉瑪依西巖體偏移了100km(Allenetal., 1995; Chouletetal., 2012)。古地磁資料顯示,西準(zhǔn)噶爾相對(duì)于西伯利亞在早二疊世經(jīng)歷了逆時(shí)針旋轉(zhuǎn),而NW向的北天山斷裂則發(fā)育右旋走滑(Wangetal., 2007; Chouletetal., 2011)。準(zhǔn)噶爾盆地地震剖面顯示,西準(zhǔn)噶爾在二疊紀(jì)被推覆到準(zhǔn)噶爾盆地基底巖石之上(Lietal., 2017)。中亞造山帶西段多數(shù)近直立的韌性剪切帶的黑云母40Ar/39Ar年齡在290~245Ma之間,如額爾齊斯斷裂在二疊紀(jì)(290~265Ma)發(fā)生了地殼尺度的逆沖(Briggsetal., 2007),北天山剪切帶北部在245Ma之前發(fā)生了區(qū)域性的右旋走滑事件(Charvetetal., 2007)。這些斷層活動(dòng)表明新疆北部在二疊至三疊紀(jì)期間發(fā)生了廣泛變形。這種變形在整個(gè)中亞地區(qū)主要表現(xiàn)為走滑作用,并伴隨著收縮變形或伸展變形(Gillespieetal., 2017, 2020)。而構(gòu)造地質(zhì)學(xué)、古地磁學(xué)、地震等地質(zhì)資料也顯示,西準(zhǔn)噶爾在中生代發(fā)生構(gòu)造活化(Chouletetal., 2013),主要表現(xiàn)為區(qū)內(nèi)大型斷裂均發(fā)生了右旋走滑作用(余養(yǎng)里, 2016)。有學(xué)者認(rèn)為蒙古-鄂霍茲克海的閉合和歐亞板塊南緣的地塊間碰撞的遠(yuǎn)程效應(yīng)波及中亞西部,包括準(zhǔn)噶爾地區(qū)(Hendrixetal., 1992; Sobel and Arnaud, 1999; Jolivetetal., 2001; Gongetal., 2021),使西準(zhǔn)噶爾地區(qū)在中生代構(gòu)造活動(dòng)非常活躍,并引發(fā)該地區(qū)整體抬升和夷平作用(Dumitruetal., 2001; 李瑋等, 2007; 李麗等, 2008; Lietal., 2014)。

    1.2 礦床地質(zhì)特征

    蘇云河斑巖型鉬礦位于巴爾魯克斷裂北側(cè)(圖1、圖2),礦區(qū)主要出露地層為中泥盆統(tǒng)巴爾魯克組,由凝灰?guī)r和安山質(zhì)凝灰?guī)r等組成(鐘世華等, 2015)。NE向斷裂構(gòu)造在礦區(qū)內(nèi)十分發(fā)育,控制著含礦帶的分布,礦化主要發(fā)生在長(zhǎng)英質(zhì)巖石及其圍巖中(圖2)。礦區(qū)內(nèi)共出露三個(gè)花崗斑巖體,分別為Ⅰ號(hào)、Ⅱ號(hào)和Ⅲ號(hào)巖體,其中Ⅱ號(hào)巖體出露面積最大。三個(gè)巖體均呈巖株?duì)钋秩氲街心嗯杞y(tǒng)地層中。礦化主要發(fā)生在斑巖體及其接觸的圍巖中,礦體呈水平-緩傾的巨厚層狀或似層狀,主要賦存在中泥盆統(tǒng)巴爾魯克組火山沉積巖中,部分賦存于花崗斑巖及花崗閃長(zhǎng)斑巖體中(圖2)。礦體與圍巖無(wú)明顯界限,呈漸變接觸關(guān)系,目前已探明的鉬儲(chǔ)量為57萬(wàn)t。

    圖2 蘇云河斑巖型鉬礦礦區(qū)地質(zhì)圖(據(jù)Shen et al., 2017修改)Fig.2 Geological map of Suyunhe porphyry molybdenum deposit (modified after Shen et al., 2017)

    礦石礦物包括輝鉬礦、白鎢礦、黃鐵礦和黃銅礦,以及少量的輝銅礦、斑銅礦、鈦鐵礦和磁黃鐵礦。脈石礦物以石英、鉀長(zhǎng)石、斜長(zhǎng)石和黑云母為主。鉬礦化主要出現(xiàn)在石英脈和細(xì)脈中(Shenetal., 2017)。鉬礦化相關(guān)的圍巖蝕變包括鉀質(zhì)、綠泥石-白云母和千枚巖蝕變。在熱液活動(dòng)的初始階段,鉀質(zhì)蝕變廣泛發(fā)育,后被綠泥石-白云母蝕變和火山碎屑圍巖中的千枚巖蝕變疊加。綠泥石-白云母蝕變具有綠泥石+白云母+伊利石+石英+鈦石的礦物組合,以黑云母部分至完全轉(zhuǎn)變?yōu)榫G泥石、斜長(zhǎng)石轉(zhuǎn)變?yōu)榧?xì)粒白云母和伊利石為代表。千枚巖蝕變具有白云母+伊利石+石英+黃鐵礦的礦物組合,廣泛存在于火山碎屑圍巖中(Shenetal., 2017)。鐘世華等(2015)將蘇云河鉬礦成礦過(guò)程劃分為3個(gè)階段。早期流體以巖漿水為主,發(fā)育石英-鉀長(zhǎng)石脈、石英-鉀長(zhǎng)石-(黃鐵礦)-輝鉬礦脈、鉀長(zhǎng)石脈和無(wú)礦石英脈;中期流體以巖漿水和大氣水的混合水為主,主要發(fā)育石英-輝鉬礦-黃鐵礦脈、石英-輝鉬礦脈和石英-黃鐵礦脈和無(wú)礦石英脈;晚期流體的混合水中,大氣水的比例增加,為中低溫、低鹽度流體,主要發(fā)育石英-碳酸鹽脈、石英-黃鐵礦脈和無(wú)礦石英脈。中期為主要成礦階段,而晚期脈體中幾乎沒(méi)有礦化。

    楊猛等(2015)和Shenetal.(2017)通過(guò)測(cè)定花崗巖樣品的鋯石U-Pb年齡限定了侵入巖的活動(dòng)年限為294~309Ma。Shenetal.(2013)和鐘世華等(2015)報(bào)道蘇云河鉬礦中輝鉬礦的Re-Os加權(quán)年齡分別為301±4Ma和295±2Ma,表明礦化發(fā)生于晚石炭世-早二疊世。對(duì)含礦斑巖及圍巖中的流體包裹體進(jìn)行成礦條件分析顯示,主要成礦階段流體均一溫度為240~443℃,壓力77.3~81.3MPa,對(duì)應(yīng)成礦深度為1.9~4.8km(楊志強(qiáng)等, 2012; 鐘世華等, 2015)。

    2 樣品與測(cè)試方法

    2.1 樣品描述

    本次研究采集了4件花崗質(zhì)巖石樣品開(kāi)展低溫?zé)崮甏鷮W(xué)研究,采樣位置見(jiàn)圖2,分別為石英閃長(zhǎng)玢巖(WJ16-02)、微細(xì)粒斑狀二長(zhǎng)花崗巖(WJ16-06)、二長(zhǎng)閃長(zhǎng)玢巖(WJ16-07)和二長(zhǎng)花崗巖(WJ17-11)。所有樣品均為斑狀或似斑狀結(jié)構(gòu),斑晶含量為35%~45%,斑晶礦物主要為斜長(zhǎng)石(20%~30%)、鉀長(zhǎng)石(少量)、石英(約10%)、角閃石(5%~10%)和黑云母(少量)。基質(zhì)主要由斜長(zhǎng)石組成(20%~50%),含少量鉀長(zhǎng)石、角閃石、石英、黑云母。樣品中的黑云母等暗色礦物普遍被白云母及綠泥石交代呈假象(圖3)。4件樣品均有不同程度礦化,其中樣品WJ16-07較為明顯,含有大量立方體晶形的黃鐵礦、輝鉬礦顆粒,部分白云母被碳酸鹽巖礦物交代,有輕微綠泥石化。樣品WJ16-06和WJ17-11中穿插有大量石英細(xì)脈,綠泥石化蝕變較強(qiáng)。樣品的基本情況見(jiàn)表1。

    表1 蘇云河鉬礦樣品位置和熱年代學(xué)分析方法匯總Table 1 Summary of sample localities and thermochronological analysis in the Suyunhe Mo deposit

    圖3 蘇云河地區(qū)花崗質(zhì)巖石的正交偏光鏡下照片(a、b)石英閃長(zhǎng)玢巖(WJ16-02);(c、d)二長(zhǎng)花崗巖(WJ16-06);(e、f)二長(zhǎng)閃長(zhǎng)玢巖(WJ16-07). Ksp-鉀長(zhǎng)石;Pl-斜長(zhǎng)石;Bi-黑云母;Chl-綠泥石;Hbl-角閃石;Mus-白云母;Cbn-碳酸鹽Fig.3 Microscopic photos of granitic rocks in the Suyunhe area (under CPL)(a, b) quartz diorite porphyrite (WJ16-02); (c, d) monzodiorite porphyrite (WJ16-06); (e, f) monzogranite (WJ16-07). Ksp-K-feldspar; Pl-plagioclase; Bi-biotite; Chl-chlorite; Hbl-hornblende; Mus-muscovite; Cbn-carbonate

    2.2 樣品測(cè)試方法

    2.2.1 鋯石和磷灰石裂變徑跡測(cè)年

    裂變徑跡分析是在美國(guó)亞利桑那大學(xué)地球科學(xué)系完成,鋯石和磷灰石顆粒通過(guò)外探測(cè)器方法進(jìn)行測(cè)年。磷灰石顆粒被安裝在環(huán)氧樹(shù)脂、 氧化鋁和金剛石中拋光,然后在20±1℃下用5.5M的HNO3蝕刻20±0.5秒,顯示出自發(fā)裂變徑跡。鋯石顆粒固定在全氟烷氧基樹(shù)脂上,使用金剛石進(jìn)行拋光,而后在鋯坩堝 220℃條件下使用KOH-NaOH共晶熔體蝕刻3~36小時(shí)(Gleadowetal., 1986)。最佳刻蝕時(shí)間取決于樣品的年齡和輻射損傷,需要通過(guò)重復(fù)刻蝕和每隔3~6小時(shí)觀察監(jiān)測(cè)。在美國(guó)俄勒岡州立大學(xué)Triga反應(yīng)堆進(jìn)行照射,中子通量使用歐洲標(biāo)準(zhǔn)物質(zhì)和測(cè)量研究所(IRMM)的鈾玻璃IRMM 540R進(jìn)行監(jiān)測(cè)。照射后通過(guò)用48% HF刻蝕20分鐘揭示出誘發(fā)徑跡。最后利用奧林巴斯BX61顯微鏡和Kinetek自動(dòng)分級(jí)系統(tǒng)計(jì)算自發(fā)和誘發(fā)的裂變徑跡密度。裂變徑跡的長(zhǎng)度和Dpar值使用FTStage軟件進(jìn)行校準(zhǔn),裂變徑跡年齡采用Hurford and Green (1983)推薦的Zeta校準(zhǔn)方法IUGS計(jì)算,以獲得±1σ誤差的中心年齡。根據(jù)Hurford (1990)的建議,通過(guò)反復(fù)校準(zhǔn)包括Durango、Fish Canyon磷灰石和Fish Canyon鋯石在內(nèi)的大量國(guó)際公認(rèn)的年齡標(biāo)準(zhǔn),獲得了磷灰石和鋯石IRMM 540R和IRMM541 Zeta校準(zhǔn)因子分別為351.9±3.8和116.0±1.3。詳細(xì)分析流程依據(jù)Thomson and Ring (2006)的標(biāo)準(zhǔn)。

    2.2.2 鋯石(U-Th)/He測(cè)年

    鋯石(U-Th)/He測(cè)年在美國(guó)科羅拉多大學(xué)波爾得分校熱年代學(xué)研究和儀器實(shí)驗(yàn)室完成。鋯石顆粒首先在顯微鏡下進(jìn)行挑選,主要觀察鋯石顆粒的尺寸,形狀和包裹體等情況。由于α粒子的射出效應(yīng)的存在,需要對(duì)實(shí)測(cè)年齡進(jìn)行校正。將挑選好的鋯石顆粒進(jìn)行拍照,并在鏡下測(cè)量鋯石顆粒的長(zhǎng)、寬和高及錐體的長(zhǎng)度,根據(jù)測(cè)量結(jié)果計(jì)算α校正系數(shù)。完成測(cè)量后鋯石顆粒裝入鈮囊中。用四極質(zhì)譜儀來(lái)測(cè)定4He/3He值,4He含量用同位素稀釋法測(cè)定。完成4He含量分析后,鋯石顆粒需要經(jīng)歷多次消解。第一次消解時(shí)樣品中加入200μL含235U-230Th-145Nd稀釋劑的HF酸。將小瓶加蓋,堆放在125mL的聚四氟乙烯瓶中,然后放置在高壓釜(Parr)中,在220℃下消解72小時(shí)。消解過(guò)的溶液在90℃的熱盤上蒸干后,加入200μL 7:1 HNO3:HF的混合液,然后加蓋,在90℃的熱盤消減4個(gè)小時(shí)。礦物溶解后用1~3mL的雙去離子水稀釋,使用Agilent 7900型四極桿電感耦合等離子體質(zhì)譜儀(Q-ICP-M_S)進(jìn)行U、Th和Sm含量分析。詳細(xì)的分析流程請(qǐng)參考Johnsonetal.(2017)。

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

    3.1 裂變徑跡測(cè)年結(jié)果

    3.1.1 磷灰石裂變徑跡測(cè)年結(jié)果

    本文采集的4件樣品均進(jìn)行了磷灰石裂變徑跡分析。在每個(gè)樣品中挑選出20個(gè)質(zhì)量良好的磷灰石顆粒進(jìn)行測(cè)試,測(cè)試結(jié)果如表2和圖4所示。4件樣品的平均裂變徑跡長(zhǎng)度為13.4~14.0μm,均顯示為較窄的單峰分布(圖4),暗示樣品快速通過(guò)部分退火帶(Gleadowetal., 1986; 呂紅華等, 2013)。所有樣品均通過(guò)卡方檢驗(yàn)(P(χ2)>5%,表2),顯示了單一的單顆粒年齡,即樣品的年齡可以反映其經(jīng)歷的構(gòu)造熱事件。樣品的中心年齡分別為190±13.8Ma(樣品WJ16-02)、126±7.3Ma(樣品WJ16-06)、186±13.9Ma(樣品WJ16-07)和146±6.1Ma(樣品WJ17-11)。4件樣品的平均Dpar值在2.02~2.45μm之間,大于標(biāo)準(zhǔn)礦物磷灰石的蝕刻半徑(Donlicketal., 1999, 2005),表明其較高的耐退火率。

    圖4 蘇云河斑巖鉬礦區(qū)磷灰石裂變徑跡年齡雷達(dá)圖(左)和徑跡長(zhǎng)度分布圖(右)Fig.4 Radial plots (left) and track length distribution (right) of apatite fission track of the Suyunhe Mo deposit

    3.1.2 鋯石裂變徑跡測(cè)年結(jié)果

    樣品WJ17-11進(jìn)行了鋯石裂變徑跡分析。在樣品中挑選14個(gè)質(zhì)量良好的鋯石顆粒進(jìn)行測(cè)試,所有單顆粒年齡分布在193~317Ma之間,且通過(guò)卡方檢驗(yàn)(P(χ2)>5%),說(shuō)明單顆粒年齡聚合性較好。中心年齡采用IUGS推薦的Zeta校準(zhǔn)方法計(jì)算(Hurford and Green, 1983),本次實(shí)驗(yàn)鋯石的Zeta值為116.0±1.3,其中心年齡為247±15.8Ma(表2)。

    3.2 鋯石(U-Th)/He測(cè)年結(jié)果

    樣品WJ16-06的鋯石(U-Th)/He分析結(jié)果見(jiàn)表3。4個(gè)鋯石顆粒的(U-Th)/He年齡分別為209±1Ma、208±2Ma、238±2Ma、239±2Ma(表3),其加權(quán)平均年齡為224±1Ma。該樣品的單顆粒ZHe年齡或者加權(quán)年齡都明顯比其鋯石U-Pb年齡更年輕,表明這些數(shù)據(jù)是可靠的,能反映巖體就位以來(lái)的冷卻過(guò)程。

    表3 蘇云河鉬礦的鋯石(U-Th)/He結(jié)果Table 3 Zircon (U-Th)/He results from the Suyunhe Mo deposit

    3.3 HeFTy軟件模擬

    本文使用HeFTy軟件對(duì)樣品進(jìn)行熱歷史模擬,以便更好地了解樣品所經(jīng)歷的時(shí)間-溫度信息。輸入數(shù)據(jù)包括磷灰石裂變徑跡的單顆粒年齡、長(zhǎng)度、蝕刻半徑以及C軸投影夾角等。擴(kuò)散動(dòng)力學(xué)模型引自Ketchametal.(2007)的多元?jiǎng)恿W(xué)模型,動(dòng)力學(xué)參數(shù)用Dpar值,初始徑跡長(zhǎng)度為16.3μm(Donelicketal., 1999)。長(zhǎng)度擬合最優(yōu)方程選用Kolmogorov-Smirnov Test,計(jì)算方法選擇Monte Carlo(Ketchametal., 2009)。起始點(diǎn)限制條件為溫度280~320℃,年齡260~280Ma(圖5),溫度的限制必須高于ZFT封閉溫度,年齡的限制范圍是ZFT和鋯石U-Pb的年齡之間。同時(shí)假定模擬終點(diǎn)的地表溫度為20±3℃。其中模擬結(jié)果分為良好熱歷史曲線(GOF>0.5)和較好熱歷史曲線(GOF>0.05)。每件樣品的模擬直到良好路徑達(dá)到100條時(shí)才會(huì)停止,模擬結(jié)果見(jiàn)圖5。其中,最佳擬合路徑和加權(quán)平均路徑分別用黑色和藍(lán)色曲線標(biāo)注。本研究熱史模擬年齡的GOF值均大于0.9,徑跡長(zhǎng)度GOF值大于0.6,說(shuō)明本次熱模擬結(jié)果可信。AFT、ZFT、ZHe的封閉溫度區(qū)間分別為60~110℃、210~270℃、160~200℃(Gleadowetal., 1986; Green, 1989; Corrigan, 1991)。由于裂變徑跡在高于封閉溫度時(shí)退火作用幾乎在瞬時(shí)發(fā)生,而低于封閉溫度時(shí)退火速率極低(Gleadowetal., 1986),因此,退火模型僅在封閉溫度(PAZ)范圍內(nèi)有效,不應(yīng)過(guò)度解釋PAZ以外獲得的冷卻曲線。本文的熱歷史模擬結(jié)果顯示蘇云河斑巖鉬礦經(jīng)歷了一期快速冷卻,即中三疊世至早白堊世(240~120Ma)。

    圖5 蘇云河鉬礦花崗質(zhì)巖石樣品的熱歷史模擬結(jié)果紫色區(qū)域代表擬合良好;綠色區(qū)域代表擬合可接受;藍(lán)色線代表加權(quán)平均路徑;黑色線代表最佳擬合曲線;方框代表模擬初始值限制;4個(gè)樣品均輸入AFT年齡作為初始條件;GOF代表擬合良好值;TL代表徑跡長(zhǎng)度.封閉溫度范圍為:AFT:60~110℃,ZHe: 160~200℃,ZFT:210~270℃Fig.5 Results of HeFTy thermal history modeling of selected granitic samples from the Suyunhe porphyry Mo deposit

    3.4 剝露速率計(jì)算

    本文使用礦物對(duì)法、年齡-封閉溫度法和熱史模擬法計(jì)算蘇云河鉬礦的剝露速率(丁汝鑫等, 2007; 常遠(yuǎn)和周祖翼, 2010; 呂紅華等, 2013)。礦物對(duì)法是指對(duì)同一樣品采用不同方法進(jìn)行測(cè)年,而后根據(jù)封閉溫度差和年齡差的比來(lái)計(jì)算剝露速率。年齡-封閉溫度法為已知樣品年齡和封閉溫度,假定區(qū)域的地溫梯度和當(dāng)今地表溫度來(lái)計(jì)算剝露速率。而熱史模擬法是根據(jù)模擬曲線及冷卻區(qū)間的劃分來(lái)計(jì)算剝露速率。因此,熱史模擬法較之其他兩種方法,更能體現(xiàn)每個(gè)階段的剝露速率和揭頂量的差異。由于AFT、ZHe和ZFT的封閉溫度分別為60~110℃、160~200℃和210~270℃(Hurford and Green, 1983; Ketchametal., 1999; Reinersetal., 2004; Enkelmann and Garver, 2016),計(jì)算時(shí),AFT、ZHe與ZFT的封閉溫度分別取85℃、180℃和240℃,計(jì)算結(jié)果見(jiàn)表4。根據(jù)熱史模擬法的計(jì)算結(jié)果,蘇云河鉬礦在中二疊世-早白堊世(240~120Ma)的剝露速率為49.0~56.7m/Myr,在早白堊世至今(120~0Ma)的剝露速率為6.7~21.7m/Myr。

    表4 基于年齡-封閉溫度法、礦物對(duì)法和熱史模擬法估算蘇云河斑巖型鉬礦不同時(shí)期的剝露速率Table 4 Denudation rate of different periods from the Suyunhe porphyry Mo deposit based on age-closure temperature, mineral pairs and thermal modeling method

    4 討論

    4.1 蘇云河斑巖型鉬礦的巖漿-熱液演化過(guò)程

    斑巖型礦床主要形成于匯聚板片邊緣,是復(fù)雜巖漿和熱液過(guò)程的產(chǎn)物,其就位古深度通常小于6km(Wilkinson and Kesler, 2007; Richards, 2009; Sillitoe, 2010; Cookeetal., 2014)。準(zhǔn)確掌握巖漿-熱液-冷卻過(guò)程是探索斑巖礦床形成、演化和保存機(jī)制的基礎(chǔ) (McInnesetal., 2005; Wilkinson and Kesler, 2009)。傳統(tǒng)鋯石U-Pb同位素體系具有較高的封閉溫度,主要記錄火成巖的侵位時(shí)代(Cherniak and Watson, 2001),輝鉬礦Re-Os年齡,鋯石裂變徑跡和(U-Th)/He年齡主要記錄成礦和熱液活動(dòng)時(shí)代(Green, 1989; Suzukietal., 1996)。楊猛等(2015)報(bào)道了蘇云河鉬礦區(qū)I#含礦花崗巖體的LA-ICP-MS鋯石U-Pb年齡為309±4Ma。而Shenetal.(2017)得到蘇云河花崗巖和花崗閃長(zhǎng)巖的SIMS鋯石U-Pb年齡為298~294Ma。這些鋯石U-Pb年齡結(jié)果表明,蘇云河鉬礦區(qū)侵入巖巖漿活動(dòng)時(shí)限為晚石炭-早二疊世。Shenetal.(2013)和鐘世華等(2015)報(bào)道蘇云河鉬礦中輝鉬礦的Re-Os加權(quán)年齡分別為301±4Ma和295±2Ma,表明礦化發(fā)生于晚石炭世晚期到早二疊世。蘇云河鉬礦區(qū)花崗巖的侵位時(shí)限與成礦年齡在誤差范圍內(nèi)一致,表明這些花崗巖可能與鉬礦化密切相關(guān)。另外,蘇云河礦床綠泥石化蝕變廣泛發(fā)育。三種類型的綠泥石化蝕變被識(shí)別出來(lái),分別為綠泥石化蝕變Ⅰ、Ⅱ和Ⅲ型。Ⅰ型綠泥石賦存于深部花崗巖中,以原生黑云母或角閃石的置換為特征。除綠泥石蝕變外,絹云母蝕變常被觀察到,并取代花崗巖中的斜長(zhǎng)石或長(zhǎng)石;Ⅱ型綠泥石主要賦存于含礦石英脈或細(xì)脈中,Ⅱ型綠泥石與輝鉬礦有密切的關(guān)系,因此它們的化學(xué)成分可以用來(lái)測(cè)定鉬礦化的物理化學(xué)條件;Ⅲ型綠泥石也伴有絹云母。但它們賦存于淺層火山沉積序列中,并形成蝕變邊,與貧瘠的石英脈或細(xì)脈接壤(鐘世華等, 2015; Caoetal., 2020)。類型Ⅱ的綠泥石化蝕變是相對(duì)高的溫度(319~353℃),與所測(cè)量的流體包裹體溫度一致(鐘世華等, 2015; Caoetal., 2020)。而類型Ⅰ和Ⅲ綠泥石化蝕變的溫度較低(247~319℃)。鋯石裂變徑跡的封閉溫度剛好可以控制晚期綠泥石化蝕變的上限,這表明該礦床的低溫蝕變不晚于240Ma。因此,蘇云河斑巖型礦床可能經(jīng)歷了長(zhǎng)期的(>55Myr)蝕變過(guò)程。

    4.2 蘇云河斑巖型鉬礦的礦床保存和剝露史

    通常,巖漿從巖漿房上侵到淺部地殼過(guò)程中,由于巖漿與圍巖巨大的溫度差,導(dǎo)致巖漿與圍巖之間快速的熱傳導(dǎo)過(guò)程而快速降溫,其冷卻速率可達(dá)1000℃/Myr,甚至更高(McInnesetal., 2005; Harrisetal., 2008; Fuetal., 2010; Lietal., 2014)。蘇云河花崗巖和花崗閃長(zhǎng)巖從鋯石U-Pb(298~294Ma)到ZFT(247±16Ma)封閉溫度的平均冷卻速率僅為10.1~19.7℃/Myr,比典型的斑巖型礦床小兩個(gè)數(shù)量級(jí)。這可能有以下三種原因:(1)花崗巖就位較深,圍巖溫度大于ZFT部分保留區(qū),導(dǎo)致ZFT系統(tǒng)長(zhǎng)時(shí)間沒(méi)有計(jì)時(shí);(2)花崗巖的ZFT年齡被成礦后巖漿熱液事件所重置;(3)成礦后接受了較厚的埋藏,礦體被埋藏至較深的位置,埋藏溫度高于ZFT的部分保留區(qū),重置了ZFT年齡。根據(jù)邵潔璉計(jì)算成礦壓力和深度的經(jīng)驗(yàn)公式(楊志強(qiáng)等, 2012),獲得蘇云河斑巖鉬礦成礦深度為1.9~4.8km(鐘世華等, 2015)。假定地表溫度為20℃,地溫梯度為25℃/km。蘇云河斑巖鉬礦最大侵位深度對(duì)應(yīng)的圍巖溫度為140℃左右?;◢弾r的形成溫度明顯高于圍巖溫度。因此,這些花崗巖只是淺部就位,第一種可能排除。蘇云河地區(qū)已報(bào)道的最年輕的鋯石U-Pb年齡為293.7±2.3Ma(Shenetal., 2017),在早二疊世之后沒(méi)有巖漿活動(dòng)的證據(jù),而且也沒(méi)有受到后期流體交代的跡象,第二種可能也排除。西準(zhǔn)噶爾在晚石炭至早二疊世主要為伸展背景(Gengetal., 2009; Yinetal., 2010; 尹繼元等, 2013),十分有利于接受沉積。由于斷裂構(gòu)造作用,西準(zhǔn)噶爾在~285~260Ma經(jīng)歷了普遍剝蝕過(guò)程(Gillespieetal., 2020),這可能為礦區(qū)之上的蓋層提供了物源。此外,位于準(zhǔn)噶爾盆地西南端的石炭紀(jì)碎屑巖樣品的熱歷史模擬也顯示,它們經(jīng)歷了二疊紀(jì)至中三疊世的埋藏過(guò)程,該樣品至少被加熱至~120℃,導(dǎo)致其AFT年齡重置(Glorieetal., 2019)。因此我們認(rèn)為,蘇云河鉬礦可能也經(jīng)歷了一期快速埋藏過(guò)程,其溫度超過(guò)ZFT的部分保留區(qū),導(dǎo)致其年齡重置。

    以地溫梯度25℃/km和ZFT封閉溫度210~270℃計(jì)算(Hurford and Green, 1983; Brandonetal., 1998),若要重置ZFT年齡,樣品的溫度至少到達(dá)270℃,對(duì)應(yīng)最大的埋藏深度為10km(圖6a, b)。根據(jù)WJ16-06和WJ17-11兩組樣品的礦物對(duì)計(jì)算結(jié)果,蘇云河鉬礦208~126Ma的冷卻速率為1.1℃/Myr,247~146Ma的冷卻速率為1.5℃/Myr,明顯低于晚石炭-中三疊世期間的冷卻速率,說(shuō)明中三疊世以來(lái)礦體和圍巖的溫度已經(jīng)基本平衡。因此,厚埋藏應(yīng)該發(fā)生于中三疊世之前。斑巖的侵位古深度為1.9~4.8km,需要在晚石炭至中三疊世期間覆蓋約5.2~8.1km的蓋層。

    圖6 蘇云河斑巖型鉬礦侵位、埋藏和剝露模型圖縱坐標(biāo)是以294Ma時(shí)古地表海拔為0km設(shè)置的參考系.(a)花崗巖質(zhì)巖漿侵入到中泥盆統(tǒng)地層中,巖漿作用和成礦作用幾乎同時(shí)形成;(b)礦床在早二疊世-中三疊世之間接受了較厚的沉積;(c)中三疊世-早白堊世期間經(jīng)歷了快速的快速剝露過(guò)程;(d)早白堊世以后經(jīng)歷緩慢剝露過(guò)程Fig.6 Schematic diagrams of emplacement, burial and exhumation of the Suyunhe porphyry Mo deposit

    自中三疊世開(kāi)始,礦床進(jìn)入快速冷卻階段。對(duì)準(zhǔn)噶爾盆地西北緣的物源研究表明,在這一階段西準(zhǔn)噶爾為盆地提供了大量的碎屑沉積(蔚遠(yuǎn)江等, 2020)。熱史模擬結(jié)果和沉積學(xué)證據(jù)均表明,礦床在三疊紀(jì)至早白堊世發(fā)生了快速冷卻。根據(jù)熱史模擬法的計(jì)算結(jié)果,蘇云河鉬礦中三疊世-早白堊世(240~120Ma)的剝露速率為49.0~56.7m/Myr,揭頂量為7.4~9.2km(圖6c)。

    早白堊世之后,礦床進(jìn)入緩慢冷卻階段,蘇云河鉬礦所受剝蝕作用明顯減緩。根據(jù)熱史模擬法的計(jì)算結(jié)果,其平均剝露速率為6.7~21.7m/Myr,較之中三疊世-早白堊世已經(jīng)明顯減緩。在緩慢冷卻階段,該地區(qū)的揭頂量為0.8~2.6km。此外,準(zhǔn)噶爾盆地的沉積學(xué)相關(guān)研究表明,白堊紀(jì)末期以來(lái),準(zhǔn)噶爾盆地受全球氣候變冷和青藏高原隆升共同影響,干旱化程度不斷加強(qiáng)(王熠哲等, 2019),這種干旱性氣候在中生代晚期及新生代一定程度上限制了西準(zhǔn)噶爾的剝蝕作用(Pullenetal., 2020)。這些證據(jù)表明蘇云河鉬礦在晚白堊世之后經(jīng)歷了緩慢剝露過(guò)程,直至裸露至地表(圖6d)。

    4.3 西準(zhǔn)噶爾中生代快速冷卻事件

    目前,一些學(xué)者認(rèn)為西準(zhǔn)噶爾地區(qū)中生代經(jīng)歷快速冷卻過(guò)程(Yinetal., 2018),也有學(xué)者認(rèn)為是緩慢冷卻過(guò)程(Gillespieetal., 2020)。本文獲得的ZFT、ZHe和AFT年齡結(jié)果表明西準(zhǔn)噶爾地區(qū)經(jīng)歷了漫長(zhǎng)的中生代冷卻過(guò)程(247~126Ma)。熱模擬也顯示該地區(qū)經(jīng)歷了一期連續(xù)的冷卻和剝露過(guò)程(圖7)。樣品WJ16-06和WJ17-11使用礦物對(duì)法計(jì)算結(jié)果表明,其在224~126Ma和246~145Ma期間的冷卻速率分別為1.1℃/Myr和1.5℃/Myr,指示了相對(duì)快速的冷卻。這些結(jié)果與準(zhǔn)噶爾盆地在中生代接受了巨厚的沉積一致。比如:準(zhǔn)噶爾盆地西北緣的紅山嘴地區(qū)在三疊紀(jì)持續(xù)拗陷并發(fā)生沉積補(bǔ)償,湖盆內(nèi)充填巨厚的三疊-白堊系沉積(高陽(yáng)等, 2016)。這些中生代沉積巖物源來(lái)自于西北地區(qū),由厚層的礫巖和沉積角礫巖組成,巖石成分成熟度差、碎屑顆粒分選差,均支持這一期的快速剝露(Hendrixetal., 1992)。另外,前期在西準(zhǔn)噶爾包古圖地區(qū)也獲得相似的中生代的ZHe和AHe熱年代學(xué)數(shù)據(jù)(Yinetal., 2018)。這表明這一時(shí)期的快速冷卻不僅僅局限于蘇云河地區(qū),而是整個(gè)西準(zhǔn)噶爾地區(qū)。前人的研究顯示,西準(zhǔn)噶爾地區(qū)中生代經(jīng)歷多期次的快速冷卻事件,可能與區(qū)域南部的遠(yuǎn)程效應(yīng)有關(guān)(如:羌塘和昆侖-柴達(dá)木碰撞或者羌塘和拉薩碰撞)(Hendrixetal., 1992; Dumitruetal., 2001)。然而,西準(zhǔn)噶爾在二疊-三疊紀(jì)呈現(xiàn)由NW向SE的逆沖、旋轉(zhuǎn)的區(qū)域構(gòu)造(何登發(fā)等, 2018),這些構(gòu)造活動(dòng)并非受到由南向北的擠壓應(yīng)力而形成的。本文的研究結(jié)果表明,單一的構(gòu)造事件的遠(yuǎn)程效應(yīng)很難解釋如此漫長(zhǎng)的冷卻過(guò)程。相比而言,區(qū)域內(nèi)走滑-逆沖斷層的活動(dòng)可能更好地解釋這期冷卻過(guò)程。

    圖7 蘇云河鉬礦綜合的熱史模擬匯編黑色圓點(diǎn)表示樣品處于封閉溫度區(qū)間中間時(shí)的年齡Fig.7 Interpreted thermal history of the Suyunhe Mo deposit derived from the modeling

    首先,本文在蘇云河地區(qū)獲得的磷灰石裂變徑跡數(shù)據(jù)(190~126Ma)和西準(zhǔn)噶爾南部地區(qū)獲得的磷灰石(U-Th)/He數(shù)據(jù)(237~120Ma)都比較分散(Yinetal., 2018)。Gillespieetal.(2020)在西準(zhǔn)噶爾地區(qū)也獲得比較分散的磷灰石裂變徑跡結(jié)果(307~177Ma)(圖8)。而熱模擬結(jié)果也顯示西準(zhǔn)噶爾地區(qū)經(jīng)歷漫長(zhǎng)的中生代快速冷卻過(guò)程(圖5、圖7)。其次,西準(zhǔn)噶爾地區(qū)在晚古生代末期到中生代由廣泛發(fā)育走滑斷層轉(zhuǎn)換為走滑逆沖斷層(Allenetal., 1995; Allen and Vincent, 1997; Laurent-Charvetetal., 2003; Chouletetal., 2012),并發(fā)育多組逆沖推覆體,由西向東依次為塔爾巴哈臺(tái)-薩吾爾推覆體、巴爾魯克-謝米斯臺(tái)推覆體、加依爾推覆體,呈前展式活動(dòng),紅車斷裂帶、克百斷裂帶和烏夏斷裂帶處于沖斷前鋒部位,同時(shí)伴有地塊旋轉(zhuǎn)(何登發(fā)等, 2018)。第三,準(zhǔn)噶爾盆地西北緣在三疊紀(jì)發(fā)生了強(qiáng)烈擠壓變形,侏羅-白堊紀(jì)雖然構(gòu)造作用有所減弱,但地層內(nèi)部存在多期局部不整合,說(shuō)明局部的斷層活動(dòng)仍然活躍(隋風(fēng)貴, 2015; 胡小文等, 2020)。最后,Yangetal.(2012)通過(guò)地震剖面分析顯示,準(zhǔn)噶爾盆地南隆起的斷裂活動(dòng)時(shí)限一直持續(xù)到白堊紀(jì)。西準(zhǔn)噶爾東南緣斷裂帶在三疊紀(jì)-侏羅紀(jì)發(fā)育右旋走滑活動(dòng)(余養(yǎng)里, 2016)。綜上所述,西準(zhǔn)噶爾地區(qū)在中生代斷裂活動(dòng)頻繁,在時(shí)間上很難與某個(gè)特定的構(gòu)造事件或者遠(yuǎn)程效應(yīng)相關(guān)。因此,我們認(rèn)為西準(zhǔn)噶爾地區(qū)中生代的快速冷卻過(guò)程可能與區(qū)內(nèi)發(fā)生的走滑-逆沖作用有關(guān)。

    圖8 西準(zhǔn)噶爾不同類型低溫年代學(xué)年齡統(tǒng)計(jì)圖圖中數(shù)據(jù)來(lái)自本研究、Yin et al. (2018)和Gillespie et al. (2020). 與圖中標(biāo)示出的冷卻區(qū)間相關(guān)的主要構(gòu)造事件有:古亞洲洋的閉合(310~250Ma)(Gillespie et al., 2020; Glorie et al., 2019);區(qū)域內(nèi)走滑-逆沖斷層活化(240~120Ma)Fig.8 Low temperature geochronological data summary of different dating systems in the West Junggar

    5 結(jié)論

    對(duì)蘇云河鉬礦花崗質(zhì)巖石進(jìn)行了低溫?zé)崮甏鷮W(xué)分析,包括ZFT、AFT和ZHe測(cè)年。綜合其它地質(zhì)背景資料和年齡數(shù)據(jù),得出以下結(jié)論:

    (1)蘇云河礦床經(jīng)歷了超過(guò)55Myr的中-低溫蝕變過(guò)程。

    (2)蘇云河斑巖型鉬礦在晚石炭-早三疊世期間接受大約5.2~8.1km的蓋層埋藏。而后經(jīng)歷了中生代快速冷卻過(guò)程,揭頂量約為7.4~9.2km,在晚白堊世之后進(jìn)入緩慢冷卻階段,揭頂量為0.8~2.6km。礦床的保存很大程度上受到后期的埋藏和抬升速率的控制。

    (3)西準(zhǔn)噶爾地區(qū)中生代的快速冷卻可能是由區(qū)內(nèi)大型斷層活化引起的。

    致謝韓春明研究員與劉麗萍博士與作者就斑巖成礦與保存機(jī)制進(jìn)行了討論;三位匿名審稿人詳細(xì)審閱了全文,并提出很多寶貴的意見(jiàn)和建議;在此一并衷心的表示感謝!

    謹(jǐn)以此文致敬李繼亮先生!李繼亮先生在多個(gè)領(lǐng)域的開(kāi)創(chuàng)性研究,引領(lǐng)和推進(jìn)了我國(guó)沉積學(xué)和造山帶大地構(gòu)造學(xué)的研究發(fā)展。他嚴(yán)謹(jǐn)?shù)目蒲袘B(tài)度、崇高的學(xué)術(shù)思想令我輩深深折服,感謝先生在中國(guó)地質(zhì)學(xué)科作出的杰出貢獻(xiàn)!

    猜你喜歡
    徑跡準(zhǔn)噶爾鉬礦
    專利名稱:一種鉬礦和鎢礦的聯(lián)合冶煉工藝
    基于蒙特卡羅模擬方法的圓筒形固體核徑跡氡探測(cè)器探測(cè)效率的研究
    四、準(zhǔn)噶爾對(duì)天山南北麓維吾爾人的統(tǒng)治
    二、準(zhǔn)噶爾汗國(guó)的建立與葉爾羌汗國(guó)
    裂變徑跡LA-ICP-MS/FT法原理、實(shí)驗(yàn)流程和應(yīng)用
    全球十大鉬礦
    基于波形特征的露天鉬礦微震事件的識(shí)別分析——以卓資山鉬礦為例
    乾隆年間清軍的第三次出征準(zhǔn)噶爾考辨
    準(zhǔn)噶爾遠(yuǎn)征軍進(jìn)藏線路考
    核乳膠中質(zhì)子徑跡計(jì)算機(jī)輔助識(shí)別方法研究
    高清欧美精品videossex| 男女高潮啪啪啪动态图| 久久精品久久精品一区二区三区| 熟女少妇亚洲综合色aaa.| 亚洲欧美成人综合另类久久久| 婷婷色麻豆天堂久久| av网站免费在线观看视频| 51午夜福利影视在线观看| 肉色欧美久久久久久久蜜桃| 51午夜福利影视在线观看| 狠狠婷婷综合久久久久久88av| 免费女性裸体啪啪无遮挡网站| 永久免费av网站大全| 亚洲av中文av极速乱| 婷婷色综合大香蕉| 一本—道久久a久久精品蜜桃钙片| 国产精品成人在线| 国产精品秋霞免费鲁丝片| 欧美另类一区| 免费黄色在线免费观看| 日韩精品免费视频一区二区三区| 成年美女黄网站色视频大全免费| 最近最新中文字幕免费大全7| 国产免费视频播放在线视频| 丝袜在线中文字幕| 在线精品无人区一区二区三| 看非洲黑人一级黄片| 日本av手机在线免费观看| 国产乱来视频区| 日韩欧美一区视频在线观看| 精品福利永久在线观看| 国产精品欧美亚洲77777| 超碰成人久久| 巨乳人妻的诱惑在线观看| 99精国产麻豆久久婷婷| 七月丁香在线播放| 天天躁日日躁夜夜躁夜夜| 久久人人97超碰香蕉20202| 亚洲久久久国产精品| 亚洲精品视频女| 日韩av免费高清视频| 色播在线永久视频| 在线亚洲精品国产二区图片欧美| 99久久人妻综合| 国产亚洲av片在线观看秒播厂| 91成人精品电影| 一级爰片在线观看| a级片在线免费高清观看视频| 国产成人一区二区在线| 国产av码专区亚洲av| 国产成人精品无人区| 国产精品久久久久久人妻精品电影 | 一级a爱视频在线免费观看| 亚洲精品日本国产第一区| 国产精品熟女久久久久浪| 99久国产av精品国产电影| 免费在线观看视频国产中文字幕亚洲 | 国产精品香港三级国产av潘金莲 | 久久天堂一区二区三区四区| 尾随美女入室| 99国产综合亚洲精品| 黄色视频不卡| 亚洲久久久国产精品| 熟女av电影| 欧美激情 高清一区二区三区| 美女午夜性视频免费| 在线观看人妻少妇| 另类精品久久| 韩国高清视频一区二区三区| 久久久久久久大尺度免费视频| 18禁动态无遮挡网站| 91成人精品电影| 日韩一区二区视频免费看| 国产一卡二卡三卡精品 | 两性夫妻黄色片| 最近最新中文字幕免费大全7| 午夜激情久久久久久久| 女人精品久久久久毛片| 中文精品一卡2卡3卡4更新| 国产高清国产精品国产三级| 久久久久人妻精品一区果冻| 一本色道久久久久久精品综合| 男女边吃奶边做爰视频| 久久热在线av| 午夜日韩欧美国产| 亚洲欧美清纯卡通| 色视频在线一区二区三区| 成人漫画全彩无遮挡| 成人18禁高潮啪啪吃奶动态图| www.熟女人妻精品国产| 热re99久久国产66热| 亚洲精品av麻豆狂野| 国产成人啪精品午夜网站| 午夜激情久久久久久久| 777米奇影视久久| 青青草视频在线视频观看| 日日爽夜夜爽网站| 亚洲精品日本国产第一区| 99热全是精品| 伊人久久国产一区二区| 天天影视国产精品| 999精品在线视频| 国产午夜精品一二区理论片| 亚洲一区中文字幕在线| 亚洲一级一片aⅴ在线观看| 一级,二级,三级黄色视频| 伦理电影大哥的女人| 久久天堂一区二区三区四区| 人体艺术视频欧美日本| 只有这里有精品99| 日韩中文字幕视频在线看片| 久久久久久久精品精品| 免费观看性生交大片5| 夜夜骑夜夜射夜夜干| 在线免费观看不下载黄p国产| 晚上一个人看的免费电影| 国产一区有黄有色的免费视频| 亚洲视频免费观看视频| 又大又黄又爽视频免费| 亚洲四区av| 美女中出高潮动态图| 美女脱内裤让男人舔精品视频| 亚洲免费av在线视频| av不卡在线播放| 国产淫语在线视频| 国产精品无大码| 亚洲,一卡二卡三卡| 女人被躁到高潮嗷嗷叫费观| 亚洲成色77777| 一区二区av电影网| 国产精品嫩草影院av在线观看| avwww免费| 亚洲精品自拍成人| 一二三四中文在线观看免费高清| 男女无遮挡免费网站观看| 久久久久久久久久久久大奶| 亚洲美女视频黄频| 色综合欧美亚洲国产小说| 校园人妻丝袜中文字幕| 国产成人免费无遮挡视频| av在线播放精品| 一本大道久久a久久精品| 性色av一级| 精品视频人人做人人爽| 伦理电影大哥的女人| 亚洲精品久久久久久婷婷小说| 99热网站在线观看| 精品一区二区三区av网在线观看 | 国产男女内射视频| 国产免费视频播放在线视频| 丰满饥渴人妻一区二区三| 久热这里只有精品99| 国产成人精品在线电影| 香蕉国产在线看| 黄色毛片三级朝国网站| 久久狼人影院| 中文字幕人妻熟女乱码| 欧美国产精品va在线观看不卡| 一区二区日韩欧美中文字幕| 成人国产麻豆网| 亚洲欧美清纯卡通| 99re6热这里在线精品视频| 久久久精品国产亚洲av高清涩受| 最近的中文字幕免费完整| 国产伦理片在线播放av一区| 女人精品久久久久毛片| 免费看av在线观看网站| 亚洲av日韩在线播放| 制服人妻中文乱码| 天天躁狠狠躁夜夜躁狠狠躁| 青草久久国产| 热re99久久国产66热| 国产色婷婷99| 国产极品粉嫩免费观看在线| 99九九在线精品视频| 亚洲欧美色中文字幕在线| a级毛片在线看网站| 日本爱情动作片www.在线观看| 18禁动态无遮挡网站| 一级毛片黄色毛片免费观看视频| 国产片内射在线| 欧美成人精品欧美一级黄| 成人亚洲欧美一区二区av| 熟女少妇亚洲综合色aaa.| 国产无遮挡羞羞视频在线观看| 久久精品人人爽人人爽视色| 高清av免费在线| 建设人人有责人人尽责人人享有的| 美女高潮到喷水免费观看| 18禁裸乳无遮挡动漫免费视频| 亚洲精品成人av观看孕妇| 久久精品国产亚洲av涩爱| 国产成人免费观看mmmm| 老汉色∧v一级毛片| 成人亚洲欧美一区二区av| 夫妻午夜视频| 精品亚洲成国产av| 国产精品无大码| 好男人视频免费观看在线| 日韩电影二区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品日本国产第一区| 午夜福利,免费看| 日日撸夜夜添| 国产成人欧美| 母亲3免费完整高清在线观看| 国产成人av激情在线播放| videos熟女内射| 午夜日本视频在线| 欧美人与性动交α欧美软件| 欧美日本中文国产一区发布| www.精华液| 纯流量卡能插随身wifi吗| 国产成人精品久久二区二区91 | 一区二区日韩欧美中文字幕| 99热网站在线观看| 亚洲三区欧美一区| 国产成人精品无人区| 啦啦啦中文免费视频观看日本| 亚洲第一av免费看| 亚洲av成人精品一二三区| 男人爽女人下面视频在线观看| 亚洲精品美女久久久久99蜜臀 | 男人爽女人下面视频在线观看| 久久久久精品人妻al黑| www.熟女人妻精品国产| 亚洲欧美色中文字幕在线| 国产精品无大码| 超碰成人久久| 性高湖久久久久久久久免费观看| 成人影院久久| 国产一区有黄有色的免费视频| 黄片无遮挡物在线观看| 曰老女人黄片| 人人妻人人澡人人爽人人夜夜| 狂野欧美激情性xxxx| 国产精品一区二区精品视频观看| 亚洲成av片中文字幕在线观看| 日日撸夜夜添| 欧美成人午夜精品| 亚洲国产精品一区二区三区在线| 中文字幕高清在线视频| 一本一本久久a久久精品综合妖精| 日本wwww免费看| 久久精品国产综合久久久| 亚洲精品成人av观看孕妇| 国产色婷婷99| 人人妻人人澡人人爽人人夜夜| 国产精品一国产av| 亚洲av电影在线观看一区二区三区| 国产黄色视频一区二区在线观看| 国产色婷婷99| 欧美日韩亚洲综合一区二区三区_| 99久久精品国产亚洲精品| 中文字幕av电影在线播放| 18禁动态无遮挡网站| 亚洲在久久综合| 国产精品国产三级专区第一集| 伦理电影免费视频| 狂野欧美激情性bbbbbb| 国产不卡av网站在线观看| 久久99精品国语久久久| 男女国产视频网站| 亚洲精品中文字幕在线视频| 久久免费观看电影| 欧美 亚洲 国产 日韩一| 捣出白浆h1v1| 丰满乱子伦码专区| 久久久国产欧美日韩av| 久久久精品国产亚洲av高清涩受| 天堂俺去俺来也www色官网| 亚洲成人av在线免费| 精品一区二区三区av网在线观看 | 91国产中文字幕| 亚洲精品日本国产第一区| 19禁男女啪啪无遮挡网站| 男女国产视频网站| 一区二区日韩欧美中文字幕| 亚洲免费av在线视频| 国产乱来视频区| 亚洲成av片中文字幕在线观看| 久久精品亚洲av国产电影网| 日韩一卡2卡3卡4卡2021年| 美女扒开内裤让男人捅视频| 欧美国产精品va在线观看不卡| 久久久欧美国产精品| 中文欧美无线码| av在线播放精品| 国产精品久久久久久久久免| 十八禁网站网址无遮挡| 久久精品亚洲av国产电影网| 日日摸夜夜添夜夜爱| 狂野欧美激情性bbbbbb| 国产精品麻豆人妻色哟哟久久| 国产精品av久久久久免费| 欧美人与性动交α欧美精品济南到| 免费av中文字幕在线| 亚洲精品在线美女| 国产一区二区激情短视频 | 高清欧美精品videossex| 午夜福利乱码中文字幕| 免费人妻精品一区二区三区视频| 日本av手机在线免费观看| 欧美日韩亚洲高清精品| 老司机影院毛片| 桃花免费在线播放| 亚洲一区中文字幕在线| 热99久久久久精品小说推荐| 欧美日韩成人在线一区二区| 日本av手机在线免费观看| 黑人巨大精品欧美一区二区蜜桃| 国产熟女欧美一区二区| 日韩av免费高清视频| 亚洲三区欧美一区| 精品人妻熟女毛片av久久网站| 欧美日韩成人在线一区二区| 亚洲色图综合在线观看| 久久久久精品人妻al黑| 搡老乐熟女国产| 精品一区二区三卡| 国产一区亚洲一区在线观看| 少妇精品久久久久久久| 爱豆传媒免费全集在线观看| 波多野结衣av一区二区av| bbb黄色大片| 美女视频免费永久观看网站| 亚洲色图 男人天堂 中文字幕| 一二三四在线观看免费中文在| 免费观看av网站的网址| 我要看黄色一级片免费的| 午夜福利在线免费观看网站| 久久韩国三级中文字幕| 亚洲中文av在线| 久久精品久久精品一区二区三区| 人妻一区二区av| 欧美日韩一级在线毛片| 黄色视频在线播放观看不卡| 亚洲av欧美aⅴ国产| av不卡在线播放| 亚洲成人av在线免费| 亚洲精品久久午夜乱码| 美女午夜性视频免费| 欧美日韩视频精品一区| 国产成人欧美| 亚洲精品日韩在线中文字幕| 天天躁日日躁夜夜躁夜夜| 欧美人与善性xxx| 国产片内射在线| avwww免费| 日韩中文字幕视频在线看片| 天天躁狠狠躁夜夜躁狠狠躁| 精品久久蜜臀av无| 国产色婷婷99| 国产亚洲午夜精品一区二区久久| 精品一品国产午夜福利视频| 一边亲一边摸免费视频| 男人添女人高潮全过程视频| 狂野欧美激情性xxxx| 91成人精品电影| 少妇人妻久久综合中文| 91成人精品电影| 操美女的视频在线观看| 女人高潮潮喷娇喘18禁视频| 成人国语在线视频| tube8黄色片| av片东京热男人的天堂| 少妇精品久久久久久久| 成人三级做爰电影| 日韩欧美一区视频在线观看| 国产视频首页在线观看| 国产伦理片在线播放av一区| 久久综合国产亚洲精品| 久久久久精品人妻al黑| 丁香六月欧美| 久久女婷五月综合色啪小说| 一区二区三区乱码不卡18| 午夜福利在线免费观看网站| 国语对白做爰xxxⅹ性视频网站| 天天添夜夜摸| 欧美人与性动交α欧美软件| 国产成人精品久久二区二区91 | 美女高潮到喷水免费观看| 免费观看性生交大片5| 国产一区二区三区综合在线观看| 国产99久久九九免费精品| 亚洲久久久国产精品| 亚洲国产精品一区三区| 深夜精品福利| 久久99热这里只频精品6学生| 老汉色av国产亚洲站长工具| 国精品久久久久久国模美| www.av在线官网国产| 亚洲,一卡二卡三卡| 国产在线免费精品| 精品国产国语对白av| 国产精品一区二区在线观看99| 91精品国产国语对白视频| 五月开心婷婷网| 国产成人免费无遮挡视频| 看十八女毛片水多多多| 欧美乱码精品一区二区三区| 狠狠婷婷综合久久久久久88av| 免费在线观看视频国产中文字幕亚洲 | 黄色视频在线播放观看不卡| 美女大奶头黄色视频| av女优亚洲男人天堂| a级毛片在线看网站| 水蜜桃什么品种好| 国产黄频视频在线观看| 久久久久久人人人人人| 日韩精品有码人妻一区| 国产乱人偷精品视频| 欧美日韩一级在线毛片| 女人高潮潮喷娇喘18禁视频| 中文精品一卡2卡3卡4更新| 精品一品国产午夜福利视频| 亚洲国产精品一区三区| 岛国毛片在线播放| 女人爽到高潮嗷嗷叫在线视频| 欧美激情高清一区二区三区 | 婷婷色麻豆天堂久久| 麻豆乱淫一区二区| 中文字幕人妻熟女乱码| 精品人妻一区二区三区麻豆| xxx大片免费视频| 在线观看免费视频网站a站| 久久久精品94久久精品| 国产一卡二卡三卡精品 | h视频一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 在线亚洲精品国产二区图片欧美| 国产精品秋霞免费鲁丝片| 国产野战对白在线观看| 啦啦啦在线观看免费高清www| 久久精品aⅴ一区二区三区四区| 国产高清国产精品国产三级| 成年人免费黄色播放视频| av卡一久久| 色网站视频免费| 热99国产精品久久久久久7| 国产成人精品无人区| 免费黄色在线免费观看| 日本av免费视频播放| 国产爽快片一区二区三区| av片东京热男人的天堂| 久久婷婷青草| 观看美女的网站| 国产精品.久久久| netflix在线观看网站| av.在线天堂| 午夜激情av网站| 免费人妻精品一区二区三区视频| 亚洲专区中文字幕在线 | 欧美黑人精品巨大| 丰满迷人的少妇在线观看| 精品人妻熟女毛片av久久网站| 热re99久久国产66热| 亚洲第一青青草原| 毛片一级片免费看久久久久| av.在线天堂| 国产黄色视频一区二区在线观看| 波多野结衣一区麻豆| av线在线观看网站| 免费久久久久久久精品成人欧美视频| kizo精华| 少妇猛男粗大的猛烈进出视频| 天天添夜夜摸| 水蜜桃什么品种好| 哪个播放器可以免费观看大片| 99久久人妻综合| 成人三级做爰电影| 亚洲,欧美,日韩| 好男人视频免费观看在线| 亚洲av日韩在线播放| 欧美成人午夜精品| 亚洲av成人精品一二三区| 在线观看一区二区三区激情| 久久影院123| 五月开心婷婷网| 2021少妇久久久久久久久久久| 亚洲精品自拍成人| 亚洲精品国产av蜜桃| 桃花免费在线播放| 久久毛片免费看一区二区三区| 国产99久久九九免费精品| 丝袜美腿诱惑在线| 国产亚洲av高清不卡| 观看av在线不卡| 91成人精品电影| 激情五月婷婷亚洲| 国产成人精品福利久久| 日本欧美国产在线视频| 国产黄频视频在线观看| 伊人久久国产一区二区| 一级毛片电影观看| 伦理电影大哥的女人| 在现免费观看毛片| 免费黄网站久久成人精品| 亚洲免费av在线视频| 久久久精品免费免费高清| 久久精品久久精品一区二区三区| 亚洲第一区二区三区不卡| 韩国av在线不卡| 在线观看免费午夜福利视频| 亚洲第一av免费看| 综合色丁香网| 黄频高清免费视频| 亚洲精品在线美女| 成人国语在线视频| 亚洲,欧美精品.| 9191精品国产免费久久| 日韩av在线免费看完整版不卡| 晚上一个人看的免费电影| 各种免费的搞黄视频| 熟女av电影| av国产久精品久网站免费入址| 伊人亚洲综合成人网| 国产精品蜜桃在线观看| 又粗又硬又长又爽又黄的视频| 亚洲国产精品999| 亚洲国产欧美网| 国产精品 国内视频| 纯流量卡能插随身wifi吗| 国产男人的电影天堂91| 精品人妻熟女毛片av久久网站| 又黄又粗又硬又大视频| 青草久久国产| 亚洲国产毛片av蜜桃av| 91老司机精品| 精品国产乱码久久久久久男人| 欧美成人午夜精品| e午夜精品久久久久久久| 久久毛片免费看一区二区三区| 毛片一级片免费看久久久久| 久久久亚洲精品成人影院| 18禁动态无遮挡网站| 久久久国产一区二区| 亚洲成av片中文字幕在线观看| 亚洲精品国产av成人精品| 只有这里有精品99| 免费黄频网站在线观看国产| 国产一区二区 视频在线| 老司机影院毛片| 亚洲五月色婷婷综合| 日韩不卡一区二区三区视频在线| 两个人免费观看高清视频| 毛片一级片免费看久久久久| 亚洲激情五月婷婷啪啪| 99久久人妻综合| 欧美精品av麻豆av| 亚洲av电影在线进入| 精品少妇一区二区三区视频日本电影 | 亚洲 欧美一区二区三区| 国产伦人伦偷精品视频| 老汉色∧v一级毛片| 在线观看三级黄色| 成年人免费黄色播放视频| 美女国产高潮福利片在线看| 久久久久精品久久久久真实原创| 成人午夜精彩视频在线观看| 亚洲国产欧美网| 丝袜人妻中文字幕| 国产日韩欧美亚洲二区| 久久综合国产亚洲精品| 亚洲欧美一区二区三区黑人| 国产精品久久久久久久久免| 亚洲国产中文字幕在线视频| 国产成人精品在线电影| 麻豆乱淫一区二区| 久久久久视频综合| 亚洲在久久综合| 国产精品久久久久久精品电影小说| 日本一区二区免费在线视频| 久久精品aⅴ一区二区三区四区| 亚洲成人手机| 欧美日韩成人在线一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 久久青草综合色| 天天操日日干夜夜撸| 日韩制服丝袜自拍偷拍| 久久ye,这里只有精品| 国产高清不卡午夜福利| 国产乱来视频区| 精品久久蜜臀av无| 国产成人精品福利久久| 少妇 在线观看| 最近最新中文字幕免费大全7| 久久精品国产亚洲av涩爱| 9热在线视频观看99| 日韩制服骚丝袜av| 宅男免费午夜| 在现免费观看毛片| 成人国产麻豆网| 1024香蕉在线观看| 不卡av一区二区三区| 国产欧美日韩综合在线一区二区| 日韩欧美一区视频在线观看| 欧美乱码精品一区二区三区| 亚洲熟女精品中文字幕| 男女无遮挡免费网站观看| 亚洲婷婷狠狠爱综合网| 欧美精品一区二区免费开放| 国产xxxxx性猛交| 中文字幕高清在线视频| 又大又爽又粗| 亚洲第一av免费看| av卡一久久| 中文欧美无线码| 亚洲精品久久成人aⅴ小说| 精品少妇内射三级| 街头女战士在线观看网站| 久久狼人影院| 久久精品国产综合久久久| 国产精品一国产av| 人妻 亚洲 视频|