Exploration Geophysics Exploration Geophysics Society
Journal of the Australian Society of Exploration Geophysicists
RESEARCH ARTICLE

Hybrid absorbing boundary condition for piecewise smooth curved boundary in 2D acoustic finite difference modelling

Bei Li 1 2 Yang Liu 1 2 3 Yumin Zhao 1 2 Xin Liu 1 2
+ Author Affiliations
- Author Affiliations

1 State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum, Beijing 102249, China.

2 CNPC Key Laboratory of Geophysical Prospecting, China University of Petroleum, Beijing 102249, China.

3 Corresponding author. Email: wliuyang@vip.sina.com

Exploration Geophysics - https://doi.org/10.1071/EG17012
Submitted: 15 January 2017  Accepted: 20 June 2017   Published online: 18 August 2017

Abstract

Flexible computational domains and their corresponding irregular absorbing boundary conditions (ABCs) have previously been shown to improve the efficiency of finite difference (FD) modelling. However, these proposed ABCs for irregular boundaries are based on irregular grid methods. Although they can improve geometric flexibility of FD modelling, irregular grid methods are still complex to implement due to computationally expensive meshing process. To avoid complex mesh generation, FDs in mesh-free discretisation have been developed for nontrivial geometric settings, where scattered nodes can be placed suitably with respect to irregular boundaries and arbitrarily shaped anomalies without coordinate mapping or mesh-element forming. Radial-basis-function-generated FD (RBF-FD) has been proven successful in modelling seismic wave propagation based on mesh-free discretisation. Using RBF-FD, we develop a hybrid ABC for piecewise smooth curved boundary in 2D acoustic wave modelling based on straightforward expanding strategies for sampled boundary and corner nodes. The whole irregular computational domain is combined by an objective zone and a transition zone separated by piecewise smooth curved boundary. Nodal distribution in the objective zone is adaptive to the model structure through special alignment and varying density of nodes, while the transition zone are designed for solving one way wave equation more easily and accurately in the presence of piecewise smooth curved boundary. Modelling examples for homogenous and heterogeneous models with differently shaped computational domains demonstrate the effectiveness of our proposed method.

Key words: acoustic wave equation, hybrid absorbing boundary condition, mesh-free method, numerical modelling, one way wave equation, piecewise smooth curved boundary, radial-basis-function-generated finite difference.


References

Bayona, V., Flyer, N., Lucas, G. M., and Baumgaertner, A. J. G., 2015, A 3-D RBF-FD solver for modeling the atmospheric global electric circuit with topography (GECRBFFD v1.0): Geoscientific Model Development, 8, 3007–3020
A 3-D RBF-FD solver for modeling the atmospheric global electric circuit with topography (GECRBFFD v1.0):CrossRef |

Bayona, V., Flyer, N., Fornberg, B., and Barnett, G. A., 2017, On the role of polynomials in RBF-FD approximations: II. Numerical solution of elliptic PDEs: Journal of Computational Physics, 332, 257–273
On the role of polynomials in RBF-FD approximations: II. Numerical solution of elliptic PDEs:CrossRef |

Bentley, J. L., 1975, Multidimensional binary search trees used for associative searching: Communications of the ACM, 18, 509–517
Multidimensional binary search trees used for associative searching:CrossRef |

Bérenger, J. P., 1994, A perfectly matched layer for the absorption of electromagnetic waves: Journal of Computational Physics, 114, 185–200
A perfectly matched layer for the absorption of electromagnetic waves:CrossRef |

Cerjan, C., Kosloff, D., Kosloff, R., and Resef, M., 1985, A nonreflecting boundary condition for discrete acoustic and elastic wave equations: Geophysics, 50, 705–708
A nonreflecting boundary condition for discrete acoustic and elastic wave equations:CrossRef |

Clayton, R. W., and Engquist, B., 1977, Absorbing boundary conditions for acoustic and elastic wave equations: Bulletin of the Seismological Society of America, 6, 1529–1540

Compani-Tabrizi, B., 1986, k-t scattering formulation of the absorptive acoustic wave equation: wraparound and edge-effect elimination: Geophysics, 51, 2185–2192
k-t scattering formulation of the absorptive acoustic wave equation: wraparound and edge-effect elimination:CrossRef |

Flyer, N., Lehto, E., Blaise, S., Wright, G. B., and St-Cyr, A., 2012, A guide to RBF-generated finite differences for nonlinear transport: shallow water simulations on a sphere: Journal of Computational Physics, 231, 4078–4095
A guide to RBF-generated finite differences for nonlinear transport: shallow water simulations on a sphere:CrossRef |

Flyer, N., Wright, G. B., and Fornberg, B., 2014, Radial basis function-generated finite differences: a mesh-free method for computational geosciences, in W. Freeden, M. Z. Nashed, and T. Sonar, eds., Handbook of Geomathematics: Springer, 1–30.

