岩石学报  2019, Vol. 35 Issue (9): 2711-2726, doi: 10.18654/1000-0569/2019.09.06   PDF    
粤北仁化棉花坑铀矿床成矿热液演化及其对成矿过程的约束
祁家明1, 朱捌1, 吴建勇1, 曹豪杰1, 刘文泉1, 徐争启2     
1. 核工业二九〇研究所, 韶关 512029;
2. 成都理工大学地球科学学院, 成都 610059
摘要: 本文以棉花坑(302)铀矿床成矿期紫黑色萤石、浅粉红色方解石及赤红色微晶石英等含铀脉石矿物及与之共生的黄铁矿为研究对象,采用流体包裹体热力学、群体成分分析及黄铁矿微量元素分析等方法,对成矿流体演化特性及对成矿过程的指示与约束开展了研究。研究表明,成矿期的萤石、方解石、微晶石英中流体包裹体类型以富液相两相NaCl-H2O型为主,平均均一温度分别为185.8℃、177.0℃、140.4℃,平均盐度分别为2.24% NaCleqv、1.36% NaCleqv、1.75% NaCleqv,矿床流体具有中低温低盐度特征;计算出平均成矿压力分别为39.5MPa、38.0MPa、30.1MPa,平均成矿深度分别为1.5km、1.3km、1.1km。流体包裹体群体成分显示成矿流体中富含K+、Na+、Ca2+等阳离子和HCO3-、F-、SO42-等阴离子及CO2、H2O等气相成分。这些脉石矿物为成矿期不同阶段沉淀的产物,随着成矿流体温度、压力逐渐降低,流体存在演化分异和不混溶现象,流体内的∑M+/∑M-逐渐升高,矿物沉淀按先析出萤石、其次方解石、最后微晶石英的顺序进行。成矿期黄铁矿Y/Ho平均比值变化显示,矿物沉淀过程逐渐改变了成矿流体性质,使得Zr/Hf、Nb/La、Co/Ni等稀土、高场强元素平均比值逐渐变小,还原性的成矿环境也会发生轻微波动;铀元素在流体演化的最晚阶段才大量与微晶石英一同沉淀,沉淀出的黄铁矿U/Th平均比值逐渐升高。
关键词: 流体包裹体    群体成分    黄铁矿    微量元素    成矿流体演化    棉花坑铀矿    广东北部    
The evolution of ore-forming fluid and its constraint on mineralization process in Mianhuakeng uranium deposit, northern Guangdong, China
QI JiaMing1, ZHU Ba1, WU JianYong1, CAO HaoJie1, LIU WenQuan1, XU ZhengQi2     
1. Research Institute No. 290, China National Nuclear Corporation, Shaoguan 512029, China;
2. College of Earth Sciences, Chengdu University of Technology, Chengdu 610059, China
Abstract: In this paper,the evolution of ore-forming fluid and its constraint on ore-mineralization process in the Mianhuakeng (No.302) uranium deposit has been studied by means of thermodynamics and composition of fluid inclusions in various syn-ore gangue minerals (such as purple black fluorite,light pink calcite and red microcrystalline quartz) and the trace elements of symbiosis pyrite. The results show that the inclusions of syn-ore fluorite,calcite and microcrystalline quartz were mainly liquid-rich two-phase NaCl-H2O with the average homogeneous temperatures of 185.8℃,177.0℃ and 140.4℃ and average salinity of 2.24% NaCleqv,1.36% NaCleqv,1.75% NaCleqv respectively. That was to say,the No.302 deposit had the characteristics of hydrothermal mineralization at low temperature and low salinity. In addition,the average ore-forming pressure was 39.5MPa,38.0MPa and 30.1MPa,and the average ore-forming depth was 1.5km,1.3km and 1.1km respectively. The composition of fluid inclusion groups has been showed that the ore-forming fluids were rich in K+,Na+,Ca2+ cations,and HCO3-,F-,SO42- anions as well as gaseous components such as CO2 and H2O. It was considered that these gangues were the products of precipitation in different stages of the syn-ore period,there were evolutionary differentiation and immiscibility in the hydrothermal syn-ore process with the decreasing of temperature and pressure of the ore-forming fluid and increasing of ∑M+/∑M- in the inclusions. Mineral precipitation was accompanied by the sequence of gangue mineral formation of fluorite,calcite and microcrystalline quartz. The change of Y/Ho average ration of pyrite associated with syn-ore period showed that the properties of the ore-forming fluid has been changed during the mineral precipitation process,which made the average ratio of REE and HFSE such as Zr/Hf,Nb/La,Co/Ni gradually decrease,and made the reductive syn-ore environment undergo slight microwave movement,uranium element would be precipitated with microcrystalline quartz in large quantities at the latest stage of fluid evolution,the U/Th average ratio of the precipitated pyrite was increased gradually.
Key words: Fluid inclusions    Groups composition    Pyrite    Trace elements    Evolution of ore-forming fluids    Mianhuakeng uranium deposit    Northern Guangdong    

华南地区是我国的重要铀矿产区,按产出围岩的不同主要分为花岗岩型、火山岩型和碳硅泥岩型等三类热液铀矿床(Hu et al., 2008, 2009)。而粤北棉花坑(302)铀矿床是华南地区最具代表性的花岗岩型铀矿床之一,也是目前粤北地区勘探、开采深度最深、研究程度最高的铀矿床。前人对该矿床的年代学、地球化学、微区矿物学以及矿床成因等方面进行了详细研究(余达淦, 1979;杜乐天和王玉明, 1984;闵茂中等, 1993;胡瑞忠等, 2004;张国全等, 2008;沈渭洲等, 2010;黄国龙等, 2010;高飞等, 2011;祁家明等, 2014, 2015;钟福军等, 2017;张龙等, 2018);尽管一些研究从流体包裹体方面对其成矿流体性质与物质来源进行了诸多讨论(倪师军等, 1994;张国全等, 2008;郭国林等, 2010;张闯等, 2016; Zhang et al., 2017),但对成矿过程中成矿热液演化特征及其对元素迁移的约束作用研究还较为缺乏,成矿过程是多重参数动态耦合的过程,其深入研究对于揭示铀成矿演化机制和成矿动力学背景具有重要意义,可以提高华南热液型铀矿床的理论研究水平,指导区内铀矿勘查工作。这也是本文撰写的初衷。

棉花坑(302)铀矿床发育大量成矿阶段形成的各类热液矿物,包括大量萤石、方解石和微晶石英等脉石矿物,是研究成矿流体性质的重要对象。而流体包裹体是获取热液矿床成矿流体信息最直接、最有效的手段。研究表明,通过不同温度下流体包裹体相变行为的观察,可以了解成矿流体的温度、压力、密度及组成等重要参数,以探讨成矿流体的演化特征(卢焕章, 2011)。

同时矿床成矿期热液脉体内还发育大量黄铁矿,是铀成矿过程中形成的重要金属矿物,黄铁矿发育多的部位矿石品位也较高,是重要的指示富铀矿体的标志矿物(张国全等, 2007)。研究表明,黄铁矿中的微量元素主要以类质同象或晶格缺陷内包体形式存在,不同环境下形成的黄铁矿的微量元素含量主要受成矿流体性质制约(Moh, 1980;毕献武等, 2004; Reich et al., 2005);同时,稀土元素可以有效指示热液体系中成矿流体来源和水-岩作用历史(Henderson, 1984)。研究不同阶段黄铁矿的微量元素,亦可反映成矿流体的演化特征而约束成矿过程。因此,本文选择棉花坑(302)铀矿床成矿期重要脉石矿物中的流体包裹体及其共生矿物黄铁矿为研究对象,重点解析铀成矿过程中成矿热液性质和组成及其演化特征,探讨铀成矿机制。

1 区域地质背景

在21世纪初,华南地区的花岗岩型铀矿储量曾占全国探明铀矿储量的近40%,而粤北地区诸广-贵东岩体的铀资源量占华南全区的1/2以上,是我国最重要的花岗岩型铀矿集区(黄国龙等, 2010;陈振宇等, 2014)。诸广-贵东地区位于南岭EW向构造-岩浆带中段,区域上位于华南板块与扬子板块接触部位(图 1),以中生代花岗岩为主,区内铀矿床多赋存于花岗岩体内部或与其围岩接触带附近的断裂构造中。

图 1 粤北区域地质简图 1-加里东期花岗岩类;2-华力西期花岗岩类;3-印支期花岗岩类;4-燕山早期花岗岩类;5-燕山晚期花岗岩类;6-震旦系-志留系;7-泥盆系-三叠系;8-下-中侏罗统;9-上白垩统;10-古近系;11-区域断裂;12-铀矿床 Fig. 1 Regional geological map in northern Guangdong 1-Caledonian granitoids; 2-Variscan granitoids; 3-Indosinian granitoids; 4-Early Yanshanian granitoids; 5-Late Yanshanian granitoids; 6-Sinian-Silurian; 7-Devonian-Triassic; 8-Lower-Middle Jurassic; 9-Upper Cretaceous; 10-Paleogene; 11-regional faults; 12-uranium deposits

棉花坑(302)铀矿床所在的诸广岩体为多期多阶段复式岩体,侵入于寒武系和泥盆系地层中,岩浆活动频繁,但以印支期和燕山期S型黑云母花岗岩为主,地球化学属性为高硅、高铝、碱性-钙碱性系列,印支期成岩年龄为245~225Ma,燕山期成岩年龄为170~150Ma(邓平等, 2003)。其岩体含铀量可达13×10-6~24×10-6,是一般花岗岩的数倍,是区内重要的产铀岩体,为铀成矿提供了丰富的铀源(邓平等, 2003)。

