基于稀疏离散τ-p变换的叠后地震道外推

这里研究的是一种基于稀疏离散τ-p变换的地震道缺失重建方法。该方法根据叠后零偏移距剖面,在局部时窗内可以看作是一系列线性同相轴的组合,使用稀疏离散τ-p和预条件双共轭梯度算法进行地震道外推,使空间方向的缺失道得到恢复。τ-p变换由于信息不足,有限的孔径和离散等因素,导致在τ-p域的结果存在假象、不准确,而在τ-p域进行稀疏对于重建缺失的信息是非常有用和必要的。同时,还可使用迭代的预条件双共轭梯度算法进行计算。     

τ-p变换的两种实现方法

τ-p变换是一种经典的函数投影变换方法,在地震数据处理中有着广泛的应用,可以根据地震数据射线参数的差异,实现信噪分离、地震道插值、平面波分解等。但是,由于τ-p变换的精度和分辨率受到相应数学反问题的限制,在保证较高变换精度的前提下,τ-p变换计算速度的加快,以及相应滤波算子的设计等都值得研究。本文提出了基于径向道变换的τ-p变换方法和基于斜率分解的τ-p变换方法,阐明了两种方法的基本原理并与常规τ-p变换方法进行对比。基于径向道变换的τ-p变换方法利用快速Fourier变换和径向道变换,能有效减少τ-p变换的耗时;基于斜率分解的τ-p变换方法运用斜率分解算法,能进行高分辨率τ-p变换,并且提供冗余变换维度（τ-x-p域）,使滤波器算子的设计更灵活。数值实验结果表明,两种实现方法分别在计算速度和重构精度上优于传统τ-p变换算法。通过面波噪声压制的实际数据处理和比较,证明本文提出的两种τ-p变换实现方法可以为实际处理提供更加有效和灵活的实施方案。     

利用τ-p变换技术实现多波波场分离

通过研究多波地震资料空间三个分量的波场特征及在τ-p域的特性,采用τ-p变换技术,编制了波场分离程序;选取适当的空间和时间τ-p参数,对模型和实际多波地震资料进行波场分离,取得了满意的效果。     

去假频高分辨率τ-p变换法波场分离

τ-p变换是消除相干噪声、波场分离的重要方法之一,但常规τ-p变换存在算子假频和端点效应,分辨率太低,无法直接使用。本文提出一种新的非迭代高分辨率τ-p变换方法,可同时压制假频和端点效应,提高了变换域的分辨率。通过理论模型记录与实际资料的试验表明了该方法的有效性。     

面向海上地震资料的τ-p域保幅处理方法

AVO分析和反演需要地震资料处理时最大限度地保持振幅相对关系,然而常规处理中球面扩散补偿、拉伸畸变和多次波均会影响振幅相对关系,海上地震资料尤为严重。提出一套面向海上地震资料的τ-p域保幅处理方法。该方法通过τ-p域无需球面扩散补偿且多次波周期性增强的特点,减少球面扩散补偿与多次波对振幅相对关系的破坏,通过τ-p域拉伸因子非时空变换的特点减少动校拉伸产生的频率和能量畸变。最后通过挪威海上实际资料的处理验证了方法的合理性与可行性。     

基于稀疏变换的地震数据重构方法

缺失地震数据重构恢复是后期地震资料处理取得良好效果的前提.笔者通过研究稀疏变换(F-K变换、Curvelet变换)与最近流行的压缩感知理论,将两者结合起来,建立基于稀疏变换的地震数据重构模型.F-K变换是将地震数据由时间—空间域变换到稀疏域频率—波数域,Curvelet变换由于其良好的方向性、局部性以及各向异性,能够将地震数据进行更优的稀疏表达.基于重构模型,分别采用这两种稀疏变换对地震数据进行重构,并且比较两者的重构效果,证实Curvelet变换重构效果优于F-K变换.最终通过Marmousi2模型以及实际地震资料处理分析,证明该重构模型的正确性和有效性.     

τ—p变换在浅层地震勘探中的应用

数据处理是地震勘探关键所在。如何压制干扰波，突出有效波是资料处理中主要课题之一，ｔ－ｘ域中各种波相互干扰，τ－ｐ域中各种波相互分离，给资料处理提取有效波压制干扰波带来极大方便。     

基于信号重构的地震道插值

在简要回顾几种插值方法的基础上 ,提出了基于信号重构的地震道插值方法。这种方法无需任何地质和地球物理假设 ,主要用于恢复缺失的地震道。本方法的基础是傅立叶变换理论和最小平方反演技术。文中给出了原理、算法、实现步骤及简要的流程框图。如果对变换域的参数进行规则采样 ,则该方法可以用高效、快速的算法来实现。本插值方法适用于规则和非规则采样的数据。文中给出了理论合成数据和实际数据的试验结果 ,并对结果进行了讨论。同时 ,本方法对于 AVO分析、DMO处理、三维叠前偏移成像等研究也将是有益的     

τ－P变换在煤田地震勘探资料处理中的应用

本文简单介绍了τ－Ｐ变换的基本原理，以及在地震勘探资料处理中的应用范围，着重介绍了利用τ－Ｐ变换来压制线性干扰波的效果。     

τ-p变换去除金属矿地震资料中的线性干扰

金属矿地震勘探资料中通常含有多种不同类型的噪声,有效地消除这些噪声,提高地震信号的信噪比是地震数据处理的主要内容之一。二维滤波(F K滤波)是消除低速干扰的一种有效手段,它可以消除视速度较低的面波及其它干扰波,但它只适用于地层倾角较缓的地区。对于复杂条件下的金属矿效果不佳。经过多次实验,用τ-p法压制线性干扰,效果显著。     

基于压缩感知理论的缺失地震数据重构方法

压缩感知是一种新的理论,它打破常规尼奎斯特-香农采样定理的制约,利用信号的稀疏特性或可压缩特性,用较少的数据即可重构恢复完整的信号。建立了基于压缩感知理论的缺失地震数据重构模型。首先在与稀疏变换不相关的测量矩阵基础上引入一种约束矩阵,使地震数据的缺失满足或接近高斯随机分布;随机缺失的地震数据变换到稀疏域会产生很多与有效信号不相干的随机噪声,接着通过一种新的自适应阈值迭代算法可以很好地消除稀疏系数中的随机噪声干扰,经过逆稀疏变换即得到重构后的地震数据。Marmousi 2模型测试及实际地震资料处理均验证了该方法的可行性和有效性。重构缺失地震数据取得了较好的效果。     

基于抛物线方程的坐标拉伸τ-p变换的P-P、P-SV波分离

分离P-P波和P-SV波最常用的是τ-p变换法, 但是在τ-p域中, P-P波和P-SV波常有重叠部分, 不能很好分离。地震勘探中, 当炮检距小于界面深度时, P-P波和P-SV波的时距曲线方程都可近似为抛物线方程。针对这种情况, 提出了一种基于抛物线方程的坐标拉伸τ-p变换的P-P、P-SV 波分离方法, 该方法首先对地震记录沿空间方向上进行坐标拉伸, 再进行τ-p正变换, 在τ-p域中分离P-P波和P-SV波, 其次对P-P波和P-SV波分别进行τ-p反变换, 实现P-P波和P-SV波的分离。与常规τ-p变换法对比分析表明, 利用该方法能很好地将合成地震记录中的P-P波和P-SV波分离。     

τ-P变换及向下延拓方法在声纳浮标折射地震测深速度反演中的应用

地震折射波场Φ（ｘ，ｔ）通过τ－ｐ变换得到倾斜迭加场φ（ｐ，τ）。从φ（ｐ，τ）中获得τ（ｐ）函数，并将其向下延拓到深度—慢度剖面，从而反演出横向均匀的速度结构。该技术能充分利用折射波场的丰富信息，在迭代反演中不仅能简单快速地得到反演结果，而且所得解具有良好的稳定性，较少受主观因素影响。将该方法应用于海洋声纳浮标地震折射测深数据处理，得到较为理想的结果     

用MATLAB 实现地震数据的小波变换

小波变换借助于时－频局部分析特性，已经成为地震资料处理中的一种重要工具。作者在阐述其基本原理后，讨论了如何选择小波基（函数），并利用MATLAB语言及小波工具箱，实现了SEG－Y格式地震数据文件的小波分析，取得了令人满意的效果。     

面波压制方法与道间距关系的数值分析

野外无检波器组合的采集方式使得面波的压制成为单点高密度地震数据处理的关键。F-K滤波、τ-p变换是两种常用的压制面波的方法,这两种方法压制面波的效果与道间距有关。对面波速度、频率、射线参数与道间距关系的数值分析发现,道间距越小,压制面波的效果越好;大道间距会产生假频而使得面波难以有效地滤除。这从另一侧面也证明这两种方法可以更好地压制高密度地震数据中的面波。     

抗噪性加权抛物线Radon变换的地震道重建

叠前CMP道集中,走时为双曲线的地震反射同相轴经过部分动校正( NMO)之后变为近似抛物线,当道集中缺失某些地震道时,可利用抛物线Radon变换通过迭代计算进行地震道的重建和恢复. 笔者基于抛物线Radon变换的基本原理,利用加权抛物线Radon变换的计算方法,进行叠前地震道的恢复和重建,并且在数据域对其加权系数进行改进和优化,模型数据和实际资料的试算结果表明,改进后的计算方法即使在较强的随机噪声下依然可以稳建同相轴,且避免了常规方法在重建过程中产生虚假同相轴的缺点,与原始的加权抛物线Radon变换的计算方法相比,新方法具有压制随机噪声的优点,即具有较好的抗噪性,计算精度也有所提高,且算法稳健、高效,具有良好的应用前景.     

基于双曲Radon变换的混采地震数据重建

采用多级中值滤波分离混采数据,可以很好地保持单炮的有效信息;对分离后的单炮数据,采用双曲Radon变换进行地震数据道重建,可以在Radon域获得极高的分辨率从而准确地重建地震数据。单炮道重建时,通过稀疏约束控制Radon变换来提高地震数据的重建精度。双曲Radon变换计算耗时较长,在算子求解时,本文采用快速迭代收缩阈值算法(FISTA),明显加快收敛速度,提高计算效率。模拟数据和实际数据表明,采用多级中值滤波和双曲Radon变换的方法重建混采地震数据可以获得较高的重建精度。