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

    貴州都勻撒朗重晶石礦床流體包裹體及H-O-C同位素地球化學(xué)研究

    2021-06-29 11:11:44何明勤李應(yīng)輝田歡歡梁坤萍
    大地構(gòu)造與成礦學(xué) 2021年3期
    關(guān)鍵詞:成礦

    何明勤, 李應(yīng)輝, 張 豐, 田歡歡, 梁坤萍

    貴州都勻撒朗重晶石礦床流體包裹體及H-O-C同位素地球化學(xué)研究

    何明勤1, 2*, 李應(yīng)輝3, 張 豐1, 田歡歡1, 梁坤萍1

    (1.貴州大學(xué) 資源與環(huán)境工程學(xué)院, 貴州 貴陽 550025; 2.貴州大學(xué) 喀斯特地質(zhì)資源與環(huán)境教育部重點(diǎn)實驗室, 貴州 貴陽 550025; 3.貴州省地質(zhì)礦產(chǎn)勘查開發(fā)局 105地質(zhì)大隊, 貴州 貴陽 550018)

    貴州都勻撒朗重晶石礦床位于江南造山帶西南段都勻陸緣盆地區(qū)之都勻南北向構(gòu)造變形區(qū), 礦體主要呈層狀、似層狀、脈狀或囊狀產(chǎn)于層間破碎帶及構(gòu)造交匯部位, 受構(gòu)造控制明顯。本文在野外礦床地質(zhì)調(diào)查的基礎(chǔ)上, 進(jìn)行了流體包裹體巖相學(xué)、顯微測溫和群體包裹體成分及H-O-C同位素研究, 結(jié)果顯示: 撒朗重晶石礦床的成礦流體為具有中低溫(95~240 ℃)、峰值溫度(140 ℃)為低溫、低鹽度(均值4.9%~9.1% NaCleqv)、中等密度(均值0.95~0.97 g/cm3)特征的Ca2+-Na+-Cl?-SO42?類型水; 重晶石流體包裹體水H同位素組成為?48‰~?38‰、O同位素組成為0.55‰~12.05‰, 顯示成礦流體水來源于大氣降水或與圍巖水巖反應(yīng)后的演化大氣降水; 礦石中與重晶石共生的方解石C同位素組成為?1.9‰~?1.0‰、O同位素組成(SMOW標(biāo)準(zhǔn))為20.40‰~21.01‰, 指示碳來自于海相碳酸鹽巖的溶解; 礦床的形成晚于圍巖, 是在圍巖地層形成之后, 為后期構(gòu)造及熱液活動使其中的成礦物質(zhì)活化遷移并在相應(yīng)構(gòu)造部位富集而形成。

    流體包裹體; H-O-C同位素; 礦床地球化學(xué); 撒朗重晶石礦床; 貴州都勻

    0 引 言

    重晶石礦床廣泛產(chǎn)于太古宇到第四系不同時代地層之中(李文炎和余洪云, 1991), 分布范圍廣, 規(guī)模大, 具有非常重要的經(jīng)濟(jì)價值和研究意義(高軍波和楊瑞東, 2015)。貴州的重晶石礦資源尤為豐富, 產(chǎn)出層位也十分多樣, 具有重要工業(yè)價值的重晶石礦床主要產(chǎn)于寒武系、泥盆系和奧陶系。其中, 寒武系代表性礦床有天柱縣大河邊超大型重晶石礦床, 其資源儲量超兩億噸(楊瑞東等, 2007a), 是世界上資源儲量最大的重晶石礦床(楊瑞東等, 2007a; Han et al., 2015), 其經(jīng)濟(jì)價值大, 研究程度高, 一般認(rèn)為是噴流熱水沉積成因(吳朝東等, 1999; 方維萱等, 2002; 夏菲等, 2004; 楊瑞東等, 2007a, 2007b; 侯東壯等, 2015; Zhou et al., 2018)。泥盆系重晶石礦床以鎮(zhèn)寧?紫云一帶的最為重要, 如樂紀(jì)大型重晶石礦床等。近些年來, 一些學(xué)者對其進(jìn)行了Sr同位素組成(高軍波等, 2013a)、礦床地球化學(xué)特征及成因(高軍波等, 2013b)和礦床地質(zhì)特征與成礦規(guī)律(張旭等, 2017)研究, 認(rèn)為礦床是斷裂控制的海底噴流熱水沉積成礦的產(chǎn)物(高軍波和楊瑞東, 2015)或甲烷滲漏成因所形成(Gao et al., 2017)。奧陶系重晶石礦床以施秉縣頂罐坡大型重晶石礦床最為典型, 根據(jù)流體包裹體測溫數(shù)據(jù), 認(rèn)為礦床屬低溫?zé)嵋旱V床(貴州省區(qū)域礦產(chǎn)志編寫組, 1986)。與主要產(chǎn)于高角度張性斷層內(nèi)的陡傾斜施秉縣頂罐坡重晶石礦床(圖1a)不同, 同樣賦存于奧陶系的都勻市王司鎮(zhèn)撒朗村撒朗重晶石礦床礦體主要呈緩傾斜的層狀、似層狀、脈狀或囊狀。同時, 該礦床與國內(nèi)外一些主要層狀的沉積重晶石礦床如天柱縣大河邊重晶石礦床(夏菲等, 2005; Zhou et al., 2018)和土耳其Dinek重晶石礦床(Elmas et al., 2012)等含硫化物礦物也明顯不同, 該礦床礦石中不含硫化物, 從而品位相對更高、純度相對更好, 因而更顯其獨(dú)特開發(fā)利用和研究價值。Maynard and Okita (1991)基于重晶石礦床產(chǎn)出的構(gòu)造背景將層狀重晶石礦床劃分為大陸邊緣型和克拉通裂谷型兩種, 前者以單一重晶石礦床為主, 后者還兼具鉛、鋅礦化, 據(jù)此可將撒朗重晶石礦床歸入大陸邊緣型層狀重晶石礦床。撒朗重晶石礦區(qū)內(nèi), 礦床(體)出露較好, 開采方便, 但現(xiàn)有的勘查和研究程度均很低, 限制了進(jìn)一步找礦和規(guī)?;_采工作。目前, 已有研究成果僅限于對礦床地質(zhì)特征和控礦因素的分析(李應(yīng)輝等, 2016), 缺乏對礦床地球化學(xué)尤其是成礦流體地球化學(xué)方面的相關(guān)研究。基于此, 本論文擬通過流體包裹體、同位素地球化學(xué)等手段對礦床成礦流體特征及其來源進(jìn)行較系統(tǒng)研究, 為礦床成因分析提供依據(jù), 并為深入總結(jié)成礦規(guī)律、找礦和開發(fā)利用工作提供思路。

    1 區(qū)域地質(zhì)及礦床地質(zhì)特征

    撒朗重晶石礦床位于江南造山帶西南段(圖1a)都勻陸緣盆地區(qū)之都勻南北向構(gòu)造變形區(qū)(貴州省地質(zhì)調(diào)查院, 2012)。區(qū)域內(nèi)出露地層主要有古生界寒武系、奧陶系、志留系和泥盆系, 其中奧陶系是區(qū)域重晶石、鉛鋅礦床等的主要賦礦層位, 巖性以白云巖和石灰?guī)r為主。區(qū)域內(nèi)構(gòu)造以近南北向為主, 褶皺構(gòu)造以王司背斜為主體, 斷裂構(gòu)造除南北向斷裂外, 還發(fā)育北東向、北西向和近東西向斷裂, 這些斷裂切割區(qū)內(nèi)褶皺與地層, 構(gòu)成復(fù)雜的構(gòu)造格局, 控制了區(qū)域多種礦產(chǎn)的形成(李應(yīng)輝等, 2016), 如都勻牛角塘鋅礦床(陳國勇等, 1992; 張碧志等, 1994)、鎘礦床(葉霖和劉鐵庚, 1997; 劉鐵庚和葉霖, 2000)或鎘鋅礦床(劉鐵庚等, 2004)、老貓石鉛鋅礦床(劉強(qiáng)和夏磊, 2013)、大梁子鋅礦床(劉成萬, 2017)和菜園河鋅礦床(王瓊等, 2017)等。

    1. 第四系; 2. 中下志留統(tǒng)翁項群; 3. 下奧陶統(tǒng)紅花園組+大灣組; 4. 下奧陶統(tǒng)桐梓組; 5. 正斷層及編號; 6. 逆斷層及編號; 7. 性質(zhì)不明斷層及編號; 8. 地層界線; 9. 角度不整合界線; 10. 礦體及編號; 11. 地層產(chǎn)狀; 12. 剖面位置及編號; 13. 頂罐坡重晶石礦床; 14. 礦體主要樣品采集位置及序號。

    礦區(qū)出露地層有下奧陶統(tǒng)桐梓組(O1)、紅花園組+大灣組(O1)、中下志留統(tǒng)翁項群(S1-2)和第四系(Q)(圖1b)。其中, 下奧陶統(tǒng)桐梓組(O1)和紅花園組+大灣組(O1)是一套主要由白云巖和石灰?guī)r組成的碳酸鹽類巖石; 中下志留統(tǒng)翁項群(S1-2)僅出露上亞群, 主要是由砂巖和頁巖組成的碎屑巖; 第四系(Q)僅出露在礦區(qū)西南部, 主要為黏土、粉砂、砂土或含礫砂土等。

    礦區(qū)構(gòu)造總體為一單斜構(gòu)造, 傾向西(260°~280°), 傾角5°~20°。因受斷裂構(gòu)造的影響, 局部地層向北、東傾斜。區(qū)內(nèi)斷裂構(gòu)造相對較發(fā)育, 可分為近南北向(F1)、北西向(F2、F3、F4)和北東向(F5)三組。

    礦體產(chǎn)于下奧陶統(tǒng)桐梓組(O1)中, 主要受層間斷裂破碎帶控制。據(jù)礦體露頭和采礦工程揭露, 目前礦區(qū)主要產(chǎn)出有Ⅰ、Ⅱ、Ⅲ共3個礦體(圖1b)。Ⅰ號礦體呈似層狀(圖2), 走向長約500 m、厚度1.93~2.08 m、傾向延伸大于100 m; Ⅱ號礦體呈層狀、似層狀, 走向長約132 m、厚度2.5~2.8 m、傾向延伸大于15 m; Ⅲ號礦體呈似層狀(圖2)、脈狀和囊狀產(chǎn)出, 走向長約100 m、厚度1~2 m、傾向延伸大于12 m。

    礦石成分相對簡單, 礦石礦物為重晶石, 脈石礦物有方解石、少量的石英和白云石; 礦石中BaSO4含量為92.50%~96.18%, 其次為SiO2, 含量為0.41%~ 1.30%。礦石構(gòu)造以條帶狀、脈狀構(gòu)造為主, 其次為花斑狀及團(tuán)塊狀構(gòu)造(圖3); 礦石結(jié)構(gòu)以放射狀、束狀和粒狀結(jié)構(gòu)為主, 其次為網(wǎng)狀、交代殘余結(jié)構(gòu)(圖4)。礦區(qū)圍巖蝕變相對較弱, 主要是方解石化, 其次為硅化。

    2 樣品采集與分析測試方法

    在野外較詳細(xì)地質(zhì)觀察的基礎(chǔ)上, 對出露或工程揭露較好的Ⅰ、Ⅱ、Ⅲ三個礦體的不同類型礦石進(jìn)行了較系統(tǒng)采樣, 共采集礦石樣品18件, 全部磨制包裹體片, 從中選擇透明度較好包裹體片進(jìn)行流體包裹體巖相學(xué)觀察和對其中的原生包裹體進(jìn)行溫度測定, 并挑選出重晶石、方解石單礦物用于流體包裹體成分和H-O-C同位素分析。

    流體包裹體巖相學(xué)和溫度測定在貴州大學(xué)資源與環(huán)境工程學(xué)院流體包裹體實驗室進(jìn)行, 測試儀器為英國產(chǎn)Linkam THMSG600型冷熱兩用臺和配套的Leica透光顯微鏡DMLP及其LinkSys溫控系統(tǒng), 儀器溫度測試范圍為?196~600 ℃, 冷凍溫度和均一溫度的測試精度分別為±0.1 ℃和±1 ℃。

    流體包裹體氣、液相成分分析在核工業(yè)北京地質(zhì)研究院分析測試研究中心完成, 操作流程是: 將富含原生流體包裹體的礦石碎至60~80目, 在雙目顯微鏡下選純重晶石礦物顆粒, 分階段爆裂, 先在較低溫度下爆裂除去可能的次生包裹體, 然后在300~350 ℃爆裂條件下爆裂。氣相成分分析采用CLarus600氣相色譜儀, 液相成分分析采用DIONEX-500離子色譜儀, 測定的室溫控制在20 ℃, 濕度控制在20%, 分析的液相陽、陰離子成分為Na+、K+、Mg2+、Ca2+、F?、Cl?、NO3?、SO42?, 氣相分析成分為H2、N2、CO、CH4、CO2、H2O, 其液相陽、陰離子和氣相成分分析嚴(yán)格按國家或行業(yè)相關(guān)標(biāo)準(zhǔn)進(jìn)行。

    1. 第四系粉砂礫石; 2. 紅花園組+大灣組泥質(zhì)灰?guī)r; 3. 桐梓組白云巖; 4. 翁項群砂巖; 5. 地層界線; 6. 逆斷層及編號; 7. 正斷層及編號; 8. 礦體及編號; 9. 剝土工程及編號。地層代號: Q. 第四系; S1-2W. 中下志留統(tǒng)翁項群; O1h+d. 下奧陶統(tǒng)紅花園組+大灣組; O1t. 下奧陶統(tǒng)桐梓組。

    (a) 條帶狀構(gòu)造; (b) 脈狀構(gòu)造; (c) 花斑狀構(gòu)造; (d) 團(tuán)塊狀構(gòu)造。礦物/巖石代號: Brt. 重晶石; Dol. 白云巖。

    H-O-C同位素分析首先將肉眼挑選出來的富含重晶石和方解石的礦石破碎至10~20目, 去除明顯的雜質(zhì)組分, 再將其置于雙目顯微鏡下, 挑選出純的重晶石和方解石單礦物, 純度大于99%。流體包裹體水H同位素分析在澳實(廣州)分析測試中心進(jìn)行, 礦物O及O-C同位素分析在中國科學(xué)院地質(zhì)與地球物理研究所穩(wěn)定同位素分析實驗室進(jìn)行。具體操作流程如下所述。

    H同位素分析流程: 先將挑選出來的重晶石單礦物置于長條狀容器中稱重, 然后待樣品自動進(jìn)入取樣器后在100 ℃下保持1小時, 其次采用爆裂法取包裹體水、鋅法制氫, 最后用同位素比值質(zhì)譜(TC-EA-IRMS)在線測定流體包裹體水的H同位素比值, 分析結(jié)果以維也納平均海水(SMOW)為標(biāo)準(zhǔn), 分析精度為3‰。

    重晶石O同位素分析: 采用BrF5法, 在真空條件下, 將重晶石與BrF5進(jìn)行高溫氧化還原反應(yīng)生成O2, 將氧氣提純, 在溫度為700 ℃的高溫條件下, O2與石墨發(fā)生反應(yīng), 最終轉(zhuǎn)化為CO2, 從而進(jìn)行O同位素質(zhì)譜分析, 分析結(jié)果采用相對國際標(biāo)準(zhǔn)SMOW之值, 分析儀器為賽默飛公司的DELTA V質(zhì)譜儀, 分析誤差在0.2‰以內(nèi)。

    方解石C-O同位素分析: 采用100%磷酸法, 將分析樣品方解石與磷酸發(fā)生化學(xué)反應(yīng), 反應(yīng)過程中會釋放CO2, 將釋放出的CO2送入MAT-253型質(zhì)譜儀上進(jìn)行C、O同位素測定, 所測數(shù)據(jù)均采用相對國際標(biāo)準(zhǔn)PDB之值, C同位素分析精度優(yōu)于0.15‰、O同位素分析精度優(yōu)于0.2‰。

    3 分析測試結(jié)果

    3.1 流體包裹體

    3.1.1 流體包裹體巖相學(xué)特征

    在透光顯微鏡下, 對磨制的18件樣品包裹體片進(jìn)行了較詳細(xì)觀察后發(fā)現(xiàn), 礦石中重晶石和方解石內(nèi)的包裹體較為發(fā)育, 以原生包裹體為主, 但個體普遍較小, 一般小于10 μm, 多數(shù)為1~6 μm, 按照室溫(25 ℃)(Mapani and Wilson, 1998)下的成分、相態(tài)特征及其相互關(guān)系(何明勤等, 2004a), 尤其是室溫下的相態(tài)關(guān)系的包裹體分類方案(何明勤等, 2004b), 將該礦區(qū)原生包裹體分為單一液相的純液體包裹體(Ⅰ型)和氣液兩相的液體包裹體(Ⅱ型)兩類(圖5)。

    (1) Ⅰ型(圖5a~b): 純液體包裹體, 該類包裹體無色透明, 由單一的液相L所組成, 氣液比為0, 形態(tài)為短透鏡狀、短柱狀或不規(guī)則狀, 呈孤立或稀疏分布, 分布稀少, 數(shù)量少于5%。

    (2) Ⅱ型(圖5c~f): 液體包裹體, 由氣相V、液相L兩相組成, 氣、液相均為無色, 氣液比3%~20%, 形態(tài)為負(fù)晶形、長條形、三角形、橢圓形、平行四邊形或不規(guī)則狀, 多呈孤立狀、零星狀或稀疏狀分布, 少數(shù)沿晶體生長方向成群分布, 是礦區(qū)最主要的包裹體類型, 分布廣泛, 數(shù)量大于95%。

    3.1.2 溫度測定

    流體包裹體溫度的測定按照先測冷凍溫度后測均一溫度的順序進(jìn)行, 由于不同礦體或不同類型的重晶石或方解石等礦物中, Ⅱ型液體包裹體是數(shù)量最多的包裹體, 因此無論是冷凍溫度還是均一溫度均是通過測定Ⅱ型液體包裹體來獲得的, 且均一時的相態(tài)均為液相。

    礦區(qū)內(nèi)不同礦體重晶石(Ⅱ號礦體還包括少量共生方解石)的流體包裹體溫度測定結(jié)果見表1, 結(jié)果顯示: Ⅰ號礦體均一溫度為95~240 ℃、均值為130~166 ℃; Ⅱ號礦體均一溫度為117~240 ℃、均值123~134 ℃; Ⅲ號礦體均一溫度為103~221 ℃, 均值為149 ℃。在均一溫度分布直方圖(圖6a)上, 均一溫度呈單峰態(tài)分布, 位于中低溫范圍, 峰值溫度為140 ℃, 表明其具較單一的中低溫成礦特點(diǎn), 且成礦高峰溫度為低溫。

    (a) 針狀重晶石由中心向外排列呈放射狀結(jié)構(gòu); (b) 針柱狀重晶石呈束狀、部分呈半自形或它形粒狀結(jié)構(gòu); (c) 后形成的重晶石沿先形成的重晶石或方解石邊緣、解理或裂開等裂隙充填呈不規(guī)則網(wǎng)狀結(jié)構(gòu); (d) 重晶石被方解石不完全交代形成交代殘余結(jié)構(gòu)。礦物代號: Brt. 重晶石; Cal. 方解石。

    (a) 重晶石中的Ⅰ型純液體包裹體; (b) 方解石中的Ⅰ型純液體包裹體; (c)~(e) 重晶石中的負(fù)晶形、長條形、三角形、橢圓形和不規(guī)則形Ⅱ型液體包裹體; (f) 方解石中的負(fù)晶形Ⅱ型液體包裹體。

    表1 撒朗重晶石礦床流體包裹體顯微測溫結(jié)果

    流體包裹體冰點(diǎn)溫度為?20.8~?0.2 ℃(表1), 采用Bodnar (1993)的NaCl-H2O二元體系鹽度計算方程:

    =0.00+1.78?0.04422+0.0005573(為冰點(diǎn)攝氏度數(shù))

    計算其鹽度(表1)。結(jié)果顯示Ⅰ號礦體鹽度范圍0.5%~19.1% NaCleqv, 均值4.9%~8.4% NaCleqv; Ⅱ號礦體鹽度范圍1.4%~22.9% NaCleqv, 均值7.0%~9.1% NaCleqv; Ⅲ號礦體鹽度范圍1.2%~19.7% NaCleqv,均值8.7% NaCleqv。各礦體的鹽度范圍和均值也較為一致, 顯示相似的流體成礦作用特點(diǎn)。在鹽度分布直方圖上, 處于中低鹽度范圍, 具多峰態(tài)分布特點(diǎn)(圖6b), 峰值鹽度分別為2% NaCleqv、12% NaCleqv和16% NaCleqv, 反映流體成礦過程中鹽度的明顯變化。

    在均一溫度與鹽度的散點(diǎn)圖上, 溫度變化范圍相對較窄, 而鹽度的變化范圍相對較寬, 鹽度與溫度之間不存在明顯的線性關(guān)系(圖7)。

    圖6 撒朗重晶石礦床均一溫度(a)和鹽度分布(b)直方圖

    圖7 撒朗重晶石礦床流體鹽度與均一溫度關(guān)系圖

    3.1.3 流體包裹體成分

    由于流體包裹體的個體較小, 氣、液相組分更小, 因此沒能進(jìn)行單個包裹體成分的原位分析, 而采用了群體包裹體成分分析方法。撒朗重晶石礦床流體包裹體液、氣相重晶石群體包裹體成分分析結(jié)果見表2。

    從表2可以看出: 流體包裹體液相成分陽離子以Ca2+和Na+為主, 陰離子以Cl?和SO42?為主, 是一種Ca2+-Na+-Cl?-SO42?類型水。其中, Ca2+含量為28.77~42.57 μg/g, Na+含量為2.612~12.550 μg/g; Cl?含量為5.484~24.730 μg/g, SO42?含量為4.356~ 10.810 μg/g。

    流體包裹體氣相成分中, H2O是最主要的成分, 含量在99.95%以上, CO2、H2、N2、CH4和CO的含量都甚微, 是一種以水為主的流體。

    3.2 H-O-C同位素組成

    礦石礦物重晶石O同位素、流體包裹體水H同位素, 以及重晶石共生方解石的C、O同位素(PDB標(biāo)準(zhǔn))分析數(shù)據(jù)見表3。其中, 方解石δ18OSMOW值計算采用Coplen et al. (1983)的計算公式δ18OSMOW= 1.03091δ18OPDB+30.91。

    根據(jù)礦物流體包裹體均一法測溫所得到的各樣品或各礦體溫度數(shù)據(jù)以及重晶石、方解石與水之間的O同位素平衡分餾方程:

    1000ln重晶石?水=3.94×106?2?5.47×103?1+1.86 (鄭永飛和陳江峰, 2000)

    1000ln方解石?水=4.01×106?2?4.66×103?1+1.71 (鄭永飛和陳江峰, 2000)

    計算出與礦物O同位素交換達(dá)到平衡的成礦流體水的O同位素組成(表3)。

    從表3可以看出, 三個礦體重晶石O同位素δ18OSMOW值變化于15.7‰~18.4‰之間, 流體包裹體水H同位素δDSMOW為?48‰~?38‰, Ⅰ號和Ⅱ號兩個礦體中與重晶石共生方解石的O同位素δ18OPDB為?10.2‰~?9.6‰, 重晶石或方解石礦物的O同位素組成以及流體包裹體水H同位素組成變化范圍均小, 組成較均一, 反映其單一流體成礦過程。

    計算出的三個礦體與重晶石O同位素交換達(dá)到平衡的流體包裹體水的O同位素(SMOW)組成分別為Ⅰ號礦體: 2.35‰~12.05‰、均值8.23‰; Ⅱ號礦體: 4.42‰~5.08‰, 均值4.68‰; Ⅲ號礦體: 0.55‰~ 8.78‰, 均值4.68‰。與方解石O同位素交換達(dá)到平衡的流體包裹體水的O同位素組成(SMOW)分別為Ⅰ號礦體: 1.78‰~12.33‰, 均值8.14‰; Ⅱ號礦體: 5.38‰~6.10‰, 均值5.67‰。無論是重晶石還是方解石其流體包裹體水的O同位素組成變化范圍都較礦物本身的O同位素組成變化范圍要大得多, 是礦物與流體O同位素交換過程中受流體溫度的變化影響所致, 表明溫度是控制流體O同位素變化的重要因素。

    表2 撒朗重晶石礦床流體包裹體液、氣相成分分析數(shù)據(jù)表

    表3 撒朗重晶石礦床重晶石H-O和方解石C-O同位素組成

    Ⅰ號和Ⅱ號兩個礦體方解石C同位素δ13CPDB組成為?1.9‰~?1.0‰, 平均?1.45‰, 變化較小, 顯示碳具相同來源。

    4 討 論

    4.1 成礦流體特征

    礦區(qū)三個礦體的均一溫度都較為接近, 變化范圍(95~240 ℃)稍大(表1), 均值(123~166 ℃)則相對集中??傮w呈單峰態(tài)分布, 位于中低溫度范圍, 峰值溫度為140 ℃, 屬于低溫區(qū)域(圖6a), 與多數(shù)重晶石礦床如陜西神河鋇礦床(吳勝華等, 2011)、伊朗Farsesh重晶石礦床(Asl et al., 2015)和Chenarvardeh重晶石礦床(Ehya and Mazraei, 2017)以及墨西哥Múzquiz重晶石礦床(González-Sánchez et al., 2017)等的形成溫度相似。各礦體的均一溫度范圍和均值也都較為接近, 應(yīng)是統(tǒng)一的熱液流體成礦活動所形成。

    三個礦體的成礦流體鹽度也很相似, 變化范圍(0.5%~22.9% NaCleqv)相對較大(表1), 系中低鹽度范圍。而均值(4.9%~9.1% NaCleqv)卻較集中, 具低鹽度特點(diǎn), 與陜西Ankang-Xunyang重晶石?毒重石鋇礦床(Wu et al., 2015)和土耳其Adana-Feke地區(qū)重晶石礦床(Ozdogan et al., 2017)等的成礦流體鹽度特點(diǎn)相似。各礦體的鹽度值范圍和均值較為一致, 顯示其相似的流體成礦作用特點(diǎn)。但鹽度具多峰分布特點(diǎn)(圖6b), 反映成礦過程中流體鹽度的明顯變化。而從均一溫度與鹽度不具明顯的線性關(guān)系(圖7)來看, 鹽度的變化可能并不僅僅是溫度變化引起, 而更可能是由于溫度降低導(dǎo)致礦石礦物等沉淀, 從而使得成礦流體中礦質(zhì)濃度的變化所致。

    根據(jù)流體包裹體的均一溫度和鹽度數(shù)據(jù), 可以求得成礦流體的密度, 其方法很多, 有公式計算法(劉斌和沈崑, 1999; Steele-MacInnis et al., 2012)、圖解法(Wilkinson, 2001; 盧煥章等, 2004)和查數(shù)法(劉斌和沈崑, 1999), 不同方法的思路總體相近, 本文采用是圖解法。利用測溫所獲得的均一溫度和鹽度數(shù)據(jù), 在Wilkinson (2001)溫度與鹽度圖上投點(diǎn)獲得流體密度(李應(yīng)輝, 2016)。結(jié)果顯示, 重晶石中流體包裹體密度為0.81~1.08 g/cm3、均值0.95 g/cm3, 方解石中流體包裹體密度為0.92~1.04 g/cm3、均值0.97 g/cm3??梢? 成礦流體均屬中等密度流體。

    4.2 成礦流體來源

    熱液流體是礦床尤其是熱液礦床形成時最主要的含礦介質(zhì), 在礦質(zhì)的活化、遷移和沉淀過程中起著極重要的作用(何明勤等, 2004a)。同時, 不同熱液成礦流體水的來源不同(White, 1974), 且具體礦床應(yīng)有其主要來源。通常, 成礦流體水的來源主要有海水、雨水和初生水, 以及由它們再循環(huán)或混合而形成的建造水、巖漿水和變質(zhì)水(Hoefs, 2009), 由于它們具明顯不同的同位素組成特征, 因此常用H、O同位素來推斷熱液流體來源及路徑(Dilles et al., 1992)。

    重晶石礦物為非含氫礦物, 不存在礦物與流體間的H同位素交換, 因此重晶石的流體包裹體水H同位素組成可作為成礦流體的H同位素組成。同時, 因與重晶石共生的方解石也是非含氫礦物, 因而形成重晶石的成礦流體H同位素組成同樣也可視作形成共生方解石的成礦流體H同位素組成。

    根據(jù)重晶石礦物O同位素值計算獲得的成礦流體O同位素平均值、最小值與最大值, 所反映的O同位素變化范圍及其對應(yīng)的H同位素值進(jìn)行δD-δ18O成圖, 結(jié)果顯示樣品H、O同位素組成及其變化除落在巖漿水或變質(zhì)水范圍內(nèi), 還更多地位于大氣降水線與變質(zhì)水之間區(qū)域(圖8)。介于大氣降水線與變質(zhì)水、甚至巖漿水之間區(qū)域的成礦流體水, 不同研究者曾有不同認(rèn)識, 有的稱為地層水(涂光熾等, 1988; Taylor, 1997), 有的稱為建造水(冉崇英, 1989), 有的稱為原生水(陳好壽和冉崇英等, 1992), 有的稱為改造水(何明勤等, 2004c), 并圈出了δD-δ18O圖上地層水的分布范圍(G?k?e and Bozkaya, 2008; Melo et al., 2019), 但因圈出的地層水范圍過大反而導(dǎo)致其含義不清而失去意義。一般地, 這種成礦流體水的來源多樣且復(fù)雜, 難以準(zhǔn)確確定, 需結(jié)合礦區(qū)具體特點(diǎn)和其區(qū)域成礦地質(zhì)背景研究后進(jìn)行綜合判斷。撒朗重晶石礦區(qū)乃至區(qū)域上缺乏巖漿巖出露或巖漿活動, 也沒有區(qū)域變質(zhì)巖或區(qū)域變質(zhì)作用。因此,成礦流體水不應(yīng)主要來源于巖漿水或變質(zhì)水, 最可能的是來源于大氣降水或演化后的大氣降水, 是下滲大氣降水被加熱并與O同位素值更大的圍巖水巖反應(yīng)發(fā)生O同位素交換的結(jié)果。且溫度不同, 其水巖反應(yīng)O同位素交換程度也不同, 從而造成流體O同位素組成也不一樣, 明顯偏離雨水線, 呈現(xiàn)δ18O的顯著右漂移現(xiàn)象, 這與其他一些熱液脈型重晶石礦床如廣西泥盆紀(jì)鉛?鋅重晶石礦床的成礦流體水最可能來自于雨水(Chen and Gao, 1988)以及一些重晶石?方鉛礦礦床如土耳其Karralar重晶石?方鉛礦礦床的成礦流體水也來源于雨水(G?k?e and Bozkaya, 2008)類似。

    δD(‰)=8δ18O+10為Craig (1961)雨水線方程; SMOW(或海水)與變質(zhì)水范圍據(jù)Taylor (1974); 原始巖漿水范圍據(jù)Campbell et al. (1984)。

    4.3 礦床成因初步分析

    撒朗重晶石礦床原生礦體主要呈層狀、似層狀、脈狀、囊狀產(chǎn)于下奧陶統(tǒng)桐梓組(O1)白云巖等碳酸鹽類巖石的層間破碎帶及其構(gòu)造交匯部位, 受區(qū)域褶皺、圍巖地層、層間斷層和巖性等因素的控制(李應(yīng)輝等, 2016)。流體包裹體測溫數(shù)據(jù)顯示, 三個礦體的均一溫度為95~240 ℃, 明顯高于正常沉積礦床的形成溫度。同時, 礦床也不具備一般熱水沉積礦床所具有的硅質(zhì)巖和電氣石巖等熱水沉積巖, 層狀和似層狀礦上下盤的明顯不對稱蝕變, 以及上部為層狀、似層狀礦而下部為細(xì)脈狀、筒狀礦的“雙層構(gòu)造”等特征(姚鳳良和孫豐月, 2006)。

    從熱液方解石C同位素組成(表3)可知, 方解石C同位素δ13CPDB組成為?1.9‰~?1.0‰, 平均?1.45‰,接近于零值。對于熱液碳酸鹽礦物來說, 其碳主要有三大來源: ①巖漿源或深部來源, δ13C均值約為?7‰; ②沉積碳酸鹽巖來源, δ13C約為0; ③地殼巖石中的有機(jī)碳來源, δ13C均值約為?25‰(何明勤等, 2004a)。撒朗重晶石礦床含碳熱液礦物僅由方解石組成或由方解石與重晶石組成的熱液體系, 其熱液方解石C同位素組成基本反映了熱液體系C同位素組成(Ohomoto, 1972)。因此, 其礦石熱液方解石接近于零值的C同位素均值(?1.45‰)反映撒朗重晶石礦床熱液成礦流體的碳主要來源于沉積碳酸鹽巖。結(jié)合熱液方解石的C同位素以及相應(yīng)的礦物O同位素特征(表3), 以及δ18OSMOW-δ13CPDB圖解(劉建明和劉家軍, 1997; 何明勤等, 2004a; 畢獻(xiàn)武等, 2019), 可以更加確定地反映熱液方解石碳的來源。且樣品點(diǎn)靠近海相碳酸鹽巖區(qū)域, 并沿碳酸鹽巖溶解作用方向分布, 進(jìn)一步表明碳來源于海相碳酸鹽巖的溶解作用, 這與礦床產(chǎn)出的圍巖地層為白云巖等海相碳酸鹽巖的實際地質(zhì)情況相吻合。

    碳來源于海相碳酸鹽巖的溶解以及重晶石H-O同位素研究顯示成礦流體水最可能來源于大氣降水或演化后的大氣降水, 而不是同生海水或深部流體水。這些特征表明, 礦床并不是同生形成的, 而是明顯晚于圍巖地層, 是圍巖形成之后, 受后期構(gòu)造作用影響, 地層圍巖發(fā)生變形形成包括層間斷層等在內(nèi)的褶皺與斷裂構(gòu)造, 為礦體的賦存提供了條件。與此同時, 大氣降水沿斷裂構(gòu)造下滲, 并逐漸被加熱形成熱液, 從而具有了更強(qiáng)的溶解能力, 并在白云巖等碳酸鹽類巖石和層間斷層等構(gòu)造的循環(huán)中, 不斷溶解巖石地層中的可溶性鹽類和成礦物質(zhì)而成為含礦鹽水熱液, 進(jìn)一步提高了其溶解萃取能力。

    富含可溶性鹽類和成礦物質(zhì)的含礦鹽水熱液, 均一溫度為95~240 ℃、峰值溫度140 ℃, 在沿地層巖石和斷裂構(gòu)造等的滲流循環(huán)過程中, 隨著溫度的降低以及與圍巖的水巖反應(yīng), 其成礦物質(zhì)濃度逐漸飽和沉淀, 并在層間破碎帶等斷裂構(gòu)造及其交匯部位富集形成層狀、似層狀、脈狀、囊狀重晶石礦體。據(jù)中國礦床編委會(1994), 貴州奧陶系碳酸鹽巖中的脈狀重晶石礦床其重晶石的S同位素δ34S值特征(δ34S=35.35‰), 與含礦地層同時期海水硫酸鹽S同位素平均值相近。從而認(rèn)為, 早奧陶世白云巖沉積時期, 就已形成礦源本底地層。并推測, 后期的燕山運(yùn)動構(gòu)造作用使賦礦地層中Ba等有關(guān)組分活化遷移并就近集中在構(gòu)造裂隙中富集形成重晶石礦床(中國礦床編委會, 1994), 屬中低溫?zé)嵋旱V床, 但主體(峰值溫度140 ℃)形成于200 ℃以下的低溫階段。綜上所述, 都勻撒朗重晶石礦床具體成因模式如圖9所示。

    5 結(jié) 論

    通過對都勻撒朗重晶石礦床流體包裹體和H-O-C同位素研究, 得到如下主要結(jié)論:

    (1) 流體包裹體均一溫度95~240 ℃, 峰值溫度140 ℃; 成礦流體鹽度0.5%~22.9% NaCleqv, 均值4.9%~9.1% NaCleqv; 流體密度0.81~1.08 g/cm3, 均值0.95~0.97 g/cm3; 流體包裹體液相成分陽離子以Ca2+和Na+為主, 陰離子以Cl?和SO42?為主。成礦流體是一種具有中低溫、峰值溫度為低溫、低鹽度、中等密度特征的Ca2+-Na+-Cl?-SO42?類型水。

    (2) 重晶石流體包裹體水δDSMOW為?48‰~ ?38‰, 成礦流體水的δ18OSMOW為0.55‰~12.05‰、均值4.86~8.23‰, 成礦流體水來源于大氣降水或是與圍巖在不同溫度條件下水巖反應(yīng)后的演化大氣降水。

    (3) 礦石中方解石的C同位素δ13CPDB為?1.9‰~?1.0‰, 接近于零值, 其相應(yīng)方解石的δ18OSMOW為20.40‰~21.01‰, 碳來自于圍巖海相碳酸鹽巖的溶解。

    (4) 礦體主要呈層狀、似層狀、脈狀、囊狀產(chǎn)于下奧陶統(tǒng)桐梓組(O1)白云巖等碳酸鹽類巖石的層間破碎帶及其構(gòu)造交匯部位, 其形成晚于圍巖地層, 是后期的構(gòu)造作用及其熱液活動使其中的成礦組分活化遷移并就近集中在相應(yīng)構(gòu)造裂隙中富集而形成, 屬中低溫?zé)嵋旱V床, 且主要形成于200 ℃以下的低溫階段。

    圖9 都勻撒朗重晶石礦床成因模式簡圖

    致謝:野外礦床地質(zhì)調(diào)查和樣品采集過程中得到貴州省黔南州三山礦業(yè)有限責(zé)任公司的配合與幫助; 兩位匿名評審專家對論文進(jìn)行了認(rèn)真細(xì)致審閱, 提出了建設(shè)性修改意見, 謹(jǐn)此一并表示感謝!

    畢獻(xiàn)武, 唐永永, 陶琰, 王長明, 胥磊落, 戚華文, 蘭青, 木蘭. 2019. 西南三江碰撞造山帶沉積巖容礦Pb-Zn- Ag-Cu賤金屬復(fù)合成礦與深部過程. 巖石學(xué)報, 35(5): 1341–1371.

    陳國勇, 鐘奕天, 黃根深. 1992. 都勻牛角塘鋅礦床地質(zhì)特征及成礦控制條件初探. 貴州地質(zhì), 9(3): 203–212.

    陳好壽, 冉崇英. 1992. 康滇地軸銅礦床同位素地球化學(xué). 北京: 地質(zhì)出版社: 68–74.

    方維萱, 胡瑞忠, 蘇文超, 漆亮, 肖加飛, 蔣國豪. 2002. 大河邊?新晃超大型重晶石礦床地球化學(xué)特征及形成的地質(zhì)背景. 巖石學(xué)報, 18(2): 247–256.

    高軍波, 楊瑞東. 2015. 水城紫云?南丹深大斷裂構(gòu)造演化與泥盆紀(jì)熱水(液)成巖、成礦效應(yīng). 貴陽: 貴州科技出版社: 102–134.

    高軍波, 楊瑞東, 陶平, 程偉, 何志威, 劉坤. 2013a. 貴州鎮(zhèn)寧泥盆系大型重晶石礦床鍶同位素組成特征研究. 地球化學(xué), 42(4): 385–392.

    高軍波, 楊瑞東, 陶平, 魏懷瑞, 劉坤, 程偉. 2013b. 貴州鎮(zhèn)寧泥盆系大型重晶石礦床地球化學(xué)特征及其成因研究. 現(xiàn)代地質(zhì), 27(1): 46–55.

    貴州省地質(zhì)調(diào)查院. 2012. 貴州省區(qū)域地質(zhì)志. 貴陽: 貴州省地質(zhì)調(diào)查院: 1035–1114.

    貴州省區(qū)域礦產(chǎn)志編寫組. 1986. 貴州省區(qū)域礦產(chǎn)志. 貴陽: 貴州省地質(zhì)礦產(chǎn)局區(qū)域地質(zhì)調(diào)查大隊: 419–435.

    何明勤, 劉家軍, 李朝陽, 李志明, 劉玉平. 2004a. 蘭坪盆地鉛鋅銅大型礦集區(qū)的流體成礦作用機(jī)制——以白秧坪銅鈷多金屬礦區(qū)為例. 北京: 地質(zhì)出版社: 43–103.

    何明勤, 劉家軍, 楊世瑜, 李朝陽. 2004b. 云南賓川雄魯摩銅多金屬礦床的成礦流體. 礦物學(xué)報, 24(3): 261– 265.

    何明勤, 楊世瑜, 陳昌勇, 馬德云, 鐘昆明. 2004c. 滇西小龍?zhí)?馬廠箐地區(qū)銅金多金屬礦床地質(zhì)地球化學(xué)及成因研究. 北京: 地質(zhì)出版社: 39–85.

    侯東壯, 吳湘濱, 李貞, 劉玉紅. 2015. 貴州省天柱大河邊重晶石礦床成礦物質(zhì)來源. 中國有色金屬學(xué)報, 25(4): 1039–1048.

    李文炎, 余洪云. 1991. 中國重晶石礦床. 北京: 地質(zhì)出版社: 1–105.

    李應(yīng)輝. 2016. 貴州都勻撒朗重晶石礦床成礦流體研究. 貴陽: 貴州大學(xué)碩士學(xué)位論文: 1–52.

    李應(yīng)輝, 何明勤, 余潤濤. 2016. 貴州都勻撤朗重晶石礦床地質(zhì)特征和控礦因素分析. 貴州大學(xué)學(xué)報(自然科學(xué)版), 33(2): 57–61, 66.

    劉斌, 沈崑. 1999. 流體包裹體熱力學(xué). 北京: 地質(zhì)出版社: 73–118.

    劉成萬. 2017. 貴州都勻大梁子鋅礦床地質(zhì)特征及成因淺析. 西部探礦工程, 29(1): 85–88.

    劉建明, 劉家軍. 1997. 滇黔桂金三角區(qū)微細(xì)浸染型金礦床的盆地流體成因模式. 礦物學(xué)報, 17(4): 448–456.

    劉強(qiáng), 夏磊. 2013. 貴州都勻老貓石鉛鋅礦地質(zhì)特征. 礦產(chǎn)與地質(zhì), 27(5): 393–395.

    劉鐵庚, 葉霖. 2000. 都勻牛角塘大型獨(dú)立鎘礦床的地質(zhì)地球化學(xué)特征. 礦物學(xué)報, 20(3): 279–285.

    劉鐵庚, 張乾, 葉霖, 邵樹勛. 2004. 貴州牛角塘鎘鋅礦床中發(fā)現(xiàn)原生硫鎘礦. 礦物學(xué)報, 24(2): 191–196.

    盧煥章, 范宏瑞, 倪培, 歐光習(xí), 沈昆, 張文淮. 2004. 流體包裹體. 北京: 地質(zhì)出版社: 201–229.

    冉崇英. 1989. 康滇地軸銅礦床的成礦機(jī)理. 北京: 地質(zhì)出版社: 33–35.

    涂光熾等. 1988. 中國層控礦床地球化學(xué)(第三卷). 北京: 科學(xué)出版社: 1–26.

    王瓊, 息朝莊, 王國坤, 程俊義. 2017. 貴州都勻菜園河鋅礦床成礦規(guī)律及找礦標(biāo)志. 現(xiàn)代礦業(yè), (6): 6–12.

    吳朝東, 楊承運(yùn), 陳其英. 1999. 新晃貢溪天柱大河邊重品石礦床熱水沉積成因探討. 北京大學(xué)學(xué)報, 35(6): 74–785.

    吳勝華, 劉家軍, 柳震江, 翟德高, 邢永亮. 2011. 陜西神河鋇礦床礦物流體包裹體特征與成礦物理化學(xué)條件. 礦床地質(zhì), 30(6): 1100–1112.

    夏菲, 馬東升, 潘家永, 孫占學(xué), 曹雙林, 陳少華, 聶文明, 吳凱, 劉莉. 2005. 天柱大河邊重晶石礦床鉛同位素特征及來源探討. 地球化學(xué), 34(5): 501–507.

    夏菲, 馬東升, 潘家永, 孫占學(xué), 曹雙林, 聶文明, 吳凱. 2004. 貴州天柱大河邊和玉屏重晶石礦床熱水沉積成因的鍶同位素證據(jù). 科學(xué)通報, 49(24): 2592–2595.

    楊瑞東, 鮑淼, 魏懷瑞, 王偉, 王強(qiáng). 2007a. 貴州天柱寒武系底部重晶石礦床中熱水生物群的發(fā)現(xiàn)及意義. 自然科學(xué)進(jìn)展, 17(9): 304–1309.

    楊瑞東, 魏懷瑞, 鮑淼, 王偉, 王強(qiáng). 2007b. 貴州天柱上公塘–大河邊寒武紀(jì)重晶石礦床海底熱水噴流沉積結(jié)構(gòu)、構(gòu)造特征. 地質(zhì)論評, 53(5): 675–680.

    姚鳳良, 孫豐月. 2006. 礦床學(xué)教程. 北京: 地質(zhì)出版社: 132–153.

    葉霖, 劉鐵庚. 1997. 都勻地區(qū)鎘(Cd)礦資源及其遠(yuǎn)景初探. 貴州地質(zhì), 14(2): 160–163.

    張碧志, 草明, 李明道. 1994 都勻牛角塘鋅礦田控礦條件及成礦模式探討. 貴州地質(zhì), 11(4): 87–293.

    張旭, 高軍波, 楊瑞東, 陳軍, 魏懷瑞, 徐世林, 鮑淼. 2017. 黔南泥盆系重晶石礦床地質(zhì)特征與成礦規(guī)律. 地球化學(xué), 46(5): 456–465.

    鄭永飛, 陳江峰. 2000. 穩(wěn)定同位素地球化學(xué). 北京: 科學(xué)出版社: 38–73.

    中國礦床編委會. 1994. 中國礦床(下冊). 北京: 地質(zhì)出版社: 281–314.

    Asl S M, Jafari M, Sahamiyeh R Z and Shahrokhi V. 2015. Geology, geochemistry, sulfur isotope composition, and fluid inclusion data of Farsesh barite deposit, Lorestan Province, Iran., 8: 7125– 7139.

    Bodnar R J. 1993. Revised equation and table for determining the freezing point depression of H2O-NaCl solutions., 57(3): 683–684.

    Campbell A, Rye D and Petersen U. 1984. A hydrogen and oxygen isotope study of the San Cristobal mine, Peru: Implications of the role of water to rock ratio for the genesis of wolframite deposits., 79(8): 1818–1832.

    Chen X P and Gao J Y. 1988. Thermal water deposition and Pb-Zn barite deposits in the Devonian System, Central Guangxi., 7(4): 321– 328.

    Coplen T B, Kendall C and Hopple J. 1983. Comparison of stable isotope reference samples., 302: 236–238.

    Craig H. 1961. Isotopic variation in meteoric waters., 133: 1702–1703.

    Dilles J H, Solomon G C, Taylor H P and Einaudi M T. 1992. Oxygen and hydrogen isotope characteristics of hydro-thermal alteration at the Ann-Mason porphyry copper deposit, Yerington, Nevada., 87(1): 44–63.

    Ehya F and Mazraei S M. 2017. Hydrothermal barite mineralization at Chenarvardeh deposit, Markazi Province, Iran: Evidences from REE geochemistry and fluid inclusions., 134: 299–307.

    Elmas N, Kumral M, Suner F and Tasdenen S. 2012. Strati-form barite deposits hosted in metamorphic assemblages of Dinek and surrounding regions, Isparta, Turkey., 48: 150–159.

    Gao J B, Yang R D, Chen J, Zheng L L, Cheng W and Wei H R. 2017. Multiple proxies indicating methane seepage as the origin of Devonian large barite deposit in Zhenning-Ziyun, Guizhou, SW China., 80: 18–26.

    G?k?e A and Bozkaya G. 2008. Fluid inclusion and stable isotope characteristics of the Karralar (Gazipasa, Antalya) barite-galena deposits, Southern Turkey., 50(2): 145–154.

    González-Sánchez F, González-Partida E, Canet C, Atudorei V, Atudorei P, Morales-Puente P, Cienfuegos-Alvarado E and González-Ruiz L. 2017. Geological setting and genesis of stratabound barite deposits at Múzquiz, Coahuila in Northeastern Mexico., 81: 1184–1192.

    Han S C, Hu K, Gao J, Pan J Y, Xia F and Wu W F. 2015. Origin of early Cambrian black-shale-hosted barite deposits in South China: Mineralogical and geochemical studies., 106: 79–94.

    Hoefs J. 2009. Stable Isotope Geochemistry. Berlin: Springer- Verlag Berlin and Heidelberg (Sixth edition): 93–227.

    Mapani B E S and Wilson C J L. 1998. Evidence for externallyderived vein forming and mineralizing fluids: An example from the Magdala gold mine, Victoria, Australia., 13(1–5): 323–343.

    Maynard J B and Okita P M. 1991. Bedded barite deposits in the United States, Canada, Germany, and China: Two major types based on tectonic setting., 86(2): 364–376.

    Melo G H C, Monteiro L V S, Xavier R P, Moreto C P N and Santiago E. 2019. Tracing fluid sources for the Salobo and Igarapé Bahia deposits: Implications for the genesis of the iron oxide copper-gold deposits in the Carajás Province, Brazil., 114(4): 697–718.

    Ohomoto H. 1972. Systematics of sulfur and carbon isotopes in hydrothermal ore deposits., 67: 551–578.

    Ozdogan A T, Uras Y and Oner F. 2017. Geochemistry of the barite deposits near Adana-Feke area (Eastern Taurides)., 58(11): 1349–1365.

    Steele-MacInnis M, Lecumberri-Sanchez P and Bodnar R J. 2012. HOKIEFLINCS_H2O-NACL: A Microsoft Excel spreadsheet for interpreting microthermometric data from fluid inclusions based on the PVTX properties of H2O-NaCl., 49: 334–37.

    Taylor H P Jr. 1974. The application of oxygen and hydrogen isotope studies to problems of hydrothermal alteration and ore deposition., 69(6): 843–883.

    Taylor H P Jr. 1997. Oxygen and hydrogen isotope relationships in hydrothermal mineral deposits // Barnes H L. Geochemistry of Hydrothermal Ore Deposits. Hoboken: John Wiley & Sons, Inc.: 229–302.

    White D E. 1974. Diverse origins of hydrothermal ore fluids., 69(6): 954–973.

    Wilkinson J J. 2001. Fluid inclusions in hydrothermal ore deposits., 55: 229–272.

    Wu S H, Liu J H and Zhai D G. 2015. Factors controlling precipitation of barite and witherite and genesis of Ankang-Xunyang barium deposits, Shanxi, China.(), 89(3): 836–851.

    Zhou Z B, Wen H J, Qin C J, Fourestier J, Liu L and Shi Q P. 2018. The genesis of the Dahebian Zn-Pb deposit and associated barite mineralization: Implications for hydrothermal fluid venting events along the Nanhua Basin, South China., 101: 785–802.

    Fluid Inclusion and H-O-C Isotope Geochemistry of the Salang Barite Deposit in Duyun, Guizhou Province, China

    HE Mingqin1, 2*, LI Yinghui3, ZHANG Feng1, TIAN Huanhuan1and LIANG Kunping1

    (1. College of Resources and Environmental Engineering, Guizhou University, Guiyang 550025, Guizhou, China; 2. MOE Key Laboratory of Karst Geological Resources and Environments, Guizhou University, Guiyang 550025, Guizhou, China; 3. 105 Geological Unit of the Bureau of Geology and Mineral Resources of Guizhou Province, Guiyang 550018, Guizhou, China)

    The Salang barite deposit is located in the Duyun north-south tectonic deformation area of the Duyun continental margin basin in the SW section of the Jiangnan orogenic belt. The ore bodies occur in the interlayer fracture zone or at the junction of structures, showing layer, layer-like, vein or capsule shapes. Based on the field geological investigation, samples were collected for laboratory study. The petrography, microthermometry, and composition of the fluid inclusions, H-O-C isotope compositions of barite and calcite of the main metallogenic epoch of the deposit were studied. The results showed that the ore-forming fluid of the Salang barite deposit is of Ca2+-Na+-Cl?-SO42?type which has the characteristics of low to medium temperature (95~240 ℃, with a peak temperature of 140 ℃), low salinity (averages=4.9%?9.1% NaCleqv), and medium density (averages=0.95?0.97 g/cm3). The fluid inclusions in barite have δD values of ?48‰ to ?38‰ and δ18O of 0.55‰ to 12.33‰, which indicate that the ore-forming fluid was originated from meteoric water or evolved meteoric water because of water-rock reaction. Carbon and oxygen isotopic compositions of calcite separated from the ores vary from ?1.9‰ to ?1.0‰ and 20.40‰ to 21.01‰ (relative to SMOW) respectively, which reflect that the carbon was likely derived from the dissolution of marine carbonate rocks. It can be inferred that the Salang barite deposit is of epigenetic, formed during the late tectonic and hydrothermal activity while ore-forming materials were derived from the surrounding rocks.

    fluid inclusions; H-O-C isotopes; ore deposit geochemistry; Salang barite deposit; Duyun of Guizhou province, China

    2019-12-20;

    2020-07-30

    國家自然科學(xué)地區(qū)基金項目(41563003)和聯(lián)合基金項目(U1812402)共同資助。

    何明勤(1966?), 男, 博士, 教授, 主要從事礦床學(xué)及礦床地球化學(xué)方面的教學(xué)與科研工作。Email: mingqinhe@126.com

    P611; P597

    A

    1001-1552(2021)03-0509-014

    10.16539/j.ddgzyckx.2021.03.004

    猜你喜歡
    成礦
    歡迎訂購2022年《大地構(gòu)造與成礦學(xué)》
    歡迎訂購2022年《大地構(gòu)造與成礦學(xué)》
    淺析淶源縣大黃峪鉬礦成礦遠(yuǎn)景
    桂西沉積型鋁土礦床成礦規(guī)律及成礦模式
    歡迎訂購2021年《大地構(gòu)造與成礦學(xué)》
    歡迎訂購2021年《大地構(gòu)造與成礦學(xué)》
    新疆寨北山銅礦成礦新認(rèn)識及找礦預(yù)測
    氧化還原作用對鈾成礦的意義——以相山西部流紋英安巖鈾礦為例
    中非銅鈷成礦帶礦床遙感特征研究
    柴達(dá)木盆地北緣鋰多金屬礦成礦條件及找礦潛力
    91字幕亚洲| 老司机在亚洲福利影院| 国产免费一级a男人的天堂| 日本 av在线| 日本免费a在线| av女优亚洲男人天堂| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 色综合婷婷激情| 亚洲一区二区三区色噜噜| 中文字幕精品亚洲无线码一区| 久久精品国产99精品国产亚洲性色| 少妇的逼水好多| 成年女人看的毛片在线观看| 国产精品一区二区三区四区久久| 禁无遮挡网站| 亚洲欧美日韩高清在线视频| 精品国产超薄肉色丝袜足j| 精品国产三级普通话版| 日本a在线网址| 51国产日韩欧美| 精品久久久久久成人av| 亚洲真实伦在线观看| 亚洲精品在线美女| 99久久99久久久精品蜜桃| 观看免费一级毛片| 两个人看的免费小视频| 日韩大尺度精品在线看网址| 成人一区二区视频在线观看| 国产伦一二天堂av在线观看| h日本视频在线播放| 精品国产超薄肉色丝袜足j| 97碰自拍视频| 色噜噜av男人的天堂激情| 在线观看一区二区三区| 精品乱码久久久久久99久播| 国产91精品成人一区二区三区| 3wmmmm亚洲av在线观看| 亚洲一区二区三区色噜噜| 久久欧美精品欧美久久欧美| av国产免费在线观看| 丰满人妻一区二区三区视频av | 露出奶头的视频| 日本a在线网址| 毛片女人毛片| 国产亚洲精品久久久com| 国产成人av教育| 国产精品久久电影中文字幕| 九色国产91popny在线| 成人欧美大片| netflix在线观看网站| 久久精品亚洲精品国产色婷小说| 国产成人a区在线观看| 精品久久久久久久久久久久久| 国产精品野战在线观看| 久久久久国产精品人妻aⅴ院| 日韩中文字幕欧美一区二区| 精品人妻1区二区| 免费av毛片视频| 丰满乱子伦码专区| 色精品久久人妻99蜜桃| 国产成人av教育| 99久久成人亚洲精品观看| 搞女人的毛片| 亚洲 欧美 日韩 在线 免费| 精品久久久久久,| 丰满人妻一区二区三区视频av | 性欧美人与动物交配| 国产麻豆成人av免费视频| 久久精品91无色码中文字幕| 一区二区三区免费毛片| 日韩有码中文字幕| 精品久久久久久成人av| 麻豆一二三区av精品| 又紧又爽又黄一区二区| 欧美日本亚洲视频在线播放| 丰满的人妻完整版| 高清日韩中文字幕在线| 久久中文看片网| 欧美一区二区亚洲| 性色avwww在线观看| av片东京热男人的天堂| 黑人欧美特级aaaaaa片| 久久久久久久亚洲中文字幕 | 97超级碰碰碰精品色视频在线观看| 午夜激情福利司机影院| 色综合欧美亚洲国产小说| 啦啦啦观看免费观看视频高清| 熟妇人妻久久中文字幕3abv| 精品一区二区三区视频在线观看免费| 在线视频色国产色| 99热只有精品国产| 日日干狠狠操夜夜爽| 久久午夜亚洲精品久久| 亚洲欧美精品综合久久99| 免费观看精品视频网站| 成年版毛片免费区| 国产伦在线观看视频一区| 成人一区二区视频在线观看| 国产在线精品亚洲第一网站| 欧美乱色亚洲激情| 欧美日韩黄片免| 精品久久久久久,| 天堂√8在线中文| 欧美成人a在线观看| 两人在一起打扑克的视频| 欧美大码av| 精品福利观看| 欧美日本亚洲视频在线播放| 国产av不卡久久| 欧美xxxx黑人xx丫x性爽| 日本黄色片子视频| 脱女人内裤的视频| 精品国产亚洲在线| 丁香欧美五月| 精品国产美女av久久久久小说| 亚洲av熟女| 九九久久精品国产亚洲av麻豆| 精品久久久久久久毛片微露脸| 99riav亚洲国产免费| 禁无遮挡网站| 国产精品美女特级片免费视频播放器| tocl精华| 国产精品美女特级片免费视频播放器| 色综合欧美亚洲国产小说| 99久久99久久久精品蜜桃| 亚洲一区二区三区色噜噜| 国产主播在线观看一区二区| 我要搜黄色片| 99riav亚洲国产免费| 国内精品久久久久久久电影| 成年女人毛片免费观看观看9| 久久久久九九精品影院| 欧美乱妇无乱码| 国产精品亚洲av一区麻豆| 免费一级毛片在线播放高清视频| 在线观看66精品国产| 搡老岳熟女国产| 19禁男女啪啪无遮挡网站| aaaaa片日本免费| 成人高潮视频无遮挡免费网站| 大型黄色视频在线免费观看| 人人妻,人人澡人人爽秒播| 国产精品av视频在线免费观看| 久久精品国产清高在天天线| 99热6这里只有精品| 成人永久免费在线观看视频| 91av网一区二区| 中国美女看黄片| 精品国内亚洲2022精品成人| av国产免费在线观看| a在线观看视频网站| 国内精品久久久久精免费| 亚洲国产精品合色在线| 熟妇人妻久久中文字幕3abv| 日韩免费av在线播放| 国产高清激情床上av| av中文乱码字幕在线| 亚洲精品国产精品久久久不卡| 中文字幕熟女人妻在线| 久久九九热精品免费| 国内久久婷婷六月综合欲色啪| 少妇熟女aⅴ在线视频| 村上凉子中文字幕在线| 亚洲狠狠婷婷综合久久图片| 嫁个100分男人电影在线观看| bbb黄色大片| 制服人妻中文乱码| 搡老熟女国产l中国老女人| 精品人妻偷拍中文字幕| 国产成年人精品一区二区| 日韩欧美精品免费久久 | 欧美3d第一页| 每晚都被弄得嗷嗷叫到高潮| www.色视频.com| 国产一区二区激情短视频| 母亲3免费完整高清在线观看| 欧美在线一区亚洲| 一个人观看的视频www高清免费观看| 精品一区二区三区av网在线观看| 日韩成人在线观看一区二区三区| 午夜精品久久久久久毛片777| 中文字幕av在线有码专区| ponron亚洲| 日本与韩国留学比较| 国产黄a三级三级三级人| 最后的刺客免费高清国语| 色尼玛亚洲综合影院| 又黄又爽又免费观看的视频| 色在线成人网| 一个人免费在线观看的高清视频| 国产成人av教育| 又黄又粗又硬又大视频| 亚洲18禁久久av| 黄片大片在线免费观看| 国产成人欧美在线观看| 国产成人av激情在线播放| 欧美日韩瑟瑟在线播放| 亚洲精品日韩av片在线观看 | 欧美日韩福利视频一区二区| 99热精品在线国产| 悠悠久久av| 国产一区二区激情短视频| 久久精品影院6| 亚洲专区中文字幕在线| 午夜a级毛片| 听说在线观看完整版免费高清| 少妇裸体淫交视频免费看高清| a在线观看视频网站| 天堂网av新在线| 身体一侧抽搐| 色综合亚洲欧美另类图片| 天天躁日日操中文字幕| 亚洲国产精品久久男人天堂| 手机成人av网站| 一区福利在线观看| 18禁裸乳无遮挡免费网站照片| 国产激情偷乱视频一区二区| 岛国视频午夜一区免费看| 国产精品亚洲av一区麻豆| 国产极品精品免费视频能看的| 免费看美女性在线毛片视频| 真实男女啪啪啪动态图| 国产高清激情床上av| 51国产日韩欧美| 欧美日韩综合久久久久久 | 小说图片视频综合网站| 欧美激情在线99| 看片在线看免费视频| 国产主播在线观看一区二区| 国产单亲对白刺激| 精品人妻1区二区| 久久性视频一级片| 又紧又爽又黄一区二区| 欧美在线一区亚洲| 内射极品少妇av片p| 亚洲熟妇熟女久久| www.999成人在线观看| 美女高潮的动态| 三级国产精品欧美在线观看| 色精品久久人妻99蜜桃| 99久久99久久久精品蜜桃| 丰满人妻一区二区三区视频av | 男人舔奶头视频| 久久久久国内视频| 国产免费男女视频| 听说在线观看完整版免费高清| 亚洲人成网站高清观看| 色综合亚洲欧美另类图片| 桃色一区二区三区在线观看| 国产精品综合久久久久久久免费| 欧美乱妇无乱码| 激情在线观看视频在线高清| 国产精品爽爽va在线观看网站| 久久人妻av系列| 色在线成人网| xxxwww97欧美| 国产爱豆传媒在线观看| 午夜精品在线福利| 国产精品98久久久久久宅男小说| av视频在线观看入口| 一区二区三区国产精品乱码| 老汉色av国产亚洲站长工具| 久久久久久久亚洲中文字幕 | 亚洲精品影视一区二区三区av| 欧美乱码精品一区二区三区| 床上黄色一级片| 高清在线国产一区| 国产不卡一卡二| 非洲黑人性xxxx精品又粗又长| 熟女电影av网| 中亚洲国语对白在线视频| 在线国产一区二区在线| 国产成人福利小说| 日韩大尺度精品在线看网址| 深夜精品福利| 两个人看的免费小视频| 久久久成人免费电影| 国产97色在线日韩免费| 大型黄色视频在线免费观看| 精品熟女少妇八av免费久了| 国产极品精品免费视频能看的| 嫩草影院入口| 免费无遮挡裸体视频| 天天添夜夜摸| 国产精品野战在线观看| 国产黄a三级三级三级人| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 香蕉av资源在线| 国产v大片淫在线免费观看| 男人舔奶头视频| 久9热在线精品视频| 波野结衣二区三区在线 | 午夜日韩欧美国产| 国产精品日韩av在线免费观看| 一本精品99久久精品77| 高清日韩中文字幕在线| 中文字幕人妻熟人妻熟丝袜美 | 中文字幕人妻熟人妻熟丝袜美 | 高潮久久久久久久久久久不卡| 国产精品av视频在线免费观看| 久久久国产成人精品二区| 国产三级中文精品| 热99re8久久精品国产| av在线蜜桃| 久久久久九九精品影院| 国产精品99久久99久久久不卡| 观看美女的网站| 18+在线观看网站| 99在线人妻在线中文字幕| 天堂影院成人在线观看| 国产欧美日韩一区二区三| 99热这里只有精品一区| 国产一区二区在线观看日韩 | 日韩欧美 国产精品| 一边摸一边抽搐一进一小说| 波野结衣二区三区在线 | 久久久久免费精品人妻一区二区| 色综合站精品国产| 国产伦精品一区二区三区四那| 欧美黄色片欧美黄色片| 亚洲在线自拍视频| 中文字幕人妻丝袜一区二区| 精品一区二区三区av网在线观看| 国产精品一区二区三区四区久久| 中文字幕久久专区| 19禁男女啪啪无遮挡网站| 亚洲av成人精品一区久久| 日韩欧美一区二区三区在线观看| 亚洲av不卡在线观看| 日韩成人在线观看一区二区三区| 久久精品国产99精品国产亚洲性色| 亚洲 国产 在线| 在线观看日韩欧美| 又粗又爽又猛毛片免费看| 一本一本综合久久| 亚洲一区二区三区色噜噜| 男人的好看免费观看在线视频| 国产精品99久久99久久久不卡| 99国产精品一区二区蜜桃av| 老司机在亚洲福利影院| 午夜福利在线观看吧| 岛国在线观看网站| 亚洲黑人精品在线| 国语自产精品视频在线第100页| 女人十人毛片免费观看3o分钟| 国产免费男女视频| 亚洲av熟女| 757午夜福利合集在线观看| 18禁美女被吸乳视频| 一个人看视频在线观看www免费 | 国产一级毛片七仙女欲春2| 久久久久九九精品影院| 国产一区二区三区视频了| 国产精华一区二区三区| 久久精品91蜜桃| 日本一二三区视频观看| 国产免费一级a男人的天堂| 亚洲成av人片免费观看| 每晚都被弄得嗷嗷叫到高潮| 国产成人av教育| 欧美+亚洲+日韩+国产| 一进一出抽搐动态| 午夜免费成人在线视频| 国产单亲对白刺激| 最新中文字幕久久久久| 国内久久婷婷六月综合欲色啪| 99久久精品热视频| 99久久成人亚洲精品观看| 丰满的人妻完整版| 亚洲av成人精品一区久久| 99精品欧美一区二区三区四区| 国产成人a区在线观看| 五月玫瑰六月丁香| 国产欧美日韩精品一区二区| 亚洲av电影在线进入| 色综合婷婷激情| 国产精品自产拍在线观看55亚洲| 欧美乱码精品一区二区三区| а√天堂www在线а√下载| 成人高潮视频无遮挡免费网站| 成人午夜高清在线视频| 波野结衣二区三区在线 | 国产精品国产高清国产av| 18禁美女被吸乳视频| 欧美日韩瑟瑟在线播放| 亚洲五月天丁香| 小说图片视频综合网站| 精品欧美国产一区二区三| 两个人的视频大全免费| 无人区码免费观看不卡| 国产野战对白在线观看| 成人亚洲精品av一区二区| 久久亚洲精品不卡| 日本一本二区三区精品| 国产一区二区亚洲精品在线观看| netflix在线观看网站| 免费大片18禁| 手机成人av网站| 久久天躁狠狠躁夜夜2o2o| 欧美最新免费一区二区三区 | 三级毛片av免费| 婷婷精品国产亚洲av在线| 亚洲一区高清亚洲精品| 国产一区二区在线观看日韩 | 一级a爱片免费观看的视频| 亚洲第一欧美日韩一区二区三区| www日本黄色视频网| 成人亚洲精品av一区二区| 国产高清激情床上av| av黄色大香蕉| 亚洲熟妇熟女久久| 99久久99久久久精品蜜桃| 香蕉丝袜av| 少妇熟女aⅴ在线视频| 母亲3免费完整高清在线观看| 亚洲精品一卡2卡三卡4卡5卡| 免费搜索国产男女视频| 三级毛片av免费| 国产精品久久久久久久电影 | 日韩精品中文字幕看吧| 亚洲av免费在线观看| 少妇的逼好多水| 色视频www国产| 亚洲国产高清在线一区二区三| 国内毛片毛片毛片毛片毛片| 一级毛片女人18水好多| 又粗又爽又猛毛片免费看| 黄色视频,在线免费观看| 小蜜桃在线观看免费完整版高清| 老汉色av国产亚洲站长工具| 51午夜福利影视在线观看| 久久久久久久亚洲中文字幕 | 国产午夜精品久久久久久一区二区三区 | 日日夜夜操网爽| 丰满的人妻完整版| 蜜桃久久精品国产亚洲av| 女人高潮潮喷娇喘18禁视频| 亚洲av电影在线进入| 天堂网av新在线| 丰满的人妻完整版| 午夜福利在线观看免费完整高清在 | 夜夜爽天天搞| 久久精品国产综合久久久| 黄片小视频在线播放| 亚洲avbb在线观看| 亚洲自拍偷在线| aaaaa片日本免费| 国产精品综合久久久久久久免费| 亚洲精品国产精品久久久不卡| 高清在线国产一区| 精华霜和精华液先用哪个| 精品无人区乱码1区二区| 免费在线观看影片大全网站| 美女cb高潮喷水在线观看| 免费av不卡在线播放| 成人午夜高清在线视频| 成人永久免费在线观看视频| 极品教师在线免费播放| 中文字幕av在线有码专区| 成年女人永久免费观看视频| 日韩国内少妇激情av| 久久久久国内视频| 大型黄色视频在线免费观看| 欧美成人免费av一区二区三区| 黄色女人牲交| 看免费av毛片| 国产伦人伦偷精品视频| 夜夜看夜夜爽夜夜摸| 无人区码免费观看不卡| 男女下面进入的视频免费午夜| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 香蕉丝袜av| 免费av不卡在线播放| 一级作爱视频免费观看| 91av网一区二区| 亚洲天堂国产精品一区在线| 欧美性猛交黑人性爽| www日本在线高清视频| 在线观看午夜福利视频| 国产爱豆传媒在线观看| 色av中文字幕| 老汉色av国产亚洲站长工具| 老司机深夜福利视频在线观看| 中亚洲国语对白在线视频| av在线蜜桃| 非洲黑人性xxxx精品又粗又长| 美女免费视频网站| 亚洲在线观看片| 香蕉av资源在线| 国产 一区 欧美 日韩| 久久久成人免费电影| 老司机午夜十八禁免费视频| ponron亚洲| 亚洲av免费高清在线观看| 亚洲专区国产一区二区| 欧美在线一区亚洲| 香蕉久久夜色| 黄色视频,在线免费观看| 免费一级毛片在线播放高清视频| 九九久久精品国产亚洲av麻豆| 亚洲人成网站在线播| 少妇丰满av| 99热这里只有精品一区| 国产精品av视频在线免费观看| 亚洲男人的天堂狠狠| 欧美bdsm另类| 亚洲av不卡在线观看| 一a级毛片在线观看| 校园春色视频在线观看| 五月伊人婷婷丁香| 欧美日韩亚洲国产一区二区在线观看| 欧美bdsm另类| 法律面前人人平等表现在哪些方面| 国产主播在线观看一区二区| 日本三级黄在线观看| 久久久成人免费电影| 99久久成人亚洲精品观看| 国产精品一及| 国产精品野战在线观看| 夜夜爽天天搞| 嫩草影视91久久| 一卡2卡三卡四卡精品乱码亚洲| 国产成年人精品一区二区| 国产伦一二天堂av在线观看| 十八禁人妻一区二区| 欧美又色又爽又黄视频| 精品久久久久久久久久久久久| 欧美一区二区国产精品久久精品| 亚洲av二区三区四区| 亚洲成人久久爱视频| 精品国产三级普通话版| 一区二区三区国产精品乱码| 亚洲精华国产精华精| 97碰自拍视频| 天堂影院成人在线观看| 深爱激情五月婷婷| 国产精品久久久久久久电影 | 国产一区在线观看成人免费| 国产色婷婷99| 岛国在线免费视频观看| 成人av一区二区三区在线看| 熟女少妇亚洲综合色aaa.| 亚洲精品国产精品久久久不卡| 欧美不卡视频在线免费观看| 无人区码免费观看不卡| 欧美午夜高清在线| 国产伦在线观看视频一区| 给我免费播放毛片高清在线观看| 日韩欧美国产一区二区入口| 久久精品国产99精品国产亚洲性色| 欧美bdsm另类| www.熟女人妻精品国产| 18禁国产床啪视频网站| 国产高清videossex| 男女床上黄色一级片免费看| 中文字幕久久专区| 国产一区二区三区在线臀色熟女| 亚洲va日本ⅴa欧美va伊人久久| 51午夜福利影视在线观看| 国语自产精品视频在线第100页| 中亚洲国语对白在线视频| av天堂中文字幕网| 欧美成人一区二区免费高清观看| 国产精品98久久久久久宅男小说| 青草久久国产| 每晚都被弄得嗷嗷叫到高潮| 少妇人妻精品综合一区二区 | 国产av一区在线观看免费| 亚洲av免费在线观看| av专区在线播放| 国产精品女同一区二区软件 | 可以在线观看的亚洲视频| 国产av在哪里看| 欧美丝袜亚洲另类 | 狂野欧美激情性xxxx| 午夜福利18| 亚洲精品乱码久久久v下载方式 | 在线播放国产精品三级| 欧美黑人欧美精品刺激| 人妻丰满熟妇av一区二区三区| av在线蜜桃| 亚洲国产中文字幕在线视频| 18禁国产床啪视频网站| 免费观看精品视频网站| a级毛片a级免费在线| 国产亚洲精品久久久久久毛片| 又爽又黄无遮挡网站| 哪里可以看免费的av片| 免费av不卡在线播放| 给我免费播放毛片高清在线观看| 午夜福利欧美成人| 小蜜桃在线观看免费完整版高清| 高清毛片免费观看视频网站| 18禁黄网站禁片免费观看直播| 欧美最黄视频在线播放免费| 亚洲五月天丁香| 午夜福利18| 国产欧美日韩精品一区二区| 欧美黑人巨大hd| 日本五十路高清| 日韩大尺度精品在线看网址| 嫩草影院入口| 亚洲自拍偷在线| 国产aⅴ精品一区二区三区波| 日本a在线网址| 国产三级中文精品| 国产aⅴ精品一区二区三区波| 国产高清有码在线观看视频| 欧美一区二区精品小视频在线| 级片在线观看| 久久精品国产清高在天天线| 亚洲中文日韩欧美视频| 成年人黄色毛片网站| 亚洲精品日韩av片在线观看 | 日韩欧美 国产精品|