Seismic modeling by Gaussian beams in the global Cartesian coordinate system
Bai Min1, Wu Juan1, Sun Zhangqing2, Yang Lei3
1. School of Resources and Environment, North China University of Water Resources and Electric Power, Zhengzhou, Henan 450046, China; 2. College of Geoexploration Science and Technology, Jilin University, Changchun, Jilin 130026, China; 3. Yellow River Institute of Hydraulic Research, Yellow River Conservancy Commission, Zhengzhou, Henan 450003, China
Abstract:For the Gaussian beam in the ray-centered coordinate system, the coordinate transformation is needed for solving wave-field values at every point in the underground.To this end, Leung et al proposed a mathematical theory method to solve Gaussian beam in a global Cartesian coordinate system using the ansatz theory, which can be understood as deriving the traveling time along the central ray and adding an imaginary part as two-order Taylor expansion of central ray in its neighborhood, to construct an approximate travel time function.The approximate travel time function of conventional Gaussian beam can be obtained from the new approximate travel time function with a local ray-centered coordinate transformation.For the first time, the mathematical theory of Gaussian beam in the global Cartesian coordinate system is applied to seismic wave propagation.One Gaussian beam, Green's function, the single-frequency wave-field of continuous medium and wave-field of complex medium are simulated respectively.And the accuracy of Gaussian beam in the ray-center coordinate system and the global Cartesian coordinate system are compared.Results show that Green's function calculated by Gaussian beams in both of the coordinate systems approximates the analytic Green's function in a homogeneous medium.For continuous medium, Gaussian beam in the global Cartesian coordinate system can handle caustics problem and simulate wavefield propagation of complex media, but the accuracy will be reduced and the lateral energy is weak.
白敏, 吴娟, 孙章庆, 杨磊. 基于全局笛卡尔坐标系的高斯束地震波场模拟[J]. 石油地球物理勘探, 2018, 53(2): 257-265.
Bai Min, Wu Juan, Sun Zhangqing, Yang Lei. Seismic modeling by Gaussian beams in the global Cartesian coordinate system. Oil Geophysical Prospecting, 2018, 53(2): 257-265.
Červený V,Popov M M,Pšeník Ⅰ.Computation of wave fields in inhomogeneous media-Gaussian beam approach.Geophysical Journal International,1982,70(1):109-128.
Konopásková J,Červený V,Pšeník R I.Numerical modelling of time-harmonic seismic wave fields in simple structures by the Gaussian beam method-Part Ⅰ.Studia Geophysica et Geodaetica,1984,28(1):19-35.
[8]
Konopásková J,Červený V,Pšeník R I.Numerical modelling of time-harmonic seismic wave fields in simple structures by the Gaussian beam method-Part Ⅱ.Studia Geophysica et Geodaetica,1984,28(2):113-128.
[9]
Nowack R,Aki K.The two-dimensional Gaussian beam synthetic method:testing and application.Journal of Geophysical Research,1984,89:7797-7819.
[10]
Popov M M.A method of summation of Gaussian beams in isotropic theory of elasticity.Physics of the Solid Earth,1984,19:699-706.
[11]
Raz S.Beam stacking:A generalized preprocessing technique.Geophysics,1987,52(9):1199-1210.
[12]
Hanyga A.Gaussian beams in anisotropic elastic media.Geophysical Journal International,1986,85(3):473-504.
[13]
Liu Y X.Gaussian-beam Modeling in Heterogeneous Anisotropic Media[D].Colorado School of Mines,2010.
[14]
吴立明,许云,乌达巴拉.高斯束射线法在二维非均匀介质复杂构造中的应用.地球物理学报,1995,38(增刊1):144-152.Wu Liming,Xu Yun,Wudabala.Applied study of Gaussian beam method in 2-D inhomogeneous media and laterally varying layered structures.Chinese Journal of Geophysics,1995,38(S1):144-152.
[15]
刘学才,周熙襄,沙椿.三维高斯射线束法合成三分量VSP记录.石油地球物理勘探,1995,30(5):669-680.Liu Xuecai,Zhou Xixiang,Sha Chun.Synthesizing three-component VSP data by 3-D Gaussian beam method.OGP,1995,30(5):669-680.
[16]
张汝杰,贺振华,王理.井间地震高斯射线束正演方法.物探化探计算技术,1997,19(2):128-137.Zhang Rujie,He Zhenhua,Wang Li.Forward modeling by Gaussian beam method in crosswell seismic.Computing Techniques for Geophysical and Geochemical Exploration,1997,19(2):128-137.
[17]
邓飞,孙祥娥.三维高斯射线束正演的研究与实现.大庆石油地质与开发,2009,28(1):127-131.Deng Fei,Sun Xiang'e.Study and realization of 3D Gaussian beam forward modeling.Petroleum Geology and Oilfield Development in Daqing,2009,28(1):127-131.
[18]
孙成禹,张文颖,倪长宽等.能量约束下的高斯射线束法地震波场正演.石油地球物理勘探,2011,46(6):856-861.Sun Chengyu,Zhang Wenying,Ni Changkuan et al.Seismic wave field modeling by energy controlled Gaussian beam method.OGP,2011,46(6):856-861.
[19]
栗学磊.基于高斯束的弹性波场计算方法研究[学位论文].吉林长春:吉林大学,2013.
[20]
刘鹏,杨长春,王彦飞.三维TTI介质格林函数的高斯束计算方法.地球物理学进展,2013,28(5):2547-2553.Liu Peng,Yang Changchun,Wang Yanfei.Three-dimensional Green's function calculation in TTI media using the Gaussian beam method.Progress in Geophysics,2013,28(5):2547-2553.
[21]
吴娟.基于衰减补偿的高斯束偏移方法研究[学位论文].北京:中国石油大学(北京),2015.
[22]
George T,Virieux J,Madariaga R.Seismic wave synthesis by Gaussian beam summation:A comparison with finite differences.Geophysics,1987,52(8):1065-1073.
[23]
Hill N R.Gaussian beam migration.Geophysics,1990,55(11):1416-1428.
[24]
Červený V.Seismic rays and ray intensities in inhomogeneous anisotropic media.Geophysical Journal of the Royal Astronomical Society,1972,29(1):1-13.
[25]
Červený V,Klimeš L.Transformation relations for second-order derivatives of travel time in anisotropic media.Studia Geophysica Et Geodaetica,2010,54(2):257-267.
[26]
Ralston J.Gaussian beams and the propagation of singularities.Mathematics Association of America Studies in Mathematics,1983,23:206-248.
[27]
Tanushev N M,Qian J,Ralston J V.Mountain waves and Gaussian beams.Multiscale Modeling & Simulation,2007,6(2):688-709.
[28]
Leung S,Qian J,Burridge R.Eulerian Gaussian beams for high-frequency wave propagation.Geophysics,2007,72(5):SM61-SM76.
[29]
Gray S and May W.Kirchhoff migration using eikonal equation traveltimes.Geophysics,1994,59(5):810-817.
[30]
Qian J and Symes W W.An adaptive finite difference method for traveltime and amplitude.Geophysics,2002,67(1):167-176.
[31]
Symes W W and Qian J.A slowness matching Eulerian method for multivalued solutions of eikonal equations.Journal of Scientific Computing.2003,19(1-3):501-526.
[32]
Wu H,Yang X.The Eulerian Gaussian beam method for high frequency wave propagation in the reduced momentum space.Wave Motion,2013,50(6):1036-1049.
[33]
Yang Xu,Lu Jianfeng,Fomel S.Seismic modeling using the frozen Gaussian approximation.SEG Technical Program Expanded Abstracts,2013,32:52-58.