A NEW NUMERICAL AND ANALYTICAL SCHEME TO SOLVE THE FULL WAVE EQUATION FOR SEISMIC MODELING BASED ON REM AND LEAPFROG METHODS

Q: How can I download an article?

To download an article from SID, first log in to the site, search for the article title, and click on the 'Download Article' option.

Q: How can I download an ISI article?

To download an ISI article on SID, enter the keyword or article title in the search bar, view the relevant results, click on the desired article, and select the 'Download Article' option.

Q: How can I access the SID database?

To access the SID database, visit SID.ir, create an account, and log in to access scientific resources.

Q: Is downloading articles from SID free?

Some articles on SID are available for free, while others require payment. Details are specified on the article's page.

MORADPOURI F.; MORADZADEH A.; PESTANA R.C.; SOLEIMANI MONFARED M.

مرکز اطلاعات علمی Scientific Information Database (SID) - Trusted Source for Research and Academic Resources

Journal Paper

Paper Information

Journal: JOURNAL OF ANALYTICAL AND NUMERICAL METHODS IN MINING ENGINEERING Year:2016 | Volume:5 | Issue:10 Page(s): 41-48

Download Full-Text

Persian Verion

View:

819

Download:

Cites:

Information Journal Paper

Title

A NEW NUMERICAL AND ANALYTICAL SCHEME TO SOLVE THE FULL WAVE EQUATION FOR SEISMIC MODELING BASED ON REM AND LEAPFROG METHODS

Author(s)

MORADPOURI F. | MORADZADEH A. | PESTANA R.C. | SOLEIMANI MONFARED M. | Issue Writer Certificate

Keywords

NUMERICAL METHODS

FULL WAVE EQUATION

WAVE-FIELD EXTRAPOLATION

SYMPLECTIC INTEGRATOR

LEAPFROG-RAPID EXPANSION METHOD

Abstract

Summary: The aim of this paper is to present a new numerical method to solve the wave equation with a good accuracy and high stability using the Leapfrog SYMPLECTIC INTEGRATOR and rapid expansion method (REM). It can be used for seismic modeling and reverse time migration (RTM). Using the REM with Fourier transform method for spatial derivative, the LEAPFROG-RAPID EXPANSION METHOD (L-REM) can be even used for larger time steps. The L-REM provides the solution of the wave equation and its first time derivative at the current time step. In addition to the very low error for the small time steps, increasing the time step also lead to more accurate results and high stability in comparison with the similar methods such as Stormer-Verlet, Leapfrog and Stormer-Verlet-rapid expansion method which will also be discussed in this paper. Introduction: WAVE-FIELD EXTRAPOLATION is implemented by solving the wave equation through various mathematical methods. The finite difference method is a well-known and popular numerical tool to discretize the wave equation, and its use has been common in the approximation of the spatial and time derivatives for a wave-field. Originally, the time operator was approximated by a second-order scheme, whereas the spatial derivatives were approximated by a fourth-order scheme. Approximating the time derivative in this way may introduce numerical error, leading to distortion of the pulse and numerical dispersion, which can be avoided with small time steps at the expense of increasing the computational time. Also some other methods such as Stormer-Verlet (SV), Leapfrog (L) and Stormer-Verlet-rapid expansion method (SV-REM) have been presented to improve the solution of wave equation. In the current study, a symplectic scheme based on Leapfrog integrator and rapid expansion method (L-REM) is proposed to extrapolate the wave-field and its first derivatives in time for the same time step which can be used to calculate Poynting vectors for wave-field separation and to calculate the reflection angles. Methodology and Approaches: In order to verify the numerical accuracy and behaviour of the error associated with Leapfrog-REM scheme a numerical example has been presented to be solved using different time sampling values. For implementation, an explosive source used in the centre of the computational domain having a Ricker wavelet with a maximum frequency of 25 Hz. Results and Conclusions: The presented L-REM scheme provides the solution of the wave equation and its first time derivative for different time steps. In addition to the very low error for the small time steps, increasing the time step also lead in the more accurate results and high stability in comparison with the similar methods.

Cites

No record.

References

No record.

Cite

APA: Copy

MORADPOURI, F., MORADZADEH, A., PESTANA, R.C., & SOLEIMANI MONFARED, M.. (2016). A NEW NUMERICAL AND ANALYTICAL SCHEME TO SOLVE THE FULL WAVE EQUATION FOR SEISMIC MODELING BASED ON REM AND LEAPFROG METHODS. JOURNAL OF ANALYTICAL AND NUMERICAL METHODS IN MINING ENGINEERING, 5(10 ), 41-48. SID. https://sid.ir/paper/266056/en

Vancouver: Copy

MORADPOURI F., MORADZADEH A., PESTANA R.C., SOLEIMANI MONFARED M.. A NEW NUMERICAL AND ANALYTICAL SCHEME TO SOLVE THE FULL WAVE EQUATION FOR SEISMIC MODELING BASED ON REM AND LEAPFROG METHODS. JOURNAL OF ANALYTICAL AND NUMERICAL METHODS IN MINING ENGINEERING[Internet]. 2016;5(10 ):41-48. Available from: https://sid.ir/paper/266056/en

IEEE: Copy

F. MORADPOURI, A. MORADZADEH, R.C. PESTANA, and M. SOLEIMANI MONFARED, “A NEW NUMERICAL AND ANALYTICAL SCHEME TO SOLVE THE FULL WAVE EQUATION FOR SEISMIC MODELING BASED ON REM AND LEAPFROG METHODS,” JOURNAL OF ANALYTICAL AND NUMERICAL METHODS IN MINING ENGINEERING, vol. 5, no. 10 , pp. 41–48, 2016, [Online]. Available: https://sid.ir/paper/266056/en

Related Journal Papers