3D reverse-time migration for pure P-wave in orthorhombic media  

作　　者：Ying-Hui Liu Jian-Ping Huang Liang Chen Qiang Mao Kun Tian 

机构地区：[1]School of Geoscience,China University of Petroleum(East China),Qingdao,266580,Shandong,China [2]Geophysical Research Institute,Shengli Oilfield Company,SINOPEC,Dongying,257022,Shandong,China

出　　处：《Petroleum Science》2024年第6期3937-3950,共14页石油科学(英文版)

基　　金：supported by the Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(No.2021QNLM020001);the Major Scientific and Technological Projects of Shandong Energy Group(No.SNKJ2022A06-R23);the Funds of Creative Research Groups of China(No.41821002);National Natural Science Foundation of China Outstanding Youth Science Fund Project(Overseas)(No.ZX20230152);the Major Scientific and Technological Projects of CNPC(No.ZD2019-183-003)。

摘　　要：Compared with the transverse isotropic(TI)medium,the orthorhombic anisotropic medium has both horizontal and vertical symmetry axes and it can be approximated as a set of vertical fissures developed in a group of horizontal strata.Although the full-elastic orthorhombic anisotropic wave equation can accurately simulate seismic wave propagation in the underground media,a huge computational cost is required in seismic modeling,migration,and inversion.The conventional coupled pseudo-acoustic wave equations based on acoustic approximation can be used to significantly reduce the cost of calculation.However,these equations usually suffer from unwanted shear wave artifacts during wave propagation,and the presence of these artifacts can significantly degrade the imaging quality.To solve these problems,we derived a new pure P-wave equation for orthorhombic media that eliminates shear wave artifacts while compromising computational efficiency and accuracy.In addition,the derived equation involves pseudo-differential operators and it must be solved by 3D FFT algorithms.In order to reduce the number of 3D FFT,we utilized the finite difference and pseudo-spectral methods to conduct 3D forward modeling.Furthermore,we simplified the equation by using elliptic approximation and implemented 3D reverse-time migration(RTM).Forward modeling tests on several homogeneous and heterogeneous models confirm that the accuracy of the new equation is better than that of conventional methods.3D RTM imaging tests on three-layer and SEG/EAGE 3D salt models confirm that the ORT media have better imaging quality.

关 键 词：Orthorho mbicmedia ANISOTROPIC Pure P-wave equation 3D reverse-time migration Graphic processing unit acceleration 

分 类 号：P631.4[天文地球—地质矿产勘探] P618.13[天文地球—地质学]