• Geophysics and Geophysical Exploration
• /
• v.4 no.3
• /
• pp.80-83
• /
• 2001

Four seismic reflection horizons in 3-D seismic data, coherence derived from the seismic data, and 38 well logs from the Boonsville Gas Filed in Texas were tried to be integrated and visualized in 3 dimensions. Time surface was constructed from pick times of the reflection horizons. Average velocities to each horizon at 38 well locations were calculated based on depth markers from the well logs and time picks from the 3-D seismic data. The time surface was transformed to depth surface through velocity interpolation. Coherence was calculated on the 3-D seismic data by semblance method. Spatial distribution of the coherence is captured easily in 3-D visualization. Comparing to a time-slice of seismic data, distinctive stratigraphic features could be correctly recognized on the 3-D visualization.

• The Journal of Engineering Geology
• /
• v.17 no.4
• /
• pp.615-621
• /
• 2007

Seismic refraction survey which is assumed to be effective for investigating of the amount of bed sediment was applied to two sites having different soil types. The ten lines of seismic surveys were carried out in the areas of the levee near the Wolchon Bridge, Yongin and the Damam levee, Yecheon. At the same time, two borehole data were obtained to be used in classifying geological formations on the inverted seismic sections. As the results of interpolating the depths of geological formations and the surface elevation values, it is possible to estimate the volume of the bed sediments at $24,000m^3\;and\;56,000m^3$ in Yongin and Yecheon, respectively.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2002.09a
• /
• pp.51-60
• /
• 2002

In this study we developed the visualization scheme of spatial ground-motion measurements in real time by using DSS data. Even though this scheme itself is useful for national earthquake mitigation plans, this scheme could be served as the crucial core for constructing rapid earthquake damage evaluation system. DSS is the abbreviation of Data Subscription Service and this is the pre-assigned request for the seismic stations to send very limited brief data with high priority and negligible transmission load. In addition to visualize the damage area with intensity, the corresponding epicenter can be estimated roughly for quick event alarm. For the interpolation of spatially irregular PGA data, the program, named as surface. of GMT was used with NetCDF grid file format. Since the grid file is similar to a postscript file, the program, called as shading, was coded with C language by using Matpak library in order to convert grid files into image files.

• Journal of the Korean Geographical Society
• /
• v.41 no.4 s.115
• /
• pp.391-403
• /
• 2006

The Hanon volcano located in the southern pan of Cheju Island, Korea has a wetland in its crater being used as a farmland. Previous researchers presumed this wetland was a maar lake in the past. Based on the seismic refraction method, the wetland sediment layer was estimated between 5 to 14 m deep, which is mostly in accordance with previous researches. However, this shows only the depths at some sites, not representing the whole spatial distribution. This study is an attempt to reconstruct the volcanic lake in Hanon crater by applying the spatial statistical techniques based on the depth information from the seismic survey and known data. The procedure of reconstruction is as follows: First, the depth information from the seismic survey and known data were collected and it was interpolated by IDW and Ordinary Kriging method. Next, with the interpolation map and the present DEM the paleo DEM was constructed. Finally, using the paleo lake level on core data, the boundary of volcanic lake was extracted from the paleo DEM. The reconstructed lake resembles a half-moon in the north of the central scoria cone. It is estimated that the lake was 5 m deep on average and 13 m deep at the deepest point. Although there are slight differences according to the interpolation techniques, it is calculated that the area of the lake was between 184,000 and $190000m^2,$ and its volume approximately $869,760m^3$. Because of the continuous deposition processes after the crater formation, the reconstructed volcanic lake would not indicate an actual lake at a specific time. Nevertheless, it offers a significant clue regarding the inner morphology and evolution of the crater.

• Economic and Environmental Geology
• /
• v.39 no.6 s.181
• /
• pp.711-717
• /
• 2006

