岩石学报  2012, Vol. 28 Issue (1): 122-132   PDF    
江西黄沙石英脉型钨矿床流体包裹体研究
王旭东1,2, 倪培2, 袁顺达1, 吴胜华1     
1. 中国地质科学院矿产资源研究所, 国土资源部成矿作用和资源评价重点实验室, 北京 100037;
2. 南京大学内生金属矿床成矿机制研究国家重点实验室, 地质流体研究所, 地球科学与工程学院, 南京 210093
摘要: 黄沙钨矿床是赣南地区一大型石英脉型钨多金属矿床。本文采用“流体包裹体组合”的研究方法,对黄沙钨矿床主成矿阶段早期的黑钨矿-石英脉和晚期的硫化物-(黑钨矿)-石英脉石英中的流体包裹体进行了显微测温和拉曼探针的分析。研究表明,黑钨矿-石英脉中包裹体主要为水溶液包裹体和含CO2水溶液包裹体,硫化物-(黑钨矿)-石英脉中主要发育水溶液包裹体。黑钨矿-石英脉中包裹体的均一温度明显高于硫化物-(黑钨矿)-石英脉中的包裹体,但两者水溶液包裹体的盐度相差不大。激光拉曼探针测试表明,两期矿脉中水溶液包裹体的组分主要为水,在黑钨矿-石英脉中的含CO2水溶液包裹体,除CO2外,还检测到CH4和N2组分。研究表明,以CO2逸失为特征的流体不混溶作用是早期黑钨矿-石英脉含矿流体中的金属络合物分解并沉淀成矿的主要机制,晚期硫化物-(黑钨矿)-石英脉中矿质的沉淀则主要是流体的混合作用导致。
关键词: CO2包裹体     流体不混溶     流体混合     石英脉型钨矿床     江西黄沙    
Fluid inclusion studies of the Huangsha quartz-vein type tungsten deposit, Jiangxi Province
WANG XuDong1,2, NI Pei2, YUAN ShunDa1, WU ShengHua1     
1. MLR Key Laboratory of Metallogeny and Mineral Assessment, Institute of Mineral Resources, Chinese Academy of Geological Sciences, Beijing 100037, China;
2. State Key Laboratory for Mineral Deposit Research, Institute of Geo-Fluids, School of Earth Science and Engineering, Nanjing University, Nanjing 210093, China
Abstract: The Huangsha tungsten deposit is an important large quartz-vein type W-polymetallic deposit in southern Jiangxi Province. Based on detailed petrographic observations and using the means of Fluid Inclusion Assemblage (FIA), we carried out microthermometric and Raman microspectroscopic studies of fluid inclusions in early wolframite-quartz veins and late sulfide-(wolframite)-quartz veins of the main ore-forming stage of Huangsha tungsten deposit. The fluid inclusions are mainly aqueous inclusions and carbon-aqueous inclusions in wolframite-quartz veins and only aqueous inclusions in sulfide-(wolframite)-quartz veins. Experimental results show that the aqueous inclusions decrease in homogenization temperature from the early veins to the late veins at constant salinity. Raman microspectroscopic studies of the fluid inclusions show that, the main component of aqueous inclusions in wolframite-quartz veins and sulfide-(wolframite)-quartz veins is H2O, whereas those in wolframite-quartz veins are major H2O, minor CO2, and rare CH4 and N2. The characteristics of fluid inclusions in the Huangsha tungsten deposit indicate that the fluid immiscibility caused by CO2 escaping possibly led to metals precipitation in the wolframite-quartz veins and the fluid mixing led to metals precipitation in the wolframite-quartz veins.
Key words: CO2-rich fluid inclusion     Fluid immiscibility     Fluid mixing     Quartz-vein type tungsten deposit     Huangsha, Jiangxi    

赣南是我国最重要的黑钨矿生产基地和资源基地,也是世界著名的钨矿资源集中区之一。黄沙钨矿是区内一大型黑钨矿石英脉型钨多金属矿床,位于江西省于都县丰田镇。前人对黄沙钨矿开展了较为详尽的研究工作(夏宏远等,19821984项尝培,1985夏宏远和梁书艺,1987陈尊达和胡力檖,1984干国梁,1988冯志文等,1989夏宏远,1989任英忱,1998於崇文,2004陈祥,2007),取得了一系列的认识和成果。但针对矿床流体包裹体的研究相对较少,相关的研究工作也限于流体包裹体的显微测温和利用包裹体群分析技术获取流体的信息方面,也没有基于流体包裹体研究的角度解释矿床的成因机制。随着测试技术和研究方法的进步,对该矿床的流体包裹体进行深入、系统的研究工作可以获得更为准确的流体信息,并有助于对矿床的成因机制进行合理的解释。

流体包裹体组合(FIA)理论是近年来流体包裹体研究的重要进展之一,FIA是指通过岩相学方法能够分辨出来的、代表了一个在时间上分得最细的包裹体捕获事件的一组包裹体(Goldstein and Reynolds 1994Goldstein,2003)。FIA的研究方法可以使测试的数据更具有效性、数据的结果更具代表性(池国祥和卢焕章,2008)。

本文以石英中流体包裹体为研究对象,在详细的岩相学观察的基础上,采用流体包裹体组合(FIA)的研究方法对黄沙钨矿床主要成矿阶段,即黑钨矿-石英脉阶段和硫化物-(黑钨矿)-石英阶段石英脉中流体包裹体进行了显微测温和激光拉曼光谱学测试,并据此对黄沙钨矿成矿流体的性质和流体中金属元素的沉淀机制进行探讨。

1 区域与矿床地质

赣南地区在大地构造上处于欧亚板块与滨西太平洋板块消减带内侧的华夏板块中,武夷山和南岭两大成矿带的交汇复合部位(朱焱龄等,1981)。区内出露自新元古界神山群至第四系地层(江西省地质矿产局,1984),其中震旦纪-寒武纪基底岩系及泥盆纪地层具有高于地壳克拉克值几倍至几十倍钨含量(徐克勤等,1984a韩久竹等,1984),被认为与区内钨的成矿作用有着密切的关系。

赣南的岩浆活动强烈,岩性以不同时代、不同类型的花岗岩类为主,其中燕山期花岗岩类与本区钨成矿关系密切,该花岗岩类主要是以高硅、富碱、富挥发分、铝过饱和为特征(徐克勤等,1984b)。

区内矿产资源丰富,尤以钨矿床(点)大量产出为特征,钨矿床(点)主要汇集于区内西南部的诸广山聚集带,中东部于山聚集带和南部的九连山聚集带。区内钨矿床类型包括石英脉型、变花岗岩型、云英岩型、矽卡岩型、蚀变破碎带型等,其中石英脉型为最重要的矿床类型(康永孚和李崇佑,1991)。黄沙钨矿处于于山聚集带中,南与盘古山、北与隘上、西与上坪等钨矿相邻。(图 1)。

