• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2000.04a
• /
• pp.55-62
• /
• 2000

In constructing a synthetic seismogram using ray-tracing method a suite of ray-code is required to obtain a realistic seismogram which is similar to the actual seismogram or earthquake record under consideration. An infinite number of rays exist for any arbitrarily located source and receiver. One select only a finite number of such rays in computing a synthetic seismogram so their selection becomes important to the validity of the seismogram being generated. Missing certain important rays or an inappropriate selection of ray-codes in tracing rays may result in wrong interpretation of the earthquake record or seismogram. Automatic ray-code generation will eliminate such problems. In this study we have developed an efficient algorithm, with which one can generate systemastically all the ray-codes connecting source and receiver arbitrarily located. The result of this work will helpful in analysing multiple reflections in seismic data processing as well as simulating Lg wave and multiply reflected or converted phases in earthquake study.

• 한국신재생에너지학회:학술대회논문집
• /
• 2005.06a
• /
• pp.643-646
• /
• 2005

The gas hydrate exploration using seismic reflection data, the detection of BSR(Bottom Simulating Reflector) on the seismic section is the most important work flow because the BSR have been interpreted as being formed at the base of a gas hydrate zone. Usually, BSR has some dominant qualitative characteristics on seismic section i.e. Wavelet phase reversal compare to sea bottom signal, Parallel layer with sea bottom, Strong amplitude, Masking phenomenon above the BSR, Cross bedding with other geological layer. Even though a BSR can be selected on seismic section with these guidance, it is not enough to conform as being true BSR. Some other available methods for verifying the BSR with reliable analysis quantitatively i.e. Interval velocity analysis, AVO(Amplitude Variation with Offset)analysis etc. Usually, AVO analysis can be divided by three main parts. The first part is AVO analysis, the second is AVO modeling and the last is AVO inversion. AVO analysis is unique method for detecting the free gas zone on seismic section directly. Therefore it can be a kind of useful analysis method for discriminating true BSR, which might arise from an Possion ratio contrast between high velocity layer, partially hydrated sediment and low velocity layer, water saturated gas sediment. During the AVO interpretation, as the AVO response can be changed depend upon the water saturation ratio, it is confused to discriminate the AVO response of gas layer from dry layer. In that case, the AVO modeling is necessary to generate synthetic seismogram comparing with real data. It can be available to make conclusions from correspondence or lack of correspondence between the two seismograms. AVO inversion process is the method for driving a geological model by iterative operation that the result ing synthetic seismogram matches to real data seismogram wi thin some tolerance level. AVO inversion is a topic of current research and for now there is no general consensus on how the process should be done or even whether is valid for standard seismic data. Unfortunately, there are no well log data acquired from gas hydrate exploration area in Korea. Instead of that data, well log data and seismic data acquired from gas sand area located nearby the gas hydrate exploration area is used to AVO analysis, As the results of AVO modeling, type III AVO anomaly confirmed on the gas sand layer. The Castagna's equation constant value for estimating the S-wave velocity are evaluated as A=0.86190, B=-3845.14431 respectively and water saturation ratio is $50\%$. To calculate the reflection coefficient of synthetic seismogram, the Zoeppritz equation is used. For AVO inversion process, the dataset provided by Hampson-Rushell CO. is used.

• Economic and Environmental Geology
• /
• v.29 no.5
• /
• pp.629-637
• /
• 1996

In this paper a numerical experiment is conducted to determine the low acoustic impedance of a thin oil or gas reservoir from a seismogram by using the generalized linear inversion method. The seismograms used are normal incident synthetic seismograms containing p-wave primary reflections, multiples, and peg-leg multiples on the layers consisting of oil-, gas-, water-filled sandstone incased in shales. In this experiment the acoustic impedance, the location of reservoir boundary, thickness, and source wavelet are assumed initially and revised iteratively by the least-squares-error technique until the difference between the seismogram and calculated one is very small. This experiment shows that the acoustic impedance and thickness, about 10 m thick, can be determined by the inversion.

• Economic and Environmental Geology
• /
• v.23 no.3
• /
• pp.353-358
• /
• 1990

This paper presents numerical experiments on migration of synthetic seismograms using by $45^{\circ}$ wave equation. The seismograms used are zero-offset seismogram (corresponding to stacked section) on point reflectors, dipping plane reflector, faulted and folded layers. The seismograms are constructed by upward continuation of seismic source wavelets, exploading on subsurface reflection interfaces, to the earth surface. The synthetic seismograms are migrated by downward continuation and imaging. The upward and downward continuations are implemented by solving the $45^{\circ}$ wave equation with the finite-difference method. Migration of the synthetic data used in this study results in relatively accurate reposition of subsurface structures while the synthetic sections are quite different from the structures.

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 1995.03a
• /
• pp.108-108
• /
• 1995

There are always some difficulties to discriminate artificial exlposions from micro-earthquakes, furthermore more difficulties to identify and determine decoupled explosions and/or multiple explosions from micro-earthquakes. In this study we use the synthetic seismogram of the in homogeneous models between the source and the observation station in order to find the source effect of the geological environment. We have found some source characteristics of the air-filled and/or water-filled cavity that we can hardly see P-n and S- waves arrivals and that the high frequency coda waves are well observed compared to the coupled explosions or earthquakes.

