• Proceedings of the KSEEG Conference
• /
• 2002.04a
• /
• pp.163-166
• /
• 2002

Initial velocity model close to real velocity structure of the subsurface governs the quality of image of prestack depth migration. Geophysicists employ velocity estimation tools such as velocity analysis (curvature method, coherency inversion), tomography and waveform inversion. We present a reflection tomography that parameterizes the subsurface into the movable blocks. By interpreting the depth-migrated section or stacked section, we can design an initial constant velocity model having only impedance boundaries. We use shooting-raytracing method that allows us to calculate the Jacobian-matrix efficiently.

• Proceedings of the KSEEG Conference
• /
• 2002.04a
• /
• pp.159-162
• /
• 2002

We tried to obtain an initial velocity model for prestack depth migration via waveform inversion. For application of any field data we chose a smooth background layered velocity model (v=v0 + k x z) as an initial velocity model. Newton type waveform inversion needs to invert huge Hessian matrix. In order to compute full Hessian matrix arising from full aperture data and full illumination zone, we meet insurmountable difficulties of paying astronomical computing cost. For the layered media, approximate Hessian emerging from single shot aperture data can be used repeatedly for split spread source configuration. In our work of using this Hessian characteristic of layered media we attempted to obtain the approximate velocity model as close as possible to the true velocity model in first iteration.

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.18-24
• /
• 2011

A multi-channel seismic reflection (MCS) survey was conducted in 2009 to explore the deep crustal structure of the Pacific Plate south of Hokkaido. The survey line happened to traverse a 250-km-wide Warm Core Ring (WCR), a current eddy that had been generated by the Kuroshio Extension. We attempted to use these MCS data to delineate the WCR fine structure. The survey line consists of two profiles: one with a shot interval of 200m and the other with a shot interval of 50 m. Records from the denser shot point line show much higher background noise than the records from the sparser shot point line. We identified the origin of this noise as acoustic reverberations between the sea surface, seafloor and subsurface discontinuities, from previous shots. Results showed that a prestack migration technique could enhance the signal buried in this background noise efficiently, if the sound speed information acquired from concurrent temperature measurements is available. The WCR is acoustically an assemblage of concave reflectors dipping inward, with steeper slopes (${\sim}2^{\circ}$) on th ocean side and gentler slopes (${\sim}1^{\circ}$) on the coastal side. Within the WCR, we recognised a 30-km-wide lens-shaped structure with reflectors on the perimeter.

• Geophysics and Geophysical Exploration
• /
• v.16 no.2
• /
• pp.71-78
• /
• 2013

The one-way wave equation migration is much more computationally efficient comparing with reverse time migration and it can provide better image than the migration algorithm based on the ray theory. We have developed the prestack depth migration module adopting (GS) propagator designed for vertical transverse isotropic media. Since GS propagator considers the higher-order term by expanding the Taylor series of the vertical slowness in the thin slab of the phase-screen propagator, the GS migration can offer more correct image for the complex subsurface with large lateral velocity variation or steep dip. To verify the validity of the developed GS migration module, we analyzed the accuracy with the order of the GS propagator for VTI media (GSVTI propagator) and confirmed that the accuracy of the wavefield propagation with the wide angles increases as the order of the GS propagator increases. Using the synthetic seismic data, we compared the migration results obtained from the isotropic GS migration module with the anisotropic GS migration module. The results show that the anisotropic GS migration provides better images and the improvement is more evident on steeply dipping structures and in a strongly anisotropic medium.

• Geophysics and Geophysical Exploration
• /
• v.14 no.3
• /
• pp.227-233
• /
• 2011

Attention has been focused on geothermal energy as an alternative energy because it is continuously operable without external supply. Most of geothermal anomalies in Korea are related to deep circulation of groundwater through a fracture system in granite area. Therefore it is very important to understand the distribution of the fracture system which is the main channel of ground water. In this research, we constructed the velocity models with a fracture system and the layered sediments, respectively, and generated synthetic data sets with them to verify the presented vertical seismic profiling (VSP) preprocessing scheme. We compared the results from conventional VSP preprocessing flow to those from VSP preprocessing flow considering fracture system. We noticed that the preprocessing flow considering fracture system retains more sufficient signal including down-going wave than conventional preprocessing. In addition, we applied 3D VSP prestack phase screen migration to the preprocessed reversed VSP (RVSP) data from Seokmo Island so that we were able to image fracture structure of the geothermal site in Seokmo Island.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2006.06a
• /
• pp.131-135
• /
• 2006

A seismic waveform inversion for prestack seismic data based on the Gauss-Newton method is presented. The Gauss-Newton method for seismic waveform inversion was proposed in the 80s but has rarely been studied. Extensive computational and memory requirements have been principal difficulties. To overcome this, we used different sizes of grids in the inversion stage from those of grids in the wave propagation simulation, temporal windowing of the simulation and approximation of virtual sources for calculating partial derivatives, and implemented this algorithm on parallel supercomputers. We show that the Gauss-Newton method has high resolving power and convergence rate, and demonstrate potential applications to real seismic data.

