DOi:10.16111/j.0258_7106.2017.01.005
西藏拉抗俄斑岩铜钼矿床流体包裹体研究 
唐攀1,陈毓川2,唐菊兴2**,郑文宝2,3,冷秋锋4,林彬2

(1 成都理工大学地球科学学院, 四川 成都610059; 2 中国地质科学院矿产资源研究 所 国土资源部成矿作用与资源评价重 点实验室, 北京100037; 3 中国地质大 学, 北京100083; 4 西藏华钰矿业股份有限公司, 西藏 拉萨850000)

第一作者简介唐攀, 男, 1989年生, 博士研究生, 矿床学专业。 Email: tangpan1 68@163.com
**通讯作者唐菊兴, 男, 1964年生, 博士生导师, 研究员, 主要从事矿床学和固体 矿产勘查与评价研究工作。 Email: tangjuxing@126.com

收稿日期2015_12_22

摘要:西藏冈底斯成矿带拉抗俄斑岩铜钼矿床位于西藏特提斯构造域拉萨地 块东段中南部 ,是近年来青藏高原地质大调查项目评价的重点矿床之一。文章在钻孔地质编录的基础上, 依据矿物组合、脉体穿切关系的不同,划分了3个成矿阶段:早阶段、中阶段以及晚阶段, 成矿中阶段为主成矿阶段。根据包裹体室温下的相态充填度特征以及是否含有子矿物,可将 其分为3大类:液相包裹体(Ⅰ)、气相包裹体(Ⅱ)和含子晶多相包裹体(Ⅲ)。成矿早 阶段流体包裹体主要为Ⅰ、Ⅱ和Ⅲ类包裹体,均一温度集中在260~400℃之间,w(N aCleq) 为2.1%~39.4%;成矿中阶段的流体包裹体主要为Ⅰ、Ⅱ和Ⅲ类包裹体,具有典型的沸腾包 裹体 组合,均一温度集中在280~360℃,w(NaCleq)为2.2%~37.1%;成矿晚 阶段流体包裹体主要为Ⅰ类包裹体,均一温度集中在220~280℃,w(NaCleq )为3.6%~5.6%,显示包裹体均一温度及盐 度呈递减趋势。成矿流体是中高温、中高盐度,且富含成矿元素的NaCl_H2O体系流体。矿 床形 成的压力为(100±10) MPa,形成深度为(3.7±0.4) km。激光拉曼探针分析结果表明 , 流体包裹体液相成分主要为H2O,气相成分含有CO2;子矿物有黄铜矿、磁铁矿、赤铁矿 、石膏、黄铁矿、金红石、长石等。成矿早阶段,岩浆房发生流体出溶;成矿中阶段,岩浆 流体 发生减压沸腾和不混溶作用。综上所述,温度降低、压力减小、pH值的增加以及氧化还原 电位的变化是造成拉抗俄矿床成矿元素沉淀的主要因素。
关键词: 地球化学;流体包裹体;岩浆_热液流体;拉抗俄斑岩铜钼矿床;冈底 斯;西藏
文章编号: 0258_7106 (2017) 01_0068_15 中图分类号: P618.41;P618.65 文献标志码:A 
A study of fluid inclusions from Lakange porphyry Cu_Mo deposit in Tibet 
TANG Pan1, CHEN YuChuan2, TANG JuXing2, ZHENG WenBao2,3, LENG QiuFen g4 and LIN Bin2

(1 College of Earth Sciences, Chengdu University of Technology, Chengdu 610059, Sichuan, China; 2 MLR Key Laboratory of Metallogeny and Mineral Assessment, In stitute of Mineral Resources, Chinese Academy of Geological Sciences, Beijing 10 0037, China; 3 China University of Geosciences, Beijing 100083, China; 4 Tibet Huayu Mining Development Co., Ltd., Lhasa 850000, Tibet, China)

Abstract:Located in central and southern part of east Lhasa block in the Tibetan Tethys t ectonic domain, the Lakange porphyry Cu_Mo deposit in the Gangdese metallogenic belt is one of the key deposits in the Tibetan Plateau for geological survey eva luation in recent years. On the basis of drilling core logging and according to the vein types and crosscutting relationships, the authors divided the mineraliz ation in the Lakange deposit into three stages: early ore_forming stage, ma in or e_formaing stage and late ore_forming stage. Fluid inclusions at Lakange can be classified into type Ⅰ(liquid phase inclusions), Ⅱ(vapour phase inclusions ) and Ⅲ(daughter mineral_bearing fluid inclusions). The early ore_forming st age has mainly type Ⅰ, Ⅱand Ⅲ inclusions, characterized by homogenization tem peratures 260~400℃, w(NaCleq) 2.1%~39.4%; the main ore_for ming stage has mainly type Ⅰ, Ⅱand Ⅲ inclusions, characterized by homogeniz ation t emperatures 280~360℃, w(NaCleq) 2.2%~37.1%; the late ore_fo rming stage has mainly type Ⅰinclus ions, characterized by homogenization temperatures 220~280℃, w(NaCl eq) 3.6%~5.6%,respectively, which suggests that the temperature and sal inity decrease from early to late stage. The ore_forming fluid is the NaCl_H2O system of rich ore_ forming elements, middle_high temperature and middle_high salinity. The formatio n pressure of deposit is (100±10) MPa, and the formation depth is (3.7±0.4) km. The laser Raman spectroscopic analysis reveals that H2O is the main compon en t in liquid, there exist minor amounts of CO2 in vapors, and daughter minerals are mainly chalcopyrite, magnetite, hematite, gypsum, pyrite, rutile and feldsp ar. The exsolution of magma chambers occurs at the early ore_forming stage; the initiation of boiling and immiscibility occurs at the main ore_forming stage . Th e major factors influencing the participation of ore_forming elements in the Lak ange deposit are decreasing temperature and pressure, increasing pH and variati on of oxidation_reduction potential.
Key words: geochemistry, fluid inclusion, magmatic_hydrothermal fluid, Lakange porphyry C u_Mo deposit, Gangdese, Tibet
        冈底斯成矿带是位于雅鲁藏布江缝合带与班公湖_怒江缝合带之间的世界级铜多金属成矿带 。近年来,冈底斯成矿带中南部找矿成果显著,已发现并评价了多个规模巨大的矿集区, 包括雄村铜金矿集区、厅宫_冲江铜矿集区、甲玛_驱龙_邦铺铜钼多金属矿集区、蒙亚啊_洞 中拉_亚贵拉铅锌矿集区等(唐菊兴等,2012)。目前,甲玛_驱龙_邦铺铜钼多金属矿集区 已探明铜资源量接近2000万吨,钼资源量超过150万吨(唐菊兴等,2012),是青藏高原研 究程度最高、探明铜多金属资源量最大的矿集区,也是中国重要的资源接续基地。拉抗俄斑 岩铜钼矿床位于甲玛_驱龙_邦铺铜钼多金属矿集区内,是其重要的组成部分(图1)。 
        拉抗俄斑岩铜钼矿床隶属于西藏达孜县,是20世纪80年代进行区域化探扫描时发现的化探异 常点。20世纪90年代到本世纪初,研究人员对该矿区进行了异常二级查证和矿点检 查工作。目 前,矿床的勘查和研究程度相对较低,仅在矿床地球化学、蚀变和成矿年龄方面开展过相关 研究(曲晓明等,2003;侯增谦等,2003;冷秋锋等,2015;2016)。成矿流体的功能是萃 取、溶解、搬运、沉淀和聚集成矿物质,沟通矿源场、运移场和储矿场的媒介与纽带(翟裕 生,1999)。成矿流体各种物理化学参数对研究岩浆流体性质、成矿过程、成矿物质来源和 成矿的物理化学条件等具有重要意义,也为判断剥蚀程度、指导找矿提供理论依据(杨志明 等,2005;李光明等,2007)。斑岩铜钼矿床的形成与浅侵位岩浆结晶作用的热液流体活动 有关,成矿流体大量出溶对斑岩矿床的形成具有重大意义(Ulrich et al., 2001;Richard s,2009)。因此,本文在详实的野外调研基础上,对矿床的流体包裹体进行了系统的岩相 学观察、激光拉曼成分测试和显微测温,确定了成矿流体的成分,初步探讨了岩浆_热液的 演化,为进一步丰富和完善整个矿床的理论研究奠定基础。
