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

    大興安嶺南段毛登Sn-Cu礦床巖漿作用對成礦制約:年代學、地球化學及Sr-Nd-Pb同位素證據(jù)

    2021-08-24 05:28:44季根源江思宏李高峰易錦俊張莉莉劉翼飛
    大地構(gòu)造與成礦學 2021年4期
    關鍵詞:成礦

    季根源, 江思宏, 李高峰, 易錦俊, 張莉莉, 劉翼飛

    大興安嶺南段毛登Sn-Cu礦床巖漿作用對成礦制約:年代學、地球化學及Sr-Nd-Pb同位素證據(jù)

    季根源1, 2, 江思宏1*, 李高峰1, 3, 易錦俊2, 張莉莉1, 劉翼飛1

    (1.中國地質(zhì)科學院 礦產(chǎn)資源研究所, 自然資源部成礦作用與資源評價重點實驗室, 北京 100037; 2.自然資源實物地質(zhì)資料中心, 河北 廊坊 065201; 3.北京大學 造山帶與地殼演化重點實驗室, 北京 100871)

    毛登Sn-Cu礦床位于大興安嶺南段內(nèi)蒙古錫林浩特境內(nèi), 其Sn、Cu儲量均達到中型。本文利用LA-ICP-MS鋯石U-Pb法測得與成礦關系密切的花崗斑巖成巖年齡為138±0.6 Ma; 利用LA-ICP-MS錫石U-Pb法測得錫石?硫化物?石英脈型礦石中錫石的Tera-Wasserburg諧和年齡為139±3.2 Ma, 表明花崗斑巖體形成年齡與毛登礦床成礦年齡在誤差范圍內(nèi)一致, 均形成于早白堊世?;◢彴邘r富硅(SiO2=73.18%~77.3%)、富堿(Na2O+K2O=8.00%~8.69%)、富鉀(K2O=4.80%~ 5.22%)、貧鋁(Al2O3=11.67%~12.83%), 貧鈣(CaO=0.49%~0.95%)、鎂(MgO=0.12%~0.46%), 具有較高的FeOT/(MgO+FeOT)值(0.74~0.91)、10000×Ga/Al值(3.64~4.28); 微量元素虧損Ba、Sr、P、Ti等元素, 富集Rb、Th、K、Hf等元素; 稀土元素配分模式圖呈典型的“海鷗型”以及顯著的Eu負異常(δEu=0.05~0.2); 鋯石飽和溫度為817~861 ℃, 以上均說明花崗斑巖具有A型花崗巖特征。全巖Sr-Nd同位素結(jié)果顯示低初始(87Sr/86Sr)i值(0.702907~0.704506)、較高的Nd()值(+0.8~+4.4)以及年輕的二階段模式年齡(DM2=572~863 Ma), 說明花崗斑巖可能為新生地殼經(jīng)部分熔融后, 在上升過程中經(jīng)歷較強烈的結(jié)晶分異演化作用而形成。Pb同位素組成指示巖體和礦石的Pb同位素組成接近, 說明成礦物質(zhì)主要來源于巖漿。微量元素特征、構(gòu)造環(huán)境判別圖解顯示花崗斑巖具有A2型造山后花崗巖特征, 成巖環(huán)境處于伸展構(gòu)造背景。結(jié)合前人研究成果以及礦區(qū)地質(zhì)特征、成巖(礦)年代學、Sr-Nd-Pb同位素特征, 認為毛登礦床成礦物質(zhì)的富集與花崗斑巖體侵入密切相關, 礦床成礦物質(zhì)的富集規(guī)模受花崗斑巖結(jié)晶分異演化程度制約。

    大興安嶺南段; 毛登錫銅礦床; 錫石U-Pb年齡; 鋯石U-Pb年齡; Sr-Nd-Pb同位素

    0 引 言

    大興安嶺南段位于中亞造山東部的興蒙造山帶內(nèi)(圖1), 行政區(qū)隸屬于內(nèi)蒙古自治區(qū), 區(qū)內(nèi)大量Sn-Ag-Cu-Pb-Zn-Mo礦成礦作用與中生代巖漿活動關系密切(毛景文等, 2013)。在20世紀末期, 區(qū)內(nèi)已經(jīng)發(fā)現(xiàn)鬧牛山、蓮花山等中小型銅礦床以及黃崗、安樂、大井、寶蓋溝、查木罕等大中型錫多金屬礦床。隨著近年來地質(zhì)找礦工作的不斷深入, 白音查干、維拉斯托、拜仁達壩、白音諾爾等一批大型錫鉛鋅銀礦床的陸續(xù)發(fā)現(xiàn), 大興安嶺南段地區(qū)已逐漸成為我國北方地區(qū)重要錫多金屬礦產(chǎn)地, 顯示著巨大的錫資源找礦潛力(陳鄭輝等, 2015; Liu et al., 2016; 王春女等, 2016; Wang et al., 2017; 姚磊等, 2017)。大興安嶺南段地區(qū)豐富的礦產(chǎn)資源引起了地質(zhì)科研人員的重點關注, 有關該區(qū)的地學研究成果日益豐富, 主要集中在巖石地球化學以及成巖(礦)年代學方面。研究表明, 區(qū)內(nèi)與錫礦床成因有關的晚中生代花崗質(zhì)巖體多為正長花崗巖、似斑狀二長花崗巖、石英斑巖等, 多數(shù)屬于鈣堿性?高鈣堿性系列的準鋁質(zhì)?弱過鋁質(zhì)的A型花崗巖(邵濟安等, 2010; 周振華等, 2010; 葉天竺, 2013; Liu et al., 2016; 劉新等, 2017), 巖體形成于巖石圈伸展減薄構(gòu)造背景(郭鋒等, 2001; Fan et al., 2003; 呂志成等, 2004; 江思宏等, 2012; Ouyang et al., 2015); 成礦年齡與成巖年齡具有一致性, 集中在130~150 Ma(王京彬等, 2005; 王長明等, 2006; 翟德高等, 2012; 要梅娟等, 2016; Zhai et al., 2017)。

    圖1 大興安嶺南段地質(zhì)簡圖及錫多金屬礦床分布圖(據(jù)徐備等, 2014; Ouyang et al., 2015; Chen et al., 2017; 姚磊等, 2017)

    毛登Sn-Cu礦床位于內(nèi)蒙古錫林浩特市北東50 km處, 是20世紀大興安嶺南段西坡發(fā)現(xiàn)的一處中型Sn-Cu礦床。近年在其深部新發(fā)現(xiàn)大型鉬礦體, 鉬金屬資源量達15萬噸。毛登Sn-Cu礦床為大興安嶺南段的典型錫礦床之一, 但地質(zhì)科研工作相對薄弱。劉玉強(1996a, 1996b)、張巧梅等(2013)、袁小平等(2016)分析了毛登礦床地質(zhì)特征, 探討了控礦因素, 認為毛登礦床是受構(gòu)造控制的高?中溫熱液脈型Sn-Cu礦床。張曉玲和張文山(2010)、劉傳寶等(2018)研究認為區(qū)內(nèi)與礦化有關的花崗斑巖與火山熔巖?火山碎屑巖是燕山晚期巖漿活動同源同期異相產(chǎn)物。程天赦等(2014)、郭碩等(2019)研究認為毛登?小狐山錫多金屬礦田的巖漿侵入與成礦作用發(fā)生在早白堊世。本文采集毛登錫銅礦區(qū)內(nèi)典型礦石及與成礦關系密切的花崗斑巖樣品, 開展巖石地球化學、Sr-Nd-Pb同位素、鋯石和錫石LA-ICP-MS U-Pb年齡研究, 厘定成巖(礦)時代, 揭示巖體起源及其形成環(huán)境, 探討巖石成因與演化及其對成礦的制約, 對認識大興安嶺南段地區(qū)燕山晚期巖漿演化與錫多金屬礦床成礦作用具有重要意義。

    1 礦區(qū)和礦床地質(zhì)

    毛登Sn-Cu礦床位于大興安嶺南段, 地處古亞洲洋、古太平洋和蒙古?鄂霍茨克洋構(gòu)造體系疊加區(qū)域, 屬于中亞造山帶東部的興蒙造山帶, 北以二連?賀根山斷裂為界, 南為西拉木倫斷裂, 東臨嫩江斷裂(圖1)。區(qū)域內(nèi)主要出露下二疊統(tǒng)大石寨組、中二疊統(tǒng)哲斯組、下侏羅統(tǒng)紅旗組、上侏羅統(tǒng)滿克頭鄂博組及和瑪尼吐組、下白堊統(tǒng)白音高老組以及第四系。斷裂構(gòu)造主要呈北東?北北東向, 北西向次之, 兩者交匯處控制著巖漿巖的空間分布(張巧梅等, 2013)。區(qū)域內(nèi)晚古生代、中生代巖漿巖沿主構(gòu)造線呈近北東向分布, 中生代侵入巖發(fā)育最廣泛, 與成礦關系也最緊密(呂志成等, 2000; Ouyang et al., 2015)。

    毛登Sn-Cu礦床主要出露地層有下二疊統(tǒng)大石寨組、下侏羅統(tǒng)紅旗組、下白堊統(tǒng)白音高老組以及第四系(圖2)。大石寨組分布在礦區(qū)中部, 巖性為火山角礫巖、含火山角礫熔巖以及變質(zhì)粉砂巖等, 毛登Sn-Cu礦體主要賦存于巖體與大石寨組的外接觸帶。下侏羅統(tǒng)紅旗組分布在礦區(qū)北部, 巖性為砂巖、砂礫巖以及泥巖等。下白堊統(tǒng)白音高老組分布在礦區(qū)西北部, 巖性為酸性熔巖、凝灰?guī)r、凝灰質(zhì)粉砂巖等。

    圖2 毛登Sn-Cu礦床地質(zhì)簡圖(據(jù)石得鳳, 2007; 張巧梅等, 2013; 郭碩等, 2019)

    礦區(qū)內(nèi)北西向斷裂及節(jié)理裂隙是區(qū)內(nèi)主要容礦構(gòu)造, 控制著礦體的產(chǎn)出位置和形態(tài)。FⅠ-2、FⅡ- 35、FⅡ-46是三條主斷裂, 斷裂總體走向315°~335°, 傾向NE, 傾角60°~85°。斷裂之間發(fā)育大量密集的節(jié)理裂隙, 眾多的斷裂與節(jié)理裂隙構(gòu)成北西向斷裂?裂隙帶, 自西向東(離巖體的由遠至近)劃分為Ⅰ、Ⅱ、Ⅲ號斷裂?裂隙帶(圖2)。礦區(qū)東部2 km處阿魯寶格山發(fā)育燕山晚期花崗斑巖, 與成礦關系密切。毛登Sn-Cu礦體以密集礦脈狀形式賦存于Ⅰ、Ⅱ、Ⅲ號斷裂?裂隙帶中(圖2、3), 礦體產(chǎn)狀與受控斷裂、節(jié)理裂隙面基本一致, 連續(xù)性較差。Ⅰ號斷裂?裂隙帶在巖體西側(cè)0.5~1 km處, 礦脈以錫石?硫化物?石英脈為主, 沿走向長48~200 m, 沿傾向深23~200 m; Ⅱ號斷裂?裂隙帶在巖體西側(cè)0.5 km以內(nèi), 賦存于巖體與圍巖的外接觸帶中, 礦脈為錫石?石英脈或錫石?輝鉬礦?石英脈型, 沿走向長28~297 m, 沿傾向深23~166 m; Ⅲ號斷裂?裂隙帶在巖體內(nèi)接觸帶中, 礦脈為錫石?石英脈型, 沿走向長26~215 m, 沿傾向深40~118 m(劉玉強, 1996a, 1996b; 石得鳳, 2007; 張巧梅等, 2013)。

    圖3 毛登Sn-Cu礦床101號勘探線剖面圖(據(jù)河南省有色金屬地質(zhì)勘查總院, 2009)

    礦床的礦石礦物主要有錫石、黃銅礦、閃鋅礦、輝鉬礦、毒砂、黃鐵礦、黑鎢礦, 含少量的黝銅礦、方鉛礦、黃錫礦、磁黃鐵礦、硫銻鉛礦等。脈石礦物主要有黃玉、絹云母、白云母、石英、螢石, 含少量的獨居石、板鈦礦等(圖4)。礦石結(jié)構(gòu)主要有自形?半自形晶結(jié)構(gòu)、它形粒狀結(jié)構(gòu)、鑲嵌結(jié)構(gòu)、交代結(jié)構(gòu)、共生結(jié)構(gòu)、交代殘余結(jié)構(gòu)等。礦石構(gòu)造主要有細脈浸染狀、塊狀、脈狀、晶簇狀構(gòu)造等。礦床近礦蝕變普遍發(fā)育, Sn礦化與云英巖化、硅化、螢石化、黃玉化、絹云母化關系密切, 而Cu礦化與硅化、螢石化、碳酸鹽化關系密切。

    根據(jù)礦石組構(gòu)、光薄片鑒定、圍巖蝕變以及脈體之間穿插關系, 將毛登Sn-Cu礦床成礦過程歸納為4個成礦階段: 分別為①錫石?黑鎢礦?黃玉?石英階段, 為氣成?高溫熱液階段, 是重要的Sn成礦階段, 金屬礦物包括錫石、黑鎢礦, 及少量黃銅礦、毒砂、輝鉬礦和閃鋅礦等, 脈石礦物主要為石英、黃玉, 及少量的獨居石、板鈦礦、白云母等(圖4a、c、d)。②黃銅礦?閃鋅礦?輝鉬礦?螢石?石英階段, 為高溫熱液階段, 是重要的Cu成礦階段, 金屬礦物包括黃銅礦、閃鋅礦、輝鉬礦、黃錫礦、黃鐵礦、毒砂、自然鉍, 及少量磁黃鐵礦等, 脈石礦物主要為石英、螢石、白云母等(圖4b、e、f)。③閃鋅礦?方鉛礦?螢石?絹云母?石英階段, 為中溫熱液階段, 金屬礦物主要為閃鋅礦、方鉛礦、磁黃鐵礦、黃鐵礦, 及少量黃錫礦、黝銅礦、毒砂、硫銻鉛礦等, 脈石礦物主要為石英、螢石和絹云母等(圖4g、h)。④黃鐵礦?石英?碳酸鹽階段, 屬低溫熱液階段, 成礦較差, 金屬礦物僅見少量黃鐵礦, 脈石礦物主要為石英、方解石, 及少量綠泥石、藍閃石、高嶺石等(圖4i)。

    (a) 錫石+黃鐵礦+石英礦石; (b) 閃鋅礦+黃銅礦+黃鐵礦礦石; (c) 黃玉+石英脈(單偏光); (d) 石英脈中獨居石+黑鎢礦+錫石(正交偏光); (e) 輝鉬礦+毒砂+自然鉍(反射光); (f) 黃銅礦和黃錫礦交代閃鋅礦(反射光); (g) 黃銅礦、閃鋅礦和磁黃鐵礦共生(反射光); (h) 硫化物礦石中金屬礦物+螢石+絹云母(單偏光); (i) 綠泥石+藍閃石(正交偏光)。

    2 樣品采集與分析方法

    2.1 樣品采集

    巖石樣品采于毛登Sn-Cu礦床1號采區(qū)至2號采區(qū)的山坡(坐標: 44°10′30.61″N, 116°34′08.75″E), 均為花崗斑巖, 風化程度較弱, 具斑狀結(jié)構(gòu)?基質(zhì)霏細結(jié)構(gòu), 塊狀構(gòu)造(圖5)。斑晶由斜長石(5%~10%)、鉀長石(15%~25%)、石英(20%~25%)及少量角閃石、輝石假象暗色礦物(2%~5%)構(gòu)成。斜長石呈半自形板狀結(jié)構(gòu), 聚片雙晶清晰, 粒徑0.2~4.0 mm, 具輕微絹云母化、高嶺土化, 局部皂石化, 少量顯破碎狀特征。鉀長石呈半自形?近半自形板狀結(jié)構(gòu), 粒徑0.5~5.0 mm, 輕微高嶺土化, 鉀長石卡式雙晶清晰。石英呈半自形粒狀?它形結(jié)構(gòu), 粒徑0.05~4.5 mm, 見熔蝕現(xiàn)象, 少量石英斑晶具破碎狀特征。暗色礦物粒徑0.2~1.2 mm, 角閃石顯褐綠色, 多色性明顯; 輝石變?yōu)槲⒘iL英質(zhì)及針柱狀角閃石等的集合體假像, 發(fā)育少量角閃石反應邊?;|(zhì)主要由長英質(zhì)、暗色礦物組成, 粒徑多小于0.02 mm。長英質(zhì)包括斜長石、鉀長石、石英, 部分具霏細?微晶結(jié)構(gòu); 暗色礦物量少, 零散可見, 偶見黑云母、角閃石等。

    礦石樣品MDX-1來自毛登Sn-Cu礦1號采區(qū)(坐標: 44°10′30.61″N, 116°34′08.75″E), 為礦區(qū)內(nèi)第①成礦階段代表性的錫石?硫化物?石英脈型礦石, 礦石礦物主要有錫石、黃銅礦、閃鋅礦、毒砂、黃鐵礦、輝鉍礦等。錫石以自形?半自形粒狀結(jié)構(gòu)零星分布, 粒徑一般小于0.4 mm, 多具環(huán)帶構(gòu)造, 顏色為黃褐色, 與石英、黃玉等緊密共生。黃銅礦、閃鋅礦呈粒狀, 一般不能分辨單晶顆粒大小, 兩者常共生發(fā)育, 交代毒砂、黃鐵礦。黃銅礦顆粒邊緣見黝銅礦; 閃鋅礦粒內(nèi)見黝銅礦、黃錫礦嵌布。

    2.2 分析方法

    2.2.1 鋯石LA-ICP-MS U-Pb測年

    將較新鮮花崗斑巖樣品MD16-13破碎, 經(jīng)人工淘洗、分選后, 在雙目鏡下挑選出自形程度較好的鋯石顆粒, 用環(huán)氧樹脂固定、制靶, 打磨鋯石靶至中心部位出露; 再進行透射光、反射光、陰極發(fā)光照相, 觀察鋯石顆粒結(jié)構(gòu)選擇合適的測年點位。

    鋯石LA-ICP-MS U-Pb定年工作在中國地質(zhì)科學院礦產(chǎn)資源研究所LA-MC-ICP-MS實驗室完成, 采用單點激光剝蝕的方式, 測試儀器為Finnigan Neptune型MC-ICP-MS及配套的Newwave UP 213激光剝蝕系統(tǒng)。數(shù)據(jù)分析前用鋯石GJ-1調(diào)試儀器, 使之達到最優(yōu)狀態(tài), U、Th含量以鋯石M127(Nasdala et al., 2008)為外標進行校正, 測試流程見侯可軍等(2009), 采用ICPMSDataCal 4.3程序(Liu et al., 2008)處理測得數(shù)據(jù), 運用Isoplot3.0程序(Ludwig, 2003)計算、繪制鋯石U-Pb年齡諧和圖。

    (a) 花崗斑巖標本; (b)輝石假象+角閃石反應邊(單偏光); (c) 石英斑晶+鉀長石斑晶+微晶長英質(zhì)基質(zhì)(正交偏光); (d) 鉀長石斑晶+斜長石斑晶+微晶長英質(zhì)基質(zhì)(正交偏光)。

    2.2.2 錫石LA-ICP-MS U-Pb測年

    采用常規(guī)重選法從錫石?硫化物礦石樣品MDX-1中分選出重砂礦物, 在雙目鏡下挑出晶形好、顆粒大的錫石單礦物制靶, 和上述鋯石制靶過程類似。在透射光、反射光和背散射顯微圖像中挑選晶形完整、裂隙少的錫石單礦物供定年測試。

    錫石LA-MC-ICP-MS U-Pb同位素年代學測試分析在中國地質(zhì)調(diào)查局天津地質(zhì)礦產(chǎn)研究所Thermo fisher Neptune型多接收電感耦合等離子體質(zhì)譜儀和ESI UP193-FXArF激光取樣系統(tǒng)完成, 測試的激光波長為193 nm, 脈沖寬度為5 ns, 激光剝蝕束斑為75 μm, 頻率為20 Hz, 能量密度為11 J/cm2。采用AY-4作為標樣(U含量約30×10?6,206Pb/238U表面年齡加權平均值為158.2±0.4 Ma, ID-TIMS; Yuan et al., 2011)。通過測試標樣的Tera-Wasserburg諧和年齡與ID-TIMS標準年齡對比, 計算出校正系數(shù)。具體實驗操作流程可參考崔玉榮等(2017)。研究表明Tera- Wasserburg諧和年齡圖解法可同時完成普通鉛和元素分餾效應的校正, 且相比206Pb/207Pb-238U/207Pb等時線年齡圖解法受高含量的普通Pb影響更小(郝爽等, 2016)。因此本次錫石測年數(shù)據(jù)采用Tera- Wasserburg諧和年齡圖解法計算毛登錫石的形成年齡, 運用Isoplot3.0程序(Ludwig, 2003)進行數(shù)據(jù)處理和圖形繪制。

    2.2.3 全巖主量、微量和稀土元素分析

    選擇MD16-13等8件花崗斑巖樣品清潔、粉碎(200目以下), 在中國核工業(yè)北京地質(zhì)研究院分析測試研究中心完成全巖主量和微量元素分析測試工作。主量元素測試采用X射線熒光光譜法(XRF), 使用儀器為Axiosm AX型X射線熒光光譜儀, 根據(jù)標準巖石樣品監(jiān)控氧化物的分析誤差小于1%。微量和稀土元素測試分析采用電感耦合等離子體質(zhì)譜法(ICP-MS), 使用的儀器為Element XR等離子體質(zhì)譜儀, 根據(jù)標準巖石樣品監(jiān)控分析誤差小于5%。

    2.2.4 全巖Sr-Nd-Pb同位素分析

    花崗斑巖的Sr-Nd-Pb同位素測試分析工作在中國核工業(yè)北京地質(zhì)研究院分析測試研究中心完成, 測試儀器為Phoenix型IsoProbe-T熱表面電離質(zhì)譜儀。Sr、Nd同位素比值分別采用86Sr/88Sr=0.1194、146Nd/144Nd=0.7219進行標準化校正; Pb同位素比值采用NBS-981進行標準化校正。Sr、Nd、Pb同位素測試分析工作具體的實驗室操作流程見趙海杰等(2010)。

    3 分析結(jié)果

    3.1 鋯石U-Pb年齡結(jié)果

    花崗斑巖(MD16-13)中鋯石陰極發(fā)光圖像(圖6a)顯示所選鋯石顆粒多為自形結(jié)構(gòu), 晶形較為完整, 多為長柱狀或短柱狀, 長約80~200 μm, 寬約30~90 μm; 無色透明鋯石顆粒居多, 少數(shù)U含量高的晶體呈暗色; 大部分具有指示巖漿成因的震蕩環(huán)帶。

    利用LA-ICP-MS分析了樣品中的21顆鋯石, Th含量為37×10?6~491×10?6, U含量為72×10?6~1400×10?6, Th/U值為0.34~0.64, 屬巖漿成因。21個分析點的206Pb/238U加權平均年齡為138±0.6 Ma(MSWD=1.02,=21)(表1, 圖7a), 代表了花崗斑巖結(jié)晶年齡。

    3.2 錫石U-Pb測年結(jié)果

    BSE圖像(圖6b)顯示錫石顆粒較大(0.2~0.4 mm), 主要為半自形長柱狀, 次為短柱狀。對礦石樣品MDX-1中錫石顆粒開展38個點位LA-MC-ICP-MS U-Pb測試, 結(jié)果見表2和圖7b。

    測得238U/206Pb值19.36~45.10,207Pb/206Pb值0.0766~0.4734。U-Pb數(shù)據(jù)在Tera-Wasserburg諧和年齡圖上構(gòu)成一條線性數(shù)組, 表示測年錫石含有一致的普通Pb成分。在Tera-Wasserburg諧和年齡圖獲得下交點的年齡為139±3.2 Ma(MSWD=2.3)。張東亮等(2011)研究認為有效擴散半徑為1 μm和1 mm的錫石, 其Pb封閉溫度分別為560 ℃和860 ℃, 本次測年的錫石顆粒多在0.2~0.4 mm之間, U-Pb封閉溫度大于560 ℃。毛登Sn-Cu礦床錫石礦化階段溫度為310~450 ℃(劉玉強, 1996a), 小于錫石的U-Pb體系封閉溫度, 因此認為本次錫石年齡測試結(jié)果139±3.2 Ma能代表錫石結(jié)晶年齡。

    3.3 全巖主量、微量和稀土元素分析結(jié)果

    毛登Sn-Cu礦床8件花崗斑巖樣品的主量、微量和稀土元素分析結(jié)果見表3。

    3.3.1 主量元素

    花崗斑巖具高硅(SiO2=73.18%~77.3%); 貧鈣、鎂、鐵、鋁(CaO=0.49%~0.97%, MgO=0.19%~0.46%, FeO=0.25%~0.92%, Fe2O3=0.44%~1.95%, Al2O3= 11.67%~12.83%); 富堿(Na2O+K2O=8.00%~8.69%, 其中Na2O=3.15%~3.51%、K2O=4.80%~5.22%), Na2O/ K2O=0.6~0.7, 低于國內(nèi)及全球同類型巖石平均值(分別為0.89和1.16; 黎彤等, 1998), 相對富K貧Na; 較高的FeOT/(MgO+FeOT)值(0.74~0.91); 分異指數(shù)DI值為93.0~96.7, 固結(jié)指數(shù)SI值為1.34~4.38, 指示巖漿上升侵位過程經(jīng)歷了明顯的結(jié)晶分異作用。

    圖6 毛登Sn-Cu礦床花崗斑巖中鋯石陰極發(fā)光圖像(a)和錫石BSE圖像(b)

    表1 毛登Sn-Cu礦床花崗斑巖(MD16-13)的LA-ICP-MS鋯石U-Pb分析數(shù)據(jù)

    圖7 毛登Sn-Cu礦床花崗斑巖中鋯石U-Pb年齡諧和圖(a)、錫石U-Pb年齡諧和圖(Tera-Wasserburg)(b)

    表2 毛登Sn-Cu礦床錫石LA-MC-ICP-MS U-Pb定年數(shù)據(jù)

    表3 毛登Sn-Cu礦床花崗斑巖主量(%)、微量和稀土元素(×10?6)組成

    續(xù)表3:

    續(xù)表3:

    注: A/CNK=Al2O3/(Na2O+K2O+CaO)(摩爾比); A/NK=Al2O3/(Na2O+K2O)(摩爾比);Zr.為鋯石飽和溫度, 據(jù)Watson and Harrison (1983)方法計算;1=4Si?11(Na+K)?2(Fe+Ti),2=6Ca+Mg+Al, 據(jù)Batchelor and Bowden (1985)方法計算; TE1, 3為稀土四分組效應指數(shù), 據(jù)Irber (1999)方法計算; DI. 分異指數(shù); SI. 固結(jié)指數(shù)。

    在TAS巖石分類圖解(圖8a)中, 樣品投點落入到花崗巖區(qū)域中, 在QAP巖石分類圖解(圖8b)中, 樣品投點落入到正長花崗巖石和堿長花崗巖區(qū)域中, 與野外觀察命名、巖礦鑒定結(jié)果一致。樣品的里特曼指數(shù)(σ)值在1.91~2.49之間, 屬鈣堿性; 在SiO2-K2O圖解(圖9a)中樣品表現(xiàn)出高鉀鈣堿性巖石的特征。在A/NK-A/CNK圖解(圖9b)中, A/CNK值為0.99~1.04, A/NK值在1.11~1.14, 投點位于準鋁質(zhì)向弱過鋁質(zhì)過渡的區(qū)域。

    在Harker圖解(圖10)中, SiO2與主量元素具有較好的相關性, Al2O3、CaO、Fe2O3T、MgO、MnO、Na2O、P2O5、TiO2的含量變化趨勢相同, 均隨SiO2含量的增加而減少, 表明花崗斑巖經(jīng)歷了分異結(jié)晶演化作用, 在演化過程中發(fā)生鐵鎂質(zhì)礦物、磷灰石、斜長石、鈦氧化物、鐵氧化物的結(jié)晶分離作用, 巖漿總體向富硅, 貧鋁、鈣、鐵、鎂、錳、鉀、鈉、磷、鈦的方向演化。

    3.3.2 稀土元素

    花崗斑巖樣品稀土總含量高, 變化范圍小(ΣREE=319.58×10?6~381.81×10?6), 相對富集輕稀土元素(LREE=291.92×10?6~354.14×10?6), 貧重稀土元素(HREE=25.41×10?6~38.29×10?6)。LREE/HREE值為8.57~12.8, (La/Yb)N值為10.6~18.3, 輕重稀土元素分異明顯。(La/Sm)N值為3.12~4.32, (Gd/Yb)N值為2.09~2.47, 表明輕稀土元素組、重稀土元素組均發(fā)生輕度分餾。稀土元素配分模式圖呈“海鷗式”分布(圖11a), 整體呈右傾, 較陡。巖石具有強烈的Eu負異常(δEu=0.05~0.2), 可能是因為斜長石的分離結(jié)晶作用或者巖漿源區(qū)熔融巖石中大量的斜長石存在; 無Ce異?;蛉魿e正異常(δCe=0.98~1.17), 可能因為巖石在弱酸性、弱氧化條件形成或經(jīng)歷了輕度風化、蝕變作用。

    圖8 毛登Sn-Cu礦花崗斑巖TAS(a, 據(jù)Middlemost, 1994; Ir界線據(jù)Irvine and Baragar, 1971)和QAP圖解(b, 據(jù)Streckeisen, 1973)

    圖9 毛登Sn-Cu礦花崗斑巖K2O-SiO2(a, 據(jù)Peccerillo and Taylor, 1976)及A/NK-A/CNK(b, 據(jù)Maniar and Piccoli, 1989; 虛線代表I型和S型花崗巖之間的邊界, 據(jù)Chappell and White, 1992)圖解

    3.3.3 微量元素

    微量元素測試結(jié)果顯示, 樣品總體上呈相對富集大離子親石元素LILE, 虧損高場強元素HFSE的特點; (Zr+Nb+Ce+Y)值(396~595)×10?6, 平均為480×10?6, 明顯高于A型花崗巖的下限值350×10?6; 10000Ga/Al值(3.64~4.28)平均為3.89, 大于A型花崗巖與其他類型花崗巖的分界值2.6(Whalen et al., 1987; 吳福元等, 2007), 顯示A型花崗巖地球化學特征。原始地幔標準化微量元素蜘蛛網(wǎng)圖(圖11b)呈現(xiàn)顯著的Ba、Sr、P、Ti、Nb、Ta等元素負異常和Rb、Th、La、Nd、Zr等元素的正異常, 顯示出大陸弧環(huán)境造山帶花崗巖微量元素特征(李昌年, 1992)。Ba、Sr、P、Ti的負異常暗示了長石、磷灰石和鈦鐵礦的結(jié)晶分離或在源區(qū)部分熔融過程中的大量殘留。毛登Sn-Cu礦床花崗斑巖整體上與大興安嶺南段錫多金屬成礦有關的燕山晚期花崗巖有相似的主、微量地球化學特征(郝立波等, 1999)。

    3.4 全巖Sr-Nd-Pb同位素分析結(jié)果

    Sr、Nd和Pb同位素研究在巖漿巖物質(zhì)來源方面有著重要的示蹤作用,Nd()、208Pb/204Pb、207Pb/204Pb、206Pb/204Pb等值可以較為準確的指示成巖物質(zhì)來源。

    3.4.1 Sr-Nd同位素

    巖石樣品Sr、Nd測試結(jié)果見表4。根據(jù)本次獲得的花崗斑巖鋯石U-Pb年齡(138 Ma)計算得出初始(87Sr/86Sr)i、Nd()。此次毛登Sn-Cu礦床花崗斑巖測得較高的Rb含量(183~226)×10?6以及較低的Sr含量(21~54)×10?6,87Rb/86Sr值(9.82~30.11)較高且范圍較大,87Sr/86Sr值(0.713429~0.780369)、(87Sr/86Sr)i(0.692360~0.721307)變化范圍較大。前人研究大興安嶺南段燕山期晚期花崗巖(87Sr/86Sr)i值多介于0.7020~0.7100之間(張德全, 1993; 劉偉等, 2007; 周振華等, 2011; 顧玉超等, 2017), 剔除毛登花崗斑巖4個(87Sr/86Sr)i值不在0.7020~0.7100之間的不合理數(shù)據(jù), 剩余樣品的(87Sr/86Sr)i值在0.702907~0.704506之間。(143Nd/144Nd)i值為0.512503~0.512686,Nd()值為+0.8~+4.4, 具有新生地殼源區(qū)特征;Sm/Nd為?0.54~?0.31(平均值?0.45), 基本介于?0.5~?0.2之間, Nd模式年齡有實際地質(zhì)意義(吳福元等, 1997), 計算獲得Nd同位素兩階段模式年齡DM2為572~863 Ma。

    3.4.2 Pb同位素

    花崗斑巖的Pb同位素測試結(jié)果見表5。樣品的208Pb/204Pb、207Pb/204Pb、206Pb/204Pb值都較高, 分別為38.572~38.831、15.572~15.601、18.642~18.851, 其初始鉛同位素比值分別為(208Pb/204Pb)i=38.222~ 38.391, (207Pb/204Pb)i=15.563~15.585, (206Pb/204Pb)i= 18.441~18.526。李睿華(2019)對毛登Sn-Cu礦床8件礦石樣品Pb同位素測試結(jié)果顯示, (208Pb/204Pb)i、(207Pb/204Pb)i、(206Pb/204Pb)i值分別為37.668~38.622、15.482~15.679、17.772~18.427。由圖12和表5可知, 花崗斑巖體的初始鉛Pb同位素比值變化范圍小, 具有相似的Pb同位素組成, 暗示Pb來源相對簡單, 而礦石Pb同位素比值變化范圍較大, 具有較好的線性分布特征, 表明為混合來源。

    在207Pb/204Pb-206Pb/204Pb圖解(圖12a)中, 花崗斑巖樣品初始Pb同位素比值均落入造山帶和地幔演化曲線之間, 且靠近造山帶區(qū)域, 具有從造山帶?地幔過渡的特征; 礦石樣品Pb同位素比值分布較分散, 多數(shù)落入造山帶和地幔演化曲線之間, 1件樣品落入上地殼演化線附近, 1件樣品緊靠地幔演化線附近。在208Pb/204Pb-206Pb/204Pb圖解(圖12b)中, 花崗斑巖樣品均落入造山帶演化曲線附近, 指示Pb可能來源于造山帶; 礦石樣品均落入造山帶演化曲線上方, 指示Pb主要來源于造山帶, 混有一定量的下地殼物質(zhì)。造山帶本身代表了殼?;旌蟻碓? 因此從鉛同位素構(gòu)造模式圖解可推斷花崗斑巖體、礦石Pb具深源特征, 可能來源于殼?;旌香U。另外, 花崗斑巖Pb同位素投點非常靠近礦石投點區(qū)域, 指示花崗斑巖體與礦石Pb的親緣性, 表明礦石Pb與巖體關系密切。

    圖10 毛登Sn-Cu礦床花崗斑巖哈克圖解

    圖11 毛登Sn-Cu礦床花崗斑巖樣品球粒隕石標準化稀土元素配分圖(a)和原始地幔標準化微量元素蜘蛛網(wǎng)圖(b)(標準化值據(jù)Sun and McDonough, 1989)

    表4 毛登Sn-Cu礦床花崗斑巖Sr-Nd同位素組成

    表5 毛登Sn-Cu礦床花崗斑巖、礦石Pb同位素組成

    注: a. 本文; b. 李睿華(2019)。

    4 討 論

    4.1 成巖(礦)時代

    毛登Sn-Cu礦床處于大興安嶺南段Pb-Zn-Ag- Cu-Sn-Fe成礦帶(劉建明等, 2004; 要梅娟等, 2012), 前人對區(qū)內(nèi)多個典型礦床的成巖(礦)年代學已開展了研究, 毛景文等(2005, 2013)認為大興安嶺南段地區(qū)大量金屬礦床聚集與區(qū)內(nèi)中生代花崗質(zhì)巖漿侵位關系密切; 周振華等(2010)、翟德高等(2012)等研究認為燕山晚期(130~140 Ma)是大興安嶺南段地區(qū)成巖(礦)的高峰期之一; 王長明等(2006)研究大興安嶺南段銅多金屬礦成礦時空結(jié)構(gòu)表明區(qū)內(nèi)錫成礦作用主要集中在燕山晚期。本文在前人統(tǒng)計的大興安嶺南段主要錫多金屬礦床的成礦年齡及與成礦關系密切巖體侵入年齡基礎上(劉瑞麟等, 2018), 補充了維拉斯托Sn-Cu-Pb-Zn-Ag多金屬礦床錫石U-Pb年齡136 Ma, 輝鉬礦Re-Os年齡136.8~138.8 Ma, 堿長花崗巖鋯石U-Pb年齡137~138 Ma(劉瑞麟, 2018); 拜仁達壩銀鉛鋅錫多金屬礦床礦石中蝕變礦物白云母Ar-Ar年齡135 Ma(常勇和賴勇, 2010); 毛登?小狐山錫多金屬礦田花崗斑巖鋯石U-Pb年齡132~133 Ma, 似斑狀花崗巖鋯石U-Pb年齡134 Ma, 錫石U-Pb年齡134.8 Ma、140 Ma(郭碩等, 2019; 李睿華, 2019); 查木罕鎢鉬錫礦床輝鉬礦Re-Os年齡140 Ma, 二長花崗巖獨居石U-Pb年齡137 Ma、138 Ma、139 Ma (Zhang et al., 2019); 白音查干錫銀多金屬礦床電氣石花崗斑巖鋯石U-Pb年齡為140.5 Ma, 花崗斑巖鋯石U-Pb年齡143.4 Ma, 錫石U-Pb年齡140 Ma、140.6 Ma(李睿華, 2019); 道倫達壩銅鎢錫銀礦床張家營子細?;◢弾r鋯石U-Pb年齡135 Ma, 礦石錫石U-Pb年齡134.7 Ma、136.8 Ma(陳公正等, 2018), 統(tǒng)計結(jié)果見圖13。從圖13可以看出大興安嶺南段錫多金屬礦床形成于126~144 Ma, 成礦時代峰值135~ 139 Ma, 與成礦關系密切的花崗質(zhì)巖體侵入時代為129~149 Ma, 成巖時代峰值137~140 Ma。綜上得出, 大興安嶺南段Sn多金屬礦床成巖(礦)時代主要集中在燕山晚期的早白堊世(135~140 Ma), 成巖時代范圍略大于成礦時代, 兩者具有一致性, Sn多金屬礦化與早白堊世花崗質(zhì)巖漿有緊密的聯(lián)系; 成礦時代晚于成巖時代, 指示花崗質(zhì)巖體成礦熱事件可能具有多期次、多階段的演化, 與區(qū)域Sn多金屬成礦系列先后發(fā)生高溫浸染狀Sn-W-Mo和云英巖化脈狀Sn-W氧化物礦化及中低溫Cu-Pb-Zn和Pb-Zn-Ag-Sb硫化物礦化也較為吻合。

    本次獲得毛登Sn-Cu礦床花崗斑巖LA-ICP-MS鋯石U-Pb年齡為138±0.6 Ma, 代表了巖漿巖的形成年齡; 錫石LA-ICP-MS U-Pb年齡為139±3.2 Ma, 代表了錫石的結(jié)晶年齡, 成巖成礦時間一致為早白堊世。毛登Sn-Cu礦床成巖(礦)年齡測試結(jié)果與上述大興安嶺南段主要錫多金屬礦床成巖(礦)年齡基本一致, 說明毛登Sn-Cu礦床與區(qū)內(nèi)錫多金屬成礦系列皆屬于大興安嶺南段燕山晚期巖漿活動的產(chǎn)物。

    圖12 毛登Sn-Cu礦床Pb同位素構(gòu)造模式圖(據(jù)Zartman and Doe, 1981)

    數(shù)據(jù)來源于: 常勇和賴勇, 2010; 陳公正等, 2018; 劉瑞麟, 2018; 劉瑞麟等, 2018; Zhang et al., 2019; 郭碩等, 2019; 李睿華, 2019。

    4.2 花崗巖成因

    4.2.1 巖石類型

    前人對大興安嶺南段廣泛發(fā)育的晚侏羅世?早白堊世花崗巖成因類型有不同的認識, 多數(shù)屬于A型花崗巖(邵濟安等, 2010), 包括白音查干礦床、巴爾哲地區(qū)、邊家大院礦床(A1型)、黃崗礦床(A1型)、德勒哈達地區(qū)(A2型)、石匠山地區(qū)(A2型)、努和特地區(qū)(A2型)等(周振華等, 2010; 楊武斌等, 2011; 王喜龍等, 2014; 王金芳等, 2017; 姚磊等, 2017; 王金芳等, 2018a, 2018b); 少數(shù)屬于I型花崗巖, 如維拉斯托礦床石英斑巖(Wang et al., 2017; 張?zhí)旄5? 2019), 白音查干花崗斑巖(李睿華, 2019)。另外楊奇荻等(2014)認為甘珠爾廟地區(qū)早白堊世花崗巖為I-A過渡型。雖然對花崗巖類型觀點不完全一致, 但研究者一致認為大興安嶺南段晚侏羅世?早白堊世花崗巖經(jīng)歷了高分異演化作用。高分異演化的I型、A型花崗巖主量元素以及礦物相趨于共熔點, 具有相似的地球化學、礦物學特征(King et al., 1997, 2001)。經(jīng)歷了高度演化結(jié)晶分異作用的殘余熔漿固相線溫度大大降低, 在最終的巖漿?熱液過渡體系中發(fā)生強烈的熔體?流體作用, 致使巖體中主(微)量元素地球化學行為發(fā)生了很大變化(吳福元等, 2007)。以上兩方面原因給大興安嶺南段燕山晚期高分異花崗巖成因類型判別造成了一定的困擾。

    毛登Sn-Cu礦床8件花崗斑巖樣品具有較高的分異指數(shù)(DI值為93.0~96.7), 固結(jié)指數(shù)較小(SI值為1.34~4.38), Ba、Sr、P、Ti、Eu顯示負異常(圖11)。微量元素Zr/Hf值可反映花崗巖結(jié)晶分異演化程度和熔體?流體相互作用情況, 隨著巖漿的演化Zr/Hf值逐漸減小(Irber, 1999; Dostal and Chatterjee, 2000)。毛登礦床花崗斑巖Zr/Hf值為25.8~33.0 (平均值30.0), 一般花崗巖Zr/Hf值為39(章榮清等, 2016)。主(微)量元素組成均表明毛登礦床花崗斑巖具有高分異演化花崗巖的特征, 表明了成巖過程經(jīng)歷了較強烈的結(jié)晶分異作用, 常用的花崗巖類型判別圖解法很可能對毛登礦區(qū)花崗斑巖失效。

    吳福元等(2017)在Whalen et al. (1987)研究A型花崗巖地球化學判別基礎上, 提出A型花崗巖在結(jié)晶分異過程中10000Ga/Al值逐漸減小, I或S型花崗巖則逐漸升高, 對區(qū)分A型花崗巖和高分異花崗巖效果較好。將毛登、白音查干、維拉斯托和黃崗等大興安嶺南段四個典型的Sn多金屬礦床的花崗質(zhì)巖體投點在A型花崗巖與高分異花崗巖判別圖解(圖14)中, 可以看出毛登、黃崗花崗巖體部分投點落入A型和I/S分異型的相交區(qū)域, 總體呈A型花崗巖的演化趨勢; 白音查干花崗巖體多數(shù)投點落入A型和I/S分異型的相交區(qū)域, 總體呈I/S分異型花崗巖演化趨勢; 維拉斯托花崗巖體投點均落入I/S分異型區(qū)域。

    黃崗、維拉斯托、白音查干巖體數(shù)據(jù)來源于周振華等, 2010; 趙輝等, 2015; 祝新友等, 2016; Wang et al., 2017; 劉新等, 2017; 姚磊等, 2017; 劉瑞麟, 2018; 張?zhí)旄5? 2019; 郭碩等, 2019; 李睿華, 2019。I、S、M和A分別代表I型、S型、M型和A型花崗巖。

    A型花崗巖形成的溫度(800~900 ℃)明顯高于I型和S型花崗巖(劉昌實等, 2003; 汪洋, 2008, 2009; 賈小輝等, 2009; 王楠等, 2017)。研究認為鋯石的飽和溫度近似代表花崗質(zhì)巖石近液相線的溫度, 可用來估算花崗質(zhì)巖漿的結(jié)晶溫度(King et al., 1997; 王濤和劉燊, 2013)。毛登Sn-Cu礦床花崗斑巖樣品的Zr含量隨著SiO2含量的增加而減少, 表明巖漿中Zr達到了飽和(圖15)。利用花崗質(zhì)熔體鋯石飽和溫度公式(Watson and Harrison, 1983), 計算得出毛登Sn-Cu礦床花崗斑巖鋯石飽和溫度為817~861 ℃, 平均834 ℃(表3), 屬高溫花崗巖, 與典型的A型花崗巖形成溫度一致。

    綜上認為毛登、黃崗礦床花崗斑巖為A型花崗巖, 成巖演化過程中經(jīng)歷了較強烈的結(jié)晶分異作用; 白音查干、維拉斯托礦床花崗質(zhì)巖體為高分異的I型花崗巖。燕山晚期造山帶垮塌和地殼逐漸減薄誘導幔源物質(zhì)底侵以及地殼減薄不均一性, 是大興南嶺南段早白堊世形成A型、I型不同類型花崗巖的主要原因(趙輝等, 2015)。

    4.2.2 巖石源區(qū)

    關于A型花崗巖的成因主要有以下3種觀點: ①地幔玄武質(zhì)巖漿的分離結(jié)晶作用(Han et al., 1997; Mushkin et al., 2003); ②幔源巖漿與地殼深熔巖漿的混合作用(Dickin, 1994; Harris et al., 1999; Mingram et al., 2000); ③下地殼源巖的部分熔融(Collins et al., 1982; Clemens et al., 1986)。Wyllie (1984)研究認為地幔物質(zhì)部分熔融形成的大量玄武質(zhì)巖漿才可能分異出少量的花崗質(zhì)巖漿, 研究區(qū)內(nèi)沒發(fā)現(xiàn)與花崗巖時空關系密切且大面積分布的玄武質(zhì)巖體出露, 因此花崗斑巖由地幔玄武質(zhì)巖漿的分離結(jié)晶作用形成的可能性不大。毛登花崗斑巖強烈虧損Eu、Ba、Sr, 表明其源區(qū)有大量長石殘留, 暗示成巖在低壓、較淺的環(huán)境(Rapp and Watson, 1995), 同樣排除為幔源巖漿分離結(jié)晶作用的成因。毛登花崗斑巖未見明顯的暗色基性包體, 不顯示幔源、殼源巖漿混合特征, 排除其為幔源巖漿與地殼深熔巖漿混合作用的成因。

    圖15 毛登Sn-Cu礦床花崗斑巖Zr-SiO2圖解

    Jahn et al. (2000)研究表明中亞造山帶是地球上顯生宙新生地殼增長最快的地區(qū), 位于其東部的大興安嶺地區(qū)顯生宙的花崗巖表現(xiàn)出正的Nd()值(平均值+2.0)和較低的DM值(多在500~1000 Ma), 暗示花崗巖源區(qū)為幔源年輕的地殼物質(zhì)(吳福元等, 1997; 洪大衛(wèi)等, 2000; Wu et al., 2000, 2002)。肖成東等(2004)測得東蒙地區(qū)燕山期花崗巖的Nd()值為+0.75~+8.12(均值+3.07),DM2為2622~930 Ma(均值705 Ma), 認為成巖物質(zhì)來源與虧損地幔有關。毛登Sn-Cu礦床花崗斑巖Sr-Nd同位素測試分析顯示其具有低的(87Sr/86Sr)i值、正的Nd()值(Nd()=+0.8~ +4.4)和年輕的Nd同位素模式年齡(DM2=572~863 Ma),與中亞造山帶顯生宙花崗巖Sr-Nd同位素特征一致(洪大衛(wèi)等, 2000; Jahn et al., 2000, 2001; Wu et al., 2003;劉偉等, 2007; 邵濟安等, 2010)。在(87Sr/86Sr)i-Nd()圖解(圖16)中, 樣品投點落入或靠近主地幔趨勢線上, 靠近虧損地幔單元, 遠離華北陸塊下地殼區(qū)域, 與大興安嶺中南段晚中生代花崗巖Sr-Nd同位素范圍(Jahn et al., 1999)基本一致, 說明成巖物質(zhì)來源于地幔。Sr-Nd同位素特征表明了毛登礦區(qū)花崗斑巖有可能來源于新元古代含有大量幔源組分的新生地殼物質(zhì), 在巖石圈伸展減薄背景下幔源玄武質(zhì)巖漿底侵作用誘使新生地殼部分熔融, 大量巖漿上升侵位并經(jīng)歷了較強烈的結(jié)晶分異作用, 形成酸性侵入巖。

    4.3 構(gòu)造背景

    關于大興安嶺中生代巖漿活動構(gòu)造背景存在著不同的觀點: 趙國龍等(1989)認為大興安嶺中生代火山巖形成于與古太平洋板塊俯沖有關的活動大陸邊緣; 張連昌等(2007)認為形成于蒙古?鄂霍茨克洋造山后伸展和古太平洋板塊俯沖的聯(lián)合作用, 反映陸內(nèi)伸展環(huán)境和陸緣弧火山巖雙重特征; 吳福元和曹林(1999)指出中生代以來東部太平洋板塊的俯沖至關重要。近來年, 越來越多的學者認為大興安嶺地區(qū)大規(guī)模的中生代巖漿活動構(gòu)造背景歸于蒙古?鄂霍茨克洋閉合后伸展構(gòu)造背景(郭鋒等, 2001; Fan et al., 2003; 呂志成等, 2004; 江思宏等, 2012; Ouyang et al., 2015), 或古太平洋板塊俯沖后撤伸展構(gòu)造背景環(huán)境(王喜龍等, 2014; 王良玉等, 2016)。

    在A型花崗巖分類判別圖解(圖17)中, 花崗斑巖樣品投點大多數(shù)落入了A2型區(qū), 代表造山后環(huán)境(如陸陸碰撞后期或島弧巖漿作用環(huán)境)。在Al2O3- SiO2、FeOT/(FeOT+MgO)-SiO2、(Y+Nb)-Rb和1-2構(gòu)造環(huán)境判別圖解(圖18、19)中, 大多數(shù)樣品落在造山期后A型花崗巖區(qū)域, 暗示其具有拉張構(gòu)造環(huán)境特征。

    大興安嶺中南段數(shù)據(jù)來源: 張德全等, 1993; 趙一鳴等, 1994; 蔡劍輝等, 2004; 肖成東等, 2004; 周振華等, 2011; 顧玉超等, 2017。

    結(jié)合區(qū)域構(gòu)造演化史(劉偉等, 2007; 周振華等, 2011; 陳公正等, 2018; 江思宏等, 2018)以及毛登Sn-Cu礦床花崗斑巖地球化學特征、成巖年代學、Sr-Nd-Pb同位素特征, 本文認為在早白堊世大興安嶺南段毛登地區(qū)可能受到蒙古?鄂霍茨克洋造山后伸展和古太平洋俯沖的雙重影響, 蒙古?鄂霍茨克洋構(gòu)造域已進入后碰撞階段, 古太平洋板塊俯沖方向由N-NNW向NW轉(zhuǎn)變(Engebretson et al., 1985), 兩大構(gòu)造域都處于構(gòu)造體系轉(zhuǎn)變期間, 區(qū)域構(gòu)造體系由擠壓環(huán)境向伸展環(huán)境轉(zhuǎn)換過渡期, 總體處于伸展構(gòu)造背景。兩大構(gòu)造體系聯(lián)合作用也被認為是大興安嶺南段地區(qū)燕山晚期成礦高峰期的主要原因(Ouyang et al., 2015; 江思宏等, 2018)。

    4.4 巖漿作用對成礦制約

    毛登Sn-Cu礦床主要為受斷裂?節(jié)理裂隙控制的石英脈型礦床, 花崗斑巖體與角礫巖等圍巖的接觸帶外側(cè)控制礦體的產(chǎn)出形態(tài), 巖體內(nèi)部礦化弱。錫石礦化和花崗斑巖形成的年齡分別為139±2.8 Ma和138±0.6 Ma, 兩者在誤差范圍內(nèi)一致, 表明毛登Sn-Cu礦床的形成與該期巖漿活動有關。李睿華(2019)研究毛登Sn-Cu礦床8件金屬硫化物S、Pb同位素特征, 發(fā)現(xiàn)δ34SV-CDT值較均一, 介于?3.3‰~ ?0.6‰之間(平均值?2.1‰), 具有巖漿成因的S同位素特征; Pb同位素組成與本文研究的花崗斑巖范圍接近(圖12), 認為金屬硫化物成礦物質(zhì)主要來源于花崗斑巖體。劉玉強(1996a)對各成礦階段流體的H-O同位素研究顯示成礦流體來自花崗斑巖的深部, 晚期蝕變和礦化流體中大氣降水逐漸增多。礦床內(nèi)黃玉、螢石等蝕變礦物共生表明, 成礦熱液具有巖漿成因的高F特征(Wang et al., 2017)。

    A1. 非造山花崗巖; A2. 造山后花崗巖。

    IAG. 島弧花崗巖類; CAG. 大陸弧花崗巖類; CCG. 大陸碰撞花崗巖類; POG. 后造山花崗巖類; RRG. 與裂谷有關的花崗巖類; CEUG. 與大陸的造陸抬升有關的花崗巖類。

    圖19 毛登Sn-Cu礦床花崗斑巖(Y+Nb)-Rb(a, Pearce et al., 1984)和R1-R2(b, 據(jù)Batchelor and Bowden, 1985)構(gòu)造背景判別圖解

    前人研究認為, 花崗質(zhì)巖漿經(jīng)歷高度結(jié)晶分異作用后, 在演化晚期的殘余熔體?熱液共存體系中富含揮發(fā)分、Na2O+K2O、REE、稀有金屬, 熔體與熱液的相互作用使得Eu的負異常增強, 以及W、Sn、Sb、REE、Nb、Ta、Zr、Be等成礦元素在殘余花崗質(zhì)巖漿中逐漸富集, 為形成W、Sn、稀土、稀有礦床提供物質(zhì)基礎(Jahn et al., 2001; Zhao et al., 2002; Wu et al., 2004)。例如華南地區(qū)西華山W-Sn礦床、千里山W-Sn礦床、喜馬拉雅造山帶扎西康Sb-Pb-Zn-Au礦床、大興安嶺南段地區(qū)維拉斯托Sn多金屬礦床、沙麥W-Mo礦床、巴爾哲REE-Nb-Zr-Be礦床(Zhao et al., 2002; Guo et al., 2012; 梁維等, 2015; Jiang et al., 2016; Wang et al., 2017)。

    前述全巖主、微量地球化學特征分析表明, 毛登礦區(qū)花崗斑巖為經(jīng)歷了較強烈的結(jié)晶分異作用的A型花崗巖, 李小偉等(2010)認為A型花崗巖的Rb含量大于(250×10?6~300×10?6)時屬于高分異型花崗巖, 毛登礦床花崗斑巖Rb含量為183×10?6~226×10?6, 顯然未達到“高”分異演化程度。研究表明高演化結(jié)晶分異的花崗質(zhì)巖石具有M型稀土四分組效應(Bau, 1996, 1997; 趙振華等, 1999), 毛登花崗斑巖稀土元素量化參數(shù)TE1, 3和關鍵元素K/Rb、Y/Ho、Zr/Hf的比值(分別為0.95~1.02, 187~237、22.4~26.9、25.8~33.0)接近或偏離具明顯四分組效應花崗巖的區(qū)間范圍(TE1, 3>1.10, K/Rb<100, Y/Ho>28, Zr/Hf<25) (Irber, 1999)。綜上認為毛登花崗斑巖為分異的A型花崗巖, 在演化過程中經(jīng)歷了較強烈的結(jié)晶分異作用, 對Sn、W等成礦金屬元素的富集起到積極作用; 但巖體并未達到維拉斯托礦床石英斑巖“高”演化結(jié)晶分異的程度(Wang et al., 2017), 這或許是毛登礦床Sn金屬資源量未達到維拉斯托礦床規(guī)模的原因之一。

    5 結(jié) 論

    (1) 毛登Sn-Cu礦床的LA-ICP-MS錫石U-Pb等時線年齡為139±3.2 Ma, 花崗斑巖的LA-ICP-MS鋯石U-Pb年齡為138±0.6 Ma, 二者在誤差范圍內(nèi)近于一致, 表明成巖成礦均發(fā)生在早白堊世。

    (2) 毛登Sn-Cu礦床花崗斑巖的形成受到蒙古?鄂霍茨克洋閉合造山后伸展和古太平洋俯沖的雙重影響, 總體處于伸展構(gòu)造背景; 巖石地球化學和Sr-Nd同位素特征表明, 毛登花崗斑巖具備A型花崗巖特點, 為新元古代新生地殼部分熔融的產(chǎn)物; Pb同位素組成特征表明, 成礦物質(zhì)主要來自深部巖漿。

    (3) 毛登Sn-Cu礦床花崗斑巖經(jīng)歷了較強烈的結(jié)晶分異演化作用, 但未達到“高”分異演化程度, 這可能是毛登礦床Sn金屬資源量未達到維拉斯托礦床規(guī)模的原因之一。

    致謝:野外地質(zhì)工作期間得到了錫林浩特市興原有色金屬有限責任公司杜祖權總經(jīng)理和史賀龍工程師的大力協(xié)助, 吉林大學孫景貴教授、中國地質(zhì)科學院礦產(chǎn)資源研究所武廣研究員對本文提出了寶貴的修改意見, 在此一并表示衷心的感謝!

    蔡劍輝, 閻國翰, 肖成東, 王關玉, 牟保磊, 張任祜. 2004. 太行山?大興安嶺構(gòu)造巖漿帶中生代侵入巖Nd、Sr、Pb同位素特征及物質(zhì)來源探討. 巖石學報, 20(5): 1225–1242.

    常勇, 賴勇. 2010. 內(nèi)蒙古銀都銀鉛鋅多金屬礦床成礦流體特征及成礦年代學研究. 北京大學學報(自然科學版), 46(4): 581–593.

    陳公正, 武廣, 李鐵剛, 劉瑞麟, 武利文, 章陪春, 張彤, 陳毓川. 2018. 內(nèi)蒙古道倫達壩銅鎢錫礦床LA-ICP- MS鋯石和錫石U-Pb年齡及其地質(zhì)意義. 礦床地質(zhì), 27(2): 225–245.

    陳鄭輝, 王登紅, 盛繼福, 應立娟, 梁停, 王成輝, 劉麗君, 王永幕. 2015. 中國錫礦成礦規(guī)律概要. 地質(zhì)學報, 89(6): 1026–1037.

    程天赦, 楊文靜, 王登紅. 2014. 內(nèi)蒙古西烏旗阿魯包格山A型花崗巖鋯石U-Pb年齡、地球化學特征及地質(zhì)意義. 大地構(gòu)造與成礦學, 38(3): 718–728.

    崔玉榮, 涂家潤, 陳楓, 郝爽, 葉麗娟, 周紅英, 李惠民. 2017. LA-(MC)-ICP-MS錫石U-Pb定年研究進展. 地質(zhì)學報, 91(6): 1386–1399.

    顧玉超, 陳仁義, 賈斌, 宋萬兵, 余昌濤, 鞠楠. 2017. 內(nèi)蒙古邊家大院鉛鋅銀礦床深部正長花崗巖年代學與形成環(huán)境研究. 中國地質(zhì), 44(1): 101–117.

    郭鋒, 范蔚茗, 王岳軍, 林舸. 2001. 大興安嶺南段中生代雙峰式火山作用. 巖石學報, 17(1): 161–168.

    郭碩, 何鵬, 張學斌, 崔玉榮, 張?zhí)旄? 張闊, 來林, 劉傳寶. 2019. 大興安嶺南段毛登?小孤山錫多金屬礦田年代學、地球化學特征及其地質(zhì)意義. 礦床地質(zhì), 38(3): 509–525.

    郝立波, 段國正, 李殿超, 呂志成, 潘軍. 1999. 大興安嶺錫多金屬成礦帶花崗巖地球化學特征. 世界地質(zhì), 18(2): 66–72.

    郝爽, 李惠民, 李國占, 耿建珍, 周紅英, 肖志斌, 崔玉榮, 涂家潤. 2016. LA-ICP-MS測定錫石U-Pb同位素年齡時兩種普通鉛扣除方法的原理及適用性比較. 地質(zhì)通報, 35(4): 622–632.

    河南省有色金屬地質(zhì)勘查總院. 2009. 內(nèi)蒙古自治區(qū)錫林浩特市毛登礦區(qū)鉬錫銅礦詳查報告: 1–179.

    洪大衛(wèi), 王式?jīng)? 謝錫林, 張季生. 2000. 興蒙造山帶正εNd(t)值花崗巖的成因和大陸地殼生長. 地學前緣, 7(2): 441–455.

    侯可軍, 李延河, 田有榮. 2009. LA-MC-ICP-MS鋯石微區(qū)原位U-Pb定年技術. 礦床地質(zhì), 28(4): 481–492.

    賈小輝, 王強, 唐功建. 2009. A型花崗巖的研究進展及意義. 大地構(gòu)造與成礦學, 33(3): 465–480.

    江思宏, 梁清玲, 劉翼飛, 劉妍. 2012. 內(nèi)蒙古大井礦區(qū)及外圍巖漿巖鋯石U-Pb年齡及其對成礦時間的約束. 巖石學報, 28(2): 495–513.

    江思宏, 張莉莉, 劉翼飛, 劉春花, 康歡, 王豐翔. 2018. 興蒙造山帶成礦規(guī)律及若干科學問題. 礦床地質(zhì), 37(4): 671–711.

    黎彤, 袁懷雨, 吳勝昔. 1998. 中國花崗巖類和世界花崗巖類平均化學成分的對比研究. 大地構(gòu)造與成礦學, 22(1): 24–29.

    李昌年. 1992. 火成巖微量元素巖石學. 武漢: 中國地質(zhì)大學出版社: 97–109.

    李睿華. 2019. 大興安嶺南段錫林浩特地區(qū)錫銅多金屬礦床的成礦作用. 北京: 北京大學博士學位論文: 1–240.

    李小偉, 莫宣學, 趙志丹, 朱弟成. 2010. 關于A型花崗巖判別過程中若干問題的討論. 地質(zhì)通報, 29(2–3): 278–285.

    梁維, 楊竹森, 鄭遠川. 2015. 藏南扎西康鉛鋅多金屬礦絹云母Ar-Ar年齡及其成礦意義. 地質(zhì)學報, 89(3): 560–568.

    劉昌實, 陳小明, 陳培榮, 王汝成, 胡歡. 2003. A型巖套的分類、判別標志和成因. 高校地質(zhì)學報, 9(4): 573– 591.

    劉傳寶, 王強, 王建忠, 郎欣欣, 張猛, 侯冠宇. 2018. 內(nèi)蒙古毛登地區(qū)阿魯包格山花崗斑巖地質(zhì)特征及U-Pb年代學意義. 礦產(chǎn)勘查, 9(7): 1316–1321.

    劉建明, 張銳, 張慶洲. 2004. 大興安嶺地區(qū)的區(qū)域成礦特征. 地學前緣, 11(1): 269–273.

    劉瑞麟, 武廣, 李鐵剛, 陳公正, 武利文, 章培春, 張彤, 江彪, 劉文元. 2018. 大興安嶺南段維拉斯托錫多金屬礦床LA-ICP-MS錫石和鋯石U-Pb年齡及其地質(zhì)意義. 地學前緣, 25(5): 1–15.

    劉瑞麟. 2018. 內(nèi)蒙古維拉斯托稀有金屬?錫多金屬成礦系統(tǒng)的成礦作用. 北京: 北京大學博士學位論文: 1–185.

    劉偉, 潘小菲, 謝烈文, 李禾. 2007. 大興安嶺南段林西地區(qū)花崗巖類的源巖: 地殼生長的時代和方式. 巖石學報, 23(2): 441–460.

    劉新, 李學剛, 祝新友, 孫雅琳, 蔣昊原, 蔣斌斌, 王海, 程細音. 2017. 內(nèi)蒙古白音查干錫多金屬礦床成礦作用研究Ⅱ: 成礦花崗斑巖年代學、地球化學特征及地質(zhì)意義. 礦產(chǎn)勘查, 8(6): 981–996.

    劉玉強. 1996a. 內(nèi)蒙古毛登錫銅礦床地質(zhì)及成因. 礦床地質(zhì), 15(2): 133–143.

    劉玉強. 1996b. 毛登錫銅礦床成礦分帶及其成因討論. 礦床地質(zhì), 15(4): 318–329.

    呂志成, 段國正, 郝立波, 李殿超, 潘軍, 吳豐昌. 2004. 大興安嶺中南段中生代中基性火山巖巖石學地球化學研究. 高校地質(zhì)學報, 10(2): 186–198.

    呂志成, 段國正, 劉叢強, 郝立波, 李殿超, 魏存第. 2000. 大興安嶺地區(qū)銀礦床類型、成礦系列及成礦地球化學特征. 礦物巖石地球化學通報, 19(4): 305–309.

    毛景文, 謝桂青, 張作衡, 李曉峰, 王義天, 張長青, 李永峰. 2005. 中國北方中生代大規(guī)模成礦作用的期次及其地球動力學背景. 巖石學報, 21(1): 169–188.

    毛景文, 周振華, 武廣, 江思宏, 劉成林, 李厚民, 歐陽荷根, 劉軍. 2013. 內(nèi)蒙古及鄰區(qū)礦床成礦規(guī)律與成礦系列. 礦床地質(zhì), 32(4): 715–729.

    邵濟安, 牟保磊, 朱慧忠, 張履橋. 2010. 大興安嶺中南段中生代成礦物質(zhì)的深部來源與背景. 巖石學報, 26(3): 649–656.

    石得鳳. 2007. 內(nèi)蒙古錫林浩特市毛登多金屬礦區(qū)地質(zhì)特征及成礦模式研究. 長沙: 中南大學碩士學位論文: 1–72.

    汪洋. 2008. 再論南嶺侏羅紀“鋁質(zhì)”A型花崗巖的成因及其對古低溫線的制約. 大地構(gòu)造與成礦學, 32(3): 365–381.

    汪洋. 2009. 北京白查A型花崗巖的地球化學特征及其成因與構(gòu)造指示意義. 巖石學報, 25(1): 13–24.

    王春女, 王全明, 于曉飛, 韓振哲. 2016. 大興安嶺南段錫礦成礦特征及找礦前景. 地質(zhì)與勘探, 52(2): 220– 227.

    王金芳, 李英杰, 李紅陽, 董培培. 2017. 內(nèi)蒙古西烏旗努和特早白堊世A型花崗巖LA-ICP-MS鋯石U-Pb年齡及其地質(zhì)意義. 地質(zhì)通報, 36(8): 1343–1358.

    王金芳, 李英杰, 李紅陽, 董培培. 2018a. 內(nèi)蒙古西烏旗德勒哈達早白堊世A型花崗巖形成時代: 鋯石U-Pb定年證據(jù). 中國地質(zhì), 45(1): 197–198.

    王金芳, 李英杰, 李紅陽, 董培培. 2018b. 內(nèi)蒙古西烏旗石匠山晚侏羅世?早白堊世A型花崗巖鋯石U-Pb年齡及構(gòu)造環(huán)境. 地質(zhì)通報, 37(2–3): 382–396.

    王京彬, 王玉往, 王莉娟. 2005. 大興安嶺南段錫多金屬成礦系列. 地質(zhì)與勘探, 41(6): 15–20.

    王良玉, 廖群安, 江云川, 羅婷, 胡朝斌, 肖典, 湯帥, 劉洋. 2016. 內(nèi)蒙古錫林浩特早白堊世晚期鉀玄質(zhì)火山巖成因及構(gòu)造環(huán)境. 地質(zhì)通報, 35(6): 919–931.

    王楠, 吳才來, 秦海鵬. 2017. 川西義敦島弧中生代典型花崗巖體礦物學、地球化學特征及巖漿來源探討. 地質(zhì)論評, 63(4): 981–1000.

    王濤, 劉燊. 2013. 膠南花崗巖鋯石飽和溫度及其地質(zhì)意義. 礦物巖石地球化學通報, 32(5): 619–624.

    王喜龍, 劉家軍, 翟德高, 楊永強, 王建平, 張琪彬, 張安立. 2014. 內(nèi)蒙古邊家大院礦區(qū)石英斑巖U-Pb年代學、巖石地球化學特征及其地質(zhì)意義. 礦物巖石地球化學通報, 33(5): 654–665.

    王長明, 張壽庭, 鄧軍. 2006. 大興安嶺南段銅多金屬成礦時空結(jié)構(gòu). 成都理工大學學報(自然科學版), 33(5): 478–484.

    吳福元, 曹林. 1999. 東北亞地區(qū)的若干重要基礎地質(zhì)問題. 世界地質(zhì), 18(2): 1–13.

    吳福元, 江博明, 林強. 1997. 中國北方造山帶造山后花崗巖的同位素特點與地殼生長意義. 科學通報, 42(20): 2188–2192.

    吳福元, 李獻華, 楊進輝, 鄭永飛. 2007. 花崗巖成因研究的若干問題. 巖石學報, 23(6): 1217–1238.

    吳福元, 劉小馳, 紀偉強, 王佳敏, 楊雷. 2017. 高分異花崗巖的識別與研究. 中國科學, 47(7): 745–765.

    肖成東, 張忠良, 趙利青. 2004. 東蒙地區(qū)燕山期花崗巖Nd、Sr、Pb同位素及其巖石成因. 中國地質(zhì), 31(1): 57–63.

    徐備, 趙盼, 鮑慶中, 周永恒, 王炎陽, 羅志文. 2014. 興蒙造山帶前中生代構(gòu)造單元劃分初探. 巖石學報, 30(7): 1841–1857.

    楊奇荻, 郭磊, 王濤, 曾濤, 張磊, 童英, 史興俊. 2014. 張建軍大興安嶺中南段甘珠爾廟地區(qū)晚中生代兩期花崗巖的時代、成因、物源及其構(gòu)造背景. 巖石學報, 30(7): 1961–1981.

    楊武斌, 單強, 趙振華, 羅勇, 于學元, 李寧波, 牛賀才. 2011. 巴爾哲地區(qū)堿性花崗巖的成巖和成礦作用: 礦化和未礦化巖體的比較. 吉林大學學報(地球科學版), 41(6): 1689–1704.

    姚磊, 呂志成, 葉天竺, 龐振山, 賈宏翔, 張志輝, 吳云峰, 李睿華. 2017. 大興安嶺南段內(nèi)蒙古白音查干Sn多金屬礦床石英斑巖的鋯石U-Pb年齡、地球化學和Nd-Hf同位素特征及地質(zhì)意義. 巖石學報, 33(10): 3183–3199.

    要梅娟, 曹燁, 劉家軍, 翟德高. 2016. 內(nèi)蒙古黃崗梁鐵錫礦床輝鉬礦Re-Os年齡及其成因意義. 礦產(chǎn)勘查, 7(3): 399–403.

    要梅娟, 劉家軍, 翟德高, 王建平, 邢永亮. 2012. 大興安嶺南段多金屬成礦帶硫、鉛同位素組成及其地質(zhì)意義. 吉林大學學報(地球科學版), 42(2): 362–373.

    葉天竺. 2013. 礦床模型綜合地質(zhì)信息預測技術方法理論框架. 吉林大學學報(地球科學版), 43(7): 1053–1072.

    袁小平, 楊瑞, 邵長亮, 薛燕萍. 2016. 內(nèi)蒙古錫林浩特毛登銅礦控礦因素及礦床成因分析. 甘肅冶金, 38(1): 81–86.

    翟德高, 劉家軍, 楊永強, 王建平, 定立, 劉星旺, 張梅, 要梅娟, 蘇犁, 張紅雨. 2012. 內(nèi)蒙古黃崗梁鐵錫礦床成巖、成礦時代與構(gòu)造背景. 巖石礦物學雜志, 31(4): 513–523.

    張德全. 1993. 大興安嶺南段不同構(gòu)造環(huán)境中的兩類花崗巖. 巖石礦物學雜志, 12(1): 1–11.

    張德全, 劉勇, 李大新. 1993. 大興安嶺地區(qū)與銅多金屬有關的侵入巖//大興安嶺及鄰區(qū)銅多金屬礦床論文集. 北京: 地震出版社: 50–64.

    張東亮, 彭建堂, 胡瑞忠, 袁順達, 鄭德順. 2011. 錫石U-Pb同位素體系的封閉性及其測年的可靠性分析. 地質(zhì)論評, 57(4): 549–554.

    張連昌, 陳志廣, 周新華, 英基豐, 王非, 張玉濤. 2007. 大興安嶺根河地區(qū)早白堊世火山巖深部源區(qū)與構(gòu)造?巖漿演化: Sr-Nd-Pb-Hf同位素地球化學制約. 巖石學報, 23(11): 2823–2835.

    張巧梅, 翟東興, 李華. 2013. 內(nèi)蒙古毛登鉬錫銅礦床地質(zhì)特征及成因探討. 礦產(chǎn)勘查, 4(3): 248–256.

    張?zhí)旄? 郭碩, 辛后田, 張云, 何鵬, 劉文剛, 張闊, 王可祥, 張超. 2019. 大興安嶺南段維拉斯托高分異花崗巖體的成因與演化及其對Sn-(Li-Rb-Nb-Ta)多金屬成礦作用的制約. 地球科學, 44(1): 248-267.

    張曉玲, 張文山. 2010. 內(nèi)蒙古毛登錫銅鉬多金屬礦床花崗巖及火山角礫熔巖地質(zhì)特征及找礦意義. 礦床地質(zhì), 29(S1): 337–338.

    章榮清, 陸建軍, 王汝成, 姚遠, 丁騰, 胡加斌, 張懷峰. 2016. 湘南王仙嶺地區(qū)中生代含鎢與含錫花崗巖的巖石成因及其成礦差異機制. 地球化學, 45(2): 105– 132.

    趙國龍, 楊桂林, 傅嘉友, 楊玉琢. 1989. 大興安嶺中南部中生代火山巖. 北京: 北京科學技術出版社: 1– 260.

    趙海杰, 毛景文, 向君峰, 周振華, 魏克濤, 柯于富. 2010. 湖北銅綠山礦床石英閃長巖石的礦物學及Sr-Nd-Pb同位素特征. 巖石學報, 26(3): 768–784.

    趙輝, 李舢, 王濤, 王文志, 焦永艷. 2015. 大興安嶺南段黃崗梁地區(qū)早白堊世巖漿作用的時代、成因及其構(gòu)造意義. 地質(zhì)通報, 34(12): 2203–2218.

    趙一鳴, 王大畏, 張德全, 傅先政, 鮑修坡, 李鶴年, 艾永富. 1994. 內(nèi)蒙古東南部銅多金屬成礦地質(zhì)條件及找礦模式. 北京: 地震出版社: 1–234.

    趙振華, 熊小林, 韓小東. 1999. 花崗巖稀土元素四分組效應形成機理探討——以千里山和巴爾哲花崗巖為例. 中國科學(D輯), 29(4): 331–338.

    周振華, 呂林素, 王挨順. 2011. 內(nèi)蒙古黃崗錫鐵礦床花崗巖深部源區(qū)特征與構(gòu)造巖漿演化: Sr-Nd-Pb-Hf多元同位素制約. 地質(zhì)科技情報, 30(1): 1–14.

    周振華, 呂林素, 楊永軍, 李濤. 2010. 內(nèi)蒙古黃崗錫鐵礦區(qū)早白堊世A型花崗巖成因: 鋯石U-Pb年代學和巖石地球化學制約. 巖石學報, 26(12): 3521–3537.

    祝新友, 張志輝, 付旭, 李柏陽, 王艷麗, 焦守濤, 孫雅琳. 2016. 內(nèi)蒙古赤峰維拉斯托大型錫多金屬礦的地質(zhì)地球化學特征. 中國地質(zhì), 43(1): 188–208.

    Batchelor R A and Bowden P. 1985. Petrogeneic interpre-tation of granitoid rock series using multicationic paramters., 48: 43–55.

    Bau M. 1996. Controls on the fractionation of isovalent trace elements in magmatic and aqueous systems: Evidence from Y/Ho, Zr/Hf, and lanthanide tetrad effect., 123: 323–333.

    Bau M. 1997. The lanthanide tetrad effect in highly evolved felsic igneous rocks—A reply to the comment by Y. Pan., 128(4): 409–412.

    Chappell B W and White A J R. 1992. I- and S-type granites in the Lachlan Fold Belt., 83(1–2): 1–26.

    Chen Y J, Zhang C, Wang P, Pirajno F and Li N. 2017. The Mo deposits of Northeast China: A powerful indicator of tectonic settings and associated evolutionary trends., 81(2): 602–640.

    Clemens J D, Holloway J R and White A J R. 1986. Origin of an A-type granite: Experimental constraints., 71(3-4): 317–324.

    Collins W J, Beams S D, White A J R and Chappell B W. 1982. Nature and origin of A-type granites with particular reference to Southeastern Australia., 80(2): 189–200.

    Dickin A P. 1994. Nd isotope chemistry of Tertiary igneous rocks from Arran, Scotland: Implications for magma evolution and crustal structure., 131(3): 329–333.

    Dostal J and Chatterjee A K. 2000. Contrasting behaviour of Nb/Ta and Zr/Hf ratios in a peraluminous granitic pluton (Nova Scotia, Canada)., 163(1): 207–218.

    Eby G N. 1992. Chemical subdivision of the A-type granitoids: Petrognetic and tectonic implications., 20(7): 641–644.

    Engebretson D C, Cox A and Gordon R G. 1985. Relative motions between oceanic and continental plates in the Pacific basin., 206: 1–59.

    Fan W M, Guo F, Wang Y J and Lin G. 2003. Late Mesozoic calcalkaline volcanism of post orogenic extension in the northern Da Hinggan Mountains, northeastern China., 121(1?2): 115–135.

    Guo C L, Chen Y C, Zeng Z L and Lou F S. 2012. Petrogenesis of the Xihuashan granites in southeastern China: Constraints from geochemistry and in-situ analyses of zircon U-Pb-Hf-O isotopes., 148: 209–227.

    Han B F, Wang S G, Jan B M, Hong D W, Kagami H and Sun Y L. 1997. Depleted-mantle source for the Ulungur River A-type granites from North Xinjiang, China: Geochemistry and Nd-Sr isotopic evidence, and impli-cations for Phanerozoic crustal growth., 138(3–4): 135–159.

    Harris C, Marsh J S and Milner S C. 1999. Petrology of the alkaline core of the Messum igneous complex, Nambibia: Evidence for the progressively decreasing effect of crustal contamination., 40(9): 1377–1397.

    Irber W. 1999. The lanthanide tetrad effect and its correlation with K/Rb, Eu/Eu*, Sr/Eu, Y/Ho, and Zu/Hf of evolving peraluminous granite suites., 63(3–4): 489–508.

    Irvine A J and Baragar W R A. 1971. A guide to the chemical classification of the common volcanic rocks., 8(5): 523–548.

    Jahn B M, Wu F Y and Chen B. 2000. Massive granitoid generation in Central Asia: Nd isotope evidence and implications for continental growth in the Phanerozoic., 23(2): 82–92.

    Jahn B M, Wu F Y, Capdevila R, Martineau F, Zhao Z H and Wang Y X. 2001. Highly evolved juvenile granites with tetrad REE patterns: the Woduhe and Baerzhe granites from the Great Xingan Mountains in NE China., 59(4): 171–198.

    Jahn B M, Wu F Y, Lo C H and Tsai C H. 1999. Crust-mantle interaction induced by deep subduction of the continental crust: Geochemical and Sr-Nd isotropic evidence from post-collosional maifc-ultramafic intrusion of the northern Dabie complex, central China., 157(1): 119–146.

    Jiang S H, Bagas L, Hu P, Han N, Chen C L, Liu Y and Kang H. 2016. Zircon U-Pb ages and Sr-Nd-Hf isotopes of the highly fractionated granite with tetrad REE patterns in the Shamai tungsten deposit in eastern Inner Mongolia, China: Implications for the timing of mineralization and ore genesis., 261: 322–339.

    King P L, Chappell B W, Allen C M and White A J R. 2001. Are A-type granites the high-temperature felsic granites? Evidence from fractionated granites of the Wangrah Suite., 48(4): 501–514.

    King P L, White A J R, Chappell B W and Allen C M. 1997. Characterization and origin of aluminous A-type granites from the Lachlan fold belt, Southeastern Australia., 38(3): 371–391.

    Le Maitre R W. 1976. The chemical variability of some common igneous rocks., 17(4): 589–598.

    Liu Y F, Jiang S H and Bagas L. 2016. The genesis of metal zonation in the Weilasituo and Bairendaba Ag-Zn-Pb- Cu-(Sn-W) deposits in the shallow part of a porphyry Sn-W-Rb system, Inner Mongolia, China., 75: 150–173.

    Liu Y S, Gao S, Hu Z C, Gao S, Gunther D, Xu J, Gao C G and Chen H H. 2008.analysis of major and trace elements of anhydrous minerals by LA-ICP-MS without applying an internal standard., 257(1?2): 34–43.

    Ludwig K R. 2003. Isoplot 3.00: A Geochronological Toolkit for Microsoft Excel. Berkeley Geochronology Center Special Publication: 1–70.

    Maniar P D and Piccoli P M. 1989. Tectonic discrimination of granitoids., 101(5): 635–643.

    Middlemost E A K. 1994. Naming materials in the magma/ igneous rock system.s, 37(3–4): 215–224.

    Mingram B, Trumbull R B, Littman S and Gertenberger H. 2000. A Petrogenetic study of anorogenic felsic magmatism in the Cretaceous Paresis ring complex, Namibia: Eviendence for mixing of crust and mantle-derived components., 54(1): 1–22.

    Mushkin A, Navon O, Halicz L, Hartmann G and Stein M. 2003. The petrogenesis of A-type magmas from the Amram Massif, Southern Israel., 44(5): 815–832.

    Nasdala L, Hofmeister W G, Norberg N, Mattinson J M, Corfu F, Dorr W, Kamo SL, Kennedy A K, Kronz A and Reiners P W. 2008. Zircon M257: A homogeneous natural reference material for the ion microprobe U-Pb analysis of zircon., 32(3): 247–265.

    Ouyang H G, Mao J W, Zhou Z H and Su H M. 2015. Late Mesozoic metallogeny and intracontinental magmatism, southern Great Xing’an Range, northeastern China., 27(3): 1153–1172.

    Pearce J A, Harris N B W and Tindle A G. 1984. Trace element discrimination diagrams for the tectonic interpretation of granitic rocks., 25(4): 956–983.

    Peccerillo A and Taylor S R. 1976. Geochemistry of Eocene calc-alkaline volcanic rocks from the Kastamonu area, Northern Turkey., 58(1): 63–81.

    Rapp R P and Watson E B. 1995. Dehydration melting of metabasalt at 8–32 kbar: Implications for continental growth and crust-mantle recycling., 36(4): 891–931.

    Streckeisen A L. 1973. Plutonic rocks: Classification and nomenclature recommended by the IUGS subcom-mission on the systematics of igneous rocks., 18(7): 26–30.

    Sun S S and McDonough W F. 1989. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes // Sauders A D and Norry M J. Magmatism in the Ocean Basins.,,, 42(1): 313–345.

    Wang F X, Bagas L, Jiang S H and Liu Y F. 2017. Geological, geochemical, and geochronological charac-teristics of Weilasituo Sn-polymetal deposit, Inner Mongolia, China., 80: 1206–1229.

    Watson E B and Harrison T M. 1983. Zircon saturation revisited: Temperature and composition effect in a variety of crustal magmas types., 64(2): 295–304.

    Whalen J B, Currie K L and Chappell B W. 1987. A-type granites: Geochemical characteristics, discrimination and petrogenesis., 95(4): 407–419.

    Wu F Y, Jahn B M, Wilde S A and Sun D Y. 2000. Phanerozoic continental crustal growth: U-Pb and Sr-Nd isotopic evidence from the granites in northeastern China., 328(1?2): 89–113.

    Wu F Y, Jahn B M, Wilde S A, Lo C H, Yui T F, Lin Q, Ge W C, Sun D Y. 2003. Highly fractionated I-type granites in NE China (II): Isotopic geochemistry and implications for crustal growth in the Phanerozoic., 67(3?4): 191–204.

    Wu F Y, Sun D Y, Li H M, Jahn B M and Wilde S A. 2002. A-type granites in northeastern China: Age and geoche-mical constraints on their petrogenesis., 187(1?2): 143–173.

    Wu F Y, Sun D Y, Jahn B M and Wilde S A. 2004. A Jurassic garnet-bearing granitic pluton from NE China showing tetrad REE patterns., 23(5): 731–744.

    Wyllie P J. 1984. Constraints imposed by experimental petrology on possible and impossible magma sources and products., 310(1514): 439–456.

    Yuan S D, Peng J T, Hao S, Li H M, Geng J Z and Zhang D L. 2011. In situ LA-MC-ICP-MS and ID-TIMS U-Pb geochronology of cassiterite in the giant Furong tin deposit, Hunan Province, South China: New constraints on the timing of tin-polymetallic mineralization., 43(1): 235–242.

    Zartman R E and Doe B R. 1981. Lead-isotope evolution. U. S. Geological Survey Professional Paper: 169–170.

    Zhai D G, Liu J J, Zhang A L and Sun Y Q. 2017. U-Pb, Re-Os, and40Ar/39Ar geochronology of porphyry Sn±Cu±Mo and polymetallic Ag-Pb-Zn-Cu vein mineralization at Bianjiadayuan, InnerMongolia, Northeast China: Implications for discrete mineralization events., 112(8): 2041–2059.

    Zhang L L, Jiang S H, Bagas L, Han N, Liu Y and Liu Y F. 2019. Element behaviour during interaction of magma and fluid: A case study of Chamuhan Granite, and implications on the genesis of W-Mo mineralization., (342?343): 31–34.

    Zhao Z H, Xiong X L, Han Z D, Wang Y X, Wang Q, Bao Z W and Jahn B M. 2002. Controls on the REE tetrad effect in granites: Evidence from the Qianlishan and Baerzhe granites, China., 36(6): 527–543.

    Zindler A and Hart S. 1986. Chemical geodynamics., 14(1): 493–571.

    Metallogenetic Control of Magmatism on the Maodeng Sn-Cu Deposit in the Southern Great Xing’an Range: Evidence from Geochronology, Geochemistry, and Sr-Nd-Pb Isotopes

    JI Genyuan1, 2, JIANG Sihong1*, LI Gaofeng1, 3, YI Jinjun2, ZHANG Lili1and LIU Yifei1

    (1.MNR Key Laboratory of Metallogeny and Mineral Assessment, Institute of Mineral Resources, Chinese Academy of Geological Sciences, Beijing 100037,China; 2.Cores and Samples Center of Natural Resources, Langfang 065201, Hebei, China; 3. Key Laboratory of Orogen and Crust Evolution, Peking University, Beijing 100871,China)

    The Maodeng Sn-Cu deposit is a medium-sized deposit in the southern segment of the Great Xing’an Range located in Xinlinhot City of Inner Mongolia. Whole rock geochemistry, LA-ICP-MS zircon and cassiterite U-Pb ages, and Sr-Nd-Pb isotopic compositions of the granite porphyry and Sn-Cu ores in the Maodeng deposit have been studied in detail and are reported in this contribution. The zircon from the ore-related granite porphyry has a weighted mean U-Pb age of 138±0.6 Ma, while the cassiterite from cassiterite-sulphide-quartz vein-type ore has a Tera-Wasserburg concordant U-Pb age of 139±3.2 Ma. These two ages are consistent within error, which indicates that the emplacement of granite porphyry and Sn-polymetallic mineralization in Maodeng took place during the Early Cretaceous. Geoche-mically, the granite porphyry is characterized by high SiO2(73.18% to 77.3%), (Na2O+K2O) (8.00% to 8.69%), and K2O (4.80% to 5.22%), with low CaO (0.49% to 0.95%), Al2O3(11.67% to 12.83%), MgO (0.12% to 0.46%), and relatively high FeOT/(MgO+FeOT) (0.74 to 0.91), and 10000×Ga/Al (3.64 to 4.28) ratios. The rock is enriched in Rb, Th, K, Hf, and depleted in Ba, Sr, P and Ti. The REE distribution patterns of the granite porphyry are of typical gull-wing shape with obviously negative Eu anomalies (δEu=0.05 to 0.2). Zircon saturation temperatures of the granite porphyry vary from 817 to 861 ℃. All these show that the granite porphyry exhibits geochemical characteristics typical for A-type granite. The granite porphyry has relatively low (87Sr/86Sr)ivalues (0.702907 to 0.704506), high εNd() values (+0.8 to +4.4) and young two-stage Nd model ages (DM2=572 to 863 Ma), suggesting that the rock was dominantly derived from partial melting of a juvenile lower crust, followed by strong fractional crystallization during magma ascent. The Pb isotope compositions of the granite porphyry and ores are similar, indicating that the Pb of the ores could be sourced from the granite porphyry. According to the trace element characteristics, chemical subdivision diagrams and tectonic environmentdiscrimination diagrams, the granite porphyry in the Maodeng belongs to A2-type granite, which formed in an extensional setting. Based on previous studies and combining geological characteristics, ages of the ores and granite porphyry, and Pb istope characteristics, it is believed that the intensely fractionated granite porphyry at the Maodeng deposit is spatiallyand temporally, and possibly genetically related to the Sn-Cu mineralization. It is proposed that the size of the Sn deposit is constrained by the degree of fractional crystallization of the magma.

    southern segment of the Great Xing’an Range; Maodeng Sn-Cu deposit; cassiterite U-Pb dating; zircon U-Pb dating; Sr-Nd-Pb isotopes

    2020-02-04;

    2020-06-19

    科技部深地項目(2017YFC0601303)、國家自然科學基金項目(41873051)、地科院基本業(yè)務費項目(YYWF201715)和地質(zhì)調(diào)查項目(DD20190437、DD20190411)聯(lián)合資助。

    季根源(1987–), 男, 工程師, 博士研究生, 礦物學、巖石學、礦床學專業(yè)。Email: csujigenyuan@163.com

    江思宏(1968–), 男, 研究員, 博士生導師, 主要從事金屬礦床成礦規(guī)律研究。Email: jiangsihong1@163.com

    P597; P612

    A

    1001-1552(2021)04-0681-024

    10.16539/j.ddgzyckx.2021.04.004

    猜你喜歡
    成礦
    歡迎訂購2022年《大地構(gòu)造與成礦學》
    歡迎訂購2022年《大地構(gòu)造與成礦學》
    淺析淶源縣大黃峪鉬礦成礦遠景
    桂西沉積型鋁土礦床成礦規(guī)律及成礦模式
    歡迎訂購2021年《大地構(gòu)造與成礦學》
    歡迎訂購2021年《大地構(gòu)造與成礦學》
    新疆寨北山銅礦成礦新認識及找礦預測
    氧化還原作用對鈾成礦的意義——以相山西部流紋英安巖鈾礦為例
    中非銅鈷成礦帶礦床遙感特征研究
    柴達木盆地北緣鋰多金屬礦成礦條件及找礦潛力
    99riav亚洲国产免费| 欧美激情久久久久久爽电影| 老熟妇乱子伦视频在线观看| av天堂在线播放| 男女床上黄色一级片免费看| 久久久久九九精品影院| 最近最新中文字幕大全电影3| 搞女人的毛片| 一级a爱片免费观看的视频| 一个人免费在线观看的高清视频| 后天国语完整版免费观看| 久久久久精品国产欧美久久久| 亚洲欧美日韩东京热| 国产男靠女视频免费网站| 国产一区二区在线av高清观看| 免费观看的影片在线观看| 国产精品永久免费网站| 小说图片视频综合网站| 亚洲,欧美精品.| 亚洲五月天丁香| 欧美3d第一页| 国产欧美日韩精品一区二区| 国产极品精品免费视频能看的| 国产精品久久久久久精品电影| 久久久久久久精品吃奶| 日本在线视频免费播放| 亚洲熟女毛片儿| 九九久久精品国产亚洲av麻豆 | 国产成人一区二区三区免费视频网站| 小蜜桃在线观看免费完整版高清| 国产av在哪里看| 噜噜噜噜噜久久久久久91| 成人av在线播放网站| 高潮久久久久久久久久久不卡| 一二三四在线观看免费中文在| 99re在线观看精品视频| 亚洲五月天丁香| 亚洲avbb在线观看| 亚洲精品一卡2卡三卡4卡5卡| 国产一区在线观看成人免费| a级毛片在线看网站| 午夜福利成人在线免费观看| 国产精品一区二区三区四区免费观看 | ponron亚洲| 五月玫瑰六月丁香| 日本一二三区视频观看| 成人一区二区视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 淫秽高清视频在线观看| 三级毛片av免费| 免费一级毛片在线播放高清视频| 99国产极品粉嫩在线观看| 日韩欧美一区二区三区在线观看| 久久中文看片网| 男女床上黄色一级片免费看| 国产主播在线观看一区二区| 欧美午夜高清在线| 久久久成人免费电影| 久久天堂一区二区三区四区| 黄片大片在线免费观看| 中文在线观看免费www的网站| 日本成人三级电影网站| 一二三四社区在线视频社区8| 午夜精品在线福利| 看片在线看免费视频| 欧美一区二区精品小视频在线| 五月玫瑰六月丁香| 久9热在线精品视频| 一级毛片女人18水好多| 亚洲精品在线观看二区| 国产精品女同一区二区软件 | 国产一区二区在线观看日韩 | 久久久色成人| 成人欧美大片| 神马国产精品三级电影在线观看| 日本a在线网址| 精品久久蜜臀av无| av福利片在线观看| 国产精品一区二区三区四区久久| 久9热在线精品视频| 日日摸夜夜添夜夜添小说| 大型黄色视频在线免费观看| 视频区欧美日本亚洲| 三级毛片av免费| 欧美日本视频| 午夜福利欧美成人| 高清毛片免费观看视频网站| 美女被艹到高潮喷水动态| 国产爱豆传媒在线观看| 亚洲无线在线观看| 99久久综合精品五月天人人| 一级黄色大片毛片| 国产一区二区在线av高清观看| 三级毛片av免费| 久久久久免费精品人妻一区二区| 日韩 欧美 亚洲 中文字幕| 亚洲国产精品成人综合色| 国产久久久一区二区三区| 免费av不卡在线播放| 国产蜜桃级精品一区二区三区| 中国美女看黄片| 制服人妻中文乱码| 男人舔奶头视频| 久久中文字幕一级| www.熟女人妻精品国产| 悠悠久久av| 国产精品日韩av在线免费观看| 国产成人av教育| 欧美乱码精品一区二区三区| 99久久精品热视频| 免费观看的影片在线观看| 看免费av毛片| ponron亚洲| 亚洲狠狠婷婷综合久久图片| 日韩欧美在线二视频| 小说图片视频综合网站| 长腿黑丝高跟| 欧美中文日本在线观看视频| 在线观看免费视频日本深夜| 视频区欧美日本亚洲| 欧美丝袜亚洲另类 | 性色avwww在线观看| 黄色视频,在线免费观看| 亚洲第一欧美日韩一区二区三区| 国产精品99久久99久久久不卡| av欧美777| 丝袜人妻中文字幕| 麻豆一二三区av精品| 日本五十路高清| av视频在线观看入口| 色尼玛亚洲综合影院| 亚洲精品在线美女| 久久中文字幕一级| 久久国产精品人妻蜜桃| 人人妻,人人澡人人爽秒播| 搡老妇女老女人老熟妇| www.熟女人妻精品国产| 波多野结衣高清无吗| 男女之事视频高清在线观看| 亚洲在线观看片| 又粗又爽又猛毛片免费看| 免费观看人在逋| 亚洲国产精品久久男人天堂| 国产爱豆传媒在线观看| 国产精品九九99| 日韩有码中文字幕| av黄色大香蕉| 国产单亲对白刺激| 真实男女啪啪啪动态图| 亚洲一区二区三区不卡视频| 老司机在亚洲福利影院| 国产乱人视频| 一区福利在线观看| 久久中文看片网| 亚洲熟妇中文字幕五十中出| 琪琪午夜伦伦电影理论片6080| 一个人免费在线观看电影 | 国产亚洲精品综合一区在线观看| 757午夜福利合集在线观看| 久久香蕉精品热| 人人妻人人澡欧美一区二区| 久久婷婷人人爽人人干人人爱| av在线蜜桃| 亚洲国产欧美人成| 给我免费播放毛片高清在线观看| 99riav亚洲国产免费| 99热精品在线国产| 国产精品亚洲一级av第二区| 亚洲 国产 在线| 国内精品美女久久久久久| 香蕉av资源在线| 日韩av在线大香蕉| 一本精品99久久精品77| 国产午夜福利久久久久久| 九色成人免费人妻av| 午夜激情欧美在线| 老熟妇乱子伦视频在线观看| 国产蜜桃级精品一区二区三区| 久久久久国内视频| 成人精品一区二区免费| 久久婷婷人人爽人人干人人爱| 国产精品 国内视频| 亚洲国产精品久久男人天堂| 精品一区二区三区视频在线观看免费| 久久中文字幕人妻熟女| 国产又黄又爽又无遮挡在线| av中文乱码字幕在线| 欧美激情在线99| 国产久久久一区二区三区| 亚洲七黄色美女视频| 国产精品亚洲av一区麻豆| 99久久99久久久精品蜜桃| 成人高潮视频无遮挡免费网站| 欧美日本视频| 欧美日韩中文字幕国产精品一区二区三区| 日韩欧美三级三区| 亚洲成人中文字幕在线播放| 色吧在线观看| 国产精品女同一区二区软件 | 亚洲午夜精品一区,二区,三区| 国产成+人综合+亚洲专区| 欧美高清成人免费视频www| 成人欧美大片| 日韩成人在线观看一区二区三区| 后天国语完整版免费观看| 亚洲av片天天在线观看| 久久精品国产清高在天天线| 99久久精品国产亚洲精品| 亚洲av五月六月丁香网| 琪琪午夜伦伦电影理论片6080| 午夜福利成人在线免费观看| 一进一出抽搐gif免费好疼| 在线观看美女被高潮喷水网站 | 欧美xxxx黑人xx丫x性爽| 夜夜躁狠狠躁天天躁| 女人高潮潮喷娇喘18禁视频| 国产亚洲av嫩草精品影院| 亚洲av片天天在线观看| 五月玫瑰六月丁香| 99久国产av精品| 成人特级黄色片久久久久久久| a级毛片a级免费在线| 精品99又大又爽又粗少妇毛片 | 国产成人欧美在线观看| 12—13女人毛片做爰片一| 日本免费a在线| 午夜精品一区二区三区免费看| 国产成人系列免费观看| 操出白浆在线播放| 在线观看免费视频日本深夜| 亚洲欧美一区二区三区黑人| 看片在线看免费视频| 在线观看一区二区三区| 很黄的视频免费| 亚洲专区国产一区二区| 国产精品女同一区二区软件 | 免费观看人在逋| 天天躁狠狠躁夜夜躁狠狠躁| 一a级毛片在线观看| 一本久久中文字幕| 亚洲av五月六月丁香网| 偷拍熟女少妇极品色| 精品久久久久久久久久久久久| 亚洲精品在线观看二区| 成年人黄色毛片网站| 亚洲无线观看免费| 国产欧美日韩精品亚洲av| 精品午夜福利视频在线观看一区| 嫩草影视91久久| 亚洲精华国产精华精| 一区二区三区国产精品乱码| 岛国视频午夜一区免费看| 美女高潮的动态| 老汉色∧v一级毛片| 亚洲国产精品成人综合色| 丁香六月欧美| 在线免费观看的www视频| 久久草成人影院| 亚洲一区二区三区色噜噜| 精品免费久久久久久久清纯| 亚洲av中文字字幕乱码综合| 日韩欧美三级三区| 国产精品 欧美亚洲| 国产三级中文精品| www.www免费av| 欧美日韩国产亚洲二区| 亚洲精品国产精品久久久不卡| 校园春色视频在线观看| 2021天堂中文幕一二区在线观| 高清毛片免费观看视频网站| 色老头精品视频在线观看| 九九热线精品视视频播放| 黄色丝袜av网址大全| 99国产极品粉嫩在线观看| 十八禁人妻一区二区| 日本熟妇午夜| 久久久精品欧美日韩精品| 最近最新中文字幕大全电影3| 少妇人妻一区二区三区视频| 美女 人体艺术 gogo| 亚洲国产高清在线一区二区三| 亚洲av第一区精品v没综合| 成人欧美大片| АⅤ资源中文在线天堂| 亚洲av成人精品一区久久| 国产不卡一卡二| 亚洲成a人片在线一区二区| 12—13女人毛片做爰片一| 一边摸一边抽搐一进一小说| 免费无遮挡裸体视频| 国产精品自产拍在线观看55亚洲| 国产精品野战在线观看| 亚洲国产欧洲综合997久久,| 最近最新中文字幕大全电影3| 中文字幕人成人乱码亚洲影| 国产单亲对白刺激| 亚洲成人精品中文字幕电影| xxxwww97欧美| 久久欧美精品欧美久久欧美| 久9热在线精品视频| 亚洲人成网站在线播放欧美日韩| 日韩av在线大香蕉| 成人国产综合亚洲| 51午夜福利影视在线观看| 真人做人爱边吃奶动态| АⅤ资源中文在线天堂| 中文字幕久久专区| 日韩欧美一区二区三区在线观看| 国产精品久久电影中文字幕| www日本黄色视频网| 日本 欧美在线| 国内精品久久久久久久电影| 欧美性猛交╳xxx乱大交人| 母亲3免费完整高清在线观看| 国产精品 欧美亚洲| 久久国产精品影院| 色综合亚洲欧美另类图片| 久久精品91无色码中文字幕| www.熟女人妻精品国产| 国产精品亚洲一级av第二区| 美女高潮喷水抽搐中文字幕| 禁无遮挡网站| 日韩欧美国产在线观看| 不卡一级毛片| 丰满人妻一区二区三区视频av | 日韩有码中文字幕| 久久久久久久久中文| avwww免费| 亚洲欧美日韩高清专用| 亚洲成av人片在线播放无| 少妇丰满av| 嫩草影院精品99| 母亲3免费完整高清在线观看| 国产综合懂色| 不卡av一区二区三区| 亚洲熟女毛片儿| 日本 av在线| 叶爱在线成人免费视频播放| 在线免费观看不下载黄p国产 | 国产一级毛片七仙女欲春2| 亚洲专区中文字幕在线| 91av网一区二区| 一进一出好大好爽视频| 欧美日韩一级在线毛片| 91麻豆av在线| 国产精品一区二区精品视频观看| 免费大片18禁| 国产精品亚洲美女久久久| 亚洲无线在线观看| 中文字幕av在线有码专区| 亚洲精品美女久久av网站| 人妻丰满熟妇av一区二区三区| 免费av毛片视频| 熟女少妇亚洲综合色aaa.| 日韩欧美精品v在线| 老司机午夜十八禁免费视频| 无遮挡黄片免费观看| 桃色一区二区三区在线观看| 亚洲国产精品999在线| 在线看三级毛片| www.自偷自拍.com| 亚洲人成伊人成综合网2020| 国产精品久久久久久亚洲av鲁大| 成年免费大片在线观看| 丁香欧美五月| 亚洲av电影不卡..在线观看| 久久亚洲精品不卡| 亚洲第一电影网av| 欧美黄色片欧美黄色片| 免费观看精品视频网站| 床上黄色一级片| 一级毛片女人18水好多| 亚洲精品粉嫩美女一区| 亚洲av中文字字幕乱码综合| 国产精品精品国产色婷婷| 看黄色毛片网站| 女生性感内裤真人,穿戴方法视频| 可以在线观看毛片的网站| 一二三四社区在线视频社区8| 99热这里只有是精品50| 老汉色av国产亚洲站长工具| 国产av麻豆久久久久久久| 琪琪午夜伦伦电影理论片6080| 麻豆一二三区av精品| 国产一区二区激情短视频| 99热精品在线国产| 亚洲国产欧美人成| 好看av亚洲va欧美ⅴa在| 免费观看人在逋| 免费看光身美女| 99国产精品一区二区蜜桃av| 精品不卡国产一区二区三区| 国产亚洲av嫩草精品影院| 999精品在线视频| 久久天躁狠狠躁夜夜2o2o| 国产精品野战在线观看| 久久久国产欧美日韩av| av福利片在线观看| 美女高潮喷水抽搐中文字幕| 欧美日韩福利视频一区二区| 国产伦精品一区二区三区四那| avwww免费| 非洲黑人性xxxx精品又粗又长| 国产精品一区二区三区四区免费观看 | 听说在线观看完整版免费高清| 日韩欧美国产在线观看| 亚洲自拍偷在线| 亚洲色图 男人天堂 中文字幕| 亚洲最大成人中文| 日韩三级视频一区二区三区| 日本三级黄在线观看| 人人妻人人澡欧美一区二区| 免费看十八禁软件| 特级一级黄色大片| 亚洲aⅴ乱码一区二区在线播放| 特级一级黄色大片| 一个人观看的视频www高清免费观看 | 日韩欧美精品v在线| 欧美成人免费av一区二区三区| 日韩 欧美 亚洲 中文字幕| 久久久久国产一级毛片高清牌| 一区福利在线观看| 最近最新中文字幕大全电影3| 亚洲精品美女久久久久99蜜臀| 脱女人内裤的视频| 亚洲,欧美精品.| 国产久久久一区二区三区| 国产精品久久电影中文字幕| 美女 人体艺术 gogo| 麻豆一二三区av精品| 久久久久久九九精品二区国产| 一个人免费在线观看电影 | www日本黄色视频网| av片东京热男人的天堂| 每晚都被弄得嗷嗷叫到高潮| 偷拍熟女少妇极品色| 亚洲男人的天堂狠狠| 在线视频色国产色| 69av精品久久久久久| 国产成+人综合+亚洲专区| 欧美激情在线99| 日本成人三级电影网站| 国产亚洲精品久久久com| 亚洲成人久久性| 热99re8久久精品国产| 最好的美女福利视频网| 国产激情久久老熟女| 日韩欧美一区二区三区在线观看| 欧美一级a爱片免费观看看| 19禁男女啪啪无遮挡网站| 亚洲无线观看免费| 免费人成视频x8x8入口观看| 日本成人三级电影网站| av女优亚洲男人天堂 | 校园春色视频在线观看| 国产精品,欧美在线| 美女被艹到高潮喷水动态| 亚洲狠狠婷婷综合久久图片| 不卡av一区二区三区| 国模一区二区三区四区视频 | 久久精品国产99精品国产亚洲性色| 欧美在线黄色| 真实男女啪啪啪动态图| 女生性感内裤真人,穿戴方法视频| 欧美乱码精品一区二区三区| 最近最新中文字幕大全免费视频| 一级黄色大片毛片| 色精品久久人妻99蜜桃| 757午夜福利合集在线观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲国产精品久久男人天堂| 男人舔奶头视频| 18禁裸乳无遮挡免费网站照片| 一级毛片高清免费大全| 男女做爰动态图高潮gif福利片| 熟女少妇亚洲综合色aaa.| 亚洲,欧美精品.| 中文字幕av在线有码专区| av视频在线观看入口| 女同久久另类99精品国产91| 女警被强在线播放| 少妇人妻一区二区三区视频| 成人18禁在线播放| 久久久久久久久久黄片| 少妇的逼水好多| 久久久久久久精品吃奶| 国内精品美女久久久久久| e午夜精品久久久久久久| 国产又色又爽无遮挡免费看| 真实男女啪啪啪动态图| or卡值多少钱| 曰老女人黄片| 麻豆成人午夜福利视频| 97人妻精品一区二区三区麻豆| 丰满人妻熟妇乱又伦精品不卡| 亚洲熟妇中文字幕五十中出| 日韩欧美在线乱码| 欧美成人性av电影在线观看| 在线观看免费视频日本深夜| 在线免费观看不下载黄p国产 | 中文在线观看免费www的网站| 国产成人福利小说| 无限看片的www在线观看| 18禁裸乳无遮挡免费网站照片| www.自偷自拍.com| 日韩欧美在线二视频| 日韩三级视频一区二区三区| 少妇熟女aⅴ在线视频| a在线观看视频网站| 亚洲 国产 在线| 精品99又大又爽又粗少妇毛片 | 天天一区二区日本电影三级| 毛片女人毛片| 免费看a级黄色片| 欧洲精品卡2卡3卡4卡5卡区| 国产一区二区在线av高清观看| 日韩国内少妇激情av| 俺也久久电影网| 嫩草影视91久久| 日本三级黄在线观看| 久久久久久国产a免费观看| 国产精品国产高清国产av| 香蕉国产在线看| 在线视频色国产色| 天堂网av新在线| 久久中文字幕一级| 亚洲 欧美 日韩 在线 免费| 99久久精品国产亚洲精品| 国产欧美日韩一区二区三| 看片在线看免费视频| 夜夜爽天天搞| 又粗又爽又猛毛片免费看| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品粉嫩美女一区| 国产乱人伦免费视频| 国产麻豆成人av免费视频| 一a级毛片在线观看| 日日摸夜夜添夜夜添小说| 欧美成人性av电影在线观看| 国产高清视频在线播放一区| 亚洲男人的天堂狠狠| 久久精品综合一区二区三区| 亚洲av五月六月丁香网| 国产成人av教育| 黄色女人牲交| 久久精品国产综合久久久| 十八禁人妻一区二区| 午夜影院日韩av| 午夜久久久久精精品| 叶爱在线成人免费视频播放| 成年女人看的毛片在线观看| 亚洲第一欧美日韩一区二区三区| 中文字幕人妻丝袜一区二区| 国产伦在线观看视频一区| 高清毛片免费观看视频网站| 久久国产乱子伦精品免费另类| 国产又黄又爽又无遮挡在线| 一个人免费在线观看电影 | 又爽又黄无遮挡网站| 五月玫瑰六月丁香| 国产免费男女视频| 麻豆av在线久日| 欧美日韩黄片免| 操出白浆在线播放| 国产伦在线观看视频一区| 免费在线观看影片大全网站| 精品欧美国产一区二区三| 三级国产精品欧美在线观看 | 国产精品99久久久久久久久| 精品国产乱码久久久久久男人| 无限看片的www在线观看| 久久久久久久精品吃奶| 九九在线视频观看精品| 精品一区二区三区视频在线 | 日本五十路高清| 人妻久久中文字幕网| 国产蜜桃级精品一区二区三区| 精品一区二区三区av网在线观看| 哪里可以看免费的av片| 真人做人爱边吃奶动态| 午夜亚洲福利在线播放| 黄色日韩在线| 国产欧美日韩一区二区精品| 午夜福利在线在线| 亚洲欧美精品综合久久99| 婷婷六月久久综合丁香| 一级毛片女人18水好多| 听说在线观看完整版免费高清| 又黄又爽又免费观看的视频| 久久天躁狠狠躁夜夜2o2o| 悠悠久久av| 日本黄色片子视频| 久久亚洲精品不卡| 久久久久久久久免费视频了| 亚洲aⅴ乱码一区二区在线播放| 国产精品久久久久久亚洲av鲁大| 黑人欧美特级aaaaaa片| 亚洲七黄色美女视频| 又粗又爽又猛毛片免费看| 真实男女啪啪啪动态图| 欧美zozozo另类| 国产av在哪里看| 中文字幕久久专区| 久久精品91无色码中文字幕| 少妇人妻一区二区三区视频| 真人做人爱边吃奶动态| 免费观看的影片在线观看| 巨乳人妻的诱惑在线观看| 国产伦精品一区二区三区四那| 桃红色精品国产亚洲av| 美女午夜性视频免费| 一级毛片高清免费大全| 国产亚洲av嫩草精品影院| 国产高清激情床上av| 特级一级黄色大片|