复杂构造区重磁电数据的三维联合反演及地质解释方法改进

摘要：针对断裂叠置与岩性横向突变并存的复杂构造区，提出重磁电三维联合反演的理论与方法体系。以共享地质界面为核心假设，耦合交叉梯度与岩性物性约束，构建方向性正则与深度衰减的多级约束；在算法上引入自适应步长与非单调线搜索，分块对角与代数多重网格混合预条件，并行化与混合精度提升规模适应性；在数据侧以小波域多尺度融合实现浅部细节与深部骨架协同。以德兴矿区东南段为例，数据拟合均方误差由一降至零点七，交叉验证残差方差由零点八四降至零点六三，界面与钻孔深度偏差中位数由五十五米降至三十八米，有效迭代次数由一百二十次降至七十二次，迭代时间由一百四十八分钟缩短至八十八分钟，钻孔命中率由百分之七十提升至百分之八十六。结果表明，该方法在保持地质合理性的同时显著提升了计算效率与解释可靠性。

关键词：重力磁力电法联合反演；交叉梯度约束；方向性正则；多尺度小波融合；复杂构造区地质解释

引言

造山带及其邻区常见多期断裂叠置与岩性横向突变，致使地球物理响应呈强非线性与显著等效性。重力与磁力对密度与磁化率变化敏感，电法对电性与含流体构造反应突出，单一物理场独立反演易将多物性耦合误判为单一异常，边界模糊与深部不确定性尤为突出。近年来的三维联合反演通过在体素层面耦合多物理量并共享地质界面，展示了整合优势，但在约束表达、尺度匹配与算力扩展方面仍存在不足，表现在地质先验难以量化注入，浅深尺度贡献失衡，以及百万级体素场景下迭代效率受限。面向矿集区勘查与深部找矿需求，亟需一种能在保持地质一致性的同时兼顾迭代效率与规模适应性的联合反演策略。本研究围绕约束条件优化、三维迭代效率提升与多尺度数据融合三个方面展开，力图在构造走向引导与岩性物性耦合下，对多物理量进行协同约束，在并行架构与自适应求解策略保障下实现大规模稳定收敛，并通过尺度域权重调度协调浅部细节与深部骨架，为复杂构造区提供更具可靠性的地下结构模型。

1 复杂构造区重磁电数据三维联合反演的理论框架

鉴于造山带断裂叠置与岩性横向突变并存的地质背景，重力与磁力观测对密度与磁化率变化敏感，电法观测对电导率与含流体结构响应突出。该类区域非线性与等效性并举，独立反演易把多物性耦合误判为单一效应，因此需把重磁电数据在三维空间内开展协同约束。联合反演以共享地质界面为核心假设，依托灵敏度互补与岩石物性关联，把多物理量在体素层面统一表达并纳入同一优化框架。

2 复杂构造区重磁电数据三维联合反演方法的改进

2.1 重磁电数据联合反演的约束条件优化

鉴于复杂构造区地质体边界受断裂叠置以及岩性横向突变的共同影响而表现模糊，本研究把地质先验信息引入约束条件体系来提升联合反演的地质合理性。围绕构造走向约束，我们把野外测绘的断裂走向、线性构造指向以及地震剖面解释出的主控边界等信息转化为方向性正则，用来引导密度、磁化率与电导率界面的梯度在给定走向内维持连续，并在走向变化显著的构造转折带设置分区权重，使约束对局部异方性保持适宜的灵活度。从物性耦合角度来看，密度‑磁性‑电性关联约束以岩石物性数据库与井孔测井资料为依据，构建岩性类别的物性联合分布，把各物理量的变化趋势按岩性域进行耦合；例如在基性与超基性岩域强调高密度与高磁化率的协同变化，同时对含流体破碎带强调高电导率与低磁化率的伴生特征。需重点关注的是，先验信息的不确定性会对反演稳定性产生影响，因此设置自适应权重与置信区间，把约束强度随数据一致性进行动态调节，避免过度收敛到先验而削弱真实异常的表达。由此，约束条件优化形成分级结构：区域尺度的走向引导提供宏观几何框架，断裂邻域的局部窗口约束刻画断裂核部与损伤带的界面细节，岩性域的物性关联把多物理量在体素层面进行协同收敛，使几何边界与物性变化在统一模型内彼此契合。与此同时，引入深度分层的约束衰减策略，凭借地表与浅层资料可信度较高而给予较强引导，向深部逐步减弱，使联合反演对深部结构的探索空间保留适度自由度并提高解释的稳健性。