1矿床地质
        拉抗俄斑岩铜钼矿床位于冈底斯成矿带南部东段。矿区出露侏罗系叶巴组一、二段,侏罗系 多底沟组及第 四系。叶巴组一岩段(J1_2y1)分布于矿区中部,呈近东西向展布,产状总体 倾向北,倾角陡立,主要岩性为片理化安山岩、中基性火山角砾岩、绢云石英片岩,绿泥石 英片岩及凝灰岩。岩石片理化发育,变形变质特征明显。叶巴组二岩段(J1_2y 2)分布于矿区北部,与一岩段呈整合接触,呈近东西向展布,主要岩性为火山角砾岩、 晶屑凝灰岩夹变质细砂岩、砂板岩等。多底沟组(J3d)分布于矿区南部,与叶巴组 呈断层接触。主要岩性为细晶灰岩、生物碎屑灰岩、大理岩化灰岩夹砂岩。岩石大理岩化及 矽卡岩化显著。第四系(Q)主要沿水系河床及阶地分布,主要为全新统洪冲积砾石、砂砾石、含亚砂土层(图2)。
图 1西藏大地构造简图(a)和冈底斯地区构造单元划分及矿产分布图(b)(据程文斌等,20 10; 唐菊兴等,2012修编)
     1—中二叠统标志层; 2—岩浆岩; 3—地质界线; 4—缝合带; 5—断层; 6—省会; 7 —县
     BASZ—班公湖_怒江缝合带; IYSZ—雅鲁藏布江缝合带; GLZCF—噶尔_隆格尔_扎日南木错 _错麦断裂带; SLYNJOMZ—狮泉河_
    拉果错_永珠_纳木错_嘉黎蛇绿岩带; SMLMF—沙莫勒 _麦拉_洛巴堆_米拉山断裂; YZSZ—雅鲁藏布缝合带
     Fig. 1Tectonic sketch map of Tibet (a) and geological unit division and minera l distribution of Gangdese region (b) 
    (modified after Cheng et al., 2010; Tang et al., 2012)
     1—Middle Permian marker bed; 2—Magmatite; 3—Geological boundary; 4—Suture zone; 5—Fault; 6—Provincial capital; 7—County
     BASZ—Bangong Co_Salween River suture zone; IYSZ—Yarlung Zangbo River suture z one; GLZCF—Gaer_Longgeer_Zharinanmucuo_Cuo 
    mai fault zone; SLYNJOMZ—Siquan he _Laguocuo_Yongzhu_Namucuo_Jiali opholite zone; SMLMF—Shamolei_Maila_Luobadui_M ilasha fault; 
    YZSZ—Yarlung_Zangbo suture zone    
        矿区断裂构造发育,主要断裂构造呈东西向展布,次级断裂呈北西向或北东向展布(图2) 。 矿区侵入岩主要分布于矿区中东部和西部(图2),分布面积约占矿区面积的20%,为喜马拉 雅期中酸性侵入岩。岩石类型主要有花岗闪长斑岩、石英二长斑岩、二长花岗斑岩、花岗闪 长岩 及花岗斑岩。火山岩主要分布于矿区中部、北部地区,出露大面积的侏罗系叶巴组火山岩, 出露约占矿区面积的60%,主要岩石类型有侏罗系叶巴组的火山熔岩,火山角砾岩、安山岩 、流纹岩及凝灰岩等。与成矿有关的岩体为花岗闪长斑岩,灰白色_白色,中_细粒斑状结构 ,块状构造,斑晶含量约20%~25%,成分主要为斜长石、石英和黑云母,少量的钾长石和角 闪石,角闪石基本都蚀变为黑云母,斜长石大部分蚀变为绢英母。基质为长英质,主 要为斜长石、石英等。副矿物主要为锆石、榍石、磷灰石、磁铁矿等。
        矿区蚀变类型丰富,主要有钾化、硅化、绢云母化、高岭土化、角岩化;次为矽卡岩化、大 理岩化。钾化主要产于花岗闪长斑岩中心部位,主要为黑云母化,黑云母呈团斑状或脉状产 出,次生黑云母呈鳞片状交代角闪石,或呈集合体产出,或与石英共生组合成细脉;少量的 钾长石化,钾化石和石英共生呈脉状产出。绢英岩化主要位于花岗闪长斑岩体的内接触带与 钾化带两侧,表现为绢云母交代斜长石、钾长 石。高岭土化、硅化绢云母带:该带为发育范围最广的蚀变带,受多种蚀变叠加,最常见的是高岭土化、绢云母化,高岭石、绢云母主要 是长石等矿物蚀变产物,生成以高岭石为主的多种次生黏土矿物呈团块状分布,该带金属硫 化物少,矿化较差。铜钼矿化与石英绢云母化及钾硅化蚀变关系密切。
