Formational Mechanisms of Homogeneous Fluid and Boiling Fluid: Evidences from Synthetic Fluid Inclusions
-
摘要: 为详细了解流体的形成机制, 对系统的流体包裹体合成实验进行研究.研究表明, 在合成流体包裹体实验中, 广泛存在流体的均一化和沸腾作用; 流体的均匀与否, 与流体p-t轨迹在TP (H2O)-CP(H2O)-CP(NaCl-H2O) 曲线的部位有密切的关系.p-t轨迹在曲线上部的流体为均匀流体, 反之则为沸腾流体.但也有例外, 如在溶解曲线上被主矿物捕获的流体.这为本次研究一定条件下流体的形成机制、探讨成矿作用提供了理论依据.Abstract: Synthetic fluid inclusions were studied in order to discuss formational mechanisms of homogeneous fluid and boiling fluid in detail.The results of experiments of fluid inclusion synthesizing show frequent occurrences of homogeneous fluid and boiling fluid, was closely in relation to the position of its p-t locus in TP (H2O)-CP(H2O)-CP(NaCl-H2O) curve of p-t phase diagram of NaCl-H2O system determined whether the fluid was homogeneous.If it was at the upper zone of TP (H2O)-CP(H2O)-CP(NaCl-H2O) curve, then the fluid was homogeneous, otherwise, the fluid was boiling.But the fluid trapped on the solution curve of NaCl in water was an exception, for example No.10 experiment.This was very beneficial for us to study on the formation mechanism of fluid and mineralization in the ore deposit.
-
Key words:
- synthetic fluid inclusion /
- homogeneous fluid /
- boiling fluid /
- formational mechanism
-
流体是地壳和上地幔地球化学系统的重要组成部分, 流体的成分及状态会因其所处的环境的不同而有较大的差异.在过去的几十年中, 通过详细研究天然流体包裹体, 地质研究者已经提供了大量有关各种地质环境中流体包裹体的宝贵资料, 而且已经可以比较可靠地测出给定的地壳或上地幔环境中的包裹体类型及其岩相学和热力学特征(刘斌和沈昆, 1995; Philippot, 1996; Bakker and Diamond, 2000; Kendrich et al., 2001).但地学界对流体包裹体的捕获机理以及被捕获以后可能发生的变化仍然存在一些疑问.另外, 关于流体在适当的PTX条件下的相平衡和体积的数据仍然缺乏, 对包裹体数据解释的一些假设尚有待于证明.因此, 很多地质学家尝试用人工生长晶体时伴生的流体包裹体(Roedder and Kopp, 1975) 或采用专门技术合成的流体包裹体研究一些相关的问题(Zhang and Frantz, 1989; Jeffery et al., 1994; Schmidt et al., 1995).在本次研究中, 作者也试图通过流体包裹体的合成实验研究, 探讨流体的捕获机制.
1. 流体包裹体的合成实验
本次实验采用的方法类似于Sterner and Bodnar (1984), 合成实验在南京大学内生金属成矿机制国家重点实验室进行.具体方法是: 将人工合成水晶(水晶中无大于1 μm的包裹体) 加工成直径为4 mm, 长为10 mm的石英柱, c轴垂直于直径4 mm的圆面, 并用超声波+去离子水反复清洗石英柱, 直到干净为止.然后放入在马氟炉中于355~375 ℃条件下加热约45 min, 取出后立即在去离子水中淬火, 使石英柱产生大量的裂隙, 但不散裂, 并再次用超声波+去离子水反复清洗.然后将淬火清洗后的石英柱放在120 ℃的烘箱中干燥24 h除去裂隙中的水分.
将处理好的石英柱样品装入直径5 mm的黄金管中, 加入不同配比的物质和溶液.用炔氧焰焊封, 称重, 放入80 ℃的烘箱中烘24 h, 再称重以检验焊封的质量.焊封合格的黄金管装入内冷淬火高压釜中进行实验.实验条件见表 1.
表 1 合成流体包裹体实验条件Table Supplementary Table The experimental conditions of synthetic fluid inclusions实验时先把高压容器的压力加到预定值(传压介质为蒸馏水), 在此压力下4 h左右后将高压容器的温度升高到设定值, 在此条件维持反应.反应完成后, 在维持高压容器压力的条件下, 用压缩空气使其冷却淬火, 待其温度降到室温后卸掉压力打开高压容器, 取出黄金管, 称重, 并与装入时比较, 以判断反应过程中是否泄露.如无泄露则说明合成实验成功, 取出石英柱, 沿圆柱长度方向将其切成厚约3 mm的薄片, 然后磨成厚约1 mm两面抛光的测温片进行分析.
2. 包裹体温度与盐度测试
合成包裹体的盐度和均一温度测试在Linkam TMS-93型冷热台和Laborlux D型显微镜上进行, 目镜倍数10×, 物镜倍数50×.均一温度的测定步骤是: 为节省时间, 先以20 ℃/min的速率将包裹体加热至300 ℃, 再以1 ℃/min的速率继续加热, 至包裹体接近均一时改用0.2 ℃/min的速率进行, 直至包裹体均一.盐度测试视流体包裹体的类型而有所差异, 气液两相包裹体的盐度测定步骤是: 将包裹体以30 ℃/min的速率快速冷却至-100 ℃以下, 恒温一段时间, 然后用1 ℃/min的速率升温, 至包裹体中的第一块冰开始融化时改用0.1 ℃/min的速率, 直至包裹体中的最后一块冰融化, 获得冰点; 上述步骤重复2~3次, 结果取平均值; 在冷冻温度-盐度表中查出包裹体盐度(张文淮和陈紫英, 1993); 或包裹体过冷却后, 测定NaCl·2H2O晶体的消失温度, 利用公式NaCl·2H2O晶体熔化温度-盐度换算表获得包裹体盐度.压力由刘斌和沈昆(1999)所提供的公式p=a+b×t+c×t2求得.而含NaCl子矿物的多相包裹体的盐度测试步骤是: 在热台上测定NaCl子矿物的消失(或溶解) 温度, 然后由饱和溶液的含盐度公式和压力公式(Bischoff, 1991) 获得NaCl的含量和包裹体形成压力; 与多相包裹体共生的气液两相包裹体由于某些参数难以确定而未参加压力计算.其测试及计算结果见表 2.
表 2 合成流体包裹体的测试结果Table Supplementary Table The measurements of synthetic fluid inclusions从表 2可以看出, 样品1和2中均存在一期测试结果(盐度、压力、温度) 与实验条件相近的流体包裹体, 且不同包裹体的形成条件(盐度、压力、温度) 相似, 流体密度相同, 为均匀流体被捕获的结果.但第二期包裹体有所差异: 样品1中的这期包裹体为气液两相包裹体, 气液比为15%~70%, 盐度与压力与实验条件有很大的差异; 且这些包裹体沿裂隙呈定向排列, 密切共生, 不同气液比的包裹体的形成条件(温度、压力) 相似, 但流体密度有所差别; 另外, 不同气液比包裹体的均一方式也有所差异, 富气相包裹体均一到气相, 而富液相包裹体均一到液相, 因此, 这期包裹体为流体不均匀捕获的结果.而样品2的第二期包裹体与样品3-样品6中的包裹体比较相似, 即由两类流体包裹体组成: 一类为气液两相包裹体; 一类为含子矿物多相包裹体.两类包裹体密切共生、温度与压力相近, 但含子矿物多相包裹体的流体密度远大于气液两相包裹体.另外, 这两类包裹体的测试结果与实验条件均有较大的差别, 如温度、压力均低于实验的温度和压力, 气液包裹体的盐度小于实验盐度, 而含子矿物多相包裹体则大于实验盐度.因此这些包裹体也为流体不均匀捕获的结果.
样品7-样品9中的包裹体均为富液的气液两相包裹体, 尽管其捕获压力与实验压力有较大的出入, 但包裹体捕获温度、盐度与实验温度和盐度非常一致; 包裹体气液比也非常稳定, 维持在10%~20%之间; 而且, 各实验形成的包裹体中的流体密度几乎相同; 说明这些样品中的包裹体是从均匀流体中捕获的.
而样品10中, 只发现含子矿物多相包裹体, 且子矿物与气泡的消失温度一致, 包裹体均一成液相, 表明流体可能是在矿物溶解曲线上被捕获(刘斌和沈昆, 1999).另外, 包裹体捕获压力远远小于实验压力, 捕获温度也与实验温度有着一定程度的不同, 但包裹体盐度与实验盐度基本一致.而且, 包裹体流体密度为1.06~1.08 g/cm3, 与350 ℃ (假设值) 时被捕获的流体密度1.07 g/cm3非常一致.这说明从实验稳定到流体被主矿物捕获的过程中, 流体都处于比较均匀的状态.因此, 本样品中的包裹体也是从均匀流体中捕获的.
3. 流体形成机制探讨
在本次实验中, 流体可分为均匀流体和沸腾流体两类.
3.1 沸腾流体
沸腾流体是不均匀流体的一种特殊类型. 张德会(1997)研究表明, 该流体被广泛发现于浅成热液矿床、斑岩铜钼矿床及多金属中脉状矿床等矿床中, 是这些矿床矿物沉淀的重要机理.前人(Spycher and Reed, 1989; 杨巍然和张文淮, 1996; 肖新建等, 2004; 张振亮等, 2005; 郑有业等, 2006) 对不同矿床(断裂) 中的沸腾流体提出了多种形成机制.但这些机制是否正确, 还有待于实验证明.
上述研究表明, 流体不均匀捕获在本次实验反映的是流体因实验条件或外界环境影响而发生的性质变化.在本次实验中, 实验流体为事先配制的不饱和盐水溶液, 实验前不存在两种或多种流体共存的现象, 因此不均匀流体的形成只能由沸腾现象引起, 因不均匀捕获而形成的包裹体为沸腾包裹体.根据具体实验情况, 本文将不均匀流体的形成又分为两种情况: 一种为由单纯压力下降而引起的流体沸腾作用.该沸腾作用的起因是实验条件稳定以后高压釜发生缓慢泄漏(发现后被纠正), 但实验温度并没有随之改变, 如样品1和样品2中的第二期包裹体.由于流体压力在缓慢降低到临界线附近后被纠正, 没有继续降低的动力, 因此流体只能发生小规模的沸腾作用.样品2中的沸腾包裹体类型之所以与样品1不同, 在于样品2的沸腾包裹体形成温度要大大高于样品1, 形成条件较样品1远离临界曲线(图 1).另一种为温度、压力均有一定幅度下降而引起的沸腾作用, 而且下降过程一直在持续.这类作用主要发生在合成包裹体实验的退火过程中.实验关闭后, 高压釜在快速冷却(即将高压釜推出高温炉膛) 的过程中, 压力也持续下降.但由于冷却时间很短(一般为2~7 h), 金管中的流体很快达到沸腾条件, 形成低密度、低盐度的气相流体和高密度、高盐度的液相流体.当主矿物捕获到不同相的流体时, 可形成不同类型的包裹体.捕获低密度的气相流体时, 所形成的是气液两相包裹体(以纯气相包裹体和富气相的气液两相包裹体为主), 包裹体中流体的盐度、密度均较小, 如样品3-样品6中的气液两相包裹体; 捕获高密度、高盐度的液相流体时, 可能形成含子矿物的多相包裹体, 包裹体中流体的盐度、密度均较大, 如样品3-样品6中的含子矿物多相包裹体.含子矿物多相包裹体按照其子晶与气泡消失温度的先后, 可以分为两个亚类: 一为子晶先消失, 气泡后消失, 通过气泡的消失而使包裹体达到均一, 密度中等, 如样品3和5;另一个亚类为气泡先消失, 子晶后消失, 通过子晶的消失而使包裹体达到均一, 密度较高, 如样品4和6.
从图 1可以看出, 所有的沸腾作用都发生在不同盐度流体气相线与TP (H2O)-CP(H2O)-CP(NaCl-H2O) 曲线的下方L+V或V+石盐相图区.由于沸腾作用, 可导致主矿物分别捕获气相和液相或液相+H相的流体; 子矿物的来源, 除直接来自高温高压流体的过饱和结晶外, 还可来自液相流体冷却到常温常压后(因为其盐度可能高于常温常压条件下饱和溶液的盐度) 的饱和或过饱和结晶作用.必须注意, 不是所有的沸腾作用都可导致含子矿物多相包裹体的产生, 要视沸腾作用发生的条件(温度、压力) 和沸腾作用的时间而定.温度越高、压力越低、沸腾作用时间越长, 液相流体的盐度就越高, 形成含子矿物多相包裹体的可能性也越大.
还需注意的是, 不是所有的降温降压都可导致流体的沸腾作用.从图 1可以看出, 如果实验1的压力没有降低到250×105 Pa左右, 而是降低到400×105 Pa左右(A点) 后就不再下降, 流体的沸腾作用就不可能发生.另外, 如果实验1按照⑦的趋势线演化, 即使温度与压力都降低, 沸腾作用也不可能发生.
因此, 本次实验形成的沸腾流体主要跟流体的降压作用或降温降压作用有关, 但并不是所有的降压作用或降温降压作用都能使流体沸腾, 所有的沸腾作用都发生在一定盐度流体气相线与TP (H2O)-CP(H2O)-CP(NaCl-H2O) 曲线的下方L+V或L+V+H (石盐) 相图区的温度与压力范围内.
3.2 均匀流体
在本次实验中, 按照包裹体形成温度与压力条件, 均匀流体可划分为两类: 一类为同一实验中各包裹体在含盐度、压力和温度等方面相同或相近, 气液比也基本相同, 密度非常一致, 如实验1形成的第一期流体.该期包裹体的形成条件(捕获温度、捕获压力和盐度) 与实验条件基本一致, 其流体代表的是实验条件(温度、压力) 稳定以后, 没有受外界环境影响而发生变化的流体, 即从实验条件稳定后到退火之前的均匀流体.第二类为同一实验中形成的各包裹体的盐度与实验盐度基本一致, 流体密度相同, 但包裹体捕获温度或捕获压力与实验温度或压力有着较大的差距(有的还相差特别大, 如样品10), 如样品7-10、样品2的第一期流体.其中, 样品2的包裹体虽在流体盐度、形成压力方面与实验的盐度、压力基本一致, 但形成温度远低于实验温度, 其流体代表的是实验条件稳定之前加热阶段的均匀流体; 样品7-9中包裹体的捕获温度与实验温度基本一致, 但压力变化较大, 其流体为实验条件稳定到退火阶段开始时的均匀流体; 而样品10中包裹体的捕获温度、捕获压力与实验温度和压力有着较大的差距, 尤其是压力远小于实验压力, 其流体为实验退火阶段的均匀流体.
上述实验是在高压(500×105~1 000×105 Pa)、低温(350~500 ℃) 的条件下进行的, 实验稳定后流体均处于均一状态, 其p-t投影点(实验10除外) 位于TP (H2O)-CP(H2O)-CP(NaCl-H2O) 曲线上部(图 2).结合沸腾流体的p-t轨迹, 可以认为流体的均匀与否与流体p-t轨迹在TP (H2O)-CP(H2O)-CP(NaCl-H2O) 曲线的部位有着密切的关系.p-t轨迹在该曲线上部的流体为均匀流体, 如实验7-实验9、实验2的第一期包裹体, 尽管其流体的捕获温度、捕获压力与实验条件有所不同, 但其p-t轨迹仍然位于曲线的上部, 流体处于均一状态, 主矿物所捕获的包裹体流体盐度与实验盐度一致.这类流体, 在矿床(郑有业等, 2006)、油藏(陈振林和许浚远, 2007) 的矿物包裹体中时有发现, 是最为常见的流体.
但也有例外, 如实验10, 流体捕获的p-t轨迹位于TP (H2O)-CP(H2O)-CP(NaCl-H2O) 曲线下部的NaCl溶解曲线上, 所形成的包裹体群中只发现含子矿物多相包裹体, 无不同气液比例的气液两相包裹体, 各包裹体含盐度变化不大, 密度也非常一致, 说明流体捕获时为均匀饱和流体, 其原因还有待于进一步研究.
可以看出, 在不考虑成分因素的情况下, 流体的性质主要受温度和压力的影响.流体的p-t演化轨迹不同, 其性质可能有所不同, 如p-t轨迹在TP (H2O)-CP(H2O)-CP(NaCl-H2O) 曲线上部的流体为均匀流体, 反之则为沸腾流体.个别情况下, 流体也可在溶解曲线上被主矿物捕获(实验10).
在此曲线附近, 即流体处于临界点时, 流体的性质是突变的.如Ikushima et al. (1998)用拉曼光谱对临界点附近水结构的分析表明: 在临界点以下, OH伸缩峰的最大频率随着温度的增加而不断地增加; 到了临界点后再升高温度, 它的变化很缓慢, 出现了一个平台; 在超临界区, 当温度固定时, 最大峰的频率在临界压力附近出现极大值.陈晋阳等(2002)的研究也表明, CO2-H2O包裹体在均一(临界) 点附近, 流体拉曼光谱的半高宽和拉曼位移变化速率均一前明显快于均一后.Hoffmann and Conradi (1997)则证实了液态水的最大峰频率与单个气相水分子频率差值△f在临界点最小, 同样反映了流体性质的变化.这些现象说明, 接近沸腾的流体中水分子氢键正不断断裂、气化作用日趋增强, 均一性减弱, 不均一性逐渐增强, 其外因为流体在临界点附近的升温或降压作用.当到临界点后, 气化作用达到最大, 不均一性也最强, 氢键作用最弱.沸腾后, 尽管继续升温或降压, 但氢键作用和气化作用均达到极值后变化缓慢的结果表明流体不均一性的继续存在.
4. 结论
(1) 在本次实验中, 被主矿物所捕获的流体可分为均匀流体和沸腾流体. (2) 不同流体的形成, 主要跟流体p-t轨迹在TP (H2O)-CP(H2O)-CP(NaCl-H2O) 曲线的部位有着密切的关系.p-t轨迹在该曲线上部的流体为均匀流体, 反之则为沸腾流体.但也有例外, 如在溶解曲线上被主矿物捕获的流体(实验10). (3) 尽管上述研究还处于初步阶段, 但仍可以为我们研究一定条件下流体的形成机制、探讨成矿作用提供理论依据.
-
表 1 合成流体包裹体实验条件
Table 1. The experimental conditions of synthetic fluid inclusions
表 2 合成流体包裹体的测试结果
Table 2. The measurements of synthetic fluid inclusions
-
[1] Bakker, R. J., Diamond, L. W., 2000. Determination of thecomposition and molar volume of H2O-CO2 fluid inclu-sions by microthermometry. Geochimicaet Cosmochimica Acta, 64 (10): 1753-1764. doi: 10.1016/S0016-7037(99)00334-8 [2] Bischoff, J. L., 1991. Densities of liquids and vapors in boil-ing NaCl-H2O solutions: A PVTX summary from300℃to500℃. Amer. J. Sci., 291: 309-338. doi: 10.2475/ajs.291.4.309 [3] Chen, J. Y., Zheng, H. F., Zeng, Y. S., 2002. Raman spectro-scopic study on hydrogen bond of water molecules insynthetic inclusions under high temperature. Rock and Mineral Analysis, 21 (3): 166-170 (in Chinese with English abstract). [4] Chen, Z. L., Xu, J. Y., 2007. Dip of the oil (gas) -water inter-faces in anticline-hydrodynamic oil (gas) pools. Earth Science—Journal of China University of Geosciences, 32 (1): 89-92 (in Chinese with English abstract). [5] Hoffmann, M. M., Conradi, S., 1997. Are there hydrogenbonds in supercritical water. J. Am. Chem. Soc., 119 (16): 3811-3817. doi: 10.1021/ja964331g [6] Ikushi ma, Y., Hatakeda, K., Soito, N., 1998. An in-situ Ra-man spectroscopy studies of subcritical and supercriticalwater: The peculiarity of hydrogen bonding near thecritical point. J. Chem. Phys., 108 (14): 5855-5860. doi: 10.1063/1.475996 [7] Kendrich, M. A., Burgess, R., Pattrich, R. A. D., et al., 2001. Fluid inclusion noble gas and halogen evidence onthe origin of Cu-porphyry mineralizingfluids. Geochimicaet Cosmochimica Acta, 65 (16): 2651-2668. doi: 10.1016/S0016-7037(01)00618-4 [8] Liu, B., Shen, K., 1995. Formulae for calculating oxygen fugacities of fluid inclusions and their applications. Acta Mineralogica Sinica, 15 (3): 291-302 (in Chinese with English abstract). [9] Liu, B., Shen, K., 1999. Thermodynamics of fluidinclusion. Geological Publishing House, Beijing, 119-140 (in Chinese). [10] Philippot, P., 1996. The chemistry of high-pressure fluids (1 to3 GPa): Natural observations vs. experi mental constraints. Earth Sciences Frontiers, 3 (3): 39-48 (in Chinese with English abstract). [11] Roedder, E., Kopp, O. C., 1975. Acheck on the validity of thepressure correction in inclusion geothermometry, using hydrothermally grown quartz. Fortschr. Mineral., 52: 431. [12] Schmidt, C., Rosso, K. M., Bodnar, R. J., 1995. Syntheticfluid inclusions: XIII. Experimental determination of thePVT properties in the system H2O+40% NaCl+5mol% CO2 at elevated temperature and pressure. Geochimicaet Cosmochimica Acta, 59: 3953-3959. doi: 10.1016/0016-7037(95)00258-2 [13] Spycher, N. F., Reed, M. H., 1989. Evolution of a broad lands-type epithermal orefluid along alternative P-Tpaths: Impli-cations for the transport and deposition of base, precious, and volatile metals. Econ. Geol., 84: 328-359. doi: 10.2113/gsecongeo.84.2.328 [14] Sterner, S. M., Bodnar, R. J., 1984. Synthetic fluid inclusionsin natural quartz: I. Compositional types synthesized andapplications to experi mental geochemistry. Geochim. Cosmochim. Acta, 48: 2659-2668. [15] Xiao, X. J., Gu, L. X., Ni, P., et al., 2004. Esti mation of gas-escaping amount during ore-formation originated fromfluid boiling at massive sulfide deposits, Tongling re-gion. Uranium Geology, 20 (4): 91-98 (in Chinese with English abstract). [16] Yang, W. R., Zhang, W. H., 1996. Character of fault proper-ty and combination of fluid inclusions. Earth Science—Journal of China University of Geosciences, 21 (3): 285-290 (in Chinese with English abstract). [17] Zhang, D. H., 1997. Some new advances in ore-forming fluidgeochemistry on boiling and mixing of fluids during theprocesses of hydrothermal deposits. Advances in Earth Science, 12 (6): 546-552 (in Chinese with English abstract). [18] Zhang, W. H., Chen, Z. Y., 1993. Geology of fluid inclusion. China University of Geosciences Press, Wuhan, 108-112 (in Chinese). [19] Zhang, Y. G., Frantz, J. D., 1989. Experimental determination of the compositional limits of immiscibility in thesystem CaCl2-H2O-CO2at high temperatures and pres-sures using synthetic fluid inclusions. Chemical Geology, 74: 289-308. doi: 10.1016/0009-2541(89)90039-9 [20] Zhang, Z. L., Huang, Z. L., Rao, B., et al., 2005. Concentration mechanism of ore-forming fluid in Huize lead-zincdeposits, Yunnan Province. Earth Science—Journal of China University of Geosciences, 30 (4): 443-450 (in Chinese with English abstract). [21] Zheng, Y. Y., Gao, S. B., Zhang, D. Q., et al., 2006. Ore-forming fluid controlling mineralization in Qulongsuper-large porphyry copper deposit, Tibet. Earth Science—Journal of China University of Geosciences, 31 (3): 349-354 (in Chinese with English abstract). [22] 陈晋阳, 郑海飞, 曾贻善, 2002. 高温下合成包裹体中流体水分子氢键的拉曼光谱分析. 岩矿测试, 21 (3): 166-170. doi: 10.3969/j.issn.0254-5357.2002.03.002 [23] 陈振林, 许浚远, 2007. 背斜-水动力复合油(气) 藏油(气) 水界面产状. 地球科学——中国地质大学学报, 32 (1): 89-92. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200701012.htm [24] 刘斌, 沈昆, 1995. 流体包裹体的氧逸度计算公式及其应用. 矿物学报, 15 (3): 291-302. doi: 10.3321/j.issn:1000-4734.1995.03.009 [25] 刘斌, 沈昆, 1999. 流体包裹体热力学. 北京: 地质出版社, 119-140. https://cdmd.cnki.com.cn/Article/CDMD-11415-1012365217.htm [26] Philippot, P., 1996. 高压流体(1-3GPa) 的化学组成: 自然观察与实验对比. 地学前缘, 3 (3): 39-48. doi: 10.3321/j.issn:1005-2321.1996.03.004 [27] 肖新建, 顾连兴, 倪培, 等, 2004. 铜陵地区金属硫化物矿床沸腾流体成矿过程中气体逸失量的估算. 铀矿地质, 20 (2): 91-98. https://www.cnki.com.cn/Article/CJFDTOTAL-YKDZ200402004.htm [28] 杨巍然, 张文淮, 1996. 断裂性质与流体包裹体组合特征. 地球科学——中国地质大学学报, 21 (3): 285-290. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX603.009.htm [29] 张德会, 1997. 流体的沸腾和混合在热液成矿中的意义. 地球科学进展, 12 (6): 546-552. https://www.cnki.com.cn/Article/CJFDTOTAL-DXJZ706.006.htm [30] 张文淮, 陈紫英, 1993. 流体包裹体地质学. 武汉: 中国地质大学出版社, 108-112. [31] 张振亮, 黄智龙, 饶冰, 等, 2005. 会泽铅锌矿床成矿流体浓缩机制探讨. 地球科学——中国地质大学学报, 30 (4): 443-450. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200504008.htm [32] 郑有业, 高顺宝, 张大权, 等, 2006. 西藏驱龙超大型斑岩铜矿床成矿流体对成矿的控制. 地球科学——中国地质大学学报, 31 (3): 349-354. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200603009.htm -