• 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.14 no.3
• /
• pp.234-241
• /
• 2011

Recently, exploration and development of oil sands are thriving due to high oil price. Because oil sands reservoir usually exists as a thin layer, multicomponent VSP, which has the advantage of the high-resolution around the borehole, is more effective than surface seismic survey in exploring oil sand reservoir. In addition, prestack phase-screen migration is effective for multicomponent seismic data because it is based on an one-way wave equation. In this study, we examined the applicability of the prestack phase-screen migration for multicomponent RVSP data to image the thin oil sand reservoir. As a preprocessing tool, we presented a method for separating P-wave and PS-wave from multicomponent RVSP data by using incidence angle and rotation matrix. To verify it, we have applied the developed wavefield separation method to synthetic data obtained from the velocity model including a horizontal layer and dipping layers. Also, we compared the migrated image by using P-wave with that by using PS-wave. As a result, the PS-wave migrated image has higher resolution and wide coverage than P-wave migrated image. Finally, we have applied the prestack phase-screen migration to the synthetic data from the velocity model simulating oil sand reservoir in Canada. The results show that the PS-wave migrated image describe the top and bottom boundaries of the thin oil sand reservoir more clearly than the P-wave migrated image.

• 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.

• Journal of the Korean Geophysical Society
• /
• v.3 no.2
• /
• pp.77-90
• /
• 2000

Five electrical resistivity dipole-dipole and two seismic reflection surveys were performed in the southeastern margin of the Yongdong basin to delineate the shallow basin architecture. To investigate the intra-basin structure, twenty four resistivity sounding points and three dipole-dipole lines were selected especially in the vicinity of volcanic masses. The basin-fault boundaries are identified in electrical dipole-dipole resistivity section as high resistivity-contrast of approximately $1,500\;{\Omega}{\cdot}m$, characterized as a band of high standard-deviation. They are also effectively clarified in the seismic reflection data: amplitude and continuity contrasts in the common shot gather, first-arrival profiles, complex attribute plots. The intra-basin resistivity structures are constructed by interpolating vertical electrical sounding data and dipole-dipole profiles. The high-resistivity anomalies most likely originate from the northsouth-trending and northeast-dipping volcanic masses, which are to be further quantitatively investigated with geomagnetic and magnetotelluric surveys.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.393-398
• /
• 2007

Seismic traveltime tomography technique was conducted at a site composed of black shale. It is well known that black shale has strong anisotropic property. Therefore, the anisotropic property of black shale has to be considered to obtain the appropriate subsurface velocity model by an inversion process. To estimate the anisotropic constant of the velocity of the black shale in the survey area, the relation between the velocity, which is calculated by the straight ray path and the first arrival time, and the angle of the ray propagation was examined. The elliptically shaped relation was found and it reveals that the black shale contains the anisotropic property of velocity. It was also noticed that the horizontal velocity is faster than the vertical velocity. When the estimated anisotropic constant was applied in the process of the velocity inversion for three sets of field data, we could obtain the appropriate velocity structures of the site that is consistent with the result of the geological survey.

• The Journal of Engineering Geology
• /
• v.16 no.4 s.50
• /
• pp.327-335
• /
• 2006

Three processing strategies of seismic refraction data are tested in terms of velocity and depth profiles or structures for mapping of geological discontinuities: GRM(generalized reciprocal method), GLI(generalized linear inversion), Tomography. The test data used in this study are the shot gathers reconstructed by numerical modeling for the structures of 3 planar layers(horizontal, inclined), the buried vertical fracture zones, and vertical fault zones. Tomography is shown to be very efficient for mapping of more complicated tone such as vertical fault and buried fracture zones, whereas GRM and GLI can be useful for horizontal and/or inclined layers, probably on their bases of analysis of first arrivals in travel time curves.

• Journal of the Korean Geophysical Society
• /
• v.2 no.1
• /
• pp.45-56
• /
• 1999

