• Spatial Information Research
• /
• v.7 no.2
• /
• pp.159-173
• /
• 1999

GIS has been regarded as one of important tools or methodologies for various geoscience applications. Recently, spatial data integration schemes for site-specific or field-specific thematic mapping are newly developed and utilized. However, these kinds of approaches are somewhat insufficient quantitative assessment of integrated layers towards known targets in-detailed . Moreover, GIS analysis scheme is rarely extended to scientific approaches. In this study, simple approach of Multi-Buffer Zone Analysis , related to GIS analystical aspect, is addressed and an actual application for predicting or favorable mapping of mineral occurrences, one of GIS-based geoscientific approaches, is performed, As for geo-processing in GIS itself, this scheme can be regarded as extension or adaptiation of cell-based buffering or proximity analysis to geoscientific data interpretation. This study is based on rationale that surface geological pattern around primitives such as a point, a line, or a polygon in GIS, representing significant geological features, can be efficiently utilized to delineate complex geological behaviors or events, especially handling multiple dta sets originated from multiple sources such as airborne geophysical/radiometric exploration, field survey, and even a classified image of remote sensing. Conclusively, this methodology associated wit GIS is though to be helpful to analyze the spatial pattern of multiple data, pointing given sources, and is expected to effectively utilize for exploratory analysis of cell-based resultant layer integrated with complex or different data sources.

• Geophysics and Geophysical Exploration
• /
• v.13 no.4
• /
• pp.407-415
• /
• 2010

MT soundings were carried out in 2008, in northern Victoria, Australia, as a continuing collaboration research of 2007 between Republic of Korea, Australia, and Japan. The main purpose of this research is to investigate electrical conductivity structure and thus help understanding of tectonic structure in central Victoria, which is believed to be closely linked to mineralization and magmatic processes of this region. The survey area is located in western Lachlan Fold Belts, which is the part of Tasman Fold Belts in southeastern Australia. An MT profile of 2008 is almost parallel to the one of 2007 and approximately 50 km away. The 2D inversion result of MT data also shows that the position of conductivity discontinuity near surface are well matched with the positions of major faults, such as Avoca Fault, which is the structural boundary between Stawell and Bendigo Zones, and Heathcote Fault Zone, which marks the boundary between Bendigo and Melbourne Zones. It is also confirmed from resistivity image that internal faults in Bendigo Zone are in listric form, which is implied to be formed by structural shortening during compressional orogenic activity in Silurian.

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.109-115
• /
• 2008

It is very important to estimate the physical properties of survey area and delineate the geological basement in marine site survey for the design of offshore structures. For the purpose of providing high quality data by means of engineering site survey, it is necessary to apply several survey techniques and carry out the integrated interpretation to each other. In this study, we applied single channel seismic reflection method and OBC (Ocean Bottom Cable) type seismic refraction method at shallow marine. We used a dual boomer-single channel streamer as a source-receiver in seismic reflection survey and airgun source-the developed OBC type streamer in seismic refraction survey. We made 24 channels OBC type streamer which has 4m channel interval and each channel is composed of single hydrophone and preamplifier. We tested the field applicability of the proposed method and applied the typical seismic data processing methods to the obtained reflection data in order to enhance the data quality and image resolution. In order to estimate the geological velocity distribution from refraction data, seismic refraction tomography technique was applied. Therefore, we could successfully perform time-depth conversion using the velocity information as an integrated interpretation. The proposed method could provide reliable geologic information such as sediment layer thickness and 3D basement depth map.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.268-274
• /
• 2007

Traditional two-dimensional (2D) interpretation of magnetotelluric (MT) data utilizes only transverse magnetic (TM)-mode data, because 2D inversion of transverse electric (TE)-mode data results in spurious features when 3D structures exist in the subsurface. The application of a 3D inversion algorithm to a single MT profile can reduce contamination due to off-profile anomalies and help us to incorporate TE-mode data in the interpretation. In this study, we conduct 2D and 3D inversions of MT data observed along two lines in Jeju Island. First, we invert apparent resistivities and phases in the TM and TE modes separately. Then, we perform 2D joint inversion of both TM- and TE-mode data and 3D inversion of both Zxy- and Zyx-mode data corresponding to TE- and TM-mode data in 2D. The resistivity images derived from all four data show that the geoelectrical structure in Jeju Island is a three-layered earth with the resistive-conductive-resistive stratigraphy within a depth of 5 km. The 3D inversion does not produce clear anomalies in the reconstructed profile image, while all of 2D do. This attributed to the possibility that 2D inversion results are distorted by exiting off-profile 3D anomalies in Jeju. With 3D inversion of 2D profile MT data, we can deduce more reliable results that are not seriously distorted by off-profile 3D anomalies.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.250-256
• /
• 2002

Resistivity tomography has become an important tool to image underground resistivity distribution. This method has been widely applied to site investigation for engineering and environmental purpose. In resistivity tomography, various electrode arrays can be used and each array has both merits and demerits. For example, the pole-pole array has high signal to noise ratio (S/N ratio), but its resolution is too low. The dipole-dipole array has low S/N ratio, but its resolution is very high. The Pole-dipole may has intermediate Snf ratio and resolution. The modified Pole-dipole array, recently proposed, shows reasonable S/N ratio and resolution, which are comparable to the pole-dipole array. These electrode arrays except the pole-pole array, however, have the problem that the apparent resistivity can diverge at some special electrode Positions. Also, the Pole-Pole array may not reflect the doe resistivity of an anomalous body. In this study, we propose a new electrode array, mixed array, where pole-dipole and modified pole-dipole ways are selectively used with the relative positions of current and potential electrodes. The mixed array has the same level of S/N ratio and resolution as the pole-dipole array and the apparent resistivity does not diverge in the receiver hole. Furthermore, the apparent resistivity using the array can reflect the true resistivity of the anomalous body.

