金成矿作用的流体-岩石相互作用模型化的研究.pdfVIP

矿床地质 1998拄 M1NERALDEPOSITS 第17卷增刊 金成矿作用的流体一岩石相互 作用模型化研究。 张哲儒 刘荣高 毛华海 (中国科学院地球化学研究所，矿床地球化学开放研究实验室，贵阳) 提要：考虑到水一岩相互作用过程中化学反应的热力学和动力学，流体流动的动力学和扩散 及弥散作用，建立了金成矿作用的流体．岩石相互作用动力学模型。以东北寨金矿床为例，设置初 始和边界条件，求得了数值解。计算了反应时间、总硫活度、热液温度、氧逸度及流体流速对金 从围岩中活化转入溶液的影响；同时还计算了当成矿溶液流经围岩发生化学反应金沉淀成矿时， 反应时间、流体流速的影响及黄铁矿和石英与金沉淀的相互关系。 关键词：金矿床 动力学水一岩相互作用 许多成矿作用与流体一岩石相互作用有关，尤其是低温成矿作用，在流体一岩石相互作用 体系中，化学反应和热质输运是关键问题。近年来国内外作了大量的工作，这方面具代表性 量迁移和化学反应的连续介质模型，提出了准稳态近似处理法和标度化。 础上，建立了微细浸染型金矿的溶液一岩石相互作用模型，以研究影响金活化和沉淀的物理 化学条件。改造分两个方面：一方面是把UNIX工作站版本的程序移植到微机上来；另一 方面是编制了输入文件编辑器，将原来程序的输入文件窗口化．并编制了输出文件的绘图接 口． 1理论基础 模型假设：①矿物粒间流体完全充满岩石iL隙并且是不可压缩的，其密度和粘度并不 随溶质含量而显著改变。②矿物粒间流体服从Darcy定律，岩石渗透率是岩石局部结构的 函数。③矿物生长和溶解的速率受表面反应控制。④溶质依靠平流、扩散和弥散作用进行 输运。⑤溶质在液相中的浓度和在矿物中的浓度相比为极小。⑥水相反应进行极快，配合 离子能迅速达到平衡。 *本研究为中科院“西南低温大面积成矿域研究”重大项目资助 张哲错，男，57岁，研究员．从事地球化学热力学与动力学研究。邮政编码：550002 1040 矿 床 地 质 1998＆ 在以上假设前提下，模拟程序充分考虑了水一岩相互作用过程中流体流动和化学反应之 下： 其中，C．代表溶液组分的摩尔浓度；Pf为流体密度；”为达西通量(U=中。)；中为孔隙 散系数D。与抑制因子F相除之和： D=Dh+等 R，为组分i的总反应速度，可分为多相反应的溶解一沉淀反应速率尺P(影响i组分在溶液 浓度的所有反应的速率总和)和单相溶液中的反应速率常数R产，即： R．=R?1”+R产； RP=一∑0imT。 其中7。为单位体积岩石内矿物m的沉淀与溶解反应速率，一般约定沉淀反应y。取正值， 中组分总数。 2模拟体系的初值和边界值条件 本计算以东北寨金矿床为例，根据东北寨金矿床的地质地球化学特征，模拟体系分两个 带，上带为泥灰岩，下带为板岩。根据原始容矿岩石的化学成分(高质彬等，1992)，计算 出模拟体系的初始矿物组成(体积百分数，表1)。 表1模拟体系矿物初始体积百分数 石英 伊利石 白云母 钾长石 白云石 方解石 黄铁矿 孔隙度 泥敏岩 28 12 9 5 20 15 1 10 板岩 35 41 3 2 25 25 25 儿5 结合矿物平衡热力学计算及流体包裹体特征，我们对微细浸染型金矿成矿动力学模拟体 系(溶液中)的初始条件和边界条件分别设定如下(表2)。 表2模拟体系初始条件与边界条件 Ⅱ 三Na ￡～ 醋 璐l ∑C = ECI =Mg ∑Ca