• Economic and Environmental Geology
• /
• v.31 no.3
• /
• pp.235-247
• /
• 1998

A new technique of simultaneons inversion for 3-D seismic velocity structure by using direct, reflected, and refracted waves is applied to the center of the Korean Peninsula including Pyongnam Basin, Kyonggi Massif, Okchon Fold Zone, Taebaeksan Fold Zone, Ryongnam Massif and Kyongsang Basin. Pg, Sg, PmP, SmS, Pn, and Sn arrival times of 32 events with 404 seismic rays are inverted for locations and crustal structure. 5 ($1^{\circ}$ along the latitude)${\times}6$ ($0.5^{\circ}$ along the longitude) ${\times}8$ block (4 km each layer) model was inverted. 3-D seismic crustal velocity tomography including eight sections from the surface to the Moho, eight profiles along latitude and longitude and the Moho depth distribution was determined. The results are as follows: (1) the average velocity and thickness of sediment are 5.15 km/sec and 3-4 km, and the velocity of basement is 6.12 km/sec. (2) the velocities fluctuate strongly in the upper crust, and the velocity distribution of the lower crust under Conrad appears basically horizontal. (3) the average depth of Moho is 29.8 km and velocity is 7.97 km/sec. (4) from the sedimentary depth and velocity, basement thickness and velocity, form of the upper crust, the Moho depth and form of the remarkable crustal velocity differences among Pyongnam Basin, Kyonggi Massif, Okchon Zone, Ryongnam Massif and Kyongsang Basin can be found. (5) The different crustal features of ocean and continent crust are obvious. (6) Some deep index of the Chugaryong Rift Zone can be located from the cross section profiles. (7) We note that there are big anisotropy bodies near north of Seoul and Hongsung in the upper crust, implying that they may be related to the Chugaryong Rift Zone and deep fault systems.

• Geophysics and Geophysical Exploration
• /
• v.12 no.2
• /
• pp.173-182
• /
• 2009

In April 2008, KIGAM carried out an ocean-bottom seismometer (OBS) survey in the central Ulleung Basin where strong bottom simulating reflectors (BSRs) were revealed from previous surveys and some gas-hydrate samples were retrieved by direct sampling. The purpose of this survey is to estimate the velocity structure near the BSR in the gas hydrate prospect area using wide-angle seismic data recorded on the ocean-bottom seismometers. Along with the OBS survey, a 2-D seismic survey was performed whereby stratigraphic and preliminary velocity information was obtained. Two methods were applied to wide-angle data for estimating P wave velocity; one is velocity analysis in the $\tau$-p domain and the other is seismic traveltime inversion. A 1-D interval velocity profile was obtained by the first method, which was refined to layered velocity structure by the latter method. A layer stripping method was adopted for modeling and inversion. All velocity profiles at each OBS site clearly show velocity reversal at BSR depths due to the presence of gas hydrates. In addition, we could confirm high velocity in the column/chimney structure.

• Geophysics and Geophysical Exploration
• /
• v.19 no.4
• /
• pp.200-211
• /
• 2016

To support cross potential array and direct potential array, the array for leakage detection of all kinds of water facilities is proposed and it is named as the D-Lux array. The D-Lux array data are arranged to a coloured matrix and it is called the D-Lux view. Low potential difference of anomalous zone shown in D-Lux view implies the indication of leakage zone. Furthermore, for an intuitive interpretation of D-Lux array, equipotential distribution map is made by using D-Lux and direct potential array data. Equipotential distribution map makes us possible to predict import point, export point and the path of water leakage that we could have not anticipated in D-Lux view and the graphs. The water tank experiment and numerical analysis were carried out as preparatory experiment and the field explorations were conducted at a concrete weir and a fill dam. As a result, effective and specific detection of leakage path was possible for the concrete weir and the fill dam.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.119-124
• /
• 2008

Due to the recent development of instruments which enabled the measurements of subtle IP effect in the ground and analysis algorithms, complex resistivity (CR) method was expanding its application to various field. In this study, we applied the CR method to the test site where the ground reinforcement had been done by injecting the cement mortar for investigating the effect of ground reinforcement. For this site, resistivity monitoring and tomography survey was carried out while the ground reinforcement had been made by the grouting. From the result, the anomalous region that was shown on the result of resistivity 4D monitoring was coincident with those of phase section in the CR method, because the cement grouting material had the strong IP effects. It might be expected that the CR method should be very powerful surveying tool for the similar purpose.

• Geophysics and Geophysical Exploration
• /
• v.4 no.3
• /
• pp.70-79
• /
• 2001

The integral equation method is a powerful tool for numerical electromagnetic modeling. But the difficulty of this technique is the size of the linear equations, which demands excessive memory and calculation time to invert. This limitation of the integral equation method becomes critical in inverse problem. The conventional Born approximation, where the electric field in the anomalous body is approximated by the background field, is very rapid and easy to compute. However, the technique is inaccurate when the conductivity contrast between the body and the background medium is large. Quasi-linear, quasi-analytical and extended Born approximations are novel approaches to 3-D EM modeling based on the linearization of the integral equations for scattered EM field. These approximation methods are much less time consuming than full integral equation method and more accurate than conventional Born approximation. They we, however, still approximate methods for 3-D EM modeling. Iterative series methods such as modified Born, quasi-linear and quasi-analytical can be used to increase the accuracy of various approximation methods. Comparisons of numerical performance against a full integral equation and various approximation codes show that the iterative series methods are very accurate and almost always converge. Furthermore, they are very fast and easy to implement on a computer. In this study, extended Born series method is developed and it shows more accurate result than that of other series methods. Therefore, Iterative series methods, including extended Born series, open principally new possibilities for fast and accurate 3-D EM modeling and inversion.

