• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.85-92
• /
• 2008

Magnetotelluric (MT) methods are widely applied as an effective exploration technique to geothermal surveys. Two-dimensional (2-D) analysis is frequently used to investigate a complicated subsurface structure in a geothermal region. A 2-D finite-element method (FEM) is usually applied to the MT analysis, but we must pay attention to the accuracy of so-called auxiliary fields. Rodi (1976) proposed an algorithm of improving the accuracy of auxiliary fields, and named it as the MOM method. Because it introduces zeros into the diagonal elements of coefficient matrix of the FEM total equation, a pivoting procedure applied to the symmetrical band matrix makes the numerical solution far less efficient. The MOM method was devised mainly for the inversion analysis, in which partial derivatives of both electric and magnetic fields with respect to model parameters are required. In the case of forward modeling, however, we do not have to resort to the MOM method; there is no need of modifying the coefficient matrix, and the auxiliary fields can be elicited from the regular FEM solution. The computational efficiency of the MOM method, however, can be greatly improved through a sophisticated rearrangement of the total equation.

• The Journal of Engineering Geology
• /
• v.14 no.1
• /
• pp.93-104
• /
• 2004

This paper intended to assess the hydro-structure characteristics of volcanic rocks based on the hydrogeological data obtained from the underground storage cavern during construction. The variation of groundwater levels was periodically measured from the 28 surface monitoring holes(NX size) and the hydraulic pressures and injection rates were daily monitored from the water curtain holes(95 horizontal holes and 63 vertical holes). The hydraulic interference tests were performed in whole water curtain holes. The distribution patterns of hydraulic pressure are closely related to the dip angles of fracture intersected to the water curtain holes. Three domains can be grouped by the distribution of hydraulic pressures in the horizontal water curtain holes. The initial hydraulic pressures measured immediately after drilling of water crutain holes are high in ascending order of the cavern C-2, C-1, and C-3. The priliminary hydrochemical data also indicate that the portions of the deep groundwater composition is relatively great in the cavern C-3 area. Some of the horizontal water curtain holes in the cavern C-3 show a steady higher groundwater pressure with the composition of shallow groundwater indicating the outer boundary as constant hydraulic boundary. The water curtain holes in the cavern C-2 is characterized as low initial hydraulic pressure and less injection rates, suggesting poor hydraulic connectivity to a shallow groundwater system. The results of the study can help to understand a hydraulic compartment concept in a fracture hydro-geology and be utilized during the surface investigation for a groundwater system.

• Economic and Environmental Geology
• /
• v.43 no.6
• /
• pp.615-623
• /
• 2010

Seismic refraction survey is a geophysical method that delineates subsurface velocity structure using direct wave and critically refracted wave. The generalized reciprocal method(GRM) is an inversion technique which uses travel-time data from several forward and reverse shots and which can provide the geometry of irregular inclined refractors and structures underlain by hidden layer such as low velocity zone and thin layer. In this study, a simple Excel-GRM routine was tested for fast mapping of the interface between weathering layer and bedrock during the survey, with employing a pair of forward and reverse shots. This routine was proved to control the maximum dip of approximately $30^{\circ}C$ and maximum velocity contrast of 0.6, based on the panel tests in terms of dipping angle and velocity contrast for the two-layer inclined models. In contrast with conventional operation of five to seven shots with sufficient offset distance and indoor data analysis thereafter, this routine was performed in the field shortly after data acquisition. Depth to the bedrock provided by Excel-GRM, during the field survey for Cheongju granite area, correlates well with the elevation of the surface of soft rock from the drill core and SPS logging data. This cost-effective routine developed for quickly delineating the bedrock surface in the field survey will be readily applicable to mapping of weathering zone in narrow zone with small variation of elevation of bedrock.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.25 no.2
• /
• pp.141-155
• /
• 2023

Monitoring technology based on electrical resistivity is widely used for non-destructive data collection and health analysis of underground structures and tunnels. Vulnerable sections such as fault zone generates many problems during construction of the tunnel. These problems cause displacement and stress changes of the ground. Therefore, it is necessary to predict the state of the fault zone section to ensure the mechanical stability of the underground structure. Monitoring the size of joints and the porosity of the fillings is essential for rocks. Previous studies have not considered the variety of fillings in rock joints. In this study, electrical resistivity tests were conducted according to the particle mixing state of the sandy fillings. When the size of fillings is decreased at the constant porosity, the electrical resistivity tends to increase. The results of this study are expected to be useful as basic electrical resistivity data for predicting the ground conditions and evaluation of the ground behavior that is containing sandy fillings in the rock joint for tunnels.

• Journal of the Korean Geotechnical Society
• /
• v.23 no.2
• /
• pp.47-60
• /
• 2007

