地球物理学报  2019, Vol. 62 Issue (9): 3580-3590   PDF    
全张量磁梯度数据的斜导数特征值边界识别方法研究
郭华1,2, 韩松1, 郑强3, 梁韧1, 段然1, 刘浩军1, 李西子3     
1. 中国自然资源航空物探遥感中心, 北京 100083;
2. 吉林大学, 长春 130026;
3. 中国地质大学(北京), 北京 100083
摘要:本文提出了一种全张量磁梯度数据的斜导数边界识别方法,该方法利用全张量磁梯度数据定义了北向斜导数、东向斜导数和垂向斜导数.模型试验表明,北向和东向斜导数无法有效识别磁性异常体边界,异常形态复杂虚假异常较多,而垂向斜导数在高纬度斜磁化或者垂直磁化条件下可以清晰而准确得识别地质体的边界,具有一定的深度识别能力.与现有磁张量梯度的边界识别方法相比,识别效果好,可以有效均衡不同深度地质体的响应.将上述方法应用于大兴安岭地区实测航磁三分量数据转换得到的全张量数据,获得了研究区浅部磁性体的分布,并与三分量数据进行联合解释,取得良好的地质效果.
关键词: 航磁三分量      全张量磁梯度      斜导数      边界识别     
Oblique derivative edge detection method research on full tensor magnetic gradient data
GUO Hua1,2, HAN Song1, ZHENG Qiang3, LIANG Ren1, DUAN Ran1, LIU HaoJun1, LI XiZi3     
1. China Aero Geophysical Survey and Remote Sensing Center for Natural Resources, Beijing 100083, China;
2. Jilin University, Changchun 130026, China;
3. China University of Geoscience(Beijing), Beijing 100083, China
Abstract: A oblique derivative edge detection method for full tensor magnetic gradient data is presented in this paper. This method uses full tensor magnetic gradient data to define northward oblique derivative, eastward oblique derivative and vertical oblique derivative. Model tests show that the derivatives of the north and east oblique derivatives cannot effectively determine the boundaries of magnetic anomalies, the anomaly pattern is complicated and there are many false anomalies. While the vertical derivatives can clearly and accurately determine the boundaries of geological bodies under high latitude oblique magnetization or vertical magnetization conditions. Compared with the existing edge detection methods of full tensor magnetic gradient data, the recognition effect is better, and the method can effectively balance the response of geological bodies with different depths. The above method is applied to the full tensor data obtained from the three-component aeromagnetic data conversion in Daxing'anling area, and the distribution of the shallow magnetic bodies in the study area is obtained. Good geological results are obtained with jointly interpretation of the three-component data.
Keywords: Three component of aeromagnetic data    Full tensor magnetic gradient    Oblique derivative    Edge detection    
0 引言

在磁法勘探中,导数被认为是确定地下异常体边界的有效方法.然而,当地下存在多个埋深不同的磁性源时,深层磁性源边界识别难度较大.因为深层磁性源受浅层磁性源的影响,异常幅度、形态会很不明显,梯度以及二阶导振幅同样会受到影响.为了解决埋深不同、倾角不同的多源场物体边界的问题,Miller和Singh(1994)首次给出了Tilt梯度的定义,这种方法求得的梯度值对于场源的深度是不敏感的,不受场源埋深的限制,可以减弱多源地质体在识别边界问题上相互影响,可以很好的探测出埋深不同的多源场物体的边界.由于斜导数方法在判断倾角不是0或π/2的地质体边界时存在一定的缺陷,Verduzco等(2004)提出了斜导数水平导数来解决这一问题.王想和李桐林(2004)理论研究了斜导数方法在重力勘探领域的应用,明确了斜导数水平导数在异常源边界的探测的有效性.Salem等(2007)提出了Tilt-depth概念,为判断地质体边界及埋深提供了有利的依据.Oruç和Keskinsezer(2008)将斜导数的方法理论应用到了全张量重力梯度数据中,在分析地质体的边界中获得了较好的效果.郭华(2009b)提出了一种非绝对值斜导数水平导数方法,在判断断裂构造的倾向上具有较好的作用.郭华(2009a, 2012)利用Tilt-depth方法,并分析了地质体埋深及与地质体边界位置之间的关联关系.Lahti和Karinen(2010)将斜导数理论和向上延拓理论进行综合运用,进而研究地质体的边界.Ferreira等(2013)将斜导数方法中的位场数据转变为水平梯度数据,达到了分析交叉地质体边界的效果.

