首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 22 毫秒
1.
加权抛物Radon变换叠前地震数据重建   总被引:10,自引:6,他引:10       下载免费PDF全文
基于部分动校正(NMO)后反射同相轴在CMP道集上的抛物线走时近似,给出了加权抛物Radon变换叠前地震数据重建方法(WPRT). WPRT通过在迭代过程中引入变化着的权系数,拓展和改进了传统抛物Radon变换方法,使其可同时完成不规则采样的规则化和空道及近偏移距道重建,且有更高的计算效率. 文中给出了应用WPRT进行近偏移距和中偏移距的空地震道重建及数据规则化的算法实现. 理论模型和实际地震资料的地震数据重建结果显示了本文算法的优点.  相似文献   

2.
The delay‐time Radon transform parametrizes coherent events in a seismic gather by the far‐offset trace delay time, instead of the conventional parabolic curvature or ray parameter. The reformulation may give a different physical insight into the aliasing effect in the Radon transformation and may also lead to a different algorithm. The delay‐time parametrization enables modelling of a seismic gather as the sum of coherent events with any form of moveout curve. For example, a parabolic curve can be used for traces within a moderate offset range and a linear moveout for far‐offset traces. When using this delay‐time Radon transform, it is the number of traces, rather than the spatial sampling, of the input gather that directly controls aliasing in the Radon transform image. A preconditioning operator that implicitly increases the number of input traces by spatial reconstruction (without physically performing the spatial resampling) may minimize aliasing noise in the Radon transform image.  相似文献   

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

4.
Radon变换在非规则观测系统中波场分离的应用   总被引:2,自引:1,他引:1       下载免费PDF全文
Radon变换是数据处理中广泛应用的一种方法技术.本文介绍了Radon变换在规则及其非规则观测系统中波场分离中应用,讨论了Radon变换的公式、实现方法和有关计算参数的选择.通过理论模型的试算,阐明了Radon变换在缺失道、道距变化情况下依然可以高精度的分离波场,并对线性、抛物Radon变换进行了比较,对隧道超前预报反射波资料进行了处理,给出了隧道超前预报处理流程中提取反射波的算例.  相似文献   

5.
逆散射级数法去除自由表面多次波   总被引:8,自引:5,他引:3       下载免费PDF全文
李翔  胡天跃 《地球物理学报》2009,52(6):1633-1640
压制多次波是地震资料处理过程中的重要环节之一,多次波去除的效果直接影响到最后成像的质量.本文研究的逆散射级数法以波动方程和Born级数为基础,构造一个与自由表面多次波相关的子序列来预测并消除多次波,该方法不需要知道地下的速度结构信息,适用于各种复杂的地下情况.在逆散射级数法的实现过程中,地震子波和近道波场是必不可少的信息,本文采用能量最低法则和带限抛物Radon变换法来弥补原始地震资料中所缺失的这些信息.合成资料和实际地震资料的处理结果表明,逆散射级数法能够在去除自由表面多次波的同时保留有效波的信息,对于地下复杂结构的情况同样有效.  相似文献   

6.
多项式Radon变换   总被引:21,自引:12,他引:21       下载免费PDF全文
Radon变换是数据处理中广泛应用的一种方法技术.本文介绍了次数为“2”的多项式Radon变换.讨论了多项式正反Radon变换的公式、实现方法和有关计算参数的选择.通过理论模型试算,对多项式、线性、抛物线三种Radon变换进行了比较.对实际资料进行了多项式Radon变换处理,给出了消除地震反射记录中线性干扰的算例.  相似文献   

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

8.
即使采用分辨率很高的双曲Radon变换,对速度各向异性发育介质及长偏移距情况下的地震数据,其Radon域内能量仍不收敛.为了克服此难题,我们在Radon变换的积分路径中考虑了非双曲走时的影响,通过引入非双曲时差公式中的各向异性非椭圆率η参数,可以准确描述出长偏移距条件下来自同一层位的时距曲线,并推导了由偏移距、慢度、非椭圆率三参数控制的积分曲线正反变换公式,我们称之为各向异性Radon变换.离散化求解时,各向异性Radon变换是时变的,频率域快速算法已不适用,本文采用了最优相似系数加权Gauss-Seidel迭代算法,保持其计算精度的同时也有较高的计算效率.将此方法应用在模型数据以及实际长偏移距海上地震数据的多次波压制处理中,收到了较好的处理效果.  相似文献   