Seismic refracrion and reflection surveys were conducted along an E-W trending track of 482 m long in Ilwall-dong, Pohang. End-on spread was employed as source-receiver configuration with 2 m for both geophone interval and offset. Seismic data were acquired using 24 channels at every shot fired every 2 m along the track. Refraction data were interpreted using equations for multi-horizontal layers. Reflection data were processed in the sequence of trace edit, gain control, CMP sorting, NMO correction, mute, common offset gathering, and filtering to produce a single fold seismic section. There are two layers in shallow subsurface of the study area. Upper layer has the P-wave velocities ranging from 267 to 566 m/s and is interpreted as a layer of unconsolidated sediments. Lower layer has P-wave velocities of 1096-3108 m/s and is interpreted as weathered rock to hard rock. Most of the lower layer classified as soft rock. Upper layer has lateral variations in both P-wave velocity and thickness. The upper layer in the eastern part of the seismic line is 3-5 m thick and has P-wave velocity of 400 m/s in average. The upper layer in the western part is 8-10 m thick and has P-wave velocity of 340 m/s in average. The eastern part is interpreted as unconsolidated beach sand, while the western part is interpreted as infilled soil to develop a construction site. Three fault systems of high angle are imaged in seismic reflection section. It is interpreted that the area between these fault systems are relatively safe. Large buildings should be located in the safe ground condition of no fault and footings should be designed to be in the basement rock of 3-10 m deep below the surface.

• Geophysics and Geophysical Exploration
• /
• v.11 no.3
• /
• pp.197-203
• /
• 2008

A shallow high-resolution seismic reflection survey was carried out at the Mineopo tidal flat on the western coast of Korea. The purpose of the survey was to investigate shallow sedimentary structure of the tidal flat associated with the recent sea level change. A total of 795 shots were generated at 1 m interval from a 5-kg hammer source and recorded on 48 channels of 100 Hz geophones along two mutually perpendicular profiles. The water-saturated ground condition resulted in suppressed ground rolls by significantly decreasing rigidity. In addition, seismic velocities over 1500 m/s provided easy segregation of reflected arrivals from lower velocity noise. As a consequence, seismic sections from the study area show significantly higher resolution and signal to noise ratio than conventional land seismic sections. The tidal flat consists of 5 sedimentary sequences above acoustic basement. The seismic sections reveal the continuous structure of the tidal flat formed in association with sea level rise during the Holocene.

• Geophysics and Geophysical Exploration
• /
• v.14 no.3
• /
• pp.220-226
• /
• 2011

In the elastic wave equations, both horizontal and vertical displacements are defined. Since we can measure both the horizontal and vertical displacements in field acquisition, these displacements compose a displacement vector. In this study, we propose a frequency-domain elastic waveform inversion technique taking advantage of the magnitudes of displacement vectors to define objective function. When we apply this displacement-vector objective function to the frequency-domain waveform inversion, the inversion process naturally incorporates the back-propagation algorithm. Through the inversion examples with the Marmousi model and the SEG/EAGE salt model, we could note that the RMS error of the solution obtained by our algorithm decreased more stably than that of the conventional method. Particularly, the density of the Marmousi model and the low-velocity sub-salt zone of the SEG/EAGE salt model were successfully recovered. Since the gradient direction obtained from the proposed objective function is numerically unstable, we need additional study to stabilize the gradient direction. In order to perform the waveform inversion using the displacementvector objective function, it is necessary to acquire multi-component data. Hence, more rigorous study should be continued for the multi-component land acquisition or OBC (Ocean Bottom Cable) multi-component survey.

• Journal of the Korean Geophysical Society
• /
• v.2 no.1
• /
• pp.75-88
• /
• 1999

A recent major subject of geophysical exploration is research into 3-D subsurface imaging with a composite information from the various geophysical data. In an attempt to interpret Schlumberger sounding data for the study area in 2-D and 3-D view, resistivity imaging was firstly performed and then pseudo-3-D resistivity volume was reconstructed by interpolating several 1-D resistivity plots. Electrical resistivity discontinuities such as fracture zone were successfully clarified in pseudo-3-D resistivity volume. The low resistivity zone mainly associated with fracture zone appears to develop down to granitic basement in the central part of the study area. Seismic velocity near the lineament is estimated to be approximately as small as 3,000 m/s, and weathering-layer for the southeastern part is interpreted to be deeper than for the northwestern part. Geophysical attributes such as electrical resistivity, seismic velocity, radioactivity for the Chojeong Area were analysed by utilizing a GIS software Arc/Info. The major fault boundaries and fracture zones were resolved through image enhancement of composite section (electrical resistivity and seismic refraction data) and were interpreted to develop in the southeastern part of the area, as characterized by low electrical resistivity and low seismic velocity. However, radioactivity attribute was found to be less sensitive to geological discontinuities, compared to resistivity and seismic velocity attributes.