首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
径向道变换压制相干噪声方法研究   总被引:6,自引:2,他引:4       下载免费PDF全文
径向道变换(Radial Trace Transform)是将地震道集振幅值从偏移距一双程旅行时坐标系变换到视速度-双程旅行时坐标系,通过这种坐标系的变换,使相干噪声与有效信号在视速度和频率方面都有效分离.本文在介绍RT变换基本原理基础上,分析了RT变换中两种常用插值方法及其特点.并利用对模拟地震资料的处理,证明了RT域模拟-相减法较其带通滤波法在相干噪声压制与反射信号保持方面具有明显优势.最后,根据噪声特点,通过选择合理RT滤波参数,对实际地震资料进行处理试验,获得了较好的去噪效果,明显提高了资料信噪比,验证了研究方法的有效性.  相似文献   

2.
基于非稳态多项式拟合的地震噪声衰减方法研究(英文)   总被引:1,自引:0,他引:1  
基于非稳态多项式拟合理论,针对地震数据中同相轴振幅变化这一特征,我们提出了一种地震噪声衰减的新方法。非稳态多项式拟合系数是时变的,通过整形正则化约束多项式拟和系数的光滑性,自适应的估计地震数据的相干分量。基于动校正后的共中心点道集(CMP)中地震信号的相干性,利用非稳态多项式拟合估计有效信号,从而衰减随机噪声。对于线性相干噪声,如地滚波,首先利用径向道变换(RadialTraceTransform,RTT)将地震数据变换到时间一视速度域,在时间—视速度域利用非稳态多项式拟合估计出相干噪声,然后减去相干噪声。该方法可以有效的估计振幅变化的相干分量,不需要相干分量振幅为常量的假设。模拟和实际资料处理结果表明,与传统的稳态多项式拟合和低切滤波相比,该方法可以更为有效的衰减地震噪声,同时保真了地震有效信号。  相似文献   

3.
Kirchhoff 3D prestack migration, as part of its execution, usually requires repeated access to a large traveltime table data base. Access to this data base implies either a memory intensive or I/O bounded solution to the storage problem. Proper compression of the traveltime table allows efficient 3D prestack migration without relying on the usually slow access to the computer hard drive. Such compression also allows for faster access to desirable parts of the traveltime table. Compression is applied to the traveltime field for each source location on the surface on a regular grid using 3D Chebyshev polynomial or cosine transforms of the traveltime field represented in the spherical coordinates or the Celerity domain. We obtain practical compression levels up to and exceeding 20 to 1. In fact, because of the smaller size traveltime table, we obtain exceptional traveltime extraction speed during migration that exceeds conventional methods. Additional features of the compression include better interpolation of traveltime tables and more stable estimates of amplitudes from traveltime curvatures. Further compression is achieved using bit encoding, by representing compression parameters values with fewer bits.  相似文献   

4.
叠前地震数据特征波场分解、偏移成像与层析反演   总被引:2,自引:2,他引:0       下载免费PDF全文
本文提出了一套叠前地震数据稀疏表达(特征波场合成)、深度偏移成像和层析成像的处理流程.不同于传统的变换域中的数据稀疏表达理论,本文利用局部平面波的传播方向(慢度矢量),在中心炮检点处同时进行波束合成,从而将地震数据投影到局部平面波域(高维空间)中.由于波束合成后的地震数据描述了局部平面波的方向特征,因此称之为特征波场.然而波束合成算法需要估计局部平面波的慢度矢量.当地震数据受噪声干扰时,难以在常规τ-p谱中自动估计局部平面波的射线参数(慢度矢量).本文提出了基于反演理论的特征波场合成方法,可以同时反演局部平面波及其传播方向,从而提高特征波合成的自动化程度并保持方法的稳健性.通过特征波场合成,可以将地震数据分解为单独的震相(波形).这样的数据可以直接用来成像及反演.在局部平面波域中,由于局部平面波的入射与出射射线参数已知,传统的Kirchhoff叠前深度偏移(PSDM)和高斯束/控制束PSDM可以实现从"沿等时面的画弧"到"向反射点(段)的直接投影"的转变,叠前偏移的效率以及成像质量可以同时提高.此外,特征波场与地下反射点(段)的一对一映射关系使得叠前深度偏移与层析成像融为一体,可以极大地提高速度反演的效率.数值试验证明了特征波场合成、叠前深度成像以及层析反演的有效性.  相似文献   

5.
多分辨率地震信号反褶积   总被引:11,自引:2,他引:9       下载免费PDF全文
基于二进小波变换提出了一种新的反褶积方法─-多分辨率地震信号反褶积.在地震信号二进小波变换域中的各尺度上分别进行其分辨率随小波尺度变化的反褶积,利用不同分辨率反褶积结果之间的相关性,以及测量噪声随尺度的衰减特性,从低分辨率反褶积结果逼近高分辨率反褶积结果.理论分析和实验表明,该方法有较高的精度,并且在较低信噪比情况下有好的效果.  相似文献   

