Please wait a minute...
E-mail Alert Rss
 
物探与化探  2025, Vol. 49 Issue (5): 1155-1163    DOI: 10.11720/wtyht.2025.0043
  方法研究信息处理仪器研制 本期目录 | 过刊浏览 | 高级检索 |
反射波时频域地震波形反演方法
汪昆(), 吴国忱(), 贾宗锋, 杨凌云
中国石油大学 (华东) 地球科学与技术学院,山东 青岛 266580
Reflection waveform inversion in the time-frequency domain
WANG Kun(), WU Guo-Chen(), JIA Zong-Feng, YANG Ling-Yun
School of Geosciences,China University of Petroleum(East China),Qingdao 266580,China
全文: PDF(2790 KB)   HTML
输出: BibTeX | EndNote (RIS)      
摘要 

反射波形反演通过交替更新模型中的中低波数组分和高波数组分为全波形反演提供良好的初始模型。但时域反射波形反演需要存储整个时间序列的震源背景波场与接收点散射波场互相关求取梯度,占用巨量的计算机存储空间。频域反射波形反演虽然有多尺度特性,但对计算机运算内存要求高。本文通过离散傅里叶变换,从时域波场提取对应的频域波场进行多尺度的时频域反演,仅需存储数片单频波场,在存储需求方面明显优于传统的时域反演,运算成本方面明显优于频率域,结合了时域反演的高效性和频域反演多尺度的特性。针对实际地震资料低频地震数据易缺失的问题,进一步结合包络反射波反演,在重构中深层地下介质低频信息的基础上进行逐级优化反演,降低对低频长偏移地震数据的要求。最后本文通过数值算例验证了基于包络数据的时频域反射波形反演方法对低频信息恢复的有效性。

服务
把本文推荐给朋友
加入引用管理器
E-mail Alert
RSS
作者相关文章
汪昆
吴国忱
贾宗锋
杨凌云
关键词 全波形反演反射波反演时频域包络反演    
Abstract

Reflection waveform inversion(RWI) provides an effective initial model for full waveform inversion(FWI) by alternately updating the low-to-intermediate and high wavenumber components in the model.However,time-domain RWI requires storing the cross-correlation between the source background wavefield and the receiver scattered wavefield over the entire time series to compute gradients,demanding substantial computational storage.Although frequency-domain RWI exhibits multi-scale properties,it also imposes high demands on computational memory.Based on the discrete Fourier transform,this study proposed a RWI method that extracts frequency-domain wavefields from corresponding time-domain wavefields for multi-scale inversion in the time-frequency domain.The proposed method requires storing only a few single-frequency wavefield snapshots,showing significantly lower storage demands compared to conventional time-domain RWI and reduced computational costs relative to frequency-domain RWI.Therefore,the proposed method effectively combines the computational efficiency of time-domain RWI with the multi-scale properties of frequency-domain RWI.Considering the frequently missing low-frequency data in actual seismic data,this study further integrated envelope-based RWI to reconstruct low-frequency information for medium-deep subsurface structures.This enables stage-wise optimization of the inversion process,reducing the dependency on low-frequency long-offset data.Finally,numerical examples validate the effectiveness of the time-frequency domain RWI method based on envelope data in recovering low-frequency information.

