延迟期望输出对有限数据最小平方滤波误差的作用

给定有限输入x和期望输出d,存在唯一的m 1长的最小平方滤波因子h~(m)=(h_0~(m),h_1~(m),…,h_m~(m)),为简单起见,记上式为q_m~2。一般q_m~2不趋于0(m→∞),只有当x为最小延迟时才有q_m~2→0(m→∞)。这样就限制了最小平方滤波的作用。但在地震数据处理中,有时允许期望输出d延迟其一定出现位置,即以x和为输入和期望输出,当最小平方滤波因子长度为m 1时,对应的误差记为q_m~2(S)。本文的结论有q_(2m)~2(m)→0,(m→∞)。更进一步有其中k>1的任何实数,[km]为km的整数部分。     

线性离散卡尔曼滤波及其在地震方法上的应用

众所周知,地震记录通常被描述为地震子波与反射系数序列的褶积。我们所说的反褶积处理就是要从地震记录中求取反射系数序列。它从图形上表示为一列脉冲,如能找到这些峰值就等于找到了它们的振幅和发生的时间,也等于提供了地下界面的几何和地质的信息。最常规的反褶积法是基于维纳滤波的理论。一般它要求时不变和平稳的条     

利用波形特征参数识别和预测储层含油性

地震反射波中隐含着大量地下地质信息，如波形的振幅、频率、速度、相关等，反映了地下和油气藏等特征。这些特在油藏抛述中是储层横向预测的重要手段。在勘探初期资料不多的情况下，进行综合研究评价，使本方法简单、实用。     

从鸣震记录估计海底反射系数的一种新方法

在利用三点式滤波或反向消去法压制浅海、深海鸣震时,海底反射系数的估计有重要意义。过去的方法是假定无鸣震记录的自相关函数γ_(yy)(·)满足γ_(yy)=0,|j|≥τ/2。τ为垂直鸣震周期。这个条件相当苛刻。本文的方法大大减弱了这个条件,同时对地震子波和地层反射序列不加限制。以γ_(yy)(·)、γ_(xx)(·)和ξ_(?)分别表示无鸣震记录自相关函数、鸣震记录自相关函数和海底反射系数,则只要存在n,使γ_(yy)(j)=0,|j|(?)n,就有由上两方程即可给出海底反射系数的估计值。但由于利用二级鸣震记录求海底反射系数是解一元二次方程的问题,一般有两个不同的解,因此,即使都是实根也难确定哪个根为所求的海底反射系数。为此,我们用满足该方程的两个方程联立求解,得到满意的结果。n的确定也是较困难的,我们介绍两种方法:一种是用统计的方法;另一种是根据鸣震记录的数学模型都为ARMA模型来确定的。     

预测反褶积用于消除深海鸣震的一种新途径

过去,在处理深海域的地震资料时,一般是采用“反向消去法”或“深水波形预测反褶积法”来消除鸣震。但是,当鸣震和反射波互相干涉时,这些方法的效果不理想。为了改善这类资料的处理效果,本文提出一种新的“深海预测反褶积法”。与上述两种方法比较,这种方法所要求的前提条件少,应用的范围更广泛。理论分析和模型试验的结果表明,这种方法的效果是令人满意的。     

最大熵谱估计

谱估计是数据处理的一种重要手段,初期谱估计主要用于维纳滤波设计滤波器的。近来谱估计已成为信息提取的一个有力工具,例如有些曲线(波形或信号)外形不同但却有相近的谱,而另一些曲线则其谱差别很大,通过谱估计我们可得到组成曲线的最佳频率成份,进而可以求出相应的振幅谱和相位谱来。这样,将各种曲线族通过谱估计转化成另一种更简明的形式,就能更便于进行分类识别。我们记录到的各种地震曲线就是这样,对有油气和无油气的地震曲线其中也需要用谱的特征加以分类识别。 另外,七十年代中期把最大熵谱估计方法引用到地震数据处理的反褶积方法中来,我们称它为最大熵反褶积,此法前提条件少,效果也不错,同时在此基础上我们能求取地震子波等。因此对最大熵谱估计方法的研究、改进和应用,我们认为是很有价值的。     

希尔伯特变换子波反褶积

在具有常规脉冲反褶积(DECON)相同条件的前提下,希尔伯特变换子波反褶积是利用多道振幅谱的几何平均值求取子波,并对所求的子波作各种整形,然后作反褶积,从而有利于提高地震资料的横向连续性和分辨率.试验结果表明,效果是显著的,并得到了有关专家的好评.同时,希尔伯特变换法利用了实因果序列的性质,所以提取子波方法简单,实现方便,且一个反褶积因子可应用于许多道,节省了求因子的时间。本文介绍了希尔伯特变换提取子波的理论推导和几个实现的方案,子波整形及反褶积的方法,并给出了实际资料处理的实例。     

最大熵谱估计中振幅谱和相位谱的求法

最大熵谱估计克服了窗函数对传统谱估计的不利影响,特别在数据短的情况下,它比传统谱估计优越得多。但最大熵谱估计所确定的信号最佳频率成分,其对应的振幅值并不反映原来振幅值,且随自回归模型的阶数的不同而变化。这是最大熵谱估计法的最大缺点,也是此法目前还不能很好地用于地震数据处理的主要原因。本文介绍两种由得到的最佳频率成分来恢复原振幅值和相位值的方法。