Efficient anti-aliasing and anti-leakage Fourier transform for high-dimensional seismic data regularization using cube removal and GPU

作     者：Lu Liu Sindi Ghada Fu-Hao Qin Youngseo Kim Vladimir Aleksic Hong-Wei Liu Lu Liu;Sindi Ghada;Fu-Hao Qin;Youngseo Kim;Vladimir Aleksic;Hong-Wei Liu

作者机构：Aramco Research Center-BeijingBeijing100102China KAUST Research Center GroupSaudi AramcoThuwal23955-6900Saudi Arabia EXPEC Advanced Research CenterSaudi AramcoDhahran31311Saudi Arabia Geophysical Imaging DepartmentSaudi AramcoDhahran31311Saudi Arabia National Key Laboratory of Deep Oil and GasChina University of Petroleum(East China)Qingdao266580ShandongChina 

出 版 物：《Petroleum Science》 (石油科学（英文版）)

年 卷 期：2024年第21卷第5期

页      面：3079-3089页

核心收录：

学科分类：081801[工学-矿产普查与勘探] 081802[工学-地球探测与信息技术] 08[工学] 0818[工学-地质资源与地质工程] 

基　　金：We appreciate the helpful discussions and comments of our colleagues Yi He  Yu-Jin Liu  Yue Ma  and Yi Luo. We also thank three anonymous reviewers for their useful suggestions 

主　　题：High-dimensional regularization GPU Anti-aliasing Anti-leakage 

摘      要：Seismic data is commonly acquired sparsely and irregularly, which necessitates the regularization of seismic data with anti-aliasing and anti-leakage methods during seismic data processing. We propose a novel method of 4D anti-aliasing and anti-leakage Fourier transform using a cube-removal strategy to address the combination of irregular sampling and aliasing in high-dimensional seismic data. We compute a weighting function by stacking the spectrum along the radial lines, apply this function to suppress the aliasing energy, and then iteratively pick the dominant amplitude cube to construct the Fourier spectrum. The proposed method is very efficient due to a cube removal strategy for accelerating the convergence of Fourier reconstruction and a well-designed parallel architecture using CPU/GPU collaborative computing. To better fill the acquisition holes from 5D seismic data and meanwhile considering the GPU memory limitation, we developed the anti-aliasing and anti-leakage Fourier transform method in 4D with the remaining spatial dimension looped. The entire workflow is composed of three steps: data splitting, 4D regularization, and data merging. Numerical tests on both synthetic and field data examples demonstrate the high efficiency and effectiveness of our approach.