图 1 于山钨矿床聚集带地质简图(据盛继福,2004) 1-元古宇变余凝灰岩、凝灰质砂岩、变余长石石英砂岩;2-下古生界变余长石石英砂岩、粉砂质板岩;3-上古生界-中三叠统石英砂岩、粉砂岩、炭质页岩;4-上三叠统-白垩系泥岩、粉砂岩、长石石英砂岩、凝灰质页岩、砾岩、凝灰质砂岩;5-燕山期花岗岩;6-地质界线;7-断裂;8-石英脉型钨矿;9-矽卡岩型钨矿 Fig. 1 Geological sketch map of the Yushan tungsten ore concentration belt(modified after Sheng,2004) 1-metamorphic tuff,tuffaceous sandstone and metamorphic feldspar-quartz sandstone of Proterozoic; 2-metamorphic feldspar-quartz sandstone and silty slate of Lower Paleozoic; 3-quartz sandstone,siltstone and carbonaceous shale of Upper Paleozoic-Middle Triassic; 4-mudstone,siltstone,feldspar-quartz sandstone,tuffaceous slate,conglomerate and tuffaceous sand of Upper Triassic-Cretaceous; 5-Yanshanian granite; 6-geological boundary; 7-fault; 8-quartz vein type tungsten deposit; 9-skarn type tungsten deposit

矿区内出露地层主要为寒武系浅变质砂岩和板岩,其次是泥盆系上统中棚组及第四系,其中寒武系的浅变质岩系是主要的赋矿围岩。对寒武系变质岩的主要成矿元素的定量分析结果表明,钨、锡、铋、铍等成矿元素含量均是克拉克值20倍以上(陈祥,2007)。

矿区内岩浆活动表现为燕山期花岗岩浆的多次侵入,形成了铁山垅复式花岗岩体和黄沙隐伏复式花岗岩体及少量呈岩脉或岩墙状产出的花岗斑岩。铁山垄花岗岩岩株出露于矿区东北部,北东向断层和北西向断裂交接复合处。平面上呈椭圆形,长轴为南北向,出露面积约24km2,该复式岩体主体为中粗粒似斑状黑云母花岗岩,夏宏远和梁书艺(1987) 曾报道其同位素K-Ar年龄为184Ma和177Ma,属燕山早期产物,在牛角山、铜岭、杨坑山等地也发育有细粒二云母花岗岩以及白云母花岗岩,呈岩滴、岩瘤状产出于主岩体之中,构成补体花岗岩。

黄沙矿区深部隐伏的花岗岩为铁山垅花岗岩的西南延伸部分,在隐伏花岩的突出部位云英岩化广泛发育,其上部强烈处形成云英岩型W(Mo)矿体。云英岩化花岗岩白云母的同位素K-Ar年龄为153.2Ma(夏宏远和梁书艺,1987)。在云英岩与围岩的接触处,局部可见似伟晶岩(长石、石英)壳发育,并有花岗岩脉产于变质岩中。往下则云英岩化减弱而钠长石化增强,铌钽矿化发育,部分可达到工业品位。

花岗斑岩呈岩脉或岩墙出露于矿区北部和南部,在14中段见早期花岗岩碎块以俘虏体的形式赋存于该期花岗斑岩之中,花岗斑岩沿EW向断层断续分布。

对区内出露岩浆岩系统的地球化学研究表明(夏宏远和梁书艺,1987李光来,2011),主体花岗岩与补体花岗岩可能为同源岩浆演化不同阶段的产物,花岗岩的分异演化可能对成矿元素的预富集发挥重要作用。而晚期就位的花岗斑岩演化程度更低,与主体花岗岩、补体花岗岩并非同源岩浆演化的产物。

矿区内断裂异常发育,以近东西向为主,北东向和南北向次之。近东西向裂隙为规模较大的剪切裂隙,北东向和南北向裂隙为张性裂隙,规模较小。前者分布于矿区中部,形成较早;后者分布于矿区南北两侧,形成较晚,以上各组裂隙与区内分布的矿脉相对应。矿区内各组构造裂隙皆遭受不同程度的矿化,其中以F3和北部的东西向断层矿化较强,与成矿作用密切相关。

黄沙钨矿床为典型的石英脉型黑钨矿床,据矿体的产出特征,可分为石英细脉带型和单脉型,其中石英细脉带型为本矿床主要的工业矿体,单脉型储量较小。石英脉(带)赋存标高在300m到800m之间,延长300m到1300m不等,倾向南或北,倾角60°~85°,走向分为NE向-NEE向、EW向、NWW 向和SN向。

自北向南分别为樟木林组、北组、中组、青山窝组、芭蕉坑组、南组、南北脉组和内带盲脉组(本组全产于花岗岩体中)(图 2图 3)。其中前五组多呈脉带产出,在剖面上具典型的“五层楼”(古菊云,1984)形态分布特征,钨矿床的主要储量存在于此;后三组多呈单脉产出,所含的储量相对较少。

图 2 黄沙钨矿地质平面略图(据夏宏远等,1984) 1-寒武系变质砂板岩;2-上泥盆统砂砾岩;3-第四系沉积物;4-黑云母花岗岩; 5-花岗斑岩墙;6-石英线脉带;7-石英线脉-细脉带;8-石英细脉-大脉带;9-断层;10-背斜轴 Fig. 2 The geological plane of the Huangsha tungsten deposit(modified after Xia et al.,1984) 1-Quaternary sediments; 2-sandy conglomerate of Upper Devonian; 3-metamorphic sandstone and slate of Cambrian; 4-biotite granite; 5-dyke of granite porphyry; 6-line quartz veins belt; 7-line and thin quartz veins belt; 8-thin and wide quartz veins belt; 9-fault; 10-anticline axis

图 3 黄沙钨矿地质剖面略图(据夏宏远,1989) 1-寒武系变质砂板岩; 2-花岗斑岩墙;3-黑云母花岗岩;4-断层;5-含矿石英脉带 Fig. 3 The geological profile of the Huangsha tungsten deposit(modified after Xia,1989) 1-metamorphic sandstone and slate of Cambrian; 2-dyke of granite porphyry; 3-biotite granite; 4-fault; 5-ore-bearing quartz veins belt

黄沙钨矿床的成矿作用具有多阶段的特征,前人根据矿物组合及矿脉之间的穿插关系,将该矿床从早到晚划分为:硅酸盐-氧化物阶段、黑钨矿-石英阶段、硫化物-(黑钨矿)-石英阶段、条带状硫化物石英阶段和萤石-碳酸盐阶段等五个矿化阶段,其中黑钨矿-石英阶段、硫化物-(黑钨矿)-石英阶段为主要的矿化阶段。在矿区内,这两个矿化阶段比较容易辨认,除在有些区段硫化物-(黑钨矿)-石英脉切割黑钨矿-石英脉外,矿脉(带)的中上部硫化物矿化明显,中下部黑钨矿的矿化发育是矿床的一个重要的特征。

矿床中辉钼矿Re-Os同位素年龄为146.4±6.1Ma,辉钼矿中的Re含量,与南岭东部石英脉型钨矿中伴生的辉钼矿的Re的较低含量接近,显示了成矿物质的壳源特征(Mao et al.,1999; 李光来,2011)。

矿床中矿石矿物以黑钨矿为主,次为白钨矿、黄铜矿、辉钼矿、辉铋矿、锡石、黄铁矿、闪锌矿和含银硫盐等;脉石矿物主要为石英,次为白云母、锂白云母、黄玉、钾长石、萤石、方解石等,其中黑钨矿主要以自形针柱状、长柱状、短柱状、以及他形不规则状产出。矿石结构有自形-半自形粒状结构和他形粒状结构、碎裂结构等。矿石构造主要为囊状构造、块状构造、条带状构造、晶洞构造等。