2.2 三维反演算法的迭代效率提升

鉴于三维联合反演在百万至千万体素规模下呈现高维非线性与强耦合特性，本研究把迭代策略与并行架构同步重构以提高整体效率。针对共轭梯度法这一用于求解大规模稀疏正则化方程组的迭代方法，构建自适应步长调整机制，把残差下降速率、搜索方向与局部曲率信息进行综合衡量，采用非单调线搜索与回溯限幅相结合的策略来抑制震荡并减少不必要的迭代。同时，把多物理量的预条件器切换为分块对角与代数多重网格混合形式，使高频误差在浅层快速衰减，低频误差在深部得到有效抑制；对交叉梯度与物性约束的权重进行迭代内自适应更新，使方向性正则与耦合项对搜索子空间保持稳定影响，从而在迭代早期加速收敛而在后期维持收敛质量。面向大规模数据的算力组织，设计多层并行计算框架，把三维模型域进行空间分解与任务图调度，依靠GPU对稀疏矩阵向量乘、向量内积与预条件平滑进行加速，并把重力、磁力与电法的正演算子常驻在设备显存以降低主机与设备的数据往返。通信层采用计算与通信重叠的策略，在梯度组装与残差规约之间进行异步流水化，减少全局同步带来的等待时间；精度层构建混合精度求解流程，把单精度用于算子应用与向量运算，双精度用于收敛判据与步长评估，并设置误差监测与动态回退以维持数值稳定。见图1，该迭代流程以数据预处理、预条件初始化、自适应步长决策、并行算子应用与收敛判据评估为主线，辅以分区权重更新与通信计算重叠的支路，使算法在复杂构造区的多物性耦合场景下获得更高的迭代效率与算力可扩展性。

图1 改进后的三维反演算法迭代流程

2.3 多尺度数据融合的反演策略改进

鉴于电法对浅部细节更敏感而重磁观测在区域尺度更稳健，本研究把离散小波分解当作融合枢纽来组织不同物理场的尺度贡献。对重力、磁力与电法观测开展统一分解，把低频近似系数表征深部与广域骨架，把多层细节系数表达断裂边界与岩性突变。从误差传播与分辨率核出发，在小波域构建尺度依赖权重，使电法主导高频细节，使重磁主导低频骨架并稳定深部趋势，同时配合软阈值与边缘保持巩固界面。融合后的多尺度观测嵌入联合反演的数据拟合项，等价于在小波域施加分层约束，使浅层高频约束更强而深部低频更稳健，并把交叉梯度与物性关联在同一尺度网格内对齐以减小物性间尺度错配。迭代过程中，权重随各尺度残差与灵敏度谱自适应更新，重建结果再投影至空间域用于体素模型更新，以在同一优化框架内并行推进细节分辨率与深部一致性。

3 改进方法的地质解释应用与效果评估

3.1 复杂构造区地质模型构建与解释流程

面向研究区的复杂构造格局，本研究把改进的三维联合反演结果当作地质解释的基础数据，围绕断层、岩体与矿化带的空间组合开展一体化建模，并把流程标准化以提高解释的可复用性与一致性。见图2，整体架构由尺度对齐与误差评估、界面抽取与几何建模、岩性域识别与物性赋值、矿化线索识别与靶区圈定、模型一致性检核与不确定性表达五个环节顺序衔接。鉴于多物理量分辨率不一，先把重磁电体素模型在统一栅格内对齐，并把浅层高频与深部低频分别强化与稳固，以减少尺度错配的解释偏差。界面抽取阶段把交叉梯度界面作为边界指示，交叉梯度界面指不同物性梯度同向集中之处，能够刻画共享地质边界，同时引入走向约束维持断裂连续性，并把断裂核部与损伤带分区表达。岩性域识别依靠岩石物性联合分布，把密度、磁化率与电导率的协同变化映射为岩性类别，并把井孔测井与露头测绘线索用于校准。矿化识别从多尺度异常叠合出发，把沿断裂的高密度高磁化率岩体与电导率增高的含流体带进行协同筛选，构建以断裂走向控制的矿化走廊。集成阶段把体素边界转化为曲面与实体，保证断层面与岩体外形在几何上互不矛盾，并开展数据拟合残差与结构逻辑的双重检核。需重点关注的是，不确定性随深度与先验置信度分层表达，把可信度以体素形式投影至三维模型，为后续勘查部署提供可视化与定量化支持。

