基于磁法数据的边缘识别方法研究——以内蒙古某地区为例

安百州， 冯志民， 李宁生， 张代磊， 白亚东1b,

(1.中国地质大学 a.地球科学与资源学院，b.地球物理与信息技术学院, 北京 100083；2.宁夏地球物理地球化学勘查院， 银川 750001；3.吉林大学 地球探测科学与技术学院， 长春 130021；)



基于磁法数据的边缘识别方法研究
——以内蒙古某地区为例

安百州1a，2， 冯志民2， 李宁生2， 张代磊3， 白亚东1b,2

(1.中国地质大学 a.地球科学与资源学院，b.地球物理与信息技术学院, 北京 100083；2.宁夏地球物理地球化学勘查院， 银川 750001；3.吉林大学 地球探测科学与技术学院， 长春 130021；)

为了准确划分内蒙古某地区的地质体边界和断裂位置，采用解析信号振幅和斜导数法对化极后的高精度磁法数据进行了处理，并将其综合解释结果与利用垂向一阶导数与总水平梯度法的综合解释结果进行了对比分析。通过解析信号振幅和斜导数法综合解释了7条断层和2条岩体分界线，解释结果与实际地质资料吻合较好，而通过垂向一阶导数与总水平梯度法综合解释了6条边界线，与实际资料具有一定误差。说明解析信号振幅和斜导数在识别一些地质体边界和确定某些断裂位置时,具有一定的优势。

高精度磁法； 边缘识别； 解析信号振幅； 斜导数

0 引言

地质体边缘和断裂构造线由于密度或磁性的差异常会导致重、磁异常变化率较大，因此所有的边缘识别方法均利用这一特点进行设计[1]。常用的边缘识别方法有数理统计类和数值计算两类方法。其中数理统计的方法有小子域滤波法和归一化标准差法等；数值计算的方法有垂向导数(VDR)、总水平导数(THDR)、解析信号振幅(ASM)和斜导数等方法[1-2]，①小子域滤波法是基于滑动平均法原理进行改进的一种低通滤波方法，但是常规的小子域滤波法容易引起异常的畸变，造成异常梯级带形态的扭曲[3-4]；②归一化标准差方法是计算一个滑动窗口内垂向坐标方向一阶导数的标准差与三个坐标方向一阶导数标准差之和的比值，并以此来识别边缘位置，它的缺点是难以比较边缘识别结果的精度[1,5]；③垂向导数方法是利用零值位置确定边界位置，但是导数的阶次和地质体的埋深都会影响垂向导数所确定的边缘位置[1,6-7]；④总水平导数是利用极大值位置来进行边界划分的，但是倾斜边界、地质体的埋深和地形起伏都会对解释结果有一定影响[8-10]；⑤解析信号振幅是利用极大值位置来确定地质体的边缘位置，其优点受磁异常分量和磁化方向的影响最小，缺点是横向分辨率能力较低[11-12]；⑥斜导数的实质是垂向导数和总水平导数的比值，由于它能很好地平衡高幅值异常和低幅值异常，所以它受地质体倾斜边界的影响比较小，其缺点是当解析信号振幅等于“0”时存在“解析奇点”，故其数值计算稳定性较差[1,13-15]。为了研究内蒙古某地区的地质体和断裂构造线位置，分别采用垂向一阶导数、总水平导数、解析信号振幅和斜导数等方法对研究区的1∶5 000 高精度磁法勘探数据进行处理，并将其解释结果进行对比分析，以此来初步确定研究区的地质体和断裂的分布情况，并作出相应的地质解释。

1 方法原理

1.1 垂向一阶导数

垂向一阶导数是利用零值位置确定地质体的边缘位置，它可直接适用于重力异常，对磁力异常而言，需先转换成假重力异常或化极磁力异常。它是沿垂直方向求一次导数的数据处理方法，它对磁场的高频成分有突出和放大作用，侧重于浅层地质体的局部异常而压制深层区的背景场，其计算公式为式(1)。

(1)

垂向一阶导数方法研究较早，方法较成熟，应用广泛。由于垂向一阶导数属高频放大因子，因此在计算中需压制高频干扰，常用的方法有向上延拓、维纳滤波、正则化滤波等。

1.2 总水平导数

总水平导数是利用其极大值位置来确定地质体的边缘位置，对磁力异常需化极后才能使用。其原理是先计算两个水平方向一阶导数，然后再计算总水平导数，计算公式为：

