重磁数据稳定向下延拓的水平导数迭代法

向下延拓是重磁数据处理的常用手段,能有效地区分叠加异常、增强浅部异常,但现有向下延拓算法的计算大多是不稳定的,且易受噪声的干扰,往往造成异常形态的畸变.提出一种基于水平导数和向上延拓联合迭代的向下延拓算法,由于向上延拓和水平导数的计算是稳定的,因此该向下延拓方法可有效地增强结果的准确性和稳定性.理论模型试验表明该方法的向下延拓结果比Fourier变换计算结果更加稳定、准确,且受噪音干扰小.将该方法应用于实际数据的处理,结果显示该方法能稳定和准确地完成异常的向下延拓任务,且有效地增强了浅部局部异常.      

王浩铁矿区重磁数据处理的效果

利用重磁垂向二次导数和向下延拓等方法,分解江苏省王浩铁矿的一个水平叠加异常,提取深部盲矿体引起的局部弱异常,经正演计算和钻探验证,在预计深度找到一个有一定厚度的铁矿体。     

向下延拓在深部矿产勘探中的应用——以青海某矿区为例

勘测深部盲矿体时,矿体距离观测点较远,磁异常的相互叠加常常在地面形成具有一个磁异常中心的宽缓磁异常,分辨不清矿致异常和背景异常。通常用向上延拓、方向导数等处理的效果差强人意。笔者采用向下延拓方法处理了青海某矿区的宽缓磁异常,向下延拓深度为30倍点距,获得了精细的各矿体异常特征。通过与大量的见矿、未见矿钻孔比较,向下延拓场反映的磁异常位置与已见矿钻孔有较好的对应关系,同时指示出全部的未见矿钻孔都处于向下延拓场反映的磁异常外围。基于当前向下延拓方法的稳定计算问题已经取得突破,笔者认为针对一些大埋藏深度的宽缓磁异常采用向下延拓处理,区分叠加异常,圈定异常位置,可以获得事半功倍的效果,为钻孔布置提供佐证。同时对勘探精度也作了一些有意义的讨论。     

位场向下延拓的相关系数法

位场向下延拓的研究不仅仅是在重磁资料的解释中,而且在导航方面也有着重要的作用。位场的向下延拓归结于求解一个第一类Fredholm型的积分方程,笔者借鉴概率成像的思想,通过计算已知场与待延拓点位置之间的相关系数来求解向下延拓的方程组,实现了位场的向下延拓。对该方法进行检验:1)理论数据向下延拓10倍点距时,均方误差为理论异常最大值的0.23%;2)理论数据含均方误差为理论数据最大值的1%与5%的白噪声时,向下延拓结果的均方误差传播系数分别为1.18与0.54;3)实际资料向下延拓10倍点距时,相对平均误差为4.9%。     

冀东南部航磁ΔT低缓异常化极下延特征及其对深部找矿的指导意义

本次研究使用河北省1∶2.5万高精度航空磁测资料。通过对航磁数据进行化极、向下延拓数据处理,对冀东地区3处典型航磁异常不同深度的下延效果进行对比分析,确定磁性体顶板位置。结合钻孔验证结果,讨论了深部矿体在航磁数据向下延拓曲线上的特征,表明了大比例尺高精度航空磁测在寻找深部铁矿上有很大作用,该研究对解决覆盖区深部找矿及危机矿山等问题具有重要意义。     

迭代Tikhonov正则化位场向下延拓方法及其在尕林格铁矿的应用

解析延拓是一种广泛应用的位场处理方法,向下延拓可以压制深部地质体的影响,突出浅部异常。但是,向下延拓滤波因子是一个高通滤波器,会造成下延结果震荡,从而限制了该方法在实际资料中的应用。文中详细介绍并实现了迭代Tikhonov正则化向下延拓方法,在理论模型上将该方法与传统频率域延拓方法进行对比,表明迭代Tikhonov正则化向下延拓方法的有效性;并将该方法应用于青海尕林格铁矿区磁测资料的处理解释中,下延结果与钻探情况相符,说明在厚覆盖层的勘查区中,运用迭代Tikhonov正则化向下延拓方法能够有效地提高资料处理解释的效果。     

航磁异常深部弱信号提取技术研究

磁异常通常是地下不同深度磁性地质体产生磁场的叠加,在规范高度的航磁测量结果中,深部磁性体所产生的异常通常表现为弱而平缓,其水平与垂向分辨率均较低,在航磁资料处理解释中难以有效捕获。因此,应用适当方法提取由深部地质体引起的弱磁信息是十分必要的。本文采用精度高且稳定的位场延拓技术将航磁异常向下延拓,可以稳定增强磁异常幅度,随着延拓面与场源之间距离的减小,浅成磁信号与深成磁信号的视深度差异将增大,在对数功率谱上可以将其区分,而后可利用匹配滤波方法将浅部信号剥离,从而得到深部弱信号,同时可计算深部弱信号的视深度。     

一种基于GCV法改进的频域内位场正则化向下延拓的参数选取方法

这里给出了一种在频域内位场正则化向下延拓方法中,用于选取最优正则化参数的极小化泛函。通过不同的正则化迭代方法进行仿真数据和实测航磁数据的向下延拓实验,来验证此泛函用于选取参数的有效性,实验结果证明了此方法用于位场向下延拓是可行和有效的。     

二度磁异常向下延拓的一种方法

从二度磁异常上延理论出发,推导出一种下延方法。该法有较高的延拓精度。为减少实测数据中高频干扰给下延结果带来的影响,还讨论了磁异常的光滑问题,并推导出一个能压制高频干扰、突出低频异常的下延公式。     

