• Geophysics and Geophysical Exploration
• /
• v.5 no.3
• /
• pp.206-216
• /
• 2002

The finite element method (FEM), a powerful numerical modeling tool for solving various engineering problems, is frequently applied to three-dimensional (3-D) modeling thanks to its capability of discretizing and simulating the shape of model with finite number of elements. Considering the accuracy of the solution and computing time in modeling of engineering problems, it is preferable to construct physical continuity and simplify mesh system. Although there exist systematic mesh generation systems for arbitrary shaped model, it is hard to model a simple cylinder in terms of 3-D coordinate system especially in the vicinity of the central axis. In this study I adopt cylindrical coordinate system for modeling the 3-D model space and define the origin of the coordinates with mathematically clear coordinate transformation. Since we can simulate the whole space with hexahedral elements, the cylindrical coordinate system is effective in handling the 3-D model structure. The 3-D do resistivity modeling scheme developed in this study provides basie principle for borehole-to-surface resistivity survey, which can be a useful tool for the application to environmental problem.

• Geophysics and Geophysical Exploration
• /
• v.8 no.1
• /
• pp.26-33
• /
• 2005

Inversion of multi-mode surface-wave phase velocity for shallow engineering site investigation has received much attention in recent years. A sensitivity analysis and inversion of both synthetic and field data demonstrates the greater effectiveness of this method over employing the fundamental mode alone. Perturbation of thickness and shear-wave velocity parameters in multi-modal Rayleigh wave phase velocities revealed that the sensitivities of higher modes: (a) concentrate in different frequency bands, and (b) are greater than the fundamental mode for deeper parameters. These observations suggest that multi-mode phase velocity inversion can provide better parameter discrimination and imaging of deep structure, especially with a velocity reversal, than can inversion of fundamental mode data alone. An inversion of the theoretical phase velocities in a model with a low velocity layer at 20 m depth can only image the soft layer when the first higher mode is incorporated. This is especially important when the lowest measurable frequency is only 6 Hz. Field tests were conducted at sites surveyed by borehole and PS logging. At the first site, an array microtremor survey, often used for deep geological surveying in Japan, was used to survey the soil down to 35 m depth. At the second site, linear multichannel spreads with a sledgehammer source were recorded, for an investigation down to 12 m depth. The f-k power spectrum method was applied for dispersion analysis, and velocities up to the second higher mode were observed in each test. The multi-mode inversion results agree well with PS logs, but models estimated from the fundamental mode alone show f large underestimation of the depth to shallow soft layers below artificial fill.

• Korean Journal of Remote Sensing
• /
• v.29 no.5
• /
• pp.443-459
• /
• 2013

Located on Java subduction zone, Merapi volcano is an active stratovolcano with a volcanic activity cycle of 1-5 years. Merapi's eruptions were relatively small with VEI 1-3. However, the most recent eruption occurred in 2010 was quite violent with VEI 4 and 386 people were killed. In this study, we have attempted to study the characteristics of Merapi's eruptions during 18 years using optical Landsat images. We have collected a total of 55 Landsat images acquired from July 6, 1994 to September 1, 2012 to identify pyroclastic flows and their temporal changes from false color images. To extract areal extents of pyroclastic flows, we have performed supervised classification after atmospheric correction by using COST model. As a result, the extracted dimensions of pyroclastic flows are nearly identical to the CVP monthly reports. We have converted the thermal band of Landsat TM and ETM+ to the surface temperature using NASA empirical formula and calculated time-series of the mean surface temperature in the area of peak temperature surrounding the crater. The mean surface temperature around the crater repeatedly showed the tendency to rapidly rise before eruptions and cool down after eruptions. Although Landsat satellite images had some limitations due to weather conditions, these images were useful tool to observe the precursor changes in surface temperature before eruptions and map the pyroclastic flow deposits after eruptions at Merapi volcano.

• Journal of the Mineralogical Society of Korea
• /
• v.17 no.3
• /
• pp.245-265
• /
• 2004