(2)

其中：ΔT0为某一平面上的数据体。

它在边缘识别技术中应用率较高，其计算方式有空间域和频率域两种，空间域是用中心差分格式计算水平方向的一阶导数，其优点是结果比较稳定，计算量较小，但计算结果易受网格间距的影响；频率域是通过频谱变换来计算水平一阶导数的，其计算结果稳定性差，计算量较大，但其受网格间距大小的影响较小。

1.3 解析信号振幅与斜导数

解析信号最初来源于信号分析。对于离散信号，不方便分析频率、相位等属性，所以需要构建一种信号的载体，再用这个载体来分析。解析信号就是其中的一种载体，在测线为剖面的情况(对应二维地质体)中，重力或磁场强度导数水平分量的希尔伯特变换就是重力或磁场强度导数的垂直分量[16-18]。所以其计算公式为式(3)。

(3)

斜导数的计算公式为式(4)。

(4)

由于解析信号受磁异常分量和磁化方向的影响较小，所以它在一些地质体的边界识别中取得了良好的应用效果[19-23]。在边界垂直的情况下，解析信号振幅的极大值对应边界，斜导数的零值线对应边界；边界倾斜时(如倾斜的板状体或倾斜的断层)，可以用斜导数的水平导数来确定上顶面的边界[24-28]。理论上，对于磁法，在二维情况下，利用解析信号(振幅ASM与斜导数TD)进行边界识别的优点在于不依赖于磁倾角和磁偏角，而水平导数和垂向一阶导数要化极后才能正确指示边界。所以解析信号振幅和斜导数在确定某些地质体边界时，比水平导数和垂向一阶导数具有更好的效果。

2 应用实例

2.1 研究区地质概况

研究区的大地构造属华北地台的阴山断隆东段北部，中元古界地层出露为阿拉善群和查尔泰群，其中阿拉善群可以分为迭布斯格下亚组第一岩段与迭布斯格下亚组第二岩段，迭布斯格下亚组第一岩段的岩性主要为斜长角闪片麻岩、斜长片麻岩、磁铁石英岩和大理岩等，迭布斯格下亚组第二岩段的岩性主要为斜长角闪片麻岩、透辉大理岩和磁铁石英岩等；查尔泰群岩性主要为千枚岩、粉砂岩和泥灰岩。新生界地层出露第三系渐新统和第四系，第三系渐新统的岩性主要为泥岩、细砂岩和砂砾岩，第四系主要为洪积物、砂砾石和砂土层。该区的岩浆活动主要为华力西期的花岗岩，主要出露于研究区的西北角，岩性主要为黑云母花岗岩，局部有高岭土化和绢云母化蚀变。下部所含磁铁石英岩矿层在研究区均有分布。

2.2 地球物理特征

本次工作在研究区展开了1:5 000的地面高精度磁测，并制作了相应的ΔT平面等值线图(图1)和相应岩性的磁性参数(表1),从图1中可以看出，研究区磁异常主要以北东向为主，这反映了本区的主构造方向为北东向，这与地质图所显示构造方向一致。本次研究区主要分为三个区域：①北部基本为分布于低值背景或负值背景上的串珠状高值磁异常，结合地质图可知，该低值异常区可能为迭布斯格下亚组第二岩段的斜长角闪片麻岩和透辉大理岩引起，而出现的串珠状高值磁异常应为迭布斯格下亚组第二岩段的磁铁石英岩引起，因此应重点勘查北部的串珠状高值磁异常；②中部为北东东向等轴状异常上叠加串珠状次级异常，该异常可能为迭布斯格下亚组第一岩段和迭布斯格下亚组第二岩段的磁铁石英岩引起，从地质图上可以看出在中部存在着北东东向的断层和岩性分界线，且它们的位置与串珠状的异常的位置重合，说明矿体可能与北东东向的断裂有关；③南部为东南与西南两平稳正值异常，磁场变化较平稳，反应了基底岩石磁性弱，盖层较厚，应与矿体关系不大。

图1 研究区ΔT平面等值线图

2.3 垂向导数与总水平梯度模解释