虽然上述文章都详细地讨论了斜导数方法及其水平导数在判断地质体的边界、产状和埋深方面的效果,但斜导数方法如何应用于全张量磁梯度数据未见文献讨论.而全张量磁梯度数据是磁矢量在不同方向的导数,可以更精确得描述磁性异常体的位置分布及形态特征,且具有较高的水平分辨率(Schmidt and Clark, 2000; 王林飞, 2016).现有的全张量磁梯度数据的边界识别技术(Schmidt and Clark, 2006; Oruç,2010ab)在实际应用中存在限定条件过多,识别效果不佳的问题.本文将针对全张量磁梯度数据,开展斜导数特征值边界识别技术研究,并将该技术应用于实测航磁三分量数据转换得到的全张量磁梯度数据,圈定研究区域内异常体的边界,取得了较好的应用效果.

1 全张量磁梯度的斜导数原理 1.1 磁三分量数据转换全张量磁梯度数据的公式推导

求垂向导数就是根据磁场F的分布来计算在该平面上的分布,当n=1时,即计算垂向一阶导数.根据场的表达式和波谱表达式可以得出垂向一阶导数及其波谱为

(1)

式中为位场的波谱表达式,同样, 式中uv分别为xy方向的圆波数.

同理,若实测磁异常的波谱仍然为,则由傅里叶变换的微分定理可以得到:

(2)

(3)

式中为异常波谱计算x方向一阶导数的波谱,异常波谱计算y方向一阶导数的波谱,式中uv分别为xy方向的圆波数.

地面上任意点地磁场总强度矢量T(即磁感应总强度矢量)通常可用直角坐标来描述,设以观测点为其坐标原点,XYZ三个轴的正向分别指向地理北、东和垂直向下,则该点的T矢量在直角坐标系内三个轴上的投影分量分别为北向分量(Hax)、东向分量(Hay)和垂直分量(Za);这三个分量也就是本文重点开展讨论的航磁三分量.由于引力位在场源体外和场源体内分别满足拉普拉斯方程和泊松方程,因此T矢量在直角坐标系内三个轴上的投影分量分别为北向分量(Hax)、东向分量(Hay)和垂直分量(Za)同样在场源体外和场源体内分别满足拉普拉斯方程和泊松方程,设实测的北向分量的磁异常波谱为,实测的东向分量的磁异常波谱为,实测的垂直分量的磁异常波谱为.

北向分量的磁异常波谱xyz方向导数的波谱的表达式为

(4)

(5)

(6)

东向分量的磁异常波谱xyz方向导数的波谱的表达式为

(7)

(8)

(9)

垂直分量的磁异常波谱xyz方向导数的波谱的表达式为

(10)

(11)

(12)

再由反傅里叶变换后,可以得出北向分量Bxxyz方向的导数为BxxBxyBxz,东向分量Byxyz方向的导数为ByxByyByz,垂直分量Bzxyz方向的导数为BzxBzyBzz.

1.2 全张量磁梯度的斜导数原理

斜导数是基于导数理论的基础上发展的新方法.是位场在x方向的导数,是位场在y方向的导数,是位场在z方向导数,是位场的水平导数,Tdr为Tilt角度.

二维复杂结构体的解析信号为

(13)

其中是解析信号,F是位场的大小,是相位.