• Geophysics and Geophysical Exploration
• /
• v.18 no.2
• /
• pp.74-84
• /
• 2015

To yield physically meaningful images through elastic reverse-time migration, the wavefield separation which extracts P- and S-waves from reconstructed vector wavefields by using elastic wave equation is prerequisite. For expanding the application of the elastic reverse-time migration to anisotropic media, not only the anisotropic modelling algorithm but also the anisotropic wavefield separation is essential. The anisotropic wavefield separation which uses pseudo-derivative filters determined according to vertical velocities and anisotropic parameters of elastic media differs from the Helmholtz decomposition which is conventionally used for the isotropic wavefield separation. Since applying these pseudo-derivative filter consumes high computational costs, we have developed the efficient anisotropic wavefield separation algorithm which has capability of parallel computing by using GPUs (Graphic Processing Units). In addition, the highly efficient anisotropic elastic reverse-time migration algorithm using MPI (Message-Passing Interface) and incorporating the developed anisotropic wavefield separation algorithm with GPUs has been developed. To verify the efficiency and the validity of the developed anisotropic elastic reverse-time migration algorithm, a VTI elastic model based on Marmousi-II was built. A synthetic multicomponent seismic data set was created using this VTI elastic model. The computational speed of migration was dramatically enhanced by using GPUs and MPI and the accuracy of image was also improved because of the adoption of the anisotropic wavefield separation.

• Geophysics and Geophysical Exploration
• /
• v.14 no.3
• /
• pp.227-233
• /
• 2011

Attention has been focused on geothermal energy as an alternative energy because it is continuously operable without external supply. Most of geothermal anomalies in Korea are related to deep circulation of groundwater through a fracture system in granite area. Therefore it is very important to understand the distribution of the fracture system which is the main channel of ground water. In this research, we constructed the velocity models with a fracture system and the layered sediments, respectively, and generated synthetic data sets with them to verify the presented vertical seismic profiling (VSP) preprocessing scheme. We compared the results from conventional VSP preprocessing flow to those from VSP preprocessing flow considering fracture system. We noticed that the preprocessing flow considering fracture system retains more sufficient signal including down-going wave than conventional preprocessing. In addition, we applied 3D VSP prestack phase screen migration to the preprocessed reversed VSP (RVSP) data from Seokmo Island so that we were able to image fracture structure of the geothermal site in Seokmo Island.

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

• Economic and Environmental Geology
• /
• v.54 no.1
• /
• pp.91-103
• /
• 2021

In order to estimate the vertical and horizontal structural in the Yangsan fault core line (Naengsuri area, Pohang), we carried out gravity field measurements and interpretation procedures such as Euler deconvolution method and curvature analysis in addition to the forward modelling technique (i.e. IGMAS+). We found a prominent gravity difference of more than 1.5 mGal across the fault core. This indicates a distinct density difference between the western and eastern crustal area across the Yangsan fault line. Comparing this gravity field interpretation with other existent geologic and geophysical survey data (e.g. LiDAR, trenching, electric resistivity measurements), It is concluded that (1) the prominent gravity difference is caused by the density difference of about 0.1 g/㎤ between the Bulguksa Granite in the west and the Cretaceous Sandstone in the east side, (2) the fault core is elongated vertically into a depth of about 2,000 meters and extended horizontally 3,000 meters to the NNE direction from Naengsuri area. Our results present that the gravity field method is a very effective tool to estimate a three -dimensional image of the active fault core.

• Geophysics and Geophysical Exploration
• /
• v.25 no.3
• /
• pp.140-161
• /
• 2022

Seismic data with missing traces are often obtained regularly or irregularly due to environmental and economic constraints in their acquisition. Accordingly, seismic data interpolation is an essential step in seismic data processing. Recently, research activity on machine learning-based seismic data interpolation has been flourishing. In particular, convolutional neural network (CNN) and generative adversarial network (GAN), which are widely used algorithms for super-resolution problem solving in the image processing field, are also used for seismic data interpolation. In this study, CNN-based algorithm, U-Net and GAN-based algorithm, and conditional Wasserstein GAN (cWGAN) were used as seismic data interpolation methods. The results and performances of the methods were evaluated thoroughly to find an optimal interpolation method, which reconstructs with high accuracy missing seismic data. The work process for model training and performance evaluation was divided into two cases (i.e., Cases I and II). In Case I, we trained the model using only the regularly sampled data with 50% missing traces. We evaluated the model performance by applying the trained model to a total of six different test datasets, which consisted of a combination of regular, irregular, and sampling ratios. In Case II, six different models were generated using the training datasets sampled in the same way as the six test datasets. The models were applied to the same test datasets used in Case I to compare the results. We found that cWGAN showed better prediction performance than U-Net with higher PSNR and SSIM. However, cWGAN generated additional noise to the prediction results; thus, an ensemble technique was performed to remove the noise and improve the accuracy. The cWGAN ensemble model removed successfully the noise and showed improved PSNR and SSIM compared with existing individual models.