• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.34-44
• /
• 2001

Forward modeling by construction of synthetic data is usually practiced in a horizontal surface and a few subsurface structures. However, in-situ surveys often take place in such topographic changes that the corrupted field data always make it difficult to interpret the right signals. To examine the propagation characteristic of elastic waves on the irregular surface, a general mesh generation code for finite element method was modified to consider the topography. By implementing this algorithm, the time domain modeling was practiced in some models with surface topography such as mound, channel, etc. The synthetic data obtained by receivers placed on surface also agreed with the analytic solution. The snapshots showing the total wave-field revealed the propagation characteristic of the elastic waves through complex subsurface structures and helped to identify the signals on the time traces. The transmission of Rayleigh waves along the surface, compressive waves, and sheer waves was observed. Moreover, it turned out that the Rayleigh waves behave like a new source at the edge.

• Economic and Environmental Geology
• /
• v.40 no.6
• /
• pp.739-749
• /
• 2007

The seismic velocity at the formation varies widely with physical properties in the layers. These features on seismic shot gathers are not capable of reproducing normally by numerical modeling of homogeneous medium, so that we need that of random inhomogeneous medium instead. In this study, we conducted Gaussian autocorrelation function (ACF), exponential autocorrelation function and von Karman autocorrelation function for getting inhomogeneous velocity model and applied a simple geological model. According to the results, von Karman autocorrelation function showed short wavelength to the inhomogeneous velocity medium. For numerical modeling for a gas hydrate, we determined a geological model based on field data set gathered in the East sea. The numerical modeling results showed that the von Karman autocorrelation function could properly describe scattering phenomena in the gas hydrate velocity model which contains an inhomogeneous layer. Besides, bottom-simulating-reflectors and scattered waves which appear at seismic shot gather of the field data showed properly in the inhomogeneous numerical modeling.

• Geophysics and Geophysical Exploration
• /
• v.15 no.3
• /
• pp.121-128
• /
• 2012

Absorbing boundary conditions are used to mitigate undesired reflections that can arise at the model's truncation boundaries. We apply a complex frequency shifted perfectly matched layer (CFS-PML) to elastic wave modeling in the frequency domain. Modeling results show that the performance of our implementation is superior to other absorbing boundaries. We consider the coefficients of CFS-PML to be optimal when the kinetic energy becomes to the minimum, and propose the modified CFS-PML that has the CFS-PML coefficient ${\alpha}_{max}$ defined as a function of frequency. Results with CFS-PML and modified CFS-PML are significantly improved compared with those of the classical PML technique suffering from large spurious reflections at grazing incidence.

• Geophysics and Geophysical Exploration
• /
• v.6 no.1
• /
• pp.13-22
• /
• 2003

TSP (Tunnel Seismic Profiling) survey is a technique for imaging and characterizing geological structures ahead of a tunnel face. The seismic modeling algorithm and the synthetic data could be helpful for TSP surveys. However, there is few algorithm to describe the propagation of the elastic waves around the tunnel. In this study, existing 2-dimensional seismic modeling algorithm using finite element method was modified to make a suitable algorithm for TSP modeling. Using this algorithm, TSP modeling was practiced in some models. And the synthetic data was analyzed to examine the propagation characteristics of the elastic waves. First of all, the modeling for the homogeneous tunnel model was practiced to examine the propagation characteristics of the direct waves in the vicinity of the tunnel. And the algorithm was applied to some models having reflector which is perpendicular or parallel to the excavation direction. From these, the propagation characteristics of the reflected waves were examined. Furthermore, two source-receiver arrays were used in respective models to investigate the properties of the two arrays. These modeling algorithm and synthetic data could be helpful in interpreting TSP survey data, developing inversion algorithm and designing new source-receiver arrays.

• Geophysics and Geophysical Exploration
• /
• v.5 no.2
• /
• pp.71-77
• /
• 2002

We developed a modeling and visualization software that can analyze 3-dimensional seismic data. The software divides 3 dimensional space into a series of vertical and horizontal polygons, and allows the various seismic attributes and other spatial information to be stored on these polygons. The program can pick a particular pattern in semi-automatic mode, and store the pattern in the spatial DB. The pattern can be modeled and visualized in 3 dimensional space.

• Economic and Environmental Geology
• /
• v.49 no.3
• /
• pp.225-242
• /
• 2016

