基于MMTs的吉林省矿区磁异常线性特征分析
2019-08-14卢昕婷
毛 贺,卢昕婷
(1.吉林省地质调查院,吉林 长春 130102;2.长春工程学院,吉林 长春 130012)
传统的磁异常处理手段包括:磁异常化极、解析延拓、磁源重力异常、磁异常分量换算、加深区域内矿石磁性研究等。这些方法将复杂的磁异常简单化,或者突出磁异常的某一特点,使得异常更加容易解释。这些方法都有自己的特点和应用条件,并且使用都受到一定的限制。比如化极需要地磁场和地质体磁化强度矢量的偏角和倾角信息。一般来说,地质体都具有显著的剩余磁化强度,其和感应磁化强度往往还不是同一方向。在这种情况下,地质体磁化强度矢量的方向可能是未知和不准确的。低磁纬度矿区化极计算同样也不稳定,其会产生平滑但部分错误的异常[1]。解析延拓使用比测量磁场高一个数量级的磁位导数,其异常图分析可能会遗漏一些深度磁源异常引起的低频信息。解析延拓产生的异常还取决于地质体磁化矢量和周围地磁场的方向。磁源重力异常的问题和类似化极处理。另外,选择磁位而不是测量磁场进行解释,会降低分辨率,特别是针对浅部位置矿体的异常,因为计算磁位会模糊高频信息。
在文献[2]中,提出了级磁变换MMTs。MMTs包括TMA本身Ta,以及基于TMA的四种变换:TMA梯度的模R、TMALaplacian算子结果L、E为TMA平方的Laplacian算子结果再开方取二分之一、Q为TMA和TMALaplacian算子结果的乘积的平方根。上述四种变换结果具有非负的分布,对磁化强度矢量方向的敏感程度比TMA自身低得多,并且它们对磁异常源能显示出比测量的磁异常场更高的中心性。本文使用MMTs对吉林省某矿区的磁异常数据进行处理分析,MMTs得到的异常更接近磁异常源的真实水平位置。
1 MMTs的方法原理
2 试验与分析
本文使用MMTs对吉林省某矿区的磁异常数据进行处理,测区ΔT平面图如图1a所示。
由图1a可以看出,磁异常主要集中在该研究矿区的南部,异常走向主要有北西向和北东向,个别异常沿近东西向展布,多呈团块状或带状分布;测区的中部偏南主要是低缓的负异常,变化梯度小;测区的北部主要以大面积低缓正异常为主,局部叠加低缓负异常。
MMTs处理结果如图1b至图1d所示。在Ta等值线平面图中,正异常向北侧移动,异常的线性特征相对突出、明显。Q、E、R等值线图的异常分布与Ta相似,异常范围边界呈递进收敛关系,位置更加清晰、明确,也更接近磁异常源的真实水平位置。L等值线图中的异常解比较分散,难以准确的确定异常矿体的边界位置,但其对浅部异常矿体有一点的指示作用。
图1
3 结论
本文使用MMTs对吉林省某矿区磁测数据进行处理,取得了较好的处理效果。
MMTs和传统磁异常处理方法相比较具有以下优势:MMTs对某些噪声的放大作用小于解析延拓;在低磁纬度矿区比化极稳定。此外,MMTs算法简单,程序实现方便,具有很高的计算效率。