围岩蚀变类型有白云母化、硅化、电气石化、绢云母化、黄铁矿化、萤石化、黄玉化及绿泥石化等,围岩蚀变在垂向上具有分带的特征,如矿体上部的白云母化较强,而硅化相对较弱;往深部则硅化、黄铁矿化增强。

2 样品特征与分析方法

本次研究的矿石样品采自黄沙钨矿北组、中组、青山窝组和芭蕉坑组矿化特征突出的黑钨矿-石英阶段和硫化物-(黑钨矿)-石英阶段的石英脉。两类矿石在矿化上有较大差别,黑钨矿-石英脉矿石,主要是黑钨矿矿化(图 4a),硫化物-(黑钨矿)石英脉矿石,黑钨矿矿化较弱,以强烈的硫化物矿化为特征(图 4b)。

图 4 黄沙钨矿黑钨矿-石英脉矿石(a)和硫化物-(黑钨矿)-石英脉矿石(b) Fig. 4 Specimen of wolframite-quartz veins(a)and sulfide-(wolframite)-quartz veins(b)in the Huangsha tungsten deposit

本次研究先将这些样品磨制成厚度约为0.3mm双面抛光的薄片,以流体包裹体组合(FIA)(Goldstein and Reynolds,1994Goldstein,2003)的方法对样品石英中流体包裹体进行详细的岩相学观察及显微测温,然后选择有代表性的包裹体进行激光拉曼探针分析。流体包裹体测试在南京大学内生金属成矿机制研究国家重点实验室包裹体室进行,显微测温所采用仪器为英国产的Linkam-THMS600冷热台(温度范围:-195~+600℃),分析精度为:±0.2℃,<30℃,±1℃,<300℃,±2℃,<600℃。激光Raman探针为英国产Renishaw RM2000型,实验条件:温度23℃,Ar离子激光器(514nm),风冷,狭缝宽50μm,光栅1800,扫描时间60s,扫描次数为1次。

3 流体包裹体研究 3.1 流体包裹体类型

黄沙钨矿床两个主要矿化阶段石英脉石英中均发育大量的流体包裹体。根据Roedder(1984) 卢焕章等(2004) 提出的流体包裹体在室温下相态分类准则及冷冻回温过程中的相态变化,可将流体包裹体划分为富液相两相水溶液包裹体(Ⅰ型)、富气相两相水溶液包裹体(Ⅱ型)、纯气相包裹体(Ⅲ型)、含液相CO2的三相水溶液包裹体(Ⅳ型)四类。

Ⅰ型:富液相两相水溶液包裹体,此类包裹体占石英中流体包裹体的总量90%以上,孤立状(图 5a)或成小群分布(图 5b)。包裹体长径长0.5~50μm,绝大多数为5~15μm,气相百分数5%~10%不等,形状一般为不规则状、椭圆形、长条形或石英负晶形。

图 5 黄沙钨矿石英中流体包裹体 Fig. 5 Fluid inclusions in the Huangsha tungsten deposit

Ⅱ型:富气相两相水溶液包裹体,此类包裹体数量非常少,个体通常较小,长径长5~10μm,形状一般为圆形、椭圆形,气相百分数在70%~90%之间,一般和Ⅰ型包裹体相伴生(图 5c)。

Ⅲ型:纯气相包裹体,同Ⅱ型包裹体一样,此类包裹体数量也非常少,形状为不规则状、长条形,长径长2~10μm,呈单一相,颜色较暗,与Ⅰ型包裹体共生或孤立状产出(图 5d)。

Ⅳ型:含液相CO2的三相水溶液包裹体,在本次研究中,此类包裹体在硫化物-(黑钨矿)-石英脉中没有发现,在黑钨矿-石英脉中则较为发育,该类包裹体中CO2相的相比例均较大,从占包裹体总体积的接近50%到80%不等,长径长8~15μm,形态呈圆形、椭圆形,空间上与Ⅰ型包裹体共生产出(图 5e,f)。

3.2 流体包裹体显微测温

在详细的岩相学观察的基础上,对流体包裹体开展了显微测温工作,因Ⅱ型包裹体数量非常少,且个体较小,在测试中难以对此类包裹体进行准确的测温工作,本次也没有获得此类包裹体准确的相变温度数据。对Ⅲ型包裹体的冷冻加热过程中没有观察到相态的变化。测试的主要对象是原生和假次生的Ⅰ型和Ⅳ型包裹体。

Ⅰ型包裹体盐度的计算采用Hall et al.(1988)的公式:WNaCl=0.00+1.78Tm-0.0442 Tm2+0.000557 Tm3,其中Tm为冰点下降温度。Ⅳ型包裹体盐度的计算据Roedder(1984) 提供的公式:WNaCl=15.52022-1.02342×T-0.05286×T2,其中T为笼形物融化温度。测温及盐度计算结果见表 1,利用测试得到的和计算得到的温度盐度数据作图见图 6图 7

表 1 黄沙钨矿流体包裹体测温结果 Table 1 Microthermometric data for fluid inclusions in the Huangsha tungsten deposit

图 6 黄沙钨矿流体包裹体均一温度直方图 1-黑钨矿-石英脉中Ⅰ型包裹体;2-Ⅳ型包裹体;3-硫化物-(黑钨矿)-石英脉中Ⅰ型包裹体 Fig. 6 Histogram of homogenization temperature in fluid inclusions from the Huangsha tungsten deposit 1-type I fluid inclusions in wolframite-quartz veins; 2-type Ⅳfluid inclusions; 3-type I fluid inclusions in sulfide-(wolframite)-quartz veins

图 7 黄沙钨矿流体包裹体盐度直方图 1-黑钨矿-石英脉中Ⅰ型包裹体;2-Ⅳ型包裹体;3-硫化物-(黑钨矿)-石英脉中Ⅰ型包裹体 Fig. 7 Histogram of salinity in fluid inclusions from the Huangsha tungsten deposit 1-type I fluid inclusions in wolframite-quartz veins; 2-type Ⅳfluid inclusions; 3-type I fluid inclusions in sulfide-(wolframite)-quartz veins

测温及计算结果表明,黑钨矿-石英阶段中Ⅰ型包裹体均一温度的分布区间为256~336℃,盐度为4%~8% NaCleqv,Ⅳ型包裹体的均一温度区间为298~351℃,盐度为1.4%~2.4% NaCleqv;硫化物-(黑钨矿)-石英阶段中Ⅰ型包裹体均一温度的主要分布区间为190~240℃,盐度主要分布区间为1%~7%NaCleqv。黑钨矿-石英阶段中Ⅳ型包裹体的三相点温度为-58.3~-64.4℃,低于纯CO2的三相点(-56.6℃),表明除CO2外,还混有其他挥发组分(Shepherd et al.,1985),这也被其后的流体包裹体激光拉曼测试所证实。在均一方式上,两个成矿阶段中Ⅰ型包裹体的均一方式都为均一到液相。一些CO2相比例较大的Ⅳ型包裹体在完全均一前爆裂,可测到完全均一温度的Ⅳ型包裹体的均一方式为均一到气相。

3.3 流体包裹体的激光拉曼探针测定

