FREQUENCY DOMAIN FINITE-DIFFERENCE WAVE PROPAGATION MODELING IN 2D ACOUSTIC MEDIA

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.

AMINI NAVID; JAVAHERIAN ABDORAHIM

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

Journal Paper

Paper Information

Journal: IRANIAN JOURNAL OF GEOPHYSICS Year:2010 | Volume:4 | Issue:1 Page(s): 1-16

Download Full-Text

Persian Verion

View:

1,404

Download:

Cites:

Information Journal Paper

Title

FREQUENCY DOMAIN FINITE-DIFFERENCE WAVE PROPAGATION MODELING IN 2D ACOUSTIC MEDIA

Author(s)

AMINI NAVID | JAVAHERIAN ABDORAHIM | Issue Writer Certificate

Keywords

FREQUENCY DOMAINQ3

FINITE-DIFFERENCEQ3

FDFDQ2

SEISMIC MODELINGQ3

PMLQ3

Abstract

Seismic wave propagation modeling is helpful in understanding waveform behavior involving velocity anomalies due to complicated geological structures. On the other hand, forward modeling is the kernel of inversion and tomography procedures, so developing forward modeling algorithms is an inevitable need. In order to model wave propagation, it is necessary to solve complete wave equation for 3D media. Because of limitations in computational resources, solving complete wave equation with all ideal considerations such as heterogeneity, anisotropy, and absorption, for 3D media is difficult. Thus by considering easier conditions we can solve wave equation with normal computers.Solving wave equation with computers needs to discretization techniques such as boundary integral, finite-element or FINITE-DIFFERENCE methods. Boundary integral methods are suitable for simple models while finite-element or FINITE-DIFFERENCE methods are predominantly used for heterogeneous models. Depending on the domain for which wave equation is going to be solved, we can categorize methods to time-space, frequencyspace, Laplace, slowness-space and etc. Recently, the FREQUENCY DOMAIN FINITE-DIFFERENCE (FDFD) method has found extensive application in multi-source experiment modeling, especially in waveform tomography. Because of the special form of wave equation in the FREQUENCY DOMAIN, the modeling of multi-source experiments is a straightforward job. On the other hand, considering absorption mechanisms is easy.This study deals with wave propagation modeling in a 2D acoustic heterogeneous media, using the second-order approximation of acoustic wave equation in the FREQUENCY DOMAIN. The acoustic wave equation is formulated as the first-order hyperbolic system involving fields of pressure and particle velocities. This system is discretized using the second-order staggered grid stencil. To avoid spurious reflections from the model boundaries, a sponge-like perfectly matched layer (PML) is implemented. Solving wave equation in the FREQUENCY DOMAIN leads to a large matrix equation. The key step in the FREQUENCY DOMAIN FINITE-DIFFERENCE modeling that controls computational efficiency is the numerical inversion of the massive matrix. The matrix structure depends on the spatial derivatives approximation. In order to solve this system, different direct solvers can be used. The UMFPACK (unsymmetric multifrontal sparse LU factorization package) solver, which is embedded in MATLAB and has acceptable performance in solving the general system of equations, was selected for this study. For each frequency component, it is necessary to solve a large system of equations to obtain a single frequency component of the pressure wavefield. In this paper we review criteria to avoid numerical dispersion and errors during the FINITE-DIFFERENCE approximation, and time aliasing during frequency sampling. Choosing appropriate spatial and frequency discretization intervals is very important in FDFD modeling. Choosing large D (spatial discretization interval) values will cause the pressure field to be inadequately sampled in space and numerical dispersion. For the scheme presented here, we should have more than 10 grid points per minimum wavelength to keep dispersion errors small. On the other hand according to the sampling theorem, if max df >1/tMAX , where tmax is maximum time and df is frequency components sampling interval, we encounter time aliasing. The present FDFD package is written in MATLAB programming language. MATLAB supports sparse algebra and makes possible the solution of large system of equations with minimum usage of memory, which is the main concern in FDFD algorithms. Examples of waveform modeling using FDFD are shown. In the first example, the headwave originating from a high velocity layer is modeled; in the second example, the wave behavior in a trap model is shown; and in the last example, a salt dome model is studied.

Cites

No record.

References

No record.

Cite

APA: Copy

AMINI, NAVID, & JAVAHERIAN, ABDORAHIM. (2010). FREQUENCY DOMAIN FINITE-DIFFERENCE WAVE PROPAGATION MODELING IN 2D ACOUSTIC MEDIA. IRANIAN JOURNAL OF GEOPHYSICS, 4(1), 1-16. SID. https://sid.ir/paper/133947/en

Vancouver: Copy

AMINI NAVID, JAVAHERIAN ABDORAHIM. FREQUENCY DOMAIN FINITE-DIFFERENCE WAVE PROPAGATION MODELING IN 2D ACOUSTIC MEDIA. IRANIAN JOURNAL OF GEOPHYSICS[Internet]. 2010;4(1):1-16. Available from: https://sid.ir/paper/133947/en

IEEE: Copy

NAVID AMINI, and ABDORAHIM JAVAHERIAN, “FREQUENCY DOMAIN FINITE-DIFFERENCE WAVE PROPAGATION MODELING IN 2D ACOUSTIC MEDIA,” IRANIAN JOURNAL OF GEOPHYSICS, vol. 4, no. 1, pp. 1–16, 2010, [Online]. Available: https://sid.ir/paper/133947/en

Related Journal Papers