快速3D抛物Radon变换地震数据保幅重建

为解决3D AVO地震数据快速保幅重建问题,在传统3D抛物Radon变换的基础上提出一种3D快速高阶抛物Radon变换方法.该方法将传统抛物Radon变换与正交多项式相结合,通过正交多项式系数描述地震数据AVO信息,确保重建后的地震数据具有良好保幅效果.同时,该方法引入新变量λ_x=q_xf和λ_y=q_yf,通过对q_xf和q_yf的整体采样,消除了3D高阶抛物Radon变换算子对频率的依赖,使变换算子的求逆过程仅需计算一次,大大节省计算时间.理论模型和实际地震资料的处理结果表明,该方法重建效率高,保幅效果良好.     

3D高阶抛物Radon变换在不规则地震数据保幅重建中的应用

3D地震数据不规则采样缺失重建是地震勘探数据处理流程中的重要问题.本文提出了一种基于具有保幅特性的非均匀高阶抛物Radon变换（NHOPRT）地震数据重建方法.在最小二乘反演方程中引入Delaunay三角网格剖分来计算空间不规则加权系数,从而获得最接近完整规则数据的高阶抛物Radon变换域系数.在用SVD求解反演方程过程中,利用高阶抛物Radon变换算子在频率域为指数函数,具有线性可分解特性,将二维空间的高阶抛物Radon变换算子分解为两个独立的一维空间变换算子,减小了变换算子的矩阵大小,从而很大程度地提高了计算效率.理论模型和实际地震数据重建测试证明了本文方法的有效性以及实用性.     

3D高阶抛物Radon变换在不规则地震数据保幅重建中的应用

3D地震数据不规则采样缺失重建是地震勘探数据处理流程中的重要问题.本文提出了一种基于具有保幅特性的非均匀高阶抛物Radon变换(NHOPRT)地震数据重建方法.在最小二乘反演方程中引入Delaunay三角网格剖分来计算空间不规则加权系数,从而获得最接近完整规则数据的高阶抛物Radon变换域系数.在用SVD求解反演方程过程中,利用高阶抛物Radon变换算子在频率域为指数函数,具有线性可分解特性,将二维空间的高阶抛物Radon变换算子分解为两个独立的一维空间变换算子,减小了变换算子的矩阵大小,从而很大程度地提高了计算效率.理论模型和实际地震数据重建测试证明了本文方法的有效性以及实用性.     

基于3D加权高阶抛物Radon变换远偏移距地震数据保幅重建

在深部地震勘探中,远偏移距数据缺失重建是非常重要的工作.利用抛物Radon变换方法重建时,经过部分动校正后的数据在近偏移距和中偏移距的同相轴近似于抛物型,但是处于远偏移距位置的同相轴偏离于抛物型,特别是浅层同相轴,会导致远偏移距缺失的数据重建效果不佳.本文基于高阶抛物Radon变换地震数据保幅重建理论,重点分析了不同偏移距以及缺失程度对重建效果的影响,提出根据数据缺失模型来确定迭代过程中加权系数的3D高阶抛物Radon变换重建方法.同时,在反演求解目标函数的最小二乘解时应用SVD方法来求解Radon域系数,得到高精度的求解结果.该方法在大道间距数据内插、远偏移距缺失数据重建中都取得满意的效果.     

加权抛物Radon变换叠前地震数据重建

基于部分动校正（NMO）后反射同相轴在CMP道集上的抛物线走时近似,给出了加权抛物Radon变换叠前地震数据重建方法（WPRT）. WPRT通过在迭代过程中引入变化着的权系数,拓展和改进了传统抛物Radon变换方法,使其可同时完成不规则采样的规则化和空道及近偏移距道重建,且有更高的计算效率. 文中给出了应用WPRT进行近偏移距和中偏移距的空地震道重建及数据规则化的算法实现. 理论模型和实际地震资料的地震数据重建结果显示了本文算法的优点.     

反假频时延Radon变换多次波压制方法研究

多次波压制是地震数据处理的重要环节之一,而Radon变换方法以其优势被广泛应用于多次波处理环节中.本文提出利用时间延迟量Γ代替传统抛物Radon变换的曲率参数q,对地震记录同相轴进行积分运算.当使用时延Radon变换时,输入数据的道数作为重要参数直接影响Radon变换域的假频,而非传统抛物Radon变换空间方向的采样间隔.使用预条件算子通过重建间接增加输入数据的道数,消除Radon域的假频问题,提高成像精度,进而更好的实现多次波压制.利用本文方法对水平层状介质和复杂模型数据进行试算,结果表明,Radon域假频得到有效抑制,多次波压制效果明显.     