This study was to compare the geological occurrences and geneses of the Myogi, Tsukinuno, Dobuyama and Kawasaki bentonite deposits distributed in the Tertiary sedimentary basins of NE Japan, and to compare the mineralogical and physicochemical properties of their bentonites. The Japanese bentonite deposits are mainly distributed in the Green-tuff region which was formed in Neogene. The shape of ore body of the Myogi, Tsukinuno and Kawasaki deposits formed by the diagenesis are layered and stratiform. In contrast to this, the Dobuyama deposit formed by hydrothermal alteration shows the cone shape. The mineralization age of four deposits are 1.8 ～ 21 Ha from Early Miocene to Pliocene. The Dobuyama bentonite with the highest montmorillonite content shows the highest surface area, CEC, MB adsorption, and strengths. The Tsukinuno bentonite with a little high montmorillonite content is characterized by strong alkalinity, high viscosity and swelling. The Kawasaki bentonite, the Na-Ca mixed type, shows higher viscosity and swelling than the Ca-type Dobuyama bentonite. The Myogi bentonite with the lowest montmorillonite content shows the properties of low viscosity, In adsorption, strengths and a little high CEC and surface area. The high CEC and surface area of this deposit is due to the sufficient occurrence of zeolite. A strong dispersion in the Na-type bentonite and a strong flocculation in the Ca-type bentonite took place, and both the types show a slow flocculation with time. The physicochemical properties of the bentonite are mainly controlled by the montmorillonite content, interlayer cations, and impurity minerals such as zeolite. But bentonites inconsistent to this factors are sometimes occurred. This is maybe due to the crystal chemistry such as layer charge of montmorillonite and crystal morphology of montmorillonite such as aspect ratio.

• Korean Journal of Remote Sensing
• /
• v.27 no.2
• /
• pp.171-180
• /
• 2011

MAI (multiple aperture SAR interferometry) method has been recently developed to improve the measurement accuracy of along-track surface deformation. By means of split-beam SAR processing, this novel technique produces forward- and backward-looking interferograms, which are combined to generate an MAI interferogram. The along-track surface deformation can then be derived from the MAI interferogram. The achieved accuracy of the along-track surface deformation is approximately 8 cm for interferograms with a coherence of 0.6. It is commonly recognized that the topographic phase on an MAI interferogram can be ignored. However, in this paper, we have generated an MAI interferogram from an ALOS P ALSAR interferometric pair spanning the 2010 Haiti earthquake, and shown that the topographic phase distortion on the MAI interferogram can reach to about $3.45{\times}10^{-4}$ rad./m. This distortion corresponds to an along-track surface deformation of about 98 cm. We have proposed an efficient method to remove the topographic phase distortion. After correcting the distortion, the topographic phase distortion on the MAI interferogram is reduced to about $7.82{\times}10^{-6}$ rad./m. This means that the proposed method can effectively remove the topographic distortion on the MAI interferogram to improve along-track surface deformation measurement.

• Geophysics and Geophysical Exploration
• /
• v.13 no.1
• /
• pp.61-68
• /
• 2010

A magnetotelluric (MT) survey has been performed to delineate deeply extended fracture systems at the geothermal field in Seokmo Island, Korea. To assist interpretation of the MT data, geological surveying and well logging of existing wells were also performed. The surface geology of the island shows Cretaceous and Jurassic granite in the north and Precambrian schist in the south. The geothermal regime has been found along the boundary between the schist and Cretaceous granite. Because of the deep circulation along the fracture system, geothermal gradient of the target area exceeds $45^{\circ}C/km$, which is much higher than the average geothermal gradient in Korea. 2D and 3D inversions of MT data clearly showed a very conductive anomaly, which is interpreted as a fracture system bearing saline water that extends at least down to 1.5 km depth and is inclined eastwards. After drilling down to the depth of 1280 m, more than 4000 tons/day of geothermal water overflowed with temperature higher than $70^{\circ}C$. This water showed very similar chemical composition and temperature to those from another existing well, so that they can be considered to have the same origin; i.e. from the same fracture system. A new geothermal project for combined heat and power generation was launched in 2009 in Seokmo Island, based on the survey. Additional geophysical investigations including MT surveys to cover a wider area, seismic reflection surveys, borehole surveys, and well logging of more than 20 existing boreholes will be conducted.

• The Journal of Engineering Geology
• /
• v.18 no.2
• /
• pp.127-136
• /
• 2008

Geophysical survey for establishing a wide site for the distribution of water content, wetting front infiltration due to the rainfall, and distribution of groundwater level has been performed by using 8round penetration radar (GPR) method, electrical resistivity method, and so on. On the other hand, a narrow area survey was performed to use a permittivity method such as time domain reflectometry, frequency domain reflectometry, and amplitude domain reflectometry methods for estimating volumetric water content, soil density, and concentration of contaminant in surface and subsurface. The permittivity methods establish more corrective physical parameters than different found survey technologies mentioned above. In this study for establishment of infiltration behaviors for wetting front in the unsaturated soil caused by an artificial rainfall, soil physical parameters for volumetric water content, pore water pressure, and pore air pressure were measured by FDR measurement device and pore water pressure meter which are installed in the unsaturated weathered granite soil with different depths. Consequently, the authors were proposed to a new establishment method for analyzing the variations of volumetric water content and wetting front infiltration from the responses of infiltrating pore water in the unsaturated soil.