Rock physics serves as a useful tool for seismic reservoir characterization and monitoring by providing quantitative relationships between rock properties and seismic data. Rock physics models can predict effective moduli for reservoirs with different mineral components and pore fluids from well-log data. The distribution of reservoirs and fluids for the entire seismic volume can also be estimated from rock physics models. The first part of this report discusses the Voigt, Reuss, and Hashin-Shtrikman bounds for effective elastic moduli and the Gassmann fluid substitution. The second part reviews various contact models for moderate- to high-porosity sands. In the third part, constant-cement model, known to work well for the sand that gradually loses porosity with deteriorating sorting, was applied to the well-log data from an oil field in the North Sea. Lastly, the rock physics template constructed from the constant-cement model and the results from the prestack inversion of 2D seismic data were combined to predict the lithology and fluid types for the sand reservoir of this oil field.

• Geophysics and Geophysical Exploration
• /
• v.9 no.2
• /
• pp.139-147
• /
• 2006

We had conducted a numerical modeling to investigate seismic properties of gas hydrate with field parameters acquired over the East sea in 1998. We used a 2-D staggered grid finite difference method to generate synthetic elastic seismograms for multi-channel seismic survey, OBC (Ocean Bottom Cable) survey and VCS (Vertical Cable Seismic) survey. The results of this study showed that the method using staggered grid yielded stable results and could be used to seismic imaging. We could find out the high amplitude anomaly and the phase reversal phenomenon of reflection wave at interface between the gas hydrate layer and free gas layer such a BSR (Bottom Simulating Reflector) which is the evidence for existence of gas hydrate in seismic reflection data. And we computed the reflection coefficients at the incident angles corresponding to offset distance with the synthetic seismograms. The reflection coefficients acquired from the numerical modeling were nearly consistent with the reflection coefficient computed by Shuey's equation.

• Geophysics and Geophysical Exploration
• /
• v.13 no.1
• /
• pp.1-8
• /
• 2010

Seismic reflection surveying is one of the most widely used and effective techniques for coal seam structure delineation and risk mitigation for underground longwall mining. However, the ability of the method can be compromised by the presence of volcanic cover. This problem arises within parts of the Bowen and Sydney Basins of Australia and seismic surveying can be unsuccessful. As a consequence, such areas are less attractive for coal mining. Techniques to improve the success of seismic surveying over basalt flows are needed. In this paper, we use elastic wave-equation-based forward modelling techniques to investigate the effects and characteristics of seismic wave propagation under different settings involving changes in basalt properties, its thickness, lateral extent, relative position to the shot position and various forms of inhomogeneity. The modelling results suggests that: 1) basalts with high impedance contrasts and multiple flows generate strong multiples and weak reflectors; 2) thin basalts have less effect than thick basalts; 3) partial basalt cover has less effect than full basalt cover; 4) low frequency seismic waves (especially at large offsets) have better penetration through the basalt than high frequency waves; and 5) the deeper the coal seams are below basalts of limited extent, the less influence the basalts will have on the wave propagation. In addition to providing insights into the issues that arise when seismic surveying under basalts, these observations suggest that careful management of seismic noise and the acquisition of long-offset seismic data with low-frequency geophones have the potential to improve the seismic results.

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.93-98
• /
• 2008

Abstract: Surface topography has a significant influence on seismic wave propagation in a reflection seismic exploration. Effects of surface topography on two-dimensional elastic wave propagation are investigated through modeling using a weighted-averaging (WA) finite-element method (FEM), which is computationally more efficient than conventional FEM. Effects of air layer on wave propagation are also investigated using flat surface models with and without air. To validate our scheme in modeling including topography, we compare WA FEM results for irregular topographic models against those derived from conventional FEM using one set of rectangular elements. For the irregular surface topography models, elastic wave propagation is simulated to show that breaks in slope act as a new source for diffracted waves, and that Rayleigh waves are more seriously distorted by surface topography than P-waves.

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.153-160
• /
• 2008

Seismic modeling is used to simulate wave propagation in the earth. Although the earth's subsurface is usually semi-infinite, we cannot handle the semi-infinite model in seismic modeling because of limited computational resources. For this reason, we usually assume a finite-sized model in seismic modeling. In that case, we need to eliminate the edge reflections arising from the artificial boundaries introducing a proper boundary condition. In this study, we changed three kinds of boundary conditions (sponge boundary condition, Clayton and Engquist's absorbing boundary condition, and Higdon's transparent boundary condition) so that they can be applied in elastic wave modeling for anisotropic media. We then apply them to several models whose Poisson's ratios are different. Clayton and Engquist's absorbing boundary condition is unstable in both isotropic and anisotropic media, when Poisson's ratio is large. This indicates that the absorbing boundary condition can be applied in anisotropic media restrictively. Although the sponge boundary condition yields good results for both isotropic and anisotropic media, it requires too much computational memory and time. On the other hand, Higdon's transparent boundary condition is not only inexpensive, but also reduce reflections over a wide range of incident angles. We think that Higdon's transparent boundary condition can be a method of choice for anisotropic media, where Poisson's ratio is large.