迭代Tikhonov向下延拓方法参数影响及应用

重磁位场向下延拓是重磁异常处理的一种重要方法。由于向下延拓具有不稳定性,一般只能向下延拓较小的深度。而向下延拓Tikhonov正则化迭代方法则改善了向下延拓的稳定性,提高向下延拓的深度。结合向下延拓Tikhonov正则化迭代方法波数域算子,分别讨论不同正则化因子、不同迭代次数以及不同延拓深度对向下延拓Tikhonov正则化迭代方法波数域算子的形态、幅度、极值点对应波数以及滤波范围的影响。通过对不同深度和不同磁化强度球体正演数据的处理,探讨正则化因子和迭代次数对延拓结果的影响。应用实际资料进行处理,评价其应用效果。     

反射地震资料全波方程的向下延拓

研制出一种压制地震记录中的多次反射的新方法。这种新方法建立在向下延拓步骤上，这种向下延拓步骤使用全波方程式（双曲形）作为向下延拓因子。我们将演示记录波场的向下延拓测绘出多次反射同相轴的射率函数的情景。本方法被成功地应用于平面波分解的（倾斜迭加）合成资料及野外资料的合个道中。     

重力归一化总梯度的稳定解法

重力归一化总梯度法是20世纪80年代兴起的一种求解重力反演的方法,其关键运算是向下延拓及求导。原有的向下延拓方法大多数是在波数域实现的,由于稳定性差,影响了归一化总梯度法的效果,从而难以推广应用。本文以原先的希尔伯特变换法求解重力归一化总梯度值为基础,采用了一种新的、稳定的向下延拓算法——迭代法,替代了原有波数域的向下延拓算子,提高了计算的精度和稳定性。模型实验和实测资料应用显示效果明显改善。     

位场向下延拓的泰勒级数迭代法

本文介绍了一种新的位场向下延拓方法—泰勒级数迭代法。首先将波数域中的向下延拓响应因子进行泰勒级数展开,取前N项和代替向下延拓响应因子进行近似向下延拓,并采用迭代法进行逐步逼近,直至观测面上的实测值与迭代m次的计算值的差值小到可以忽略。由此推导出了位场向下延拓的泰勒级数迭代法的迭代通式,证明了该方法的收敛性,分析了该方法的收敛速度。同时指出当通式中所取项数N=0时,与文献[11]的延拓模式一致。模型试验对比分析表明,当向下延拓深度较大且项数N=1和N=0迭代得到的延拓误差相等时,N=1的迭代次数远小于N=0的,且迭代的结果具有更好的保幅性。     

位场向下延拓的波数域正则—积分迭代法

本文将正则方法与积分迭代法相结合,提出波数域正则—积分迭代法用于位场向下延拓。该迭代法以位场正则向下延拓值为迭代初值,并用正则向下延拓算子处理迭代过程中的剩余谱,使得在加快迭代收敛的同时,压制噪声谱,以提高延拓的稳定性;计算了改进迭代法的最优步长,对比分析了改进迭代法的收敛性和滤波特性,同时提出采用L-曲线准则计算正则向下延拓的正则参数,并将其扩展到波数域以提高计算速度和适应波数域迭代法。基于理论模型和航磁数据的对比试验分析表明,改进迭代法位场向下延拓的精度、稳定性及收敛速度都有提高,其延拓效果优于现有广义逆法、积分迭代法、泰勒级数迭代法和Landweber迭代法。     

位场向下延拓的边界单元法

引言众所周知,位场向下延拓可以突出局部异常,放大某些在低缓异常中不够明显的异常特征,从而提高了异常解释的可靠性。因此,它在金属矿勘探和石油勘探中得到了越来越广泛的应用。但是,由于位场向下延拓的不适定性,寻找精确而稳定的算法一直是理论和实际工作中的难题之一。近些年来,不少学者从事这方面的研究,提出了一些新的算法或改进措施。候重初提出了频率域位场逐次向下延拓的方法;王延忠、熊光楚讨论了频率域利用组合滤波器实现位场向下延拓,取得了较好的效果。本文利用边界单元法原理,把位场向上延拓的泊松积分转化为位场向下延拓的第一类积分方程,然后求解。模型试算和实际应用结果都表明,这种方法有着较高的精度和稳定性。     

地震反射数据的全波动方程向下延拓

本文提出一种压制地震记录中多次反射的新方法。它是建立在应用全声波方程(双曲线形式)作为向下延拓算子的向下延拓方法基础上的。我们证明记录波场的向下延拓可绘成反射率函数曲线,而没有多次反射波。该方法成功地应用于单道平面波以分解(倾斜叠加)合成记录和野外数据。     

用克希霍夫积分作偏移速度分析

引言反射地震资料可看作由来自地下构造的源脉冲散射产生的一组波前。为了从地震资料进行地下成像,栽们采用了波动方程偏移。偏移包括向下延拓和成像两部分:向下延拓使这些波前收敛到它们的散射点上,成像使这些聚焦的波前在这些点上成像。当用于向下延拓的速度与由波前决定的地下速度一致时,波前就收敛到它们的散射点上。然而,地下速度往往是不可能预先得到的,因此通常我们只能根据地震资料估计地下速度。目前,已经研究出一些根据地震资料决定偏移速度的技术。应用若干常速加上偏移空间速度分析进行f-k域叠前偏移(Shurtleff,1984)就是其中一例。虽然f-k方法计算速度快,但反     

