应用计算机“推导”位场圆滑滤波及二次导数公式的方法

根据台劳级数展开式,将窗口内网格结点位场数据在计算点展开,且对应点两两相加,由最小二乘法求得窗口内网格结点的圆滑滤波算子和二次导数算子,将此算子与观测的位场数据褶积,就可得到位场圆滑后的数据及水平、垂直二次导数。     

基于最小平方滤波的面波压制方法

在地震勘探中,强振幅、低频率的面波对于后续地震资料的处理和解释会产生重要影响。基于最小平方滤波的面波压制方法是一种新的面波压制方法,它与一般的频率域滤波方法不同。该方法的主要原理是利用扫描信号与滤波因子的褶积来预测面波,其中滤波因子的求取需要建立一个目标泛函。当目标泛函取得最小值时可从地震数据中求取滤波因子进而预测出面波,再将预测出的面波作为输入进行第二次面波预测,从地震记录中减去提取的面波,就得到去除面波后的地震数据。采用非线性的扫描信号提取面波,并对扫描信号起止频率进行了分析;通过对合成地震记录与煤田实际地震资料的处理,充分证明该方法能有效地压制面波。     

蓝色滤波及其应用

反褶积是通过压缩地震子波以提高地震数据纵向（时间）分辨率的过程。地震数据可以用褶积模型来表示,即地层反射脉冲与地震子波的褶积。而地震脉冲（反射系数）包含了石油勘探和地震研究中必不可少的有关大地地球物理性质的信息,因此,从地震数据中得到反射系数就变得非常有意义。蓝色滤波通过对反射系数有色成分的模拟,得到与反射系数有色成分有关的蓝色滤波算子。然后对经反褶积处理后的地震数据进行褶积运算,就能获得比常规反褶积方法好得多的结果。通过对合成地震记录和实际地震数据的试验,证明了本方法在高分辨率处理方面具有较好的效果。     

水平层电测深正演计算的改进

自1971年Ghosh首先用线性滤波法计算电测深正演曲线以来,电测深曲线的计算大为简便。我们不再重复用线性滤波法计算电测深曲线的一般原理,仅针对本文给出的两个程序作一些必要的解释和说明。线性滤波法计算电测深曲线的一般公式据文献资料(1)、(2)、(4),用线性滤波法计算电测深曲线的褶积公式为考虑到滤波器的有限长度及位移因子,则有     

一种非线性滤波方法的试验效果

在重磁资料解释中,往往需要对异常曲线进行圆滑处理或分离后再作进一步的解释。通常采用的是所谓线性滤波,即设计一种响应或在空间域作褶积或在频率域作乘积。由于使用这类滤波器,要消去的干扰不能全部去掉,而保留的有用信号受到畸变,滤波有一定的局限性。为此,有人提出一种非线性滤波方法及位变滤     

频率空间域显式偏移

频率空间域显式偏移利用有限长度的褶积算子，通过对不同的频率进行空间褶积来实现偏移成像，且具有简洁、高效、便于扩展等特点。褶积算子一般通过泰勒展开或非线性优化来求取。针对以往方法进行分析的基础上，提出了一种新的褶积延拓算子计算方法，通过对模型资料和实际资料的处理，表明该方法是可行的。     

浅谈物探问题的TI-59计算器程序编制方法

本文以四个典型的问题:规则形体正演、褶积、匹配滤波因子求解、物性参数的统计分析来介绍TI—59计算器的编程方法。一、规则形体理论场的正演计算物探电、磁、重理论模型位场正演计算的共同特点是,需要用同一个公式反复计算不同座标点上的理论场值,这在程序设计上是属于循环问题。TI—59拥有两种(组)组织这类循环的指令,一为计数转移指令,一为比较转移指令。 1.用计数转移指令DSZ编正演计算程序     

多次波贡献道集优化在波场延拓压制多次波中的应用

基于波场延拓的多次波预测和减去技术是压制水层多次波的有效手段。常规波场延拓方法中,由已知的原始波场通过与水底各点的格林函数褶积来预测地震数据中的多次波,但在进行褶积相加时只有稳相点附近的能量才对求和结果有影响。提出了一种基于高精度拉东变换的稳相点拾取方法,是利用拉东域中的低曲率信息提取稳相点附近的同相轴,利用KL变换突出水平同相轴,得到优化的多次波贡献道集,对其进行叠加,能够很好地解决波场延拓预测的多次波存在假象的问题,并采用贝叶斯匹配减去技术将预测的多次波去除,提高了预测的多次波的精度,在同相轴两端出现的假象得到了很好的消除,理论和实际数据的应用结果验证了该方法的有效性和可行性。     

用于位场处理中的优先操作滤波器设计

基于匹配滤波和优选延拓的思想，提出了一种改进的优先延拓滤波器，与Ｐａｗｌｏｗｓｋｉ的优先延拓因子相比，它解决了对位场信号的中等频段成分优选延拓的问题。在此基础上，抽象出一种优先操作因子，这样就可以对位场信号的不同频段成分同时进行不同的处理。     

应用计算机“推导”位场水平导数公式的方法及程序

根据台劳级数展开式，将窗口内网格结点位场数据在计算点展开，且对应点两两相减，由最小二乘法求得窗口内网格结点的水平导数算子，将此算子与观测的位场数据褶积，就可求得位场的水平导数。     

含透射衰减的地震反射信号褶积模型在波阻抗反演中的应用

研究了地震波在层状地质模型内透射衰减对地震反射信号振幅的影响,给出了含透射衰减的褶积模型,为基于模型的波阻抗反演提供了一种新的正演方法。采用含透射衰减的褶积模型做正演运算,能够校正地层透射对波阻抗反演结果的影响,尤其能够消除地层的屏蔽现象。实际的资料处理表明含透射衰减的褶积模型有较好的应用价值。     

