• Geophysics and Geophysical Exploration
• /
• v.9 no.4
• /
• pp.271-278
• /
• 2006

Investigations of underground cavities are required to provide useful information for the reinforcement design and monitoring of the ground subsidence areas. It is, therefore, necessary to develop integrated geophysical techniques incorporating different geophysical methods in order to accurately image and to map underground cavities in the ground subsidence areas. In this study, we conducted geophysical investigations for development of integrated geophysical techniques to detect underground cavities at the field test site in the ground subsidence area, located at Yongweol-ri, Muan-eup, Muan-gun, Jeollanam-do. We examined the applicability of geophysical methods such as electrical resistivity, electromagnetic, and microgravity to cavity detection with the aid of borehole survey results. The underground cavities are widely present within the limestone bedrock overlain by the alluvial deposits in the test site where the ground subsidences have occurred in the past. The limestone cavities are mostly filled with groundwater or clays saturated with water in the site. The cavities, thus, have low electrical resistivity and density compared to the surrounding host bedrock. The results of the study have shown that the zones of low resistivity and density correspond to the zones of the cavities identified in the boreholes at the site, and that the geophysical methods used are very effective to detect the underground cavities. Furthermore, we could map the distribution of cavities more precisely with the study results incorporated from the various geophysical methods. It is also important to notice that the microgravity method, which has rarely used in Korea, is a very promising tool to detect underground cavities.

• Geophysics and Geophysical Exploration
• /
• v.10 no.3
• /
• pp.173-182
• /
• 2007

It is well known that manganese nodules enriched with valuable metals are abundantly distributed in the abyssal plain area in the Clarion-Clipperton (C-C) fracture zone of the northeast Pacific. Previous studies using deep-sea camera (DSC) system reported different observations about the relation of seafloor topographic change and nodule abundance, and they were sometimes contradictory. Moreover, proper foundation on the estimation of DSC underwater position, was not introduced clearly. The variability of the mining condition of manganese nodule according to seafloor topography was examined in the Korea Deep Ocean Study (KODOS) area, located in the C-C zone. In this paper, it is suggested that the utilization of deep towing system such as DSC is very useful approach to whom are interested in analysing the distributional characteristics of manganese nodule filed and in selecting promising minable area. To this purpose, nodule abundance and detailed bathymetry were acquired using deep-sea camera system and multi-beam echo sounder, respectively on the seamount free abyssal hill area of southern part ($132^{\circ}10'W$, $9^{\circ}45'N$) in KODOS regime. Some reasonable assumptions were introduced to enhance the accuracy of estimated DSC sampling position. The accuracy in the result of estimated underwater position was verified indirectly through the comparison of measured abundances on the crossing point of neighboring DSC tracks. From the recorded seafloor images, not only nodules and sediments but cracks and cliffs could be also found frequently. The positions of these probable unminable area were calculated by use of the recorded time being encountered with them from the seafloor images of DSC. The results suggest that the unminable areas are mostly distributed on the slope sides and hill tops, where nodule collector can not travel over.

• Geophysics and Geophysical Exploration
• /
• v.10 no.3
• /
• pp.183-190
• /
• 2007

The site effects of seismic stations were evaluated by conducting a simultaneous inversion of the stochastic point-source ground-motion model (STGM model; Boore, 2003) parameters based on the accumulated dataset of horizontal shear-wave Fourier spectra. A model parameter $K_0$ and frequency-dependent site amplification function A(f) were used to express the site effects. Once after a H/V ratio of the Fourier spectra was used as an initial estimate of A(f) for the inversion, the final A(f) which is considered to be the result of combined effect of the crustal amplification and loca lsite effects was calculated by averaging the log residuals at the site from the inversion and adding the mean log residual to the H/V ratio. The seismic stations were classified into five classes according to $logA_{1-10}^{max}$(f), the maximum level of the site amplification function in the range of 1 Hz < f < 10 Hz, i.e., A: $logA_{1-10}^{max}$(f) < 0.2, B: 0.2 $\leq$ $logA_{1-10}^{max}$(f) < 0.4, C: 0.4 $\leq$ $logA_{1-10}^{max}$(f) < 0.6, D: 0.6 $\leq$ $logA_{1-10}^{max}$(f) < 0.8, E: 0.8 $\leq$ $logA_{1-10}^{max}$(f). Implication of the classified result was supported by observing a shift of the dominant frequency of average A(f) for each classified stations as the class changes. Change of site classes after moving seismic stations to a better site condition was successfully described by the result of the station classification. In addition, the observed PGA (Peak Ground Acceleration)-values for two recent moderate earthquakes were well classified according to the proposed station classes.

• Geophysics and Geophysical Exploration
• /
• v.6 no.3
• /
• pp.126-133
• /
• 2003