论电磁场的向下延拓

刀丁得ylJ特定深度上电磁场的表达式、表近水平场作向下延拓计算.需要对地明推导过程‘设,本文对此给出了一种简分析中采用了均质板状地层掉划的假论电磁场的向下延拓@J.R.Wait @挥之     

Geophysics,Vol.55,No.11,1990论文文摘

正如可以通过表示来自震源的波场辐射为一组高斯射束来产生合成地震数据那样,将记录波场表示为在地表面出射的高斯射束来进行记录数据的向下延拓。在这种情况下,当地震速度有横向变化时用高斯射束来描述地震波的传播是很有利的。高斯射束向下延拓与波动方程计算地震波的传播相适应,并且保持该传播射线路径的解释描述特征,本文描述了使用高斯束对记录波场进行向下延拓的零偏移距深度偏移方     

基于向下延拓Milne法的重力归一化总梯度法

重力归一化总梯度法需进行向下延拓和垂向导数的计算,这两个过程都会放大高波数成分,影响计算结果的稳定性。为此,将向下延拓Milne法及积分垂向二阶导数法引入重力归一化总梯度的计算。将该方法应用于估算无限长水平圆柱体模型的中心位置,计算结果准确。然后针对实测重力数据,利用本方法得到的矿洞中心埋深较基于泰勒级数展开的重力归一化总梯度法更准确,因此本文方法具有一定的实用性。     

基于凸集投影的重力数据扩充下延一体化方法

重力勘探的工区经常是不规则的，有时重力数据还存在空白区。在进行波数域向下延拓处理前，需对数据空白区进行填充。同时，为了提高处理精度，还需要将数据预扩边至快速Fourier变换所需要的长度。常规重力向下延拓方法一般仅考虑向下延拓本身，或将数据填充、扩边与向下延拓这三个步骤独立进行。本文统一考虑这三个不适定问题，提出一种基于凸集投影原理的重力数据填充、扩边、下延一体化方法，即采用Gerchberg-Saxton迭代法对数据进行插值和扩边，直至达到预定的迭代次数，然后下延。利用该方法和经典的插值、扩边、下延组合方法对理论模型和航空实测重力数据进行对比实验，结果表明该方法原理简单、操作方便，扩边结果光滑无畸变，插值和下延结果精度较高，优于经典组合方法。     

复杂地形条件下高精度重力勘探的地质效果

分别用斜顶三棱柱体及方柱体模型分近，中、远区对某山区１：５万高精度重力勘探成果进行地形改正。用常规方法对地改后的布格重力异常进行处理效果不佳，本文提出了一种二维视深度滤波法，其实现过程为：首先利用解析法对观测数据向上延拓，然后利用差分法向下延拓，重复进行上、下延拓一定的次数，再用差分法向下延拓一次，达到滤波的目的。     

山形重力异常的成因机制及消除方法

在地形起伏大的地区进行重力测量，所得重力异常中包含了和地形线性相关与非线性相关的干扰异常。后者来源于布格改正前最大地改半径之外的物质，它们是引起山形异常的主要原因。     

滤波因子V_(zz)在重力资料褶积处理中的应用

重力异常的振幅谱是主频为零并随频率的增大而呈指数衰减的曲线。为了获得反映局部构造的高、中频成份,选用具有带通滤波特性的重力位二阶导数作滤波因子对重力异常进行滤波,用比较法、切线法和对数谱法反演界面的埋藏深度,该方法具有降低区域背景对重力测量值的影响、突出微弱局部异常的作用。通过理论模型和珠江口某工区实际资料计算验证,所求界面埋深精度有明显提高。图5表2参3     

基于径向谱的位场向下延拓正则参数选取方法

Tikhonov正则算法是解决位场向下延拓不适定性广泛采用的一种有效方法,最优正则参数的选取对下延精度和计算时间都有影响.为了快速、有效地选取正则参数,提出一种基于位场径向平均功率谱的Tikhonov正则参数选取方法.该方法通过建立由径向平均功率谱确定的截止波数和Tikhonov正则低通滤波函数的截止波数的关系来确定正则参数.基于理论重力模型数据及航磁实测数据将新方法同以往普遍采用的L-曲线法、C-范数法进行了正则向下延拓效果的对比实验.实验结果表明,新方法实现简单、物理意义明确,且选取的正则参数优于另外两种经典方法.     

均衡异常的初步认识

重力勘探范畴中的校正很多,均衡异常值的正、负及高、低首先与地壳的浅部构造——山区、沉积盆地以及盆地中的隆起、凹陷相对应。均衡校正实际上是地壳厚薄不均的校正。本文以艾里假说为依据,参考地质部物探所编制的ALGOL 60地形-均衡校正程序,编制了求均衡异常的标准FORTRAN程序,计算了新疆天山及准噶尔盆地的地形-均衡校正值,成果与手算查表法作了比较,精度是满意的,通过对计算成果的分析,得出了对均衡异常认识的若干结论。     

频率域位场下延的振荡机制及消除方法

频率域中位场延拓及转换已在重磁场处理中得到了广泛的应用。但位场下延，求高阶导数对出现的振荡问题一直没有得到很好的解决，目前比较一致的认识是浅层异常源是导致下延振荡的唯一原因。