本文属于类型a，即单一原创研究的学术报告。以下是根据文档内容撰写的综合报告。

研究报告：集成真实构造与岩石物理信息的三维大地电磁约束反演——以铜陵矿区为例

一、作者、机构与发表信息

本研究由来自多个科研机构的学者合作完成。第一作者为RuiJin Kong（孔瑞金），通讯作者为Zhengyong Ren（任政勇），二人均来自中南大学（Central South University）地球科学与信息物理学院及关键矿产勘查与开发国家重点实验室。其他合作者包括来自重庆大学（Chongqing University）资源与安全工程学院的Huang Chen（陈煌）、中南大学数学与统计学院的Zhengguang Liu（刘正光），以及瑞士洛桑大学（University of Lausanne）地球科学研究所的Linan Xu（徐立南）。该研究成果发表于勘探地球物理学家学会（Society of Exploration Geophysicists）主办的期刊《Geophysics》第91卷第5期（2026年9–10月号），页码为B185–B210，文献编号为10.1190/geo-2024-0797，于2026年8月12日在线先行发表，2026年9月8日正式在线出版。

二、研究背景与目的

大地电磁测深法（magnetotellurics, MT）是一种通过观测地表或海底天然时变电磁场信号来探测地球内部电性结构的地球物理方法。由于其设备轻便、探测深度范围广，已被广泛应用于矿产资源勘探、地热资源探测、火山及俯冲带研究等领域。然而，在复杂地质环境中解释MT数据一直是地球电磁感应研究领域的难点。复杂地质环境主要包括起伏的地表地形和复杂的地下电性结构。传统的一维或二维反演方法往往忽略地形影响，导致解释结果不准确甚至错误。尽管近年来基于非结构化四面体网格的三维MT反演方法得到了快速发展，但现有方法大多缺乏同时整合真实地质构造信息和岩石物理约束的能力。这一局限主要源于在非结构化网格上构建包含真实地质与岩石物理信息的约束模型存在较大困难。

针对这一问题，本研究旨在开发一种能够同时整合详细构造信息和岩石物理信息的三维MT约束反演算法。具体目标包括：在正演计算引擎中实现基于非结构化网格的嵌套网格策略，以精确计算具有起伏地形的复杂地下结构的电磁响应；在反演目标函数中引入交叉梯度（cross-gradient）约束和引导模糊C均值（guided fuzzy c-means, GFCM）聚类约束，分别整合来自地质与地球物理数据的界面信息和来自岩石样品的电性参数信息；利用高斯–牛顿法（Gauss–Newton method）最小化目标函数以保证算法的快速收敛；并通过合成数据测试和铜陵矿区实际资料应用验证算法的有效性与可靠性。

三、研究流程与详细方法

本研究的方法流程可分为以下几个主要环节。

第一，正演建模与嵌套网格策略。为确保反演中的正演引擎能精确计算具有任意地形的复杂地下电性结构的电磁响应函数，研究团队在正演建模中实现了一种基于非结构化四面体网格的嵌套网格策略。首先构建一个相对稀疏的反演网格（inversion mesh），其单元数量通常限制在30万以内；然后以该稀疏网格为基础，在数据点附近采用二分加密方法生成更精细的正演网格（forward modeling mesh），其单元数量一般不超过50万。这种嵌套网格策略实现了反演网格与正演网格之间电导率的精确映射，消除了因网格间电性参数不匹配引起的数值误差。此外，使用较稀疏的反演网格也有助于缓解反演问题的不适定性。在正演计算中，采用了研究团队此前开发的基于棱边的有限元求解器（edge-based finite element solver）来计算MT响应函数。