6.
复杂介质下保真振幅Kirchhoff深度偏移   总被引:18,自引:2,他引:18       下载免费PDF全文
本文详细讨论地下复杂介质地震成像问题.当速度模型含有强横向变化时,射线场会折叠形成多次走时,在这种情况下,文中证明在保真振幅叠前深度偏移时,应当考虑所有的到时才能得到定量的像.通过实例证明,在多次走时情况下,常规的共炮集、共检波器集都存在强的假像;由此提出共角度偏移的概念.最后,在多次走时下,讨论多次偏移算子(走时、相移、振幅等)的高效重构,给出了三维叠前深度偏移的快速算法.  相似文献   

7.
基于Dreamlet变换的地震数据压缩理论与方法   总被引:2,自引:1,他引:1       下载免费PDF全文
为达到更加有效地表示地震数据的目的,仅仅将地震数据当作普通的图像数据处理是远远不够的,地震数据中蕴含的地震波的运动学特性也应作为重要因素而被考虑到.本文讨论了利用Dreamlet变换方法实现地震数据压缩的方法,并针对地震数据本身所蕴含的频散关系特性进一步提出了多尺度Dreamlet变换压缩方法.Dreamlet变换由2个一维局部谐波变换的张量积构成,它在提供地震波场时间-空间局部化性质的同时可以保留波场的运动学特性.通过对二维SEG/EAGE叠前、叠后数据的算例说明了Dreamlet变换用于地震数据压缩的有效性.利用压缩后的数据进行成像的结果更表明,与Curvelet变换方法相比,Dreamlet与多尺度Dreamlet方法可以提供更高的压缩比;在相同压缩比的条件下,使用Dreamlet与多尺度Dreamlet方法压缩重建后的数据进行成像能更好地保留成像结果中的重要结构.  相似文献   

8.
In log time-frequency spectra, the nonstationary convolution model is a linear equation and thus we improved the Gabor deconvolution by employing a log hyperbolic smoothing scheme which can be implemented as an iteration process. Numerical tests and practical applications demonstrate that improved Gabor deconvolution can further broaden frequency bandwidth with less computational expenses than the ordinary method. Moreover, we attempt to enlarge this method’s application value by addressing nonstationary and evaluating Q values. In fact, energy relationship of each hyperbolic bin (i.e., attenuation curve) can be taken as a quantitative indicator in balancing nonstationarity and conditioning seismic traces to the assumption of unchanging wavelet, which resultantly reveals more useful information for constrained reflectivity inversion. Meanwhile, a statistical method on Q-value estimation is also proposed by utilizing this linear model’s gradient. In practice, not only estimations well agree with geologic settings, but also applications on Q-compensation migration are favorable in characterizing deep geologic structures, such as the pinch-out boundary and water channel.  相似文献   

9.
2D共炮时间域高斯波束偏移   总被引:1,自引:1,他引:0       下载免费PDF全文
针对传统射线方法在奇异区成像精度不高,而2D频率域高斯波束叠前深度偏移需要计算成像点处每个频率的格林函数,影响计算效率的问题,本文通过使用复走时代替实走时,改变频率域下成像公式的积分顺序,给出了在时间域下进行高斯波束偏移的方法和计算公式.本文使用复杂数值模型验证了2D时间域高斯波束叠前偏移方法的正确性,并同传统射线偏移成像结果做了对比.对比结果表明时间域高斯波束偏移在成像精度上优于传统射线偏移.  相似文献   

10.
冯波  王华忠  冯伟 《地球物理学报》2019,62(4):1471-1479
地震波的运动学信息(走时、斜率等)通常用于宏观速度建模.针对走时反演方法,一个基本问题是走时拾取或反射时差的估计.对于成像域反演方法,可以通过成像道集的剩余深度差近似计算反射波时差.在数据域中,反射地震观测数据是有限频带信号,如果不能准确地确定子波的起跳时间,难以精确地确定反射波的到达时间.另一方面,如果缺乏关于模型的先验信息,则很难精确测量自地下同一个反射界面的观测数据同相轴和模拟数据同相轴之间的时差.针对走时定义及时差测量问题,首先从叠前地震数据的稀疏表达出发,利用特征波场分解方法,提取反射子波并估计局部平面波的入射和出射射线参数.进一步,为了实现自动和稳定的走时拾取,用震相的包络极值对应的时间定义反射波的到达时,实现了立体数据中间的自动生成.理论上讲,利用包络极值定义的走时大于真实的反射波走时,除非观测信号具有无限带宽(即delta脉冲).然而,走时反演的目的是估计中-大尺度的背景速度结构,因此走时误差导致的速度误差仍然在可以接受的误差范围内.利用局部化传播算子及特征波聚焦成像条件将特征波数据直接投影到地下虚拟反射点,提出了一种新的反射时差估计方法.既避免了周期跳跃现象以及串层等可能性,又消除了振幅因素对时差测量的影响.最后,在上述工作基础之上,提出了一种基于特征波场分解的新型全自动反射走时反演方法(CWRTI).通过对泛函梯度的线性化近似,并用全变差正则化方法提取梯度的低波数部分,实现了背景速度迭代反演.在理论上,无需长偏移距观测数据或低频信息、对初始模型依赖性低且计算效率高,可以为后续的全波形反演提供可靠的初始速度模型.理论和实际资料的测试结果证明了本文方法的有效性.  相似文献   

