用等效磁源法进行磁异常转换

用等效磁源法实现常规磁异常换算(延拓、求导、假重力异常换算、化磁极)的模型计算算例显示了等效源法在使用上的灵活性和方便性。与频率域换算方法相比,说明了等效源法换算结果更优于频率域换算结果。     

磁源重力技术及其应用

磁异常的复杂性对资料处理和解释带来了一定的困难,而重力异常相对更为简单、更易处理和解释。通过运用空间域的泊松公式,可以在频率域内推导出由磁异常换算重力异常的计算公式,这种由磁异常换算所得到的重力异常称为磁源重力异常。这里对理论模型的磁异常数据进行了计算与实验分析,并对实际航测的磁异常资料进行了计算与处理,取得了良好的应用效果,表明了磁源重力技术的可行性、正确性和精确性。     

起伏观测面上位场导数换算的B样条法

曲面位场导数的换算方法大多基于调和函数的积分方程数值解。本文提出了基于B样条函数的曲面上位场导数直接计算法。引入按弧长S作变量的三次样条插值,得出了位场水平导数的计算式。该方法还有一个特点,那就是可以用一元样条来实现曲面上的求导。这些,都决定了该方法的简便性和精确性。利用曲线上和曲面上的两种理论模型进行试算,分别计算了Zα/x和Δg/Z~2计算和理论的结果对比表明,符合很好。误差主要取决于曲面起伏大小,异常变化陡度和取样的间隔。     

用等效源法进行重磁异常换算的效果

在磁力、重力异常的解释中,常常要进行异常场的换算,例如向上延拓、向下延拓,曲化平(将起伏地形上的观测值换算到某一水平面上的场值),磁异常不同分量之间的换算,求变换磁化方向后的异常(包括计算异常的化到磁极),以及异常的一次、二次垂向导数等等。这些换算已有一些方法可以进行,在上延、不同分量的换算和化到磁极等方面均取得了好的效果,但是在许多时候,特别对于观测面起伏高差较大的情况,常常效果不佳。而且有的方法丢掉边部测点的信息较多。在这里我们介绍重磁异常换算的等效源法,这种方法可以将上述不同方面场的换算问     

Za_⊥和△g_m的八方位换算方法

本文讨论Za_⊥及△g m八方位计算方法,Za_⊥为垂直磁化下的垂直磁异常,△g_m为磁源重力异常。这种计算方法也统称为化到地磁极。经常使用的换算方法有频率域和空间域两种。 本文提出的八方位方法是空间域换算方法之一。同时还存在六方位和十二方位的方法。以上提到的三种空间域换算方法的基本原理都是相同的,并在文献中已详细的讨论。     

估计磁场源位置的倾斜角梯度算法改进

倾斜角梯度法在未知场源物性先验知识和构造指数情况下,可快速估计出场源的边界位置和深度分布,对于大面积实测磁异常的快速解释具有应用前景。笔者在Salem的倾斜角梯度法基础上,提出了一些算法改进。在导数换算上采用空间域有限差分法换算磁异常的水平导数,采用ISVD算法换算磁异常的高阶垂直导数。在剔除坏解上,采用磁异常广义倾斜角梯度的总水平梯度滤波法和场源深度有效性筛选法相结合。理论模型数据试验表明本方法简单快速,效果明显。     

关于位场垂直二次导数换算公式的频率响应

在重、磁资料数据处理中,经常进行垂直二次导数的换算,以消除区域异常干扰,突出局部异常。位场垂直二次导数换算公式的种类繁多,且使用效果不甚一样。为了评价各公式的计算精度,正确地选用公式及其环半径,已有一些地球物理工作者从各公式的频率响应特性出发讨论了这些问题[1－3],     

起伏地形二维位场上延与换算的边界单元法

本文用边界单元法进行起伏地形二维位场向上延拓与换算,介绍了方法的原理和具体的计算公式。以地形上的磁异常Z为例,换算出地形上的H,延拓出地形上部的Z及导数δz/δx、δx/δy。计算结果与精确解符合得很好。在有效异常范围内,误差基本小于1％。     

不损失点的二度磁异常延拓方法及精度

在磁异常的解释推断中,较广泛地应用了解析延拓方法。人们为了了解剖面上磁异常的空间分布特点,经常把实测场值换算到上半空间或下半空间。然而,当位场在上延或下延过程中,边缘点将会受到损失,尤其在多次延拓后,边缘点损失更为严重。     

磁源重力异常及其导数的计算方法

引言 当矿体斜磁化时,对实测的磁异常曲线进行定量解释就会发生很大困难。克服这个困难的有效办法就是把斜磁化的实测数据换算为垂直磁化的异常数据,这种换算要借助于磁源重力异常及其导数。此外,在磁场及其导数的研究中,在重磁资料的对比中,都有赖于对磁源重力异常及其导数的计算。磁源重力异常的概念是由Baranov[1] 提出来的,并由他给出了一个在殊特情况下的六方位计算方法。使用这个方法时,必须假定     

均质多面体重磁异常正演计算表达式的一致性

根据重力位、重力异常、磁位和磁异常公式,应用泊松关系,通过严格的公式推导,给出了重力位、重力异常、磁位和磁异常之间的数学关系,从而使重磁异常正演计算表达式包含于简单的规则之中,揭示出均质多面体重磁异常正演理论的一致性。     

一种二度磁异常分量换算的数字计算方法