选择黄沙钨矿黑钨矿-石英脉和硫化物-(黑钨矿)-石英脉中具有代表性的、不同类型的流体包裹体进行了拉曼探针分析。 测试结果显示,在黑钨矿-石英脉Ⅰ型包裹体中主要检测到宽泛的液相H2O的包络峰(图 8a),该阶段石英脉中的Ⅳ型包裹体除检测到CO2外,均不同程度的检测到CH4和N2,相对于CH4,表现在拉曼谱图上出现典型的CO2谱峰,典型的N2谱峰,以及CH4谱峰(图 8b,c)从拉曼谱图上上看,个别包裹体中CH4含量很高,但包裹体中N2的含量均较少。

图 8 黄沙钨矿黑钨矿-石英脉流体包裹体拉曼图谱 Fig. 8 Raman spectra of fluid inclusions in wolframite-quartz veins of Huangsha tungsten deposit

在硫化物-(黑钨矿)-石英脉Ⅰ型包裹体中主要检测到宽泛的液相H2O的包络峰(图 9a,b),没有检测到CO2组分,在极少数Ⅰ型包裹体中检测到了少量CH4组分(图 9c)。

图 9 黄沙钨矿硫化物-(黑钨矿)-石英脉流体包裹体拉曼图谱 Fig. 9 Raman spectra of fluid inclusions in sulfide-(wolframite)-quartz veins of the Huangsha tungsten deposit
4 讨论 4.1 成矿流体的性质

流体包裹体测试结果表明,黄沙钨矿主要成矿阶段中早期的黑钨矿-石英脉的成矿流体为中-高温、低盐度的NaCl-H2O-CO2-CH4流体体系;晚阶段的硫化物-(黑钨矿)-石英脉的成矿流体则为中-低温、低盐度的NaCl-H2O流体体系。从成矿早阶段到晚阶段,流体的温度降低明显,盐度变化不大。

因黑钨矿-石英脉中Ⅳ型包裹体中不同程度的检测到CH4组分,利用该类包裹体CO2的熔化温度和部分均一温度作图(图 10),从图上可以看出,该成矿阶段Ⅳ型包裹体的摩尔体积除个别大于70cm3/mol外,主要在55~60cm3/mol,部分包裹体中CH4的摩尔分数介于0.07~0.12之间,但也有部分点投在了CH4-CO2的临界曲线之上,CH4的摩尔分数的变化,可能与流体的不混溶作用有关(范宏瑞等,2005)。本文没有对流体中CH4的来源开展进一步的工作,CH4可能来自初始的岩浆流体,也可能是在复杂的流体演化过程中,与H2渗透进入到流体包裹体内发生再平衡作用有关(胡芳芳等,2007)。

图 10 黄沙钨矿黑钨矿-石英脉阶段Ⅳ型包裹体CO2部分均一温度(ThCO2)和熔化温度(TmCO2)关系图(底图据 van den Kerkhof,1990; Thiéry et al.,1994) TcCO2-纯CO2临界温度;Tp CO2-纯CO2三相点 Fig. 10 Homogenization(ThCO2)vs. melting temperature(TmCO2)of type Ⅳfluid inclusions in wolframite-quartz veins of Huangsha tungsten deposit(after van den Kerkhof,1990; Thiéry et al.,1994) TcCO2-critical temperature of CO2Tp CO2-triple point temperature of CO2

根据流体包裹体的显微测温数据,利用Flincor流体包裹体计算软件(Brown,1989)对黄沙钨矿成矿流体的密度和压力进行了计算。结果表明,黑钨矿-石英脉Ⅰ型包裹体流体密度介于0.73~0.82g/cm3之间,Ⅳ型流体包裹体CO2的摩尔分数为0.29~0.49,密度为0.68~0.77g/cm3,采用Ⅳ型流体包裹体对该阶段的流体压力进行了估算,该阶段的流体压力介于97.3~156.6MPa之间;硫化物-(黑钨矿)-石英脉中其Ⅰ型包裹体流体密度介于0.86~0.92g/cm之间,因该阶段含矿石英脉中不发育Ⅳ型流体包裹体,采用Ⅰ型包裹体流体对流体压力进行了参考性的估算,结果表明,该阶段的流体压力介于6.0~28.0MPa之间。可以看出,两个矿化阶段的流体压力差别较大,矿化晚阶段的流体压力较低。根据矿床的产出特征,采用相关公式(Shepherd et al.,1985)对成矿深度进行了估算,其相应两个成矿阶段的成矿深度分别为3.7~5.9km和0.2 ~1.1km。两个成矿阶段的流体压力和估算的成矿深度存在较大的差别,硫化物-(黑钨矿)-石英阶段流体压力偏低可能是计算压力时采用流体体系的不同所致。黑钨矿-石英脉阶段的成矿深度大于硫化物-(黑钨矿)-石英脉与本矿床硫化物矿化发生在浅部,黑钨矿矿化发生在深部的地质事实是相符的。同时作者也考虑,因存在早期矿脉被晚期矿脉截然切割的事实,两阶段成矿作用是否是两期相对独立的热液活动的结果,而这两期流体在温压地球化学特征上存在较为明显的不同?这有待于开展进一步的工作去研究。

4.2 流体的不混溶作用和流体混合作用

黄沙钨矿主成矿阶段较早的黑钨矿-石英脉中Ⅳ型包裹体较为发育,且该类包裹体和Ⅰ型包裹体共生于同一流体包裹体组合(FIA)的现象较为常见,显示了两者同时捕获的特征。在显微测温过程中,Ⅳ型包裹体和Ⅰ型包裹体表现出不同的均一方式,Ⅰ型包裹体均一到液相,Ⅳ型包裹体均一到气相,后者的均一温度略高、而盐度低于前者,以上现象表明其捕获后流体发生了不混溶作用(Shepherd et al.,1985)。造成富液相两相水溶液包裹体盐度高于含CO2水溶液包裹体的原因可解释为在发生流体不混溶时,由于压力和温度的降低,使得在较高压力和温度条件下溶解于流体中CO2相分离出来,并由于气体的逸失,导致剩余流体中的盐度的升高。

在流体包裹体均一温度-盐度相关图(图 11)上,可以明显的看出,黄沙钨矿黑钨矿-石英脉中的流体在演化过程中经历了不混溶的作用。

图 11 黄沙钨矿流体包裹体均一温度-盐度图 1-黑钨矿-石英脉中Ⅰ型包裹体;2-Ⅳ型包裹体;3-硫化物-(黑钨矿)-石英脉中Ⅰ型包裹体 Fig. 11 Plot of homogenization temperatures vs. salinities of fluid inclusion in Huangsha tungsten deposit 1-type I fluid inclusions in wolframite-quartz veins; 2-type Ⅳfluid inclusions; 3-type I fluid inclusions in sulfide-(wolframite)-quartz veins

本次研究在黑钨矿-石英脉中Ⅳ型包裹体还不同程度的检测到CH4组分。已有研究表明,CH4组分的加入可以使NaCl-H2O-CO2流体在更深的部位发生不混溶作用(Naden and Shepherd,1989)。

