基于反射地震数据的时频域包络反演  

Reflection seismic data based envelope inversion in the time-frequency domain

作　　者：胡勇 黄兴国 陈玉婷 陈雪 张银涛[5] 万效国[5] 许永忠[1] 臧丽媛 杨磊[1] HU Yong;HUANG XingGuo;CHEN YuTing;CHEN Xue;ZHANG YinTao;WAN XiaoGuo;XU YongZhong;ZANG LiYuan;YANG Lei(School of Resources and Geosciences,China University of Mining and Technology,Xuzhou 221008,China;College of Instrumentation and Electrical Engineering,Jilin University,Changchun 130026,China;Exploration and Development Research Institute,PetroChina Huabei Oilfield Company,Renqiu 062552,China;Hebei GEO University,Hebei Key Laboratory of Strategic Critical Mineral Resources,Shijiazhuang 050031,China;Institute of China Petroleum Tarim Oilfield Company,Korla 841000,China)

机构地区：[1]中国矿业大学,资源与地球科学学院,徐州221008 [2]吉林大学,仪器科学与电气工程学院,长春130026 [3]中国石油华北油田,勘探开发研究院,任丘062552 [4]河北地质大学,河北省战略性关键矿产资源重点实验室,石家庄050031 [5]中国石油塔里木油田勘探开发研究院,库尔勒841000

出　　处：《地球物理学进展》2021年第5期1988-1994,共7页Progress in Geophysics

基　　金：江苏省自然科学基金(BK20200644);河北省自然科学基金资助(D2019403082)联合资助。

摘　　要：包络信号含有丰富的低频分量,即使在地震数据缺失低频条件下,包络目标函数也能有效缓解全波形反演的周期跳跃现象.但是,当初始速度模型较为平滑时,观测数据中的反射地震事件在模拟数据中没有与之相对应的波形,导致包络反演初期无法很好地利用反射波信号进行速度建模.本文提出基于反射地震数据的时频域包络反演方法,通过结合反射波全波形反演理论,构建反射波时频域包络目标函数,来提高包络反演的速度建模精度.本文首先利用Gabor变换获取时频域地震数据,并提取振幅信息,即为时频域包络信号.然后,推导反射波时频域包络反演的伴随震源和梯度算子.Marmousi模型数据测试结果表明,基于反射地震数据的时频域包络反演方法可以为全波形反演提供一个较好的初始速度模型.Envelope signals contains abundant low-frequency components, even if the seismic data lack low-frequencies, envelope misfit can effectively mitigate the cycle skipping problem of Full Waveform Inversion(FWI). However, when the initial velocity model is smooth, the reflected seismic event in the observed data does not have a corresponding waveform in the synthetic data, which results in the inability to make good use of the reflection data for velocity modeling in the early stage of envelope inversion. This paper proposes to use the reflection seismic data to conduct the envelope inversion in the time-frequency domain. By combining the theory of reflection FWI, the time-frequency domain envelope misfit function is established to improve the velocity modeling accuracy. This article first uses Gabor transform to obtain time-frequency domain seismic data, then extracts the amplitude information, which is regarded as time-frequency domain envelope signals. After that, the gradient operator and adjoint source of the time-frequency reflection envelope inversion are derived. Finally, the Marmousi model tests show that the reflection seismic data based envelope inversion can provide a better initial velocity model for the FWI method.

关 键 词：反射波全波形反演 包络反演 目标函数 时频域 速度模型 

分 类 号：P631[天文地球—地质矿产勘探]