Key wordsfull waveform inversion(FWI)    reflection waveform inversion(RWI)    time-frequency domain    envelope inversion
收稿日期: 2025-02-25      修回日期: 2025-07-03      出版日期: 2025-10-20
ZTFLH:  P631.4  
基金资助:国家自然科学基金重点项目“深部裂缝储层地震预测理论方法及关键技术研究”(42430809)
通讯作者: 吴国忱(1965-),男,教授,博士生导师,现从事地球物理方面的教学和科研工作。Email :guochenwu@upc.edu.cn
作者简介: 汪昆(1998-),男,硕士研究生,2020年毕业于长安大学,获地球物理学学士学位;现于中国石油大学(华东)攻读地质工程专业硕士学位,主要从事地震波正演模拟方面的学习和研究工作。Email:1514073707@qq.com
引用本文:   
汪昆, 吴国忱, 贾宗锋, 杨凌云. 反射波时频域地震波形反演方法[J]. 物探与化探, 2025, 49(5): 1155-1163.
WANG Kun, WU Guo-Chen, JIA Zong-Feng, YANG Ling-Yun. Reflection waveform inversion in the time-frequency domain. Geophysical and Geochemical Exploration, 2025, 49(5): 1155-1163.
链接本文:  
https://www.wutanyuhuatan.com/CN/10.11720/wtyht.2025.0043      或      https://www.wutanyuhuatan.com/CN/Y2025/V49/I5/1155
Fig.1  反射波路径示意
Fig.2  层状模型
Fig.3  时频域与频率域单频波场对比
Fig.4  时频域与频率域单频波场抽道对比
Fig.5  平层速度模型
Fig.6  Sigebee2A速度模型
Fig.7  单道速度反演结果对比
单核/min 十核/min 存储/GB
时间域 185.3 16.3 85.8
时频域 197.9 17.5 14.1
Table 1  效率测试
Fig.8  复杂速度模型
Fig.9  单道速度反演结果对比
[23] Wu R S, Luo J R, Wu B Y. Seismic envelope inversion and modulation signal model[J]. Geophysics, 2014, 79(3):WA13-WA24.
[24] Chi B X, Dong L G, Liu Y Z. Full waveform inversion method using envelope objective function without low frequency data[J]. Journal of Applied Geophysics, 2014,109:36-46.
[25] 梁展源, 吴国忱, 张晓语. 基于频移目标函数的包络反演方法[J]. 地球物理学进展, 2019, 34(4):1481-1488.
[25] Liang Z Y, Wu G C, Zhang X Y. Envelope inversion method based on frequency-shifted objective function[J]. Progress in Geophysics, 2019, 34(4):1481-1488.
[1] Lailly P. The seismic inverse problem as a sequence of before stack migrations[C]// Conference on Inverse Scattering,Theory and Applications, Society for Industrial and Applied Mathematics,1983:206-220.
[2] Tarantola A. Inversion of seismic reflection data in the acoustic approximation[J]. Geophysics, 1984, 49(8):1259-1266.
[3] 刘畅, 李振春, 曲英铭, 等. 地震层析成像方法综述[J]. 物探与化探, 2020, 44(2):227-234.
[3] Liu C, Li Z C, Qu Y M, et al. A review of seismic tomography methods[J]. Geophysical and Geochemical Exploration, 2020, 44(2):227-234.
[4] 国运东. 基于组合震源编码的多尺度全波形反演方法[J]. 物探与化探, 2022, 46(3):729-736.
[4] Guo Y D. Multi-scale full waveform inversion method using combined source encoding[J]. Geophysical and Geochemical Exploration, 2022, 46(3):729-736.
[5] Xu S, Wang D L, Chen F, et al. Inversion on reflected seismic wave[C]// SEG Technical Program Expanded Abstracts 2012, Society of Exploration Geophysicists,2012:1-7.
[6] Wang S, Chen F, Zhang H Z, et al. Reflection-based full waveform inversion (RFWI) in the frequency domain[C]// SEG Technical Program Expanded Abstracts 2013, Society of Exploration Geophysicists,2013:877-881.
[7] Chi B X, Dong L G, Liu Y Z. Correlation-based reflection full-waveform inversion[J]. Geophysics, 2015, 80(4):R189-R202.
[8] Luo Y, Ma Y, Wu Y, et al. Full-traveltime inversion[J]. Geophysics, 2016, 81(5):R261-R274.
[9] 姚刚, 吴迪. 反射波全波形反演[J]. 中国科学:地球科学, 2017, 47(10):1220-1232.
[9] Yao G, Wu D. Reflection full waveform inversion[J]. Scientia Sinica:Terrae, 2017, 47(10):1220-1232.
[10] Guo Q, Alkhalifah T. Elastic reflection-based waveform inversion with a nonlinear approach[J]. Geophysics, 2017, 82(6):R309-R321.
[11] Wang T F, Cheng J B, Guo Q, et al. Elastic wave-equation-based reflection kernel analysis and traveltime inversion using wave mode decomposition[J]. Geophysical Journal International, 2018, 215(1):450-470.
[12] Li Y Y, Guo Q, Li Z C, et al. Elastic reflection waveform inversion with variable density[J]. Geophysics, 2019, 84(4):R553-R567.
[13] 李青阳, 吴国忱, 段沛然, 等. 基于互相关目标函数的反射波波形反演[J]. 石油地球物理勘探, 2020, 55(4):754-765,700.
[13] Li Q Y, Wu G C, Duan P R, et al. Reflection waveform inversion based on cross-correlation misfit function[J]. Oil Geophysical Prospecting, 2020, 55(4):754-765,700.
[14] Bunks C, Saleck F M, Zaleski S, et al. Multiscale seismic waveform inversion[J]. Geophysics, 1995, 60(5):1457-1473.
[15] Boonyasiriwat C, Valasek P, Routh P, et al. An efficient multiscale method for time-domain waveform tomography[J]. Geophysics, 2009, 74(6):WCC59-WCC68.
[16] Pratt R G. Seismic waveform inversion in the frequency domain,Part 1:Theory and verification in a physical scale model[J]. Geophysics, 1999, 64(3):888-901.
[17] Pratt R G, Shipp R M. Seismic waveform inversion in the frequency domain,Part 2:Fault delineation in sediments using crosshole data[J]. Geophysics, 1999, 64(3):902-914.
[18] Operto S, Virieux J, Amestoy P, et al. 3D finite-difference frequency-domain modeling of visco-acoustic wave propagation using a massively parallel direct solver:A feasibility study[J]. Geophysics, 2007, 72(5):SM195-SM211.
[19] Fichtner A, Kennett B L N, Igel H, et al. Theoretical background for continental- and global-scale full-waveform inversion in the time-frequency domain[J]. Geophysical Journal International, 2008, 175(2):665-685.
[20] 王官超. 基于时域正演的弹性波频域全波形反演方法[D]. 东营: 中国石油大学(华东), 2016.
[20] Wang G C. Full waveform inversion method of elastic wave in frequency domain based on time domain forward modeling[D]. Dongying: China University of Petroleum(Huadong), 2016.
[21] 胡勇, 韩立国, 于江龙, 等. 基于自适应非稳态相位校正的时频域多尺度全波形反演[J]. 地球物理学报, 2018, 61(7):2969-2988.
doi: 10.6038/cjg2018L0421
[21] Hu Y, Han L G, Yu J L, et al. Time-frequency domain multi-scale full waveform inversion based on adaptive non-stationary phase correction[J]. Chinese J. Geophys., 2018, 61(7):2969-2988.
[22] Wang J H, Dong L G, Yang J Z, et al. Monochromatic wave equation reflection traveltime inversion[J]. Geophysics, 2024, 89(1):S83-S97.
[26] Luo J R, Wu R S, Hu Y, et al. Strong scattering elastic full waveform inversion with the envelope Fréchet derivative[J]. IEEE Geoscience and Remote Sensing Letters, 2021,19:8008805.
[1] 国运东. 川东北探区层析与全波形反演联合建模方法及应用[J]. 物探与化探, 2025, 49(5): 1090-1098.
[2] 陈子龙, 王海燕, 郭华, 王光文, 赵玉莲. 地震全波形反演研究进展与应用现状综述[J]. 物探与化探, 2023, 47(3): 628-637.
[3] 王莉利, 杜功鑫, 高新成, 王宁, 王维红. 基于U-Net网络的FWI地震低频恢复方法[J]. 物探与化探, 2023, 47(2): 391-400.
[4] 石战战, 庞溯, 王元君, 池跃龙, 周强. 基于f-x域时频非凸正则化低秩矩阵近似的共偏移距道集去噪方法[J]. 物探与化探, 2022, 46(6): 1444-1453.
[5] 姚含, 徐海. 基于梯度投影法的全变差正则化全波形反演[J]. 物探与化探, 2022, 46(4): 977-981.
[6] 国运东. 基于组合震源编码的多尺度全波形反演方法[J]. 物探与化探, 2022, 46(3): 729-736.
[7] 周斯琛, 李振春, 张敏, 张凯. 基于截断牛顿法的频率域全波形反演方法[J]. 物探与化探, 2017, 41(1): 147-152.
[8] 苗永康. 叠前地震全波形反演实际应用的关键影响因素[J]. 物探与化探, 2016, 40(5): 947-954.
[9] 尚新民, 刁瑞, 冯玉苹, 赵翠霞. 谱模拟方法在高分辨率地震资料处理中的应用[J]. 物探与化探, 2014, 38(1): 75-80.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
京ICP备05055290号-3
版权所有 © 2021《物探与化探》编辑部
通讯地址:北京市学院路29号航遥中心 邮编:100083
电话:010-62060192;62060193 E-mail:whtbjb@sina.com , whtbjb@163.com