岩浆铜镍硫化物矿床是指与镁铁-超镁铁岩相关的富含Cu、Ni、PGE的硫化物矿床，占世界上镍资源量的40%，铜资源量的5.5%，以及铂族元素资源量的96%，具有重要的经济意义和科学价值。针对该类矿床的研究和勘探工作已开展百年有余，直至上世纪80-90年代，人们陆续发现了加拿大Voisey’s Bay（图1a）、俄罗斯Noril’sk、以及中亚造山带内一系列赋存于小型侵入体中的大型-超大型铜镍硫化物矿床，岩浆通道成矿的模式才逐步得到广泛认可。但时至今日，人们对这一经典模型的解释却始终停留在定性演绎阶段，仍迟迟无法回答“上升岩浆能够携带多少硫化物液滴”以及“为什么会在岩浆通道的出口形成大规模块状矿化”等关键科学问题。
很多学者认为镁铁质-超镁铁质岩浆密度（～2700 kg/m3）本就大于地壳围岩，而硫化物的密度更甚（4000-5200 kg/m3）,推测上升岩浆很难携带硫化物液滴，而通道出口处的矿化则是硫化物渗滤回流所致。但是不可否认的是，夏威夷的镁铁质-超镁铁质火山熔岩中仍携带着大量的橄榄石斑晶，其密度也高于岩浆和围岩。因此，占决定性控制因素的应为深部中间岩浆房超压，而非显而易见的密度差异。当多期岩浆注入等因素使得岩浆房形成超压时，会使得围岩压裂，并驱使岩浆裹挟矿物斑晶及硫化物液滴向上运移。在此过程中，岩浆的超压不断被自身重力作用及粘滞阻力所消耗（图1b），若通道顶部的岩浆超压始终能够压裂上部围岩，则最终会喷出地表，否则侵位于浅部地层的薄弱带，从而形成含矿侵入体。
图1 (a)加拿大Voisey’s Bay 矿床剖面图； (b)岩浆通道不同区域内的压强平衡
中国科学院地质与地球物理研究所矿产资源研究院重点实验室姚卓森博士（现为加拿大卡尔顿大学博士后）在秦克章研究员和卡尔顿大学James Mungall教授的指导下，以垂直岩浆通道中的压力平衡为基础（图1b），细致厘定了单一硫化物液滴上升的动力学过程、岩浆对其最大荷载量的计算以及挥发份的影响和作用，初步建立了通道式岩浆成矿的定量化模型。
首先，通过理论计算证实，对于单个硫化物液滴而言，其在上侵镁铁质岩浆中粒度和成分基本保持不变，而免于受到破碎或元素扩散的影响，同时因密度差异而产生的相对沉降速度远小于岩浆上侵速度，故硫化物能够以相对稳定速度向上运移。随着所携带硫化物含量的增加，岩浆整体密度升高、粘度增加，使得上侵过程中岩浆超压的损耗增加，直至无法继续向上运移，而该情况所对应的硫化物比例即定义为岩浆的最大荷载量。对于镁铁质岩浆而言，其对硫化物的最大荷载量至少高于3%-5%左右。
另一方面，岩浆中若存在挥发份，会显著降低其密度，并适当降低整体粘度（图2），从而减小了重力和粘滞阻力在上侵过程中对岩浆超压的消耗，提升对硫化物的荷载能力。同时，上侵减压过程中，挥发份很可能会触发出溶作用，形成大量气泡，其与硫化物易于耦合形成混合液滴，密度小于岩浆，进一步促进了硫化物的向上运移。
图2 岩浆挥发份含量对(a)硫化物携带量(若只考虑密度差异)和(b)粘度的影响，其中灰色区域代表岩浆密度大于围岩密度
当岩浆裹挟着硫化物上侵至通道出口，则类似于流体力学中两相突扩流，已有的实验和数值模拟表明，岩浆在此区域内会形成再循环-回流区域(图3)，流速减弱，硫化物含量升高，非常有利于其沉淀和汇聚作用，从而形成最初的富集矿化。
相关工作首次建立了定量化模型，以期能够细致刻画开放通道体系下岩浆铜镍硫化物矿床的动力学运移过程，同时首次强调了岩浆挥发份对于成矿作用的重要意义，对我国金川、夏日哈木及中亚造山带内一系列铜镍矿床的勘探和研究均大有裨益。
图3 (a)两相突扩流及(b)两相缓扩流的示意图
研究成果发表于国际岩石学领域的权威期刊Journal of Petrology。（Yao Z*, Mungall J E, Qin K. A preliminary model for the migration of sulfide droplets in a magmatic conduit and the significance of volatiles[J]. Journal of Petrology, 2019, 60(12): 2281-2316. DOI: 10.1093/petrology/egaa005）（原文链接）成果受国家自然科学基金重点项目和国家重点科技研发项目资助。
