基于方向tilt-Euler的三维磁数据快速反演

罗潇， 王彦国， 葛坤朋， 邓居智， 杨亚新

1 东华理工大学地球物理与测控技术学院， 南昌 330013 2 东华理工大学核资源与环境国家重点实验室， 南昌 330013 3 核工业230研究所， 长沙 410011

0 引言

磁法数据处理与解释是磁法勘探工作的核心环节，选取适当的方法技术是获得良好地质解释的关键.受磁化方向影响，磁异常往往与地质体分布之间存在一定偏差，直接进行异常解释难度较大.化极处理可以将感应磁化磁异常转换为垂直磁化异常，在一定程度上有利于异常的解释和推断.另外，许多数据处理方法建立在化极磁异常基础上，如归一化标准差(Cooper and Cowan, 2008)和 Tilt-depth法(Salem et al., 2007).然而，当磁化角度变化较大，研究区位于低纬度，或工区存在强剩磁时，化极很难获得满意的效果.解析信号是磁异常解释的一类有效方法(Bastani and Pedersen, 2001; Srivastava and Agarwal, 2010; Cooper, 2015)，二维解析信号不受磁化方向影响(Atchuta et al., 1981)，三维磁异常解析信号受磁化方向影响较小，但仍不能无视斜磁化的影响(Huang and Guan, 1998; Li, 2006).Stavrev和Gerovsk(2000)在磁三分量和梯度张量基础上给出了R、L、E、Q四个模量转换量，并证实了这些转换量在三维中对磁化方向敏感度较小，二维中不受磁化方向影响.Verduzco等(2004)在tilt梯度的基础上提出了tilt梯度水平导数模，同样在二维磁异常处理中不受斜磁化影响，在三维中受磁化角度影响较小，也被广泛运用于磁异常识别中，不过该方法使用了磁异常二阶导数，易受噪声干扰影响，同时还会产生假异常.Wijns等(2005)提出了Theta map，这种方法可以有效识别低纬度磁异常的边界信息，但在中高纬度地区的……