• Journal of Soil and Groundwater Environment
• /
• v.9 no.1
• /
• pp.79-86
• /
• 2004

Hydraulic conductivity along rock fracture is mainly dependent on fracture geometries such as orientation, aperture, roughness and connectivity. Therefore, it needs to consider fracture geometries sufficiently on a fracture model for a numerical analysis to calculate permeability coefficient in a fracture. This study performed new type of numerical analysis using a homogenization analysis method to calculate permeability coefficient accurately along single fractures with several fracture models that were considered fracture geometries as much as possible. First of all, fracture roughness and aperture variation due to normal stress applied on a fracture were directly measured under a confocal laser scaning microscope (CLSM). The acquired geometric data were used as input data to construct fracture models for the homogenization analysis (HA). Using the constructed fracture models, the homogenization analysis method can compute permeability coefficient with consideration of material properties both in microscale and in macroscale. The HA is a new type of perturbation theory developed to characterize the behavior of a micro inhomogeneous material with a periodic microstructure. It calculates micro scale permeability coefficient at homogeneous microscale, and then, computes a homogenized permeability coefficient (C-permeability coefficient) at macro scale. Therefore, it is possible to analyze accurate characteristics of permeability reflected with local effect of facture geometry. Several computations of the HA were conducted to prove validity of the HA results compared with the empirical equations of permeability in the previous studies using the constructed 2-D fracture models. The model can be classified into a parallel plate model that has fracture roughness and identical aperture along a fracture. According to the computation results, the conventional C-permeability coefficients have values in the range of the same order or difference of one order from the permeability coefficients calculated by an empirical equation. It means that the HA result is valid to calculate permeability coefficient along a fracture. However, it should be noted that C-permeability coefficient is more accurate result than the preexisting equations of permeability calculation, because the HA considers permeability characteristics of locally inhomogeneous fracture geometries and material properties both in microscale and macroscale.

• Journal of the Korea Institute of Information and Communication Engineering
• /
• v.7 no.3
• /
• pp.478-484
• /
• 2003

The electromagnetic wave scattering problems from inhomogeneous material bodies are considered. The formulation is made in terms of mixed potentials for the moment methods (MoM). The surfaces of a three-dimensional inhomogeneous scatterer of arbitrary shape are divide into triangular patches for descretization. Application of the boundary conditions leads to the coupled surface integral equations to be satisfied for the unknown surface equivalent electric and magnetic currents. The radar cross-section (RCS) for some structures is computed and the results are compared with the reported data.

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

Most of seismic reflection prospecting assumes subsurface formation to be homogeneous media. These models are not capable of estimating small scale heterogeneity which is verified by well log data or drilling core. And those synthetic seismograms by homogeneous media are limited to explain various changes at field data. So we developed a inhomogeneous velocity model which can estimate inhomogeneity of background medium to implement numerical modeling from homogeneous medium and inhomogeneous medium on the model. Background medium using three autocorrelation functions in order to generate inhomogeneous velocity media was according to dominant wavelength of background medium and correlation length of random medium. And then we compared shot gathers. The results show that numerical modeling implemented at inhomogeneous medium depicts complex wave propagation of field data.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.515-515
• /
• 2017

하천에서 취수원을 개발하는데 있어 지층의 구성으로 물리, 화학적 여과, 흡착 등을 통해 자연 정화되는 간접 취수 방식이 활발히 도입되고 있다. 양질의 취수원을 공급할 수 있는 간접 취수 방식은 수량 확보 측면에서의 불확실성과 유지관리상의 어려움 때문에 많은 시행착오가 발생된다. 이와 같은 단점을 개선하기 위해 하상을 개착하여 불균질한 대수층을 치환하고 스크린을 통해 간접 취수원을 개발하는 하상여과 방식이 도입되고 있다. 대수층을 치환하여 여재를 구성함에 있어 오염물질 및 탁도의 저감 효율을 극대화하기 위해 다양한 연구가 진행되고 있다. 본 연구에서는 개착식 하상여과에서 치환하는 여재의 구성에 따른 탁도의 저감 효율을 분석하기 위해 축소 모형실험으로 구성하였다. 각각의 여재의 구성에 대해서는 상수도 시설기준을 통해 축소된 입경의 매질로 구성하였다. 실험실 규모의 모형 수조($1500mm{\times}500mm{\times}1700mm$)를 구성하고 하부에는 내경 80mm이고 길이 1300mm인 기능성 스크린이 부착된 취수관을 설치하였다. 모형 수조에서 여재의 두께는 총 1000mm로 구성하였고, 각각의 층에 대해서는 250mm로 하여 4개의 층을 구성할 수 있도록 하였다. 치환하는 여재의 매질에는 자갈, 왕사, 중사, 화산석을 사용하였고, 각각의 입경은 5-10mm, 2-5mm, 1-2mm, 2-5mm이고, 탁도를 유발하는 물질로는 입경이 $20{\mu}m$인 황토를 사용하였다. 단일매질 구성을 통해 각각의 여재 종류에 따른 탁도 저감 효과에 대해서 분석하였고, 세 가지의 혼합매질 구성을 통해 치환층 여재의 배치에 따른 탁도 저감 효과를 분석하였고, 각각의 구성은 중사-왕사-자갈-화산석, 화산석-중사-왕사-자갈, 중사-왕사-화산석-자갈로 하였다. 주입수는 30-50NTU를 유지하였으며 유출수의 탁도를 통해 저감 효율을 분석하였다. 분석된 결과를 통해 개착식 하상여과 방식의 여재 구성에 대해서 탁도 저감에 효율을 극대화할 수 있을 것으로 보여진다.

