基于TOUGHREACT的成矿过程化学反应数值模拟——以虎头崖铅锌多金属矿床为例

虎头崖铅锌多金属矿床是青海省祁漫塔格成矿带中典型的矽卡岩型矿床。本文基于前人对该矿区地质特征、矿物学和地球化学特征等方面的研究成果，利用TOUGHREACT数值模拟软件，针对石英-晚期硫化物成矿阶段进行成矿过程化学反应数值模拟研究。通过建立地质概念模型和网格化物化参数的设置，在设置的时间步长内对多相流体流动和热对流传导、溶质运移以及地球化学反应过程进行模拟计算，研究了矿区内不同温压条件下成矿元素化学平衡浓度变化的总体趋势以及对方铅矿和闪锌矿矿物溶解沉淀情况的影响。结果表明：当温度从250℃降至155℃时，Pb²⁺的化学平衡浓度出现骤减；当温度从300℃降至100℃时，Zn²⁺的化学平衡浓度出现骤减；由于矿物的沉淀需要消耗大量相关的成矿元素，方铅矿和闪锌矿最佳成矿温度范围分别为155℃~250℃和100℃~300℃；不同温度和压力对成矿元素化学平衡浓度影响不同，温度是控制方铅矿和闪锌矿沉淀化学反应的核心因素，压力的影响较小。模拟结果与传统地质解译结果相对符合，证明了模拟的有效性。