视磁化率的计算及其应用

引言 航磁资料数据处理是从航磁资料获取尽可能多的地质信息的手段。诸如为了消除地表局部磁性体的干扰,区分局部场和区域场而采用的向上延拓;为了区分复杂场,了解场源分布的向下延拓;为了区分同级相互干扰的异常,圈定磁性体边界的垂向导数计算;以及为了排除斜磁化影响,简化异常解释的化磁极运算等等,都已为广大物探和地质工作者所熟悉和应用,并在航磁资料的地质解释中发挥了一定作用。视磁化率换算则     

波数域重力归一化总梯度法中的圆滑滤波因子

讨论了波数域重力归一化总梯度法,提出在位场转换的向下延拓和求导过程中分别引入圆滑滤波因子,该方法增强了计算的可靠性和稳定性,扩大了可延拓深度;通过模型实验和数据分析,确定了向下延拓和导数因子中圆滑滤波因子指数参数的取值规律,指出了欲求地质体越浅(深),采用的谐波数越大(小)的地质体质心埋深与谐波数N之间的变化规律;发现了均匀密度三度背斜模型和非均匀密度储油背斜模型与二度模型相似的G^H场等值线特征,为地质体质心埋深计算及储油构造研究提供了思路。     

随机误差在频率域的传播

随机误差在重磁数据的测量、各种数学变换处理以及正反演计算中普遍存在。重磁反演解的不稳定其原因之一是由于存在随机干扰,在向下延拓、求高次微商等的计算中,常将误差过份放大,使有用异常被严重干扰,无法辨认。作者为研究随机误差对重磁变换的影响,对随机误差在频率域中传播的情况作了深入探讨。     

计算磁性体上顶与下底深度的功率谱法

我们曾从导数的概念出发,导出了利用磁异常的垂向导数或磁异常本身的功率谱直接计算三度磁性体下界深度的方法。本文则以同样途径,除了给出直接计算二度磁性体下界深度的方法外,还导出计算二度和三度磁性体上界埋深的方法。理论模型试验结果表明,计算深度具有相当高的精度,特别是下界埋深。     

确定水平圆柱体中心埋深和走向长度的方法

在磁异常解释工作中,磁性体的中心埋深和走向长度是确定磁性体的空间位置和几何尺度的主要参量,因而是异常解释的重要内容.在用似二度体模型对三度体磁异常进行选择法计算解释时,也首先遇到如何根据异常形态合理选取似二度体模型的中心埋深和走向长度的问题.对任意形状的三度体和任意截面的似二度体进行严格的理论计算将遇到数学上的困难.作为对这一问题的初步探讨,本文对有限长的水平圆柱体进行了讨论.现提出一种根据异常的特征点——零值点的位置来确定磁性体中心埋深和走向长度     

二维位场向上延拓与向下延拓的样条函数法

本文用样条函数法实现了二维位场的向上延拓和向下延拓，向上延拓用样条函数法计算泊松积分，向下延拓用样条函数法求解第一类Ｆｒｅｄｈｏｌｍｔｍ积分方程。模型和实例计算表明，本文方法正确可行且精度高。     

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

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

将三度体磁异常转化为二度体磁异常的问题

三度体磁异常的定量解释,由于缺乏合适的方法,或是为简化计算,人们曾设想,沿观测平面上一定的方向,对实测的三度体磁异常进行线积分,而获得具有一定横截面(其形状与原三度体有关)的二度体的磁异常,从而采用简单的二度异常反演方法或利用二度量板的选择法,间接地对原来的三度体磁异常作出推断解释.这种方法的原理、具体做法、误差和实际应用,文献中已有记载.但对该法的实质还有含混不清之处,误差的分析比较片面,因此,妨碍了对该法作     

重磁位场的正则化向下延拓

传统的位场向下延拓方法,在下延至场源体顶面深度附近区域,位场将发生强烈的振荡效应,无法继续向下延拓,难以用空间位场特征确定场源体的深度等特征,当然更无法划分出垂向迭加体。本文讨论了一种新的正则化方法：利用引入位场频率,埋深与正则化因子等有关的校正函数。使得向下延拓过场源体时,场值不奇异。因此,可以达到场源体以下任意需要的深度。理论场和实测资料的正则化向下延拓成果表明,本方法能有效地反映出场源体的深度等特征,也能较准确地划分出具有垂向与水平迭加的单个场源体。     

视磁化率(SUSC)填图的理论与实践

视磁化率(SUSC)填图是对总磁场强度数据进行高分辨率增强滤波的一种方法.虽然也可以对地面磁测或一维数据进行视磁化率的计算,但用于二维航磁图的分析,更能发挥其巨大的潜力.计算视磁化率最简单直接的方法,是在波数域中进行分支线性滤波(component linear filter),需要作以下四种滤波处理:低通稳定滤波、向下延拓、化到磁极、消去几何因子.方法所依据的统计前提是:假定地     

关于不规则二度体磁异常的计算

前言地球物理文献中有许多计算不规则体的重磁异常的方法.其中大部分方法不外乎是以一系列平面所限制的模型来逼近物体,计算每个面的效应,或者是把物体划分成水平的薄板,计算这些板的效应,再在其延深的范围内求其积分.还有人在计算三度体磁异常时,将物体切成垂直于磁化方向的板,计算各垂直磁化板所产生的异常.本文所介绍的方法是将二度体沿磁化方