| Peer-Reviewed

Second-order Central Difference Finite-difference Wave Equation Datuming Formulation

Received: 13 June 2022     Accepted: 15 August 2022     Published: 8 October 2022
Views:       Downloads:
Abstract

Several methods of referencing seismic reflection data to a new datum have been researched and presented in the literature, that differ in the approach used to implement the calculation of datumed wavefields. The implementation time determines, in part, the efficiency of such a method. In this paper, a simple finite-difference approach for extrapolating seismic data from one datum to another for seismic near-surface irregularity correction is derived. Kirchhoff’s integral forms the basis on which the formulation is based. The Kirchhoff summation method is used, in a novel way, in agreement with summation theory by considering a single trace as an input wavefield to the algorithm. The approach is considered time-saving in computing the output wavefields as only a single trace serves as an input in any occasion of producing every output trace. For prestack data, the method consists of two steps; calculating the wavefields referenced to the output datum for a source and receivers and time- shifting the resulting data along the mean travel-time ray paths. The algorithm is successfully applied to real data and synthetic shot gathers. For a horizontally layered earth model with a planar surface, mathematical analysis shows that the gradient of the locus of the reflection travel time increases with offset, for a shot-on-end reflection profiling, in conformity with the synthetic data obtained for such a model. In the said analysis, considering parabolic, hyperbolic, and nonhyperbolic travel-time approximations, the gradient forms a useful tool for normal moveout velocity determination. The implementation of the algorithm on field data samples suppressed coherent noise, including the long-wavelength type such as ground roll, in contrast with the static correction method. The normal moveout corrections show that shift along ray-path produced less normal moveout stretch than shift along vertical for near-surface events for the same normal moveout parameters.

Published in Earth Sciences (Volume 11, Issue 5)
DOI 10.11648/j.earth.20221105.16
Page(s) 289-306
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), 2022. Published by Science Publishing Group

Keywords

Datuming, Wave-equation, Central Difference, Finite-difference, Extrapolation, Second-order, Static, NMO Stretch

