二维数字滤波问题讨论——介绍一个垂直二阶导数公式

物探重磁数据处理工作中的垂直二阶导数法,由于具有较高的分辨力以及对区域场的压制作用,使它在实际工作中得到了越来越广泛的应用。前人曾通过不同途径求出了各种垂直二阶导数近似计算公式。本文根据二维数字滤波原理,通过计算垂直二阶导数滤波运算频率响应的反富立叶变换,求得一组权系数,由此归纳出一个新的近似计算公式。最后以实例检验本公式的适用性。     

富里叶变换解释倾斜厚板状体磁异常的应用

迄今,国内外对磁异常进行富里叶变换求其频谱以解释磁测资料方面的文章发表得很多,V.L.S.伯赫美树卡雷姆等人的《用富里叶变换解释岩脉磁异常》一文即是其中之一。该论文在频率域给出了岩脉顶宽、顶部埋深及磁化方向与岩脉之间夹角等参数的准确解法,但尚存在着如下缺陷:     

重磁异常在频率域的解释方法(下) 四、振幅谱的解释

四、振幅谱的解释上一期已经提到,重磁异常的频谱可表示为某些因子之乘积(85)。而且对于振幅谱,有关因子还要少,因为某些因子仅与相位谱有关。对于二维异常和三维场在频率轴上的振幅谱,有关的因子为四个,即     

重磁异常在频率域的解释方法(上)

频率域解释方法,指的是利用异常的频谱特性进行反演解释的方法。异常的频谱分为振幅谱和相位谱两部分,它们是作为空间座标函数的位场的付立叶变换。在国外,异常频率域解释的研究兴起于六十年代,虽然经过了一段相对冷落的时期,但近年来又重新活跃起来。在这篇文章中,我们拟讨论以下内容:(1)频谱的计算,(2)位场频谱的基本规律,(3)规则物体异常的频谱,(4)振幅谱的解释方法,(5)相位谱的解释方法。最后我们通过一些实例来看看频率域解释方法的特点。     

用磁异常的谱及其导数求场源埋深和磁矩

随着电算的普及,物探重磁数据的处理已广泛使用了数字滤波技术,旨在提取信息、压低干扰的一类线性滤波,例如:根据信号实际输出与希望输出的方差最小为原则设计的维纳滤波和卡曼滤波;根据输出信号的瞬时功率与干扰的平均功率之比最大为原则设计的最大讯噪比滤波(D.O.诺斯,当干扰为白噪声时称为匹配滤波);根据:①在信号失真最小的同时将干扰(包括随机干扰和规则干扰)压低到一定程度;②在信号保持一定程度失真的同时,将干扰压低到最小为原则设计的最佳线性滤波(B.H.斯特拉霍夫-空间域,王继伦-频率域);还有根据反馈原理,采用Z-变换设计的旨在节省计算时间和减少边缘损失的递归滤波(B.K.布哈塔恰赖阿),等等.上述滤波只有事先知道输入信号和干扰的谱的统计表达式时,才能设计出一个可供实用的数字滤波器.但在实际问题中,这些谱的准确表达式当然是不知道的.常用的处理办法是将场源近似看作是形状规则的物体,只考虑场源埋深h和     

用磁异常二维振幅谱确定磁性体埋深

目前电子计算机已经广泛应用于重磁异常的数据处理,借助于电子计算机利用磁异常的二维振幅谱来确定磁性体的埋深,是一个可行的方法。现对该方法作一简单介绍。一个周期函数,只要满足一些并不很严格的条件,就可以展成富里叶级数:     

重磁异常的方向维纳串联滤波

本文从频谱分析和最佳线性滤波理论出发,讨论了重磁方向滤波和维纳滤波的应用范围和局限性,提出了重磁异常的方向维纳串联滤波方法,从理论上进行了论证,通过垂直和水平迭加模型试算,对三种滤波结果作了对比和误差分析,讨论了异常交角、埋深、点距、地质体的形态和空间位置对滤波效果的影响,验证了方法的可靠性和适用性。     

重磁异常—维频率域位场变换与滤波

重磁异常一维频率域位场变换与滤波是二度体重磁异常推断解释的基础工作之一。编制在PC—1500机上实现的程序系统,大大有利于剖面重磁数据的现场处理。本文提供的程序系统具有功能齐全,使用灵活,处理数据量大的优点。     

第三讲 频率域中的位场变换

为了突出有用异常、压制各种干扰,在重磁异常解释工作中经常要进行延拓、求导、化向磁极、磁场各分量之间的转换以及分离区域场和局部场等多种处理。通常把上述处理统称为位场变换。在第一讲中我们已经提到,从频率域看,位场变换就是某种滤波。这一讲的中心内容就是具     

三度重磁异常在波数域的变换及滤波

对重磁异常进行变换及滤波,可以在空间域实现,也可以在波数域实现。用专门设计的算子(权系数组)对位场数据作卷积,是在空间域实现这类变换或滤波的一般形式。由富氏变换的卷积定理可知,二函数卷积的富氏变换等于两函数各自的富氏变换的乘积。在波数域实现上述位场     

Tau-p变换的频率域算法及Tau-p域的偏移处理