High-resolution shallow marine seismic surveys have been carried out for the resources exploration, engineering applications and Quaternary mapping. To improve the resolution of subsurface structure image, multichannel digital technique has been applied. The quality of the image depends on the vertical and horizontal resolution and signal to noise (S/N) ratio which are associated with the data acquisition parameters such as sample interval, common midpoint (CMP) interval and CMP fold. To understand the effect of the acquisition parameters, a test survey was carried out off Yeosu and the acquired data were analyzed. A 30 $in^3$ small air gun was used as a seismic source and 8 channel streamer cable with a 5 m group interval was used as a receiver. The data were digitally recorded with a shot interval of 2 s and sample interval of 0.1 ms. The acquired data were resampled with various sample intervals, CMP intervals and CMP folds. The resampled data were processed, plotted as seismic sections and compared each other. The analysis results show that thin bed structure with ${\~}1m$ thickness and ${\~}6^{\circ}$ slope can be imaged with good resolution and continuity and low noise using the acquisition parameters with a sample interval shorter than 0.2 ms, CMP interval shorter than 2.5 m and CMP fold more than 4. Because seismic resolution is associated with the acquisition parameters, the quality of the subsurface structure can be imaged successfully using suitable and optimum acquisition parameters.

• Geophysics and Geophysical Exploration
• /
• v.13 no.4
• /
• pp.357-363
• /
• 2010

In this study, we investigated the geologic structure of the basement and overlying sediments of the construction site of the dinosaur egg fossil museum in Hwasung, Gyeonggi Province through refraction seismology, drilling, and downward seismic velocity measurements in the drill holes. The construction site ($350{\times}750\;m^2$) is located in the reclaimed area south of Sihwa Lake, Gojeong-ri. About 6,950 m of seismic refraction data consisting of 11 lines were acquired using a sledge hammer source. Drilling to the basement was performed at five sites. Sediment samples from drilling were analysed for grain-size distribution and age dating. At two drill holes, seismic velocity was measured with depth using a hammer as a seismic source. The geological structure of the study area consists of, from top to bottom, a tidal flat layer (5 ~ 12 m thick), a weathered soil layer (2 ~ 8 m thick), and the basement. The basement is interpreted as Cretaceous sedimentary rocks that tend to be shallow eastward. The volume of the tidal flat sediments and weathered soil in the study area is estimated as $1.4{\times}10^6\;m^3$, weighing $3.5{\times}10^6$ tons. The rate of sea level rise since 8,000 yrs BP is estimated to be 0.1 ~ 0.15 cm/yr.

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.127-136
• /
• 2008

To identify the potential area of gas hydrate in the Ulleung Basin, 2-D and 3-D seismic surveys using R/V Tamhae II were conducted in 2005 and 2006. Seismic survey equipment consisted of navigation system, recording system, streamer cable and air-gun source. For reliable velocity analysis in a deep sea area where water depths are mostly greater than 1,000 m and the target depth is up to about 500 msec interval below the seafloor, 3-km-long streamer and 1,035 $in^3$ tuned air-gun array were used. During the survey, a suite of quality control operations including source signature analysis, 2-D brute stack, RMS noise analysis and FK analysis were performed. The source signature was calculated to verify its conformity to quality specification and the gun dropout test was carried out to examine signature changes due to a single air gun's failure. From the online quality analysis, we could conclude that the overall data quality was very good even though some seismic data were affected by swell noise, parity error, spike noise and current rip noise. Especially, by checking the result of data quality enhancement using FK filtering and missing trace restoration technique for the 3-D seismic data inevitably contaminated with current rip noises, the acquired data were accepted and the field survey could be conducted continuously. Even in survey areas where the acquired data would be unsuitable for quality specification, the marine seismic survey efficiency could be improved by showing the possibility of noise suppression through onboard data processing.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.383-392
• /
• 2007

Karstic features and mining-related cavities not only lead to severe restrictions in land utilizations, but also constitute serious concern about geohazard and groundwater contamination. A microgravity survey was applied for detecting, mapping and monitoring karstic cavities in the test site at Muan prepared by KIGAM. The gravity data were collected using an AutoGrav CG-3 gravimeter at about 800 stations by 5 m interval along paddy paths. The density distribution beneath the profiles was drawn by two dimensional inversion based on the minimum support stabilizing functional, which generated better focused images of density discontinuities. We also imaged three dimensional density distribution by growing body inversion with solution from Euler deconvolution as a priori information. The density image showed that the cavities were dissolved, enlarged and connected into a cavity network system, which was supported by drill hole logs. A time-lapse microgravity was executed on the road in the test site for monitoring the change of the subsurface density distribution before and after grouting. The data were adjusted for reducing the effects due to the different condition of each survey, and inverted to density distributions. They show the change of density structure during the lapsed time, which implies the effects of grouting. This case history at the Muan test site showed that the microgravity with accuracy and precision of ${\mu}Gal$ is an effective and practical tool for detecting, mapping and monitoring the subsurface cavities.

