• International Union of Geodesy and Geophysics Korean Journal of Geophysical Research
• /
• v.25 no.1
• /
• pp.23-34
• /
• 1997

Seismic tomography using the graph theoretical method of ray tracing is performed in two synthetic data sets with laterally varying velocity structures. The straight-ray tomography shows so poor results in imaging the laterally varying velocity structure that the ray-traced tomographic techniques should be used. Conventional ray tracing methods have serious drawbacks, i.e. problems of convergence and local minima, when they are applied to seismic tomography. The graph theretical method finds good approximated raypaths in rapidly varying media even in shadow zones, where shooting methods meet with convergence problems. The graph theoretical method ensures the globally minimal traveltime raypath while bending methods often cause local minima problems. Especially, the graph theoretical method is efficient in case that many sources and receivers exist, since it can find the traveltimes and corresponding raypaths to all receivers from a specific source at one time. Moreover, the algorithm of graph theoretical method is easily applicable to the ray tracing in anisotropic media, and even to the three dimensional case. Among the row-active inversion techniques, the conjugate gradient (CG) method is used because of fast convergence and high efficiency. The iterative sequence of the ray tracing by the graph theoretical method and the inversion by the CG method is an efficient and robust algorithm for seismic tomography in laterally varying velocity structures.

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

The estimation of seawater intrusion into deep aquifers has been becoming an important subject in terms of site characterization for geological disposal of radioactive waste. Conventional direct-current resistivity methods have been used for ground water explorations and recently have been applied to environmental problems. However, electromagnetic methods are more practical and useful for such a deep investigation. We consider audio-frequency magnetotelluric (AMT) and surface-to-borehole electromagnetic (EM) tomography methods as promising tools for the investigation of deep aquifer. These methods were tested in the Hasunuma area, Chiba Prefecture, Japan. Although the study area is in an urban area, high-quality AMT data were acquired, which was mainly accomplished by night-time data recording and remote-reference data processing. One-dimensional inversion results of the AMT data revealed two extremely conductive zones, which is consistent with the electrical conductivity profile of pore water in core samples. It can be interpreted as the seawater intrusions into both zones. However, the chemical analysis of the groundwater sampled in the deep zone suggests that this groundwater must be fossil seawater that had been confined during sedimentation processes. In addition, the permeability coefficient of the deep layer is very low. Thus the deep conductive zone corresponds to the fossil seawater regarded as being difficult to flow.

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

In order to increase the signal-to-noise ratio in multi-component seismic data, we developed new polarization filters based on the method of multicomponent complex trace analysis. Unlike the previous polarization filters, the present filters separately compute linear and elliptic components at each time sample using amplitude ratio of horizontal and vertical components of body waves and ellipticity of Rayleigh waves. The polarization filters work ideally even with low S/N data. Application of the filters to both synthetic and real seismic data shows that Rayleigh waves of elliptic motions are effectively eliminated and both P and S waves of linear motions are well separated each other.

• Proceedings of the KSRS Conference
• /
• 2008.10a
• /
• pp.414-417
• /
• 2008

We have extracted a surface deformation map of a part of Canisteo Peninsula on Amundsen Sea in West Antarctica by applying 4-pass DInSAR technique to two ERS-1/2 tandem pairs obtained on October 21-22, 1995 (diff-pair) and March 9-10, 1996 (topo-pair), and analyzed changes of glaciers, sea ice, ice shelves, and their kinematic interactions. We observed fast motion of glaciers pushing the adjoining sea ice. Some interferometric phases indicate the up-rise of sea ice of which type is thought to be land-fast ice to exert repulsive force against the pushing glacier. There were other glaciers and sea ice that moved to the same direction, suggesting that the sea ice in these regions was land-fast ice weakly harnessed to sea bottom or pack ice not harnessed at all. Several small circular fringes in ice shelves suggested that islands or seamounts on the bottom of ice shelves deterred the movement of ice shelves, resulting in the rise of ice surface.

• Proceedings of the KSRS Conference
• /
• 2008.10a
• /
• pp.422-425
• /
• 2008

