• Geophysics and Geophysical Exploration
• /
• v.9 no.3
• /
• pp.207-224
• /
• 2006

It has been widely known that the seismic piezo-cone penetration test (SCPTu) is one of the most useful techniques for investigating the geotechnical characteristics such as static and dynamic soil properties. As practical applications in Korea, SCPTu was carried out at two sites in Busan and four sites in Incheon, which are mainly composed of alluvial or marine soil deposits. From the SCPTu waveform data obtained from the testing sites, the first arrival times of shear waves and the corresponding time differences with depth were determined using the cross-over method, and the shear wave velocity $(V_S)$ profiles with depth were derived based on the refracted ray path method based on Snell's law. Comparing the determined $V_S$ profile with the cone tip resistance $(q_t)$ profile, both trends of profiles with depth were similar. For the application of the conventional CPTu to earthquake engineering practices, the correlations between $V_S$ and CPTu data were deduced based on the SCPTu results. For the empirical evaluation of $V_S$ for all soils together with clays and sands which are classified unambiguously in this study by the soil behavior type classification index $(I_C)$, the authors suggested the $V_S-CPTu$ data correlations expressed as a function of four parameters, $q_t,\;f_s,\;\sigma'_{v0}$ and $B_q$, determined by multiple statistical regression modeling. Despite the incompatible strain levels of the downhole seismic test during SCPTu and the conventional CPTu, it is shown that the $V_S-CPTu$ data correlations for all soils, clays and sands suggested in this study is applicable to the preliminary estimation of $V_S$ for the soil deposits at a part in Korea and is more reliable than the previous correlations proposed by other researchers.

• 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.10 no.3
• /
• pp.191-202
• /
• 2007

To identify the detailed attenuation structure in the southern Korean Peninsula, a numerical test was conducted for the Q tomography inversion to be applied to the accumulated dataset until 2005. In particular, the stochastic pointsource ground-motion model (STGM model; Boore, 2003) was adopted for the 2D Q tomography inversion for direct application to simulating the strong ground-motion. Simultaneous inversion of the STGM model parameters with a regional single Q model was performed to evaluate the source and site effects which were necessary to generate an artificial dataset for the numerical test. The artificial dataset consists of simulated Fourier spectra that resemble the real data in the magnitude-distance-frequency-error distribution except replacement of the regional single Q model with a checkerboard type of high and low values of laterally varying Q models. The total number of Q blocks used for the checkerboard test was 75 (grid size of $35{\times}44km^2$ for Q blocks); Q functional form of $Q_0f^{\eta}$ ($Q_0$=100 or 500, 0.0 < ${\eta}$ < 1.0) was assigned to each Q block for the checkerboard test. The checkerboard test has been implemented in three steps. At the first step, the initial values of Q-values for 75 blocks were estimated. At the second step, the site amplification function was estimated by using the initial guess of A(f) which is the mean site amplification functions (Yun and Suh, 2007) for the site class. The last step is to invert the tomographic Q-values of 75 blocks based on the results of the first and second steps. As a result of the checkerboard test, it was demonstrated that Q-values could be robustly estimated by using the 2D Q tomography inversion method even in the presence of perturbed source and site effects from the true input model.

• 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.15 no.1
• /
• pp.39-46
• /
• 2012

Recently, stratigraphic reservoirs are getting more attention than structural reservoirs which have mostly developed. However, recognizing stratigraphic thin gas reservoirs in a stacked section is usually difficult because of tuning effects. Moreover, if the reflections from the brine-saturated region of a thin layer have the same polarity with those from the gas-saturated region, we could not easily identify the gas reservoir with conventional data processing technique. In this study, we introduced a way to delineate the gas-saturated region in a thin layer reservoir using a spectral decomposition method. First of all, amplitude spectrum with the variation of the frequency and the incident angle was investigated for the medium which represents property of Class 3, Class 1 or Class 4 AVO response. The results show that the maximum difference in the amplitude spectra between brine and gas-saturated thin layers occurs around the peak frequency independent of the incident angle and the type of AVO responses. In addition, the amplitude spectra of the gas-saturated zone are greater than those of brine-saturated one in Class 3 and Class 4 at the peak frequency while those of phenomenon occur oppositely in Class 1. Based on the results, we applied spectral decomposition method to the stacked section in order to distinguish the gas-saturated zone from the brine-saturated zone in a thin layer reservoir. To verify our new method, we constructed a thin-layer velocity model which contains both gas and brine-saturated zones which have the same reflection polarities. As a result, in the spectral decomposed sections near the peak frequency obtained by Wigner-Ville Distribution (WVD), we could identify the difference between reflections from gas- and brinesaturated region in the thin layer reservoir, which was hardly distinguishable in the stacked section.

• Geophysics and Geophysical Exploration
• /
• v.1 no.1
• /
• pp.25-30
• /
• 1998

Time-lapse crosswell seismic data, recorded before and after the cavity filling, showed that the filling increased the velocity at a known cavity zone in an old mine site in Inchon area. The seismic response depicted on the tomogram and in conjunction with the geologic data from drillings imply that the size of the cavity may be either small or filled by debris. In this study, I attempted to evaluate the filling effect by analyzing velocity measured from the time-lapse tomograms. The data acquired by a downhole airgun and 24-channel hydrophone system revealed that there exists measurable amounts of source statics. I presented a methodology to estimate the source statics. The procedure for this method is: 1) examine the source firing-time for each source, and remove the effect of irregular firing time, and 2) estimate the residual statics caused by inaccurate source positioning. This proposed multi-step inversion may reduce high frequency numerical noise and enhance the resolution at the zone of interest. The multi-step inversion with different starting models successfully shows the subtle velocity changes at the small cavity zone. The inversion procedure is: 1) conduct an inversion using regular sized cells, and generate an image of gross velocity structure by applying a 2-D median filter on the resulting tomogram, and 2) construct the starting velocity model by modifying the final velocity model from the first phase. The model was modified so that the zone of interest consists of small-sized grids. The final velocity model developed from the baseline survey was as a starting velocity model on the monitor inversion. Since we expected a velocity change only in the cavity zone, in the monitor inversion, we can significantly reduce the number of model parameters by fixing the model out-side the cavity zone equal to the baseline model.

• Geophysics and Geophysical Exploration
• /
• v.23 no.1
• /
• pp.22-37
• /
• 2020

We analyze the comprehensive three-dimensional (3D) magnetic structure characteristics from the seafloor to the deep layer of the Tofua Arc (TA) 12 seamount in the Tonga Arc, Southwestern Pacific, using bathymetric and geomagnetic data, and magnetization vector inversion (MVI) results. The seafloor features surrounding TA 12 highlight a NW-SE-oriented elliptical caldera at the summit of the seamount, two small cones in the depressed area. A large-scale sea valley is present on the western flank of the seamount, extending from these caldera cones to the southwestern base of the seamount. TA 12 seamount exhibits a low magnetic anomaly in the caldera depression, whereas a high magnetic anomaly is observed surrounding the low magnetic anomaly across the caldera summit. It is therefore presumed that there may be a strong magnetic material distribution or magma intrusion in the caldera. The 3D MVI results show that the high anomaly zones are mainly present in the surrounding slopes of the seamount from the seafloor to the -3,000 m (below the seafloor) level, with the magnetic susceptibility intensity increasing as the seafloor level increases at the caldera depression. However, small high anomaly zones are present across the study area near the seafloor level. Therefore, we expect that the magma ascent in TA 12 seamount migrated from the flanks to the depression area. Furthermore, we assume that the complex magnetic distribution near the seafloor is due to the remnant magnetization.