The finite difference modelling of wave equation are generally performed in Cartesian coordinates which is relatively easy to implement, but when dealing 3D axisymmetric media, even for simple layered medium, the simulation can be very time-consuming. Fortunately, for 3D (approximate) axisymmetric model, a more economical approach is to solve the wave equation in cylindrical coordinates, which can reduce the computation time to nearly as short as that of the 2D calculations. In this paper, the unsplit convolutional PML absorbing boundary is first implemented in 3D acoustic wave equation modelling in cylindrical coordinates. Numerical simulation results show that it’s an accurate and efficient algorithm and exhibits good performance for absorbing artificial reflection. Thus our algorithm can be a tantalizing choice to efficiently simulate axisymmetric models.


Article metrics loading...

Loading full text...

Full text loading...


  1. Dong, C.H., Wang, S.X., Zhao, J.G. and Tang, G.Y.
    [2013] Numerical experiment and analysis of the differential acoustic resonance spectroscopy for elastic property measurements. Journal of Geophysics and Engineering, 10, 054002.
    [Google Scholar]
  2. Komatitsch, D. and Martin, R.
    [2007] An unsplit convolutional perfectly matched layer improved at grazing incidence forthe seismic wave equation. Geophysics, 72(5), SM155–SM167.
    [Google Scholar]
  3. Takenaka, H. and Tanaka, H.
    [2003] Quasi-cylindrical 2.5D wave modeling for large-scale seismic surveys. Geophysical Research Letters, 30(21), 2086, doi:10.1029/2003GL018068.
    https://doi.org/10.1029/2003GL018068 [Google Scholar]

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error