9.
多次波偏移中的假象主要来自于不同地震事件之间的互相关,由于这种互相关满足成像条件,很难直接在偏移过程中去除.但是对于准确的速度模型,真实的成像结果在角度域内应该是平直的.根据这个判断准则,可以在角度域内移除多次波偏移中的假象.本文以数据自相关偏移为例,提出了在单程波多次波偏移中移除假象的主要流程:首先在在单程波偏移过程中高效地提取角度域共成像点道集,然后对角度域共成像点道集应用高分辨率的抛物线型Radon变换,用合适的切除函数处理后,反变换回到角度域,最后叠加各个角度成分,得到偏移结果.Marmousi模型的合成数据测试表明,这种方法可以很好地压制多次波偏移过程中产生的假象,有效地提高成像结果的信噪比.  相似文献   

10.
A conventional velocity-stack gather consists of constant-velocity CMP-stacked traces. It emphasizes the energy associated with the events that follow hyperbolic traveltime trajectories in the CMP gather. Amplitudes along a hyperbola on a CMP gather ideally map onto a point on a velocity-stack gather. Because a CMP gather only includes a cable-length portion of a hyperbolic traveltime trajectory, this mapping is not exact. The finite cable length, discrete sampling along the offset axis and the closeness of hyperbolic summation paths at near-offsets cause smearing of the stacked amplitudes along the velocity axis. Unless this smearing is removed, inverse mapping from velocity space (the plane of stacking velocity versus two-way zero-offset time) back to offset space (the plane of offset versus two-way traveltime) does not reproduce the amplitudes in the original CMP gather. The gather resulting from the inverse mapping can be considered as the model CMP gather that contains only the hyperbolic events from the actual CMP gather. A least-squares minimization of the energy contained in the difference between the actual CMP gather and the model CMP gather removes smearing of amplitudes on the velocity-stack gather and increases velocity resolution. A practical application of this procedure is in separation of multiples from primaries. A method is described to obtain proper velocity-stack gathers with reduced amplitude smearing. The method involves a t2-stretching in the offset space. This stretching maps reflection amplitudes along hyperbolic moveout curves to those along parabolic moveout curves. The CMP gather is Fourier transformed along the stretched axis. Each Fourier component is then used in the least-squares minimization to compute the corresponding Fourier component of the proper velocity-stack gather. Finally, inverse transforming and undoing the stretching yield the proper velocity-stack gather, which can then be inverse mapped back to the offset space. During this inverse mapping, multiples, primaries or all of the hyperbolic events can be modelled. An application of velocity-stack processing to multiple suppression is demonstrated with a field data example.  相似文献   

11.
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.  相似文献   

12.
3D高阶抛物Radon变换地震数据保幅重建   总被引:1,自引:2,他引:1       下载免费PDF全文
本文结合传统3D抛物Radon变换(PRT)和AVO数据正交多项式拟合,给出了3D高阶抛物Radon变换方法(HOPRT).该变换增加了描述AVO数据变化的梯度信息和曲率信息,拓展了传统3D抛物Radon变换方法,使其在具有AVO特征的数据重建中具有更高的准确度,从而提高AVO分析的可靠性.文中给出了3D高阶抛物Radon变换进行地震数据保幅重建的流程.理论模型和实际地震资料的重建结果显示了本文方法的优点.  相似文献   

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

14.
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.  相似文献   

15.
动校正拉伸是地震资料处理的一个基本问题,解决拉伸问题的处理方法是切除.现代地震数据大多为长排列采集,动校正拉伸更为严重.依据褶积模型和Fourier变换的基本性质,本文给出频谱代换无拉伸动校正方法.算法实现就是将CMP道集变换到频率域,取参考道的相位谱替换其它偏移距道的相位,同时保持其振幅谱不变,再做Fourier反变换就得到动校正后的地震剖面.通过其实现过程可知该方法不需要地下介质的速度信息,算法可完全自动实现,且具有较高的计算效率.频谱代换无拉伸动校正可适用于任何偏移距的地震资料,而且还可有效保持地震资料的AVO效应.理论模拟数据及其叠加结果显示频谱代换法的有效性和实用性,同时该方法具有较强的抗随机噪音能力.  相似文献   

