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.
| Published in | Applied and Computational Mathematics (Volume 5, Issue 2) |
| DOI | 10.11648/j.acm.20160502.14 |
| 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), 2016. Published by Science Publishing Group |
Finite Difference Method, Forward Modeling, Synthetic Seismogram, Stability, Wavefield, Perfectly Matched Layer
| [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. |
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
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
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
Share This Article
Twitter
Linked In
Facebook
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://article.sciencepublishinggroup.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}
}
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 -
Email: