吉林大学学报(地球科学版) ›› 2019, Vol. 49 ›› Issue (6): 1780-1787.doi: 10.13278/j.cnki.jjuese.20180299

• 地球探测与信息技术 • 上一篇    下一篇

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

郑植升1, 刘洋1, 刘财1, 张亮2   

  1. 1. 吉林大学地球探测科学与技术学院, 长春 130026;
    2. 中国石油吉林油田公司乾安采油厂, 吉林 松原 138000
  • 收稿日期:2018-11-07 发布日期:2019-11-30
  • 通讯作者: 刘洋(1979-),男,教授,博士生导师,主要从事地震数据非平稳模式分析、地震数据和海洋电磁数据处理方法研究,以及地质-地球物理综合研究工作,E-mail:yangliu1979@jlu.edu.cn E-mail:yangliu1979@jlu.edu.cn
  • 作者简介:郑植升(1993-),男,硕士研究生,主要从事地震数据处理方法研究,E-mail:zhengzs17@mails.jlu.edu.cn
  • 基金资助:
    国家自然科学基金项目(41774127,41522404);吉林大学高层次科技创新团队建设项目(2017TD-14)

Two Implementations for τ-p Transforms

Zheng Zhisheng1, Liu Yang1, Liu Cai1, Zhang Liang2   

  1. 1. College of GeoExploration Science and Technology, Jilin University, Changchun 130026, China;
    2. Qian'an Oil Factory, Jilin Oilfield, PetroChina, Songyuan 138000, Jilin, China
  • Received:2018-11-07 Published:2019-11-30
  • Supported by:
    Supported by National Natural Science Foundation of China (41774127,41522404) and High-Level Science and Technology Innovation Team Building Project of Jilin University (2017TD-14)

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

关键词: τ-p变换, 投影切片定理, 径向道变换, 斜率分解

Abstract: τ-p transform is a classic projection transform method, which is widely used in seismic data processing such as signal-noise separation, seismic trace interpolation, and plane wave decomposition based on different ray parameters. However, the precision and resolution of τ-p transform are limited by mathematical inverse problems, and it is necessary to study the methods to accelerate the calculation speed of τ-p transform or to design relevant filter, while maintaining high precision. We propose two τ-p transform methods based on radial trace transform and slope decomposition respectively. In this paper,we demonstrate the fundamental theory of the two methods and compare them with conventional τ-p transform. Using fast Fourier transform and radial trace transform, the τ-p transform method can efficiently reduce the time consumption. Utilizing stable slope decomposition algorithm, the τ-p transform method can perform high resolution τ-p transform,and provide redundant transformation dimension (τ-x-p domain),which enable the design of filter operator more flexible. The numerical experiments show that these two proposed methods are superior to the conventional τ-p transform in calculation speed and reconstruction accuracy. The field-data examples for ground-roll noise suppression confirm that the two methods can provide more effective and flexible implementations for seismic data processing.

Key words: τ-p transform, projection-slice theorem, radial trace transform, slope decomposition

中图分类号: 

  • P631.4
[1] Stoffa P L, Buhl P, Diebold J B, et al. Direct Mapping of Seismic Data to the Domain of Intercept Time and Ray Parameter:A Plane-Wave Decomposition[J]. Geophysics, 1981, 46(3):255-267.
[2] Durrani T S, Bisset D. The Radon Transform and Its Properties[J]. Geophysics, 1984, 49(8):1180-1187.
[3] Turner G. Aliasing in the Tau-p Transform and the Removal of Spatially Aliased Coherent Noise[J]. Geophysics, 1990, 55(11):1496-503.
[4] Yilmaz Öz. Seismic Data Analysis:Processing, Inversion, and Interpretation of Seismic Data[M]. Katras:Society of Exploration Geophysicists, 2001.
[5] Wilson C, Guitton A. Teleseismic Wavefield Interpolation and Signal Extraction Using High-Resolution Linear Radon Transforms[J]. Geophysical Journal International, 2007, 168(1):171-181.
[6] Levin S A. Frequency-Dip Formulation of Wave-Theoretic Migration in Stratified Media:Acoustical Imaging[M]. Boston:Springer, 1980.
[7] Forbriger T. Inversion of Shallow-Seismic Wavefields:I:Wavefield Transformation[J]. Geophysical Journal International, 2003, 153(3):719-734.
[8] Dev A, McMechan G A. Spatial Antialias Filtering in the Slowness-Frequency Domain[J]. Geophysics, 2009, 74(2):V35-V42.
[9] McMechan G A. Px Imaging by Localized Slant Stacks of Tx Data[J]. Geophysical Journal International, 1983, 72(1):213-221.
[10] Tang R, Carswell A, Moon W. Velocity Analysis in the p-x Plane from a Slant Stack Wavefield[J]. Geophysical Prospecting, 1984, 32(6):1016-1032.
[11] Ottolini R. Signal/Noise Separation in Dip Space[J]. SEP Report, 1983, 37:143-149.
[12] Ghosh S, Fomel S. Multiple Suppression in the t-x-p Domain[C]//SEG Technical Program Expanded Abstracts 2012. Las Vegas:Society of Exploration Geophysicists, 2012:1-6.
[13] Radon J. Vber die Bestimmung von Funktionen durch ihre Integralwerte längs gewisser Mannigfaltigkeiten, Ber Verh Sächs Akad Wiss Leipzig[J]. Math-Nat kl, 1917, 69:262-277.
[14] Claerbout J F, Johnson A G. Extrapolation of Time-Dependent Waveforms Along Their Path of Propagation[J]. Geophysical Journal International, 1971, 26(1):285-293.
[15] Claerbout J F. Earth Soundings Analysis:Processing Versus Inversion[M]. Boston:Blackwell Scientific Publications, 1992.
[16] 刘清林. Tau-p变换的频率域算法及Tau-p域的偏移处理[J]. 长春地质学院学报, 1986, 16(4):77-80. Liu Qinglin.The Calculating Method of Tau-p Transform in F-K Domain and Migration Method in Tau-p Domain[J]. Journal of Changchun University of Earth Science, 1986, 16(4):77-80.
[17] Averbuch A, Coifman R R, Donoho D L, et al. Fast Slant Stack:A Notion of Radon Transform for Data in a Cartesian Grid Which Is Rapidly Computible, Algebraically Exact, Geometrically Faithful and Invertible[J]. SIAM Scientific Computing, 2001, 37(3):192-206.
[18] Claerbout J F. Slant-Stacks and Radial Traces[R]. Stanford:Stanford Exploration Project, 1975:1-12.
[19] Fomel S. Inverse B-Spline Interpolation[R]. Stanford:Stanford Exploration Project, 2000.
[20] Liu Y, Fomel S. Seismic Data Analysis Using Local Time-Frequency Decomposition[J]. Geophysical Prospecting, 2012, 61(3):516-525.
[21] 张雪冰,刘财,刘洋,等.基于局部均值分解的地震信号时频分解方法[J].吉林大学学报(地球科学版), 2017, 47(5):1562-1571. Zhang Xuebing, Liu Cai, Liu Yang, et al. Seismic Data Time-Frequency Decomposition Based on Local Mean Decomposition[J]. Journal of Jilin University (Earth Science Edition), 2017, 47(5):1562-1571.
[22] Fomel S. Shaping Regularization in Geophysical-Estimation Problems[J]. Geophysics, 2007, 72(2):R29-R36.
[1] 丁梦颜, 冯晅, 刘财. 裂缝对岩石非线性弹性特征的影响[J]. 吉林大学学报(地球科学版), 2026, 56(3): 1051-1061.
[2] 范佳奇, 吴秋莹, 王典, 李鹏.

 Radon变换的贪婪-快速迭代收缩阈值算法实现及多次波压制应用 [J]. 吉林大学学报(地球科学版), 2026, 56(2): 684-693.

[3] 赵天硕, 宋超, , 刘财, , , , 徐雨歆.

物理信息神经网络地震走时层析成像程函方程因式分解方法 [J]. 吉林大学学报(地球科学版), 2026, 56(2): 694-702.