λ-f域加权抛物Radon变换地震数据重建方法研究

大多数地震处理技术都要求地震数据具备完整性和规则性,然而,由于诸多因素的影响,地震勘探所采集到的资料普遍存在数据缺失,需要对地震数据进行重建.本文在传统Radon变换重建的基础上提出一种λ-f域加权抛物Radon变换的地震数据重建方法.通过引入新变量λ,消除了Radon变换算子对频率的依赖,使得Radon变换算子及算子的逆只需计算一次,显著提高了计算效率.同时,在λ-f域抛物Radon变换迭代计算过程中引入变化的权系数,更好地实现了λ-f域的能量聚焦.理论模型及实际数据试算表明,文中方法对地震数据重建精度较高,单道对比吻合较好.     

抛物Radon变换法近偏移距波场外推

本文给出了抛物Radon变换的基本原理,以及部分动校正后的CMP道集抛物线近似有效性的证明,基于带限正反最小平方抛物Radon变换的Levinson递推算法,对缺失的近偏移距地震波场进行叠前重建和外推.给出了抛物Radon变换法地震道重建外推的基本原理和叠前地震数据规则化的处理流程,另外对于Radon域均匀采样的情形,本文给出了均匀层状介质和Marmousi模型的近偏移距外推结果,计算结果验证了算法的稳定性和适用性.     

基于SLA范数的高阶高分辨率λ-f域Radon变换多次波压制方法（英文）

传统的Radon变换是利用一次波与多次波的速度差异将其聚焦在Radon域内的不同"点"或"直线"上。然而有限的"偏移孔径"、"离散采样"和"地震数据的AVO特性"会使Radon变换的聚焦变差,降低了该方法的分辨率,影响多次波压制精度。此外传统Radon变换没有考虑地震数据的AVO特性,使用L_0范数的近似形式L1范数来提高Radon域的聚焦特性,计算量极大。本文就上述问题,本文将SL_0范数和正交多项式同时引入λ-f域Radon变换,提出基于SL_0范数的高阶高分辨率λ-f域Radon变换。在考虑地震数据的AVO特性基础上,将正交多项式和Radon变换相结合。首先推导出时间域高阶Radon正反变换形式,随后使用傅里叶变换将其变换到频率域,再引入变量,将频率和曲率相结合,对其进行整体采样,具备了较快的计算效率,最后使用平滑高斯函数替代L_1范数来近似L_0范数(SL_0范数),通过最速下降和梯度投影原理求得Radon变换在SL_0范数下的稀疏解。理论模型和实际资料试算表明,同基于L_1范数的Radon变换比较,本文方法既在Radon域内具有更好的聚焦特性,同时在有效压制多次波(NMO)后存在剩余时差)的同时能更好地保存一次波的AVO特性。     

混合域高分辨率抛物Radon变换及在衰减多次波中的应用

高分辨率Radon变换存在计算效率和分辨率不能兼得的困境.时间域算法可以获得很高的分辨率,但计算效率非常低;频率域算法具有良好计算效率,但分辨率不理想.为此发展了混合域高分辨率抛物Radon变换,即对频率域抛物Radon变换引入时变的稀疏权.本文给出了一种新的混合域高分辨率抛物Radon变换实现方法,并将该算法应用于叠前数据衰减多次波.文中给出了Radon变换和衰减多次波的流程.理论和实际数据算例表明本文方法既有较高的分辨率又有很高的计算效率.     

多方向正交多项式变换压制多次波

提出一种基于Radon 变换和正交多项式变换的多方向正交多项式变换压制多次波方法.抛物Radon变换对不同曲率方向的同相轴叠加,根据速度差异区分一次波和多次波,但Radon反变换会损伤振幅特性,不利于AVO分析.多方向正交多项式变换在Radon变换(某一曲率方向的零阶特性)的基础上,利用正交多项式变换进一步分析同相轴的高阶多项式分布特性,用正交多项式谱表征同相轴AVO特性;根据一次波和多次波速度差异和同相轴能量分布特征实现多次波压制.该方法的优点是仅用一个曲率参数就可描述同相轴剩余时差参数,提高了一次波和多次波的剩余时差分辨率.实验结果表明,该方法可以有效压制多次波并保留一次波AVO特性.     

基于保幅拉东变换的多次波衰减

