• 한국지구물리탐사학회:학술대회논문집
• /
• 2005.09a
• /
• pp.21-36
• /
• 2005

Although the main geology of Korea consists of granite and gneiss, it Is not uncommon to encounter anisotropy Phenomena in crosshole radar tomography even when the basement is crystalline rock. To solve the anisotropy Problem, we have developed and continuously upgraded an anisotropic inversion algorithm assuming a heterogeneous elliptic anisotropy to reconstruct three kinds of tomograms: tomograms of maximum and minimum velocities, and of the direction of the symmetry axis. In this paper, we discuss the developed algorithm and introduce some case histories on the application of anisotropic radar tomography in Korea. The first two case histories were conducted for the construction of infrastructure, and their main objective was to locate cavities in limestone. The last two were performed In a granite and gneiss area. The anisotropy in the granite area was caused by fine fissures aligned in the same direction, while that in the gneiss and limestone area by the alignment of the constituent minerals. Through these case histories we showed that the anisotropic characteristic itself gives us additional important information for understanding the internal status of basement rock. In particular, the anisotropy ratio defined by the normalized difference between maximum and minimum velocities as well as the direction of maximum velocity are helpful to interpret the borehole radar tomogram.

• Geophysics and Geophysical Exploration
• /
• v.11 no.3
• /
• pp.250-262
• /
• 2008

In order to assess the feasibility of the dipole-dipole electric method to the investigation of metallic ore deposit, both field data simulation and inversion are carried out for several simplified ore deposit models. Our interest is in a vein-type model, because most of the ore deposits (more than 70%) exist in a vein type in Korea. Based on the fact that the width of the vein-type ore deposits ranges from tens of centimeters to 2m, we change the width and the material property of the vein, and we use 40m-electrode spacing for our test. For the vein-type model with too small width, the low resistivity zone is not detected, even though the resistivity of the vein amounts to 1/300 of that of the surrounding rock. Considering a wide electrode interval and cell size used in the inversion, it is natural that the size of the low resistivity zone is overestimated. We also perform field data simulation and inversion for a vein-type model with surrounding hydrothermal alteration zones, which is a typical structure in an epithermal ore deposits. In the model, the material properties are assumed on the basis of resistivity values directly observed in a mine originated from an epithermal ore deposits. From this simulation, we can also note that the high resistivity value of the vein does not affect the results when the width of the vein is narrow. This indicates that our main target should be surrounding hydrothermal alteration zones rather than veins in field survey. From these results, we can summarize that when the vein is placed at the deep part and the difference of resistivity values between the vein and the surrounding rock is not large enough, we cannot detect low resistivity zone and interpret the subsurface structures incorrectly using the electric method performed at the surface. Although this work is a little simple, it can be used as references for field survey design and field data Interpretation. If we perform field data simulation and inversion for a number of models and provide some references, they will be helpful in real field survey and interpretation.

• Geophysics and Geophysical Exploration
• /
• v.11 no.3
• /
• pp.242-249
• /
• 2008

Recently, AVO analysis has been widely used in oil exploration with seismic subsurface section as a direct indicator of the existence of the gas. In the case of the deep reservoirs like the gas reservoirs in the East-sea, it is often difficult to observe AVO responses in CMP gathers even though the bright spots are shown in the stacked section. Because the reservoir becomes more consolidated as its depth deepens, P-wave velocity does not decrease significantly when the pore fluid is replaced by the gas. Thus the difference in Poisson's ratio, which is a key factor for AVO response, between the reservoir and the layer above it does not increase significantly. In this study, we analyzed the effects of Poisson's ratio difference on AVO response with a variety of Poisson's ratios for the upper and lower layers. The results show that, as the difference in Poisson's ratio between the upper and lower layers decreases, the change in the reflection amplitude with incidence angle decreases and AVO responses become insignificant. To consider the limitation of AVO responses shown in the gas reservoir in East-sea, the velocity model was made by simulation Gorae V structure with seismic data and well logs. The results of comparing AVO responses observed from the synthetic data with theoretical AVO responses calculated by using material properties show that the amount of the change in reflection amplitude with increasing incident angle is very small when the difference in Poisson's ratio between the upper and lower layers is small. In addition, the characteristics of AVO responses were concealed by noise or amplitude distortion arisen during preprocessing. To overcome such limitations of AVO analysis of the data from deep reservoirs, we need to acquire precisely reflection amplltudes In data acquisition stage and use processing tools which preserve reflection amplitude in data processing stage.

• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.89-97
• /
• 2007