图 2拉抗俄矿区地质简图(据冷秋锋等,2015修编)
     1—第四系残坡积物、冲洪积物; 2—上侏罗统多底沟组一段; 3—中-下侏罗统叶巴组二 段; 4—中-下侏罗统叶巴组一段; 5—花岗斑岩; 
    6—花岗闪长斑岩; 7—地质界线; 8—实测断层; 9—矿体; 10—钻孔及编号; 11—勘探线及编号
     Fig. 2Geological map of the Lakange mining area (modified after Leng et al., 2015)
     1—Quaternary residual, slope, alluvial and diluvial materials; 2—1st member of Upper Jurassic Duodigou Formation; 3—2nd member 
    of Middle-Lower Jurassic Ye ba Formation; 4—1st member of Middle-Lower Jurassic Yeba Formation; 5—Granite porphyry; 6—Granodiorite 
    porphyry; 7—Geological boundary; 8—Measured fault ; 9—Ore body; 10—Drill hole and its serial number; 10—Exploration line 
    and it s serial number    
        脉体十分发育,类型丰富,贯穿成矿流体演化的各个阶段,是矿化和蚀变的之间记录。根据 矿物组合、脉体直接的穿插关系,将脉体划分为成矿早、中、晚阶段(表1)。
成矿早阶段无矿石英脉(图3a~b)、黑云母±石英脉(图3b),还可见少量的 磁铁矿脉(图3c)。矿物组合为石英、黑云母、钾长石及少量的磁铁矿。
成矿中阶段脉体产出连续,硫化物呈 线状分布,石英颗粒垂直于脉壁呈长柱状或对称 生长,主要包括5种类型:辉钼矿脉(图3d )、石英+辉钼矿±黄铜矿脉(图3e~g)、黄铜矿脉(图3h)、石英+黄铜矿±黄铁矿脉( 图 3j)、石英+黄铜矿+辉钼矿脉。矿物组合主要为石英、辉钼矿、黄铜矿、黄铁矿、绢云母以 及少量的绿泥石、绿帘石、黑云母。其中,石英+辉钼矿±黄铜矿脉、石英+黄铜矿±黄铁矿 脉分别为拉抗俄矿床钼、铜矿化最主要的产出形式。
        成矿晚阶段 多为近直立、高角度产出,主要包括3种类型:石英±黄铁矿脉、石膏±石英脉、方解石± 石英脉(图3k~m);矿物组合主要为石英、黄铁矿、方解石、石膏以及少量的方铅矿、辉 钼矿。
 表 1拉抗俄斑岩铜钼矿床主要脉体类型及其特征
     Table 1Paragenetic sequence of vein types and its characteristics in Lakange Cu_Mo deposit      
2研究方法
        本次研究样品主要采自拉抗俄矿床钻孔岩芯,岩性为与矿化有关的花岗闪长斑岩,采样时兼 顾不同矿化阶段、不同矿化及蚀变类型等方面的因素。将上述样品制成厚0.3 mm双面抛光 的 包裹体片,然后在偏光显微镜下观察,记录不同成矿阶段、不同类型包裹体的分布密度、大 小、气相分数等岩相学特征。
        流体包裹体显微测温工作在成都理工大学能源学院国家油气藏重点实验室完成,首先在光学 显微镜下辨识各种含包裹体矿物、包裹体类型并选出适合测温的流体包裹体,然后采用英国 Linkam THMSG 600型冷热台进行显微温度测试,温度控制范围-196~600℃,0~-196℃精度 为±0.5℃,0~600℃的精度为±2℃。测试过程中,升/降温速率为0.50~25.0℃,相变 点速率降为0.5~1.0℃。
        单个流体包裹体激光拉曼光谱显微探针分析在国土资源部成都地质调查中心实验室完成,所 用仪器为英国产的Renishaw RM2000型显微激光拉曼光谱仪,使用波长514 nm Ar+激光器 ,激光功率40 mW,曝光时间10 s/3次叠加,光谱仪狭缝10 μm。
