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

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 for the travel time computation commonly used, however, may not give sufficiently precise results 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. 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.

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.335-358
• /
• 2020

We present a numerical model to simulate coupled hydro-mechanical behavior of fault using TOUGH-FLAC simulator. This study aims to develop a numerical method to estimate fluid injection-induced fault reactivation in low permeability rock and to access the relevant hydro-mechanical stability in rock as part of DECOVALEX-2019 Task B. A coupled fluid flow and mechanical interface model to explicitly represent a fault was suggested and validated from the applications to benchmark simulations and the field experiment at Mont Terri underground laboratory in Switzerland. The pressure build-up, hydraulic aperture evolution, displacement, and stress responses matched those obtained at the site, which indicates the capability of the model to appropriately capture the hydro-mechanical processes in rock fault.

• 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.

• Tunnel and Underground Space
• /
• v.23 no.4
• /
• pp.261-270
• /
• 2013

As part of the development of the 3-D geologic modeling software, this study addresses on new development of software modules that can perform the analysis and visualization of the fracture network system in 3-D. The developed software modules, such as BOUNDARY, DISK3D, FNTWK3D, CSECT and BDM, are coded on Microsoft Visual Studio platform using the MFC and OpenGL library supported by C++ program language. Each module plays a role in construction of analysis domain, visualization of fracture geometry in 3-D, calculation of equivalent pipes, production of cross-section map and management of borehole data, respectively. The developed software modules for analysis and visualization of the 3-D fracture network system can be used to tackle the geomechanical problems related to strength, deformability and hydraulic behaviors of the fractured rock masses. All these benefits will further enhance the economic competitiveness of the domestic software industry.

• Journal of the Korean earth science society
• /
• v.26 no.7
• /
• pp.691-697
• /
• 2005

Recently, the development of accurate gravity-meter and GPS make it possible to obtain high resolution gravity data. Though gravity data interpretation like modeling and inversion has significantly improved, gravity data processing itself has improved very little. Conventional gravity data processing removes gravity effects due to mass and height difference between base and measurement level. But, it would be a biased density model when some or whole part of anomalous bodies exist above the base level. We attempted to make a multiquadric surface of the survey area from topography with DEM (Digital Elevation Map) data. Then we constituted rectangular blocks which reflect real topography of the survey area by the multiquadric surface. Thus, we were able to carry out 3-D inversions which include information of topography. We named this technique, 3-D Gravity Terrain Inversion (3DGTI). The model test showed that the inversion model from 3DGTI made better results than conventional methods. Furthermore, the 3-dimensional model from the 3DGTI method could maintain topography and as a result, it showed more realistic geologic model. This method was also applied on real field data in Masan-Changwon area. Granitic intrusion is an important geologic characteristic in this area. This method showed more critical geological boundaries than other conventional methods. Therefore, we concluded that in the case of various rocks and rugged terrain, this new method will make better model than convention ones.

• Economic and Environmental Geology
• /
• v.49 no.4
• /
• pp.291-299
• /
• 2016

This study examined stratigraphic research containing extreme climate event during Quaternary period in Saban-ri, Haeri-myeon, Gochang by constructing 3D topographic model and 3D geological model. As a result of 3D topographic model and subsurface geological model, the geology of study area accumulated bedrock, Pleistocene series, and Holocene series chronologically. Most of the study area consist of bedrock on basement and Holocene series on upper layer. Additionally, Pleistocene series are presented as lens-shaped deposit on eastern part, and wedge-shaped deposit on northeastern part. Holocene layers consist of sand and clay-silt layer deposited sequentially where implies fluvial deposits on transgression environment. Distinctively, Pleistocene clayey silt layer and Holocene sand layer on eastern are observed as pond shape deposits that are considered as storm-related deposits originated from overwash system caused by extreme paleoclimate.

• Geophysics and Geophysical Exploration
• /
• v.27 no.1
• /
• pp.37-50
• /
• 2024

The magnetotelluric (MT) survey can be affected by external environmental factors. In particular, when acquiring MT data in islands, it is essential to consider the combined effect of topography and sea to understand the results and make accurate interpretations. To analyze the MT response (apparent resistivity, phase) with consideration of the effect of topography and sea, a small cone-shaped island model surrounded by deep sea was created. Two-dimensional (2-D) and three-dimensional (3-D) forward modeling were performed on the terrain model considering topography and the island model considering both topography and sea. The 2-D MT response did not reflect the topographic and sea effect of the direction orthogonal to the 2-D profile. The 3-D MT response included topographic and sea effects in all directions. The XY and YX components of the apparent resistivity were separated on undulating topography, such as a hill. A conductor at 1 km below sea level could be distinguished from topographic and sea effects in the MT response, and low resistivity anomaly was attenuated at greater depths. This study will facilitate understanding of field data measured on small islands.