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

    巖漿鎳銅鉑族礦床成礦過程中流體的作用
    ——對(duì)小巖體超大型礦床的啟示

    2015-01-04 02:25:21張銘杰湯慶艷李文淵胡沛青李建平
    中國工程科學(xué) 2015年2期
    關(guān)鍵詞:硫化物巖漿礦床

    張銘杰,湯慶艷,李文淵,2,余 明,胡沛青,李建平

    (1.蘭州大學(xué)地質(zhì)科學(xué)與礦產(chǎn)資源學(xué)院,蘭州730000;2.中國地質(zhì)調(diào)查局西安地質(zhì)調(diào)查中心,西安710054)

    巖漿鎳銅鉑族礦床成礦過程中流體的作用
    ——對(duì)小巖體超大型礦床的啟示

    張銘杰1,湯慶艷1,李文淵1,2,余 明1,胡沛青1,李建平1

    (1.蘭州大學(xué)地質(zhì)科學(xué)與礦產(chǎn)資源學(xué)院,蘭州730000;2.中國地質(zhì)調(diào)查局西安地質(zhì)調(diào)查中心,西安710054)

    超大型巖漿鎳、鉑族元素(PGE)硫化物礦床控制了Ni-PGE的資源量,本文根據(jù)地幔巖漿事件中Ni、Cu和PGE等在部分熔融、巖漿結(jié)晶和硫化物熔離過程的行為,討論了巖漿Ni-Cu-PGE硫化物礦床的成礦巖漿類型與規(guī)模、成礦金屬元素聚集的方式、巖漿演化過程中硫飽和度與硫化物液相不混溶的控制因素,探討了流體對(duì)成礦金屬的運(yùn)移聚集作用等。認(rèn)為地幔大規(guī)模的高度部分熔融(即要求高M(jìn)g巖漿)促使地幔硫化物及橄欖石中巨量的Ni和PGE進(jìn)入巖漿,硫化物熔離富集成礦金屬形成大型層狀(或通道)巖體賦存的超大型Ni-Cu-PGE硫化物礦床;地幔超臨界流體可能大范圍內(nèi)搬運(yùn)聚集大量Ni-PGE金屬元素,對(duì)巖漿硫化物礦床成礦可能有重要貢獻(xiàn),是小規(guī)模巖漿(小巖體)成礦的可能機(jī)制。

    流體;小巖體;幔源巖漿;超大型巖漿礦床;成礦作用

    1 前言

    鎳、銅和鉑族元素(PGE)由于特殊的用途是十分重要的戰(zhàn)略資源,Ni和PGE資源主要來源于與鎂鐵質(zhì)-超鎂鐵質(zhì)巖相關(guān)的礦床,如巖漿硫化物型鎳礦、大型層狀雜巖體鉑族礦床及紅土型鎳礦床。巖漿硫化物型鎳礦床占全球鎳儲(chǔ)量(約1.16×1010t)的28%,紅土型鎳礦床占72%[1]。全球具有經(jīng)濟(jì)價(jià)值的約99%的PGE資源和40%的Ni資源都賦存于與幔源巖漿有關(guān)的巖漿硫化物礦床中[2,3]。紅土型鎳礦床中約85%與增生地體內(nèi)的純橄巖、方輝橄欖巖和橄欖巖有關(guān),其余15%與克拉通內(nèi)科馬提巖和層狀鎂鐵質(zhì)-超鎂鐵質(zhì)巖漿巖有關(guān)。因此,巖漿型Ni-Cu-PGE硫化物礦床是Ni與PGE資源最終的重要成礦類型。

    巖漿Ni-Cu-PGE硫化物礦床品位較高、雜質(zhì)較少、選冶容易,是目前全球60%以上Ni金屬的來源礦床。世界范圍內(nèi)Ni和PGE資源的空間分布極不均勻,巖漿硫化物型鎳礦主要集中在加拿大(Sudbury和Voisey’s Bay礦床)、俄羅斯(Noril’sk礦床)和中國(金川礦床)等國,紅土型鎳礦主要集中于新喀里多尼亞、印度尼西亞和古巴等國,大型層狀雜巖鉑礦主要集中于南非(Bushveld礦床)和美國(Stillwater礦床)?,F(xiàn)有勘探實(shí)踐表明全球共發(fā)現(xiàn)5個(gè)超大型鎳銅硫化物礦床,Ni、Cu和PGE資源量主要集中在這幾個(gè)超大型巖漿礦床中[3],如Noril’sk礦床、Sudbury礦床、Voisey’s Bay礦床和金川礦床等超大型Ni-Cu-PGE硫化物礦床控制著該類礦床的資源量,鎳資源總量估計(jì)高達(dá)3.8×107t,占世界上已知該類礦床資源總量(含PGE)的73%[3,4]。可見Ni-Cu-PGE硫化物礦床超大規(guī)模成礦作用具有重要意義。

    Sudbury Ni-Cu-PGE硫化物礦床由隕石沖擊誘發(fā)的地殼重熔形成[5],其他超大型巖漿Ni-Cu-PGE硫化物礦床是地幔(地幔柱)巖漿作用及相關(guān)的殼幔相互作用的產(chǎn)物[3,4,6]。這些超大型礦床的形成需要巨量的成礦金屬,要求成礦巖漿在地球內(nèi)部較大范圍內(nèi)富集成礦金屬物質(zhì),因此Ni-PGE硫化物礦床常賦礦于大型巖體中,或伴生有同期的巖漿巖。Noril’sk礦床周圍有同時(shí)代的噴出巖和侵入巖。賦礦巖體較小的Voisey’s Bay礦床被認(rèn)為是幔源巖漿體系巖漿通道賦礦[7]。

    金川Ni-Cu-PGE硫化物礦床賦存巖體極小,且區(qū)域內(nèi)沒有發(fā)現(xiàn)同時(shí)代的幔源噴出巖和侵入巖,金川礦床區(qū)域內(nèi)碎屑鋯石也沒有顯示同期的巖漿過程[8],因此難以用大規(guī)模巖漿作用(大巖體)或通道成礦模式解釋,成礦巖漿深部熔離、分期灌入時(shí),硅酸鹽熔體由于密度較小,比硫化物熔體優(yōu)先被擠出巖漿房侵位。碎屑鋯石指示金川巖體應(yīng)該屬于華北克拉通[8],并非構(gòu)造推覆地質(zhì)體,因此推斷金川礦床成礦巖漿的規(guī)??赡懿淮?。小規(guī)模的巖漿作用(小巖體)是否可以形成超大型Ni-PGE硫化物礦床?巨量成礦金屬富集的機(jī)理是什么?

    地質(zhì)觀測(cè)和實(shí)驗(yàn)數(shù)據(jù)表明,流體組分(特別是地幔條件下超臨界態(tài)流體)對(duì)幔源硅酸鹽巖漿形成及成礦元素的遷移和富集具有不可忽視的作用,在Ni-PGE成礦金屬富集過程中可能具有重要作用[9~13]。本文從幔源成礦巖漿作用體系、硫飽和與硫化物熔離等成礦控制因素、流體的作用等方面總結(jié)了Ni-Cu-PGE硫化物礦床形成的機(jī)制,探討小規(guī)模巖漿形成超大型巖漿礦床的可能性。

    2 成礦巖漿作用

    2.1 成礦巖漿類型

    Ni和PGE等為親鐵元素,在地球早期化學(xué)分異演化時(shí)期,原始地球核-幔分異過程中,Ni和PGE高度富集于地核,致使地幔的Ni和PGE豐度低,比地球平均值和球粒隕石低兩個(gè)數(shù)量級(jí);地殼的Ni和PGE豐度更低[14~16]。Cu是典型的親硫元素,在地幔中Cu的虧損明顯弱于Ni和PGE,其地殼豐度高于地幔豐度。但部分幔源鎂鐵質(zhì)火成巖中Ni、Cu和PGE的豐度卻能達(dá)到地殼平均值的數(shù)倍至數(shù)十倍,形成巖漿礦床[17,18],因此,幔源巖漿系統(tǒng)是Ni-Cu-PGE硫化物礦床形成的重要體系。

    幔源巖漿的Ni和PGE組成與含量主要取決于地幔源區(qū)的成分、Ni和PGE的賦存狀態(tài)以及部分熔融程度[4,19]。Ni在硫化物與硅酸鹽巖漿的分配系數(shù)遠(yuǎn)高于在橄欖石與硅酸鹽巖漿的分配系數(shù),Ni、Cu和PGE優(yōu)先賦存于地幔源區(qū)硫化物中,幔源巖漿中Ni、Cu和PGE的含量取決于硫化物是否被全部熔融出來進(jìn)入成礦巖漿中[19,20]。但地幔中的硫化物含量有限,因而地幔中大量的Ni主要賦存于橄欖石中,成礦巖漿中Ni的含量受控于難熔礦物橄欖石的熔融量,即部分熔融程度[2,20,21]。部分熔融程度高的科馬提巖和苦橄巖具有較高的Ni和PGE含量,而地核或核幔邊界起源的地幔柱巖漿作用中Ni和PGE的含量較高,如Bushveld、Noril’sk和Stillwater礦床。源區(qū)Ni含量低且部分熔融程度較低的大洋中脊玄武巖(MORB)具有非常低的PGE含量[19]。因此高鎂巖漿是Ni-Cu-PGE硫化物礦床成礦巖漿的主要類型[3,4,22,23]。

    巖漿Ni-Cu-PGE硫化物礦床成礦巖漿類型主要有M科馬提質(zhì)巖漿和玄武質(zhì)巖漿[21],成分變化較大[3,4,24]??岂R提質(zhì)巖漿是一些重要硫化物礦床的母巖漿,相關(guān)礦床形成局限于太古代(如Kambalda礦床)和古元古代(如Thompson和Raglan等礦床)[3],硫化物中Ni含量最高,具有最高的Ni/Cu值和最低的Pd/Ir值。拉斑玄武質(zhì)巖漿形成的礦床中的Ni-PGE資源量明顯高于科馬提質(zhì)巖漿(見圖1),主要發(fā)育于克拉通地區(qū),礦化不如科馬提質(zhì)巖漿形成的礦床那么普遍,但卻形成了重要的礦床[3,4]。

    圖1 不同巖漿類型的Ni-Cu-PGE硫化物礦床的PGE、Co、Cu與Ni礦石量Fig.1 Ore tonnages of PGE,Co,Cu and Ni from different types of magma in magmatic Ni-Cu-PGE sulfide deposits

    科馬提質(zhì)巖漿和拉斑玄武質(zhì)巖漿形成的Ni-Cu-PGE硫化物礦床中礦石的PGE配分型式不同,(Pt+Pd)/(Os+Ir+Ru)值有所差異,科馬提質(zhì)巖漿礦床具平坦的PGE配分型式,(Pt+Pd)/(Os+Ir+Ru)= 3~3.5。拉斑玄武質(zhì)巖漿礦床具有Pt-Pd富集的PGE配分型式,(Pt+Pd)/(Os+Ir+Ru)=5.6~55.6[11,25]。

    2.2 成礦巖漿(巖體)規(guī)模

    幔源硅酸質(zhì)巖漿(科馬提質(zhì)巖漿和玄武質(zhì)巖漿)中Ni-PGE的溶解度有限[14,21],因此超大型Ni-PGE硫化物礦床的形成一般需要大規(guī)模的巖漿作用以聚集巨量的成礦金屬作為成礦物質(zhì)。賦礦巖體有大型(層狀)巖體、通道及小巖體等不同的類型[23,24,26],成礦巖漿(巖體)規(guī)模有較大差異。

    2.2.1 大型(層狀)巖體中的礦床

    大型(層狀)巖體產(chǎn)出礦床的成礦巖漿規(guī)模巨大,常由多個(gè)小礦床組成礦群。例如,Noril’sk礦床和Sudbury礦床是世界上兩個(gè)最大的Ni富集區(qū),是由多個(gè)小礦床組成的礦群[4]。南非的Bushveld礦床(2.06 Ga)[27]、美國的Stillwater礦床(2.7 Ga)[28]和Duluth礦床(1.1 Ga)[29]賦存于大型層狀巖體中,賦存鉑族金屬礦床,同時(shí)蘊(yùn)含有巨量的鎳硫化物礦石儲(chǔ)量。

    加拿大安大略湖北部山脈Sudbury Ni-Cu硫化物礦床賦存于1.85 Ga的層狀火成雜巖體[3]中,成礦巖體侵位于云英閃長片麻巖與石英二長巖間的接觸帶上,是隕石巨大撞擊導(dǎo)致的大規(guī)模巖漿活動(dòng)的產(chǎn)物,具獨(dú)特的巖石組合:深部為基性-超基性巖石,基底為石英蘇長巖,中間為石英輝長巖,近地表由長英質(zhì)巖石組成,頂部為花崗巖帽,賦礦巖石為蘇長巖-輝長巖[3]。Sudbury礦床由多個(gè)小的礦床群組成,不同礦床中Ni、Cu和PGE等金屬元素含量差異大,具有垂直分帶特征。

    Noril’sk-Talnakh Ni-PGE硫化物礦床位于西伯利亞地臺(tái)西北部,賦存在二疊紀(jì)末期西伯利亞地幔柱巖漿作用相關(guān)的大型超鎂鐵質(zhì)巖體中。玄武巖鋯石U-Pb年齡為251.2±0.3 Ma[30]。鎂鐵-超鎂鐵質(zhì)侵入巖主要由橄欖輝長巖、磁鐵輝長巖、斑雜輝長巖和輝長閃長巖組成,侵位于泥盆紀(jì)沉積巖和大火成巖省最底部的大陸溢流玄武巖中[3]。巨大致密的硫化物礦體賦存于苦橄巖-輝長巖-輝綠巖巖床,輝長巖-輝綠巖中具有圓形硫化物析離體的滴狀礦石(下部為鎳黃鐵礦-黃鐵礦帶,上部為黃銅礦帶),礦石具海綿隕鐵結(jié)構(gòu)和填隙結(jié)構(gòu),在輝長巖-輝綠巖和接觸變質(zhì)圍巖中疊加了硫化物礦化作用形成的角礫狀和細(xì)脈狀礦石。Noril’sk礦床的資源量可能需要多達(dá)200 km3的成礦巖漿[24]。

    2.2.2 ?。ㄍǖ佬停r體賦存礦床

    小巖體及通道型巖體產(chǎn)出礦床的賦礦巖體較小,如Voisey’s Bay礦床和金川礦床。Voisey’s Bay礦床是20世紀(jì)末發(fā)現(xiàn)的超大型巖漿Ni-PGE硫化物礦床,為巖漿通道賦礦,賦礦巖體規(guī)模小,成礦巖漿規(guī)模可能巨大[7,31]。金川Ni-Cu-PGE硫化物礦床由一個(gè)單獨(dú)的大礦床組成,賦存世界第三的Ni-Cu資源量,同時(shí)富集PGE,但賦礦巖體極小[32]。

    Voisey’s Bay Ni-Cu硫化物礦床產(chǎn)于Nain深成巖系(NPS,1350~1290Ma)的橄長巖席(1334Ma)中,侵位于東太古代Nain省和NNW向Torngat造山帶(1 850~1 730 Ma)交界部位[33]。Voisey’s Bay礦床成礦巖漿在早期遭受中下地殼太古宙巖石的混染,晚期則遭受古元古代上地殼Tasiuyak片麻巖的混染[34]。

    金川Ni-Cu-PGE硫化物礦床位于阿拉善地塊西南緣的龍首山隆起,呈NW-SE向不整合侵位于前長城系白家嘴子組大理巖、片麻巖及少量混合花崗巖中。含礦超鎂鐵質(zhì)巖體主要由富橄欖石和鉻鐵礦的二輝橄欖巖、斜長二輝橄欖巖、橄長巖、橄欖二輝巖和二輝巖組成[32]。超鎂鐵質(zhì)侵入體鋯石-斜鋯石U-Pb年齡為831 Ma[35]。金川Ni-Cu硫化物礦床賦礦巖體極小,硫化物與硅酸鹽比例極高[36]。

    2.3 成礦巖漿過程

    地幔高程度部分熔融形成的成礦巖漿演化過程中,經(jīng)歷巖漿結(jié)晶分異、不混熔硫化物熔漿與硅酸鹽巖漿熔離等巖漿過程,Ni及PGE硫化物得以富集成礦。賦存于鎂鐵-超鎂鐵質(zhì)雜巖中的Ni-PGE硫化物礦床主要形成于上述兩種巖漿作用過程。

    1)結(jié)晶分異作用:巖漿就地結(jié)晶過程中堆積作用形成了鉑族金屬礦床及鉻鐵礦、鈦鐵礦和磁鐵礦礦床,如Bushveld礦床、Stillwater礦床和Duluth鉑族金屬礦床的大型層狀雜巖具層狀堆積特征。南非Bushveld雜巖體中Merensky礦層是一個(gè)PGE富集層位(占總資源量的30%),硫化物是PGE的主要聚集體。該礦層由兩個(gè)或兩個(gè)以上含硫化物的基性巖漿上升匯聚形成50 cm至數(shù)米厚的超鎂鐵質(zhì)堆晶巖,即Merensky輝石巖,主要是斜方輝石巖,局部包括橄欖巖。如果沒有硫化物的早期熔離,則巖漿分異結(jié)晶過程使成礦金屬分散進(jìn)入巖漿礦物,并發(fā)生分異。

    2)硫化物熔離作用:由于Ni、Cu、PGE等具有很高的硫化物-硅酸鹽巖漿分配系數(shù),硫化物的熔離使成礦物質(zhì)進(jìn)入硫化物富集[19,20]。地幔條件下幔源成礦巖漿中的硫不飽和,隨巖漿侵位上升過程中,壓力減小促使巖漿中硫的溶解度增加,因而巖漿中的硫仍處于不飽和狀態(tài)。因此要使硫化物從硅酸鹽巖漿中熔離出來形成Ni-Cu-PGE礦床,需要其他因素促使硫飽和[4,24,36],如Noril’sk、金川和Voisey’s Bay礦床成礦過程中地殼物質(zhì)混染促使巖漿中硫達(dá)到飽和。因此,硫飽和是巖漿中成礦物質(zhì)富集、硫化物礦床形成的主要控制因素。

    3 成礦作用控制因素

    在Ni-PGE硫化物礦床超大規(guī)模成礦過程中,成礦母巖漿的組成中成礦金屬的大量聚集及巖漿演化系統(tǒng)中硫飽和、硫化物熔離、Ni-Cu-PGE元素富集是礦床形成的兩個(gè)重要過程,具不同的控制因素。

    3.1 成礦物質(zhì)聚集

    巖漿型Ni、Cu和PGE硫化物礦床形成的物質(zhì)基礎(chǔ)是成礦金屬元素及硫。與硅酸鹽巖漿平衡的不混熔硫化物熔漿富集Ni、Cu和PGE金屬,Ni-Cu-PGE硫化物礦床的形成是成礦巖漿演化過程中不混溶硫化物的熔離聚集、在特定部位(底部、通道等)濃集、堆積的過程[3,4,21,31]。

    3.1.1 成礦金屬

    在地幔巖漿作用體系中Ni、Cu和PGE的地球化學(xué)行為在硫含量飽和與否的條件下有所差異。在巖漿體系硫不飽和的條件下,Ni、Os、Ir和Ru為相容元素,而Cu和Pd是強(qiáng)不相容元素,巖漿分離結(jié)晶過程中二者發(fā)生分異。在體系硫飽和的條件下,Ni、Cu和PGE為強(qiáng)不相容元素,硅酸鹽巖漿與硫化物分配系數(shù)低[19,25,37]、特別是PGE具有極高的硫化物熔漿/硅酸鹽熔漿分配系數(shù),且受R因子(硅酸鹽熔漿與硫化物熔漿的質(zhì)量比值)的控制[38,39],相關(guān)元素進(jìn)入硫化物熔漿。

    Ni、Cu和PGE在地幔部分熔融和幔源巖漿分離結(jié)晶過程中地球化學(xué)性質(zhì)差異導(dǎo)致其發(fā)生分異。Ni、Ir和Ru對(duì)于橄欖石、輝石和尖晶石是相容元素,Rh是橄欖石的相容元素,Pt是尖晶石和輝石的相容元素,特別是Ir、Ru和Rh的DSp/Sil可達(dá)102~103數(shù)量級(jí),而Pd和Cu在這些礦物中為不相容元素,分離結(jié)晶過程中總是殘留在巖漿中[2]。在某些層狀巖體的鉻鐵礦層(如南非Bushveld巖體的UG1和UG2層)或蛇綠巖帶的鉻鐵礦礦床(如西藏羅布莎)中可能形成銥鉑族(IPGE:Ir、Os、Ru)的富集。

    除隕石沖擊成因的Sudbury超大型Ni-Cu硫化物礦床為陸殼變質(zhì)火山巖和沉積巖熔融物提供成礦物質(zhì)外,其他巖漿Ni-Cu硫化物礦床的成礦物質(zhì)主要來源于地幔。地幔是地殼PGE富集的主要源庫,大陸拉斑玄武巖中的Ni含量為99 ppm(1 ppm= 1×10-6),Cu為90 ppm。PGE遷移主要有兩個(gè)途徑:地幔部分熔融產(chǎn)生的巖漿侵入地殼,或地幔板片就位于俯沖/碰撞帶。

    3.1.2 硫

    地球演化過程的化學(xué)分異產(chǎn)生地幔不均一性,不同地幔源區(qū)的硫含量也有很大差異,地幔中的硫含量一般僅為5~400 ppm,平均為200 ppm[40]。 在地幔壓力下,科馬提質(zhì)與玄武質(zhì)巖漿中硫的溶解度為500~1 000 ppm[41,42],是巖漿溫度、壓力、氧逸度(fO2)和硫逸度(fS2)的函數(shù)[43,44]。

    1)溫度:在地幔條件下,巖漿中硫的飽和度與溫度呈正相關(guān)。硫不飽和巖漿體系中,隨溫度的升高,硅酸鹽熔體中硫的溶解度增大,在巖漿溫度為1 000℃時(shí),溫度每增加100℃,硫溶解度增長6倍,當(dāng)巖漿溫度為1 400℃時(shí),硫溶解度增長3倍[45]。在溫度為1 200℃、fO2=10-10.5時(shí),鐵硫化物達(dá)到飽和,含10%(質(zhì)量分?jǐn)?shù))FeO的拉斑玄武質(zhì)巖漿可以溶解大約0.04%的硫[42],溫度每增加100℃就可以增加溶解0.05%~0.09%的硫。

    2)壓力:巖漿中硫的飽和度隨壓力的升高呈指數(shù)減小,硫化物和硅酸鹽熔體間的可混溶性隨著壓力的升高而增大[41]。

    3)氧逸度:氧逸度對(duì)硫飽和度的影響相對(duì)小,主要影響硫的存在形式。在地幔巖漿溫度(1 400~1 500℃)范圍內(nèi),巖漿的氧逸度大于10-6大氣壓時(shí),硫以SO42-的形式溶解;氧逸度小于10-6大氣壓時(shí),硫則以硫化物的形式存在。

    4)巖漿成分:硫飽和度隨著熔體成分而變化,F(xiàn)eO、MgO及CaO組分的增加促使硫溶解度增加,而SiO2和Al2O3的作用相反。

    巖漿上升侵位、演化分異、與地殼圍巖同化混染等物理化學(xué)條件的變化都會(huì)引起巖漿中硫溶解度的改變。例如,層狀侵入體結(jié)晶分異過程中,在常壓下,巖漿溫度、FeO含量和氧逸度對(duì)硅酸鹽巖漿中硫化物的溶解度至關(guān)重要。在封閉條件下,當(dāng)尖晶石沒有大量結(jié)晶時(shí),巖漿冷卻和結(jié)晶遵從緩沖曲線,此時(shí)氧逸度不是硫化物溶解度的主要控制因素,巖漿溫度和FeO含量的下降會(huì)導(dǎo)致硫化物溶解度降低。橄欖石結(jié)晶過程中,巖漿溫度和FeO含量的下降會(huì)使硫化物溶解度降低。古銅輝石結(jié)晶的過程中,溫度連續(xù)快速下降,但FeO含量相對(duì)穩(wěn)定,以后隨FeO/MgO值的增加而增加。隨斜長石的結(jié)晶,巖漿的FeO含量開始增加,影響硫化物溶解度連續(xù)下降的速度。

    3.2 硫化物熔離作用的控制因素

    硫化物熔體的熔離是Ni-Cu-PGE硫化物礦床形成的關(guān)鍵過程,而硫化物熔體熔離的關(guān)鍵在于成礦巖漿中硫的過飽和,受巖漿體系中硫含量的影響[21,42]。Vogt[46]提出巖漿硫化物礦床的不混熔成因機(jī)制——巖漿熔離說,即原始含礦巖漿在結(jié)晶重力分異之前發(fā)生硫化物液態(tài)分異作用。硫化物的熔離取決于巖漿中硫的溶解度、氧化還原電位以及巖漿的溫度和粘度等眾多因素。

    地幔硅酸鹽巖漿中硫的飽和度與溫度正相關(guān),與壓力負(fù)相關(guān)。成礦巖漿上侵至地殼巖漿房時(shí),壓力下降,硫的溶解度增高,大多數(shù)巖漿源區(qū)硫不飽和,即使硫飽和的成礦巖漿上升后也將不再飽和[41,47],因此促使巖漿中硫達(dá)到飽和需要其他因素,包括巖漿快速冷卻、巖漿混合、結(jié)晶分異作用、外來硫的加入或地殼混染等因素[3,4,24,48,49]。

    1)巖漿快速冷卻。隨溫度的降低巖漿中硫的溶解度降低,巖漿快速冷卻可促使硫達(dá)到飽和,如Sudbury礦床[50]。

    2)巖漿混合。巖漿混合可降低巖漿的液相線溫度,成分差異較大的巖漿混合可改變硫的飽和曲線,使硫進(jìn)入飽和區(qū)間[24]。

    3)巖漿分離結(jié)晶。分離結(jié)晶作用的不斷進(jìn)行促使幔源巖漿中硫富集,隨橄欖石、輝石等造巖礦物的不斷結(jié)晶,成礦巖漿中不相容元素硫的濃度不斷增高達(dá)到飽和。巖漿經(jīng)歷60%以上的分離結(jié)晶作用后達(dá)到硫飽和[24,41,51]。

    4)外來硫的加入。富硫圍巖被巖漿同化,增加硫含量促使硫飽和;外來硫的補(bǔ)充是超大型Ni-Cu-PGE硫化物礦床形成的關(guān)鍵因素[49]。例如,加拿大的Thompson礦石的硫同位素組成指示硫來源于海水碳酸鹽,俄羅斯Noril’sk硫化物礦床形成中硫飽和是圍巖硫加入巖漿[49,52]。源于圍巖硫的同化作用促使幔源巖漿硫飽和[41,51]。

    5)富Si圍巖混染。SiO2的加入可降低巖漿中硫的溶解度。地殼物質(zhì)的污染、高硅物質(zhì)加入鎂鐵質(zhì)巖漿,改變了巖漿中四面體和八面體的比例,也改變了巖漿中硫的溶解度,硫飽和促使硫與金屬結(jié)合形成硫化物[53]。Bushveld、Voisey’s Bay Ni-Cu硫化物礦床中O、Os、Nd、Sr和Pb同位素體系均表明強(qiáng)烈的地殼物質(zhì)的混染促使硫化物的熔離[53,54]。

    Ni-Cu-PGE硫化物礦床超大規(guī)模成礦作用中,需要大量的硫促使硫化物從硅酸鹽熔體中熔離。巖漿快速冷卻和巖漿混合在理論上是完全可行的,但相關(guān)巖漿過程的地質(zhì)實(shí)際較少。巖漿分離結(jié)晶作用可促使硫飽和,但易使成礦金屬分散,難以形成品位較高的礦床,如美國的Stillwater礦床和Duluth鉑族金屬礦床。

    3.3 超大規(guī)模成礦的構(gòu)造環(huán)境因素

    超大型Ni-Cu-PGE硫化物礦床的形成一般需要下列關(guān)鍵因素:富Mg巖漿,大陸地殼環(huán)境,圍巖中的硫源,同期巖體存在親銅元素虧損,與圍巖相互反應(yīng)和巖漿通道[3,4]。在構(gòu)造環(huán)境上超大型Ni-Cu硫化物礦床主要形成于裂谷、地幔柱等張性環(huán)境中[6,55]。成礦巖漿通道系統(tǒng)是Ni-Cu-PGE硫化物礦床形成聚集的重要環(huán)境。例如,Noril’sk-Talnakh、Voisey’s Bay等礦床均發(fā)現(xiàn)于較小的通道型巖體中[7,49]。在不同構(gòu)造控制的部位形成的巖體形態(tài)各異[31]。

    4 巖漿流體揮發(fā)分與成礦作用

    現(xiàn)有巖漿硫化物礦床成礦理論建立在硅酸鹽巖漿溶解成礦金屬、硫化物熔離富集成礦的基礎(chǔ)上,受成礦金屬元素和硫在幔源硅酸鹽巖漿中溶解度的限制[14~16,21,42],需要大規(guī)模的巖漿才足以聚集巨量的成礦物質(zhì),即超大型巖漿硫化物礦床的形成要求大型(層狀)巖體或通道巖體。由于成礦巖漿硫飽和、硫化物熔離富集成礦金屬在地殼巖漿房不難實(shí)現(xiàn),所以小規(guī)模巖漿作用(即小巖體)能否形成超大型礦床的關(guān)鍵問題和必要條件是巨量成礦金屬元素聚集的可能性。地幔條件下超臨界流體(特別是CO2)對(duì)金屬元素(特別是納米級(jí)金屬元素)的超強(qiáng)搬運(yùn)能力提供了成礦潛在的可能性[10~12,56]。

    地幔部分熔融形成巖漿的過程中流體揮發(fā)分的助熔作用得到模擬實(shí)驗(yàn)的證實(shí)和廣泛認(rèn)同,近年來的實(shí)踐觀察和實(shí)驗(yàn)確認(rèn)流體揮發(fā)分對(duì)成礦金屬元素具有淋濾、遷移和聚集作用[10~13,56]。

    4.1 巖漿流體揮發(fā)分對(duì)金屬元素的遷移富集

    現(xiàn)有成礦金屬元素氣液遷移實(shí)驗(yàn)研究側(cè)重于水熱溶液及固-液反應(yīng)[57],對(duì)氣相遷移的研究較少。許多金屬元素的遷移、富集和成礦過程中氣相(特別是超臨界流體)起著重要作用。地幔排氣、巖漿去氣、火山噴氣及海底噴流中的各種地球內(nèi)部流體揮發(fā)分是元素遷移和富集的重要介質(zhì)[58~60]。在火山噴氣和巖漿礦床的氣體包裹體中保留著金屬元素氣相遷移的痕跡。

    很多金屬和金屬絡(luò)合物在巖漿條件下呈高揮發(fā)性氣相運(yùn)移,親銅元素及Au、Ir和Pt等貴金屬元素高度富集于火山氣體和海底噴流中[10,11,61~63]。隨著巖漿成分與溫度的不同,火山噴氣中金屬元素的濃度變化很大,約n×10-9~n×10-6[57,62,64,65]。玄武質(zhì)巖漿(尤其是堿性玄武質(zhì)巖漿)噴氣中Cu、Zn、Pb、As、Ag和Au的濃度較高。深部來源的高溫火山氣體攜帶的大量金屬元素,隨溫度降低在火山噴氣孔及附近析凝為固態(tài)的礦物集合體(升華殼)和凝結(jié)水,其中富集Cu、Pb、Ag和Au等多種成礦金屬[59,62,64,66]。例如,從夏威夷Kilauea火山噴氣的凝結(jié)物中檢測(cè)到Co、V、Cr、Cu、Au、Ag、In、Re和Ir等元素,Kilauea火山氣相中PGE Ir的高富集與該區(qū)玄武質(zhì)巖漿的深部來源及噴氣的高F、Cl含量有關(guān)[61,67]。凝結(jié)物中Cu等元素相對(duì)于母巖漿具很高的富集因子(即富集系數(shù)=凝結(jié)物元素濃度/母巖漿元素濃度),如Cu為10~102,Zn、W、Ir和Ag為102~105,Au、Cd、Re和S則高于105,且最高富集因子見于最富氣體的噴發(fā)期內(nèi)[67]。印度尼西亞Merapi火山液態(tài)的凝結(jié)水中Ni、Cr、Sb、Cu、Pb的富集因子為5~35[68]。

    金屬呈簡(jiǎn)單氯化物、氟化物、硫化物、含氧酸鹽及原子遷移的理論假設(shè)低估了氣體搬運(yùn)金屬的能力,基于實(shí)驗(yàn)數(shù)據(jù)的熱力學(xué)模擬計(jì)算表明,水熱氣相可運(yùn)移足夠的成礦金屬以形成具有經(jīng)濟(jì)價(jià)值的淺成熱液型和斑巖型礦床[56,69,70]。某些斑巖型Cu-Au礦床的富氣體包裹體比共存的鹽水包裹體更富集Cu和Au[69,71]。如,Grasberg斑巖銅礦床中富氣相包裹體的Cu和Au含量平均為1.2%和10.17 ppm,而共存的鹽水(液相)包裹體的Cu和Au平均含量為0.3%和0.26 ppm[71]。

    流體揮發(fā)分在地球分異過程中對(duì)親銅元素在硅酸鹽和硫化物熔體間的分異具有重要影響[16,72]。金屬元素在蒸汽相中的溶解實(shí)驗(yàn)研究表明Cu、Au和As等元素在富含硫的巖漿流體相分離時(shí)優(yōu)先進(jìn)入蒸汽相[70]。

    火山噴氣的凝結(jié)物中高濃度的Cu、Zn、Pb、As、Ag和Au,以及斑巖型礦床中低密度流體包裹體(蒸氣相)中硫化物礦物的存在,預(yù)示著這些金屬元素可能是以氣相搬運(yùn)的。

    4.2 地幔環(huán)境流體揮發(fā)分對(duì)金屬的搬運(yùn)作用

    地幔條件下超臨界流體對(duì)Cu、Ni和PGE(特別是納米級(jí)PGE)的搬運(yùn)聚集具有重要的作用。超臨界CO2流體對(duì)Cu、Ni和Au等金屬的溶解、搬運(yùn)能力巨大,并在適當(dāng)?shù)沫h(huán)境下富集[9,12,13,59,70,73~75]。例如,Cu在巖漿氣相中的濃度比共存熔體中高幾百倍[67,73]。

    地幔捕虜體普遍存在Cu、Ni和PGE合金或微粒,南非Bushveld雜巖Merenskey礦層橄欖石、鉻鐵礦及金屬硫化物中發(fā)現(xiàn)直徑<100 nm的Cu、Ni和PGE的微粒包體[76,77],揭示金屬元素以微?;蚴╟lusters)的形式通過幔源巖漿而遷移[78],表明Cu、Ni和PGE在地球深部超臨界態(tài)地幔流體上涌過程中被淋濾、搬運(yùn)大量注入成礦系統(tǒng),為成礦作用的持續(xù)進(jìn)行和巨量成礦元素堆積提供了有利條件。由此可見超臨界流體對(duì)金屬成礦物質(zhì)的搬運(yùn)與聚集對(duì)超大型礦床的形成具有重要的意義。

    流體化學(xué)組成、硫的含量對(duì)巖漿的氧化還原狀態(tài)和硫飽和具有重要的影響[79]。Cu、Ag和Au在硅酸鹽熔體+蒸汽+鹵水±礦物(磁鐵礦±磁黃鐵礦)體系中的分配系數(shù)在無硫和含硫體系中有較大差別,含硫體系明顯較高[80]。Lowenstern[74]認(rèn)為含有CO2的巖漿流體在中-上地殼環(huán)境下普遍存在不混溶現(xiàn)象。Au在富硫巖漿流體相分離時(shí)優(yōu)先進(jìn)入蒸汽相,因此CO2的存在對(duì)成礦流體中Au的再分配與富集具有重要的作用。

    4.3 巖漿Ni-Cu硫化物礦床中流體的作用

    Ni-Cu-PGE硫化物礦床成礦巖漿體系中流體的化學(xué)組成隨構(gòu)造環(huán)境不同有所差異[81~84]??死ōh(huán)境的金川Ni-Cu-PGE硫化物礦床成礦巖漿流體主要由CO2和SO2組成,稀有氣體同位素表明成礦巖漿起源于地幔來源的流體與殼源流體混合的環(huán)境,演化過程中混入大氣來源再循環(huán)組分[85]。中亞造山帶新疆北部喀拉通克和東天山Ni-Cu-PGE硫化物礦床成礦巖漿流體揮發(fā)分以H2O為主(平均為92%),其次為CO2、H2、H2S和SO2,成礦巖漿流體主要以H2(黃山)、H2O(喀拉通克)為主,橄欖石和輝石中高的H2和H2O揮發(fā)分相表明大量流體可能與再循環(huán)板片脫水相關(guān)[81,83]。地幔柱大火成巖省巖漿Ni-Cu硫化物礦床成礦巖漿流體主要由H2O、CO2和SO2組成[82,84]。

    中亞造山帶喀拉通克硫化物礦床中CO2與成礦金屬Cu和Ni含量呈正比[81](圖2a、圖2b),地幔柱大火成巖省朱布Ni-Cu硫化物礦床CO2與成礦金屬Cu和Ni含量沒有相關(guān)性(圖2c、圖2d)。結(jié)果表明造山帶環(huán)境成礦巖漿體系中高溫高壓下超臨界態(tài)的CO2參與搬運(yùn)與聚集金屬成礦物質(zhì),可能是大量的Cu和Ni等金屬聚集的關(guān)鍵因素;但地幔柱大規(guī)模成礦巖漿體系中Cu和Ni等金屬聚集可能與流體揮發(fā)分關(guān)系不明顯。

    圖2 喀拉通克和朱布Ni-Cu硫化物礦床CO2與Cu和Ni含量相關(guān)圖(數(shù)據(jù)來自文獻(xiàn)[83,84])Fig.2 The correlation diagram of CO2vs.Cu and Ni contents in Karatongke and Zhubu Ni-Cu sulfide deposits (data from Ref.[83,84])

    俄羅斯Dovyren層狀雜巖體中含有浸染狀硫化物,流體包裹體中CH4和H2等含量峰與PGE和Au層位一致[86]。Sudbury雜巖體中接觸帶型Cu-PGEAu礦床以CH4為主要成分的碳?xì)浠衔飳?duì)成礦金屬元素起搬運(yùn)作用,Cu、Au、B和Ag在CH4包裹體與鹵水包裹體中的分配系數(shù)遠(yuǎn)大于Na、Ca、Fe、Mn、Zn和Pb在兩相中的分配系數(shù),兩種流體發(fā)生不混溶時(shí),碳?xì)浠衔锪黧w中更富集成礦金屬。碳?xì)浠衔锏蓉毸畵]發(fā)分從硫化物熔體中的出溶最終決定了Au和PGE在成礦系統(tǒng)中的分配。足夠證據(jù)表明巖漿礦床成礦巖漿起源于流體富集的地幔源區(qū)[23]。

    5 結(jié)語

    1)超大型巖漿Ni-Cu-PGE硫化物礦床的形成取決于地幔巖漿作用中Ni、Cu和PGE元素在成礦巖漿中的聚集數(shù)量以及巖漿演化過程中硫飽和與硫化物熔離的程度。高度部分熔融和流體對(duì)成礦金屬的運(yùn)移富集具有重要作用。

    2)大規(guī)模的高度部分熔融(高M(jìn)g)促使地幔硫化物及橄欖石中巨量的Ni和PGE進(jìn)入巖漿、硫化物熔離聚集形成大型層狀(或通道)巖體賦存的超大型Ni-Cu-PGE硫化物礦床。

    3)地幔超臨界流體可能大范圍內(nèi)搬運(yùn)聚集大量Ni-PGE金屬元素,對(duì)巖漿硫化物礦床成礦可能有重要貢獻(xiàn),是小規(guī)模巖漿(小巖體)成礦的可能機(jī)制。

    致謝:特將本文獻(xiàn)給湯中立教授從事地質(zhì)工作60年暨80華誕。李文淵和張銘杰攜同事和弟子祝賀湯中立院士80華誕,感謝湯老師多年來的指導(dǎo)與關(guān)懷。在本文完成過程中Chusi Li、焦建剛、張照偉及匿名評(píng)審專家等給予了指導(dǎo)和幫助,在此表示衷心的感謝。

    [1] 李文淵.巖漿Cu-Ni-PGE礦床研究現(xiàn)狀及發(fā)展趨勢(shì)[J].西北地質(zhì),2007,40(2):1-18.

    [2] 宋謝炎,胡瑞忠,陳列錳.銅、鎳、鉑族元素地球化學(xué)性質(zhì)及其在幔源巖漿起源、演化和巖漿硫化物礦床研究中的意義[J].地學(xué)前緣,2009,16(4):287-305.

    [3] Naldrett A J.Magmatic Sulfide Deposits:Geology,Geochemistry and Exploration[M].Berlin:Springer Verlag,Heidelberg,2004.

    [4] Naldrett A J.World-class Ni-Cu-PGE deposits:Key factors in their genesis[J].Mineralium Deposita,1999,34(3):227-240.

    [5] Dietz R S.Sudbury structure as an astrobleme[J].Journal of Geology,1964,72(4):412-434.

    [6] Lambert D D,F(xiàn)oster J G,F(xiàn)rick L R,et al.Geodynamics of magmatic Cu-Ni-PGE sulfide deposits;new insights from the Re-Os isotope system[J].Economic Geology,1998,93(2):121-136.

    [7] Evans-Lamswood D M,Butt D P,Jackson R S,et al.Physical controls associated with the distribution of sulfides in the Voisey’s Bay Ni-Cu-Co deposit,Labrador[J].Economic Geology,2000,95(4):749-769.

    [8] Tang Qingyan,Li Chusi,Zhang Mingjie,et al.Detrital zircon constraint on the timing of amalgamation between Alxa and Ordos,with exploration implications for Jinchuan-type Ni-Cu ore deposit in China[J].Precambrian Research,2014,255:748-755.

    [9] Boudreau A E,McCallum I S.Concentration of platinum-group elements by magmatic fluids in layered intrusions[J].Economic Geology,1992,87(7):1830-1848.

    [10] Fleet M E,Wu T W.Volatile transport of precious metals at 1 000℃:Speciation,fractionation,and effect of base-metal sulfide[J].Geochimica et Cosmochimica Acta,1995,59(3):487-495.

    [11] Fleet M E,Crocket J H,Stone W E.Partitioning of platinumgroup elements(Os,Ir,Ru,Pt,Pd)and gold between sulfide liquid and basalt melt[J].Geochimica et Cosmochimica Acta,1996,60(13):2397-2412.

    [12] Blackburn J M,Long D P,Cabanas A,et al.Deposition of conformal copper and nickel films from supercritical carbon dioxide[J].Science,2001,294(5540):141-145.

    [13] Fiorentini M L,Beresford S W.Role of volatiles and metasomatized subcontinental lithospheric mantle in the genesis of magmatic Ni-Cu-PGE mineralization[J].Terra Nova,2008,20(5):333-340.

    [14] Barnes S J,Picard C P.The behavior of platinum-group elements during partial melting,crystal fractionation,and sulfide segregation[J].Geochimica et Cosmochimica Acta,1993,57(1):79-87.

    [15] Snow J E,Schmidt G.Constraints on Earth accretion deduced from noble metals in the oceanic mantle[J].Nature,1998,391 (6663):166-169.

    [16] Day J M D,Pearson D G,Taylor L A.Highly siderophile element constraints on accretion and differentiation of the Earth-Moon System[J].Science,2007,315(5809):217-219.

    [17] Yuan Feng,Zhou Taofa,Zhang Dayu,et al.Siderophile and chalcophile metal variations in basalts:Implications for the sulfide saturation history and Ni-Cu-PGE mineralization potential of the Tarim continental flood basalt province,Xinjiang Province,China[J].Ore Geology Reviews,2012,45:5-15.

    [18] Kamenetsky V S,Maas R,F(xiàn)onseca R O C,et al.Noble metals potential of sulfide-saturated melts from the subcontinental lithosphere[J].Geology,2013,41(5):575-578.

    [19] Peach C L,Mathez E A,Keays R R.Sulfide melt silicate melt distribution coefficients for noble metals and other chalcophile elements as deduced from MORB[J].Geochimica et Cosmochimica Acta,1990,54(12):3379-3389.

    [20] Bezmen N I,Asif M,Brugmann G E,et al.Distribution of Pd,Rh,Ru,Ir,Os and Au between sulfide and silicate melts [J].Geochimica et Cosmochimica Acta,1994,58(4):1251-1260.

    [21] Keays R R.The role of komatiitic and picritic magmatism and S-saturation in the formation of ore-deposits[J].Lithos,1995,34(1-3):1-18.

    [22] Chai G,Naldrett A J.The Jinchuan ultramafic intrusion:Cumulate of a high-Mg basaltic magma[J].Journal of Petrology,1992,33(2):277-303.

    [23] Arndt N.The lithospheric mantle plays no active role in the formation of orthomagmatic ore deposits[J].Economic Geology,2013,108(8):1953-1970.

    [24] Ripley E M,Li C.Sulfide saturation in mafic magmas is external sulfur required for magmatic Ni-Cu-(PGE)ore genesis[J]. Economic Geology,2013,108(1):45-58.

    [25] Barnes S J,Makovicky E,Makovicky M,et al.Partition coefficients for Ni,Cu,Pd,Pt,Rh and Ir between monosulfide solid solution and sulfide liquid and the formation of compositionally zoned Ni-Cu sulfide bodies by fractional crystallization of sulfide liquid[J].Canadian Journal of Earth Sciences,1997,34 (4):366-374.

    [26] Maier W D,Groves D I.Temporal and spatial controls on the formation of magmatic PGE and Ni-Cu deposits[J].Mineralium Deposita,2011,46(8):841-857.

    [27] Scoates J S,F(xiàn)riedman R M.Precise age of the platiniferous Merensky reef,Bushveld complex,south Africa,by the U-Pb zircon chemical abrasion ID-TIMS technique[J].Economic Geology,2008,103(3):465-471.

    [28] Premo W R,Premo W R,Helz R T,et al.U-Pb and Sm-Nd ages for the Stillwater Complex and its associated sills and dikes,Beartooth Mountains,Montana[J].Geology,1990,18(11):1065-1068.

    [29] Paces James B,Miller Jr James D.Precise U-Pb ages of Duluth Complex and related mafic intrusions,northeastern Minnesota:Geochronological insights to physical,petrogenetic,paleomagnetic,and tectonomagmatic processes associated with the 1.1 Ga Midcontinent Rift System[J].Journal of Geophysical Research:Solid Earth(1978—2012),2012,98(B8):13997-14013.

    [30] Kamo L,Czamanske G K,Krogh T E.A minimum U-Pb age for Siberian flood-basalt volcanism[J].Geochimica et Cosmochimica Acta,1996,60(18):3505-3511.

    [31] Lightfoot P C,Evans-Lamswood D.Structural controls on the primary distribution of mafic-ultramafic intrusions containing Ni-Cu-Co-(PGE)sulfide mineralization in the roots of large igneous provinces[J].Ore Geology Reviews,2015,64:354-386.

    [32] 湯中立,李文淵.金川銅鎳硫化物(含鉑)礦床成礦模式及地質(zhì)對(duì)比[M].北京:地質(zhì)出版社,1995.

    [33] Nicholson S W,Cannon W F,Schulz K J.Metallogeny of the Midcontinent rift system of North America[J].Precambrian Research,1992,58(1-4):355-386.

    [34] Li C,Naldrett A J.Melting ions of gneissic inclusions with enclosing magma at Voisey’s Bay,Labrador,Canada[J].Economic Geology,2000,95(4):801-814.

    [35] Zhang Mingjie,Kamo S L,Li C,et al.Precise U-Pb zirconbaddeleyite age of the Jinchuan sulfide ore-bearing ultramafic intrusion,western China[J].Mineralium Deposita,2010,45(1):3-9.

    [36] Li C,Ripley E M.The Giant Jinchuan Ni-Cu-(PGE)deposit:Tectonic setting,magma evolution,ore genesis and exploration implications[M]//Li C,Ripley E M,eds.Magmatic Ni-Cu and PGE Deposits:Geology,Geochemistry and Genesis:Reviews in Economic Geology(Volume 17).Denver,Colorado:Society of Economic Geologists,Inc,2011:163-180.

    [37] Li C,Barnes S J,Makovicky E,et al.Partitioning of nickel,copper,iridium,rhenium,platinum,and palladium between monosulfide solid solution and sulfide liquid[J].Geochimica et Cosmochimica Acta,1996,60(7):1231-1238.

    [38] Campbell I H,Naldrett A J.The influence of silicate:Sulfide ratio on the geochemistry of the magmatic sulfides[J].Economic Geology,1979,74(6):1503-1506.

    [39] Campbell I H,Naldrett A J,Barnes S J.A model for the origin of the platinum-rich sulfide horizons in the Bushveld and Stillwater complexes[J].Journal of Petrology,1983,24(2):133-165.

    [40] David H E,Lorand J P.Mantle sulfide geobarometry[J].Geochimica et Cosmochimica Acta,1993,57(10):2213-2222.

    [41] Mavrogenes J A,O’Neill H S.The relative effect s of pressure,temperature and oxygen fugacity on the solubility of sulfide in mafic magmas[J].Geochimica et Cosmochimica Acta,1999,63(7/8):1173-1180.

    [42] Haughton D R,Roeder P L,Skinner B J.Solubility of sulfur in mafic magmas[J].Economic Geology,1974,69(4):451-467.

    [43] Ballhaus C,Tredoux M,Spaeth A.Phase relations in the Fe-Ni-Cu-PGE-S system at magmatic temperature and application to massive sulfide ores of the Sudbury Igneous Complex[J].Journal of Petrology,2001,42(10):1911-1926.

    [44] Moretti R,Papale P,Ottonello G.A model for the saturation of C-O-H-S fluids in silicate melts[J].Geological Society,London,Special Publications,2003,213:81-101.

    [45] Buchanan D L,Nolan J.Solubility of sulfur and sulfide immiscibility in synthetic tholeiitic melts and their relevance to Bushveld-Complex[J].The Canadian Mineralogist,1979,17(2):483-494.

    [46] Vogt J H L.Beitrage zur genetischen classification der durch magmatische differentiations processe und der durch pneumatolyse entstandenen erzvorkomme[J].Zeitschrift Prakt Geol,1894,2:391-399.

    [47] Wendlandt R F.Sulfide saturation of basalt and andesite melts at high pressures and temperatures[J].American Mineralogist,1982,67(9210):877-885.

    [48] Li C,Ripley E M.Empirical equations to predict the sulfur content of mafic magma at sulfide saturation and applications to magmatic sulfide deposits[J].Mineralium Deposita,2005,40 (2):218-230.

    [49] Keays R R,Lightfoot P C.Crustal sulfur is required to form magmatic Ni-Cu sulfide deposits:Evidence from chalcophile element signatures of Siberian and Deccan Trap basalts[J].Mineralium Deposita,2010,45(3):241-257.

    [50] Maier W D,Barnes S J,De Waal S A.Exploration for magmatic Cu-Ni-PGE sulphide deposits;a review of recent advances in the use of geochemical tools,and their application to some South African ores[J].South African Journal of Geology,1998,101(3):237-253.

    [51] Ripley E M,Li C,Shin D.Paragneiss assimilation in the genesis of magmatic Ni-Cu-Co sulfide mineralization at Voisey’s Bay,Labrador:δ34S,δ13C,and Se/S Evidence[J].Economic Geology,2002,97(6):1307-1318.

    [52] Ripley E M,Lightfoot P C,Li C,et al.Sulfur isotopic studies of continental flood basalts in the Noril’sk region:Implications for the association between lavas and ore-bearing intrusions[J]. Geochimica et Cosmochimica Acta,2003,67(15):2805-2817.

    [53] Maier W D,Arndt N T,Curl E A.Progressive crustal contamination of the Bushveld Complex[J].Contributions to Mineralogy and Petrology,2000,140(3):316-327.

    [54] Lambert D D,F(xiàn)rick L R,F(xiàn)oster J G,et al.Re-Os isotopic systematics of the Voisey’s Bay Ni-Cu-Co magmatic sulfide system,Labrador,Canada:II.Implications for parental magma chemistry,ore genesis,and metal redistribution[J].Economic Geology,2000,95(4):867-888.

    [55] Li X H,Su L,Chung S L,et al.Formation of the Jinchuan ultramafic intrusion and the world’s third largest Ni-Cu sulfide deposit:Associated with the~825 Ma south China mantle plume?[J].Geochemistry,Geophysics,Geosystems,2005,6 (11):16-32.

    [56] Philips G N,Evans K A.Role of CO2in the formation of gold deposits[J].Nature,2004,429(6994):860-863.

    [57] Hanley J J,Mungall J E,Pettke T,et al.Ore metal redistribution by hydrocarbon-brine and hydrocarbon-halide melt phases,North Range footwall of the Sudbury Igneous Complex,Ontario,Canada[J].Mineralium Deposita,2005,40(3):237-256.

    [58] Hedenquist J W,Aoki M,Shinohara H.Flux of volatiles and ore-forming metals from the magmatic-hydrothermal system of Satsuma Iwojima volcano[J].Geology,1994,22(7):585-588.

    [59] Simmons S F,Brown K L.Gold in magmatic hydrothermal solutions and the rapid formation of a giant ore deposit[J].Science,2006,314(5797):288-291.

    [60] 張 生,熊小林,Seward T M.成礦元素的氣相遷移與實(shí)驗(yàn)研究[J].地學(xué)前緣,2009,16(1):68-75.

    [61] Olmez I,F(xiàn)innegan D L,Zoller W H.Iridium emissions from Kilauea Volcano[J].Journal of Geophysical Research:Solid Earth,1986,91(B1):653-663.

    [62] Gemmell B J.Geochemistry of metallic trace elements in furmarolic condensates from Nicaragua and Costa Rican volcanos[J]. Journal of Volcanology and Geothermal Research,1987,33(1-3):161-181.

    [63] Rubin K H.Degassing of metals and metalloids from erupting seamount and mid-ocean ridge volcanoes[J].Geochimica et Cosmochimica Acta,1997,61(17):3525-3542.

    [64] Giggenbach W F.Magma degassing and mineral deposition inhydrothermal systems along convergent plate boundaries[J]. Economic Geology,1992,87(7):1927-1944.

    [65] Taran Y A,Bernard A,Gavilanes J C,et al.Native gold in mineral precipitates from high-temperature volcanic gases of Colima volcano,Mexico[J].Applied Geochemistry,2000,15(3):337-346.

    [66] Quisefit J P,Toutain J P,Bergametti G,et al.Evolution versus cooling of gaseous volcanic emissions from Momtombo volcano,Nicaragua[J].Geochimica et Cosmochimica Acta,1989,53(10):2591-2608.

    [67] Crowe B R,F(xiàn)innegan D L,Zoller W H,et al.Trace element geochemistry of volcanic gases and particles from 1983—1984 eruptive episodes of Kilauea Volcano[J].Journal of Geophysical Research:Solid Earth,1987,92(B13):13708-13714.

    [68] Symonds R B,Rose W,Reed M H,et al.Volatilization,transport and sublimation of metallic and nonmetallic elements in high temperature gases at Merapi Volcano,Indonesia[J].Geochimica et Cosmochimica Acta,1987,51(8):2083-2101.

    [69] Williams-Jones A E,Heinrich C A.Vapor transport of metals and the formation of magmatic-hydrothermal ore deposits[J]. Economic Geology,2005,100(7):1287-1312.

    [70] Heinrich C A,Driesner T,Stefansson A,et al.Magmatic vapor contraction and the transport of gold from the porphyry environment to epithermal ore deposits[J].Geology,2004,32(9):761-764.

    [71] Ulrich T,Guenther D,Heinrich C A.Gold concentrations of magmatic brines and the metal budget of porphyry copper deposits[J].Nature,1999,399(6737):676-679.

    [72] Jana D,Walker D.Core formation in the presence of various C—H—O volatile species[J].Geochimica et Cosmochimica Acta,1999,63(15):2299-2310.

    [73] Lowenstern J B,Mahood G A,Rivers M L.Evidence for extreme partitioning of copper into a magmatic vapor phase[J]. Science,1991,252(5011):1405-1409.

    [74] Lowenstern J B.Carbon dioxide in magmas and implications for hydrothermal systems[J].Mineralium Deposita,2001,36(6):490-502.

    [75] Yang Kaihui,Scott S D.Possible contribution of a metal-rich magmatic fluid to a sea-floor hydrothermal system[J].Nature,1996,383(6599):420-423.

    [76] Ballhaus C,Sylvester P.Noble metal enrichment processes in the merensky reef,Bushveld complex[J].Journal of Petrology,2000,41(4):545-561.

    [77] Wirth R,Reid D,Schreiber A.Nanometer-sized platinumgroup minerals(PGM)in base metal sulfides:New evidence for an orthomagmatic origin of the Merensky Reef PGE ore deposit,Bushveld Complex,South Africa[J].The Canadian Mineralogist,2013,51(1):143-155.

    [78] Malitch K N,Auge T,Badanina I Yu,et al.Os-rich nuggets from Au-PGE placers of the Maimecha-Kotui Province,Russia:A multi-disciplinary study[J].Mineralogy and Petrology,2002,76(1/2):121-148.

    [79] Burgisser A,Scaillet B,Harshvarhan H.Chemical patterns of erupting silicic magmas and their influence on the amount of degassing during ascent[J].Journal of Geophysical Research:Solid Earth,2008,113(B12):204-217.

    [80] Simon A C,Pettke T,Candela P A,et al.Copper partitioning in a melt-vapor-brine-magnetite-pyrrhotite assemblage[J].Geochimica et Cosmochimica Acta,2006,70(22):5583-5600.

    [81] Fu Piaoer,Tang Qingyan,Zhang Mingjie,et al.The ore genesis of Kalatongke Cu-Ni Sulfide Deposit,west China[J].Acta Geologica Sinica,2012,86(3):568-578.

    [82] Tang Qingyan,Zhang Mingjie,Li C,et al.The chemical compositions and abundances of volatiles in the Siberian large igneous province:Constraints on magmatic CO2and SO2emissions into the atmosphere[J].Chemical Geology,2013,339:84-91.

    [83] 傅飄兒.新疆北部晚古生代巖漿銅鎳硫化物礦床成因:巖石及流體地球化學(xué)研究[D].蘭州:蘭州大學(xué),2012.

    [84] 湯慶艷.峨眉山二疊紀(jì)地幔柱巖漿銅鎳鉑族硫化物礦床成礦體系對(duì)比[D].蘭州:蘭州大學(xué),2013.

    [85] Zhang Mingjie,Tang Qingyan,Hu Peiqing,et al.Noble gas isotopic constraints on the origin and evolution of the Jinchuan Ni-Cu-(PGE)sulfide ore-bearing ultramafic intrusion,Western China[J].Chemical Geology,2013,339:301-312.

    [86] Konnikov E G,Meurer W P,Neruchev S S,et al.Fluid regime of platinum group elements(PGE)and gold-bearing reef formationintheDovyrenmafic-ultramaficlayeredcomplex,easternSiberia,Russia[J].MineraliumDeposita,2000,35(6):526-532.

    The roles of volatiles in mineralizations of magmatic Ni-Cu-PGE sulfide deposits—Implications for potential metallogenic mechanism of super-large scale magmatic deposits in small magma

    Zhang Mingjie1,Tang Qingyan1,Li Wenyuan1,2,Yu Ming1,Hu Peiqing1,Li Jianping1
    (1.School of Earth Sciences,Lanzhou University,Lanzhou 730000,China;2.Xi’an Institute of Geology and Mineral Resources,Xi’an 710054,China)

    The resources of nickel and platinum group element(PGE)are contributed by super-large scale magmatic Ni-PGE sulfide deposits in the world.In this paper,on the basis of the behaviors of Ni,Cu and PGE in the processes of partial melting,magma crystallization and sulfide segregation etc.mantle magmatic events,we summarized the types and scales of oreforming magma,the gathering way of ore-forming metal elements,the controlling factors of sulfur saturation and the immiscibility of the sulfide liquid phase during the magmatism of Ni-Cu-PGE sulfide deposits,and discussed the role of volatiles in the transport and enrichment of ore-forming metals.It is concluded that the high degree of partial melting(i.e.,magma with high Mg content)in mantle is crucial to gather the huge amounts of Ni and PGE from mantle sulfide and olivine into ore-forming magma,and the sulfide segregation will form super-large scale Ni-PGE sulfide deposits in a large scale layered(or conduit)intrusion.It is inferred that the supercritical fluid(volatile)in mantle can carry and aggregate a large number of Ni-PGE metals from a large volume of mantle,and make a great contribution to the formation of Ni-PGE sulfide deposits.It is possible to be the metallogenic mechanism of magmatic Ni-PGE sulfide deposits by small mantle derived magma(i.e.,small intrusion).

    volatile;small intrusion;mantle-derived magma;super-large scale magmatic deposits;mineralizations

    O35;P611

    A

    1009-1742(2015)02-0040-10

    2014-11-20

    國家自然科學(xué)基金項(xiàng)目(41372095,41472070,41072056);高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金項(xiàng)目(20120211110023);中國地質(zhì)調(diào)查局項(xiàng)目(12120114044401)

    張銘杰,1965年出生,男,甘肅漳縣人,教授,博士生導(dǎo)師,主要從事巖石圈演化與成礦研究工作;E-mail:mjzhang@lzu.edu.cn

    猜你喜歡
    硫化物巖漿礦床
    構(gòu)造疊加暈法在深部找礦中的應(yīng)用——以河南小秦嶺楊砦峪金礦床S60號(hào)礦脈為例
    黑龍江省林口縣三合村探明超大型石墨礦床
    巖漿里可以開采出礦物質(zhì)嗎?
    火山冬天——巖漿帶來的寒冷
    西昆侖新發(fā)現(xiàn)鹽湖型鹵水硼鋰礦床
    大洋多金屬硫化物自然氧化行為研究
    連續(xù)流動(dòng)法測(cè)定沉積物中的酸揮發(fā)性硫化物
    Li2S-P2S5及Li2S-SiS2基硫化物固體電解質(zhì)研究進(jìn)展
    鎢酸錳催化氧化脫除模擬油硫化物
    遼南分水金礦床鉛同位素特征及礦床成因
    av在线天堂中文字幕| 国产成人免费观看mmmm| 性插视频无遮挡在线免费观看| 色噜噜av男人的天堂激情| 久久精品人妻少妇| 国产免费又黄又爽又色| 狂野欧美白嫩少妇大欣赏| 噜噜噜噜噜久久久久久91| 最近中文字幕2019免费版| 岛国毛片在线播放| 久久精品国产亚洲av天美| 亚洲国产成人一精品久久久| 一个人观看的视频www高清免费观看| 色综合站精品国产| 亚洲欧美日韩无卡精品| 高清午夜精品一区二区三区| 国产国拍精品亚洲av在线观看| 麻豆成人av视频| 国产精品麻豆人妻色哟哟久久 | 91狼人影院| 成人亚洲精品av一区二区| 国产精品美女特级片免费视频播放器| 一区二区三区乱码不卡18| 中文字幕熟女人妻在线| 91aial.com中文字幕在线观看| 精品不卡国产一区二区三区| 天天一区二区日本电影三级| 可以在线观看毛片的网站| 亚洲高清免费不卡视频| 精品免费久久久久久久清纯| 婷婷色麻豆天堂久久 | 黑人高潮一二区| 22中文网久久字幕| 精品酒店卫生间| 中文字幕制服av| 禁无遮挡网站| 国产精品久久久久久精品电影| 亚洲成人av在线免费| 春色校园在线视频观看| 日韩大片免费观看网站 | 99热这里只有是精品50| 日日摸夜夜添夜夜爱| 日韩三级伦理在线观看| 日本五十路高清| 亚洲欧美精品自产自拍| 欧美日韩综合久久久久久| 麻豆成人av视频| 中文资源天堂在线| 国产精品电影一区二区三区| 国产精品不卡视频一区二区| 精品久久久久久成人av| 午夜激情欧美在线| 一级毛片我不卡| 青春草亚洲视频在线观看| 久久久成人免费电影| 蜜臀久久99精品久久宅男| 国产老妇女一区| 亚洲欧美日韩东京热| 可以在线观看毛片的网站| 中文字幕久久专区| 男女那种视频在线观看| 91精品国产九色| 国产成人aa在线观看| 亚洲av免费在线观看| 精品少妇黑人巨大在线播放 | 久久草成人影院| 久久久久久九九精品二区国产| 国产成人精品一,二区| 欧美3d第一页| 99久国产av精品| 国产单亲对白刺激| 国产亚洲一区二区精品| 如何舔出高潮| 亚洲欧美一区二区三区国产| 嫩草影院入口| 两个人视频免费观看高清| 日韩强制内射视频| 国产精品国产三级专区第一集| 欧美激情久久久久久爽电影| 国产乱来视频区| 国产亚洲精品久久久com| 好男人视频免费观看在线| 麻豆国产97在线/欧美| 亚洲精品乱久久久久久| 亚洲av免费在线观看| 亚洲av免费在线观看| 亚洲av免费在线观看| 亚洲综合色惰| 精品午夜福利在线看| 欧美一区二区国产精品久久精品| 国产精品人妻久久久影院| 亚洲精品乱久久久久久| 美女黄网站色视频| 人妻制服诱惑在线中文字幕| 网址你懂的国产日韩在线| 久久精品国产亚洲av涩爱| 欧美日韩国产亚洲二区| 免费观看在线日韩| 国内少妇人妻偷人精品xxx网站| 国产高清不卡午夜福利| 69av精品久久久久久| 少妇裸体淫交视频免费看高清| 综合色丁香网| 观看美女的网站| 国产精品人妻久久久影院| 亚洲精品色激情综合| 亚洲精品色激情综合| 久久婷婷人人爽人人干人人爱| 男女下面进入的视频免费午夜| 久久久亚洲精品成人影院| 日本-黄色视频高清免费观看| a级毛色黄片| 成人午夜高清在线视频| 色综合站精品国产| 日韩强制内射视频| 国产精品女同一区二区软件| 最近中文字幕高清免费大全6| 久久6这里有精品| 亚洲精品国产av成人精品| 搡老妇女老女人老熟妇| 免费电影在线观看免费观看| 国产精品久久久久久精品电影| 26uuu在线亚洲综合色| 国产三级中文精品| www.色视频.com| 波多野结衣高清无吗| 波多野结衣高清无吗| 免费电影在线观看免费观看| 日韩欧美三级三区| 国产成人免费观看mmmm| 最近最新中文字幕免费大全7| 成人国产麻豆网| 日韩中字成人| 97人妻精品一区二区三区麻豆| 自拍偷自拍亚洲精品老妇| 国产 一区精品| 国产真实伦视频高清在线观看| 精品不卡国产一区二区三区| 国产激情偷乱视频一区二区| 国产午夜精品一二区理论片| 久久久a久久爽久久v久久| 亚洲人成网站在线观看播放| 国产毛片a区久久久久| 国产白丝娇喘喷水9色精品| 国产综合懂色| 亚洲在线自拍视频| 赤兔流量卡办理| 国产亚洲精品av在线| 不卡视频在线观看欧美| 精品国产露脸久久av麻豆 | 亚洲av一区综合| 色播亚洲综合网| 在现免费观看毛片| 亚洲国产精品成人久久小说| 中文字幕精品亚洲无线码一区| 中文字幕亚洲精品专区| 国产在线一区二区三区精 | 日本黄色视频三级网站网址| 国产真实伦视频高清在线观看| 美女黄网站色视频| 精品久久久久久久久av| 天堂√8在线中文| 淫秽高清视频在线观看| 国产精品日韩av在线免费观看| 激情 狠狠 欧美| 精品久久久久久久末码| 久久久久网色| 韩国av在线不卡| 别揉我奶头 嗯啊视频| 国产又黄又爽又无遮挡在线| 又黄又爽又刺激的免费视频.| 联通29元200g的流量卡| 国产精品电影一区二区三区| 中文字幕av在线有码专区| 国产精品久久久久久久电影| 国产 一区精品| 男女边吃奶边做爰视频| 国产黄色小视频在线观看| 久久韩国三级中文字幕| 能在线免费观看的黄片| 蜜桃亚洲精品一区二区三区| 又爽又黄a免费视频| 国产探花在线观看一区二区| 国产免费男女视频| 在线天堂最新版资源| 国产免费一级a男人的天堂| 日韩av不卡免费在线播放| 97人妻精品一区二区三区麻豆| 亚洲精品日韩av片在线观看| 边亲边吃奶的免费视频| 汤姆久久久久久久影院中文字幕 | 秋霞伦理黄片| 国产伦一二天堂av在线观看| 伊人久久精品亚洲午夜| 最近中文字幕2019免费版| 久久热精品热| 26uuu在线亚洲综合色| 尾随美女入室| 丝袜美腿在线中文| 深爱激情五月婷婷| 亚洲国产精品成人综合色| 日韩强制内射视频| 免费搜索国产男女视频| 少妇熟女欧美另类| 国产大屁股一区二区在线视频| 91狼人影院| 两个人的视频大全免费| 3wmmmm亚洲av在线观看| 久久久久久国产a免费观看| 亚洲高清免费不卡视频| 国产精品国产三级国产专区5o | 大话2 男鬼变身卡| 最新中文字幕久久久久| 综合色av麻豆| 亚洲五月天丁香| 久久精品国产亚洲网站| 亚洲三级黄色毛片| 伦理电影大哥的女人| 亚洲怡红院男人天堂| 国产一区二区亚洲精品在线观看| 97热精品久久久久久| 亚洲人成网站在线播| 精品国内亚洲2022精品成人| 国产精品乱码一区二三区的特点| 国产高潮美女av| 国产亚洲91精品色在线| 久久国内精品自在自线图片| 亚洲欧美日韩高清专用| 99热精品在线国产| 色综合色国产| 国产精品久久久久久久电影| 神马国产精品三级电影在线观看| 久久精品夜夜夜夜夜久久蜜豆| 床上黄色一级片| 亚洲av熟女| 99热全是精品| 天堂影院成人在线观看| 国产伦一二天堂av在线观看| 久久久久国产网址| 久久久午夜欧美精品| 久久婷婷人人爽人人干人人爱| 91在线精品国自产拍蜜月| 丰满少妇做爰视频| 99久久精品国产国产毛片| 午夜福利成人在线免费观看| 大又大粗又爽又黄少妇毛片口| 伦精品一区二区三区| 永久网站在线| 精品一区二区免费观看| 69av精品久久久久久| 直男gayav资源| 特级一级黄色大片| 在线a可以看的网站| 亚洲av中文av极速乱| 麻豆国产97在线/欧美| 成人av在线播放网站| 秋霞伦理黄片| 亚洲激情五月婷婷啪啪| 97在线视频观看| 日韩亚洲欧美综合| 国产色婷婷99| 久久精品91蜜桃| 日韩制服骚丝袜av| 中国国产av一级| 日韩欧美精品免费久久| 久久人人爽人人片av| 欧美另类亚洲清纯唯美| 久久久久久久久中文| 日本黄色视频三级网站网址| 美女国产视频在线观看| 日本欧美国产在线视频| 麻豆成人av视频| 婷婷色综合大香蕉| 亚洲av一区综合| 国产精品国产三级国产专区5o | 欧美最新免费一区二区三区| 1000部很黄的大片| 久久99热这里只有精品18| 国产黄色小视频在线观看| 国产亚洲av嫩草精品影院| 久久精品国产鲁丝片午夜精品| 两性午夜刺激爽爽歪歪视频在线观看| 国内少妇人妻偷人精品xxx网站| 最近最新中文字幕免费大全7| 亚洲va在线va天堂va国产| videos熟女内射| 看免费成人av毛片| 久久久亚洲精品成人影院| 精品久久久久久电影网 | 99久久成人亚洲精品观看| 日本熟妇午夜| 精华霜和精华液先用哪个| 少妇裸体淫交视频免费看高清| 欧美xxxx性猛交bbbb| kizo精华| 精品人妻视频免费看| 一个人免费在线观看电影| 日本黄色片子视频| 欧美潮喷喷水| 亚洲激情五月婷婷啪啪| 黑人高潮一二区| 天美传媒精品一区二区| 欧美zozozo另类| 亚洲av免费高清在线观看| 国产精品蜜桃在线观看| 亚洲国产日韩欧美精品在线观看| 免费观看在线日韩| 国产精品av视频在线免费观看| 日韩三级伦理在线观看| 国语对白做爰xxxⅹ性视频网站| 国产又黄又爽又无遮挡在线| 最近中文字幕2019免费版| 亚洲av男天堂| 日韩高清综合在线| 少妇丰满av| 成人午夜高清在线视频| 99久久人妻综合| 内射极品少妇av片p| 久久久亚洲精品成人影院| 97超视频在线观看视频| 日本免费a在线| 亚洲四区av| 九九爱精品视频在线观看| 久久这里有精品视频免费| 在线观看66精品国产| 亚洲精品,欧美精品| 亚洲av电影在线观看一区二区三区 | 九九在线视频观看精品| 在线观看av片永久免费下载| 九九久久精品国产亚洲av麻豆| 久久精品久久精品一区二区三区| 午夜福利视频1000在线观看| 蜜桃亚洲精品一区二区三区| 亚洲av成人精品一区久久| 国产极品精品免费视频能看的| 国产一区二区在线观看日韩| 成人特级av手机在线观看| 亚洲va在线va天堂va国产| 国产亚洲av片在线观看秒播厂 | 日韩一区二区三区影片| 国产在线男女| 日本与韩国留学比较| 卡戴珊不雅视频在线播放| 麻豆av噜噜一区二区三区| 成人亚洲欧美一区二区av| 最近最新中文字幕免费大全7| 国产精品精品国产色婷婷| 午夜免费激情av| 国产在视频线在精品| 日本免费在线观看一区| 长腿黑丝高跟| 非洲黑人性xxxx精品又粗又长| 天天一区二区日本电影三级| 久久这里只有精品中国| 一本久久精品| 中文字幕久久专区| 日韩一本色道免费dvd| 国产精品爽爽va在线观看网站| 免费在线观看成人毛片| 欧美日本视频| 久久欧美精品欧美久久欧美| 免费观看人在逋| 中文字幕人妻熟人妻熟丝袜美| 看免费成人av毛片| 久久久久久久久大av| 国产淫片久久久久久久久| 久久亚洲国产成人精品v| 日韩精品有码人妻一区| 国产黄片视频在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 黄色日韩在线| 亚洲国产精品专区欧美| 精品一区二区三区视频在线| 91精品一卡2卡3卡4卡| 亚洲av成人精品一二三区| 日韩亚洲欧美综合| 少妇熟女aⅴ在线视频| 久久这里只有精品中国| 人妻少妇偷人精品九色| 看十八女毛片水多多多| 久久久久久久久久黄片| 亚洲婷婷狠狠爱综合网| 午夜福利成人在线免费观看| 午夜福利网站1000一区二区三区| 国产精品日韩av在线免费观看| 婷婷色av中文字幕| 熟女电影av网| 国产精品三级大全| 欧美成人免费av一区二区三区| 日本黄色视频三级网站网址| 最近视频中文字幕2019在线8| 晚上一个人看的免费电影| 插逼视频在线观看| 色综合站精品国产| 麻豆精品久久久久久蜜桃| 亚洲欧美日韩东京热| 日韩中字成人| 好男人视频免费观看在线| 国产在线男女| 日产精品乱码卡一卡2卡三| 听说在线观看完整版免费高清| 色尼玛亚洲综合影院| 午夜福利高清视频| 久久精品久久精品一区二区三区| 久久精品国产亚洲av天美| 插逼视频在线观看| 日韩制服骚丝袜av| 亚洲av免费高清在线观看| 国产精品人妻久久久久久| av国产免费在线观看| 美女内射精品一级片tv| 国产探花在线观看一区二区| 国产av在哪里看| 国产精品精品国产色婷婷| 国语对白做爰xxxⅹ性视频网站| 99久国产av精品| 麻豆成人av视频| 中文亚洲av片在线观看爽| 变态另类丝袜制服| 又爽又黄无遮挡网站| 欧美最新免费一区二区三区| 国产亚洲91精品色在线| 国产亚洲一区二区精品| 欧美区成人在线视频| 最近视频中文字幕2019在线8| 久久人人爽人人片av| 亚洲在久久综合| 七月丁香在线播放| 国产美女午夜福利| 国产精品乱码一区二三区的特点| 国产精品久久久久久精品电影| 亚洲va在线va天堂va国产| a级一级毛片免费在线观看| 深夜a级毛片| 一边摸一边抽搐一进一小说| 青春草视频在线免费观看| 国产亚洲5aaaaa淫片| 六月丁香七月| 亚洲三级黄色毛片| 亚洲真实伦在线观看| av卡一久久| 成人欧美大片| 亚洲综合色惰| 久久久国产成人精品二区| 日本黄大片高清| 观看美女的网站| 国产乱人偷精品视频| 欧美一区二区精品小视频在线| 国产高清三级在线| 国产av码专区亚洲av| 亚洲av成人av| 国产单亲对白刺激| 日韩欧美 国产精品| 国产黄色小视频在线观看| 日本黄色片子视频| 亚洲电影在线观看av| 久久久久性生活片| 99热全是精品| 国产高清不卡午夜福利| 九色成人免费人妻av| 一卡2卡三卡四卡精品乱码亚洲| 天堂中文最新版在线下载 | 亚洲丝袜综合中文字幕| av国产免费在线观看| 国产精品一二三区在线看| 免费观看在线日韩| 少妇猛男粗大的猛烈进出视频 | 日韩三级伦理在线观看| 十八禁国产超污无遮挡网站| 日日撸夜夜添| 亚洲精品乱码久久久v下载方式| 深夜a级毛片| 亚洲一区高清亚洲精品| 亚洲精品日韩在线中文字幕| 久久久久网色| ponron亚洲| 麻豆成人午夜福利视频| 国产极品精品免费视频能看的| 少妇猛男粗大的猛烈进出视频 | 国产精品嫩草影院av在线观看| 色综合色国产| 国产视频首页在线观看| 高清在线视频一区二区三区 | a级毛片免费高清观看在线播放| 高清毛片免费看| 亚洲精品亚洲一区二区| 在现免费观看毛片| 精品酒店卫生间| 亚洲精品国产成人久久av| 久久久精品大字幕| 国内揄拍国产精品人妻在线| 熟女人妻精品中文字幕| 99视频精品全部免费 在线| 一夜夜www| 国产日韩欧美在线精品| 亚洲精品日韩av片在线观看| 日韩 亚洲 欧美在线| 九九在线视频观看精品| 搡老妇女老女人老熟妇| 亚洲无线观看免费| 少妇的逼好多水| 天堂影院成人在线观看| 亚洲自拍偷在线| 久久人人爽人人片av| 亚洲欧美精品自产自拍| 卡戴珊不雅视频在线播放| 日韩在线高清观看一区二区三区| 亚洲国产欧美在线一区| 中文字幕免费在线视频6| 欧美日韩在线观看h| 少妇人妻精品综合一区二区| 亚洲电影在线观看av| 亚洲国产精品成人综合色| 91狼人影院| 欧美一区二区国产精品久久精品| 日本一本二区三区精品| 婷婷六月久久综合丁香| 国产v大片淫在线免费观看| 亚洲av中文av极速乱| 欧美色视频一区免费| 在线a可以看的网站| 狠狠狠狠99中文字幕| 中文字幕av在线有码专区| 亚洲国产欧洲综合997久久,| 久久久久久久久久久丰满| 精品人妻视频免费看| 久久精品久久久久久久性| 日韩一本色道免费dvd| 免费观看的影片在线观看| 亚洲欧美日韩东京热| 色综合站精品国产| 男女那种视频在线观看| 在线观看一区二区三区| 国产69精品久久久久777片| 亚洲精品自拍成人| 国产亚洲精品久久久com| 最近中文字幕2019免费版| 亚洲欧美日韩卡通动漫| 在线播放无遮挡| 久久人人爽人人爽人人片va| 国产色爽女视频免费观看| www.色视频.com| 亚洲美女搞黄在线观看| 黄片无遮挡物在线观看| 成人鲁丝片一二三区免费| 日韩av不卡免费在线播放| 好男人在线观看高清免费视频| 精品国产一区二区三区久久久樱花 | 国产高清三级在线| 国产伦精品一区二区三区视频9| 蜜桃亚洲精品一区二区三区| 亚洲真实伦在线观看| 51国产日韩欧美| ponron亚洲| 99国产精品一区二区蜜桃av| 伦精品一区二区三区| 18禁在线播放成人免费| 亚洲精华国产精华液的使用体验| 亚洲欧洲日产国产| 一级毛片aaaaaa免费看小| 亚洲18禁久久av| 在线播放国产精品三级| 亚洲怡红院男人天堂| 级片在线观看| 国产成人91sexporn| 高清毛片免费看| 草草在线视频免费看| 亚洲精品一区蜜桃| 亚洲怡红院男人天堂| 高清毛片免费看| 草草在线视频免费看| av黄色大香蕉| 久久久成人免费电影| 久久精品国产鲁丝片午夜精品| 精品国产三级普通话版| 亚洲国产成人一精品久久久| 日韩中字成人| 特大巨黑吊av在线直播| 3wmmmm亚洲av在线观看| 亚洲国产成人一精品久久久| 国产不卡一卡二| 中文字幕制服av| 久久精品夜夜夜夜夜久久蜜豆| 国产精品久久电影中文字幕| 在线观看一区二区三区| 中文字幕免费在线视频6| 成人三级黄色视频| 99久久中文字幕三级久久日本| 高清毛片免费看| 美女黄网站色视频| 夜夜爽夜夜爽视频| a级毛片免费高清观看在线播放| 久久久久久久久久成人| 六月丁香七月| 天堂中文最新版在线下载 | 亚洲丝袜综合中文字幕| 麻豆国产97在线/欧美| 欧美高清成人免费视频www| 91久久精品国产一区二区成人| 成人亚洲欧美一区二区av| 男人的好看免费观看在线视频| 亚洲欧美精品自产自拍| 午夜精品在线福利| 午夜日本视频在线| 又爽又黄无遮挡网站| 成人国产麻豆网| 啦啦啦观看免费观看视频高清| 久久精品综合一区二区三区| 日韩欧美在线乱码| 两个人视频免费观看高清| 在线天堂最新版资源| 亚洲人与动物交配视频| 深夜a级毛片| 最近最新中文字幕大全电影3| 国产黄a三级三级三级人|