区内构造活动强烈,形成了一系列断裂构造,晚侏罗世-早白垩世多为控岩断裂构造,而晚白垩世-新近纪形成的断裂多为控盆断裂构造。断裂展布方向有SN、EW、NW、NE等,这些断裂构造多为深大断裂,具有多次活动的特征。总体趋势是SN向断裂较早,EW向、NW向断裂次之,最晚为NE向断裂。深大断裂控制了区内中-新生代以来的岩浆活动和盆地形成与发展,同时也控制了铀矿和多金属矿产的形成。

区内铀矿多产于花岗岩区域内NE向主干断裂附近及与其伴生的二、三级硅化碎裂带中,铀矿类型以微晶石英型、碎裂蚀变带型和“交点”型为主,成矿温度集中于250~150℃,区内主要铀成矿呈多期多阶段性,其年龄集中于140Ma、120Ma、90Ma、70Ma和50Ma等阶段,成矿高峰期为90~70Ma,与花岗岩的年龄相差较大,而与区内的6次地壳伸展运动的时间则是一一对应的(张国全等, 2007)。可认为该些地壳伸展运动引起的深大断裂活动,是驱动区内铀成矿的动力学背景。

2 矿床地质

棉花坑(302)铀矿床位于诸广岩体东南部印支期中细粒二云母花岗岩(锆石U-Pb年龄232.0±4.0Ma; 黄国龙等, 2012)与燕山期中粒黑云母花岗岩(锆石U-Pb年龄157.2±1.7Ma; 黄国龙等, 2014)之接触部位,且由NE向棉花坑断裂与NW向油洞断裂所夹持的区域(图 2)。

图 2 棉花坑(302)铀矿床地质简图 1-碱交代岩;2-中粒黑云母花岗岩;3-中粗粒黑云母花岗岩;4-中细粒二云母花岗岩;5-中粗粒二云母;6-辉绿岩岩;7-断裂;8-铀矿床 Fig. 2 Geological sketch map of Mianhuakeng (No.302) uranium deposit 1-alkali metasomatic rock; 2-medium grain biotite granite; 3-medium-coarse grained biotite granite; 4-medium-fine grained two mica granite; 5-coarse-medium grain two mica granite; 6-diabase; 7-fault; 8-uranium deposit

铀矿体主要呈似脉状、扁豆状和透镜状产于上述夹持区域中的近SN向断裂蚀变破碎带内(图 3)。铀矿化垂幅大,自地表(海拔500m左右)至深部(海拔-647m)均见有工业铀矿体分布。矿石成分相对简单,矿石矿物以沥青铀矿为主,呈肾状、葡萄状、分散球粒状等形式产出,沥青铀矿形成温度介于250~150℃之间;少量次生铀矿物;铀成矿呈多期多阶段性,沥青铀矿年龄为~120Ma(张国全等, 2008;张龙等, 2018)、~102Ma(张龙等, 2018)、93±15Ma(Bonnetti et al., 2018)、~92Ma(张龙等, 2018)、70±11Ma(黄国龙等, 2010;沈渭洲等, 2010)、~68Ma(张龙等, 2018)、60.0±0.5Ma(Zhong et al., 2018)、~54Ma(张国全等, 2008)等系列成矿年龄。

图 3 棉花坑(302)铀矿床9号带剖面及样品采集示意图 1-细粒黑云母花岗岩;2-花岗斑岩;3-中粒二云母花岗岩;4-绢云母化碎裂花岗岩;5-硅化带;6-矿体;7-含矿构造蚀变带编号;8-采样点位 Fig. 3 Schematic map of samples collected from section of silicification zone No.9 in the Mianhuakeng (No.302) uranium deposit 1-fine-grained biotite granite; 2-granite porphyry; 3-medium-grained mica granite; 4-sericite cataclastic granite; 5-silicification zone; 6-uranium ore body; 7-ore-bearing structural alteration zone number; 8-samples site

矿体围岩蚀变发育,受断裂构造带控制,分布于构造破碎带及围岩蚀变带内,且具有明显的分带现象,距矿体由近及远可划分为6条带,即:矿化中心带(强赤铁矿化硅化带,褐红色微晶石英)、矿体旁强蚀变带(强绢云母、绿泥石化碎裂花岗岩)、矿体蚀变带(赤铁矿化、绢云母化黑云母花岗岩)、近矿端蚀变带(赤铁矿化黑云母花岗岩)、远矿端蚀变带(弱赤铁矿化黑云母花岗岩)和正常花岗岩带。

与铀矿化相关的伴生脉石矿物有褐红色微晶石英(图 4a, d-f)、紫黑色萤石(图 4b)、浅粉红色方解石(图 4c)和少量绿泥石、绢云母等。矿物组合通常为紫黑色萤石-黄铁矿-沥青铀矿、粉红色方解石-黄铁矿-沥青铀矿和赤红色微晶石英-赤铁矿-黄铁矿-沥青铀矿。与沥青铀矿伴生的金属矿物主要是胶状黄铁矿(图 4g, h),其次为赤铁矿(图 4g)、少量褐铁矿(图 4g)和方铅矿等。

图 4 棉花坑(302)铀矿床成矿期各类脉石特征照片 (a)赤红色微晶石英内沥青铀矿与胶状黄铁矿共生;(b)紫黑色萤石与灰色石英共生,萤石内含少量胶状黄铁矿;(c)粉红色方解石裂隙内含沥青铀矿脉,含少量胶状黄铁矿;(d)紫黑色萤石、方解石与微晶石英共生,内含胶状黄铁矿;(e)方解石与赤红色微晶石英共生,含少量胶状黄铁矿;(f)微晶石英与方解石共生;(g)样品F150内黄铁矿、褐铁矿及赤铁矿共生,反射光(-);(h)样品F41-3内浸染状黄铁矿,反射光(-). Py-黄铁矿;Lim-褐铁矿;Hem-赤铁矿;U-沥青铀矿;F-萤石;Si-微晶石英;Cal-方解石 Fig. 4 Photographs of metallogenic epoch gangues in the Mianhuakeng (No.302) uranium deposit (a) symbiosis of pitchblende and colloidal pyrite in red microcrystalline quartz; (b) purple-black fluorite coexists with grey quartz, and fluorite contains a small amount of colloidal pyrite; (c) pink calcite fissures contain pitch-uranium veins with a small amount of colloidal pyrite; (d) purple-black fluorite, calcite and microcrystalline quartz coexist, containing colloidal pyrite; (e) calcite coexists with red microcrystalline quartz and contains a small amount of colloidal pyrite; (f) symbiosis of microcrystalline quartz and calcite; (g) pyrite, limonite and hematite coexist in sample F150, reflecting light (-); (h) impregnated pyrite in sample F41-3, reflecting light (-). Py-pyrite; Lim-limonite; Hem-hematite; U-pitchblende; F-fluorite; Si-micro quartz; Cal-calcite

根据棉花坑(302)铀矿床的矿石矿物共生组合特征及其穿插关系,其热液活动(脉体)可划分为成矿前、成矿期和成矿后3个阶段,矿石矿物生成大致归纳为(表 1)。

表 1 棉花坑(302)矿床矿石矿物生成顺序表 Table 1 Mineral arisen sequence of the Mianhuakeng (No.302) uranium deposit

(1) 成矿前阶段:为白色高温石英脉(Ⅰ)形成阶段,为铀成矿热液活动的前奏。

主要以白色高温石英脉为主,偶见少量黄铁矿、黄铜矿,辉铋矿、闪锌矿等金属硫化物,脉体两侧围岩发育云英岩化及少量绿泥石化,该阶段的石英脉体已破碎成角砾,后被成矿阶段微晶石英脉所胶结。

(2) 成矿阶段:为铀成矿热液活动的主阶段,分为紫黑色含铀萤石脉(Ⅱ)、粉红色含铀方解石脉(Ⅲ)和赤红色含铀微晶石英(Ⅳ)等三个亚阶段,沥青铀矿主要是在赤红色微晶石英形成的后阶段析出的。

紫黑色含铀萤石脉多呈脉状、不规则块状或角砾状产出,脉宽仅数厘米,偶见少量胶状黄铁矿,部分萤石脉角砾被微晶石英或方解石胶结(图 4b),紫黑色萤石分布于矿床中、下部范围内,其铀含量较高,而矿床中、上部范围发育紫色、浅紫色萤石,其铀含量较低。

粉红色含铀方解石脉多与紫黑色萤石共生,呈脉状或不规则状,被微晶石英或灰白色方解石胶结,脉宽数厘米,内含数毫米宽的沥青铀矿脉(图 4c),偶见少量胶状黄铁矿。

赤红色含铀微晶石英脉是矿床最重要的含铀热液脉体,脉体组成矿物以隐晶、微晶石英和沥青铀矿为主(图 4a, d-f),其次有少量水(绢)云母、绿泥石、赤铁矿、黄铁矿、方铅矿和黄铜矿,偶尔可见少量的钠长石、萤石、高岭石和碳酸盐矿物。部分赤红色微晶石英呈角砾状,被后期灰白色微晶石英胶结。

(3) 成矿后阶段:为灰白色微晶石英、浅色萤石、灰白色方解石和石英晶洞(Ⅴ)阶段,是热液活动的尾期。