A three-dimensional (3D) magnetotelluric (MT) survey has been carried out to delineate subsurface structures and possible fractures, for development of low-temperature geothermal resources in Pohang, Korea. Quite good quality MT data could be obtained throughout the survey region by locating the remote reference in Kyushu, Japan, which is ${\sim}480\;km$ from the centre of the field site. 3D modelling and inversion are performed taking into account the sea effect in MT measurements near the seashore. The nearby sea in the Pohang area affects MT data at frequencies below $1\;Hz{\sim}0.2\;Hz$, depending on the distance from the seashore. The most severe sea effects were observed in the south-east parts of the survey area, closer to Youngil Bay. 3D inversion with and without the seawater constraint showed very similar results at shallow depths, roughly down to 2 km. At greater depths, however, a strong sea effect seems to form a fictitious conductive structure in ordinary 3D inversion, especially in the south-eastern part of the survey region. Comparison between drilling results and the resistivity profiles from inversions showed that five layered structures can be distinguished the subsurface beneath the target area. They are: (a) semi-consolidated mudstones with resistivity less than $10\;{\Omega}m$, which are ${\sim}300\;m$ thick in the northern part and ${\sim}600\;m$ thick in the southern part of the survey area; (b) occasional occurrence of trachybasalt and lapilli tuff within the mudstone layer has resistivity of a few tens of${\Omega}m$, (c) intrusive rhyolite ${\sim}400\;m$ thick has resistivity of several hundreds of ${\Omega}m$, (d) alternating sandstone and mudstone down to 1.5 km depth shows resistivity of ${\sim}100\;{\Omega}m$, (e) a conductive structure was found at a depth of ${\sim}3\;km$, but more geological and geophysical study should be carried out to identify this structure.

• Geophysics and Geophysical Exploration
• /
• v.6 no.3
• /
• pp.126-133
• /
• 2003

High-resolution shallow marine seismic surveys have been carried out for the resources exploration, engineering applications and Quaternary mapping. To improve the resolution of subsurface structure image, multichannel digital technique has been applied. The quality of the image depends on the vertical and horizontal resolution and signal to noise (S/N) ratio which are associated with the data acquisition parameters such as sample interval, common midpoint (CMP) interval and CMP fold. To understand the effect of the acquisition parameters, a test survey was carried out off Yeosu and the acquired data were analyzed. A 30 $in^3$ small air gun was used as a seismic source and 8 channel streamer cable with a 5 m group interval was used as a receiver. The data were digitally recorded with a shot interval of 2 s and sample interval of 0.1 ms. The acquired data were resampled with various sample intervals, CMP intervals and CMP folds. The resampled data were processed, plotted as seismic sections and compared each other. The analysis results show that thin bed structure with ${\~}1m$ thickness and ${\~}6^{\circ}$ slope can be imaged with good resolution and continuity and low noise using the acquisition parameters with a sample interval shorter than 0.2 ms, CMP interval shorter than 2.5 m and CMP fold more than 4. Because seismic resolution is associated with the acquisition parameters, the quality of the subsurface structure can be imaged successfully using suitable and optimum acquisition parameters.

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.127-136
• /
• 2008

To identify the potential area of gas hydrate in the Ulleung Basin, 2-D and 3-D seismic surveys using R/V Tamhae II were conducted in 2005 and 2006. Seismic survey equipment consisted of navigation system, recording system, streamer cable and air-gun source. For reliable velocity analysis in a deep sea area where water depths are mostly greater than 1,000 m and the target depth is up to about 500 msec interval below the seafloor, 3-km-long streamer and 1,035 $in^3$ tuned air-gun array were used. During the survey, a suite of quality control operations including source signature analysis, 2-D brute stack, RMS noise analysis and FK analysis were performed. The source signature was calculated to verify its conformity to quality specification and the gun dropout test was carried out to examine signature changes due to a single air gun's failure. From the online quality analysis, we could conclude that the overall data quality was very good even though some seismic data were affected by swell noise, parity error, spike noise and current rip noise. Especially, by checking the result of data quality enhancement using FK filtering and missing trace restoration technique for the 3-D seismic data inevitably contaminated with current rip noises, the acquired data were accepted and the field survey could be conducted continuously. Even in survey areas where the acquired data would be unsuitable for quality specification, the marine seismic survey efficiency could be improved by showing the possibility of noise suppression through onboard data processing.

• Journal of the Korean Geophysical Society
• /
• v.5 no.4
• /
• pp.283-292
• /
• 2002