In order to study gas hydrate, potential future energy resources, Korea Institute of Geoscience and Mineral Resources has conducted seismic reflection survey in the East Sea since 1997. one of evidence for presence of gas hydrate in seismic reflection data is a bottom simulating reflector (BSR). The BSR occurs at the interface between overlaying higher velocity, hydrate-bearing sediment and underlying lower velocity, free gas-bearing sediment. That is often characterized by large reflection coefficient and reflection polarity reverse to that of seafloor reflection. In order to apply depth migration to seismic reflection data. we need high performance computers and a parallelizing technique because of huge data volume and computation. Phase shift plus interpolation (PSPI) is a useful method for migration due to less computing time and computational efficiency. PSPI is intrinsically parallelizing characteristic in the frequency domain. We conducted conventional data processing for the gas hydrate data of the Ease Sea and then applied prestack depth migration using message-passing-interface PSPI (MPI_PSPI) that was parallelized by MPI local-area-multi-computer (MPI_LAM). Velocity model was made using the stack velocities after we had picked horizons on the stack image with in-house processing tool, Geobit. We could find the BSRs on the migrated stack section were about at SP 3555-4162 and two way travel time around 2,950 ms in time domain. In depth domain such BSRs appear at 6-17 km distance and 2.1 km depth from the seafloor. Since energy concentrated subsurface was well imaged we have to choose acquisition parameters suited for transmitting seismic energy to target area.

• Geophysics and Geophysical Exploration
• /
• v.9 no.3
• /
• pp.241-249
• /
• 2006

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms commonly used, however, may not give sufficiently precise computational results of traveltime data particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. Considering the complex geology of Korea, we assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution. The performance of the algorithm developed in this study is demonstrated by the comparison of the analytic and numerical solutions for the homogeneous anisotropic earth as well as through the numerical experiment for the two layer model whose anisotropic properties are greatly different each other. We expect that the developed modeling algorithm can be used in the development of processing and inversion schemes of seismic data acquired in strongly anisotropic environment, such as migration, velocity analysis, cross-well tomography and so on.

• Journal of the Korean Geotechnical Society
• /
• v.34 no.5
• /
• pp.19-35
• /
• 2018

The necessity of predicting the spatial information of the site-specific seismic response, which is essential information for the comprehensive earthquake disaster countermeasures, is increasing for the mid-west urban areas where the earthquake-induced damages can be increased due to frequent occurrence of mid-scale earthquake such as 2016 Gyeongju Earthquake and 2017 Pohang Earthquake. Especially, researches on strategic securing of site survey datasets and understanding the site-specific site response characteristics were conducted for Gyeonggi-do, South Korea. In this study, a GIS-based framework for site-specific assessment of site response and induced vulnerable area in Gyeonggi-do, South Korea was proposed. Geo-Data based on GIS platform was constructed for regional estimation of geotechnical characteristics by collecting borehole and land coverage datasets. And the geo-spatial grid information was developed for deriving spatial distribution of geotechnical layer and site response parameters based on the optimization of the geostatistical interpolation method. Accordingly, base information for Improving earthquake preparedness measures was derived as seismic zonation map with administrative sub-units considering the quantitative site effect of Gyeonggi-do.

• Journal of the Korean GEO-environmental Society
• /
• v.15 no.3
• /
• pp.65-74
• /
• 2014

Borehole investigation which is mainly used to figure out geotechnical characterizations at construction work has the benefit that it provides a clear and convincing geotechnical information. But it has limitations to get the overall information of the construction site because it is performed at point location. In contrast, geophysical measurements like seismic survey has the advantage that the geological stratum information of a large area can be characterized in a continuous cross-section but the result from geophysics survey has wide range of values and is not suitable to determine the geotechnical design values directly. Therefore it is essential to combine borehole data and geophysics data complementally. Accordingly, in this study, a three-dimensional spatial interpolation of the cross-sectional distribution of seismic refraction was performed using digitizing and geostatistical method (krigring). In the process, digital map were used to increase the trustworthiness of method. Using this map, errors of ground height which are broken out in measurement from boring investigation and geophysical measurements can be revised. After that, average seismic velocity are derived by comparing borehole data with geophysical speed distribution data of each soil layer. During this process, outlier analysis is adapted. On the basis of the average seismic velocity, integrated analysis techniques to determine the three-dimensional geological stratum information is established. Finally, this analysis system is applied to dam construction field.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.26-32
• /
• 1999