针对时间—空间域倾斜迭加算法的一些缺点,本支根据倾斜迭加运算与富里叶变换的关系,提出在频率一波数域实现τ—p变换。在频率一波数域进行τ—p变换,还可与视速度滤波相结合,以便更有效地消除干扰波,提高处理质量。根据τ—p变换频率一波数域算法原理,提出了在τ—p域直接实现偏移处理,同时指出,为了避免假同相轴的干扰,必须进行插值处理。     

解斜磁化水平圆柱体磁异常的JU变换法

目前，求解各种形体磁异常的富里叶变换法，特别是快速富里叶变换法，已被广泛应用；希尔伯特变换方法也已引入了磁异常的解释领域。笔者在学习特殊函数论时发现，建立一种正交于1     

对二度重磁资料数字滤波的一些讨论

数字滤波是抑制干扰、提高信噪比的一种手段,在物化探数据处理工作中获得了日益广泛的应用.在重磁资料处理和推断解释工作中,对异常进行光滑、求导数或求剩余异常,目的也是压低干扰突出有用信号.解析延拓是异常解释方法的一种,它实际上是对异常的频谱进行某种特定的改造.对于二度体,数字滤波的运算可用下面的褶积分方程来描述:     

重磁异常的直接推断解释方法——(在异常体为水平柱状体的情况下)

前言在任何物探异常定量解释的实践中,通常假定干扰体具有规则的形状或是已知干扰平面的深度.偏离了这种假设,将导致解释结果含糊不清.到目前为止所研究成功的各种简单几何形体的重磁异常的解释方法,可大致划分为下列四类:1.经验法;2.特征曲线法;3.曲线拟合     

综合物探在四川会理拉拉M163-1航磁异常查证中的应用

为了调查拉拉M163-1航磁异常的精细磁异常特征与深部电性结构特征,采用地面高精度磁法与音频大地电磁测深法,并结合区域重磁资料与地质资料对异常区开展了综合物探查证工作。通过区域重磁、地面高精度磁法与音频大地电磁测深资料的综合分析,得出了与前人不同的认识,并初步确定了M163-1航磁异常的精细磁异常特征、异常体的电性结构特征。这为该异常的找矿远景调查评价提供了重要的物探依据。     

频谱分析的一般知识(第一讲)

由于物探工作中广泛地应用了电子计算机,因此大大提高了数据处理与解释工作的速度与精度。目前,由于航空物探的发展,获得了大量的数据,对于这样大量的数据若只用以往空间域的处理方法与解释方法,虽然用高速电子计算机进行处理仍感不便,这就提出了怎样提高方法本身的速度。而频谱分析方法在作多种数据处理中比空间域方法快得多。那么如何进行频谱分析、空间域资料转换成频率域资料,频率域资料的解释等问题就是本讲座的主要内容。本讲座共分四讲:一、频谱分析的一般知识,二、快速富氏变换,三、频谱分析在位场变换中之应用,四、频谱分析在正反问题中之应用。     

介绍一种重磁异常的快速正演方法

前言最近文献中刊载了一种快速计算任意形状物体重磁异常的富氏变换法,其基本原理如下: 具有均匀密度或均匀磁化强度分布的任意形状物体的重力与磁力异常均可表示为场源几何形状函数与格林函数的褶积。因此,它们的频谱可以表示为这两个函数的富氏变换的乘积。要计算物体引起的重磁异常,关键在于给出格林函数与场源几何形状的富氏变换,然后再对这两个函数的富氏变换的乘积作反富氏变换。场源几何形状由描述物体形状的解析式给出,格林函数则仅与观测点的位置和磁化强度矢量或密度有关,而与场源的几何形状无关。所以,对于某种场位来说,只要一次     

在第四系覆盖区如何辨别基性岩体与铁矿引起的重磁异常

铁矿勘查工作采用重、磁手段一般非常有效,但在第四系覆盖区有基性岩体存在情况下,就会产生较强的重磁干扰,往往造成误判,因此,区分基性岩体与铁矿引起的重磁异常成为勘查工作的关键。以滦南杜蒿坨铁矿勘查为例,通过对磁法资料和重力资料的处理、解释,分析了基性岩体重磁异常特征及对铁矿异常的掩盖作用,总结了辨别铁矿异常的方法。并经钻探验证这些方法可靠和适用,在今后铁矿勘查工作中有一定的指导意义。     

重磁数据的线性滤波与反演失真问题

文中对国内研制的4种主要线性滤波方法进行了对比分析。选出一种具代表性的滤波方法,用计算机模拟观测场,探讨了在复杂干扰下重磁异常因滤波引起的失真,并就重磁异常在因滤波而引起失真的情况下如何反演解释,提出了可行的办法和建议。     

综合物探在会东航磁异常查证中的应用

根据会东航磁异常的基本特征,判断其异常为深部多个磁性体的综合反映.采用地面高精度磁测和音频大地电磁测深法,并结合区域重力勘探资料对异常区开展了综合物探调查.通过对区域重力、高精度磁法资料和音频大地电磁测深数据开展综合解释,初步确定了异常体在深部的空间展布特征,为在该区实施进一步的勘查工作提供了重要找矿线索.