灰白色微晶石英-浅色萤石脉规模大小不等,宽度为数厘米至数米,主要由浅色萤石、灰白色微晶石英和一些粘土矿物(以高岭石为主)组成。脉体中的石英、萤石和粘土矿物常构成韵律条带。胶结成矿期形成的脉石角砾。

灰白色方解石脉呈细脉状和不规则脉体穿插到成矿阶段的热液脉体中,个别大脉可长达几米。

3 样品采集及测试方法

在棉花坑(302)铀矿床-150m中段以下共采集与成矿相关的脉石矿物6件,主要为微晶石英、方解石、萤石,在空间上均与铀矿紧密共生,样品描述见表 2,采样分布示意见图 3。样品经过手工挑选,磨制成双面抛光薄片14片,然后分别进行岩相学观察、显微测温,由于部分微晶石英薄片在显微镜下较难寻找到包裹体,因而能找到包裹体的薄片只10片。对含黄铁矿的样品经粉碎、过筛,在双目镜下挑选40~60目、纯度大于99%的黄铁矿单矿物。

表 2 棉花坑(302)铀矿床含矿脉石样品特征 Table 2 Characteristics of metallogenic epoch gangue samples in the Mianhuakeng (No.302) uranium deposit

包裹体岩相学观察和显微测温均在长江大学包裹体实验室完成,其中显微测温使用仪器为Lin Kam THMS600型冷热台,理论上可达到实验温度范围为-196~+600℃,温度显示0.1℃,控制稳定温度±0.1℃,光控直径1.13mm,样品轴向移动16mm,加热/冷冻速率0.01~130℃/min。仪器标定采用标准物质为KNO3、K2CrO3、CCl4及人工配制的NaCl标准溶液。当实验测度低于31℃时,误差为±0.2℃;介于31~300℃之间时,误差为±1℃;高于300℃时,误差为±2℃。包裹体盐度计算采用Hall et al. (1988)公式,密度计算采用刘斌和沈昆(1999)经验公式。分析结果见表 3

表 3 棉花坑(302)铀矿床成矿期各类脉石流体包裹体特征表 Table 3 Characteristics of fluid inclusions from syn-ore gangues of the Mianhuakeng (No.302) uranium deposit

包裹体的群体成分测试采用Shimadzu HICSP离子色谱仪和GC2010气相色谱仪,氦气作为载气,包裹体爆裂采用SGE公司热爆裂炉,其爆裂取样温度范围在100~500℃之间,阳离子最低检出限为×10-6,阴离子最低检出限位×10-9,分析结果见表 4

表 4 棉花坑(302)铀矿床成矿期各类脉石流体包裹体气相、液相成分特征表 Table 4 Characteristics of gas and liquid composition of liquid inclusions of various syn-ore gangues in the Mianhuakeng (No.302) uranium deposit

黄铁矿单矿物样品送至核工业北京地质研究院进行测试。对黄铁矿溶解后再ICP-MS上采用在线内标(Rh)法进行测试,测试仪器为PerkinElmer Elan DCR-e等离子质谱仪,主要实验条件为:仪器功率1400W,冷却气16L/min,进样冲洗时间50s,单个元素积分时间0.001s。测试依据标准为GB/T 14506《硅酸盐岩石化学分析方法》第29、30部分。各元素的检出限为0.002×10-6。分析结果见表 5

表 5 棉花坑(302)铀矿床成矿期各类脉石共生黄铁矿微量元素含量(×10-6) Table 5 Contents (×10-6) of trace elements in pyrite associated with various syn-ore gangues in the Mianhuakeng (No.302) uranium deposit
4 流体包裹体特征 4.1 包裹体岩相学特征

棉花坑(302)铀矿床成矿期脉体较为丰富,微晶石英、紫黑色萤石、浅肉红色方解石脉岩均发育大量包裹体,以原生为主,次生包裹体很少(图 5)。包裹体气液比变化较大,从10%~55%均有产出,但以10%~25%之间为主;长轴直径分布范围为5~68μm,但绝大部分分布在5~20μm之间;大部分包裹体形状以椭圆形、四边形为主,少数呈不规则状、长条状等。

图 5 棉花坑(302)铀矿床成矿期各脉石包裹体镜下照片 (a)方解石包裹体;(b)微晶石英包裹体;(c、d)萤石包裹体 Fig. 5 Microphotographs of fluid inclusions of various syn-ore gangues in the Mianhuakeng (No.302) uranium deposit (a) inclusions in calcite; (b) inclusions in microcrystalline quartz; (c, d) inclusions in fluorite

根据室温条件下流体包裹体的物理相态和成分,矿床内原生包裹体主要为H2O-NaCl包裹体,该类型亦可细分为以下两亚型:(1)富气相L+V两相水溶液包裹体(Ⅰ型);(2)富液相L+V两相水溶液包裹体(Ⅱ型)。其中L主要为水溶液,V主要为水蒸汽。各脉体的包裹体类型均较为单一,体系类型为NaCl-H2O,以Ⅱ型气液两相包裹体为主,只有方解石中含少量Ⅰ型两相气液包裹体。

萤石中全发育Ⅱ型气液两相包裹体,其气液比为10%~40%,包裹体一般呈椭圆形、四边形和不规则状,呈小群体分布或自由状态分布,包裹体大小为4~42μm。

方解石Ⅱ型包裹体气液比一般为10%~35%,包裹体主要呈椭圆形、不规则状和四边形,呈群状分布或自由状态分布,包裹体大小为4~38μm。

微晶石英Ⅱ型包裹体气液比为3%~40%,包裹体一般呈椭圆形、四边形或不规则状,呈群体分布或自由状态分布,包裹体大小为4~62μm。

4.2 包裹体显微热力学特征

对棉花坑(302)铀矿床各成矿期脉石矿物的流体包裹体进行冰点温度、均一温度等测试,数据见表 3,统计结果见图 6

图 6 棉花坑(302)铀矿床成矿期各脉石包裹体测温特征直方图 Fig. 6 Histograms of fluid inclusions temperature measurement of various syn-ore gangues in the Mianhuakeng (No.302) uranium deposit

从包裹体测温分析来看,紫黑色萤石的包裹体冰点温度为-3.5~-0.3℃(平均值-1.3℃);均一温度为92.3~302.1℃(平均值185.8℃)。计算得到的盐度为0.53%~5.74% NaCleqv(平均值2.24% NaCleqv),密度为0.727~0.973g/cm3(平均值0.892g/cm3)。

浅粉红色方解石主要以Ⅱ型包裹体为主,并含少量Ⅰ型包裹体。其中Ⅱ型包裹体的冰点温度为-2.3~-0.2℃(平均值-0.8℃);均一温度为65.4~254.7℃(平均值177.0℃)。计算得到的盐度为0.35%~3.88% NaCleqv(平均值1.36% NaCleqv);密度为0.779~1.003g/cm3(平均值0.902g/cm3)。

微晶石英包裹体冰点温度为-3.3~-0.3℃(平均值-1.0℃);均一温度为84.7~298.8℃(平均值140.4℃);计算得到的盐度为0.53%~5.43% NaCleqv(平均值1.75% NaCleqv);密度为0.747~0.993g/cm3(平均值0.937g/cm3)。

总体来看(图 6),萤石、方解石、微晶石英等三类脉石矿物包裹体测温参数的峰值分布非常紧密,同时在空间关系上紧密共生,表明它们应该是成矿期不同阶段沉淀的产物。此外,萤石、方解石、微晶石英等脉石包裹体平均的均一温度分别为185.8℃、177.0℃、140.4℃,平均盐度则为2.24% NaCleqv、1.36% NaCleqv、1.75% NaCleqv,表明该矿床为中低温热液成矿,这与前人的研究一致(张国全等, 2008;郭国林等, 2010;张闯等, 2016);同时,在成矿演化过程中流体温度、盐度在逐渐降低。

4.3 包裹体群成分特征

本文有代表性的分别选取了成矿期的萤石、方解石、微晶石英流体包裹体进行了气液成分分析(表 4),三类脉石矿物包裹体液相阳离子均以K+、Na+、Ca2+为主,Mg2+次之,液相阴离子以HCO3-、F-、SO42-为主,Cl-次之,成分上具有相似性。但三类包裹体内离子含量比值又有所不同,如萤石、方解石、微晶石英包裹体液相总阴阳离子比∑M+/∑M-分别为0.592、0.882、1.650,呈逐渐增大的趋势;而液相∑U含量分别为0.0034mol/L、0.0015mol/L、0.0005mol/L,呈逐渐降低的趋势。

三类包裹体气相成分主要为CO2、H2O,萤石、方解石包裹体内还含有少量还原性气体CH4,而H2含量则均低于检出限。F-、HCO3-是成矿流体中重要的矿化剂,当流体pH < 5时,氧化态的铀(U6+)与F-组成铀酰氟化物络合物主导形式运移,当pH更高时,氧化态的铀(U6+)与HCO3-组成铀碳酸盐络合物主导形式运移(凌洪飞, 2011),不同类型包裹体中CO2的含量变化基本可以代表成矿流体中CO2的含量变化,成矿期萤石、方解石、微晶石英流体包裹体CO2/H2O比值分别为0.0159、0.0198、0.0018,微晶石英包裹体CO2/H2O比值最低,液相总阴阳离子比∑M+/∑M-最高,液相∑U含量最低,均一温度最低,存在成矿流体的演化分异现象,表明该期成矿流体发生过沸腾作用或不混溶现象,导致CO2发生逃逸。

5 黄铁矿地球化学特征 5.1 微量元素特征