Among the various seismic data processing sequences, the velocity analysis is the most time consuming and man-hour intensive processing steps. For the production seismic data processing, a good velocity analysis tool as well as the high performance computer is required. The tool must give fast and accurate velocity analysis. There are two different approches in the velocity analysis, batch and interactive. In the batch processing, a velocity plot is made at every analysis point. Generally, the plot consisted of a semblance contour, super gather, and a stack pannel. The interpreter chooses the velocity function by analyzing the velocity plot. The technique is highly dependent on the interpreters skill and requires human efforts. As the high speed graphic workstations are becoming more popular, various interactive velocity analysis programs are developed. Although, the programs enabled faster picking of the velocity nodes using mouse, the main improvement of these programs is simply the replacement of the paper plot by the graphic screen. The velocity spectrum is highly sensitive to the presence of the noise, especially the coherent noise often found in the shallow region of the marine seismic data. For the accurate velocity analysis, these noise must be removed before the spectrum is computed. Also, the velocity analysis must be carried out by carefully choosing the location of the analysis point and accuarate computation of the spectrum. The analyzed velocity function must be verified by the mute and stack, and the sequence must be repeated most time. Therefore an iterative, interactive, and unified velocity analysis tool is highly required. An interactive velocity analysis program, xva(X-Window based Velocity Analysis) was invented. The program handles all processes required in the velocity analysis such as composing the super gather, computing the velocity spectrum, NMO correction, mute, and stack. Most of the parameter changes give the final stack via a few mouse clicks thereby enabling the iterative and interactive processing. A simple trace indexing scheme is introduced and a program to nike the index of the Geobit seismic disk file was invented. The index is used to reference the original input, i.e., CDP sort, directly A transformation techinique of the mute function between the T-X domain and NMOC domain is introduced and adopted to the program. The result of the transform is simliar to the remove-NMO technique in suppressing the shallow noise such as direct wave and refracted wave. However, it has two improvements, i.e., no interpolation error and very high speed computing time. By the introduction of the technique, the mute times can be easily designed from the NMOC domain and applied to the super gather in the T-X domain, thereby producing more accurate velocity spectrum interactively. The xva program consists of 28 files, 12,029 lines, 34,990 words and 304,073 characters. The program references Geobit utility libraries and can be installed under Geobit preinstalled environment. The program runs on X-Window/Motif environment. The program menu is designed according to the Motif style guide. A brief usage of the program has been discussed. The program allows fast and accurate seismic velocity analysis, which is necessary computing the AVO (Amplitude Versus Offset) based DHI (Direct Hydrocarn Indicator), and making the high quality seismic sections.

• Geophysics and Geophysical Exploration
• /
• v.20 no.2
• /
• pp.88-95
• /
• 2017

Chogye basin, which is surrounded by country rock, has a closed-basin form. In such a basin, incident seismic energy can form multiply reflected waves, thus causing energy concentration to occur at this closed-basin area. Microtremor measurement survey was performed at the Chogye basin, which is located in Chogye-myeon and Jeokjungmyeon, Hapcheon-gun, Gyeongsangnam-do, Republic of Korea. Microtremor data were transformed into the frequency domain, and then the horizontal-to-vertical spectral ratios (HVSR) were calculated. Fundamental resonance frequencies were estimated from the HVSR results for every observation point. Using the empirical relationship between site period and thickness for sediment sites in Korea known from the previous study, the distribution of sediment thickness of the Chogye basin was estimated from the fundamental resonance frequencies. Being compared with the mountainous rim with steep slope, the measurement points inside the basin have low values of the fundamental resonance frequency with the minimum of 1.03 Hz, which corresponds to the thickness of sedimentary layer with the maximum depth of about 100 m. A three-dimensional basin model was constructed for bedrock topography of the Chogye basin by an interpolation of basin depths estimated at each measurement site.