Flyer, N., Fornberg, B., Bayona, V., and Barnett, G. A., 2016, On the role of polynomials in RBF-FD approximations: I. Interpolation and accuracy: Journal of Computational Physics, 321, 21–38
On the role of polynomials in RBF-FD approximations: I. Interpolation and accuracy:CrossRef |

Fornberg, B., and Flyer, N., 2015a, Solving PDEs with radial basis functions: Acta Numerica, 24, 215–258
Solving PDEs with radial basis functions:CrossRef |

Fornberg, B., and Flyer, N., 2015b, Fast generation of 2-D node distributions for mesh-free PDE discretizations: Computers & Mathematics with Applications, 69, 531–544
Fast generation of 2-D node distributions for mesh-free PDE discretizations:CrossRef |

Fornberg, B., and Piret, C., 2008, A stable algorithm for flat radial basis functions on a sphere: SIAM Journal on Scientific Computing, 30, 60–80
A stable algorithm for flat radial basis functions on a sphere:CrossRef |

Fornberg, B., and Wright, G., 2004, Stable computation of multiquadric interpolants for all values of the shape parameters: Computers & Mathematics with Applications, 48, 853–867
Stable computation of multiquadric interpolants for all values of the shape parameters:CrossRef |

Fries, T., and Matthies, H., 2004, Classification and overview of meshfree methods: Informatikbericht Nr. 2003–3, Institute of Scientific Computing.

Gao, H., and Zhang, J., 2008, Implementation of perfectly matched layers in an arbitrary geometrical boundary for elastic wave modeling: Geophysical Journal International, 174, 1029–1036
Implementation of perfectly matched layers in an arbitrary geometrical boundary for elastic wave modeling:CrossRef |

Gao, Y., Song, H., Zhang, J., and Yao, Z., 2016, Comparison of artificial absorbing boundaries for acoustic wave equation modeling: Exploration Geophysics, 48, 76–93
Comparison of artificial absorbing boundaries for acoustic wave equation modeling:CrossRef |

Gedney, S. D., and Zhao, B., 2010, An auxiliary differential equation formulation for the complex-frequency shifted PML: IEEE Transactions on Antennas and Propagation, 58, 838–847
An auxiliary differential equation formulation for the complex-frequency shifted PML:CrossRef |

Heidari, A. H., and Guddati, M. N., 2006, Highly accurate absorbing boundary conditions for wide-angle wave equations: Geophysics, 71, S85–S97
Highly accurate absorbing boundary conditions for wide-angle wave equations:CrossRef |

Higdon, R. L., 1991, Absorbing boundary conditions for elastic waves: Geophysics, 56, 231–241
Absorbing boundary conditions for elastic waves:CrossRef |

Hu, W., Abubakar, A., and Habashy, T. M., 2007, Application of the nearly perfectly matched layer in acoustic wave modeling: Geophysics, 72, SM169–SM175
Application of the nearly perfectly matched layer in acoustic wave modeling:CrossRef |

Kindelan, M., Moscoso, M., and González-Rodríguez, P., 2016, Radial basis function interpolation in the limit of increasingly flat basis functions: Journal of Computational Physics, 307, 225–242
Radial basis function interpolation in the limit of increasingly flat basis functions:CrossRef |

Komatitsch, D., and Martin, R., 2007, An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation: Geophysics, 72, SM155–SM167
An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation:CrossRef |

Kosloff, R., and Kosloff, D., 1986, Absorbing boundaries for wave propagation problems: Journal of Computational Physics, 63, 363–376
Absorbing boundaries for wave propagation problems:CrossRef |

Lee, G. H., Chung, H. J., and Choi, C. K., 2003, Adaptive crack propagation analysis with the element-free Galerkin method: International Journal for Numerical Methods in Engineering, 56, 331–350
Adaptive crack propagation analysis with the element-free Galerkin method:CrossRef |

Liu, Y., and Sen, M. K., 2009, A new time-space domain high-order finite-difference method for the acoustic wave equation: Journal of Computational Physics, 228, 8779–8806
A new time-space domain high-order finite-difference method for the acoustic wave equation:CrossRef | 1:CAS:528:DC%2BD1MXht1GnsLvL&md5=7ce9f8d8c31ba8af95d6592ae7b10951CAS |

Liu, Y., and Sen, M. K., 2010, A hybrid scheme for absorbing edge reflections in numerical modeling of wave propagation: Geophysics, 75, A1–A6
A hybrid scheme for absorbing edge reflections in numerical modeling of wave propagation:CrossRef |

Liu, Y., and Sen, M. K., 2012, A hybrid absorbing boundary condition for elastic staggered-grid modeling: Geophysical Prospecting, 60, 1114–1132
A hybrid absorbing boundary condition for elastic staggered-grid modeling:CrossRef |

Marcinkovich, C., and Olsen, K., 2003, On the implementation of perfectly matched layers in a three-dimensional fourth-order velocity-stress finite difference scheme: Journal of Geophysical Research, 108, 1978–2012
On the implementation of perfectly matched layers in a three-dimensional fourth-order velocity-stress finite difference scheme:CrossRef |

