Seismic anisotropy has been extensively acknowledged as a crucial element that influences the wave propagation characteristic during wavefield simulation,inversion and imaging.Transversely isotropy(TI)and orthorhombic...Seismic anisotropy has been extensively acknowledged as a crucial element that influences the wave propagation characteristic during wavefield simulation,inversion and imaging.Transversely isotropy(TI)and orthorhombic anisotropy(OA)are two typical categories of anisotropic media in exploration geophysics.In comparison of the elastic wave equations in both TI and OA media,pseudo-acoustic wave equations(PWEs)based on the acoustic assumption can markedly reduce computational cost and complexity.However,the presently available PWEs may experience SV-wave contamination and instability when anisotropic parameters cannot satisfy the approximated condition.Exploiting pure-mode wave equations can effectively resolve the above-mentioned issues and generate pure P-wave events without any artifacts.To further improve the computational accuracy and efficiency,we develop two novel pure qP-wave equations(PPEs)and illustrate the corresponding numerical solutions in the timespace domain for 3D tilted TI(TTI)and tilted OA(TOA)media.First,the rational polynomials are adopted to estimate the exact pure qP-wave dispersion relations,which contain complicated pseudo-differential operators with irrational forms.The polynomial coefficients are produced by applying a linear optimization algorithm to minimize the objective function difference between the expansion formula and the exact one.Then,the developed optimized PPEs are efficiently implemented using the finite-difference(FD)method in the time-space domain by introducing a scalar operator,which can help avoid the problem of spectral-based algorithms and other calculation burdens.Structures of the new equations are concise and corresponding implementation processes are straightforward.Phase velocity analyses indicate that our proposed optimized equations can lead to reliable approximation results.3D synthetic examples demonstrate that our proposed FD-based PPEs can produce accurate and stable P-wave responses,and effectively describe the wavefield features in complicated TTI and TOA media.展开更多
The anisotropic properties of subsurface media cause waveform distortions in seismic wave propagation,resulting in a negative infl uence on seismic imaging.In addition,wavefields simulated by the conventional coupled ...The anisotropic properties of subsurface media cause waveform distortions in seismic wave propagation,resulting in a negative infl uence on seismic imaging.In addition,wavefields simulated by the conventional coupled pseudo-acoustic equation are not only aff ected by SV-wave artifacts but are also limited by anisotropic parameters.We propose a least-squares reverse time migration(LSRTM)method based on the pure q P-wave equation in vertically transverse isotropic media.A fi nite diff erence and fast Fourier transform method,which can improve the effi ciency of the numerical simulation compared to a pseudo-spectral method,is used to solve the pure q P-wave equation.We derive the corresponding demigration operator,migration operator,and gradient updating formula to implement the LSRTM.Numerical tests on the Hess model and field data confirm that the proposed method has a good correction eff ect for the travel time deviation caused by underground anisotropic media.Further,it signifi cantly suppresses the migration noise,balances the imaging amplitude,and improves the imaging resolution.展开更多
基金supported by the National Key R&D Program of China(2021YFA0716902)National Natural Science Foundation of China(NSFC)under contract number 42374149 and 42004119National Science and Technology Major Project(2024ZD1002907)。
文摘Seismic anisotropy has been extensively acknowledged as a crucial element that influences the wave propagation characteristic during wavefield simulation,inversion and imaging.Transversely isotropy(TI)and orthorhombic anisotropy(OA)are two typical categories of anisotropic media in exploration geophysics.In comparison of the elastic wave equations in both TI and OA media,pseudo-acoustic wave equations(PWEs)based on the acoustic assumption can markedly reduce computational cost and complexity.However,the presently available PWEs may experience SV-wave contamination and instability when anisotropic parameters cannot satisfy the approximated condition.Exploiting pure-mode wave equations can effectively resolve the above-mentioned issues and generate pure P-wave events without any artifacts.To further improve the computational accuracy and efficiency,we develop two novel pure qP-wave equations(PPEs)and illustrate the corresponding numerical solutions in the timespace domain for 3D tilted TI(TTI)and tilted OA(TOA)media.First,the rational polynomials are adopted to estimate the exact pure qP-wave dispersion relations,which contain complicated pseudo-differential operators with irrational forms.The polynomial coefficients are produced by applying a linear optimization algorithm to minimize the objective function difference between the expansion formula and the exact one.Then,the developed optimized PPEs are efficiently implemented using the finite-difference(FD)method in the time-space domain by introducing a scalar operator,which can help avoid the problem of spectral-based algorithms and other calculation burdens.Structures of the new equations are concise and corresponding implementation processes are straightforward.Phase velocity analyses indicate that our proposed optimized equations can lead to reliable approximation results.3D synthetic examples demonstrate that our proposed FD-based PPEs can produce accurate and stable P-wave responses,and effectively describe the wavefield features in complicated TTI and TOA media.
基金financially supported by the National Key R&D Program of China (No. 2019YFC0605503)the Major Scientific and Technological Projects of CNPC (No. ZD2019-183-003)the National Natural Science Foundation of China (No. 41922028,41874149)。
文摘The anisotropic properties of subsurface media cause waveform distortions in seismic wave propagation,resulting in a negative infl uence on seismic imaging.In addition,wavefields simulated by the conventional coupled pseudo-acoustic equation are not only aff ected by SV-wave artifacts but are also limited by anisotropic parameters.We propose a least-squares reverse time migration(LSRTM)method based on the pure q P-wave equation in vertically transverse isotropic media.A fi nite diff erence and fast Fourier transform method,which can improve the effi ciency of the numerical simulation compared to a pseudo-spectral method,is used to solve the pure q P-wave equation.We derive the corresponding demigration operator,migration operator,and gradient updating formula to implement the LSRTM.Numerical tests on the Hess model and field data confirm that the proposed method has a good correction eff ect for the travel time deviation caused by underground anisotropic media.Further,it signifi cantly suppresses the migration noise,balances the imaging amplitude,and improves the imaging resolution.
