地震勘探数据处理的一种算法——分布参数系统最优理论的应用

本文将根据地震数据识别地层结构问题归结为分布参数系统的最优化问题,利用脉冲变分原理给出了一个新的算法,此算法与文献[4]中算法相比较,计算效率要高。     

地震勘探问题的一种新模型及其算法

本文针对地震勘探中对地层介质的识别问题,提出了一种新的数学模型:用双曲型方程Cauchy问题代替混合问题。得到一种逐层递推的快速反演算法,取得了较好的计算效果。     

地震勘探问题的一种新模型及其算法

本文针对地震勘探中对地层介质的识别问题,提出了一种新的数学模型:用双曲型方程Cauchy问题代替混合问题。得到一种逐层递推的快速反演算法,取得了较好的计算效果。     

一种适于近场理论地震图宽频域计算的矩阵分解算法

理论地震图方法在地震震源过程的研究中得到了广泛的应用。为了研究震源过程的细节,必须利用近场地震资料的高频信息,但是在应用Haskell矩阵法计算近场理论地震图时,格林函数的高频成分的数值不稳定性是计算宽频域理论地震图的一个基本困难。因此,目前仅有不超过10Hz的算例。本研究采用Haskell矩阵的一种新的分解组合形式,在数值计算时有效地避免了在计算过程中数值结果的溢出,实现了近场理论地震图的宽频域计算。同时,由于利用了矩阵运算中的解析关系,减少了运算的次数,从而提高了计算的速度。根据本算法建立的Fortran程序,在Univoc-1100计算机上进行了数值检验。结果表明,计算得到的格林函数至少在0—40Hz频率域范围内仍保持良好的数值稳定性。本文给出的这种Haskell矩阵分解组合形式的矩阵元素简单,且具有一定的对称性。对于这种分解所包含的物理意义,将有待进一步深入研究。     

一种适于近场理论地震图宽频域计算的矩阵分解算法

理论地震图方法在地震震源过程的研究中得到了广泛的应用。为了研究震源过程的细节,必须利用近场地震资料的高频信息,但是在应用Haskell矩阵法计算近场理论地震图时,格林函数的高频成分的数值不稳定性是计算宽频域理论地震图的一个基本困难。因此,目前仅有不超过10Hz的算例。本研究采用Haskell矩阵的一种新的分解组合形式,在数值计算时有效地避免了在计算过程中数值结果的溢出,实现了近场理论地震图的宽频域计算。同时,由于利用了矩阵运算中的解析关系,减少了运算的次数,从而提高了计算的速度。根据本算法建立的Fortran程序,在Univoc-1100计算机上进行了数值检验。结果表明,计算得到的格林函数至少在0-40Hz频率域范围内仍保持良好的数值稳定性。本文给出的这种Haskell矩阵分解组合形式的矩阵元素简单,且具有一定的对称性。对于这种分解所包含的物理意义,将有待进一步深入研究。     

地震资料处理中的聚束滤波方法

聚束滤波方法的基本原理,包括信号与噪音模型、设计准则和所获得的聚束滤波器．本文所用的自适应聚束滤波器是根据具有参数化的动校正量、振幅及相位随偏移距变化（ＭＶＯ、ＡＶＯ及ＰＶＯ）的一次反射和多次反射信号模型而设计的．但在实际应用中ＰＶＯ通常不能得到证实．结果分析能为进一步的从ＡＶＯ及ＰＶＯ获取岩石类型信息的专题研究提供资料．人工合成资料的例子给出了参数化的自适应聚束滤波的实施细节和所设计滤波器的响应特性．实际资料的例子表明数据自适应聚束滤波在叠前共中点道集地震多次波消除上比Ｒａｄｏｎ变换方法更灵活、有效．类似于其它叠前处理过程,自适应聚束滤波的优越性在信噪比较高的资料上体现得最为明显．     

地震资料处理中小波函数的选取研究

本文给出了常见地震子波的一个模拟公式,可以很好地模拟零相位及混合相位子波,在一定意义上也可以近似模拟最大相位及最小相位子波．模拟出的子波加上适当的修正项后满足允许条件,可用作小波函数．与Ｍｏｒｌｅｔ小波类似,在实际应用中这些修正项在一定条件下可以略去,文中对Ｍｏｒｌｅｔ小波作了改造,使其能更好地适应于地震资料处理．研究了反射波能量及噪声等干扰波在时间－尺度域的分布特征与所选基本小波的关系．提出用地震子波（或与地震子波相近的函数）作为基本小波,对地震资料进行去噪及分频解释的方法．最后用实例证明方法的有效性．     