垂向导数与总水平导数的计算一般在频率域中通过频谱变换的方法实现，由于垂向导数与总水平导数的频率响应因子属高频放大因子，在计算时需要压制高频干扰的影响。这里通过不同高度的向上延拓，发现向上延拓30m能取得较好的效果。因ΔT磁异常需经化极后才能应用垂向导数和总水平导数来识别磁性体的边缘位置。故对ΔT磁异常先做上延后化极，再用化极磁异常分别作垂向一阶导数(图2(a))与总水平导数(图2(b))，最后进行解释。解释结果如图2(c)所示，其中红色线段为解释的构造线，黑色线段为填图的成果图。

表1 研究区岩矿石磁性参数统计表

由图2可看出，垂向一阶导数图的零值线与总水平导数的极值所确定的地质体边界，与地质图对应情况总体效果较好，但是对局部的一些构造或岩体分界线反映不够清楚，可能是由于受地质体边界的影响较大造成的。从图2可以看出，华力西期花岗岩与迭布斯格下亚组第二岩段边界位置与地质图基本吻合，且在迭布斯格下亚组第二岩段中解释的断层与地质图中北东向断层比较吻合，这说明解释效果较好。从图1可以看出，迭布斯格下亚组第二岩段显示磁性分布不均匀，前人资料表明，当断裂带为一构造破碎带时，岩层或岩体的磁场因为断裂破碎带的影响而降低，磁场在断裂带分布范围内显示降低的特点，往往形成磁场的间断。据此认为迭布斯格下亚组第一岩段分布范围应多发育次级断层，这在我们所解释的断裂图上也有所反映。但是该解释图(图2(c))对迭布斯格下亚组第一岩段与迭布斯格下亚组第二岩段的岩性分界不明显，且对西北角上的华力西期花岗岩中断裂没有反应，说明该解释图(图2(c))的精确度还有待提高。

图2 研究区垂向导数与总水平梯度模解释图Fig.2 The map of veritical derivative and the map of total horizontal derivative for the reasearch area

图3 解析信号振幅和斜导数解释图

2.4 解析信号振幅和斜导数解释

解析信号振幅(ASM)是利用极大值位置来确定地质体的边缘位置，其优点是受磁异常分量和磁化方向的影响小。斜导数(TD)的实质是垂向导数和总水平导数的比值，由于斜导数为一阶导数的比值，所以它能将高幅值异常和低幅值异常很好地平衡，从而起到增强边缘的效果。通过运用解析信号振幅和斜导数对研究区的化极磁异常处理后得到的等值线图如图3(a)、3(b)所示，将两者的等值线进行综合解释，得到的解释图如图3(c)所示，其中红色线段为解释的构造线，黑色线段为地质填图结果。从图3(c)中可以看出，迭布斯格下亚组第一岩段与迭布斯格下亚组第二岩段的岩性分界在图中有较好地反映，华力西期的玄武岩与迭布斯格下亚组第二岩段的岩性分界也在图中表现明显，而且对迭布斯格下亚组第一岩段与第四系的岩性分界有一定的显示，对西北角上的华力西期花岗岩中的断层也具有反映。这说明了斜导数解释与地质图吻合良好，具有更好的解释效果。

2.5 解析信号振幅、斜导数与垂向导数、总水平 导数法对比

通过对比解析信号振幅、斜导数与垂向导数、总水平导数，我们可以发现，地质图中西北角华力西期花岗岩的北东向断层在解析信号振幅和斜导数中反应较好，而在垂向一阶导数和总水平导数中均反应不明显或无反应；研究区中部的迭布斯格下亚组第一岩段与迭布斯格下亚组第二岩段的边界在解析信号振幅、斜导数中反映更加明显，而垂向导数和总水平导数对该条岩性分界线的解释效果较差；对于迭布斯格下亚组第一岩段中的次级断裂，斜导数、垂向导数和总水平导数均有一定程度的显示；迭布斯格下亚组第一岩段与第四系的岩性分界在解析信号振幅、斜导数解释中反映较明显，而在垂向导数和总水平导数中无反应。综上所述，解析信号振幅、斜导数比垂向导数、总水平导数具有更好的解释效果，与地质资料更为吻合。

2.6 地质解释

根据解析信号振幅和斜导数的解释结果，对该区构造进行划分，推断断层7条(F1、F2、F3、F4、F5、F7、F8)，其中F1、F2和F3为北东向断层，F4、F5、F7和F8为北西向断层；推断岩体分界线2条(F6、F9)。从图3(a)可以看出，本区的异常极大值基本接近于断层两侧位置，显示出高磁由断层控制的特点，认为断层为高磁的形成提供了条件。其中F1为华力西期花岗岩中的北东向断层，F2断层控制着华力西期花岗岩与迭布斯格下亚组第二岩段的分布，该断层西北为华力西期花岗岩，从图1可看出，整套花岗岩整体显示负值，磁场变化平稳，说明花岗岩基本无磁或者弱磁。沿该断层两侧展布数串珠状高值，其极值分别达到542 nT、655 nT、1 005 nT、1 084 nT、524 nT。经过不同高度的上延，当上延高度达到100 m时，这些次级异常基本全部消失，根据经验公式，即三度体上延高度的三分之一为异常源埋藏深度，推测高磁体深度在地表下30 m左右，说明异常由浅部高磁体引起。综合地质图分析，认为沿断裂分布的串珠状异常可能与磁铁矿体有关，后经钻探工程验证，在地下发现了磁铁矿，说明解释效果较好。F3、F4和F5断层以不同角度错段迭布斯格下亚组第二岩段，其中F4和F5两断裂北西部均出露华力西期花岗岩体，沿F3、F5展布的串珠状高值，其极值可达到 1 394 nT、1 626 nT。经上延后，F3西南端及至中断的异常全部消失，说明F3西南端及至中断高磁由浅源异常引起。F3东段异常经上延后异常依然明显。综合地质图，以上异常可能与迭布斯格下亚组第二岩段磁铁石英岩有关。沿F5断裂展布的带状异常，走向北西，由两个次级异常组成，极值分别达到669 nT、1 153 nT，上延后异常明显，综合地质图，异常可能由矿体引起。F3、F5断裂交夹位置展布一近南北向的等轴状高值异常，该异常形态较不甚规则，极值略低于沿断层展布的极值，可能为中磁性变质岩引起，是否与含矿体岩石有关需进一步工作验证。F6为迭布斯格下亚组第一岩段与迭布斯格下亚组第二岩段的岩性分界线，从图3(c)中可以看出与地质图的吻合情况较好；F7和F8为迭布斯格下亚组第一岩段的北西向断裂，沿这两条断裂均有高值分布，其最高值可达1 788 nT；在研究区东南角的F9为迭布斯格下亚组第一岩段与第四系的分界线，该解释线与实际地质比较吻合，解释效果较好。

3 结论与建议

解析信号振幅和斜导数可有效识别断裂和地质体边界，与垂向导数、总水平梯度模等方法相比具有较明显的优势，这两种方法在解决一些倾角不同、埋深不同的多源场地质体边界和产状，以及划分一些大地构造单元，确定某些断裂构造带的位置，进行构造分区等问题上具有一定的优势。

在研究区的F2断裂上，经过钻孔验证发现了磁铁矿，可以初步判断北东向的断裂为控矿断裂。矿体可能是由于华力西期花岗岩侵入带来的岩浆热液，使得研究区的含铁硅质岩发生变质作用，从而形成含铁石英岩矿床。

[1] 王万银，邱之云，杨永，等.位场边缘识别方法研究进展[J].地球物理学进展，2010，25(1)：196-210. WANG W Y,QIU Z Y,YANG Y,et al. Some advances in the edge recognition of the potential field.Progress in Geophysics,2010,25(1):196-210.(In Chinese)

[2] 郭华，吴燕冈，高铁. 重力斜导数方法在时间域中的理论模型与研究[J]. 吉林大学学报：地球科学版，2006,36:9-14. GUO H,WU Y G,GAO T. The research of theories model in time area with the method of tilt derivative in gravity.Journal of Jilin University:Earth Science Edition,2006,36:9-14.(In Chinese)

[3] 张凤旭，张凤琴，刘财，等. 裂构造精细解释技术—三方向小子域滤波法[J]. 物理学报，2007，50(5)：1543-1550. ZHANG F X ,ZHANG F Q ,LIU C,et al. A technique for elaborate explanation of faulted structures : three-directional small subdomain filtering. Chinese Journal of Geophysics,2007 ,50 (5) :1543-1550.(In Chinese)

[4] 马国庆，杜晓娟，李丽丽.优化小子域滤波方法研究及其应用[J]，石油地球物理勘探，2013,48(4)：658-662. MA G Q，DU X J,LI L L. An optimized small sub-domain filtering method,Oil Geophysical Prospecting.2013,48(4):658-662.(In Chinese)

[5] COOPER G.RJ, COWAN DR. Edge enhancement of potential—field data using normalized statistics[J]. Geophysics,2008,73(3):H1-H4.