16.
A new seismic interpolation and denoising method with a curvelet transform matching filter, employing the fast iterative shrinkage thresholding algorithm (FISTA), is proposed. The approach treats the matching filter, seismic interpolation, and denoising all as the same inverse problem using an inversion iteration algorithm. The curvelet transform has a high sparseness and is useful for separating signal from noise, meaning that it can accurately solve the matching problem using FISTA. When applying the new method to a synthetic noisy data sets and a data sets with missing traces, the optimum matching result is obtained, noise is greatly suppressed, missing seismic data are filled by interpolation, and the waveform is highly consistent. We then verified the method by applying it to real data, yielding satisfactory results. The results show that the method can reconstruct missing traces in the case of low SNR (signal-to-noise ratio). The above three problems can be simultaneously solved via FISTA algorithm, and it will not only increase the processing efficiency but also improve SNR of the seismic data.  相似文献   

17.
熊登  赵伟  张剑锋 《地球物理学报》2009,52(4):1068-1077
高分辨率Radon变换存在计算效率和分辨率不能兼得的困境.时间域算法可以获得很高的分辨率,但计算效率非常低;频率域算法具有良好计算效率,但分辨率不理想.为此发展了混合域高分辨率抛物Radon变换,即对频率域抛物Radon变换引入时变的稀疏权.本文给出了一种新的混合域高分辨率抛物Radon变换实现方法,并将该算法应用于叠前数据衰减多次波.文中给出了Radon变换和衰减多次波的流程.理论和实际数据算例表明本文方法既有较高的分辨率又有很高的计算效率.  相似文献   

18.
在野外数据采集过程中,空间非均匀采样下的地震道缺失现象经常出现,为了不影响后续资料处理,必须进行高精度数据重建.然而大多数常规方法只能对空间均匀采样下的地震缺失道进行重建,而对于非均匀采样的地震数据则无能为力.为此本文在以往多尺度多方向二维曲波变换的基础上,首先引入非均匀快速傅里叶变换,建立均匀曲波系数与空间非均匀采样下地震缺失道数据之间的规则化反演算子,在L1最小范数约束下,使用线性Bregman方法进行反演计算得到均匀曲波系数,最后再进行均匀快速离散曲波反变换,从而形成基于非均匀曲波变换的高精度地震数据重建方法.该方法不仅可以重建非均匀带假频的缺失数据,而且具有较强的抗噪声能力,同时也可以将非均匀网格数据归为到任意指定的均匀采样网格.理论与实际数据的处理表明了该方法重建效果远优于非均匀傅里叶变换方法,可以有效地指导复杂地区数据采集设计及重建.  相似文献   

19.
Receiver ghost reflections adversely affect variable-depth streamer (VDS) data acquisition. In addition, the frequency notches caused by the interference between receiver ghosts and primary waves strongly affect seismic data processing and imaging. We developed a high-resolution Radon transform algorithm and used it to predict receiver ghosts from VDS data. The receiver ghost reflections are subtracted and removed from the raw data. We propose a forward Radon transform operator of VDS data in the frequency domain and, based on the ray paths of the receiver ghosts, we propose an inverse Radon transform operator. We apply the proposed methodology to model and field data with good results. We use matching and subtracting modules of commercially available seismic data processing software to remove the receiver ghosts. The frequency notches are compensated and the effective frequency bandwidth of the seismic data broadens.  相似文献   

20.
多次波的存在会降低地震资料的信噪比,影响地震资料处理效果和后续的地质解释精度,压制多次波干扰是地震资料处理中的一个重要环节。利用有限差分方法正演模拟几个典型地质模型的地震波场,之后进行动校正,再通过抛物线Radon变换方法把t-x域一次反射波和多次波变换到τ-q域,由于存在速度的差异,Radon域记录中的一次反射波和多次波的能量互相分开,在Radon域对多次波的能量切除,保留反射波的有效信息,达到压制多次波的效果,提高资料的信噪比。通过对实际地震数据进行抛物线Radon变换结果也进一步验证了该方法在压制多次波中的应用效果。   相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号