3流体包裹体研究
3.1流体包裹体岩相学特征
        拉抗俄矿床各成矿阶段不同脉体石英中广泛发育流体包裹体,类型丰富(图4a~k)。根据 包 裹体成分及室温下的相态特征、充填度特征以及是否含有子矿物(卢焕章等,2004;Rusk e t al., 2008),可将其分为3大类:
        液相包裹体(Ⅰ)室温下由液相和气相组成,气相分数在10%~45%之间,加热后 均一为液 相,是矿区最发育的一类包裹体;多呈乳滴状、椭圆形、负晶形和不规则状等,粒径2~40 μm。与气相包裹体、含子晶三相包裹体密切共生。
        气相包裹体(Ⅱ)室温下由气相和液相组成,气相分数在45%~100%之间,可见 少量的 纯气相包裹体,加热后均一为气相;多呈椭圆形、乳滴状和负晶形等,粒径5~35 μm。
        含子晶多相包裹体(Ⅲ)室温下由气相、液相和子晶相组成,加热时部分为子晶 先消失, 部分为气泡先消失,最后都均一为液相,分布较广,各个阶段均发育;子晶以石盐为主,还 有钾盐子晶、石膏子晶和不透明子矿物等,呈负晶形、椭圆形和不规则状等,粒径4~30 μ m。石盐子晶多呈立方体,淡绿色,大小一般<10 μm;钾盐子晶呈浑圆状、长方体状;石膏 子晶呈板状。不透明矿物主要为黄铜矿、磁铁矿、赤铁矿和黄铁矿等,以黄铜矿、磁铁矿为 主;黄铜矿呈三角形、不规则状;磁铁矿呈不规则状,多数粒度较小。
   图 3拉抗俄斑岩铜钼矿床各种脉体及其穿切关系
     a. Q脉,具有Ksp蚀变晕,切穿斜长石斑晶; b. Q+Bt脉,切穿Q脉; c. Mt脉; d. Mol脉 ,被晚期的Q脉切穿; e. Q+Mol脉,切穿早期的两组Q脉; f. Q+Mol脉,切穿早期的Bt细脉 ; g. 关系明显的三期脉,产状较陡的Q2脉切穿Q+Mol脉和早期产状较缓无矿Q1脉,Q+Mol脉 切穿早期产状较缓Q1脉; h. Cp细脉,被早期的Mol细脉切穿; j. Q+Cp脉; k. 晚期Q脉; l. Gp+Q脉; m. Cal脉,切穿早期的
    Q_Mol脉
    Q—石英;Mt—磁铁矿;Bt—黑云母;Mol—辉钼矿;Cp—黄铜矿;Cal—方解石;Gp—石膏 
    Fig. 3Vein types and crosscutting relationships in Lakange Cu_Mo deposit
     a. Quartz vein with K_feldspar alteration, cutting a plagioclase; b. Quartz + bi otite vein, cutting quartz vein; c. Magnetite vein; d. Molybdenite vein, cut by quartz vein in later period; e. Quartz + molybdenite vein, cutting quartz vein; f. Quartz + molybdenite vein, cutting biotite vein; g. Quartz + molybdenite vein , cut by quartz vein of later period, but cutting quartz vein of early period; h . Chalcopyrite vein, cut by molybdenite vein of early period; j. Quartz + chalco pyrite vein; k. Quartz vein; l. Gypsum + quartz vein; m. Calcite vein, cutting q uartz + molybdenite vein 
    of early period
    Q—Quartz; Mt—Magnetite; Bt—Biotite; Mol—Molybdenite; Cp—Chalcopyrite; Cal— Calcite; Gp—Gypsum 
    图 4拉抗俄斑岩铜钼矿床中流体包裹体的显微照片
     a. 液相包裹体; b. 气相包裹体; c~i. 含子晶包裹体; j~k. 为不同种类的包裹体组 合
     Q—石英; Mol—辉钼矿; Cp—黄铜矿; Mag—磁铁矿; Hem—赤铁矿; Py—黄铁矿; Rt —金红石; Gp—石膏; Fl—长石; Hl—石盐; K—钾盐
     L—液相; V—气相; Ⅰ—液相包裹体; Ⅱ—气相包裹体; Ⅲ—含子晶多相包裹体 
     Fig. 4Micrographs of fluid inclusions of Lakange porphyry Cu_Mo deposit 
     a. Liquid phase inclusion; b. Vapor phase inclusion; c~i. Daughter mineral_bear ing fluid inclusion; j~k. Different types of inclusionsQ—Quartz; Mol—Molybdenite; Cp—Chalcopyrite; Mag—Magnetite; Hem—Hematite; P y—Pyrite; Rt—Rutile; Gp—Gypsum; Fl—Felspar; 
    Hl—Halite; K—Sylvite
     L—Liquid phase; V—Vapor phase; Ⅰ—Liquid phase inclusion; Ⅱ—Vapor phase inclusion; Ⅲ—Daughter mineral_bearing fluid inclusion     
        成矿早阶段发育Ⅰ、Ⅱ和Ⅲ类包裹体。中阶段包裹体发育Ⅰ、Ⅱ和Ⅲ类包裹体,Ⅲ类包裹体 类型发育较多,类型也较复杂,这些包裹体从含有一个石盐子矿物到含有几个透明子矿物( 石盐、石膏等)和不透明的子矿物(黄铜矿、磁铁矿、赤铁矿等)。晚阶段,发育 较多的Ⅰ类包裹体,少量的Ⅱ和Ⅲ类包裹体。
3.2流体包裹体激光拉曼分析
        拉曼探针是针对单个流体包裹体进行原位无损分析测定的最为有效的方法之一(张敏等,20 07)。
        对拉抗俄斑岩铜钼矿床中典型的Ⅰ、Ⅱ和Ⅲ类包裹体气、液、固三相做拉曼光谱分析 ,结果表明除寄主矿物石英的特征峰外,包裹体气相成分的拉曼谱图含宽泛的H 2O峰(3310~3610 cm-1),还有CO2特征峰(1386 cm-1和1283 cm -1);液相成分含宽泛的H2O峰(3310~3610 cm-1);子矿物有黄铜矿(290 cm -1和354 cm-1)、磁铁矿(534 cm-1和663 cm-1)、赤铁矿(60 9 cm-1和1314 cm-1)、石膏(1015 cm-1)、金红石(608 cm-1 )、长石(512 cm-1)和黄铁矿(340 cm-1和377 cm-1)等(图5a~k )。
3.3流体包裹体显微测温
        本文对拉抗俄矿床各成矿阶段295个流体包裹体进行了显微测温,测试的包裹体 均为原生包裹体,详细测温数据列于表2,相关图解见图6。
        成矿早阶段包裹体发育Ⅰ、Ⅱ和Ⅲ类包裹体。Ⅰ类包裹体均一温度在271~376℃ 之间,平均 值316℃;冰点温度为-1.2~-10.4℃,对应w(NaCleq)为2.1% ~14.4%。Ⅱ类包裹体均一温度在323~436℃之间,平均值389℃。Ⅲ类包裹体均一温度在2 51~ 315℃之间,平均值286℃;石盐子晶溶化温度为206~315℃(部分子晶先消失,部分为气泡 先消失),对应w(NaCleq)为32.2%~39.4%。
        成矿中阶段包裹体发育Ⅰ、Ⅱ和Ⅲ类包裹体。Ⅰ类包裹体均一温度在256~441℃ 之间,平均 值为325℃;冰点温度为-1.3~-5.4℃,对应w(NaCleq)为2.2%~8.4% 。Ⅱ类包裹体均一温度在300~440℃之间,平均值346℃;冰点温度为-5.5~-8.4℃ ,对应 w(NaCleq)为8.6%~12.2%。Ⅲ类包裹体均一温度在243~326℃,平均值 29 2℃;石盐子晶溶化温度为162~286℃(部分子晶先消失,部分为气泡先消失),对应w (NaCleq)为30.2%~37.1%。
        成矿晚阶段主要为Ⅰ类包裹体发育,少量的Ⅱ类包裹体。Ⅰ类包裹体均一温度在 185~342℃之 间,平均值250℃;冰点温度为-2.1~-3.3℃,对应w(NaCleq)为3 .6%~5.4%。Ⅱ类包裹体均一温度在210~326℃,平均值260℃;冰点温度为-2.2~ -3.4℃,对应w(NaCleq)为3.7%~5.6%。
