文献类型：期刊文献

英文题名：One-way propagators based on matrix multiplication in arbitrarily lateral varying media with GPU implementation

作者：You, Jiachun[1,2];Liu, Zhege[1];Liu, Jianli[1];Li, Canping[3]

机构：[1]Chengdu Univ Technol, State Key Lab Oil & Gas Reservoir Geol & Exploita, Chengdu 610059, Sichuan, Peoples R China;[2]Shandong Univ Sci & Technol, Shandong Key Lab Deposit Mineralizat & Sedimentar, Qingdao 266590, Shandong, Peoples R China;[3]Guangdong Ocean Univ, Lab Ocean Remote Sensing & Informat Technol, Zhanjiang 524088, Peoples R China

年份：2019

卷号：130

起止页码：32

外文期刊名：COMPUTERS & GEOSCIENCES

收录：SCI-EXPANDED(收录号：WOS:000478712700004)、、EI(收录号：20192407047091)、Scopus(收录号：2-s2.0-85067042601)、WOS

基金：This work was supported by China Postdoctoral Science Foundation (Grant No.2019M653357) and Shandong Provincial Key Laboratory of Depositional Mineralization & Sedimentary Mineral in Shandong University of Science and Technology (Grant No. DMSM2018015) and the National Natural Science Foundation of China (Grant No. 41306050).

语种：英文

外文关键词：One-way wave equation; Matrix multiplication; Helmholtz operator; Square root operator; Parallel implementation

外文摘要：Conventional one-way migration methods have some shortcomings, such as their limitations in imaging large angles. A novel strategy is presented here to solve these issues. The Helmholtz operator of the two-way wave equation is a power of the square root operator of the one-way wave equation. Based on this connectivity, the computation of the square root of the Helmholtz operator based on matrix multiplication and the calculation of a one-way propagator with a Taylor series expansion, is studied. The one-way propagators used in this paper involve complete matrix multiplication, which is a distinct feature of this method and is suitable for parallel implementation. Analysing the values of matrix square roots generated by the classical eigenvalues decomposition method and an algorithm of computing matrix square roots with matrix multiplication, we prove that our algorithm can calculate a stable result, which lays a solid foundation for further study. The impulse response test proves the advantage of our proposed method over conventional one-way migration in describing the wavefield propagation of large angles in a medium with a strong lateral velocity variation. Moreover, we implement GPU and CPU versions of the algorithm and compared their efficiencies. The post-stack migration result of the Salt model and the pre-stack migration result of the Marmousi model further illustrate the superiority of our proposed migration method in imaging the complex and fine-scale structures compared to the conventional one-way migration method. This provides a promising practical application in seismic exploration.