结构尺度对等间距热液成矿的控制及基于人工智能的预测方法

《Ore Geology Reviews》：Structural-Scale Controls on the Evenly-Spaced Hydrothermal Mineralization and AI-Based Predictive Methods

　　

受多尺度构造控制的等间距热液成矿是全球主要成矿带中常见的空间分布格局，已在成矿带、矿集区、矿床和矿体尺度上得到证实。然而，对其形成机制、主导控制因素、跨尺度关系及定量预测模型的系统综合研究仍然有限。本研究系统梳理了从成矿带到矿体/脉体尺度受构造控制的等间距成矿的代表性实例，总结了主要控制因素，并完善了构造控矿与空间预测的理论框架。研究结果揭示了主导控制机制随尺度变化的明显规律。随着观测尺度的减小，主要控制因素从区域应力场和构造背景转变为岩石力学性质和岩性非均质性。具体而言，在成矿带尺度上，间距主要受板块构造、岩浆作用和先存基底构造控制；在矿集区尺度上，受脆性盖层厚度控制；在矿床尺度上，受岩石力学性质和岩性非均质性控制；在矿体-脉体尺度上，受局部断层运动学、倾角和空间排列控制。以大梁子铅锌矿床为案例，通过PFC-3D模拟评估了断层间距对脆性盖层抗拉强度的敏感性。此外，研究人员通过整合野外边界条件与模拟获得的抗拉强度参数，开发了数据驱动的多层感知器回归模型。该人工智能模型成功预测了隐伏矿体的位置，与构造模拟结果和F2氟地球化学异常靶区高度吻合。未来工作应重点关注小尺度构造控矿机制、多源数据集成以及数据驱动的人工智能模型，以实现深部隐伏矿床的精准预测。

**论文解读：结构尺度对等间距热液成矿的控制及基于人工智能的预测方法**

**一、研究背景与科学问题**

受多尺度构造控制的等间距热液成矿是全球主要成矿带中普遍存在的空间分布格局，其间距范围从数米到数百公里不等，已在成矿带、矿集区、矿床和矿体等多个尺度上被广泛记录。这种空间有序性受多尺度构造相互作用和力学耦合机制的制约，为分析构造对成矿的控制作用以及深部隐伏矿床的定位预测提供了重要基础。

在力学机制研究方面，应力波干涉模型和应力阴影模型为均质脆性介质中周期性断裂排列提供了经典解释。Bai和Pollard进一步量化了断裂间距与层厚之间的关系，定义了标志着断裂饱和的临界间距-厚度比。然而，与纯构造间距相比，热液矿床的分布反映了构造变形、流体汇聚和矿物沉淀的综合效应。尽管研究已认识到等间距成矿是构造、岩性和力学因素耦合作用的结果，但仍存在三方面不足：第一，多数研究仅涉及单一尺度的机制，缺乏从成矿带到矿体尺度的统一认识；第二，不同尺度的主导控制因素各异，对各层次一级控制因素的认识存在不确定性；第三，传统力学模型和模拟实验多为定性或半定量，难以转化为隐伏矿体的预测工具。

近年来，机器学习和深度学习方法已被用于整合多源异构数据以识别复杂空间模式，但多数人工智能模型纯属数据驱动，缺乏明确的力学约束，很少纳入构造变形、断裂间距或岩石力学规律，限制了其在深部、勘探程度低地区的可解释性和可靠性。因此，迫切需要建立一个将多尺度构造机制与定量预测相结合的数据驱动人工智能框架。

**二、主要技术方法**

本研究采用多方法综合的研究路线。首先，通过系统收集和综合分析全球多个典型成矿带、矿集区、矿床和矿体尺度的等间距成矿实例，进行跨尺度构造控矿因素的综合归纳。其次，以四川西南部大梁子铅锌矿床为案例，采用基于三维离散元法（PFC-3D）的数值模拟方法，构建了2.7 km × 0.9 km × 0.45 km规模的砂箱式模型，力学参数依据矿区野外和室内试验标定，系统模拟不同抗拉强度条件下的Riedel剪切间距变化。再次，基于数值模拟产生的200个样本数据，构建了以岩石抗拉强度为单一输入特征、Riedel剪切带平均间距为目标变量的多层感知器回归模型，采用均方误差损失函数、ReLU激活函数、Dropout正则化、学习率调度器和早停机制进行训练。最后，将模型输出投影至矿区地质构造图，与三维离散元模拟结果、F2氟地球化学异常区和构造解译靶区进行对比验证。

**三、研究结果**

**3.1 跨尺度控制因素**

通过对不同尺度等间距热液矿床控制因素的综合研究，系统揭示了从成矿带到矿体尺度随尺度变化的耦合机制。在成矿带尺度，等间距分布主要受先存构造格架和板块俯冲所决定的沿弧岩浆产能的规律性波动控制，形成周期性分布的热源节点。在矿集区尺度，基底走滑断裂形成的平行等间距分布主要受脆性盖层厚度控制。在矿床尺度，Riedel剪切的间距不仅受脆性层厚度影响，还受岩石力学参数制约，其中抗拉强度变化对间距的影响最为显著。在矿体/矿段尺度，断裂和矿段间距对局部应力场变化表现出强烈依赖性，不同产状和构造性质的控矿断裂形成的雁列间距往往差异显著。

**3.2 成矿带尺度等间距分布特征与机制**