4讨论
4.1成矿流体特征
        成矿早阶段发育Ⅰ、Ⅱ和Ⅲ类包裹体,均一温度出现2个峰值段(图6),分别集中于260~3 20℃和340~400℃,变化范围较大,在邦铺、朱诺和驱龙矿床也观察到同样的现象(罗茂澄 等, 2012;李淼等,2015;杨志明等,2009);w(NaCleq)为2.1%~1.64%和32 .2%~39.3%。造成该现 象的原因可能是受晚期流体叠加改造(李淼等,2015);Ⅲ类包裹体既可通过石盐溶化而 均一,也可通过气泡的消失而均一,许多这类包裹体的石盐子晶溶化温度小于气液相均一温 度,并且部分包裹体的石盐子晶溶化温度和气液相均一温度相差50℃以上,可能是由于捕获 不均匀流体形成(Rusk et al., 2008)。中阶段包裹体发育Ⅰ、Ⅱ和Ⅲ类包裹体,均一温 度集中于280~360℃之间,少量的包裹体均一温度在380~460℃(图6),可能是由于包裹 体 泄露导致均一温度偏高;w(NaCleq)为2.2%~12.2%和30.02%~37.1%。 该阶段发育较多的多相 含子矿物包裹体,这些包裹体从含有一个石盐子矿物到含有几个透明子矿物 (石盐、石膏等)和不透明的子矿物(黄铜矿、磁铁矿、赤铁矿等)。晚阶段,发育较多的Ⅰ类包裹体 , 少量的Ⅱ和Ⅲ类包裹体,均一温度集中于220~280℃,w(NaCleq)为3.6%~5 .6%。拉抗俄矿床 流体包裹体中的子矿物绝大部分是石盐子晶,少量钾盐子晶(图4);另外,激光拉曼测 试结果显示,气相成分主要为H2O,少量CO2;因此,本矿床的成矿流体主要为含NaCl 的水溶液,其他盐类组分很小,可以近似地作为NaCl_H2O体系(卢焕章等,2004)。
   图 5拉抗俄斑岩铜钼矿床流体包裹体激光拉曼图谱
     a. 液相和气相; b. 子晶; c. 气相; d~k. 子晶Cp—黄铜矿; Mag—磁铁矿; Hem—赤铁矿; Py—黄铁矿; Rt—金红石; Gp—石膏; Fl —长石 
     Fig. 5Transmitted light photomicrographs of Lakange fluid inclusions
     a. Liquid phase and vapor phase; b. Daughter mineral;c. Vapor phase; d~k. Daug hter crystal
     Cp—Chalcopyrite; Mag—Magnetite; Hem—Hematite; Py—Pyrite; Rt—Rutile; Gp—Gy psum;Fl—Felspar 
   表 2拉抗俄铜钼矿床不同成矿阶段流体包裹体显微测温结果
     Table 2Microthermometric properties of fluid inclusions from different ore_for ming stages at Lakange 
        综上所述,拉抗俄矿床的成矿流体为中高温、中高盐度,且富含成矿元素的NaCl_H2O体系 流体,从早期到晚期,温度和盐度逐渐减低。
4.2成矿压力及深
        流体发生沸腾作用,流体内部饱和蒸汽压和外部压力相等,故可用沸腾包裹体组合较精确的 估计斑岩矿床形成时的压力(卢焕章等,2004;Roedder et al., 1980)。拉抗俄矿床成矿 中阶段发育典型的沸腾包裹体组合(下文详细描述),因此,用此阶段的压力代表成矿的压 力 。本次压力估计,对w(NaCleq)为0~30%的低盐度流体的压力采用Knigh t等(1989)的压力公式计算,对w(NaCleq)>30%的中高盐度流体采用Beck er等(2008)的方法计算,获得成矿的压力为(100±10) MPa。
        若采用27 MPa/km的静岩压力,则成矿时的古深度约为(3.7±0.4) km。驱龙矿床含矿斑 岩 上部的叶巴组厚度在3~5 km之间(阴家润,2006;杨志明等,2008);流体包裹体估 算 获得的压力为(105±15) MPa~(90±20) MPa,形成深度在(4.2±0.6) km~(3.6 ±0.8) km(杨志明等,2009)。拉抗俄矿床位于驱龙矿床向西20 km处,含矿斑岩上部都 为叶巴组。 因此,本次估算的成矿压力和形成深度与地质观察基本一致。
4.3岩浆_热液过程
    4.3.1岩浆出溶流体
        岩浆出溶热液是斑岩成矿系统中的一个基本事实。Redmond等(2004)研究表明,在Bingham 斑岩矿床约3.5 km的古深度,岩浆分离出临界流体,并用之代表深部岩浆房中出溶形成的 流 体。Rusk等(2008)得出相似的结论,在Butte斑岩铜矿约9 km的古深度(目前矿体下部1.4 ~2.6 km处), 岩浆分离出一种单相低盐度岩浆流体(w(NaCleq)为2%~ 10%),形成石英_辉钼矿_黄铜矿脉,后期发生钾硅化 蚀变。Richards(2005)研究认为深部岩浆房出溶的流体为中等盐度的超临界流体,w(NaCleq)约为10%;张德会等(2001)研 究认为岩浆出溶压力≤100 MPa,出溶深度大致为3~4 km。
 图 6拉抗俄矿床不同阶段矿流体包裹体均一温度及盐度直方图
     Fig. 6Histograms showing homogenization temperature and salinity of different stage inclusions at Lakange   
        由激光拉曼成分测试可知,拉抗俄矿床的成矿早阶段,部分包裹体中的子矿物为长石、金红 石。该子矿物可能为机械捕获,也可能是流体直接沉淀形成。如果是流体直接沉淀形成,表 明流体继承了岩浆的某些成分,是岩浆出溶流体(岩浆_热液演化)的一个证据。另外,对 成矿早阶段少量的Ⅱ类包裹体加热可观察到临界均一,代表从岩浆房出溶后尚未发生相分离 的流体。
4.3.2不混溶作用
        Sillitoe(2010)研究认为,斑岩铜矿中的高盐度流体主要来自岩浆结晶晚期的二次沸腾, 即熔体中出溶的富含挥发分的低盐度岩浆流体上升至近地表时,向岩株或岩枝顶部聚集,内 压力不断增大,导致岩体顶部破裂,引起压力突然释放,造成岩浆流体发生减压沸腾作用和 不混溶作用,形成低密度的蒸汽和高盐度流体。
        拉抗俄矿床中,成矿中阶段Ⅰ、Ⅱ和Ⅲ类包裹体共存(图7),且具有相似的均一温度(表2 、 图6),说明拉抗俄矿床岩浆流体在上升过程中经历了沸腾作用和不混溶作用。从图8还可看 到,拉抗俄矿床包裹体盐度基本处于0~15%与30%~40%两个区间内,缺少中间的过渡盐度区 ,这也是由岩浆流体在温压条件下不稳定发生相分离导致(简伟等,2010)。