斜导数(Tilt Derivative)与相位θ是类似的,是用垂直梯度Vdr与水平梯度Thdr比值的绝对值的反正切值,即:

(14)

其中VdrThdr分别是位场F在垂直方向和水平方向上的一阶导数,位场在垂直方向上的一阶导数Vdr可能为正值也可能为负值, 在地质体区域上Vdr为正, 在地质体区域外为负, 位场在水平方向上的一阶导数总是正值;对于二维空间来说有:

(15)

xy是水平方向相互垂直的两个方向.在边界位置处斜导数值为零;Thdr是位场水平方向的一阶导数的绝对值,总为正值,在物体边界位置处取最大值;所以对于深部场源,即使它的垂直梯度和水平梯度都很小,两者的比值仍然会很大,地质体的边界基本在零位置处取得,不受埋深深浅的影响.

二维空间的全张量磁梯度均满足位场的概念,可以利用斜导数方法来提取地质体的边界,即:

(16)

(17)

(18)

Tdr_Bx为北向斜导数值,Tdr_By为东向斜导数值,Tdr_Bz为垂向斜导数值,同样BxzByzBzz可能为正值也可能为负值, 但是这和磁异常总强度磁异常不同的是,仅仅Bzz在地质体区域上为正, 在地质体区域外为负,Tdr_Bz值的零值位置基本对应地质体的边界,而BxzByz却有本质的不同,在地质体区域上既有正也有负,如果简单的对BxzByz取绝对值,对于利用斜导数来判断地质体的边界基本在零值位置处取得的目的是不一致的,经过深入研究,本文也针对全张量磁梯度数据提出了斜导数特征值的概念,为判断地质体的边界带来便利.

1.3 斜导数特征值

斜导数方法是对地质体边界反映较理想的一种方法,但是利用零值线来判断地质体的边界往往带来较大的误差,能否在正值区间范围内确定一个特征值,而仅仅绘制大于该特征值的等值线图,达到更精确地判断地质体的边界的目的呢?本文针对总该问题开展了研究.

斜导数方法的主要作用是突出地质体的边界,而地下地质体多以有限延伸板状体为主,因此本文设计了有限延伸板状体的理论模型并开展深入研究,图 1所示折线为有限延伸板状体中心剖面总场异常值斜导数Tdr的剖面曲线图,根据模型体边界所对应的斜导数数值来确定地质体边界的斜导数特征值.当Tdr的数值位于45°±5°之间时,对应于模型体的边界,经过多次理论模型研究得知,Tdr大约大于40°时,异常曲线基本上位于地质体上的位置,将Tdr大于40°的等值线绘出即可,40°转换为弧度后的值为0.7,即认为斜导数只要绘制弧度值大于0.7的等值线图就可以简单明了的判断出地质体的边界,而弧度值为0.7时本文就确定为斜导数的特征值.此时需要将北向和东向斜导数的计算公式进行适当地修改,公式为

图 1 有限延伸板状体总场磁异常斜导数Tdr剖面图 顶面埋深500 m,底面埋深900 m,模型体宽度为400 m,中心位置为(2000 m,-700 m),磁化强度为50 A·m-1,其磁化方向与地磁场方向始终保持一致. Fig. 1 Oblique derivative Tdr profile of total field magnetic anomaly of a limited extension plate The depth of the top surface is 500 m, the depth of the bottom is 900 m, the width of the model body is 400 m, the center position coordinate is (2000 m, -700 m), and the magnetization intensity is 50 A·m-1. The direction of magnetization is consistent with the direction of the geomagnetic field.

(19)

(20)

此时,再根据公式(18)、(19)、(20)就可以分别利用斜导数特征值方法来更加精确地判断地质体的边界了,因为对于垂直斜导数来说,正值仍然在地质体区域上,而且更加接近地质体的边界;对于北向和东向斜导数来说,所有正值既有在地质体区域上的,也有不在地质体区域上的,而特征值大大缩小了地质体的范围.