安第斯中部斑岩铜矿带中，Domeyko成矿带内相邻矿集区间距约100–134 km，总体空间周期约120 km；Pocuro成矿带内相邻矿集区间距约75–101 km，总体空间周期约85–90 km。两个成矿带均表现出双周期性特征：主要矿集区间半间距位置上发育次级成矿中心，间距约为主要矿集区的一半。这种等间距串珠状分布与先存基底断裂的规则间距密切相关，基底断裂在晚期挤压-张扭应力体制下重新活化，形成高角度相交的高渗透性流体通道，深部成矿流体沿交汇节点周期性上升并形成等间距成矿中心。

**3.3 矿集区尺度等间距分布特征与机制**

滇东北矿集区受小江断裂带及其派生的Riedel剪切系统控制，次级断裂平均间距约66 km，矿床呈带状或串珠状分布，相邻矿床间距一般为数十公里。砂箱实验表明，脆性层厚度是Riedel剪切间距最直接的控制因素，间距与厚度呈正相关，S/T归一化值约为1.32。内摩擦角对间距-厚度线性函数斜率有显著影响，而内聚力和基底摩擦对间距的影响可忽略。先存基底构造实验揭示了基底几何形态对Riedel剪切分布模式的显著控制：楔形基底对应弱到强递进型Riedel剪切，阶梯状基底对应突变型，地堑/地垒型基底对应间歇性和局部连续型。

**3.4 矿床尺度等间距分布特征与机制**

大梁子铅锌矿床属层控密西西比河谷型，矿体空间分布受右行张扭性F15断层及其派生的Riedel剪切断层控制，次级张扭-剪切Riedel破裂自西向东排列为F5、F6(N)、F6(S)、F8和F9，相邻矿体间距约150 m。F15断层在成矿期兼具导矿和容矿功能，形成负花状构造。三维离散元数值模拟表明：当侧向拉张占主导时，Riedel破裂趋于与主断层方向一致，应力分布变得更加密集，应力阴影区减小，断裂间距减小；抗拉强度增大时，断裂形成需要更多应力积累，Riedel剪切间距近似按幂函数增大。

**3.5 矿体尺度等间距分布特征与机制**

毛坪超大型铅锌矿床的矿体受NE向毛坪压扭性断层及其上盘毛毛山背斜构成的压扭-褶皱构造控制，矿体间距平均为数米至数十米。矿体在剖面上呈“宽-窄”分布，空间上呈倾斜排列。矿体侧伏方向与断层上盘运动方向一致，左行剪切导致左行侧伏，右行剪切导致右行侧伏。矿体的延伸和厚度受局部主压应力倾向和倾角变化共同控制，层间断裂带与其低序次成矿构造之间的夹角越小，矿体的延伸和厚度越大。

**四、人工智能预测模型**

**4.1 模型构建与训练**

基于三维离散元模拟获得的200个样本，构建了“输入层-隐藏层1-Dropout-隐藏层2-输出层”架构的多层感知器回归模型。模型以岩石抗拉强度为核心输入特征，以Riedel剪切带平均间距为目标变量。输入特征经标准化处理，采用均方误差作为损失函数，在隐藏层2中引入ReLU激活函数增强非线性表示能力，通过Dropout正则化减少对局部样本模式的依赖。训练过程中采用学习率调度器和早停机制控制训练进程，最大训练迭代次数设为5000。

**4.2 模型评估与应用**

模型预测曲线与离散元模拟样本的趋势高度一致，整体形态与回归拟合结果接近，表明模型能够捕捉抗拉强度变化对Riedel剪切间距的控制作用。损失函数在迭代过程中呈持续下降趋势并在后期趋于稳定，权重分布直方图显示参数主要集中在有限范围内，无极端发散或异常偏差。将模型输出投影至大梁子矿区地质构造图后，模型预测的深部隐伏矿体区域与三维离散元模拟结果及F2地球化学异常区高度吻合，基本满足深部矿体优选和精细靶区圈定的实际找矿要求。

**五、讨论与结论**

本研究通过综合野外构造分析、物理模拟和数值模拟，系统阐明了热液矿床从成矿带到矿体尺度等间距分布的尺度依赖性机制。主导控制因素表现出明显的尺度效应：成矿带尺度主要受先存基底构造和规律分布的岩浆-热液系统控制；矿集区尺度以脆性盖层厚度为主导；矿床尺度上岩石力学性质对Riedel剪切间距影响显著；矿体尺度则以局部断层的运动学特征、倾角和空间排列以及赋矿围岩的岩性非均质性为主要控制因素。

从成矿带到矿集区尺度，等间距分布主要受区域构造和流变学格架控制，受多因素耦合影响强烈；随着尺度缩小至矿床和矿体水平，区域构造应力场的影响明显减弱，而局部构造特征和岩性非均质性的控制作用显著增强。这一尺度转换模式为不同尺度的勘探策略制定提供了重要依据。

以大梁子铅锌矿床为案例，将三维离散元模拟获得的Riedel剪切间距数据与多层感知器回归模型相结合，构建了数据驱动的人工智能预测模型。将构造-力学机制纳入人工智能模型有效提高了预测结果的可解释性和可靠性。模型成功预测了隐伏矿体的空间位置，与构造模拟实验和F2地球化学异常圈定的靶区高度一致。

需要指出的是，本研究中的实验和数值模型均基于均质脆性层的假设，该假设适用于脆性盖层几乎全由岩性均匀的白云岩组成的大梁子铅锌矿床。然而，许多自然地质环境中脆性盖层常表现出明显的岩性非均质性、层理各向异性和先存构造不连续面，这些因素可能显著改变应力阴影模式，从而改变断裂间距。此外，当前模型仅输入了单一岩石力学参数（抗拉强度），无法完全捕捉自然成矿系统的复杂性。未来工作应纳入更多地质、构造和力学变量，以更好地再现地质背景和成矿过程，提高模型在更复杂地质环境中的预测能力和适用性。