西藏多龙斑岩型Cu-Au矿区深部至浅部构造框架与区域勘探靶区筛选：频谱约束的二维多尺度小波分析

《Minerals》：Deep-to-Shallow Structural Architecture and Regional Exploration Target Screening in the Duolong Porphyry Cu-Au District, Tibet: Spectrum-Constrained 2D Multiscale Wavelet Analysis

　　

位于班公湖-怒江缝合带（BNSZ）中部的多龙矿区是西藏主要的斑岩-浅成低温热液型Cu-Au矿集区。研究人员将径向平均对数功率谱与五级db4二维离散小波变换（2D-DWT）相结合，以比较深度敏感频带在完整布格重力和化极（RTP）航磁网格中的平面视图表达。五级变换产生从细尺度到粗尺度排序的有符号细节重建d1–d5，以及最终的累积低通近似a5。重力/磁法谱估计值：浅部集合为1.02/1.09 km，中间集合为5.12/4.85 km，最深分辨集合为18.60/21.58 km。这些值描述的是统计源深度范围，而非离散界面。深部源估计值之间的2.98 km差异（配对均值的14.8%）被视为方法学不确定性和/或密度与磁化强度转换深度的差异，而非定量一致性。最粗的细节d5和低通近似a5表达了宽泛的上部至中部地壳背景转换，并未成像莫霍面或地幔源。共定位的中尺度重力和磁异常与复合侵入-基底源区一致，但缺乏特定地点的密度和磁化率测量数据，无法进行唯一的岩性归属。因此，一个以约84.0°E、31.7°N为中心的东南异常区被提出作为宽3–8 km源深度带内的一级区域筛选假设，而非定量封闭的钻探目标。该解释区分了地球物理观测与地质假设，并明确纳入了谱和小波分离的局限性。

**论文解读：西藏多龙斑岩型Cu-Au矿区深部至浅部构造框架与区域勘探靶区筛选**

**研究背景、问题与意义**

斑公湖-怒江缝合带（Bangong-Nujiang Suture Zone，BNSZ）是青藏高原中部记录特提斯洋演化和羌塘-拉萨地体碰撞的关键构造带。位于该缝合带中段的多龙矿区是一个超大型斑岩-浅成低温热液Cu-Au（铜-金）矿集区，已探明铜资源超过1000万吨，具有重要的战略和经济价值。解析其深部至浅部地壳构造、以及侵入体、构造和热液系统之间的空间关系，对于评估造山背景下岩浆和流体运移模型至关重要。然而，多期叠加的构造-岩浆事件、广泛的地表覆盖以及强烈变化的围岩性质，使得常规地球物理解释面临挑战。完整布格重力异常和总场航磁数据包含来自不同深度源的叠加信号，从区域地壳对比到浅部侵入体和构造特征，使得深部地壳性质转换、中尺度异常源和浅部侵入-构造系统之间的空间关系尚未被充分约束。为此，研究人员开展了本研究，旨在通过频谱约束的二维多尺度小波分析（2D-DWT），建立从深部到浅部的统计源深度范围与平面尺度重建之间的关联，并筛选区域勘探靶区。该研究发表于《Minerals》，其意义在于提供了一种透明、可重复的尺度分离和假设构建工具，明确区分地球物理观测与地质解释，并系统纳入不确定性，为后续勘探工作提供一级筛选假设。

**主要技术方法**

本研究未使用样本队列，而是基于提供者提供的完整布格重力异常网格（mGal）和化极（RTP）航磁网格（nT）进行分析。关键技术方法包括：（1）径向平均对数功率谱（radially averaged log-power spectrum）用于统计源深度估计，通过分段线性拟合从不同波数窗口提取浅部、中间和深部源集合的深度范围；（2）五级db4二维离散小波变换（2D-DWT）用于非冗余多尺度滤波，重建从细尺度到粗尺度的有符号细节（d1–d5）及累积低通近似（a5），并将这些尺度与谱估计的深度窗口进行频谱引导的尺度耦合。此外，通过连续谱混合物和基于小波基的基准测试验证了方法的稳定性和分组能量的有序性。所有解释均基于谱和小波分离的局限性，并明确排除了柱状体或莫霍面等刚性界面的成像。