相对于较早的黑钨矿-石英脉,本次研究没有在硫化物-(黑钨矿)-石英脉中发现Ⅳ型包裹体,该类脉中Ⅰ型包裹体的盐度变化范围较大,且均一温度与盐度之间有较为明显的线性关系,从高温区间到低温区间,随着均一温度的降低,盐度也表现出相应的降低趋势(图 10),显示了一定的流体混合特征。表明了在流体演化过程中可能发生了较高温度、盐度的流体与较低温度、盐度流体的混合作用。该阶段流体混合作用也得到了H、O同位素研究的证实,冯志文等(1989) 干国梁(1988) 的研究表明,黄沙钨矿黑钨矿-石英脉中的成矿流体主要是岩浆水,在硫化物-(黑钨矿)-石英脉阶段,有向大气水过渡的趋势。李光来(2011) 通过开展该矿床流体包裹体He-Ar同位素的研究工作,也证实了在硫化物-(黑钨矿)-石英脉矿化阶段有大气水的加入。本次流体包裹体研究工作进一步证实了该成矿阶段流体的混合作用。

研究表明,钨在成矿流体中呈多种络合物的形式迁移(Wood,1992Wood and Samson,2000),热液型钨矿床的成矿作用,本质上就是流体中钨的络合物的分解、沉淀作用。

钨的络合物在热液中分解、沉淀的机制一般被解释成流体体系的自然冷却(Ramboz et al.,1985Seal et al.,1987Samson,1990O’Reilly et al.,1997)、不同流体的混合(Landis and Rye,1974Ramboz et al.,1985Jackson et al.,1989Heinrich,1990Bailly et al.,2002Beuchat et al.,2004Wang et al.,2010)、流体的不混溶(沸腾或CO2的泡腾)(Higgins and Kerrich,1982Ramboz et al.,1982Seal et al.,1987Lynch,1989Polya,1989Giamello et al.,1992So and Yun,1994Graupner et al.,1999)、pH值升高、水岩反应(Kelly and Rye,1979Seal et al.,1987Polya,1989Clark et al.,1990Catallani and Williams,1991)、压力的降低(Polya,19891990)及围岩中非极性挥发份的加入(O’Reilly et al.,1997)等。大量研究表明,不同流体的混合和流体的不混溶作用是含矿流体中钨的络合物的分解、沉淀的主要机制。

本次对黄沙钨矿床主成矿阶段的流体包裹体研究工作表明,其早阶段黑钨矿-石英脉的成矿流体在演化过程中发生了以CO2逸失为特征的流体不混溶作用,流体不混溶作用是该阶段含矿流体中络合物分解并沉淀成矿的主要因素。晚阶段硫化物-(黑钨矿)-石英脉中成矿元素沉淀的主导因素则是流体的混合作用。

4.3 CO2在成矿中作用的讨论

研究表明,不同程度的含有CO2组分是大多数钨矿床成矿流体的普遍特征(Landis and Rye,1974Kelly and Rye,1979Higgins,1980Bussink and Kreulen,1984Ramboz et al.,1985Kamilli et al.,1992Noronha et al.,1992Graupner et al.,1999Rios et al.,2003王旭东等,2010)但也有部分矿床不含或仅含微量的CO2(Rankin and Alderton,1983Zaw,1984Shelton et al.,1986Thorn,1988Vallance et al.,2001刘若兰等,1985Giuliani et al.,1988刘若兰和慕纪录,1990王旭东等,2008胡东泉等,2011)。就赣南地区而言,崇-犹-余矿集区的钨矿床含矿石英脉石英中一般不含或仅含微量的CO2等挥发分,而不同程度的含有CO2组分是九连山矿集区及于山矿集区中的钨矿床的显著特征,对于本次研究的黄沙钨矿床,主成矿阶段黑钨矿-石英脉流体中含有CO2组分,硫化物-(黑钨矿)-石英脉中则没有发现含CO2的包裹体。是什么样的地质因素造成这样的差别,仍需进一步的研究。

Higgins(1980) 在对前人工作总结的基础上,认为在高温高压的条件下,在富CO2流体中钨可能以碳酸盐、重碳酸盐的形式迁移,Ramboz et al.(1985) 则认为没有确定的证据表明钨是在含CO2的流体中迁移或是在水溶液中迁移,Wood and Samson(2000) 的研究表明,在一定的温压条件下,在NaCl-HCl-H2O溶液中,钨主要以简单钨酸(如H2WO4,HWO44-,WO2-和KHWO40)和碱性钨酸盐离子对(如NaHWO40,KWO4-和NaWO4-)形式迁移,钨的浓度可达数千×10-6的浓度,足以形成钨的热液矿床。

尽管CO2组分在钨在成矿流体中迁移所起到的作用尚需进一步的研究,但大多数钨矿床成矿流体中均不同程度的含CO2组分表明CO2是热液钨矿床形成的有利因素之一。

4.4 矿床的成因机制

黄沙钨矿主成矿阶段含矿石英脉石英中流体包裹体研究表明,黑钨矿-石英脉与硫化物-(黑钨矿)-石英脉的含矿流体具有不同的演化过程,并由此引起成矿元素沉淀成矿的方式的不同。早阶段的黑钨矿-石英脉中成矿元素的沉淀主要是含矿流体发生以CO2逸失为特征的不混溶作用导致,其可能的过程为:岩浆作用形成的并具有较大内压的含矿流体,在岩浆期后热动力的驱动下,沿矿区内的断裂向裂隙发育的低压地带运移,在此过程中,随着温度、压力等条件的改变,流体发生了以CO2等挥发分逸出为特征的不混溶作用,流体不混溶使含矿流体的酸碱性、氧化还原等条件发生改变,致使含矿流体中的络合物分解,并且由于CO2等挥发分的逸失导致残余流体浓度升高,甚至造成流体的瞬时过饱和。在这个过程中,金属络合物分解,WO42-与流体中的Fe2+、Mn2+等金属阳离子结合而沉淀成矿。晚阶段的硫化物-(黑钨矿)-石英脉中成矿元素的沉淀则主要由混合作用引起,其可能的过程为:上升的、富含成矿元素的岩浆流体与通过矿区内断裂和裂隙下渗的、被加热并萃取了围岩中一定量金属元素的大气水混合,混合作用引起流体体系物理化学条件的不断改变,在此过程中,金属络合物分解,在有利的构造地段沉淀成矿。

5 结论

(1) 黄沙钨矿主要矿化阶段中早阶段黑钨矿-石英脉成矿流体为中-高温、低盐度的NaCl-H2O-CO2-CH4体系,晚阶段硫化物-(黑钨矿)-石英脉的成矿流体为中-低温、低盐度的NaCl-H2O体系,从早阶段到晚阶段,流体的温度降低明显,盐度变化不大,两者的成矿压力及由此计算的成矿深度存在较大的差别,表明两者成矿条件的不同。

(2) 矿床主成矿阶段中早阶段的黑钨矿-石英脉的成矿流体在演化过程中,经历了以CO2逸失为特征的流体不混溶作用,流体不混容作用是导致含矿流体中的金属络合物分解并沉淀成矿的主要机制。晚阶段硫化物-(黑钨矿)-石英脉中矿质的沉淀则主要是流体的混合作用导致。