• Journal of Korea Soil Environment Society
• /
• v.1 no.2
• /
• pp.15-21
• /
• 1996

Effects of acid deposition on the acidification of surface water over a short term period were studied during June~October, 1995. A reservoir located in Chunchon, Kangwondo was selected for this study because this reservoir had received little man-made pollutants. Geological survey within the study area was carried out and the pH of a small and shallow stream flowing from the mountain valley down to the reservoir was monitored. Rainfall was measured by using a rain gauge and the precipitation was collected by automatic precipitation sampler. During the study period, rain pH was in the range of 3.81~5.77 with an average of 4.8. The EC was in the range of 5~189$\mu$S/cm with an average of 10.6$\mu$S/cm. ionic composition of the rainfall showed that of the anions, deposition of ${SO_4}^{2-}$, was highest with 3,119.7 kg/$\textrm{km}^2$, and it was $NH^{4+}$ with 1,053.2kg/$\textrm{km}^2$ for cations. The surface water pHs were maintained neutral or weak basic representing little evidence of acidification regardless of the acidic rainfall pH. Every time of sampling, however, the pHs exhibited increasing tendency as elevation of measuring site of the stream became lowered, which indicated the impacts of acid-base reactions by acidic or basic substances during travel of water downstream. This result suggested the necessity of further research to determine the respective buffering capacity of soil and underlying rocks.

• The Journal of Engineering Geology
• /
• v.23 no.2
• /
• pp.149-160
• /
• 2013

Volcanic activity commonly creates a highly complicated volcanic complex due to the admixture of lava flow and sedimentation of volcanic ash. The Song-Gok site is composed of volcanic rocks that collapsed at the lower part of the slope, in combination with several discontinuities in and around a fault. The results of projection analysis indicated the possibility of plane, wedge, and toppling failure in the failure section. The results of discontinuity modeling using the Distinct Element Method (DEM) revealed a total displacement of 207 mm and a joint shear displacement of 114 mm. The yield surface zone was verified at the fault plane of the failure section. In geotechnical terms, volcanic rock slopes are characteristically vulnerable to failure because of differential weathering among the various rock types, the effect of groundwater based on the permeability of the rocks, and the presence of systematic joints generated by the cooling and contraction of lava. When considering the stability of a volcanic rock slope, it is necessary to consider data such as the geological features of the rock, as obtained through detailed geological survey, and variations in discontinuities and rock blocks.

• 한국석유지질학회:학술대회논문집
• /
• autumn
• /
• pp.28-46
• /
• 2001

Natural gas is a mixture of hydrocarbon gases and impurities such as nitrogen, hydrogen sulfide, and carbon dioxide and a clean energy producing no pollution materials for combustion. Currently, the demand of the natural gas is rapidly increasing due to worldwide environmental problems. According to Hubbert's study in the past, the natural gas was predicted as rapidly depleted resources, and then the results led to high gas price and limitation of usage during 1980s. Afterward, the study of natural gas resources based on geology identified the additional natural gas resources that were not considered in Hubbert's study. They are unconventional gas, additional resources in the existed reservoirs, and natural gas in deep subsurface areas. Such additional resouces made the future of natural gas bright and pormised low and stable gas price in the future. Deep natural gas is defined as the gas existing at or below 15,000ft$(4,752{\cal}m)$ in depth from the surface. According to the study from the U.S. Geological Survey(USGS) in 1995, 1,412 TCF of technically recoverable natural gas was remained to be discovered or developed in the onshore of United States. A significant part of that resource base, 114 TCF, exists at deep sedimentary basins, and it shows wide distribution with various geological environments. In 1995, the deep gas contributed to $6.7\% of total supply amount of natural gas in the United States and is expected to be $18.7\% by 201.5. However, the development of the deep gas is a high risky business due to expensive investment and high portion of dry holes, although it is developed. Thus, for developing the deep gas economically, it is necessary to overcome many technical challenges. In this paper, for increasing success rate of the deep gas, 1) geologic and compositional characteristics, and production cost have been analyzed according to depth, 2) technical problems related to deep gas production have been summarized, and 3) finally future study areas for increasing application of the deep gas have been suggested. For reference, this paper was written based on the study results from USGS and Gas Research Institute(GRI), for the United States is doing the most active R&D in the deep gas area, and thus, has many reliable data.