**研究结果**

**4.1 谱源深度估计（分段回归）**
通过分段最小二乘拟合，从径向平均对数功率谱中识别出三个近似斜率区间，分别对应三个统计源集合：
- 4.1.1 布格重力谱响应：深部源集合（0 < k ≤ 0.10 rad/km）深度为18.60 km；中间源集合（k = 0.10–0.42 rad/km）深度为5.17 km；浅部源集合（k = 0.42–1.65 rad/km）深度为1.02 km。
- 4.1.2 航磁RTP谱响应：深部源集合（0 < k ≤ 0.08 rad/km）深度为21.58 km；中间源集合（k = 0.08–0.35 rad/km）深度为4.85 km；浅部源集合（k = 0.35–1.48 rad/km）深度为1.09 km。
重力与磁法深度配对均值分别为1.055 km（浅部）、4.985 km（中间）和20.09 km（深部），深部差异14.8%被归因于方法学不确定性和/或密度与磁化强度转换深度的差异。

**4.2 小波分解与多尺度平面视图分析**
基于谱窗口顺序和相邻重建的持续性空间模式，将场分为三个操作尺度范围：浅部（d1–d2，对应约1 km谱集合）、中间（d3–d4，对应约5 km谱集合）和长波长（d5 + a5，对应18.60–21.58 km谱集合）。
- 4.2.2 布格重力多尺度平面视图特征：d1–d2显示局部梯度与短波长异常，与已知断裂吻合；d3–d4定义了多龙矿区及其周围的正负残差异常源区；d5和a5表达了区域地壳梯度，而非尖界面。
- 4.2.3 航磁RTP多尺度平面视图特征：d1–d2的尖锐短波长磁异常与已知斑岩中心局部重合；d3–d4的单个异常合并为更宽的高磁区，与磁化侵入系统或基底源一致；d5和a5的中低磁响应可能反映长英质/蚀变/剪切物质或居里温度效应，不具诊断性。

**结论总结与讨论**

讨论部分强调了多尺度异常的性质对应性和跨方法深度验证：浅部重力与磁法深度差异仅0.07 km，d1–d2异常与已知斑岩中心和构造重合，支持浅部侵入-构造关联；中间深度差异0.27 km，共定位的d3–d4残差暗示重叠的密度和磁化强度对比，但未证明均匀高密度酸性岩基或AFC（同化混染-分离结晶）驱动的磁铁矿富集；最深集合差异2.98 km，可能反映密度和磁化强度边界的解耦。研究还比较了向上延拓、连续小波变换和欧拉反褶积等替代方法，指出其局限性，并明确将未包括欧拉解作为主要解释约束，建议后续进行系统网格源位置分析、约束建模、物性采样和钻探。

研究结论部分翻译如下：
**6. 结论**
**6.1 方法学有效性**
结合谱与2D-DWT的工作流程将统计源深度窗口与有符号平面视图尺度重建相关联。谱截断为二进尺度的分组提供信息，但并非客观或唯一决定。重力和磁法估计在浅部和中间集合上高度一致；14.8%的深部差异被报告为不确定性及可能存在的物性解耦。该工作流程是一种尺度分离和假设构建工具，而非对唯一地质体的反演。
**6.2 地质发现与勘探意义总结**
结果区分了三个操作尺度范围并阐明了其证据限制：
- 深部范围（18.60–21.58 km）：从重力和磁法谱统计推断出宽泛的上部至中部地壳源深度范围，未成像莫霍面、地幔源、唯一断层倾角或直接地幔通道。
- 中间范围（4.85–5.12 km）：重叠的d3–d4异常定义了解释性异常源区，复合岩浆储层在物性上非唯一，且未得到物性证据支持的均匀高密度酸性岩基或AFC驱动的磁铁矿富集。
- 浅部范围（1.02–1.09 km）：d1–d2异常在空间上与已知斑岩中心和构造相关，但未唯一分辨单个矿体。
- 东南d3扰动定义了一个一级区域筛选靶区，位于约84.0°E、31.7°N，解释性3–8 km源深度带内，其顶部深度和体积未知，需通过系统数值网格源位置分析、约束建模、物性采样和钻探进行验证。