Applied and Computational Mathematics

| Peer-Reviewed |

Effect of Grid Step Sizes on Computational Time Using Finite-Difference Method for Seismic Wave Modeling

Received: 05 September 2015    Accepted: 21 September 2015    Published: 15 April 2016
Views:       Downloads:

Share This Article

Abstract

Finite-difference method, a popular seismic forward modeling method is a technique which allows us to numerically solve partial differential equations like the wave equation solved in this paper. Beyond its use in standard data acquisition, it is a very instructive tool to understand how waves propagate in the earth's subsurface. Since the accuracy obtainable by using the finite difference scheme lies solely on its stability and ability to handle grid dispersion, this is achievable by applying appropriate grid step sizes. The developed finite-difference method was employed to generate snapshots of synthetic seismograms to highlight the effect of grid step sizes on computational time while ensuring numerical stability of the scheme used through accurate discretization of the medium and adopting Perfectly Matched Layer (PML) absorbing boundary conditions. Results shows that for a grid size of 5m x 5m x 5m having 260 x 260 x 100 grid points and time step of 100 - 500, the wavefield propagating is horizontally symmetric. From the results, the importance of grid step sizes on computational time is re-emphasized.

DOI 10.11648/j.acm.20160502.14
Published in Applied and Computational Mathematics (Volume 5, Issue 2, April 2016)
Page(s) 56-63
Creative Commons

This is an Open Access article, distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution and reproduction in any medium or format, provided the original work is properly cited.

Copyright

Copyright © The Author(s), 2024. Published by Science Publishing Group

Keywords

Finite Difference Method, Forward Modeling, Synthetic Seismogram, Stability, Wavefield, Perfectly Matched Layer

References
[1] Bohlen, T. & E. H. Saenger., 2006: Accuracy of heterogeneous staggered-grid finite-difference modeling of Rayleigh waves, Geophysics, 71, No. 4, Pages T109-T115.
[2] Collinos, F. and Tsogka, C., 2001, Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media. Geophysics 66, No 1, 294 – 307.
[3] Higham N. J., 1996. Accuracy and stability of numerical algorithms, Society for Industrial and applied Mathematics, Philadelphia, ISBN 0-8987 - 1355-2.
[4] Kriiger, O. S., Saenger E. H., Oates, S. J and Shapiro, S. A., 2007, A numerical study on reflection coefficients of fractured media. Geophysics, 72, No. 3. D61 – D67.
[5] Kurzmann, A. Przebindowska, A, Köhn, D. and Bohlen, T., 2013: Acoustic full waveform tomography in the presence of attenuation: a sensitivity analysis, Geophysical Journal International, 195 (2), 985 – 1000.
[6] Moczo, P., Jozef, K and Halada, L., (2004), The finite-difference method for seismologist.
[7] Moczo, P., J. Kristek, M. Galis., P. Pazak., M. Balazovjech., 2007, The Finite-Difference and Finite-Element Modeling of Seismic Wave Propagation and Earthquake motion. Acta Physica Slovaca, 57, No. 2, 177-406.
[8] Olowofela, J. A., Akinyemi, O. D., and Ajani, O. O., (2015): Stability analysis for finite-difference schemeused for seismic imaging using amplitude and phase portrait. Applied and Computational Mathematics 4(1): 1-6.
[9] Saenger, E. & T. Bohlen., 2004: Finite-difference modeling of viscoelastic and anisotropic wave propagation using the rotated staggered grid, Geophysics, 69, No. 2, 583-591.
[10] Virieux, J., 1986, P – SV wave propagation in heterogeneous media: Velocity – Stress finite-difference method: Geophysics, 51, 889-901.
Author Information
  • Department of Physics, Federal University of Agriculture (FUNAAB), Abeokuta, Nigeria

  • Department of Physics, Federal University of Agriculture (FUNAAB), Abeokuta, Nigeria

  • Department of Physics and Solar Energy, Bowen University, Iwo, Nigeria

