POPs in sediments and soil in the PRD are comparable to or much higher than those reported in other regions. Some sites may be classified as POPs- polluted with high ecological risks. Large-scale land transform in the process of regional urbanization may facilitate the transfer of POPs in the soil to the sedimentary system by enhancing the soil run-off. Urban atmospheric PCBs in PRD are found to be less than some of the North American or European urbans, but PAHs are significantly higher. The center of the PRD has been the major source area of PAHs and organochlorine pesticides in the PRD. The northern part of the PRD serves as a regional sink for the air particulates and affiliated POPs.

Reverse-time migration has no dip limitations and one of the most promising methods to preserve true amplitudes. We applied 3-D prestack reverse time migration based on a pseudo-spectral implementation of the acoustic wave equation to the SEG/EAGE salt dome synthetic data set. We were able to illuminate sub salt reflectors of the SEG/EAGE salt model that were barely observable in the Kirchhoff migration images. Using the pseudo-spectral modeling technique, we could implement reverse-time migration within the core memory, which could be equipped to a personal computer.

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.

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.

The extended Born, or localized nonlinear approximation of integral equation (IE) solution has been applied to inverting single-hole electromagnetic (EM) data using a cylindrically symmetric model. The extended Born approximation is less accurate than a full solution but much superior to the simple Born approximation. When applied to the cylindrically symmetric model with a vertical magnetic dipole source, however, the accuracy of the extended Born approximation is greatly improved because the electric field is scalar and continuous everywhere. One of the most important steps in the inversion is the selection of a proper regularization parameter for stability. Occam's inversion (Constable et al., 1987) is an excellent method for obtaining a stable inverse solution. It is extremely slow when combined with a differential equation method because many forward simulations are needed but suitable for the extended Born solution because the Green's functions, the most time consuming part in IE methods, are repeatedly re-usable throughout the inversion. In addition, the If formulation also readily contains a sensitivity matrix, which can be revised at each iteration at little expense. The inversion algorithm developed in this study is quite stable and fast even if the optimum regularization parameter Is sought at each iteration step. Tn this paper we show inversion results using synthetic data obtained from a finite-element method and field data as well.

경제적이며 정확한 3차원 전자탐사 모델링을 위해 위해 Habashy et al. (1993)에 의해 제안된 국소 비선형 근사(localized nonlinear approximation)를 이용하여 전자탐사 모델링 알고리듬을 개발하였다. 전자탐사 수치모델링시 많은 계산시간 및 기억용량을 필요로 하는 Green 텐서 적분을 정확하고 빠르게 계산하기 위해, 단일 미소체를 이용한 공간파수 영역에서의 Green 텐서 적분 알고리듬을 제안하였다. 더욱이 Green 텐서의 송수신 방향 및 상반성을 고려하여 각각의 미소체에 의한 전체 미소체에의 Green 텐서 적분을 한 개의 미소체에 의한 전체 미소체에의 Green 텐서 적분 값으로 구하게 하므로 매우 적은 기억용량 만으로 Green 텐서 적분 행렬을 구성할 수 있어, 역산법에 효과적으로 적용할 수 있다. 이 수치 모델링 알고리듬을 기본으로 하여 평활화 제한을 가한 최소자승 역산 알고리듬을 개발하였다. 이 역산 알고리듬을 지표 전자탐사 및 시추공-지표 전자탐사 등에 적용하여 PC에서도 빠르게 3차원 전자탐사 역산이 수행됨을 보였다.

We evaluated the potential of the limestone deposits distributed in the Uljin area, the north-eastern part of Gyungsang province, based on the results of the surface and underground geological survey. The geology of the Uljin area consists mainly of the Cambro-Ordovician Joseon Supergroup correlated to the Duwibong sequence of the south-eastern part of Gangwon province. The Joseon Supergroup rested unconformably on the Precambrian schists is divided into six formations , namely, Maesan, Gusanri, Namsusan, Songchon, Mongcheon, and Keummaeri formations in ascending order. Main geological structures are two thrust faults with NS trending, some normal and reverse faults, and folds.

연구지역은 좁은 면적에 비하여 800m 이상의 지형기복을 보이며, 계곡과 하천 충적층을 따라 마을이 발달된 전형적인 수계 상류의 농촌지역이다. 주요 상수원은 천부지하수이며, 농업용으로는 하천수를 이용한다. 지표수는 연구지역의 소유역(A, B)간 계절별 수질차이를 보이나, 지하수의 경우 유역간 및 계절별 변동이 상대적으로 작게 나타난다. B 소유역의 하류부에서는 지질매체의 영향으로 보이는 Fe, Mn, Zn 농도의 부화현상이 나타났으며, 이는 방이리층과 길왕리층을 기반암으로 하는 지역에 대한 전체적인 지하수질 조사가 필요함을 지시한다. 국지적인 수질오염원으로 과거에 개발되었던 폐금광이 존재하며, 이로 인한 산성수의 발생도 인지되었다.