• <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í)別方法研究
    av有码第一页| 免费在线观看完整版高清| 热99国产精品久久久久久7| 精品一区二区三区四区五区乱码 | 国产精品久久久久久精品电影小说| 日日摸夜夜添夜夜爱| 大片电影免费在线观看免费| 在线观看三级黄色| 国产成人免费观看mmmm| 欧美另类一区| 一区二区三区激情视频| 九九爱精品视频在线观看| 熟女av电影| 国产成人精品久久二区二区91 | 亚洲国产精品成人久久小说| 亚洲五月色婷婷综合| 一级黄片播放器| 久久久精品区二区三区| av不卡在线播放| 婷婷色综合www| 亚洲 欧美一区二区三区| www.熟女人妻精品国产| www.精华液| 伦理电影免费视频| 我的亚洲天堂| 最近手机中文字幕大全| 亚洲熟女毛片儿| 五月天丁香电影| 无遮挡黄片免费观看| 晚上一个人看的免费电影| 乱人伦中国视频| 天天影视国产精品| 欧美变态另类bdsm刘玥| 国产一区二区激情短视频 | 亚洲一码二码三码区别大吗| 男的添女的下面高潮视频| 久久天躁狠狠躁夜夜2o2o | 精品午夜福利在线看| 人人澡人人妻人| 亚洲 欧美一区二区三区| 午夜福利视频精品| 日本爱情动作片www.在线观看| 亚洲国产精品一区三区| www.熟女人妻精品国产| √禁漫天堂资源中文www| 亚洲国产中文字幕在线视频| 午夜老司机福利片| 制服诱惑二区| 永久免费av网站大全| 欧美xxⅹ黑人| 天天操日日干夜夜撸| 男女之事视频高清在线观看 | avwww免费| 精品一品国产午夜福利视频| 亚洲精品一二三| 性高湖久久久久久久久免费观看| 国产黄频视频在线观看| 国产精品亚洲av一区麻豆 | 精品一区二区三卡| 国产成人精品久久久久久| 亚洲av男天堂| 中文乱码字字幕精品一区二区三区| 日本91视频免费播放| 国产成人精品福利久久| 99热网站在线观看| 99久久综合免费| 九九爱精品视频在线观看| 女人久久www免费人成看片| 欧美人与性动交α欧美精品济南到| 久久久久久久久免费视频了| 亚洲国产成人一精品久久久| 美女脱内裤让男人舔精品视频| 极品人妻少妇av视频| 少妇猛男粗大的猛烈进出视频| 女的被弄到高潮叫床怎么办| 久久久国产精品麻豆| 在线观看免费午夜福利视频| 热re99久久国产66热| 18禁观看日本| 在线天堂最新版资源| 中国三级夫妇交换| 黄色视频在线播放观看不卡| 亚洲国产精品999| 日韩大码丰满熟妇| 亚洲精品久久午夜乱码| 国产欧美日韩一区二区三区在线| 美女脱内裤让男人舔精品视频| 一级毛片 在线播放| 午夜福利影视在线免费观看| 精品视频人人做人人爽| 超色免费av| 青春草亚洲视频在线观看| 国产黄色视频一区二区在线观看| 人人妻人人爽人人添夜夜欢视频| 一区二区av电影网| 欧美人与性动交α欧美软件| 国产日韩欧美在线精品| a级毛片在线看网站| 激情五月婷婷亚洲| 久热爱精品视频在线9| 考比视频在线观看| 一区在线观看完整版| 免费观看性生交大片5| 悠悠久久av| 国产欧美亚洲国产| 午夜日本视频在线| 中文精品一卡2卡3卡4更新| av不卡在线播放| 欧美 亚洲 国产 日韩一| 搡老乐熟女国产| 人人妻人人添人人爽欧美一区卜| 大片免费播放器 马上看| 欧美 亚洲 国产 日韩一| 老司机在亚洲福利影院| 久久久久国产精品人妻一区二区| 性色av一级| videos熟女内射| 国产精品久久久久久精品古装| 欧美日韩一级在线毛片| 国产成人精品无人区| 亚洲国产欧美一区二区综合| 老鸭窝网址在线观看| 精品国产露脸久久av麻豆| 亚洲欧美一区二区三区黑人| 丝袜人妻中文字幕| 婷婷色综合大香蕉| 国产男人的电影天堂91| 丰满迷人的少妇在线观看| 久久韩国三级中文字幕| 日韩欧美精品免费久久| av不卡在线播放| 妹子高潮喷水视频| 自线自在国产av| 亚洲欧美色中文字幕在线| 亚洲,欧美,日韩| 欧美激情高清一区二区三区 | 美女国产高潮福利片在线看| 久久狼人影院| 男女下面插进去视频免费观看| 久久精品亚洲av国产电影网| 国产97色在线日韩免费| 人人妻人人添人人爽欧美一区卜| 中文字幕人妻熟女乱码| 性高湖久久久久久久久免费观看| 成人免费观看视频高清| 午夜福利免费观看在线| 国产福利在线免费观看视频| 中文字幕人妻熟女乱码| 欧美国产精品一级二级三级| av天堂久久9| 日日爽夜夜爽网站| 国产成人免费无遮挡视频| 国产av国产精品国产| 午夜免费鲁丝| 亚洲国产av影院在线观看| 中国三级夫妇交换| 欧美日韩成人在线一区二区| 大香蕉久久成人网| 男人操女人黄网站| 中文字幕最新亚洲高清| 一本久久精品| 乱人伦中国视频| 久久鲁丝午夜福利片| www.av在线官网国产| 免费女性裸体啪啪无遮挡网站| 如日韩欧美国产精品一区二区三区| 人人妻人人澡人人看| 国产欧美日韩综合在线一区二区| 高清av免费在线| 亚洲婷婷狠狠爱综合网| 亚洲五月色婷婷综合| 晚上一个人看的免费电影| 国产欧美亚洲国产| 国产国语露脸激情在线看| netflix在线观看网站| 午夜福利视频在线观看免费| 肉色欧美久久久久久久蜜桃| 欧美精品一区二区免费开放| 不卡av一区二区三区| 少妇被粗大的猛进出69影院| 国产麻豆69| 欧美国产精品va在线观看不卡| 秋霞伦理黄片| 午夜av观看不卡| 亚洲婷婷狠狠爱综合网| 久久精品aⅴ一区二区三区四区| 久久精品亚洲av国产电影网| 久久精品亚洲熟妇少妇任你| 成年人免费黄色播放视频| 操美女的视频在线观看| 精品国产超薄肉色丝袜足j| 不卡视频在线观看欧美| 国产精品欧美亚洲77777| 黄色 视频免费看| 中文字幕人妻丝袜制服| 久久青草综合色| 少妇精品久久久久久久| 日本色播在线视频| 国产 精品1| 啦啦啦 在线观看视频| 18禁裸乳无遮挡动漫免费视频| 亚洲欧美日韩另类电影网站| 欧美97在线视频| 黄片无遮挡物在线观看| 99久久综合免费| 日韩欧美一区视频在线观看| 精品一区二区免费观看| 久久久精品区二区三区| 交换朋友夫妻互换小说| 亚洲情色 制服丝袜| 国产日韩欧美亚洲二区| 午夜av观看不卡| 2021少妇久久久久久久久久久| 人妻 亚洲 视频| 免费观看a级毛片全部| 又黄又粗又硬又大视频| 五月开心婷婷网| 亚洲国产中文字幕在线视频| 极品少妇高潮喷水抽搐| 如日韩欧美国产精品一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 国产一卡二卡三卡精品 | 好男人视频免费观看在线| 日韩欧美一区视频在线观看| e午夜精品久久久久久久| 麻豆乱淫一区二区| 国产视频首页在线观看| 下体分泌物呈黄色| 99国产综合亚洲精品| 夫妻午夜视频| 久久久久久人妻| 亚洲国产欧美一区二区综合| 国产精品久久久久成人av| 中文字幕人妻丝袜一区二区 | 亚洲视频免费观看视频| 桃花免费在线播放| 女的被弄到高潮叫床怎么办| 欧美精品人与动牲交sv欧美| 日韩欧美精品免费久久| 亚洲,欧美精品.| 成人国语在线视频| 九九爱精品视频在线观看| 亚洲精品自拍成人| 亚洲精品日本国产第一区| 日本黄色日本黄色录像| av天堂久久9| videos熟女内射| 亚洲美女黄色视频免费看| 免费在线观看视频国产中文字幕亚洲 | 美女大奶头黄色视频| 久久国产亚洲av麻豆专区| 午夜福利一区二区在线看| 亚洲伊人久久精品综合| 国产一级毛片在线| 久久99精品国语久久久| 国产精品99久久99久久久不卡 | 亚洲四区av| 久久鲁丝午夜福利片| 国产亚洲午夜精品一区二区久久| av在线老鸭窝| 69精品国产乱码久久久| 免费日韩欧美在线观看| 夫妻性生交免费视频一级片| 中文字幕av电影在线播放| 午夜福利一区二区在线看| 国产成人91sexporn| 天堂俺去俺来也www色官网| 国产免费福利视频在线观看| 久久精品国产综合久久久| 久久av网站| 日韩中文字幕视频在线看片| 国产精品99久久99久久久不卡 | 啦啦啦中文免费视频观看日本| 亚洲精品国产av蜜桃| 色网站视频免费| 99精品久久久久人妻精品| 欧美久久黑人一区二区| 麻豆精品久久久久久蜜桃| 熟女av电影| 99国产综合亚洲精品| 操出白浆在线播放| 丝袜美足系列| 精品亚洲成国产av| 波多野结衣av一区二区av| 久久女婷五月综合色啪小说| 成人手机av| 国产激情久久老熟女| 一边摸一边做爽爽视频免费| 欧美黄色片欧美黄色片| 成人黄色视频免费在线看| 在线观看免费高清a一片| 亚洲av在线观看美女高潮| 精品免费久久久久久久清纯 | 汤姆久久久久久久影院中文字幕| 国产精品女同一区二区软件| 国产在线一区二区三区精| 在线天堂最新版资源| svipshipincom国产片| 超碰97精品在线观看| 亚洲av在线观看美女高潮| 搡老乐熟女国产| 激情视频va一区二区三区| 赤兔流量卡办理| 国产极品天堂在线| 久久综合国产亚洲精品| av在线app专区| 欧美在线一区亚洲| av天堂久久9| 一边摸一边抽搐一进一出视频| av有码第一页| 中国国产av一级| 高清在线视频一区二区三区| 大陆偷拍与自拍| 日韩不卡一区二区三区视频在线| 一区二区三区乱码不卡18| 在线观看人妻少妇| 日韩人妻精品一区2区三区| 制服人妻中文乱码| 热re99久久精品国产66热6| 免费看不卡的av| 男女边摸边吃奶| 亚洲av电影在线观看一区二区三区| 亚洲欧美激情在线| 欧美日本中文国产一区发布| 色婷婷av一区二区三区视频| 曰老女人黄片| 黄片小视频在线播放| 日本色播在线视频| 日韩免费高清中文字幕av| 赤兔流量卡办理| av福利片在线| 老司机深夜福利视频在线观看 | 大码成人一级视频| 欧美成人午夜精品| 波野结衣二区三区在线| 91精品伊人久久大香线蕉| 国产精品二区激情视频| 久久人人爽av亚洲精品天堂| 人妻一区二区av| 久久久久精品国产欧美久久久 | 国产精品熟女久久久久浪| 久久女婷五月综合色啪小说| 色播在线永久视频| 老司机影院毛片| 91精品国产国语对白视频| 人妻一区二区av| 国产视频首页在线观看| 国产一卡二卡三卡精品 | 无遮挡黄片免费观看| 亚洲av成人精品一二三区| 亚洲精品国产区一区二| 亚洲美女黄色视频免费看| 韩国av在线不卡| 一级片免费观看大全| 成年人午夜在线观看视频| 国产伦人伦偷精品视频| 老司机深夜福利视频在线观看 | 久久久久精品人妻al黑| 啦啦啦中文免费视频观看日本| 国产毛片在线视频| 青青草视频在线视频观看| 国产不卡av网站在线观看| 欧美中文综合在线视频| 精品人妻在线不人妻| 大码成人一级视频| 日韩精品免费视频一区二区三区| 日本猛色少妇xxxxx猛交久久| 日韩制服丝袜自拍偷拍| 久久精品亚洲av国产电影网| 亚洲精品中文字幕在线视频| 国产熟女午夜一区二区三区| 晚上一个人看的免费电影| 国产精品一国产av| 黄色 视频免费看| 亚洲熟女精品中文字幕| 国产欧美日韩一区二区三区在线| 亚洲欧洲国产日韩| 国产一级毛片在线| 久久精品久久久久久噜噜老黄| 亚洲国产成人一精品久久久| 亚洲av成人精品一二三区| 欧美成人精品欧美一级黄| 久久性视频一级片| 悠悠久久av| 69精品国产乱码久久久| 免费黄网站久久成人精品| 成人国产麻豆网| 三上悠亚av全集在线观看| 波多野结衣一区麻豆| 国产精品久久久久久精品古装| 中文欧美无线码| 韩国av在线不卡| 超色免费av| 激情视频va一区二区三区| 久久性视频一级片| 久久久久久人人人人人| 欧美日韩一级在线毛片| 国产在线一区二区三区精| 精品一区二区免费观看| 欧美精品高潮呻吟av久久| 可以免费在线观看a视频的电影网站 | 麻豆av在线久日| 侵犯人妻中文字幕一二三四区| 亚洲av欧美aⅴ国产| 黄色视频在线播放观看不卡| 国产精品国产三级专区第一集| 麻豆乱淫一区二区| 制服诱惑二区| 嫩草影视91久久| 日本vs欧美在线观看视频| 在线亚洲精品国产二区图片欧美| 性少妇av在线| 国产一区有黄有色的免费视频| 亚洲国产欧美在线一区| 欧美黑人精品巨大| 国产精品一区二区在线不卡| 老汉色av国产亚洲站长工具| 中文字幕高清在线视频| 中文字幕人妻丝袜制服| 午夜激情av网站| 国产黄色免费在线视频| 欧美亚洲日本最大视频资源| 久久久久精品久久久久真实原创| 91成人精品电影| 777米奇影视久久| 欧美97在线视频| √禁漫天堂资源中文www| 国产不卡av网站在线观看| 91老司机精品| 午夜福利免费观看在线| 2021少妇久久久久久久久久久| 纵有疾风起免费观看全集完整版| 欧美成人午夜精品| 最黄视频免费看| 韩国精品一区二区三区| 国产精品人妻久久久影院| 日韩电影二区| 亚洲欧洲精品一区二区精品久久久 | 久久99一区二区三区| 国产精品一区二区精品视频观看| 这个男人来自地球电影免费观看 | 日韩一卡2卡3卡4卡2021年| av在线app专区| 777米奇影视久久| 成人午夜精彩视频在线观看| 视频区图区小说| 亚洲国产精品国产精品| 一级片'在线观看视频| 赤兔流量卡办理| 青青草视频在线视频观看| 国产一区二区激情短视频 | 欧美日韩视频高清一区二区三区二| 熟女av电影| 一区二区三区四区激情视频| 国产av码专区亚洲av| 亚洲综合精品二区| 老鸭窝网址在线观看| 一区二区三区激情视频| 少妇人妻 视频| 亚洲国产毛片av蜜桃av| 国产片内射在线| 中文字幕人妻丝袜一区二区 | 永久免费av网站大全| 一级a爱视频在线免费观看| 超色免费av| 欧美在线一区亚洲| 国产乱来视频区| 男女之事视频高清在线观看 | 久久精品国产亚洲av高清一级| 啦啦啦视频在线资源免费观看| 赤兔流量卡办理| 久久久久国产一级毛片高清牌| 国产精品蜜桃在线观看| 国产福利在线免费观看视频| 国产成人a∨麻豆精品| 在线精品无人区一区二区三| 激情五月婷婷亚洲| 日韩电影二区| 中文字幕色久视频| www.自偷自拍.com| 黄片无遮挡物在线观看| 亚洲精品第二区| 国产免费现黄频在线看| 天堂8中文在线网| 人妻 亚洲 视频| 大码成人一级视频| 亚洲精品一区蜜桃| 香蕉国产在线看| 国产乱来视频区| 满18在线观看网站| 欧美日韩国产mv在线观看视频| 国产精品 国内视频| 亚洲成人一二三区av| 亚洲伊人久久精品综合| 久热这里只有精品99| 精品少妇久久久久久888优播| 色网站视频免费| 国产精品一区二区精品视频观看| 别揉我奶头~嗯~啊~动态视频 | 男人添女人高潮全过程视频| 亚洲色图 男人天堂 中文字幕| 天堂俺去俺来也www色官网| 51午夜福利影视在线观看| 日日爽夜夜爽网站| 午夜福利,免费看| 日韩中文字幕视频在线看片| 新久久久久国产一级毛片| 久久久久久人人人人人| 国产精品国产三级国产专区5o| 99国产综合亚洲精品| 99九九在线精品视频| 成人亚洲精品一区在线观看| 中文字幕av电影在线播放| 成人亚洲欧美一区二区av| 黄色毛片三级朝国网站| 免费高清在线观看日韩| 91精品国产国语对白视频| 老汉色av国产亚洲站长工具| 国精品久久久久久国模美| 成年人免费黄色播放视频| 中文字幕av电影在线播放| 成人免费观看视频高清| 午夜日韩欧美国产| 大陆偷拍与自拍| 亚洲四区av| 亚洲av欧美aⅴ国产| 老司机亚洲免费影院| 波野结衣二区三区在线| 无限看片的www在线观看| 天堂中文最新版在线下载| 又黄又粗又硬又大视频| 久久精品国产亚洲av高清一级| 五月天丁香电影| 精品福利永久在线观看| 咕卡用的链子| 亚洲七黄色美女视频| 777久久人妻少妇嫩草av网站| 一级,二级,三级黄色视频| 亚洲精品成人av观看孕妇| 亚洲国产欧美一区二区综合| 一级黄片播放器| 国产精品成人在线| 纯流量卡能插随身wifi吗| 国产麻豆69| 免费少妇av软件| 免费av中文字幕在线| 久久天堂一区二区三区四区| 亚洲精品国产色婷婷电影| 欧美黑人精品巨大| 国产片内射在线| 精品人妻熟女毛片av久久网站| 中文字幕av电影在线播放| 久久精品人人爽人人爽视色| 99久久综合免费| 国产av精品麻豆| 老司机在亚洲福利影院| 成人手机av| 精品亚洲乱码少妇综合久久| 色婷婷av一区二区三区视频| 午夜福利乱码中文字幕| 亚洲国产精品一区三区| 国产成人精品福利久久| 最近手机中文字幕大全| 满18在线观看网站| 色婷婷久久久亚洲欧美| 母亲3免费完整高清在线观看| 黄色视频在线播放观看不卡| 丝袜人妻中文字幕| 久久久久久久大尺度免费视频| 黄片无遮挡物在线观看| 亚洲精品在线美女| 日本一区二区免费在线视频| 在线观看人妻少妇| 一级片'在线观看视频| 午夜福利一区二区在线看| 亚洲情色 制服丝袜| 精品人妻熟女毛片av久久网站| 亚洲免费av在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲美女视频黄频| 国产精品一二三区在线看| 亚洲国产精品一区二区三区在线| 亚洲精品日本国产第一区| 老鸭窝网址在线观看| 亚洲成人一二三区av| 男女边摸边吃奶| 国产精品久久久久久精品电影小说| 永久免费av网站大全| 亚洲成人一二三区av| 国产精品成人在线| 女性被躁到高潮视频| 美女大奶头黄色视频| 一区二区三区激情视频| 亚洲av成人精品一二三区| 久久精品熟女亚洲av麻豆精品| 亚洲国产精品一区三区| 亚洲在久久综合| xxxhd国产人妻xxx| 夫妻性生交免费视频一级片| 热99国产精品久久久久久7| 日韩av不卡免费在线播放| 青春草国产在线视频| av免费观看日本| 国产深夜福利视频在线观看| 又粗又硬又长又爽又黄的视频| 国产精品国产av在线观看| 一区二区三区四区激情视频| 国产av精品麻豆| 一级a爱视频在线免费观看| 操出白浆在线播放| 美女午夜性视频免费| 精品一区在线观看国产| 精品人妻熟女毛片av久久网站| 天天躁日日躁夜夜躁夜夜| 亚洲av日韩在线播放|