表 5可以看出,成矿期各类脉石的黄铁矿中异常富集U、Pb等大离子亲石元素,其中U含量为88.60×10-6~854.0×10-6,Pb含量为498.0×10-6~1148×10-6;微量元素标准化配分曲线整体向右倾斜(图 7),除U、Pb异常富集外,其余微量元素均出现不同程度的亏损,这可能与成矿流体性质有关,U、Pb由于离子半径不同很难以类质同象方式代替黄铁矿中的Fe,更有可能以包体或晶格缺陷方式存在于黄铁矿中,受成矿流体成分的不同而发生变化(周家云等, 2008)。

图 7 棉花坑(302)铀矿床成矿期脉石矿物内黄铁矿原始地幔标准化微量元素蛛网图(标准化值据Sun and McDonough, 1989) Fig. 7 Primitive mantle-normalized trace element spider diagrams of pyrite associated with various syn-ore gangues in the Mianhuakeng (No.302) uranium deposit (normalization values after Sun and McDonough, 1989)

不同样品的Cs、Rb、Pb、Ba、Th、U等大离子亲石元素以及Cu、Sc、Co、Cr、Ni等过渡族元素的分布特征大致相同,但Nb、Zr、Hf等高场强元素含量差异明显,在曲线上表现出明显的分异特征,可能是黄铁矿在随脉石矿物沉淀析出的过程中其成分发生了演化分异效应。

5.2 稀土元素特征

表 5可看出,三类脉石矿物内黄铁矿的稀土元素含量较低,∑REE含量为3.07×10-6~8.73×10-6,LREE含量为2.54×10-6~6.53×10-6,HREE含量为0.53×10-6~2.20×10-6。从稀土元素标准化配分图(图 8)可见,三类脉石矿物的黄铁矿均稍微富集轻稀土,稍亏损重稀土,Eu无较明显的负异常,配分曲线略呈右倾斜。δEu值为0.70~0.75,δCe值为1.02~1.47,说明该些黄铁矿的形成具有大致相同的环境。

图 8 棉花坑(302)铀矿床成矿期脉石矿物内黄铁矿球粒陨石标准化稀土元素配分曲线(标准化值据Sun and McDonough, 1989) Fig. 8 Chondrite-normalized REE patterns of pyrite associated with various syn-ore gangues in the Mianhuakeng (No.302) uranium deposit (normalization values after Sun and McDonough, 1989)

同时,三类脉石矿物的黄铁矿稀土配分模式又有所差别,表现于各黄铁矿稀土元素含量的系统差别(图 8),如与粉红色方解石内共生的黄铁矿稀土元素含量最高,与赤红色微晶石英共生的黄铁矿稀土元素含量最低,而与紫黑色萤石共生的黄铁矿各稀土元素含量位于二者之间,表现出黄铁矿同各自脉石矿物在沉淀析出过程中的热液演化差异。

6 讨论 6.1 成矿压力及深度演化

根据流体包裹体均一温度和盐度,利用公式(邵洁涟, 1988)可计算成矿流体压力:P=P0Th/T0(式中P0=219+2620ωT0=374+920ωP(105Pa)为成矿压力,Th为均一温度,ω为成矿流体盐度);同时,根据包裹体流体压力,利用适用于中低温流体成矿深度经验公式(张长青等, 2007;薛伟等, 2010)可计算对于的成矿深度:H=P/270(km)。

经计算,棉花坑矿床成矿期紫色萤石为代表的成矿压力为20.1~61.5MPa,平均值为39.4MPa,成矿深度为0.8~2.3km,平均值为1.5km;以成矿期粉红色方解石代表的成矿压力为14.4~52.8MPa,平均值为38.0MPa,成矿深度为0.5~2.0km,平均值为1.4km;以成矿期赤红色微晶石英为代表的成矿压力为18.5~61.2MPa,平均值为30.2MPa,成矿深度为0.7~2.3km,平均值为1.1km。

可见,棉花坑矿床成矿压力平均值为37.0MPa,成矿深度平均值为1.4km。同时,以成矿期萤石、方解石、微晶石英代表的成矿压力、成矿深度总体上呈递减趋势,表现出以不同脉石标型矿物为代表的成矿流体在成矿过程中的演化分异现象。

6.2 成矿流体演化分异效应

成矿流体的温度、盐度呈降低的趋势,流体各参数的改变会发生沸腾作用,流体的沸腾作用是引起矿床成矿物质沉淀的重要机制(Cox et al., 2001; Wilkinson, 2001; Hagemann and Lüders, 2003;陈衍景等, 2007;倪培等, 2005, 2018;邓明国等, 2018),会造成部分气液分离,该过程中包裹体捕获不同密度与盐度流体成分,当矿物晶体捕获流体后,捕获的流体盐度随温度降低而保持不变(侯明兰等, 2007;池国祥和卢焕章, 2008;池国祥和赖健清, 2009)。从棉花坑矿床成矿期各脉体Th-ωTh-ρ图解(图 9a, b)可以看出,若成矿流体中同时包含石英、方解石、萤石等三种脉石矿物的成分,在沿裂隙向上运移过程中,流体均一温度逐渐降低,流体的总体盐度是逐渐降低的,而密度则逐渐升高,呈现出流体的演化特性,但当包裹体缺陷捕获流体后,其均一温度则不再发生变化;铀矿在沉淀过程中流体发生∑CO2沸腾去气作用(商朋强等, 2006;张国全等, 2007),沸腾状态下捕获的包裹体,其均一压力等于捕获压力,可代表矿物形成时的压力(Roedder and Bodnar, 1980),因此,流体按形成温度高低依次析出萤石、方解石和微晶石英。在某种程度上,热液流体的矿物沉淀析出顺序与反映岩浆流体的“鲍文反应序列”是类似的。

图 9 棉花坑(302)铀矿床成矿期脉石矿物Th-ω (a)和Th-ρ (b)演化趋势图 Fig. 9 Th vs. ω (a) and Th vs. ρ (b) evolutionary trends of fluid inclusions in various syn-ore gangues in the Mianhuakeng (No.302) uranium deposit

同时,成矿时部分流体发生不混溶分异演化现象。某一种体系中具有两种或更多流体相互组合而共存的现象,称为流体的不混溶性(卢焕章, 2011)。一般可将流体相的不混溶性分为三类(Hollister and Crawford, 1981;赵甫峰等, 2011):(1)以水流体为主的沸腾;(2)从富水流体中分离出来的富CO2流体;(3)从某种硅酸盐熔体中分离出来的低密度流体。从对棉花坑矿床成矿流体包裹体测温数据来看,石英或方解石Ⅱ型包裹体中存在盐度、密度不同的流体,这说明在捕获时发生流体沸腾而呈现不混溶状态(宋玉财等, 2006;卢焕章, 2011)。从棉花坑矿床各脉石包裹体Th-ω(图 9a)来看,中低温成矿流体沸腾作用并不强烈,但随着温度的逐渐下降,流体不混溶分异还是存在一定的差异,即表现出一种为低盐度、低密度流体,另一种为较高盐度、较高密度流体,这是流体不混溶分异演化的直接证据。

至于为什么沥青铀矿主要产于微晶石英中,可能与铀元素的地球化学性质和成矿流体因发生不混溶而析出矿物的顺序有关,由上所知,萤石、方解石、微晶石英代表的平均成矿压力分别为39.5MPa、38.0MPa、30.1MPa,成矿流体温度、压力也逐渐降低,成分也相应发生变化,存在流体演化分异效应(魏浩等, 2011;王琳琳等, 2015)。成矿流体以含铀的碳酸盐-氟化物-硫化物-硅酸盐热液的形式存在,当流体250℃时溶解铀的能力最强,之后适度的CO2逸出有利于铀矿的沉淀,大量的CO2逸出反而会使沉淀的铀溶解(刘正义和刘红旭, 2009)。当成矿流体温度逐渐降低至300~210℃时,铀酰氟化物络合物分解,萤石逐渐沉淀形成,其晶体缺陷就能捕获该温度区间的相对高温、相对高盐度流体,以及捕获沸腾后的低温低盐度流体,由于该温度区间成矿流体具有最强的铀溶解能力,因而沉淀出的铀矿较少;当成矿流体温度降至210~180℃时,由于压力的降低,铀酰碳酸盐络合物分解,成矿流体开始发生不混溶现象逸出CO2,方解石逐渐形成,其晶体缺陷开始捕获该温度区间的流体及沸腾后的较低温流体,这一过程中由于CO2逸出较完全,沉淀的铀大部分会溶解重新进入流体(刘正义和刘红旭, 2009);由于U属于高度不相容大离子亲石元素,前期并不随流体的分异而大量进入析出的矿物中,而是在流体分异的最晚阶段,才与矿物一同沉淀(凌洪飞, 2011),因此,当流体温度降至180~160℃时,微晶石英与铀开始析出,其缺陷开始捕获该温度区间及低于该温度的流体。不同类型的矿物析出序列在图 9中得到很好的展现,同时也很好的解释了沥青铀矿主要赋存于赤红色微晶石英脉中的原因。

同时,热液中的大量存在的S2-、Fe2+是铀沉淀极其重要的还原剂,该些强还原性离子会与铀酰络合物破裂分解出的U6+发生氧化还原反应(张国全等, 2007;凌洪飞, 2011),U6+被还原成沥青铀矿(UO2),而S2-、Fe2+等则被氧化成胶状黄铁矿(FeS2)和浸染状赤铁矿(Fe2O3),所以,黄铁矿与铀矿共生,是热液成因产物,而不是先存在于围岩中提供还原环境的。