11.
地层品质因子Q的可用于地震资料高分辨率处理,而从VSP资料下行直达波更容易获取准确的地层品质因子。通过对零偏移距VSP资料的监控子波和下行初至波的频谱进行综合分析,仿照Ricker子波频谱的表达式,本文提出了震源子波频谱新的表达式。在震源子波频谱新的表达式基础上,我们介绍了改进的频谱拟合法和改进的谱比法的层Q值反演方法及相应的处理流程。基于本文提出的层Q值反演方法,利用实际的零偏移距VSP资料的下行直达波,反演稳定的层Q值,并用于零偏移距VSP资料及井旁地面地震资料的反Q滤波振幅补偿处理,提高了地震资料的分辨率。  相似文献   

12.
在地震弹性矢量波场框架下,推导了多波联合层析速度反演方程以及走时残差与角道集剩余曲率的转换关系式,提出了一种利用成像域角道集更新P波、S波速度的走时层析反演方法.其实现过程可以概括为:将弹性波多分量数据作为输入,基于高斯束实现矢量波场成像并提取角道集,利用层析反演方程求解慢度更新量,最终获得多波联合反演结果.模型试算和实际资料处理验证了该方法的反演效果,能够为弹性矢量波联合深度偏移提供高质量的叠前速度场.  相似文献   

13.
The resolution of seismic data is critical to seismic data processing and the subsequent interpretation of fine structures. In conventional resolution improvement methods, the seismic data is assumed stationary and the noise level not changes with space, whereas the actual situation does not satisfy this assumption, so that results after resolution improvement processing is not up to the expected effect. To solve these problems, we propose a seismic resolution improvement method based on the secondary time–frequency spectrum. First, we propose the secondary time-frequency spectrum based on S transform (ST) and discuss the reflection coefficient sequence and time-dependent wavelet in the secondary time–frequency spectrum. Second, using the secondary time–frequency spectrum, we design a twodimensional filter to extract the amplitude spectrum of the time-dependent wavelet. Then, we discuss the improvement of the resolution operator in noisy environments and propose a novel approach for determining the broad frequency range of the resolution operator in the time–frequency–space domain. Finally, we apply the proposed method to synthetic and real data and compare the results of the traditional spectrum-modeling deconvolution and Q compensation method. The results suggest that the proposed method does not need to estimate the Q value and the resolution is not limited by the bandwidth of the source. Thus, the resolution of the seismic data is improved sufficiently based on the signal-to-noise ratio (SNR).  相似文献   

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.
目前对多次波的有效利用仅围绕多次波成像技术展开,通过成像多次波试图获取更丰富的地下构造信息.不同于该思路,本文另辟蹊径,从利用多次波提高地震数据分辨率角度出发,对多次波的有效利用进行了深入挖掘.首先基于聚焦变换思想在聚焦域内实现多次波的降阶,通过理论推导得出聚焦域内多次波表现为原始数据的多维子波反褶积这一重要结论,从理论上证明了本文方法提高地震数据分辨率的可行性;然后采用引入整形正则化的非稳态回归自适应匹配滤波方法将聚焦域内由多次波构建的高分辨率数据分离出来,实现原始数据的高分辨率转换.与常规反褶积模型不同,该方法基于波动理论推导得出,可以适用于任意复杂情况;每一道输出结果中所有炮记录都参与了运算,从空间上加以约束,在提高纵向分辨率的同时可以改善数据的横向分辨率.最后通过模型试算和实际资料处理对本文方法的有效性、适应性和实用性进行了验证.  相似文献   

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

