• Geophysics and Geophysical Exploration
• /
• v.22 no.2
• /
• pp.56-61
• /
• 2019

In this paper, we studied efficient modeling strategies when we simulate the 3D time-domain acoustic wave propagation using a cell-based finite difference method which can handle the variations of both P-wave velocity and density. The standard finite difference method assigns physical properties such as velocities of elastic waves and density to grid points; on the other hand, the cell-based finite difference method assigns physical properties to cells between grid points. The cell-based finite difference method uses average physical properties of adjacent cells to calculate the finite difference equation centered at a grid point. This feature increases the computational cost of the cell-based finite difference method compared to the standard finite different method. In this study, we used additional memory to mitigate the computational overburden and thus reduced the calculation time by more than 30 %. Furthermore, we were able to enhance the performance of the modeling on several media with limited density variations by using the cell-based and standard finite difference methods together.

• Geophysics and Geophysical Exploration
• /
• v.6 no.2
• /
• pp.77-86
• /
• 2003

We designed a new time-domain, finite-difference, elastic wave modeling technique, based on a displacement formulation. which yields nearly correct solutions to Lamb's problem. Unlike the conventional, displacement-based, finite-difference method using a node-based grid set (where both displacements and material properties such as density and Lame constants are assigned to nodal points), in our new finite-difference method, we use a cell-based grid set (where displacements are still defined at nodal points but material properties within cells). In the case of using the cell-based grid set, stress-free conditions at the free surface are naturally described by the changes in the material properties without any additional free-surface boundary condition. Through numerical tests, we confirmed that the new second-order finite differences formulated in the cell-based grid let generate numerical solutions compatible with analytic solutions unlike the old second-order finite-differences formulated in the node-based grid set.

• Geophysics and Geophysical Exploration
• /
• v.19 no.2
• /
• pp.76-83
• /
• 2016

Migration image can be distorted due to reflected waves in the source and receiver wavefields when discontinuities of input velocity model exist in seismic imaging. To remove reflected waves coming from layer interfaces, it is a common practice to smooth the velocity model for migration. If the velocity model is smoothed, however, the subsurface image can be distorted because the velocity changes around interfaces. In this paper, we attempt to minimize the distortion by reducing reflection energy in the source and receiver wavefields through acoustic impedance homogenization. To make acoustic impedance constant, we define fake density model and use it for migration. When the acoustic impedance is constant over all layers, the reflection coefficient at normal incidence becomes zero and the minimized reflection energy results in the improvement of migration result. To verify our algorithm, we implement the reverse-time migration using cell-based finite-difference method. Through numerical examples, we can note that the migration image is improved at the layer interfaces with high velocity contrast, and it shows the marked improvement particularly in the shallow part.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.51-56
• /
• 2008

The waveform inversion for isotropic media has ever been studied since the 1980s, but there has been few studies for anisotropic media. We present a seismic waveform inversion algorithm for 2-D heterogeneous transversely isotropic structures. A cell-based finite difference algorithm for anisotropic media in time domain is adopted. The steepest descent during the non-linear iterative inversion approach is obtained by backpropagating residual errors using a reverse time migration technique. For scaling the gradient of a misfit function, we use the pseudo Hessian matrix which is assumed to neglect the zero-lag auto-correlation terms of impulse responses in the approximate Hessian matrix of the Gauss-Newton method. We demonstrate the use of these waveform inversion algorithm by applying them to a two layer model and the anisotropic Marmousi model data. With numerical examples, we show that it's difficult to converge to the true model when we assumed that anisotropic media are isotropic. Therefore, it is expected that our waveform inversion algorithm for anisotropic media is adequate to interpret real seismic exploration data.