基于磁异常模量的三维磁化率成像研究应用.pdfVIP

·696· 中国地球物理2011 (19)地质调查和矿产勘查地球物理 基于磁异常模量的三维磁化率成像研究及应用 3D basedonthe magneticsusceptibilityimaging of anomalies ampfitudemagnetic 刘圣博1‘陈超2崔放1 Cui Liu Cha02 Fan91 Shengb01+Chen 1．中国地质调查局武汉地质调查中心 武汉430205； 2．中国地质大学(武汉)地球物理与空间信息学院武汉430074 1．引言 以往的磁异常资料解释一般假设地下磁性体的磁化方向与地磁场一致，它所具有的剩磁和退磁效应很 小以致不影响磁性体的磁化方向，然而在许多实际问题中地下磁性体具有很强的剩余磁化强度。由于所有 的磁异常分量以及AT的特征都强烈依赖于磁性体的磁化方向，当磁性体的剩磁很强且使实际磁化方向偏 离原地磁场方向时，就不满足前面的假设条件，由此将给磁异常解释工作带来困难。理论上可以证明，磁 异常模量基本不依赖或弱依赖磁性体的磁化方向。因此，利用磁异常模量进行磁法解释，能压制强剩磁的 影响。同时，三维磁化率成像技术是一种定量的磁法解释方法，能较形象的体现出磁性异常体的分布特征。 2．磁异常模量的特征 通常磁测数据是A丁，可以在频率域将△丁数据转换成磁场的三分量Ek、H。和乙。如果是三分量磁力 测量的话，可以直接利用(1)式计算磁异常模量死。理论上可以证明，对于二度体而言磁异常模量形态 完全不依赖于异常体的磁化方向，而三度体的磁异常模量形态微弱地受异常体磁化方向影响。磁异常模量 具有这个特征，是由于其相对磁异常△r而言融合了磁异常三分量的信息。 瓦=(^乙2+^b2+z。2)。／2 (1) 3．基于磁异常模量的三维磁化率成像原理 Occam最小构造反演方法，其过程是将地下三维模型空间剖分成几何大小相等或者不等的直立长方 体，并假设每个模型单元是均质均匀磁化的。反演时，假设地下模型单元的磁化方向为正常场的方向，然 后基于观测的磁异常数据实现重构地下每个单元磁化率的数值，这便实现了磁化率三维成像。磁化率成像 是一个反演问题，利用最小模型和光滑模型进行约束以降低反演中的自由度获得唯一解。构建的反演目标 函数中有数据拟合目标函数和模型目标函数两项，之间的权重关系依据Tikhonov曲线的原则来确定。文中 采用共轭梯度法进行求解。对磁化率范围约束反演时，在迭代过程中采用非线性变换法逐步实现对磁化率 进行约束。磁异常模量反演的目标函数形式如下式： 矽(七)=8黟0(d咖一瓦∥耐(七))0‘+∥』阡乙(七一t订)旷 (2) 4．结论 磁异常模量包含了磁异常三分量的信息，弱化了磁化方向的影响，具有弱依赖异常体磁化方向的特性。 因此，基于磁异常模量的磁化率成像能压制强剩磁的影响。虽然磁异常模量融合了磁异常三个方向上分量 信息，但由于它是标量，只体现了各分量信息中的幅值信息，而缺失了相位信息。所以磁异常模量又会模 糊部分信息，导致磁化率成像结果垂向上分辨力严重受损。同时，磁异常模量数据是通过磁测△r数据转 换获得，若实际磁性体存在剩磁或退磁时，其转换结果自然与实际不同。所以，利用三分量磁测数据进行 磁异常模量解释或反演，才能充分体现该方法的特点。 参考文献 of S ShearerS．Three—dimensionalinversion datainthe ofremancnt thesis，Col