• Geophysics and Geophysical Exploration
• /
• v.5 no.1
• /
• pp.18-22
• /
• 2002

As a single arrival traveltime, maximum energy arrival traveltime has been known as the most proper operator for Kirchhoff migration. In case of the model having the simple structure, both the first arrival traveltime and the maximum energy arrival traveltime can be used as the correct operators for Kirchhoff migration. However for some model having the complex and high velocity contrast structure, the migration using the first arrival traveltime can't give the correct depth section. That is, traveltime to be required in Kirchhoff migration is the maximum energy traveltime, but, needs considerably more calculation time than that of first arrival. In this paper, we propose the method for calculating the traveltime approximated to the maximum energy arrival using one-way wave equation. After defining the WAS(Wrap Around Suppression) factor to be used for calculating the first arrival traveltime using one-way wave equation as the function of lateral grid interval and depth and considering the delay time of source wavelet. we calculate the traveltime approximated to the maximum energy arrival. to verify the validity of this traveltime, we applied this to the migraion for simple structure and complex structure and compared the depth section with that obtained by using the first arrival traveltime.

• Geophysics and Geophysical Exploration
• /
• v.6 no.2
• /
• pp.57-63
• /
• 2003

For the dectection of small cavity in the hard rock, we investigated the feasibility of crosswell travel-time tomography and Kirchhoff migration technique. In travel-time tomography, first arrival anomaly caused by small cavity was investigated by numerical modeling based on the knowledge of actual field information. First arrival delay was very small (<0.125 msec) and detectable receiver offset range was limited to 4m with respect to $1\%$ normalized first arrival anomaly. As a consequence, it was turned out that carefully designed survey array with both sufficient narrow spatial spacing and temporal (<0.03125 msec) sampling were required for small cavity detection. Also, crosswell Kirchhoff migration technique was investigated with both numerical and real data. Stack section obtained by numerical data shows the good cavity image. In crosswell seismic data, various unwanted seismic events such as direct wave and various mode converted waves were alto recorded. To remove these noises und to enhance the diffraction signal, combination of median and bandpass filtering was applied and prestack and stacked migration images were created. From this, we viewed the crosswell migration technique as one of the adoptable method for small cavity detection.

• Geophysics and Geophysical Exploration
• /
• v.1 no.1
• /
• pp.8-18
• /
• 1998

In exploration seismology, the Kirchhoff hyperbola has been successfully used to migrate reflection seismo-grams. The mathematical basis of Kirchhoff hyperbola has not been clearly defined and understood for the application of prestack or poststack migration. The travel time from the scatterer in the subsurface to the receivers (exploding reflector model) on the surface can be a kinematic approximation of Green's function when the source is excited at position of the scatterer. If we add the travel time from the source to the scatterer in the subsurface to the travel time of exploding reflector model, we can view this travel time as a kinematic approximation of the partial derivative wavefield with respect to the velocity or the density in the subsurface. The summation of reflection seismogram along the Kirchhoff hyperbola can be evaluated as an inner product between the partial derivative wavefield and the field reflection seismogram. In addition to this kinematic interpretation of Kirchhoff hyperbola, when we extend this concept to shallow refraction seismic data, the stacking of refraction data along the straight line can be interpreted as a measurement of an inner product between the first arrival waveform of the partial derivative wavefield and the field refraction data. We evaluated the Kirchhoff hyperbola and the straight line for stacking the refraction data in terms of the first arrival waveform of the partial derivative wavefield with respect to the velocity or the density in the subsurface. This evaluation provides a firm and solid basis for the conventional Kirchhoff migration and the straight line stacking of the refraction data.

• Geophysics and Geophysical Exploration
• /
• v.23 no.1
• /
• pp.13-21
• /
• 2020

The acquisition and processing of long-offset data are essential for imaging deep geological structures in marine seismic surveys. It is challenging to derive an accurate subsurface image by employing conventional data processing to long-offset data owing to the normal moveout (NMO) stretch and non-hyperbolic moveout phenomena induced by seismic anisotropy. In 2017, the Korea Institute of Geoscience and Mineral Resources conducted a simultaneous two-dimensional multichannel streamer and ocean-bottom seismic survey using a 5.7-km streamer and an ocean-bottom seismometer to identify the deep geological structure of the Ulleung Basin. Herein, the actual geological subsurface structure was obtained via the sequential iterative updating of the velocity and anisotropic parameters of the long-offset data obtained using a multichannel streamer, and anisotropic prestack Kirchhoff migration was performed using the updated velocity and anisotropic parameters as input parameters. As a result, the reflection energy in the long-offset traces, which showed non-hyperbolic moveout owing to seismic anisotropy, was well aligned horizontally and NMO stretches were also reduced. Thus, a more precise and accurate migrated image was obtained, minimizing the distortion of reflectors and mispositioned reflection energy.