4.4成矿物质沉淀机制
        许多学者对斑岩矿床中金属特别是Cu的沉淀机制开展了详细的实验研究工作,研究表明,在 成矿过程中温度降低、压力降低、pH值增加或氧逸度增加,均可促进Cu的沉淀(Hezarkhani e t al., 1999; Ulrich et al., 2001; Redmond et al., 2004; Landtwing et al., 2005; 卢焕章,2011)。
        研究表明,不混溶作用可以破坏成矿流体的相平衡,相态变化导致溶液中金属络合物发生分 解并沉淀出金属矿物,是矿质沉淀的重要机制(Roedder, 1984; Reed et al., 2006; Klemm et al., 2008)。不混熔分离(或沸腾作用)即原始均 匀的流体,在地质演化过程中由于 多种原因,导致原始均匀流体发生不混溶分离,分成物理或化学性质不一致的两个相(卢焕 章等,2004)。拉抗俄矿床(图6)从成矿早阶段到晚阶段,温度和盐度逐渐降低, 说明温度降低导致成矿元素的沉淀,流体包裹体盐度降低。成矿中阶段,大量的Ⅰ、Ⅱ和Ⅲ 类包裹体共存,且具有相似的均一温度,说明由于压力降低发生了不混熔作用,导致流体的 相分离,从而导致CO2等挥发分的逸出,而这些挥发分的逸出能有效升高流体的pH值,降 低 矿物溶解度,促使黄铜矿、辉钼矿等矿石矿物沉淀(Drummond et al., 1985;张德会,199 7;卢焕章,2011)。
 图 7拉抗俄矿床中沸腾包裹体组合
     Fig. 7Boiling inclusion composition in Lakange图 8拉抗俄矿床流体包裹体均一温度度_盐度图解
     (底图据Calagari, 2004; Wilkinso n, 2001)
     Fig. 8Homogenization temperature versus salinity of 
    inclusions in Lakange ( base map after Calagari, 
    2004; Wilkinson, 2001)         
        流体包裹体中子矿物的组合常见2种类型:一种是金属子矿物多与石盐、钾盐等共存, 另外一种则是石盐、钾盐与石膏等共存;从子矿物氧化还原角度看,硬石膏—石盐—钾盐— 气相代表了氧化性质的流体,而石盐—钾盐—金属子矿物—气相代表了还原性质的流体。 而拉抗俄矿床成矿中阶段存在石盐、钾盐与氧化物石膏等共存(图4d、h、i),金属 子矿物与石盐、钾盐等共存(图4e~h)以及金属子矿物(黄铜矿、磁铁矿等)、氧化物( 赤 铁矿、石膏等)与石盐、钾盐共存的流体包裹体(图4d~f),说明流体的氧化还原电位发 生了变化,导致成矿元素沉淀。
因此,温度降低、压力减小、pH值的增加以及氧化还原电位的变化是影响拉抗俄矿床成矿元 素沉淀的主要因素。
5结论
        (1) 拉抗俄斑岩型钼铜矿床的形成划分为成矿早阶段、成矿中阶段以及 成矿晚阶段3个成矿阶段。流体包裹体分为Ⅰ(液相包裹体)、Ⅱ(气相包裹体)和 Ⅲ(含子晶多相包裹体)3类。成矿早阶段发育Ⅰ、Ⅱ和Ⅲ类包裹体;中阶段发育Ⅰ、Ⅱ和 Ⅲ类包裹体,Ⅲ类包裹体较多,类型也比较复杂;晚阶段主要发育较多的Ⅰ类包裹体,少 量的Ⅱ和Ⅲ类包裹体。
        (2) 激光拉曼探针分析结果表明,流体包裹体液相成分主要为H2O,气相成分含有CO2 ;子矿 物有黄铜矿、磁铁矿、赤铁矿、石膏、黄铁矿、金红石和长石等。成矿早阶段流体主要为Ⅰ 、Ⅱ和Ⅲ类包裹体,均一温度集中在260~400℃,w(NaCleq)为2.1%~39 .4%;成矿中阶段流体主要为Ⅰ、Ⅱ和Ⅲ类包裹体,具有典型的沸腾包裹体组合,均一温度 集中在280~360℃,w(NaCleq)为2.2%~37.1%;成矿晚阶段流体主要为 Ⅰ类包裹体,均一温度集中在220~280℃,w(NaCleq)为3.6%~5.6%; 从成矿早阶段到晚阶段 ,包裹体均一温度及盐度呈递减趋势。矿床形成的压力为(100±10) MPa,形成深度为(3.7 ±0.4) km。成矿流体是中高温、中高盐度,且富含成矿元素的NaCl_H2O体系流体。
        (3) 成矿早阶段,岩浆房发生出溶。成矿中阶段,岩浆流体发生减压沸腾和不混溶作用。 温度降低、压力减小、pH值的增加以及氧化还原电位的变化是影响拉抗俄矿床成矿元素沉淀 的主要因素。
    
        志谢感谢西藏地质六队张焕彬队长、西藏地调院刘鸿飞院长等为笔者野外工作 的大力支持;感谢成都地质矿产研究所包裹体实验室潘忠习老师的帮助;感谢审稿专家对本 文提出的宝贵修改意见。 