图 2 复杂构造区地质解释流程架构图

3.2 改进方法的效果验证与对比分析

鉴于江西德兴铜矿东南段断裂叠置与岩性横向突变共存的构造背景，本研究把改进前后联合反演结果与ZK系列钻孔揭示界面进行对照，并把数据拟合、几何一致性以及计算代价三类指标进行量化评估。综合对比表明，交叉梯度与岩性耦合约束在多尺度融合框架中的联动作用，使浅层高频异常得到更清晰的边界表达，同时维持深部低频骨架的稳定，表现为拟合误差下降与交叉验证残差收敛同步发生，断裂核部与岩性界面的深度偏差显著缩小。迭代层面，分块对角与代数多重网格混合预条件配合自适应步长与非单调线搜索，把有效迭代步数与总耗时同步压缩，峰值内存占用随之降低，满足百万体素规模计算的可扩展性需求。结合矿区24口钻孔的路标信息，含磁性斑岩体与沿NE向主断裂的导电带在改进结果中呈现连续配置，钻孔命中率的提升与异常体几何一致性相互印证，靶区圈定更加集中且与已知矿体的空间关系更为协调。关键量化指标的对比如表1所列，其变化幅度与前述机理改进保持一致性，支撑了地质解释可靠性的整体提升。

表1 改进前后反演方法效果对比表

结语

本研究构建了面向复杂构造区的重磁电三维联合反演框架，突出共享地质界面思想，以交叉梯度与岩性物性耦合为主线，配合走向引导与深度衰减形成分级约束体系。在算法层面，通过自适应步长与非单调线搜索、混合预条件与并行加速，显著压缩有效迭代与总耗时；在数据层面，通过小波域多尺度融合，实现电法主导的高频细节与重磁主导的低频骨架协同。德兴矿区实例验证了方法的有效性，拟合误差、界面深度偏差与计算代价均明显下降，钻孔命中率显著提升，联合反演结果更好地支撑了断层几何与岩性分区的综合解释。本文还给出面向地质建模的标准化工作流，通过界面抽取、岩性域识别与不确定性表达，增强了成果的可用性与可复用性。后续工作将面向各向异性与诱导极化等扩展物理场，完善非线性岩石物性关系与时变约束表达，并引入更严格的贝叶斯不确定性量化与先验自适应，以进一步提升复杂构造区深部成像与找矿预测能力。

参考文献

[1]高秀鹤,熊盛青,孙思源,曾昭发,于长春.构建组合模型加权函数并应用于航磁数据三维反演建模（英文）[J].Applied Geophysics,2025,22(02):342-353+556.

[2]张旭,徐得忠,韩媛,任华宁,戚志鹏,李貅.多源地空瞬变电磁一维反演及其在三维模型数据解释中的应用分析[J].西北地质,2025,58(03):33-40.

[3]罗姣,陈炳锦,郭文波,熊伯全,王备战,王伟,陈靖,冯凡,徐毅.基于小波多尺度分析的磁法三维可视化建模研究[J].西北地质,2025,58(03):41-49.

[4]刘晓,郑芳文.基于主轴各向异性的张量CSAMT数据空间域Occam三维反演（英文）[J].Applied Geophysics,2025,22(02):252-263+554.

[5]周官群,赵宏达,张维鑫,吴小平,梁希凯,王宗涛,王亚飞.基于跨巷观测的三维直流电法底板水害精细化探测[J].煤矿安全,2025,56(06):192-201.