Cite This Article
  • APA Style

    Olowofela Joseph. A., Akinyemi Olukayode. D., Ajani Olumide. O. (2016). Effect of Grid Step Sizes on Computational Time Using Finite-Difference Method for Seismic Wave Modeling. Applied and Computational Mathematics, 5(2), 56-63. https://doi.org/10.11648/j.acm.20160502.14

    Copy | Download

    ACS Style

    Olowofela Joseph. A.; Akinyemi Olukayode. D.; Ajani Olumide. O. Effect of Grid Step Sizes on Computational Time Using Finite-Difference Method for Seismic Wave Modeling. Appl. Comput. Math. 2016, 5(2), 56-63. doi: 10.11648/j.acm.20160502.14

    Copy | Download

    AMA Style

    Olowofela Joseph. A., Akinyemi Olukayode. D., Ajani Olumide. O. Effect of Grid Step Sizes on Computational Time Using Finite-Difference Method for Seismic Wave Modeling. Appl Comput Math. 2016;5(2):56-63. doi: 10.11648/j.acm.20160502.14

    Copy | Download

  • @article{10.11648/j.acm.20160502.14,
      author = {Olowofela Joseph. A. and Akinyemi Olukayode. D. and Ajani Olumide. O.},
      title = {Effect of Grid Step Sizes on Computational Time Using Finite-Difference Method for Seismic Wave Modeling},
      journal = {Applied and Computational Mathematics},
      volume = {5},
      number = {2},
      pages = {56-63},
      doi = {10.11648/j.acm.20160502.14},
      url = {https://doi.org/10.11648/j.acm.20160502.14},
      eprint = {https://download.sciencepg.com/pdf/10.11648.j.acm.20160502.14},
      abstract = {Finite-difference method, a popular seismic forward modeling method is a technique which allows us to numerically solve partial differential equations like the wave equation solved in this paper. Beyond its use in standard data acquisition, it is a very instructive tool to understand how waves propagate in the earth's subsurface. Since the accuracy obtainable by using the finite difference scheme lies solely on its stability and ability to handle grid dispersion, this is achievable by applying appropriate grid step sizes. The developed finite-difference method was employed to generate snapshots of synthetic seismograms to highlight the effect of grid step sizes on computational time while ensuring numerical stability of the scheme used through accurate discretization of the medium and adopting Perfectly Matched Layer (PML) absorbing boundary conditions. Results shows that for a grid size of 5m x 5m x 5m having 260 x 260 x 100 grid points and time step of 100 - 500, the wavefield propagating is horizontally symmetric. From the results, the importance of grid step sizes on computational time is re-emphasized.},
     year = {2016}
    }
    

    Copy | Download

  • TY  - JOUR
    T1  - Effect of Grid Step Sizes on Computational Time Using Finite-Difference Method for Seismic Wave Modeling
    AU  - Olowofela Joseph. A.
    AU  - Akinyemi Olukayode. D.
    AU  - Ajani Olumide. O.
    Y1  - 2016/04/15
    PY  - 2016
    N1  - https://doi.org/10.11648/j.acm.20160502.14
    DO  - 10.11648/j.acm.20160502.14
    T2  - Applied and Computational Mathematics
    JF  - Applied and Computational Mathematics
    JO  - Applied and Computational Mathematics
    SP  - 56
    EP  - 63
    PB  - Science Publishing Group
    SN  - 2328-5613
    UR  - https://doi.org/10.11648/j.acm.20160502.14
    AB  - Finite-difference method, a popular seismic forward modeling method is a technique which allows us to numerically solve partial differential equations like the wave equation solved in this paper. Beyond its use in standard data acquisition, it is a very instructive tool to understand how waves propagate in the earth's subsurface. Since the accuracy obtainable by using the finite difference scheme lies solely on its stability and ability to handle grid dispersion, this is achievable by applying appropriate grid step sizes. The developed finite-difference method was employed to generate snapshots of synthetic seismograms to highlight the effect of grid step sizes on computational time while ensuring numerical stability of the scheme used through accurate discretization of the medium and adopting Perfectly Matched Layer (PML) absorbing boundary conditions. Results shows that for a grid size of 5m x 5m x 5m having 260 x 260 x 100 grid points and time step of 100 - 500, the wavefield propagating is horizontally symmetric. From the results, the importance of grid step sizes on computational time is re-emphasized.
    VL  - 5
    IS  - 2
    ER  - 

    Copy | Download

  • Sections