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

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

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

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

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

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

利用正则化方法计算重力归一化总梯度

针对重力归一化总梯度计算过程中向下延拓的不稳定性,将正则化因子引入到向下延拓中,提出了利用正则化方法计算重力归一化总梯度的方法。在模型试验中,对非均匀球冠模型中不同厚度油气藏进行试验分析,结果表明:选取合适的谐波数和正则化因子,可以较好识别出"两高夹一低"的储油特征,且识别能力强于常规的Fourier变换法,这说明本文方法能够提高重力归一化总梯度计算结果的稳定性和分辨率。在实例应用中,利用此方法不但清晰地显示了已有油气藏的位置,还推测了5个油气构造区。     

位场向下延拓的导数迭代法

针对向下延拓的不适定问题,提出了一种新的位场向下延拓方法-导数迭代法。从位场垂向一阶导数的定义出发,将观测面和向上延拓、向下延拓同等高度平面上的位场值近似联系起来,采用迭代法进行逐次逼近,推导出空间域位场向下延拓的导数迭代法的递推公式。考虑到空间域向上延拓积分方程实现的复杂性,对递推公式进行快速傅里叶变换,整理得到波数域中的导数迭代公式,同时证明了该方法的收敛性。模型检验和实例分析均表明迭代法相对直接FFT法具有稳定性强和下延深度大的优点。     

位场正则化下延在自然电位概率成像中的应用

通过位场正则化向下延拓方法,弥补自然电位法对于地下全空间信息反映不足的客观缺陷,改进仅通过地表数据成像而导致的深度不灵敏状况。在模型分析中,新方案的成像结果和真实模型非常吻合。对于实测资料解释,新方案的处理结果在横向上符合实测数据的异常反映,同时还能够准确的给出地下半空间的纵向信息。     

重力向下延拓的迭代法对比分析研究

向下延拓是位场数据处理和解释的基本方法,但其不适定性导致计算结果稳定性差。近年来,迭代法在提高向下延拓计算稳定性和计算精度方面上的优势备受重视。为比较各种迭代法的优缺点,对积分迭代法、泰勒级数迭代法和补偿法进行了对比分析研究。在模型试算的基础上,将三种方法应用到了江西相山铀矿田重力数据处理之中。研究表明,三种迭代法的迭代通式具有较强的相似性,只是选取的初始滤波算子不同,但初值滤波算子选取的合理性直接决定着迭代法对低频成分的转换能力和高频成分的压制能力。补偿法在向下延拓计算中具有更高的计算精度和更强的计算稳定性,因此选择补偿法具有更好的实际应用价值。     

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

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

对迭代法位场向下延拓方法的剖析

位场向下延拓迭代法的实质内容是"向上延拓而不是向下延拓"和以"迭代的结果趋近于观测值"为标准的操作过程。根据数据操作流程剖析了位场向下延拓迭代法的运行机制,得到了迭代数据在空间域的变化规律,即用迭代法将观测高度的位场向下延拓一个深度h。这实际上是通过不同高度的向上延拓来实现的。也就是说,迭代次数增加一次,涉及的上延平面就增大一个h的高度。一般地,迭代次数n与上延高度h的关系为n~(n+1)h。在空间域中,初值、上延结果、差以及每一次校正后的结果都能用满足莱布尼兹定理的交错级数表示,从而得出了迭代法能够收敛的结论;或者,以"观测高度上的实测值与计算值的差值小到可以忽略"为标准,从数学上也能证明迭代法能够收敛。数学推论和模型试验结果说明了迭代的位场初值可以任意给定。在实际操作中,迭代误差标准的影响和由于迭代误差标准不恰当可能出现不能达到迭代标准的情况,需引起注意,也值得进一步研究。     

基于层源位场的重力及其梯度数据联合相关成像

为提高重力与重力梯度数据联合成像的纵向分辨率,笔者提出基于层源位场的重力及其梯度联合相关成像方法.为体现不同深度位置地质目标体的空间分布特征,采用均值滤波法和插值-切割分离法分离位场异常,并使用向下延拓方法得到不同深度层的层源重力及其梯度数据,在此基础上实现分层重力及其梯度的联合相关成像.将本文方法应用到多个立方体的组...     