• Geophysics and Geophysical Exploration
• /
• v.21 no.3
• /
• pp.171-182
• /
• 2018

We performed the seismic field data processing using an open-source software (Madagascar) to verify if it is applicable to processing of field data, which has low signal-to-noise ratio and high uncertainties in velocities. The Madagascar, based on Python, is usually supposed to be better in the development of processing technologies due to its capabilities of multidimensional data analysis and reproducibility. However, this open-source software has not been widely used so far for field data processing because of complicated interfaces and data structure system. To verify the effectiveness of the Madagascar software on field data, we applied it to a typical seismic data processing flow including data loading, geometry build-up, F-K filter, predictive deconvolution, velocity analysis, normal moveout correction, stack, and migration. The field data for the test were acquired in Gunsan Basin, Yellow Sea using a streamer consisting of 480 channels and 4 arrays of air-guns. The results at all processing step are compared with those processed with Landmark's ProMAX (SeisSpace R5000) which is a commercial processing software. Madagascar shows relatively high efficiencies in data IO and management as well as reproducibility. Additionally, it shows quick and exact calculations in some automated procedures such as stacking velocity analysis. There were no remarkable differences in the results after applying the signal enhancement flows of both software. For the deeper part of the substructure image, however, the commercial software shows better results than the open-source software. This is simply because the commercial software has various flows for de-multiple and provides interactive processing environments for delicate processing works compared to Madagascar. Considering that many researchers around the world are developing various data processing algorithms for Madagascar, we can expect that the open-source software such as Madagascar can be widely used for commercial-level processing with the strength of expandability, cost effectiveness and reproducibility.

• Geophysics and Geophysical Exploration
• /
• v.21 no.3
• /
• pp.150-161
• /
• 2018

Recently, unmanned aircraft EM (electromagnetic) survey based on ICT (Information and Communication Technology) has been widely utilized because of the efficiency in regional survey. We performed the theoretical study on the unmanned airship EM system developed by KIGAM (Korea Institute of Geoscience and Mineral resources) as part of the practical application of unmanned aircraft EM survey. Since this system has different configurations of transmitting and receiving loops compared to the conventional aircraft EM systems, a new technique is required for the appropriate interpretation of measured responses. Therefore, we proposed a method to calculate the EM field for the arbitrary shaped transmitter and verified its validity through the comparison with analytic solution for circular loop. In addition, to simulate the magnetic responses by three-dimensionally (3D) distributed anomalies, we have adapted our algorithm to 3D frequency-domain EM modeling algorithm based on the edge-FEM (finite element method). Though the analysis on magnetic field responses from a subsurface anomaly, it was found that the response decreases as the depth of the anomaly increases or the flight altitude increases. Also, it was confirmed that the response became smaller as the resistivity of the anomaly increases. However, a nonlinear trend of the out-of-phase component is shown depending on the depth of the anomaly and the used frequency, that makes it difficult to apply simple analysis based on the mapping of the magnitude of the responses and can cause the non-uniqueness problem in calculating the apparent resistivity. Thus, it is a prerequisite to analyze the appropriate frequency band and flight altitude considering the purpose of the survey and the site conditions when conducting a survey using the unmanned aircraft EM system.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.280-290
• /
• 2002

Three-dimensional (3-D) electromagnetic (EM) modeling algorithm has been developed using finite element method (FEM) to acquire more efficient interpretation techniques of EM data. When FEM based on nodal elements is applied to EM problem, spurious solutions, so called 'vector parasite', are occurred due to the discontinuity of normal electric fields and may lead the completely erroneous results. Among the methods curing the spurious problem, this study adopts vector element of which basis function has the amplitude and direction. To reduce computational cost and required core memory, complex bi-conjugate gradient (CBCG) method is applied to solving complex symmetric matrix of FEM and point Jacobi method is used to accelerate convergence rate. To verify the developed 3-D EM modeling algorithm, its electric and magnetic field for a layered-earth model are compared with those of layered-earth solution. As we expected, the vector based FEM developed in this study does not cause ny vector parasite problem, while conventional nodal based FEM causes lots of errors due to the discontinuity of field variables. For testing the applicability to high frequencies 100 MHz is used as an operating frequency for the layer structure. Modeled fields calculated from developed code are also well matched with the layered-earth ones for a model with dielectric anomaly as well as conductive anomaly. In a vertical electric dipole source case, however, the discontinuity of field variables causes the conventional nodal based FEM to include a lot of errors due to the vector parasite. Even for the case, the vector based FEM gave almost the same results as the layered-earth solution. The magnetic fields induced by a dielectric anomaly at high frequencies show unique behaviors different from those by a conductive anomaly. Since our 3-D EM modeling code can reflect the effect from a dielectric anomaly as well as a conductive anomaly, it may be a groundwork not only to apply high frequency EM method to the field survey but also to analyze the fold data obtained by high frequency EM method.