References
[1] Abedi, M. M., and Riahi, M. A, 2016 Nonhyperbolic stretch-free normal moveout correction. Geophysics 81, NO. 6, P. u87-u95, Doi: 10.1190/GEO2016-0078.1.
[2] Adegoke, K., 2016. Interpreting the summation notation when the lower limit is greater than the upper limit. Downloaded from https://vixra.org/pdf/1601.0207v2.pdf
[3] Alkhalifah, T., and Tsvankin, I., 1995, Velocity analysis for transversely isotropic media: Geophysics, 60, 1550- 1566.
[4] Barison, E.; Brancatelli, G.; Nicolich, R.; Accaino, F.; Giustiniani, M.; Tinivell, U. Wave equation datuming applied to marine OBS data and to land high-resolution seismic profiling Journal of Applied Geophysics 73, 267- 277.
[5] Berryhill, J. R 1979. Wave-equation datuming. Geophysics 44, 1329-1344.
[6] Berryhill, J. R 1984. Wave-equation datuming before stack. Geophysics 49, 2064-2066.
[7] Bevc, D 1997. Flooding the topography: Wave equation datuming of land data with rugged acquisition topography. Geophysics 62, 1558-1569.
[8] Causon, D. M and Mingham, C. G. 2010. Introductory Finite difference Methods for PDEs. Ventus Publishing ApS, ISBN 978-87-7681-642-1.
[9] Claerbout, J. F., and Doherty, S., 1972, Downward continuation of Moveout-corrected seismograms: Geophysics, 37, no. 5, 741-768.
[10] Jimenez-Tejero, C. E., Ranero, C. R., Sallares, V, and Gras, C., 2022. Downward continuation of marine seismic reflection data: an undervalued tool to improve velocity models. Geophys. J. Int. 230, 831-848.
[11] Ke, B., Zhao, B., Liu, C and Chen, B., 1999 Improve image quality by wave-equation datuming based on a single shot gather SEG Technical Program Expanded Abstracts, DOI: 10.1190/1.1845210.
[12] Martini, F. and Bean, C. J, 2002. Application of pre-stack wave equation datuming to remove interface scattering in sub-salt imaging Firstbreak. Vol. 20, No 6, 395-403.
[13] Mcdonald, M. 2012. Numerical Methods in Seismic Wave Propagation. A thesis submitted to the faculty of graduate studies in partial fulfillment of the requirements for the degree of M.sc. in applied mathematics, Department of mathematics and statistics Calgary, Alberta.
[14] Nino, C . L . B and Vides, L . A . M, 2005. Wave equation datuming to correct topography effect on foothill seismic data. Earth Sci. Res. J. Vol. 9, No. 2, 132-138.
[15] Perea Pineda. C. A.; Sanabria Gomez, J. D. and Gonzalez, C. A, 2016 “Static Corrections for Weathering Layer using Wave Equation Datuming and Delay Time Techniques”. Athens: ATINER’S Conference Paper Series, No: ERT2016-1996.
[16] Rasheed, M. 2004. Fifty years of stacking Acta Geophysica vol. 62, no. 3, pp. 505-528, DOI: 10.2478/s11600-013-0191-4.
[17] Schneider, S. A, 1978. Integral formulation for migration in two and three dimensions. Geophysics 43, 49-76.
[18] Shearer, P. M. 2009. Introduction to SEISMOLOGY Second Edition. Cambridge University Press. ISBN-13 978-0-511-58010-9.
[19] Shtivelman, V., 1984. A hybrid method for wavefield computation. Geophysical prospecting 32, 236-257.
[20] Shtivelman, V,1985. Two-dimensionalacousticmodeling by a hybrid method. Geophysics 50, 1273-1284.
[21] Shtivelman, V and Canning, A, 1988. Datum correction by wave-equation extrapolation. Geophysics 53, 1311- 1322.
[22] Tegtmeier-Last, S. 2007. Redatuming of Sparse 3D Seismic Data. THESIS to obtain the degree of doctor at Delft University of Technology. Downloaded from https://repository.tudelft.nl/islandora/object/uuid: 28edc70e-85fe-4351-9c22-e753fe9b0975/datastream/OBJ
[23] Tinivella, U.; Giustiniani, M. and Nicolich, R. 2017 Wave Equation Datuming Applied to S-wave reflection seismic data Journal of Applied Geophysics; https://doi.org/10.1016/j.jappgeo.2018.03.015.
[24] Tinivella, U.; Giustiniani, M. and Vargas-Cordero, I. 2017 Wave Equation Datuming Applied to Seismic Data in Shallow Water Environment and Post-Critical Water Bottom Reflection Energies, 10 (9), 1414; https://doi.org/10.3390/en10091414
[25] Wiggins, J. W, 1984. Kichhoff integral extrapolation and migration of nonplanar data. Geophysics 49, 1239-1248.
[26] Yilmaz, O. and Lucas, D, 1986. Prestack layer replacement. Geophysics 51, 1355-1369.
Cite This Article
  • APA Style

    Daisi Israel Komolafe. (2022). Second-order Central Difference Finite-difference Wave Equation Datuming Formulation. Earth Sciences, 11(5), 289-306. https://doi.org/10.11648/j.earth.20221105.16

    Copy | Download

    ACS Style

    Daisi Israel Komolafe. Second-order Central Difference Finite-difference Wave Equation Datuming Formulation. Earth Sci. 2022, 11(5), 289-306. doi: 10.11648/j.earth.20221105.16

    Copy | Download

    AMA Style

    Daisi Israel Komolafe. Second-order Central Difference Finite-difference Wave Equation Datuming Formulation. Earth Sci. 2022;11(5):289-306. doi: 10.11648/j.earth.20221105.16

    Copy | Download

  • @article{10.11648/j.earth.20221105.16,
      author = {Daisi Israel Komolafe},
      title = {Second-order Central Difference Finite-difference Wave Equation Datuming Formulation},
      journal = {Earth Sciences},
      volume = {11},
      number = {5},
      pages = {289-306},
      doi = {10.11648/j.earth.20221105.16},
      url = {https://doi.org/10.11648/j.earth.20221105.16},
      eprint = {https://article.sciencepublishinggroup.com/pdf/10.11648.j.earth.20221105.16},
      abstract = {Several methods of referencing seismic reflection data to a new datum have been researched and presented in the literature, that differ in the approach used to implement the calculation of datumed wavefields. The implementation time determines, in part, the efficiency of such a method. In this paper, a simple finite-difference approach for extrapolating seismic data from one datum to another for seismic near-surface irregularity correction is derived. Kirchhoff’s integral forms the basis on which the formulation is based. The Kirchhoff summation method is used, in a novel way, in agreement with summation theory by considering a single trace as an input wavefield to the algorithm. The approach is considered time-saving in computing the output wavefields as only a single trace serves as an input in any occasion of producing every output trace. For prestack data, the method consists of two steps; calculating the wavefields referenced to the output datum for a source and receivers and time- shifting the resulting data along the mean travel-time ray paths. The algorithm is successfully applied to real data and synthetic shot gathers. For a horizontally layered earth model with a planar surface, mathematical analysis shows that the gradient of the locus of the reflection travel time increases with offset, for a shot-on-end reflection profiling, in conformity with the synthetic data obtained for such a model. In the said analysis, considering parabolic, hyperbolic, and nonhyperbolic travel-time approximations, the gradient forms a useful tool for normal moveout velocity determination. The implementation of the algorithm on field data samples suppressed coherent noise, including the long-wavelength type such as ground roll, in contrast with the static correction method. The normal moveout corrections show that shift along ray-path produced less normal moveout stretch than shift along vertical for near-surface events for the same normal moveout parameters.},
     year = {2022}
    }
    

    Copy | Download

  • TY  - JOUR
    T1  - Second-order Central Difference Finite-difference Wave Equation Datuming Formulation
    AU  - Daisi Israel Komolafe
    Y1  - 2022/10/08
    PY  - 2022
    N1  - https://doi.org/10.11648/j.earth.20221105.16
    DO  - 10.11648/j.earth.20221105.16
    T2  - Earth Sciences
    JF  - Earth Sciences
    JO  - Earth Sciences
    SP  - 289
    EP  - 306
    PB  - Science Publishing Group
    SN  - 2328-5982
    UR  - https://doi.org/10.11648/j.earth.20221105.16
    AB  - Several methods of referencing seismic reflection data to a new datum have been researched and presented in the literature, that differ in the approach used to implement the calculation of datumed wavefields. The implementation time determines, in part, the efficiency of such a method. In this paper, a simple finite-difference approach for extrapolating seismic data from one datum to another for seismic near-surface irregularity correction is derived. Kirchhoff’s integral forms the basis on which the formulation is based. The Kirchhoff summation method is used, in a novel way, in agreement with summation theory by considering a single trace as an input wavefield to the algorithm. The approach is considered time-saving in computing the output wavefields as only a single trace serves as an input in any occasion of producing every output trace. For prestack data, the method consists of two steps; calculating the wavefields referenced to the output datum for a source and receivers and time- shifting the resulting data along the mean travel-time ray paths. The algorithm is successfully applied to real data and synthetic shot gathers. For a horizontally layered earth model with a planar surface, mathematical analysis shows that the gradient of the locus of the reflection travel time increases with offset, for a shot-on-end reflection profiling, in conformity with the synthetic data obtained for such a model. In the said analysis, considering parabolic, hyperbolic, and nonhyperbolic travel-time approximations, the gradient forms a useful tool for normal moveout velocity determination. The implementation of the algorithm on field data samples suppressed coherent noise, including the long-wavelength type such as ground roll, in contrast with the static correction method. The normal moveout corrections show that shift along ray-path produced less normal moveout stretch than shift along vertical for near-surface events for the same normal moveout parameters.
    VL  - 11
    IS  - 5
    ER  - 

    Copy | Download

Author Information
  • Science Department, Greensprings School, Lagos, Nigeria

  • Sections