为在去除多次波时有效保护地震一次反射波数据的AVO现象,给后续反演、解释提供准确的地震数据,本文提出了一种基于保幅拉东变换的多次波衰减方法,该方法是对常规抛物拉东变换的修改,把常规的稀疏拉东变换在拉东域分成两部分：一部分用于模拟零偏移距处的反射波能量,增加的另一部分用于模拟反射波振幅的AVO特性.该方法不仅考虑了反射波同相轴的形状,还考虑了反射波同相轴振幅幅度的变化,从而可把反射波信息进行有效转换,进而有利于多次波的消除,更好地恢复有效波的能量.在把地震数据由时间域转换到拉东域时,本文采用了IRLS算法实现保幅拉东算子的反演.模型数据和实际地震道集的试算分析表明,与常规拉东变换相比,保幅拉东变换在去除多次波的同时可有效保护一次反射波的AVO现象.     

基于稀疏反演三维表面多次波压制方法

三维表面多次波压制是海洋地震资料预处理中的重要研究课题,基于波动理论的三维表面多次波压制方法(3DSRME)是数据驱动的方法,理论上来说,可有效压制复杂构造地震数据表面多次波.但该方法因对原始地震数据采集要求高而很难在实际资料处理中广泛应用.本文基于贡献道集的概念,将稀疏反演方法引入到表面多次波压制中,应用稀疏反演代替横测线积分求和,无需对横测线进行大规模重建,进而完成三维表面多次波预测,这样可有效解决实际三维地震数据横测线方向稀疏的问题.基于纵测线多次波积分道集为抛物线的假设,为保证预测后三维表面多次波和全三维数据预测的多次波在运动学和动力学特征上基本一致,文中对预测数据实施基于稳相原理的相位校正.理论模型和实际数据的测试结果表明,本文基于稀疏反演三维表面多次波压制方法可在横测线稀疏的情况下,有效压制三维复杂介质地震资料中的表面多次波,从而更好地提高海洋地震资料的信噪比,为高分辨率地震成像提供可靠的预处理数据保障.     

Multiple attenuation using λ-f domain high-resolution Radon transform

The parabolic Radon transform has been widely used in multiple attenuation. To further improve the accuracy and efficiency of the Radon transform, we developed the 2- fdomain high-resolution Radon transform based on the fast and modified parabolic Radon transform presented by Abbad. The introduction of a new variable 2 makes the transform operator frequency-independent. Thus, we need to calculate the transform operator and its inverse operator only once, which greatly improves the computational efficiency. Besides, because the primaries and multiples are distributed on straight lines with different slopes in the 2-fdomain, we can easily choose the filtering operator to suppress the multiples. At the same time, the proposed method offers the advantage of high-resolution Radon transform, which can greatly improve the precision of attenuating the multiples. Numerical experiments suggest that the multiples are well suppressed and the amplitude versus offset characteristics of the primaries are well maintained. Real data processing results further verify the effectiveness and feasibility of the method.     

三维Radon变换的一种快速解析方法

3D Radon变换及其反变换是X-CT三维图像重建理论的核心,它在其他许多学科领域也有广泛应用。3D Radon变换的表达式是一个三重积分,按照定义直接计算相当费时。为此,研究一种新的快速的方法实现3D Radon变换,对X-CT图像重建理论及相关领域的发展有重要意义。本文以算法仿真常用椭球模型为基础,通过求解椭球模型与空间任意平面的面积,实现了用解析的方法快速得到模型的Radon变换,进一步比较了它与传统方法的优缺点,最后根据Radon反变换重建出原物体模型;计算机仿真结果验证了这种方法的正确。     

Time-domain hyperbolic Radon transform for separation of P-P and P-SV wavefields

A time-domain hyperbolic Radon transform based method for separating multicomponent seismic data into P-P and P-SV wavefields is presented. This wavefield separation method isolates P-P and P-SV wavefields in the Radon panel due to their differences in slowness, and an inverse transform of only part of the data leads to separated wavefields. A problem of hyperbolic Radon transform is that it works in the time domain entailing the inversion of large operators which is prohibitively time-consuming. By applying the conjugate gradient algorithm during the inversion of hyperbolic Radon transform, the computational cost can be kept reasonably low for practical application. Synthetic data examples prove that P-P and P-SV wavefield separation by hyperbolic Radon transform produces more accurate separated wavefields compared with separation by high-resolution parabolic Radon transform, and the feasibility of the proposed separation scheme is also verified by a real field data example.