Trapezoid-Grid Finite-Difference Time-Domain Method for 3D Seismic Wavefield Modeling Using CPML Absorbing Boundary Condition

被引:2
作者
Wu, Bangyu [1 ,2 ]
Tan, Wenzhuo [1 ,2 ]
Xu, Wenhao [3 ]
Li, Bo [2 ]
机构
[1] Xi An Jiao Tong Univ, Sch Math & Stat, Xian, Peoples R China
[2] Sinopec Geophys Res Inst, Nanjing, Peoples R China
[3] Xi An Jiao Tong Univ, Sch Informat & Commun Engn, Xian, Peoples R China
关键词
finite difference; trapezoid-grid method; seismic wave simulation; 3D; time-domain method; COMPLEX-FREQUENCY; SYMPLECTIC METHODS; INTERPOLATION; EQUATIONS; INVERSION;
D O I
10.3389/feart.2021.777200
中图分类号
P [天文学、地球科学];
学科分类号
07 ;
摘要
The large computational memory requirement is an important issue in 3D large-scale wave modeling, especially for GPU calculation. Based on the observation that wave propagation velocity tends to gradually increase with depth, we propose a 3D trapezoid-grid finite-difference time-domain (FDTD) method to achieve the reduction of memory usage without a significant increase of computational time or a decrease of modeling accuracy. It adopts the size-increasing trapezoid-grid mesh to fit the increasing trend of seismic wave velocity in depth, which can significantly reduce the oversampling in the high-velocity region. The trapezoid coordinate transformation is used to alleviate the difficulty of processing ununiform grids. We derive the 3D acoustic equation in the new trapezoid coordinate system and adopt the corresponding trapezoid-grid convolutional perfectly matched layer (CPML) absorbing boundary condition to eliminate the artificial boundary reflection. Stability analysis is given to generate stable modeling results. Numerical tests on the 3D homogenous model verify the effectiveness of our method and the trapezoid-grid CPML absorbing boundary condition, while numerical tests on the SEG/EAGE overthrust model indicate that for comparable computational time and accuracy, our method can achieve about 50% reduction on memory usage compared with those on the uniform-grid FDTD method.
引用
收藏
页数:10
相关论文
共 50 条
[41]   A Hybrid Chebyshev Pseudo-Spectral Finite-Difference Time-Domain Method for Numerical Simulation of 2D Acoustic Wave Propagation [J].
Tong, Xiaozhong ;
Sun, Ya .
MATHEMATICS, 2024, 12 (01)
[42]   A time-space domain stereo finite difference method for 3D scalar wave propagation [J].
Chen, Yushu ;
Yang, Guangwen ;
Ma, Xiao ;
He, Conghui ;
Song, Guojie .
COMPUTERS & GEOSCIENCES, 2016, 96 :218-235
[43]   3D electromagnetic modeling of graphitic faults in the Athabasca Basin using a finite-volume time-domain approach with unstructured grids [J].
Lu, Xushan ;
Farquharson, Colin G. ;
Miehe, Jean-Marc ;
Harrison, Grant .
GEOPHYSICS, 2021, 86 (06) :B349-B367
[44]   Time-space domain scalar wave modeling by a novel hybrid staggered-grid finite-difference method with high temporal and spatial accuracies [J].
Zhou, Hongyu ;
Liu, Yang ;
Wang, Jing .
JOURNAL OF COMPUTATIONAL PHYSICS, 2022, 455
[45]   3D Laplace-domain waveform inversion using a low-frequency time-domain modeling algorithm [J].
Ha, Wansoo ;
Kang, Seung-Goo ;
Shin, Changsoo .
GEOPHYSICS, 2015, 80 (01) :R1-R13
[46]   Adaptive finite element for 3D time-domain airborne electromagnetic modeling based on hybrid posterior error estimation [J].
Zhang, Bo ;
Yin, Changchun ;
Ren, Xiuyan ;
Liu, Yunhe ;
Qi, Yanfu .
GEOPHYSICS, 2018, 83 (02) :WB71-WB79
[47]   3D frequency-domain finite-difference viscoelastic-wave modeling using weighted average 27-point operators with optimal coefficients [J].
Gosselin-Cliche, Benjamin ;
Giroux, Bernard .
GEOPHYSICS, 2014, 79 (03) :T169-T188
[48]   Fast 3D time-domain airborne EM forward modeling using random under-sampling [J].
Wang, Haoman ;
Liu, Yunhe ;
Yin, Changchun ;
Ren, Xiuyan ;
Cao, Jin ;
Su, Yang ;
Xiong, Bin .
JOURNAL OF APPLIED GEOPHYSICS, 2021, 191
[49]   Time-Domain Solver for 3D Electromagnetic Problems Using the Method of Moments and the Fast Inverse Laplace Transform [J].
Ohnuki, Shinichiro ;
Kitaoka, Yuya ;
Takeuchi, Takashi .
IEICE TRANSACTIONS ON ELECTRONICS, 2016, E99C (07) :797-800
[50]   Time-space-domain temporal high-order staggered-grid finite-difference schemes by combining orthogonality and pyramid stencils for 3D elastic-wave propagation [J].
Xu, Shigang ;
Liu, Yang ;
Ren, Zhiming ;
Zhou, Hongyu .
GEOPHYSICS, 2019, 84 (04) :T259-T282