萤石、方解石、微晶石英等脉石矿物包裹体内∑U含量呈逐渐降低的趋势(表 4),成矿热液由于温度压力的逐渐降低,流体中的铀会达到饱和而逐渐沉淀,因而成矿流体内的铀含量是逐渐降低的,而矿石的铀矿品位则是逐渐升高的。

6.3 成矿过程中微量元素的演化与约束

由于稀土元素(REE)和高场强元素(HFSE)主要存在于黄铁矿内包裹体或晶体缺陷中,受黄铁矿晶体结构影响不大,而主要受形成黄铁矿的成矿流体介质特征控制(毕献武等, 2004),因此,棉花坑铀矿床黄铁矿的REE和HFSE特征可反应其成矿流体的某些特征。从成矿期与萤石、方解石、微晶石英等三类脉石矿物共生黄铁矿微量、稀土配分模式图(图 7图 8)来看,其配分模式及形态基本相同,表明它们具有相同的物质来源或相似的形成环境,但随着流体的分异演化,与矿物同时沉淀析出的微量元素也会发生分异演化现象,表现出各种元素含量或比值的逐渐演化。

富F的热液可同时有效富集HFSE、LREE,Zr/Hf、Nb/La等比值一般大于1(Oreskes and Einaudi, 1990;毕献武等, 2004);在热液富F、HCO3-络合物时,Y/Ho < 28(Bau and Dulski, 1999);热液成因黄铁矿Co/Ni值介于1~5之间(Monteiro et al., 2008;单强等, 2009)。棉花坑矿床与萤石、方解石、微晶石英等三类脉石矿物共生黄铁矿Zr/Hf(28.8~40.4)、Nb/La(1.02~1.49)、Y/Ho(21.54~40.0)、Co/Ni(1.17~1.71)、LREE/HREE(2.97~5.34)等均大于1,且各数值变化范围较小,说明成矿过程中并未受到其它热液活动的干扰。与萤石、方解石共生黄铁矿Y/Ho比值均小于28、而与微晶石英共生黄铁矿Y/Ho平均值为38.8,这表明成矿流体随着温度降低,在大量析出萤石、方解石后,流体的成分与性质发生了相应改变,使得Y/Ho比值发生变化,而Zr/Hf、Nb/La、Co/Ni平均比值逐渐变小;U/Th平均比值分别为175、422、609,呈现出逐渐增大的趋势(表 5)。

从成矿环境演化方面来看,δEu和δCe可作为氧化还原环境的指示标志(涂光炽, 1998),Eu2+易于在高温、还原性质的热液中存在,导致Eu正异常,但不易存在于低温还原的热液中;氧化条件下,Ce3+氧化为Ce4+,与其它元素分离,导致Ce异常(Deng et al., 2015;陈炳翰等, 2014;郭林楠等, 2019)。棉花坑矿床成矿期与萤石、方解石、微晶石英共生黄铁矿的平均δEu(0.75、0.71、0.72) < 1,平均δCe(1.16、1.05、1.39)>1 (表 5),说明其形成于中低温还原性的环境中,这与包裹体研究的结论相同;同时,在三类脉石矿物按生成温度依次沉淀演化过程中δEu、δCe值也在发生相应微小变化,结合包裹体测温,成矿流体温度降至300~210℃,F-析出形成萤石,此时δEu=0.75、δCe=1.16,流体为中温还原环境;当流体温度降至210~180℃,HCO3-析出方解石,此时δEu=0.71、δCe=1.05,流体中因HCO3-的大量减少呈低温较强还原环境;当流体温度降至180~160℃时,可能由于地下水的逐渐混合,此时δEu=0.72、δCe=1.39,流体呈低温弱还原环境。

6.4 铀成矿作用过程

在华南地区印支-燕山期复式花岗岩体形成之后,自中生代以来,该区共经历了6次构造伸展运动:(1)与晚侏罗-早白垩世岩浆断陷盆岭阶段:145~135Ma、120~115Ma和110~100Ma,(2)与晚白垩世-古近纪的陆内伸展断陷沉积盆地阶段:95~85Ma、75~70Ma和55~45Ma(胡瑞忠等, 2004;张国全等, 2007)。同时,粤北地区的辉绿岩Ar-Ar年龄主要分为~140Ma、~105Ma、~90Ma等3个阶段(李献华等, 1997),也与区内的该6次地壳伸展运动相吻合。而区内铀矿的主成矿期为95~65Ma,以95Ma陆内伸展断陷沉积活动为华南铀成矿大爆发的标志。因此,棉花坑铀矿床巨大的成岩、成矿时差,远大于岩浆活动的余热10Myr左右的最大时限,铀矿成因几乎不可能来自岩浆热液分异;而棉花坑(302)铀矿床的一系列成矿年龄看似杂乱无章,其实是与华南的6次地壳伸展运动一一对应(胡瑞忠等, 2004),非常符合华南断陷伸展活动年龄制约(表 6)。

表 6 棉花坑铀矿年龄与华南区域构造运动关系 Table 6 Relationship between uranium age of Mianhuakeng (No.302) uranium deposit and tectonic movement in South China

事实上,断陷红盆边缘控盆深大断裂、产铀花岗岩的分布明显控制了铀矿田的空间定位(林锦荣等, 2016)。因此,从华南铀成矿动力学背景出发,可对棉花坑(302)铀矿床及华南地区的其它花岗岩型铀矿床的成矿作用过程做如下厘定:

在区内印支-燕山产铀花岗岩体形成后,区内又经历了3期岩浆断陷盆岭伸展运动,使得区内花岗岩体的铀得以预活化,并发育宽广的绢云母化(胡瑞忠等, 2004)。在95Ma左右开始的持续3期的陆内断陷盆地伸展运动,华南地区由NW-SE向挤压环境转变为拉张环境,沿NE向断裂带张裂形成断陷盆地,其边缘的深大断裂带切穿地壳深部,此时少量基性岩浆上涌并快速侵位形成辉绿岩脉。之后深大断裂带的持续活动导致地壳深部热液上涌,但此时的热液氧逸度低(凌洪飞, 2011),在高温高压下沿断裂向上运移的过程中,不断的与岩体内下渗的地热裂隙水(大气降水)不断混合,伴随着∑CO2的加入(胡瑞忠等, 2004),形成高氧逸度、高萃取性的临界流体(凌洪飞, 2011),与所经岩体发生水-岩反应,使岩体内以晶质铀矿形式存在的预活化的U4+被氧化成U6+而形成成矿流体,并以铀酰络离子F-、UO2(CO3)n2(1-n)的形式迁移(石少华等, 2010)。随着成矿流体的不断上升,流体温度、压力的变化、氧化还原电位的变化、氧逸度的降低等多种因素耦合而发生铀的沉淀,沉淀过程中成矿热液会发生演化分异现象,按不同温度析出相应脉石矿物,而热液中的S2-、Fe2+等强还原性离子会与铀酰络合物破裂分解出的U6+发生氧化还原反应形成沥青铀矿。

7 结论

(1) 流体包裹体热力学研究表明,棉花坑铀矿床为中低温低盐度热液成矿。均一温度、盐度、成矿压力及成矿深度的变化显示,铀成矿过程存在热液的演化分异和不混溶现象,铀矿沉淀伴随先析出萤石、其次方解石、最后微晶石英的脉石矿物生成顺序,受铀元素的高度不相容性和成矿流体性质的制约,铀在最晚阶段才大量与微晶石英一同沉淀。

(2) 流体包裹体群体成分研究表明,成矿流体中富含K+、Na+、Ca2+等阳离子、HCO3-、F-、SO42-等阴离子及CO2、H2O等气相成分,随着成矿流体的逐渐演化分异,流体内的∑M+/∑M-逐渐升高。

(3) 成矿期黄铁矿微量元素研究表明,成矿热液的演化分异过程伴随着微量元素的演化现象,Y/Ho平均比值显示成矿物质的逐渐析出改变成矿流体性质,使得Zr/Hf、Nb/La、Co/Ni等稀土、高场强元素平均比值逐渐变小,而U/Th平均比值逐渐增高,同时还原性的成矿环境也会发生轻微波动,对成矿演化过程有较好的指示与约束作用。

致谢      本文的包裹体测试得到了路远发老师的指导和帮助;两位匿名审稿专家对本文提出了宝贵的修改意见;在此一并表示衷心的感谢!