致谢 论文的写作得到了毛景文研究员的悉心指导;两位审稿人对论文初稿提出了宝贵的修改意见;在此一并表示衷心的感谢。
参考文献
[] Bailly L, Grancea L, Kouzmanov K. 2002. Infrared microthermometry and chemistry of wolframite from the baia sprie epithermal deposit, Romania. Economic Geology, 97: 415–421. DOI:10.2113/gsecongeo.97.2.415
[] Beuchat S, Moritza R, Pettkeb T. 2004. Fluid evolution in the W-Cu-Zn-Pb San Cristobal vein, Peru: Fluid inclusion and stable isotope evidence. Chem. Geol., 210: 201–224. DOI:10.1016/j.chemgeo.2004.06.008
[] Brown PE. 1989. FLINCOR: A microcomputer program for the reduction and investigation of fluid-inclusion data. Am. Mineral., 74: 1390–1393.
[] Bureau of Geology and Mineral Resources of Jiangxi Province. 1984. Regional Geology of Jiangxi Province. Beijing: Geological Publishing House : 1 -921.
[] Bussink RW, Kreulen R. 1984. Gas analyses fluid inclusion and stable isotopes of Panasqueira W-Sn deposits, Portugal. Bull. Mineral, 107: 703–713.
[] Cattalani S, Williams-Jones AE. 1991. C-O-H-N fluid evolution at Saint-Robert, Quebec: Implications for W-Bi-Ag mineral deposition. Canadian Mineralogist, 29: 435–452.
[] Chen X. 2007. Origin of ore-host fault system in Huangsha tungsten vein deposit, southern Jiangxi Province. Ph. D. Dissertation. Beijing: China University of Geosciences, 1-110(in Chinese with English summary)
[] Chen ZD, Hu LS. 1984. The geological characteristics and primary zoning of Huangsha vein type tungsten deposit. In: Yu HZ (ed.). Proceedings of Symposium on Tungsten Geology. Nanchang, 1981. Beijing: Geological Publishing House : 25 -34.
[] Chi GX, Lu HZ. 2008. Validation and representation of fluid inclusion microthermometric data using the fluid inclusion assemblage (FIA) concept. Acta Petrologica Sinica, 24(9): 1945–1953.
[] Clark AH, Kontak DJ, Farrar E. 1990. The San Judas Tadeo W (-Mo, Au) deposit: Permian lithophile mineralization in southeastern Peru. Economic Geology, 85: 1651–1668. DOI:10.2113/gsecongeo.85.7.1651
[] Fan HR, Hu FF, Yang JH, Shen K, Zhai MG. 2005. luid evolution and large-scale gold metallogeny during Mesozoic tectonic transition in the eastern Shandong Province. Acta Petrologica Sinica, 5(21): 1317–1328.
[] Feng ZW, Xia WH, Zhang JT, Chen ZY. 1989. Geological characteristics of the Huangsha vein-type tungsten deposit and a discussion on the property of its ore-forming fluid, Jiangxi. Earth Science-Journal of China University of Geoscience, 14(4): 423–431.
[] Gan GL. 1988. The characteristics of inclusions in minerals of the Huangsha vein-type tungsten deposit, Jiangxi. Journal of Guilin University of Technology, 8(2): 133–139.
[] Giamello M, Protano G, Riccobono F, et al. 1992. The W-Mo deposit of Perda Majori (SE Sardinia, Italy): A fluid inclusion study of ore and gangue minerals. Eur. J. Mineral., 4: 1079–1084. DOI:10.1127/ejm/4/5/1079
[] Giuliani G, Li YD, Sheng TF. 1988. Fluid inclusion study of Xihuashan tungsten deposit in the southern Jiangxi Province, China. Mineralium Deposita, 23: 24–33.
[] Goldstein RH, Reynolds TJ. 1994. Systematics of fluid inclusions in diagenetic minerals. SEPM Short Course, 31: 199.
[] Goldstein RH. 2003. Petrographic analysis of fluid inclusions. In: Samson I, Anderson A and Marshall D (eds.). Fluid Inclusions-Analysis and Interpretation. Mineralogical Association of Canada, Short Course Series, 32: 9–53.
[] Graupner T, Kempe U, Dombon E. 1999. Fluid regime and ore formation in the tungsten (-yttrium) deposits of Kyzyltau (Mongolian Altai) evidence for fluid variability in tungsten-tin ore systems. Chem. Geol., 154: 29–40.
[] Gu JY. 1984. Morphological zoning of the vein-type tungsten deposits in south China. In: Yu HZ (ed.). Proceedings of Symposium on Tungsten Geology. Nanchang, 1981. Beijing: Geological Publishing House : 35 -45.
[] Han JZ, Hu XM, Liu SL. 1984. The congtrols of strata over tungsten ore deposits as viewed from tungsten abundances in the strata of southern Jiangxi. Geochimica(2): 176–179.
[] Heinrich CA. 1990. The chemistry of hydrothermal tin (-tungsten) ore deposition. Economic Geology, 85(3): 457–481. DOI:10.2113/gsecongeo.85.3.457
[] Higgins NC. 1980. Fluid inclusion evidence for the transport of tungsten by carbonate complexes in hydrothermal solutions. Can. J. Earth Sci., 17: 823–830. DOI:10.1139/e80-082
[] Higgins NC, Kerrich R. 1982. Progressive 18O depletion during CO2 separation from a carbon dioxide-rich hydrothermal fluid: Evidence from the Grey River tungsten deposit, Newfoundland: Can. J. Earth Sci., 19: 2247–2257.
[] Hu DQ, Hua RM, Li GL, Wei XL, Huang XE. 2011. Study on the fluid inclusions of Maoping tungsten deposit, southern Jiangxi Province. Geological Journal of China Universities, 17(2): 327–336.
[] Hall DL, Sterner SM, Bodnar RJ. 1988. Freezing point depression of NaCl-KCl-H2O solutions. Economic Geology, 83: 197–202. DOI:10.2113/gsecongeo.83.1.197
[] Hu FF, Fan HR, Yang KF, Shen K, Zhai MG, Jin CW. 2007. Fluid inclusions in the Denggezhuang lode gold deposit at Muping, Jiaodong Peninsula. Acta Petrologica Sinica, 23(9): 2155–2164.
[] Jackson NJ, Willis-Richards J, Manning DAC, et al. 1989. Evolution of the Cornubian ore field, Southwest England: Part II. Mineral deposits and ore-forming processes. Economic Geology, 84: 1101–1133.
[] Kamilli RJ, Cole JC, Elliott JE. 1992. Geology and genesis of the Baid Al Jimalah tungsten deposit, Kingdom of Saudi Arabia. Economic Geology, 88: 1743–1767.
[] Kang YF, Li CY. 1991. Geological characteristics, types and distribution of tungsren deposits in China. Mineral Deposits, 10(1): 19–26.
[] Kelly WC, Rye RO. 1979. Geologic fluid inclusion and stable isotope studies of the tin-tungsten deposits of Panasqueira, Portual. Economic Geology, 74: 1721–1822. DOI:10.2113/gsecongeo.74.8.1721
[] Landis GP, Rye RO. 1974. Geologic, fluid inclusion and stable isotope studies of the Pasto Bueno tungsten base metal deposit, Northern Peru. Economic Geology, 69: 1025–1059. DOI:10.2113/gsecongeo.69.7.1025
[] Li GL. 2011. The evolution of Yanshanian Granite and Tungsten Mineralization in Southern Jiangxi Province and adjacent region. Ph. D. Dissertation. Nanjing: Nanjing University, 1-189 (in Chinese with English summary)
[] Liu RL, Mu JL, Xu R. 1985. A study of Dalongshan tungsten-molybdenum deposit, Jiangxi Province. Journal of Chengdu University of Technology(2): 25–31.
[] Liu RL, Mu JL. 1990. Research on fluid inclusions and its applicationg to the study of ore-forming stages in the Muziyuan W-Mo deposit, Jiangxi Province. Journal of Chengdu University of Technology, 17(3): 18–28.
[] Lu HZ, Fan HR, Ni P, Ou GX, Shen K, Zhang WX. 2004. Fluid Inclusions. Beijing: Science Press : 406 -419.
[] Lynch JVG. 1989. Hydrothermal alteration, veining, and fluid inclusion characteristics of the Kalzas wolframite deposit, Yukon: Can. J. Earth Sci., 26: 2106–2115.
[] Mao JW, Zhang ZC, Zhang ZH, Du AD. 1999. Re-Os isotopic dating of molybdenites in the Xiaoliugou W (Mo) deposit in the Northern Qilian Mountains and its geological significance. Geochimica et Cosmochimica Acta, 63(11-12): 1815–1818. DOI:10.1016/S0016-7037(99)00165-9
[] Naden J, Shepherd TJ. 1989. Role of methane and carbon dioxide in gold deposition. Nature, 342: 793–795. DOI:10.1038/342793a0
[] Noronha F, Doria A, Dubessy J. 1992. Characterization and timing of the different types of fluids present in the barren and ore-veins of the W-Sn deposit of Panasqueira. Central Portugal. Mineralium Deposita, 27: 72–79.
[] O’Reilly C, Gallagher V, Feely M. 1997. Fluid inclusion study of the Ballinglen W-Sn-sulphide mineralization, SE Ireland. Mineralium Deposita, 32: 569–580. DOI:10.1007/s001260050123
[] Polya DA. 1989. Chemistry of the main-stage ore-forming fluids of the Panasqueira W-Cu(Ag)-Sn deposit. Portugal: Implications for models of ore genesis. Economic Geology, 84: 1134–1152.
[] Polya DA. 1990. Pressure-dependence of wolframite solubility for hydrothermal vein formation. Transactions of the Institute of Mining and Metallurgy, 99: B120–B124.
[] Ramboz C, Pichavant M, Weisbrod A. 1982. Fluid immiscibility in natural processes: Use and misuse of fluid inclusion data. II. Chem. Geol., 37: 29–48. DOI:10.1016/0009-2541(82)90065-1
[] Ramboz C, Schnapper D, Dubessy J. 1985. The P-V-T-X-fO2 evolution of H2O-CO2-CH4-bearing fluids n a wolframite vein: Reconstruction from fluid inclusion studies. Geochim. Cosmochim. Acta, 49: 205–219. DOI:10.1016/0016-7037(85)90205-4
[] Rankin AH, Alderton DHM. 1983. Fluid Inclusion Petrography of SW England granites and its potential in mineral exploration. Mineralium Deposita, 18: 335–347.
[] Ren YC. 1998. Restudy of bismuth sulfosalt minerals in Pangushan-Huangsha wolfram-quartz vein deposit, Jiangxi Province. Contributions to Geology and Mineral Resources Research, 13(4): 1–17.
[] Rios FJ, Villas RN, Fuzikawa K, et al. 2003. Fluid evolution in the Pedra Preta wolframite ore deposit Paleoproterozoic Musa Granite, eastern Amazon craton, Brazil. South American Earth Sciences, 15: 787–802. DOI:10.1016/S0895-9811(02)00128-1
[] Roedder E. 1984. Fluid inclusions. Mineralogical Society of America, reviews in Mineralogy, 12: 644.
[] Samson IM. 1990. Fluid evolution and mineralization in a subvolcanic granite stock: The Mount Pleasant W-Mo-Sn deposits, New Brunswick, Canada. Economic Geology, 85: 145–163.
[] Seal RR II, Clark AH, Morrissey CJ. 1987. Stockwork tungsten (scheelite)-molybdenum mineralization, Lake George, Southwestern NewBrunswick. Economic Geology, 82: 1259–1282. DOI:10.2113/gsecongeo.82.5.1259
[] Shelton KL, Supso C, Rye DM, et al. 1986. Geologic, sulfur isotope, and fluid inclusion studies of the Sannae W-Mo mine, Republic of Korea: Comparison of sulfur isotope systematics in Korean W deposits. Economic Geology, 81: 430–446. DOI:10.2113/gsecongeo.81.2.430
[] Sheng JF. 2004. Metallogeny of the Major Metallic Ore Deposits in China. Beijing: Geological Publishing House : 146 .
[] Shepherd TJ, RaKin A, Alderton DHM. 1985. A Practical Guide to Fluid Inclusion Studies. Blackie and Son Limited: 1–154.
[] So CS, Yun ST. 1994. Origin and evolution of W-Mo-producing fluids in a granitic hydrothermal system: Geochemical studies of quartz vein deposits around the Susan granite, Hwanggangri district, Republic of Korea. Economic Geology, 89: 246–267. DOI:10.2113/gsecongeo.89.2.246
[] Thiéry R, van der Kerkhof AM, Dubessy F. 1994. VX properties of CH4-CO2 and CO2-N2 fluid inclusions: modeling for T<31℃ and P<400bar. European Journal of Mineralogy, 6: 753–771. DOI:10.1127/ejm/6/6/0753
[] Thorn PG. 1988. Fluid inclusion and stable isotope studies at the Chicote tungsten deposit, Bolivia. Economic Geology, 83: 62–68. DOI:10.2113/gsecongeo.83.1.62
[] Vallance J, Cathelineau M, Mariganc C, et al. 2001. Microfracturing and fluid mixing in granites: W-(Sn)ore deposition at Vaulry (NE French Massif Central). Tectonophysics, 336: 43–61. DOI:10.1016/S0040-1951(01)00093-2
[] van den Kerkhof A. 1990. Isochoric phase diagram in the system CO2-CH2 and CO2-N2: Application to fluid inclusions. Geochim. Cosmochim. Acta, 54: 621–629. DOI:10.1016/0016-7037(90)90358-R
[] Wang XD, Ni P, Jiang SY, Huang JB. 2008. Fluid inclusion study on the Piaotang tungsten deposit, southern Jiangxi Province, China. Acta Petrologica Sinica, 24(9): 2163–2170.
[] Wang XD, NI P, Jiang SY, et al. 2010. Origin of ore-forming fluids of the Piaotang tungsten deposit in Jiangxi Province: Evidence from Helium and argon isotopes. Chinese Science Bulletin, 55(7): 628–634. DOI:10.1007/s11434-009-0416-8
[] Wang XD, Ni P, Zhang PS, Wang TG. 2011. Fluid inclusion studies of the Pangushan quartz-vein type tungsten deposit in southern Jiangxi Province. Acta Petrologica Et Mineralogica, 29(5): 539–550.
[] Wood SA. 1992. Experimental determination of the solubility of W03(s) and the thermodynamic properties of H2WO4(aq) in the range 300~600℃ at 1kbar: Calculation of scheelite salability. Geochim. Cosmochim. Acta, 56: 1827–1836. DOI:10.1016/0016-7037(92)90313-8
[] Wood SA, Samson IM. 2000. The hydrothermal geochemistry of tungsten in granitoid environments: I. Relative solubilities of ferberite and scheelite as a function of T, P, pH and mNaCl. Economic Geology, 95: 143–182.
[] Xia HY, Liang SY, Shuai DQ, Xie WX. 1982. The existing form and enrichment regularity of Huangsha wolframite-quartz vein type deposit, Jiangxi Province. Chinese Science Bulletin(3): 164–167.
[] Xia HY, Liang SY, Xie WX, et al. 1984. The primary zoning and its genesis of Huangsha tungsten deposit, Jixangxi Province. In: Yu HZ (ed.). Proceedings of Symposium on Tungsten Geology. Nanchang, 1981. Beijing: Geological Publishing House : 221 -230.
[] Xia HY, Liang SY. 1987. The evolution of Huangsha-Tieshanlong ore-bearing granite and its REE geochemistry. Geochimica(4): 330–340.
[] Xia HY. 1989. The tungsten ore deposits in Huangsha, Jiangxi Province. In: Chen YC, Pei RF, Zhong HL et al. (eds.). The Geology of Nonferrous and Rare Metal Deposits Related to Mesozoic Granitoids in Nanling Region. Beijing: Geological Publishing House : 87 -101.
[] Xiang CP. 1985. Geological features and genesis of the Huangsha vein tungsten (copper) deposit, Jiangxi. Geological Review, 31(5): 422–428.
[] Xu KQ, Hu SX, Sun MZ, et al. 1984a. Regional factors controlling the formation of tungsten deposits in south China. In: Yu HZ (ed.). Proceedings of Symposium on Tungsten Geology. Nanchang. Beijing: Geological Publishing House : 245 -249.
[] Xu KQ, Sun N, Wang DZ, et al. 1984b. Genesis and metallogenesis of granites of South China. In: Xu KQ and Tu GC (eds.). Granite and Related of Metallogenesis (Collected Papers of International Symposium). Nanjing, 1982. Nanjing: Jiangsu Science and Technology Press: 1–20.
[] Yu CW. 2004. Fractal dilatation of multiple Hydraulic fracturing. Earth Science Frontiers, 11(1): 11–44.
[] Zaw UK. 1984. Geology and geothermometry of vein-type W-Sn deposits at Pennaichaung and Yetkanzintaung prospects, Tavoy Township, Tennasserim Division, Southern Burma. Mineralium Deposita, 19: 138–144.
[] Zhu YL, Li CY, Lin YH. 1981. Geology of Tungsten Deposit of Southern Jixangxi Province. Nanchang: Jiangxi People's Publishing House: 1–419.
[] 陈祥. 2007. 赣南黄沙脉状钨矿床容矿断裂构造系统成因探讨. 博士学位论文. 北京:中国地质大学,1-110
[] 陈尊达, 胡立檖. 1984. 1984黄沙脉钨矿床地质特征及原生分带. 北京: 地质出版社: 25-34.
[] 池国祥, 卢焕章.2008. 流体包裹体组合对测温数据有效性的制约及数据表达方法. 岩石学报, 24(9): 1945–1953.
[] 冯志文, 夏卫华, 章锦统, 陈紫英.1989. 江西黄沙脉钨矿床特征及成矿流体性质讨论. 地球科学—中国地质大学学报, 4(14): 423–431.
[] 范宏瑞, 胡芳芳, 杨进辉, 沈昆, 翟明国.2005. 胶东中生代构造体制转折过程中流体演化和金的大规模成矿. 岩石学报, 5(21): 1317–1328.
[] 干国梁.1988. 江西黄沙钨矿床的矿物流体包裹体特征. 桂林冶金地质学院学报, 8(2): 133–139.
[] 古菊云. 1984. 华南脉钨矿床的形态分带. 北京: 地质出版社: 35-45.
[] 韩久竹, 胡心铭, 刘世濂.1984. 从赣南地层含钨丰度看地层对钨矿的控制. 地球化学, 2: 176–179.
[] 胡东泉, 华仁民, 李光来, 韦星林, 黄小娥.2011. 赣南茅坪钨矿流体包裹体研究. 高校地质学报, 17(2): 327–336.
[] 胡芳芳, 范宏瑞, 杨奎锋, 沈昆, 翟明国, 金成伟.2007. 胶东牟平邓格庄金矿床流体包裹体研究. 岩石学报, 23(9): 2155–2164.
[] 江西省地质矿产局. 1984. 江西省区域地质志. 北京: 地质出版社: 1-921.
[] 康永孚, 李崇佑.1991. 中国钨矿地质特征、类型及其分布. 矿床地质, 10(1): 19–26.
[] 李光来. 2011. 赣南及邻区燕山期花岗岩演化与钨成矿作用. 博士学位论文. 南京:南京大学,1-189
[] 刘若兰, 慕纪录, 徐荣.1985. 江西大龙山钨钼矿床流体包裹体的研究. 成都地质学院学报(2): 23–31.
[] 刘若兰, 慕纪录.1990. 木梓园钨钼矿床流体包裹体研究及在成矿阶段划分中的应用. 成都地质学院学报, 17(3): 18–28.
[] 卢焕章, 范宏瑞, 倪培, 欧光习, 沈昆, 张文淮. 2004. 流体包裹体. 北京: 科学出版社: 406-419.
[] 任英忱.1998. 江西盘古山-黄沙黑钨矿石英脉矿床铋硫盐矿物再研究. 地质找矿论丛, 13(4): 1–17.
[] 盛继福. 2004. 钨矿床. 见: 赵一呜, 吴良士, 白鸽等编. 中国主要金属矿床成矿规律. 北京: 地质出版社: 1-146.
[] 王旭东, 倪培, 蒋少涌, 黄建宝.2008. 赣南漂塘钨矿流体包裹体研究. 岩石学报, 24(9): 2163–2170.
[] 王旭东, 倪培, 张伯声, 王天刚.2011. 江西盘古山石英脉型钨矿床流体包裹体研究. 岩石矿物学杂志, 29(5): 539–550.
[] 夏宏远, 梁书艺, 帅德权, 谢为鑫.1982. 江西黄沙黑钨矿石英脉中银的存在形式和富集规律. 科学通报(3): 164–167.
[] 夏宏远, 梁书艺, 谢为鑫, 等.1984. 江西黄沙钨矿床的原生分带及其成因. 北京: 地质出版社: 221–230.
[] 夏宏远, 梁书艺.1987. 黄沙-铁山垄含矿花岗岩演化和稀土元素地球化学. 地球化学(4): 330–340.
[] 夏宏远. 1989. 江西黄沙钨矿床.见:陈毓川, 裴荣富, 张宏良等编.南岭地区与中生代花岗岩有关的有色及稀有金属矿床地质. 北京: 地质出版社: 87-101.
[] 项尝培.1985. 江西黄沙脉状钨(铜)矿床地质特征及成因探讨. 地质论评, 31(5): 422–428.
[] 徐克勤, 胡受奚, 孙明志, 等. 1984a. 华南钨矿床的区域成矿条件分析. 北京: 地质出版社: 245-249.
[] 徐克勤, 孙鼐, 王德滋, 等. 1984b. 华南花岗岩成因与成矿. 南京: 江苏科学技术出版社: 1-20.
[] 於崇文.2004. 多重水力断裂的分形扩张. 地学前缘, 11(1): 11–44.
[] 朱焱龄, 李崇佑, 林运淮. 1981. 赣南钨矿地质. 南昌: 江西人民出版社: 1-419.