• Geophysics and Geophysical Exploration
• /
• v.9 no.3
• /
• pp.241-249
• /
• 2006

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms commonly used, however, may not give sufficiently precise computational results of traveltime data particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. Considering the complex geology of Korea, we assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution. The performance of the algorithm developed in this study is demonstrated by the comparison of the analytic and numerical solutions for the homogeneous anisotropic earth as well as through the numerical experiment for the two layer model whose anisotropic properties are greatly different each other. We expect that the developed modeling algorithm can be used in the development of processing and inversion schemes of seismic data acquired in strongly anisotropic environment, such as migration, velocity analysis, cross-well tomography and so on.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.229-240
• /
• 2007

In this study, seismic anisotropic tomography algorithm was developed for imaging the seismic velocity anisotropy of the subsurface. This algorithm includes several inversion schemes in order to make the inversion process stable and robust. First of all, the set of the inversion parameters is limited to one slowness, two ratios of slowness and one direction of the anisotropy symmetric axis. The ranges of the inversion parameters are localized by the pseudobeta transform to obtain the reasonable inversion results and the inversion constraints are controlled efficiently by ACB(Active Constraint Balancing) method. Especially, the inversion using the Fresnel volume is applied to the anisotropic tomography and it can make the anisotropic tomography more stable than ray tomography as it widens the propagation angle coverage. The algorithm of anisotropic tomography is verified through the numerical experiments. And, it is applied to the real field data measured at limestone region and the results are discussed with the drill log and geological survey data. The anisotropic tomography algorithm will be able to provide the useful tool to evaluate and understand the geological structure of the subsurface more reasonably with the anisotropic characteristics.

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

In Korean geology that crystalline rock is dominant, the properties of subsurface including the anisotropy are distributed complexly and changed abruptly. Because of such geological environments, cross-hole seismic traveltime tomography is widely used to obtain the high resolution image of the subsurface for the engineering purposes in the geotechnical sites. However, because the cross-hole tomography has a wide propagation angle coverage relatively, its data tend to include the seismic velocity anisotropy comparing with the surface seismic methods. It can cause the misinterpretation that the cross-hole seismic data including the anisotropic effects are analyzed and treated with the general processing techniques assuming the isotropy. Therefore, we need to consider the seismic anisotropy in cross-hole seismic traveltime tomography. The seismic anisotropic tomography algorithm, which is developed for evaluation of the velocity anisotropy, includes several inversion schemes in order to make the inversion process stable and robust. First of all, the set of the inversion parameters is limited to one slowness, two ratios of slowness and one direction of the anisotropy symmetric axis. The ranges of the inversion parameters are localized by the pseudo-beta transform to obtain the reasonable inversion results and the inversion constraints are controlled efficiently by ACB(Active Constraint Balancing) method. Especially, the inversion using the Fresnel volume is applied to the anisotropic tomography and it can make the anisotropic tomography more stable than ray tomography as it widens the propagation angle coverage.

• Journal of the Korean Geophysical Society
• /
• v.6 no.4
• /
• pp.231-244
• /
• 2003

