Research on Wave Field Numerical Simulation of High Order Finite Difference in Multi-Scale Grid Wave Equations
[1]
Zhang Xiaodan, College of Electronic and Information, Xi`an Polytechnic University, Xi’an, China; College of Electronics and Information Engineering (College of Microelectronics), Xi`an Jiaotong University, Xi’an, China.
[2]
She Yichong, College of Electronic and Information, Xi`an Polytechnic University, Xi’an, China.
[3]
Liu Guizhong, College of Electronics and Information Engineering (College of Microelectronics), Xi`an Jiaotong University, Xi’an, China.
[4]
Zhang Zhiyu, College of Automation and Information engineering, Xi'an University of Technology, Xi’an, China.
[5]
Zhu Lei, College of Electronic and Information, Xi`an Polytechnic University, Xi’an, China.
In the numerical simulation of seismic wave field, the problem of how to ensure both high efficiency and precision has always been one of the hot spots of seismic exploration scholars. The traditional method used the constant small step length in finite difference, which greatly reduces the calculation efficiency. A method of adopt different scale grids according to the characteristics of the geological model and optimize the transition zone has been proposed. firstly, analysis the speed model of the research object to determine the scale of grid; secondly, determine the scope of the transition zone; finally, calculate the coefficient and the differential points of both inside and outside the transition zone, gain the wave field value of every grid point of the model. According to the experimental results in the paper, the calculation efficiency of multi-scale grid method can be improved obviously, and the case’s results of this article can as high as 25.16% in average.
Multi-scale Grid, Step Length, Transition Zone, Wave Field Simulation, Computational Efficiency
[1]
Tian Renfei, Cao Junxing, The efficiency improvement of the pseudo-spectral method and the simulation of the cross-well seismic wave modeling [J]. Geophysical and Geochemical Exploration: 2008, 4: 430-433.
[2]
Li Lin, Liu Tao, Hu Tianyue, Spectral element method with triangular mesh and its application in seismic modeling [J]. Chinese Journal Of Geophysics: 2014, 4: 1224-1234.
[3]
DU Huakun, FENG Deshan, GPR simulation by finite element method of unstructured grid based on Delaunay triangulation [J]. Journal of Central South University (Science and Technology), 2015, 46 (4): 1326-1334.
[4]
Lian Ximeng, Shan Lianyu, Sui Zhiqiang, An overview of research on perfectly matched layers absorbing boundary condition of seismic forward numerical simulation [J]. Progress In Geophysics: 2015, 30 (4): 1725-1733 1004-2903.
[5]
Zhang Geng, Tuo Xianguo, Zhang Yi, Forward modeling study on data acquisition parameters of shallow seismic reflection method [J]. Geotechnical Investigation&Surveying: 2015, (1): 93-98, 1000-1433.
[6]
Huang Jianping, Yang Jidong, Li Zhenchun, Gaussian beam seismic forward modeling for 3 dimensional undulating surface under the Effective neighborhood wave field approximation framework [J]. Oil Geophysical Prospecting: 2015, 50 (5) P896-904, 804 1000-7210.
[7]
Jian-Ping Huang; Ying-Ming Qu, Variable-coordinate forward modeling of irregular surface based on dual-variable grid [J]. Applied Geophysics: 2015, 12 (1): 101-110.
[8]
Bai, Chao-ying; Hu, Guang-yi, Seismic wavefield propagation in 2D anisotropic media: Ray theory versus wave-equation simulation [J]. Journal of Applied Geophysics: 2014, 104: 167-171.
[9]
Cai Xiao-Hui, Liu Yang, Ren Zhi-Ming, Three—dimensional acoustic wave equation modeling based on the optimal finite-difference scheme [J]. APPLIED GEOPHYSICS: 2015, 12 (3): 409—420.
[10]
Li Bin, Wen Mingming, Mu Zelin, Staggered-grid Finite-difference Method with High-order Accuracy in Time-space Domains for Acoustic Forward Modeling [J]. Computerized Tomography Theory And Application: 2017, 3: 259-266.
[11]
Yue Xiaopeng, Bai Chaoyin, Yue Chongwang, Accuracy analysis of elastic wave field simulation based on high-order staggered grid finite difference scheme [J]. Coal Geology & Exploration: 2017, 1, 125-130.
[12]
X Yue; C Yue, Simulation of acoustic wave propagation in a borehole surrounded by cracked media using a finite difference method based on Hudson’s approach [J]. Journal of Geophysics and Engineering. 2017, 14: 633–640.
[13]
Wei-Zhong Wang, Tian-Yue Hu, x, Varia-ble-order rotated staggered-grid method for elastic-wave forward modeling [J]. Applied Geophysics: 2015, 12 (3): 389-400.
[14]
Hongyong Yan, Lei Yang, Seismic modeling with an optimal staggered-grid finite-difference scheme based on combining Taylor-series expansion and minimax approximation [J]. Studia Geophysica et Geodaetica: 2017, 61 (3): 560-574.
[15]
Wang Jian, Meng Xiaohong, Liu Hong, Cosine-modulated window function-based staggered-grid finite-difference forward modeling [J]. Applied Geophysics: 2017, 1: 115-124.
[16]
Wang Jian; Meng Xiaohong; Liu Hong; Optimization of finite difference forward modeling for elastic waves based on optimum combined window functions [J]. Journal of Applied Geophysics, 2017, 138: 62-71.
[17]
Feng Jihao, Cao Danping, Qin Haixu, Analysis of reflection characteristics of multi-scale seismic data based on the wave equation forward modeling [J]. Progress In Geophysics: 2016 (3): 1058-1065 1004-2903.
[18]
Li Yuesheng, Wu Guochen, Seismic modeling in three-dimensional monoclinic fractured media [J]. Progress In Geophysics: 2016, 2: 525-536 1004-2903.
[19]
Liang Wenquan, Liao Yongming, Jiang Lin, Acoustic Wave Equation Modeling with a New Time-Space Domain Finite Difference Stencil and an Improved Linear Algorithm [J]. China Earthquake Engineering Journal: 2016, 5: 815-821.
[20]
Wang, Yanfei; Liang, Wenquan; Nashed, Zuhair; Determination of finite difference coefficients for the acoustic wave equation using regularized least-squares inversion [J]. J. Inverse Ill-Posed Probl, 2016, 24 (6): 743-760.
[21]
Peter Moczo, Finite-difference technique for SH-waves in 2-D media using irregular grids-application to the seismic response problem [J]. Geophysical Journal International: 1989, 99: 321-329.
[22]
Emmanuel Chaljub; Emeline Maufroy; Peter Moczo; 3-D numerical simulations of earthquake ground motion in sedimentary basins: testing accuracy through stringent models [J]. Geophysical Journal International, 2015, 201 (1): 90-111.
[23]
Jastram C, Tessmer E. Elastic modeling on a grid with vertically varying spacing [J]. Geophysics Prospecting: 1994, 42 (4): 357-370.
[24]
Wang Y, and Schuster G T. Finite difference variable grid scheme for acoustic and elastic wave equation modeling [C]. 1996 SEG annual meeting., Soc. Expl. Geophysics Expanded Abstracted: 1996, 674-677.
[25]
Zhang Jianfeng, Non-regular grid difference method for elastic wave numerical simulation. [J]. Chinese Journal Of Geophysics 1998, 41 (supplementary issue): 357-366.
[26]
Dong Liangguo, Li Peiming, Dispersive Problem in Seismic Wave Propagation Numerical Modeling [J]. Nature Gas Industry, 2004, 24 (6): 53-56.
[27]
Bérenger JP. A perfectly matched layer for the absorption of electromagnetic waves [J]. Journal of computational physics, 1994, 114 (2): 185. 200.
[28]
Haiqiang Lan, Application of a perfectly matched layer in seismic wave field simulation with an irregular free surface [J]. Geophysical Prospecting, 2016, 64 (1): 112-128.
[29]
Lan, Haiqiang; Zhang, Zhongjie, Comparative study of the free-surface boundary condition in two-dimensional finite-difference elastic wave field simulation [J]. Journal of Geophysics and Engineering, 2011, 8 (2): 275-286.
[30]
Haiqiang Lan and Ling Chen, An upwind fast sweeping scheme for calculating seismic wave first-arrival traveltimes for models with an irregular free surface [J]. Geophysical Prospecting, 2017, DOI: 10.1111/1365-2478.12513.
[31]
Yingjie Gao; Hanjie Song; Jinhai Zhang, Comparison of artificial absorbing boundaries for acoustic wave equation modelling [J]. Exploration Geophysics, 2017, 48 (1) 76: 93.
[32]
Weijuan Meng and Li-Yun Fu, Seismic wavefield simulation by a modified finite element method with a perfectly matched layer absorbing boundary [J]. Journal of Geophysics and Engineering, 2017, 14 (4): 852.
[33]
Michael Brun; Eliass Zafati; Irini Djeran-Maigre; Florent Prunier, Hybrid Asynchronous Perfectly Matched Layer for seismic wave propagation in unbounded domains [J]. Finite Elements in Analysis and Design, 2016, 122: 1-15.
[34]
Michael Brun Jensen; José Alberto Manresa; Ken Steffen Frahm; Analysis of muscle fiber conduction velocity enables reliable detection of surface EMG crosstalk during detection of nociceptive withdrawal reflexes [J]. BMC neuroscience, 2013, 14: 39, doi: 10.1186/1471-2202-14-39.