[6] 崔瑞华，谷社峰，李成立.位场垂向导数的稳定算法[J].物探化探计算技术，2009，31(5)：426-430. CUI R H,GU SH F,LI CH L.A stable algorithm for calculating the vertical derivatives of potential field,Computing Techniques for Geophysical and Geochemical Exploration,2009,31(5):426-430.(In Chinese)

[7] 韩利，袁炳强，张春灌，等.位场各阶垂向导数的ISVD算法及其应用[J].物探化探计算技术，2013，35(6)：668-674. HAN L,YUAN B Q,ZHANG CH G,et al. Isvd algorithm and its application for calculating each order vertical derivative of potential field, Computing Techniques for Geophysical and Geochemical Exploration ,2013，35(6)：668-674.(In Chinese)

[8] 钟清.区域重力资料在地质填图中的边界定位问题研究[D].北京：中国地质大学，2006. ZHONG Q.A study on Bounding problem using Gravity data in the Geological Mapping.Beijing:China University of Geosciences，2006.(In Chinese)

[9] 钟清，孟小红，刘士毅.重力资料地位地质体边界问题的探讨[J].物探化探计算技术，2007，29(增刊)：35-38. ZHONG Q,MENG X H,LIU SH Y. A discussion of Bounding Mapping by Gravity data. Computing Techniques for Geophysical and Geochemical Exploration，2007，29(supp1)：35-38.(In Chinese)

[10]WANG W Y, PAN Y, QIU Z Y. A new edge recognition technology based on the normalized for potential field data[J].Applied Geophysics,2009,6(3):226-233.

[11]HSU S K, SIBUET J C, SHYU C T. High-resolution detection of geologic boundaries from pitentital-field anomalies: An enhanced analytic signal technique[J]. Geophysics, 1996,61(2):373-386.

[12]DEBEGLIA N, CORPEL J. Automatic 3-D interpretation of potential field data using analytic signal derivatives[J]. Geophysics,1997,62(1):87-96.

[13]王明, 郭志宏, 何辉, 等.基于反双曲正切的位场边界识别技术[J].物探与化探，2013，37 (4): 655-663. WANG M,GUO ZH H,HE H,et al. Edge detection of potential field data using inverse hyperbolic tangent.Geophysical and Geochemical Exploration，2013，37 (4): 655-663.(In Chinese)

[14]王明, 郭志宏, 骆遥. Tilt-Euler方法在位场数据处理及解释中的应用[J].物探与化探. 2012，36 (1): 126-132. WANG M,GUO ZH H,LUO Y. The application of Tilt-Euler deconvolution method to potential field data processing and interpretation. Geophysical and Geochemical Exploration. 2012，36 (1): 126-132.(In Chinese)

[15]李卓聪.斜导数方法研究及其在道郎和都格矿区重力资料解释中的应用[D].长春：吉林大学，2009. LI ZH C. Reaserching of Tilt derivative and application in Daolangheduge mine.Changchun:Jilin University，2009.(In Chinese)

[16]胡中栋, 余钦范,楼海. 三维解析信号法[J]. 物探化探计算技术, 1995, 17(3)：36-42. HU ZH D,YU Q F,LOU H. 3-D analytic signal method. Computing Techniques for Geophysical and Geochemical Exploration, 1995, 17(3)：36-42. (In Chinese)

[17]罗利春. 用希尔伯特变换构造解析信号进行时频分析[J]. 航天电子对抗, 2003 (3): 26-29. LUO L CH. Time-frequency analysis of analytic signal based on Hilbert transform. Aerospace Electronic Warfare, 2003 (3): 26-29.(In Chinese)

[18]骆遥, 王明, 罗锋, 等. 重磁场二维希尔伯特变换—直接解析信号解释方法[J]. 地球物理学报, 2011, 54(7): 1912-1920. LUO Y,WANG M,LUO F,et al. Direct analytic signal interpretation of potential field data using 2-D Hilbert transform.Chinese Journal of Geophysics, 2011, 54(7): 1912-1920.(In Chinese)

[19]刘义波, 曾昭发, 陈雄, 等. 解析信号和水平导数法超贫铁矿磁探测—以内蒙古凉城铁矿为例[J]. 世界地质, 2013, 32(4): 839-846. LIU Y B,ZENG SH F,CHEN X,et al. Method of analytic signal and horizontal derivative in magnetic prospecting lean iron deposit:a case study in Liangcheng area,Inner Mongolia. Global Geology, 2013, 32(4): 839-846.(In Chinese)