2 理论模型研究

在公式推导和理论分析的基础上,我们设计简单模型和复杂模型,分别赋予不同的地磁场参数与物性参数,来研究斜导数方法对于全张量磁梯度数据的边界识别效果.

2.1 简单模型

先从简单模型入手,在无磁背景中设计单一立方体的模型,该立方体的边长为400 m,中心位置埋深为400 m,磁化强度为10 A·m-1,模型几何特征如图 2所示.

图 2 简单模型XY平面及XZ平面剖面图 Fig. 2 XY and XZ plane profiles of simple model

假定磁化方向与地磁场方向一致,分别改变地磁场倾角和偏角形成不同的5种地磁参数与物性参数组合,如表 1所示.

表 1 不同磁化情况下简单模型参数表 Table 1 Parameter table of simple model under different magnetization conditions

分别计算简单模型在上述5种参数组合下的全张量磁梯度数据,并根据公式(18)、(19)、(20)得到全张量磁梯度的斜导数值Tdr_BxTdr_ByTdr_Bz.根据1.3节中斜导数特征值的分析,0.7为斜导数的特征值,揭示地质体的边界特征.绘制简单模型5种参数组合下Tdr (弧度值)大于0.7的等值线图,如图 34567所示.

图 3 垂直磁化1下的北向斜导数Tdr_Bx (a)、东向斜导数Tdr_By (b)、垂向斜导数Tdr_Bz (c)等值线图 Fig. 3 Contour maps of northward oblique derivative Tdr_Bx (a), eastward oblique derivative Tdr_Bx (b) and vertical oblique derivative Tdr_Bz (c) under vertical magnetization 1
图 4 垂直磁化2下北向斜导数Tdr_Bx (a)、东向斜导数Tdr_By (b)、垂向斜导数Tdr_Bz (c)等值线图 Fig. 4 Contour maps of northward oblique derivative Tdr_Bx (a), eastward oblique derivative Tdr_Bx (b) and vertical oblique derivative Tdr_Bz (c) under vertical magnetization 2
图 5 高纬度斜磁化下北向斜导数Tdr_Bx (a)、东向斜导数Tdr_By (b)、垂向斜导数Tdr_Bz (c)等值线图 Fig. 5 Contour maps of northward oblique derivative Tdr_Bx (a), eastward oblique derivative Tdr_Bx (b) and vertical oblique derivative Tdr_Bz (c) under high latitude oblique magnetization
图 6 低纬度斜磁化下北向斜导数Tdr_Bx (a)、东向斜导数Tdr_By (b)、垂向斜导数Tdr_Bz (c)等值线图 Fig. 6 Contour maps of northward oblique derivative Tdr_Bx (a), eastward oblique derivative Tdr_Bx (b) and vertical oblique derivative Tdr_Bz (c) under low latitude oblique magnetization
图 7 水平磁化下北向斜导数Tdr_Bx (a)、东向斜导数Tdr_By (b)、垂向斜导数Tdr_Bz (c)等值线图 Fig. 7 Contour maps of northward oblique derivative Tdr_Bx (a), eastward oblique derivative Tdr_Bx (b) and vertical oblique derivative Tdr_Bz (c) under horizontal magnetization