[4] 谭晓淼, 周建波, 饶莹, 王海燕, 侯贺晟, 李明芮, 高锐. 中亚造山带中南部索伦缝合带岩石圈结构及其对古亚洲洋演化的制约[J]. 吉林大学学报(地球科学版), 2026, 56(1): 209-218.
[5] 毛子雄, 侯贺晟, 周建波, 符伟. 佳木斯地块—那丹哈达地体近地表速度结构[J]. 吉林大学学报(地球科学版), 2026, 56(1): 219-228.
[6] 韩复兴, 王源, 高正辉, 常志邈, 马飞, 秦昊, 尚浩. 基于NEWUOA的CRS叠加成像技术[J]. 吉林大学学报(地球科学版), 2026, 56(1): 386-396.
[7] 黄兴国, 翁央央, 韩丽. 基于广义递归卷积的孔隙黏弹地震波正演模拟[J]. 吉林大学学报(地球科学版), 2026, 56(1): 377-385.
[8] 韩复兴, 刘水源, 高正辉, 韩江涛, 张涛, 尚浩. 基于机器学习的微动HVSR数据干扰信号压制方法[J]. 吉林大学学报(地球科学版), 2025, 55(6): 2153-2163.
[9] 刘财, 张焱喆, 刘洋. 地震数据非平稳特征分析综述[J]. 吉林大学学报(地球科学版), 2025, 55(6): 2132-2152.
[10] 齐娇, 曹思远. 基于先验知识的深度学习表面相关多次波压制方法[J]. 吉林大学学报(地球科学版), 2025, 55(5): 1702-1714.
[11] 李丛, 张栋, 董果果, 袁青松, 许军, 朱德胜, 代磊, 李鹏飞, 焦通, 郑玉生, 魏俏巧, 刘家橘. 连续随机离散缝网表征技术在中牟凹陷裂缝发育区中的应用[J]. 吉林大学学报(地球科学版), 2025, 55(5): 1715-1727.
[12] 田广, 赵岩. 信噪比约束的可调节振幅补偿算子反Q滤波方法[J]. 吉林大学学报(地球科学版), 2025, 55(4): 1351-1360.
[13] 王志勇, 刘国昌, 王梓旭, 郭严粮, 秦晨.

基于振幅一致性残差卷积编码-解码器的不规则缺失数据重建 [J]. 吉林大学学报(地球科学版), 2025, 55(4): 1336-1350.

[14] 孙敬雯, 吕子强, 孔庆翰, 唐泽豪, 邱俊辉, 刘珈君. 基于背景噪声成像的沂沭断裂带及邻区波速变化[J]. 吉林大学学报(地球科学版), 2025, 55(4): 1361-1371.
[15] 丁乾龙, 沈金松, 陈双全, 冉尚, 龙刚. 频率域传播矩阵法地震AVA多参数反演[J]. 吉林大学学报(地球科学版), 2025, 55(3): 957-969.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 吴远坤, 刘成林, 于春勇. 松辽盆地双城断陷深层原油成藏模式[J]. 吉林大学学报(地球科学版), 2024, 54(5): 1443 -1456 .
[2] 马荣,石建省,刘继朝. 人工内分泌网络模型在水文地质参数研究中的应用[J]. 吉林大学学报(地球科学版), 2013, 43(3): 914 -921 .
[3] 李宁, 王成文. 东北及邻区晚古生代地层接触关系与佳-蒙地块的形成和演化[J]. 吉林大学学报(地球科学版), 2017, 47(5): 1331 -1340 .
[4] 胡大千,初凤友,姚 杰. 中太平洋YJA海山富钴结壳矿物组成与元素地球化学[J]. J4, 2006, 36(01): 32 -0037 .
[5] 姜 雪, 程日辉,于民凤. 裂谷地层的气候和构造控制:Zscape模型分析与在松辽盆地北安断陷的应用[J]. J4, 2006, 36(01): 54 -0059 .
[6] 李春柏,张新涛,刘 立,任延广,孟 鹏. 布达特群热流体活动及其对火山碎屑岩的改造作用--以海拉尔盆地贝尔凹陷为例[J]. J4, 2006, 36(02): 221 -0226 .
[7] 孟宪纲,薄万举,刘志广,刘勇,畅柳,李朝柱,王子平. 芦山7.0级地震与巴颜喀拉块体中东段的活动性[J]. 吉林大学学报(地球科学版), 2014, 44(5): 1705 -1711 .
[8] 陈欢庆, 梁淑贤, 舒治睿, 邓晓娟, 彭寿昌. 冲积扇砾岩储层构型特征及其对储层开发的控制作用——以准噶尔盆地西北缘某区克下组冲积扇储层为例[J]. 吉林大学学报(地球科学版), 2015, 45(1): 13 -24 .
[9] 贾大成,邢立新, 潘 军, M. J. van Bergen, H. van Roermund. 伊通上地幔剪切带捕虏体中富铝尖晶石的地球化学特征[J]. J4, 2006, 36(04): 497 -502 .
[10] 谢忠雷,杨佰玲,包国章,董德明. 茶园土壤不同形态镍的含量及其影响因素[J]. J4, 2006, 36(04): 599 -604 .