Conventional active surface wave measurements performed using a transient or continuous source are often limited in the maximum depth of penetration due to the difficulty of generating low-frequency energy with reasonably portable sources. This limitation may inhibit accurate seismic site response calculations because of the inability to define deeper subsurface structure. By measuring surface wave generated by passive sources including microtremors and cultural noise, it is possible to overcome this problem and develop soil stiffness profiles to much larger depth. Reliability of dispersion estimates from the passive surface wave measurements is critical to present reliable shear wave velocity profiles and can be improved by the measurements and analyses of passive surface waves based on correct understanding of systematic errors included in passive dispersion data. In this study, the systematic errors caused by poor wavenumber resolution and energy leakage into sidelobes in passive tests are mainly explored. Recommendations for reliable passive surface wave measurements and dispersion estimates are presented and illustrated at a site in San Jose, California, U.S.

• The Journal of Engineering Geology
• /
• v.32 no.1
• /
• pp.157-171
• /
• 2022

Electrical resistivity tomography (ERT) is a traditional and representative geophysical method for determining the resistivity distributions of surrounding soil and rock volumes. Depth resolution profiles and sensitivity distribution sections of the resistivities with respect to various electrode configurations are calculated and investigated using numerical model data. Shallow vertical resolution decreases in the order of Wenner, Schlumberger, and dipole-dipole arrays. A high investigable depth in homogeneous medium is calculated to be 0.11-0.19 times the active electrode spacing, but is counterbalanced by a low vertical resolution. For the application of ERT depth resolution profiles and sensitivity distributions, we provide subsurface structure models for two types of land-creping failure (planar and curved), subvertical fracture, and weathered layer over felsic and mafic igneous rocks. The dipole-dipole configuration appears to be most effective for mapping land-creeping failure planes (especially for curved planes), whereas the Wenner array gives the best resolution of soil horizons and shallow structures in the weathered zone.

• Journal of Korean Society of Disaster and Security
• /
• v.16 no.4
• /
• pp.157-164
• /
• 2023

The estuary dam is a structure installed and operated in a closed state except when flood event occurs to prevent inland saltwater intrusion and secure freshwater supply. However, the closed state of dam leads to issues such as eutrophication, so it is necessary to examine the extent of saltwater intrusion resulting from the opening of sluice gates. Groundwater, due to its subsurface conditions and slow flow characteristics, is widely analyzed using numerical models. OpenGeoSys, an open-source software capable of simulating Thermal- Hydraulic- Mechanical- Chemical phenomena, was adopted for this study. Simulations were conducted assuming natural flow conditions without dam and operating considering busy farming season, mostly from March to September. Verification of the model through analytical solutions showed error of 3.7%, confirming that OpenGeoSys is capable of simulating saltwater intrusion for these cases. From results simulated for 10 years, considering for the busy farming season, resulted in about 46% reduction in saltwater intrusion length compared to natural flow conditions, approximately 74.36 m. It may be helpful to make choices to use groundwater as a water resource.

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.80-87
• /
• 2011

Autonomous underwater vehicles (AUVs) present the important advantage of being able to approach the seafloor more closely than surface vessel surveys can. To collect bathymetric data, bottom material information, and sub-surface images, multibeam echosounder, sidescan sonar (SSS) and subbottom profiler (SBP) equipment mounted on an AUV are powerful tools. The 3000m class AUV URASHIMA was developed by the Japan Agency for Marine-Earth Science and Technology (JAMSTEC). After finishing the engineering development and examination phase of a fuel-cell system used for the vehicle's power supply system, a renovated lithium-ion battery power system was installed in URASHIMA. The AUV was redeployed from its prior engineering tasks to scientific use. Various scientific instruments were loaded on the vehicle, and experimental dives for science-oriented missions conducted from 2006. During the experimental cruise of 2007, high-resolution acoustic images were obtained by SSS and SBP on the URASHIMA around the northern Kumano Basin off Japan's Kii Peninsula. The map of backscatter intensity data revealed many debris objects, and SBP images revealed the subsurface structure around the north-eastern end of our study area. These features suggest a structure related to the formation of the latest submarine fan. However, a strong reflection layer exists below ~20 ms below the seafloor in the south-western area, which we interpret as a denudation feature, now covered with younger surface sediments. We continue to improve the vehicle's performance, and expect that many fruitful results will be obtained using URASHIMA.

• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.31 no.2
• /
• pp.259-266
• /
• 1998