图 3为垂直磁化下的模型斜导数等值线图,由图可知,Tdr_Bx斜导数等值线图东西对称,在研究区域东西两侧出现与模型边界无关的等值线,而研究区域中心位置存在两个近似为椭圆的小圈闭,两圈闭的长轴与模型的东西边界相对应,但并不能准确圈定其东西边界. Tdr_Bx斜导数等值线没有实现对模型的南北边界的控制.与Tdr_Bx特征相似,Tdr_By斜导数等值线图南北对称,在研究区域南北两侧出现与模型边界无关的等值线,而研究区域中心位置也存在两个近似为椭圆的小圈闭,两圈闭的长轴与模型的南北边界相对应,同样不能准确圈定其南北边界. Tdr_Bz斜导数等值线圈定的范围与模型体的实际位置相吻合,边界识别效果较好,研究区域边界处没有出现假异常.图 4图 3存在着6°的偏角差异,但Tdr_BxTdr_ByTdr_Bz斜导数等值线图的特征与图 3一致.推测可能由于偏角较小,偏角的差异并未影响斜导数对于磁异常体边界的圈定.同时,模型这两种参数组合的斜导数等值线特征也表明0.7的斜导数特征的选取是适合的.

图 5为高纬度斜磁化下的模型斜导数等值线图,图 6为低纬度斜磁化下的模型斜导数等值线图.与垂直磁化相比,斜磁化下斜导数Tdr_BxTdr_ByTdr_Bz斜导数等值线图不再是轴对称,Tdr_BxTdr_ByTdr_Bz斜导数等值线均向南部偏移.而且,Tdr_BxTdr_By均在研究区域北部出现假异常,异常范围相较于垂直磁化变大,南部边界处假异常消失. Tdr_Bz受高纬度斜磁化的影响较小.图 7为水平磁化下的斜导数等值线图,Tdr_BxTdr_By斜导数等值线图基本可视为轴对称,对称轴分别为图框对角线和Y=2000 m,异常的中心点与模型的中心点位置相一致,但Tdr_BxTdr_By斜导数等值线图无法确定异常体的边界. Tdr_Bz斜导数等值线图在北部出现假异常,异常体位置向南偏移,中心偏移量约250 m,偏移量较斜磁化更大.

综上,根据单一立方体斜导数处理结果的分析,在地磁参数与物性参数一致时,磁化方向对于斜导数的边界识别效果影响较大,越接近垂直磁化,识别效果越好;垂直磁化条件下,Tdr_BxTdr_ByTdr_Bz斜导数等值线可以较为准确地识别磁性异常体边界,Tdr_Bx仅能识别南北边界,Tdr_By仅能识别东西边界,Tdr_Bz相对于其他两个方向的识别效果更好,可以较为准确地识别磁性异常体边界,且假异常较少;水平磁化下Tdr_BxTdr_By斜导数等值线图的中心可以指示异常体的中心点,但无法进行边界识别.

2.2 复杂模型

为了能够更加有效的对地下地质体的分布情况进行解释,本文设计了地磁场参数与物性参数一致情况下的复杂模型体,分析异常体复杂组合情况下,全张量磁梯度斜导数方法对于地质体边界的识别效果.从简单模型分析可知,垂向斜导数特征值Tdr_Bz对异常体边界识别效果较好,因此对于复杂模型的分析仅针对垂向斜导数特征值Tdr_Bz.复杂模型由4个位置不同的立方体组成,位置如图 8所示.

图 8 复杂理论模型XY平面及XZ平面剖面图 Fig. 8 XY and XZ plane profiles of complex theoretical model

模型1的坐标为(1250 m, 2850 m,200 m),其顶面埋藏深度为100 m,底面埋藏深度为300 m,其地质体规模为200 m×200 m×200 m,磁化强度为10 A·m-1.

模型2的质心坐标为(2000 m,2000 m,400 m),其顶面埋藏深度为200 m,底面埋藏深度为600 m,其地质体规模为400 m×400 m×400 m,磁化强度为10 A·m-1.

模型3的质心坐标为(3300 m,1600 m,700 m),其顶面埋藏深度为500 m,底面埋藏深度为900 m,其地质体规模为400 m×400 m×400 m,磁化强度为20 A·m-1.

模型4的质心坐标为(2000 m, 2000 m,1750 m),其顶面埋藏深度为1000 m,底面埋藏深度为2500 m,其地质体规模为1500 m×1500 m×1500 m,磁化强度为10 A·m-1.