17.
基于单程波深度延拓算法,得到震源波场在成像点的入射角度,结合对地层倾角的估计,获得入射地震波与界面法线的夹角.通过运用"保幅"的反褶积成像条件和考虑累加的炮数,解决了炮点覆盖不均匀导致的成像幅值误差问题,进而建议了炮域波动方程叠前深度偏移直接产生角道集的方法和流程.与基于空间移动或时移成像条件的波动方程叠前深度偏移提取角道集的方法相比,本文建议的方法只需少量额外的存储空间,又可补偿观测系统非均匀覆盖对成像幅值的影响;其增加的计算量与炮域偏移算法相比几乎可以忽略.文中算例表明,本文方法提取的角道集可为叠前反演提供较精确的AVO振幅特性.此外,就改善地震成像效果本身而言,提取角道集使得可在波动方程叠前深度偏移中应用剩余动校和拉伸切除技术,从而可更好地保持高频成分并提高成像的信噪比.  相似文献   

18.
Radon transform is a powerful tool with many applications in different stages of seismic data processing, because of its capability to focus seismic events in the transform domain. Three-parameter Radon transform can optimally focus and separate different seismic events, if its basis functions accurately match the events. In anisotropic media, the conventional hyperbolic or shifted hyperbolic basis functions lose their accuracy and cannot preserve data fidelity, especially at large offsets. To address this issue, we propose an accurate traveltime approximation for transversely isotropic media with vertical symmetry axis, and derive two versions of Radon basis functions, time-variant and time-invariant. A time-variant basis function can be used in time domain Radon transform algorithms while a time-invariant version can be used in, generally more efficient, frequency domain algorithms. Comparing the time-variant and time-invariant Radon transform by the proposed basis functions, the time-invariant version can better focus different seismic events; it is also more accurate, especially in presence of vertical heterogeneity. However, the proposed time-invariant basis functions are suitable for a specific type of layered anisotropic media, known as factorized media. We test the proposed methods and illustrate successful applications of them for trace interpolation and coherent noise attenuation.  相似文献   

19.
张雅晨  刘洋  刘财  武尚 《地球物理学报》2019,62(3):1181-1192
地震数据本质上是时变的,不仅有效同相轴表现出确定性信号的时变特征,而且复杂地表和构造条件以及深部探测环境总是引入时变的非平稳随机噪声.标准的频率-空间域预测滤波只适合压制平面波信号假设下的平稳随机噪声,而处理非平稳地震随机噪声时,需要将数据体分割为小窗口进行分析,但效果不够理想,而传统非预测类随机噪声压制方法往往适应性不高,因此开发能够保护地震信号时变特征的随机噪声压制方法具有重要的工业价值.压缩感知是近年出现的一个新的采样理论,通过开发信号的稀疏特性,已经在地震数据处理中的数据插值以及噪声压制中得到了应用.本文系统地分析了压缩感知理论框架下的地震随机噪声压制问题,建立了阈值消噪的数学反演目标函数;针对时变有效信息具有的可压缩性,利用有限差分算法求解炮检距连续方程,构建有限差分炮检距连续预测算子(FDOC),在seislet变换框架下,提出一种新的快速稀疏变换域———FDOC-seislet变换,实现地震数据的高度稀疏表征;结合非平稳随机噪声不可压缩的特征,提出了一种整形迭代消噪方法,该方法是一种广义的迭代收缩阈值(IST)算法,在无法计算稀疏变换伴随算子的条件下,仍然能够对强噪声环境中的时变有效信息进行有效恢复.通过对模型数据和实际数据的处理,验证了FDOC-seislet稀疏变换域随机噪声迭代压制方法能够在保护复杂构造地震波信息的前提下,有效地衰减原始数据中的强振幅随机噪声干扰.  相似文献   

20.
非稳态地震稀疏约束反褶积研究(英文)   总被引:1,自引:1,他引:0  
传统Robinson褶积模型主要受缚于三种不合理的假设,即白噪反射系数、最小相位地震子波与稳态假设,而现代反射系数反演方法(如稀疏约束反褶积等)均在前两个假设上寻求突破的同时却忽视了一个重要事实:实际地震信号具有典型的非稳态特征,这直接冲击着反射系数反演中地震子波不随时间变化的这一基础性假设。本文首先通过实际反射系数测试证实,非稳态效应造成重要信息无法得到有效展现,且对深层影响尤为严重。为校正非稳态影响,本文从描述非稳态方面具有普适性的非稳态褶积模型出发,借助对数域的衰减曲线指导检测非稳态影响并以此实现对非稳态均衡与校正。与常规不同,本文利用对数域Gabor反褶积仅移除非稳态影响,而将分离震源子波和反射系数的任务交给具有更符合实际条件的稀疏约束反褶积处理,因此结合两种反褶积技术即可有效解决非稳态特征影响,又能避免反射系数和地震子波理想化假设的不利影响。海上地震资料的应用实际表明,校正非稳态影响有助于恢复更丰富的反射系数信息,使得与地质沉积和构造相关的细节特征得到更加清晰的展现。  相似文献   

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

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