-1-航磁数据位场转换处理及效果航磁T测量数据是不同深度、不同形态、规模的磁性地质体磁场信息在观测面上的综合反映。由于场的叠加效应,使得某些具有一定地质意义的异常变得复杂,在原始图件上很难识别,给地质解释工作带来了难度。为了提高对航磁异常的分辨能力,突出更多有用信息,根据测区航磁异常特征和地质解释需要,对原始测量数据进行了原平面化极、上延、垂向一阶导数以及剩余异常提取等几种位场转换处理。第一节位场转换处理及效果航磁平面网格数据位场转换处理采用表达式简单、运算速度快捷的频率域算法,进行化极、导数换算、解析延拓等处理。频率域转换的过程是:首先对异常资料进行傅立叶正变换,以得到异常资料的频谱;而后把异常的频谱和与转换相应的频率相应函数点积,得到处理后异常的频谱;最后对处理后异常的频谱进行傅立叶反变换,从而得到处理后的异常。位场转换处理使用的软件是中国国土资源航空物探遥感中心自主开发的WINDOWS系统下地球物理数据处理解释软件(GeoProbeMager)及航空物探彩色矢量成图系统(AgsMGis)。一、原平面化极处理化极,即化磁极,就是把斜磁化异常转变为垂直磁化异常,相当于在磁北极观测异常。测区处于中纬度地区,由于倾斜磁化的影响,造成磁异常中心不是正好对应在地质体的正上方,而是相对于地质体的中心向南部产生一定的偏移。这对于确定磁性地质体的空间位置、形态、分布范围以及对磁异常的定性定量解释均带来一定的困难。化极可用于消除由于非垂直磁化引起的异常不对称性,在剩磁很小或感磁远大于剩磁且两者方向一致的情况下,将实测的斜磁化异常转化为垂直磁化异常,这样可以较为准确的确定异常的场源位置,提高异常解释的定位精度。从而使异常形态简化,并与磁性体位置保持一致,有利于圈定磁性体边界和走向。作化极处理时要注意剩磁的影响,化极处理一般都假定磁化方向与地磁场方向一致,对于那些剩磁远远大于感磁且剩磁方向与地磁场方向不一致的磁性体就不符合这一假设条件,特别是测区中的火山岩分布区,由于剩磁较大会出现磁场畸变现象,使用时应注意甄别。从项目组野外物性测量结果看,区内多数岩石以感磁为主,剩磁方向与感磁方向接近,符合化极的前提条件。全区采用"频率域偶层位变倾角磁方向转换方法"实现磁场全变倾角化极。在观测面上建立笛卡尔直角坐标系,使x轴志向磁北,z轴垂直向下。假设观测场T是-2-一分布在观测面下方z=h平面上的偶层磁荷面引起的。它在观测点P(x,y,z)处产生的磁位U与磁场T分别为114UMddr(1)0TutU(2)222MMlxyzh(3)设,,xyz00ttt(4)00lll(5)式中M表示偶层磁荷面的磁化强度矢量;0t、0l表示均为常矢量研究区t、l的平均值;0t、0l表示他们的变化值。将(4)式代入(2)式并进行傅氏变换得:000FTugtFUFtU(6)移项得00011FUFTFtUgtu(7)式中222,2,2,giuivffuv对(1)式两端作傅氏变换得:21/4fhzFUgFMef(8)同样的,将(5)式代入(3)式进行傅氏变换后代入(8)式得:20014fhzFUglFMgFMlef移项得20014fhzFMfeFUgFMlgl(9)(7)式、(9)式即为频率域变倾角磁方向转换的两个基本公式。已知观测场T,可应用(7)式计算FU,再将FUFU代入(9)式中计算FM。再把FM代入地磁极处的FT极,即可实现变倾角化极处理。FT极与FM的关系如下。在地磁北极有:00000,0,1,0tltl,由(7)式、(9)式可得:20fhzFMufeFM极(10)-3-经傅立叶反变换后可得化磁极长Z210fhzZFufeFM(11)假设要换算的场为pT,其磁场方向单位矢为t,磁化方向单位矢为l,则只要把他们代入(7),令0t、0l偏差为零,即得FU的初值0FU:00011pFUFTugt对0FU反变换求得初值0U,把0U代入(7)式得1FU,如此反复迭代,直到求得的U值之差小于给定的标准为止。求得FU后,类似地取0FU为2004fhzfFMFUegl反变换求得0M,把0M代入(9)式,求得1FM,如此反复可最后求得FM,并代入(11)式求得化磁极磁场。由于本区处于xxx°xxx′~xxx°xxx′,属于中纬度地区,斜磁化能够产生一定影响,对原磁场数据进行化极处理后,在垂直磁化的条件下,磁异常的形态以及磁异常与磁性体的关系都比较简单,便于进行地质解释。对比航磁T等值线平面图和航磁T化极等值线平面图,航磁化极处理作用非常明显(图1...