为了检验垂向斜导数特征值边界识别的应用效果,应用张量不变量法、分量组合法两种现有张量边界识别方法与垂向斜导数方法进行模型对比试验.

张量不变量法的表达式为

(21)

分量组合法的表达式为

(22)

三种方法在不同磁化方向下的边界识别效果如图 91011所示.

图 9 低纬度斜磁化下复杂模型边界识别结果图 (a)分量组合法;(b)张量不变量法;(c) Tdr_Bz斜导数法. Fig. 9 Edge detection result diagram of complex model under low latitude oblique magnetization (a) Components combination method; (b) Tensor invariant method; (c) Tdr_Bz oblique derivative method.
图 10 高纬度斜磁化下复杂模型边界识别结果图 (a)分量组合法;(b)张量不变量法;(c) Tdr_Bz斜导数法. Fig. 10 Edge detection result diagram of complex model under high latitude oblique magnetization (a) Components combination method; (b) Tensor invariant method; (c) Tdr_Bz oblique derivative method.
图 11 垂直磁化下复杂模型边界识别结果图 (a)分量组合法;(b)张量不变量法;(c) Tdr_Bz斜导数法. Fig. 11 Edge detection result diagram of complex model under vertical oblique magnetization (a) Components combination method; (b) Tensor invariant method; (c) Tdr_Bz oblique derivative method.

低纬度斜磁化条件下,张量不变量法、分量组合法与Tdr_Bz斜导数法所识别出的边界均受到倾斜磁化的影响存在一定程度的变形,存在不同程度的偏移.张量不变量法受斜磁化影响向北偏移,但边界范围偏大,不能有效反映模型3的边界,未能反映模型4的存在.分量组合法受到一定程度斜磁化的影响,异常中心向南偏移,模型1和模型2的识别效果优于张量不变量法,未能反映模型3、4的存在;Tdr_Bz斜导数法受斜磁化影响最大,模型1、2、3向南偏移量较大,模型3有较为明确的反映,未能反映模型4的存在.低纬度斜磁化下,三种方法各有优劣,Tdr_Bz斜导数法优势不明显.

高纬度斜磁化条件下,张量不变量法、分量组合法与Tdr_Bz斜导数法所识别出的边界受倾斜磁化影响较小,整体偏移量较小,不影响整体边界范围的识别.张量不变量法与分量组合法的识别效果与低纬度斜磁化相似,识别的模型1和模型2边界范围均偏大,均不能有效反映模型3和模型4的边界,未能反映模型4的存在.分量组合法受到一定程度斜磁化的影响,异常中心向南偏移. Tdr_Bz斜导数法可以较为有效地识别模型1、模型2和模型3,其中模型1和模型2边界识别准确,与地质体实际边界吻合,模型3虽然由于深部异常体模型4的影响有假异常存在,但异常主体位置依然可识别.

表 2 低纬度斜磁化下复杂模型参数表 Table 2 Parameter table of complex model under low latitude oblique magnetization
表 3 高纬度斜磁化下复杂模型参数表 Table 3 Parameter table of complex model under high latitude oblique magnetization
表 4 垂直磁化下复杂模型参数表 Table 4 Parameter table of complex model under vertical magnetization

垂直磁化条件下,张量不变量法与分量组合法仅能有效识别模型1和模型2的存在,但对模型1、2的边界识别效果不好,边界范围偏大.而Tdr_Bz斜导数法可以较为有效地识别模型1、模型2和模型3,其中模型1和模型2边界识别准确,与地质体实际边界吻合,模型3虽然由于深部异常体模型4的影响有假异常存在,但异常主体位置依然可识别.

根据低纬度斜磁化、高纬度斜磁化和垂直磁化的模型试验结果可知,在高纬度斜磁化及垂直磁化下,Tdr_Bz斜导数法可以清晰的显示浅部地质体的边界,与地质体实际地质体边界相吻合,同时该方法有一定的深度识别能力,可以有效辨识不同深度的地质体,但深度越深边界识别效果越差.