A study on the biological and chemical characteristics in the middle last Sea of Korea was carried out at 31 stations in October $11\~18$, 1995 on board the R/V Tam-Yang. The chlorophyll a concentration, new and regenerated production, and the vertical diffusion of nitrate from the thermocline structure were investigated. From the vertical distribution of chlorophyll a, subsurface maxima were observed near the thermorline at most stations including the frontal zone, except at the southern stations where the maximum chloropyll a concentration occurred at the surface, The nanophytoplankton was the most dominant fraction comprising $83.5\%$ of total phytoplankton cell numbers, but netphytoplankton were common at the southern stations where the dominant species were Rhizosolenia sp. Nitrogenous new production and regenerated productions were measured using the stable isotope $^{15}N$ nitrate and ammonia uptake method. The vertically integrated nitrogen production varied between 8.470 and $72.945\;mg\;N\;m^{-2}\;d^{-1}$. The f-ratio, which is the traction of new production from primary production, waried between 0.03 and 0.72, indicating that $3\%$ to $72\%$ of primary production was supported by the input of nutrients from below the euphotic zone and the rest are supported by ammonia recycled within the euphotic layer. This range of f-ratio encompasses from extremely oligotrophic to eutrophic area characteristics. The differences in productivity and f-ratio among stations were related to frontal structure and the bottom topography. The values were high near the frontal zone and low outside of it, and the station near Ulleng Island showed the highest f-ratio. Vertical diffusion coefficients were calculated from both the water column stability (Kz-1) of King and Devol's equation (1979) and new nitrogen requirement (Kz-2). The values of Kz-2 ($0.11\~0.55\;cm^2/s$) were relatively low compared to the values reported previously.

• Economic and Environmental Geology
• /
• v.12 no.2
• /
• pp.95-104
• /
• 1979

The position of any point on the earth's surface can be. represented in the spherical coordinates by surface spherical harmonics. Since geomagnetic field is a function of position on the earth, it can be also expressed by spherical harmonic analysis as spherical harmonics of trigonometric series of $a_m({\theta})$ cos $m{\phi}$ and $b_m({\theta})$ sin $m{\phi}$. Coefficients of surface spherical harmonics, $a_m({\theta})$ and $b_m({\theta})$, can be drawn from the components of the geomagnetic field, declination and inclination, and vice versa. In this paper, components of geomagnetic field, declination and inclination in the Korean peninsula are obtained by spherical harmonic analysis using the Gauss coefficients calculated from the world-wide magnetic charts of 1960. These components correspond to the values of normal geomagnetic field having no disturbances of subsurface mass, structure, and so on. The vertical and total components offer the zero level for the interpretation of geomagnetic data obtained by magnetic measurement in the Korean peninsula. Using this zero level, magnetic anomaly map is obtained from the data of airborne magnetic. prospecting carried out during 1958 to 1960. The conclusions of this study are as follows; (1) The intensity of horizontal component of normal geomagnetic field in Korean peninsula ranges from $2{\times}10^4$ gammas to $2.45{\times}10^4$ gammas. It decreases about 500 with the increment of $1^{\circ}$ in latitude. Along the same. latitude, it increases 250 gammas with the increment of $1^{\circ}$ in longitude. (2) Intensity of vertical component ranges from $3.85{\times}10^4$ gammas to $5.15{\times}10^4$ gammas. It increases. about 1000 gammas with the increment of $1^{\circ}$ in latitude. Along the same latitude, it decreases. 150~240 gammas with the increment of $1^{\circ}$ in longitude. Decreasing rate is considerably larger in higher latitude than in lower latitude. (3) Total intensity ranges from $4.55{\times}10^4$ gammas to $5.15{\times}10^4$ gammas. It increases 600~700 gammas with the increament of $1^{\circ}$ in latitude. Along the same latitude, it decreases 10~90 gammas. with the increment of $1^{\circ}$ in longitude. Decreasing rate is considerably larger in higher latitude as the case of vertical component. (4) The declination ranges from $-3.8^{\circ}$ to $-11.5^{\circ}$. It increases $0.6^{\circ}$ with the increment of $1^{\circ}$ in latitude. Along the same latutude, it increases $0.6^{\circ}$ with the increment of l O in longitude. Unlike the cases of vertical and total component, the rate of change is considerably larger in lower latitude than in higher latitude. (5) The inclination ranges from $57.8^{\circ}$ to $66.8^{\circ}$. It increases about $1^{\circ}$ with 'the increment of $1^{\circ}$ in latitude Along the same latitude, it dereases $0.4^{\circ}$ with the increment of $1^{\circ}$ in longitude. (6) The Boundaries of 5 anomaly zones classified on the basis of the trend and shape of anomaly curves correspond to the geologic boundaries. (7) The trend of anomaly curves in each anomaly zone is closely related to the geologic structure developed in the corresponding zone. That is, it relates to the fault in the 3rd zone, the intrusion. of granite in the 1st and 5th zones, and mountains in the 2nd and 4th zones.