Martin, R., Komatitsch, D., Gedney, S. D., and Bruthiaux, E., 2010, A high-order time and space formulation of the unsplit perfectly matched layer for the seismic wave equation using Auxiliary Differential Equations (ADE-PML): Computer Modeling in Engineering & Sciences, 56, 17–41

Martin, B., Flyer, N., Fornberg, B., and St-Cyr, A., 2013, Development of mesh-less computational algorithms for seismic exploration: 83rd Annual International Meeting, SEG, Expanded Abstracts, 3543–3547.

Martin, B., Fornberg, B., and St-Cyr, A., 2015, Seismic modeling with radial-basis-function-generated finite differences: Geophysics, 80, T137–T146
Seismic modeling with radial-basis-function-generated finite differences:CrossRef |

Rao, Y., and Wang, Y., 2013, Seismic waveform simulation with pseudo-orthogonal grids for irregular topographic models: Geophysical Journal International, 194, 1778–1788
Seismic waveform simulation with pseudo-orthogonal grids for irregular topographic models:CrossRef |

Ren, Z., and Liu, Y., 2013, A hybrid absorbing boundary condition for frequency-domain finite-difference modeling: Journal of Geophysics and Engineering, 10, 054003
A hybrid absorbing boundary condition for frequency-domain finite-difference modeling:CrossRef |

Reynolds, A. C., 1978, Boundary conditions for the numerical solution of wave propagation problems: Geophysics, 43, 1099–1110
Boundary conditions for the numerical solution of wave propagation problems:CrossRef |

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, 334–339
Convolutional PML (CPML): an efficient FDTD implementation of the CFS-PML for arbitrary media:CrossRef |

Sochacki, J., Kubichek, R., George, J., Fletcher, W., and Smithson, S., 1987, Absorbing boundary conditions and surface waves: Geophysics, 52, 60–71
Absorbing boundary conditions and surface waves:CrossRef |

Takekawa, J., and Mikada, H., 2016, An absorbing boundary condition for acoustic wave propagation using a mesh-free method: Geophysics, 81, T145–T154
An absorbing boundary condition for acoustic wave propagation using a mesh-free method:CrossRef |

Takekawa, J., Mikada, H., and Imamura, N., 2015, A mesh-free method with arbitrary-order accuracy for acoustic wave propagation: Computers & Geosciences, 78, 15–25
A mesh-free method with arbitrary-order accuracy for acoustic wave propagation:CrossRef |

Virieux, J., Calandra, H., and Plessix, R. É., 2011, A review of the spectral, pseudo-spectral, finite-difference and finite-element modelling techniques for geophysical imaging: Geophysical Prospecting, 59, 794–813
A review of the spectral, pseudo-spectral, finite-difference and finite-element modelling techniques for geophysical imaging:CrossRef |

Wenterodt, C., and von Estorff, O., 2009, Dispersion analysis of the meshfree radial point interpolation method for the Helmholtz equation: International Journal for Numerical Methods in Engineering, 77, 1670–1689
Dispersion analysis of the meshfree radial point interpolation method for the Helmholtz equation:CrossRef |

Zeng, Y. Q., He, J. Q., and Liu, Q. H., 2001, The application of the perfectly matched layer in numerical modeling of wave propagation in poroelastic media: Geophysics, 66, 1258–1266
The application of the perfectly matched layer in numerical modeling of wave propagation in poroelastic media:CrossRef |

Zhang, J., and Gao, H., 2011, Irregular perfectly matched layers for 3D elastic wave modeling: Geophysics, 76, T27–T36
Irregular perfectly matched layers for 3D elastic wave modeling:CrossRef |

Zhang, J., and Liu, T., 1999, P-SV-wave propagation in heterogeneous media: grid method: Geophysical Journal International, 136, 431–438
P-SV-wave propagation in heterogeneous media: grid method:CrossRef |

Zhang, J., and Liu, T., 2002, Elastic wave modelling in 3D heterogeneous media: 3D grid method: Geophysical Journal International, 150, 780–799
Elastic wave modelling in 3D heterogeneous media: 3D grid method:CrossRef |

Zhang, W., and Shen, Y., 2010, Unsplit complex frequency-shifted PML implementation using auxiliary differential equations for seismic wave modeling: Geophysics, 75, T141–T154
Unsplit complex frequency-shifted PML implementation using auxiliary differential equations for seismic wave modeling:CrossRef |

Zhang, X., Liu, Y., Cai, X., and Ren, Z., 2015, The AWWE-based hybrid absorbing boundary condition for finite-difference modeling and its application in reverse-time migration: Journal of Applied Geophysics, 123, 93–101
The AWWE-based hybrid absorbing boundary condition for finite-difference modeling and its application in reverse-time migration:CrossRef |



Export Citation