高精度双吸收边界条件在地震波场模拟中的应用(英文)
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
边界反射处理是地震数值模拟中需要研究的重要问题之一。本文将数学科学计算中提出的双吸收边界条件(DABC)(Hagstrom et al.,2014)应用到地震波场数值模拟中。该方法将一种局部高精度吸收边界(ABC)应用于两条平行的人工边界上,从而实现对边界反射的双吸收。本文以二维声波方程为例,设计一个基于DABC的有限差分(FD)正演方案,给出详细推导过程和实现步骤,相对于完全匹配层(PML),其理论分析和实现难度大大降低,稳定性和灵活性亦得到一定程度的增强。文章最后对典型的均质模型和SEG盐丘模型进行有限差分数值模拟实验,实验结果表明,DABC能够得到远远优于Clayton-Engquist边界条件(CEBC)的吸收效果,并与PML吸收效果相当。
We apply the newly proposed double absorbing boundary condition(DABC)(Hagstrom et al., 2014) to solve the boundary reflection problem in seismic finite-difference(FD) modeling. In the DABC scheme, the local high-order absorbing boundary condition is used on two parallel artificial boundaries, and thus double absorption is achieved. Using the general 2D acoustic wave propagation equations as an example, we use the DABC in seismic FD modeling, and discuss the derivation and implementation steps in detail. Compared with the perfectly matched layer(PML), the complexity decreases, and the stability and fl exibility improve. A homogeneous model and the SEG salt model are selected for numerical experiments. The results show that absorption using the DABC is considerably improved relative to the Clayton–Engquist boundary condition and nearly the same as that in the PML.
引文
Becache,E.,Givoli,D.,and Hagstrom,T.,2010,Highorder absorbing boundary conditions for anisotropic and convective wave equations:Journal of Computational Physics,229(4),1099-1129.
    Berenger,J.P.,1994,A perfectly matched layer for the absorption of electromagnetic waves:Journal of Computational Physics,114,185-200.
    Becache,E.,Fauqueux,S.,and Joly,P.,2003,Stability of perfectly matched layers,group velocities and anisotropic waves:Journal of Computational Physics,188(2),399-433.
    Cerjan,C,Kosloff,D.,Kosloff,R.,and Moshe,R.,1985,A non-reflection boundary condition for discrete acoustic and elastic wave equation:Geophysics,50(4),705-708.
    Clayton,R.,and Engquist,B.,1977,Absorbing boundary conditions for acoustic and elastic wave equations:Bulletin of the Seismological Society of America,67(6),1529-1540.
    Collino,F.,1993,High order absorbing boundary conditions for wave propagation models:Straight line boundary and corner cases:Proceedings of the 2nd International Conference on Mathematical and Numerical Aspects of Wave Propagation,161-171.
    Collino,F.,and Tsogka,C,2001,Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media:Geophysics,66(1),294-307.
    Gedney,S.D.,1996,An anisotropic perfectly matched layer-absorbing medium for the truncation of FDTD lattices:Antennas and Propagation,IEEE Transactions on,44(12),1630-1639.
    Givoli,D.,Neta,B.,2004,High-order non-reflecting boundary scheme for time-dependent waves:Journal of Computational Physics,186(1),24-46.
    Hagstrom,T.,Givoli,D.,Rabinovich,D.,and Bielak,J.,2014,The Double Absorbing Boundary method:Journal of Computational Physics,259,220-241.
    Higdon,R.L.,1987,Numerical absorbing boundary conditions for the wave equation:Mathematics of Computation,49(179),65-90.
    Komatitsch,D.and Tromp,J.,2003,A perfectly matched layer absorbing boundary condition for the secondorder seismic wave equation:Geophysical Journal International,154(1),146-153.
    Liu,Y.,and Sen,M.K.,2010,A hybrid scheme for absorbing edge reflections in numerical modeling of wave propagation:Geophysics,75(2),A1-A6.
    Liu,Y.,and Sen,M.K.,2012,A hybrid absorbing boundary condition for elastic staggered-grid modeling:Geophysical Prospecting,60(6),1114-1132.
    Rabinovich,D.,Givoli,D.,and Becache,E.,2010,Comparison of high-order absorbing boundary conditions and perfectly matched layers in the frequency domain:International Journal for Numerical Methods in Biomedical Engineering,26(10),1351-1369.
    Roden,J.A.,and Gedney S.D.,2000,Convolutional PML(CPML):An efficient FDTD implementation of the CFS-PML for arbitrary media:Microwave and Optical Technology Letters,27(5),334-339.
    Song J.Y,Zheng X.D.,Zhang Y,Xu J.Q.,Qin Z.,and Song X.J.,2011,Frequency domain wave equation forward modeling using gaussian elimination with static pivoting:Applied Geophysics,8(1),60-68.
    Yan,H.Y,and Liu Y,2013a,Visco-acoustic prestack reverse-time migration based on the time-space domain adaptive high-order finite-difference method:Geophysical Prospecting,61,941-954.
    Yan,H.Y,Liu Y,2013b,Acoustic VTI modeling and prestack reverse-time migration based on the time-space domain staggered-grid finite-difference method:Journal of Applied Geophysics,90,41-52.
    Zhao,J.G.,and Shi,R.Q.,2013,Perfectly matched layerabsorbing boundary condition for finite-element timedomain modeling of elastic wave equations:Applied Geophysics,10(3),323-336.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心