• Economic and Environmental Geology
• /
• v.28 no.5
• /
• pp.481-485
• /
• 1995

This paper presents a new approaching method to determine the velocity and geometry of shallow subsurface from seismic refraction events. After picking the first breaks from seismic refraction data, we assume that field refraction seismogram can be replaced by the unit delta function having time shift of first break. Time curves are generated by shooting ray tracing. The partial derivatives seismogram for a damped least squares method is computed analytically at each step of the forward ray tracing. The technique is successfully tested on synthetic and real data. It has the advantage of real full waveform inversion, which is robust at low frequency band even if the initial guess is far from the true model.

• The Journal of Engineering Geology
• /
• v.23 no.4
• /
• pp.457-465
• /
• 2013

We integrated and correlated datasets from surface and subsurface geophysics, drilling cores, and engineering geology to identify geological interfaces and characterize the joints and fracture zones within the rock mass. The regional geometry of a geologically weak zone was investigated via a fence projection of electrical resistivity data and a borehole image-processing system. Subsurface discontinuities and intensive fracture zones within the rock mass are delineated by cross-hole seismic tomography and analyses of dip directions in rose diagrams. The dynamic elastic modulus is studied in terms of the P-wave velocity and Poisson's ratio. Subsurface discontinuities, which are conventionally identified using the N value and from core samples, can now be identified from anomalous reflection coefficients (i.e., acoustic impedance contrast) calculated using a pair of well logs, comprising seismic velocity from suspension-PS logging and density from logging. Intensive fracture zones identified in the synthetic seismogram are matched to core loss zones in the drilling core data and to a high concentration of joints in the borehole imaging system. The upper boundaries of fracture zones are correlated to strongly negative amplitude in the synthetic trace, which is constructed by convolution of the optimal Ricker wavelet with a reflection coefficient. The standard deviations of dynamic elastic moduli are higher for fracture zones than for acompact rock mass, due to the wide range of velocities resulting from the large numbers of joints and fractures within the zone.

• Geophysics and Geophysical Exploration
• /
• v.11 no.4
• /
• pp.361-367
• /
• 2008

When constructing a synthetic seismogram in the earthquake study or in seismic data interpretation by using a ray-tracing technique, the most troublesome and error-prone task is to define a suite of ray codes for the corresponding rays to trace in advance. An infinite number of rays exist for any arbitrarily located source and receiver in a medium. Missing certain important rays or an inappropriate selection of ray codes in tracing rays may result in wrong interpretation of the earthquake record or seismogram. Automatic ray code generation could be able to eliminate those problems. In this study we have developed an efficient algorithm with which one can generate systematically all the ray codes for the source(s) and receiver(s) arbitrarily located in a model. The result of this work could be used not only in analysing multiples in seismic data processing and interpretation, but also in coda wave study, study on the amplification effects in a basin and phase identification of the waves multiply reflected/refracted in earthquake study.

• 한국석유지질학회:학술대회논문집
• /
• autumn
• /
• pp.2-16
• /
• 1999

This study includes structural analysis of the northern Pattani Basin, areal description of depositional facies, and their spatial relationships using 3-D seismic and well data. Well log data indicate that the representative depositional facies of the studied intervals are sandy, fluvial, channel-fill facies encased in shaly floodplain deposits. Seismic responses were predicted from a synthetic seismogram using a model of dominant depositional facies. Peak-to-trough amplitude and instantaneous frequency seismic attributes are used in depositional facies interpretation. Three Intervals A, B and C are interpreted on the successive stratal surfaces. The shallowest interval, A, is the Quaternary transgressive succession. Each stratal surface showed flow pattern variation of fluvial channel facies. Two transgressive cycles were identified in interval A. Interval B also indicated fluvial facies. Depositional facies architectures are described by interpreting seismic attributes on the successive stratal surfaces.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2003.11a
• /
• pp.723-730
• /
• 2003

In general, well log and core data have been utilized for reservoir characterization. These well data can provide valuable information on reservoir properties with high vertical resolution at well locations. While the seismic surveys cover large areas of field but give only indirect features about reservoir properties. Therefore it is possible to estimate the reservoir properties guided by seismic data on entire area if a relationship of seismic data and well data can be defined. Seismic attributes calculated from seismic surveys contain the particular reservoir features, so that they should be extracted and used properly according to the purpose of study. The method to select the suitable seismic attributes among enormous ones is needed. The stepwise regression and fuzzy curve analysis based on fuzzy logics are used for selecting the best attributes. The relationship can be utilized to estimate reservoir properties derived from seismic attributes. This methodology is applied to a synthetic seismogram and a sonic log acquired from velocity model. Seismic attributes calculated from the seismic data are reflection strength, instantaneous phase, instantaneous frequency and pseudo sonic logging data as well as seismic trace. The fuzzy curve analysis is used for choosing the best seismic attributes compared to sonic log as well data, so that seismic trace, reflection strength, instantaneous frequency, and pseudo sonic logging data are selected. The relationship between the seismic attribute and well data is found out by the statistical regression method and estimates the reliable well data at a specific field location derived from only seismic attributes. For a future work in this study, the methodology should be checked an applicability of the real fields with more complex and various reservoir features.