[20]秦葆瑚. 用解析信号计算倾角异常[J]. 物探化探计算技术, 1998, 20(3): 261-266. QIN B H. Calculation of Tilt by analytic signal. Computing Techniques for Geophysical and Geochemical Exploration, 1998, 20(3): 261-266.(In Chinese)

[21]张季生. 用三维解析信号技术处理化南航磁异常[J]. 物探与化探, 2000, 24(3): 190-196. ZHANG J SH. The application of 3-D analytic signal technique to the treatment of aeromagnetic anomalies in south China. Geophysical and Geochemical Exploration, 2000, 24(3): 190-196.(In Chinese)

[22]王万银. 位场解析信号振幅极值位置空间变化规律研究[J]. 地球物理学报, 2012, 55(4): 1288-1299. WANG W Y. Spatial variation law of the extreme value position of analytical signal amplitude for potential field data. Chinese Journal of Geophysics, 2012, 55(4): 1288-1299.(In Chinese)

[23]张季生. 增强解析信号技术的应用条件[J]. 物探与化探, 1999, 23(4): 296-300. ZHANG J SH. The application conditions of the technique for the enhancement of analytic signals. Geophysical and Geochemical Exploration, 1999, 23(4): 296-300.(In Chinese)

[24]郭华, 于长春, 吴燕冈. 改进的斜导数方法及应用[J]. 物探与化探, 2009, 33(2): 212-216. GUO H,YU CH CH,WU Y G. An improved Tilt derivative method and its application. Geophysical and Geochemical Exploration, 2009, 33(2): 212-216.(In Chinese)

[25]孙晖, 刘万崧, 王洪昌,等. 基于重磁场特征的敦化盆地构造格架研究[J].世界地质,2012,31(1):171-178. SUN H,LIU W S,WANG H CH, et al. Tectonic framework of Dunhua Basin based on gravity and magnetic characteristics.Global Geology,2012,31(1):171-178.(In Chinese)

[26]郭华, 王明, 谢汝宽, 等. 向上延拓对斜导数方法处理航磁异常的影响[J]. 物探与化探, 2012, 36(1): 137-143. GUO H,WANG M,XIE R K,et al. The effects of the upward continuation on the application of Tilt derivative method to aeromagnetic data processing. Geophysical and Geochemical Exploration, 2012, 36(1): 137-143.(In Chinese)

[27]戴明刚, 曲寿利. 位场梯度模法及其在碳酸盐岩地区断裂识别中的应用[J]. 地球物理学进展, 2009 (3): 951-958. DAI M G,QU S L. Determination of faults using gradient modules of potential field anomalies and its application to carbonate rock areas of south China. Progress in Geophysics, 2009 (3): 951-958.(In Chinese)

Case study: A research on edge detection methods based on magneticdata in Inner Mongolia

AN Bai-zhou1a,2， FENG Zhi-min2, LI Ning-sheng2， ZHANG Dai-lei3， BAI Ya-dong1b,2

(1.China University of Geosciences a.School of Earth Sciences and Resources，b.School of Geophysics and Information Technology，Beijing 100083,China；2.Ningxia Geophysical and Geochemical Exploration Institute，Yinchuan 750001,China；3. The College of Geoexploration Science and technology，Jilin universiy，Changchun 130021,China)

To recognize the location of geology body and faults in someplace in Inner Mongolia, we adopted analytical single amplitude and tilt derivative method to process high resolution magnetic data, and the data have been reduced to pole. Then we adopted vertical derivative and total horizontal gradient method to process the same data. The result shows that by analytical single amplitude and tilt derivative method, 7 faults and 2 boundaries of rock mass can be found, and for derivative and total horizontal gradient method, only 6 boundaries can be found. The former matches better with geology data. It shows that analytical single amplitude and tilt derivative method have some advantages when it was used to recognize geology bodies and locate the position of faults.

high resolution magnetic; edge detection; analytic signal amplitude; tilt derivative

2014-11-03改回日期：2015-01-18

宁夏回族自治区地勘基金项目(2014-地勘基金07)

安百州(1990-)，男，硕士，主要从事重磁资料的应用研究工作， E-mail：941868522@qq.com。

1001-1749(2015)05-0599-07

P 631.2

A

10.3969/j.issn.1001-1749.2015.05.09