波场延拓反Q 滤波的正则化方法

基于波场向下延拓的反Q 滤波是补偿大地衰减效应的一种常用方法,但其振幅补偿的不稳定性限制了该方法在中、深层地震资料处理中的应用。针对于此,提出波场延拓反Q 滤波的正则化方法,该方法不仅较好地解决了地震波能量补偿的不稳定问题,同时对经过品质因子较低地层的深层地震波有着较好的能量补偿效果。模型试验和实例应用均表明: 与常规稳定的反Q 滤波方法相比,基于正则化法的反Q 滤波不仅具有较强的抗干扰能力和较好的振幅补偿能力,而且有效地提高了中、深层地震资料的分辨能力。     

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

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

重力场向下延拓的三阶Adams-Bashforth公式法

重力场向上延拓是稳定且收敛的过程,而向下延拓是不稳定且发散的过程。为此,本文提出一种重力场向下延拓新方法。首先,对重力场及其垂向一阶导数向上延拓,得到不同高度的重力场垂向导数;然后,基于求解微分方程的三阶Adams-Bashforth多步法,推导出稳定的向下延拓公式;最后,为验证本文方法,将其分别应用于模型数据和实际数据。理论模型试验及误差曲线表明,相对于经典下延方法——傅里叶变换下延法和积分迭代下延法,新方法三阶Adams-Bashforth公式法下延过程稳定,边界效应不明显,下延深度可达5倍点距,下延结果与真实值的相对误差更小,结果更准确。将本文方法应用于加拿大某区实测航空重力数据,得到有效且准确的下延结果,能够识别和圈定一些细小异常特征。     

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

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

位场的频率域向下延拓方法

众所周知,位场频率域向下延拓的理论滤滤算子是exp〔2πh(u²+v²)^1/2〕,其中u和v分别为沿x和y方向的频率(波数),h是向下延拓的深度。     

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

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

利用重力归一化总梯度及相位法研究断裂构造

用波数域的位场转换和快速傅立叶变换理论研究重力归一化总梯度和归一化相位;在波数域向下延拓和导数计算的滤波算子中,分别引入圆滑滤波因子,抑制了向下延拓计算中对噪声干扰的放大作用,增加了计算的稳定性;利用改进的方法,计算并阐述了4种断裂构造GH场等值线和相位曲线的变化特征;别列兹金重力归一化总梯度法中极大值理论在研究断裂构造中具有局限性,对于无限延伸断裂构造,应根据GH场等值线的走向、形状的特点,结合归一化相位曲线的转折点确定断裂空间位置.     

位场垂向高阶导数的Tikhonov正则化迭代法

在位场数据处理中,垂向导数具有重要的物理意义。其在一定程度上可以划分不同深度和大小异常源产生的叠加异常,且导数的阶次越高,这种分辨能力就越强,但通常认为高阶导数的换算是不稳定的。本文在Tikhonov正则化求位场垂向高阶导数的基础上,结合迭代法进行逐次逼近,提出了位场高阶导数的Tikhonov正则化迭代法,并且得到Tikhonov正则化迭代法的递推公式。通过对该方法的滤波特性分析可以看出,该方法计算的位场垂向高阶导数具有一定的稳定性及保幅性。模型试验和实际数据的处理表明,该方法计算结果较常规FFT求导法有更高的稳定性和实用价值。     

确定接触面的位场矢量方法

本文叙述了一种定量计算接触面参数的位场矢量方法.该方法以位场的频谱表示为基础,得出重磁异常函数向下延拓的表达式;再借助于归一化全梯度的非线性变换函数,可以定量地计算接触面的空间位置、埋藏深度及倾斜方向.为了克服位场函数在向下半空间延拓过程中产生的不稳定性,采用了径向频率随深度增加而发生位移的圆滑因子.该方法的优点在于能够得到稳定的计算结果、精度高、能直接地得出接触面在深部的空间参数.本文用几个模型试算和实例证明了这些优点.     

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

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