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

    贛南地區(qū)石英脈型鎢礦成礦流體特征*

    2016-08-02 08:37:27周龍全李光來
    華東地質(zhì) 2016年2期

    周龍全,李光來,唐 傲,蘇 曄

    (東華理工大學核資源與環(huán)境省部共建國家重點實驗室培育基地,南昌 330013)

    ?

    贛南地區(qū)石英脈型鎢礦成礦流體特征*

    周龍全,李光來,唐傲,蘇曄

    (東華理工大學核資源與環(huán)境省部共建國家重點實驗室培育基地,南昌 330013)

    摘要:贛南是我國鎢礦床最密集的地區(qū),尤以石英脈型鎢礦最為發(fā)育。本文通過分析近年來該區(qū)石英脈型鎢礦流體包裹體類型、流體包裹體特征、顯微測溫、激光拉曼光譜等方面的最新成果,結(jié)合碳、氫、氧及鍶同位素的研究成果,探討贛南石英脈型鎢礦的流體特征,重點探討石英脈型鎢礦形成過程中的流體演化。認為贛南石英脈型鎢礦成礦流體主要來源于巖漿水,流體演化始于高溫高鹽度的巖漿—熱液過渡階段,與黑鎢礦沉淀密切相關(guān)的流體溫度主要集中于260~360 ℃,鹽度主要集中于4~9 wt% NaCl eq.,屬中—低鹽度、富含SiO2、揮發(fā)組分及多種成礦元素的熱液體系;礦質(zhì)主要以流體沸騰和混合作用為主,自然冷卻僅為少數(shù)礦床的主要礦石沉淀機制。

    關(guān)鍵詞:石英脈型鎢礦;成礦流體;贛南地區(qū)

    南嶺是中國乃至世界上鎢錫稀有金屬的重要產(chǎn)地,空間上具有“東鎢西錫”的特點,即南嶺東段以鎢礦為主,中段鎢錫礦并重,西段以錫礦為主[1]。其中贛南地區(qū)以石英脈型鎢礦的密集產(chǎn)出為特征,且還發(fā)育矽卡巖型、云英巖型、偉晶巖型、破碎帶型等多種類型的鎢礦床,是世界著名的鎢礦集區(qū),素有“世界鎢都”的美譽,該區(qū)發(fā)育的石英脈型鎢礦在中國乃至世界范圍內(nèi)具有代表性。

    對熱液礦床而言,成礦流體是成礦元素活化、遷移乃至沉淀的重要媒介[2-3],長期以來其性質(zhì)、來源及演化歷史一直是礦床學研究的重要內(nèi)容[2-9]。隨著現(xiàn)代測試技術(shù)的不斷發(fā)展,一些學者已不滿足于利用透明礦物的顯微測溫、爆裂法測溫、C-H-O穩(wěn)定同位素、顯微激光拉曼光譜分析等常規(guī)手段研究成礦流體,不透明礦物流體包裹體的紅外顯微鏡研究成了研究成礦流體不可或缺的技術(shù)手段[10-13]。在新技術(shù)條件下,贛南石英脈型鎢礦成礦流體方面的研究取得了突破性進展,積累了大量研究資料[12-34]。本文收集近年來贛南石英脈型鎢床成礦流體的研究成果,總結(jié)該區(qū)石英脈型鎢礦成礦流體的物理化學特征,探討該類礦床成礦的流體演化特征。

    1地質(zhì)背景及礦床分布

    贛南位于歐亞板塊與濱西太平洋板塊消減帶內(nèi)側(cè)華夏板塊的羅霄褶皺帶中部,晚古生代諸廣—武夷隆起與粵北拗陷過渡帶,屬濱太平洋構(gòu)造域中生代構(gòu)造帶的南東部,次級構(gòu)造單元為南嶺緯向構(gòu)造帶東段與武夷山北東—北北東向構(gòu)造帶南段的復合部位[35]。該區(qū)構(gòu)造變形強烈,褶皺、斷裂發(fā)育,特別是中生代以來發(fā)生了多期次、多性質(zhì)、多方向的擠壓和拉張,形成了武夷、羅霄、萬洋三條規(guī)模巨大的北東—北北東向隆起帶和桂東—興國—石城、郴州—崇義—會昌、韶關(guān)—三南—尋烏三條東西向構(gòu)造—巖漿—成礦帶。研究區(qū)基底為震旦—志留紀碎屑巖和變碎屑巖,蓋層為從泥盆紀到三疊紀的碳酸鹽巖、泥灰?guī)r夾碎屑巖,而侏羅—白惡紀發(fā)育的地層以斷裂盆地中的碎屑巖、火山巖和紅層為代表[36]。晚古生代、寒武系、泥盆系地層是鎢的重要礦源層,多時代鎢的礦源層的存在是控制該區(qū)鎢礦床形成的主要因素之一[37]。

    該區(qū)花崗質(zhì)巖體出露眾多,以燕山期花崗巖最為發(fā)育,具有多期多階段成巖特點。石英脈型鎢礦的形成時代多集中于160~150 Ma[1,36,38-39]。與鎢礦關(guān)系密切的花崗巖類主要為改造型,少數(shù)為同熔型[40],而石英脈型鎢礦多與燕山早期隱伏的高演化花崗巖體有密切成因聯(lián)系,這類花崗巖的形成時代多集中于165~150 Ma[41-42],石英脈型鎢礦床主要分布于崇—猶—余礦集區(qū)、于山礦集區(qū)以及三南礦集區(qū),著名的石英脈型鎢礦有西華山、漂塘、茅坪、蕩坪、盤古山、黃沙、大吉山等鎢礦床。

    根據(jù)脈幅大小石英脈型鎢礦床分為石英大脈型和細脈帶型,根據(jù)礦體產(chǎn)出位置又可分為內(nèi)帶型、內(nèi)外帶型及外帶型三種類型。石英脈型鎢礦的礦體一般呈脈狀產(chǎn)于震旦系、寒武系、泥盆系地層和燕山期花崗巖體中,圍繞成礦巖體的內(nèi)外接觸帶成群成帶產(chǎn)出,礦脈常見波狀彎曲、膨大縮小、分支復合等現(xiàn)象,由于受不同構(gòu)造控制,礦脈常呈疊瓦狀、雁行、人字型等成群分布[38]。大部分礦脈有“上部擴散,下部收斂”的特征,剖面上形如一把打開的紙扇,“扇根”多插入花崗巖中,并逐漸尖滅。發(fā)育完全的礦床在垂向上具有“五層樓”礦化模式,如漂塘鎢礦、木梓園鎢礦、茅坪鎢礦、鋸板坑鎢礦,在一些石英脈型鎢礦床的“五層樓”下部常有“地下室”鎢礦床發(fā)育于隱伏巖體的頂部,典型的礦例有茅坪鎢礦、鐵山垅鎢礦、樟東坑鎢礦、淘錫坑鎢礦等。

    2流體包裹體研究

    2.1流體包裹體特征

    根據(jù)流體包裹體在室溫下相態(tài)分類準則及冷凍回溫過程中的相態(tài)變化等研究成果[2,6],贛南石英脈型鎢礦床的流體包裹體可分為五種類型:氣液兩相水溶液包裹體(Ⅰ型)、單相包裹體(Ⅱ型)、富CO2的三相包裹體(Ⅲ型)、含子礦物包裹體(Ⅳ型)、熔融包裹體(Ⅴ型)(表1)。比較而言,石英、黃玉中的流體包裹體類型較豐富,Ⅰ型和Ⅱ型包裹體普遍較發(fā)育,部分發(fā)育Ⅲ型和Ⅳ型包裹體,而綠柱石、螢石、錫石、黑鎢礦中的流體包裹體類型相對較為單一,以發(fā)育Ⅰ型包裹體為主,僅少數(shù)礦床流體包裹體中發(fā)現(xiàn)子礦物,且發(fā)育數(shù)量較少,主要為石鹽、方解石和含鎢酸根礦物[12,17,23,26]。研究表明,綠柱石除了大量發(fā)育Ⅰ型包裹體外,還發(fā)育Ⅴ型包裹體[14]。Ⅲ型包裹體在不同礦區(qū)的發(fā)育程度不同,主要發(fā)育于石英和黃玉中。各寄主礦物中的包裹體粒徑變化范圍為0.5~64 μm,絕大多數(shù)為5~25 μm,絕大多數(shù)氣相百分數(shù)為5%~45%,包裹體形狀多呈不規(guī)則形、圓形、橢圓形、長管狀,少數(shù)呈寄主礦物的負晶形,多呈孤立、成群、沿生長環(huán)帶生長、串珠狀等分布,往往多世代的包裹體共存于同一礦脈中。

    2.2流體包裹體均一溫度

    贛南部分石英脈型鎢礦床流體包裹體均一溫度如表1。對不同或同一礦床的礦化階段劃分總體上可歸納為硅酸鹽—氧化物石英脈階段、氧化物—硫化物—石英脈階段、碳酸鹽階段。 贛南地區(qū)石英脈型鎢礦不同礦物中包裹體均一溫度范圍為70~435 ℃[12-34],其中,硅酸鹽—氧化物階段的包裹體均一溫度主要集中于220~360 ℃之間;除新安子鎢礦氧化物—硫化物階段包裹體均一溫度處于較高溫度(260~300 ℃)外,大多數(shù)礦床主要為150~260 ℃;而碳酸鹽階段包裹體均一溫度相對較低,一般在230 ℃以下。在均一溫度直方圖中(圖1),與黑鎢礦共生的石英中的流體包裹體均一溫度變化較大,部分礦床集中在較明顯的兩個及以上的溫度區(qū)間。如木梓園鎢礦,均一溫度集中于180~240 ℃和260~300 ℃;盤古山鎢礦均一溫度集中于160~220 ℃和300~320 ℃;茅坪鎢礦均一溫度集中于190~250 ℃和300~360 ℃;新安子鎢礦均一溫度集中于180~230 ℃、260~300 ℃和320~360 ℃,其中,均一溫度的高溫區(qū)基本上處于硅酸鹽—氧化物階段,低溫間處于氧化物—硫化物階段。

    黑鎢礦中流體包裹體的均一溫度變化范圍為185~412 ℃,除盤古山均一溫度峰值集中于240~280 ℃外,大部分均一溫度峰值主要位于260~360 ℃,且黑鎢礦、錫石、黃玉中包裹體的均一溫度總體高于石英的均一溫度(圖1),在較高溫度(260~412 ℃)階段可見黑鎢礦、錫石、黃玉多在硅酸鹽—氧化物—石英脈階段大規(guī)模結(jié)晶,且黑鎢礦、錫石、黃玉中流體包裹體均一溫度相對石英變化較小,其中錫石、黃玉中流體包裹體的均一溫度處于相對狹窄范圍內(nèi)(290~360 ℃),說明錫石、黃玉相對黑鎢礦、石英而言結(jié)晶持續(xù)的時間較短,而不同礦物中流體包裹體均一溫度存在差異,這種差異主要可能由于所處的“結(jié)晶時段”不同,而非礦物形成后的影響[2,43]。石英是長期結(jié)晶的“貫穿性”礦物,其包裹體記錄了流體較長時間的演化歷史,同時與其它礦物相比,石英在應力作用和流體改造下更易遭受破壞,其流體包裹體的均一溫度顯示較大的變化范圍,因此其形成溫度并不能準確代表黑鎢礦形成溫度[10,28],相比之下,黑鎢礦、錫石、黃玉流體包裹體中的均一溫度處于相對狹窄且溫度較高的范圍內(nèi),較石英而言遭受改造程度小并能較多地保存原生包裹體[28],且在主成礦期間大量形成的礦石礦物(多為不透明礦物)記載了成礦流體的最佳成礦信息[10]。因此,該區(qū)黑鎢礦中的包裹體更真實地反應了礦床中鎢的形成條件,說明贛南地區(qū)石英脈型鎢礦床鎢成礦的主要溫度集中于260~360 ℃。

    表1 贛南石英脈鎢礦床流體包裹體特征

    圖1 贛南石英脈型鎢礦流體包裹體均一溫度直方圖Fig.1    Homogenization temperature histograms of fluid inclusions of quartz-vein type tungsten deposits in southern Jiangxi Province西華山數(shù)據(jù)[12];蕩坪、盤古山及大吉山數(shù)據(jù)[13];漂塘數(shù)據(jù)[13,27,28];茅坪數(shù)據(jù)[25];淘錫坑數(shù)據(jù)[26];新安子數(shù)據(jù)[29];大龍山數(shù)據(jù)[32];木梓園數(shù)據(jù)[21];黃沙數(shù)據(jù)[22]。

    2.3爆裂法測溫

    一般認為爆裂溫度代表成礦溫度的上限,是獲取爆裂溫度最快捷的方法。盤古山黑鎢礦的爆裂溫度為240~330 ℃,石英的爆裂溫度為290~350 ℃,較礦脈石英中的均一溫度185~325℃略高[44];大吉山鎢礦共測得238個爆裂溫度,溫度范圍為258~494 ℃,主要集中于350~400 ℃,此外還有28個數(shù)據(jù)>500 ℃,相比于石英中包裹體的均一溫度范圍(84~324 ℃)高很多,其爆裂曲線顯示多個爆裂縫[17-18]。爆裂溫度指示石英脈型鎢礦具有較高的溫度上限,多個爆裂縫的存在則可能反映了流體活動的多期多階段性特征。

    2.4流體包裹體鹽度、密度及壓力

    石英脈型鎢礦成礦流體的鹽度范圍為0.18~14.2 wt%NaCl eq.(表1),顯示其具有中—低鹽度的特點。其中石英中流體包裹體的鹽度峰值主要集中于4~8 wt%NaCl eq.,黃玉中流體包裹體的鹽度峰值主要集中于6~11 wt%NaCl eq.,綠柱石中流體包裹體的鹽度峰值主要集中于1~3 wt%NaCl eq.,黑鎢礦中流體包裹體的鹽度峰值主要集中于4~9 wt%NaCl eq.,錫石中流體包裹體的鹽度峰值主要集中于8~9 wt%NaCl eq.(圖2),說明黃玉、錫石、黑鎢礦中流體的鹽度略高于石英、綠柱石中流體的鹽度,且黑鎢礦成礦流體鹽度范圍較石英中流體的鹽度范圍窄。該區(qū)包裹體多為NaCl-H2O±CO2體系,利用溫度—密度相圖[45]、溫度—鹽度—密度相圖[46]及Flincor流體包裹體計算軟件[47]等不同方法確定的密度為0.55~1.03 g/cm3,成礦流體密度主要集中于0.7~0.99 g/cm3[15,23,20-21,26,28],獲得的成礦流體均一壓力介于6~156.6 MPa[14,15,20-21,23,26],除西華山鎢礦成礦深度為1.0~3.3 km外[12],成礦深度位于3.1~5.9 km[17,20,22-23]。

    3同位素地球化學

    在贛南石英脈型鎢礦研究中,碳氧同位素、氫氧同位素以及銣鍶同位素被廣泛應用于成礦流體來源的示蹤,并取得一系列成果[12,16,26,32-33,48-54]。

    圖2 贛南石英脈型鎢礦流體包裹體鹽度直方圖Fig.2 Salinity histograms for fluid inclusions of quartz-vein type tungsten deposits in Southern Jiangxi Province西華山數(shù)據(jù)[12];蕩坪、盤古山及大吉山數(shù)據(jù)[13];漂塘數(shù)據(jù)[13,23,30];茅坪數(shù)據(jù)[25];新安子數(shù)據(jù)[29];木梓園數(shù)據(jù)[21];黃沙數(shù)據(jù)[22]

    3.1碳氧同位素特征

    漂塘鎢礦礦脈中3個層解石樣品的δ13C值為-9.03‰~-5.86‰(平均為-7.06‰)[48],西華山鎢礦礦脈中6個層解石的δ13C值為-7.53‰~-4.51‰(平均為-6.50‰)[48];大吉山鎢礦礦脈中2個方解石的δ13C值為-8.2‰~-8.0‰[49]。漂塘、西華山鎢礦礦脈中層解石和大吉山鎢礦礦脈中方解石的碳同位素組成與巖漿碳的同位素組成(-8.0‰~-5.0‰)[35,49]基本一致,說明它們屬于巖漿—熱液成因。漂塘鎢礦礦脈中11個方解石樣品δ13C值為-6.31‰~4.4‰(其中3個δ13C值<-5‰)[48],δ13C值的變程較寬,且出現(xiàn)正值,較層解石更富13C,說明成礦流體可能因與圍巖發(fā)生水巖反應而提高了13C比例;也可能在成礦流體演化的晚期階段有大氣降水的兌入,將鈣質(zhì)圍巖淋濾而獲得13C,使稍晚階段的方解石δ13C值有所增高。

    3.2氫氧同位素特征

    氫氧同位素組成很大程度上可反應成礦流體中水的來源[55-56],因此被廣泛應用于與巖漿熱液有關(guān)的成礦流體的研究。贛南石英脈型鎢礦石英中氫氧同位素組成為:δD為-98‰~-45‰,δ18O為+7.3‰~+14.05‰。根據(jù)石英—水平衡分餾方程計算獲得了與石英呈平衡的礦液的δ18OH2O值變化范圍為-3.1‰~+10.8‰;黑鎢礦氫氧同位素組成為:δD為-135‰~-50.4‰,δ18O為+3.67‰~+7.61‰,根據(jù)黑鎢礦—水平衡分餾方程計算獲得了與黑鎢礦呈平衡的礦液的δ18OH2O值變化范圍為+4.90‰~+8.86‰(表2)。石英屬含氧礦物,幾乎不含氫原子數(shù),雖易與它所含的水發(fā)生同位素平衡再交換反應,但對流體包裹體的氫同位素組成所造成的影響很小[57],可以認為氫同位素組成代表了原始溶液的組成[26]。西華山、漂塘、茅坪、黃沙、盤古山、大吉山等礦床的石英及黑鎢礦流體中氫氧同位素數(shù)據(jù)大多數(shù)落在巖漿水區(qū)域,且礦液中的δD值基本位于-80‰~-50‰范圍(圖3),說明原始熱液源于巖漿水。少數(shù)氫氧同位素數(shù)據(jù)落于巖漿水附近區(qū)域,應是成礦晚階段的流體開始有大氣降水加入所致。而僅淘錫坑、大吉山鎢礦礦液中δ18OH2O值偏離正常巖漿水區(qū)域,顯示具有大氣降水的特征,原因可能為研究者較多地使用了成礦晚階段的樣品。

    表2 贛南石英脈型鎢礦氫氧同位素組成

    圖3 贛南石英脈型鎢礦成礦流體氫氧同位素特征圖Fig.3 Hydrogen and oxygen isotopic compositions of ore-forming fluids from quartz-vein type tungsten deposits in southern Jiangxi Province底圖[35];西華山、漂塘數(shù)據(jù)[48];大吉山數(shù)據(jù)[49];黃沙數(shù)據(jù)[51];茅坪數(shù)據(jù)[34];淘錫坑數(shù)據(jù)[26];盤古山數(shù)據(jù)[31]。

    3.3鍶同位素特征

    鍶同位素的初始比值(ISr)常被用來示蹤成礦物質(zhì)或成礦流體的來源。該區(qū)流體包裹體Rb-Sr法定年獲得的初始鍶比值(ISr)>0.714(大吉山除外),雖然淘錫坑鎢礦三組石英脈流體包裹體ISr值差異較大,卻也同樣>0.714(圖4),顯示成礦流體可能主要源于地殼。

    圖4 贛南成礦流體初始鍶比值Fig.4 Initial strontium ratios of ore-forming fluids in southern Jiangxi Province底圖、漂塘、柯樹嶺、牛嶺數(shù)據(jù)[58];西華山、蕩坪數(shù)據(jù)[16];淘錫坑數(shù)據(jù)[59]

    4討論

    4.1成礦流體來源

    根據(jù)贛南石英脈型鎢礦床同位素研究,除漂塘鎢礦礦脈中方解石δ13C值>-5‰外,西華山、漂塘、大吉山方解石、層解石δ13C值多數(shù)為-8‰~-5‰[35,48-49],與巖漿碳的同位素組成基本一致,顯示在成礦流體演化早階段碳主要來源于巖漿,而晚期階段由于大氣降水的兌入將鈣質(zhì)圍巖淋濾的13C帶入流體中,使稍晚階段的方解石δ13C值增高。石英的δ18O值變化范圍為+7.3‰~+14.1‰,黑鎢礦變化范圍為+4.90‰~+8.86‰,除淘錫坑外,所有石英的δ18O值變化為+10.58‰~+14.1‰(表2)。其中西華山、漂塘鎢礦全巖δ18O值變化為+8.23‰~+13.04‰[48,50,52],無論全巖還是石英的δ18O值均比“正?!睅r漿水(+6‰~+9‰)[55]高,而花崗巖較高的δ18O值與來源物質(zhì)的成因密切相關(guān)[52],由含高δ18O的侵入巖漿固結(jié)成巖后形成[50],且礦液中δD值基本為-80‰~-50‰,說明礦液主要來源于巖漿。另外,少數(shù)氫氧同位素數(shù)據(jù)落于巖漿水附近(圖3),一方面可能研究者較多使用了成礦晚階段的樣品所致,另一方可能是來自巖漿的原始流體與大氣降水不同程度混合導致。西華山鎢礦成礦流體在硅酸鹽—氧化物階段主要以巖漿水為主[15,48,50,52],魏文鳳等[15]認為西華山成礦流體δ18O值從硅酸鹽—氧化物階段到晚期碳酸鹽階段逐漸降低,并向雨水線方向漂移。漂塘、大吉山、淘錫坑鎢礦成礦流體在硅酸鹽—氧化物階段雖以巖漿水為主,但有相當數(shù)量的大氣降水參與,晚期大氣降水更顯著[15,26,49]。綜上所述,認為該區(qū)石英脈型鎢礦早期成礦流體主要來源于巖漿水,有不同程度大氣降水參與,到成礦晚階段大氣降水呈增多趨勢。

    4.2成礦流體性質(zhì)

    贛南石英脈型鎢礦成礦流體在演化中具多階段性和復雜性,不同成礦階段具有不同特征。同一礦床從早階段的硅酸鹽—氧化物階段至碳酸鹽階段,成礦溶液的鹽度大體呈降低趨勢,揮發(fā)組分明顯減少。如黃沙鎢礦在氧化物—石英階段石英中氣液兩相和含CO2三相包裹體均一溫度及鹽度范圍分別為256~336 ℃、4~9 wt%NaCl eq.和298~351℃、1.4~2.4wt%NaCl eq.,氣相組分檢測到H2O、CO2、CH4、N2等,硫化物—石英階段石英中氣液兩相包裹體均一溫度主要分布于190~240 ℃,鹽度為1~7 wt%NaCl eq.,氣相組分只檢測到H2O和少量CH4,顯示流體從高溫階段到低溫階段鹽度呈下降趨勢,氣相組分不斷減少[22];盤古山鎢礦主要礦化階段中早階段輝鉍礦—黑鎢礦—石英脈石英中的氣液兩相和含CO2三相包裹體的均一溫度分別集中于250~330℃和299~353 ℃,明顯高于晚階段(輝鉍礦)—黑鎢礦—石英脈石英中的氣液兩相和含CO2三相包裹體的均一溫度(150~220℃,213~299 ℃),從早階段到晚階段流體鹽度和密度相差不大,含CO2三相包裹體數(shù)量呈減少趨勢,而CO2單相包裹體數(shù)量則呈增多趨勢[20];淘錫坑鎢礦硅酸鹽-氧化物階段石英中流體包裹體均一溫度范圍集中在180~260 ℃和280~400 ℃兩個區(qū)間,鹽度集中于0~4.5 wt%NaCl eq.,較氧化物—硫化物階段溫度(160~260 ℃)高很多,鹽度(3~7 wt%NaCl eq.)相對增高,到碳酸鹽階段均一溫度(120~200 ℃)和鹽度(0~4.73 wt%NaCl eq.)較氧化物—硫化物階段呈降低趨勢,且流體中揮發(fā)組分不斷減少[26];西華山鎢礦黑鎢礦中流體包裹體均一溫度集中于320~370 ℃,鹽度集中于4~9 wt%NaCl eq.,硅酸鹽—氧化物階段的石英中包裹體的均一溫度、鹽度分別集中于240~320 ℃和1~6 wt%NaCl eq.,到硫化物階段則集中于180~230 ℃和1~5 wt%NaCl eq.,顯示黑鎢礦在較早階段形成,鹽度呈降低趨勢[15]。贛南石英脈型鎢礦成礦流體在硅酸鹽-氧化物階段成礦流體溫度主要集中于220~360 ℃[12-34],多數(shù)礦床富含H2O、CO2、CH4、N2等揮發(fā)組分,屬于高—中溫富含揮發(fā)分的熱液,且鎢大規(guī)模沉淀主要發(fā)生在該階段,溫度和鹽度主要集中于260~360 ℃和4~9 wt%NaCl eq.[13,15,25,26-29];氧化物—硫化物階段成礦流體溫度主要集中于150~260 ℃,鹽度較硅酸鹽—氧化物階段略低,揮發(fā)組分較硅酸鹽-氧化物階段明顯減少;碳酸鹽階段流體溫度多<230 ℃,鹽度相對較低。

    4.3成礦流體的流體演化

    研究表明南嶺地區(qū)與鎢錫礦化有關(guān)的巖漿巖及其成礦流體是巖漿—熱液過渡性流體[14,68-74]。西華山礦脈中綠柱石熔融包裹體、熔—流包裹體的發(fā)現(xiàn)及其與流體包裹體的共存[14,16,27,73],也證實了脈鎢礦床的成礦流體早期具有巖漿—熱液過渡性流體的特征,且富含SiO2、揮發(fā)組分及成礦元素[70],其演化始于巖漿—熱液過度階段,溫度可達650~820 ℃[16,27,73]。伴隨溫度、壓力降低及物化條件的改變,熔漿與熱液發(fā)生液態(tài)分離[74],經(jīng)液態(tài)分離作用分異的原始高溫高鹽度的巖漿熱液,在溫度差、壓力差和濃度差等驅(qū)動下運移。一方面,原始流體可與圍巖發(fā)生水巖反映,同時可能萃取含鎢圍巖中部分鎢元素;另一方面,運移的含礦溶液與流經(jīng)圍巖裂隙的大氣降水發(fā)生混合而促使鎢礦石沉淀。研究顯示贛南石英脈型鎢礦與黑鎢礦共生的石英流體包裹體均一溫度變化范圍為115~435 ℃,鹽度范圍為0.18~14.2 wt%NaCl eq.,且同種或不同種寄主礦物的包裹體類型、組成、均一溫度以及鹽度等既有相似也存在差異,說明成礦流體經(jīng)歷了復雜的演化過程。目前的研究成果記錄了72~435 ℃溫度階段的流體演化史[12-34],但卻缺少435~650 ℃溫度階段的數(shù)據(jù)(有待進一步研究)。從現(xiàn)有數(shù)據(jù)看,黑鎢礦大規(guī)模結(jié)晶的溫度主要發(fā)生在260~360 ℃階段,說明該區(qū)石英脈型鎢礦成礦流體經(jīng)液態(tài)分離后隨溫度、壓力降低以及物理化學條件的改變而經(jīng)歷了復雜的演化過程,在260~360 ℃溫度階段發(fā)生黑鎢礦沉淀。

    熱液型鎢礦礦石沉淀機制有流體不混溶作用[3,12, 63-64 ,75-76]、流體混合作用[3,11-12,62,66]和自然冷卻[62]等。由于礦床地質(zhì)環(huán)境、成礦流體物理化學條件等因素的異同,鎢沉淀的主導因素在不同礦床中具有差異。研究顯示大吉山、木梓園、西華山、盤古山、黃沙、茅坪、新安子等鎢礦床成礦流體經(jīng)歷了不同規(guī)模的沸騰作用[12-13,17,20-22,24-25,29],且多發(fā)生于早階段,溫度多介于250~340℃[17,20,21,22,25],與黑鎢礦主要成礦溫度(260~360℃)基本一致,說明流體的沸騰作用與鎢成礦關(guān)系密切。早期高溫階段流體灌入裂隙構(gòu)造,隨壓力驟然下降使流體產(chǎn)生減壓沸騰,同時大量CO2、CH4、H2S、N2等氣體揮發(fā)分從熱液中溢出,促使流體溶液中的酸性揮發(fā)性組分進入氣相,相分離的同時可使液相pH值升高,使鎢的絡合物的穩(wěn)定性降低,從而導致鎢礦物的結(jié)晶。伴隨成礦流體的演化,富含成礦元素的巖漿流體與通過礦區(qū)內(nèi)裂隙、斷裂的大氣降水發(fā)生混合,如西華山、淘錫坑、漂塘、黃沙、大吉山、茅坪等鎢礦床顯示了在成礦流體演化過程中發(fā)生了高溫、高鹽度的流體與低溫度、低鹽度流體不同程度的混合作用[12-13,22-26,34],引起流體體系物理化學條件不斷改變,使金屬絡合物分解并在有利的構(gòu)造地段沉淀成礦[12-13,22-23,26]。氫、氧同位素研究顯示也該區(qū)成礦流體主要源于巖漿水,到成礦晚階段逐漸有大氣降水參與。漂塘鎢礦床錫石中流體包裹體的均一溫度介于320~370℃,與Ni et al.[13]獲得的黑鎢礦的成礦溫度(320~360℃)基本吻合,王旭東等[23]認為與錫成礦相關(guān)的流體在成礦早期經(jīng)歷了自然冷卻過程,在成礦溫度上基本一致,說明與鎢成礦相關(guān)的成礦流體在高溫階段亦經(jīng)歷了自然冷卻的過程。自然冷卻的過程僅是該區(qū)少數(shù)鎢礦床主要的礦石沉淀機制,例如僅蕩坪、漂塘、大吉山鎢礦床的流體包裹體的鹽度變化不大且均一溫度有逐漸降低的趨勢[13,17,23],說明成礦流體經(jīng)歷自然冷卻的過程,導致成礦物質(zhì)溶解度隨溫度下降而降低,發(fā)生黑鎢礦沉淀。

    5結(jié)論

    贛南石英脈型鎢礦硅酸鹽—氧化物階段的流體與鎢成礦最密切的溫度集中于260~360 ℃,鹽度集中于4~9 wt%NaCl eq,壓力介于27~159 MPa,屬于高—中溫、中—低鹽度、富含揮發(fā)分的熱液體系。氧化物—硫化物階段屬于中—低溫、中—低鹽度熱液體系,主要為成礦晚期階段流體。碳、氫、氧同位素資料表明該區(qū)石英脈型鎢礦早期成礦流體主要來源于巖漿水,有不同程度大氣降水參與,到成礦晚階段大氣降水呈增多趨勢。贛南石英脈型鎢礦流體演化始于巖漿—熱液過渡階段,經(jīng)液態(tài)分離作用后隨溫度、壓力降低及物理化學條件的改變而經(jīng)歷復雜的演化過程,主要的礦石沉淀機制可能為流體沸騰和混合作用,自然冷卻是該區(qū)少數(shù)鎢礦床的主要礦石沉淀機制。

    參考文獻

    [1]華仁民,李光來,張文蘭,等.華南鎢和錫大規(guī)模成礦作用的差異及其原因初探[J].礦床地質(zhì),2010,29(1):9-23.

    [2]Roedder E. Fluid inclusions[J]. Mineralogical Society of America, Reviews in Mineralogy, 1984,12: 644.

    [3]Wilkinson J J. Fluid inclusions in hydrothermal ore deposits[J]. Lithos, 2001, 55(1/4): 229-272.

    [4]張文淮,陳紫英.流體包裹體地質(zhì)學[M].北京:中國地質(zhì)大學出版社,1993:1-242.

    [5]池國祥,賴健清.流體包裹體在礦床研究中的作用[J].礦床地質(zhì),2009,28(6):850-855.

    [6]盧煥章,范宏瑞,倪培,等.流體包裹體[M].北京:科學出版社,2004:406-419.

    [7]Naumov V B, Dorofeev V A, Mironova O F. Physicochemical parameters of the formation of hydrothermal deposits:a fluid inclusion study. I.Tin and Tungsten Deposits[J]. Geochemistry International, 2011, 49(10): 1002-1021.

    [8]Bodnar R J, Lecumberri-Sanchez P, Moncada D, et al. Fluid inclusions in hydrothermal ore deposits[C]//Treatise on Geochemistry (Second Edition), Geochemistry of Mineral Deposits, 2014: 119-142.

    [9]倪培,范宏瑞,丁俊英.流體包裹體研究進展[J].礦物巖石地球化學通報,2014,33(1):1-5.

    [10]Moritz R. Fluid salinities obtained by infrared microthemometry of opaque minerals Implications for ore deposit modeling-A note of caution[J]. Journal of Geochemical Exploration, 2006, 89: 284-287.

    [11]ZHU Yanan, PENG Jiantang. Infrared microthermometric and noble gas isotope study of fluid inclusions in ore minerals at the Woxi orogenic Au-Sb-W deposit, western Hunan, South China[J]. Ore Geology Reviews, 2015, 65(1): 55-69.

    [12]WEI Wenfeng, HU Ruizhong, BI Xianwu, et al. Infrared microthermometric and stable isotopic study of fluid inclusions in wolframite at the Xihuashan Tungsten deposit, Jiangxi province, China[J]. Mineralium Deposita, 2012, 47(6): 589-605.

    [13]NI Pei, WANG Xudong, WANG Guoguang, et al. An infrared microthermometric study of fluid inclusions in coexisting quartz and wolframite from Late Mesozoic Tungsten deposits in the Gannan metallogenic belt, South China[J]. Ore Geology Reviews, 2015, 65(4): 1062-1077.

    [14]常海亮,黃惠蘭.西華山鎢礦床中熔融包裹體的初步研究與礦床成因探討[J].巖石礦物學雜志,2002,21(2):143-150.

    [15]魏文鳳.贛南西華山鎢礦床成礦流體及成礦機制研究[D].北京:中國科學院研究生院,2011:1-122.

    [16]劉家齊,汪雄武,曾貽善,等.西華山花崗巖及鎢錫鈹?shù)V田成礦流體演化[J].華南地質(zhì)與礦產(chǎn),2002(3):91-96.

    [17]席斌斌,張德會,周利敏.江西省全南縣大吉山鎢礦成礦流體演化特征[J].地質(zhì)學報,2008,82(7):956-966.

    [18]徐文剛,張德會,席斌斌,等.流體包裹體爆裂法測溫技術(shù)可靠性討論——以江西大吉山鎢礦為例[J].現(xiàn)代地質(zhì),2008,22(5):757-765.

    [19]王旭東,倪培,蔣少涌,等.贛南漂塘鎢礦流體包裹體研究[J].巖石學報,2008,24(9):2163-2170.

    [20]王旭東,倪培,張伯聲,等.江西盤古山石英脈型鎢礦床流體包裹體研究[J].巖石礦物學雜志,2010,29(5):539-550.

    [21]王旭東,倪培,袁順達,等.贛南木梓園鎢礦流體包裹體特征及其地質(zhì)意義[J].中國地質(zhì),2012,39(6):1790-1797.

    [22]王旭東,倪培,袁順達,等.江西黃沙石英脈型鎢礦床流體包裹體研究[J].巖石學報,2012,28(1):122-132.

    [23]王旭東,倪培,袁順達,等.贛南漂塘鎢礦錫石及共生石英中流體包裹體研究[J].地質(zhì)學報,2013,87(6):850-859.

    [24]王旭東,倪培,袁順達,等.江西大吉山鎢多金屬礦床流體包裹體研究[J].礦床地質(zhì),2013,32(2):308-322.

    [25]胡東泉,華仁民,李光來,等.贛南茅坪鎢礦流體包裹體研究[J].高校地質(zhì)學報,2011,17(2):327-336.

    [26]宋生瓊.贛南淘錫坑鎢礦床流體地球化學和成礦機理研究[D].貴陽:中國科學院地球化學研究所,2011:1-115.

    [27]黃惠蘭,李芳,譚靖,等.贛南西華山黑鎢礦中熔融包裹體的發(fā)現(xiàn)和初步研究[J].華南地質(zhì)與礦產(chǎn),2012,28(2):181-183.

    [28]黃惠蘭,常海亮,李芳,等.西華山鎢礦床晶洞中水晶與黑鎢礦流體包裹體顯微測溫與特征元素測定[J].地學前緣,2013,20(2):205-212.

    [29]吳開興,張戀,陳陵康,等.贛南新安子鎢錫礦床流體包裹體地球化學研究[J].有色金屬科學與工程,2013,4(5):70-78.

    [30]趙波.江西漂塘石英脈型黑鎢礦床成礦深度估算[D].北京:中國地質(zhì)大學,2013:1-77.

    [31]李光來,華仁民,王旭東,等.江西南部盤古山鎢礦的氫氧同位素研究[J].東華理工大學學報:自然科學版,2014,37(2):164-169.

    [32]劉若蘭,慕紀錄,徐榮.江西大龍山鎢鉬礦床流體包裹體的研究[J].成都地質(zhì)學院學報,1985(2):25-31.

    [33]劉若蘭,慕紀錄.木梓園鎢鉬礦床流體包裹體研究及在成礦階段劃分中的應用[J].成都地質(zhì)學院學報,1990,17(3):18-28.

    [34]鄭躍鵬,喻鐵階,吳開華.茅坪鎢錫多金屬礦床流體包裹體特征及地質(zhì)意義[J].礦產(chǎn)與地質(zhì),1991,5(4):311-317.

    [35]朱焱齡,李崇佑,林運淮.江西南部脈鎢礦床主要礦化特征[C]//余鴻彰.鎢礦地質(zhì)討論會論文集.北京:地質(zhì)出版社,1981:339-351.

    [36]毛景文,謝桂青,郭春麗,等.南嶺地區(qū)大規(guī)模鎢錫多金屬成礦作用:成礦時限及地球動力學背景[J].巖石學報,2007,23(10):2329-2338.

    [37]徐克勤,胡受奚,孫明志,等.華南鎢礦床的區(qū)域成礦條件分析[C]//余鴻彰.鎢礦地質(zhì)討論會論文集.北京:地質(zhì)出版社,1981:243-258.

    [38]柳志青.脈狀鎢礦床成礦預測理論[M].北京:科學出版社,1980:1-7.

    [39]俞受鋆.南嶺地區(qū)燕山早期多階段花崗巖體與鎢礦成礦關(guān)系研究質(zhì)出版社[C]//余鴻彰.鎢礦地質(zhì)討論會論文集.北京:地質(zhì)出版社,1981:305-313.

    [40]徐克勤,程海.中國鎢礦形成的大地構(gòu)造背景[J].地質(zhì)找礦論叢,1987(3):1-7.

    [41]華仁民,陳培榮,張文蘭,等.華南中、新生代與花崗巖類有關(guān)的成礦系統(tǒng)[J].中國科學D輯,2003,33(4):335-343.

    [42]華仁民,陳培榮,張文蘭,等.論華南地區(qū)中生代3次大規(guī)模成礦作用[J].礦床地質(zhì),2005,24(2):99-107.

    [43]Campbell A R, Panter K S. Comparison of fluid inclusions in coexisting wolframite, cassiterite, and quartz from St.Michael’s mount and cligga head, Cornwall, England[J]. Geochimica Et Cosmochimica Acta, 1990, 54(3): 673-681.

    [44]蔡建明,劉若蘭,曾廣勝. 江西盤古山鎢礦流體包裹體及其與成礦關(guān)系的研究[C]//余鴻彰. 鎢礦地質(zhì)討論會論文集(中文版). 北京:地質(zhì)出版社,1981:1-11.

    [45]Bischoff J L. Densities of liquids and vapors in boiling NaCl-H2O solutions:A PVTX summary from 300℃ to 500℃[J]. American Journal of Science, 1991, 291(4): 309-338.

    [46]Bodnar R J. A method of calculateing fluid inclusion volumes based on vppor bubble diameters and PVTX properties of inclusion fluids[J]. Economic Geology, 1983, 78(3): 535-542.

    [47]Flincor B P. FLINCOR: A microcomputer program for the reduction and investigation of fluid-inclusiondata[J]. American Mineralogist, 1989, 74(11): 1390-1393.

    [48]穆治國,黃福生,陳成業(yè),等. 漂塘—西華山石英脈型鎢礦床碳、氫、氧穩(wěn)定同位素研究[C]//余鴻彰.鎢礦地質(zhì)討論會論文集(中文版). 北京:地質(zhì)出版社,1981,153-168.

    [49]莊龍池,林偉圣,謝霆煥.大吉山鎢礦的穩(wěn)定同位素地球化學[J].中國地質(zhì)科學院宜昌地質(zhì)礦產(chǎn)研究所所刊,1991(16):109-120.

    [50]吳永樂,梅勇文. 西華山鎢礦田多次成巖成礦及其演化規(guī)律[C]//余鴻彰.鎢礦地質(zhì)討論會論文集(中文版). 北京: 地質(zhì)出版社,1981,207-220.

    [51]馮志文,夏衛(wèi)華,章錦統(tǒng),等.江西黃沙脈鎢礦床特征及成礦流體性質(zhì)討論[J].地球科學:中國地質(zhì)大學學報,1989,14(4):423-432.

    [52]張理剛,莊龍池,錢雅倩,等.江西西華山-漂塘地區(qū)花崗巖及其鎢錫礦床的穩(wěn)定同位素地球化學[C]//鎢礦地質(zhì)討論會論文集,北京:地質(zhì)出版社,1981:325-338.

    [53]劉家齊.西華山花崗巖及其成礦作用[J].中國地質(zhì)科學院院報,1989,19:83-105.

    [54]張理剛.成巖成礦理論與找礦—中國主要類型礦床及花崗巖類巖石的穩(wěn)定同位素地質(zhì)學[M].北京:北京工業(yè)大學出版社,1989:1-200.

    [55]Taylor H P. The application of Oxygen and Hydrogen isotope studies to problems of hydrothermal alteration and ore depositions[J]. Economic Geology, 1974, 69(6): 843-883.

    [56]White D E. Diverse origins of hydrothermal ore Fluids[J]. Economic Geology, 1974, 69(6): 954-973.

    [57]丁悌平.氫氧同位素地球化學[M].北京:地質(zhì)出版社,1980:1-184.

    [58]王登紅,陳富文,張永忠,等.南嶺有色-貴金屬成礦潛力及綜合探測技術(shù)研究[M].北京:地質(zhì)出版社,2010:1-472.

    [59]郭春麗,王登紅,陳毓川,等.贛南中生代淘錫坑鎢礦區(qū)花崗巖鋯石SHRIMP年齡及石英脈Rb-Sr年齡測定[J].礦床地質(zhì),2007,26(4):432-442.

    [60]Giuliani G. Li Y D, Sheng T F. Fluid inclusion study of xihuashan Tungsten deposit in the southern Jiangxi province,China[J]. Mineralium Deposita, 1988, 23(1): 24-33.

    [61]Rios F J, Villas R N, Fuzikawa K. Fluid evolution in the Pedra Preta wolframite ore deposit, Paleoproterozoic Musa granite, eastern Amazon craton, Brazil[J]. Journal of South American Earth Sciences, 2003, 15(7): 787-802.

    [62]Beuchat S. Moritz R, Pettke T. Fluid evolution in the W-Cu-Zn-Pb San cristobal vein,Peru:fluid inclusion and stable isotope evidence[J]. Chemical Geology, 2004, 210: 201-224.

    [64]YANG Yongfei, LI Nuo, CHEN Yanjing. Fluid inclusion study of the Nannihu giant porphyry Mo-W deposit, Henan Province, China: Implications for the Nature of porphyry ore-fluid systems formed in a continental collision setting[J]. Ore Geology Reviews, 2012, 46: 83-94.

    [65]Mohamed A M. Evolution of mineralizing fluids of cassiterite-wolframite and fluorite deposits from Mueilha Tin mine area, Eastern Desert of Egypt, evidence from fluid inclusion[J]. Arabian Journal of Geosciences, 2013, 6(3): 775-782.

    [66] Moura A, Dória A , Neiva A.M.R, et al. Metallogenesis at the Carris W-Mo-Sn deposit (Gerês, Portugal):Constraints from fluid inclusions, mineral geochemistry, Re-Os and He-Ar isotopes[J]. Ore Geology Reviews, 2014,56(56):73-93.

    [67] Higgins N C. Fluid inclusion evidence for the transport of tungsten by carbonate complexes in hydrothermal solutions[J]. Canadian Journal of Earth Science. 2011, 17(7):823-830.

    [68]朱焱齡,李崇佑,林運淮.贛南鎢礦地質(zhì)[M].南昌:江西人民出版社,1981:324-337.

    [69]林新多,張德會,章傳玲.湖南宜章瑤崗仙黑鎢礦石英脈成礦流體性質(zhì)的探討[J].地球科學,1986,11(2):153-160.

    [70]張德會.石英脈型黑鎢礦床成礦流體性質(zhì)的進一步探討[J].地球科學,1987,12(2):185-192.

    [71]余行禎,李佩蘭.西華山鎢礦田成礦熱流體性質(zhì)分析[J].礦產(chǎn)與地質(zhì),1988(1):81-88.

    [72]干國良.江西黃沙脈鎢礦床成礦流體性質(zhì)的初步探討[J].宜昌地礦所???,1991:83-94.

    [73]常海亮,汪雄武,王曉地,等.西華山黑鎢礦-石英脈綠柱石中熔融包裹體的成分[J].巖石礦物學雜志,2007,26(3):259-268.

    [74]張德會.試論石英脈型黑鎢礦床的液態(tài)分離成因[J].地質(zhì)與勘探,1988(7):15-20.

    [75]Giamello M, Protano G, Riccobono F, et al. The W-Mo deposit of Perda Majori (SE Sardinia, Italy): a fluid inclusion study of ore and gangue minerals[J]. European Journal of Mineralogy, 1992, 4(5): 1079-1084.

    [76]Smith M P. Fluid inclusion and stable isotope constraints on the Genesis of the Cligga Head Sn-W deposit, S.W.England[J]. European Journal of Mineralogy, 1996, 8(5): 961-997.

    DOI:10.16788/j.hddz.32-1865/P.2016.02.008

    * 收稿日期:2015-07-11改回日期:2015-12-30責任編輯:汪建寧

    基金項目:國家自然科學基金項目(41302053),博士啟動資金(DHBK201120)資助。

    第一作者簡介:周龍全,1987年生,男,研究生,礦產(chǎn)普查與勘探專業(yè)。

    通訊作者簡介:李光來,1983年生,講師,碩士生導師,主要從事花崗巖與成礦作用研究。

    中圖分類號:P611

    文獻標識碼:A

    文章編號:2096-1871(2016)02-136-11

    Ore-forming fluids characteristics of quartz-vein tungsten deposits in southern Jiangxi Province

    ZHOU Long-quan, LI Guang-lai, TANG Ao, SU Ye

    (StateKeyLaboratoryBreedingBaseofNuclearResourcesandEnvironment,EastChinaInstituteofTechnology,Nanchang330013,China)

    Abstract:Southern Jiangxi is one important cluster of tungsten deposits in China, especially for hosting well developed quartz-vein tungsten deposits. Through analyzing the latest research achievements of fluid inclusion types, ore-forming fluids characteristics of the quartz-veins tungsten deposit, microthermometry and Raman spectroscopy, and combined with research results of C, H, O and Sr isotope compositions, this study discussed the characteristics of ore-forming fluids of quartz-vein tungsten deposits in southern Jiangxi Province, with an emphasis on the fluid evolution during the mineralization of quartz-vein tungsten deposits. The ore-forming fluids indicated by this study would originate from magmatic waters with the ore-forming fluid evolution beginning at magmatic hydrothermal transitional stage of high temperature and high salinity. The temperature and salinity of fluids related to wolframite mineralization are around 260~360℃ and 4~9 wt%NaCl eq., respectively, suggesting that the ore-forming fluids belonged to a mid-low salinity hydrothermal system enriched in SiO2, volatile components and ore-forming elements. The ore-forming mechanism in this area was dominated by fluid boiling and mixing, with minor deposits resulting from the ore precipitation due to fluid cooling.

    Key words:quartz-vein type tungsten deposit; ore-forming fluid; southern Jiangxi Province

    亚洲aⅴ乱码一区二区在线播放 | 国语自产精品视频在线第100页| 国产亚洲精品久久久久久毛片| 午夜福利18| 免费电影在线观看免费观看| 欧美黑人巨大hd| 久久草成人影院| 天堂av国产一区二区熟女人妻 | 岛国在线免费视频观看| 91麻豆av在线| 亚洲最大成人中文| 日韩精品免费视频一区二区三区| 免费观看人在逋| 成人手机av| 丰满人妻熟妇乱又伦精品不卡| 欧美乱色亚洲激情| 午夜福利18| 中文字幕熟女人妻在线| xxxwww97欧美| 人人妻人人看人人澡| 成人av一区二区三区在线看| 日韩精品免费视频一区二区三区| 日本精品一区二区三区蜜桃| 视频区欧美日本亚洲| 极品教师在线免费播放| 精品国产乱码久久久久久男人| 亚洲五月天丁香| 午夜精品久久久久久毛片777| 99热6这里只有精品| 国产成人啪精品午夜网站| 亚洲va日本ⅴa欧美va伊人久久| 国产av在哪里看| 亚洲av日韩精品久久久久久密| 两性夫妻黄色片| 亚洲欧美日韩高清专用| 黄色a级毛片大全视频| 亚洲精品av麻豆狂野| 成年版毛片免费区| 男女做爰动态图高潮gif福利片| 亚洲av熟女| www国产在线视频色| 欧美日韩国产亚洲二区| 五月伊人婷婷丁香| 久久久久久久精品吃奶| 校园春色视频在线观看| 国产精品1区2区在线观看.| 国产真人三级小视频在线观看| 国产精品免费一区二区三区在线| 在线免费观看的www视频| 精品久久久久久久毛片微露脸| 亚洲一码二码三码区别大吗| 久久精品成人免费网站| 成人手机av| 久久这里只有精品19| 欧美日韩亚洲国产一区二区在线观看| 人妻丰满熟妇av一区二区三区| 欧美乱色亚洲激情| 黄色 视频免费看| 精华霜和精华液先用哪个| 欧美大码av| 亚洲色图 男人天堂 中文字幕| 国产爱豆传媒在线观看 | 久久国产乱子伦精品免费另类| 午夜福利18| 国产精品免费一区二区三区在线| 中文字幕熟女人妻在线| 国产精品久久视频播放| 很黄的视频免费| 老司机福利观看| 中文字幕高清在线视频| 亚洲国产日韩欧美精品在线观看 | 欧美一区二区精品小视频在线| 婷婷精品国产亚洲av在线| 亚洲成人久久性| 国产精品免费视频内射| 国产视频一区二区在线看| 亚洲一区二区三区色噜噜| 美女大奶头视频| 亚洲熟妇中文字幕五十中出| 国产高清有码在线观看视频 | 中文亚洲av片在线观看爽| 国产私拍福利视频在线观看| 欧美人与性动交α欧美精品济南到| 一级毛片女人18水好多| 欧美又色又爽又黄视频| 69av精品久久久久久| 美女扒开内裤让男人捅视频| 全区人妻精品视频| 90打野战视频偷拍视频| 成熟少妇高潮喷水视频| 精品不卡国产一区二区三区| 亚洲熟妇熟女久久| 欧美3d第一页| 小说图片视频综合网站| 久久久久久久午夜电影| 老司机福利观看| 亚洲精品一卡2卡三卡4卡5卡| 两个人免费观看高清视频| 99riav亚洲国产免费| 最新美女视频免费是黄的| 久久婷婷成人综合色麻豆| 变态另类丝袜制服| 欧美3d第一页| 婷婷精品国产亚洲av| 精品国产美女av久久久久小说| 18禁裸乳无遮挡免费网站照片| xxx96com| 又粗又爽又猛毛片免费看| 久久久精品欧美日韩精品| 亚洲免费av在线视频| 欧美绝顶高潮抽搐喷水| 在线永久观看黄色视频| 制服诱惑二区| 久久久久免费精品人妻一区二区| 亚洲精品中文字幕在线视频| 欧美日韩黄片免| 精品欧美国产一区二区三| 亚洲欧美激情综合另类| 久久国产精品人妻蜜桃| 日韩国内少妇激情av| 国内少妇人妻偷人精品xxx网站 | 午夜两性在线视频| 亚洲av成人不卡在线观看播放网| 国产免费av片在线观看野外av| 少妇裸体淫交视频免费看高清 | 看片在线看免费视频| 一边摸一边做爽爽视频免费| 久久精品影院6| 99热6这里只有精品| 91九色精品人成在线观看| 一级a爱片免费观看的视频| 久久天躁狠狠躁夜夜2o2o| 国内精品久久久久精免费| 亚洲全国av大片| 久久人妻福利社区极品人妻图片| 欧美国产日韩亚洲一区| 妹子高潮喷水视频| 在线观看美女被高潮喷水网站 | 亚洲欧美日韩无卡精品| 久久久久国内视频| 亚洲av中文字字幕乱码综合| 国产精品av视频在线免费观看| 亚洲成av人片免费观看| 欧美日韩国产亚洲二区| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲va日本ⅴa欧美va伊人久久| 在线永久观看黄色视频| 国产日本99.免费观看| 又黄又爽又免费观看的视频| 日本精品一区二区三区蜜桃| 老司机午夜福利在线观看视频| 国产乱人伦免费视频| 免费在线观看完整版高清| 欧美日韩福利视频一区二区| 久久精品亚洲精品国产色婷小说| 久久婷婷人人爽人人干人人爱| 两性夫妻黄色片| 亚洲午夜理论影院| 久久天堂一区二区三区四区| 很黄的视频免费| 老熟妇乱子伦视频在线观看| 88av欧美| 国产成人欧美在线观看| 不卡一级毛片| 亚洲欧美一区二区三区黑人| 日韩欧美国产在线观看| 国产三级在线视频| 日韩高清综合在线| 又黄又粗又硬又大视频| 精品国产美女av久久久久小说| 免费在线观看黄色视频的| 亚洲一区二区三区不卡视频| 在线永久观看黄色视频| 18禁黄网站禁片午夜丰满| 日本熟妇午夜| 亚洲精品美女久久久久99蜜臀| 最新在线观看一区二区三区| 国产午夜福利久久久久久| 淫秽高清视频在线观看| 蜜桃久久精品国产亚洲av| 欧美日韩黄片免| 我要搜黄色片| 91大片在线观看| 久久久久亚洲av毛片大全| cao死你这个sao货| 天堂影院成人在线观看| 妹子高潮喷水视频| 老司机福利观看| 男人的好看免费观看在线视频 | 中文字幕av在线有码专区| 白带黄色成豆腐渣| 欧美av亚洲av综合av国产av| 最近视频中文字幕2019在线8| 狂野欧美白嫩少妇大欣赏| 淫妇啪啪啪对白视频| 黄色 视频免费看| 国产私拍福利视频在线观看| 波多野结衣高清作品| 精品久久久久久,| 日本一本二区三区精品| 在线观看舔阴道视频| 日韩欧美 国产精品| 天堂√8在线中文| 制服诱惑二区| 怎么达到女性高潮| 老汉色∧v一级毛片| 天堂动漫精品| 欧美黄色淫秽网站| 我的老师免费观看完整版| 国产精品自产拍在线观看55亚洲| 九色国产91popny在线| 午夜福利在线在线| tocl精华| 777久久人妻少妇嫩草av网站| 色综合站精品国产| 国产精品 欧美亚洲| www.自偷自拍.com| 毛片女人毛片| 午夜激情av网站| 床上黄色一级片| 欧美成人免费av一区二区三区| 老汉色∧v一级毛片| 亚洲一码二码三码区别大吗| 一级毛片女人18水好多| 国产又色又爽无遮挡免费看| aaaaa片日本免费| 全区人妻精品视频| 母亲3免费完整高清在线观看| 久久 成人 亚洲| 两个人免费观看高清视频| 精品一区二区三区四区五区乱码| 亚洲电影在线观看av| 午夜亚洲福利在线播放| 亚洲精品av麻豆狂野| 婷婷六月久久综合丁香| 欧美中文日本在线观看视频| ponron亚洲| 18美女黄网站色大片免费观看| 午夜日韩欧美国产| 国产av不卡久久| 校园春色视频在线观看| 亚洲美女黄片视频| 久久这里只有精品19| 亚洲国产高清在线一区二区三| 少妇裸体淫交视频免费看高清 | 成人国产一区最新在线观看| 午夜福利在线观看吧| 国产午夜精品论理片| 亚洲精品色激情综合| 成熟少妇高潮喷水视频| 亚洲自拍偷在线| 最近在线观看免费完整版| 亚洲av五月六月丁香网| 国产伦人伦偷精品视频| 丁香六月欧美| 亚洲男人天堂网一区| 成熟少妇高潮喷水视频| 天天躁狠狠躁夜夜躁狠狠躁| 免费在线观看成人毛片| 男人舔女人的私密视频| 欧美人与性动交α欧美精品济南到| 午夜激情av网站| av在线天堂中文字幕| 日韩大尺度精品在线看网址| 免费在线观看完整版高清| 色综合欧美亚洲国产小说| 成人国产一区最新在线观看| 亚洲欧美日韩高清专用| 久热爱精品视频在线9| 国产精品自产拍在线观看55亚洲| 精品久久久久久久人妻蜜臀av| 丰满的人妻完整版| 国产精品久久电影中文字幕| 精华霜和精华液先用哪个| 三级毛片av免费| 亚洲精品一卡2卡三卡4卡5卡| 99在线人妻在线中文字幕| e午夜精品久久久久久久| 欧美性长视频在线观看| 美女午夜性视频免费| 精品免费久久久久久久清纯| 亚洲第一欧美日韩一区二区三区| 国语自产精品视频在线第100页| 精品国产美女av久久久久小说| 亚洲欧美日韩高清专用| 波多野结衣巨乳人妻| 免费看a级黄色片| 91国产中文字幕| 天天一区二区日本电影三级| 少妇裸体淫交视频免费看高清 | 亚洲,欧美精品.| 成熟少妇高潮喷水视频| 欧美日本亚洲视频在线播放| 久久久精品国产亚洲av高清涩受| 国产午夜福利久久久久久| 男女视频在线观看网站免费 | 久久久久国内视频| 亚洲欧美精品综合一区二区三区| 久久久久久久精品吃奶| 久久久久国产精品人妻aⅴ院| 精品久久久久久久毛片微露脸| 亚洲 欧美一区二区三区| 宅男免费午夜| 亚洲av电影不卡..在线观看| 国产一区二区在线观看日韩 | 国产伦人伦偷精品视频| 亚洲专区字幕在线| 国产亚洲精品久久久久5区| 麻豆一二三区av精品| 五月玫瑰六月丁香| 九色国产91popny在线| 99久久精品热视频| www.熟女人妻精品国产| 国产又色又爽无遮挡免费看| 亚洲avbb在线观看| 色噜噜av男人的天堂激情| 成年免费大片在线观看| 亚洲美女黄片视频| 国产精品亚洲av一区麻豆| 男人舔女人的私密视频| 三级毛片av免费| 国产高清视频在线播放一区| 久99久视频精品免费| 老汉色∧v一级毛片| 男人的好看免费观看在线视频 | 狂野欧美白嫩少妇大欣赏| 欧美日韩黄片免| 日韩大码丰满熟妇| 精品不卡国产一区二区三区| a在线观看视频网站| 亚洲欧美精品综合一区二区三区| 久久天躁狠狠躁夜夜2o2o| 婷婷亚洲欧美| 免费人成视频x8x8入口观看| 午夜激情av网站| 69av精品久久久久久| 不卡av一区二区三区| 日韩欧美免费精品| 无限看片的www在线观看| 听说在线观看完整版免费高清| 麻豆一二三区av精品| 九色成人免费人妻av| 精品久久久久久久久久免费视频| 天堂动漫精品| 国产精品电影一区二区三区| 久久精品综合一区二区三区| 亚洲国产欧美网| 亚洲精品国产精品久久久不卡| 亚洲av成人精品一区久久| 丰满人妻一区二区三区视频av | 国产亚洲av嫩草精品影院| 女同久久另类99精品国产91| 午夜福利视频1000在线观看| 日本黄色视频三级网站网址| 两性夫妻黄色片| 精品国产亚洲在线| 麻豆成人午夜福利视频| 一区二区三区高清视频在线| 三级毛片av免费| 好男人电影高清在线观看| 欧美日韩一级在线毛片| 亚洲av成人不卡在线观看播放网| 成年人黄色毛片网站| 中文在线观看免费www的网站 | 久久午夜亚洲精品久久| 午夜日韩欧美国产| 丰满的人妻完整版| 亚洲成a人片在线一区二区| 一级黄色大片毛片| 露出奶头的视频| 在线播放国产精品三级| 啦啦啦免费观看视频1| 婷婷丁香在线五月| 久久99热这里只有精品18| 日本在线视频免费播放| 精品乱码久久久久久99久播| 久久久精品大字幕| 国产一区二区在线av高清观看| 日日爽夜夜爽网站| 亚洲性夜色夜夜综合| 亚洲国产日韩欧美精品在线观看 | 亚洲国产欧洲综合997久久,| 久久久久久久久免费视频了| 麻豆国产97在线/欧美 | 变态另类丝袜制服| 日韩欧美在线乱码| 91av网站免费观看| 白带黄色成豆腐渣| 国产成人影院久久av| 全区人妻精品视频| 久久久国产成人免费| 少妇粗大呻吟视频| 亚洲欧美激情综合另类| 丁香欧美五月| 欧美中文综合在线视频| 免费在线观看完整版高清| 精品久久久久久久久久免费视频| 国产精品九九99| 99国产精品99久久久久| 免费看美女性在线毛片视频| 在线十欧美十亚洲十日本专区| 亚洲aⅴ乱码一区二区在线播放 | 99久久国产精品久久久| 精品人妻1区二区| 国产97色在线日韩免费| 国产精品香港三级国产av潘金莲| 日韩欧美在线二视频| 精品国产美女av久久久久小说| 在线播放国产精品三级| 一二三四社区在线视频社区8| 精品国产乱子伦一区二区三区| 国产午夜精品久久久久久| 亚洲av中文字字幕乱码综合| 免费看美女性在线毛片视频| 曰老女人黄片| 18禁裸乳无遮挡免费网站照片| 999久久久精品免费观看国产| 国产av又大| 999久久久国产精品视频| 国产成人一区二区三区免费视频网站| 天天躁狠狠躁夜夜躁狠狠躁| 国产久久久一区二区三区| 久久婷婷人人爽人人干人人爱| 欧美大码av| 亚洲av成人av| 亚洲午夜理论影院| 久久久久久久精品吃奶| 国产精品久久久久久久电影 | 国产成人影院久久av| 国产99白浆流出| 在线观看日韩欧美| 久久人人精品亚洲av| 午夜福利成人在线免费观看| 18禁国产床啪视频网站| 久久香蕉国产精品| 亚洲国产中文字幕在线视频| 亚洲自拍偷在线| 国产麻豆成人av免费视频| 亚洲精品在线美女| av超薄肉色丝袜交足视频| 欧美绝顶高潮抽搐喷水| 床上黄色一级片| 国产又黄又爽又无遮挡在线| 在线观看一区二区三区| 国产日本99.免费观看| 一区二区三区高清视频在线| 黄片大片在线免费观看| 老司机在亚洲福利影院| 露出奶头的视频| 日日夜夜操网爽| 美女黄网站色视频| 男女那种视频在线观看| 一二三四社区在线视频社区8| 麻豆一二三区av精品| 免费观看人在逋| 国产精品99久久99久久久不卡| 桃红色精品国产亚洲av| 天天躁夜夜躁狠狠躁躁| 91麻豆av在线| 宅男免费午夜| 在线国产一区二区在线| 国产高清有码在线观看视频 | 日日干狠狠操夜夜爽| 精品不卡国产一区二区三区| 可以在线观看的亚洲视频| 国产三级在线视频| 色尼玛亚洲综合影院| 亚洲国产欧洲综合997久久,| 女人爽到高潮嗷嗷叫在线视频| 久久精品综合一区二区三区| 国产精品野战在线观看| 亚洲av成人av| 亚洲一码二码三码区别大吗| 淫妇啪啪啪对白视频| 色噜噜av男人的天堂激情| 夜夜夜夜夜久久久久| 嫩草影视91久久| 久久人人精品亚洲av| 亚洲成人久久爱视频| 两个人的视频大全免费| 国产成人影院久久av| 国产乱人伦免费视频| 露出奶头的视频| 成人精品一区二区免费| 长腿黑丝高跟| 啦啦啦观看免费观看视频高清| 黄片小视频在线播放| 国产精品永久免费网站| 色噜噜av男人的天堂激情| 老司机福利观看| 一a级毛片在线观看| 制服人妻中文乱码| or卡值多少钱| 精品久久久久久成人av| 欧美黑人欧美精品刺激| 欧美日本视频| 国产精品 国内视频| 国产精品亚洲一级av第二区| 欧美黑人巨大hd| 欧美日韩精品网址| 亚洲七黄色美女视频| 日本熟妇午夜| 一个人免费在线观看电影 | 九色国产91popny在线| 午夜免费观看网址| 午夜久久久久精精品| 国产精品免费一区二区三区在线| 男女那种视频在线观看| 亚洲精品色激情综合| 国内精品久久久久精免费| 国产三级在线视频| 国产一区二区在线av高清观看| 级片在线观看| www日本黄色视频网| 久久久久性生活片| 一区二区三区高清视频在线| 18禁裸乳无遮挡免费网站照片| www日本在线高清视频| 黄色视频,在线免费观看| 日日夜夜操网爽| 亚洲av成人av| 国产探花在线观看一区二区| 88av欧美| 国内精品久久久久精免费| 长腿黑丝高跟| 黄色毛片三级朝国网站| 日韩三级视频一区二区三区| 久久 成人 亚洲| 日本精品一区二区三区蜜桃| 欧美绝顶高潮抽搐喷水| 制服诱惑二区| 美女午夜性视频免费| 国产探花在线观看一区二区| 大型av网站在线播放| 老汉色∧v一级毛片| 91麻豆av在线| 19禁男女啪啪无遮挡网站| 久久久久国产精品人妻aⅴ院| 久久久国产成人免费| 日韩 欧美 亚洲 中文字幕| 怎么达到女性高潮| 亚洲aⅴ乱码一区二区在线播放 | 国产欧美日韩一区二区精品| 美女免费视频网站| 国产亚洲av嫩草精品影院| 久久久国产欧美日韩av| 亚洲av成人精品一区久久| 国产亚洲精品久久久久5区| 精品国产超薄肉色丝袜足j| 国产单亲对白刺激| 长腿黑丝高跟| 成年版毛片免费区| 久久这里只有精品19| 国产欧美日韩一区二区三| 韩国av一区二区三区四区| 日本一二三区视频观看| 亚洲成人精品中文字幕电影| 国产在线精品亚洲第一网站| 搡老妇女老女人老熟妇| 999精品在线视频| 亚洲欧美精品综合久久99| 床上黄色一级片| 欧美日韩精品网址| 亚洲欧美激情综合另类| 黄片小视频在线播放| 欧美 亚洲 国产 日韩一| 日韩高清综合在线| 日韩欧美 国产精品| 国产日本99.免费观看| 黄色a级毛片大全视频| 亚洲精品久久成人aⅴ小说| 色尼玛亚洲综合影院| 18禁国产床啪视频网站| 91老司机精品| 视频区欧美日本亚洲| 国产伦在线观看视频一区| 大型av网站在线播放| 久久香蕉精品热| 亚洲人成电影免费在线| www.自偷自拍.com| 亚洲第一欧美日韩一区二区三区| 99久久精品热视频| 精品久久久久久久久久免费视频| 亚洲精品粉嫩美女一区| 国产亚洲精品av在线| 国内毛片毛片毛片毛片毛片| 亚洲欧美日韩东京热| 两个人视频免费观看高清| 精品少妇一区二区三区视频日本电影| 欧美成狂野欧美在线观看| 嫩草影视91久久| 欧美成人免费av一区二区三区| 色综合亚洲欧美另类图片| 免费在线观看日本一区| 国产精品av久久久久免费| 成人国产综合亚洲| 男人的好看免费观看在线视频 | 亚洲欧美激情综合另类| 99精品欧美一区二区三区四区| 少妇熟女aⅴ在线视频| 无人区码免费观看不卡| 免费在线观看完整版高清| 国产精品一区二区三区四区免费观看 | 伊人久久大香线蕉亚洲五| 午夜久久久久精精品| 亚洲精品中文字幕一二三四区| 好男人电影高清在线观看| xxx96com| 99久久无色码亚洲精品果冻| 久久久久亚洲av毛片大全| 亚洲欧美日韩高清在线视频| av超薄肉色丝袜交足视频| 亚洲欧美日韩无卡精品| 亚洲av中文字字幕乱码综合| 午夜福利视频1000在线观看|