局部尺度频率波数域相位反演  

Local-scale frequency-wavenumber domain phase inversion

作　　者：康鹏飞 胡勇 刘瑞东 孙冲 赵卓 袁萍 郑明君 蒙毅凡 许永忠[2] KANG PengFei;HU Yong;LIU RuiDong;SUN Chong;ZHAO Zhuo;YUAN Ping;Zhen MingJun;MENG YiFan;XU YongZhong(Institute of China Petroleum Tarim Oilfield Company,Korla 841000,China;School of Resources and Geosciences,China University of Mining and Technology,Xuzhou 221008,China;Shandong Provincial Key Laboratory of Deep Oil and Gas,China University of Petroleum(East China),Qingdao 266580,China;R&D Center for Ultra-Deep Complex Reservior Exploration and Development,CNPC,Korla 841000,China)

机构地区：[1]中石油塔里木油田勘探开发研究院,库尔勒841000 [2]中国矿业大学资源与地球科学学院,徐州221008 [3]中国石油大学(华东),山东省深层油气重点实验室,青岛266580 [4]中国石油天然气集团有限公司超深层复杂油气藏勘探开发技术研发中心,库尔勒841000

出　　处：《地球物理学进展》2025年第1期155-165,共11页Progress in Geophysics

基　　金：国家自然科学基金项目(42104116);中国博士后科学基金(2021M703580);企业重大技术攻关项目(T202201.04,041022070109)联合资助.

摘　　要：全波形反演(FWI)是获得高分辨率速度模型的一种重要方法.但是当地震数据缺失低频分量时,常规FWI结果会出现严重的周期跳跃问题,影响最终速度建模的精度.为此,本文提出了局部尺度频率波数域相位反演方法,充分考虑了地震数据在时间和偏移距方向上的局部特征,并利用地震数据局部尺度分解策略和二维傅里叶变换,来构建局部尺度频率波数域相位目标函数,恢复速度模型的低波数分量,为常规FWI方法提供一个较好的初始速度模型.本文首先利用二维滑动窗函数来提取局部尺度地震数据,结合二维傅里叶变换建立局部尺度频率波数域指数相位差目标函数.然后,推导局部尺度频率波数域相位反演对应的伴随震源和梯度算子.最后,Marmousi模型和火成岩-碳酸盐模型测试结果表明,局部尺度频率波数域相位反演方法可以为FWI提供一个较好的初始速度模型,实现缓解FWI周期跳跃的目标.Full Waveform Inversion(FWI)is an important method to obtain high-resolution velocity models.However,the seismic data are absence of low-frequency components,the conventional FWI results will have serious cycle skipping problem,which affects the accuracy of the final velocity modeling.For this reason,this paper proposes a local-scale frequency-wavenumber domain phase inversion method,which fully considers the local features of seismic data in the time and offset directions,and utilizes the local-scale decomposition strategy of seismic data and 2D Fourier transform to construct a local-scale frequency-wavenumber domain phase misfit function,to recover the low-wavenumber components of the velocity model and to provide a better initial velocity model for the FWI method.In this paper,we first utilize a 2D sliding window function to extract the local-scale seismic data,and combine it with the 2D Fourier transform to establish the local-scale frequency-wavenumber domain phase information based misfit function.Then,we derived the adjoint-source and gradient operator corresponding to the phase inversion in the local-scale frequency-wavenumber domain.Finally,the test results of the Marmousi model and the igneous-carbonate model show that the local-scale frequency-wavenumber domain phase inversion can provide a better initial velocity model for the conventional FWI method and mitigate the FWI cycle skipping.

关 键 词：频率波数域 全波形反演 相位反演 局部尺度 目标函数 

分 类 号：TP3[自动化与计算机技术—计算机科学与技术]