Perfectly matched layer on curvilinear grid for the second-order seismic acoustic wave equation

A curvilinear-grid perfectly matched layer (PML) absorbing boundary condition for the second-order seismic acoustic wave equation is presented in this paper. The rectangular grids are transformed into curvilinear grids by using a mathematical mapping to fit the curvilinear boundary, and the original wave equation is reformulated under the curvilinear coordinate system. Based on the reformulated wave equation, theoretical expressions and analysis of the curvilinear-grid PML are given. Furthermore, PML model 1 with symmetric form and PML model 2 with asymmetric form are derived from the same acoustic wave equation. By combination with the finite difference (FD) method, these two models are applied to seismic wave modelling with surface topography. The results show that the absorption effect of these two models discretised by the same second-order time difference and second-order space difference are different, and the symmetric-form PML yields better modelling results than the asymmetric-form.

[1]  M. Piggott,et al.  Three‐dimensional tsunami propagation simulations using an unstructured mesh finite element model , 2013 .

[2]  Bengt Fornberg,et al.  The pseudospectral method; accurate representation of interfaces in elastic wave calculations , 1988 .

[3]  Tsili Wang,et al.  Finite-difference modeling of elastic wave propagation: A nonsplitting perfectly matched layer approach , 2003 .

[4]  Luc Giraud,et al.  New curvilinear scheme for elastic wave propagation in presence of curved topography , 2011 .

[5]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[6]  Zhang Jianfeng,et al.  P–SV-wave propagation in heterogeneous media: grid method , 1999 .

[7]  Bent O. Ruud,et al.  3-D finite-difference elastic wave modeling including surface topography , 1998 .

[8]  Jeroen Tromp,et al.  A perfectly matched layer absorbing boundary condition for the second-order seismic wave equation , 2003 .

[9]  Jingyi Chen.  Application of the nearly perfectly matched layer for seismic wave propagation in 2D homogeneous isotropic media , 2011 .

[10]  D. Komatitsch,et al.  An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation , 2007 .

[11]  René Matzen.  An efficient finite element time‐domain formulation for the elastic second‐order wave equation: A non‐split complex frequency shifted convolutional PML , 2011 .

[12]  J. Bérenger.  Three-Dimensional Perfectly Matched Layer for the Absorption of Electromagnetic Waves , 1996 .

[13]  Hongwei Gao,et al.  Irregular perfectly matched layers for 3D elastic wave modeling , 2011 .

[14]  S. Cummer,et al.  A simple, nearly perfectly matched layer for general electromagnetic media , 2003, IEEE Microwave and Wireless Components Letters.

[15]  Ushnish Basu,et al.  Explicit finite element perfectly matched layer for transient three‐dimensional elastic waves , 2009 .

[16]  Weng Cho Chew,et al.  A 3D perfectly matched medium from modified maxwell's equations with stretched coordinates , 1994 .

[17]  John B. Schneider,et al.  Application of the perfectly matched layer (PML) absorbing boundary condition to elastic wave propagation , 1996 .

[18]  Aria Abubakar,et al.  Application of the Nearly Perfectly Matched Layer In Acoustic Wave Modeling , 2006 .

[19]  Toshifumi Matsuoka,et al.  Acceleration of computation speed for elastic wave simulation using a Graphic Processing Unit , 2011 .

[20]  Qing Huo Liu,et al.  PERFECTLY MATCHED LAYERS FOR ELASTODYNAMICS: A NEW ABSORBING BOUNDARY CONDITION , 1996 .

[21]  David Trubitt,et al.  Physical Modeling , 1991, Advanced Series on Ocean Engineering.

[22]  D. Appelö,et al.  A stable finite difference method for the elastic wave equation on complex geometries with free surfaces , 2007 .

[23]  Hossein Memarian,et al.  イラン,ガッチサラン蒸発形成層の無水石膏,泥灰土,塩の動的特性 , 2013 .

[24]  Hongwei Gao,et al.  Implementation of perfectly matched layers in an arbitrary geometrical boundary for elastic wave modelling , 2008 .

[25]  Gunilla Kreiss,et al.  A new absorbing layer for elastic waves , 2006, J. Comput. Phys..

[26]  Wei Zhang,et al.  Unsplit complex frequency-shifted PML implementation using auxiliary differential equations for seismic wave modeling , 2010 .

[27]  Kim B. Olsen,et al.  On the implementation of perfectly matched layers in a three‐dimensional fourth‐order velocity‐stress finite difference scheme , 2003 .

[28]  Olivier Bou Matar,et al.  Convolutional perfectly matched layer for elastic second-order wave equation. , 2010, The Journal of the Acoustical Society of America.

[29]  J. Vilotte,et al.  The Newmark scheme as velocity–stress time-staggering: an efficient PML implementation for spectral element simulations of elastodynamics , 2005 .

[30]  Jianghai Xia,et al.  Application of the multiaxial perfectly matched layer (M-PML) to near-surface seismic modeling with Rayleigh waves , 2011 .

[31]  Qing Huo Liu,et al.  The application of the perfectly matched layer in numerical modeling of wave propagation in poroelastic media , 2001 .

[32]  Stephen D. Gedney,et al.  Convolution PML (CPML): An efficient FDTD implementation of the CFS–PML for arbitrary media , 2000 .

[33]  Shangxu Wang,et al.  Analysis of complicated structure seismic wave fields , 2010 .

[34]  Ekkehart Tessmer,et al.  3-D elastic modeling with surface topography by a Chebychev spectral method , 1994 .

[35]  C. Tsogka,et al.  Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media , 2001 .