Lots of various utilities are buried under the surface of the earth. The effective handling of the underground utilities is becoming the big subject and project for the harmonious management and administration of the city. To detect the position and depth of buried underground utilities, GPR and Induced EM surveys are commonly used. However, they have limitations, such as shallow skin-depth and non-availability in the areas where subsurface materials are not homogeneous and are compose of clays and/or salts and gravels. The aim of this study is to find the efficient geophysical method which can overcome these limitations. For this purpose, various geophysical mehods were applied in the site of poor geotechnical environment.

• Economic and Environmental Geology
• /
• v.53 no.6
• /
• pp.703-713
• /
• 2020

In many geological fields, including hydrogeology, it is of great importance to determine the heterogeneity of the subsurface media. This study briefly introduces the concept and theory of the method that can estimate the hydraulic properties of the media constituting the aquifer, which was recently introduced by Park (2020). After the introduction, the method was applied to the synthetic aquifer to demonstrate the practicality, from which various implications were drawn. The introduced technique uses a global optimization technique called the covariance matrix adaptation evolution strategy (CMA-ES). Conceptually, it is a methodology to characterize the aquifer heterogeneity by assimilating the groundwater level time-series data due to the imposed hydraulic stress. As a result of applying the developed technique to estimate the hydraulic conductivity of a hypothetical aquifer, it was confirmed that a total of 40000 unknown values were estimated in an affordable computational time. In addition, the results of the estimates showed a close numerical and structural similarity to the reference hydraulic conductivity field, confirming that the quality of the estimation by the proposed method is high. In this study, the developed method was applied to a limited case, but it is expected that it can be applied to a wider variety of cases through additional development of the method. The development technique has the potential to be applied not only to the field of hydrogeology, but also to various fields of geology and geophysics. Further development of the method is currently underway.

• Proceedings of the International Union of Geodesy And Geophysics Korea Journal of Geophysical Research Conference
• /
• 2003.05a
• /
• pp.16-16
• /
• 2003

본 연구는 ILP(International Lithosphere Project) Task Group II-4가 진행하고 있는 상부맨틀에 대한 3차원 구조도 작성 연구의 일환으로 수행되어졌으며 구조도 작성을 위한 데이터 베이스의 구조는 task group의 표준안을 따랐다. 기존 문헌과 기존의 데이터 베이스를 통해서 획득된 자료를 이용해 동해와 그 주변을 대상으로 하는 지역의 ($32-45^{\circ}$E, $122-148^{\circ}$N) 상부 670km까지의 3차원 구조도 작성을 위한 초기 모델을 구축하였으며, 이 절차를 최대한 자동화하는 프로그램을 포트란을 이용해 만들어보았다. 연구 지역에 대한 곡율을 계산하기 위해 표준타원체 모델인 WGS84과 geoid undulation 모델인 EGM96을 사용했으며 지형 고도 자료는 GTOPO30, GLOBE 1.0, 그리고 Smith and Sandwell 데이터베이스를 사용하여 지구 중심으로부터 지표까지의 거리를 구하였다. 연구지역은 $0.25^{\circ}$간격으로 나누었으며 총 5777개의 격자점을 정의하였으며 각각의 격자점에 1차원 수직구조를 부여함으로써 3차원 모델을 구축하였다. 그리고 지형적으로나 지질학적으로 유사한 지역을 하나의 구역으로 정의하고 동일한 수직구조를 부여함으로써 모든 격자점에 1차원 수직구조를 정의하지 않도록 하였다. 본 연구에서는 지표 지질은 모델에 고려하지 않았지만 지형학적으로 의미가 있는 분지나 수평적으로 불균질성이 뚜렷한 지역을 중심으로 연구 지역의 리젼을 정의하였다. 중요 리젼에 대한 지각구조에 대해서는 기존의 문헌을 통해 모델치를 정의하였으며 지각 하부부터 상부 670km에 대한 부분은 Task Group에서 제시한 표준 모델을 이용했다. 모델을 정의하기 위해 주어진 격자점에 대한 지구 중심으로부터 지오이드까지의 거리, 지오이드로부터 지표까지의 거리를 정의해주었으며, 각 격자점의 수직구조를 정의하기 위해 깊이에 따른 각 매질의 밀도, P파의 속도, S파의 속도, P파에 대한 Q값, S파에 대한 Q값을 정의 해주었다. S파의 속도를 구하기 위해서 지구 내부 물질을 포아송 매질이라는 가정 하에, 관계식을 $Vp{\;}={\;}SQRT(3){\;}{\times}{\;}Vs$ 이용하였다. 획득한 모델치들을 이용해 동해와 동해 인근 지역에 대한 초기모델을 구축하였다.