3 实际数据应用

2017中国成功研制了具有国际领先水平的航磁三分量测量系统,获得了国内首个直接测量的航磁三分量数据.测区位于大兴安岭北段,该区域地磁倾角约为66°,磁偏角约为-11°,平均飞行高度200 m,数据测量精度优于30 nT,主要为晚古生代中-酸性侵入岩及中生代中酸性火山岩(佘宏全等,2012),该区域主要发育NE-NNE向断裂以及NW向、近EW向断裂.截取部分实测航磁三分量数据,实测三分量数据如图 12所示.实测航磁三分量数据转换得到全张量磁梯度数据,应用垂向斜导数方法对其进行处理,处理的结果如图 13所示.

图 12 测区实测航磁三分量异常图 (a)北向分量;(b)东向分量;(c)垂向分量;(d)矢量箭头图. Fig. 12 Aeromagnetic three-component anomaly map in the survey area (a) North component; (b) East somponent; (c) Vertical component; (d) Vector map.
图 13 (a) Tdr_Bz斜导数法边界识别结果图; (b)边界识别结果与垂向分量对比图 Fig. 13 (a) Tdr_Bz oblique derivative edge detection result diagram; (b) Edge detection result and vertical component comparison diagram

图 12可知,区域内主要的磁性异常走向为北东向、北西向、近南北向、近东西向,区域内高磁异常主要出现在东部测区边界.采用Tdr_Bz斜导数法边界识别方法得到的地质体水平分布如图 13a所示,磁性地质体主要为条带状分布,呈现北东、北西、近南北和近东西向,与航磁三分量实测数据的分析结果一致,且存在明显的嵌套特征,即条带状磁性体环绕多个不规则区域.将Tdr_Bz斜导数法边界识别结果图透明化处理后叠置在航磁三分量垂向分量图上,发现识别到的磁性异常体主要分布在垂向分量高磁异常和低磁异常的边界上,环形异常包裹主要异常圈闭.图 13b中东部边界高磁异常并没有识别到磁性体存在,推测垂向分量的高磁异常位于深部,超出了边界识别的深度,因而未能识别到异常体的存在.而多处地磁异常无磁性异常体分布可能表明在本方法的深度识别范围内没有磁性异常分布.

4 结论

本文基于常规的斜导数方法,推导了全张量磁梯度数据的斜导数张量识别方法,对斜导数特征值Tdr_BxTdr_ByTdr_Bz进行模型试验.模型试验及方法对比结果表明,Tdr_Bz斜导数具有较好的边界识别效果,可以清晰的显示浅部地质体的边界,同时该方法有一定的深度识别能力,可以有效辨识不同深度的地质体.但倾斜磁化会影响该方法的使用,在中低纬度应用时需要先进行化极处理.应用上述方法对大兴安岭地区的试验测区进行处理解释,取得了较好的应用效果.