第二，目标函数的构建。反演目标函数由四项组成：数据拟合项（data misfit term）、模型粗糙度项（model roughness term）、交叉梯度项（cross-gradient term）和GFCM聚类项（GFCM clustering term）。数据拟合项通过数据加权矩阵对正演响应与观测数据之间的残差进行加权。模型粗糙度项采用参考模型约束，参考模型为均匀半空间。交叉梯度项用于度量电阻率模型与先验速度模型之间的结构相似性，其定义为两个模型梯度的叉积。当两个模型的结构一致时，其梯度方向平行，叉积为零。为在四面体网格上生成空间梯度算子，采用了Lelièvre和Farquharson（2013）提出的在每个单元的小邻域内拟合线性趋势的方法。GFCM聚类项则用于将先验岩石物理信息引入反演。假设研究区内存在C类地质单元，其电阻率分别为t₁, t₂, …, t_C，GFCM聚类项通过隶属度值（membership value）uᵢⱼ和聚类中心vⱼ来约束模型参数向先验聚类中心靠拢。隶属度值在反演迭代过程中不断更新，而聚类中心则根据隶属度值和先验聚类中心的加权平均来计算。当某个先验聚类中心可靠时，赋予较大的权重因子ηⱼ；若不可靠或未知，则设ηⱼ为零。

第三，目标函数的最小化。采用高斯–牛顿法最小化目标函数。该方法需要求解一个大型线性方程组，其中涉及灵敏度矩阵（sensitivity matrix）的显式构建与存储、数据协方差矩阵、模型协方差矩阵以及由交叉梯度项和GFCM项导出的矩阵。模型更新向量的求解采用预处理共轭梯度法（preconditioned conjugate gradient, PCG），步长则通过线搜索方法确定。权衡参数（trade-off parameters）的计算采用了基于随机列向量的初始化方法，通过用户设定的λ和c参数来放大或抑制各约束项的影响。

第四，合成数据测试。研究设计了两组合成模型测试。第一组为带地形的双块体模型，包含一个低阻体（10 Ω·m）和一个高阻体（1000 Ω·m），背景电阻率为100 Ω·m，地表起伏在0至−2 km之间，共77个测点，频率范围为0.001–10 Hz的16个频点。第二组为带地形的层状模型，包含四层电阻率分别为10、100、1000和10000 Ω·m的介质层，地表起伏在0至−3.5 km之间，共171个测点，频率范围为0.0001–0.1 Hz的16个频点。合成数据中加入了5%的高斯随机噪声。为评估非理想先验信息对反演结果的影响，还设计了不完整和不正确的速度模型以及不完整和不正确的岩石物理信息等测试情景。

第五，铜陵矿区实际资料应用。铜陵矿区位于中国安徽省，是长江中下游成矿带的重要组成部分。该区域在中生代早期经历了强烈的构造运动和岩浆活动，成矿作用与岩浆活动密切相关。SinoProbe项目在该区域采集了6条MT剖面共51个测点的数据，频率范围为0.00055–320 Hz，共24个频点。由于研究区内缺乏完整的三维地震速度模型，研究团队利用已有的地震反射数据和区域地质信息构建了一个受地质约束的三维速度模型。该模型基于叠后时间偏移（poststack time-migrated, PSTM）地震剖面解释的关键构造特征（包括主要断裂和地壳界面），并采用平均地壳速度6000 m/s进行时深转换。岩石物理先验信息则通过对MT数据进行一维反演分析获得，确定了7个电阻率聚类中心（约1、3、15、75、175、750和3000 Ω·m）。

四、主要研究结果

在合成数据测试中，不同约束条件的反演方法表现出明显差异。第一组双块体模型的测试结果表明，无约束反演虽然恢复了两个异常体的基本特征，但高阻体的电阻率仅恢复到300 Ω·m（真实值为1000 Ω·m），且边界非常模糊。加入交叉梯度约束后，低阻体的边界清晰度显著提高，高阻体的电阻率恢复到400 Ω·m。加入GFCM约束后，两个异常体的边界都变得非常清晰，低阻体恢复到约10 Ω·m，高阻体恢复到约1000 Ω·m，与真实模型高度吻合。同时施加两种约束时，反演结果在结构和物性两方面均达到最佳。而未考虑地形的反演结果在地表产生了明显的虚假异常，尤其是在地形起伏较大的区域。模型均方根误差（model RMS）的统计进一步证实了这些观察：单独MT反演的模型RMS为113.6，加入速度模型后降至110.4，加入岩石物理信息后降至108.9，同时加入两种约束后进一步降至103.6，而未考虑地形的模型RMS最高，达125.4。