本文从二度磁异常水平分量与垂直分量之间的基本关系式——希尔伯特积分出发,应用数字积分的梯形公式,提出了由二度磁异常垂直分量值换算水平分量的新方法。其换算公式为:Ha〔K〕=1/2π(Za〔K 1〕-Za〔K-1〕) 1/π(?)m〔j〕Za〔i〕/i-ki≠K式中 m〔i〕(?)(i=1或 n)式中 (i≠1或 n)Ha〔K〕和 Za〔i〕分别是第 K 个观测点上的磁异常水平分量值和第 i 个观测点上的磁异常垂直分量观测值;n 为磁异常剖面上观测点的个数。本文所推荐的方法具有边缘损失小、计算精度高和便于在电子计算器上实现数字计算等优点。文中还给出了水平圆柱体和有限延深薄板状体理论模型上的算例及换算结果的反演数据,以验证方法的可靠性。     

用等效源法进行磁异常换算的效果(二)

在物化探电子计算技术79年第2期中,我们讨论了用等效源法进行重磁异常换算的效果[1]。方法的基本思想是,采用最优化方法求出观测异常的场源分布,根据这个场源分布,再来计算任意点的各种场值,如向上延拓,向下延拓,曲化平(将起伏地形上的观测值换算到某一水平面上的场值),磁异常不同分量之间的换     

利用归赤道换算法解释磁法资料

磁场总强度异常的形状受磁化方向(磁倾角和磁偏角)和地磁场矢量的影响.除非磁倾角是90或是0.否则磁异常的中心就不在异常体的正上方通常利用归极换算法消除这种磁倾角和磁偏角的影响使异常中心移至异常体上方许多年前就认识到归极换算法不适用于在低磁纬度地区观察到的异常.这种方法使换算得到的异常有误差,走向也发生偏差.本文提出另一种处理方法来消除磁倾角和磁偏角的影响.这种新方法是将异常换算到磁赤道(磁倾角为0).这是一种滤波方法,重新计算磁场总强度异常,引起异常的物体有如是在水平方向被磁化.经滤波后,异常低值将位于引起异常的物体正上方为了便于解释,使倒相180,则低值异常变成     

■换算成■的效果

本文叙述了应用位场转换方法,对M512—2地磁异常垂直分量的(?)a转换为水平分量(?)a的计算,并与实测(?)a对比。讨论了有关误差问题,结果表明,在一定条件下,根据实测(?)a换算(?)a是可行的,能够代替野外实测的(?)a。一、问题的提出如果知道磁异常的水平分量(?)a,可以使一些需要(?)a的解释方法成为现实,从而提高磁法勘探的解释水平,取得更好的地质效果。     

关于磁异常化到地磁极问题的几点认识

文中阐述了磁异常化到地磁极（简称化极）的物理意义,指出了化极和换算为垂直磁化条件下的垂直磁异常的差别。在做化极处理时,实测的磁异常要完整,否则将会引起化极结果的畸变。利用教科书中给出的公式对二度体磁异常做化极处理,对化极结果解释所得的磁化强度是磁性体的总磁化强度,一般情况下不是原始的ΔＴ（或Ｚａ）异常所对应的有效磁化强度。     

基于二维经验模态分解的重磁异常分离

由于地质过程的复杂性及成矿过程的多期次叠加性,原始重磁异常往往是多种地质因素的混合信息,既包含区域背景异常信息,也包含与矿床(体)、矿化蚀变带以及隐伏岩体等与找矿密切相关的地质要素所引起的局部重磁异常.如何从复杂的叠加重磁异常中分离出具有找矿意义的局部异常,是当前矿产勘查和资源潜力评价工作中面临的难题之一.采用经验模态分解(EMD)方法来分解重磁异常,为提高分解的稳健性提出了用双调和样条插值(BSI)进行包罗面插值的新方法,并以云南个旧地区重磁数据为例,对其进行非线性多尺度分解,实现对区域异常与局部重磁异常的分离,揭示了深层次找矿信息并拓宽了经验模态分解方法的应用领域.      

地球磁异常(EMAG2)数据中海域资料质量评估——以北极地区Kolbeinsey脊南段为例

为了评估地球磁异常(EMAG2)数据中海域资料的质量问题,选择Kolbeinsey脊南段地区约193 500 km2的1∶500 000及1∶1 000 000的航磁数据进行对比分析。基于地球磁异常(EMAG2)数据,利用解析延拓方法获得了EMAG2(向下延拓4 km)及航磁异常(向上延拓4 km)。利用相关分析法对EMAG2与航磁异常(向上延拓4 km)、EMAG2(向下延拓4 km)与航磁异常进行相关分析,得到了相应的相关系数。通过对这些数据之间相关系数及差值特征的综合分析,评估了EMAG2数据中Kolbeinsey脊南段地区磁力数据的质量。研究结果表明,EMAG2数据融合了大量海域航空磁测及海洋磁测资料,在测线较密的海域,其数据质量相对较高。然而若将4 km高度的EMAG2数据换算到海平面高度的异常数据,其数据质量较低。     

卡乌留克塔格地区地质及磁异常特征

卡乌留克塔格C-12磁异常位于兴地大断裂南侧,航磁为低缓异常,地磁异常则极值较大.通过地面高精度磁法工作,确定了航磁异常在地面的具体位置及形态、范围及变化特征.根据磁异常的推断解释,C-12异常为超基性岩体引起.由此可见,磁法测量对寻找基性、超基性含铁磷岩体具有重要的意义.     

时变场中导电导磁体的边界元素法

该文研究了低频人工时变场（似稳场）中任意导电导磁矿体的边界元素法。通过对矢量位的边界积分方程，再以样条边界元素法进行离散，得到实分量的常系数方程组和虚分量的二阶微分方程组。此方程可利用时间积分的样条函数方法求解，解出适量位后，利用磁场强度与矢量位Ａ的相应关系而求得二次磁场的水平分量和垂直分量。