航磁数据处理资料

航磁T测量数据是不同深度、不同形态、规模的磁性地质体磁场信息在观测面上的综合反映。由于场的叠加效应，使得某些具有一定地质意义的异常变得复杂，在原始图件上很难识别，给地质解释工作带来了难度。为了提高对航磁异常的分辨能力，突出更多有用信息，根据测区航磁异常特征和地质解释需要，对原始测量数据进行了原平面化极、上延、垂向一阶导数以及剩余异常提取等几种位场转换处理。第一节位场转换处理及效果航磁平面网格数据位场转换处理采用表达式简单、运算速度快捷的频率域算法，进行化极、导数换算、解析延拓等处理。频率域转换的过程是：首先对异常资料进行傅立叶正变换，以得到异常资料的频谱；而后把异常的频谱和与转换相应的频率相应函数点积，得到处理后异常的频谱；最后对处理后异常的频谱进行傅立叶反变换，从而得到处理后的异常。位场转换处理使用的软件是中国国土资源航空物探遥感中心自主开发的WINDOWS系统下地球物理数据处理解释软件(GeoProbeMage)及航空物探彩色矢量成图系统化极，即化磁极，就是把斜磁化异常转变为垂直磁化异常，相当于在磁北极观测异常。测区处于中纬度地区，由于倾斜磁化的影响，造成磁异常中心不是正好对应在地质体的正上方，而是相对于地质体的中心向南部产生一定的偏移。这对于确定磁性地质体的空间位置、形态、分布范围以及对磁异常的定性定量解释均带来一定的困难。化极可用于消除由于非垂直磁化引起的异常不对称性，在剩磁很小或感磁远大于剩磁且两者方向一致的情况下，将实测的斜磁化异常转化为垂直磁化异常，这样可以较为准确的确定异常的场源位置，提高异常解释的定位精度。从而使异常形态简化，并与磁性体位置保持一致，有利于圈定磁性体边对于那些剩磁远远大于感磁且剩磁方向与地磁场方向不一致的磁性体就不符合这一假设条件，特别是测区中的火山岩分布区，由于剩磁较大会出现磁场畸变现象，使用时应注意甄别。从项目组野外物性测量结果看，区内多数岩石以感磁为主，剩磁方向与感磁方向接近，符合化极的前提条件。全区采用"频率域偶层位变倾角磁方向转换方法"实现磁场全变倾角化极。在观一分布在观测面下方z=h平面上的偶层磁荷面引起的。它在观测点P(x,y,z处产生的磁位U与磁场T分别为T1MgddrX'yoo式中M表示偶层磁荷面的磁化强度矢量;均o表示他们的变化值。将(4)式代入(2)式并进行傅氏变换得：移项得uo对(1)式两端作傅氏变换得：to、lo表示均为常矢量研究区t、I的平⑺°°同样的，将(5)式代入⑶式进行傅氏变换后代入(8)式得:oo移项得1ogo(7)式、(9)式即为频率域变倾角磁方向转换的两个基本公式。已知观测场极极o经傅立叶反变换后可得化磁极长Zofe2fhzFM假设要换算的场为T，其磁场方向单位矢为t,磁化方向单位矢为I，则只要把p他们代入⑺,令to、Io偏差为零，即得FU的初值FoU:o差小于给定的标准为止。求得FU后，类似地取FU为ooo反变换求得M，把M代入(9)式，求得FM，如此反复可最后求得FM，并代入(11)oo磁场数据进行化极处理后，在垂直磁化的条件下，磁异常的形态以及磁异常与磁性体的关系都比较简单，便于进行地质解释。对比航磁T等值线平面图和航磁T化极等值线平面图，航磁化极处理作用非常明显(图1):局部异常整体向北偏移，表明通过化极处理，使异常回归到磁性地质体上方；减小或消除了由于斜磁化而引起的多数局部异常正负异常伴生现象，为进一步圈定岩体边界创造条件；使异常带及梯度带更加明显，有利于揭示出不同地质体的分布与形态，对圈定各种不同类型的断裂、确定磁性体的性质及边界航磁局部异常通常是叠加在区域背景场上的次级异等基础图件中表现并不明显，需要通过一定的数学处理手段来突出其特征。垂向导数处理是解决这个问题的一种有效手段，它反映了磁场在垂直方向上的梯度变化，在增强由浅部磁性体引起的局部异常、压制长波区域场有很强的功能，可以突出在总场图上不明显的细节，并能分解横向叠加异常，理论上导数的次数越高，这种分辨能力就越强。磁异常垂向导数换算公式如下:女口果令SzxX,y,z、Szyx,y,z、爲x,y,z及2zx,y,z、Szyyx,y,z、Szzzx,y,z分别为Zax,y,z对X、y、z的一阶导数及二阶导数的频谱，则有微分定力易于得到:zxzy22e2ee2e2u2v22u2v2221/2gzu2v221/22u2v21/2gz同理，可以写出:SS2222由此可知，求磁场的n阶垂向导数的频谱，应乘上的导数因子为而求磁场沿x方向或y方向的n阶水平导数的频谱，应乘上导数因子为频谱（即求x,yxnyzmx,yxnyzm航磁T垂向一阶导数已经广泛地应用于磁异常的解释，它能区分相邻磁性体异常，减少其相互叠加的影响，并把叠加在背景场中的局部异常分离出来，是压制区域场，圈定局部异常，分离叠加异常的常用方法。在实际磁场转换处理中，由于垂向一次导数相波器，在突出高频异常的同时，也突出了测量、磁场调平等干扰误差。对本区化极场的数据后的图件与原磁场图相比（图2a、b）,突出了浅部磁性体信息，而压制了深层区域场的影响。该处理也消除或减弱了局部异常之间的叠加和干扰现象。因此，航磁T化极垂向一阶导数处理在提取强背景场中的弱缓异常，圈定局部异常、火山构造、划分构造边界等方面具有根据厚板状磁性体异常公式，垂向二阶导数的零值线为磁性体边界位置。因此，航磁T化极垂向二阶导数处理的主要目的是利用航磁异常垂向二次变换率来圈定磁性体的范围和边界。本区航磁T垂向二阶导数处理是在化极处理的基础上，对化极后的网格数据采用频率域位场转换方法求取磁异常沿垂直方向上的二次变换率，并编制了航磁T垂向二阶导数等值线平面图（图2c）。在理论上，经垂向二阶导数处理后，区域场得到了进一步的压制，很大程度上消除了深部磁性体的影响，使得磁性体的范围和边界更加明显，仅供参考使用。三、化极0°方向水平一阶导数处理化极0°方向水平一阶导数处理的目的是突出异常在东西向的线性特征，分辨东西方向上构造线的展布，以准确的划定浅层构造、断裂构造，以便推断区内的构造格架。磁异常水平导数换算公式如下：女口果令SzxX,y,z、Szyx,y,z、Sx,y,z及Szxxx,y,z、Szyyx,y,z、Szzzx,y,z分别为Zzx,y,z磁T化极垂向二阶导数等值线平面图；对X、y、z的一阶导数及二阶导数的频谱，则有微分定力易于得到：SzxSzy22e2222uv2e1/2gzgSzz同理，可以写出:2u2v2u1/2gzSS22222222nn而求磁场沿x方向或y方向的n阶水平导数的频谱，应乘上导数因子为nI是实测平面上任一方向，它与X轴的夹角为a，则有:Ixy两边作傅氏变换并应用微分定理，得知利用(13)式即可实现磁场的频率域方向导数计算，当a0，代入(13)式即可求得0航磁T化极0°方向水平导数处理结果显示(图3)，局部域近东西向的线性异常特征及弧形异常特征都非常明显，为该区划分浅层构造、近东西向断裂构造等提供依据。力ee磁场向上延拓处理就是将原观测面上的磁场值向上换算到另一个高度面上。随着上延高度的增加，磁性体引起的异常幅度按指数规律衰减。衰减最快的为浅部局部磁性体引起的高频异常成分，而具有一定延伸的大规模磁性体引起的低频异常成分衰减较慢。可见，向上延拓处理起到压制浅部小规模磁性体异常而突出深大地质体异常的作用。设场源位于z=H平面一下(H>0)，则磁场在z=H平面以上是对x、y、z的连续函数。若z=0观测平面上的磁场Tx,y,0为已知，可以得到向上延拓公式为y/32由褶积积分公式可知，上式为Tx,y,0与1zy2z23/2关于变量x,y3/2二维褶积。空间域的褶积与频率域的乘积相对应。下面分别求Tx,y,0及223/2傅立叶变换，设Tx,y,z对于变量x,y的傅立叶变换为Spu,v,z，有则利用上式可以由已知的Tx,y,0求出其频谱Sru,v,0。进一步求z的2x2y2z23/2傅立叶变换，应用Erdelyi(1954)给出的积分变换表可以得到：z2u2v21/2gzu2v21/2gzr(19)式即为向上延拓的频谱表达式。通过磁场向上换算，相当于加大了观测面与场源的距离，可以使局部小规模异常随换算高度的增加而减小，而深部规模较大的磁性体所产生的异常更加凸出。为了了解深源磁性体的特征及航磁异常随高度衰减变化特征，判断磁性体的埋深及延伸情况，在化极基础上进行了0.5km、1.0km、3.0km、5.0km四种不同高度的向上延拓处理。通过不同高度的向上延拓，消除了高频磁异常的干扰，使得磁场面貌逐渐单调，达到了突出低频区域异常的目的，对了解深源磁性体的特征和基底构造具有一定的对比测区上延高度磁场图可以看出(图4)：化极上延0.5km后，高频干扰异常被压制，有意义的局部异常基本保留；化极上延1.0km后，规模较小的局部异常衰减得很快，中等规模的异常明显突出，区域磁场面貌反映得更加清晰；化极上延3km后，由于测区覆盖较浅，而引起局部异常的磁性地质体延伸有限，高频异常几乎消失；化极上延5km后，有效的压制了浅部磁性体引起的异常，突出了深源低频磁异常。因此，化极上延0.5km或1.0km磁场图，对研究本区区域构造、划分隐伏岩体非常有效；上延3.0km后仍然存在的磁异常则反映出了规模较大、延伸较深或埋深较大的磁性地质体；上延5.0km后反映的深源低频磁异常，对于确定磁性基岩、深大断裂及区域构造格架有着重要意义。五、化极匹配滤波求取局部及区域场区域场和局部场的分离问题是航磁数据处理的主要内容之一，对实际资料的解释有重要意义。利用匹配滤波算法可对航磁数据进行区域场和局部场的分离，进一步达到突出浅部异常或突出深部异常的目的。知，这类磁性体可以用随机分布的偶极源组成的等效层代替。设偶极源等效层的深度为d2，源的偶极磁距为mdp,,d2，其傅氏正变换为mdpu,v,d2；并设观测面高度为z，则以求出偶极等效层所产生的场的振幅谱为(为简洁式中略去u0/4，并不影响结果)Adpu,v,z42.u2v2Mdpu,v,d2e2d2zuv显然，振幅Adpu,v,d2是和•.u2—v2Mdpu,v,d2e2八丹成正比的。在假定区域场是由许多下延伸很大的磁性体所引起，它们可以用由随机分布的电极等效层所代替。设点极等效层深度为di,电极源的磁荷密度为mp,,di，其傅氏正变换为叫u,v,di，则由点极等效层引起场的振幅谱为当不考虑干扰场时，设实测场为区域场与局部场之和，它的振幅谱为Au,v,z2Mpu,v,die2dl42Ju2v2gMdpu’vQe2d2z^^为了从实测场中提取区域场，可以设计一个只对区域场响应的滤波器，即与区域场匹配的滤波器。令滤波器的频率响应函数为)将(21)、(22)式代入(23)式可得:u,v,z，为使它仅对区域场响应,p1e2MdpVuv2TUTVd2dip实际上此时u,v,z与z无关，应为u,v。求出频率相应u,v后，由(23)式求出Au,v,z，用无相移的傅氏反变换(即相位谱不变)，就获得了区域场。同样，如果从实测场中提取局部场，也可以设计个只对局部异常匹配的滤波器，用同样的方法也可以得到它的频率响应函数为2Su^V2d2d式中，X表示点距，k表示直线斜率，h为区域场与局部场分离的最佳深度。航磁数据处理中，匹配滤波法是分离局部与区域异常场的重要手段之一，因为匹配滤波器是已知相关滤波器，它要求二者具有明显差异的波数成分，该工区的航磁异常特征满足这个条件（图5浅部大多为带状、块状或零散分布的火山岩、侵入岩体，提取区域场时，它是一个低通滤波器，这与数学解析向上延拓不同之处在于它有一个较为复杂的类似于汉宁滤波器的窗函数。在提取高频成分时，它不会放大导致高频成分的振荡效应，因为高通时的滤波器渐近线为1。与向上延拓相比,该方法简单易行，而且还能够获取分析局部场与区域场的相关窗口。通过匹配滤波算法可以分离浅部异常信息，通过不同的匹配因子，可以逐层剥离异常，为地下异常区域的分层解析提供了技术手段。六、剩余异常剖面平面图对于那些频率高、幅值低的航磁异常，在航磁T等值线平面图上受网格化取数圆滑滤波影响以及成图精度的制约，往往显示不清或被漏掉；而在航磁T平面剖面图上往往叠加在较大异常的背景场上，不容易识别出来。为了从剖面上将一些局部异常和较微弱的短波异常从区域背景中分离提取出来，采用空间域非线性滤波法对测区剖面数据进行非线性滤波处理计算，提取航磁剩余异常。非线性滤波方法是根据剖面异常曲线的曲率变化，按一定的异常宽度窗口拟合计算出剖面的区域背景异常值，然后将其从原始剖面异常中减去，所获得的剖面异常即为相应的剩余异常。该方法可以有效的提取指定波长的剩余磁异常，具体计算公式如下和迭代步骤如下(图6)：步骤i用剖面实测异常作为第一次圆滑的原始数据，用(26)、(27)、(28)式估计给定滤波窗口宽度WD两倍内异常曲线弯曲方向的特征值。SiTFKmFTFKS<3TF0.5gTFK2步骤2如果Ski、Sk2、Sk3同时满足下面条件:Q点距(30)、(31)式计算出Kmi、K、Kpi三点区域异常值，然后向前移动窗口，回到步骤1。TFKmTFKmiSki.0.5TpiSk3.0.5T步骤3如果Ski、Sk2、Sk3不同时买足步骤2的条件，则向前移动窗口，步长等于滤波窗口宽度，回到步骤i,指导整条剖面计算完成为止。步骤4以点距为单位，逐步减少窗口宽度，用新得到剖面异常代替i2、3,共重复WD次，使滤波窗口内每一个异常点都得到处理，提高滤波效果。步骤5为保证计算精度，重复上述步骤i2、3、4过程2~5次，获得剖面区域异常。为了获得较好的滤波效果，根据本区地质解释的需要，并分析对比不同滤波异常窗口处理结果信息，经过使用对比，本区最后选用滤波异常半宽度羽000m,沿剖面逐线完成对航磁弱缓异常的提取。从航磁弱缓异常处理效果看，消除了背景场的影响，分离了叠加异常，局部异常的形态更加完整清晰，并对弱小异常起到增强作用。如图7中，经过弱缓异常处理后，弱小磁异常在数条测线上反映明显，异常边界也更加清晰。目前对起伏地形条件下航磁数据进行定量正反演计算主要采用最优化算法，原起伏测量面上直接进行解释。本次研究以最优化理论为基础，使用二维半模型在为了提供模型正演计算的速度