致谢  感谢段树岭老师提供实际数据支持,感谢骆遥高级工程师和谢汝宽工程师对于文章写作的帮助,感谢审稿专家的宝贵建议!
References
Ferreira F J F, De Souza J, Bongiolo A D B E S, et al. 2013. Enhancement of the total horizontal gradient of magnetic anomalies using the tilt angle. Geophysics, 78(3): J33-J41. DOI:10.1190/geo2011-0441.1
Guo H, Wang P, Xie R K. 2012. The relationship of the boundary and the depthcomputed with the method of tilt-depth. Progress in Geophysics (in Chinese), 27(5): 2151-2159. DOI:10.6038/j.issn.1004-2903.2012.05.041
Guo H, Xiong S Q, Yu C C. 2009a. The application of tilt-depth technology. Geophysical and Geochemical Exploration (in Chinese), 33(5): 583-586.
Guo H, Yu C C, Wu Y G. 2009b. An improved tilt derivative method and its application. Geophysical and Geochemical Exploration (in Chinese), 33(2): 212-216.
Lahti I, Karinen T. 2010. Tilt derivative multiscale edges of magnetic data. The Leading Edge, 29(1): 24-29. DOI:10.1190/1.3284049
Miller H G, Singh V. 1994. Potential field tilt-a new concept for location of potential field sources. Journal of Applied Geophysics, 32(2-3): 213-217. DOI:10.1016/0926-9851(94)90022-1
Oruç B. 2010a. Depth estimation of simple causative sources from gravity gradient tensor invariants and vertical component. Pure and Applied Geophysics, 167(10): 1259-1272. DOI:10.1007/s00024-009-0021-4
Oruç B. 2010b. Location and depth estimation of point-dipole and line of dipoles using analytic signals of the magnetic gradient tensor and magnitude of vector components. Journal of Applied Geophysics, 70(1): 27-37. DOI:10.1016/j.jappgeo.2009.10.002
Oruç B, Keskinsezer A. 2008. Structural setting of the Northeastern Biga Peninsula (Turkey) from tilt derivatives of gravity gradient tensors and magnitude of horizontal gravity components. Pure and Applied Geophysics, 165(9-10): 1913-1927. DOI:10.1007/s00024-008-0407-8
Salem A S, Williams S, Fairhead J D, et al. 2007. Tilt-depth method:A simple depth estimation method using first-order magnetic derivatives. The Leading Edge, 2007(26): 1502-1505.
Schmidt P W, Clark D A. 2006. The magnetic gradient tensor:Its properties and uses in source characterization. The Leading Edge, 25(1): 75-78. DOI:10.1190/1.2164759
Schmidt P W, Clark DA. 2000. Advantages of measuring the magnetic gradient tensor. Preview, 85: 26-30.
She H Q, Li J W, Xiang A P, et al. 2012. U-Pb ages of the Zircons from primary rocks in middle-northern Daxinganling and its implications to geotectonic evolution. Acta Petrologica Sinica (in Chinese), 28(2): 571-594.
Verduzco B, Fairhead J D, Green C M, et al. 2004. New insights into magnetic derivatives for structural mapping. The Leading Edge, 23(2): 116-119. DOI:10.1190/1.1651454
Wang L F, Guo C C, Xue D J, et al. 2016. Analytic signals of magnetic gradient tensor and their application to estimate source location. Progress in Geophysics (in Chinese), 31(3): 1164-1172. DOI:10.6038/pg20160333
Wang X, Li T L. 2004. Locating the boundaries of magnetic or gravity sources with Tdr and Tdr-Thdr method. Progress in Geophysics (in Chinese), 19(3): 625-630.
郭华, 王平, 谢汝宽. 2012. 利用Tilt-depth求取埋藏深度及其与边界位置的关系研究. 地球物理学进展, 27(5): 2151-2159. DOI:10.6038/j.issn.1004-2903.2012.05.041
郭华, 熊盛青, 于长春. 2009a. Tilt-depth方法的应用. 物探与化探, 33(5): 583-586.
郭华, 于长春, 吴燕冈. 2009b. 改进的斜导数方法及应用. 物探与化探, 33(2): 212-216.
佘宏全, 李进文, 向安平, 等. 2012. 大兴安岭中北段原岩锆石U-Pb测年及其与区域构造演化关系. 岩石学报, 28(2): 571-594.
王林飞, 郭灿灿, 薛典军, 等. 2016. 磁梯度张量解析信号分析法及其在场源位置识别中的应用. 地球物理学进展, 31(3): 1164-1172. DOI:10.6038/pg20160333
王想, 李桐林. 2004. Tilt梯度及其水平导数提取重磁源边界位置. 地球物理学进展, 19(3): 625-630. DOI:10.3969/j.issn.1004-2903.2004.03.022