TerraSAR-X is new radar satellite operated at X-band, multi polarization, and multi beam mode. Compared with C-band or L-band SAR, the X-band system inherently suffers from more temporal decorrelation, but is more sensitive to surface deformation monitoring due to short wavelength (3.1 cm) and high spatial resolution (1m-3m). It is generally expected that sensitivity to estimate surface movement using TerraSAR-X will be increased by the factor of 10, compared to current C-band system with low spatial resolution such as ERS-2, Envisat. Many urban areas are experiencing land subsidence due to water, oil and natural gas withdrawal, underground excavation, sediment compaction, and so on. Monitoring of surface deformation is valuable for effectively limiting damage areas. In addition high accuracy and spatially dense subsidence map can be achieved by X-band InSAR observation, promoting identification and separation of various subsidence processes and leading to enhanced understanding via mechanical modeling. In this study we will introduce some initial InSAR results using new TerraSAR-X SAR data for surface deformation monitoring.

• International Union of Geodesy and Geophysics Korean Journal of Geophysical Research
• /
• v.26 no.1
• /
• pp.1-14
• /
• 1998

Various numerical methods for the two dimensional shallow water equations have been applied to the problems of flood routing, tidal circulation, storm surges, and atmospheric circulation. These methods are often based on the Alternating Direction Implicity(ADI) method. However, the ADI method results in inaccuracies for large time steps when dealing with a complex geometry or bathymetry. Since this method reduces the performance considerably, a fully implicit method developed by Wilders et al. (1998) is used to improve the accuracy for a large time step. Finite Difference Methods are defined on a rectangular grid. Two drawbacks of this type of grid are that grid refinement is not possibile locally and that the physical boundary is sometimes poorly represented by the numerical model boundary. Because of the second deficiency several purely numerical boundary effects can be involved. A boundary fitted curvilinear coordinate transformation is used to reduce these difficulties. It the curvilinear coordinate transformation is used to reduce these difficulties. If the coordinate transformation is orthogonal then the transformed shallow water equations are similar to the original equations. Therefore, an orthogonal coorinate transformation is used for defining coordinate system. A multigrid (MG) method is widely used to accelerate the convergence in the numerical methods. In this study, a technique using a MG method is proposed to reduce the computing time and to improve the accuracy for the orthogonal to reduce the computing time and to improve the accuracy for the orthogonal grid generation and the solutions of the shallow water equations.

• Geophysics and Geophysical Exploration
• /
• v.1 no.3
• /
• pp.155-160
• /
• 1998

EM tomography technique has been developed. The algorithm used the extended Born approximations for forward modeling and reconstructed a conductivity image by a smoothness constraint least squares inversion method. Observed data, the vertical components of secondary magnetic fields, were simulated with the 3-D integral equation code. The results showed that the location of anomalous body could be imaged very well, but conductivity of the body was lower than real one and the vertical resolution was much higher than the horizontal resolution.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2007.06a
• /
• pp.139-144
• /
• 2007

To delineate subglacial and englacial images of the western Fourcade Glacier in Antarctica, GPR data were acquired along a 470 m profile in November 2006. Signature deconvolution, migration velocity analysis, and finite difference depth migration were effective in increasing vertical resolution, obtaining the velocity function, and yielding clear depth images of the mixed-phase single-channel data, respectively. The source signature was extracted through stacking reflection signals from the basement. The deconvolution successfully compressed the source signatures to increase temporal resolution and attenuated reverberations. Medium velocities were analyzed with 112 diffraction events by the constant-velocity migration method based on the Stolt algorithm. Finally, depth migration based on the finite difference algorithm provided clear englacial and subglacial images in the laterally-varying velocity field.

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

A preliminary study was made for density correction for the effect of separation between sonde and borehole wall due to the existence of casing in a borehole. Firstly, series of gamma-gamma measurements were obtained with 3 different separations in 3 different density models constructed in Kangwon National University. Secondly, density correction charts were made based on the 'spine and ribs' technique in mud cake correction. Using these charts, we could determine effectively the true density from the measurement in the thinly cased borehole.

• The Journal of Engineering Geology
• /
• v.18 no.2
• /
• pp.137-144
• /
• 2008

The stand-off errors due to the different separations between the sonde and the borehole wall were measured and analyzed in 4 physical borehole models located in Kangwon National University having different densities with cylindrical and half cylindrical PVC and/or acrylic casings. The analysis of the stand-off error data based on the "spine and ribs" technique suggests a well defined rib line for each model irrespectively of the types and thicknesses of the casing, and that the gradients of the ribs are proportional to the densities of the models. By using these characteristics successful density correction could be made for the plastic casings in NX sized boreholes.