第二组层状模型的测试结果显示出类似的规律。无约束反演得到的各层电阻率值与真实值偏差较大（如第一层为5000 Ω·m而非10000 Ω·m，第四层为30 Ω·m而非10 Ω·m），且层界面模糊。交叉梯度约束显著增强了电性界面的清晰度，特别是第二和第三层之间的界面。GFCM约束使四层的电阻率均接近真实值，边界清晰。同时使用两种约束时结果最优。模型RMS从单独MT反演的2573.03降至同时使用两种约束时的2240.44，而未考虑地形的结果RMS高达2978.91。

关于非理想先验信息的测试表明，即使结构信息不完整（速度模型中只有低阻体对应的低速度异常而无高阻体对应的高速度异常），交叉梯度约束仍能改善低阻体的反演结果，使其边界更接近真实模型。不正确的结构信息则会显著扭曲反演结果。岩石物理信息不完整或存在一定误差时，反演结果的劣化程度相对较小，模型RMS仍低于无约束反演。

在铜陵矿区实际资料应用中，四种反演方法的数据拟合RMS分别为：无约束反演2.02，交叉梯度约束反演2.57，GFCM约束反演2.03，联合约束反演2.19。所有测点的观测数据均得到了较好的拟合。从水平切片和垂直剖面来看，交叉梯度约束使电性界面与已知断裂的空间展布更加吻合，GFCM约束使电性体更加聚焦，边界清晰。联合约束反演结果兼具二者的优势。研究识别出五个低阻异常体（C1–C5），并通过敏感性分析验证了这些异常体的可靠性。C1位于下地壳23 km深度，电阻率约10 Ω·m，被解释为含矿元素的下地壳部分熔融带。C2、C3和C4位于中地壳10–15 km深度，空间上与C1相连，被解释为热液流体上升的通道。C5与阳新–常州断裂（YCF）空间对应，电阻率低至1 Ω·m，深度约8 km。这些低阻通道与已知的铜陵中央断裂（TCF）和丁桥–戴汇断裂（DDF）相连，热液流体沿断裂运移并最终成矿，解释了地表矿床主要分布在断裂附近的观测事实。

五、结论与意义

本研究成功开发了一种能够整合真实构造信息和岩石物理约束的三维MT反演方法。该方法利用非结构化四面体网格精确模拟复杂地形，通过嵌套网格策略实现了正演与反演网格的独立构建和精确映射，通过交叉梯度方法整合地震速度模型提供的构造信息，通过GFCM聚类方法整合岩石样品提供的电性参数信息，并采用高斯–牛顿法实现反演的快速收敛。合成测试证明了该方法的有效性，铜陵矿区的实际应用则展示了其在复杂地质环境中的实用价值。

本研究的科学价值在于：提出了将地质构造约束与岩石物理约束同时纳入三维MT反演的完整技术方案，解决了现有基于非结构化网格的三维MT反演软件缺乏此类功能的局限。在实际应用方面，本研究为铜陵矿区提供了更为详细和可靠的三维地电模型，揭示了下地壳可能的部分熔融带和地壳内的岩浆通道，提出深部热液流体沿通道和断裂上升驱动成矿的解释，为理解长江中下游成矿带的成矿机制提供了新的地球物理证据。此外，本研究还讨论了今后将方法推广到音频大地电磁（AMT）数据和开展地震–MT联合反演的前景。

六、研究亮点

本研究的亮点主要体现在三个方面。第一，方法创新：将交叉梯度约束与GFCM聚类约束同时集成于基于非结构化四面体网格的三维MT反演框架中，并配套开发了嵌套网格策略，这是现有文献中尚未见报道的能力。第二，系统性的合成测试设计：不仅测试了理想先验信息条件下的反演效果，还系统地评估了不完整和不正确先验信息对反演结果的影响，为实际应用中先验信息的可靠性评估提供了参考。第三，实际应用价值：铜陵矿区的案例研究表明，整合多种先验信息后获得的三维地电模型不仅拟合了MT观测数据，还与已知地质构造和岩石物理信息一致，成功揭示了与成矿作用相关的深部结构和流体通道，体现了该方法在矿产资源勘查领域的应用前景。