断裂系统的双扩散对流与成矿元素富集研讨.pdfVIP

500 以科学发展蕊促进科技创新(上 断裂系统的双扩散对流与成矿元素富集研究 杨瑞琰h2马东升2 1．中国地质大学(武汉)数理学院，430074；2．南京大学地球科学系，210093 摘要由于双扩散对流可驱动成矿流体的运移，因此建立了在双扩散对流作用下成矿 流体富集机制的数学模型，并针对不同的断裂带规模和断裂带的倾向，以及在有外力作用流体 发生流动的情况时，对成矿元素的富集范围进行了数值模拟。模拟结果显示：断裂带的规模越 大在断裂带上部、对流层的下部区域成矿元素富集强度也越高；同时若流体的流动方向与断裂 的倾向一致将加速成矿元素的富集，若流体流动的方向与断裂的倾向相反将不利于元素的 富集。 关键词 地球化学 双扩散对流 盆地流体戍矿模式 盆地流体由于温度梯度而发生热扩散，浓度梯度而产生溶质扩散。当这两种扩散同时发生在盆地 流体系统时将导致盆地流体的双扩散对流(简称DDC)发生。将DDC引入到多孔介质流体领域是在20 世纪60年代末，主要研究热液系统和陆地地热系统。由于受实验条件和手段的限制，关于此项研究的 实验报告非常少，目前通行的做法是进行数值模拟研究。 本文作者在文献中，通过系统分析沉积盆地热液改造矿床的成矿特点，认为呈分散状态存在的成矿 元素从矿源层到赋矿层的漫长时期内，至少发生过一次富集过程。驱动成矿流体运移的驱动力，除了重 力驱动、压力驱动外。还可能存在双扩散对流驱动模式。其次利用双扩散对流方程建立了在双扩散对流 素的富集范围进行了数值模拟。数值模拟显示对不同的Ra取值，在沉积盆地基底隆起的上方或隆起 的斜坡部位，均发生了富集作用。模拟还证实若隆起部位存在垂直断裂，同时有盆地流体流向此隆起和 断裂区域，那么双扩散对流系统更有利于成矿元素的富集。 本文的研究是文献的继续，现在的断裂被研究的不再是垂直的，而是有不同倾向的。因此本文的研 究目的就是讨论不同倾向的断裂和流体的流动对成矿元素富集的影响。 一、双扩散对流的数学模型 多孔介质是指流体流动的介质环境(岩体和围岩)，尽管其孔隙度和渗透率有不同的变化，但流体的 流动基本是属于渗流，满足达西定律。在中低温矿床中，成矿作用一般总是与盆地流体的输运有关。根 据研究问题，我们考虑一个温度、溶质浓度梯度并存的盆地流体对流系统。在这相互藕合的盆地流体对 流系统内的盆地流体运动规律涉及到流体连续性方程、达西定律、能量守恒方程和溶质守恒方程。这些 方程经过变量无量缌化后，将有量缨的瀛函数口、时间c、温度T、浓度s以及坐标分量z和z化为无量 纲单位的流函数p、时间r、温度o、浓度s以及坐标分量x和z。从而得到均质多孔介质的无量纲 DDC控制方程组。(具体推导过程见文献的附录1)。 咀科学发展观促进科技创新(二) 501 丁。圳} 图1研究区域示意图 警+K磐一R。(裘咕孺OS) ㈩ 踟．a尘甜 a审瑚 O‘睁．O‘0 … 4磊+瓦葫一丽瑟2霜+万 L2， 、 1、 cs，… aXaz e箜Or+警筹一器荔=击(赛+雾)一R?azaX Le＼aXzaz2} 度的函数，如果在某时刻某处的浓度高于它的温度，我们就假设该时刻在该处溶质发生结晶沉淀，其沉 淀强度等于浓度与温度的差；如果该处的浓度小于其区域的温度，就认为不发生结晶沉淀(或溶解已结 晶的物质，本文没有考虑)。因此有： R=㈦oi ss圣80≤ ㈤ 如果多