小波滤波方法应用效果浅析

为降低利用NaI(Tl)闪烁探测器得到的γ能谱数据因统计涨落而产生误差,又保留谱峰的全部重要特征,以便谱的分析,必须对实测伽马能谱数据进行光滑。文章研究了变换能谱光滑的各种方法:算术滑动平均法、重心法、多项式最小二乘拟合法、离散函数褶积滑动变换和小波滤波法等。采用理论计算,结合实测谱线,通过峰面积标准偏差,与H/FWHM的比较,对各种方法进行评价。结果表明采用小波滤波方法实现伽玛能谱的滤波处理,优于其它滤波效果。     

基于修正褶积模型的地震反演方法

在地震反演中,地震信号可以认为反射系数与地震子波的褶积,在此过程地震子波振幅保持不变。这里提出一种修正的褶积模型,将透射损失融入褶积过程中,表现为子波振幅的衰减,对于砂泥岩互层介质,修正褶积模型更加真实地刻画地震信号形成过程。基于修正褶积模型建立新的反演目标函数,反演得到的反射系数消除了透射损失的影响。与基于常规褶积反演得到的反射系数相比较,此方法保证了更多有效信息的体现。实际资料应用表明,基于修正褶积模型的地震反演方法提高了反演结果的合理性。     

利用Matlab函数识别沉积物中的米兰柯维奇旋回信号

针对旋回分析中需要对数据进行的各种处理,利用Matlab中集成的滤波（filter）函数、插值（interp1）函数、求均值（mean）函数、去线性趋势（detrend）函数、求方差（variance）函数及谱分析（pwelch）函数等的组合,编制了识别沉积物中米兰柯维奇旋回信号的程序。选择窗口宽度(window)约等于有用信号中最低频率信号的波长、交叠点数(noverlap)为窗口宽度的50％、傅里叶变换的数据点数(nfft)大于全部采样数据的长度且为2的最小整数次幂时,可获得较好的谱估计结果。其不足之处是不能估计背景的红噪声谱及给出谱峰值的置信水平。     

一种三维重、磁异常的向上延拓的近似计算

本文利用立体角求和的方法对向上延拓进行近似计算。方法假设边界平面上第一小元域内的场值对延拓值的贡献与小元域对延拓点所张的立体角有关。求出立体角就能够求出延拓系数。实例计算表明减少了边缘的损失,其结果可与原有的积分延拓方法相比拟。     

区域重力资料收录、整理、成图自动化程序简介

在1/20万区域重力测量工作中,宁夏物探队利用PC-1500机建立了一套适合野外现场进行区域重力资料收录、整理、处理以及成图自动化的生产程序。一,实测重力值的记录和计算程序携带装有该程序的PC-1500主机和重力仪到野外现场,一切观测要素直接输入到机内记录保存,对重复观测的读数,机内可自动求取平均值。工作日期和观测时间可不输入,由主机时钟控制自动给定。每当读完一个观测点后,主机自动关闭,待下一个点读     

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

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

最基本的三种混合信号分离法的综述

一、信号(波形)的混合形式经常遇到的信号混合形式,最基本的有下列三种形式:(a)加性:如 x(t)=s(t)+n(t)(b)乘性:如 x(t)=a(t)×c(t)(c)褶积性:如 x(t)=b(t)*r(t)当然,更复杂的混合形式可以通过以上三种运算的复合得到。滤波含有把信号分离的意思,即在混合信号中分离出所需要的信号分量,而同时压制其它分量。当然,这里所说的需要的分量也是相对的,由实际问题而定。     

叠前非局部平均滤波压制随机噪音

非局部平均滤波方法的去噪性能优异,但其在地震资料处理中的应用刚刚起步。该方法利用数据具有的结构冗余,以包含局部结构的小窗口或邻域为单元,利用局部结构相似性进行加权运算,增强有效信号,压制随机噪音。针对叠前地震资料数据量大、噪音背景强、局部结构简单;原始非局部平均算法对每一点滤波,需要对数据体内所有点计算权系数后进行加权计算,计算量大,对强噪音背景适用性差等不足,对原始非局部平均算法进行了改进,主要包括:基于速度谱的搜索窗口分割;基于梯度域奇异值分解的局部结构相似集选择方法;基于相似集大小的自适应滤波参数选择方法。试验结果表明,该方法改进后对于叠前地震数据的随机噪声具有较好的压制作用。      

边部无损最小二乘圆滑问题——关于地质上应用程序优化问题的讨论之一

一、数据圆滑处理工作中存在的问题 在地质工作中,特别在物化探工作中,为了消去局部干挠,突出有用异常,对观测数据进行分块最小二乘圆滑处理的方法,早已被广泛采用,并且收到了较好的效果。但是直到目前,在进行数据圆滑处理时,不管是手算还是电算,都是只计算中心区的圆滑值,边部要损失相当一部分（见图1）。为了说明边部损失问题的严重性,可做如下统计:设有M=100×100个数据的正方形网格点。用不同圆滑半径N（以点距为单位）的圆滑方法,造成边部损失的点数m和损失率δ=m/M如下表。从表中可以看出,随着圆滑半径的增加,边部损失很快增大,当N=3时,δ已达11.64%,若增大到N=6,δ将达到22.56%。另外,统计还可表明,如果原始数据的总数小于100×100,损失率还会增大。由于边部损失,不但使资料不能充分利用,造成浪费,而且还可能使资料残缺不全,给解释推断带来困难。