参考文献
Bau M and Dulski P. 1999. Comparing yttrium and rare earths in hydrothermal fluids from the Mid-Atlantic ridge:Implications for Y and REE behaviour during near-vent mixing and for the Y/Ho ratio of Proterozoic seawater. Chemical Geology, 155(1-2): 77-90 DOI:10.1016/S0009-2541(98)00142-9
Bi XW, Hu RZ, Peng JT and Wu KX. 2004. REE and HFSE geochemical characteristics of pyrites in Yao'an gold deposit:Tracing ore forming fluid signatures. Bulletin of Mineralogy, Petrology and Geochemistry, 23(1): 1-4 (in Chinese with English abstract)
Bonnetti C, Liu XD, Mercadier J, Cuney M, Deloule E, Villeneuve J and Liu WQ. 2018. The genesis of granite-related hydrothermal uranium deposits in the Xiazhuang and Zhuguang ore fields, North Guangdong Province, SE China:Insights from mineralogical, trace elements and U-Pb isotopes signatures of the U mineralisation. Ore Geology Reviews, 92: 588-612 DOI:10.1016/j.oregeorev.2017.12.010
Chen BH, Wang ZL, Li HL, Li JK, Li JL and Wang GQ. 2014. Evolution of ore fluid of the Taishang gold deposit, Jiaodong:Constraints on REE and trace element component of auriferous pyrite. Acta Petrologica Sinica, 30(9): 2518-2532 (in Chinese with English abstract)
Chen YJ, Ni P, Fan HR, Pirajno F, Lai Y, Su WC and Zhang H. 2007. Diagnostic fluid inclusions of different types hydrothermal gold deposits. Acta Petrologica Sinica, 23(9): 2085-2108 (in Chinese with English abstract)
Chen ZY, Huang GL, Zhu B, Chen ZH, Huang F, Zhao Z and Tian ZJ. 2014. The characteristics and metallogenic specialization of granite-hosted uranium deposits in the Nanling region. Geotectonica et Metallogenia, 38(2): 264-275 (in Chinese with English abstract)
Chi GX and 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 (in Chinese with English abstract)
Chi GX and Lai JQ. 2009. Roles of fluid inclusions in study of mineral deposits. Mineral Deposits, 28(6): 850-855 (in Chinese with English abstract)
Cox SF, Knackstedt MA and Braun J. 2001. Principles of structural control on permeability and fluid flow in hydrothermal systems. In: Richards JP and Tosdal RM (eds.). Structural Controls on Ore Genesis, Structural Controls on Ore Genesis. Reviews in Economic Geology. Littleton: Society of Economic Geologists, 1-24
Deng J, Liu XF, Wang QF and Pan RG. 2015. Origin of the Jiaodong-type Xinli gold deposit, Jiaodong Peninsula, China:Constraints from fluid inclusion and C-D-O-S-Sr isotope compositions. Ore Geology Reviews, 65: 674-686 DOI:10.1016/j.oregeorev.2014.04.018
Deng MG, Chen W, Wang XW, Liu FX, Guan SJ, Lu YX, Yu HJ and Zhao F. 2018. Fluid inclusion and ore genesis of the Luziyuan distal skarn Pb-Zn-Fe (-Cu) poly-metallic deposit, West Yunnan, SW China. Acta Petrologica Sinica, 34(5): 1239-1257 (in Chinese with English abstract)
Deng P, Shu LS and Tan ZZ. 2003. The geological setting for the formation of rich uranium ores in Zhuguang-Guidong large-scale uranium metallogenetic area. Geological Review, 49(5): 486-494 (in Chinese with English abstract)
Du LT and Wang YM. 1984. Unification of metallogenic mechanism of granite-type, volcanic-type, carbosiliceous-mudstone-type and sandstone-type uranium deposits in South China. Radioactive Geology, (3): 1-10 (in Chinese)
Gao F, Lin JR, Zhong QL, Guo SY, Pang YQ, Rong JS and Hu ZH. 2011. The wall rock alteration and its geochemical characteristics of uranium deposit No.302. Uranium Geology, 27(5): 274-281 (in Chinese with English abstract)
Guo GL, Liu XD, Pan JY, Liu CD, Yan ZB and Chen YP. 2010. Study of fluid inclusion from uranium deposit No.302 in North Guangdong. Uranium Geology, 26(6): 350-354, 368 (in Chinese with English abstract)
Guo LN, Huang CM, Zhang L, Chen BH, Li RH and Liu Y. 2019. Source of ore-forming fluids in the Luoshan gold deposit, Jiaodong:Constrains from REE and trace element features of auriferous pyrite in the altered-rock type and auriferous quartz vein type ores. Geoscience, 33(1): 121-136 (in Chinese with English abstract)
Hagemann SG and Lüders V. 2003. P-T-X conditions of hydrothermal fluids and precipitation mechanism of stibnite-gold mineralization at the Wiluna lode-gold deposits, Western Australia:Conventional and infrared microthermometric constraints. Mineralium Deposita, 38(8): 936-952 DOI:10.1007/s00126-003-0351-6
Hall DL, Sterner SM and Bodnar RJ. 1988. Freezing point depression of NaCl-KCl-H2O solutions. Economic Geology, 83(1): 197-202 DOI:10.2113/gsecongeo.83.1.197
Henderson P. 1984. Rare Earth Element Geochemistry. Amsterdam: Elsevier, 1-510
Hollister LS and Crawford ML. 1981. Fluid Inclusions: Applications to Petrology. Mineralog: Mineralogical Association of Canada, 1-304
Hou ML, Jiang SY, Shen K, Lian GJ, Liu QC and Xiao FL. 2007. Fluid inclusion and H-O isotope study of gold mineralization in the Penglai gold field, eastern Shandong. Acta Petrologica Sinica, 23(9): 2241-2256 (in Chinese with English abstract)
Hu RZ, Bi XW, Su WC, Peng JT and Li CY. 2004. The relationship between uranium metallogenesis and crustal extension during the Cretaceous-Tertiary in South China. Earth Science Frontiers, 11(1): 153-160 (in Chinese with English abstract)
Hu RZ, Bi XW, Zhou MF, Peng JT, Su WC, Liu S and Qi HW. 2008. Uranium metallogenesis in South China and its relationship to crustal extension during the Cretaceous to Tertiary. Economic Geology, 103(3): 583-598 DOI:10.2113/gsecongeo.103.3.583
Hu RZ, Burnard PG, Bi XW, Zhou MF, Peng JT, Su WC and Zhao JH. 2009. Mantle-derived gaseous components in ore-forming fluids of the Xiangshan uranium deposit, Jiangxi province, China:Evidence from He, Ar and C isotopes. Chemical Geology, 266(1-2): 86-95 DOI:10.1016/j.chemgeo.2008.07.017
Huang GL, Yin ZP, Lin HF, Deng P, Zhu B and Shen WZ. 2010. Formation age, geochemical characteristics and genesis of pitchblende from No.302 uranium deposit in northern Guangdong. Mineral Deposits, 29(2): 352-360 (in Chinese with English abstract)
Huang GL, Cao HJ, Ling HF, Shen WZ, Wang XD and Fu SC. 2012. Zircon SHRIMP U-Pb age, geochemistry and genesis of the Youdong granite in northern Guangdong. Acta Geologica Sinica, 86(4): 577-586 (in Chinese with English abstract)
Huang GL, Liu XY, Sun LQ, Li ZS and Zhang SJ. 2014. Zircon U-Pb dating, geochemical characteristic and genesis of the Changjiang granite in northern Guangdong. Acta Geologica Sinica, 88(5): 836-849 (in Chinese with English abstract)
Li XH, Hu RZ and Rao B. 1997. Geochronology and geochemistry of Cretaceous mafic dikes from northern Guangdong, SE China. Geochimica, 26(2): 14-31 (in Chinese with English abstract)
Lin JR, Li ZY, Hu ZH, Wang YJ and Wang F. 2016. Controlling factors for the spatial positioning of hydrothermal uranium orefield. Uranium Geology, 32(6): 333-339 (in Chinese with English abstract)
Ling HF. 2011. Origin of hydrothermal fluids of granite-type uranium deposits:Constraints from redox conditions. Geological Review, 57(2): 193-206 (in Chinese with English abstract)
Liu B and Shen K. 1999. Thermodynamics of Fluid Inclusions. Beijing: Geological Publishing House, 1-290 (in Chinese)
Liu ZY and Liu HX. 2009. The simulation test of granite uranium mineralization. Earth Science Frontiers, 16(1): 99-113 (in Chinese with English abstract)
Lu HZ. 2011. Fluids immiscibility and fluid inclusions. Acta Petrologica Sinica, 27(5): 1253-1261 (in Chinese with English abstract)
Min MZ, Zhang ZH and Zhang BT. 1993. Mineralization features of the granite-type uranium deposits in southern China. Chinese Journal of Nuclear Science and Engineering, 13(1): 67-76 (in Chinese with English abstract)
Moh GH. 1980. Ore syntheses, phase equilibria studies and applications. Abhandlungen: Neues Jahrbuch für Mineralogie, 114-184
Monteiro LVS, Xavier RP, de Carvalho ER, Hitzman MW, Johnson CA, de Souza Filho CR and Torresi I. 2008. Spatial and temporal zoning of hydrothermal alteration and mineralization in the Sossego iron oxide-copper-gold deposit, Carajas Mineral Province, Brazil:Paragenesis and stable isotope constraints. Mineralium Deposita, 43(2): 129-159 DOI:10.1007/s00126-006-0121-3
Ni P, Tian JH, Zhu XT, Ling HF, Jiang SY and Gu LX. 2005. Fluid inclusion studies on footwall stringer system mineralization of Yongping massive copper deposit, Jiangxi Province, China. Acta Petrologica Sinica, 21(5): 1339-1346 (in Chinese with English abstract)
Ni P, Chi Z, Pan JY, Wang GG, Chen H and Ding JY. 2018. The characteristics of ore-forming fluids and mineralization mechanism in hydrothermal deposits:A case study of some typical deposits in China. Bulletin of Mineralogy, Petrology and Geochemistry, 37(3): 369-394 (in Chinese with English abstract)
Ni SJ, Hu RZ and Jin JF. 1994. A vertical zoning model generated by the mixing and boiling of hydrothermal solution for uranium deposit No.302. Uranium Geology, 10(2): 70-77 (in Chinese with English abstract)
Oreskes N and Einaudi MT. 1990. Origin of rare earth element-enriched hematite breccias at the Olympic Dam Cu-U-Au-Ag deposit, Roxby Downs, South Australia. Economic Geology, 85(1): 1-28
Qi JM, Huang GL, Zhu B, Fu SC, Xu Y and Ye SX. 2014. Compositions study of auxiliary minerals in altered granitic rocks of the Mianhuakeng uranium deposit in northern Guangdong. Acta Geologica Sinica, 88(9): 1691-1704 (in Chinese with English abstract)
Qi JM, Luo CW, Huang GL, Cao HJ and Xu LL. 2015. Geochemical character of pyrite and its tracing effects to ore forming fluid for granite type uranium deposit in north Guangdong. Uranium Geology, 31(2): 73-80 (in Chinese with English abstract)
Reich M, Kesler SE, Utsunomiya S, Palenik CS, Chryssoulis SL and Ewing RC. 2005. Solubility of gold in arsenian pyrite. Geochimica et Cosmochimica Acta, 69(11): 2781-2796 DOI:10.1016/j.gca.2005.01.011
Roedder E and Bodnar RJ. 1980. Geologic pressure determinations from fluid inclusion studies. Annual Review of Earth and Planetary Sciences, 8: 263-301 DOI:10.1146/annurev.ea.08.050180.001403
Shan Q, Zhang B, Luo Y, Zhou CP, Yu XY, Zeng QS, Yang WB and Niu HC. 2009. Characteristics and trace element geochemistry of pyrite from the Songhu iron deposit, Nilek County, Xinjiang, China. Acta Petrologica Sinica, 25(6): 1456-1464 (in Chinese with English abstract)
Shang PQ, Hu RZ, Bi XW, Liu L and Zhang GQ. 2006. Study on carbon and oxygen isotopes in granite-type hydrothermal uranium deposits:A case of the Xiazhuang uranium ore field in northern Guangdong. Journal of Mineralogy and Petrology, 26(3): 71-76 (in Chinese with English abstract)
Shao JL. 1988. Prospecting Mineralogy of Gold Deposit. Wuhan: China University of Geosciences Press, 38-45 (in Chinese)
Shen WZ, Ling HF, Deng P, Zhu B, Huang GL and Tan ZZ. 2010. Study on isotope geochemistry of uranium deposit 302 in northern Guangdong Province. Uranium Geology, 26(2): 80-87 (in Chinese with English abstract)
Shi SH, Hu RZ, Wen HJ, Sun RL, Wang JS and Chen H. 2010. Geochronology of the Shazijiang uranium ore deposit, Northern Guangxi, China:U-Pb ages of pitchblende and their geological significance. Acta Geologica Sinica, 84(8): 1175-1182 (in Chinese with English abstract)
Song YC, Hu WX, Jin ZJ and Chen Y. 2006. Fluid and melt inclusions and their fluid species in corundum megacrysts from the basalts in Changle, Shandong Province, eastern China. Geochimica, 35(4): 377-387 (in Chinese with English abstract)
Sun SS and McDonough WF. 1989. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes. In: Sanders AD and Norry MJ (eds.). Magmatism in Ocean Basins. Geological Society, London, Special Publication, 42(1): 313-345
Tu GC. 1998. Cryogenic Geochemistry. Beijing:Science Press: 93-106 (in Chinese)
Wang LL, Sun FY, Xu JH, Wang KY, Wang L and Xu QL. 2015. Ore-forming fluid and evolution of Luming molybdenum deposit, Heilongjiang Province. Acta Petrologica Sinica, 31(4): 979-990 (in Chinese with English abstract)
Wei H, Xu JH, Zeng QD, Wang YH, Liu JM and Chu SX. 2011. Fluid evolution of alteration and mineralization at the Duobaoshan porphyry Cu (Mo) deposit, Heilongjiang province. Acta Petrologica Sinica, 27(5): 1361-1374 (in Chinese with English abstract)
Wilkinson JJ. 2001. Fluid inclusions in hydrothermal ore deposits. Lithos, 55(1-4): 229-272 DOI:10.1016/S0024-4937(00)00047-5
Xue W, Xue CJ, Chi GX, Shi HG, Gao BY and Yang SF. 2010. Study on the fluid inclusions of Baiyangping poly-metallic deposit in Lanping basin, northwestern Yunnan, China. Acta Petrologica Sinica, 26(6): 1773-1784 (in Chinese with English abstract)
Yu DG. 1979. Evolution of Yanshan granite and uranium mineralization enrichment:Discussion on diagenesis and mineralization of some uranium-bearing rocks in South China. Uranium Geology, (4): 13-25 (in Chinese)
Zhang C, Cai YQ, Xu H and Liu JL. 2016. Mineralization mechanism of 302 uranium deposit, North Guangdong Province:Evidence from fluid inclusions. Journal of East China University of Technology (Natural Science), 39(2): 156-164 (in Chinese with English abstract)
Zhang C, Cai YQ, Xu H, Dong Q, Liu JL and Hao RX. 2017. Mechanism of mineralization in the Changjiang uranium ore field, South China:Evidence from fluid inclusions, hydrothermal alteration, and H-O isotopes. Ore Geology Reviews, 86: 225-253 DOI:10.1016/j.oregeorev.2017.01.013
Zhang CQ, Mao JW, Yu JJ and Li HM. 2007. Study on fluid inclusion and the metallogenetic mechanism of Chipu Pb-Zn deposit in Sichuan, China. Acta Petrologica Sinica, 23(10): 2541-2552 (in Chinese with English abstract)
Zhang GQ, Hu RZ, Shang PQ, Liu L and Yang SF. 2007. An overview on the ore-forming mechanism of the granite-type uranium deposit in South China. Bulletin of Mineralogy, Petrology and Geochemistry, 26(4): 399-404 (in Chinese with English abstract)
Zhang GQ, Hu RZ, Shang PQ, Tian JJ and Shuang Y. 2008. Study on the C-O isotopic composition of calcites and metallogenic dynamics background in the No.302 uranium deposi. Acta Mineralogica Sinica, 28(4): 413-420 (in Chinese with English abstract)
Zhang L, Chen ZY, Li SR and Huang GL. 2018. Characteristics of uranium minerals in wail-rock alteration zones of the Mianhuakeng (No.302) uranium deposit, northern Guangdong, South China. Acta Petrologica Sinica, 34(9): 2657-2670 (in Chinese with English abstract)
Zhao FF, Liu XF, Chu YT, Li CH, Song XF, Wu R and Xiao JX. 2011. Immiscible characteristics of mantle-derived fluid inclusions in special xenoliths from Cenozoic alkalic-rich porphyry in West Yunnan. Geochimica, 40(4): 305-323 (in Chinese with English abstract)
Zhong FJ, Pan JY, Xia F, Liu GQ, Qi JM, Guo XW and Li HD. 2017. Geochemical characteristics of rare earth elements in mineralization process in the Changjiang uranium ore field, northern Guangdong, China. Journal of Mineralogy and Petrology, 37(2): 63-73 (in Chinese with English abstract)
Zhong FJ, Pan JY, Qi JM, Yan J, Liu WQ and Li HD. 2018. New in-situ LA-ICP-MS U-Pb ages of uraninite from the Mianhuakeng uranium deposit, northern Guangdong Province, China:Constraint on the metallogenic mechanism. Acta Geologica Sinica, 92(2): 852-854 DOI:10.1111/1755-6724.13558
Zhou JY, Zheng RC, Zhu ZM, Chen JB, Shen B, Li XY and Luo LP. 2008. Geochemical characteristics of trace elements of pyrite and its implications to the metallogenesis in the Lala copper deposit. Journal of Mineralogy and Petrology, 28(3): 64-71 (in Chinese with English abstract)
毕献武, 胡瑞忠, 彭建堂, 吴开兴. 2004. 黄铁矿微量元素地球化学特征及其对成矿流体性质的指示. 矿物岩石地球化学通报, 23(1): 1-4. DOI:10.3969/j.issn.1007-2802.2004.01.001
陈炳翰, 王中亮, 李海林, 李金奎, 李京濂, 王国强. 2014. 胶东台上金矿床成矿流体演化:载金黄铁矿稀土元素和微量元素组成约束. 岩石学报, 30(9): 2518-2532.
陈衍景, 倪培, 范宏瑞, Pirajno F, 赖勇, 苏文超, 张辉. 2007. 不同类型热液金矿系统的流体包裹体特征. 岩石学报, 23(9): 2085-2108. DOI:10.3969/j.issn.1000-0569.2007.09.009
陈振宇, 黄国龙, 朱捌, 陈郑辉, 黄凡, 赵正, 田泽瑾. 2014. 南岭地区花岗岩型铀矿的特征及其成矿专属性. 大地构造与成矿学, 38(2): 264-275.
池国祥, 卢焕章. 2008. 流体包裹体组合对测温数据有效性的制约及数据表达方法. 岩石学报, 24(9): 1945-1953.
池国祥, 赖健清. 2009. 流体包裹体在矿床研究中的作用. 矿床地质, 28(6): 850-855. DOI:10.3969/j.issn.0258-7106.2009.06.013
邓明国, 陈伟, 王学武, 刘凤祥, 管申进, 卢映祥, 余海军, 赵发. 2018. 滇西芦子园远程矽卡岩Pb-Zn-Fe(Cu)多金属矿床流体包裹体初探及矿床成因探讨. 岩石学报, 34(5): 1239-1257.
邓平, 舒良树, 谭正中. 2003. 诸广-贵东大型铀矿聚集区富铀矿成矿地质条件. 地质论评, 49(5): 486-494. DOI:10.3321/j.issn:0371-5736.2003.05.005
杜乐天, 王玉明. 1984. 华南花岗岩型、火山岩型、碳硅泥岩型、砂岩型铀矿成矿机理的统一性. 放射性地质, (3): 1-10.
高飞, 林锦荣, 钟启龙, 郭树英, 庞雅庆, 戎家树, 胡志华. 2011. 302铀矿床围岩蚀变分带性及地球化学特征. 铀矿地质, 27(5): 274-281. DOI:10.3969/j.issn.1000-0658.2011.05.003
郭国林, 刘晓东, 潘家永, 刘成东, 严兆彬, 陈益平. 2010. 粤北302铀矿床流体包裹体研究. 铀矿地质, 26(6): 350-354, 368. DOI:10.3969/j.issn.1000-0658.2010.06.005
郭林楠, 黄春梅, 张良, 陈炳翰, 李瑞红, 刘跃. 2019. 胶东罗山金矿床成矿流体来源:蚀变岩型和石英脉型矿石载金黄铁矿稀土与微量元素特征约束. 现代地质, 33(1): 121-136.
侯明兰, 蒋少涌, 沈昆, 连国建, 刘其臣, 肖风利. 2007. 胶东蓬莱金矿区流体包裹体和氢氧同位素地球化学研究. 岩石学报, 23(9): 2241-2256. DOI:10.3969/j.issn.1000-0569.2007.09.022
胡瑞忠, 毕献武, 苏文超, 彭建堂, 李朝阳. 2004. 华南白垩-第三纪地壳拉张与铀成矿的关系. 地学前缘, 11(1): 153-160. DOI:10.3321/j.issn:1005-2321.2004.01.012
黄国龙, 尹征平, 凌洪飞, 邓平, 朱捌, 沈渭洲. 2010. 粤北地区302矿床沥青铀矿的形成时代、地球化学特征及其成因研究. 矿床地质, 29(2): 352-360. DOI:10.3969/j.issn.0258-7106.2010.02.017
黄国龙, 曹豪杰, 凌洪飞, 沈渭洲, 王小冬, 伏顺成. 2012. 粤北油洞岩体SHRIMP锆石U-Pb年龄、地球化学特征及其成因研究. 地质学报, 86(4): 577-586. DOI:10.3969/j.issn.0001-5717.2012.04.004
黄国龙, 刘鑫扬, 孙立强, 李钟枢, 张世佳. 2014. 粤北长江岩体的锆石U-Pb定年、地球化学特征及其成因研究. 地质学报, 88(5): 836-849.
李献华, 胡瑞忠, 饶冰. 1997. 粤北白垩纪基性岩脉的年代学和地球化学. 地球化学, 26(2): 14-31. DOI:10.3321/j.issn:0379-1726.1997.02.004
林锦荣, 李子颖, 胡志华, 王勇剑, 王峰. 2016. 热液型铀矿空间定位的控制因素. 铀矿地质, 32(6): 333-339. DOI:10.3969/j.issn.1000-0658.2016.06.002
凌洪飞. 2011. 论花岗岩型铀矿床热液来源——来自氧逸度条件的制约. 地质论评, 57(2): 193-206.
刘斌, 沈昆. 1999. 流体包裹体热力学. 北京: 地质出版社, 1-290.
刘正义, 刘红旭. 2009. 花岗岩铀成矿作用的模拟实验. 地学前缘, 16(1): 99-113. DOI:10.3321/j.issn:1005-2321.2009.01.012
卢焕章. 2011. 流体不混溶性和流体包裹体. 岩石学报, 27(5): 1253-1261.
闵茂中, 张祖还, 章邦桐. 1993. 华南花岗岩型铀矿床的成矿特征. 核科学与工程, 13(1): 67-76.
倪培, 田京辉, 朱筱婷, 凌洪飞, 蒋少涌, 顾连兴. 2005. 江西永平铜矿下盘网脉状矿化的流体包裹体研究. 岩石学报, 21(5): 1339-1346.
倪培, 迟哲, 潘君屹, 王国光, 陈辉, 丁俊英. 2018. 热液矿床的成矿流体与成矿机制——以中国若干典型矿床为例. 矿物岩石地球化学通报, 37(3): 369-394.
倪师军, 胡瑞忠, 金景福. 1994. 302铀矿床热液的混合和沸腾垂直分带模式. 铀矿地质, 10(2): 70-77.
祁家明, 黄国龙, 朱捌, 伏顺成, 许幼, 叶松鑫. 2014. 粤北棉花坑铀矿床蚀变花岗岩副矿物特征研究. 地质学报, 88(9): 1691-1704.
祁家明, 罗春梧, 黄国龙, 曹豪杰, 许丽丽. 2015. 粤北花岗岩型铀矿黄铁矿地球化学特征及对成矿流体的指示作用. 铀矿地质, 31(2): 73-80. DOI:10.3969/j.issn.1000-0658.2015.02.002
单强, 张兵, 罗勇, 周昌平, 于学元, 曾乔松, 杨武斌, 牛贺才. 2009. 新疆尼勒克县松湖铁矿床黄铁矿的特征和微量元素地球化学. 岩石学报, 25(6): 1456-1464.
商朋强, 胡瑞忠, 毕献武, 刘雷, 张国全. 2006. 花岗岩型热液铀矿床C、O同位素研究——以粤北下庄铀矿田为例. 矿物岩石, 26(3): 71-76. DOI:10.3969/j.issn.1001-6872.2006.03.012
邵洁涟. 1988. 金矿找矿矿物学. 武汉:中国地质大学出版社: 38-45.
沈渭洲, 凌洪飞, 邓平, 朱捌, 黄国龙, 谭正中. 2010. 粤北302铀矿床同位素地球化学研究. 铀矿地质, 26(2): 80-87. DOI:10.3969/j.issn.1000-0658.2010.02.003
石少华, 胡瑞忠, 温汉捷, 孙如良, 王加昇, 陈恒. 2010. 桂北沙子江铀矿床成矿年代学研究:沥青铀矿U-Pb同位素年龄及其地质意义. 地质学报, 84(8): 1175-1182.
宋玉财, 胡文瑄, 金之钧, 陈勇. 2006. 山东昌乐刚玉巨晶中的流体和熔融包裹体及其流体组分特征. 地球化学, 35(4): 377-387. DOI:10.3321/j.issn:0379-1726.2006.04.005
涂光炽. 1998. 低温地球化学. 北京: 科学出版社, 93-106.
王琳琳, 孙丰月, 徐九华, 王可勇, 王力, 许庆林. 2015. 黑龙江鹿鸣钼矿床成矿流体及演化. 岩石学报, 31(4): 979-990.
魏浩, 徐九华, 曾庆栋, 王燕海, 刘建明, 褚少雄. 2011. 黑龙江多宝山斑岩铜(钼)矿床蚀变-矿化阶段及其流体演化. 岩石学报, 27(5): 1361-1374.
薛伟, 薛春纪, 池国祥, 石海岗, 高炳宇, 杨寿发. 2010. 滇西北兰坪盆地白秧坪多金属矿床流体包裹体研究. 岩石学报, 26(6): 1773-1784.
余达淦. 1979. 燕山花岗岩演化与铀矿化富集——对华南一些含铀岩体成岩成矿关系探讨. 铀矿地质, (4): 13-25.
张闯, 蔡煜琦, 徐浩, 刘佳林. 2016. 粤北302铀矿床成矿机制探讨——来自流体包裹体的证据. 东华理工大学学报(自然科学版), 39(2): 156-164. DOI:10.3969/j.issn.1674-3504.2016.02.009
张长青, 毛景文, 余金杰, 李厚民. 2007. 四川甘洛赤普铅锌矿床流体包裹体特征及成矿机制初步探讨. 岩石学报, 23(10): 2541-2552. DOI:10.3969/j.issn.1000-0569.2007.10.022
张国全, 胡瑞忠, 商朋强, 刘雷, 杨社锋. 2007. 华南花岗岩型铀矿床成矿机理研究进展. 矿物岩石地球化学通报, 26(4): 399-404. DOI:10.3969/j.issn.1007-2802.2007.04.013
张国全, 胡瑞忠, 商朋强, 田建吉, 双燕. 2008. 302铀矿床方解石C-O同位素组成与成矿动力学背景研究. 矿物学报, 28(4): 413-420. DOI:10.3321/j.issn:1000-4734.2008.04.013
张龙, 陈振宇, 李胜荣, 黄国龙. 2018. 粤北棉花坑(302)铀矿床围岩蚀变分带的铀矿物研究. 岩石学报, 34(9): 2657-2670.
赵甫峰, 刘显凡, 楚亚婷, 李春辉, 宋祥峰, 吴冉, 肖继雄. 2011. 滇西富碱斑岩中特殊包体岩石的流体包裹体幔源不混溶特征. 地球化学, 40(4): 305-323.
钟福军, 潘家永, 夏菲, 刘国奇, 祁家明, 郭新文, 李海东. 2017. 粤北长江铀矿田成矿过程中稀土元素地球化学特征. 矿物岩石, 37(2): 63-73.
周家云, 郑荣才, 朱志敏, 陈家彪, 沈冰, 李潇雨, 罗丽萍. 2008. 拉拉铜矿黄铁矿微量元素地球化学特征及其成因意义. 矿物岩石, 28(3): 64-71. DOI:10.3969/j.issn.1001-6872.2008.03.011