地震资料处理中的聚束滤波方法

聚束滤波方法的基本原理,包括信号与噪音模型、设计准则和所获得的聚束滤波器．本文所用的自适应聚束滤波器是根据具有参数化的动校正量、振幅及相位随偏移距变化（ＭＶＯ、ＡＶＯ及ＰＶＯ）的一次反射和多次反射信号模型而设计的．但在实际应用中ＰＶＯ通常不能得到证实．结果分析能为进一步的从ＡＶＯ及ＰＶＯ获取岩石类型信息的专题研究提供资料．人工合成资料的例子给出了参数化的自适应聚束滤波的实施细节和所设计滤波器的响应特性．实际资料的例子表明数据自适应聚束滤波在叠前共中点道集地震多次波消除上比Ｒａｄｏｎ变换方法更灵活、有效．类似于其它叠前处理过程,自适应聚束滤波的优越性在信噪比较高的资料上体现得最为明显．     

地震资料处理中小波函数的选取研究

本文给出了常见地震子波的一个模拟公式，可以很好地模拟零相位及混合相位子波，在一定意义上也可以近似模拟最大相位及最小相位子波．模拟出的子波加上适当的修正项后满足允许条件，可用作小波函数．与Ｍｏｒｌｅｔ小波类似，在实际应用中这些修正项在一定条件下可以略去，文中对Ｍｏｒｌｅｔ小波作了改造，使其能更好地适应于地震资料处理．研究了反射波能量及噪声等干扰波在时间－尺度域的分布特征与所选基本小波的关系．提出用地震子波（或与地震子波相近的函数）作为基本小波，对地震资料进行去噪及分频解释的方法．最后用实例证明方法的有效性．     

浅层地震勘探资料处理中的速度分析参数选取

浅层人工地震勘探是当前城市地震活断层探测的主要方法,勘探精度不仅受激发、接收等因素的影响,还与地震数据处理的精度密切相关。速度分析是地震勘探资料处理中的关键环节,它的准确与否对动校正、共深度点叠加、叠后偏移及时深转换等都将产生影响。针对浅层地震资料目的层较浅、叠加次数少和城市背景干扰大等特点,文中对速度分析的参数选取问题进行了研究。首先通过对相邻CDP道集内各接收道炮检距的分析,以及单个和多个CDP道集的对比,提出了用于速度分析的CDP大道集的抽取原则。并通过对不同的CDP大道集、速度扫描间隔、计算时窗长度、动校切除比例等所获得的速度谱的对比分析,对速度谱扫描过程中各参数的选取提出了建议,并对速度分析及时深转换过程中值得注意的问题进行了讨论     

地震台站观测资料微机管理实用程序

本文简要地论述了台站观测资料微机管理实用程序的设计原则、结构及其特点，给出了操作流程示意图和简要说明。该程序使用效果较好，使观测资料处理及时，月报输出快速准确，减少了工作人员的工作量。     

反假频非均匀地震数据重建方法研究

研究基于Fourier变换的数据重建方法,既能进行非均匀采样数据重建,又可以去除空间假频. 将不规则采样数据重建问题归结为信息重建的地球物理反演问题,采用最小二乘方法从观测的稀疏或不规则数据反演模型空间完全信息. 在求解信息重建反演问题时,引入DFT 加权范数规则化策略,采用预条件共轭梯度法（PCG）求解,保证解的稳定性和收敛速度. 处理线性同相轴假频问题时,根据采样定理,引入线性预测方法,采用Yule Walker方程由带限信号的无假频低频功率谱预测高频功率谱,达到反假频目的. 本文研究了均匀采样数据内插,非均匀采样数据重建,非均匀分布高频信息重建等方面问题,数值试验取得较好效果.     

TRACE INTERPOLATION IN SEISMIC DATA PROCESSING *

Trace interpolation of spatially aliased seismic data is used to improve the quality of multi-trace processing, especially that of migration. We consider non-rectangular sampling lattices. After stating the real sampling requirements for (two-dimensional or 3-D) digitized signals, the paper explains geometrically the effect of migration and of the 2-D dip filtering upon the aliased events, as well as our ability to interpolate traces from undersampled seismic data.     

MATLAB在地震信号处理中的应用实例

实现数据的可视化是地震监测和研究工作中最重要的问题之一。本文介绍了MATLAB在信号处理工作中的主要特点,结合实例以代码的形式介绍了几个台站工作中最常遇到的数据可视化编程方法。通过分析指出MATLAB实现地震数据可视化和其它实现数据可视化方法相比具有方便灵活,图形质量高的特点。