• Geophysics and Geophysical Exploration
• /
• v.8 no.2
• /
• pp.145-155
• /
• 2005

A two-dimensional (2-D) interpretation of MT data has been performed for the purpose of fracture detection for geothermal development. Remote stations have been operated in Kyushu, Japan (480 km apart) as well as in Korea (60 km and 165 km apart in 2002 and 2003 data set, respectively). Apparent resistivity and phase curves calculated by remote processing with the Japan remote data showed enough quality for 2-D inversion for the whole frequency range. Remote reference processing with Korea remote reference data also showed quite good continuity in apparent resistivity and phase curves except some noisy frequency bands; around the power frequency, 60 Hz, and around the dead band $10^{-1}Hz\;Hz\;\~1\;Hz$, where the natural EM signal is known to be very weak. Even though the subsurface showed severe three-dimensional (3-D) characteristics in the survey area so that 2-D inversion by itself could not give enough information for deep geological structures, the 2-D inversion for the 5 survey lines showed several common features. The conductive semi-consolidate mudstone layer is dipping from north to south (about 500 m depth on the south and 200 m on the north most part of the survey area). The boundary between the low (L-2) and high (H-2) resistivity anomalies can be thought as a major fault with strike $N15^{\circ}E$, passing through the sites 206, 112 and 414. The shallow (< 1 km) conductive anomalies (L-4) seem to be fracture zones having strike E-W (at site 105) and $N60^{\circ}W$ (at site 434). And there exists a conductive layer in the western and west-southern part of the survey area in the depth below $2\~3\;km$, for which further investigation is to be needed.

• Geophysics and Geophysical Exploration
• /
• v.12 no.4
• /
• pp.361-366
• /
• 2009

Layered-earth Green's functions in electormagnetic (EM) surveys play a key role in modeling the response of exploration targets. They are computed through the Hankel transforms of analytic kernels. Computational precision depends upon the choice of algebraically equivalent forms by which these kemels are expressed. Since three-dimensional (3D) modeling can require a huge number of Green's function evaluations, total computational time can be influenced by computational time for the Hankel transform evaluations. Linear digital filters have proven to be a fast and accurate method of computing these Hankel transforms. In EM modeling for 3D inversion, electric fields are generally evaluated by the secondary field formulation to avoid the singularity problem. In this study, three components of electric fields for five different sources on the surface of homogeneous half-space were derived as primary field solutions. Moreover, reflection coefficients in TE and TM modes were produced to calculate EM responses accurately for a two-layered model having a sea layer. Accurate primary fields should substantially improve accuracy and decrease computation times for Green's function-based problems like MT problems and marine EM surveys.

• 한국해양학회지
• /
• v.29 no.3
• /
• pp.228-238
• /
• 1994

The change of sea-level is a good indicator of the change of climate during the Quaternary period. The sea-levels in the world have been changing very irregularly during that time. The pattern of the Quaternary sea-level change was assumed to be a stochastic fractal in this study. We measured fractal dimensions of the Holocene sea-levels of the Hudson river estuary and the Delaware coast. A box counting method gave almost the same values. i.e., D=1.358 for the Hudson sea-level changes and D+1.346 for the Delaware sea-level changes. the ability of the inverse method of fractal interposea-levels. IFIF reproduction the realistic sea-levels for the both of them. The delaware sea-level data made less statistical errors for the interpolation of IFIF than the Hudson and the Delaware sea-levels. IFIF reproduction the realistic sea-levels for the both of them. The Delaware sea-level data made less statistical errors for the interpolation of IFIF than the Hudson sea-level data. This suggests that the Delaware sea-level data are more reliable than the Hudson sea-level data was calculated from the fractal dimension of the Delaware sea-level data. Fractal interpolation functions (FIF) was used to reconstruct the peleosea-levels of the Korean coasts and the Atlantic Ocean coasts of the United States. The Korean Peleosea-level change generacted by FIF is different from the peleosea-level change of the eastern U.S.. The Korean peleosea-levels are much higher than the eastern U.S. Paleosea-levels, comparing to each other from the present to 8,000 BP.

• Economic and Environmental Geology
• /
• v.27 no.5
• /
• pp.479-489
• /
• 1994

In this paper, the iterative least-squares inversion method is used to determine shapes and density contrasts of 2-D structures from the gravity data. The 2-D structures are represented by their cross-sections of N-sided polygons with density contrasts which are constant or varying with depth. Gravity data are calculated by theoretical formulas for the above structure models. The data are considered as observed ones and used for inversions. The inversions are performed by the following processes: I) polygon's vertices and density contrast are initially assumed, 2) gravity are calculated for the assumed model and error between the true (observed) and calculated gravity are determined, 3) new vertices and density contrast are determined from the error by using the damped least-squares inversion method, and 4) final model is determined when the error is very small. Results of this study show that the shape and density contrast of each model are accurately determined when the density contrast is constant or vertical density gradient is known. In case where the density gradient is unknown, the inversion gives incorrect results. But the shape and density gradient of the model are determined when the surface density contrast is known.

• Geophysics and Geophysical Exploration
• /
• v.13 no.4
• /
• pp.397-406
• /
• 2010

Resistivity method is a practical and effective geophysical technique to detect leakage zones in embankment dams. Generally, resistivity survey conducted along the crest assumes that the embankment dam has a 2D structure. However, the 3D topography of the embankment distorts significantly resistivity data measured on anywhere of the dam. This study evaluates the influence from 3D effects created by specific dam geometry and effects of water level fluctuations through the 3D finite element modeling technique. Also, a comparison between different locations of survey line are carried out, and topographic correction technique is developed for the resistivity data obtained along the embankment dam. Furthermore, using synthetic resistivity data for an embankment dam model with leakage zone, detectability of leakage zones is estimated through 2.5D inversion.