The Korean peninsula is located at the edge of the East Asian active margin. The seismic activity in the Korean Peninsula is relatively low compared with the neighboring countries China and Japan. According to the available Seismic information, the Korean Peninsula is not totally safe from the Earthquake disaster. Moreover, the area is surrounded by varies tectonic forces which is resulted from the relative movements of the surrounding tectonic plates "Pacific, Philippine Sea, Eurasian and South China". Nowadays South Korea has 65 GPS stations belong to 5 governmental organizations "each organization figure out its own GPS stations for different requirements" In order to minimize the seismic hazard in the Korean Peninsula a program for monitoring the recent crustal movement has been designed considering the uses of the available GPS station "some selected stations from the previously mentioned stations" and the tectonic settings in and around the Korean Peninsula. This program is composed of two main parts, the first part to monitor the crustal deformation around the Korean Peninsula with the collaboration of the surrounding countries "China and Japan" this part is composed of two phases "East Sea Phase and Yellow Sea Phase". These phases will be helpful in determining the deformation parameters in the East Sea and the Yellow Sea respectively While the Second part of this program, is designed to determine the deformation parameters id and around the main faults in the Korean Peninsula and the relative movement between the Korean Peninsula and the Cheju Island. Through out this study the needs of crustal movement center rose up to collect the data from the previously mentioned stations and Organizations in order to use such reliable data in different geodynamical application.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2005.05a
• /
• pp.103-107
• /
• 2005

In order to investigate the velocity structure of the southern part of the Korean peninsula, exploded seismic signals were recorded for 120 s along a 294-km WNW-ESE line and 150 s along a 335-km NNW-SSE line in 2002 and 2004, respectively. Velocity tomograms were derived from inverting first arrival times. One-dimensional velocity models derived by joint analyses of teleseismic receiver functions and surface wave dispersion at several stations near the profiles were uesd to build initial models. The raypaths indicate several midcrust interfaces including ones at approximate depths of 2.0 and 14.9 km with refraction velocities of approximately 6.0 and 7.1 km/s, respectively. The deepest significant interface varies in depth from 30.8 km to 36.1 km. The critically refracting velocity varies from 7.8 to 8.1 km/s along this interface which may correspond to the Moho discontinuity. The velocity tomograms show (1) existence of a low-velocity zone centered at 6-7 km depth under the Okchon fold belt, (2) extension of the Yeongdon fault down to greater than 10 km, and (3) existence of high-velocity materials under the Gyeongsan basin whose thickness is less than 4.2 km.

• Journal of the Korean Geophysical Society
• /
• v.3 no.1
• /
• pp.37-48
• /
• 2000

In order to investigate in-line ground motions caused by earthquakes, we examine the multicomponent complex trace analysis method (MCTAM) for the synthetic data and apply it to real earthquake data. An experimental result for synthetic data gives correct information on the arrival times, duration of individual phases, and approaching angles for body waves. Rayleigh waves are also easily identified with the MCTAM. A deep earthquake with magnitude of 7.3 was chosen to test various polarization attributes of ground motions. For P waves, instantaneous phase difference between the vertical and the in-line horizontal components ${\phi}(t)$, instantaneous reciprocal ellipticity ${\rho}(t)$, and approaching angle ${\tau}(t)$ are computed to be ${\pm}180^{\circ},\;0{\sim}0.25,\;and\;-30^{\circ}{\sim}-45^{\circ}$, respectively. For S waves, ${\phi}(t)$ tends to vary while ${\rho}(t)$ have values of $0{\sim}0.3\;and\;{\tau}(t)$ remains near vertical, respectively. A relatively low frequency signal registered just prior to the S wave event is interpreted as a P-wave phase based on its polarization characteristics. Velocities of P and S waves are computed to be 8.633 km/s and 4.762 km/s, and their raypath parameters 0.074 s/km and 0.197 s/km. Dynamic Poisson's ratio is obtained as 0.281 from the velocities of P and S waves.

• Journal of the Korean Geophysical Society
• /
• v.3 no.4
• /
• pp.215-226
• /
• 2000

In order to map the acoustic basement and to locate fracture zones in the Galgok fault, seismic refraction data were acquired near the Chonbuk ranch in Gyeongju. Along three profiles of 72m(Line 1), 72m(Line 2), and 36m(Line 3) long, seismic signals were generated by a 5kg hammer. The refraction data were collected by employing twelve 8 Hz geophones at an interval of 3m and recording time of 192ms at a sampling rate of 0.2ms. The data are interpreted using GRM method. The top layer (Layer 1) is characterized as the velocity of approximately250 m/s and thickness of approximately 2.1m. This layer is regarded as a soil layer. Underneath Layer 1 lies unconsolidated layer (Layer 2) whose refraction velocity is determined to be $1,030{\sim}1,400m/s$. Layer 2 is approximately 4.6m thick and is regarded as a Quaternary gravel layer. The third layer (Layer 3) has the mean refraction velocity of $2,100{\sim}2,200m/s$ and is interpreted to be the acoustic basement. In some parts of Lines 1 and 3, the difference in depth to the top of Layer 2 is greater than 20 cm indicating the possibility of existence of Quaternary faults. Along Line 3 and the eastern part of Line 1, refracted energy from the acoustic basement was not recorded. This may highly indicate that a relatively large scale fault exists under the western part of Line 1.