参考文献
 References   
     Becker S P, Fall A and Bodnar R J. Synthetic fluid inclusions. ⅩⅦ. PVT X proper ties of high salinity H2O_NaCl solutions (>30 wt% NaCl): Application to fluid in clusions that homogenize by halite disappearance from porphyry copper and other hydrothermal ore deposits[J]. Econ. Geol., 103(3): 539_554.
     Calagari A A. 2004. Fluid inclusion studies in quartz veinlets in the porphyry c opper deposit at Sungun, East-Azarbaidjan, Iran[J]. Journal of Asian Earth Sc iences, 23(2): 179-189.
     Cheng W B, Gu X X, Tang J X, Wang Li Q, Lv P R, Zhong K H, Liu X J and Gao Y M. 2010. Lead isotope characteristics of ore sulfides from typical deposits in the Gangdese_Nyainqentanglha metallogenic belt implications for the zonation of ore_ forming elements[J]. Acta Petrologica Sinica, 26(11): 3350_3362(in Chinese wit h English abstract).
     Drummond S E and Ohmoto H. 1985. Chemical evolution and mineral deposition in bo iling hydrothermal systems[J]. Econ. Geol., 80(1): 126_147.
     Hall D L, Sterner S M and Bodnar R J. 1988. Freezing point depression of NaCl_KC l_H2O solutions[J]. Econ. Geol., 83(1): 197_202.
     Hezarkhani A, Williams_Jones A E and Gammons C H. 1999. Factors controlling copp er solubility and chalcopyrite deposition in the Sungun porphyry copper deposit, Iran[J]. Mineralium Deposita, 34(8): 770_783.
     Hou Z Q, Qu X M, Wang S X, Gao Y F, Du A D and Huang W. 2003. Re_Os dating of mo lybdenite from the Gangdese metallogenic belt, Tibet: Applications of mineraliza tion time and dynamic background[J]. Science in China(Series D), 33(7):609_618 (in Chinese with English abstract).
     Jian W, Liu W and Shi L H. 2010. Advance in study of porphyry molybdenum deposit s[J]. Mineral Deposits, 29(2): 308_316(in Chinese with English abstract).
     Klemm L M, Pettke T and Heinrich C A. 2008. Fluid and source magma evolution of the Questa porphyry Mo deposit, New Mexico, USA [J]. Mineralium Deposita, 43(5 ): 533_552.
     Knight C L and Bodnar R J. 1989. Synthetic fluid inclusions: Ⅳ. Critical PVTX p roperties of NaCI_H2O solutions[J]. Geochimica et Cosmochimica Acta, 53: 3_8 .
     Landtwing M R, Pettke T, Halter W E, Heinrich C A, Redmond P B, Einaudi M T and Kunze K. 2005. Copper deposition during quartz dissolution by cooling magmatic_h ydrothermal fluids: The Bingham porphyry[J]. Earth and Planetary Science Lette rs, 235(1_2): 229_243.
     Leng Q F, Tang J X, Zheng W B, Zhang J S, Tang P, Yan G and Dong Y. 2015. Re_Os dating of molybdenite from the Lakange porphyry Cu_Mo deposit in Tibet and its g eological significance[J]. Geology in China, (2): 570_584(in Chinese with English abstract).
     Leng Q F, Tang J X, Zheng W B, Wang B H, Tang P and Wang H. 2016. Geochronology, geochemistry and zircon Hf isotopic compositions of the ore_bearing porphyry in the Lakange porphyry Cu_Mo deposit, Tibet[J]. Earth Science, 41(6): 999_1015 (in Chinese with English abstract).
     Li G M, Qin K Z, Zhang T P and Xiao B. 2007. High temperature, salinity and stro ng oxidation ore_forming fluid at Duobuza gold_rich porphyry copper deposit in t he Bangonghu tectonic belt, Tibet: Evidence from fluid inclusions[J]. Acta Petrologica Sinica, 23(5): 935_952(in Chinese with English abstract).
     Li M, Sun X, Zheng Y Y and Guo F. 2015. Characteristic of fluid inclusions of th e Zhunuo porphyry copper deposit in the Gangdese belt, Tibet[J]. Acta Petrolog ica Sinica, 31(5): 1355_1347(in Chinese with English abstract).
     Li Y Q, Rui Z Y and Cheng L X. 1981. Fluid inclusions and mineralization of the Yulong porphyry copper (molybdenum) deposit[J]. Acta Geologica Sinica, (3): 21 6_232(in Chinese with English abstract).
     Lu H Z, Fan H R, Ou X X, Shen K and Zhang W H. 2004. Fluid inclusions[M]. Beij ing: Science Press. 1_485(in Chinese).
     Lu H Z. 2011. Fluid in the earth[M]. Beijing: Higher Education Press. 329_330 (in Chinese). 
     Luo M C, Mao J W, Wang L Q, Leng Q F and Chen W. 2012. Fluid inclusion evidence for magmatic_hydrothermal evolution in the Bangpu porphyry molybdenum_copper dep osit, Tibet[J]. Acta Geoscientica Sinica, 33(4): 471_481(in Chinese with Engli sh abstract).
     Qu X M, Hou Z Q and Li Z Q. 2003. 40Ar/39Ar ages of the ore_bearing porphyries o f the Gangdese porphyry copper belt and their geological significances[J]. Act a Geologica Sinica, 77(2):245_252(in Chinese with English abstract).
     Redmond P B, Einaudi M T, Inan E E, Landtwing M R and Heinrch C A. 2004. Copper deposition by fluid cooling in intrusion_centered systems: New insights from the Bingham porphyry ore deposit, Utah[J]. Geology, 32(3): 217_220.
     Reed M H and Palandri J. 2006. Sulfide mineral precipitation from hydrothermal f luids[J]. Reviews in Mineralogy and Geochemistry, 61(1): 609_631.
     Richards J P. 2005. Cumulative factors in the generation of giant calc_alkaline porphyry Cu deposits[C].In:Porter T M,ed. Super Porphyry Copper & Gold Depos its PGC Publishing. 7_25.
     Richards J P. 2009. Postsubduction porphyry Cu_Au and epithermal Au deposits: Pr oducts of remelting of subduction_modified lithosphere[J]. Geology, 37(3): 247_250.
     Roedder E and Bodnar R J. 1980. Geologic pressure determinations from fluid incl usion studies[J]. Annual Review of Earth and Planetary Sciences, 8: 263_30l.
     Roedder E. 1984. Fluid inclusions[J]. Reviews in Mineralogy, 12: 644_668.
     Rusk B G, Reed M H and Dilles J H. 2008. Fluid inclusion evidence for magmatic_h ydrothermal fluid evolution in the porphyry copper_molybdenum deposit at Butte, Montana[J]. Econ. Geol., 103(2): 307_334.
     Sillitoe R H. 2010. Porphyry copper systems[J]. Econ. Geol., 105(1): 3_41.
     Tang J X, Duo J, Liu H F, Lang X H, Zhang J S, Zheng W B and Ying L J. 2012. Min erogenetic series of ore deposits in the east part of the Gangdise metallogenic belt[J]. Acta Geoscientica Sinica, 33(4): 393_410(in Chinese with English abst ract).
     Ulrich T, Gunther D and Heinrich C A. 2001. The evolution of a porphyry Cu_Au de posit, based on LA_ICP_MS analysis of fluid inclusions: Bajo de la Alumbrea, Arg entina[J]. Econ. Geol., 96: 1743_1774.
     Wilkinson J J. 2001. Fluid inclusions in hydrothermal ore de_posits[J]. Lithos , 55(1_4): 229_272.
     Yang Z M, Xie Y L, Li G M and Xu J H. 2005. Characteristics and forming process of ore_ forming fluids at Qulong copper deposit in Gangdise porphyry copper belt , Tibet[J]. Geology and Prospecting, 41(3): 21_26(in Chinese with English abstract).
     Yang Z M, Hou Z Q, Song Y C, Li Z Q and Xia D X. 2008. Qulong superlarge porphyr y Cu deposit in Tibet:Geology, alteration and mineralization[J]. Mineral Depo sits, 27(3): 279_318(in Chinese with English abstract).
     Yang Z M and Hou Z Q. 2009. Genesis of giant porphyry Cu deposit at Qulong, Tibe t: Constraints from fluid inclusions and H_O isotopes[J]. Acta Geologica S inica, 83(12): 1838_1859(in Chinese with English abstract).
     Yin J R, Cai H W, Zhou Z G, Zhang Y Y, Duan X and Xie Y W. 2006. Study of marine Triassic/Jurassic boundary stratigraphy and the latest Triassic mass extinction in Tibet[J]. Earth Science Frontier, 13(4): 244_254(in Chinese with Eng lish abstract).
     Zhai Y S. 1999. On the metallogenic system[J]. Earth Science Frontiers, 6(1): 13_27(in Chinese with English abstract).
     Zhang D H. 1997. Some new advances in ore forming fluid geochemistry on boiling and mixing of fluids during the processes of hydrothermal deposits[J]. Advance in Earth Sciences, 12(6): 546_552(in Chinese with English abstract). 
     Zhang D H, Zhang W H and Xu G J. 2001. Exsolution and evolution of magmatic hydr othermal fluids and their constraints on the porphyry ore_forming system[J]. E arth Science Frontier, 8(3): 194_202(in Chinese with English abstract).
     Zhang M,Zhang J F,Li L Q and Qiu L F. 2007. The application of laser Raman mic roprobe to the study of fluid inclusion[J]. World Nuclear Geoscience, 24(4): 2 4(4): 238_244(in Chinese with English abstract).
    
     附中文参考文献
    
     程文斌, 顾雪祥, 唐菊兴, 王立强, 吕鹏瑞, 钟康惠, 刘晓吉, 高一鸣. 2010. 西 藏冈底斯_念青唐古拉成矿带典型矿床硫化物Pb同位素特征_对成矿元素组合分带性的指示 [J]. 岩石学报, 26(11):3350_3362.
     侯增谦, 曲晓明, 王淑贤, 高永丰, 杜安道, 黄卫. 2003. 西藏高原冈底斯斑岩铜矿带辉钼 矿Re_Os年龄:成矿作用时限与动力学背景应用[J].中国科学(D辑), 33(7):609_618.
     简伟, 柳维, 石黎红. 2010. 斑岩型钼矿床研究进展[J]. 矿床地质, 29(2): 308_316.
     冷秋锋, 陈毓川, 唐菊兴, 郑文宝, 张金树, 唐攀, 严刚, 董宇. 西藏拉抗俄斑岩铜钼矿床 辉钼矿Re_Os同位素测年及其地质意义[J]. 中国地质, 2015,(2):570_584.
     冷秋锋, 唐菊兴, 郑文宝, 王保宏, 唐攀, 王豪. 2016. 西藏拉抗俄斑岩Cu_Mo矿床含矿斑 岩地球化学、锆石U-Pb年代学及Hf同位素组成[J]. 地球科学, 41(6): 999_1015.
     李光明, 李金祥, 秦克章, 张天平, 肖波. 2007. 西藏班公湖带多不杂超大型富金斑岩铜矿 的高温高盐度高氧化成矿流体:流体包裹体证据[J]. 岩石学报, 23(5): 935_952.
     李淼, 孙祥, 郑有业, 郭峰. 2015. 西藏冈底斯朱诺斑岩型铜矿床流体包裹体特征[J]. 岩石学报, 31(5): 1335_1347.
     李荫清, 芮宗瑶, 程莱仙. 1981. 玉龙斑岩铜(钼)矿床的流体包裹体及成矿作用研究[J ]. 地质学报, (3): 216_232.
     卢焕章, 范宏瑞, 倪培, 欧先习, 沈昆, 张文淮. 2004. 流体包裹体[M]. 北京: 科 学出版社. 1_485
     卢焕章. 2011. 地球中的流体[M]. 北京: 高等教育出版社. 329_330.
     罗茂澄, 毛景文, 王立强, 冷秋锋, 陈伟. 2012. 西藏邦铺斑岩钼铜矿床岩浆_热液流体演 化: 流体包裹体研究[J]. 地球学报, 33(4): 471_481.
     曲晓明, 侯增谦, 李振清. 2003. 冈底斯铜矿带含矿斑岩的40Ar/39 Ar年龄及地质意义[J].地质学报, 77(2):245_252.
     唐菊兴, 多吉, 刘鸿飞, 郎兴海, 张金树, 郑文宝, 应立娟. 2012. 冈底斯成矿带东段矿床 成矿系列及找矿突破的关键问题研究[J]. 地球学报, 33(4): 393_410.
     杨志明, 谢玉玲, 李光明, 徐九华. 2005. 西藏冈底斯斑岩铜矿带驱龙铜矿成矿流体特征及 其演化[J]. 地质与勘探, 41(2): 21_26.
     杨志明, 侯增谦, 宋玉财, 李振清, 夏代详, 潘凤雏. 2008. 西藏驱龙超大型斑岩铜矿床: 地质、蚀变与矿化[J]. 矿床地质, 27(3): 279_318.
     杨志明, 侯增谦. 2009. 西藏驱龙超大型斑岩铜矿的成因: 流体包裹体及H_O同位素证据[J ]. 地质学报, 83(12): 1838_1859.
     阴家润, 蔡华伟, 周志光, 张翼翼, 段翔, 谢尧武. 2006. 西藏海相三叠系—侏罗系界限及 晚三叠世生物灭绝事件研究[J]. 地学前缘, 13(4): 244-254.
     翟裕生. 1999. 论成矿系统[J]. 地学前缘, 6(1): 13_27.
     张德会. 1997. 流体的沸腾和混合在热液成矿中的意义[J]. 地球科学进展, 12(6): 546_ 552.
     张德会, 张文淮, 许国建. 2001. 岩浆热液出溶和演化对斑岩成矿系统金属成矿的制约[J ]. 地学前缘, 8(3): 194_202.
     张敏, 张建锋, 李林强, 邱林飞. 2007. 激光拉曼探针在流体包裹体研究中的应用[J]. 世界核地质科学, 24(4): 238_244. 矿床地质 2017 年  36 卷第 1 期36 (1): 68~82